Ab initio Description of Bond Breaking in Large Electric Fields

Dielectric Properties of Nanoconfined Water from Ab Initio Thermopotentiostat Molecular Dynamics

We discuss how to include our recently proposed thermopotentiostat technique [Deissenbeck et al. Phys. Rev. Lett. 2021, 126, 136803] into any existing ab initio molecular dynamics (AIMD) package. Using thermopotentiostat AIMD simulations in the canonical NVTΦ ensemble at a constant electrode potential, we compute the polarization bound charge and dielectric response of interfacial water from first principles.

Spontaneous Water Radical Cation Oxidation at Double Bonds in Microdroplets

Spontaneous oxidation of compounds containing diverse X=Y moieties (e.g., sulfonamides, ketones, esters, sulfones) occurs readily in organic-solvent microdroplets. This surprising phenomenon is proposed to be driven by the generation of an intermediate species [M+H₂O]^+·: a covalent adduct of water radical cation (H₂O^+·) with the reactant molecule (M). The adduct is observed in the positive ion mass spectrum while its formation in the interfacial region of the microdroplet (i.e., at the air-droplet interface) is indicated by the strong dependence of the oxidation product formation on the spray distance (which reflects the droplet size and consequently the surface-to-volume ratio) and the solvent composition. Importantly, based on the screening of a ca. 21,000-compound library and the detailed consideration of six functional groups, the formation of a molecular adduct with the water radical cation is a significant route to ionization in positive ion mode electrospray, where it is favored in those compounds with X=Y moieties which lack basic groups. A set of model monofunctional systems was studied and in one case, benzyl benzoate, evidence was found for oxidation driven by hydroxyl radical adduct formation followed by protonation in addition to the dominant water radical cation addition process. Significant implications of molecular ionization by water radical cations for oxidation processes in atmospheric aerosols, analytical mass spectrometry and small-scale synthesis are noted.

Status and Direction of Atom Probe Analysis of Frozen Liquids

Imaging of liquids and cryogenic biological materials by electron microscopy has been recently enabled by innovative approaches for specimen preparation and the fast development of optimized instruments for cryo-enabled electron microscopy (cryo-EM). Yet, cryo-EM typically lacks advanced analytical capabilities, in particular for light elements. With the development of protocols for frozen wet specimen preparation, atom probe tomography (APT) could advantageously complement insights gained by cryo-EM. Here, we report on different approaches that have been recently proposed to enable the analysis of relatively large volumes of frozen liquids from either a flat substrate or the fractured surface of a wire. Both allowed for analyzing water ice layers which are several micrometers thick consisting of pure water, pure heavy water, and aqueous solutions. We discuss the merits of both approaches and prospects for further developments in this area. Preliminary results raise numerous questions, in part concerning the physics underpinning field evaporation. We discuss these aspects and lay out some of the challenges regarding the APT analysis of frozen liquids.

Can electric fields drive chemistry for an aqueous microdroplet?

Reaction rates of common organic reactions have been reported to increase by one to six orders of magnitude in aqueous microdroplets compared to bulk solution, but the reasons for the rate acceleration are poorly understood. Using a coarse-grained electron model that describes structural organization and electron densities for water droplets without the expense of ab initio methods, we investigate the electric field distributions at the air-water interface to understand the origin of surface reactivity. We find that electric field alignments along free O-H bonds at the surface are ~16 MV/cm larger on average than that found for O-H bonds in the interior of the water droplet. Furthermore, electric field distributions can be an order of magnitude larger than the average due to non-linear coupling of intramolecular solvent polarization with intermolecular solvent modes which may contribute to even greater surface reactivity for weakening or breaking chemical bonds at the droplet surface.

Reflections on the Spatial Performance of Atom Probe Tomography in the Analysis of Atomic Neighborhoods

Atom probe tomography (APT) is often introduced as providing “atomic-scale” mapping of the composition of materials and as such is often exploited to analyze atomic neighborhoods within a material. Yet quantifying the actual spatial performance of the technique in a general case remains challenging, as it depends on the material system being investigated as well as on the specimen's geometry. Here, by using comparisons with field-ion microscopy experiments, field-ion imaging and field evaporation simulations, we provide the basis for a critical reflection on the spatial performance of APT in the analysis of pure metals, low alloyed systems and concentrated solid solutions (i.e., akin to high-entropy alloys). The spatial resolution imposes strong limitations on the possible interpretation of measured atomic neighborhoods, and directional neighborhood analyses restricted to the depth are expected to be more robust. We hope this work gets the community to reflect on its practices, in the same way, it got us to reflect on our work.

Aluminum depletion induced by co-segregation of carbon and boron in a bcc-iron grain boundary

The local variation of grain boundary atomic structure and chemistry caused by segregation of impurities influences the macroscopic properties of polycrystalline materials. Here, the effect of co-segregation of carbon and boron on the depletion of aluminum at a Σ5 (3 1 0 )[0 0 1] tilt grain boundary in a α - Fe-4 at%Al bicrystal is studied by combining atomic resolution scanning transmission electron microscopy, atom probe tomography and density functional theory calculations. The atomic grain boundary structural units mostly resemble kite-type motifs and the structure appears disrupted by atomic scale defects. Atom probe tomography reveals that carbon and boron impurities are co-segregating to the grain boundary reaching levels of >1.5 at%, whereas aluminum is locally depleted by approx. 2 at.%. First-principles calculations indicate that carbon and boron exhibit the strongest segregation tendency and their repulsive interaction with aluminum promotes its depletion from the grain boundary. It is also predicted that substitutional segregation of boron atoms may contribute to local distortions of the kite-type structural units. These results suggest that the co-segregation and interaction of interstitial impurities with substitutional solutes strongly influences grain boundary composition and with this the properties of the interface.

Open and strong-scaling tools for atom-probe crystallography: high-throughput methods for indexing crystal structure and orientation

Volumetric crystal structure indexing and orientation mapping are key data processing steps for virtually any quantitative study of spatial correlations between the local chemical composition features and the microstructure of a material. For electron and X-ray diffraction methods it is possible to develop indexing tools which compare measured and analytically computed patterns to decode the structure and relative orientation within local regions of interest. Consequently, a number of numerically efficient and automated software tools exist to solve the above characterization tasks. For atom-probe tomography (APT) experiments, however, the strategy of making comparisons between measured and analytically computed patterns is less robust because many APT data sets contain substantial noise. Given that sufficiently general predictive models for such noise remain elusive, crystallography tools for APT face several limitations: their robustness to noise is limited, and therefore so too is their capability to identify and distinguish different crystal structures and orientations. In addition, the tools are sequential and demand substantial manual interaction. In combination, this makes robust uncertainty quantification with automated high-throughput studies of the latent crystallographic information a difficult task with APT data. To improve the situation, the existing methods are reviewed and how they link to the methods currently used by the electron and X-ray diffraction communities is discussed. As a result of this, some of the APT methods are modified to yield more robust descriptors of the atomic arrangement. Also reported is how this enables the development of an open-source software tool for strong scaling and automated identification of a crystal structure, and the mapping of crystal orientation in nanocrystalline APT data sets with multiple phases.

Dielectric Properties of Nanoconfined Water: A Canonical Thermopotentiostat Approach

We introduce a novel approach to sample the canonical ensemble at constant temperature and applied electric potential. Our approach can be straightforwardly implemented into any density-functional theory code. Using thermopotentiostat molecular dynamics simulations allows us to compute the dielectric constant of nanoconfined water without any assumptions for the dielectric volume. Compared to the commonly used approach of calculating dielectric properties from polarization fluctuations, our thermopotentiostat technique reduces the required computational time by 2 orders of magnitude.

First-Principles Calculation of the Evaporation Field and Roll-up Effect of M (M = Fe, Cu, Si, and Mn) on the Fe (001) and Fe Step Structure

First-principles calculations were performed on the evaporation field of Fe, Cu, Mn, and Si in Fe (001) and on the evaporation field and roll-up effect of Fe, Cu, and Mn in the Fe (001) step structure. The larger the evaporation barrier energy tendency, at an electric field of 0 V/nm (absorption energy), the larger was the evaporation field. Electric field evaporation calculation results indicate that the order in which the electric field is easily evaporated is Mn > Cu > Fe > Si. The tendency that Mn and Cu evaporate more easily than does Fe and that the evaporation of Si is less probable is consistent with the experiment of a dilute element in steel. In the Fe (001) step structure, when the electric field is low, the roll-up effect where the evaporated atoms move on the step is large, and when the electric field is large, the roll-up effect is small. The roll-up effect of Cu was almost the same as that of Fe, and the roll-up effect of Mn was small because the chemical bond between Mn and Fe was weak.

Atom probe tomography

Atom probe tomography (APT) provides three-dimensional compositional mapping with sub-nanometre resolution. The sensitivity of APT is in the range of parts per million for all elements, including light elements such as hydrogen, carbon or lithium, enabling unique insights into the composition of performance-enhancing or lifetime-limiting microstructural features and making APT ideally suited to complement electron-based or X-ray-based microscopies and spectroscopies. Here, we provide an introductory overview of APT ranging from its inception as an evolution of field ion microscopy to the most recent developments in specimen preparation, including for nanomaterials. We touch on data reconstruction, analysis and various applications, including in the geosciences and the burgeoning biological sciences. We review the underpinnings of APT performance and discuss both strengths and limitations of APT, including how the community can improve on current shortcomings. Finally, we look forwards to true atomic-scale tomography with the ability to measure the isotopic identity and spatial coordinates of every atom in an ever wider range of materials through new specimen preparation routes, novel laser pulsing and detector technologies, and full interoperability with complementary microscopy techniques.

Evaporation and Fragmentation of Organic Molecules in Strong Electric Fields Simulated with DFT

Atom probe tomography allows us to measure the three-dimensional composition of materials with up to atomic resolution by evaporating the material using high electric fields. Initially developed for metals, it is increasingly used for covalently bound structures. To aid the interpretation of the obtained fragmentation pattern, we modeled the fragmentation and desorption of self-assembled monolayers of thiolate molecules on a gold surface in strong electrostatic fields using density functional theory. We used a cluster model and a periodic model of amino-undecanethiolate, NH₂(CH₂)₁₁S, and fluoro-decanethiolate, CF₃(CF₂)₇(CH₂)₂S. In the former molecule, the fragment CH₂NH₂⁺ was found to evaporate at fields of 5.4-7.7 V/nm. It was followed by different hydrocarbon fragments. Fluoro-decanethiolate evaporates CF₃⁺ at fields of 5.7-6.7 V/nm in the cluster model and at 15.4-23.1 V/nm in the periodic model, followed by CF₂⁺ and C₂F₄²⁺. Detailed analysis of the electronic structure during the evaporation process revealed a stepwise accumulation of the charge in the head groups exposed to the strongest fields, followed by dissociation of covalent bonds. These observations will facilitate the analysis of atom probe experiments of covalently bound structures.

Preferential Evaporation in Atom Probe Tomography: An Analytical Approach

Atom probe tomography (APT) analysis conditions play a major role in the composition measurement accuracy. Preferential evaporation (PE), which significantly biases the apparent composition, more than other well-known phenomena in APT, is strongly connected to those analysis conditions. One way to optimize them, in order to have the most accurate measurement, is therefore to be able to predict and then to estimate their influence on the apparent composition. An analytical model is proposed to quantify the PE. This model is applied to three different alloys such as NiCu, FeCrNi, and FeCu. The model explains not only the analysis temperature dependence, as in an already existing model, but also the dependence to the pulse fraction and the pulse frequency. Moreover, the model can also provide an energetic constant directly linked to the energy barrier required to field evaporate atom from the sample surface.