Journal of Chemical Theory and Computation

We present a generalization of the connectivity-based hierarchy (CBH) of isodesmic-based correction schemes to a multilayered fragmentation platform for overall cost reduction while retaining high accuracy. The newly developed multilayered CBH approach, called stepping-stone CBH (SSCBH), is benchmarked on a diverse set of 959 medium-sized organic molecules. Applying SSCBH corrections to the PBEh-D3 density functional resulted in an average error of 0.76 kcal/mol for the full test set compared to accurate CCSD(T)-quality enthalpies and an even lower error of 0...

Efficient ways to prepare Fermionic ground states on quantum computers are in high demand, and different techniques have been developed over the past few years. Despite having a vast set of methods, it is still unclear which method performs well for which system. In this work, we combine interpretable circuit designs with an effective basis approach in order to optimize a multiconfigurational valence bond wave function. Based on selected model systems, we show how this leads to an explainable performance. We demonstrate that the developed methodology outperforms related methods in terms of the size of the effective basis, as well as individual quantum resources for the involved circuits...

We present a linear scaling atomic orbital based algorithm for the computation of the most expensive exchange-type RI-MP2-F12 term by employing numerical quadrature in combination with CABS-RI to avoid six-center-three-electron integrals. Furthermore, a robust distance-dependent integral screening scheme, based on integral partition bounds [Thompson, T. H.; Ochsenfeld, C. J. Chem. Phys. 2019, 150, 044101], is used to drastically reduce the number of the required three-center-one-electron integrals substantially...

The projected atomic orbital (PAO) technique is presented for the construction of virtual orbital spaces in projection-based embedding (PbE) applications. The proposed straightforward procedure produces a set of virtual orbitals that are used in the final, high-level calculation of the embedded active subsystem. The PAO scheme is demonstrated on intermolecular potentials of bimolecular complexes in ground and excited states, including Rydberg excitations. The results show the outstanding performance of the PbE method when used with PAO virtual orbitals compared with those produced using common orbital localization techniques...

We present an efficient analytical energy gradient algorithm for the cluster-in-molecule resolution-of-identity second-order Møller-Plesset perturbation (CIM-RI-MP2) method based on the Lagrange multiplier method. Our algorithm independently constructs the Lagrangian formalism within each cluster, avoiding the solution of the coupled-perturbed Hartree-Fock (CPHF) equation for the whole system. Due to this feature, the computational cost of the CIM-RI-MP2 gradients is much lower than that of other local MP2 algorithms...

A specific checkpoint between target DNA binding and cleavage primarily governs the precision of Cas9 gene editing. Although various CRISPR-Cas9 variants have been developed to improve DNA cleavage accuracy, we still lack a comprehensive understanding of how they work at the molecular level. Herein, we have focused on studying the late-stage conformational transitions of Cas9 and an evolved Cas9 mutant (evoCas9) that start from the precleavage state. Our submilliseconds of dynamic simulations reveal that the presence of base mismatches leads the HNH nuclease domain of Cas9 to alter its principal functional modes of motion, thereby impairing its conformational activation...

Identifying optimal collective variables to model transformations using atomic-scale simulations is a long-standing challenge. We propose a new method for the generation, optimization, and comparison of collective variables that can be thought of as a data-driven generalization of the path collective variable concept. It consists of a kernel ridge regression of the committor probability, which encodes a transformation's progress. The resulting collective variable is one-dimensional, interpretable, and differentiable, making it appropriate for enhanced sampling simulations requiring biasing...

It has been experimentally demonstrated that molecular-vibration polaritons formed by strong coupling of a molecular vibration to an infrared cavity mode can significantly modify the physical properties and chemical reactivities of various molecular systems. However, a complete theoretical understanding of the underlying mechanisms of the modifications remains elusive due to the complexity of the hybrid system, especially the collective nature of polaritonic states in systems containing many molecules. We develop here the semiclassical theory of molecular vibration-polariton dynamics based on the truncated Wigner approximation (TWA) that is tractable in large molecular systems and simultaneously captures the quantum character of photons in the optical cavity...

State preparation for quantum algorithms is crucial for achieving high accuracy in quantum chemistry and competing with classical algorithms. The localized active space-unitary coupled cluster (LAS-UCC) algorithm iteratively loads a fragment-based multireference wave function onto a quantum computer. In this study, we compare two state preparation methods, quantum phase estimation (QPE) and direct initialization (DI), for each fragment. We test the two state preparation methods on three systems, ranging from a model system, a set of interacting hydrogen molecules, to more realistic chemical problems, like the C-C double bond breaking in transbutadiene and the spin ladder in a bimetallic system...

The photochemistry of nitrophenols is a source of smog as nitrous acid is formed from their photolysis. Nevertheless, computational studies of the photochemistry of these widespread toxic molecules are scarce. In this work, the initial photodeactivation of ortho -nitrophenol and para -nitrophenol is modeled, both in gas phase and in aqueous solution to simulate atmospheric and aerosol environments. A large number of excited states, six for ortho -nitrophenol and 11 for para -nitrophenol, have been included and were all populated during the decay...

Over the years, Hedin's GW self-energy has been proven to be a rather accurate and simple approximation to evaluate electronic quasiparticle energies in solids and in molecules. Attempts to improve over the simple GW approximation, the so-called vertex corrections, have been constantly proposed in the literature. Here, we derive, analyze, and benchmark the complete second-order term in the screened Coulomb interaction W for finite systems. This self-energy named G 3 W 2 contains all the possible time orderings that combine 3 Green's functions G and 2 dynamic W ...

On the one hand, much of computational chemistry is concerned with "bottom-up" calculations which elucidate observable behavior starting from exact or approximated physical laws, a paradigm exemplified by typical quantum mechanical calculations and molecular dynamics simulations. On the other hand, "top down" computations aiming to formulate mathematical models consistent with observed data, e.g., parametrizing force fields, binding or kinetic models, have been of interest for decades but recently have grown in sophistication with the use of Bayesian inference (BI)...

Memory effects emerge as a fundamental consequence of dimensionality reduction when low-dimensional observables are used to describe the dynamics of complex many-body systems. In the context of molecular dynamics (MD) data analysis, accounting for memory effects using the framework of the generalized Langevin equation (GLE) has proven efficient, accurate, and insightful, particularly when working with high-resolution time series data. However, in experimental systems, high-resolution data are often unavailable, raising questions about the impact of the data resolution on the estimated GLE parameters...

The Protein Structure Transformer (PeSTo), a geometric transformer, has exhibited exceptional performance in predicting protein-protein binding interfaces and distinguishing interfaces with nucleic acids, lipids, small molecules, and ions. In this study, we introduce PeSTo-Carbs, an extension of PeSTo specifically engineered to predict protein-carbohydrate binding interfaces. We evaluate the performance of this approach using independent test sets and compare them with those of previous methods. Furthermore, we highlight the model's capability to specialize in predicting interfaces involving cyclodextrins, a biologically and pharmaceutically significant class of carbohydrates...

We show that chemically accurate potential energy surfaces (PESs) can be generated from quantum computers by measuring only the density along an adiabatic transition between different molecular geometries. In lieu of using phase estimation, the energy is evaluated by performing line-integration using the inverted real-space time-dependent density functional theory Kohn-Sham (KS) potential obtained from the geometry-varying densities of the full wave function. The accuracy of this method depends on the validity of the adiabatic evolution itself and the potential inversion process (which is theoretically exact but can be numerically unstable), whereas the total evolution time is the defining factor for the precision of phase estimation...

We propose a state-averaged orbital optimization scheme for improving the accuracy of excited states of the electronic structure Hamiltonian for use on near-term quantum computers. Instead of parameterizing the orbital rotation operator in the conventional fashion as an exponential of an antihermitian matrix, we parameterize the orbital rotation as a general partial unitary matrix. Whereas conventional orbital optimization methods minimize the state-averaged energy using successive Newton steps of the second-order Taylor expansion of the energy, the method presented here optimizes the state-averaged energy using an orthogonally constrained gradient projection method that does not require any expansion approximations...

A major obstacle hindering the application of orbital-free density functional theory (OF-DFT) to all metals is the lack of accurate local pseudopotentials (LPSs), especially for transition metals. In this work, we developed high-quality LPSs for all simple and transition metals by fitting the atomic eigenvalues and orbital norms beyond the cutoff radii. Due to the lack of nonlocality in LPSs, it is very challenging to simultaneously fit the semicore and outermost valence orbitals of transition metals. We overcame this issue by excluding the semicore orbitals from the LPS optimizations...

We present low-scaling algorithms for GW and constrained random phase approximation based on a symmetry-adapted interpolative separable density fitting (ISDF) procedure that incorporates the space-group symmetries of crystalline systems. The resulting formulations scale cubically, with respect to system size, and linearly with the number of k -points, regardless of the choice of single-particle basis and whether a quasiparticle approximation is employed. We validate these methods through comparisons with published literature and demonstrate their efficiency in treating large-scale systems through the construction of downfolded many-body Hamiltonians for carbon dimer defects embedded in hexagonal boron nitride supercells...

The human L-type amino acid transporter 1 (LAT1; SLC7A5), is an amino acid exchanger protein, primarily found in the blood-brain barrier, placenta, and testis, where it plays a key role in amino acid homeostasis. Cholesterol is an essential lipid that has been highlighted to play a role in regulating the activity of membrane transporters, such as LAT1, yet little is known about the molecular mechanisms driving this phenomenon. Here we perform a comprehensive computational analysis to investigate cholesterol's role in LAT1 structure and function, focusing on four cholesterol-binding sites (CHOL1-4) identified in a recent LAT1-apo inward-open conformation cryo-EM structure...

The characterization of solvation shells of atoms, ions, and molecules in solution is essential to relate solvation properties to chemical phenomena such as complex formation and reactivity. Different definitions of the first-shell coordination sphere from simulation data can lead to potentially conflicting data on the structural properties and associated ligand exchange dynamics. The definition of a solvation shell is typically based on a given threshold distance determined from the respective solute-solvent pair distribution function g ( r ) (i...