Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ILs”

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 1,549 records · Page 86
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Covalent Adaptable Networks: Reprocessable Cross-Linked Polymers

Thermoset polymers have desirable properties, such as excellent thermal and mechanical stability, but their covalent cross-links typically prevent repair or recycling. By enabling and controlling dynamic exchange reactions within polymer networks, their covalent bonds rearrange and allow the polymer to be reshaped. These viscoelastic polymer networks, now known as covalent adaptable networks (CANs), are an important frontier for improving plastic circularity, as well as for designing valuable stimuli-responsive materials. This Review describes the history of CANs, dating back to the early days of polymer science, and the evolution of their classification and nomenclature. A comprehensive survey of dynamic reactions and linkage chemistries is provided, as well as methods to characterize and reprocess CANs. Beyond straightforward reprocessing, many advanced applications of CANs and their composites are now emerging. Lastly, we provide perspective on how the development of new chemistries, strategies to control stimuli-responsive bond exchange and mechanical properties, and a deep understanding of exchange reactions will advance this field toward scalable, sustainable, and high-value materials.

Sun, Molly [Northwestern University, Evanston, IL ↗
DOE OSTI2025

Alizarin Multilayers Adsorbed onto Glassy Carbon Electrodes for Electrochemical Sequestration of Manganese, Sodium, and Lithium Cations

The development of efficient methods for metal ion recovery and water remediation is critical for addressing the nation’s urgent need for secure domestic supply chains, overcoming critical materials challenges, and ensuring resilient manufacturing. Here, we present a novel approach for synthesizing stable Alizarin (Alz, 1,2-dihydroxy anthraquinone) multilayers on glassy carbon electrodes (GCEs) for electrochemical sequestration of manganese, sodium, and lithium cations. Alz-GCEs exhibit negatively shifted reduction potentials, indicating strong interactions between Alz and metal cations through metal-coupled electron transfer (MCET) mechanisms. The cation binding interactions and redox behavior of these electrodes were investigated using cyclic voltammetry (CV) and density functional theory (DFT) calculations. Our results demonstrate that Alz forms multilayer structures on GCEs with redox properties that are modulated by the presence of metal cations in the electrolyte. DFT calculations provide insights into the electrochemical mechanism, indicating both stepwise and concerted pathways for metal binding. The findings highlight the potential of Alz-GCEs for efficient and selective metal ion capture, which could be useful for developing sustainable materials for critical metal recovery and water remediation. This work suggests that redox-mediated organic adlayers are promising candidates for advancing electrochemical separation technologies of metal ions.

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

GPU-Accelerated Solution of the Bethe–Salpeter Equation for Large and Heterogeneous Systems

We present a massively parallel GPU-accelerated implementation of the Bethe–Salpeter equation (BSE) for the calculation of the vertical excitation energies (VEEs) and optical absorption spectra of condensed and molecular systems, starting from single-particle eigenvalues and eigenvectors obtained with density functional theory. The algorithms adopted here circumvent the slowly converging sums over empty and occupied states and the inversion of large dielectric matrices through a density matrix perturbation theory approach and a low-rank decomposition of the screened Coulomb interaction, respectively. Further computational savings are achieved by exploiting the nearsightedness of the density matrix of semiconductors and insulators to reduce the number of screened Coulomb integrals. We scale our calculations to thousands of GPUs with a hierarchical loop and data distribution strategy. The efficacy of our method is demonstrated by computing the VEEs of several spin defects in wide-band-gap materials, showing that supercells with up to 1000 atoms are necessary to obtain converged results. We discuss the validity of the common approximation that solves the BSE with truncated sums over empty and occupied states. In conclusion, we then apply our GW-BSE implementation to a diamond lattice with 1727 atoms to study the symmetry breaking of triplet states caused by the interaction of a point defect with an extended line defect.

Absorption spectra↗
DOE OSTI2025

Polynomial Scaling Localized Active Space Unitary Selective Coupled Cluster Singles and Doubles

We present a polynomial-scaling algorithm for the localized active space unitary selective coupled cluster singles and doubles (LAS-USCCSD) method. In this approach, cluster excitations are selected based on a threshold ϵ determined by the absolute gradients of the LAS-UCCSD energy with respect to cluster amplitudes. Using the generalized Wick’s theorem for multireference wave functions, we derive the gradient expression as a polynomial function of one-, two-, and three-body reduced density matrices and 1- and 2-electron integrals, valid for any multireference wave function. The resulting gradient implementation exhibits a memory scaling of 𝒪(N 6 ), with N spin orbitals in the combined active space of all fragments. The variational quantum eigensolver is used to optimize the selected cluster excitations on a quantum simulator. Furthermore, by plotting the energy error, defined as the difference between the LAS-USCCSD and corresponding CASCI energies, against the inverse cluster amplitude selection threshold (ϵ –1 ) for polyene chains containing 2 to 5 π-bond units, we establish a relationship between the energy error and the threshold. To further validate the accuracy of LAS-USCCSD, we computed the cis–trans isomerization energy of stilbene (a 20-qubit system) and the magnetic coupling constant of the tris-hydroxo-bridged chromium dimer [Cr 2 (OH) 3 (NH 3 ) 6 ] 3+ (evaluated as both 12- and 20-qubit systems) using the Qiskit-Qulacs simulator. Assessing such examples is important to determine the practical feasibility of quantum simulations for chemically realistic systems. Toward this goal, with the LAS-USCCSD algorithm we estimated the quantum resources required for simulating an active space of (30e,22o) in [Cr 2 (OH) 3 (NH 3 ) 6 ] 3+ , a size that remains beyond the reach of current quantum simulators for accurate treatment.

Algorithms↗
DOE OSTI2026

Characterizing Defect Dynamics in Silicon Carbide Using Symmetry-Adapted Collective Variables and Machine Learning Interatomic Potentials

Silicon carbide (SiC) divacancies are attractive candidates for spin-defect qubits possessing long coherence times and optical addressability. The high activation barriers associated with SiC defect formation and motion pose challenges for their study by first-principles molecular dynamics. In this work, we develop and deploy machine learning interatomic potentials (MLIPs) to accelerate defect dynamics simulations while retaining ab initio accuracy. We employ an active learning strategy comprising symmetry-adapted collective variable discovery and enhanced sampling to compile configurationally diverse training data, calculation of energies and forces using density functional theory (DFT), and training of an E(3)-equivariant MLIP based on the Allegro model. Here, the trained MLIP reproduces DFT-level accuracy in defect transition activation free energy barriers, enables the efficient and stable simulation of multidefect 216-atom supercells, and permits an analysis of the temperature dependence of defect thermodynamic stability and formation/annihilation kinetics to propose an optimal annealing temperature to maximally stabilize VV divacancies.

Computer simulations↗
DOE OSTI2026

Encumbering the Copper(I)-Diimine Coordination Sphere with a Benzyl Group: Impact on the Stability and Photoredox Properties

Copper(I)-bisdiimine complexes, [Cu(NN) 2 ] + (where NN represents a diimine such as phenanthroline or bipyridine derivatives), are considered promising photosensitizers for various photochemical applications. However, their effectiveness is subject to several key challenges. In particular, controlling steric strain around the copper(I) center by introducing appropriate substituents (R) at the α-position of the nitrogen atoms is crucial for optimizing the excited-state properties of the complex. In brief, increasing the size of R leads to longer emission lifetimes and higher quantum yields. Additionally, the energy of the singlet excited state rises with increasing steric bulk, enhancing photoinduced reactivity. However, excessive steric strain from bulky substituents can significantly destabilize the coordination sphere. To balance complex stability with increased steric bulk around the metal center, we have developed two novel non-symmetrical ligands featuring branched alkyl chains and benzyl groups at the α-position of the nitrogen atoms. Here, our findings demonstrate that intramolecular π-stacking interactions between the benzyl group and the opposing phenanthroline ligand contribute to stabilizing the coordination sphere. Furthermore, the flexibility of the benzyl group reinforces the tetrahedral geometry around copper(I), resulting in an increased singlet excitedstate energy compared to benchmark complexes. Notably, we show that this enhancement in excited-state energy translates into greater excited-state reactivity.

Copper(I)↗
DOE OSTI2025

Coupling of Charge Regulation and Geometry in Soft Ionizable Molecular Assemblies

The size, shape, and charge of structures, such as proteins and amphiphile assemblies, respond in an interconnected manner to solution ionic conditions. Here, we analyze assemblies of an amphiphile (C 16 K 2 ), with two ionizable amino acids [lysine (K)] coupled to a 16-carbon alkyl tail, via small-angle X-ray scattering (SAXS), nonlinear Poisson–Boltzmann theory (nl-PB), and hybrid Monte Carlo-molecular dynamics (MC-MD) simulations. SAXS revealed structural transitions from spherical micelles to cylindrical micelles to bilayers with increasing pH. By combining SAXS-determined structural information and nl-PB, we derived the molecular degree of ionization as a function of pH. The back-calculated titration curves matched the experimental data over an extended pH range, without adjustable parameters. Similarly, the SAXS data on the evolution of spherical micelle structure with ionic strength were combined with nl-PB and MC-MD to derive the bare and effective charges. MC-MD, which considered finite ion sizes, showed that bare and effective charges saturate quickly with increasing salt concentration. Furthermore, the calculated effective charges closely matched results from Zeta-potential measurements. The presented approach has advantages over customary methods for charge regulation, such as the Henderson–Hasselbalch (HH) or Hill models, where molecular ionization/deionization in assemblies is described by effective pKs that are distinct from the pK for isolated molecules. However, these models lack a physical explanation for these pK shifts. By contrast, our approach of combining structural details with an electrostatic model and simulations provides a more intuitive understanding of structure-charge coupling and a framework for understanding charge regulation in many synthetic and biological systems.

ionization↗
DOE OSTI2025

Computational Study of Indium Oxide Photoelectrodes

Using a combination of first-principles molecular dynamics simulations and electronic structure calculations, we characterize the atomistic structure and vibrational properties of a photocatalytic surface of In 2 O 3 , a promising photoelectrode for the production of hydrogen peroxide. We then investigate the surface in contact with water and show that the electronic states of In 2 O 3 are appropriately positioned in energy to facilitate the two-electron water oxidation reaction (WOR) over the competing four-electron oxygen evolution reaction. In conclusion, we further propose that the use of strained thin films interfaced with water is beneficial in decreasing the optical gap of In 2 O 3 and thus in utilizing a wider portion of the solar spectrum for the WOR.

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

Spectra-Orthogonal Optical Anisotropy in Wafer-Scale Molecular Crystal Monolayers

Controlling the spectral and polarization response of two-dimensional (2D) crystals is vital for developing ultrathin platforms for compact optoelectronic devices. However, independently tuning optical anisotropy and spectral response remains challenging in conventional semiconductors due to the intertwined nature of their lattice and electronic structures. Here, we report spectra-orthogonal optical anisotropy─where polarization anisotropy is tuned independently of spectral response─in wafer-scale, one-atom-thick 2D molecular crystal (2DMC) monolayers synthesized on monolayer transition-metal dichalcogenide (TMD) crystals. Utilizing the concomitant spectral consistency and structural tunability of perylene derivatives, we demonstrate tunable optical polarization anisotropy in 2DMCs with similar spectral profiles, as confirmed by room-temperature scanning tunneling microscopy and cross-polarized reflectance microscopy. Additional angle-dependent analysis of the single-crystal and polycrystalline molecular domains reveals an epitaxial relationship between the 2DMC and TMD. In conclusion, our results establish a scalable, molecule-based 2D crystalline platform for unique and tunable functionalities unattainable in covalent 2D solids.

2D materials↗
DOE OSTI2025

Supported Organoiridium-Pincer Catalysts for the Nonoxidative Dehydrogenation of High-Density Polyethylene

The iridium-pincer complex {p-OP( t Bu) 2 -C 6 H 2 -2,6-OP( t Bu) 2 ] 2 }Ir(C 2 H 4 ) ( P [Ir]) has been reported as a stable and active catalyst toward alkane dehydrogenation in homogeneous and supported heterogeneous systems. Dehydrogenation has been shown as a practical method toward functional polyolefins, with dehydrogenated high-density polyethylene (deHDPE) demonstrated as a valuable synthon for upcycling, as orthogonal C-H strategies are key to end-of-life upcycling. The heterogenization of P [Ir] on oxides (SiO 2 , Al 2 O 3 , and TiO 2 ; P [Ir]/E y O x ) yields a mixture of organometallic Ir-fragments whose catalytic nonoxidative dehydrogenation activity is modulated by the binding modes of the active metal on the surface. The binding mode was elucidated by a combination of solid-state NMR and XAFS analyses and supported by DFT calculations. Surface binding through the ligand enables active organoiridium that catalyzes internal olefination of deHDPE up to 1.23 mol % at 200 °C under dynamic vacuum. Alternatively, when the organoiridium is bonded though the metal center (Ir-O SiO ), catalyst activity is negligible. Furthermore, the catalytic activity of P [Ir]/SiO 2 showed comparable reactivity with the homogeneous analogue under the same catalytic conditions, and the heterogenized catalyst can be reused up to three cycles. In conclusion, this work highlights the importance of understanding how organometallic precursors react with hydroxylated metal oxide surfaces to establish structure-property relationships.

Hunt, Samuel B. [Argonne National Laboratory (ANL)↗
DOE OSTI2025

Minimization of Cathode|Solid-Electrolyte Interfacial Delamination through the Application of Interphase Layers

Next-generation lithium-ion batteries are expected to use solid electrolytes (SEs) to enable higher energy density and extreme fast-charge capabilities. One major mode of degradation at the cathode|SE interface is delamination between the cathode active materials and SEs, which leads to performance decay. Experimental observations indicate that implementation of interphase layers can minimize the cathode|SE delamination induced capacity fade. A multiscale computational methodology is developed here to investigate the applicability of boron substituted lithium carbonate (Li 2+x B x C 1–x O 3 , x = 0.5, or LBCO) to minimize the delamination at the cathode|SE interface. Atomistic simulations indicate that the fracture energies at both the cathode|LBCO and LBCO|SE interfaces are higher than those at the cathode|SE interface, which reduces the extent of delamination. Mesoscale simulations indicate that, apart from increasing the fracture energy, decreasing the evolution of strain energy by lowering the elastic modulus of the interphase layer can also minimize the extent of delamination at the cathode|SE interface. However, the adoption of an interphase layer with high ionic conductivity is necessary to minimize the ohmic losses during operation at higher current densities. This study provides guidance on selecting interphase layers with specific properties and thicknesses to minimize both interfacial delamination and impedance growth.

LBCO↗
DOE OSTI2025

Molecular Level Understanding of Polyethylene Terephthalate (PET) Depolymerization in Base/Alcohol Hybrid Systems

Polyethylene terephthalate (PET) depolymerization in base/alcohol hybrid systems represents a promising low-energy approach for chemically recycling PET waste into valuable monomers. This study investigates the mechanistic pathways of PET depolymerization in NaOH/alcohol solutions, emphasizing the competing roles of hydroxide and alkoxide species. Utilizing a combination of experimental techniques, density functional theory (DFT) calculations, and molecular dynamics (MD) simulations, we explore how factors such as base concentration, alcohol chain length, and pKa values of alcohols influence PET depolymerization efficiency and pathways. Our findings indicate that alkoxide ions (RO⁻) exhibit notably higher reactivity than hydroxide ions (HO⁻), favoring an alcoholysis pathway in the base/alcohol hybrid system. Experimental results across a series of C1 to C5 alcohols show that longer-chain alcohols, particularly 1-butanol, achieve higher PET conversion, although this does not align solely with simple nucleophilicity trends of alkoxides. While DFT calculations reveal comparable activation energies for various alkoxides in PET depolymerization, MD simulations underscore the significant role of alcohol chain length, with longer-chain alcohols forming more stable or frequent interactions with PET. Additionally, the alkoxide concentration, influenced by the alcohol’s pK a , directly impacts PET conversion. These suggest that PET depolymerization is governed by a balance between alkoxide concentration and alkoxide-PET interactions, rather than activation energies or nucleophilicity alone. From a practical perspective, incorporating long-chain alcohols as cosolvents may enhance process efficiency but increases raw material costs by approximately 30%. However, long-chain alcohols present a safer and more sustainable alternative to hazardous cosolvents such as dichloromethane. This work offers a molecular-level understanding of PET depolymerization in base/alcohol systems and provides insights into optimizing these systems for more efficient and sustainable PET recycling processes.

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

Boosting the Oxygen Evolution Reaction by Tuning the Interfacial Iron Adsorption on Layered Double Hydroxide

Understanding the interaction between ions in the electrolyte and electrode materials plays an important role in optimizing the water electrolysis performance for hydrogen production. Herein, the synergistic effect of iron (Fe) in the electrolyte and interlayer anions within the layered structure on the oxygen evolution reaction (OER) has been investigated by combining material synthesis with controlled structure, multiple characterization techniques, and first-principles calculations. Nickel aluminum layered double hydroxides (NiAl-LDHs) with different interlayer anions (CO 3 2– , Cl – , and Br – ) show similar oxygen evolution activity in the absence of Fe species in the electrolyte. The addition of Fe into the electrolyte results in improved performance for all of the NiAl-LDHs, following the rank LDH-Br > LDH-Cl > LDH-CO 3 , under all of the conditions with varied concentration of Fe. X-ray absorption spectroscopy and identical location electron microscopy analyses show that the LDH structure remains unchanged after the OER activity test, while in situ stationary probe rotating disk electrode inductively coupled plasma mass spectrometry (SPRDE-ICP-MS) measurements show partial dissolution of the intercalating halide ions during cycling, with less dissolution for Br-intercalated materials. Insights from theoretical calculations demonstrate the thermodynamic preference of Br – to remain intercalated in the presence of Fe, while the stronger adsorption of Fe(OH) 3 species on the LDH-Br sample promotes the OER activity. In conclusion, these results provide mechanistic insights into the rational design of active layered materials with an enhanced OER performance for efficient water electrolysis.

58 GEOSCIENCES↗
DOE OSTI2026

“Antenna-like” Light-Harvesting in Partially Metalated Porphyrin Hydrogen-Bonded Organic Frameworks

Arrays of light harvesting molecules with long-range order hold promise as materials capable of efficient collection of photons as well as rapid transport of captured solar energy to chemical catalysts. Hydrogen-bonded organic frameworks (HOFs) could be promising artificial light-harvesting platforms, due to their crystalline nature and their comparative ease of synthesis. In this work, mixed-linkers of either zinc-metallated or free base tetrakis(4-carboxyphenyl)porphyrin HOFs were examined and found to demonstrate ‘antenna-like’ light harvesting behavior. Using time-resolved emission spectroscopy, Förster resonance energy transfer from zinc porphyrin to free-base porphyrin was assessed. Based on changes in the emission peak profile as the ratio of zinc porphyrin to free-base porphyrin is varied, we find that, in the limit of 100% zinc porphyrin, a photogenerated exciton can sample ~48 chromophores. These findings suggest that stacked linkers within HOFs can act as light harvesting ‘antennae’ and with encapsulation of suitable catalysts, may offer prospective application as photochemical energy conversion systems.

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

Direct Determination of 3C/4H Silicon Carbide Heterophase Interfaces by Electron Ptychography

Silicon carbide is a technologically relevant material that can exist in various structural polymorphs, each with its own electronic properties. Through the controlled growth of abrupt interfaces between the polymorphs, systems with emergent properties arise from differences in local bonding and stacking sequences. Determining the structure of these interfaces provides important insights into predicting and ultimately controlling their properties. Here, in this study, we use multislice electron ptychography (MEP) to qualitatively and quantitatively analyze interfaces between 3C and 4H silicon carbide in three dimensions. Two distinct interfaces are investigated: a coherent interface between $(1\bar{1}1)$ (3C) and $(0001)$ (4H) , and an incoherent interface between $(1\bar{1}2)$ (3C) and $(1\bar{1}00)$ (4H) . At the coherent interface, MEP enables direct quantification of a sharp boundary with an inclined step, whereas at the incoherent interface, local atomic displacements are determined along dislocation cores. By quantifying the spatial distribution and atomic-scale structure of these interfaces in three dimensions, this study provides insight into the complexity of silicon carbide heteropolytype interfaces, paving the way for theoretical modeling and the development of design rules to achieve desired properties.

heterophase interface↗
DOE OSTI2025

Free-Energy Landscapes and Surface Dynamics in Methane Activation on Ni(511) via Machine Learning and Enhanced Sampling

Methane activation on stepped Ni(511) surfaces involves the rearrangement of surface atoms as the chemical reaction proceeds. This process is particularly sensitive to temperature. Using machine-learned interatomic potentials (MLIPs) coupled with enhanced sampling techniques, we investigate the activation of methane under realistic operando conditions. Our analysis reveals that methane dissociation occurs predominantly at step-edge nickel atoms. As CH x (where x = 3 or 4) species bind to additional surface nickel atoms, their reduced mobility leads to entropic penalties that suppress certain configurations and transition states. This is reflected in the underlying free energy surfaces, where configurations such as methyl binding to hollow sites and activation routes involving two nickel atoms become unfavorable as temperature increases. At elevated temperatures, methane activation extends from step-edge sites to terrace regions because of reduced free-energy barriers and enhanced surface dynamics. By decomposing the free-energy into enthalpic and entropic contributions, we uncover temperature-dependent shifts in the preferences of methane for the relevant active sites and arrive at a detailed molecular picture of methane activation.

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

Role of the Mobile Active Site Flap in IMP Dehydrogenase Inhibitor Binding

Inosine 5′-monophosphate dehydrogenase (IMPDH) is a promising antibiotic target. This enzyme catalyzes the NAD-dependent oxidation of inosine 5′-monophosphate (IMP) to xanthosine 5′-monophosphate (XMP), which is the rate-limiting step in guanine nucleotide biosynthesis. Bacterial IMPDH-specific inhibitors have been developed that bind to the NAD + site. These inhibitors display varied affinities to different bacterial IMPDHs that are not easily rationalized by X-ray crystal structures of enzyme–inhibitor complexes. Inspection of X-ray crystal structures of 25 enzyme–inhibitor complexes, including 10 newly described, suggested that a mobile active site flap may be a structural determinant of inhibitor potency. Saturation transfer difference NMR experiments also suggested that the flap may contact the inhibitors to varying extents in different IMPDHs. Flap residue Leu413 contacted some inhibitors but was not structured in the crystal structures of other inhibitor complexes. The substitution of Leu413 with Phe or Ala in Bacillus anthracis IMPDH had inhibitor-selective effects, suggesting residue 413 could be a structural determinant of affinity. Curiously, the Ala substitution increased the potency of most inhibitors, even those that contacted Leu413 in the crystal structures. Presteady-state and steady-state kinetics experiments showed that the Leu413Ala substitution had comparable effects on inhibitor binding to the noncovalent E·IMP complex and the covalent intermediate E-XMP*, suggesting that the flap had similar interactions in both complexes. These results demonstrate that contacts do not necessarily indicate favorable interactions, and poorly structured mobile regions should not be discounted when assessing binding determinants.

Conformation↗
DOE OSTI2025

Advanced Polymer Characterization: Modular Operations for Spectral Alignment by Iterative Compression (MOSAIC)

Matrix-assisted laser desorption/ionization (MALDI) mass spectrometry encodes structural information across diverse homo- and copolymer ensembles, yet decrypting these spectra requires a systematic analytical approach. We introduce Modular Operations for Spectral Alignment by Iterative Compression (MOSAIC)─a general cipher algorithm that applies modular arithmetic to filter monomer-derived mass contributions and cluster MALDI peaks by nonconstitutional repeating units (non-CRUs). MOSAIC performs sequential modular operations using monomer mass differences as base units to compress complex spectral data, revealing end-group distributions and comonomer incorporation. As a demonstration, we applied MOSAIC to five copolymers formed by two different polymerization mechanisms. Furthermore, the resulting remainder–mass plots clearly resolve polymer homologs with distinct non-CRUs into visually apparent clusters, enabling intuitive assignment of mass spectral features.

Wang, Hanlin M. [University of Illinois at Urbana−↗