Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Dynamic range”

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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Design of nematic liquid crystals to control microscale dynamics

The dynamics of small particles, both living such as swimming bacteria and inanimate, such as colloidal spheres, has fascinated scientists for centuries. If one could learn how to control and streamline their chaotic motion, that would open technological opportunities in the transformation of stored or environmental energy into systematic motion, with applications in micro-robotics, transport of matter, guided morphogenesis. This review presents an approach to command microscale dynamics by replacing an isotropic medium with a liquid crystal. Orientational order and associated properties, such as elasticity, surface anchoring, and bulk anisotropy, enable new dynamic effects, ranging from the appearance and propagation of particle-like solitary waves to self-locomotion of an active droplet. By using photoalignment, the liquid crystal can be patterned into predesigned structures. In the presence of the electric field, these patterns enable the transport of solid and fluid particles through nonlinear electrokinetics rooted in anisotropy of conductivity and permittivity. Director patterns command the dynamics of swimming bacteria, guiding their trajectories, polarity of swimming, and distribution in space. This guidance is of a higher level of complexity than a simple following of the director by rod-like microorganisms. Namely, the director gradients mediate hydrodynamic interactions of bacteria to produce an active force and collective polar modes of swimming. The patterned director could also be engraved in a liquid crystal elastomer. When an elastomer coating is activated by heat or light, these patterns produce a deterministic surface topography. The director gradients define an activation force that shapes the elastomer in a manner similar to the active stresses triggering flows in active nematics. Finally, the patterned elastomer substrates could be used to define the orientation of cells in living tissues. The liquid-crystal guidance holds a major promise in achieving the goal of commanding microscale active flows.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Examining the vertical heterogeneity of aerosols over the Southern Great Plains

Abstract. Atmospheric aerosols affect the global energy budget by scattering and absorbing sunlight (direct effects) and by changing the microphysical structure, lifetime, and coverage of clouds (indirect effects). Both aerosol direct and indirect effects are affected by the vertical distribution of aerosols in the atmosphere, which is further influenced by a range of processes, such as aerosol dynamics, long-range transport, and entrainment. However, many observations of these processes are based on ground measurements, limiting our ability to understand the vertical distribution of aerosols and simulate their impact on clouds and climate. In this work, we examined the vertical heterogeneity of aerosols over the Atmospheric Radiation Measurement (ARM) Southern Great Plains (SGP) using data collected from the Holistic Interactions of Shallow Clouds, Aerosols and Land Ecosystems (HI-SCALE) campaign. The vertical profiles of meteorological and aerosol physiochemical properties up to 2500 m above are examined based on the 38 flights conducted during the HI-SCALE campaign. The aerosol properties over the SGP show strong vertical heterogeneity and seasonal variabilities. The aerosol concentrations at the surface are the highest due to strong emissions at ground level. In general, the mode diameter of aerosols during summer (∼ 100 nm) is larger than that during spring (∼ 30 nm), a result of enhanced condensational growth due to enriched volatile organic compounds in summer. The concentration of aerosols below 30 nm in the boundary layer (BL) (e.g., below 1000 m) during spring is higher than that during summer, a result of the stronger new particle formation (NPF) events due to the reduced condensation sink in spring. In the BL, the size of the aerosols gradually increases with altitude due to condensational growth and cloud processing. However, the chemical composition of the aerosols remained similar, with organics and sulfates representing 59.8 ± 2.2 % and 22.7 ± 2.1 %, respectively, of the total mass in the BL. Through the vertical profiles of aerosol properties, we observed NPF events in the upper BL during 7 out of 38 research flights, where the newly formed particles continue to grow as they are mixed down to the surface. There is also an indication that deep convection brings aerosols from the free troposphere (FT) to the surface, where they grow to contribute to the cloud condensation nuclei (CCN). Overall, the vertical heterogeneity of aerosols over the SGP is influenced by aerosol dynamics (new particle formation, growth, and cloud processing) and transport processes (mixing in the BL, long-range transport, entrainment, and convective downward transport). Case studies showing the influence of these factors are discussed.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

The Solar Neighborhood. XLVII. Comparing M-dwarf Models with Hubble Space Telescope Dynamical Masses and Spectroscopy

We use HST/STIS optical spectroscopy of 10 M dwarfs in five closely separated binary systems to test models of M-dwarf structure and evolution. Individual dynamical masses ranging from 0.083 to 0.405 M {sub ⊙} for all stars are known from previous work. We first derive temperature, radius, luminosity, surface gravity, and metallicity by fitting the BT-Settl atmospheric models. We verify that our methodology agrees with empirical results from long-baseline optical interferometry for stars of similar spectral types. We then test whether or not evolutionary models can predict those quantities given the stars’ known dynamical masses and the conditions of coevality and equal metallicity within each binary system. We apply this test to five different evolutionary model sets: the Dartmouth models, the MESA/MIST models, the models of Baraffe et al., the PARSEC models, and the YaPSI models. We find marginal agreement between evolutionary model predictions and observations, with few cases where the models respect the condition of coevality in a self-consistent manner. We discuss the pros and cons of each family of models and compare their predictive power.

47 OTHER INSTRUMENTATION↗
DOE OSTI2024

IRIS: Exploring Performance Scaling of the Intelligent Runtime System and its Dynamic Scheduling Policies

High-Performance Computing is becoming increasingly heterogeneous, relying on a diverse mix of hardware to achieve good performance. Paradoxically, current drivers and frameworks for these devices typically require separate languages and implementations for each vendor. Furthermore, there are few tools and little support to schedule codes between these devices in a truly heterogeneous manner-partly because of this fragmentation between vendors and the languages each supports. To overcome both limitations, the Intelligent Runtime System (IRIS) was developed. It allows a common task abstraction to automatically be shared among contemporary vendors and is run from a single host-side API. At runtime, IRIS queries the host system and registers which frameworks and drivers are available, these determine which kernels can be used by the scheduler-CPUs via OpenMP, Nvidia GPUs (CUDA), AMD GPUs (HIP), and Intel and Xilinx FPGAs with OpenCL. IRIS enables tasks to be scheduled to any heterogeneous device and resolves to the appropriate kernel binary at runtimeit only uses the devices supported by the system on which it is run. IRIS supports single-task and graph-based expressions of dependencies of tasks. Additionally, IRIS features a range of dynamic scheduling policies, allowing complex chains of tasks and interactions to be executed, relieving the programmer/user from considering the system to assign tasks to devices optimally. This paper presents the peak performance attainable by IRIS over a range of systems-each with different numbers and types of accelerator devices, it highlights the flexibility of IRIS since these devices are truly heterogeneous, relying on different backends (drivers, frameworks, and languages) which historically required unique implementations to utilize them. We then use this peak performance as a baseline to compare increasingly complex chains of tasks (with increasingly complex task dependencies) and evaluate how IRIS copes. Finally, we consider the performance of different IRIS scheduling policies on this range of task graphs.

Johnston, Beau↗
DOE OSTI2022

Modelling of edge plasma dynamics with active wall boundary conditions

A self-consistent 2D model is presented for transport in boundary plasma and plasma-facing material walls. Plasma dynamics in the domain is represented by a 2D collisional plasma fluid model in the edge-plasma code UEDGE, and transport of hydrogen and heat in the wall is represented by a system of reaction–diffusion equations in the 1D wall code FACE. To account for variation of parameters along the wall, in the coupled model multiple instances of the FACE code run in parallel. Here, the coupled model provides a tool for investigating a range of dynamic plasma–material interactions phenomena in 2D. For demonstration of its capability, one application of particular interest is the role of active wall in tokamak strike point sweeping proposed for mitigation of divertor heat loads. In the present study, the coupled calculations are applied to investigation of the impact of heat and hydrogen transport in the material wall on the divertor plasma and target heat load during sweeping of the target strike point for parameters of a high-power tokamak.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Ultra-Fast Non-Volatile Resistive Switching Devices with Over 512 Distinct and Stable Levels for Memory and Neuromorphic Computing

Low-current multilevel programmability with inherent non-volatility and high stability of resistance states is required for both multi-bit memory storage and deep learning accelerators but is difficult to achieve. Here, in a resistive switching system, this work realizes >512 (>9 bits) distinct non-volatile conductance levels with stable retention for each state with current levels down to the nanoampere range, highly promising for potential integration with small processing nodes with ultra-low power consumption requirements. This is achieved by demonstrating a new thin film design concept that encompasses three key features: an ultra-thin epitaxial oxygen ionic switching layer that provides a tunable energy barrier at the bottom electrode, an overcoat amorphous layer that acts as an ion migration barrier for stable state retention, and a partial conductive filament as a localized electronic transport channel to the epitaxial switching layer. A large dynamic resistance range of up to seven orders of magnitude is achieved with reset-free transitions among intermediate states, and programmability is demonstrated with ultra-fast (20 ns) pulses. Artificial neural network (ANN) simulations, based on the experimental performance and its non-idealities, demonstrate close-to-ideal inference accuracies for various Modified National Institute of Standards and Technology (MNIST) data sets.

36 MATERIALS SCIENCE↗
DOE OSTI2020

X-Ray Afterglows from the Gamma-Ray Burst “Large-angle” Emission

In this work, we derive basic analytical results for the timing and decay of the gamma-ray burst (GRB) counterpart and delayed afterglow light curves for a brief emission episode from a relativistic surface endowed with angular structure, consisting of a uniform core of size ${\theta }_{c}$ (Lorentz factor ${{\rm{\Gamma }}}_{c}$ and surface emissivity ${i}_{\nu ^{\prime} }^{{\prime} }$ are angle independent) and an axially symmetric power-law envelope (${\rm{\Gamma }}\sim {\theta }^{-g}$). In this "large-angle emission" model, radiation produced during the prompt emission phase (GRB) at angles $\theta \gt {\theta }_{c}$ arrives at the observer well after the burst (delayed emission). The dynamical time range of the very fast decaying GRB "tail" and of the flat afterglow "plateau" and the morphology of the GRB counterpart/afterglow are all determined by two parameters: the core's parameter ${{\rm{\Gamma }}}_{c}{\theta }_{c}$ and the envelope's Lorentz factor index g, leading to three types of light curves that display three post-GRB phases (type 1: tail, plateau/slow decay, post-plateau/normal decay), two post-GRB phases (type 2: tail and fast decay), or just one (type 3: normal decay). We show how X-ray light-curve features can be used to determine core and envelope dynamical and spectral parameters. Testing of the large-angle emission model is done using the Swift/XRT X-ray emission of two afterglows of type 1 (GRB 060607A, GRB 061121), one of type 2 (GRB 061110A), and one of type 3 (GRB 061007). We find that the X-ray afterglows with plateaus require an envelope Lorentz factor ${\rm{\Gamma }}\sim {\theta }^{-2}$ and a comoving-frame emissivity ${i}_{\nu ^{\prime} }^{{\prime} }\sim {\theta }^{2}$; thus, for a typical afterglow spectrum ${F}_{\nu }\sim {\nu }^{-1}$, the lab-frame energy release is uniform over the emitting surface.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

NESSi : The N on- E quilibrium S ystems S imulation package

The nonequilibrium dynamics of correlated many-particle systems is of interest in connection with pump–probe experiments on molecular systems and solids, as well as theoretical investigations of transport properties and relaxation processes. Nonequilibrium Green’s functions are a powerful tool to study interaction effects in quantum many-particle systems out of equilibrium, and to extract physically relevant information for the interpretation of experiments. Here, we present the open-source software package NESSi (The Non-Equilibrium Systems Simulation package) which allows to perform many-body dynamics simulations based on Green’s functions on the L-shaped Kadanoff–Baym contour. NESSi contains the library libcntr which implements tools for basic operations on these nonequilibrium Green’s functions, for constructing Feynman diagrams, and for the solution of integral and integro-differential equations involving contour Green’s functions. The library employs a discretization of the Kadanoff–Baym contour into time points and a high-order implementation of integration routines. The total integrated error scales up to $\mathcal{O}(N^{-7})$, which is important since the numerical effort increases at least cubically with the simulation time. A distributed-memory parallelization over reciprocal space allows large-scale simulations of lattice systems. We provide a collection of example programs ranging from dynamics in simple two-level systems to problems relevant in contemporary condensed matter physics, including Hubbard clusters and Hubbard or Holstein lattice models. The libcntr library is the basis of a follow-up software package for nonequilibrium dynamical mean-field theory calculations based on strong-coupling perturbative impurity solvers.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

CERF: IM3 Projected Western US Power Plant Locations

Overview The Capacity Expansion Regional Feasibility (CERF) model is an open-source geospatial python package that provides new power plant locations at a 1km resolution. The model ingests U.S. state or regional-scale electricity system capacity expansion plans, such as those produced by the Global Change Analysis Model (GCAM-USA), and identifies feasible, site-specific locations for individual new power plants (renewable and non-renewable). CERF combines high-resolution geospatial suitability analyses with an economic algorithm that selects individual plant siting locations based on grid interconnection costs and the locational marginal value of new generation. The model incorporates a wide range of dynamic constraints and opportunities, such as protected lands, population density, existing infrastructure, and water availability. This dataset provides CERF power plant siting results for IM3 Phase 2 simulations across eight different scenarios for the Western US through 2055. The scenarios include combinations of two Shared Socioeconomic Pathways (SSP3 and SSP5) with four high-resolution climate projections specific to the United States (see, https://tgw-data.msdlive.org/). These climate projections include "hotter" and "cooler" variants for two Representative Concentration Pathways (RCP4.5 and RCP8.5). The resulting eight simulations are: rcp45cooler_ssp3 rcp45cooler_ssp5 rcp45hotter_ssp3 rcp45hotter_ssp5 rcp85cooler_ssp3 rcp85cooler_ssp5 rcp85hotter_ssp3 rcp85hotter_ssp5 CERF siting results in this dataset correspond to capacity expansion plans in the GCAM-USA IM3 Phase 2 simulation data and are available for each of the above scenarios. Data Details Temporal Range: 2015-2055 in 5-year timesteps. Note that 2015 is the experiment base year and 2020 and beyond represent model simulation years. Spatial Range: Plant locations are provided for the eleven states in the Western US including Arizona, California, Colorado, Idaho, Montana, New Mexico, Nevada, Oregon, Utah, Washington, and Wyoming. Spatial Resolution: 1 km-squared, provided in x and y coordinates Geospatial Projection: Albers Equal Area Conic (ESRI:102003) File Type: csv The dataset contains subdirectories for each of the eight scenarios described in the overview. Each scenario folder contains two subfolders with the following information: 1. Power Plant Data This directory contains a single .csv file of power plant locations for both pre-existing (non-CERF sited plants in operation in 2015) and new (CERF-sited) power plants across the temporal range along with additional CERF model output parameters for CERF-sited plants. Plant with a siting year earlier than 2020 correspond to facilities that are operational leading into the first timestep CERF simulation. For a more detailed description of CERF model output parameters, see the CERF model documentation. Note that the cerf_plant_id parameter is unique within each scenario file but not across scenario files. Parameter Descriptions scenario - Name of scenario cerf_plant_id - Unique siting identifier cerf_sited - If True, indicates that plant was sited by CERF model. If False, indicates pre-existing facility region_name - Name of region (state) tech_id - Technology ID tech_name - Full generation technology name inclusive of cooling type (if applicable) and additional characteristics tech_simple - Simplified generation technology type unit_size_mw - Power plant unit size (MW) xcoord - X coordinate in the default CRS (meters) ycoord - Y coordinate in the default CRS (meters) index - Index position in the flattend 2D array buffer_in_km - Exclusion buffer around site (km) sited_year - Year of siting retirement_year - Year of retirement lmp_zone - Locational marginal price (LMP) zone ID locational_marginal_price_usd_per_mwh - Locational marginal price ($/MWh) generation_mwh_per_year - Generation output (MWh/yr) operating_cost_usd_per_year - Cost of plant operations ($/yr) net_operational_value - Net operational value based on LMP and and operating costs ($/yr) interconnection_cost - Cost of interconnection for transmission & gas pipeline (if applicable) net_locational_cost -- Difference of interconnection cost and operating value ($/yr) capacity_factor_fraction - Capacity factor (fraction) carbon_capture_rate_fraction - Carbon capture rate (fraction) fuel_co2_content_tons_per_btu - Fuel CO2 content (tons/Btu) fuel_price_usd_per_mmbtu - Fuel price ($/MMBtu) fuel_price_esc_rate_fraction - Fuel price escalation rate (fraction) heat_rate_btu_per_kWh - Heat rate (Btu/kWh) lifetime_yrs - Technology lifetime for annuity (years) operational_life_yrs - Operational lifetime for retirement (years) variable_om_usd_per_mwh - Variable operation and maintenance costs of yearly capacity use ($/MWh) variable_om_esc_rate_fraction - Variable operation and maintenance costs escalation rate (fraction) carbon_tax_usd_per_ton - Carbon tax ($/ton) carbon_tax_esc_rate_fraction - Carbon tax escalation rate (fraction) 2. Storage Data This directory contains information on new and pre-existing energy storage facilities operational in each timestep along with various storage operational parameters. The 2015 timestep provides pre-existing energy storage data and corresponds with facilities that are operational leading into the first model simulation timestep. Note that coordinates in the storage files correspond to the interconnection point on the grid (substation location), not individual energy storage locations. Energy storage is added in a cumulative process at each given interconnection point. That is, each individual file provides the total operational storage capacity interconnected to the specified substation for the given timestep, inclusive of previously installed storage at that location and new storage installed in that timestep at that location. Parameters scenario - Name of scenario timestep - Simulation timestep name - Unique storage identifier s_typ - Type of energy storage technology (battery or pumped storage hydro) s_node - Node ID of interconnecting substation xcoord - X coordinate in the default CRS (meters) ycoord - Y coordinate in the default CRS (meters) charge_rate - Maximum charge rate (power capacity) of storage system (MW) discharge_rate - Maximum discharge rate (power capacity) of storage system (MW) duration - Duration of storage system (hours) max_SoC - Allowed maximum state of charge (energy capacity) of storage system (MWh) min_SoC -Allowed minimum state of charge (energy capacity) of storage system (MWh) charge_eff - Efficiency of charge (fraction between 0 and 1) discharge_eff - Efficiency of discharge (fraction between 0 and 1) Acknowledgment IM3 is a multi-institutional effort led by Pacific Northwest National Laboratory and supported by the U.S. Department of Energy's Office of Science as part of research in MultiSector Dynamics, Earth and Environmental Systems Modeling Program.

CERF↗
DOE OSTI2024

Double diffraction imaging of x-ray induced structural dynamics in single free nanoparticles

Abstract Because of their high photon flux, x-ray free-electron lasers (FEL) allow to resolve the structure of individual nanoparticles via coherent diffractive imaging (CDI) within a single x-ray pulse. Since the inevitable rapid destruction of the sample limits the achievable resolution, a thorough understanding of the spatiotemporal evolution of matter on the nanoscale following the irradiation is crucial. We present a technique to track x-ray induced structural changes in time and space by recording two consecutive diffraction patterns of the same single, free-flying nanoparticle, acquired separately on two large-area detectors opposite to each other, thus examining both the initial and evolved particle structure. We demonstrate the method at the extreme ultraviolet (XUV) and soft x-ray Free-electron LASer in Hamburg (FLASH), investigating xenon clusters as model systems. By splitting a single XUV pulse, two diffraction patterns from the same particle can be obtained. For focus intensities of about 2 × 10 12 W cm −2 we observe still largely intact clusters even at the longest delays of up to 650 picoseconds of the second pulse, indicating that in the highly absorbing systems the damage remains confined to one side of the cluster. Instead, in case of five times higher flux, the diffraction patterns show clear signatures of disintegration, namely increased diameters and density fluctuations in the fragmenting clusters. Future improvements to the accessible range of dynamics and time resolution of the approach are discussed.

Sauppe, M. (ORCID:0000000163123102)↗
DOE OSTI2020

Three-Dimensional Imaging Lidar for Characterizing Particle Fields and Organisms in the Mesopelagic Zone

The ocean’s mesopelagic zone is largely uncharacterized despite its vital role in sustaining ocean ecosystems. The composition, cycling, and fate of particle fields in the mesopelagic lacks an integrative multi-scale understanding of organism migration patterns, distribution, and diversity. This problem is addressed by combining complementary technologies with overlapping size spectra, including profiler mounted optical scattering sensors, profiler, and ship mounted acoustic devices, and a custom Unobtrusive Multi-Static Lidar Imager (UMSLI). This unique sensor suite can observe distributions of particles including organisms over a six order of magnitude dynamic size range, from microns to meters. Overlapping size ranges between different methods allows for cross-validation. This work focuses on the lidar imaging measurements and optical backscattering and attenuation, covering a combined particle size range of 0.1 mm to several cm. Particles at the small end of this range are sized using an existing backscattering time series inversion method after Briggs et al. (2013). Larger particles are resolved with UMSLI over an expanding volume using three-dimensional photo-realistic laser serial imaging. UMSLI’s image rectifying ability over time allows for derivation of particle concentration, size, and spatial distribution. Technical details on the development and post-processing methods for the novel UMSLI system are provided. Image resolved particle size distributions (PSDs) revealed a size shift from smaller to larger particles (>0.5 mm) as indicated by flatter slopes from dawn (slope = 2.6) to dusk (slope = 3.0). PSD trends are supported by an optical backscatter and transmissometer time series inversion analysis. Size shifts in the particle field are largely attributed to aggregation effects. Images support evidence of temporal variation between dusk and dawn stations through statistical analysis of particle concentrations for particle sizes 0.50–5.41 mm. Spatial analysis of the particle field revealed a dominantly uniform distributed marine snow background. The importance and potential of integrated approaches to studying particle and organism dynamics in ocean environments are discussed.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

A constitutive model for elastomers tailored by ionic bonds and entanglements

Over the past decade or two, the concept has emerged of using multiple types of weak interactions simultaneously to enhance the mechanical properties of elastomers. These weak interactions include physical entanglements, hydrogen bonds, metal-coordination bonds, dynamic covalent bonds, and ionic bonds. The combination of entanglements and ionic bonding has been minimally explored and is particularly exciting because of the broad application space for polyelectrolytes. In this work, a constitutive model framework is developed to describe the response of elastomers with both ionic bonds and entanglements. Here we formulate a micromechanical model that couples together chain stretching, ionic bond slipping, and entanglement evolution. The ionic bonds provide toughness by enabling plastic deformation in comparison to covalently crosslinked material and add strength compared to a linear polymer. Evolution of the entanglement density is taken as a key mechanism that can govern stiffness, toughness, and self-recovery in elastomers. The model is used to match bulk polyelectrolytes with different fractions of ionic components under a variety of loading histories. The variations in material parameters are then used to help understand the relative importance of different governing mechanisms in the bulk polymers. We show that the theoretical framework can explain our experimental uniaxial tensile experimental results for polyelectrolytes. This model can help to design better material with high stiffness and toughness. We expect that our model can be extended to explain the mechanical behavior of other polyelectrolytes and other soft materials with a wide range of dynamic bonds.

42 ENGINEERING↗
DOE OSTI2020

NaGaS 2 : An Elusive Layered Compound with Dynamic Water Absorption and Wide‐Ranging Ion‐Exchange Properties

Abstract Most ternary sulfides belonging to the MGaS 2 structure‐type have been known for many years and are well‐characterized. Surprisingly, there have been no reports of the NaGaS 2 composition, which contains Na, a monovalent cation slightly larger in size than Li, found in LiGaS 2 , a compound known for its non‐linear optical properties. Now it is demonstrated for the first time that the unique reversible water absorption in NaGaS 2 has resulted in its absence from previous reports owing to difficulties encountered when characterizing this compound by SC XRD. The layered structure of this compound coupled with uniquely easy migration of water molecules between the layers allows for ion exchange with 3d and 5f metal cations. Some cations, for example, Ni 2+ , facilitate exfoliation of the layers, providing a facile synthetic route to a new class of 2D chalcogenide materials and furthermore demonstrating that NaGaS 2 can readily uptake uranyl species from aqueous solutions.

Klepov, Vladislav V.↗
DOE OSTI2021

Protein Assembly by Design

Proteins are nature’s primary building blocks for the construction of sophisticated molecular machines and dynamic materials, ranging from protein complexes such as photosystem II and nitrogenase that drive biogeochemical cycles to cytoskeletal assemblies and muscle fibers for motion. Such natural systems have inspired extensive efforts in the rational design of artificial protein assemblies in the last two decades. As molecular building blocks, proteins are highly complex, in terms of both their three-dimensional structures and chemical compositions. To enable control over the self-assembly of such complex molecules, scientists have devised many creative strategies by combining tools and principles of experimental and computational biophysics, supramolecular chemistry, inorganic chemistry, materials science, and polymer chemistry, among others. Owing to these innovative strategies, what started as a purely structure-building exercise two decades ago has, in short order, led to artificial protein assemblies with unprecedented structures and functions and protein-based materials with unusual properties. Furthermore, our goal in this review is to give an overview of this exciting and highly interdisciplinary area of research, first outlining the design strategies and tools that have been devised for controlling protein self-assembly, then describing the diverse structures of artificial protein assemblies, and finally highlighting the emergent properties and functions of these assemblies.

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

Perspective on Many-Body Methods for Molecular Polaritonic Systems

Recent advances in strong light–matter interactions have revealed a wealth of new physical phenomena in molecules embedded in optical cavities, including modified chemical reactivity, altered excitation spectra, and novel quantum correlations. To describe these effects from first-principles, the field of ab initio quantum electrodynamics (QED) has emerged as a compelling extension of quantum chemistry that treats electronic and photonic degrees of freedom on equal footing. In this Perspective, we review the growing landscape of many-body QED methods, including Hartree–Fock, density functional theory (QEDFT), time-dependent DFT (QED-TDDFT), configuration interaction (QED-CI), complete active space (QED-CASSCF), coupled cluster (QED-CC), quantum Monte Carlo (QED-QMC), and density matrix renormalization group (QED-DMRG), highlighting recent developments and implementations. We further explore real-time methods, gradient and Hessian formalisms, and the integration of nonadiabatic nuclear dynamics. Applications range from benchmark simulations of polaritonic chemistry to quantum simulations on emerging quantum hardware. We conclude by outlining future directions for theory development and interdisciplinary efforts at the interface of quantum chemistry, condensed matter, and quantum optics.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Quasi-one-dimensional hydrogen bonding in nanoconfined ice

Abstract The Bernal-Fowler ice rules stipulate that each water molecule in an ice crystal should form four hydrogen bonds. However, in extreme or constrained conditions, the arrangement of water molecules deviates from conventional ice rules, resulting in properties significantly different from bulk water. In this study, we employ machine learning-driven first-principles simulations to identify a new stabilization mechanism in nanoconfined ice phases. Instead of forming four hydrogen bonds, nanoconfined crystalline ice can form a quasi-one-dimensional hydrogen-bonded structure that exhibits only two hydrogen bonds per water molecule. These structures consist of strongly hydrogen-bonded linear chains of water molecules that zig-zag along one dimension, stabilized by van der Waals interactions that stack these chains along the other dimension. The unusual interplay of hydrogen bonding and van der Waals interactions in nanoconfined ice results in atypical proton behavior such as potential ferroelectric behavior, low dielectric response, and long-range proton dynamics.

Science & Technology - Other Topics↗
DOE OSTI2026

Mapping the Perseus galaxy cluster with XRISM: Gas kinematic features and their implications for turbulence

We present extended gas kinematic maps of the Perseus cluster based on a combination of five new XRISM/Resolve pointings observed in 2025 with four performance verification datasets from 2024, totaling a net exposure of 745 ks. To date, Perseus remains the only cluster that has been extensively mapped out to ≃0.7 r 2500 by XRISM/Resolve, while simultaneously offering sufficient spatial resolution to resolve gaseous substructures driven by mergers and active galactic nucleus (AGN) feedback. Our observations cover multiple radial directions and a broad range of dynamical scales, enabling us to characterize the kinematic properties of the intracluster medium up to a scale of ∼500 kpc. In the measurements, we detected high-velocity dispersions (≃300km s −1 ) in the eastern region of the cluster that are spatially coincident with the extended X-ray surface brightness excess and correspond to a nonthermal pressure fraction of ≃7 − 13%. The velocity field outside the AGN-dominant region can be effectively described by a single, large-scale kinematic driver based on the velocity structure function, which statistically favors an energy injection scale of at least a few hundred kpc. The estimated turbulent dissipation energy is comparable to the gravitational potential energy released by a recent merger, implying a significant role of turbulent cascade in the merger energy conversion. In the bulk velocity field, we observed a dipole-like pattern along the east-west direction with an amplitude of ≃ ± 200 − 300 km s −1 , indicating rotational motions induced by the recent merger event. This feature constrains the viewing direction to ≃30° −50° relative to the normal of the merger plane. Our hydrodynamic simulations suggest that Perseus has experienced at least two energetic mergers since redshift z ∼ 1, the most recent of which is associated with the radio galaxy IC310, in agreement with recent SRG/eROSITA findings. This study showcases exciting scientific opportunities for future missions with high-resolution spectroscopic capabilities (e.g., HUBS, LEM, and NewAthena).

79 ASTRONOMY AND ASTROPHYSICS↗