@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