Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “strain energy density”

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 523 records · Page 29
All sourcesNASA NTRSDOE OSTI
NASA NTRS2003

Water Impact Test and Simulation of a Composite Energy Absorbing Fuselage Section

In March 2002, a 25-ft/s vertical drop test of a composite fuselage section was conducted onto water. The purpose of the test was to obtain experimental data characterizing the structural response of the fuselage section during water impact for comparison with two previous drop tests that were performed onto a rigid surface and soft soil. For the drop test, the fuselage section was configured with ten 100-lb. lead masses, five per side, that were attached to seat rails mounted to the floor. The fuselage section was raised to a height of 10-ft. and dropped vertically into a 15-ft. diameter pool filled to a depth of 3.5-ft. with water. Approximately 70 channels of data were collected during the drop test at a 10-kHz sampling rate. The test data were used to validate crash simulations of the water impact that were developed using the nonlinear, explicit transient dynamic codes, MSC.Dytran and LS-DYNA. The fuselage structure was modeled using shell and solid elements with a Lagrangian mesh, and the water was modeled with both Eulerian and Lagrangian techniques. The fluid-structure interactions were executed using the fast general coupling in MSC.Dytran and the Arbitrary Lagrange-Euler (ALE) coupling in LS-DYNA. Additionally, the smooth particle hydrodynamics (SPH) meshless Lagrangian technique was used in LS-DYNA to represent the fluid. The simulation results were correlated with the test data to validate the modeling approach. Additional simulation studies were performed to determine how changes in mesh density, mesh uniformity, fluid viscosity, and failure strain influence the test-analysis correlation.

Fasanella, Edwin L.↗
DOE OSTI2022

Ab initio calculations of third-order elastic coefficients

Third-order elasticity (TOE) theory predicts strain-induced changes in second-order elastic coefficients (SOECs) and can model elastic wave propagation in stressed media. Although third-order elastic tensors have been determined based on first principles in previous studies, their current definition is based on an expansion of thermodynamic energy in terms of the Lagrangian strain near the natural, or zero pressure, reference state. This definition is inconvenient for predictions of SOECs under significant initial stresses. Therefore, when TOE theory is necessary to study the strain dependence of elasticity, the seismological community has resorted to an empirical version of the theory. This study reviews the thermodynamic definition of the third-order elastic tensor and proposes using an “effective” third-order elastic tensor. An explicit expression for the effective third-order elastic tensor is given and verified. Additionally, we extend the ab initio approach to calculate third-order elastic tensors under finite pressure and apply it to two cubic systems, namely, NaCl and MgO. As applications and validations, we evaluate (a) strain-induced changes in SOECs and (b) pressure derivatives of SOECs based on ab initio calculations. Good agreement between third-order elasticity-based predictions and numerically calculated values confirms the validity of our theory.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Informed unsupervised machine learning analysis of dislocation microstructure from high-resolution differential aperture X-ray structural microscopy data

This study leverages high-resolution differential-aperture X-ray structural microscopy (DAXM) to probe the local dislocation structure in deformed 304L-stainless steel at small strain, by measuring the lattice rotation and deviatoric elastic strain with a sub-micron resolution. For a single grain in a polycrystalline specimen, the measured lattice rotation field over the measured volume exhibited a multimodal distribution while the deviatoric elastic strain showed a single-mode distribution. An unsupervised Cauchy mixture machine learning model was developed to resolve the multimodal distribution of the lattice rotation. By mapping the lattice rotation data associated with each Cauchy peak in the model back onto the measured volume, we identify contiguous regions of the crystal rotated near the average values corresponding to the peaks of the overall rotation distribution. These regions represent the grain subdivision in the microstructure. Finally, the dislocation density tensor was also computed and its norm was laid over the rotation field to detect the subgrain boundaries. This step provided a validation of the Cauchy mixture model for the analysis of the lattice rotation distribution. The current study highlights the integration of advanced X-ray microscopy techniques with data-driven analysis methods to uncover detailed microstructure scales in deformed crystals.

Machine learning; Lattice rotation; High-energy X-↗
DOE OSTI2022

Nanoscale Control of Polar Surface Phases in Layered van der Waals CuInP 2 S 6

Antiferroelectric (AFE) materials, in which alternating dipole moments cancel out to a zero net macroscopic polarization, can be used for high-density energy storage and memory applications. The AFE phase can exist in bulk CuInP 2 Se 6 , CuBiP 2 S 6 , and a few other transition-metal thiophosphates below 200 K. The required low temperature poses challenges for practical applications. In this work, we report the coexistence of ferrielectric (FE) states and a stable surface phase that does not show piezoelectric response (“zero-response phase”) in bulk CuInP 2 S 6 at room temperature. Using piezoresponse force microscopy (PFM) tomographic imaging together with density functional theory, we find that direct and alternating voltages can locally and stably convert FE to zero-response phases and vice versa. While PFM loops show pinched hystereses reminiscent of antiferroelectricity, PFM tomography reveals that the zero-response areas form only on top of the FE phase in which the polarization vector is pointing up. Theoretical calculations suggest that the zero-response phase may correspond to AFE ordering where stacked CuInP 2 S 6 layers have alternating polarization orientations thereby leading to a net-zero polarization. Consistent with experimental findings, theory predicts that the FE polarization pointing down is robust up to the top surface, whereas FE polarization pointing up energetically favors the formation of an AFE surface layer, whose thickness is likely to be sensitive to local strains. AFE order is likely to be more robust against detrimental size effects than polar order, therefore providing additional opportunities to create multifunctional heterostructures with 2D electronic materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2003

Foam on Tile Impact Modeling for the Space Shuttle Program

Following the breakup of the Space Shuttle Columbia during reentry a NASA-wide investigation team was formed to examine the probable damage inflicted on Orbiter Thermal Protection System (TPS) elements by impact of External Tank insulating foam projectiles. Our team was to apply rigorous, physics-based analysis techniques to help determine parameters of interest for an experimental test program, utilize validated codes to investigate the full range of impact scenarios, and use analysis derived models to predict aero-thermal-structural responses to entry conditions. We were to operate on a non-interference basis with the j Team, and were to supply significant findings to that team and to the Orbiter Vehicle Engineering Working Group, being responsive to any solicitations for support from these entities. The authors formed a working sub-group within the larger team to apply the Smooth Particle Hydrodynamics code SPHC to the damage estimation problem. Numerical models of the LI-900 TPS tiles and of the BX-250 foam were constructed and used as inputs into the code. Material properties needed to properly model the tiles and foam were obtained from other working sub-groups who performed tests on these items for this purpose. Two- and three- dimensional models of the tiles were constructed, including the glass outer layer, the densified lower layer of LI-900 insulation, the Nomex felt Strain Isolation Pad (SIP) mounting layer, and the underlying aluminum 2024 vehicle skin. A model for the BX-250 foam including porous compression, elastic rebound, and surface erosion was developed. Code results for the tile damage and foam behavior were extensively validated through comparison with the Southwest Research Institute (SwRI) foam-on-tile impact experiments carried out in 1999. These tests involved small projectiles striking individual tiles and small tile arrays. Following code and model validation we simulated impacts of larger ET foam projectiles on the TPS tile systems used on the wings of the orbiter. Tiles used on the Wing Acreage, the Main Landing Gear Door, and the Carrier Panels near the front edge of the wing were modeled. Foam impacts shot for the CAB investigation were modeled, as well as impacts at larger angles, including rapid rotation of the projectile, and with varying foam properties. General results suggest that foam impacts on tiles at about 500 mph could cause appreciable damage if the impact angle is greater than about 20 degrees. Some variations of the foam properties, such as increased brittleness or increased density could increase damage in some cases. Rapid (17 rps) rotation failed to increase the damage for the two cases considered. This does not rule out other cases in which the rotational energy might lead to an increase in tile damage, but suggests that in most cases rotation will not be an important factor. Similar models will be applied for other impacting materials, other velocities, and other geometries as part of the Return to Flight process.

Stellingwerf, R. F.↗
DOE OSTI2020

Initial development of an RIA envelope for dispersed nuclear fuel

A reactivity insertion accident (RIA) is a design basis accident in which reactivity is rapidly injected as a result of a control rod ejection or blade drop scenario. The resultant power increase can result in fuel rod failure and subsequent release of radioactive material into the reactor’s primary system. For light water reactors, failure under RIAs is typically a result of fuel melt, pellet-cladding interaction, or rod over-pressurization. Fuel melt occurs when the fuel system is unable to transport heat out of the fuel system. Pellet-cladding interaction occurs when an aggressive fuel pellet expansion pushes on an embrittled cladding beyond its strain limit. Rod over-pressurization occurs when a fuel rod exceeds the critical heat flux and departs from nucleate boiling, causing the rod internal pressure to rapidly increase, exceeding the systems pressure, and balloon until it bursts. To prevent failure, regulators have developed safety requirements that limit the injected enthalpy based on the state of the fuel system. However, dispersed fuel could enable the removal or relaxation of regulator-imposed safety criteria. Dispersed fuel embeds fuel particles in a highly conductive metal. Dispersed particles can have a particle radius of 1 µm to greater than 100 µm. Traditional fuel systems are plagued by poor thermal conductivity, especially at higher burnups, thus reducing the ability of the fuel system to respond to an RIA. However, dispersed fuel could improve the ability of the fuel system to transport injected energy away from the fissile material and into the coolant more efficiently. Furthermore, dispersing the fuel particles mitigates hard contact that typically occurs in the traditional zirconium-uranium dioxide fuel system, thus removing pellet-cladding mechanical interaction as a potential failure mechanism. This paper evaluates fuel particles dispersed in a metal matrix using the BISON fuel performance code to develop an initial failure threshold based on melt temperature of the fuel particle and/or cladding material as a function of fission density for beginning-of-life and end-of-life conditions. BISON results show that (1) reducing the size of the particles allows fission density to increase as a function of time, (2) particle-to-particle proximity must be considered to evaluate the limiting conditions, and (3) improving the fuel particle thermal conductivity improves thermal performance.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Solute effects upon dislocation motion and recovery in Mg alloys (Final Report)

The objective the research was to develop a firmer understanding of the interactions between substitutional solute atoms and dislocations within Mg alloys. These interactions govern the absolute and relative mobilities of various types of dislocations in Mg (e.g., basal < a>, non-basal < a>, and pyramidal < c+a>). Furthermore, they have an impact on dislocation recovery processes (e.g., cross-glide, climb, rearrangement, and annihilation). Ultimately, solute-dislocation interactions strongly impact a) strain hardening, b) strain rate sensitivity, c) plastic anisotropy of textured polycrystals, d) texture evolution, e) dislocation substructures evolution, and even f) recrystallization behavior. It is broadly known that Mg alloys exhibit poor low temperature formability. Since formability is largely governed by properties (a) – (c) in this list, it is critical to develop a better understanding of solute-dislocation interactions, if one hopes to improve the situation. While the theory of static solute strengthening is well developed, especially for alloys with face centered cubic (FCC) crystal structures, there are outstanding questions related to applications to hexagonal close packed (HCP) crystal structures and dynamic strain aging (DSA) in materials of various structures. DSA has far-reaching implications for metal formability and, in the case of Mg alloys, appears to correlate with the so-called rare-earth (RE) texture which has been shown to benefit formability. The Portevin-Le Chatelier (PLC) effect, and associated negative strain rate sensitivity, occur at higher temperatures (>100℃) in Mg alloys as compared with similar Al alloys, even though they have similar melting points and solute diffusivities. Our preliminary research has shown that modern, physics-based models of DSA can be tuned to describe the behavior of Mg alloys if a rather higher activation enthalpy is assumed for cross-core diffusivity. While this partially explains the delay in DSA to higher temperatures, it is also hypothesized that this delay is due in part to the intrinsically more thermally activated (rate sensitive) nature of non-basal < a> dislocation motion which is required for macroscopic flow of Mg alloys, whereas octahedral slip in many FCC metals like aluminum is essentially athermal at room temperature. It was originally proposed to employ a combination of in-situ diffraction-based experimental characterization to validate existing theory. It was envisioned to perform in-situ transmission electron microscopy (TEM) to assess individual dislocation behavior and in-situ high-energy X-ray diffraction (HEXRD) techniques which were showing great promise for elucidating collective dislocation behavior, 2 including recovery, especially if the contributions to various forms of diffraction peak broadening (𝜂𝜂,𝜔𝜔,and 2𝜃𝜃) can be effectively integrated. Finally, it was envisioned to perform discrete dislocation dynamics (DDD) modeling approaches to aide in the interpretation of both TEM and HEXRD experiments. In the end, mechanical tests were performed on more complex Mg alloys which exhibited evidence of dynamic strain aging, and this led to the establishment of another project. Mechanical test data obtained at McMaster University served as the basis of an assessment of the applicability of Bazinski’s “stress equivalence” theory of solute strengthening to polycrystalline alloys of Mg. Although we did not succeed in applying the approach to Mg alloys, we did develop expertise with the HEXRD approach using a BCC, β-Ti alloy and demonstrated numerous new capabilities that may be applied to any polycrystalline material in collaboration with researchers at CHESS and around the world: (1) assessment of details of the elastoplastic transition (yielding) using a combination of HEXRD and full-field polycrystal plasticity modeling, (2) the first-ever experimental observation of strong stress rotation within the individual grains of a polycrystalline material, and (3) a comprehensive analysis of the grain-level dislocation density evolution based upon diffraction peak broadening along 𝜂𝜂,𝜔𝜔,and 2𝜃𝜃 directions. This final aspect allowed us to confirm that dislocations were gliding on multiple plane types and not restricted to {110} type planes, and it also provided clues as to why some grains were unloading during straining, with surprising implications for our understanding of the effects of geometrically necessary dislocations (GNDs). Finally, graduate student, Mohammed Shabana, developed a MATLAB code which confirmed the conclusions of Prof. Catalin Picu (Rensallear Polytechnic Institute, RPI) regarding the effect of solute-trapped, forest dislocations on the breaking stress of Lomer lock junctions in FCC metal alloys. He applied the same anisotropic line-tension model to a variety of dislocation junction configurations and found an inconsistency in the widely cited results of Dupuy and Fivel regarding the Hirth Lock, and he outlined an approach to extend these finding to HCP Mg alloys that we are still pursuing with discretionary fundings at UVA.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Effects of chemical segregation on ductility-anisotropy in high strength Fe-Mn-Al-C lightweight austenitic steels

Ductile Fe-Mn-Al-C lightweight steels offer great potential as high specific strength materials for weight critical applications. In the present work, very high yield (~1 GPa) and ultimate tensile strength levels (~1.5 GPa) were obtained along the rolling and transverse directions of Fe-30Mn-8.5Al-0.9Si-0.9C-0.5Mo (wt.%), with high tensile elongation to failure of ~35%. This excellent mechanical behavior was attributed to the nanoscale, ordered κ-carbide precipitates formed during aging after rolling. However, tensile ductility was not isotropic through all directions of the rolled plate. In particular, the ductility was consistently lower than 5% true strain along the normal direction. Upon detailed microstructural investigations, micro-scale lamellar chemical micro-segregation bands were detected, which caused the observed anisotropic embrittlement along certain plate directions. These bands resulted in low failure strains along the normal direction since it is oriented perpendicular to these chemical segregation bands, resulting in a failure type reminiscent of delamination-induced cleavage fracture. Chemical segregation causes an inhomogeneous κ-carbide precipitate distribution upon aging, producing more pronounced mechanical anisotropy after peak-aging as compared to the solution heat treated samples. Surprisingly, deformation twinning was observed after tensile deformation of solution-heat treated specimens via electron-backscatter diffraction experiments, which is anomalous for a high stacking faulty energy (SFE) material. Further, it was proposed that chemical micro-segregation promotes locally different SFE regions, leading to deformation twinning in this nominally high SFE material. Overall, this work demonstrates the effects of chemical segregation on mechanical properties and microstructural evolution in a high performance, lightweight Fe-Mn-Al-C steel.

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

A phase-field model for non-isothermal phase transformation and plasticity in polycrystalline yttria-stabilized tetragonal zirconia

Here we propose an elastoplastic phase-field (PF) model to investigate the mechanics of tetragonal-to-monoclinic phase transformation (TMPT) and elastoplastic deformation of polycrystalline yttria-stabilized tetragonal zirconia (YSTZ). A Landau polynomial with non-vanishing chemical energy at the equilibrium temperature is introduced to account for the actual formation energies of the phases. The effects of different grain orientations, latent heat, and temperature on TMPT and deformation mechanisms are considered. The suppressive transformation effects of the grain boundaries (GBs) is modeled using an inhomogeneous kinetic coefficient in the bulk and GBs. The simulation results for single crystals demonstrate the capability of the model to reproduce the orientation-dependent compressive deformation of YSTZ similar to atomistic simulations and micropillar experiments. The single crystal with [100] crystallographic orientation along the loading direction (SC[100]) displays both TMPT and plasticity, SC[101] experiences only phase transformation, while SC[001] undergoes only plastic yielding. The TMPT induced by compressive loading exhibits shape memory effect (SME) below the equilibrium transformation temperature and pseudoelasticity above it, while the critical transformation stress increases with increasing loading temperature. The irrecoverable plastic strain is found to trap a part of the monoclinic phase, which prevents a complete reverse transformation. The polycrystalline cases also display SME and PE at low and high temperatures, respectively. Due to the orientation differences between grains and the stress concentrations at geometric nonlinearities, plastic deformation occurs in polycrystalline YSTZ for an applied load less than the yield stress. The results suggest a possible limitation of plasticity and an improvement of the shape recovery of YSTZ if one can control the orientation of the grains and/or increase the density of stacking faults at the GBs during material processing.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Factors Affecting the Electron Conductivity in Single Crystal Li 7 La 3 Zr 2 O 12 and Li 7 P 3 S 11

One of the serious challenges in all solid-state Li ion batteries is neutral Li intrusion into the solid-state electrolyte that can ultimately cause catastrophic failure. One possibility for this is due to n-type electron conductivity that induces the reaction Li + + e – → Li 0 at sites where the potential is less than the Li + /Li potential. This paper reports hybrid density functional theory calculations of the electronic conductivity in two prototype single crystalline solid-state electrolytes, cubic Li 7 La 3 Zr 2 O 12 (c-LLZO) and Li 7 P 3 S 11 (LPS). The formation energies of important point defects that can affect electron conductivity are determined, and we find that the mechanism of n-type electron conductivity for both solid-state electrolytes is via “small” electron polaron hopping, where the quotes signify that substantial Li ion rearrangement is associated with the polaron formation and its migration. In both electrolytes, the formation energies for the small polarons at the Fermi energy are too high to generate measurable electron conductivity at room temperature. For c-LLZO, the concentration of electron polarons necessary to ensure charge neutrality from positively charged oxygen vacancies formed in synthesis can be significantly higher. Hence, the electron conductivity could be significant when measured with ion-blocking metal electrodes, and we discuss how the synthesis conditions could affect this magnitude. However, in the solid-state battery, these polarons are replaced by negatively charged Li vacancies so that the electron conductivity should remain minimal. For LPS single crystals, the inherent minimal electron conductivity is independent of synthesis conditions. Here we also show that the cost of forming Li 0 in bulk c-LLZO is enormous due to strain effects so that it could only potentially form at voids, grain boundaries, or around vacancy defects which relax the lattice strain.

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

Laser Powder Directed Energy Deposition of Steels for Nuclear Applications

This comprehensive investigation examines the structure–property relationships in two nuclear alloy systems—Alloy 709 (A709) austenitic stainless steel and Grade 92 (G-92) ferritic/martensitic (F/M) steel—manufactured via directed energy deposition (DED) for sodium-cooled fast reactor applications. This study establishes the fundamental mechanisms for controlling microstructures for optimizing the mechanical performance of additively manufactured nuclear materials through systematic heat treatment optimization and multiscale characterization. As-deposited A709 steel develops a complex multiscale strengthening architecture consisting of a fine cellular solidification structure with diameter of 2-3 µm within10–50 µm grains, elevated dislocation densities from rapid thermal cycling, and grain boundary precipitates that activate concurrent Hall–Petch, dislocation, and precipitation hardening mechanisms to achieve exceptional properties [yield strength (YS): 603 MPa, ultimate tensile strength (UTS): 844 MPa, Vickers hardness: 220 HV] that achieve a 44% superior strength compared to that of the wrought material. Heat treatments produce different results. Solution annealing (SA) dissolves the cellular structure and reduces the hardness to 190 HV. Precipitation treatment (PT) keeps the cellular structure but adds carbides, allowing the hardness to reach 205 HV. The best approach combines both treatments (SA+PT) and creates uniform precipitate distributions with M 23 C 6 carbides at the grain boundaries and MX carbonitrides in the matrix, achieving a hardness of 195 HV. However, directional differences persist, with a 12%–15% strength variation between orientations due to the inherited layered microstructural architecture that survives aggressive heat treatment. While tensile testing at 550°C demonstrates 40%–50% thermal softening with dynamic strain aging, DED A709 steel still maintains a 71% higher YS than that of the wrought material. Ion irradiation studies (100–400 dpa) of DED A709 steel reveal progressive radiation damage with increasing void density and radiation-induced segregation causing nickel enrichment and chromium depletion, which will ultimately compromise mechanical properties. As-deposited G-92 exhibits exceptional strength (UTS: 1650–1700 MPa, 430 HV) through a complex microstructure containing both ferrite and martensite phases, a high geometrically necessary dislocation (GND) density (17.04×10 14 /m 2 ), and fine carbides. Heat treatments create distinct changes. Normalizing produces fresh martensite with the highest hardness (460 HV) and an increased GND density (20.23×10 14 /m 2 ). Tempering develops dual precipitation systems and reduces the hardness to 290 HV. The optimal approach uses sequential normalizing plus tempering, achieving balanced properties with the lowest hardness (250 HV) and a reduced GND density (11.01×10 14 /m 2 ). A processing-dependent anisotropy is observed: horizontal specimens achieve superior ductile behavior, while vertical specimens exhibit brittle failure. A tempering heat treatment successfully mitigates this anisotropic behavior by transforming the hard martensitic as-deposited structure into tempered martensite enabling both horizontal and vertical specimens to exhibit similar stress–strain characteristics with visible necking behavior. Remarkably, testing at 550°C reveals a reversal in the anisotropy, where as-deposited specimens achieve near isotropy with superior thermal stability (a 15%–20% strength reduction), while tempered specimens develop an orientation dependence with a 25%–30% strength reduction. Both alloy systems demonstrate that DED processing creates specimens with a superior strength through refined microstructural features, though with distinct strengthening mechanisms—austenitic through cellular structures and precipitates versus F/M through phase transformations and precipitates. Heat treatment optimization requires alloy-specific approaches, with A709 benefiting from controlled precipitation while G-92 requires careful phase transformation control. The results show that DED manufacturing can produce nuclear materials with exceptional performance, but directional effects and temperature-dependent behavior must be carefully considered for reactor component design and qualification.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

Engineering the fracture resistance of 2H-transition metal dichalcogenides using vacancies: An in-silico investigation based on HRTEM images

Vacancy engineering of 2H-transition metal dichalcogenides (2H-TMDs) has recently attracted great attention due to its potential to fine-tune the phonon and opto-electric properties of these materials. From a mechanical perspective, this symmetry-breaking process typically reduces the overall crack resistance of the material and adversely affects its reliability. However, vacancies can trigger the formation of heterogeneous phases that synergistically improve fracture properties. In this study, using MoSe2 as an example, we characterize the types and density of vacancies that can emerge under electron irradiation and quantify their effect on fracture. Molecular dynamic (MD) simulations, employing a re-parameterized Tersoff potential capable of accurately capturing bond dissociation and structural phase changes, reveal that isolated transition metal monovacancies or chalcogenide divacancies tend to arrest the crack tip and hence enhance the monolayer toughness. In contrast, isolated chalcogenide monovacancies do not significantly affect toughness. The investigation further reveals that selenium vacancy lines, formed by high electron dose rates, alter the crack propagating direction and lead to multiple crack kinking. Using atomic displacements and virial stresses together with a continuum mapping, displacement, strain, and stress fields are computed to extract mechanistic information, e.g., conditions for crack kinking and size effects in fracture events. The study also reveals the potential of specific defect patterns, “vacancy engineering,” to improve the toughness of 2H-TMDs materials.

36 MATERIALS SCIENCE↗
DOE OSTI2025

All 2D Material Printed Diodes and Circuits on Paper for Sustainable Electronics

Sustainable electronics aim to reduce environmental impact by using ecofriendly materials, energy-efficient manufacturing, and recyclable components. However, existing approaches rely on complex, resource-intensive methods, rare metals, or nanomaterials with limited stability as well as plastic substrates, raising sustainability issues. Solution-processed two-dimensional (2D) materials offer a promising alternative: water-based and biocompatible conductive, semiconductive, and insulating 2D material inks can be produced with scalable techniques and are suitable for the fabrication of fully printed devices on low-cost and biodegradable paper substrates. However, 2D material only and fully printed diodes on paper have not yet been reported. Here, we demonstrate fully inkjet-printed 2D material-based diodes on paper using metal-insulator-semiconductor and metal-insulator-metal-semiconductor architectures. Water-based graphene and MoS2 inks, prepared by liquid-phase exfoliation, are used for the metallic and insulating films, while electrochemical exfoliation is used to produce the semiconducting MoS2 ink. The highest forward-to-reverse current ratio obtained is 330 (at ±2 V), while the forward current density is 1 mA/cm2 (at 1 V), making the diode performance comparable to the best solution-processed diodes reported so far. However, in contrast to previous works, fabrication occurs entirely at room temperature and ambient pressure, without using high-pressure sputtering, thermal evaporation, and any precious metal ink. The devices maintain stable performance under bending up to strain of 4% over 10,000 cycles. Finally, the diodes are successfully integrated with other 2D-material based electrical components to realize fully printed RC circuits, differentiators, integrators, and AC-to-DC converters onto paper, hence demonstrating the suitability of our approach for sustainable and disposable integrated circuits.

2D materials↗
DOE OSTI2024

Development of “GaSb-on-silicon” metamorphic substrates for optoelectronic device growth

The epitaxial development and characterization of metamorphic “GaSb-on-silicon” buffers as substrates for antimonide devices is presented. The approach involves the growth of a spontaneously and fully relaxed GaSb metamorphic buffer in a primary epitaxial reactor, and use of the resulting “GaSb-on-silicon” wafer to grow subsequent layers in a secondary epitaxial reactor. The buffer growth involves four steps—silicon substrate preparation for oxide removal, nucleation of AlSb on silicon, growth of the GaSb buffer, and finally capping of the buffer to prevent oxidation. This approach on miscut silicon substrates leads to a buffer with negligible antiphase domain density. The growth of this buffer is based on inducing interfacial misfit dislocations between an AlSb nucleation layer and the underlying silicon substrate, which results in a fully relaxed GaSb buffer. A 1 μm thick GaSb layer buffer grown on silicon has ~9.2 × 10 7 dislocations/cm 2 . The complete lack of strain in the epitaxial structure allows subsequent growths to be accurately lattice matched, thus making the approach ideal for use as a substrate. Here we characterize the GaSb-on-silicon wafer using high-resolution x-ray diffraction and transmission electron microscopy. The concept’s feasibility is demonstrated by growing interband cascade light emitting devices on the GaSb-on-silicon wafer. The performance of the resulting LEDs on silicon approaches that of counterparts grown lattice matched on GaSb.

36 MATERIALS SCIENCE↗
DOE OSTI2024

GeSn-on-GaAs with photoconductive carrier lifetime >400 ns: role of substrate orientation and atomistic simulation

Group IV GeSn quantum material finds application in electronics and silicon-compatible photonics. Synthesizing these materials with low defect density and high carrier lifetime is a potential challenge due to lattice mismatch induced defects and tin segregation at higher growth temperature. Recent advancements in the growth of these GeSn materials on Si, Ge, GaAs, and with substrate orientations, demonstrated different properties using epitaxial and chemical deposition methods. This article addresses the effect of GaAs substrate orientation and misorientation on the materials’ properties and carrier lifetimes in epitaxial Ge 0.94 Sn 0.06 layers. With starting GaAs substrates of (100)/2°, (100)/6°, (110) and (111)A orientations, Ge 0.94 Sn 0.06 epitaxial layers were grown with an intermediate Ge buffer layer by molecular beam epitaxy and analyzed by several analytical tools. X-ray analysis displayed good crystalline quality, and Raman spectroscopy measurements showed blue shifts in phonon wavenumber due to biaxial compressive strain in Ge 0.94 Sn 0.06 epilayers. Cross-sectional transmission electron microscopy analysis confirmed the defect-free heterointerface of Ge 0.94 Sn 0.06 /Ge/GaAs heterostructure. Minority carrier lifetimes of the unintentionally doped n-type Ge 0.94 Sn 0.06 epilayers displayed photoconductive carrier lifetimes of >400 ns on (100)/6°, 319 ns on (100)/2°, and 434 ns on (110) GaAs substrate at 1500 nm excitation wavelength. On the other hand, Ge 0.94 Sn 0.06 layer showed poor carrier lifetime on (111)A GaAs substrate. The observed differences in carrier lifetimes were correlated with the formation energy of the Ge on (100)/6° and (100)/2° GaAs heterointerface using Stillinger–Weber interatomic potential model-based atomistic simulation with different heterointerfacial bonding by Synopsys QuantumATK tool. Total energy computation of 6280-atom Ge/GaAs supercell on (100)/6° leads to lower formation energy than (100)/2°, making it more thermodynamically stable. Hence, the growth of the GeSn/III–V material system using misoriented (100) substrates that are more thermodynamically stable will enhance the performances of optoelectronic devices.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1999

Production and Characterization of Carbon Nanotubes and Nanotube-Based Composites

The Nobel Prize winning discovery of the Buckuball (C60) in 1985 at Rice University by a group including Dr. Richard Smalley led to the whole new class of carbon allotropes including fullerenes and nanotubes. Especially interesting from many viewpoints are single-walled carbon nanotubes, which structurally are like a single graphitic sheet wrapped around a cylinder and capped at the ends. This cylinders have diameter as small as 0.5 - 2 nm (1/100,000th the diameter of a human hair) and are as long as 0.1 - 1 mm. Nanotubes are really individual molecules and believed to be defect-free, leading to high tensile strength despite their low density. Additionally, these fibers exhibit electrical conductivity as high as copper, thermal conductivity as high as diamond, strength 100 times higher than steel at one-sixth the weight, and high strain to failure. Thus it is believed that developments in the field of nanotechnology will lead to stronger and lighter composite materials for next generation spacecraft. Lack of a bulk method of production is the primary reason nanotubes are not used widely today. Toward this goal JSC nanotube team is exploring three distinct production techniques: laser ablation, arc discharge and chemical vapor deposition (CVD, in collaboration with Rice University). In laser ablation technique high-power laser impinges on the piece of carbon containing small amount of catalyst, and nanotubes self-assemble from the resulting carbon vapor. In arc generator similar vapor is created in arc discharge between carbon electrodes with catalyst. In CVD method nanotubes grow at much lower temperature on small catalyst particles from carbon-containing feedstock gas (methane or carbon monoxide). As of now, laser ablation produces cleanest material, but mass yield is rather small. Arc discharge produces grams of material, but purity is low. CVD technique is still in baby steps, but preliminary results look promising, as well as perspective of scaling the process up. We are also working on necessary purification of nanotubes. Applications of nanotubes are in such various fields as lightweight composites, molecular electronics, energy storage (electrodes in Li ion batteries), flat panel displays, conductive polymers, etc. JSC nanotube team is focused on development of lightweight materials. We work on the injection thermoset epoxies reinforced with nanotubes. Early results show good wetting of nanotube surface with epoxy, which is very important. More research will be possible as more nanotubes become available.

Nikolaev, Pavel↗
DOE OSTI2024

A Self-Healing, Flowable, Yet Solid Electrolyte Suppresses Li-Metal Morphological Instabilities

In this article, lithium metal (Li 0 ) solid-state batteries encounter implementation challenges due to dendrite formation, side reactions, and movement of the electrode–electrolyte interface in cycling. Notably, voids and cracks formed during battery fabrication/operation are hot spots for failure. Here, a self-healing, flowable yet solid electrolyte composed of mobile ceramic crystals embedded in a reconfigurable polymer network is reported. This electrolyte can auto-repair voids and cracks through a two-step self-healing process that occurs at a fast rate of 5.6 µm h -1 . A dynamical phase diagram is generated, showing the material can switch between liquid and solid forms in response to external strain rates. The flowability of the electrolyte allows it to accommodate the electrode volume change during Li 0 stripping. Simultaneously, the electrolyte maintains a solid form with high tensile strength (0.28 MPa), facilitating the regulation of mossy Li 0 deposition. The chemistries and kinetics are studied by operando synchrotron X-ray and in situ transmission electron microscopy (TEM). Solid-state NMR reveals a dual-phase ion conduction pathway and rapid Li + diffusion through the stable polymer-ceramic interphase. This designed electrolyte exhibits extended cycling life in Li 0 –Li 0 cells, reaching 12 000 h at 0.2 mA cm -2 and 5000 h at 0.5 mA cm -2 . Furthermore, owing to its high critical current density of 9 mA cm -2 , the Li 0 –LiNi 0.8 Mn 0.1 Co 0.1 O 2 (NMC811) full cell demonstrates stable cycling at 5 mA cm -2 for 1100 cycles, retaining 88% of its capacity, even under near-zero stack pressure conditions.

25 ENERGY STORAGE↗
DOE OSTI2025

Experimental and theoretical investigation of the crystalline surface, film, and interface properties of antiperovskite Mn 3 GaN grown by molecular beam epitaxy on MgO(001)

Here, we present a study of the epitaxial growth, characterization, and theoretical modeling of thin film antiperovskite Mn 3 GaN, an antiferromagnetic material with kagome structure which is grown on MgO (001) substrates using N-plasma-assisted molecular beam epitaxy. Reflection high energy electron diffraction is used to assess the in-plane evolution of the film structure during growth, and the surface is investigated in-situ using scanning tunneling microscopy and Auger electron spectroscopy. These results are combined with precision measurements done ex-situ determining the film lattice constants using a combination of x-ray diffraction with reciprocal space mapping and scanning transmission electron microscopy. Overall, a uniform, homogeneous film with an atomically smooth vacuum surface and atomically sharp substrate interface is found having very small in-plane tensile strain and mild out-of-plane compressive strain. First-principles theoretical calculations are applied in order to ascertain the lowest energy models for both the Mn 3 GaN surface and the Mn 3 GaN/MgO film/substrate interface. Models including MnGa versus MnN surface layers and MnGa versus MnN interfacial layers are considered as functions of both the Mn and Ga chemical potentials. The predictions are discussed in comparison to the experimental results. The overall findings suggest that Mn 3 GaN on MgO(001) is a viable epitaxial film which can be further explored in connection with antiferromagnetic spintronics.

Density functional theory↗