<div dir="ltr"><div>Hi, <br>I am facing a problem regarding the phase transition for PbS, moreover I am following the research article: "<u><span style="font-family:-apple-system,BlinkMacSystemFont,"Segoe UI",Roboto,Oxygen-Sans,Ubuntu,Cantarell,"Helvetica Neue",sans-serif;font-size:16px">Liu, S., Guo, M., Yu, S. </span><i style="box-sizing:inherit;font-family:-apple-system,BlinkMacSystemFont,"Segoe UI",Roboto,Oxygen-Sans,Ubuntu,Cantarell,"Helvetica Neue",sans-serif;font-size:16px">et al.</i><span style="font-family:-apple-system,BlinkMacSystemFont,"Segoe UI",Roboto,Oxygen-Sans,Ubuntu,Cantarell,"Helvetica Neue",sans-serif;font-size:16px"> Hidden pressure-stabilized lead reservoirs in Earth’s mantle. </span><i style="box-sizing:inherit;font-family:-apple-system,BlinkMacSystemFont,"Segoe UI",Roboto,Oxygen-Sans,Ubuntu,Cantarell,"Helvetica Neue",sans-serif;font-size:16px">Nat Commun</i><span style="font-family:-apple-system,BlinkMacSystemFont,"Segoe UI",Roboto,Oxygen-Sans,Ubuntu,Cantarell,"Helvetica Neue",sans-serif;font-size:16px"> </span><span style="font-weight:bolder;box-sizing:inherit;font-family:-apple-system,BlinkMacSystemFont,"Segoe UI",Roboto,Oxygen-Sans,Ubuntu,Cantarell,"Helvetica Neue",sans-serif;font-size:16px">17</span></u><span style="font-family:-apple-system,BlinkMacSystemFont,"Segoe UI",Roboto,Oxygen-Sans,Ubuntu,Cantarell,"Helvetica Neue",sans-serif;font-size:16px"><u>, 2913 (2026). <a href="https://doi.org/10.1038/s41467-026-69772-8">https://doi.org/10.1038/s41467-026-69772-8</a></u> "</span> to check the initial correctness, i am using same pseudopotential and also got same values of all the phases for atomic positions and lattice parameters as they got at high pressures,  and used the phonopy-qha similar to what they used, so my question is that what i might be doing? Because the PM3M phase is metallic, is that the reason why I am getting an incorrect value?  At 20 GPa, they report a phase transition at 2600 K, while I get it at 380 K, which is a huge difference. </div><div>Here I am attaching the plot for the same (literature and ours), which I got, and the input file starting for PM3M:<br><br>&CONTROL<br>  calculation = 'scf'<br>  prefix = 'pbs'<br>  outdir = './out_super/'<br>  pseudo_dir = './pseudo'<br>  tprnfor = .true.<br>  tstress = .true.<br>  disk_io = 'low'<br>  wf_collect = .false.<br>/<br>&SYSTEM<br>  ibrav = 0<br>  nat = 128<br>  ntyp = 2<br>  ecutwfc = 100<br>  ecutrho = 800<br>  occupations = 'smearing'<br>  smearing = 'mv'<br>  degauss = 0.01<br>  nosym = .true.<br>/<br>&ELECTRONS<br>  conv_thr = 1.0d-9<br>  mixing_beta = 0.3<br>/<br>K_POINTS automatic<br>2 2 2 0 0 0</div></div>