Search NASA⌕ Search

Engineering topics

Li, Wan-Lu

Publications and source records attributed to Li, Wan-Lu.

11 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Mechanistic origin of solvent-dependent thermal stability in sodiated Sn anodes for sodium-ion batteries

Understanding the thermal stability of high-energy density alloy anodes is critical for the safe deployment of sodium-ion batteries (SIBs). Here, accelerating rate calorimetry (ARC), post-mortem characterizations, and density functional theory (DFT) calculations are combined to understand the thermal reactivity of fully sodiated Sn, Sn-hard carbon (HC) blends, and HC anodes in carbonate- and ether-based electrolytes. ARC measurements show that propylene carbonate (PC) causes earlier self-heating rate (SHR) onset and higher reactivity than tetraethylene glycol dimethyl ether (TEGDME), indicating inferior thermal stability. Sodiated Sn exhibits better thermal stability than sodiated HC, while Sn-HC blends show intermediate behavior that improves with increasing Sn content. Post-ARC analyses reveal desodiation of Na15Sn4 to metallic Sn with particle coalescence, whereas Sn-HC blends and HC retain finer morphologies. PC promotes Sn oxidation to SnO, while TEGDME suppresses oxide formation; NaPF6-containing electrolytes additionally form NaF. DFT calculations show that PC adsorption lowers Na extraction energy and enhances interfacial electronic interactions, facilitating Na release and reductive decomposition. These results establish a direct correlation between solvent-dependent reaction pathways and thermal stability in SIB alloy anodes.

Accelerating rate calorimetry↗
DOE OSTI2023

How thermal fluctuations influence the function of the FeMo cofactor in nitrogenase enzymes

The catalytic mechanism of N 2 fixation by nitrogenase remains unresolved in how the strong N≡N bond is activated and why the reductive elimination of H 2 is required. Here, we use density functional theory and physiologically relevant thermal simulations to elucidate the mechanism of the complete nitrogenase catalytic cycle. Over the accumulation of four reducing equivalents, we find that protons and electrons transfer to the FeMo cofactor to weaken and break its bridge Fe–S bond, leading to temporary H 2 S formation that exposes the Fe sites to weakly bind N 2 . Remarkably, we find that subsequent H 2 formation is responsible for chemical activation to an N=N double bond accompanied by a low barrier for H 2 release. We emphasize that finite temperature effects smooth out mechanistic differences between DFT functionals observed at 0 K, thus leading to a consistent understanding as to why H formation is an obligatory step in N 2 adsorption and activation.

DFT↗
DOE OSTI2022

Optimizing the Solvent Reorganization Free Energy by Metal Substitution for Nanocage Catalysis

The supramolecular capsule Ga 4 L 6 12– has been found to catalyze a number of important chemical reactions, but the close proximity of a poorly organized external solvent environment can serve as an impediment to the reaction chemistry. Here, using ab initio molecular dynamic calculations and electric field analyses, we find that metal substitution of indium for gallium lowers the total and electrostatic activation barriers by ~3–4 kcal/mol in water solvent. Using an energy decomposition analysis of the interaction of water solvent with the metal vertices, we find that Pauli repulsion between the metal and water decreases, allowing the water coordination with the metal to increase, upon substitution of In for Ga. We therefore determine that the stabilization of the transition state is due to a better arrangement of water molecules around the metal vertices of In 4 L 6 12– and provides a proof of concept of how to redesign the nanocage scaffold in order to reduce the solvent reorganization energy.

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

Linear Combination of Atomic Dipoles to Calculate the Bond and Molecular Dipole Moments of Molecules and Molecular Liquids

Here we report a linear combination of atomic dipole (LCAD) method for calculating the bond dipole moments of molecules. We show that the LCAD method reproduces the known molecular dipole moments of small to large molecules with a small error with respect to experimental and benchmark ab initio calculations, and molecular dipole distributions of bulk water that agree with maximally localized Wannier functions. The bond dipole moments derived from LCAD are also chemically interpretable in terms of the trend in bond ionicity in going from neutral to charged molecules. Moreover, the LCAD method accurately captures the influence of electric fields, supported by the correct trend in the change of the dipole moment under a uniform external electric field. The better grounding of bond dipole calculations indicates that it should also serve as a useful approach to bond dipole-field models used in catalysis or to reconstruct the small dipole of a H-terminated graphene flake.

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

Challenges for density functional theory: calculation of CO adsorption on electrocatalytically relevant metals

Density Functional Theory (DFT) is currently the most tractable choice of theoretical model used to understand the mechanistic pathways for electrocatalytic processes such as CO 2 or CO reduction. In this work, we assess the performance of two DFT functionals designed specifically to describe surface interactions, RTPSS and RPBE, as well as two popular meta-GGA functionals, SCAN and B97M-rV, that have not been a priori optimized for better interfacial properties. We assess all four functionals against available experimental data for prediction of bulk and bare surface properties on four electrocatalytically relevant metals, Au, Ag, Cu, and Pt, and for binding CO to surfaces of these metals. To partially mitigate issues such as thermal and anharmonic corrections associated with comparing computations with experiments, molecular benchmarks against high level quantum chemistry are reported for CO complexes with Au, Ag, Cu and Pt atoms, as well as the CO–water complex and the water dimer. Overall, we find that the surface modified RPBE functional performs reliably for many of the benchmarks examined here, and the meta-GGA functionals also show promising results. Specifically B97M-rV predicts the correct site preference for CO binding on Ag and Au (the only functional tested here to do so), while RTPSS performs well for surface relaxations and binding of CO on Pt and Cu.

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

An isolated water droplet in the aqueous solution of a supramolecular tetrahedral cage

Water under nanoconfinement at ambient conditions has exhibited low-dimensional ice formation and liquid–solid phase transitions, but with structural and dynamical signatures that map onto known regions of water’s phase diagram. Using terahertz (THz) absorption spectroscopy and ab initio molecular dynamics, we have investigated the ambient water confined in a supramolecular tetrahedral assembly, and determined that a dynamically distinct network of 9 ± 1 water molecules is present within the nanocavity of the host. The low-frequency absorption spectrum and theoretical analysis of the water in the Ga 4 L 6 12– host demonstrate that the structure and dynamics of the encapsulated droplet is distinct from any known phase of water. A further inference is that the release of the highly unusual encapsulated water droplet creates a strong thermodynamic driver for the high-affinity binding of guests in aqueous solution for the Ga 4 L 6 12– supramolecular construct.

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

Spherical trihedral metallo-borospherenes

The discovery of borospherenes unveiled the capacity of boron to form fullerene-like cage structures. While fullerenes are known to entrap metal atoms to form endohedral metallofullerenes, few metal atoms have been observed to be part of the fullerene cages. Here we report the observation of a class of remarkable metallo-borospherenes, where metal atoms are integral parts of the cage surface. We have produced La 3 B 18 - and Tb 3 B 18 - and probed their structures and bonding using photoelectron spectroscopy and theoretical calculations. Global minimum searches revealed that the most stable structures of Ln 3 B 18 - are hollow cages with D 3h symmetry. The B 18 -framework in the Ln 3 B 18 - cages can be viewed as consisting of two triangular B6 motifs connected by three B 2 units, forming three shared B 10 rings which are coordinated to the three Ln atoms on the cage surface. These metallo-borospherenes represent a new class of unusual geometry that has not been observed in chemistry heretofore.

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

Interplay of water and a supramolecular capsule for catalysis of reductive elimination reaction from gold

Abstract Supramolecular assemblies have gained tremendous attention due to their ability to catalyze reactions with the efficiencies of natural enzymes. Using ab initio molecular dynamics, we identify the origin of the catalysis by the supramolecular capsule Ga 4 L 6 12− on the reductive elimination reaction from gold complexes and assess their similarity to natural enzymes. By comparing the free energies of the reactants and transition states for the catalyzed and uncatalyzed reactions, we determine that an encapsulated water molecule generates electric fields that contributes the most to the reduction in the activation free energy. Although this is unlike the biomimetic scenario of catalysis through direct host-guest interactions, the electric fields from the nanocage also supports the transition state to complete the reductive elimination reaction with greater catalytic efficiency. However it is also shown that the nanocage poorly organizes the interfacial water, which in turn creates electric fields that misalign with the breaking bonds of the substrate, thus identifying new opportunities for catalytic design improvements in nanocage assemblies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗