Low-temperature heat capacities and thermodynamic properties of rare-earth triisothiocyanate hydrates (Ⅲ)——Sm(NCS)_3·6H_2O,Gd(NCS)_3·6H_2O,Yb(NCS)_3·6H_2O and Y(NCS)_3·6H_2O

The heat capacities of four RE isothiocyanate hydrates,Sm( NCS)3 6H2O,Gd( NCS)3 6H2O,Yb(NCS)3 6H2O and Y( NCS)3 6H2O,have been measured from 13 to 300 K with a fully-automated adiabatic calorimeter No obvious thermal anomaly was observed for the above-mentioned compounds in the experimental tem-peiatnre ranges.The polynomial equations for calculating the heat capacities of the four compounds in the range of 13-300K were obtained by the least-squares fitting based on the experimental Cp data.The Cp values below 13 K were estimated by using the Debye-Einstem heat capacity functions.The standard molar thermodynamic functions were calculated from 0 to 300 K.Gibbs energies of formation were also calculated.     

Isotope effects in ice Ih: a path-integral simulation

Ice Ih has been studied by path-integral molecular dynamics simulations, using the effective q-TIP4P/F potential model for flexible water. This has allowed us to analyze finite-temperature quantum effects in this solid phase from 25 to 300 K at ambient pressure. Among these effects we find a negative thermal expansion of ice at low temperatures, which does not appear in classical molecular dynamics simulations. The compressibility derived from volume fluctuations gives results in line with experimental data. We have analyzed isotope effects in ice Ih by considering normal, heavy, and tritiated water. In particular, we studied the effect of changing the isotopic mass of hydrogen on the kinetic energy and atomic delocalization in the crystal as well as on structural properties such as interatomic distances and molar volume. For D(2)O ice Ih at 100 K we obtained a decrease in molar volume and intramolecular O-H distance of 0.6% and 0.4%, respectively, as compared to H(2)O ice.     

尿嘧啶和5-溴尿嘧啶的低温热容

采用综合物性测量系统(PPMS)的热容测量模块在1.9-300 K温度区间内对两种药物中间体(尿嘧啶和5-溴尿嘧啶)的低温热容进行了测量与研究. 结果表明, 在测量温区内两种化合物的低温热容随温度的上升而逐步增加, 无任何热异常现象产生; 在相同温度下, 5-溴尿嘧啶的热容数值始终高于尿嘧啶. 利用低温热容理论模型对热容数据进行了拟合, 并计算得到了0-300 K温区的摩尔熵变、焓变等热力学函数. 此外, 通过热容拟合数据计算得到的尿嘧啶和5-溴尿嘧啶在298.15 K的标准摩尔规定熵分别为(132.48±1.32)和(165.39±1.65) J·K^-1·mol^-1.     

Thermodynamics of Sodium Uranoniobate and Its Monohydrate

Reaction calorimetry was used to determine standard enthalpies of formation at 298.15 K for crystalline NaNbUO6 (-2580.0±2.0 kJ/mol) and NaNbUO₆·H₂O (-2876.5±1.5 kJ/mol). The heat capacities of these compounds were studied in the range 80-300 K by adiabatic vacuum calorimetry, and their thermodynamic functions were calculated. Standard entropies (-540.5±4.1 and -730.6±4.1 J mol^{- 1} K^{- 1}) and Gibbs functions of formation at 298.15 K (-2419.0±2.0 and -2658.5±2.5 kJ/mol) for NaNbUO₆ and NaNbUO₆. H2O, respectively, were calculated. Thermodynamic functions for a number of reactions yielding these compounds were calculated and examined.     

NMR investigation of methyl-2,4-dimethoxysalicylate: effect of solvent and temperature

(1)H NMR and (13)C NMR of methyl-2,4-dimethoxysalicylate 2 was measured in chloroform-d at the temperature range of 220-330 K, in dimethyl sulfoxide-d(6) at the temperature range of 300-400 K and in a polar protic solvent (CD(3)OD) at 300 K. The structure of 2 in liquid phase (solvent) is compared with those in solid phase (X-ray) and in the gas phase (quantum mechanical calculations). The relationship between molecular geometry, (1)H NMR chemical shift and W coupling of involved protons has a complex nature, but hydrogen bonds [C=O...H-O and C=O...H-CH(2)O] strength is the principle factor.     

Stereographic projection path-integral simulations of (HF)n clusters

We perform several quantum canonical ensemble simulations of (HF)(n) clusters. The HF stretches are rigid, and the stereographic projection path-integral method is employed for the simulation in the resulting curved configuration space. We make use of the reweighted random series techniques to accelerate the convergence of the path-integral simulation with respect to the number of path coefficients. We develop and test estimators for the total energy and heat capacity based on a finite difference approach for non-Euclidean spaces. The quantum effects at temperatures below 400 K are substantial for all sizes. We observe interesting thermodynamic behaviors in the quantum simulations of the octamer and the heptamer.     

液态水热容的量子校正

在不同温度下对液态水进行分子动力学模拟，研究各温度下液态水中各个原子的速度自相关函数密度谱，以考察液态水热容的量子校正随温度的变化规律．研究结果表明，水分子的三个内部振动模式对热容的量子校正不随温度变化，而转动和分子距平衡位置的摆动运动模式的量子校正随温度升高而逐步减小．对于分子动力学模拟结果经温度涨落计算所得的热容进行了量子校正，校正结果与实验值能符合．     

Thermodynamics and molecular dynamics of some ferrocene derivatives

The heat capacities of acetylferrocene, 1,1′-diacetylferrocene, and 1,1′-diethylferrocene were investigated by low-temperature adiabatic calorimetry in the temperature range from 5 to 300 K and their thermodynamic functions were calculated. The enthalpies of combustion of the substances were determined by calorimetry of combustion, and the thermodynamic functions of their formation were calculated by quantum chemistry methods. Inter- and intramolecular interactions of the ferrocene derivatives were also studied by the methods of molecular mechanics and molecular dynamics. Translated fromIzvestiya Akademii Nauk. Seriya Khimicheskaya, No. 9, pp. 1668–1676, September, 1999.     

Oxygen vacancy formation for transient structures on the CeO2(110) surface at 300 and 750 K

Ab initio embedded-cluster calculations have been performed for the CeO2(110) surface using temperature induced structures from molecular dynamics (MD) snapshots. As a first step towards understanding how temperature induced distortions of the surface structure influence the surface oxygen reactivity, the energy cost of removing an O atom from the surface was calculated for 41 snapshots from the MD simulation at 300 K. The quantum mechanical embedded-cluster calculations show that already at 300 K the dynamics causes significant fluctuations (root mean square of 0.37 eV) in the O vacancy formation energy (Evac) while the distribution of the two excess electrons associated with the vacancy is virtually unaffected by the surface dynamics and remains localized on the two Ce ions close to the vacancy. It is also found that the quantum mechanical Evac fluctuations can be reproduced by oxygen vacancy calculations using only the relaxed shell-model force field (FF) itself and the MD geometries. Using the FF as the interaction model, the effect of raising the temperature to 750 K and the effect of doping with Ca were investigated for the oxygen vacancy formation.     

Towards an assessment of the accuracy of density functional theory for first principles simulations of water. II

A series of 20 ps ab initio molecular dynamics simulations of water at ambient density and temperatures ranging from 300 to 450 K are presented. Car-Parrinello (CP) and Born-Oppenheimer (BO) molecular dynamics techniques are compared for systems containing 54 and 64 water molecules. At 300 K, an excellent agreement is found between radial distribution functions (RDFs) obtained with BO and CP dynamics, provided an appropriately small value of the fictitious mass parameter is used in the CP simulation. However, we find that the diffusion coefficients computed from CP dynamics are approximately two times larger than those obtained with BO simulations for T>400 K, where statistically meaningful comparisons can be made. Overall, both BO and CP dynamics at 300 K yield overstructured RDFs and slow diffusion as compared to experiment. In order to understand these discrepancies, the effect of proton quantum motion is investigated with the use of empirical interaction potentials. We find that proton quantum effects may have a larger impact than previously thought on structure and diffusion of the liquid.     

5-氨基间苯二甲酸钠和5-羟基间苯二甲酸钠的热力学性质

在水溶液中合成了5-氨基间苯二甲酸钠(1)和5-羟基间苯二甲酸钠(2)固态样品,元素分析和TG-DTG确定其组成符合C8H5O4NNa2·H2O(1)和C8H4O5Na2·H2O(2).用精密自动绝热热量计测定了它们在78～400K温区的低温热容,将实验值用最小二乘法拟合,得到热容随温度变化的多项式方程,用此方程进行数值积分,得到该温区内每隔5K的舒平热容值和各种热力学函数值.用RD496-2000型微热量计测定了样品在298.15K时的标准摩尔溶解焓分别为ΔsolHmθ(1,s)=-44.552±0.164kJmol-1和θΔsolHm(2,s)=-36.055±0.154kJmol-1,计算了其水合阴离子标准摩尔生成焓分别为θΔfHm(C8H5O4N2-,aq)=-684.56±1.67kJmol-1和ΔfHmθ(C8H4O52-,aq)=-1263.43±2.13kJmol-1.用RBC-II型精密转动弹热量计测定了样品的恒容燃烧热分别为ΔcU(1,s)=-13382.14±5.28Jg-1和ΔcU(2,s)=-10339.15±4.15Jg-1,计算了它们的标准摩尔燃烧焓和标准摩尔生成焓分别为ΔcHmθ(1,s)=-3252.90±1.28kJmol-1和θΔcHm(2,s)=-2522.64±1.01kJmol-1,ΔfHmθ(1,s)=-1406.46±1.66kJmol-1,θΔfHm(2,s)=-1993.79±1.46kJmol-1.     

First principle and ReaxFF molecular dynamics investigations of formaldehyde dissociation on Fe(100) surface

Detailed formaldehyde adsorption and dissociation reactions on Fe(100) surface were studied using first principle calculations and molecular dynamics (MD) simulations, and results were compared with available experimental data. The study includes formaldehyde, formyl radical (HCO), and CO adsorption and dissociation energy calculations on the surface, adsorbate vibrational frequency calculations, density of states analysis of clean and adsorbed surfaces, complete potential energy diagram construction from formaldehyde to atomic carbon (C), hydrogen (H), and oxygen (O), simulation of formaldehyde adsorption and dissociation reaction on the surface using reactive force field, ReaxFF MD, and reaction rate calculations of adsorbates using transition state theory (TST). Formaldehyde and HCO were adsorbed most strongly at the hollow (fourfold) site. Adsorption energies ranged from ?22.9 to ?33.9 kcal/mol for formaldehyde, and from ?44.3 to ?66.3 kcal/mol for HCO, depending on adsorption sites and molecular direction. The dissociation energies were investigated for the dissociation paths: formaldehyde → HCO + H, HCO → H + CO, and CO → C + O, and the calculated energies were 11.0, 4.1, and 26.3 kcal/mol, respectively. ReaxFF MD simulation results were compared with experimental surface analysis using high resolution electron energy loss spectrometry (HREELS) and TST based reaction rates. ReaxFF simulation showed less reactivity than HREELS observation at 310 and 523 K. ReaxFF simulation showed more reactivity than the TST based rate for formaldehyde dissociation and less reactivity than TST based rate for HCO dissociation at 523 K. TST‐based rates are consistent with HREELS observation. © 2013 Wiley Periodicals, Inc.     

Investigating the relationship between infrared spectra of shared protons in different chemical environments: a comparison of protonated diglyme and protonated water dimer

The energetics, dynamics, and infrared spectroscopy of the shared proton in different chemical environments is investigated using molecular dynamics simulations. A three-dimensional potential energy surface (PES) suitable for describing proton transfer between an acceptor and a donor oxygen atom is combined with an all-atom force field to carry out reactive molecular dynamics simulations. The construction of the fully dimensional PES is inspired from the established mixed quantum mechanics/molecular mechanics treatment of larger systems. The "morphing potential" method is used to transform the generic PES for proton transfer along an O...H+...O motif into a three-dimensional PES for proton transfer in protonated diglyme. Using molecular dynamics simulations at finite temperature, the gas phase infrared spectra are calculated for both species from the Fourier transform of the dipole moment autocorrelation function. For protonated diglyme the modes involving the H+ motion are strongly mixed with other degrees of freedom. At low temperature, the O...H+...O asymmetric stretching vibration is found at 870 cm-1, whereas for H5O2+ this band is at 724 cm-1. As expected, the vibrational bands of protonated diglyme show no temperature dependence whereas for H5O2+ at T = 100 K the proton transfer mode is found at 830 cm-1, in good agreement with 861 cm-1 from very recent molecular dynamics simulations.     

Communication: A concerted mechanism between proton transfer of Zundel anion and displacement of counter cation

Ab initio path integral molecular dynamics simulation of M(+)(H(3)O(2)(-)) (M = Li, Na, and K) has been carried out to analyze how the structure and dynamics of a low-barrier hydrogen-bonded Zundel anion, H(3)O(2)(-), can be affected by the counter alkali metal cation, M(+). Our simulation predicts that the quantum proton transfer in Zundel anion can be strongly coupled to the motion of counter cation located nearby. A smaller cation can induce larger structural distortion of the Zundel anion fragment making the proton transfer barrier higher, and hence, lower the vibrational excitation energy. It is also argued that a large H∕D isotope effect is present.     

Comparison of approximate quantum simulation methods applied to normal liquid helium at 4 K

The Feynman-Kleinert linearized path integral molecular dynamics (FK-LPI), ring polymer molecular dynamics (RPMD), and centroid molecular dynamics (CMD) methods are applied to the simulation of normal liquid helium. Comparisons of the simulation results at the T = 4 K and rho = 0.01873 A-3 state point are presented. The calculated quantum correlation functions for the three methods show significant differences, both in the short time and in the intermediate regions of the spectrum. Our simulation results are also compared to the recent results of other approximate quantum simulation methods. We find that FK-LPI qualitatively agrees with other approximate quantum simulation results while CMD and RPMD predict a qualitatively different impulsive rebound in the velocity autocorrelation function. Frequency space analysis reveals that RPMD exhibits a broad high-frequency tail similar to that from quantum mode coupling theory and numerical analytic continuation approaches, while FK-LPI provides a somewhat more rapid decay at high frequency than any of these three methods. CMD manifests a high-frequency component that is greatly reduced compared with the other methods.     

Hydrogen bond dynamics in heavy water studied with quantum dynamical simulations

The structure and dynamics of the hydrogen-bond network in heavy water (D(2)O) is studied as a function of the temperature using quantum dynamical simulations. Our approach combines an ab initio-based representation of the water interactions with an explicit quantum treatment of the molecular motion. A direct connection between the calculated linear and nonlinear vibrational spectra and the underlying molecular dynamics is made, which provides new insights into the rearrangement of the hydrogen-bond network in heavy water. A comparison with previous calculations on liquid H(2)O suggests that tunneling does not effectively contribute to the dynamics of the water hydrogen-bond network above the melting point. However, the effects of nuclear quantization are not negligible at all temperatures and become increasingly important near the melting point, which is in agreement with recent experimental analysis of the structural properties of liquid water as well as of the proton momentum distribution in supercooled water.     

Molecular dynamics simulations of atmospheric oxidants at the air-water interface: solvation and accommodation of OH and O3

A comparative study of OH, O3, and H2O equilibrium aqueous solvation and gas-phase accommodation on liquid water at 300 K is performed using a combination of ab initio calculations and molecular dynamics simulations. Polarizable force fields are developed for the interaction potential of OH and O3 with water. The free energy profiles for transfer of OH and O3 from the gas phase to the bulk liquid exhibit a pronounced minimum at the surface, but no barrier to solvation in the bulk liquid. The calculated surface excess of each oxidant is comparable to calculated and experimental values for short chain, aliphatic alcohols. Driving forces for the surface activity are discussed in terms of the radial distribution functions and dipole orientation distributions for each molecule in the bulk liquid and at the surface. Simulations of OH, O3, and H2O impinging on liquid water with a thermal impact velocity are used to calculate thermal accommodation (S) and mass accommodation (alpha) coefficients. The values of S for OH, O3, and H2O are 0.95, 0.90, and 0.99, respectively. The approaching molecules are accelerated toward the liquid surface when they are approximately 5 angstroms above it. The molecules that reach thermal equilibrium with the surface do so within 2 ps of striking the surface, while those that do not scatter into the gas phase with excess translational kinetic energy in the direction perpendicular to the surface. The time constants for absorption and desorption range from approximately 35 to 140 ps, and the values of alpha for OH, O3, and H2O are 0.83, 0.047, and 0.99, respectively. The results are consistent with previous formulations of gas-phase accommodation from simulations, in which the process occurs by rapid thermal and structural equilibration followed by diffusion on the free energy profile. The implications of these results with respect to atmospheric chemistry are discussed.     

Isotope quantum effects in water around the freezing point

We have measured the difference in electronic structure factors between liquid H(2)O and D(2)O at temperatures of 268 and 273 K with high energy x-ray diffraction. These are compared to our previously published data measured from 279 to 318 K. We find that the total structural isotope effect increases by a factor of 3.5 over the entire range, as the temperature is decreased. Structural isochoric temperature differential and isothermal density differential functions have been used to compare these data to a thermodynamic model based upon a simple offset in the state function. The model works well in describing the magnitude of the structural differences above approximately 310 K, but fails at lower temperatures. The experimental results are discussed in light of several quantum molecular dynamics simulations and are in good qualitative agreement with recent temperature dependent, rotationally quantized rigid molecule simulations.     

Thermodynamics of Lithium Uranoniobate and Its Monohydrate

Standard enthalpies of formation for crystalline LiNbUO₆ (-2619.5±1.5 KJ/mol) and LiNbUO₆· 2H2O (-3251.0±3.0 KJ/mol) at 298.15 K were determined by reaction calorimetry. The heat capacity of these compounds was studied in the range 80-300 K by adiabatic vacuum calorimetry, and their thermodynamic functions were calculated. Standard entropies (-528.5±4.1 and -976.7±4.1 J mol^{- 1} K^{- 1}) and Gibbs functions of formation at 298.15 K (-2462.0±2.5 and -2960.0±4.0 kJ/mol) for LiNbUO₆ and LiNbUO₆·2H₂O, respectively, were calculated. Thermodynamic functions for a number of reactions yielding these compounds were calculated and examined.