Spin 1/2 example¶

In [1]:

# Import Package
from spinfinity import *

import numpy as np
from matplotlib import pyplot as plt

# Import Package from spinfinity import * import numpy as np from matplotlib import pyplot as plt

In [2]:

# Create Spin System
spinsys = SpinSys([0.5])

# Modify Spin Hamiltonian Parameters
spinsys.B = [0, 0, 0.5] # Tesla

# Create Spin System spinsys = SpinSys([0.5]) # Modify Spin Hamiltonian Parameters spinsys.B = [0, 0, 0.5] # Tesla

ESR Frequency Sweep¶

In [3]:

f, I = ESR.calc_FreqSweep(spinsys, freq_array=np.linspace(13, 15, 1000));

plt.figure()
plt.plot(f, I)
plt.xlabel('Frequency (GHz)')
plt.ylabel('Intensity (a.u.)')
plt.title('ESR Spectrum')
plt.show()

f, I = ESR.calc_FreqSweep(spinsys, freq_array=np.linspace(13, 15, 1000)); plt.figure() plt.plot(f, I) plt.xlabel('Frequency (GHz)') plt.ylabel('Intensity (a.u.)') plt.title('ESR Spectrum') plt.show()

Found 1 potential transitions:
|0> -> |1>; f = 13.9953 GHz