<div dir="ltr">Hi all,<div><br></div><div>I'm having a lot of troubles in modeling the adsorption of molecules on surfaces. I've quite a lot of experience with QE when working with crystals made of open-shell molecules but I'm not familiar with the description of metal surfaces.</div><div><br></div><div>First of all, I'm using Qespresso v.5.1.1<br><br>I succeed in simulating 4 layers of a 111 Gold surface with Quantum Espresso (QE). However, using the same surface and unit cell, when I include a magnetic (open-shell) molecule on top of that surface, even if its far away from it, I cannot reach convergence. I'm 99% sure that the input is fine in terms of the unit cell definition (see input file attached) and the problem is just the convergence of the electronic structure. <br><br>I believe the problem is the fact that the system is highly inhomogeneous as it has a metal surface and an open-shell molecule. I take as an example the smearing, since, although it is mandatory when modeling the metal alone, it seems to hinder the description of an open-shell molecule. Originally, I used methfessel-paxton smearing with a 'degauss' of 0.2 eV. The convergence seems to be even worse as the value is increased, and changing the type of smearing doesn't help neither.</div><div><br>I've tried several tricks to try to improve electronic structure's convergence. Changing the 'mixing_beta', the 'type of mixing', removing the smearing, switching to 'cg' algorith, increasing number K points ... all of them with identical output (take Summary.out as an example of the typical evolution of the SCF energy, you'll see that is awful). <br><br>Anyone has experience on modeling such king of systems?<br><br>Any help is appreciated, thanks in advance<br></div><div><br></div><div>Dr. Sergi Vela, </div><div>University of Strasbourg, France.</div><div><br></div><div><br></div></div>