Introduction to Real Time propagation in Yambo

Introduction

A new feature in yambo is the numerical integration of a time-dependent (TD) equation of motion (EOM), able to describe the evolution of the electronic system under the action of an external laser pulse. Similarly to the equilibrium case, the most diffuse ab initio approaches to real-time propagation are based on TD-DFT and there exist a number of GPL codes available to this end. On the contrary the implementation of real time propagation within MBPT is an almost unique feature of the yambo code. Two different schemes are available. In one case, the density matrix of the system, is propagated in time, as described in the first section. In the second case an effective time dependent Schrodinger equation for the occupied orbitals is propagated, as described in the second section.

Time Dependent Equation for the Reduced One--Body Density--Matrix

The first scheme is a Time Dependent equation for the one body reduced density matrix. The EOM for the density matrix projected in the space of the single particle wave-functions is derived from non-equilibrium (NEQ) many-body perturbation theory and reads

Here we underline quantities which are vectors in the trans­ition space. The Hamiltonian contains the equilibrium eigenvalues plus the variation of the self-energy. More details can be found in Ref.

represents the external potential written in the length gauge; shape, polarization, intensity (and eventually frequency) of the field E can be selected in input. r is the position operator. The coupling to the external field is exact up to first order.

Linear Response

From the knowledge of the density matrix, the first order polarization is computed at each time step. The spectrum of the system can then be obtained by the Fourier transform of the polarization, which can be done as a post-processing step. Absorption is thus obtained via the dipole-dipole response function (equivalent to the length gauge in linear response). [math]\displaystyle{ \chi(\omega) = \frac{P(\omega)}{E(\omega)} }[/math]

More details will be given during the tutorials:

• Real time approach to linear response
• Real time Bethe-Salpeter Equation (TDSE)

Time Dependent Effective Schrödinger Equation

The second scheme is a Time Dependent Effective Schrödinger Equation(TD-ESE), where the coupling between electrons and the external field is described by means of the Modern Theory of Polarization:

Where [math]\displaystyle{ | v_{i \mathbf{k}} \rangle }[/math] are the time-dependent valence bands, Hsysk is the Hamiltonian, [math]\displaystyle{ \mathcal E(t) }[/math] is the external field and the term [math]\displaystyle{ i e \partial k }[/math] corresponds to the dipole operator in periodic systems. For more details on this last term see Ref. and . Differently from other codes, the equations of motion of Yambo are in the length gauge and not in the velocity one. If you want to know more about the advantages and disadvantages of the two gauges read section 2.7 of Ref. . The Hamiltonian and the initial wave-functions are obtained from DFT. Then in order to obtain linear response, we probe the system with a delta function field, that excites all the frequencies at the same footing, and Fourier transforms. The Hamiltonian is obtained from Green's function theory and is written as:

where [math]\displaystyle{ \rho_0 }[/math] and [math]\displaystyle{ \gamma_0 }[/math] are the density and the single-particle density matrix at equilibrium, [math]\displaystyle{ \Sigma^{COHSEX} }[/math] is the COHSEX self-energy and [math]\displaystyle{ V^H }[/math] is the Hartree potential. This Hamiltonian corresponds to a real-time version of the Bethe-Salpeter Equation. Then the polarization is calculated by means of Berry's phase

Linear Response

In general the induced polarization is the sum of the responses of the system to all orders in the external field: [math]\displaystyle{ \bf{P}(t) = \sum_{n=-\infty}^{+\infty} \bf{P}_n (t) }[/math] If the perturbation is weak enought the response is dominated by the first order term and from its Fourier transform we get the response function. Indeed

[math]\displaystyle{ \chi(\omega) = \frac{P(\omega)}{E(\omega)} }[/math]

that is related to the dielectric constant through the relation [math]\displaystyle{ \epsilon(\omega) = 1 + 4 \pi \chi(\omega) }[/math]. See also Ref.

More details will be given during the tutorials:

• Real time approach to linear response
• Real time Bethe-Salpeter Equation (TDSE)

Non Linear Response

Non-linear response calculations are similar to the linear response case, with the only difference that we excite the system with a sinusoidal external field. A perturbation with a specific frequency allows to expand the polarization in the form [math]\displaystyle{ \bf{P}(t) = \sum_{n=-\infty}^{+\infty} \bf{p}_n e^{-i\omega_n t} }[/math] where the coefficient [math]\displaystyle{ \bf{p}_1,...,\bf{p}_n }[/math] are related to [math]\displaystyle{ \chi^{(1)},...,\chi^{(n)} }[/math]. In practice we use the scheme described in Ref. and show below to extract the non-linear coefficients:

More details will be given during the tutorials:

• Real time approach to non-linear response
• Correlation effects in the non-linear response

References

Prev: ICTP Tutorials

Now: ICTP Tutorials --> Linear Response --> Theory

Next: Prerequisites