Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “interfacial forces”

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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Finite Element Modeling of the Fiber-Matrix Interface in Polymer Composites

Polymer composites are used in numerous industries due to their high specific strength and high specific stiffness. Composites have markedly different properties than both the reinforcement and the matrix. Of the several factors that govern the final properties of the composite, the interface is an important factor that influences the stress transfer between the fiber and matrix. The present study is an effort to characterize and model the fiber-matrix interface in polymer matrix composites. Finite element models were developed to study the interfacial behavior during pull-out of a single fiber in continuous fiber-reinforced polymer composites. A three-dimensional (3D) unit-cell cohesive damage model (CDM) for the fiber/matrix interface debonding was employed to investigate the effect of interface/sizing coverage on the fiber. Furthermore, a two-dimensional (2D) axisymmetric model was used to (a) analyze the sensitivity of interface stiffness, interface strength, friction coefficient, and fiber length via a parametric study; and (b) study the shear stress distribution across the fiber-interface-matrix zone. It was determined that the force required to debond a single fiber from the matrix is three times higher if there is adequate distribution of the sizing on the fiber. The parametric study indicated that cohesive strength was the most influential factor in debonding. Moreover, the stress distribution model showed the debonding mechanism of the interface. It was observed that the interface debonded first from the matrix and remained in contact with the fiber even when the fiber was completely pulled out.

36 MATERIALS SCIENCE↗
NASA NTRS2013

Effect of Residual Noncondensables on Pressurization and Pressure Control of a Zero-Boil-Off Tank in Microgravity

The Zero-Boil-Off Tank (ZBOT) Experiment is a small-scale experiment that uses a transparent ventless Dewar and a transparent simulant phase-change fluid to study sealed tank pressurization and pressure control with applications to on-surface and in-orbit storage of propellant cryogens. The experiment will be carried out under microgravity conditions aboard the International Space Station in the 2014 timeframe. This paper presents preliminary results from ZBOT's ground-based research that focuses on the effects of residual noncondensable gases in the ullage on both pressurization and pressure reduction trends in the sealed Dewar. Tank pressurization is accomplished through heating of the test cell wall in the wetted and un-wetted regions simultaneously or separately. Pressure control is established through mixing and destratification of the bulk liquid using a temperature controlled forced jet flow with different degrees of liquid jet subcooling. A Two-Dimensional axisymmetric two-phase CFD model for tank pressurization and pressure control is also presented. Numerical prediction of the model are compared to experimental 1g results to both validate the model and also indicate the effect of the noncondensable gas on evolution of pressure and temperature distributions in the ullage during pressurization and pressure control. Microgravity simulations case studies are also performed using the validated model to underscore and delineate the profound effect of the noncondensables on condensation rates and interfacial temperature distributions with serious implications for tank pressure control in reduced gravity.

Marangoni and Natural Convection↗
NASA NTRS1998

Oscillatory/Chaotic Thermocapillary Flow Induced by Radiant Heating

There is a continuing need to understand the fluid physics occurring under low gravity conditions in processes such as crystal growth, materials processing, and the movement of bubbles or droplets. The fluid flow in such situations is often caused by a gradient in interfacial tension. If a temperature gradient is created due to a heat source, the resulting flow is called thermocapillary flow, a special case of Marangoni Convection. In this study, an experimental investigation was conducted using silicone oil in cylindrical containers with a laser heat source at the free surface. It was desired to determine the conditions under which steady, axisymmetrical thermocapillary flow becomes unstable and oscillatory three-dimensional flow states develop. The critical Marangoni number for each observed oscillatory state was measured as a function of the container aspect ratio and the dynamic Bond number, a measure of buoyant force versus ii thermocapillary force. Various oscillatory modes were observed during three- dimensional convection, and chaotic flow was reached in one test condition. The critical Marangoni numbers are compared with those measured in previous studies, and the power spectra and phase trajectories of the instantaneous surface temperature distributions are used to characterize the routes of transitions to the chaotic flow state. Results show that only superharmonic modes appear in the routes to chaos while infinite number of subharmonic modes occur in flow transitions for pure Rayleigh convection.

DeWitt, Kenneth J.↗
DOE OSTI2021

Process Robustness of Friction Stir Dovetailing of AA7099 to Steel with In Situ AA6061 Interlayer Linking

Rolled homogeneous armor (RHA) steel was joined to aluminum alloy AA7099-T7451 with an interlayer of AA6061-T651 by single-pass friction stir dovetailing (FSD). The evolution of the AA7099/AA6061 and AA6061/steel interfacial microstructures and mechanical properties of the joints were examined. The AA6061 interlayer prevented Zn migration from the AA7099 to the steel during FSD and inhibited formation of brittle, Zn-rich, Fe-Al intermetallic compounds within the joint. The flow of the dissimilar aluminum alloys at the interface was controlled with a partially threaded friction stir welding pin and a dovetail geometric configuration optimized in previous work. The process robustness was explored to establish parameter windows by elucidating parametric effects on the relationship between process conditions and the joint’s microstructure and mechanical properties. Ranges of process forge force (50–70 kN) and interface temperature (475°C–515°C) consistently yielded joints with adequate performance. The load-carrying capacity in lap shear tensile tests was 2025 ? 100 N/mm.

Friction Stir Dovetailing, FSW, AA7099, AA6061, Ro↗
DOE OSTI2023

Microscale Electrochemical Corrosion of Uranium Oxide Particles

Understanding the corrosion of spent nuclear fuel is important for the development of long-term storage solutions. However, the risk of radiation contamination presents challenges for experimental analysis. Adapted from the system for analysis at the liquid–vacuum interface (SALVI), we developed a miniaturized uranium oxide (UO 2 )-attached working electrode (WE) to reduce contamination risk. To protect UO 2 particles in a miniatured electrochemical cell, a thin layer of Nafion was formed on the surface. Atomic force microscopy (AFM) shows a dense layer of UO 2 particles and indicates their participation in electrochemical reactions. Particles remain intact on the electrode surface with slight redistribution. X-ray photoelectron spectroscopy (XPS) reveals a difference in the distribution of U(IV), U(V), and U(VI) between pristine and corroded UO 2 electrodes. The presence of U(V)/U(VI) on the corroded electrode surface demonstrates that electrochemically driven UO 2 oxidation can be studied using these cells. Our observations of U(V) in the micro-electrode due to the selective semi-permeability of Nafion suggest that interfacial water plays a key role, potentially simulating a water-lean scenario in fuel storage conditions. This novel approach offers analytical reproducibility, design flexibility, a small footprint, and a low irradiation dose, while separating the α-effect. This approach provides a valuable microscale electrochemical platform for spent fuel corrosion studies with minimal radiological materials and the potential for diverse configurations.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
NASA NTRS1999

Development of a 12-Thrust Chamber Kerosene /Oxygen Primary Rocket Sub-System for an Early (1964) Air-Augmented Rocket Ground-Test System

Airbreathing/Rocket combined-cycle, and specifically rocket-based combined- cycle (RBCC), propulsion systems, typically employ an internal engine flow-path installed primary rocket subsystem. To achieve acceptably short mixing lengths in effecting the "air augmentation" process, a large rocket-exhaust/air interfacial mixing surface is needed. This leads, in some engine design concepts, to a "cluster" of small rocket units, suitably arrayed in the flowpath. To support an early (1964) subscale ground-test of a specific RBCC concept, such a 12-rocket cluster was developed by NASA's Marshall Space Flight Center (MSFC). The small primary rockets used in the cluster assembly were modified versions of an existing small kerosene/oxygen water-cooled rocket engine unit routinely tested at MSFC. Following individual thrust-chamber tests and overall subsystem qualification testing, the cluster assembly was installed at the U. S. Air Force's Arnold Engineering Development Center (AEDC) for RBCC systems testing. (The results of the special air-augmented rocket testing are not covered here.) While this project was eventually successfully completed, a number of hardware integration problems were met, leading to catastrophic thrust chamber failures. The principal "lessons learned" in conducting this early primary rocket subsystem experimental effort are documented here as a basic knowledge-base contribution for the benefit of today's RBCC research and development community.

Pryor, D.↗
DOE OSTI2020

An in situ shearing x-ray measurement system for exploring structures and dynamics at the solid–liquid interface

Revealing interfacial structure and dynamics has been one of the essential thematic topics in material science and condensed matter physics. Synchrotron-based x-ray scattering techniques can deliver unique and insightful probing of interfacial structures and dynamics, in particular, in reflection geometries with higher surface and interfacial sensitivity than transmission geometries. We demonstrate the design and implementation of an in situ shearing x-ray measurement system, equipped with both inline parallel-plate and cone-and-plate shearing setups and operated at the advanced photon source at Argonne National Laboratory, to investigate the structures and dynamics of end-tethered polymers at the solid-liquid interface. With a precise lifting motor, a micrometer-scale gap can be produced by aligning two surfaces of a rotating upper shaft and a lower sample substrate. A torsional shear flow forms in the gap and applies tangential shear forces on the sample surface. The technical combination with nanoscale rheology and the utilization of in situ x-ray scattering allow us to gain fundamental insights into the complex dynamics in soft interfaces under shearing. In this work, we demonstrate the technical scope and experimental capability of the in situ shearing x-ray system through the measurements of charged polymers at both flat and curved interfaces upon shearing. Furthermore, through the in situ shearing x-ray scattering experiments integrated with theoretical simulations, we aim to develop a detailed understanding of the short-range molecular structure and mesoscale ionic aggregate morphology, as well as ion transport and dynamics in soft interfaces, thereby providing fundamental insight into a long-standing challenge in ionic polymer brushes with a significant technological impact.

36 MATERIALS SCIENCE↗
NASA NTRS1996

The Micromechanics of the Moving Contact Line

The proposed research is divided into three components concerned with molecular structure, molecular orientation, and continuum averages of discrete systems. In the experimental program, we propose exploring how changes in interfacial molecular structure generate contact line motion. Rather than rely on the electrostatic and electrokinetic fields arising from the molecules themselves, we augment their interactions by an imposed field at the solid/liquid interface. By controling the field, we can manipulate the molecular structure at the solid/liquid interface. In response to controlled changes in molecular structure, we observe the resultant contact line motion. In the analytical portion of the proposed research we seek to formulate a system of equations governing fluid motion which accounts for the orientation of fluid molecules. In preliminary work, we have focused on describing how molecular orientation affects the forces generated at the moving contact line. Ideally, as assumed above, the discrete behavior of molecules can be averaged into a continuum theory. In the numerical portion of the proposed research, we inquire whether the contact line region is, in fact, large enough to possess a well-defined average. Additionally, we ask what types of behavior distinguish discrete systems from continuum systems. Might the smallness of the contact line region, in itself, lead to behavior different from that in the bulk? Taken together, our proposed research seeks to identify and accurately account for some of the molecular dynamics of the moving contact line, and attempts to formulate a description from which one can compute the forces at the moving contact line.

Han, Minsub↗
NASA NTRS1986

Ice banding as a response of the coupled ice-ocean system to temporally varying winds

This study models formation of ice bands in the marginal ice zones. A one-dimensional coupled ice-ocean model is used in which the ice model is coupled to a reduced gravity ocean model through interfacial stresses. The internal ice stresses are important only at high ice concentrations (90-100 percent); otherwise, the main balance for the ice motion is between the air-ice stress and the ice-water stress, i.e., free drift. The drag coefficients were chosen so that the air-ice momentum flux is 3 times greater than the air-ocean momentum flux. Thus the Ekman transport is larger under the ice than in the open water, so that winds parallel to the ice edge, with the ice on the right, produce upwelling. The upwelling simulation was extended to include temporally varying forcing, which was chosen to vary sinusoidally with a 4-day period. This forcing resembles successive cyclone passings perpendicular to the ice edge. When the oceanic upper layer was thin, which means that the dynamics are strongly nonlinear, the ice bands were formed. The up/downwelling signals do not disappear in wind reversals because of nonlinear advection. This leads to convergences and divergences in oceanic and ice velocities that manifest themselves as ice banding. At least one wind reversal is needed to produce one ice band.

Hakkinen, S.↗
DOE OSTI2020

Visualizing pectin polymer-polymer entanglement produced by interfacial water movement

In this report, we investigated the physical conditions for creating pectin polymer-polymer (homopolymer) entanglement. The potential role of water movement in creating pectin entanglement was investigated by placing water droplets-equivalent to the water content of two gel phase films-between two glass phase films and compressing the films at variable probe velocities. Slow probe velocity (0.5 mm/sec) demonstrated no significant debonding. Furthermore, corresponding videomicroscopy demonstrated an occasional water bridge, but no evidence of stranding or polymer entanglement. In contrast, fast probe velocity (5 mm/sec) resulted in 1) an increase in peak adhesion strength, 2) a progressive debonding curve, and 3) increased work of cohesion (p < .001). Corresponding videomicroscopy demonstrated pectin stranding and delamination between pectin films. Scanning electron microscopy images obtained during pectin debonding provided additional evidence of both stranding and delamination. We conclude that water movement can supply the motive force for the rapid chain entanglement between pectin films.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Mathematical modeling of novel porous transport layer architectures for proton exchange membrane electrolysis cells

Thin foil based porous transport layers (PTLs) that contain highly structured pore arrays have shown promise as anode PTLs in proton exchange membrane electrolysis cells. These novel PTLs, fabricated with advanced manufacturing techniques, produce thin, tunable, multifunctional layers with reduced flow and interfacial resistances and high thermal and electric conductivities. To further optimize their design, it is important to understand their fundamental impact on the transport of protons, electrons, and liquid/vapor mixtures in the electrode. In this work, we develop a two-dimensional multiphysics model to simulate the coupled electrochemistry and multiphase transport in an electrolysis cell operated with the novel PTL architecture. The results show that larger pores improve access of water to the anode catalyst layer, which is beneficial for both the oxygen evolution reaction and membrane hydration. Larger pore sizes also improve oxygen gas transport from the catalyst layer, because generated oxygen gas is forced to travel in-plane through the anode catalyst layer until it reaches a pore opening that is connected to a channel. The discussed results confirm that the proposed thin foil based PTLs are fundamentally different from conventional PTLs, such as felts or layered meshes. The model developed in this work also provides generalizable insight into fundamental PEMEC phenomena, such as the competition between liquid and gas phase transport, membrane hydration and water management, and nonuniform electrochemical reactions, which are processes relevant to all PEMEC designs.

25 ENERGY STORAGE↗
DOE OSTI2025

Tailoring the Physicochemical Properties of Nb Thin Films via Surface Engineering Methods

The modification of surface oxide layers formed on niobium (Nb) thin films via chemical mechanical planarization (CMP) and accelerated neutral atom beam (ANAB) processing provides a promising route toward tailoring their emergent properties and performance when used as superconducting qubits. Here, in this study, we show that CMP- and ANAB-formed Nb oxides are significantly thinner and smoother than the native oxide, as revealed by transmission electron microscopy (TEM) and atomic force microscopy. Scanning TEM and energy-dispersive X-ray spectroscopy along with X-ray photoelectron spectroscopy identified an oxidation gradient within the native and surface-engineered oxides. The topside layer is dominated by Nb 5+ (Nb 2 O 5 ), with various Nb suboxides present closer to the oxide/metal interface. Time-of-flight secondary ion mass spectrometry (ToF-SIMS) depth profiling confirmed the presence of an oxygen content gradient and demonstrated the enhanced resistance of the CMP- and ANAB-formed oxides to oxygen surface exchange and subsequent diffusion via 18 O 2 isotopic labeling experiments. ToF-SIMS also identified an interfacial layer containing trapped hydrogen (H)-containing species at the Nb oxide/metal interface. In situ ToF-SIMS and TEM revealed migration of the H/OH interfacial layer coinciding with decomposition of the surface oxide. Furthermore, our density functional theory calculations indicated that both H from moisture present in ambient air and bulk H in Nb films tend to segregate at the interface. These findings underscore the importance of understanding surface oxidation mechanisms, hydrogen incorporation, and their impact on the designed functionalities of Nb-based devices.

36 MATERIALS SCIENCE↗
NASA NTRS1994

Two-layer viscous instability in a rotating couette device

A novel experiment to study the interfacial shear instability between two liquids is described. Density-matched immiscible liquids are confined between concentric cylinders such that the interface is parallel to the cylinder walls. Interfacial waves that develop because of viscosity differences between the shearing fluids are studied as a function of rotation rate and depth ratio using optical techniques. Conditions neutral stability and the most unstable wavenumber agree reasonably well with predictions from linear stability analysis of the Navier-Stokes equations. Illumination using laser sheets allows precise measurement of the interface shape. Future experiments will verify the correctness of weakly nonlinear theories that describe energy transfer and saturation of wave growth by nonlinear effects. Measurements of solitary wave shapes, that occur far above neutral stability, will be compared to similar measurements for systems that have gravity as an important force to determine how gravity effects large disturbances. These results will be used to interpret slug and annular flow data that have been obtained in other mu g studies.

Gallagher, Christopher↗
NASA NTRS2000

Measurement of the Surface Dilatational Viscosity of an Insoluble Surfactant Monolayer at the Air/Water Interface Using a Pendant Drop Apparatus

When a fluid interface with surfactants is at rest, the interfacial stress is isotropic (as given by the equilibrium interfacial tension), and is described by the equation of state which relates the surface tension to the surfactant surface concentration. When surfactants are subjected to shear and dilatational flows, flow induced interaction of the surfactants; can create interfacial stresses apart from the equilibrium surface tension. The simplest relationship between surface strain rate and surface stress is the Boussinesq-Scriven constitutive equation completely characterized by three coefficients: equilibrium interfacial tension, surface shear viscosity, and surface dilatational viscosity Equilibrium interfacial tension and surface shear viscosity measurements are very well established. On the other hand, surface dilatational viscosity measurements are difficult because a flow which change the surface area also changes the surfactant surface concentration creating changes in the equilibrium interfacial tension that must be also taken into account. Surface dilatational viscosity measurements of existing techniques differ by five orders of magnitude and use spatially damped surface waves and rapidly expanding bubbles. In this presentation we introduce a new technique for measuring the surface dilatational viscosity by contracting an aqueous pendant drop attached to a needle tip and having and insoluble surfactant monolayer at the air-water interface. The isotropic total tension on the surface consists of the equilibrium surface tension and the tension due to the dilation. Compression rates are undertaken slow enough so that bulk hydrodynamic stresses are small compared to the surface tension force. Under these conditions we show that the total tension is uniform along the surface and that the Young-Laplace equation governs the drop shape with the equilibrium surface tension replaced by the constant surface isotropic stress. We illustrate this technique using DPPC as the insoluble surfacant monolayer and measured for it a surface dilatational viscosity in the LE phase that is 20 surface poise.

Lorenzo, Jose↗
DOE OSTI2021

Enhanced van der Waals epitaxy of germanium by out-of-plane dipole moment induced from transferred graphene on TiN/AlN multilayers

We report recent advances in 3D/2D heterostructures have opened up tremendous opportunities in building highly flexible and durable optoelectronic devices. However, the inherit lack of interfacial bonding and low surface energy of van der Waals surfaces limit the nucleation and growth of 3D materials. Enhancing wettability by providing a porous buffer is effective in growing compound semiconductors on graphene while van der Waals epitaxy of Ge remains challenging. Here, the nucleation of Ge has been significantly improved from an islanded mode to granular modes by using a TiN/AlN multilayered buffer prior to Ge/graphene integration. Highly textured Ge growth with dominating (111), (220), and (311) peaks are identified by x-ray diffraction. The microstructure of the buffer TiN/AlN demonstrates a polycrystalline quality with clean interfaces between each interlayer and the substrate. Kelvin probe force microscopy measurements along the lateral TiN/AlN interface identify a potential drop corresponding to the AlN phase. This contact potential difference between TiN and AlN is the key in generating the out-of-plane dipole moment that modifies the surface energy of the monolayer graphene, resulting in enhanced wettability of the Ge adatoms nucleated on top. Surface dipole induced nucleation of 3D semiconductor thin films on 2D materials via the proper design of buffer layer is fundamentally important to enhance the 3D/2D growth toward flexible optoelectronic applications.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Mesocrystal growth through oriented sliding and attachment of nanoplates

Oriented attachment is a critical, yet poorly understood, crystal growth pathway based on the self-assembly of nanocrystals. During oriented attachment, solvent-separated particles align and coalesce through forces that enable precise rotation and translation. While prior studies emphasized intragap forces driving crystallographic alignment, the forces enabling uniform stacking and superlattice formation remain unclear. Here, we demonstrate how macroscopic gibbsite mesocrystals emerge from nanoplates guided into staggered positions by directional sliding. Electron microscopy and X-ray scattering reveal the monoclinic superlattice structure, based on nanoplate stacking with a uniform ≈50° stagger along the gibbsite [010] direction. In situ liquid-cell TEM captures preferential sliding along the gibbsite [010] direction, decelerating with increasing particle overlap. Molecular dynamics simulations reveal that this staggered arrangement corresponds to a global free-energy minimum, rather than full alignment. The simulations also confirm that sliding along the [010] direction is energetically favored and provide insight into the role of interfacial water in achieving long-range ordered assemblies. These insights highlight the energy landscape’s role in oriented attachment, with implications for material synthesis and hierarchical structures in nature.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Electrolyte–Electrocatalyst Interfacial Effects of Polymeric Materials for Tandem CO 2 Capture and Conversion Elucidated Using In Situ Electrochemical AFM

Integrating CO 2 capture and electrochemical conversion has been proposed as a strategy to reduce the net energy required for CO 2 regeneration in traditional CO 2 capture and conversion schemes and can be coupled with carbon-free renewable electricity. Polyethylenimine (PEI)-based materials have been previously studied as CO 2 capture materials and can be integrated in these reactive capture processes. PEI-based electrolytes have been found to significantly increase the CO 2 loading, and impact selectivity and rate of product formation when compared to the conventional aqueous electrolytes. However, the influence of these materials at the catalyst–electrode interface is currently not well understood. In this study, PEI-based electrolytes were prepared and their impact on the morphology of a silver electrode performing electrochemical CO 2 reduction (CO 2 R) was studied using in situ electrochemical atomic force microscopy (EC-AFM). The presence of PEI on the electrode surface could be distinguished based on nanomechanical properties (DMT modulus), and changes were observed as negative polarization was applied, revealing a reorganization of the PEI chains due to electrostatic interactions. These changes were impacted by the electrolyte composition, including the addition of supporting electrolyte KHCO 3 salt, as well as CO 2 captured by the PEI-based electrolyte, which minimized the change in surface mechanical properties and degree of PEI alignment on the electrode surface. The changes in surface mechanical properties were also dependent on the PEI polymer length, with higher molecular weight PEI showing different reconfiguration than the shorter polymer brushes. The study highlights that the choice of polymer material, the electrolyte composition, and CO 2 captured impact the nearelectrode environment, which has implications for CO 2 R, and presents EC-AFM as a new tool that can be used to probe the dynamic behavior of these interfaces during electrocatalysis.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Modeling Equilibrium Solid–Liquid Interfaces under Effective Constant Chemical Potential Using Machine Learning Interatomic Potentials

The chemical potential (μ) of species in solution is essential for understanding various chemical processes at interfaces. Molecular dynamics (MD) simulations, constrained by fixed compositions, cannot maintain constant chemical potential with reference to a targeted concentration or chemical potential under nonequilibrium or dynamic conditions, as solute species can migrate to the interface and deplete (or enrich) the bulk due to solute-interface interactions. In this study, we introduce a simple and computationally efficient approach named iterative quasi-constant chemical potential molecular dynamics (iqCμMD) simulation, which helps simulate targeted molar concentrations of species in solution. iqCμMD overcomes the limitations of conventional MD by adjusting the number of species in the solution to reach a target bulk concentration (chemical potential), which allows simulation of the interface under the bulk conditions comparable to experiment. We demonstrate our approach using machine learning interatomic potential (MLIP)-based MD simulations of the Na 2 SO 4,aq –graphene interface, and to show the transferability of our approach, we also perform classical force field-based MD simulations of NaCl aq –air and NaCl aq –graphite interfaces, which produce comparable results to previous CμMD simulations. Our results also show that the iqCμMD approach efficiently achieves the desired bulk ion concentration within two iterations, and by utilizing MLIPs, we can achieve converged results using relatively small-scale simulations compared to previous CμMD simulations. By combining iqCμMD with MLIP-driven simulations, solid–liquid interfaces can be modeled under an effective constant chemical potential with DFT-level accuracy. Here, we show that iqCμMD offers a robust and simple computational framework for constant chemical potential simulations, as its only requirement is to be able to converge interfacial simulations with a measurable bulk region.

Chemical structure↗