Skip to Content.
Sympa Menu

forum - Re: Re: [abinit-forum] TiO2 lattice constant optimization

forum@abinit.org

Subject: The ABINIT Users Mailing List ( CLOSED )

List archive

Re: Re: [abinit-forum] TiO2 lattice constant optimization


Chronological Thread 
  • From: "Anglade Pierre-Matthieu" <anglade@gmail.com>
  • To: forum@abinit.org
  • Subject: Re: Re: [abinit-forum] TiO2 lattice constant optimization
  • Date: Thu, 26 Apr 2007 13:57:18 +0200
  • Domainkey-signature: a=rsa-sha1; c=nofws; d=gmail.com; s=beta; h=received:message-id:date:from:to:subject:in-reply-to:mime-version:content-type:content-transfer-encoding:content-disposition:references; b=je+C5iLL0gMTS+thwMnc64KMJ6MKk+Oru5pA81XCr0Z5d9rWHb4p0iKXsrWbwgo8WUVa6dDJdwdeoCG9cK/euEWA1NXjYoebAMzMGb8ye6meP9C1o0pBg+UwGH9PIbTq0e6IswET7t3YBhgD6BSYD3+eaX/i2ajILeDhee1rrKU=

Hi,

Here are a few comment and pieces of advice. Comment 3 is certainly
most important for your present problem.

1) I can't tell you wether the pseudopotentials are OK or not. This is
you to discover. Usualy, to do that people compare, for each elements
the properties obtained with the pseudo and some experimental or full
electron results. When you download pseudopotentials from the abinit
site you are warn that you use them at your own risk. Then you must
pay attention to all necessary pseuopotentials checkings.

2) ecut=120 hartree is anyway too much for your calculation. You waste
your time with that

3) You must be more carefull with the way you set your kpt mesh. It is
quite certainly not correct now. Try having a look at the abinit
tutorials or at the description of Abinit input variables about kpt.
"prtkpt" will probably be very usefull for you.

4) Be certain that 30 setps are enough to get an accurate convergence
of SCF cycles. Otherwise mind having a look at variables like iprcel
and diemac.

Regards

PMA

On 4/26/07, taoj0003@ntu.edu.sg <taoj0003@ntu.edu.sg> wrote:
My input file is shown below.

the pseduopotetial i am using are download from the Abinit website. Both the
Ti and O are Troullier-Martins pseudopotentials. is there anything wrong?
thanks


# Crystalline TiO2



kptopt 1 # Option for the automatic generation of k points,
ecut 120.0 #

toldfe 1.0d-10

acell 4.55 4.55 2.95812 angstrom
optcell 6

ntypat 2
znucl 22 8

natom 6
typat 1 1 2 2 2 2
xred
0.00000000 0.00000000 0.00000000 #Ti
0.50000000 0.50000000 0.50000000 #Ti
0.30530000 0.30530000 0.00000000 #O
-0.30530000 -0.30530000 0.00000000 #O
0.80530000 0.19470000 0.50000000 #O
0.19470000 0.80530000 0.50000000 #O


# 0.00000000 0.00000000 0.00000000 #Ti
# 2.29686500 2.29686500 1.47906000 #Ti
# 1.40246577 1.40246577 0.00000000 #O
# -1.40246577 -1.40246577 0.00000000 #O
# 3.69933077 0.89439923 1.47906000 #O
# 0.89439923 3.69933077 1.47906000 #O


nstep 30



--
Pierre-Matthieu Anglade



Archive powered by MHonArc 2.6.16.

Top of Page