<div dir="ltr">Hi Sergio,<div><br><div>Is your system metallic or having very small band gap? You may have to do BOMD. If nothing else works you may use our development version of doing BMD which uses damped dynamics to optimize wavefunction. It is quite efficient and cheaper than BOMD in PW, at least for the systems we tested so far. We plan to put this in QE in the next version. We can send you our version and you may give us your feedback. Please send us an e-mail if you wish to try it.</div><div><br></div><div>Best wishes,<br><div class="gmail_extra"><br clear="all"><div><div class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div>Biswajit Santra</div><div>Mobile: +1-609-227-9202</div><div><a href="http://www.princeton.edu/~bsantra/" style="font-size:12.8px" target="_blank">http://www.princeton.edu/~bsantra/</a><br></div></div></div></div></div>
<br><div class="gmail_quote">On Mon, Apr 24, 2017 at 9:04 PM, Sergio Castaneda Ramirez <span dir="ltr"><<a href="mailto:scastanr@unal.edu.co" target="_blank">scastanr@unal.edu.co</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div>Dear users</div><div><br></div><div>I'm carrying out a dynamical simulation for a hydroxide ion in a hydrated segment of a conductive polymer using the cp.x code of Quantum Espresso V6.1. </div><div><br></div><div>This is the general procedure I have been followed: Firstly, I perform an electronic minimization with fixed ions and cell to reach an electronic ground state. Secondly, I carry out an ionic relaxation run using damped molecular dynamics. After that, I add a random displacement of 0.01 a.u. to take out the ions from equilibrium to perform equilibration (NVT ensemble) and production (NVE ensemble) runs using verlet algorithm for both electrons and ions. </div><div><br></div><div>However, I haven't be able to keep the adiabaticity of the system when I perform NVE and NVT simulations. In both cases, kinetic electronic energy always increases from magnitudes of 1e-6 eV to orders of 0.1-10 after few time steps, which is totally incorrect. I have tried with different electronic fictitious masses, time step sizes and pseudopotentials. I also performed a relaxation of the cell volume using vc-cp, add and thermostat for the electronic kinectic energy during the NVT runs, but nothing seems to work. What could you suggest to me to do?</div><div><br></div><div>Very thanks in Advance.</div><div><br></div><div>Sergio Castaneda</div><div>Universidad Nacional de Colombia</div><div>Medellin, Colombia</div><div><br></div><div><br><div><br></div><div><br><div><br></div><div><br></div></div></div></div>
<br>______________________________<wbr>_________________<br>
Pw_forum mailing list<br>
<a href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a><br>
<a href="http://pwscf.org/mailman/listinfo/pw_forum" rel="noreferrer" target="_blank">http://pwscf.org/mailman/<wbr>listinfo/pw_forum</a><br></blockquote></div><br></div></div></div></div>