Excited state absorption of DNA bases in the gas phase and in chloroform solution: a comparative quantum mechanical study

Excited-State Absorption: Reference Oscillator Strengths, Wave Function, and TDDFT Benchmarks

Excited-state absorption (ESA) corresponds to the transition between two electronic excited states and is a fundamental process for probing and understanding light-matter interactions. Accurate modeling of ESA is indeed often required to interpret time-resolved experiments. In this contribution, we present a dataset of 53 ESA oscillator strengths in three different gauges and the associated vertical transition energies between 71 excited states of 21 small- and medium-sized molecules from the QUEST database. In a few cases, we additionally investigated the effect of geometry relaxation on excited-state geometries. The reference values were obtained within the quadratic response (QR) CC3 formalism using eight different Dunning basis sets. We found that the d-aug-cc-pVTZ basis set is always adequate while its more compact double-ζ counterpart, d-aug-cc-pVDZ, performs well in most cases. These QR-CC3 data allow us to assess the performance of QR-TDDFT, with and without applying the Tamm-Dancoff approximation, using a panel of global and range-separated hybrids (B3LYP, BH&HLYP, CAM-B3LYP, LC-BLYP33, and LC-BLYP47), as well as several lower-order wave function methods, i.e., QR-CCSD, QR-CC2, EOM-CCSD, ISR-ADC(2), and ISR-ADC(3). We show that QR-TDDFT delivers acceptable errors for ESA oscillator strengths with CAM-B3LYP showing particular promise, especially for the largest molecules of our set, and in the Franck-Condon (FC) region. We also find that ISR-ADC(3) exhibits excellent performance in this region. When using excited-state optimal geometries, the relative performance of wave function-based approaches remains consistent with trends observed in the Franck-Condon region. However, for TD(A)-DFT, the accuracy varies more significantly, as the performance of different exchange-correlation functionals significantly depends on the chosen geometry.

Theoretical aspects of the adsorption of normal and modified base pairs of DNA on graphene models toward DNA sequencing

A theoretical understanding of the adsorption of DNA base pairs (GC, AT, CAF-T and CAF-C) on the graphene models (Gr, SiGr and SiGr-COOH) is investigated. Among the complexes, SiGr-COOH_AT is found to have the highest adsorption energies of -202.83 kcal/mol. The strong adsorption between DNA base pairs and the SiGr-COOH model leads to concomitant charge transfer responsible for the stability of the corresponding models and is verified with NBO analysis. AIM analysis discloses the high orbital overlap that signifies the strong interaction. Closed-shell interactions are observed through the positive values of total electron density, and it is also observed that Si-O(N) interaction has both covalent and electrostatic characteristics. This is the first theoretical attempt to investigate the adsorption of DNA base pairs on SiGr-COOH, which is more favourable than other models and may call for further experimental studies, which is crucial in developing new bio-sensors.Communicated by Ramaswamy H. Sarma.

Assessment of State-Averaged Driven Similarity Renormalization Group on Vertical Excitation Energies: Optimal Flow Parameters and Applications to Nucleobases

We present a comprehensive excited-state benchmark for the state-averaged (SA) driven similarity renormalization group (DSRG) [Li, C.; Evangelista, F. A. J. Chem. Phys. 2018, 148, 124106]. Following the QUEST database [Véril, M.; Scemama, A.; Caffarel, M.; Lipparini, F.; Boggio-Pasqua, M.; Jacquemin, D.; Loos, P.-F. Wiley Interdiscip. Rev. Comput. Mol. Sci. 2021, 11, e1517], 280 vertical transition energies of 35 medium-sized molecules are computed using the SA-DSRG derived second- and third-order perturbation theories (PT2/PT3) along with a nonperturbative approach [sq-LDSRG(2)]. Comparing to the theoretical best estimates, the optimal flow parameter is found to be 0.35 and 2.0 Eh-2 for SA-DSRG-PT2 and SA-DSRG-PT3, respectively. For SA-sq-LDSRG(2), a flow parameter of 1.5 Eh-2 provides converged equations without compromising the accuracy. We then assess the accuracy of the SA-DSRG hierarchy using these parameters. The SA-DSRG-PT2 scheme outperforms the level-shifted CASPT2 by 0.10 eV in mean absolute error (MAE), yet this accuracy is slightly inferior than that of CASPT2 with the ionization-potential-electron-affinity shift. Both SA-DSRG-PT3 and SA-sq-LDSRG(2) yield a MAE of 0.10 eV, which is comparable to that of CASPT3 (0.09 eV). Finally, we compute vertical excitation energies of several low-lying singlet states of nucleobases. The SA-sq-LDSRG(2) approach provides highly accurate results for π → π* excitations, while n → π* transitions are better described by SA-DSRG-PT3.

Nonadiabatic Vibrational Resonance Raman Spectra from Quantum Dynamics Propagations with LVC Models. Application to Thymine

We present a viable protocol to compute vibrational resonance Raman (vRR) spectra for systems with several close-lying and potentially coupled electronic states. It is based on the parametrization of linear vibronic coupling (LVC) models from time-dependent density functional theory (TD-DFT) calculations and quantum dynamics propagations of vibronic wavepackets with the multilayer version of the multiconfiguration time-dependent Hartree (ML-MCTDH) method. Our approach is applied to thymine considering seven coupled electronic states, comprising the three lowest bright states, and all vibrational coordinates. Computed vRR at different excitation wavelengths are in good agreement with the available experimental data. Up to 250 nm the signal is dominated by the lowest HOMO → LUMO transition, whereas at 233 nm, in the valley between the two lowest energy absorption bands, the contributions of all the three bright states, and their interferences and couplings, are important. Inclusion of solvent (water) effects improves the agreement with experiment, reproducing the coalescence of vibrational bands due to CC and C═O stretchings. With our approach we disentangle and assess the effect of interferences between the contribution of different quasi-resonant states to the transition polarizability and the effect of interstate couplings. Our findings strongly suggest that in cases of close-lying and potentially coupled states a simple inclusion of interference effects is not sufficient, and a fully nonadiabatic computation should instead be performed. We also document that for systems with strong couplings and quasi-degenerate states, the use of HT perturbative approach, not designed for these cases, may lead to large artifacts.

Automatized protocol and interface to simulate QM/MM time-resolved transient absorption at TD-DFT level with COBRAMM

We present a series of new implementations that we recently introduced in COBRAMM, the open-source academic software developed in our group. The goal of these implementations is to offer an automatized workflow and interface to simulate time-resolved transient absorption (TA) spectra of medium-to-big chromophore embedded in a complex environment. Therefore, the excited states absorption and the stimulated emission are simulated along nonadiabatic dynamics performed with trajectory surface hopping. The possibility of treating systems from medium to big size is given by the use of time-dependent density functional theory (TD-DFT) and the presence of the environment is taken into account employing a hybrid quantum mechanics/molecular mechanics (QM/MM) scheme. The full implementation includes a series of auxiliary scripts to properly setup the QM/MM system, the calculation of the wavefunction overlap along the dynamics for the propagation, the evaluation of the transition dipole moment at linear response TD-DFT level, and scripts to setup, run and analyze the TA from an ensemble of trajectories. Altogether, we believe that our implementation will open the door to the easily simulate the time-resolved TA of systems so far computationally inaccessible.

Quantum dynamics of the photoinduced charge separation in a symmetric donor–acceptor–donor triad: The role of vibronic couplings, symmetry and temperature

The photoinduced charge separation in a symmetric donor–acceptor–donor (D–A–D) triad is studied quantum mechanically using a realistic diabatic vibronic coupling model. The model includes a locally excited DA*D state and two charge-transfer states D+A−D and DA−D+ and is constructed according to a procedure generally applicable to semirigid D–A–D structures and based on energies, forces, and force constants obtained by quantum chemical calculations. In this case, the electronic structure is described by time-dependent density functional theory, and the corrected linear response is used in conjunction with the polarizable continuum model to account for state-specific solvent effects. The multimode dynamics following the photoexcitation to the locally excited state are simulated by the hybrid Gaussian-multiconfigurational time-dependent Hartree method, and temperature effects are included using thermo field theory. The dynamics are connected to the transient absorption spectrum obtained in recent experiments, which is simulated and fully assigned from first principles. It is found that the charge separation is mediated by symmetry-breaking vibrations of relatively low frequency, which implies that temperature should be accounted for to obtain reliable estimates of the charge transfer rate.