def figure3():
U = 10
se = chain.UniformChain(1, [0, 0], 0, [], U)
Es = np.linspace(-5, 15, 1024)
taus = np.linspace(0, 2, 1024)
figures.g2_coherent_tau(se, (0, 0), (1, 1), Es, taus, directory=se.directory(), ticks=[0, 0.5, 1, 1.5, 2], title="for $U = 10 \Gamma$")
def dimer_tau(js, t=1, U=1, taus=None):
se = chain.UniformChain(N=2, js=js, E=0, t=t, U=U)
se.label = r'$U={0}\Gamma$'.format(U)
omegas = np.linspace(-5 * t, 15 * t, 1024)
if taus is None:
taus = np.linspace(0, 4, 1024)
figures.g2_coherent_tau(se, (0, 0), (1, 1), omegas, taus, directory=se.directory(), ticks=[0.5, 0.75, 1, 1.25, 1.5], title="for $U = {}\Gamma$".format(U))
def ring_tau(t=1, U=1, taus=None, zticks=None):
se = ring.UniformRing(N=6, js=(0, 2), E=0, t=t, U=U)
se.label = r'$U={0}\Gamma$'.format(U)
omegas = np.linspace(-5 * t, 15 * t, 1024)
if taus is None:
taus = np.linspace(0, 4, 1024)
if zticks is None:
zticks = [0.5, 0.75, 1, 1.25, 1.5]
figures.g2_coherent_tau(se, (0, 0), (1, 1), omegas, taus, directory=se.directory(), ticks=zticks, title="for $U = {}\Gamma$".format(U))
def chain_tau(t=1, U=1, parasite=0, taus=None, zticks=None, logscale=False):
se = chain.UniformChain(N=10, js=(0, 9), E=0, t=t, U=U, parasite=parasite)
se.label = r'$U={0}\Gamma$'.format(U)
omegas = np.linspace(-5 * t, 5 * t, 1024)
if taus is None:
taus = np.linspace(0, 4, 1024)
if zticks is None:
zticks = [0.5, 0.75, 1, 1.25, 1.5]
figures.g2_coherent_tau(se, (0, 0), (1, 1),
omegas,
taus,
directory=se.directory(),
ticks=zticks,
title="for $U = {}\Gamma$".format(U),
logscale=logscale)
def cube_tau(t=1, U=1, taus=None, zticks=None, logscale=False):
se = cube.UniformCube(L=4, W=4, H=4, js=(0, 3), E=0, t=t, U=U)
se.label = r'$U={0}\Gamma$'.format(U)
omegas = np.linspace(-5 * t, 15 * t, 1024)
if taus is None:
taus = np.linspace(0, 4, 1024)
if zticks is None:
zticks = [0.5, 0.75, 1, 1.25, 1.5]
figures.g2_coherent_tau(se, (0, 0), (1, 1),
omegas,
taus,
directory=se.directory(),
ticks=zticks,
logscale=logscale,
title="for $U = {}\Gamma$".format(U))
def plane_tau(t=1,
U=1,
parasite=0,
phase=0,
taus=None,
zticks=None,
logscale=False):
"""tau dependence"""
se = plane.UniformPlane(L=8,
W=8,
js=(0, 7),
E=0,
t=t,
U=U,
phase=phase,
parasite=parasite)
se.label = r'$U={0}\Gamma$'.format(U)
omegas = np.linspace(-8 * t, 8 * t, 1024)
if taus is None:
taus = np.linspace(0, 4, 1024)
if zticks is None:
zticks = [0, 0.5, 1, 1.5, 2]
extra = None
if parasite > 0:
extra = 'parasite-{}'.format(parasite)
if phase > 0:
extra += 'phase-{}'.format(phase)
figures.g2_coherent_tau(se, (0, 0), (1, 1),
omegas,
taus,
directory=se.directory(extra),
ticks=zticks,
logscale=logscale,
title="for $U = {}\Gamma$".format(U),
verbose=True)