<html dir="ltr">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<style id="owaParaStyle">P {
        MARGIN-TOP: 0px; MARGIN-BOTTOM: 0px
}
</style>
</head>
<body fPStyle="1" ocsi="0">
<div style="FONT-SIZE: 13px; COLOR: #000000; DIRECTION: ltr; FONT-FAMILY: Tahoma">
<div>Dear pwscf users</div>
<div style="FONT-SIZE: 16px; COLOR: #000000; FONT-FAMILY: Times New Roman">
<div>
<div style="FONT-SIZE: 13px; COLOR: #000000; DIRECTION: ltr; FONT-FAMILY: Tahoma">
<div>&nbsp;</div>
<div>I've been trying to do phonon calculations, for elemental and compound semiconductors(Si, Ge, InP and SiC). Si and Ge both have exactly the same structure and space group, so symmetry operations are also same. Below is&nbsp;the self consistent input file, i'm
 using.&nbsp;In case of Si and&nbsp;Ge, the scf.out file works with all the 48 symmetry operations but in case of Ge it rejects all the inversion symmetry opertors. Also, for further calculations i'm using exactly the procedure given in example07 directory.
</div>
<div>&nbsp;</div>
<div>Also, for InP, I used the same procedure, but i'm missing the splitting of phonon spectra at higher frequency, as reported in publications. After going through the forums, I got to know that I should use zasr=.true. and epsil=.true. for calculations in
 InP, since it is a polar semiconductor. I tried this and the phonon calculation job keeps on crashing. I don't understand what should I do to solve this problem?</div>
<div>&nbsp;</div>
<div><br>
&nbsp;<em>&amp;control<br>
&nbsp;&nbsp;&nbsp; calculation='scf'<br>
&nbsp;&nbsp;&nbsp; restart_mode='from_scratch',<br>
&nbsp;&nbsp;&nbsp; prefix='ge',<br>
&nbsp;&nbsp;&nbsp; pseudo_dir = '/nfs/18/osu5465/espresso/espresso-3.2/pseudo/',<br>
&nbsp;&nbsp;&nbsp; outdir='/nfs/18/osu5465/espresso/tmp/'<br>
&nbsp;/<br>
&nbsp;&amp;system<br>
&nbsp;&nbsp;&nbsp; ibrav= 2, celldm(1) =10.691, nat=2 , ntyp= 1,<br>
&nbsp;&nbsp;&nbsp; ecutwfc =25.0,&nbsp;&nbsp; <br>
&nbsp;&nbsp;&nbsp; occupations='smearing', smearing='methfessel-paxton', degauss=0.05,<br>
&nbsp;/<br>
&nbsp;&amp;electrons<br>
&nbsp;&nbsp;&nbsp; conv_thr =&nbsp; 1.0d-8<br>
&nbsp;&nbsp;&nbsp; mixing_beta = 0.7 <br>
&nbsp;/<br>
ATOMIC_SPECIES<br>
&nbsp;Ge&nbsp; 72.61&nbsp;&nbsp; Ge.pz-bhs.UPF<br>
ATOMIC_POSITIONS {crystal}<br>
Ge&nbsp;&nbsp;&nbsp;&nbsp; 0.250000&nbsp; 0.250000&nbsp;&nbsp;&nbsp; 0.25000000<br>
Ge&nbsp;&nbsp;&nbsp;&nbsp; 0.000000&nbsp; 0.000000&nbsp;&nbsp;&nbsp; 0.00000000<br>
K_POINTS {automatic}<br>
&nbsp;8 8 8 0 0 0<br>
</em></div>
<div><em>&nbsp;</div>
</em>
<div>--</div>
<div class="BodyFragment"><font size="2"><span style="FONT-SIZE: 10pt">
<div class="PlainText">Ashutosh Kumar<br>
Grad Student<br>
The Ohio State University<br>
</div>
</span></font></div>
</div>
</div>
</div>
</div>
</body>
</html>