Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “QTAIM”

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.

12 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Synergistic Density Functional Theory and Molecular Dynamics Approach to Elucidate PNIPAM–Water Interaction Mechanisms

This study employs Density Functional Theory (DFT) and Molecular Dynamics (MD) simulations to investigate interactions between water molecules and Poly(N-isopropylacrylamide) (PNIPAM). DFT reveals preferential water binding sites, with enhanced binding energy observed in the linker zone. Quantum Theory of Atoms in Molecules (QTAIM) and electron localization function (ELF) analyses highlight the roles of hydrogen bonding and steric hindrance. MD simulations unveil temperature-dependent hydration dynamics, with structural transitions marked by changes in the radius of gyration (Rg) and the radial distribution function (RDF), aligning with DFT findings. Our work goes beyond prior studies by combining a DFT, QTAIM and MD simulations approach across different PNIPAM monomer-to-30mer structures. It introduces a systematic quantification of pseudo-saturation thresholds and explores water clustering dynamics with structural specificity, which have not been previously reported in the literature. These novel insights establish a more complete molecular-level picture of PNIPAM hydration behavior and temperature responsiveness, emphasizing the importance of amide hydrogen and carbonyl oxygen sites in hydrogen bonding, which weakens above the lower critical solution temperature (LCST), resulting in increased hydrophobicity and paving the way for understanding water sorption mechanisms, offering guidance for future applications such as dehumidification and atmospheric water harvesting.

08 HYDROGEN↗
DOE OSTI2026

Cosolvent-tuned interactions in ionic liquids: A vibrational and quantum-chemical study of ethylene glycol ratio effects

Ionic liquids (ILs) are attractive media for CO 2 capture but remain limited by viscosity and cost. Blending ILs with ethylene glycol (EG) is a practical route to mitigate these constraints, yet the molecular origins of cosolvent effects and their dependence on composition are not well resolved. We combine Fourier-transform infrared (FT-IR) spectroscopy with quantum-chemical (DFT) analysis to elucidate how the IL:EG molar ratio modulates intermolecular interactions and electronic structure. Computed vibrational frequencies enable mode assignment and deconvolution of overlapping bands, revealing systematic, ratio-dependent shifts and broadenings in (i) EG O–H stretching, (ii) cation and EG C–H stretchings (imidazolium C2–H, C4–H, C5–H, methyl and ethyl groups, -CH2 of EG), (iii) anion signature modes (e.g., CN motifs), and (iv) EG C–O and C–C stretchings, consistent with the redistribution of hydrogen-bonding networks. Molecular electrostatic potential (MESP) maps quantify attenuation of extreme potential regions with increasing EG, indicating progressive screening of cation–anion electrostatic interactions. Quantum Theory of Atoms in Molecules (QTAIM) identifies emergent bond critical points between EG and the IL ions, while Reduced Density Gradient–Noncovalent Interaction (RDG–NCI) analysis differentiates strong directional hydrogen bonds from dispersive contacts across compositions. Together, these results show that EG fraction controls a switch from predominantly ion–ion to mixed ion–EG coordination, altering local polarity and polarizability that underlie the observed FT-IR trends. The framework provides composition–structure–spectrum relationships that can guide rational selection of IL:EG ratios to balance favorable molecular interactions with practical performance targets in scalable CO 2 capture systems.

DAC↗
DOE OSTI2026

Pyridine Complexes of Iodobismuthate(III) Anions

We report a rare family of pyridine-coordinated iodobismuthate(III) salts supported by alkyltriphenylphosphonium and tetraphenylphosphonium cations. Reactions of BiI 3 with Ph 3 PR + I − (R = Me, Et, n Pr, n Bu, Ph) in neat pyridine, followed by crystallization, yield structurally tunable bismuth-halide-pyridine anions dictated by reagent stoichiometry. Combination of BiI 3 and Ph 3 PR + I − in 2:1 ratio produced [Ph 3 PR] 2 [BiI 5 Py], 1 (R = Me, Et, n Pr, Ph), while combination in 1:1 ratio resulted in three compounds: [Ph 3 PR][cis-BiI 4 Py 2 ], 2 (R = n Pr, Ph), [Ph 3 PR][trans-BiI 4 Py 2 ], 3 (R = Me, Et, Ph), and [Ph 3 PR] 2 [transoid-Bi 2 I 8 Py 2 ], 4 (R = Me, Et, n Pr, n Bu, Ph). In many cases, the compounds were isolated as Py or Et 2 O solvates, and in some cases, multiple degrees of solvation or polymorphism were encountered. Hirshfeld analysis of 1–4 showed the major anion–cation/anion/solvent interactions to be H⋯I, H⋯H, and C⋯H. Diffuse reflectance measurements of representative compounds, all of which were yellow-orange to red-orange, revealed bandgaps in the range of 1.9–2.2 eV, where density-of-states KS-DFT calculations attribute the absorption to metal-centered charge transfer within the anionic unit. NLMO and QTAIM analyses further indicate predominantly ionic Bi(III)–I/pyridine bonding with robust inner-sphere coordination that is insensitive to anion speciation.

Bismuth Compounds↗
DOE OSTI2024

Ternary Complexes of BiI 3 /CuI and SbI 3 /CuI with Tetrahydrothiophene

Reactions of BiI 3 /CuI mixtures with tetrahydrothiophene (THT) in toluene produce 2-D sheet networks BiCu 3 I 6 (THT) n (n = 2, 3, or 4), depending on reaction conditions. All three structures are based on BiI 6 octahedra, which share pairs of (μ 2 -I) 2 with Cu 3 (THT) n units. BiCu 3 I 6 (THT) 2 features Cu 2 (μ 2 -I) 2 rhombs with close Cu···Cu interactions and is accompanied by formation of the very complex HBi 3 Cu 12 I 22 (THT) 8 . Reactions of SbI 3 /CuI with THT in toluene produced a SbCu 3 I 6 (THT) 2 network shows Cu 3 (μ 2 -THT) 2 units, like its Bi congener, but Cu 6 (μ 2 -I) 6 barrels rather than rhombs. Isolated SbI 3 units are stacked above the Cu 6 I 6 barrels. A molecular compound, Sb 3 Cu 3 I 12 (THT) 6 consists of a face-sharing Sb 3 I 12 stack, in which the Cu-THT units are bonded in asymmetric fashion about the central SbI 6 . Metal-halide bonds were investigated via QTAIM and NLMO analyses, demonstrating that these bonds are largely ionic and occur between the Bi/Sb and I p orbitals. Hirshfeld analysis shows significant H···H and H···I interactions. Diffuse reflectance spectroscopy (DRS) reveals band edges for the Bi species of 1.71–1.82 eV, while those for the neutral Sb complexes are in the range of 1.94–2.06 eV. Mapping of the electronic structure via density of state calculations indicates population of antibonding Bi/Sb–I orbitals in the excited state.

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

Elucidation of an Unusually Long Pu–N Bond in a Plutonium(III)-Tetrazolate Complex

Four trivalent, f-element tetrazolate hydrate complexes [M(H 2 O) 9 ](Hdtb) 3 ·nH 2 O (Nd1, n = 7 and Pu1, n = 9; dtb 2– = 1,3-di(tetrazolate-5-yl)benzene) and [M(Hdtb)(H 2 O) 8 ](dtb)·11H 2 O (Nd2 and Pu2) were prepared using metathesis reactions. These complexes contain hydrated M(III) cations, but in the latter complexes, Nd2 and Pu2, one of the water molecules has been displaced by a long interaction between the M(III) cation and a Hdtb – anion. Notably, the Pu(III)–N bond in Pu2, representing the longest IX Pu(III)–N ( IX = nine coordinate) bond reported has a length of 2.8338(15) Å and is slightly shorter than the Nd(III)–N bond length of 2.8425(13) Å in Nd2. Analysis of bond lengths, Wiberg bond indices (WBI), natural localized molecular orbitals (NLMOs), and quantum theory of atoms in molecules (QTAIM) reveals that the metal contribution to the Pu(III)–N bond is marginally greater than that of the Pu(III)–OH 2 bonds in Pu2 and the Nd(III)–N bond in Nd2. Thus, this rather long M–N interaction provides an example where the expectation that An(III) compounds exhibit greater covalency with soft donor ligands compared to harder ligands fails. Furthermore, the absorption spectra of Pu1 and Pu2 further support this observation, highlighting a surprising degree of similarity in their electronic structures.

Covalent bonding↗
DOE OSTI2024

Role of Intermolecular Interactions in Deep Eutectic Solvents for CO 2 Capture: Vibrational Spectroscopy and Quantum Chemical Studies

Recent research and reviews on CO 2 capture methods, along with advancements in industry, have highlighted high costs and energy-intensive nature as the primary limitations of conventional direct air capture and storage (DACS) methods. In response to these challenges, deep eutectic solvents (DESs) have emerged as promising absorbents due to their scalability, selectivity, and lower environmental impact compared to other absorbents. However, the molecular origins of their enhanced thermal stability and selectivity for DAC applications have not been explored before. Therefore, the current study focuses on a comprehensive investigation into the molecular interactions within an alkaline DES composed of potassium hydroxide (KOH) and ethylene glycol (EG). Combining Fourier transform infrared (FT-IR) and quantum chemical calculations, the study reports structural changes and intermolecular interactions induced in EG upon addition of KOH and its implications on CO 2 capture. Experimental and computational spectroscopic studies confirm the presence of noncovalent interactions (hydrogen bonds) within both EG and the KOH-EG system and point to the aggregation of ions at higher KOH concentrations. Additionally, molecular electrostatic potential (MESP) surface analysis, natural bond orbital (NBO) analysis, quantum theory of atoms-in-molecules (QTAIM) analysis, and reduced density gradient-noncovalent interaction (RDG-NCI) plot analysis elucidate changes in polarizability, charge distribution, hydrogen bond types, noncovalent interactions, and interaction strengths, respectively. Evaluation of explicit and hybrid models assesses their effectiveness in representing intermolecular interactions. This research enhances our understanding of molecular interactions in the KOH-EG system, which are essential for both the absorption and desorption of CO 2 . The study also aids in predicting and selecting DES components, optimizing their ratios with salts, and fine-tuning the properties of similar solvents and salts for enhanced CO 2 capture efficiency.

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

Reactions of Np(VI) and Pu(VI) with Phenanthroline Result in Bending the [NpVIO2]2+ Unit and a Reduced Pu(V) Species

Despite decades of actinyl chemistry, genuinely bent actinyl structures remain rare beyond uranium. We report the first crystallographic characterization of a bent neptunyl(VI) complex, NpO 2 Cl 2 (phen) 2, along with a linear plutonyl species that, unexpectedly, adopts the +V oxidation state. Coordination of two 1,10-phenanthroline ligands to NpO 2 2+ enforces pronounced bending of the O yl –Np–O yl unit to 162° through steric clashes with the axial phenathroline ligand, representing the sharpest angle reported for any neptunyl(VI) complex. In contrast, synthesis with PuO 2 2+ produces a reduced linear Pu(V) species, PuO 2 Cl(phen) 2, underscoring the redox lability and the reluctance of plutonium to distort. Raman and IR spectra yield the first experimental stretching and interaction force constants for plutonyl(V). Comparison with values for plutonyl(VI) and bent and linear uranyl(VI) and neptunyl(VI) compounds shows that reduction perturbs the actinyl bond more than bending within one oxidation state. Electrochemical and spectroscopic studies clarify phenanthroline binding to Pu(VI) and the system’s reactivity. Quantum chemical calculations indicate that bending is energetically favored in the U–Np–Pu series, but the stabilization energy decreases as the actinide atomic number increases. NBO and QTAIM analyses reveal systematic trends: increasing 5f occupancy and decreasing An–O yl bond covalency from U to Pu. These results demonstrate that ligand-induced bending in neptunyl(VI) species is not only electronically feasible but also synthetically achievable, though plutonyl(VI) is constrained by redox reactivity. This work expands the frontier of nonlinear actinyl chemistry and illuminates how structure, oxidation state, and electronic configuration interrelate across the actinyl series.

actinides↗
DOE OSTI2026

Covalency of M–N Bonds in Isomorphous Lanthanide and Actinide 5-(2-Pyridyl)-1H-tetrazolate Complexes

Experimental and computational analyses of [M(pdtz) 3 (H 2 O) 3 ]·3.5H 2 O (M 3+ = Pu 3+ −Cm 3+ , La 3+ −Nd 3+ , and Sm 3+ −Ho 3+ , pdtz− = 5-(2-pyridyl)-1H-tetrazolate) were conducted to understand potential differences in bonding between lanthanide and actinide complexes with a N-donor ligand. Structural analyses show that the An−N bond distances in the Pu 3+ , Am 3+ , and Cm 3+ complexes are within error of one another. Whereas in the lanthanide series, there is a nearly linear decrease in the Ln−N bond lengths from La 3+ to Ho 3+ (excluding Pm 3+ ). The An−N bond lengths are ∼0.015 Å shorter than their similarly-sized lanthanide analogs, in agreement with computational results that suggest greater covalent character in these bonds versus those with lanthanides. QTAIM analysis indicates that the An−N orbital mixing remains essentially unchanged from Pu 3+ to Cm 3+ , consistent with the nearly identical An−N bond lengths. However, upon deconvolution of the NLMOs into orbital compositions, the metal orbital contributions to An−N bonding decreases slightly overall wherein the 6d involvement remains constant, 7s involvement slightly increases, and 5f participation decreases. The molecular orbital energy diagram indicates that energy degeneracy between the 5f metal and 2p ligand orbitals increases from Pu 3+ to Cm 3+ and counteracts the contraction of the 5f orbtials. Together with prior reports of decreasing energy degeneracy between 5f and 3p orbitals from Np 3+ to Cf 3+ , these observations provide guidance on understanding how chemical bonding evolves in the actinide series.

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

Non–linear bonding trends in maleonitrile-1,2–dithiolate complexes of the transuranium actinides

The trivalent actinides are produced in the nuclear fuel cycle during power production and provide the largest long-term radiation dose in used nuclear fuel. It is ideal for these elements to be removed from used nuclear fuel for disposal and a necessity for fuel recycling. A key challenge to this is the similarity of chemical behavior of the trivalent actinides to the lanthanides that are also present as fission products in used fuel. Thus far, some of the most effective separations of actinides from lanthanides utilise chelating agents containing sulfur moieties such as dithiophosphinates that selectively bind to actinide ions because of a greater bond covalency relative to lanthanide ions. Typically, greater differences between actinide and lanthanide ions are observable the more ligands and chelators bonds have a covalent character. Here, a series of complexes of the trivalent actinides Np(III) through Cf(III) (excluding Bk(III)) with maleonitrile-1,2-dithiolate (mnt 2– ) are synthesized along with their lanthanide counterparts (La(III) – Nd(III), Sm(III) – Gd(III), Dy(III)), in order to characterize the nature of chemical bonds with these metal ions and a polarizable, non-innocent, sulfur-donor ligand. The metal-sulfur bonds in these complexes trend shorter than measured for lanthanides with equivalent ionic radii. However, particularly large deviations are observed in the neptunium and plutonium complexes in both structure and bonding, resulting in a nonlinear bond length trendline for the actinide series. Density Functional Theory (DFT) calculations with Quantum Theory of Atoms in Molecules (QTAIM) and Natural Bond Order (NBO) analyses indicate that for the neptunium and plutonium complexes, the presence of increased 5f-orbital participation, energy degeneracy of the metal and ligand orbitals, and the structure packing result in shortened M–S bonds. The stabilization of the energy of the 5f-orbitals and the decrease in f-contribution to bonding orbitals in the later actinides results in structural properties more similar to the lanthanide complexes.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2024

Volatile fatty acid extraction from fermentation broth using a hydrophobic ionic liquid and in situ enzymatic esterification

Efficient recovery of volatile fatty acids (VFAs) from fermentation broth is a challenge due to low VFA titers and thus limits the commercialization of VFA production using biological routes. Liquid–liquid extraction using hydrophobic ionic liquids (ILs) shows great promise for the extraction and esterification of hydrophilic VFAs. In this study, several ILs were evaluated to select a water-immiscible and efficient extraction solvent. The selected IL, trihexyltetradecyl phosphonium dibutylphosphate ([P 666,14 ][DBP]), gave a cumulative VFA extraction of around 842.8 mg per g IL. The predicted excess enthalpy (H E ) and logarithmic activity coefficients ln(γ) using the COSMO-RS model were validated with the experimentally obtained VFA recovery from fermentation broth. To understand the extraction mechanism of VFAs, quantum theory of atoms in molecules (QTAIM) and noncovalent interaction (NCI) were performed. The results suggest that long chain fatty acids exhibit strong van der Waals interaction with the DBP anion leading to higher VFA extraction. The enzymatic esterification of VFAs with ethanol in [P 666,14 ][DBP] was optimized using the Box–Behnken response surface design of experiment. Under the optimized conditions, up to 83.7% of hexanoic acid was converted to ethyl esters, while other shorter chain VFAs have lower conversion efficiency (38.3–63.2%).

09 BIOMASS FUELS↗
DOE OSTI2025

Unveiling the covalency of versatile Pu(iii)-N bonds in a unique plutonium(iii) complex

A trivalent plutonium–pyrazinyl–tetrazolate complex Na 2 [Pu(Hdtp)(dtp) 2 (H 2 O) 4 ]·9H 2 O (Pu_dtp, H 2 dtp = 2,3-di-1H-tetrazol-5-ylpyrazine) was synthesized through metathesis reaction of plutonium bromide and Na 2 (dtp)·2H 2 O in water. This structure is particularly notable among complexes formed by trivalent f-elements and the dtp 2− ligand in aqueous media. In contrast to other trivalent f-elements, including all Ln 3+ (with the exception of Pm 3+ ) and Cm 3+ , preferentially coordinated with eight water molecules rather than the nitrogen donors of the dtp 2− ligand, Pu 3+ exhibits a distinct affinity for nitrogen coordination. This observation provides strong evidence that the 5f electrons in Pu 3+ are more delocalized than other studied trivalent f-elements. In Pu_dtp, three distinct Pu(III)–N bonds are present: one Pu(III)–N5 from pyrazinyl, one Pu(III)–N4 from the least electronegative nitrogen in the tetrazolate, and three Pu(III)–N1/N2/N3 from the most electronegative nitrogens in the tetrazolate. Experimental Pu(III)–N bond lengths, Wiberg bond indices (WBI), natural localized molecular orbitals (NLMO), quantum theory of atoms in molecules (QTAIM), and energy decomposition analysis (EDA), reveal a covalency trend: Pu(III)–N from the most electronegative nitrogen in tetrazolate > Pu(III)–N from the least electronegative nitrogen in tetrazolate > Pu(III)–N from pyrazinyl. This trend arises from the increased negative charge on the most electronegative nitrogen atoms in the tetrazolate ring, enhancing electrostatic Pu–N1/N2/N3 interactions. These stronger electrostatic interactions lead to shorter bond lengths, thereby enhancing orbital overlap and greater covalency, compared to the less electronegative nitrogen in tetrazolate (Pu–N4) and the neutral pyrazinyl nitrogen (Pu–N5).

Bai, Zhuanling [Colorado School of Mines, Golden, ↗
DOE OSTI2026

An–imidophosphorane (An = U–Pu) bond covalency and proton-coupled electron transfer thermodynamics driven by orbital energy matching

A series of mid-actinide (An = U–Pu) tetrahomoleptic complexes supported by highly electron-donating imidophosphorane ligands, NPC ([NP t Bu(pyrr) 2 ] − , where t Bu = C(CH 3 ) 3 ; pyrr = pyrrolidinyl = N(C 4 H 8 )), are systematically investigated computationally and experimentally to elucidate the nature of actinide–ligand (An–L) covalency across the An 3+/4+/5+ oxidation states. Trends in An–L bonding and redox properties for these complexes, together with their protonated counterparts, are examined using orbital-, electron density-, and energy-decomposition-based methods. This integrated approach reveals progressively improved energy matching between α-spin An 5f and N im 2p orbitals with increasing atomic number and oxidation state, becoming particularly pronounced in the ligand-dominant π-bonding orbitals of An 4+ and An 5+ . In contrast to the An 3+ species, the enhanced An 5f π contributions in the higher-valent counterparts drive the increase in An–N im covalency for later An, thereby inverting the covalency trend to U < Np < Pu. Redistribution of electron density towards the An and N im atomic basins due to the growing energy-matching assisted covalency correlates with higher pKa values and increased N im –H bond dissociation free energies in protonated An 4+ complexes. Electron density at Nim in An 4+ shows a linear correlation with the p K a values calculated via the Bordwell equation. Calculations predict a cathodic shift of 0.84–1.00 V in the redox couples upon protonation, a trend validated when experimentally accessible. These findings demonstrate an increasing role of covalency driven by orbital energy matching from U to Pu in tuning the thermodynamic driving force for proton-coupled electron transfer in the An 5+ species.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗