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)