def anharm_plot2(qdt, fig_width=9.0, fig_height=6.0, ymin=-1.5, ymax=1.0):
"""Lamb shifted anharmonicity plot"""
pl = Plotter(fig_width=fig_width, fig_height=fig_height)
fw0 = linspace(2e9, 7e9, 2000)
lsfw0 = array([sqrt(f * (f - 2 * qdt._get_Lamb_shift(f=f))) for f in fw0])
Ej = qdt._get_Ej_get_fq(fq=lsfw0)
E0, E1, E2 = qdt._get_transmon_energy_levels(Ej=Ej, n_energy=3)
anharm = (E2 - E1) - (E1 - E0)
E0p, E1p, E2p = qdt._get_lamb_shifted_transmon_energy_levels(Ej=Ej,
n_energy=3)
anharmp = (E2p - E1p) - (E1p - E0p)
line(lsfw0 / 1e9,
anharm / h / 1e9,
plotter=pl,
linewidth=0.5,
color="purple",
label=r"anharm")
line(lsfw0 / 1e9,
anharmp / h / 1e9,
plotter=pl,
linewidth=0.5,
color="black",
label=r"ls anharm")
line(fw0 / 1e9,
qdt._get_coupling(fw0) / qdt.max_coupling,
label=r"$G_a/2C$",
color="blue",
plotter=pl)
pl.xlabel = r"$E_J/E_C$"
pl.ylabel = r"$\Delta$ (GHz)"
pl.legend(loc='lower left')
return pl
def coupling_plot():
pl=Plotter(fig_width=6.0, fig_height=4.0)
fw0=linspace(4e9, 7e9, 2000)
line(fw0/1e9, qdt._get_coupling(fw0)/1e9, label=r"$G_a/2C$", color="blue", plotter=pl)
line(fw0/1e9, qdt._get_Lamb_shift(fw0)/1e9, label=r"$-B_a/2C$", color="red", plotter=pl)
line(fw0/1e9, idt._get_coupling(fw0)/4/1e9, label=r"$G_a^{IDT}/2C/4$", color="green", linewidth=1.0, plotter=pl)
pl.set_ylim(-1.0, 1.5)
pl.legend(loc='lower right')
return pl
def anharm_plot2(qdt, fig_width=9.0, fig_height=6.0, ymin=-1.5, ymax=1.0):
"""Lamb shifted anharmonicity plot"""
pl=Plotter(fig_width=fig_width, fig_height=fig_height)
fw0=linspace(2e9, 7e9, 2000)
lsfw0=array([sqrt(f*(f-2*qdt._get_Lamb_shift(f=f))) for f in fw0])
Ej=qdt._get_Ej_get_fq(fq=lsfw0)
E0, E1, E2=qdt._get_transmon_energy_levels(Ej=Ej, n_energy=3)
anharm=(E2-E1)-(E1-E0)
E0p, E1p, E2p=qdt._get_lamb_shifted_transmon_energy_levels(Ej=Ej, n_energy=3)
anharmp=(E2p-E1p)-(E1p-E0p)
line(lsfw0/1e9, anharm/h/1e9, plotter=pl, linewidth=0.5, color="purple", label=r"anharm")
line(lsfw0/1e9, anharmp/h/1e9, plotter=pl, linewidth=0.5, color="black", label=r"ls anharm")
line(fw0/1e9, qdt._get_coupling(fw0)/qdt.max_coupling, label=r"$G_a/2C$", color="blue", plotter=pl)
pl.xlabel=r"$E_J/E_C$"
pl.ylabel=r"$\Delta$ (GHz)"
pl.legend(loc='lower left')
return pl
