Computing linear optical spectra in the presence of nonadiabatic effects on Graphics Processing Units using molecular dynamics and tensor-network approaches

Evan Lambertson,Dayana Bashirova,Kye E. Hunter,Benhardt Hansen,Tim J. Zuehlsdorff

2024-08-13

Abstract:We compare two recently developed strategies, implemented in open source software packages, for computing linear optical spectra in condensed phase environments in the presence of nonadiabatic effects. Both approaches rely on computing excitation energy and transition dipole fluctuations along molecular dynamics (MD) trajectories, treating molecular and environmental degrees of freedom on the same footing. Spectra are then generated in two ways: In the recently developed Gaussian Non-Condon Theory (GNCT), the linear response functions are computed in terms of independent adiabatic excited states, with non-Condon effects described through spectral densities of transition dipole fluctuations. For strongly coupled excited states, we instead parameterize a linear vibronic coupling (LVC) Hamiltonian directly from spectral densities of energy fluctuations and diabatic couplings computed along the MD trajectory. The optical spectrum is then calculated using powerful, numerically exact tensor-network approaches. Both the electronic structure calculations to sample system fluctuations and the quantum dynamics simulations using tensor-network methods are carried out on graphics processing units (GPUs), enabling rapid calculations on complex condensed phase systems. We assess the performance of the approaches using model systems in the presence of a conical intersection (CI), and the pyrazine molecule in different solvent environments.

Chemical Physics

What problem does this paper attempt to address?

The paper is primarily dedicated to addressing the efficient and accurate calculation of linear optical spectra in complex condensed phase environments in the presence of non-adiabatic effects. Specifically, the study focuses on the following key points:

1. **Non-Adiabatic Effects**: When the coupling between electronic states within a molecule becomes significant, the traditional adiabatic approximation (i.e., the Born-Oppenheimer approximation) is no longer applicable. This non-adiabatic effect is particularly pronounced near conical intersections (CIs), posing challenges for spectral simulation.

2. **Environmental Coupling**: The interaction between a molecule and its surrounding environment (such as a solvent) has a significant impact on spectral properties. Traditional methods often handle these effects through simplified means such as Polarizable Continuum Models (PCMs), but this approach may not fully capture the complex solvent-solute dynamics.

To address the above challenges, the paper proposes and compares two strategies, both of which combine Molecular Dynamics (MD) sampling with tensor network methods to efficiently calculate linear optical spectra. Specifically:

- **Gaussian Non-Adiabatic Coupling Theory (GNCT)**: Suitable for cases of weak non-adiabatic interactions, capable of describing weak coupling effects between different electronic excited states, such as intensity borrowing. This method is based on the framework of open quantum systems, retaining the full coupling between electronic states and environmental degrees of freedom without the need for assuming a frozen solvent environment.

- **Tensor Network Methods**: Targeting strong coupling scenarios, particularly when there is significant interaction between electronic excited states. By directly extracting energy fluctuations and dihedral couplings from MD trajectories, parameterizing the Linear Vibronic Coupling (LVC) Hamiltonian, and utilizing powerful tensor network techniques to compute optical spectra.

Both methods implement electronic structure calculations and quantum dynamics simulations on Graphics Processing Units (GPUs), greatly enhancing the computational speed for complex condensed phase systems. Additionally, the effectiveness and advantages of these methods are validated through tests on model systems and pyrazine molecules in different solvent environments.