Explanations of input files (dielectric function of crystalline silicon)

From salmon
Revision as of 14:57, 23 November 2017 by Yabana (talk | contribs) (Created page with "=== &calculation === Mandatory: calc_mode &calculation calc_mode = 'GS_RT' / This indicates that the ground state (GS) and the real time (RT) calculations are carried...")
(diff) ← Older revision | Latest revision (diff) | Newer revision → (diff)
Jump to: navigation, search

&calculation

Mandatory: calc_mode

&calculation
  calc_mode = 'GS_RT'
/

This indicates that the ground state (GS) and the real time (RT) calculations are carried out sequentially in the present job. See Tutorial-v.1.0.0#&calculation for detail.

&control

Mandatory: none

&control
  sysname = 'Si'
/

'Si' defined by sysname = 'C2H2' will be used in the filenames of output files.

&system

Mandatory: periodic, al, state, nelem, natom

&system
  iperiodic = 3
  al = 10.26d0,10.26d0,10.26d0
  nstate = 32
  nelec = 32
  nelem = 1
  natom = 8
/

iperiodic = 3 indicates that three dimensional periodic boundary condition (bulk crystal) is assumed. al = 10.26d0, 10.26d0, 10.26d0 specifies the lattice constans of the unit cell. nstate = 32 indicates the number of Kohn-Sham orbitals to be solved. nelec = 32 indicate the number of valence electrons in the system. nelem = 1 and natom = 8 indicate the number of elements and the number of atoms in the system, respectively. See Tutorial-v.1.0.0#&system for more information.

&pseudo

&pseudo
  izatom(1)=14
  pseudo_file(1) = './Si_rps.dat'
  lloc_ps(1)=2
/

izatom(1) = 14 indicates the atomic number of the element #1. pseudo_file(1) = 'Si_rps.dat' indicates the pseudopotential filename of element #1. lloc_ps(1) = 2 indicate the angular momentum of the pseudopotential that will be treated as local.

&functional

&functional
  xc = 'PZ'
/

This indicates that the adiabatic local density approximation with the Perdew-Zunger functional is used. We note that meta-GGA functionals that reasonably reproduce the band gap of various insulators may also be used in the calculation of periodic systems. See Tutorial-v.1.0.0#&functional for detail.

&rgrid

Mandatory: dl or num_rgrid

&rgrid
  num_rgrid = 12,12,12
/

num_rgrid=12,12,12 specifies the number of the grids for each Cartesian direction. See Tutorial-v.1.0.0#&rgrid for more information.

&kgrid

Mandatory: none

This namelist provides grid spacing of k-space for periodic systems.

&kgrid
  num_kgrid = 4,4,4
/

&tgrid

&tgrid
 nt=3000
 dt=0.16  
/

dt=0.16 specifies the time step of the time evolution calculation. nt=3000 specifies the number of time steps in the calculation.

&propagation

&propagation
  propagator='etrs'
/

propagator = 'etrs' indicates the use of enforced time-reversal symmetry propagator. See Tutorial-v.1.0.0#&propagation for more information.

&scf

Mandatory: nscf

This namelists specify parameters related to the self-consistent field calculation.

&scf
  ncg = 5
  nscf = 120
/

ncg = 5 is the number of conjugate-gradient iterations in solving the Kohn-Sham equation. Usually this value should be 4 or 5. nscf = 120 is the number of scf iterations.

&emfield

Mandatory:ae_shape1

&emfield
  trans_longi = 'tr'
  ae_shape1 = 'impulse'
  epdir_re1 = 0.,0.,1.
/

as_shape1 = 'impulse' indicates that a weak impulsive field is applied to all electrons at t=0

epdir_re1(3) specify a unit vector that indicates the direction of the impulse.

trans_longi = 'tr' specifies the treatment of the polarization in the time evolution calculation, transverse for 'tr' and longitudinal for 'lo'.

See Tutorial-v.1.0.0#&emfield for detail.

&analysis

&analysis
 nenergy=1000
 de=0.001
/

nenergy=1000 specifies the number of energy steps, and de=0.001 specifies the energy spacing in the time-frequency Fourier transformation.

&atomic_red_coor

Mandatory: atomic_coor or atomic_red_coor (they may be provided as a separate file)

&atomic_red_coor
 'Si'    .0      .0      .0      1
 'Si'    .25    .25    .25    1
 'Si'    .5      .0      .5      1
 'Si'    .0      .5      .5      1
 'Si'    .5      .5      .0      1
 'Si'    .75    .25   .75     1
 'Si'    .25    .75   .75     1
 'Si'    .75    .75   .25     1
/

Cartesian coordinates of atoms are specified in a reduced coordinate system. First column indicates the element, next three columns specify reduced Cartesian coordinates of the atoms, and the last column labels the element.