Search

The phenomenon of underscreening in concentrated electrolyte solutions leads to a larger decay length of the charge-charge correlation than the prediction of Debye-Hückel (DH) theory and has found a resurgence of both theoretical and experimental interest in the chemical physics community. To systematically understand and investigate this phenomenon in electrolytes requires a theory of concentrated electrolytes to describe charge-charge correlations beyond the DH theory. We review the theories of electrolytes that can transition from the DH limit to concentrations where charge correlations dominate, giving rise to underscreening and the associated Kirkwood Transitions (KTs). In this perspective, we provide a conceptual approach to a theoretical formulation of electrolyte solutions that exploits the competition between molecular-informed short-range (SR) and long-range interactions. We demonstrate that all deviations from the DH limit for real electrolyte solutions can be expressed through a single function ΣQ that can be determined both theoretically and numerically. Importantly, ΣQ can be directly related to the details of SR interactions and, therefore, can be used as a tool to understand how differences in representations of interaction can influence collective effects. The precise function form of ΣQ can be inferred through a Gaussian field theory of both the number and charge densities. The resulting formulation is validated by experiment and can accurately describe the collective phenomenon of screening in concentrated bulk electrolytes. Importantly, the Gaussian field theory predictions of the screening lengths appear to be less than ∼1 nm at concentrations above KTs., (© 2024 Author(s). All article content, except where otherwise noted, is licensed under a Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).)

We establish the connection between the measured small angle x-ray scattering signal and the charge-charge correlations underlying Kirkwood transitions (KTs) in 1:1, 2:1, and 3:1 aqueous electrolytes. These measurements allow us to obtain underscreening lengths for bulk electrolytes independently verified by theory and simulations. Furthermore, we generalize the concept of KTs beyond those theoretically predicted for 1:1 electrolytes, which involves the inverse screening length, a0, and the inverse periodicity length, Q0. Above the KTs, we find a universal scaling of a0∝c-ζ/3 and Q0 ∝ c1/3 for the studied electrolyte solutions, where ζ is the ionic strength factor., (© 2024 Author(s). All article content, except where otherwise noted, is licensed under a Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).)

We study the prototypical S N 2 reaction Cl - + CH 3 Cl → CH 3 Cl + Cl - in water using quantum mechanics/molecular mechanics (QM/MM) computer simulations with transition path sampling and inertial likelihood maximization. We have identified a new solvent coordinate to complement the original atom-exchange coordinate used in the classic analysis by Chandrasekhar, Smith, and Jorgensen [J. Am. Chem. Soc. 107, 154 (1985)]. The new solvent coordinate quantifies instantaneous solvent-induced polarization relative to the equilibrium average charge density at each point along the reaction pathway. On the basis of likelihood scores and committor distributions, the new solvent coordinate improves upon the description of solvent dynamical effects relative to previously proposed solvent coordinates. However, it does not increase the transmission coefficient or the accuracy of a transition state theory rate calculation.

CP2K is an open source electronic structure and molecular dynamics software package to perform atomistic simulations of solid-state, liquid, molecular, and biological systems. It is especially aimed at massively parallel and linear-scaling electronic structure methods and state-of-the-art ab initio molecular dynamics simulations. Excellent performance for electronic structure calculations is achieved using novel algorithms implemented for modern high-performance computing systems. This review revisits the main capabilities of CP2K to perform efficient and accurate electronic structure simulations. The emphasis is put on density functional theory and multiple post-Hartree-Fock methods using the Gaussian and plane wave approach and its augmented all-electron extension.

Predicting accurate nuclear magnetic resonance chemical shieldings relies upon cancellation of different types of errors between the theoretically calculated shielding constant of the analyte of interest and the reference. Often, the intrinsic error in computed shieldings due to basis sets, approximations in the Hamiltonian, description of the wave function, and dynamic effects is nearly identical between the analyte and reference, yet if the electronic structure or sensitivity to local environment differs dramatically, this cannot be taken for granted. Detailed prior work has examined the octahedral trivalent cation Al(H 2 O) 6 3+ , accounting for ab initio intrinsic errors. However, the use of this species as a reference for the chemically distinct tetrahedral anion Al(OH) 4 - requires an understanding of how these errors cancel in order to define the limits of accurately predicting Al27 chemical shielding in Al(OH) 4 - . In this work, we estimate the absolute shielding of the Al27 nucleus in Al(OH) 4 - at the coupled cluster level (515.1 ± 5.3 ppm). Shielding sensitivity to the choice of method approximation and atomic basis sets used has been evaluated. Solvent and thermal effects are assessed through ensemble averaging techniques using ab initio molecular dynamics. The contribution of each type of intrinsic error is assessed for the Al(H 2 O) 6 3+ and Al(OH) 4 - ions, revealing significant differences that fundamentally hamper the ability to accurately calculate the Al27 chemical shift of Al(OH) 4 - from first principles.

Confinement of water in sub-nanometer pores strongly alters its vibrational dynamics from that of bulk water. The effect of confinement can, furthermore, be finely tuned by small changes in the size and symmetry of the confining pore. Using inelastic neutron scattering (INS), we recently studied the dynamics of water confined in the channels of beryl and cordierite in which, at low temperatures, water shows similar behavior, indicating an absence of hydrogen bonds acting on the water molecule and a shallow water potential in the direction perpendicular to the channels. In addition, we observed multiple tunneling modes (between 0.66 and 14.7 meV) in the INS spectra of beryl due to transitions between the split ground-state of the water protons. Here, we present a study of (i) the effect of pressure on the dynamics of water in beryl, (ii) the dynamics of water in beryl containing alkali metals (which results in changing the orientation of the water molecule in the crystal), and (iii) the dynamics of water in cordierite at low energies. We found a shift in the tunneling and vibrational modes of water in beryl to higher energies at 22 kbar relative to 1 bar. No tunneling modes were observed for water in cordierite and type-II water in beryl. Therefore, we conclude that very small differences in the size and structure of the pores and the orientation of the water molecule in these minerals result in changes in the potential of the water protons and drastic changes in the confined water dynamics.

The aqueous solvation structure of the Na + ion is studied both by molecular dynamics (MD) simulations and K-edge X-ray absorption near edge spectroscopy (XANES). We present a systematic study contrasting the differences in the predicted XANES spectra, using molecular configurations generated from classical and quantum MD. Using a set of suitable order parameters, we elucidate how the spectroscopic features are influenced by the local solvation structure around the Na + ion. Our findings suggest that XANES is sensitive to fluctuations between the first and second solvation shells and can distinguish between classical and quantum interaction representations used to generate molecular ensembles.

Studies on ion pairing at interfaces have been intensified recently because of their importance in many chemical reactive phenomena, such as ion-ion interactions that are affected by interfaces and their influence on kinetic processes. In this study, we performed simulations to examine the thermodynamics and kinetics of small polarizable sodium iodide ions in the bulk and near the water liquid-vapor interface. Using classical transition state theory, we calculated the dissociation rates and corrected them with transmission coefficients obtained from the reactive flux formalism and Grote-Hynes theory. Our results show that in addition to affecting the free energy of ions in solution, the interfacial environments significantly influence the kinetics of ion pairing. The results on the relaxation time obtained using the reactive flux formalism and Grote-Hynes theory present an unequivocal picture that the interface suppresses ion dissociation. The effects of the use of molecular models on the ion interactions as well as the ion-pair configurations at the interface are also quantified and discussed.

Determining the solvation free energies of single ions in water is one of the most fundamental problems in physical chemistry and yet many unresolved questions remain. In particular, the ability to decompose the solvation free energy into simple and intuitive contributions will have important implications for models of electrolyte solution. Here, we provide definitions of the various types of single ion solvation free energies based on different simulation protocols. We calculate solvation free energies of charged hard spheres using density functional theory interaction potentials with molecular dynamics simulation and isolate the effects of charge and cavitation, comparing to the Born (linear response) model. We show that using uncorrected Ewald summation leads to unphysical values for the single ion solvation free energy and that charging free energies for cations are approximately linear as a function of charge but that there is a small non-linearity for small anions. The charge hydration asymmetry for hard spheres, determined with quantum mechanics, is much larger than for the analogous real ions. This suggests that real ions, particularly anions, are significantly more complex than simple charged hard spheres, a commonly employed representation.

First principles molecular dynamics simulation protocol is established using revised functional of Perdew-Burke-Ernzerhof (revPBE) in conjunction with Grimme's third generation of dispersion (D3) correction to describe the properties of water at ambient conditions. This study also demonstrates the consistency of the structure of water across both isobaric (NpT) and isothermal (NVT) ensembles. Going beyond the standard structural benchmarks for liquid water, we compute properties that are connected to both local structure and mass density fluctuations that are related to concepts of solvation and hydrophobicity. We directly compare our revPBE results to the Becke-Lee-Yang-Parr (BLYP) plus Grimme dispersion corrections (D2) and both the empirical fixed charged model (SPC/E) and many body interaction potential model (MB-pol) to further our understanding of how the computed properties herein depend on the form of the interaction potential.

A combination of theory, X-ray diffraction (XRD) and extended x-ray absorption fine structure (EXAFS) are used to probe the hydration structure of aqueous Na + . The high spatial resolution of the XRD measurements corresponds to Q max = 24 Å -1 while the first-reported Na K-edge EXAFS measurements have a spatial resolution corresponding to 2k = Q max = 16 Å -1 . Both provide an accurate measure of the shape and position of the first peak in the Na-O pair distribution function, g NaO (r). The measured Na-O distances of 2.384 ± 0.003 Å (XRD) and 2.37 ± 0.024 Å (EXAFS) are in excellent agreement. These measurements show a much shorter Na-O distance than generally reported in the experimental literature (Na-O avg ∼ 2.44 Å) although the current measurements are in agreement with recent neutron diffraction measurements. The measured Na-O coordination number from XRD is 5.5 ± 0.3. The measured structure is compared with both classical and first-principles density functional theory (DFT) simulations. Both of the DFT-based methods, revPBE and BLYP, predict a Na-O distance that is too long by about 0.05 Å with respect to the experimental data (EXAFS and XRD). The inclusion of dispersion interactions (-D3 and -D2) significantly worsens the agreement with experiment by further increasing the Na-O distance by 0.07 Å. In contrast, the use of a classical Na-O Lennard-Jones potential with SPC/E water accurately predicts the Na-O distance as 2.39 Å although the Na-O peak is over-structured with respect to experiment.

X-ray diffraction measurements of liquid water are reported at pressures up to 360 MPa corresponding to a density of 0.0373 molecules per Å(3). The measurements were conducted at a spatial resolution corresponding to Q(max) = 16 Å(-1). The method of data analysis and measurement in this study follows the earlier benchmark results reported for water under ambient conditions having a density of 0.0333 molecules per Å(3) and Q(max) = 20 Å(-1) [J. Chem. Phys. 138, 074506 (2013)] and at 70 °C having a density of 0.0327 molecules per Å(3) and Q(max) = 20 Å(-1) [J. Chem. Phys. 141, 214507 (2014)]. The structure of water is very different at these three different T and P state points and thus they provide the basis for evaluating the fidelity of molecular simulation. Measurements show that at 360 MPa, the 4 waters residing in the region between 2.3 and 3 Å are nearly unchanged: the peak position, shape, and coordination number are nearly identical to their values under ambient conditions. However, in the region above 3 Å, large structural changes occur with the collapse of the well-defined 2nd shell and shifting of higher shells to shorter distances. The measured structure is compared to simulated structure using intermolecular potentials described by both first-principles methods (revPBE-D3) and classical potentials (TIP4P/2005, MB-pol, and mW). The DFT-based, revPBE-D3, method and the many-body empirical potential model, MB-pol, provide the best overall representation of the ambient, high-temperature, and high-pressure data. The revPBE-D3, MB-pol, and the TIP4P/2005 models capture the densification mechanism, whereby the non-bonded 5th nearest neighbor molecule, which partially encroaches the 1st shell at ambient pressure, is pushed further into the local tetrahedral arrangement at higher pressures by the more distant molecules filling the void space in the network between the 1st and 2nd shells.

Thermal fluctuation and hydrophobicity are two hallmarks of fluid hydrodynamics on the nano-scale. It is a challenge to consistently couple the small length and time scale phenomena associated with molecular interaction with larger scale phenomena. The development of this consistency is the essence of mesoscale science. In this study, we use a nanoscale fluid model based on smoothed dissipative particle dynamics that accounts for the phenomena associated with density fluctuations and hydrophobicity. We show consistency in the fluctuation spectrum across scales. In doing so, it is necessary to account for finite fluid particle size. Furthermore, we demonstrate that the present model can capture the void probability and solvation free energy of nonpolar hard particles of different sizes. The present fluid model is well suited for an understanding of emergent phenomena in nano-scale fluid systems.

Results of ab initio molecular dynamics (AIMD) simulations (density functional theory+PBE96) of the dynamics of waters in the hydration shells surrounding the Zn(2+) ion (T approximately 300 K, rho approximately 1 gm/cm(3)) are compared to simulations using a combined quantum and classical molecular dynamics [AIMD/molecular mechanical (MM)] approach. Both classes of simulations were performed with 64 solvating water molecules ( approximately 15 ps) and used the same methods in the electronic structure calculation (plane-wave basis set, time steps, effective mass, etc.). In the AIMD/MM calculation, only six waters of hydration were included in the quantum mechanical (QM) region. The remaining 58 waters were treated with a published flexible water-water interaction potential. No reparametrization of the water-water potential was attempted. Additional AIMD/MM simulations were performed with 256 water molecules. The hydration structures predicted from the AIMD and AIMD/MM simulations are found to agree in detail with each other and with the structural results from x-ray data despite the very limited QM region in the AIMD/MM simulation. To further evaluate the agreement of these parameter-free simulations, predicted extended x-ray absorption fine structure (EXAFS) spectra were compared directly to the recently obtained EXAFS data and they agree in remarkable detail with the experimental observations. The first hydration shell contains six water molecules in a highly symmetric octahedral structure is (maximally located at 2.13-2.15 A versus 2.072 A EXAFS experiment). The widths of the peak of the simulated EXAFS spectra agree well with the data (8.4 A(2) versus 8.9 A(2) in experiment). Analysis of the H-bond structure of the hydration region shows that the second hydration shell is trigonally bound to the first shell water with a high degree of agreement between the AIMD and AIMD/MM calculations. Beyond the second shell, the bonding pattern returns to the tetrahedral structure of bulk water. The AIMD/MM results emphasize the importance of a quantum description of the first hydration shell to correctly describe the hydration region. In these calculations the full d(10) electronic structure of the valence shell of the Zn(2+) ion is retained. The simulations show substantial and complex charge relocation on both the Zn(2+) ion and the first hydration shell. The dipole moment of the waters in the first hydration shell is 3.4 D (3.3 D AIMD/MM) versus 2.73 D bulk. Little polarization is found for the waters in the second hydration shell (2.8 D). No exchanges were seen between the first and the second hydrations shells; however, many water transfers between the second hydration shell and the bulk were observed. For 64 waters, the AIMD and AIMD/MM simulations give nearly identical results for exchange dynamics. However, in the larger particle simulations (256 waters) there is a significant reduction in the second shell to bulk exchanges.

We applied the self-consistent polarization density functional theory (SCP-DFT) to water. SCP-DFT requires only minimal parametrization, self-consistently includes the dispersion interaction neglected by standard DFT functionals, and has a cost similar to standard DFT despite its improved performance. Compared to the DFT functionals BLYP and BLYP-D (where the latter contains a simple dispersion correction), SCP-DFT yields interaction energies per molecule and harmonic frequencies of clusters in better agreement with experiment, with errors in the former of only a few tenths of a kcal/mol. BLYP and BLYP-D underbind and overbind the clusters, respectively, by up to about 1 kcal/mol. For liquid water, both BLYP and SCP-DFT predict radial distribution functions that are similar and overstructured compared to experiment. However, SCP-DFT improves over BLYP in predicting the experimental enthalpy of vaporization. A decomposition of the dimer interaction energy attempts to rationalize the performance of SCP-DFT. The SCP-DFT approach holds promise as an efficient and accurate method for describing large hydrogen-bonded systems, and has the potential to model complex systems with minimal parametrization.

Incoherent inelastic neutron scattering is used to probe the effects of dihydrogen bonding on the vibrational dynamics in the molecular crystal of ammonia borane. The thermal neutron energy loss spectra of (11)B enriched ammonia borane isotopomers ((11)BH(3)NH(3), (11)BD(3)NH(3), and (11)BH(3)ND(3)) are presented and compared to the vibrational power spectrum calculated using ab initio molecular dynamics. A harmonic vibrational analysis on NH(3)BH(3) clusters was also explored to check for consistency with experiment and the power spectrum. The measured neutron spectra and computed ab initio power spectrum compare extremely well (50-500 cm(-1)). Some assignment of modes to simple harmonic motion, e.g., NH(3) and BH(3) torsion in the molecular crystal is possible, and it is confirmed that the lowest modes are dominated by collective motion. We show that the vibrational dynamics as modeled with ab initio molecular dynamics provides a more complete description of anharmonic and collective dynamics in the low frequency region of the inelastic incoherent neutron scattering spectra when compared to the conventional harmonic approach.

Semiempirical self-consistent field (SCF) methods based on the neglect of diatomic differential overlap (NDDO) formalism have the ability to treat the formation and breaking of chemical bonds but have been found to poorly describe hydrogen bonding and weak electrostatic complexes. In contrast, most empirical potentials are not able to describe bond breaking and formation but have the ability to add missing elements of hydrogen bonding by using classical electrostatic interactions. We present a new method which combines aspects of both NDDO-based SCF techniques and classical descriptions of polarization to describe the diffuse nature of the electronic wavefunction in a self-consistent manner. We develop the "self-consistent polarization neglect of diatomic differential overlap" (SCP-NDDO) theory with the additional description of molecular dispersion developed as a second-order perturbation theory expression. The current study seeks to model water-water interactions as a test case. To this end, we have parametrized the method to accurate ab initio complete basis set limit estimates of small water cluster binding energies of Xantheas and co-workers [J. Chem. Phys. 116, 1493 (2002); 120, 823 (2004)]. Overall agreement with the ab initio binding energies (n=2-6, and 8) is achieved with a rms error of 0.19 kcal/mol. We achieve noticeable improvements in the structure, vibrational frequencies, and energetic predictions of water clusters (n< or =21) relative to standard NDDO-based methods.

First principles molecular dynamics simulations of the hydration shells surrounding UO(2)(2+) ions are reported for temperatures near 300 K. Most of the simulations were done with 64 solvating water molecules (22 ps). Simulations with 122 water molecules (9 ps) were also carried out. The hydration structure predicted from the simulations was found to agree with very well-known results from x-ray data. The average U=O bond length was found to be 1.77 A. The first hydration shell contained five trigonally coordinated water molecules that were equatorially oriented about the O-U-O axis with the hydrogen atoms oriented away from the uranium atom. The five waters in the first shell were located at an average distance of 2.44 A (2.46 A, 122 water simulation). The second hydration shell was composed of distinct equatorial and apical regions resulting in a peak in the U-O radial distribution function at 4.59 A. The equatorial second shell contained ten water molecules hydrogen bonded to the five first shell molecules. Above and below the UO(2)(2+) ion, the water molecules were found to be significantly less structured. In these apical regions, water molecules were found to sporadically hydrogen bond to the oxygen atoms of the UO(2)(2+), oriented in such a way as to have their protons pointed toward the cation. While the number of apical waters varied greatly, an average of five to six waters was found in this region. Many water transfers into and out of the equatorial and apical second solvation shells were observed to occur on a picosecond time scale via dissociative mechanisms. Beyond these shells, the bonding pattern substantially returned to the tetrahedral structure of bulk water.

We use ab initio electronic structure methods to calculate the many-body decomposition of the binding energies of the OH.(H2O)n (n=2,3) complexes. We employ MP2 and CCSD(T) levels of theory with aug-cc-pVDZ and aug-cc-pVTZ basis sets and analyze the significance of the nonpairwise interactions between OH radical and the surrounding water molecules. We also evaluate the accuracy of our newly developed potential function, the modified Thole-type model, for predicting the many-body terms in these complexes. Our analysis of the many-body contributions to the OH.(H2O)n binding energies clearly shows that they are just as important in the OH interactions with water as they are for interactions in pure water systems.

Activation energies for water cluster evaporation are of interest in many areas of chemical physics. We present the first computation of activation energies for monomer evaporation of small water clusters using the formalism of dynamical nucleation theory (DNT). To this end, individual evaporation rate constants are computed for water clusters (H(2)O)(i), where i=2-10 for temperatures ranging from 243 to 333 K. These calculations employ a parallel sampling technique utilizing a Global Arrays toolkit. The resulting evaporation rate constants for each cluster are then fitted to Arrhenius equations to obtain activation energies. We discuss DNT evaporation rate constants and their relation to potentials of mean force, activation energies, and how to account for nonseparability of the reaction coordinate in the reactant state partition function.

Raman spectra of single crystal ammonia borane, NH3BH3, were recorded as a function of temperature from 88 to 300 K using Raman microscopy and a variable temperature stage. The orthorhombic to orientationally disordered tetragonal phase transition at 225 K was clearly evident from the decrease in the number of vibrational modes. However, some of the modes in the orthorhombic phase appeared to merge 10-12 K below the phase transition perhaps suggesting the presence of an intermediate phase. Factor group analysis of vibrational spectra for both orthorhombic and tetragonal phase is provided. In addition, electronic structure calculations are used to assist in the interpretation and assignment of the normal modes.

We present an approach to calculate the free energy profile along a condensed-phase reaction path based on high-level electronic structure methods for the reactive region. The bulk of statistical averaging is shifted toward less expensive descriptions by using a hierarchy of representations that includes molecular mechanics, density functional theory, and coupled cluster theories. As an application of this approach we study the reaction of CHCl3 with OH- in aqueous solution.

In recent years, the ClO free radical has been found to play an important role in the ozone removal processes in the atmosphere. In this work, the authors present a potential energy surface scan of the ClO.H2O system with high-level ab initio methods. Because of the existence of low-lying excited states of the ClO.H2O complex and their potential impact on the chemical behavior of the ClO radical in the atmosphere, the authors perform the potential energy surface scan at the CCSD(T)/aug-cc-pVTZ level of theory of both the first excited and ground states. Analytical potentials for both ground and excited states, with the ClO and H2O units held fixed at their optimized geometries and with anisotropic atomic polarizabilities modeling the physics of the unpaired electron in the ClO radical, were built based on a Thole-type model. The two minima of the ClO.H2O complex are recovered by the analytical potential.

We report converged quantum statistical mechanical simulations of liquid water with the Thole-type Model (version 2.1), Flexible, polarizable (TTM2.1-F) interaction potential for water. Simulations of total length of 600 ps with a 0.05 fs time step for a periodic unit cell of 256 molecules with up to 32 replicas per atom suggest that the quantum effects contribute 1.01+/-0.02 kcal/mol to the liquid enthalpy of formation at 298.15 K. They furthermore demonstrate for the first time a quantitative agreement with experiment for the heights and broadening of the intramolecular OH and HH peaks in the radial distribution functions.

The OH radical is one of the most important oxidants in the atmosphere due to its high reactivity. The study of hydrogen-bonded complexes of OH with the water molecules is a topic of significant current interest. In this work, we present the development of a new analytical functional form for the interaction potential between the rigid OH radical and H(2)O molecules. To do this we fit a selected functional form to a set of high level ab initio data. Since there is a low-lying excited state for the H(2)O.OH complex, the impact of the excited state on the chemical behavior of the OH radical can be very important. We perform a potential energy surface scan using the CCSD(T)/aug-cc-pVTZ level of electronic structure theory for both excited and ground states. To model the physics of the unpaired electron in the OH radical, we develop a tensor polarizability generalization of the Thole-type all-atom polarizable rigid potential for the OH radical, which effectively describes the interaction of OH with H(2)O for both ground and excited states. The stationary points of (H(2)O)(n)OH clusters were identified as a benchmark of the potential.

The Gibbs ensemble Monte Carlo simulation technique was used to compare the phase equilibria of the rigid TIP4P water model [Jorgensen et al., J. Chem. Phys. 79, 926 (1983)] utilizing classical and quantum statistical mechanics. The quantum statistical mechanical treatment generally resulted in lower liquid densities and higher vapor densities, narrowing the phase envelope. As a result, the calculated critical temperatures and normal boiling points were lower from the quantum simulations than the classical by 22 and 17 K, respectively, but the critical densities were equal within the estimated uncertainties. When the phase diagram from the quantum statistical mechanical treatment was increased by 22 K, it agreed with the classical results quite well throughout the entire simulated temperature range. A semiclassical treatment, involving a low order expansion in Planck's constant, resulted in good agreement with the path integral results for second virial coefficients, but gave densities and vapor pressures that fluctuated between the values for the classical and quantum statistical mechanics values, with no definite agreement with either.

Vapor to liquid multicomponent nucleation is a dynamical process governed by a delicate interplay between condensation and evaporation. Since the population of the vapor phase is dominated by monomers at reasonable supersaturations, the formation of clusters is governed by monomer association and dissociation reactions. Although there is no intrinsic barrier in the interaction potential along the minimum energy path for the association process, the formation of a cluster is impeded by a free energy barrier. Dynamical nucleation theory provides a framework in which equilibrium evaporation rate constants can be calculated and the corresponding condensation rate constants determined from detailed balance. The nucleation rate can then be obtained by solving the kinetic equations. The rate constants governing the multistep kinetics of multicomponent nucleation including sensitivity analysis and the potential influence of contaminants will be presented and discussed.

We construct a rigid-body (five-dimensional) potential-energy surface for the water-hydrogen complex using scaled perturbation theory (SPT). An analytic fit of this surface is obtained, and, using this, two minima are found. The global minimum has C2v symmetry, with the hydrogen molecule acting as a proton donor to the oxygen atom on water. A local minimum with Cs symmetry has the hydrogen molecule acting as a proton acceptor to one of the hydrogen atoms on water, where the OH bond and H2 are in a T-shaped configuration. The SPT global minimum is bound by 1097 microEh (Eh approximately 4.359744 x 10(-18) J). Our best estimate of the binding energy, from a complete basis set extrapolation of coupled-cluster calculations, is 1076.1 microEh. The fitted surface is used to calculate the second cross virial coefficient over a wide temperature range (100-3000 K). Three complementary methods are used to quantify quantum statistical mechanical effects that become significant at low temperatures. We compare our results with experimental data, which are available over a smaller temperature range (230-700 K). Generally good agreement is found, but the experimental data are subject to larger uncertainties., ((c) 2004 American Institute of Physics)