(You should click this button, you dont need to see the code).
In [1]:
%pylab inline
import sys
sys.path.append('../')
from mpl_toolkits.axes_grid1 import Grid
from Icarus.Classes.QuantumDot import QuantumDot
from constants import Constants
from Experiments.utils.plot_fidelity_curves import *
In [2]:
constants = Constants()
qd = QuantumDot(constants.xtau, constants.xxtau, constants.ptau, constants.FSS, constants.crosstau)
In [3]:
qd.crosstau = 0.
qd.bg_emission_rate = 0.
fss = np.linspace(-10, 10, 500)*1e-6
fidelities = [qd.ideal_fidelity_lorentzian(f)[0] for f in fss]
ghv = qd.ideal_fidelity_lorentzian(f)[1]
numerical_data_no_dec = np.loadtxt('../out/fidelity_vs_fss/xtau1/fss_v_fidelity.txt', delimiter = ',')
plt.figure(figsize = (16./1.5,9./1.5))
plt.plot(numerical_data_no_dec[:,0], numerical_data_no_dec[:,1], 'r-', fss/1e-6, fidelities, 'b--')
plt.legend(['Monte carlo', 'Theoretical'])
plt.text(6, 0.9, ' coherence ' + np.array(ghv).astype('|S4').tostring())
plt.xlim([-10, 10]) ; plt.ylim([0.45, 1.05])
plt.xticks(np.linspace(-10, 10, 11))
plt.xlabel('Fine structure splitting $eV$') ; plt.ylabel('Fidelity')
Out[3]:
In [4]:
qd.crosstau = 2.5
fss = np.linspace(-10, 10, 500)*1e-6
fidelities = [qd.ideal_fidelity_lorentzian(f)[0] for f in fss]
ghv = qd.ideal_fidelity_lorentzian(f)[1]
numerical_data_with_dec = np.loadtxt('../out/fidelity_vs_fss/cross_dephasing2.5ns/fss_v_fidelity.txt', delimiter = ',')
plt.figure(figsize = (16./1.5,9./1.5))
plt.plot(numerical_data_with_dec[:,0], numerical_data_with_dec[:,1], 'r-', fss/1e-6, fidelities, 'b--')
plt.legend(['Monte carlo', 'Theoretical'])
plt.text(6, 0.9, ' coherence ' + np.array(ghv).astype('|S4').tostring())
plt.xlim([-10, 10]) ; plt.ylim([0.45, 1.05])
plt.xticks(np.linspace(-10, 10, 11))
plt.xlabel('Fine structure splitting $eV$') ; plt.ylabel('Fidelity')
Out[4]:
In [6]:
qd.crosstau = 2.5
qd.bg_emission_rate = 0.05
fss = np.linspace(-10, 10, 500)*1e-6
fidelities = [qd.ideal_fidelity_lorentzian(f)[0] for f in fss]
ghv = qd.ideal_fidelity_lorentzian(f)[1]
numerical_data_dec_and_bg = np.loadtxt('../out/2013-08-07/fid_v_fss/fid_v_fss.txt', delimiter = ',')
plt.figure(figsize = (16./1.5,9./1.5))
plt.plot(numerical_data_dec_and_bg[:,0]/1e-6, numerical_data_dec_and_bg[:,1], 'r-', fss/1e-6, fidelities, 'b--')
plt.legend(['Monte carlo', 'Theoretical'])
plt.text(6, 0.9, ' coherence ' + np.array(ghv).astype('|S4').tostring())
plt.xlim([-10, 10]) ; plt.ylim([0.45, 1.05])
plt.xticks(np.linspace(-10, 10, 11))
plt.xlabel('Fine structure splitting $eV$') ; plt.ylabel('Fidelity')
Out[6]:
In [ ]: