# Difference between revisions of "Introduction to Real Time propagation in Yambo"

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 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 $| v_{i \mathbf{k}} \rangle$ are the time-dependent valence bands, Hsysk is the Hamiltonian, $\mathcal E(t)$ is the external field and the term $i e \partial k$ corresponds to the dipole operator in periodic systems. For more details on this last term see Ref. [1] and [2]. 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. [3]. 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. We calculate the Berry polarization

### Linear Response

In general the induced polarization is the sum of the responses of the system to all orders in the external field: $\bf{P}(t) = \sum_{n=-\infty}^{+\infty} \bf{P}_n (t)$ 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

$\chi(\omega) = \frac{P(\omega)}{E(\omega)}$

that is related to the dielectric constant through the relation $\epsilon(\omega) = 1 + 4 \pi \chi(\omega)$. See also Ref. [2] More details will be given during the tutorial

### 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 $\bf{P}(t) = \sum_{n=-\infty}^{+\infty} \bf{p}_n e^{-i\omega_n t}$ where the coefficient $\bf{p}_1,...,\bf{p}_n$ are related to $\chi^{(1)},...,\chi^{(n)}$. In practice we use the scheme described in Ref. [4] and show below to extract the non-linear coefficients:

More details will be given during the tutorials Real time approach to non-linear response and Correlation effects in the non-linear response

## References

1. PRB 69, 085106 (2004)
2. C. Attaccalite and M. Gruning Rev. B, 88, 235113 (2013)
3. Non-linear response in extended systems: a real-time approach
4. C. Attaccalite and M. Gruning, Rev. B, 88, 235113 (2013)