Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Solutions”

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

On the molecular correlations that result in field-dependent conductivities in electrolyte solutions

Employing recent advances in response theory and nonequilibrium ensemble reweighting, we study the dynamic and static correlations that give rise to an electric field-dependent ionic conductivity in electrolyte solutions. We consider solutions modeled with both implicit and explicit solvents, with different dielectric properties, and at multiple concentrations. Implicit solvent models at low concentrations and small dielectric constants exhibit strongly field-dependent conductivities. We compare these results to Onsager-Wilson theory of the Wien effect, which provides a qualitatively consistent prediction at low concentrations and high static dielectric constants but is inconsistent away from these regimes. The origin of the discrepancy is found to be increased ion correlations under these conditions. Explicit solvent effects act to suppress nonlinear responses, yielding a weakly field-dependent conductivity over the range of physically realizable field strengths. By decomposing the relevant time correlation functions, we find that the insensitivity of the conductivity to the field results from the persistent frictional forces on the ions from the solvent. Our findings illustrate the utility of nonequilibrium response theory in rationalizing nonlinear transport behavior.

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

The “good,” the “bad,” and the “hidden” in neutron scattering and molecular dynamics of ionic aqueous solutions

We characterize a concentrated 7.3 m CaCl 2 solution, combining neutron diffraction with chloride isotopic substitution (Cl-NDIS) in null water and molecular dynamics (MD) simulations. We elucidate the solution structure, thermodynamic properties, and extent of ion pairing previously suggested as concentration-dependent and often not observed at lower concentrations. Our Cl-NDIS measurements designate the solvent-shared ion pairing as dominant and the contact ion pairing (CIP) as insignificant even under conditions close to the solubility limit. The MD models parameterized against neutron diffraction with calcium isotopic substitution (Ca-NDIS) overestimate CIP despite successfully reproducing most of the Cl-NDIS signal. This drawback originates from the fact that Ca 2+ –Cl - interactions were primarily “hidden” in the Ca-NDIS signal due to overlapping with Ca 2+ –O w and Ca 2+ –H w contributions to the total scattering. Contrary, MD models with moderate CIP and possessing generally good performance at high concentrations fail to reproduce the NDIS measurements accurately. Therefore, the electronic polarization, introduced in most of the recent MD models via scaling ionic charges, resolves some but not all parameterization drawbacks. We conclude that despite improving the quality of MD models “on average,” the question “which model is the best” has not been answered but replaced by the question “which model is better for a given research.” An overall “good” model can still be inappropriate or, in some instances, “bad” and, unfortunately, produce erroneous results. Additionally, the accurate interpretation of several NDIS datasets, complemented by MD simulations, can prevent such mistakes and help identify the strengths, weaknesses, and convenient applications for corresponding computational models.

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

Structure–transport correlations in Na 11 Sn 2 SbSe 12 and its sulfide solid solutions

We report a new ion conductor—Na 11 Sn 2 SbSe 12 —as a possible candidate for a solid catholyte in composite cathodes for all-solid state Na-batteries, which exhibits a room temperature ionic conductivity of 0.15 ± 0.03 mS cm −1 and an activation energy of 0.39 ± 0.02 eV. The sulfide solid solutions of Na 11 Sn 2 SbSe 12 , namely, Na 11 Sn 2 SbS 12−x Se x (x = 1 and 6), were also investigated through a combination of Rietveld refinement against powder x-ray diffraction data and electrochemical impedance spectroscopy to reveal the complex structure–property relationships governing ion transport in this class of materials. Meanwhile, broadening of the Na-ion diffusion pathways in Na 11 Sn 2 SbSe 12 is expected to facilitate Na-ion transport compared to the sulfide-rich member of the solid solution, the opposite holds: increasing the Se fraction in Na 11 Sn 2 SbS 12−x Se x leads to a little change in the activation energy but a reduction in the ionic conductivity. We ascribe this to the lowering of the prefactor, σ 0 , in the Arrhenius relationship with increasing lattice “softening” as a function of higher Se content.

Ramos, Erika P.↗
DOE OSTI2023

Investigation of defects in melt and solution grown perovskite CsPbBr 3 single crystals

Recently, there has been considerable interest in x-ray and gamma ray detectors with large volume and high energy resolution that operate at room temperature. To improve detector energy resolution, the carrier mobility-lifetime product needs to be increased, and the electronic trap state concentration needs to be minimized. Defect concentrations in the part per billion range can alter the charge transport and carrier recombination lifetime. In this work, thermally stimulated current spectroscopy measurements were systematically carried out in bulk halide perovskite single crystals of CsPbBr 3 over a temperature range of 80–320 K. The origins and trap parameters of CsPbBr 3 crystals from the solution growth and melt growth procedures were determined and compared. Trap concentrations were ranged from 1 × 10 11 to 1 × 10 16 cm –3 . Appreciable detector performance was observed for CsPbBr3 crystals with trap concentrations less than 1 × 10 14 cm –3 . The comparison of spectral responses of crystal samples grown using two different methods shows that, after purification, solution-grown crystals are comparable to melt-grown crystals in terms of low defect concentration and improved detector performance. Finally, for an improved mobility-lifetime product and enhanced spectral response to high energy radiation from fissile materials, trap states in either type of a crystal ingot must be reduced closer to 10 11 cm –3 .

36 MATERIALS SCIENCE↗
DOE OSTI2024

Training self-learning circuits for power-efficient solutions

As the size and ubiquity of artificial intelligence and computational machine learning models grow, the energy required to train and use them is rapidly becoming economically and environmentally unsustainable. Recent laboratory prototypes of self-learning electronic circuits, such as “physical learning machines,” open the door to analog hardware that directly employs physics to learn desired functions from examples at a low energy cost. In this work, we show that this hardware platform allows for an even further reduction in energy consumption by using good initial conditions and a new learning algorithm. Using analytical calculations, simulations, and experiments, we show that a trade-off emerges when learning dynamics attempt to minimize both the error and the power consumption of the solution—greater power reductions can be achieved at the cost of decreasing solution accuracy. Finally, we demonstrate a practical procedure to weigh the relative importance of error and power minimization, improving the power efficiency given a specific tolerance to error.

Stern, Menachem (ORCID:0000000302158082)↗
DOE OSTI2024

Dynamic and relativistic effects on Pt–Pt indirect spin–spin coupling in aqueous solution studied by ab initio molecular dynamics and two- vs four-component density functional NMR calculations

Treating 195Pt nuclear magnetic resonance parameters in solution remains a considerable challenge from a quantum chemistry point of view, requiring a high level of theory that simultaneously takes into account the relativistic effects, the dynamic treatment of the solvent–solute system, and the dynamic electron correlation. A combination of Car-Parrinello molecular dynamics (CPMD) and relativistic calculations based on two-component zeroth order regular approximation spin–orbit Kohn–Sham (2c-ZKS) and four-component Dirac–Kohn–Sham (4c-DKS) Hamiltonians is performed to address the solvent effect (water) on the conformational changes and JPtPt1 coupling. A series of bridged PtIII dinuclear complexes [L1–Pt2(NH3)4(Am)2–L2]n+ (Am = α–pyrrolidonate and pivalamidate; L = H2O, Cl−, and Br−) are studied. The computed Pt–Pt coupling is strongly dependent on the conformational dynamics of the complexes, which, in turn, is correlated with the trans influence among axial ligands and with the angle N–C–O from the bridging ligands. The J-coupling is decomposed in terms of dynamic contributions. The decomposition reveals that the vibrational and explicit solvation contributions reduce JPtPt1 of diaquo complexes (L1 = L2 = H2O) in comparison to the static gas-phase magnitude, whereas the implicit solvation and bulk contributions correspond to an increase in JPtPt1 in dihalo (L1 = L2 = X−) and aquahalo (L1 = H2O; L2 = X−) complexes. Relativistic treatment combined with CPMD shows that the 2c-ZKS Hamiltonian performs as well as 4c-DKS for the JPtPt1 coupling.

Chemistry↗
DOE OSTI2024

Nonlinear Poisson–Boltzmann solutions for charged parallel plates: When opposite charges repel

I present an exact solution of the Poisson–Boltzmann equation for two parallel plates and discuss the solution properties. I discuss in more detail plates with opposite charges: In this case, there are two critical separations, L c,1 < L c,2 . For separations less than L c,1 , the force between plates is repulsive. It switches to attractive at L c,1 , but with the electric potential having the same sign on both plates. For L > L c,2 , the force remains attractive, and the potential at the plates has the same sign as the charge on each plate. I also describe charge regulation, determined by pK a , and provide formulas for both the critical distance where oppositely charged plates repel and their charging process. Finally, the implications of these results for the nanoparticle assembly, as driven by electrostatic interactions, are also discussed.

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

Modeling ethanol/water adsorption in all-silica zeolites using the real adsorbed solution theory

A comprehensive set of single-component and binary isotherms were collected for ethanol/water adsorption into the siliceous forms of 185 known zeolites using grand-canonical Monte Carlo simulations. Using these data, a systematic analysis of ideal/real adsorbed-solution theory (IAST/RAST) was conducted and activity coefficients were derived for ethanol/water mixtures adsorbed in different zeolites based on RAST. It was found that activity coefficients of ethanol are close to unity while activity coefficients of water are larger in most zeolites, indicating a positive excess free energy of the mixture. This observation can be attributed to water/ethanol interactions being less favorable than water/water interactions in the single-component adsorption of water at comparable loadings. The deviation from ideal behavior can be highly structure-dependent but no clear correlation with pore diameters was identified. Furthermore, our analysis also demonstrates the following: (1) accurate unary isotherms in the low-loading regime are critical for obtaining physically sensible activity coefficients; (2) the global regression scheme to solve for activity model parameters performs better than fitting activity models to activity coefficients calculated locally at each binary state point; and (3) including the dependence on adsorption potential offers only a minor benefit for describing binary adsorption at the lowest fugacities. Finally, the Margules activity model was found incapable of capturing the non-ideal adsorption behavior over the entire range of fugacities and compositions in all zeolites, but for conditions typical of solution-phase adsorption, RAST predictions using zeolite-specific or even bulk Margules parameters provide an improved description compared to IAST.

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

Cesium exhibits different mesoscale segregation and ion pairing than lithium, sodium, or potassium in concentrated alkaline aqueous nitrite solutions

Understanding ion pairing in concentrated alkaline electrolytes facilitates the prediction and control of chemical processes in nuclear waste. While sodium is the dominant alkali metal in such systems, cesium often exhibits distinct bulk behavior relative to lighter alkali cations. Aqueous mixtures of cesium hydroxide and sodium nitrite were compared to alkali hydroxide analogs using small-angle x-ray scattering, revealing that cesium disrupts the sodium nitrite electrolyte structure, forming cesium-rich domains. This mesoscale segregation contrasts with the near-ideal mixing observed in other mixed hydroxide–nitrite systems. Raman spectroscopy indicates cesium–nitrite ion pairing, evidenced by vibrational shifts distinct from those associated with sodium. Multinuclear magnetic resonance spectroscopy further supports cesium–nitrite and cesium–hydroxide association, revealing a distinct local environment for nitrite in cesium-containing solutions. Together, these findings show that cesium promotes a unique solution structure dominated by specific ion pairing within segregated domains. This structural organization may influence radical generation pathways in high-ionic-strength alkaline media relevant to nuclear waste processing and management.

Nienhuis, Emily T. [Pacific Northwest National Lab↗
DOE OSTI2026

Metal–bicarbonate ion pairing in alkaline aqueous solutions from multilevel embedded correlated wavefunction theory and molecular dynamics

In this work, we examine ion-pairing mechanisms of bicarbonates in alkaline aqueous solutions with the divalent metal ions most abundantly present in seawater, namely, Ca 2+ and Mg 2+ . We employ a rare-event enhanced sampling approach within first-principles molecular dynamics to explore regions of phase space spanning solvent-shared to contact ion pairs. Second-order Møller–Plesset perturbation theory (MP2) corrections are subsequently applied in an embedding framework (EMB) to refine the electronic structure of stationary states and associated reaction barriers along the free-energy profiles while retaining the extended solvent effects at the density functional theory (DFT) level. Ca 2+ –HCO 3 − was previously hypothesized to exist in a solvent-shared ion pair (SSHIP) by DFT studies with an endergonic contact ion pair (CIP) formation; however, our EMB-MP2 refinement of the DFT ion-pairing pathways reveals that Ca 2+ and HCO 3 − form a virtually barrier-free CIP in alkaline solutions, with even more energetic ease than the widely studied Ca–CO 3 ion pair. We find qualitative agreement between DFT and EMB-MP2 for Mg 2+ —unlike Ca 2+ , Mg 2+ refuses to shed its strong hydration shell, thereby preferring a SSHIP state with a significant activation barrier to crossover to the CIP forms—a trait reminiscent of ion pairing in Mg–CO 3 and closely related to the kinetic limitations underlying the famous subject of the dolomite problem. Our study highlights the importance of improved electronic structure descriptions of liquids, modeled as a condensed phase of matter lacking in long-range crystalline order. It also strongly suggests that Ca 2+ –HCO 3 − CIPs are likely precursors involved in prenucleation of CaCO 3 mineral formation in seawater.

Sharma, Vidushi [Princeton Plasma Physics Laborato↗
DOE OSTI2020

Experimental and theoretical evidence for hydrogen doping in polymer solution-processed indium gallium oxide

The field-effect electron mobility of aqueous solution-processed indium gallium oxide (IGO) thin-film transistors (TFTs) is significantly enhanced by polyvinyl alcohol (PVA) addition to the precursor solution, a >70-fold increase to 7.9 cm 2 /Vs. To understand the origin of this remarkable phenomenon, microstructure, electronic structure, and charge transport of IGO:PVA film are investigated by a battery of experimental and theoretical techniques, including In K-edge and Ga K-edge extended X-ray absorption fine structure (EXAFS); resonant soft X-ray scattering (R-SoXS); ultraviolet photoelectron spectroscopy (UPS); Fourier transform-infrared (FT-IR) spectroscopy; time-of-flight secondary-ion mass spectrometry (ToF-SIMS); composition-/processing-dependent TFT properties; high-resolution solid-state 1 H, 71 Ga, and 115 In NMR spectroscopy; and discrete Fourier transform (DFT) analysis with ab initio molecular dynamics (MD) liquid-quench simulations. The 71 Ga{ 1 H} rotational-echo double-resonance (REDOR) NMR and other data indicate that PVA achieves optimal H doping with a Ga···H distance of ~3.4 Å and conversion from six- to four-coordinate Ga, which together suppress deep trap defect localization. This reduces metal-oxide polyhedral distortion, thereby increasing the electron mobility. Finally, hydroxyl polymer doping thus offers a pathway for efficient H doping in green solvent-processed metal oxide films and the promise of high-performance, ultra-stable metal oxide semiconductor electronics with simple binary compositions.

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

Enhanced polymorph metastability drives glycine nucleation in aqueous salt solutions

Crystal nucleation from aqueous solutions influences countless geological, biochemical, astrophysical, environmental, and materials science–related phenomena, including ice formation, the manufacturing of active pharmaceutical ingredients, development of diseases such as Alzheimer’s and the origin of life itself. Understanding and controlling nucleation is essential for designing materials with specific properties, developing strategies to inhibit or promote crystallization in various contexts and preventing pathological aggregation in neurodegenerative diseases. Similar to the protein structure prediction problem—where a single amino acid sequence can in theory adopt one most stable conformation but in practice may sample multiple competing conformations—crystal nucleation faces a parallel challenge: the same chemical species can form diverse polymorphs under different environmental conditions (e.g., temperature, pressure, solvent). Each polymorph presents its own set of physical and chemical properties, highlighting the importance of understanding and controlling polymorph selection in fields ranging from pharmaceuticals to materials design. Despite advances in experimental and computational methods for studying phase transitions and polymorph stability, nucleation remains challenging due to its nanoscale nature. Furthermore, in practical settings, salts and impurities can further influence crystal nucleation in diverse contexts, from scaling in pipelines and desalination plants to the durability of concrete and the efficiency of battery materials. This can lead to the formation of polymorphs that may differ from the most stable phase in pure solutions. Or, even though the final structure might appear same irrespective of whether the environment contained impurities or not, the mechanism through which it was formed might be completely different and not intuitive.

Wang, Ruiyu [University of Maryland, College Park,↗
DOE OSTI2025

Separation of dilute plutonium(IV) from nitric acid solutions via Reillex HPQ anion exchange resin

Anion exchange is a common method for the separation of plutonium (Pu) in aqueous systems, providing high recovery and decontamination from impurities. Reillex HPQ anion exchange resin has been used around the Department of Energy complex in this application for decades but has typically been used for Pu feed concentrations greater than 1 g/L. It is desirable to recover Pu at more dilute concentrations from a dissolved used fuel or target element, for example. In this work, the performance of Reillex HPQ anion exchange resin was characterized with feed solutions of ca. 0.24 g/L Pu. A significant reduction in the performance of the resin was realized with this dilute concentration relative to a feed solution of 2.35 g/L Pu with similar feed flow rates and the same resin column, likely due to the increased driving force for mass transfer with a higher feed concentration. The maximum capacity of Reillex HPQ resin for Pu can be improved by reducing the flow rate and increasing the feed Pu concentration. Furthermore, the breakthrough point of Pu can be improved in the same manner and by increasing the resin bed height. Diffusion of Pu into the resin was illustrated to be the rate limiting step.

Plutonium separations↗
DOE OSTI2022

Benchmark of few-level quantum theory vs ab initio numerical solutions for the strong-field Autler–Townes effect in photoionization of hydrogen

Abstract The temporal and spectral consequences of an intermediate resonance en route to photoionization are investigated theoretically in two ways: by solving few-level model equations and by ab initio numerical solution of the time-dependent Schrödinger equation, in both cases for hydrogen in three dimensions. The model consists of atomic states resonantly field-dressed in a three-level reduction of the hydrogen atom that consists of the 2 p –3 d (Balmer) transition and one energetically-distant continuum state. The model’s level occupation probabilities are derived from three Schrödinger amplitude equations and are benchmarked against an ab initio numerical solution for the hydrogen electron’s wavefunction under the same field. We examine contrasts between the results of the two approaches with a particular focus on Autler–Townes doublets that appear in the photoelectron spectrum.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

Mathematical solutions in internal dose assessment: A comparison of Python-based differential equation solvers in biokinetic modeling

Abstract In biokinetic modeling systems employed for radiation protection, biological retention and excretion have been modeled as a series of discretized compartments representing the organs and tissues of the human body. Fractional retention and excretion in these organ and tissue systems have been mathematically governed by a series of coupled first-order ordinary differential equations (ODEs). The coupled ODE systems comprising the biokinetic models are usually stiff due to the severe difference between rapid and slow transfers between compartments. In this study, the capabilities of solving a complex coupled system of ODEs for biokinetic modeling were evaluated by comparing different Python programming language solvers and solving methods with the motivation of establishing a framework that enables multi-level analysis. The stability of the solvers was analyzed to select the best performers for solving the biokinetic problems. A Python-based linear algebraic method was also explored to examine how the numerical methods deviated from an analytical or semi-analytical method. Results demonstrated that customized implicit methods resulted in an enhanced stable solution for the inhaled 60 Co (Type M) and 131 I (Type F) exposure scenarios for the inhalation pathway of the International Commission on Radiological Protection (ICRP) Publication 130 Human Respiratory Tract Model (HRTM). The customized implementation of the Python-based implicit solvers resulted in approximately consistent solutions with the Python-based matrix exponential method ( expm ). The differences generally observed between the implicit solvers and expm are attributable to numerical precision and the order of numerical approximation of the numerical solvers. This study provides the first analysis of a list of Python ODE solvers and methods by comparing their usage for solving biokinetic models using the ICRP Publication 130 HRTM and provides a framework for the selection of the most appropriate ODE solvers and methods in Python language to implement for modeling the distribution of internal radioactivity.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2020

An efficient approach to explore the solution space of a wind turbine rotor design process

This study proposes a novel approach to explore the solution space of a wind turbine rotor design process. The goal is to offer to blade designers the possibility to select an optimal rotor for given market conditions, assessing trade-offs and alternatives in a matter of minutes. The design process consists of sequential aerodynamic-structural optimizations, where the two loops are linked via a blade-loading parameter that can be varied by the user. A design study is performed starting from the IEA Wind Task 37 land-based reference wind turbine. The solution space is characterized for rotor diameters in the range of 130–160 m, rated generator power values in the range of 3.0–6.0 MW, and tip-speed ratios in the range of 7.5–12.5. Results are discussed highlighting optimal design decisions and trade-offs.

17 WIND ENERGY↗
DOE OSTI2020

Stability and response of trapped solitary wave solutions of coupled nonlinear Schrödinger equations in an external, $\mathcal{PT}$- and supersymmetric potential

In this work, we present trapped solitary wave solutions of a coupled nonlinear Schrödinger (NLS) system in 1 + 1 dimensions in the presence of an external, supersymmetric and complex $\mathcal{PT}$-symmetric potential. The Schrödinger system this work focuses on possesses exact solutions whose existence, stability, and spatio-temporal dynamics are investigated by means of analytical and numerical methods. Two different variational approximations are considered where the stability and dynamics of the solitary waves are explored in terms of eight and twelve time-dependent collective coordinates (CCs). We find regions of stability for specific potential choices as well as analytic expressions for the small oscillation frequencies in the CC approximation. Our findings are further supported by performing systematic numerical simulations of the NLS system.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Atomistic structure of the SARS-CoV-2 pseudoknot in solution from SAXS-driven molecular dynamics

SARS-CoV-2 depends on -1 programmed ribosomal frameshifting (-1 PRF) to express proteins essential for its replication. The RNA pseudoknot stimulating -1 PRF is thus an attractive drug target. However, the structural models of this pseudoknot obtained from cryo-EM and crystallography differ in some important features, leaving the pseudoknot structure unclear. We measured the solution structure of the pseudoknot using small-angle X-ray scattering (SAXS). The measured profile did not agree with profiles computed from the previously solved structures. Beginning with each of these solved structures, we used the SAXS data to direct all atom molecular dynamics (MD) simulations to improve the agreement in profiles. In all cases, this refinement resulted in a bent conformation that more closely resembled the cryo-EM structures than the crystal structure. Applying the same approach to a point mutant abolishing -1 PRF revealed a notably more bent structure with reoriented helices. This work clarifies the dynamic structures of the SARS-CoV-2 pseudoknot in solution.

60 APPLIED LIFE SCIENCES↗