Excited state gradients within a polarizable QM/MM formulation

Multiscale approaches that partition the system into an active site (where the electronic process under study occurs) and a remaining region, the environment, have proven to be good strategies for the computation of electronic excitations in complex systems. In this work the implementation of a polarizable QM/MM scheme for the computation of excited state gradients is presented and are applied to a test case

Deciphering the Chemical Basis of Fluorescence of a Selenium-Labeled Uracil Probe when Bound at the Bacterial Ribosomal A-Site

We unveil in this work the main factors that govern the turn-on/off fluorescence of a Se-modified uracil probe at the ribosomal RNA A-site. Whereas the constraint into an “in-plane” conformation of the two rings of the fluorophore is the main driver for the observed turn-on fluorescence emission in the presence of the antibiotic paromomycin, the electrostatics of the environment plays a minor role during the emission process. Our computational strategy clearly indicates that, in the absence of paromomycin, the probe prefers conformations that show a dark S1 electronic state with participation of nπ* electronic transition contributions between the selenium atom and the π-system of the uracil moiety

Surface hopping within an exciton picture. An electrostatic embedding scheme

We report the development and the implementation of an exciton approach that allows ab initio nonadiabatic dynamics simulations of electronic excitation energy transfer in multichromophoric systems. For the dynamics, a trajectory-based strategy is used within the surface hopping formulation. The approach features a consistent hybrid formulation that allows the construction of potential energy surfaces and gradients by combining quantum mechanics and molecular mechanics within an electrostatic embedding scheme. As an application, the study of a molecular dyad consisting of a covalently bound BODIPY moiety and a tetrathiophene group is presented using time-dependent density functional theory (TDDFT). The results obtained with the exciton model are compared to previously performed full TDDFT dynamics of the same system. Our results show excellent agreement with the full TDDFT results, indicating that the couplings that lead to excitation energy transfer (EET) are dominated by Coulomb interaction terms and that charge-transfer states are not necessary to properly describe the nonadiabatic dynamics of the system. The exciton model also reveals ultrafast coherent oscillations of the excitation between the two units in the dyad, which occur during the first 50 fs

Highly efficient surface hopping dynamics using a linear vibronic coupling model.

We report an implementation of the linear vibronic coupling (LVC) model within the surface hopping dynamics approach and present utilities for parameterizing this model in a blackbox fashion. This results in an extremely efficient method to obtain qualitative and even semi-quantitative information about the photodynamical behavior of a molecule, and provides a new route toward benchmarking the results of surface hopping computations. The merits and applicability of the method are demonstrated in a number of applications. First, the method is applied to the SO2 molecule showing that it is possible to compute its absorption spectrum beyond the Condon approximation, and that all the main features and timescales of previous on-the-fly dynamics simulations of intersystem crossing are reproduced while reducing the computational effort by three orders of magnitude. The dynamics results are benchmarked against exact wavepacket propagations on the same LVC potentials and against a variation of the electronic structure level. Four additional test cases are presented to exemplify the broader applicability of the model. The photodynamics of the isomeric adenine and 2-aminopurine molecules are studied and it is shown that the LVC model correctly predicts ultrafast decay in the former and an extended excited-state lifetime in the latter. Futhermore, the method correctly predicts ultrafast intersystem crossing in the modified nucleobase 2-thiocytosine and its absence in 5-azacytosine while it fails to describe the ultrafast internal conversion to the ground state in the latter

Solvent effects on electronically excited states: QM/Continuum versus QM/Explicit models

The inclusion of solvent effects in the calculation of excited states is vital to obtain reliable absorption spectra and density of states of solvated chromophores. Here we analyze the performance of three classical approaches to describe aqueous solvent in the calculation of the absorption spectra and density of states of pyridine, tropone, and tropothione. Specifically, we compare the results obtained from quantum mechanics/polarizable continuum model (QM/PCM) versus quantum mechanics/molecular mechanics (QM/MM) in its electrostatic-embedding (QM/MMee) and polarizable-embedding (QM/MMpol) fashions, against full-QM computations, in which the solvent is described at the same level of theory as the chromophore. We show that QM/PCM provides very accurate results describing the excitation energies of ππ∗ and nπ∗ transitions, the last ones dominated by strong hydrogen-bonding effects, for the three chromophores. The QM/MMee approach also performs very well for both types of electronic transitions, although the description of the ππ∗ ones is slightly worse than that obtained from QM/PCM. The QM/MMpol approach performs as well as QM/PCM for describing the energy of ππ∗ states, but it is not able to provide a satisfactory description of hydrogen-bonding effects on the nπ∗ states of pyridine and tropone. The relative intensity of the absorption bands is better accounted for by the explicit-solvent models than by the continuum-solvent approach.LG and JJN further acknowledge the University of Vienna for financial support, while MDV and MFSJM thank the Marie Curie Actions, within the Innovative Training Network-European Join Doctorate in Theoretical Chemistry and Computational Modelling TCCM-ITN-EJD-642294, for their respective PhD grant

HOAX: a hyperparameter optimisation algorithm explorer for neural networks

Computational chemistry has become an important tool to predict and understand molecular properties and reactions. Even though recent years have seen a significant growth in new algorithms and computational methods that speed up quantum chemical calculations, the bottleneck for trajectory-based methods to study photo-induced processes is still the huge number of electronic structure calculations. In this work, we present an innovative solution, in which the amount of electronic structure calculations is drastically reduced, by employing machine learning algorithms and methods borrowed from the realm of artificial intelligence. However, applying these algorithms effectively requires finding optimal hyperparameters, which remains a challenge itself. Here we present an automated user-friendly framework, HOAX, to perform the hyperparameter optimisation for neural networks, which bypasses the need for a lengthy manual process. The neural network-generated potential energy surfaces (PESs) reduce the computational costs compared to the ab initio-based PESs. We perform a comparative investigation on the performance of different hyperparameter optimiziation algorithms, namely grid search, simulated annealing, genetic algorithm, and Bayesian optimizer in finding the optimal hyperparameters necessary for constructing the well-performing neural network in order to fit the PESs of small organic molecules. Our results show that this automated toolkit not only facilitate a straightforward way to perform the hyperparameter optimisation but also the resulting neural networks-based generated PESs are in reasonable agreement with the ab initio-based PESs

Methods to Calculate Electronic Excited-State Dynamics for Molecules on Large Metal Clusters with Many States:Ensuring Fast Overlap Calculations and a Robust Choice of Phase

We present an efficient set of methods for propagating excited-state dynamics involving a large number of configuration interaction singles (CIS) or Tamm-Dancoff approximation (TDA) single-reference excited states. Specifically, (i) following Head-Gordon et al., we implement an exact evaluation of the overlap of singly-excited CIS/TDA electronic states at different nuclear geometries using a biorthogonal basis and (ii) we employ a unified protocol for choosing the correct phase for each adiabat at each geometry. For many-electron systems, the combination of these techniques significantly reduces the computational cost of integrating the electronic Schrodinger equation and imposes minimal overhead on top of the underlying electronic structure calculation. As a demonstration, we calculate the electronic excited-state dynamics for a hydrogen molecule scattering off a silver metal cluster, focusing on high-lying excited states, where many electrons can be excited collectively and crossings are plentiful. Interestingly, we find that the high-lying, plasmon-like collective excitation spectrum changes with nuclear dynamics, highlighting the need to simulate non-adiabatic nuclear dynamics and plasmonic excitations simultaneously. In the future, the combination of methods presented here should help theorists build a mechanistic understanding of plasmon-assisted charge transfer and excitation energy relaxation processes near a nanoparticle or metal surface