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')