<div dir="ltr">Dear Chiara,<div><br></div><div>I've applied acoustic sum rules with my <a href="http://matdyn.in">matdyn.in</a> files in this way:</div><div><br></div><div>#############</div><div><br></div><div>&input<br>asr = 'simple'<br>flfrc = 'MoS2_BI_Dipole.fc'<br>flfrq = 'MoS2_BI_Dipole.freq'<br>q_in_band_form = .true.<br>q_in_cryst_coord = .true.<br>/<br>5<br>0.00 0.00 0.00 100 #Gamma<br>0.500 0.000 0.000 50 #M<br>0.3333 0.333 0.000 50 #K<br>0.000 0.000 0.000 100 #Gamma<br>0.500 0.000 0.000 1 #M<br></div><div><br></div><div>###############</div><div>and attached there's my matdyn.x file, my scf output file (MoS2_BI.out) and my ph.x output file(MoS2_BI_ph333.out). You can also find  a picture of the phonon dispersion relations.</div><div>kind regards,</div><div>Giuseppe</div><div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">Il giorno gio 3 ago 2023 alle ore 09:03 Cignarella Chiara via users <<a href="mailto:users@lists.quantum-espresso.org">users@lists.quantum-espresso.org</a>> ha scritto:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Dear Giuseppe,<br>
<br>
Can you include also the output files?<br>
<br>
Did you apply the acoustic sun rules with matdyn.x (or dynmat)?<br>
<br>
Best<br>
<br>
Chiara Cignarella<br>
PhD student<br>
THEOS, EPFL Lausanne<br>
<br>
> On 3 Aug 2023, at 00:26, Giuseppe Scala <<a href="mailto:g.scala2@campus.unimib.it" target="_blank">g.scala2@campus.unimib.it</a>> wrote:<br>
> <br>
> <br>
> I get negative frequencies at gamma points for Bilayer MoS2. I relaxed my structure and checked convergence with respect to energy cutoff. That's my input file for the scf calculation:<br>
> <br>
> ##############<br>
> <br>
> &CONTROL<br>
> calculation ='scf'<br>
> prefix = 'MoS2_BI'<br>
> pseudo_dir = '/home/Giuseppe'<br>
> outdir = './'<br>
> tstress = .true.<br>
> tprnfor = .true.<br>
> forc_conv_thr = 1.0d-3<br>
> nstep = 200<br>
> /<br>
> <br>
>      &SYSTEM<br>
> ibrav = 4<br>
> a = 3.1429635960300000000<br>
> c = 20.006663000000000000<br>
> nat = 6<br>
> ntyp = 2<br>
> ecutwfc = 60<br>
> ecutrho = 240<br>
> occupations = 'smearing'<br>
> degauss = 0.01<br>
> /<br>
> <br>
>     &ELECTRONS<br>
> mixing_beta = 0.7<br>
> conv_thr = 1.0d-8<br>
> diagonalization = 'david'<br>
> electron_maxstep=200<br>
> /<br>
> <br>
>     &IONS<br>
> trust_radius_max = 0.2<br>
> /<br>
>     &CELL<br>
> cell_dynamics = "bfgs"<br>
> cell_dofree = '2Dxy'<br>
> /<br>
> <br>
> ATOMIC_SPECIES<br>
> Mo 95.9400 Mo.pbesol-spn-kjpaw_psl.1.0.0.UPF<br>
> S 32.0660 S.pbesol-n-kjpaw_psl.1.0.0.UPF<br>
> ATOMIC_POSITIONS {angstrom}<br>
> Mo 1.5714815237 0.9072955423 6.6890512735<br>
> Mo 0.0000000000 1.8145901001 13.377562735<br>
>  S 1.5714815237 0.9072955423 11.825385701<br>
>  S 1.5714815237 0.9072955423 14.928508325<br>
>  S 0.0000000000 1.8145901001 8.2412092563<br>
>  S 0.0000000000 1.8145901001 5.1381015776<br>
> <br>
> K_POINTS {automatic}<br>
> 4 4 4 0 0 0<br>
> <br>
> ######### <br>
> <br>
> <br>
> And here my input file for ph.x<br>
> <br>
> <br>
> #######<br>
> <br>
> &inputph<br>
> prefix = 'MoS2_BI'<br>
> fildvscf = 'dvscf'<br>
> ldisp = .true.<br>
> fildyn = 'MoS2_BI.dyn'<br>
> alpha_mix(1) = 0.3<br>
> nq1 = 2,<br>
> nq2 = 2,<br>
> nq3 = 2,<br>
> tr2_ph = 1.0d-12<br>
> /<br>
> ########<br>
> <br>
> Is my bilayer wrong? I've come the exact same calculations for Bulk and Monolayer MoS2 and everything went smoothly<br>
> <br>
> _______________________________________________<br>
> The Quantum ESPRESSO community stands by the Ukrainian<br>
> people and expresses its concerns about the devastating<br>
> effects that the Russian military offensive has on their<br>
> country and on the free and peaceful scientific, cultural,<br>
> and economic cooperation amongst peoples<br>
> _______________________________________________<br>
> Quantum ESPRESSO is supported by MaX (<a href="http://www.max-centre.eu" rel="noreferrer" target="_blank">www.max-centre.eu</a>)<br>
> users mailing list <a href="mailto:users@lists.quantum-espresso.org" target="_blank">users@lists.quantum-espresso.org</a><br>
> <a href="https://lists.quantum-espresso.org/mailman/listinfo/users" rel="noreferrer" target="_blank">https://lists.quantum-espresso.org/mailman/listinfo/users</a><br>
_______________________________________________<br>
The Quantum ESPRESSO community stands by the Ukrainian<br>
people and expresses its concerns about the devastating<br>
effects that the Russian military offensive has on their<br>
country and on the free and peaceful scientific, cultural,<br>
and economic cooperation amongst peoples<br>
_______________________________________________<br>
Quantum ESPRESSO is supported by MaX (<a href="http://www.max-centre.eu" rel="noreferrer" target="_blank">www.max-centre.eu</a>)<br>
users mailing list <a href="mailto:users@lists.quantum-espresso.org" target="_blank">users@lists.quantum-espresso.org</a><br>
<a href="https://lists.quantum-espresso.org/mailman/listinfo/users" rel="noreferrer" target="_blank">https://lists.quantum-espresso.org/mailman/listinfo/users</a></blockquote></div>