Density scaling of the transport properties of molecular and ionic liquids

Casalini and Roland [Phys. Rev. E 69, 062501 (2004); J. Non-Cryst. Solids 353, 3936 (2007)] and other authors have found that both the dielectric relaxation times and the viscosity, η, of liquids can be expressed solely as functions of the group (TV?(γ)), where T is the temperature, V is the molar volume, and γ a state-independent scaling exponent. Here we report scaling exponents γ, for the viscosities of 46 compounds, including 11 ionic liquids. A generalization of this thermodynamic scaling to other transport properties, namely, the self-diffusion coefficients for ionic and molecular liquids and the electrical conductivity for ionic liquids is examined. Scaling exponents, γ, for the electrical conductivities of six ionic liquids for which viscosity data are available, are found to be quite close to those obtained from viscosities. Using the scaling exponents obtained from viscosities it was possible to correlate molar conductivity over broad ranges of temperature and pressure. However, application of the same procedures to the self-diffusion coefficients, D, of six ionic and 13 molecular liquids leads to superpositioning of poorer quality, as the scaling yields different exponents from those obtained with viscosities and, in the case of the ionic liquids, slightly different values for the anion and the cation. This situation can be improved by using the ratio (D∕T), consistent with the Stokes-Einstein relation, yielding γ values closer to those of viscosity.     

Communication: Thermodynamic scaling of the Debye process in primary alcohols

The molecular dynamics of hydrogen-bonded liquids usually does not satisfy the thermodynamic scaling. However, very recently, two opposite conclusions about validity of thermodynamical scaling in monohydroxy alcohol, 2-ethyl-1-hexanol, were presented by Reiser et al. [J. Chem. Phys. 132, 181101 (2010)] and Fragiadakis et al. [J. Chem. Phys. 132, 144505 (2010)]. In this communication we present new experimental results that can explain this ostensible contradiction.     

Thermodynamical and structural properties of imidazolium based ionic liquids from molecular simulation

We have performed molecular dynamics simulations to determine the densities and heat of vaporization as well as structural information for the 1-alkyl-3-methyl-imidazolium based ionic liquids [amim][Cl] and [amim][BF(4)] in the temperature range from 298 to 363 K. In this simulation study, we used an united atom model of Liu et al. [Phys. Chem. Chem. Phys. 8, 1096 (2006)] for the [emim(+)] and [bmim(+)] cations, which we have extended for simulation in [hmim]-ILs and combined with parameters of Canongia Lopes et al. [J. Phys. Chem. B 108, 2038 (2004)] for the [Cl(-)] anion. Our simulation results prove that both the original united atoms approach by Liu et al. and our extension yield reasonable predictions for the ionic liquid with a considerably reduced computational expense than that required for all atoms models. Radial distribution functions and spatial distribution functions where employed to analyze the local structure of this ionic liquid, and in which way it is influenced by the type of the anion, the size of the cation, and the temperature. Our simulations give evidence for the occurrence of tail aggregations in these ionic liquids with increasing length of the side chain and also increasing temperature.     

Optimization of the anisotropic united atoms intermolecular potential for n-alkanes: improvement of transport properties

The parameters of the anisotropic united atom (AUA) intermolecular potential for n-alkanes originally proposed by Toxvaerd [J. Chem. Phys. 93, 4290 (1990)] [AUA(3)] was optimized by Ungerer et al. [J. Chem. Phys. 112, 5499 (2000)] [AUA(4)] on the basis of equilibrium properties (vapor pressures, vaporization enthalpies, and liquid densities). In this work we analyze the influence of the torsion potential in the internal and collective dynamics of the AUA model. The modified potential [AUA(4m)] preserves all the intermolecular parameters and only explores an increment in the trans-gauche and gauche(+)-gauche(-) transition barrier of the torsion potential. This modification better reproduce different transport properties (shear viscosity, self-diffusion coefficient, and internal relaxation times), keeping the accuracy achieved in our previous work for equilibrium properties. An extensive investigation of the shear viscosity of ethane, n-pentane, n-dodecane, and n-eicosane in a wide range of pressures and temperatures shows that the AUA(4m) improves the accuracy of the original AUA(4), reducing the absolute average deviation from 30% to 14.5%. Finally, the self-diffusion coefficient of n-hexane computed with the new model in the range of 223-333 K and from 0.1 to 295 MPa is in better agreement with respect to the experimental data than the original model.     

Polarizable and flexible model for ethanol

A polarizable, flexible model for ethanol is obtained based on an extensive series of B3LYP/6-311++G(d,p) calculations and molecular dynamics simulations. The ethanol model includes electric-field dependence in both the atomic charges and the intramolecular degrees of freedom. Field-dependent intramolecular potentials have been attempted only once previously, for OH and HH stretches in water [P. Cicu et al., J. Chem. Phys. 112, 8267 (2000)]. The torsional potential involving the hydrogen-bonding hydrogen in ethanol is found to be particularly field sensitive. The methodology for developing field-dependent potentials can be readily generalized to other molecules and is discussed in detail. Molecular dynamics simulations of bulk ethanol are performed and the results are assessed based on comparisons with the self-diffusion coefficient [N. Karger et al., J. Chem. Phys. 93, 3437 (1990)], dielectric constant [J. T. Kindt and C. A. Schmuttenmaer, J. Phys. Chem. 100, 10373 (1996)], enthalpy of vaporization [R. C. Wilhoit and B. J. Zwolinski, J. Phys. Chem. Ref. Data, Suppl. 2, 2 (1973)], and experimental interatomic distributions [C. J. Benmore and Y. L. Loh, J. Chem. Phys. 112, 5877 (2000)]. The simultaneous variation of the atomic charges and the intramolecular potentials requires modified equations of motion and a multiple time step algorithm has been implemented to solve these equations. The article concludes with a discussion of the bulk structure and properties with an emphasis on the hydrogen bonding network.     

Prevalence for the universal distribution of relaxation times near the glass transitions in experimental model systems: rodlike liquid crystals and orientationally disordered crystals

Recently, Nielsen et al. [J. Chem. Phys. 130, 154508 (2009); Philos. Mag. 88, 4101 (2008)] demonstrated a universal pattern for the high frequency wing of the loss curve for primary relaxation time on approaching the glass transition for organic liquids. In this contribution it is presented that a similar universality occurs for glass-forming liquid crystals and orientationally disordered crystals (plastic crystals). Empirical correlations of the found behavior are also briefly discussed.     

Higher excitations for an exponential multireference wavefunction Ansatz and single-reference based multireference coupled cluster Ansatz: application to model systems H(4), P(4), and BeH(2)

This article reports on the convergence of the exponential multireference wavefunction Ansatz (MRexpT) [J. Chem. Phys. 123, 84102 (2005)] and the single-reference based multireference coupled cluster Ansatz [J. Chem. Phys. 94, 1229 (1991)] with respect to higher cluster excitations. The approaches are applied to the H(4), P(4), and BeH(2) model systems according to the recently published analysis by Evangelista et al. [J. Chem. Phys. 125, 154113 (2006)]. The results show both MRexpT and SRMRCC to be highly accurate although SRMRCC shows problems due to its lack of Fermi vacuum invariance (symmetry breaking).     

A fully simultaneously optimizing genetic approach to the highly excited coupled-cluster factorization problem

In this article we report on the coupled-cluster factorization problem. We describe the first implementation that optimizes (i) the contraction order for each term, (ii) the identification of reusable intermediates, (iii) the selection and factoring out of common factors simultaneously, considering all projection levels in a single step. The optimization is achieved by means of a genetic algorithm. Taking a one-term-at-a-time strategy as reference our factorization yields speedups of up to 4 (for intermediate excitation levels, smaller basis sets). We derive a theoretical lower bound for the highest order scaling cost and show that it is met by our implementation. Additionally, we report on the performance of the resulting highly excited coupled-cluster algorithms and find significant improvements with respect to the implementation of Kállay and Surján [J. Chem. Phys. 115, 2945 (2001)] and comparable performance with respect to MOLPRO's handwritten and dedicated open shell coupled cluster with singles and doubles substitutions implementation [P. J. Knowles, C. Hampel, and H.-J. Werner, J. Chem. Phys. 99, 5219 (1993)].     

Close-coupling study of rotational energy transfer of CO (upsilon=2) by collisions with He atoms

Quantum close-coupling scattering calculations of rotational energy transfer in the vibrationally excited CO due to collisions with He atom are presented for collision energies between 10(-5) and approximately 1000 cm-1 with CO being initially in the vibrational level upsilon=2 and rotational levels j=0,1,4, and 6. The He-CO interaction potential of Heijmen et al. [J. Chem. Phys. 107, 9921 (1997)] was adopted for the calculations. Cross sections for rovibrational transitions and state-to-state rotational energy transfer from selected initial rotational levels were computed and compared with recent measurements of Carty et al. [J. Chem. Phys. 121, 4671 (2004)] and available theoretical results. Comparison in all cases is found to be excellent, providing a stringent test for the scattering calculations as well as the reliability of the He-CO interaction potential by Heijmen et al.     

Comment on "Study of dielectric relaxations of anhydrous trehalose and maltose glasses" [J. Chem. Phys. 134, 014508 (2011)

Very recently Kwon et al. [H.-J. Kwon, J.-A. Seo, H. K. Kim, and Y. H. Hwang, J. Chem. Phys. 134, 014508 (2011)] published an article on the study of dielectric relaxation in trehalose and maltose glasses. They carried out broadband dielectric measurements at very wide range of temperatures covering supercooled liquid as well as glassy state of both saccharides. It is worth to mention that authors have also applied a new method for obtaining anhydrous glasses of trehalose and maltose that enables avoiding their caramelization. Four relaxation processes were identified in dielectric spectra of both saccharides. The slower one was identified as structural relaxation process the next one, not observed by the others, was assigned as Johari-Goldstein (JG) β-relaxation, while the last two secondary modes were of the same nature as found by Kaminski et al. [K. Kaminski, E. Kaminska, P. Wlodarczyk, S. Pawlus, D. Kimla, A. Kasprzycka, M. Paluch, J. Ziolo, W. Szeja, and K. L. Ngai, J. Phys. Chem. B 112, 12816 (2008)]. In this comment we show that the authors mistakenly assigned the slowest relaxation process as structural mode of disaccharides. We have proven that this relaxation process is an effect of formation of thin layer of air or water between plate of capacitor and sample. The same effect can be observed if plates of capacitor are oxidized. Thus, we concluded that their slowest mode is connected to the dc conduction process while their β JG process is primary relaxation of trehalose and maltose.     

Comment on: "Disentangling density and temperature effects in the viscous slowing down of glass forming liquids" [J. Chem. Phys. 120, 6135 (2004)

Recently, Tarjus et al. [G. Tarjus, D. Kivelson, S. Mossa, and C. Alba-Simionesco, J. Chem. Phys. 120, 6135 (2004)] concluded from a review of data for a variety of glass formers that the supercooled dynamics are almost invariably dominated by temperature T, rather than by density rho. By including additional published data into such a compilation, we show that for van der Waals molecular liquids, the dynamics near T(g) are in fact governed as much by density as by temperature. Moreover, relaxation times measured at various temperatures and pressures can be superimposed by plotting as a function rho(gamma)/T. This scaling form can arise from an assumed inverse power law for the intermolecular repulsive potential, with gamma a material constant. (c) 2004 American Institute of Physics.     

Reactivity enhancement of ultracold O(3P)+H2 collisions by van der Waals interactions

The role of van der Waals forces in O((3)P)+H(2)(upsilon=1,j=0) collisions is investigated theoretically at low and ultralow temperatures. Quantum scattering calculations have been performed for zero total angular momentum using the lowest London-Eyring-Polanyi-Sato double-polynomial (3)A(") potential-energy surface reported by [Rogers et al., J. Phys. Chem. A 104, 2308 (2000)] and its recent BMS1 and BMS2 extensions developed by [Brandao et al., J. Chem. Phys. 121, 8861 (2004)] which provide a more accurate treatment of the van der Waals interaction. Our calculations show that van der Waals forces strongly influence chemical reactivity at ultracold translational energies. The presence of a zero-energy resonance for the BMS1 surface is found to enhance reactivity in the ultracold regime and shift the Wigner threshold to lower temperatures.     

Cross sections for electron impact excitation of the vibrationally resolved A 1Pi electronic state of carbon monoxide

The authors report new differential cross section measurements for electron impact excitation of the A (1)Pi(v(')) states of carbon monoxide. The energy range is 20-200 eV. They also reanalyze the A (1)Pi(v(')) manifold cross sections of Middleton et al. [J. Phys. B 26, 1743 (1993)] in order to provide a basis for comparison with our new vibrationally resolved differential cross sections. Excellent agreement is found between the two sets of measurements at all common energies. From 20 to 200 eV the present differential cross sections are extrapolated and integrated, and the corresponding integral excitation cross sections determined. New scaled Born integral cross sections, calculated as a part of the present study, are compared against these experimental integral cross sections, with excellent agreement being found for all the A (1)Pi(v(')=0-7)<--X (1)Sigma(g) (+)(v(")=0) transitions. In addition our scaled Born integral cross sections are found to be in excellent agreement between 300 and 1500 eV with those derived from the previous experiments of Lassettre and Skerbele [J. Chem. Phys. 54, 1597 (1971)] and of Zhong et al. [Phys. Rev. A 55, 1799 (1997)] and from near threshold to 15 eV with those derived from Zobel et al. [J. Phys. B 29, 813 (1996)] and Zetner et al. (J. Phys. B 31, 2395 (1998)].     

Comparison of different methods for calculating the paramagnetic relaxation enhancement of nuclear spins as a function of the magnetic field

The enhancement of the spin-lattice relaxation rate for nuclear spins in a ligand bound to a paramagnetic metal ion [known as the paramagnetic relaxation enhancement (PRE)] arises primarily through the dipole-dipole (DD) interaction between the nuclear spins and the electron spins. In solution, the DD interaction is modulated mostly by reorientation of the nuclear spin-electron spin axis and by electron spin relaxation. Calculations of the PRE are in general complicated, mainly because the electron spin interacts so strongly with the other degrees of freedom that its relaxation cannot be described by second-order perturbation theory or the Redfield theory. Three approaches to resolve this problem exist in the literature: The so-called slow-motion theory, originating from Swedish groups [Benetis et al., Mol. Phys. 48, 329 (1983); Kowalewski et al., Adv. Inorg. Chem. 57, (2005); Larsson et al., J. Chem. Phys. 101, 1116 (1994); T. Nilsson et al., J. Magn. Reson. 154, 269 (2002)] and two different methods based on simulations of the dynamics of electron spin in time domain, developed in Grenoble [Fries and Belorizky, J. Chem. Phys. 126, 204503 (2007); Rast et al., ibid. 115, 7554 (2001)] and Ann Arbor [Abernathy and Sharp, J. Chem. Phys. 106, 9032 (1997); Schaefle and Sharp, ibid. 121, 5387 (2004); Schaefle and Sharp, J. Magn. Reson. 176, 160 (2005)], respectively. In this paper, we report a numerical comparison of the three methods for a large variety of parameter sets, meant to correspond to large and small complexes of gadolinium(III) and of nickel(II). It is found that the agreement between the Swedish and the Grenoble approaches is very good for practically all parameter sets, while the predictions of the Ann Arbor model are similar in a number of the calculations but deviate significantly in others, reflecting in part differences in the treatment of electron spin relaxation. The origins of the discrepancies are discussed briefly.     

Phase behavior and critical properties of size-asymmetric, primitive-model electrolytes

The theory of J. Jiang et al. [J. Chem. Phys. 116, 7977 (2002)] for size-symmetric electrolytes is extended to size-asymmetric electrolytes. When compared to molecular-simulation results, this extension gives the correct trend of critical properties with size asymmetry.     

Recovering the Crooks equation for dynamical systems in the isothermal-isobaric ensemble: a strategy based on the equations of motion

The Crooks equation [Eq. (10) in J. Stat. Phys. 90, 1481 (1998)] relates the work done on a system during a nonequilibrium transformation to the free energy difference between the final and the initial state of the transformation. Recently, the authors have derived the Crooks equation for systems in the canonical ensemble thermostatted by the Nose-Hoover or Nose-Hoover chain method [P. Procacci et al., J. Chem. Phys. 125, 164101 (2006)]. That proof is essentially based on the fluctuation theorem by Evans and Searles [Adv. Phys. 51, 1529 (2002)] and on the equations of motion. Following an analogous approach, the authors derive here the Crooks equation in the context of molecular dynamics simulations of systems in the isothermal-isobaric (NPT) ensemble, whose dynamics is regulated by the Martyna-Tobias-Klein algorithm [J. Chem. Phys. 101, 4177 (1994)]. Their present derivation of the Crooks equation correlates to the demonstration of the Jarzynski identity for NPT systems recently proposed by Cuendet [J. Chem. Phys. 125, 144109 (2006)].     

Quantum Monte Carlo calculations of the dissociation energy of the water dimer

We report diffusion quantum Monte Carlo (DMC) calculations of the equilibrium dissociation energy D(e) of the water dimer. The dissociation energy measured experimentally, D(0), can be estimated from D(e) by adding a correction for vibrational effects. Using the measured dissociation energy and the modern value of the vibrational energy Mas et al., [J. Chem. Phys. 113, 6687 (2000)] leads to D(e)=5.00+/-0.7 kcal mol(-1), although the result Curtiss et al., [J. Chem. Phys. 71, 2703 (1979)] D(e)=5.44+/-0.7 kcal mol(-1), which uses an earlier estimate of the vibrational energy, has been widely quoted. High-level coupled cluster calculations Klopper et al., [Phys. Chem. Chem. Phys. 2, 2227 (2000)] have yielded D(e)=5.02+/-0.05 kcal mol(-1). In an attempt to shed new light on this old problem, we have performed all-electron DMC calculations on the water monomer and dimer using Slater-Jastrow wave functions with both Hartree-Fock approximation (HF) and B3LYP density functional theory single-particle orbitals. We obtain equilibrium dissociation energies for the dimer of 5.02+/-0.18 kcal mol(-1) (HF orbitals) and 5.21+/-0.18 kcal mol(-1) (B3LYP orbitals), in good agreement with the coupled cluster results.     

Diffusion in linear porous media with periodic entropy barriers: A tube formed by contacting spheres

The problem of transport in quasi-one-dimensional periodic structures has been studied recently by several groups [D. Reguera et al., Phys. Rev. Lett.96, 130603 (2006); P. S. Burada et al., Phys. Rev. E75, 051111 (2007); B. Q. Ai and L. G. Liu, ibid.74, 051114 (2006); B. Q. Ai et al., ibid.75, 061126 (2007); B. Q. Ai and L. G. Liu, J. Chem. Phys.126, 204706 (2007); 128, 024706 (2008); E. Yariv and K. D. Dorfman, Phys. Fluids19, 037101 (2007); N. Laachi et al., Europhys. Lett.80, 50009 (2007); A. M. Berezhkovskii et al., J. Chem. Phys.118, 7146 (2003); 119, 6991 (2003)]. Using the concept of "entropy barrier" [R. Zwanzig, J. Phys. Chem.96, 3926 (1992)] one can classify such structures based on the height of the entropy barrier. Structures with high barriers are formed by chambers, which are weakly connected with each other because they are connected by small apertures. To escape from such a chamber a diffusing particle has to climb a high entropy barrier to find an exit that takes a lot of time [I. V. Grigoriev et al., J. Chem. Phys.116, 9574 (2002)]. As a consequence, the particle intrachamber lifetime tau(esc) is much larger than its intrachamber equilibration time, tau(rel), tau(esc)>tau(rel). When the aperture is not small enough, the intrachamber escape and relaxation times are of the same order and the hierarchy fails. This is the case of low entropy barriers. Transport in this case is analyzed in the works of Schmid and co-workers, Liu and co-workers, and Dorfman and co-workers, while the work of Berezhkovskii et al. is devoted to diffusion in the case of high entropy barriers.     

Long-range-corrected hybrids using a range-separated Perdew-Burke-Ernzerhof functional and random phase approximation correlation

We build on methods combining a short-range density functional approximation with a long-range random phase approximation [B. G. Janesko, T. M. Henderson, and G. E. Scuseria, J. Chem. Phys. 130, 081105 (2009)] or second-order screened exchange [J. Paier et al., J. Chem. Phys. 132, 094103 (2010)] by replacing the range-separated local density approximation functional with a range-separated generalized gradient approximation functional in the short range. We present benchmark results that show a marked improvement in the thermodynamic tests over the previous local density approximation-based methods while retaining those methods' excellent performance in van der Waals interactions.     

Sound velocity dispersion in room temperature ionic liquids studied using the transient grating method

Sound velocity is determined by the transient grating method in a range from 10(6) to 10(10) Hz in three room temperature ionic liquids, 1-butyl-3-methylimidazolium bis(trifluoromethanesulfonyl)imide, 1-butyl-3-methylimidazolium hexafluorophosphate, and N,N,N-trimethyl-N-propylammonium bis(trifluoromethanesulfonyl)imide. In all room temperature ionic liquids studied, the sound velocity increased with increasing frequency. The cause of this change is posited to be structural relaxation in the room temperature ionic liquids. Frequency dependence of the sound velocity is not reproduced by a simple Debye relaxation model. The sound velocity dispersion relation in 1-butyl-3-methylimidazolium hexafluorophosphate matches a Cole-Davidson function with parameters determined by a dielectric relaxation [C. Daguenet et al., J. Phys. Chem. B 110, 12682 (2006)], indicating that structural and reorientational relaxations are strongly coupled. Conversely, the sound velocity dispersions of the other two ionic liquids measured do not match those measured for dielectric relaxation, implying that structural relaxation is much faster than the reorientational relaxation. This difference is discussed in relation to the motilities of anions and cations.