<div>Dear users:</div> <div> i want to calculte the phonons of SrTiO3 which Sr atom is fixed. i do it by using "nat_todo".</div> <div>there is my inputfile of ph.x:<BR> &INPUTPH<BR> amass(1) = 87.62000,<BR> amass(2) = 47.88000,<BR> amass(3) = 15.99960,<BR> outdir = 'tmp/' ,<BR> prefix = 'st'
,<BR> tr2_ph = 1.0d-14 ,<BR> fildyn = 'st.dyn5' ,<BR> nat_todo= 4<BR> epsil = .true.,<BR> trans = .true.,<BR> elph =
.false.,<BR> zue = .true.,<BR>/<BR>0 0 0<BR>2 3 4 5</div> <div>but in dyn file there is not phonons. by dynamat.x i deal with the dyn file.the results is as follows.</div> <div> omega( 1) = -2.746240 [THz] = -91.605330 [cm-1]<BR> ( -0.229443 0.000000 -0.099229 0.000000 0.000000 0.000000 )<BR> ( -0.089673 0.000000 -0.038782 0.000000 0.000000 0.000000 )<BR> ( 0.541628 0.000000 0.234241 0.000000 0.000000 0.000000 )<BR> ( 0.541628
0.000000 0.190990 0.000000 0.000000 0.000000 )<BR> ( 0.441620 0.000000 0.234241 0.000000 0.000000 0.000000 )<BR> omega( 2) = -2.746240 [THz] = -91.605330 [cm-1]<BR> ( 0.099229 0.000000 -0.229443 0.000000 0.000000 0.000000 )<BR> ( 0.038782 0.000000 -0.089673 0.000000 0.000000 0.000000 )<BR> ( -0.234241 0.000000 0.541628 0.000000 0.000000 0.000000 )<BR> ( -0.234241 0.000000
0.441620 0.000000 0.000000 0.000000 )<BR> ( -0.190990 0.000000 0.541628 0.000000 0.000000 0.000000 )<BR> omega( 3) = 0.000000 [THz] = -0.000009 [cm-1]<BR> ( -0.037890 0.000000 -0.025352 0.000000 0.444884 0.000000 )<BR> ( -0.037890 0.000000 -0.025352 0.000000 0.444884 0.000000 )<BR> ( -0.037890 0.000000 -0.025352 0.000000 0.444884 0.000000 )<BR> ( -0.037890 0.000000 -0.025352
0.000000 0.444884 0.000000 )<BR> ( -0.037890 0.000000 -0.025352 0.000000 0.444884 0.000000 )<BR> omega( 4) = 0.000000 [THz] = 0.000005 [cm-1]<BR> ( -0.180783 0.000000 0.408968 0.000000 0.007909 0.000000 )<BR> ( -0.180783 0.000000 0.408968 0.000000 0.007909 0.000000 )<BR> ( -0.180783 0.000000 0.408968 0.000000 0.007909 0.000000 )<BR> ( -0.180783 0.000000 0.408968
0.000000 0.007909 0.000000 )<BR> ( -0.180783 0.000000 0.408968 0.000000 0.007909 0.000000 )<BR> omega( 5) = 0.000000 [THz] = 0.000005 [cm-1]<BR> ( 0.407286 0.000000 0.179171 0.000000 0.044898 0.000000 )<BR> ( 0.407286 0.000000 0.179171 0.000000 0.044898 0.000000 )<BR> ( 0.407286 0.000000 0.179171 0.000000 0.044898 0.000000 )<BR> ( 0.407286 0.000000 0.179171
0.000000 0.044898 0.000000 )<BR> ( 0.407286 0.000000 0.179171 0.000000 0.044898 0.000000 )<BR> omega( 6) = 3.509935 [THz] = 117.079598 [cm-1]<BR> ( 0.000000 0.000000 0.000000 0.000000 -0.513490 0.000000 )<BR> ( 0.000000 0.000000 0.000000 0.000000 0.574186 0.000000 )<BR> ( 0.000000 0.000000 0.000000 0.000000 0.292212 0.000000 )<BR> ( 0.000000 0.000000 0.000000
0.000000 0.400781 0.000000 )<BR> ( 0.000000 0.000000 0.000000 0.000000 0.400781 0.000000 )<BR> omega( 7) = 4.714264 [THz] = 157.251966 [cm-1]<BR> ( 0.384775 0.000000 0.055007 0.000000 0.000000 0.000000 )<BR> ( -0.866552 0.000000 -0.123881 0.000000 0.000000 0.000000 )<BR> ( 0.181436 0.000000 0.025938 0.000000 0.000000 0.000000 )<BR> ( 0.181436 0.000000 0.017608
0.000000 0.000000 0.000000 )<BR> ( 0.123170 0.000000 0.025938 0.000000 0.000000 0.000000 )<BR> omega( 8) = 4.714264 [THz] = 157.251966 [cm-1]<BR> ( 0.055007 0.000000 -0.384775 0.000000 0.000000 0.000000 )<BR> ( -0.123881 0.000000 0.866552 0.000000 0.000000 0.000000 )<BR> ( 0.025938 0.000000 -0.181436 0.000000 0.000000 0.000000 )<BR> ( 0.025938 0.000000 -0.123170
0.000000 0.000000 0.000000 )<BR> ( 0.017608 0.000000 -0.181436 0.000000 0.000000 0.000000 )<BR> omega( 9) = 6.891958 [THz] = 229.892506 [cm-1]<BR> ( 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 )<BR> ( 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 )<BR> ( 0.000905 0.000000 0.005360 0.000000 0.000000 0.000000 )<BR> ( -0.000905 0.000000 0.000000
0.000000 0.707086 0.000000 )<BR> ( 0.000000 0.000000 -0.005360 0.000000 -0.707086 0.000000 )<BR> omega(10) = 6.891958 [THz] = 229.892506 [cm-1]<BR> ( 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 )<BR> ( 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 )<BR> ( -0.703970 0.000000 0.066527 0.000000 0.000000 0.000000 )<BR> ( 0.703970 0.000000 0.000000
0.000000 0.000397 0.000000 )<BR> ( 0.000000 0.000000 -0.066527 0.000000 -0.000397 0.000000 )<BR> omega(11) = 6.891958 [THz] = 229.892506 [cm-1]<BR> ( 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 )<BR> ( 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 )<BR> ( 0.066522 0.000000 0.703950 0.000000 0.000000 0.000000 )<BR> ( -0.066522 0.000000 0.000000
0.000000 -0.005421 0.000000 )<BR> ( 0.000000 0.000000 -0.703950 0.000000 0.005421 0.000000 )<BR> omega(12) = 13.570374 [THz] = 452.661953 [cm-1]<BR> ( 0.000000 0.000000 0.000000 0.000000 -0.020309 0.000000 )<BR> ( 0.000000 0.000000 0.000000 0.000000 -0.269344 0.000000 )<BR> ( 0.000000 0.000000 0.000000 0.000000 -0.350802 0.000000 )<BR> ( 0.000000 0.000000 0.000000 0.000000
0.634026 0.000000 )<BR> ( 0.000000 0.000000 0.000000 0.000000 0.634026 0.000000 )<BR> omega(13) = 16.727499 [THz] = 557.973034 [cm-1]<BR> ( -0.011300 0.000000 0.000932 0.000000 0.000000 0.000000 )<BR> ( -0.019840 0.000000 0.001636 0.000000 0.000000 0.000000 )<BR> ( -0.365336 0.000000 0.030127 0.000000 0.000000 0.000000 )<BR> ( -0.365336 0.000000 -0.070254 0.000000 0.000000
0.000000 )<BR> ( 0.851929 0.000000 0.030127 0.000000 0.000000 0.000000 )<BR> omega(14) = 16.727499 [THz] = 557.973034 [cm-1]<BR> ( -0.000932 0.000000 -0.011300 0.000000 0.000000 0.000000 )<BR> ( -0.001636 0.000000 -0.019840 0.000000 0.000000 0.000000 )<BR> ( -0.030127 0.000000 -0.365336 0.000000 0.000000 0.000000 )<BR> ( -0.030127 0.000000 0.851929 0.000000 0.000000 0.000000 )<BR> (
0.070254 0.000000 -0.365336 0.000000 0.000000 0.000000 )<BR> omega(15) = 23.660653 [THz] = 789.239679 [cm-1]<BR> ( 0.000000 0.000000 0.000000 0.000000 -0.049737 0.000000 )<BR> ( 0.000000 0.000000 0.000000 0.000000 -0.280773 0.000000 )<BR> ( 0.000000 0.000000 0.000000 0.000000 0.951708 0.000000 )<BR> ( 0.000000 0.000000 0.000000 0.000000 0.080452 0.000000 )<BR> ( 0.000000
0.000000 0.000000 0.000000 0.080452 0.000000 )<BR>in the results, i found eignvector of Sr atom is not zero.but i found the eignvector of Sr atom in SrTiO3 is zero when Sr atom is fixed.please give me some suggestion how to modify the inputfile of ph.x. </div> <div> </div> <div> </div> <div>Regards</div> <div><BR> </div> <div> zhuzhenye</div><p>
<hr size=1><a href="http://cn.mail.yahoo.com/" target=blank>
ÑÅ»¢Ãâ·ÑÓÊÏä-3.5GÈÝÁ¿£¬20M¸½¼þ</a>