Hi Michele and Venkat,
Thank you both for your answers!
Michele, I thought about faking a 1-D simulation by choosing initial conditions so there is no momenta in y and z and using, for example the 1D harmonic potential driver for x. For an NVE simulation that should work.
However, that will not work with a thermostat, right? Is that what you meant?
Thanks,
Barak