<div dir="ltr">Thank you for your response.</div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Wed, Oct 18, 2023 at 6:59 PM Kazume NISHIDATE <<a href="mailto:nisidate@iwate-u.ac.jp">nisidate@iwate-u.ac.jp</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Dear Jayraj;<br>
<br>
> the temperature fluctuations are too high compared to the other<br>
>  literature.<br>
<br>
Because you are using "thousands of atoms" in LAMMPS while you are using<br>
only the 128 atoms in QE.<br>
<br>
> MD simulation of thousands of atoms in LAMMPS<br>
…….<br>
> On my metallic system, which has 128 atoms,<br>
<br>
<br>
You must understand the concept of ensemble average.<br>
A good text was introduced by Nicola:<br>
<br>
> 3) not sure where to start - you need to learn and understand MD - a<br>
>  good start is if you google "ercolessi md primer"<br>
<br>
<br>
In addition to his elemental comments, I want to add some trivial<br>
things.<br>
<br>
> "conv_thr" (i.e. ecutwfc=5 or 10) below a minimum cutoff of<br>
>  pseudopotential? So that I can get the next ionic position very<br>
<br>
Suppose you have a potential curvature having a shape like a cup and <br>
a ball is falling along its curvature from its edge. The ball will<br>
perform oscillating behavior continually if there is no energy loss by<br>
a friction (energy conservation).<br>
<br>
In a primitive molecular dynamics calculation (MD), total energy<br>
(kinetic+potential) must be preserved in the long time scale (say<br>
10000 time steps) when the system does not have a restriction on its<br>
kinetics (like temperature regulation). <br>
<br>
The "ecutwfc=5" apparently insufficient to evaluate the potential<br>
curvature and hence your MD calculation will collapse. It is <br>
equivalent to a MD in LAMMPS with random interaction potentials.<br>
<br>
Use the default cutoff energy and monitor the fluctuation of total<br>
energies in the MD runs with changing the dt (time increment value) to<br>
find the appropriate (largest) dt with preserving the total energies.<br>
<br>
And you can also consider a gamma-point only calculation if your <br>
system is large enough.<br>
<br>
<br>
<br>
西館数芽<br>
Kazume NISHIDATE Ph.D<br>
<br>
Department of Systems Innovation Engineering,<br>
Graduate School of Science and Engineering, Iwate University<br>
4-3-5 Ueda, Morioka, Iwate, 020-8551 JAPAN<br>
Phone:+81-19-621-6391<br>
<a href="mailto:kazume.nishidate@gmail.com" target="_blank">kazume.nishidate@gmail.com</a>, <a href="mailto:nisidate@iwate-u.ac.jp" target="_blank">nisidate@iwate-u.ac.jp</a><br>
<a href="https://sites.google.com/site/nisidatelab/" rel="noreferrer" target="_blank">https://sites.google.com/site/nisidatelab/</a><br>
<br>
_______________________________________________<br>
The Quantum ESPRESSO community stands by the Ukrainian<br>
people and expresses its concerns about the devastating<br>
effects that the Russian military offensive has on their<br>
country and on the free and peaceful scientific, cultural,<br>
and economic cooperation amongst peoples<br>
_______________________________________________<br>
Quantum ESPRESSO is supported by MaX (<a href="http://www.max-centre.eu" rel="noreferrer" target="_blank">www.max-centre.eu</a>)<br>
users mailing list <a href="mailto:users@lists.quantum-espresso.org" target="_blank">users@lists.quantum-espresso.org</a><br>
<a href="https://lists.quantum-espresso.org/mailman/listinfo/users" rel="noreferrer" target="_blank">https://lists.quantum-espresso.org/mailman/listinfo/users</a></blockquote></div>