<div dir="auto">There is some issue in your bilayer </div><div dir="auto">Try to make a new bilayer with CRYSTAL SOFTWARE </div><div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Thu, 3 Aug 2023 at 03:55, Giuseppe Scala <g.scala2@campus.unimib.it> wrote:<br></div><blockquote class="gmail_quote" style="margin: 0px 0px 0px 0.8ex; border-left-width: 1px; border-left-style: solid; padding-left: 1ex; border-left-color: rgb(204, 204, 204);"><div dir="auto">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:<div dir="auto"><br><div dir="auto">##############</div><div dir="auto"><br></div><div dir="auto">&CONTROL</div><div dir="auto">calculation ='scf'</div><div dir="auto">prefix = 'MoS2_BI'</div><div dir="auto">pseudo_dir = '/home/Giuseppe'</div><div dir="auto">outdir = './'</div><div dir="auto">tstress  = .true.</div><div dir="auto">tprnfor = .true.</div><div dir="auto">forc_conv_thr = 1.0d-3</div><div dir="auto">nstep = 200</div><div dir="auto">/</div><div dir="auto"><br></div><div dir="auto">     &SYSTEM</div><div dir="auto">ibrav = 4</div><div dir="auto">a = 3.1429635960300000000</div><div dir="auto">c = 20.006663000000000000</div><div dir="auto">nat = 6</div><div dir="auto">ntyp = 2</div><div dir="auto">ecutwfc = 60</div><div dir="auto">ecutrho = 240</div><div dir="auto">occupations = 'smearing'</div><div dir="auto">degauss = 0.01</div><div dir="auto">/</div><div dir="auto"><br></div><div dir="auto">    &ELECTRONS</div><div dir="auto">mixing_beta = 0.7</div><div dir="auto">conv_thr = 1.0d-8</div><div dir="auto">diagonalization =  'david'</div><div dir="auto">electron_maxstep=200</div><div dir="auto">/</div><div dir="auto"><br></div><div dir="auto">    &IONS</div><div dir="auto">trust_radius_max = 0.2</div><div dir="auto">/</div><div dir="auto">    &CELL</div><div dir="auto">cell_dynamics = "bfgs"</div><div dir="auto">cell_dofree = '2Dxy'</div><div dir="auto">/</div><div dir="auto"><br></div><div dir="auto">ATOMIC_SPECIES</div><div dir="auto">Mo 95.9400 Mo.pbesol-spn-kjpaw_psl.1.0.0.UPF</div><div dir="auto">S 32.0660 S.pbesol-n-kjpaw_psl.1.0.0.UPF</div><div dir="auto">ATOMIC_POSITIONS {angstrom}</div><div dir="auto">Mo   1.5714815237    0.9072955423   6.6890512735</div><div dir="auto">Mo   0.0000000000    1.8145901001   13.377562735</div><div dir="auto"> S   1.5714815237    0.9072955423   11.825385701</div><div dir="auto"> S   1.5714815237    0.9072955423   14.928508325</div><div dir="auto"> S   0.0000000000    1.8145901001   8.2412092563</div><div dir="auto"> S   0.0000000000    1.8145901001   5.1381015776</div><div dir="auto"><br></div><div dir="auto">K_POINTS {automatic}</div><div dir="auto">4 4 4 0 0 0</div><div dir="auto"><br></div><div dir="auto">######### </div><div dir="auto"><br></div><div dir="auto"><br></div><div dir="auto">And here my input file for ph.x</div><div dir="auto"><br></div><div dir="auto"><br></div><div dir="auto">#######</div><div dir="auto"><br></div><div dir="auto">&inputph</div><div dir="auto">prefix = 'MoS2_BI'</div><div dir="auto">fildvscf = 'dvscf'</div><div dir="auto">ldisp = .true.</div><div dir="auto">fildyn = 'MoS2_BI.dyn'</div><div dir="auto">alpha_mix(1) = 0.3</div><div dir="auto">nq1 = 2,</div><div dir="auto">nq2 = 2,</div><div dir="auto">nq3 = 2,</div><div dir="auto">tr2_ph = 1.0d-12</div><div dir="auto">/</div><div dir="auto">########</div><div dir="auto"><br></div><div dir="auto">Is my bilayer wrong? I've come the exact same calculations for Bulk and Monolayer MoS2 and everything went smoothly</div><div dir="auto"><br></div></div></div>________________________________<wbr>_______________
<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>________________________________<wbr>_______________
<br>Quantum ESPRESSO is supported by MaX (www.max-centre.eu)
<br>users mailing list users@lists.quantum-espresso.org
<br>https://lists.quantum-espresso.org/<wbr>mailman/listinfo/users</blockquote></div></div>