Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “atomic ordering”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 721 records · Page 40
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Artificial neural network potentials for mechanics and fracture dynamics of two-dimensional crystals **

Understanding the mechanics and failure of materials at the nanoscale is critical for their engineering and applications. The accurate atomistic modeling of brittle failure with crack propagation in covalent crystals requires a quantum mechanics-based description of individual bond-breaking events. Artificial neural network potentials (NNPs) have emerged to overcome the traditional, physics-based modeling tradeoff between accuracy and accessible time and length scales. Previous studies have shown successful applications of NNPs for describing the structure and dynamics of molecular systems and amorphous or liquid phases of materials. However, their application to deformation and failure processes in materials is still uncommon. In this study, we discuss the apparent limitations of NNPs for the description of deformation and fracture under loadings and propose a way to generate and select training data for their employment in simulations of deformation and fracture simulations of crystals. We applied the proposed approach to 2D crystalline graphene, utilizing the density-functional tight-binding method for more efficient and extensive data generation in place of density functional theory. Then, we explored how the data selection affects the accuracy of the developed artificial NNPs. It revealed that NNP’s reliability should not only be measured based on the total energy and atomic force comparisons for reference structures but also utilize comparisons for physical properties, e.g. stress–strain curves and geometric deformation. In sharp contrast to popular reactive bond order potentials, our optimized NNP predicts straight crack propagation in graphene along both armchair and zigzag (ZZ) lattice directions, as well as higher fracture toughness of ZZ edge direction. Our study provides significant insight into crack propagation mechanisms on atomic scales and highlights strategies for NNP developments of broader materials.

2D materials↗
DOE OSTI2020

Why Is the Range of Timescale So Wide in Glass-Forming Liquid?

The viscosity and the relaxation time of a glass-forming liquid vary over 15 orders of magnitude before the liquid freezes into a glass. The rate of the change with temperature is characterized by liquid fragility. The mechanism of such a spectacular behavior and the origin of fragility have long been discussed, but it remains unresolved because of the difficulty of carrying out experiments and constructing theories that bridge over a wide timescale from atomic (ps) to bulk (minutes). Through the x-ray diffraction measurement and molecular dynamics simulation for metallic liquids we suggest that large changes in viscosity can be caused by relatively small changes in the structural coherence which characterizes the medium-range order. Here the structural coherence does not imply that of atomic-scale structure, but it relates to the coarse-grained density fluctuations represented by the peaks in the pair-distribution function (PDF) beyond the nearest neighbors. The coherence length is related to fragility and increases with decreasing temperature, and it diverges only at a negative temperature. This analysis is compared with several current theories which predict a phase transition near the glass transition temperature.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Probing electron-hole coherence in strongly driven 2D materials using high-harmonic generation

Under strong electromagnetic excitation, electron–hole (e-h) pairs may be generated in solids, which are subsequently driven to high energy and high momentum, producing high harmonics (HH) of the driving field. The HH efficiency depends on the degree of coherence between the driven electron and hole created by the laser field. Here, we disrupt this e-h coherence in an atomically thin semiconductor by photodoping via incoherent e-h pairs. We observe a strong, systematic harmonic order-dependent intensity reduction. This trend is explained by an exponential decay of the inter-band polarization, proportional to the sub-cycle excursion time of the e-h pair. Our study provides a platform to probe the importance of many-body effects, such as excitation density-dependent decoherence time for strongly driven electrons without the need of ultrashort laser pulses.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

Core-level binding energy shifts between interior, terrace and edge atoms in MnO(001) thin films

Understanding the physical origin of core-level photoemission line shapes can offer valuable information about the chemical and physical properties of surfaces. For instance, in a large number of transition metals oxides, changes in the line shape allow the accurate determination of their oxidation states and cation site symmetry. Yet, despite this importance, experimental investigations on core-level shifts have been much neglected in recent years. In order to provide further evidence of the physical relevance, we have, in this contribution, introduced a new aspect of interior-, terrace- and edge- atom core-level binding energy shifts to describe monolayers of MnO(001) films grown on an Au(111) substrate. By this means we were able to distinguish the line shape contributions related to different types of atomic sites. Here we show that their relative intensities and energy shifts are able to provide information about the relative amount of under coordinated atoms on the surface and, thus give insights into their catalytic properties. Our findings reveal the importance of a detailed surface science characterization to provide the correct interpretation of distinct photoemission line shapes when considering thin film samples as well as nanostructures in general.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Vacancy-mediated complex phase selection in high entropy alloys

Phase selection in Ti-Zr-Hf-Al high-entropy alloys was investigated by in-situ high-energy x-ray diffraction, single-crystal x-ray diffraction, and KKR-CPA electronic-structure methods that address disorder and vacancies, predicting formation enthalpy and chemical short-range order (SRO). Samples with varying Al content were synthesized, characterized, and computationally assessed to ascertain the composition-dependent phase selection, as increased Al content often acts as a stabilizer of a body-centered-cubic structure. Equiatomic TiZrHfAl was especially interesting due to its observed bcc superstructure – a variant of γ-brass with 4 vacancies per cell (not 2 as in γ-brass). We highlight how vacancy ordering mediates selection of this variant of γ-brass, which is driven by vacancy-atom SRO that dramatically suppress all atomic SRO. As vacancies are inherent in processing refractory systems, we expect that similar discoveries await in other high entropy alloys or in revisiting older experimental data.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Equation of state of neutron-rich matter in d -dimensions

Nuclear systems under constraints, with high degrees of symmetries and/or collectivities may be considered as moving effectively in spaces with reduced spatial dimensions. We first derive analytical expressions for the nucleon specific energy E 0 (ρ), pressure P 0 (ρ), incompressibility coefficient K 0 (ρ) and skewness coefficient J 0 (ρ) of symmetric nucleonic matter (SNM), the quadratic symmetry energy E sym (ρ), its slope parameter L(ρ) and curvature coefficient K sym (ρ) as well as the fourth-order symmetry energy E sym,4 (ρ) of neutron-rich matter in general d spatial dimensions (abbreviated as "dD") in terms of the isoscalar and isovector parts of the isospin-dependent single-nucleon potential according to the generalized Hugenholtz-Van Hove (HVH) theorem. The equation of state (EOS) of nuclear matter in dD can be linked to that in the conventional 3-dimensional (3D) space by the ε-expansion which is a perturbative approach successfully used previously in treating second-order phase transitions and related critical phenomena and more recently in studying the EOS of cold atoms. The ε-expansion of nuclear EOS in dD based on a reference dimension d f = d – ε is shown to be effective with –1 ≲ ε ≲ 1 starting from 1 ≲ d f ≲ 3 in comparison with the exact expressions derived using the HVH theorem. Moreover, the EOS of SNM (with/without considering its potential part) is found to be reduced (enhanced) in lower (higher) dimensions, indicating in particular that the many-nucleon system tends to be deeper bounded but saturate at higher densities in spaces with lower dimensions. Here, the links between the EOSs in 3D and dD spaces from the ε-expansion provide new perspectives to the EOS of neutron-rich matter.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Equilibrium carbon isotope fractionation factors of hydrocarbons: Semi-empirical force-field method

Here, we have calculated the reduced partition function ratios for carbon isotopes (β-factor) of 67 hydrocarbons (alkanes, alkenes, alkynes, cycloalkanes, and aromatics), including their 267 single-substituted isotopomers. The calculations were performed using the harmonic oscillator – rigid rotator model and the Urey / Bigeleisen-Mayer approach. Normal frequencies of molecular vibrations of the isotopologues were calculated from the molecular structures data and valence field force constants stored in the Light-handled Elucidation of Vibrations (LEV) database elaborated by Dr. Gribov and his colleagues in Russia. The LEV database is constructed from experimental spectroscopic and structural data by solving inverse problems of molecular vibrations. The LEV is internally consistent and the most comprehensive database available to date. The β-factor were calculated in the temperature range of 200–800 K with a 10 K step. Our calculations predict that the β-factors increase with an increasing number of C atoms within the same groups of hydrocarbons (e.g., alkanes). Our calculations also show a general descending order of 13 C enrichments among the different groups of hydrocarbons: cycloalkanes, aromatics, alkenes (double bonds) and isoalkanes, alkanes, alkynes (triple bonds). Position-specific, intramolecular isotope effects within hydrocarbons are determined by the β-factors of C in the different functional groups in the order: quaternary (C), tertiary (methine, CH), secondary (methylene, CH 2 ), primary (methyl, CH 3 ) and double bond (C=C), saturated bond (C - C), triple bond (HC≡). Our calculations on bulk and position-specific carbon isotope β-factors of the hydrocarbons, which are generally consistent with very limited ab initio calculations in the literature, are internally consistent and the most comprehensive to date for future applications to position-specific isotope geochemistry of hydrocarbons.

58 GEOSCIENCES↗
DOE OSTI2022

Diabatic Valence-Hole States in the C 2 Molecule: “Putting Humpty Dumpty Together Again”

Despite the long history of spectroscopic studies of the C$_2$ molecule, fundamental questions about its chemical bonding are still being hotly debated. The complex electronic structure of C$_2$ is a consequence of its dense manifold of near-degenerate, low-lying electronic states. A global multi-state diabatic model is proposed here to disentangle the numerous configuration interactions within four symmetry manifolds of C$_2$ ($^{1}\Pi_g$, $^{3}\Pi_g$, $^{1}\Sigma_u^+$, and $^{3}\Sigma_u^+$). The key concept of our model is the existence of two "valence-hole" configurations, $2\sigma_g^22\sigma_u^11\pi_{u}^33\sigma_g^2$ for $^{1,3}\Pi_g$ states and $2\sigma_g^22\sigma_u^11\pi_{u}^43\sigma_g^1$ for $^{1,3}\Sigma_u^+$ states that derive from $3\sigma_g\leftarrow2\sigma_u$ electron promotion. The lowest-energy state from each of the four C$_2$ symmetry species is dominated by this type of valence-hole configuration at its equilibrium internuclear separation. Finally, as a result of their large binding energy (nominal bond order of 3) and correlation with the 2s$^2$2p$^2$+2s2p$^3$ separated-atom configurations, the presence of these valence-hole configurations has a profound impact on the $global$ electronic structure and unimolecular dynamics of C$_2$.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Is the Local Ion Density Sufficient to Drive NaCl Nucleation from the Melt and Aqueous Solution?

Even though nucleation is ubiquitous in different science and engineering problems, investigating nucleation is extremely difficult due to the complicated ranges of time and length scales involved. In this work, we simulate NaCl nucleation in both molten and aqueous environments using enhanced sampling of all-atom molecular dynamics with deep-learning-based estimation of reaction coordinates. By incorporating various structural order parameters and learning the reaction coordinate as a function thereof, we achieve significantly improved sampling relative to traditional ad hoc descriptions of what drives nucleation, particularly in an aqueous medium. Our results reveal a one-step nucleation mechanism in both environments, with reaction coordinate analysis highlighting the importance of local ion density in distinguishing solid and liquid states. However, although fluctuations in the local ion density are necessary to drive nucleation, they are not sufficient. Our analysis shows that near the transition states, descriptors such as enthalpy and local structure become crucial. Furthermore, our protocol proposed here enables robust nucleation analysis and phase sampling and could offer insights into nucleation mechanisms for generic small molecules in different environments.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Multiscale Reactive Model for 1,3,5-Triamino-2,4,6-trinitrobenzene Inferred by Reactive MD Simulations and Unsupervised Learning

When high-energy-density materials are subjected to thermal or mechanical insults at extreme conditions (shock loading), a coupled response between the thermo-mechanical and chemical behaviors is systematically induced. Herein we develop a reaction model for the fast chemistry of 1,3,5-triamino-2,4,6-trinitrobenzene (TATB) at the mesoscopic scale, where the chemical behavior is determined by underlying microscopic reactive simulations. The slow carbon cluster formation is not discussed in the present work. All-atom reactive molecular dynamics (MD) simulations are performed with the ReaxFF potential, and a reduced-order chemical kinetics model for TATB is fitted to isothermal and adiabatic simulations of single crystal chemical decomposition. Unsupervised machine learning techniques based on non-negative matrix factorization are applied to MD trajectories to model the decomposition kinetics of TATB in terms of a four-component model. The associated heats of reaction are fit to the temperature evolution from adiabatic decomposition trajectories. Using a chemical species analysis, we show that non-negative matrix factorization captures the main chemical decomposition steps of TATB and provides an accurate estimation of their evolution with temperature. The final analytical formulation, coupled to a diffusion term, is incorporated into a continuum formalism, and simulation results are compared one-to-one against MD simulations of 1D reaction propagation along different crystallographic directions and with different initial temperatures. A good agreement is found for both the temporal and spatial evolution of the temperature field.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Superior Piezo-/Ferro-Electricity in Antiferroelectric Ag x NbO 3-δ Thin Films by Nanopillar Local Structure Design

Antiferroelectrics are fundamental mother compounds critical in developing innovative lead-free piezoelectrics and ferroelectrics and hold great promise for wide-ranging applications in energy conversion and electronic devices. However, harnessing their superior properties presents a significant challenge due to the delicate balance required between their various states. Here, in this study, through the unique design of nanopillar structures to alleviate the local polar heterogeneity, we have achieved significantly improved piezo-/ferro-electricity in classic lead-free antiferroelectric Ag x NbO 3-δ (x = 1, 0.9, and 0.8) epitaxial thin films. The effective piezoelectric coefficient reaches 440 pm V –1 , 1 order of magnitude larger than the stoichiometric AgNbO 3 , rivaling classic lead zirconate titanate piezoelectrics. Atomic-scale electron microscopy investigations unravel the underlying mechanisms. The nanopillars, characterized by antisite occupancy of both Ag and Nb atoms and forming out-of-phase boundaries with the matrix, reduce the local crystal symmetry via interphase strain. This leads to the creation of flexible multinanodomain structures that significantly facilitate polarization rotation, thus substantially enhancing the piezoelectric performance. This study demonstrates the feasibility of engineering local heterogeneity through nanopillar design, offering a generally applicable method for property improvement of a wide range of antiferroelectrics.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Unveiling the complexity of nanodiamond structures

Understanding nanodiamond structures is of great scientific and practical interest. It has been a long-standing challenge to unravel the complexity underlying nanodiamond structures and to resolve the controversies surrounding their polymorphic forms. Here, we use transmission electron microscopy with high-resolution imaging, electron diffraction, multislice simulations, and other supplementary techniques to study the impacts of small sizes and defects on cubic diamond nanostructures. The experimental results show that common cubic diamond nanoparticles display the (200) forbidden reflections in their electron diffraction patterns, which makes them indistinguishable from new diamond (n-diamond). The multislice simulations demonstrate that cubic nanodiamonds smaller than 5 nm can present the d -spacing at 1.78 Å corresponding to the (200) forbidden reflections, and the relative intensity of these reflections increases as the particle size decreases. Our simulation results also reveal that defects, such as surface distortions, internal dislocations, and grain boundaries can also make the (200) forbidden reflections visible. These findings provide valuable insights into the diamond structural complexity at nanoscale, the impact of defects on nanodiamond structures, and the discovery of novel diamond structures.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Magic-angle semimetals with chiral symmetry

Here, we construct and solve a two-dimensional, chirally symmetric model of Dirac cones subjected to a quasiperiodic modulation. In real space, this is realized with a quasiperiodic hopping term. This hopping model, as we show, at the Dirac node energy has a rich phase diagram with a semimetal-to-metal phase transition at intermediate amplitude of the quasiperiodic modulation, and a transition to a phase with a diverging density of states (DOS) and subdiffusive transport when the quasiperiodic hopping is strongest. We further demonstrate that the semimetal-to-metal phase transition can be characterized by the multifractal structure of eigenstates in momentum space and can be considered as a unique “unfreezing” transition. This unfreezing transition in momentum space generates flat bands with a dramatically renormalized bandwidth in the metallic phase similar to the phenomena of the band structure of twisted bilayer graphene at the magic angle. We characterize the nature of this transition numerically as well as analytically in terms of the formation of a band of topological zero modes. For pure quasiperiodic hopping, we provide strong numerical evidence that the low-energy DOS develops a divergence and the eigenstates exhibit Chalker (quantum-critical) scaling despite the model not being random. At particular commensurate limits the model realizes higher-order topological insulating phases. We discuss how these systems can be realized in experiments on ultracold atoms and metamaterials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Electronic structure of α - RuCl 3 by fixed-node and fixed-phase diffusion Monte Carlo methods

Layered material α-RuCl 3 has caught wide attention due to its possible realization of Kitaev's spin liquid and its electronic structure that involves the interplay of electron-electron correlations and spin-orbit effects. Several DFT+U studies have suggested that both electron-electron correlations and spin-orbit effects are crucial for accurately describing the band gap. This work studies the importance of these two effects using fixed-node and fixed-phase diffusion Monte Carlo calculations both in spin-averaged and explicit spin-orbit formalisms. In the latter, the Slater-Jastrow trial function is constructed from two-component spin orbitals using our recent quantum Monte Carlo (QMC) developments and thoroughly tested effective core potentials. Our results show that the gap in the ideal crystal is already accurately described by the spin-averaged case, with the dominant role being played by the magnetic ground state with significant exchange and electron correlation effects. We find qualitative agreement between hybrid DFT, DFT+U, and QMC. In addition, QMC results agree very well with available experiments, and we identify the values of exact Fock exchange mixing that provide comparable gaps. Explicit spin-orbit QMC calculations reveal that the effect of spin-orbit coupling on the gap is minor, of the order of 0.2 eV, which corresponds to the strength of the spin orbit of the Ru atom.

36 MATERIALS SCIENCE↗
DOE OSTI2020

K-shell internal conversion coefficient for M4 decay of the 30.8 keV isomer in 93 Nb

The 30.8 keV (1/2) – isomeric state in $^{93}_{41}$Nb decays by an M4 transition to the (9/2) + ground state. Here, we have measured the K-shell internal conversion coefficient α K for this transition in order to test the validity of current methods of calculation for a case with relatively low atomic number, low transition energy, and a high value of α K . Taking the fluorescence yield of niobium to be 0.751(4), we obtain α K = 2.56(9) × 10 4 , a result that agrees with the Dirac-Fock calculations that take into account the presence of a K vacancy in the final state, and disagrees by 1.7 standard deviations with calculations that ignore the vacancy. We also determine the energy of the isomeric state to be 30.760(5) keV, a fourfold improvement in precision over the currently accepted value.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Learning from metastable symmetric-tilt grain boundaries using physics-based descriptors

Grain boundaries (GBs) govern critical properties of polycrystalline materials. Although significant advancements have been made in characterizing minimum energy and ordered GBs, real GBs are seldom found in such well-defined states. This diversity of atomic arrangements in metastable states makes it challenging to establish structure-property relationships with physical insights. Here, to address this challenge, we use data-driven methods to explore these relationships and examine the underlying physics. In this study, we utilize a large atomistic database (~5000) of minimum energy and metastable states of symmetric-tilt copper GBs, combined with physically motivated local atomic environment (LAE) descriptors [strain functional descriptors (SFDs)], to predict GB properties and gain physical insights. Our regression models exhibit robust predictive capabilities using only 19 descriptors, generalizing to atomic environments in nanocrystals. A significant highlight of our work is the integration of an unsupervised method with SFDs to elucidate LAEs at GBs and their role in determining properties. The model, trained on these minimum energy and metastable GBs using SFDs, predicts the properties of unseen nanocrystals with good accuracy. Our research underscores the role of a physics-based representation of LAEs and the efficacy of data-driven methods in establishing GB structure-property relationships.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Spontaneous exciton dissociation in transition metal dichalcogenide monolayers

Since the seminal work on MoS 2 , photoexcitation in atomically thin transition metal dichalcogenides (TMDCs) has been assumed to result in excitons, with binding energies order of magnitude larger than thermal energy at room temperature. Here, we reexamine this foundational assumption and show that photoexcitation of TMDC monolayers can result in a substantial population of free charges. Performing ultrafast terahertz spectroscopy on large-area, single-crystal TMDC monolayers, we find that up to ~10% of excitons spontaneously dissociate into charge carriers with lifetimes exceeding 0.2 ns. Scanning tunneling microscopy reveals that photocarrier generation is intimately related to mid-gap defects, likely via trap-mediated Auger scattering. Only in state-of-the-art quality monolayers, with mid-gap trap densities as low as 10 9 cm –2 , does intrinsic exciton physics start to dominate the terahertz response. Our findings reveal the necessity of knowing the defect density in understanding photophysics of TMDCs.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Visualizing broken symmetry and topological defects in a quantum Hall ferromagnet

The interaction between electrons in graphene under high magnetic fields drives the formation of a rich set of quantum Hall ferromagnetic (QHFM) phases with broken spin or valley symmetry. Visualizing atomic-scale electronic wave functions with scanning tunneling spectroscopy (STS), we resolved microscopic signatures of valley ordering in QHFM phases and spectral features of fractional quantum Hall phases of graphene. At charge neutrality, we observed a field-tuned continuous quantum phase transition from a valley-polarized state to an intervalley coherent state, with a Kekulé distortion of its electronic density. Mapping the valley texture extracted from STS measurements of the Kekulé phase, we could visualize valley skyrmion excitations localized near charged defects. Our techniques can be applied to examine valley-ordered phases and their topological excitations in a wide range of materials.

valley symmetry↗