Search

FEASST is an open-source Monte Carlo software package for particle-based simulations. This software, which was released in 2017, has been used to study phase equilibrium, self-assembly, aggregation or gelation in biological materials, colloids, polymers, ionic liquids, and adsorption in porous networks. We highlight some of the unique features available in FEASST, such as flat-histogram grand canonical ensemble, Gibbs ensemble, and Mayer-sampling simulations with support for anisotropic models and parallelization with flat-histogram and prefetching. We also discuss how the challenges of supporting a variety of Monte Carlo algorithms were overcome by an object-oriented design. This also allows others to extend classes, which improves software interoperability, as inspired by LAMMPS classes and user packages. This article describes version 0.25.1 with benchmarks, compilation instructions, and introductory tutorials for running, restarting, and testing simulations, user guidelines, software design strategies, alternative interfaces, and the test-driven development strategy., (© 2024 Author(s).)

We develop a multiscale coarse-grain model of the NIST Monoclonal Antibody Reference Material 8671 (NISTmAb) to enable systematic computational investigations of high-concentration physical instabilities such as phase separation, clustering, and aggregation. Our multiscale coarse-graining strategy captures atomic-resolution interactions with a computational approach that is orders of magnitude more efficient than atomistic models, assuming the biomolecule can be decomposed into one or more rigid bodies with known, fixed structures. This method reduces interactions between tens of thousands of atoms to a single anisotropic interaction site. The anisotropic interaction between unique pairs of rigid bodies is precomputed over a discrete set of relative orientations and stored, allowing interactions between arbitrarily oriented rigid bodies to be interpolated from the precomputed table during coarse-grained Monte Carlo simulations. We present this approach for lysozyme and lactoferrin as a single rigid body and for the NISTmAb as three rigid bodies bound by a flexible hinge with an implicit solvent model. This coarse-graining strategy predicts experimentally measured radius of gyration and second osmotic virial coefficient data, enabling routine Monte Carlo simulation of medically relevant concentrations of interacting proteins while retaining atomistic detail. All methodologies used in this work are available in the open-source software Free Energy and Advanced Sampling Simulation Toolkit., (© 2024 Author(s).)

The solution viscosity and protein-protein interactions (PPIs) as a function of temperature (4-40 °C) were measured at a series of protein concentrations for a monoclonal antibody (mAb) with different formulation conditions, which include NaCl and sucrose. The flow activation energy ( E η ) was extracted from the temperature dependence of solution viscosity using the Arrhenius equation. PPIs were quantified via the protein diffusion interaction parameter ( k D ) measured by dynamic light scattering, together with the osmotic second virial coefficient and the structure factor obtained through small-angle X-ray scattering. Both viscosity and PPIs were found to vary with the formulation conditions. Adding NaCl introduces an attractive interaction but leads to a significant reduction in the viscosity. However, adding sucrose enhances an overall repulsive effect and leads to a slight decrease in viscosity. Thus, the averaged (attractive or repulsive) PPI information is not a good indicator of viscosity at high protein concentrations for the mAb studied here. Instead, a correlation based on the temperature dependence of viscosity (i.e., E η ) and the temperature sensitivity in PPIs was observed for this specific mAb. When k D is more sensitive to the temperature variation, it corresponds to a larger value of E η and thus a higher viscosity in concentrated protein solutions. When k D is less sensitive to temperature change, it corresponds to a smaller value of E η and thus a lower viscosity at high protein concentrations. Rather than the absolute value of PPIs at a given temperature, our results show that the temperature sensitivity of PPIs may be a more useful metric for predicting issues with high viscosity of concentrated solutions. In addition, we also demonstrate that caution is required in choosing a proper protein concentration range to extract k D . In some excipient conditions studied here, the appropriate protein concentration range needs to be less than 4 mg/mL, remarkably lower than the typical concentration range used in the literature.

Theories of small systems play an important role in the fundamental understanding of finite size effects in statistical mechanics, as well as the validation of molecular simulation results as no computer can simulate fluids in the thermodynamic limit. Previously, a shell particle was included in the isothermal-isobaric ensemble in order to resolve an ambiguity in the resulting partition function. The shell particle removed either redundant volume states or redundant translational degrees of freedom of the system and yielded quantitative differences from traditional simulations in this ensemble. In this work, we investigate the effect of including a shell particle in the canonical, grand canonical, and Gibbs ensembles. For systems comprised of a pure component ideal gas, analytical expressions for various thermodynamic properties are obtained. We also derive the Metropolis Monte Carlo simulation acceptance criteria for these ensembles with shell particles, and the results of the simulations of an ideal gas are in excellent agreement with the theoretical predictions. The system size dependence of various important ensemble averages is also analyzed., (© 2024 Author(s). Published under an exclusive license by AIP Publishing.)

Molecular simulations of water adsorption in porous materials often converge slowly due to sampling bottlenecks that follow from hydrogen bonding and, in many cases, the formation of water clusters. These effects may be exacerbated in metal-organic framework (MOF) adsorbents, due to the presence of pore spaces (cages) that promote the formation of discrete-size clusters and hydrophobic effects (if present), among other reasons. In Grand Canonical Monte Carlo (MC) simulations, these sampling challenges are typically manifested by low MC acceptance ratios, a tendency for the simulation to become stuck in a particular loading state (i.e., macrostates ), and the persistence of specific clusters for long periods of the simulation. We present simulation strategies to address these sampling challenges, by applying flat-histogram MC (FHMC) methods and specialized MC move types to simulations of water adsorption. FHMC, in both Transition-matrix and Wang-Landau forms, drives the simulation to sample relevant macrostates by incorporating weights that are self-consistently adjusted throughout the simulation and generate the macrostate probability distribution (MPD). Specialized MC moves, based on aggregation-volume bias and configurational bias methods, separately address low acceptance ratios for basic MC trial moves and specifically target water molecules in clusters; in turn, the specialized MC moves improve the efficiency of generating new configurations which is ultimately reflected in improved statistics collected by FHMC. The combined strategies are applied to study the adsorption of water in CuBTC and ZIF-8 at 300 K, through examination of the MPD and the adsorption isotherm generated by histogram reweighting. A key result is the appearance of nontrivial oscillations in the MPD, which we show to be associated with water clusters in the adsorption system. Additionally, we show that the probabilities of certain clusters become similar in value near the boundaries of the isotherm hysteresis loop, indicating a strong connection between cluster formation/destruction and the thermodynamic limits of stability. For a hydrophobic MOF, the FHMC results show that the phase transition from low density to high density is suppressed to water pressure far above the bulk-fluid saturation pressure; this is consistent with results presented elsewhere. We also compare our FHMC simulation isotherm to one measured by a different technique but with ostensibly the same molecular interactions and comment on observed differences and the need for follow-up work. The simulation strategies presented here can be applied to the simulation of water in other MOFs using heuristic guidelines laid out in our text, which should facilitate the more consistent and efficient simulation of water adsorption in porous materials in future applications.

We simultaneously designed the porosity and plane symmetry of self-assembling colloidal films by using isohedral tiles to determine the location and shape of enthalpically interacting surface patches on motifs being functionalized. The symmetries of both the tile and motif determine the plane symmetry group of the final assembly. Previous work has either ignored symmetry considerations altogether or accounted for only the tile's properties, applicable only when the motif is asymmetric; this approach provides a complete account and enables the design of symmetric colloids using this tile-based approach, which are often more practical to manufacture. We present the methodology, based on the type of the tile, and provide computational tools that enable the automatic classification of all tiles for a given motif and the optimization of the tile to fit the motif, sometimes referred to as "Escherization".

Monoclonal antibodies (mAbs) make up a major class of biotherapeutics with a wide range of clinical applications. Their physical stability can be affected by various environmental factors. For instance, an acidic pH can be encountered during different stages of the mAb manufacturing process, including purification and storage. Therefore, understanding the behavior of flexible mAb molecules in acidic solution environments will benefit the development of stable mAb products. This study used small-angle X-ray scattering (SAXS) and complementary biophysical characterization techniques to investigate the conformational flexibility and protein-protein interactions (PPI) of a model mAb molecule under near-neutral and acidic conditions. The study also characterized the interactions between Fab and Fc fragments under the same buffer conditions to identify domain-domain interactions. The results suggest that solution pH significantly influences mAb flexibility and thus could help mAbs remain physically stable by maximizing local electrostatic repulsions when mAbs become crowded in solution. Under acidic buffer conditions, both Fab and Fc contribute to the repulsive PPI observed among the full mAb at a low ionic strength. However, as ionic strength increases, hydrophobic interactions lead to the self-association of Fc fragments and, subsequently, could affect the aggregation state of the mAb.

We introduce Gaussian Process Regression (GPR) as an enhanced method of thermodynamic extrapolation and interpolation. The heteroscedastic GPR models that we introduce automatically weight provided information by its estimated uncertainty, allowing for the incorporation of highly uncertain, high-order derivative information. By the linearity of the derivative operator, GPR models naturally handle derivative information and, with appropriate likelihood models that incorporate heterogeneous uncertainties, are able to identify estimates of functions for which the provided observations and derivatives are inconsistent due to the sampling bias that is common in molecular simulations. Since we utilize kernels that form complete bases on the function space to be learned, the estimated uncertainty in the model takes into account that of the functional form itself, in contrast to polynomial interpolation, which explicitly assumes the functional form to be fixed. We apply GPR models to a variety of data sources and assess various active learning strategies, identifying when specific options will be most useful. Our active-learning data collection based on GPR models incorporating derivative information is finally applied to tracing vapor-liquid equilibrium for a single-component Lennard-Jones fluid, which we show represents a powerful generalization to previous extrapolation strategies and Gibbs-Duhem integration. A suite of tools implementing these methods is provided at https://github.com/usnistgov/thermo-extrap., (© 2023 Public Domain.)

Recent interest in parallelizing flat-histogram transition-matrix Monte Carlo simulations in the grand canonical ensemble, due to its demonstrated effectiveness in studying phase behavior, self-assembly and adsorption, has led to the most extreme case of single-macrostate simulations, where each macrostate is simulated independently with ghost particle insertions and deletions. Despite their use in several studies, no efficiency comparisons of these single-macrostate simulations have been made with multiple-macrostate simulations. We show that multiple-macrostate simulations are up to 3 orders of magnitude more efficient than single-macrostate simulations, which demonstrates the remarkable efficiency of flat-histogram biased insertions and deletions, even with low acceptance probabilities. Efficiency comparisons were made for supercritical fluids and vapor-liquid equilibrium of bulk Lennard-Jones and a three-site water model, self-assembling patchy trimer particles and adsorption of a Lennard-Jones fluid confined in a purely repulsive porous network, using the open source simulation toolkit FEASST. By directly comparing with a variety of Monte Carlo trial move sets, this efficiency loss in single-macrostate simulations is attributed to three related reasons. First, ghost particle insertions and deletions in single-macrostate simulations incur the same computational expense as grand canonical ensemble trials in multiple-macrostate simulations, yet ghost trials do not reap the sampling benefit from propagating the Markov chain to a new microstate. Second, single-macrostate simulations lack macrostate change trials that are biased by the self-consistently converging relative macrostate probability, which is a major component of flat histogram simulations. Third, limiting a Markov chain to a single macrostate reduces sampling possibilities. Existing parallelization methods for multiple-macrostate flat-histogram simulations are shown to be more efficient than parallel single-macrostate simulations by approximately an order of magnitude or more in all systems investigated.

Marine environmental monitoring efforts often rely on the bioaccumulation of persistent anthropogenic contaminants in organisms to create a spatiotemporal record of the ecosystem. Intercorrelation results from the origin, uptake, and transport of these contaminants throughout the ecosystem and may be affected by organism-specific processes such as biotransformation. Here, we explore trends that machine learning tools reveal about a large, recently released environmental chemistry data set of common anthropogenic pollutants measured in the eggs of five seabird species from the North Pacific Ocean. We modeled these data with a variety of machine learning approaches and found models that could accurately determine a range of taxonomic and spatiotemporal trends. We illustrate a general workflow and set of analysis tools that can be used to identify interpretable models which perform nearly as well as state-of-the-art "black boxes." For example, we found shallow decision trees that could resolve genus with greater than 96% accuracy using as few as two analytes and a k -nearest neighbor classifier that could resolve species differences with more than 94% accuracy using only five analytes. The benefits of interpretability outweighed the marginally improved accuracy of more complex models. This demonstrates how machine learning may be used to discover rational, quantitative trends in these systems.

Computational screening of adsorbent materials often uses the Henry's law constant ( K H ) (at a particular temperature) as a first discriminator metric due to its relative ease of calculation. The isosteric heat of adsorption in the limit of zero pressure ( q st ∞ ) is often calculated along with the Henry's law constant, and both properties are informative metrics of adsorbent material performance at low-pressure conditions. In this article, we introduce a method for extrapolating K H as a function of temperature, using series-expansion coefficients that are easily computed at the same time as K H itself; the extrapolation function also yields q st ∞ . The extrapolation is highly accurate over a wide range of temperatures when the basis temperature is sufficiently high, for a wide range of adsorbent materials and adsorbate gases. Various results suggest that the extrapolation is accurate when the extrapolation range in inverse-temperature space is limited to |β - β 0 | < 0.5 mol/kJ. Application of the extrapolation to a large set of materials is shown to be successful provided that K H is not extremely large and/or the extrapolation coefficients converge satisfactorily. The extrapolation is also able to predict q st ∞ for a system that shows an unusually large temperature dependence. The work provides a robust method for predicting K H and q st ∞ over a wide range of industrially relevant temperatures with minimal effort beyond that necessary to compute those properties at a single temperature, which facilitates the addition of practical operating (or processing) conditions to computational screening exercises.

We describe a method for deriving surface functionalization patterns for colloidal systems that can induce self-assembly into any chosen periodic symmetry at a planar interface. The result is a sequence of letters, s ∈ {A,T,C,G}, or a gene, that describes the perimeter of the colloidal object and programs its self-assembly. This represents a genome that is finite and can be exhaustively enumerated. These genes derive from symmetry, which may be topologically represented by two-dimensional parabolic orbifolds; since these orbifolds are surfaces that may be derived from first principles, this represents an ab initio route to colloid functionality. The genes are human readable and can be employed to easily design colloidal units. We employ a biological (genetic) analogy to demonstrate this and illustrate their connection to the designs of Maurits Cornelis (M. C.) Escher.

Variational autoencoders (VAEs) are rapidly gaining popularity within molecular simulation for discovering low-dimensional, or latent, representations, which are critical for both analyzing and accelerating simulations. However, it remains unclear how the information a VAE learns is connected to its probabilistic structure and, in turn, its loss function. Previous studies have focused on feature engineering, ad hoc modifications to loss functions, or adjustment of the prior to enforce desirable latent space properties. By applying effectively arbitrarily flexible priors via normalizing flows, we focus instead on how adjusting the structure of the decoding model impacts the learned latent coordinate. We systematically adjust the power and flexibility of the decoding distribution, observing that this has a significant impact on the structure of the latent space as measured by a suite of metrics developed in this work. By also varying weights on separate terms within each VAE loss function, we show that the level of detail encoded can be further tuned. This provides practical guidance for utilizing VAEs to extract varying resolutions of low-dimensional information from molecular dynamics and Monte Carlo simulations.

Discovering meaningful collective variables for enhancing sampling, via applied biasing potentials or tailored MC move sets, remains a major challenge within molecular simulation. While recent studies identifying collective variables with variational autoencoders (VAEs) have focused on the encoding and latent space discovered by a VAE, the impact of the decoding and its ability to act as a generative model remains unexplored. We demonstrate how VAEs may be used to learn (on-the-fly and with minimal human intervention) highly efficient, collective Monte Carlo moves that accelerate sampling along the learned collective variable. In contrast to many machine learning-based efforts to bias sampling and generate novel configurations, our methods result in exact sampling in the ensemble of interest and do not require reweighting. In fact, we show that the acceptance rates of our moves approach unity for a perfect VAE model. While this is never observed in practice, VAE-based Monte Carlo moves still enhance sampling of new configurations. We demonstrate, however, that the form of the encoding and decoding distributions, in particular the extent to which the decoder reflects the underlying physics, greatly impacts the performance of the trained VAE.

We derive properties of self-assembling rings which can template the organization of an arbitrary colloid into any periodic symmetry in two Euclidean dimensions. By viewing this as a tiling problem, we illustrate how the shape and chemical patterning of these rings are derivable, and are explicitly reflected by the symmetry group's orbifold symbol. We performed molecular dynamics simulations to observe their self-assembly and found 5 different characteristics which could be easily rationalized on the basis of this symbol. These include systems which undergo chiral phase separation, are addressably complex, exhibit self-limiting growth into clusters, form ordered "rods" in only one-dimension akin to a smectic phase, and those from symmetry groups which are pluripotent and allow one to select rings which exhibit different behaviors. We discuss how the curvature of the ring's edges plays an integral role in achieving correct self-assembly, and illustrate how to obtain these shapes. This provides a method for patterning colloidal systems at interfaces without explicitly programming this information onto the colloid itself.

The most direct approach to determining if two aqueous solutions will phase-separate upon mixing is to exhaustively screen them in a pair-wise fashion. This is a time-consuming process that involves preparation of numerous stock solutions, precise transfer of highly concentrated and often viscous solutions, exhaustive agitation to ensure thorough mixing, and time-sensitive monitoring to observe the presence of emulsion characteristics indicative of phase separation. Here, we examined the pair-wise mixing behavior of 68 water-soluble compounds by observing the formation of microscopic phase boundaries and droplets of 2278 unique 2-component solutions. A series of machine learning classifiers (artificial neural network, random forest, k-nearest neighbors, and support vector classifier) were then trained on physicochemical property data associated with the 68 compounds and used to predict their miscibility upon mixing. Miscibility predictions were then compared to the experimental observations. The random forest classifier was the most successful classifier of those tested, displaying an average receiver operator characteristic area under the curve of 0.74. The random forest classifier was validated by removing either one or two compounds from the input data, training the classifier on the remaining data and then predicting the miscibility of solutions involving the removed compound(s) using the classifier. The accuracy, specificity, and sensitivity of the random forest classifier were 0.74, 0.80, and 0.51, respectively, when one of the two compounds to be examined was not represented in the training data. When asked to predict the miscibility of two compounds, neither of which were represented in the training data, the accuracy, specificity, and sensitivity values for the random forest classifier were 0.70, 0.82 and 0.29, respectively. Thus, there is potential for this machine learning approach to improve the design of screening experiments to accelerate the discovery of aqueous two-phase systems for numerous scientific and industrial applications.

We investigate the thermodynamic properties of various super-critical model adsorptive systems with different fluid-solid attractive strengths using the confined-density virial expansion, with coefficients calculated using the Mayer-sampling Monte Carlo method up to fifth order. We find that the virial expansion converges for adsorptive systems over a density range corresponding approximately to the film-formation regime. Beyond this regime, higher order effects become increasingly important. The virial expansion of the density profile is also investigated. It is determined that this expansion gives insight into the structure associated with adsorption. We also find that weakly attractive systems have a more negative second virial coefficient than strongly attractive systems. This runs counter to the usual interpretation of bulk fluid virial coefficients. This is due to the infinite-dilution limit being very different for adsorbed fluids compared to bulk fluids.

Thermodynamic extrapolation has previously been used to predict arbitrary structural observables in molecular simulations at temperatures (or relative chemical potentials in open-system mixtures) different from those at which the simulation was performed. This greatly reduces the computational cost in mapping out phase and structural transitions. In this work, we explore the limitations and accuracy of thermodynamic extrapolation applied to water, where qualitative shifts from anomalous to simple-fluid-like behavior are manifested through shifts in the liquid structure that occur as a function of both temperature and density. We present formulas for extrapolating in volume for canonical ensembles and demonstrate that linear extrapolations of water's structural properties are only accurate over a limited density range. On the other hand, linear extrapolation in temperature can be accurate across the entire liquid state. We contrast these extrapolations with classical perturbation theory techniques, which are more conservative and slowly converging. Indeed, we show that such behavior is expected by demonstrating exact relationships between extrapolation of free energies and well-known techniques to predict free energy differences. An ideal gas in an external field is also studied to more clearly explain these results for a toy system with fully analytical solutions. We also present a recursive interpolation strategy for predicting arbitrary structural properties of molecular fluids over a predefined range of state conditions, demonstrating its success in mapping qualitative shifts in water structure with density.

The accurate prediction of stable crystalline phases is a long-standing problem encountered in the study of conventional atomic and molecular solids as well as soft materials. One possible solution involves enumerating a reasonable set of candidate structures and then screening them to identify the one(s) with the lowest (free) energy. Candidate structures in this set can also serve as starting points for other routines, such as genetic algorithms, which search via optimization. Here, we present a framework for crystal structure enumeration of two-dimensional systems that utilizes a combination of symmetry- and stoichiometry-imposed constraints to compute valid configurations of particles that tile Euclidean space. With mild assumptions, this produces a computationally tractable total number of proposed candidates, enabling multicomponent systems to be screened by direct enumeration of possible crystalline ground states. The python code that enables these calculations is available at https://github.com/usnistgov/PACCS.

Inverse design methods are powerful computational approaches for creating colloidal systems which self-assemble into a target morphology by reverse engineering the Hamiltonian of the system. Despite this, these optimization procedures tend to yield Hamiltonians which are too complex to be experimentally realized. An alternative route to complex structures involves the use of several different components, however, conventional inverse design methods do not explicitly account for the possibility of phase separation into compositionally distinct structures. Here, we present an inverse design scheme for multicomponent colloidal systems by combining active learning with a method to directly compute their ground state phase diagrams. This explicitly accounts for phase separation and can locate stable regions of Hamiltonian parameter space which grid-based surveys are prone to miss. Using this we design low-density, binary structures with Lennard-Jones-like pairwise interactions that are simpler than in the single component case and potentially realizable in an experimental setting. This reinforces the concept that ground states of simple, multicomponent systems might be rich with previously unappreciated diversity, enabling the assembly of non-trivial structures with only few simple components instead of a single complex one.

The phenomenon of dynamic arrest, more commonly referred to as gel and glass formation, originates as particle motion slows significantly. Current understanding of gels and glasses stems primarily from dispersions of spherical particles, but much less is known about how particle shape affects dynamic arrest transitions. To better understand the effects of particle shape anisotropy on gel and glass formation, we systematically measure the rheology, particle dynamics, and static microstructure of thermoreversible colloidal dispersions of adhesive hard rods (AHR). First, the dynamic arrest transitions are mapped as a function of temperature T, aspect ratio L/D≈ 3 to 7, and volume fraction φ≈ 0.1 to 0.5. The critical gel temperature T gel and glass volume fraction φ g are determined from the particle dynamics and rheology. Second, an effective orientation-averaged, short-range attraction between rods is quantified from small-angle scattering measurements and characterized by a reduced temperature τ. Similar τ is found at low rod concentrations, indicating that rod gelation occurs at similar effective attraction strength independent of L/D. Monte Carlo simulations reveal a similar convergence in τ when rods cluster and percolate with an average bond coordination number 〈n c 〉≈ 2.4, supporting the link between physical gelation and rigidity percolation. Lastly, AHR results are mapped onto a dimensionless state diagram to compare with previous predictions of attraction-driven gels, repulsion-driven glasses, and liquid crystal phases.

Protein therapeutics are vitally important clinically and commercially, with monoclonal antibody (mAb) therapeutic sales alone accounting for $115 billion in revenue for 2018.[1] In order for these therapeutics to be safe and efficacious, their protein components must maintain their high order structure (HOS), which includes retaining their three-dimensional fold and not forming aggregates. As demonstrated in the recent NISTmAb Interlaboratory nuclear magnetic resonance (NMR) Study[2], NMR spectroscopy is a robust and precise approach to address this HOS measurement need. Using the NISTmAb study data, we benchmark a procedure for automated outlier detection used to identify spectra that are not of sufficient quality for further automated analysis. When applied to a diverse collection of all 252 1 H, 13 C gHSQC spectra from the study, a recursive version of the automated procedure performed comparably to visual analysis, and identified three outlier cases that were missed by the human analyst. In total, this method represents a distinct advance in chemometric detection of outliers due to variation in both measurement and sample., Competing Interests: Declaration of Interest The authors have no competing interests to declare.

Vast numbers of unstudied hypothetical porous frameworks continue to spark interest in optimizing adsorption and catalytic processes. Evaluating the use of such materials depends on the accessibility of thermodynamic metrics such as the free energy, which, in turn, depend on the satisfactory estimation or calculation of the adsorption entropy, which often remains elusive. Previous works using simulations and experimental data have demonstrated relationships between the entropy and system descriptors, allowing for sensible predictions based on more-easily obtained physical parameters. However, the resultant conclusions were either based on experimental data for industrially relevant alkanes or lacked a significant sample size. In this paper, we evaluate correlations between gas-phase and adsorbed-phase entropies for a larger and more chemically diverse set of adsorbate molecules by using force fields and statistical mechanical expressions to calculate those entropies. In total, we perform calculations for 37 molecules across 10 chemical categories available in the TraPPE force field set, as adsorbed in five siliceous zeolites. Our results show that linear correlations between the gas- and adsorbed-phase entropies persist for the larger and diverse set of adsorbate molecules studied here, proving a broader applicability and justifying the use of simple correlations for many adsorbates and, presumably, adsorbent materials.

While ionic liquids have promising applications as industrial solvents, predicting their fluid phase properties and coexistence remains a challenge. Grand canonical Monte Carlo simulation is an effective method for such predictions, but equilibration is hampered by the apparent requirement to insert and delete neutral sets of ions simultaneously in order to maintain charge neutrality. For relatively high densities and low temperatures, previously developed methods have been shown to be essential in improving equilibration by gradual insertion and deletion of these neutral sets of ions. We introduce an expanded ensemble approach which may be used in conjunction with these existing methods to further improve efficiency. Individual ions are inserted or deleted in one Monte Carlo trial rather than simultaneous insertion/deletion of neutral sets. We show how charge neutrality is maintained and show rigorous quantitative agreement between the conventional and the proposed expanded ensemble approaches, but with up to an order of magnitude increase in efficiency at high densities. The expanded ensemble approach is also more straightforward to implement than simultaneous insertion/deletion of neutral sets, and its implementation is demonstrated within open source software.

We demonstrate a method based on symmetry to predict the structure of self-assembling, multicomponent colloidal mixtures. This method allows us to feasibly enumerate candidate structures from all symmetry groups and is many orders of magnitude more computationally efficient than combinatorial enumeration of these candidates. In turn, this permits us to compute ground-state phase diagrams for multicomponent systems. While tuning the interparticle potentials to produce potentially complex interactions represents the conventional route to designing exotic lattices, we use this scheme to demonstrate that simple potentials can also give rise to such structures which are thermodynamically stable at moderate to low temperatures. Furthermore, for a model two-dimensional colloidal system, we illustrate that lattices forming a complete set of 2-, 3-, 4-, and 6-fold rotational symmetries can be rationally designed from certain systems by tuning the mixture composition alone, demonstrating that stoichiometric control can be a tool as powerful as directly tuning the interparticle potentials themselves.

Employing 2 different coarse-grained models, we evaluated the effect of intramolecular domain-domain distances and hinge flexibility on the general solution structure of monoclonal antibodies (mAbs), within the context of protein-protein steric repulsion. These models explicitly account for the hinge region, and represent antibodies at either domain or subdomain levels (i.e., 4-bead and 7-bead representations, respectively). Additionally, different levels of mAb flexibility are also considered. When evaluating mAbs as rigid structures, analysis of small-angle scattering profiles showed that changes in the relative internal distances between Fc and Fab domains significantly alter the local arrangement of neighboring molecules, as well as the molecular packing of the concentrated mAb solutions. Likewise, enabling hinge flexibility in either of the mAb models led to qualitatively similar results, where flexibility increases the spatial molecular arrangement at elevated concentrations. This occurs because fluctuations in mAb quaternary structure are modulated by the close proximity between molecules at elevated concentrations (>50 mg mL -1 ), yielding an increased molecular packing and osmotic compressibility. However, our results also showed that the mechanism behind this synergy between flexibility and packing strongly depends on both the level of structural detail and the number of degrees-of-freedom considered in the coarse-grained model., (Copyright © 2019 American Pharmacists Association®. All rights reserved.)

We present analytic expressions for the second virial coefficient of hard particles in slit and cylindrical pores with arbitrary square-well fluid-fluid and fluid-solid interactions, number of fluid components, and pore sizes. We derive multiple methods to predict the isotherms in confined materials based on the virial expansion and find that the thermodynamic and structural properties calculated from the virial expansion to second order in the dilute limit display excellent qualitative agreement with previous simulation results.

The creation of molecular or colloidal building blocks which can self-assemble into complex, ordered porous structures has been long sought-after, and so are the guiding principles behind this creation. The pursuit of this goal has led to the creation of novel classes of materials like metal organic frameworks (MOFs) and covalent organic frameworks (COFs). In theory, a tremendous number of structures can be formed by these materials due to the variety of geometries available to their building blocks. However, most realized crystal structures tend to be simple or homoporous and typically assemble from building blocks with high degrees of symmetry. Building blocks with low degrees of symmetry suitable for assembly into the more complex structures tend to assemble into polymorphous or disordered structures instead. In this work, we use Monte Carlo simulations of patchy vertex-like building blocks to show how the addition of chemical specificity via orthogonally reacting functional sites can allow vertex-like building blocks with even asymmetric geometries to self-assemble into ordered crystallites of various complex structures. In addition to demonstrating the utility of such a strategy in creating ordered, heteroporous structures, we also demonstrate that it can be used as a means for tuning specific features of the crystal structure, accomplishing such aims as the control of relative pore sizes. We also discuss heuristics for properly designing molecules so that they can assemble into target structures.

Cylindrical or rod-like particles are promising materials for the applications of fillers in nanocomposite materials and additives to control rheological properties of colloidal suspensions. Recent advances in particle synthesis allows for cylinders to be manufactured with short-ranged attractions to study the gelation as a function of packing fraction, aspect ratio and attraction strength. In order to aid in the analysis of small-angle scattering experiments of rod-like particles, computer simulation methods were used to model these particles with specialized Monte Carlo algorithms and tabular superquadric potentials. The attractive interaction between neighboring rods increases with the amount of locally-accessible surface area, thus leading to patchy-like interactions. We characterize the clustering and percolation of cylinders as the attractive interaction increases from the homogenous fluid at relatively low attraction strength, for a variety of aspect ratios and packing fractions. Comparisons with the experimental scattering results are also presented, which are in agreement.

The theoretical framework to evaluate small-angle scattering (SAS) profiles for multi-component macromolecular solutions is re-examined from the standpoint of molecular simulations in the grand-canonical ensemble, where the chemical potentials of all species in solution are fixed. This statistical mechanical ensemble resembles more closely scattering experiments, capturing concentration fluctuations that arise from the exchange of molecules between the scattering volume and the bulk solution. The resulting grand-canonical expression relates scattering intensities to the different intra- and intermolecular pair distribution functions, as well as to the distribution of molecular concentrations on the scattering volume. This formulation represents a generalized expression that encompasses most of the existing methods to evaluate SAS profiles from molecular simulations. The grand-canonical SAS methodology is probed for a series of different implicit-solvent, homogeneous systems at conditions ranging from dilute to concentrated. These systems consist of spherical colloids, dumbbell particles, and highly flexible polymer chains. Comparison of the resulting SAS curves against classical methodologies based on either theoretical approaches or canonical simulations (i.e., at a fixed number of molecules) shows equivalence between the different scattering intensities so long as interactions between molecules are net repulsive or weakly attractive. On the other hand, for strongly attractive interactions, grand-canonical SAS profiles deviate in the low- and intermediate- q range from those calculated in a canonical ensemble. Such differences are due to the distribution of molecules becoming asymmetric, which yields a higher contribution from configurations with molecular concentrations larger than the nominal value. Additionally, for flexible systems, explicit discrimination between intra- and inter-molecular SAS contributions permits the implementation of model-free, structural analysis such as Guinier's plots at high molecular concentrations, beyond what the traditional limits are for such analysis.

Binary superlattices constructed from nano- or micron-sized colloidal particles have a wide variety of applications, including the design of advanced materials. Self-assembly of such crystals from their constituent colloids can be achieved in practice by, among other means, the functionalization of colloid surfaces with single-stranded DNA sequences. However, when driven by DNA, this assembly is traditionally premised on the pairwise interaction between a single DNA sequence and its complement, and often relies on particle size asymmetry to entropically control the crystalline arrangement of its constituents. The recently proposed "multi-flavoring" motif for DNA functionalization, wherein multiple distinct strands of DNA are grafted in different ratios to different colloids, can be used to experimentally realize a binary mixture in which all pairwise interactions are independently controllable. In this work, we use various computational methods, including molecular dynamics and Wang-Landau Monte Carlo simulations, to study a multi-flavored binary system of micron-sized DNA-functionalized particles modeled implicitly by Fermi-Jagla pairwise interactions. We show how self-assembly of such systems can be controlled in a purely enthalpic manner, and by tuning only the interactions between like particles, demonstrate assembly into various morphologies. Although polymorphism is present over a wide range of pairwise interaction strengths, we show that careful selection of interactions can lead to the generation of pure compositionally ordered crystals. Additionally, we show how the crystal composition changes with the like-pair interaction strengths, and how the solution stoichiometry affects the assembled structures.

We describe a methodology for extrapolating the structural properties of multicomponent fluids from one thermodynamic state to another. These properties generally include features of a system that may be computed from an individual configuration such as radial distribution functions, cluster size distributions, or a polymer's radius of gyration. This approach is based on the principle of using fluctuations in a system's extensive thermodynamic variables, such as energy, to construct an appropriate Taylor series expansion for these structural properties in terms of intensive conjugate variables, such as temperature. Thus, one may extrapolate these properties from one state to another when the series is truncated to some finite order. We demonstrate this extrapolation for simple and coarse-grained fluids in both the canonical and grand canonical ensembles, in terms of both temperatures and the chemical potentials of different components. The results show that this method is able to reasonably approximate structural properties of such fluids over a broad range of conditions. Consequently, this methodology may be employed to increase the computational efficiency of molecular simulations used to measure the structural properties of certain fluid systems, especially those used in high-throughput or data-driven investigations.

We used flat-histogram sampling Monte Carlo to study capillary phase transitions in deformable adsorbent materials. Specifically, we considered a pure adsorbate fluid below its bulk critical temperature within a slit pore of variable pore width. The instantaneous pore width is dictated by a number of factors, such as adsorbate loading, reservoir pressure, fluid-wall interaction, and bare adsorbent properties. In the slit pores studied here, the bare adsorbent free energy was assumed to be biparabolic, consisting of two preferential pore configurations, namely, the narrow pore and the large pore configurations. Four distinct phases could be found in the adsorption isotherms. We found a low-pressure phase transition, driven primarily by capillary condensation/evaporation and accompanied by adsorbent deformation in response. The deformation can be a relatively small contraction/expansion as seen in elastic materials, or a large-scale structural transformation of the adsorbent. We also found a high-pressure transition driven by excluded volume effects, which tends to expand the material and thus results in a large-scale structural transformation of the adsorbent. The adsorption isotherms and osmotic free energies can be rationalized by considering the relative free energy differences between the basins of the bare adsorbent free energy.

The Ninth Industrial Fluid Properties Simulation Challenge aimed to test the ability of molecular modeling approaches to predict water/oil interfacial tension (IFT) at conditions of high temperature and pressure. In particular, the challenge featured water/oil IFT where the oil was n-dodecane, toluene, or a 50:50 n-dodecane/toluene blend at 1.825 MPa and temperatures in the range of 383 K to 443 K. Seven entries were received including approaches such as molecular dynamics (MD) and Monte Carlo (MC) simulations, COSMO-RS, and iSAFT, and they were judged by comparison to pendant drop tensiometer benchmark data. The quality of predictions varied among the entries between approximately 20 % and 70 % of the total points possible with the entries based on MD and MC having the highest scores in most cases. As is often the case in molecular modeling, predictions of the relative trends tended to be reliable even if the absolute values were not.

We derive an approach for extrapolating the free energy landscape of multicomponent systems in the grand canonical ensemble, obtained from flat-histogram Monte Carlo simulations, from one set of temperature and chemical potentials to another. This is accomplished by expanding the landscape in a Taylor series at each value of the order parameter which defines its macrostate phase space. The coefficients in each Taylor polynomial are known exactly from fluctuation formulas, which may be computed by measuring the appropriate moments of extensive variables that fluctuate in this ensemble. Here we derive the expressions necessary to define these coefficients up to arbitrary order. In principle, this enables a single flat-histogram simulation to provide complete thermodynamic information over a broad range of temperatures and chemical potentials. Using this, we also show how to combine a small number of simulations, each performed at different conditions, in a thermodynamically consistent fashion to accurately compute properties at arbitrary temperatures and chemical potentials. This method may significantly increase the computational efficiency of biased grand canonical Monte Carlo simulations, especially for multicomponent mixtures. Although approximate, this approach is amenable to high-throughput and data-intensive investigations where it is preferable to have a large quantity of reasonably accurate simulation data, rather than a smaller amount with a higher accuracy.

Virial coefficients are predicted over a large range of both temperatures and model parameter values (i.e., alchemical transformation) from an individual Mayer-sampling Monte Carlo simulation by statistical mechanical extrapolation with minimal increase in computational cost. With this extrapolation method, a Mayer-sampling Monte Carlo simulation of the SPC/E (extended simple point charge) water model quantitatively predicted the second virial coefficient as a continuous function spanning over four orders of magnitude in value and over three orders of magnitude in temperature with less than a 2% deviation. In addition, the same simulation predicted the second virial coefficient if the site charges were scaled by a constant factor, from an increase of 40% down to zero charge. This method is also shown to perform well for the third virial coefficient and the exponential parameter for a Lennard-Jones fluid.

Using Monte Carlo and molecular dynamics simulations, we examine the adsorption of methane in cylindrical silica mesopores in an effort to understand a possible phase transition of adsorbed methane in MCM-41 and SBA-15 silica that was previously identified by an unexpected increase in the adsorbed fluid density following capillary condensation, as measured by small-angle neutron scattering (SANS) [Chiang, W-S., et al., Langmuir 2016, 32, 8849]. Our initial simulation results identify a roughly 10 % increase in the density of the liquidlike adsorbed phase for either an isotherm with increasing pressure or an isobar with decreasing temperature and that this densification is associated with a local maximum in the isosteric enthalpy of adsorption. Subsequent analysis of the simulated fluid, via computation of bond-orientational order parameters of specific annular layers of the adsorbed fluid, showed that the layers undergo an ordering transition from a disordered, amorphous state to one with two-dimensional hexagonal structure. Furthermore, this two-dimensional restructuring of the fluid occurs at the same thermodynamic state points as the aforementioned densification and local maximum in the isosteric enthalpy of adsorption. We thus conclude that the densification of the fluid is the result of structural reorganization, which is signaled by the maximum in the isosteric enthalpy. Owing to the qualitative similarity of the structural transitions in the simulated and experimental methane fluids, we propose this hexagonal reorganization as a plausible explanation of the densification observed in SANS measurements. Lastly, we speculate how this structural transition may impact the transport properties of the adsorbed fluid.

Using molecular simulations, we investigate the relationship between the pore-averaged and position-dependent self-diffusivity of a fluid adsorbed in a strongly attractive pore as a function of loading. Previous work (Krekelberg, W. P.; Siderius, D. W.; Shen, V. K.; Truskett, T. M.; Errington, J. R. Connection between thermodynamics and dynamics of simple fluids in highly attractive pores. Langmuir 2013, 29, 14527-14535, doi: 10.1021/la4037327) established that pore-averaged self-diffusivity in the multilayer adsorption regime, where the fluid exhibits a dense film at the pore surface and a lower density interior pore region, is nearly constant as a function of loading. Here we show that this puzzling behavior can be understood in terms of how loading affects the fraction of particles that reside in the film and interior pore regions as well as their distinct dynamics. Specifically, the insensitivity of pore-averaged diffusivity to loading arises from the approximate cancellation of two factors: an increase in the fraction of particles in the higher diffusivity interior pore region with loading and a corresponding decrease in the particle diffusivity in that region. We also find that the position-dependent self-diffusivities scale with the position-dependent density. We present a model for predicting the pore-average self-diffusivity based on the position-dependent self-diffusivity, which captures the unusual characteristics of pore-averaged self-diffusivity in strongly attractive pores over several orders of magnitude.

We systematically investigate the assembly of binary multi-flavored colloidal mixtures in two dimensions. In these mixtures all pairwise interactions between species may be tuned independently. This introduces an additional degree of freedom over more traditional binary mixtures with fixed mixing rules, which is anticipated to open new avenues for directed self-assembly. At present, colloidal self-assembly into non-trivial lattices tends to require either high pressures for isotropically interacting particles, or the introduction of directionally anisotropic interactions. Here we demonstrate tunable assembly into a plethora of structures which requires neither of these conditions. We develop a minimal model that defines a three-dimensional phase space containing one dimension for each pairwise interaction, then employ various computational techniques to map out regions of this phase space in which the system self-assembles into these different morphologies. We then present a mean-field model that is capable of reproducing these results for size-symmetric mixtures, which reveals how to target different structures by tuning pairwise interactions, solution stoichiometry, or both. Concerning particle size asymmetry, we find that domains in this model's phase space, corresponding to different morphologies, tend to undergo a continuous "rotation" whose magnitude is proportional to the size asymmetry. Such continuity enables one to estimate the relative stability of different lattices for arbitrary size asymmetries. Owing to its simplicity and accuracy, we expect this model to serve as a valuable design tool for engineering binary colloidal crystals from multi-flavored components.

We derive a method for extrapolating the grand canonical free energy landscape of a multicomponent fluid system from one temperature to another. Previously, we introduced this statistical mechanical framework for the case where kinetic energy contributions to the classical partition function were neglected for simplicity [N. A. Mahynski et al., J. Chem. Phys. 146, 074101 (2017)]. Here, we generalize the derivation to admit these contributions in order to explicitly illustrate the differences that result. Specifically, we show how factoring out kinetic energy effects a priori, in order to consider only the configurational partition function, leads to simpler mathematical expressions that tend to produce more accurate extrapolations than when these effects are included. We demonstrate this by comparing and contrasting these two approaches for the simple cases of an ideal gas and a non-ideal, square-well fluid.

Using molecular simulations, we investigate how the range of fluid-fluid (adsorbate-adsorbate) interactions and the strength of fluid-solid (adsorbate-adsorbent) interactions impact the strong connection between distinct adsorptive regimes and distinct self-diffusivity regimes reported in [Krekelberg, W. P.; Siderius, D. W.; Shen, V. K.; Truskett, T. M.; Errington, J. R. Langmuir 2013 , 29 , 14527-14535]. Although increasing the fluid-fluid interaction range changes both the thermodynamics and the dynamic properties of adsorbed fluids, the previously reported connection between adsorptive filling regimes and self-diffusivity regimes remains. Increasing the fluid-fluid interaction range leads to enhanced layering and decreased self-diffusivity in the multilayer-formation regime but has little effect on the properties within film-formation and pore-filling regimes. We also find that weakly attractive adsorbents, which do not display distinct multilayer formation, are hard-sphere-like at super- and subcritical temperatures. In this case, the self-diffusivity of the confined and bulk fluid has a nearly identical scaling-relationship with effective density.

We use Monte Carlo simulations to investigate the thermodynamic behavior of soft porous crystal (SPC) adsorbents under the influence of an external barostat. We consider SPCs that naturally exhibit polymorphism between crystal forms of two distinct pore sizes. In the absence of barostatting, these crystals may be naturally divided into two categories depending on their response to stress applied by the adsorbate fluid: those which macroscopically deform and change the volume of their unit cell ("breathing") and those which instead undergo internal rearrangements that change the adsorbate-accessible volume without modifying the unit cell volume ("gate-opening"). When breathing SPCs have a constant external pressure applied, in addition to the thermodynamic pressure of the adsorbate fluid, we find that the free energy difference between the crystal polymorphs is shifted by a constant amount over the entire course of adsorption. Thus, their relative stability may be easily controlled by the barostat. However, when the crystal is held at a fixed overall pressure, changes to the relative stability of the polymorphs tend to be more complex. We demonstrate a thermodynamic analogy between breathing SPCs held at a fixed pressure and macroscopically rigid gate-opening ones which explains this behavior. Furthermore, we illustrate how this implies that external mechanical forces may be employed to tune the effective free energy profile of an empty SPC, which may open new avenues to engineer the thermodynamic properties of these polymorphic adsorbents, such as selectivity.

The self-assembly of patchy trimer particles consisting of one attractive site and two repulsive sites is investigated with nonequilibrium molecular dynamics simulations in the presence of a velocity gradient, as would be produced by the application of a shear stress on the system. As shear is increased, globular-shaped micellar clusters increase in size and become more elongated. The globular clusters are also more stable at higher temperatures in the presence of shear than at equilibrium. These results help to increase our understanding of the effect of shear on self-assembly for a variety of applications.

Measurement of the pore-size distribution (PSD) via gas adsorption and the so-called "kernel method" is a widely used characterization technique for rigid adsorbents. Yet, standard techniques and analytical equipment are not appropriate to characterize the emerging class of flexible adsorbents that deform in response to the stress imparted by an adsorbate gas, as the PSD is a characteristic of the material that varies with the gas pressure and any other external stresses. Here, we derive the PSD for a flexible adsorbent using statistical mechanics in the osmotic ensemble to draw analogy to the kernel method for rigid materials. The resultant PSD is a function of the ensemble constraints including all imposed stresses and, most importantly, the deformation free energy of the adsorbent material. Consequently, a pressure-dependent PSD is a descriptor of the deformation characteristics of an adsorbent and may be the basis of future material characterization techniques. We discuss how, given a technique for resolving pressure-dependent PSDs, the present statistical mechanical theory could enable a new generation of analytical tools that measure and characterize certain intrinsic material properties of flexible adsorbents via otherwise simple adsorption experiments.

We present a method for predicting the free energy landscape of fluids at low temperatures from flat-histogram grand canonical Monte Carlo simulations performed at higher ones. We illustrate our approach for both pure and multicomponent systems using two different sampling methods as a demonstration. This allows us to predict the thermodynamic behavior of systems which undergo both first order and continuous phase transitions upon cooling using simulations performed only at higher temperatures. After surveying a variety of different systems, we identify a range of temperature differences over which the extrapolation of high temperature simulations tends to quantitatively predict the thermodynamic properties of fluids at lower ones. Beyond this range, extrapolation still provides a reasonably well-informed estimate of the free energy landscape; this prediction then requires less computational effort to refine with an additional simulation at the desired temperature than reconstruction of the surface without any initial estimate. In either case, this method significantly increases the computational efficiency of these flat-histogram methods when investigating thermodynamic properties of fluids over a wide range of temperatures. For example, we demonstrate how a binary fluid phase diagram may be quantitatively predicted for many temperatures using only information obtained from a single supercritical state.

We use flat-histogram Monte Carlo simulations to study how changing the flexibility of soft porous crystals (SPCs) affects their selective adsorption of a binary, size-asymmetric supercritical fluid. Specifically, we consider mesoporous SPCs which have multiple minima in their free energy profiles as a function of pore size such that they are capable of exhibiting polymorphism between a narrow and large pore phase. While specific fluid-pore interactions determine the shape of both pores' selectivity curve as a function of adsorbate pressure, an individual pore tends to selectively adsorb a species based on the size of the adsorbate molecule relative to itself, thereby shifting the pore's selectivity curve relative to its polymorph. By controlling the flexibility of a SPC, the relative thermodynamic stability of the two pore phases may be varied, thereby changing the overall selectivity of the SPC during adsorbate loading. We investigate this for two classes of SPCs: one representative of "gate-opening" materials and another of "breathing" materials. For gate-opening materials, this control is much more salient than in breathing ones. However, for the latter, we illustrate how to tune the free energy profile to create materials which breathe multiple times during adsorption/desorption.

We demonstrate an extensible flat-histogram Monte Carlo simulation methodology for studying the adsorption of multicomponent fluids in flexible porous solids. This methodology allows us to easily obtain the complete free energy landscape for the confined fluid-solid system in equilibrium with a bulk fluid of any arbitrary composition. We use this approach to study the adsorption of a prototypical coarse-grained binary fluid in "Hookean" solids, where the free energy of the solid may be described as a simple spring. However, our approach is fully extensible to solids with arbitrarily complex free energy profiles. We demonstrate that by tuning the fluid-solid interaction ranges, the inhomogeneous fluid structure inside the pore can give rise to enhanced selective capture of a larger species through cooperative adsorption with a smaller one. The maximum enhancement in selectivity is observed at low to intermediate pressures and is especially pronounced when the larger species is very dilute in the bulk. This suggest a mechanism by which the selective capture of a minor component from a bulk fluid may be enhanced.

Ultrasonic experiments allow one to measure the elastic modulus of bulk solid or fluid samples. Recently such experiments have been carried out on fluid-saturated nanoporous glass to probe the modulus of a confined fluid. In our previous work [G. Y. Gor et al., J. Chem. Phys., 143, 194506 (2015)], using Monte Carlo simulations we showed that the elastic modulus K of a fluid confined in a mesopore is a function of the pore size. Here we focus on the modulus-pressure dependence K(P), which is linear for bulk materials, a relation known as the Tait-Murnaghan equation. Using transition-matrix Monte Carlo simulations we calculated the elastic modulus of bulk argon as a function of pressure and argon confined in silica mesopores as a function of Laplace pressure. Our calculations show that while the elastic modulus is strongly affected by confinement and temperature, the slope of the modulus versus pressure is not. Moreover, the calculated slope is in a good agreement with the reference data for bulk argon and experimental data for confined argon derived from ultrasonic experiments. We propose to use the value of the slope of K(P) to estimate the elastic moduli of an unknown porous medium.