Modeling the polychromism of oxide minerals: The case of alexandrite and cordierite.

In this work, we investigate the spectroscopic properties of photochromic alexandrite and cordierite by TD-DFT. The objective is to assess the TD-DFT for the simulation of pleochroism (change of color depending on the crystallographic direction of the observation) and the change of color as a function of the light source. For these simulations, we compared an embedding where dangling bonds are saturated by hydrogen atoms and an electrostatic embedding. The electrostatic embedding provided numerically more stable results and allowed a good reproduction of the pleochroism of cordierite, based on a Fe2+-Fe3+ intervalence charge transfer transition. However, the pleochroism of alexandrite is not as well reproduced, suggesting that TD-DFT has some difficulties to reproduce the anisotropy of the transition dipole moment, an aspect that is not deeply documented in the literature.

Unraveling the Oxidation Kinetics Through Electronic Structure Regulation of MnCo2O4.5@Ni3S2 p-n Junction for Urea-Assisted Electrocatalytic Activity.

A promising strategy to boost electrocatalytic performance is via assembly of hetero-nanostructured electrocatalysts that delivers the essential specific surface area and also active sites by lowering the reaction barrier. However, the challenges associated with the intricate designs and mechanisms remain underexplored. Therefore, the present study constructs a p-n junction in a free-standing MnCo2O4.5@Ni3S2 on Ni-Foam. The space-charge region's electrical characteristics is dramatically altered by the formed p-n junction, which enhances the electron transfer process for urea-assisted electrocatalytic water splitting (UOR). The optimal MnCo2O4.5@Ni3S2 electrocatalyst results in greater oxygen evolution reactivity (OER) than pure systems, delivering an overpotential of only 240 mV. Remarkably, upon employing as UOR electrode the required potential decreases to 30 mV. The impressive performance of the designed catalyst is attributed to the enhanced electrical conductivity, greater number of electrochemical active sites, and improved redox activity due to the junction interface formed between p-MnCo2O4.5 and n-Ni3S2. There are strong indications that the in situ formed extreme-surface NiOOH, starting from Ni3S2, boosts the electrocatalytic activity, i.e., the electrochemical  surface reconstruction generates the active species. In conclusion, this work presents a high-performance p-n junction design for broad use, together with a viable and affordable UOR electrocatalyst.

One-Step Electrocatalytic Approach Applied to the Synthesis of ß-Propiolactones from CO2 and Dienes.

ß-Lactones are common substructures in a variety of natural products and drugs, and they serve as versatile synthetic intermediates in the production of valuable chemical derivatives. Traditional ß-lactone synthesis relies on laborious multi-step synthetic methods that use toxic compounds, sophisticated catalysts, expensive, and/or reactive chemicals. Based on the in situ electrochemical formation of metal-based nanoclusters, this paper describes the development of a one-step, room temperature electrocatalytic method for the formation of stable ß-lactone from CO2 and dienes. This one-step "electrosynthesis" method results in the formation of a new class of ß-lactone with high selectivity (up to 100%) and activity (up to 80% yields with respect to the reacted diene) by regulating the applied potential and current density. This work paves the way for more sustainable and environmentally friendly reaction pathways based on the in situ formation of nanoclusters as organic electrosynthesis catalysts.

Gaussian attractive potential for carboxylate/cobalt surface interactions.

Ligand-decorated metal surfaces play a pivotal role in various areas of chemistry, particularly in selective catalysis. Molecular dynamics simulations at the molecular mechanics level of theory are best adapted to gain complementary insights to experiments regarding the structure and dynamics of such organic films. However, standard force fields tend to capture only weak physisorption interactions. This is inadequate for ligands that are strongly adsorbed such as carboxylates on metal surfaces. To address this limitation, we employ the Gaussian Lennard-Jones (GLJ) potential, which incorporates an attractive Gaussian potential between the surface and ligand atoms. Here, we develop this approach for the interaction between cobalt surfaces and carboxylate ligands. The accuracy of the GLJ approach is validated through the analysis of the interaction of oxygen with two distinct cobalt surfaces. The accuracy of this method reaches a root mean square deviation (RMSD) of about 3 kcal/mol across all probed configurations, which corresponds to a percentage error of roughly 4%. Application of the GLJ force field to the dynamics of the organic layer on these surfaces reveals how the ligand concentration influences the film order, and highlights differing mobility in the x and y directions, attributable to surface corrugation on Co(112̄0). GLJ is versatile, suitable for a broad range of metal/ligand systems, and can, subsequently, be utilized to study the organic film on the adsorption/desorption of reactants and products during a catalytic process.

Beyond single-crystal surfaces: The GAL21 water/metal force field.

Solvent effects are notoriously difficult to describe for metallic nanoparticles (NPs). Here, we introduce GAL21 which is the first pairwise additive force field that is specifically designed to modulate the near chemisorption energy of water as a function of the coordination numbers of the metallic atoms. We find a quadratic dependence to be most suitable for capturing the dependence of the adsorption energy of water on the generalized coordination number (GCN) of the metal atoms. GAL21 has been fitted against DFT adsorption energies for Cu, Ag, Au, Ni, Pd, Pt, and Co on 500 configurations and validated on about 3000 configurations for each metal, constructed on five surfaces with GCNs varying from 2.5 to 11.25. Depending on the metals, the root mean square deviation is found between 0.7 kcal mol-1 (Au) to 1.6 kcal mol-1 (Ni). Using GAL21, as implemented in the open-source code CP2K, we then evaluate the solvation energy of Au55 and Pt55 NPs in water using thermodynamic integration. The solvation free energy is found to be larger for Pt than for Au and systematically larger than 200 kcal mol-1, demonstrating the large impact of solvent on the surface energetics of NPs. Still, given that the amorphous NPs are both, the most stable and the most solvated ones, we do not predict a change in the preferred morphology between the gas-phase and in water. Finally, based on a linear regression on three sizes of NPs (from 38 to 147), the solvation energy for Au and Pt surface atoms is found to be -5.2 and -9.9 kcal mol-1, respectively.

DockOnSurf: A Python Code for the High-Throughput Screening of Flexible Molecules Adsorbed on Surfaces.

We present the open-source python package DockOnSurf which automates the generation and optimization of low-energy adsorption configurations of molecules on extended surfaces and nanoparticles. DockOnSurf is especially geared toward handling polyfunctional flexible adsorbates. The use of this high-throughput workflow allows us to carry out the screening of adsorbate-surface configurations in a systematic, customizable, and traceable way, while keeping the focus on the chemically relevant structures. The screening strategy consists in splitting the exploration of the adsorbate-surface configurational space into chemically meaningful domains, that is, by choosing among different conformers to adsorb, surface adsorption sites, adsorbate anchoring points, and orientations and allowing dissociation of (acidic) protons. We demonstrate the performance of the main features based on varying examples, ranging from CO adsorption on a gold nanoparticle to sorbitol adsorption on hematite. Through the use of the presented program, we aim to foster efficiency, traceability, and ease of use in research within tribology, catalysis, nanoscience, and surface science in general.

What does graphitic carbon nitride really look like?

Graphitic carbon nitrides (g-CNs) have become popular light absorbers in photocatalytic water splitting cells. Early theoretical work on these structures focused on fully polymerized g-C3N4. Experimentally, it is known that the typically employed melamine polycondensation does not go toward completion, yielding structures with ∼15 at% hydrogen. Here, we study the conformational stability of "melon", with the [C6N9H3]n structural formula using DFT. Referencing to a 2D melon sheet, B3LYP-dDsC and PBE-MBD computations revealed the same qualitative trend in stability of the 3D structures, with several of them within 5 kJ mol-1 per tecton. Fina's orthorhombic melon is the most stable of the studied conformers, with Lotsch' monoclinic melon taking an intermediate value. Invoking a simple Wannier-Mott-type approach, Fina's and Lotsch' structures exhibited the lowest optical gaps (2.8 eV), within the error margin of the experimental value (2.7 eV). All conformers yielded gaps below that of the monolayer's (3.2 eV), suggesting Jelley-type ("J") aggregation effects.

(Dis)Similarities of adsorption of diverse functional groups over alumina and hematite depending on the surface state.

To accelerate the conversion to more sustainable lubricants, there is a need for an improved understanding of the adsorption at the solid/liquid interface. As a first step, the density functional theory computed adsorption energies can be used to screen the ability of additives to cover a surface. Analogously to what has been found in catalysis with the universal scaling relations, we investigate here if a general universal ranking of additives can be found, independently of the surface considered. We divided our set of 25 diverse representative molecules into aprotic and protic molecules. We compared their adsorption over alumina and hematite, which are models of surface oxidized aluminum and steel, respectively. The adsorption energy ranking of our set is not strongly affected by alumina hydration. In contrast, adsorption on hematite is more strongly affected by hydration since all exposed Fe Lewis acid sites are converted into hydroxylated Brønsted basic sites. However, the ranking obtained on hydrated hematite is close to the one obtained on dry alumina, paving the road to a fast screening of additives. In our library, protic molecules are more strongly adsorbed than non-protic molecules. In particular, methyl and dimethyl phosphates are the most strongly adsorbed ones, followed by N-methyldiethanolamine, succinimide, and ethanoic acid. Additives combining these functional groups are expected to strongly adsorb at the solid/liquid interface and, therefore, likely to be relevant components of lubricant formulations.

Revisiting the Atomistic Structures at the Interface of Au(111) Electrode-Sulfuric Acid Solution.

Knowledge of atomistic structures at solid/liquid interfaces is essential to elucidate interfacial processes in chemistry, physics, and materials sciences. The (√3 × âˆš7) structure associated with a pair of sharp reversible current spikes in the cyclic voltammogram on a Au(111) electrode in sulfuric acid solution represents one of the most classical ordered structures at electrode/electrolyte interfaces. Although more than 10 adsorption configurations have been proposed in the past four decades, the atomistic structure remains ambiguous and is consequently an open problem in electrochemistry and surface science. Herein, by combining high-resolution electrochemical scanning tuning microscopy, electrochemical infrared and Raman spectroscopies, and, in particular, the newly developed quantitative computational method for electrochemical infrared and Raman spectra, we unambiguously reveal that the adstructure is Au(111)(√3 × âˆš7)-(SO4···w2) with a sulfate anion (SO4*) and two structured water molecules (w2*) in a unit cell, and the crisscrossed [w···SO4···w]n and [w···w···]n hydrogen-bonding network comprises the symmetric adstructure. We further elucidate that the electrostatic potential energy dictates the proton affinity of sulfate anions, leading to the potential-tuned structural transformations. Our work enlightens the structural details of the inner Helmholtz plane and thus advances our fundamental understanding of the processes at electrochemical interfaces.

Two-sites are better than one: revisiting the OER mechanism on CoOOH by DFT with electrode polarization.

We uncover the existence of several competitive mechanisms of water oxidation on the ß-CoOOH (10-14) surface by going beyond the classical 4-step mechanism frequently used to study this reaction at the DFT level. Our results demonstrate the importance of two-site reactivity and of purely chemical steps with the associated activation energies. Taking the electrochemical potential explicitly into account leads to modifications of the reaction energy profiles finally leading to the proposition of a new family of mechanisms involving tetraoxidane intermediates. The two-site mechanisms revealed in this work are of key importance to rationalize and predict the impact of dopants in the design of future catalysts.

Parameter-free coordination numbers for solutions and interfaces.

Coordination numbers are among the central quantities to describe the local environment of atoms and are thus used in various applications such as structure analysis, fingerprints, and parameters. Yet, there is no consensus regarding a practical algorithm, and many proposed methods are designed for specific systems. In this work, we propose a scale-free and parameter-free algorithm for nearest neighbor identification. This algorithm extends the powerful Solid-Angle based Nearest-Neighbor (SANN) framework to explicitly include local anisotropy. As such, our Anisotropically corrected SANN (ASANN) algorithm provides with a fast, robust, and adaptive method for computing coordination numbers. The ASANN algorithm is applied to flat and corrugated metallic surfaces to demonstrate that the expected coordination numbers are retrieved without the need for any system-specific adjustments. The same applies to the description of the coordination numbers of metal atoms in AuCu nanoparticles, and we show that ASANN based coordination numbers are well adapted for automatically counting neighbors and the establishment of cluster expansions. Analysis of classical molecular dynamics simulations of an electrified graphite electrode reveals a strong link between the coordination number of Cs+ ions and their position within the double layer, a relation that is absent for Na+, which keeps its first solvation shell even close to the electrode.

Water adlayers on noble metal surfaces: Insights from energy decomposition analysis.

Water molecules adsorbed on noble metal surfaces are of fundamental interest in surface science, in heterogeneous catalysis, and as a model for the metal/water interface. Herein, we analyze 28 water structures adsorbed on five noble metal surfaces (Cu, Ag, Au, Pd, and Pt) via density functional theory and energy decomposition analysis based on the block localized wave function technique. Structures, ranging from monomers to ice adlayers, reveal that the charge transfer from water to the surface is nearly independent from the charge transfer between the water molecules, while the polarization energies are cooperative. Dense water-water networks with small surface dipoles, such as the 39×39 unit cell [experimentally observed on Pt(111)], are favored compared to the highly ordered and popular Hup and Hdown phases. The second main result of our study is that the many-body interactions, which stabilize the water assemblies on the metal surfaces, are dominated by the polarization energies, with the charge transfer scaling with the polarization energies. Hence, if an empirical model could be found that reproduces the polarization energies, the charge transfer could be predicted as well, opening exciting perspectives for force field development.

Can microsolvation effects be estimated from vacuum computations? A case-study of alcohol decomposition at the H2O/Pt(111) interface.

Converting biomass into sustainable chemicals and energy feed-stocks requires innovative heterogeneous catalysts, which are able to efficiently work under aqueous conditions. Computational chemistry is a key asset in the design of these novel catalysts, but it has to face two challenges: the large reaction networks and the potential role of hydration. They can be addressed using scaling relations such as Brønsted-Evans-Polanyi (BEP) and solvation models, respectively. In this study, we show that typical reaction and activation energies of alcohol decomposition on Pt(111) are not strongly modified by the inclusion of the water solvent as a continuum model. In contrast, adding a single water molecule strongly favours O-H and C-OH scission and it prevents C-O and to a lesser extent C-C scissions. The resulting BEP relationships partially reflect these changes induced by the solvent. Predicting Pt-catalysed alcohol decomposition in water should thus account for the influence of the solvent on thermodynamics and kinetics. In addition, we found that the reaction energies obtained in the presence of an explicit water molecule scale with the ones obtained in a vacuum. Hence, we reveal that vacuum computations in combination with corrections based on our linear regressions are able to capture the important H-bonding effect.

Implicit self-consistent electrolyte model in plane-wave density-functional theory.

The ab initio computational treatment of electrochemical systems requires an appropriate treatment of the solid/liquid interfaces. A fully quantum mechanical treatment of the interface is computationally demanding due to the large number of degrees of freedom involved. In this work, we develop a computationally efficient model where the electrode part of the interface is described at the density-functional theory (DFT) level, and the electrolyte part is represented through an implicit solvation model based on the Poisson-Boltzmann equation. We describe the implementation of the linearized Poisson-Boltzmann equation into the Vienna Ab initio Simulation Package, a widely used DFT code, followed by validation and benchmarking of the method. To demonstrate the utility of the implicit electrolyte model, we apply it to study the surface energy of Cu crystal facets in an aqueous electrolyte as a function of applied electric potential. We show that the applied potential enables the control of the shape of nanocrystals from an octahedral to a truncated octahedral morphology with increasing potential.

Molecular mechanics models for the image charge, a comment on "including image charge effects in the molecular dynamics simulations of molecules on metal surfaces".

We re-investigate the image charge model of Iori and Corni (Iori and Corni, J. Comput. Chem. 2008, 29, 1656). We find that a simple symmetrization of their model allows to obtain quantitatively correct results for the electrostatic interaction of a water molecule with a metallic surface. This symmetrization reduces the magnitude of the electrostatic interaction to less than 10% of the total interaction energy. © 2017 Wiley Periodicals, Inc.

A machine learning approach to graph-theoretical cluster expansions of the energy of adsorbate layers.

The accurate description of the energy of adsorbate layers is crucial for the understanding of chemistry at interfaces. For heterogeneous catalysis, not only the interaction of the adsorbate with the surface but also the adsorbate-adsorbate lateral interactions significantly affect the activation energies of reactions. Modeling the interactions of the adsorbates with the catalyst surface and with each other can be efficiently achieved in the cluster expansion Hamiltonian formalism, which has recently been implemented in a graph-theoretical kinetic Monte Carlo (kMC) scheme to describe multi-dentate species. Automating the development of the cluster expansion Hamiltonians for catalytic systems is challenging and requires the mapping of adsorbate configurations for extended adsorbates onto a graphical lattice. The current work adopts machine learning methods to reach this goal. Clusters are automatically detected based on formalized, but intuitive chemical concepts. The corresponding energy coefficients for the cluster expansion are calculated by an inversion scheme. The potential of this method is demonstrated for the example of ethylene adsorption on Pd(111), for which we propose several expansions, depending on the graphical lattice. It turns out that for this system, the best description is obtained as a combination of single molecule patterns and a few coupling terms accounting for lateral interactions.

Solvation free energies for periodic surfaces: comparison of implicit and explicit solvation models.

The evaluation of solvation energies is a great challenge. We focus here on an organic molecule chemisorbed at a metal-liquid interface, as a prototypical system, essential in tribology, electrochemistry and heterogeneous catalysis. We compare an established implicit solvation scheme with a strategy based on molecular mechanics (MM) free energy perturbation (FEP) seeded by QM computations. First, we benchmark the approaches against experimental hydration energies of standard (organic) molecules and find acceptable errors in the order of 0.06 eV (1.3 kcal mol-1). Then, the impact of various parameters on the solvation energy of an adsorbate have been assessed on a typical system of interest, levulinic acid adsorbed at a Ru(0001)/water interface. We identify the need for dipole corrections or symmetric slabs when including solvation effects on metallic surfaces. The MM-FEP scheme is revealed to be as reliable as the implicit solvent for water. In the case of levulinic acid, both PCM and MM-FEP agree that the bulk solvation effect is not sufficient to change the adsorption mode from bidentate to mono-dentate, despite the fact that the COOH group is desolvated in the bidentate case. MM-FEP has the great advantage of being more easily generalized to other solvents and to be further improved which will be particularly useful to study solvent and (counter-)ion effects on interfacial reactions.

Modeling the HCOOH/CO2 Electrocatalytic Reaction: When Details Are Key.

Our first principles simulations of the electrooxidation of formic acid over nickel identify the reorientation of the formate intermediate and the desorption of CO2 as the rate-limiting steps. Although they are not associated with an electron transfer, these barriers are strongly modified when the electrochemical potential is explicitly accounted for and when modeling the influence of the solvent. Hence, such a level of modeling is key to understand the kinetic limitations that penalize the reaction.

Impacts of electrode potentials and solvents on the electroreduction of CO2: a comparison of theoretical approaches.

Since CO2 is a readily available feedstock throughout the world, the utilization of CO2 as a C1 building block for the synthesis of valuable chemicals is a highly attractive concept. However, due to its very nature of energy depleted "carbon sink", CO2 has a very low reactivity. Electrocatalysis offers the most attractive means to activate CO2 through reduction: the electron is the "cleanest" reducing agent whose energy can be tuned to the thermodynamic optimum. Under protic conditions, the reduction of CO2 over many metal electrodes results in formic acid. Thus, to open the road to its utilization as a C1 building block, the presence of water should be avoided to allow a more diverse chemistry, in particular for C-C bond formation with alkenes. Under those conditions, the intrinsic reactivity of CO2 can generate carbonates and oxalates by C-O and C-C bond formation, respectively. On Ni(111), almost exclusively carbonates and carbon monoxide are evidenced experimentally. Despite recent progress in modelling electrocatalytic reactions, determining the actual mechanism and selectivities between competing reaction pathways is still not straight forward. As a simple but important example of the intrinsic reactivity of CO2 under aprotic conditions, we highlight the shortcomings of the popular linear free energy relationship for electrode potentials (LFER-EP). Going beyond this zeroth order approximation by charging the surface and thus explicitly including the electrochemical potential into the electronic structure computations allows us to access more detailed insights, shedding light on coverage effects and on the influence of counterions.