Two-photon absorption

Two-photon absorption#

In experimental studies of two- and multi-photon absorption processes, the multiphoton transition strength, a function of all the frequencies of all photons absorbed, is analyzed. For two laser sources of circular frequency \(\omega_1\) and \(\omega_2\) with associated wave lengths \(\lambda_1\) and \(\lambda_2\), the two-photon transition strength \(\delta (\omega_1,\omega_2)\) for the transition between states \(\mid 0 \rangle \) and \( \mid n \rangle \) in isotropic samples is given by:

\[\begin{eqnarray*} \delta^{\rm TPA}_{0n}(\omega_1,\omega_2) & = & {\sf F} {\mathcal{S}}_{\alpha\alpha}^{0n}(\omega_1,\omega_2) {\mathcal{S}}_{\beta\beta}^{on, *}(\omega_1,\omega_2) + {\sf G} {\mathcal{S}}_{\alpha\beta}^{0n}(\omega_1,\omega_2) {\mathcal{S}}_{\alpha\beta}^{0n, *}(\omega_1,\omega_2)+ \\ & & + {\sf H} {\mathcal{S}}_{\alpha\beta}^{0n}(\omega_1,\omega_2) {\mathcal{S}}_{\beta\alpha}^{0n, *}(\omega_1,\omega_2), \end{eqnarray*}\]

where \({\sf F}\), \({\sf G}\) and \({\sf H}\) are numbers depending on the polarization state of the two photons and on the geometrical set-up (mutual direction of the laser beams) and, in the dipole approximation,

\[ {\mathcal{S}}_{\alpha\beta}^{0n}(\omega_1,\omega_2) = \frac{1}{\hbar} \sum_{m} \Big \{ \frac{(\hat{\mu}_{\alpha})_{0m} (\hat{\mu}_{\beta})_{mn} } {\omega_{m0} - \omega_1} + \frac{(\hat{\mu}_{\beta})_{0m} (\hat{\mu}_{\alpha})_{mn} } {\omega_{m0} - \omega_2} \Big \} \]

is the second-rank, two-photon tensor. In the equation for \(\delta^{\rm TPA}_{0n}(\omega_1,\omega_2)\) the summation runs over the whole set of excited states, the energy conservation relation \(\omega_1+\omega_2=\omega_{n0}\) applies, and off-resonance conditions are implied - that is, the frequencies \(\omega_1\) and \(\omega_2\) are sufficiently far off the values at which the denominators vanish. The tensor is non-symmetric in the exchange of the two frequencies except for \(\omega_1=\omega_2\). For the special case of a one-colour beam - that is, a monochromatic light source the transition matrix is symmetric and (using \(\omega=\omega_1=\omega_2\))

\[ \delta^{\rm TPA}_{0n}(\omega) = {\mathcal{S}}_{\lambda\lambda}^{0n}(\omega) {\mathcal{S}}_{\mu\mu}^{0n, *}(\omega) + 2 {\mathcal{S}}_{\lambda\mu}^{0n}(\omega) {\mathcal{S}}_{\lambda\mu}^{0n, *}(\omega) \]

for linear polarization of the incident light and

\[ \delta^{\rm TPA}_{0n}(\omega) = - {\mathcal{S}}_{\lambda\lambda}^{0n}(\omega) {\mathcal{S}}_{\mu\mu}^{0n, *}(\omega) + 3 {\mathcal{S}}_{\lambda\mu}^{0n}(\omega) {\mathcal{S}}_{\lambda\mu}^{0n, *}(\omega) \]

for circular polarization, respectively.

It can be shown that the two-photon absorption transition rate (cross section) can be obtained from the single residue of the cubic response function.

Two-photon absorption transition amplitudes \({\mathcal{S}}_{\alpha\beta}^{0n}(\omega_1,\omega_2)\), can also be extracted from the single residue of a quadratic response function.

The following relations hold for the matrix elements of the two-photon transition tensor

\[\begin{eqnarray*} S_{\alpha\beta}^{0n}(\omega) & = & S_{\beta\alpha}^{0n}(\omega_{n0} - \omega), \\ S_{\alpha\beta}^{0n}(\omega) & = & S_{\alpha\beta}^{n0}(- \omega)^\ast. \end{eqnarray*}\]

For methods which do not fulfill this last equation (such as the coupled cluster approach), one can instead use the two-photon transition strength

\[ {F}^{0n}_{\alpha\beta,\gamma\delta}(\omega) = \frac{1}{2} \{ S_{\alpha\beta}^{0n}(-\omega) S_{\gamma\delta}^{n0}(\omega) +S_{\gamma\delta}^{0n, *}(-\omega) S_{\alpha\beta}^{0n, *}(\omega)\}, \]

where \(S_{\alpha\beta}^{0n}(\omega)\) and \(S_{\alpha\beta}^{n0}(\omega)\) are the left and right transition moments.

Three-photon absorption has been described through the single residue of the cubic response function. Similarly to two-photon absorption, one can discuss the third-order (left and right) transition moments and strengths and define the three-photon transition strengths in this manner for nonvariational wave functions.

In general, the computational requirements for the calculation of molecular properties from the residues of the response functions inherit the requirements from the response functions themselves; that is, the selection of the basis sets has to be done considering the operators appearing in the expression for the transition moments and excited-state properties. However, as the residues are connected to specific excited states, the nature of the probed excited state also needs to considered