Asymptotic trends in thermodynamic perturbation theory

The development of transferable force fields for n-alkanes has enabled molecular-dynamics simulation of the reference (A0) and perturbation (A1,A2) terms in thermodynamic perturbation theory (TPT) over a broad range of chain length. The implied equations of state yield 9.1% average error in vapor pressure and 4.7% error in liquid density for compounds ranging from propane to triacontane. Further simulations extend to nC80, but there are no experimental data to which comparisons can be made. With reliable TPT terms from molecular simulation, it is possible to analyze the trends with respect to molecular weight. Each TPT contribution is shown to approach an asymptote in the long chain limit. The asymptotes and the approaches to them are quantitatively characterized. A0 and A1 approach their asymptotes at relatively short chain lengths (nC30). A2, on the other hand, approaches its asymptote slowly (nC80). Simulation-based TPT terms also permit unambiguous interpretation of the number of coarse-grained segments relative to the number of carbons in the chain. Previous attempts have relied on characterizations that included the repulsive and attractive contributions simultaneously in a manner susceptible to a cancellation of errors. In this work, the reference fluid alone provides the characterization and the result is shown to be consistent with expectations for the A1 term. The conclusion is that the number of carbons per segment approaches roughly 10 in the long chain limit, much larger than previously reported. A small adjustment to the chain contribution from Wertheim's [J. Stat. Phys. 42, 477 (1986)] TPT1 model is sufficient to provide quantitative accuracy for A0.     

An improved thermodynamic perturbation theory for Mercedes-Benz water

We previously applied Wertheim's thermodynamic perturbation theory for associative fluids to the simple Mercedes-Benz model of water. We found that the theory reproduced well the physical properties of hot water, but was less successful in capturing the more structured hydrogen bonding that occurs in cold water. Here, we propose an improved version of the thermodynamic perturbation theory in which the effective density of the reference system is calculated self-consistently. The new theory is a significant improvement, giving good agreement with Monte Carlo simulations of the model, and predicting key anomalies of cold water, such as minima in the molar volume and large heat capacity, in addition to giving good agreement with the isothermal compressibility and thermal expansion coefficient.     

Hard sphere perturbation theory for fluids with soft-repulsive-core potentials

The thermodynamic properties of fluids with very soft repulsive-core potentials, resembling those of some liquid metals, are predicted with unprecedented accuracy using a new first-order thermodynamic perturbation theory. This theory is an extension of Mansoori-Canfield/Rasaiah-Stell (MCRS) perturbation theory, obtained by including a configuration integral correction recently identified by Mon, who evaluated it by computer simulation. In this work we derive an analytic expression for Mon's correction in terms of the radial distribution function of the soft-core fluid, g(0)(r), approximated using Lado's self-consistent extension of Weeks-Chandler-Andersen (WCA) theory. Comparisons with WCA and MCRS predictions show that our new extended-MCRS theory outperforms other first-order theories when applied to fluids with very soft inverse-power potentials (n< or =6), and predicts free energies that are within 0.3 kT of simulation results up to the fluid freezing point.     

Calculating thermodynamic properties from perturbation theory: I. An analytic representation of square-well potential hard-sphere perturbation theory

The Barker–Henderson macroscopic compressibility approximation of the second-order perturbation term is improved by assuming that the numbers of molecules in every two neighbour shells are correlated, based upon the original assumptions. The results are better than those for the original macroscopic compressibility and local compressibility approximation, especially at high densities. A simple analytic representation of square-well potential hard-sphere perturbation theory is derived based upon this improvement. The method is tested by calculating thermodynamic properties with the four-term truncated form, and the results are in good agreement with those of Monte Carlo and Molecular Dynamics simulation.     

Thermodynamic perturbation theory for fused sphere hard chain fluids using nonadditive interactions

A model is developed for the equation of state of fused chains based on Wertheim thermodynamic perturbation theory and nonadditive size interactions. The model also assumes that the structure (represented by the radial distribution function) of the fused chain fluid is the same as that of the touching hard sphere chain fluid. The model is completely based on spherical additive and nonadditive size interactions. The model has the advantage of offering good agreement with simulation data while at the same time being independent of fitted parameters. The model is most accurate for short chains, small values of Delta (slightly fused spheres) and at intermediate (liquidlike) densities.     

Calculation of free energy surfaces using the methods of thermodynamic perturbation theory

A molecular dynamics method for determining the free energy difference between systems separated in configuration space has been developed. With this new approach, which is based on thermodynamic perturbation techniques, potentials of mean force for conformational changes may be calculated. As a test of the method, the potential of mean force and radial distribution function for liquid argon have been computed. The results are in good agreement with those obtained from an ordinary simulation.     

Resummed thermodynamic perturbation theory for central force associating potential. Multi-patch models

A resummed thermodynamic perturbation theory for associating fluids with multiply bondable central force associating potential is extended for the fluid with multiple number of multiply bondable associating sites. We consider a multi-patch hard-sphere model for associating fluids. The model is represented by the hard-sphere fluid system with several spherical attractive patches on the surface of each hard sphere. Resummation is carried out to account for blocking effects, i.e., when the bonding of a particle restricts (blocks) its ability to bond with other particles. Closed form analytical expressions for thermodynamical properties (Helmholtz free energy, pressure, internal energy, and chemical potential) of the models with arbitrary number of doubly bondable patches at all degrees of the blockage are presented. In the limiting case of total blockage, when the patches become only singly bondable, our theory reduces to Wertheim's thermodynamic perturbation theory (TPT) for polymerizing fluids. To validate the accuracy of the theory we compare to exact values, for the thermodynamical properties of the system, as determined by Monte Carlo computer simulations. In addition we compare the fraction of multiply bonded particles at different values of the density and temperature. In general, predictions of the present theory are in good agreement with values for the model calculated using Monte Carlo simulations, i.e., the accuracy of our theory in the case of the models with multiply bondable sites is similar to that of Wertheim's TPT in the case of the models with singly bondable sites.     

Monte carlo simulation and thermodynamic perturbation theory for mixtures of diatomic molecules

The thermodynamic properties and site—site distribution functions of mixtures of non-spherical molecules are obtained by Monte Carlo simulation. A non-spherical reference-system perturbation theory based on the RISM equation is developed to predict these results. The agreement between theory and simulation for the thermodynamic properties is encouraging. Important differences in the relative peak heights of the site—site distribution functions from theory and simulation are attributed to the role of attractive forces in determining local structure in the fluid mixtures, where the volumes of the components are similar but the well depths differ.     

Mixing rules for van-der-Waals type equation of state based on thermodynamic perturbation theory

A concept based on the thermodynamic perturbation theory for a `simple fluid' has been applied to the attractive term of a van-der-Waals type equation of state (EOS) to derive a simple mixing rule for the a parameter. The new mixing rule is a small correction to the original one-fluid approximation to account for the influence of particles of j-type on the correlation function of ii-type in a mixture consisting of particles of i and j types. The importance of the correction has been shown by comparison of the calculated results for binary mixtures of Lennard–Jones fluids with the data obtained by numerical method (Monte-Carlo simulation). The new mixing rules can be considered as a flexible generalization of the conventional mixing rules and can be reduced to the original v-d-W mixing rules by defaulting the extra binary parameters to zero. In this way the binary parameters already available in the literature for many systems can be used without any additional regression work. Extension of the new mixing rules to a multicomponent system do not suffer from `Michelsen–Kistenmacher syndrome' and provide the correct limit for the composition dependence of second virial coefficients. Their applicability has been illustrated by various examples of vapor–liquid and liquid–liquid equilibria using a modified Patel–Teja EOS. The new mixing rules can be applied to any EOS of van-der-Waals type, i.e., EOS containing two terms which reflect the contributions of repulsive and attractive intermolecular forces.     

Fluid-fluid and fluid-solid transitions in the Kern-Frenkel model from Barker-Henderson thermodynamic perturbation theory

We study the Kern-Frenkel model for patchy colloids using Barker-Henderson second-order thermodynamic perturbation theory. The model describes a fluid where hard sphere particles are decorated with one patch, so that they interact via a square-well potential if they are sufficiently close one another, and if patches on each particle are properly aligned. Both the gas-liquid and fluid-solid phase coexistences are computed and contrasted against corresponding Monte Carlo simulations results. We find that the perturbation theory describes rather accurately numerical simulations all the way from a fully covered square-well potential down to the Janus limit (half coverage). In the region where numerical data are not available (from Janus to hard-spheres), the method provides estimates of the location of the critical lines that could serve as a guideline for further efficient numerical work at these low coverages. A comparison with other techniques, such as integral equation theory, highlights the important aspect of this methodology in the present context.     

Solid phase thermodynamic perturbation theory: test and application to multiple solid phases

A simple procedure for the determination of hard sphere (HS) solid phase radial distribution function (rdf) is proposed, which, thanks to its physical foundation, allows for extension to other crystal structures besides the fcc structure. The validity of the procedure is confirmed by comparing (1) the predicted HS solid phase rdf's with corresponding simulation data and (2) the predicted non-HS solid phase Helmholtz free energy by the present solid phase first-order thermodynamic perturbation theory (TPT) whose numerical implementation depends on the HS solid phase rdf's as input, with the corresponding predictions also by the first-order TPT but the required HS solid phase rdf is given by an "exact" empirical simulation-fitted formula. The present solid phase first-order TPT predicts isostructural fcc-fcc transition of a hard core attractive Yukawa fluid, in very satisfactory agreement with the corresponding simulation data and is far more accurate than a recent thermodynamically consistent density functional perturbation theory. The present solid phase first-order TPT is employed to investigate multiple solid phases. It is found that a short-ranged potential, even if it is continuous and differentiable or is superimposed over a long-ranged potential, is sufficient to induce the multiple solid phases. When the potential range is short enough, not only isostructural fcc-fcc transition but also isostructural bcc-bcc transition, simple cubic (sc)-sc transition, or even fcc-bcc, fcc-sc, and bcc-sc transitions can be induced. Even triple point involving three solid phases becomes possible. The multiple solid phases can be stable or metastable depending on the potential parameters.     

Performance evaluation of third-order thermodynamic perturbation theory and comparison with existing liquid state theories

To evaluate the performance of a recently proposed third-order thermodynamic perturbation theory (TPT), we employ the third TPT for calculation of thermodynamic properties such as compressibility factor, internal energy, excess chemical potential, gas-liquid coexistence curve, and critical properties of several fluids. By comparing the third-order TPT results with corresponding simulation data available in literature and supplied in the present report and theoretical results from several other theoretical approaches, one concludes that the third-order TPT is, in general, more accurate than other approaches such as Barker-Henderson second-order TPT using a macroscopic compressibility approximation (MCA-TPT), self-consistent Ornstein-Zernike approach, Monte Carlo perturbation theory, and a specially devised equation of state. Specifically, the third-order TPT can predict quantitatively a double critical phenomena of gas-liquid transition and a low-density liquid (LDL)-high-density liquid (HDL) transition associated with a soft core (SC) potential fluid very satisfactorily, but the predictions for the LDL-HDL transition based on the second-order MCA-TPT are quantitatively very bad or qualitatively incorrect. The failure of the second-order MCA-TPT for the SC fluid can be ascribed to the facts that for the SC potential the second-order and third-order terms of the perturbation expansion are not small quantities and that the second-order term is underestimated by the MCA. It is concluded that the present third-order version of the TPT is reliable for varying model fluids.     

Outer sphere perturbation of delocalized mixed-valence complexes

The complexes [[Ru(ttp)(bpy)](2)(micro-adpc)][PF(6)](2) and [[Ru(ttp)(bpy)](2)(micro-dicyd)][PF(6)](2), where ttp is 4-toluene-2,2':6',2' '-terpyridine, bpy is 2,2'-bipyridine, adpc(2)(-) is azodi(phenylcyanamide), and dicyd(2)(-) is 1,4-dicyanamidebenzene, were prepared and characterized by IR and NIR, vis spectroelectrochemistry, and cyclic voltammetry. The crystal structure of the complex, [[Ru(ttp)(bpy)](2)(micro-adpc)][PF(6)](2).6DMF, revealed a planar bridging adpc(2)(-) ligand with the cyanamide groups adopting an anti configuration. IR and comproportionation data are consistent with delocalized mixed-valence complexes, and a spectroscopic analysis assuming C(2)(h) microsymmetry leads to a prediction of multiple MMCT transitions with the lowest energy transition equal to the resonance exchange integral for the mixing of ruthenium donor and acceptor orbitals with a bridging ligand orbital (the preferred superexchange pathway). The solvent dependence of the MMCT band energy that is seen for [[Ru(ttp)(bpy)](2)(micro-adpc)](3+) is due to a ground state weakening of metal-metal coupling because of solvent donor interactions with the acceptor azo group of the bridging ligand.     

Stationary perturbation theory

Summary After a short recapitulation of the basic concepts of stationary perturbation theory, this is applied to a many-electron Hamiltonian, with or without an external field, given in a Fock space formulation in terms of a finite basis, the exact eigenfunctions of which are the full-CI wave functions. The Lie algebra _c ⁿ of the variational group corresponding to this problem is presented. It has an important subalgebra _c ⁽¹⁾ of one-particle transformations. Hartree-Fock and coupled Hartree-Fock (also uncoupled Hartree-Fock) as well as MC-SCF and coupled MC-SCF are outlined in this framework. Many-body perturbation theory and Møller-Plesset perturbation theory are derived from the same kind of stationarity condition and a new non-perturbative iteration construction of the full-CI wave function is proposed, the first Newton-Raphson iteration cycle of which is CEPA-0. For the treatment of electron correlation for properties two variants of Møller-Plesset theory referred to as coupled (CMP) and uncoupled (UCMP) are defined, neither of which is fully satisfactory. While CMP satisfies a Brillouin condition, which implies that first order correlation corrections to first- and second-order properties vanish, it does not satisfy a Hellmann-Feynman theorem, i.e. a first order property isnot the expectation value of the operator associated with the property. Conversely UCMP satisfies a Hellmann-Feynman theorem but no Brillouin theorem. The incompatibility of the two theorems is related to an unbalanced treatment of one-particle- and higher excitations in MP theory. CMP, which is based on coupled Hartree-Fock as uncorrelated reference, appears to have slight advantages over UCMP, but neither variant looks very promising for the evaluation of 2nd order correlation corrections to 2nd-order properties. Then four variants of the perturbation theory of properties with a nonperturbative treatment of electron correlation on CEPA-0 level (but extendable to a higher level) are discussed. While those variants which are the direct counterpart of UCMP and CMP must be discarded, the perturbative CEPA-0 derived from a perturbative treatment on full-CI level appears to satisfy all important criteria, in particular it satisfies a Brillouin-Brueckner condition and a Hellmann-Feynman theorem. A simplified version, the coupled Brillouin-Brueckner CEPA-0 appears to have essentially the same qualities. It is important to replace the Brillouin condition of MP theory by the Brillouin-Brueckner condition in non-perturbative approaches, especially if one is interested in properties.     

Large-Order perturbation theory

The original motivation for studying the asymptotic behavior of the coefficients of perturbation series came from quantum field theory. An overview is given of some of the attempts to understand quantum field theory beyond finite-order perturbation series. At least in the case of the Thirring model and probably in general, the full content of a relativistic quantum field theory cannot be recovered from its perturbation series. This difficulty, however, does not occur in quantum mechanics, and the anharmonic oscillator is used to illustrate the methods used in large-order perturbation theory. Two completely different methods are discussed, the first one using the WKB approximation, and a second one involving the statistical analysis of Feynman diagrams. The first one is well developed and gives detailed information about the desired asymptotic behavior, while the second one is still in its infancy and gives instead information about the distribution of vertices of the Feynman diagrams.     

Calculating thermodynamic properties from perturbation theory: Part III. The mixtures of square-well chain fluid

A completely analytic perturbation theory has been developed to calculate the Helmholtz energy, compressibility factor, internal energy and constant-volume heat capacity for square-well chain fluid mixtures. This theory is based on the improved Barker–Henderson macroscopic compressibility (mc) approximation proposed by Zhang, the first-order perturbation theory of Wertheim in which Zhang’s analytic monomer radial distribution function as the function of temperature and monomer density is used, and a simple mixing rule similar to that of Hino–Prausnitz. The validity of the perturbation theory is evaluated by comparing the calculated compressibility factor, internal energy and constant-volume heat capacity for the freely jointed square-well chain mixtures from the theory to MC simulation data. The results show that the theory predicts results in good agreement with simulation results.     

On perturbation graph theory

A perturbation theoretical formalism is developed which enables the calculation of the (topological) resonance energy of arbitrary heteroconjugated π-electron systems. The previous method of Herndon and Párkányi is thus generalized.     

Auxiliary density perturbation theory

A new approach, named auxiliary density perturbation theory, for the calculation of second energy derivatives is presented. It is based on auxiliary density functional theory in which the Coulomb and exchange-correlation potentials are expressed by auxiliary function densities. Different to conventional coupled perturbed Kohn-Sham equations the perturbed density matrix is obtained noniteratively by solving an inhomogeneous equation system with the dimension of the auxiliary function set used to expand the auxiliary function density. A prototype implementation for the analytic calculation of molecular polarizabilities is presented. It is shown that the polarizabilities obtained with the newly developed auxiliary density perturbation approach match quantitative with the ones from standard density functional theory if augmented auxiliary function sets are used. The computational advantages of auxiliary density perturbation theory are discussed, too.