<div dir="ltr">Dear Slava,<div><br></div><div>Can I ask you a question? </div><div>I want to applied DFT+U method too, but I don't kown how to measure U_MINUS_J value , I find in you case the value is only  0.1234777, So what dose 'J' mean and how to calculate?</div><div>Thank you very much.</div><div><br></div><div>Xiaomei </div><div><br>在 2019年4月24日星期三 UTC-4下午12:22:27，Vyacheslav Bryantsev写道：<blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;"><div dir="ltr"><div>Dear CP2K Group,</div><div><br></div><div>There are a couple of ways to setup antiferromagnetic calculations for 6 Cr3+ ions (3d3 state) dissolved in MgCl2-KCl melt. <br></div><div><br></div><div>One way is to use MAGNETIZATION keyword. Another way is to use the BS section (commented below)</div><div><br></div><div>&KIND Cr1<br>      ELEMENT Cr<br>      BASIS_SET DZVP-MOLOPT-SR-GTH<br>      POTENTIAL GTH-PBE-q14<br>      <b>MAGNETIZATION 4.0</b><br>      &DFT_PLUS_U<br>        L 2<br>        U_MINUS_J 0.1234777<br>        # U = 3.36 eV<br>        U_RAMPING 0.01234777<br>       &END DFT_PLUS_U</div><div> <b>!       &BS<br> !        &ALPHA<br> !         NEL +1 -1<br> !         L 2 0<br> !         N 3 4<br> !       &END ALPHA<br> !       &BETA<br> !        NEL -5 -1<br> !        L 2 0<br> !        N 3 4<br> !      &END BETA</b><br>     &END BS<br>    &END KIND<br>    &KIND Cr2<br>      ELEMENT Cr<br>      BASIS_SET DZVP-MOLOPT-SR-GTH<br>      POTENTIAL GTH-PBE-q14<br>      <b>MAGNETIZATION -4.0</b><br>      &DFT_PLUS_U<br>        L 2<br>        U_MINUS_J 0.1234777<br>        # U = 3.36 eV<br>        U_RAMPING 0.01234777<br>       &END DFT_PLUS_U<br></div><div><b>!       &BS<br>!        &ALPHA<br>!         NEL -5 -1<br>!         L 2 0<br>!         N 3 4<br>!       &END ALPHA<br>!       &BETA<br>!        NEL +1 -1<br>!        L 2 0<br>!        N 3 4<br>!      &END BETA</b></div><div></div><div>    &END KIND<br></div><div><br></div><div>The problem I have is that the total E from two setups (MAGNETIZATION and BS) can be somewhat different:</div><div>For example, the Etot of a system (mol fractions are shown) is lower when setting up calculations using BS<br></div><div>50-MgCl2-50-KCl-5-CrCl3</div><div>Etot = -8997.64290729994173 using BS section<br></div><div>Etot = -8997.62650294793639 using MAGNETIZATION keyword<br></div><div> </div><div><br></div><div>However, for a different composition, the Etot is lower using MAGNETIZATION <br></div><div><div>15-MgCl2-85-KCl-5-CrCl3</div><div>Etot = -6895.30145964661006 using BS section<br></div><div>Etot = -6895.31001188029859 using MAGNETIZATION keyword</div><div><br></div><div>The resulting spin moment on each Cr3+ in both cases are close to +3 and -3, as expected.<br></div><div><br></div><div><b>What are the general recommendations to correctly find a ground electronic state using DFT+U for transition metal ions in solution/melts?<br></b></div><div><b>Are there any additional checks one can do to make sure the ground electronic state is found?<br></b></div><div>After the ground electronic state is found, the next step will be to run AIMD simulations. <b><br></b></div><div> </div><div>The input and output files for these examples are given as attachments.</div><div><br></div><div>Thank you,</div><div>Slava</div><div><br></div></div><div><br></div><div><br></div><div><br></div></div></blockquote></div></div>