<div dir="ltr"><div><div><div>Dear Marcella,<br><br></div>Thank you very much for you suggestion. I will read these literatures carefully.<br><br></div>Best regards,<br></div>Huan<br><br><br><div><div><div><div><div class="gmail_extra">
<br><br><div class="gmail_quote">On Tue, Oct 29, 2013 at 7:59 AM, Marcella Iannuzzi <span dir="ltr"><<a href="mailto:marc...@pci.uzh.ch" target="_blank">marc...@pci.uzh.ch</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div style="word-wrap:break-word"><div><span style="text-indent:0px;letter-spacing:normal;font-variant:normal;text-align:-webkit-auto;font-style:normal;font-weight:normal;line-height:normal;border-collapse:separate;text-transform:none;font-size:medium;white-space:normal;font-family:Helvetica;word-spacing:0px"><span style="border-collapse:separate;font-variant:normal;letter-spacing:normal;line-height:normal;text-align:-webkit-auto;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px"><div style="word-wrap:break-word">
<div><br></div></div></span></span>Dear Huan,
</div>
<br><div><div>I think you should not use methods blindly, and read first the related literature. </div><div>In the specific case, you should read more on the Martyna-Tuckerman Poisson solver and more on the GPW method, </div>
<div>in order to understand how to prepare your setup for the simulations. </div><div><div>The correct references are even given in the cp2k reference manual:</div><div><a name="142031673fd70063_14202cbc0480d3e5_reference_46" style="text-indent:0px;letter-spacing:normal;font-variant:normal;text-align:start;font-style:normal;font-weight:normal;line-height:normal;text-transform:none;font-size:medium;white-space:normal;font-family:Times;word-spacing:0px">Martyna, GJ; Tuckerman, ME. </a><a href="http://dx.doi.org/10.1063/1.477923" style="font-family:Times;font-size:medium;font-style:normal;font-variant:normal;font-weight:normal;letter-spacing:normal;line-height:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px" target="_blank">JOURNAL OF CHEMICAL PHYSICS, 110 (6), 2810-2821 (1999).</a></div>
<div><span style="font-family:Times">A reciprocal space based method for treating long range interactions in ab initio and force-field-based calculations in clusters.</span></div><div><span style="font-family:Times"><br>
</span></div>
<div><a name="142031673fd70063_14202cbc0480d3e5_reference_31" style="text-indent:0px;letter-spacing:normal;font-variant:normal;text-align:start;font-style:normal;font-weight:normal;line-height:normal;text-transform:none;font-size:medium;white-space:normal;font-family:Times;word-spacing:0px">VandeVondele, J; Krack, M; Mohamed, F; Parrinello, M; Chassaing, T; Hutter, J.</a><a href="http://dx.doi.org/10.1016/j.cpc.2004.12.014" style="font-family:Times;font-size:medium;font-style:normal;font-variant:normal;font-weight:normal;letter-spacing:normal;line-height:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px" target="_blank">COMPUTER PHYSICS COMMUNICATIONS, 167 (2), 103-128 (2005).</a><br style="text-indent:0px;letter-spacing:normal;font-variant:normal;text-align:start;font-style:normal;font-weight:normal;line-height:normal;text-transform:none;font-size:medium;white-space:normal;font-family:Times;word-spacing:0px">
<span style="text-indent:0px;letter-spacing:normal;font-variant:normal;text-align:start;font-style:normal;display:inline!important;font-weight:normal;float:none;line-height:normal;text-transform:none;font-size:medium;white-space:normal;font-family:Times;word-spacing:0px">QUICKSTEP: Fast and accurate density functional calculations using a mixed Gaussian and plane waves approach.</span></div>
</div><div><span style="text-indent:0px;letter-spacing:normal;font-variant:normal;text-align:start;font-style:normal;display:inline!important;font-weight:normal;float:none;line-height:normal;text-transform:none;font-size:medium;white-space:normal;font-family:Times;word-spacing:0px"><br>
</span></div><div>In particular, there you can find how to choose the cell size for the MT poisson solver, </div><div>and how multi grids and relative cutoff are used within the GPW scheme. </div><div><br></div><div>regards</div>
<div>Marcella</div><div><div><font face="Times"><br></font></div><div><font face="Times"><br></font></div><blockquote type="cite"><div dir="ltr">The structure of my benzene-water cluster was optimized by Gaussian09 and can be put into a cell with size as 5*5*5. I tried 7*7*7, 10*10*10, 15*15*15 and 20*20*20 cells for NVE calculation in CP2K with PBE functional. However, only cell with 20*20*20 can maintain the structure of benzene-water cluster. For the rest of results with smaller cell, all the structures break due to the dramatic increasing of temperature from 50 K to 10000 K after 30 steps.<br>
<br>I have no idea about these testing results. Why could not use smaller cell?<br><br>I consulted the CP2K reference and found a description on &MT section. <b><i>Sets up parameters of Martyna-Tuckerman poisson solver. Note that exact results are only guaranteed if the unit cell is twice
as large as charge density (and serious artefacts can result if the cell is much smaller). <br></i></b><br></div></blockquote><div><br></div></div><blockquote type="cite"><div><div dir="ltr">In addition, I noticed that the REL_CUTOFF value is related to the &MGRID section. <b><i>Specify the multiplicative factor for the CUTOFF keyword in MULTI_GRID section. The result gives the cutoff at which the
1/r non-periodic FFT3D is evaluated.Default is 2.0</i></b><br><b><i><br></i></b><i></i>However, both of these lines do not have any detail explanations or examples. The CP2K reference is really too poor for a beginner to learn how to use CP2K, comparing the Gaussian reference. I do not understand what is the meaning of "if the cell is twice as large as charge density". How to check the charge density of my system? And also, is there any rule of thumb to choose the value of REL_CUTOFF?<br>
<br>Best regards,<br>Huan<br><br><br><br><br><b><i></i></b></div><div><br></div>
-- <br></div>
You received this message because you are subscribed to the Google Groups "cp2k" group.<br>
To unsubscribe from this group and stop receiving emails from it, send an email to <a href="mailto:cp2k+uns...@googlegroups.com" target="_blank">cp2k+uns...@googlegroups.com</a>.<div><br>
To post to this group, send email to <a href="mailto:cp...@googlegroups.com" target="_blank">cp...@googlegroups.com</a>.<br>
Visit this group at <a href="http://groups.google.com/group/cp2k" target="_blank">http://groups.google.com/group/cp2k</a>.<br>
For more options, visit <a href="https://groups.google.com/groups/opt_out" target="_blank">https://groups.google.com/groups/opt_out</a>.<br>
</div></blockquote></div><br></div><div><div>
<p></p>
-- <br>
You received this message because you are subscribed to a topic in the Google Groups "cp2k" group.<br>
To unsubscribe from this topic, visit <a href="https://groups.google.com/d/topic/cp2k/Rt1JAuuZKxg/unsubscribe" target="_blank">https://groups.google.com/d/topic/cp2k/Rt1JAuuZKxg/unsubscribe</a>.<br>
To unsubscribe from this group and all its topics, send an email to <a href="mailto:cp2k%2Bun...@googlegroups.com" target="_blank">cp2k+uns...@googlegroups.com</a>.<br>
To post to this group, send email to <a href="mailto:cp...@googlegroups.com" target="_blank">cp...@googlegroups.com</a>.<br>
Visit this group at <a href="http://groups.google.com/group/cp2k" target="_blank">http://groups.google.com/group/cp2k</a>.<br>
For more options, visit <a href="https://groups.google.com/groups/opt_out" target="_blank">https://groups.google.com/groups/opt_out</a>.<br>
</div></div></blockquote></div><br></div></div></div></div></div></div>