<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-15">
</head>
<body text="#000000" bgcolor="#FFFFFF">
<p>Hello, i'm Francesco Moriggi, a master student at the Milan University, i have a problem with ph.x. I am calculating NO frequency and i'm not reaching a convergence in values compared to cutoff energy, using norm conserving and ultrasoft pseudopotentials.
Values i get are (i've put an image too in the email): </p>
<div class="moz-forward-container">
<p><b>Cutoff freq Norm freq US</b></p>
<b></b>
<p><b>75 1879.0 1844.4</b></p>
<b></b>
<p><b>80 1875.5 1858.2</b></p>
<b></b>
<p><b>85 1890.6 1907.1</b></p>
<b></b>
<p><b>90 1923.5 1895.2</b></p>
<b></b>
<p><b>100 1889.6 1909.3</b></p>
<b></b>
<p><b>110 1914.2 1893.1</b></p>
<b></b>
<p><b>120 1890.5 1900.4</b></p>
<b></b>
<p><b>130 1890.7 1901.0</b></p>
<b></b>
<p><b>140 1914.7 1890.9</b></p>
<b></b>
<p><b>150 1886.9 ////</b></p>
<p>With Gaussian (aug-cc-pv5Z) i've got <b>1889.21</b> cm-1.<br>
</p>
<p>An example of the input of geometry optimization is:</p>
<p><b>&CONTROL</b><b><br>
</b><b> calculation = "relax",</b><b><br>
</b><b> prefix = "NO",</b><b><br>
</b><b> pseudo_dir = "/galileo/home/userexternal/fmoriggi/QE/NO/pseudoNO/ultrasoft",</b><b><br>
</b><b> outdir = "./out",</b><b><br>
</b><b>/</b><b><br>
</b><b>&SYSTEM</b><b><br>
</b><b> ibrav = 1,</b><b><br>
</b><b> a = 20,</b><b><br>
</b><b> nat = 2,</b><b><br>
</b><b> ntyp = 2,</b><b><br>
</b><b> nspin = 2,</b><b><br>
</b><b> ecutwfc = 150.D0,</b><b><br>
</b><b> ecutrho = 1200.D0,</b><b><br>
</b><b> occupations = 'smearing', smearing = 'gauss', degauss=0.02,</b><b><br>
</b><b> starting_magnetization(1)=-0.9,</b><b><br>
</b><b> starting_magnetization(2)=1.5,</b><b><br>
</b><b>/</b><b><br>
</b><b>&ELECTRONS</b><b><br>
</b><b> electron_maxstep = 500,</b><b><br>
</b><b> conv_thr = 1.D-7,</b><b><br>
</b><b> mixing_beta = 0.7D0,</b><b><br>
</b><b>/</b><b><br>
</b><b>&IONS</b><b><br>
</b><b>/</b><b><br>
</b><b>ATOMIC_SPECIES</b><b><br>
</b><b>O 16.00 O.pbe-van_ak.UPF</b><b><br>
</b><b>N 14.00 N.pbe-van_ak.UPF</b><b><br>
</b><b>ATOMIC_POSITIONS {angstrom}</b><b><br>
</b><b>N 1.16 0.0 0.0</b><b><br>
</b><b>O 0.000 0.0 0.0 0 0 0</b><b><br>
</b><b>K_POINTS</b><b><br>
</b><b>1</b><b><br>
</b><b>0.0 0.0 0.0 1.0</b></p>
<p>Later i've done an scf with a better conv_thr. This is the input.<br>
</p>
<p><b>&CONTROL</b><b><br>
</b><b> calculation = "scf",</b><b><br>
</b><b> prefix = "150",</b><b><br>
</b><b> pseudo_dir = "/galileo/home/userexternal/fmoriggi/QE/NO/pseudoNO/ultrasoft",</b><b><br>
</b><b> outdir = "./out150",</b><b><br>
</b><b>/</b><b><br>
</b><b>&SYSTEM</b><b><br>
</b><b> ibrav = 1,</b><b><br>
</b><b> a = 20,</b><b><br>
</b><b> nat = 2,</b><b><br>
</b><b> ntyp = 2,</b><b><br>
</b><b> nspin = 2,</b><b><br>
</b><b> ecutwfc = 150.D0,</b><b><br>
</b><b> ecutrho = 1200.D0,</b><b><br>
</b><b> occupations = 'smearing', smearing = 'gauss', degauss=0.02,</b><b><br>
</b><b> starting_magnetization(1)=-0.9,</b><b><br>
</b><b> starting_magnetization(2)=1.5,</b><b><br>
</b><b>/</b><b><br>
</b><b>&ELECTRONS</b><b><br>
</b><b> electron_maxstep = 500,</b><b><br>
</b><b> conv_thr = 1.D-9,</b><b><br>
</b><b> mixing_beta = 0.7D0,</b><b><br>
</b><b>/</b><b><br>
</b><b>&IONS</b><b><br>
</b><b>/</b><b><br>
</b><b>ATOMIC_SPECIES</b><b><br>
</b><b>O 16.00 O.pbe-van_ak.UPF</b><b><br>
</b><b>N 14.00 N.pbe-van_ak.UPF</b><b><br>
</b><b>ATOMIC_POSITIONS {angstrom}</b><b><br>
</b><b>N 1.166278 0.0 0.0</b><b><br>
</b><b>O 0.000 0.0 0.0 0 0 0</b><b><br>
</b><b>K_POINTS</b><b><br>
</b><b>1</b><b><br>
</b><b>0.0 0.0 0.0 1.0</b></p>
<p>And then i run the ph.x with this input:</p>
<p><b>phonon calculation - NO</b><b><br>
</b><b> &inputph</b><b><br>
</b><b> tr2_ph=4.0d-12,</b><b><br>
</b><b> prefix='150'</b><b><br>
</b><b> amass(1)=14.007,</b><b><br>
</b><b> amass(2)=15.9994,</b><b><br>
</b><b> outdir='/galileo/home/userexternal/fmoriggi/QE/NO/relax/ultra/out150'</b><b><br>
</b><b> /</b><b><br>
</b><b> 0.0 0.0 0.0</b></p>
<p>If you can help me it would be very useful to my future TiO2 based calculations.<br>
<br>
<br>
Best regards,</p>
<p>Francesco Moriggi<br>
</p>
<p>Webpage: <a class="moz-txt-link-freetext" href="http://users.unimi.it/ceotto/about_FrancescoMoriggi.html">
http://users.unimi.it/ceotto/about_FrancescoMoriggi.html</a><br>
</p>
</div>
</body>
</html>