@author: jornrietveld """ from FunctionsSpinHam import Adatom, Structure, TemperatureVariation, Measurement, getHigestProbabillity from time import time import matplotlib.pyplot as plt from matplotlib import transforms import numpy as np start = time() T_rangeShort = np.logspace(0.5,-0.5,100) V_bias = 0.002 Fe1 = Adatom(s=2,D=-2.1,E=0.31,g=2.11) Fe2 = Adatom(s=2,D=-3.6,E=0.31,g=2.11) Fe3 = Adatom(s=2,D=-2.1,E=0.311,g=2.11) JList = np.linspace(0.05,0.1,2) simDict = {'tipPos':0,'Gs':2.56e-7,'G':1.5e-9,'eta':0.5,'u':0.4} fig,ax = plt.subplots(1,3,figsize=(12,4.5)) scaledtrans = transforms.ScaledTranslation(-0.55, -0.15, fig.dpi_scale_trans) letters = ['a', 'b', 'c', 'd', 'e', 'f', 'g', 'h', 'i', 'j', 'k', 'l'] for i,axis in enumerate(ax): axis.text(0, 1, letters[i], fontsize=14, fontweight="bold", va="bottom", ha="left", transform=axis.transAxes + scaledtrans)
import matplotlib.pyplot as plt from mpl_toolkits.axes_grid1.inset_locator import inset_axes from matplotlib import transforms import numpy as np start = time() T_rangeShort = np.logspace(np.sqrt(2), 0, 20) V_bias = 0.001 #V_range = np.linspace(0.005,0.020,200) T = 1 FielDir = 'By' # 'Bx', 'By' or 'Bz' Fe1 = Adatom(s=2, D=-2.1, E=0.31, g=2.11) Fe2 = Adatom(s=2, D=-3.6, E=0.31, g=2.11) Fe3 = Adatom(s=2, D=-3.7, E=0.31, g=2.11) Fe4 = Adatom(s=2, D=-2.3, E=0.31, g=2.11) n = 625 BList = np.linspace(0, 11, 34) J = 0.7 JList = np.linspace(0.4, 2.2, 10) EList = np.linspace(0.16, 0.61, 10) simDict = { 'tipPos': 0,
for i, axis in enumerate(ax): axis.text(0, 1, letters[i], fontsize=14, fontweight="bold", va="bottom", ha="left", transform=axis.transAxes + scaledtrans) axis.set_ylabel(r'Switching rate ($s^{-1}$)') axis.set_xlabel(r'$\frac{1}{T}$ $(\frac{1}{k})$', fontsize=14) axis.grid() linestArray = ['k-', '--sC0', '--sC1', '--sC2', '--sC3'] Fe1 = Adatom(s=2, D=-2.1, E=0.31, g=2.11) Fe2 = Adatom(s=2, D=-3.6, E=0.31, g=2.11) Fe3 = Adatom(s=2, D=-2.1, E=0.311, g=2.11) Fe1D = Adatom(s=2, D=-3.1, E=0.31, g=2.11) Fe2D = Adatom(s=2, D=-4.6, E=0.31, g=2.11) Fe3D = Adatom(s=2, D=-3.1, E=0.311, g=2.11) Fe1E = Adatom(s=2, D=-2.1, E=0.61, g=2.11) Fe2E = Adatom(s=2, D=-3.6, E=0.61, g=2.11) Fe3E = Adatom(s=2, D=-2.1, E=0.611, g=2.11) J = 0.7 J2 = 1.2
from FunctionsSpinHam import Adatom, Structure, BigStructure, TemperatureVariation from time import time import matplotlib.pyplot as plt import numpy as np start = time() T_rangeShort = np.logspace(1, 0, 16) V_bias = 0.002 #Dividing by meV because I am stupid #Btip = 0.0 #Jtip = 2.11*Btip*muB/meV Fe1 = Adatom(s=2, D=-2.1, E=0.31, g=2.11) Fe2 = Adatom(s=2, D=-3.6, E=0.31, g=2.11) Fe3 = Adatom(s=2, D=-2.1, E=0.311, g=2.11) Fe4 = Adatom(s=2, D=-2.1, E=0.31, g=2.11) Fe5 = Adatom(s=2, D=-3.6, E=0.31, g=2.11) Fe6 = Adatom(s=2, D=-2.1, E=0.311, g=2.11) Fe7 = Adatom(s=2, D=-1.5, E=0.31, g=2.11) Jin = 0.65 Jout = 0.75 J = 0.7 J4 = J #0.43 J6 = J #0.65 J8 = J #1.03