Bulk Water Molecules Research Articles

Kinetics and mechanistic interpretation of fluoride removal by nanocrystalline hydroxyapatite derived from Limacine artica shells

The present investigation describes the synthesis and characterization of nanocrystalline adsorbent, hydroxyapatite (HAp) from Limacine artica shells and its application towards fluoride removal from aqueous solution. The high content of calcium carbonate (92.9%) analyzed by inductively coupled plasma optical emission spectrometry (ICP-OES) confirmed the snail shell to be a good precursor for HAp synthesis. The X-ray diffraction (XRD) and Braunauer-Emmet-Teller (BET) analyses supported the formation of a nanocrystalline material with a surface area and pore diameter of 152.45m2/g and 4.2nm, respectively. The observed decrease in crystallite size after fluoride ion adsorption was attributed to the dispersion of agglomerated HAp particles due to H-bonding between adsorbed fluoride ions with bulk water molecules. The good fitting of Lagergren pseudo second-order kinetics and Elovich model indicated the rate limiting step to be of chemisorptions type. The adsorption data followed Langmuir, Freundlich and Dubinin-Radushkevich (D-R) isotherms satisfactorily supporting both monolayer and multilayer physisorption at HAp-water interface. Thermodynamic study revealed the adsorption process to be spontaneous and endothermic in nature. Electrostatic attraction, ion exchange and hydrogen bonding are mainly responsible for the adsorption mechanism. The used HAp could be reused for fourth cycle without a significant loss of adsorption ability.

Elucidation of the Catalytic Mechanism of a Miniature Zinc Finger Hydrolase.

To improve our mechanistic understanding of zinc metalloenzymes, we report a joint computational and experimental study of a minimal carbonic anhydrase (CA) mimic, a 22-residue Zn-finger hydrolase. We combine classical molecular dynamics (MD) simulations, quantum mechanics/molecular mechanics (QM/MM) geometry optimizations, and QM/MM free energy simulations with ambient and high-pressure kinetic measurements to investigate the mechanism of the hydrolysis of the substrate p-nitrophenylacetate (pNPA). The zinc center of the hydrolase prefers a pentacoordinated geometry, as found in most naturally occurring CAs and CA-like enzymes. Two possible mechanisms for the catalytic reaction are investigated. The first one is analogous to the commonly accepted mechanism for CA-like enzymes: a sequential pathway, in which a Zn2+-bound hydroxide acts as a nucleophile and the hydrolysis proceeds through a tetrahedral intermediate. The initial rate-limiting step of this reaction is the nucleophilic attack of the hydroxide on pNPA to form the tetrahedral intermediate. The computed free energy barrier of 18.5 kcal/mol is consistent with the experimental value of 20.5 kcal/mol obtained from our kinetics experiments. We also explore an alternative reverse protonation pathway for the hydrolase, in which a nearby hydroxide ion from the bulk acts as the nucleophile (instead of a zinc-bound hydroxide). According to QM/MM MD simulations, hydrolysis occurs spontaneously along this pathway. However, this second scenario is not viable in our system, as the tertiary structure of the hydrolase lacks a suitably positioned residue that would act as a general base and generate a hydroxide ion from a nearby bulk water molecule. Hence, our combined theoretical and experimental study indicates that the investigated minimal CA mimic retains the essential mechanistic features of CA-like enzyme catalysis. The high-pressure experiments show that its catalytic efficiency can be enhanced by applying hydrostatic pressure. According to the simulations, more drastic improvements might be afforded by mutations that make the reverse protonation pathway accessible.

Origins of Water Molecules in the Photosystem II Crystal Structure.

The cyanobacterial photosystem II (PSII) crystal structure includes more than 1300 water molecules in each monomer unit; however, their precise roles in water oxidation are unclear. To understand the origins of water molecules in the PSII crystal structure, the accessibility of bulk water molecules to channel inner spaces in PSII was investigated using the water-removed PSII structure and molecular dynamics (MD) simulations. The inner space of the channel that proceeds toward the D1-Glu65/D2-Glu312 pair (E65/E312 channel) was entirely filled with water molecules from the bulk region. In the same channel, a diamond-shaped cluster of water molecules formed near redox-active TyrZ in MD simulations. Reorientation of the D2-Leu352 side chain resulted in formation of a hexagonal water network at the Cl-2 binding site. Water molecules could not enter the main region of the O4-water chain, which proceeds from the O4 site of the Mn4CaO5 cluster. However, in the O4-water chain, the two water binding sites that are most distant from the protein bulk surface were occupied by water molecules that approached along the E65/E312 channel, one of which formed an H-bond with the O4 site. These findings provide key insights into the significance of the channel ends, which may utilize water molecules during the PSII photocycle.

Raman scattering-based multiconformational analysis for probing the structural differences between acetylcholine and acetylthiocholine

Acetylcholine is the first discovered neurotransmitter that has received a great attention regarding its capability of binding to several cellular targets. The chemical composition of acetylcholine, including a positively charged trimethylammonium and a carbonyl group, as well as its conformational flexibility was pointed out as the key factors in the stabilization of its interactions. Here, the possibilities offered by a Raman scattering-based multiconformatioal analysis to access the most stable conformers of acetylcholine, is discussed. To control the validity of this protocol, acetylcholine and one of its closely structured analogues, acetylthiocholine, were simultaneously analyzed. Solution Raman spectra revealed distinct and well resolved strong markers for each molecule. Density functional theory calculations were consistent with the fact that the energy order of the low energy conformers is considerably affected by the acyloxy oxygen→sulfur atom substitution. Raman spectra were calculated on the basis of the thermal average of the spectra arising from the low energy conformers. It has been evidenced that the carbonyl and trimethylammonium groups are the most favorable hydration sites in aqueous environment. Taking into account the large gap between the carbonyl bond-stretch and aliphatic bending bands, Raman spectra also allowed separation of the HOH bending vibrations arising from the bound and bulk water molecules.

A hydrated ion model of [UO2]2+ in water: Structure, dynamics, and spectroscopy from classical molecular dynamics.

A new ab initio interaction potential based on the hydrated ion concept has been developed to obtain the structure, energetics, and dynamics of the hydration of uranyl in aqueous solution. It is the first force field that explicitly parameterizes the interaction of the uranyl hydrate with bulk water molecules to accurately define the second-shell behavior. The [UO2(H2O)5]2+ presents a first hydration shell U-O average distance of 2.46 Å and a second hydration shell peak at 4.61 Å corresponding to 22 molecules using a coordination number definition based on a multisite solute cavity. The second shell solvent molecules have longer mean residence times than those corresponding to the divalent monatomic cations. The axial regions are relatively de-populated, lacking direct hydrogen bonding to apical oxygens. Angle-solved radial distribution functions as well as the spatial distribution functions show a strong anisotropy in the ion hydration. The [UO2(H2O)5]2+ solvent structure may be regarded as a combination of a conventional second hydration shell in the equatorial and bridge regions, and a clathrate-like low density region in the axial region. Translational diffusion coefficient, hydration enthalpy, power spectra of the main vibrational modes, and the EXAFS spectrum simulated from molecular dynamics trajectories agree fairly well with the experiment.

On the solvation properties of supercritical electrolyte solutions

A series of hydrothermal diamond anvil cell experiments was conducted to constrain the effect of anions (Cl−, F−) and cations (Mg+2, Na+) on the molecular structure of high temperature/-pressure fluids, and to express the extent of hydration as a function of solute's speciation and concentration at temperatures ranging from 400 to 800°C and pressures of 1.7 to 30.1kbar. Results reveal that hydrogen-bonding environments for H2O molecules in the first hydration shell of anions evolve from pure OH⋯O to dynamic mixtures with OH⋯X− (F−, Cl−). The distribution of the “OH⋯X−” structures (OH⋯O, OH⋯X−) is proportional to the concentration of dissolved anions; strongly dependent, however, on ionic speciation. To this end, the OH⋯O structures inside the hydration shell of F− are weaker than the OH⋯Cl− and OH⋯O in Cl−bearing solutions. On the contrary, cations (Na+, Mg2+) appear to impose a minimal effect on the Raman spectra of the OH⋯O structures. The thermodynamic stability of “OH⋯Cl−” environments shows a strong correlation with the activity of H2O, indicative of the localized effect of the Cl− electric field that decouples the hydration from the bulk water molecules. The enthalpy change required to rupture these H-bonding environments (ΔHOH⋯Cl−) provides a measure for the extent of ion hydration and constrains the solvation properties of supercritical electrolyte solutions. For example, the dominant presence of “OH⋯Cl− “ environments lowers the concentration of isolated water dipoles and contributes to the decrease in the dielectric constant of supercritical electrolyte solutions, which in turn, could suppress the solubility of hydrophobic non-polar aqueous species, such as SiO2(aq).

Effect of polyvinylpyrrolidone at methane hydrate-liquid water interface. Application in flow assurance and natural gas hydrate exploitation

Polyvinylpyrrolidone (PVP) is one of the most studied low dosages hydrate inhibitor (LDHI), while its hydrate inhibiting effect is well known, its surface active properties by which PVP molecules alters the hydrate-liquid water interface has not been understood properly. In the present work, influence of PVP molecules at the hydrate-liquid water interface was studied at molecular level using molecular dynamics (MD) simulation. In addition, impact of various molecular weights (or chain length) of PVP molecules was also investigated. The force field parameters for PVP monomer (F4) and polymers were developed and validated against PVP physical properties. The free energy of binding of PVP with methane hydrate and methane hydrate decomposition kinetics was studied in presence of PVP at hydrate-water interface. Structural properties of hydrate were analyzed using four body order parameter and mutually coordinated guest order parameter (MCG-OP). The decomposition rate of methane hydrate in presence of PVP molecules (in bulk water near hydrate interface) was studied, it was observed that PVP polymer changes the surface properties of hydrate and enhances the hydrate decomposition rate.

Exploration of Binding Interactions of Cu(2+) with d-Penicillamine and its O- and Se- Analogues in Both Gas and Aqueous Phases: A Theoretical Approach.

We have theoretically explored the entire binding phenomena of d-penicillamine and its O- and Se-analogues with Cu(2+) in both gas and aqueous phases. At first, a brief conformational analysis has been performed via -XH and -COOH rotations to investigate such conformers that are suitable for binding in both bidentate as well as tridentate fashions. The stability of each bidentate and tridentate complex is determined on the basis of relative energy (ΔE) and gas phase metal ion affinity (MIA) along with the bonding analysis by using atoms in molecule theory. The effect of conformational change on the stability of the complexes is also examined thoroughly. By analyzing the MIA values, we have shown that the side chain substitution makes an impact on the binding process. To delve into the binding phenomena in aqueous phase, we have introduced both the first and second hydration sphere models. In first hydration sphere model, to realize the precise effect of water molecules we have considered stable octahedral hexa-aqua copper complex, [Cu(H2O)6](+2) and accordingly substituted water molecules depending on the bidentate or tridentate nature of the chelating agents. The influence of bulk water molecules on the energetics and geometries of the first hydrated sphere complexes have also been investigated by employing second hydration sphere model assuming physiological pH through the implementation of implicit COSMO and polarizable continuum models, respectively. In the second hydration sphere model, the zwitterionic structures of the amino acids and their side chain deprotonated forms are also included to study the binding phenomena with Cu(2+). The complete work furnishes both the binding properties and the energetics of the copper-artificial amino acid complexes in both gas and aqueous phases that will reflect a realistic overview of the entire binding phenomena.

Label-Free Fluctuation Spectroscopy Based on Coherent Anti-Stokes Raman Scattering from Bulk Water Molecules.

Nanoparticles (NPs) and molecules can be analyzed by inverse fluorescence correlation spectroscopy (iFCS) as they pass through an open detection volume, displacing fractions of the fluorescence‐emitting solution in which they are dissolved. iFCS does not require the NPs or molecules to be labeled. However, fluorophores in μm–mm concentrations are needed for the solution signal. Here, we instead use coherent anti‐Stokes Raman scattering (CARS) from plain water molecules as the signal from the solution. By this fully label‐free approach, termed inverse CARS‐based correlation spectroscopy (iCARS‐CS), NPs that are a few tenths of nm in diameter and at pM concentrations can be analyzed, and their absolute volumes/concentrations can be determined. Likewise, lipid vesicles can be analyzed as they diffuse/flow through the detection volume by using CARS fluctuations from the surrounding water molecules. iCARS–CS could likely offer a broadly applicable, label‐free characterization technique of, for example, NPs, small lipid exosomes, or microparticles in biomolecular diagnostics and screening, and can also utilize CARS signals from biologically relevant media other than water.

Identification of Multiple Water-Iodide Species in Concentrated NaI Solutions Based on the Raman Bending Vibration of Water.

The influence of aqueous electrolytes on the water bending vibration was studied with Raman spectroscopy. For all salts investigated (NaI, NaBr, NaCl, and NaSCN), we observed a nonlinear intensity increase of the water bending vibration with increasing concentration. Different lasers and a tunable frequency-doubled optical parametric oscillator system were used to achieve excitation wavelengths between 785 and 374 nm. Focusing on NaI solutions, the relative enhancement of the water bending vibration was found to increase strongly with excitation photon energy, in line with a preresonance effect from the iodide-water charge-transfer transition. We used multivariate curve resolution (MCR) to decompose the measured Raman spectra of NaI solutions into three interconverting spectral components assigned to bulk water and water molecules interacting with one (X···H-O-H···O) and two (X···H-O-H···X) iodide ions (X = I(-)). The Raman spectrum of solid sodium iodide dihydrate supports the assignment of the latter. Using the MCR results, relative Raman scattering cross sections of 4.0 ± 0.6 and 14.0 ± 0.1 were calculated for the mono- and di-iodide species, respectively (compared to that of bulk water set to unity). In addition, it was found that at relatively low concentrations each iodide ion affects the Raman spectrum of roughly 22 surrounding water molecules, indicating that the influence of iodide extends beyond the first solvation shell. Our results demonstrate that the Raman bending vibration of water is a sensitive probe, providing new insights into anion solvation in aqueous environments.

Water Exchange at a Hydrated Platinum Electrode is Rare and Collective

We use molecular dynamics simulations to study the exchange kinetics of water molecules at a model metal electrode surface–exchange between water molecules in the bulk liquid and water molecules bound to the metal. This process is a rare event, with a mean residence time of a bound water of ∼40 ns for the model we consider. With analysis borrowed from the techniques of rare-event sampling, we show how this exchange or desorption is controlled by (1) reorganization of the hydrogen bond network within the adlayer of bound water molecules and (2) interfacial density fluctuations of the bulk liquid adjacent to the adlayer. We define collective coordinates that describe the desorption mechanism. Spatial and temporal correlations associated with a single event extend over nanometers and tens of picoseconds.

Solvation thermodynamics of amino acid side chains on a short peptide backbone.

The hydration process of side chain analogue molecules differs from that of the actual amino acid side chains in peptides and proteins owing to the effects of the peptide backbone on the aqueous solvent environment. A recent molecular simulation study has provided evidence that all nonpolar side chains, attached to a short peptide backbone, are considerably less hydrophobic than the free side chain analogue molecules. In contrast to this, the hydrophilicity of the polar side chains is hardly affected by the backbone. To analyze the origin of these observations, we here present a molecular simulation study on temperature dependent solvation free energies of nonpolar and polar side chains attached to a short peptide backbone. The estimated solvation entropies and enthalpies of the various amino acid side chains are compared with existing side chain analogue data. The solvation entropies and enthalpies of the polar side chains are negative, but in absolute magnitude smaller compared with the corresponding analogue data. The observed differences are large; however, owing to a nearly perfect enthalpy-entropy compensation, the solvation free energies of polar side chains remain largely unaffected by the peptide backbone. We find that a similar compensation does not apply to the nonpolar side chains; while the backbone greatly reduces the unfavorable solvation entropies, the solvation enthalpies are either more favorable or only marginally affected. This results in a very small unfavorable free energy cost, or even free energy gain, of solvating the nonpolar side chains in strong contrast to solvation of small hydrophobic or nonpolar molecules in bulk water. The solvation free energies of nonpolar side chains have been furthermore decomposed into a repulsive cavity formation contribution and an attractive dispersion free energy contribution. We find that cavity formation next to the peptide backbone is entropically favored over formation of similar sized nonpolar side chain cavities in bulk water, in agreement with earlier work in the literature on analysis of cavity fluctuations at nonpolar molecular surfaces. The cavity and dispersion interaction contributions correlate quite well with the solvent accessible surface area of the nonpolar side chains attached to the backbone. This correlation however is weak for the overall solvation free energies owing to the fact that the cavity and dispersion free energy contributions are almost exactly cancelling each other.

Mechanism of Meropenem Hydrolysis by New Delhi Metallo β-Lactamase

New Delhi metallo β-lactamase (NDM-1) is a recent addition to the metallo-β-lactamases family that is capable of hydrolyzing most of the available antibiotics, including the new generation carbapenems. Here, we report the mechanism of Meropenem hydrolysis catalyzed by NDM-1 based on hybrid quantum-mechanical/molecular-mechanical metadynamics simulations. Our work elicits the molecular details of the catalytic mechanism and free energy profiles along the reaction pathway. We identified the ring opening step involving the nucleophilic addition of the bridging hydroxyl group on the β-lactam ring of the drug as the rate-determining step. Subsequent protonation of β-lactam nitrogen occurs from a bulk water molecule that diffuses into the active site and is preferred over proton transfer from the bridging hydroxyl group or from the protonated Asp124. The roles of important active site residues of NDM-1 and change in the coordination environment of Zn ions during the hydrolysis are also scrutinized.

Crystal structure of Legionella pneumophila dephospho-CoA kinase reveals a non-canonical conformation of P-loop.

Dephospho-CoA kinase (DPCK; EC 2.7.1.24) catalyzes the final step in the coenzyme A biosynthetic pathway. DPCK transfers a phosphate group from ATP to the 3-hydroxyl group of the ribose of dephosphocoenzyme A (dCoA) to yield CoA and ADP. Upon the binding of ligands, large conformational changes is induced in DPCKs, as well as in many other kinases, to shield the bound ATP in their catalytic site from the futile hydrolysis by bulk water molecules. To investigate the molecular mechanisms underlying the phosphoryl transfer during DPCK catalytic cycle, we determined the crystal structures of the Legionellapneumophila DPCK (LpDPCK) both in its apo-form and in complex with ATP. The structures reveal that LpDPCK comprises of three domains, the classical core domain, the CoA domain, and the LID domain, which are packed together to create a central cavity for substrate-binding and enzymatic catalysis. The binding of ATP induces large conformational changes, including a hinge-bending motion of the CoA binding domain and the "helix to loop" conformational change of the P-loop. Finally, modeling of a dCoA molecule to the enzyme provides insights into the catalytic mechanism of DPCK.

Vibrational dephasing and frequency shifts of hydrogen-bonded pyridine–water complexes

In this paper we present the picosecond vibrational dynamics and Raman shifts of hydrogen-bonded pyridine–water complexes present in aqueous solutions in a wide concentration range from dense to extreme dilute solutions. We studied the vibrational dephasing and vibrational frequency modulation by calculating time correlation functions of vibrational relaxation by fits in the frequency domain. The concentration induced variations in bandwidths, band frequencies and characteristic dephasing times have been estimated and interpreted as effects due to solute–solvent interactions. The time-correlation functions of vibrational dephasing were obtained for the ring breathing mode of both “free” and hydrogen-bonded pyridine molecules and it was found that sufficiently deviate from the Kubo model. There is a general agreement in the whole concentration range with the modeling proposed by the Rothschild approach, which applies to complex liquids. The results have shown that the reorientation of pyridine aqueous solutions is very slow and hence in both scattering geometries only vibrational dephasing is probed. It is proposed that the spectral changes depend on the perturbations induced by the dynamics of the water molecules in the first hydration cell and water in bulk, while at extreme dilution conditions, the number of bulk water molecules increases and the interchange between molecules belonging to the first hydration cell may not be the predominant modulation mechanism. The evolution of several parameters, such as the characteristic times, the percentage of Gaussian character in the peak shape and the a parameter are indicative of drastic variations at extreme dilution revealing changes in the vibrational relaxation of the pyridine complexes in the aqueous environment. The higher dilution is correlated to diffusion of water molecules into the reference pyridine system in agreement with the jump diffusion model, while at extreme dilutions, almost all pyridine molecules are elaborated in hydrogen bonding. The results are discussed in the framework of the current phenomenological status of the field.

Hydrodynamic and nonhydrodynamic contributions to the bimolecular collision rates of solute molecules in supercooled bulk water.

Bimolecular collision rate constants of a model solute are measured in water at T = 259–303 K, a range encompassing both normal and supercooled water. A stable, spherical nitroxide spin probe, perdeuterated 2,2,6,6-tetramethyl-4-oxopiperidine-1-oxyl, is studied using electron paramagnetic resonance spectroscopy (EPR), taking advantage of the fact that the rotational correlation time, τR, the mean time between successive spin exchanges within a cage, τRE, and the long-time-averaged spin exchange rate constants, Kex, of the same solute molecule may be measured independently. Thus, long- and short-time translational diffusion behavior may be inferred from Kex and τRE, respectively. In order to measure Kex, the effects of dipole–dipole interactions (DD) on the EPR spectra must be separated, yielding as a bonus the DD broadening rate constants that are related to the dephasing rate constant due to DD, Wdd. We find that both Kex and Wdd behave hydrodynamically; that is to say they vary monotonically with T/η or η/T, respectively, where η is the shear viscosity, as predicted by the Stokes–Einstein equation. The same is true of the self-diffusion of water. In contrast, τRE does not follow hydrodynamic behavior, varying rather as a linear function of the density reaching a maximum at 276 ± 2 K near where water displays a maximum density.

Molecular Dynamics Study of Conserved Water Molecules in PrPc and Pathological Point Mutation T188R

The secondary conformational change of cellular prion protein (PrPc) into scrapie prion protein (PrPsc) is critical important for prion disease including bovine spongiform encephalopathy and Creutzfeldt-Jakob disease.This study performed molecular dynamics simulations of both human prion protein (WT-PrPc) and pathogenic mutant prion protein (T188R) and analyze its conserved water molecules in equilibrium state.Although the secondary structure is not dramatically different, increase of Cα- RMSF occurs not only around mutation point but also in a wide range. This increase of intensity may result from the instability caused by point mutation.Next, we calculated mean residence time for each protein atoms to find conserved water molecules. Residence time is defined that the time from when the water molecule get into the sphere whose center is a protein atom to the evacuation of water molecule from the sphere. We found four conserved water molecule sites in WT-PrPc. Conserved water molecules also keep the higher rotational correlation than water molecules in bulk. These results shows that conserved water molecules are not only wrapped but also fixed by hydrogen bond between conserved water molecule and amino acids.On the other hand, every residence time of T188R decrease compared by WT-PrPc and the rotational relaxation is also faster than WT-PrPc although secondary structure is almost not different. These results show that conserved water molecules cannot survive related to the wide range fluctuation caused by point mutation.

Proton Transfer in Wild-Type GFP and S205V Mutant Is Reduced by Conformational Changes of Residues in the Proton Wire

To study the dynamics and mechanisms of the proton wire of wild type green fluorescence protein (wt-GFP) and its S205V mutant, we applied molecular dynamics (MD) simulations and compared the results with the X-ray structures of both proteins and with the proton transfer kinetics of these proteins studied by the time-resolved emission technique. The MD simulations for the wt-GFP show that, in the proton wire, the Ser205 exists in two conformations with similar probabilities. One conformation supports the proton transfer, and the other does not. The fluctuation between the two conformers is relatively slow. This result may explain the time-resolved emission spectrum's long-time fluorescence tail of the wt-GFP chromophore's protonated form. The MD simulations of the S205V mutant show that the water molecule in the proton wire is replaced by other bulk water molecules along the simulations of 60 ns. Furthermore, as in the wt-GFP, the Thr203 also exists in two conformations in which only one conformation supports the proton transfer. These two findings give an insight into the relatively slow proton transfer rate in the S205V mutant in comparison to the wt-GFP.

Role of Water in Mukaiyama–Aldol Reaction Catalyzed by Lanthanide Lewis Acid: A Computational Study

Carbon-carbon bond formations, such as Kobayashi modification of the Mukaiyama-Aldol reaction, catalyzed by lanthanide (Ln) Lewis acid in aqueous solution comprise one of the most attractive types of reactions in terms of green chemistry. However, their detailed mechanisms and the role of water molecules remained unclear. In order to explore complex potential energy surfaces for the water and substrate coordination around Eu(3+) as well as the detailed mechanism of the Mukaiyama-Aldol reaction between trimethylsilyl (TMS) cylcohexenolate and benzaldehyde (BA) catalyzed by Eu(3+), the recently developed anharmonic downward distortion following (ADDF) and artificial force-induced reaction (AFIR) methods were used with the B3LYP-D3 theory. The most favorable water coordination structures are Eu(3+)(H2O)8 and Eu(3+)(H2O)9; they are comparable in free energy and are likely to coexist, with an effective coordination number of 8.3. Eu(3+)(H2O)8(BA) is the best aldehyde coordinated structure. Starting with this complex, the Mukaiyama-Aldol reaction proceeds via a stepwise mechanism, first C-C bond formation between the substrates, followed by proton transfer from water to BA and then TMS dissociation caused by nucleophilic attack by bulk water molecules. Why did the yield of the Mukaiyama-Aldol reaction catalyzed by Ln(3+) in organic solvent dramatically increase upon addition of water? Without water, the reverse reaction (C-C cleavage) takes place easily. Why did this reaction show syn-preference in water? The anti transition state for C-C formation in water is entropically less favored relative to the syn transition state because of the existence of a rigid hydrogen bond between the TMS part and coordination water around Eu(3+) in the former.

Nonadiabatic Dynamics of Positive Charge during Photocatalytic Water Splitting on GaN(10-10) Surface: Charge Localization Governs Splitting Efficiency

Photochemical water splitting is a promising avenue to sustainable, clean energy and fuel production. Gallium nitride (GaN) and its solid solutions are excellent photocatalytic materials; however, the efficiency of the process is low on pure GaN, and cocatalysts are required to increase the yields. We present the first time-domain theoretical study of the initial steps of photocatalytic water splitting on a GaN surface. Our state-of-the-art simulation technique, combining nonadiabatic molecular dynamics and time-dependent density functional theory, allows us to characterize the mechanisms and time scales of the evolution of the photogenerated positive charge (hole) and the subsequent proton transfer at the GaN/water interface. The calculations show that the hole loses its excess energy within 100 fs and localizes primarily on the nitrogen atoms of the GaN surface, initiating a sequence of proton-transfer events from the surface N-H group to the nearby OH groups and bulk water molecules. Water splitting requires hole localization on oxygen rather than nitrogen, necessitating nonadiabatic transitions uphill in energy on pure GaN. Such transitions happen rarely, resulting in low yields of the photocatalytic water splitting observed experimentally. We conclude that efficient cocatalysts should favor localization of the photogenerated hole on oxygen-containing species at the semiconductor/water interface.