The Influence of Ni Addition in the Mechanism of CO2 Electroreduction on Cu Crystals-Mechanistic Insight from DFT Simulations.

We present a DFT analysis of the role of the Cu-Ni synergistic effect for the CO2 reduction to C2H4, in comparison to the pure Cu catalyst. The analysis is focused on the thermodynamic stability of reactive intermediates along the proposed pathway of C2 species formation. We have observed that the potential needed for the reaction decreases with the addition of Ni in the investigated model. In addition, we have observed the differences in the preferred pathway based on the significant differences in stability of the reactive intermediates depending on th Cu:Ni ratio. The results suggest that despite the fact the Cu surface is always exposed, and it is the only one that is able to directly interact with the intermediates, the presence of the Ni in the underlying sections of the crystal is significant enough to change the mechanism of the reaction.

Reactive Grand-Canonical Monte Carlo Simulations for Modeling Hydration of MgCl2.

Thermochemical heat-storage applications, based on the reversible endo-/exothermic hydration reaction of salts, are intensively investigated to search for compact heat-storage devices. To achieve a truly valuable storage system, progressively complex salts are investigated. For these salts, the equilibrium temperature and pressure conditions are not always easy to predict. However, these conditions are crucial for the design of thermochemical heat-storage systems. A biased grand-canonical Monte Carlo (GCMC) tool is developed, enabling the study of equilibrium conditions at the molecular level. The GCMC algorithm is combined with reactive force field molecular dynamics (ReaxFF), which allows bond formation within the simulation. The Weeks-Chandler-Andersen (WCA) potential is used to scan multiple trial positions for the GCMC algorithm at a small cost. The most promising trial positions can be selected for recomputation with the more expensive ReaxFF. The developed WCA-ReaxFF-GCMC tool was used to study the hydration of MgCl2·nH2O. The simulation results show a good agreement with experimental and thermodynamic equilibriums for multiple hydration levels. The hydration shows that water, present at the surface of crystalline salt, deforms the surface layers and promotes further hydration of these deformed layers. Additionally, the WCA-ReaxFF-GCMC algorithm can be used to study other, non-TCM-related, reactive sorption processes.

Gibbs Ensemble Monte Carlo for Reactive Force Fields to Determine the Vapor-Liquid Equilibrium of CO2 and H2O.

Absorption and reactive properties of fluids in porous media are key to the design and improvement of numerous energy related applications. Molecular simulations of these systems require accurate force fields that capture the involved chemical reactions and have the ability to describe the vapor-liquid equilibrium (VLE). Two new reactive force fields (ReaxFF) for CO2 and H2O are developed, which are capable of not only modeling bond breaking and formation in reactive environments but also predicting their VLEs at saturation conditions. These new force fields include extra terms (ReaxFF-lg) to improve the long-range interactions between the molecules. For validation, we have developed a new Gibbs ensemble Monte Carlo (GEMC-ReaxFF) approach to predict the VLE. Computed VLE data show good agreement with National Institute of Standards and Technology reference data as well as existing nonreactive force fields. This validation proves the applicability of the GEMC-ReaxFF method to test new reactive force fields, and simultaneously it proves the applicability to extend newly developed ReaxFF force fields to other more complex reactive systems.

Absolute energy level positions in tin- and lead-based halide perovskites.

Metal halide perovskites are promising materials for future optoelectronic applications. One intriguing property, important for many applications, is the tunability of the band gap via compositional engineering. While experimental reports on changes in absorption or photoluminescence show rather good agreement for different compounds, the physical origins of these changes, namely the variations in valence and conduction band positions, are not well characterized. Here, we determine ionization energy and electron affinity values of all primary tin- and lead-based perovskites using photoelectron spectroscopy data, supported by first-principles calculations and a tight-binding analysis. We demonstrate energy level variations are primarily determined by the relative positions of the atomic energy levels of metal cations and halide anions and secondarily influenced by the cation-anion interaction strength. These results mark a significant step towards understanding the electronic structure of this material class and provides the basis for rational design rules regarding the energetics in perovskite optoelectronics.

Correlations between Density-Based Bond Orders and Orbital-Based Bond Energies for Chemical Bonding Analysis.

Quantum chemistry-based codes and methods provide valuable computational tools to estimate reaction energetics and elucidate reaction mechanisms. Electronic structure methods allow directly studying the chemical transformations in molecular systems involving breaking and making of chemical bonds and the associated changes in the electronic structure. The link between the electronic structure and chemical bonding can be provided through the crystal orbital Hamilton population (COHP) analysis that allows quantifying the bond strength by computing Hamilton-weighted populations of localized atomic orbitals. Another important parameter reflecting the nature and strength of a chemical bond is the bond order that can be assessed by the density derived electrostatic and chemical (DDEC6) method which relies on an electron and spin density-partitioning scheme. Herein, we describe a linear correlation that can be established between the DDEC6-derived bond orders and the bond strengths computed with the COHP formalism. We demonstrate that within defined boundaries, the COHP-derived bond strengths can be consistently compared among each other and linked to the DDEC6-derived bond orders independent of the used model. The validity of these correlations and the effective model independence of the electronic descriptors are demonstrated for a variety of gas-phase chemical systems, featuring different types of chemical bonds. Furthermore, the applicability of the derived correlations to the description of complex reaction paths in periodic systems is demonstrated by considering the zeolite-catalyzed Diels-Alder cycloaddition reaction between 2,5-dimethylfuran and ethylene.

Mechanistic Insight into the [4 + 2] Diels-Alder Cycloaddition over First Row d-Block Cation-Exchanged Faujasites.

The Diels-Alder cycloaddition (DAC) is a powerful tool to construct C-C bonds. The DAC reaction can be accelerated in several ways, one of which is reactant confinement as observed in supramolecular complexes and Diels-Alderases. Another method is altering the frontier molecular orbitals (FMOs) of the reactants by using homogeneous transition-metal complexes whose active sites exhibit d-orbitals suitable for net-bonding orbital interactions with the substrates. Both features can be combined in first row d-block (TM) exchanged faujasite catalysts where the zeolite framework acts as a stabilizing ligand for the active site while confining the reactants. Herein, we report on a mechanistic and periodic DFT study on TM-(Cu(I), Cu(II), Zn(II), Ni(II), Cr(III), Sc(III), V(V))exchanged faujasites to elucidate the effect of d-shell filling on the DAC reaction between 2,5-dimethylfuran and ethylene. Two pathways were found: one being the concerted one-step and the other being the stepwise two-step pathway. A decrease in d-shell filling results in a concomitant increase in reactant activation as evidenced by increasingly narrow energy gaps and lower activation barriers. For models holding relatively small d-block cations, the zeolite framework was found to bias the DAC reaction toward an asynchronous one-step pathway instead of the two-step pathway. This work is an example of how the active site properties and the surrounding chemical environment influence the reaction mechanism of chemical transformations.

Tin, The Enabler-Hydrogen Diffusion into Ruthenium.

Hydrogen interaction with ruthenium is of particular importance for the ruthenium-capped multilayer reflectors used in extreme ultraviolet (EUV) lithography. Hydrogen causes blistering, which leads to a loss of reflectivity. This problem is aggravated by tin. This study aims to uncover the mechanism via which tin affects the hydrogen uptake, with a view to mitigation. We report here the results of a study of hydrogen interaction with the ruthenium surface in the presence of tin using Density Functional Theory and charge density analyses. Our calculations show a significant drop in the energy barrier to hydrogen penetration when a tin atom or a tin hydride molecule (SnHx) is adsorbed on the ruthenium surface; the barrier has been found to drop in all tested cases with tin, from 1.06 eV to as low as 0.28 eV in the case of stannane (SnH4). Analyses show that, due to charge transfer from the less electronegative tin to hydrogen and ruthenium, charge accumulates around the diffusing hydrogen atom and near the ruthenium surface atoms. The reduced atomic volume of hydrogen, together with the effect of electron⁻electron repulsion from the ruthenium surface charge, facilitates subsurface penetration. Understanding the nature of tin's influence on hydrogen penetration will guide efforts to mitigate blistering damage of EUV optics. It also holds great interest for applications where hydrogen penetration is desirable, such as hydrogen storage.

Electronic Structure Analysis of the Diels-Alder Cycloaddition Catalyzed by Alkali-Exchanged Faujasites.

The Diels-Alder cycloaddition (DAC) reaction is a commonly employed reaction for the formation of C-C bonds. DAC catalysis can be achieved by using Lewis acids and via reactant confinement in aqueous nanocages. Low-silica alkali-exchanged faujasite catalysts combine these two factors in one material. They can be used in the tandem DAC/dehydration reaction of biomass-derived 2,5-dimethylfuran (DMF) with ethylene toward p-xylene, in which the DAC reaction step initiates the overall reaction cycle. In this work, we performed periodic density functional theory (DFT) calculations on the DAC reaction between DMF and C2H4 in low-silica alkali(M)-exchanged faujasites (MY; Si/Al = 2.4; M = Li+, Na+, K+, Rb+, Cs+). The aim was to investigate how confinement of reactants in MY catalysts changed their electronic structure and the DAC-reactivity trend among the evaluated MY zeolites. The conventional high-silica alkali-exchanged isolated site model (MFAU; Si/Al = 47) served as a reference. The results show that confinement leads to initial-state (IS) destabilization and transition-state (TS) stabilization. Among the tested MY, most significant IS destabilization is found in RbY. Only antibonding orbital interactions between the reactants/reactive complex and cations were found, indicating that TS stabilization arises from ionic interactions. Additionally, in RbY the geometry of the transition state is geometrically most similar to that of the initial and final state. RbY also exhibits an optimal combination of the confinement-effects, resulting in having the lowest computed DAC-activation energy. The overall effect is a DAC-reactivity trend inversion in MY as compared to the trend found in MFAU where the activation energy correlates with the Lewis acidity of the exchangeable cations.

An Active Alkali-Exchanged Faujasite Catalyst for p-Xylene Production via the One-Pot Diels-Alder Cycloaddition/Dehydration Reaction of 2,5-Dimethylfuran with Ethylene.

The one-pot Diels-Alder cycloaddition (DAC)/dehydration (D) tandem reaction between 2,5-dimethylfuran and ethylene is a potent pathway toward biomass-derived p-xylene. In this work, we present a cheap and active low-silica potassium-exchanged faujasite (KY, Si/Al = 2.6) catalyst. Catalyst optimization was guided by a computational study of the DAC/D reaction mechanism over different alkali-exchanged faujasites using periodic density functional theory calculations complemented by microkinetic modeling. Two types of faujasite models were compared, i.e., a high-silica alkali-exchanged faujasite model representing isolated active cation sites and a low-silica alkali-exchanged faujasite in which the reaction involves several cations in the proximity. The mechanistic study points to a significant synergetic cooperative effect of the ensemble of cations in the faujasite supercage on the DAC/D reaction. Alignment of the reactants by their interactions with the cationic sites and stabilization of reaction intermediates contribute to the high catalytic performance. Experiments confirmed the prediction that KY is the most active catalyst among low-silica alkali-exchanged faujasites. This work is an example of how the catalytic reactivity of zeolites depends on multiple interactions between the zeolite and reagents.

Elucidation of the origin of chiral amplification in discrete molecular polyhedra.

Chiral amplification in molecular self-assembly has profound impact on the recognition and separation of chiroptical materials, biomolecules, and pharmaceuticals. An understanding of how to control this phenomenon is nonetheless restricted by the structural complexity in multicomponent self-assembling systems. Here, we create chiral octahedra incorporating a combination of chiral and achiral vertices and show that their discrete nature makes these octahedra an ideal platform for in-depth investigation of chiral transfer. Through the construction of dynamic combinatorial libraries, the unique possibility to separate and characterise each individual assembly type, density functional theory calculations, and a theoretical equilibrium model, we elucidate that a single chiral unit suffices to control all other units in an octahedron and how this local amplification combined with the distribution of distinct assembly types culminates in the observed overall chiral amplification in the system. Our combined experimental and theoretical strategy can be applied generally to quantify discrete multi-component self-assembling systems.

Photoisomerization induced scission of rod-like micelles unravelled with multiscale modeling.

HYPOTHESIS: In photorheological fluids, subtle molecular changes caused by light lead to abrupt macroscopic alterations. Upon UV irradiation of an aqueous cetyltrimethylammonium bromide (CTAB) and trans-ortho-methoxycinnamic acid (trans-OMCA) solution, for instance, the viscosity drops over orders of magnitude. Multiscale modeling allows to elucidate the mechanisms behind these photorheological effects. EXPERIMENTS: We use time-dependent DFT calculations to study the photoisomerization, and a combination of atomistic molecular dynamics (MD) and DFT to probe the influence of both OMCA isomers on the micellar solutions. FINDINGS: The time-dependent DFT calculations show that the isomerization pathway occurs in the first triplet excited state with a minimum energy conformation closest to the after photoisomerization predominant cis configuration. In the MD simulations, with sub-microsecond timescales much shorter than the experimental morphological transition, already a clear difference is observed in the packing of the two OMCA isomers: contrary to trans-OMCA, cis-OMCA exposes notable part of its hydrophobic aromatic rings at the micelle surface. This can explain why trans-OMCA adopts rod-like micellar packing (high viscosity) while cis-OMCA spherical micellar packing (low viscosity). Moreover, lowering of the OMCA co-solute concentration allowed us to perform full simulation of the breakup process of the rod-like micelles which are stable prior to isomerization.

Scaling Relations for Acidity and Reactivity of Zeolites.

Zeolites are widely applied as solid acid catalysts in various technological processes. In this work we have computationally investigated how catalytic reactivity scales with acidity for a range of zeolites with different topologies and chemical compositions. We found that straightforward correlations are limited to zeolites with the same topology. The adsorption energies of bases such as carbon monoxide (CO), acetonitrile (CH3CN), ammonia (NH3), trimethylamine (N(CH3)3), and pyridine (C5H5N) give the same trend of acid strength for FAU zeolites with varying composition. Crystal orbital Hamilton populations (COHP) analysis provides a detailed molecular orbital picture of adsorbed base molecules on the Brønsted acid sites (BAS). Bonding is dominated by strong σ donation from guest molecules to the BAS for the adsorbed CO and CH3CN complexes. An electronic descriptor of acid strength is constructed based on the bond order calculations, which is an intrinsic parameter rather than adsorption energy that contains additional contributions due to secondary effects such as van der Waals interactions with the zeolite walls. The bond order parameter derived for the CH3CN adsorption complex represents a useful descriptor for the intrinsic acid strength of FAU zeolites. For FAU zeolites the activation energy for the conversion of π-adsorbed isobutene into alkoxy species correlates well with the acid strength determined by the NH3 adsorption energies. Other zeolites such as MFI and CHA do not follow the scaling relations obtained for FAU; we ascribe this to the different van der Waals interactions and steric effects induced by zeolite framework topology.

Role of Adsorbed Water on Charge Carrier Dynamics in Photoexcited TiO2.

Overall photocatalytic water splitting is one of the most sought after processes for sustainable solar-to-chemical energy conversion. The efficiency of this process strongly depends on charge carrier recombination and interaction with surface adsorbates at different time scales. Here, we investigated how hydration of TiO2 P25 affects dynamics of photogenerated electrons at the millisecond to minute time scale characteristic for chemical reactions. We used rapid scan diffuse-reflectance infrared Fourier transform spectroscopy (DRIFTS). The decay of photogenerated electron absorption was substantially slower in the presence of associated water. For hydrated samples, the charge carrier recombination rates followed an Arrhenius-type behavior in the temperature range of 273-423 K; these became temperature-independent when the material was dehydrated at temperatures above 423 K or cooled below 273 K. A DFT+U analysis revealed that hydrogen bonding with adsorbed water stabilizes surface-trapped holes at anatase TiO2(101) facet and lowers the barriers for hole migration. Hence, hole mobility should be higher in the hydrated material than in the dehydrated system. This demonstrates that adsorbed associated water can efficiently stabilize photogenerated charge carriers in nanocrystalline TiO2 and suppress their recombination at the time scale up to minutes.

Structure of palladium nanoparticles under oxidative conditions.

Using density functional theory (DFT) and thermodynamic considerations we study the shape and stability of Pd nanoparticles in oxygen-lean and oxygen-rich atmospheres. We find that at very high oxygen coverage cubes exposing (100) faces will form, which are stabilized due to the formation of a O/(√5 × âˆš5)R27° overlayer. The shape of oxygen-covered Pd and Pt nanoparticles is compared in this study.

Unraveling the self-assembly of the (S)-glutamic acid "flower" structure on Ag(100).

(S)-Glutamic acid adsorbed on Ag(100) organizes in different self-assembled structures depending on surface temperature [Smerieri, M.; Vattuone, L.; Kravchuk, T.; Costa, D.; Savio, L. (S)-Glutamic Acid on Ag(100): Self-Assembly in the Nonzwitterionic Form. Langmuir2011, 27, 2393-2404]. In particular, two of these structures, referred to as "square" and "flower" geometries, are found to coexist on the surface upon deposition at T = 350 K. The former assembly was fully resolved at the atomic level in the work of Smerieri et al. [Smerieri, M.; Vattuone, L.; Costa, D.; Tielens, F.; Savio, L. Self-Assembly of (S)-Glutamic Acid on Ag(100): A Combined LT-STM and Ab Initio Investigation. Langmuir2010, 26, 7208-7215], in which we proved that the driving force for adsorption is the van der Waals interactions between the molecules and the Ag surface, that is, that molecules are in a physisorbed state. In this paper, we complete our work by presenting the characterization of the "flower" structure. In contrast to the case of the "square" assembly, a strong chemical bond between glutamic acid radicals and the surface is at the basis of the "flowers" geometry. Whereas the chemisorbed central GLU tetramer interacts strongly with the surface, the physisorbed surrounding GLU molecules conserve some degree of freedom in the layer which counterbalances the weak adsorption energy. The "flower" and the "square" assemblies have similar dispersion energy and H-bond interaction energy; as a consequence of the different chemical state of the GLU molecules, however, such contributions have a very different relative weight in the stabilization of the two structures.