Silicon Radical-Induced CH4 Dissociation for Uniform Graphene Coating on Silica Surface.

Due to the manufacturability of highly well-defined structures and wide-range versatility in its microstructure, SiO2 is an attractive template for synthesizing graphene frameworks with the desired pore structure. However, its intrinsic inertness constrains the graphene formation via methane chemical vapor deposition. This work overcomes this challenge by successfully achieving uniform graphene coating on a trimethylsilyl-modified SiO2 (denote TMS-MPS). Remarkably, the onset temperature for graphene growth dropped to 720 °C for the TMS-MPS, as compared to the 885 °C of the pristine SiO2. This is found to be mainly from the Si radicals formed from the decomposition of the surface TMS groups. Both experimental and computational results suggest a strong catalytic effect of the Si radicals on the CH4 dissociation. The surface engineering of SiO2 templates facilitates the synthesis of high-quality graphene sheets. As a result, the graphene-coated SiO2 composite exhibits a high electrical conductivity of 0.25 S cm-1. Moreover, the removal of the TMP-MPS template has released a graphene framework that replicates the parental TMS-MPS template on both micro- and nano- scales. This study provides tremendous insights into graphene growth chemistries as well as establishes a promising methodology for synthesizing graphene-based materials with pre-designed microstructures and porosity.

The impact of stoichiometry on the initial steps of crystal formation: Stability and lifetime of charged triple-ion complexes.

Minerals form in natural systems from solutions with varying ratios of their lattice ions, yet non-stoichiometric conditions have generally been overlooked in investigations of new formation (nucleation) of ionic crystals. Here, we investigated the influence of cation:anion ratio in the solution on the initial steps of nucleation by studying positively and negatively charged triple ion complexes and subsequent particle size evolution. Our model systems are carbonates and sulfates of calcium and barium, as it was recently shown that solution stoichiometry affects the timing and rate of their nucleation. Molecular dynamics (MD) simulations and dynamic light scattering (DLS) flow experiments show that nucleation correlates with the stability and lifetime of the initial complexes, which were significantly impacted by the cation:anion stoichiometry and ion type. Specifically, B a S O 4 2 2 - ${{\rm B}{\rm a}{\left({{\rm S}{\rm O}}_{4}\right)}_{2}^{2-}}$ was found to have higher association constants and its lifetime was twofold longer than B a 2 S O 4 2 + ${{{\rm B}{\rm a}}_{2}{{\rm S}{\rm O}}_{4}^{2+}}$ . Similar trends were observed for B a C O 3 ${{{\rm B}{\rm a}{\rm C}{\rm O}}_{3}}$ and C a S O 4 ${{{\rm C}{\rm a}{\rm S}{\rm O}}_{4}}$ . Contrastingly, for C a C O 3 ${{{\rm C}{\rm a}{\rm C}{\rm O}}_{3}}$ , C a C O 3 2 2 - ${{\rm C}{\rm a}{\left({{\rm C}{\rm O}}_{3}\right)}_{2}^{2-}}$ was found to have lower association constants and its lifetime was shorter than C a 2 C O 3 2 + ${{{\rm C}{\rm a}}_{2}{{\rm C}{\rm O}}_{3}^{2+}}$ . These trends in stability and lifetime follow the same asymmetrical behaviour as observed experimentally for particle formation using techniques like DLS. This suggests a causal relationship between the stability and lifetime of the initial charged complexes and the nucleation under non-stoichiometric conditions.

Single-Atom Iridium on Hematite Photoanodes for Solar Water Splitting: Catalyst or Spectator?

Single-atom catalysts (SACs) on hematite photoanodes are efficient cocatalysts to boost photoelectrochemical performance. They feature high atom utilization, remarkable activity, and distinct active sites. However, the specific role of SACs on hematite photoanodes is not fully understood yet: Do SACs behave as a catalytic site or as a spectator? By combining spectroscopic experiments and computer simulations, we demonstrate that single-atom iridium (sIr) catalysts on hematite (α-Fe2O3/sIr) photoanodes act as a true catalyst by trapping holes from hematite and providing active sites for the water oxidation reaction. In situ transient absorption spectroscopy showed a reduced number of holes and shortened hole lifetime in the presence of sIr. This was particularly evident on the second timescale, indicative of fast hole transfer and depletion toward water oxidation. Intensity-modulated photocurrent spectroscopy evidenced a faster hole transfer at the α-Fe2O3/sIr/electrolyte interface compared to that at bare α-Fe2O3. Density functional theory calculations revealed the mechanism for water oxidation using sIr as a catalytic center to be the preferred pathway as it displayed a lower onset potential than the Fe sites. X-ray photoelectron spectroscopy demonstrated that sIr introduced a mid-gap of 4d state, key to the fast hole transfer and hole depletion. These combined results provide new insights into the processes controlling solar water oxidation and the role of SACs in enhancing the catalytic performance of semiconductors in photo-assisted reactions.

Molecular Insights into the Stability of Titanium in Electrolytes Containing Chlorine and Fluorine Ions.

Titanium and its alloys are protected by a compact and stable passive film, which confers resistance to corrosion by the primary halogen chloride (Cl-) while being less effective against fluoride (F-). Although researchers have recognized different macroscopic corrosion effects of these halide ions on titanium, the underlying mechanisms remain largely unexplored. In this work, the bonding of Cl-/F- with stable passive films was studied in neutral and acidic (pH = 2.3) conditions. The synergistic effect between the interfacial hydrogen bond (HB) structure and halogens on titanium corrosion was first revealed using first-principles calculation and Raman spectroscopy. F- forms more stable halogen-Ti bonds than Cl-, resulting in titanium degradation. The proton combined with F- exhibits a specific synergistic effect, causing corrosion of the passive film. The water hydrogen bond transformation index (HBTI) at the titanium/aqueous interface was 1.88 in an acidic solution containing F-, significantly higher than that in neutral/acid solutions containing Cl- (1.80/1.81) and a neutral solution containing F- (1.81). This work clarifies the structure-activity relationship between HBTI and the destruction of titanium passive films. We propose that the microstructure of the interfacial HB is an undeniable factor in the corrosion of titanium.

Adsorption, activation, and conversion of carbon dioxide on small copper-tin nanoclusters.

Carbon dioxide (CO2) conversion to value-added chemicals is an attractive solution to reduce globally accelerating CO2 emissions. Among the non-precious and abundant metals tested so far, copper (Cu) is one of the best electrocatalysts to convert CO2 into more than thirty different hydrocarbons and alcohols. However, the selectivity for desired products is often too low. We present a computational investigation of the effects of nanostructuring, doping, and support on the activity and selectivity of Cu-Sn catalysts. Density functional theory calculations were conducted to explore the possibility of using small Cu-Sn clusters, Cu4-nSnn (n = 0-4), isolated or supported on graphene and γ-Al2O3, to activate CO2 and convert it to carbon monoxide (CO) and formic acid (HCOOH). First, a detailed analysis of the structure, stability, and electronic properties of Cu4-nSnn clusters and their ability to absorb and activate CO2 was considered. Then, the kinetics of the gas phase CO2 direct dissociation on Cu4-nSnn to generate CO was determined. Finally, the mechanism of electrocatalytic CO2 reduction to CO and HCOOH on Cu4-nSnn, Cu4-nSnn/graphene and Cu4-nSnn/γ-Al2O3 was computed. The selectivity towards the competitive electrochemical hydrogen evolution reaction on these catalysts was also considered. The Cu2Sn2 cluster suppresses the hydrogen evolution reaction and is highly selective towards CO, if unsupported, or HCOOH if supported on graphene. This study demonstrates that the Cu2Sn2 cluster is a potential candidate for the electrocatalytic conversion of the CO2 molecule. Moreover, it identifies insightful structure-property relationships in Cu-based nanocatalysts, highlighting the influence of composition and catalyst support on CO2 activation.

The carbon chain growth during the onset of CVD graphene formation on γ-Al2O3 is promoted by unsaturated CH2 ends.

Chemical vapor deposition of methane onto a template of alumina (Al2O3) nanoparticles is a prominent synthetic strategy of graphene meso-sponge, a new class of nano porous carbon materials consisting of single-layer graphene walls. However, the elementary steps controlling the early stages of graphene growth on Al2O3 surfaces are still not well understood. In this study, density functional theory calculations provide insights into the initial stages of graphene growth. We have modelled the mechanism of CH4 dissociation on the (111), (110), (100), and (001) γ-Al2O3 surfaces. Subsequently, we have considered the reaction pathway leading to the formation of a C6 ring. The γ-Al2O3(110) and γ-Al2O3(100) are both active for CH4 dissociation, but the (100) surface has higher catalytic activity towards the carbon growth reaction. The overall mechanism involves the formation of the reactive intermediate CH2* that then can couple to form CnH2n* (n = 2-6) intermediates with unsaturated CH2 ends. The formation of these species, which are not bound to the surface-active sites, promotes the sustained carbon growth in a nearly barrierless process. Also, the short distance between terminal carbon atoms leads to strong interactions, which might lead to the high activity between unsaturated CH2* of the hydrocarbon chain. Analysis of the electron localization and geometries of the carbon chains reveals the formation of C-Al-σ bonds with the chain growing towards the vacuum rather than C-Al-π bonds covering the γ-Al2O3(100) surface. This growth behaviour prevents catalyst poisoning during the initial stage of graphene nucleation.

Bridging atomistic simulations and thermodynamic hydration models of aqueous electrolyte solutions.

Chemical thermodynamic models of solvent and solute activities predict the equilibrium behavior of aqueous solutions. However, these models are semi-empirical. They represent micro-scale ion and solvent behaviors controlling the macroscopic properties using small numbers of parameters whose values are obtained by fitting to activities and other partial derivatives of the Gibbs energy measured for the bulk solutions. We have conducted atomistic simulations of aqueous electrolyte solutions (MgCl2 and CaCl2) to determine the parameters of thermodynamic hydration models. We have implemented a cooperative hydration model to categorize the water molecules in electrolyte solutions into different subpopulations. The value of the electrolyte-specific parameter, k, was determined from the ion-affected subpopulation with the lowest absolute value of the free energy of removing the water molecule. The other equilibrium constant parameter, K1, associated with the first degree of hydration, was computed from the free energy of hydration of hydrated clusters. The hydration number, h, was determined from a reorientation dynamic analysis of the water subpopulations compared to bulk-like behavior. The reparameterized models [R. H. Stokes and R. H. Robinson, J. Solution Chem. 2, 173 (1973) and Balomenos et al., Fluid Phase Equilib. 243, 29 (2006)] using the computed values of the parameters lead to the osmotic coefficients of MgCl2 solutions that are consistent with measurements. Such an approach removes the dependence on the availability of experimental data and could lead to aqueous thermodynamic models capable of estimating the values of solute and solvent activities as well as thermal and volumetric properties for a wide range of compositions and concentrations.

Hydrogen-Bond Structure and Low-Frequency Dynamics of Electrolyte Solutions: Hydration Numbers from ab Initio Water Reorientation Dynamics and Dielectric Relaxation Spectroscopy.

We present an atomistic simulation scheme for the determination of the hydration number (h) of aqueous electrolyte solutions based on the calculation of the water dipole reorientation dynamics. In this methodology, the time evolution of an aqueous electrolyte solution generated from ab initio molecular dynamics simulations is used to compute the reorientation time of different water subpopulations. The value of h is determined by considering whether the reorientation time of the water subpopulations is retarded with respect to bulk-like behavior. The application of this computational protocol to magnesium chloride (MgCl2 ) solutions at different concentrations (0.6-2.8 mol kg-1 ) gives h values in excellent agreement with experimental hydration numbers obtained using GHz-to-THz dielectric relaxation spectroscopy. This methodology is attractive because it is based on a well-defined criterion for the definition of hydration number and provides a link with the molecular-level processes responsible for affecting bulk solution behavior. Analysis of the ab initio molecular dynamics trajectories using radial distribution functions, hydrogen bonding statistics, vibrational density of states, water-water hydrogen bonding lifetimes, and water dipole reorientation reveals that MgCl2 has a considerable influence on the hydrogen bond network compared with bulk water. These effects have been assigned to the specific strong Mg-water interaction rather than the Cl-water interaction.

Density functional theory based molecular dynamics study of solution composition effects on the solvation shell of metal ions.

We present an ab initio molecular dynamics study of the alkali metal ions Li+, Na+, K+ and Cs+, and of the alkaline earth metal ions Mg2+ and Ca2+ in both pure water and electrolyte solutions containing the counterions Cl- and SO42-. Simulations were conducted using different density functional theory methods (PBE, BLYP and revPBE), with and without the inclusion of dispersion interactions (-D3). Analysis of the ion-water structure and interaction strength, water exchange between the first and second hydration shell, and hydrogen bond network and low-frequency reorientation dynamics around the metal ions have been used to characterise the influence of solution composition on the ionic solvation shell. Counterions affect the properties of the hydration shell not only when they are directly coordinated to the metal ion, but also when they are at the second coordination shell. Chloride ions reduce the sodium hydration shell and expand the calcium hydration shell by stabilizing under-coordinated hydrated Na(H2O)5+ complexes and over-coordinated Ca(H2O)72+. The same behaviour is observed in CaSO4(aq), where Ca2+ and SO42- form almost exclusively solvent-shared ion pairs. Water exchange between the first and second hydration shell around Ca2+ in CaSO4(aq) is drastically decelerated compared with the simulations of the hydrated metal ion (single Ca2+, no counterions). Velocity autocorrelation function analysis, used to probe the strength of the local ion-water interaction, shows a smoother decay of Mg2+ in MgCl2(aq), which is a clear indication of a looser inter-hexahedral vibration in the presence of chloride ions located in the second coordination shell of Mg2+. The hydrogen bond statistics and orientational dynamics in the ionic solvation shell show that the influence on the water-water network cannot only be ascribed to the specific cation-water interaction, but also to the subtle interplay between the level of hydration of the ions, and the interactions between ions, especially those of opposite charge. As many reactive processes involving solvated metal ions occur in environments that are far from pure water but rich in ions, this computational study shows how the solution composition can result in significant differences in behaviour and function of the ionic solvation shell.

Solid and Aqueous Speciation of Yttrium in Passive Remediation Systems of Acid Mine Drainage.

Yttrium belongs to the rare earth elements (REEs) together with lanthanides and scandium. REEs are commonly used in modern technologies, and their limited supply has made it necessary to look for new alternative resources. Acid mine drainage (AMD) is a potential resource since it is moderately enriched in REEs. In fact, in passive remediation systems, which are implemented to minimize the environmental impacts of AMD, REEs are mainly retained in basaluminite, an aluminum hydroxysulfate precipitate. In this study, the solid and liquid speciation and the local structure of yttrium are studied in high-sulfate aqueous solutions, basaluminite standards, and samples from remediation columns using synchrotron-based techniques and molecular modeling. Pair distribution function (PDF) analyses and ab initio molecular dynamics density functional theory models of the yttrium sulfate solution show that the YSO4+ ion pair forms a monodentate inner-sphere complex. Extended X-ray absorption fine structure (EXAFS) and PDF analyses show that Y is retained by basaluminite, forming a monodentate inner-sphere surface complex on the aluminum hydroxide surface. EXAFS of the column samples shows that more than 72% of their signal is represented by the signal of basaluminite with which YSO4+ forms an inner-sphere complex. The atomic view of the REE configuration in AMD environments could facilitate a deeper research of REE recovery from waste generated in AMD remediation systems.

Prediction of self-assembly of adenosine analogues in solution: a computational approach validated by isothermal titration calorimetry.

The recent discovery of the role of adenosine-analogues as neuroprotectants and cognitive enhancers has sparked interest in these molecules as new therapeutic drugs. Understanding the behavior of these molecules in solution and predicting their ability to self-assemble will accelerate new discoveries. We propose a computational approach based on density functional theory, a polarizable continuum solvation description of the aqueous environment, and an efficient search procedure to probe the potential energy surface, to determine the structure and thermodynamic stability of molecular clusters of adenosine analogues in solution, using caffeine as a model. The method was validated as a tool for the prediction of the impact of small structural variations on self-assembly using paraxanthine. The computational results were supported by isothermal titration calorimetry experiments. The thermodynamic parameters enabled the quantification of the actual percentage of dimer present in solution as a function of concentration. The data suggest that both caffeine and paraxanthine are present at concentrations comparable to the ones found in biological samples.

Simulations reveal the role of composition into the atomic-level flexibility of bioactive glass cements.

Bioactive glass ionomer cements (GICs), the reaction product of a fluoro-alumino-silicate glass and polyacrylic acid, have been in effective use in dentistry for over 40 years and more recently in orthopaedics and medical implantation. Their desirable properties have affirmed GIC's place in the medical materials community, yet are limited to non-load bearing applications due to the brittle nature of the hardened composite cement, thought to arise from the glass component and the interfaces it forms. Towards helping resolve the fundamental bases of the mechanical shortcomings of GICs, we report the 1st ever computational models of a GIC-relevant component. Ab initio molecular dynamics simulations were employed to generate and characterise three fluoro-alumino-silicate glasses of differing compositions with focus on resolving the atomic scale structural and dynamic contributions of aluminium, phosphorous and fluorine. Analyses of the glasses revealed rising F-content leading to the expansion of the glass network, compression of Al-F bonding, angular constraint at Al-pivots, localisation of alumino-phosphates and increased fluorine diffusion. Together, these changes to the structure, speciation and dynamics with raised fluorine content impart an overall rigidifying effect on the glass network, and suggest a predisposition to atomic-level inflexibility, which could manifest in the ionomer cements they form.

Modelling the effects of salt solutions on the hydration of calcium ions.

Classical molecular dynamics simulations of several aqueous alkali halide salt solutions have been used to determine the effect of electrolytes on the structure of water and the hydration properties of calcium ions. Compared with the simulations of Ca(2+) ions in pure liquid water, the frequency of water exchange in the first hydration shell of calcium, which is a fundamental process in controlling the reactivity of calcium(ii) aqua-ions, is drastically reduced in the presence of other electrolytes in solution. The strong stabilization of the hydration shell of Ca(2+) occurs not only when the halide anions are directly coordinated to calcium, but also when the alkali and halide ions are placed at or outside the second coordination shell of Ca(2+), suggesting that the reactivity of the first solvation shell of the calcium ion can be influenced by the specific affinity of other ions in solution for the water molecules coordinated to Ca(2+). Analysis of the hydrogen-bonded structure of water in the vicinity of the calcium ion shows that the average number of hydrogen bonds per water molecules, which is 1.8 in pure liquid water, decreases as the concentration of alkali-halide salts in solution increases, and that the temporal fluctuations of hydrogen bonds are significantly larger than those obtained for Ca(2+) in pure liquid water. This effect has been explained in terms of the dynamics of reorganization of the O-H···X(-) (X = F, Cl and Br) hydrogen bond. This work shows the importance of solution composition in determining the hydrogen-bonding network and ligand-exchange dynamics around metal ions, both in solution and at the mineral-water interfaces, which in turn has implications for interactions occurring at the mineral-water interface, ultimately controlling the mobilization of ions in the environment as well as in industrial processes.

Density functional theory study of the oligomerization of carboxylic acids.

We present a density functional theory [M06-2X/6-31+G(d,p)] study of the structures and free energies of formation of oligomers of four carboxylic acids (formic acid, acetic acid, tetrolic acid, and benzoic acid) in water, chloroform, and carbon tetrachloride. Solvation effects were treated using the SMD continuum solvation model. The low-lying energy structures of molecular complexes were located by adopting an efficient search procedure to probe the potential energy surfaces of the oligomers of carboxylic acids (CA)n (n = 2-6). The free energies of the isomers of (CA)n in solution were determined as the sum of the electronic energy, vibrational-rotational-translational gas-phase contribution, and solvation free energy. The assessment of the computational protocol adopted in this study with respect to the dimerization of acetic acid, (AA)2, and formic acid, (FA)2, located new isomers of (AA)2 and (FA)2 and gave dimerization constants in good agreement with the experimental values. The calculation of the self-association of acetic acid, tetrolic acid, and benzoic acid shows the following: (i) Classic carboxylic dimers are the most stable isomer of (CA)2 in both the gas phase and solution. (ii) Trimers of carboxylic acid are stable in apolar aprotic solvents. (iii) Molecular clusters consisting of two interacting classic carboxylic dimers (CA)4,(D+D) are the most stable type of tetramers, but their formation from the self-association of classic carboxylic dimers is highly unfavorable. (iv) For acetic acid and tetrolic acid the reactions (CA)2 + 2CA → (CA)4,(D+D) and (CA)3 + CA → (CA)4,(D+D) are exoergonic, but these aggregation pathways go through unstable clusters that could hinder the formation of tetrameric species. (v) For tetrolic acid the prenucleation species that are more likely to form in solution are dimeric and trimeric structures that have encoded structural motifs resembling the α and ß solid forms of tetrolic acid. (vi) Stable tetramers of benzoic acid could form in carbon tetrachloride from the aggregation of trimers and monomers. (vii) Higher order clusters such as acetic acid pentamers and tetrolic acid hexamers are highly unstable in all solvation environments.

Comparative Investigation of the Microstructure of MgCl2 Aqueous Solutions Using Different X-ray Scattering Sources, Raman Spectroscopy, and Atomistic Simulations.

Aqueous solutions of magnesium chloride (MgCl2(aq)) are often used to test advances in the theory of electrolyte solutions because they are considered an ideal strong 2:1 electrolyte. However, there is evidence that some ion association occurs in these solutions, even at low concentrations. Even a small ion-pairing constant can have a significant impact on the chemical speciation of ions, so it is important to determine whether ion pairing actually occurs. In this study, MgCl2(aq) with concentrations ranging from 1 to 35% was studied using three methods: X-ray scattering (XRS) with the Shanghai Synchrotron Radiation Facility (SSRF) and silver-anode laboratory sources, Raman spectroscopy, and molecular dynamics (MD) simulations with the COMPASS-II and Madrid force fields. XRS results were analyzed in the framework of PDF theory to obtain the reduced structure function F(Q) and the reduced pair distribution function G(r). The F(Q) values from synchrotron radiation and laboratory sources both showed that the tetrahedral hydrogen bonds in bulk water were destroyed with the increased MgCl2 concentration. The results of G(r) indicated that the main peaks centered at 2.05 and 2.80 Å can be ascribed to the interactions of Mg-O and O-O, respectively. The peak at 3.10 Å is attributed to the combined effect of O-O and Cl-O. By comparing the structural information on MgCl2 solution obtained from the two light sources, it was found that both SSRF and silver-anode laboratory sources can reflect the above-mentioned structural information on MgCl2 solution. The radial distribution function (RDF) obtained from MD simulations of MgCl2 solutions assigned the peaks at 2.0, 2.8, and 3.2 Å to the Mg-O, O-O, and Cl-O interatomic pairs, respectively. The decrease in the O-O coordination number confirms that the hydrogen-bonding network of water is disrupted by increasing MgCl2 observed by X-ray scattering. The proportion of Mg-Cl contact ion pairs gradually increases with MgCl2 concentration as does the coordination number. Raman spectroscopy results show that the bond type changes from double donor double acceptor (DDAA) to single donor-single acceptor (DA) with increasing concentration, providing explicit details of the hydrogen-bond evolution in the aqueous solution.

Density functional theory simulations of the structure, stability and dynamics of iron sulphide clusters in water.

Density Functional Theory-based calculations have been employed to investigate the structure, stability and dynamics of iron sulphide clusters, FexSy (x, y ≤ 4), in water. Car-Parrinello molecular dynamics simulations of the building unit FeS in explicit water show that the iron is only four-coordinated, which indicates that the effect of sulphur is to significantly reduce the coordination shell of iron compared with the typical octahedral arrangement of hexa-aqua iron complexes in water. The molecular dynamics simulations of FexSy particles (x, y ≥ 2) in explicit water reveal that these clusters are highly unstable as they dissociate after a few picoseconds. The Gibbs free energies to form the FeS and Fe2S2 species have been evaluated in a simulated aqueous environment, using the mPW1B95 density functional theory level for the gas-phase component and the UAHF-CPCM solvation model for the hydration contribution, and the results indicate that while FeS is thermodynamically stable in aqueous solution, the formation of a Fe2S2 cluster is endergonic, and dissociation is preferred under natural water conditions.

Advanced Catalyst Design and Reactor Configuration Upgrade in Electrochemical Carbon Dioxide Conversion.

Electrochemical carbon dioxide reduction reaction (CO2 RR) driven by renewable energy shows great promise in mitigating and potentially reversing the devastating effects of anthropogenic climate change and environmental degradation. The simultaneous synthesis of energy-dense chemicals can meet global energy demand while decoupling emissions from economic growth. However, the development of CO2 RR technology faces challenges in catalyst discovery and device optimization that hinder their industrial implementation. In this contribution, a comprehensive overview of the current state of CO2 RR research is provided, starting with the background and motivation for this technology, followed by the fundamentals and evaluated metrics. Then the underlying design principles of electrocatalysts are discussed, emphasizing their structure-performance correlations and advanced electrochemical assembly cells that can increase CO2 RR selectivity and throughput. Finally, the review looks to the future and identifies opportunities for innovation in mechanism discovery, material screening strategies, and device assemblies to move toward a carbon-neutral society.

CO2-mineralization and carbonation reactor rig: Design and validation for in situ neutron scattering experiments-Engineering and lessons learned.

CO2 mineralization via aqueous Mg/Ca/Na-carbonate (MgCO3/CaCO3/Na2CO3) formation represents a huge opportunity for the utilization of captured CO2. However, large-scale mineralization is hindered by slow kinetics due to the highly hydrated character of the cations in aqueous solutions (Mg2+ in particular). Reaction conditions can be optimized to accelerate carbonation kinetics, for example, by the inclusion of additives that promote competitive dehydration of Mg2+ and subsequent agglomeration, nucleation, and crystallization. For tracking mineralization and these reaction steps, neutron scattering presents unprecedented advantages over traditional techniques for time-resolved in situ measurements. However, a setup providing continuous solution circulation to ensure reactant system homogeneity for industrially relevant CO2-mineralization is currently not available for use on neutron beamlines. We, therefore, undertook the design, construction, testing and implementation of such a self-contained reactor rig for use on selected neutron beamlines at the ISIS Neutron and Muon Source (Harwell, UK). The design ensured robust attachment via suspension from the covering Tomkinson flange to stabilize the reactor assembly and all fittings (~25 kg), as well as facilitating precise alignment of the entire reactor and sample (test) cell with respect to beam dimension and direction. The assembly successfully accomplished the principal tasks of providing a continuous flow of the reaction mixture (~500 mL) for homogeneity, quantitative control of CO2 flux into the mixture, and temperature and pressure regulation throughout the reaction and measurements. The design is discussed, with emphasis placed on the reactor, including its geometry, components, and all technical specifications. Descriptions of the off-beamline bench tests, safety, and functionality, as well as the installation on beamlines and trial experimental procedure, are provided, together with representative raw neutron scattering results.

Polarizable force field development and molecular dynamics study of phosphate-based glasses.

Molecular dynamics simulations of phosphate-based glasses P(2)O(5)-CaO-Na(2)O have been carried out, using an interatomic force field that has been parameterized to reproduce the structural and mechanical properties of crystalline phosphorus pentoxide, o(')(P(2)O(5))(∞) orthorhombic phase. Polarization effects have been included through the shell-model potential and formal charges have been used to aid transferability. A modification to the DL_POLY code (version 2.20) was used to model the high temperature shell dynamics. Structural characterizations of three biomedically applicative molar compositions, (P(2)O(5))(0.45)(CaO)(x)(Na(2)O)(0.55-x) (x = 0.30, 0.35, and 0.40), have been undertaken. Good agreement with available experimental and ab initio data is obtained. The simulations show that, dependent on composition, the phosphorus atoms are primarily bonded to two or three oxygens that in turn bridge to neighbouring phosphorus atoms. Na(+) and Ca(2+) modifiers are found to occupy a pseudo-octahedral bonding environment with mean oxygen coordination numbers of 6.55 and 6.85, respectively, across all compositions studied.

Optimal Icosahedral Copper-Based Bimetallic Clusters for the Selective Electrocatalytic CO2 Conversion to One Carbon Products.

Electrochemical CO2 reduction reactions can lead to high value-added chemical and materials production while helping decrease anthropogenic CO2 emissions. Copper metal clusters can reduce CO2 to more than thirty different hydrocarbons and oxygenates yet they lack the required selectivity. We present a computational characterization of the role of nano-structuring and alloying in Cu-based catalysts on the activity and selectivity of CO2 reduction to generate the following one-carbon products: carbon monoxide (CO), formic acid (HCOOH), formaldehyde (H2C=O), methanol (CH3OH) and methane (CH4). The structures and energetics were determined for the adsorption, activation, and conversion of CO2 on monometallic and bimetallic (decorated and core@shell) 55-atom Cu-based clusters. The dopant metals considered were Ag, Cd, Pd, Pt, and Zn, located at different coordination sites. The relative binding strength of the intermediates were used to identify the optimal catalyst for the selective CO2 conversion to one-carbon products. It was discovered that single atom Cd or Zn doping is optimal for the conversion of CO2 to CO. The core@shell models with Ag, Pd and Pt provided higher selectivity for formic acid and formaldehyde. The Cu-Pt and Cu-Pd showed lowest overpotential for methane formation.