Electrostatic calculation of linear and non-linear optical properties of ice Ih,II, IX and VIII

Macroscopic first- and third-order susceptibilities of ice Ih, ice II, ice IX and ice VIII are calculated using static and frequency-dependent electronic and static vibrational molecular (hyper)polarizabilities at the MP2 level. The molecular properties are in good agreement with experiment and with high-level ab initio calculations. Intermolecular electrostatic and polarization effects due to induced dipoles are taken into account using a rigorous local-field theory. The electric field due to permanent dipoles is used to calculate effective in-crystal (hyper)polarizabilities. The polarizability depends only weakly on the permanent field, but the dipole moment and the hyperpolarizabilities are strongly affected. The calculated linear susceptibility is in good agreement with available experimental data for ice Ih, and the third-order susceptibility for a third harmonic generation experiment is in reasonable agreement with experimental values for liquid water. The molecular vibrational contributions have a small effect on the susceptibilities. The electric properties of a water tetramer are calculated and used to estimate the effect of non-dipolar interactions on the susceptibilities of ice Ih, which are found to be small.     

Molecular-dynamics studies of surface of ice Ih

We performed molecular dynamics calculations of surface of ice Ih in order to investigate formation mechanism of melting layer on the surface. The results showed that the vibrational amplitude of the atoms in the surface layer greatly depends on the crystal orientation, whereas that in the ice bulk is isotropic. The anisotropy of the vibration is due to a dangling motion of the free O-H bonds exist at the surface layer. The dangling motion enhances the rotational motion of the water molecules. The vibrational density of state showed a coupling between the rotational vibration and the lattice vibration of the water molecules in the surface layer. The coupling of the vibrations causes a distortion of ice lattice. Through the hydrogen-bonding network, the distortion transmits to the interior of the crystal. We conclude that the dangling motion of the free O-H bonds exist at the surface layer is one of the dominant factors governing the surface melting of ice crystal.     

Two-dimensional infrared spectroscopy of neat ice Ih

The OH stretch line shape of ice Ih exhibits distinct peaks, the assignment of which remains controversial. We address this longstanding question using two dimensional infrared (2D IR) spectroscopy of the OH stretch of H(2)O and the OD stretch of D(2)O of ice Ih at T = 80 K. The isotropic response is dominated by a 2D line shape component which does not depend on the pump pulse frequency. The decay time of the component that does depend on the pump frequency is calculated using singular value decomposition (bi-exponential decay H(2)O: 30 fs, 490 fs; D(2)O: 40 fs, 690 fs). The anisotropic contribution exhibits on-diagonal peaks, which decay on a very fast timescale (H(2)O: 85 fs; D(2)O: 65 fs), with no corresponding anisotropic cross-peaks. Both isotropic and anisotropic results indicate that randomization of excited dipoles occurs with a very rapid rate, just like in neat liquid water. We conclude that the underlying mechanism relates to the complex interplay between exciton migration and exciton-phonon coupling.     

Two-dimensional infrared spectroscopy of isotope-diluted ice Ih

We present experimental 2D IR spectra of isotope diluted ice Ih (i.e., the OH stretch mode of HOD in D(2)O and the OD stretch mode of HOD in H(2)O) at T = 80 K. The main spectral features are the extremely broad 1-2 excited state transition, much broader than the corresponding 0-1 groundstate transition, as well as the presence of quantum beats. We do not observe any inhomogeneous broadening that might be expected due to proton disorder in ice Ih. Complementary, we perform simulations in the framework of the Lippincott-Schroeder model, which qualitatively reproduce the experimental observations. We conclude that the origin of the observed line shape features is the coupling of the OH-vibrational coordinate with crystal phonons and explain the beatings as a coherent oscillation of the O···O hydrogen bond degree of freedom.     

Transient absorption of vibrationally excited ice Ih

The ultrafast dynamics of HDO:D2O ice Ih at 180 K is studied by midinfrared ultrafast pump-probe spectroscopy. The vibrational relaxation of HDO:D2O ice is observed to proceed via an intermediate state, which has a blueshifted absorption spectrum. Polarization resolved measurements reveal that the intermediate state is part of the intramolecular relaxation pathway of the HDO molecule. In addition, slow dynamics on a time scale of the order of 10-100 ps is observed, related to thermally induced collective reorganizations of the ice lattice. The transient absorption line shape is analyzed within a Lippincott-Schroeder model for the OH-stretch potential. This analysis identifies the main mechanism behind the strong spectral broadening of the v(OH)=1-->2 transition.     

The structure,thermodynamic properties and infrared spectra of liquid water and ice

Monte Carlo simulations are used, together with models of the intramolecular and intermolecular potential surfaces, to model liquid water and several phases of ice. Intramolecular relaxation makes important contributions to both thermodynamic and structural properties. A quantum local mode analysis of the Monte Carlo configurations is used to predict the density of states and infrared absorption intensities for the intramolecular bending and stretching vibrations. The large shifts from the gas phase OH stretch frequencies observed experimentally in the liquid and solid phases are due to anharmonic terms in the intramolecular surface rather than to harmonic intermolecular coupling. A significant contribution to observed changes in IR intensity on condensation arises from the large molecular polarisability.     

Simulations of proton order and disorder in ice Ih

Computer simulations of ice Ih with different proton orientations are presented. Simulations of proton disordered ice are carried out using a Monte Carlo method which samples over proton degree of freedom, allowing for the calculation of the dielectric constant and for the examination of the degree of proton disorder. Simulations are also presented for two proton ordered structures of ice Ih, the ferroelectric Cmc2(1) structure or ice XI and the antiferroelectric Pna2(1) structure. These simulations indicate that a transition to a proton ordered phase occurs at low temperatures (below 80 K). The symmetry of the ordered phase is found to be dependent on the water potential. The stability of the two proton ordered structures is due to a balance of short-ranged interactions which tend to stabilize the Pna2(1) structure and longer-range interactions which stabilize the Cmc2(1) structure.     

A van der Waals density functional study of ice Ih

Density functional theory with the van der Waals density functional (vdW-DF) is used to calculate equilibrium crystal structure, binding energy, and bulk modulus of ice Ih. It is found that although it overestimates the equilibrium volume, vdW-DF predicts accurate binding energy of ice Ih, as compared with high level quantum chemistry calculations and experiment. Inclusion of the nonlocal correlation, i.e., van der Waals interaction, leads to an overall improvement over the standard generalized gradient approximation in describing water ice.     

Melting temperature of ice Ih calculated from coexisting solid-liquid phases

We carried out molecular-dynamics simulations by using the two-phase coexistence method with the constant pressure, particle number, and enthalpy ensemble to compute the melting temperature of proton-disordered hexagonal ice I(h) at 1-bar pressure. Four models of water were considered, including the widely used TIP4P [W. L. Jorgensen, J. Chandrasekha, J. D. Madura, R. W. Impey, and M. L. Klein, J. Chem. Phys.79, 926 (1983)] and TIP5P [M. W. Mahoney and W. L. Jorgensen J. Chem. Phys.112, 8910 (2000)] models, as well as recently improved TIP4P and TIP5P models for use with Ewald techniques-the TIP4P-Ew [W. Horn, W. C. Swope, J. W. Pitera, J. C. Madura, T. J. Dick, G. L. Hura, and T. Head-Gordon, J. Chem. Phys.120, 9665 (2004)] and TIP5P-Ew [S. W. Rick, J. Chem. Phys.120, 6085 (2004)] models. The calculated melting temperature at 1 bar is T(m) = 229 +/- 1 K for the TIP4P and T(m) = 272.0 +/- 0.6 K for the TIP5P ice I(h), both are consistent with previous simulations based on free-energy methods. For the TIP4P-Ew and TIP5P-Ew models, the calculated melting temperature is T(m) = 257.0 +/- 1.1 K and T(m) = 253.9 +/- 1.1 K, respectively.     

Infrared and Raman line shapes for ice Ih. II. H2O and D2O

We present a theoretical study of infrared and Raman line shapes of polycrystalline and single crystal ice Ih, for both water and heavy water, at 1, 125, and 245 K. Our calculations involve a mixed quantum/classical approach, a new water simulation model with explicit three-body interactions, transition frequency and dipole maps, and intramolecular and intermolecular vibrational coupling maps. Our theoretical spectra are in reasonable agreement with experimental spectra (available only near the two higher temperatures). We trace the origins of the different spectral peaks to weak and strong intermolecular couplings. We also discuss the delocalization of the vibrational eigenstates in terms of the competing effects of disorder and coupling.     

Identification,composition, thermodynamic and structural properties of a pyroaurite-like iron(II)-iron (III) hydroxy-oxalate green rust

The aerial oxidation of aqueous suspensions of ferrous hydroxide precipitated from ferrous oxalate and caustic soda can lead to an iron (II)-iron (III) hydroxy-oxalate of the pyroaurite group, a GR(C₂O₄^2-) Green Rust. As other GR compounds, it is unstable with respect to the action of oxygen and oxidises later on. Its chemical composition was determined to be [Fe^II₆ Fe^III₂(OH)₁₆]²⁺[C₂O₄^2- · nH₂O], with n more likely equal to 3 on the basis of structural considerations. The composition does not vary and the Fe (II) / Fe (III) ratio in the compound is measured by means of transmission Mössbauer spectroscopy at 78 K and 20 K in the range from 2.8 to 3.2 for various samples at various stages of the reaction. GR(C₂O₄^2-) is paramagnetic at both temperatures and is unambiguously distinguished from ferrous hydroxide, the initial reactant, and magnetite, the main final product, which are magnetically ordered at 20 K. The spectrum of the GR compound is composed of three quadrupole doublets, one due to the Fe(III) cations characterised by a small quadrupole splitting ΔE_Q of 0.40 mm s⁻¹, and two due to the Fe(II) cations, characterised by larger ΔE_Q values of about 2.55 and 2.85 mm s⁻¹. Finally, from the observed equilibrium conditions between ferrous hydroxide and GR(C₂O₄^2-), the standard free enthalpy of formation of GR(C₂O₄^2-) was computed to be : ΔG°_f[Fe^II₆ Fe^III₂ (OH)₁₆]²⁺[C₂O₄^2- · 3H₂O] = −5383 ± 3 kJ mol⁻¹.     

Electrochemical and thermodynamic properties of hexacyanoferrate(II)/(III) redox system on multi-walled carbon nanotubes

Novel films consist of multi-walled carbon nanotubes (MWCNT) were fabricated by means of catalytic chemical vapor deposition (CVD) technique with decomposition of either acetonitrile (ACN) or benzene (BZ) using ferrocene (FeCp₂) as catalyst. The electrochemical and thermodynamic behavior of the ferrocyanide/ferricyanide, [Fe(CN)₆]^3−/4− redox couple on synthesized MWCNT-based films was investigated by means of cyclic voltammetry (CV) and electrochemical impedance spectroscopy (EIS) techniques at T = (278.15, 283.15, 293.15, and 303.15) K. The redox couple [Fe(CN)₆]^3−/4− behaves quasi-reversibly on fabricated MWCNT-based films and its reversibility is enhanced upon increasing temperature. Namely, the findings establish that with the rise in temperature the barrier for interfacial electron transfer decreases, leading, consequently, to an enhancement of the kinetics of the charge transfer process. According to thermodynamics the equilibrium of the redox process is shifted towards the formation of [Fe(CN)₆]³⁻ at elevated temperatures.     

Direct dynamics study of ultrafast vibrational energy relaxation in ice Ih

Car-Parrinello molecular dynamics (CPMD) and a previously developed wave packet model are used to study ultrafast relaxation in water clusters. Water clusters of 15 water molecules are used to represent ice Ih. The relaxation is studied by exciting a symmetric or an asymmetric stretch mode of the central water molecule. The CPMD results suggest that relaxation occurs within 100 fs. This is in agreement with experimental work by Woutersen and Bakker and the earlier wave packet calculations. The CPMD results further indicate that the excitation energy is transferred both intramolecularly and intermolecularly on roughly the same time scale. The intramolecular energy transfer occurs predominantly between the symmetric and asymmetric modes while the bend mode is largely left unexcited on the short time scale studied here.     

Some thermodynamic properties of carboxylato-pentammine cobalt(III) ions

Measurements are reported on the thermodynamic properties for the acid ionisation of various carboxylato-pentammine cobalt(III) ions, where the organic ligand is a dicarboxylate ion ^?OCORCOOH. It is shown that the results can be correlated with the ionisation of the free acids, on the basis of a simple electrostatic model. Results are also reported on the enthalpies of hydrolysis of the same ions, obtained indirectly from the enthalpies of reaction with sodium sulphide. The results for the complexes of the dicarboxylic acids which were studied (maleato, phthalato, fumarato, succinato and malonato) are not very different, while the acetato complex is somewhat more unstable. There is some indication that the complexes which hydrolyse faster are also thermodynamically less stable.     

Modeling equilibrium concentrations of Bjerrum and molecular point defects and their complexes in ice Ih

We present a model for the determination of the thermal equilibrium concentrations of Bjerrum defects, molecular point defects, and their aggregates in ice I(h). First, using a procedure which minimizes the free energy of an ice crystal with respect to the numbers of defect species, we derive a set of equations for the equilibrium concentrations of free Bjerrum and point defects, as well their complexes. Using density-functional-theory calculations, we then evaluate the binding energies of Bjerrum-defect/vacancy and Bjerrum-defect/interstitial complexes. In contrast to the complexes which involve the molecular vacancy, the results suggest that the molecular interstitial binds preferentially to the D-type Bjerrum defect. Using both theoretical binding and formation free energies as well as the available experimental data, we find that the preferential binding and the substantial presence of the interstitial as the predominant point defect in ice I(h) may lead to conditions in which the number of free D defects becomes considerably smaller than that of free L defects. Such a scenario could possibly be involved in the experimentally observed inactivity of D-type Bjerrum defects in the electrical properties of ice I(h).     

Crystal structure and thermodynamic properties of dipotassium diiron(III) hexatitanium oxide

In the present work, the temperature dependence of heat capacity of dipotassium diiron(III) hexatitanium oxide has been measured for the first time in the range from 10 to 300 K by means of precision adiabatic vacuum calorimetry. The experimental data were used to calculate standard thermodynamic functions, namely the heat capacity $ C_{p}^{ \circ } (T) $ , enthalpy $ H^{ \circ } (T) - H^{ \circ } (0) $ , entropy $ S^{ \circ } (T) - S^{ \circ } (0), $ and Gibbs function $ G^{ \circ } (T) - H^{ \circ } (0) $ for the range from T → 0 to 300 K. The structure of K₂Fe₂Ti₆O₁₆ is refined by the Rietveld method: space group I4/m, Z = 1, a = 10.1344(2) Å, c = 2.97567(4) Å, V = 305.618(7) Å³. The high-temperature X-ray diffraction was used for the determination of coefficients of thermal expansion.     

Investigation of the hydrogen bonding in ice Ih by first-principles density function methods

It is a well recognized difficult task to simulate the vibrational dynamics of ices using the density functional theory (DFT), and there has thus been rather limited success in modelling the inelastic neutron scattering (INS) spectra for even the simplest structure of ice, ice Ih, particularly in the translational region below 400 cm(-1). The reason is partly due to the complex nature of hydrogen bonding (H-bond) among water-water molecules which require considerable improvement of the quantum mechanical simulation methods, and partly owing to the randomness of protons in ice structures which often requires simulation of large super-lattices. In this report, we present the first series of successful simulation results for ice Ih using DFT methods. On the basis of the recent advancement in the DFT programs, we have achieved for the first time theoretical outcomes that not only reproduce the rotational frequencies between 500 to 1200 cm(-1) for ice Ih, but also the two optic peaks at ～240 and 320 cm(-1) in the translational region of the INS spectra [J. C. Li, J. Chem. Phys 105, 6733 (1996)]. Besides, we have also investigated the impact of pairwise configurations of H(2)O molecules on the H-bond and found that different proton arrangements of pairwise H(2)O in the ice Ih crystal lattice could not alter the nature of H-bond as significantly as suggested in an early paper [J. C. Li and D. K. Ross, Nature (London) 365, 327 (1993)], i.e., reproducing the two experimental optic peaks do not need to invoke the two H-bonds as proposed in the previous model which led to considerable debates. The results of this work suggest that the observed optic peaks may be attributed to the coupling between the two bands of H-O stretching modes in H(2)O. The current computational work is expected to shed new light on the nature of the H-bonds in water, and in addition to offer a new approach towards probing the interaction between water and biomaterials for which H-bond is essential.     

Quantum path integral simulation of isotope effects in the melting temperature of ice Ih

The isotope effect in the melting temperature of ice Ih has been studied by free energy calculations within the path integral formulation of statistical mechanics. Free energy differences between isotopes are related to the dependence of their kinetic energy on the isotope mass. The water simulations were performed by using the q-TIP4P/F model, a point charge empirical potential that includes molecular flexibility and anharmonicity in the OH stretch of the water molecule. The reported melting temperature at ambient pressure of this model (T=251?K) increases by 6.5±0.5 and 8.2±0.5?K upon isotopic substitution of hydrogen by deuterium and tritium, respectively. These temperature shifts are larger than the experimental ones (3.8 and 4.5 K, respectively). In the classical limit, the melting temperature is nearly the same as that for tritiated ice. This unexpected behavior is rationalized by the coupling between intermolecular interactions and molecular flexibility. This coupling makes the kinetic energy of the OH stretching modes larger in the liquid than in the solid phase. However, the opposite behavior is found for intramolecular modes, which display larger kinetic energy in ice than in liquid water.     

Ruthenium(II) thiol and H2S complexes: synthesis, characterization, and thermodynamic properties

The known, green, five-coordinate species trans-RuCl(2)(P-N)(PPh(3)) react with R'SH thiols to give yellow cis-RuCl(2)(P-N)(PPh(3))(R'SH) products (P-N = o-diphenylphosphino-N,N'-dimethylaniline; R' = alkyl). The MeSH and EtSH compounds are structurally characterized, with the former being the first reported for a transition metal-MeSH complex, while the thiol complexes with R' = (n)Pr, (i)Pr, (n)Pn (pentyl), (n)Hx (hexyl), and Bn (benzyl) are synthesized in situ. Other trans-RuX(2)(P-N)(PR(3)) complexes (X = Br, I; R = Ph, p-tolyl) are synthesized, and their H(2)S adducts, of a type reported earlier by our group, are also prepared. Thermodynamic data are presented for the reversible formation of the MeSH and EtSH complexes and the H(2)S analogues. The Ru(II)Cl(2)(P-N)(PPh(3)) complex in solution decomposes under O(2) to form [Ru(III)Cl(P-N)](2)(μ-O)(μ-Cl)(2).