Molecular-scale mechanisms of CO2 mineralization in nanoscale interfacial water films.

The calamitous impacts of unabated carbon emission from fossil-fuel-burning energy infrastructure call for accelerated development of large-scale CO2 capture, utilization and storage technologies that are underpinned by a fundamental understanding of the chemical processes at a molecular level. In the subsurface, rocks rich in divalent metals can react with CO2, permanently sequestering it in the form of stable metal carbonate minerals, with the CO2-H2O composition of the post-injection pore fluid acting as a primary control variable. In this Review, we discuss mechanistic reaction pathways for aqueous-mediated carbonation with carbon mineralization occurring in nanoscale adsorbed water films. In the extreme of pores filled with a CO2-dominant fluid, carbonation reactions are confined to angstrom to nanometre-thick water films coating mineral surfaces, which enable metal cation release, transport, nucleation and crystallization of metal carbonate minerals. Although seemingly counterintuitive, laboratory studies have demonstrated facile carbonation rates in these low-water environments, for which a better mechanistic understanding has come to light in recent years. The overarching objective of this Review is to delineate the unique underlying molecular-scale reaction mechanisms that govern CO2 mineralization in these reactive and dynamic quasi-2D interfaces. We highlight the importance of understanding unique properties in thin water films, such as how water dielectric properties, and consequently ion solvation and hydration behaviour, can change under nanoconfinement. We conclude by identifying important frontiers for future work and opportunities to exploit these fundamental chemical insights for decarbonization technologies in the twenty-first century.

Connecting particle interactions to agglomerate morphology and rheology of boehmite nanocrystal suspensions.

HYPOTHESIS: The rheology of complex suspensions, such as nuclear waste slurries at the Hanford and Savannah River sites, imposes significant challenges on industrial-scale processing. Investigating the rheology and connecting it to the agglomerate morphology and underlying particle interactions in slurries will provide important fundamental knowledge, as well as prescriptive data for practical applications. Here, we use suspensions of nano-scale aluminum oxyhydroxide minerals in the form of boehmite as an analog of the radioactive waste slurry to investigate the correlation between particle interactions, agglomerate morphology, and slurry rheology. EXPERIMENTS: A combination of Couette rheometry and small-angle scattering techniques (independently and simultaneously) were used to understand how agglomerate structure of slurry changes under flow and how these structural changes manifest themselves in the bulk rheology of the suspensions. FINDINGS: Our experiments show that the boehmite slurries are thixotropic, with the rheology and structure of the suspensions changing with increasing exposure to flow. In the slurries, particle agglomerates begin as loose, system-spanning clusters, but exposure to moderate shear rates causes the agglomerates to irreversibly consolidate into denser clusters of finite size. The structural changes directly influence the rheological properties of the slurries such as viscosity and viscoelasticity. Our study shows that solution pH affects the amount of structural rearrangement and the kinetics of the rearrangement process, with an increase in pH leading to faster and more dramatic changes in bulk rheology, which can be understood via correlations between particle interactions and the strength of particle network. Nearly identical structural changes were also observed in Poiseuille flow geometries, implying that the observed changes are relevant in pipe flow as well.

Effects of Ionic Strength, Salt, and pH on Aggregation of Boehmite Nanocrystals: Tumbler Small-Angle Neutron and X-ray Scattering and Imaging Analysis.

The US government currently spends significant resources managing the legacies of the Cold War, including 300 million liters of highly radioactive wastes stored in hundreds of tanks at the Hanford (WA) and Savannah River (SC) sites. The materials in these tanks consist of highly radioactive slurries and sludges at very high pH and salt concentrations. The solid particles primarily consist of aluminum hydroxides and oxyhydroxides (gibbsite and boehmite), although many other materials are present. These form complex aggregates that dramatically affect the rheology of the solutions and, therefore, efforts to recover and treat these wastes. In this paper, we have used a combination of transmission and cryo-transmission electron microscopy, dynamic light scattering, and X-ray and neutron small and ultrasmall-angle scattering to study the aggregation of synthetic nanoboehmite particles at pH 9 (approximately the point of zero charge) and 12, and sodium nitrate and calcium nitrate concentrations up to 1 m. Although the initial particles form individual rhombohedral platelets, once placed in solution they quickly form well-bonded stacks, primary aggregates, up to ∼1500 Å long. These are more prevalent at pH = 12. Addition of calcium nitrate or sodium nitrate has a similar effect as lowering pH, but approximately 100 times less calcium than sodium is needed to observe this effect. These aggregates have fractal dimension between 2.5 and 2.6 that are relatively unaffected by salt concentration for calcium nitrate at high pH. Larger aggregates (>∼4000 Å) are also formed, but their size distributions are discrete rather than continuous. The fractal dimensions of these aggregates are strongly pH-dependent, but only become dependent on solute at high concentrations.

Direction-specific interaction forces underlying zinc oxide crystal growth by oriented attachment.

Crystallization by particle attachment is impacting our understanding of natural mineralization processes and holds promise for novel materials design. When particles assemble in crystallographic alignment, expulsion of the intervening solvent and particle coalescence are enabled by near-perfect co-alignment via interparticle forces that remain poorly quantified. Here we report measurement and simulation of these nanoscale aligning forces for the ZnO(0001)-ZnO(000[Formula: see text]) system in aqueous solution. Dynamic force spectroscopy using nanoengineered single crystal probes reveals an attractive force with 60o rotational periodicity. Calculated distance and orientation-dependent potentials of mean force show several attractive free energy wells distinguished by numbers of intervening water layers, which reach a minimum when aligned. The calculated activation energy to separate the attractively bound solvated interfaces perfectly reproduces the measured 60o periodicity, revealing the key role of intervening water structuring as a basis to generate the interparticle torque that completes alignment and enables coalescence.Crystal growth is a fundamental process, important in a wide range of fields, but the interparticle forces responsible for molecule alignment are not well understood. Here, the authors measure the alignment forces in ZnO using dynamic force spectroscopy, highlighting the role of intervening water molecules.

The origin of facet selectivity and alignment in anatase TiO2 nanoparticles in electrolyte solutions: implications for oriented attachment in metal oxides.

Oriented attachment (OA) is an important nonclassical pathway for crystal growth from solution, occurring by the self-assembly of nanoparticles and often leading to highly organized three-dimensional crystal morphologies. The forces that drive nanocrystal reorientation for face-selective attachment and exclude improperly aligned particles have remained unknown. Here we report evidence at the microscopic level that ion correlation forces arising from dynamically interacting electrical double layers are responsible for face-selective attraction and particle rotation into lattice co-alignment as particles interact at long range. Atomic-to-mesoscale simulations developed and performed for the archetype OA system of anatase TiO2 nanoparticles in aqueous HCl solutions show that face-selective attraction from ion correlation forces outcompetes electrostatic repulsion at several nanometers apart, drawing particle face pairs into a metastable solvent-separated captured state. The analysis of the facet and pH dependence of interparticle interactions is in quantitative agreement with the observed decreasing frequency of attachment between the (112), (001), and (101) face pairs, revealing an adhesion barrier that is largely due to steric hydration forces from structured intervening solvents. This finding helps open new avenues for controlling crystal growth pathways leading to highly ordered three-dimensional nanomaterials.

Spontaneous water oxidation at hematite (α-Fe2O3) crystal faces.

Hematite (α-Fe2O3) persists as a promising candidate for photoelectrochemical water splitting, but a slow oxygen evolution reaction (OER) at its surfaces remains a limitation. Here we extend a series of studies that examine pH-dependent surface potentials and electron-transfer properties of effectively perfect low-index crystal faces of hematite in contact with simple electrolyte. Zero-resistance amperometry (ZRA) was performed in a two electrode configuration to quantify spontaneous dark current between hematite crystal face pairs (001)/(012), (001)/(113), and (012)/(113) at pH 3. Exponentially decaying currents initially of up to 200 nA were reported between faces over 4 min experiments. Fourth-order ZRA kinetics indicated rate limitation by the OER for current that flows between (001)/(012) and (001)/(113) face pairs, with the (012) and (113) faces serving as the anodes when paired with (001). The cathodic partner reaction is reductive dissolution of the (001) face, converting surface Fe(3+) to solubilized aqueous Fe(2+), at a rate maintained by the OER at the anode. In contrast, OER rate limitation does not manifest for the (012)/(113) pair. The uniqueness of the (001) face is established in terms of a faster intrinsic ability to accept the protons required for the reductive dissolution reaction. OER rate limitation inversely may thus arise from sluggish kinetics of hematite surfaces to dispense with the protons that accompany the four-electron OER. The results are explained in terms of semiquantitative energy band diagrams. The finding may be useful as a consideration for tailoring the design of polycrystalline hematite photoanodes that present multiple terminations to the interface with electrolyte.

Effect of surface site interactions on potentiometric titration of hematite (α-Fe2O3) crystal faces.

Time dependent potentiometric pH titrations were used to study the effect of atomic scale surface structure on the protonation behavior of the structurally well-defined hematite/electrolyte interfaces. Our recently proposed thermodynamic model [1,25] was applied to measured acidimetric and alkalimetric titration hysteresis loops, collected from highly organized (001), (012), and (113) crystal face terminations using pH equilibration times ranging from 15 to 30 min. Hysteresis loop areas indicate that (001) faces equilibrate faster than the (012) and (113) faces, consistent with the different expected ensembles of singly-, doubly-, and triply-coordinated surface sites on each face. Strongly non-linear hysteretic pH-potential relationships were found, with slopes exceeding Nernstian, collectively indicating that protonation and deprotonation is much more complex than embodied in present day surface complexation models. The asymmetrical shape of the acidimetric and alkalimetric titration branches were used to illustrate a proposed steric "leaky screen" repulsion/trapping interaction mechanism that stems from high affinity singly-coordinated sites electrostatically and sterically screening lower affinity doubly- and triply-coordinated sites. Our data indicate that site interaction is the dominant phenomenon defining surface potential accumulation behavior on single crystal faces of metal oxide minerals.

Synthesis and properties of titanomagnetite (Fe(3-x)Ti(x)O4) nanoparticles: a tunable solid-state Fe(II/III) redox system.

Titanomagnetite (Fe(3-x)Ti(x)O(4)) nanoparticles were synthesized by room temperature aqueous precipitation, in which Ti(IV) replaces Fe(III) and is charge compensated by conversion of Fe(III) to Fe(II) in the unit cell. A comprehensive suite of tools was used to probe composition, structure, and magnetic properties down to site-occupancy level, emphasizing distribution and accessibility of Fe(II) as a function of x. Synthesis of nanoparticles in the range 0≤x≤0.6 was attempted; Ti, total Fe and Fe(II) content were verified by chemical analysis. TEM indicated homogeneous spherical 9-12 nm particles. µ-XRD and Mössbauer spectroscopy on anoxic aqueous suspensions verified the inverse spinel structure and Ti(IV) incorporation in the unit cell up to x≤0.38, based on Fe(II)/Fe(III) ratio deduced from the unit cell edge and Mössbauer spectra. Nanoparticles with a higher value of x possessed a minor amorphous secondary Fe(II)/Ti(IV) phase. XANES/EXAFS indicated Ti(IV) incorporation in the octahedral sublattice (B-site) and proportional increases in Fe(II)/Fe(III) ratio. XA/XMCD indicated that increases arise from increasing B-site Fe(II), and that these charge-balancing equivalents segregate to those B-sites near particle surfaces. Dissolution studies showed that this segregation persists after release of Fe(II) into solution, in amounts systematically proportional to x and thus the Fe(II)/Fe(III) ratio. A mechanistic reaction model was developed entailing mobile B-site Fe(II) supplying a highly interactive surface phase that undergoes interfacial electron transfer with oxidants in solution, sustained by outward Fe(II) migration from particle interiors and concurrent inward migration of charge-balancing cationic vacancies in a ratio of 3:1.

Electrostatic potential of specific mineral faces.

Reaction rates of environmental processes occurring at hydrated mineral surfaces are in part controlled by the electrostatic potential that develops at the interface. This potential depends on the structure of exposed crystal faces as well as the pH and the type of ions and their interactions with these faces. Despite its importance, experimental methods for determining fundamental electrostatic properties of specific crystal faces such as the point of zero charge are few. Here we show that this information may be obtained from simple, cyclic potentiometric titration using a well-characterized single-crystal electrode exposing the face of interest. The method exploits the presence of a hysteresis loop in the titration measurements that allows the extraction of key electrostatic descriptors using the Maxwell construction. The approach is demonstrated for hematite (α-Fe(2)O(3)) (001), and thermodynamic proof is provided for the resulting estimate of its point of zero charge. Insight gained from this method will aid in predicting the fate of migrating contaminants, mineral growth/dissolution processes, and mineral-microbiological interactions and in testing surface complexation theories.

A high-pressure atomic force microscope for imaging in supercritical carbon dioxide.

A high-pressure atomic force microscope (AFM) that enables in situ, atomic scale measurements of topography of solid surfaces in contact with supercritical CO(2) (scCO(2)) fluids has been developed. This apparatus overcomes the pressure limitations of the hydrothermal AFM and is designed to handle pressures up to 100 atm at temperatures up to ∼350 K. A standard optically-based cantilever deflection detection system was chosen. When imaging in compressible supercritical fluids such as scCO(2), precise control of pressure and temperature in the fluid cell is the primary technical challenge. Noise levels and imaging resolution depend on minimization of fluid density fluctuations that change the fluid refractive index and hence the laser path. We demonstrate with our apparatus in situ atomic scale imaging of a calcite (CaCO(3)) mineral surface in scCO(2); both single, monatomic steps and dynamic processes occurring on the (1014) surface are presented. This new AFM provides unprecedented in situ access to interfacial phenomena at solid-fluid interfaces under pressure.

Bond paths and van der Waals interactions in orpiment, As2S3.

The calculated electron density distribution for orpiment, As(2)S(3), reveals that As-S, S-S, and As-As bond paths are associated with the experimental interlayer directed bonded interactions detected in a combined infrared and Raman study. The successful modeling of the infrared- and Raman-determined interlayer bonded interactions together with bond paths and the structuralization of a variety of inorganic molecules in terms of "key-lock" bond path mainstays support the argument that van der Waals forces in inorganic molecular crystals are directional.

Role of directed van der Waals bonded interactions in the determination of the structures of molecular arsenate solids.

Bond paths, local energy density properties, and Laplacian, L(r) = -wedge(2)rho(r), composite isosurfaces of the electron density distributions were calculated for the intramolecular and intermolecular bonded interactions for molecular solids of As(2)O(3) and AsO(2) composition, an As(2)O(5) crystal, a number of arsenate molecules, and the arsenic metalloid, arsenolamprite. The directed intermolecular van der Waals As-O, O-O, and As-As bonded interactions are believed to serve as mainstays between the individual molecules in each of the molecular solids. As-O bond paths between the bonded atoms connect Lewis base charge concentrations and Lewis acid charge depletion domains, whereas the O-O and As-As paths connect Lewis base pair and Lewis acid pair domains, respectively, giving rise to sets of intermolecular directed bond paths. The alignment of the directed bond paths results in the periodic structures adopted by the arsenates. The arrangements of the As atoms in the claudetite polymorphs of As(2)O(3) and the As atoms in arsenolamprite are similar. Like the As(2)O(3) polymorphs, arsenolamprite is a molecular solid connected by relatively weak As-As intermolecular directed van der Waals bond paths between the layers of stronger As-As intramolecular bonded interactions. The bond critical point and local energy density properties of the intermolecular As-As bonded interactions in arsenolamprite are comparable with the As-As interactions in claudetite I. As such, the structure of claudetite I can be viewed as a stuffed derivative of the arsenolamprite structure with O atoms between pairs of As atoms comprising the layers of the structure. The cubic structure adopted by the arsenolite polymorph can be understood in terms of sets of directed acid-base As-O and base-base O-O pair domains and bond paths that radiate from the tetrahedral faces of its constituent molecules, serving as face-to-face key-lock mainstays in forming a periodic tetrahedral array of molecules rather than one based on some variant of close packing. The relatively dense structure and the corrugation of the layers in claudetite I can also be understood in terms of directed van der Waals As-O bonded interactions. Our study not only provides a new basis for understanding the crystal chemistry and the structures of the arsenates, but it also calls for a reappraisal of the concept of van der Waals bonded interactions, how the structures of molecular solids are viewed, and how the molecules in these solids are bonded in a periodic structure.

Experimental bond critical point and local energy density properties determined for Mn-O, Fe-O, and Co-O bonded interactions for tephroite, Mn2SiO4, fayalite, Fe2SiO4, and Co2SiO4 olivine and selected organic metal complexes: comparison with properties calculated for non-transition and transition metal M-O bonded interactions for silicates and oxides.

Bond critical point (bcp) and local energy density properties for the electron density (ED) distributions, calculated with first-principle quantum mechanical methods for divalent transition metal Mn-, Co-, and Fe-containing silicates and oxides are compared with experimental model ED properties for tephroite, Mn 2SiO 4, fayalite, Fe 2SiO 4, and Co 2SiO 4 olivine, each determined with high-energy synchrotron single-crystal X-ray diffraction data. Trends between the experimental bond lengths, R(M-O), (M = Mn, Fe, Co), and the calculated bcp properties are comparable with those observed for non-transition M-O bonded interactions. The bcp properties, local total energy density, H( r c), and bond length trends determined for the Mn-O, Co-O, and Fe-O interactions are also comparable. A comparison is also made with model experimental bcp properties determined for several Mn-O, Fe-O, and Co-O bonded interactions for selected organometallic complexes and several oxides. Despite the complexities of the structures of the organometallic complexes, the agreement between the calculated and model experimental bcp properties is fair to good in several cases. The G( r c)/rho( r c) versus R(M-O) trends established for non-transition metal M-O bonded interactions hold for the transition metal M-O bonded interactions with G( r c)/rho( r c) increasing in value as H( r c) becomes progressively more negative in value, indicating an increasing shared character of the interaction as G( r c)/rho( r c) increases in value. As observed for the non-transition metal M-O bonded interactions, the Laplacian, nabla (2)rho( r c), increases in value as rho( r c) increases and as H( r c) decreases and becomes progressive more negative in value. The Mn-O, Fe-O, and Co-O bonded interactions are indicated to be of intermediate character with a substantial component of closed-shell character compared with Fe-S and Ni-S bonded interactions, which show greater shared character based on the | V( r c)|/ G( r c) bond character indicator. The atomic charges conferred on the transition metal atoms for the three olivines decrease with increasing atomic number from Mn to Fe to Co as the average M-O bond lengths decrease from 2.219 to 2.168 to 2.128 A, respectively.

Shared and closed-shell O-O interactions in silicates.

Bond paths of maximum electron density spanning O-O edges shared between equivalent or quasiequivalent MOn (n > 4) coordination polyhedra are not uncommon electron density features displayed by silicates. On the basis of the positive values for the local electronic energy density, H(rc), at the bond critical points, rc, they qualify as weak "closed-shell" interactions. As observed for M-O bonded interactions (M = first and second row metal atoms), the electron density, rho(rc), and the Laplacian of the electron density increase in a regular way as the separation between the O atoms, R(O-O), decreases. A simple model, based on R(O-O) and the distances of the Si atoms from the midpoint between adjacent pairs of O atoms, partitions the O-O bond paths in the high-pressure silica polymorph coesite into two largely disjoint domains, one with and one without bond paths. The occurrence of O-O bond paths shared in common between equivalent coordination polyhedra suggests that they may be grounded in some cases on factors other than bonded interactions, particularly since they are often displayed by inert procrystal representations of the electron density. In these cases, it can be argued that the accumulation of the electron density along the paths has its origin, at least in part, in the superposition of the peripheral electron density distributions of the metal M atoms occupying the edge-sharing polyhedra. On the other hand, the accumulation of electron density along the paths may stabilize a structure by shielding the adjacent M atoms in the edge-sharing polyhedra. For closed-shell Li-O, Na-O, and Mg-O interactions, H(rc) is positive and increases as the value of rho(rc) increases, unlike the "shared" Be-O, B-O, C-O, Al-O, Si-O, P-O, and S-O interactions, where H(rc) is negative and decreases as rho(rc) increases. The H(rc) values for the weak closed-shell O-O interactions also increase as rho(rc) increases, as observed for the closed-shell M-O interactions. On the basis of the bond critical point properties and the negative H(rc) value, the O-O interaction comprising the O2 molecule in silica III qualifies as a shared interaction.

Theoretical electron density distributions for Fe- and Cu-sulfide earth materials: a connection between bond length, bond critical point properties, local energy densities, and bonded interactions.

Bond critical point and local energy density properties together with net atomic charges were calculated for theoretical electron density distributions, rho(r), generated for a variety of Fe and Cu metal-sulfide materials with high- and low-spin Fe atoms in octahedral coordination and high-spin Fe atoms in tetrahedral coordination. The electron density, rho(rc), the Laplacian, triangle down2rho(rc), the local kinetic energy, G(rc), and the oxidation state of Fe increase as the local potential energy density, V(rc), the Fe-S bond lengths, and the coordination numbers of the Fe atoms decrease. The properties of the bonded interactions for the octahedrally coordinated low-spin Fe atoms for pyrite and marcasite are distinct from those for high-spin Fe atoms for troilite, smythite, and greigite. The Fe-S bond lengths are shorter and the values of rho(rc) and triangle down2rho(rc) are larger for pyrite and marcasite, indicating that the accumulation and local concentration of rho(r) in the internuclear region are greater than those involving the longer, high-spin Fe-S bonded interactions. The net atomic charges and the bonded radii calculated for the Fe and S atoms in pyrite and marcasite are also smaller than those for sulfides with high-spin octahedrally coordinated Fe atoms. Collectively, the Fe-S interactions are indicated to be intermediate in character with the low-spin Fe-S interactions having greater shared character than the high-spin interactions. The bond lengths observed for chalcopyrite together with the calculated bond critical point properties are consistent with the formula Cu+Fe3+S2. The bond length is shorter and the rho(rc) value is larger for the FeS4 tetrahedron displayed by metastable greigite than those displayed by chalcopyrite and cubanite, consistent with a proposal that the Fe atom in greigite is tetravalent. S-S bond paths exist between each of the surface S atoms of adjacent slabs of FeS6 octahedra comprising the layer sulfide smythite, suggesting that the neutral Fe3S4 slabs are linked together and stabilized by the pathways of electron density comprising S-S bonded interactions. Such interactions not only exist between the S atoms for adjacent S8 rings in native sulfur, but their bond critical point properties are similar to those displayed by the metal sulfides.

Si-O bonded interactions in silicate crystals and molecules: a comparison.

Bond critical point, local kinetic energy density, G(rc), and local potential energy density, V(rc), properties of the electron density distributions, rho(r), calculated for silicates such as quartz and gas-phase molecules such as disiloxane are similar, indicating that the forces that govern the Si-O bonded interactions in silica are short-ranged and molecular-like. Using the G(rc)/rho(rc) ratio as a measure of bond character, the ratio increases as the Si-O bond length, the local electronic energy density, H(rc)= G(rc) + V(rc), and the coordination number of the Si atom decrease and as the accumulation of the electron density at the bond critical point, rho(rc), and the Laplacian, inverted Delta2 rho(rc), increase. The G(rc)/rho(rc) and H(rc)/rho(rc) ratios categorize the bonded interaction as observed for other second row atom M-O bonds into discrete categories with the covalent character of each of the M-O bonds increasing with the H(rc)/rho(rc) ratio. The character of the bond is examined in terms of the large net atomic charges conferred on the Si atoms comprising disiloxane, stishovite, quartz, and forsterite and the domains of localized electron density along the Si-O bond vectors and on the reflex side of the Si-O-Si angle together with the close similarity of the Si-O bonded interactions observed for a variety of hydroxyacid silicate molecules and a large number of silicate crystals. The bond critical point and local energy density properties of the electron density distribution indicate that the bond is an intermediate interaction between Al-O and P-O bonded interactions rather than being a closed-shell or a shared interaction.

Bond length and local energy density property connections for non-transition-metal oxide-bonded interactions.

For a variety of molecules and earth materials, the theoretical local kinetic energy density, G(r(c)), increases and the local potential energy density, V(r(c)), decreases as the M-O bond lengths (M = first- and second-row metal atoms bonded to O) decrease and the electron density, rho(r(c)), accumulates at the bond critical points, r(c). Despite the claim that the local kinetic energy density per electronic charge, G(r(c))/rho(r(c)), classifies bonded interactions as shared interactions when less than unity and closed-shell when greater, the ratio was found to increase from 0.5 to 2.5 au as the local electronic energy density, H(r(c)) = G(r(c)) + V(r(c)), decreases and becomes progressively more negative. The ratio appears to be a measure of the character of a given M-O bonded interaction, the greater the ratio, the larger the value of rho(r(c)), the smaller the coordination number of the M atom and the more shared the bonded interaction. H(r(c))/rho(r(c)) versus G(r(c))/rho(r(c)) scatter diagrams categorize the M-O bonded interactions into domains with the local electronic energy density per electron charge, H(r(c))/rho(r(c)), tending to decrease as the electronegativity differences for the bonded pairs of atoms decrease. The values of G(r(c)) and V(r(c)), estimated with a gradient-corrected electron gas theory expression and the local virial theorem, are in good agreement with theoretical values, particularly for the bonded interactions involving second-row M atoms. The agreement is poorer for shared C-O and N-O bonded interactions.

Classification of metal-oxide bonded interactions based on local potential- and kinetic-energy densities.

A classification of the hydrogen fluoride H-F-bonded interactions comprising a large number of molecules has been proposed by Espinosa et al. [J. Chem. Phys. 117, 5529 (2002)] based on the ratio /Vr(c)/ / Gr(c) where /Vr(c)/ is the magnitude of the local potential-energy density and Gr(c) is the local kinetic-energy density, each evaluated at a bond critical point r(c). A calculation of the ratio for the M-O bonded interactions comprising a relatively large number of oxide molecules and earth materials, together with the constraints imposed by the values of inverted Delta2rho r(c) and the local electronic energy density, Hr(c) = Gr(c) + Vr(c), in the H-F study, yielded practically the same classification for the oxides. This is true despite the different trends that hold between the bond critical point and local energy density properties with the bond lengths displayed by the H-F and M-O bonded interactions. On the basis of the ratio, Li-O, Na-O, and Mg-O bonded interactions classify as closed-shell ionic bonds, Be-O, Al-O, Si-O, B-O, and P-O interactions classify as bonds of intermediate character with the covalent character increasing from Be-O to P-O. N-O interactions classify as shared covalent bonds. C-O and S-O bonded interactions classify as both intermediate and covalent bonded interactions. The C-O double- and triple-bonded interactions classify as intermediate-bonded interactions, each with a substantial component of covalent character and the C-O single-bonded interaction classifies as a covalent bond whereas their local electronic energy density values indicate that they are each covalent bonded interactions. The ratios for the Be-O, Al-O, and Si-O bonded interactions indicate that they have a substantial component of ionic character despite their classification as bonds of intermediate character. The trend between the ratio and the character of the bonded interactions is consistent with trends expected from electronegativity considerations. The ratio increases as the net charges and the coordination numbers for the atoms for several Ni-sulfides decrease. On the contrary, the ratio for the Si-O bonded interactions for the orthosilicate, forsterite, Mg2SiO4, and the high-pressure silica polymorph, stishovite, decreases as the observed net atomic charges and the coordination numbers of Si and O increase in value. The ratio for the Ni-Ni bonded interactions for the Ni-sulfides and bulk Ni metal indicate that the interactions are intermediate in character with a substantial component of ionic character.

Theoretical characterization of charge transport in chromia (alpha-Cr2O3).

Transport of conduction electrons and holes through the lattice of alpha-Cr2O3 (chromia) is modeled as a valence alternation of chromium cations using ab initio electronic structure calculations and electron-transfer theory. In the context of the small polaron model, a cluster approach was used to compute quantities controlling the mobility of localized electrons and holes, i.e., the reorganization energy and the electronic coupling matrix element that enter Marcus' theory. The calculation of the electronic coupling followed the generalized Mulliken-Hush approach using the complete active space self-consistent-field (CASSCF) method and the quasidiabatic method. Our findings indicate that hole mobility is more than three orders of magnitude larger than electron mobility in both (001) and [001] lattice directions. The difference arises mainly from the larger internal reorganization energy calculated for electron-transport relative to hole-transport processes while electronic couplings have similar magnitudes. The much larger hole mobility versus electron mobility in alpha-Cr2O3 is in contrast to similar hole and electron mobilities in hematite alpha-Fe2O3 previously calculated. Our calculations also indicate that the electronic coupling for all charge-transfer processes of interest is smaller than for the corresponding processes in hematite. This variation is attributed to the weaker interaction between the metal 3d states and the O(2p) states in chromia than in hematite, leading to a smaller overlap between the charge-transfer donor and acceptor wave functions and smaller superexchange coupling in chromia. Nevertheless, the weaker coupling in chromia is still sufficiently large to suggest that charge-transport processes in chromia are adiabatic in nature. The electronic coupling is found to depend on both the superexchange interaction through the bridging oxygen atoms and the d-shell electron-spin coupling within the Cr-Cr donor-acceptor pair, while the reorganization energy is essentially independent of the electron-spin coupling.

Charge transport in metal oxides: a theoretical study of hematite alpha-Fe2O3.

Transport of conduction electrons and holes through the lattice of alpha-Fe(2)O(3) (hematite) is modeled as a valence alternation of iron cations using ab initio electronic structure calculations and electron transfer theory. Experimental studies have shown that the conductivity along the (001) basal plane is four orders of magnitude larger than the conductivity along the [001] direction. In the context of the small polaron model, a cluster approach was used to compute quantities controlling the mobility of localized electrons and holes, i.e., the reorganization energy and the electronic coupling matrix element that enter Marcus' theory. The calculation of the electronic coupling followed the generalized Mulliken-Hush approach using the complete active space self-consistent field method. Our findings demonstrate an approximately three orders of magnitude anisotropy in both electron and hole mobility between directions perpendicular and parallel to the c axis, in good accord with experimental data. The anisotropy arises from the slowness of both electron and hole mobilities across basal oxygen planes relative to that within iron bilayers between basal oxygen planes. Interestingly, for elementary reaction steps along either of the directions considered, there is only less than one order of magnitude difference in mobility between electrons and holes, in contrast to accepted classical arguments. Our findings indicate that the most important quantity underlying mobility differences is the electronic coupling, albeit the reorganization energy contributes as well. The large values computed for the electronic coupling suggest that charge transport reactions in hematite are adiabatic in nature. The electronic coupling is found to depend on both the superexchange interaction through the bridging oxygen atoms and the d-shell electron spin coupling within the Fe-Fe donor-acceptor pair, while the reorganization energy is essentially independent of the electron spin coupling.