Dear all gromacs users,I'm new to gromacs, and I have some questions about bond constraints.Now I have a system with methane, hydrogen, and TIP4P/Ew water. In special, topology of hydrogen molecule is including [ constraints ], for constructing a virtual site in the center of hydrogen molecule.By including "constraint-algorithm = lincs" and "constraints = h-bonds", I found that all bonds including H were constrained.So I changed *.mdp file with "constraint-algorithm = lincs" and "constraints = none", and I found that LINCS are applied to hydrogen, but it seemed that SETTLE is applied to water instead of LINCS.If I want to apply LINCS both on hydrogen molecules and waters, and not on methane or other hydrocarbons including hydrogen atoms, how can I do that?Any advices will be very helpful.Dong Woo Kang -- Gromacs Users mailing list
* Please search the archive at http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting! * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists * For (un)subscribe requests visit https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a mail to gmx-users-requ...@gromacs.org.