Engineering Papers⌕ Search

Engineering topics

Wang, Xiaoyu

Publications and source records attributed to Wang, Xiaoyu.

pH-Driven Restructuring of Hydration Layers and Cation Ad-sorption at the Alumina-Water Interface

Oxide-water interfaces underpin ion separation, catalysis, and electrochemical energy technologies, where the electrical double layer (EDL) controls adsorption, transport, and reactivity. Yet, the molecular-scale link between pH-dependent surface protonation, hydration-layer structure, and counter-ion adsorption remains poorly defined. Here, we combine in situ crystal truncation rod (CTR) and resonant anomalous X-ray reflectivity (RAXR) with streaming potential measurements and ab initio molecular dynamics (AIMD) simulations to resolve the chemical and structural evolution of the EDL at the single-crystal alumina (012)-water interface in 10 mM Rb+ over pH 3-12. CTR measurements reveal two distinct adsorbed water layers at ~2.2 and ~3.5 Å above the surface that each shift toward the substrate at transition pHs near 6.5 and 10.6, respectively, directly reflecting changes in primary hydration layer structure in response to the deprotonation of bridging and terminal aluminol groups. RAXR shows a 10-fold increase in Rb+ coverage and a decrease in mean adsorption height from ~3.5 to ~2.7 Å with increasing pH, indicating enhanced counter-ion binding accompanied by Stern layer contraction. Streaming potential measurements demonstrate that the zeta potential, i.e., potential at the hydrodynamic shear plane, is positive at pH 3 and becomes negative at pH ≥3.5. This negative charge magnitude increases with increasing pH, consistent with progressive surface deprotonation at higher pH. AIMD identifies inner- and outer-sphere Rb+ complexes whose adsorption heights and coordination geometries depend sensitively on the protonation state of surface oxygens, providing atomistic support for the experimentally inferred trends. These measurements establish two discrete, site-specific pH transitions in hydration-layer structure that track aluminol (de)protonation and quantitatively link them to a pH-driven contraction of the Stern layer (increasing Rb+ coverage and decreasing adsorption height). This provides a direct structural basis for connecting surface acid-base chemistry to ion binding distances at an oxide-water interface.

Electrical double layer (EDL), Surface protonation↗

High Pressure Synthesis of Rubidium Superhydrides

Through laser-heated diamond anvil cell experiments, we synthesize a series of rubidium superhydrides and explore their properties with synchrotron x-ray powder diffraction and Raman spectroscopy measurements, combined with density functional theory calculations. Upon heating rubidium monohydride embedded in H 2 at a pressure of 18 GPa, we form RbH 9 − I , which is stable upon decompression down to 8.7 GPa, the lowest stability pressure of any known superhydride. At 22 GPa, another polymorph, RbH 9 − II is synthesised at high temperature. Unique to the Rb-H system among binary metal hydrides is that further compression does not promote the formation of polyhydrides with higher hydrogen content. Instead, heating above 87 GPa yields RbH 5 , which exhibits two polymorphs ( RbH 5 − I and RbH 5 − II ). All of the crystal structures comprise a complex network of quasimolecular H 2 units and H − anions, with RbH 5 providing the first experimental evidence of linear H 3 − anions. Published by the American Physical Society 2025

Kuzovnikov, Mikhail A.↗

Metadynamics investigation of lanthanide solvation free energy landscapes and insights into separations energetics

Lanthanide ion solvation chemistry in nonaqueous phases is key to understanding and developing effective separation processes for these critical materials. Due to the complexity and inherent disorder of the solution phase, a comprehensive picture of the solvated metal ion is often difficult to generate solely from conventional spectroscopic approaches and electronic structure calculations, particularly in the extractant phase. In this work, we use classical molecular dynamics (MD) simulation with an advanced sampling technique, metadynamics, supplemented by experimental spectroscopy and speciation analysis, to measure lanthanide solvation free energy landscapes. We define coordination-based collective variables to probe the entire range of solvation configurations in the organic phase of lanthanum (La), europium (Eu), and lutetium (Lu) nitrate salts bound with a commonly used extractant, N,N′-dimethyl, N,N′-dioctylhexylethoxymalonamide (DMDOHEMA). The known lanthanide extraction trend of La ≈ Eu > Lu is readily explained by the measured free energy surfaces, which show consistent DMDOHEMA coordination from La to Eu, followed by loss of DMDOHEMA coordination from Eu to Lu. These simulations suggest how ligand crowding at the metal center can control selectivity, in this case resulting in the opposite extraction trend as observed with other conventional extractants, where the enthalpic contribution from increasing lanthanide charge density across the series dominates the extraction energetics. We also find that the presence of inner-sphere water, verified by time-resolved fluorescence, diversifies the accessible solvation structures. As a result, understanding solvation requires consideration of an entire thermodynamic ensemble, rather than the single dominant lowest-energy structure, as is often considered out of necessity in interpretation of spectroscopic data or in electronic structure-based ligand design approaches. In general, we demonstrate how metadynamics uniquely enables investigation of complex, multidimensional solvation energetic landscapes, and how it can explain selectivity trends where extraction is controlled by more complex mechanisms than simple charge density-based selectivity.

Wang, Xiaoyu↗

Adsorption and Solvation Modulate Rhodamine B Diffusion in Ethanol/Water-Filled Anodic Alumina Nanopores

Confinement of solvents and solutes within nanoporous materials frequently leads to the emergence of unique mass transport behaviors that, once fully understood, may lead to improved chemical separations. Here, the diffusion of Rhodamine B (RhB) dye within 10 and 20 nm diameter anodic aluminum oxide (AAO) nanopores filled with binary ethanol/water mixtures is investigated. Mixture compositions spanning from pure ethanol to pure water are employed. The results of confocal fluorescence correlation spectroscopy studies reveal that RhB diffusion occurs by a two-component mechanism comprising composition-dependent fast and slow motions, characterized by diffusion coefficients D f and D s . The results are consistent with those of previous studies performed under more limited conditions [J. Phys. Chem. C, 2023, 127, 411-420]. The fast component scales with mixture viscosity and is assigned to hindered bulk-like diffusion in central pore regions. Slow diffusion likely involves adsorption of RhB to the pore surface and may be described by a desorption mediated mechanism. The occurrence of RhB adsorption to the AAO surface is verified at the single-molecule level by wide-field fluorescence imaging of membrane cross-sectional surfaces. Unique composition dependent trends in the autocorrelation amplitude and in D s that mimic bulk RhB solubility are revealed. D s is found to be smallest in pure ethanol and pure water and largest in intermediate mixtures. These results suggest that RhB surface adsorption is strongest in the pure liquids and weakest in mixtures of intermediate composition, where the dye is least soluble, and most soluble, respectively. As a result, molecular dynamics simulations reveal that a water layer appears on the pore surface under most conditions, while RhB is solvated primarily by ethanol. The composition dependence of RhB diffusion is concluded to reflect its solvation dependent interactions with the pore walls.

36 MATERIALS SCIENCE↗

Influence of Aqueous Phase Acidity on Ln(III) Coordination by N , N , N ', N '-Tetraoctyldiglycolamide

Here, this study highlights the importance of combining distribution ratio measurements with multiple spectroscopic techniques to provide a more comprehensive understanding of organic phase Ln coordination chemistry. Solvent extraction investigations with N,N,N',N'-tetraoctyldiglycolamide (TODGA) in n-heptane reveal the sensitivity of Ln complexation to the HNO 3 concentration. Distribution ratio measurements in tandem with UV–Vis demonstrated that increasing the concentration of HNO3 above 0.5 M with a constant NO 3 – of 1 M increases the number of coordinating TODGA molecules, from a 1:2 to a 1:3 Ln:TODGA complex. At each concentration of HNO 3 considered herein (from 0.01 to 1 M), Eu lifetime analysis demonstrated no evidence of H 2 O coordination. Results from Fourier transform infrared investigations suggest the presence of inner-sphere NO 3 – under low concentrations of HNO 3 when the 1:2 Ln:TODGA complex is present. Increasing the HNO 3 concentration above 0.5 M increases the propensity for outer-sphere interactions by removing the coordinated NO 3 – and saturating the Ln coordination sphere with three TODGA molecules, resulting in the well-established cationic, trischelate homoleptic [Ln(TODGA) 3 ] 3+ complex. This work demonstrates the importance in considering the NO 3 – source for solvent extraction systems. In particular, for systems with an affinity for outer-sphere interactions with molar concentrations of HNO 3 , changing the NO 3 – source can change the inner-sphere coordination of the Ln complex, which, in turn, affects the separation efficacy.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Realization of Organocerium-Based Fullerene Molecular Materials Showing Mott Insulator-Type Behavior

Abstract Electron-rich organocerium complexes (C 5 Me 4 H) 3 Ce and [(C 5 Me 5 ) 2 Ce(ortho-oxa)], with redox potentials E 1/2 = –0.82 V and E 1/2 = –0.86 V versus Fc/Fc + respectively, were reacted with fullerene (C 60 ) in different stoichiometries to obtain molecular materials. Structurally characterized co-crystals: [(C 5 Me 4 H) 3 Ce] 2 •C 60 (1) and [(C 5 Me 5 ) 2 Ce(ortho-oxa)] 3 •C 60 (2) of C 60 with cerium-based rare earth molecular precursors are reported for the first time. The extent of charge transfer in 1 and 2 was evaluated using a series of physical measurements: FT-IR, Raman, solidstate UV-vis-NIR spectroscopy, X-ray absorption near edge structure (XANES) spectroscopy, and magnetic susceptibility measurements. The physical measurements indicate that 1 and 2 comprise the cerium(III) oxidation state with formally neutral C 60 as a co-crystal in both cases. Pressure-dependent periodic density functional theory calculations were performed to study the electronic structure of 1. Inclusion of a Hubbard-U parameter removes Ce f states from the Fermi level, opens up a band gap, and stabilizes FM/AFM magnetic solutions that are isoenergetic because of the large distances between the Ce(III) cations. Furthermore, the electronic structure of this strongly correlated Mott insulator-type system is reminiscent of the well-studied Ce 2 O 3 .

36 MATERIALS SCIENCE↗

Twisted-layer boron nitride ceramic with high deformability and strength

Moiré superlattices formed by twisted stacking in van der Waals materials have emerged as a new platform for exploring the physics of strongly correlated materials and other emergent phenomena. However, there remains a lack of research on the mechanical properties of twisted-layer van der Waals materials, owing to a lack of suitable strategies for making three-dimensional bulk materials. Here we report the successful synthesis of a polycrystalline boron nitride bulk ceramic with high room-temperature deformability and strength. This ceramic, synthesized from an onion-like boron nitride nanoprecursor with conventional spark plasma sintering and hot-pressing sintering, consists of interlocked laminated nanoplates in which parallel laminae are stacked with varying twist angles. The compressive strain of this bulk ceramic can reach 14% before fracture, about one order of magnitude higher compared with traditional ceramics (less than 1% in general), whereas the compressive strength is about six times that of ordinary hexagonal boron nitride layered ceramics. The exceptional mechanical properties are due to a combination of the elevated intrinsic deformability of the twisted layering in the nanoplates and the three-dimensional interlocked architecture that restricts deformation from propagating across individual nanoplates. The advent of this twisted-layer boron nitride bulk ceramic opens a gate to the fabrication of highly deformable bulk ceramics.

36 MATERIALS SCIENCE↗

Structural Trends and Vibrational Analysis of N,N,N′,N′−Tetramethylmalonamide Complexes Across the Lanthanide Series

Fundamental understanding of coordination chemistry across the lanthanide series is essential for explaining chemical behavior of rare-earth metals in complex liquid-liquid extraction processes, which in turn affects the distribution ratios and efficacy of separations as a whole. In this work, we explore the structural trends between the lanthanides and a neutral N,N,N',N'-tetramethylmalonamide (TMMA) ligand within four isolated families of solid-state compounds: Ln(trans-TMMA) 2 (NO 3 ) 3 Ln=La-Nd, Sm; Ln(cis-TMMA) 2 (NO 3 ) 3 Ln=Eu-Tb, Er; [Ln(TMMA) 3 (NO 3 ) 2 ][Ln(TMMA)(NO 3 ) 4 ] Ln=Dy-Tm; Ln(Κ 2 -TMMA)(iPrOH)(NO 3 ) 3 and Ln(Κ 1 -TMMA)(Κ 2 -TMMA)(NO 3 ) 3 Ln=Yb, Lu. Moving across the lanthanide series, we note the formation of both discrete charge-neutral complexes, as well as charged molecular anion-cation pairs, with variations in spatial ligand arrangement, coordination numbers, and ligand denticities. IR and Raman spectroscopy paired with DFT frequency calculations were used for an in-depth investigation of vibrational modes unique to each structural family. The collection of isolated model compounds was also discussed in the context of liquid-liquid separations based on reported distribution ratios from malonamide extraction.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Using Metadynamics to Reveal Extractant Conformational Free Energy Landscapes

Understanding the impact of extractant functionalization on metal-binding energetics in liquid-liquid extraction is essential to guide the development of better separation processes. Traditionally, computational extractant design uses electronic structure calculations on metal-ligand clusters to determine the metal-binding energy of the lowest energy state. Although highly accurate, this approach does not account for all of the relevant physics encountered under experimental conditions. Such methodologies often neglect entropic contributions such as temperature effects and ligand flexibility, in addition to approximating solvent-extractant interactions with implicit solvent models. In this study, we use classical molecular dynamics simulations with an advanced sampling method, metadynamics, to map out extractant molecule conformational free energies in the condensed phase. Here we generate the complete conformational landscape in solution for a family of bidentate malonamide-based extractants with different functionalizations of the headgroup and the side chains. In particular, we show how such alkyl functionalization reshapes the free energy landscape, affecting the free energy penalty of organizing the extractant into the cis-like metal-binding conformation from the trans-like conformation of the free extractant in solution. Specifically, functionalizing alkyl tails to the center of the headgroup has a greater influence on increasing molecular rigidity and disfavoring the binding conformation than functionalizing side chains. These findings are consistent with trends in metal-binding energetics based on experimentally reported distribution ratios. We also consider a different bidentate extractant molecule, carbamoylmethylphosphine oxide, and show how the choice of solvent can further reshape the conformational energetic landscape. This study demonstrates the feasibility of using molecular dynamics simulations with advanced sampling techniques to investigate extractant conformational energetics in solution, which, more broadly, will enable extractant design that accounts for entropic effects and explicit solvation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Solvent effects on extractant conformational energetics in liquid–liquid extraction: a simulation study of molecular solvents and ionic liquids

Extractant design in liquid–liquid extraction (LLE) is a research frontier of metal ion separations that typically focuses on the direct extractant–metal interactions. However, a more detailed understanding of energetic drivers of separations beyond primary metal coordination is often lacking, including the role of solvent in the extractant phase. In this work, we propose a new mechanism for enhancing metal-complexant energetics with nanostructured solvents. Using molecular dynamics simulations with umbrella sampling, we find that the organic solvent can reshape the energetics of the extractant's intramolecular conformational landscape. Here, we calculate free energy profiles of different conformations of a representative bidentate extractant, n-octyl(phenyl)-N,N-diisobutyl carbamoyl methyl phosphinoxide (CMPO), in four different solvents: dodecane, tributyl phosphate (TBP), and dry and wet ionic liquid (IL) 1-ethyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide ([EMIM][Tf 2 N]). By promoting reorganization of the extractant molecule into its binding conformation, our findings reveal how particular solvents can ameliorate this unfavorable step of the metal separation process. In particular, the charge alternating nanodomains formed in ILs substantially reduce the free energy penalty associated with extractant reorganization. Importantly, using alchemical free energy calculations, we find that this stabilization persists even when we explicitly include the extracted cation. These findings provide insight into the energetic drivers of metal ion separations and potentially suggest a new approach to designing effective separations using a molecular-level understanding of solvent effects.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Xe Recovery from Nuclear Power Plants Off-Gas Streams: Molecular Simulations of Gas Permeation through DD3R Zeolite Membrane

Recent experimental work has shown zeolite membrane-based separation as a promising potential technology for Kr/Xe gas mixtures due to its much lower energy requirements in comparison to cryogenic distillation, the conventional separation method for such mixtures. Such a separation is also economically rewarding because Xe is in high demand, as a valuable product for many applications/processes. In this work, we have used Molecular Dynamics (MD) simulations to study the effects of different conditions, i.e., temperature, pressure, and gas feed composition, on Kr/Xe separation performance via DD3R zeolite membranes. We provide a comprehensive study of the permeation of the different gas species, density profiles, and diffusion coefficients. Molecular simulations show that if the feed is changed from pure Kr/Xe to an equimolar mixture, the Kr/Xe separation factor increases, which agrees with experiments. In addition, when Ar is introduced as a sweep gas, the adsorption of both Kr and Xe increases, while the permeation of pure Kr increases. A similar behavior is observed with equimolar mixtures of Kr/Xe with Ar as the sweep gas. High-separation Kr/Xe selectivity is observed at 50 atm and 425 K but with low total permeation rates. Changing pressure and temperature are found to have profound effects on optimizing the separation selectivity and the permeation throughput.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Structure, stability, and superconductivity of N-doped lutetium hydrides at kbar pressures

Here, the structure of the material responsible for the room temperature and near ambient pressure superconductivity reported in an N-doped lutetium hydride has not been conclusively determined. Herein, density functional theory calculations are performed in an attempt to uncover what it might be. Guided by a range of strategies including crystal structure prediction and modifications of existing structure types, we present an array of Lu-N-H phases that are dynamically stable at experimentally relevant pressures. Although none of the structures found are thermodynamically stable, and none are expected to remain superconducting above ~17 K at 10 kbar, a number of metallic compounds with fcc Lu lattices – as suggested by the experimental X-ray diffraction measurements of the majority phase – are identified. The system whose calculated equation of states matches best with that measured for the majority phase is fluorite-type LuH 2 , whose 10 kbar superconducting critical temperature was estimated to be 0.1 K using the Allen-Dynes modified McMillan equation.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Reinforcement learning based hybrid bond-order coarse-grained interatomic potentials for exploring mesoscale aggregation in liquid–liquid mixtures

Exploring mesoscopic physical phenomena has always been a challenge for brute-force all-atom molecular dynamics simulations. Although recent advances in computing hardware have improved the accessible length scales, reaching mesoscopic timescales is still a significant bottleneck. Coarse-graining of all-atom models allows robust investigation of mesoscale physics with a reduced spatial and temporal resolution but preserves desired structural features of molecules, unlike continuum-based methods. Here, we present a hybrid bond-order coarse-grained forcefield (HyCG) for modeling mesoscale aggregation phenomena in liquid–liquid mixtures. The intuitive hybrid functional form of the potential offers interpretability to our model, unlike many machine learning based interatomic potentials. We parameterize the potential with the continuous action Monte Carlo Tree Search (cMCTS) algorithm, a reinforcement learning (RL) based global optimizing scheme, using training data from all-atom simulations. The resulting RL-HyCG correctly describes mesoscale critical fluctuations in binary liquid–liquid extraction systems. cMCTS, the RL algorithm, accurately captures the mean behavior of various geometrical properties of the molecule of interest, which were excluded from the training set. The developed potential model along with the RL-based training workflow could be applied to explore a variety of other mesoscale physical phenomena that are typically inaccessible to all-atom molecular dynamics simulations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Topological electride phase of sodium at high pressures and temperatures

Sodium, a textbook example of a nearly-free-electron metal, exhibits unforeseen pressure-induced behavior including the stabilization of numerous polymorphs—some possessing extremely complex unit cells—as well as a metal-to-insulator transition to the iconic hP4 phase. However, until recently, most of the experimental and theoretical studies on solid sodium have been restricted to the low-temperature regime. In this work, using ab initio evolutionary structure searches coupled with quasiharmonic calculations, we discover seven new phases of sodium that are more stable than the known hP4 phase at pressure-temperature conditions that were recently attained in ramp-compression experiments. From these, our calculations suggest that P6 3 /m Na is the ground state between ≈250 GPa at 710 K and ≈350 GPa at 1900 K. Electronic structure calculations show that this phase is a topological semimetal with a Dirac nodal surface that is protected by a nonsymmorphic symmetry S 2z and an electride owing to its non-nuclear charge localized within one-dimensional honeycomb channels and zero-dimensional cages. Our results highlight the complexity of dense sodium's free-energy landscape and intricate electronic structure that emerges at finite temperatures and conditions where ionic cores overlap.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Transient Gene and miRNA Expression Profile Changes of Confluent Human Fibroblast Cells in Space

Microgravity or an altered gravity environment from the static 1 gravitational constant has been shown to influence global gene expression patterns and protein levels in cultured cells. However, most of the reported studies conducted in space or using simulated microgravity on the ground have focused on the growth or differentiation of the cells. Whether non-dividing cultured cells will sense the presence of microgravity in space has not been specifically addressed. In an experiment conducted on the International Space Station, confluent human fibroblast cells were fixed after being cultured in space for 3 and 14 days for investigations of gene and miRNA (microRNA) expression profile changes in these cells. A fibroblast is a type of cell that synthesizes the extracellular matrix and collagen, the structural framework for tissues, and plays a critical role in wound healing and other functions. Results of the experiment showed that on Day 3, both the flown and ground cells were still proliferating slowly even though they were confluent, as measured by the expression of the protein Ki-67 positive cells, and the cells in space grew slightly faster. Gene and miRNA expression data indicated activation of NF(sub kappa)B (nuclear factor kappa-light-chain-enhancer of activated B cells) and other growth related pathways involving HGF and VEGF in the flown cells. On Day 14 when the cells were mostly non-dividing, the gene and miRNA expression profiles between the flight and ground samples were indistinguishable. Comparison of gene and miRNA expressions in the Day 3 samples in respect to Day 14 revealed that most of the changes observed on Day 3 were related to cell growth for both the flown and ground cells. Analysis of cytoskeleton changes by immunohistochemistry staining of the cells with antibodies for alpha-tubulin showed no difference between the flight and ground samples. Results of our study suggest that in true non-dividing human fibroblast cells, microgravity in space has little effect on the gene and miRNA expression. Gene and miRNA expression changes were observed in cells that were confluent, but still proliferating slowly. The faster growth in the flown cells was associated with the activation of NF(sub kappa)B pathways which triggers the expression of several growth factors and the suppression of the cell cycle checkpoint.

Zhang, Ye↗