SEARCH · Engineering Papers

Results for “Functional priors”

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 685 records · Page 38

All sources NASA NTRS DOE OSTI

DOE OSTI2026

Excipient screening by lyophilization provides insights into spray drying formulations for nanoparticle vaccines

Nanoparticles have shown great promise as delivery platforms in the development of tunable and safe vaccines. Nanolipoprotein particles (NLPs), also known as nanodiscs, are discoidal nanoparticles composed of a lipid bilayer stabilized at their periphery by apolipoproteins. Under the right conditions, the NLP self-assembly process is highly customizable in terms of lipids and apolipoprotein constituents, allowing for tunable physical and chemical characteristics. This flexibility allows a wide range of vaccine antigens and adjuvants to be incorporated onto the NLP platform for tailored vaccine design. The stability of NLPs during long term storage is a very important factor in developing a vaccine delivery platform suitable for widespread global use. When stored in a solution for extended periods of time, NLPs dissociate into their corresponding lipids and protein constituents, leading to particle degradation. Proper stabilization of NLPs can often be achieved by lyophilization (i.e. freeze-drying), a method widely used for various applications including pharmaceuticals. This process, however, can be damaging to particles without the presence of lyoprotectants or excipients that help maintain particle stability during lyophilization. Another method used to stabilize vaccines and pharmaceuticals is spray drying, a process that converts liquid formulations into dry powders through controlled heating and airflow. While spray drying is rapid, scalable, and cost-effective, lyophilization is typically a gentler process that better retains biomolecule structure and function. Both processes use excipients for particle stabilization, so lyophilization can be used as a surrogate to test stability of NLPs, to down-select formulations that may withstand the harsher conditions of spray drying. To screen different formulations, NLPs were synthesized and purified to homogeneity by size exclusion chromatography (SEC) and samples were prepared with a wide range of lyoprotectants and/or excipients. To assess the protective effects of excipients on NLPs upon spray drying, both pre- and post-lyophilized samples were analyzed by SEC. To assess the protective effects upon heating (encountered during the spray drying process), NLP samples were incubated at elevated temperatures prior to SEC analysis. The lyoprotectants and excipients evaluated in this study had different efficiencies in protecting NLPs during lyophilization and heating tests. Trehalose, for example, exhibits stabilization on NLPs both upon lyophilization and heating whereas leucine accelerated NLP dissociation. Although some lyoprotectants are effective by themselves, different combinations can decrease the stabilization of NLPs. Shelf-stable vaccines that do not require cold-chain storage are essential for global accessibility and our findings provide fundamental insight into how to advance NLP-based vaccines for these applications.

Serrano, Litzay J [Lawrence Livermore National Lab↗

DOE OSTI2020

Discovery of Self-Assembling π-Conjugated Peptides by Active Learning-Directed Coarse-Grained Molecular Simulation

Electronically active organic molecules have demonstrated great promise as novel soft materials for energy harvesting and transport. Self-assembled nanoaggregates formed from π-conjugated oligopeptides composed of an aromatic core flanked by oligopeptide wings offer emergent optoelectronic properties within a water-soluble and biocompatible substrate. Nanoaggregate properties can be controlled by tuning core chemistry and peptide composition, but the sequence–structure–function relations remain poorly characterized. Here, we employ coarse-grained molecular dynamics simulations within an active learning protocol employing deep representational learning and Bayesian optimization to efficiently identify molecules capable of assembling pseudo-1D nanoaggregates with good stacking of the electronically active π-cores. We consider the DXXX-OPV3-XXXD oligopeptide family, where D is an Asp residue and OPV3 is an oligophenylenevinylene oligomer (1,4-distyrylbenzene), to identify the top performing XXX tripeptides within all 20 3 = 8000 possible sequences. By direct simulation of only 2.3% of this space, we identify molecules predicted to exhibit superior assembly relative to those reported in prior work. Spectral clustering of the top candidates reveals new design rules governing assembly. This work establishes new understanding of DXXX-OPV3-XXXD assembly, identifies promising new candidates for experimental testing, and presents a computational design platform that can be generically extended to other peptide-based and peptide-like systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Dark Energy Survey Year 1 Results: Wide-field mass maps via forward fitting in harmonic space

Herein we present new wide-field weak lensing mass maps for the Year 1 Dark Energy Survey (DES) data, generated via a forward fitting approach. This method of producing maps does not impose any prior constraints on the mass distribution to be reconstructed. The technique is found to improve the map reconstruction on the edges of the field compared to the conventional Kaiser–Squires method, which applies a direct inversion on the data; our approach is in good agreement with the previous direct approach in the central regions of the footprint. The mapping technique is assessed and verified with tests on simulations; together with the Kaiser–Squires method, the technique is then applied to data from the DES Year 1 data and the differences between the two methods are compared. We also produce the first DES measurements of the convergence Minkowski functionals and compare them to those measured in simulations.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

Modeling Savings for ENERGY STAR Smart Home Energy Management Systems

The objective of this study was to develop a repeatable and defensible methodology to analyze the energy savings for Home Energy Management Systems (HEMS) that meets the minimum requirements for certification under ENERGY STAR ® Smart Home Energy Management System (SHEMS) Version 1. Mandatory connected loads include a smart thermostat, two smart lights, and one smart power strip or smart outlet. Control strategies must include feedback to occupants through an in-home display, user programming, occupancy sensor-based controls, and responsiveness to utility signals such as demand response programs. Several occupant behavior patterns were selected to quantify the range of energy savings potential for a HEMS with this basic functionality. A literature review was conducted to establish realistic room-by-room occupancy levels and usage patterns for connected devices. A series of event-driven hourly profiles were created, followed by adjustments based on application of HEMS control strategies to thermostats, interior lighting, and plug load schedules. EnergyPlus modeling was performed using these hourly schedules in three locations (Boston, Houston, and Phoenix) to examine climate dependence of energy savings. Total site energy savings ranged from 4.3 to 27.1 MBtu/year (7%-35%), and utility bill savings ranged from $\$$123 to $\$$670/year (6%-29%). The highest predicted savings was realized by occupants that were not energy conscious prior to HEMS installation, but highly engaged with the HEMS controls once the system was installed. The smart thermostat accounted for most of the savings, followed by the smart power strip. Smart lighting did not save a significant amount of energy in our analysis, based on an assumption that efficient LEDs with no standby power would normally be installed anyway.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2021

Radiative dynamics and delayed emission in pure and doped InP/ZnSe/ZnS quantum dots

We have used time-correlated single photon counting to elucidate the radiative dynamics of InP/ZnSe/ZnS core/shell/shell quantum dots (QDs) that differ in the amount and distribution of excess indium. Stoichiometric QDs having an In:P atom ratio very near unity exhibit simple luminescence kinetics. The photoluminescence (PL) rises with the 40 ps instrument response function and exhibits a decay that is close to a single exponential with a time constant that decreases from 32 to 28 ns with increasing shell thickness. QDs having excess indium (In:P ratio of 1.15–1.63) show a significant component of a slower rise time assigned to transient population of indium-based hole traps in the ZnSe shell. They also have a slower PL decay, attributed to an equilibrium between these traps, which are optically dark, and the emissive valence-band state. This results in a radiative lifetime that increases from 32 to 48 ns with increasing shell thickness. Different treatments of the InP cores prior to shell deposition result in different core/shell interfaces as indicated by resonance Raman spectroscopy, as well as differences in the amplitude and timescale of the slow PL rise and the PL decay time. These are interpreted in terms of different radial distributions of the indium-based hole traps, which can be related to differences in the interfacial lattice strain.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

MultiGreen: A multiplexing architecture for GreenGate cloning

Genetic modification of plants fundamentally relies upon customized vector designs. The ever-increasing complexity of transgenic constructs has led to increased adoption of modular cloning systems for their ease of use, cost effectiveness, and rapid prototyping. GreenGate is a modular cloning system catered specifically to designing bespoke, single transcriptional unit vectors for plant transformation—which is also its greatest flaw. MultiGreen seeks to address GreenGate’s limitations while maintaining the syntax of the original GreenGate kit. The primary limitations MultiGreen addresses are 1) multiplexing in series, 2) multiplexing in parallel, and 3) repeated cycling of transcriptional unit assembly through binary intermediates. MultiGreen efficiently concatenates bespoke transcriptional units using an additional suite of level 1acceptor vectors which serve as an assembly point for individual transcriptional units prior to final, level 2, condensation of multiple transcriptional units. Assembly with MultiGreen level 1 vectors scales at a maximal rate of 2*⌈ log 6 n ⌉+3 days per assembly, where n represents the number of transcriptional units. Further, MultiGreen level 1 acceptor vectors are binary vectors and can be used directly for plant transformation to further maximize prototyping speed. MultiGreen is a 1:1 expansion of the original GreenGate architecture’s grammar and has been demonstrated to efficiently assemble plasmids with multiple transcriptional units. MultiGreen has been validated by using a truncated violacein operon from Chromobacterium violaceum in bacteria and by deconstructing the RUBY reporter for in planta functional validation. MultiGreen currently supports many of our in-house multi transcriptional unit assemblies and will be a valuable strategy for more complex cloning projects.

Science & Technology - Other Topics↗

DOE OSTI2022

Hicks Pries, Sulman, et al 2018. In situ incubation of 13C-labeled litter in three soil pits at Blodgett Forest, CA from 2013 to 2016

Even though over half of the world's soil organic carbon (SOC) is stored in subsoils (>20 cm deep), and the old ages of subsoil OC indicate its cycling differs from surface SOC, there are few studies examining in situ decomposition processes in deep soils. The purpose of this dataset is to elucidate these processes. In this study we added 13C-labeled fine roots to 15, 55, and 95 cm depths of a well-characterized coniferous forest Alfisol and monitored the amount of root-derived C remaining over 6, 12, and 30 months. We recovered the root-derived C in microbial phospholipid fatty acids (PLFAs) after 6 months and in coarse (>2 mm) particulate, fine (<2 mm) particulate, and dense, mineral-associated pools after 6, 12, and 30 months. Overall, root decomposition in the first 6 months was similar among all depths but significantly diverged at 30 months with faster decomposition at 15 cm than at 95 cm. There were more fungal and Gram negative-associated PLFAs at 15 cm than at 95 cm, and 13C analysis revealed those microbial groups preferred the added root carbon to native SOC. Mineral-associations were not the cause of slower decomposition at depth because similar amounts of applied root C was recovered in the dense fraction at all depths. The largest difference among depths was in the amount of root C recovered in the coarse particulate fraction, which was greater at 95 cm (50%) than at 15 cm (15%). Slower decomposition of the particulate pool at depth likely contributed to the increase in C:N ratios and depletion of δ13C values below 60 cm depth in our soil profiles. Simulations of these soils using the CORPSE model, which incorporates microbial priming effects and mineral stabilization of SOC, reproduced patterns of particulate and mineral-associated SOC over both time and depth and suggested that a lack of priming by root exudates at depth could account for the slower decomposition rate of particulate root material. Decomposition of deep particulate SOC may increase if root exudation or dissolved OC transport to depth increases.This dataset includes (1) characterization of the soil pits prior to litter addition, (2) characterization of how the litter changed after 6, 12, and 30 months, and (3) the CORPSE model structure, parameterization, and output.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2026

Acclimation to high and low diurnal light is flexible in Chlamydomonas reinhardtii

Chlamydomonas acclimates to repeated low (LL) or high light (HL) days by changing the abundance of photosynthetic complexes and the ultrastructure of its thylakoid membranes. These phenotypes persist through the night phases, suggesting a readiness for the daylight environment that is routinely experienced despite the intervening dark periods [S. Dupuis et al., Plant Cell 37, koaf086 (2025), 10.1093/plcell/koaf086]. Here, we investigate how prior acclimation impacts algal fitness upon a change in daylight intensity and how quickly Chlamydomonas can reprogram its photoprotective strategy in a diurnal context. We performed a systems analysis of synchronized populations acclimated to diurnal LL when subjected to HL days and of populations acclimated to diurnal HL when subjected to LL days. In the latter case, diurnal photoacclimation decreased fitness during the first day at a new light intensity: HL-acclimated cells barely increased in size over the first LL period, and they failed to complete a cell cycle. However, although LL-acclimated cells showed severe photodamage after 6 h of HL, they recovered chloroplast form and function later that afternoon and successfully divided at nightfall. These cells rapidly altered their thylakoid membrane ultrastructure, increased their photoprotective quenching capacity, and decreased their inventory of photosystem and antenna proteins by the end of the first HL day. Transcriptomic and proteomic analyses revealed rapid induction of thousands of genes, including those encoding proteases, chaperones, and other proteins involved in the chloroplast unfolded protein response. These results show that the alga is highly flexible and competent to rapidly acclimate to changes in diurnal light intensity.

Diel↗

DOE OSTI2025

Low Peripheral Blood Counts and Elevated Proinflammatory Cytokines Signal a Poor CD19 Chimeric Antigen Receptor T-cell Response in Acute Lymphoblastic Leukemia

CD19 chimeric antigen receptor T-cell (CAR-T) therapy has significantly improved outcomes for patients with relapsed/refractory B-cell acute lymphoblastic leukemia (R/R B-ALL). However, approximately 20% of patients fail to achieve a complete remission (CR), and some develop severe, life-threatening toxicities. Understanding the biological mechanisms underlying both dysfunctional responses and severe toxicity is essential for optimizing patient management and improving therapeutic efficacy. This study aimed to (1) characterize cytokine profiles associated with dysfunctional responses and severe toxicity following CAR-T infusion, (2) examine the timing and trajectory of cytokine changes in relation to treatment outcomes, and evaluate potential strategies for mitigating toxicity and treatment failure. We conducted a comprehensive analysis of serum cytokine profiles in 86 adult and pediatric patients undergoing autologous CD19 CAR-T therapy for B-ALL. Patients were categorized into three groups: (1) Dysfunctional response—Patients who failed to achieve a minimal residual disease-negative CR (MRD-CR) by Day 63 or who experienced recurrence of CD19+ disease in the setting ongoing CAR-T cell detection before Day 63. (2) Functional response with severe cytokine release syndrome (CRS) and/or neurotoxicity (NTX)—Patients with best response of MRD-CR by Day 63 who experienced grade 3 or higher CRS or NTX. (3) Functional response without severe CRS or NTX—Patients with best response of MRD-CR by Day 63 who did not experience grade =3 CRS or NTX. Cytokine levels were measured during the first-week postinfusion and correlated with treatment efficacy, toxicity outcomes, complete blood counts, and CAR-T expansion dynamics. This analysis aimed to better understand how cytokine profiles relate to patient outcomes and immune responses in CAR-T therapy. Patients with dysfunctional response exhibited decreased neutrophils, platelets, and levels of granulocytic cytokines (suggestive of low bone marrow reserve) alongside elevated pro-inflammatory cytokines by Day 1. Functional response with severe toxicity patients showed a progressive rise in proinflammatory cytokines, reaching similar levels to dysfunctional response patients by Day 7. We observed that high cytokines at both the Day 1 and Day 7 time points were associated with poor survival. These findings remained significant when adjusting for high disease burden, a known predictor of severe inflammatory toxicity and lack of response. Early post-CAR-T infusion inflammation is associated with both dysfunctional response and severe toxicity—even after adjusting for disease burden. This suggests that inflammation, in addition to disease burden, plays a role in determining patient outcome. Therefore, strategies aimed at reducing the pro-inflammatory state prior to or early after CAR-T cell infusion may improve outcomes for R/R B-ALL patients.

Serum cytokines↗

DOE OSTI2022

Chemical and Structural Alterations in the Amorphous Structure of Obsidian due to Nanolites

Obsidian is volcanic glass that results from the rapid cooling of silica-rich melt. Nanoscale crystallites precipitate out of the melt prior to solidification and remain embedded in the amorphous matrix. These crystallites provide information on the flow kinetics and composition of the melt. Due to the sparsity and size of nanolites, studies often focus on supramicron crystallites. This research takes advantage of the conchoidal fracture of obsidian by knapping samples with nanometer-thin edges for transmission electron microscopy characterization. Nanolites in the amorphous matrix are studied using energy-dispersive spectroscopy (EDS) and electron diffraction. Certain alkali and alkaline-earth cations exhibit patterns of depletion near Fe-oxide nanolites. EDS is used to identify nanolites and variations in the composition of the matrix. Parallel beam diffraction and radial distribution function analysis of nearest-neighbor distances determine average bond lengths in the matrix near nanolites, showing that nanolites influence the nearby short-range ordering and atomic character of the matrix. Analysis reveals decreased mean nearest-neighbor distances in the matrix adjacent to nanolites compared to the bulk. Our methods exhibit the required sensitivity to detect variations in the composition and structure near nanolites, and our findings indicate that obsidian nanolites contribute to quantifiable localized changes in the amorphous structure.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Strain-induced majority carrier inversion in ferromagnetic epitaxial LaCo O 3 - δ thin films

Tensile-strained LaCoO 3-δ thin films are ferromagnetic, in sharp contrast to the zero-spin bulk, although no clear consensus has emerged as to the origin of this phenomenon. While magnetism has been heavily studied, relatively little attention has been paid to electronic transport, due to the insulating nature of the strain-stabilized ferromagnetic state. Here, structure, magnetism, and transport are studied in epitaxial LaCoO 3-δ films (10–22-nm thick) on various substrates (from 1.4% compressive to 2.5% tensile strain), using synchrotron x-ray diffraction, scanning probe and transmission electron microscopy, magnetometry, polarized neutron reflectometry, resistivity, and Hall effect. High quality, smooth films are obtained, exhibiting superstructures associated with both oxygen vacancy ordering and periodic in-plane ferroelastic domains. Consistent with prior work, ferromagnetism with an approximately 80–85 K Curie temperature is observed under tension; polarized neutron reflectometry confirms a relatively uniform magnetization depth profile, albeit with interfacial dead layer formation. Electrical transport is found to have similar semiconducting nature to bulk, but with reduced resistivity and activation energy. Hall effect measurements, however, reveal a striking inversion of the majority carrier type, from p-type in the bulk and under compression to n-type under tension. While thus far overlooked, ferromagnetism in epitaxial LaCoO 3-δ films is thus directly correlated with n-type behavior, providing important insight into the ferromagnetic state in this system. Here, aided by density functional theory calculations, these results are interpreted in terms of tensile-strain-induced orbital occupation and band structure changes, including a rapid decrease in effective mass at the e g -derived conduction band minimum, and corresponding increase at the valence band maximum.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Reducing the matrix effect in mass spectral imaging of biofilms using flow-cell culture

The interactions between soil microorganisms and soil minerals play a crucial role in the formation and evolution of minerals and the stability of soil aggregates. Due to the heterogeneity and diversity of the soil environment, the under-standing of the functions of bacterial biofilms in soil minerals at the microscale is limited. A soil mineral-bacterial biofilm system was used as a model in this study, and it was analyzed by time-of-flight secondary ion mass spectrometry (ToF-SIMS) to acquire molecular level information. Static culture in multi-wells and dynamic flow-cell culture in microfluidics of biofilms were investigated. Our results show that more characteristic molecules of biofilms can be observed in SIMS spectra of the flow-cell culture. In contrast, biofilm signature peaks are buried under the mineral components in SIMS spectra in the static culture case. Spectral overlay was used in peak selection prior to performing Principal component analysis (PCA). Comparisons of the PCA results between the static and flow-cell culture show more pronounced molecular features and higher loadings of organic peaks of the dynamic cultured specimens. For example, fatty acids secreted from bacterial biofilm extracellular polymeric substance are likely to be responsible for biofilm dispersal due to mineral treatment up to 48 h. Such findings suggest that the use of microfluidic cells to dynamically culture biofilms be a more suitable method for reducing the matrix effect arisen from the growth medium and minerals as a perturbation fac-tor for improved spectral and multivariate analysis of complex mass spectral data in ToF-SIMS. These results show that the interaction mechanism between biofilms and soil minerals at the molecular level can be better studied using the flow-cell culture and advanced mass spectral imaging techniques like ToF-SIMS.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Templated synthesis enhances the cobalt adsorption capacity of a porous organic polymer

Divalent transition metals such as Co(II) are important targets for removal from water sources, due to their potential toxicity as well as their high value. In this study, we found that a series of porous organic polymers based on amide-linked tetraphenylmethane units are effective Co(II) ion adsorbents in aqueous solution. To increase the density of Co(II) binding sites, we then developed a templated synthesis in which the branched, rigid monomers are pre-assembled around Co(II) ions prior to polymerization. Furthermore, after polymer formation, the Co(II) template ions are removed to yield a material rich in Co(II) binding sites. Ion adsorption isotherms show that the Co(II)-templated material has an ion adsorption capacity significantly greater than those of the non-templated materials, highlighting the utility of a templated synthetic route. SEM and TEM images show the morphology of the templated polymer to be dramatically different from the non-templated polymers and to be similar in size and shape to the Co(II)–monomer precursors, emphasizing the role of the template ions in directing the formation of the resulting polymer. This guest-templated approach requires no functionalization of the generic monomer and represents a promising synthetic route to high-capacity ion adsorbents for water purification and aqueous separations.

36 MATERIALS SCIENCE↗

DOE OSTI2021

From noise to information: The transfer function formalism for uncertainty quantification in reconstructing the nuclear density

The neutron distribution of neutron-rich nuclei provides critical information on the structure of finite nuclei and neutron stars. Parity violating experiments—such as PREX and CREX—provide a clean and largely model-independent determination of neutron densities. Such experiments, however, are challenging and expensive, which is why sound statistical arguments are required to maximize the information gained. We introduce a new framework, the “transfer function formalism,” aimed at uncertainty quantification, model selection, and experimental design in the context of neutron densities. The transfer functions (TFs) are built analytically by expressing the linear response of the objective function (e.g., χ2) to small perturbations of the data. Using the TF formalism, we are able to analyze the expected overall uncertainty—quantified in terms of bias and variance—of the mean square radius and interior density of 48 Ca and 208 Pb. Using relativistic mean field models as a proxy for the weak-charge density—and assuming that a total of five measurements could be performed on the weak form factor of 48 Ca and 208 Pb—we identify the optimal models and experimental locations that minimize the uncertainty in the extraction of the radius and interior density. We also explore the use of the TF formalism to understand the influence of prior distributions for the model parameters, as well as the optimization of model hyperparameters not constrained by the data. Here,we establish how the choice of experimental locations and the model that is used can have a significant impact on the final uncertainties of the extracted quantities of interest. For challenging experiments such as CREX and PREX, a proper quantification of such uncertainties is critical. We have demonstrated how the TF formalism provides several advantages for this type of analysis.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2020

Solvent depression of transition temperature to selectively stimulate actuation of shape memory polymer foams

An embodiment of the invention is a shape memory polymer (SMP) foam designed to recover its original shape through exposure to a solvent. Thermo-responsive SMPs are polymers designed to maintain a programmed secondary shape until heated above their transition temperature, upon which the polymer recovers its original, or primary, shape. The thermo-responsive SMP foam is programmed to its secondary shape prior to use, typically compression of the foam to a small volume, and remains in this programmed shape until exposed to a selected solvent such as dimethyl sulfoxide or ethyl alcohol. Upon exposure to the solvent, the transition temperature of the SMP foam decreases below the temperature of the environment and the SMP foam actuates to its primary shape. The SMP foam is tailored to actuate upon exposure to specific solvents while minimizing or preventing actuation when exposed to water or other solvents. This selective solvent actuation can be used to increase working time of a SMP foam device, that is, the time allowed for use of a device without undesired actuation, while maintaining functional SMP actuation. Solvent actuated SMP foams can be used in various applications including, but not limited to, treatment of aneurysms and arterio-venous malformations, tissue engineering, and wound healing.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Solvent depression of transition temperature to selectively stimulate actuation of shape memory polymer foams

Boyle, Anthony↗

DOE OSTI2020

Density-Matrix Based Extended Lagrangian Born–Oppenheimer Molecular Dynamics

Extended Lagrangian Born–Oppenheimer molecular dynamics [ Phys. Rev. Lett. 2008, 100, 123004] is presented for Hartree–Fock theory, where the extended electronic degrees of freedom are represented by a density matrix, including fractional occupation numbers at elevated electronic temperatures. In contrast to regular direct Born–Oppenheimer molecular dynamics simulations, no iterative self-consistent field optimization is required prior to the force evaluations. To sample regions of the potential energy landscape where the gap is small or vanishing, which leads to particular convergence problems in regular direct Born–Oppenheimer molecular dynamics simulations, an adaptive integration scheme for the extended electronic degrees of freedom is presented. The integration scheme is based on a tunable, low-rank approximation of a fourth-order kernel, K, that determines the metric tensor, T ≡ K T K, used in the extended harmonic oscillator of the Lagrangian that generates the dynamics of the electronic degrees of freedom. Here, the formulation and algorithms provide a general guide to implement extended Lagrangian Born–Oppenheimer molecular dynamics for quantum chemistry, density functional theory, and semiempirical methods using a density matrix formalism.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

All-Atom Molecular Dynamics Simulations of the Temperature Response of Densely Grafted Polyelectrolyte Brushes

Water acts a “good” solvent for most polyelectrolyte (PE) systems due to its high dielectric strength. Theoretical and experimental studies have highlighted some remarkably unconventional properties of strongly water-swollen PE brushes such as a reduction in brush height with increase in temperature. However, an understanding of the atomistic-scale response of the grafted PE chains as well as the local arrangement and organization of the brush-supported counterions and water molecules to changes in temperature is missing. Here, we conduct an all-atom molecular dynamics (MD) study to probe the influence of temperature on the “microstructure” of densely grafted, fully ionized polyacrylic acid (PAA) brushes. The atomistic insights obtained from our study shed light on prior experimental observations and elucidate the key role played by changes in hydrogen bonding network in dictating the hydrophilicity of the brushes. In addition to the temperature-mediated reduction in height of the PE brush layer, we study the influence of temperature on several properties of the brush-supported water molecules such as the number distribution and kinetics of water-water and water-PE hydrogen bonds (HBs), mass density, orientational tetrahedral order parameter, and the activation energy associated with the self-diffusion of water. The effect of temperature on the properties of the brush-trapped counterions is also quantified via changes in their solvation structure, counterion-water radial distribution function (RDF), and the translational mobility. Our findings unravel several hitherto unknown phenomena such as the enhanced propensity for attaining PE brush-induced “water-in-salt”-like scenarios at elevated temperatures, and an increase in the activation energy for the self-diffusion of water with increase in the degree of PE brush-induced nanoconfinement (from ‘bulk water’ to ‘sparsely grafted PE brushes’ to ‘densely grafted PE brushes’). Further, we anticipate that this work will be helpful in improving our predictive capabilities of the temperature response of water-swollen PE brushes, which would be instrumental for better design of several PE brush systems used in a multitude of applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗