Hi,
I am using AmberTools22, which was installed on our supercomputer last week along with Amber22 (no issues during installation and/or testing). So it should have all the necessary patches.
I just tested, the HMR is not the issue, similar behavior happens also without it at 2fs timestep.
Here is the tleap:
source leaprc.protein.ff19SB
source leaprc.RNA.OL3
WAT = O3P
source leaprc.water.opc3pol
loadAmberParams frcmod.opc3pol_HMR4fs
loadAmberParams frcmod.ionslm_126_opc
molc = loadPdb "input.pdb"
solvatebox molc OPC3BOX 10.0
addIonsRand molc Na+ 88
addIonsRand molc Cl- 92
saveAmberParm molc structure.prmtop structure.inpcrd
savePdb molc structure.pdb
quit
Cheers,
Piia
--
Piia Barto¹, Ph.D. (Pharmacy)
Project Researcher
University of Eastern Finland | UEF | School of Pharmacy
Kuopio Campus | Canthia
Yliopistonranta 1 C | P.O. BOX 1627 | 70211 Kuopio
+358 50 548 8376 | piia.bartos.uef.fi<mailto:piia.bartos.uef.fi>
From: Carlos Simmerling <carlos.simmerling.gmail.com<mailto:carlos.simmerling.gmail.com>>
Sent: maanantai 21. marraskuuta 2022 18.13
To: Piia Bartos <piia.bartos.uef.fi<mailto:piia.bartos.uef.fi>>; AMBER Mailing List <amber.ambermd.org<mailto:amber.ambermd.org>>
Subject: Re: [AMBER] OPC3-pol
can you share your leap input, and which version of AmberTools you are using? is it fully up to date with patches?
if you suspect HMR, did you try a run without HMR?
On Mon, Nov 21, 2022 at 6:21 AM Piia Bartos via AMBER <amber.ambermd.org<mailto:amber.ambermd.org>> wrote:
Hi Amber people,
I tried to test the new OPC3-pol water model with my protein system.
It seems to be incompatible with some of my simulation parameter, as the simulations with the OPC3-pol water appear to be "frozen", where atoms remain almost still thorough the simulation (500 ns). The same production script works fine for OPC water model. Any ideas what I need to change? The issue might be with the Parmed, HMR and using 4fs timestep, but I think I built the system correctly...
Here is the production script:
# 500 ns production Langevin thermostat + Monte carlo barostat dynamics
&cntrl
imin=0,
irest=1, ntx=5,
ntpr=25000, ntwx=25000, ntwr=25000,
ntb=2, barostat=2, cut=10, iwrap=1,
nstlim=125000000, nscm=0, dt=0.004,
ntt=3, gamma_ln = 1, ig=-1, temp0=310, tempi=310,
ntp=1,
ntc=2, ntf=2,
ioutfm=1,
/
&ewald
netfrc = 0,
skin_permit = 0.75,
&end
/
Cheers,
Piia
--
Piia Barto¹, Ph.D. (Pharmacy)
Project Researcher
University of Eastern Finland | UEF | School of Pharmacy
Kuopio Campus | Canthia
Yliopistonranta 1 C | P.O. BOX 1627 | 70211 Kuopio
+358 50 548 8376 | piia.bartos.uef.fi<mailto:piia.bartos.uef.fi><mailto:piia.bartos.uef.fi<mailto:piia.bartos.uef.fi>>
_______________________________________________
AMBER mailing list
AMBER.ambermd.org<mailto:AMBER.ambermd.org>
http://lists.ambermd.org/mailman/listinfo/amber<https://eur03.safelinks.protection.outlook.com/?url=http%3A%2F%2Flists.ambermd.org%2Fmailman%2Flistinfo%2Famber&data=05%7C01%7C%7C898a3ec4ab854e835bfd08dacbdb4cad%7C87879f2e73044bf2baf263e7f83f3c34%7C0%7C0%7C638046440022550254%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000%7C%7C%7C&sdata=i2y3ygYfEftmnTB%2BFUv2BEnSH13yScUZAMbnfJNKFdA%3D&reserved=0>
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Tue Nov 29 2022 - 01:00:02 PST