Search NASA⌕ Search

SEARCH · Search NASA

Results for “extensive and intensive behavior”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Aggregate residential demand flexibility behavior: A novel assessment framework

Residential demand flexibility (DF) could save the U.S. electric grid up to 10 GW of peak demand while supporting increased amounts of renewable generation. However, less than 40% of the estimated DF peak reduction capacity is currently realized, and less than 8% of American households are enrolled. These low participation rates are combined with high rates of "overriding" a DF event and eventual opt-outs among enrolled customers. There is still not a comprehensive understanding of the drivers of DF participation and occupant interaction with DF program signals. We, therefore, present a novel survey processing framework to assess occupant DF-relevant behavior from the American Time Use Survey (ATUS). Using the framework, we summarize both the extensive and intensive behavior of more than 200,000 ATUS respondents (2003-2018 data) and provide insights on the DF-relevant behavior of residential occupants, which is generally overlooked in the literature. Here we use the framework to identify the activity priorities of residential occupants in the United States during different DF-relevant periods (critical peak, peak, and off-peak). These preferred activities capture overlooked routine behaviors that could be barriers to DF participation, if ignored, and might explain the high levels of overrides often exhibited by participants of demand response.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2025

Three-dimensional modeling of hyphal fusion, branching, and nutrient transport in filamentous fungi

Fungi exhibit behaviors distinct from other microbes. Filamentous fungi grow by extending complex networks of branched filaments collectively referred to as the mycelium. These networks can expand over large distances and traverse low-nutrient areas by translocating nutrients through the filament network. This spatial characteristic makes filamentous fungi crucial for soil ecosystems, supporting stable microbial communities and promoting plant growth. However, simulating these behaviors is complex. The elongated nature of fungal compartments results in different mechanical interactions compared to the commonly modeled spherical bacteria. These detailed hyphal mechanics require specialized consideration and are often excluded from conventional fungal simulation packages. Additionally, the extensive fungal networks in nature demand computationally intensive simulations, necessitating high-performance algorithms. Therefore, realistic fungi simulations require specialized software. Here, we introduce a fungal modeling expansion to the high-performance biological modelling and interface exchange (bmx) software suite. bmx leverages adaptive mesh refinement in AMReX for chemical diffusion and incorporates a full mechanical model for bacterial cells, accelerated by GPUs. By extending bmx to model filamentous particles, we demonstrate the formation of complex filament networks through interactions like hyphal branching and fusion (anastomosis). We show that the networks produced match real-world fungal structures through various metrics. This work supports computational studies of fungal growth dynamics and can be adapted to investigate the growth of other filamentous structures in biology or materials science. The expanded-BMX package is open-sourced and is available online.

Cell mechanics↗
DOE OSTI2025

Reorganization of Water at Aqueous Aluminum Chloride (AlCl 3 ) Interfaces: Vibrational Sum Frequency Generation and Molecular Dynamics Simulations

AlCl 3 hydration states and complexation are not well understood both in solutions and at the air–aqueous interface despite their potential significance in natural waters and their industrial and energy-related applications. Here, we investigated Al 3+ and Cl – ion behaviors in an AlCl 3 aqueous bulk solution and at the air–aqueous interface using interface-selective vibrational sum frequency generation (SFG), Raman and infrared spectroscopies, molecular dynamics (MD) simulation, as well as molecular-informed reduced modeling. Our reduced modeling reveals relatively long-range effects for Al 3+ as compared to monovalent ions such as Na + indicating that the interfacial depth of trivalent ions can be significantly larger than that of monovalent ions at the air–water interface. MD simulations reveal interfacial stratification and multiple layering of the ions. Compression of the Al 3+ and Cl – distributions with increasing concentrations from 0.5 to 2.5 m is also observed in the subsurface regions. Significant SSP- and PPP-polarized SFG OH spectral intensity increases are observed from 0.5 to 1.5 m and 0.5 to 2.5 m, respectively, indicative of interfacial depth increases and a change in average orientation above 1.5 m. Extensive evaluation of SFG spectra, Fresnel-corrected using several approaches, shows the same trends. The nonmonotonic trend points to a changing structure in surface and subsurface water orientation and hydrogen bonding environment generally consistent with the MD simulation of stratification and water orientation changes. Furthermore, solvent-shared ion pairing is implicated with MD simulation radial distribution analysis and consistent with infrared spectral identification of the hexaaqua aluminum ion in the solution phase. Spectral evidence of a strong Al 3+ hydration shell and the acidic behavior of the Al 3+ ions is obvious in the Raman and infrared spectra of the bulk solution. In conclusion, we show that the MD dipole potential is directly related to the MD second-order susceptibility of the interface, χ SFG–MD (2) , both of which correlate up to ∼35 Å with the spectral observations of increasing and then saturating intensities, suggesting that both ion stratification and interfacial depth determine the water orientations at an air–water interface of 1-3 electrolyte solutions.

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

Computational modeling of microalgal biofilm growth in heterogeneous rotating algal biofilm reactors (RABRs) for wastewater treatment

Rotating algal biofilm reactors (RABRs) are innovative systems designed to cultivate microalgae biofilms efficiently. In this paper, we have developed a novel mathematical model to accurately capture the growth dynamics of algae biofilms within RABR. By considering the spatial heterogeneity of the RABR, we introduce a PDE-based model that addresses the spatial variations across the substratum, enabling a more accurate simulation of biofilm growth in RABRs. The photosynthesis process is modeled through reactive kinetics, driving the growth of the algae biofilm. To analyze the system's behavior, we employ finite difference numerical methods to solve the complex PDE model. We then conduct extensive numerical simulations to understand algae biofilm growth in the RABR environment under various operational factors and environmental conditions. One primary focus in these simulations is to investigate the impact of various harvesting strategies, harvesting frequencies, light intensity, and light exposure on the overall biomass productivity of the algae biofilm. The numerical results provide valuable insights into optimizing algae biofilm growth and designing harvesting techniques in RABR systems. Our proposed novel mathematical model provides an effective platform for the theoretical investigation and design of RABRs for wastewater treatment.

09 BIOMASS FUELS↗
DOE OSTI2024

Reinforcement expectation in the honeybee ( Apis mellifera ): Can downshifts in reinforcement show conditioned inhibition?

When animals learn the association of a conditioned stimulus (CS) with an unconditioned stimulus (US), later presentation of the CS invokes a representation of the US. When the expected US fails to occur, theoretical accounts predict that conditioned inhibition can accrue to any other stimuli that are associated with this change in the US. Empirical work with mammals has confirmed the existence of conditioned inhibition. But the way it is manifested, the conditions that produce it, and determining whether it is the opposite of excitatory conditioning are important considerations. Invertebrates can make valuable contributions to this literature because of the well-established conditioning protocols and access to the central nervous system (CNS) for studying neural underpinnings of behavior. Nevertheless, although conditioned inhibition has been reported, it has yet to be thoroughly investigated in invertebrates. Here, we evaluate the role of the US in producing conditioned inhibition by using proboscis extension response conditioning of the honeybee (Apis mellifera). Specifically, using variations of a “feature-negative” experimental design, we use downshifts in US intensity relative to US intensity used during initial excitatory conditioning to show that an odorant in an odor–odor mixture can become a conditioned inhibitor. We argue that some alternative interpretations to conditioned inhibition are unlikely. However, we show variation across individuals in how strongly they show conditioned inhibition, with some individuals possibly revealing a different means of learning about changes in reinforcement. We discuss how the resolution of these differences is needed to fully understand whether and how conditioned inhibition is manifested in the honeybee, and whether it can be extended to investigate how it is encoded in the CNS. It is also important for extension to other insect models. In particular, work like this will be important as more is revealed of the complexity of the insect brain from connectome projects.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2024

Selective pseudocapacitive immobilization of REE elements on carbon based electrodes

Current extraction technologies for rare earth elements at ambient conditions are reagent and energy intensive, giving rise to significant quantities of secondary waste streams. Exploration of the removal of common cations from aqueous streams has been undertaken with carbon materials and have shown promising results. However, those promising materials have not been explored extensively for rare earth element capture from aqueous process streams. In this study, a carbon electrode was investigated for the adsorption and immobilization of REE elements from aqueous solutions. Cyclic voltammetry studies of the carbon electrode displayed a pseudocapacitive behavior where it was verified that cation adsorption is accompanied by an electron transfer process. Preliminary tests showed that a current density of 89.1 mA. g -1 , allowed for the pseudocapacitive adsorption (PSA) of Nd 3+ cation without the formation of rare earth hydroxides. Hence, the selective PSA of Nd 3+ was verified in an electrolyte solution with equimolar concentrations of Nd 3+ , Mg 2+ , Li + , Na + and K + , achieving separation factors of 8.6, 1.1, and 10.9 for Nd 3+ /Li + , Nd 3+ /Na + , and Nd 3+ /Mg 2+ , respectively. Analysis of the potential-time curves for the various cations suggests that storage of ions in the electrode involved pore-spacings rather than interlayer spacings and was corroborated by XRD analysis. Specific capacitance as high as >640 F.g -1 for Nd 3+ was also observed for the carbon electrode, with Faradaic efficiencies (FE) of > 20% for Nd 3+ and Mg 2+ , and > 9% for Na + . Furthermore, adsorption capacity of >125 mg g -1 after 4hrs of electrosorption was observed for Nd 3+ . In the presence of the Nd 3+ electrolyte, the electrode achieved increasing storage of Nd 3+ with >94% retention of Nd 3+ with a FE > 20% over 6 cycles between loading and releases in 1 M KCl solution. In conclusion, these preliminary results provide information and some guidance on the selective recovery of rare earth elements, such as Nd, in aqueous streams in the presence of competitive cations employing the use of carbon-based materials.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2026

From molecular to macroscopic: predicting liquid–liquid phase equilibria and small-angle scattering of mixtures of organic liquids from atomistic simulation using Kirkwood–Buff theory

Macroscopic phase equilibria between solutions define the functionality of many biological and industrial processes, yet they are challenging to predict due to the inherent complexity of liquids containing large molecules. This work introduces an approach for the purely predictive calculation of such phase equilibria in temperature-composition space from molecular dynamics (MD) simulations at one temperature in the single-phase region. We use an approach developed previously to obtain the entropic and enthalpic contributions to the free energy of mixing from the atomic-scale information given by MD simulations via Kirkwood–Buff theory. This allows us to accurately estimate the free energy of mixing as a function of temperature, and thus obtain liquid–liquid phase equilibria, including liquid–liquid critical points, associated binodal and spinodal lines, and composition fluctuations across a region of temperature and composition. Results for binary malonamide–alkane systems are validated by comparison to a direct experimental probe of the fluctuations: the small angle X-ray scattering intensity near zero wavenumber. The MDKB → Phase method demonstrated here provides a significant improvement in predicting liquid–liquid equilibria and free energy as a function of temperature for our systems of interest compared to conventional thermodynamic models. The accurate performance of this purely predictive approach lies in its preservation of atomistic details when determining thermodynamic properties. Furthermore, its inherent extensibility to multi-component systems will likely make the MDKB → Phase approach a valuable general tool for connecting molecular interactions to macroscopic phase equilibria and for the computational screening of materials for targeted thermodynamic behavior.

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

Understanding the High-Temperature Deformation Behaviors in Additively Manufactured Al6061+TiC Composites via In Situ Neutron Diffraction

Aluminum matrix composites (AMCs) are designed to enhance the performance of conventional aluminum alloys for engineering applications at both room and elevated temperatures. However, the dynamic phase-specific deformation behavior and load-sharing mechanisms of AMCs at elevated temperatures have not been extensively studied and remain unclear. Here, in situ neutron diffraction experiments are employed to reveal the phase-specific structure evolution of additively manufactured Al6061+TiC composites under compressive loading at 250 °C. It is found that the addition of a small amount of nano-size TiC significantly alters the deformation behavior and increases the strength at 250 °C in comparison to the as-printed Al6061. Unlike the two-stage behavior observed in Al6061, the Al6061+TiC composites exhibit three stages during compression triggered by changes in the interphase stress states. Further analysis of Bragg peak intensity and broadening reveals that the presence of TiC alters the dislocation activity during deformation at 250 °C by influencing dislocation slip planes and promoting dislocation accumulation. These findings provide direct experimental observations of the phase-specific dynamic process in AMCs under deformation at an elevated temperature. The revealed mechanisms provide insights for the future design and optimization of high-performance AMCs.

36 MATERIALS SCIENCE↗
DOE OSTI2025

AUTOMATIC GENERATION OF EVENT TREES AND FAULT TREES: A MODEL-BASED APPROACH

In the past few decades, increasing complexity in modern engineering systems has been driven by the integration of a large number of components and by the fact that the system operations involve many disciplines (e.g., thermal-hydraulics, plant operations, cyber-security). Current safety/reliability modeling approaches to such systems are labor intensive, difficult to learn, and rely heavily on simplistic Boolean logic to depict failure propagation and accident progression. While these methods serve well for simple systems (i.e., linear causal systems with limited small inter- and intra-system interactions), their results are difficult to verify when modeling complex systems (typically performed through the extensive use of modeling assumptions). The development of new methods is addressed to meet these challenges through a model-based system engineering (MBSE) lens. Under MBSE philosophy, every aspect of the system (form or function) is represented by a model that completely characterizes its architecture or behavior. MBSE approach greatly improves the management of design, analysis and verification of complex systems. An integration of Dynamic Probabilistic Risk Assessment (DPRA) methods with MBSE models is proposed to perform safety/reliability analyses of engineering systems. In particular, MBSE representation of the system (performed using Systems Modeling Language [SysML]) is coupled with DPRA methods to automatically generate event trees and fault trees.

97 - MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Magnetically propagating Hund’s exciton in van der Waals antiferromagnet NiPS3

Abstract Magnetic van der Waals (vdW) materials have opened new frontiers for realizing novel many-body phenomena. Recently NiPS 3 has received intense interest since it hosts an excitonic quasiparticle whose properties appear to be intimately linked to the magnetic state of the lattice. Despite extensive studies, the electronic character, mobility, and magnetic interactions of the exciton remain unresolved. Here we address these issues by measuring NiPS 3 with ultra-high energy resolution resonant inelastic x-ray scattering (RIXS). We find that Hund’s exchange interactions are primarily responsible for the energy of formation of the exciton. Measuring the dispersion of the Hund’s exciton reveals that it propagates in a way that is analogous to a double-magnon. We trace this unique behavior to fundamental similarities between the NiPS 3 exciton hopping and spin exchange processes, underlining the unique magnetic characteristics of this novel quasiparticle.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Understanding High-Voltage Behavior of Sodium-Ion Battery Cathode Materials Using Synchrotron X-ray and Neutron Techniques: A Review

Despite substantial research efforts in developing high-voltage sodium-ion batteries (SIBs) as high-energy-density alternatives to complement lithium-ion-based energy storage technologies, the lifetime of high-voltage SIBs is still associated with many fundamental scientific questions. In particular, the structure phase transition, oxygen loss, and cathode–electrolyte interphase (CEI) decay are intensely discussed in the field. Synchrotron X-ray and neutron scattering characterization techniques offer unique capabilities for investigating the complex structure and dynamics of high-voltage cathode behavior. In this review, to accelerate the development of stable high-voltage SIBs, we provide a comprehensive and thorough overview of the use of synchrotron X-ray and neutron scattering in studying SIB cathode materials with an emphasis on high-voltage layered transition metal oxide cathodes. We then discuss these characterizations in relation to polyanion-type cathodes, Prussian blue analogues, and organic cathode materials. Finally, future directions of these techniques in high-voltage SIB research are proposed, including CEI studies for polyanion-type cathodes and the extension of neutron scattering techniques, as well as the integration of morphology and phase characterizations.

25 ENERGY STORAGE↗
DOE OSTI2024

Airborne LiDAR to Improve Canopy Fuels Mapping for Wildfire Modeling

Increasing conflict between wildfire and the built environment has increased the need for more up-to-date and finer resolution canopy fuels data to improve wildfire modeling and associated risk forecasts. The US Forest Service and US Department of the Interior’s LANDFIRE product, which provides 30-m resolution canopy fuels data for the entire US, is one of the most widely used sources of fuels data. However, the last complete mapping effort for LANDFIRE is based on 2016 conditions, and subsequent updates reflect disturbances 1-2 years behind the release year. Airborne systems equipped with Light Detection and Ranging (LiDAR) sensors can be deployed to actively sense canopy structure and estimate canopy fuels data (cover, height, base height, bulk density) at finer resolutions. Canopy base height (CBH) and canopy bulk density (CBD) are difficult to measure both in the field and in LiDAR point clouds. Still, they are important for accurately modeling crown fires, which are often intense and difficult to contain. Additionally, point cloud datasets are large, and calculations require efficient utilization of computational resources. To address these challenges, we are working on an approach that uses openly available National Ecological Observatory Network (NEON) airborne LiDAR data, with calculations processed in the R programming language and parallelized through the lidR package. CBH and CBD are often derived from tree height, diameter at breast height, and species-specific allometries using the Fire and Fuels Extension of the Forest Vegetation Simulator (FFE-FVS). We aim to test if airborne LiDAR can estimate CBH and CBD without the use of empirical equations. Reliable estimates of canopy fuels data directly from airborne LiDAR could streamline quick, fine-resolution updates for use in wildfire behavior models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Rigorous Oxidation State Assignments for Supported Ga-Containing Catalysts Using Theory-Informed X-ray Absorption Spectroscopy Signatures from Well-Defined Ga(I) and Ga(III) Compounds

Gallium-based heterogeneous catalysts originally developed for commercial use in propane dehydroaromatization have been explored extensively as potential replacements for Pt- and Cr-based catalysts used in propane dehydrogenation for on-demand propylene production. In a large number of experimental and theoretical studies, active sites with a variety of Ga nuclearities, coordination environments, and oxidation states have been proposed. Isolated Ga(I) ions are often invoked, despite the scarcity of well-defined molecular species and their well-documented instability. In this study, we investigate the appearance of a high-intensity, low-energy white line at the Ga K-edge upon reduction of Ga/HZSM-5 by H 2 at a temperature of ca. 500 °C, accompanied by a dramatic reduction in the extended X-ray absorption fine structure (EXAFS) intensity. In contrast, Ga/γ-Al 2 O 3 does not show such behavior. In order to lay a rigorous foundation for characterizing these types of systems and to establish experimental signatures for the elusive Ga(I) oxidation state, we recorded Ga K-edge X-ray absorption spectra [including high-energy-resolution fluorescence detection-X-ray absorption near-edge spectroscopy (HERFD-XANES)] for several well-defined molecular and crystalline Ga(I) compounds. XANES is essential to establishing the presence of Ga(I), despite the overlap in edge energies with organoGa(III) compounds, because Ga(I)-containing oxide materials show very weak EXAFS scattering. Compared to the XANES of trigonal Ga(III)-containing materials, Ga(I) spectra display a significantly more intense white line feature. Theoretical simulations agree well with this experimental observation and reveal that the strong XANES intensity originates from the superposition of transitions to several empty, nearly degenerate p-like states. Finally, these signatures provide compelling evidence for assigning the intense white line and dramatic loss of EXAFS intensity in Ga/HZSM-5 to the near-quantitative reduction of Ga(III) to Ga(I), while the weaker white line and conventional EXAFS signal of Ga/γ-Al 2 O 3 point to, at most, a minor fraction of Ga(I) sites.

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

Tip‐Enhanced Imaging and Control of Infrared Strong Light‐Matter Interaction

Optical antenna resonators enable control of light‐matter interactions on the nano‐scale via electron–photon hybrid states in strong coupling. Specifically, mid‐infrared (MIR) nano‐antennas coupled to saturable intersubband transitions in multi‐quantum‐well (MQW) semiconductor heterostructures allow for the coupling strength to be tuned through antenna resonance and field intensity. Here, in this study, tip‐enhanced nano‐scale variation of antenna‐MQW coupling across the antenna is demonstrated, with a spatially‐dependent coupling strength $g_{\textrm{aq}}$ varying from 73 (strong coupling) to 24 cm -1 (weak coupling). This behavior is modeled based on the spatially dependent local constructive and destructive interference between tip and antenna fields. Using a quantum‐mechanical density‐matrix model of the MQW system with its designed values of transition dipole moment, doping density, and population decay time, the picosecond IR pulse coupling to intersubband transitions and the associated tip induced strong‐field saturation effects are described. These results present a new regime of nonlinear IR light‐matter control based on the dynamic manipulation of quantum hybrid states on the nanoscale and in the infrared, with a perspective regarding extension to molecular vibrations.

Wang, Yueying↗
DOE OSTI2024

Mechanically activated and deactivated ion transport across nanopores with heterogeneous surface charge distributions

To mimic the intricate and adaptive functionalities of biological ion channels, electrohydrodynamic ion transport has been studied extensively, albeit mostly, across uniformly charged nanochannels. Here, we analyze the ion transport under coupled electric field and pressure across heterogeneously charged nanopores with oppositely charged sections on their lateral surface. We only consider such pores with symmetric hourglass-like and cylindrical shapes to focus on the effects of the non-uniform surface charge distribution. Finite-element simulations of a continuum model demonstrate that a pressure applied in either direction of the pore-axis equally suppresses or amplifies the ionic conductance, depending on the electric field polarity, by distorting the quasi-static distribution of ions in the pore. The resulting anomalous mechanical deactivation and activation of ionic current under opposite voltage biases exhibit the functional modularity of our setup, while their intensities are highly tunable, substantially greater than those of analogous behaviors in other nanochannels, and fundamentally correlated to ionic current rectification (ICR) in our pores. A detailed study of ICR subsequently reveals its counterintuitive non-monotonous variations, in the pores, with the magnitude of applied voltage and the pore length, that can help optimize their diode-like behavior. We further illustrate that while the hourglass-shaped nanopores yield the more efficient mechanical suppressors of ion transport, their cylindrical analogs are the superior rectifiers and mechanical amplifiers of ion conduction. Therefore, this article provides a blueprint for the strategic design of nanofluidic circuits to attain a robust, modular, and tunable control of ion transport under external electrical and mechanical stimuli.

Physics↗
DOE OSTI2026

The role of focusing geometry in MeV x-ray production from petawatt laser–solid interaction

Relativistic laser–plasma interactions provide a compact and flexible route to generating bright, ultrashort pulses of MeV x rays, with applications in high-energy density science, nuclear physics, and radiography. Despite extensive study of intensity scaling in laser-driven electron acceleration, the role of focusing geometry and focal-volume effects in MeV radiation production remains insufficiently understood. Here, we present experimental and kinetic simulation results from the Texas Petawatt Laser (120 J, 140 fs) in which the focusing geometry (f/3 or f/1.5) is varied, while the laser energy and pulse duration are fixed. Experimentally, the f/3 geometry produces approximately three times more MeV radiation than the f/1.5 geometry, despite its twofold lower nominal vacuum intensity. This result is based on the time-integrated radiation per unit solid angle, as measured along the diagnostic line of sight. Three-dimensional particle-in-cell simulations reproduce this trend when a modest (10s of μm) effective focal plane shift is introduced, demonstrating that relativistic laser–plasma coupling is highly sensitive to focal geometry in the presence of a preplasma. This behavior is consistent with differences in interaction length and effective intensity at the critical surface between the two configurations. The sensitivity of the tightly focused f/1.5 configuration ($z_R ≈$⁠ 5 μ m) reflects the combined influence of thermal lensing and an extended preplasma, while the f/3 geometry (⁠$z_R ≈$ 80 μ m) remains comparatively robust. These results demonstrate a breakdown of conventional intensity scaling and identify focusing geometry as a critical control parameter for MeV electron and x-ray generation at petawatt powers. Implications for laser–plasma-based radiographic facilities are also discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Molecular concentration field design using closed-form steady-state solutions

Control over spatial concentration fields represents a fundamental challenge in designing synthetic biological systems and programmable soft materials. While nature creates morphogen gradients that orchestrate complex developmental processes, synthetic approaches have largely relied on empirical optimization and computationally intensive simulations. Here, we present an analytical framework for steady-state concentration fields generated by finite-sized localized sources in diffusion–degradation systems and derive closed-form solutions for one-, two-, and three-dimensional geometries. By expressing these solutions in dimensionless form, we show that gradient steepness and spatial structure are organized by the Thiele modulus, which captures the competition between diffusion and degradation length scales. The analysis reveals distinct design regimes: in degradation-dominated systems, gradient shape is governed by exponential decay and becomes dimension-independent, whereas in diffusion-dominated systems, gradient magnitude and extent follow dimension-dependent power-law scaling. Building on these results, we introduce a quantitative design strategy that uses threshold-based criteria to program concentration ranges by tuning physically accessible parameters, most directly the production rate, while holding transport and degradation properties fixed. Comparisons with numerical solutions and reported experimental systems demonstrate consistency with the predicted scaling behavior. Together, this work provides a generalizable and physically transparent framework for designing steady-state concentration fields in synthetic biological and soft matter systems, enabling predictive control of gradient-mediated organization without reliance on extensive numerical optimization.

Kim, Dong Woo [Johns Hopkins University, Baltimore↗
DOE OSTI2024

Metal DiSulfide thin films synthesized by pulsed ion beam ablation: Self-assembled structures with improved mechanical properties in humid conditions

Thin films composed of transition metal dichalcogenides such as molybdenum disulfide (MoS 2 ) or tungsten disulfide (WS 2 ) are well known for their solid lubricating behavior in dry environments such as space, but deteriorate rapidly in environments where oxygen is present. There have been many attempts over recent decades to improve MoS 2 lubricating performance in humid environments by addition of metals such as titanium to the film. The additions have come in the form of co-sputtering with the MoS 2 , or formation of nano multilayer coatings. We have previously reported on the synthesis of MoS 2 films by high-power ion beam ablation [Prasad et al, Matls. Lett. 65 (2010) 4–6]. The resulting films are observed to contain spheres of pure Mo of 10–100 nm size evenly dispersed throughout the ~1 μm-thick film, a form of self-assembly not before seen in films produced by sputtering or by pulsed laser deposition (PLD). The remaining Mosingle bondS matrix was observed to crystalize in situ under frictional contact into MoS 2 with the basal planes orienting themselves in the direction of sliding underneath the wear surface. The present study concerns metal additions in the form of alternating layers with MoS 2 deposited in between from an ablation target. The addition of Ti resulted in the extension of good tribological performance to humid air conditions, but the path to good performance appears to be relatively narrow. The same distributed Mo spheres were present in the new film, and similar self-assembled morphologies when V was substituted for the Ti, and when multilayers were formed of WS 2 with Zr additions. Furthermore, it appears that addition of Ti resulted in improvements in mechanical film strength that contributed to the wear ability in humid conditions.

Intense ion beams↗