
<html dir="ltr">

<head>

<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">

<style id="owaParaStyle" type="text/css">P {margin-top:0;margin-bottom:0;}</style>

</head>

<body ocsi="0" fpstyle="1">

<div style="direction: ltr;font-family: Tahoma;color: #000000;font-size: 10pt;">Dear All,<br>

<br>

I am doing a scf calculation for a particular non-collinear spin configuration. A am using LDA&#43;U method. When I use the following case.inorb and case.indmc files:<br>

case.inorb <br>

&nbsp;1&nbsp; 3&nbsp; 0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; nmod, natorb, ipr<br>

PRATT&nbsp; 1.0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; BROYD/PRATT, mixing<br>

&nbsp; 1 1 3<br>

&nbsp; 2 1 3<br>

&nbsp; 3 1 2<br>

&nbsp; 1&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; nsic 0..AMF, 1..SIC, 2..HFM<br>

&nbsp;&nbsp; 0.65 0.07&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; U J (Ry)<br>

&nbsp;&nbsp; 0.65 0.07<br>

&nbsp;&nbsp; 0.45 0.07<br>

case.indmc<br>

-9.&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Emin cutoff energy<br>

&nbsp;3&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; number of atoms for which density matrix is calculated<br>

&nbsp;1&nbsp; 1&nbsp; 3&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; index of 1st atom, number of L's, L1<br>

&nbsp;2&nbsp; 1&nbsp; 3&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; dtto for 2nd atom, repeat NATOM times<br>

&nbsp;3&nbsp; 1&nbsp; 2&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; index of 1st atom, number of L's, L1<br>

&nbsp;0 0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; r-index, (l,s)index<br>

<br>

the program executes properly without any error. However when I use the other set of case.inorb and case.indmc files:<br>

case.inorb<br>

&nbsp; 1&nbsp; 3&nbsp; 0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; nmod, natorb, ipr<br>

PRATT&nbsp; 1.0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; BROYD/PRATT, mixing<br>

&nbsp; 1 2 3 2<br>

&nbsp; 2 2 3 2<br>

&nbsp; 3 1 2<br>

&nbsp; 1&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; nsic 0..AMF, 1..SIC, 2..HFM<br>

&nbsp;&nbsp; 0.65 0.07<br>

&nbsp;&nbsp; 0.45 0.07&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; U J (Ry)<br>

&nbsp;&nbsp; 0.65 0.07<br>

&nbsp;&nbsp; 0.45 0.07<br>

&nbsp;&nbsp; 0.45 0.07<br>

case.indmc<br>

-9.&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Emin cutoff energy<br>

&nbsp;3&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; number of atoms for which density matrix is calculated<br>

&nbsp;1&nbsp; 2&nbsp; 3&nbsp; 2 &nbsp;&nbsp;&nbsp;&nbsp; index of 1st atom, number of L's, L1<br>

&nbsp;2&nbsp; 2&nbsp; 3&nbsp; 2 &nbsp;&nbsp;&nbsp;&nbsp; dtto for 2nd atom, repeat NATOM times<br>

&nbsp;3&nbsp; 1&nbsp; 2&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; index of 1st atom, number of L's, L1<br>

&nbsp;0 0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; r-index, (l,s)index<br>

the program stops by giving the error: error in Vorb. It is not writing case.vorbup/dn files. In the next set I just use U at the d states as well as at the f states for atom 1 and 2. Am I giving any wrong input? Help me to solve this problem.<br>

<br>

Waiting for suggestions.<br>

<br>

Best,<br>

Hena<br>

<br>

<br>

</div>

</body>

</html>