Advances in machine learned interatomic potentials (MLIPs), such as those using neural networks, have resulted in short-range models that can infer interaction energies with near ab initio accuracy and orders of magnitude reduced computational cost. For many atom systems, including macromolecules, biomolecules, and condensed matter, model accuracy can become reliant on the description of short- and long-range physical interactions. The latter terms can be difficult to incorporate into an MLIP framework. Recent research has produced numerous models with considerations for nonlocal electrostatic and dispersion interactions, leading to a large range of applications that can be addressed using MLIPs. In light of this, we present a Perspective focused on key methodologies and models being used where the presence of nonlocal physics and chemistry are crucial for describing system properties. The strategies covered include MLIPs augmented with dispersion corrections, electrostatics calculated with charges predicted from atomic environment descriptors, the use of self-consistency and message passing iterations to propagated nonlocal system information, and charges obtained via equilibration schemes. We aim to provide a pointed discussion to support the development of machine learning-based interatomic potentials for systems where contributions from only nearsighted terms are deficient.

Traditional computational approaches to design chemical species are limited by the need to compute properties for a vast number of candidates, e.g., by discriminative modeling. Therefore, inverse design methods aim to start from the desired property and optimize a corresponding chemical structure. From a machine learning viewpoint, the inverse design problem can be addressed through so-called generative modeling. Mathematically, discriminative models are defined by learning the probability distribution function of properties given the molecular or material structure. In contrast, a generative model seeks to exploit the joint probability of a chemical species with target characteristics. The overarching idea of generative modeling is to implement a system that produces novel compounds that are expected to have a desired set of chemical features, effectively sidestepping issues found in the forward design process. In this contribution, we overview and critically analyze popular generative algorithms like generative adversarial networks, variational autoencoders, flow, and diffusion models. We highlight key differences between each of the models, provide insights into recent success stories, and discuss outstanding challenges for realizing generative modeling discovered solutions in chemical applications.

A molecular dynamics study of water-soluble polymers: analysis of force fields from atomistic simulationsA force field (FF) analysis was performed on three water-soluble polymers (PAM, PNIPAAm, and PEO), that are of significant interest for biomedical applications, by measuring the polymer radius of gyration (Rg), solvent accessible surface area (SASA), radial distribution functions (g(r)), and relative shape anisotropy (κ 2 ) in dilute conditions. Three generalized FFs were utilized to model PAM and PNIPAAm (DREIDING, GAFF, and GAFF2) in conjunction with five water models: SPC, SPC/E, TIP3P, TIP4P, and TIP4P/2005. It was found that the DREIDING FF showed better agreement with PAM experimental data reported for Rg, irrespective of the water model. For PNIPAAm, the DREIDING FF was also the best performing among the FFs studied; however, the choice of water model played an important role in the predicted properties. Although the trends in SASA and κ 2 are consistent with Rg, g(r) of hydrophilic polymer atoms with oxygen of water molecules shows only slight changes with polymer-water force field combinations.For PEO modelled with the CHARMM C35r FF, all water models except TIP3P resulted in good agreement with experimental and previous simulation data. These results highlight the considerable impact that polymer and water force fields can have on sampling appropriate polymer conformations in solvated systems.

High-throughput calculations based on molecular simulations to predict the adsorption of molecules inside metal–organic frameworks (MOFs) have become a useful complement to experimental efforts to identify promising adsorbents for chemical separations and storage. For computational convenience, all existing efforts of this kind have relied on simulations in which the MOF is approximated as rigid. In this paper, we use extensive adsorption–relaxation simulations that fully include MOF flexibility effects to explore the validity of the rigid framework approximation. We also examine the accuracy of several approximate methods to incorporate framework flexibility that are more computationally efficient than adsorption–relaxation calculations. We first benchmark various models of MOF flexibility for four MOFs with well-established CO2 experimental consensus isotherms. We then consider a range of adsorption properties, including Henry’s constants, nondilute loadings, and adsorption selectivity, for seven adsorbates in 15 MOFs randomly selected from the CoRE MOF database. Our results indicate that in many MOFs adsorption–relaxation simulations are necessary to make quantitative predictions of adsorption, particularly for adsorption at dilute concentrations, although more standard calculations based on rigid structures can provide useful information. Finally, we investigate whether a correlation exists between the elastic properties of empty MOFs and the importance of including framework flexibility in making accurate predictions of molecular adsorption. Our results did not identify a simple correlation of this type.

The enormous number of combinations of adsorbing molecules and porous materials that exist is known as adsorption space. The adsorption space for microporous polymers has not yet been systematically explored, especially when compared with efforts for crystalline adsorbents. We report molecular simulation data for the adsorptive and structural properties of polymers of intrinsic microporosity with a diverse set of adsorbate species with 345 distinct adsorption isotherms and over 240,000 fresh and swollen structures. These structures and isotherms were obtained using a sorption-relaxation technique that accounts for the critical role of flexibility of the polymeric adsorbents. This enables us to introduce a set of correlations that can estimate adsorbent swelling and fractional free volume dilation as a function of adsorbate uptake based on readily characterized properties. The separation selectivity of the 276 distinct binary molecular pairs in our data is reported and high-performing adsorbent systems are identified.