<html>
<head>
<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>
</head>
<body dir="ltr">
<div style="font-family: Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
I am using v.6.3MaX</div>
<div style="font-family: Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<br>
</div>
<div style="font-family: Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
Casey</div>
<div id="appendonsend"></div>
<hr style="display:inline-block;width:98%" tabindex="-1">
<div id="divRplyFwdMsg" dir="ltr"><font face="Calibri, sans-serif" style="font-size:11pt" color="#000000"><b>From:</b> users <users-bounces@lists.quantum-espresso.org> on behalf of Pietro Davide Delugas <pdelugas@sissa.it><br>
<b>Sent:</b> Wednesday, March 27, 2019 11:20 AM<br>
<b>To:</b> users@lists.quantum-espresso.org<br>
<b>Subject:</b> Re: [QE-users] incorrect LO-TO splitting for CaO and MgO</font>
<div> </div>
</div>
<div style="background-color:#FFFFFF">
<div class="x_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"><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="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="font-family:Helvetica; font-size:12px; font-weight:normal; orphans:auto; widows:auto">
<br>
</div>
<div style="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="font-family:Helvetica; font-size:12px; font-weight:normal; orphans:auto; widows:auto">
<b><br>
</b></div>
<div style="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="font-family:Helvetica; font-size:12px; font-weight:normal; orphans:auto; widows:auto">
<br>
</div>
<div style="font-family:Helvetica; font-size:12px; font-weight:normal; orphans:auto; widows:auto">
Thanks in advance,</div>
<div style="font-family:Helvetica; font-size:12px; font-weight:normal; orphans:auto; widows:auto">
Casey Brock</div>
<div style="font-family:Helvetica; font-size:12px; font-weight:normal; orphans:auto; widows:auto">
<br>
</div>
<div style="font-family:Helvetica; font-size:12px; font-weight:normal; orphans:auto; widows:auto">
input to pw.x</div>
<div style="font-family:Helvetica; font-size:12px; font-weight:normal; orphans:auto; widows:auto">
-------------------------------------------</div>
<div style="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="font-family:Helvetica; font-size:12px; font-weight:normal; orphans:auto; widows:auto">
<br>
</div>
<div style="font-family:Helvetica; font-size:12px; font-weight:normal; orphans:auto; widows:auto">
input to ph.x</div>
<div style="font-family:Helvetica; font-size:12px; font-weight:normal; orphans:auto; widows:auto">
---------------------------------------</div>
<div style="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="font-family:Helvetica; font-size:12px; font-weight:normal; orphans:auto; widows:auto">
<br>
</div>
<div style="font-family:Helvetica; font-size:12px; font-weight:normal; orphans:auto; widows:auto">
input to dynmat.x<br>
---------------------------------------</div>
<div style="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="x_mimeAttachmentHeader"></fieldset> <br>
<pre>_______________________________________________
users mailing list
<a class="x_moz-txt-link-abbreviated" href="mailto:users@lists.quantum-espresso.org">users@lists.quantum-espresso.org</a>
<a class="x_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>
</div>
</body>
</html>