<html>
<head>
<style>
.hmmessage P
{
margin:0px;
padding:0px
}
body.hmmessage
{
FONT-SIZE: 9pt;
FONT-FAMILY:Tahoma
}
</style>
</head>
<body class='hmmessage'>
<FONT color=#000000>Dear users:</FONT><BR>
<FONT color=#000000></FONT>&nbsp;<BR>
<FONT color=#000000>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;I have calculated the phonon dispersion of zinc-blende CrSb. Before calculating phonon dispersion, I relaxed the CrSb lattice and got the optimized lattice parameters,which agreed well with that in previous literature, and&nbsp;after ph.x calculation I found that the two lower frequences at X(-1.0 0.0 0.0) are imaginary,which indicates the lattice instability of two modes.The frequencies are shown as following</FONT><BR>
<FONT color=#000000></FONT>&nbsp;<BR>
<FONT color=#000000>&nbsp;&nbsp;  q =&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; -1.0000&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.0000&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.0000<BR>&nbsp;**************************************************************************<BR>&nbsp;&nbsp;&nbsp;&nbsp; omega( 1) =&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; -2.790113 [THz] =&nbsp;&nbsp;&nbsp;&nbsp; -93.068763 [cm-1]<BR>&nbsp;(&nbsp; 0.000000&nbsp;&nbsp; 0.000000&nbsp;&nbsp;&nbsp; -0.403189&nbsp;&nbsp; 0.000000&nbsp;&nbsp;&nbsp;&nbsp; 0.754599&nbsp;&nbsp; 0.000000&nbsp;&nbsp; )<BR>&nbsp;(&nbsp; 0.000000&nbsp;&nbsp; 0.000000&nbsp;&nbsp;&nbsp;&nbsp; 0.456614&nbsp;&nbsp; 0.000000&nbsp;&nbsp;&nbsp; -0.243973&nbsp;&nbsp; 0.000000&nbsp;&nbsp; )<BR>&nbsp;&nbsp;&nbsp;&nbsp; omega( 2) =&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; -2.790113 [THz] =&nbsp;&nbsp;&nbsp;&nbsp; -93.068763 [cm-1] <BR>&nbsp;(&nbsp; 0.000000&nbsp;&nbsp; 0.000000&nbsp;&nbsp;&nbsp; -0.754599&nbsp;&nbsp; 0.000000&nbsp;&nbsp;&nbsp; -0.403189&nbsp;&nbsp; 0.000000&nbsp;&nbsp; )<BR>&nbsp;(&nbsp; 0.000000&nbsp;&nbsp; 0.000000&nbsp;&nbsp;&nbsp; -0.243973&nbsp;&nbsp; 0.000000&nbsp;&nbsp;&nbsp; -0.456614&nbsp;&nbsp; 0.000000&nbsp;&nbsp; )<BR>&nbsp;&nbsp;&nbsp;&nbsp; omega( 3) =&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 3.418041 [THz] =&nbsp;&nbsp;&nbsp;&nbsp; 114.014317 [cm-1] <BR>&nbsp;(&nbsp; 0.000000&nbsp;&nbsp; 0.000000&nbsp;&nbsp;&nbsp;&nbsp; 0.000000&nbsp;&nbsp; 0.000000&nbsp;&nbsp;&nbsp;&nbsp; 0.000000&nbsp;&nbsp; 0.000000&nbsp;&nbsp; )<BR>&nbsp;( -1.000000&nbsp;&nbsp; 0.000000&nbsp;&nbsp;&nbsp;&nbsp; 0.000000&nbsp;&nbsp; 0.000000&nbsp;&nbsp;&nbsp;&nbsp; 0.000000&nbsp;&nbsp; 0.000000&nbsp;&nbsp; )<BR>&nbsp;&nbsp;&nbsp;&nbsp; omega( 4) =&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 3.580039 [THz] =&nbsp;&nbsp;&nbsp;&nbsp; 119.418057 [cm-1]<BR>&nbsp;(&nbsp; 1.000000&nbsp;&nbsp; 0.000000&nbsp;&nbsp;&nbsp;&nbsp; 0.000000&nbsp;&nbsp; 0.000000&nbsp;&nbsp;&nbsp;&nbsp; 0.000000&nbsp;&nbsp; 0.000000&nbsp;&nbsp; )<BR>&nbsp;(&nbsp; 0.000000&nbsp;&nbsp; 0.000000&nbsp;&nbsp;&nbsp;&nbsp; 0.000000&nbsp;&nbsp; 0.000000&nbsp;&nbsp;&nbsp;&nbsp; 0.000000&nbsp;&nbsp; 0.000000&nbsp;&nbsp; )<BR>&nbsp;&nbsp;&nbsp;&nbsp; omega( 5) =&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 5.338839 [THz] =&nbsp;&nbsp;&nbsp;&nbsp; 178.085673 [cm-1]<BR>&nbsp;(&nbsp; 0.000000&nbsp;&nbsp; 0.000000&nbsp;&nbsp;&nbsp; -0.817099&nbsp;&nbsp; 0.000000&nbsp;&nbsp;&nbsp;&nbsp; 0.000000&nbsp;&nbsp; 0.000000&nbsp;&nbsp; )<BR>&nbsp;(&nbsp; 0.000000&nbsp;&nbsp; 0.000000&nbsp;&nbsp;&nbsp;&nbsp; 0.000000&nbsp;&nbsp; 0.000000&nbsp;&nbsp;&nbsp;&nbsp; 0.576498&nbsp;&nbsp; 0.000000&nbsp;&nbsp; )<BR>&nbsp;&nbsp;&nbsp;&nbsp; omega( 6) =&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 5.338839 [THz] =&nbsp;&nbsp;&nbsp;&nbsp; 178.085673 [cm-1]<BR>&nbsp;(&nbsp; 0.000000&nbsp;&nbsp; 0.000000&nbsp;&nbsp;&nbsp;&nbsp; 0.000000&nbsp;&nbsp; 0.000000&nbsp;&nbsp;&nbsp; -0.817099&nbsp;&nbsp; 0.000000&nbsp;&nbsp; )<BR>&nbsp;(&nbsp; 0.000000&nbsp;&nbsp; 0.000000&nbsp;&nbsp;&nbsp;&nbsp; 0.576498&nbsp;&nbsp; 0.000000&nbsp;&nbsp;&nbsp;&nbsp; 0.000000&nbsp;&nbsp; 0.000000&nbsp;&nbsp; )<BR>&nbsp;**************************************************************************</FONT><BR>
<FONT color=#000000></FONT>&nbsp;<BR>
<FONT color=#000000>Usuallly, <FONT color=#ff0000>an instability implys that the atoms&nbsp;are not located &nbsp;in local minimum energy postion</FONT>,thus next I displaced atoms in the cell accroding to the eigen-displacements shown above,and got a cubic cell containg 8 atoms. </FONT><FONT color=#ff0000>But I found the energy rised with the increase of atom displacements, which was totally contrary to what I expected. In fact, only the decrease of energy with atom displacments will be consistent the results from phonon self-consitsten calculation.</FONT><BR>
<FONT color=#000000></FONT>&nbsp;<BR>
<FONT color=#000000>The electronic structure calculation input file and phonon input file are list as : </FONT><BR>
<FONT color=#000000></FONT>&nbsp;<BR>
&nbsp;<BR>
<FONT color=#000000></FONT>&nbsp;<BR>
<FONT color=#000000></FONT>&nbsp;<BR>
<FONT color=#000000>###&nbsp;&nbsp;&nbsp; electronic structure calculation </FONT><BR>
<FONT color=#000000>####</FONT><BR>
<FONT color=#000000>&nbsp;&amp;CONTROL<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; title = CrSb ,<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; calculation = 'scf' ,<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; restart_mode = 'from_scratch' ,<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; prefix='CrSb',<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; pseudo_dir = '/home/phymwm/PWSCF/espresso-3.2.3/pseudo' ,<BR>&nbsp;/<BR>&nbsp;&amp;SYSTEM<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ibrav = 2,<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; celldm(1) = 11.55,<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; nat = 2,<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ntyp = 2,<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ecutwfc=60,<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ecutrho=600,<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; occupations = 'smearing' ,<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; degauss = 0.025 ,<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; smearing = 'methfessel-paxton' ,<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; nspin = 2 ,<BR>&nbsp;&nbsp; starting_magnetization(1) = 3.0,<BR>&nbsp;&nbsp; starting_magnetization(2) = 1.0,<BR>&nbsp;/<BR>&nbsp;&amp;ELECTRONS<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; electron_maxstep = 200,<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; conv_thr = 1.0D-8 ,<BR>&nbsp;/<BR>ATOMIC_SPECIES<BR>&nbsp;&nbsp; Cr&nbsp;&nbsp; 52.00000&nbsp; Cr.pw91-sp-van.UPF<BR>&nbsp;&nbsp; Sb&nbsp; 121.80000&nbsp; Sb.pw91-n-van.UPF<BR>ATOMIC_POSITIONS alat<BR>&nbsp;&nbsp; Cr&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.000000000&nbsp;&nbsp;&nbsp; 0.000000000&nbsp;&nbsp;&nbsp; 0.000000000<BR>&nbsp;&nbsp; Sb&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.250000000&nbsp;&nbsp;&nbsp; 0.250000000&nbsp;&nbsp;&nbsp; 0.250000000<BR>K_POINTS automatic<BR>&nbsp; 8 8 8&nbsp;&nbsp; 0 0 0<BR>&nbsp;</FONT><BR>
<FONT color=#000000></FONT>&nbsp;<BR>
<FONT color=#000000>#####</FONT><BR>
<FONT color=#000000>###</FONT><BR>
phonon dispersion for CrSb<BR>&nbsp;&amp;inputph<BR>&nbsp; tr2_ph=1.0d-16,<BR>&nbsp; prefix='CrSb',<BR>&nbsp; fildvscf='dv',<BR>&nbsp; amass(1)=52.00,<BR>&nbsp; amass(2)=121.8,<BR>&nbsp; outdir='./',<BR>&nbsp; fildyn='CrSb.dyn',<BR>&nbsp; trans=.true.,<BR>&nbsp; ldisp=.true.<BR>&nbsp; nq1=4, nq2=4, nq3=4<BR>&nbsp;/<BR>
&nbsp;<BR>
&nbsp;<BR>
&nbsp;<BR>
Thanks in advance !<BR>
&nbsp;<BR>
Wenmei Ming <BR>
&nbsp;<BR><br /><hr />轻松把Hotmail下载到本地,试试 Windows Live Mail。 <a href='http://get.live.com/wl/all' target='_new'>立即尝试!</a></body>
</html>