[Wien] [EXTERNAL] Re: YCo5 magnetic anisotropy calculations with GGA+U on Co

Samolyuk, German D. samolyukgd at ornl.gov
Fri Dec 1 19:19:21 CET 2023


>This depends on how you start the EIGN calculation. From a scalar
>relativistic GGA+U calculation or just from GGA ??

I started from GGA+U scalar relativistic

>Not a big surprise. During the scf cycle the potential changes.
>Remember, by adding a constant potential, the EIGEN values would change
>by shift*NE, but the potential energy would also shift in opposite
>direction by -shift*NE, so Etot remains constant, but EIGEN not. (NE ..
>number of electrons)

Is it possible to output the shift value? It will allow to shift ebnd by the same value in both magnetization direction calculations. I don't expect to have the same MAE energy values, but if EBND gives main contribution to the value obtained from total energy difference it will be meaningful  to discuss.

Thank you detailed explanation. I really appreciate it

German

Dr. German D Samolyuk
Materials Theory Group
Materials Science & Technology Division
Oak Ridge National Laboratory
Post Office Box 2008
Oak Ridge, TN 37831-6138
(865) 241-5394
(865) 241-3650 (FAX)
________________________________
From: Wien <wien-bounces at zeus.theochem.tuwien.ac.at> on behalf of Peter Blaha <peter.blaha at tuwien.ac.at>
Sent: Friday, December 1, 2023 2:27 AM
To: A Mailing list for WIEN2k users <wien at zeus.theochem.tuwien.ac.at>
Subject: [Wien] [EXTERNAL] Re: YCo5 magnetic anisotropy calculations with GGA+U on Co


Initially I executed calculation in the "spirit of force theorem" with
just one iteration and obtained well converged K value calculated from
EIGEN. However reviewer mentioned that force theorem is not obeyed in
case of +U calculation - corrections are not in second order. He is
right, the K value obtained from full scf calculations and total energy
difference is ~40% higher. For just GGA both approaches gives the same
result.
-------------------------------------------

This depends on how you start the EIGN calculation. From a scalar
relativistic GGA+U calculation or just from GGA ??

And in addition it is probably again a problem of symmetry. Reliable
values can only be obtained with identical symmetry in the scalarel and
SO calculation. So you have to perform the scalarel. calc with the same
struct file (symmetry operations) and case.in2c (LM combinations) as the
2 SO calculations.

It is NOT only a matter of k-points and cannot be "fixed by x kgen -fbz
---------------------------------------------------
---------------------------------------------------
What surprises me is the non-monotonic dependence of K value as a number
of k-points. At the same time K calculated from the total energy
converged. While, the  band structure energy is part of the total energy.
------------------------------
Not a big surprise. During the scf cycle the potential changes.
Remember, by adding a constant potential, the EIGEN values would change
by shift*NE, but the potential energy would also shift in opposite
direction by -shift*NE, so Etot remains constant, but EIGEN not. (NE ..
number of electrons)

--
--------------------------------------------------------------------------
Peter BLAHA, Inst.f. Materials Chemistry, TU Vienna, A-1060 Vienna
Phone: +43-1-58801-165300
Email: peter.blaha at tuwien.ac.at    WIEN2k: https://urldefense.us/v2/url?u=http-3A__www.wien2k.at&d=DwICAg&c=v4IIwRuZAmwupIjowmMWUmLasxPEgYsgNI-O7C4ViYc&r=lTD7jQRfMViWEsN8TZ1wLmkMhVe4MCRH76GADmGBpD4&m=2rJoswOBQ92PKOEckFLkl8FPEmeAYw5-vJ6W5FNVUsIRT9dNeM3rsFcJOTPpNFJb&s=8Bc_0Dag-FLMvY2FdwzRLEjWbmXdtmAew2_77xLIJkg&e=
WWW:   https://urldefense.us/v2/url?u=http-3A__www.imc.tuwien.ac.at&d=DwICAg&c=v4IIwRuZAmwupIjowmMWUmLasxPEgYsgNI-O7C4ViYc&r=lTD7jQRfMViWEsN8TZ1wLmkMhVe4MCRH76GADmGBpD4&m=2rJoswOBQ92PKOEckFLkl8FPEmeAYw5-vJ6W5FNVUsIRT9dNeM3rsFcJOTPpNFJb&s=DzH1tAa_e6jpwNi0n5EBke6n2M4jAHbZpKYBS3ISpC8&e=
-------------------------------------------------------------------------
_______________________________________________
Wien mailing list
Wien at zeus.theochem.tuwien.ac.at
https://urldefense.us/v2/url?u=http-3A__zeus.theochem.tuwien.ac.at_mailman_listinfo_wien&d=DwICAg&c=v4IIwRuZAmwupIjowmMWUmLasxPEgYsgNI-O7C4ViYc&r=lTD7jQRfMViWEsN8TZ1wLmkMhVe4MCRH76GADmGBpD4&m=2rJoswOBQ92PKOEckFLkl8FPEmeAYw5-vJ6W5FNVUsIRT9dNeM3rsFcJOTPpNFJb&s=3kFudcBlmqk8rCLFIS8Q6qbM3pxWK9eXVxgLnGDlinE&e=
SEARCH the MAILING-LIST at:  https://urldefense.us/v2/url?u=http-3A__www.mail-2Darchive.com_wien-40zeus.theochem.tuwien.ac.at_index.html&d=DwICAg&c=v4IIwRuZAmwupIjowmMWUmLasxPEgYsgNI-O7C4ViYc&r=lTD7jQRfMViWEsN8TZ1wLmkMhVe4MCRH76GADmGBpD4&m=2rJoswOBQ92PKOEckFLkl8FPEmeAYw5-vJ6W5FNVUsIRT9dNeM3rsFcJOTPpNFJb&s=sqNi0wE0N-5TvU2CVqKzeYBfBb4jjrcIwejWUItyZYw&e=
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://zeus.theochem.tuwien.ac.at/pipermail/wien/attachments/20231201/89c1d66e/attachment.htm>


More information about the Wien mailing list