<html>
<head>
<meta http-equiv="Content-Type" content="text/html;
charset=windows-1252">
</head>
<body text="#000000" bgcolor="#FFFFFF">
<div class="moz-cite-prefix">Hi Casey <br>
using your inputs cut and paste I get the right values. <br>
<br>
# mode [cm-1] [THz] IR<br>
1 -0.00 -0.0000 0.0000<br>
2 -0.00 -0.0000 0.0000<br>
3 0.00 0.0000 0.0000<br>
4 291.10 8.7268 11.3884<br>
5 291.10 8.7268 11.3884<br>
6 569.57 17.0753 11.3884<br>
<br>
What version are you using ? <br>
<br>
<br>
On 03/27/2019 03:19 PM, Casey Brock wrote:<br>
</div>
<blockquote type="cite"
cite="mid:DM6PR18MB2713B0BD9E2C78437436EC8AFB580@DM6PR18MB2713.namprd18.prod.outlook.com">
<meta http-equiv="Content-Type" content="text/html;
charset=windows-1252">
<style type="text/css" style="display:none;"> P {margin-top:0;margin-bottom:0;} </style>
<div style="font-family: Calibri, Helvetica, sans-serif;
font-size: 12pt; color: rgb(0, 0, 0);">
<span style="caret-color: rgb(0, 0, 0); font-family: Helvetica;
font-size: 12px; font-weight: normal; orphans: auto; widows:
auto; display: inline !important">I am using pw.x -> ph.x
-> dynmat.x to get the LO and TO phonon frequencies at
gamma for rock-salt CaO. The value for the LO phonon I am
getting differs from experiment and other DFPT calculations in
the literature by 5 THz (167 cm^-1). Perhaps this is a problem
with Quantum Espresso, but it’s very possible I’m doing
something wrong. I’ve included results and input files below,
and any advice would be greatly appreciated.</span>
<div style="caret-color: rgb(0, 0, 0); font-family: Helvetica;
font-size: 12px; font-weight: normal; orphans: auto; widows:
auto">
<br>
</div>
<div style="caret-color: rgb(0, 0, 0); font-family: Helvetica;
font-size: 12px; font-weight: normal; orphans: auto; widows:
auto">
<div>The TO and LO frequencies in THz are given below. Notice
that Quantum Espresso using DFPT gives a value that is too
low.</div>
<div>experiment: 8.9, 17.1</div>
<div>frozen phonon (Phonopy + Quantum Espresso): 8.71, 17.1</div>
<div>DFPT (Quantum Espresso, pw.x->ph.x->dynmat.x):
8.53, <b>12.2</b></div>
</div>
<div style="caret-color: rgb(0, 0, 0); font-family: Helvetica;
font-size: 12px; font-weight: normal; orphans: auto; widows:
auto">
<b><br>
</b></div>
<div style="caret-color: rgb(0, 0, 0); font-family: Helvetica;
font-size: 12px; font-weight: normal; orphans: auto; widows:
auto">
If I use Phonopy (frozen phonon technique) combined with
Quantum Espresso (pw.x and ph.x) I get the correct splitting
at gamma. I have tried the same calculation with MgO and I see
the same problem, where the frozen phonon technique works
while DFPT does not. Here are some things I have tried:
different PAW pseudopotentials, increasing ecutwfc and
ecutrho, increasing k-points, reducing conv_thr, and reducing
tr2_ph.</div>
<div style="caret-color: rgb(0, 0, 0); font-family: Helvetica;
font-size: 12px; font-weight: normal; orphans: auto; widows:
auto">
<br>
</div>
<div style="caret-color: rgb(0, 0, 0); font-family: Helvetica;
font-size: 12px; font-weight: normal; orphans: auto; widows:
auto">
Thanks in advance,</div>
<div style="caret-color: rgb(0, 0, 0); font-family: Helvetica;
font-size: 12px; font-weight: normal; orphans: auto; widows:
auto">
Casey Brock</div>
<div style="caret-color: rgb(0, 0, 0); font-family: Helvetica;
font-size: 12px; font-weight: normal; orphans: auto; widows:
auto">
<br>
</div>
<div style="caret-color: rgb(0, 0, 0); font-family: Helvetica;
font-size: 12px; font-weight: normal; orphans: auto; widows:
auto">
input to pw.x</div>
<div style="caret-color: rgb(0, 0, 0); font-family: Helvetica;
font-size: 12px; font-weight: normal; orphans: auto; widows:
auto">
-------------------------------------------</div>
<div style="caret-color: rgb(0, 0, 0); font-family: Helvetica;
font-size: 12px; font-weight: normal; orphans: auto; widows:
auto">
&control<br>
calculation='scf'<br>
restart_mode='from_scratch'<br>
prefix='cao'<br>
pseudo_dir = './'<br>
outdir='./output/'<br>
tprnfor = .true.<br>
tstress = .true.<br>
/<br>
&system<br>
ibrav = 0<br>
nat = 2<br>
ntyp = 2<br>
ecutwfc = 70.0<br>
ecutrho = 280.0<br>
/<br>
&electrons<br>
conv_thr = 1.0d-10<br>
mixing_beta = 0.7<br>
/<br>
K_POINTS {automatic}<br>
12 12 12 0 0 0<br>
<br>
CELL_PARAMETERS bohr<br>
4.45 4.45 0<br>
0 4.45 4.45<br>
4.45 0 4.45<br>
ATOMIC_SPECIES<br>
O 16.00 O.pz-n-kjpaw_psl.0.1.UPF<br>
Ca 40.08 Ca.pz-spn-kjpaw_psl.1.0.0.UPF<br>
ATOMIC_POSITIONS crystal<br>
Ca 0.0 0.0 0.0<br>
O 0.5 0.5 0.5</div>
<div style="caret-color: rgb(0, 0, 0); font-family: Helvetica;
font-size: 12px; font-weight: normal; orphans: auto; widows:
auto">
<br>
</div>
<div style="caret-color: rgb(0, 0, 0); font-family: Helvetica;
font-size: 12px; font-weight: normal; orphans: auto; widows:
auto">
input to ph.x</div>
<div style="caret-color: rgb(0, 0, 0); font-family: Helvetica;
font-size: 12px; font-weight: normal; orphans: auto; widows:
auto">
---------------------------------------</div>
<div style="caret-color: rgb(0, 0, 0); font-family: Helvetica;
font-size: 12px; font-weight: normal; orphans: auto; widows:
auto">
Phonon for CaO<br>
&inputph<br>
prefix='cao',<br>
outdir='./output/',<br>
fildyn='dyn.G'<br>
epsil = .true.<br>
tr2_ph=1.0d-14,<br>
/<br>
0 0 0</div>
<div style="caret-color: rgb(0, 0, 0); font-family: Helvetica;
font-size: 12px; font-weight: normal; orphans: auto; widows:
auto">
<br>
</div>
<div style="caret-color: rgb(0, 0, 0); font-family: Helvetica;
font-size: 12px; font-weight: normal; orphans: auto; widows:
auto">
input to dynmat.x<br>
---------------------------------------</div>
<div style="caret-color: rgb(0, 0, 0); font-family: Helvetica;
font-size: 12px; font-weight: normal; orphans: auto; widows:
auto">
&input<br>
fildyn = 'dyn.G'<br>
asr = 'simple'<br>
q(1)=0.0, q(2)=0.0, q(3)=1.0<br>
/</div>
<br>
</div>
<br>
<fieldset class="mimeAttachmentHeader"></fieldset>
<br>
<pre wrap="">_______________________________________________
users mailing list
<a class="moz-txt-link-abbreviated" href="mailto:users@lists.quantum-espresso.org">users@lists.quantum-espresso.org</a>
<a class="moz-txt-link-freetext" href="https://lists.quantum-espresso.org/mailman/listinfo/users">https://lists.quantum-espresso.org/mailman/listinfo/users</a></pre>
</blockquote>
<p><br>
</p>
</body>
</html>