QM/MM Nonadiabatic Dynamics: the SHARC/COBRAMM Approach

Nonadiabatic Excited-State Molecular Dynamics with an Explicit Solvent: NEXMD-SANDER Implementation

In this article, the nonadiabatic excited-state Molecular dynamics (NEXMD) package is linked with the SANDER package, provided by AMBERTOOLS. The combination of these software packages enables the simulation of photoinduced dynamics of large multichromophoric conjugated molecules involving several coupled electronic excited states embedded in an explicit solvent by using the quantum/mechanics/molecular mechanics (QM/MM) methodology. The fewest switches surface hopping algorithm, as implemented in NEXMD, is used to account for quantum transitions among the adiabatic excited-state simulations of the photoexcitation and subsequent nonadiabatic electronic transitions, and vibrational energy relaxation of a substituted polyphenylenevinylene oligomer (PPV3-NO2) in vacuum and methanol as an explicit solvent has been used as a test case. The impact of including specific solvent molecules in the QM region is also analyzed. Our NEXMD-SANDER QM/MM implementation provides a useful computational tool to simulate qualitatively solvent-dependent effects, like electron transfer, stabilization of charge-separated excited states, and the role of solvent reorganization in the molecular optical properties, observed in solution-based spectroscopic experiments.

Dissociation Energies via Embedding Techniques

Due to the large number of interactions, evaluating interaction energies for large or periodic systems results in time-consuming calculations. Prime examples are liquids, adsorbates, and molecular crystals. Thus, there is a high demand for a cheap but still accurate method to determine interaction energies and gradients. One approach to counteract the computational cost is to fragment a large cluster into smaller subsystems, sometimes called many-body expansion, with the fragments being molecules or parts thereof. These subsystems can then be embedded into larger entities, representing the bigger system. In this work, we test several subsystem approaches and explore their limits and behaviors, determined by calculations of trimer interaction energies. The methods presented here encompass mechanical embedding, point charges, polarizable embedding, polarizable density embedding, and density embedding. We evaluate nonembedded fragmentation, QM/MM (quantum mechanics/molecular mechanics), and QM/QM (quantum mechanics/quantum mechanics) embedding theories. Finally, we make use of symmetry-adapted perturbation theory utilizing density functional theory for the monomers to interpret the results. Depending on the strength of the interaction, different embedding methods and schemes prove favorable to accurately describe a system. The embedding approaches presented here are able to decrease the interaction energy errors with respect to full system calculations by a factor of up to 20 in comparison to simple/unembedded approaches, leading to errors below 0.1 kJ/mol.

Ultrafast dynamics in spatially confined photoisomerization: accelerated simulations through machine learning models

This study sheds light on the exploration of photoresponsive host-guest systems, highlighting the intricate interplay between confined spaces and photosensitive guest molecules. Conducting nonadiabatic molecular dynamics (NAMD) simulations based on electronic structure calculations for such large systems remains a formidable challenge. By leveraging machine learning (ML) as an accelerator for NAMD simulations, we analytically constructed excited-state potential energy surfaces along relevant collective variables to investigate photoisomerization processes efficiently. Combining the quantum mechanics/molecular mechanics (QM/MM) methodology with ML-based NAMD simulations, we elucidated the reaction pathways and identified the key degrees of freedom as reaction coordinates leading to conical intersections. A machine learning-based nonadiabatic dynamics model has been developed to compare the excited-state dynamics of the guest molecule, benzopyran, in both the gas phase and its behavior within the confined space of cucurbit[5]uril. This comparative analysis was designed to determine the influence of the environment on the photoisomerization rate of the guest molecule. The results underscore the effectiveness of ML models in simulating trajectory evolution in a cost-effective manner. This research offers a practical approach to accelerate NAMD simulations in large-scale systems of photochemical reactions, with potential applications in other host-guest complex systems.

Simulation of the non-adiabatic dynamics of an enone-Lewis acid complex in an explicit solvent

Unlocking the full potential of Lewis acid catalysis for photochemical transformations requires a comprehensive understanding of the ultrafast dynamics of substrate-Lewis acid complexes. In a previous article [Peschel et al., Angew. Chem. Int. Ed., 2021, 60, 10155], time-resolved spectroscopy supported by static calculations revealed that the Lewis acid remains attached during the relaxation of the model complex cyclohexenone-BF3. In contrast to the experimental observation, surface-hopping dynamics in the gas phase predicted ultrafast heterolytic dissociation. We attributed the discrepancy to missing solvent interactions. Thus, in this work, we present an interface between the SHARC and FermiONs++ program packages, which enables us to investigate the ultrafast dynamics of cyclohexenone-BF3 in an explicit solvent environment. Our simulations demonstrate that the solvent prevents the dissociation of the complex, leading to an intriguing dissociation-reassociation mechanism. Comparing the dynamics with and without triplet states highlights their role in the relaxation process and shows that the Lewis acid inhibits intersystem crossing. These findings provide a clear picture of the relaxation process, which may aid in designing future Lewis acid catalysts for photochemical applications. They underscore that an explicit solvent model is required to describe relaxation processes in weakly bound states, as energy transfer to the solvent is crucial for the system to reach its minimum geometries.

Multiconfigurational Surface Hopping: a Time-Dependent Variational Approach with Momentum-Jump Trajectories

The Ehrenfest mean field dynamics and trajectory surface hopping have been widely used in nonadiabatic dynamics simulations. Based on the time-dependent variational principle (TDVP), the multiconfigurational Ehrenfest (MCE) method has also been developed and can be regarded as a multiconfigurational extension of the traditional Ehrenfest dynamics. However, it is not straightforward to apply the TDVP to surface hopping trajectories because there exists momentum jump during surface hops. To solve this problem, we here propose a multiconfigurational surface hopping (MCSH) method, where continuous momenta are obtained by linear interpolation and the interpolated trajectories are used to construct the basis functions for TDVP in a postprocessing manner. As demonstrated in a series of representative spin-boson models, MCSH achieves high accuracy with only several hundred trajectory bases and can uniformly improve the performance of surface hopping. In principle, MCSH can be combined with all kinds of mixed quantum-classical trajectories and thus has the potential to properly describe general nonadiabatic dynamics.

The development of the QM/MM interface and its application for the on-the-fly QM/MM nonadiabatic dynamics in JADE package: Theory, implementation, and applications

Understanding the nonadiabatic dynamics of complex systems is a challenging task in computational photochemistry. Herein, we present an efficient and user-friendly quantum mechanics/molecular mechanics (QM/MM) interface to run on-the-fly nonadiabatic dynamics. Currently, this interface consists of an independent set of codes designed for general-purpose use. Herein, we demonstrate the ability and feasibility of the QM/MM interface by integrating it with our long-term developed JADE package. Tailored to handle nonadiabatic processes in various complex systems, especially condensed phases and protein environments, we delve into the theories, implementations, and applications of on-the-fly QM/MM nonadiabatic dynamics. The QM/MM approach is established within the framework of the additive QM/MM scheme, employing electrostatic embedding, link-atom inclusion, and charge-redistribution schemes to treat the QM/MM boundary. Trajectory surface-hopping dynamics are facilitated using the fewest switches algorithm, encompassing classical and quantum treatments for nuclear and electronic motions, respectively. Finally, we report simulations of nonadiabatic dynamics for two typical systems: azomethane in water and the retinal chromophore PSB3 in a protein environment. Our results not only illustrate the power of the QM/MM program but also reveal the important roles of environmental factors in nonadiabatic processes.

Resolving Photoinduced Femtosecond Three-Dimensional Solute-Solvent Dynamics through Surface Hopping Simulations

Photoinduced dynamics in solution is governed by mutual solute-solvent interactions, which give rise to phenomena like solvatochromism, the Stokes shift, dual fluorescence, or charge transfer. Understanding these phenomena requires simulating the solute's photoinduced dynamics and simultaneously resolving the three-dimensional solvent distribution dynamics. If using trajectory surface hopping (TSH) to this aim, thousands of trajectories are required to adequately sample the time-dependent three-dimensional solvent distribution functions, and thus resolve the solvent dynamics with sub-Ångstrom and femtosecond accuracy and sufficiently low noise levels. Unfortunately, simulating thousands of trajectories with TSH in the framework of hybrid quantum mechanical/molecular mechanical (QM/MM) can be prohibitively expensive when employing ab initio electronic structure methods. To tackle this challenge, we recently introduced a computationally efficient approach that combines efficient linear vibronic coupling models with molecular mechanics (LVC/MM) via electrostatic embedding [Polonius et al., JCTC 2023, 19, 7171-7186]. This method provides solvent-embedded, nonadiabatically coupled potential energy surfaces while scaling similarly to MM force fields. Here, we employ TSH with LVC/MM to unravel the photoinduced dynamics of two small thiocarbonyl compounds solvated in water. We describe how to estimate the number of trajectories required to produce nearly noise-free three-dimensional solvent distribution functions and present an analysis based on approximately 10,000 trajectories propagated for 3 ps. In the electronic ground state, both molecules exhibit in-plane hydrogen bonds to the sulfur atom. Shortly after excitation, these bonds are broken and reform perpendicular to the molecular plane on timescales that differ by an order of magnitude due to steric effects. We also show that the solvent relaxation dynamics is coupled to the electronic dynamics, including intersystem crossing. These findings are relevant to advance the understanding of the coupled solute-solvent dynamics of solvated photoexcited molecules, e.g., biologically relevant thio-nucleobases.

A Portrait of the Chromophore as a Young System-Quantum-Derived Force Field Unraveling Solvent Reorganization upon Optical Excitation of Cyclocurcumin Derivatives

The study of fast non-equilibrium solvent relaxation in organic chromophores is still challenging for molecular modeling and simulation approaches, and is often overlooked, even in the case of non-adiabatic dynamics simulations. Yet, especially in the case of photoswitches, the interaction with the environment can strongly modulate the photophysical outcomes. To unravel such a delicate interplay, in the present contribution we resorted to a mixed quantum-classical approach, based on quantum mechanically derived force fields. The main task is to rationalize the solvent reorganization pathways in chromophores derived from cyclocurcumin, which are suitable for light-activated chemotherapy to destabilize cellular lipid membranes. The accurate and reliable decryption delivered by the quantum-derived force fields points to important differences in the solvent's reorganization, in terms of both structure and time scale evolution.

Potential energy interpolation with target-customized weighting coordinates: application to excited-state dynamics of photoactive yellow protein chromophore in water

Interpolation of potential energy surfaces (PESs) can provide a practical route to performing molecular dynamics simulations with a reliability matching a high-level quantum chemical calculation. An obstacle to its widespread use is perhaps the lack of general and optimal interpolation settings that can be applied in a black-box manner for any given molecular system. How to set up the weights for interpolation is one such task, and we still need to diversify the approaches in order to treat various systems. Here, we develop a new interpolation weighting scheme, which allows us to choose the weighting coordinates in a system-specific manner, by amplifying the contribution from specific internal coordinates. The new weighting scheme with an appropriate selection of coordinates is proved to be effective in reducing the interpolation error along the reaction pathway. As a demonstration, we consider the photoactive yellow protein chromophore system, as it constitutes itself as an interesting target that bears long-standing questions related to excited-state dynamics inside protein environments. We build its two-state diabatic interpolated PES with the new weighting scheme. We indeed see the utility of our scheme by conducting nonadiabatic molecular dynamics simulations with the required semi-global PES based on a limited number of data points.

Galván I, Alavi A, Aleotti F, Aquilante F, Autschbach J, Avagliano D, Baiardi A, Bao JJ, Battaglia S, Birnoschi L, Blanco-González A, Bokarev SI, Broer R, Cacciari R, Calio PB, Carlson RK, Carvalho Couto R, Cerdán L, Chibotaru LF, Chilton NF, Church JR, Conti I, Coriani S, Cuéllar-Zuquin J, Daoud RE, Dattani N, Decleva P, de Graaf C, Delcey M, De Vico L, Dobrautz W, Dong SS, Feng R, Ferré N, Filatov(Gulak) M, Gagliardi L, Garavelli M, González L, Guan Y, Guo M, Hennefarth MR, Hermes MR, Hoyer CE, Huix-Rotllant M, Jaiswal VK, Kaiser A, Kaliakin DS, Khamesian M, King DS, Kochetov V, Krośnicki M, Kumaar AA, Larsson ED, Lehtola S, Lepetit MB, Lischka H, López Ríos P, Lundberg M, Ma D, Mai S, Marquetand P, Merritt ICD, Montorsi F, Mörchen M, Nenov A, Nguyen VHA, Nishimoto Y, Oakley MS, Olivucci M, Oppel M, Padula D, Pandharkar R, Phung QM, Plasser F, Raggi G, Rebolini E, Reiher M, Rivalta I, Roca-Sanjuán D, Romig T, Safari AA, Sánchez-Mansilla A, Sand AM, Schapiro I, Scott TR, Segarra-Martí J, Segatta F, Sergentu DC, Sharma P, Shepard R, Shu Y, Staab JK, Straatsma TP, Sørensen LK, Tenorio BNC, Truhlar DG, Ungur L, Vacher M, Veryazov V, Voß TA, Weser O, Wu D, Yang X, Yarkony D, Zhou C, Zobel JP, Lindh R. The OpenMolcas Web: A Community-Driven Approach to Advancing Computational Chemistry

The developments of the open-source OpenMolcas chemistry software environment since spring 2020 are described, with a focus on novel functionalities accessible in the stable branch of the package or via interfaces with other packages. These developments span a wide range of topics in computational chemistry and are presented in thematic sections: electronic structure theory, electronic spectroscopy simulations, analytic gradients and molecular structure optimizations, ab initio molecular dynamics, and other new features. This report offers an overview of the chemical phenomena and processes OpenMolcas can address, while showing that OpenMolcas is an attractive platform for state-of-the-art atomistic computer simulations.

LVC/MM: A Hybrid Linear Vibronic Coupling/Molecular Mechanics Model with Distributed Multipole-Based Electrostatic Embedding for Highly Efficient Surface Hopping Dynamics in Solution

We present a theoretical framework for a hybrid linear vibronic coupling model electrostatically embedded into a molecular mechanics environment, termed the linear vibronic coupling/molecular mechanics (LVC/MM) method, for the surface hopping including arbitrary coupling (SHARC) molecular dynamics package. Electrostatic embedding is realized through the computation of interactions between environment point charges and distributed multipole expansions (DMEs, up to quadrupoles) that represent each electronic state and transition densities in the diabatic basis. The DME parameters are obtained through a restrained electrostatic potential (RESP) fit, which we extended to yield higher-order multipoles. We also implemented in SHARC a scheme for achieving roto-translational invariance of LVC models as well as a general quantum mechanics/molecular mechanics (QM/MM) interface, an OpenMM interface, and restraining potentials for simulating liquid droplets. Using thioformaldehyde in water as a test case, we demonstrate that LVC/MM can accurately reproduce the solvation structure and energetics of rigid solutes, with errors on the order of 1-2 kcal/mol compared to a BP86/MM reference. The implementation in SHARC is shown to be very efficient, enabling the simulation of trajectories on the nanosecond time scale in a matter of days.

Nature of Hops, Coordinates, and Detailed Balance for Nonadiabatic Simulations in the Condensed Phase

Photoinduced processes play a crucial role in a multitude of important molecular phenomena. Accurately modeling these processes in an environment other than a vacuum requires a detailed description of the electronic states involved as well as how energy flows are coupled to the surroundings. Nonadiabatic effects must also be included in order to describe the exchange of energy between electronic and nuclear degrees of freedom correctly. In this work, we revisit the ring-opening reaction 1,3-cylohexadiene (CHD) in a solvent environment. Using our newly developed Interface for Non-Adiabatic Quantum mechanics/molecular mechanics in Solvent (INAQS) we trace the evolution of the reaction via hybrid quantum mechanics/molecular mechanics (QM/MM) surface hopping with a focus on the solvent's participation in the nonadiabatic relaxation process and the long-time approach to equilibrium. We explicitly include the MM solvent contribution to the nonadiabatic coupling vector─enabling an accurate approach to equilibrium at long times─and find that in highly multidimensional systems gradients can have little or nothing to do with the nonadiabatic couplings.

NEXMD v2.0 Software Package for Nonadiabatic Excited State Molecular Dynamics Simulations

We present NEXMD version 2.0, the second release of the NEXMD (Nonadiabatic EXcited-state Molecular Dynamics) software package. Across a variety of new features, NEXMD v2.0 incorporates new implementations of two hybrid quantum-classical dynamics methods, namely, Ehrenfest dynamics (EHR) and the Ab-Initio Multiple Cloning sampling technique for Multiconfigurational Ehrenfest quantum dynamics (MCE-AIMC or simply AIMC), which are alternative options to the previously implemented trajectory surface hopping (TSH) method. To illustrate these methodologies, we outline a direct comparison of these three hybrid quantum-classical dynamics methods as implemented in the same NEXMD framework, discussing their weaknesses and strengths, using the modeled photodynamics of a polyphenylene ethylene dendrimer building block as a representative example. We also describe the expanded normal-mode analysis and constraints for both the ground and excited states, newly implemented in the NEXMD v2.0 framework, which allow for a deeper analysis of the main vibrational motions involved in vibronic dynamics. Overall, NEXMD v2.0 expands the range of applications of NEXMD to a larger variety of multichromophore organic molecules and photophysical processes involving quantum coherences and persistent couplings between electronic excited states and nuclear velocity.

Deciphering Methylation Effects on S2(ππ*) Internal Conversion in the Simplest Linear α,β-Unsaturated Carbonyl

Chemical substituents can influence photodynamics by altering the location of critical points and the topography of the potential energy surfaces (electronic effect) and by selectively modifying the inertia of specific nuclear modes (inertial effects). Using nonadiabatic dynamics simulations, we investigate the impact of methylation on S₂(ππ*) internal conversion in acrolein, the simplest linear α,β-unsaturated carbonyl. Consistent with time constants reported in a previous time-resolved photoelectron spectroscopy study, S₂ → S₁ deactivation occurs on an ultrafast time scale (∼50 fs). However, our simulations do not corroborate the sequential decay model used to fit the experiment. Instead, upon reaching the S₁ state, the wavepacket bifurcates: a portion undergoes ballistic S₁ → S₀ deactivation (∼90 fs) mediated by fast bond-length alternation motion, while the remaining decays on the picosecond time scale. Our analysis reveals that methyl substitution, generally assumed to mainly exert inertial influence, is also manifested in important electronic effects due to its weak electron-donating ability. While methylation at the β C atom gives rise to effects principally of an inertial nature, such as retarding the twisting motion of the terminal -CHCH₃ group and increasing its coupling with pyramidalization, methylation at the α or carbonyl C atom modifies the potential energy surfaces in a way that also contributes to altering the late S₁-decay behavior. Specifically, our results suggest that the observed slowing of the picosecond component upon α-methylation is a consequence of a tighter surface and reduced amplitude along the central pyramidalization, effectively restricting the access to the S₁/S₀-intersection seam. Our work offers new insight into the S₂(ππ*) internal conversion mechanisms in acrolein and its methylated derivatives and highlights site-selective methylation as a tuning knob to manipulate photochemical reactions.

Structural Sampling and Solvation Models for the Simulation of Electronic Spectra: Pyrazine as a Case Study

The impact of sampling methods on spectral broadening in the gas phase and on the convergence of spectra in aqueous solution when using microsolvation, continuum solvation, and hybrid models is studied using pyrazine as a test case. For the sake of comparing classical Maxwell-Boltzmann and Wigner samplings in the gas phase, static and time-resolved X-ray absorption spectra after photoexcitation to the lowest ¹B_2u(ππ*) state, as well as the static UV-vis absorption spectrum, are considered. In addition, the UV-vis absorption spectrum of pyrazine in aqueous solution is also computed in order to systematically investigate its convergence with the number of explicitly included solvent shells with and without taking bulk solvation effects into account with the conductor-like screening model to represent implicit water beyond such explicit solute complexes. Concerning the static and time-resolved X-ray absorption spectra of pyrazine at the carbon K-edge as well as its UV-vis absorption spectrum in the gas phase, we find that these spectra obtained with Wigner and Maxwell-Boltzmann samplings substantially agree. For the UV-vis absorption spectrum in the aqueous solution, only the first two energetically low-lying bands converge quickly with the size of the explicitly included solvation shells, either with or without an additional continuum solvation taken into account. In stark contrast, calculations of the higher-lying excitations relying on finite microsolvated clusters without additional continuum solvation severely suffer from unphysical charge-transfer excitations into Rydberg-like orbitals at the cluster/vacuum interface. This finding indicates that computational UV-vis absorption spectra covering sufficiently high-lying states converge only if continuum solvation of the explicitly microsolvated solutes is included in the models.

Solvent Dynamics of Aqueous Halides before and after Photoionization

Electron transfer reactions can be strongly influenced by solvent dynamics. We study the photoionization of halides in water as a model system for such reactions. There are no internal nuclear degrees of freedom in the solute, allowing the dynamics of the solvent to be uniquely identified. We simulate the equilibrium solvent dynamics for Cl^-, Br^-, I^-, and their respective neutral atoms in water, comparing quantum mechanical/molecular mechanical (QM/MM) and classical molecular dynamics (MD) methods. On the basis of the obtained configurations, we calculate the extended X-ray absorption fine structure (EXAFS) spectra rigorously based on the MD snapshots and compare them in detail with other theoretical and experimental results available in the literature. We find our EXAFS spectra based on QM/MM MD simulations in good agreement with their experimental counterparts for the ions. Classical MD simulations for the ions lead to EXAFS spectra that agree equally well with the experiment when it comes to the oscillatory period of the signal, even though they differ from the QM/MM radial distribution functions extracted from the MD. The amplitude is, however, considerably overestimated. This suggests that to judge the reliability of theoretical simulation methods or to elucidate fine details of the atomistic dynamics of the solvent based on EXAFS spectra, the amplitude as well as the oscillatory period need to be considered. If simulations fail qualitatively, as does the classical MD for the aqueous neutral halogen atoms, the resulting EXAFS will also be strongly affected in both oscillatory period and amplitude. The good reliability of QM/MM-based EXAFS simulations, together with clear qualitative differences in the EXAFS spectra found between halides and their atomic counterparts, suggests that a combined theory and experimental EXAFS approach is suitable for elucidating the nonequilibrium solvent dynamics in the photoionization of halides and possibly also for electron transfer reactions in more complex systems.

MoBioTools: A toolkit to setup quantum mechanics/molecular mechanics calculations

We present a toolkit that allows for the preparation of QM/MM input files from a conformational ensemble of molecular geometries. The package is currently compatible with trajectory and topology files in Amber, CHARMM, GROMACS and NAMD formats, and has the possibility to generate QM/MM input files for Gaussian (09 and 16), Orca (≥4.0), NWChem and (Open)Molcas. The toolkit can be used in command line, so that no programming experience is required, although it presents some features that can also be employed as a python application programming interface. We apply the toolkit in four situations in which different electronic-structure properties of organic molecules in the presence of a solvent or a complex biological environment are computed: the reduction potential of the nucleobases in acetonitrile, an energy decomposition analysis of tyrosine interacting with water, the absorption spectrum of an azobenzene derivative integrated into a voltage-gated ion channel, and the absorption and emission spectra of the luciferine/luciferase complex. These examples show that the toolkit can be employed in a manifold of situations for both the electronic ground state and electronically excited states. It also allows for the automatic correction of the active space in the case of CASSCF calculations on an ensemble of geometries, as it is shown for the azobenzene derivative photoswitch case.

State averaged CASSCF in AMOEBA polarizable water model for simulating nonadiabatic molecular dynamics with nonequilibrium solvation effects

This paper presents a state-averaged complete active space self-consistent field (SA-CASSCF) in the atomic multipole optimized energetics for biomolecular application (AMOEBA) polarizable water model, which enables rigorous simulation of non-adiabatic molecular dynamics with nonequilibrium solvation effects. The molecular orbital and configuration interaction coefficients of the solute wavefunction, and the induced dipoles on solvent atoms, are solved by minimizing the state averaged energy variationally. In particular, by formulating AMOEBA water models and the polarizable continuum model (PCM) in a unified way, the algorithms developed for computing SA-CASSCF/PCM energies, analytical gradients, and non-adiabatic couplings in our previous work can be generalized to SA-CASSCF/AMOEBA by properly substituting a specific list of variables. Implementation of this method will be discussed with the emphasis on how the calculations of different terms are partitioned between the quantum chemistry and molecular mechanics codes. We will present and discuss results that demonstrate the accuracy and performance of the implementation. Next, we will discuss results that compare three solvent models that work with SA-CASSCF, i.e., PCM, fixed-charge force fields, and the newly implemented AMOEBA. Finally, the new SA-CASSCF/AMOEBA method has been interfaced with the ab initio multiple spawning method to carry out non-adiabatic molecular dynamics simulations. This method is demonstrated by simulating the photodynamics of the model retinal protonated Schiff base molecule in water.

Exciton-vibrational dynamics induces efficient self-trapping in a substituted nanoring

Cycloparaphenylenes, being the smallest segments of carbon nanotubes, have emerged as prototypes of the simplest carbon nanohoops. Their unique structure-dynamics-optical properties relationships have motivated a wide variety of synthesis of new related nanohoop species. Studies of how chemical changes, introduced in these new materials, lead to systems with new structural, dynamics and optical properties, expand their functionalities for optoelectronics applications. Herein, we study the effect that conjugation extension of a cycloparaphenylene through the introduction of a satellite tetraphenyl substitution has on its structural and dynamical properties. Our non-adiabatic excited state molecular dynamics simulations suggest that this substitution accelerates the electronic relaxation from the high-energy band to the lowest excited state. This is partially due to efficient conjugation achieved between specific phenyl units as introduced by the tetraphenyl substitution. We observe a particular exciton redistribution during relaxation, in which the tetraphenyl substitution plays a significant role. As a result, an efficient inter-band energy transfer takes place. Besides, the observed phonon-exciton interplay induces a significant exciton self-trapping. Our results encourage and guide the future studies of new phenyl substitutions in carbon nanorings with desired optoelectronic properties.

Infinitene: Computational Insights from Nonadiabatic Excited State Dynamics

Progress in organic synthesis opens exploration of a rich diversity of molecules with interesting new structural topologies. This is the case of a recently synthesized helically twisted figure-eight molecule coined infinitene. The molecule belongs to a numerous family of looped polyarenes, where the degree of π-conjugation is controlled by high strain energies and steric hindrances. A particular balance of these ingredients leads to unusual optoelectronic properties potentially suitable for a range of applications in nanoelectronics and photonics. Due to its recent discovery, the photophysical properties of infinitene remain unexplored. In this Letter, atomistic nonadiabatic excited state molecular dynamics modeling unveils unique features of intramolecular electronic and vibrational energy relaxation and redistribution that take place after molecular photoexcitation. Our results detail relationships between optical and electronic properties providing useful knowledge for future molecular designs related to infinitene.

Adiabatic Excited-State Molecular Dynamics with an Explicit Solvent: NEXMD-SANDER Implementation

We present a method to link the Nonadiabatic EXcited-state Molecular Dynamics (NEXMD) package to the SANDER package supplied by AMBERTOOLS to provide excited-state adiabatic quantum mechanics/molecular mechanics (QM/MM) simulations. NEXMD is a computational package particularly developed to perform simulations of the photoexcitation and subsequent nonadiabatic electronic and vibrational energy relaxation in large multichromophoric conjugated molecules involving several coupled electronic excited states. The NEXMD-SANDER exchange has been optimized in order to achieve excited-state adiabatic dynamics simulations of large conjugated materials in a QM/MM environment, such as an explicit solvent. Dynamics of a substituted polyphenylene vinylene oligomer (PPV3-NO₂) in vacuum and different explicit solvents has been used as a test case by performing comparative analysis of changes in its optical spectrum, state-dependent conformational changes, and quantum bond orderings. The method has been tested and compared with respect to previous implicit solvent implementations. Also, the impact on the expansion of the QM region by including a variable number of solvent molecules has been analyzed. Altogether, these results encourage future implementations of NEXMD simulations using the same combination of methods.

INAQS, a Generic Interface for Nonadiabatic QM/MM Dynamics: Design, Implementation, and Validation for GROMACS/Q-CHEM simulations

The accurate description of large molecular systems in complex environments remains an ongoing challenge for the field of computational chemistry. This problem is even more pronounced for photoinduced processes, as multiple excited electronic states and their corresponding nonadiabatic couplings must be taken into account. Multiscale approaches such as hybrid quantum mechanics/molecular mechanics (QM/MM) offer a balanced compromise between accuracy and computational burden. Here, we introduce an open-source software package (INAQS) for nonadiabatic QM/MM simulations that bridges the sampling capabilities of the GROMACS MD package and the excited-state infrastructure of the Q-CHEM electronic structure software. The interface is simple and can be adapted easily to other MD codes. The code supports a variety of different trajectory-based molecular dynamics, ranging from Born-Oppenheimer to surface hopping dynamics. To illustrate the power of this combination, we simulate electronic absorption spectra, free-energy surfaces along a reaction coordinate, and the excited-state dynamics of 1,3-cyclohexadiene in solution.

Ultrafast Intersystem Crossing Dynamics of 6-Selenoguanine in Water

Rationalizing the photochemistry of nucleobases where an oxygen is replaced by a heavier atom is essential for applications that exploit near-unity triplet quantum yields. Herein, we report on the ultrafast excited-state deactivation mechanism of 6-selenoguanine (6SeGua) in water by combining nonadiabatic trajectory surface-hopping dynamics with an electrostatic embedding quantum mechanics/molecular mechanics (QM/MM) scheme. We find that the predominant relaxation mechanism after irradiation starts on the bright singlet S₂ state that converts internally to the dark S₁ state, from which the population is transferred to the triplet T₂ state via intersystem crossing and finally to the lowest T₁ state. This S₂ → S₁ → T₂ → T₁ deactivation pathway is similar to that observed for the lighter 6-thioguanine (6tGua) analogue, but counterintuitively, the T₁ lifetime of the heavier 6SeGua is shorter than that of 6tGua. This fact is explained by the smaller activation barrier to reach the T₁/S₀ crossing point and the larger spin-orbit couplings of 6SeGua compared to 6tGua. From the dynamical simulations, we also calculate transient absorption spectra (TAS), which provide two time constants (τ₁ = 131 fs and τ₂ = 191 fs) that are in excellent agreement with the experimentally reported value (τ_exp = 130 ± 50 fs) (Farrel et al. J. Am. Chem. Soc. 2018, 140, 11214). Intersystem crossing itself is calculated to occur with a time scale of 452 ± 38 fs, highlighting that the TAS is the result of a complex average of signals coming from different nonradiative processes and not intersystem crossing alone.

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.

Sampling effects in quantum mechanical/molecular mechanics trajectory surface hopping non-adiabatic dynamics

The impact of different initial conditions in non-adiabatic trajectory surface hopping dynamics within a hybrid quantum mechanical/molecular mechanics scheme is investigated. The influence of a quantum sampling, based on a Wigner distribution, a fully thermal sampling, based on classical molecular dynamics, and a quantum sampled system, but thermally equilibrated with the environment, is investigated on the relaxation dynamics of solvated fulvene after light irradiation. We find that the decay from the first singlet excited state to the ground state shows high dependency on the initial condition and simulation parameters. The three sampling methods lead to different distributions of initial geometries and momenta, which then affect the fate of the excited state dynamics. We evaluated both the effect of sampling geometries and momenta, analysing how the ultrafast decay of fulvene changes accordingly. The results are expected to be of interest to decide how to initialize non-adiabatic dynamics in the presence of the environment. This article is part of the theme issue 'Chemistry without the Born-Oppenheimer approximation'.

Nonadiabatic Dynamics Algorithms with Only Potential Energies and Gradients: Curvature-Driven Coherent Switching with Decay of Mixing and Curvature-Driven Trajectory Surface Hopping

Direct dynamics by mixed quantum-classical nonadiabatic methods is an important tool for understanding processes involving multiple electronic states. Very often, the computational bottleneck of such direct simulation comes from electronic structure theory. For example, at every time step of a trajectory, nonadiabatic dynamics requires potential energy surfaces, their gradients, and the matrix elements coupling the surfaces. The need for the couplings can be alleviated by employing the time derivatives of the wave functions, which can be evaluated from overlaps of electronic wave functions at successive time steps. However, evaluation of overlap integrals is still expensive for large systems. In addition, for electronic structure methods for which the wave functions or the coupling matrix elements are not available, nonadiabatic dynamics algorithms become inapplicable. In this work, building on recent work by Baeck and An, we propose new nonadiabatic dynamics algorithms that only require adiabatic potential energies and their gradients. The new methods are named curvature-driven coherent switching with decay of mixing (κCSDM) and curvature-driven trajectory surface hopping (κTSH). We show how powerful these new methods are in terms of computation time and accuracy as compared to previous mixed quantum-classical nonadiabatic dynamics algorithms. The lowering of the computational cost will allow longer nonadiabatic trajectories and greater ensemble averaging to be affordable, and the ability to calculate the dynamics without electronic structure coupling matrix elements extends the dynamics capability to new classes of electronic structure methods.

Scalable Ehrenfest Molecular Dynamics Exploiting the Locality of Density-Functional Tight-Binding Hamiltonian

To explore the science behind excited-state dynamics in high-complexity chemical systems, a scalable nonadiabatic molecular dynamics (MD) technique is indispensable. In this study, by treating the electronic degrees of freedom at the density-functional tight-binding level, we developed and implemented a reduced scaling and multinode-parallelizable Ehrenfest MD method. To achieve this goal, we introduced a concept called patchwork approximation (PA), where the effective Hamiltonian for real-time propagation of the electronic density matrix is partitioned into a set of local parts. Numerical results for giant icosahedral fullerenes, which comprise up to 6000 atoms, suggest that the scaling of the present PA-based method is less than quadratic, which yields a significant advantage over the conventional cubic scaling method in terms of computational time. The acceleration by the parallelization on multiple nodes was also assessed. Furthermore, the electronic and structural dynamics resulting from the perturbation by the external electric field were accurately reproduced with the PA, even when the electronic excitation was spatially delocalized.