← Back to branch summary

~albertog/siesta/efield-2.5-jjunquera

« back to all changes in this revision

Viewing changes to Tests/Reference/force_2.out

  • Committer: Alberto Garcia
  • Date: 2007-11-21 11:45:49 UTC
  • mfrom: (192.1.67)
  • Revision ID: Arch-1:siesta@uam.es--2006%siesta-devel--reference--2.3--patch-1
Direct merge into master branch of initial BSC changes
The changes along the BSC branches, up to the end of the restructuring
of siesta.F and associated changes, have been merged into a direct
descendant of the main development line. The BSC work originally
started as a branch of siesta-devel--reference--2.1--patch-29. Development
along 2.1 continued, and a new continuation branch 2.3 has been created
specifically for this merge.

Main patches applied:

 * ref@bsc--2007/siesta-bsc--master--2.1--base-0
   tag of siesta@uam.es--2006/siesta-bsc--reference--2.1--patch-7

 * ref@bsc--2007/siesta-bsc--master--2.1--patch-1
   Use of xalast in analysis routines. Exit of geometry loop

 * ref@bsc--2007/siesta-bsc--master--2.1--patch-2
   Creation of a module to hold the siesta options

 * ref@bsc--2007/siesta-bsc--master--2.1--patch-3
   New geometry module

 * ref@bsc--2007/siesta-bsc--master--2.1--patch-4
   New stub module for sparse matrices

 * ref@bsc--2007/siesta-bsc--master--2.1--patch-5
   More options for running tests

 * ref@bsc--2007/siesta-bsc--master--2.1--patch-6
   Encapsulation of k-point handling

 * ref@bsc--2007/siesta-bsc--master--2.1--patch-7
   Initialize iza in struct_init

 * ref@bsc--2007/siesta-bsc--master--2.1--patch-8
   Siesta_todo slimming by M. Quero

 * ref@bsc--2007/siesta-bsc--master--2.1--patch-9
   Fix import of no_l in born_charge

 * ref@bsc--2007/siesta-bsc--master--2.1--patch-10
   Some conversions to Fortran90 by M. Quero

 * ref@bsc--2007/siesta-bsc--master--2.1--patch-11
   Fourth session at the BSC

 * ref@bsc--2007/siesta-bsc--master--2.1--patch-12
   Clarification of the scope of the stress variables

 * ref@bsc--2007/siesta-bsc--master--2.1--patch-13
   Creation of siesta_forces

 * ref@bsc--2007/siesta-bsc--master--2.1--patch-14
   Replacement of some allocatables by pointers and automatics

 * ref@bsc--2007/siesta-bsc--master--2.1--patch-15
   New m_energies and m_steps modules. Back to old k-point behavior

 * ref@bsc--2007/siesta-bsc--master--2.1--patch-16
   Merge of removal of integer and real variables from siesta_todo

 * ref@bsc--2007/siesta-bsc--master--2.1--patch-17
   Fixes for troublesome bugs in reference code

 * ref@bsc--2007/siesta-bsc--master--2.1--patch-18
   Final cleanup of siesta_todo

 * ref@bsc--2007/siesta-bsc--master--2.1--patch-19
   Add character(len=*) routines to alloc.F90

 * ref@bsc--2007/siesta-bsc--master--2.1--patch-20
   Pointers in fixed and setspatial. si2x1h test added to bsc-Makefile

 * ref@bsc--2007/siesta-bsc--master--2.1--patch-21
   Re-organization of pulay module

 * ref@bsc--2007/siesta-bsc--master--2.1--patch-22
   Reorganization of hsparse/xijorb calls with new neighbor module

 * ref@bsc--2007/siesta-bsc--master--2.1--patch-23
   VPATH-aware compilation for multiple executable versions

 * ref@bsc--2007/siesta-bsc--master--2.1--patch-24
   Allocatables to pointers  I

 * ref@bsc--2007/siesta-bsc--master--2.1--patch-25
   Allocatables to pointers  II

 * ref@bsc--2007/siesta-bsc--master--2.1--patch-26
   Allocatables to pointers III

 * ref@bsc--2007/siesta-bsc--master--2.1--patch-27
   Allocatables to pointers IV -- new neighbor code + nspecies fix

 * ref@bsc--2007/siesta-bsc--master--2.1--patch-28
   Explicit array-ness in calls in initatom and cellxc

 * ref@bsc--2007/siesta-bsc--master--2.1--patch-29
   Fix wrong allocations in cellxc.F

 * ref@bsc--2007/siesta-bsc--master--2.1--patch-30
   Explicit array extents in initatom.f

 * ref@bsc--2007/siesta-bsc--master--2.1--patch-31
   Avoid shrinking of density-matrix arrays for extrapol.

 * ref@bsc--2007/siesta-bsc--master--2.1--patch-32
   Fix typo in state_init.F

 * ref@bsc--2007/siesta-bsc--master--2.1--patch-33
   Execute SCF loop when nscf=1

 * ref@bsc--2007/siesta-bsc--master--2.1--patch-34
   Clarify bounds of SCF loop in siesta_forces.

 * ref@bsc--2007/siesta-bsc--master--2.2--base-0
   tag of ref@bsc--2007/siesta-bsc--master--2.1--patch-31

 * ref@bsc--2007/siesta-bsc--master--2.3--base-0
   tag of ref@bsc--2007/siesta-bsc--master--2.2--base-0

 * ref@bsc--2007/siesta-bsc--master--2.3--patch-1
   Prepare CHANGES file

 * ref@bsc--2007/siesta-bsc--master--2.3--patch-2
   Merge patch-log for patch-30 from 2.1

 * ref@bsc--2007/siesta-bsc--master--2.3--patch-3
   New treatment of fractional atoms in VCA. Bug fix in lmxo

 * ref@bsc--2007/siesta-bsc--master--2.3--patch-4
   Re-enabling of kgrid update in variable-cell calculations

 * ref@bsc--2007/siesta-bsc--master--2.3--patch-5
   Fix typo in state_init.F

 * ref@bsc--2007/siesta-bsc--master--2.3--patch-6
   Merge filtering package by Eduardo Anglada

 * ref@bsc--2007/siesta-bsc--master--2.3--patch-7
   Add graphite_c6_full test for more realistic vdW test

 * ref@bsc--2007/siesta-bsc--master--2.3--patch-8
   Add patchlog for k-point fix

 * ref@bsc--2007/siesta-bsc--master--2.3--patch-9
   Fix of zmatrix code to deal with degenerate case

 * ref@bsc--2007/siesta-bsc--master--2.3--patch-10
   Units conversion in Util/Optical/optical.f. Scripts. Cosmetics

 * ref@bsc--2007/siesta-bsc--master--2.3--patch-11
   Implementation of basis_enthalpy calculation. Zmatrix dependency

 * ref@bsc--2007/siesta-bsc--master--2.3--patch-12
   Option to use fractional rc's for multiple zeta

 * ref@bsc--2007/siesta-bsc--master--2.3--patch-13
   Sync vpath changes. Update siesta.tex

 * ref@bsc--2007/siesta-bsc--master--2.3--patch-14
   Fix tag in MPI send/receive in mulliken

 * ref@bsc--2007/siesta-bsc--master--2.3--patch-15
   Option to use fractional rc's for multiple zeta

 * ref@bsc--2007/siesta-bsc--master--2.3--patch-16
   Merge of XML tester by Eduardo Anglada. Portability fixes

 * ref@bsc--2007/siesta-bsc--master--2.3--patch-17
   Fix marenostrum-mpi.make. Syntax in compare_m.f90

 * ref@bsc--2007/siesta-bsc--master--2.3--patch-18
   MareNostrum fixes. zdrot to blas. obj_setup. compare_m syntax

 * ref@bsc--2007/siesta-bsc--master--2.3--patch-19
   Sync to bsc--2.1. Change banner in CHANGES to bsc--2.3

 * ref@bsc--2007/siesta-bsc--master--2.3--patch-20
   Execute SCF loop when nscf=1

 * ref@bsc--2007/siesta-bsc--master--2.3--patch-21
   Refinements of XML tester. New ioncat program

 * ref@bsc--2007/siesta-bsc--master--2.3--patch-22
   Undef var in kgridinit, iohs MPI write, pdosg array bound, Origin shift

 * ref@bsc--2007/siesta-bsc--master--2.3--patch-23
   Zmatrix optimization enhancements. Sign change in MM stress.

 * ref@bsc--2007/siesta-bsc--master--2.3--patch-24
   Clarify bounds of SCF loop in siesta_forces.

 * ref@bsc--2007/siesta-bsc--master--2.3--patch-25
   Avoid MP-grid permutations in trivial gamma case

 * ref@bsc--2007/siesta-bsc--master--2.3--patch-26
   Sync to reference--2.1

 * siesta@uam.es--2006/siesta-bsc--reference--2.1--base-0
   tag of siesta@uam.es--2006/siesta-devel--reference--2.1--patch-29

 * siesta@uam.es--2006/siesta-bsc--reference--2.1--patch-1
   First stage of siesta.F splitting at BSC

 * siesta@uam.es--2006/siesta-bsc--reference--2.1--patch-2
   Created struct_init for initial geometry setup. New test force_2

 * siesta@uam.es--2006/siesta-bsc--reference--2.1--patch-3
   Consolidate geometry updates at the end of loop

 * siesta@uam.es--2006/siesta-bsc--reference--2.1--patch-4
   Initialize vol2 correctly in m_check_supercell.f

 * siesta@uam.es--2006/siesta-bsc--reference--2.1--patch-5
   Kgrid setup streamlined. Bands. Proximity check. Hsparse allocation 

 * siesta@uam.es--2006/siesta-bsc--reference--2.1--patch-6
   Work by Manuel Quero before the meeting on Feb 7th

 * siesta@uam.es--2006/siesta-bsc--reference--2.1--patch-7
   Moved the Born-effective-charge code




Show diffs side-by-side

added added

removed removed

Lines of Context:
Tests/Reference/force_2.out
 
1
Siesta Version: siesta@uam.es--2006/siesta-bsc--reference--2.1--patch-1
 
2
Architecture  : i686-pc-linux-gnu--unknown
 
3
Compiler flags: g95 -g -O2
 
4
SERIAL version
 
5
 
 
6
* Running in serial mode
 
7
>> Start of run:  29-JAN-2007  17:58:13
 
8
 
 
9
                           ***********************       
 
10
                           *  WELCOME TO SIESTA  *       
 
11
                           ***********************       
 
12
 
 
13
reinit: Reading from standard input
 
14
************************** Dump of input data file ****************************
 
15
Systemlabel  force_2
 
16
Systemname   si2 Force constants
 
17
NumberOfSpecies 1
 
18
%block chemicalspecieslabel
 
19
1 14 Si
 
20
%endblock chemicalspecieslabel
 
21
PAO.BasisSize SZ
 
22
MeshCutoff 130.0 Ry
 
23
# In a real calculation, a high degree of self-consistency should be achieved
 
24
#DM.MixingWeight  0.3
 
25
#DM.NumberPulay    3
 
26
#DM.Tolerance     1.0d-5
 
27
# Here we use Harris since we are just interested in testing the displacement
 
28
# logic
 
29
HarrisFunctional T
 
30
%include FC.fdf
 
31
************************** End of input data file *****************************
 
32
 
 
33
reinit: -----------------------------------------------------------------------
 
34
reinit: System Name: si2 Force constants                                        
 
35
reinit: -----------------------------------------------------------------------
 
36
reinit: System Label: force_2                                                     
 
37
reinit: -----------------------------------------------------------------------
 
38
 
 
39
initatom: Reading input for the pseudopotentials and atomic orbitals ----------
 
40
 Species number:  1  Label: Si Atomic number: 14
 
41
Ground state valence configuration:   3s02  3p02
 
42
Reading pseudopotential information in formatted form from Si.psf
 
43
 
 
44
Valence configuration for pseudopotential generation:
 
45
3s( 2.00) rc: 1.89
 
46
3p( 2.00) rc: 1.89
 
47
3d( 0.00) rc: 1.89
 
48
4f( 0.00) rc: 1.89
 
49
Dumping pseudopotential information in formatted form in Si.psdump
 
50
For Si, standard SIESTA heuristics set lmxkb to 2
 
51
 (one more than the basis l, including polarization orbitals).
 
52
Use PS.lmax or PS.KBprojectors blocks to override.
 
53
 
 
54
<basis_specs>
 
55
===============================================================================
 
56
Si                   Z=  14    Mass=  28.090        Charge= 0.17977+309
 
57
Lmxo=1 Lmxkb=2     BasisType=split      Semic=F
 
58
L=0  Nsemic=0  Cnfigmx=3
 
59
          n=1  nzeta=1  polorb=0
 
60
            splnorm:    .15000    
 
61
               vcte:    .00000    
 
62
               rinn:    .00000    
 
63
                rcs:    .00000    
 
64
            lambdas:    1.0000    
 
65
L=1  Nsemic=0  Cnfigmx=3
 
66
          n=1  nzeta=1  polorb=0
 
67
            splnorm:    .15000    
 
68
               vcte:    .00000    
 
69
               rinn:    .00000    
 
70
                rcs:    .00000    
 
71
            lambdas:    1.0000    
 
72
-------------------------------------------------------------------------------
 
73
L=0  Nkbl=1  erefs: 0.17977+309
 
74
L=1  Nkbl=1  erefs: 0.17977+309
 
75
L=2  Nkbl=1  erefs: 0.17977+309
 
76
===============================================================================
 
77
</basis_specs>
 
78
 
 
79
atom: Called for Si                    (Z =  14)
 
80
 
 
81
read_vps: Pseudopotential generation method:
 
82
read_vps: ATM3      Troullier-Martins                       
 
83
Total valence charge:    4.00000
 
84
 
 
85
xc_check: Exchange-correlation functional:
 
86
xc_check: Ceperley-Alder
 
87
V l=0 = -2*Zval/r beyond r=  2.5494
 
88
V l=1 = -2*Zval/r beyond r=  2.5494
 
89
V l=2 = -2*Zval/r beyond r=  2.5494
 
90
All V_l potentials equal beyond r=  1.8652
 
91
This should be close to max(r_c) in ps generation
 
92
All pots = -2*Zval/r beyond r=  2.5494
 
93
Using large-core scheme for Vlocal
 
94
 
 
95
atom: Estimated core radius    2.54944
 
96
 
 
97
atom: Including non-local core corrections could be a good idea
 
98
atom: Maximum radius for 4*pi*r*r*local-pseudopot. charge    2.85303
 
99
atom: Maximum radius for r*vlocal+2*Zval:    2.58151
 
100
GHOST: No ghost state for L =  0
 
101
GHOST: No ghost state for L =  1
 
102
GHOST: No ghost state for L =  2
 
103
 
 
104
KBgen: Kleinman-Bylander projectors: 
 
105
   l= 0   rc=  1.936440   el= -0.796617   Ekb=  4.661340   kbcos=  0.299756
 
106
   l= 1   rc=  1.936440   el= -0.307040   Ekb=  1.494238   kbcos=  0.301471
 
107
   l= 2   rc=  1.936440   el=  0.002313   Ekb= -2.808672   kbcos= -0.054903
 
108
 
 
109
KBgen: Total number of  Kleinman-Bylander projectors:    9
 
110
atom: -------------------------------------------------------------------------
 
111
 
 
112
atom: SANKEY-TYPE ORBITALS:
 
113
 
 
114
SPLIT: Orbitals with angular momentum L= 0
 
115
 
 
116
SPLIT: Basis orbitals for state 3s
 
117
 
 
118
SPLIT: PAO cut-off radius determined from an
 
119
SPLIT: energy shift=  0.020000 Ry
 
120
 
 
121
   izeta = 1
 
122
                 lambda =    1.000000
 
123
                     rc =    5.007352
 
124
                 energy =   -0.777669
 
125
                kinetic =    0.573829
 
126
    potential(screened) =   -1.351499
 
127
       potential(ionic) =   -3.827441
 
128
 
 
129
SPLIT: Orbitals with angular momentum L= 1
 
130
 
 
131
SPLIT: Basis orbitals for state 3p
 
132
 
 
133
SPLIT: PAO cut-off radius determined from an
 
134
SPLIT: energy shift=  0.020000 Ry
 
135
 
 
136
   izeta = 1
 
137
                 lambda =    1.000000
 
138
                     rc =    6.270866
 
139
                 energy =   -0.288955
 
140
                kinetic =    0.877930
 
141
    potential(screened) =   -1.166885
 
142
       potential(ionic) =   -3.426998
 
143
atom: Total number of Sankey-type orbitals:  4
 
144
 
 
145
atm_pop: Valence configuration (for local Pseudopot. screening):
 
146
 3s( 2.00)                                                            
 
147
 3p( 2.00)                                                            
 
148
Vna: chval, zval:    4.00000   4.00000
 
149
 
 
150
Vna:  Cut-off radius for the neutral-atom potential:   6.270866
 
151
 
 
152
atom: _________________________________________________________________________
 
153
 
 
154
prinput: Basis input ----------------------------------------------------------
 
155
 
 
156
PAO.BasisType split     
 
157
 
 
158
%block ChemicalSpeciesLabel
 
159
    1   14 Si                      # Species index, atomic number, species label
 
160
%endblock ChemicalSpeciesLabel
 
161
 
 
162
%block PAO.Basis                 # Define Basis set
 
163
Si                    2                    # Species label, number of l-shells
 
164
 n=3   0   1                         # n, l, Nzeta 
 
165
   5.007   
 
166
   1.000   
 
167
 n=3   1   1                         # n, l, Nzeta 
 
168
   6.271   
 
169
   1.000   
 
170
%endblock PAO.Basis
 
171
 
 
172
prinput: ----------------------------------------------------------------------
 
173
 
 
174
coor:   Atomic-coordinates input format  =     Fractional
 
175
 
 
176
siesta: Atomic coordinates (Bohr) and species
 
177
siesta:      0.00000   0.00000   0.00000  1        1
 
178
siesta:      2.56530   2.56530   2.56530  1        2
 
179
 
 
180
siesta: System type = bulk      
 
181
 
 
182
siesta: ******************** Simulation parameters ****************************
 
183
siesta:
 
184
siesta: The following are some of the parameters of the simulation.
 
185
siesta: A complete list of the parameters used, including default values,
 
186
siesta: can be found in file out.fdf
 
187
siesta:
 
188
redata: Number of spin components        =     1
 
189
redata: Long output                      =     F
 
190
redata: Number of Atomic Species         =     1
 
191
redata: Charge density info will appear in .RHO file
 
192
redata: Write Mulliken Pop.              =     NO
 
193
redata: Mesh Cutoff                      =   130.0000  Ry
 
194
redata: Net charge of the system         =     0.0000 |e|
 
195
redata: Max. number of SCF Iter          =    50
 
196
redata: Mixing is linear
 
197
redata: Mix DM in first SCF step ?       =     F
 
198
redata: Write Pulay info on disk?        =     F
 
199
redata: New DM Mixing Weight             =     0.2500
 
200
redata: New DM Occupancy tolerance       = 0.000000000001
 
201
redata: No kicks to SCF
 
202
redata: DM Mixing Weight for Kicks       =     0.5000
 
203
redata: DM Tolerance for SCF             =     0.000100
 
204
redata: Require Energy convergence for SCF =     F
 
205
redata: DM Energy tolerance for SCF      =     0.000100 eV
 
206
redata: Using Saved Data (generic)   =     F
 
207
redata: Use continuation files for DM    =     F
 
208
redata: Neglect nonoverlap interactions  =     F
 
209
redata: Method of Calculation            =     Diagonalization
 
210
redata: Divide and Conquer               =     T
 
211
redata: Electronic Temperature           =     0.0019  Ry
 
212
redata: Fix the spin of the system       =     F
 
213
redata: Dynamics option                  =     Force Constants Matrix Calculation
 
214
redata: Atomic displ for force constants  =     0.0400  Bohr
 
215
redata: First atom to move               =        1
 
216
redata: Last atom to move                =        2
 
217
redata: ***********************************************************************
 
218
 
 
219
initatomlists: Number of atoms, orbitals, and projectors:      2     8    18
 
220
Total number of electrons:     8.000000
 
221
Total ionic charge:     8.000000
 
222
 
 
223
siesta: k-grid: Number of k-points =     1
 
224
siesta: k-grid: Cutoff (effective) =     1.920 Ang
 
225
siesta: k-grid: Supercell and displacements
 
226
siesta: k-grid:    0   1   0      0.000
 
227
siesta: k-grid:    0   0   1      0.000
 
228
siesta: k-grid:    1   0   0      0.000
 
229
 
 
230
* Maximum dynamic memory allocated =     1 MB
 
231
 
 
232
siesta:                 ==============================
 
233
                            Begin FC step =      0
 
234
                            Undisplaced coordinates
 
235
                        ==============================
 
236
 
 
237
outcell: Unit cell vectors (Ang):
 
238
        0.000000    2.715000    2.715000
 
239
        2.715000    0.000000    2.715000
 
240
        2.715000    2.715000    0.000000
 
241
 
 
242
outcell: Cell vector modules (Ang)   :    3.839590    3.839590    3.839590
 
243
outcell: Cell angles (23,13,12) (deg):     60.0000     60.0000     60.0000
 
244
outcell: Cell volume (Ang**3)        :     40.0258
 
245
 
 
246
InitMesh: MESH =    24 x    24 x    24 =       13824
 
247
InitMesh: Mesh cutoff (required, used) =   130.000   161.974 Ry
 
248
 
 
249
* Maximum dynamic memory allocated =     4 MB
 
250
 
 
251
stepf: Fermi-Dirac step function
 
252
 
 
253
siesta: Program's energy decomposition (eV):
 
254
siesta: Eions   =       380.802124
 
255
siesta: Ena     =       115.502327
 
256
siesta: Ekin    =        82.397263
 
257
siesta: Enl     =        39.713031
 
258
siesta: DEna    =        -8.221426
 
259
siesta: DUscf   =         0.353217
 
260
siesta: DUext   =         0.000000
 
261
siesta: Exc     =       -66.294639
 
262
siesta: eta*DQ  =         0.000000
 
263
siesta: Emadel  =         0.000000
 
264
siesta: Emeta   =         0.000000
 
265
siesta: Emolmec =         0.000000
 
266
siesta: Ekinion =         0.000000
 
267
siesta: Eharris =      -194.939312
 
268
siesta: Etot    =      -217.352349
 
269
siesta: FreeEng =      -217.352349
 
270
 
 
271
siesta: Eharris(eV) =    -194.939312
 
272
 
 
273
 
 
274
siesta: Eharris(eV) =    -194.939312
 
275
 
 
276
timer: Routine,Calls,Time,% = IterSCF        1       0.694  64.62
 
277
elaps: Routine,Calls,Wall,% = IterSCF        1       0.694  64.64
 
278
 
 
279
siesta: E_KS - E_eggbox =      -194.6964
 
280
 
 
281
siesta: Atomic forces (eV/Ang):
 
282
----------------------------------------
 
283
   Tot    0.000080    0.000082    0.000082
 
284
----------------------------------------
 
285
   Max    0.000082
 
286
   Res    0.000058    sqrt( Sum f_i^2 / 3N )
 
287
----------------------------------------
 
288
   Max    0.000082    constrained
 
289
 
 
290
Stress-tensor-Voigt (kbar):     -476.61     -476.61     -476.61       -0.00        0.00        0.00
 
291
Target enthalpy (eV/cell)     -194.6964
 
292
 
 
293
* Maximum dynamic memory allocated =     4 MB
 
294
 
 
295
siesta:                 ==============================
 
296
                            Begin FC step =      1
 
297
                            displace atom        1
 
298
                            in direction         1
 
299
                            by        -0.0400 Bohr
 
300
                        ==============================
 
301
 
 
302
outcell: Unit cell vectors (Ang):
 
303
        0.000000    2.715000    2.715000
 
304
        2.715000    0.000000    2.715000
 
305
        2.715000    2.715000    0.000000
 
306
 
 
307
outcell: Cell vector modules (Ang)   :    3.839590    3.839590    3.839590
 
308
outcell: Cell angles (23,13,12) (deg):     60.0000     60.0000     60.0000
 
309
outcell: Cell volume (Ang**3)        :     40.0258
 
310
 
 
311
siesta: Eharris(eV) =    -194.932993
 
312
 
 
313
 
 
314
siesta: Eharris(eV) =    -194.932993
 
315
 
 
316
 
 
317
siesta: Atomic forces (eV/Ang):
 
318
----------------------------------------
 
319
   Tot   -0.010059    0.000081    0.000081
 
320
----------------------------------------
 
321
   Max    0.782261
 
322
   Res    0.448744    sqrt( Sum f_i^2 / 3N )
 
323
----------------------------------------
 
324
   Max    0.782261    constrained
 
325
 
 
326
Stress-tensor-Voigt (kbar):     -479.67     -479.81     -479.81        0.00       22.98       -0.00
 
327
Target enthalpy (eV/cell)     -194.6894
 
328
 
 
329
* Maximum dynamic memory allocated =     4 MB
 
330
 
 
331
siesta:                 ==============================
 
332
                            Begin FC step =      2
 
333
                            displace atom        1
 
334
                            in direction         1
 
335
                            by         0.0400 Bohr
 
336
                        ==============================
 
337
 
 
338
outcell: Unit cell vectors (Ang):
 
339
        0.000000    2.715000    2.715000
 
340
        2.715000    0.000000    2.715000
 
341
        2.715000    2.715000    0.000000
 
342
 
 
343
outcell: Cell vector modules (Ang)   :    3.839590    3.839590    3.839590
 
344
outcell: Cell angles (23,13,12) (deg):     60.0000     60.0000     60.0000
 
345
outcell: Cell volume (Ang**3)        :     40.0258
 
346
 
 
347
siesta: Eharris(eV) =    -194.932993
 
348
 
 
349
 
 
350
siesta: Eharris(eV) =    -194.932993
 
351
 
 
352
 
 
353
siesta: Atomic forces (eV/Ang):
 
354
----------------------------------------
 
355
   Tot    0.010222    0.000081    0.000081
 
356
----------------------------------------
 
357
   Max    0.782423
 
358
   Res    0.448791    sqrt( Sum f_i^2 / 3N )
 
359
----------------------------------------
 
360
   Max    0.782423    constrained
 
361
 
 
362
Stress-tensor-Voigt (kbar):     -479.67     -479.81     -479.81       -0.00      -22.98       -0.00
 
363
Target enthalpy (eV/cell)     -194.6894
 
364
 
 
365
* Maximum dynamic memory allocated =     4 MB
 
366
 
 
367
siesta:                 ==============================
 
368
                            Begin FC step =      3
 
369
                            displace atom        1
 
370
                            in direction         2
 
371
                            by        -0.0400 Bohr
 
372
                        ==============================
 
373
 
 
374
outcell: Unit cell vectors (Ang):
 
375
        0.000000    2.715000    2.715000
 
376
        2.715000    0.000000    2.715000
 
377
        2.715000    2.715000    0.000000
 
378
 
 
379
outcell: Cell vector modules (Ang)   :    3.839590    3.839590    3.839590
 
380
outcell: Cell angles (23,13,12) (deg):     60.0000     60.0000     60.0000
 
381
outcell: Cell volume (Ang**3)        :     40.0258
 
382
 
 
383
siesta: Eharris(eV) =    -194.932993
 
384
 
 
385
 
 
386
siesta: Eharris(eV) =    -194.932993
 
387
 
 
388
 
 
389
siesta: Atomic forces (eV/Ang):
 
390
----------------------------------------
 
391
   Tot    0.000081   -0.010059    0.000082
 
392
----------------------------------------
 
393
   Max    0.782261
 
394
   Res    0.448744    sqrt( Sum f_i^2 / 3N )
 
395
----------------------------------------
 
396
   Max    0.782261    constrained
 
397
 
 
398
Stress-tensor-Voigt (kbar):     -479.81     -479.67     -479.81       -0.00        0.00       22.98
 
399
Target enthalpy (eV/cell)     -194.6894
 
400
 
 
401
* Maximum dynamic memory allocated =     4 MB
 
402
 
 
403
siesta:                 ==============================
 
404
                            Begin FC step =      4
 
405
                            displace atom        1
 
406
                            in direction         2
 
407
                            by         0.0400 Bohr
 
408
                        ==============================
 
409
 
 
410
outcell: Unit cell vectors (Ang):
 
411
        0.000000    2.715000    2.715000
 
412
        2.715000    0.000000    2.715000
 
413
        2.715000    2.715000    0.000000
 
414
 
 
415
outcell: Cell vector modules (Ang)   :    3.839590    3.839590    3.839590
 
416
outcell: Cell angles (23,13,12) (deg):     60.0000     60.0000     60.0000
 
417
outcell: Cell volume (Ang**3)        :     40.0258
 
418
 
 
419
siesta: Eharris(eV) =    -194.932993
 
420
 
 
421
 
 
422
siesta: Eharris(eV) =    -194.932993
 
423
 
 
424
 
 
425
siesta: Atomic forces (eV/Ang):
 
426
----------------------------------------
 
427
   Tot    0.000080    0.010223    0.000081
 
428
----------------------------------------
 
429
   Max    0.782424
 
430
   Res    0.448791    sqrt( Sum f_i^2 / 3N )
 
431
----------------------------------------
 
432
   Max    0.782424    constrained
 
433
 
 
434
Stress-tensor-Voigt (kbar):     -479.81     -479.67     -479.81       -0.00        0.00      -22.98
 
435
Target enthalpy (eV/cell)     -194.6894
 
436
 
 
437
* Maximum dynamic memory allocated =     4 MB
 
438
 
 
439
siesta:                 ==============================
 
440
                            Begin FC step =      5
 
441
                            displace atom        1
 
442
                            in direction         3
 
443
                            by        -0.0400 Bohr
 
444
                        ==============================
 
445
 
 
446
outcell: Unit cell vectors (Ang):
 
447
        0.000000    2.715000    2.715000
 
448
        2.715000    0.000000    2.715000
 
449
        2.715000    2.715000    0.000000
 
450
 
 
451
outcell: Cell vector modules (Ang)   :    3.839590    3.839590    3.839590
 
452
outcell: Cell angles (23,13,12) (deg):     60.0000     60.0000     60.0000
 
453
outcell: Cell volume (Ang**3)        :     40.0258
 
454
 
 
455
siesta: Eharris(eV) =    -194.932993
 
456
 
 
457
 
 
458
siesta: Eharris(eV) =    -194.932993
 
459
 
 
460
 
 
461
siesta: Atomic forces (eV/Ang):
 
462
----------------------------------------
 
463
   Tot    0.000081    0.000081   -0.010059
 
464
----------------------------------------
 
465
   Max    0.782260
 
466
   Res    0.448744    sqrt( Sum f_i^2 / 3N )
 
467
----------------------------------------
 
468
   Max    0.782260    constrained
 
469
 
 
470
Stress-tensor-Voigt (kbar):     -479.81     -479.81     -479.67       22.98       -0.00       -0.00
 
471
Target enthalpy (eV/cell)     -194.6894
 
472
 
 
473
* Maximum dynamic memory allocated =     4 MB
 
474
 
 
475
siesta:                 ==============================
 
476
                            Begin FC step =      6
 
477
                            displace atom        1
 
478
                            in direction         3
 
479
                            by         0.0400 Bohr
 
480
                        ==============================
 
481
 
 
482
outcell: Unit cell vectors (Ang):
 
483
        0.000000    2.715000    2.715000
 
484
        2.715000    0.000000    2.715000
 
485
        2.715000    2.715000    0.000000
 
486
 
 
487
outcell: Cell vector modules (Ang)   :    3.839590    3.839590    3.839590
 
488
outcell: Cell angles (23,13,12) (deg):     60.0000     60.0000     60.0000
 
489
outcell: Cell volume (Ang**3)        :     40.0258
 
490
 
 
491
siesta: Eharris(eV) =    -194.932993
 
492
 
 
493
 
 
494
siesta: Eharris(eV) =    -194.932993
 
495
 
 
496
 
 
497
siesta: Atomic forces (eV/Ang):
 
498
----------------------------------------
 
499
   Tot    0.000081    0.000081    0.010222
 
500
----------------------------------------
 
501
   Max    0.782424
 
502
   Res    0.448792    sqrt( Sum f_i^2 / 3N )
 
503
----------------------------------------
 
504
   Max    0.782424    constrained
 
505
 
 
506
Stress-tensor-Voigt (kbar):     -479.81     -479.81     -479.67      -22.98        0.00       -0.00
 
507
Target enthalpy (eV/cell)     -194.6894
 
508
 
 
509
* Maximum dynamic memory allocated =     4 MB
 
510
 
 
511
siesta:                 ==============================
 
512
                            Begin FC step =      7
 
513
                            displace atom        2
 
514
                            in direction         1
 
515
                            by        -0.0400 Bohr
 
516
                        ==============================
 
517
 
 
518
outcell: Unit cell vectors (Ang):
 
519
        0.000000    2.715000    2.715000
 
520
        2.715000    0.000000    2.715000
 
521
        2.715000    2.715000    0.000000
 
522
 
 
523
outcell: Cell vector modules (Ang)   :    3.839590    3.839590    3.839590
 
524
outcell: Cell angles (23,13,12) (deg):     60.0000     60.0000     60.0000
 
525
outcell: Cell volume (Ang**3)        :     40.0258
 
526
 
 
527
siesta: Eharris(eV) =    -194.932994
 
528
 
 
529
 
 
530
siesta: Eharris(eV) =    -194.932994
 
531
 
 
532
 
 
533
siesta: Atomic forces (eV/Ang):
 
534
----------------------------------------
 
535
   Tot   -0.010141   -0.000000    0.000000
 
536
----------------------------------------
 
537
   Max    0.782341
 
538
   Res    0.448767    sqrt( Sum f_i^2 / 3N )
 
539
----------------------------------------
 
540
   Max    0.782341    constrained
 
541
 
 
542
Stress-tensor-Voigt (kbar):     -479.67     -479.81     -479.81       -0.00      -22.98       -0.00
 
543
Target enthalpy (eV/cell)     -194.6894
 
544
 
 
545
* Maximum dynamic memory allocated =     4 MB
 
546
 
 
547
siesta:                 ==============================
 
548
                            Begin FC step =      8
 
549
                            displace atom        2
 
550
                            in direction         1
 
551
                            by         0.0400 Bohr
 
552
                        ==============================
 
553
 
 
554
outcell: Unit cell vectors (Ang):
 
555
        0.000000    2.715000    2.715000
 
556
        2.715000    0.000000    2.715000
 
557
        2.715000    2.715000    0.000000
 
558
 
 
559
outcell: Cell vector modules (Ang)   :    3.839590    3.839590    3.839590
 
560
outcell: Cell angles (23,13,12) (deg):     60.0000     60.0000     60.0000
 
561
outcell: Cell volume (Ang**3)        :     40.0258
 
562
 
 
563
siesta: Eharris(eV) =    -194.932993
 
564
 
 
565
 
 
566
siesta: Eharris(eV) =    -194.932993
 
567
 
 
568
 
 
569
siesta: Atomic forces (eV/Ang):
 
570
----------------------------------------
 
571
   Tot    0.010141    0.000001   -0.000000
 
572
----------------------------------------
 
573
   Max    0.782341
 
574
   Res    0.448767    sqrt( Sum f_i^2 / 3N )
 
575
----------------------------------------
 
576
   Max    0.782341    constrained
 
577
 
 
578
Stress-tensor-Voigt (kbar):     -479.67     -479.81     -479.81       -0.00       22.98       -0.00
 
579
Target enthalpy (eV/cell)     -194.6894
 
580
 
 
581
* Maximum dynamic memory allocated =     4 MB
 
582
 
 
583
siesta:                 ==============================
 
584
                            Begin FC step =      9
 
585
                            displace atom        2
 
586
                            in direction         2
 
587
                            by        -0.0400 Bohr
 
588
                        ==============================
 
589
 
 
590
outcell: Unit cell vectors (Ang):
 
591
        0.000000    2.715000    2.715000
 
592
        2.715000    0.000000    2.715000
 
593
        2.715000    2.715000    0.000000
 
594
 
 
595
outcell: Cell vector modules (Ang)   :    3.839590    3.839590    3.839590
 
596
outcell: Cell angles (23,13,12) (deg):     60.0000     60.0000     60.0000
 
597
outcell: Cell volume (Ang**3)        :     40.0258
 
598
 
 
599
siesta: Eharris(eV) =    -194.932993
 
600
 
 
601
 
 
602
siesta: Eharris(eV) =    -194.932993
 
603
 
 
604
 
 
605
siesta: Atomic forces (eV/Ang):
 
606
----------------------------------------
 
607
   Tot    0.000000   -0.010141   -0.000000
 
608
----------------------------------------
 
609
   Max    0.782341
 
610
   Res    0.448767    sqrt( Sum f_i^2 / 3N )
 
611
----------------------------------------
 
612
   Max    0.782341    constrained
 
613
 
 
614
Stress-tensor-Voigt (kbar):     -479.81     -479.67     -479.81       -0.00       -0.00      -22.98
 
615
Target enthalpy (eV/cell)     -194.6894
 
616
 
 
617
* Maximum dynamic memory allocated =     4 MB
 
618
 
 
619
siesta:                 ==============================
 
620
                            Begin FC step =     10
 
621
                            displace atom        2
 
622
                            in direction         2
 
623
                            by         0.0400 Bohr
 
624
                        ==============================
 
625
 
 
626
outcell: Unit cell vectors (Ang):
 
627
        0.000000    2.715000    2.715000
 
628
        2.715000    0.000000    2.715000
 
629
        2.715000    2.715000    0.000000
 
630
 
 
631
outcell: Cell vector modules (Ang)   :    3.839590    3.839590    3.839590
 
632
outcell: Cell angles (23,13,12) (deg):     60.0000     60.0000     60.0000
 
633
outcell: Cell volume (Ang**3)        :     40.0258
 
634
 
 
635
siesta: Eharris(eV) =    -194.932993
 
636
 
 
637
 
 
638
siesta: Eharris(eV) =    -194.932993
 
639
 
 
640
 
 
641
siesta: Atomic forces (eV/Ang):
 
642
----------------------------------------
 
643
   Tot    0.000000    0.010141    0.000001
 
644
----------------------------------------
 
645
   Max    0.782341
 
646
   Res    0.448767    sqrt( Sum f_i^2 / 3N )
 
647
----------------------------------------
 
648
   Max    0.782341    constrained
 
649
 
 
650
Stress-tensor-Voigt (kbar):     -479.81     -479.67     -479.81       -0.00        0.00       22.98
 
651
Target enthalpy (eV/cell)     -194.6894
 
652
 
 
653
* Maximum dynamic memory allocated =     4 MB
 
654
 
 
655
siesta:                 ==============================
 
656
                            Begin FC step =     11
 
657
                            displace atom        2
 
658
                            in direction         3
 
659
                            by        -0.0400 Bohr
 
660
                        ==============================
 
661
 
 
662
outcell: Unit cell vectors (Ang):
 
663
        0.000000    2.715000    2.715000
 
664
        2.715000    0.000000    2.715000
 
665
        2.715000    2.715000    0.000000
 
666
 
 
667
outcell: Cell vector modules (Ang)   :    3.839590    3.839590    3.839590
 
668
outcell: Cell angles (23,13,12) (deg):     60.0000     60.0000     60.0000
 
669
outcell: Cell volume (Ang**3)        :     40.0258
 
670
 
 
671
siesta: Eharris(eV) =    -194.932993
 
672
 
 
673
 
 
674
siesta: Eharris(eV) =    -194.932993
 
675
 
 
676
 
 
677
siesta: Atomic forces (eV/Ang):
 
678
----------------------------------------
 
679
   Tot    0.000001   -0.000000   -0.010142
 
680
----------------------------------------
 
681
   Max    0.782342
 
682
   Res    0.448767    sqrt( Sum f_i^2 / 3N )
 
683
----------------------------------------
 
684
   Max    0.782342    constrained
 
685
 
 
686
Stress-tensor-Voigt (kbar):     -479.81     -479.81     -479.67      -22.98       -0.00        0.00
 
687
Target enthalpy (eV/cell)     -194.6894
 
688
 
 
689
* Maximum dynamic memory allocated =     4 MB
 
690
 
 
691
siesta:                 ==============================
 
692
                            Begin FC step =     12
 
693
                            displace atom        2
 
694
                            in direction         3
 
695
                            by         0.0400 Bohr
 
696
                        ==============================
 
697
 
 
698
outcell: Unit cell vectors (Ang):
 
699
        0.000000    2.715000    2.715000
 
700
        2.715000    0.000000    2.715000
 
701
        2.715000    2.715000    0.000000
 
702
 
 
703
outcell: Cell vector modules (Ang)   :    3.839590    3.839590    3.839590
 
704
outcell: Cell angles (23,13,12) (deg):     60.0000     60.0000     60.0000
 
705
outcell: Cell volume (Ang**3)        :     40.0258
 
706
 
 
707
siesta: Eharris(eV) =    -194.932993
 
708
 
 
709
 
 
710
siesta: Eharris(eV) =    -194.932993
 
711
 
 
712
 
 
713
siesta: Atomic forces (eV/Ang):
 
714
----------------------------------------
 
715
   Tot    0.000001    0.000001    0.010141
 
716
----------------------------------------
 
717
   Max    0.782341
 
718
   Res    0.448767    sqrt( Sum f_i^2 / 3N )
 
719
----------------------------------------
 
720
   Max    0.782341    constrained
 
721
 
 
722
Stress-tensor-Voigt (kbar):     -479.81     -479.81     -479.67       22.98        0.00       -0.00
 
723
Target enthalpy (eV/cell)     -194.6894
 
724
 
 
725
* Maximum dynamic memory allocated =     4 MB
 
726
 
 
727
siesta: Program's energy decomposition (eV):
 
728
siesta: Eions   =       380.802124
 
729
siesta: Ena     =       115.501483
 
730
siesta: Ekin    =       115.577902
 
731
siesta: Enl     =        29.935837
 
732
siesta: DEna    =        -8.120549
 
733
siesta: DUscf   =         1.765896
 
734
siesta: DUext   =         0.000000
 
735
siesta: Exc     =       -68.547879
 
736
siesta: eta*DQ  =         0.000000
 
737
siesta: Emadel  =         0.000000
 
738
siesta: Emeta   =         0.000000
 
739
siesta: Emolmec =         0.000000
 
740
siesta: Ekinion =         0.000000
 
741
siesta: Eharris =      -194.932993
 
742
siesta: Etot    =      -194.689433
 
743
siesta: FreeEng =      -194.689433
 
744
 
 
745
siesta: Final energy (eV):
 
746
siesta:       Kinetic =     115.577902
 
747
siesta:       Hartree =      23.678156
 
748
siesta:    Ext. field =       0.000000
 
749
siesta:   Exch.-corr. =     -68.547879
 
750
siesta:  Ion-electron =    -120.258615
 
751
siesta:       Ion-ion =    -145.138997
 
752
siesta:       Ekinion =       0.000000
 
753
siesta:         Total =    -194.689433
 
754
 
 
755
siesta: Atomic forces (eV/Ang):
 
756
siesta:      1    0.000001    0.000000    0.782341
 
757
siesta:      2   -0.000000    0.000001   -0.772200
 
758
siesta: ----------------------------------------
 
759
siesta:    Tot    0.000001    0.000001    0.010141
 
760
 
 
761
siesta: Stress tensor (static) (eV/Ang**3):
 
762
siesta:    -0.299468    0.014341   -0.000000
 
763
siesta:     0.014341   -0.299468    0.000000
 
764
siesta:    -0.000000    0.000000   -0.299383
 
765
 
 
766
siesta: Cell volume =         40.025752 Ang**3
 
767
 
 
768
siesta: Pressure (static):
 
769
siesta:                Solid            Molecule  Units
 
770
siesta:           0.00326128          0.00335847  Ry/Bohr**3
 
771
siesta:           0.29943952          0.30836375  eV/Ang**3
 
772
siesta:         479.76018441        494.05853840  kBar
 
773
 
 
774
* Maximum dynamic memory allocated : Node    0 =     4 MB
 
775
 
 
776
* Maximum memory occured during ranger                        
 
777
 
 
778
timer: CPU execution times:
 
779
timer:  Routine       Calls   Time/call    Tot.time        %
 
780
timer:  siesta            1      23.294      23.294   100.00
 
781
timer:  Setup             1       0.338       0.338     1.45
 
782
timer:  bands             1       0.001       0.001     0.00
 
783
timer:  writewave         1       0.003       0.003     0.01
 
784
timer:  KSV_init          1       0.000       0.000     0.00
 
785
timer:  IterMD           13       1.765      22.941    98.48
 
786
timer:  hsparse          14       0.003       0.049     0.21
 
787
timer:  overfsm          26       0.003       0.090     0.39
 
788
timer:  IterSCF          13       0.423       5.496    23.59
 
789
timer:  kinefsm          26       0.003       0.084     0.36
 
790
timer:  nlefsm           26       0.016       0.409     1.76
 
791
timer:  DHSCF            26       0.853      22.184    95.23
 
792
timer:  DHSCF1            1       0.130       0.130     0.56
 
793
timer:  DHSCF2           13       0.249       3.242    13.91
 
794
timer:  REORD           234       0.000       0.022     0.09
 
795
timer:  POISON           39       0.004       0.159     0.68
 
796
timer:  DHSCF3           26       0.138       3.581    15.37
 
797
timer:  rhoofd           26       0.065       1.702     7.31
 
798
timer:  cellXC           39       0.032       1.247     5.35
 
799
timer:  vmat             26       0.034       0.885     3.80
 
800
timer:  diagon           13       0.001       0.010     0.04
 
801
timer:  rdiag            13       0.000       0.006     0.03
 
802
timer:  rdiag1           13       0.000       0.001     0.00
 
803
timer:  rdiag2           13       0.000       0.002     0.01
 
804
timer:  rdiag3           13       0.000       0.000     0.00
 
805
timer:  rdiag4           13       0.000       0.000     0.00
 
806
timer:  MolMec           26       0.000       0.000     0.00
 
807
timer:  DHSCF4           13       1.171      15.223    65.35
 
808
timer:  dfscf            13       0.908      11.800    50.65
 
809
timer:  optical           1       0.000       0.000     0.00
 
810
  
 
811
 
 
812
elaps: ELAPSED times:
 
813
elaps:  Routine       Calls   Time/call    Tot.time        %
 
814
elaps:  siesta            1      23.299      23.299   100.00
 
815
elaps:  Setup             1       0.339       0.339     1.46
 
816
elaps:  bands             1       0.001       0.001     0.01
 
817
elaps:  writewave         1       0.002       0.002     0.01
 
818
elaps:  KSV_init          1       0.001       0.001     0.00
 
819
elaps:  IterMD           13       1.765      22.943    98.47
 
820
elaps:  hsparse          14       0.003       0.047     0.20
 
821
elaps:  overfsm          26       0.004       0.092     0.39
 
822
elaps:  IterSCF          13       0.423       5.495    23.59
 
823
elaps:  kinefsm          26       0.003       0.082     0.35
 
824
elaps:  nlefsm           26       0.016       0.410     1.76
 
825
elaps:  DHSCF            26       0.853      22.189    95.23
 
826
elaps:  DHSCF1            1       0.130       0.130     0.56
 
827
elaps:  DHSCF2           13       0.250       3.244    13.92
 
828
elaps:  REORD           234       0.000       0.025     0.11
 
829
elaps:  POISON           39       0.004       0.162     0.70
 
830
elaps:  DHSCF3           26       0.138       3.582    15.38
 
831
elaps:  rhoofd           26       0.065       1.699     7.29
 
832
elaps:  cellXC           39       0.032       1.248     5.36
 
833
elaps:  vmat             26       0.034       0.888     3.81
 
834
elaps:  diagon           13       0.001       0.009     0.04
 
835
elaps:  rdiag            13       0.000       0.006     0.03
 
836
elaps:  rdiag1           13       0.000       0.000     0.00
 
837
elaps:  rdiag2           13       0.000       0.000     0.00
 
838
elaps:  rdiag3           13       0.000       0.001     0.01
 
839
elaps:  rdiag4           13       0.000       0.000     0.00
 
840
elaps:  MolMec           26       0.000       0.000     0.00
 
841
elaps:  DHSCF4           13       1.171      15.223    65.34
 
842
elaps:  dfscf            13       0.908      11.800    50.65
 
843
elaps:  optical           1       0.001       0.001     0.00
 
844
  
 
845
>> End of run:  29-JAN-2007  17:58:36