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