Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “continuous improvement”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

High yield, large-scale synthesis of calcium-based microporous metal-organic framework and examination of the long-term stability for xenon adsorption applications

Here, scale-up synthesis of calcium-based microporous metal-organic framework (SBMOF-1) up to ~400 g in a batch with a yield of >90 % was achieved by a solvothermal reaction of sulfonyldibenzoic acid (SDB) with an excess of calcium chloride. Here, we observed that recrystallization of unreacted SDB at a solvothermal condition caused a moderate reaction yield (40–50 %) at the reference condition of $C_{CaCl_2}$/$C_{SDB}$=1 and $C_{SDB}$/$C_{SDB(ref)}$=1. Simply adding more reagents to the reactor did not increase the mass of product formed per unit volume due to a more pronounced loss of the yield at those conditions. By simultaneously changing the molar ratio of CaCl 2 to SDB, $C_{CaCl_2}$/$C_{SDB}$, and the molar concentration of the SDB reagent, $C_{SDB}$/$C_{SDB(ref)}$, we explored %yield of the reaction. Interestingly, a linear improvement in the yield was observed from 21% (at $C_{CaCl_2}$/$C_{SDB}$=0.5) up to 78% (at $C_{CaCl_2}$/$C_{SDB}$=6) at a fixed ratio of $C_{SDB}$/$C_{SDB(ref)}$=2 and the yield leveled off after further addition. Unlike those at $C_{CaCl_2}$/$C_{SDB}$=1, the yields at a high CaCl 2 excess continued to improve with increasing the $C_{SDB}$/$C_{SDB(ref)}$. When a large pressure vessel (2500 mL EtOH, $C_{CaCl_2}$/$C_{SDB}$=6, $C_{SDB}$/$C_{SDB(ref)}$=8) was used, about 415 g of SBMOF-1 with a yield of 92.3% was produced, indicating 16 × the space yield improvement. The ability to synthesize SBMOF-1 on a large scale allowed us to examine the long-term stability of SBMOF-1 for almost 200 days in the presence of varying levels of relative humidity.

36 MATERIALS SCIENCE↗
DOE OSTI2023

The Dark Energy Survey Year 3 and eBOSS: constraining galaxy intrinsic alignments across luminosity and colour space

ABSTRACT We present direct constraints on galaxy intrinsic alignments (IAs) using the Dark Energy Survey Year 3 (DES Y3), the Extended Baryon Oscillation Spectroscopic Survey (eBOSS), and its precursor, the Baryon Oscillation Spectroscopic Survey (BOSS). Our measurements incorporate photometric red sequence (redMaGiC) galaxies from DES with median redshift z ∼ 0.2–1.0, luminous red galaxies from eBOSS at z ∼ 0.8, and also an SDSS-III BOSS CMASS sample at z ∼ 0.5. We measure two-point IA correlations, which we fit using a model that includes lensing, magnification, and photometric redshift error. Fitting on scales 6 Mpc h−1 < rp < 70 Mpc h−1, we make a detection of IAs in each sample, at 5σ–22σ (assuming a simple one-parameter model for IAs). Using these red samples, we measure the IA–luminosity relation. Our results are statistically consistent with previous results, but offer a significant improvement in constraining power, particularly at low luminosity. With this improved precision, we see detectable dependence on colour between broadly defined red samples. It is likely that a more sophisticated approach than a binary red/blue split, which jointly considers colour and luminosity dependence in the IA signal, will be needed in future. We also compare the various signal components at the best-fitting point in parameter space for each sample, and find that magnification and lensing contribute $\sim 2\!-\!18~{{\ \rm per\ cent}}$ of the total signal. As precision continues to improve, it will certainly be necessary to account for these effects in future direct IA measurements. Finally, we make equivalent measurements on a sample of emission-line galaxies from eBOSS at z ∼ 0.8. We constrain the non-linear alignment amplitude to be $A_1=0.07^{+0.32}_{-0.42}$ (|A1| < 0.78 at 95 per cent CL).

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Annual Stormwater Report for the Y-12 National Security Complex Oak Ridge, Tennessee

This is the eighth annual stormwater report prepared in accordance with the modified National Pollutant Discharge Elimination System permit issued to the Y-12 National Security Complex (Y-12) on April 1, 2014. This modified permit removed the authorization to augment flow in East Fork Poplar Creek (EFPC) and thus removed the requirement to sample this effluent in conjunction with other major outfalls on EFPC. The general condition of the site continues to improve. The repair of Outfall 014 is complete and sampling was performed in 2020 from the new Outfall 014 structure. This repair marks a significant improvement in the storm drain network infrastructure. There are several construction and demolition activities currently in progress at Y-12. These include the construction of the Uranium Processing Facility, the Mercury Treatment Facility on EFPC, the site-wide power grid upgrade, and the demolition of the Biology Complex. Current stormwater protections have proven effective in 2020; however, moving forward, these projects will present challenges for the protection of stormwater. Therefore, oversight of these activities, and all other construction and demolition projects, will continue to be a priority in 2021. Y-12 will continue to implement appropriate best management practices and reduce outside material storage areas where possible.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Digital image correlation and infrared thermography data for seven unique geometries of 304L stainless steel

Material Testing 2.0 (MT2.0) is a paradigm that advocates for the use of rich, full-field data, such as from digital image correlation and infrared thermography, for material identification. By employing heterogeneous, multi-axial data in conjunction with sophisticated inverse calibration techniques such as finite element model updating and the virtual fields method, MT2.0 aims to reduce the number of specimens needed for material identification and to increase confidence in the calibration results. To support continued development, improvement, and validation of such inverse methods—specifically for rate-dependent, temperature-dependent, and anisotropic metal plasticity models—we provide here a thorough experimental data set for 304L stainless steel sheet metal. The data set includes full-field displacement, strain, and temperature data for seven unique specimen geometries tested at different strain rates and in different material orientations. Commensurate extensometer strain data from tensile dog bones is provided as well for comparison. We believe this complete data set will be a valuable contribution to the experimental and computational mechanics communities, supporting continued advances in material identification methods.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Analysis of Lobe Power Calculator and Indication System with Physics and Cycle Based Models

The Advanced Test Reactor (ATR) at INL measures reactor power through various methods, two of them being thermal and Nitrogen-16 (N-16) activity. Water power calculator (WPC) is a thermal power system that measures flow and temperature to determine thermal quadrant powers. The N-16 system utilizes a beta detector that outputs Nitrogen activation levels to calculate lobe power through an algorithm called Lobe Power Calculation and Indication system (LPCIS). The LPCIS utilizes the N-16 system and the WPC system to determine reactor core power levels. The WPC provides accurate calculations of quadrant and total reactor thermal power. With the use of WPC measurements, thermal-to-N-16 (T2N) power ratios are produced to determine if the two indication systems agree on core power. Relative magnitude equations are used to utilize N-16 coefficients and multipliers to improve the indications of the LPCIS. These correct indications are crucial for maintaining safety limits because operators rely on this information for decision making. Currently, the LPCIS system uses linear equations and matrices to calculate lobe power through multipliers and coefficients. Advancements in technology and system upgrades have increased the accuracy of power readings by making the system more dynamic. The new proposed coefficient and multiplier method implements a cycle specific and physics-based model intended for changing coefficients during operation. This calibration experiment focused on power splitting, outer shim and neck shim, as well as fuel burning into the reactor digital acquisition system (RDAS) weighting factors. Results demonstrated that the physics learning method yields a smaller error margin inside of the desired power range for the data set 166-A. Continuing to improve the physics-based model will help improve the power accuracy of the LPCIS system.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2021

Implications of an improved water equation of state for water-rich planets

Water (H 2 O), in all forms, is an important constituent in planetary bodies, controlling habitability and influencing geological activity. Under conditions found in the interior of many planets, as the pressure increases, the H-bonds in water gradually weaken and are replaced by ionic bonds. Recent experimental measurements of the water equation of state (EOS) showed both a new phase of H-bonded water ice, ice-VIIt, and a relatively low transition pressure just above 30 GPa to ionic bonded ice-X, which has a bulk modulus 2.5 times larger. The higher bulk modulus of ice-X produces larger planets for a given mass, thereby either reducing the atmospheric contribution to the volume of many exoplanets or limiting their water content. We investigate the impact of the new EOS measurements on the planetary mass–radius relation and interior structure for water-rich planets. We find that the change in the planet mass–radius relation caused by the systematic differences between previous and new experimental EOS measurements is comparable to the observational uncertainties in some planet sizes – an issue that will become more important as observations continue to improve.

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

Analysis of Lobe Power Calculator and Indication System with Physics and Cycle Based Models

The Advanced Test Reactor (ATR) at INL measures reactor power through two methods, thermal and Nitrogen-16 (N-16) activity. Water power calculator (WPC) is a thermal power system that measures flow and temperature to determine the thermal quadrant powers. The N-16 system utilizes a beta chamber detector that outputs reactivity levels to calculate lobe power through an algorithm called lobe power calculation and indication system (LPCIS). The LPCIS utilizes the N-16 system and the WPC system to determine reactor core power levels. The WPC provides accurate calculations of quadrant and total reactor thermal power. With the use of WPC measurements, thermal-to-N-16 (T2N) power ratios are produced to determine if the two indication systems agree on core power. Relative magnitude equations are used to utilize N-16 coefficients and multipliers to improve the indications of the LPCIS. This is crucial for maintaining safety limits. Currently, the LPCIS system uses linear equations and matrices to calculate lobe power through multipliers and coefficients. Advancements in technology and system upgrades have assisted system engineers at INL with the objective to reach a more dynamic system. In return the system demonstrates an increase in the accuracy of power reading while maintaining safety margins. The new proposed coefficient and multiple method implements a cycle specific, and a physics-based model intended for changing coefficients during operation. This calibration experiment focuses on power splitting, outer shim and neck shim, as well as fuel burning into the RDAS weighting factors. Results show that the Physics learning method yields a smaller error margin inside of the desired power range for the data set 166-A. This proved true for both constrained and unconstrained testing. This is most likely due to the physics data fitting approach that resulted in favorable coefficients and multipliers for 166-A. Continuing to improve the physics-based model will help improve the power accuracy of the LPCIS system.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2023

Field‐scale analysis of miscanthus production indicates climate change may increase the opportunity for water quality improvement in a key Iowa watershed

Abstract The Raccoon River Basin is the primary source for drinking water in Iowa's largest city and plays a major role in the Mississippi River Basin's high nutrient exports. Future climate change may have major impacts on the biological, physiological, and agronomic processes imposing a threat to ecosystem services. Efforts to reduce nitrogen (N) loads within this basin have included local litigation and the implementation of the Iowa Nutrient Reduction Strategy, which suggest incorporating bioenergy crops (i.e., miscanthus) within the current corn–soybean landscape to reach a 41% reduction in nitrate loads. This study focuses on simulating N export for historical and future land use scenarios by using an agroecosystem model (Agro‐IBIS) and a hydrology model (THMB) at the 500‐m resolution, similar to the scale of agricultural fields. Model simulations are driven by CMIP5 climate data for historical, mid‐century, and late‐century under the RCP 4.5 and 8.5 warming projections. Using recent crop profit analyses for the state of Iowa, profitability maps were generated and nitrogen leaching thresholds were used to determine where miscanthus should replace corn–soybean area to maximize reductions in N pollution. Our results show that miscanthus inclusion on low profit and high N leaching areas can result in a 4% reduction of N loss under current climate conditions and may reduce N loss by 21%–26% under future climate conditions, implying that water quality has the potential continue to improve under future climate conditions when strategically implemented conservation practices are included in future farm management plans.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Vertical z-axis discontinuous carbon fibers for improved lightning strike performance of continuous fiber-reinforced polymer composites

Effective lightning strike protection for critical aerospace and wind applications requires high electrical conductivity to dissipate current efficiently. However, polymer matrix composites face a challenge due to their inherently insulating nature. While conventional carbon fiber-reinforced composites (CFRP) exhibit electrical conductivity in the planar direction, achieving through-thickness conductivity remains an ongoing challenge. In this work, we have undertaken the fabrication of CFRP interleaved with vertically oriented carbon fibers (Z-fiber) to impart higher electrical conductivity along the thickness direction. Two Z-fiber composite variations are prepared: Z-1 with a single layer of Z-fiber and Z-5 with five interleaved layers and compared with no Z-fiber layer (Z-0) composite. The composite panels were subjected to lab-scale lightning strike tests with a current magnitude of 100 kA. To emulate real-world service conditions, an aerospace-grade paint coating was applied to the composite laminates. Comparative analysis shows Z-1 reduces damage diameter to ∼22 mm compared to Z-0 (∼26 mm), while Z-5 exhibits the least damage (∼16.7 mm), confirmed by optical microscopy. Z-5 demonstrates nine times higher through-thickness electrical conductivity than Z-0, reducing electrical anisotropy substantially. Thermal-electric finite element damage modeling predicts surface damage within 6% of experimental values for both Z-0 and Z-5 composites. Flexural tests post-lightning reveal Z-5 retains 66% flexural strength and 86% modulus, significantly better than Z-0, which retains less than 40% for both properties. This study highlights the efficacy of Z-fiber composites in lightning strike protection, offering improved through-thickness conductivity and mechanical property retention.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Evaluation of H-Canyon Ventilation Exhaust Tunnel Inspection at the Savannah River Site - 20052

The Department of Energy (DOE) H-Canyon Facility at the Savannah River Site (SRS) is operated by Savannah River Nuclear Solutions (SRNS). The Facility Structural Integrity Program [1] provides a process to inspect, evaluate, and document the conditions of passive safety related structural systems and components, their degradation mechanisms and their impact on the safety envelope. The Structural Integrity Program performs periodic visual inspections to confirm the systems and components can perform their safety function and recommend corrective actions before the function is compromised. The H-Canyon Exhaust Tunnel is part of the H-Canyon Ventilation System and is periodically inspected under the Structural Integrity Program. The Tunnel performs a Safety Class passive design feature function that is available 100% of the time. The underground reinforced concrete Tunnel directs Canyon process air from the Canyon to the Sand Filter System. The radiological and chemical airborne activity coupled with the physical configuration of the tunnel precludes a manned entry into the tunnel to perform the inspection. Unlike active systems that can be subjected to physical testing to confirm capability, qualification calculations, including modeling, is performed to assess passive design feature capabilities relative to operational, accident and natural phenomena hazards events. Input for these calculations come from knowledge of field conditions, design drawings and conservative assumptions. The knowledge of field conditions comes from concrete and soil testing activities as well as inspection data. The use of cameras on crawlers to obtain inspection data and how that data is used in qualification calculations is the focus of this discussion. Remote Tunnel inspections have been performed using cameras on a stick, referred to as pole cameras, prior to Calendar Year 2003 (CY2003). Since CY2003, six camera equipped crawlers have been used. Crawlers have evolved and improved by incorporating lessons learned from the previous inspections. The latest 2019 crawler successfully traveled the entire length of the tunnel between the Canyon and Sand Filter System and captured images that were previously unseen, thus providing key visual data for the Structural Integrity Program report [2]. The inspection data serves as 1) input to qualification calculations, 2) validates current conditions remain bounded by qualification calculations, and 3) documents a record of change over time. This data is vital for the effective evaluation of the tunnel structure and aids with predicting the life cycle service life safety envelope to protect the public, environment, and facility worker. The goal of the 2019 deployment was to use a more robust vehicle to ensure travel over obstacles, perform the inspection with higher resolution cameras, and have the ability to elevate the cameras to view surfaces previously obscured by abandon-in-place ducts in the tunnel. The inspection information will serve to strengthen the input basis for qualification calculations, versus using excessive conservative calculation assumptions. This 2019 inspection was compared with previous crawler visual evidence to identify change over time to project remaining service life and confirm field conditions remain bounded by qualification calculations. The results of the 2019 Exhaust Tunnel inspection continue to improve and provide a firm basis to make sound engineering judgements regarding tunnel capability and approaches in projecting remaining service life with an aging system and component. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2022

Advances in SnO 2 for Efficient and Stable n-i-p Perovskite Solar Cells

Perovskite solar cells (PSCs) based on the regular n-i-p device architecture have reached above 25% certified efficiency with continuously reported improvements in recent years. A key common factor for these recent breakthroughs is the development of SnO 2 as an effective electron transport layer in these devices. In this review, we discuss the key advances in SnO 2 development, including various deposition approaches and surface treatment strategies, to enhance the bulk and interface properties of SnO 2 for highly efficient and stable n-i-p PSCs. We also discuss the general materials chemistry associated with SnO 2 along with the corresponding materials challenges and improvement strategies, focusing on defects, intrinsic properties, and impact on device characteristics. Finally, we highlight some SnO2 implementations related to scalable processes and flexible devices, and we also provide our perspective on the future development of efficient and stable large-scale perovskite solar modules.

14 SOLAR ENERGY↗
DOE OSTI2026

MFC 5.0: An exascale many-physics flow solver

Many problems of interest in engineering, medicine, and the fundamental sciences rely on high-fidelity flow simulation, making performant computational fluid dynamics solvers a mainstay of the open-source software community. Previous work MFC 3.0 was made a published, documented, and open-source solver via Bryngelson et al. Comp. Phys. Comm. (2021) with numerous physical features, numerical methods, and scalable infrastructure. MFC 5.0 is a significant update to MFC 3.0, featuring a broad set of well-established and novel physical models and numerical methods, as well as the introduction of GPU and APU (or superchip) acceleration. Here, we exhibit state-of-the-art performance and ideal scaling on the first two exascale supercomputers, OLCF Frontier and LLNL El Capitan. Combined with MFC’s single-accelerator performance, MFC achieves exascale computation in practice, and achieved the largest-to-date public CFD simulation at 200 trillion grid points as a 2025 ACM Gordon Bell Prize finalist. New physical features include the immersed boundary method, N-fluid phase change, Euler–Euler and Euler–Lagrange sub-grid bubble models, fluid-structure interaction, hypo- and hyper-elastic materials, chemically reacting flow, two-material surface tension, magnetohydrodynamics (MHD), and more. Numerical techniques now represent the current state-of-the-art, including general relaxation characteristic boundary conditions, WENO variants, Strang splitting for stiff sub-grid flow features, and low Mach number treatments. Weak scaling to tens of thousands of GPUs on OLCF Summit and Frontier and LLNL El Capitan achieves efficiencies within 5% of ideal to over 90% of their respective system sizes. Strong scaling results for a 16-times increase in device count show parallel efficiencies over 90% on OLCF Frontier. MFC’s software stack has undergone further improvements, including continuous integration, which ensures code resilience and correctness through over 300 regression tests; metaprogramming, which reduces code length while maintaining performance portability; and code generation for computing chemical reactions

Computational fluid dynamics↗
DOE OSTI2023

Neural network-based control of an ultrafast laser

With the recent advances in machine learning (ML) and data science (DS), the control, modeling, and analysis of these complex systems continues to improve. In this work, we report on the optimization of the intensity of a femtosecond laser using feedforward neural networks (FFNN) that model the input–output relationships of the data. The input parameters of the system were optimized to achieve the required performance of the femtosecond laser. We propose a neural network-based control system to model the relationship between the spectral amplitude and phase of the input laser pulse at the amplifier input and the shape of the output pulse. Low-jitter laser parameter inputs and the resulting laser pulse duration were modeled, and the resulting correlation between the input and output data was used to optimize the laser pulse. Here, we demonstrate improved processing and laser control performance.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2020

Rheology and stability of nanoparticle-stabilized CO 2 foam under reservoir conditions

This work investigates the rheology and stability of nanoparticle-stabilized CO 2 foams under reservoir conditions (high temperature and high pressure) for fracturing applications. The effects of different parameters on foam apparent viscosity and foam stability were experimentally investigated including the effects of nanoparticle concentration, salinity, foam quality (Γ), shear rate and temperature. The power law model was applied to calculate foam apparent viscosity owing to its pseudo-plastic behavior, and the changes of foam heights over time were used to evaluate foam stability. Results showed that the CO 2 foam apparent viscosity featured a mountain-shaped curve versus Γ, with the peak apparent viscosity obtained at 70% foam quality. The increase of salinity (up to 11%) in the continuous phase improved both foam stability and rheology. Higher nanoparticle concentration could contribute to better foam stability, but there was a threshold concentration, above which the foam apparent viscosity remained stabilized. The CO 2 foam stabilized by nanoparticles displayed a shear-thickening behavior as the experiment flow rate increased from 6 mL/min to 18 mL/min. Further studies showed that as the total flow rate increased, the CO 2 foam became finer-textured with better stability. Elevated temperatures could undermine foam apparent viscosity and long-term stability. The results of this study could provide guidelines as to the design of foam fracturing systems for field applications.

03 NATURAL GAS↗
DOE OSTI2025

Accelerating room air conditioner efficiency in India: Grid, economic, and policy implications through 2035

India is poised for a rapid surge in space cooling demand, driven by rising incomes, urbanization, and intensifying heat. Between 2025 and 2035, the country is expected to add 130–150 million new room air conditioners (ACs). If Minimum Energy Performance Standards (MEPS) continue to improve at the historical rate of 2–3 % annually, room ACs alone could contribute over 180 GW to peak electricity demand by 2035-nearly 30 % of the projected national total. This study evaluates the impact of an accelerated MEPS trajectory, proposing to raise the 1-star threshold to ISEER 5.0 by 2027, ISEER 6.3 by 2030, and ISEER 7.4 by 2033. Drawing on engineering cost analysis, stock turnover modeling, and retail pricing data, we find that this pathway could reduce peak demand by over 60 GW, save 118 TWh of electricity annually, avoid 49 MtCO₂ of electricity-related emissions per year, avert ₹7.5 trillion (∼US$85 billion) in power system investments, and yield ₹0.7–2.3 trillion (∼US$8–26 billion) in net consumer savings by 2035. Contrary to affordability concerns, empirical trends show that higher efficiency does not increase AC prices. These results highlight the value of ambitious MEPS as a cost-effective strategy for improving grid reliability, reducing emissions, and advancing consumer welfare in emerging economies.

Abhyankar, Nikit↗
DOE OSTI2026

Technology Impact and Resource Assessment of Existing and Planned U.S. Biofuel Production: Life Cycle Water Consumption, Water Stress, Land Use, and Criteria Air Pollutants

Biofuels have the potential to strengthen the U.S. energy supply, enhance energy security, and promote economic development. As the United States continues to expand biofuel production, quantifying resource requirements and location-specific constraints is crucial for planning, siting, and technology development to support long-term viability. Accordingly, this work assesses the life cycle resource consumption (water consumption, and land use), water stress and criteria air pollutants associated with expanded U.S. biofuel production over the 2020–2035 period, based on producers’ plans. We perform a bottom-up technology impact assessment and resource assessment by integrating facility-level production statistics with Argonne’s Research and Development (R&D) GREET model and county-level water-stress characterization factors from Available Water Remaining for the United States (AWARE-US) model. Results suggest that by 2035, biofuels could meaningfully contribute to U.S. energy demand, driven primarily by first-generation and waste-based feedstocks with plans for substantial capacity expansion, although cellulosic and e-fuel technologies remain limited. However, growth must be managed to minimize impacts on water resources and land use. These impacts vary by fuel type and facility location, specifically, projected expansion increases water consumption and can elevate water-stress impacts in certain regions like Nebraska, Kansas, Colorado, Idaho, California and North Texas. Direct land use also increases overall, particularly for first-generation feedstocks such as corn and soybeans. These findings underscore the need for continued technological improvements and innovative strategies to manage resource demands as the industry scales and to support complementary deployment within the evolving U.S. energy system.

09 BIOMASS FUELS↗
DOE OSTI2024

Role of Dynamic Polarization Interactions in the Electrical Double Layer at Calcite (104) Interfaces with Aqueous Solutions

Reactions at mineral interfaces with aqueous solutions control many geochemical and biogeochemical processes in the Earth’s critical zone. At the molecular level, insights into important properties such as the structure of the electrical double layer (EDL) at specific mineral interfaces continue to improve due to the increasing fidelity of laboratory instrumentation and computational approaches. However, molecular simulation approaches suffer from limited reach into relevant scales of time, length, and system complexity. To span this gap, a novel hybrid approach that couples first principle plane-wave density functional theory (DFT) with classical DFT (cDFT) is demonstrated and applied to calcite (104) interfaces with various electrolytes. In this approach, a region of interest described using DFT interacts with the surrounding medium described using cDFT to arrive at a self-consistent ground state. Benchmarking against experimental observations and entirely first principle DFT simulations demonstrates that this hybrid model efficiently encompasses the key short-range and collective interactions in the EDL. Simulations of calcite (104)/solution interfaces reveal the key static and dynamic polarization interactions that give rise to structuring of ions and water. Ion hydration interactions have the strongest effect on the depth of the first minimum in the density distribution of counterions at the surface, and the position and width of the first density peak is largely determined by the strength of ion-correlation forces. Finer details of ion distributions are controlled by mutual polarization of the calcite surface and interfacial electrolyte. Finally, this new ability to efficiently and rigorously predict EDL structure at mineral surfaces in contact with complex solutions paves the way to accurately modeling sorption, nucleation, dissolution, and growth in realistic systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗