<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40"> <head> <meta http-equiv="Content-Type" content="text/html; charset=utf-8"> <meta name="Generator" content="Microsoft Word 15 (filtered medium)"> <style></style> </head> <body lang="EN-US" link="blue" vlink="purple"> <div class="WordSection1"> Dear Zuoheng<o:p></o:p> <o:p> </o:p> Though the regression test inputs may serve as a kind of initial input template for an application case, the settings for these tests are rather sloppy allowing for a quick processing in order to check a specific program feature or functionality. The core-shell model requires smaller time steps than rigid-ion force fields because of the usually faster motion of the light shell particles. At ambient conditions 0.5 fs might be a good guess, but smaller and smaller time steps are required for increasing temperatures until the break-down of the approach (“polarisation catastrophe”) at high temperatures will occur. The shell particles should always stay close to their core particles. Core-shell distances of more than 0.5A indicate that the approach is going to break-down.<o:p></o:p> There is the option to minimise (relax) the core-shell distances (<a href="https://manual.cp2k.org/trunk/CP2K_INPUT/MOTION/SHELL_OPT.html">MOTION/SHELL_OPT</a> section) if you specify for MAX_DISTANCE a value greater that zero and that value is surpassed. That is useful for initialisation or structure optimisation runs, but the energy conservation is affected when such relaxation steps are randomly performed during an MD run. MAX_DISTANCE has only an effect if the SHELL_OPT section is explicitly specified.<o:p></o:p> <o:p> </o:p> HTH<o:p></o:p> <o:p> </o:p> Matthias<o:p></o:p> <o:p> </o:p> From: cp...@googlegroups.com <...@googlegroups.com> On Behalf Of zuohen...@gmail.com Sent: Donnerstag, 20. September 2018 10:10 To: cp2k <...@googlegroups.com> Subject: Re: [CP2K:10764] Why doesn't the Cons Qty include the potential energy of shell (shell model) in NVE ensemble?<o:p></o:p> <o:p> </o:p> <div> Dear Matthias,<o:p></o:p> <div> <o:p> </o:p> </div> <div> Thank you for your reply.<o:p></o:p> </div> <div> Could I ask an another question? <o:p></o:p> </div> <div> I found core-shell model should use a smaller time step (0.1 - 0.2fs). However the time step of some U2O simulation in cp2k test files is larger (0.5fs even 1fs). <o:p></o:p> </div> <div> Is it because you use a heavier shell mass (MASS_FRACTION of O is 0.1) [1] and use a MAX_DISTANCE = 0.5 angstrom?<o:p></o:p> </div> <div> How do I determine the MAX_DISTANCE in my simulation of O2- and OH- ions with shell model?<o:p></o:p> </div> <div> <o:p> </o:p> </div> <div> Thank you very much~<o:p></o:p> </div> <div> <o:p> </o:p> </div> <div> Best wishes,<o:p></o:p> </div> <div> Zuoheng<o:p></o:p> </div> <div> <o:p> </o:p> </div> <div> [1] Devynck F, Iannuzzi M, Krack M. Frenkel pair recombinations in UO 2: importance of explicit description of polarizability in core-shell molecular dynamics simulations. Physical Review B, 2012, 85(18): 184103.<o:p></o:p> </div> 在 2018年9月19日星期三 UTC+2下午7:55:28，Matthias Krack写道：<o:p></o:p> <blockquote style="border:none;border-left:solid #CCCCCC 1.0pt;padding:0cm 0cm 0cm 6.0pt;margin-left:4.8pt;margin-top:5.0pt;margin-right:0cm;margin-bottom:5.0pt"> <div> <div> Dear Zuoheng<o:p></o:p> <o:p></o:p> md_ener%epot contains the potential energy of all particles including the shell particles.<o:p></o:p> <o:p></o:p> Best regards<o:p></o:p> <o:p></o:p> Matthias<o:p></o:p> <o:p></o:p> Von: <a href="javascript:" target="_blank">cp...@googlegroups.com</a> <<a href="javascript:" target="_blank">cp...@googlegroups.com</a>> Im Auftrag von <a href="javascript:" target="_blank">zuoh...@gmail.com</a> Gesendet: Mittwoch, 19. September 2018 17:15 An: cp2k <<a href="javascript:" target="_blank">cp...@googlegroups.com</a>> Betreff: [CP2K:10755] Why doesn't the Cons Qty include the potential energy of shell (shell model) in NVE ensemble?<o:p></o:p> <o:p></o:p> <div> <code>Dear CP2K developers,</code> <code>I am sorry to disturb you. I have some questions.</code> <code>Recently I want to analyze whether the total energy is converged with core-shell model in NVE ensemble. I didn't use extra thermostat for shell.</code> <code>I don't know why the Cons Qty energy = md_ener%ekin+md_ener%epot+md_ener%ekin_shell but not including the potential of shell? Is the Cons Qty the total energy? In NVE ensemble, which energy should be constant? ekin+epot+ekin_shell+epot_shell or ekin+epot+ekin_shell?</code> <code>Thank you very much~</code> <code>Best wishes,</code> <code>Zuoheng</code> <code>PS:The whole code of get_econs_nve function is following:</code><o:p></o:p> <div style="border:solid #BBBBBB 1.0pt;padding:0cm 0cm 0cm 0cm"> <div> SUBROUTINE get_econs_nve(md_env, md_ener, para_env) TYPE(md_environment_type), POINTER :: md_env TYPE(md_ener_type), INTENT(inout) :: md_ener TYPE(cp_para_env_type), POINTER :: para_env CHARACTER(LEN=*), PARAMETER :: routineN = 'get_econs_nve', routineP = moduleN//':'//routineN TYPE(force_env_type), POINTER :: force_env TYPE(thermostat_type), POINTER :: thermostat_coeff, thermostat_shell NULLIFY (force_env, thermostat_coeff, thermostat_shell) CALL get_md_env(md_env, force_env=force_env, thermostat_coeff=thermostat_coeff, & thermostat_shell=thermostat_shell) md_ener%constant = md_ener%ekin+md_ener%epot+md_ener%ekin_shell CALL get_thermostat_energies(thermostat_shell, md_ener%thermostat_shell_pot, & md_ener%thermostat_shell_kin, para_env) md_ener%constant = md_ener%constant+md_ener%thermostat_shell_kin+md_ener%thermostat_shell_pot END SUBROUTINE get_econs_nve<o:p></o:p> </div> </div> <o:p></o:p> </div> -- 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 <a href="javascript:" target="_blank">cp2k+u...@googlegroups.com</a>. To post to this group, send email to <a href="javascript:" target="_blank">cp...@googlegroups.com</a>. Visit this group at <a href="https://groups.google.com/group/cp2k" target="_blank"> https://groups.google.com/group/cp2k</a>. For more options, visit <a href="https://groups.google.com/d/optout" target="_blank"> https://groups.google.com/d/optout</a>.<o:p></o:p> </div> </div> </blockquote> </div> -- 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 <a href="mailto:cp2k+unsu...@googlegroups.com">cp2k+unsu...@googlegroups.com</a>. To post to this group, send email to <a href="mailto:cp...@googlegroups.com">cp...@googlegroups.com</a>. Visit this group at <a href="https://groups.google.com/group/cp2k">https://groups.google.com/group/cp2k</a>. For more options, visit <a href="https://groups.google.com/d/optout">https://groups.google.com/d/optout</a>.<o:p></o:p> </div> </body> </html>