Excited states dynamics in time-dependent density functional theory: high-field molecular dissociation and harmonic generation

Alberto Castro,M. A. L. Marques,Julio A. Alonso,George F. Bertsch,Angel Rubio

DOI: https://doi.org/10.48550/arXiv.cond-mat/0206307

2002-06-17

Abstract:We present a theoretical description of femtosecond laser induced dynamics of the hydrogen molecule and of singly ionised sodium dimers, based on a real-space, real-time, implementation of time-dependent density functional theory (TDDFT). High harmonic generation, Coulomb explosion and laser induced photo-dissociation are observed. The scheme also describes non-adiabatic effects, such as the appearance of even harmonics for homopolar but isotopically asymmetric dimers, even if the ions were treated classically. This TDDFT-based method is reliable, scalable, and extensible to other phenomena such as photoisomerization, molecular transport and chemical reactivity.

Condensed Matter

What problem does this paper attempt to address?

The problem that this paper attempts to solve is: how to use time - dependent density functional theory (TDDFT) to describe and simulate the excited - state dynamic behavior of molecules under the action of a strong laser field. Specifically, the research focuses on:

1. **High - Harmonic Generation (HHG)**:
- When a molecule is irradiated by a strong laser pulse, high - order harmonics will be generated. These harmonics can be used as soft X - ray sources and have important application prospects.

2. **Coulomb Explosion**:
- Under the action of a strong laser field, the molecule may dissociate, causing the Coulomb repulsion between charged particles to make the molecular fragments fly apart at high speed.

3. **Laser - Induced Photodissociation**:
- A laser pulse can trigger the photodissociation of a molecule, that is, after the molecule absorbs a photon, the chemical bond breaks, and finally the molecule decomposes.

To achieve these goals, the author proposes a method based on real - space, real - time time - dependent density functional theory (TDDFT), which combines the quantum - mechanical propagation of the electronic wave function and the classical treatment of nuclear motion. This method can effectively simulate the dynamic processes of multi - electron systems in a strong laser field and can be extended to larger systems, such as photoisomerization, molecular transport, and chemical reactivity problems.

### Formula Explanation

The formulas involved in the paper mainly include the Lagrangian and equations of motion:

1. **Lagrangian**:
\[
L=\sum_{\alpha}\left[\frac{1}{2m_{\alpha}}\dot{\vec{R}}_{\alpha}^{2}+Z_{\alpha}\vec{R}_{\alpha}\cdot\vec{E}(t)\right]-\sum_{\alpha < \alpha'}\frac{Z_{\alpha}Z_{\alpha'}}{\vert\vec{R}_{\alpha}-\vec{R}_{\alpha'}\vert}+\sum_{i}\langle\phi_{i}\vert i\hbar\frac{\partial}{\partial t}-e\vec{x}\cdot\vec{E}(t)\vert\phi_{i}\rangle - E_{DFT}(\{\phi\}, \{\vec{R}\})
\]
where \(E_{DFT}\) is the Kohn - Sham density functional, which depends on the electronic orbitals \(\{\phi\}\) and nuclear coordinates \(\{\vec{R}\}\), and \(\vec{E}(t)\) is the time - varying electric field.

2. **Newton's Equation of Motion**:
\[
m_{\alpha}\frac{d^{2}\vec{R}_{\alpha}}{dt^{2}}=-\nabla_{\vec{R}_{\alpha}}E_{DFT}-\sum_{\alpha'\neq\alpha}\frac{(\vec{R}_{\alpha}-\vec{R}_{\alpha'})Z_{\alpha}Z_{\alpha'}}{\vert\vec{R}_{\alpha}-\vec{R}_{\alpha'}\vert^{3}}+Z_{\alpha}\vec{E}(t)
\]

3. **Generalized Hellmann - Feynman Theorem**:
\[
-\nabla_{\vec{R}_{\alpha}}E_{DFT}=-\sum_{i}\langle\phi_{i}\vert\nabla_{\alpha}H_{KS}\vert\phi_{i}\rangle
\]
where \(H_{KS}=\frac{\delta E_{DFT}}{\delta n(\vec{r}, t)}\) is the Kohn - Sham Hamiltonian.