def setup(U_input): system = harmonic.hhg(cycles=cycles, delta=delta, libN=libN, field=field, nup=number, ndown=number, nx=nx, ny=0, U=U_input, t=t, F0=F0, a=a, bc='pbc') return system
a = 4 scalefactor = 1 scalefactor2 = 1 ascale = 1.001 ascale2 = 1 Jscale = 1 cutoff = 100 cutoff2 = cutoff Tracking = True prop = hams.hhg(field=field, nup=number, ndown=number, nx=nx, ny=0, U=U, t=t, F0=F0, a=a, bc='pbc') prop2 = hams.hhg(field=field, nup=number, ndown=number, nx=nx, ny=0, U=U2, t=t2, F0=F0, a=a, bc='pbc') print(prop.field)
U = t * f delta = 0.05 cycles = 10 # field= 32.9 field = 32.9 F0 = 10 a = 4 parameternames = '-%s-nsites-%s-cycles-%s-U-%s-t-%s-n-%s-delta-%s-field-%s-amplitude' % ( nx, cycles, U, t, number, delta, field, F0) lat = harmonic.hhg(field=field, nup=number, ndown=number, nx=nx, ny=0, U=U, t=t, F0=F0, a=a, bc='pbc') print('\n') print(vars(lat)) time = cycles psi_temp = harmonic.hubbard(lat)[1].astype(complex) init = psi_temp h = hub.create_1e_ham(lat, True) N = int(time / (lat.freq * delta)) + 1 print(N) times = np.linspace(0.0, cycles, N) prop = lat r = ode(evolve.integrate_f).set_integrator('zvode', method='bdf')