Pump and Probe

This tutorial works only with Yambo version > 5.2, that will be released soon.

In this tutorial we will show you how to setup two external fields in yambo_nl, to perform pump and probe simulation. At present this tutorial explain only how to run simulation and obtain the final polarization. We are working on python script to analyze results that will be make available soon.

We start from the AlAs wave-functions of the Real time approach to non-linear response tutorial and remove symmetries according to the Pump and Probe external fields as explained in the Prerequisites for Real Time propagation with Yambo.

In our example we choose direction [1,1,0] for the pump and [0,0,1] for the probe, the corresponding ypp.in file will be:

fixsyms                          # [R] Remove symmetries not consistent with an external perturbation
% Efield1
 1.000000 | 1.000000 | 0.000000 |        # First external Electric Field
%
% Efield2
 0.000000 | 0.000000 | 1.000000 |        # Additional external Electric Field
%
BField= 0.000000           T     # [MAG] Magnetic field modulus
Bpsi= 0.000000             deg   # [MAG] Magnetic field psi angle [degree]
Btheta= 0.000000           deg   # [MAG] Magnetic field theta angle [degree]
#RmAllSymm                     # Remove all symmetries 
RmTimeRev                     # Remove Time Reversal
#RmSpaceInv                    # Remove Spatial Inversion

then we generate the Pump&Probe input with the command yambo_nl -u p -V rt:

nloptics                         # [R] Non-linear spectroscopy
DIP_Threads=0                    # [OPENMP/X] Number of threads for dipoles
NL_Threads=0                     # [OPENMP/NL] Number of threads for nl-optics
% NLBands
   3 |  6 |                           # [NL] Bands range
%
NLverbosity= "high"              # [NL] Verbosity level (low | high)
NLstep= 0.010000           fs    # [NL] Time step length
 NLtime= 75.0000         fs      # [NL] Simulation Time
NLintegrator= "INVINT"           # [NL] Integrator ("EULEREXP/RK2/RK4/RK2EXP/HEUN/INVINT/CRANKNIC")
NLCorrelation= "IPA"             # [NL] Correlation ("IPA/HARTREE/TDDFT/LRC/LRW/JGM/SEX")
NLLrcAlpha= 0.000000             # [NL] Long Range Correction
 NLDamping= 0.100000        eV      # [NL] Damping (or dephasing)
 RADLifeTime= 15.00000      fs      # [RT] Radiative life-time (if negative Yambo sets it equal to Phase_LifeTime in NL)
#UseDipoles                    # [NL] Use Covariant Dipoles (just for test purpose)
#FrSndOrd                      # [NL] Force second order in Covariant Dipoles
#EvalCurrent                   # [NL] Evaluate the current
HARRLvcs= 10417            RL    # [HA] Hartree     RL components
EXXRLvcs= 10417            RL    # [XX] Exchange    RL components
% ExtF_Dir
  1.000000 | 1.000000 | 0.000000 |        # [NL ExtF] Field Versor
%
% ExtF_Freq
  1.500000 | 1.500000 |         eV    # [NL ExtF] Field Frequency
%
 ExtF_Int=  1000.00         kWLm2  # [NL ExtF] Field Intensity
ExtF_Width= 5.00000       fs    # [NL ExtF] Field Width
 ExtF_kind= "QSSIN"               # [NL ExtF] Kind(SIN|SOFTSIN|RES|ANTIRES|GAUSS|DELTA|QSSIN)
ExtF_Tstart= 0.010000      fs    # [NL ExtF] Initial Time 
% ExtF2_Dir
  0.000000 | 0.000000 | 1.000000 |        # [NL ExtF] Field Versor
%
% ExtF2_Freq
 0.100000 | 0.100000 |         eV    # [NL ExtF] Field Frequency
%
 ExtF2_Int= 10.0000        kWLm2 # [NL ExtF] Field Intensity
ExtF2_Width= 0.000000      fs    # [NL ExtF] Field Width
 ExtF2_kind= "DELTA"              # [NL ExtF] Kind(SIN|SOFTSIN|RES|ANTIRES|GAUSS|DELTA|QSSIN)
 ExtF2_Tstart= 30.00000     fs    # [NL ExtF] Initial Time