Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “RESPONSE FUNCTIONS”

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

Graph-Based Modeling for the Detection and Tracking of Sarin-Surrogate-Induced Neurotoxicity Using a Human-Relevant, In-Vitro Brain Model

Organophosphorus (OP) nerve agents are a chemical threat to the United States, to the civilian population (e.g., pesticides) and historically weaponized (e.g., sarin) as chemical warfare agents. The unprecedented, accelerated process from “bench-to-bedside” during the SARSCov2 pandemic has made it clear that technology and tools need to be readily available for immediate response. Advances in human organ tissue mimetic systems are a promising technology to evaluate the human-relevant response in vitro for basic and applied research and drug screening. In particular, current brain microphysiological systems (MPS) have the capability to monitor and detect changes in engineered human neural circuit activity. However, current data analytics approaches for these systems lack the granularity to functionally detect and distinguish the different mechanisms that occur in the brain following neurotoxicity, injury, and disease. The goal of this project was to advance the computational analytical capabilities of the brain MPS to detect functional changes in neural circuit structure at different stages of Sarin surrogate-induced neurotoxicity. We developed graph-based models to (1) identify the composition of the neural circuit structure; (2) detect and monitor how this structure changes following sarin-induced neurotoxicity; and (3) evaluate the analytical pipeline using known/promising oxime reactivators. Through experiments on the bMPS where in vitro neuronal cultures were exposed to a sarin surrogate, we demonstrated the capabilities of our computational pipeline to identify different responses in the functional networks of brain cells exposed to low and high concentrations of the nerve agent. We identified a biphasic response of human neural network activity following exposure to a sarin-surrogate that had not been reported in the literature before. The graph-based models and software developed in this project can be used for future studies that leverage the brain MPS technology, such as treatment efficacy assessment.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Giant Doping Response of Magnetic Anisotropy in MnTe

Developing simple ways to control spin states in spintronic devices is a crucial step towards increasing their functionality. MnTe is a room-temperature antiferromagnet with promising spintronic properties, including for thermospintronics and magnon-based devices. Here, we show that, in MnTe, less than 1% Li is sufficient to produce a dramatic spin reorientation as observed by neutron diffraction. The behavior of the 0001 magnetic Bragg peak reveals a significant reorientation of the Mn spins from planar in the pure material to almost completely axial with minimal Li doping. Temperature dependence of the magnetic peaks in Li-doped samples indicates that axial spins shift back to planar suddenly upon approaching the Néel temperature (T N = 307 K). Density functional theory calculations support the idea that a shift in the Fermi level caused by doping is responsible for switching the material between two competing magnetic ground states. These results pave the way for developing easy switching of magnetic states in functional materials such as spintronic devices and topological insulators.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Perennial grass root system specializes for multiple resource acquisitions with differential elongation and branching patterns

Roots optimize the acquisition of limited soil resources, but relationships between root forms and functions have often been assumed rather than demonstrated. Furthermore, how root systems co-specialize for multiple resource acquisitions is unclear. Theory suggests that trade-offs exist for the acquisition of different resource types, such as water and certain nutrients. Measurements used to describe the acquisition of different resources should then account for differential root responses within a single system. To demonstrate this, we grew Panicum virgatum in split-root systems that vertically partitioned high water availability from nutrient availability so that root systems must absorb the resources separately to fully meet plant demands. We evaluated root elongation, surface area, and branching, and we characterized traits using an order-based classification scheme. Plants allocated approximately 3/4th of primary root length towards water acquisition, whereas lateral branches were progressively allocated towards nutrients. However, root elongation rates, specific root length, and mass fraction were similar. Our results support the existence of differential root functioning within perennial grasses. Similar responses have been recorded in many plant functional types suggesting a fundamental relationship. Root responses to resource availability can be incorporated into root growth models via maximum root length and branching interval parameters.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

CISME report on experiments with HMX powder, PBX 9501 prills, and heated PBX 9501

The ability to predict post-ignition explosive response is relevant to explosives safety, particularly weapons assembly/disassembly operations that must operate under a sufficiently conservative safety basis to meet risk mitigation criteria for inadvertent nuclear detonation (IND). In particular, hazard determinations hinge on whether the PBX 9501 main charge can ultimately attain detonation from an initial, non-shock, ignition stimulus. The CISME experiment uses a spherical explosive charge, which is thermally ignited in the center using a fiber-delivered infrared laser pulse. This ignition mechanism - i.e. a localized thermal hotspot - is representative of the ignition-causing insults typical in accident scenarios, including low-velocity impact, drop, skid, and punch events. In fiscal year 2018 CISME experiments were conducted with pristine PBX 9501 spheres, confined and unconfined. Sub-detonative, low-violence reactions were observed. Here we report experiments conducted since 2018. The explosive charge was varied: experiments were conducted with room temperature pour-density HMX powder and prills of PBX 9501, and heated, thermally damaged PBX 9501. Experiments were conducted in a confined configuration, with two 4 mm-thick stainless-steel hemispherical shells bolted around the sphere. Detonation, via DDT (deflagration-to-detonation transition) was observed in the 3-inch diameter powdered HMX and the 6-inch diameter heated and thermally-damaged PBX 9501 tests. No transition to DDT was observed in smaller diameters for those configurations, and no DDT was observed for pour-density prills of PBX 9501 of any diameter. The original test series demonstrated that the response of pristine, undamaged PBX 9501 to a central ignition is relatively benign. This new test series shows that increased porosity, lack of polymer binders and thermal damage can greatly increase the violence of the post-ignition response, up to and including detonation. These results bridge the explosive response of PBX 9501 from sub-detonative through detonation as a function of diameter and provide valuable test cases for modeling efforts that strive to capture the DDT threshold. By exploring explosive response in a more representative 3D geometry, we can also benchmark these new results to more plentiful DDT datasets generated from testing in 1D tube confinement (the historical test configuration)–the configuration upon which almost all DDT models were developed.

42 ENGINEERING↗
DOE OSTI2020

Enzalutamide response in a panel of prostate cancer cell lines reveals a role for glucocorticoid receptor in enzalutamide resistant disease

Representative in vitro model systems that accurately model response to therapy and allow the identification of new targets are important for improving our treatment of prostate cancer. Here we describe molecular characterization and drug testing in a panel of 20 prostate cancer cell lines. The cell lines cluster into distinct subsets based on RNA expression, which is largely driven by functional Androgen Receptor (AR) expression. KLK3 , the AR-responsive gene that encodes prostate specific antigen, shows the greatest variability in expression across the cell line panel. Other common prostate cancer associated genes such as TMPRSS2 and ERG show similar expression patterns. Copy number analysis demonstrates that many of the most commonly gained (including regions containing TERC and MYC ) and lost regions (including regions containing TP53 and PTEN) that were identified in patient samples by the TCGA are mirrored in the prostate cancer cell lines. Assessment of response to the anti-androgen enzalutamide shows a distinct separation of responders and non-responders, predominantly related to status of wild-type AR. Surprisingly, several AR-null lines responded to enzalutamide. These AR-null, enzalutamide-responsive cells were characterized by high levels of expression of glucocorticoid receptor (GR) encoded by NR3C1 . Treatment of these cells with the anti-GR agent mifepristone showed that they were more sensitive to this drug than enzalutamide, as were several of the enzalutamide non-responsive lines. This is consistent with several recent reports that suggest that GR expression is an alternative signaling mechanism that can bypass AR blockade. This study reinforces the utility of large cell line panels for the study of cancer and identifies several cell lines that represent ideal models to study AR-null cells that have upregulated GR to sustain growth.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Oxygen Vacancy Injection as a Pathway to Enhancing Electromechanical Response in Ferroelectrics

Since their discovery in late 1940s, perovskite ferroelectric materials have become one of the central objects of condensed matter physics and materials science due to the broad spectrum of functional behaviors they exhibit, including electro-optical phenomena and strong electromechanical coupling. In such disordered materials, the static properties of defects such as oxygen vacancies are well explored but the dynamic effects are less understood. In this work, the first observation of enhanced electromechanical response in BaTiO3 thin films is reported driven via dynamic local oxygen vacancy control in piezoresponse force microscopy (PFM). A persistence in peizoelectricity past the bulk Curie temperature and an enhanced electromechanical response due to a created internal electric field that further enhances the intrinsic electrostriction are explicitly demonstrated. The findings are supported by a series of temperature dependent band excitation PFM in ultrahigh vacuum and a combination of modeling techniques including finite element modeling, reactive force field, and density functional theory. Furthermore, this study shows the pivotal role that dynamics of vacancies in complex oxides can play in determining functional properties and thus provides a new route toward– achieving enhanced ferroic response with higher functional temperature windows in ferroelectrics and other ferroic materials.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Biphasic response of human iPSC-derived neural network activity following exposure to a sarin-surrogate nerve agent

Organophosphorus nerve agents (OPNA) are hazardous environmental exposures to the civilian population and have been historically weaponized as chemical warfare agents (CWA). OPNA exposure can lead to several neurological, sensory, and motor symptoms that can manifest into chronic neurological illnesses later in life. There is still a large need for technological advancement to better understand changes in brain function following OPNA exposure. The human-relevant in vitro multi-electrode array (MEA) system, which combines the MEA technology with human stem cell technology, has the potential to monitor the acute, sub-chronic, and chronic consequences of OPNA exposure on brain activity. However, the application of this system to assess OPNA hazards and risks to human brain function remains to be investigated. In a concentration-response study, we have employed a human-relevant MEA system to monitor and detect changes in the electrical activity of engineered neural networks to increasing concentrations of the sarin surrogate 4-nitrophenyl isopropyl methylphosphonate (NIMP). We report a biphasic response in the spiking (but not bursting) activity of neurons exposed to low (i.e., 0.4 and 4 μM) versus high concentrations (i.e., 40 and 100 μM) of NIMP, which was monitored during the exposure period and up to 6 days post-exposure. Regardless of the NIMP concentration, at a network level, communication or coordination of neuronal activity decreased as early as 60 min and persisted at 24 h of NIMP exposure. Once NIMP was removed, coordinated activity was no different than control (0 μM of NIMP). Interestingly, only in the high concentration of NIMP did coordination of activity at a network level begin to decrease again at 2 days post-exposure and persisted on day 6 post-exposure. Notably, cell viability was not affected during or after NIMP exposure. Also, while the catalytic activity of AChE decreased during NIMP exposure, its activity recovered once NIMP was removed. Gene expression analysis suggests that human iPSC-derived neurons and primary human astrocytes resulted in altered genes related to the cell’s interaction with the extracellular environment, its intracellular calcium signaling pathways, and inflammation, which could have contributed to how neurons communicated at a network level.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Stoichiometric Thiol Redox Proteomics for Quantifying Cellular Responses to Perturbations

Post-translational modifications regulate the structure and function of proteins that can result in changes to the activity of different pathways. These include modifications altering the redox state of thiol groups on protein cysteine residues, which are sensitive to oxidative environments. While mass spectrometry has advanced the identification of protein thiol modifications and expanded our knowledge of redox-sensitive pathways, the quantitative aspect of this technique is critical for the field of redox proteomics. In this review, we describe how mass spectrometry-based redox proteomics has enabled researchers to accurately quantify the stoichiometry of reversible oxidative modifications on specific cysteine residues of proteins. We will describe advancements in the methodology that allow for the absolute quantitation of thiol modifications, as well as recent reports that have implemented this approach. We will also highlight the significance and application of such measurements and why they are informative for the field of redox biology.

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

EcoSpec: Highly Equipped Tower-Based Hyperspectral and Thermal Infrared Automatic Remote Sensing System for Investigating Plant Responses to Environmental Changes

Despite an advanced ability to forecast ecosystem functions and climate at regional and global scales, little is known about relationships between local variations in water and carbon fluxes and large-scale phenomena. To enable data collection of local-scale ecosystem functions to support such investigations, we developed the EcoSpec system, a highly equipped remote sensing system that houses a hyperspectral radiometer (350–2500 nm) and five optical and infrared sensors in a compact tower. Its custom software controls the sequence and timing of movement of the sensors and system components and collects measurements at 12 locations around the tower. The data collected using the system was processed to remove sun-angle effects, and spectral vegetation indices computed from the data (i.e., the Normalized Difference Vegetation Index (NDVI), Normalized Difference Water Index (NDWI), Photochemical Reflectance Index (PRI), and Moisture Stress Index (MSI)) were compared with the fraction of photochemically active radiation (fPAR) and canopy temperature. The results showed that the NDVI, NDWI, and PRI were strongly correlated with fPAR; the MSI was correlated with canopy temperature at the diurnal scale. These correlations suggest that this type of near-surface remote sensing system would complement existing observatories to validate satellite remote sensing observations and link local and large-scale phenomena to improve our ability to forecast ecosystem functions and climate. The system is also relevant for precision agriculture to study crop growth, detect disease and pests, and compare traits of cultivars.

47 OTHER INSTRUMENTATION↗
DOE OSTI2021

Kinetic Properties of Microbial Exoenzymes Vary With Soil Depth but Have Similar Temperature Sensitivities Through the Soil Profile

Current knowledge of the mechanisms driving soil organic matter (SOM) turnover and responses to warming is mainly limited to surface soils, although over 50% of global soil carbon is contained in subsoils. Deep soils have different physicochemical properties, nutrient inputs, and microbiomes, which may harbor distinct functional traits and lead to different SOM dynamics and temperature responses. We hypothesized that kinetic and thermal properties of soil exoenzymes, which mediate SOM depolymerization, vary with soil depth, reflecting microbial adaptation to distinct substrate and temperature regimes. We determined the Michaelis-Menten (MM) kinetics of three ubiquitous enzymes involved in carbon (C), nitrogen (N) and phosphorus (P) acquisition at six soil depths down to 90 cm at a temperate forest, and their temperature sensitivity based on Arrhenius/ Q 10 and Macromolecular Rate Theory (MMRT) models over six temperatures between 4–50°C. Maximal enzyme velocity ( V max ) decreased strongly with depth for all enzymes, both on a dry soil mass and a microbial biomass C basis, whereas their affinities increased, indicating adaptation to lower substrate availability. Surprisingly, microbial biomass-specific catalytic efficiencies also decreased with depth, except for the P-acquiring enzyme, indicating distinct nutrient demands at depth relative to microbial abundance. These results suggested that deep soil microbiomes encode enzymes with intrinsically lower turnover and/or produce less enzymes per cell, reflecting distinct life strategies. The relative kinetics between different enzymes also varied with depth, suggesting an increase in relative P demand with depth, or that phosphatases may be involved in C acquisition. V max and catalytic efficiency increased consistently with temperature for all enzymes, leading to overall higher SOM-decomposition potential, but enzyme temperature sensitivity was similar at all depths and between enzymes, based on both Arrhenius/ Q 10 and MMRT models. In a few cases, however, temperature affected differently the kinetic properties of distinct enzymes at discrete depths, suggesting that it may alter the relative depolymerization of different compounds. We show that soil exoenzyme kinetics may reflect intrinsic traits of microbiomes adapted to distinct soil depths, although their temperature sensitivity is remarkably uniform. These results improve our understanding of critical mechanisms underlying SOM dynamics and responses to changing temperatures through the soil profile.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Temperature dependence of viscoelastic Poisson’s ratio of cement mortar

Abstract Concrete creep research has focused primarily on uniaxial response. However, biaxially prestressed concrete structures are common, resulting in a multiaxial stress state that can complicate the behavior of a viscoelastic material like concrete. Significant creep strains may be induced in directions transverse to each principle stress due to Poisson’s effect. Past research is unclear regarding the viscoelastic or viscoplastic properties of concrete outside of uniaxial response. It has been reported in separate studies that concrete viscoelastic/viscoplastic Poisson’s ratio (VPR) is an increasing, decreasing and constant function with time, with all reported measurements performed at room temperature. In this paper, the 3D basic creep response of mature cement mortar is examined using a confined compression experiment that allows direct determination of the full stress and infinitesimal strain tensors in a single test, which enables the determination of VPR under a multiaxial stress state. For this purpose, a unique, miniature version of the standardized concrete creep frame is designed that is amenable to placing in climate chambers and temperature ovens. The experimental results indicate that the VPR of sealed, mature cement mortar is nearly constant and equal to the elastic at room temperature, while the VPR gradually increases with time when measured at 60 °C.

Construction & Building Technology↗
DOE OSTI2020

High broadband photoconductivity of few-layered MoS 2 field-effect transistors measured using multi-terminal methods: effects of contact resistance

Among the layered two dimensional semiconductors, molybdenum disulfide (MoS 2 ) is considered to be an excellent candidate for applications in optoelectronics and integrated circuits due to its layer-dependent tunable bandgap in the visible region, high ON/OFF current ratio in field-effect transistors (FET) and strong light–matter interaction properties. In this study, using multi-terminal measurements, we report high broadband photocurrent response ( R ) and external quantum efficiency (EQE) of few-atomic layered MoS 2 phototransistors fabricated on a SiO 2 dielectric substrate and encapsulated with a thin transparent polymer film of Cytop. The photocurrent response was measured using a white light source as well as a monochromatic light of wavelength λ = 400 nm–900 nm. We measured responsivity using a 2-terminal configuration as high as R = 1 × 10 3 A W -1 under white light illumination with an optical power P opt = 0.02 nW. The R value increased to 3.5 × 10 3 A W -1 when measured using a 4-terminal configuration. Using monochromatic light on the same device, the measured values of R were 10 3 and 6 × 10 3 A W -1 under illumination of λ = 400 nm when measured using 2- and 4-terminal methods, respectively. The highest EQE values obtained using λ = 400 nm were 10 5 % and 10 6 % measured using 2- and 4-terminal configurations, respectively. The wavelength dependent responsivity decreased from 400 nm to the near-IR region at 900 nm. The observed photoresponse, photocurrent–dark current ratio (PDCR), detectivity as a function of applied gate voltage, optical power, contact resistances and wavelength were measured and are discussed in detail. The observed responsivity is also thoroughly studied as a function of contact resistance of the device.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Response of salt intrusion in a tidal estuary to regional climatic forcing

Abstract Salinity distribution in a large tidal estuary is subject to estuarine adjustment under the influences of multiple physical drivers such as freshwater pulses and sea level rise, and is crucial to upstream water quality, aquaculture, and ecosystem functions of the estuary. To better understand the estuarine salinity response to climate change, the unstructured-grid Finite Volume Community Ocean Model was implemented to simulate the salt intrusion in the Delaware Bay Estuary. The model was first validated by multiple observational data sets and subsequently applied in an idealized setting to examine the response of salt front to freshwater pulses in high flow conditions, followed by a long-term drought condition supported by a multi-decadal streamflow drought analysis in the estuary. The model results showed that after the freshwater pulses the salt front location moved further upstream with sea level rise. Under the simulated long-term drought condition, the adjustment timescale of salt intrusion varies nonlinearly with sea level rise. With a significant increase in sea level rise, the adjustment timescale starts to decrease. This shift suggests a transition into a different regime where the estuary becomes more stratified, as indicated by an increasing bulk Simpson number with rising sea levels.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Maximizing the Radar Generalized Image Quality Equation for Bistatic SAR Using Waveform Frequency Agility

The radar generalized image quality equation (RGIQE) is a metric used to measure both monostatic and bistatic synthetic aperture radar (BSAR) image quality, it is a function of signal-to-noise ratio (SNR) and 2-D bandwidth. The 2-D bandwidth is equal to the area of the transfer function’s (TF) passband region. With the exception of side-looking monostatic geometries, almost all monostatic and bistatic geometries have skewed passband shapes when waveform frequency parameters remain unchanged from pulse to pulse. Most synthetic aperture radar (SAR) applications require a rectangular-shaped passband region, this is achieved by inscribing a rectangular region within the skewed intrinsic passband region. Increasing skewness results in less inscription area reducing 2-D bandwidth, image SNR, and thus RGIQE capacity. In this article, a waveform with frequency agility is used to rectify the skewness that degrades RGIQE capacity. By changing the waveform’s center frequency and instantaneous bandwidth from pulse to pulse in a particular manner, the intrinsic passband region can be de-skewed. The de-skewed shape maximizes the inscription area thus maximizing 2-D bandwidth, image SNR, and RGIQE capacity. Here, three examples are given in this article, one monostatic geometry, and two bistatic geometries. RGIQE capacity is increased by 52.02%, 44.42%, and 79.09% for the three examples.

47 OTHER INSTRUMENTATION↗
DOE OSTI2022

Investigation of plasmon relaxation mechanisms using nonadiabatic molecular dynamics

Hot carriers generated from the decay of plasmon excitation can be harvested to drive a wide range of physical or chemical processes. However, their generation efficiency is limited by the concomitant phonon-induced relaxation processes by which the energy in excited carriers is transformed into heat. However, simulations of dynamics of nanoscale clusters are challenging due to the computational complexity involved. Here, in this paper, we adopt our newly developed Trajectory Surface Hopping (TSH) nonadiabatic molecular dynamics algorithm to simulate plasmon relaxation in Au 20 clusters, taking the atomistic details into account. The electronic properties are treated within the Linear Response Time-Dependent Tight-binding Density Functional Theory (LR-TDDFTB) framework. The relaxation of plasmon due to coupling to phonon modes in Au 20 beyond the Born–Oppenheimer approximation is described by the TSH algorithm. The numerically efficient LR-TDDFTB method allows us to address a dense manifold of excited states to ensure the inclusion of plasmon excitation. Starting from the photoexcited plasmon states in Au 20 cluster, we find that the time constant for relaxation from plasmon excited states to the lowest excited states is about 2.7 ps, mainly resulting from a stepwise decay process caused by low-frequency phonons of the Au 20 cluster. Furthermore, our simulations show that the lifetime of the phonon-induced plasmon dephasing process is ~10.4 fs and that such a swift process can be attributed to the strong nonadiabatic effect in small clusters. Our simulations demonstrate a detailed description of the dynamic processes in nanoclusters, including plasmon excitation, hot carrier generation from plasmon excitation dephasing, and the subsequent phonon-induced relaxation process.

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

Changes in polarization dictate necessary approximations for modeling electronic deexcitation intensity: Application to x-ray emission

Accurate simulation of electronic excitations and deexcitations are critical for complementing complex spectroscopic experiments and can provide validation to theoretical approaches. Here, using a generalized framework, we contrast the accuracy and validity of orbital-constrained and linear-response approaches that build upon Kohn-Sham density functional theory (DFT) to simulate emission spectra of electronic origin and propose an efficient approximation, named many-body x-ray emission spectroscopy or MBXES, for simulating such processes. We show analytically as well as with computed examples that for electronic (de)excitation leading to an appreciable change in polarization (i.e., density rearrangement), the adiabatic approximation in a response-based formalism will be inadequate for the calculation of oscillator strength. Thus, such a change (e.g., in the net electrostatic dipole moment of a finite system) can be used as a metric for evaluating the applicability of the adiabatic response-based approach and can be particularly valuable in x-ray emission spectroscopy. On the other hand, MBXES, the flexible method introduced in this paper, can compute oscillator strengths accurately at a much lower computational expense on the basis of two DFT-based self-consistent field calculations. Using illustrative examples of emission spectra, the efficacy of the MBXES method is demonstrated by comparison with its parent theory, orbital-optimized DFT, and with experiments.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Too big, too small, or just right? A benchmark assessment of density functional theory for predicting the spatial extent of the electron density of small chemical systems

Multipole moments are the first-order responses of the energy to spatial derivatives of the electric field strength. The quality of density functional theory prediction of molecular multipole moments thus characterizes errors in modeling the electron density itself, as well as the performance in describing molecules interacting with external electric fields. However, only the lowest non-zero moment is translationally invariant, making the higher-order moments origin-dependent. Therefore, instead of using the 3 × 3 quadrupole moment matrix, we utilize the translationally invariant 3 × 3 matrix of second cumulants (or spatial variances) of the electron density as the quantity of interest (denoted by K). The principal components of K are the square of the spatial extent of the electron density along each axis. A benchmark dataset of the principal components of K for 100 small molecules at the coupled cluster singles and doubles with perturbative triples at the complete basis set limit is developed, resulting in 213 independent K components. The performance of 47 popular and recent density functionals is assessed against this Var213 dataset. Several functionals, especially double hybrids, and also SCAN and SCAN0 predict reliable second cumulants, although some modern, empirically parameterized functionals yield more disappointing performance. The H, Li, and Be atoms, in particular, are challenging for nearly all methods, indicating that future functional development could benefit from the inclusion of their density information in training or testing protocols.

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

Intracellular Biomacromolecule Delivery by Stimuli-Responsive Protein Vesicles Loaded by Hydrophobic Ion Pairing

Proteins can perform ideal therapeutic functions. However, their large size and significant surface hydrophilicity and charge prohibit them from reaching intracellular targets. These chemical features also render them poorly encapsulated by nanoparticles used for intracellular delivery. In this work, a novel combination of protein vesicles and hydrophobic ion pairing (HIP) was used to load protein cargo and achieve cytosolic delivery to overcome the limitations of previous protein vesicle properties. Protein vesicles are thermally self-assembling nanoparticles made from elastin-like polypeptide (ELP) fused to an arginine-rich leucine zipper and a globular protein fused to a glutamate-rich leucine zipper. To impart stimuli-responsive disassembly, physiological stability, and small size, the ELP sequence was modified to include histidine and tyrosine residues. HIP was used to load and release protein cargo requiring endosomal escape for cytosolic function. HIP vesicles enabled delivery of cytochrome c, a cytosolically active protein, and a significant reduction in viability in both a traditional two-dimensional (2D) human cancer cell line culture and a biomimetic three-dimensional (3D) organoid model of acute myeloid leukemia. By examining the uptake of positively and negatively charged fluorescent protein cargos loaded by HIP, this work revealed the necessity of HIP for cytosolic cargo delivery and how HIP loading influences protein vesicle self-assembly and disassembly using microscopy, small-angle X-ray scattering, and nanoparticle tracking analysis. HIP protein vesicles have the potential to broaden the use of intracellular proteins as therapeutics for various diseases and extend protein vesicles to deliver other biomacromolecules, as the strategy developed here resulted in the first cytosolic protein cargo delivery using protein vesicles.

59 BASIC BIOLOGICAL SCIENCES↗