Optimizing Active Sites for High CO Selectivity during CO2 Hydrogenation over Supported Nickel Catalysts.

Controlling the selectivity of CO2 hydrogenation catalysts is a fundamental challenge. In this study, the selectivity of supported Ni catalysts prepared by the traditional impregnation method was found to change after a first CO2 hydrogenation reaction cycle from 100 to 800 °C. The usually high CH4 formation was suppressed leading to full selectivity toward CO. This behavior was also observed after the catalyst was treated under methane or propane atmospheres at elevated temperatures. In situ spectroscopic studies revealed that the accumulation of carbon species on the catalyst surface at high temperatures leads to a nickel carbide-like phase. The catalyst regains its high selectivity to CH4 production after carbon depletion from the surface of the Ni particles by oxidation. However, the selectivity readily shifts back toward CO formation after exposing the catalysts to a new temperature-programmed CO2 hydrogenation cycle. The fraction of weakly adsorbed CO species increases on the carbide-like surface when compared to a clean nickel surface, explaining the higher selectivity to CO. This easy protocol of changing the surface of a common Ni catalyst to gain selectivity represents an important step for the commercial use of CO2 hydrogenation to CO processes toward high-added-value products.

Multiple pathways in pressure-induced phase transition of coesite.

High-pressure single-crystal X-ray diffraction method with precise control of hydrostatic conditions, typically with helium or neon as the pressure-transmitting medium, has significantly changed our view on what happens with low-density silica phases under pressure. Coesite is a prototype material for pressure-induced amorphization. However, it was found to transform into a high-pressure octahedral (HPO) phase, or coesite-II and coesite-III. Given that the pressure is believed to be hydrostatic in two recent experiments, the different transformation pathways are striking. Based on molecular dynamic simulations with an ab initio parameterized potential, we reproduced all of the above experiments in three transformation pathways, including the one leading to an HPO phase. This octahedral phase has an oxygen hcp sublattice featuring 2 × 2 zigzag octahedral edge-sharing chains, however with some broken points (i.e., point defects). It transforms into α-PbO2 phase when it is relaxed under further compression. We show that the HPO phase forms through a continuous rearrangement of the oxygen sublattice toward hcp arrangement. The high-pressure amorphous phases can be described by an fcc and hcp sublattice mixture.

Retention of contaminants Cd and Hg adsorbed and intercalated in aluminosilicate clays: A first principles study.

Layered clay materials have been used to incorporate transition metal (TM) contaminants. Based on first-principles calculations, we have examined the energetic stability and the electronic properties due to the incorporation of Cd and Hg in layered clay materials, kaolinite (KAO) and pyrophyllite (PYR). The TM can be (i) adsorbed on the clay surface as well as (ii) intercalated between the clay layers. For the intercalated case, the contaminant incorporation rate can be optimized by controlling the interlayer spacing of the clay, namely, pillared clays. Our total energy results reveal that the incorporation of the TMs can be maximized through a suitable tuning of vertical distance between the clay layers. Based on the calculated TM/clay binding energies and the Langmuir absorption model, we estimate the concentrations of the TMs. Further kinetic properties have been examined by calculating the activation energies, where we found energy barriers of ∼20 and ∼130 meV for adsorbed and intercalated cases, respectively. The adsorption and intercalation of ionized TM adatoms were also considered within the deprotonated KAO surface. This also leads to an optimal interlayer distance which maximizes the TM incorporation rate. By mapping the total charge transfers at the TM/clay interface, we identify a net electronic charge transfer from the TM adatoms to the topmost clay surface layer. The effect of such a charge transfer on the electronic structure of the clay (host) has been examined through a set of X-ray absorption near edge structure (XANES) simulations, characterizing the changes of the XANES spectra upon the presence of the contaminants. Finally, for the pillared clays, we quantify the Cd and Hg K-edge energy shifts of the TMs as a function of the interlayer distance between the clay layers and the Al K-edge spectra for the pristine and pillared clays.

First principles characterization of silicate sites in clay surfaces.

Aluminosilicate clays like Montmorillonite (MMT) and Muscovite Mica (MT) have siloxane cavities on the basal plane. The hydroxyl groups localized in these cavities and van der Waals (vdW) forces contribute significantly to adsorption processes. However, the basal sites are found to be difficult to characterize experimentally. Here, (001) surfaces of MMT and MT clays were investigated using first-principles calculations to understand how these silicate surface sites are influenced by hydroxyl groups and the effective role of inner layer vdW interactions. Based on density-functional theory (DFT) within the generalized gradient approximation (GGA), different types of exchange-correlation functionals were tested to check the effect of vdW dispersion correction. Noncontact atomic force microscopy (nc-AFM), X-ray absorption spectroscopy (XAS) in the near-edge region and solid-state nuclear magnetic resonance (SS-NMR) spectroscopy were simulated. In both clays, the oxygen surface sites are directly affected by the intralayer interaction through hydroxyl groups. Our results indicated that the chemical environment of the hydroxyl groups is distinct in the MMT and MT structures. The vdW correction was essential for a better description of the surface oxygen sites and correctly describes the similarity between both clays. Particularly, the bulk apical oxygen sites in the MT structure are less influenced by vdW interaction. Compared to MMT, the silicon surface sites of MT are more sensitive to the intralayer changes in Si-Oapical-Al and with less effect of the hydroxyl groups. These results provide a clear understanding of influence of the siloxane cavity on the oxygen and silicon surface sites in aluminosilicates.

NMR characterization of hydrocarbon adsorption on calcite surfaces: a first principles study.

The electronic and coordination environment of minerals surfaces, as calcite, are very difficult to characterize experimentally. This is mainly due to the fact that there are relatively few spectroscopic techniques able to detect Ca(2+). Since calcite is a major constituent of sedimentary rocks in oil reservoir, a more detailed characterization of the interaction between hydrocarbon molecules and mineral surfaces is highly desirable. Here we perform a first principles study on the adsorption of hydrocarbon molecules on calcite surface (CaCO3 (101¯4)). The simulations were based on Density Functional Theory with Solid State Nuclear Magnetic Resonance (SS-NMR) calculations. The Gauge-Including Projector Augmented Wave method was used to compute mainly SS-NMR parameters for (43)Ca, (13)C, and (17)O in calcite surface. It was possible to assign the peaks in the theoretical NMR spectra for all structures studied. Besides showing different chemical shifts for atoms located on different environments (bulk and surface) for calcite, the results also display changes on the chemical shift, mainly for Ca sites, when the hydrocarbon molecules are present. Even though the interaction of the benzene molecule with the calcite surface is weak, there is a clearly distinguishable displacement of the signal of the Ca sites over which the hydrocarbon molecule is located. A similar effect is also observed for hexane adsorption. Through NMR spectroscopy, we show that aromatic and alkane hydrocarbon molecules adsorbed on carbonate surfaces can be differentiated.

First-Principles Prediction of Amorphous Silica Nanoparticle Surface Charge: Effect of Size, pH, and Ionic Strength.

The quantification of surface charge properties of silica nanoparticles is essential for several applications. To determine these properties, many experimental and theoretical methods have been introduced, which are time-consuming and/or challenging to use. In this study, a first-principles approach is developed to determine the surface charge properties of amorphous silica nanoparticles against the nanoparticle size, pH, and ionic strength without relying on experimental data. An amorphous silica nanoparticle of 1.34 nm diameter is simulated by using integrated molecular dynamics and Monte Carlo methods. A detailed analysis of the nanoparticle structure is provided by analyzing the types of silanol groups on the surface. Moreover, a model is developed to estimate the probability distribution of the surface silanol groups based on the nearest neighbor distances and the diameter of the nanoparticle to determine the number of surface silanols on larger nanoparticles. Thereafter, a computational chemistry approach is used to calculate the acid dissociation constants of the corresponding deprotonation reactions. The calculated constants and the point of zero charge value are in excellent agreement with experiments. The surface charge properties of the nanoparticle with various diameters are then estimated by using a mean-field model at different pH and ionic strength values. The results of the developed model are compared to the Poisson-Boltzmann equation as a reference model. The developed model predictions agree well with the reference model for low and mid-electrolyte concentrations (1 and 10 mM) and small nanoparticles (smaller than 100 nm). However, the developed model seems to qualitatively predict the surface charge properties more accurately than the Poisson-Boltzmann model for high electrolyte concentrations.

Comparison of thermodynamic stabilities and mechanical properties of CO2, SiO2, and GeO2 polymorphs by first-principles calculations.

The recent discovery that molecular CO(2) transforms under compression into carbon four-coordinated, 3-dimensional network solid phases has generated considerable interests on possible new phases in the fourth-main-group elemental oxides. Based on density-functional theory calculations, we have investigated the thermodynamic stability, mechanical properties and electronic structure of proposed guest-free clathrates, quartz and cristobalite phases for CO(2), SiO(2), and GeO(2), and the dry ice phase for CO(2). It was predicted that a GeO(2) clathrate, likely a semiconductor, could be synthesized presumably with some suitable guest molecules. The hypothetical CO(2) guest-free clathrate phase was found hardly to be formed due to the large energy difference with respect to the other polymorphs. This phase is unstable at all pressures, which is also implied by its different electronic structure in comparison with SiO(2) and GeO(2). Finally, the SiO(2) clathrate presents a uniquely high bulk modulus, which is higher than that of quartz and three times of the experimental data, might not be a weak point of ab-initio calculations such as pseudopotentials, correlation functional etc., instead it can be readily understood by the constraint as imposed by the high symmetry. Either temperature or an "exhausted" relaxation (without any symmetry constraint) can remedy this problem.

Interface tension of silica hydroxylated nanoparticle with brine: a combined experimental and molecular dynamics study.

We have used molecular dynamics simulations to calculate the interfacial tension of hydroxylated SiO(2) nanoparticles under different temperatures and solutions (helium and brine with monovalent and divalent salts). In order to benchmark the atomistic model, quartz SiO(2) interfacial tension was measured based on inverse gas chromatography under He atmosphere. The experimental interfacial tension values for quartz were found between 0.512 and 0.617 N/m. Our calculated results for the interfacial tension of silica nanoparticles within helium atmosphere was 0.676 N/m, which is higher than the value found for the system containing He∕α-quartz (0.478 N/m), but it is similar to the one found for amorphous silica surface. We have also studied the interfacial tension of the nanoparticles in electrolyte aqueous solution for different types and salts concentrations (NaCl, CaCl(2), and MgCl(2)). Our calculations indicate that adsorption properties and salt solutions greatly influence the interfacial tension in an order of CaCl(2) > MgCl(2) > NaCl. This effect is due to the difference in distribution of ions in solution, which modifies the hydration and electrostatic potential of those ions near the nanoparticle.

Structural properties and phase transitions in a silica clathrate.

Melanophlogite, a low-pressure silica polymorph, has been extensively studied at different temperatures and pressures by molecular dynamics simulations. While the high-temperature form is confirmed as cubic, the low-temperature phase is found to be slightly distorted, in agreement with experiments. With increasing pressure, the crystalline character is gradually lost. At 8 GPa, the radial distribution function is consistent with an amorphous state. Like pristine glass, the topology changes, plastic behavior, and permanent densification appear above ∼12 GPa, triggered by Si coordination number changes. We predict that a partial crystalline and amorphous sample can be obtained by recovering the sample from a pressure of ∼12-16 GPa.

Self-accumulation of aromatics at the oil-water interface through weak hydrogen bonding.

It is well-known that the amphiphilic solutes are surface-active and can accumulate at the oil-water interface. Here, we have investigated the water and a light-oil model interface by using molecular dynamic simulations. It was found that aromatics concentrated in the interfacial region, whereas the other hydrocarbons were uniformly distributed throughout the oil phase. Similar to previous studies, such concentrations were not observed at pure aromatics-water interfaces. We show that the self-accumulation of aromatics at the oil-water interface is driven by differences in the interfacial tension, which is lower for aromatics-water than between the others. The weak hydrogen bonding between the aromatic rings and the water protons provides the mechanism for lowering the interfacial tension.

Discovery of Low-Modulus Ti-Nb-Zr Alloys Based on Machine Learning and First-Principles Calculations.

The discovery of low-modulus Ti alloys for biomedical applications is challenging due to a vast number of compositions and available solute contents. In this work, machine learning (ML) methods are employed for the prediction of the bulk modulus (K) and the shear modulus (G) of optimized ternary alloys. As a starting point, the elasticity data of more than 1800 compounds from the Materials Project fed linear models, random forest regressors, and artificial neural networks (NN), with the aims of training predictive models for K and G based on compositional features. The models were then used to predict the resultant Young modulus (E) for all possible compositions in the Ti-Nb-Zr system, with variations in the composition of 2 at. %. Random forest (RF) predictions of E deviate from the NN predictions by less than 4 GPa, which is within the expected variance from the ML training phase. RF regressors seem to generate the most reliable models, given the selected target variables and descriptors. Optimal compositions identified by the ML models were later investigated with the aid of special quasi-random structures (SQSs) and density functional theory (DFT). According to a combined analysis, alloys with 22 Zr (at. %) are promising structural materials to the biomedical field, given their low elastic modulus and elevated beta-phase stability. In alloys with Nb content higher than 14.8 (at. %), the beta phase has lower energy than omega, which may be enough to avoid the formation of omega, a high-modulus phase, during manufacturing.

Brine-Oil Interfacial Tension Modeling: Assessment of Machine Learning Techniques Combined with Molecular Dynamics.

The physical chemistry mechanisms behind the oil-brine interface phenomena are not yet fully clarified. The knowledge of the relation between brine composition and concentration for a given oil may lead to the ionic tuning of the injected solution on geochemical and enhanced oil recovery processes. Thus, it is worth examining the parameters influencing the interfacial properties. In this context, we have combined machine learning (ML) techniques with classical molecular dynamics simulations (MD) to predict oil/brine interfacial tensions (IFT) effectively and compared this process to a linear regression (LR) method. To diversify our data set, we have introduced a new atomistic crude oil model (medium) with 36 different types of hydrocarbon molecules. The MD simulations were performed for mono- and multicomponent (toluene, heptane, Heptol, light, and medium) oil systems interfaced with sulfate and chloride brines with varying cations (Na+, K+, Ca2+, and Mg2+) and salinity concentration. Thus, a consistent IFT data set was built for the ML training and LR fitting at room temperature and pressure conditions, over the feature space considering oil density, oil composition, salinity, and ionic concentrations. On the basis of gradient boosted (GB) algorithms, we have observed that the dominant quantities affecting the IFT are related to the oil attributes and the salinity concentration, and no specific ion dominates the IFT changes. When the obtained LR model was validated against MD and experimental data from the literature, the error varied up to 2% and 9%, respectively, showing a robust and consistent transferability. The combination of MD simulations and ML techniques may provide a fast and cost-effective IFT determination over multiple and complex fluid-fluid and fluid-solid interfaces.

Infrared absorption of MgO at high pressures and temperatures: a molecular dynamic study.

We calculate by molecular dynamics the optical functions of MgO in the far infrared region 100-1000 cm(-1), for pressures up to 40 GPa and temperatures up to 4000 K. An ab initio parametrized many-body force field is used to generate the trajectories. Infrared spectra are obtained from the time correlation of the polarization, and from Kramers-Kronig relations. The calculated spectra agree well with experimental data at ambient pressure. We find that the infrared absorption of MgO at CO(2) laser frequencies increases substantially with both pressure and temperature and we argue that this may explain the underestimation, with respect to theoretical calculations, of the high-pressure melting temperature of MgO determined in CO(2) laser-heated diamond-anvil cell experiments.

Multilevel Molecular Modeling Approach for a Rational Design of Ionic Current Sensors for Nanofluidics.

The complexity displayed by nanofluidic-based systems involves electronic and dynamic aspects occurring across different size and time scales. To properly model such kind of system, we introduced a top-down multilevel approach, combining molecular dynamics simulations (MD) with first-principles electronic transport calculations. The potential of this technique was demonstrated by investigating how the water and ionic flow through a (6,6) carbon nanotube (CNT) influences its electronic transport properties. We showed that the confinement on the CNT favors the partially hydrated Na, Cl, and Li ions to exchange charge with the nanotube. This leads to a change in the electronic transmittance, allowing for the distinguishing of cations from anions. Such an ionic trace may handle an indirect measurement of the ionic current that is recorded as a sensing output. With this case study, we are able to show the potential of this top-down multilevel approach, to be applied on the design of novel nanofluidic devices.

Molecular dynamics studies of aqueous silica nanoparticle dispersions: salt effects on the double layer formation.

The ion distribution around hydroxylated silica nanoparticles (NP-H) dispersed in brine was investigated by fully atomistic molecular dynamics. The NP-H dispersions in aqueous electrolyte media are simulated in solutions of varying salinity (NaCl, CaCl2, and MgCl2), salt concentration (0.06 × 10(-3) to 3.00 × 10(-3) mole fraction [Formula: see text]), and temperature (300 and 350 K) at 1 atm. The NP-H models reproduce the observed experimental concentration of silanol and geminal surface sites, which are responsible for local charge variations on the nanoparticles' surface. Interestingly, under certain salt concentration conditions, the formation of an electrical double layer (DL) around the overall neutral NP-H occurs. The resulting DLs are attenuated with increasing temperature for all evaluated salts. With increasing salt concentration, a sign inversion of the effective charge at the first ion layer is observed, which modifies the electrostatic environment around the nanoparticle. The minimum salt concentration that leads to a DL formation at 300 K is 1.05 × 10(-3), 0.37 × 10(-3), and 0.06 × 10(-3) χs for NaCl, CaCl2, and MgCl2, respectively. The width of the DL decreases sequentially in ionic strength from NaCl to CaCl2 to MgCl2, which is similar to that found for highly charged surfaces. These results are in line with our previous experimental data for negative charged NP-H. All together, these observations suggest an interplay mechanism between the formation and narrowing of electric double layers on the stability of NP dispersions in both neutral and negatively charged NP-H.

First-principles calculations of H, O and OH adsorption on metallic layered supported thin films.

In this work, the adsorption of hydrogen, oxygen and hydroxyl on metallic thin films is studied through first-principles calculations. We explore how the structural and electronic properties of palladium, platinum and gold thin films change with respect to the type of substrate. As a major result we find that Pd/Au(111) and Pt/Au(111) thin films present enhanced adsorption properties for H, O and OH. This improvement is a result of the induced tensile strain on the film due to the misfit between the lattice parameters of the film and the substrate. For these systems, the tensile strain results in a shift of the d-band center position towards to the Fermi level, with implications for the enhancement of adsorption properties. Our results suggest that the location of the unadsorbed d-band center for Pd/Au(111), Pt/Au(111) and Au thin films is a good parameter to predict the reactivity between these surfaces and H, O and OH. However, when considering different numbers of atomic monolayers, changes in adsorption energy are observed and there is no correlation for Pd/Au(111) and Au/Pt(111) films. For Pd/Pt(111) and Pt/Pd(111) films the difference between lattice parameters is relatively small, and no correlation is found, since no considerable strain is induced. In addition, our results support that a compressive strain will always lead to weaker adsorption. We also observe that the work function is strongly affected by adsorption. In particular, H adsorption results in an expansion of the interlayer distance between the topmost layers of the film. Furthermore, after atomic insertion, the interlayer distance of Pd/Pt(111) films is similar to the interlayer distance for bulk PdH0.6, which indicates that these thin films can act as precursor states for hydride formation.

Molecular dynamics studies of fluid/oil interfaces for improved oil recovery processes.

In our paper, we study the interface wettability, diffusivity, and molecular orientation between crude oil and different fluids for applications in improved oil recovery (IOR) processes through atomistic molecular dynamics (MD). The salt concentration, temperature, and pressure effects on the physical chemistry properties of different interfaces between IOR agents [brine (H(2)O + % NaCl), CO(2), N(2), and CH(4)] and crude oil have been determined. From the interfacial density profiles, an accumulation of aromatic molecules near the interface has been observed. In the case of brine interfaced with crude oil, our calculations indicate an increase in the interfacial tension with increasing pressure and salt concentration, which favors oil displacement. On the other hand, with the other fluids studied (CO(2), N(2), and CH(4)), the interfacial tension decreases with increasing pressure and temperature. With interfacial tension reduction, an increase in fluid diffusivity in the oil phase is observed. We also studied the molecular orientation properties of the hydrocarbon and fluids molecules in the interface region. We perceived that the molecular orientation could be affected by changes in the interfacial tension and diffusivity of the molecules in the interface region with the increased pressure and temperature: pressure (increasing) → interfacial tension (decreasing) → diffusion (increasing) → molecular ordering. From a molecular point of view, the combination of low interfacial tension and high diffusion of molecules in the oil phase gives the CO(2) molecules unique properties as an IOR fluid compared with other fluids studied here.

Ab initio investigation of ammonia-borane complexes for hydrogen storage.

The structural, electronic, and thermodynamic properties of ammonia-borane complexes with varying amounts of hydrogen have been characterized by first principles calculations within density functional theory. The calculated structural parameters and thermodynamic functions (free energy, enthalpy and entropy) were found to be in good agreement with experimental and quantum chemistry data for the crystals, dimers, and molecules. The authors find that zero-point energies change several H2 release reactions from endothermic to exothermic. Both the ammonia-borane polymeric and borazine-cyclotriborazane cycles show a strong exothermic decomposition character (approximately -10 kcal/mol), implying that rehydrogenation may be difficult to moderate H2 pressures. Hydrogen bonding in these systems has been characterized and they find the N-H bond to be more covalent than the more ionic B-H bond.

Tuning oxygen packing in silica by nonhydrostatic pressure.

The transformation of SiO2 from low pressure tetrahedral phases into denser octahedral phases takes place via the collapse of the oxygen sublattice into a close-packed arrangement. The transition paths and the resulting products are known to be affected by the presence of anisotropic stresses, which are difficult to control, so interpretation of the experimental results is problematic. Based on nonhydrostatic molecular dynamics simulations, we show that the collapse of the oxygen sublattice in the specific case of cristobalite is concomitant with the disappearance of tetrahedral units and that non hydrostatic stresses can be tuned to yield phases with different oxygen close-packed sublattices, including the alpha-PbO2-like phase, for which we provide a microscopic formation path, and phases with a cubic close packing, like anatase, not seen in experiments yet.

Infrared and Raman spectra of silica polymorphs from an ab initio parametrized polarizable force field.

The general aim of this study is to test the reliability of polarizable model potentials for the prediction of vibrational (infrared and Raman) spectra in highly anharmonic systems such as high temperature crystalline phases. By using an ab initio parametrized interatomic potential for SiO2 and molecular dynamics simulations, we calculate the infrared and Raman spectra for quartz, cristobalite, and stishovite at various thermodynamic conditions. The model is found to perform very well in the prediction of infrared spectra. Raman peak positions are also reproduced very well by the model; however, Raman intensities calculated by explicitly taking the derivative of the polarizability with respect to the atomic displacements are found to be in poorer agreement than intensities calculated using a parametrized "bond polarizability" model. Calculated spectra for the high temperature beta phases, where the role of dynamical disorder and anharmonicities is predominant, are found to be in excellent agreement with experiments. For the octahedral phases, our simulations are able to reproduce changes in the Raman spectra across the rutile-to-CaCl2 transition around 50 GPa, including the observed phonon softening.