[gmx-developers] Thermostatting in selected directions

Milan Predota predota at prf.jcu.cz
Thu Feb 2 13:54:31 CET 2017


I need to do NEMD simulations of flows and for correct thermostatting 
I'd need to select only 1 or 2 velocity directions used for setting the 

I need  e.g. to calculate  kinetic energy only from v_x and v_y and 
avoid v_z where there is a nonzero and a priory unknown streaming velocity.

The temperature coupling would act of course only on the v_x and v_y 
components, leaving the v_z to equilibrate due to intermolecular 

I have done this in my own MD code (only water + ions) with success, 

but I'm now using Gromacs for more complex molecules.

I don't find this feature even in Gromacs 2016. *Is there a plan to 
incorporate this feature in Gromacs?*

I guess to modify the source code by myself would be a tedious and 
dangerous task.

I have found out that LAMMPS probably has this feature

but I'd prefer to continue using Gromacs.
Many thanks,

Institute of Physics and Biophysics
Faculty of Science, University of South Bohemia
Branisovska 1760, 370 05 Ceske Budejovice
Czech Republic
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://maillist.sys.kth.se/pipermail/gromacs.org_gmx-developers/attachments/20170202/0a1bf32f/attachment.html>

More information about the gromacs.org_gmx-developers mailing list