Communication: Analytic continuation of the virial series through the critical point using parametric approximants

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.

Analytic solution of the SEIR epidemic model via asymptotic approximant

An analytic solution is obtained to the SEIR Epidemic Model. The solution is created by constructing a single second-order nonlinear differential equation in ln S and analytically continuing its divergent power series solution such that it matches the correct long-time exponential damping of the epidemic model. This is achieved through an asymptotic approximant (Barlow et al., 2017) in the form of a modified symmetric Padé approximant that incorporates this damping. The utility of the analytical form is demonstrated through its application to the COVID-19 pandemic.

Accurate closed-form solution of the SIR epidemic model

An accurate closed-form solution is obtained to the SIR Epidemic Model through the use of Asymptotic Approximants (Barlow et al., 2017). The solution is created by analytically continuing the divergent power series solution such that it matches the long-time asymptotic behavior of the epidemic model. The utility of the analytical form is demonstrated through its application to the COVID-19 pandemic.

Cluster integrals and virial coefficients for realistic molecular models

We present a concise, general, and efficient procedure for calculating the cluster integrals that relate thermodynamic virial coefficients to molecular interactions. The approach encompasses nonpairwise intermolecular potentials generated from quantum chemistry or other sources; a simple extension permits efficient evaluation of temperature and other derivatives of the virial coefficients. We demonstrate with a polarizable model of water. We argue that cluster-integral methods are a potent yet underutilized instrument for the development and application of first-principles molecular models and methods.

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.