SEARCH · Search NASA

Results for “GROUP DYNAMICS”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 37 records · Page 2

All sources NASA NTRS DOE OSTI

DOE OSTI2026

Hyporheic‐Zone Processes and Stream Oxygen Dynamics: Insights From a Multiscale Reactive Transport Model

Aquatic ecosystem metabolism encapsulates the daily fixation (gross primary production, GPP d ) and mineralization (ecosystem respiration, ER d ) of organic carbon. In fluvial systems, these are commonly estimated by inverse solutions to field observations using a model that describes oxygen concentrations varying in the water column in response to metabolic fluxes and air‐water gas exchange controlled by a rate coefficient (K 600 ). The most common conceptual model is the single‐station metabolism (SSM) model. The simplicity and flexibility of this conceptualization make it attractive; however, it implicitly assumes that all the processes that consume oxygen in fluvial systems can be lumped into a bulk estimate of respiration with poorly understood consequences for estimates of GPP d , ER d , and K 600 . Here, we focus on the implications of using SSM conceptualization when estimating metabolic fluxes from oxygen dynamics in channels where hyporheic exchange occurs. We use a new multiscale numerical model for reactive transport in streams that represents hyporheic exchange and streambed heterotrophic respiration. Nondimensionalization of this model reveals dimensionless groups that collectively control oxygen dynamics. Numerical experiments offer a mechanistic understanding of the impacts of hyporheic exchange on diel oxygen dynamics revealing that potential biases arise from neglecting mass transfer limitations. Specifically, we found that hyporheic exchange significantly affects diel oxygen dynamics, even for nonreactive streambed sediments. Moreover, while the SSM performs well in many situations, we find conditions where significant bias is produced by hyporheic exchange, even when oxygen data are well‐fitted. These situations pose a major challenge in the interpretation of metabolism assessment estimates.

Gomez‐Velez, Jesus D. [Oak Ridge National Laborato↗

DOE OSTI2026

Strain-Driven Evolution of Structural and Carrier Recombination Dynamics in Germanium

Strain engineering in group-IV semiconductors, in particular germanium (Ge), offers a promising route for advancing next-generation optoelectronic device performance. High minority carrier lifetime serves as a direct indicator of superior material quality and device efficiency. Several factors govern the minority carrier lifetime in both strained Ge and strained GeSn, such as strain level, epilayer thickness, defects, surface roughness, and Sn alloying. In this study, the impact of biaxial strain on effective minority carrier lifetimes of tensile-strained Ge (0% ..epsilon..-Ge/Al?As, 0.78% ..epsilon..-Ge/In0.115?Ga0.885?As, and 0.93% ..epsilon..-Ge/In0.16?Ga0.84?As) and compressively strained GeSn (-0.53% c-Ge0.97?Sn0.03/Al?As) epilayers, grown by molecular beam epitaxy, was experimentally investigated via the microwave-reflection photoconductive decay (..mu..-PCD) technique as a simultaneous function of thickness and strain. The structural and strain-relaxation properties were analyzed by high-resolution x-ray diffraction and Raman spectroscopy, while surface morphology analysis using atomic force microscopy provided insights into the surface roughness of strained Ge-based epilayers. Defect characterization by cross-section transmission electron microscopy revealed superior crystalline quality of ..epsilon..-Ge in one epilayer, while twin boundary formation was observed in the other. We investigated the carrier recombination dynamics in epitaxially strained Ge by extracting the bulk recombination lifetime and surface recombination velocity using the thickness-dependent minority carrier lifetime. High effective minority carrier lifetimes of approximately 60-340 ns were achieved at high strain levels with minimal dependence on the thickness of the strained Ge-based epilayers by means of the ..mu..-PCD technique. Minority carrier lifetimes in Ge-based epilayers exhibit a quadratic dependence on the biaxial strain, ranging from -0.88% (compressive) to +1.6% (tensile). The analysis of the carrier recombination dynamics in strained Ge-based heterostructures offers a pathway to designing group-IV based optoelectronic devices.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Tough and circular glass fiber composites via a tailored dynamic boronic ester interface

Glass fiber reinforced polymer (GFRP) composites are valued for their strength and cost-effectiveness. However, traditional GFRPs often face challenges for end-of-life recycling due to their non-depolymerizable thermoset matrices, and long-term performance due to inadequate interfacial adhesion, which can lead to fiber–matrix delamination. Here, in this work, we have designed dynamic fiber–matrix interfaces to allow tough and closed-loop recyclable GFRPs by utilizing a vitrimer, derived from upcycled polystyrene-b-poly(ethylene-co-butylene)-b-polystyrene (SEBS) with boronic ester (S-Bpin) and amine-based diol crosslinker. The boronic ester groups in S-Bpin form dynamic covalent bonds with the naturally present hydroxyl groups on the unsized GF surface, which eliminates the need for fiber sizing and enables facile closed-loop recyclability of both the fibers and the vitrimer matrix. The resulting strong fiber–matrix interface, depicted by the Raman mapping, leads to a 552% increase in-plane shear toughness (6.2 ± 0.3 MJ m -3 ) and 27% ultimate tensile strength (361 ± 89.2 MPa) compared to those of the conventional epoxy-based matrix (0.95 ± 0.05 MJ m -3 and 264 ± 59.7 MPa, respectively). The network rearrangement through dynamic boronic ester exchange enables fast thermoformability and repairability of micro-cracks at elevated temperatures. Additionally, both the matrix and composite demonstrate strong adhesion to various surfaces including steel and glasses exhibiting ≥6 MPa lap shear strength, which expands their suitability for diverse industrial applications. The readily created dynamic interface between boronic ester functionalized vitrimer and neat GFs presents a promising strategy for developing closed-loop recyclable, multifunctional structural materials, offering a sustainable alternative to non-recyclable thermoset GFRPs and contributes to a circular economy in composite materials.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Quantifying the Resolution Sensitivity of the Kain–Fritsch Scheme Across the Gray Zone by Isolating Interactions: A TWP‐ICE Case Study

The resolution sensitivity of the Kain–Fritsch (KF) convection scheme and the role of interactions between the physics and dynamics within the gray zone (<10 km) were investigated using the Separate Physics and Dynamics Experiment (SPADE) framework. Two groups of experiments were conducted using the Weather Research and Forecasting (WRF) model via traditional (Tradition) runs and SPADE runs with resolutions of 1, 2, 4, and 8 km during the wet period of the Tropical Warm Pool–International Cloud Experiment (TWP‐ICE). Results show that the KF scheme simulates the weakened convective processes well as the resolution increases in both groups, and the changes in the convective variables with resolution in SPADE are smaller than in the Tradition group. This indicates the important effects of interactions between model components on convection parameterizations as the resolution changes. Additionally, the microphysics variables remain nearly unchanged with resolution in SPADE and weaken slightly in Tradition as the resolution decreases, suggesting the relatively weaker influences of model interactions for the resolved‐cloud parameterization. Therefore, the scale‐aware behavior of KF scheme is further strengthened in Tradition runs, primarily through inhibiting the strength of stratiform processes through physics–dynamics interactions and physical components.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Influence of substitution pattern on the dynamics of internal conversion and intersystem crossing in thiopyridone isomers

We report a combined experimental and theoretical investigation of the ultrafast internal conversion (IC) and intersystem crossing (ISC) dynamics of two thiopyridone (TP) isomers in solution. Our study used ultrafast transient X-ray absorption spectroscopy (XAS) at the sulfur K-edge, in conjunction with electronic excited state surface hopping molecular dynamics and simulations of the excited state XAS, to investigate the impact of the functional group substitution pattern and solvent on the dynamics of IC and ISC. The combination of the localized X-ray probe and the simulation results enables, in part, the differentiation between ππ* and nπ* character excited states, as well as singlet and triplet states. Access to nπ* character excitations has particular value since they often prove challenging to assess with optical spectroscopy. For 2-TP, the photoexcited S 2 (ππ*) state rapidly undergoes IC to the S 1 (nπ*) state below the instrument response time, followed by ISC to the T 1 (ππ*) state on a timescale of 600 fs in acetonitrile. For 4-TP, the timescale of S 2 to S 1 IC increases to 330 fs and the timescale of ISC increases to more than 10 ps. The differences between isomers are rationalized by considering the key role of the, nπ* intermediates in mediating the intersystem crossing of these systems. Varying the substitution pattern of the molecule can stabilize or destabilize these intermediates leading to the increase in ISC rate in the ortho isomer as compared to the para isomer, while changing the solvent from acetonitrile to water had minimal effect on the electronic excited state relaxation mechanism.

Garratt, Douglas [SLAC National Accelerator Labora↗

DOE OSTI2026

Charge Dependence of Local Hydration Dynamics in Poly(Acrylic Acid) Solutions

Localized measurements of hydration dynamics surrounding poly(acrylic acid) (pAA) chains are executed using Overhauser dynamic nuclear polarization relaxometry to determine the role of polymer protonation state on the polymer solvation environment. Alterations in polymer concentration and solution pH uncovered three different regimes of local water dynamics ranging from bulk-water-like to subdued water diffusivity. These experiments, combined with molecular dynamics simulations, reveal that the presence of deprotonated carboxylic acid groups promotes a hydrophilic local environment and extended polymer chain configurations that together incentivize the incorporation of bulk-like water near the polymer. This results in a stretched polymer conformation that leads to faster observed hydration dynamics in the local hydration environment. Meanwhile, protonated carboxylic acid groups are expected to engage in intrapolymer hydrogen bonding between monomer groups, resulting in a more collapsed conformation and slower observed hydration dynamics. The importance of the carboxylic acid groups on mediating hydration behavior is further illustrated by their role in dictating the solvation environment around poly(acrylic acid-stat-(poly(ethylene glycol) methyl ether acrylate) (p(AA-stat-PEGMEA)) copolymers with varying monomer ratios. Furthermore, these local hydration measurements establish a connection between local water dynamics and polymer chain conformation and may provide additional molecular perspective into poly(acrylic acid)’s classification as a superabsorbent polymer.

Conformation↗

DOE OSTI2026

Semi-Analytical Hierarchical Bayesian Inference of Nonlinear Model Structure in Stochastic Dynamics: Applied to Compartmental Models of Infectious Diseases

A Bayesian computational framework for parsimonious inference in stochastic nonlinear dynamical systems is presented. This framework enables the concurrent estimation of system states, time-varying parameters, time-invariant parameters, and the optimal sparsity structure of the model parameters. Because differential equation-based models are often simplified mechanistic or phenomenological representations, robust inference from noisy measurement data requires explicit treatment of model error and uncertainty. Model error and time-varying parameters can be represented as random processes, enabling inference while making minimal assumptions about the underlying sources of discrepancy and variability. Adopting stochastic differential equation representations affords the model significant flexibility, but can also render it susceptible to overfitting during statistical inversion, where the inferred model may track noise rather than the underlying signal. To alleviate the effects of overfitting and to enable the discovery of the optimal sparse representation of the time-invariant parameters, a Bayesian sparse learning algorithm is embedded within the framework. This sparse learning framework adopts an approximate hierarchical Bayesian setting defined by a series of semi-analytical expressions. The model structure inference framework is validated using a stochastic compartmental model for tracking and forecasting active cases of an infectious disease. Compartmental models describe population-level infectious disease dynamics through interactions among population fractions grouped by disease state. Mathematically, such models consist of a system of coupled ordinary differential equations. This example adopts an expressive compartmental model that includes multiple possible interactions between disease states, motivated by early uncertainty surrounding COVID-19 reinfection dynamics and their implications for long-term epidemic forecasting. The sparse learning exercise permits the inference of a priori unknown epidemiological dynamics from simulated public health data, discovering the nested compartmental model that optimizes the trade-off between average data-fit and model complexity. It is shown that inducing sparsity among the model parameters eliminates redundant interactions between compartments, equivalently revealing the optimal coupling structure between differential equations.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2025

Introducing Kynema, an Open-Source Performance-Portable Flexible-Multibody-Dynamics Solver

In this talk we introduce Kynema, an open-source general flexible-multibody-dynamics solver that is well suited for simulating wind turbine structural dynamics. Kynema uses a Lie-group time integrator for constrained systems and runs on both CPUs and GPUs. Timing results for simulations are presented for the IEA 15-MW turbine with and without aerodynamic forces.

17 WIND ENERGY↗

DOE OSTI2024

A spatially resolved spectral analysis of giant radio galaxies with MeerKAT

ABSTRACT In this study we report the spatially resolved, wideband spectral properties of three giant radio galaxies (GRGs) in the COSMOS field: MGTC J095959.63+024608.6, MGTC J100016.84+015133.0, and MGTC J100022.85+031520.4. One of these galaxies, MGTC J100022.85+031520.4, is reported here for the first time, with a projected linear size of 1.29 Mpc at a redshift of 0.1034. Unlike the other two, it is associated with a brightest cluster galaxy (BCG), making it one of the few GRGs known to inhabit cluster environments. We examine the spectral age distributions of the three GRGs using new MeerKAT UHF-band (544–1088 MHz) observations, and L-band (900–1670 MHz) data from the MeerKAT International GHz Tiered Extragalactic Exploration (MIGHTEE) survey. We test two models of spectral ageing, the Jaffe–Perola and Tribble models, using the Broadband Radio Astronomy Tools (brats) software, and find that they agree well with each other. We estimate the Tribble spectral age for MGTC J095959.63+024608.6 as 68 Myr, for MGTC J100016.84+015133.0 as 47 Myr, and for MGTC J100022.85+031520.4 as 67 Myr. We find significant disagreements between these spectral age estimates and the estimates of the dynamical ages of these GRGs, modelled in cluster and group environments. Our results highlight the need for additional processes that are not accounted for in either the dynamic age or the spectral age estimations.

Charlton, K. K. L. (ORCID:0000000229252047)↗

DOE OSTI2025

Preliminary modeling of triply periodic minimal surface (TPMS) structures using RELAP5-3D

With the United States Department of Energy (DOE)’s goal of quadrupling the nation’s nuclear energy supply by 2050, and with the Advanced Fuels Campaign pushing for new types of advanced reactor fuels and geometries, the need has arisen for new nuclear fuel designs. One such design is to swap out current nuclear fuel geometries in exchange for another type of geometry, called a Triply Periodic Minimal Surface (TPMS). TPMSs are self-supporting, infinitely repeating lattices—attributes that lend themselves well to additive manufacturing. These surfaces also possess enhanced heat transfer properties thanks to their internal area changes and large surface-area-to-volume ratios. Their drawback, however, is an increased pressure drop. Given the small amount of correlations and data (Reynolds numbers in the 2,000–8,000 range), and the minimal amount of experience so far obtained by modeling TPMS structures using 1D systems codes such as the Reactor Excursion and Leak Analysis Program (RELAP5-3D), further research into this topic was needed. Using data from the University of Wisconsin - Madison (UW), curve fits were created for both a Heat Transfer Coefficient (HTC) correlation and a Darcy friction factor empirical coefficient correlation. The curves’ coefficients and multipliers were then output and utilized in RELAP5-3D models of two upcoming experiments—Flow Loop for INFLUX Pressure drop (FLIP) and Microreactor Agile Non-nuclear Experimental Test (MAGNET)—aimed at increasing the available data for Reynolds numbers to the 16,000–36,000 range for TPMS structures. The models were run under the conditions utilized by a Computational Fluid Dynamics (CFD) analysis performed by another group at Idaho National Laboratory. Only CFD pressure drop values were obtained from the FLIP test, and those values showed that the RELAP5-3D models had a lower rate of pressure increase in comparison to the CFD values. In addition, there seemed to be a vertical shift upward in the pressure drop for both models whenever the TPMS porosity decreased, and the RELAP5-3D models showed a higher vertical shift in comparison to the CFD values. The MAGNET results did not correspond to any CFD or experimental results against which they could be compared, so they were instead compared against the proposed CFD input conditions. These values were then compared with each other to make sure the model seemed to be performing as expected, paving the way for future tests that can be run for the purpose of further analyses and comparisons. The pressure drop increased with temperature and mass flow rate independently. The temperature change would decrease with increasing mass flow rate and temperature, which was just as we expected based on the fact that the lower viscosity and decreased density would result in higher friction and churning losses. The last metric that was assessed was the enthalpy flow change, which increased with increasing mass flow rate and decreasing temperature.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2025

Engineering micromotion in Floquet prethermalization via space-time symmetries

We present a systematic framework for Floquet prethermalization under strong resonant driving, emphasizing the pivotal role of dynamical space-time symmetries. Our approach demonstrates how dynamical space-time symmetries map onto the projective static symmetry group of the prethermal Hamiltonian governing the prethermal regime. We introduce techniques for detecting dynamical symmetries through the time evolution of local observables, facilitating a detailed analysis of micromotion within each period and surpassing the limitations of conventional stroboscopic Floquet prethermal dynamics. To implement this framework, we present a prethermal protocol that preserves order-2 dynamical symmetry in a spin-ladder model, confirming the predicted relationships between the expectation values of local observables at distinct temporal points in the Floquet cycle, linked by this symmetry.

floquet systems↗

DOE OSTI2026

The Future Polarized Target Program at Jefferson Lab

Polarized targets have played a crucial role in Jefferson Lab's exploration of nuclear structure over the past four decades. The three original experimental halls have seen 19 separate installations of polarized solid or gas targets for use in the particle physics scattering experiments, and this trend will continue in the next decade. Five polarized target systems are in preparation for use at JLab in the coming years. Hall B will see the use of two solid polarized targets, one longitudinally polarized to the beam, the other transversely, as well as a novel 3He gas polarized target. In Hall C, new experiments will augment the tensor polarization in dynamically polarized solids. Plans are under development to bring a polarized solid target to Jefferson Lab's photon beam hall, Hall D, for the first time. To support these efforts, the JLab polarized target group is building a test laboratory to develop dynamic nuclear polarization techniques, as well as an apparatus to irradiate target material using electrons from JLab's injector test facility. In this talk, we will explore the development progress and plans for each of these efforts.

Maxwell, James [Thomas Jefferson National Accelera↗

DOE OSTI2025

Attosecond X-Ray Core-Level Chronoscopy of Aromatic Molecules

Attosecond photoemission or photoionization delays are a unique probe of the structure and the electronic dynamics of matter. However, the spectral congestion of valence photoelectron spectra sets fundamental limits to the complexity of systems that can be studied, and the delocalization of valence electron wave functions blurs the spatial origin of the photoelectron wave packet. Using attosecond x-ray pulses from LCLS, we demonstrate the key advantages of measuring core-level delays: The photoelectron spectra remain atomlike, the measurements become element specific, and the observed scattering dynamics originate from a pointlike source when multicenter interference effects are negligible. We exploit these unique features to reveal the effects of changing functional groups (C-H vs N) and symmetry on attosecond scattering dynamics by measuring and calculating the photoionization delays between N−1⁢𝑠 and C−1⁢𝑠 core shells of a series of aromatic azabenzene molecules. Remarkably, the delays increase with the number of nitrogen atoms in the molecule and reveal multiple resonances. We identify two previously unknown mechanisms regulating the associated attosecond dynamics, namely the enhanced confinement of the trapped wave function with the replacement of C-H groups by N atoms and the decrease of the coupling strength among the photoemitted partial waves with increasing symmetry. This study demonstrates the unique opportunities opened by measurements of core-level photoionization delays for unraveling attosecond electron dynamics in complex matter.

Ji, Jia-Bao [Eidgenoessische Technische Hochschule↗

DOE OSTI2025

Characterizing Hydrated Polymers via Dielectric Relaxation Spectroscopy: Connecting Relative Permittivity, State of Water, and Salt Transport Properties of Sulfonated Polysulfones

Sulfonated polysulfone is a promising membrane material for separation and energy generation processes that rely on membranes to control the rates of small-molecule (e.g., water and ions) transport. The interactions among water molecules, ions, and the sulfonate groups in these polymers play a key role in controlling these rates of transport, but much remains unknown about these fundamental interactions in sulfonated polymers. In this study, we used dielectric relaxation spectroscopy to characterize water molecule dynamics in sulfonated polysulfone and Nafion. We found that the charged sulfonate groups contribute to a restriction of water molecule dynamics (i.e., a reduction in the characteristic time scale of dipolar motions) in a manner that is governed by the concentration and nature (i.e., conjugate base strength) of the sulfonate group. Additionally, we develop strategies to use these data to aid in modeling ion transport in sulfonated polysulfone. These results may be useful to guide engineering strategies for polymeric membranes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Interpreting the spatial distribution of soil properties with a physically-based distributed hydrological model

Digital soil maps are commonly data-driven as the development of physically-based models for soil mapping is difficult due to the complexity of soils. However, physically-based hydrologic models have been successful in simulating water dynamics. Since water movement is a major driver of pedogenesis, the physical rules that govern water movement might help explain and predict the spatial variation of soil properties. Here, we demonstrate the novel use of a physically-based, distributed hydrologic model to inform the spatial distribution of soil properties. The Distributed Hydrology Soil Vegetation Model (DHSVM) was utilized to simulate soil moisture content (SM) and water table depth (WTD) in two hillslope catchments under pasture and forest management wherein hydrologic model outputs were then compared with soil properties measured in situ. SM sensors and wells were installed in both catchments to validate simulations of soil water movement via Nash-Sutcliffe Efficiency (E). In-situ observations were made at 87 sites within both catchments to study the connection between simulated water movement (SM and WTD) and observed soil properties, namely the depth and thickness of the argillic (Bt), fragic (Btx), and C horizons, and the depth of redoximorphic features. The simulated time series of SM and WTD were also clustered per season using Dynamic Time Warping (DTW), which identified similarity among time series at varying timescales. Model validation suggested that simulations of surficial SM (0–20 cm) were reasonable (E = 0.45), however, simulated subsurface SM (45–60 cm) and WTD were not sufficiently accurate. The thickness of Btx horizons were spatially grouped into different populations by SM clusters from every season except spring. For the other properties, only SM dynamics of specific seasons grouped into significantly different populations, suggesting that the explanatory power of simulated water movement varies seasonally and was greater during winter. Here, we show clusters of simulated SM separated soil properties into statistically different populations, showing that hydrologic models could inform areas that followed different water dynamics related to pedogenic trajectories and related biogeochemical processes not necessarily simulated by the model. As such, physically-based modeling of water dynamics can, therefore, inform and advance digital soil mapping by linking water movement patterns stemming from hydrologic model outputs to spatial patterns of soil properties and pedogenesis.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

The Fluid Dynamics Uncertainty Quantification Challenge Problem: XFOIL vs. MFOIL

Uncertainty quantification (UQ) has become more critical in aerospace engineering due to the growing dependence on computational tools for design optimization and performance analyses of aerospace vehicles. Even though the significance of UQ in assessing the credibility of computational analyses is well recognized, its costs and complexity impede its integration into standard practices, particularly in computational fluid dynamics (CFD) and other fluid analyses. This paper presents a UQ study for low-fidelity computational aerodynamics analyses with XFOIL and mfoil (i.e., the MATLAB version of XFOIL with several implementation modifications); these tools are utilized widely in both research and education. The main contributions of this paper are as follows: 1) improved precision in quantifying the uncertainty of the baseline Monte Carlo results used to benchmark surrogate modeling techniques for UQ, 2) quantification of the effect of the implementation differences between XFOIL and mfoil on solution quantities of interest (QoIs), such as lift and pitching moment coefficients, and 3) development of an open-source UQ library for use with XFOIL and mfoil, which has educational values and helps promote UQ for fluid analyses with aerospace applications. Results and discussions revolve around cases 1-4 of the challenge problem posed by the AIAA Fluid Dynamics Technical Committee’s Uncertainty Quantification Discussion Group (UQDG). In case 3, this work employs CFDverify, an open-source solution verification software, to quantify the discretization error and evaluate the extrapolated QoIs based on the grid convergence index (GCI). This UQ study differentiates itself from previous studies in the rigor of handling baseline Monte Carlo uncertainty and in including mfoil, which is a more accessible alternative to XFOIL. Finally, despite the growing computing power, low-fidelity computational tools remain valuable, such as for aerodynamic shape optimization at Mach numbers below 0.65 and low-to-mid Reynolds numbers.

Lay, Aidan S [University of Tennessee, Knoxville (↗

DOE OSTI2025

Towards excitations and dynamical quantities in correlated lattices with density matrix embedding theory

Density matrix embedding theory (DMET) provides a framework to describe ground-state expectation values in strongly correlated systems, but its extension to dynamical quantities is still an open problem. We show one route to obtaining excitations and dynamical spectral functions by using the techniques of DMET to approximate the matrix elements that arise in a single-mode inspired excitation ansatz. We demonstrate this approach in the one-dimensional Hubbard model, comparing the neutral excitations, single-particle density of states, charge, and spin dynamical structure factors to benchmarks from the Bethe ansatz and density matrix renormalization group. Finally, our work highlights the potential of these ideas in building computationally efficient approaches for dynamical quantities.

1-dimensional systems↗

DOE OSTI2024

Computational multiphysics modeling of radioactive aerosol deposition in diverse human respiratory tract geometries

The evaluation of aerosol exposure relies on generic mathematical models that assume uniform particle deposition profiles over the human respiratory tract and do not account for subject-specific characteristics. Here we introduce a hybrid-automated computational workflow that generates personalized particle deposition profiles in 3D reconstructed human airways from computed tomography scans using Computational Fluid and Particle Dynamics simulations. This is the first large-scale study to consider realistic airways variability, where 380 lower and 40 upper human respiratory tract 3D geometries are reconstructed and parameterized. The data is clustered into nine groups using random forest regression. Computational fluid and particle dynamics simulations are conducted on these representative geometries using a realistic heavy-breathing respiratory cycle and radioactive iodine-131 as a source term. Monte Carlo radiation transport simulations are performed to obtain detailed energy deposition maps. Our findings emphasize the importance of personalized studies, as minor respiratory tract variations notably influence deposition patterns rather than global parameters of the lower airways, observing more than 30% variance in the mass deposition fraction.

62 RADIOLOGY AND NUCLEAR MEDICINE↗