How spin state and oxidation number of transition metal atoms determine molecular adsorption: a first-principles case study for NH3.

Understanding how the electronic state of transition metal atoms can influence molecular adsorption on a substrate is of great importance for many applications. Choosing NH3 as a model molecule, its adsorption behavior on defected SnS2 monolayers is investigated. The number of valence electrons n is controlled by decorating the monolayer with different transition metal atoms, ranging from Sc to Zn. Density-Functional Theory based calculations show that the adsorption energy of NH3 molecules oscillates with n and shows a clear odd-even pattern. There is also a mirror symmetry of the adsorption energies for large and low electron numbers. This unique behavior is mainly governed by the oxidation state of the TM ions. We trace back the observed trends of the adsorption energy to the orbital symmetries and ligand effects which affect the interaction between the 3σ orbitals (NH3) and the 3d orbitals of the transition metals. This result unravels the role which the spin state of TM ions plays in different crystal fields for the adsorption behavior of molecules. This new understanding of the role of the electronic structure on molecular adsorption can be useful for the design of high efficiency nanodevices in areas such as sensing and photocatalysis.

Amorphous MoS2 from a machine learning inter-atomic potential.

Amorphous molybdenum disulfide has shown potential as a hydrogen evolution catalyst, but the origin of its high activity is unclear, as is its atomic structure. Here, we have developed a classical inter-atomic potential using the charge equilibration neural network method, and we have employed it to generate atomic models of amorphous MoS2 by melting and quenching processes. The amorphous phase contains an abundance of molybdenum and sulfur atoms in low coordination. Besides the 6-coordinated molybdenum typical of the crystalline phases, a substantial fraction displays coordinations 4 and 5. The amorphous phase is also characterized by the appearance of direct S-S bonds. Density functional theory shows that the amorphous phase is metallic, with a considerable contribution of the 4-coordinated molybdenum to the density of states at the Fermi level. S-S bonds are related to the reduction of sulfur, with the excess electrons spread over several molybdenum atoms. Moreover, S-S bond formation is associated with a distinctive broadening of the 3s states, which could be exploited for experimental characterization of the amorphous phases. The large variety of local environments and the high density of electronic states at the Fermi level may play a positive role in increasing the electrocatalytic activity of this compound.

DFT insights into electrocatalytic CO2 reduction to methanol on α-Fe2O3(0001) surfaces.

Electrocatalytic reduction of CO2 to manufacture fuels and other useful chemicals is one of the appealing methods to reuse CO2. Herein, electrocatalytic CO2 reduction on a model α-Fe2O3(0001) surface catalyst has been investigated by means of density functional theory. This systematic study, involving 20 reaction intermediates and 63 distinct elementary reaction steps, has allowed the identification of a novel mechanism for the decomposition of the key intermediate *COOH. Methanol is the preferred product, with an overpotential of 0.8 V, over carbon monoxide (CO), formic acid (HCOOH), and formaldehyde (CH2O). Formaldehyde formed on the surface will be converted into methanol. This work demonstrates the need for a complete investigation of possible pathways to find the most favourable one, beyond chemical intuition. Moreover, it suggests that hematite could be an interesting material for CO2 reduction.

Catalytic properties of α-MnO2 for Li-air battery cathodes: a density functional investigation.

Details of the formation and dissociation of the first layer of Li2O2 on the α-MnO2(100) surface as the cathode in Li-air batteries have been studied using first principles density functional theory. The bias dependence of the electrochemical steps of charge (Li2O2 dissociation) and discharge (Li2O2 formation) via two different mechanisms has been studied. Discharge potential is found to be 2.94 V for the mechanism in which O2 adsorption is followed by lithiation. Charging potential for the reverse process is 3.37 V, giving an overpotential of 0.43 V, which is much lower than that on carbon electrodes. This is also in good agreement with experiments on α-MnO2 cathodes. In Li2O2 formation via the disproportionation of two LiO2 adsorbates, a maximum discharge potential of 2.61 V and a minimum charging potential of 3.48 V are obtained. The minimum energy pathway in this mechanism has a moderate kinetic barrier of 0.57 eV. Charging potentials of 3.37 V and 3.48 V imply that the typical charging potentials applied in the experiments (∼3.8 V) will dissociate the entire Li2O2 layer. These findings explain why α-MnO2 performs so well as a catalyst in Li-air battery cathodes, and suggest that a larger area of α-MnO2(100) can help reduce capacity loss.

Understanding the electrochemical double layer at the hematite/water interface: A first principles molecular dynamics study.

Using first principles molecular dynamics simulations, we probe the electrochemical double layer formed at the interface between the hematite surface and water. We consider two terminations of the (001) surface, viz., the fully hydroxylated (OH) and the stoichiometric (FeO3Fe) termination. We explicitly incorporate the counterions (Na+ and F-) in the solution, and model both specific and nonspecific adsorption of F- ions. We find that F- ions prefer to bind directly to the Fe ions (specific adsorption), with a substantial energy gain (0.75 eV/ion). We investigate the effect of the interface and the counterions on the dipole of individual water molecules. We find significant deviations of +0.2/-0.15 D for dipoles of the first solvation shell water molecules of F-/Na+ ions, respectively. Additionally, the hydration layers at the interface show an enhancement in the dipole moment resulting from stronger hydrogen bonding interactions between the water molecules and surface charged species. Furthermore, we analyze the electrostatic potential profile at the solid/liquid interface as a function of the kind of counterion present in the double layer and compute the capacitance of the compact (Helmholtz) layer. We find that our results (40.3 ± 3.5 µF/cm2 for the OH termination and 51 ± 5 µF/cm2 for the FeO3Fe termination) compare favorably with values reported by potentiometric titration based experimental studies (10-100 µF/cm2).

Magnetoelectric ϵ-Fe2O3: DFT study of a potential candidate for electrode material in photoelectrochemical cells.

The metastable iron oxide ϵ-Fe2O3 is rare but known for its magnetoelectric properties. While the more common alpha phase has been recognized for a long time as a suitable material for photoelectrochemical cells, its use is limited because of the electron-hole recombination problem when exposed to light. The indirect bandgap of the epsilon phase with its spontaneous polarization may offer a better potential for the application in photoelectrochemistry. Here, we report a detailed study of the electronic and structural features of the epsilon phase of iron oxide, its stability in thin films, and possible water dissociation reactions. Our studies are performed using density functional theory with a Hubbard-U correction. We observe that the stable ϵ-Fe2O3 surfaces favor the dissociation of water. The average difference in the energies of the states when water is adsorbed and when it is dissociated is roughly found to be -0.40 eV. Our results compare with the available experimental results where the epsilon phase is reported to be more efficient for the release of hydrogen from renewable oxygenates when exposed to sunlight.

Passivation of surface states of α-Fe2O3(0001) surface by deposition of Ga2O3 overlayers: A density functional theory study.

There is a big debate in the community regarding the role of surface states of hematite in the photoelectrochemical water splitting. Experimental studies on non-catalytic overlayers passivating the hematite surface states claim a favorable reduction in the overpotential for the water splitting reaction. As a first step towards understanding the effect of these overlayers, we have studied the system Ga2O3 overlayers on hematite (0001) surfaces using first principles computations in the PBE+U framework. Our computations suggest that stoichiometric terminations of Ga2O3 overlayers are energetically more favored than the bare surface, at ambient oxygen chemical potentials. Energetics suggest that the overlayers prefer to grow via a layer-plus-island (Stranski-Krastanov) growth mode with a critical layer thickness of 1-2 layers. Thus, a complete wetting of the hematite surface by an overlayer of gallium oxide is thermodynamically favored. We establish that the effect of deposition of the Ga2O3 overlayers on the bare hematite surface is to passivate the surface states for the stoichiometric termination. For the oxygen terminated surface which is the most stable termination under photoelectrochemical conditions, the effect of deposition of the Ga2O3 overlayer is to passivate the hole-trapping surface state.

CO2 conversion to methanol on Cu(I) oxide nanolayers and clusters: an electronic structure insight into the reaction mechanism.

The mechanism of carbon dioxide reduction to methanol on Cu(I) oxide nanolayers and clusters using water as the source of hydrogen was traced using density functional theory. The nature of the active sites is revealed, namely the role of surface copper dimers, which are present on the Cu2O(001) surface and in the nanoclusters of size Cu32O16 and Cu14O7. The major difference between metal catalysts and Cu2O is outlined: the CO2 molecule interacts strongly with the oxide and undergoes bending prior to hydrogenation. The first step of CO2 hydrogenation results in the formation of a stable carboxyl intermediate, -CO(OH), which in the following steps is converted to methanol via formic acid and formaldehyde intermediates. The consumption of hydrogen from water leaves surface peroxo- and hydroperoxo-species. The peroxides easily desorb molecular oxygen, while for hydroperoxides the reaction of oxygen evolution requires an activation energy of 130 kJ mol(-1). The maxima in the absorption spectra correspond well with the required activation energies in the elementary steps.

Nitrogen electrochemically reduced to ammonia with hematite: density-functional insights.

Licht et al. (Science, 2014, 345, 637) recently proposed a procedure to synthesize NH3 from N2 and by steam electrolysis in molten hydroxide suspensions of nano-Fe2O3. This highly exciting investigation undoubtedly boosts the hope of the CO2-free and low-cost ammonia industry. To provide insights at the atomistic level into the reduction process of N2, we have carried out a density-functional study on the electrochemical formation of NH3 molecules on hematite(0001) surfaces. By considering associative and dissociative mechanisms, we have identified a reaction path that requires an applied bias of -1.1 V to allow the proton transfer processes to occur downhill. The most energy-demanding step is the addition of the first proton to the adsorbed molecular nitrogen. The computed bias is in good agreement with experimental electrolysis potentials that activate the electric current.

Defective α-Fe2O3(0001): an ab initio study.

By using density functional theory calculations at the PBE+U level, we investigated the properties of hematite (0001) surfaces decorated with adatoms/vacancies/substituents. For the most stable surface termination over a large range of oxygen chemical potentials (muO), the vacancy formation and adsorption energies were determined as a function of muO. Under oxygen-rich conditions, all defects are metastable with respect to the ideal surface. Under oxygen-poor conditions, O vacancies and Fe adatoms become stable. Under ambient conditions, all defects are metastable; in the bulk, O vacancies form more easily than Fe vacancies, whereas at the surface the opposite is true. All defects, that is, O and Fe vacancies, Fe and Al adatoms, and Al substituents, induce important modifications to the geometry of the surface in their vicinity. Dissociative adsorption of molecular oxygen is likely to be exothermic on surfaces with Fe/Al adatoms or O vacancies.

Photo-driven oxidation of water on α-Fe2O3 surfaces: an ab initio study.

Adopting the theoretical scheme developed by the Nørskov group [see, for example, Nørskov et al., J. Phys. Chem. B 108, 17886 (2004)], we conducted a density functional theory study of photo-driven oxidation processes of water on various terminations of the clean hematite (α-Fe2O3) (0001) surface, explicitly taking into account the strong correlation among the 3d states of iron through the Hubbard U parameter. Six best-known terminations, namely, Fe−Fe−O3− (we call S1), O−Fe−Fe−(S2), O2−Fe−Fe−(S3), O3−Fe−Fe− (S4), Fe−O3−Fe− (S5), and O−Fe−O3−(S6), are first exposed to water, the stability of resulting surfaces is investigated under photoelectrochemical conditions by considering different chemical reactions (and their reaction free energies) that lead to surfaces covered by O atoms or/and OH groups. Assuming that the water splitting reaction is driven by the redox potential for photogenerated holes with respect to the normal hydrogen electrode, UVB, at voltage larger than UVB, most 3-oxygen terminated substrates are stable. These results thus suggest that the surface, hydroxylated in the dark, should release protons under illumination. Considering the surface free energy of all the possible terminations shows that O3­S5 and O3­S1 are the most thermodynamically stable. While water oxidation process on the former requires an overpotential of 1.22 V, only 0.84 V is needed on the latter.

The Role of Water Molecules on Polaron Behavior at Rutile (110) Surface: A Constrained Density Functional Theory Study.

Understanding the behavior of a polaron in contact with water is of significant importance for many photocatalytic applications. We investigated the influence of water on the localization and transport properties of polarons at the rutile (110) surface by constrained density functional theory. An excess electron at a dry surface favors the formation of a small polaron at the subsurface Ti site, with a preferred transport direction along the [001] axis. As the surface is covered by water, the preferred spatial localization of the polarons is moved from the subsurface to the surface. When the water coverage exceeds half a monolayer, the preferred direction of polaron hopping is changed to the [110] direction toward the surface. This characteristic behavior is related to the Ti3d-orbital occupations and crystal field splitting induced by different distorted structures under water coverage. Our work describes the reduced sites that might eventually play a role in photocatalysis for rutile (110) surfaces in a water environment.

Water adsorption and dissociation on α-Fe2O3(0001): PBE+U calculations.

Adsorption and dissociation of water on different oxygen- and iron-terminated hematite(0001) surfaces at monolayer coverage have been studied by density-functional theory calculations, including a Hubbard-like+U correction. We considered six possible surface terminations, including four oxygen- and two iron-terminations. Binding energy of water on these terminations can be as large as 1.0 eV. On these terminations the energy barrier for the dissociation of the molecularly adsorbed water is less than 0.3 eV, and in few cases the dissociation is even spontaneous, i.e., without any detectable barrier. Our results thus suggest that water can be adsorbed on the α-Fe2O3(0001) surface dissociatively at room temperature, as previously found by experiment. This study also presents a very first theoretical insight into the adsorption and dissociation of water on all known terminations of the hematite(0001) surface.

High-Throughput Screening of Promising Redox-Active Molecules with MolGAT.

Redox flow batteries (RFBs) have emerged as a promising option for large-scale energy storage, owing to their high energy density, low cost, and environmental benefits. However, the identification of organic compounds with high redox activity, aqueous solubility, stability, and fast redox kinetics is a crucial and challenging step in developing an RFB technology. Density functional theory-based computational materials prediction and screening is a time-consuming and computationally expensive technique, yet it has a high success rate. To speed up the discovery of new materials with desired properties, machine-learning-based models can be trained on large data sets. Graph neural networks (GNNs) are particularly well-suited for non-Euclidean data and can model complex relationships, making them ideal for accelerating the discovery of novel materials. In this study, a GNN-based model called MolGAT was developed to predict the redox potential of organic molecules using molecular structures, atomic properties, and bond attributes. The model was trained on a data set of over 15,000 compounds with redox potentials ranging from -4.11 to 2.56. MolGAT outperformed other GNN variants, such as the Graph Attention Network, Graph Convolution Network, and AttentiveFP models. The trained model was used to screen a vast chemical data set comprising 581,014 molecules, namely OMDB, QM9, ZINC, CHEMBL, and DELANEY, and identified 23,467 potential redox-active compounds for use in redox flow batteries. Of those, 20,716 molecules were identified as potential catholytes with predicted redox potentials up to 2.87 V, while 2,751 molecules were deemed potential anolytes with predicted redox potentials as low as -2.88 V. This work demonstrates the capabilities of graph neural networks in condensed matter physics and materials science to screen promising redox-active species for further electronic structure calculations and experimental testing.

The carbonyl-lock mechanism underlying non-aromatic fluorescence in biological matter.

Challenging the basis of our chemical intuition, recent experimental evidence reveals the presence of a new type of intrinsic fluorescence in biomolecules that exists even in the absence of aromatic or electronically conjugated chemical compounds. The origin of this phenomenon has remained elusive so far. In the present study, we identify a mechanism underlying this new type of fluorescence in different biological aggregates. By employing non-adiabatic ab initio molecular dynamics simulations combined with a data-driven approach, we characterize the typical ultrafast non-radiative relaxation pathways active in non-fluorescent peptides. We show that the key vibrational mode for the non-radiative decay towards the ground state is the carbonyl elongation. Non-aromatic fluorescence appears to emerge from blocking this mode with strong local interactions such as hydrogen bonds. While we cannot rule out the existence of alternative non-aromatic fluorescence mechanisms in other systems, we demonstrate that this carbonyl-lock mechanism for trapping the excited state leads to the fluorescence yield increase observed experimentally, and set the stage for design principles to realize novel non-invasive biocompatible probes with applications in bioimaging, sensing, and biophotonics.

Ab initio parameterization of an all-atom polarizable and dissociable force field for water.

A novel all-atom, dissociative, and polarizable force field for water is presented. The force field is parameterized based on forces, stresses, and energies obtained form ab initio calculations of liquid water at ambient conditions. The accuracy of the force field is tested by calculating structural and dynamical properties of liquid water and the energetics of small water clusters. The transferability of the force field to dissociated states is studied by considering the solvation of a proton and the ionization of water at extreme conditions of pressure and temperature. In the case of the solvated proton, the force field properly describes the presence of both Eigen and Zundel configurations. In the case of the pressure-induced ice VIII/ice X transition and the temperature-induced transition to a superionic phase, the force field is found to describe accurately the proton symmetrization and the melting of the proton sublattice, respectively.

Carbon in palladium catalysts: A metastable carbide.

The catalytic activity of palladium toward selective hydrogenation of hydrocarbons depends on the partial pressure of hydrogen. It has been suggested that the reaction proceeds selectively toward partial hydrogenation only when a carbon-rich film is present at the metal surface. On the basis of first-principles simulations, we show that carbon can dissolve into the metal because graphite formation is delayed by the large critical nucleus necessary for graphite nucleation. A bulk carbide Pd(6)C with a hexagonal six-layer fcc-like supercell forms. The structure is characterized by core level shifts of 0.66-0.70 eV in the core states of Pd, in agreement with experimental x-ray photoemission spectra. Moreover, this phase traps bulk-dissolved hydrogen, suppressing the total hydrogenation reaction channel and fostering partial hydrogenation.

Heterostructures of ε-Fe2O3 and α-Fe2O3: insights from density functional theory.

Many materials used in energy devices or applications suffer from the problem of electron-hole pair recombination. One promising way to overcome this problem is the use of heterostructures in place of a single material. If an electric dipole forms at the interface, such a structure can lead to a more efficient electron-hole pair separation and thus prevent recombination. Here we model and study a heterostructure comprised of two polymorphs of Fe2O3. Each one of the two polymorphs, α-Fe2O3 and ε-Fe2O3, individually shows promise for applications in photoelectrochemical cells. The heterostructure of these two materials is modeled by means of density functional theory. We consider both ferromagnetic as well as anti-ferromagnetic couplings at the interface between the two systems. Both individual oxides are insulating in nature and have an anti-ferromagnetic spin arrangement in their ground state. The same properties are found also in their heterostructure. The highest occupied electronic orbitals of the combined system are localized at the interface between the two iron-oxides. The localization of charges at the interface is characterized by electrons residing close to the oxygen atoms of ε-Fe2O3 and electron-holes localized on the iron atoms of α-Fe2O3, just around the interface. The band alignment at the interface of the two oxides shows a type-III broken band-gap heterostructure. The band edges of α-Fe2O3 are higher in energy than those of ε-Fe2O3. This band alignment favours a spontaneous transfer of excited photo-electrons from the conduction band of α- to the conduction band of ε-Fe2O3. Similarly, photo-generated holes are transferred from the valence band of ε- to the valence band of α-Fe2O3. Thus, the interface favours a spontaneous separation of electrons and holes in space. The conduction band of ε-Fe2O3, lying close to the valence band of α-Fe2O3, can result in band-to-band tunneling of electrons which is a characteristic property of such type-III broken band-gap heterostructures and has potential applications in tunnel field-effect transistors.

Tuning optoelectronic properties of triphenylamine based dyes through variation of pi-conjugated units and anchoring groups: A DFT/TD-DFT investigation.

Dye-sensitized solar cells (DSSCs) have attracted widespread attention due to their unique features. In the present work, molecular engineered triphenylamine based dyes featuring donor-bridge-acceptor architecture have been considered and investigated for suitable properties for DSSCs applications. Hydantoin anchoring group has been introduced replacing the commonly used cyanoacrylic acid to improve the long-term stability of the device. Results on the effects of varied anchoring groups and pi-spacers have been interpreted from the viewpoint of DFT/TD-DFT calculations. Designed sensitizers exhibit suitable light-harvesting efficiencies, excited-state lifetimes, electron injection and regeneration abilities. Red-shifted electronic spectra are observed for three hydantoin dyes compared to others in the same family. Further analysis of chemical descriptors and observation from full-electron donor-acceptor map reveal that the three dyes among nine are potential materials with promising properties towards improving DSSCs performance.

Ab initio thermodynamics of lithium oxides: from bulk phases to nanoparticles.

Lithium ion batteries are nowadays key devices for energy storage, and a great research effort is under way to develop and apply new materials. Recently, new approaches have been proposed that rely on the reversible formation of either Li2O or Li2O2 at the electrodes. The details of their formation and dissolution are, however, still unclear. As a first step towards the understanding of these processes, bulk lithium oxides, their surfaces and their nanoparticles have been here investigated by density functional theory and ab initio thermodynamics. At a pressure of 1 atmosphere of oxygen, Li2O2 is the stable bulk phase below 5 K, where a transition to Li2O takes place. Wulff's construction predicts an octahedral shape for Li2O nanoparticles and the form of a hexagonal prism for Li2O2. By taking into account the effect of the surfaces, a size-dependent phase diagram is calculated. At an oxygen pressure of 1 atmosphere and a temperature of 300 K, Li2O is the stable phase for particles with a diameter larger than approximately 2.5 nm. This size-dependent oxidation behavior of lithium should be taken into account in the design of nanostructured oxygen cathodes.