Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Correlation coefficient”

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 559 records · Page 31
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

The calculation of light element impurity (α, n) yield curves in a PuO 2 matrix and associated specific yield coefficients: Influence of the reaction cross sections

Most of the Pu separated from irradiated commercial nuclear fuel is stored as PuO 2 . The primary quantitative nondestructive measurement technique used to verify the amount of Pu in storage containers is passive neutron correlation counting. An important physical property of the oxide material is the ratio, α, of the rate of (α, n) neutrons produced inside the item to the rate of neutrons produced by spontaneous fission. This ratio influences the precision of the correlated counting method and affects the interpretation of the data because of how it changes both the primary total neutron production rate and the rate of induced fission events taking place inside the item. In addition to the main O(α, n) contribution, additional contributions come from α-particle interactions with light element impurities that are inevitably present. In this work, we calculate specific (α, n) yield coefficients, expressed in units of neutrons per second per gram of α-emitting nuclide per part per million by mass of the specified impurity element distributed in a pure PuO 2 matrix, for some key α-emitting actinides commonly present in reprocessed Pu ( 238–242 Pu+ 241 Am). These coefficients are directly applicable to nuclear safeguards verification work in which the α ratio is often calculated from the Pu-isotopic composition and chemical information obtained by other means. Further, they also provide a convenient up-to-date reference set against which values generated by other methods can be compared. Results are presented for impurities with atomic number from 3 to 17 inclusive, plus K and Fe. In most cases, these coefficients are not expected to change by more than 5%–10% at any time in the future. However, as new data become available, changes as large as 20% may be needed for some targets (e.g., F). The present yield calculations are limited by the general shortage of quality experimental total (α, n) reaction cross section data, which, together with unexplained variation between determinations, means that an objective and coherent evaluation is not possible. The situation is even less satisfactory for the partial differential cross section needed to calculate neutron spectra.

(𝛼, n) reactions↗
DOE OSTI2025

Pool boiling heat transfer evaluation of next-generation dielectric fluid: Opteon™ 2P50

The growing use of artificial intelligence has led to heavy thermal loads and high heat dissipation rates in data centers. Conventional air-cooled technologies are not able to fulfill these requirements. To overcome these challenges, two-phase immersion cooling (2PIC) has emerged as one of the leading technologies for high power-density chips. 2PIC increases the heat dissipation rate and efficiency of the system while reducing the footprint of the cooling equipment. A fluid with adequate dielectric properties, a suitable normal boiling temperature to maintain chip temperatures, and good material compatibility, is desired for 2PIC system. In this study, the pool boiling heat transfer of a new developmental dielectric fluid, Opteon™ 2P50, was experimentally investigated. The heat transfer coefficients at various heat fluxes (20–150 kW/m 2 ) and the critical heat flux were measured using a smooth aluminum surface. Compared with HFE-7100, Opteon™ 2P50 shows higher heat transfer coefficient (up to 59% higher) and a slightly lower value of critical heat flux (around 5.9% lower). The modified Cooper correlation with the optimized leading constant resulted in reliable prediction accuracy with a 5.3% mean absolute error percentage. Overall, these results indicate that the new dielectric fluid provides similar thermal performance to some legacy fluids.

2P50↗
DOE OSTI2023

Separation of volatile organic contaminants from water using a direct-contact dehumidifier: An experimental study and modeling

A packed-bed, direct-contact dehumidifier (PB-DCD) has been developed to selectively condense water from a gaseous mixture of organic contaminants. In this process, hot air is humidified by mixing with wastewater vapor in a static mixer and is dehumidified by subcooled fresh water passing through the PB-DCD in a counter-current configuration. The goal of this study is to understand the relationship between various operating parameters and how they affect separation of isopropyl alcohol (IPA) from the gaseous mixture with the purpose of expanding the investigation to other contaminants in the future. The process has been successfully used to separate IPA from wastewater. Different parameters including the concentration of IPA in water, feed temperature, air-to-vapor ratio, and cooling-water-to-vapor mass ratio, as well as residence time were examined to determine their influence on contaminant separation and clean water recovery rate. The Buckingham Pi theorem was applied to reduce the number of studying variables by generating non-dimensional groups. Experiments were performed to determine the relationship between the non-dimensional parameters and model mass transfer in the system. Mass transfer coefficients for multicomponent gaseous mixture condensation in a packed-bed direct contact condenser column are calculated and a modified correlation of mass transfer is developed. The parameters of the mass transfer model are obtained from operating conditions of the PB-DCD, stream conditions and thermodynamic properties. A modified Sherwood correlation is developed to predict IPA separation from a gaseous mixture consisting of air, water, and IPA. The modified Sherwood correlation is applied to successfully predict the Sherwood number with a mean absolute error equal to 6.6 %. According to the experimental results, the maximum IPA separation achieved was 76.3 %, which corresponds to a water recovery rate of 58.7 %.

42 ENGINEERING↗
DOE OSTI2020

Heat Transfer Coefficients of Additively Manufactured Tubes with Internal Pin Fins for Supercritical Carbon Dioxide Cycle Recuperators

This paper describes the measurement of convective heat transfer coefficients and friction factors for sCO 2 flowing in additively manufactured tubes with internal pin fins at the US DOE’s National Energy Technology Laboratory in Morgantown, WV. The measurement procedures were validated by conducting benchmark tests with smooth stainless-steel tube and comparing the results with published correlations for Nusselt number (Nu) and friction factor. Over Reynolds numbers (ReD) ranging from 5 × 10 4 to 2.5 × 10 5 , measured Nu was within 5% of the Dittus-Boelter correlation and measured friction factors were within 5% of the McAdams correlation for smooth tube flow. The candidate pin fin patterned pipes were additively manufactured (AM) at the Oak Ridge National Laboratory. The pins were circular or elliptical in cross-section and were printed at a 30⁰ angle relative to the inner wall (to meet AM constraints). The pin arrangement was helical to promote enhanced heat transfer due to swirl flow. Pin length to diameter aspect ratio was 1.33, 2, and 8, while the pin diameter to tube diameter ratio was 0.188, 0.125, and 0.063. Tests were performed for ReD varying from 6.9 × 10 4 to 2.2 × 10 5 and at conditions equivalent to the low pressure (LP) outlet (8.69 MPa, 361 K) and the high pressure (HP) inlet (20.7 MPa, 350 K) of the low temperature recuperator (LTR) in an indirect sCO 2 cycle. The Wilson plot technique was utilized to measure the bulk heat transfer coefficients. For the best performing design (tube A, pin length to tube diameter ratio: 1.33, pin diameter to tube diameter ratio: 0.19), the local heat transfer coefficient increased by 136% relative to the Dittus-Boelter correlation at the LTR low pressure outlet and 194% at the LTR high pressure inlet. These correspond to a 282% and a 271% increase in the product of the heat transfer coefficient and surface area (adjusted for fin efficiency) product, respectively. Large pressure drops across the test articles were observed. For Tube Design A, the average friction factor, across the range of ReD considered, was significantly larger than the McAdams correlation at both the LTR LP outlet and the LTR HP inlet. A thermal performance factor was utilized to express the ratio of material required to build a finned heat exchanger relative to a finless heat exchanger with the same heat duty and pumping power. Tube Design A was estimated to decrease the required heat exchanger material by 13%.

Searle, Matthew↗
DOE OSTI2021

Correlation between antiferromagnetic and Mott states in spin orbit coupled Sr 2 Ir O 4 : A study of Sr 2 Ir 1 – x M x O 4 ( M = Fe or Co)

Sr 2 IrO 4 is an archetypal spin-orbit-coupled Mott insulator with an antiferromagnetic state below 240 K. Here we report results of our study on single crystals of Sr 2 Ir 1–x Fe x O 4 (0 ≤ x < 0.32) and Sr 2 Ir 1–x Co x O 4 (0 ≤ x < 0.22). Fe doping retains the antiferromagnetic state but simultaneously precipitates an emergent metallic state whereas Co doping causes a rapid collapse of both the antiferromagnetic and Mott states, giving rise to a confined metallic state featuring a pronounced linearity of the basal-plane resistivity up to 700 K. The results indicate tetravalent Fe 4+ (3d 4 ) ions in the intermediate spin state with S = 1 and Co 4+ (3d 5 ) ions in the high spin state with S = 5/2 substituting for Ir 4+ (5d 5 ) ions in Sr 2 IrO 4 , respectively. The effective magnetic moment closely tracks the Néel temperature as doping increases, suggesting that the spin state of the dopant predominately determines the magnetic properties in doped Sr 2 IrO 4 . Furthermore, all relevant properties including charge-carrier density (e.g., 10 28 /m 3 ), Sommerfeld coefficient (e.g., 19 mJ / mole K 2 ), and Wilson ratio (e.g., 2.6) consistently demonstrate an emergent metallic state that is both robust and highly correlated in the two systems, arising from the percolation of bound states and the weakening of structural distortions. Furthermore, this study strongly indicates that the antiferromagnetic and Mott states are correlated in an unconventional manner in Sr 2 IrO 4 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

In‐mold rheology and automated process control for injection molding of recycled polypropylene

Abstract Manufacturing plastic parts with secondary feedstocks has risen to the forefront of importance in recent years. However, the variation in molecular weight and rheology of secondary feedstock can lead to inconsistent part quality. This work evaluates the effectiveness of a novel closed‐loop adaptive process control system that adjusts nozzle pressure in response to in‐mold pressure data. Five different recycled polypropylene blends, with a broad distribution of flow properties, were evaluated to determine the effectiveness of the control system at reducing processing variation. The experimental results show that the process control strategy reduced the variation within the mold, as seen by in‐mold pressure curves and calculated in‐mold viscosity values. Additionally, the parameters that control the automated process adjustments were investigated, showing the importance of optimization. The analysis of the correlation between in‐mold rheology and mechanical properties showed a slight variation in the mechanical properties and parts weight with a coefficient of variation of under 5%. Overall, the results demonstrate the ability of pressure‐controlled molding and automated viscosity adjustment to reduce the variability when molding a secondary feedstock. Highlights Pressure‐controlled injection molding of recycled polypropylene. Automated closed‐loop adaptive process control methodology. Methodology resulted in a reduction in pressure variation during molding. Changes in mechanical properties and in‐mold viscosity were investigated. Results show the potential of pressure‐controlled molding at reducing variation.

Krantz, Joshua↗
DOE OSTI2023

Quantitatively controlled electrophoretic deposition of nanocrystal films from non-aqueous suspensions

This study presents a novel method to correlate the mass and charge transfer kinetics during the electrophoretic deposition of nanocrystal films by using a purpose-built double quartz crystal microbalance combined with simultaneous current-measurement. Our data support a multistep process for film formation: generation of charged nanocrystal flux, charge transfer at the electrode, and polarization of neutral nanocrystals near the electrode surface. The polarized particles are then subject to dielectrophoretic forces that reduce diffusion away from the interface, generating a sufficiently high neutral particle concentration at the interface to form a film. The correlation of mass and charge transfer enables quantification of the nanocrystal charge, the fraction of charged nanocrystals, and the initial sticking coefficient of the particles. These quantities permit calculation of the film thickness, providing a theoretical basis for using concentration and voltage as process parameters to grow films of targeted thicknesses.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2020

Carrier recombination and diffusion in high-purity diamond after electron irradiation and annealing

Carrier-transport mechanisms are studied in high-purity diamond irradiated with 6 MeV electrons in the dose range of 10 12 –10 16 cm -2 and annealed at different temperatures up to 1450°C. Lifetimes and diffusion coefficients are extracted using two pump–probe techniques based on free-carrier absorption and transient-grating principles and then correlated with the corresponding defect evolution from spectroscopic measurements. Furthermore, the neutral monovacancy is revealed as the main carrier recombination center in the as-irradiated diamond, providing bipolar carrier lifetimes of a few nanoseconds at the highest irradiation dose. Carrier-capture cross sections are reduced during annealing as vacancies aggregate into divacancies at ≤1000°C and extended vacancy clusters at 1450°C.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Correlated Displacement of Dynamic Elastic Dipoles Produces Nonclassical Electrostriction in Zr-Doped Ceria

By combining experimental data with density functional theory-based ab initio molecular dynamics modeling, this work provides evidence that nonclassical electrostriction in isovalent Zr-doped ceria is due to the correlated anharmonic motion of dynamic elastic dipoles associated with multiple [ZrO 8 ]-local bonding units with a high Zr concentration (Zr 0.1 Ce 0.9 O 2 ). Introduction of 0.5 mol % trivalent or divalent codopants (Sc, Yb, La, or Ca) reduces the longitudinal electrostriction strain coefficient by more than a factor of 10, produces a 3-fold decrease in the relative dielectric permittivity, and increases the elastic modulus. Since these changes depend neither on the radius nor on the valency of the codopant, we conclude that the responsible species are charge-compensating oxygen vacancies (VO). For trivalent dopants (Do 0.005 Zr 0.1 Ce 0.895 O 1.9975 ), oxygen vacancies are present at a concentration ratio 1:40 with respect to Zr, giving, for random distribution, a characteristic interaction distance of ≤2.3 unit cells (1.2 nm). Oxygen vacancies participate in [ZrO 7 -V O ] local bonding units, disrupting the correlated dynamic displacements of the connected [ZrO 8 ]-local bonding units. Finally, such correlated motion of dynamic elastic dipoles may also explain the exponential increase in the longitudinal electrostriction strain coefficient with an increase in Zr concentration to <0.2 mole fraction and must be taken into account for further development of nonclassical electrostrictors based on Zr-doped ceria.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Relaxor behavior and electrothermal properties of Sn- and Nb-modified (Ba,Ca)TiO 3 Pb-free ferroelectric

Relaxor ferroelectrics have drawn attention for possible applications in solid-state cooling and thermal energy harvesting, owing to their electrothermal energy conversion properties. Here, we have synthesized and characterized the structure–property correlations of a new Sn- and Nb-doped (Ba,Ca)TiO 3 relaxor ferroelectric with large pyroelectric and electrocaloric effects over a broad temperature range. We observed two peaks for the temperature-dependent pyroelectric coefficient: (i) -(∂P/∂T) ~ 563 µC/(m 2 K) at T~ 270 K and (ii) -(∂P/∂T) ~ 1021 µC/(m 2 K) at T~ 320 K. In addition, a broad peak for electrocaloric temperature change is observed near 320 K with a relative cooling power of ~17 J/kg. These properties could be correlated to structural changes observed using X-ray diffraction at two different temperature ranges in the material. Analysis of high-energy X-ray scattering and specific heat capacity data revealed a transition from the cubic to tetragonal phase near T m ~ 320 K, whereas an additional increase in the tetragonality (c/a) of the polar phase is observed below T s ~ 270 K.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Enhanced Low-Temperature Thermoelectric Performance in (PbSe) 1+δ (VSe 2 ) 1 Heterostructures due to Highly Correlated Electrons in Charge Density Waves

We explore the effect of charge density wave (CDW) on the in-plane thermoelectric transport properties of (PbSe) 1+δ (VSe 2 ) 1 and (PbSe) 1+δ (VSe 2 ) 2 heterostructures. In (PbSe) 1+δ (VSe 2 ) 1 we observe an abrupt 86% increase in the Seebeck coefficient, 245% increase in the power factor, and a slight decrease in resistivity over the CDW transition. This behavior is not observed in (PbSe) 1+δ (VSe 2 ) 2 and is rather unusual compared to the general trend observed in other materials. The abrupt transition causes a deviation from the Mott relationship through correlated electron states. Raman spectra of the (PbSe) 1+δ (VSe 2 ) 1 material show the emergence of additional peaks below the CDW transition temperature associated with VSe 2 material. Temperature-dependent in-plane X-ray diffraction (XRD) spectra show a change in the in-plane thermal expansion of VSe 2 in (PbSe) 1+δ (VSe 2 ) 1 due to lattice distortion. Here, the increase in the power factor and decrease in the resistivity due to CDW suggest a potential mechanism for enhancing the thermoelectric performance at the low temperature region.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Simplified calculation of the area specific impedance for solid-state battery design

Simplified algebraic area specific impedance (ASI) correlations have been developed for solid-state composite battery electrodes made of a single ion conducting electrolyte, conductive additive, and intercalation active material. Two ASI expressions were developed, one for short times ( i.e. , pulsed power operation) and another for the pseudo steady state operation (i.e., sustained discharge for energy estimation). A full electrochemical model based on porous electrode theory was developed to examine the accuracy of the simplified ASI expressions. The simplified expressions agree favorably with full model results over a wide range of parameters (i.e. , electrode thicknesses, electrolyte conductivities, solid-state diffusion coefficients, specific surface areas, etc.) and conditions (i.e. , C-rates, states of charge, and pulse times). Under most conditions, the error between the full model and the correlations is well below 7 %. Higher errors were observed for the pseudo steady state expression at high/ low states of charge where the assumption of uniform reaction distributions loses validity. Here, the short time ASI has higher error at low states of charge due to the nonlinearity of the open circuit voltage equation, which is assumed linear in the formulation of the simplified algebraic expression.

25 ENERGY STORAGE↗
DOE OSTI2024

Improving the modeling of near-wall interphase heat transfer in porous media models of Pebble Bed Reactors

Here, this work aims to improve capabilities for modeling localized effects in porous media models of Pebble Bed Reactors. The wall-channeling effect is the primary local phenomenon of interest in a PBR, where the presence of the reflector wall disrupts the pebble packing, causing the pebbles near the wall to pack less efficiently and creating large void regions. Accurate modeling of the near-wall region is important as it will affect core bypass flow and temperature predictions. Porous media models are commonly used for design scoping and plant-level simulations of PBRs. Although these models have some capabilities to model the near-wall region, the correlations that are available in porous media codes are often inaccurate when a multi-region model is used to discretize the near-wall region. This work employs a high-to-low analysis to study the accuracy of available interphase heat transfer closures. NekRS, a spectral element computational fluid dynamics code, is used to perform Large Eddy Simulations. These LES simulation results are compared to porous media model results from the Pronghorn porous media code. The friction term of the KTA drag closure is first improved, reducing the error in the prediction of the near-wall velocity from over 50% to less than 5%. This is combined with improvements to the form term from previous works to produce a drag closure that is capable of accurately modeling the wall-channeling effect across a variety of flow conditions. The Nusselt number predictions of several heat transfer correlations are compared to the high-fidelity results where it is found that the KTA heat transfer correlation is capable of accurately predicting the local Nusselt numbers that were determined in the high-fidelity simulation. Comparison of the radial solid temperature profiles, however, reveal discrepancies between NekRS and Pronghorn. It is discovered that the implementation of the interphase heat transfer coefficient that exists in many current porous media codes is not valid when local porosities are modeled. Instead, it is suggested that the interphase heat transfer coefficient should be dependent on the local porosity, the Nusselt number, and the local solid surface-to-volume ratio. Implementation of this change produces improvement in the agreement between the results obtained by NekRS and Pronghorn while using the KTA heat transfer correlation.

interphase heat transfer↗
DOE OSTI2023

Particle Receiver Models for Systems Analysis

Particle receivers are gaining importance in the field of Concentrating Solar Power (CSP) due to the high temperature that particles can achieve without degradation. Several researchers are studying the potential of this technology by means of system analyses, which need simple and light models of the system. This study presents two simple models for the particle receiver. The simplest model is a correlation obtained by fitting the results calculated with a more complex receiver model simulated in CFD. The other model is a 1D model, which is benchmarked against the same CFD results. Although both models achieve high coefficient of determination, R2, when compared to CFD results, the 1D model seems to provide more accurate results (especially during sunsets and sunrises). Both models are integrated into a tecno-economic model developed in previous work. The LCOE obtained with the 1D model is between 7% and 10% greater than the one obtained with the correlation.

González-Portillo, Luis F. (ORCID:0000000348643825↗
DOE OSTI2023

Gallium isotopic constraints for the origin of the Earth-Moon system

Comparatively heavy isotopic compositions of moderately volatile elements (MVE) in lunar rocks have been advocated to reflect the loss of light isotopes during devolatilization processes from the Moon. In this study we present new gallium (Ga) isotope data for lunar highland rocks, with a focus on the Ferroan Anorthosite Suite (FAS). These are commonly thought to be direct crystallization products from the late lunar magma ocean (LMO) and should contain the majority of the Ga inventory of the Moon. As such, FAS rocks are crucial for identifying the processes that drove Ga isotope fractionation as well as for inferring the Ga isotopic composition of the bulk Moon. Our data reveal that FAS samples have a range in δ 71 Ga from -0.27 to 0.22‰ and are generally isotopically light in Ga compared to other lunar lithologies, but straddle values typical of terrestrial rocks. Although Ga is defined as an MVE, these Ga isotope variations do not correspond with concentrations of more volatile elements, indicating that Ga isotope variations in the FAS are not primarily controlled by devolatilization processes. Instead, the Ga isotopic compositions of bulk FAS rocks broadly correlate with the composition of plagioclase, with the calcium content of plagioclase decreasing as Ga becomes isotopically heavier. This suggests that fractionation of Ga isotopes in FAS rocks was caused by the preferential incorporation of isotopically light Ga into plagioclase during the later solidification stages of the LMO. The progressive crystallization and extraction of plagioclase forces the residual melt towards increasingly heavier Ga isotope ratios, corroborating similar conclusions derived from correlations between δ 71 and Eu* in the mare basalt suite. Using Ga isotope partitioning calculations, we demonstrate that an isotope fractionation coefficient between plagioclase and coexisting melt of -0.3 to -0.4‰ could explain the observed range of δ 71 values in FAS, mare basalt suite rocks, and KREEP. These calculations allow for a first order estimate of the Ga isotopic composition of the bulk silicate Moon prior to plagioclase fractionation and suggest it was close to the composition of the bulk silicate Earth. This would imply that the Moon did not lose a substantial fraction of its Ga inventory during accretion, consistent with new constraints from Rb-Sr isotope systematics that indicate the Moon's volatile deficit was primarily inherited from Theia. In conjunction with the overlap in non-mass-dependent isotope ratios, these collective observations could be reconciled if Theia and the proto-Earth formed in similar regions of the inner Solar System that were already volatile-depleted.

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

Pore connectivity influences mass transport in natural rocks: Pore structure, gas diffusion and batch sorption studies

For this work, six rocks (one granodiorite, one limestone, two chalks, one mudstone, and one dolostone) with different extents of heterogeneity at six different particle sizes (from 75 to 8000 μm) were studied to describe the effects of pore connectivity on mass transport. The methods applied were (i) porosity measurement of granular rocks, (ii) analyses of gas-phase diffusive transport in a bed of packed particles, along with a solid quartz method at these six particle sizes being developed to identify the contribution of intraparticle diffusion, and (iii) batch sorption tests of multiple ions (anions and cations) with subsequent analyses of inductively coupled plasma-mass spectrometry. Granular porosity measurement results reveal that with decreasing particle sizes, the effective porosities for the “heterogenous” group of rocks (Grimsel granodiorite and Edwards limestone) increase, whereas the porosities of another “homogeneous” group (two Israel chalk samples, Japan mudstone, and Wyoming dolostone) remain constant. Gas diffusion results show that the intraparticle gas diffusion coefficient among these two sample groups, varying in the magnitude of 10 -8 to 10 -6 m 2 /s, are not directly correlated to the porosity differences. Moreover, the batch sorption work displays a different affinity of rocks for various tracers. For Grimsel granodiorite, Japan mudstone, and Wyoming dolostone, the adsorption capacity of Sm 3+ and Eu 3+ increases as the particle size decreases. In general, this integrated research of grain size distribution, granular rock porosity, intraparticle diffusivity, and ionic sorption capacity gives insights into the pore connectivity effect on both physical and chemical transport behaviors for different lithologies and/or different particle sizes.

58 GEOSCIENCES↗
DOE OSTI2021

Toward a First-Principles Framework for Predicting Collective Properties of Electrolytes

Conspectus Given the universal importance of electrolyte solutions, it is natural to expect that we have a nearly complete understanding of the fundamental properties of these solutions (e.g., the chemical potential) and that we can therefore explain, predict, and control the phenomena occurring in them. In fact, reality falls short of these expectations. But, recent advances in the simulation and modeling of electrolyte solutions indicate that it should soon be possible to make progress toward these goals. In this Account, we will discuss the use of first-principles interaction potentials based in quantum mechanics (QM) to enhance our understanding of electrolyte solutions. Specifically, we will focus on the use of quantum density functional theory (DFT) combined with molecular dynamics simulation (DFT-MD) as the foundation for our approach. The overarching concept is to understand and accurately reproduce the balance between local or short-ranged (SR) structural details and long-range (LR) correlations, allowing the prediction of the thermodynamics of both single ions in solution as well as the collective interactions characterized by activity/osmotic coefficients. In doing so, relevant collective motions and driving forces characterized by chemical potentials can be determined. Here, in this Account, we will make the case that understanding electrolyte solutions requires a faithful QM representation of the SR nature of the ion–ion, ion–water, and water–water interactions. However, the number of molecules that is required for collective behavior makes the direct application of high-level QM methods that contain the best SR physics untenable, making methods that balance accuracy and efficiency a practical goal. Alternatives such as continuum solvent models (CSMs) and empirically based classical molecular dynamics have been extensively employed to resolve this problem but without yet overcoming the fundamental issue of SR accuracy. We will demonstrate that accurately describing the SR interaction is imperative for predicting both intrinsic properties, namely, at infinite dilution, and collective properties of electrolyte solutions. DFT has played an important role in our understanding of condensed phase systems, e.g., bulk liquid water, the air–water interface, ions in bulk, and at the air–water interface. This approach holds huge promise to provide benchmark calculations of electrolyte solution properties that will allow for the development and improvement of more efficient methods, as well as an enhanced understanding of fundamental phenomena. However, the standard protocol using the generalized gradient approximation with van der Waals (vdW) correction requires improvement in order to achieve a high level of quantitative accuracy. Simply simulating with higher level DFT functionals may not be the best route considering the significant computational cost. Alternative methods of incorporating information from higher levels of QM should be explored; e.g., using force matching techniques on small clusters, where high level benchmark calculations are possible, to develop ideal correction terms to the DFT functional is a promising possibility. We argue that DFT with statistical mechanics is becoming an increasingly useful framework enabling the prediction of collective electrolyte properties.

Duignan, Timothy T.↗
DOE OSTI2023

Tuning the Radius Ratio to Enhance Thermoelectric Properties in the Zintl Compounds AM 2 Sb 2 (A = Ba, Sr; M = Zn, Cd)

Five novel Zintl phase solid solutions in the Ba 1–x Sr x Zn 2–y Cd y Sb 2 (0 ≤ x ≤ 0.13(1); 0 ≤ y ≤ 0.32(2)) system were successfully synthesized by the molten Pb metal-flux method, and the powder X-ray diffraction and single-crystal X-ray diffraction analyses proved that all five title compounds adopted the BaCu 2 S 2 -type phase having the orthorhombic Pnma space group (Z = 4, Pearson code oP20) with five crystallographically independent atomic sites. The previously studied BaCu 2 S 2 -type antimonides demonstrated a limited tolerance for doping in contrast to the CaAl 2 Si 2 -type antimonides. To understand the relatively narrower phase width and limited dopability of the title BaCu 2 S 2 -type phase than the CaAl 2 Si 2 -type phase in the overall Ba 1–x Sr x Zn 2–y Cd y Sb 2 system, the radius ratio of cations and anionic elements r + /r – for two structure types were thoroughly investigated. For the first time, the r + /r – ratio was identified as a critical factor for the phase selectivity: (1) r + /r – > 1 favored the BaCu 2 S 2 -type phase, and (2) r + /r – < 1 favored the CaAl 2 Si 2 -type phase. Further, we also revealed the structural transformation mechanism from the more widely observed CaAl 2 Si 2 -type phase to the title BaCu 2 S 2 -type phase as the relatively larger cationic elements were introduced to the system. A series of DFT calculations using the three hypothetical models indicated that a resonance peak near EF in the density of states curves was descended from the relatively flat band structure at several special symmetry points rationalizing the enhanced Seebeck coefficients of Ba 0.94(1) Sr 0.06 Zn 1.86(3) Cd 0.14 Sb 2 and Ba 0.96(1) Sr 0.04 Zn 1.68(2) Cd 0.32 Sb 2 . Electron localization function analysis rationalized the correlation between the polarity change of anionic Zn/Cd–Sb bonds and the charge carrier mobility on the anionic frameworks. Temperature-dependent thermoelectric properties were studied for the four title compounds, and the results proved that the Sr and Cd doping in the title Ba 1–x Sr x Zn 2–y Cd y Sb 2 system successfully enhanced the ZT values through the increased Seebeck coefficients and the reduced total thermal conductivities.

36 MATERIALS SCIENCE↗