Optimized theory for simple and molecular fluids

An optimized closure approximation for both simple and molecular fluids is presented. A smooth interpolation between Perkus-Yevick and hypernetted chain closures is optimized by minimizing the free energy self-consistently with respect to the interpolation parameter(s). The molecular version is derived from a refinement of the method for simple fluids. In doing so, a method is proposed which appropriately couples an optimized closure with the variant of the diagrammatically proper integral equation recently introduced by this laboratory [K. M. Dyer et al., J. Chem. Phys. 123, 204512 (2005)]. The simplicity of the expressions involved in this proposed theory has allowed the authors to obtain an analytic expression for the approximate excess chemical potential. This is shown to be an efficient tool to estimate, from first principles, the numerical value of the interpolation parameters defining the aforementioned closure. As a preliminary test, representative models for simple fluids and homonuclear diatomic Lennard-Jones fluids were analyzed, obtaining site-site correlation functions in excellent agreement with simulation data.     

Perturbation theory without wave function for multidimensional systems

A new method to obtain perturbation corrections to the eigenvalues of multidimensional quantummechanical models is developed. It consists of rearranging the Rayleigh-Schrodinger perturbation theory so that any coefficient of the perturbation series is obtained from a simple and compact recursion relationship. The Zeeman effect in hydrogen and the hydrogen molecule-ion are used to illustrate the procedure.     

On the limiting radial distribution function for hydrogenic orbitals

An exact reduced limiting expression for the generalized radial distribution functionD _n (r) is derived and compared with quantum distributions for various degrees of excitation. It represents the quantum result at large quantum numbers significantly better than a prior empirical representation of the universal reduced distribution and gives a somewhat larger electronic partition function for the hydrogen atom than that based on the previous distribution.     

Hard-sphere radial distribution function again

A theoretically based closed-form analytical equation for the radial distribution function, g(r), of a fluid of hard spheres is presented and used to obtain an accurate analytic representation. The method makes use of an analytic expression for the short- and long-range behaviors of g(r), both obtained from the Percus-Yevick equation, in combination with the thermodynamic consistency constraint. Physical arguments then leave only three parameters in the equation of g(r) that are to be solved numerically, whereas all remaining ones are taken from the analytical solution of the Percus-Yevick equation.     

Approximative "one particle" bridge function B(1)(r) for the theory of simple fluids

New properties for the one particle bridge function B(1)(r), which are necessary to the calculation of the excess chemical potential betamue), are derived for the hard sphere fluid. The method, which only requires the knowledge of the bridge function B(2)(r), is based on an investigation of the correlation function dependence on the Kirkwood charging parameter. In this framework, the unavoidable question of topological homotopy is addressed. As far as B(2)(r) is considered as exact, this work provides useful information on B(1)(r) in the well identified dynamical regimes of the hard sphere fluid. Signatures of the transitions between these regimes are identified on the trends of B(1)(r). This approach provides self-consistent results for betamue) that agree very well with simulation data.     

An accurate expression for radial distribution function of the Lennard-Jones fluid

A simple and accurate expression for radial distribution function (RDF) of the Lennard-Jones fluid is presented. The expression explicitly states the RDF as a continuous function of reduced interparticle distance, temperature, and density. It satisfies the limiting conditions of zero density and infinite distance imposed by statistical thermodynamics. The distance dependence of this expression is expressed by an equation which contains 11 adjustable parameters. These parameters are fitted to 353 RDF data, obtained by molecular dynamics calculations, and then expressed as functions of reduced distance, temperature and density. This expression, having a total of 65 constants, reproduces the RDF data with an average root-mean-squared deviation of 0.0152 for the range of state variables of 0.5  T^*  5.1 and 0.35  ρ^*  1.1 (T^*=kT/ε and ρ^* = ρσ³ are reduced temperature and density, respectively). The expression predicts the pressure and the internal energy of the Lennard-Jones fluid with an uncertainty that is comparable to that obtained directly from the molecular dynamics simulations.     

Second-order perturbation theory for the angular pair correlation function in molecular fluids

Second-order perturbation theory is used to calculate spherical harmonic coefficients of the angular pair correlation function g(rω₁ω₂) for a liquid in which the molecules interact with a pair potential that is the sum of Lennard-Jones and quadrupole-quadrupole parts. The theory is compared with both molecular dynamics results and with the predictions of the GMF ≡ LHNC, QHNC and first-order perturbation theories. Second-order perturbation theory gives excellent results for the harmonic coefficient g(224,r), but is poorer for g(222,r) and g(202,r).     

Perturbation theory for confined systems

We discuss the application of perturbation theory to a system of particles confined in a spherical box. A simple argument shows that the particles behave almost independently in sufficiently strong confinement. We choose the helium atom with a moving nucleus as a particular example and compare results of first order with those for the nucleus clamped at the center of the box. We provide a suitable explanation for some numerical results obtained recently by other authors.     

Generalized Langevin theory on the dynamics of simple fluids under external fields

A theory on the time development of the density and current fields of simple fluids under an external field is formulated through the generalized Langevin formalism. The theory is applied to the linear solvation dynamics of a fixed solute regarding the solute as the external field on the solvent. The solute-solvent-solvent three-body correlation function is taken into account through the hypernetted-chain integral equation theory, and the time correlation function of the random force is approximated by that in the absence of the solute. The theoretical results are compared with those of molecular-dynamics (MD) simulation and the surrogate theory. As for the transient response of the density field, our theory is shown to be free from the artifact of the surrogate theory that the solvent can penetrate into the repulsive core of the solute during the relaxation. We have also found a large quantitative improvement of the solvation correlation function compared with the surrogate theory. In particular, the short-time part of the solvation correlation function is in almost perfect agreement with that from the MD simulation, reflecting that the short-time expansion of the theoretical solvation correlation function is exact up to t(2) with the exact three-body correlation function. A quantitative improvement is found in the long-time region, too. Our theory is also applied to the force-force time correlation function of a fixed solute, and similar improvement is obtained, which suggests that our present theory can be a basis to improve the mode-coupling theory on the solute diffusion.     

Perturbation theory for multiconfiguration reference states

Vector method procedures are adapted to evaluate Rayleigh-Schrödinger perturbation corrections to a multiconfiguration zeroth order function. If this function is sufficiently flexible, this perturbation theory can be applied to low lying excited states. The effectiveness of our theory is demonstrated on the ground state of F₂ and the low lying excited states of Mg₂. Energies calculated through fourth order are compared with appropriate CI results.     

Symmetrical radial distribution functions in the potentially theory of electrolyte solutions

A symmetrical radial distribution function g_ij in electric double layer theory has been proposed by Féat and Levine. We adopt their work to the electrolyte bulk to show how a symmetrical g_ij may be obtained in the potential theory of electrolyte solutions.     

Parametrization of the radial distribution function of a Yukawa fluid

We give semi-empirical expressions for five terms in an inverse temperature expansion of the radial distribution function of a Yukawa fluid. The parameters in these terms are chosen to fit what we believe should be an accurate equation of state when either the energy or pressure routes is used. Thus, a measure of self consistency is achieved. The equation of state which is the basis of our fit is the inverse temperature expansion of the Yukawa fluid free energy, obtained from the mean spherical approximation but modified to give reasonable results at low densities.     

An extended rism equation for molecular polar fluids

The RISM integral equation is extended to molecules with charged sites via a renormalization of the Coulomb potentials and the introduction of appropriate closure relations. For a fluid of diatomics with atomic charges of ±0.2 e the equation yields site-site correlation functions in qualitative agreement with those from computer simulation.     

Blip function calculation of the radial distribution functions for a binary liquid mixture

It is shown that the use of the high temperature approximation to blip function theory together with a two-component hard sphere reference fluid leads to reasonably accurate predictions for the radial distribution functions in a moderate temperature, moderate density, binary Lennard-Jones mixture.     

A simple molecular theory of double reentrance exhibited by highly polar compounds

We show by simple model calculations that dipole-induced dipole and interchain interactions lead to a parallel alignment of molecules with strong longitudinal dipoles for intermolecular separations below a certain value. On the other hand an antiparallel alignment is favoured for larger separations because of dipole-dipole interactions. We incorporate this change in the intermolecular configuration, which naturally leads to two 'lengths' in the problem, in developing a simple McMillan type molecular theory of the double reentrance phenomenon in which the following phase sequence is observed on cooling the sample: isotropic-nematic-smectic A_d-reentrant nematic-smectic A₁. The calculated properties including the phase diagram are in broad agreement with experimental trends.