Michael Ashton

Behavior of molecules and molecular ions near a field emitter

The cold emission of particles from surfaces under intense electric fields is a process which underpins a variety of applications including atom probe tomography (APT), an analytical microscopy technique with near-atomic spatial resolution. Increasingly relying on fast laser pulsing to trigger the emission, APT experiments often incorporate the detection of molecular ions emitted from the specimen, in particular from covalently or ionically bonded materials. Notably, it has been proposed that neutral molecules can also be emitted during this process. However, this remains a contentious issue. To investigate the validity of this hypothesis, a careful review of the literature is combined with the development of new methods to treat experimental APT data, the modelling of ion trajectories, and the application of density-functional theory (DFT) simulations to derive molecular ion energetics. It is shown that the direct thermal emission of neutral molecules is extremely unlikely. However, neutrals can still be formed in the course of an APT experiment by dissociation of metastable molecular ions.

Topology-Scaling Identification of Layered Solids and Stable Exfoliated 2D Materials

The Materials Project crystal structure database has been searched for materials possessing layered motifs in their crystal structures using a topology-scaling algorithm. The algorithm identifies and measures the sizes of bonded atomic clusters in a structures unit cell, and determines their scaling with cell size. The search yielded 826 stable layered materials that are considered as candidates for the formation of two-dimensional monolayers via exfoliation. Density-functional theory was used to calculate the exfoliation energy of each material and 680 monolayers emerge with exfoliation energies below those of already-existent two-dimensional materials. The crystal structures of these two-dimensional materials provide templates for future theoretical searches of stable two-dimensional materials. The optimized structures and other calculated data for all 826 monolayers are provided at our database (https://materialsweb.org).

Computational characterization of lightweight multilayer MXene Li-ion battery anodes

MXenes, a class of two-dimensional transition metal carbides and nitrides, have shown promise experimentally and computationally for use in energy storage applications. In particular, the most lightweight members of the monolayer MXene family (M = Sc, Ti, V, or Cr) are predicted to have gravimetric capacities above 400 mAh/g, higher than graphite. Additionally, intercalation of ions into multilayer MXenes can be accomplished at low voltages, and low diffusion barriers exist for Li diffusing across monolayer MXenes. However, large discrepancies have been observed between the calculated and experimental reversible capacities of MXenes. Here, dispersion-corrected density functional theory calculations are employed to predict reversible capacities and other battery-related properties for six of the most promising members of the MXene family (O-functionalized Ti- and V-based carbide MXenes) as bilayer structures. The calculated reversible capacities of the V2CO2 and Ti2CO2 bilayers agree more closely with experiment than do previous calculations for monolayers. Additionally, the minimum energy paths and corresponding energy barriers along the in-plane [1000] and [0100] directions for Li travelling between neighboring MXene layers are determined. V4C3O2 exhibits the lowest diffusion barrier of the compositions considered, at 0.42 eV, but its reversible capacity (148 mAh/g) is dragged down by its heavy formula unit. Conversely, the V2CO2 MXene shows good reversible capacity (276 mAh/g), but a high diffusion barrier (0.82 eV). We show that the diffusion barriers of all bilayer structures are significantly higher than those calculated for the corresponding monolayers, advocating the use of dispersed monolayer MXenes instead of multilayers in high performance anodes.

Two-Dimensional Intrinsic Half-Metals With Large Spin Gaps

Through a systematic search of all layered bulk compounds combined with density functional calculations employing hybrid exchange-correlation functionals, we predict a family of three magnetic two-dimensional (2D) materials with half-metallic band structures. The 2D materials, FeCl2, FeBr2, and FeI2, are all sufficiently stable to be exfoliated from bulk layered compounds. The Fe2+ ions in these materials are in a high-spin octahedral d6 configuration leading to a large magnetic moment of 4 μB. Calculations of the magnetic anisotropy show an easy-plane for the magnetic moment. A classical XY model with nearest neighbor coupling estimates critical temperatures, Tc, for the Berezinskii-Kosterlitz-Thouless transition ranging from 122 K for FeI2 to 210 K for FeBr2. The quantum confinement of these 2D materials results in unusually large spin gaps, ranging from 4.0 eV for FeI2 to 6.4 eV for FeCl2, which should defend against spin current leakage even at small device length scales. Their purely spin-polarized currents and dispersive interlayer interactions should make these materials useful for 2D spin valves and other spintronic applications.

Computational discovery and characterization of polymorphic two-dimensional IV–V materials

First-principles calculations predict the stability and properties of two-dimensional (2D) group IV–V MX compounds (M = Si, Ge, Sn, Pb; X = P, As, Sb, Bi) in two distinct crystallographic configurations: a low symmetry Cm layer inspired by the layered bulk structures of SiP, SiAs, GeP, and GeAs, and a high symmetry P6¯m2 structure previously predicted for 2D SiP. The calculations predict that the Cm structure is more stable for X = As, Sb, and Bi, and the P6¯m2 structure for X = P. The electronic properties of the two structures are quite different. Specifically, the Cm band gaps are systematically larger by about 15% for most compounds, and the band offsets indicate that all of these compounds in the Cm crystal structure will form type II heterojunctions in contact with their P6¯m2 polymorphs. Pourbaix diagrams predict that a few of the 2D IV–V compounds are stable when exposed to water.

Computational Study of Low Interlayer Friction in Tin+1Cn (n = 1, 2, and 3) MXene

The friction of adjacent Tin+1Cn (n = 1, 2, and 3) MXene layers is investigated using density functional theory (DFT) calculations and classical molecular dynamics simulations with ReaxFF potentials. The calculations reveal the sliding pathways in all three MXene systems with low energy barriers. The friction coefficients for interlayer sliding are evaluated using static calculations. Both DFT and ReaxFF methods predict friction coefficients between 0.24 and 0.27 for normal loads less than 1.2 GPa. The effect of titanium (Ti) vacancies in sublayers and terminal oxygen (O) vacancies at surfaces on the interlayer friction is further investigated using the ReaxFF potential. These defects are found to increase the friction coefficients by increasing surface roughness and creating additional attractive forces between adjacent layers. However, these defective MXenes still maintain friction coefficients below 0.31. We also consider functionalized Ti3C2 MXene terminated with -OH and -OCH3 and find that compared to the -O-terminated surface both groups further reduce the interlayer friction coefficient to 0.10-0.14.