forum@abinit.org

Subject: The ABINIT Users Mailing List ( CLOSED )

List archive

Re: [abinit-forum] ecut for the RF calculation of BaTiO3

From: Xavier Gonze <xavier.gonze@uclouvain.be>
To: forum@abinit.org
Subject: Re: [abinit-forum] ecut for the RF calculation of BaTiO3
Date: Wed, 4 Mar 2009 18:30:04 +0100
Domainkey-signature: a=rsa-sha1; c=nofws; d=uclouvain.be; h=message-id:from: to:in-reply-to:content-type:content-transfer-encoding: mime-version:subject:date:references; q=dns; s=selucl; b=aqU60dv TZ7gEcPcxTLQOqqmolMcjQ+/T3dL2wiQ205hDO+NAYe3bPU6EC91CkxS3RlphF+4 oNAKnLJ6gA4dM3bYmC1lCvwiS4Xj3CXo25oHo96HPeZj4XJ8k/+484dLIbd/kgMA PPBvkjK8abXr1RZArCvGp35u3sLUOirAeUQM=

Dear Takeshi Nishimatsu,

Please, see Philippe Ghosez's thesis, precisely on BaTiO3 :
http://www.phythema.ulg.ac.be/Download/download.php
You will find interesting information about the cut-off energy convergence in an appendix
of this thesis.

Also, please note that you should compare RF and frozen-phonon
computations at exactly the same cut off energy, and you will find
that, PROVIDED you do properly your frozen-phonon finite-difference computation, that they agree
with more than 5 digits ...

Wish you the best,
Xavier Gonze

On 04 Mar 2009, at 03:06, Takeshi NISHIMATSU wrote:

Hi,

I'm testing the RF calculations of BaTiO3 with abinit-5.6.4,
Teter's extended norm-conserving LDA pseudopotentials,
8x8x8 k-points and input files in
https://loto.svn.sourceforge.net/svnroot/loto/feram/trunk/src/10example-PbTiO3-DDB2short/
(though they are for PbTiO3).

Below, I compared eigenvalues [eV/Angstrom^2] and eigenvectors
of force constant matrices at the Gamma and X points
between the frozen-phonon calculations and the RF calculations
with various ecut values, 35-80 [Hartree]. It seems that
ecut > 80 [Hartree] is required for accurate RF calculations
of (the soft mode of) BaTiO3.

Am I right?
Or should I use the other pseudopotentials?
Please give me your advice.

I am most grateful to the developers of ABINIT for their efforts!
-----Takeshi Nishimatsu

BaTiO3 (a_0 = 3.93889 [Angstrom])
eigenvalues [eV/Angstrom^2] and eigenvectors of the force constant matrix at the
Gamma point
a frozen-phonon calculation
ecut = 43 [Hartree]
eigenvalue Ba Ti O_I O_II O_III
1 -0.602341 0.1989 0.7306 -0.2908 -0.2908 -0.5074
2 0.038567 -0.4653 -0.4992 -0.4273 -0.4273 -0.4113
3 4.906813 -0.0000 -0.0000 0.7071 -0.7071 -0.0000
4 9.567793 0.8522 -0.4080 -0.2314 -0.2314 0.0118
5 14.546525 -0.1328 0.2249 -0.4235 -0.4235 0.7571

RF calculations
ecut = 35 [Hartree]
1 -0.456122 -0.128 0.000 -0.001 -0.635 0.000 -0.004 0.573 -0.000 0.002
0.355 -0.000 0.002 0.355 -0.000 0.004
2 -0.456122 -0.001 -0.000 0.128 -0.004 -0.000 0.635 0.004 0.000 -0.355
0.002 0.000 -0.355 0.002 0.000 -0.573
3 -0.456096 -0.000 -0.128 -0.000 -0.000 -0.635 -0.000 0.000 0.355 0.000
0.000 0.573 0.000 0.000 0.355 0.000
4 0.195207 -0.000 0.502 -0.000 -0.000 0.606 -0.000 -0.000 0.374 -0.000
-0.000 0.320 -0.000 -0.000 0.374 -0.000
5 0.195223 -0.502 -0.000 0.002 -0.606 -0.000 0.003 -0.320 -0.000 0.002
-0.374 -0.000 0.002 -0.374 -0.000 0.002
6 0.195223 -0.002 -0.000 -0.502 -0.003 -0.000 -0.606 -0.002 -0.000 -0.374
-0.002 -0.000 -0.374 -0.002 -0.000 -0.320
7 5.003634 -0.000 -0.000 0.000 -0.000 -0.000 -0.000 0.000 0.000 -0.707
-0.000 0.000 0.707 0.000 -0.000 0.000
8 5.003634 -0.000 -0.000 -0.000 0.000 -0.000 -0.000 -0.000 0.000 -0.000
0.707 0.000 0.000 -0.707 -0.000 0.000
9 5.003634 -0.000 0.000 -0.000 0.000 0.000 0.000 0.000 -0.707 -0.000
0.000 -0.000 0.000 -0.000 0.707 -0.000
10 9.692868 0.000 -0.845 0.000 -0.000 0.417 -0.000 0.000 0.237 -0.000
-0.000 -0.019 -0.000 -0.000 0.237 0.000
11 9.692879 0.845 0.000 -0.003 -0.417 -0.000 0.002 0.019 -0.000 0.001
-0.237 0.000 0.001 -0.237 -0.000 -0.000
12 9.692879 -0.003 -0.000 -0.845 0.002 0.000 0.417 -0.000 0.000 0.237
0.001 -0.000 0.237 0.001 0.000 -0.019
13 14.637381 -0.137 0.000 0.002 0.236 -0.000 -0.003 0.754 0.000 0.005
-0.422 -0.000 0.005 -0.422 0.000 -0.009
14 14.637381 -0.002 0.000 -0.137 0.003 -0.000 0.236 0.009 0.000 -0.422
-0.005 -0.000 -0.422 -0.005 0.000 0.754
15 14.637381 0.000 0.137 0.000 -0.000 -0.236 -0.000 -0.000 0.422 0.000
0.000 -0.754 0.000 0.000 0.422 -0.000

ecut = 40 [Hartree]
1 -0.994634 0.237 -0.000 0.001 0.788 -0.000 0.002 -0.460 0.000 -0.001
-0.237 0.000 -0.001 -0.237 0.000 -0.001
2 -0.994634 -0.001 -0.000 0.237 -0.002 -0.000 0.788 0.001 0.000 -0.237
0.001 0.000 -0.237 0.001 0.000 -0.460
3 -0.994627 -0.000 -0.237 -0.000 -0.000 -0.788 -0.000 0.000 0.237 0.000
0.000 0.460 0.000 0.000 0.237 0.000
4 -0.028809 0.436 -0.000 -0.005 0.416 -0.000 -0.005 0.466 -0.000 -0.006
0.458 -0.000 -0.006 0.458 -0.000 -0.006
5 -0.028809 -0.005 0.000 -0.436 -0.005 0.000 -0.416 -0.006 0.000 -0.458
-0.006 0.000 -0.458 -0.006 0.000 -0.466
6 -0.028807 0.000 0.436 0.000 0.000 0.416 0.000 0.000 0.458 0.000
0.000 0.466 0.000 0.000 0.458 0.000
7 4.869130 0.000 0.000 -0.000 0.000 0.000 -0.000 -0.000 -0.707 0.003
-0.000 -0.000 -0.003 0.000 0.707 -0.000
8 4.869130 0.000 -0.000 -0.000 -0.000 -0.000 0.000 -0.000 0.003 0.707
0.000 0.000 -0.707 -0.000 -0.003 -0.000
9 4.869130 0.000 0.000 -0.000 -0.000 0.000 0.000 0.000 -0.000 -0.000
0.707 -0.000 0.000 -0.707 0.000 0.000
10 9.492791 -0.000 -0.858 0.000 0.000 0.393 -0.000 -0.000 0.234 -0.000
0.000 -0.009 -0.000 0.000 0.234 0.000
11 9.492800 -0.858 0.000 0.002 0.393 -0.000 -0.001 -0.009 -0.000 -0.000
0.234 0.000 -0.000 0.234 -0.000 0.000
12 9.492800 0.002 0.000 0.858 -0.001 -0.000 -0.393 0.000 -0.000 -0.234
-0.000 0.000 -0.234 -0.000 -0.000 0.009
13 14.494830 0.000 0.135 -0.000 -0.000 -0.227 0.000 -0.000 0.424 -0.000
0.000 -0.756 -0.000 0.000 0.424 0.000
14 14.494831 -0.003 -0.000 -0.135 0.005 0.000 0.227 0.015 -0.000 -0.424
-0.009 0.000 -0.424 -0.009 -0.000 0.756
15 14.494831 -0.135 0.000 0.003 0.227 -0.000 -0.005 0.756 0.000 0.009
-0.424 -0.000 0.009 -0.424 0.000 -0.015

ecut = 45 [Hartree]
1 -0.400629 -0.000 -0.185 -0.000 -0.000 -0.709 -0.000 0.000 0.311 0.000
0.000 0.520 0.000 0.000 0.311 0.000
2 -0.400612 0.185 -0.000 0.008 0.708 -0.000 0.029 -0.519 0.000 -0.013
-0.311 0.000 -0.013 -0.311 0.000 -0.021
3 -0.400612 -0.008 -0.000 0.185 -0.029 -0.000 0.708 0.021 0.000 -0.311
0.013 0.000 -0.311 0.013 0.000 -0.519
4 0.046698 -0.474 -0.000 -0.001 -0.527 -0.000 -0.001 -0.391 -0.000 -0.001
-0.415 -0.000 -0.001 -0.415 -0.000 -0.001
5 0.046698 -0.001 0.000 0.474 -0.001 0.000 0.527 -0.001 0.000 0.415
-0.001 0.000 0.415 -0.001 0.000 0.391
6 0.046717 0.000 -0.474 0.000 0.000 -0.527 0.000 0.000 -0.415 0.000
0.000 -0.391 0.000 0.000 -0.415 0.000
7 4.906948 -0.000 0.000 0.000 0.000 0.000 0.000 -0.000 0.000 -0.000
-0.707 -0.000 0.000 0.707 -0.000 -0.000
8 4.906948 -0.000 0.000 -0.000 -0.000 0.000 0.000 0.000 0.000 0.707
-0.000 -0.000 -0.707 0.000 -0.000 -0.000
9 4.906948 -0.000 0.000 -0.000 -0.000 0.000 -0.000 0.000 0.707 -0.000
0.000 -0.000 0.000 -0.000 -0.707 -0.000
10 9.613126 -0.851 -0.000 0.003 0.414 0.000 -0.002 -0.012 0.000 -0.001
0.228 -0.000 -0.001 0.228 0.000 0.000
11 9.613126 0.003 0.000 0.851 -0.002 -0.000 -0.414 0.000 -0.000 -0.228
-0.001 0.000 -0.228 -0.001 -0.000 0.012
12 9.613127 -0.000 0.851 -0.000 0.000 -0.414 0.000 -0.000 -0.228 0.000
0.000 0.012 0.000 0.000 -0.228 -0.000
13 14.561896 0.000 -0.130 -0.000 -0.000 0.220 0.000 -0.000 -0.423 -0.000
0.000 0.759 -0.000 0.000 -0.423 0.000
14 14.561899 -0.130 -0.000 0.002 0.220 0.000 -0.003 0.759 -0.000 0.006
-0.423 0.000 0.006 -0.423 -0.000 -0.012
15 14.561899 0.002 -0.000 0.130 -0.003 0.000 -0.220 -0.012 -0.000 0.423
0.006 0.000 0.423 0.006 -0.000 -0.759

ecut = 50 [Hartree]
1 -0.616130 0.203 -0.000 -0.004 0.735 -0.000 -0.015 -0.505 0.000 0.006
-0.286 0.000 0.006 -0.286 0.000 0.011
2 -0.616130 -0.004 -0.000 -0.203 -0.015 -0.000 -0.735 0.011 0.000 0.286
0.006 0.000 0.286 0.006 0.000 0.505
3 -0.616122 0.000 0.203 -0.000 0.000 0.735 -0.000 -0.000 -0.286 0.000
-0.000 -0.505 0.000 -0.000 -0.286 0.000
4 0.037663 -0.463 -0.000 0.001 -0.493 -0.000 0.002 -0.416 -0.000 0.001
-0.430 -0.000 0.001 -0.430 -0.000 0.001
5 0.037663 0.001 0.000 0.463 0.002 0.000 0.493 0.001 0.000 0.430
0.001 0.000 0.430 0.001 0.000 0.416
6 0.037667 0.000 -0.463 0.000 0.000 -0.493 0.000 0.000 -0.430 0.000
0.000 -0.416 0.000 0.000 -0.430 0.000
7 4.925377 -0.000 -0.000 -0.000 -0.000 0.000 0.000 0.000 0.707 0.000
0.000 0.000 -0.000 -0.000 -0.707 -0.000
8 4.925377 -0.000 -0.000 0.000 -0.000 0.000 0.000 -0.000 0.000 -0.000
-0.707 0.000 0.000 0.707 -0.000 -0.000
9 4.925377 0.000 -0.000 0.000 0.000 0.000 0.000 -0.000 0.000 -0.707
0.000 0.000 0.707 -0.000 -0.000 0.000
10 9.575858 -0.000 0.852 0.000 0.000 -0.407 -0.000 -0.000 -0.232 -0.000
0.000 0.012 -0.000 0.000 -0.232 0.000
11 9.575867 0.023 0.000 -0.852 -0.011 -0.000 0.407 0.000 -0.000 0.232
-0.006 0.000 0.232 -0.006 -0.000 -0.012
12 9.575867 -0.852 -0.000 -0.023 0.407 0.000 0.011 -0.012 0.000 0.006
0.232 -0.000 0.006 0.232 0.000 -0.000
13 14.572737 0.000 -0.133 -0.000 -0.000 0.226 0.000 -0.000 -0.424 -0.000
0.000 0.756 -0.000 0.000 -0.424 0.000
14 14.572739 -0.133 -0.000 0.000 0.226 0.000 -0.000 0.756 -0.000 0.000
-0.424 0.000 0.000 -0.424 -0.000 -0.000
15 14.572739 0.000 -0.000 0.133 -0.000 0.000 -0.226 -0.000 -0.000 0.424
0.000 0.000 0.424 0.000 -0.000 -0.756

ecut = 55 [Hartree]
1 -0.638695 -0.001 -0.000 0.210 -0.003 -0.000 0.745 0.002 0.000 -0.278
0.001 0.000 -0.278 0.001 0.000 -0.497
2 -0.638695 0.210 -0.000 0.001 0.745 -0.000 0.003 -0.497 0.000 -0.001
-0.278 0.000 -0.001 -0.278 0.000 -0.002
3 -0.638630 -0.000 -0.210 -0.000 -0.000 -0.745 -0.000 0.000 0.278 0.000
0.000 0.497 0.000 0.000 0.278 0.000
4 0.025415 0.000 -0.459 0.000 0.000 -0.479 0.000 0.000 -0.435 0.000
0.000 -0.425 0.000 0.000 -0.435 0.000
5 0.025427 0.459 0.000 -0.002 0.480 0.000 -0.002 0.425 0.000 -0.002
0.435 0.000 -0.002 0.435 0.000 -0.001
6 0.025427 -0.002 -0.000 -0.459 -0.002 -0.000 -0.480 -0.001 -0.000 -0.435
-0.002 -0.000 -0.435 -0.002 -0.000 -0.425
7 4.917949 0.000 0.000 0.000 -0.000 0.000 0.000 0.000 -0.000 -0.000
-0.707 -0.000 0.000 0.707 0.000 -0.000
8 4.917949 -0.000 -0.000 0.000 -0.000 -0.000 -0.000 0.000 0.000 -0.707
0.000 0.000 0.707 -0.000 -0.000 -0.000
9 4.917949 -0.000 0.000 0.000 0.000 0.000 0.000 -0.000 -0.707 -0.000
0.000 -0.000 0.000 -0.000 0.707 -0.000
10 9.574672 0.000 -0.853 -0.000 -0.000 0.406 0.000 0.000 0.232 0.000
-0.000 -0.012 0.000 -0.000 0.232 -0.000
11 9.574693 0.000 -0.000 0.853 -0.000 0.000 -0.406 0.000 0.000 -0.232
-0.000 -0.000 -0.232 -0.000 0.000 0.012
12 9.574693 0.853 0.000 -0.000 -0.406 -0.000 0.000 0.012 -0.000 0.000
-0.232 0.000 0.000 -0.232 -0.000 -0.000
13 14.549221 0.000 -0.133 0.000 -0.000 0.226 -0.000 -0.000 -0.424 0.000
0.000 0.756 0.000 0.000 -0.424 -0.000
14 14.549253 -0.133 -0.000 -0.000 0.226 0.000 0.000 0.756 -0.000 -0.001
-0.424 0.000 -0.001 -0.424 -0.000 0.001
15 14.549253 0.000 -0.000 -0.133 -0.000 0.000 0.226 -0.001 -0.000 -0.424
0.001 0.000 -0.424 0.001 -0.000 0.756

ecut = 57 [Hartree]
1 -0.652137 -0.000 -0.000 0.217 -0.001 -0.000 0.753 0.001 0.000 -0.271
0.000 0.000 -0.271 0.000 0.000 -0.489
2 -0.652137 -0.217 0.000 -0.000 -0.753 0.000 -0.001 0.489 -0.000 0.000
0.271 -0.000 0.000 0.271 -0.000 0.001
3 -0.652112 0.000 0.217 0.000 0.000 0.753 0.000 -0.000 -0.271 -0.000
-0.000 -0.489 -0.000 -0.000 -0.271 -0.000
4 0.013472 0.000 -0.454 0.000 0.000 -0.468 0.000 0.000 -0.440 0.000
0.000 -0.434 0.000 0.000 -0.440 0.000
5 0.013494 0.454 0.000 0.008 0.468 0.000 0.008 0.434 0.000 0.007
0.440 0.000 0.007 0.440 0.000 0.007
6 0.013494 0.008 -0.000 -0.454 0.008 -0.000 -0.468 0.007 -0.000 -0.440
0.007 -0.000 -0.440 0.007 -0.000 -0.434
7 4.910288 -0.000 -0.000 0.000 -0.000 -0.000 -0.000 0.000 0.000 -0.707
-0.000 0.000 0.707 0.000 -0.000 0.000
8 4.910288 -0.000 -0.000 -0.000 -0.000 -0.000 0.000 0.000 0.000 0.000
-0.707 0.000 -0.000 0.707 -0.000 -0.000
9 4.910288 -0.000 0.000 0.000 0.000 0.000 0.000 0.000 -0.707 -0.000
-0.000 -0.000 0.000 0.000 0.707 -0.000
10 9.570446 0.000 -0.854 -0.000 -0.000 0.405 0.000 0.000 0.231 0.000
-0.000 -0.011 0.000 -0.000 0.231 -0.000
11 9.570457 0.854 0.000 -0.007 -0.405 -0.000 0.004 0.011 -0.000 0.002
-0.231 0.000 0.002 -0.231 -0.000 -0.000
12 9.570457 -0.007 0.000 -0.854 0.004 -0.000 0.405 -0.000 -0.000 0.231
0.002 0.000 0.231 0.002 -0.000 -0.011
13 14.542755 0.067 0.000 0.115 -0.114 -0.000 -0.194 -0.383 0.000 0.366
0.214 -0.000 0.366 0.214 0.000 -0.653
14 14.542755 0.115 0.000 -0.067 -0.194 -0.000 0.114 -0.653 0.000 -0.214
0.366 -0.000 -0.214 0.366 0.000 0.383
15 14.542764 -0.000 0.133 0.000 0.000 -0.225 -0.000 0.000 0.424 0.000
-0.000 -0.757 0.000 -0.000 0.424 -0.000

ecut = 60 [Hartree]
1 -0.644218 0.217 -0.000 -0.003 0.753 -0.000 -0.009 -0.489 0.000 0.003
-0.271 0.000 0.003 -0.271 0.000 0.006
2 -0.644218 0.003 0.000 0.217 0.009 0.000 0.753 -0.006 -0.000 -0.271
-0.003 -0.000 -0.271 -0.003 -0.000 -0.489
3 -0.644203 -0.000 -0.217 0.000 -0.000 -0.753 0.000 0.000 0.271 -0.000
0.000 0.489 -0.000 0.000 0.271 -0.000
4 0.014311 0.000 0.454 -0.000 0.000 0.468 -0.000 0.000 0.440 -0.000
0.000 0.434 -0.000 0.000 0.440 -0.000
5 0.014318 -0.454 0.000 0.019 -0.467 0.000 0.020 -0.433 0.000 0.019
-0.439 0.000 0.019 -0.439 0.000 0.018
6 0.014318 0.019 0.000 0.454 0.020 0.000 0.467 0.018 0.000 0.439
0.019 0.000 0.439 0.019 0.000 0.433
7 4.912821 0.000 0.000 0.000 0.000 0.000 -0.000 -0.000 -0.001 0.707
0.000 -0.000 -0.707 -0.000 0.001 0.000
8 4.912821 0.000 -0.000 -0.000 -0.000 -0.000 0.000 -0.000 0.000 0.000
-0.707 0.000 -0.000 0.707 -0.000 -0.000
9 4.912821 -0.000 -0.000 0.000 -0.000 -0.000 0.000 0.000 0.707 0.001
0.000 0.000 -0.001 -0.000 -0.707 0.000
10 9.568857 -0.000 0.854 0.000 0.000 -0.405 -0.000 -0.000 -0.231 -0.000
0.000 0.011 -0.000 0.000 -0.231 0.000
11 9.568865 0.001 0.000 -0.854 -0.001 -0.000 0.405 0.000 -0.000 0.231
-0.000 0.000 0.231 -0.000 -0.000 -0.011
12 9.568865 -0.854 -0.000 -0.001 0.405 0.000 0.001 -0.011 0.000 0.000
0.231 -0.000 0.000 0.231 0.000 -0.000
13 14.545352 -0.133 0.000 0.000 0.225 -0.000 -0.000 0.757 0.000 0.001
-0.424 -0.000 0.001 -0.424 0.000 -0.002
14 14.545352 0.000 -0.000 0.133 -0.000 0.000 -0.225 -0.002 -0.000 0.424
0.001 0.000 0.424 0.001 -0.000 -0.757
15 14.545354 0.000 0.133 0.000 -0.000 -0.225 -0.000 -0.000 0.424 0.000
0.000 -0.757 0.000 0.000 0.424 -0.000

ecut = 61 [Hartree]
1 -0.643811 0.217 0.000 0.005 0.753 0.000 0.017 -0.489 -0.000 -0.006
-0.271 -0.000 -0.006 -0.271 -0.000 -0.011
2 -0.643811 0.005 -0.000 -0.217 0.017 -0.000 -0.753 -0.011 0.000 0.271
-0.006 0.000 0.271 -0.006 0.000 0.489
3 -0.643778 0.000 -0.217 0.000 0.000 -0.753 0.000 -0.000 0.271 -0.000
-0.000 0.489 -0.000 -0.000 0.271 -0.000
4 0.014193 -0.000 0.454 -0.000 -0.000 0.467 -0.000 -0.000 0.440 -0.000
-0.000 0.434 -0.000 -0.000 0.440 -0.000
5 0.014197 -0.454 -0.000 0.001 -0.468 -0.000 0.001 -0.434 -0.000 0.001
-0.440 -0.000 0.001 -0.440 -0.000 0.001
6 0.014197 -0.001 -0.000 -0.454 -0.001 -0.000 -0.468 -0.001 -0.000 -0.440
-0.001 -0.000 -0.440 -0.001 -0.000 -0.434
7 4.912773 -0.000 -0.000 -0.000 -0.000 -0.000 -0.000 0.000 0.000 0.000
0.707 0.000 -0.000 -0.707 -0.000 0.000
8 4.912773 -0.000 -0.000 0.000 -0.000 -0.000 -0.000 0.000 0.000 -0.707
0.000 0.000 0.707 -0.000 -0.000 0.000
9 4.912773 -0.000 -0.000 0.000 -0.000 -0.000 0.000 0.000 0.707 0.000
-0.000 0.000 -0.000 0.000 -0.707 0.000
10 9.568942 0.000 -0.854 0.000 -0.000 0.405 -0.000 0.000 0.231 -0.000
-0.000 -0.011 -0.000 -0.000 0.231 0.000
11 9.568962 0.854 0.000 -0.004 -0.405 -0.000 0.002 0.011 -0.000 0.001
-0.231 0.000 0.001 -0.231 -0.000 -0.000
12 9.568962 -0.004 -0.000 -0.854 0.002 0.000 0.405 -0.000 0.000 0.231
0.001 -0.000 0.231 0.001 0.000 -0.011
13 14.545186 0.000 0.133 -0.000 -0.000 -0.225 0.000 -0.000 0.424 -0.000
0.000 -0.757 -0.000 0.000 0.424 0.000
14 14.545195 -0.133 0.000 0.000 0.225 -0.000 -0.000 0.757 0.000 0.001
-0.424 -0.000 0.001 -0.424 0.000 -0.001
15 14.545195 0.000 0.000 0.133 -0.000 -0.000 -0.225 -0.001 0.000 0.424
0.001 -0.000 0.424 0.001 0.000 -0.757

ecut = 62 [Hartree]
1 -0.642569 0.011 0.000 -0.217 0.037 0.000 -0.752 -0.024 -0.000 0.271
-0.013 -0.000 0.271 -0.013 -0.000 0.489
2 -0.642569 -0.217 0.000 -0.011 -0.752 0.000 -0.037 0.489 -0.000 0.013
0.271 -0.000 0.013 0.271 -0.000 0.024
3 -0.642548 0.000 0.217 0.000 0.000 0.753 0.000 -0.000 -0.271 -0.000
-0.000 -0.489 -0.000 -0.000 -0.271 -0.000
4 0.014140 -0.000 0.454 -0.000 -0.000 0.467 -0.000 -0.000 0.440 -0.000
-0.000 0.434 -0.000 -0.000 0.440 -0.000
5 0.014148 0.454 0.000 -0.004 0.467 0.000 -0.004 0.434 0.000 -0.004
0.440 0.000 -0.004 0.440 0.000 -0.004
6 0.014148 -0.004 -0.000 -0.454 -0.004 -0.000 -0.467 -0.004 -0.000 -0.440
-0.004 -0.000 -0.440 -0.004 -0.000 -0.434
7 4.912736 -0.000 -0.000 0.000 -0.000 -0.000 0.000 0.000 0.000 -0.000
-0.707 0.000 0.000 0.707 -0.000 -0.000
8 4.912736 -0.000 0.000 0.000 -0.000 0.000 -0.000 0.000 -0.000 0.707
-0.000 -0.000 -0.707 0.000 0.000 0.000
9 4.912736 0.000 -0.000 0.000 -0.000 -0.000 -0.000 0.000 0.707 0.000
0.000 0.000 -0.000 -0.000 -0.707 0.000
10 9.568814 -0.000 -0.854 0.000 0.000 0.405 -0.000 -0.000 0.231 -0.000
0.000 -0.011 -0.000 0.000 0.231 0.000
11 9.568822 -0.854 0.000 0.012 0.405 -0.000 -0.006 -0.011 -0.000 -0.003
0.231 0.000 -0.003 0.231 -0.000 0.000
12 9.568822 0.012 0.000 0.854 -0.006 -0.000 -0.405 0.000 -0.000 -0.231
-0.003 0.000 -0.231 -0.003 -0.000 0.011
13 14.545466 0.000 0.133 -0.000 -0.000 -0.225 0.000 -0.000 0.424 -0.000
0.000 -0.757 -0.000 0.000 0.424 0.000
14 14.545472 -0.001 0.000 0.133 0.002 -0.000 -0.225 0.005 0.000 0.424
-0.003 -0.000 0.424 -0.003 0.000 -0.757
15 14.545472 -0.133 0.000 -0.001 0.225 -0.000 0.002 0.757 0.000 -0.003
-0.424 -0.000 -0.003 -0.424 0.000 0.005

ecut = 65 [Hartree]
1 -0.656336 -0.000 0.222 -0.000 -0.000 0.759 -0.000 0.000 -0.266 0.000
0.000 -0.484 0.000 0.000 -0.266 0.000
2 -0.656328 -0.217 -0.000 0.045 -0.743 -0.000 0.153 0.474 0.000 -0.054
0.260 0.000 -0.054 0.260 0.000 -0.098
3 -0.656328 0.045 0.000 0.217 0.153 0.000 0.743 -0.098 -0.000 -0.260
-0.054 -0.000 -0.260 -0.054 -0.000 -0.474
4 0.007038 -0.451 0.000 0.003 -0.458 0.000 0.003 -0.440 0.000 0.003
-0.443 0.000 0.003 -0.443 0.000 0.003
5 0.007038 -0.003 0.000 -0.451 -0.003 0.000 -0.458 -0.003 0.000 -0.443
-0.003 0.000 -0.443 -0.003 0.000 -0.440
6 0.007044 0.000 0.451 0.000 0.000 0.458 0.000 0.000 0.443 0.000
0.000 0.440 0.000 0.000 0.443 0.000
7 4.910474 0.000 0.000 -0.000 -0.000 0.000 0.000 -0.000 0.000 -0.000
-0.707 -0.000 0.000 0.707 -0.000 -0.000
8 4.910474 0.000 0.000 -0.000 0.000 0.000 0.000 -0.000 0.002 -0.707
0.000 -0.000 0.707 -0.000 -0.002 -0.000
9 4.910474 -0.000 0.000 0.000 -0.000 0.000 -0.000 0.000 0.707 0.002
0.000 -0.000 -0.002 -0.000 -0.707 0.000
10 9.563921 -0.001 0.000 -0.854 0.000 -0.000 0.404 -0.000 -0.000 0.231
0.000 0.000 0.231 0.000 -0.000 -0.011
11 9.563921 0.854 -0.000 -0.001 -0.404 0.000 0.000 0.011 0.000 0.000
-0.231 -0.000 0.000 -0.231 0.000 -0.000
12 9.563922 -0.000 -0.854 -0.000 0.000 0.404 0.000 -0.000 0.231 0.000
0.000 -0.011 0.000 0.000 0.231 -0.000
13 14.543022 -0.133 -0.000 -0.003 0.224 0.000 0.004 0.757 -0.000 -0.008
-0.424 0.000 -0.008 -0.424 -0.000 0.015
14 14.543022 -0.003 -0.000 0.133 0.004 0.000 -0.224 0.015 -0.000 0.424
-0.008 0.000 0.424 -0.008 -0.000 -0.757
15 14.543025 0.000 -0.133 -0.000 -0.000 0.224 0.000 -0.000 -0.424 -0.000
0.000 0.757 -0.000 0.000 -0.424 0.000

ecut = 70 [Hartree]
1 -0.685880 0.000 -0.221 -0.000 0.000 -0.759 -0.000 -0.000 0.265 0.000
-0.000 0.484 0.000 -0.000 0.265 0.000
2 -0.685874 -0.000 0.000 -0.221 -0.001 0.000 -0.759 0.001 -0.000 0.265
0.000 -0.000 0.265 0.000 -0.000 0.484
3 -0.685874 -0.221 -0.000 0.000 -0.759 -0.000 0.001 0.484 0.000 -0.000
0.265 0.000 -0.000 0.265 0.000 -0.001
4 0.007253 -0.001 -0.000 -0.451 -0.001 -0.000 -0.458 -0.001 -0.000 -0.443
-0.001 -0.000 -0.443 -0.001 -0.000 -0.440
5 0.007253 0.451 -0.000 -0.001 0.458 -0.000 -0.001 0.440 -0.000 -0.001
0.443 -0.000 -0.001 0.443 -0.000 -0.001
6 0.007259 0.000 0.451 -0.000 0.000 0.458 -0.000 0.000 0.443 -0.000
0.000 0.440 -0.000 0.000 0.443 -0.000
7 4.908248 -0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 -0.000
-0.707 -0.000 0.000 0.707 -0.000 -0.000
8 4.908248 -0.000 0.000 -0.000 0.000 0.000 0.000 0.000 0.001 0.707
-0.000 -0.000 -0.707 0.000 -0.001 0.000
9 4.908248 -0.000 0.000 -0.000 -0.000 0.000 0.000 0.000 0.707 -0.001
0.000 -0.000 0.001 -0.000 -0.707 -0.000
10 9.559438 -0.854 -0.000 0.018 0.403 0.000 -0.008 -0.011 0.000 -0.005
0.232 -0.000 -0.005 0.232 0.000 0.000
11 9.559438 -0.018 -0.000 -0.854 0.008 0.000 0.403 -0.000 0.000 0.232
0.005 -0.000 0.232 0.005 0.000 -0.011
12 9.559441 -0.000 0.854 -0.000 0.000 -0.404 0.000 -0.000 -0.232 0.000
0.000 0.011 0.000 0.000 -0.232 -0.000
13 14.532712 0.000 0.133 0.000 -0.000 -0.225 -0.000 -0.000 0.424 0.000
0.000 -0.756 0.000 0.000 0.424 -0.000
14 14.532715 -0.002 -0.000 0.133 0.003 0.000 -0.225 0.010 -0.000 0.424
-0.006 0.000 0.424 -0.006 -0.000 -0.756
15 14.532715 -0.133 0.000 -0.002 0.225 -0.000 0.003 0.756 0.000 -0.006
-0.424 -0.000 -0.006 -0.424 0.000 0.010

ecut = 75 [Hartree]
1 -0.685909 0.000 -0.221 -0.000 0.000 -0.759 -0.000 -0.000 0.265 0.000
-0.000 0.484 0.000 -0.000 0.265 0.000
2 -0.685903 0.001 0.000 -0.221 0.003 0.000 -0.759 -0.002 -0.000 0.265
-0.001 -0.000 0.265 -0.001 -0.000 0.484
3 -0.685903 -0.221 -0.000 -0.001 -0.759 -0.000 -0.003 0.484 0.000 0.001
0.265 0.000 0.001 0.265 0.000 0.002
4 0.007671 -0.002 0.000 -0.451 -0.002 0.000 -0.459 -0.002 0.000 -0.443
-0.002 0.000 -0.443 -0.002 0.000 -0.439
5 0.007671 -0.451 0.000 0.002 -0.459 0.000 0.002 -0.439 0.000 0.002
-0.443 0.000 0.002 -0.443 0.000 0.002
6 0.007677 0.000 0.451 0.000 0.000 0.459 0.000 0.000 0.443 0.000
0.000 0.439 0.000 0.000 0.443 0.000
7 4.909921 0.000 0.000 -0.000 0.000 0.000 0.000 -0.000 0.008 -0.707
-0.000 -0.000 0.707 0.000 -0.008 -0.000
8 4.909921 -0.000 -0.000 -0.000 0.000 -0.000 0.000 0.000 -0.000 0.000
-0.707 0.000 -0.000 0.707 0.000 -0.000
9 4.909921 -0.000 0.000 0.000 -0.000 0.000 0.000 0.000 0.707 0.008
-0.000 -0.000 -0.008 0.000 -0.707 0.000
10 9.553733 -0.854 0.000 -0.030 0.404 -0.000 0.014 -0.011 -0.000 0.008
0.231 0.000 0.008 0.231 -0.000 -0.000
11 9.553733 -0.030 -0.000 0.854 0.014 0.000 -0.404 -0.000 0.000 -0.231
0.008 -0.000 -0.231 0.008 0.000 0.011
12 9.553736 -0.000 -0.854 -0.000 0.000 0.404 0.000 -0.000 0.232 0.000
0.000 -0.011 0.000 0.000 0.232 -0.000
13 14.542708 0.000 -0.133 0.000 -0.000 0.225 -0.000 -0.000 -0.424 0.000
0.000 0.756 0.000 0.000 -0.424 -0.000
14 14.542712 -0.133 -0.000 -0.002 0.225 0.000 0.003 0.756 -0.000 -0.005
-0.424 0.000 -0.005 -0.424 -0.000 0.009
15 14.542712 0.002 -0.000 -0.133 -0.003 0.000 0.225 -0.009 -0.000 -0.424
0.005 0.000 -0.424 0.005 -0.000 0.756

ecut = 80 [Hartree]
eigenvalues [eV/Angstrom^2] and eigenvectors of the force constant matrix at the
Gamma point
1 -0.648047 -0.221 0.000 -0.001 -0.758 0.000 -0.003 0.485 -0.000 0.001
0.267 -0.000 0.001 0.267 -0.000 0.002
2 -0.648047 -0.001 -0.000 0.221 -0.003 -0.000 0.758 0.002 0.000 -0.267
0.001 0.000 -0.267 0.001 0.000 -0.485
3 -0.648020 0.000 0.221 0.000 0.000 0.758 0.000 -0.000 -0.267 -0.000
-0.000 -0.485 -0.000 -0.000 -0.267 -0.000
4 0.008204 -0.452 0.000 -0.001 -0.460 0.000 -0.001 -0.438 0.000 -0.001
-0.442 0.000 -0.001 -0.442 0.000 -0.001
5 0.008204 -0.001 -0.000 0.452 -0.001 -0.000 0.460 -0.001 -0.000 0.442
-0.001 -0.000 0.442 -0.001 -0.000 0.438
6 0.008213 0.000 0.452 0.000 0.000 0.460 0.000 0.000 0.442 0.000
0.000 0.439 0.000 0.000 0.442 0.000
7 4.911302 0.000 0.000 -0.000 0.000 0.000 0.000 -0.000 -0.000 0.707
0.001 -0.000 -0.707 -0.001 0.000 -0.000
8 4.911302 0.000 0.000 0.000 0.000 0.000 -0.000 0.000 -0.000 -0.001
0.707 -0.000 0.001 -0.707 0.000 0.000
9 4.911302 -0.000 0.000 -0.000 0.000 0.000 -0.000 -0.000 -0.707 -0.000
-0.000 -0.000 0.000 0.000 0.707 -0.000
10 9.563699 0.000 -0.854 -0.000 -0.000 0.405 0.000 0.000 0.231 0.000
-0.000 -0.011 0.000 -0.000 0.231 -0.000
11 9.563710 0.854 0.000 0.003 -0.405 -0.000 -0.002 0.011 -0.000 -0.001
-0.231 0.000 -0.001 -0.231 -0.000 0.000
12 9.563710 -0.003 -0.000 0.854 0.002 0.000 -0.405 -0.000 0.000 -0.231
0.001 -0.000 -0.231 0.001 0.000 0.011
13 14.547952 0.000 -0.133 0.000 -0.000 0.224 -0.000 -0.000 -0.424 0.000
0.000 0.757 0.000 0.000 -0.424 -0.000
14 14.547976 -0.133 -0.000 0.001 0.224 0.000 -0.002 0.757 -0.000 0.004
-0.424 0.000 0.004 -0.424 -0.000 -0.006
15 14.547976 0.001 0.000 0.133 -0.002 -0.000 -0.224 -0.006 0.000 0.424
0.004 -0.000 0.424 0.004 0.000 -0.757

BaTiO3 (a_0 = 3.93889 [Angstrom])
eigenvalues [eV/Angstrom^2] and eigenvectors of the force constant matrix at the X
point
a frozen-phonon calculation
ecut = 43 [Hartree]
1 0.664710 -0.000 0.000 0.000 0.000 0.788 -0.000 -0.000 0.000 0.000
0.000 -0.488 0.000 -0.000 -0.375 0.000
2 0.664773 -0.000 0.000 -0.000 0.000 -0.000 -0.788 -0.000 0.000 0.000
-0.000 0.000 0.375 0.000 0.000 0.488
3 3.434234 0.000 -0.000 0.000 -0.000 -0.582 0.000 0.000 -0.000 -0.000
-0.000 -0.392 0.000 0.000 -0.713 0.000
4 3.434345 0.000 -0.000 0.000 -0.000 0.000 0.582 0.000 -0.000 0.000
0.000 0.000 0.713 -0.000 0.000 0.392
5 5.983030 -0.000 0.709 0.000 0.000 -0.000 0.000 -0.000 0.705 0.000
0.000 0.000 0.000 0.000 -0.000 -0.000
6 5.983127 0.000 -0.000 0.709 0.000 -0.000 -0.000 0.000 -0.000 0.705
-0.000 0.000 -0.000 0.000 0.000 -0.000
7 6.588437 -0.000 0.000 0.000 -0.000 -0.000 -0.000 -0.000 0.000 0.000
0.707 0.000 -0.000 -0.707 -0.000 0.000
8 6.952857 -0.000 0.000 0.000 -0.358 -0.000 -0.000 -0.000 -0.000 -0.000
-0.660 -0.000 -0.000 -0.660 0.000 0.000
9 7.809087 -0.000 0.000 -0.705 -0.000 0.000 0.000 -0.000 -0.000 0.709
-0.000 0.000 -0.000 0.000 -0.000 -0.000
10 7.809342 0.000 -0.705 -0.000 -0.000 -0.000 -0.000 0.000 0.709 0.000
-0.000 -0.000 -0.000 -0.000 0.000 0.000
11 11.984757 -0.930 -0.000 -0.000 -0.000 -0.000 -0.000 -0.368 -0.000 -0.000
-0.000 -0.000 -0.000 -0.000 -0.000 -0.000
12 12.742579 -0.000 0.000 0.000 0.000 -0.000 0.202 -0.000 -0.000 0.000
-0.000 -0.000 -0.593 0.000 0.000 0.780
13 12.742890 0.000 -0.000 0.000 -0.000 0.202 0.000 0.000 0.000 -0.000
-0.000 0.780 -0.000 0.000 -0.593 0.000
14 30.316858 -0.368 0.000 0.000 0.000 0.000 0.000 0.930 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000
15 37.059905 0.000 -0.000 -0.000 0.934 -0.000 -0.000 0.000 -0.000 0.000
-0.253 0.000 0.000 -0.253 -0.000 -0.000

RF calculations
ecut = 35
1 0.856958 -0.000 0.000 0.000 0.000 -0.754 0.000 -0.000 0.000 0.000
0.000 0.516 -0.000 0.000 0.406 -0.000
2 0.856998 -0.000 -0.000 -0.000 -0.000 -0.000 -0.754 -0.000 -0.000 0.000
-0.000 0.000 0.406 -0.000 0.000 0.516
3 3.645806 0.000 0.000 -0.000 0.000 -0.000 0.622 0.000 0.000 0.000
-0.000 -0.000 0.695 0.000 -0.000 0.361
4 3.645884 -0.000 0.000 0.000 0.000 -0.622 -0.000 -0.000 0.000 0.000
0.000 -0.361 -0.000 0.000 -0.695 -0.000
5 6.043002 0.000 0.011 0.720 -0.000 0.000 -0.000 0.000 0.011 0.694
-0.000 -0.000 -0.000 -0.000 0.000 -0.000
6 6.043002 0.000 -0.720 0.011 0.000 -0.000 0.000 0.000 -0.694 0.011
0.000 -0.000 0.000 0.000 -0.000 -0.000
7 6.682457 -0.000 0.000 0.000 0.000 0.000 0.000 -0.000 -0.000 -0.000
0.707 -0.000 0.000 -0.707 0.000 0.000
8 7.097130 0.000 -0.000 -0.000 -0.354 -0.000 0.000 0.000 -0.000 0.000
-0.661 0.000 -0.000 -0.661 -0.000 -0.000
9 7.887216 0.000 -0.005 -0.694 0.000 0.000 0.000 0.000 0.005 0.720
0.000 -0.000 -0.000 0.000 0.000 -0.000
10 7.887216 -0.000 0.694 -0.005 0.000 -0.000 0.000 -0.000 -0.720 0.005
-0.000 0.000 -0.000 0.000 -0.000 -0.000
11 12.030936 -0.930 -0.000 -0.000 0.000 0.000 0.000 -0.368 0.000 0.000
-0.000 0.000 -0.000 -0.000 -0.000 0.000
12 12.825515 -0.000 0.000 0.000 -0.000 -0.213 0.000 -0.000 -0.000 -0.000
-0.000 -0.777 -0.000 -0.000 0.593 0.000
13 12.825554 -0.000 0.000 0.000 0.000 -0.000 -0.213 -0.000 0.000 -0.000
0.000 -0.000 0.593 0.000 0.000 -0.777
14 30.407870 -0.368 0.000 0.000 -0.000 -0.000 -0.000 0.930 -0.000 -0.000
0.000 0.000 -0.000 0.000 -0.000 0.000
15 37.590543 -0.000 -0.000 0.000 0.935 -0.000 -0.000 0.000 0.000 0.000
-0.250 -0.000 -0.000 -0.250 -0.000 0.000

ecut = 40
1 0.243667 0.000 -0.000 0.000 -0.000 0.815 0.000 0.000 -0.000 0.000
-0.000 -0.471 -0.000 -0.000 -0.336 -0.000
2 0.243688 -0.000 0.000 -0.000 0.000 0.000 -0.815 -0.000 0.000 0.000
0.000 -0.000 0.336 0.000 -0.000 0.471
3 3.324505 -0.000 0.000 -0.000 -0.000 -0.000 -0.541 -0.000 0.000 -0.000
-0.000 -0.000 -0.731 -0.000 -0.000 -0.415
4 3.324513 0.000 -0.000 0.000 -0.000 0.541 -0.000 0.000 -0.000 0.000
-0.000 0.415 -0.000 -0.000 0.731 -0.000
5 5.963268 -0.000 0.002 0.702 -0.000 -0.000 -0.000 0.000 0.002 0.712
-0.000 -0.000 -0.000 -0.000 -0.000 -0.000
6 5.963268 -0.000 -0.702 0.002 0.000 -0.000 -0.000 -0.000 -0.712 0.002
0.000 -0.000 -0.000 0.000 -0.000 0.000
7 6.541910 -0.000 -0.000 -0.000 0.000 -0.000 -0.000 -0.000 -0.000 -0.000
-0.707 0.000 -0.000 0.707 -0.000 -0.000
8 6.852537 -0.000 -0.000 -0.000 -0.362 -0.000 0.000 0.000 -0.000 0.000
-0.659 0.000 0.000 -0.659 -0.000 0.000
9 7.788136 0.000 0.001 -0.712 0.000 -0.000 0.000 0.000 -0.001 0.702
0.000 0.000 -0.000 0.000 0.000 -0.000
10 7.788136 -0.000 -0.712 -0.001 -0.000 0.000 0.000 -0.000 0.702 0.001
-0.000 -0.000 0.000 -0.000 0.000 -0.000
11 11.970059 0.929 -0.000 0.000 -0.000 -0.000 -0.000 0.370 0.000 -0.000
-0.000 0.000 -0.000 -0.000 -0.000 0.000
12 12.684698 0.000 -0.000 0.000 -0.000 0.000 -0.205 0.000 -0.000 -0.000
-0.000 0.000 0.594 -0.000 -0.000 -0.778
13 12.684712 -0.000 -0.000 -0.000 0.000 0.205 0.000 -0.000 0.000 0.000
0.000 0.778 -0.000 0.000 -0.594 0.000
14 30.239420 -0.370 -0.000 -0.000 -0.000 -0.000 -0.000 0.929 0.000 -0.000
0.000 -0.000 0.000 0.000 0.000 -0.000
15 36.693310 0.000 -0.000 0.000 -0.932 0.000 -0.000 -0.000 -0.000 -0.000
0.256 0.000 0.000 0.256 -0.000 -0.000

ecut = 45 [Hartree]
1 0.881766 -0.000 0.000 -0.000 -0.000 -0.780 0.000 -0.000 0.000 -0.000
-0.000 0.489 -0.000 -0.000 0.390 -0.000
2 0.881826 0.000 0.000 0.000 -0.000 0.000 0.780 0.000 0.000 -0.000
0.000 -0.000 -0.390 -0.000 -0.000 -0.489
3 3.448139 -0.000 -0.000 -0.000 -0.000 -0.000 -0.594 -0.000 -0.000 -0.000
0.000 -0.000 -0.706 -0.000 -0.000 -0.385
4 3.448168 0.000 -0.000 -0.000 -0.000 0.594 -0.000 0.000 -0.000 -0.000
0.000 0.385 -0.000 -0.000 0.706 -0.000
5 5.985984 0.000 0.004 -0.706 0.000 0.000 0.000 0.000 0.004 -0.709
0.000 -0.000 0.000 0.000 -0.000 0.000
6 5.985984 -0.000 -0.706 -0.004 0.000 -0.000 0.000 -0.000 -0.709 -0.004
0.000 -0.000 0.000 0.000 -0.000 -0.000
7 6.583952 -0.000 -0.000 -0.000 -0.000 -0.000 -0.000 -0.000 0.000 -0.000
-0.707 -0.000 -0.000 0.707 0.000 -0.000
8 6.951193 -0.000 -0.000 -0.000 -0.358 0.000 0.000 -0.000 -0.000 0.000
-0.660 0.000 0.000 -0.660 -0.000 0.000
9 7.820245 -0.000 0.002 -0.709 0.000 -0.000 0.000 -0.000 -0.002 0.706
0.000 0.000 -0.000 0.000 -0.000 0.000
10 7.820245 0.000 0.709 0.002 0.000 -0.000 0.000 0.000 -0.706 -0.002
0.000 0.000 0.000 0.000 0.000 -0.000
11 12.010457 -0.930 0.000 0.000 0.000 0.000 0.000 -0.367 -0.000 -0.000
0.000 -0.000 0.000 0.000 0.000 -0.000
12 12.765853 -0.000 -0.000 0.000 0.000 -0.000 0.195 -0.000 -0.000 -0.000
-0.000 -0.000 -0.591 0.000 0.000 0.783
13 12.765884 -0.000 -0.000 0.000 0.000 0.195 0.000 -0.000 -0.000 -0.000
-0.000 0.783 -0.000 0.000 -0.591 0.000
14 30.411583 -0.367 -0.000 0.000 0.000 -0.000 -0.000 0.930 0.000 0.000
-0.000 -0.000 -0.000 -0.000 -0.000 -0.000
15 37.133106 -0.000 -0.000 0.000 -0.934 -0.000 -0.000 0.000 -0.000 -0.000
0.253 0.000 -0.000 0.253 -0.000 -0.000

ecut = 50 [Hartree]
1 0.650240 0.000 0.000 -0.000 -0.000 -0.789 0.000 0.000 0.000 0.000
-0.000 0.489 -0.000 -0.000 0.372 -0.000
2 0.650363 0.000 -0.000 -0.000 -0.000 -0.000 -0.789 -0.000 -0.000 0.000
-0.000 0.000 0.372 0.000 0.000 0.489
3 3.437387 -0.000 0.000 0.000 -0.000 0.000 -0.580 -0.000 -0.000 -0.000
0.000 0.000 -0.714 -0.000 0.000 -0.393
4 3.437478 0.000 -0.000 0.000 0.000 0.580 0.000 0.000 -0.000 -0.000
0.000 0.393 0.000 0.000 0.714 0.000
5 5.995243 -0.000 -0.711 -0.001 -0.000 -0.000 0.000 -0.000 -0.703 -0.001
-0.000 -0.000 -0.000 -0.000 -0.000 -0.000
6 5.995243 0.000 0.001 -0.711 0.000 -0.000 -0.000 0.000 0.001 -0.703
0.000 0.000 -0.000 0.000 -0.000 -0.000
7 6.604551 -0.000 0.000 0.000 0.000 -0.000 0.000 -0.000 -0.000 0.000
0.707 -0.000 0.000 -0.707 0.000 0.000
8 6.949866 -0.000 -0.000 0.000 0.359 -0.000 -0.000 -0.000 0.000 -0.000
0.660 0.000 -0.000 0.660 -0.000 -0.000
9 7.822770 0.000 0.703 0.000 0.000 -0.000 0.000 0.000 -0.711 -0.000
0.000 0.000 0.000 0.000 0.000 -0.000
10 7.822770 0.000 0.000 -0.703 0.000 0.000 0.000 0.000 -0.000 0.711
0.000 0.000 -0.000 0.000 0.000 -0.000
11 11.989834 0.930 -0.000 0.000 0.000 -0.000 -0.000 0.369 0.000 -0.000
0.000 -0.000 0.000 0.000 0.000 -0.000
12 12.762793 0.000 0.000 -0.000 0.000 -0.000 0.203 0.000 -0.000 0.000
-0.000 -0.000 -0.593 0.000 0.000 0.779
13 12.762825 -0.000 0.000 -0.000 0.000 -0.203 -0.000 -0.000 -0.000 0.000
0.000 -0.779 0.000 0.000 0.593 -0.000
14 30.326878 -0.369 -0.000 -0.000 -0.000 -0.000 -0.000 0.930 0.000 -0.000
0.000 -0.000 -0.000 0.000 -0.000 -0.000
15 37.050130 0.000 -0.000 0.000 -0.933 -0.000 0.000 -0.000 -0.000 -0.000
0.254 -0.000 -0.000 0.254 0.000 -0.000

ecut = 55 [Hartree]
1 0.622722 0.000 0.000 0.000 0.000 -0.000 0.793 0.000 0.000 0.000
0.000 0.000 -0.368 0.000 0.000 -0.486
2 0.622731 -0.000 -0.000 -0.000 -0.000 0.793 0.000 -0.000 -0.000 0.000
-0.000 -0.486 -0.000 -0.000 -0.368 -0.000
3 3.421091 0.000 0.000 0.000 0.000 -0.575 0.000 0.000 0.000 0.000
0.000 -0.396 0.000 0.000 -0.716 0.000
4 3.421116 0.000 0.000 -0.000 0.000 0.000 0.575 0.000 0.000 -0.000
0.000 0.000 0.716 -0.000 0.000 0.396
5 5.991277 -0.000 -0.007 -0.708 0.000 -0.000 -0.000 -0.000 -0.007 -0.706
0.000 -0.000 -0.000 0.000 -0.000 -0.000
6 5.991277 0.000 -0.708 0.007 0.000 -0.000 0.000 0.000 -0.706 0.007
0.000 -0.000 0.000 0.000 -0.000 0.000
7 6.596521 -0.000 0.000 -0.000 0.000 -0.000 0.000 -0.000 0.000 0.000
-0.707 -0.000 0.000 0.707 0.000 0.000
8 6.936373 0.000 0.000 0.000 0.360 0.000 -0.000 -0.000 0.000 -0.000
0.660 -0.000 -0.000 0.660 0.000 0.000
9 7.824611 -0.000 0.005 0.706 -0.000 0.000 0.000 -0.000 -0.005 -0.708
-0.000 -0.000 0.000 -0.000 0.000 0.000
10 7.824611 -0.000 0.706 -0.005 0.000 -0.000 0.000 -0.000 -0.708 0.005
0.000 0.000 0.000 -0.000 -0.000 -0.000
11 11.993944 0.929 0.000 0.000 -0.000 -0.000 -0.000 0.369 -0.000 -0.000
-0.000 -0.000 -0.000 -0.000 0.000 0.000
12 12.742406 -0.000 -0.000 -0.000 -0.000 0.000 0.202 -0.000 -0.000 -0.000
0.000 0.000 -0.593 -0.000 -0.000 0.779
13 12.742423 0.000 -0.000 0.000 0.000 0.202 -0.000 0.000 0.000 -0.000
0.000 0.779 0.000 0.000 -0.593 -0.000
14 30.318716 -0.369 0.000 0.000 -0.000 0.000 -0.000 0.929 -0.000 -0.000
0.000 0.000 -0.000 0.000 -0.000 -0.000
15 37.007008 0.000 -0.000 0.000 -0.933 0.000 0.000 -0.000 -0.000 -0.000
0.254 0.000 -0.000 0.254 -0.000 0.000

ecut = 60 [Hartree]
1 0.615824 0.000 -0.000 -0.000 0.000 0.000 -0.795 0.000 0.000 0.000
-0.000 -0.000 0.365 0.000 -0.000 0.484
2 0.615852 0.000 -0.000 0.000 -0.000 0.795 0.000 0.000 -0.000 0.000
-0.000 -0.484 -0.000 -0.000 -0.365 -0.000
3 3.405616 -0.000 0.000 0.000 0.000 -0.572 0.000 -0.000 0.000 -0.000
0.000 -0.399 0.000 -0.000 -0.717 0.000
4 3.405647 -0.000 0.000 -0.000 0.000 0.000 0.572 -0.000 0.000 -0.000
-0.000 0.000 0.717 0.000 0.000 0.399
5 5.988686 0.000 -0.019 0.707 -0.000 -0.000 0.000 0.000 -0.019 0.707
-0.000 -0.000 0.000 -0.000 -0.000 0.000
6 5.988686 0.000 -0.707 -0.019 0.000 -0.000 0.000 0.000 -0.707 -0.019
0.000 0.000 0.000 0.000 -0.000 -0.000
7 6.590920 0.000 0.000 0.000 -0.000 -0.000 0.000 0.000 -0.000 0.000
0.707 -0.000 0.000 -0.707 0.000 0.000
8 6.927831 -0.000 -0.000 0.000 0.360 0.000 -0.000 -0.000 0.000 0.000
0.660 -0.000 -0.000 0.660 -0.000 -0.000
9 7.820878 -0.000 0.018 -0.707 0.000 -0.000 0.000 -0.000 -0.018 0.707
0.000 0.000 -0.000 0.000 -0.000 -0.000
10 7.820878 0.000 -0.707 -0.018 -0.000 0.000 0.000 0.000 0.707 0.018
-0.000 -0.000 -0.000 -0.000 0.000 -0.000
11 11.993530 0.929 0.000 -0.000 -0.000 -0.000 -0.000 0.369 -0.000 -0.000
0.000 -0.000 0.000 0.000 0.000 -0.000
12 12.739142 0.000 -0.000 -0.000 0.000 0.000 0.201 0.000 0.000 0.000
-0.000 0.000 -0.593 0.000 -0.000 0.779
13 12.739144 0.000 -0.000 -0.000 -0.000 0.201 -0.000 0.000 0.000 0.000
0.000 0.779 0.000 -0.000 -0.593 -0.000
14 30.321260 -0.369 0.000 0.000 -0.000 0.000 -0.000 0.929 -0.000 -0.000
0.000 0.000 0.000 0.000 -0.000 -0.000
15 36.981755 0.000 -0.000 0.000 -0.933 0.000 0.000 -0.000 -0.000 -0.000
0.255 -0.000 -0.000 0.255 0.000 0.000

ecut = 65 [Hartree]
1 0.602196 0.000 0.000 0.000 -0.000 -0.798 -0.000 0.000 0.000 0.000
-0.000 0.482 0.000 -0.000 0.363 0.000
2 0.602245 0.000 -0.000 0.000 -0.000 -0.000 0.798 0.000 0.000 -0.000
0.000 0.000 -0.363 -0.000 0.000 -0.482
3 3.395252 -0.000 0.000 -0.000 0.000 0.000 0.569 -0.000 0.000 -0.000
-0.000 0.000 0.718 0.000 0.000 0.401
4 3.395265 -0.000 -0.000 -0.000 0.000 0.569 -0.000 -0.000 -0.000 -0.000
0.000 0.401 -0.000 0.000 0.718 -0.000
5 5.987714 -0.000 0.014 0.707 -0.000 0.000 0.000 -0.000 0.014 0.707
0.000 0.000 0.000 -0.000 0.000 0.000
6 5.987714 0.000 -0.707 0.014 0.000 -0.000 0.000 0.000 -0.707 0.014
0.000 -0.000 0.000 0.000 -0.000 0.000
7 6.588096 0.000 -0.000 -0.000 0.000 -0.000 0.000 -0.000 -0.000 -0.000
0.707 0.000 0.000 -0.707 0.000 0.000
8 6.922382 0.000 0.000 0.000 0.360 -0.000 -0.000 0.000 0.000 -0.000
0.660 0.000 -0.000 0.660 -0.000 -0.000
9 7.818962 0.000 0.012 0.707 -0.000 -0.000 -0.000 -0.000 -0.012 -0.707
-0.000 0.000 0.000 -0.000 0.000 0.000
10 7.818962 0.000 -0.707 0.012 -0.000 0.000 -0.000 0.000 0.707 -0.012
-0.000 -0.000 -0.000 -0.000 -0.000 0.000
11 11.993228 0.929 0.000 -0.000 -0.000 0.000 -0.000 0.369 -0.000 0.000
-0.000 -0.000 0.000 -0.000 0.000 -0.000
12 12.736876 0.000 0.000 -0.000 0.000 -0.000 0.201 0.000 -0.000 0.000
0.000 -0.000 -0.593 0.000 0.000 0.779
13 12.736912 0.000 -0.000 0.000 -0.000 0.201 0.000 0.000 0.000 -0.000
-0.000 0.779 -0.000 -0.000 -0.593 0.000
14 30.322459 -0.369 0.000 0.000 -0.000 0.000 -0.000 0.929 -0.000 0.000
0.000 0.000 -0.000 0.000 -0.000 -0.000
15 36.961780 0.000 -0.000 0.000 -0.933 -0.000 0.000 -0.000 0.000 -0.000
0.255 -0.000 -0.000 0.255 0.000 -0.000

ecut = 70 [Hartree]
1 0.570639 -0.000 -0.000 -0.000 -0.000 -0.000 -0.798 -0.000 -0.000 0.000
0.000 0.000 0.361 -0.000 0.000 0.482
2 0.570645 -0.000 0.000 0.000 0.000 -0.798 0.000 -0.000 0.000 -0.000
0.000 0.482 -0.000 0.000 0.361 -0.000
3 3.394085 -0.000 0.000 0.000 -0.000 0.000 0.568 -0.000 0.000 0.000
0.000 0.000 0.719 -0.000 0.000 0.401
4 3.394087 -0.000 0.000 0.000 0.000 -0.568 0.000 -0.000 0.000 -0.000
0.000 -0.401 0.000 -0.000 -0.719 0.000
5 5.986515 -0.000 -0.707 -0.005 0.000 -0.000 0.000 -0.000 -0.708 -0.005
0.000 -0.000 0.000 0.000 -0.000 -0.000
6 5.986515 -0.000 -0.005 0.707 0.000 -0.000 -0.000 -0.000 -0.005 0.708
0.000 -0.000 -0.000 -0.000 -0.000 -0.000
7 6.587487 0.000 -0.000 -0.000 -0.000 -0.000 -0.000 0.000 -0.000 -0.000
0.707 -0.000 -0.000 -0.707 0.000 -0.000
8 6.918812 -0.000 0.000 -0.000 -0.360 -0.000 -0.000 0.000 -0.000 0.000
-0.660 -0.000 -0.000 -0.660 0.000 -0.000
9 7.817227 -0.000 -0.708 -0.002 -0.000 0.000 -0.000 -0.000 0.707 0.002
-0.000 -0.000 -0.000 -0.000 0.000 0.000
10 7.817227 0.000 0.002 -0.708 0.000 -0.000 0.000 0.000 -0.002 0.707
0.000 0.000 -0.000 0.000 0.000 -0.000
11 11.989303 0.929 -0.000 0.000 -0.000 -0.000 0.000 0.369 0.000 -0.000
0.000 0.000 -0.000 0.000 -0.000 0.000
12 12.725944 -0.000 -0.000 0.000 -0.000 0.202 -0.000 -0.000 0.000 0.000
0.000 0.779 0.000 -0.000 -0.594 -0.000
13 12.725952 0.000 -0.000 -0.000 0.000 -0.000 -0.202 0.000 0.000 -0.000
-0.000 -0.000 0.594 0.000 0.000 -0.779
14 30.302732 -0.369 -0.000 -0.000 -0.000 -0.000 0.000 0.929 0.000 -0.000
0.000 -0.000 -0.000 0.000 0.000 0.000
15 36.948464 0.000 -0.000 -0.000 -0.933 0.000 -0.000 -0.000 -0.000 -0.000
0.255 -0.000 0.000 0.255 0.000 -0.000

ecut = 75 and 80 for the X point are not calculated becasue of memory shortage.

[abinit-forum] ecut for the RF calculation of BaTiO3, Takeshi NISHIMATSU, 03/04/2009
- Re: [abinit-forum] ecut for the RF calculation of BaTiO3, Xavier Gonze, 03/04/2009
  - Re: [abinit-forum] ecut for the RF calculation of BaTiO3, Takeshi NISHIMATSU, 03/05/2009
    - Re: [abinit-forum] ecut for the RF calculation of BaTiO3, Eric Bousquet, 03/05/2009
      - Re: [abinit-forum] ecut for the RF calculation of BaTiO3, Takeshi NISHIMATSU, 03/07/2009
      - Re: [abinit-forum] ecut for the RF calculation of BaTiO3, Takeshi NISHIMATSU, 03/10/2009
  - [abinit-forum] nband value and occopt, Souraya Goumri-Said, 03/05/2009