Search NASA⌕ Search

SEARCH · Search NASA

Results for “HAND”

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 199 records · Page 11

Scanning Electrochemical Microscopy for Kinetic Investigations in Viscous Deep Eutectic Solvents: Identifying Practical Approach Curves and Deviations from Electron Transfer Models

Determining heterogeneous electrochemical electron transfer (ET) kinetics in electrolytes with a wide range of physical properties is of great interest for achieving high-performance redox flow batteries. Among such electrolytes, concentrated hydrogen-bonded electrolytes (CoHBEs), including deep eutectic solvents (DESs), have recently garnered significant attention. Unfortunately, traditional Tafel analysis using macroelectrodes often encounters issues with mass transfer limitations in CoHBEs with high viscosities, thereby restricting kinetic analysis to a narrow potential window. Here, in this work, we introduce a methodology for evaluating ET kinetics in viscous DES using the scanning electrochemical microscopy (SECM). We first determined practical solutions to SECM tip positioning in ethaline DES, which yield pseudopositive feedback responses. Lattice Boltzmann method (LBM) simulations helped us rationalize the impact of the fluid and concentration fields, as well as tip geometry, tip approach velocity v, and the solvent viscosity ηs, on the shape of the approach curves. In addition to successfully recreating approach curves over a variety of conditions, we found that approaching a conductor ensured a practical point where the normalized tip response (Ni T = 2) converged at L = 0.7 within ∼10% error regardless of tip velocity. With positioning capabilities at hand, we investigated the kinetics of Fe 3+ /Fe 2+ redox couple in aqueous and the ethaline media. The experimental kinetic results were interpreted using the Butler–Volmer (BV) and Marcus–Hush–Chidsey (MHC) models. For ethaline, a nonideal kinetic behavior was observed, potentially attributed to solvent dynamics within DESs or to the interplay of chloride anions in the charge transfer process.

electrodes↗

Revealing Structural Evolution of Nickel Phosphide-Iron Oxide Core–Shell Nanocatalysts in Alkaline Medium for the Oxygen Evolution Reaction

Metal phosphide-containing materials have emerged as a potential candidate of non-precious metal-based catalysts for alkaline oxygen evolution reaction (OER). While it is known that metal phosphide undergoes structural evolution, considerable debate persists regarding the effects of dynamics on the surface activation and morphological stability of the catalysts. In this study, we synthesize NiP x -FeO x core-shell nanocatalysts with an amorphous NiP x core designed for enhanced OER activity. Using ex-situ X-ray absorption spectroscopy, we elucidate the local structural changes as a function of cyclic voltammetry cycles. Our studies suggest that the presence of corner-sharing octahedra in the FeO x shell improves structural rigidity through interlayer cross-linking, thereby inhibiting the diffusion of OH - /H 2 O. Thus, the FeO x shell preserves the amorphous NiP x core from rapid oxidation to Ni 3 (PO 4 ) 2 and Ni(OH) 2 . On the other hand, the incorporation of Ni from the core into the FeO x shell facilitates absorption of hydroxide ions for OER. As a result, the Ni/Fe(OH) x at the surface oxidizes to the active γ-(oxy)hydroxide phase under the applied potentials, promoting OER. This intriguing synergistic behavior holds significance as such synthetic route involving the FeO x shell can be extended to other systems, enabling manipulation of surface adsorption and diffusion of hydroxide ions. These findings also demonstrate that nanomaterials with core-shell morphology can be tuned to leverage the strength of each metallic component for improved electrochemical activities.

25 ENERGY STORAGE↗

Defect Diffusion Graph Neural Networks for Materials Discovery in High-Temperature Energy Applications

Here, the migration of crystallographic defects dictates material properties and performance for a plethora of technological applications. Density functional theory (DFT)-based nudged elastic band (NEB) calculations are a powerful computational technique for predicting defect migration activation energy barriers, yet they become prohibitively expensive for high-throughput screening of defect diffusivities. Without introducing hand-crafted (i.e., chemistry- or structure-specific) descriptors, we propose a generalized deep learning approach to train surrogate models for NEB energies of vacancy migration by hybridizing graph neural networks with transformer encoders and simply using pristine host structures as input. With sufficient training data, computationally efficient and simultaneous inference of vacancy defect thermodynamics and migration activation energies can be obtained to compute temperature-dependent vacancy diffusivities and to down-select candidates for more thorough DFT analysis or experiments. Thus, as we specifically demonstrate for potential water-splitting materials, candidates with desired defect thermodynamics, kinetics, and host stability properties can be more rapidly targeted from open-source databases of experimentally validated or hypothetical materials.

14 SOLAR ENERGY↗

Ion Size Effects on the Thermodynamic, Kinetic, and Mechanical Properties during Ion Exchange in Solid-State Electrolytes

Ion exchange offers a pathway to impose residual compressive stresses to mitigate the electro-chemo-mechanical cracking of solid-state electrolytes such as lithium lanthanum zirconium oxide. This study uses a coupled multiscale framework (integrating density functional theory (DFT), molecular dynamics (MD), and continuum modeling) to examine how exchange ion size influences stress, diffusion, fracture toughness, and electronic properties. Larger isovalent ions (Na + , Ag + , K + ) were exchanged with Li + , with DFT confirming their preference for octahedral 96h sites and a linear relationship between ion size and chemical free expansion coefficient. MD simulations reveal stress and concentration effects on exchange ion diffusivity at elevated temperatures, with Na + and Ag + maintaining favorable mobility while K + showing concentration-dependent clustering. Continuum modeling predicts the range of fracture strength improvements and the required ion exchange concentration profile. It was shown that a 5% surface exchange concentration can induce ∼0.6 GPa of surface compressive stress using Na + and ∼1.0 GPa of surface compressive stress using Ag + . On the other hand, larger ion exchange species may penalize Li + diffusivity by increasing the activation volume and activation energy. Interestingly, Na + has a negligible penalty on Li-ion diffusivity. The room temperature Li + ion diffusivity is reduced by ∼40% with Ag + ion exchange. Electronic band structure analysis shows no size-dependent change in the bandgap, though Ag + introduces localized defect states near the valence band maximum. This study highlights ion size as a key factor in optimizing LLZO properties, offering a framework to improve the solid-state battery performance.

Jagad, Harsh D. [Brown Univ., Providence, RI (Unit↗

Ultrafast Sintering and Dopant Effects in Garnet LLZO Solid Electrolytes

High-throughput, low-cost manufacturing, and optimization of solid electrolytes are necessary for the adoption of solid-state batteries. In this work, garnet-type Li 7 La 3 Zr 2 O 12 (LLZO) with different aliovalent dopants, 𝑇𝑎$_{^·_{𝑍𝑟}}$, 𝐴𝑙$^{··}_{𝐿𝑖}$, and 𝐺𝑎$^{··}_{𝐿𝑖}$, have been ultrafast-sintered with different temperature ramping rates. The densification behavior, phases, their evolution, and surface chemistry of different LLZO have been investigated and linked to their electrochemical performances. It has been shown that LLZO with 𝑇𝑎$_{^·_{𝑍𝑟}}$ dopant demonstrates the highest garnet phase purity and overall best electrochemical performances, and ultrafast sintering further improves densification, ionic conductivity, and electrochemical stability. On the other hand, LLZO doped with 𝐴𝑙$^{··}_{𝐿𝑖}$ and 𝐺𝑎$^{··}_{𝐿𝑖}$ are reaching higher cubic phase purities and ionic conductivities via conventional sintering, indicating undesirable dopant migration and segregation during the ultrafast sintering process. In conclusion, these findings provide insights into the manufacturing of solid electrolyte materials.

36 MATERIALS SCIENCE↗

Variable Surface Termination and Ligand Passivation of Lead Sulfide Nanocrystals Synthesized with Excess Lead Chloride

The surface termination and ligand passivation of semiconducting nanocrystals (NCs) impact the stability, optical properties, and self-assembly of NCs. In this work, we definitively characterize the surface of lead sulfide (PbS) NCs synthesized from excess PbCl 2 . With a combination of small-angle neutron scattering (SANS), photoluminescence, and 1D and 2D NMR experiments, we show that the surface termination of PbS NCs depends on the presence of PbCl 2 during ligand exchange. When excess PbCl 2 is removed prior to ligand exchange, PbS[RNH 3 + Cl – ] NCs are obtained, which are terminated by a monolayer of PbCl x on the {100} PbS facet and passivated by oleylammonium chloride ligands. On the other hand, when excess PbCl 2 remains in solution during ligand exchange, lead oleate forms and attaches to the {111} PbS facets of PbS[RNH 3 + Cl – ] NCs, creating PbS@PbCl x NCs. PbS@PbCl x NCs are coated in an epitaxial layer of PbCl x on both the {100} PbS and {111} PbS facets, making them 0.3–0.4 nm larger on average than PbS[RNH 3 + Cl – ] NCs with identical absorption wavelengths. Additionally, PbS@PbCl x NCs have a consistently higher photoluminescence quantum yield and longer photoluminescence lifetimes than PbS[RNH 3 + Cl – ] NCs. This study clarifies the surface structure of PbS NCs synthesized from excess PbCl 2 , highlighting ligand exchange strategies and reconciling observations from across the literature.

Anions↗

Interfacial Hydrogen-Bond Dynamics in Transition Metal Compounds

Understanding how water behaves when confined within atomic layers of active transition-metal carbides, nitrides, and carbonitrides is essential for uncovering the fundamental principles needed to engineer solid–liquid interfaces at the atomic scale. Yet, how lattice element chemistry and surface termination groups collectively regulate the structure and mobility of such interlayer water remains poorly understood. Here, we present a composition-controlled investigation of interlayer water dynamics in layered transition-metal nitride, carbide, and carbonitride systems using a systematic integration of quasi-elastic neutron scattering (QENS), ab initio molecular dynamics (AIMD) simulations, and density functional theory (DFT) calculations. QENS measurements show that nitride-rich systems host mobile, translationally diffusing water with thermally activated self-diffusion coefficients on the order of 10 –10 m 2 s –1 , whereas mixed C/N lattices confine water to localized, nontranslational motion that is insensitive to temperature. AIMD and DFT reveal that lattice C/N chemistry and surface functional group composition reshape the first hydration layer by modulating the surface electronic structure and termination-dependent hydrogen-bond networks, leading to pronounced differences in water ordering and thermal resilience. On the other hand, fully carbide systems exhibit intermediate behavior, highlighting that water mobility is not primarily controlled by the hydration level alone but by the coupling between lattice composition and surface chemistry. Overall, this study establishes how surface chemistry and lattice composition jointly control interfacial hydrogen bond dynamics, offering a mechanistic framework for designing transition-metal layered materials with tailored interfacial transport properties.

Hydration↗

Multireference Methods for Chemistry and Materials Science: Automated Active Spaces, Efficient Dynamic Correlation, and Extended Systems

While multiconfigurational approaches have long been relegated to expert practitioners working on a case-by-case basis, recent developments have increasingly made these methods more routine and applicable to broader sets of systems. This article outlines the state-of-the-art in multiconfigurational approaches, with an emphasis on moving from delicate hand-selected pathways through configuration space toward more robust and efficient approaches to treating a host of challenging chemical systems accurately. First, we overview recent work in automated active-space selection, which has enabled increasingly large-scale applications of multireference methods to modeling vertical excitations and reactivity. Second, we highlight the increasingly efficient methods for recovering correlation energy beyond the active space, as headlined by extensions of pair-density functional theory and its role in accurate and efficient treatment of excited-state dynamics and its utilization to train machine-learned potentials. Finally, we highlight recent efforts to treat extended systems that until recently have lied beyond the traditional limits of active-space methods, giving center stage to product-form wave functions of the localized active space family of methods that allow for the computation of multiconfigurational band structures. These recent advancements point to a broader use of multireference approaches for high-impact chemical and materials science applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Prioritization of Early-Stage Research and Development of a Hydrogel-Encapsulated Anaerobic Technology for Distributed Treatment of High Strength Organic Wastewater

This study aims to support the prioritization of research and development (R&D) pathways of an anaerobic technology leveraging hydrogel-encapsulated biomass to treat high-strength organic industrial wastewaters, enabling decentralized energy recovery and treatment to reduce organic loading on centralized treatment facilities. To characterize the sustainability implications of early-stage design decisions and to delineate R&D targets, an encapsulated anaerobic process model was developed and coupled with design algorithms for integrated process simulation, techno-economic analysis, and life cycle assessment under uncertainty. Across the design space, a single-stage configuration with passive biogas collection was found to have the greatest potential for financial viability and the lowest life cycle carbon emission. Through robust uncertainty and sensitivity analyses, we found technology performance was driven by a handful of design and technological factors despite uncertainty surrounding many others. Hydraulic retention time and encapsulant volume were identified as the most impactful design decisions for the levelized cost and carbon intensity of chemical oxygen demand (COD) removal. Encapsulant longevity, a technological parameter, was the dominant driver of system sustainability and thus a clear R&D priority. Ultimately, we found encapsulated anaerobic systems with optimized fluidized bed design have significant potential to provide affordable, carbon-negative, and distributed COD removal from high strength organic wastewaters if encapsulant longevity can be maintained at 5 years or above.

Anaerobic Treatment↗

Design of Experiments for Dynamic Test Runs in Solvent-Based CO 2 Capture Pilot Plants

Test runs in the pilot plants consume significant resources, and therefore, the learning from test runs should be maximized. Test runs conducted in the pilot plants are often steady state. It takes several hours for reaching steady-state in the pilot plants, and thus, the duration of the test runs needs to be long even for collecting few steady-state data points. On the other hand, a large number of measurements can be collected through dynamic test runs in a short span of time. This paper presents a systematic design of dynamic experiments (DoDEs) for identifiability of model parameters, which is achieved by persistently exciting the inputs signals. A pseudorandom binary sequence (PRBS) is designed as the input signal for DoDE due to its efficiency in obtaining sufficient spectral content. However, due to the long sequence size of the PRBS signal, a Schroeder-phase input signal, which is a multisine signal, is also designed. Tests for both types of signals are run in the Pilot Solvent Test Unit (PSTU) at the National Carbon Capture Center in Wilsonville, Alabama. The transient data are used to solve dynamic data reconciliation and parameter estimation problem. The estimated parameters are found to be not only superior to those estimated from using data collected from hundreds of steady-state test runs in a nonreactive (air–water) system, but the parameters could be estimated by using the dynamic data collected for about 24 h from the pilot plant for the MEA-H 2 O–CO 2 system.

CO2 capture↗

Organic Acid-Assisted Thermal Dehalogenation of Halide Salt Nuclear Wastes: From Waste Salts to Borosilicate Glass

Only a handful of high-halide salt waste forms have been demonstrated for vitrification-based immobilization strategies for halide-salt nuclear waste streams (e.g., pyroprocessing wastes, molten salt reactor wastes) and they all have low waste loading potential and most have low chemical durabilities for high-alkali streams. An alternative approach to direct salt immobilization is salt partitioning prior to waste form fabrication and one option for partitioning is halide removal (called dehalogenation). Removing the halogen fraction through dehalogenation can significantly reduce the waste volume required for disposal in the primary waste form. Furthermore, when dehalogenation is performed using organic acids, the dehalogenation reagent can decompose during high-temperature vitrification, reducing waste loading limitations in the waste form. In the current work, different organic acids (i.e., oxalic, formic, acetic, oxamic, and citric) were evaluated for dehalogenation efficiency of a simple chloride salt simulant (7.19% LaCl 3 , 53.77% LiCl, and 39.04% KCl, by mole) and a more complex chloride salt simulant called ERV3 (electrorefiner version 3) at 150 °C–300 °C and using H + /Cl – molar ratios of 1:1, 2:1, and 3:1. Additionally, a borosilicate glass waste form called TARS (or the average of refined specifications) was formulated, produced, and characterized for dehalogenated ERV3.

Amorphous materials↗

Extraction, Characterization, and Stability Studies of Bistriazinyl-Derived Carboxylic Acids

Separation of An(III) and Ln(III) ions will benefit the recycling of used nuclear fuel (UNF). For this purpose, many ligands have been tested over the years, and several separation processes have been successfully demonstrated on the laboratory scale. Current research aims at the development of new ligands that are built only with carbon, hydrogen, oxygen, and nitrogen (CHON), as they can be incinerated completely without secondary waste production. Here, we tested a new class of water-soluble ligands, the bistriazinyl-octa-carboxylic acids. One member, in particular, 2,6-bis-[5,6-di(3,4-dicarboxyphenyl)-1,2,4- triazin-3-yl]-pyridine (BTPOA), was found to be suitable for the selective separation of Am(III) and Cm(III) ions from Ln(III) ions and may act as a CHON alternative to its sulfonated analogue (SO 3 -Ph-BTP). BTPOA exhibited good extraction results and a high selectivity for Am(III) over Eu(III) ions. This ligand’s complexation of metal ions was further studied using potentiometric spectroscopy, as well as time-resolved laser-induced spectroscopy with Cm(III) and Eu(III) in aqueous HClO 4 and HNO 3 media. Conditional stability constants of each formed species were determined. In the HClO 4 system, Cm(III) formed three species (1:1, 1:2, and 1:3) through the stepwise addition of a single BTPOA molecule. On the other hand, in HNO 3 , Cm(III) formed two 1:2 complexes and one 1:3 complex, while the stepwise formation of three species was observed for Eu(III). The stability constants are comparable to the values for SO 3 -Ph-BTP. The radiolytic behavior of BTPOA was also investigated using electron pulse irradiation measurements to determine absolute rate coefficients (k) under ambient temperature conditions for the reaction of BTPOA with typical UNF reprocessing radical radiolysis products the hydrated electron (e aq − , k = (1.60 ± 0.02) × 10 10 M −1 s −1 ), the hydrogen atom (H • , k = (2.17 ± 0.03) × 10 9 M −1 s −1 ), and hydroxyl ( • OH, k = (6.95 ± 0.06) × 10 9 M −1 s −1 ) and nitrate (NO 3 • , k = (0.37 ± 0.02) × 10 7 M −1 s−1) radicals. These rate coefficients indicate that the radiolytic longevity of BTPOA should increase with HNO 3 concentration, owing to the consumption of e aq − /H • , by nitrate anions, and the replacement of • OH by the less reactive NO 3 • .

Diaz Gomez, Laura [Forschungszentrum Juelich (Germ↗

Computational Investigation of the Chemical Bond between An(III) Ions and Soft-Donor Ligands

The chemical bonding of actinide ions with arene and borohydride ligands is explored via quantum chemical methods to understand how the transuranium elements interact with softdonor ligands. Specifically, the [An(C 6 Me 6 )(BH 4 ) 3 ] complexes (An = U, Np, and Pu) and their reduced congeners are studied. Density functional theory (DFT) shows that the metal–ligand interactions in the neutral complexes are governed by electrostatic interactions. Both DFT and complete active space (CASSCF) results show that as one moves from U to Pu, the 5f-orbitals are stabilized leading to a poorer energy match with the ligand orbitals. This contributes to progressively weaker metal-arene and metal-borohydride interactions across the series due to a decrease in energy-driven covalency. A reduction in orbital contributions to bonding is obtained for the transuranium-arene interactions as well. Upon reduction, the arene is reduced, forming a δ-bond. This causes the An–arene distances to contract by 0.1–0.2 Å compared to the neutral complexes. The ground state is assigned as the intermediate-spin state where the arene radical is antiferromagnetically coupled to the metal-centered f-electrons in Np and Pu. On the other hand, the ferromagnetically and antiferromagnetically coupled states are close in energy in the uranium complex, but do not mix when spin– orbit coupling is included using a state-interaction approach (SO-CASPT2). The population of the CASSCF δ*-antibonding natural orbital increases from U to Pu consistent with the increased An–arene distances, weaker interactions, and decreasing covalency across the series. Although the An–B distance increases by ca. 0.06 Å upon reduction, both the neutral and reduced species involve an An(III)–borohydride bond and as such are qualitatively similar. The Np complexes can be assigned to have slightly weaker bonding than the uranium analogs but are overall “uranium-like”. The Pu complexes are predicted to have less covalent contributions to bonding in both the Pu–arene and Pu–borohydride interactions; however, the Pu–arene interaction is predicted to be particularly weak.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Rhenium-Sulfido and -Dithiolato Corroles: Reflections on Chalcophilicity

The high-temperature (~180 °C) reaction between free-base meso-triarylcorroles and Re 2 (CO) 10 , followed by exposure to PCl 3 and thiols (or elemental sulfur), affords rhenium-sulfido (ReS) corroles in 67-76% yields. The use of shorter reaction times, lower temperatures (~130 °C), and a dithiol (e.g., ethane-1,2-dithiol) also allows the isolation of rhenium-dithiolato corroles, presumptive intermediates on the path to ReS corroles. The ReS corroles exhibit high thermal stability and two reversible oxidations and reductions in their cyclic voltammograms, with redox potentials nearly identical to those observed for analogous ReO corroles. The electrochemical HOMO-LUMO gaps of the complexes, at 2.2 eV, are consistent with ligand-centered oxidation and reduction. The UV-vis spectra of ReS corroles, on the other hand, differ significantly from those of their ReO counterparts. Scalar-relativistic DFT calculations suggest that this difference reflects low-energy LUMO+2 and LUMO+3 levels, consisting of Re-S π-antibonding interactions; the ReO corroles, in contrast, exhibit a larger LUMO+1/LUMO+2 gap, as expected for a relatively classical Gouterman-type metalloporphyrin analogue. The high stability of ReS corroles is consistent with geochemists' view of rhenium as a moderately chalcophilic element (i.e., one that partitions into sulfide melts) as well as with a recent quantitative analysis of thiophilicity, which indicates that rhenium's oxophilicity and thiophilicity are essentially evenly balanced.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

In the Mix : A Workshop Merging Computational Chemistry and Electrochemistry Alongside Data Science

As chemistry expands to more complex and interdisciplinary areas, a new generation of diverse researchers must engage with science and learn effective cross-disciplinary collaboration and communication. To these ends, we designed and implemented In the Mix, a graduate student-led, two-day workshop for undergraduate students promoting collaborative science in the context of energy storage innovations. Here, the interactive workshop was designed for future and emerging researchers to gain hands-on experience with data science, computational chemistry, and electrochemistry techniques that are critical for developing materials for battery technologies. Participants also visited commercial renewable energy facilities to help them connect discovery-based research with industry and broader societal considerations. The workshop content and structure ensured that participants experienced the interrelatedness of the fields and understood the importance of collaborative research to yield scientific advances with real-world applications. An external team evaluated the workshop and participants’ perceptions of their experiences. While our research context was energy storage, the workshop goals and outcomes are applicable to other contexts. Interdisciplinary, experiential workshops are a key avenue to broadening participation in science and research, and the ideas presented here can be readily modified for other scientific contexts and/or incorporated as broader impact activities.

25 ENERGY STORAGE↗

Developing and Demonstrating a Lab Method for Quantifying Hair Exposure to Environmental Tobacco Smoke with a Forensic Perspective

Forensic science and chemistry curricula often lack hands-on experimental designs that effectively simulate the impact of environmental tobacco smoke (ETS) or other common elements found at crime scenes, such as marijuana, on trace forensic evidence. Hair, a critical form of trace evidence, offers unique advantages for assessing long-term exposure to environmental pollutants, including ETS. This study presents a novel, noninvasive forensic laboratory module designed to evaluate ETS exposure on various human hair types (untreated, dyed, and bleached). The experimental procedure involved controlled cigarette smoke exposure, followed by analysis using UV−visible spectroscopy, FTIR spectroscopy, and zeta potential measurements. Thirteen students participated in the three-week lab module (three sessions per week). Pre- and postlab assessments were conducted to evaluate learning outcomes. The prelab assessment focused on baseline knowledge of forensic hair analysis, as well as student expectations and confidence. The postlab assessment evaluated knowledge gained, technical insights, application of techniques, self-reflection, conceptual understanding, and practical skill development. This design helped students comprehend the effects of chemical treatments that significantly influence hair’s capacity to adsorb ETS residues by altering its physical and chemical properties. Integration of this experiment into the forensic chemistry curriculum led to measurable gains in student understanding, technical competency, and appreciation for realworld forensic applications. This method offers a valuable teaching and investigative tool for assessing individual ETS exposure in forensic contexts.

FTIR↗

A Theory of Ultrafast Charge Transfer Relaxation with Non-Innocent Solvent Molecules

In this work we revisit the photodynamics of tetracyanoethylene-hexamethylbenzene (TCNE- HMB), the molecular complex studied by Hochstrasser et al. [J. Chem. Phys. 100, 4797–4810, 1994] that has long challenged the applicability of Marcus’s theory of elec- tron transfer for predicting photochemical reactions. Using a novel black-box electronic structure algorithm (time-dependent density functional theory with one double, TD- DFT-1D) to efficiently run molecular dynamics that can treat charge recombination, we run ab initio surface hopping molecular dynamics and confirm that, for a polar solvent, charge recombination rates can be incredibly fast (indeed faster than the sol- vent relaxation time); for non-polar solvents, the rate is much slower. We demonstrate that, although Marcus theory cannot be directly applied, these nonequilibrium (and sometimes incredibly fast) photoexcited dynamics can be effectively explained within a two-state model without any evidence of a transition through a conical intersection. Most importantly, for this paradigmatic model system, we are able to identify two nuclear coordinates of interest (rather than the single coordinate predicted by Marcus or a full set of internal quantum modes studied by Bixon and Jortner): the solvent relaxation in the first shell (that strongly modulates the energies of the charge trans- fer state and differentiates time scales for relaxation) and a nuclear displacement in the TCNE-HMB complex arising from a handful of vibrations that induces non-Born Oppenheimer motion and eventually facilitates an abrupt electronic transition to the ground state. Altogether, these findings suggest a tractable generalization of Marcus theory for future simulations of photochemistry with non-innocent solvent environ- ments in the spirit of a Hamiltonian suggested by Stuchebrukhov (J. Chem. Phys. 107, 3821, 1997).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Reducing the Cost of CCSD Basis Set Extrapolation in Ab Initio Computational Thermochemistry

Here, a series of approximations to CCSD contributions in computational model chemistries is presented in the context of kcal mol –1 , kJ mol –1 , and 20 cm –1 theoretical predictions of total atomization energies, benchmarked within the HEAT+CH 4 test suite. A specific set of circumstances where MP2, without empirical scaling, may be used as an effective intermediate in the first two of these accuracy ranges was determined. However, SDQ-MP4, a method long used in pursuit of kcal mol –1 accuracy but relatively unstudied in the subchemical accuracy community, offers significant improvement over the quality of MP2 as a basis-set intermediate at significantly reduced cost compared to CCSD. Given this, we argue for SDQ-MP4 as the de facto CCSD basis-set intermediate in sub-chemical accuracy calculations when CCSD in a desired basis set becomes unaffordable. We additionally report on a “CBS-like” scheme, where MP2 and SDQ-MP4 are used in conjunction to create a “cheap” three-part approximation of large CCSD basis set limits. The data for the CCSD approximation schemes are organized in such a way that model chemistry developers can locate an analog of their current approach for the CCSD basis set limit and explore alternative intermediates that either decrease computational cost or increase computational accuracy. We also show, for a handful of molecules, that SDQ-MP4 shows promise as an effective basis-set intermediate for harmonic and fundamental frequency computations, allowing for zero-point corrections of nearly CCSD(T)/ANO1 quality using simple composite methods that only require CCSD(T)/ANO0.

Thorpe, James H. [Argonne National Laboratory (ANL↗