Deep learning study of tyrosine reveals that roaming can lead to photodamage

Surprising Dynamics Phenomena in the Diels-Alder Reaction of C60 Uncovered with AI

We performed an extensive artificial intelligence-accelerated quasi-classical molecular dynamics investigation of the time-resolved mechanism of the Diels-Alder reaction of fullerene C60 with 2,3-dimethyl-1,3-butadiene. In a substantial fraction (10%) of reactive trajectories, the larger C60 noncovalently attracts the 2,3-dimethyl-1,3-butadiene long before the barrier so that the diene undergoes the series of complex motions including roaming, somersaults, twisting, and twisting somersaults around the fullerene until it aligns itself to pass over the barrier. These complicated processes could be easily missed in typically performed quantum chemical simulations with shorter and fewer trajectories. After the barrier is passed, the bonds take longer to form compared to the simplest prototypical Diels-Alder reaction of ethene with 1,3-butadiene despite high similarities in transition states and barrier widths evaluated with intrinsic reaction coordinate (IRC) calculations. C60 is mainly responsible for these differences as its reaction with 1,3-butadiene is similar to the reaction with 2,3-dimethyl-1,3-butadiene: the only substantial difference being that the extra methyl groups double the probability of the prolonged alignment phase in dynamics. These additional calculations of C60 with 1,3-butadiene could be performed via active learning more easily by reusing the data generated for the other two reactions, showing the potential for larger-scale exploration of the effects of different substrates in the same types of reactions.

On the Simulation of Photoreactions Using Restricted Open-Shell Kohn-Sham Theory

It is a well-established standard to describe ground-state chemical reactions at an ab initio level of multi-electron theory. Fast reactions can be directly simulated. The most widely used approach is density functional theory for the electronic structure in combination with molecular dynamics for the nuclear motion. This approach is known as ab initio molecular dynamics. In contrast, the simulation of excited-state reactions at this level of theory is significantly more difficult. It turns out that the self-consistent solution of the Kohn-Sham equations is not easily reached in excited-state simulations. The first program that solved this problem was the Car-Parrinello molecular dynamics code, using restricted open-shell Kohn-Sham theory. Meanwhile, there are alternatives, most prominently the Q-Chem code, which widens the range of applications. The present study investigates the suitability of both codes for the molecular dynamics simulation of excited-state motion and presents applications to photoreactions.

SpaiNN: equivariant message passing for excited-state nonadiabatic molecular dynamics

Excited-state molecular dynamics simulations are crucial for understanding processes like photosynthesis, vision, and radiation damage. However, the computational complexity of quantum chemical calculations restricts their scope. Machine learning offers a solution by delivering high-accuracy properties at lower computational costs. We present SpaiNN, an open-source Python software for ML-driven surface hopping nonadiabatic molecular dynamics simulations. SpaiNN combines the invariant and equivariant neural network architectures of SchNetPack with SHARC for surface hopping dynamics. Its modular design allows users to implement and adapt modules easily. We compare rotationally-invariant and equivariant representations in fitting potential energy surfaces of multiple electronic states and properties arising from the interaction of two electronic states. Simulations of the methyleneimmonium cation and various alkenes demonstrate the superior performance of equivariant SpaiNN models, improving accuracy, generalization, and efficiency in both training and inference.

Recent Advances in 3D Printing of Smart Scaffolds for Bone Tissue Engineering and Regeneration

The repair and functional reconstruction of bone defects resulting from severe trauma, surgical resection, degenerative disease, and congenital malformation pose significant clinical challenges. Bone tissue engineering (BTE) holds immense potential in treating these severe bone defects, without incurring prevalent complications associated with conventional autologous or allogeneic bone grafts. 3D printing technology enables control over architectural structures at multiple length scales and has been extensively employed to process biomimetic scaffolds for BTE. In contrast to inert and functional bone grafts, next-generation smart scaffolds possess a remarkable ability to mimic the dynamic nature of native extracellular matrix (ECM), thereby facilitating bone repair and regeneration. Additionally, they can generate tailored and controllable therapeutic effects, such as antibacterial or antitumor properties, in response to exogenous and/or endogenous stimuli. This review provides a comprehensive assessment of the progress of 3D-printed smart scaffolds for BTE applications. It begins with an introduction to bone physiology, followed by an overview of 3D printing technologies utilized for smart scaffolds. Notable advances in various stimuli-responsive strategies, therapeutic efficacy, and applications of 3D-printed smart scaffolds are discussed. Finally, the review highlights the existing challenges in the development and clinical implementation of smart scaffolds, as well as emerging technologies in this field.

In Silico Chemical Experiments in the Age of AI: From Quantum Chemistry to Machine Learning and Back

Computational chemistry is an indispensable tool for understanding molecules and predicting chemical properties. However, traditional computational methods face significant challenges due to the difficulty of solving the Schrödinger equations and the increasing computational cost with the size of the molecular system. In response, there has been a surge of interest in leveraging artificial intelligence (AI) and machine learning (ML) techniques to in silico experiments. Integrating AI and ML into computational chemistry increases the scalability and speed of the exploration of chemical space. However, challenges remain, particularly regarding the reproducibility and transferability of ML models. This review highlights the evolution of ML in learning from, complementing, or replacing traditional computational chemistry for energy and property predictions. Starting from models trained entirely on numerical data, a journey set forth toward the ideal model incorporating or learning the physical laws of quantum mechanics. This paper also reviews existing computational methods and ML models and their intertwining, outlines a roadmap for future research, and identifies areas for improvement and innovation. Ultimately, the goal is to develop AI architectures capable of predicting accurate and transferable solutions to the Schrödinger equation, thereby revolutionizing in silico experiments within chemistry and materials science.

MLatom Software Ecosystem for Surface Hopping Dynamics in Python with Quantum Mechanical and Machine Learning Methods

We present an open-source MLatom@XACS software ecosystem for on-the-fly surface hopping nonadiabatic dynamics based on the Landau-Zener-Belyaev-Lebedev algorithm. The dynamics can be performed via Python API with a wide range of quantum mechanical (QM) and machine learning (ML) methods, including ab initio QM (CASSCF and ADC(2)), semiempirical QM methods (e.g., AM1, PM3, OMx, and ODMx), and many types of ML potentials (e.g., KREG, ANI, and MACE). Combinations of QM and ML methods can also be used. While the user can build their own combinations, we provide AIQM1, which is based on Δ-learning and can be used out-of-the-box. We showcase how AIQM1 reproduces the isomerization quantum yield of trans-azobenzene at a low cost. We provide example scripts that, in dozens of lines, enable the user to obtain the final population plots by simply providing the initial geometry of a molecule. Thus, those scripts perform geometry optimization, normal mode calculations, initial condition sampling, parallel trajectories propagation, population analysis, and final result plotting. Given the capabilities of MLatom to be used for training different ML models, this ecosystem can be seamlessly integrated into the protocols building ML models for nonadiabatic dynamics. In the future, a deeper and more efficient integration of MLatom with Newton-X will enable a vast range of functionalities for surface hopping dynamics, such as fewest-switches surface hopping, to facilitate similar workflows via the Python API.

Machine Learning of Reactive Potentials

In the past two decades, machine learning potentials (MLPs) have driven significant developments in chemical, biological, and material sciences. The construction and training of MLPs enable fast and accurate simulations and analysis of thermodynamic and kinetic properties. This review focuses on the application of MLPs to reaction systems with consideration of bond breaking and formation. We review the development of MLP models, primarily with neural network and kernel-based algorithms, and recent applications of reactive MLPs (RMLPs) to systems at different scales. We show how RMLPs are constructed, how they speed up the calculation of reactive dynamics, and how they facilitate the study of reaction trajectories, reaction rates, free energy calculations, and many other calculations. Different data sampling strategies applied in building RMLPs are also discussed with a focus on how to collect structures for rare events and how to further improve their performance with active learning.

Artificial-Intelligence-Enhanced On-the-Fly Simulation of Nonlinear Time-Resolved Spectra

Time-resolved spectroscopy is an important tool for unraveling the minute details of structural changes in molecules of biological and technological significance. The nonlinear femtosecond signals detected for such systems must be interpreted, but it is a challenging task for which theoretical simulations are often indispensable. Accurate simulations of transient absorption or two-dimensional electronic spectra are, however, computationally very expensive, prohibiting the wider adoption of existing first-principles methods. Here, we report an artificial-intelligence-enhanced protocol to drastically reduce the computational cost of simulating nonlinear time-resolved electronic spectra, which makes such simulations affordable for polyatomic molecules of increasing size. The protocol is based on the doorway-window approach for the on-the-fly surface-hopping simulations. We show its applicability for the prototypical molecule of pyrazine for which it produces spectra with high precision with respect to ab initio reference while cutting the computational cost by at least 95% compared to pure first-principles simulations.

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.

Ultrafast Photocontrolled Rotation in a Molecular Motor Investigated by Machine Learning-Based Nonadiabatic Dynamics Simulations

The thermal helix inversion (THI) of the overcrowded alkene-based molecular motors determines the speed of the unidirectional rotation due to the high reaction barrier in the ground state, in comparison with the ultrafast photoreaction process. Recently, a phosphine-based motor has achieved all-photochemical rotation experimentally, promising to be controlled without a thermal step. However, the mechanism of this photochemical reaction has not yet been fully revealed. The comprehensive computational studies on photoisomerization still resort to nonadiabatic molecular dynamics (NAMD) simulations based on electronic structure calculations, which remains a high computational cost for large systems such as molecular motors. Machine learning (ML) has become an accelerating tool in NAMD simulations recently, where excited-state potential energy surfaces (PESs) are constructed analytically with high accuracy, providing an efficient approach for simulations in photochemistry. Herein the reaction pathway is explored by a spin-flip time-dependent density functional theory (SF-TDDFT) approach in combination with ML-based NAMD simulations. According to our computational simulations, we notice that one of the key factors of fulfilling all-photochemical rotation in the phosphine-based motor is that the excitation energies of four isomers are similar. Additionally, a shortcut photoinduced transformation between unstable isomers replaces the THI step, which shares the conical intersection (CI) with photoisomerization. In this study, we provide a practical approach to speed up the NAMD simulations in photochemical reactions for a large system that could be extended to other complex systems.

Machine Learning Seams of Conical Intersection: A Characteristic Polynomial Approach

The machine learning of potential energy surfaces (PESs) has undergone rapid progress in recent years. The vast majority of this work, however, has been focused on the learning of ground state PESs. To reliably extend machine learning protocols to excited state PESs, the occurrence of seams of conical intersections between adiabatic electronic states must be correctly accounted for. This introduces a serious problem, for at such points, the adiabatic potentials are not differentiable to any order, complicating the application of standard machine learning methods. We show that this issue may be overcome by instead learning the coordinate-dependent coefficients of the characteristic polynomial of a simple decomposition of the potential matrix. We demonstrate that, through this approach, quantitatively accurate machine learning models of seams of conical intersection may be constructed.

Dynamics of Hydrogen Bond Breaking Induced by Outer-Valence Intermolecular Coulombic Decay

The photoexcitation of weakly bound complexes can lead to several decay pathways, depending on the nature of the potential energy surfaces. Upon excitation of a chromophore in a weakly bound complex, ionization of its neighbor upon energy transfer can occur due to a unique relaxation process known as intermolecular Coulombic decay (ICD), a phenomenon of renewed focus owing to its relevance in biological systems. Herein, we report the evidence for outer-valence ICD induced by multiphoton excitation by near-ultraviolet radiation of 4.4 eV photons, hitherto unknown in molecular systems. In the binary complexes of 2,6-difluorophenylacetylene with aliphatic amines, a resonant two-photon excitation localized on the 2,6-difluorophenylacetylene chromophore results in the formation of an amine cation following an outer-valence ICD process. The unique trends in experimentally observed translational energy distribution profiles of the amine cations following hydrogen bond dissociation, analyzed with the help of electronic structure and ab initio molecular dynamics calculations, revealed the presence of a delicate interplay of roaming dynamics, methyl-rotor dynamics, and binding energy.

Neural network potentials for chemistry: concepts, applications and prospects

Artificial Neural Networks (NN) are already heavily involved in methods and applications for frequent tasks in the field of computational chemistry such as representation of potential energy surfaces (PES) and spectroscopic predictions. This perspective provides an overview of the foundations of neural network-based full-dimensional potential energy surfaces, their architectures, underlying concepts, their representation and applications to chemical systems. Methods for data generation and training procedures for PES construction are discussed and means for error assessment and refinement through transfer learning are presented. A selection of recent results illustrates the latest improvements regarding accuracy of PES representations and system size limitations in dynamics simulations, but also NN application enabling direct prediction of physical results without dynamics simulations. The aim is to provide an overview for the current state-of-the-art NN approaches in computational chemistry and also to point out the current challenges in enhancing reliability and applicability of NN methods on a larger scale.

Photodegradation Kinetics and Deep Learning-Based Intelligent Colorimetric Method for Bioavailability-Based Dissolved Iron Speciation

Via the photodegradation of dissolved iron (dFe) complexes in the euphotic zone, released free Fe(III) is the most important source of bioavailable iron for eukaryotic phytoplankton. There is an urgent need to establish bioavailability-based dissolved iron speciation (BDIS) methods. Herein, an intelligent system with dFe pretreatment and a colorimetric sensor is developed for real-time monitoring of newly generated Fe(III) ions. According to the photodegradation kinetics of dFe, including kinetic constant and photogenerated time of free Fe(III) ions, 3 sources, 6 kinds, and 12 species of dFe are determined by our photocatalytic-assisted colorimetric sensor and deep learning model within 20.0 min. The algal dFe-uptake for 4 days can be predicted by BDIS with correlation coefficient 0.85, which could be explained by the hard and soft acids and bases theory (HSAB) and density functional theory (DFT). These results successfully demonstrate the proof-of-concept for photodegradation kinetics-based speciation and bioavailability assessments of dissolved metals.

A Look Inside the Black Box of Machine Learning Photodynamics Simulations

ConspectusPhotochemical reactions are of great importance in chemistry, biology, and materials science because they take advantage of a renewable energy source, mild reaction conditions, and high atom economy. Light absorption can excite molecules to a higher energy electronic state of the same spin multiplicity. The following nonadiabatic processes induce molecular transformations that afford exotic molecular architectures and high-energy-isomers that are inaccessible by thermal means. Computational simulations now complement time-resolved instrumentation to reveal ultrafast excited-state mechanistic information for photochemical reactions that is essential in disentangling elusive spectroscopic features, excited-state lifetimes, and excited-state mechanistic critical points. Nonadiabatic molecular dynamics (NAMD), powered by surface hopping techniques, is among the most widely applied techniques to model the photochemical reactions of medium-sized molecules. However, the computational efficiency is limited because of the requisite thousands of multiconfigurational quantum-chemical calculations multiplied by hundreds of trajectories. Machine learning (ML) has emerged as a revolutionary force in computational chemistry to predict the outcome of the resource-intensive multiconfigurational calculations on the fly. An ML potential trained with a substantial set of quantum-chemical calculations can predict the energies and forces with errors under chemical accuracy at a negligible cost. The integration of ML potentials in NAMD dramatically extends the maximum simulation time scale by ∼10 000-fold to the nanosecond regime.In this Account, we present a comprehensive demonstration of ML photodynamics simulations and summarize our most recent applications in resolving complex photochemical reactions. First, we address three fundamental components of ML techniques for photodynamics simulations: the quantum-chemical data set, the ML potential, and NAMD. Second, we describe best practices in building training data and our procedure toward training the ML photodynamics model with our recent literature contributions. We introduce a convenient training data generation scheme combining Wigner sampling and geometrical interpolation. It trains reliable and effective ML potentials suitable for subsequent active learning to detect undersampled data. We demonstrate how active learning automatically discovers new mechanistic pathways and reproduces experimental results. We point out that atomic permutation is an essential data augmentation approach to improve the learnability of distance-based molecular descriptors for highly symmetric molecules. Third, we demonstrate the utility of ML-photodynamics by showing the results of ML photodynamics simulations of (1) photo-torquoselective 4π disrotatory electrocyclic ring closing of norbornyl cyclohexadiene, which reveals a thermal conversion from experimentally unobserved intermediates to the reactant in 1 ns; (2) [2 + 2] photocycloaddition of substituted [3]-syn-ladderdienes in competition with 4π and 6π electrocyclic ring-opening reactions, uncovering substituent effects to explain the reported increased quantum yield of substituted cubane precursors; and (3) photochemical 4π disrotatory electrocyclic reactions of fluorobenzenes in nanoseconds with XMS-CASPT2-level training data. We expect this Account to broaden understanding of ML photodynamics and inspire future developments and applications to increasingly large molecules within complex environments on long time scales.