Search NASA⌕ Search

SEARCH · Search NASA

Results for “surface stability”

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 1,315 records · Page 73
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Hybrid Doping Strategy with High‐Entropy Cu/Fe Surface Modification and Zr Bulk Incorporation for Ni‐Rich Cathodes

A hybrid doping strategy combining Zr 4+ bulk doping with high-entropy Cu 2+ /Fe 3+ surface doping is developed to enhance the structural and interfacial stability of Ni-rich layered oxide cathodes. Cu and Fe are selectively introduced at the particle surface via a surface-selective ion-exchange process, forming a ≈15 nm Fe-rich layer while preserving the layered framework. Compared to the pristine cathode, the hybrid sample exhibits significantly improved electrochemical performance in both half-cell and full-cell configurations. In half-cells, the hybrid retains 88.5% and 90.2% after 100 cycles at 1C under 4.6 and 4.5 V, respectively. During high-voltage full-cell cycling, the hybrid cathode maintains over 80% capacity retention, whereas the pristine counterpart retains less than 10% under identical conditions over the same cycling period. XPS, EELS, and DEMS analyses confirm improved oxygen retention, suppressed gas evolution, and stable surface chemistry, while DFT calculations indicate enhanced Me–O bonding in the selected Fe 0.75 Cu 0.25 (Mn 1/16 Co 2/16 Ni 13/16 )O 2 surface composition, which is identified through DFT-calculated mixing energy reaching a minimum at this ratio, indicating the most thermodynamically favorable configuration. In conclusion, these results demonstrate the effectiveness of this hybrid doping strategy in mitigating coupled degradation pathways in Ni-rich cathodes.

15 GEOTHERMAL ENERGY↗
NASA NTRS2011

Effects of Cavities and Protuberances on Transition over Hypersonic Vehicles

Surface protuberances and cavities on a hypersonic vehicle are known to cause several aerodynamic or aerothermodynamic issues. Most important of all, premature transition due to these surface irregularities can lead to a significant rise in surface heating. To help understand laminar-turbulent transition induced by protuberances or cavities on a Crew Exploration Vehicle (CEV) surface, high-fidelity numerical simulations are carried out for both types of trips on a CEV wind tunnel model. Due to the large bluntness, these surface irregularities reside in an accelerating subsonic boundary layer. For the Mach 6 wind tunnel conditions with a roughness Reynolds number Re(sub kk) of 800, it was found that a protuberance with a height to boundary layer thickness ratio of 0.73 leads to strong wake instability and spontaneous vortex shedding, while a cavity with identical geometry only causes a rather weak flow unsteadiness. The same cavity with a larger Reynolds number also leads to similar spontaneous vortex shedding and wake instability. The wake development and the formation of hairpin vortices for both protuberance and cavity were found to be qualitatively similar to that observed for an isolated hemisphere submerged in a subsonic, low speed flat-plate boundary layer. However, the shed vortices and their accompanying instability waves were found to be slightly stabilized downstream by the accelerating boundary layer along the CEV surface. Despite this stabilizing influence, it was found that the wake instability spreads substantially in both wall-normal and azimuthal directions as the flow is evolving towards a transitional state. Similarities and differences between the wake instability behind a protuberance and a cavity are investigated. Computations for the Mach 6 boundary layer over a slender cylindrical roughness element with a height to the boundary layer thickness of about 1.1 also shows spontaneous vortex shedding and strong wake instability. Comparisons of detailed flow structures associated with protuberances at subsonic and supersonic edge Mach numbers indicate distinctively different instability mechanisms.

Chang, Chau-Lyan↗
DOE OSTI2024

Insights into the Electrochemical Oxidation and Reduction of Nickel Oxide Surfaces

Surface oxidation/reduction processes, driven by varying electrochemical potentials, can substantially impact catalyst effectiveness and, consequently, electrolyzer performance. Here, this study combines theoretical and experimental approaches to explore the surface redox behavior of nickel oxides, which are cost-effective and efficient catalysts for many electrochemical reactions. Surface Pourbaix diagrams for three different phases of nickel oxides, i.e., nickel hydroxide (Ni(OH) 2 ), nickel oxyhydroxide (NiOOH), and nickel dioxide (NiO 2 ), were constructed using density functional theory-based simulations. Various experimental methods, including cyclic voltammetry, in situ Raman spectroscopy, and electrochemical titration, were employed to probe the surface redox processes of nickel oxide thin films. Our findings indicate that the ABAB stacking sequence of Ni(OH) 2 lacks stability under oxidizing conditions to host the surface oxidation (deprotonation) events, while the AABBCC stacking sequence of NiOOH is energetically favorable due to the presence of interlayer hydrogen bonding. Rapid charge transfer facilitated by interlayer hydrogen bonding accounts for the higher reactivity of partially oxidized/reduced NiOOH (001) surfaces compared to Ni(OH) 2 (001) and NiO 2 (001) surfaces with the same stoichiometry, where interlayer hydrogen bonding is absent. Insights presented in this work can offer guidelines for optimizing operational conditions and tailoring the surface structures and oxidation states of nickel oxides to enhance performance in applications such as electrocatalysis and supercapacitors.

36 MATERIALS SCIENCE↗
NASA NTRS1976

Atomic oxygen-metal surface studies as applied to mass spectrometer measurements of upper planetary atmospheres

The problem of atomic oxygen loss in mass spectrometer ion sources can be reduced to an understanding of the possible surface interactions between oxygen atoms and the metal surface of the ion source. Results are presented for an experimental study in which an atomic oxygen beam apparatus and a mass spectrometer were used to measure the oxygen atom reflection, recombination, general surface reaction, and occlusion probabilities on six different engineering surfaces as a function of atomic oxygen exposure. The materials studied are gold, Nichrome V, aluminum, titanium, silver, and platinum. The variation in measured reflection probability seems to occur with metals that form oxides, Nichrome V being stable in terms of reflection stability. Recombination is observed an all surfaces except aluminum and platinum. Variation in the complete set of measurements in a single experiment is the result of varying surface conditions.

Sjolander, G. W.↗
NASA NTRS1986

Boundary layer transition

The boundary layer stability, its active control by sound and surface heating and the effect of curvature are studied numerically and experimentally for subsonic flow. In addition, the experimental and flight test data are correlated using the stability theory for supersonic Mach numbers. Active transition fixing and feedback control of boundary layer by sound interactions are experimentally investigated at low speed over an airfoil. Numerical simulation of active control by surface heating and cooling in air shows that by appropriate phase adjustment a reduction in the level of perturbation can be obtained. This simulation is based on the solution of two-dimensional compressible Navier-Stokes equations for a flat plate. Goertler vortices are studied experimentally on an airfoil in the Low Turbulence Pressure Tunnel (LTPT). The flow pattern was visualized using the sublimating chemical technique and data were obtained using a three component laser velocimeter. The effect of curvature on swept leading-edge stability on a cylinder was numerically studied. The results suggest that transition is dominated by traveling disturbance waves and that the waves with the greatest total amplification has an amplitude ratio of e sup 11. Experimental data from the quiet supersonic tunnel and flight tests are analyzed using linear compressible stability theory.

Maestrello, L.↗
DOE OSTI2022

How the Bioinspired Fe 2 Mo 6 S 8 Chevrel Breaks Electrocatalytic Nitrogen Reduction Scaling Relations

The nitrogen reduction reaction (NRR) is a renewable alternative to the energy- and CO 2 -intensive Haber–Bosch NH 3 synthesis process but is severely limited by the low activity and selectivity of studied electrocatalysts. The Chevrel phase Fe 2 Mo 6 S 8 has a surface Fe–S–Mo coordination environment that mimics the nitrogenase FeMo-cofactor and was recently shown to provide state-of-the-art activity and selectivity for NRR. Here, we elucidate the previously unknown NRR mechanism on Fe 2 Mo 6 S 8 via grand-canonical density functional theory (GC-DFT) that realistically models solvated and biased surfaces. Fe sites of Fe 2 Mo 6 S 8 selectively stabilize the key *NNH intermediate via a narrow band of free-atom-like surface d-states that selectively hybridize with p-states of *NNH, which results in Fe sites breaking NRR scaling relationships. These sharp d-states arise from an Fe–S bond dissociation during N 2 adsorption that mimics the mechanism of the nitrogenase FeMo-cofactor. Furthermore, we developed a new GC-DFT-based approach for calculating transition states as a function of bias (GC-NEB) and applied it to produce a microkinetic model for NRR at Fe 2 Mo 6 S 8 that predicts high activity and selectivity, in close agreement with experiments. Furthermore, our results suggest new design principles that may identify effective NRR electrocatalysts that minimize the barriers for *N 2 protonation and *NH 3 desorption and that may be broadly applied to the rational discovery of stable, multinary electrocatalysts for other reactions where narrow bands of surface d-states can be tuned to selectively stabilize key reaction intermediates and guide selectivity toward a target product. Furthermore, our results highlight the importance of using GC-DFT and GC-NEB to accurately model electrocatalytic reactions.

10 SYNTHETIC FUELS↗
DOE OSTI2025

Stabilization of magnetic bubbles in [Ni/Co]n multilayers on an oxygen-reconstructed Nb(110) surface via an ultra-thin Cu interlayer

Magnetic thin films hosting topological spin textures, such as magnetic skyrmions, hold high potential for breakthroughs in the field of spintronics, due to good scalability and energy efficiency. Novel computational architectures such as memory-in-logic devices rely on material platforms able to host those topological spin textures. Furthermore, recently proposed designs of novel quantum information technologies are based on heterostructures where topological spin textures are in direct proximity to a superconducting layer. Here, we demonstrate the stabilization of out-of-plane magnetic bubbles in highly ordered [Ni/Co]n multilayers on a Nb(110) single crystal. This is achieved without the need for the removal of the well-known Nb(110)-oxide surface reconstruction, due to the introduction of an ultra-thin Cu interlayer in between the Nb substrate and the magnetic multilayer. The Cu interlayer generates a well-ordered hexagonal surface, which is key for the epitaxial growth of the [Ni/Co]n multilayers hosting the desired out-of-plane anisotropy. The magnetic ground state of the prepared material stacks is directly imaged via spin-polarized low-energy electron microscopy, revealing the presence of magnetic bubble domains with lateral sizes as small as 450 nm.

Dibajeh, Ahmad↗
NASA NTRS2002

Factors Affecting the Latitudinal Location of the Intertropical Convergence Zone in a GCM

The dominant role of the latitudinal peak of the sea surface temperature (SST) in determining the latitudinal location of the intertropical convergence zone (ITCZ) is well-known. However, the roles of the other factors are less well-known and are the topic of this study. These other factors include the inertial stability, the interaction between convection and surface fluxes and the interaction between convection and radiation. Since these interactions involve convection, in a model they involve the cumulus parameterization scheme. These factors are studied with a general circulation model with uniform SST and solar angle. Under the aforementioned model settings, the latitudinal location of the ITCZ is the latitude where the balance of two types of attraction on the ITCZ, both due to earth's rotation, exists. Directly related to the Coriolis parameter, the first type pulls the ITCZ toward the equator and is not sensitive to model design changes. Related to the convective circulation, the second type pulls the ITCZ poleward and is sensitive to model design changes. Due to the shape and the magnitude of the attractors, the balance of the two types of attractions is reached either at the equator or more than 10 degrees away from the equator. The former case results in a single ITCZ over the equator and the latter case a double ITCZ straddling the equator.

Chao, Winston C.↗
DOE OSTI2021

Vapor-phase self-assembly for generating thermally stable single-atom catalysts

Preparation of thermally stable metal single-atom catalysts (SACs) is a challenge in heterogeneous catalysis, especially on conventional supports that provide a weak metal-support interaction. In this work, we report that a modified support MgAl 2 O 4 can stabilize Pt single atoms by a mechanism of vapor-phase self-assembly in a high-temperature treatment (800°C, air). The experimental results on the formation mechanism and the structure are validated by DFT and ab initio molecular dynamics simulations. We infer that stable triangular K 3 O 3 structures help stabilize Pt single atoms at high temperatures in oxidizing conditions, exhibiting excellent reactivity for methane oxidation. The obtained Pt/K/MgAl 2 O 4 SAC presents excellent stability in methane oxidation after steam treatment at elevated temperatures, whereas the Pt/MgAl 2 O 4 nanocatalyst suffers from rapid deactivation due to Pt nanoparticle growth. Now this work paves the way for preparing thermally stable and highly active SACs using conventional high-surface-area supports, despite the weak metal-support interaction.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1995

Vestibular plasticity following orbital spaceflight: recovery from postflight postural instability

Results of previous studies suggested that the vestibular mediated postural instability observed in astronauts upon return to earth from orbital spaceflight may be exacerbated by an increased weighting of visual inputs for spatial orientation and control of movement. This study was performed to better understand the roles of visual and somatosensory contributions to recovery of normal sensori-motor postural control in returning astronauts. Preflight and postflight, 23 astronaut volunteers were presented randomly with three trials of six sensory organization test (SOT) conditions in the EquiTest system test battery. Sagittal plane center-of-gravity (COG) excursions computed from ground reaction forces were significantly higher on landing day than preflight for those test conditions presenting sway-referenced visual and/or somatosensory orientation cues. The ratio of summed peak-to-peak COG sway amplitudes on the two sway-referenced vision tests (SOTs 3 + 6) compared to the two eyes closed tests (SOTs 2 + 5) was increased on landing day, indicating an increased reliance on visual orientation cues for postural control. The ratio of peak-to-peak COG excursions on sway-referenced surfaces (SOTs 4, 5 & 6) to an earth fixed support surfaces (SOTs 1, 2 & 3) increased even more after landing suggesting primary reliance on somatosensory orientation cues for recovery of postflight postural stability. Readaptation to sway-referenced support surfaces took longer than readaptation to sway-referenced vision. The increased reliance on visual and somatosensory inputs disappeared in all astronauts 4-8 days following return to earth.

NASA Program Space Physiology and Countermeasures↗
DOE OSTI2020

A short tutorial on uranium electropolishing

Electropolishing of uranium metal is best performed in a simple aqueous solution of phosphoric acid. With +3 to +10 volts of bias potential applied to the uranium (anode) relative to a Pt counter electrode (cathode). As summarized in a Pourbaix stability diagram, these conditions yield dissolution of the U surface as UOX containing solution species without the possibility of injecting hydrogen into the metal surface. However, if the uranium is biased to negative potentials relative to the counter electrode then water may be reduced at the U surface and hydrogen gas becomes available for reaction with the U metal. The stability diagram may be modified if there are other complexing species present in solution and the potential scale shifted if a different metal is use for the counter electrode. These effects can be understood in a known way as long as solution and electrochemical cell conditions are recognized.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
NASA NTRS1947

Fundamental Aerodynamic Investigations for Development of Arrow-Stabilized Projectiles

The numerous patent applications on arrow-stabilized projectiles indicate that the idea of projectiles without spin is not new, but has appeared in various proposals throughout the last decades. As far as projectiles for subsonic speeds are concerned, suitable shapes have been developed for sometime, for example, numerous grenades. Most of the patent applications, though, are not practicable particularly for projectiles with supersonic speed. This is because the inventor usually does not have any knowledge of aerodynamic flow around the projectile nor any particular understanding of the practical solution. The lack of wind tunnels for the development of projectiles made it necessary to use firing tests for development. These are obviously extremely tedious or expensive and lead almost always to failures. The often expressed opinion that arrow-stabilized projectiles cannot fly supersonically can be traced to this condition. That this is not the case has been shown for the first time by Roechling on long projectiles with foldable fins. Since no aerodynamic investigations were made for the development of these projectiles, only tedious series of firing tests with systematic variation of the fins could lead to satisfactory results. These particular projectiles though have a disadvantage which lies in the nature cf foldable fins. They occasionally do not open uniformly in flight, thus causing unsymmetry in flow and greater scatter. The junctions of fins and body are very bad aerodynamically and increase the drag. It must be possible to develop high-performance arrow-stabilized projectiles based on the aerodynamic research conducted during the last few years at Peenemuende and new construction ideas. Thus the final shape, ready for operational use, could be developed in the wind tunnel without loss of expensive time in firing tests. The principle of arrow-stabilized performance has been applied to a large number of caliburs which were stabilized by various means Most promising was the development of a subcaliber wing-stabilized projectile with driving disc (Treibspiegel) where rigid control surfaces extend beyond the caliber of the projectile into the free stream. The stabilized projectiles of full-caliber, wing-stabilized projectiles with fins within the caliber is considerably more difficult. A completely satisfactory solution for the latter has not been found yet.

Kurzweg, Hermann↗
DOE OSTI2021

Electronic, magnetic, and optical properties of bulk and (1 1 1)-surfaces of CoMnZnSi quaternary Heusler alloy

The structural stability, electronic structure, and magnetic and optical properties of CoMnZnSi quaternary Heusler alloy (QHA) in the bulk and (1 1 1)-slab forms have been investigated by performing density functional theory (DFT) calculations. Formation and cohesive energies, and elastic constants confirm that the bulk CoMnZnSi is chemically and mechanically stable at an equilibrium lattice parameter of 5.81 Å. The calculated elastic constants also indicate that this QHA has ductile and anisotropic features. We investigate and discuss the bonding behavior from charge density distribution and density of states. Furthermore, these calculations in the bulk phase show a perfect half-metallic behavior with an integer value of magnetic moment (4 µB) and a large spin-flip gap of 0.39 eV. On the other hand, the Co (1 1 1) and Si (1 1 1)-slabs exhibit semi-metallic nature at the Fermi level. These results indicate that CoMnZnSi qualifies for spintronic applications. In the bulk and Co (1 1 1), Zn (1 1 1), and Si (1 1 1)-surfaces, the magnetic moment of Co and Mn atoms are ferromagnetically aligned, while the ferrimagnetic alignment of Co with Mn has been found in Mn (1 1 1)-surface. We also investigate the optimal conductivity, dielectric functions, reflectivity, absorption coefficient, refractive index, and loss function to understand the underlying optical properties of this alloy.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Evolution of Active Sites in Pt-Based Nanoalloy Catalysts for the Oxidation of Carbonaceous Species by Combined in Situ Infrared Spectroscopy and Total X-ray Scattering

Here, we present results from combined in situ infrared spectroscopy and total X-ray scattering studies on the evolution of catalytically active sites in exemplary binary and ternary Pt-based nanoalloys during a sequence of CO oxidation–reactivation–CO oxidation reactions. We find that when within a particular compositional range, the fresh nanoalloys may exhibit high catalytic activity for low-temperature CO oxidation. Using surface-specific atomic pair distribution functions (PDFs) extracted from the in situ total X-ray scattering data, we find that, regardless of their chemical composition and initial catalytic activity, the fresh nanoalloys suffer a significant surface structural disorder during CO oxidation. Upon reactivation in oxygen atmosphere, the surface of used nanoalloy catalysts both partially oxidizes and orders. Remarkably, it largely retains its structural state when the nanoalloys are reused as CO oxidation catalysts. The seemingly inverse structural changes of studied nanoalloy catalysts occurring under CO oxidation and reactivation conditions affect the active sites on their surface significantly. In particular, through different mechanisms, both appear to reduce the CO binding strength to the nanoalloy’s surface and thus increase the catalytic stability of the nanoalloys. The findings provide clues for further optimization of nanoalloy catalysts for the oxidation of carbonaceous species through optimizing their composition, activation, and reactivation. Besides, the findings demonstrate the usefulness of combined in situ infrared spectroscopy and total X-ray scattering coupled to surface-specific atomic PDF analysis to the ongoing effort to produce advanced catalysts for environmentally and technologically important applications.

36 MATERIALS SCIENCE↗
NASA NTRS1947

Wind-Tunnel Tests of the 1/25-Scale Powered Model of the Martin JRM-1 Airplane. IV - Tests with Ground Board and with Modified Wing and Hull - TED No. NACA 232: Tests with Ground Board and with Modified Wing and Hull, TED No. NACA 232 - Part 4

Wind-tunnel tests were made of a 1/25 scale model of the Martin JRM-1 airplane to determine: (1) The longitudinal stability and control characteristics of the JRM-1 model near the water and lateral and directional stability characteristics with power while moving on the surface of the water, the latter being useful for the design of tip floats; (2) The stability and stalling characteristics of the wing with a modified airfoil contour; (3) Stability characteristics of a hull of larger design gross weight; The test results indicated that the elevator was powerful enough to trim the original model in a landing configuration at any lift coefficient within the specified range of centers of gravity. The ground-board tests for evaluating the aerodynamic forces and moments on an airplane in a simulated cross wind indicate a high dihedral effect in the presence of the ground board and, consequently, during low-speed taxying and take-off, large overturning moments would result which would have to be overcome by the tip floats.

Lockwood, Vernard E.↗
NASA NTRS2006

Three Years of Atmospheric Infrared Sounder Radiometric Calibration Validation using Sea Surface Temperatures

This paper evaluates the absolute accuracy and stability of the radiometric calibration of the Atmospheric Infrared Sounder (AIRS) by analyzing the difference between the brightness temperatures measured at 2616 cm(exp -1) and those calculated at the top of the atmosphere (TOA), using the Real-Time Global Sea Surface Temperature (RTGSST) for cloud-free night tropical oceans between +/- 30 degrees latitude. The TOA correction is based on radiative transfer. The analysis of the first 3 years of AIRS radiances verifies the absolute calibration at 2616 cm(exp -1) to better than 200 mK, with better than 16 mK/yr stability. The AIRS radiometric calibration uses an internal full aperture wedge blackbody with the National Institute of Standards and Technology (NIST) traceable prelaunch calibration coefficients. The calibration coefficients have been unchanged since launch. The analysis uses very tight cloud filtering, which selects about 7000 cloud-free tropical ocean spectra per day, about 0.5% of the data. The absolute accuracy and stability of the radiometry demonstrated at 2616 cm(sup -1) are direct consequences of the implementation of AIRS as a thermally controlled, cooled grating-array spectrometer and meticulous attention to details. Comparable radiometric performance is inferred from the AIRS design for all 2378 channels. AIRS performance sets the benchmark for what can be achieved with a state-of-the-art hyperspectral radiometer from polar orbit and what is expected from future hyperspectral sounders. AIRS was launched into a 705 km altitude polar orbit on NASA's Earth Observation System (EOS) Aqua spacecraft on 4 May 2002. AIRS covers the 3.7-15.4 micron region of the thermal infrared spectrum with a spectral resolution of nu/Delta nu = 1200 and has returned 3.7 million spectra of the upwelling radiance each day since the start of routine data gathering in September 2002.

climate↗
DOE OSTI2021

A Statistical Evaluation of WRF-LES Trace Gas Dispersion Using Project Prairie Grass Measurements

In recent years, new measurement systems have been deployed to monitor and quantify methane emissions from the natural gas sector. Large-eddy simulation (LES) has complemented measurement campaigns by serving as a controlled environment in which to study plume dynamics and sampling strategies. However, with few comparisons with controlled-release experiments, the accuracy of LES for modeling natural gas emissions is poorly characterized. In this paper, we evaluate LES from the Weather Research and Forecasting (WRF) Model against Project Prairie Grass campaign measurements and surface layer similarity theory. Using WRF-LES, we simulate continuous emissions from 30 near-surface trace gas sources in two stability regimes: strong convection and weak convection. We examine the impact of grid resolutions ranging from 6.25 to 52 m in the horizontal dimension on model results. We evaluate performance in a statistical framework, calculating fractional bias and conducting Welch’s t tests. WRF-LES accurately simulates observed surface concentrations at 100 m and beyond under strong convection; simulated concentrations pass t tests in this region irrespective of grid resolution. However, in weakly convective conditions with strong winds, WRF-LES substantially overpredicts concentrations—the magnitude of fractional bias often exceeds 30%, and all but one t test fails. The good performance of WRF-LES under strong convection correlates with agreement with local free convection theory and a minimal amount of parameterized turbulent kinetic energy. The poor performance under weak convection corresponds to misalignment with Monin–Obukhov similarity theory and a significant amount of parameterized turbulent kinetic energy.

17 WIND ENERGY↗
DOE OSTI2017

Na-Montmorillonite Edge Structure and Surface Complexes: An Atomistic Perspective

The edges of montmorillonite (MMT) react strongly with metals and organic matter, but the atomic structure of the edge and its surface complexes are not unambiguous since the experimental isolation of the edge is challenging. In this study, we introduce an atomistic model of a Na MMT edge that is suitable for classical molecular dynamics (MD) simulations, in particular for the B edge, a representative edge surface of 2:1 phyllosilicates. Our model possesses the surface groups identified through density functional theory (DFT) geometry optimizations performed with variation in the structural charge deficit and Mg substitution sites. The edge structure of the classical MD simulations agreed well with previous DFT-based MD simulation results. Our MD simulations revealed an extensive H-bond network stabilizing the Na-MMT edge surface, which required an extensive simulation trajectory. Some Na counter ions formed inner-sphere complexes at two edge sites. The stronger edge site coincided with the exposed vacancy in the dioctahedral sheet; a weaker site was associated with the cleaved hexagonal cavity of the tetrahedral sheet. The six-coordinate Na complexes were not directly associated with the Mg edge site. Our simulations have demonstrated the heterogeneous surface structures, the distribution of edge surface groups, and the reactivity of the MMT edge.

58 GEOSCIENCES↗