<div dir="ltr"><div dir="ltr">On Wed, Mar 10, 2021 at 5:10 PM Kevin Kayang <<a href="mailto:kwkayang@crimson.ua.edu">kwkayang@crimson.ua.edu</a>> wrote:<br></div><div class="gmail_quote"><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"> ibrav     = 0,<br><div> celldm(1)= 1.88972687,<br></div></div></blockquote><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div>[...] <br></div></div></blockquote><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div>CELL_PARAMETERS cubic <br>   20.000000000   0.0000000000   0.0000000000<br>   0.0000000000   20.000000000   0.0000000000<br>   0.0000000000   0.0000000000   20.000000000<br></div></div></blockquote><div><br></div><div>this is an exceptionally obscure way to provide input data. A better one: ibrav=1, a=20.0</div><br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div> nosym = .FALSE.<br></div></div></blockquote><div><br></div><div>why? do you know what it means?<br></div><div> <br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div> input_dft = vdw-df2<br></div></div></blockquote><div><br></div><div>correct syntax: input_dft='vdw-df2'</div><div> <br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div> one_atom_occupations = .TRUE.,<br></div></div></blockquote><div><br></div><div>why? do you know what it means?<br></div><div><br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div> starting_magnetization(1) = 1,<br></div></div></blockquote><div><br></div><div>this means 'all electrons with spin up'. You may use a smaller value to break the symmetry<br></div><div><br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div> starting_magnetization(2) = 1,<br></div></div></blockquote><div><br></div><div>this is wrong: there is just one kind of atom, not two, and may lead to unpredictable behavior.</div><div><br></div>The energy vs distance curve of a dimer is a simple but tricky problem. What you observe might be the result of the crossing of two different curves, one with no spin and the other with spin polarization.</div><div class="gmail_quote"><br></div><div class="gmail_quote">Paolo<br clear="all"></div>-- <br><div dir="ltr" class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div>Paolo Giannozzi, Dip. Scienze Matematiche Informatiche e Fisiche,<br>Univ. Udine, via delle Scienze 206, 33100 Udine, Italy<br>Phone +39-0432-558216, fax +39-0432-558222<br><br></div></div></div></div></div></div>