Search

The relationship between the thermodynamic and computational properties of physical systems has been a major theoretical interest since at least the 19th century. It has also become of increasing practical importance over the last half-century as the energetic cost of digital devices has exploded. Importantly, real-world computers obey multiple physical constraints on how they work, which affects their thermodynamic properties. Moreover, many of these constraints apply to both naturally occurring computers, like brains or Eukaryotic cells, and digital systems. Most obviously, all such systems must finish their computation quickly, using as few degrees of freedom as possible. This means that they operate far from thermal equilibrium. Furthermore, manycomputers, both digital and biological, are modular, hierarchical systems with strong constraints on the connectivity among their subsystems. Yet another example is that to simplify their design, digital computers are required to be periodic processes governed by a global clock. None of these constraints were considered in 20th-century analyses of the thermodynamics of computation. The new field of stochastic thermodynamics provides formal tools for analyzing systems subject to all of these constraints. We argue here that these tools may help us understand at a far deeper level just how the fundamental thermodynamic properties of physical systems are related to the computation they perform. [ABSTRACT FROM AUTHOR]

The Faraday Discussions article discusses the general discussion on stochastic and low-scaling techniques for extended systems. Researchers explored various kernels and global descriptors to predict energy accurately. They found that global descriptors performed better for energy prediction, especially in homogeneous systems. The methodology was tested on hydrogen chains and showed promise for accurately representing wavefunctions and potential energy surfaces. The approach was found to be data-efficient and could potentially be applied to more complex systems in the future. [Extracted from the article]

Key to being able to accurately model the properties of realistic materials is being able to predict their properties in the thermodynamic limit. Nevertheless, because most many-body electronic structure methods scale as a high-order polynomial, or even exponentially, with system size, directly simulating large systems in their thermodynamic limit rapidly becomes computationally intractable. As a result, researchers typically estimate the properties of large systems that approach the thermodynamic limit by extrapolating the properties of smaller, computationally-accessible systems based on relatively simple scaling expressions. In this work, we employ Gaussian processes to more accurately and efficiently extrapolate many-body simulations to their thermodynamic limit. We train our Gaussian processes on Smooth Overlap of Atomic Positions (SOAP) descriptors to extrapolate the energies of one-dimensional hydrogen chains obtained using two high-accuracy many-body methods: coupled cluster theory and Auxiliary Field Quantum Monte Carlo (AFQMC). In so doing, we show that Gaussian processes trained on relatively short 10–30-atom chains can predict the energies of both homogeneous and inhomogeneous hydrogen chains in their thermodynamic limit with sub-milliHartree accuracy. Unlike standard scaling expressions, our GPR-based approach is highly generalizable given representative training data and is not dependent on systems' geometries or dimensionality. This work highlights the potential for machine learning to correct for the finite size effects that routinely complicate the interpretation of finite size many-body simulations. [ABSTRACT FROM AUTHOR]

The first magnetic 2D material discovered, monolayer (ML) CrI3, is particularly fascinating due to its ground state ferromagnetism. However, because ML materials are difficult to probe experimentally, much remains unresolved about ML CrI3's structural, electronic, and magnetic properties. Here, we leverage Density Functional Theory (DFT) and high-accuracy Diffusion Monte Carlo (DMC) simulations to predict lattice parameters, magnetic moments, and spin–phonon and spin–lattice coupling of ML CrI3. We exploit a recently developed surrogate Hessian DMC line search technique to determine CrI3's ML geometry with DMC accuracy, yielding lattice parameters in good agreement with recently published STM measurements—an accomplishment given the ∼10% variability in previous DFT-derived estimates depending upon the functional. Strikingly, we find that previous DFT predictions of ML CrI3's magnetic spin moments are correct on average across a unit cell but miss critical local spatial fluctuations in the spin density revealed by more accurate DMC. DMC predicts that magnetic moments in ML CrI3 are 3.62 μB per chromium and −0.145 μB per iodine, both larger than previous DFT predictions. The large disparate moments together with the large spin–orbit coupling of CrI3's I-p orbital suggest a ligand superexchange-dominated magnetic anisotropy in ML CrI3, corroborating recent observations of magnons in its 2D limit. We also find that ML CrI3 exhibits a substantial spin–phonon coupling of ∼3.32 cm−1. Our work, thus, establishes many of ML CrI3's key properties, while also continuing to demonstrate the pivotal role that DMC can assume in the study of magnetic and other 2D materials. [ABSTRACT FROM AUTHOR]

Spurred by recent technological advances, there is a growing demand for computational methods that can accurately predict the dynamics of correlated electrons. Such methods can provide much-needed theoretical insights into the electron dynamics probed via time-resolved spectroscopy experiments and observed in non-equilibrium ultracold atom experiments. In this article, we develop and benchmark a numerically exact Auxiliary Field Quantum Monte Carlo (AFQMC) method for modeling the dynamics of correlated electrons in real time. AFQMC has become a powerful method for predicting the ground state and finite temperature properties of strongly correlated systems mostly by employing constraints to control the sign problem. Our initial goal in this work is to determine how well AFQMC generalizes to real-time electron dynamics problems without constraints. By modeling the repulsive Hubbard model on different lattices and with differing initial electronic configurations, we show that real-time AFQMC is capable of accurately capturing long-lived electronic coherences beyond the reach of mean field techniques. While the times to which we can meaningfully model decrease with increasing correlation strength and system size as a result of the exponential growth of the dynamical phase problem, we show that our technique can model the short-time behavior of strongly correlated systems to very high accuracy. Crucially, we find that importance sampling, combined with a novel adaptive active space sampling technique, can substantially lengthen the times to which we can simulate. These results establish real-time AFQMC as a viable technique for modeling the dynamics of correlated electron systems and serve as a basis for future sampling advances that will further mitigate the dynamical phase problem. [ABSTRACT FROM AUTHOR]

Finite temperature auxiliary field-based quantum Monte Carlo methods, including determinant quantum Monte Carlo and Auxiliary Field Quantum Monte Carlo (AFQMC), have historically assumed pivotal roles in the investigation of the finite temperature phase diagrams of a wide variety of multidimensional lattice models and materials. Despite their utility, however, these techniques are typically formulated in the grand canonical ensemble, which makes them difficult to apply to condensates such as superfluids and difficult to benchmark against alternative methods that are formulated in the canonical ensemble. Working in the grand canonical ensemble is furthermore accompanied by the increased overhead associated with having to determine the chemical potentials that produce desired fillings. Given this backdrop, in this work, we present a new recursive approach for performing AFQMC simulations in the canonical ensemble that does not require knowledge of chemical potentials. To derive this approach, we exploit the convenient fact that AFQMC solves the many-body problem by decoupling many-body propagators into integrals over one-body problems to which non-interacting theories can be applied. We benchmark the accuracy of our technique on illustrative Bose and Fermi--Hubbard models and demonstrate that it can converge more quickly to the ground state than grand canonical AFQMC simulations. We believe that our novel use of HS-transformed operators to implement algorithms originally derived for non-interacting systems will motivate the development of a variety of other methods and anticipate that our technique will enable direct performance comparisons against other many-body approaches formulated in the canonical ensemble. [ABSTRACT FROM AUTHOR]

We review recent advances in the capabilities of the open source ab initio Quantum Monte Carlo (QMC) package QMCPACK and the workflow tool Nexus used for greater efficiency and reproducibility. The auxiliary field QMC (AFQMC) implementation has been greatly expanded to include k-point symmetries, tensor-hypercontraction, and accelerated graphical processing unit (GPU) support. These scaling and memory reductions greatly increase the number of orbitals that can practically be included in AFQMC calculations, increasing the accuracy. Advances in real space methods include techniques for accurate computation of bandgaps and for systematically improving the nodal surface of ground state wavefunctions. Results of these calculations can be used to validate application of more approximate electronic structure methods, including GW and density functional based techniques. To provide an improved foundation for these calculations, we utilize a new set of correlation-consistent effective core potentials (pseudopotentials) that are more accurate than previous sets; these can also be applied in quantum-chemical and other many-body applications, not only QMC. These advances increase the efficiency, accuracy, and range of properties that can be studied in both molecules and materials with QMC and QMCPACK. [ABSTRACT FROM AUTHOR]

The rate of modern drug discovery using experimental screening methods still lags behind the rate at which pathogens mutate, underscoring the need for fast and accurate predictive simulations of protein evolution. Multidrug-resistant bacteria evade our defenses by expressing a series of proteins, the most famous of which is the 29-kilodalton enzyme, TEM β-lactamase. Considering these challenges, we applied a covalent docking heuristic to measure the effects of all possible alanine 237 substitutions in TEM due to this codon's importance for catalysis and effects on the binding affinities of commercially-available β-lactam compounds. In addition to the usual mutations that reduce substrate binding due to steric hindrance, we identified two distinctive specificity-shifting TEM mutations, Ala237Arg and Ala237Lys, and their respective modes of action. Notably, we discovered and verified through minimum inhibitory concentration assays that, while these mutations and their bulkier side chains lead to steric clashes that curtail ampicillin binding, these same groups foster salt bridges with the negatively-charged side-chain of the cephalosporin cefixime, widely used in the clinic to treat multi-resistant bacterial infections. To measure the stability of these unexpected interactions, we used molecular dynamics simulations and found the binding modes to be stable despite the application of biasing forces. Finally, we found that both TEM mutants also bind strongly to other drugs containing negatively-charged R-groups, such as carumonam and ceftibuten. As with cefixime, this increased binding affinity stems from a salt bridge between the compounds' negative moieties and the positively-charged side chain of the arginine or lysine, suggesting a shared mechanism. In addition to reaffirming the power of using simulations as molecular microscopes, our results can guide the rational design of next-generation β-lactam antibiotics and bring the community closer to retaking the lead against the recurrent threat of multidrug-resistant pathogens. Author summary: Resistance to antibiotics is a major public health threat. Microorganisms are able to resist commonly used drugs by evolving and expressing enzymes capable of neutralizing antibiotics. Understanding the relationships between structural elements in these enzymes and their drug-clearing functions can lead to crucial insights for the discovery of next-generation antibiotics. In this study, we have used cutting-edge computational modeling methods to predict the effects of all naturally-occurring variations of an important region of the binding site of TEM β-lactamase, one of the major resistance-granting enzymes in bacteria. In an effort to identify patterns that could be useful for drug discovery, our simulations sought to understand how chemical changes in the tested region can affect resistance against a collection of over 90 widely used antibiotics. Crucially, through our simulations, we have identified a pathway for bacterial resistance against β-lactam antibiotics containing a negatively-charged moiety. We have also elucidated the mechanism behind the gain of resistance, which involves strong interactions between the drug's negative moieties and the positively-charged chemical shifts resulting from the mutation. Finally, we validated our predictions against fitness experiments for two commonly used antibiotics, which qualitatively corroborated our most important findings. [ABSTRACT FROM AUTHOR]

Identifying structural differences among proteins can be a non-trivial task. When contrasting ensembles of protein structures obtained from molecular dynamics simulations, biologically-relevant features can be easily overshadowed by spurious fluctuations. Here, we present SINATRA Pro, a computational pipeline designed to robustly identify topological differences between two sets of protein structures. Algorithmically, SINATRA Pro works by first taking in the 3D atomic coordinates for each protein snapshot and summarizing them according to their underlying topology. Statistically significant topological features are then projected back onto a user-selected representative protein structure, thus facilitating the visual identification of biophysical signatures of different protein ensembles. We assess the ability of SINATRA Pro to detect minute conformational changes in five independent protein systems of varying complexities. In all test cases, SINATRA Pro identifies known structural features that have been validated by previous experimental and computational studies, as well as novel features that are also likely to be biologically-relevant according to the literature. These results highlight SINATRA Pro as a promising method for facilitating the non-trivial task of pattern recognition in trajectories resulting from molecular dynamics simulations, with substantially increased resolution. Author summary: Structural features of proteins often serve as signatures of their biological function and molecular binding activity. Elucidating these structural features is essential for a full understanding of underlying biophysical mechanisms. While there are existing methods aimed at identifying structural differences between protein variants, such methods do not have the capability to jointly infer both geometric and dynamic changes, simultaneously. In this paper, we propose SINATRA Pro, a computational framework for extracting key structural features between two sets of proteins. SINATRA Pro robustly outperforms standard techniques in pinpointing the physical locations of both static and dynamic signatures across various types of protein ensembles, and it does so with improved resolution. [ABSTRACT FROM AUTHOR]

Summary Single amino acid variations (SAVs) are a primary contributor to variations in the human genome. Identifying pathogenic SAVs can provide insights to the genetic architecture of complex diseases. Most approaches for predicting the functional effects or pathogenicity of SAVs rely on either sequence or structural information. This study presents 〈Lai Yang Rubenstein Uzun Sarkar〉 (LYRUS), a machine learning method that uses an XGBoost classifier to predict the pathogenicity of SAVs. LYRUS incorporates five sequence-based, six structure-based and four dynamics-based features. Uniquely, LYRUS includes a newly proposed sequence co-evolution feature called the variation number. LYRUS was trained using a dataset that contains 4363 protein structures corresponding to 22 639 SAVs from the ClinVar database, and tested using the VariBench testing dataset. Performance analysis showed that LYRUS achieved comparable performance to current variant effect predictors. LYRUS's performance was also benchmarked against six Deep Mutational Scanning datasets for PTEN and TP53. Availability and implementation LYRUS is freely available and the source code can be found at https://github.com/jiaying2508/LYRUS. Supplementary information Supplementary data are available at Bioinformatics Advances online. [ABSTRACT FROM AUTHOR]

Data encoded in molecules offers opportunities for secret messaging and extreme information density. Here, we explore how the same chemical and physical dimensions used to encode molecular information can expose molecular messages to detection and manipulation. To address these vulnerabilities, we write data using an object's pre-existing surface chemistry in ways that are indistinguishable from the original substrate. While it is simple to embed chemical information onto common objects (covers) using routine steganographic permutation, chemically embedded covers are found to be resistant to detection by sophisticated analytical tools. Using Turbo codes for efficient digital error correction, we demonstrate recovery of secret keys hidden in the pre-existing chemistry of American one dollar bills. These demonstrations highlight ways to improve security in other molecular domains, and show how the chemical fingerprints of common objects can be harnessed for data storage and communication. [ABSTRACT FROM AUTHOR]

Chemical mixtures can be leveraged to store large amounts of data in a highly compact form and have the potential for massive scalability owing to the use of large-scale molecular libraries. With the parallelism that comes from having many species available, chemical-based memory can also provide the physical substrate for computation with increased throughput. Here, we represent non-binary matrices in chemical solutions and perform multiple matrix multiplications and additions, in parallel, using chemical reactions. As a case study, we demonstrate image processing, in which small greyscale images are encoded in chemical mixtures and kernel-based convolutions are performed using phenol acetylation reactions. In these experiments, we use the measured concentrations of reaction products (phenyl acetates) to reconstruct the output image. In addition, we establish the chemical criteria required to realize chemical image processing and validate reaction-based multiplication. Most importantly, this work shows that fundamental arithmetic operations can be reliably carried out with chemical reactions. Our approach could serve as a basis for developing more advanced chemical computing architectures. [ABSTRACT FROM AUTHOR]

Molecular data systems have the potential to store information at dramatically higher density than existing electronic media. Some of the first experimental demonstrations of this idea have used DNA, but nature also uses a wide diversity of smaller non-polymeric molecules to preserve, process, and transmit information. In this paper, we present a general framework for quantifying chemical memory, which is not limited to polymers and extends to mixtures of molecules of all types. We show that the theoretical limit for molecular information is two orders of magnitude denser by mass than DNA, although this comes with different practical constraints on total capacity. We experimentally demonstrate kilobyte-scale information storage in mixtures of small synthetic molecules, and we consider some of the new perspectives that will be necessary to harness the information capacity available from the vast non-genomic chemical space. [ABSTRACT FROM AUTHOR]

One of the long-standing holy grails of molecular evolution has been the ability to predict an organism's fitness directly from its genotype. With such predictive abilities in hand, researchers would be able to more accurately forecast how organisms will evolve and how proteins with novel functions could be engineered, leading to revolutionary advances in medicine and biotechnology. In this work, we assemble the largest reported set of experimental TEM-1 β-lactamase folding free energies and use this data in conjunction with previously acquired fitness data and computational free energy predictions to determine how much of the fitness of β-lactamase can be directly predicted by thermodynamic folding and binding free energies. We focus upon β-lactamase because of its long history as a model enzyme and its central role in antibiotic resistance. Based upon a set of 21 β-lactamase single and double mutants expressly designed to influence protein folding, we first demonstrate that modeling software designed to compute folding free energies such as FoldX and PyRosetta can meaningfully, although not perfectly, predict the experimental folding free energies of single mutants. Interestingly, while these techniques also yield sensible double mutant free energies, we show that they do so for the wrong physical reasons. We then go on to assess how well both experimental and computational folding free energies explain single mutant fitness. We find that folding free energies account for, at most, 24% of the variance in β-lactamase fitness values according to linear models and, somewhat surprisingly, complementing folding free energies with computationally-predicted binding free energies of residues near the active site only increases the folding-only figure by a few percent. This strongly suggests that the majority of β-lactamase's fitness is controlled by factors other than free energies. Overall, our results shed a bright light on to what extent the community is justified in using thermodynamic measures to infer protein fitness as well as how applicable modern computational techniques for predicting free energies will be to the large data sets of multiply-mutated proteins forthcoming. [ABSTRACT FROM AUTHOR]

Biomolecular information systems offer exciting potential advantages and opportunities to complement conventional semiconductor technologies. Much attention has been paid to information-encoding polymers, but small molecules also play important roles in biochemical information systems. Downstream from DNA, the metabolome is an information-rich molecular system with diverse chemical dimensions which could be harnessed for information storage and processing. As a proof of principle of small-molecule postgenomic data storage, here we demonstrate a workflow for representing abstract data in synthetic mixtures of metabolites. Our approach leverages robotic liquid handling for writing digital information into chemical mixtures, and mass spectrometry for extracting the data. We present several kilobyte-scale image datasets stored in synthetic metabolomes, which can be decoded with accuracy exceeding 99% using multi-mass logistic regression. Cumulatively, >100,000 bits of digital image data was written into metabolomes. These early demonstrations provide insight into some of the benefits and limitations of small-molecule chemical information systems. [ABSTRACT FROM AUTHOR]

The extent of coupling between the folding of a protein and its binding to a substrate varies from protein to protein. Some proteins have highly structured native states in solution, while others are natively disordered and only fold fully upon binding. In this Letter, we use Monte Carlo simulations to investigate how disordered polymer chains grafted around a binding site affect the folding and binding of three model proteins. The protein that approaches the substrate fully folded is more hindered during the binding process than those whose folding and binding are cooperative. The polymer chains act as localized crowding agents and can select correctly folded and bound configurations in favor of nonspecifically adsorbed states. The free energy change for forming all intraprotein and protein-substrate contacts can depend nonmonotonically on the polymer length. [ABSTRACT FROM AUTHOR]

The first magnetic 2D material discovered, monolayer (ML) CrI3, is particularly fascinating due to its ground state ferromagnetism. However, because ML materials are difficult to probe experimentally, much remains unresolved about ML CrI3's structural, electronic, and magnetic properties. Here, we leverage Density Functional Theory (DFT) and high-accuracy Diffusion Monte Carlo (DMC) simulations to predict lattice parameters, magnetic moments, and spin–phonon and spin–lattice coupling of ML CrI3. We exploit a recently developed surrogate Hessian DMC line search technique to determine CrI3's ML geometry with DMC accuracy, yielding lattice parameters in good agreement with recently published STM measurements—an accomplishment given the ∼10% variability in previous DFT-derived estimates depending upon the functional. Strikingly, we find that previous DFT predictions of ML CrI3's magnetic spin moments are correct on average across a unit cell but miss critical local spatial fluctuations in the spin density revealed by more accurate DMC. DMC predicts that magnetic moments in ML CrI3 are 3.62 μB per chromium and −0.145 μB per iodine, both larger than previous DFT predictions. The large disparate moments together with the large spin–orbit coupling of CrI3's I-p orbital suggest a ligand superexchange-dominated magnetic anisotropy in ML CrI3, corroborating recent observations of magnons in its 2D limit. We also find that ML CrI3 exhibits a substantial spin–phonon coupling of ∼3.32 cm−1. Our work, thus, establishes many of ML CrI3's key properties, while also continuing to demonstrate the pivotal role that DMC can assume in the study of magnetic and other 2D materials. [ABSTRACT FROM AUTHOR]

Quantum Monte Carlo (QMC) algorithms have long relied on Jastrow factors to incorporate dynamic correlation into trial wave functions. While Jastrow-type wave functions have been widely employed in real-space algorithms, they have seen limited use in second-quantized QMC methods, particularly in projection methods that involve a stochastic evolution of the wave function in imaginary time. Here we propose a scheme for generating Jastrow-type correlated trial wave functions for auxiliary-field QMC methods. The method is based on decoupling the two-body Jastrow into one-body projectors coupled to auxiliary fields, which then operate on a single determinant to produce a multideterminant trial wave function. We demonstrate that intelligent sampling of the most significant determinants in this expansion can produce compact trial wave functions that reduce errors in the calculated energies. Our technique may be readily generalized to accommodate a wide range of two-body Jastrow factors and applied to a variety of model and chemical systems. [ABSTRACT FROM AUTHOR]

We report the observation of a π-type dipole-bound state (π-DBS) in cryogenically cooled deprotonated 9-anthrol molecular anions (9AT-) by resonant two-photon photoelectron imaging. A DBS is observed 191 cm-1 (0.0237 eV) below the detachment threshold, and the existence of the π-DBS is revealed by a distinct (s+d)-wave photoelectron angular distribution. The π-DBS is stabilized by the large anisotropic in-plane polarizability of 9AT. The population of the dipole-forbidden π-DBS is proposed to be via a nonadiabatic coupling with the dipole-allowed σ-type DBS mediated by molecular rotations. [ABSTRACT FROM AUTHOR]

The accurate mapping of potential energy surfaces (PESs) is crucial to our understanding of the numerous physical and chemical processes mediated by atomic rearrangements, such as conformational changes and chemical reactions, and the thermodynamic and kinetic feasibility of these processes. Stochastic electronic structure theories, e.g., Quantum Monte Carlo (QMC) methods, enable highly accurate total energy calculations that in principle can be used to construct the PES. However, their stochastic nature poses a challenge to the computation and use of forces and Hessians, which are typically required in algorithms for minimum-energy pathway (MEP) and transition state (TS) identification, such as the nudged elastic band (NEB) algorithm and its climbing image formulation. Here, we present strategies that utilize the surrogate Hessian line-search method, previously developed for QMC structural optimization, to efficiently identify MEP and TS structures without requiring force calculations at the level of the stochastic electronic structure theory. By modifying the surrogate Hessian algorithm to operate in path-orthogonal subspaces and at saddle points, we show that it is possible to identify MEPs and TSs by using a force-free QMC approach. We demonstrate these strategies via two examples, the inversion of the ammonia (NH 3 ) molecule and the nucleophilic substitution (S N 2) reaction F - + CH 3 F → FCH 3 + F - . We validate our results using Density Functional Theory (DFT)- and Coupled Cluster (CCSD, CCSD(T))-based NEB calculations. We then introduce a hybrid DFT-QMC approach to compute thermodynamic and kinetic quantities, free energy differences, rate constants, and equilibrium constants that incorporates stochastically optimized structures and their energies, and show that this scheme improves upon DFT accuracy. Our methods generalize straightforwardly to other systems and other high-accuracy theories that similarly face challenges computing energy gradients, paving the way for highly accurate PES mapping, transition state determination, and thermodynamic and kinetic calculations at significantly reduced computational expense.