Bulk water MD, PM6-FM and full PBC

214 views
Skip to first unread message

Ivan Gladich

unread,
Sep 30, 2021, 7:19:32 AM9/30/21
to cp2k

Dear CP2K users

 I would like to perform a bulk water simulation using PM6-FM semiempirical method. My goal is to have a flavor of a chemical reaction between organics in bulk water and, afterward, refine it at DFT-MD level.

I am a newbie on semiempirical methods. I used PM6-FM because it seems appropriate for liquid system

(https://groups.google.com/g/cp2k/c/BjBMWoq0doE/m/QADgmxMhAQAJ )

 I started inspired by some examples in CP2k github. I made two attempts, but both give me unreasonable results

 1)    A water cluster in vacuum at 300 K. After a short while, the system "kind of freezes" and nothing is moving

 2)    A bulk water box in PBC: after a short while all hydrogens detach from the water oxygen

 Thus, Is it possible to perform semiempirical  (PM6-FM or other) MD in full PBC? Does anyone have a successful input to share or point that  I can use as a starting point?

 I attach my inputs if someone is interested to check.

The initial configuration was taken from an equilibrated water box of 1.9 nm X 1.9nm X 1.9nm equilibrated at  classical MD at 1 bar and 300 K

 Thank you very much for any possible help or suggestion

Best regards

Ivan

pm6-FM-MT.inp
wat-box.xyz
pm6-FM-PBC.inp

Ruel

unread,
Jun 12, 2022, 4:27:39 AM6/12/22
to cp2k
Hi Ivan, did you manage to get PM6-FM method to work with PBC? I have also the same problem of broken O..H bonds.
Ruel

Thomas Kühne

unread,
Jun 21, 2022, 5:10:34 AM6/21/22
to 'Dorothea Golze' via cp2k
Dear Ivan and Ruel, 

please find a sample input that we have used for a related system attached - maybe it’s helpful. 

Best, 
Thomas
PM6FM.inp

Ivan Gladich

unread,
Jun 21, 2022, 5:22:29 PM6/21/22
to cp...@googlegroups.com
Dear Thomas

Thank you very much
I am going to test it

Best and thank you again

Ivan


-- 
You received this message because you are subscribed to a topic in the Google Groups "cp2k" group.
To unsubscribe from this topic, visit https://groups.google.com/d/topic/cp2k/C5LWpGsPgEg/unsubscribe.
To unsubscribe from this group and all its topics, send an email to cp2k+uns...@googlegroups.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/cp2k/A6F6A878-7529-4505-8DB1-583797622CEF%40gmail.com.
<PM6FM.inp>

-- 
You received this message because you are subscribed to the Google Groups "cp2k" group.
To unsubscribe from this group and stop receiving emails from it, send an email to cp2k+uns...@googlegroups.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/cp2k/b8cda270-55da-4387-9831-1e07199f91b8n%40googlegroups.com.



==============================
Thomas D. Kühne
Dynamics of Condensed Matter
Chair of Theoretical Chemistry
University of Paderborn
Warburger Str. 100
D-33098 Paderborn
Germany


-- 
You received this message because you are subscribed to a topic in the Google Groups "cp2k" group.
To unsubscribe from this topic, visit https://groups.google.com/d/topic/cp2k/C5LWpGsPgEg/unsubscribe.
To unsubscribe from this group and all its topics, send an email to cp2k+uns...@googlegroups.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/cp2k/A6F6A878-7529-4505-8DB1-583797622CEF%40gmail.com.

Reply all
Reply to author
Forward
0 new messages