~nickpapior/siesta/trunk-kovalp

« back to all changes in this revision

Viewing changes to Src/Tests/var_cell/var_cell.fdf

  • Committer: Alberto Garcia
  • Date: 2005-10-04 13:49:31 UTC
  • mfrom: (unknown (missing))
  • Revision ID: Arch-1:siesta@uam.es--2005%siesta-devel--reference--1.4--patch-13
New kind of structure input -- constant-volume variable-cell - no sig
* Using the option UseStructFile, the structural information will be
read from SiestaLabel.STRUCT_IN. This option is incompatible with
UseSaveXV, and will prevail.

* Siesta now always produces a SiestaLabel.STRUCT_OUT file.

The .STRUCT_XX files have the cell vectors in Ang and the atomic
positions in atomic coordinates.

* The option MD.ConstantVolume will result in a constant-volume
variable-cell simulation (i.e., only the cell shape and the atomic
positions change).

* After every MD or relaxation step, Siesta now writes the "enthalpy"
E+pV, where p is actually the target pressure.

* New tests: var_cell and constant_volume.

* Changed si test to use the new UseStructFile feature.

* Removed "sig" fossils from source.


Show diffs side-by-side

added added

removed removed

Lines of Context:
 
1
#
 
2
# MgCO3 in primitive cell. SZ. Variable cell
 
3
#
 
4
SystemName      MgCo3 R-3c test -- SZ, 100 Ry -- variable cell at 50 Gpa
 
5
SystemLabel     var_cell
 
6
NumberOfSpecies         3
 
7
NumberOfAtoms          10
 
8
%block ChemicalSpeciesLabel
 
9
      1      12     Mg
 
10
      2       6     C
 
11
      3       8     O
 
12
%endblock ChemicalSpeciesLabel
 
13
 
 
14
PAO.BasisSize  SZ
 
15
 
 
16
# Rhombohedral primitive cell
 
17
# (HEX 3-fold cell had A=4.635 and C=15.023, so alpha=48.179 degrees
 
18
#  and a=5.67783 Ang)
 
19
 
 
20
LatticeConstant     5.67783 Ang
 
21
%block LatticeParameters   
 
22
1.0 1.0 1.0 48.179 48.179 48.179
 
23
%endblock LatticeParameters   
 
24
 
 
25
AtomicCoordinatesFormat    Fractional
 
26
 
 
27
%block AtomicCoordinatesAndAtomicSpecies
 
28
  0.0         0.0         0.0      1
 
29
  0.5         0.5         0.5      1
 
30
  0.25        0.25        0.25     2
 
31
 -0.25       -0.25       -0.25     2
 
32
  0.5274     -0.0274      0.25     3
 
33
  0.25        0.5274     -0.0274   3
 
34
 -0.0274      0.25        0.5274   3
 
35
 -0.5274      0.0274     -0.25     3
 
36
 -0.25       -0.5274      0.0274   3
 
37
  0.0274     -0.25       -0.5274   3
 
38
%endblock AtomicCoordinatesAndAtomicSpecies
 
39
 
 
40
kgrid_cutoff 15.0 Ang
 
41
 
 
42
Solution.Method       diagon
 
43
MeshCutoff             100 Ry
 
44
DM.NumberBroyden       4
 
45
DM.UseSaveDM           T
 
46
DM.MixingWeight        0.1         # New DM amount for next SCF cycle
 
47
DM.Tolerance           1.d-4       # Tolerance in maximum difference
 
48
                                   # between input and output DM
 
49
MaxSCFIterations       20
 
50
 
 
51
WriteCoorStep      .true.
 
52
WriteForces        .true.
 
53
 
 
54
MD.TypeOfRun         CG
 
55
MD.Variable-Cell      T
 
56
MD.Target-pressure    50.0 GPa
 
57
MD.Num-CG-steps      5
 
58
MD.Max-Stress-Tol    2.0 GPa
 
59
 
 
60