CalculationMode = td
FromScratch = yes

ExperimentalFeatures = yes
PeriodicDimensions = 1

Spacing = 0.4

%LatticeParameters
2.0*angstrom | 3.5*angstrom | 3.5*angstrom
%

PseudopotentialSet = pseudodojo_pbe
%Coordinates
"H" | -0.3707*angstrom | 0.0 | 0.0
"H" | 0.3707*angstrom | 0.0 | 0.0
%

%KPointsGrid
25 | 1 | 1
%
KPointsUseSymmetries = yes
%SymmetryBreakDir
1 | 0 | 0
%

omega = 1.5*eV
period = 2*pi/omega
stime = 4*period

TDPropagationTime = stime+period
TDPropagator = exp_mid
TDExponentialMethod = lanczos
TDTimeStep = 0.1

%TDExternalFields
vector_potential | 1 | 0 | 0 | omega | "envelope_function"
%

electric_amplitude = -c/omega*(sqrt(1*10^10)/sqrt(3.509470*10^16)) # in atomic units
%TDFunctions
"envelope_function" | tdf_from_expr | "electric_amplitude*(sin(pi/stime*t))^2*(1-step(t-stime))"
%

RestartWriteInterval = 1000

%TDOutput
laser
total_current
%
