Search

The theory on electric conductivity of electrolyte solutions we have developed [T. Yamaguchi, T. Matsuoka, and S. Koda, J. Chem. Phys. 127, 064508 (2007)] is applied to a model electrolyte solution that shows a minimum of equivalent conductivity as the function of concentration [T. Yamaguchi, T. Akatsuka, and S. Koda, J. Chem. Phys. 134, 244506 (2011)]. The theory succeeds in reproducing the equivalent conductivity minimum, whereas the mode-coupling theory (MCT) underestimates the conductivity in the low-concentration regime. The theory can also reproduce the decrease in the relaxation time of conductivity with increasing the concentration we have demonstrated with a Brownian dynamics simulation. A detailed analysis shows that the relaxation of the conductivity occurs through two processes. The faster one corresponds to the collision between a cation and an anion, and the slower one does to the polarization of the ionic atmosphere. The increase in the equivalent conductivity with concentration is attributed to the decrease in the effect of the ionic atmosphere, which is in turn explained by the fact that the counter ion cannot penetrate into the repulsive core when the Debye screening length is compatible or smaller than the ionic diameter. The same mechanism is also observed in MCT calculation with static structure factor determined by mean-spherical approximation. [ABSTRACT FROM AUTHOR]

The theory on the frequency-dependent electric conductivity of electrolyte solutions proposed previously by Yamaguchi et al. [J. Chem. Phys. 127, 234501 (2007)] is extended to include the hydrodynamic interaction between ions. The theory is applied to the aqueous solution of NaCl and the concentration dependence of the conductivity agrees well with that determined by experiments. The effects of the hydrodynamic and relaxation effects are highly nonadditive in the concentrated solution, because the hydrodynamic interaction between ions affects the time-dependent response of the ionic atmosphere. The decrease in the electric conductivity is divided into the contributions of ion pair distribution at various distances. The long-range ionic atmosphere plays a major role at the concentration as low as 0.01 mol/kg, whereas the contribution of the contact ion pair region is important at 1 mol/kg. The magnitude of the contribution of the contact ion pair region is scarcely dependent on the presence of the hydrodynamic interaction. The transport number of cation is calculated to be a decreasing function of concentration as is observed in experiments. [ABSTRACT FROM AUTHOR]

The theory on the ultrasonic absorption of electrolyte solutions we have proposed previously [T. Yamaguchi et al., J. Chem. Phys. 126, 144505 (2007)] is extended to calculate the frequency-dependent electric conductivity of the solution. The ionic contribution of the dielectric relaxation spectrum is obtained at the same time. The theory is able to handle the contributions of both the ion-pair dynamics and the relaxation of ionic atmosphere, as is the case of ultrasonic absorption. The effect of the barrier height between the contact and solvent-separated ion pairs is investigated in detail. It is clarified that the competition between the dissociation and reorientational relaxation rates of the contact ion pair is an important factor for the ion pair to be regarded as the ion pair in terms of ionic conductivity. [ABSTRACT FROM AUTHOR]

The theory on the ultrasonic absorption spectrum of electrolytic solutions recently proposed by us is applied to the model system that resembles to the aqueous solution of MgSO4. The charges on ions are reduced to ±1.5e in order to obtain the equilibrium structure by the integral equation theory. The theory reproduces the existence of two relaxations around 100 kHz and 1 GHz. The physical origin of the relaxation is analyzed based on the theoretical expression. The slower relaxation is shown to originate in the formation of contact ion pair, in harmony with the conventional assignment. The amplitude of this relaxation agrees with the experimental one fairly well. The absorption cross section is a weakly increasing function of the concentration of the salt in theory, whereas it depends little on the concentration in experiment, which is ascribed to the weaker association of the pair in the theory. The deviation from the Debye relaxation is found for the faster process, and the concentration dependence is small. The analysis shows that this relaxation stems from the coupling between the pressure and the long-range concentration fluctuation, and the concentration independence and the non-Debye relaxation are explained based on the theoretical analysis. In particular, the theory demonstrates that this process has the t-3/2 tail in the time domain, which is confirmed by numerical calculation. The deviation of the theoretical relaxation amplitude from the experimental one is elucidated in terms of the theoretical expression of the coefficient. [ABSTRACT FROM AUTHOR]

A theory is formulated that describes the sound absorption of electrolytic solutions due to the relative motion of ions, including the formation of ion pairs. The theory is based on the Kubo-Green formula for the bulk viscosity. The time correlation function of the pressure is projected onto the bilinear product of the density modes of ions. The time development of the product of density modes is described by the diffusive limit of the generalized Langevin equation, and approximate expressions for the three- and four-body correlation functions required are given with the hypernetted-chain integral equation theory. Calculations on the aqueous solutions of model electrolytes are performed. It is demonstrated that the theory describes both the activated barrier crossing between contact and solvent-separated ion pairs and the Coulombic correlation between ions. [ABSTRACT FROM AUTHOR]

A theory on the time development of the density and current fields of simple fluids under an external field is formulated through the generalized Langevin formalism. The theory is applied to the linear solvation dynamics of a fixed solute regarding the solute as the external field on the solvent. The solute-solvent-solvent three-body correlation function is taken into account through the hypernetted-chain integral equation theory, and the time correlation function of the random force is approximated by that in the absence of the solute. The theoretical results are compared with those of molecular-dynamics (MD) simulation and the surrogate theory. As for the transient response of the density field, our theory is shown to be free from the artifact of the surrogate theory that the solvent can penetrate into the repulsive core of the solute during the relaxation. We have also found a large quantitative improvement of the solvation correlation function compared with the surrogate theory. In particular, the short-time part of the solvation correlation function is in almost perfect agreement with that from the MD simulation, reflecting that the short-time expansion of the theoretical solvation correlation function is exact up to t2 with the exact three-body correlation function. A quantitative improvement is found in the long-time region, too. Our theory is also applied to the force-force time correlation function of a fixed solute, and similar improvement is obtained, which suggests that our present theory can be a basis to improve the mode-coupling theory on the solute diffusion. [ABSTRACT FROM AUTHOR]

The transient response function of the density profile of the solvent around a solute during the translational diffusion of the solute is formulated based on the generalized Langevin formalism. The resultant theory is applied to both neat Lennard-Jones fluids and cations in liquid water, and the response functions are obtained from the analysis of the molecular dynamics simulations. In the case of the self-diffusion of Lennard-Jones fluids, the responses of the solvation structures are in harmony with conventional pictures based on the mode-coupling theory, that is, the binary collision in the low-density fluids, the backflow effect from medium to high density fluids, and the backscatter effect in the liquids near the triple point. In the case of cations in water, the qualitative behavior is strongly dependent on the size of cations. The pictures similar to simple dense liquids are obtained for the large ion and the neutral molecule, while the solvent waters within the first solvation shell of small ions show an oscillatory response in the short-time region. In particular, the oscillation is remarkably underdumped for lithium ion. The origin of the oscillation is discussed in relation to the theoretical treatment of the translational diffusion of ions in water. © 2005 American Institute of Physics. [ABSTRACT FROM AUTHOR]

The molecular motion of water in water–hydrophobic solute mixtures was investigated by the mode-coupling theory for molecular liquids based on the interaction-site description. When the model Lennard-Jones solute was mixed with water, both the translational and reorientational motions of solvent water become slower, in harmony with various experiments and molecular dynamics simulations. We compared the mechanism of the slowing down with that of the pressure dependence of the molecular motion of neat water [T. Yamaguchi, S.-H. Chong, and F. Hirata, J. Chem. Phys. 119, 1021 (2003)]. We found that the decrease in the solvent mobility caused by the solute can essentially be elucidated by the same mechanism: That is, the fluctuation of the number density of solvent due to the cavity formation by the solute strengthens the friction on the collective polarization through the dielectric friction mechanism: We also employed the solute molecule that is the same as solvent water except for the amount of partial charges, in order to alter the strength of the solute–solvent interaction continuously. The mobility of the solvent water was reduced both by the hydrophobic and strongly hydrophilic solutes, but it was enhanced in the intermediate case. Such a behavior was discussed in connection with the concept of positive and negative hydrations. © 2004 American Institute of Physics. [ABSTRACT FROM AUTHOR]

The transport coefficients of three ionic liquids, lithium bromide (LiBr), rubidium bromide (RbBr), and molten silica (SiO[sub 2]) are calculated by the mixture mode-coupling theory. The static partial structure factors required are obtained from the interionic interaction potential by the Ornstein–Zernike/hypernetted-chain integral equation. The anomalous pressure dependence of the transport properties, the increase in the molar ionic conductivity of LiBr and the fluidity of SiO[sub 2] in the low-pressure region, is reproduced qualitatively well by our theoretical calculation. The calculated results are analyzed in the similar way as that for water performed by Yamaguchi et al. [T. Yamaguchi, S.-H. Chong, and F. Hirata, J. Chem. Phys. 119, 1021 (2003)], and we found that the friction on the electric current caused by the coupling between the charge- and number-density modes is effective to the increase of the transport coefficient with pressure, as is the case of water. Comparing the results for LiBr and RbBr, the contribution of the electrostatic friction is smaller for RbBr, which leads to the normal pressure dependence of its molar ionic conductivity. The negative values of the Nernst–Einstein deviation parameter for the ionic conductivity of LiBr and SiO[sub 2] reported by previous MD simulations are also explained consistently. Furthermore, it is shown that the mechanism for the anomalous pressure dependence of the fluidity of molten SiO[sub 2] demonstrated in this work is consistent with a conventional picture that the five-coordinated silicon atom is important to enhance the ionic mobility. © 2003 American Institute of Physics. [ABSTRACT FROM AUTHOR]

The ultrasonic vibration potential (UVP) and the electrokinetic sonic amplitude (ESA) are described in terms of equilibrium time-correlation functions using the linear response theory. The reciprocal relationship between UVP and ESA is shown based on the formulation. By introducing the generalized Langevin theory and taking the hydrodynamic limit, it is discussed how the effective volume of ions in the UVP measurement is related to their partial molar thermodynamic quantities. The effective volume is proven exactly equal to the isothermal partial molar volume in the isothermal formulation. The effect of the adiabaticity of sound wave is also investigated. © 2003 American Institute of Physics. [ABSTRACT FROM AUTHOR]

Equilibrium and non-equilibrium molecular dynamics simulations are performed on liquid hexane in order to clarify the origin of the Cox-Merz rule of liquids composed of chain-like molecules. The relation between the frequency-dependent complex shear viscosity and the shear-rate dependent nonlinear viscosity follows the Cox-Merz rule as expected. The slowest viscoelastic relaxation mode is explained by the translation-orientation coupling mechanism, and the saturation of the shear-induced orientational order is observed in the non-equilibrium simulation at the onset of the shear thinning. The origin of the Cox-Merz rule is discussed in terms of the translation-orientation coupling.

A theoretical method for real-time dynamics of nonadiabatic reorganization of electronic configurations in molecules is developed, with dual aim that the intramolecular electron dynamics can be probed by means of direct and/or indirect photoionizations and that the physical origins behind photoionization signals attained in the time domain can be identified in terms of the language of time-dependent quantum chemistry. In doing so, we first formulate and implement a new computational scheme for nonadiabatic electron dynamics associated with molecular ionization, which well fits in the general theory of nonadiabatic electron dynamics. In this method, the total nonadiabatic electron wavepackets are propagated in time directly with complex natural orbitals without referring to Hartree-Fock molecular orbitals, and the amount of electron flux from a molecular region leading to ionization is evaluated in terms of the relevant complex natural orbitals. In the second half of this paper, we apply the method to electron dynamics in the elementary processes consisting of the Auger decay to demonstrate the methodological significance. An illustrative example is taken from an Auger decay starting from the 2a 1 orbital hole-state of H 2 O + . The roles of nuclear momentum (kinetic) couplings in electronic-state mixing during the decay process are analyzed in terms of complex natural orbitals, which are schematically represented in the conventional language of molecular symmetry of the Hartree-Fock orbitals.

A theory for dynamics of molecular photoionization from nonadiabatic electron wavepackets driven by intense pulse lasers is proposed. Time evolution of photoelectron distribution is evaluated in terms of out-going electron flux (current of the probability density of electrons) that has kinetic energy high enough to recede from the molecular system. The relevant electron flux is in turn evaluated with the complex-valued electronic wavefunctions that are time evolved in nonadiabatic electron wavepacket dynamics in laser fields. To uniquely rebuild such wavefunctions with its electronic population being lost by ionization, we adopt the complex-valued natural orbitals emerging from the electron density as building blocks of the total wavefunction. The method has been implemented into a quantum chemistry code, which is based on configuration state mixing for polyatomic molecules. Some of the practical aspects needed for its application will be presented. As a first illustrative example, we show the results of hydrogen molecule and its isotope substitutes (HD and DD), which are photoionized by a two-cycle pulse laser. Photon emission spectrum associated with above threshold ionization is also shown. Another example is taken from photoionization dynamics from an excited state of a water molecule. Qualitatively significant effects of nonadiabatic interaction on the photoelectron spectrum are demonstrated.

The NaCl salt-solution interface often serves as an example of an uncharged surface. However, recent laser-Doppler electrophoresis has shown some evidence that the NaCl crystal is positively charged in its saturated solution. Using molecular dynamics (MD) simulations, we have investigated the NaCl salt-solution interface system, and calculated the solubility of the salt using the direct method and free energy calculations, which are kinetic and thermodynamic approaches, respectively. The direct method calculation uses a salt-solution combined system. When the system is equilibrated, the concentration in the solution area is the solubility. In the free energy calculation, we separately calculate the chemical potential of NaCl in two systems, the solid and the solution, using thermodynamic integration with MD simulations. When the chemical potential of NaCl in the solution phase is equal to the chemical potential of the solid phase, the concentration of the solution system is the solubility. The advantage of using two different methods is that the computational methods can be mutually verified. We found that a relatively good estimate of the solubility of the system can be obtained through comparison of the two methods. Furthermore, we found using microsecond time-scale MD simulations that the positively charged NaCl surface was induced by a combination of a sodium-rich surface and the orientation of the interfacial water molecules.

Low-temperature and high-pressure experiments were performed with filled ice Ih structure of methane hydrate under 2.0-77.0 GPa and 30-300 K using diamond anvil cells and a helium-refrigeration cryostat. In situ X-ray diffractometry revealed distinct changes in the compressibility of the axial ratios of the host framework with pressure. Raman spectroscopy showed a split in the C-H vibration modes of the guest methane molecules, which was previously explained by the orientational ordering of the guest molecules. The pressure and temperature conditions at the split of the vibration modes agreed well with those of the compressibility change. The results indicate the following: (i) the orientational ordering of the guest methane molecules from an orientationally disordered state occurred at high pressures and low temperatures; and (ii) this guest ordering led to anisotropic contraction in the host framework. Such guest orientational ordering and subsequent anisotropic contraction of the host framework were similar to that reported previously for filled ice Ic hydrogen hydrate. Since phases with different guest-ordering manners were regarded as different phases, existing regions of the guest disordered-phase and the guest ordered-phase were roughly estimated by the X-ray study. In addition, above the pressure of the guest-ordered phase, another high-pressure phase developed in the low-temperature region. The deuterated-water host samples were also examined, and the influence of isotopic effects on guest ordering and phase transformation was observed.

Low-temperature and high-pressure experiments were performed on the filled ice Ic structure of hydrogen hydrate at previously unexplored conditions of 5-50 GPa and 30-300 K using diamond anvil cells and a helium-refrigeration cryostat. In situ x-ray diffractometry revealed that the cubic filled ice Ic structure transformed to tetragonal at low temperatures and high pressures; the axis ratio of the tetragonal phase changed depending on the pressure and temperature. These results were consistent with theoretical predictions performed via first principle calculations. The tetragonal phase was determined to be stable above 20 GPa at 300 K, above 15 GPa at 200 K, and above 10 GPa at 100 K. Further changes in the lattice parameters were observed from about 45-50 GPa throughout the temperature region examined, which suggests the transformation to another high-pressure phase above 50 GPa. In our previous x-ray study that was performed up to 80 GPa at room temperature, a similar transformation was observed above 50 GPa. In this study, the observed change in the lattice parameters corresponds to the beginning of that transformation. The reasons for the transformation to the tetragonal structure are briefly discussed: the tetragonal structure might be induced due to changes in the vibrational or rotational modes of the hydrogen molecules under low temperature and high pressure.

The interfacial structure of heptane and toluene at oil-silica interfaces has previously been studied by sum frequency generation [Z. Yang et al., J. Phys. Chem. C. 113, 20355 (2009)]. It was found that the toluene molecule is almost perpendicular to the silica surface with a tilt angle of about 25°. Here, we have investigated the structural properties of toluene and heptane at oil-silica interfaces using molecular dynamics simulations for two different surfaces: the oxygen-bridging (hydrophobic) and hydroxyl-terminated (hydrophilic) surfaces of quartz (silica). Based on the density profile, it was found that both heptane and toluene oscillate on silica surfaces, with heptane showing more oscillation peaks. Furthermore, the toluene molecules of the first layer were found to have an asymmetric distribution of orientations, with more CH(3) groups pointed away from the silica surface than towards the silica surface. These findings are generally consistent with previous experiments, and reveal enhanced molecular structures of liquids at oil-silica interfaces.

The recent discovery that molecular CO(2) transforms under compression into carbon four-coordinated, 3-dimensional network solid phases has generated considerable interests on possible new phases in the fourth-main-group elemental oxides. Based on density-functional theory calculations, we have investigated the thermodynamic stability, mechanical properties and electronic structure of proposed guest-free clathrates, quartz and cristobalite phases for CO(2), SiO(2), and GeO(2), and the dry ice phase for CO(2). It was predicted that a GeO(2) clathrate, likely a semiconductor, could be synthesized presumably with some suitable guest molecules. The hypothetical CO(2) guest-free clathrate phase was found hardly to be formed due to the large energy difference with respect to the other polymorphs. This phase is unstable at all pressures, which is also implied by its different electronic structure in comparison with SiO(2) and GeO(2). Finally, the SiO(2) clathrate presents a uniquely high bulk modulus, which is higher than that of quartz and three times of the experimental data, might not be a weak point of ab-initio calculations such as pseudopotentials, correlation functional etc., instead it can be readily understood by the constraint as imposed by the high symmetry. Either temperature or an "exhausted" relaxation (without any symmetry constraint) can remedy this problem.