Equation of state for mercury: revisited

An analytical equation of state by Song and Mason is developed to calculate the PVT properties of mercury. The equation of state is based on the statistical-mechanical perturbation theory of hard convex bodies and can be written as a fifth-order polynomial in the density. There exists three temperature-dependent parameters in the equation of state; the second virial coefficient, an effective molecular volume, and a scaling factor for the average contact pair distribution function of hard convex bodies. The temperature-dependant parameters have been calculated using corresponding-states correlations based on the surface tension and the liquid density at the normal boiling point. Employing the present equation of state, we have calculated the PVT properties of mercury over a wide range of temperatures and pressures. The average absolute deviation for the calculated density of mercury in the saturation and compressed states is 1.09%.     

基于化学缔合统计理论的链状流体状态方程

基于化学缔合统计理论的链状流体状态方程(EOS)能够反映实际分子的形状、链节成链、缔合等具体信息,在实际流体热力学性质计算中有着广泛应用.一般的链状流体EOS仅考虑相邻链节间的相关性,我们则借助统计力学和计算机模拟结果在模型中纳入了相间链节间的相关性,获得的硬球链流体(HSCF)模型能够更好地预测模型流体的压缩因子和第二维里系数.以HSCF为参考,引入方阱色散微扰项获得了实际方阱链流体(SWCF)EOS;结合根据黏滞球模型导得的缔合项,进一步构建了缔合流体EOS.最近,我们根据微扰理论和积分方程方法又开发了一新的变阱宽方阱链流体(SWCF-VR)模型.SWCF和SWCF-VREOSs可很好地用于计算小分子、聚合物、离子液体等纯流体及混合物的相行为、热焓、表面张力、黏度等热力学及传递性质,显示了模型良好的工程应用价值.本文就本课题组多年来在自由空间范畴内基于化学缔合统计理论开发链状流体EOS及其实际应用作系统的总结.     

Excited state coupling in low order many-body theory

By using a matrix analogue to the standard summation of diagonal ladder diagrams, it is shown that one can sum an infinite glass of diagrams and obtain excited state coupling in low order diagrams.     

Quantum variational transition state theory revisited

A quantum version of classical variational transition state theory suggested by McLafferty and Pechukas is refined. In this new quantum version, the variational property of the theory leads to the identification of an optimal smeared dividing surface. This optimal function is shown to be the eigenfunction associated with the lowest eigenvalue of a positive quantum transition state theory operator. The lowest eigenvalue is the optimal bound on the quantum rate. Application of the theory to the parabolic barrier provides better bounds but does not give an essential improvement when compared to previous quantum transition state theories.     

Self-consistent fluid variational theory for the dissociation of dense nitrogen

The self-consistent fluid variational theory is used to calculate the pressure dissociation of dense nitrogen at high temperatures. The accurate high-pressure and high-temperature effective pair potentials are adopted to describe the intermolecular interactions, which are made to consider molecular dissociation. This paper focuses on a mixture of nitrogen atoms and molecules and is devoted to the study of the phenomenon of pressure dissociation at finite temperature. The equation of state and dissociation degree are calculated from the free-energy functions in the range of temperature of 2000-15 000 K and density of 0.2-3.0 gcm(3), which can be compared with other approaches and experiments.     

Third-order many-body perturbation theory for the ground state of the carbon monoxide molecule

Many-body perturbation calculations have been performed for the ground state of the carbon monoxide molecule at its equilibrium internuclear separation. The calculations are complete through third order within the algebraic approximation; i.e., the state functions are parameterized by expansion in a finite basis set. All two-, three-, and four-body terms are rigorously determined, and many-body effects are found to be very important. A detailed comparison is made with a previously reported configuration interaction study. Padé approximants to the energy expansion are constructed. The many-body perturbative wave function is used in the Rayleigh quotient to produce upper bounds to the electronic energy.     

Isotropic hard core dumbell fluid: Equation of state

A monte-Carlo calculation has been made for 108 hard core dumbbells of the compressibility factor in the isotropic density range. These compressibility factors are compared to several approximate theories, and comments are made on the utility of these approximate theories.     

Statistical theory for hydrogen bonding fluid system of AaDd type (Ⅲ): Equation of state and fluctuations

The equation of the state of the hydrogen bonding fluid system of AaDd type is studied by the principle of statistical mechanics. The influences of hydrogen bonds on the equation of state of the system are obtained based on the change in volume due to hydrogen bonds. Moreover, the number density fluctuations of both molecules and hydrogen bonds as well as their spatial correlation property are investigated. Furthermore, an equation describing relation between the number density correlation function of "molecules-hydrogen bonds" and that of molecules and hydrogen bonds is derived. As application,taking the van der Waals hydrogen bonding fluid as an example, we considered the effect of hydrogen bonds on its relevant statistical properties.     

Equation of state for hard-sphere fluid in restricted geometry

Many practical applications require the knowledge of the equation of state of fluids in restricted geometry. We study a hard-sphere fluid at equilibrium in a narrow cylindrical pore with hard walls for pore radii R<((square root 3)+2)/4 (in units of the hard sphere diameter). In this case each particle can interact only with its nearest neighbors, which makes possible the use of analytical methods to study the thermodynamics of the system. Using a transfer operator formalism and expanding in low- and high-pressure regions, we can obtain a simple analytical equation of state for almost all ranges of pressure. The results agree with Monte Carlo simulations. Additionally, it is shown that a convenient analytical representation can be chosen to accurately describe the equation of state within the error of the Monte Carlo simulation.     

Equation of state of a seven-dimensional hard-sphere fluid. Percus-Yevick theory and molecular-dynamics simulations

Following the work of Leutheusser [Physica A 127, 667 (1984)], the solution to the Percus-Yevick equation for a seven-dimensional hard-sphere fluid is explicitly found. This allows the derivation of the equation of state for the fluid taking both the virial and the compressibility routes. An analysis of the virial coefficients and the determination of the radius of convergence of the virial series are carried out. Molecular-dynamics simulations of the same system are also performed and a comparison between the simulation results for the compressibility factor and theoretical expressions for the same quantity is presented.     

Modern state of intermolecular interaction theory

The concept of potential surfaces and classification of various types of intermolecular forces are given. The possibilities and the criteria of applicability of modern methods for potential surface calculation at short, intermediate, and long distances are discussed. Special attention is paid to the methods for calculating interactions between large molecules.     

Comparison of high-order many-body perturbation theory and configuration interaction for H2O

Diagrammatic many-body perturbation theory, coupled with a recursive computational procedure, is employed to obtain the correlation energy of H₂O within a 39-STO basis set by evaluating all double-excitation diagrams through twelfth order without any approximations. This provides, for the first time, the complete double-excitation diagrams contributions to the correlation energy, which is ?0.28826 hartree, compared with a correlation energy of ?0.27402 hartree obtained from a configuration interaction calculation which includes all double excitations. The difference of 0.0142 hartree includes the “size consistency” correction to the all-double-excitations CI energy, due to the “pathological” unliked-diagram terms remaining in that result, but also involves certain fourth- and higher-order rearrangement diagrams. Contrary to conventional belief, the unshifted, or Møller-Plesset partitioning of the hamiltonian provides a much more rapid convergence of the perturbation series that does the shifted, or Epstein-Nesbet partitioning. In both cases. Padé approximants enhance the convergence of the series considerably. A simple variation-perturbation scheme based on the first-order MBPT wavefunction is sufficient to provide 97.5% of the all-doubles CI correlation energy.     

Calculating rate constants with updated Hessians using variational transition state theory with multidimensional tunneling

Variational transition state theory with multidimensional tunneling (VTST/MT) has been used for calculating the rate constants of reactions. The updated Hessians have been used to reduce the computational costs for both geometry optimization and trajectory following procedures. In this paper, updated Hessians are used to reduce the computational costs while calculating the rate constants applying VTST/MT. Although we found that directly applying the updated Hessians will not generate good vibrational frequencies along the minimum energy path (MEP), however, we can either re-compute the full Hessian matrices at fixed intervals or calculate the Block Hessians, which is constructed by numerical one-side difference for the Hessian elements in the "critical" region and Bofill updating scheme for the rest of the Hessian elements. Due to the numerical instability of the Bofill update method near the saddle point region, we have suggested a simple strategy in which we follow the MEP until certain percentage of the classical barrier height from the barrier top with full Hessians computed and then performing rate constant calculation with the extended MEP using Block Hessians. This strategy results a mean unsigned percentage deviation (MUPD) around 10% with full Hessians computed till the point with 80% classical barrier height for four studied reactions. This proposed strategy is attractive not only it can be implemented as an automatic procedure but also speeds up the VTST/MT calculation via embarrassingly parallelization to a personal computer cluster.     

Multi-component molecular orbital theory for electrons and nuclei including many-body effect with full configuration interaction treatment: isotope effects on hydrogen molecules

A multi-component molecular orbital (MC_MO) theory is developed for a combined quantum system of electrons and nuclei with the full configuration interaction (CI) scheme of Cartesian Gaussian-type functions. The technique of graphical unitary group approach (GUGA) is modified to obtain the CI matrix elements for many kinds of quantum particles efficiently. The optimum basis sets for quantum nuclei are proposed with the fully variational procedure. The average internuclear distances, dipole polarizabilities, and nuclear vibrational excitation energies for isotopic hydrogen molecules calculated with optimized basis sets are found to adequately reproduce the corresponding experimental values.     

Transition state theory: variational formulation, dynamical corrections, and error estimates

Transition state theory (TST) is revisited, as well as evolutions upon TST such as variational TST in which the TST dividing surface is optimized so as to minimize the rate of recrossing through this surface and methods which aim at computing dynamical corrections to the TST transition rate constant. The theory is discussed from an original viewpoint. It is shown how to compute exactly the mean frequency of transition between two predefined sets which either partition phase space (as in TST) or are taken to be well-separated metastable sets corresponding to long-lived conformation states (as necessary to obtain the actual transition rate constants between these states). Exact and approximate criterions for the optimal TST dividing surface with minimum recrossing rate are derived. Some issues about the definition and meaning of the free energy in the context of TST are also discussed. Finally precise error estimates for the numerical procedure to evaluate the transmission coefficient kappaS of the TST dividing surface are given, and it is shown that the relative error on kappaS scales as 1/square root(kappaS) when kappaS is small. This implies that dynamical corrections to the TST rate constant can be computed efficiently if and only if the TST dividing surface has a transmission coefficient kappaS which is not too small. In particular, the TST dividing surface must be optimized upon (for otherwise kappaS is generally very small), but this may not be sufficient to make the procedure numerically efficient (because the optimal dividing surface has maximum kappaS, but this coefficient may still be very small).     

The accuracy of transition state theory in its absolute rate theory and variational formulations

The accuracy of the absolute rate theory and variational formulations of transition state theory have been determined for a series of collinear reactions A + BC → AB + C on the Porter—Karplus potential energy surface by comparing theoretical values of classical mean reaction probabilities with classical mechanical trajectory results.     

A comparison of different many-body perturbation theory calculations of the ground state of SiS

A comparison of different many-body perturbation theory (MBPT ) calculations of the ground state rotational and vibrational constants of SiS is made. The calculations are performed up to the complete fourth-order MBPT level, and in all cases two basis sets are utilized. The results of the third-order and some incomplete fourth-order calculations are in good agreement, but the complete fourth-order is among the worst as compared with the experimental data. Analysis of the different contributions to the calculated correlation eneriges points towards the necessity of including even higher-order terms of the(MBPT ) expansion.     

Equation of state based on the weeks-chandler-andersen perturbation theory

An equation of state based on the Weeks-Chandler-Andersen (WCA) perturbation theory utilizing the simplified random-phase approximation (RPA) term is p of state is applied successfully to calculate the P-V-T relationship for a spherical molecule, argon.For nonspherical molecules, the effects of anisotropic interactions are treated empirically. The calculated P-V-T relationships and saturated properties for nonspherical and nonpolar molecules agree well with experimental data. The potential parameters for nonpolar substances are well correlated with the acentric factors.     

Equation of state and structural properties of the Weeks-Chandler-Andersen fluid

Molecular dynamics simulations have been carried out for the equation of state and percolation properties of the Weeks-Chandler-Andersen (WCA) system in its fluid phase as functions of density and temperature. The compressibility factor Z collapses well for the various isotherms, using an effective particle diameter for the WCA particle which is (in the usual WCA reduced units) sigma(e)=2(16)(1+T)(16), where T is the temperature. A corresponding "effective" packing fraction is zeta(e)=pisigma(e) (3)N6V, for N particles in volume V, which therefore scales out the effects of temperature. Using zeta(e) the simulation derived Z can be fitted to a simple analytic form which is similar to the Carnahan-Starling hard sphere equation of state and which is valid at all temperatures and densities where the WCA fluid is thermodynamically stable. The data, however, are not scalable onto the hard sphere equation of state for the complete packing fraction range. We explored the continuum percolation behavior of the WCA fluids. The percolation distance sigma(p) for the various states collapses well onto a single curve when plotted as sigma(p)sigma(e) against zeta(e). The ratio sigma(p)sigma(e) exhibits a monotonic decrease with increasing zeta(e) between the percolation line for permeable spheres and the glass transition limit, where sigma(p)sigma(e) approximately 1. The percolation packing fraction was calculated as a function of effective packing fraction and fitted to an empirical expression. The local coordination number at the percolation threshold showed a transition between the soft core and hard core limits from ca. 2:74 to 1:5, as previously demonstrated in the literature for true hard spheres. A number of simple analytic expressions that represent quite well the percolation characteristics of the WCA system are proposed.     

A variational transition state theory description of periselectivity effects on cycloadditions of ketenes with cyclopentadiene

The reaction of cycloaddition of ketene and cyclopentadiene presents experimentally a competing mechanism where the branching ratio between the Woodward?CHoffmann allowed [4+2] and forbidden [2+2] cycloaddition product is 4.56?at ?20?°C, but because the minimum energy path misses the [2+2] product altogether, it has been claimed to lie beyond the scope of transition state theory. In this paper, a variational transition state theory study on this reaction is presented. It is found that the minimum energy path affording the [4+2] product travels through a potential energy plateau very close to the minimum energy path that describes the interconversion between both cycloaddition products, allowing the transfer between both and the direct formation of the forbidden [2+2] product, in this way acting as a means to circumvent the Woodward?CHooffmann rules. Within the domain of the competitive canonical unified statistical theory, a value for the branching ratio in very good agreement with experiment is computed.