Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “control layers”

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

Identifying higher oxygenate synthesis sites in Cu catalysts promoted and stabilized by atomic layer deposited Fe 2 O 3

The conversion of syngas to higher oxygenates is an important step in the pathway towards the long-term sustainable production of fuels and chemicals. Mixed Cu/Fe catalysts combine metals with different CO activation characteristics, a property that makes them promising candidates for higher oxygenate production. However, Cu/Fe catalysts suffer from instability due to sintering and phase separation, and the active sites for higher oxygenate formation on these materials are not well understood. In this work, we use atomic layer deposition (ALD) to modify silica-supported Cu nanoparticles with monolayer precise amounts of Fe 2 O 3 and study the effects on catalyst structure, reactivity, and stability. We demonstrate that the Fe 2 O 3 ALD process is inherently selective towards deposition on the SiO 2 surface versus the Cu nanoparticles. As a result, a sufficiently thick Fe 2 O 3 layer surrounding the Cu nanoparticles acts as a barrier against Cu migration, preventing catalyst deactivation through sintering without covering the Cu. In syngas conversion reactions, we find that higher oxygenate selectivity increases with Fe 2 O 3 loading at low and intermediate ALD cycle numbers. Furthermore, using in situ X-ray absorption spectroscopy, we demonstrate that the Fe near the Cu nanoparticles is reduced to a metallic state under reaction conditions. From these measurements, we identify these metallic Cu/Fe interfaces as the active sites for higher oxygenate formation, which are most abundant at intermediate Fe 2 O 3 ALD cycle numbers. As the Fe 2 O 3 loading increases, higher oxygenate selectivity falls, resulting from encapsulation of the Cu nanoparticles by thick Fe 2 O 3 layers. Overall, the results demonstrate how thickness-controlled ALD Fe 2 O 3 can be used as both a promoter and stabilizer of Cu-based higher oxygenate synthesis catalysts.

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

Nanoscale electrostatic control in ultraclean van der Waals heterostructures by local anodic oxidation of graphite gates

In an all-van der Waals heterostructure, the active layer, gate dielectrics and gate electrodes are assembled from two-dimensional crystals that have a low density of atomic defects. This design allows two-dimensional electron systems with very low disorder to be created, particularly in heterostructures where the active layer also has intrinsically low disorder, such as crystalline graphene layers or metal dichalcogenide heterobilayers. A key missing ingredient has been nanoscale electrostatic control, with existing methods for fabricated local gates typically introducing unwanted contamination. Here we describe a resist-free local anodic oxidation process for patterning sub-100 nm features in graphite gates, and their subsequent integration into an all-van der Waals heterostructure. We define a quantum point contact in the fractional quantum Hall regime as a benchmark device and observe signatures of chiral Luttinger liquid behaviour, indicating an absence of extrinsic scattering centres in the vicinity of the point contact. In the integer quantum Hall regime, we demonstrate in situ control of the edge confinement potential, a key requirement for the precision control of chiral edge states. In conclusion, this technique may enable the fabrication of devices capable of single anyon control and coherent edge-state interferometry in the fractional quantum Hall regime.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Impact of Flow Configurations on Response Time and Data Quality in Real-Time, In-Line Fourier Transform Infrared (FTIR) Monitoring of Viscous Flows

The real-time, in-line monitoring of continuous flow concentrations is widely conducted via infrared (IR) spectrometry by using a flow cell connected to a reactive flow stream. For protective purposes, the IR sensor tip is typically offset from the flow. This offset can cause the formation of a stagnant boundary layer above the sensor, especially when dealing with high-viscosity fluids. As a result, the IR signal response time is often controlled by the slow diffusional exchange of fluid in the boundary layer, as confirmed via 2D computational fluid dynamics (CFD) simulations. We evaluated several flow configuration modifications in a typical IR flow cell in order to identify the changes to the flow dynamics that enable improved response times with minimal changes to the cell configuration: the use of (i) vertical flow, where the standard horizontal flow over the sensor is redirected to contact vertically with the sensor, (ii) a static mixer to create radial flow momentum above the IR sensor, and (iii) horizontal or vertical nozzles to direct the flow toward the IR sensor. The vertical flow configuration did not show any significant improvement over the standard horizontal flow configuration. However, the static mixer, horizontal nozzle, and vertical nozzle configurations all resulted in markedly improved response times and signal quality, albeit at the expense of a higher pressure drop across the flow cell. These results point toward straightforward, user-accessible modifications of in-line IR flow cells that result in significant improvements in signal stability and acquisition times.

42 ENGINEERING↗
DOE OSTI2024

Aerosol‐Correlated Cloud Activation for Clean Conditions in the Tropical Atlantic Boundary Layer During LASIC

Abstract Aerosol measurements during the DOE ARM Layered Atlantic Smoke Interactions with Clouds (LASIC) campaign were used to quantify the differences between clean and smoky cloud condensation nuclei (CCN) budgets. Accumulation‐mode particles accounted for ∼70% of CCN at supersaturations <0.3% in clean and smoky conditions. Aitken‐mode particles contributed <20% and sea‐spray‐mode particles <10% at supersaturations <0.3%, but at supersaturations >0.3% Aitken particles contributions increased to 30%–40% of clean CCN. For clean conditions, the Hoppel minimum diameter was correlated to the accumulation‐mode number concentration, indicating aerosol‐correlated cloud activation was controlling the lower diameter cutoff for which particles serve as CCN. For smoky conditions, the contributions of Aitken particles increase and the correlation of cloud activation to accumulation‐mode particles is masked by the lower‐hygroscopicity smoke. These results provide the first multi‐month in situ quantitative constraints on the role of aerosol number size distributions in controlling cloud activation in the tropical Atlantic boundary layer.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Local cation ordering in compositionally complex Ruddlesden–Popper n = 1 oxides

The Ruddlesden–Popper (RP) layered perovskite structure is of great interest due to its inherent tunability, and the emergence and growth of the compositionally complex oxide (CCO) concept endows the RP family with further possibilities. Here, a comprehensive assessment of thermodynamic stabilization, local order/disorder, and lattice distortion was performed in the first two reported examples of lanthanum-deficient La n+1 B n O 3n+1 (n = 1, B = Mg, Co, Ni, Cu, Zn) obtained via various processing conditions. Chemical short-range order (CSRO) at the B-site and the controllable excess interstitial oxygen (δ) in RP-CCOs are uncovered by neutron pair distribution function analysis. Reverse Monte Carlo analysis of the data, Metropolis Monte Carlo simulations, and extended x-ray absorption fine structure analysis implies a modest degree of magnetic element segregation on the local scale. Further, ab initio molecular dynamics simulations results obtained from special quasirandom structure disagree with experimentally observed CSRO but confirm Jahn–Teller distortion of CuO 6 octahedra. These findings highlight potential opportunities to control local order/disorder and excess interstitial oxygen in layered RP-CCOs and demonstrate a high degree of freedom for tailoring application-specific properties. They also suggest a need for expansion of theoretical and data modeling approaches in order to meet the innate challenges of CCO and related high-entropy phases.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Experimental validation of a high fidelity Monte Carlo neutron transport model of the MIT graphite exponential pile

High-fidelity modeling and simulation were performed for the MIT graphite exponential pile (MGEP) using Monte Carlo neutron transport codes OpenMC and MCNP, and the results were validated by experimental data. The MGEP is being used as the test bed for the design of an autonomous control system for the pile's neutron flux distribution. The main contribution of this work is to generate the training data sets of neutron flux distributions with different locations of control rods that perturb the neutron flux profiles. First, code -to-code cross verification between OpenMC and MCNP was performed to ensure consistency of the numerical modeling within statistical uncertainties. To validate the accuracy of this high-fidelity model, a series of neutron flux measurements were conducted using a Helium-3 (He-3) neutron detector on a mobile platform that is placed inside the pile. Second, the neutron flux profiles were measured in four vertical layers of interest, and compared to the corresponding simulation results. The comparison results shows that the root mean square error is less than 2.5% in the two upper layers, and less than 4.5% in all four measured layers. Here the results validated the accuracy of the modeling and simulation. Finally, the relative change of the neutron flux profiles from moving control rods was analyzed, which identified the layer that has the best sensitivity regarding the control rods movements. Thus, this work identified and provided training data sets of both simulated and experimental neutron flux profiles in the most sensitive layer, paving the path forward to the real-time experimental demonstration of the autonomous control system.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Tuning 3-D Nanomaterial Architectures Using Atomic Layer Deposition to Direct Solution Synthesis

The ability to synthesize nanoarchitected materials with tunable geometries provides a means to control their functional properties, with applications in biological, environmental, and energy fields. To this end, various bottom-up and top-down synthesis processes have been developed. However, many of these processes require prepatterning or etching steps, making them challenging to scale-up to complex, nonplanar substrates. Furthermore, the ability to integrate nanomaterials into hierarchical arrays with precise control of feature spacing and orientation remains a challenge. One approach to overcome these patterning challenges is the use of surface modification layers to guide the resulting geometry of nanomaterial architectures grown from the substrate. A powerful strategy to accomplish this is what we will refer to as “surface-directed assembly,” where the resulting geometric parameters (feature size, shape, orientation) are predetermined by the initial surface layer. In particular, the use of Atomic Layer Deposition (ALD) to form a surface layer, followed by solution-based growth processes, has the ability to synthesize architected structures with tunable geometries on complex, nonplanar surfaces. Over the past decade, we have reported a series of studies where surface-directed assembly is used to synthesize ZnO nanowires (NWs) on top of a variety of substrates. In this case, a thin film of ZnO is deposited onto the substrate using ALD, which can guide the NW diameter, spacing, and angular orientation with respect to the substrate by controlling epitaxial relationships. Furthermore, we have shown that by depositing a submonolayer overcoat of a secondary material (e.g., amorphous TiO 2 ), nucleation sites are partially blocked, which can further tune the spacing between nanowires while minimizing changes to their other geometric properties. This approach can be used to generate multilevel hierarchical structures, such as hyperbranched NW arrays with tunable control of each level of hierarchy using ALD. Finally, we have demonstrated that the tunable control of geometric parameters can be scaled-up to curved, nonplanar substrates. This highlights the power of ALD to conformally and uniformly deposit the seed layers on complex substrates with subnanometer precision. To complement these seeded hydrothermal approaches, we expanded this strategy to include conversion chemistry of the initial ALD seed layers. For example, by replacing ZnO with Al 2 O 3 as the seed layer without changing the hydrothermal growth conditions, Al–Zn layered-double hydroxide nanosheets can be formed instead of nanowires. In another example of conversion chemistry, a solution anion-exchange process was used to incorporate sulfur into ALD metal oxide films. In both of these conversion processes, the properties of the initial ALD film enabled tuning of the resulting nanostructure geometry. In this Account, we describe the use of ALD to guide the growth of diverse nanomaterial systems, with tunable control over their geometry and composition. We further show how these approaches can be used to tune functional properties for a range of applications, including superomniphobic surfaces, antibiofouling coatings, and photocatalysis. In conclusion, we conclude with an outlook on how the combination of ALD and solution synthesis can enable future directions in scalable nanomanufacturing to overcome the limitations of traditional top-down and bottom-up approaches.

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

Phase transition mechanism and bandgap engineering of Sb2S3 at gigapascal pressures

Abstract Earth-abundant antimony trisulfide (Sb 2 S 3 ), or simply antimonite, is a promising material for capturing natural energies like solar power and heat flux. The layered structure, held up by weak van-der Waals forces, induces anisotropic behaviors in carrier transportation and thermal expansion. Here, we used stress as mechanical stimuli to destabilize the layered structure and observed the structural phase transition to a three-dimensional (3D) structure. We combined in situ x-ray diffraction (XRD), Raman spectroscopy, ultraviolet-visible spectroscopy, and first-principles calculations to study the evolution of structure and bandgap width up to 20.1 GPa. The optical band gap energy of Sb 2 S 3 followed a two-step hierarchical sequence at approximately 4 and 11 GPa. We also revealed that the first step of change is mainly caused by the redistribution of band states near the conduction band maximum. The second transition is controlled by an isostructural phase transition, with collapsed layers and the formation of a higher coordinated bulky structure. The band gap reduced from 1.73 eV at ambient to 0.68 eV at 15 GPa, making it a promising thermoelectric material under high pressure.

Cui, Zhongxun (ORCID:0000000275044177)↗
DOE OSTI2022

Steering CO 2 hydrogenation toward C–C coupling to hydrocarbons using porous organic polymer/metal interfaces

The conversion of CO 2 into fuels and chemicals is an attractive option for mitigating CO 2 emissions. Controlling the selectivity of this process is beneficial to produce desirable liquid fuels, but C–C coupling is a limiting step in the reaction that requires high pressures. Here, we propose a strategy to favor C–C coupling on a supported Ru/TiO 2 catalyst by encapsulating it within the polymer layers of an imine-based porous organic polymer that controls its selectivity. Such polymer confinement modifies the CO 2 hydrogenation behavior of the Ru surface, significantly enhancing the C 2+ production turnover frequency by 10-fold. We demonstrate that the polymer layers affect the adsorption of reactants and intermediates while being stable under the demanding reaction conditions. Our findings highlight the promising opportunity of using polymer/metal interfaces for the rational engineering of active sites and as a general tool for controlling selective transformations in supported catalyst systems.

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

Cyber–Physical System Security of Distribution Systems

The Information and Communications Technology (ICT) for control and monitoring of power systems is a layer on top of the physical power system infrastructure. The cyber system and physical power system components form a tightly coupled Cyber–Physical System (CPS). Sources of vulnerabilities arise from the computing and communication systems of the cyber–power grid. Cyber intrusions targeting the power grid are serious threats to the reliability of electricity supply that is critical to society and the economy. In a typical Information Technology environment, numerous attack scenarios have shown how unauthorized users can access and manipulate protected information from a network domain. The need for cyber security has led to industry standards that power grids must meet to ensure that the monitoring, operation, and control functions are not disrupted by cyber intrusions. Cyber security technologies such as encryption and authentication have been deployed on the CPS. Intrusion or anomaly detection and mitigation tools developed for power grids are emerging. Furthermore, this survey paper provides the basic concepts of cyber vulnerabilities of distribution systems and CPS security. The important ICT subjects for distribution systems covered in this paper include Supervisory Control And Data Acquisition, Distributed Energy Resources, including renewable energy and smart meters.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Nanostructured polyelectrolytes for ion-selective membranes

Nanostructured polyelectrolyte bilayers deposited by Layer-by-Layer deposition on nanoporous membranes can be selectively crosslinked to modify the polyelectrolyte charge density and control ionic selectivity independent of ionic conductivity. For example, the polyelectrolyte bilayer can comprise a cationic polymer layer, such as poly(ethyleneimine), and an anionic polymer layer, such as poly(acrylic acid). Increasing the number of bilayers increases the cation selectivity when the poly(ethyleneimine) layer is crosslinked with glutaraldehyde. Crosslinking the membranes also increases the chemical and mechanical strength of the polyelectrolyte films. This controllable and inexpensive method can be used to create ion-selective and mechanically robust membranes on porous supports for a wide range of applications.

25 ENERGY STORAGE↗
DOE OSTI2021

Molecular dopants: Tools to control the electronic structure of metal halide perovskite interfaces

In the standard configurations of metal halide perovskite solar cell, the active layer, or absorber, follows a p-i-n or n-i-p electronic structure that is designed to enhance the separation and extraction of photo-induced charge carriers. The control of the Fermi level position across the film, between electron and hole transport layers, is therefore of paramount importance. Direct localized doping in metal halide perovskites being still elusive, the design of n-i-p and p-i-n structures has so far relied predominantly on surface and interface doping of the perovskite as well as on the control of the work function of the substrate and transport layers on which, or between which, the absorber is being placed. We provide here a short review of that work, emphasizing the fundamental studies of electronic structure performed on systems modified with organic molecular dopants. The review starts with a justification for the effectiveness of interface doping, based on the ability to move the Fermi level across the gap of the perovskite. We then review work done on the deposition of molecular oxidants and reductants on perovskite surfaces, including the mitigation of the surface states, and the impact of these dopants on energy level alignment with substrate and charge transport layers. The second part of the review focuses on the use of molecular dopants to either modify the work function of electron or hole transport layers to establish the boundary conditions for a p-i-n or n-i-p structure, or to enhance the conductivity of these layers in order to facilitate charge carrier extraction. Final considerations are also given on recent work on bulk doping of the perovskite layer with molecular dopants.

14 SOLAR ENERGY↗
DOE OSTI2021

Tuning electrical and interfacial thermal properties of bilayer MoS2 via electrochemical intercalation

Layered two-dimensional (2D) materials such as MoS 2 have attracted much attention for nano- and opto-electronics. Recently, intercalation (e.g. of ions, atoms, or molecules) has emerged as an effective technique to modulate material properties of such layered 2D films reversibly. We probe both the electrical and thermal properties of Li-intercalated bilayer MoS 2 nanosheets by combining electrical measurements and Raman spectroscopy. Furthermore, we demonstrate reversible modulation of carrier density over more than two orders of magnitude (from 0.8 × 10 12 to 1.5 × 10 14 cm -2 ), and we simultaneously obtain the thermal boundary conductance between the bilayer and its supporting SiO 2 substrate for an intercalated system for the first time. This thermal coupling can be reversibly modulated by nearly a factor of eight, from 14 ± 4.0 MW m -2 K -1 before intercalation to 1.8 ± 0.9 MW m -2 K -1 when the MoS 2 is fully lithiated. These results reveal electrochemical intercalation as a reversible tool to modulate and control both electrical and thermal properties of 2D layers.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2025

Synthesis of superconducting freestanding infinite-layer nickelate heterostructures on the millimetre scale

The development of superconductivity in infinite-layer nickelates through hole-doping relies on the controlled synthesis of Ni in a high oxidation state, followed by topotactic reduction to a very low oxidation state. So far, superconductivity has been realized only in epitaxial thin films. Here, in this work, we integrate these techniques with heterostructures that include an epitaxial soluble buffer layer, enabling the release of freestanding (Nd,Sr)NiO 2 heterostructures. The released heterostructures exhibit comparable structural and electronic properties to those of optimized thin films, with lateral dimensions ranging from millimetres to ~100 μm, depending on the degree of strain released with respect to the initial substrate. The changes in the superconducting transition temperature upon release of the heterostructure mirror those reported for variations in substrate and pressure, suggesting a common underlying response to strain in the infinite-layer nickelate superconductivity. These freestanding structures will facilitate a range of experimental studies without the constraints of substrates.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Growth of Large‐Sized 2D Ultrathin SnSe Crystals with In‐Plane Ferroelectricity

Abstract Tin (II) selenide (SnSe) is an emerging 2D material with many intriguing properties, such as record‐high thermoelectric figure of merit (ZT), purely in‐plane ferroelectricity, and excellent nonlinear optical properties. To explore these functional properties and related applications, a crucial step is to develop controllable routes to synthesize large‐area, ultrathin, and high‐quality SnSe crystals. Physical vapor deposition (PVD) constitutes a reliable method to synthesize 2D SnSe, however, effects of various growth parameters have not yet been systematically investigated, and current PVD‐synthesized flakes are often thick (>10 nm) with small lateral sizes (<10 µm). In this work, high‐quality 2D SnSe crystals are synthesized via low‐pressure PVD, which display in‐plane ferroelectric domains observed by piezoresponse force microscopy and polarization‐dependent reflection spectroscopy. Detailed studies regarding the roles of various parameters are further carried out, including substrate pre‐annealing, growth duration, temperature, and pressure, which enable to rationally optimize the growth and obtain 2D SnSe crystals with lateral sizes up to ≈23.0 µm and thicknesses down to ≈2.0 nm (3–4 layers). This work paves the way for the controlled growth of large‐area 2D SnSe, facilitating the future exploration of many interesting multiferroic properties and applications with atomic thickness.

2D materials↗
DOE OSTI2023

Density variance dynamics in disparate shock tubes

This report discusses two experiments which investigate thin layer, heavy curtain fragmentation from the perspective of a Reynolds-Averaged mix model, in drastically disparate experimental regimes. The first, the centimeter/millisecond-scale “Horizontal Shock Tube” (HST) is a compressed-gas piston-driven shock tube experiment. The second, the “Multishock thin layer” (Mshock) experiment performed at the National Ignition Facility, is a micrometer/nanosecond-scale laser-driven shock tube experiment. Both are situations in which a heavy plane layer (a ‘curtain’) is initially suspended in a lighter medium. After being shocked from at least one side, the layer translates while its interfaces evolve due to the excitation of the Richtmyer-Meshkov instability at its surfaces. The evolution of density variance, which initially exists only on the surface of the layer, as it comes to encompass the whole layer interior is used as a description of layer fragmentation and dissolution. These experiments have each been simulated in the Los Alamos National Laboratory multi-physics code xRAGE, which includes fundamental hydrodynamics, extended plasma physics and radiation effects which are important to drive the high-energy density experiment, and the Besnard-Harlow-Rauenzahn (BHR) turbulence model. In each, the principal diagnostic for comparison is an experimental metric for the density (co)variance, b, which tracks the moments of the density field at the curtain interfaces and body. Due to experimental constraints in different regimes (i.e. optical diagnostics can be deployed on conventional shock tubes, while the plasma shock tubes must be imaged by x-rays; interfaces can be imposed to specification on laser-driven experiments, which are stored in the solid phase, while conventional experiments have imperfect control of the flow fields which separate the layer, etc.) the experiments are not designed to be perfect scaled cognates of one another. However, despite the separation of six orders of magnitude of scaling in time, and four in space, we are able to demonstrate that the same turbulence model, operating in the same fashion in the same computer code, is able to reproduce results in each experiment, by tracking evolution due to common relevant physics. Additionally, we will present preliminary work toward density variance comparisons in a single-interface Richtmyer-Meshkov configuration, the conventional fluid “Vertical Shock Tube” (VST) experiment, and the Modal Initial Conditions (ModCons) campaign fielded at the OMEGA-EP laser facility.

42 ENGINEERING↗
DOE OSTI2021

SIRIUS: Science-Driven Data Management for Multi-Tiered Storage

The data sets being generated by large applications on very large-scale systems are increasing in both size and complexity. At the same time, there are new ways available to store and access these data sets. The goal in this project is to develop software that applications can use to make use of new and existing storage technologies in more sophisticated ways. One challenge in scientific data management is handling ‘hot’ vs ‘cold’ data. Data that is hot is data that is needed (or will be needed soon) in order for the program to continue progressing, while cold data is either output (and so will not be need further during the life of the program) or will not be needed until significantly later in the program’s run. Hot data should be stored in a way that allows fast access. On most systems, economic factors lead to an inverse relationship between storage performance and storage capacity and so fast access storage is limited. This makes it important to correctly place hot and cold data and avoid cold data unnecessarily consuming precious resources. In this reporting period, we addressed this challenge in various ways and at various levels. Data management frameworks offer only limited control to applications in how data is stored. We have added software capabilities for seamlessly moving data between layers of the storage technology using promote and demote functions to existing software frameworks. This gives direct control to applications in deciding what priority data receives. Additionally, we integrated different storage layer management frameworks in order to allow data to be exchanged and moved between storage layers in a consistent way across the application. Further, applications are not always able to directly decide what storage level makes sense for a given piece of data without an understanding of the underlying storage technologies. Data storage frameworks are often positioned to make these sorts of decisions in service of the application. We have added machine-learning based capabilities to data staging frameworks in order to make intelligent decisions about where data should be stored given learning about patterns in previous usage of similar data.

97 MATHEMATICS AND COMPUTING↗
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↗