<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=utf-8">
<style type="text/css" style="display:none;"><!-- P {margin-top:0;margin-bottom:0;} --></style>
</head>
<body dir="ltr">
<div id="divtagdefaultwrapper" style="font-size:14pt;color:#000000;font-family:OpenSans;" dir="ltr">
<p>Thanks Vahid! This is actually a nice hint. I didn't relax the lattice with QE since I got the relaxed structure from someone else (collaboration) and minor differences between codes had - in my experience (3D solids) - never such a big influence on the
phonons modes. But I will have a look at the stresses...</p>
<p><br>
</p>
<p>Anyways, meanwhile the phonon dispersion was also calculated using phonopy, i.e., the supercell method, and there is no imaginary mode. Dunno if the code interpolates differently or enforces the acoustic sum rule in a different way, but it is weird that
the collaborator uses phonopy with QE and gets positive modes and I don't...<br>
</p>
<p><br>
</p>
<p>Thanks again for the hint!</p>
<p><br>
</p>
<p>Thomas</p>
<p><br>
</p>
<div id="Signature">
<div name="divtagdefaultwrapper" style="font-family:Calibri,Arial,Helvetica,sans-serif; font-size:; margin:0">
<pre class="moz-signature" cols="72"><font size="3">--<br>Dr. rer. nat. Thomas Brumme<br>Theoretical chemistry<br>TU Dresden - BAR / II49<br>Helmholtzstr. 18<br>01069 Dresden<br><br>Tel: +49 (0)351 463 40844<br><br><font color="000000">email: thomas.brumme@tu-dresden.de</font></font><br></pre>
</div>
</div>
<br>
<br>
<div style="color: rgb(0, 0, 0);">
<div>
<hr tabindex="-1" style="display:inline-block; width:98%">
<div id="x_divRplyFwdMsg" dir="ltr"><font style="font-size:11pt" face="Calibri, sans-serif" color="#000000"><b>Von:</b> users <users-bounces@lists.quantum-espresso.org> im Auftrag von Vahid Askarpour <vh261281@dal.ca><br>
<b>Gesendet:</b> Mittwoch, 14. Juli 2021 17:29<br>
<b>An:</b> Quantum ESPRESSO users Forum<br>
<b>Betreff:</b> Re: [QE-users] Interpolation for 2D phonons results in one imaginary branch near Gamma</font>
<div> </div>
</div>
</div>
<font size="2"><span style="font-size:10pt;">
<div class="PlainText">Hi Thomas,<br>
<br>
This may not be directly related to your question and you may have already tried it but I found that for 2D materials, if in the vc-relax step, I ensure that the in-plane stresses are almost or exactly zero, I can make the flexural mode energies entirely positive
along the various BZ segments. <br>
<br>
Cheers,<br>
Vahid<br>
<br>
Vahid Askarpour<br>
Department of physics and atmospheric science<br>
Dalhousie University<br>
Halifax, NS<br>
CANADA<br>
<br>
> On Jul 14, 2021, at 12:00 PM, Thomas Brumme <tbrumme@msx.tu-dresden.de> wrote:<br>
> <br>
> CAUTION: The Sender of this email is not from within Dalhousie.<br>
> <br>
> Did anyone try recently the 2D cutoff technique for phonons? Is there<br>
> maybe an error in my input?<br>
> <br>
> Regards<br>
> <br>
> Thomas<br>
> <br>
> On 7/7/21 11:37 AM, Stefano Baroni wrote:<br>
>> Oh, certo: scusa, dimenticavo. S.<br>
>> <br>
>> --<br>
>> Stefano Baroni, Trieste -- swift message written and sent on the go<br>
>> <br>
>>> Il giorno 7 lug 2021, alle ore 10:42, Thomas Brumme <tbrumme@msx.tu-dresden.de> ha scritto:<br>
>>> <br>
>>> Dear all,<br>
>>> <br>
>>> I think I don't see the wood for the trees and hope that someone can help.<br>
>>> <br>
>>> I want to calculate the phonons for a 2D material using the 2D cutoff technique. The scf calculations for the well-relaxed system work without a problem and I also get no error calculating the phonons on a 5x5x1 q-point grid. The relevant portions of the
scf and phonon input are attached below. All frequencies are positive at every single q point, except the 3 acoustic branches at the Gamma point which are "only" -10 cm^-1 and -15 cm^-1. However, the interpolation, in order to get the dispersion, results in
one branch which is negative for quite some portion near Gamma - for the first half along G-K I have up to -130 cm^-1 and along G-M up to -60 cm^-1. I guess that this is the out-of-plane acoustic branch and now I wonder if this is real (i.e., the system is
really unstable) or if I make some obvious mistake. I will also check by explicitly calculating few q points in the relevant region but before I wanted to ask here.<br>
>>> <br>
>>> Any help would be appreciated!<br>
>>> <br>
>>> Regards<br>
>>> <br>
>>> Thomas Brumme<br>
>>> <br>
>>> --------------------------------------------<br>
>>> scf:<br>
>>> <br>
>>> &control<br>
>>> calculation = 'scf'<br>
>>> prefix = 'pasc12',<br>
>>> pseudo_dir = '/home/tbrumme/pseudos/pslibrary.1.0.0/rel-pbe/',<br>
>>> outdir = '/scratch/ws/0/tbrumme-pasc/',<br>
>>> verbosity = 'high',<br>
>>> wf_collect = .true.,<br>
>>> /<br>
>>> <br>
>>> &system<br>
>>> ibrav = 0,<br>
>>> nat = 14,<br>
>>> ntyp = 3,<br>
>>> ecutwfc = 70.0,<br>
>>> ecutrho = 560,<br>
>>> noncolin = .true.,<br>
>>> lspinorb = .true.,<br>
>>> vdw_corr = 'DFT-D',<br>
>>> nbnd = 96,<br>
>>> assume_isolated = '2D',<br>
>>> /<br>
>>> <br>
>>> &electrons<br>
>>> conv_thr = 1.D-14,<br>
>>> /<br>
>>> <br>
>>> ATOMIC_SPECIES<br>
>>> C 12.0 'C.rel-pbe-n-kjpaw_psl.1.0.0.UPF'<br>
>>> P 30.973761 'P.rel-pbe-n-kjpaw_psl.1.0.0.UPF'<br>
>>> As 74.921595 'As.rel-pbe-dn-kjpaw_psl.1.0.0.UPF'<br>
>>> <br>
>>> ATOMIC_POSITIONS crystal<br>
>>> ...<br>
>>> <br>
>>> CELL_PARAMETERS angstrom<br>
>>> 6.69668 0.0 0.0<br>
>>> -3.348340 5.799495 0.0<br>
>>> 0.0 0.0 15.0<br>
>>> <br>
>>> K_POINTS automatic<br>
>>> 25 25 1 0 0 0<br>
>>> <br>
>>> --------------------------------------------<br>
>>> phonons:<br>
>>> <br>
>>> Phonon dispersions for pasc12<br>
>>> &inputph<br>
>>> tr2_ph=1.0d-18,<br>
>>> prefix='pasc12',<br>
>>> ldisp=.true.,<br>
>>> nq1=5, nq2=5, nq3=1,<br>
>>> outdir='/scratch/ws/0/tbrumme-pasc/',<br>
>>> fildyn='pasc12.dyn.xml',<br>
>>> max_seconds = 83700,<br>
>>> verbosity='high',<br>
>>> recover=.true.<br>
>>> alpha_mix(1)=0.4,<br>
>>> /<br>
>>> <br>
>>> --------------------------------------------<br>
>>> q2r and matdyn:<br>
>>> <br>
>>> &input<br>
>>> fildyn='pasc12.dyn.xml', flfrc='pasc12.551.fc', loto_2d=.false., zasr='simple',<br>
>>> /<br>
>>> <br>
>>> &input<br>
>>> flfrc='pasc12.551.fc.xml', flfrq='pasc12.freq', q_in_band_form=.true., loto_2d=.false., asr='simple',<br>
>>> fldyn='pasc12.dyn.xml'<br>
>>> /<br>
>>> 4<br>
>>> 0.00000000 0.00000000 0.0 60<br>
>>> 0.33333333 0.57735027 0.0 30<br>
>>> 0.00000000 0.57735027 0.0 45<br>
>>> 0.00000000 0.00000000 0.0 1<br>
>>> <br>
>>> --<br>
>>> Dr. rer. nat. Thomas Brumme<br>
>>> Theoretical chemistry<br>
>>> TU Dresden - BAR / II49<br>
>>> Helmholtzstr. 18<br>
>>> 01069 Dresden<br>
>>> <br>
>>> Tel: +49 (0)351 463 40844<br>
>>> <br>
>>> email: thomas.brumme@tu-dresden.de<br>
>>> <br>
>>> _______________________________________________<br>
>>> Quantum ESPRESSO is supported by MaX (<a href="http://www.max-centre.eu" id="LPlnk208982" previewremoved="true">www.max-centre.eu</a>)<br>
>>> users mailing list users@lists.quantum-espresso.org<br>
>>> <a href="https://lists.quantum-espresso.org/mailman/listinfo/users" id="LPlnk454782" previewremoved="true">
https://lists.quantum-espresso.org/mailman/listinfo/users</a><br>
>> _______________________________________________<br>
>> Quantum ESPRESSO is supported by MaX (<a href="http://www.max-centre.eu" id="LPlnk917922" previewremoved="true">www.max-centre.eu</a>)<br>
>> users mailing list users@lists.quantum-espresso.org<br>
>> <a href="https://lists.quantum-espresso.org/mailman/listinfo/users" id="LPlnk381279" previewremoved="true">
https://lists.quantum-espresso.org/mailman/listinfo/users</a><br>
> <br>
> --<br>
> Dr. rer. nat. Thomas Brumme<br>
> Theoretical chemistry<br>
> TU Dresden - BAR / II49<br>
> Helmholtzstr. 18<br>
> 01069 Dresden<br>
> <br>
> Tel: +49 (0)351 463 40844<br>
> <br>
> email: thomas.brumme@tu-dresden.de<br>
> <br>
> _______________________________________________<br>
> Quantum ESPRESSO is supported by MaX (<a href="http://www.max-centre.eu" id="LPlnk77439" previewremoved="true">www.max-centre.eu</a>)<br>
> users mailing list users@lists.quantum-espresso.org<br>
> <a href="https://lists.quantum-espresso.org/mailman/listinfo/users" id="LPlnk199620" previewremoved="true">
https://lists.quantum-espresso.org/mailman/listinfo/users</a><br>
<br>
_______________________________________________<br>
Quantum ESPRESSO is supported by MaX (<a href="http://www.max-centre.eu" id="LPlnk541386" previewremoved="true">www.max-centre.eu</a>)<br>
users mailing list users@lists.quantum-espresso.org<br>
<a href="https://lists.quantum-espresso.org/mailman/listinfo/users" id="LPlnk666092" previewremoved="true">https://lists.quantum-espresso.org/mailman/listinfo/users</a></div>
</span></font></div>
</div>
</body>
</html>