Dear users,
I am a Physics Masters student and am using tkwant to simulate braiding
(adiabatic exchanging) of Majorana fermions in a 2D system. Currently, even for
the smallest systems, although these do work using Jupyter, I get a memory
error when using my university's HPC supercomputer and MPI ranks. Also, when
running simulations on Jupyter notebook (supercomputer client) any large-ish
system never actually completes.
More specifically, my issue seems to lie with 'wave_function.evolve' and
'wave_function.evaluate' for the time for loop.
density_operator = kwant.operator.Density(syst, np.eye(8)).bind(params=params)
wave_function = manybody.WaveFunction(psi_init, tasks)
for time in times:
wave_function.evolve(time)
density = wave_function.evaluate(density_operator)
densities.append(density)
plt.plot(sites, density - density0, label='time={}'.format(time))
In my system, I have 4 time-dependent gates and they are supposed to turn on
and off according to what time the system is in. I have used 'if' statements to
specify when the gates should be on or off.
If anyone can help it would be of great help. Thank you.
