Surface structure of cubic aluminum oxide

Atomistic Simulations of the Crystalline-to-Amorphous Transformation of γ-Al2O3 Nanoparticles: Delicate Interplay between Lattice Distortions, Stresses, and Space Charges

The size-dependent phase stability of γ-Al₂O₃ was studied by large-scale molecular dynamics simulations over a wide temperature range from 300 to 900 K. For the γ-Al₂O₃ crystal, a bulk transformation to α-Al₂O₃ by an FCC-to-HCP transition of the O sublattice is still kinetically hindered at 900 K. However, local distortions of the FCC O-sublattice by the formation of quasi-octahedral Al local coordination spheres become thermally activated, as driven by the partial covalency of the Al-O bond. On the contrary, spherical γ-Al₂O₃ nanoparticles (NPs) (with sizes of 6 and 10 nm) undergo a crystalline-to-amorphous transformation at 900 K, which starts at the reconstructed surface and propagates into the core through collective displacements of anions and cations, resulting in the formation of 7- and 8-fold local coordination spheres of Al. In parallel, the reconstructed Al-enriched surface is separated from the stoichiometric core by a diffuse Al-depleted transition region. This compositional heterogeneity creates an imbalance of charges inside the NP, which induces a net attractive Coulombic force that is strong enough to reverse the initial stress state in the NP core from compressive to tensile. These findings disclose the delicate interplay between lattice distortions, stresses, and space-charge regions in oxide nanosystems. A fundamental explanation for the reported expansion of metal-oxide NPs with decreasing size is provided, which has significant implications for, e.g., heterogeneous catalysis, NP sintering, and additive manufacturing of NP-reinforced metal matrix composites.

Super-Hydrophilic Microporous Ni(OH)x/Ni3 S2 Heterostructure Electrocatalyst for Large-Current-Density Hydrogen Evolution

Exploiting active and stable non-precious metal electrocatalysts for alkaline hydrogen evolution reaction (HER) at large current density plays a key role in realizing large-scale industrial hydrogen generation. Herein, a self-supported microporous Ni(OH)x/Ni₃ S₂ heterostructure electrocatalyst on nickel foam (Ni(OH)x/Ni₃ S₂ /NF) that possesses super-hydrophilic property through an electrochemical process is rationally designed and fabricated. Benefiting from the super-hydrophilic property, microporous feature, and self-supported structure, the electrocatalyst exhibits an exceptional HER performance at large current density in 1.0 M KOH, only requiring low overpotential of 126, 193, and 238 mV to reach a current density of 100, 500, and 1000 mA cm^-2 , respectively, and displaying a long-term durability up to 1000 h, which is among the state-of-the-art non-precious metal electrocatalysts. Combining hard X-rays absorption spectroscopy and first-principles calculation, it also reveals that the strong electronic coupling at the interface of the heterostructure facilitates the dissociation of H₂ O molecular, accelerating the HER kinetics in alkaline electrolyte. This work sheds a light on developing advanced non-precious metal electrocatalysts for industrial hydrogen production by means of constructing a super-hydrophilic microporous heterostructure.

Hydration Structures on γ-Alumina Surfaces With and Without Electrolytes Probed by Atomistic Molecular Dynamics Simulations

A wide range of systems, both engineered and natural, feature aqueous electrolyte solutions at interfaces. In this study, the structure and dynamics of water at the two prevalent crystallographic terminations of gamma-alumina, [110] and [100], and the influence of salts─sodium chloride, ammonium acetate, barium acetate, and barium nitrate on such properties─were investigated using equilibrium molecular dynamics simulations. The resulting interfacial phenomena were quantified from simulation trajectories via atomic density profiles, angle probability distributions, residence times, 2-D density distributions within the hydration layers, and hydrogen bond density profiles. Analysis and interpretation of the results are supported by simulation snapshots. Taken together, our results show stronger interaction and closer association of water with the [110] surface, compared to [100], while ion-induced disruption of interfacial water structure was more prevalent at the [100] surface. For the latter, a stronger association of cations is observed, namely sodium and ammonium, and ion adsorption appears determined by their size. The differences in surface-water interactions between the two terminations are linked to their respective surface features and distributions of surface groups, with atomistic-scale roughness of the [110] surface promoting closer association of interfacial water. The results highlight the fundamental role of surface characteristics in determining surface-water interactions, and the resulting effects on ion-surface and ion-water interactions. Since the two terminations of gamma-alumina considered represent interfaces of significance to numerous industrial applications, the results provide insights relevant for catalyst preparation and adsorption-based water treatment, among other applications.

Physical Properties of γ Alumina Surface Hydroxyls Revisited through a Large Scale Periodic Quantum-Chemistry Approach

We have studied surface hydroxyls adsorbed onto (001), (011), and (111) gamma alumina surfaces using a quantum-chemistry approach in order to compare with empirical models proposed in the literature. Local electronic structures and geometries in the low OH coverage limit have been evaluated for both ideal and relaxed surfaces with the help of a large scale periodic quantum-chemical code. Hydroxyl groups are adsorbed onto surfaces, and a study of their local electronic structure, vibrational frequencies, charges, and adsorption energies is performed and analyzed as a function of their adsorption site geometry. Our results show that, even on ideal (nonrelaxed) surfaces, OH local environments are more complicated than those stated by empirical models and strongly influence the hydroxyl stretching vibrational mode. Large scale simulation shows that disorder takes place even at 0 K, and the analysis of the vibrational frequencies leads to a revision of Knözinger's empirical model. Cationic vacancies in the first surface layers have also been taken into account; they have a significant influence on the surface atomic and electronic structures, modifying the physical properties of adsorbed OH entities. This work emphasizes the necessity to perform an electronic structure calculation to better understand adsorbed OH properties on gamma alumina surfaces and reveals the difficulty to make a one-to-one correspondence between OH stretching frequencies and their other physical properties. Finally, we show that these results agree with some available experimental studies.

Nonequilibrium energy dissipation at the interface of sliding model hydroxylated α-alumina surfaces

Nonequilibrium molecular dynamics simulations were performed to study the dynamics of energy transfer at the interface of a small nanoscale hydroxylated alpha-alumina surface sliding across a much larger surface of the same material. Sliding velocities of 0.05, 0.5, 5, and 50 ms and loads of 0, 0.0625, 5, 15, 25, and 100 nN were considered. Nonequilibrium energy distributions were found at the interface for each of these conditions. The velocity distribution P(v) for the atoms in a sublayer of the smaller surface oscillates during the sliding, reflecting the periodicity of the interfacial intermolecular potential. When averaged over the sliding, this P(v) for each of the sublayers is bimodal with Boltzmann and non-Boltzmann components. The non-Boltzmann component, with temperatures in excess of 1000 K and as high as 2500 K, is most important for the interfacial H-atom sublayer and becomes less important in moving to a sublayer further from the interface. Similarly, the temperature of the Boltzmann component decreases for sublayers further from the interface and approaches the 300 K temperature of the boundary. The temperature of the Boltzmann component decreases, but the importance of the non-Boltzmann component increases, as the sliding velocity is decreased. The temperature of the non-Boltzmann component is relatively insensitive to the sliding velocity. Friction forces are determined by calculating the energy dissipation during the sliding, and different regimes are found for variation in the friction force versus sliding velocity v(s) and applied load. For v(s) of 0.05, 0.5, and 5 ms, the friction force is inversely proportional to v(s) reflecting the increased time for energy dissipation as v(s) is decreased.