Improved description of hematite surfaces by the SCAN functional

Density Functional Theory Study of the Spin-Orbit Insulating Phase in SnTe Cubic Nanowires: Implications for Topological Electronics

We investigate the electronic, structural, and topological properties of the SnTe and PbTe cubic nanowires using ab initio calculations. Using standard and linear-scale density functional theory, we go from the ultrathin limit up to the nanowire thicknesses observed experimentally. Finite-size effects in the ultrathin limit produce an electric quadrupole and associated structural distortions; these distortions increase the band gap, but they get reduced with the size of the nanowires and become less and less relevant. Ultrathin SnTe cubic nanowires are trivial band gap insulators; we demonstrate that by increasing the thickness, there is an electronic transition to a spin-orbit insulating phase due to trivial surface states in the regime of thin nanowires. These trivial surface states with a spin-orbit gap of a few meV appear at the same k-point of the topological surface states. Going to the limit of thick nanowires, we should observe the transition to the topological crystalline insulator phase with the presence of two massive surface Dirac fermions hybridized with the persistent trivial surface states. Therefore, we have the copresence of massive Dirac surface states and trivial surface states close to the Fermi level in the same region of the k-space. According to our estimation, the cubic SnTe nanowires are trivial insulators below the critical thickness tc1 = 10 nm, and they become spin-orbit insulators between tc1 = 10 nm and tc2 = 17 nm, while they transit to the topological phase above the critical thickness of tc2 = 17 nm. These critical thickness values are in the range of typical experimental thicknesses, making the thickness a relevant parameter for the synthesis of topological cubic nanowires. Pb1-xSnxTe nanowires would have both these critical thicknesses tc1 and tc2 at larger values depending on the doping concentration. We discuss the limitations of density functional theory in the context of topological nanowires and the consequences of our results on topological electronics.

Structural, electronic and magnetic properties of greigite Fe3S4by GGA and GGA+Uversus SCAN meta-GGA density functionals

The performance of exchange-correlation functional of density functional theory represented in generalized gradient approximation (GGA) and in the strongly constrained and appropriately normed (SCAN) meta-GGA scheme to study structural, electronic, and magnetic properties of greigite (Fe3S4) was investigated. The effects of inclusion of strong electron correlations represented by on-site Hubbard correctionU, and nonlocality of the long-range van der Waals (vdW) interactions were also considered. Geometry optimization yielded the inverse spinel structure and lattice parameter of greigite in good agreement with experimental data. Calculated electronic structure revealed a half-metallic nature of the greigite bands for the applied functionals except for GGA, which predicts metallic behavior. Antiferromagnetic coupling of iron ions in tetrahedral and octahedral coordinations makes the overall crystal structure ferrimagnetic. In general the GGA+Uand SCAN show comparable performance in prediction physical properties of greigite. Inclusion of the vdW correction does not change the character of the bands.

Calculating shear viscosity with confined non-equilibrium molecular dynamics: a case study on hematite - PAO-2 lubricant

The behaviour of confined lubricants at the atomic scale as affected by the interactions at the surface-lubricant interface is relevant in a range of technological applications in areas such as the automotive industry. In this paper, by performing fully atomistic molecular dynamics, we investigate the regime where the viscosity starts to deviate from the bulk behaviour, a topic of great practical and scientific relevance. The simulations consist of setting up a shear flow by confining the lubricant between iron oxide surfaces. By using confined Non-Equilibrium Molecular Dynamics (NEMD) simulations at a pressure range of 0.1-1.0 GPa at 100 °C, we demonstrate that the film thickness of the fluid affects the behaviour of viscosity. We find that by increasing the number of lubricant molecules, we approach the viscosity value of the bulk fluid derived from previously published NEMD simulations for the same system. These changes in viscosity occurred at film thicknesses ranging from 10.12 to 55.93 Å. The viscosity deviations at different pressures between the system with the greatest number of lubricant molecules and the bulk simulations varied from -16% to 41%. The choice of the utilized force field for treating the atomic interactions was also investigated.

Excitation-Wavelength-Dependent Charge-Carrier Lifetime in Hematite: An Insight from Nonadiabatic Molecular Dynamics

Experiments have reported that the photoexcited carrier lifetime in α-Fe₂O₃ has a significant excitation-wavelength dependence but leave the physical mechanism unresolved. In this work, we rationalize the puzzling excitation-wavelength dependence of the photoexcited carrier dynamics in Fe₂O₃ by performing nonadiabatic molecular dynamics simulation based on the strongly constrained and appropriately normed functional, which accurately describes the electronic structure of Fe₂O₃. Photogenerated electrons with lower-energy excitation relax fast in the t_2g conduction band within about 100 fs, while the photogenerated electrons with higher-energy excitation undergo first a slower interband relaxation from the e_g lower state to the t_2g upper state on a time scale of 135 ps, followed by the much faster t_2g intraband relaxation. This study provides insight into the experimentally reported excitation-wavelength dependence of the carrier lifetime in Fe₂O₃ and a reference for regulating photogenerated carrier dynamics in transition-metal oxides through the light excitation wavelength.

Passivation of Hematite by a Semiconducting Overlayer Reduces Charge Recombination: An Insight from Nonadiabatic Molecular Dynamics

Hematite (α-Fe₂O₃) is a promising photoanode material for photoelectrochemical water splitting. Surface-passivating layers are effective in improving water oxidation kinetics; however, the passivation mechanism is not fully understood due to the complexity of interfacial reactions. Focusing on the Fe-terminated Fe₂O₃ (0001) surface that exhibits surface states in the band gap, we perform ab initio quantum dynamics simulations to study the effect of an α-Ga₂O₃ overlayer on charge recombination. The overlayer eliminates surface states and suppresses charge recombination 4-fold. This explains in part the observed cathodic shift in the onset potential for water oxidation. The increased charge carrier lifetime is an outcome of two factors, energy gap and electron-vibrational coupling, with a positive contribution from the former but a negative contribution from the latter. This work presents an advance in the atomistic time-domain understanding of the influence of surface passivation on charge recombination dynamics and provides guidance for designing novel α-Fe₂O₃ photoanodes.

Adsorption of small organic acids and polyphenols on hematite surfaces: Density Functional Theory + thermodynamics analysis

HYPOTHESIS

The interactions of organic molecules with mineral surfaces are influenced by several factors such as adsorbate speciation, surface atomic and electronic structure, and environmental conditions. When coupled with thermodynamic techniques, energetics from atomistic modeling can provide a molecular-level picture of which factors determine reactivity. This is paramount for evaluating the chemical processes which control the fate of these species in the environment.

EXPERIMENTS

Inner-sphere adsorption of oxalate and pyrocatechol on (001), (110), and (012) α-Fe₂O₃ surfaces was modeled using Density Functional Theory (DFT). Unique bidentate binding modes were sampled along each facet to study how different adsorbate and surface factors govern site preference. Adsorption energetics were then calculated using a DFT + thermodynamics approach which combines DFT energies with tabulated data and Nernst-based corrective terms to incorporate different experimental parameters.

FINDINGS

Instead of a universal trend, each facet displays a unique factor that dominates site preference based on either strain (001), functional groups (110), or topography (012). Adsorption energies predict favorable inner-sphere adsorption for both molecules but opposite energetic trends with varying pH. Additionally, vibrational analysis was conducted for each system and compared to experimental IR data. The work presented here provides an effective, computational methodology to study numerous adsorption processes occurring at the surface-aqueous interface.

Intrinsic defect formation and the effect of transition metal doping on transport properties in a ductile thermoelectric material α-Ag2S: a first-principles study

In this paper, the electronic structure and transport properties of a ductile thermoelectric material α-Ag2S are examined using first-principles calculations combined with the Boltzmann transport equation within a constant relaxation-time approximation. The use of the exchange-correlation functional SCAN + rVV10 successfully describes the geometric and electronic structure of α-Ag2S with a direct bandgap value of 0.99 eV, which is consistent with the previous experimental observations. Based on the calculations of the formation energy of typical intrinsic defects, it is found that intrinsic defect formation greatly affects the conductivity of the system where silver vacancy and interstitial silver act as p-type and n-type defects, respectively. Large Seebeck coefficients at room-temperature, of around -760 μV K-1 for n-type and 1400 μV K-1 for p-type, are realized. It is also suggested that the doping of fully filled d-block elements such as Cu and Au not only maintained the Seebeck coefficients at high values but also improved electrical conductivity by more than 1.4 times, leading to the improvement of the power factor by up to 40% compared to the non-doped sample at low carrier concentration.

Polaronic structure of excess electrons and holes for a series of bulk iron oxides

With the use of a gap-optimized hybrid functional and large supercells, it is found that while the electron hole polaron generally localises onto a single iron site, the electron polaron localises across two iron sites of the same spin layer.