An official website of the United States government
Free energy evaluation in molecular simulations of both classical and quantum systems is computationally intensive and requires sophisticated algorithms. This is because free energy depends on the volume of accessible phase space, a quantity that is inextricably linked to the integration measure in a coordinate representation of a many-body problem. In contrast, the same problem expressed as a field theory (auxiliary field or coherent states) isolates the particle number as a simple parameter in the Hamiltonian or action functional and enables the identification of a chemical potential field operator. We show that this feature leads a “direct” method of free energy evaluation, in which a particle model is converted to a field theory and appropriate field operators are averaged using a field-theoretic simulation conducted with complex Langevin sampling. These averages provide an immediate estimate of the Helmholtz free energy in the canonical ensemble and the entropy in the microcanonical ensemble. The method is illustrated for a classical polymer solution, a block copolymer melt exhibiting liquid crystalline and solid mesophases, and a quantum fluid of interacting bosons.
more »
« less
This content will become publicly available on December 28, 2025
Preserving positivity in density-explicit field-theoretic simulations
Field-theoretic simulations are numerical methods for polymer field theory, which include fluctuation corrections beyond the mean-field level, successfully capturing various mesoscopic phenomena. Most field-theoretic simulations of polymeric fluids use the auxiliary field (AF) theory framework, which employs Hubbard–Stratonovich transformations for the particle-to-field conversion. Nonetheless, the Hubbard–Stratonovich transformation imposes significant limitations on the functional form of the non-bonded potentials. Removing this restriction on the non-bonded potentials will enable studies of a wide range of systems that require multi-body or more complex potentials. An alternative representation is the hybrid density-explicit auxiliary field theory (DE-AF), which retains both a density field and a conjugate auxiliary field for each species. While the DE-AF representation is not new, density-explicit field-theoretic simulations have yet to be developed. A major challenge is preserving the real and non-negative nature of the density field during stochastic evolution. To address this, we introduce positivity-preserving schemes that enable the first stable and efficient density-explicit field-theoretic simulations (DE-AF FTS). By applying the new method to a simple fluid, we find thermodynamically correct results at high densities, but the algorithm fails in the dilute regime. Nonetheless, DE-AF FTS is shown to be broadly applicable to dense fluid systems including a simple fluid with a three-body non-bonded potential, a homopolymer solution, and a diblock copolymer melt.
more »
« less
- Award ID(s):
- 2104255
- PAR ID:
- 10633198
- Publisher / Repository:
- American Institute of Physics
- Date Published:
- Journal Name:
- The Journal of Chemical Physics
- Volume:
- 161
- Issue:
- 24
- ISSN:
- 0021-9606
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
null (Ed.)Dinitrogen pentoxide (N2O5) is an important intermediate in the atmospheric chemistry of nitrogen oxides. Although there has been much research, the processes that govern the physical interactions between N2O5 and water are still not fully understood at a molecular level. Gaining a quantitative insight from computer simulations requires going beyond the accuracy of classical force fields while accessing length scales and time scales that are out of reach for high-level quantum-chemical approaches. To this end, we present the development of MB-nrg many-body potential energy functions for nonreactive simulations of N2O5 in water. This MB-nrg model is based on electronic structure calculations at the coupled cluster level of theory and is compatible with the successful MB-pol model for water. It provides a physically correct description of long-range many-body interactions in combination with an explicit representation of up to three-body short-range interactions in terms of multidimensional permutationally invariant polynomials. In order to further investigate the importance of the underlying interactions in the model, a TTM-nrg model was also devised. TTM-nrg is a more simplistic representation that contains only two-body short-range interactions represented through Born–Mayer functions. In this work, an active learning approach was employed to efficiently build representative training sets of monomer, dimer, and trimer structures, and benchmarks are presented to determine the accuracy of our new models in comparison to a range of density functional theory methods. By assessing the binding curves, distortion energies of N2O5, and interaction energies in clusters of N2O5 and water, we evaluate the importance of two-body and three-body short-range potentials. The results demonstrate that our MB-nrg model has high accuracy with respect to the coupled cluster reference, outperforms current density functional theory models, and thus enables highly accurate simulations of N2O5 in aqueous environments.more » « less
-
Density functional theory (DFT) has been applied to modeling molecular interactions in water for over three decades. The ubiquity of water in chemical and biological processes demands a unified understanding of its physics, from the single molecule to the thermodynamic limit and everything in between. Recent advances in the development of data-driven and machine-learning potentials have accelerated simulation of water and aqueous systems with DFT accuracy. However, anomalous properties of water in the condensed phase, where a rigorous treatment of both local and non-local many-body (MB) interactions is in order, are often unsatisfactory or partially missing in DFT models of water. In this review, we discuss the modeling of water and aqueous systems based on DFT and provide a comprehensive description of a general theoretical/computational framework for the development of data-driven many-body potentials from DFT reference data. This framework, coined MB-DFT, readily enables efficient many-body molecular dynamics (MD) simulations of small molecules, in both gas and condensed phases, while preserving the accuracy of the underlying DFT model. Theoretical considerations are emphasized, including the role that the delocalization error plays in MB-DFT potentials of water and the possibility to elevate DFT and MB-DFT to near-chemical-accuracy through a density-corrected formalism. The development of the MB-DFT framework is described in detail, along with its application in MB-MD simulations and recent extension to the modeling of reactive processes in solution within a quantum mechanics/MB molecular mechanics (QM/MB-MM) scheme, using water as a prototypical solvent. Finally, we identify open challenges and discuss future directions for MB-DFT and QM/MB-MM simulations in condensed phases.more » « less
-
A<sc>bstract</sc> We investigate finite-temperature observables in three-dimensional largeNcritical vector models taking into account the effects suppressed by$$ \frac{1}{N} $$ . Such subleading contributions are captured by the fluctuations of the Hubbard-Stratonovich auxiliary field which need to be handled with care due to a subtle divergence structure which we clarify. The examples we consider include the scalarO(N) model, the Gross-Neveu model, the Nambu-Jona-Lasinio model and the massless Chern-Simons Quantum Electrodynamics. We present explicit results for the free energy density to the subleading order in$$ \frac{1}{N} $$ , which captures the thermal one-point function of the stress-energy tensor to this order. We also include the dependence on a chemical potential. We determine the Wilson coefficient in the thermal effective action that is sensitive to global symmetry for the first time directly in interacting CFTs, which produces a symmetry-resolved asymptotic density of states. We further provide a formula from diagrammatics for the one-point functions of general single-trace higher-spin currents. We observe that in most cases considered, these subleading effects lift the apparent degeneracies between observables in different models at infiniteN, while in special cases the discrepancies only start to appear at the next-to-subleading order.more » « less
-
Simulations of anharmonic vibrational motion rely on computationally expedient representations of the governing potential energy surface. The n-mode representation (n-MR)—effectively a many-body expansion in the space of molecular vibrations—is a general and efficient approach that is often used for this purpose in vibrational self-consistent field (VSCF) calculations and correlated analogues thereof. In the present analysis, a lack of convergence in many VSCF calculations is shown to originate from negative and unbound potentials at truncated orders of the n-MR expansion. For cases of strong anharmonic coupling between modes, the n-MR can both dip below the true global minimum of the potential surface and lead to effective single-mode potentials in VSCF that do not correspond to bound vibrational problems, even for bound total potentials. The present analysis serves mainly as a pathology report of this issue. Furthermore, this insight into the origin of VSCF non-convergence provides a simple, albeit ad hoc, route to correct the problem by “painting in” the full representation of groups of modes that exhibit these negative potentials at little additional computational cost. Somewhat surprisingly, this approach also reasonably approximates the results of the next-higher n-MR order and identifies groups of modes with particularly strong coupling. The method is shown to identify and correct problematic triples of modes—and restore SCF convergence—in two-mode representations of challenging test systems, including the water dimer and trimer, as well as protonated tropine.more » « less
