Proton solvation in methanol solutions of HCl

The MFTIR IR spectra of solutions of HCl in methanol were obtained in the 900–4000 cm^–1 frequency range. It was found that each proton binds two molecules of methanol. The spectra exhibit intense, continuous absorption (CA) with an intensity coefficient at 2000 cm^–1 of 174±10 liter/(mole·cm), which is in agreement with the corresponding coefficient for H₅O ₂ ⁺ . The optical densities of CA are linear functions of the concentration of HCl at 900–1600 cm^–1; there is no linearity at higher frequencies for C_HCl>4 M, and there are less than two molecules of MeOH for each (MeOH)₂H⁺ ion. The results obtained are in agreement with the model in which CA arises in solutions of strong acids because of the interaction of proton vibrations in a strong symmetric H bond with the vibrations of other groups of the proton disolvate.N. N. Semenov Institute of Chemical Physics, Russian Academy of Sciences, 117977 Moscow. Translated from Izvestiya Akademii Nauk, Seriya Khimicheskaya, No. 10, pp. 2261–2268, October, 1992.     

HCl solvation at the surface and within methanol clusters/nanoparticles II: Evidence for molecular wires

Condensed-phase solvation of HCl on and within methanol nanoparticles was investigated by Fourier transform infrared (FTIR) spectroscopy, on-the-fly molecular dynamics as implemented in the density functional code Quickstep (which is part of the CP2K package), and ab initio calculations. Adsorption and solvation stages are identified and assigned with the help of calculated infrared spectra obtained from the simulations. The results have been further checked with MP2-level ab initio calculations. The range of acid solvation states extends from the single-coordinated slightly stretched HCl to proton-sharing with Zundel-like methanol O...H+...X- states, and finally to MeOH2+...Cl- units with full proton transfer. Furthermore, once the proton moves to methanol, it is mobilized along methanol molecular chains. Since the proton dynamics reflects the evolving local structures, the "proton" spectra display broad bands usually with underlying continua.     

Tailoring the mobility of a 3D molecule adsorbed on a metal surface

Because of its tetrahedral structure, spirobifluorene is an innovative molecule for molecular mechanics studies by means of scanning tunneling microscopy. On Cu(100), it was observed only anchored at defects because of its mobility at room temperature. To frustrate its diffusion, it was functionalized with phenyl and thiophene groups. Tetraphenylspirobifluorene is also mobile on Cu(100), whereas tetrathienylspirobifluorene is fixed in the middle of the terraces. This very different behavior is an original and unexpected result because both benzene and thiophene are reported to be weakly bound to Cu(100) with almost the same adsorption energy.     

Liquid-Vapor Equilibrium and Ionization of HCl in the System HCl-H2SO4-H2O at 298 K

The equilibrium pressures and compositions of the vapor over HCl-H₂SO₄-H₂O solutions with a low content of HCl and several fixed concentrations of H₂SO₄ were measured by the flow method at 298 K. The semiconcentration values of the mixed dissociation constants of HCl were calculated and extrapolated to infinite dilution.     

Structure of CO2 adsorbed on the KCl(100) surface

The structure and dynamics of the adsorbate CO(2)/KCl(100) from a diluted phase to a saturated monolayer have been investigated with He atom scattering (HAS), low-energy electron diffraction (LEED), and polarization dependent infrared spectroscopy (PIRS). Two adsorbate phases with different CO(2) coverage have been found. The low-coverage phase is disordered at temperatures near 80 K and becomes at least partially ordered at lower temperatures, characterized by a (2√2×√2)R45° diffraction pattern. The saturated 2D phase has a high long-range order and exhibits (6√2×√2)R45° symmetry. Its isosteric heat of adsorption is 26 ± 4 kJ mol(-1). According to PIRS, the molecules are oriented nearly parallel to the surface, the average tilt angle in the saturated monolayer phase is 10° with respect to the surface plane. For both phases, structure models are proposed by means of potential calculations. For the saturated monolayer phase, a striped herringbone structure with 12 inequivalent molecules is deduced. The simulation of infrared spectra based on the proposed structures and the vibrational exciton approach gives reasonable agreement between experimental and simulated infrared spectra.     

(HCl)2 and (HF)2 in small helium clusters: quantum solvation of hydrogen-bonded dimers

We present a rigorous theoretical study of the solvation of (HCl)(2) and (HF)(2) by small ((4)He)(n) clusters, with n=1-14 and 30. Pairwise-additive potential-energy surfaces of He(n)(HX)(2) (X=Cl and F) clusters are constructed from highly accurate four-dimensional (rigid monomer) HX-HX and two-dimensional (rigid monomer) He-HX potentials and a one-dimensional He-He potential. The minimum-energy geometries of these clusters, for n=1-6 in the case of (HCl)(2) and n=1-5 for (HF)(2), correspond to the He atoms in a ring perpendicular to and bisecting the HX-HX axis. The quantum-mechanical ground-state energies and vibrationally averaged structures of He(n)(HCl)(2) (n=1-14 and 30) and He(n)(HF)(2) (n=1-10) clusters are calculated exactly using the diffusion Monte Carlo (DMC) method. In addition, the interchange-tunneling splittings of He(n)(HCl)(2) clusters with n=1-14 are determined using the fixed-node DMC approach, which was employed by us previously to calculate the tunneling splittings for He(n)(HF)(2) clusters, n=1-10 [A. Sarsa et al., Phys. Rev. Lett. 88, 123401 (2002)]. The vibrationally averaged structures of He(n)(HX)(2) clusters with n=1-6 for (HCl)(2) and n=1-5 for (HF)(2) have the helium density localized in an effectively one-dimensional ring, or doughnut, perpendicular to and at the midpoint of the HX-HX axis. The rigidity of the solvent ring varies with n and reaches its maximum for the cluster size at which the ring is filled, n=6 and n=5 for (HCl)(2) and (HF)(2), respectively. Once the equatorial ring is full, the helium density spreads along the HX-HX axis, eventually solvating the entire HX dimer. The interchange-tunneling splitting of He(n)(HCl)(2) clusters hardly varies at all over the cluster size range considered, n=1-14, and is virtually identical to that of the free HCl dimer. This absence of the solvent effect is in sharp contrast with our earlier results for He(n)(HF)(2) clusters, which show a approximately 30% reduction of the tunneling splitting for n=4. A tentative explanation for this difference is proposed. The implications of our results for the interchange-tunneling dynamics of (HCl)(2) in helium nanodroplets are discussed.     

Reactions of HCl and D2O with molten alkali carbonates

The acidic oxide SO? and protic acid HCl are among the gases released in the combustion of coal and the incineration of municipal waste. They are typically removed by wet or dry scrubbing involving calcium carbonate or calcium hydroxide. The molten alkali carbonate eutectic provides a liquid-state alternative that readily absorbs SO? and HCl and does not become covered with a passivating layer. Gas-liquid scattering experiments utilizing the eutectic mixture (44 mol % Li?CO?, 31 mol % Na?CO?, 25 mol % K?CO?) reveal that the reaction probability for HCl(g) + CO?2? → CO?(g) + OH? + Cl? is 0.31 ± 0.02 at 683 K and rises to 0.39 at 783 K. Gaseous CO? is formed within 10?? s or less, implying that the reaction takes place in a liquid depth of less than 1000 ?. When the melt is exposed to D?O, the analogous reaction D?O(g) + CO?2? → CO?(g) + 2OD? occurs too slowly to measure and no water uptake is observed. Together with previous studies of SO?(g) + CO?2? → CO?(g) + SO?2?, these results demonstrate that molten carbonates efficiently remove both gaseous HCl and SO? while reacting at most weakly with water vapor. The experiments further highlight the remarkable ability of hot CO?2? ions to behave as a base in reactions with protic and Lewis acids.     

Ionization of H2 Rydberg molecules at a metal surface

The ionization of a beam of H2 Rydberg molecules in collision with a metal surface (evaporated Au or Al) is studied. The Rydberg states are excited in an ultraviolet-vacuum ultraviolet double-resonant process and are state selected with a core rotational quantum number N+=0 or 2 and principal quantum numbers n=17-22 (N+=2) or n=41-45 (N+=0). It is found that the N+=0 states behave in a very similar manner to previous studies with atomic xenon Rydberg states, the distance of ionization from the surface scaling with n2. The N+=2 states, however, undergo a process of surface-induced rotational autoionization in which the core rotational energy transfers to the Rydberg electron. In this case the ionization distance scales approximately with nu0(2), the effective principal quantum number with respect to the adiabatic threshold. This process illustrates the close similarity between field ionization in the gas phase and the surface ionization process which is induced by the field due to image charges in the metal surface. The surface ionization rate is enhanced at certain specific values of the field, which is applied in the time interval between excitation and surface interaction. It is proposed here that these fields correspond to level crossings between the N+=0 and N+=2 Stark manifolds. The population of individual states of the N+=2, n=18 Stark manifold in the presence of a field shows that the surface-induced rotational autoionization is more facile for the blueshifted states, whose wave function is oriented away from the surface, than for the redshifted states. The observed processes appear to show little dependence on the chemical nature of the metallic surface, but a significant change occurs when the surface roughness becomes comparable to the Rydberg orbit dimensions.     

Measurements of the heats of dilution and description of the system H2O/H2SO4/HCl with a solvation model

Heat of dilution measurements have been performed on the system H₂O/H₂SO₄/HCl in the concentration range of 0.3 mol% to 20 mol% electrolyte and in the temperature range of 20°C to 60°C. The constituent binary systems H₂O/H₂SO₄ and H₂O/HCl could be described using the solvation model introduced by Engels [H. Engels, Anwendung des modells der lokalen zusammensetzung auf elektroytlösungen, Dissertation, RWTH Aachen, 1985]. The vapor–liquid-equilibria (VLE) and the thermal properties are described simultaneously by one set of parameters. Applying the solvation model to the ternary system H₂O/H₂SO₄/HCl, the measured heats of dilution as well as VLE-data taken from the literature could be predicted. The heat of dilution data are presented with a mean relative deviation of 3.5%. The boiling temperatures of the three sources of ternary VLE data considered in this paper are predicted with a mean absolute deviation of less than 1.2 K. The corresponding mean absolute deviations of the vapor mole fractions of HCl are within 0.06.     

Efficient treatment of solvation shells in 3D molecular theory of solvation

We developed a technique to decrease memory requirements when solving the integral equations of three‐dimensional (3D) molecular theory of solvation, a.k.a. 3D reference interaction site model (3D‐RISM), using the modified direct inversion in the iterative subspace (MDIIS) numerical method of generalized minimal residual type. The latter provides robust convergence, in particular, for charged systems and electrolyte solutions with strong associative effects for which damped iterations do not converge. The MDIIS solver (typically, with 2 × 10 iterative vectors of argument and residual for fast convergence) treats the solute excluded volume (core), while handling the solvation shells in the 3D box with two vectors coupled with MDIIS iteratively and incorporating the electrostatic asymptotics outside the box analytically. For solvated systems from small to large macromolecules and solid–liquid interfaces, this results in 6‐ to 16‐fold memory reduction and corresponding CPU load decrease in MDIIS. We illustrated the new technique on solvated systems of chemical and biomolecular relevance with different dimensionality, both in ambient water and aqueous electrolyte solution, by solving the 3D‐RISM equations with the Kovalenko–Hirata (KH) closure, and the hypernetted chain (HNC) closure where convergent. This core–shell‐asymptotics technique coupling MDIIS for the excluded volume core with iteration of the solvation shells converges as efficiently as MDIIS for the whole 3D box and yields the solvation structure and thermodynamics without loss of accuracy. Although being of benefit for solutes of any size, this memory reduction becomes critical in 3D‐RISM calculations for large solvated systems, such as macromolecules in solution with ions, ligands, and other cofactors. © 2012 Wiley Periodicals, Inc.     

Ionization condition of lithium ionic liquid electrolytes under the solvation effect of liquid and solid solvents

Ionization condition and ionic structures of the lithium ionic liquid electrolytes, LiTFSI/EMI-TFSI/(PEG or silica), were investigated through the measurements of ionic conductivity and diffusion coefficient. The size of the hydrodynamic lithium species (rLi) evaluated from the Stokes-Einstein equation was 0.90 nm before gelation with the PEG or silica. This reveals that the TFSI- anions from the solvent are coordinated on Li+ for solvation, forming, for example, Li(TFSI)4(3-) and Li(TFSI)2- in the electrolyte solution. By the dispersion of PEG for gelation, rLi increased up to 1.8 nm with the 10 wt % of PEG. This indicates that the lithium species is directly interacted with the oxygen sites on the polymer chains and the lithium species migrate, reflecting the polymer by hopping from site to site. In case of the silica dispersion, rLi decreased to 0.7 nm at 10 wt % silica. Although the silica surface with silanol groups fundamentally attracts Li+, the lithium does not migrate from site to site on the silica surface as in the gel of the polymer and follows random walk behavior in the network of the liquid-phase pathways in the two-phase gel. In the process, that solvated TFSI- anions are partially removed may be due to the attractive effect of H+, which was dissociated from the silanol group. It is concluded that the dispersed silica is effective to modify the hydrodynamic lithium species to be appropriate for charge transport as reducing the size and anionic charge of Li(TFSI)4(3-) by removing one or two TFSI- anions.     

吸附氢气后的钯团簇在氧化亚铜表面上的稳定性研究

为了探究吸附H2后的Pdn团簇在Cu2O（111）完整表面和铜缺陷表面上的稳定性,计算了负载在Cu2O（111）完整表面和铜缺陷表面上的Pdn（n=1-4）对H2分子的最稳定吸附结构;利用在给定H2压力和温度下Pdn / Cu2O表面吸附H2的相图揭示了Pdn团簇在Cu2O（111）两个表面的变化情况。结果表明,在吸附了H2分子以后,Pdn团簇更倾向于保持原有的结构,且随着Pd团簇的增大,吸附H2的数量也逐渐增长。     

Prediction of the orientations of adsorbed protein using an empirical energy function with implicit solvation

When simulating protein adsorption behavior, decisions must first be made regarding how the protein should be oriented on the surface. To address this problem, we have developed a molecular simulation program that combines an empirical adsorption free energy function with an efficient configurational search method to calculate orientation-dependent adsorption free energies between proteins and functionalized surfaces. The configuration space is searched systematically using a quaternion rotation technique, and the adsorption free energy is evaluated using an empirical energy function with an efficient grid-based calculational method. In this paper, the developed method is applied to analyze the preferred orientations of a model protein, lysozyme, on various functionalized alkanethiol self-assembled monolayer (SAM) surfaces by the generation of contour graphs that relate adsorption free energy to adsorbed orientation, and the results are compared with experimental observations. As anticipated, the adsorbed orientation of lysozyme is predicted to be dependent on the discrete organization of the functional groups presented by the surface. Lysozyme, which is a positively charged protein, is predicted to adsorb on its 'side' on both hydrophobic and negatively charged surfaces. On surfaces with discrete positively charged sites, attractive interaction energies can also be obtained due to the presence of discrete local negative charges present on the lysozyme surface. In this case, 'end-on' orientations are preferred. Additionally, SAM surface models with mixed functionality suggest that the interactions between lysozyme and surfaces could be greatly enhanced if individual surface functional groups are able to access the catalytic cleft region of lysozyme, similar to ligand-receptor interactions. The contour graphs generated by this method can be used to identify low-energy orientations that can then be used as starting points for further simulations to investigate conformational changes induced in protein structure following initial adsorption.     

Complex investigation of adsorbed solvation layers in polymer-containing aqueous suspensions of zinc oxide

Theoretical and Experimental Chemistry -     

O（~1D）与CF_2HCl反应的理论研究

用量子化学密度泛函理论（DFT）和G3B3方法，对为（~1D）与CF_2HCl的反应 进行了研究，在B3LYP/6-311+G(d)，B3LYP/6-311+G(2df,2pd)和G3B3计算水平上， 优化了反应热能面上各驻点的几何结构，通过内禀反应坐标（IRC）计算和振动分 析，对反应过渡态进行了确认，并确定了反应机理。     

Dynamics of D2 released from the dissociation of D2O on a zirconium surface

Hydrogen is efficiently released during water dissociation on zirconium (Zr), while even very rapid temperature programmed heating of a hydrogen covered Zr surface predominantly leads to dissolution (approximately 99% dissolution). To help resolve these apparently contradictory observations, we have studied the dynamics of water (D2O) dissociation on a crystalline Zr surface by probing the rotational and vibrational energy distributions of the D2 produced using resonant enhanced multiphoton ionization spectroscopy. The internal-state energy distribution of the D2 product was found to be rotationally cold and vibrationally hot with respect to the temperature of the surface. The rotational distribution shows slight deviations from Boltzmann's law, with a mean rotational temperature of 426 K while the surface is at 800 K. The population of the nu"=1 vibration is at least four times higher than a 800 K temperature would allow, this corresponding to a vibrational temperature of 1100 K. Information on the translational energy of the D2 product have also been obtained by time-of-flight spectroscopy and it is found to be nearly thermally equilibrated with the surface temperature. Similar results were obtained from studies of D2 scattered from a clean Zr surface, and of D2 released by a slow thermal desorption process which involves dissolved hydrogen as the source. The reconciliation of the present results with those for thermal desorption of preadsorbed hydrogen implies a role for both surface and subsurface adsorption sites on the Zr surface and clearly demonstrates that at high temperatures, the release of D2 arises from the recombinative desorption of adsorbed hydrogen formed by the complete dissociation of D2O.     

The states of hydrogen adsorbed on the surface of skeleton nickel

The experimental data on hydrogen adsorption on skeleton nickel were used to construct the pK spectrum of atomic hydrogen on this surface. The pK spectrum obtained corresponded to the presence of three adsorption centers on the surface of skeleton nickel. A chemical model of hydrogen adsorption on the surface of skeleton nickel was suggested.