Thanks Duy, <br>I appreciate it very much. I did that. I wanted to see how the units were converted. <br><br>Best Regards,<br>Vikas<br><br><div class="gmail_quote">On Tue, Jun 21, 2011 at 5:52 PM, Duy Le <span dir="ltr">&lt;<a href="mailto:ttduyle@gmail.com">ttduyle@gmail.com</a>&gt;</span> wrote:<br>
<blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">I don&#39;t remember. You can always calculate it easily:<br>
the ratio between C6 (or R0) of Cacbon in mm_dispersion.f90 and in<br>
dft-d3 code will give you an idea how to convert them.<br>
<div class="im">--------------------------------------------------<br>
Duy Le<br>
PhD Candidate<br>
Department of Physics<br>
University of Central Florida.<br>
<br>
&quot;Men don&#39;t need hand to do things&quot;<br>
<br>
<br>
<br>
</div><div><div></div><div class="h5">On Tue, Jun 21, 2011 at 3:19 PM, Vikas Varshney &lt;<a href="mailto:vv0210@gmail.com">vv0210@gmail.com</a>&gt; wrote:<br>
&gt; Dear Duy, All,<br>
&gt; I have converted the R0 from the dftd3.f code to QE code (from Angs to Bohr<br>
&gt; with 1.1 further scaling factor). However, I dont know what is the scaling<br>
&gt; factor that is used for C6 in mm_dispersion.f90. Any insight will be highly<br>
&gt; appreciated.<br>
&gt;<br>
&gt; Best Regards,<br>
&gt;<br>
&gt; Vikas Varshney,<br>
&gt; Computational Materials Scientist,<br>
&gt; Wright-Patterson Air Force Base, Dayton-OH (USA).<br>
&gt;<br>
&gt; On Tue, Jun 21, 2011 at 11:27 AM, Vikas Varshney &lt;<a href="mailto:vv0210@gmail.com">vv0210@gmail.com</a>&gt; wrote:<br>
&gt;&gt;<br>
&gt;&gt; Dear Duy,<br>
&gt;&gt; Thank you very much for the input. I will look into that today itself.<br>
&gt;&gt;<br>
&gt;&gt; Best Regards,<br>
&gt;&gt; Vikas<br>
&gt;&gt;<br>
&gt;&gt; On Tue, Jun 21, 2011 at 10:51 AM, Duy Le &lt;<a href="mailto:ttduyle@gmail.com">ttduyle@gmail.com</a>&gt; wrote:<br>
&gt;&gt;&gt;<br>
&gt;&gt;&gt; Dear Vikas,<br>
&gt;&gt;&gt; One simplest thing you can do is to read the DFT-D2 part of DFT-D3<br>
&gt;&gt;&gt; code provided by Grimme&#39;s group (<a href="http://toc.uni-muenster.de/DFTD3/" target="_blank">http://toc.uni-muenster.de/DFTD3/</a>)<br>
&gt;&gt;&gt; You will be able to find all necessary coefficients for all atoms with<br>
&gt;&gt;&gt; Z&lt;=86. Then you have to add them to mm_dispersion.f90. Recompile QE,<br>
&gt;&gt;&gt; then you are good to go.<br>
&gt;&gt;&gt; --------------------------------------------------<br>
&gt;&gt;&gt; Duy Le<br>
&gt;&gt;&gt; PhD Candidate<br>
&gt;&gt;&gt; Department of Physics<br>
&gt;&gt;&gt; University of Central Florida.<br>
&gt;&gt;&gt;<br>
&gt;&gt;&gt; &quot;Men don&#39;t need hand to do things&quot;<br>
&gt;&gt;&gt;<br>
&gt;&gt;&gt;<br>
&gt;&gt;&gt;<br>
&gt;&gt;&gt; On Tue, Jun 21, 2011 at 10:40 AM, Vikas Varshney &lt;<a href="mailto:vv0210@gmail.com">vv0210@gmail.com</a>&gt;<br>
&gt;&gt;&gt; wrote:<br>
&gt;&gt;&gt; &gt; Dear All,<br>
&gt;&gt;&gt; &gt;<br>
&gt;&gt;&gt; &gt; I am trying to relax my WS2 unit cell structure using vc-relax (using<br>
&gt;&gt;&gt; &gt; PBE<br>
&gt;&gt;&gt; &gt; pps) and I want to include disperson corrections to it. So, I tried<br>
&gt;&gt;&gt; &gt; using<br>
&gt;&gt;&gt; &gt; the london flag = .true. However, on starting my job I got an error<br>
&gt;&gt;&gt; &gt;<br>
&gt;&gt;&gt; &gt;<br>
&gt;&gt;&gt; &gt;  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%<br>
&gt;&gt;&gt; &gt;      from  init_london  : error #         1<br>
&gt;&gt;&gt; &gt;      atom W   not found<br>
&gt;&gt;&gt; &gt;<br>
&gt;&gt;&gt; &gt;  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%<br>
&gt;&gt;&gt; &gt;<br>
&gt;&gt;&gt; &gt; After going through mm_dispersion.f90, I found that the vdw_radius and<br>
&gt;&gt;&gt; &gt; c_6<br>
&gt;&gt;&gt; &gt; parameters are given for first 54 atoms only (which corresponds to<br>
&gt;&gt;&gt; &gt; first 5<br>
&gt;&gt;&gt; &gt; periods of periodic table). Tungsten (W), with atomic number of 74 is<br>
&gt;&gt;&gt; &gt; in 6th<br>
&gt;&gt;&gt; &gt; period.<br>
&gt;&gt;&gt; &gt;<br>
&gt;&gt;&gt; &gt; I emailed Prof. Stefan Grimme (author of the references that are<br>
&gt;&gt;&gt; &gt; mentioned<br>
&gt;&gt;&gt; &gt; in mm_dispersion.f90) if anything was further pursued. He replied<br>
&gt;&gt;&gt; &gt; mentioning<br>
&gt;&gt;&gt; &gt; that they have approached this problem further and referred me to his<br>
&gt;&gt;&gt; &gt; website,<br>
&gt;&gt;&gt; &gt;<br>
&gt;&gt;&gt; &gt; <a href="http://toc.uni-muenster.de/DFTD3/" target="_blank">http://toc.uni-muenster.de/DFTD3/</a><br>
&gt;&gt;&gt; &gt;<br>
&gt;&gt;&gt; &gt; and told me that this has already been implemented into major QM codes.<br>
&gt;&gt;&gt; &gt;<br>
&gt;&gt;&gt; &gt; Mentioning that, I have few questions.<br>
&gt;&gt;&gt; &gt; 1. I wanted to ask if QE has incorporated the additions DFT-D,<br>
&gt;&gt;&gt; &gt; recently.<br>
&gt;&gt;&gt; &gt; 2. Are there any other ways to incorporate dispersion interactions in<br>
&gt;&gt;&gt; &gt; calculations in DFT.<br>
&gt;&gt;&gt; &gt;<br>
&gt;&gt;&gt; &gt; Currently, I am running my simulations on 4.2 version of the code that<br>
&gt;&gt;&gt; &gt; I<br>
&gt;&gt;&gt; &gt; downloaded few months ago. I would high appreciate your insights in how<br>
&gt;&gt;&gt; &gt; should I go about solving this problem regarding dispersion<br>
&gt;&gt;&gt; &gt; interactions in<br>
&gt;&gt;&gt; &gt; Tungsten.<br>
&gt;&gt;&gt; &gt;<br>
&gt;&gt;&gt; &gt; Looking forward for a positive reply.<br>
&gt;&gt;&gt; &gt;<br>
&gt;&gt;&gt; &gt; Best Regards,<br>
&gt;&gt;&gt; &gt;<br>
&gt;&gt;&gt; &gt; Vikas Varshney,<br>
&gt;&gt;&gt; &gt; Computational Materials Scientist,<br>
&gt;&gt;&gt; &gt; Wright-Patterson Air Force Base, Dayton-OH (USA).<br>
&gt;&gt;&gt; &gt;<br>
&gt;&gt;&gt; &gt; _______________________________________________<br>
&gt;&gt;&gt; &gt; Pw_forum mailing list<br>
&gt;&gt;&gt; &gt; <a href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a><br>
&gt;&gt;&gt; &gt; <a href="http://www.democritos.it/mailman/listinfo/pw_forum" target="_blank">http://www.democritos.it/mailman/listinfo/pw_forum</a><br>
&gt;&gt;&gt; &gt;<br>
&gt;&gt;&gt; &gt;<br>
&gt;&gt;&gt; _______________________________________________<br>
&gt;&gt;&gt; Pw_forum mailing list<br>
&gt;&gt;&gt; <a href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a><br>
&gt;&gt;&gt; <a href="http://www.democritos.it/mailman/listinfo/pw_forum" target="_blank">http://www.democritos.it/mailman/listinfo/pw_forum</a><br>
&gt;&gt;<br>
&gt;<br>
&gt;<br>
&gt; _______________________________________________<br>
&gt; Pw_forum mailing list<br>
&gt; <a href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a><br>
&gt; <a href="http://www.democritos.it/mailman/listinfo/pw_forum" target="_blank">http://www.democritos.it/mailman/listinfo/pw_forum</a><br>
&gt;<br>
&gt;<br>
_______________________________________________<br>
Pw_forum mailing list<br>
<a href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a><br>
<a href="http://www.democritos.it/mailman/listinfo/pw_forum" target="_blank">http://www.democritos.it/mailman/listinfo/pw_forum</a><br>
</div></div></blockquote></div><br>