<div dir="ltr"><div>  Dear Wien2k community and developers,<div dir="auto"><br></div><div dir="auto">I am encountering an issue with LMBJ calculations for 2D BaSe monolayer </div><div dir="auto">using Wien2k 23.2. The nlvdw step crashes immediately with "Error in NLVDW", </div><div dir="auto">while the identical workflow works perfectly for CaSe and SrSe 2D monolayers.</div><div dir="auto"><br></div><div dir="auto">== System Details ==</div><div dir="auto"><br></div><div dir="auto">Material: 2D BaSe monolayer (PbO-type structure)</div><div dir="auto">Space group: 123 (P4/mmm)</div><div dir="auto">Lattice: a = b = 8.587093 Bohr, c = 42.038689 Bohr (~22 Å vacuum)</div><div dir="auto">Atoms: Ba at (0, 0, 0.5), Se at (0.5, 0.5, 0.5)</div><div dir="auto">RMT: Ba = 2.5, Se = 2.3</div><div dir="auto">Spin-unpolarized calculation</div><div dir="auto"><br></div><div dir="auto">== Key Parameters ==</div><div dir="auto"><br></div><div dir="auto">case.in1:</div><div dir="auto">  RKmax = 7.0, MAX L IN WF = 10, V-NMT = 4</div><div dir="auto">  Local orbitals included for Ba 4d (-6.21 Ry), 5s (-2.17 Ry), </div><div dir="auto">  5p (-1.03 Ry), and Se 3d (-3.57 Ry)</div><div dir="auto">  Energy cutoff used in lstart: -8.0 Ry</div><div dir="auto"><br></div><div dir="auto">case.in2:</div><div dir="auto">  EMIN = -14.0 Ry, NE = 36</div><div dir="auto">  Lattice harmonics: max L = 6</div><div dir="auto">  GMAX = 12.00</div><div dir="auto"><br></div><div dir="auto">case.innlvdw (identical to working CaSe):</div><div dir="auto">  kernel type: 1</div><div dir="auto">  kernel parameter: -0.8491</div><div dir="auto">  GMAX: 25</div><div dir="auto">  rhoc: 0.3</div><div dir="auto">  potential: T</div><div dir="auto">  GMAXpot: 10</div><div dir="auto">  sigma: 3.78 bohr</div><div dir="auto"><br></div><div dir="auto">== What Works ==</div><div dir="auto"><br></div><div dir="auto">PBE SCF: Converges cleanly in 11 iterations.</div><div dir="auto">  - :NEC = 90.00000 / 90.00000000 (zero deviation)</div><div dir="auto">  - :ENE change < 1e-7 Ry between last iterations</div><div dir="auto">  - :DIS = 0.0000276 (well below threshold)</div><div dir="auto">  - :GAP (global) = 0.148086 Ry = 2.015 eV</div><div dir="auto">  - No warnings</div><div dir="auto"><br></div><div dir="auto">CaSe and SrSe 2D monolayers: LMBJ works without any issues</div><div dir="auto">(same workflow, same parameter ranges, same Wien2k version).</div><div dir="auto"><br></div><div dir="auto">== The Problem ==</div><div dir="auto"><br></div><div dir="auto">After:</div><div dir="auto">  save_lapw BaSe_2D_pbe</div><div dir="auto">  init_mbj_lapw (chose: 2 for lmBJ, then 0 for default Koller 2012)</div><div dir="auto">  run_lapw -i 1 -NI -p</div><div dir="auto">  save_lapw PBE</div><div dir="auto">  init_mbj_lapw (chose: 2 for lmBJ, then 0)</div><div dir="auto">  run_lapw -lmbj -i 100 -ec 0.0001 -cc 0.0001 -p</div><div dir="auto"><br></div><div dir="auto">The dayfile shows:</div><div dir="auto">  cycle 1 (...) (100/99 to go)</div><div dir="auto">  > lapw0 -grr 5.144u 0.121s 0:01.71 307.6%</div><div dir="auto">  > nlvdw -lmbj 0.007u 0.008s 0:00.01 0.0% <-- crashes immediately</div><div dir="auto">  > stop error</div><div dir="auto"><br></div><div dir="auto">The nlvdw.error file contains only: "Error in NLVDW"</div><div dir="auto"><br></div><div dir="auto">== Diagnostic Information ==</div><div dir="auto"><br></div><div dir="auto">- nlvdw binary exists at /opt/WIEN_23/nlvdw</div><div dir="auto">- All required input files exist with correct sizes:</div><div dir="auto">  BaSe_2D.r2v: 565,653 bytes</div><div dir="auto">  BaSe_2D.r2vdn: 565,653 bytes (copy of r2v for spin-unpolarized)</div><div dir="auto">  BaSe_2D.tausum: 598,913 bytes</div><div dir="auto">  BaSe_2D.innlvdw: parameters identical to working CaSe/SrSe</div><div dir="auto">- Removed all empty/stale tauup, taudn files</div><div dir="auto">- Tested both serial and parallel runs - same error</div><div dir="auto">- case.in0 has XC_LMBJ correctly set</div><div dir="auto">- 2*lnsmax (=8) >= lmax in case.in2 (=6), so this is not the </div><div dir="auto">  lmax-in-in2 vs. lnsmax issue from the older nlvdw thread (msg16382, 2017)</div><div dir="auto"><br></div><div dir="auto">== What I Have Tried ==</div><div dir="auto"><br></div><div dir="auto">1. Compared every parameter with working CaSe directory - identical except </div><div dir="auto">   for atom-specific values (EMIN, NE, GMAX which is 12 vs 16).</div><div dir="auto">2. Tried with -p and without -p (parallel/serial)</div><div dir="auto">3. Verified all files have non-zero sizes</div><div dir="auto">4. Manually copied r2v_grr -> r2v (Wien2k 23.2 uses _grr suffix)</div><div dir="auto">5. Restored tausum from saved PBE state</div><div dir="auto"><br></div><div dir="auto">== Question ==</div><div dir="auto"><br></div><div dir="auto">Could this be a numerical issue specific to heavy elements with deep </div><div dir="auto">semicore d-states (Ba 4d at -6.21 Ry)? The fact that CaSe/SrSe work </div><div dir="auto">but BaSe doesn't, with the same workflow and parameter ranges, </div><div dir="auto">suggests a Ba-specific limitation in nlvdw -lmbj.</div><div dir="auto"><br></div><div dir="auto">Is there any known workaround for LMBJ with heavy alkaline earth </div><div dir="auto">elements like Ba? Or should I expect regular mBJ to give comparable </div><div dir="auto">results given my 22 Å vacuum thickness?</div><div dir="auto"><br></div><div dir="auto">I am happy to provide case.struct, case.in0, case.in1, case.in2, </div><div dir="auto">case.innlvdw or any other files if needed for diagnosis.</div><div dir="auto"><br></div><div dir="auto">Thank you for any guidance!</div><div dir="auto"><br></div></div><div><div dir="ltr" class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><p class="MsoNormal" style="margin-bottom:0.0001pt;text-align:justify;line-height:14.65pt"><b><i><span style="font-family:Georgia,serif"><font color="#20124d">Thanks and Regards</font></span></i></b></p><p class="MsoNormal" style="margin-bottom:0.0001pt;text-align:justify;line-height:14.65pt"><b><i><span style="font-family:Georgia,serif"><font color="#4c1130">Dr. Karin</font><font color="#741b47">a</font><font color="#4c1130"> </font></span></i></b></p><p class="MsoNormal" style="margin-bottom:0.0001pt;line-height:normal"><b style="color:rgb(102,0,0)"><i><span style="font-family:Georgia,serif">| </span></i></b><b><i><span style="font-family:Georgia,serif"><a href="https://scholar.google.com/citations?user=2_K6k2MAAAAJ&hl=en" target="_blank"><font color="#660000">Google Scholar</font></a><font color="#660000"> |</font></span></i></b><b style="color:rgb(102,0,0)"><i><span style="font-family:Georgia,serif"> </span></i></b></p></div><div style="margin:0px;padding:0px 0px 20px;width:1142px;font-family:Roboto,RobotoDraft,Helvetica,Arial,sans-serif;font-size:medium"><div><div style="font-size:12.8px;margin:8px 0px 0px;padding:0px"><div><div dir="ltr"><div><div dir="ltr"><span><div><br><span></span><span></span></div><b><font color="#0b5394"></font></b></span></div></div><div></div></div><div></div></div></div><div></div></div></div></div></div></div></div></div></div>