<div dir="ltr">Replying to
Dr. Lorenzo Paulatto<div><br></div><div>Sir as you say that "algorithm fails because the cell constraint does not allow the cell</div><div>pressure to get under the threshold, you may need more vacuum or less<br>tight constraint." My current c value is 37.789979707283219 Boher. Its an hexagonal unit cell. What do you suggest should have been the c vale and convergence criterion.</div><div> </div><div>As you said " but there will be no net magnetization as you are not breaking the spin symmetry in the initial conditions." To observe that if my structure is showing magnetic behaviour or spin splitting in the bandstructure, what should have I given the initial magentization ?</div></div><br><div class="gmail_quote gmail_quote_container"><div dir="ltr" class="gmail_attr">On Sat, Apr 26, 2025 at 3:38 PM <<a href="mailto:users-request@lists.quantum-espresso.org">users-request@lists.quantum-espresso.org</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">Send users mailing list submissions to<br>
<a href="mailto:users@lists.quantum-espresso.org" target="_blank">users@lists.quantum-espresso.org</a><br>
<br>
To subscribe or unsubscribe via the World Wide Web, visit<br>
<a href="https://lists.quantum-espresso.org/mailman/listinfo/users" rel="noreferrer" target="_blank">https://lists.quantum-espresso.org/mailman/listinfo/users</a><br>
or, via email, send a message with subject or body 'help' to<br>
<a href="mailto:users-request@lists.quantum-espresso.org" target="_blank">users-request@lists.quantum-espresso.org</a><br>
<br>
You can reach the person managing the list at<br>
<a href="mailto:users-owner@lists.quantum-espresso.org" target="_blank">users-owner@lists.quantum-espresso.org</a><br>
<br>
When replying, please edit your Subject line so it is more specific<br>
than "Re: Contents of users digest..."<br>
<br>
<br>
Today's Topics:<br>
<br>
1. Incorrect Visualization of Ba?ZrS? Body-Centered Tetragonal<br>
Structure in XCrySDen from Quantum ESPRESSO Input (Asma Khatun)<br>
2. Re: [SPAM] Incorrect Visualization of Ba?ZrS?<br>
Body-Centered Tetragonal Structure in XCrySDen from Quantum<br>
ESPRESSO Input (Husak Michal)<br>
3. bfgs failed after 25 scf cycles and 23 bfgs steps,<br>
convergence not achieved (Barsha Pal)<br>
4. Re: bfgs failed after 25 scf cycles and 23 bfgs steps,<br>
convergence not achieved (Lorenzo Paulatto)<br>
5. Re: bfgs failed after 25 scf cycles and 23 bfgs steps,<br>
convergence not achieved (Paolo Giannozzi)<br>
<br>
<br>
----------------------------------------------------------------------<br>
<br>
Message: 1<br>
Date: Fri, 25 Apr 2025 22:54:27 +0600<br>
From: Asma Khatun <<a href="mailto:asma@tch.hstu.ac.bd" target="_blank">asma@tch.hstu.ac.bd</a>><br>
To: <a href="mailto:users@lists.quantum-espresso.org" target="_blank">users@lists.quantum-espresso.org</a><br>
Subject: [QE-users] Incorrect Visualization of Ba?ZrS? Body-Centered<br>
Tetragonal Structure in XCrySDen from Quantum ESPRESSO Input<br>
Message-ID:<br>
<CAOqCWc7sT0x7N=7q=81wOX=+CMeqv8t9BYHP0UTpcwbMDq=<a href="mailto:ucw@mail.gmail.com" target="_blank">ucw@mail.gmail.com</a>><br>
Content-Type: text/plain; charset="utf-8"<br>
<br>
Dear Quantum ESPRESSO Community,<br>
<br>
I am trying to model the Ba?ZrS? compound with a body-centered tetragonal<br>
(bct) structure using Quantum ESPRESSO. However, when I visualize the<br>
structure using XCrySDen from my SCF input file (<a href="http://ba2zrs4.scf.in" rel="noreferrer" target="_blank">ba2zrs4.scf.in</a>), the<br>
crystal does not appear as expected. The atomic positions and symmetry seem<br>
incorrect even before running any calculation.<br>
To help identify the issue, I have attached:<br>
<br>
a) <a href="http://ba2zrs4.scf.in" rel="noreferrer" target="_blank">ba2zrs4.scf.in</a>: My input file used for visualization.<br>
<br>
b) unitcell.png: The structure as visualized in XCrySDen.<br>
<br>
Could anyone kindly review my input file and advise if there is any mistake<br>
in the lattice vectors, atomic positions, or crystal symmetry definition?<br>
Any help would be sincerely appreciated.<br>
<br>
Thank you very much for your time and support.<br>
<br>
Best regards,<br>
Mst. Asma Khatun<br>
PhD Fellow, Department of Physics<br>
University of Rajshahi, Bangladesh<br>
<br>
whatsapp: +8801723674001<br>
<br>
e-mail: <a href="mailto:asma@tch.hstu.ac.bd" target="_blank">asma@tch.hstu.ac.bd</a><br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <<a href="http://lists.quantum-espresso.org/pipermail/users/attachments/20250425/50706c30/attachment-0001.html" rel="noreferrer" target="_blank">http://lists.quantum-espresso.org/pipermail/users/attachments/20250425/50706c30/attachment-0001.html</a>><br>
-------------- next part --------------<br>
A non-text attachment was scrubbed...<br>
Name: <a href="http://ba2zrs4.scf.in" rel="noreferrer" target="_blank">ba2zrs4.scf.in</a><br>
Type: application/octet-stream<br>
Size: 1574 bytes<br>
Desc: not available<br>
URL: <<a href="http://lists.quantum-espresso.org/pipermail/users/attachments/20250425/50706c30/attachment-0001.obj" rel="noreferrer" target="_blank">http://lists.quantum-espresso.org/pipermail/users/attachments/20250425/50706c30/attachment-0001.obj</a>><br>
-------------- next part --------------<br>
A non-text attachment was scrubbed...<br>
Name: unitcell.png.png<br>
Type: image/png<br>
Size: 123464 bytes<br>
Desc: not available<br>
URL: <<a href="http://lists.quantum-espresso.org/pipermail/users/attachments/20250425/50706c30/attachment-0001.png" rel="noreferrer" target="_blank">http://lists.quantum-espresso.org/pipermail/users/attachments/20250425/50706c30/attachment-0001.png</a>><br>
<br>
------------------------------<br>
<br>
Message: 2<br>
Date: Fri, 25 Apr 2025 21:37:19 +0000<br>
From: Husak Michal <<a href="mailto:Michal.Husak@vscht.cz" target="_blank">Michal.Husak@vscht.cz</a>><br>
To: Quantum ESPRESSO users Forum <<a href="mailto:users@lists.quantum-espresso.org" target="_blank">users@lists.quantum-espresso.org</a>><br>
Subject: Re: [QE-users] [SPAM] Incorrect Visualization of Ba?ZrS?<br>
Body-Centered Tetragonal Structure in XCrySDen from Quantum ESPRESSO<br>
Input<br>
Message-ID: <<a href="mailto:8d53e7a0da3b4123b2ad4d19ef2e482a@vscht.cz" target="_blank">8d53e7a0da3b4123b2ad4d19ef2e482a@vscht.cz</a>><br>
Content-Type: text/plain; charset="euc-kr"<br>
<br>
Can you supply a CIF file with space group, so we can see what you expect ?<br>
Can you check the CIF yourself by e.g. Mercury ?<br>
Automatic conversion from CIF is most reliable method for most DFT programs<br>
input creation ...<br>
<br>
Michal<br>
________________________________<br>
From: users <<a href="mailto:users-bounces@lists.quantum-espresso.org" target="_blank">users-bounces@lists.quantum-espresso.org</a>> on behalf of Asma Khatun <<a href="mailto:asma@tch.hstu.ac.bd" target="_blank">asma@tch.hstu.ac.bd</a>><br>
Sent: Friday, April 25, 2025 6:54:27 PM<br>
To: <a href="mailto:users@lists.quantum-espresso.org" target="_blank">users@lists.quantum-espresso.org</a><br>
Subject: [SPAM] [QE-users] Incorrect Visualization of Ba?ZrS? Body-Centered Tetragonal Structure in XCrySDen from Quantum ESPRESSO Input<br>
<br>
<br>
Dear Quantum ESPRESSO Community,<br>
<br>
I am trying to model the Ba?ZrS? compound with a body-centered tetragonal (bct) structure using Quantum ESPRESSO. However, when I visualize the structure using XCrySDen from my SCF input file (<a href="http://ba2zrs4.scf.in" rel="noreferrer" target="_blank">ba2zrs4.scf.in</a><<a href="http://ba2zrs4.scf.in" rel="noreferrer" target="_blank">http://ba2zrs4.scf.in</a>>), the crystal does not appear as expected. The atomic positions and symmetry seem incorrect even before running any calculation.<br>
To help identify the issue, I have attached:<br>
<br>
a) <a href="http://ba2zrs4.scf.in" rel="noreferrer" target="_blank">ba2zrs4.scf.in</a><<a href="http://ba2zrs4.scf.in" rel="noreferrer" target="_blank">http://ba2zrs4.scf.in</a>>: My input file used for visualization.<br>
<br>
b) unitcell.png: The structure as visualized in XCrySDen.<br>
<br>
Could anyone kindly review my input file and advise if there is any mistake in the lattice vectors, atomic positions, or crystal symmetry definition? Any help would be sincerely appreciated.<br>
<br>
Thank you very much for your time and support.<br>
<br>
Best regards,<br>
Mst. Asma Khatun<br>
PhD Fellow, Department of Physics<br>
University of Rajshahi, Bangladesh<br>
<br>
whatsapp: +8801723674001<br>
<br>
e-mail: <a href="mailto:asma@tch.hstu.ac.bd" target="_blank">asma@tch.hstu.ac.bd</a><mailto:<a href="mailto:asma@tch.hstu.ac.bd" target="_blank">asma@tch.hstu.ac.bd</a>><br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <<a href="http://lists.quantum-espresso.org/pipermail/users/attachments/20250425/82cc8f8e/attachment-0001.html" rel="noreferrer" target="_blank">http://lists.quantum-espresso.org/pipermail/users/attachments/20250425/82cc8f8e/attachment-0001.html</a>><br>
<br>
------------------------------<br>
<br>
Message: 3<br>
Date: Sat, 26 Apr 2025 13:32:46 +0530<br>
From: Barsha Pal <<a href="mailto:barsha.pal@aus.ac.in" target="_blank">barsha.pal@aus.ac.in</a>><br>
To: <a href="mailto:users@lists.quantum-espresso.org" target="_blank">users@lists.quantum-espresso.org</a><br>
Subject: [QE-users] bfgs failed after 25 scf cycles and 23 bfgs steps,<br>
convergence not achieved<br>
Message-ID:<br>
<<a href="mailto:CAE7kQ%2BxT%2BXPHU2hWUzucAg2pU3wZgZidc-z-A0fAHCn996X9wQ@mail.gmail.com" target="_blank">CAE7kQ+xT+XPHU2hWUzucAg2pU3wZgZidc-z-A0fAHCn996X9wQ@mail.gmail.com</a>><br>
Content-Type: text/plain; charset="utf-8"<br>
<br>
The error I have got is:<br>
<br>
Message from routine bfgs:<br>
history already reset at previous step: exiting<br>
Energy error = 3.4E-09 Ry<br>
Gradient error = 4.0E-05 Ry/Bohr<br>
Cell gradient error = 2.8E-03 kbar<br>
<br>
bfgs failed after 25 scf cycles and 23 bfgs steps, convergence not<br>
achieved<br>
(criteria: energy < 1.0E-06 Ry, force < 1.0E-05 Ry/Bohr, cell <<br>
5.0E-03 kbar)<br>
<br>
End of BFGS Geometry Optimization<br>
<br>
Final enthalpy = - ** Ry<br>
<br>
File ./vc/g.bfgs deleted, as requested<br>
Begin final coordinates<br>
new unit-cell volume = ** a.u.^3 ( **Ang^3 )<br>
density = ** g/cm^3<br>
<br>
CELL_PARAMETERS (alat= )<br>
<br>
<br>
ATOMIC_POSITIONS (crystal)<br>
<br>
End final coordinates<br>
---------<br>
<br>
I get this message for a spin orbit coupling calculations, I had used all<br>
relativistic pseudos. Input file description is as follows:<br>
<br>
<br>
&CONTROL<br>
calculation = "vc-relax"<br>
etot_conv_thr = 1.00e-06<br>
forc_conv_thr = 1.00e-05<br>
nstep = 100<br>
outdir = "./vc"<br>
prefix = "g"<br>
pseudo_dir = "./"<br>
verbosity = "high"<br>
/<br>
<br>
&SYSTEM<br>
celldm(1) =<br>
degauss = 0.000735<br>
ecutrho = 440<br>
ecutwfc = 110<br>
ibrav = 0<br>
nat =<br>
ntyp = 2<br>
occupations = "smearing"<br>
smearing = "gaussian"<br>
starting_magnetization(1) = 0.00000<br>
starting_magnetization(2) = 1.00000<br>
noncolin = .true.<br>
lspinorb = .true.<br>
/<br>
<br>
&ELECTRONS<br>
conv_thr = 1.00000e-08<br>
electron_maxstep = 200<br>
mixing_beta = 7.00000e-01<br>
/<br>
<br>
&IONS<br>
ion_dynamics = "bfgs"<br>
/<br>
<br>
&CELL<br>
cell_dofree = "2Dxy"<br>
cell_dynamics = "bfgs"<br>
press = 0.00000e+00<br>
press_conv_thr = 5.00000e-03<br>
/<br>
<br>
K_POINTS {automatic}<br>
8 8 1 0 0 0<br>
<br>
CELL_PARAMETERS {alat}<br>
<br>
<br>
ATOMIC_SPECIES<br>
A ******* A.upf !using norm-conserving and full relativistic<br>
with pbe<br>
B ******* B.upf !using norm-conserving and full relativistic<br>
with pbe<br>
<br>
ATOMIC_POSITIONS {crystal}<br>
<br>
<br>
What should I do next for optimization to be converged?<br>
Also, please provide some guidelines for applying the magnetization value.<br>
How much should it be used?<br>
If I have two types of atoms and put zero for starting_magnetization for<br>
both, how is it taken? If magnetization is given for one type of atom then<br>
how is it calculated? And if startine_magnetization is given for both, then<br>
what does the code consider?<br>
<br>
<br>
Thank you.<br>
Ashley Cooper<br>
Phd<br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <<a href="http://lists.quantum-espresso.org/pipermail/users/attachments/20250426/eaaed2d7/attachment-0001.html" rel="noreferrer" target="_blank">http://lists.quantum-espresso.org/pipermail/users/attachments/20250426/eaaed2d7/attachment-0001.html</a>><br>
<br>
------------------------------<br>
<br>
Message: 4<br>
Date: Sat, 26 Apr 2025 10:27:54 +0200<br>
From: Lorenzo Paulatto <<a href="mailto:lorenzo.paulatto@cnrs.fr" target="_blank">lorenzo.paulatto@cnrs.fr</a>><br>
To: <<a href="mailto:users@lists.quantum-espresso.org" target="_blank">users@lists.quantum-espresso.org</a>><br>
Subject: Re: [QE-users] bfgs failed after 25 scf cycles and 23 bfgs<br>
steps, convergence not achieved<br>
Message-ID: <<a href="mailto:01e45f58-dd29-4b2b-94d5-833dc37403d1@cnrs.fr" target="_blank">01e45f58-dd29-4b2b-94d5-833dc37403d1@cnrs.fr</a>><br>
Content-Type: text/plain; charset="utf-8"; Format="flowed"<br>
<br>
> ? ? cell_dofree ? ?= "2Dxy"<br>
<br>
Hello Barsha,<br>
<br>
you give very little detail about your system, but I would guess that <br>
the algorithm fails because the cell constraint does not allow the cell <br>
pressure to get under the threshold, you may need more vacuum or less <br>
tight constraint.<br>
<br>
/If I have two types of atoms?and put zero for starting_magnetization <br>
for both, how is it taken? If magnetization is given for one type of <br>
atom then how is it calculated? And if startine_magnetization is given <br>
for both, then what does the code consider?/<br>
<br>
Because you have set lspinorb and noncolin=.true., the calculation will <br>
include spin-orbit coupling (if at least one of your pseudopotentials is <br>
full relativistic), but there will be no net magnetization as you are <br>
not breaking the spin symmetry in the initial conditions.<br>
hth<br>
<br>
<br>
-- <br>
Dr. Lorenzo Paulatto<br>
IdR @ IMPMC - CNRS UMR 7590 & Sorbonne Universit?<br>
phone: +33 (0)1 442 79822 / telegram: lpaulatto<br>
<a href="http://www.impmc.upmc.fr/~paulatto/" rel="noreferrer" target="_blank">http://www.impmc.upmc.fr/~paulatto/</a> - <a href="https://anharmonic.github.io/23-24/423" rel="noreferrer" target="_blank">https://anharmonic.github.io/<br>
23-24/423</a> B115, 4 place Jussieu 75252 Paris CX 05<br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <<a href="http://lists.quantum-espresso.org/pipermail/users/attachments/20250426/b527a6ab/attachment-0001.html" rel="noreferrer" target="_blank">http://lists.quantum-espresso.org/pipermail/users/attachments/20250426/b527a6ab/attachment-0001.html</a>><br>
<br>
------------------------------<br>
<br>
Message: 5<br>
Date: Sat, 26 Apr 2025 10:31:02 +0200<br>
From: Paolo Giannozzi <<a href="mailto:paolo.giannozzi@uniud.it" target="_blank">paolo.giannozzi@uniud.it</a>><br>
To: Quantum ESPRESSO users Forum <<a href="mailto:users@lists.quantum-espresso.org" target="_blank">users@lists.quantum-espresso.org</a>>,<br>
Barsha Pal <<a href="mailto:barsha.pal@aus.ac.in" target="_blank">barsha.pal@aus.ac.in</a>><br>
Subject: Re: [QE-users] bfgs failed after 25 scf cycles and 23 bfgs<br>
steps, convergence not achieved<br>
Message-ID: <<a href="mailto:c96cece2-4179-4802-b19c-007e1757d797@uniud.it" target="_blank">c96cece2-4179-4802-b19c-007e1757d797@uniud.it</a>><br>
Content-Type: text/plain; charset=UTF-8; format=flowed<br>
<br>
On 26/04/2025 10:02, Barsha Pal wrote:<br>
<br>
> ? ? ?Energy error ? ? ? ? ? ?= ? ? ?3.4E-09 Ry<br>
> ? ? ?Gradient error ? ? ? ? ?= ? ? ?4.0E-05 Ry/Bohr<br>
> ? ? ?Cell gradient error ? ? = ? ? ?2.8E-03 kbar<br>
<br>
these values are already exceedingly small. There is no point in pushing <br>
the optimization further in my opinion.<br>
<br>
Paolo<br>
-- <br>
Paolo Giannozzi, Dip. Scienze Matematiche Informatiche e Fisiche,<br>
Univ. Udine, via delle Scienze 206, 33100 Udine Italy, +39-0432-558216<br>
<br>
<br>
<br>
------------------------------<br>
<br>
Subject: Digest Footer<br>
<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><br>
<br>
------------------------------<br>
<br>
End of users Digest, Vol 213, Issue 23<br>
**************************************<br>
</blockquote></div>