<html>
  <head>
    <meta http-equiv="Content-Type" content="text/html; charset=UTF-8">
  </head>
  <body>
    <p>Hi </p>
    <p> the starting_magnetization is specified per species not per
      site, this may be sometime not very practical, but in your case
      warrants  that the supercell is initialized in the same way as
      the  bulk one. Consider also that if the ferromagnetic solution is
      stable, the final density and magnetization are not very sensitive
      to the specific initialization, the starting_magnetization only
      breaks the symmetry between the two polarization and defines which
      is the majority spin. <br>
    </p>
    <p>Contrains instead affect the energy with an additive term that
      depends on the deviation from the constaint, and as you are going
      to do finite differences between displaced systems, where you
      might actually expect some small variation of the local
      magnetizations this could influence you results. <br>
    </p>
    <p>hope this helps - regards <br>
    </p>
    <p>Pietro<br>
    </p>
    <div class="moz-cite-prefix">On 28/07/20 00:35, Hien Vo wrote:<br>
    </div>
    <blockquote type="cite"
cite="mid:BYAPR11MB3766ED92504D4C80496EC361CC720@BYAPR11MB3766.namprd11.prod.outlook.com">
      <meta http-equiv="Content-Type" content="text/html; charset=UTF-8">
      <meta name="Generator" content="Microsoft Word 15 (filtered
        medium)">
      <style><!--
/* Font Definitions */
@font-face
        {font-family:"Cambria Math";
        panose-1:2 4 5 3 5 4 6 3 2 4;}
@font-face
        {font-family:Calibri;
        panose-1:2 15 5 2 2 2 4 3 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
        {margin:0in;
        margin-bottom:.0001pt;
        font-size:11.0pt;
        font-family:"Calibri",sans-serif;}
.MsoChpDefault
        {mso-style-type:export-only;}
@page WordSection1
        {size:8.5in 11.0in;
        margin:1.0in 1.0in 1.0in 1.0in;}
div.WordSection1
        {page:WordSection1;}
--></style>
      <div class="WordSection1">
        <p class="MsoNormal">Hi all,<o:p></o:p></p>
        <p class="MsoNormal"><o:p> </o:p></p>
        <p class="MsoNormal">I’m trying to compute phonons by the frozen
          phonon approach of a ferromagnetic system. I’m not exactly
          sure how to go about it. <o:p></o:p></p>
        <p class="MsoNormal"><o:p> </o:p></p>
        <p class="MsoNormal">I’m planning to converge a unit cell of 20
          atoms using the starting_magnetization option. Then I will use
          Phonopy to generate a set of 80-atom supercells with finite
          displacements and run pwscf on those. How should I treat
          magnetization in these supercells? I’m concerned that if I use
          starting_magnetization again, the systems might converge to
          different ‘magnetic moment per site’ values. <o:p></o:p></p>
        <p class="MsoNormal"><o:p> </o:p></p>
        <p class="MsoNormal">I’m thinking of trying the
          constrained_magnetization=’atomic’ option and constrain each
          atomic magnetic moment to the value achieved in the 20-atom
          unit cell. Is this the way to go? <o:p></o:p></p>
        <p class="MsoNormal"><o:p> </o:p></p>
        <p class="MsoNormal">Best,<o:p></o:p></p>
        <p class="MsoNormal">Hien Vo<o:p></o:p></p>
        <p class="MsoNormal">PhD Student<o:p></o:p></p>
        <p class="MsoNormal">Chemistry Department | The University of
          Chicago <o:p></o:p></p>
        <p class="MsoNormal"><o:p> </o:p></p>
        <p class="MsoNormal"><o:p> </o:p></p>
        <p class="MsoNormal">Best,<o:p></o:p></p>
        <p class="MsoNormal"><o:p> </o:p></p>
        <p class="MsoNormal">Hien Vo<o:p></o:p></p>
        <p class="MsoNormal"><o:p> </o:p></p>
      </div>
      <br>
      <fieldset class="mimeAttachmentHeader"></fieldset>
      <pre class="moz-quote-pre" wrap="">_______________________________________________
Quantum ESPRESSO is supported by MaX (<a class="moz-txt-link-abbreviated" href="http://www.max-centre.eu/quantum-espresso">www.max-centre.eu/quantum-espresso</a>)
users mailing list <a class="moz-txt-link-abbreviated" href="mailto:users@lists.quantum-espresso.org">users@lists.quantum-espresso.org</a>
<a class="moz-txt-link-freetext" href="https://lists.quantum-espresso.org/mailman/listinfo/users">https://lists.quantum-espresso.org/mailman/listinfo/users</a></pre>
    </blockquote>
  </body>
</html>