This article summarizes technical advances contained in the fifth major release of the Q-Chem quantum chemistry program package, covering developments since 2015. A comprehensive library of exchange–correlation functionals, along with a suite of correlated many-body methods, continues to be a hallmark of the Q-Chem software. The many-body methods include novel variants of both coupled-cluster and configuration-interaction approaches along with methods based on the algebraic diagrammatic construction and variational reduced density-matrix methods. Methods highlighted in Q-Chem 5 include a suite of tools for modeling core-level spectroscopy, methods for describing metastable resonances, methods for computing vibronic spectra, the nuclear–electronic orbital method, and several different energy decomposition analysis techniques. High-performance capabilities including multithreaded parallelism and support for calculations on graphics processing units are described. Q-Chem boasts a community of well over 100 active academic developers, and the continuing evolution of the software is supported by an “open teamware” model and an increasingly modular design.

Photo-induced processes are fundamental in nature but accurate simulations are seriously limited by the cost of the underlying quantum chemical calculations, hampering their application for long time scales. Here we introduce a method based on machine learning to overcome this bottleneck and enable accurate photodynamics on nanosecond time scales, which are otherwise out of reach with contemporary approaches. Instead of expensive quantum chemistry during molecular dynamics simulations, we use deep neural networks to learn the relationship between a molecular geometry and its high-dimensional electronic properties. As an example, the time evolution of the methylenimmonium cation for one nanosecond is used to demonstrate that machine learning algorithms can outperform standard excited-state molecular dynamics approaches in their computational efficiency while delivering the same accuracy.

A highly efficient protocol for performing nonadiabatic dynamics simulations is implemented and applied to ultrafast internal conversion and intersystem crossing in various molecules.

We present a computational strategy to simulate the absorption lineshape of a molecule embedded in a complex environment by using a polarizable QM/MM approach. This strategy is presented in two alternative formulations, one based on a molecular dynamics simulation of the structural fluctuations of the system and the other using normal modes and harmonic frequencies calculated on optimized geometries. The comparison for the case of a chromophore within a strongly inhomogeneous and structured environment, namely the intercalation pocket of DNA, shows that the MD-based approach is able to reproduce the experimental spectral bandshape. In contrast, the static approach overestimates the vibronic coupling, resulting in a much broader band.

We report the development and the implementation of an exciton approach that allows the computing of ab initio non-adiabatic 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 also 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 non-adiabatic 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 rst 50 fs.

In a systematic study of the Au‐catalyzed reaction of o‐alkynylphenols with aryldiazonium salts, we find that essentially the same reaction conditions lead to a change in mechanism when a light source is applied. If the reaction is carried out at room temperature using a AuI catalyst, the diazonium salt undergoes electrophilic deauration of a vinyl AuI intermediate and provides access to substituted azobenzofurans. If the reaction mixture is irradiated with blue LED light, C−C bond formation due to N2‐extrusion from the diazonium salt is realized selectively, using the same starting materials without the need for an additional photo(redox) catalyst under aerobic conditions. We report a series of experiments demonstrating that the same vinyl AuI intermediate is capable of producing the observed products under photolytic and thermal conditions. The finding that a vinyl AuI complex can directly, without the need for an additional photo(redox) catalyst, result in C−C bond formation under photolytic conditions is contrary to the proposed mechanistic pathways suggested in the literature till date and highlights that the role of oxidation state changes in photoredox catalysis involving Au is thus far only poorly understood and may hold surprises for the future. Computational results indicate that photochemical activation can occur directly from a donor–acceptor complex formed between the vinyl AuI intermediate and the diazonium salt.

No abstract

Charge and structural relaxation of electronically excited states in embedded systems are strongly affected by the environment. It is known that the largest part of environment effects comes from electrostatics. However, polarization can also play a role by tuning the electronic and geometrical properties of the states, finally modifying the fluorescence. Here we present the formulation of analytical excited-state gradients within a polarizable QM/MM approach and their implementation within the ONIOM framework. A time-dependent DFT level of theory is used in combination with an induced dipole formulation of the polarizable embedding. The formation and relaxation of the bright excited state of an organic dye (DAPI) intercalated in a DNA pocket is used to quantify the role played by the mutual polarization between the QM subsystem and the embedding and also to investigate the onset of overpolarization, which is a known limit of the model with potentially detrimental effects. On the one hand, the results indicate the robustness of the QM-classical interface and, on the other hand, show the non-negligible effect of polarization between DAPI and a DNA pocket in determining the fluorescence properties of the embedded dye.