{"title":"The evolution of the Amber additive protein force field: History, current status, and future.","authors":"Xianwei Wang, Danyang Xiong, Yueqing Zhang, Jihang Zhai, Yu-Cheng Gu, Xiao He","doi":"10.1063/5.0227517","DOIUrl":"https://doi.org/10.1063/5.0227517","url":null,"abstract":"<p><p>Molecular dynamics simulations are pivotal in elucidating the intricate properties of biological molecules. Nonetheless, the reliability of their outcomes hinges on the precision of the molecular force field utilized. In this perspective, we present a comprehensive review of the developmental trajectory of the Amber additive protein force field, delving into researchers' persistent quest for higher precision force fields and the prevailing challenges. We detail the parameterization process of the Amber protein force fields, emphasizing the specific improvements and retained features in each version compared to their predecessors. Furthermore, we discuss the challenges that current force fields encounter in balancing the interactions of protein-protein, protein-water, and water-water in molecular dynamics simulations, as well as potential solutions to overcome these issues.</p>","PeriodicalId":15313,"journal":{"name":"Journal of Chemical Physics","volume":"162 3","pages":""},"PeriodicalIF":3.1,"publicationDate":"2025-01-21","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143006055","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"化学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Kohn-Sham inversion for open-shell systems.","authors":"Jannis Erhard, Egor Trushin, Andreas Görling","doi":"10.1063/5.0239422","DOIUrl":"https://doi.org/10.1063/5.0239422","url":null,"abstract":"<p><p>Methods based on density-functional theory usually treat open-shell atoms and molecules within the spin-unrestricted Kohn-Sham (KS) formalism, which breaks symmetries in real and spin space. Symmetry breaking is possible because the KS Hamiltonian operator does not need to exhibit the full symmetry of the physical Hamiltonian operator, but only the symmetry of the spin density, which is generally lower. Symmetry breaking leads to spin contamination and prevents a proper classification of the KS wave function with respect to the symmetries of the physical electron system. Formally well-justified variants of the KS formalism that restore symmetries in real space, in spin space, or in both have been introduced long ago, but have rarely been used in practice. Here, we introduce numerically stable KS inversion methods to construct reference KS potentials from reference spin-densities for all four possibilities to treat open shell systems, non-symmetrized, spin-symmetrized, space-symmetrized, and fully-symmetrized. The reference spin-densities are obtained by full configuration interaction and high-level coupled cluster methods for the considered atoms and diatomic molecules. The decomposition of the total energy in contributions such as the non-interacting kinetic, the exchange, and the correlation energy is different in the four KS formalisms. Reference values for these differences are provided for the considered atoms and molecules. All KS inversions, except the fully symmetrized one, lead in some cases to solutions violating the Aufbau principle. In the purely spin-symmetrized KS formalism, this represents a violation of the KS v-representability condition, i.e., no proper KS wave functions exist in those cases.</p>","PeriodicalId":15313,"journal":{"name":"Journal of Chemical Physics","volume":"162 3","pages":""},"PeriodicalIF":3.1,"publicationDate":"2025-01-21","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143006131","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"化学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Magnus Bukhave Johansen, Andreas Erbs Hillers-Bendtsen, Hector H Corzo, Ashleigh Barnes, Kurt V Mikkelsen, Dmytro Bykov
{"title":"The divide expand consolidate scheme for unrestricted second order Møller-Plesset perturbation theory ground state energies.","authors":"Magnus Bukhave Johansen, Andreas Erbs Hillers-Bendtsen, Hector H Corzo, Ashleigh Barnes, Kurt V Mikkelsen, Dmytro Bykov","doi":"10.1063/5.0228963","DOIUrl":"https://doi.org/10.1063/5.0228963","url":null,"abstract":"<p><p>The linear scaling divide-expand-consolidate (DEC) framework is expanded to include unrestricted Hartree-Fock references. By partitioning the orbital space and employing local molecular orbitals, the full molecular calculation can be performed as independent calculations on individual fragments, making the method well-suited for massively parallel implementations. This approach also incorporates error control through the fragment optimization threshold (FOT), which maintains precision and consistency throughout the calculations. A benchmark was conducted for correlation energies of open-shell systems and the relative energies of both open- and closed-shell molecules at the MP2 level of theory. The full calculation result is achieved as the FOT approaches zero. For correlation energies, an FOT of 10-3 is sufficient to recover over 98% of the full result in all cases. However, for relative energies and the electronic energy component of oxidation potentials, a tighter FOT of 10-4 is required to keep the DEC error within 10% for both open- and closed-shell molecules. This is likely due to a lack of systematic error cancellation for the molecules with vastly different chemical natures. Therefore, for accurate relative energies, the FOT should be an order of magnitude lower, and additional caution is needed, particularly for large systems. The DEC method extension to unrestricted references maintains favorable features of linear scaling and can be implemented in a massively parallel algorithm to calculate correlation energies for large open-shell systems.</p>","PeriodicalId":15313,"journal":{"name":"Journal of Chemical Physics","volume":"162 3","pages":""},"PeriodicalIF":3.1,"publicationDate":"2025-01-21","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143006054","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"化学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Ab initio study on the dynamics and spectroscopy of collective rovibrational polaritons.","authors":"Tamás Szidarovszky","doi":"10.1063/5.0244977","DOIUrl":"https://doi.org/10.1063/5.0244977","url":null,"abstract":"<p><p>Accurate rovibrational molecular models are employed to gain insight in high-resolution into the collective effects and intermolecular processes arising when molecules in the gas phase interact with a resonant infrared (IR) radiation mode. An efficient theoretical approach is detailed, and numerical results are presented for the HCl, H2O, and CH4 molecules confined in an IR cavity. It is shown that by employing a rotationally resolved model for the molecules, revealing the various cavity-mediated interactions between the field-free molecular eigenstates, it is possible to obtain a detailed understanding of the physical processes governing the energy level structure, absorption spectra, and dynamic behavior of the confined systems. Collective effects, arising due to the cavity-mediated interaction between molecules, are identified in energy level shifts, in intensity borrowing effects in the absorption spectra, and in the intermolecular energy transfer occurring during Hermitian or non-Hermitian time propagation.</p>","PeriodicalId":15313,"journal":{"name":"Journal of Chemical Physics","volume":"162 3","pages":""},"PeriodicalIF":3.1,"publicationDate":"2025-01-21","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143006128","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"化学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Alessandro Simon, Luc Belloni, Daniel Borgis, Martin Oettel
{"title":"The orientational structure of a model patchy particle fluid: Simulations, integral equations, density functional theory, and machine learning.","authors":"Alessandro Simon, Luc Belloni, Daniel Borgis, Martin Oettel","doi":"10.1063/5.0248694","DOIUrl":"https://doi.org/10.1063/5.0248694","url":null,"abstract":"<p><p>We investigate the orientational properties of a homogeneous and inhomogeneous tetrahedral four-patch fluid (Bol-Kern-Frenkel model). Using integral equations, either (i) HNC or (ii) a modified HNC scheme with a simulation input, the full orientational dependence of pair and direct correlation functions is determined. Density functionals for the inhomogeneous problem are constructed via two different methods. The first, molecular density functional theory, utilizes the full direct correlation function and an isotropic hard-sphere bridge functional. The second method, a machine learning approach, uses a decomposition of the functional into an isotropic reference part and a mean-field orientational part, where both parts are improved by machine learning techniques. A comparison with the simulation data at hard walls and around hard tracers shows a similar performance of the two functionals. Machine learning strategies are discussed to eliminate residual differences, with the goal of obtaining machine-learning enhanced functionals for the general anisotropic fluid.</p>","PeriodicalId":15313,"journal":{"name":"Journal of Chemical Physics","volume":"162 3","pages":""},"PeriodicalIF":3.1,"publicationDate":"2025-01-21","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143006057","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"化学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Lixue Cheng, P Bernát Szabó, Zeno Schätzle, Derk P Kooi, Jonas Köhler, Klaas J H Giesbertz, Frank Noé, Jan Hermann, Paola Gori-Giorgi, Adam Foster
{"title":"Highly accurate real-space electron densities with neural networks.","authors":"Lixue Cheng, P Bernát Szabó, Zeno Schätzle, Derk P Kooi, Jonas Köhler, Klaas J H Giesbertz, Frank Noé, Jan Hermann, Paola Gori-Giorgi, Adam Foster","doi":"10.1063/5.0236919","DOIUrl":"https://doi.org/10.1063/5.0236919","url":null,"abstract":"<p><p>Variational ab initio methods in quantum chemistry stand out among other methods in providing direct access to the wave function. This allows, in principle, straightforward extraction of any other observable of interest, besides the energy, but, in practice, this extraction is often technically difficult and computationally impractical. Here, we consider the electron density as a central observable in quantum chemistry and introduce a novel method to obtain accurate densities from real-space many-electron wave functions by representing the density with a neural network that captures known asymptotic properties and is trained from the wave function by score matching and noise-contrastive estimation. We use variational quantum Monte Carlo with deep-learning Ansätze to obtain highly accurate wave functions free of basis set errors and from them, using our novel method, correspondingly accurate electron densities, which we demonstrate by calculating dipole moments, nuclear forces, contact densities, and other density-based properties.</p>","PeriodicalId":15313,"journal":{"name":"Journal of Chemical Physics","volume":"162 3","pages":""},"PeriodicalIF":3.1,"publicationDate":"2025-01-21","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143006129","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"化学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"pyVPT2: Interoperable software for anharmonic vibrational frequency calculations.","authors":"Philip M Nelson, C David Sherrill","doi":"10.1063/5.0251445","DOIUrl":"https://doi.org/10.1063/5.0251445","url":null,"abstract":"<p><p>We present pyVPT2, a program to perform second-order vibrational perturbation theory (VPT2) computations to obtain anharmonic vibrational frequencies. This program is written in Python and can utilize any of the several quantum chemistry programs that have been interfaced to the QCEngine project of the Molecular Sciences Software Institute (MolSSI). The requisite single point energy, gradient, or Hessian computations can be automatically performed in a distributed-parallel fashion by optionally using the MolSSI's QCFractal software. With pyVPT2, VPT2 anharmonic frequencies can now be computed using quantum chemistry programs that lack their own VPT2 capabilities.</p>","PeriodicalId":15313,"journal":{"name":"Journal of Chemical Physics","volume":"162 3","pages":""},"PeriodicalIF":3.1,"publicationDate":"2025-01-21","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143006053","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"化学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Mpemba effect in the relaxation of an active Brownian particle in a trap without metastable states.","authors":"Apurba Biswas, R Rajesh","doi":"10.1063/5.0246857","DOIUrl":"https://doi.org/10.1063/5.0246857","url":null,"abstract":"<p><p>We explore the role of activity in the occurrence of the Mpemba effect within a system of an active colloid diffusing in a potential landscape devoid of metastable minimum. The Mpemba effect is characterized by a phenomenon where a hotter system reaches equilibrium quicker than a colder one when both are rapidly cooled to the same low temperature. While a minimal asymmetry in the potential landscape is crucial for observing this effect in passive colloidal systems, the introduction of activity can either amplify or reduce the threshold of this minimal asymmetry, resulting in the activity-induced and suppressed Mpemba effect. We attribute these variations in the Mpemba effect to the effective translational shift in the phase boundaries, which occurs as activity is changed.</p>","PeriodicalId":15313,"journal":{"name":"Journal of Chemical Physics","volume":"162 3","pages":""},"PeriodicalIF":3.1,"publicationDate":"2025-01-21","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143006133","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"化学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Daniel Gómez Bustos, Sreeprasad Sreenivasan, Balazs Pinter
{"title":"A computational study on the effect of structural isomerism on the excited state lifetime and redox energetics of archetype iridium photoredox catalyst platforms [Ir(ppy)2(bpy)]+ and Ir(ppy)3.","authors":"Daniel Gómez Bustos, Sreeprasad Sreenivasan, Balazs Pinter","doi":"10.1063/5.0239293","DOIUrl":"https://doi.org/10.1063/5.0239293","url":null,"abstract":"<p><p>This study investigates the impact of structural isomerism on the excited state lifetime and redox energetics of heteroleptic [Ir(ppy)2(bpy)]+ and homoleptic Ir(ppy)3 photoredox catalysts using ground-state and time-dependent density functional theory methods. While the ground- and excited-state reduction potentials differ only slightly among the isomers of these complexes, our findings reveal significant variations in the radiative and non-radiative decay rates of the reactivity-controlling triplet 3MLCT states of these closely related species. The observed differences in radiative decay rates could be traced back to variations in the transition dipole moment, vertical energy gaps, and spin-orbit coupling of the isomers. In [Ir(ppy)2(bpy)]+, transition dipole moment differences play a significant role in controlling the relative lifetime of the triplet states, which we rationalized by a vectorial analysis of permanent dipole moments of the ground and excited states. Regarding the two isomers of Ir(ppy)3, changes in radiative decay rates were primarily attributed to variations in vertical energy gaps and intensity borrowing from other singlet-singlet transitions driven by spin-orbit coupling. Non-radiative decay variations were assessed in terms of differences in reorganization energies, adiabatic energy gap, and spin-orbit coupling. For both complexes, reorganization energies associated with low-energy molecular vibrations and metal-ligand bond length changes following the de-excitation process were major contributors. These insights provide a deeper understanding of how molecular design can be leveraged to optimize the performance of iridium-based photoredox catalysts, potentially guiding the development of more efficient catalytic systems for future applications.</p>","PeriodicalId":15313,"journal":{"name":"Journal of Chemical Physics","volume":"162 2","pages":""},"PeriodicalIF":3.1,"publicationDate":"2025-01-14","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142978569","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"化学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Darin E Mihalik, R Wang, B H Yang, P C Stancil, T J Price, R C Forrey, N Balakrishnan, R V Krems
{"title":"Accurate machine learning of rate coefficients for state-to-state transitions in molecular collisions.","authors":"Darin E Mihalik, R Wang, B H Yang, P C Stancil, T J Price, R C Forrey, N Balakrishnan, R V Krems","doi":"10.1063/5.0242182","DOIUrl":"https://doi.org/10.1063/5.0242182","url":null,"abstract":"<p><p>We present an algorithm that combines quantum scattering calculations with probabilistic machine-learning models to predict quantum dynamics rate coefficients for a large number of state-to-state transitions in molecule-molecule collisions much faster than with direct solutions of the Schrödinger equation. By utilizing the predictive power of Gaussian process regression with kernels, optimized to make accurate predictions outside of the input parameter space, the present strategy reduces the computational cost by about 75%, with an accuracy within 5%. Our method uses temperature dependences of rate coefficients for transitions from the isolated states of initial rotational angular momentum j, determined via explicit calculations, to predict the temperature dependences of rate coefficients for other values of j. The approach, demonstrated here for rovibrational transitions of SiO due to thermal collisions with H2, uses different prediction models and is thus adaptive to various time and accuracy requirements. The procedure outlined in this work can be used to extend multiple inelastic molecular collision databases without exponentially large computational resources required for conventional rigorous quantum dynamics calculations.</p>","PeriodicalId":15313,"journal":{"name":"Journal of Chemical Physics","volume":"162 2","pages":""},"PeriodicalIF":3.1,"publicationDate":"2025-01-14","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142978581","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"化学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}