> However, suppose the Hubbard atom is Ce, "[Xe]6s2 4f1 5d1" with two localized 
> incomplete orbitals of Ce-4f and Ce-5d. Then how can one specify the 
> starting_ns_eigenvalue for both of them with l=3 and l=2?


This is currently not possible to do (not implemented). starting_ns_eigenvalue 
changes only the eigenvalues of the occupation matrix corresponding to the 
first Hubbard channel.


Note that applying +U to two Hubbard channels of the same atomic type is a 
relatively new feature, and hence many things are not supported for this case 
(like the issue mentioned about).


> It seems that the second localized orbital is treated as "background part" 
> and accidentally(?!) the eigenvectors of its nod-degenerate eigenvalues are 
> set to 1.0 !


The second Hubbard channel can be called "background", but it is nothing but 
just the second channel. The populations of the first and second Hubbard 
channels are initialized in the code (PW/src/init_ns.f90) and I do not really 
understand what is "accidental". If you want to have more precise answers, you 
should be more specific.


HTH


Iurii


--
Dr. Iurii TIMROV
Senior Research Scientist
Theory and Simulation of Materials (THEOS)
Swiss Federal Institute of Technology Lausanne (EPFL)
CH-1015 Lausanne, Switzerland
+41 21 69 34 881
http://people.epfl.ch/265334
________________________________
From: users <users-boun...@lists.quantum-espresso.org> on behalf of Mahmoud 
Payami Shabestari via users <users@lists.quantum-espresso.org>
Sent: Friday, August 12, 2022 7:48:32 AM
To: Mahmoud Payami Shabestari; Quantum ESPRESSO users Forum
Subject: Re: [QE-users] syntax of using starting_ns_eigenvalue for two 
different localized orbitals of the same atom

It seems that the second localized orbital is treated as "background part" and 
accidentally(?!) the eigenvectors of its nod-degenerate eigenvalues are set to 
1.0 !

mahmoud




From: "Mahmoud Payami Shabestari" <mpay...@aeoi.org.ir>
To: "Quantum ESPRESSO users Forum" <users@lists.quantum-espresso.org>
Date: Fri, 12 Aug 2022 07:32:05 +0430
Subject: Re: [QE-users] syntax of using starting_ns_eigenvalue for two 
different localized orbitals of the same atom

Dear Stefano,

Hi.
Thank you for your detailed explanation for the Fe case.
However, suppose the Hubbard atom is Ce, "[Xe]6s2 4f1 5d1" with two localized 
incomplete orbitals of Ce-4f and Ce-5d. Then how can one specify the 
starting_ns_eigenvalue for both of them with l=3 and l=2?
Thank you in advance.

Best regards,
mahmoud


From: Stefano de Gironcoli <degir...@sissa.it>
To: Mahmoud Payami Shabestari via users <users@lists.quantum-espresso.org>
Date: Thu, 11 Aug 2022 21:19:56 +0200
Subject: Re: [QE-users] syntax of using starting_ns_eigenvalue for two 
different localized orbitals of the same atom

suppose your Hubbard atom of type 1 is an Fe ( U on 3d orbitals)..., then
starting_ns_eigenvalue(3,1,1)=1.0  overwrites the third eigenvalue of the 
occupations with spin=1 (up) to 1.0
starting_ns_eigenvalue(4,1,1)=0.3 overwrites the forth eigenvalue of the 
occupations with spin=1 (up) to 0.3
starting_ns_eigenvalue(5,2,1)=0.0 overwrites the fifth eigenvalue of the 
occupations with spin=2 (dw) to 0.0
the other elements are left unchanged.
In general at the first iteration the code computes the occupation matrix ns, 
diagonalizes it, overwrites the eigenvalues if required in input and then 
calculates the output ns matrix from the (overwritten) eigenvalue decomposition.
if you specify more starting_ns_eigenvalue lines for a given atomic type the 
code will apply all of them independently
stefano
On 11/08/22 21:00, Mahmoud Payami Shabestari via users wrote:
Dear QE-7.1 Developers,

Hi. I would like to use the starting_ns_eigenvalue()=1.0 for two different 
localized orbitals of the same atom. I do not have problem in the HUBBARD card, 
but I do not know how to use starting_ns_eigenvalue()=1.0 for different orbital 
angular momenta of the same atom.
Any comments is highly appreciated.

Best regards,
Mahmoud Payami


---------------------------------------
Mahmoud Payami
NSTRI, AEOI, Tehran, Iran

Email: mpay...@aeoi.org.ir<mailto:mpay...@aeoi.org.ir>
Phone: +98 (0)21 82066504
--------------------------------------------------------






_______________________________________________
The Quantum ESPRESSO community stands by the Ukrainian
people and expresses its concerns about the devastating
effects that the Russian military offensive has on their
country and on the free and peaceful scientific, cultural,
and economic cooperation amongst peoples
_______________________________________________
Quantum ESPRESSO is supported by MaX 
(www.max-centre.eu<http://www.max-centre.eu>)
users mailing list 
users@lists.quantum-espresso.org<mailto:users@lists.quantum-espresso.org>
https://lists.quantum-espresso.org/mailman/listinfo/users
_______________________________________________
The Quantum ESPRESSO community stands by the Ukrainian
people and expresses its concerns about the devastating
effects that the Russian military offensive has on their
country and on the free and peaceful scientific, cultural,
and economic cooperation amongst peoples
_______________________________________________
Quantum ESPRESSO is supported by MaX (www.max-centre.eu)
users mailing list users@lists.quantum-espresso.org
https://lists.quantum-espresso.org/mailman/listinfo/users

Reply via email to