<html>
<head>
<meta name="generator" content="Windows Mail 17.5.9600.20856">
<style data-externalstyle="true"><!--
p.MsoListParagraph, li.MsoListParagraph, div.MsoListParagraph {
margin-top:0in;
margin-right:0in;
margin-bottom:0in;
margin-left:.5in;
margin-bottom:.0001pt;
}
p.MsoNormal, li.MsoNormal, div.MsoNormal {
margin:0in;
margin-bottom:.0001pt;
}
p.MsoListParagraphCxSpFirst, li.MsoListParagraphCxSpFirst, div.MsoListParagraphCxSpFirst, 
p.MsoListParagraphCxSpMiddle, li.MsoListParagraphCxSpMiddle, div.MsoListParagraphCxSpMiddle, 
p.MsoListParagraphCxSpLast, li.MsoListParagraphCxSpLast, div.MsoListParagraphCxSpLast {
margin-top:0in;
margin-right:0in;
margin-bottom:0in;
margin-left:.5in;
margin-bottom:.0001pt;
line-height:115%;
}
--></style></head>
<body dir="ltr">
<div data-externalstyle="false" dir="ltr" style="font-family: 'Meiryo', 'Calibri', 'Segoe UI', 'Microsoft YaHei UI', 'Microsoft JhengHei UI', 'Malgun Gothic', 'sans-serif';font-size:11.4975pt;">

<div>Dear Dr. Paolo Giannozzi,</div><div><br></div><div>I am sorry for the delay in my reply.</div><div><br></div><div><div>> My suggestion for the next release is to align 3D-RISM-SCF to the QE version that will be released (this should be easy); </div><div>> to release 3D-RISM-SCF as a modified QE (not yet sure how in practice, if as a patch, or as a bunch of modified routines,</div><div>> or as new PW/src and Modules/ directories); then, to start seriously to think at how to merge it into QE. </div><div>> Unfortunately, merging takes a lot of time and nobody is really in charge of it (in the future this may change, hopefully).</div><div><br></div></div><div>As you point out, the 3D-RISM-SCF works inside the self-consistent loop.</div><div>And it may require careful consideration and enough time to merge the code.</div><div><div>I accept your suggestion to avoid the integration of the 3D-RISM-SCF for the next release of Q.E.</div></div><div>I hope it will be integrated in the future.</div><div><br></div><div><div><div>> About "custom" FFTs: these are used only in hybrid calculations, to speed up the computation of (V_x\psi)(r) using FFTs</div><div>> with a grid corresponding to a reduced cutoff (this is 1-2 times ecutwfc, the cutoff for plane waves used </div><div>> in the expansion of \psi(r)=\sum_G\psi(G) e^{iGr}, instead of 4*ecutwfc, needed in principle to perform the exact computation).</div><div>> I am surprised to hear that this is used in 3D-RISM-SCF (or are you planning to use them in the future)?</div></div><br></div><div>I know well the calculation of exact-exchange with reduced FFT’s mesh,</div><div>because I have implemented Hybrid-GGA to another first-principles software.</div><div>In the 3D-RISM-SCF, the purpose to use “Custom” FFTs is to speed up the computation  (same as exact-exchange).</div><div>Solvent’s distributions have to be calculated by “Dense” FFTs, strictly.</div><div>However, “Custom” FFTs with small cutoff yield little decreasing of precision and a lot of speed-up.<br></div><div data-signatureblock="true"><div>Thus, I have used “Custom” FFTs, and I will use them in the future.<br></div><div>If you change or repeal the implementation of “Custom” FFTs,</div><div>I would like for you to remain their versatility or </div><div>to prepare another generic modules of FFT with arbitrary cutoff.</div><div><br></div><div>Regards,</div><div>Satomichi Nishihara</div><div><br></div><div><br></div></div><div style="padding-top: 5px; border-top-color: rgb(229, 229, 229); border-top-width: 1px; border-top-style: solid;"><div><font face=" 'Meiryo', 'Calibri', 'Segoe UI', 'Microsoft YaHei UI', 'Microsoft JhengHei UI', 'Malgun Gothic', 'sans-serif'" style='line-height: 15pt; letter-spacing: 0.02em; font-family: "Meiryo", "Calibri", "Segoe UI", "Microsoft YaHei UI", "Microsoft JhengHei UI", "Malgun Gothic", "sans-serif"; font-size: 12pt;'><b>差出人:</b> <a href="mailto:paolo.giannozzi@uniud.it" target="_parent">Paolo Giannozzi</a><br><b>送信日時:</b> ‎2015‎年‎6‎月‎9‎日 ‎火曜日 ‎0‎:‎56<br><b>宛先:</b> <a href="mailto:nisihara225@gmail.com" target="_parent">nisihara225@gmail.com</a><br><b>CC:</b> <a href="mailto:q-e-developers@qe-forge.org" target="_parent">General discussion list for Quantum ESPRESSO developers</a>, <a href="mailto:spiga.filippo@gmail.com" target="_parent">Filippo S</a></font></div></div><div><br></div><div dir=""><div dir="ltr"><div><div>This is a rather complex issue and I am not sure it can be solved in the next few days. If I understand correctly, <br>3D-RISM-SCF modifies the way the self-consistent code works internally, by adding computation of additional <br>quantities inside the self-consistent loop. Until now this kind of changes has been added inside QE directly by <br>the author, but after many uncoordinated additions the result has become a mess. One idea that has been around<br>for a while is to locate some "hooks" in the code where additional computations can be confined; then release the<br>implementation of additional methods (like 3D-RISM-SCF) as "add-ons" that modify QE in those specific points only<br></div><div>(notice the plugins_* subroutines). Unfortunately this is still ongoing (or maybe "notgoing") work. <br><br>My suggestion for the next release is to align 3D-RISM-SCF to the QE version that will be released (this should be easy); <br>to release 3D-RISM-SCF as a modified QE (not yet sure how in practice, if as a patch, or as a bunch of modified routines,<br></div><div>or as new PW/src and Modules/ directories); then, to start seriously to think at how to merge it into QE. <br>Unfortunately, merging takes a lot of time and nobody is really in charge of it (in the future this may change, hopefully).<br><br>About "custom" FFTs: these are used only in hybrid calculations, to speed up the computation of (V_x\psi)(r) using FFTs<br></div>with a grid corresponding to a reduced cutoff (this is 1-2 times ecutwfc, the cutoff for plane waves used in the expansion of \psi(r)=\sum_G\psi(G) e^{iGr}, instead of 4*ecutwfc, needed in principle to perform the exact computation). I am surprised<br></div><div>to hear that this is used in 3D-RISM-SCF (or are you planning to use them in the future)?<br><br></div><div>Paolo<br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Jun 8, 2015 at 3:38 PM, Filippo Spiga <span dir="ltr"><<a href="mailto:spiga.filippo@gmail.com" target="_parent">spiga.filippo@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin: 0px 0px 0px 0.8ex; padding-left: 1ex; border-left-color: rgb(204, 204, 204); border-left-width: 1px; border-left-style: solid;">Dear Satomichi,<br>
<br>
thank you very much for your quick feedback. I will attempt myself to merge your branch into a local copy of current QE codebase and run the tests. If I get something weird, I will contact you directly. I suggest to sync by email again beginning next week, I will get in touch privately. From my perspective, once the code compiles and is aligned and test are provided then it can be released. I let senior QE developers and PIs comment further.<br>
<br>
Thank you to contributing to Quantum ESPRESSO, your effort is really appreciated!<br>
<br>
Regards,<br>
Filippo<br>
<div class="HOEnZb"><div class="h5"><br>
On Jun 8, 2015, at 1:02 PM, <a href="mailto:nisihara225@gmail.com" target="_parent">nisihara225@gmail.com</a> wrote:<br>
> Dear Dr. Filippo Spiga,<br>
><br>
> the answers to your questions are as follows.<br>
><br>
> > Which QE version is compatible with it?<br>
> The branch of 3D-RISM-SCF is compatible with the Q.E. 5.1.2.<br>
> In more detail, it has branched from the 11516th revision of SVN.<br>
> The URL of the branch is <a href="http://qeforge.qe-forge.org/svn/q-e/branches/espresso-3drismscf" target="_parent">http://qeforge.qe-forge.org/svn/q-e/branches/espresso-3drismscf</a> .<br>
><br>
> > Did you develop 3D-RISM-SCF as external plug-in or directly inside the QE code?<br>
> The 3D-RISM-SCF was implemented directly inside the Q.E. code.<br>
> Source files in Modules/ and PW/ were mainly modified.<br>
><br>
> > Do you have extensive tests to validate your implementation?<br>
> Yes, I prepare some tests or examples:<br>
> /branches/espresso-3drismscf/espresso/PW/examples/RISM_example .<br>
> And I confirm that these examples provide reasonable results.<br>
><br>
> > Do you have one or more reference papers that explains the method utilized?<br>
> Yes, main reference papers are<br>
> H.Sato et al., J. Chem. Phys. 2000, 112, 9463 , and<br>
> A.Kovalenko, F.Hirata, Chem. Phys. Lett. 1998, 290, 237 .<br>
> Some formulas are utilized in the 3D-RISM-SCF, and<br>
> their references are written in comments of corresponding source codes or manuals.<br>
><br>
> > Is your code suitable for general public release under GPL license?<br>
> Yes, my code is based on GPL license (<a href="http://www.gnu.org/copyleft/gpl.txt" target="_parent">http://www.gnu.org/copyleft/gpl.txt</a>).<br>
> I have wrriten so, in header comments of the source codes.<br>
><br>
> > Are you willingness to maintain and take care of bug fixes and improvements of this particular piece of code for at least one year starting from the release date?<br>
> Of course, I will maintain the code.<br>
><br>
> Regards,<br>
> Satomichi Nishihara<br>
<br>
</div></div><div class="HOEnZb"><div class="h5">--<br>
Mr. Filippo SPIGA, M.Sc.<br>
Quantum ESPRESSO Foundation<br>
<a href="http://fspiga.github.io" target="_parent">http://fspiga.github.io</a> ~ skype: filippo.spiga<br>
<br>
«Nobody will drive us out of Cantor's paradise.» ~ David Hilbert<br>
<br>
*****<br>
Disclaimer: "Please note this message and any attachments are CONFIDENTIAL and may be privileged or otherwise protected from disclosure. The contents are not to be disclosed to anyone other than the addressee. Unauthorized recipients are requested to preserve this confidentiality and to advise the sender immediately of any error in transmission."<br>
<br>
</div></div></blockquote></div><br><br clear="all"><br>-- <br><div class="gmail_signature"><div dir="ltr"><span><font color="#888888">Paolo Giannozzi, Dept. Chemistry&Physics&Environment,<br>
Univ. Udine, via delle Scienze 208, 33100 Udine, Italy<br>
hone <a href="tel:%2B39-0432-558216" target="_parent">+39-0432-558216</a>, fax <a href="tel:%2B39-0432-558222" target="_parent">+39-0432-558222</a><br>
</font></span></div></div>
</div>
</div>





</div>
</body>
</html>