Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Liquid solid interfaces”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Final Technical Report for DE-SC0021049: Manipulating interfacial reactivity with atomically layered heterostructures

This final technical report summarizes the work accomplished in this DOE Early Career Research Program project that has established moiré superlattice materials and two-dimensional (2D) heterostructures as a highly tunable platform for controlling heterogeneous charge transfer (ET) kinetics at solid-liquid interfaces. By precisely engineering van der Waals heterostructures of atomically thin 2D materials, particularly bilayer and trilayer graphene with controlled twist angles, this project demonstrated systematic control of interfacial charge transfer rates spanning three orders of magnitude. This research addresses fundamental questions about how electronic structure, charge localization, and atomic layer-dependent properties govern charge transfer at electrochemical interfaces, with broad implications for energy conversion, electrocatalysis, and next-generation electrochemical devices.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Designing carbon wetting layers with inorganic salts as additives for enhanced adhesion and molten sodium wettability in electrochemical application

ß” alumina solid-state electrolyte (BASE) is considered one of the most promising materials for sodium batteries due to its high ionic conductivity and exceptional stability with molten sodium. However, at lower operating temperatures, the poor physical contact at the liquid-solid interface between molten sodium and BASE results in high interfacial resistance. In this study, we explore the optimization of carbon wetting layers with various salt additives to enhance adhesion and improve molten sodium wettability. Using a spray-coating method, carbon layers incorporating sodium hexametaphosphate (Na6(PO3)6, SHMP) as an additive demonstrate superior adhesion, mechanical stability, and enhanced wettability with molten sodium. Scanning Electron Microscopy (SEM) analysis reveals that SHMP preserves the porous carbon network required for efficient sodium transport while minimizing structural defects such as cracks and voids. Electrochemical testing using Na symmetric cells confirms that SHMP-modified layers significantly reduce interfacial resistance and overpotential, outperforming No salts counterparts. These findings highlight that the strategic incorporation of suitable additives, such as SHMP, in the design of microscale carbon layers could significantly enhance the performance of molten sodium batteries at lower operating temperatures.

Han, Henry H.↗
DOE OSTI2025

In-situ/operando study of Cu-based nanocatalysts for CO 2 electroreduction using electrochemical liquid cell TEM

The structure of a nanocatalyst during electrocatalytic reactions often deviates from its pristine structure due to intrinsic properties, or physical and chemical adsorption at the catalytic surfaces. Taking Cu-based catalysts for CO 2 electroreduction reactions (CO 2 RR) as an example, they often experience segregation, leaching, and alloying during reactions. With the recent breakthrough development of high-resolution polymer electrochemical liquid cells, in-situ electrochemical liquid cell transmission electron microscopy (EC-TEM) alongside other advanced microscopy techniques, has become a powerful platform for revealing electrocatalysts restructuring at the atomic level. Considering the complex reactions involving electrified solid-liquid interfaces and catalyst structural evolution with intermediates, systematic studies with multimodal approaches are crucial. In this article, we demonstrate a research protocol for the study of electrocatalysts structural evolution during reactions using the in-situ EC-TEM platform. Using Cu and CuAg nanowire catalysts for CO 2 RR as model systems, we describe the experimental procedures and findings. We highlight the platform’s crucial role in elucidating atomic-scale pathways of nanocatalyst restructuring and identifying catalytic active sites, as well as avoiding potential artifacts to ensure unbiased conclusions. Using the multimodal characterization toolbox, we provide the opportunity to correlate the structure of a working catalyst with its performance. Finally, we discuss advancements as well as the remaining gap in elucidating the structural-performance relationship of working catalysts. We expect this article will assist in establishing guidelines for future investigations of complex electrochemical reactions, such as CO 2 RR and other catalytic processes, using the in-situ EC-TEM platform.

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

Liquid Core Detection and Strand Condition Monitoring in a Continuous Caster Using Optical Fiber

Real-time monitoring of the liquid core position during the continuous casting of steel has been demonstrated using low-cost distributed optical-fiber-based strain sensors. These sensors were installed on the containment roll support structures in the segments of a production continuous caster to detect the position of the solid–liquid interface and monitor the strand condition during the continuous casting. Distributed Fiber Bragg Grating sensors (FBGs) were used in this work to monitor strain at six roll positions in the caster. The sensor performance was first validated by comparing optical strain measurements with conventional strain gauge measurements in the lab. Next, optical strain measurements were performed on an isolated caster segment in a segment maintenance facility using hydraulic jacks to simulate the presence of a liquid core under the roll. Finally, the sensors were evaluated during caster operation. The sensors successfully detected the load increase associated with the presence of a liquid core under each instrumented roll location. Incidents of bulging and roll eccentricity were also detected using frequency analysis of the optical strain signal. The liquid core position measurements were compared using predictions from computer models (digital twins) in use at the production site. The measurements were in good agreement with the model predictions, with a few exceptions. Under certain transient caster operating conditions, such as spraying practice changes and SEN exchanges, the model predictions deviated slightly from the liquid core position determined from strain measurements.

47 OTHER INSTRUMENTATION↗
DOE OSTI2020

Efficient Mesofluidic Separation of Large Particles in nuclear slurries - 20408

Efficient and effective particle separation is essential to cleanup of many nuclear wastes. For example, particle separation may be used to accelerate settle-decant bottlenecks that limit the throughput necessary to achieve the U.S. Department of Energy waste cleanup mission in a timely manner. Particle separation may be used to protect downstream processing equipment from pulses of solids that may be encountered near (within a few feet of) liquid-solid interfaces during waste transfer operations, thereby preventing work stoppages. Effective particle separation may permit efficient sludge washing. These are only a few of the ways in which particle separation is important. Yet, efficient and effective particle separation in nuclear processing environments remains challenging. For example, settle-decant operations permit larger and heavier particles to settle, leaving smaller, lighter particles suspended. However, settle-decant operations are slow, convection currents may resuspend solids, hindered settling of multicomponent slurries remains incompletely understood, and predictive models for settling in graduated cylinders fail to match limited observations of settling in large waste tanks. Additionally, pumping operations can cause turbulent resuspension of particles when the pump intake is close to the settled solids layer. Filtration techniques that use membranes or partially permeable barriers retain larger particles, permitting smaller particles to permeate. However, filters, including dead-end filtration, are prone to clogging and caking, operate at elevated pressures due to minimal void volume, and increase in pressure during operations, which require more control systems. Therefore, the need for high throughput particle separation techniques that operate with modest pressure drops persists. A novel mesofluidic separator presents the opportunity to effectively and efficiently accelerate the waste cleanup mission. The separator separates large particles from process streams across a broad range of particle sizes and has no moving parts or media to replace, regenerate, or clean. This separator design has an unusually large void volume, permitting operation at much higher flow rates (and lower pressures) than traditional filtration (e.g., dead-end filtration). Industrial-scale flow rates have been demonstrated. In performance testing, mesofluidic separators operate at might flow rates (>90 gpm (0.0006 m{sup 3}/s) in piping 3 inches (0.08 m) in diameter; Re>10{sup 5}), with modest pressure drop (∼25 psi (170 KPa) in testing). In complex, aggregating waste simulants, the separator loses <25% of flow projected over months without back pulsing or chemical cleaning. Mesofluidic separation presents infrastructure advantages, reduces risk, and provides mission impact. Infrastructure advantages include implementation within existing transfer systems, plug and play without facility modifications to safety systems, and operation at low pressures. The separator is flexible in location, may be positioned within or outside of waste tanks, and may replace or augment dead-end and cross-flow filters. Mesofluidic separation reduces risk by minimizing waste (no media to replace or regenerate and no cleaning chemicals to handle and dispose) and by minimizing filter change-out consistent with as low as reasonably achievable (ALARA) exposure to workers. The potential mission impact of these separators is substantial. The separator has potential to support, simplify, and accelerate in-farm transfers and waste feed delivery. Furthermore, the separator can free up settling and holding tanks in the Direct Feed Low-Activity Waste (DFLAW) mission by reducing or eliminating post-transfer settling and wait times with the potential to free up 1 Mgal (4000 m{sup 3}) of double-shell tank space. This paper quantitatively compares mesofluidic separation to dead-end filtration, discusses scale-up results, and considers the separator's potential to efficiently and effectively reduce the long-term environmental ability of particulate-rich nuclear wastes. (authors)

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

Liquid Interfacial Electron Microscopy Identifies Nanogalvanic Corrosion in Pearlitic Steel

The nanoscale mechanisms of localized corrosion in low carbon steels have remained elusive due to the complexity of studying the degradative material behavior at nanoscale solid-liquid interfaces. We identified various steps in the nanogalvanic corrosion processes using in-situ liquid-cell scanning transmission electron microscopy (STEM) using a microfluidic holder by Hummingbird Scientific. Initial work, performed at low magnification, identified the initiation point on a 1018 low-carbon steel surface. This initiation point was determined to be a triple junction of two ferrite grains bridging a cementite grain in contact with a baseline electrolyte of 6 uM CO2 dissolved in a buffered (2.78 uM Na2SO4) aqueous solution, pH 6.1. The pre-etched low-carbon steel surface was prepared using focused ion beam lift-out procedures to extract a cross-section of the low-carbon steel surface, which then was thinned to about 150 nm and transferred to a SiN membrane microfluidic window. The transfer was made using a lift-out needle to attach the low-carbon steel lamella to the corner of the SiN window, and then Pt/C deposition held the lamella in contact with the window while it was released from the lift out needle. To identify the triple point on the low carbon steel lamella, prior to attachment on the SiN window, the sample was characterized for compositional variations with energy dispersive x-ray spectroscopy mapping, grain orientation and phase mapping with precession electron diffraction, and thickness mapping with energy filtered transmission electron microscopy. This pre-characterization prior to the in-situ experiment provided a map of the multiphase and multigrain structure, where the in-situ liquid cell imaging provided a clear understanding of the initiation point on the sample. These data were cross-correlated to paint a holistic picture of the triple junction site, enabling low electron-fluence in-situ snapshot imaging to avoid dominating the native corrosion reactions with effects from the incident electron beam. This initial result identified that localized, nanogalvanic corrosion at the phase interface was the dominant corrosion process in the low-carbon steel, so we next targeted the observation of an array of these nanogalvanic features phase boundaries in a pearlite grain. Near-surface ferrite/cementite phase interfaces that typify pearlitic low-carbon steel were extracted, pre-characterized, and imaged for the in-situ corrosion processes. The sample was a cross-section from a pearlite grain, with alternating ferrite and cementite grains that extended microns down from the pre-etched low-carbon steel pipe surface. After contact with a buffered aqueous solution, the phase boundaries between the ferrite and cementite began to dissolve, with observable material loss and thickness changes in the dark-field and bright-field STEM images. Within minutes, the corrosion front proceeded deeper into the material, claiming a thin layer of ferrite around all exposed phase boundaries before progressing laterally into the ferrite matrix, converting the ferrite to corrosion product normal to each buried cementite grain. Formation of the corrosion product causes a volumetric expansion, creating a lateral wedging force that mechanically ejects the cementite grains from their grooves and leaves behind percolation channels into the steel substructure. Rapid and deleterious, this nanogalvanic corrosion pathway represents an important target for understanding and preventing run-away degradation in this common building material. Observation of this corrosion mechanism was enabled by the combination of pre-characterization using standard structural, grain, and compositional analysis in the TEM, which provides maps for understanding the reaction propagation captured in low-dose, in-situ, liquid-cell STEM.

corrosion↗
DOE OSTI2024

Characterization of compositional variation and solidification condition on rapid solidification behavior of beta stabilized titanium

Rapid solidification (RS) is a field of study that has increased in importance due to the rise of additive manufacturing (AM) techniques. The breakdown of equilibrium at the solid-liquid interface results in an inability to use traditional solidification models that have been applied to great success in casting. As such, there is a current need for a high-throughput testing method to screen alloys at rapid solidification conditions without the significant cost for powder atomization. In this work, twin hammer splat quenching (SQ) is utilized and theoretically developed as a technique uniquely placed to test alloy systems for AM applications. Using analytical thermal spray models and numerical heat transfer simulations, the cooling rate and solidification conditions of splat quenching are compared across several parameters to evaluate the key parameters for generating cooling rates relevant to AM process. Once a framework for the cooling rates and relevancy of SQ are developed, SQ was used to process alloy compositions in the Ti5553 alloy space (AMS 7026) to study the formation of detrimental secondary phases occurring upon rapid quenching from the melt. A martensitic phase is found to occur at the lower end of the standard, but the composition listed in the standard were found to be safe. Finally, the rapid solidification behavior of Ti-5553 is evaluated by comparing splat quenching, laser welding, and plasma atomization that offer a wide window of solidification conditions. The segregation and solidification structure are analyzed and compared to models to provide a deeper understanding of the solidification behavior of beta stabilized titanium alloys.

Williamson, Christopher J.↗
DOE OSTI2022

Manipulating Interfacial Dissolution–Redeposition Dynamics to Resynthesize Electrode Surface Chemistry

The solid-liquid electrochemical interface offers a two-dimensional environment for geometrically confined interfacial reactions to tailor electrode surface chemistry under operating conditions. In this study, we demonstrate that the dissolution and redeposition kinetics of transition metal cations, a ubiquitous phenomenon at the electrochemical interface, can be manipulated to regulate the chemical composition and crystal structure of the electrode surface as well as the overall electrochemical performance. Foreign cations (e.g., Ti 4+ ), either added as electrolyte additives or dissolved from surface coatings, can rapidly participate in the electrode dissolution-redeposition process, and facilitate the establishment of the dissolution-redeposition equilibrium. Our work expands the control over the electrochemical reactions at the solid-liquid interface and provides new insights into interfacial studies in electrochemistry and surface science.

25 ENERGY STORAGE↗
DOE OSTI2022

Drag reduction and the Vogel exponent of a flexible beam in transient shear flows

Interactions between a flexible beam and a fluid in a channel are of great relevance to biological hairy surfaces, aquatic vegetation, marine life (e.g., fish gills), and many industrial systems alike. While steady state response of a beam to such flows is fairly well-explored, their behavior in the transient regime is not fully understood. Here, a series of numerical simulations are performed to study the laminar Couette flow of an incompressible viscous fluid past an elastic beam in a two-dimensional channel. The flexible beam is perpendicular to the direction of flow, and its base is fixed to the stationary bottom of the channel. We measure the evolution of the Vogel exponent, drag reduction, and reconfiguration number during the transient and steady-state response of the fluid–structure system for different geometrical and physical properties. Our benchmark shows a good agreement between numerical and experimental observations. Our results show that the system's steady-state response at different bulk-fluid velocities can be reproduced by investigating the shear flow response during the transient regime. We define a new variable that characterizes the evolution of the local velocity profile in the proximity of the free end of the beam and use that to characterize the transient-regime response. The analysis yields insight into the competing effects of elasticity of the beam and non-linear flow response.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

In situ x-ray imaging to understand subsurface behavior during continuous wave laser drilling

A limited understanding regarding the underlying dynamics and mechanisms of material removal during continuous wave laser drilling has presented significant challenges in achieving precision and process control. Here, to address this, we employed high-fidelity, in situ synchrotron x-ray imaging to reveal previously unknown material behaviors during continuous wave laser drilling with power modulation. Our findings highlight that high-aspect ratio drill holes are achieved when the laser modulation frequency falls within the range of 8–12 kHz, provided that the laser average power and modulation amplitude levels meet the specified limits. Under these conditions, we identified a material removal mechanism driven by incremental accumulation of recoil pressure that gradually pushes material upward from deep within the substrate to the surface. This mechanism manifested as a low-frequency fluctuation in the vapor depression depth, resulting in periodic instances of material ejection. Furthermore, our study underscores that rapid expansion of the melt pool and the widening of the drill hole opening can impede effective material removal by redirecting energy from material ejection to increasing the melt pool size. This investigation contributes essential insights into the subsurface dynamics involved in the drilling of high-aspect ratio holes, furthering our fundamental understanding of this intricate process.

47 OTHER INSTRUMENTATION↗
DOE OSTI2022

Alcove formation in dissolving cliffs driven by density inversion instability

We demonstrate conditions that give rise to cave-like features commonly found in dissolving cliffsides with a minimal two-phase physical model. Alcoves that are wider at the top and tapered at the bottom, with sharp-edged ceilings and sloping floors, are shown to develop on vertical solid surfaces dissolving in aqueous solvents. As evident from descending plumes, sufficiently large indentations evolve into alcoves as a result of the faster dissolution of the ceiling due to a solutal Rayleigh–Bénard density inversion instability. In contrast, defects of size below the boundary layer thickness set by the critical Rayleigh number smooth out, leading to stable planar interfaces. Furthermore, the ceiling recession rate and the alcove opening area evolution are shown to be given to first-order by the critical Rayleigh number. By tracking passive tracers in the fluid phase, we show that the alcoves are shaped by the detachment of the boundary layer flow and the appearance of a pinned vortex at the leading edge of the indentations. The attached boundary layer past the developing alcove is then found to lead to rounding of the other sides and the gradual sloping of the floor.

58 GEOSCIENCES↗
DOE OSTI2021

Simulations of the IR and Raman spectra of water confined in amorphous silica slit pores

Water in nano-scale confining environments is a key element in many biological, material, and geological systems. The structure and dynamics of the liquid can be dramatically modified under these conditions. Probing these changes can be challenging, but vibrational spectroscopy has emerged as a powerful tool for investigating their behavior. A critical, evolving component of this approach is a detailed understanding of the connection between spectroscopic features and molecular-level details. In this paper, this issue is addressed by using molecular dynamics simulations to simulate the linear infrared (IR) and Raman spectra for isotopically dilute HOD in D 2 O confined in hydroxylated amorphous silica slit pores. The effect of slit-pore width and hydroxyl density on the silica surface on the vibrational spectra is also investigated. The primary effect of confinement is a blueshift in the frequency of OH groups donating a hydrogen bond to the silica surface. This appears as a slight shift in the total (measurable) spectra but is clearly seen in the distance-based IR and Raman spectra. Analysis indicates that these changes upon confinement are associated with the weaker hydrogen-bond accepting properties of silica oxygens compared to water molecules.

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

Liquid–Solid Boundaries Dominate Activity of CO 2 Reduction on Gas-Diffusion Electrodes

Electrochemical CO 2 electrolysis to produce hydrocarbon fuels or material feedstocks offers a renewable alternative to fossilized carbon sources. Gas-diffusion electrodes (GDEs), composed of solid electrocatalysts on porous supports positioned near the interface of a conducting electrolyte and CO 2 gas, have been able to demonstrate the substantial current densities needed for future commercialization. These higher reaction rates have often been ascribed to the presence of a three-phase interface, where solid, liquid, and gas provide electrons, water, and CO 2 , respectively. Conversely, mechanistic work on electrochemical reactions implicates a fully two-phase reaction interface, where gas molecules reach the electrocatalyst’s surface by dissolution and diffusion through the electrolyte. Because the discrepancy between an atomistic three-phase versus two-phase reaction has substantial implications for the design of catalysts, gas-diffusion layers, and cell architectures, the nuances of nomenclatures and governing phenomena surrounding the three-phase-region require clarification. In this paper we outline the macro, micro, and atomistic phenomena occurring within a gas-diffusion electrode to provide a focused discussion on the architecture of the often-discussed three-phase region for CO 2 electrolysis. From this information, we comment on the outlook for the broader CO 2 electroreduction GDE cell architecture.

CO2 reduction↗
DOE OSTI2021

Molecular Properties and Chemical Transformations Near Interfaces

The properties of bulk water and aqueous solutions are known to change in the vicinity of an interface and/or in a confined environment, including the thermodynamics of ion selectivity at interfaces, transition states and pathways of chemical reactions, and nucleation events and phase growth. Here we describe joint progress in identifying unifying concepts about how air, liquid, and solid interfaces can alter molecular properties and chemical reactivity compared to bulk water and multicomponent solutions. Furthermore, we also discuss progress made in interfacial chemistry through advancements in new theory, molecular simulation, and experiments.

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

Validating first-principles molecular dynamics calculations of oxide/water interfaces with x-ray reflectivity data

Metal oxide/water interfaces play a crucial role in many electrochemical and photocatalytic processes, such as photoelectrochemical water splitting, the creation of fuel from sunlight, and electrochemical CO 2 reduction. First-principles electronic structure calculations can reveal unique insights into these processes, such as the role of the alignment of the oxide electronic energy levels with those of liquid water. An essential prerequisite for the success of such calculations is the ability to predict accurate structural models of these interfaces, which in turn requires careful experimental validation. Here we report a general, quantitative validation protocol for first-principles molecular dynamics simulations of oxide/aqueous interfaces. The approach makes direct comparisons of interfacial x-ray reflectivity (XR) signals from experimental measurements and those obtained from ab initio simulations with semilocal and van der Waals functionals. The protocol is demonstrated here for the case of the Al 2 O 3 (001)/water interface, one of the simplest oxide/water interfaces. We discuss the technical requirements needed for validation, including the choice of the density functional, the simulation cell size, and the optimal choice of the thermodynamic ensemble. Our results establish a general paradigm for the validation of structural models and interactions at solid/water interfaces derived from first-principles simulations. Furthermore, while there is qualitative agreement between the simulated structures and the experimental best-fit structure, direct comparisons of simulated and measured XR intensities show quantitative discrepancies that derive from both bulk regions (i.e., alumina and water) as well as the interfacial region, highlighting the need for accurate density functionals to properly describe interfacial interactions. Our results show that XR data are sensitive not only to the atomic structure (i.e., the atom locations) but also to the electron-density distributions in both the substrate and at the interface.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Lattice Boltzmann model for conjugate heat transfer across thin walls

A lattice Boltzmann (LB) model with an efficient and accurate interface treatment for conjugate heat transfer across a thin wall between two different media is developed. The proposed interface treatment avoids fine meshing and computation within the thin layer; instead, the energy balance within the thin layer and the conjugate conditions on each interface are utilized to construct explicit updating schemes for the microscopic distribution functions of the LB model at the interior lattice nodes of the two media next to the thin layer. The proposed interface schemes reduce to the standard interface scheme for conjugate conditions in the literature in the limit of zero thickness of the thin layer, and thus it can be considered a more general interface treatment. A simplified version of the interface treatment is also proposed when the heat flux variation along the tangential direction of the thin layer is negligible. Three representative numerical tests are conducted to verify the applicability and accuracy of the proposed interface schemes. The results demonstrate that the intrinsic second-order accuracy of the LB model is preserved with the proposed interface schemes for thin layers with constant tangential fluxes, while for general situations with varying tangential fluxes, first-order accuracy is obtained. Finally, this interface treatment within the LB framework is attractive in conjugate heat transfer modeling involving thin layers for its simplicity, accuracy, and significant reduction in computational resources.

42 ENGINEERING↗
DOE OSTI2020

Effect of Collective Dynamics and Anharmonicity on Entropy in Heterogenous Catalysis: Building the Case for Advanced Molecular Simulations

In this paper, we present a perspective on the computational determination of entropy and its effects and consequences on heterogeneous catalysis. Special attention is paid to the role of anharmonicity and the deviations from the standard harmonic oscillator approximations which can fail to provide a reliable assessment of entropy. To address these challenges, advanced methodologies are needed that can explicitly account for these thermodynamic drivers through the appropriate statistical sampling of reactive free energy surfaces. We discuss where anharmonicity should be expected, where it has been observed from a theoretical perspective, and the methods currently employed to address it. We concentrate on three types of systems where we have observed major, nonnegligible anharmonic effects: (1) supported nanoparticles, where the migration of metal atoms, complexes, and entire clusters exhibit anharmonic behavior in their dynamic motion; (2) porous solids where confinement effects distort potential energy surfaces and hinder molecular motions resulting in large entropic terms; and (3) solid/liquid interfaces, where interactions between solvent molecules and adsorbed species can result in large solvent organization free energy and unique reactivity.

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

Validated Reactive Force Field Quantifies MXene Interfacial Properties, Mechanics, and Thermal Transport

MXenes combine rich surface chemistry, mechanical strength, and high conductivity for a multitude of emerging applications. Predictive modeling supports accelerated materials designs and has been limited by the absence of validated and transferable force fields. Here, we introduce an interpretable, reactive INTERFACE force field (IFF and IFF-R) for Ti 3 C 2 T x MXenes that is trained based on chemical knowledge and achieves quantitative agreement with experiments across lattice parameters (<0.5%), density (<0.2%), liquid contact angles, Raman spectra, and the in-plane elastic modulus (∼320 GPa). The models cover surface terminations from hydroxyl (−OH) to fluorine (−F) groups and are extensible to other chemistries. We introduce pH-resolved surface chemistry and identify dopamine adsorption mechanisms at MXene–aqueous interfaces supported by QCM-D and UV–Vis experiments. The data reveal coplanar and perpendicular binding modes and concentration-dependent multilayer assembly. We predict previously inaccessible properties, including termination-dependent cleavage energies, interlayer shear moduli and dynamic shear failure, nanoindentation and brittle fracture, anisotropic in-plane and out-of-plane thermal conductivities, including the role of defects. Agreement with available experimental data is consistently close and exceeds DFT accuracy across the benchmark properties examined. The IFF/IFF-R model is compatible with CHARMM, AMBER, OPLS, and CVFF force fields for simulations of MXenes with diverse surface terminations, electrolyte interfaces, biointerfaces, and polymer composites without additional parameters. Parameter sets, 3D models, and analysis scripts are provided for community use. The validated, reactive, and transferable IFF framework facilitates predictive design of MXene-based films, membranes, sensing interfaces, and composites.

MXene↗