<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=utf-8">
<style type="text/css" style="display:none;"><!-- P {margin-top:0;margin-bottom:0;} --></style>
</head>
<body dir="ltr">
<div id="divtagdefaultwrapper" style="font-size:12pt;color:#000000;font-family:Calibri,Arial,Helvetica,sans-serif;" dir="ltr">
<p><span style="font-size: 12pt;">I would like to perform a spin-polarized md-simulation for a system containing among other things Sulfur atoms </span><span style="font-size: 12pt;">using the DFTB+ package. </span><br>
</p>
<p><span style="font-size: 12pt;"><br>
</span></p>
<p><span style="font-size: 12pt;">I have a few doubts that I would very much appreciate if you could help me resolve:</span></p>
<p><span style="font-size: 12pt;"><br>
</span></p>
<p><span style="font-size: 12pt;">1) <span style="font-family: Calibri, Arial, Helvetica, sans-serif, "Apple Color Emoji", "Segoe UI Emoji", NotoColorEmoji, "Segoe UI Symbol", "Android Emoji", EmojiSymbols; font-size: 16px;">The basis set for Sulfur atom in
 DFTB3 includes d-orbitals. I have read in the literature (</p>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif, "Apple Color Emoji", "Segoe UI Emoji", NotoColorEmoji, "Segoe UI Symbol", "Android Emoji", EmojiSymbols; font-size: 16px;">
"Zoltan Bodrog* and Balint Aradi Phys. Status Solidi B 249, No. 2, 259–269 (2012) / DOI 10.1002/pssb.201100524<span style="font-size: 12pt;">) the following: "</span><span style="font-size: 12pt;">The simplistic, </span><span style="font-size: 12pt;">‘semidiagonal’
 U matrix is a good approximation for </span>elements up to the third period, but it is significantly wrong</div>
</span></span><span style="font-size: 12pt;"><span style="font-family: Calibri, Arial, Helvetica, sans-serif, "Apple Color Emoji", "Segoe UI Emoji", NotoColorEmoji, "Segoe UI Symbol", "Android Emoji", EmojiSymbols; font-size: 16px;">
<div>above (subshell hardness errors going well beyond 10%, see<br>
Table 1) if d orbitals are present in the atomic basis."</div>
<div><br>
</div>
<div>I understand that S is not third period but certainly polarizable. Is it OK to therefore use an atomic based Hubbard matrix in this particular case for open shell situations for example including an excess electron? In other words, since <span style="font-size: 12pt;">I
 would like to use DFTB3  and the 3ob-3-1 S-K parameters and orbital resolved Hubbard derivatives are not available for most of the atoms </span><span style="font-size: 12pt;">is it reasonable to use the atomically-resolved methodology with the option </span><span style="font-size: 12pt;">ShellResolvedSpin=Yes
 in DFTB+</span><span style="font-size: 12pt;">? </span><span style="font-size: 12pt;">I notice that the spin parameters exist for all the orbitals in my system.</span></div>
</span></span>
<p></p>
<p><br>
</p>
<p>2) My second question is related to Wpp vs Wdd in this particular case. In </p>
<p><span>Patrick Melix, Augusto Faria Oliveira, Robert Rüger, Thomas Heine, <span>Theor Chem Acc (2016) 135:232</span></span><span style="font-size: 12pt;"> spin polarized SCC-DFTB using 3ob-3-1 is used. Only one spin parameter is provided per atom. If I look
 at this spin parameter in the case of the S atom it appears to be that corresponding to Wpp. However, 3ob-3-1 includes d orbitals for S.</span></p>
<p><span style="font-size: 12pt;"><br>
</span></p>
<p><span style="font-size: 12pt;">How should I understand this (beyond the fact that d orbitals are in principle empty but are there for polarization effects)? How should one decide about using Wpp or Wdd in this case, particularly considering that the condensed
 phase system includes an excess electron (not necessarily localized on S but who knows)?</span><span style="font-size: 12pt;"></span></p>
<p><br>
</p>
<p>Sorry that my email is verbose, I am trying to understand as much as I can how to properly use DFTB without trivial mistakes. I have found almost no calculations of open shell systems in the literature and I would like to know if this is for some good reason.</p>
<p><br>
</p>
<p>Best Regards,</p>
<p>Kamal</p>
</div>
</body>
</html>