Job control

| prefix | calc_mode | fklist | fqlist | ftemper | debug | yaml_fname | hole | tmp_dir | load_scatter_eph | find_efermi | sampling | cauchy_scale | scat_impl | nsamples | use_mem | cum_inner_emin | cum_inner_emax | cum_outer_emin | cum_outer_emax | cum_de | spectral_emin | spectral_emax | cum_outer_np | spectral_np |

Boltzmann Transport Equation

| boltz_kdim | boltz_qdim | band_min | band_max | boltz_emin | boltz_emax | boltz_nstep | boltz_de | delta_smear | full_ite | phfreq_cutoff | trans_thr | lmagsym |

Polar correction

| polar_split |

Ultra-fast dynamics

| time_step | output_nstep | boltz_init_dist | pump_pulse | pump_pulse_fname | boltz_init_e0 | boltz_init_smear | boltz_init_ampl | solver | boltz_efield | boltz_norm_dist | boltz_acc_thr | boltz_nstep_min | ph_mode_exclude_ranges |


 Job control

JOB CONTROL
prefix

Variable type: string
Default value: ''

Job name prefix. It should be the same as the prefix used in QE.

Mandatory for all calculation types.

JOB CONTROL
calc_mode

Variable type: string
Default value: ''
Options:
'bands'
interpolated electronic band structure (see in the section)
'phdisp'
interpolated phonon dispersions along a given crystal momentum path (see in the section)
'ephmat'
absolute values of the e-ph matrix elements (see in the section)
'setup'
set up transport property calculations (see in the section)
'imsigma'
imaginary part of the lowest-order e-ph self-energy (see in the section)
'meanfp'
e-ph mean free paths for electronic states (see in the section)
'trans-rta'
electrical conductivity and carrier mobility tensors using relaxation time approximation (see in the section)
'trans-ita'
electrical conductivity and carrier mobility tensors using iterative method (see in the section)
'trans-mag-rta'
phonon-limited conductivity and carrier mobility using relaxation time approximation in finite magnetic field (see in the section)
'trans-mag-ita'
phonon-limited conductivity and carrier mobility using full Boltzmann equation in a finite magnetic field (see in the section)
'trans-pp'
computes Seebeck coefficient (see in the section)
'dynamics-run'
ultrafast hot carrier dynamics (see in the section)
'dynamics-pp'
computes the energy-dependent carrier population as a function of time (see in the section)
'ephmat_spin'
absolute values of the e-ph spin-flip matrix elements (see in the section)
'spins'
interpolated spin texture (see in the section)
'imsigma_spin'
imaginary part of the lowest-order e-ph spin-flip self-energy (see in the section)
'spinlifetime'
spin relaxation times (see in the section)
'spectral-re'
on-shell real-part of the Fan-Migdal self-energy (see in the section)
'spectral-se'
frequency-dependent imaginary part of the Fan-Migdal self-energy (see in the section)
'spectral-cum'
cumulant spectral function (see in the section)

Calculation mode. To see the typical input files for different calculation modes, click here.

Mandatory for all calculation types.

JOB CONTROL
fklist

Variable type: string
Default value: ''
Typical value: 'prefix_tet.kpt'

Name of the file containing the k-point list (in crystal coordiates).

Mandatory for calc_mode:  bandsephmatimsigmameanfpephmat_spinspinsimsigma_spin
spectral-respectral-sespectral-cum

JOB CONTROL
fqlist

Variable type: string
Default value: ''
Typical value: 'prefix_phdisp.qpt'

Name of the file containing the q-point list (in crystal coordiates).

Mandatory for calc_mode:  phdispephmatephmat_spin Optional for calc_mode:  imsigmaimsigma_spin

JOB CONTROL
ftemper

Variable type: string
Default value: ''
Typical value: 'prefix.temper'

Name of the file containing values for the temperature (K), chemical potential (eV), and carrier concentration (cm-2 or cm-3).

Mandatory for calc_mode:  setupimsigmameanfptrans-rtatrans-itatrans-mag-rtatrans-mag-ita
trans-ppdynamics-rundynamics-ppimsigma_spinspinlifetimespectral-respectral-sespectral-cum

JOB CONTROL
debug

Variable type: logical
Default value: .false.

Debug mode.

Optional for calc_mode:  qe2pert

JOB CONTROL
yaml_fname

Variable type: string
Default value: '.yml'

Name of the YAML output file. If the parameter (yaml_fname) is not changed by the user, the output name of the YAML file will be {prefix}_{calc_mode}.yml. For example, for Si bands calculation, the YAML file would be called si_bands.yml.

Optional for all calculation types.

JOB CONTROL
hole

Variable type: logical
Default value: .false.

Set to .true. for calculations on hole carriers.

Optional for calc_mode:  setuptrans-rtatrans-itatrans-mag-rtatrans-mag-itatrans-ppdynamics-run
dynamics-pp

JOB CONTROL
tmp_dir

Variable type: string
Default value: '.'
Typical value: './tmp'

The directory where the e-ph matrix elements are stored when calc_mode='trans'.

Optional for calc_mode:  trans-rtatrans-itatrans-mag-rtatrans-mag-itadynamics-run

JOB CONTROL
load_scatter_eph

Variable type: logical
Default value: .false.

Read the e-ph matrix elements from the files in tmp_dir. Used for calc_mode='trans'.

Optional for calc_mode:  trans-rtatrans-itatrans-mag-rtatrans-mag-itadynamics-run

JOB CONTROL
find_efermi

Variable type: logical
Default value: .false.
Typical value: .false.

Flag to indicate whether to find the Fermi energy, which means fix the Fermi energy or occupancy.

Optional for calc_mode:  setup

JOB CONTROL
sampling

Variable type: string
Default value: 'uniform'
Options:
'uniform'
random sampling from an uniform distribution.
'cauchy'
random sampling from a Cauchy disrtibution.

Random q points sampling method.

Optional for calc_mode:  imsigmaimsigma_spinspectral-respectral-sespectral-cum

JOB CONTROL
cauchy_scale

Variable type: real
Default value: 0.05
Typical value: 1.0

Scale parameter gamma for the Cauchy distribution; used when sampling='cauchy'.

Optional for calc_mode:  imsigmaimsigma_spinspectral-respectral-sespectral-cum

JOB CONTROL
scat_impl

Variable type: string
Default value: 'std'
Options:
'std'
Standard CPU-friendly carrier dynamics implementation; reduced memory usage requirements
'tgt'
\"Target-oriented\" GPU-friendly carrier dynamics implementation; increased memory usage requirements

Specify which carrier-dynamics implementation to use, either 'std' for the standard CPU-friendly implementation, or 'tgt' for the GPU-friendly target-oriented implementation.

Mandatory for calc_mode:  trans-itatrans-mag-rtatrans-mag-itadynamics-run

JOB CONTROL
scat_impl

Variable type: string
Default value: 'std'
Options:
'std'
Standard CPU-friendly carrier dynamics implementation; reduced memory usage requirements
'tgt'
\"Target-oriented\" GPU-friendly carrier dynamics implementation; increased memory usage requirements

Specify which carrier-dynamics implementation to use, either 'std' for the standard CPU-friendly implementation, or 'tgt' for the GPU-friendly target-oriented implementation.

Mandatory for calc_mode:  transdynamics-run

JOB CONTROL
nsamples

Variable type: integer
Default value: 100000

Number of q-points for the summation over the q-points in imsigma calculation.

Optional for calc_mode:  imsigmaimsigma_spinspectral-respectral-sespectral-cum

JOB CONTROL
use_mem

Variable type: logical
Default value: .true.
Typical value: .true.

Flag for using memory or not and only used in boltz_scatter


JOB CONTROL
cum_inner_emin

Variable type: real
Default value: -0.1
Typical value: -0.1

Lower boundry for inner window when do cumulant calculation

Mandatory for calc_mode:  spectral-sespectral-cum

JOB CONTROL
cum_inner_emax

Variable type: real
Default value: 0.1
Typical value: 0.1

Upper boundry for inner window when do cumulant calculation

Mandatory for calc_mode:  spectral-sespectral-cum

JOB CONTROL
cum_outer_emin

Variable type: real
Default value: -0.1
Typical value: -0.1

Lower boundry for outer window when do cumulant calculation

Mandatory for calc_mode:  spectral-sespectral-cum

JOB CONTROL
cum_outer_emax

Variable type: real
Default value: 0.1
Typical value: 0.1

Upper boundry for outer window when do cumulant calculation

Mandatory for calc_mode:  spectral-sespectral-cum

JOB CONTROL
cum_de

Variable type: real
Default value: 0.1
Typical value: 0.1

Energy difference. Decide output spectral energy step with spectral : cum_de / spectral_np

Optional for calc_mode:  spectral-sespectral-cum

JOB CONTROL
spectral_emin

Variable type: real
Default value: -0.1
Typical value: -0.1

lower boundry of spectral energy scale

Mandatory for calc_mode:  spectral-cum

JOB CONTROL
spectral_emax

Variable type: real
Default value: 0.1
Typical value: 0.1

upper boundry of spectral energy scale

Mandatory for calc_mode:  spectral-cum

JOB CONTROL
cum_outer_np

Variable type: integer
Default value: 1
Typical value: 1

Number of points of outer window

Optional for calc_mode:  spectral-sespectral-cum

JOB CONTROL
spectral_np

Variable type: integer
Default value: 1
Typical value: 1

Number of points in the energy window. Decide output spectral energy step with cum_de: cum_de / spectral_np




 Boltzmann Transport Equation

BOLTZMANN TRANSPORT EQUATION
boltz_kdim

Variable type: integer
Default value: (1,1,1)
Dimensions: (3)

Number of k points along each dimension for the Boltzmann equation.

Mandatory for calc_mode:  setuptrans-rtatrans-itatrans-mag-rtatrans-mag-itatrans-ppdynamics-run
dynamics-ppspinlifetime

BOLTZMANN TRANSPORT EQUATION
boltz_qdim

Variable type: integer
Default value: (0,0,0)
Dimensions: (3)

Number of q points along each dimension for the Boltzmann equation. It should be same with boltz_kdim.

Optional for calc_mode:  trans-rtatrans-itatrans-mag-rtatrans-mag-itadynamics-run

BOLTZMANN TRANSPORT EQUATION
band_min

Variable type: integer
Default value: 1

Lowest band included.

Mandatory for calc_mode:  spinlifetimespectral-respectral-sespectral-cum Optional for calc_mode:  ephmatsetupimsigmameanfptrans-rtatrans-itatrans-mag-rta
trans-mag-itatrans-ppdynamics-rundynamics-ppephmat_spinimsigma_spin

BOLTZMANN TRANSPORT EQUATION
band_max

Variable type: integer
Default value: 9999999

Highest band included.

Mandatory for calc_mode:  spinlifetimespectral-respectral-sespectral-cum Optional for calc_mode:  ephmatsetupimsigmameanfptrans-rtatrans-itatrans-mag-rta
trans-mag-itatrans-ppdynamics-rundynamics-ppephmat_spinimsigma_spin

BOLTZMANN TRANSPORT EQUATION
boltz_emin

Variable type: real
Default value: -9999.0
Units: eV

Bottom of the energy window for the Boltzmann equation.

Optional for calc_mode:  setuptrans-rtatrans-itatrans-mag-rtatrans-mag-itatrans-ppdynamics-run
dynamics-ppspinlifetime

BOLTZMANN TRANSPORT EQUATION
boltz_emax

Variable type: real
Default value: 9999.0
Units: eV

Top of the energy window for the Boltzmann equation.

Optional for calc_mode:  setuptrans-rtatrans-itatrans-mag-rtatrans-mag-itatrans-ppdynamics-run
dynamics-ppspinlifetime

BOLTZMANN TRANSPORT EQUATION
boltz_nstep

Variable type: integer
Default value: 0
Typical value: 50

Number of iterations for solving the Boltzmann transport equation.

Mandatory for calc_mode:  trans-itatrans-mag-rtatrans-mag-itadynamics-run Optional for calc_mode:  trans-rta

BOLTZMANN TRANSPORT EQUATION
boltz_de

Variable type: real
Default value: 1.0
Units: meV

Energy step for the integrals in the Boltzmann equation.

Optional for calc_mode:  trans-rtatrans-itatrans-mag-rtatrans-mag-itatrans-ppdynamics-pp

BOLTZMANN TRANSPORT EQUATION
delta_smear

Variable type: real
Default value: 10.0
Units: meV

Smearing for the Dirac delta function.

Optional for calc_mode:  imsigmatrans-rtatrans-itatrans-mag-rtatrans-mag-itadynamics-runimsigma_spin
spectral-respectral-sespectral-cum

BOLTZMANN TRANSPORT EQUATION
full_ite

Variable type: logical
Default value: .true.

Solve BTE with both E- and T-fields iteratively.

Optional for calc_mode:  trans-rtatrans-itatrans-mag-rtatrans-mag-ita

BOLTZMANN TRANSPORT EQUATION
phfreq_cutoff

Variable type: real
Default value: 1.0
Typical value: 1.0
Units: meV

Phonon energy threshold. Phonons with energy smaller than phfreq_cutoff will be excluded.

Optional for calc_mode:  ephmatimsigmatrans-rtatrans-itatrans-mag-rtatrans-mag-itadynamics-run
ephmat_spinimsigma_spinspectral-respectral-sespectral-cum

BOLTZMANN TRANSPORT EQUATION
trans_thr

Variable type: real
Default value: 0.002

Threshold for the iterative procedure.

Optional for calc_mode:  trans-rtatrans-itatrans-mag-rtatrans-mag-ita

BOLTZMANN TRANSPORT EQUATION
lmagsym

Variable type: logical
Default value: .false.

Use magnetic symmetry with axis along z.




 Polar correction

POLAR CORRECTION
polar_split

Variable type: string
Default value: ''
Options:
''
(leave blank) both the polar and nonpolar parts, no split
'polar'
polar part only
'rmpol'
remainder part

Polar correction mode.

Optional for calc_mode:  imsigmaimsigma_spin



 Ultra-fast dynamics

ULTRA-FAST DYNAMICS
time_step

Variable type: real
Default value: 1.0
Typical value: 1.0
Units: fs

Time step for the carrier dynamics.

Mandatory for calc_mode:  dynamics-run

ULTRA-FAST DYNAMICS
output_nstep

Variable type: integer
Default value: 1

Print out the results every output_nstep time steps.

Mandatory for calc_mode:  dynamics-run

ULTRA-FAST DYNAMICS
boltz_init_dist

Variable type: string
Default value: ''
Typical value: 'gaussian'
Options:
'restart'
restart from the previous calculations.
'lorentz'
Loretnzian distribution
'fermi'
Fermi-Dirac distribution
'gaussian'
Gaussian distribution

Initial electron distribution at time zero.

Mandatory for calc_mode:  dynamics-run

ULTRA-FAST DYNAMICS
pump_pulse

Variable type: logical
Default value: .false.
Typical value: .false.

During first pulse_num_steps of dynamics, add extra occupations (hence the name, 'pump'). A file pump_pulse.h5 (or specified name with pump_pulse_fname) must be present in the directory and it must have the following datasets: pulse_num_steps, pulse_snap_t_[...], time_step.

Optional for calc_mode:  dynamics-run

ULTRA-FAST DYNAMICS
pump_pulse_fname

Variable type: string
Default value: 'pump_pulse.h5'

Name of HDF5 file containing pump pulse information, refer to pump_pulse.

Optional for calc_mode:  dynamics-run

ULTRA-FAST DYNAMICS
boltz_init_e0

Variable type: real
Default value: -9999.0
Typical value: 1.0
Units: eV

Energy parameter used to generate initial distribution. Needs to be specified for boltz_init_dist='lorentz' (center), 'gaussian' (center), or 'fermi' (chemical potential).

Optional for calc_mode:  dynamics-run

ULTRA-FAST DYNAMICS
boltz_init_smear

Variable type: real
Default value: 20.0
Typical value: 1.0
Units: meV

The broadening or width of the initial distribution for boltz_init_dist='lorentz' or 'gaussian', or temperature (in meV) for 'fermi'.

Optional for calc_mode:  dynamics-run

ULTRA-FAST DYNAMICS
boltz_init_ampl

Variable type: real
Default value: 0.1
Typical value: 0.1
Units: arbitrary

Amplitude of the initial distribution, currently implemented for boltz_init_dist='gaussian'.

Optional for calc_mode:  dynamics-run

ULTRA-FAST DYNAMICS
solver

Variable type: string
Default value: 'rk4'
Options:
'euler'
Euler method (first-order)
'rk4'
fourth-order Runge-Kutta method

Solver type for the Boltzmann transport equation.

Optional for calc_mode:  dynamics-run

ULTRA-FAST DYNAMICS
boltz_efield

Variable type: real
Default value: (0.0, 0.0, 0.0)
Dimensions: (3)
Units: V/cm

External electric field for the ultrafast real-time dynamics.

Optional for calc_mode:  dynamics-rundynamics-pp

ULTRA-FAST DYNAMICS
boltz_norm_dist

Variable type: logical
Default value: .false.

Normalize the distribution function at each step of the real-time Boltzmann Transport Equation simulation.

Optional for calc_mode:  dynamics-run

ULTRA-FAST DYNAMICS
boltz_acc_thr

Variable type: real
Default value: 0.0
Typical value: 1.0
Units: cm/s2

Drift acceleration threshold. If specified, the real-time simulation will stop when at least 10 last iterations had a drift acceleration lower than the threshold.

Optional for calc_mode:  dynamics-run

ULTRA-FAST DYNAMICS
boltz_nstep_min

Variable type: integer
Default value: 1
Typical value: 100

Minimum number of iterations of real-time dynamics. Used with boltz_acc_thr. Will not be applied if boltz_acc_thr is not specified, in this case, use boltz_nstep instead.

Optional for calc_mode:  dynamics-run

ULTRA-FAST DYNAMICS
ph_mode_exclude_ranges

Variable type: integer
Default value: -1
Dimensions: (30, 2)

Specify the phonon modes to exclude in rt-BTE: the g2 factor will be set to zero for these modes. The modes to exclude are specified with ranges, e.g. to exclude modes from 2 to 5, write in the input file: ph_mode_exclude_ranges(1,:)=2,5 . To add another range to exclude, specify ph_mode_exclude_ranges(2,:)=7,10 . One can use up to 30 ranges. To exclude one mode, set first and second number to the phonon mode index. The indices of the phonon modes to exclude must be smaller than the total number of modes in a system.

Optional for calc_mode:  dynamics-run