Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ATOMIC DISPLACEMENTS”

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

Modeling the effect of MHD activity on runaway electron generation during SPARC disruptions

Magnetohydrodynamic (MHD) instabilities and runaway electrons (REs) interact in several ways, making it important to self-consistently model these interactions for accurate predictions of RE generation and the design of mitigation strategies, such as massive gas injection (MGI). Using M3D-C1 – an extended MHD code with a RE fluid model – we investigate the effects of 3-D nonlinear MHD activity, material injection, and 2-D axisymmetric vertical displacement events (VDEs) on RE evolution during disruptions on SPARC – a high-field, high-current tokamak designed to achieve a fusion gain Q > 1. Several cases, comprising different combinations of neon (Ne) and deuterium (D 2 ) injection, are considered. Here, our results demonstrate key effects that arise from the self-consistent RE + MHD coupling, such as an initial increase in RE generation due to MHD instability growth, decreased saturation energies of the m/n = 1/1 mode driving sawteeth-like activity, RE losses in stochastic magnetic fields, and subsequent RE confinement and plateau formation due to re-healing of flux surfaces. Large RE plateaus (>5 MA) are obtained with Ne-only injection (2-5 × 10 21 atoms), while combined D 2 + Ne injection (2 × 10 21 Ne atoms; 1.8 × 10 22 D 2 molecules) produces a lower RE current (<2 MA). With D 2 + Ne injection, a post thermal quench "cold" VDE terminates the RE beam, preventing a steady plateau. These simulations couple REs, 3-D MHD instabilities, MGI, and axisymmetric VDEs for the first time in SPARC disruption simulations and represent a crucial step in understanding RE generation and mitigation in high-current devices like SPARC.

Datta, Rishabh [Massachusetts Inst. of Technology ↗
DOE OSTI2020

Low-frequency Raman spectrum of 2D layered perovskites: Local atomistic motion or superlattice modes?

We report the low-frequency Raman spectrum (ω = 10 cm -1 –150 cm -1 ) of a wide variety of alkylammonium iodide based 2D lead halide perovskites (2D LHPs) as a function of A-site cation (MA = methylammonium and FA = formamidinium), octahedral layer thickness (n = 2–4), organic spacer chain length (butyl-, pentyl-, hexyl-), and sample temperature (T = 77 K–293 K). Using density functional theory calculations under the harmonic approximation for n = 2 BA:MAPbI, we assign several longitudinal/transverse optical phonon modes between 30 cm -1 and 100 cm -1 , the eigen displacements of which are analogous to that observed previously for octahedral twists/distortions in bulk MAPbI. Additionally, we propose an alternative assignment for low-frequency modes below this band (<30 cm -1 ) as zone-folded longitudinal acoustic phonons corresponding to the periodicity of the entire layered structure. Furthermore, we compare measured spectra to predictions of the Rytov elastic continuum model for zone-folded dispersion in layered structures. Our results are consistent across the various 2D LHPs studied herein, with energetic shifts of optical phonons corresponding to microscopic structural differences between materials and energetic shifts of acoustic phonons according to changes in the periodicity and elastic properties of the perovskite/organic subphases. This study highlights the importance of both the local atomic order and the superlattice structure on the vibrational properties of layered 2D materials.

2D materials↗
DOE OSTI2024

Teaching Nonradiative Transitions with MATLAB and Python

Nonradiative transitions are changes in energy states of atoms, ions, or molecules that do not involve the emission or absorption of photons. Despite their importance in understanding luminescent properties and photochemical reaction mechanisms, nonradiative transitions are rarely given more than a qualitative overview in undergraduate and even graduate physical chemistry curricula. To supplement the coverage of nonradiative transition topics, we provide here a set of active learning exercises to help students develop an intuitive understanding of the factors that determine the rate of nonradiative transitions. Here, we start by outlining the theoretical background through the formulations of the Franck–Condon factor and its relation to the rate of nonradiative transition. We then introduce three teaching modules, with associated MATLAB and Python codes, to explore how (1) the excited state nuclear displacement, (2) the electronic energy gap between excited and ground state and (3) the excited/ground state vibrational mode frequencies affect the magnitude of the Franck–Condon factor and thereby the rate of nonradiative transitions. The wave function overlap plots that accompany all teaching modules provide direct visualization of the effect of input parameters on the magnitude of Franck–Condon overlap integral.

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

Structure and phase transitions in niobium and tantalum derived nanoscale transition metal perovskites, Ba(Ti,MV)O3, M=Nb,Ta

The prospect of creating ferroelectric or high permittivity nanomaterials provides motivation for investigating complex transition metal oxides of the form Ba(Ti, MV)O3, where M = Nb or Ta. Solid state processing typically produces mixtures of crystalline phases, rarely beyond minimally doped Nb/Ta. Using a modified sol-gel method, we prepared single phase nanocrystals of Ba(Ti, M)O3. Compositional and elemental analysis puts the empirical formulas close to BaTi0.5Nb0.5O3−δ and BaTi0.5Ta0.5O3−δ. For both materials, a reversible temperature dependent phase transition (non-centrosymmetric to symmetric) is observed in the Raman spectrum in the region 533–583 K (260–310 °C); for Ba(Ti, Nb)O3, the onset is at 543 K (270 °C); and for Ba(Ti, Ta)O3, the onset is at 533 K (260 °C), which are comparable with 390–393 K (117–120 °C) for bulk BaTiO3. The crystal structure was resolved by examination of the powder x-ray diffraction and atomic pair distribution function (PDF) analysis of synchrotron total scattering data. It was postulated whether the structure adopted at the nanoscale was single or double perovskite. Double perovskites (A2B′B″O6) are characterized by the type and extent of cation ordering, which gives rise to higher symmetry crystal structures. PDF analysis was used to examine all likely candidate structures and to look for evidence of higher symmetry. The feasible phase space that evolves includes the ordered double perovskite structure Ba2(Ti, MV)O6 (M = Nb, Ta) Fm-3m, a disordered cubic structure, as a suitable high temperature analog, Ba(Ti, MV)O3Pm-3m, and an orthorhombic Ba(Ti, MV)O3Amm2, a room temperature structure that presents an unusually high level of lattice displacement, possibly due to octahedral tilting, and indication of a highly polarized crystal.

Chemistry↗
DOE OSTI2020

Multilayered microstructures with shape memory effects for vertical deployment

This paper presents a fabrication and characterization of multilayered microstructures with shape memory effects enabling large vertical deployment under electro-thermal actuation. Our previous research demonstrated vertical deployment of such microstructures by the effect of thermal mismatch. Development of equiatomic NiTi layers in the multilayered microstructure is investigated further for shape memory effects. Multilayered microstructures are built by sputtered NiTi layers and a lift-off process. Negative photoresist ma-N1420 enables clean lift-off of 500 nm thick NiTi layers by forming a significant undercut profile after development. The parametric study on co-sputter powers for Ti and Ni 50 Ti 50 targets suggests that 100 W RF on Ti target and 200 W DC on Ni 50 Ti 50 target can deposit Ni 49.62 Ti 50.38 layers. X-Ray Diffraction (XRD) and Atomic Force Microscopy (AFM) were used to study the crystal structures and surface topography of NiTi layers. XRD results of post-annealed Ni 49.62 Ti 50.38 layers show coexistence of austenite and martensitic phases at room temperature, suggesting that the transformation temperature of such NiTi layers should be approximate 20 °C. The surface topography of Ni 49.62 Ti 50.38 layers reveals substantial increase of surface roughness at ambient conditions after the annealing. Experimental verification of the multilayered microstructure for vertical deployment was carried out by Signatone Probe Station and Dual Scanning Electron Microscope/Focused Ion Beam (SEM/FIB) system. Finally, a vertical deployment of the two-dimensional (2D) multilayered microstructures for three-dimensional (3D) can be detected by applying a constant voltage of 0.04 V, and the expected 3D deployment displacement is enlarged from 2 μm to 10 μm by introducing the shape memory effect.

42 ENGINEERING↗
DOE OSTI2024

Sub‐Diffraction Correlation of Quantum Emitters and Local Strain Fields in Strain‐Engineered WSe 2 Monolayers

Strain-engineering in atomically thin metal dichalcogenides is a useful method for realizing single-photon emitters (SPEs) for quantum technologies. Correlating SPE position with local strain topography is challenging due to localization inaccuracies from the diffraction limit. Currently, SPEs are assumed to be positioned at the highest strained location and are typically identified by randomly screening narrow-linewidth emitters, of which only a few are spectrally pure. In this work, hyperspectral quantum emitter localization microscopy is used to locate 33 SPEs in nanoparticle-strained WSe 2 monolayers with sub-diffraction-limit resolution (≈30 nm) and correlate their positions with the underlying strain field via image registration. In this system, spectrally pure emitters are not concentrated at the highest strain location due to spectral contamination; instead, isolable SPEs are distributed away from points of peak strain with an average displacement of 240 nm. These observations point toward a need for a change in the design rules for strain-engineered SPEs and constitute a key step toward realizing next-generation quantum optical architectures.

2D material↗
DOE OSTI2025

Status of the development and testing of in-vessel and ECH-protection components for the ITER low-field side reflectometer

The ITER Low-Field Side Reflectometer (LFSR) is a critical diagnostic system designed to measure edge electron density profiles, fluctuations, and plasma rotation in ITER. This paper presents the latest developments in the design, testing, and validation of key in-vessel and Electron Cyclotron Heating (ECH) protection components. The LFSR antenna array has been optimized to provide robust coverage over expected plasma vertical displacements and a Doppler measurement for plasma rotation. Further, a comprehensive assessment of a 170-GHz diffraction grating and a novel stray-ECH power monitor demonstrates effectiveness in mitigating the impact of stray ECH power on sensitive microwave electronics. Additionally, an advancement in anti-reflective technology for millimeter waves significantly improves transmission of LFSR's vacuum window while meeting ITER's stringent safety and operational requirements. These results support the continued integration of LFSR into ITER, ensuring its diagnostic capabilities remain resilient under reactor-relevant conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Inferring relative dose-dependent color center populations in proton irradiated thoria single crystals using optical spectroscopy

In this work, we have utilized photoluminescence spectroscopy and optical ellipsometry to characterize the dose-dependence of the photoluminescence emission intensity and changes in optical absorption of thoria single crystals subject to irradiation with energetic protons at room- and elevated-temperatures. The photoluminescence peaks and the optical absorption bands are attributed to creation of new electronic states emerging from defects resulting from displacement damage. These bands are most likely associated with electrons trapped at the oxygen vacancy sites similar to color centers formed in other irradiated oxides and halides. Our experimental observations are supported by a standard density functional theory calculation of the electronic structure in pristine and oxygen vacancy-bearing thoria crystals. The dose-dependence of the intensity of the photoluminescence peaks is used to parameterize a rate theory model that estimates the concentration of color centers in the irradiated crystals. This parameterization provides optimized migration barrier parameters for oxygen interstitials and vacancies that simultaneously capture the optical response of the crystals irradiated at room- and elevated-temperature. These optical spectroscopy techniques offer a promising pathway to characterize the population of color centers formed at the sites of oxygen anion vacancies, particularly in irradiated nuclear fuels, where atomic-level defects cannot be readily imaged using electron microscopy. When combined with other direct and indirect characterization tools, our approach can provide new insight into defect formation and accumulation in energy materials over single atomic to extended length scales.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Molecular dynamics study of interstitial He clusters in nickel

This study presents a molecular dynamics analysis focusing on the behavior of interstitial helium (He) clusters in nickel (Ni), examining their formation, stability, and migration energetics. Consistent with previous research, we found that the binding energies of interstitial helium within a helium cluster are positive and increase with the cluster size, indicating a preference for helium atoms to cluster together. However, our findings also reveal that while the formation energy increases monotonically with cluster size, the increase in binding energy is non-monotonic. Importantly, small He clusters were observed to be thermally unstable at reactor operational temperatures (approximately 600 K), with the He 2 cluster exhibiting instability even at room temperature. With a binding energy of 0.49 eV for a He 4 cluster, we hypothesize that for helium bubbles to form via homogeneous nucleation (i.e., through trap mutation) at reactor operating temperatures, the helium concentration must be high enough to facilitate the formation of helium clusters of at least size 4 or larger. As expected, interstitial helium and small helium clusters are highly mobile. This mobility was observed not only at room temperature but also at temperatures as low as approximately 200 K. Furthermore, the mean squared displacement method has been utilized to determine the migration barriers and the corresponding prefactors for clusters ranging from He 1 to He 6

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Electrochemical and In Situ Spectroscopic Study of the Effect of Prussian Blue as a Mediator in a Solid-State Supercapacitor

In situ X-ray absorption spectroscopy (XAS) along with electrochemical measurements were used to explore the redox activity of Prussian blue (PB) and its contribution to pseudocapacitance in an asymmetrical solid-state supercapacitor (SC). The SC, with a configuration of (activated carbon)/polymer electrolyte membrane/(activated carbon + PB), was fabricated with an embedded KCl saturated Ag/AgCl reference electrode. The potentials of the two electrodes were generally displaced proportionally with increase in the cell voltage. Upon negative polarization from the open circuit voltage (OCV) condition to -3 V, the charge capacity was greater than that upon positive polarization from the OCV condition to +3 V, in agreement with the higher fraction of Fe(III) than that of Fe(II) and the low concentration of K ions in PB. Detailed analysis of cyclic voltammetry data in the range between 0 and -3 V indicates that the battery-like contribution to the charge current was 96% that of the capacitor-like contribution at a scanning rate of 2.5 mV s -1 but only 15% at 100 mV s -1 . The in situ XAS characterization results indicate that the majority of Fe was reversibly charged and discharged between Fe(III) and Fe(II) during positive and negative polarizations with minimal changes in local atomic structure.

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

Ab Initio Polariton Transport Dynamics with the Classical Path Approximation

We present an ab initio framework for simulating polariton transport dynamics based on the classical path approximation (CPA). The quantum dynamics of polariton transport involves simulating many electronic degrees of freedom, making a fully ab initio dynamics simulation computationally expensive. We demonstrate that the CPA, which removes the need for excited-state nuclear gradients, is well-suited for polaritonic systems because collective light–matter coupling leads to vanishing excited-state forces. Benchmark comparisons between CPA and full evaluation of the excited-state forces show excellent agreement for polariton transport results in model light–matter systems such as polariton group velocities and mean-squared displacements. Ab initio simulations of polariton transport using CPA reproduce key physical trends that are observed in experiments with BODIPY molecules. Our work establishes the CPA as a highly efficient tool for ab initio investigations of transport and energy flow in hybrid light–matter systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Real-time observation of toroidal current redistributions induced by three-dimensional MHD phenomena triggering vertical displacement events in tokamak plasmas

Three-dimensional MHD instabilities, including edge-localized modes (ELMs) and internal reconnection events (IREs), have been observed to precipitate loss of vertical stability in tokamak plasmas, resulting in vertical displacement events (VDEs). This vertical destabilization can occur due to toroidal current redistributions and/or shape changes resulting from these phenomena. Using a recently introduced method for rapidly reconstructing the two-dimensional toroidal plasma current density profile in real-time, results are presented that demonstrate the specific current distribution changes that occur during ELMs (on KSTAR) and IREs (on MAST-U) that lead to loss of vertical control. The method most efficiently reconstructs the toroidal current density profile by doing so on a basis of principal components of historical profiles. These principal components isolate dominant current profile dynamics, improving interpretability, increasing speed, and reducing dimensionality of the profile computation. On KSTAR, this computation is executed in the real-time plasma control system at a rate of 10 kHz (limited by available CPU cycle times), allowing the current profile evolution to be assessed at several times over the course of each ELM event. Further, by incorporating the reconstructions into a novel vertical stability metric, the contribution of specific current profile dynamics to the loss of vertical stability can be assessed in real-time for VDE avoidance and improved understanding of the causal relationship between three-dimensional MHD phenomena and VDEs. The success of this method in approximating toroidal current density profiles from kinetic equilibrium reconstructions is also presented ($R^2=0.990$), along with its capability to produce other equilibrium quantities of interest in real-time at high time resolution.

edge-localized modes↗
DOE OSTI2020

Viscoelastic Response of Dispersed Entangled Polymer Melts

Any polymer synthesis routes result in molecular weight dispersity DAI, depending on the polymerization mechanism. The lowest dispersity polymers are those made by anionic and atom-transfer radical polymerization, which exhibit relatively narrow distributions DM=g,„1A172 — 1.02-1.04. This small divergence from monodispersity results in significant number of molecules that differ in their molecular weight from the average and their impact on the viscoelastic response remains an open question. Here the effects of low dispersity on stress relaxation and shear viscosity of entangled melts are studied using a coarse-grained model for polyethylene. Polymer melts with chain length dispersity set to follow a Schulz-Zimm distribution of Dm = 1.0 — 1.16 for an entangled polyethylene melt using molecular dynamic simulations. The systems were studied to times up to 800 ,us which is beyond the terminal time. These systems are compared to those with binary and ternary distributions. The stress relaxation functions were extracted from the Green-Kubo relation and from stress-strain relaxation following a uniaxial extension. We find on the entanglement time scale, both the mean squared displacement and the stress relaxation are independent of D m . The tube diameter and the rubbery plateau for entangled dynamics do not depend on D m . At longer times, the terminal relaxation time decreases as the faster motion of the shorter chains results in constraint release for the longer chains. Probing this low dispersity regime opens the way to evaluate the degree of dispersity that affects mechanical properties.

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

Phase segregation and miscibility of TiO x nanocomposites in Gd-doped ceria solid electrolyte material

Electro-chemo-mechanical (ECM) coupling refers to mechanical deformation due to electrochemically driven compositional change in a solid. An ECM actuator producing micrometre-size displacements and long-term stability at room temperature was recently reported, comprising a 20 mol% Gd-doped ceria (20GDC), a solid electrolyte membrane, placed between two working bodies made of TiO x /20GDC (Ti-GDC) nanocomposites with Ti concentration of 38 mol%. The volumetric changes originating from oxidation or reduction in the local TiO x units are hypothesized to be the origin of mechanical deformation in the ECM actuator. Studying the Ti concentration-dependent structural changes in the Ti-GDC nanocomposites is therefore required for (i) understanding the mechanism of dimensional changes in the ECM actuator and (ii) maximizing the ECM response. Here, the systematic investigation of the local structure of the Ti and Ce ions in Ti-GDC over a broad range of Ti concentrations using synchrotron X-ray absorption spectroscopy and X-ray diffraction is reported. The main finding is that, depending on the Ti concentration, Ti atoms either form a cerium titanate or segregate into a TiO 2 anatase-like phase. The transition region between these two regimes with Ti(IV) concentration between 19% and 57% contained strongly disordered TiO x units dispersed in 20GDC containing Ce(III) and Ce(IV) and hence rich with oxygen vacancies. As a result, this transition region is proposed to be the most advantageous for developing ECM-active materials.

47 OTHER INSTRUMENTATION↗
DOE OSTI2021

Nearly hyperuniform, nonhyperuniform, and antihyperuniform density fluctuations in two-dimensional transition metal dichalcogenides with defects

Hyperuniform many-body systems in d-dimensional Euclidean space R d are characterized by completely suppressed (normalized) infinite-wavelength density fluctuations, and appear to be endowed with novel exotic physical properties. Recently, hyperuniform systems of disordered varieties have been observed in the context of various atomic-scale two-dimensional (2D) materials. Here, we analyze the effects of localized defects on the density fluctuations across length scales and on the hyperuniformity property of experimental samples of two-dimensional transition metal dichalcogenides. In particular, we extract atomic coordinates from time series annular dark field-scanning transmission electron microscopy (ADF-STEM) imaging data of 2D tungsten chalcogenides with the 2H structure (Te-doped 2H-WSe 2 ) showing continuous development and evolution of electron-beam induced defects, and construct the corresponding chemical-bonding informed coordination networks between the atoms. We then compute a variety of pair statistics and bond-orientational statistics to characterize the samples. At low defect concentrations, the corresponding materials are nearly hyperuniform, characterized by significantly suppressed scattering at the zero wave-number limit (omitting forward/ballistic scattering). As more defects are introduced, the (approximate) hyperuniformity of the materials is gradually destroyed, and the system becomes non-hyperuniform even when the material still contains a significant amount of crystalline regions. At high defect concentrations, the structures become antihyperuniform with diverging (normalized) large-scale density fluctuations, mimicking those typically observed at the thermal critical points associated with phase transitions. Overall, the defected materials possess varying degrees of orientation order, and there is apparently no intermediate hexatic phase emerging. To understand the observed nearly hyperuniform density fluctuations in the slightly defected materials, we construct a minimalist structural model and demonstrate that the experimental samples can be essentially viewed as perturbed honeycomb crystals with small correlated displacements and double chalcogen vacancies. Moreover, the small correlated displacements alone can significantly degrade hyperuniformity of the perfect honeycomb structure. Therefore, even a small amount of vacancies, when coupled with correlated displacements, can completely destroy hyperuniformity of the system.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2023

Prediction and Control of Atomic Ordering in Electrodeposited Binary Alloy Films: Direct Synthesis of L10 Magnetic Phases

Permanent magnet materials are essential components in renewable energy technology, underlying, electric vehicles, hybrids, wind turbines and more. Today commercial permanent magnets based on Nd-Fe-B exhibit an energy product of ~ 400 kJ/m 3 , capable to run wind turbines efficiently. Unfortunately, the supplies of rare-earth elements are centralized, and their availability may see volatility over time. This state of affairs stimulates the search for the search for rare-earth-free permanent magnetic materials. The L1 0 crystal structure underlines an important class of chemical ordered alloys that exhibit uniaxial magnetocrystalline anisotropy. This material has been considered an excellent candidate for rare-earth-free permanent magnets. In this program, a series of systematic experiments were performed and analysed in order to understand the effect of each. In order to synthesize L 10-Fe-Ni, one has to circumvent the sluggish kinetic and the limited driving force resulted from the low order-disorder temperature (~320°C) of the system. A strategy was found to trigger a martensitic transformation from the FCC phase to one of the non-cubic phases by high-strain methods. The body-centered tetragonal Fe-Ni obtained by this strategy was considered as a precursor of L10 Fe-Ni. The formation of BCT Fe-Ni via a displacive process was also implemented at the nanoscale by using tailored AuNi@FeNi core@shell.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Intramolecular 1,2 C-H Addition of o -Methyl Groups to Form Unique Ruthenium Pincer Tuck-in Complexes

New RPN H P ligands containing 2,4-xylyl (4mXPN H P) and mesityl (MesPN H P) groups on the phosphorus atoms were synthesized. 4mXPN H P reacts with [(cymene)RuCl 2 ] 2 followed by PMe 3 to produce κ 3 -4mXPN H PRu(PMe 3 )Cl 2 . MesPN H P reacts with [(cymene)RuCl 2 ] 2 to produce monomeric κ 3 -MesPN H PRuCl 2 that reacts with CO forming κ 3 -MesPN H PRu(CO)Cl 2 . Surprisingly, dehydrohalogenation of these complexes results in the activation of ortho methyl groups of the pincer ligands, rather than formation of κ 3 -RPNPRuLCl complexes. This results from transient κ 3 -RPNPRuLCl formation followed by 1,2-addition of an ortho C-H bond across the Ru-amide bond. The transient amide complex of κ 4 -4mXPN H PRu(PMe 3 )Cl was trapped with CO forming κ 3 -4mXPNPRu(PMe 3 )(CO)Cl. In contrast, κ 4 -MesPNHPRu(CO)Cl does not react with ligands to trap the expected amide complex of reverse C-H addition. Instead, CO and PMe 3 displace the chloride ligand forming cationic complexes. In both cases, hemi-lability of the pincer ligand was observed spectroscopically. In conclusion, the new complexes serve as precursors to moderately active catalysts for the acceptorless dehydrogenative coupling of n-butanol.

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

Effect of Solvents on Lignin–Surface Interactions via Molecular Dynamics Simulations

Lignin, an essential building block of lignocellulosic biomass, is a potential abundant source of aromatic monomers for the polymer and chemical industry. Reductive catalytic fractionation (RCF) is one promising process that can produce high yields of phenolic monomers and oligomers from lignin under different catalytic conditions. An important choice in optimizing RCF is the selection of solvent; however, detailed insights into solvent effects on lignin behaviors and interactions remain limited. Here, in this work, we perform all-atom molecular dynamics simulations to study the solvation of lignin, solvent-mediated conformational changes, and the interaction of solvated lignin oligomers with model surfaces. We focus on the behavior of an oligomeric lignin model compound in methanol, ethanol, a binary mixture of ethanol and water, and water at both the RCF reaction temperature (473 K) and room temperature. Analysis of structural features of lignin suggests that these three organic solvent systems favorably solvate lignin, resulting in a more extended conformation suitable for catalytic conversion to valuable chemicals. We further introduce model palladium (Pd) and carbon (C) surfaces to understand how solvent choice impacts adsorption onto a representative catalytic surface and support, and to quantify the competition among the reactant and solvent molecules for the surface. Unbiased simulations suggest strong adsorption of lignin on both Pd and C surfaces at 473 K, with notable solvent-mediated differences in adsorption energies. Additionally, our findings indicate that lignin adsorption is promoted by the entropy change resulting from the displacement of solvent molecules from the surface. This study provides a molecular perspective of adsorption of lignin onto varying surfaces, which is a step towards understanding and optimizing the catalytic conversion of lignin into valuable chemicals.

adsorption↗