[Pw_forum] The total forces of TiO2(101) surface do not converge

Giuseppe Mattioli giuseppe.mattioli at ism.cnr.it
Fri Sep 8 11:20:04 CEST 2017


Dear Andrey
I would be very surprised if your structure converged. :-)
I've opened your file with xcrysden and the slab seems anything but  
anatase (101).
Check the structure, first. Then

ecutwfc = 85,
ecutrho=850

It is not clear whether your PPs are norm conserving or ultrasoft  
ones. If they are NC, then you do not need to specify ecutrho, because  
it is automatically fixed at the ecutwfc*4 value. If they are US, then  
ecutwfc=85.0 is an extremely high cutoff for wavefunction, likely  
leading to a huge waste of time.

  cell_dofree = 'xy',

It is my opinion, but I suppose that when you want to simulate a  
surface with a slab, you should not optimize the xy lattice  
parameters, because in "real systems" they are generally kept fixed at  
the underlying bulk values.

HTH
Giuseppe

Quoting Andrey Chibisov <andreichibisov at yandex.ru>:

> Dear Colleagues.
> I am studing the TiO2(101) surface. But during the structure relaxation the
> total energy do not converge.
> I attach the input file.
>
>  
> -- 
> Best regards,
> Andrey Chibisov. Ph.D.
> Numerical method of mathematical physics Laboratory,
> Computational Center, Russian Academy of Sciences.
> Khabarovsk, Russia
> Web page: https://www.researchgate.net/profile/A_Chibisov
> http://ru.linkedin.com/pub/andrey-chibisov/55/253/986/en
>  


-- 
********************************************************
- Article premier - Les hommes naissent et demeurent
libres et égaux en droits. Les distinctions sociales
ne peuvent être fondées que sur l'utilité commune
- Article 2 - Le but de toute association politique
est la conservation des droits naturels et
imprescriptibles de l'homme. Ces droits sont la liberté,
la propriété, la sûreté et la résistance à l'oppression.
********************************************************

    Giuseppe Mattioli
    CNR - ISTITUTO DI STRUTTURA DELLA MATERIA
    v. Salaria Km 29,300 - C.P. 10
    I 00015 - Monterotondo Stazione (RM)
    Tel + 39 06 90672836 - Fax +39 06 90672316
    E-mail: <giuseppe.mattioli at ism.cnr.it>
    http://www.ism.cnr.it/en/staff/giuseppe-mattioli/
    ResearcherID: F-6308-2012





More information about the users mailing list