Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Single crystal model”

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 1,045 records · Page 58
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

A super-resolution technique to analyze single-crystal inelastic neutron scattering measurements using direct-geometry chopper spectrometers

Direct-geometry time-of-flight chopper neutron spectroscopy is instrumental in studying dynamics in liquid, powder, and single crystal systems. We report here that real-space techniques in optical imagery can be adapted to obtain reciprocal-space super resolution dispersion for phonon or magnetic excitations from single-crystal neutron spectroscopy measurements. The procedure to reconstruct super-resolution energy dispersion of excitations relies on an accurate determination of the momentum and energy-dependent point spread function and a dispersion correction technique inspired by an image disparity calculation technique commonly used in stereo imaging. Here by, applying these methods to spinwave dispersion data from a virtual neutron experiment demonstrates ~5-fold improvement over nominal energy resolution.

47 OTHER INSTRUMENTATION↗
DOE OSTI2020

Electron Paramagnetic Resonance Characteristics and Crystal Structure of a Tutton Salt Analogue: Copper-Doped Cadmium Creatininium Sulfate

Electron paramagnetic resonance and crystallographic studies on copper-doped cadmium creatininium sulfate (CdCrnS) were undertaken to study the characteristics of a copper-hexahydrate complex in an organic analogue of Tutton's salt. X-ray diffraction experiments determined the crystal structure of CdCrnS at both 100 and 298 K. CdCrnS, like Tutton salt, crystallizes in the monoclinic space group P2 1 /n. The unit cell contains two cadmium hexahydrate complexes, four creatininium ions, four sulfates, and four additional solvation waters. Both crystallography and EPR find that the doped copper replaces the cadmium in the structure. Single-crystal EPR measurements at room temperature determined the g and copper hyperfine (A Cu ) tensors (principal values: g = 2.437, 2.134, and 2.080 and A Cu = -327, -84.8, and 7.33 MHz). EPR spectra of the powder at room temperature gave g = 2.448, 2.125, and 2.085 and A Cu = -315, -75.0, and 35.0 MHz and at 110 K gave g = 2.462, 2.116, and 2.077 and A Cu = -340, -30.0, and 35.0 MHz. The room-temperature tensors are close to the "rigid lattice limit" values found in copper-doped Tutton salts but with a higher gmin and weaker A Cu x coupling than average. A small but measurable temperature dependency of the tensors indicated the presence of a dynamic Jahn-Teller (JT) effect. In addition, the EPR line width changed dramatically with temperature, which is like that found in all copper-doped Tutton crystals. Utilizing the model of Silver-Getz for the g-value variation gave an estimate for the energy difference (δ 12 = 640 cm -1 ) between the ground and next highest JT configurations. An empirical correlation appears to exist between δ 12 and g min and A Cu x for the copper hexahydrates studied in similar crystals. Finally, this suggests a relationship between the amount of unpaired spins in the copper d-orbital x lobe and the gap between wells of the adiabatic potential surface.

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

Atomistic Simulations of the Elastic Compression of Platinum Nanoparticles

Abstract The elastic behavior of nanoparticles depends strongly on particle shape, size, and crystallographic orientation. Many prior investigations have characterized the elastic modulus of nanoscale particles using experiments or simulations; however their reported values vary widely depending on the methods for measurement and calculation. To understand these discrepancies, we used classical molecular dynamics simulation to model the compression of platinum nanoparticles with two different polyhedral shapes and a range of sizes from 4 to 20 nm, loaded in two different crystal orientations. Multiple standard methods were used to calculate the elastic modulus from stress-vs-strain data for each nanoparticle. The magnitudes and particle-size dependence of the resulting moduli varied with calculation method and, even for larger nanoparticles where bulk-like behavior may be expected, the effective elastic modulus depended strongly on shape and orientation. Analysis of per-atom stress distributions indicated that the shape- and orientation-dependence arise due to stress triaxiality and inhomogeneity across the particle. When the effective elastic modulus was recalculated using a representative volume element in the center of a large nanoparticle, the elastic modulus had the expected value for each orientation and was shape independent. It is only for single-digit nanoparticles that meaningful differences emerged, where even the very center of the particle had a lower modulus due to the effect of the surface. These findings provide better understanding of the elastic properties of nanoparticles and disentangle geometric contributions (such as stress triaxiality and spatial inhomogeneity) from true changes in elastic properties of the nanoscale material.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Emergent quasi-one-dimensional antiferromagnetism in the distorted kagome magnet CePtPb

CePtPb hosts a distorted kagome lattice of Ce 3+ ions, providing a clean platform to investigate how reduced local symmetry and strong spin-orbit coupling reshape frustrated magnetism. Magnetization, specific heat, and magnetocaloric effect measurements, combined with a symmetry analysis of the single-ion anisotropy, demonstrate that the local 𝑚⁢2⁢𝑚 site symmetry selects a nearly Ising-like Kramers doublet with easy axes lying within the 𝑎⁢𝑏 plane. This results in three distinct in-plane Ising directions and an overall easy-plane anisotropy. The low-energy magnetic response is well captured by a three-sublattice Ising model, which quantitatively reproduces the saturation magnetization for arbitrary in-plane field orientations, including $[110]$ and $[1\bar{⁢1}⁢0]$, as well as the ratio of the field-induced critical fields. For 𝐵∥$[110]$, the phase diagram exhibits two quantum critical points at 𝐵 c⁢1 = 0.25T and 𝐵 c⁢2 = 0.55T, arising from the sequential polarization of the three Ising sublattices. In conclusion, these results reveal that the system develops quasi-one-dimensional spin chains along the 𝑐 axis, emerging from the nominally three-dimensional crystal structure composed of stacked kagome layers, and illustrate how reduced local symmetry can drive effective dimensional reduction in rare-earth Ising magnets.

Li, Fangli [Southern University of Science and Tec↗
DOE OSTI2026

MnRhBi3: A Cleavable Antiferromagnetic Metal

This dataset contains DFT input and output files supporting the theoretical modeling in the associated publication (Chem. Mater. 2024, 36, 11306-11316). The calculations characterize MnRhBi3, an orthorhombic (Cmmm) van der Waals-layered intermetallic compound that cleaves easily between neighboring Bi layers. The dataset is organized into three calculation types: (i) Bulk: Structural relaxations of the periodic MnRhBi3 crystal in antiferromagnetic (AFM) and ferromagnetic (FM) configurations, using the vdW-DF-optB86b functional. These provide the equilibrium lattice constants, magnetic energy differences (AFM is 0.5 meV/f.u. lower than FM), and magnetic moments (4.4 µB/Mn, 0.17 µB/Rh, 0.18 µB/Bi) reported in Table 1 of the main text. (ii) Slab: Same magnetic configurations computed with an 18 Ang vacuum layer introduced between Bi layers, used to calculate the cleavage energy Ec = 0.56 J/m2 (AFM) and 0.57 J/m2 (FM), establishing MnRhBi3 as a van der Waals-layered material comparable to graphite, MoS2, and CrI3. (iii) ELF: Single-point calculation on the relaxed bulk AFM geometry with LELF=.TRUE., producing the ELFCAR file used to generate electron localization function isosurfaces and contour maps (Fig. 2, main text) showing Bi lone pairs directed into the van der Waals gaps. All folders contain CONTCAR, INCAR, KPOINTS, OUTCAR, and POSCAR. The ELF/ folder additionally contains ELFCAR. Calculations were performed using VASP 6.3.2 with PBE + vdW-DF-optB86b, PAW potentials, and an energy cutoff of 800 eV.

36 MATERIALS SCIENCE↗
DOE OSTI2024

INSPIRED: Inelastic neutron scattering prediction for instantaneous results and experimental design

Inelastic neutron scattering (INS) has unique advantages in probing how atoms vibrate and how the vibrations propagate and interact. Such dynamic information is crucial in understanding various material properties, from heat capacity, thermal conductivity, phase transitions, and chemical reactions to more exotic quantum behavior. The analysis and interpretation of the INS spectra often start from a model structure of the sample, followed by a series of calculations to obtain the simulated spectra to compare with experiments. The conventional way to perform such calculations usually requires significant time, computing resources, and specialized expertise. Here, we present a new program named INSPIRED (Inelastic Neutron Scattering Prediction for Instantaneous Results and Experimental Design), which enables users to perform rapid INS simulations in several different ways on their personal computers in just a few clicks, with the crystal structure as the only input file. Specifically, the users can choose a pre-trained symmetry-aware neural network (coupled with an autoencoder) to predict the phonon density of states (DOS), 1D S(E) and 2D S(|Q|,E) spectra for any given structure. One can also choose an existing density functional theory (DFT) calculation from a database (containing over 12,000 crystals), and quickly obtain the simulated INS spectra for single crystals and powders. It is also possible to use pre-trained universal machine learning force fields to relax a given crystal structure, calculate the phonon dispersion and DOS, and, subsequently, the INS spectra. All these functions are implemented with a PyQt graphic user interface. Finally, we expect these new tools will benefit broad user communities and significantly improve the efficiency of experiment design, execution, and data analysis for INS.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Shock Hugoniot measurements of single-crystal 1,3,5-triamino-2,4,6-trinitrobenzene (TATB) compressed to 83 GPa

We present laser-driven shock Hugoniot measurements of single-crystal (SC) 1,3,5-triamino-2,4,6-trinitrobenzene (TATB) between 15 and 83 GPa, spanning pressures below and well above the Chapman–Jouguet pressure of ~28 GPa for TATB formulations (TATB grains mixed with plastic binders at 5–10 wt. %). The new SC data are generally ~3% more compressible than previously published data on neat and formulated TATB measured in gas-gun and explosive-driven experiments. An exception is at compressions in the density of ~1.5 (~30–40 GPa), where our new SC data exhibit significantly lower pressures than previous results on overdriven TATB formulations, suggesting that our SC samples remain largely unreacted below 35 GPa over the short nanosecond-time scales inherent to our laser-driven experiments. We report these novel equation-of-state measurements are a critical step toward understanding TATB in its most fundamental form and improving predictive modeling of TATB-based explosives.

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

Noncollinear antiferromagnetic order in the buckled honeycomb lattice of magnetoelectric Co 4 Ta 2 O 9 determined by single-crystal neutron diffraction

Co 4 Ta 2 O 9 exhibits a three-dimensional magnetic lattice based on the buckled honeycomb motif. It shows unusual magnetoelectric effects, including the sign change and nonlinearity. These effects cannot be understood without the detailed knowledge of the magnetic structure. Herein, we report neutron diffraction and direction-dependent magnetic susceptibility measurements on Co 4 Ta 2 O 9 single crystals. In this work, below 20.3 K, we find a long-range antiferromagnetic order in the alternating buckled and flat honeycomb layers of Co 2+ ions stacked along the c axis. Within experimental accuracy, the magnetic moments lie in the ab plane. They form a canted antiferromagnetic structure with a tilt angle of ~14° at 15 K in the buckled layers, while the magnetic moments in each flat layer are collinear. This is directly evidenced by a finite (0, 0, 3) magnetic Bragg peak intensity, which would be absent in the collinear magnetic order. The magnetic space group is C2/c'. It is different from the previously reported C2/c' group, also found in the isostructural Co 4 Nb 2 O 9 . The revised magnetic structure successfully explains the major features of the magnetoelectric tensor of Co 4 Ta 2 O 9 within the framework of the spin-flop model.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Structure-function analysis of microRNA 3'-end trimming by Nibbler

Nibbler (Nbr) is a 3'-to-5' exoribonuclease whose catalytic 3'-end trimming activity impacts microRNA (miRNA) and PIWI-interacting RNA (piRNA) biogenesis. Here, we report on structural and functional studies to decipher the contributions of Nbr’s N-terminal domain (NTD) and exonucleolytic domain (EXO) in miRNA 3'-end trimming. We have solved the crystal structures of the NTD core and EXO domains of Nbr, both in the apo-state. The NTD-core domain of Aedes aegypti Nbr adopts a HEAT-like repeat scaffold with basic patches constituting an RNA-binding surface exhibiting a preference for binding double-strand RNA (dsRNA) over single-strand RNA (ssRNA). Structure-guided functional assays in Drosophila S2 cells confirmed a principal role of the NTD in exonucleolytic miRNA trimming, which depends on basic surface patches. Gain-of-function experiments revealed a potential role of the NTD in recruiting Nbr to Argonaute-bound small RNA substrates. The EXO domain of A. aegypti and Drosophila melanogaster Nbr adopt a mixed α/β-scaffold with a deep pocket lined by a DEDDy catalytic cleavage motif. We demonstrate that Nbr’s EXO domain exhibits Mn2+-dependent ssRNA-specific 3'-to-5' exoribonuclease activity. Modeling of a 3' terminal Uridine into the catalytic pocket of Nbr EXO indicates that 2'-O-methylation of the 3'-U would result in a steric clash with a tryptophan side chain, suggesting that 2'-O-methylation protects small RNAs from Nbr-mediated trimming. Overall, our data establish that Nbr requires its NTD as a substrate recruitment platform to execute exonucleolytic miRNA maturation, catalyzed by the ribonuclease EXO domain.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Uniaxial stress study of spin and charge stripes in La1.875Ba0.125CuO4 by 139La NMR and 63Cu NQR

We study the response of spin and charge order in single crystals of La1.875Ba0.125CuO4 to uniaxial stress, through 139La nuclear magnetic resonance (NMR) and 63Cu nuclear quadrupole resonance (NQR), respectively. In unstressed La1.875Ba0.125CuO4, the low-temperature tetragonal structure onsets below TLTT = 57 K, while the charge order and the spin order transition temperatures are TCO = 54 K and TSO = 37 K, respectively. We find that uniaxial stress along the [110] lattice direction strongly suppresses TCO and TSO, but has little effect on TLTT. In other words, under stress along [110] a large splitting (≈ 21 K) opens between TCO and TLTT, showing that these transitions are not tightly linked. On the other hand, stress along [100] causes a slight suppression of TLTT but has essentially no effect on TCO and TSO. Magnetic field H along [110] stabilizes the spin order: the suppression of TSO under stress along [110] is slower under H ∥ [110] than H ∥ [001]. We develop a Landau free energy model and interpret our findings as an interplay of symmetry breaking terms driven by the orientation of spins.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Supramolecular Metal‐Organic Framework Electrocatalysts With Hydration‐Adaptive Gallery Expansion and Linker‐Mediated Distal Ni‐Co Cooperation Enable Framework‐Retentive Oxygen Evolution

Atomic‑level tailoring of electronic communication between spatially separated metal sites offers a route to accelerate multielectron electrocatalysis. We report Ni-TPTC, Co-TPTC, and heterobimetallic Ni-Co-TPTC supramolecular metal-organic frameworks (SMOFs) for the oxygen evolution reaction (OER) built from a terphenyl‑tetracarboxylate (TPTC) linker. Single-crystal X-ray diffraction of heterobimetallic Ni–Co–TPTC reveals trans-bis-aqua octahedral chains assembled into a π-stacked, hydrogen-bonded lattice and confirms retention of the parent architecture, while powder X-ray diffraction supports an isostructural Co analog. In 1.0 M KOH, Ni-Co-TPTC reaches 20 mA cm −2 at 350 mV and maintains operation at ∼60 mA cm −2 for 45 h with modest decay of the initial current. Time‑dependent ex situ diffraction shows dominant framework reflections with low‑angle shifts consistent with gallery dilation (d‑spacing ∼7.7 → ∼9.2 Å) during early operation, while Co K‑edge X‑ray absorption and microscopy reveal Co‑rich oxide/oxyhydroxide‑like domains at longer times that correlate with decay. X‑ray photoelectron spectroscopy, quantified by the Remote Binding‑Energy Differential, tracks composition‑dependent Ni-Co polarization, and DFT free‑energy analysis suggests that a model Ni-Co environment lowers the Co‑centered potential‑limiting step by ∼0.24 eV relative to the Co‑only framework. Furthermore, these results suggest that framework‑retentive expansion and distal Ni-Co coupling can drive OER activity before gradual deligation and oxide formation mark longer‑time degradation.

Supramolecular metal-organic frameworks (SMOFs)↗
DOE OSTI2020

An updated estimate of the Mu2e experiment sensitivity

The Mu2e experiment at Fermilab will search for the conversion of a negative muon into an electron inside the field of a nucleus. This process does not conserve charged-lepton flavour and is heavily suppressed in the Standard Model (SM), with a branching ratio < 10-50. Any evidence of it would be an undeniable evidence of new physics beyond the SM. The project sets out to achieve a single event sensitivity of $\sim 3 × 10^{-17}$ on the ratio between the probability for a conversion of a negative muon into an electron and the one for a muon capture by the nucleus. Such a sensitivity would represent a 4 orders of magnitude improvement on the previous upper limit for the process, making possible to test predictions of different extensions of the SM. Mu2e uses three superconducting solenoids to produce and measure the muon conversions. In the first solenoid, the Production Solenoid, pions and kaons are produced, together with other secondary products, by 8 GeV kinetic energy proton interactions in a tungsten target. A gradient magnetic field is specifically designed to direct low momentum particles into the Transport Solenoid, an S-shaped magnet that filters out particles with unwanted charge and momentum. Muons, produced by pion and kaon decays, finally reach the aluminum Stopping Target in the Detector Solenoid, where they eventually stop and convert. The result of the conversion process is a monochromatic electron of ~105 MeV/c momentum. The Detector Solenoid also hosts the two main detectors: a straw tube tracker and two CsI calorimeter disks, both providing measurement of event kinematics. A germanium detector and a LaBr crystal are located downstream of the Detector Solenoid to measure the X and gamma rays produced by the muon captures in the Stopping Target. A veto system of scintillators covering the Detector Solenoid and half of the Transport Solenoid is used to identify and reject cosmic rays interactions. With respect to the initial project, the Mu2e running plan has evolved to a staged configuration with 2 years at reduced intensity before the 2025 accelerator shutdown for the neutrino beam upgrade and 2 or 3 years at full intensity after that. This, together with geometry changes and a better knowledge of detector performances obtained by the first slice tests, has required a full revision of the signal over background selection that is the subject of this thesis. In order to achieve a new estimate of Mu2e sensitivity, the simulation of the data corresponding to the first 2 years of data acquisition has been performed. This includes both conversion electrons (CE) and the main sources of background: cosmics, decay-in orbits (DIO), radiative pion captures (RPC) and antiprotons. The characteristics of the signal and of the main backgrounds have been studied to define the best selection variables for CE. A special effort has been devoted to the evaluation of the antiproton background. The lack of experimental data for antiproton production cross section makes the systematic uncertainty on this background significant. A new parameterization of the cross section has been developed to fit the existing data and to provide a more reliable estimate of the systematic uncertainty by comparing the results of the old and the new model. A special effort has been devoted to the optimization of the antiprotons Monte Carlo generator. This study has also revealed that the dominant component within this background is represented by antiprotons produced in the opposite direction with respect to the Transport Solenoid entrance and then redirected to it by back-scattering processes in the Production Target. This ultimately highlights the sensitivity of the background estimate to the G4 handling of antiproton interactions in the 1 to 3 GeV/c momentum range. Finally, the final experimental sensitivity has been studied. The momentum and time selection have been optimized to obtain the 5-sigma discovery reach or, in case of no signal, the upper limit on conversion probability. The results confirm that, in the firs t two years of data taking, Mu2e will be able to improve the current experimental sensitivity for muon-to-electron conversion in an atomic field by more than 3 order of magnitudes.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
NASA NTRS2014

Amorphous and Crystalline H20 Ice at Rhea's Inktomi Crater

We present the analysis of Cassini spectral data from spectral mapping of Saturnian icy moons Dione and Rhea, to investigate possible effects of impact crater formation on the relative abundances of crystalline and amorphous water ice in the moons' ice crusts. Both moons display morphologically young ray craters as well as older craters. Possible changes in ice properties due to crater formation are conjectured to be more visible in younger craters, and as such Rhea's well imaged ray crater Inktomi is analysed, as are older craters for comparison. We used data from Cassini's Visual and Infrared Mapping Spectrometer (VIMS). For each pixel in the VIMS maps, spectral data were extracted in the near-infrared range (1.75 micrometers less than lambda less than 2.45 micrometers). Analysis was begun by fitting a single Gaussian to the peak in absorption at 2.0 micrometers, which was then subtracted from the data, leaving residuals with a minimum on either side of the original 2.0-micrometers band. The spectra of the individual spatial pixels were then clustered by the differences between these minima, which are sensitive to changes in both ice grain size and crystallinity. This yielded preliminary maps which approximated the physical characteristics of the landscape and were used to identify candidates for further analysis. Spectra were then clustered by the properties of the 1.5-micrometers band, to divide the map into regions based on inferred grain size. For each region, the predicted differences in minima from the Gaussian residuals, over a range of crystallinities, were calculated based on the found grain sizes. This model was used to find the crystallinity of each pixel via grain size and characteristics of the residual function. Preliminary results show a greater degree of crystallization of young crater interiors, particularly in Rhea's ray crater Inktomi, where ice showed crystalline ice abundances between 33 percent and 61 percent. These patterns in ice crystallization are possibly attributable to increased heat generated during crater formation.

Saturn Satellites↗
DOE OSTI2022

High‐Pressure and High‐Temperature Single‐Crystal Elasticity of Cr‐Pyrope: Implications for the Density and Seismic Velocity of Subcontinental Lithospheric Mantle

Abstract Single‐crystal X‐ray diffraction and Brillouin spectroscopy experiments were performed on a natural Cr‐pyrope (Prp 71.0 Alm 12.6 Sps 0.7 Grs 3.5 Uvr 12.2 ) at high pressure and high temperature up to 11.0 GPa and 800 K. Fitting the collected data to the third‐order finite strain equation yields bulk modulus ( K S 0 ), shear modulus ( G 0 ), their pressure ((∂ K S /∂ P ) T and (∂ G /∂ P ) T ) and temperature (( ∂K S /∂T ) P and ( ∂G/∂T ) P ) derivatives, K S 0 = 167.7(8) GPa, G 0 = 91.5(5) GPa, (∂ K S /∂ P ) T = 4.3(1), (∂ G /∂ P ) T = 1.4(1), ( ∂K S /∂T ) P = − 0.0175(1) GPa/K and ( ∂G/∂T ) P = − 0.0073(1) GPa/K. Using the obtained results, we examined whether the elastic properties of the Cr‐pyrope can be accurately calculated from those of endmembers including pyrope, almandine, grossular, and uvarovite assuming a linear relationship between elastic properties and composition (end‐member model). The results indicate that the end‐member model provides a sufficient approximation for the elastic properties of Cr‐pyrope in calculating the density and velocity of the subcontinental lithospheric mantle (SCLM). We modeled the densities and velocities of three typical types of SCLM (Archon, Proton, and Tecton) in order to investigate how the variation of chemical composition influences the SCLM. We obtained that the compositional change from the Archon to the Tecton increases the density of the SCLM significantly, which can be an important prerequisite for SCLM delamination. However, the compositional variation only slightly changes the velocity of the SCLM and the change is within the uncertainty of the calculation. Moreover, in comparison to the velocity, ρ / V P and ρ / V S are much more sensitive to the compositional change of the SCLM.

Xu, Jingui↗
DOE OSTI2024

Characterization of electromigration-induced short-range stress development in Al(0.25 at. % Cu) conductor line

Scanning x-ray microbeam topography and fluorescence experiments were conducted in situ to study the electromigration behavior of a 0.5 μm thick, 10 μm wide, and 200 μm long Al(0.25 at. % Cu) conductor line with 1.5 μm-thick SiO2 passivation on a single crystal Si substrate. The strain sensitivity of x-ray topography measurement allowed detailed examination of the electromigration-induced stress distribution and evolution in the conductor line in response to the depletion of Cu solute early in the electromigration process. Upon electromigration at 0.4 MA/cm2 and 303 °C, a short-range stress gradient was quickly induced by Al migration in the Cu-depleted cathode region to counteract further Al flow. The stress gradient was fully developed during the 5.3 h incubation time, extending over the critical Blech length of about 66 μm from the cathode end. Plastic deformation then occurred at the downstream end of the Cu-depleted region. The preferential electromigration of Cu did not cause detectable stress change outside the Cu-depleted region, except for the significant stress development from the Al2Cu precipitation at the anode end which appeared to initiate the fracture in the passivation. Preliminary finite difference modeling was undertaken to simulate the experimental observations, from which important parameters dictating electromigration in Al(Cu) line were extracted: an apparent effective valence of −5.6 and −1.9 for Cu and Al in Al(Cu), respectively, and a critical Cu concentration of 0.16 at. % above which Al grain boundary diffusion is effectively blocked.

Physics↗
DOE OSTI2021

Computational study of Li3BO3 and Li3BN2 I: Electrolyte properties of pure and doped crystals

Both Li3BO3 and Li3BN2 materials have promising properties for use in all-solid-state batteries and other technologies dependent on electrolytes with significant ionic conductivity. As the first of a two-part study, this paper reports the analysis of detailed simulations of Li ion diffusion in the monoclinic forms of these materials. Using both NEB and MD methods, it is clear that Li ion migration via vacancy mechanisms provides the most efficient ion transport in each material. While the results suggest that interstitial defects in these materials do not play a direct role in Li ion migration, their relative stability seems to enhance vacancy production via the formation of Frenkel-type defects. This may partially explain why the Li ion conductivities computed from MD simulations of samples initially containing a single Li ion vacancy are in reasonable agreement with measured values of this work for Li3BO3 and those reported in the literature for poorly crystalline samples of both materials. The possibility of increasing vacancy concentrations by substitutional doping (F for O in Li3BO3 and C for B in Li3BN2) is also examined, finding simulated conductivities comparable to those of the ideal vacancy model.

Li, Yan↗
DOE OSTI2022

Integrated Computational Materials and Mechanical Modeling for Additive Manufacturing of Alloys with Graded Structure Used in Fossil Fuel Power Plants

Wire-arc additive manufacturing (WAAM) has demonstrated its unique capability of producing large-size alloy components with a significantly reduced fabrication time and enhanced geometry design freedom. In this project, the team has developed an ICME (Integrated Computational Materials Engineering) modeling framework, which supports the WAAM of the AUSC (Advanced Ultra-Supercritical) power plant components. The manufacturing design has been applied to Inconel 740H, steel P91, as well as the dissimilar alloy components between steel P91 and Inconel 740H. The ICME model framework is developed by considering two types of modeling. First, mechanistic modeling has been applied to control the printing quality and understand the sequence of the dissimilar printing of the wall structure. The following models have been included in the developed ICME framework: finite element thermal model, grain structure model, residual stress simulation, crystal plasticity model, CALPHAD-based precipitation kinetic model, phase stability prediction, thermal expansion predictive model, and heuristic creep model. Secondary, a physics-based machine learning model has also been developed based on the ICME model structure. The machine learning model development is based on the ICME model prediction with calibration of the experiments. In addition, the WAAM has been utilized as a high-throughput experimental tool rapidly generating a gradient of alloy composition to facilitate experimental database generation for process-structure-property relationships. Such a database directly supported the ICME-enhanced machine learning, which further assisted in intermediate composition block design between P91 and 740H. A high-throughput screening study of the oxidation resistance has been performed based on such high-throughput experimentation. Based on the computational design, several dissimilar alloy manufacturing with post-heat treatment have been performed with a comprehensive evaluation of mechanical performance, including hardness mapping, yield strength, creep resistance. In this project, the single component of P91 and 740H processed by WAAM after heat treatment designed by ICME has demonstrated higher performance in yield strength and creep resistance than the wrought materials. The P91 sample prepared by WAAM with ICME-designed heat treatment performs better than P92 in creep resistance. The designed graded alloy printing with intermediate block shows a promising performance that exceeds the traditional welding. Moreover, the current research indicates the high need for location-specific design analysis with uncertainty quantification, an important topic that deserves more dedicated research. The achievement of this project demonstrated the promising future of WAAM in structural alloy manufacturing for energy power plant development. Successful printing requires synergetic efforts made by manufacturing, mechanical, and materials sciences.

20 FOSSIL-FUELED POWER PLANTS↗
NASA NTRS1990

A numerical and analytical study of nonlinear bifurcations associated with the morphological stability of two-dimensional single crystals

The nonlinear stability of a two-dimensional single crystal of pure material in an undercooled melt is studied both analytically and numerically. The quasi-steady state approximation is used for the thermal fields, and the effects of different solid and liquid thermal conductivities, isotropic interfacial growth kinetics, and isotropic surface tension are included. The bifurcation analysis is performed by calculating the instantaneous value of the fundamental component of the local normal growth speed for an interface perturbed by a single Fourier shape component. Numerically, the fundamental component of the interfacial growth speed is found by Fourier analysis of the solution to an integrodifferential equation obeyed at the interface. Analytically, an expansion technique is used to derive a solvability condition defining each of these bifurcation points. The analytical and numerical results are in very close agreement. Almost all of the bifurcations are subcritical, and the results are presented by giving values of the Landau coefficient as a function of the different dimensionless parameters used in the model.

Brush, L. N.↗