mirror of https://github.com/abinit/abinit.git
97 lines
3.9 KiB
Plaintext
97 lines
3.9 KiB
Plaintext
# Input for RT-TDDFT tutorial
|
|
# RT-TDDFT calculation
|
|
# Diamond at equilibrium volume
|
|
# First dataset: ground state calculation
|
|
# Second dataset: RT-TDDFT with impulsive electric field
|
|
|
|
ndtset 2
|
|
|
|
# RT-TDDFT section
|
|
optdriver2 9
|
|
getwfk2 -1 # Read initial orbitals from the first dataset
|
|
dtele 0.10 # Electronic time step
|
|
ntime 100 # Maximum number or time step
|
|
td_ef_type 1 # Impulse electric field
|
|
td_ef_tzero 1.0 # applied at time tzero
|
|
|
|
#Definition of the unit cell
|
|
acell 3*3.53 angstrom # Lengths of the primitive vectors (equilibrium volume)
|
|
rprim # 3 orthogonal primitive vectors (FCC lattice)
|
|
0.0 1/2 1/2
|
|
1/2 0.0 1/2
|
|
1/2 1/2 0.0
|
|
nsym 0 # Automatic detection of symetries
|
|
|
|
#Definition of the atom types and pseudopotentials
|
|
ntypat 1 # There is only one type of atom
|
|
znucl 6 # Atomic number of the possible type(s) of atom. Here carbon.
|
|
pp_dirpath "$ABI_PSPDIR" # Path to the directory where
|
|
# pseudopotentials for tests are stored
|
|
pseudos "Psdj_paw_pw_std/C.xml" # PAW atomic data for Carbon
|
|
|
|
#Definition of the atoms
|
|
natom 2 # There are two atoms
|
|
typat 1 1 # They both are of type 1, that is, Carbon
|
|
xred # Location of the atoms: Diamond structure
|
|
0.0 0.0 0.0 # Triplet giving the reduced coordinates of atom 1
|
|
1/4 1/4 1/4 # Triplet giving the reduced coordinates of atom 2
|
|
|
|
#Definition of bands and occupation numbers
|
|
nband 8 # Compute 8 bands (4 occupied, 4 empty)
|
|
nbdbuf 2 # No need to converge the last two empty bands
|
|
occopt 1 # Automatic generation of occupation numbers, as a semiconductor
|
|
|
|
#Numerical parameters of the calculation : planewave basis set and k point grid
|
|
ecut 18 # Maximal plane-wave kinetic energy cut-off, in Hartree
|
|
pawecutdg 36 # Max. plane-wave kinetic energy cut-off, in Ha, for the PAW double grid
|
|
kptopt 3 # Automatic generation of k points, no taking symmetry into account
|
|
ngkpt 2 2 2 # k-points grid based on the primitive vectors
|
|
nshiftk 4 # of the reciprocal space, repeated four times,
|
|
shiftk # with different shifts:
|
|
0.5 0.5 0.5
|
|
0.5 0.0 0.0
|
|
0.0 0.5 0.0
|
|
0.0 0.0 0.5
|
|
istwfk *1 # Wavefunctions are complex for all k-points
|
|
|
|
#Parameters for the SCF procedure
|
|
nstep 20 # Maximal number of SCF cycles
|
|
tolvrs 1.0d-10 # Will stop when, twice in a row, the difference
|
|
# between two consecutive evaluations of potential residual
|
|
# differ by less than tolvrs
|
|
tolwfr 1.0d-20 # Additional convergence criterion on the largest residual
|
|
# to ensure that wavefunctions are also tighlty converged
|
|
# as they will be used as initial conditions for the following
|
|
# RT-TDDFT run
|
|
|
|
#Miscelaneous parameters
|
|
prtwf1 1 # Print wavefunctions (required for following RT-TDDFT runs)
|
|
prtwf2 0 # Do not print wavefunctions at each step during RT-TDDFT
|
|
prtden 0 # Do not print density
|
|
prteig 0 # Do not print eigenvalues
|
|
prtebands 0 # Do not print ebands
|
|
prtcurrent 1 # Print current
|
|
|
|
|
|
##############################################################
|
|
# This section is used only for regression testing of ABINIT #
|
|
##############################################################
|
|
#%%<BEGIN TEST_INFO>
|
|
#%% [setup]
|
|
#%% executable = abinit
|
|
#%% [files]
|
|
#%% files_to_test =
|
|
#%% trttddft_4.abo, tolnlines= 0, tolabs= 0.0, tolrel= 0.0, fld_options= -medium ;
|
|
#%% trttddft_4o_DS2_TDENER, tolnlines= 0, tolabs= 0.0, tolrel= 0.0, fld_options= -medium ;
|
|
#%% trttddft_4o_DS2_TDCURRENT, tolnlines= 0, tolabs= 0.0, tolrel= 0.0, fld_options -=medium ;
|
|
#%% [paral_info]
|
|
#%% max_nprocs = 8
|
|
#%% [extra_info]
|
|
#%% authors = F. Brieuc
|
|
#%% keywords = RTTDDFT, PAW
|
|
#%% description =
|
|
#%% Input for RT-TDDFT tutorial
|
|
#%% Diamond at equilibrium volume
|
|
#%% RT-TDDFT calculation with TD electric field
|
|
#%%<END TEST_INFO>
|