Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Diffuse interface methods”

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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Low-energy radiative backgrounds in CCD-based dark-matter detectors

The reach of sub-GeV dark-matter detectors is at present severely affected by low-energy events from various origins. We present the theoretical methods to compute the single- and few-electron events that arise from secondary radiation emitted by high-energy particles as they pass through detector materials and perform a detailed simulation to quantify them at (Skipper) CCD-based experiments, focusing on the SENSEI data collected at Fermilab near the MINOS cavern. The simulations account for the generation of secondaries from Cherenkov and luminescent recombination radiation; photo-absorption in the bulk, backside layer, pitch adapter, and epoxy; the photon reflection and refraction at interfaces; thin-film interference; the roughness of the interfaces; the dynamics of charges produced in the highly doped CCD-backside-layers; and the partial charge collection on the CCD backside. We consider several systematic uncertainties, notably those stemming from the backside modeling, which we estimate with a “fiducial” and an “extreme” charge-diffusion model, with the former model being preferred due to better agreement with partial-charge collection data. We find that Cherenkov photons constitute about 30% of the observed single-electron events for both diffusion models; radiative recombination contributes negligibly to the event rate for the fiducial model, although it can dominate over Cherenkov for the extreme model. We also estimate the fraction of 2-electron events that arise from 1-electron event coincidences in the same pixel, finding that the entire 2-electron rate can be explained by coincidences of radiative events and spurious charge. Accounting for both radiative and non-radiative backgrounds, we project the sensitivity of future Skipper-CCD-based experiments to different dark-matter models. For light-mediator models with dark-matter masses of 1, 5, and 10 MeV, we find that future experiments with 10-kg-year exposures and successful background mitigation could have a sensitivity that is larger by 9, 3, and 2 orders of magnitude, respectively, when compared to an experiment without background improvements.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Boosting the Low-Temperature Performance of Graphite Anodes by Creating an Electrochemically Active Interface

Graphite is the major anode material used in commercial lithium-ion batteries (LIBs). However, the sluggish ion-transfer kinetics associated with graphite anodes significantly restrict the operation of LIBs over a wide temperature. This is primarily due to their low reversible capacity and the substantial overpotential exhibited under low-temperature conditions. To address this limitation, we demonstrate herein an approach that involves grafting an electrochemically active lithium benzenesulfonate layer onto a graphite surface through a typical reduction reaction of diazonium cations, followed by ion exchange process. This surface modification reduces the charge transfer resistance of graphite anodes, leading to an excellent reversible capacity of ~150 mAh g –1 at low-temperatures (-20 °C, 0.1C). Electrochemical impedance spectroscopy indicates that both desolvation of the lithium ions outside the graphite, and lithium diffusion within the solid electrolyte interphase and graphite lattice are two crucial rate-limiting steps during the Li (de)lithiation, with the latter dominating during the low-temperature operation. In conclusion, these findings demonstrate a facile method for enhancing the low-temperature performance of graphite through surface modification and provide valuable insights into fundamental understandings that can guide the future design of better -low-temperature graphite anodes.

25 ENERGY STORAGE↗
DOE OSTI2023

Neural Network Enhanced RKPM for Electrochemical-Mechanical Coupled Damage Modeling of Energy Storage Materials

Energy storage materials undergo significant charge cycling, which makes understanding their reliability and durability fundamental in predicting performance and service life. Strong electrochemical-mechanical coupling and highly anisotropic material properties contribute to the formation and propagation of micro-cracking, largely along material interfaces and grain boundaries. For Li-ion batteries, for example, lithium moving between electrodes during charging and discharging process causes expansion and contraction of grains, and the strongly anisotropic and nonlinearly [Li]-dependent grain material properties can cause grains to expand into and contract away from each other, leading to chemo-mechanical cracking. In the first part of this work, a RKPM based computational framework for solving the coupled solid-phase lithium conservation with Fickian diffusion and the lithium concentration dependent anisotropic mechanical problem subjected to a highly nonlinear Butler-Volmer boundary condition is introduced. The choice of RKPM completeness conditions for lithium concentration and mechanical deformation fields, and the variational consistency condition for the domain integration of the coupled problem is first determined. In the second part of this work, a neural network-enhanced reproducing kernel particle method (NN-RKPM) [1] is leveraged to accurately capture damage and crack propagation throughout the material, by learning the location, orientation, and sharpness of discontinuity while allowing for a coarser nodal distribution than that is necessary for capturing sharp solution transitions using traditional mesh-based methods. NN-RKPM is used to inform how crack opening and closure in turn affect the coupled chemical equations and material microstructure.

damage modeling↗
DOE OSTI2021

Analysis of a new multispecies tumor growth model coupling 3D phase-fields with a 1D vascular network

In this study, we present and analyze a mathematical model for tumor growth incorporating ECM erosion, interstitial flow, and the effect of vascular flow and nutrient transport. The model is of phase-field or diffused-interface type in which multiple phases of cell species and other constituents are separated by smooth evolving interfaces. The model involves a mesoscale version of Darcy’s law to capture the flow mechanism in the tissue matrix. Modeling flow and transport processes in the vasculature supplying the healthy and cancerous tissue, one-dimensional (1D) equations are considered. Since the models governing the transport and flow processes are defined together with cell species models on a three-dimensional (3D) domain, we obtain a 3D–1D coupled model.

3D–1D coupled blood flow models↗
DOE OSTI2021

Toward predictive permeabilities: Experimental measurements and multiscale simulation of methanol transport in Nafion

A polymer membrane's permeability to solutes determines its suitability for various applications: a permeability value is essential for predicting performance in diverse contexts. Using aqueous methanol permeation through Nafion as an example, we describe a methodology for determining membrane permeability that accounts for boundary layer effects and the possibility of swelling. For the materials and apparatus used herein, analysis of a permeance measurement and computational fluid dynamics simulations show that the mass transfer boundary layer is on the order of ones to tens of microns. Additionally, the data are used to develop and validate a multiscale model describing solute permeation through a hydrated membrane as a series of physical mechanistic steps: reversible adsorption from solution at the membrane interface, diffusion driven by a concentration gradient within the membrane, and reversible desorption into solution at the opposite membrane interface. The validated model is used to predict methanol transport across a solar-driven CO 2 reduction device and to assess the impact of polymer changes on the measured value. The approach of combining experimental data, computational fluid dynamics, and the mechanistic multiscale model is expected to provide more accurate analysis of membrane permeation data in cases with polymer swelling or unusual device geometries, among others.

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

Regulating hydrogel mechanical properties with an electric field

Stimuli-responsive polymeric materials have attracted significant attention due to their ability to change properties in response to various external stimuli. Using an electric field as the stimulus is of particular interest as it possesses the potential for seamless integration of materials with electronic systems. While many materials with electric field responsive actuation have an associated mechanical property change, it is beneficial to develop materials that exhibit mechanical property changes without accompanying significant shape deformation. To address this challenge, here we designed a semi-interpenetrating polymer network (semi-IPN) hydrogel system containing both polyelectrolytes and salt ions, which enables electric field induced changes in mechanical properties while minimizing actuation. We first successfully verified the viability of our design by removing salt ions through a diffusion-only method where we witnessed the stiffness increased to 4.5 times the initial value while still being highly deformable. After this, we applied an electric field to transport the salt ions out of the hydrogel, as shown by both Raman spectroscopy and scanning electron microscopy. We were able to show a time-dependent stiffness increase, the maximum of which was 5 times the original stiffness. We quantified ion transport and water-splitting in the hydrogel by both experiments and simulations. Following this, we showed functional system reversibility by reversing the direction of the current to reinject salt ions into the semi-IPN hydrogel and reducing its stiffness to below the initial value. It's worth noting that our simulations enable us to understand the governing mechanisms behind ion generation and salt transport that leads to mechanical property changes. Finally, we were able to fabricate a spatially variable stiffness haptic interface with our hydrogel, with demonstrated reversibility and cyclability. This research can possibly find applications in soft robotics and haptics and also inspire the development of bio-compatible electronics related devices.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Two-dimensional perovskite heterostructures for single crystal semiconductor devices

Two-dimensional (2D) perovskites have gained much attention lately owing to their excellent optoelectronic properties, chemical tunability, and environmental stability. Multiple methods have been devised to synthesize high quality 2D perovskite single crystals, and recent progress in fabricating its heterostructures is notable as well. In particular, with growing interest in 2D van der Waals heterostructures, 2D perovskites have become a strong candidate as a new building block for heterostructures to reveal unique physical properties across different interfaces. Until now, various heterostructure devices of 2D perovskite single crystals with other types of 2D materials such as transition metal dichalcogenides (TMDs) and graphene have been studied, which have shown intriguing results including interlayer excitons and enhanced electronic properties. Here, we introduce various synthetic approaches to realize 2D perovskite single crystals and unique characteristics of their single crystal heterostructures fabricated with precision, possessing sharp interfaces. Moreover, recent studies of semiconductor devices based on 2D perovskite single crystal heterostructures are discussed in-depth. New perspectives to further the horizon in the field of 2D perovskite heterostructures are suggested in this work including the consideration of metal–2D material van der Waals contact, application of dry transfer techniques, electric bias driven ion diffusion studies, and nanocrystal array fabrication. 2D perovskite heterostructure single crystal devices factoring in these novel perspectives will further uncover the true potential of these materials for highly efficient and stable semiconductor devices.

Physics↗
DOE OSTI2025

Mesoporous Thin Film Architectures: Addressing Material Demands through Molecular Self-Assembly

Mesoporous thin films spark interest across a wide range of disciplines due to their tunable nanostructures, large internal surface areas, and strong compatibility with planar optical, electronic, and microfluidic devices. While attention in the porous materials community has shifted toward macroporous or disordered nanoporous systems, a resurgence in mesoporous thin film research is underway, driven by new molecular self-assembly methods, advanced materials chemistry, and improved characterization techniques. The integration of high-χN block copolymer design, kinetically persistent micelle templating, and postdeposition processing protocols now allows control over structural parameters such as pore size, wall thickness, porosity, and connectivity. These advances have overcome many of the thermodynamic and processing constraints that previously limited widespread adoption. Rather than serving only as high-surface-area supports, mesoporous thin films are engineered as active interfaces where responsive chemistries and nanoscale confinement act in tandem. Embedding switchable ligands, thermoresponsive polymers, redox mediators, or ion-selective groups directly within the pore walls enables real-time control over transport, optical, and electrochemical properties. These capabilities open up new directions in adaptive coatings, gated membranes, and fast-response biosensors. To further expand their functional scope, mesoporous films are integrated into hierarchical and multicomponent architectures. Techniques such as triblock terpolymer templating, crack-directed assembly, and nanoimprint lithography allow for control over spatial organization on the micron and submicron scale and pore system orientation. This enables programmable anisotropy, enhanced molecular diffusion, and wavelength-selective photonic behavior, essential for next-generation sensing, catalysis, and energy applications. Such structural and functional complexity requires equally sophisticated characterization. Multimodal and in situ techniques can track material dynamics under operational conditions. Recent progress includes extended-range ellipsometric porosimetry (EP) for hierarchical architectures, vacuum EP for interface energetics, time-resolved EP for diffusion kinetics, and correlative AFM-SAXS mapping. The introduction of advanced neutron-based spectroscopies, particularly quasielastic neutron scattering (QENS), promises to provide real-time access to ion transport dynamics and segmental motion under nanoscale confinement, offering a path toward deeper mechanistic understanding of structure-performance correlations in mesoporous systems. This Account reflects the technical advances made and the interdisciplinary collaborations that have shaped our collective vision. The particular dimensions of mesopores enable us to subtly tune interactions at the molecular, interfacial, and mesoscopic levels that permit us to harness nanoconfinement. What emerges is a versatile, modular platform capable of chemical gating, energy transduction, and sensing with a level of tunability unmatched by other porous materials. We highlight critical challenges including the need for more robust large-area processing, a deeper understanding of dynamic behavior under cycling, and better integration with device-level architectures. Our strategies support the transition of mesoporous thin films into active high-performance components in next-generation energy, environmental, and biomedical systems.

oxides↗
DOE OSTI2024

Optical and microstructural characterization of Er 3+ doped epitaxial cerium oxide on silicon

Rare-earth ion dopants in solid-state hosts are ideal candidates for quantum communication technologies, such as quantum memories, due to the intrinsic spin–photon interface of the rare-earth ion combined with the integration methods available in the solid state. Erbium-doped cerium oxide (Er:CeO 2 ) is a particularly promising host material platform for such a quantum memory, as it combines the telecom-wavelength (~ 1.5 μm) 4f–4f transition of erbium, a predicted long electron spin coherence time when embedded in CeO 2 , and a small lattice mismatch with silicon. In this work, we report on the epitaxial growth of Er:CeO 2 thin films on silicon using molecular beam epitaxy, with controlled erbium concentration between 2 and 130 parts per million (ppm). We carry out a detailed microstructural study to verify the CeO 2 host structure and characterize the spin and optical properties of the embedded Er 3+ ions as a function of doping density. In as-grown Er:CeO 2 in the 2–3 ppm regime, we identify an EPR linewidth of 245(1) MHz, an optical inhomogeneous linewidth of 9.5(2) GHz, an optical excited state lifetime of 3.5(1) ms, and a spectral diffusion-limited homogeneous linewidth as narrow as 4.8(3) MHz. We test the annealing of Er:CeO 2 films up to 900 °C, which yields narrowing of the inhomogeneous linewidth by 20% and extension of the excited state lifetime by 40%.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Leveraging Artificial Intelligence to Expedite Antibody Design and Enhance Antibody–Antigen Interactions

This perspective sheds light on the transformative impact of recent computational advancements in the field of protein therapeutics, with a particular focus on the design and development of antibodies. Cutting-edge computational methods have revolutionized our understanding of protein–protein interactions (PPIs), enhancing the efficacy of protein therapeutics in preclinical and clinical settings. Central to these advancements is the application of machine learning and deep learning, which offers unprecedented insights into the intricate mechanisms of PPIs and facilitates precise control over protein functions. Despite these advancements, the complex structural nuances of antibodies pose ongoing challenges in their design and optimization. Our review provides a comprehensive exploration of the latest deep learning approaches, including language models and diffusion techniques, and their role in surmounting these challenges. We also present a critical analysis of these methods, offering insights to drive further progress in this rapidly evolving field. The paper includes practical recommendations for the application of these computational techniques, supplemented with independent benchmark studies. These studies focus on key performance metrics such as accuracy and the ease of program execution, providing a valuable resource for researchers engaged in antibody design and development. Through this detailed perspective, we aim to contribute to the advancement of antibody design, equipping researchers with the tools and knowledge to navigate the complexities of this field.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

The Application of Differential Scanning Calorimetry to Investigate Precipitation Behavior in Nickel-Base Superalloys Under Continuous Cooling and Heating Conditions

In this study, a suite of experimental tools and fast-acting, numerical-simulation techniques was used to quantify the precipitation behavior of three nickel-base superalloys: IN-100, LSHR, and 718. Experimental methods comprised differential scanning calorimetry (DSC) to establish the specific heat as a function of temperature and selected direct-resistance heating trials (using a Gleeble ® machine) to obtain samples for microstructural analysis. For the DSC experiments, each alloy was cooled at a prescribed constant rate (between 5 and 20 K/min) after an initial soak/equilibration in the high-temperature, single-phase (supersolvus) temperature regime. On-heating DSC trials beginning at ambient temperature were also performed on alloy 718 in three different starting conditions: super- δ -solvus solution treated and water quenched (denoted as ST), solution treated and aged (STA), and solution treated and overaged (STOA). DSC results, revealing the thermal signatures associated with the kinetics of precipitation of γ ' (IN-100, LSHR) or γ ' and γ " (718), were interpreted using a previously-developed fast-acting routine that treats concurrent nucleation, growth, coarsening, and dissolution. For these simulations, special attention was paid to various thermo-kinetic input parameters including equilibrium solvus-approach curves, bulk free energies of transformation, matrix-precipitate interface energies, and effective diffusivities. For the γ - γ ' superalloys (IN-100 and LSHR), estimates of precipitate volume fraction as a function of temperature from the specific-heat data revealed semi-quantitative agreement with simulation predictions. For the γ - γ '- γ " superalloy (718), simulation predictions of precipitate volume fractions were converted to specific heat as a function of temperature and showed semi-quantitative agreement with the direct measurements.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Deposition pressure dependent structural and optoelectronic properties of ex-situ boron-doped poly-Si/SiO x passivating contacts based on sputtered silicon

Among common methods to form polycrystalline silicon (poly-Si) films for passivating-contact solar cells, physical vapor deposition, in particular sputtering, is the safest one as it does not require any toxic gaseous precursors. One of the critical parameters to control the properties of sputtered silicon films is their deposition pressure. In this work, structural and optoelectronic characteristics of ex-situ boron-doped poly-Si/SiO x passivating contacts, formed from sputtered intrinsic amorphous silicon (a-Si) deposited at different pressures on top of SiO x /c-Si substrates and subjected to a high-temperature boron diffusion step, are investigated. The deposition rate and density of the as-deposited a-Si films increase with reducing pressure. Low-temperature photoluminescence spectra captured from the as-deposited samples at different pressures do not show typical emissions from hydrogenated a-Si. Meanwhile, their Fourier-transform infrared absorption spectra all show Si–H stretching modes, indicating that hydrogen had been initially incorporated into the chemical SiO x layers and eventually hydrogenated the a-Si/SiO x interfaces during the sputtering process. After the high-temperature boron-diffusion step, all hydrogen-related peaks disappear. Lower pressure films (1.5 and 2.5 mTorr) show more consistent improved performance after hydrogen treatments, compared to higher pressure films (4 and 5 mTorr). Overall, the resultant passivating contacts at 2.5 mTorr achieve a low single-side recombination current density Jo of ~9 fA/cm 2 , whereas their contact resistivity is still low at 15 mO cm 2 .

14 SOLAR ENERGY↗
DOE OSTI2022

Pauling’s rules for oxide-based minerals: A re-examination based on quantum mechanical constraints and modern applications of bond-valence theory to Earth materials

Since their introduction in 1929, Pauling’s five rules have been used by scientists from many disciplines to rationalize and predict stable arrangements of atoms and coordination polyhedra in crystalline solids; amorphous materials such as silicate glasses and melts; nanomaterials, poorly crystalline solids; aqueous cation and anion complexes; and sorption complexes at mineral-aqueous solution interfaces. The predictive power of these simple yet powerful rules was challenged recently by George et al. (2020), who performed a statistical analysis of the performance of Pauling’s five rules for about 5000 oxide crystal structures. They concluded that only 13% of the oxides satisfy the last four rules simultaneously and that the second rule has the most exceptions. They also found that Pauling’s first rule is satisfied for only 66% of the coordination environments tested and concluded that no simple rule linking ionic radius to coordination environment will be predictive due to the variable quality of univalent radii. We address these concerns and discuss quantum mechanical calculations that complement Pauling’s rules, particularly his first (radius sum and radius ratio rule) and second (electrostatic valence rule) rules. We also present a more realistic view of the bonded radii of atoms, derived by determining the local minimum in the electron density distribution measured along trajectories between bonded atoms known as bond paths, i.e., the bond critical point (rc). Electron density at the bond critical point is a quantum mechanical observable that correlates well with Pauling bond strength. Moreover, a metal atom in a polyhedron has as many bonded radii as it has bonded interactions, resulting in metal and O atoms that may not be spherical. O atoms, for example, are not spherical in many oxide-based crystal structures. Instead, the electron density of a bonded oxygen is often highly distorted or polarized, with its bonded radius decreasing systematically from ~1.38 Å when bonded to highly electropositive atoms like sodium to 0.64 Å when bonded to highly electronegative atoms like nitrogen. Bonded radii determined for metal atoms match the Shannon (1976) radii for more electropositive atoms, but the match decreases systematically as the electronegativities of the M atoms increase. As a result, significant departures from the radius ratio rule in the analysis by George et al. (2020) is not surprising. We offer a modified, more fundamental version of Pauling’s first rule and demonstrate that the second rule has a one-to-one connection between the electron density accumulated between the bonded atoms at the bond critical point and the Pauling bond strength of the bonded interaction. Pauling’s second rule implicitly assumes that bond strength is invariant with bond length for a given pair of bonded atoms. Many studies have since shown that this is not the case, and Brown and Shannon (1973) developed an equation and a set of parameters to describe the relation between bond length and bond strength, now redefined as bond valence to avoid confusion with Pauling bond-strength. Brown (1980) used the valence-sum rule, together with the path rule and the valence-matching principle, as the three axioms of bond-valence theory (BVT), a powerful method for understanding many otherwise elusive aspects of crystals and also their participation in dynamic processes. We show how a priori bond-valence calculations can predict unstrained bond-lengths and how bond-valence mapping can locate low-Z atoms in a crystal structure (e.g., Li) or examine possible diffusion pathways for atoms through crystal structures. In addition, we briefly discuss Pauling’s third, fourth, and fifth rules, the first two of which concern the sharing of polyhedron elements (edges and faces) and the common instability associated with structures in which a polyhedron shares an edge or face with another polyhedron and contains high-valence cations. The olivine [α-(MgxFe1–x)2SiO4] crystal structure is used to illustrate the distortions from hexagonal close-packing of O atoms caused by metal-metal repulsion across shared polyhedron edges. We conclude by discussing several applications of BVT to Earth materials, including the use of BVT to: (1) locate H+ ions in crystal structures, including the location of protons in the crystal structures of nominally anhydrous minerals in Earth’s mantle; (2) determine how strongly bonded (usually anionic) structural units interact with weakly bonded (usually cationic) interstitial complexes in complex uranyl-oxide and uranyl-oxysalt minerals using the valence-matching principle; (3) calculate Lewis acid strengths of cations and Lewis base strengths of anions; (4) determine how (H2O) groups can function as bond-valence transformers by dividing one bond into two bonds of half the bond valence; (5) help characterize products of sorption reactions of aqueous cations (e.g., Co2+ and Pb2+) and oxyanions [e.g., selenate (Se6+O4)2- and selenite (Se4+O3)2-] at mineral-aqueous solution interfaces and the important role of protons in these reactions; and (6) help characterize the local coordination environments of highly charged cations (e.g., Zr4+, Ti4+, U4+, U5+, and U6+) in silicate glasses and melts.

Geochemistry & Geophysics↗
DOE OSTI2020

Mechanisms of mechanochemical synthesis of cesium lead halides: pathways toward stabilization of α-CsPbI 3

Cesium lead iodide with cubic perovskite structure (α-CsPbI 3 ) is gaining significant interest in photovoltaic applications due to its excellent absorbance of the visible solar light and other attractive optoelectronic properties. However, the synthesis of stable α-CsPbI 3 poses a significant challenge. Mechanochemical synthesis is emerging as a suitable method for the preparation of cesium lead halides. This work investigates the ball milling-induced synthesis of cesium lead halides perovskite phase using halide mixing or doping approaches. The synthesis in the CsI + PbI 2 , CsBr + PbBr 2 , CsBr + PbI 2 , and CsI + PbI 2 + NdI 3 mixtures and halide exchange reactions in the CsPbBr 3 + 3KI and CsBr + PbBr 2 + 3KI systems are investigated to elucidate the mechanism of this process. Then, CsPb(I 1–x Br x ) 3 and CsPb( 1–y )Nd y I 3 materials with different x and y ratios are prepared, and their stability is probed in the air using light absorption spectroscopy. These results suggest that Nd doping is more efficient in the stabilization of the perovskite structure than partial replacement of iodine with bromine. Microstructure observations reveal the existence of two different product formation mechanisms depending on the mechanical properties of reactants. The results reveal that the milling temperature has a significant impact on the reaction kinetics. Here, the produced particles nucleate and grow at the reactant interface and retard the synthesis reaction by creating a diffusion barrier. Extended milling reduces the product particle size and creates fresh contact between reactants, thus facilitating reaction completion.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Development and Applications of an eReaxFF Force Field for Graphitic Anodes of Lithium-Ion Batteries

Graphene is one of the most promising materials for lithium-ion battery anodes due to its superior electronic conductivity, high surface area for lithium intercalation, fast ionic diffusivity and enhanced specific capacity. A reliable description of many battery processes requires an explicit description of electrochemical interactions involving electrons. A detailed atomistic modeling of electronic conduction and non-zero voltage simulations of graphitic materials require the inclusion of an explicit electronic degree of freedom. To enable large length- and time-scale simulations of electron conduction in graphitic anodes, we developed an eReaxFF force field concept describing graphitic materials with an explicit electron. The newly developed force field, verified against quantum chemistry-based data describing, amongst others, electron affinities and equation of states, reproduces the qualitative behavior of electron conductivity in pristine and imperfect graphitic materials at different applied temperatures and voltages. In addition, excess electron localization near a defect site estimated from eReaxFF simulations agree quite well with the corresponding density functional theory calculations. Here, our eReaxFF simulations show the initiation of lithium-metal-plating driven by electron transfer from the graphene surface to the exposed lithium ions demonstrating the method’s potential for studying lithium-graphene interactions with explicit electrons and explain many unresolved electrode and electrode-electrolyte interface processes.

25 ENERGY STORAGE↗
DOE OSTI2020

Converting 2D Nanofiber Membranes to 3D Hierarchical Assemblies with Structural and Compositional Gradients Regulates Cell Behavior

Abstract New methods are described for converting 2D electrospun nanofiber membranes to 3D hierarchical assemblies with structural and compositional gradients. Pore‐size gradients are generated by tuning the expansion of 2D membranes in different layers with incorporation of various amounts of a surfactant during the gas‐foaming process. The gradient in fiber organizations is formed by expanding 2D nanofiber membranes composed of multiple regions collected by varying rotating speeds of mandrel. A compositional gradient on 3D assemblies consisting of radially aligned nanofibers is prepared by dripping, diffusion, and crosslinking. Bone mesenchymal stem cells (BMSCs) on the 3D nanofiber assemblies with smaller pore size show significantly higher expression of hypoxia‐related markers and enhanced chondrogenic differentiation compared to BMSCs cultured on the assemblies with larger pore size. The basic fibroblast growth factor gradient can accelerate fibroblast migration from the surrounding area to the center in an in vitro wound healing model. Taken together, 3D nanofiber assemblies with gradients in pore sizes, fiber organizations, and contents of signaling molecules can be used to engineer tissue constructs for tissue repair and build biomimetic disease models for studying disease biology and screening drugs, in particular, for interface tissue engineering and modeling.

Chen, Shixuan↗
DOE OSTI2023

Modeling laser interactions with aluminum and tantalum targets using a hybrid atomistic-continuum model

A hybrid atomistic-continuum method can model the microstructure evolution of metals subjected to laser irradiation. This method combines classical molecular dynamics (MD) simulations with the two-temperature model (TTM) to account for the laser energy absorption and heat diffusion behavior. Accurate prediction of the temperature evolution in the combined MD-TTM method requires reliable accuracy in electron heat capacity, electron thermal conductivity, and electron–phonon coupling factor across the temperatures generated. This study uses the electronic density of states (DOS) obtained from first-principle calculations. The calculated electron temperature-dependent parameters are used in MD-TTM simulations to study the laser metal interactions in FCC and BCC metals and the phenomenon of laser shock loading and melting. This study uses FCC Al and BCC Ta as model systems to demonstrate this capability. When subjected to short pulsed laser shocks, the dynamic failure behavior predicted using temperature-dependent parameters is compared with the experimentally reported single-crystal and nanocrystalline Al and Ta systems. The MD-TTM simulations also investigate laser ablation and melting behavior of Ta to compare with the ablation threshold reported experimentally. In conclusion, this manuscript demonstrates that integrating the temperature-dependent parameters into MD-TTM simulations leads to the accurate modeling of the laser–metal interaction and allows the prediction of the kinetics of the solid–liquid interface.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Leak-tight crack repair for 304L stainless steel using friction surfacing

Here, friction surfacing (FS) was investigated as a method for sealing 50 μm wide thru-cracks in 3-mm-thick 304L stainless steel (304L) plate. Friction surfacing was performed using two sizes of 304L consumable rod, with diameters of 9.52 mm and 12.7 mm, on two substrate conditions (clean and oxidized). Friction surfacing was able to deposit a 600 μm thick coating and repair cracks to a depth of 100 μm - 200 μm below the original surface, when using the 12.7-mm-diameter consumable rod for both substrate conditions. Helium leak rates of 10 -10 atm-cc/s were achieved on crack repairs, designating them as leak-tight by the ANSI N14.5 standard. Optical and scanning electron microscopy were used to investigate the coating microstructure and bond interfaces. In this study, the 12.7-mm-diameter rod performed better with disruption of the oxide layer on the oxidized substrate and providing a more homogenous coating. Mechanical properties of the coated samples were evaluated by performing micro-indentation, tensile, bending, and adhesion testing. Higher hardness in the coatings was observed due to the fine equiaxed grains resulting from dynamic recrystallization. Complete consolidation of plastically deformed material on the substrate and a strong diffusion bond across the interface was observed. This was reflected in tensile and bending properties of the FS coatings being comparable to those of uncoated specimen. No delamination of the FS coating was observed during adhesion testing. The results demonstrate that friction surfacing is a viable option to seal cracks in stainless steels.

36 MATERIALS SCIENCE↗