<div>Dear all</div>  <div> </div>  <div>I met a problem when I used ph.x to calculate Dielectric constant of BaTiO3. The Dielectric constants are abnormally large, and so do the effective charges. The output and input files are attached below. I have searched the forum and find someone fixed this problem by using denser k-grid. However, in my case this does not work. Even I use a 16x16x16 MP k-point, the problem also exsits. <BR>I also calculate effective charges of BaTiO3 using Berry Phase method, and the result is reasonable. It is not likely pseudopotentials are the problem, because Linear Response calculations of PbTiO3, BaZrO3, PbZrO3 all give reasonable results, and agree well with Berry Phase results.</div>  <div> </div>  <div>Really appreciate your help, thanks!<BR>Ruizhi Zhang</div>  <div><BR>Here are the input and output files</div>  <div>bto.ph.out<BR>.......<BR>.......<BR> Dielectric constant in cartesian axis</div>  <div> ( 
 -10960.282004778       0.000000000       0.000000000 )<BR> (       0.000000000  -10960.282004778       0.000000000 )<BR> (       0.000000000       0.000000000  -10960.282004778 )</div>  <div> Effective charges E-U in cartesian axis</div>  <div>  atom      1<BR> (      -30.52734        0.00000        0.00000 )<BR> (        0.00000      -30.52734        0.00000 )<BR> (        0.00000        0.00000      -30.52734 )<BR> 
 atom      2<BR> (      -48.86837        0.00000        0.00000 )<BR> (        0.00000      -48.86837        0.00000 )<BR> (        0.00000        0.00000      -48.86837 )<BR>  atom      3<BR> (       25.10465        0.00000        0.00000 )<BR> (        0.00000       27.75341        0.00000 )<BR> (        0.00000       
 0.00000       27.75341 )<BR>  atom      4<BR> (       27.75341        0.00000        0.00000 )<BR> (        0.00000       27.75341        0.00000 )<BR> (        0.00000        0.00000       25.10465 )<BR>  atom      5<BR> (       27.75341        0.00000        0.00000 )<BR> (        0.00000       25.10465        0.00000
 )<BR> (        0.00000        0.00000       27.75341 )<BR>.......<BR>.......</div>  <div>bto.scf.in<BR>&control<BR>    calculation  = 'scf'<BR>    restart_mode = 'from_scratch'<BR>    prefix = 'bto'<BR>    pseudo_dir   = '/home/zrz/espresso-3.1.1/pseudo/'<BR>    outdir       = '/home/zrz/phonon/bto.cub/tmp/'<BR> /<BR> &system<BR>    ibrav=1<BR>    celldm(1)=7.56<BR>    nat=5<BR>    ntyp=3<BR>    ecutwfc=30<BR>  /<BR> &electrons<BR>    conv_thr = 1e-7<BR>    mixing_beta=0.7<BR> /<BR>&IONS<BR> /<BR>ATOMIC_SPECIES<BR>  Ba    137.33     Ba.vdb.UPF<BR> 
 Ti    47.90      Ti.vdb.UPF<BR>  O     15.9994    O.vdb.UPF<BR>ATOMIC_POSITIONS<BR>  Ba    0.000    0.000    0.000<BR>  Ti    0.500    0.500    0.500<BR>  O     0.000    0.500    0.500<BR>  O     0.500    0.500    0.000<BR>  O     0.500    0.000    0.500<BR>K_POINTS {automatic}<BR>4 4 4 0 0 0</div>  <div><BR>bto.ph.in<BR>phonons of bto<BR>&inputph<BR>  tr2_ph=1.0d-12,<BR>  prefix='bto',<BR>  epsil = .true.<BR>  trans = .false.<BR>  amass(1)=137.33,<BR>  amass(2)=47.90,<BR>  amass(3)=15.9996<BR>  outdir='/home/zrz/phonon/bto.cub/tmp/',<BR>  fildyn='bto.dyn',<BR> /<BR>0.0 0.0
 0.0</div>  <div> </div><p>

<hr size=1>Finding fabulous fares is fun.<br><a href="http://farechase.yahoo.com/promo-generic-14795097;_ylc=X3oDMTFtNW45amVpBF9TAzk3NDA3NTg5BF9zAzI3MTk0ODEEcG9zAzEEc2VjA21haWx0YWdsaW5lBHNsawNxMS0wNw--
">Let Yahoo! FareChase search your favorite travel sites</a> to find flight and hotel bargains.