Fifth to eleventh virial coefficients of hard spheres

Origins of the Failure of the Activity Virial Series

We examine the development of the virial equation of state when expressed as a series in the activity with coefficients labeled b_n. Using the one-dimensional hard-rod model as a prototype, we consider steps in its development that introduce inaccuracy that leads it to form a divergent series. We discuss the role of volume dependence of the virial coefficients and present expressions and calculations for volume-dependent coefficients b_n(V) for the hard-rod model up to n = 200. We examine alternative methods for computing properties from the b_n. We recommend that further efforts be made to compute volume-dependent virial coefficients as a means to understand better the virial equation of state and to make it more robust in applications.

Virial coefficients of hard, homonuclear dumbbells in two- to four-dimensional Euclidean spaces

We calculated virial coefficients up to the eighth order for hard dumbbells in two-, three-, and four-dimensional Euclidean spaces employing Mayer-sampling Monte Carlo simulations. We improved and extended available data in two dimensions, provide virial coefficients in R^{4} in dependence on their aspect ratio, and recalculated virial coefficients for three-dimensional dumbbells. Highly accurate, semianalytical values for the second virial coefficient of homonuclear, four-dimensional dumbbells are provided. We compare the influence of the aspect ratio and the dimensionality to the virial series for this concave geometry. Lower-order reduced virial coefficients B[over ̃]_{i}=B_{i}/B_{2}^{i-1} depend in first approximation linearly from the inverse excess part of their mutual excluded volume.

Virial coefficients of hard hyperspherocylinders in R^{4}: Influence of the aspect ratio

We provide second- to sixth-order virial coefficients of hard hyperspherocylinders in dependence on their aspect ratio ν. Virial coefficients of an anisotropic geometry in four dimensions are calculated employing an optimized Mayer-sampling algorithm. As the second virial coefficient of a hard particle is identical to its excluded hypervolume, the numerically obtained second virial coefficients can be compared to analytical relations for the excluded hypervolume based on geometric measures of the respective, convex geometry in dependence on its aspect ratio ν.

Calculation of third to eighth virial coefficients of hard lenses and hard, oblate ellipsoids of revolution employing an efficient algorithm

We provide third to eighth virial coefficients of oblate, hard ellipsoids of revolution and hard lenses in dependence on their aspect ratio ν. Employing an algorithm optimized for hard anisotropic shapes, highly accurate data are accessible with comparatively small numerical effort. For both geometries, reduced virial coefficients B[over ̃]_{i}(ν)=B_{i}(ν)/B_{2}^{i-1}(ν) are in first approximation proportional to the inverse excess contribution α^{-1} of their excluded volume. The latter quantity is directly accessible from second virial coefficients and analytically known for convex bodies.

Evaluation of Osmotic Virial Coefficients via Restricted Gibbs Ensemble Simulations, with Support from Gas-Phase Mixture Coefficients

We present a method for computing osmotic virial coefficients in explicit solvent via simulation in a restricted Gibbs ensemble. Two equivalent phases are simulated at once, each in a separate box at constant volume and temperature and each in equilibrium with a solvent reservoir. For osmotic coefficient BN, a total of N solutes are individually exchanged back and forth between the boxes, and the average distribution of solute numbers between the boxes provides the key information needed to compute BN. Separately, expressions are developed for BN as a series in solvent reservoir density ρ₁, with the coefficients of the series expressed in terms of the usual gas-phase mixture coefficients B_ij. Normally, the B_ij are defined for an infinite volume, but we suggest that the observed dependence of B_ij on system size L can be used to estimate L dependence of the BN, allowing them to be computed accurately at L → ∞ while simulating much smaller system sizes than otherwise possible. The methods for N = 2 and 3 are demonstrated for two-component mixtures of size-asymmetric additive hard spheres. The proposed methods are demonstrated to have greater precision than established techniques, for a given amount of computational effort. The ρ₁ series for BN when applied by itself is (for this noncondensing model) found to be the most efficient in computing accurate osmotic coefficients for the solvent densities considered here.

Algebraic second virial coefficient of the Mie m - 6 intermolecular potential based on perturbation theory

We propose several simple algebraic approximations for the second virial coefficient of fluids whose molecules interact by a generic Mie m - 6 intermolecular pair potential. In line with a perturbation theory, the parametric equations are formulated as the sum of a contribution due to a reference part of the intermolecular potential and a perturbation. Thereby, the equations provide a convenient (low-density) starting point for developing equation-of-state models of fluids or for developing similar approximations for the virial coefficient of (polymeric-)chain fluids. The choice of Barker and Henderson [J. Chem. Phys. 47, 4714 (1967)] and Weeks, Chandler, and Andersen [Phys. Rev. Lett. 25, 149 (1970); J. Chem. Phys. 54, 5237 (1971); and Phys. Rev. A 4, 1597 (1971)] for the reference part of the potential is considered. Our analytic approximations correctly recover the virial coefficient of the inverse-power potential of exponent m in the high-temperature limit and provide accurate estimates of the temperatures for which the virial coefficient equals zero or takes on its maximum value. Our description of the reference contribution to the second virial coefficient follows from an exact mapping onto the second virial coefficient of hard spheres; we propose a simple algebraic equation for the corresponding effective diameter of the hard spheres, which correctly recovers the low- and high-temperature scaling and limits of the reference fluid's second virial coefficient.

Toward a Suite of Coarse-Grained Models for Molecular Simulation of Monoclonal Antibodies and Therapeutic Proteins

A series of coarse-grained models for molecular simulation of proteins are considered, with emphasis on the application of predicting protein-protein self-interactions for monoclonal antibodies (MAbs). As an illustrative example and for quantitative comparison, the models are used to predict osmotic virial coefficients over a broad range of attractive and repulsive self-interactions and solution conditions for a series of MAbs where the second osmotic virial coefficient has been experimentally determined in prior work. The models are compared based on how well they can predict experimental behavior, their computational burdens, and scalability. An intermediate-resolution model is also introduced that can capture specific electrostatic interactions with improved efficiency and similar or improved accuracy when compared to the previously published models. Guidance is included for the selection of coarse-grained models more generally for capturing a balance of electrostatic, steric, and short-ranged nonelectrostatic interactions for proteins from low to high concentrations.

Performance of the asymptotic expansion method to derive equations of state for hard polyhedron fluids

The asymptotic expansion method is used to derive analytical expressions for the equations of state of 14 hard polyhedron fluids such as cube, octahedron, rhombic dodecahedron, etc., by knowing the values of only the first eight virial coefficients.

Predicting maximally random jammed packing density of non-spherical hard particles via analytical continuation of fluid equation of state

We developed a formalism for accurately predicting the density of MRJ packing state of a wide spectrum of congruent non-spherical hard particles in 3D via analytical fluid EOS.

Predicting High-Concentration Interactions of Monoclonal Antibody Solutions: Comparison of Theoretical Approaches for Strongly Attractive Versus Repulsive Conditions

Nonspecific protein-protein interactions of a monoclonal antibody were quantified experimentally using light scattering from low to high protein concentrations (c₂) and compared with prior work for a different antibody that yielded qualitatively different behavior. The c₂ dependence of the excess Rayleigh ratio (R^ex) provided the osmotic second virial coefficient (B₂₂) at low c₂ and the static structure factor (S_q=0) at high c₂, as a function of solution pH, total ionic strength (TIS), and sucrose concentration. Net repulsive interactions were observed at pH 5, with weaker repulsions at higher TIS. Conversely, attractive electrostatic interactions were observed at pH 6.5, with weaker attractions at higher TIS. Refined coarse-grained models were used to fit model parameters using experimental B₂₂ versus TIS data. The parameters were used to predict high-c₂ R^ex values via Monte Carlo simulations and separately with Mayer-sampling calculations of higher-order virial coefficients. For both methods, predictions for repulsive to mildly attractive conditions were quantitatively accurate. However, only qualitatively accurate predictions were practical for strongly attractive conditions. An alternative, higher resolution model was used to show semiquantitatively and quantitatively accurate predictions of strong electrostatic attractions at low c₂ and low ionic strength.

Thermodynamic and dynamical properties of the hard sphere system revisited by molecular dynamics simulation

Revised thermodynamic and dynamical properties of the hard sphere (HS) system are obtained from extensive molecular dynamics calculations carried out with large system sizes (number of particles, N) and long times. Accurate formulas for the compressibility factor of the HS solid and fluid branches are proposed, which represent the metastable region and take into account its divergence at close packing. Some basic second-order thermodynamic properties are obtained and a maximum in some of their derivatives in the metastable fluid region is found. The thermodynamic parameters associated with the melting-freezing transition have been determined to four digit accuracy, which generates accurate new values for the coexistence properties of the HS system. For the self-diffusion coefficient, D, it is shown that relatively large systems (N > 104) are required to achieve an accurate linear extrapolation of D to the infinite size limit with a D vs. N-1/3 plot. Moreover, it is found that there is a density dependence of the value of the slope in the linear regime. The density dependent correction becomes practically insignificant at higher densities and the hydrodynamic formula found in the literature is still accurate. However, with decreasing density the density dependence of the size correction cannot be neglected, which indicates that other sources of N-dependence, apart from those derived on purely hydrodynamic grounds, may also be important (and as yet unaccounted for). A detailed analytic representation of the density dependence of the HS self-diffusion coefficient and the HS viscosity, η, is given. It is shown that the HS viscosity near freezing and in the metastable region can be described well by the Krieger-Dougherty equation. Both D and η start to scale at high densities and in the metastable region in such a way that Dηp = const, where p ≃ 0.97, and D → 0 and η → ∞ at a packing fraction of 0.58, which coincides with some previous predictions of the HS glass transition density.

Equations of the state of hard sphere fluids based on recent accurate virial coefficients B5–B12

A review on the numerical virial coefficients, compressibility factor, fluid–solid phase transition point and equations of the state of hard sphere fluids.

Virial Coefficients and Equations of State for Hard Polyhedron Fluids

Hard polyhedra are a natural extension of the hard sphere model for simple fluids, but there is no general scheme for predicting the effect of shape on thermodynamic properties, even in moderate-density fluids. Only the second virial coefficient is known analytically for general convex shapes, so higher-order equations of state have been elusive. Here we investigate high-precision state functions in the fluid phase of 14 representative polyhedra with different assembly behaviors. We discuss historic efforts in analytically approximating virial coefficients up to B₄ and numerically evaluating them to B₈. Using virial coefficients as inputs, we show the convergence properties for four equations of state for hard convex bodies. In particular, the exponential approximant of Barlow et al. (J. Chem. Phys. 2012, 137, 204102) is found to be useful up to the first ordering transition for most polyhedra. The convergence behavior we explore can guide choices in expending additional resources for improved estimates. Fluids of arbitrary hard convex bodies are too complicated to be described in a general way at high densities, so the high-precision state data we provide can serve as a reference for future work in calculating state data or as a basis for thermodynamic integration.

The Uhlenbeck-Ford model: Exact virial coefficients and application as a reference system in fluid-phase free-energy calculations

The Uhlenbeck-Ford (UF) model was originally proposed for the theoretical study of imperfect gases, given that all its virial coefficients can be evaluated exactly, in principle. Here, in addition to computing the previously unknown coefficients B₁₁ through B₁₃, we assess its applicability as a reference system in fluid-phase free-energy calculations using molecular simulation techniques. Our results demonstrate that, although the UF model itself is too soft, appropriately scaled Uhlenbeck-Ford (sUF) models provide robust reference systems that allow accurate fluid-phase free-energy calculations without the need for an intermediate reference model. Indeed, in addition to the accuracy with which their free energies are known and their convenient scaling properties, the fluid is the only thermodynamically stable phase for a wide range of sUF models. This set of favorable properties may potentially put the sUF fluid-phase reference systems on par with the standard role that harmonic and Einstein solids play as reference systems for solid-phase free-energy calculations.

Calculation of high-order virial coefficients for the square-well potential

Accurate virial coefficients B_{N}(λ,ɛ) (where ɛ is the well depth) for the three-dimensional square-well and square-step potentials are calculated for orders N=5-9 and well widths λ=1.1-2.0 using a very fast recursive method. The efficiency of the algorithm is enhanced significantly by exploiting permutation symmetry and by storing integrands for reuse during the calculation. For N=9 the storage requirements become sufficiently large that a parallel algorithm is developed. The methodology is general and is applicable to other discrete potentials. The computed coefficients are precise even near the critical temperature, and thus open up possibilities for analysis of criticality of the system, which is currently not accessible by any other means.

Computation of virial coefficients from integral equations

A polynomial-time method of computing the virial coefficients from an integral equation framework is presented. The method computes the truncated density expansions of the correlation functions by series transformations, and then extracts the virial coefficients from the density components. As an application, the method was used in a hybrid-closure integral equation with a set of self-consistent conditions, which produced reasonably accurate virial coefficients for the hard-sphere fluid and Gaussian model in high dimensions.