Difference between revisions of "Explanations for C2H2 rt pulse input files"
| Line 12: | Line 12: | ||
If you do not specify the units for some physical quantities, atomic unit will be used for those quantities. | If you do not specify the units for some physical quantities, atomic unit will be used for those quantities. | ||
| − | == &calculation | + | == &calculation == |
| + | |||
| + | Mandatory: calc_mode | ||
&calculation | &calculation | ||
| Line 22: | Line 24: | ||
For periodic systems (specified by <code>iperiodic = 3</code> in <code>&system</code>), both ground state and real time calculations should be carried out as a single task (<code>calc_mode = 'GS_RT'</code>). | For periodic systems (specified by <code>iperiodic = 3</code> in <code>&system</code>), both ground state and real time calculations should be carried out as a single task (<code>calc_mode = 'GS_RT'</code>). | ||
| − | == &control | + | == &control == |
| + | |||
| + | Mandatory: none | ||
&control | &control | ||
| Line 30: | Line 34: | ||
'C2H2' defined by <code>surname = 'C2H2'</code> will be used in the filenames of output files. | 'C2H2' defined by <code>surname = 'C2H2'</code> will be used in the filenames of output files. | ||
| − | == ¶llel | + | == ¶llel == |
| + | |||
| + | Mandatory: none? | ||
¶llel | ¶llel | ||
| Line 44: | Line 50: | ||
<code>nproc_domain_s = 1,1,1</code> indicates the spatial divisions for Hartree potential in x,y,z directions. | <code>nproc_domain_s = 1,1,1</code> indicates the spatial divisions for Hartree potential in x,y,z directions. | ||
| − | == &tgrid | + | == &tgrid == |
| + | |||
| + | Mandatory: dt, Nt | ||
&tgrid | &tgrid | ||
| Line 54: | Line 62: | ||
<code>Nt=24000</code> specifies the number of time steps in the calculation. | <code>Nt=24000</code> specifies the number of time steps in the calculation. | ||
| − | == &emfield | + | == &emfield == |
| + | |||
| + | Mandatory: ae_shape1, epdir_re1, {rlaser_int1 or amplitude1}, omega1, pulse_tw1, phi_cep1 | ||
&emfield | &emfield | ||
| Line 79: | Line 89: | ||
The time delay can be specified using the variable 't1_t2'. | The time delay can be specified using the variable 't1_t2'. | ||
| − | == &hartree | + | == &hartree == |
| + | |||
| + | Mandatory: none | ||
&hartree | &hartree | ||
Revision as of 13:49, 11 June 2017
&unit
Mandatory: none
&units unit_length='Angstrom' unit_energy='eV' unit_time='fs' /
This namelist specifies the unit system used in the input and the output files. If you do not specify the units for some physical quantities, atomic unit will be used for those quantities.
&calculation
Mandatory: calc_mode
&calculation calc_mode = 'RT' /
The variable calc_mode should be one of 'GS', 'RT', and 'GS-RT'.
Note that the ground state ('GS') and real time ('RT') calculations should be done separately and sequentially for isolated systems (specified by iperiodic = 0 in &system).
For periodic systems (specified by iperiodic = 3 in &system), both ground state and real time calculations should be carried out as a single task (calc_mode = 'GS_RT').
&control
Mandatory: none
&control sysname = 'C2H2' /
'C2H2' defined by surname = 'C2H2' will be used in the filenames of output files.
¶llel
Mandatory: none?
¶llel domain_parallel = 'y' nproc_ob = 1 nproc_domain = 1,1,1 nproc_domain_s = 1,1,1 /
domain_parallel = 'y' indicates that the spatial grid is divided and parallely executed.
nproc_ob = 1 indicates the number of MPI parallelization for orbitals.
nproc_domain = 1,1,1 indicates the spatial division for orbitals in x,y,z directions.
nproc_domain_s = 1,1,1 indicates the spatial divisions for Hartree potential in x,y,z directions.
&tgrid
Mandatory: dt, Nt
&tgrid dt=1.25d-3 Nt=24000
/
dt=1.25d-3 specifies the time step of the time evolution calculations.
Nt=24000 specifies the number of time steps in the calculation.
&emfield
Mandatory: ae_shape1, epdir_re1, {rlaser_int1 or amplitude1}, omega1, pulse_tw1, phi_cep1
&emfield ae_shape1 = 'esin2cos' epdir_re1 = 0.d0,0.d0,1.d0 rlaser_int1 = 1.d8 omega1=9.26d0 pulse_tw1=30.d0 phi_cep1=0.75d0 /
ae_shape1 = 'esin2cos' specifies the pulse shape of the electric field, having sin-square envelope.
epdir_re1 = 0.d0,0.d0,1.d0 specifies the unit vector of polarization direction.
Specifying real part of the polarization vector by 'epdir_re1', linear polarization is assumed.
Using both real ('epdir_re1') and imaginary ('epdir_im1') parts of the polarization vector, circularly (and general ellipsoidary) polarized pulses may also be described.
laser_int1 = 1.d8 specifies maximum intensity of the applied electric field in unit of W/cm^2.
omega1=9.26d0 specifies photon energy (frequency multiplied with hbar).
pulse_tw1=30.d0 specifies the pulse duration.
Note that it is not FWHM but a full duration of the sin-square envelope.
phi_cep1=0.75d0 specifies the carrier envelope phase of the pulse.
Note that 'phi_cep' must be 0.75 (or 0.25) if one employs 'esin2cos' pulse shape, since otherwise the time integral of the electric field does not vanish.
It is possible to use two pulses simultaneously to simulate pump-probe experiments, adding information for two pulses. The time delay can be specified using the variable 't1_t2'.
&hartree
Mandatory: none
&hartree meo = 3 num_pole_xyz = 2,2,2 /
meo = 3 specifies the order of multipole expansion of electron density that is used to prepare boundary condition for the Hartree potential.
