Abstract:Theoretical prediction of electronic absorption spectra without input from experiments is no easy feat, as it requires addressing all of the factors that affect line shapes. In practice, however, the methodologies are limited to treat these ingredients only to a certain extent. Here, we present a multiscale protocol that addresses the temperature, solvent, and nuclear quantum effects as well as anharmonicity and the reconstruction of the final spectra from individual transitions. First, quantum mechanics/molecular mechanics (QM/MM) molecular dynamics is conducted to obtain trajectories of solute–solvent configurations, from which the corresponding quantum-corrected ensembles are generated through the generalized smoothed trajectory analysis (GSTA). The optical spectra of the ensembles are then produced by calculating vertical transitions using time-dependent density-functional theory (TDDFT) with implicit solvation. To obtain the final spectral shapes, the stick spectra from TDDFT are convoluted with Gaussian kernels where the half-widths are determined by a statistically motivated strategy. We have tested our method by calculating the UV–vis spectra of a recently discovered acridine photocatalyst in two redox states. Vibronic progressions and broadenings due to the finite lifetime of the excited states are not included in the methodology yet. Nuclear quantization affects the relative peak intensities and widths, which is necessary to reproduce the experimental spectrum. We have also found that using only the optimized geometry of each molecule works surprisingly well if a proper empirical broadening factor is applied. This is explained by the rigidity of the conjugated chromophore moieties of the selected molecules, which are mainly responsible for the excitations in the spectra. In contrast, we have also shown that other parts of the molecules are flexible enough to feature anharmonicities that impair the use of other techniques such as Wigner sampling.This article has not yet been cited by other publications.

On-the-fly ab initio semiclassical evaluation of third-order response functions for two-dimensional electronic spectroscopy

On-the-fly ab initio semiclassical evaluation of vibronic spectra at finite temperature

On-the-fly ab initio semiclassical evaluation of absorption spectra of polyatomic molecules beyond the Condon approximation

On-the-fly ab initio three thawed Gaussians approximation: A semiclassical approach to Herzberg-Teller spectra

Applicability of the thawed Gaussian wavepacket dynamics to the calculation of vibronic spectra of molecules with double-well potential energy surfaces

Two-dimensional electronic spectroscopy from first principles

Multiscale Modeling of Electronic Spectra Including Nuclear Quantum Effects

Two-dimensional Electronic Spectra from the Hierarchical Equations of Motion Method: Application to Model Dimers.

On-the-Fly Ab Initio Hagedorn Wavepacket Dynamics: Single Vibronic Level Fluorescence Spectra of Difluorocarbene

Exploring the Influence of Approximations for Simulating Valence Excited X-ray Spectra

Optical spectra in the condensed phase: Capturing anharmonic and vibronic features using dynamic and static approaches

Toward the improvement of vibronic spectra and non-radiative rate constants using the vertical Hessian method

The {\it ab initio} calculation of spectra of open shell diatomic molecules

Accelerating Anharmonic Spectroscopy Simulations via Local-Mode, Multilevel Methods.

Ab initio simulation of single vibronic level fluorescence spectra of anthracene using Hagedorn wavepackets

Which Electronic Structure Method to Choose in Trajectory Surface Hopping Dynamics Simulations? Azomethane as a Case Study

The symmetric quasi-classical model using on-the-fly time-dependent density functional theory within the Tamm-Dancoff approximation

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