<div dir="auto">I haven't used CASTEP, but from the description you mentioned, trust_radius_max/min is same. BFGS relaxation converge/finish when total force and change in total energy is below forc_conv_thr and etot_conv_thr. Through trust radius you can adjust what can be the max, min and initial atomic displacement system will make in each BFGS cycle to achieve converge.</div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Sun, 15 Nov 2020, 1:37 am Husak Michal, <<a href="mailto:Michal.Husak@vscht.cz">Michal.Husak@vscht.cz</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">> First of all, you can find all the information about units of data in pw.x input manual.<br>
<br>
For conv_thr the units are missing in the manual (but is is probably (a.u.)<br>
<br>
<br>
>Regarding the ecutrho, it is a very important factor for any DFT calculation,<br>
<br>
ecutrho equivalent is impossible to find in any CASTEP manual mot in related materials.<br>
<br>
I will ask CASTP developers directly ...<br>
<br>
<br>
According trust_radius_max - there is not clear description in manual how it works,  In manual is only written<br>
<br>
"<br>
<br>
Maximum ionic displacement in the structural relaxation.<br>
(bfgs only)<br>
"<br>
<br>
Does it mean  BFGS runs until all shifts are less than this criteria ?<br>
<br>
________________________________<br>
From: users <<a href="mailto:users-bounces@lists.quantum-espresso.org" target="_blank" rel="noreferrer">users-bounces@lists.quantum-espresso.org</a>> on behalf of Yuvam Bhateja <<a href="mailto:yuvamb16@gmail.com" target="_blank" rel="noreferrer">yuvamb16@gmail.com</a>><br>
Sent: Saturday, November 14, 2020 8:44:25 PM<br>
To: Quantum ESPRESSO users Forum<br>
Subject: Re: [QE-users] CASTEP to QE setup parameters conversions<br>
<br>
Hey Michal,<br>
<br>
First of all, you can find all the information about units of data in pw.x input manual. And yes, you can change the units and use, keep in mind about convergence criteria, QE uses total force where as some DFT codes uses max force present on a single atom.<br>
<br>
Regarding the ecutrho, it is a very important factor for any DFT calculation, some programs have different name, for example in SIESTA it is meshcuttoff. Please refer to input manual for detailed info of all parameters, it can be helpful.<br>
<br>
For your geom_disp_tol in CASTEP, you can check trust_radius_(min/max/ini) which offers the same purpose.<br>
<br>
All the energy units are in Rydberg, as mentioned in manual.<br>
<br>
Hope it helps.<br>
I am also facing similar problems by going from QE to SIESTA.<br>
<br>
Keep in touch<br>
<br>
Regards<br>
Yuvam<br>
<br>
On Sun, 15 Nov 2020, 12:43 am Husak Michal, <<a href="mailto:Michal.Husak@vscht.cz" target="_blank" rel="noreferrer">Michal.Husak@vscht.cz</a><mailto:<a href="mailto:Michal.Husak@vscht.cz" target="_blank" rel="noreferrer">Michal.Husak@vscht.cz</a>>> wrote:<br>
Hi all<br>
<br>
<br>
I was working a lot with CASTEP last years ...<br>
<br>
I want to switch some my projects to  QE + do head to head comparison with CASTEP.<br>
<br>
I am confused by some selected corresponding setup parameters conversion<br>
<br>
(I had handled the rest including input formats conversion myself):<br>
<br>
<br>
QE: etot_conv_thr  (a.u.]<br>
<br>
CASTEP: geom_energy_tol (eV/atom)<br>
<br>
Should I simply  recalculate eV to a.u. and multiply by number of atoms in unit cell ?<br>
<br>
Why the definition in QE/CASTEP differs ?<br>
<br>
<br>
QE: forc_conv_thr  (a.u.)<br>
<br>
CASTEP: geom_force_tol (eV/angstrem)<br>
<br>
I have no idea how to convert this  parameter ....<br>
<br>
<br>
QE: non existing criteria<br>
<br>
CASTEP: geom_disp_tol (angstrem) , Maximal atom displacement during a cycle ..<br>
<br>
Sound like QE does not check this criteria ?<br>
<br>
<br>
QE: ecutrho (Ry), default 4 * ecutwfc<br>
<br>
CASTEP: non existing parameter<br>
<br>
Any idea why CASTEP does not use this parameter ?<br>
<br>
<br>
QE: conv_thr  , units not found/missing in manual, probably (a.u.) ?<br>
<br>
CASTEP: elec_energy_tol (eV/atom)<br>
<br>
Should I simply  recalculate eV to a.u. and multiply by number of atoms in unit cell ?<br>
<br>
Why the definition differs ?<br>
<br>
<br>
Michal Husak<br>
<br>
UCT Prague, Department of Solid State Chemistry<br>
<br>
<br>
<br>
<br>
_______________________________________________<br>
Quantum ESPRESSO is supported by MaX (<a href="http://www.max-centre.eu" rel="noreferrer noreferrer" target="_blank">www.max-centre.eu</a><<a href="http://www.max-centre.eu" rel="noreferrer noreferrer" target="_blank">http://www.max-centre.eu</a>>)<br>
users mailing list <a href="mailto:users@lists.quantum-espresso.org" target="_blank" rel="noreferrer">users@lists.quantum-espresso.org</a><mailto:<a href="mailto:users@lists.quantum-espresso.org" target="_blank" rel="noreferrer">users@lists.quantum-espresso.org</a>><br>
<a href="https://lists.quantum-espresso.org/mailman/listinfo/users" rel="noreferrer noreferrer" target="_blank">https://lists.quantum-espresso.org/mailman/listinfo/users</a><br>
_______________________________________________<br>
Quantum ESPRESSO is supported by MaX (<a href="http://www.max-centre.eu" rel="noreferrer noreferrer" target="_blank">www.max-centre.eu</a>)<br>
users mailing list <a href="mailto:users@lists.quantum-espresso.org" target="_blank" rel="noreferrer">users@lists.quantum-espresso.org</a><br>
<a href="https://lists.quantum-espresso.org/mailman/listinfo/users" rel="noreferrer noreferrer" target="_blank">https://lists.quantum-espresso.org/mailman/listinfo/users</a><br>
</blockquote></div>