An efficient way to model complex magnetite: Assessment of SCC-DFTB against DFT

Computational Exploration of Phenolic Compounds in Corrosion Inhibition: A Case Study of Hydroxytyrosol and Tyrosol

The corrosion of materials remains a critical challenge with significant economic and infrastructural impacts. A comprehensive understanding of adsorption characteristics of phytochemicals can facilitate the effective design of high-performing environmentally friendly inhibitors. This study conducted a computational exploration of hydroxytyrosol (HTR) and tyrosol (TRS) (potent phenolic compounds found in olive leaf extracts), focusing on their adsorption and reactivity on iron surfaces. Utilizing self-consistent-charge density-functional tight-binding (SCC-DFTB) simulations, molecular dynamics (MD) simulations, and quantum chemical calculations (QCCs), we investigated the molecules' structural and electronic attributes and interactions with iron surfaces. The SCC-DFTB results highlighted that HTR and TRS coordinated with iron atoms when adsorbed individually, but only HTR maintained bonding when adsorbed alongside TRS. At their individual adsorption, HTR and TRS had interaction energies of -1.874 and -1.598 eV, which became more negative when put together (-1.976 eV). The MD simulations revealed parallel adsorption under aqueous and vacuum conditions, with HTR demonstrating higher adsorption energy. The analysis of quantum chemical parameters, including global and local reactivity descriptors, offered crucial insights into molecular reactivity, stability, and interaction-prone atomic sites. QCCs revealed that the fraction of transferred electron ∆N aligned with SCC-DFTB results, while other parameters of purely isolated molecules failed to predict the same. These findings pave the way for potential advancements in anticorrosion strategies leveraging phenolic compounds.

Unraveling Bonding Mechanisms and Electronic Structure of Pyridine Oximes on Fe(110) Surface: Deeper Insights from DFT, Molecular Dynamics and SCC-DFT Tight Binding Simulations

The development of corrosion inhibitors with outstanding performance is a never-ending and complex process engaged in by researchers, engineers and practitioners. The computational assessment of organic corrosion inhibitors' performance is a crucial step towards the design of new task-specific materials. Herein, the electronic features, adsorption characteristics and bonding mechanisms of two pyridine oximes, namely 2-pyridylaldoxime (2POH) and 3-pyridylaldoxime (3POH), with the iron surface were investigated using molecular dynamics (MD), and self-consistent-charge density-functional tight-binding (SCC-DFTB) simulations. SCC-DFTB simulations revealed that the 3POH molecule can form covalent bonds with iron atoms in its neutral and protonated states, while the 2POH molecule can only bond with iron through its protonated form, resulting in interaction energies of -2.534, -2.007, -1.897, and -0.007 eV for 3POH, 3POH⁺, 2POH⁺, and 2POH, respectively. Projected density of states (PDOSs) analysis of pyridines-Fe(110) interactions indicated that pyridine molecules were chemically adsorbed on the iron surface. Quantum chemical calculations (QCCs) revealed that the energy gap and Hard and Soft Acids and Bases (HSAB) principles were efficient in predicting the bonding trend of the molecules investigated with an iron surface. 3POH had the lowest energy gap of 1.706 eV, followed by 3POH⁺ (2.806 eV), 2POH⁺ (3.121 eV), and 2POH (3.431 eV). In the presence of a simulated solution, MD simulation showed that the neutral and protonated forms of molecules exhibited a parallel adsorption mode on an iron surface. The excellent adsorption properties and corrosion inhibition performance of 3POH may be attributed to its low stability compared to 2POH molecules.

Improving the Oxygen Evolution Reaction on Fe3O4(001) with Single-Atom Catalysts

Doping magnetite surfaces with transition-metal atoms is a promising strategy to improve the catalytic performance toward the oxygen evolution reaction (OER), which governs the overall efficiency of water electrolysis and hydrogen production. In this work, we investigated the Fe₃O₄(001) surface as a support material for single-atom catalysts of the OER. First, we prepared and optimized models of inexpensive and abundant transition-metal atoms, such as Ti, Co, Ni, and Cu, trapped in various configurations on the Fe₃O₄(001) surface. Then, we studied their structural, electronic, and magnetic properties through HSE06 hybrid functional calculations. As a further step, we investigated the performance of these model electrocatalysts toward the OER, considering different possible mechanisms, in comparison with the pristine magnetite surface, on the basis of the computational hydrogen electrode model developed by Nørskov and co-workers. Cobalt-doped systems were found to be the most promising electrocatalytic systems among those considered in this work. Overpotential values (∼0.35 V) were in the range of those experimentally reported for mixed Co/Fe oxide (0.2-0.5 V).

An Efficient Way to Model Complex Iron Carbides: A Benchmark Study of DFTB2 against DFT

Iron carbides have attracted increasing attention in recent years due to their enormous potential in catalytic fields, such as Fischer-Tropsch synthesis and the growth of carbon nanotubes. Theoretical calculations can provide a more thorough understanding of these reactions at the atomic scale. However, due to the extreme complexity of the active phases and surface structures of iron carbides at the operando conditions, calculations based on density functional theory (DFT) are too costly for realistically large models of iron carbide particles. Therefore, a cheap and efficient quantum mechanical simulation method with accuracy comparable to DFT is desired. In this work, we adopt the spin-polarized self-consistent charge density functional tight-binding (DFTB2) method for iron carbides by reparametrization of the repulsive part of the Fe-C interactions. To assess the performance of the improved parameters, the structural and electronic properties of iron carbide bulks and clusters obtained with DFTB2 method are compared with the previous experimental values and the results obtained with DFT approach. Calculated lattice parameters and density of states are close to DFT predictions. The benchmark results show that the proposed parametrization of Fe-C interactions provides transferable and balanced description of iron carbide systems. Therefore, spin-polarized DFTB2 is valued as an efficient and reliable method for the description of iron carbide systems.

Multiscale simulations of the hydration shells surrounding spherical Fe3O4 nanoparticles and effect on magnetic properties

Iron oxide magnetic nanoparticles (NPs) are excellent systems in catalysis and in nanomedicine, where they are mostly immersed in aqueous media. Even though the NP solvation by water is expected to play an active role, the detailed structural insight at the nanostructure oxide/water interface is still missing. Here, based on our previous efforts to obtain accurate models of dehydrated Fe3O4 NPs and of their magnetic properties and through multiscale molecular dynamics simulations combining the density functional tight binding method and force field, we unravel the atomistic details of the short range (chemical) and long range (physical) interfacial effects when magnetite nanoparticles are immersed in water. The influence of the first hydration shell on the structural, electronic and magnetic properties of Fe3O4 NPs is revealed by high-level hybrid density functional calculations. Hydrated Fe3O4 NPs possess larger magnetic moment than dehydrated ones. This work bridges the large gap between experimental studies on solvated Fe3O4 NPs and theoretical investigations on flat Fe3O4 surfaces covered with water and paves the way for further study of Fe3O4 NPs in biological environments.

Parametrization of the Fe-Owater cross-interaction for a more accurate Fe3O4/water interface model and its application to a spherical Fe3O4 nanoparticle of realistic size

The accurate description of iron oxides/water interfaces requires reliable force field parameters that can be developed through comparison with sophisticated quantum mechanical calculations. Here, a set of CLASS2 force field parameters is optimized to describe the Fe-O_water cross-interaction through comparison with hybrid density functional theory (HSE06) calculations of the potential energy function for a single water molecule adsorbed on the Fe₃O₄ (001) surface and with density functional tight binding (DFTB+U) molecular dynamics simulations for a water trilayer on the same surface. The performance of the new parameters is assessed through the analysis of the number density profile of a water bulk (12 nm) sandwiched between two magnetite slabs of large surface area. Their transferability is tested for water adsorption on the curved surface of a spherical Fe₃O₄ nanoparticle of realistic size (2.5 nm).

DFTB+, a software package for efficient approximate density functional theory based atomistic simulations

DFTB+ is a versatile community developed open source software package offering fast and efficient methods for carrying out atomistic quantum mechanical simulations. By implementing various methods approximating density functional theory (DFT), such as the density functional based tight binding (DFTB) and the extended tight binding method, it enables simulations of large systems and long timescales with reasonable accuracy while being considerably faster for typical simulations than the respective ab initio methods. Based on the DFTB framework, it additionally offers approximated versions of various DFT extensions including hybrid functionals, time dependent formalism for treating excited systems, electron transport using non-equilibrium Green's functions, and many more. DFTB+ can be used as a user-friendly standalone application in addition to being embedded into other software packages as a library or acting as a calculation-server accessed by socket communication. We give an overview of the recently developed capabilities of the DFTB+ code, demonstrating with a few use case examples, discuss the strengths and weaknesses of the various features, and also discuss on-going developments and possible future perspectives.

Insight into the interface between Fe3O4 (001) surface and water overlayers through multiscale molecular dynamics simulations

In this work, we investigate the Fe₃O₄ (001) surface/water interface by combining several theoretical approaches, ranging from a hybrid functional method (HSE06) to density-functional tight-binding (DFTB) to molecular mechanics (MM). First, we assess the accuracy of the DFTB method to correctly reproduce HSE06 results on structural details and energetics and available experimental data for adsorption of isolated water, dimers, and trimers up to a water monolayer. Second, we build two possible configurations of a second and a third overlayer and perform molecular dynamics simulations with DFTB, monitoring the water orientation, the H-bond network, and the ordered water structure formation. To make our models more realistic, we then build a 12 nm-thick water multilayer on top of the Fe₃O₄ (001) surface slab model, which we investigate through MM-molecular dynamics (MD). The water layer structuring, revealed by the analysis of the atomic positions from a long MM-MD run for this large MM model, extends up to about 6-7 Å and nicely compares with that observed for a water trilayer model. However, MM and DFTB MD simulations show some discrepancy due to the poor description of the Fe⋯OH₂ distance in MM that calls for further work in the parameterization of the model.

Photoelectrochemistry of Ferrites: Theoretical Predictions vs. Experimental Results

This paper gives an overview about recent theoretical and experimental work on electronic and optical properties of spinel ferrites MFe2O4. These compounds have come into focus of research due to their possible application as photocatalyst material for photoelectrochemical water splitting. The theoretical background of state-of-the-art quantum-chemical approaches applied for predicting electronic and optical band gaps, absolute band positions, optical absorption spectra, dielectric functions and Raman spectra, is briefly reviewed. Recent applications of first-principles methods on magnetic and electronic properties of ferrites with M = Mg and the first row of subgroup elements Sc to Zn are presented, where it is shown that the fundamental band gap is strongly dependent on the spin state and the degree of inversion of the spinel structure. The observed variation of electronic properties may serve as an explanation for the large scattering of experimental results. The exchange of M and Fe cations has also a pronounced effect on the Raman spectra of ferrites, which is analyzed at atomic scale from first principles. Calculated optical absorption spectra of ferrites are compared to experimental spectra. The electronic nature of the first excitations and the role of oxygen vacancies are discussed. For the calculation of absolute band positions, which have a significant impact on the photoelectrochemical activity of the ferrites, models of the most stable ferrite surfaces are developed that take into account their polar nature and the interaction with the solvent. Theoretically predicted valence and conduction band edges are compared to results from electrochemical measurements. The role of cation exchange on the surface electronic structure is investigated both theoretically and experimentally.

Shaping Magnetite Nanoparticles from First Principles

Iron oxide magnetic nanoparticles (NPs) are stimuli-responsive materials at the forefront of nanomedicine. Their realistic finite temperature simulations are a formidable challenge for first-principles methods. Here, we use density functional tight binding to open up the required time and length scales and obtain global minimum structures of Fe_{3}O_{4} NPs of realistic size (1400 atoms, 2.5 nm) and of different shapes, which we then refine with hybrid density functional theory methods to accomplish proper electronic and magnetic properties, which have never been accurately described in simulations. On this basis, we develop a general empirical formula and prove its predictive power for the evaluation of the total magnetic moment of Fe_{3}O_{4} NPs. By converting the total magnetic moment into the macroscopic saturation magnetization, we rationalize the experimentally observed dependence with shape. We also reveal interesting reconstruction mechanisms and unexpected patterns of charge ordering.