Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Atomistic models”

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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Subtle metastability of the layered magnetic topological insulator MnBi 2 Te 4 from weak interactions

Abstract Layered quantum materials can host interesting properties, including magnetic and topological, for which enormous computational predictions have been done. Their thermodynamic stability is much less visited computationally, which however determines the existence of materials and can be used to guide experimental synthesis. MnBi 2 Te 4 is one of such layered quantum materials that was predicted to be an intrinsic antiferromagnetic topological insulator, and later experimentally realized but in a thermodynamically metastable state. Here, using a combined first-principles-based approach that considers lattice, charge, and spin degrees of freedom, we investigate the metastability of MnBi 2 Te 4 by calculating the Helmholtz free energy for the reaction Bi 2 Te 3 + MnTe → MnBi 2 Te 4 . We identify a temperature range (~500–873 K) in which the compound is stable with respect to the competing binary phases, consistent with experimental observation. We validate the predictions by comparing the calculated specific heats contributed from different degrees of freedom with experimental results. Our findings indicate that the degrees of freedom responsible for the van der Waals interaction, lattice vibration, magnetic coupling, and nontrivial band topology in MnBi 2 Te 4 not only enable emergent phenomena but also play a crucial role in determining its thermodynamic stability. This conclusion lays the foundation for the future computational material synthesis of novel layered systems.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Microscopic mechanism of unusual lattice thermal transport in TlInTe 2

We investigate the microscopic mechanism of ultralow lattice thermal conductivity (κ l ) of TlInTe 2 and its weak temperature dependence using a unified theory of lattice heat transport, that considers contributions arising from the particle-like propagation as well as wave-like tunneling of phonons. While we use the Peierls–Boltzmann transport equation (PBTE) to calculate the particle-like contributions (κ l (PBTE)), we explicitly calculate the off-diagonal (OD) components of the heat-flux operator within a first-principles density functional theory framework to determine the contributions (κ l (OD)) arising from the wave-like tunneling of phonons. At each temperature, T, we anharmonically renormalize the phonon frequencies using the self-consistent phonon theory including quartic anharmonicity, and utilize them to calculate κ l (PBTE) and κ l (OD). With the combined inclusion of κ l (PBTE), κ l (OD), and additional grain-boundary scatterings, our calculations successfully reproduce the experimental results. Our analysis shows that large quartic anharmonicity of TlInTe 2 (a) strongly hardens the low-energy phonon branches, (b) diminishes the three-phonon scattering processes at finite T, and (c) recovers the weaker than T –1 decay of the measured κ l .

36 MATERIALS SCIENCE↗
DOE OSTI2021

Performant implementation of the atomic cluster expansion (PACE) and application to copper and silicon

The atomic cluster expansion is a general polynomial expansion of the atomic energy in multi-atom basis functions. Here we implement the atomic cluster expansion in the performant C++ code that is suitable for use in large-scale atomistic simulations. We briefly review the atomic cluster expansion and give detailed expressions for energies and forces as well as efficient algorithms for their evaluation. We demonstrate that the atomic cluster expansion as implemented in shifts a previously established Pareto front for machine learning interatomic potentials toward faster and more accurate calculations. Moreover, general purpose parameterizations are presented for copper and silicon and evaluated in detail. We show that the Cu and Si potentials significantly improve on the best available potentials for highly accurate large-scale atomistic simulations.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Efficient and universal characterization of atomic structures through a topological graph order parameter

Abstract A graph-based order parameter, based on the topology of the graph itself, is introduced for the characterization of atomistic structures. The order parameter is universal to any material/chemical system and is transferable to all structural geometries. Four sets of data are used to validate both the generalizability and accuracy of the algorithm: (1) liquid lithium configurations spanning up to 300 GPa, (2) condensed phases of carbon along with nanotubes and buckyballs at ambient and high temperature, (3) a diverse set of aluminum configurations including surfaces, compressed and expanded lattices, point defects, grain boundaries, liquids, nanoparticles, all at nonzero temperatures, and (4) eleven niobium oxide crystal phases generated with ab initio molecular dynamics. We compare our proposed method to existing, state-of-the-art methods for the cases of aluminum and niobium oxide. Our order parameter uniquely classifies every configuration and outperforms all studied existing methods, opening the door for its use in a multitude of complex application spaces that can require fine structure-level characterization of atomistic graphs.

36 MATERIALS SCIENCE↗
DOE OSTI2022

A flexible and scalable scheme for mixing computed formation energies from different levels of theory

Abstract Computational materials discovery efforts are enabled by large databases of properties derived from high-throughput density functional theory (DFT), which now contain millions of calculations at the generalized gradient approximation (GGA) level of theory. It is now feasible to carry out high-throughput calculations using more accurate methods, such as meta-GGA DFT; however recomputing an entire database with a higher-fidelity method would not effectively leverage the enormous investment of computational resources embodied in existing (GGA) calculations. Instead, we propose here a general procedure by which higher-fidelity, low-coverage calculations (e.g., meta-GGA calculations for selected chemical systems) can be combined with lower-fidelity, high-coverage calculations (e.g., an existing database of GGA calculations) in a robust and scalable manner. We then use legacy PBE(+ U ) GGA calculations and new r 2 SCAN meta-GGA calculations from the Materials Project database to demonstrate that our scheme improves solid and aqueous phase stability predictions, and discuss practical considerations for its implementation.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Ultra-fast interpretable machine-learning potentials

Abstract All-atom dynamics simulations are an indispensable quantitative tool in physics, chemistry, and materials science, but large systems and long simulation times remain challenging due to the trade-off between computational efficiency and predictive accuracy. To address this challenge, we combine effective two- and three-body potentials in a cubic B-spline basis with regularized linear regression to obtain machine-learning potentials that are physically interpretable, sufficiently accurate for applications, as fast as the fastest traditional empirical potentials, and two to four orders of magnitude faster than state-of-the-art machine-learning potentials. For data from empirical potentials, we demonstrate the exact retrieval of the potential. For data from density functional theory, the predicted energies, forces, and derived properties, including phonon spectra, elastic constants, and melting points, closely match those of the reference method. The introduced potentials might contribute towards accurate all-atom dynamics simulations of large atomistic systems over long-time scales.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Circumventing the phonon bottleneck by multiphonon-mediated hot exciton cooling at the nanoscale

Abstract Nonradiative processes govern efficiencies of semiconductor nanocrystal (NC)-based devices. A central process is hot exciton cooling, or the nonradiative relaxation of a highly excited electron/hole pair to form a band-edge exciton. Due to quantum confinement effects, the timescale and mechanism of cooling are not well understood. A mismatch between electronic energy gaps and phonon frequencies has led to the hypothesis of a phonon bottleneck and extremely slow cooling, while enhanced electron-hole interactions have suggested ultrafast cooling. Experimental measurements of the cooling timescale range six orders of magnitude. Here, we develop an atomistic approach to describe phonon-mediated exciton dynamics and simulate cooling in NCs of experimentally relevant sizes. We find that cooling occurs on ~30 fs timescales in CdSe NCs, in agreement with the most recent measurements, and that the phonon bottleneck is circumvented through a cascade of multiphonon-mediated relaxation events. Furthermore, we identify NC handles for tuning the cooling timescale.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Bayesian blacksmithing: discovering thermomechanical properties and deformation mechanisms in high-entropy refractory alloys

Finding alloys with specific design properties is challenging due to the large number of possible compositions and the complex interactions between elements. This study introduces a multi-objective Bayesian optimization approach guiding molecular dynamics simulations for discovering high-performance refractory alloys with both targeted intrinsic static thermomechanical properties and also deformation mechanisms occurring during dynamic loading. The objective functions are aiming for excellent thermomechanical stability via a high bulk modulus, a low thermal expansion, a high heat capacity, and for a resilient deformation mechanism maximizing the retention of the BCC phase after shock loading. Contrasting two optimization procedures, we show that the Pareto-optimal solutions are confined to a small performance space when the property objectives display a cooperative relationship. Conversely, the Pareto front is much broader in the performance space when these properties have antagonistic relationships. Density functional theory simulations validate these findings and unveil underlying atomic-bond changes driving property improvements.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Feature engineering descriptors, transforms, and machine learning for grain boundaries and variable-sized atom clusters

Abstract Obtaining microscopic structure-property relationships for grain boundaries is challenging due to their complex atomic structures. Recent efforts use machine learning to derive these relationships, but the way the atomic grain boundary structure is represented can have a significant impact on the predictions. Key steps for property prediction common to grain boundaries and other variable-sized atom clustered structures include: (1) describing the atomic structure as a feature matrix, (2) transforming the variable-sized feature matrix to a fixed length common to all structures, and (3) applying a machine learning algorithm to predict properties from the transformed matrices. We examine how these steps and different combinations of engineered features impact the accuracy of grain boundary energy predictions using a database of over 7000 grain boundaries. Additionally, we assess how different engineered features support interpretability, offering insights into the physics of the structure-property relationships.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Anisotropic epitaxial stabilization of a low-symmetry ferroelectric with enhanced electromechanical response

Here piezoelectrics interconvert mechanical energy and electric charge and are widely used in actuators and sensors. The best performing materials are ferroelectrics at a morphotropic phase boundary, where several phases coexist. Switching between these phases by electric field produces a large electromechanical response. In ferroelectric BiFeO3, strain can create a morphotropic-phase-boundary-like phase mixture and thus generate large electric-field-dependent strains. However, this enhanced response occurs at localized, randomly positioned regions of the film. Here, we use epitaxial strain and orientation engineering in tandem-anisotropic epitaxy-to craft a low-symmetry phase of BiFeO3 that acts as a structural bridge between the rhombohedral-like and tetragonal-like polymorphs. Interferometric displacement sensor measurements reveal that this phase has an enhanced piezoelectric coefficient of x2.4 compared with typical rhombohedral-like BiFeO3. Band-excitation frequency response measurements and first-principles calculations provide evidence that this phase undergoes a transition to the tetragonal-like polymorph under electric field, generating an enhanced piezoelectric response throughout the film and associated field-induced reversible strains. These results offer a route to engineer thin-film piezoelectrics with improved functionalities, with broader perspectives for other functional oxides.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Electrochemically induced amorphous-to-rock-salt phase transformation in niobium oxide electrode for Li-ion batteries

Intercalation-type metal oxides are promising negative electrode materials for safe rechargeable lithium-ion batteries due to the reduced risk of Li plating at low voltages. Nevertheless, their lower energy and power density along with cycling instability remain bottlenecks for their implementation, especially for fast-charging applications. Here, we report a nanostructured rock-salt Nb 2 O 5 electrode formed through an amorphous-to-crystalline transformation during repeated electrochemical cycling with Li + . This electrode can reversibly cycle three lithiums per Nb 2 O 5 , corresponding to a capacity of 269 mAh g –1 at 20 mA g –1 , and retains a capacity of 191 mAh g –1 at a high rate of 1 A g –1 . Further, it exhibits superb cycling stability with a capacity of 225 mAh g –1 at 200 mA g –1 for 400 cycles, and a Coulombic efficiency of 99.93%. We attribute the enhanced performance to the cubic rock-salt framework, which promotes low-energy migration paths. Our work suggests that inducing crystallization of amorphous nanomaterials through electrochemical cycling is a promising avenue for creating unconventional high-performance metal oxide electrode materials.

25 ENERGY STORAGE↗
DOE OSTI2024

Substantial oxygen loss and chemical expansion in lithium-rich layered oxides at moderate delithiation

Delithiation of layered oxide electrodes triggers irreversible oxygen loss, one of the primary degradation modes in lithium-ion batteries. However, the delithiation-dependent mechanisms of oxygen loss remain poorly understood. Here we investigate the oxygen non-stoichiometry in Li 1.18–x Ni 0.21 Mn 0.53 Co 0.08 O 2–δ electrodes as a function of Li content by using cycling protocols with long open-circuit voltage steps at varying states of charge. Surprisingly, we observe substantial oxygen loss even at moderate delithiation, corresponding to 2.5, 4.0 and 7.6 ml O 2 per gram of Li 1.18–x Ni 0.21 Mn 0.53 Co 0.08 O 2–δ after resting at upper capacity cut-offs of 135, 200 and 265 mAh g −1 for 100 h. Our observations suggest an intrinsic oxygen instability consistent with predictions of high oxygen activity at intermediate potentials versus Li/Li + . In addition, we observe a large chemical expansion coefficient with respect to oxygen non-stoichiometry, which is about three times greater than those of classical oxygen-deficient materials such as fluorite and perovskite oxides. Furthermore, our work challenges the conventional wisdom that deep delithiation is a necessary condition for oxygen loss in layered oxide electrodes and highlights the importance of calendar ageing for investigating oxygen stability.

Atomistic models↗
DOE OSTI2020

Observation of single-defect memristor in an MoS2 atomic sheet

Non-volatile resistive switching, also known as memristor effect, where an electric field switches the resistance states of a two-terminal device, has emerged as an important concept in the development of high-density information storage, computing and reconfigurable systems. The past decade has witnessed substantial advances in non-volatile resistive switching materials such as metal oxides and solid electrolytes. It was long believed that leakage currents would prevent the observation of this phenomenon for nanometre-thin insulating layers. However, the recent discovery of non-volatile resistive switching in two-dimensional monolayers of transition metal dichalcogenide and hexagonal boron nitride sandwich structures (also known as atomristors) has refuted this belief and added a new materials dimension owing to the benefits of size scaling. Here we elucidate the origin of the switching mechanism in atomic sheets using monolayer MoS 2 as a model system. Atomistic imaging and spectroscopy reveal that metal substitution into a sulfur vacancy results in a non-volatile change in the resistance, which is corroborated by computational studies of defect structures and electronic states. These findings provide an atomistic understanding of non-volatile switching and open a new direction in precision defect engineering, down to a single defect, towards achieving the smallest memristor for applications in ultra-dense memory, neuromorphic computing and radio-frequency communication systems.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2020

Observations of grain-boundary phase transformations in an elemental metal

The theory of grain boundary (the interface between crystallites, GB) structure has a long history 1 and the concept of GBs undergoing phase transformations was proposed 50 years ago 2,3 . The underlying assumption was that multiple stable and metastable states exist for different GB orientations 4,5,6 . The terminology ‘complexion’ was recently proposed to distinguish between interfacial states that differ in any equilibrium thermodynamic property 7 . Different types of complexion and transitions between complexions have been characterized, mostly in binary or multicomponent systems 8,9,10,11,12,13,14,15,16,17,18,19 . Simulations have provided insight into the phase behaviour of interfaces and shown that GB transitions can occur in many material systems 20,21,22,23,24 . However, the direct experimental observation and transformation kinetics of GBs in an elemental metal have remained elusive. In this work we demonstrate atomic-scale GB phase coexistence and transformations at symmetric and asymmetric $[11\bar{1}]$ tilt GBs in elemental copper. Atomic-resolution imaging reveals the coexistence of two different structures at Σ19b GBs (where Σ19 is the density of coincident sites and b is a GB variant), in agreement with evolutionary GB structure search and clustering analysis 21,25,26 . We also use finite-temperature molecular dynamics simulations to explore the coexistence and transformation kinetics of these GB phases. Our results demonstrate how GB phases can be kinetically trapped, enabling atomic-scale room-temperature observations. Our work paves the way for atomic-scale in situ studies of metallic GB phase transformations, which were previously detected only indirectly 9,15,27,28,29 , through their influence on abnormal grain growth, non-Arrhenius-type diffusion or liquid metal embrittlement.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Diversity-oriented synthesis of polymer membranes with ion solvation cages

Microporous polymers feature shape-persistent free volume elements (FVEs), which are permeated by small molecules and ions when used as membranes for chemical separations, water purification, fuel cells, and batteries. It remains a significant challenge to identify FVEs with analyte specificity, due to difficulties in generating microporous polymer libraries with sufficient diversity for screening their properties. Here, we describe a diversity-oriented synthetic (DOS) strategy for microporous polymer membranes from which we identified those whose FVEs serve as solid solvation cages for lithium ions (Li + ). Furthermore, key elements of our strategy included diversification of bis(catechol)-type monomers via multi-component Mannich reactions to introduce Li+-coordinating functionality within individual FVEs, topology-enforcing polymerizations for generating macromolecular skeletal diversity for networking FVEs into different pore architectures, and several classes of on-polymer reactions for diversifying pore geometries and dielectric properties. Lead candidate polymer membranes featuring explicit ion solvation cages exhibited both higher ionic conductivity and higher cation transference number than control membranes where FVEs were aspecific, which indicates conventional bounds for membrane permeability and selectivity for ion transport can be overcome.4 These advantages are tied to enhanced Li + partitioning from the electrolyte when the cages are present, higher diffusion barriers for anions within the pores, and network-enforced restrictions on the number of solvent molecules bound to Li+ by comparison to the bulk electrolyte, which reduces the effective mass of the working ion. Such membranes show promise as anode-stabilizing interlayers in high-voltage lithium-metal batteries for electric mobility.

42 ENGINEERING↗
DOE OSTI2024

The persistence of memory in ionic conduction probed by nonlinear optics

Predicting practical rates of transport in condensed phases enables the rational design of materials, devices and processes. This is especially critical to developing low-carbon energy technologies such as rechargeable batteries. For ionic conduction, the collective mechanisms, variation of conductivity with timescales and confinement, and ambiguity in the phononic origin of translation, call for a direct probe of the fundamental steps of ionic diffusion: ion hops. However, such hops are rare-event large-amplitude translations, and are challenging to excite and detect. Here we use single-cycle terahertz pumps to impulsively trigger ionic hopping in battery solid electrolytes. This is visualized by an induced transient birefringence, enabling direct probing of anisotropy in ionic hopping on the picosecond timescale. The relaxation of the transient signal measures the decay of orientational memory, and the production of entropy in diffusion. We extend experimental results using in silico transient birefringence to identify vibrational attempt frequencies for ion hopping. Using nonlinear optical methods, we probe ion transport at its fastest limit, distinguish correlated conduction mechanisms from a true random walk at the atomic scale, and demonstrate the connection between activated transport and the thermodynamics of information.

25 ENERGY STORAGE↗
DOE OSTI2022

A representation-independent electronic charge density database for crystalline materials

Abstract In addition to being the core quantity in density-functional theory, the charge density can be used in many tertiary analyses in materials sciences from bonding to assigning charge to specific atoms. The charge density is data-rich since it contains information about all the electrons in the system. With the increasing prevalence of machine-learning tools in materials sciences, a data-rich object like the charge density can be utilized in a wide range of applications. The database presented here provides a modern and user-friendly interface for a large and continuously updated collection of charge densities as part of the Materials Project. In addition to the charge density data, we provide the theory and code for changing the representation of the charge density which should enable more advanced machine-learning studies for the broader community.

36 MATERIALS SCIENCE↗