Rotational excitation of water by hydrogen molecules: comparison of results from classical and quantum mechanics

Quasiclassical trajectory calculations are carried out for rotational excitation of water by hydrogen molecules. State-to-state rate coefficients are determined at 100 K and are compared to available quantum results. A good agreement between classical and quantum rates is observed for downward transitions, with an average accuracy of classical results better than a factor of 2. It is thus found that the ambiguities described by Faure and Wiesenfeld [J. Chem. Phys. 121, 6771 (2004)] can be solved in the particular case of waterlike asymmetric-top molecules.     

Liquid-vapor and liquid-liquid phase equilibria of the Brodholt-Sampoli-Vallauri polarizable water model

Liquid-vapor and liquid-liquid phase equilibria of the polarizable Brodholt-Sampoli-Vallauri water model have been investigated by Gibbs ensemble Monte Carlo computer simulations. The coexisting liquid and vapor densities and energy of vaporization of the model is found to be in a reasonable agreement with experimental data in the entire temperature range of liquid-vapor coexistence. The critical temperature and density of the model are found to be 615 K and 0.278 gcm(3), respectively, close to the experimental values of 647.1 K and 0.322 gcm(3). In the supercooled state two distinct liquid-liquid coexistence regions are observed. The existence of liquid-liquid phase separation of a polarizable water model is demonstrated for the first time.     

Theory and phenomenology for a variety of classical and quantum phase transitions

In this review, we first introduce recent progress in the mathematical structure of the three-dimensional Ising model, from the points of view of topologic, algebraic and geometric aspects. Then we discuss in turn Anderson localization due to disorder and then first- and second-order metal-insulator transitions, depending on electron correlation, with and without a magnetic field. Finally, we make intimate contact with the phase diagram showing the equilibrium between low temperature regimes of the magnetically induced Wigner electron solid and the so-called Laughlin electron liquid in the two-dimensional case.     

Influence of bond flexibility on the vapor-liquid phase equilibria of water

The authors performed Gibbs ensemble simulations on the vapor-liquid equilibrium of water to investigate the influence of incorporating intramolecular degrees of freedom in the simple point charge (SPC) water model. Results for vapor pressures, saturation densities, heats of vaporization, and the critical point for two different flexible models are compared with data for the corresponding rigid SPC and SPC/E models. They found that the introduction of internal vibrations, and also their parametrization, has an observable effect on the prediction of the vapor-liquid coexistence curve. The flexible SPC/Fw model, although optimized to describe bulk diffusion and dielectric constants at ambient conditions, gives the best prediction of saturation densities and the critical point of the examined models.     

Ground-state properties of the retinal molecule: from quantum mechanical to classical mechanical computations of retinal proteins

Retinal proteins are excellent systems for understanding essential physiological processes such as signal transduction and ion pumping. Although the conjugated polyene system of the retinal chromophore is best described with quantum mechanics, simulations of the long-timescale dynamics of a retinal protein in its physiological, flexible, lipid-membrane environment can only be performed at the classical mechanical level. Torsional energy barriers are a critical ingredient of the classical force-field parameters. Here we review briefly current retinal force fields and discuss new quantum mechanical computations to assess how the retinal Schiff base model and the approach used to derive the force-field parameters may influence the torsional potentials.     

An NMR and quantum mechanical investigation of solvent effects on conformational equilibria of butanedinitrile

Vicinal proton-proton NMR couplings and ab initio quantum mechanics have been used to investigate solvent effects on conformational equilibria of butanedinitrile. The trans and gauche conformations are about equally favored at room temperature in solvents of low dielectric constant while the equilibrium is essentially the statistical proportions of one-third trans and two-thirds gauche in water with a high dielectric constant. The coupling assignments were confirmed with the aid of stereospecific deuterium-labeled (R,R or S,S)-1,2-dideuteriobutanedinitrile. The calculations support the observed trends. Similar results were observed for 1,2-dibromo- and dichloroethanes.     

Fluid phase equilibria in the system n-butane + water

Although the phase equilibria in the system n-butane + water have been studied frequently, a review of the experimental results has revealed serious disagreement among the various investigators. In this work, the data from the literature are supplemented with some new solubility data. These data are then used construct a model, based on Henry's law, for the phase equilibria.     

Effect of noise on the classical and quantum mechanical nonlinear response of resonantly coupled anharmonic oscillators

Multidimensional infrared spectroscopy probes coupled molecular vibrations in complex, condensed phase systems. Recent theoretical studies have focused on the analytic structure of the nonlinear response functions required to calculate experimental observables in a perturbative treatment of the radiation-matter interaction. Classical mechanical nonlinear response functions have been shown to exhibit unbounded growth for anharmonic, integrable systems, as a consequence of the nonlinearity of classical mechanics, a feature that is absent in a quantum mechanical treatment. We explore the analytic structure of the third-order vibrational response function for an exactly solvable quantum mechanical model that includes some of the important and theoretically challenging aspects of realistic models of condensed phase systems: anharmonicity, resonant coupling, fluctuations, and a well-defined classical mechanical limit.     

The phase diagram of water from quantum simulations

The phase diagram of water has been calculated from the TIP4PQ/2005 model, an empirical rigid non-polarisable model. The path integral Monte Carlo technique was used, permitting the incorporation of nuclear quantum effects. The coexistence lines were traced out using the Gibbs-Duhem integration method, once having calculated the free energies of the liquid and solid phases in the quantum limit, which were obtained via thermodynamic integration from the classical value by scaling the mass of the water molecule. The resulting phase diagram is qualitatively correct, being displaced to lower temperatures by 15-20 K. It is found that the influence of nuclear quantum effects is correlated to the tetrahedral order parameter.     

Finite element interpolation for combined classical/quantum mechanical molecular dynamics simulations

A method is presented to interpolate the potential energy function for a part of a system consisting of a few degrees of freedom, such as a molecule in solution. The method is based on a modified finite element (FE) interpolation scheme. The aim is to save computer time when expensive methods such as quantum-chemical calculations are used to determine the potential energy function. The expensive calculations are only carried out if the molecule explores new unknown regions of the conformation space. If the molecule resides in regions previously explored, a cheap interpolation is performed instead of an expensive calculation, using known neighboring points. We report the interpolation techniques for the energies and the forces of the molecule, the handling of the FE mesh, and an application to a simple test example in molecular dynamics (MD) simulations. Good performance of the method was obtained (especially for MD simulations with a preceding Monte Carlo mesh generation) without losing accuracy. © 1997 John Wiley & Sons, Inc. J Comput Chem 18 : 1484–1495, 1997     

Critical examination of Ga2Se3 phase equilibria

The Ga₂Se₃ phase equilibria were examined by DTA, X-ray diffractometry, visual observation, and isothermal crystal growth to resolve discrepancies in the reported binary phase diagrams. The Se-rich solid solubility extends to less than 60.5 at. % Se. There is no peritectic transformation at 880°C, nor is there a liquid immiscibility from 75 to 85 at. % Se. The rate of crystal growth influences the incorporation of lattice defects and thus determines which polymorphic form of Ga₂Se₃ crystallizes from the melt. The preferred liquidus curve from 60 to 100 at. % Se is presented and compared to those reported in the literature.     

Environmental effects on the spectroscopic properties of gallic acid: a combined classical and quantum mechanical study

The solvation of gallic acid (in water and acetonitrile) is studied by means of its spectroscopic properties. IR, UV, and NMR spectra are predicted by using various solvation models obtained in terms of both purely classical and density functional approaches. Comparison with experiments is used to validate solvation models. Hydrogen-bond and long-range (or bulk) effects are evaluated by comparing different solvation models. A continuum-only approach, a purely discrete, and a mixed continuum/discrete approach based on quantum-mechanical and classical molecular-dynamics solute-solvent clusters are tested.     

Single molecule charge transport: from a quantum mechanical to a classical description

This paper explores charge transport at the single molecule level. The conductive properties of both small organic molecules and conjugated polymers (molecular wires) are considered. In particular, the reasons for the transition from fully coherent to incoherent charge transport and the approaches that can be taken to describe this transition are addressed in some detail. The effects of molecular orbital symmetry, quantum interference, static disorder and molecular vibrations on charge transport are discussed. All of these effects must be taken into account (and may be used in a functional way) in the design of molecular electronic devices. An overview of the theoretical models employed when studying charge transport in small organic molecules and molecular wires is presented.     

Inelastic collisions of molecular hydrogen: a comparison of results from quantum and classical mechanics

Full-dimensional quantum and classical calculations have been carried out for inelastic (nonreactive) energy transfer in H2+H2 on the ab initio potential energy surface of Boothroyd et al. [J. Chem. Phys. 116, 666 (2002)]. State-to-state cross sections are determined and compared for transitions from H2(0,j(ab))+H2(1,j(cd)). While there is excellent agreement for transitions involving small Deltaj, for larger Deltaj and for vibrational relaxation, significant differences are observed which exhibit no systematic trends. Reasons for this disagreement are discussed.     

Cumulative reaction probabilities: A comparison between quasiclassical and quantum mechanical results

This article presents a quasiclassical trajectory (QCT) method for determining the cumulative reaction probability (CRP) as a function of the total energy. The method proposed is based on a discrete sampling using integer values of the total and orbital angular momentum quantum numbers for each trajectory and on the development of equations that have a clear counterpart in the quantum mechanical (QM) case. The calculations comprise cumulative reaction probabilities at a given total angular momentum J, as well as those summed over J. The latter are used to compute QCT rate constants. The method is illustrated by comparing QCT and exact QM results for the H+H2, H+D2, D+H2, and H+HD reactions. The agreement between QCT and QM results is very good, with small discrepancies between the two data sets indicating some genuine quantum effects. The most important of these involves the value of the CRP at low energies which, due to the absence of tunneling, is lower in the QCT calculations, causing the corresponding rate constants to be smaller. The second is the steplike structure that is clearly displayed in the QM CRP for J = 0, which is much smoother in the corresponding QCT results. However, when the QCT density of reactive states, i.e., the derivatives of the QCT CRP with respect to the energy, is calculated, a succession of maxima and minima is obtained which roughly resembles those found in the QM calculations, although the latter are considerably sharper. The analysis of the broad peaks in the QCT density of reactive states indicates that the distributions of collision times associated with the maxima are somewhat broader, with a tail extending to larger collision times, than those associated with the minima. In addition, the QM and QCT dynamics of the isotopic variants mentioned above are compared in the light of their CRPs. Issues such as the compliance of the QCT CRP with the law of microscopic reversibility, as well as the similarity between the CRPs for ortho and para species in the QM and QCT cases, are also addressed.     

Structure and phase equilibria of mixtures of the complex salt hexadecyltrimethylammonium polymethacrylate, water and different oils

This work reports on phase diagrams for mixtures of a complex salt formed by a cationic surfactant and an oppositely charged polyelectrolyte, hexadecyltrimethylammonium polymethacrylate, in binary mixtures with water and in ternary mixtures containing water and organic solvents of different polarity ('oils'): decanol, octanol, p-xylene and cyclohexane. The liquid crystalline structures formed were identified by small angle X-ray scattering measurements, which also provided information about changes in the size of the aggregates as a function of the system composition. These results are analysed in comparison with others previously reported [Bernardes et al., J. Phys. Chem. B 110 (2006) 10332-10340] for the analog complex formed with polyacrylate and, in general, reveal that the presence of an extra methylene group in the polymer chain does not produce significant changes in the complex phase diagrams nor in the structure of the liquid crystalline phases formed. Additionally, the obtained results confirm once more the approach used to analyze these kinds of systems formed by polymer and oppositely charged surfactant.     

Geometrical interpretation of crystal-vapor phase equilibria and equilibria in helium systems

A geometrical interpretation of phase equilibria was suggested and used to determine several rules governing crystal-vapor equilibria, analytically describe the transition of helium systems into the superfluid state, and reveal the general properties of the low-temperature state of matter.     

Solid and liquid phase equilibria and solid-hydrate formation in binary mixtures of water with amines

Solid and liquid phase diagrams have been constructed for {water+triethylamine,or+N,N-dimethylformamide(DMF) or+N,N-dimethlacetamide (DMA)} Solid-hydrates form with the empirical formulae N(C2H5)3 3H2O,DMF 3H2O,DMF 2H2O,DMA 3H2O and (DMA)2 3H2O.All are congruently melting except the first which melts incongruently.The solid-hydrate formation is attributed to hydrogen bond.The results are compared with the references     

Dimerisation of urea in water solution: a quantum mechanical investigation

The effect of water solvation on the structure and stability of cyclic dimers of urea has been investigated with the aid of density functional theory at the B3LYP/6-311++G** level. Several hydration models have been discussed. Specific solvent effects have been simulated through single and multiple water-urea interactions involving all the hydration sites of urea. The bulk solvent effects have been estimated through polarised continuum models. Under all the hydration patterns cyclic dimers continue to be stable structures although the solvent weakens the urea-urea interaction. Single and multiple specific urea-water interactions are competitive with urea dimerisation. The anticooperative nature of the two intermolecular interactions is largely due to the changes on sigma- and pi-electron density of urea caused by hydrogen bonding with water. The stability of the dimer is however, lost within a few ps when the hydrated dimer is described by a quantum mechanical molecular dynamics approach (ADMP). The cyclic dimer evolves towards structures where urea molecules are linked not more directly but through water molecules which have a bridge function.