Effective density terms in proper integral equations

Two complementary routes to a new integral equation theory for site-site molecular fluids are presented. First, a simple approximation to a subset of the atomic site bridge functions in the diagrammatically proper integral equation theory is presented. This in turn leads to a form analogous to the reactive fluid theory, in which the normalization of the intramolecular distribution function and the value of the off-diagonal elements in the density matrix of the proper integral equations are the means of propagating the bridge function approximation. Second, a derivation from a topological expansion of a model for the single-site activity followed by a topological reduction and low-order truncation is given. This leads to an approximate numerical value for the new density coefficient. The resulting equations give a substantial improvement over the standard construction as shown with a series of simple diatomic model calculations.     

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.     

Towards a more accurate reference interaction site model integral equation theory for molecular liquids

A systematic approach for increasing the accuracy of the reference interaction site model (RISM) theory is introduced that uses input from simulation results to produce very accurate site-site pair correlation functions for single component molecular liquids. The methodology allows the computation of the "RISM bridge function." Realistic molecular liquids such as water, alcohols, amides, and others are investigated, and the merits and limitations of the method for each of these liquids are examined in relation to the known deficiencies of the RISM theory.     

Self-consistent field equations in the distribution function theory of polymeric liquids

In the distribution function approach to the conformational and thermodynamic properties of polymeric liquids site-site (pair) distribution functions are essential components of the theory. These site-site pair distribution functions are basically mean fields obeying integral equations. In our recent works, a set of self-consistent field equations has been proposed for site-site pair correlation functions which allow us to study conformational and thermodynamic properties of polymeric liquids. In this article, we present a short review of the theory and its applications to a number of aspects of polymeric liquids we have made until now. We also present a self-consistent version of the polymer reference interaction site model where the integral equations for the intramolecular site-site correlation functions are obtained from the Kirkwood hierarchy on the basis of the present theory. The present theory is shown to predict correctly the scaling properties associated with swollen and collapsed polymers in good and poor solvents, respectively. At finite densities, self-consistent solutions of the intra- and intermolecular equations yield the structures and thermodynamics of polymer melts which are favorably compared with Monte Carlo simulation results. Self-consistent theory results are found to be more accurate than the non-self-consistent approaches that use an ideal Gaussian chain conformation distribution function. © 1995 John Wiley & Sons, Inc.     

An interaction site model integral equation study of molecular fluids explicitly considering the molecular orientation

We implemented an interaction site model integral equation for rigid molecules based on a density-functional theory where the molecular orientation is explicitly considered. In this implementation of the integral equation, multiple integral of the degree of freedom of the molecular orientation is performed using efficient quadrature methods, so that the site-site pair correlation functions are evaluated exactly in the limit of low density. We apply this method to Cl(2), HCl, and H(2)O molecular fluids that have been investigated by several integral equation studies using various models. The site-site pair correlation functions obtained from the integral equation are in good agreement with the one from a simulation of these molecules. Rotational invariant coefficients, which characterize the microscopic structure of molecular fluids, are determined from the integral equation and the simulation in order to investigate the accuracy of the integral equation.     

Second virial coefficients of asymmetric top molecules

A short self-contained derivation is given for the second virial coefficient B2(T) of a gas consisting of identical interacting asymmetric rigid rotors. The resulting expression is correct through variant Planck's over h2. First, the canonical partition function is derived by means of an variant Planck's over h expansion of exp[-H/(k(B)T)] due to Friedmann [Adv. Chem. Phys. 4, 225 (1962)]. The present work applies angular momentum operators and known facts from angular momentum theory. It is considerably more accessible than Friedmann's exposition, which is not based on angular momentum operators, but instead on explicit derivatives with respect to Euler angles. The partition function obtained from the variant Planck's over h expansion is applied to the derivation of an expression for B2(T) that is identical in appearance to the expression for symmetric rotors of T Pack [J. Chem. Phys. 78, 7217 (1983)]. The final equation in this work is valid for rigid rotors of any symmetry.     

On the exchange-hole model of London dispersion forces

First-principles derivation is given for the heuristic exchange-hole model of London dispersion forces by Becke and Johnson [J. Chem. Phys. 122, 154104 (2005)]. A one-term approximation is used for the dynamic charge density response function, and it is shown that a central nonempirical ingredient of the approximate nonexpanded dispersion energy is the charge density autocorrelation function, a two-particle property, related to the exchange-correlation hole. In the framework of a dipolar approximation of the Coulomb interaction around the molecular origin, one obtains the so-called Salem-Tang-Karplus approximation to the C(6) dispersion coefficient. Alternatively, by expanding the Coulomb interaction around the center of charge (centroid) of the exchange-correlation hole associated with each point in the molecular volume, a multicenter expansion is obtained around the centroids of electron localization domains, always in terms of the exchange-correlation hole. In order to get a formula analogous to that of Becke and Johnson, which involves the exchange-hole only, further assumptions are needed, related to the difficulties of obtaining the expectation value of a two-electron operator from a single determinant. Thus a connection could be established between the conventional fluctuating charge density model of London dispersion forces and the notion of the "exchange-hole dipole moment" shedding some light on the true nature of the approximations implicit in the Becke-Johnson model.     

Theoretical and computational investigations on thermodynamic properties, effective site diameters, and molecular free volume of carbon disulfide fluid

A newly proposed theory [R. Laghaei et al., J. Chem. Phys. 124, 154502 (2006)] was extended to polyatomics and applied to compute the density and temperature dependence of the effective site diameters of carbon disulfide fluids. The generic van der Waals (GvdW) theory was also extended to polyatomics in order to calculate the GvdW parameters and the molecular free volume using the effective site diameters as the repulsion-attraction separation distance. A three-site Lennard-Jones potential available in the literature was slightly modified and used in Monte Carlo simulations to obtain the functions appearing in the effective site diameter and GvdW expressions. The interaction potential was examined to reproduce the fluid phase thermodynamic properties using Gibbs ensemble Monte Carlo simulations and also the equation of state in the liquid phase using NVT Monte Carlo (NVT-MC) simulations. Comparison between the simulation results and experimental data shows excellent agreement for the densities of the coexisting phases, the vapor pressure, properties of the predicted critical point, and the equation of state. NVT-MC simulations were performed over a wide range of densities and temperatures in sub- and supercritical regions to compute the effective site diameters, the GvdW parameters, and the molecular free volume. The molecular structure in terms of the site-site pair correlation functions, the density dependence of the effective site diameters, and the density and temperature dependence of the GvdW parameters and molecular free volume were studied and discussed. The GvdW parameters were fitted to empirical expressions as a function of density and temperature. The computed molecular free volume will be used in future investigations to study the transport properties of carbon disulfide.     

Wavelet algorithm for solving integral equations of molecular liquids. A test for the reference interaction site model

A new efficient method is developed for solving integral equations based on the reference interaction site model (RISM) of molecular liquids. The method proposes the expansion of site-site correlation functions into the wavelet series and further calculations of the approximating coefficients. To solve the integral equations we have applied the hybrid scheme in which the coarse part of the solution is calculated by wavelets with the use of the Newton-Raphson procedure, while the fine part is evaluated by the direct iterations. The Coifman 2 basis set is employed for the wavelet treatment of the coarse solution. This wavelet basis set provides compact and accurate approximation of site-site correlation functions so that the number of basis functions and the amplitude of the fine part of solution decrease sufficiently with respect to those obtained by the conventional scheme. The efficiency of the method is tested by calculations of SPC/E model of water. The results indicated that the total CPU time to obtain solution by the proposed procedure reduces to 20% of that required for the conventional hybrid method.     

An analysis of the solution of groups method for component activity coefficients

The solution of groups method is analyzed via a statistical mechanical formalism in order to uncover its inherent physical assumptions. Viewing the molecules in solution as a physical ‘cluster’ of groups, or interaction sites, provides a framework through which the implicit assumptions may be cast in standard molecular terminology. We find that, in addition to the assumed site-site form of the molecular pair potential, asymmetric distribution of sites within each molecule and use of an ‘effective’ angle average of the pair correlation function are also necessary in order to arrive at the commonly used solution of groups formula. Alternatively, the latter approximations can be interpreted as the site analogue of the n-fluid theory for simple mixtures. Attempts to improve upon the present empirical models may be made by relaxing some of the more restrictive assumptions.     

Modified free volume theory of self-diffusion and molecular theory of shear viscosity of liquid carbon dioxide

In previous work on the density fluctuation theory of transport coefficients of liquids, it was necessary to use empirical self-diffusion coefficients to calculate the transport coefficients (e.g., shear viscosity of carbon dioxide). In this work, the necessity of empirical input of the self-diffusion coefficients in the calculation of shear viscosity is removed, and the theory is thus made a self-contained molecular theory of transport coefficients of liquids, albeit it contains an empirical parameter in the subcritical regime. The required self-diffusion coefficients of liquid carbon dioxide are calculated by using the modified free volume theory for which the generic van der Waals equation of state and Monte Carlo simulations are combined to accurately compute the mean free volume by means of statistical mechanics. They have been computed as a function of density along four different isotherms and isobars. A Lennard-Jones site-site interaction potential was used to model the molecular carbon dioxide interaction. The density and temperature dependence of the theoretical self-diffusion coefficients are shown to be in excellent agreement with experimental data when the minimum critical free volume is identified with the molecular volume. The self-diffusion coefficients thus computed are then used to compute the density and temperature dependence of the shear viscosity of liquid carbon dioxide by employing the density fluctuation theory formula for shear viscosity as reported in an earlier paper (J. Chem. Phys. 2000, 112, 7118). The theoretical shear viscosity is shown to be robust and yields excellent density and temperature dependence for carbon dioxide. The pair correlation function appearing in the theory has been computed by Monte Carlo simulations.     

Parameterized local hybrid functionals from density-matrix similarity metrics

We recently proposed a real-space similarity metric comparing the Kohn-Sham one-particle density matrix to the local spin-density approximation model density matrix [Janesko and Scuseria, J. Chem. Phys. 127, 164117 (2007)]. This metric provides a useful ingredient for constructing local hybrid density functionals that locally mix exact exchange and semilocal density functional theory exchange. Here we present two lines of inquiry: An approximate similarity metric comparing exact versus generalized gradient approximation (GGA), exchange and parameterized mixing functions using these similarity metrics. This approach yields significantly improved thermochemistry, including GGA local hybrids whose thermochemical performance approaches GGA global hybrids.     

Exactly Solvable Model Mimicking the H<Subscript>2</Subscript> Molecule in the Limit of Large Nuclear Masses

Avoiding the limitations of the Born–Oppenheimer approximation remains an issue of continuing importance in molecular physics. If we take the H₂ molecule as a prototype, then one uses, in some sense, the ratio of electron mass m to nuclear mass M as an expansion parameter, sometimes cited to be explicitly (m/M)^1/4. Here, we employ a model of such a two-electron diatomic molecule set up by Makarewicz [Am. J. Phys. 54 (1986) 178] to study the exact ground-state wave function obtained by treating electrons and nuclei on the same footing, Of course, to obtain such a wave function analytically means adopting model force laws, both for confinement and for the interelectronic interaction. From the exact wave function and ground-state energy, contact is finally established with the large nuclear mass limit.     

Reactive and nonreactive scattering of N2 from Ru(0001): a six-dimensional adiabatic study

We have studied the dissociative chemisorption and scattering of N(2) on and from Ru(0001), using a six-dimensional quasiclassical trajectory method. The potential energy surface, which depends on all the molecular degrees of freedom, has been built applying a modified Shepard interpolation method to a data set of results from density functional theory, employing the RPBE generalized gradient approximation. The frozen surface and Born-Oppenheimer [Ann. Phys. (Leipzig) 84, 457 (1927)] approximations were used, neglecting phonons and electron-hole pair excitations. Dissociative chemisorption probabilities are found to be very small even for translational energies much higher than the minimum reaction barrier, in good agreement with experiment. A comparison to previous low dimensional calculations shows the importance of taking into account the multidimensional effects of N(2) rotation and translation parallel to the surface. The new calculations strongly suggest a much smaller role of nonadiabatic effects than previously assumed on the basis of a comparison between low dimensional results and experiments [J. Chem. Phys. 115, 9028 (2001)]. Also in agreement with experiment, our theoretical results show a strong dependence of reaction on the initial vibrational state. Computed angular scattering distributions and parallel translation energy distributions are in good agreement with experiments on scattering, but the theory overestimates vibrational and rotational excitations in scattering.     

An integral equation theory for inhomogeneous molecular fluids: the reference interaction site model approach

An integral equation theory which is applicable to inhomogeneous molecular liquids is proposed. The "inhomogeneous reference interaction site model (RISM)" equation derived here is a natural extension of the RISM equation to inhomogeneous systems. This theory makes it possible to calculate the pair correlation function between two molecules which are located at different density regions. We also propose approximations concerning the closure relation and the intramolecular susceptibility of inhomogeneous molecular liquids. As a preliminary application of the theory, the hydration structure around an ion is investigated. Lithium, sodium, and potassium cations are chosen as the solute. Using the Percus trick, the local density of solvent around an ion is expressed in terms of the solute-solvent pair correlation function calculated from the RISM theory. We then analyze the hydration structure around an ion through the triplet correlation function which is defined with the inhomogeneous pair correlation function and the local density of the solvent. The results of the triplet correlation functions for cations indicate that the thermal fluctuation of the hydration shell is closely related to the size of the solute ion. The triplet correlation function from the present theory is also compared with that from the Kirkwood superposition approximation, which substitutes the inhomogeneous pair correlation by the homogeneous one. For the lithium ion, the behavior of the triplet correlation functions from the present theory shows marked differences from the one calculated within the Kirkwood approximation.     

Site-site memory equation approach in study of density/pressure dependence of translational diffusion coefficient and rotational relaxation time of polar molecular solutions: acetonitrile in water, methanol in water, and methanol in acetonitrile

We present results of the theoretical study and numerical calculation of the dynamics of molecular liquids based on the combination of the memory equation formalism and the reference interaction site model (RISM). Memory equations for the site-site intermediate scattering functions are studied in the mode-coupling approximation for the first-order memory kernels, while equilibrium properties such as site-site static structure factors are deduced from RISM. The results include the temperature-density (pressure) dependence of translational diffusion coefficients D and orientational relaxation times tau for acetonitrile in water, methanol in water, and methanol in acetonitrile--all in the limit of infinite dilution. Calculations are performed over the range of temperatures and densities employing the extended simple point charge model for water and optimized site-site potentials for acetonitrile and methanol. The theory is able to reproduce qualitatively all main features of temperature and density dependences of D and tau observed in real and computer experiments. In particular, anomalous behavior, i.e, the increase in mobility with density, is observed for D and tau of methanol in water, while acetonitrile in water and methanol in acetonitrile do not show deviations from the ordinary behavior. The variety exhibited by the different solute-solvent systems in the density dependence of the mobility is interpreted in terms of the two competing origins of friction, which interplay with each other as density increases: the collisional and dielectric frictions which, respectively, increase and decrease with increasing density.     

Improvement on macroscopic compressibility approximation and beyond

A numerical procedure is proposed to extend the thermodynamic perturbation expansion (TPE) to a higher order. It is shown that the present second order term is superior to that due to a macroscopic compressibility approximation (MCA), a local compressibility approximation, and a superposition approximation by Barker and Henderson [Rev. Mod. Phys. 48, 587 (1976)]. Extensive model calculation and comparison with simulation data available in literature and supplied in the present report indicate that the present third order TPE is superior to a previous second order TPE based on the MCA, two previous perturbation theories, which are respectively based on an analytical mean spherical approximation for an Ornstein-Zernike equation, and an assumed explicit functional form for the Laplace transform of radial distribution function multiplied by radial distance, and a recent generalized van der Waals theory. The present critical temperature for a hard core attractive Yukawa fluid of varying range is in very good agreement with that due to a hierarchical reference theory. The present third order TPE is computationally far more modest than the self-consistent integral equation theory, and therefore is a viable alternative to use of the latter.     

Quantum theory of molecular collisions in a magnetic field: efficient calculations based on the total angular momentum representation

An efficient method is presented for rigorous quantum calculations of atom-molecule and molecule-molecule collisions in a magnetic field. The method is based on the expansion of the wave function of the collision complex in basis functions with well-defined total angular momentum in the body-fixed coordinate frame. We outline the general theory of the method for collisions of diatomic molecules in the (2)Σ and (3)Σ electronic states with structureless atoms and with unlike (2)Σ and (3)Σ molecules. The cross sections for elastic scattering and Zeeman relaxation in low-temperature collisions of CaH((2)Σ(+)) and NH((3)Σ(-)) molecules with (3)He atoms converge quickly with respect to the number of total angular momentum states included in the basis set, leading to a dramatic (>10-fold) enhancement in computational efficiency compared to the previously used methods [A. Volpi and J. L. Bohn, Phys. Rev. A 65, 052712 (2002); R. V. Krems and A. Dalgarno, J. Chem. Phys. 120, 2296 (2004)]. Our approach is thus well suited for theoretical studies of strongly anisotropic molecular collisions in the presence of external electromagnetic fields.     

Molecular nitrogen-N2 properties: the intermolecular potential and the equation of state

Quantum mechanical (QM) high precision calculations were used to determine N(2)-N(2) intermolecular interaction potential. Using QM numerical data the anisotropic potential energy surface was obtained for all orientations of the pair of the nitrogen molecules in the rotation invariant form. The new N(2)-N(2) potential is in reasonably good agreement with the scaled potential obtained by van der Avoird et al. using the results of Hartree-Fock calculations [J. Chem. Phys. 84, 1629 (1986)]. The molecular dynamics (MD) of the N(2) molecules has been used to determine nitrogen equation of state. The classical motion of N(2) molecules was integrated in rigid rotor approximation, i.e., it accounted only translational and rotational degrees of freedom. Fincham [Mol. Simul. 11, 79 (1993)] algorithm was shown to be superior in terms of precision and energy stability to other algorithms, including Singer [Mol. Phys. 33, 1757 (1977)], fifth order predictor-corrector, or Runge-Kutta, and was therefore used in the MD modeling of the nitrogen pressure [S. Krukowski and P. Strak, J. Chem. Phys. 124, 134501 (2006)]. Nitrogen equation of state at pressures up to 30 GPa (300 kbars) and temperatures from the room temperature to 2000 K was obtained using MD simulation results. Results of MD simulations are in very good agreement (the error below 1%) with the experimental data on nitrogen equation of state at pressures below 1 GPa (10 kbars) for temperatures below 1800 K [R. T. Jacobsen et al., J. Phys. Chem. Ref. Data 15, 735 (1986)]. For higher temperatures, the deviation is slightly larger, about 2.5% which still is a very good agreement. The slightly larger difference may be attributed to the vibrational motion not accounted explicitly by rigid rotor approximation, which may be especially important at high temperatures. These results allow to obtain reliable equation of state of nitrogen for pressures up to 30 GPa (300 kbars), i.e., close to molecular nitrogen stability limit, determined by Nellis et al. [Phys. Rev. Lett. 53, 1661 (1984)].