next up previous contents index
Next: Keywords Up: Input file Previous: Input file   Contents   Index

An example: methane molecule

An input file Methane.dat in the directory 'work' is shown below. This input file has a flexible data format, in which a parameter is given behind a keyword, the order of keywords is arbitrary, and a blank and a comment can be also described freely. For the keywords and options, both capital, small letters, and the mixture are acceptable, although these options in below example are written in a specific form.

  #
  #  SCF calculation of a methane molecule by the LDA
  #  and the cluster method    
  #

  #
  # File Name      
  #

  System.CurrrentDirectory         ./    # default=./
  System.Name                     met
  level.of.stdout                   1    # default=1 (1-3)
  level.of.fileout                  1    # default=1 (0-2)

  #
  # Definition of Atomic Species  
  #

  Species.Number       2
  <Definition.of.Atomic.Species
   H   H4.0-s1          H_TM
   C   C4.5-s1p1        C_TM_PCC
  Definition.of.Atomic.Species>

  #
  # Atoms
  #

  Atoms.Number         5
  Atoms.SpeciesAndCoordinates.Unit   Ang # Ang|AU|FRAC
  <Atoms.SpeciesAndCoordinates
   1  C      0.000000    0.000000    0.000000     2.0  2.0 
   2  H     -0.889981   -0.629312    0.000000     0.5  0.5
   3  H      0.000000    0.629312   -0.889981     0.5  0.5
   4  H      0.000000    0.629312    0.889981     0.5  0.5
   5  H      0.889981   -0.629312    0.000000     0.5  0.5
  Atoms.SpeciesAndCoordinates>
  Atoms.UnitVectors.Unit             Ang # Ang|AU
  #<Atoms.UnitVectors                     
  #  10.0   0.0   0.0  
  #   0.0  10.0   0.0  
  #   0.0   0.0  10.0  
  #Atoms.UnitVectors>

  #
  # SCF or Electronic System
  #

  scf.XcType                  LDA        # LDA|LSDA-CA|LSDA-PW|GGA-PBE
  scf.SpinPolarization        off        # On|Off|NC
  scf.ElectronicTemperature  100.0       # default=300 (K)
  scf.energycutoff           120.0       # default=150 (Ry)
  scf.maxIter                 60         # default=40
  scf.EigenvalueSolver        cluster    # DC|GDC|Cluster|Band
  scf.Kgrid                  1 1 1       # means n1 x n2 x n3
  scf.Mixing.Type           rmm-diis     # Simple|Rmm-Diis|Gr-Pulay|Kerker|Rmm-Diisk
  scf.Init.Mixing.Weight     0.30        # default=0.30 
  scf.Min.Mixing.Weight      0.001       # default=0.001 
  scf.Max.Mixing.Weight      0.400       # default=0.40 
  scf.Mixing.History          7          # default=5
  scf.Mixing.StartPulay       4          # default=6
  scf.criterion             1.0e-8       # default=1.0e-6 (Hartree) 
  scf.lapack.dste            dstevx      # dstevx|dstedc|dstegr,default=dstevx

  #
  # 1D FFT  
  #

  1DFFT.NumGridK             900         # default=900
  1DFFT.NumGridR             900         # default=900
  1DFFT.EnergyCutoff        2500.0       # default=3600 (Ry)

  #
  # Orbital Optimization
  #

  orbitalOpt.Method           off        # Off|Unrestricted|Restricted
  orbitalOpt.InitCoes     Symmetrical    # Symmetrical|Free
  orbitalOpt.initPrefactor   0.1         # default=0.1
  orbitalOpt.scf.maxIter      15         # default=12
  orbitalOpt.MD.maxIter        7         # default=5
  orbitalOpt.per.MDIter       20         # default=1000000
  orbitalOpt.criterion      1.0e-4       # default=1.0e-4 (Hartree/borh)

  # 
  # output of contracted orbitals
  #

  CntOrb.fileout               off       # on|off, default=off
  Num.CntOrb.Atoms             1         # default=1
  <Atoms.Cont.Orbitals
   1
  Atoms.Cont.Orbitals>
 
  #
  # SCF Order-N
  #

  orderN.HoppingRanges        6.0        # default=5.0 (Ang) 
  orderN.NumHoppings           2         # default=2

  #
  # MD or Geometry Optimization
  #

  MD.Type                     nomd       # Nomd|Opt|NVE|NVT_VS|NVT_NH
                                       # Constraint_Opt|DIIS
  MD.Opt.DIIS.History          4
  MD.Opt.StartDIIS             5         # default=5

  MD.maxIter                   1         # default=1
  MD.TimeStep                 1.0        # default=0.5 (fs)
  MD.Opt.criterion         1.0e-5        # default=1.0e-4 (Hartree/bohr)

  #
  # restart using a restart file, *.rst  
  #

  scf.restart                 off        # on|off,default=off

  #
  # MO output
  #

  MO.fileout                   off       # on|off
  num.HOMOs                     1        # default=1
  num.LUMOs                     1        # default=1
  MO.Nkpoint                    1        # default=1 
  <MO.kpoint
    0.0  0.0  0.0
  MO.kpoint>

  #
  # DOS and PDOS
  #

  Dos.fileout                  off       # on|off, default=off
  Dos.Erange               -10.0  10.0   # default = -20 20 
  Dos.Kgrid                  1  1  1     # default = Kgrid1 Kgrid2 Kgrid3

  #
  # output Hamiltonian and overlap
  #  

  HS.fileout                   off       # on|off, default=off



2009-08-28