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.