UV-induced protonation of molecules adsorbed on ice surfaces at low temperature

UV irradiation of ice films adsorbed with methylamine molecules induces protonation of the adsorbate molecules at low temperature (50-130 K). The observation indicates that long-lived protonic defects are created in the ice film by UV light, and they transfer protons to the adsorbate molecules via tunneling mechanism at low temperature. The methylammonium ion formed by proton transfer remains to be stable at the ice surface. It is suggested that this solid-phase protonation might play a significant role in the production of molecular ions in interstellar clouds.     

Exceptional Isotopic‐Substitution Effect: Breakdown of Collective Proton Tunneling in Hexagonal Ice due to Partial Deuteration

Multiple proton transfer controls many chemical reactions in hydrogen‐bonded networks. However, in contrast to well‐understood single proton transfer, the mechanisms of correlated proton transfer and of correlated proton tunneling in particular have remained largely elusive. Herein, fully quantized ab initio simulations are used to investigate H/D isotopic‐substitution effects on the mechanism of the collective tunneling of six protons within proton‐ordered cyclic water hexamers that are contained in proton‐disordered ice, a prototypical hydrogen‐bonded network. At the transition state, isotopic substitution leads to a Zundel‐like complex, [HO???D???OH], which localizes ionic defects and thus inhibits perfectly correlated proton tunneling. These insights into fundamental aspects of collective proton tunneling not only rationalize recent neutron‐scattering experiments, but also stimulate investigations into multiple proton transfer in hydrogen‐bonded networks much beyond ice.     

Electrical processes upon the evaporation and condensation of water and ice

A mathematical model of electrical processes that take place upon the evaporation of water and sublimation of ice, as well as the condensation growth of water and ice phases from vapor, is proposed. The transfer of the main charge carriers, such as (i) protons and hydroxide ions (in ice, water, and vapor and (ii) orientational defects (in ice and water) is taken into account. Upon the evaporation of water and the sublimation of ice, the first carriers are accumulated before the phase front and cause positive charges in the surface of the water and ice, while the second carriers are depleted (their concentration becomes lower than the thermodynamic value) and impart a negative charge to water and ice. The contribution of protons and hydroxide ions dominates at a low rate of evaporation. In the condensation of vapor and relevant growth of water and ice phases, the polarity of surface charge is opposite to that observed upon evaporation. The values of interfacial current and signs of phase charges upon sublimation, evaporation, and condensation that are predicted in the model comply with experimental data.     

Evidence for the surface origin of point defects in ice: control of interior proton activity by adsorbates

Spectroscopic studies are presented of H-D isotopic exchange in the interior of ice nanocrystals. The exchange process is dominated by ionic and orientational defects long viewed as governing the electrical properties of ice. A new finding that interior exchange rates can be controlled by acidic and basic adsorbates is evidence that the defects originate at the ice surface. In particular, it is argued that interior isotopic exchange is a reflection of proton concentrations equilibrated at the ice surface.     

Mechanistic study of proton transfer and HD exchange in ice films at low temperatures (100-140 K)

We have examined the elementary molecular processes responsible for proton transfer and HD exchange in thin ice films for the temperature range of 100-140 K. The ice films are made to have a structure of a bottom D(2)O layer and an upper H(2)O layer, with excess protons generated from HCl ionization trapped at the D(2)OH(2)O interface. The transport behavior of excess protons from the interfacial layer to the ice film surface and the progress of the HD exchange reaction in water molecules are examined with the techniques of low energy sputtering and Cs(+) reactive ion scattering. Three major processes are identified: the proton hopping relay, the hop-and-turn process, and molecular diffusion. The proton hopping relay can occur even at low temperatures (<120 K), and it transports a specific portion of embedded protons to the surface. The hop-and-turn mechanism, which involves the coupling of proton hopping and molecule reorientation, increases the proton transfer rate and causes the HD exchange of water molecules. The hop-and-turn mechanism is activated at temperatures above 125 K in the surface region. Diffusional mixing of H(2)O and D(2)O molecules additionally contributes to the HD exchange reaction at temperatures above 130 K. The hop-and-turn and molecular diffusion processes are activated at higher temperatures in the deeper region of ice films. The relative speeds of these processes are in the following order: hopping relay>hop and turn>molecule diffusion.     

Proton mobility in thin ice films: a revisit

We have examined proton transport through an ice film in the temperature range 73-140 K by initially adding hydronium ions into the interior of the film and then monitoring the build-up of hydronium ion population at the film surface. The result confirms that the proton exhibits limited mobility in the ice film at low temperature, but it becomes highly mobile at temperature above 130 K. Based on this result we suggest an explanation of the anomalous experimental observations in the literature for the proton mobility in ice films.     

Mechanistic aspects of proton chain transfer: a computational study for the green fluorescent protein chromophore

We explore several models for the ground-state proton chain transfer pathway between the green fluorescent protein chromophore and its surrounding protein matrix, with a view to elucidating mechanistic aspects of this process. We have computed quantum chemically the minimum energy pathways (MEPs) in the ground electronic state for one-, two-, and three-proton models of the chain transfer. There are no stable intermediates for our models, indicating that the proton chain transfer is likely to be a single, concerted kinetic step. However, despite the concerted nature of the overall energy profile, a more detailed analysis of the MEPs reveals clear evidence of sequential movement of protons in the chain. The ground-state proton chain transfer does not appear to be driven by the movement of the phenolic proton off the chromophore onto the neutral water bridge. Rather, this proton is the last of the three protons in the chain to move. We find that the first proton movement is from the bridging Ser205 moiety to the accepting Glu222 group. This is followed by the second proton moving from the bridging water to the Ser205--for our model this is where the barrier occurs. The phenolic proton on the chromophore is hence the last in the chain to move, transferring to a bridging "water" that already has substantial negative charge.     

H/D isotopic exchange between water molecules at ice surfaces

H/D isotopic exchange between H(2)O and D(2)O molecules was studied at the surface of ice films at 90-140 K by the technique of Cs(+) reactive ion scattering. Ice films were deposited on a Ru(0001) substrate in different compositions of H(2)O and D(2)O and in various structures to study the kinetics of isotopic exchange. H/D exchange was very slow on an ice film at 95-100 K, even when H(2)O and D(2)O were uniformly mixed in the film. At 140 K, H/D exchange occurred in a time scale of several minutes on the uniform mixture film. Kinetic measurement gave the rate coefficient for the exchange reaction, k(140 K)=1.6(+/-0.3) x 10(-19) cm(2) molecule(-1) s(-1) and k(100 K)< or =5.7(+/-0.5) x 10(-21) cm(2) molecule(-1) s(-1) and the Arrhenius activation energy, E(a)> or =9.8 kJ mol(-1). Addition of HCl on the film to provide excess protons greatly accelerated the isotopic exchange reaction such that it went to completion very quickly at the surface. The rapid reaction, however, was confined within the first bilayer (BL) of the surface and did not readily propagate to the underlying sublayer. The isotopic exchange in the vertical direction was almost completely blocked at 95 K, and it slowly occurred only to a depth of 3 BLs from the surface at 140 K. Thus, the proton transfer was highly directional. The lateral proton transfer at the surface was attributed to the increased mobility of protonic defects at the molecularly disordered and activated surface. The slow, vertical proton transfer was probably assisted by self-diffusion of water molecules.     

Sticking of CO to crystalline and amorphous ice surfaces

We present results of classical trajectory calculations on the sticking of hyperthermal CO to the basal plane (0001) face of crystalline ice Ih and to the surface of amorphous ice Ia. The calculations were performed for normal incidence at a surface temperature Ts = 90 K for ice Ia, and at Ts = 90 and 150 K for ice Ih. For both surfaces, the sticking probability can be fitted to a simple exponentially decaying function of the incidence energy, Ei: Ps = 1.0e(-Ei(kJ/mol)/90(kJ/mol)) at Ts = 90 K. The energy transfer from the impinging molecule to the crystalline and the amorphous surface is found to be quite efficient, in agreement with the results of molecular beam experiments on the scattering of the similar molecule, N2, from crystalline and amorphous ice. However, the energy transfer is less efficient for amorphous than for crystalline ice. Our calculations predict that the sticking probability decreases with Ts for CO scattering from crystalline ice, as the energy transfer from the impinging molecule to the warmer surfaces becomes less efficient. At high Ei (up to 193 kJ/mol), no surface penetration occurs in the case of crystalline ice. However, for CO colliding with the amorphous surface, a penetrating trajectory was observed to occur into a large water pore. The molecular dynamics calculations predict that the average potential energy of CO adsorbed to ice Ih is -10.1 +/- 0.2 and -8.4 +/- 0.2 kJ/mol for CO adsorbed to ice Ia. These values are in agreement with previous experimental and theoretical data. The distribution of the potential energy of CO adsorbed to ice Ia was found to be wider (with a standard deviation sigma of 2.4 kJ/mol) than that of CO interacting with ice Ih (sigma = 2.0 kJ/mol). In collisions with ice Ia, the CO molecules scatter at larger angles and over a wider distribution of angles than in collisions with ice Ih.     

Non-joule heating of ice in an electric field

We theoretically predict and calculate non-Joule heating/cooling caused by a direct electric current in ordinary crystalline ice Ih. The cause of this effect is related to partial ordering/disordering occurring in the proton subsystem of ice when protons either drift or diffuse in the ice. Depending on relative directions of the electric current and the configuration vector of ice, the non-Joule effect can be either positive, that is, heat generation, or negative, that is, heat absorption, and its absolute magnitude is usually comparable with that of normal Joule heating. The magnitude of this phenomenon is also approximately inversely proportional to the ice temperature and, thus, is more pronounced at low temperatures.     

Combinatorial entropy and phase diagram of partially ordered ice phases

A close analytical estimate for the combinatorial entropy of partially ordered ice phases is presented. The expression obtained is very general, as it can be used for any ice phase obeying the Bernal-Fowler rules. The only input required is a number of crystallographic parameters, and the experimentally observed proton site occupancies. For fully disordered phases such as hexagonal ice, it recovers the result deduced by Pauling, while for fully ordered ice it is found to vanish. Although the space groups determined for ice I, VI, and VII require random proton site occupancies, it is found that such random allocation of protons does not necessarily imply random orientational disorder. The theoretical estimate for the combinatorial entropy is employed together with free energy calculations in order to obtain the phase diagram of ice from 0 to 10 GPa. Overall qualitative agreement with experiment is found for the TIP4P model of water. An accurate estimate of the combinatorial entropy is found to play an important role in determining the stability of partially ordered ice phases, such as ice III and ice V.     

Guanidinium groups act as general-acid catalysts in phosphoryl transfer reactions: a two-proton inventory on a model system

Cleavage/transesterification of phosphodiesters is catalyzed by various acidic groups in solution and with enzymes. General-acid catalysts can transfer protons to the developing phosphorane intermediate, resulting in a monoprotic-monoanionic intermediate, giving the so-called "triester mechanism". Using a proton inventory on a model compound (1) possessing an intramolecular hydrogen bond between a phosphodiester and a guanidinium group, we find that two protons move in the rate-determining step for cleavage/transesterification. In contrast, HPNP shows a single-proton inventory and is a substrate well accepted to react with the movement of only one proton at the transition state. We therefore propose a mechanism for 1 that involves general-acid catalysis by the guanidinium group. This leads one to conclude that other, more acidic groups, such as ammonium and imidazolium, would also act as general-acid catalysts.     

Voltammetry of azobenzene microcrystals

Azobenzene solid particles have been mechanically attached to a graphite electrode and measured by cyclic staircase voltammetry. Well-developed and widely separated cathodic and anodic peaks were observed. Redox reaction is controlled by both the heterogeneous charge transfer kinetics and the mass transfer. Its mechanism is explained by the surface diffusion model. Reaction starts at the three-phase boundary, where electrons are transferred from the electrode surface to the attached azobenzene molecules. The electrons are then propagated over the surface of microcrystals by a series of exchange reactions between hydrazobenzene and azobenzene molecules, with the participation of proton donors from the solution. The apparent mass transfer occurs in this way. In the crystal lattice the transmission of protons is not possible, and consequently there is no propagation of electrons. Received: 31 January 1997 / Accepted: 21 February 1997     

Interaction of low-energy ions and atoms of light elements with a fluorinated carbon molecular lattice

The mechanism of interaction of low-energy atoms and ions of light elements (H, H+, He, Li, the kinetic energy of the particles 2-40 eV) with C6H6, C6F12, C60, and C60F48 molecules was studied by ab initio MD simulations and quantum-chemical calculations. It was shown that starting from 6 A from the carbon skeleton for the "C6H6 + proton" and "C60 + proton" systems, the electronic charge transfer from the aromatic molecule to H+ occurs with a probability close to 1. The process transforms the H+ to a hydrogen atom and the neutral C6H6 and C60 molecules to cation radicals. The mechanism of interaction of low-energy protons with C6F12 and C60F48 molecules has a substantially different character and can be considered qualitatively as the interaction between a neutral molecule and a point charge. The Coulomb perturbation of the system arising from the interaction of the uncompensated proton charge with the Mulliken charges of fluorine atoms results in an inversion of the energies of the electronic states localized on the proton and on the C6F12 and C60F48 molecules and makes the electronic charge transfer energetically unfavorable. On the different levels of theory, the barriers of the proton penetration for the C6F12 and C60F48 molecules are from two to four times lower than those for the corresponding parent systems (C6H6 and C60). The penetration barriers of the He atom and Li+ ion depend mainly on the effective radii of the bombarding particles. The theoretical penetration and escaped barriers for the "Li+ + C60" process qualitatively explain the experimental conditions of synthesis of the Li@C60 complex.     

470—The proton pump at work: Part I. The basic concept of excess proton/defect proton translocation

The Nagle-Morowitz proton pump, which is based on proton transport in water and ice, is shown to be inapplicable to weakly H-bonded proton transport systems. It is demonstrated that in weakly or non H-bonded systems protons can migrate either as excess protons H using a high-lying proton conduction band or as defect protons H′ using intermediate energy levels. As H? the protons act as positive charge carriers endowed with a very high mobility. As H′ they act as negative charge carriers. The H′ transport mechanism involves thermally activated proton tunneling. Owing to the very large mobility differences between H? and H′, > 10⁶, very high potentials can be self-generated in any situation which creates a concentration gradient. Proton conductivity data on inorganic model compounds are presented. Applying these results to proton transport across biomembranes, transmembrane potentials U_m, acidification Δ pH and transient phenomena can be explained as result of H? and H′ translocation.     

Theory of the dielectric constant of ice

The first result of this paper is to show that the Onsager—Slater theory of the dielectric constant is consistent for some reasoable model of ice in the limit of no intrinsic defects. Accordingly, a model is presented, called the unit model, which has unit cells with no dipole moments for which the Onsager—Slater theory is exact. The second result of this paper is to show that the unit model is physically and chemically realistic. Bjerrum defects are introduced into the model and the relation between the dielectric constant and the Bjerrum defect charge found by Onsager and Dupuis for a less realistic model continues to hold and is satisfied by the experimental values. In a simple point charge approximation the charge distribution determined by the model requirements and the experimentally determined Bjerrum fault charge are found and seem reasonble. Higher order multipole interaction energies are consistent with eviations from pure 1/T dependence of the dielectric constant of real ice with intrinsic defects, can be derived in the context of the unit model. This formula interpolates between the Onsager—Slater formula in the limit of no intrinsic defects and the Kirkwood—Frohlich formula in the limit of many intrinsic defects. For the concentration of defects in real ice the interpolation formula is practically the same as the Onsager—Slater formula and differs from the Kirkwood—Frohlich formula by a factor of nearly $\text{3}\text{2}$.     

Adsorption of humic acids onto goethite: effects of molar mass, pH and ionic strength

In this paper, the LCD (ligand charge distribution) model is applied to describe the adsorption of (Tongbersven) humic acid (HA) to goethite. The model considers both electrostatic interactions and chemical binding between HA and goethite. The large size of HA particles limits their close access to the surface. Part of the adsorbed HA particles is located in the compact part at the goethite surface (Stern layers) and the rest in the less structured diffuse double layer (DDL). The model can describe the effects of pH, ionic strength, and loading on the adsorption. Compared to fulvic acid (FA), adsorption of HA is stronger and more pH- and ionic-strength-dependent. The larger number of reactive groups on each HA particle than on a FA particle results in the stronger HA adsorption observed. The stronger pH dependency in HA adsorption is related to the larger number of protons that are coadsorbed with HA due to the higher charge carried by a HA particle than by a FA particle. The positive ionic-strength dependency of HA adsorption can be explained by the conformational change of HA particles with ionic strength. At a higher ionic strength, the decrease of the particle size favors closer contact between the particles and the surface, leading to stronger competition with electrolyte ions for surface charge neutralization and therefore leading to more HA adsorption.     

Anomalous temperature-isotope dependence in proton-coupled electron transfer

Motivated by the experiments of Hodgkiss et al. [J. Phys. Chem. (submitted)] on electron transfer (ET) through a H-bonding interface, we present a new theoretical model for proton-coupled electron transfer (PCET) in the condensed phase, that does not involve real proton transfer. These experiments, which directly probe the joint T-isotope effects in coupled charge transfer reactions, show anomalous T dependence in k(H)k(D), where k(H) and k(D) are the ET rates through the H-bonding interface with H-bonded protons and deuterons, respectively. We address the anomalous T dependence of the k(H)k(D) in our model by attributing the modulation of the electron tunneling dynamics to bath-induced fluctuations in the proton coordinate, so that the mechanism for coupled charge transfer might be better termed vibrationally assisted ET rather than PCET. We argue that such a mechanism may be relevant to understanding traditional PCET processes, i.e., those in which protons undergo a transfer from donor to acceptor during the course of ET, provided there is an appropriate time scale separating both coupled charge transfers. Likewise, it may also be useful in understanding long-range ET in proteins, where tunneling pathways between redox cofactors often pass through H-bonded amino acid residues, or other systems with sufficiently decoupled proton and electron donating functionalities.     

Dynamics of colloidal particles in ice

We use x-ray photon correlation spectroscopy (XPCS) to probe the dynamics of colloidal particles in polycrystalline ice. During freezing, the dendritic ice morphology and rejection of particles from the ice created regions of high particle density, where some of the colloids were forced into contact and formed disordered aggregates. The particles in these high density regions underwent ballistic motion, with a characteristic velocity that increased with temperature. This ballistic motion is coupled with both stretched and compressed exponential decays of the intensity autocorrelation function. We suggest that this behavior could result from ice grain boundary migration.