diff --git a/nanowire-conductance.ipynb b/nanowire-conductance.ipynb index 69fd7e8d480bf2a7551d6028f0325f5917c703f9..f439bd3b102903cd161a0e5f8a3bb9358c8081a3 100755 --- a/nanowire-conductance.ipynb +++ b/nanowire-conductance.ipynb @@ -1,5 +1,12 @@ { "cells": [ + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "# Majorana nanowire conductance" + ] + }, { "cell_type": "code", "execution_count": null, @@ -10,7 +17,6 @@ "import kwant\n", "from scipy.constants import hbar, m_e, eV, physical_constants\n", "import matplotlib.pyplot as plt\n", - "%matplotlib inline\n", "\n", "meV = eV * 1e-3 # one meV in J\n", "mu_B = physical_constants['Bohr magneton'][0] / meV # Bohr magneton in meV\n", @@ -31,13 +37,17 @@ "\n", " # We apply a magnetic field in all parts of the system\n", " def onsite_sc(site, mu_l, mu_sc, B_x, alpha, delta, Vbarrier):\n", - " return (2 * t - mu_sc) * np.kron(s0, sz) + (0.5 * B_x * g * mu_B) * np.kron(sx, s0) + delta * np.kron(s0, sx)\n", + " return ((2 * t - mu_sc) * np.kron(s0, sz) \n", + " + (0.5 * B_x * g * mu_B) * np.kron(sx, s0) \n", + " + delta * np.kron(s0, sx))\n", "\n", " def onsite_normal(site, mu_l, mu_sc, B_x, alpha, delta, Vbarrier):\n", - " return (2 * t - mu_l) * np.kron(s0, sz) + (0.5 * B_x * g * mu_B) * np.kron(sx, s0)\n", + " return ((2 * t - mu_l) * np.kron(s0, sz) \n", + " + (0.5 * B_x * g * mu_B) * np.kron(sx, s0))\n", "\n", " def onsite_barrier(site, mu_l, mu_sc, B_x, alpha, delta, Vbarrier):\n", - " return (2 * t - mu_l + Vbarrier) * np.kron(s0, sz) + (0.5 * B_x * g * mu_B) * np.kron(sx, s0)\n", + " return ((2 * t - mu_l + Vbarrier) * np.kron(s0, sz) \n", + " + (0.5 * B_x * g * mu_B) * np.kron(sx, s0))\n", "\n", " # The hopping is the same in all subsystems. There is normal hopping and\n", " # spin orbit interaction.\n",