<br><br><div><span class="gmail_quote">   Dear pwscf users,<br>
<br>
</span>   My  calculated gamma phonon frequency for Si (spin-less)  is<br>
    = 510 cm-1    When  I switch on the spin  tag to specify the spin polarized<br>
   calculation for   the same sample, then I got   corresponding<br>
   gamma phonon frequency = 214 cm-1  <br>
   while all the other input parameters are exactly same.<div><br>
   <br>
   Paolo wrote:  </div>   please provide a job demostrating such behaviour<br>
<br>
   Dear Paolo, <br>
   I'm also sending the following input file what is demostrating such behaviour<br>
   by attachment.<br>
   <br>
 **************************************************************************************<br>
   #!/bin/sh<br>
 executable_dir=./<br>
 tmp_dir=./tmp<br>
 pseudo_dir=./pseudo<br>
cat > <a href="http://si.scf.in">si.scf.in</a> << EOF<br>
 &control<br>
    calculation='scf',<br>
    restart_mode='from_scratch',<br>
    prefix='si'<br>
    pseudo_dir = '$pseudo_dir/',<br>
    outdir='$tmp_dir/'<br>
 /<br>
 &system<br>
    ibrav = 2, celldm(1) =10.20, nat=  2, ntyp= 1,<br>
    ecutwfc = 18.0, nelec  =   8.00, nbnd =10<br>
    nspin=2, tot_magnetization=0,<br>
    occupations='smearing', smearing='fd', degauss=0.005,<br>
/<br>
 &electrons<br>
    mixing_beta = 0.7<br>
    conv_thr =  1.0d-8<br>
 /<br>
ATOMIC_SPECIES<br>
 Si  28.086  Si.vbc.UPF<br>
ATOMIC_POSITIONS<br>
 Si 0.00 0.00 0.00<br>
 Si 0.25 0.25 0.25<br>
K_POINTS<br>
  10<br>
   0.1250000  0.1250000  0.1250000   1.00<br>
   0.1250000  0.1250000  0.3750000   3.00<br>
   0.1250000  0.1250000  0.6250000   3.00<br>
   0.1250000  0.1250000  0.8750000   3.00<br>
   0.1250000  0.3750000  0.3750000   3.00<br>
   0.1250000  0.3750000  0.6250000   6.00<br>
   0.1250000  0.3750000  0.8750000   6.00<br>
   0.1250000  0.6250000  0.6250000   3.00<br>
   0.3750000  0.3750000  0.3750000   1.00<br>
   0.3750000  0.3750000  0.6250000   3.00<br>
EOF<br>
$ECHO "  running the scf calculation...\c"<br>
$executable_dir/pw.x  <  <a href="http://si.scf.in">si.scf.in</a> > si.scf.out<br>
$ECHO " done"<br>
<br>
# phonon calculation at Gamma<br>
cat > <a href="http://si.phG.in">si.phG.in</a> << EOF<br>
phonons of Si at Gamma<br>
 &inputph<br>
  tr2_ph=1.0d-14,<br>
  prefix='si',<br>
  amass(1)=28.08,<br>
  outdir='$tmp_dir/',<br>
  fildyn='si.dynG',<br>
 /<br>
0.0 0.0 0.0<br>
EOF<br>
$ECHO "  running the phonon calculation at Gamma...\c"<br>
$executable_dir/ph.x < <a href="http://si.phG.in">si.phG.in</a> > si.phG.out<br>
$ECHO " done"<br>
</div><br>
********************************************************************************<br>
  To run the spin less calculation I just remove tot_magnetization tag <br>
  and reset nspin=1 and I get correct Gamma point<br>
  phonon frequency (510 cm-1) whereas above input file for spin calculation<br>
  will give us a completely wrong result (214 cm-1).<br>
 <br>
  Looking forward to your valuable suggestions.<br>
 <br>
  With best regards,<br>
  Srijan<br>