Search

Interconversion of the oxidation states of uranium enables separations and reactivity schemes involving this element and contributes to technologies for recycling of spent nuclear fuels. The redox behaviors of uranium species impact these processes, but use of electrochemical methods to drive reactions of molecular uranium complexes and to obtain molecular insights into the outcomes of electrode-driven reactions has received far less attention than it deserves. Here, we show that electro-reduction of the uranyl ion (UO 2 2+ ) can be used to promote stepwise functionalization of the typically unreactive oxo groups with exogenous triphenylborane (BPh 3 ) serving as a moderate electrophile, avoiding the conventional requirement for a chemical reductant. Parallel electroanalytical, spectrochemical, and chemical reactivity studies, supported by spectroscopic findings and structural data from X-ray diffraction analysis on key reduced and borylated products, demonstrate that our electrochemical approach largely avoids undesired cross reactions and disproportionation pathways; these usually impact the multicomponent systems needed for uranyl functionalization chemistry. Joint computational studies have been used to map the changes associated with U-O activation and to quantify the free energy differences related to key reactions. Taken together, the results suggest that electrochemical methods can be used for selective interconversion of molecular actinide species, reminiscent of methods commonly employed in transition metal redox catalysis.

Ensemble -average sampling of structures from ab initio molecular dynamics (AIMD) simulations can be used to predict theoretical extended X-ray absorption fine structure (EXAFS) signals that closely match experimental spectra. However, AIMD simulations are time-consuming and resource-intensive, particularly for solvated lanthanide ions, which often form multiple nonrigid geometries with high coordination numbers. To accelerate the characterization of lanthanide structures in solution, we employed the Northwest Potential Energy Surface Search Engine (NWPEsSe), an adaptive-learning global optimization algorithm, to efficiently screen first-shell structures. As case studies, we examine two systems: Eu(NO 3 ) 3 dissolved in acetonitrile with a terpyridine ligand (terpyNO 2 ), and Nd(NO 3 ) 3 dissolved in acetonitrile. The theoretical spectra for structures identified by NWPEsSe were compared to both experimental and AIMD-derived EXAFS spectra. The NWPEsSe algorithm successfully identified the proper solvation structure for both Eu(NO 3 ) 3 (terpyNO 2 ) and Nd(NO 3 )(acetonitrile) 3 , with the calculated EXAFS signals closely matching the experimental spectra for the Eu-ligand complex and showing good similarity for the Nd salt; the better agreement with the ligand-containing structure is attributed to a less dynamic coordination environment due to the rigid ligand. The key advantage of the global optimization algorithm lies in its ability to sample the coordination environment across the potential energy surface and reduce the time required to identify structures from generally a month to within a week. Additionally, this approach is versatile and can be adapted to characterize main-group metal complexes.

Although crown ethers can selectively bind many metal cations, little is known regarding the solution properties of crown ether complexes of the uranyl dication, UO 2 2+ . Here, the synthesis and characterization of isolable complexes in which the uranyl dication is bound in an 18-crown-6-like moiety are reported. A tailored macrocyclic ligand, templated with a Pt(II) center, captures UO 2 2+ in the crown moiety, as demonstrated by results from single-crystal X-ray diffraction analysis. The U(V) oxidation state becomes accessible at a quite positive potential ( E 1/2 ) of -0.18 V vs Fc +/0 upon complexation, representing the most positive U VI /U V potential yet reported for the UO 2 n + core. Isolation and characterization of the U(V) form of the crown complex are also reported here; there are no prior reports of reduced uranyl crown ether complexes, but U(V) is clearly stabilized by crown chelation. Joint computational studies show that the electronic structure of the U(V) form results in significant weakening of U-O oxo bonding despite the quite positive reduction potential at which this species can be accessed, underscoring that crown-ligated uranyl species could demonstrate unique reactivity under only modestly reducing conditions.

Single-site copper-based catalysts have shown remarkable activity and selectivity for a variety of reactions. However, deactivation by sintering in high-temperature reducing environments remains a challenge and often limits their use due to irreversible structural changes to the catalyst. Here, we report zeolite-based copper catalysts in which copper oxide agglomerates formed after reaction can be repeatedly redispersed back to single sites using an oxidative treatment in air at 550 °C. Under different environments, single-site copper in Cu-Zn-Y/deAlBeta undergoes dynamic changes in structure and oxidation state that can be tuned to promote the formation of key active sites while minimizing deactivation through Cu sintering. For example, single-site Cu 2+ reduces to Cu 1+ after catalyst pretreatment (270 °C, 101 kPa H 2 ) and further to Cu 0 nanoparticles under reaction conditions (270-350 °C, 7 kPa EtOH, 94 kPa H 2 ) or accelerated aging (400-450 °C, 101 kPa H 2 ). After regeneration at 550 °C in air, agglomerated CuO was dispersed back to single sites in the presence and absence of Zn and Y, which was verified by imaging, in situ spectroscopy, and catalytic rate measurements. Ab initio molecular dynamics simulations show that solvation of CuO monomers by water facilitates their transport through the zeolite pore, and condensation of the CuO monomer with a fully protonated silanol nest entraps copper and reforms the single-site structure. The capability of silanol nests to trap and stabilize copper single sites under oxidizing conditions could extend the use of single-site copper catalysts to a wider variety of reactions and allows for a simple regeneration strategy for copper single-site catalysts.

Solvation and ion valency effects on selectivity of metal oxyanions at redox-polymer interfaces are explored through in situ spatial-temporally resolved neutron reflectometry combined with large scale ab initio molecular dynamics. The selectivity of ReO 4 - vs MoO 4 2- for two redox-metallopolymers, poly(vinyl ferrocene) (PVFc) and poly(3-ferrocenylpropyl methacrylamide) (PFPMAm) is evaluated. PVFc has a higher Re/Mo separation factor compared to PFPMAm at 0.6 V vs Ag/AgCl. In situ techniques show that both PVFc and PFPMAm swell in the presence of ReO 4 - (having higher solvation with PFPMAm), but do not swell in contact with MoO 4 2- . Ab initio molecular simulations suggest that MoO 4 2- maintains a well-defined double solvation shell compared to ReO 4 - . The more loosely solvated anion (ReO 4 - ) is preferably adsorbed by the more hydrophobic redox polymer (PVFc), and electrostatic cross-linking driven by divalent anionic interactions could impair film swelling. Thus, the in-depth understanding of selectivity mechanisms can accelerate the design of ion-selective redox-mediated separation systems for transition metal recovery and recycling., Competing Interests: The authors declare no competing financial interest., (© 2024 The Authors. Published by American Chemical Society.)

Proton transfer is critically important to many electrocatalytic reactions, and directed proton delivery could open new avenues for the design of electrocatalysts. However, although this approach has been successful in molecular electrocatalysis, proton transfer has not received the same attention in heterogeneous electrocatalyst design. Here, we report that a metal oxide proton relay can be built within heterogeneous electrocatalyst architectures and improves the kinetics of electrochemical hydrogen evolution and oxidation reactions. The volcano-type relationship between activity enhancement and p K a of amine additives confirms this improvement; we observe maximum rate enhancement when the p K a of a proton relay matches the pH of the electrolyte solution. Density-functional-theory-based reactivity studies reveal a decreased proton transfer energy barrier with a metal oxide proton relay. These findings demonstrate the possibility of controlling the proton delivery and enhancing the reaction kinetics by tuning the chemical properties and structures at heterogeneous interfaces.

Cement and concrete are vital materials used to construct durable habitats and infrastructure that withstand natural and human-caused disasters. Still, concrete cracking imposes enormous repair costs on societies, and excessive cement consumption for repairs contributes to climate change. Therefore, the need for more durable cementitious materials, such as those with self-healing capabilities, has become more urgent. In this review, we present the functioning mechanisms of five different strategies for implementing self-healing capability into cement based materials: (1) autogenous self-healing from ordinary portland cement and supplementary cementitious materials and geopolymers in which defects and cracks are repaired through intrinsic carbonation and crystallization; (2) autonomous self-healing by (a) biomineralization wherein bacteria within the cement produce carbonates, silicates, or phosphates to heal damage, (b) polymer-cement composites in which autonomous self-healing occurs both within the polymer and at the polymer-cement interface, and (c) fibers that inhibit crack propagation, thus allowing autogenous healing mechanisms to be more effective. In all cases, we discuss the self-healing agent and synthesize the state of knowledge on the self-healing mechanism(s). In this review article, the state of computational modeling across nano- to macroscales developed based on experimental data is presented for each self-healing approach. We conclude the review by noting that, although autogenous reactions help repair small cracks, the most fruitful opportunities lay within design strategies for additional components that can migrate into cracks and initiate chemistries that retard crack propagation and generate repair of the cement matrix.

In situ Al K-edge X-ray absorption near edge structure (XANES) and Extended X-ray absorption fine structure (EXAFS) spectroscopy in conjunction with ab initio molecular dynamics (AIMD) simulations show that adsorption of 1-propanol alters the structure of the Brønsted acid site through changes in the associated aluminum-oxygen tetrahedron in zeolite H-MFI. The decreasing intensity of the pre-edge signal of the in situ Al K-edge XANES spectra with increasing 1-propanol coverage shows that Al T-sites become more symmetric as the sorbed alcohol molecules form monomers, dimers, and trimers. The adsorption of monomeric 1-propanol on Brønsted acid sites reduces the distortion of the associated Al T-site, shortens the Al-O distance, and causes the formation of a Zundel-like structure. With dimeric and trimeric alcohol clusters, the zeolite proton is fully transferred to the alcohols and the aluminum-oxygen tetrahedron becomes fully symmetric. The subtle changes in Al-K-edge XANES in the presence of sorbate structures, with the use of theory, are used to probe the local zeolite structures and provide a basis to predict the population and chemical state of the sorbed species., Competing Interests: The authors declare no competing financial interest., (© 2023 The Authors. Published by American Chemical Society.)

Efficient direct air capture (DAC) of CO 2 will require strategies to deal with the relatively low concentration in the atmosphere. One such strategy is to employ the combination of a CO 2 -selective membrane coupled with a CO 2 capture solvent acting as a draw solution. Here, the interactions between a leading water-lean carbon-capture solvent, a polyether ether ketone (PEEK)-ionene membrane, CO 2 , and combinations were probed using advanced NMR techniques coupled with advanced simulations. We identify the speciation and dynamics of the solvent, membrane, and CO 2 , presenting spectroscopic evidence of CO 2 diffusion through benzylic regions within the PEEK-ionene membrane, not spaces in the ionic lattice as expected. Our results demonstrate that water-lean capture solvents provide a thermodynamic and kinetic funnel to draw CO 2 from the air through the membrane and into the bulk solvent, thus enhancing the performance of the membrane. The reaction between the carbon-capture solvent and CO 2 produces carbamic acid, disrupting interactions between the imidazolium (Im + ) cations and the bistriflimide anions within the PEEK-ionene membrane, thereby creating structural changes through which CO 2 can diffuse more readily. Consequently, this restructuring results in CO 2 diffusion at the interface that is faster than CO 2 diffusion in the bulk carbon-capture solvent., (© 2023 Battelle Memorial Institute, Oak Ridge National Laboratory and The Authors. ChemSusChem published by Chemistry Europe and Wiley-VCH GmbH.)

Research interest in single-atom catalysts (SACs) has been continuously increasing. However, the lack of understanding of the dynamic behaviors of SACs during applications hinders catalyst development and mechanistic understanding. Herein, we report on the evolution of active sites over Pd/TiO 2 -anatase SAC (Pd 1 /TiO 2 ) in the reverse water-gas shift (rWGS) reaction. Combining kinetics, in situ characterization, and theory, we show that at T ≥ 350 °C, the reduction of TiO 2 by H 2 alters the coordination environment of Pd, creating Pd sites with partially cleaved Pd-O interfacial bonds and a unique electronic structure that exhibit high intrinsic rWGS activity through the carboxyl pathway. The activation by H 2 is accompanied by the partial sintering of single Pd atoms (Pd 1 ) into disordered, flat, ∼1 nm diameter clusters (Pd n ). The highly active Pd sites in the new coordination environment under H 2 are eliminated by oxidation, which, when performed at a high temperature, also redisperses Pd n and facilitates the reduction of TiO 2 . In contrast, Pd 1 sinters into crystalline, ∼5 nm particles (Pd NP ) during CO treatment, deactivating Pd 1 /TiO 2 . During the rWGS reaction, the two Pd evolution pathways coexist. The activation by H 2 dominates, leading to the increasing rate with time-on-stream, and steady-state Pd active sites similar to the ones formed under H 2 . This work demonstrates how the coordination environment and nuclearity of metal sites on a SAC evolve during catalysis and pretreatments and how their activity is modulated by these behaviors. These insights on SAC dynamics and the structure-function relationship are valuable to mechanistic understanding and catalyst design.

Polyoxometalates (POMs) with localized radical or open-shell metal sites have the potential to be used as transformative electronic spin based molecular qubits (MQs) for quantum computing (QC). For practical applications, MQs have to be immobilized in electronically or optically addressable arrays which introduces interactions with supports as well as neighboring POMs. Herein, we synthesized Keggin POMs with both tungsten (W) and vanadium (V) addenda atoms. Ion soft landing, a highly-controlled surface modification technique, was used to deliver mass-selected V-doped POMs to different self-assembled monolayer surfaces on gold (SAMs) without the solvent, counterions, and contaminants that normally accompany deposition from solution. Alkylthiol, perfluorinated, and carboxylic-acid terminated monolayers were employed as representative model supports on which different POM-surface and POM-POM interactions were characterized. We obtained insights into the vibrational properties of supported V-doped POMs and how they are perturbed by interactions with specific surface functional groups using infrared reflection absorption and scattering-type scanning near-field optical microscopy, as well as tip enhanced Raman spectroscopy. Different functional groups on SAMs and nanoscale heterogeneity are both shown to modulate the observed spectroscopic signatures. Spectral shifts are also found to be dependent on POM-POM interactions. The electronic structure of the V-doped POMs was determined in the gas phase using negative ion photoelectron spectroscopy and on surfaces with scanning Kelvin probe microscopy. The chemical functionality and charge transfer properties of the SAMs are demonstrated to exert an influence on the charge state and electronic configuration of supported V-doped POMs. The geometric and electronic structure of the POMs were also calculated using density functional theory. Our joint experimental and theoretical findings provide insight into how V substitution as well as POM-surface and POM-POM interactions influence the vibrational properties of POMs.

Absolute thermodynamic quantities for critical chemical reactions are needed to determine the role of solvents and reactive environments in catalysis and electrocatalysis. Theoretical methods can provide such quantification but are often hindered by the innate complexity of electron correlation and dynamic relaxation of solvent environments. We present and validate a protocol for calculating the redox potentials of the ferrocene/ferrocenium redox pair in acetonitrile. Equation-of-motion and effective fragment potential (EFP) methods are used to characterize the adiabatic and vertical ionization potentials as well as the electron affinity processes. We benchmark molecular mechanics against the EFP model to show the differences in the ferrocene electronic polarizability in two redox states. Our best estimate of the redox potential (4.94 eV) agrees well with the experimental value (4.93 eV). This demonstrates the ability of modern computational methods to predict absolute redox potentials quantitatively and to quantify the correlation of dynamic effects, which underlie their origin.

A mechanistic understanding of metal-organic framework (MOF) synthesis and scale-up remains underexplored due to the complex nature of the interactions of their building blocks. In this work, we investigate the collective assembly of building units at the early stages of MOF nucleation, using MIL-101(Cr) as a prototypical example. Using large-scale molecular dynamics simulations, we observe that the choice of solvent (water and N , N -dimethylformamide), the introduction of ions (Na + and F - ) and the relative populations of MIL-101(Cr) half-secondary building unit (half-SBU) isomers have a strong influence on the cluster formation process. Additionally, the shape, size, nucleation and growth rates, crystallinity, and short and long-range order largely vary depending on the synthesis conditions. We evaluate these properties as they naturally emerge when interpreting the self-assembly of MOF nuclei as the time evolution of an undirected graph. Solution-induced conformational complexity and ionic concentration have a dramatic effect on the morphology of clusters emerging during assembly. While pure solvents lead to the rapid formation of a small number of large clusters, the presence of ions in aqueous solutions results in smaller clusters and slower nucleation. This diversity is captured by the key features of the graph representation. Principle component analysis on graph properties reveals that only a small number of molecular descriptors is needed to deconvolute MOF self-assembly. Descriptors such as the average coordination number between half-SBUs and fractal dimension are of particulalr interest as they can be can be followed experimentally by techniques like by time-resolved spectroscopy. Ultimately, graph theory emerges as an approach that can be used to understand complex processes revealing molecular descriptors accessible by both simulation and experiment.

Functionalization of graphene oxide (GO) membranes is generally achieved using carboxyl groups as binding sites for ligands. Herein, by taking advantage of the ability of imidazolium-based ionic liquids (ILs) to undergo an epoxide ring-opening reaction, a new approach of GO modification was established, in which ILs were bonded to the abundant epoxides on GO without sacrificing the carboxyl groups. Computational methods confirmed this unique configuration of ILs on GO, which enabled the dispersion of IL/GO flakes in water for facile casting into laminate membranes. Compared with neat GO, the ILs in IL/GO membranes served as spacers that substantially reduced the multi-valent cation mobility, simultaneously facilitated ion desolvation, and increased the water flux across the membrane. Our studies found that the higher separation efficiency of IL/GO membranes may be attributed to the synergistic modification of the hydrophobicity and surface charge. Specifically, the protonated nitrogen of the imidazolium cations altered the surface charge of GO, thereby generating electrostatic repulsion that enhanced the selectivity of cation rejection. On the other hand, the increased length of the alkyl chains bound to the imidazolium rings was found to increase the hydrophobicity of GO, which, in turn, aided the fine-tuning of the water desolvation/transport dynamics at the GO/IL interface to achieve a high water flux. Additionally, the water retention was reduced on the hydrophobic planes, which inhibited GO swelling during aqueous separations. Molecular dynamics simulations revealed increased water diffusivity when ILs were intercalated within GO layers. We establish that without requiring a high energy input, functionalization of GO membranes with ILs may be a promising approach to achieve efficient ion separation and critical material recovery.

Increasing atmospheric concentrations of greenhouse gases due to industrial activity have led to concerning levels of global warming. Reducing carbon dioxide (CO 2 ) emissions, one of the main contributors to the greenhouse effect, is key to mitigating further warming and its negative effects on the planet. CO 2 capture solvent systems are currently the only available technology deployable at scales commensurate with industrial processes. Nonetheless, designing these solvents for a given application is a daunting task requiring the optimization of both thermodynamic and transport properties. Here, we discuss the use of atomic scale modeling for computing reaction energetics and transport properties of these chemically complex solvents. Theoretical studies have shown that in many cases, one is dealing with a rich ensemble of chemical species in a coupled equilibrium that is often difficult to characterize and quantify by experiment alone. As a result, solvent design is a balancing act between multiple parameters which have optimal zones of effectiveness depending on the operating conditions of the application. Simulation of reaction mechanisms has shown that CO 2 binding and proton transfer reactions create chemical equilibrium between multiple species and that the agglomeration of resulting ions and zwitterions can have profound effects on bulk solvent properties such as viscosity. This is balanced against the solvent systems needing to perform different functions (e.g., CO 2 uptake and release) depending on the thermodynamic conditions (e.g., temperature and pressure swings). The latter constraint imposes a "Goldilocks" range of effective parameters, such as binding enthalpy and p K a , which need to be tuned at the molecular level. The resulting picture is that solvent development requires an integrated approach where theory and simulation can provide the necessary ingredients to balance competing factors., Competing Interests: The authors declare no competing financial interest., (© 2022 The Authors. Published by American Chemical Society.)

Adsorption of organics in the aqueous phase is an area which is experimentally difficult to measure, while computational techniques require extensive configurational sampling of the solvent and adsorbate. This is exceedingly computationally demanding, which excludes its routine use. If implicit solvent could be applied instead, this would dramatically reduce the computational cost as configurational sampling of solvent is not needed. Here, using statistical thermodynamic arguments and DFT calculations with implicit solvent models, we show that semiquantitative values for the free energy and entropy change of adsorption in the aqueous phase (Δ G ads solv and Δ S ads solv ) for small organics can be calculated, for a range of coverages. We parametrize the soft sphere based solute dielectric cavity to an approximated free energy of solvation for a single Pt atom at the (111) facet, forming upper and lower bounds based on the entropy of water at the aqueous metal interface (Δ G solv (Pt) = -4.35 to -7.18 kJ mol -1 ). This captures the decrease in Δ G ads solv compared to the free energy of adsorption in the vacuum phase (Δ G ads vac ), while solvent models with electron density based cavities fail to do so. For a range of oxygenated aromatics, the adsorption energetics using horizontal gas phase geometries significantly overestimate Δ G ads solv compared to experiment by ∼100 kJ mol -1 , but they agree with ab initio MD simulations using similar geometries. This suggests oxygenated aromatic compounds adsorb perpendicular to the metallic surface, while the Δ G ads solv for vertical geometries of furfural and cyclohexanol agree to within 20 kJ mol -1 of experimental studies. The proposed techniques provide an inexpensive toolset for validation and prediction of adsorption energetics on solvated metallic surfaces, which could be further validated by the future availability of more experimental measurements for the aqueous entropy/free energy of adsorption.

A combined experimental and theoretical study has been carried out on the wetting and reactivity of water-lean carbon capture solvents on the surface of common column packing materials. Paradoxically, these solvents are found to be equally able to wet hydrophobic and hydrophilic surfaces. The solvents are amphiphilic and can adapt to any interfacial environment, owing to their inherent heterogeneous (nonionic/ionic) molecular structure. Ab initio molecular dynamics indicates that these structures enable the formation of a strong adlayer on the surface of hydrophilic surfaces like oxidized steel which promotes solvent decomposition akin to hydrolysis from surface oxides and hydroxides. This decomposition passivates the surface, making it effectively hydrophobic, and the decomposed solvent promotes leaching of the iron into the bulk fluid. This study links the wetting behavior to the observed corrosion of the steels by decomposition of solvent at steel interfaces. The overall affect is strongly dependent on the chemical composition of the solvent in that amines are stable, whereas imines and alcohols are not. Moreover, plastic packing shows little to no solvent degradation, but an equal degree of wetting., (© 2021 Wiley-VCH GmbH.)

Actinide molten salts represent a class of important materials in nuclear energy. Understanding them at a molecular level is critical for the proper and optimal design of relevant technological applications. Yet, owing to the complexity of electronic structure due to the 5f orbitals, computational studies of heavy elements in condensed phases using ab initio potentials to study the structure and dynamics of these elements embedded in molten salts are difficult. This lack of efficient computational protocols makes it difficult to obtain information on properties that require extensive statistical sampling like transport properties. To tackle this problem, we adopted a machine-learning approach to study ThCl 4 -NaCl and UCl 3 -NaCl binary systems. The machine-learning potential with the density functional theory accuracy allows us to obtain long molecular dynamics trajectories (ns) for large systems (10 3 atoms) at a considerably low computing cost, thereby efficiently gaining information about their bonding structures, thermodynamics, and dynamics at a range of temperatures. We observed a considerable change in the coordination environments of actinide elements and their characteristic coordination sphere lifetime. Our study also suggests that actinides in molten salts may not follow well-known entropy-scaling laws.

The analogy between single-atom catalysts (SACs) and molecular catalysts predicts that the specific catalytic activity of these systems is constant. We provide evidence that this prediction is not necessarily true. As a case in point, we show that the specific activity over ceria-supported single Pd atoms linearly increases with metal atom density, originating from the cumulative enhancement of CeO 2 reducibility. The long-range electrostatic footprints (≈1.5 nm) around each Pd site overlap with each other as surface Pd density increases, resulting in an observed deviation from constant specific activity. These cooperative effects exhaust previously active O atoms above a certain Pd density, leading to their permanent removal and a consequent drop in reaction rate. The findings of our combined experimental and computational study show that the specific catalytic activity of reducible oxide-supported single-atom catalysts can be tuned by varying the surface density of single metal atoms., (© 2021 Wiley-VCH GmbH.)

The catalytic reduction in carbon dioxide is a crucial step in many chemical industrial reactions, such as methanol synthesis, the reverse water-gas shift reaction, and formic acid synthesis. Here, we investigate the role of bulk hydrogen, where hydrogen atoms are found deep inside a metal surface as opposed to subsurface ones, upon CO 2 reduction over a Ni(110) surface using density functional theory and ab initio molecular dynamics simulations. While it has previously been shown that subsurface hydrogen stabilizes CO 2 and can aid in overcoming reaction barriers, the role of bulk hydrogen is less studied and thus unknown with regard to CO 2 reduction. We find that the presence of bulk hydrogen can significantly alter the electronic structure of the Ni(110) surface, particularly the work function and d-band center, such that CO 2 adsorbs more strongly to the surface and is more easily reduced. Our results show an enhanced CO 2 dissociation in the presence of bulk hydrogen, shedding light on a hitherto underappreciated mechanistic pathway for CO 2 reduction on metal surfaces.

Cu-zeolites are able to directly convert methane to methanol via a three-step process using O 2 as oxidant. Among the different zeolite topologies, Cu-exchanged mordenite (MOR) shows the highest methanol yields, attributed to a preferential formation of active Cu-oxo species in its 8-MR pores. The presence of extra-framework or partially detached Al species entrained in the micropores of MOR leads to the formation of nearly homotopic redox active Cu-Al-oxo nanoclusters with the ability to activate CH 4 . Studies of the activity of these sites together with characterization by 27 Al NMR and IR spectroscopy leads to the conclusion that the active species are located in the 8-MR side pockets of MOR, and it consists of two Cu ions and one Al linked by O. This Cu-Al-oxo cluster shows an activity per Cu in methane oxidation significantly higher than of any previously reported active Cu-oxo species. In order to determine unambiguously the structure of the active Cu-Al-oxo cluster, we combine experimental XANES of Cu K- and L-edges, Cu K-edge HERFD-XANES, and Cu K-edge EXAFS with TDDFT and AIMD-assisted simulations. Our results provide evidence of a [Cu 2 AlO 3 ] 2+ cluster exchanged on MOR Al pairs that is able to oxidize up to two methane molecules per cluster at ambient pressure., Competing Interests: The authors declare no competing financial interest., (© 2021 The Authors. Published by American Chemical Society.)

We have developed a new set of norm-conserving pseudopotentials and companion Gaussian basis sets for the actinide (An) series (Ac-Lr) using the Goedecker, Teter, and Hutter (GTH) formalism with the Perdew, Burke, and Ernzerhof (PBE) exchange-correlation functional of generalized gradient approximation. To test the accuracy and reliability of the newly parameterized An-GTH pseudopotentials and basis sets, a variety of benchmarks on actinide-containing molecules were carried out and compared to all-electron and available experimental results. The new pseudopotentials include both medium- ([Xe]4f 14 ) and large-core ([Xe]4f 14 5d 10 ) options that successfully reproduce the structures and energetics, particularly redox processes. The medium-core size set, in particular, reproduces all-electron calculations over multiple oxidation states from 0 to VII, whereas the large-core set is suitable only for the early series elements and low oxidation states. The underlying reason for these transferability issues is discussed in detail. This work fills a critical void in the literature for studying the chemistry of 5f-block elements in the condensed phase.

In catalysis, MgO is often used to modify the acid-base properties of support oxides and to stabilize supported metal atoms and particles on oxides. In this study, we show how the sublimation of MgO powder can be used to deposit MgO monomers, hither on anatase TiO 2 (101). A combination of x-ray electron spectroscopy, high-resolution scanning tunneling microscopy, and density functional theory is employed to gain insight into the MgO monomer binding, electronic and vibrational properties, and thermal stability. In the most stable configuration, the Mg and O of the MgO monomer bind to two surface oxygens and one undercoordinated surface titanium, respectively. The additional binding weakens the Mg-O monomer bond and makes Mg more ionic. The monomers are thermally stable up to 600 K, where the onset of diffusion into the TiO 2 bulk is observed. The monomeric MgO species on TiO 2 (101) represent an ideal atomically precise system with modified acid-base properties and will be employed in our future catalytic studies.

Chemical reactivity and sorption in zeolites are coupled to confinement and-to a lesser extent-to the acid strength of Brønsted acid sites (BAS). In presence of water the zeolite Brønsted acid sites eventually convert into hydronium ions. The gradual transition from zeolite Brønsted acid sites to hydronium ions in zeolites of varying pore size is examined by ab initio molecular dynamics combined with enhanced sampling based on Well-Tempered Metadynamics and a recently developed set of collective variables. While at low water content (1-2 water/BAS) the acidic protons prefer to be shared between zeolites and water, higher water contents (n > 2) invariably lead to solvation of the protons within a localized water cluster adjacent to the BAS. At low water loadings the standard free energy of the formed complexes is dominated by enthalpy and is associated with the acid strength of the BAS and the space around the site. Conversely, the entropy increases linearly with the concentration of waters in the pores, favors proton solvation and is independent of the pore size/shape.

Understanding how the ligand shell controls low-frequency electron paramagnetic resonance (EPR) spectroscopic properties of metal ions is essential if they are to be used in EPR-based bioimaging schemes. In this work, we probe how specific variations in the ligand structure impact L-band (ca. 1.3 GHz) EPR spectroscopic linewidths in the trichloride salts of five Cr(iii) complexes: [Cr(RR-dphen)3]3+ (RR-dphen = (1R,2R)-(+)-diphenylethylenediamine, 1), [Cr(en)3]3+ (en = ethylenediamine, 2), [Cr(me-en)3]3+ (me-en = 1,2-diaminopropane, 3), [Cr(tn)3]3+ (tn = 1,3-diaminopropane, 4) [Cr(trans-chxn)3]3+ (trans-chxn = trans-(±)-1,2-diaminocyclohexane, 5). Spectral broadening varies in a nonintuitive manner across the series, showing the sharpest peaks for 1 and broadest for 5. Molecular dynamics simulations provide evidence that the broadening is correlated to rigidity in the inner coordination sphere and reflected in ligand-dependent distribution of Cr-N bond distances that can be found in frozen solution.

Single-atom catalysts are often reported to have catalytic properties that surpass those of nanoparticles, while a direct comparison of sites common and different for both is lacking. Here we show that single atoms of Pt-group metals embedded into the surface of Fe 3 O 4 have a greatly enhanced interaction strength with CO 2 compared with the Fe 3 O 4 surface. The strong CO 2 adsorption on single Rh atoms and corresponding low activation energies lead to 2 orders of magnitude higher conversion rates of CO 2 compared to Rh nanoparticles. This high activity of single atoms stems from the partially oxidic state imposed by their coordination to the support. Fe 3 O 4 -supported Rh nanoparticles follow the behavior of single atoms for CO 2 interaction and reduction, which is attributed to the dominating role of partially oxidic sites at the Fe 3 O 4 -Rh interface. Thus, we show a likely common catalytic chemistry for two kinds of materials thought to be different, and we show that single atoms of Pt-group metals on Fe 3 O 4 are especially successful materials for catalyzed reactions that depend primarily upon sites with the metal-O-Fe environment.

A key problem associated with the design of graphene oxide (GO) materials and their tuning for nanoscale separations is how specific functional groups influence the competitive adsorption of solvated ions and water at liquid/graphene interfaces. Computation accompanied by experiment shows that OH and COOH exert an influence on water adsorption properties stronger than that of O and H functional groups. The COO - anions, following COOH deprotonation, stabilize Pb(II) through strong electrostatic interactions. This suggests that, among the functional groups under study, COOH offers the best Pb(II) adsorption capacity and the ability to regenerate the sorbent through a pH swing. In line with computation, striking experimental observations revealed that a substantial increase in Pb(II) adsorption occurs with increasing pH. Our findings provide a systematic framework for controlled design and implementation of regenerable C-based sorbents used in separations and desalination., Competing Interests: The authors declare no competing financial interest., (© 2021 The Authors. Published by American Chemical Society.)

To resolve the fleeting structures of lanthanide Ln 3+ aqua ions in solution, we (i) performed the first ab initio molecular dynamics (AIMD) simulations of the entire series of Ln 3+ aqua ions in explicit water solvent using pseudopotentials and basis sets recently optimized for lanthanides and (ii) measured the symmetry of the hydrating waters about Ln 3+ ions (Nd 3+ , Dy 3+ , Er 3+ , Lu 3+ ) for the first time with extended X-ray absorption fine structure (EXAFS). EXAFS spectra were measured experimentally and generated from AIMD trajectories to directly compare simulation, which concurrently considers the electronic structure and the atomic dynamics in solution, with experiment. We performed a comprehensive evaluation of EXAFS multiple-scattering analysis (up to 6.5 Å) to measure Ln-O distances and angular correlations (i.e., symmetry) and elucidate the molecular geometry of the first hydration shell. This evaluation, in combination with symmetry-dependent L 3 - and L 1 -edge spectral analysis, shows that the AIMD simulations remarkably reproduces the experimental EXAFS data. The error in the predicted Ln-O distances is less than 0.07 Å for the later lanthanides, while we observed excellent agreement with predicted distances within experimental uncertainty for the early lanthanides. Our analysis revealed a dynamic, symmetrically disordered first coordination shell, which does not conform to a single molecular geometry for most lanthanides. This work sheds critical light on the highly elusive coordination geometry of the Ln 3+ aqua ions.

Hierarchically ordered oxides are of critical importance in material science and catalysis. Unfortunately, the design and synthesis of such systems remains a key challenge to realizing their potential. In this study, we demonstrate how the deposition of small oligomeric (MoO 3 ) 1-6 clusters-formed by the facile sublimation of MoO 3 powders-leads to the self-assembly of locally ordered arrays of immobilized mono-oxo (MoO 3 ) 1 species on anatase TiO 2 (101). Using both high-resolution imaging and theoretical calculations, we reveal the dynamic behavior of the oligomers as they spontaneously decompose at room temperature, with the TiO 2 surface acting as a template for the growth of this hierarchically structured oxide. Transient mobility of the oligomers on both bare and (MoO 3 ) 1 -covered TiO 2 (101) areas is identified as key to the formation of a complete (MoO 3 ) 1 overlayer with a saturation coverage of one (MoO 3 ) 1 per two undercoordinated surface Ti sites. Simulations reveal a dynamic coupling of the reaction steps to the TiO 2 lattice fluctuations, the absence of which kinetically prevents decomposition. Further experimental and theoretical characterizations demonstrate that (MoO 3 ) 1 within this material are thermally stable up to 500 K and remain chemically identical with a single empty gap state produced within the TiO 2 band structure. Finally, we see that the constituent (MoO 3 ) 1 of this material show no proclivity for step and defect sites, suggesting they can reliably be grown on the (101) facet of TiO 2 nanoparticles without compromising their chemistry., Competing Interests: The authors declare no competing interest.

The hydrogenation of benzaldehyde to benzyl alcohol on carbon-supported metals in water, enabled by an external potential, is markedly promoted by polarization of the functional groups. The presence of polar co-adsorbates, such as substituted phenols, enhances the hydrogenation rate of the aldehyde by two effects, that is, polarizing the carbonyl group and increasing the probability of forming a transition state for H addition. These two effects enable a hydrogenation route, in which phenol acts as a conduit for proton addition, with a higher rate than the direct proton transfer from hydronium ions. The fast hydrogenation enabled by the presence of phenol and applied potential overcompensates for the decrease in coverage of benzaldehyde caused by competitive adsorption. A higher acid strength of the co-adsorbate increases the intensity of interactions and the rates of selective carbonyl reduction., (© 2020 The Authors. Published by Wiley-VCH GmbH.)

Sustainable energy generation calls for a shift away from centralized, high-temperature, energy-intensive processes to decentralized, low-temperature conversions that can be powered by electricity produced from renewable sources. Electrocatalytic conversion of biomass-derived feedstocks would allow carbon recycling of distributed, energy-poor resources in the absence of sinks and sources of high-grade heat. Selective, efficient electrocatalysts that operate at low temperatures are needed for electrocatalytic hydrogenation (ECH) to upgrade the feedstocks. For effective generation of energy-dense chemicals and fuels, two design criteria must be met: (i) a high H:C ratio via ECH to allow for high-quality fuels and blends and (ii) a lower O:C ratio in the target molecules via electrochemical decarboxylation/deoxygenation to improve the stability of fuels and chemicals. The goal of this review is to determine whether the following questions have been sufficiently answered in the open literature, and if not, what additional information is required:(1)What organic functionalities are accessible for electrocatalytic hydrogenation under a set of reaction conditions? How do substitutions and functionalities impact the activity and selectivity of ECH?(2)What material properties cause an electrocatalyst to be active for ECH? Can general trends in ECH be formulated based on the type of electrocatalyst?(3)What are the impacts of reaction conditions (electrolyte concentration, pH, operating potential) and reactor types?

Water-lean CO 2 capture solvents show promise for more efficient and cost-effective CO 2 capture, although their long-term behavior in operation has yet to be well studied. New observations of extended structure solvent behavior show that some solvent formulations transform into a glass-like phase upon aging at operating temperatures after contact with CO 2 . The glassification of a solvent would be detrimental to a carbon-capture process due to plugging of infrastructure, introducing a critical need to decipher the underlying principles of this phenomenon to prevent it from happening. We present the first integrated theoretical and experimental study to characterize the nano-structure of metastable and glassy states of an archetypal single-component alkanolguanidine carbon-capture solvent and assess how minute changes in atomic-level interactions convert the solvent between metastable and glass-like states. Small-angle neutron scattering and neutron diffraction coupled with small- and wide-angle X-ray scattering analysis demonstrate that minute structural changes in solution precipitae reversible aggregation of zwitterionic alkylcarbonate clusters in solution. Our findings indicate that our test system, an alkanolguanidine, exhibits a first-order phase transition, similar to a glass transition, at approximately 40 °C-close to the operating absorption temperature for post-combustion CO 2 capture processes. We anticipate that these phenomena are not specific to this system, but are present in other classes of colvents as well. We discuss how molecular-level interactions can have vast implications for solvent-based carbon-capture technologies, concluding that fortunately in this case, glassification of water-lean solvents can be avoided as long as the solvent is run above its glass transition temperature.

An accurate prediction of NMR chemical shifts at affordable computational cost is very important for different types of structural assignments in experimental studies. Density functional theory (DFT) and gauge-including atomic orbital (GIAO) are two of the most popular computational methods for NMR calculation, yet they often fail to resolve ambiguities in structural assignments. Here, we present a new method that uses machine learning (ML) techniques (DFT + ML) that significantly increases the accuracy of 13 C/ 1 H NMR chemical shift prediction for a variety of organic molecules. The input of the generalizable DFT + ML model contains two critical parts: one is a vector providing insights into chemical environments, which can be evaluated without knowing the exact geometry of the molecule; the other one is the DFT-calculated isotropic shielding constant. The DFT + ML model was trained with a data set containing 476 13 C and 270 1 H experimental chemical shifts. For the DFT methods used here, the root mean square deviations (RMSDs) for the errors between predicted and experimental 13 C/ 1 H chemical shifts can be as small as 2.10/0.18 ppm, which is much lower than those from simple DFT (5.54/0.25 ppm), or DFT + linear regression (LR) (4.77/0.23 ppm) approaches. It also has a smaller maximum absolute error than two previously proposed NMR-predicting ML models. The robustness of the DFT + ML model is tested on two classes of organic molecules (TIC10 and hyacinthacines), where the correct isomers were unambiguously assigned to the experimental ones. Overall, the DFT + ML model shows promise for structural assignments in a variety of systems, including stereoisomers, that are often challenging to determine experimentally.

Precise functionalization of electrodes with size-selected ionic liquid (IL) clusters may improve the application of ILs in electrochemical separations. Herein we report our combined experimental and theoretical investigation of the IL clusters 1-ethyl-3-methylimidazolium tetrafluoroborate [EMIM] n [BF 4 ] n +1 - ( n = 1-9) and demonstrate their selectivity and efficiency toward targeted adsorption of ions from solution. The structures and energies of the IL clusters, predicted with global optimization, agree with and help interpret the ion abundances and stabilities measured by high-mass-resolution electrospray ionization mass spectrometry and collision-induced dissociation experiments. The [EMIM][BF 4 ] 2 - cluster, which was identified as the most stable IL cluster, was selectively soft-landed onto a working electrode. Electrochemical impedance spectroscopy revealed a lower charge transfer resistance on the soft-landed electrode containing [EMIM][BF 4 ] 2 - compared with an electrode prepared by drop-casting of an IL solution containing the full range of IL clusters. Our findings indicate that specific IL clusters may be used to increase the efficiency of electrochemical separations by lowering the overpotentials involved.

Capturing carbon dioxide from post-combustion gas streams is an energy-intensive process that is required prior to either converting or sequestering CO 2 . Although a few commercial 1st and 2nd generation aqueous amine technologies have been proposed, the cost of capturing CO 2 with these technologies remains high. One approach to decrease costs of capture has been the development of water-lean solvents that aim to increase efficiency by reducing the water content in solution. Water-lean solvents, such as γ-aminopropyl aminosilicone/triethylene glycol (GAP/TEG), are promising technologies, with the potential to halve the parasitic load to a coal-fired power plant, albeit only if high solution viscosities and hydrolysis of the siloxane moieties can be mitigated. This study concerns an integrated multidisciplinary approach to overhaul the GAP/TEG solvent system at the molecular level to mitigate hydrolysis while also reducing viscosity. Cosolvents and diluents are found to have negligible effects on viscosity and are not needed. This finding allows for the design of single-component siloxane-free diamine derivatives with site-specific incorporation of selective chemical moieties for direct placement and orientation of hydrogen bonding to reduce viscosity. Ultimately, these new formulations are less susceptible to hydrolysis and exhibit up to a 98 % reduction in viscosity compared to the initial GAP/TEG formulation., (© 2020 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim.)

Global optimization constitutes an important and fundamental problem in theoretical studies in many chemical fields, such as catalysis, materials, or separations problems. In this paper, a novel algorithm has been developed for the global optimization of large systems including neat and ligated clusters in the gas phase and supported clusters in periodic boundary conditions. The method is based on an updated artificial bee colony (ABC) algorithm method, that allows for adaptive-learning during the search process. The new algorithm is tested against four classes of systems of diverse chemical nature: gas phase Au 55 , ligated Au 8 2+ , Au 8 supported on graphene oxide and defected rutile, and a large cluster assembly [Co 6 Te 8 (PEt 3 ) 6 ][C 60 ] n , with sizes ranging between 1 and 3 nm and containing up to 1300 atoms. Reliable global minima (GMs) are obtained for all cases, either confirming published data or reporting new lower energy structures. The algorithm and interface to other codes in the form of an independent program, Northwest Potential Energy Search Engine (NWPEsSe), is freely available, and it provides a powerful and efficient approach for global optimization of nanosized cluster systems.

CO 2 has attracted considerable attention in recent years due to its role in the greenhouse effect and environmental management. While its reaction with water has been studied extensively, the same cannot be said for reactivity in the supercritical CO 2 phase, where the conjugate acid/base equilibria proceed through different mechanisms and activation barriers. In spite of the apparent simplicity of the CO 2 + H 2 O reaction, the collective effect of different environments has a drastic influence on the free energy profile. Enhanced sampling techniques and well-tailored collective variables provide a detailed picture of the enthalpic and entropic drivers underscoring the differences in the formation mechanism of carbonic acid in the gas, aqueous, and supercritical CO 2 phases.

Supported single-atom catalysts (SACs) have gained increasing attention for improved catalytic activity and selectivity for industrially relevant reactions. In this study, we explore the hydrogenation of acetylene over single Pt, Ru, Rh, Pd, and Ir atoms supported on the Fe 3 O 4 (001) surface using density functional theory calculations. The thermodynamic profile of H diffusion is significantly modified by the type of single metal atoms used, suggesting that H spillover from the single atom dopant to the Fe 3 O 4 (001) surface is favored and will likely lead to high H coverages of the functioning catalyst. Correspondingly, as the surface H coverage increases, the important desorption step of ethylene becomes energetically competitive against the detrimental hydrogenation steps of ethylene to ethane. A kinetic model is employed to explore how the activity and selectivity of SACs toward ethylene production change as a function of mass of the catalyst loaded into a flow reactor. Overall, we show that the selectivity of SACs toward ethylene production can be tuned by considering the proper type of metal and controlling the redox state of the support.

Implicit solvent models are a computationally efficient method of representing solid/liquid interfaces prevalent in electrocatalysis, energy storage, and materials science. However, electronic structure changes induced at the metallic surface by the dielectric continuum are not fully understood. To address this, we perform DFT calculations for the Pt(111)/water interface, in order to compare Poisson-Boltzmann continuum solvation methods with ab initio molecular dynamics (AIMD) simulations of explicit solvent. We show that the implicit solvent cavity can be parametrized in terms of the electric dipole moment change at the equilibrated explicit Pt/water interface to obtain the potential of zero charge (PZC). We also compare the accuracy of aqueous enthalpies of adsorption of phenol on Pt(111) using geometry and charge density based dielectric cavitation methods. The ability to parametrize the cavity according to individual atoms, as afforded in the geometry based approach, is key to obtaining accurate enthalpy changes of adsorption under aqueous conditions. We also show that the electronic structure changes induced by explicit solvent and our proposed implicit solvent parametrization scheme yield comparable density difference profiles and d-band projected density of states. We therefore demonstrate the capability of implicit solvent approaches to capture both the energetics of adsorption processes and the main electronic effects of aqueous solvent on the metallic surface. This work therefore provides a scheme for computationally efficient simulations of interfacial processes for applications in areas such as heterogeneous catalysis and electrochemistry.

The effect of co-mingled dopants, Co(II) and Cr(III), on Tc(IV) incorporation and retention in magnetite under varying temperatures (75-700 °C) was explored using ab initio molecular dynamics simulations, batch experiments, and solid phase characterization. Tc(IV) stabilization was achieved with a magnetite surface oversaturated with or containing an equal number of Tc and Cr. Under oversaturation conditions, the forced formation of a Cr 2 O 3 phase on the magnetite surface may help prevent Tc release. Upon Co addition, and depending on the relative concentration of Tc, Cr, and Co at the magnetite surface, Co was found to preferentially stabilize Cr rather than Tc and suppress the formation of the protective Cr 2 O 3 surface phase. Only systems with similar Cr/Co concentrations or relatively high Cr concentrations stabilized Tc within magnetite. As such, the relative concentration of Tc, Cr, and Co was identified as a critical parameter for maximizing dopant efficacy towards Tc stabilization in magnetite., (Published by Elsevier B.V.)

Acid functionalization of a carbon support allows to enhance the electrocatalytic activity of Pd to hydrogenate benzaldehyde to benzyl alcohol proportional to the concentration of Brønsted-acid sites. In contrast, the hydrogenation rate is not affected when H 2 is used as a reduction equivalent. The different responses to the catalyst properties are shown to be caused by differences in the hydrogenation mechanism between the electrochemical and the H 2 -induced hydrogenation pathways. The enhancement of electrocatalytic reduction is realized by the participation of support-generated hydronium ions in the proximity of the metal particles., (© 2020 The Authors. Published by Wiley-VCH Verlag GmbH & Co. KGaA.)

A complete set of pseudopotentials and accompanying basis sets for all lanthanide elements are presented based on the relativistic, norm-conserving, separable, dual-space Gaussian-type pseudopotential protocol of Goedecker, Teter, and Hutter (GTH) within the generalized gradient approximation (GGA) and the exchange-correlation functional of Perdew, Burke, and Ernzerhof (PBE). The corresponding basis sets have been molecularly optimized (MOLOPT) using a contracted form with a single set of Gaussian exponents for the s, p, and d states. The f states are uncontracted explicitly with Gaussian exponents. Moreover, the Hubbard U values for each lanthanide element, for DFT+ U calculations, are also tabulated, allowing for the proper treatment of the strong on-site Coulomb interactions of localized 4f electrons. The accuracy and reliability of our GTH pseudopotentials and companion basis sets optimized for lanthanides is illustrated by a series of test calculations on lanthanide-centered molecules, and solid-state systems, with the most common oxidation states. We anticipate that these pseudopotentials and basis sets will enable larger-scale density functional theory calculations and ab initio molecular dynamics simulations of lanthanide molecules in either gas or condensed phases, as well as of solid state lanthanide-containing materials, allowing further exploration of the chemical and physical properties of lanthanide systems.

The assembly mechanism of  Metal-Organic Frameworks (MOFs) is controlled by the choice of solvent and the presence of spectator ions. In this paper, we apply  enhanced sampling molecular dynamics methods to investigate the role of solvent and ions in the early stages of the synthesis of MIL-101(Cr). Microsecond-long well-tempered metadynamics simulations uncover a  rich  structural free energy landscape, with secondary building units (SBUs) adopting distinct crystal and noncrystal like configurations. In the presence of ions (Na + , F - ), we observe a complex effect on the crystallinity of SBUs. By  modulating the interactions between terephthalate linkers and Cr atoms, ions affect the abundance of crystal-like SBUs, consequently controlling the percentage of defects. Solvent effects are assessed by comparing water with   N, N-dimethylformamide, in which SBU adducts are appreciably more stable and compact. These results shed light on how solvent and ionic strength impact the free energy of the assembly phenomena that ultimately control material synthesis.

The thermodynamic state of H 2 adsorbed on Pt in the aqueous phase was determined by kinetic analysis of H 2 reacting with D 2 O to HDO, HD, and D 2 , and by DFT-based ab initio molecular dynamics simulations of H 2 adsorption on Pt(111), Pt(110), and Pt nanoparticles. Dissociative adsorption of H 2 on Pt is significantly weakened in the aqueous phase compared to adsorption at gas-solid interfaces. Water destabilizes the adsorbed H atoms, decreasing the heat of adsorption by 19-22 kJ  m o l H 2 - 1 while inducing an additional entropy loss of 50-70 J  m o l H 2 - 1  K -1 . Upon dissociative adsorption of H 2 , the average distance of water from the Pt surface increases and the liquid adopts a structure that is more ordered than before close to the Pt surface, which limits the translation mobility of the adsorbed H atoms. The presence of hydrated hydronium ions next to the Pt surface further lowers the H-Pt bond strength., (© 2019 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim.)

Efficient regeneration of organolithium compounds is a challenging aspect in the process of novel organometathetical catalytic cycles. One of these catalytic cycles is a newly suggested method for Mg production from seawater that capitalizes on the rich chemistry of Grignard reagents. The proposed three-step catalytic cycle with Cp 2 MCl L catalyst ( M = Ti, Zr; L = select organic ligands) requires the regeneration of n BuLi from Li(s), butene, and H 2 . The potential of this approach is evaluated with density functional theory-based molecular simulations. The results reveal that the high affinity of Li toward Cl and N results in the formation of alkanes, and the strong coupling between the catalyst and BuLi leads to catalyst deactivation. To improve its catalytic performance, we proposed the use of a diamine cocatalyst and a modified catalyst with a ligand that does not contain N, which would help release BuLi from the vicinity of the catalytic center. Ab initio molecular dynamics simulations at 298 K in explicit solvent (THF) were used to estimate the Gibbs free energetics and equilibrium constants obtained from the vibrational density of states using velocity autocorrelation functions. The results show a marked improvement in the free energetics with lower barriers toward the completion of the catalytic cycle and suppression of deactivation channels.