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)