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

Investigating the origin of the far-field reflection interference fringe (RIF) of microdroplets

We show that the reflection interference fringe (RIF) is formed on a screen far away from the microdroplets placed on a prism-based substrate, which have low contact angles and thin droplet heights, caused by the dual convex–concave profile of the droplet, not a pure convex profile. The geometric formulation shows that the interference fringes are caused by the optical path difference when the reflected rays from the upper convex profile at the droplet–air interface interfere with reflection from the lower concave profile at oblique angles lower than the critical angle. Analytic solutions are obtained for the droplet height and the contact angle out of the fringe number and the fringe radius in RIF from the geometric formulation. Furthermore, the ray tracing simulation is conducted using the custom-designed code. The geometric formulation and the ray tracing show excellent agreement with the experimental observation in the relation between the droplet height and the fringe number and the relation between the contact angle and the fringe radius. This study is remarkable as the droplet's dual profile cannot be easily observed with the existing techniques. However, the RIF technique can effectively verify the existence of a dual profile of the microdroplets in a simple setup. In this work, the RIF technique is successfully developed as a new optical diagnostic technique to determine the microdroplet features, such as the dual profile, the height, the contact angle, the inflection point, and the precursor film thickness, by simply measuring the RIF patterns on the far-field screen.

42 ENGINEERING↗
DOE OSTI2022

Molecular Driving Force for Facet Selectivity of Sequence-Defined Amphiphilic Peptoids at Au–Water Interfaces

Shape-controlled colloidal nanocrystal syntheses often require facet-selective solution-phase chemical additives to regulate surface free energy, atom addition/migration fluxes or particle attachment rates. Because of their highly tunable properties and robustness to a wide range of experimental conditions, peptoids represent a very promising class of next generation functional additives for control over nanocrystal growth. However, understanding the origin of facet selectivity at the molecular level is critical to generalizing their design. In this work, we employ molecular dynamics simulations and biased sampling methods and report stronger selectivity to Au(111) than to Au(100) for Nce3Ncp6, a peptoid that has been shown to assist the formation of five-fold twinned Au nanostars. We find that facet-selectivity is achieved through synergistic effects of both peptoid-surface and solvent-surface interactions. Moreover, the amphiphilic nature of Nce3Ncp6 together with the order of peptoid-peptoid and peptoid-surface binding energies, i.e., peptoid-Au(100) < peptoid-peptoid < peptoid-Au(111), further amplifies its distinct collective behavior on different Au surfaces. Our studies provide a fundamental understanding of the molecular origin of facet-selective adsorption and highlight the possibility of future designs and uses of sequence-defined peptoids for predictive syntheses of nanocrystals with designed shapes and properties.

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

Y(III) Sorption at the Orthoclase (001) Surface Measured by X-ray Reflectivity

Interactions of heavy metals with charged mineral surfaces control their mobility in the environment. Here, we investigate the adsorption of Y(III) onto the orthoclase (001) basal plane, the former as a representative of rare earth elements and an analogue of trivalent actinides and the latter as a representative of naturally abundant K-feldspar minerals. In this study, we apply in situ high-resolution X-ray reflectivity to determine the sorption capacity and molecular distribution of adsorbed Y species as a function of the Y 3+ concentration, [Y 3+ ], at pH 7 and 5. With [Y 3+ ] ≥ 1 mM at pH 7, we observe an inner-sphere (IS) sorption complex at a distance of ~1.5 Å from the surface and an outer-sphere (OS) complex at 3–4 Å. Based on the adsorption height of the IS complex, a bidentate, binuclear binding mode, in which Y 3+ binds to two terminal oxygens, is proposed. In contrast, mostly OS sorption is observed at pH 5. The observed maximum Y coverage is ~1.3 Y 3+ /A UC (A UC : area of the unit cell = 111.4 Å 2 ) for all the investigated pH values and Y concentrations, which is in the expected range based on the estimated surface charge of orthoclase (001).

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

Differences in Mechanism and Rate of Zeolite-Catalyzed Cyclohexanol Dehydration in Apolar and Aqueous Phase

The rate of acid-base catalyzed dehydration of alcohols strongly depends on the solvent and the environment of the acid sites. In this work, we find that Brønsted acidic sites in large-pore zeolites, but not in medium-pore zeolites, catalyze cyclohexanol dehydration in decalin at significantly higher rates than hydrated hydronium ions in aqueous phase. Specifically, the difference in turnover rates between the two solvents amounts to two to three orders of magnitude on H-BEA and H-FAU, while being very modest (within a factor of 2) for H-MFI. Combining kinetic, isotopic tracer and 2H NMR measurements, it is established that cyclohexanol dehydration generally follows an E1-elimination pathway in decalin. A notable exception is the monomer dehydration route on H-MFI, which exhibits a much lower activation energy and a substantially negative activation entropy that appear to be associated with an E2-type mechanism. The C-O bond cleavage displays a dominant degree of rate control in decalin, which stands in contrast to deprotonation (C-H cleavage) being rate-limiting in aqueous-phase dehydration.

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

High-speed impact of micron-sized diesel drop trains—Splashing dynamics, secondary droplet formation, and effects of pre-existing film thickness

This work investigates surface impingement of monodispersed diesel drop trains using computational simulations. The three-dimensional simulations include high-speed impact of micrometer-sized fuel drops onto initially dry and wetted stainless steel substrates. The drop size and impact velocity are representative of fuel injection conditions in internal combustion engines (ICEs). The drop trains serve as a simplified representation of fuel spray. To accurately capture the film that naturally develops on initially dry substrates, a dynamic diesel-stainless steel contact angle model was developed using experiments of single drop impact under ICE representative temperature. Simulations are presented for a highly splashing case and the effects of a pre-existing film on splashing dynamics were investigated, including the temporal evolution of splashed mass and film thickness. It was concluded that for thin films, the effect of pre-existing film thickness is less significant than in single drop impingement. Secondary droplet characterization was performed on simulation results, leading to instantaneous and time-averaged distributions of secondary droplet size, velocity magnitude, and trajectory angle. It was found that for each drop impingement, approximately 58% of the splashed mass is from the impinging drop itself, while the remainder of the splashed mass is composed of film liquid. From a high resolution case, details of secondary droplet formation are observed and three distinct phases of secondary droplet formation are identified. The detailed analysis of drop train impingement under engine-relevant conditions serves as a first step toward a robust understanding of fuel impingement and the development of cleaner and more efficient ICEs.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Less can be more: Insights on the role of electrode microstructure in redox flow batteries from two-dimensional direct numerical simulations

Understanding how to structure a porous electrode to facilitate fluid, mass, and charge transport is key to enhancing the performance of electrochemical devices, such as fuel cells, electrolyzers, and redox flow batteries (RFBs). Here, using a parallel computational framework, direct numerical simulations are carried out on idealized porous electrode microstructures for RFBs. Strategies to improve an electrode design starting from a regular lattice are explored. By introducing vacancies in the ordered arrangement, it is possible to achieve higher voltage efficiency at a given current density, thanks to improved mixing of reactive species, despite reducing the total reactive surface. Careful engineering of the location of vacancies, resulting in a density gradient, outperforms disordered configurations. Our simulation framework is a new tool to explore transport phenomena in RFBs, and our findings suggest new ways to design performant electrodes.

25 ENERGY STORAGE↗
DOE OSTI2022

Phase diagram of the TIP4P/Ice water model by enhanced sampling simulations

In this work, we studied the phase diagram for the TIP4P/Ice water model using enhanced sampling molecular dynamics simulations. Our approach is based on the calculation of ice–liquid free energy differences from biased coexistence simulations that reversibly sample the melting and growth of layers of ice. We computed a total of 19 melting points for five different ice polymorphs, which are in excellent agreement with the melting lines obtained from the integration of the Clausius–Clapeyron equation. For proton-ordered and fully proton-disordered ice phases, the results are in very good agreement with previous calculations based on thermodynamic integration. For the partially proton-disordered ice III, we find a large increase in stability that is in line with previous observations using direct coexistence simulations for the TIP4P/2005 model. This issue highlights the robustness of the approach employed here for ice polymorphs with diverse degrees of proton disorder. Our approach is general and can be applied to the calculation of other complex phase diagrams.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

A multimodal flow reactor for photocatalysis under atmospheric conditions

Photocatalysis is a promising concept for the direct conversion of solar energy into fuels and chemicals. The design, experimental protocol, and performance of a multimodal and versatile flow reactor for the characterization of powdered and immobilized photocatalysts are presented in this report. Ultimately, this instrument enables rigorous evaluation of photocatalysis performance metrics. The apparatus quantifies transient gas-phase reaction products via online real-time gas analyzer mass spectrometry (RTGA-MS). For H 2 , the most challenging gas, the photocatalytic system’s RTGA-MS gas detection sensitivity spans over three orders of magnitude and can detect down to tens of parts per million under atmospheric conditions. Using Pt nanoparticles supported on anatase TiO 2 photocatalyst via wet impregnation, the instrument’s capability for the characterization of photocatalytic H 2 evolution is demonstrated, resulting in an apparent quantum yield (AQY) of 48.1% ± 0.9% at 320 nm, 45.7% ± 0.3% at 340 nm and 31% ± 1% at 360 nm. The photodeposition of Pt on anatase TiO 2 was employed to demonstrate the instrument’s capability to track the transient behavior of photocatalysts, resulting in an improved 55% ± 2% AQY for H 2 evolution at 340 nm from aqueous methanol. This photocatalytic instrument enables systematic study of a wide variety of photocatalytic reactions such as water splitting and CO 2 reduction to valuable C2+ fuels and chemicals.

14 SOLAR ENERGY↗
DOE OSTI2025

Incorporating the molecular-scale into a hydrodynamic description of confined aqueous systems

Hydrodynamics provides a continuum-level description of fluid motion, but its applicability at the nanoscale becomes uncertain due to the emerging importance of molecular-level effects such as spatial heterogeneity. Hydrodynamic boundary conditions that incorporate molecular details allow us to partition the system into a near-wall region and a bulk fluid region. We identify a hydrodynamic wall located inside the fluid that determines where slip begins. By extending the hydrodynamic wall with the slip length, the position of the extrapolated wall is established. This offers a unified description of both slip and stagnant flow behaviors, with wall hydrophobicity characterized by the relative location of the extrapolated wall with respect to the physical wall. Employing this concept in analyses of equilibrium molecular dynamics (MD) and non-equilibrium MD simulations of Couette and Poiseuille flows, our results demonstrate consistency between equilibrium and non-equilibrium approaches across different flow types and confinement levels. This demonstrates the robust nature of linear response theory. We then explore the effects of fluid-wall and bulk fluid interactions on the hydrodynamic properties. These findings enhance the effectiveness of molecular-based simulations for investigating complex confined systems in nanofluidics, biology, and colloidal science, offering a complementary molecular-scale perspective to traditional continuum approaches.

flow simulations↗
DOE OSTI2025

Structure, dynamics, and electrochemistry of choline chloride/ethylene glycol eutectic solvents at an electrode surface explored by molecular dynamics simulations

Choline chloride and ethylene glycol mixtures with 1:2, 1:4, and 1:6 molar ratios on the surfaces of graphite and gold electrodes were studied using classical molecular dynamics simulations. Both neutral and charged electrodes were considered. The liquid composition, solvation structure, molecular orientation, and dynamics at the electrode surface are significantly different from those of the bulk liquid. These properties strongly depend on the electrode material and charge density, whereas they are less sensitive to the overall solvent composition. The effect of the electrode on the composition, structure, and orientation of the liquid fades beyond ∼10 Å from the surface of the electrode. This distance corresponds to about two layers of the structured solvent, despite the fact that the layered structure extends to at least five layers or about 25 Å. However, the electrode influences solvent dynamics over a longer distance. The electrochemical properties of the eutectic solvent at both electrode surfaces were also studied. The simulations captured the experimental differential capacitance shapes for both electrode systems, although the magnitudes and exact shapes differ. The simulations further revealed that the solvent in the first solvation layer can both overscreen and underscreen the electrode charges depending on the electrode material and electrode potential.

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

Advances in Operando Electrocatalysis with High-Resolution Hard X-Ray Spectroscopy

Electrocatalysis unfolds at the interface—where solids, liquids, and gases meet to exchange charge, transform molecules, and set the foundation for technologies like CO 2 conversion, and water splitting. These interfacial regions are inherently dynamic, chemically diverse, and structurally heterogeneous. How atoms rearrange, oxidize, coordinate ligands, or adsorb reaction intermediates under electro-chemical potential defines a catalyst’s performance. Furthermore, capturing these changes with element specificity and electronic sensitivity under operando conditions remains one of the most critical challenges in the field.

Garcia-Esparza, Angel T. [SLAC National Accelerato↗
DOE OSTI2022

Machine learning approaches to streamline and enhance the analysis of multiscale imaging data for bioaerosol and soil particles

Bioaerosol and soil particles are ubiquitous in the environment. They are multicomponent and complex in nature displaying mixed inorganic and organic components. The way components are mixed in a bioaerosol sample is referred to as its mixing state. Soil particles are also a mixture of inorganic (mineral) and organic (soil organic matter) components. Bioaerosol particles contribute to a major fraction of coarse mode atmospheric particles, especially in the tropical areas, contributing up to 80 % of the particle mass concentration. The mixing state of particles is crucial to evaluate because it impacts several important environmental processes such as warm and cold cloud formation and radiation budget. Mixing states in aerosols are accompanied by chemical reactions across solid-liquid-gas interfaces. In this study, we utilized elemental compositions and microcopy images of thousands of atmospheric particles acquired by computer-controlled scanning electron microscope equipped with an energy-dispersive x-ray spectrometer to compute the mixing state of atmospheric particles. A 2D convolutional neural network (CNN), also known as convnet, was used to model the relationship between low resolution imaging data and higher resolution spectroscopy data, with the former as training input and the latter as target output. Two types of CNNs were implemented and tested; a basic CNN and an Inception-v3 network. For binary classification, the basic CNN achieved an accuracy of 84.29 % across all atom types, and the Inception-v3-like network achieved an accuracy 85.51 %. This study demonstrates the applicability of deep learning to handle large amounts of imaging/chemical spectroscopy data efficiently and evaluate particle mixing state from a range of environmental samples.

54 ENVIRONMENTAL SCIENCES↗