The intermolecular potential energy functions of krypton and xenon

The intermolecular potential energy functions for krypton and xenon have been determined using new semi-inversion techniques. These techniques, which have previously been applied to the data for argon, enable information about intermolecular forces to be obtained directly from second virial coefficient and gas viscosities measurements.     

An intermolecular pair potential energy function for methane

A numerical spherically averaged intermolecular pair potential energy function is presented which describes the interactions of methane molecules. The function incorporates accurate estimates of the long-range dispersion force coefficients, and closely reproduces shear viscosity and second virial coefficient data over a wide temperature range. Agreement with other experimental data which is generally satisfactory is also discussed. The potential is shown to be superior to a large number of analytical potentials previously proposed for this gas.     

A reliable single parameter interatomic potential for argon

A simple semiempirical approximation previously proposed for the isotropic intermolecular forces between two closed shell systems is tested in detail for the argon-argon interaction. The potential is based on the knowledge of the first-order coulomb interaction energy, a suitably damped three term long range asymptotic expansion of the second order coulomb energy, and a semiempirical representation of the exchange interaction energy which contains one adjustable parameter. The single adjustable parameter can be reliably determined by fitting the second virial coefficient for argon in the 130–773 K temperature range with the long range interaction coefficients being constrained within the theoretical bounds specified by Tang, Norbeck and Certain. The reliability of the potential is compared with that of several literature potentials by comparing the theoretical predictions obtained from the potentials with experimental results for the second virial coefficient, viscosity, thermal conductivity and thermal diffusion ratios for dilute argon gas, and with spectroscopic data for the dimer, and with SCF calculations of the Ar-Ar potential at small interatomic separations. Our best potential predicts these properties with a precision as good as or better than other recent potentials which generally contain more adjustable parameters and/or involve more input data. The results confirm earlier work that suggested that the scheme tested is capable of yielding reliable isotropic potentials for the interaction of closed shell systems for 0·3 ? R/R_m ? ∞ where R_m is the intermolecular distance at the van der Waals minimum. The scheme appears to offer a method for obtaining reliable potentials while avoiding problems associated with optimizing many parameters with respect to fitting experimental constraints.     

The effect of steepness of soft-core square-well potential model on some fluid properties

The effect of repulsive steepness of the soft-core square well (SCSW) potential model on the second virial coefficient, critical behaviour (two- phase region and the position of critical point), and coordination number are investigated. The soft-core thermodynamic perturbation theory (TPT) presented by Weeks-Chandler-Anderson (WCA) recently developed by Ben-Amotz and Stell (BAS) has been used for the reference system, and the Barker-Henderson TPT for the perturbed system. The Barker-Henderson macroscopic compressibility approximation has been used for all order perturbation terms in which the second-order one is improved by assuming that the molecules in every two neighbouring shells are correlated upon the original assumption. By using the hard-sphere isothermal compressibility consistency for the radial distribution function (RDF), an analytical closed expression has been derived for the Helmholtz free energy function contained effective hard-sphere diameter. The accuracy of the model has been examined for the hard-core system, and an appropriate range found for the attractive width of the potential well (R), then the effect of steepness parameter on the critical quantities, coordination number, and the inversion temperature of the second virial coefficient, has been investigated qualitatively. The predicted results are in good agreement with the computer simulation data for the critical constants, and coordination number at the limit of the hard-core square-well potential model at least qualitatively, and for the attractive range 1.55 ≤ R ≤ 1.7, quantitatively. It was found that the steepness of the potential model has a marginal effect on the critical behaviour, and also every thermodynamic quantity at low and medium temperatures for which the molecular penetration is negligible, but since the penetration at high temperatures is significant, the role of the steepness of potential on the inversion temperature of the second virial coefficient and coordination number is highlighted.     

Intermolecular potential energy functions from gaseous viscosity data. Choice of the well depth

Dilute gas viscosity data may be inverted directly to give the intermolecular potential energy function if the well depth is known. The consequences of using different values of the well depth are studied, and it is concluded that the correct value may be distinguished by using second virial coefficient data.     

Inter-DNA attraction mediated by divalent counterions

Can nonspecifically bound divalent counterions induce attraction between DNA strands? Here, we present experimental evidence demonstrating attraction between short DNA strands mediated by Mg2+ ions. Solution small angle x-ray scattering data collected as a function of DNA concentration enable model independent extraction of the second virial coefficient. As the [Mg2+] increases, this coefficient turns from positive to negative reflecting the transition from repulsive to attractive inter-DNA interaction. This surprising observation is corroborated by independent light scattering experiments. The dependence of the observed attraction on experimental parameters including DNA length provides valuable clues to its origin.     

Magnetic anisotropies of some aromatic molecules

A common method for the estimation of uncertainties introduced by surface and impurity effects into experimental measurements of virial coefficients is described. The sign and the amplitude of the second virial coefficient response to perturbation caused by adsorption of molecules on the internal surface of the vessel have been determined. It has been shown that the magnitude of the second virial coefficient distortion depends on such competing factors as adsorption-impurity perturbation parameter, mixture composition which has been corrected taking into account this perturbation, and the nature of the impurity expressed in terms of its second virial coefficient and of the solvent-impurity cross second virial coefficient. The character of the Lennard-Jones 12–6 potential parameters perturbation, caused by the adsorption-impurity effects, is determined using second virial coefficient data inversion technique. Numerical estimates are made for nitrogen, helium, argon, xenon, their binary mixtures, and also for krypton-sulphur hexafluoride gaseous mixtures.     

H-H and C-H coupling constants in neutral and protonated azines

New measurements are reported of the density dependent depolarization ratio for argon, krypton, xenon, methane and sulphur hexafluoride, and the results are analysed to provide values for the second and third depolarization virial coefficients. The relationships between the second depolarization virial coefficient, the zeroth moment of the two-body Rayleigh spectrum and the second Kerr virial coefficient are considered, and it is shown that they now provide consistent results for the collision-induced pair polarizability anisotropy. Former inconsistencies are attributed to insufficient allowance for the effects of three-body interactions. Calculations of the second and third depolarization virial coefficients based on the DID model and using the Maitland-Smith potential are in excellent agreement with the experimental results for argon, krypton and xenon.     

Atomic interactions in neon and helium

The first two quantum corrections to the second virial coefficients of the Smith-Thakkar potential are calculated. Parameters for neon and helium, gases in which quantum effects are important, are then determined by fitting to semiempirical dispersion coefficients and experimental second virial coefficients. Viscosity coefficients for both gases and vibrational energy level spacings for the neon dimer are calculated as independent tests of the potentials. Overall agreement with experiment is excellent for neon and moderate for helium. The previously determined parameters for argon are found to be only very slightly perturbed by the inclusion of quantum corrections in the calculated second virial coefficients.     

Simple pair potential model for real fluids

Simple pair potential accounting realistically for both the short-range repulsions and long-range attractions is proposed and its main properties are discussed. It is shown that for simple fluids (i) its attractive force constant, C₆, corresponds to that found experimentally and (ii) it gives perfect agreement of the second virial coefficient over the entire temperature ranges for which experimental data are known.     

基于分子动力学模拟流体输运性质的稳定性分析

利用线性响应理论对Ar流体输运参数进行了分子动力学模拟,结果发现:导热系数和黏度会随着自相关积分函数积分时间的增加而产生剧烈波动,而扩散系数却相对稳定. 针对积分稳定性这一问题,对导热系数和黏度中的热流密度和应力张量进行了分解分析,发现含分子间作用力项是影响稳定性的最大因素. 从牛顿力学出发对作用力项的影响机理进行了分析,指明减小这种影响的最主要方法是使在体系进行统计输运参数前达到稳定平衡状态,即最小的预平衡步数应该满足使体系达到该状态下熵最大或者能量最低,并尽量减小温度对体系的影响. 同时,还对模拟盒尺寸、统计步长等因素对积分稳定性的影响进行了分析,给出了保持稳定性的建议. 关键词： 分子动力学 输运性质 自相关函数 稳定性     

Ab initio calculation of the refractivity and hyperpolarizability second virial coefficients of neon gas

Second virial coefficients for the density dependence of a number of electric properties are calculated for neon gas. Employing an accurate CCSD(T) potential for the Ne² van der Waals dimer and interaction-induced electric dipole polarizabilities and hyperpolarizabilities obtained from CCSD response theory, we evaluated the dielectric, refractivity, Kerr and ESHG second virial coefficients using both a semiclassical and a quantum statistical approach. The results cover a wide range of temperatures and are expected to be more reliable than the available experimental and empirical data. Quantum effects are found to be important only for temperatures below 100 K. The frequency-dependence of the refractivity virial coefficient is found to be small, but not negligible. For frequencies in the visible region it accounts for a few percent of the final results. For the ESHG virial coefficient of neon, frequency dependence is found to be very important, accounting for 20–25% of the second virial coefficient at the typical frequencies employed in experiments.     

Ab initio pair potential energy curve for the argon atom pair and thermophysical properties for the dilute argon gas. II. Thermophysical properties for low-density argon

A recent argon–argon interatomic potential energy curve determined from quantum-mechanical ab initio calculations and described with an analytical representation [B. Jäger, R. Hellmann, E. Bich, and E. Vogel, Mol. Phys. 107, 2181 (2009); 108, 105 (2010)] was used to calculate the most important thermophysical properties of argon governed by two-body interactions. Second pressure, acoustic, and dielectric virial coefficients as well as viscosity and thermal conductivity in the limit of zero density were computed for natural argon from 83 to 10,000?K. The calculated values for the different thermophysical properties are compared with available experimental data and values computed for other argon–argon potentials. This extensive analysis shows that the proposed potential is superior to all previous ones and that the calculated thermophysical property values are accurate enough to be applied as standard values for the complete temperature range of the calculations.     

Editorial

Using coupled cluster singles and doubles linear response theory and the d-aug-cc-pVTZ basis set extended with a 3s3p2d1f1g set of midbond functions, the interaction induced electric dipole polarisability surface of the CO–Ar van der Waals complex is computed. Combining this surface with accurate intermolecular potential energy and electric dipole surfaces, the pressure and dielectric second virial coefficients of the complex are calculated by a classical statistical approach. Excellent agreement with experimental results (to within the experimental error bars) is obtained for the pressure second virial coefficient over a range of temperatures. No previous experimental or theoretical investigations have been carried out for the dielectric second virial coefficient, B _ε(T), which is estimated to be about 1.9 cm⁶ mol^??1 at room temperature. This value results from a balance of terms due to the interaction induced electric dipole polarisability (predominant at high temperatures) and orientational electric dipole contributions.     

Ab initio intermolecular potential energy surfaces for the Ar–NCCN van der Waals complexes

The intermolecular potential energy surface of complex pairing argon with cyanogen molecule (NCCN) was calculated using the coupled cluster with single and double and perturbative triple excitations (CCSD(T)) with aug-cc-pvdz basis set extended with a set of mid-bond (3s3p2d1f1g) functions. The interaction energies were calculated by the supermolecular approach with the full counterpoise correction for the basis set superposition error. The calculated potential energies were fitted to an analytical expression. The calculated Ar–NCCN potential energy surface shows a global minimum at 3.35 Å, the distance between argon and centre of mass of cyanogen, for the T-shaped geometry and two local minimum at distance of 5.54 Å for the linear geometry on one side of cyanogen. Finally, the interaction second virial coefficients were calculated using the fitted potential energy surface and were compared with those obtained by the parameters of the Beattie–Bridgeman equation of states of pure argon and cyanogens fluids, approximately.     

Accurate second Kerr virial coefficient of rare gases from the state-of-the-art ab initio potentials and (hyper)polarizabilities

The second Kerr virial coefficient of rare gases is studied in this work using the best ab initio potentials and (hyper)polarizabilities in the literature. The second Kerr virial coefficient of helium-4, helium-3, neon, argon, and krypton and its polarizability component of xenon are computed by the semi-classical method together with the Padé approximant over a wide temperature range. In addition, the uncertainty of second Kerr virial coefficient is estimated from the uncertainties of the ab initio interaction-induced properties. The experimental and theoretical data in the literature is compared with our calculated values to examine the quality of this work. It is shown that our computed values in the supplementary materials are as accurate as the literature data at medium and high temperatures and are more reliable at low temperatures.     

The influence of intermolecular interactions on the Kerr effect in gases

Expressions are derived for the second and third Kerr virial coefficients, B _K and C _K, of spherical top molecules in terms of irreducible cluster polarizabilities, and values are calculated using the dipole-induced dipole model for argon, krypton, xenon, methane, tetrafluoromethane, neopentane and sulphur hexafluoride. For mixtures of rare gases it is shown that the collision-induced dipole moment makes a negligible contribution to B _K. The effect of the choice of intermolecular potential function on the calculated second Kerr virial coefficients is also demonstrated. It is found that the predominant contributions to C _K arise from the pair polarizability, and that the triplet polarizability is only of minor importance.     

Molecular dynamics studies of thermal conductivity time correlation functions

The Green–Kubo time correlation function for the thermal conductivity in liquid argon is studied for a thermodynamic state close to the triple point by standard molecular dynamics simulations. The collective heat flux vector has been separated into contributions originated at the kinetic energy, the intermolecular potential and the pair virial function. Furthermore, the Green–Kubo time correlation functions have been broken down into partial n-body terms (n=1,2,3,4). The most important contribution to the thermal conductivity is represented by the auto correlation of the virial term. In contrast to other collective phenomena described by time correlation functions involving n-body terms, the partial Green–Kubo time correlation functions for the thermal conductivity are not affected by exponential long-time tails.