def chain_tau0(t=1,
ylim1=None,
ylim2=None,
Es=None,
directory="",
yticks2=None,
parasite=0):
Us = [0.5, 1, 2, 10]
ses = [
chain.UniformChain(N=10,
js=(0, 9),
E=0,
t=t,
U=Ui * t,
parasite=parasite) for Ui in Us
]
for i in range(0, len(Us)):
ses[i].label = r'${}$'.format(ses[i].info['Us'][0])
if Es is None:
Es = np.linspace(-5 / 2 * t, 15 / 2 * t, 1024)
figures.g2_coherent(ses, (0, 0), (1, 1),
Es,
2 * Es,
ses[0].directory() + directory,
offset=0,
ylims=[ylim1, ylim2],
yticks2=yticks2)
def cube_tau0(t=1,
ylim1=None,
ylim2=None,
Es=None,
directory="",
yticks2=None):
Us = [0.5, 4, 20]
ses = [
cube.UniformCube(L=4, W=4, H=4, js=(0, 3), E=0, t=t, U=Ui * t)
for Ui in Us
]
for i in xrange(0, len(Us)):
ses[i].label = r'${}$'.format(ses[i].info['Us'][0] * t)
if Es is None:
Es = np.linspace(-5 / 2 * t, 25 / 2 * t, 1024)
figures.g2_coherent(ses, (0, 0), (1, 1),
Es,
2 * Es,
ses[0].directory() + directory,
offset=0,
ylims=[ylim1, ylim2],
yticks2=yticks2)
def figure1():
Us = [.5, 1, 2, 10]
ses = [chain.UniformChain(1, [0, 0], 0, [], Ui) for Ui in Us]
for i in range(len(Us)):
ses[i].label = r'${}$'.format(Us[i])
Es = np.linspace(-5 / 2, 15 / 2, 1024)
figures.g2_coherent(ses, (0, 0), (1, 1), Es, 2 * Es, ses[0].directory(), offset=0, ylims=[[1e-1, 2], [1e-1, 1e1]])
def dimer_tau0(js, t=1, ylim1=None, ylim2=None):
Us = [0.5, 1, 2, 10]
ses = [chain.UniformChain(N=2, js=js, E=0, t=t, U=Ui * t) for Ui in Us]
for i in range(0, len(Us)):
ses[i].label = r'${}$'.format(ses[i].info['Us'][0])
Es = np.linspace(-5 / 2 * t, 15 / 2 * t, 1024)
figures.g2_coherent(ses, (0, 0), (1, 1), Es, 2 * Es, ses[0].directory(), offset=0, ylims=[ylim1, ylim2])
def ring_tau0(t=1, ylim1=None, ylim2=None, Es=None, directory="", yticks2=None):
Us = [0.5, 1, 2, 10]
ses = [ring.UniformRing(N=6, js=(0, 2), E=0, t=t, U=Ui * t) for Ui in Us]
for i in range(0, len(Us)):
ses[i].label = r'${}$'.format(ses[i].info['Us'][0] * t)
if Es is None:
Es = np.linspace(-5 / 2 * t, 15 / 2 * t, 1024)
figures.g2_coherent(ses, (0, 0), (1, 1), Es, 2 * Es, ses[0].directory() + directory, offset=0, ylims=[ylim1, ylim2], yticks2=yticks2)
def plane_tau0(t=1,
parasite=0,
phase=0,
ylim1=None,
ylim2=None,
Es=None,
directory="",
yticks2=None):
"""generate g1 and g2 for plane"""
Us = [0.5, 1, 2, 10]
ses = [
plane.UniformPlane(L=8,
W=8,
js=(56, 0),
E=0,
t=t,
U=Ui * t,
phase=phase,
parasite=parasite) for Ui in Us
]
for i in range(0, len(Us)):
ses[i].label = r'${}$'.format(ses[i].info['Us'][0] * t)
if Es is None:
Es = np.linspace(-8 / 2 * t, 8 / 2 * t, 512)
extra = None
if parasite > 0:
extra = 'parasite-{}'.format(parasite)
if phase > 0:
extra += 'phase-{}'.format(phase)
figures.g2_coherent(ses, (0, 0), (1, 1),
Es,
2 * Es,
ses[0].directory(extra) + directory,
offset=0,
ylims=[ylim1, ylim2],
yticks2=yticks2)