Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Interface properties”

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

Preparation of Ultrathin Gold Films with Subatomic Surface Roughness

Nanoscale device fabrication requires control over film growth at the atomic scale. Growth conditions must be tuned in consideration of interface parameters like chemical bonding, surface free energy, and lattice matching. In metals, electronic properties may also be utilized for control of physical parameters. Quantum size effects can induce metals to spontaneously form specific shapes and sizes according to their electronic structure. Unfortunately, such electronic growth is generally known only for a few systems and is typically only stable under cryogenic conditions. In this work, we explore a recently discovered class of electronic growth systems in which metal films are grown upon the relatively inert surfaces of van der Waals crystals. In this class of materials, the electronic growth is highly stable at room temperature and actually requires higher temperature annealing to achieve proper equilibrium. We work with the Au/MoS 2 system, which shows excellent stability and can readily form discrete and atomically flat nanostructures. Here, we show how the electronic growth modes facilitate the formation of atomically flat films with nanometer scale thickness. Here, the surface roughness of these films was found to be less than a single atom over several square microns, creating nearly perfect surfaces for studies of self-assembled monolayers or other applications.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Signatures of Coherent Phonon Transport in Ultralow Thermal Conductivity Two-Dimensional Ruddlesden–Popper Phase Perovskites

An emerging class of methylammonium lead iodide (MAPbI 3 )-based Ruddlesden-Popper (RP) phase perovskites, BA 2 MA n-1 Pb n I 3n+1 (n = 1-7), exhibit enhanced stability to environmental conditions relative to MAPbI 3 , yet still degrade at elevated temperatures. Here, we experimentally determine the thermal conductivities of these layered RP phases for n = 1-6, where n defines the number of repeated perovskite octahedra per layer. We measure thermal conductivities of 0.37 ± 0.13/0.12, 0.17 ± 0.08/0.07, 0.21 ± 0.05/0.04, and 0.19 ± 0.04/0.03 W/m∙K in thin films of n = 1-4 and 0.08 ± 0.06/0.04, 0.06 ± 0.04/0.03, 0.06 ± 0.03/0.03, and 0.08 ± 0.07/0.04 W/m∙K in single crystals of n = 3-6. With the exception of n = 1, these thermal conductivities are lower than the range of 0.34-0.50 W/m∙K reported for single-crystal MAPbI 3 . Reduced-order lattice dynamics modeling suggests that the initially decreasing trend of thermal conductivity in similarly oriented perovskites with increasing n may result from the transport properties of coherent phonons, emergent from the superstructure, that do not scatter at the interfaces of organic butylammonium chains and perovskite octahedra. Reduced group velocity of coherent phonons in n = 3-6, a consequence of band flattening in the phonon dispersion, is primarily responsible for their ultralow thermal conductivities. Similar effects on thermal conductivity have been experimentally demonstrated in deposited superlattices, but never in naturally defined materials such as RP phases. GIWAXS measurements reveal that higher n RP phase thin films are less orientationally controlled and therefore possess apparently elevated thermal conductivities relative to single crystals of the same n.

2D perovskites↗
DOE OSTI2022

Understanding the Impacts of Support–Polymer Interactions on the Dynamics of Poly(ethyleneimine) Confined in Mesoporous SBA-15

Supported amines are a promising class of CO 2 sorbents offering large uptake capacities and fast uptake rates. Additionally, among supported amines, poly(ethyleneimine) (PEI) physically impregnated in the mesopores of SBA-15 silica is widely used. Within these composite materials, the chain dynamics and morphologies of PEI strongly influence the CO 2 capture performance, yet little is known about chain and macromolecule mobility in confined pores. Here, we probe the impact of the support–PEI interactions on the dynamics and structures of PEI at the support interface and the corresponding impact on CO 2 uptake performance, which yields critical structure–property relationships. The pore walls of the support are grafted with organosilanes with different chemical end groups to differentiate interaction modes (spanning from strong attraction to repulsion) between the pore surface and PEI. Combinations of techniques, such as quasi-elastic neutron scattering (QENS), 1 H T 1 –T 2 relaxation correlation solid-state NMR, and molecular dynamics (MD) simulations, are used to comprehensively assess the physical properties of confined PEI. We hypothesized that PEI would have faster dynamics when subjected to less attractive or repulsive interactions. However, we discover that complex interfacial interactions resulted in complex structure–property relationships. Indeed, both the chain conformation of the surface-grafted chains and of the PEI around the surface influenced the chain mobility and CO 2 uptake performance. By coupling knowledge of the dynamics and distributions of PEI with CO 2 sorption performance and other characteristics, we determine that the macroscopic structures of the hybrid materials dictate the first rapid CO 2 uptake, and the rate of CO 2 sorption during the subsequent gradual uptake stage is determined by PEI chain motions that promote diffusive jumps of CO 2 through PEI-packed domains.

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

Design Principles for Using Amphiphilic Polymers To Create Microporous Water

Aqueous dispersions of microporous nanocrystals with dry, gas-accessible pores—referred to as “microporous water”—enable high densities of gas molecules to be transported through water. For many applications of microporous water, generalizable strategies are required to functionalize the external surface of microporous particles to control their dispersibility, stability, and interactions with other solution-phase components—including catalysts, proteins, and cells—while retaining as much of their internal pore volume as possible. Here, in this work, we establish design principles for the noncovalent surface functionalization of hydrophobic metal–organic frameworks with amphiphilic polymers that render the particles dispersible in water and enhance their hydrolytic stability. Specifically, we show that block co-polymers with persistence lengths that exceed the micropore aperture size of zeolitic imidazolate frameworks (ZIFs) can dramatically enhance ZIF particle dispersibility and stability while preserving porosity and >80% of the theoretical O 2 carrying capacity. Moreover, enhancements in hydrolytic stability are greatest when the polymer can form strong bonds to exposed metal sites on the external particle surface. More broadly, our insights provide guidelines for controlling the interface between polymers and metal–organic framework particles in aqueous environments to augment the properties of microporous water.

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

Enabling the Application of Large Footprint Open-Bottom Permeameters through New Shape Factors

Open-bottom permeameters (OBP) are commonly used for estimating the hydraulic properties of sediment beds under water bodies and fluxes across the sediment-water interface. The specific shape factor for an OBP flow situation is necessary for data interpretation. It is determined by the ratio of the permeameter tube radius to the tube's penetration depth within the sediment. Previous permeameter designs favored small radii which allowed for deeper penetration and stability of the tube. The shape factor is provided here for large footprint OBPs that allows for shallow penetration and greater radius. Through a numerical and analytical modeling study of the OBP hydrodynamics, previous knowledge on the OBP shape factor is generalized to the entire range of radius-to-penetration ratios that might be practically applied. The results reveal the advantages of the large footprint design, including better sensitivity due to the ability to impose larger discharge over a shorter test duration.

42 ENGINEERING↗
DOE OSTI2021

New material platform for superconducting transmon qubits with coherence times exceeding 0.3 milliseconds

The superconducting transmon qubit is a leading platform for quantum computing and quantum science. Building large, useful quantum systems based on transmon qubits will require significant improvements in qubit relaxation and coherence times, which are orders of magnitude shorter than limits imposed by bulk properties of the constituent materials. This indicates that relaxation likely originates from uncontrolled surfaces, interfaces, and contaminants. Previous efforts to improve qubit lifetimes have focused primarily on designs that minimize contributions from surfaces. However, significant improvements in the lifetime of two-dimensional transmon qubits have remained elusive for several years. Here, we fabricate two-dimensional transmon qubits that have both lifetimes and coherence times with dynamical decoupling exceeding 0.3 milliseconds by replacing niobium with tantalum in the device. We have observed increased lifetimes for seventeen devices, indicating that these material improvements are robust, paving the way for higher gate fidelities in multi-qubit processors.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Structural and dynamic mechanisms of GABA A receptor modulators with opposing activities

γ-Aminobutyric acid type A (GABA A ) receptors are pentameric ligand-gated ion channels abundant in the central nervous system and are prolific drug targets for treating anxiety, sleep disorders and epilepsy. Diverse small molecules exert a spectrum of effects on γ-aminobutyric acid type A (GABA A ) receptors by acting at the classical benzodiazepine site. They can potentiate the response to GABA, attenuate channel activity, or counteract modulation by other ligands. Structural mechanisms underlying the actions of these drugs are not fully understood. Here we present two high-resolution structures of GABA A receptors in complex with zolpidem, a positive allosteric modulator and heavily prescribed hypnotic, and DMCM, a negative allosteric modulator with convulsant and anxiogenic properties. These two drugs share the extracellular benzodiazepine site at the α/γ subunit interface and two transmembrane sites at β/α interfaces. Structural analyses reveal a basis for the subtype selectivity of zolpidem that underlies its clinical success. Molecular dynamics simulations provide insight into how DMCM switches from a negative to a positive modulator as a function of binding site occupancy. Together, these findings expand our understanding of how GABA A receptor allosteric modulators acting through a common site can have diverging activities.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Interplay of negative electronic compressibility and capacitance enhancement in lightly-doped metal oxide Bi 0.95 La 0.05 FeO 3 by quantum capacitance model

Light-sensitive capacitance variation of Bi 0.95 La 0.05 FeO 3 (BLFO) ceramics has been studied under violet to UV irradiation. The reversible capacitance enhancement up to 21% under 405 nm violet laser irradiation has been observed, suggesting a possible degree of freedom to dynamically control this in high dielectric materials for light-sensitive capacitance applications. By using ultraviolet photoemission spectroscopy (UPS), we show here that exposure of BLFO surfaces to UV light induces a counterintuitive shift of the O 2p valence state to lower binding energy of up to 243 meV which is a direct signature of negative electronic compressibility (NEC). A decrease of BLFO electrical resistance agrees strongly with the UPS data suggesting the creation of a thin conductive layer on its insulating bulk under light irradiation. By exploiting the quantum capacitance model, we find that the negative quantum capacitance due to this NEC effect plays an important role in this capacitance enhancement.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Strong Fermi-level pinning at metal contacts to halide perovskites

The performance of halide perovskite-based electronic and optoelectronic devices is often related to interfacial charge transport. To shed light on the underlying physical and chemical properties of CH 3 NH 3 PbI 3 (MAPbI 3 ) in direct contact with common electrodes Al, Ti, Cr, Ag, and Au, the evolution of interfacial properties and Fermi level pinning is systematically studied. Given a unique experimental facility, pristine interfaces without any exposure to ambient air were prepared. We observe aggregation of substantial amounts of metallic lead (Pb 0 ) at the metal/MAPbI 3 interface, resulting from the interfacial reaction between the deposited metal and iodine ions from MAPbI 3 . It is found that the Schottky barrier height at the metal/MAPbI 3 interface is independent of the metal work function due to strong Fermi level pinning, possibly due to the metallic Pb 0 aggregates, which act as interfacial trap sites. The charge neutrality level of MAPbI 3 is consistent with the energy level of Pb 0 -related defects, indicating that Pb 0 interfacial trap states can be nonradiative recombination sites. Here this work underlines that control of chemical bonding at interfaces is a key factor for designing future halide perovskite-based devices.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Modeled sensitivity of multi-MA accelerator performance to electrode contaminant inventory

Significant particle-in-cell code development has enabled simulations of power flow in multi-MA accelerators to include the desorption of surface contaminants, their ionization into surface plasmas, and the impact of these plasmas on efficiency. The simulations base desorption on an Arrhenius equation, whose most significant unknown is the surface contaminant inventory. The sensitivity of power-flow simulations to this inventory is studied here using Sandia National Laboratories' Z accelerator with a 7-nH MagLIF load [Phys. Plasmas 17, 056303 (2010)]. Simulations are conducted in 3D cylindrical coordinates for the current-adder, or “convolute,” region of Z and in 2D for the final feed only. Simulated contaminant inventories are varied from 1 to 32 monolayers (MLs) in 2D, and 2 to 4 ML in 3D. The results reveal sensitivities to the local ratio of E/B⁠. The high B-field, low E-field region near the short-circuit load is insensitive to the contaminant inventory, where assumed values of 4–32 ML change the load current by ≤ 2%, and agree with experiment to within 2% at peak current. A 1-ML value is the outlier, increasing the load current by 5%, but still within measurement uncertainty. In contrast, the relatively higher E-field, lower B-field convolute region has slower contaminant desorption and higher-magnitude E-field penetration of the surface plasmas. The current loss in the convolute region does increase with contaminant inventory. The loss assuming 4 ML is 12% larger than for 2 ML, with 4 ML being the better match to experiment.

Arrhenius equation↗
DOE OSTI2022

The effect of metal–insulator interface interactions on electrical transport in granular metals

Here, we present an in-depth study of metal–insulator interfaces within granular metal (GM) films and correlate their interfacial interactions with structural and electrical transport properties. Nominally 100 nm thick GM films of Co and Mo dispersed within yttria-stabilized zirconia (YSZ), with volumetric metal fractions (φ) from 0.2–0.8, were grown by radio frequency co-sputtering from individual metal and YSZ targets. Scanning transmission electron microscopy and DC transport measurements find that the resulting metal islands are well-defined with 1.7–2.6 nm average diameters and percolation thresholds between φ = 0.4–0.5. The room temperature conductivities for the φ = 0.2 samples are several orders of magnitude larger than previously-reported for GMs. X-ray photoemission spectroscopy indicates both oxygen vacancy formation within the YSZ and band-bending at metal–insulator interfaces. The higher-than-predicted conductivity is largely attributed to these interface interactions. In agreement with recent theory, interactions that reduce the change in conductivity across the metal–insulator interface are seen to prevent sharp conductivity drops when the metal concentration decreases below the percolation threshold. These interface interactions help interpret the broad range of conductivities reported throughout the literature and can be used to tune the conductivities of future GMs.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Effect of buffer termination on intermixing and conductivity in LaTiO 3 /SrTiO 3 heterostructures integrated on Si(100)

Here, the control of chemical exchange across heterointerfaces formed between ultrathin functional transition-metal oxide layers provides an effective route to manipulate the electronic properties of these systems. By determining the layer-resolved structural profile across the interface between the Mott insulator, LaTiO 3 (LTO) grown epitaxially on SrTiO 3 (STO)-buffered silicon by molecular beam epitaxy, we find that interfacial cationic exchange depends on the surface termination of the strained STO buffer. Using a combination of temperature-dependent transport and synchrotron x-ray crystal truncation rods and reciprocal space mapping, an enhanced conductivity in STO/LTO/SrO-terminated STO buffers compared to heterostructures with TiO 2 -terminated STO buffers is correlated with La/Sr exchange and the formation of metallic La 1-x Sr x TiO 3 . La/Sr exchange effectively reduces the strain energy of the system due to the large lattice mismatch between the nominal oxide layers and the Si substrate.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Modeling diffusion and depletion in high-aspect-ratio atomic layer deposition processes: Process parameters and manufacturing impacts

Atomic layer deposition (ALD) is a powerful technique for modifying the surface chemistry and properties of substrates with complex and nonplanar topologies. However, achieving uniform and conformal deposition on ultrahigh-aspect-ratio substrates remains challenging, typically requiring large quantities of precursors and long exposure times. Furthermore, process optimization is often performed empirically and involves substantial trial and error. In this work, we perform a combined experimental and computational study of ALD Al 2 O 3 infiltration into silica aerogel monoliths (aspect ratio >10 5 ). A reaction-diffusion model is used to explore the effects of key processing parameters, namely, exposure time per dose, precursor source temperature, number of aerogels in the reactor, and reactor volume. The model is based on quasi-static mode ALD, where the dosed precursor is held in the chamber for a fixed period of time before purging. We analyze the trade-offs between process throughput and precursor utilization for each of these parameters. Furthermore, we investigate the co-optimization and interactions between multiple process parameters, demonstrating the potential for further improvements. Furthermore, this physics-based model can be used to identify a set of process parameters for high-aspect-ratio ALD that meet specific manufacturing objective functions, including throughput, cost, and sustainability.

Aerogel↗
DOE OSTI2024

Fluidized-Bed Gasification of Coal-Biomass-Plastics for Hydrogen Production

Coal is one of the most abundant fossil energy resources in the United States and in the world. The recoverable reserves in the United States are estimated to be about 252 billion tons – more than 350 years of supply at current rates of usage. However, the share of coal in total primary energy consumption in the US has been declining over the years. The decline of coal is mainly attributed to cheap natural gas and precipitous declines in the cost of electricity production from renewable technologies such as wind and solar. Coal can potentially be used if it is coupled with carbon-neutral feedstock such as biomass-agricultural residues, forest biomass, and forest residues. Co-gasification of coal and biomass can become a negative carbon emission technology if the carbon dioxide (CO 2 ) is captured and sequestered. Although the biomass gasification process has a lot of similarities to coal gasification, the large-scale adaptation of power production sourced from biomass has not come to fruition. The main reason is that the power production from biomass is still expensive when compared with natural gas or coal power technologies. To address the feedstock cost, one approach is to use low-cost feedstocks, such as municipal solid wastes (MSW) or plastics, for gasification. Gasification involves the partial oxidation of coal and/or biomass feedstocks to produce a combustible fuel called synthesis gas (syngas) which is composed of carbon monoxide (CO), hydrogen (H 2 ), CO 2 , methane (CH 4 ), nitrogen (N 2 ), water (H 2 O), and other compounds that might be considered as contaminants. Raw syngas from gasification must go through multiple steps to produce high-purity hydrogen. The specific steps depend upon the quality (gas composition, contaminants, and their concentration) and condition (pressure and temperature) of syngas. The long-term goal of the project was to utilize coal and plastics together with biomass to produce energy and fuels using a gasification platform while reducing greenhouse gas emissions. The main objective of this research was to examine gasification performance in a laboratory-scale fluidized-bed gasifier for hydrogen production. The specific objectives of the research were to: (i) study coal-plastic-biomass mixture flowability for consistent feeding in the gasifier; (ii) understand the gasification behavior of the mixture in steam and oxygen environments; (iii) perform thermal property characterization of ash and slag from the mixture feedstock and refractory-ash interface of the mixture under gasification conditions; and (iv) develop process models to determine the technology needed for syngas cleanup and contaminants. The study found that there was no apparent segregation when biomass, coal, and waste plastics were mixed together during feeding. Although there were differences in hydrogen production when individual feedstock were fed, the hydrogen concentration remained almost constant with various blends. Therefore, blending waste plastics with biomass and coal, which are all abundant, is a better approach for energy production. Results of the techno-economic analysis suggested that integration of advanced gasification (GTI’s R-GAS™) and syngas cleanup and conditioning technologies (RTI’s WDP and AFWGS) for clean hydrogen production resulted in substantial benefits, including significant capital cost and operating cost reductions. Advanced technologies resulted in 16% reduction in the hydrogen production cost (COH) from 2.94 $\$$/kg to 2.47 $\$$/g, with further scope for optimization and cost reduction. These advanced technologies also result in lower emissions, and improved energy efficiency.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2023

Neural Network-Enhanced Reproducing Kernel Particle Method for Image-Based Multiphysics Damage Modeling of Energy Storage Materials

Energy storage materials undergo significant stresses during charge/discharge cycling, which makes understanding their reliability and durability fundamental in predicting performance and service life. Strong electrochemical-mechanical coupling and highly anisotropic material properties contribute to the formation and propagation of micro-cracking, largely along material interfaces and grain boundaries. With microstructural images supplied by the National Renewable Energy Laboratory (NREL), image-based modeling techniques are used to represent the complex material microstructures that dictate the coupled physics of these systems. Traditional electrochemical-mechanical models rely on mesh-based finite element methods, which can lead to difficulties in capturing crack propagation due to mesh dependency. Additionally, commonly used damage models, such as the continuous damage model and the cohesive zone model, often have steep tradeoffs between discontinuous field accuracy and computational expense. In this work, a neural network-enhanced reproducing kernel particle method (NN-RKPM) [1] is leveraged to accurately capture damage and crack propagation throughout the material by learning the location, orientation, and sharpness of discontinuity while allowing for a coarser nodal distribution than that necessary for capturing sharp solution transitions using traditional mesh-based methods. NN-RKPM is used to inform how crack opening and closure in turn affect the coupled chemical equations and material microstructure. Reference: [1] Baek, J., Chen, J. S., Susuki, K., "Neural Network enhanced Reproducing Kernel Particle Method for Modeling Localizations," International Journal for Numerical Methods in Engineering, Vol. 123, pp 4422-4454, https://doi.org/10.1002/nme.7040, 2022.

damage modeling↗
DOE OSTI2025

Entrapment Behavior of Solid Surrogate Fission Products at Engineered UN Nano‐Hetero‐Interfaces Within Metallic Nuclear Fuels

Nanometric hetero-interfaces provide a wealth of scientific and engineering opportunities due to their complex and often misunderstood properties that can differ from their respective bulk constituents. In this work, the ability for engineered nanostructures within a bulk U─Mo alloy to arrest simulant fission products is investigated experimentally and computationally. Nanostructured 90 wt% U/ 10 wt% Mo (U-10Mo) with 7.1 at% Nd is consolidated using spark-plasma- sintering (SPS) techniques and is heat-treated at 500 °C under vacuum for 24, 100, 500, and 1000 h. Analysis on the sintered and heat-treated U-10Mo reveals rapid kinetics in Nd diffusion to nanocluster sites, with evidence of Nd diffusion occurring during sintering and during the following heat-treatment. The segregation behavior of Nd at two different U─Mo/UN interfaces is computationally verified using density functional theory (DFT) to reinforce experimental data.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Beyond Local Solvation Structure: Nanometric Aggregates in Battery Electrolytes and their Effect on Electrolyte Properties

Electrolytes are an essential component of all electrochemical storage and conversion devices, such as batteries. In the history of battery development, the complex nature of electrolytes has often been a bottleneck. Fundamental knowledge of electrolyte systems encompasses elucidation of structure-property relationships of the solution species. Recently, nanometric aggregates have been observed in several classes of electrolytes, including super-concentrated, redox-flow, multivalent, polymer, and ionic liquid-based electrolytes. Compared with the well-studied local solvation structures such as contact ion pairs and solvent-separated ions, these aggregates impose unique effects on the ion distribution and transport both within bulk electrolytes and at electrode/electrolyte interfaces. This Perspective highlights the discovery of the aggregates in various battery electrolytes and their impact on electrolyte properties. We also present an outlook for future studies of this emerging field of nanometric aggregates and the need for the development of new experimental and computational tools to study their properties.

Yu, Zhou↗
DOE OSTI2023

Distance-dependent dielectric constant at the calcite/electrolyte interface: Implication for surface complexation modeling

Hypothesis: The electrical double layer formed at the mineral/electrolyte interface is often modeled using mean-field approaches based on a continuum description of the solvent whose dielectric constant is assumed to decrease monotonically with decreasing distance to the surface. In contrast, molecular simulations show that the solvent polarizability oscillates near the surface similar to the water density profile - as shown previously, for example, by Bonthuis et al. (D.J. Bonthuis, S. Gekle, R.R. Netz, Dielectric Profile of Interfacial Water and its Effect on Double-Layer Capacitance, Phys Rev Lett 107(16) (2011) 166102). Here, we showed that molecular and mesoscale pictures agree by spatially averaging the dielectric constant obtained from molecular dynamics simulations over the distances relevant to the mean-field representation. In addition, the values of capacitances used to describe the electrical double layer in Surface Complexation Models (SCMs) of the mineral/electrolyte interface can be estimated using molecularly informed spatially averaged dielectric constants and positions of hydration layers. Experiments: First, we used molecular dynamics simulations to model the calcite 101¯4/electrolyte interface. Next, by using atomistic trajectories, we calculated the distance-dependent static dielectric constant and water density in the direction normal to the. Finally, we applied spatial compartmentalization consistent with the model of parallel-plate capacitors connected in series to estimate SCM capacitances. Findings: Computationally expensive simulations are required to determine the dielectric constant profile of interfacial water near the mineral surface. On the other hand, water density profiles are readily assessable from much shorter simulation trajectories. Our simulations confirmed that dielectric and water density oscillations at the interface are correlated. Here, we parametrized linear regression models to estimate the dielectric constant directly from the local water density. This is a significant computational shortcut compared to slowly converging calculations relying on total dipole moment fluctuations. The amplitude of the interfacial dielectric constant oscillation can exceed the dielectric constant of the bulk water, suggesting an ice-like frozen state, but only if there are no electrolyte ions. The interfacial accumulation of electrolyte ions causes a decrease in the dielectric constant due to the reduction of water density and re-orientation of water dipoles in ion hydration shells. Finally, we show how to use the computed dielectric properties to estimate SCM's capacitances.

58 GEOSCIENCES↗