Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “kinetic parameters”

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 865 records · Page 48
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Lagrangian particle model for 3D simulation of pellets and SPI fragments in tokamaks

A 3D numerical model for the ablation of pellets and shattered pellet injection (SPI) fragments in tokamaks in the plasma disruption mitigation and fueling parameter space has been developed based on the Lagrangian particle code [R. Samulyak, X. Wang, H.-S. Chen, Lagrangian Particle Method for Compressible Fluid Dynamics, J. Comput. Phys., 362 (2018), 1-19]. The pellet code implements the low magnetic Reynolds number MHD equations, kinetic models for the electronic heating, a pellet surface ablation model, an equation of state that supports multiple ionization states, radiation, and a model for grad-B drift of the ablated material across the magnetic field. The Lagrangian particle algorithm is highly adaptive, capable of simulating a large number of fragments in 3D while eliminating numerical difficulties of dealing with the tokamak background plasma. The code has achieved good agreement with theory for spherically symmetric ablation flows. Axisymmetric simulations of neon and deuterium pellets in magnetic fields ranging from 1 to 6 Tesla have been compared with previous simulations using the FronTier code, and very good agreement has also been obtained. Furthermore, the main physics contribution of the paper is a detailed study of the influence of 3D effects, in particular grad-B drift, on pellet ablation rates and properties of ablation clouds. Smaller reductions of ablation rates in magnetic fields compared to axially symmetric simulations have been demonstrated because the ablated material is not confined to narrowing channels in the presence of grad-B drift. Contribution of various factors in the grad-B drift model has also been quantified.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Resistive wall mode stability and resonant field amplification in MAST high beta plasma

The n=1 (n is the toroidal mode number) resistive wall mode (RWM) instability and the resonant field amplification (RFA) due to a stable RWM response are numerically investigated for a MAST high-pressure plasma scenario, utilizing the MARS-F (Liu Y.Q. et al 2000 Phys. Plasmas 7 3681) and MARS-K (Liu Y.Q. et al 2008 Phys. Plasmas 15 112503) codes, for the purpose of understanding the potential RWM behavior in MAST-U plasmas. Ideal kink stability analysis identifies a target plasma, with parameters similar to that of the reference equilibrium reconstructed from a MAST high-pressure discharge, that accesses the RWM regime. The unstable n=1 RWM for the target equilibrium is subject to strong damping by the plasma toroidal flow and/or the drift kinetic effects from thermal particles. As a result, the mode is found to be stable under the experimental flow conditions, or even without flow stabilization if drift kinetic stabilization is included. The stability prediction is robust against variation of the assumed resistive wall minor radius. Active magneto-hydrodynamic (MHD) spectroscopy modeling, using the magnetic coils designed for controlling the edge localized modes (ELMs) in MAST-U as the antenna, shows strong resonant field amplification due to a stable RWM response in the target plasma. Maximal amplification, of 6-7 times larger than the applied vacuum field within the plasma, is obtained assuming a coil phasing of 120-180 degrees between the upper and lower rows of the ELM control coils. The MHD-kinetic hybrid model predicts generally higher RFA than the fluid model, in particular near the low-field and high-field sides of the plasma boundary surface. Furthermore, the MARS-F/K modeling thus shows that the RWM in MAST-U is subject to strong flow and kinetic stabilization, but can nevertheless be detected by active MHD spectroscopy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Quantitative surface studies of moiré intercalated Pb structures under Gr/SiC

Quantitative growth studies using surface diffraction and scanning tunneling microscopy are carried out to determine the intercalation conditions for Pb under Gr/SiC, with initial phase a mixture of buffer-layer and single-layer graphene. Intercalation is a kinetically complex, multistage process with intercalated Pb moving between different subsurface locations with annealing, eventually to the location of higher stability. The transfer is determined by monitoring the difference of a given spot intensity from its clean value, as a function of growth parameters. In addition, a broad, low-intensity (10×10) gr moiré phase is observed exclusively originating from intercalated Pb. Well-defined growth conditions are identified for the moiré spot to be observed reversibly. The (10×10) gr phase follows changes of the Pb(10) spot during growth and annealing experiments, instead of changes of the other Gr/SiC spots indicating preferred nucleation of Pb islands on the (10×10) gr domains. These structural results are discussed in connection to different models proposed for the (10×10) gr supercell and the subsurface Pb coverage. In conclusion, mass transport of subsurface Pb is unusually high, possibly because of the very high compression of the intercalated Pb, when compared to metallic Pb(111).

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Zero-field spin resonance in graphene with proximity-induced spin-orbit coupling

In this work we investigate collective spin excitations in graphene with proximity-induced spin-orbit coupling (SOC) of the Rashba and valley-Zeeman types, as it is the case, e.g., for graphene on transition-metal-dichalcogenide substrates. It is shown that, even in the absence of an external magnetic field, such a system supports collective spin modes, which correspond to coupled oscillations of the uniform and valley-staggered magnetizations. These modes can be detected via both zero-field electron spin resonance (ESR) and zero-field electric-dipole spin resonance (EDSR), with EDSR response arising solely due to Rashba SOC. We analyze the effect of electron-electron interaction within the Fermi-liquid kinetic equation and show that the interaction splits both the ESR and EDSR peaks into two. The magnitude of splitting and the relative weights of the resonances can be used to extract the spin-orbit coupling constants and many-body interaction parameters that may not be accessible by other methods.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Computational Modeling of Molten Salt Infiltration and Oxidation in Nuclear Graphite

Graphite is utilized as a moderator and reflector in advanced nuclear reactor designs due to its high thermal conductivity, neutron moderation properties, and resistance to radiation damage. However, its longterm performance and reliability are challenged by degradation mechanisms such as molten salt infiltration in molten salt reactors (MSRs) and oxidation in gas-cooled reactors (GCRs). These mechanisms can compromise the structural integrity and operational lifetime of graphite components, necessitating a more detailed assessment of their physical behavior. This report focuses on the development of computational models for molten salt infiltration and oxidation of graphite to aid the design and performance analysis of graphite components. For molten salt infiltration, a computational framework is developed that couples incompressible Navier-Stokes and phase-field model to simulate the penetration of molten salt into graphite?s interconnected pore structure. Initial model verification is performed using two-phase flows in two dimensions, demonstrating the models ability to capture fundamental physical behavior and agree with analytical solution. This framework is then applied to a realistic IG110 nuclear graphite , where a computed tomography extracted pore geometry is used to analyse the infiltration behavior of FLiNaK molten salt. This model provides insights into how the microstructure and other relevant parameters influence the transport pathways of molten salt into graphite, potentially offering a means to rapidly evaluate a graphite grade?s resistance to infiltration. For oxidation, the report details pore-scale mass and heat transport models, describing the diffusion of gases, reaction kinetics, and thermal effects. Additionally, this report highlights inconsistencies in the existing volume-averaged macroscopic model, particularly in upscaling of reaction kinetics and flux terms, and surface to volume transformations. These inconsistencies suggest that current formulations may not accurately capture the experimentally observed graphite oxidation process, highlighting the need for improved model development. This work advances the development of physics-based computational models for graphite degradation, contributing to improved predictive models for next-generation nuclear reactor designs. Future efforts will focus on refining the infiltration model to address non-physical behaviors and enhance its robustness. Additionally, for oxidation, further studies will employ the principles of volume averaging to rigorously derive the upscaled equations, potentially in collaboration with subject matter experts.

Computational Modeling of Molten Salt Infiltration↗
DOE OSTI2025

Delving into the depths of NGC 3783 with XRISM

We present our study of the X-Ray Imaging and Spectroscopy Mission (XRISM) observation of the Seyfert-1 galaxy NGC 3783. XRISM’s Resolve microcalorimeter has enabled, for the first time, a detailed characterization of the highly ionized outflows in this active galactic nucleus. Our analysis constrains their outflow and turbulent velocities, along with their ionization parameter (ξ) and column density (N H ). The high-resolution Resolve spectrum reveals a distinct series of Fe absorption lines between 6.4 and 7.8 keV, ranging from Fe XVIII to Fe XXVI. At lower energies (1.8−3.3 keV), absorption features from Si, S, and Ar are also detected. Our spectroscopy and photoionization modeling of the time-averaged Resolve spectrum uncovers six outflow components, five of which exhibit relatively narrow absorption lines with outflow velocities ranging from 560 to 1170 km s −1 . In addition, a broad absorption feature is detected, which is consistent with Fe XXVI outflowing at 14 300 km s −1 (0.05 c). The kinetic luminosity of this component is 0.8−3% of the bolometric luminosity. Our analysis of the Resolve spectrum shows that more highly ionized absorption lines are intrinsically broader than those of lower-ionization species, indicating that the turbulent velocity of the six outflow components (ranging from 0 to 3500 km s −1 ) increases with ξ. Furthermore, we find that the column density (N H ) of the outflows generally declines with the ionization parameter up to log ξ = 3.2 but rises beyond this point, suggesting a complex ionization structure. The absorption profile of the Fe XXV resonance line is intriguingly similar to UV absorption lines (Lyα and C IV) observed by the Hubble Space Telescope, from which we infer that the outflows are clumpy in nature. Our XRISM/Resolve results from lower- and higher-ionization regimes support a “hybrid wind” scenario in which the observed outflows have multiple origins and driving mechanisms. We explore various interpretations of our findings within active galactic nucleus wind models.

X-rays: galaxies↗
DOE OSTI2022

A workflow for automatic generation and efficient refinement of individual pressure-dependent networks

Manual chemical kinetic model construction often requires modelists to at least implicitly guess all possible decay pathways and their relative fluxes for all important species. We present a workflow and associated tool for automatic generation and re finement of pressure-dependent networks that should enable modelists to efficiently and comprehensively identify decay pathways and estimate respective parameters for chemical species. This tool combines the capabilities of the Reaction Mechanism Generator (RMG) software to generate possible reaction paths, determine thermochemistry, approximate rate coefficients and estimate frequencies with the capabilities of the Arkane software to use quantum chemical parameters to compute thermochemistry and pressure dependent rates. A flux-based algorithm is used to decide which isomers to add to the network. Isomers added to the network are reacted to form other channels. When enough of the flux is accounted for by the isomers and bimolecular product channels, the generation process terminates to yield a comprehensive network. Network sensitivity analysis is applied to the network to identify the most important wells and barriers that can be refined using quantum chemistry calculations. Iterative refinement can then be used to achieve accuracy. The comprehensive network can be easily reduced using energy and flux-based algorithms to the most important channels and isomers in the network.

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

A Bezier Curve Informed Melt Pool Geometry to Model Additive Manufacturing Microstructures Using SPPARKS

Additive manufacturing is a transformative technology with the potential to manufacture designs which traditional subtractive machining methods cannot. Additive manufacturing offers fast builds at near final desired geometry; however, material properties and variability from part to part remain a challenge for certification and qualification of metallic components. AM induced metallic microstructures are spatially heterogeneous and highly process dependent. Engineering properties such as strength and toughness are significantly affected by microstructure morphologies resulting from the manufacturing process Linking process parameters to microstructures and ultimately to the dynamic response of AM materials is critical to certifying and qualifying AM built parts and components and improving the performance of AM materials. The AM fabrication process is characterized by building parts layer by layer using a selective laser melt process guided by a computer. A laser selectively scans and melts metal according to a designated geometry. As the laser scans, metal melts, fuses, and solidifies forming the final geometry in a layerwise fashion. As the laser heat source moves away, the metal cools and solidifies forming metallic microstructures. This work describes a microstructure modeling application implemented in the SPPARKS kinetic Monte Carlo computational framework for simulating the resulting microstructures. The application uses Bzier curves and surfaces to model the melt pool surface and spatial temperature profile induced by moving the laser heat source; it simulates the melting and fusing of metal at the laser hot spot and microstructure formation and evolution when the laser moves away. The geometry of the melt pool is quite flexible and we explore effects of variances in model parameters on simulated microstructures.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Combining Computational Modeling with Reaction Kinetics Experiments for Elucidating the In Situ Nature of the Active Site in Catalysis

Microkinetic modeling based on density functional theory (DFT) derived energetics is important for addressing fundamental questions in catalysis. The quantitative fidelity of microkinetic models (MKMs), however, is often insufficient to conclusively infer the mechanistic details of a specific catalytic system. Here, this can be attributed to a number of factors such as an incorrect model of the active site for which DFT calculations are performed, deficiencies in the hypothesized reaction mechanism, inadequate consideration of the surface environment under reaction conditions, and intrinsic errors in the DFT exchange-correlation functional. Despite these limitations, we aim at developing a rigorous understanding of the reaction mechanism and of the nature of the active site for heterogeneous catalytic chemistries under reaction conditions. By achieving parity between experimental and modeling outcomes through robust parameter estimation and by ensuring coverage-consistency between DFT calculations and MKM predictions, it is possible to systematically refine the mechanistic model and, thereby, our understanding of the catalytic active site in situ.

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

Modeling discontinuous dynamic recrystallization containing second phase particles in magnesium alloys utilizing phase field method

Pre-existing second phase particles can induce particle stimulated nucleation (PSN) and pinning effect on grain boundaries during DRX process of magnesium alloys. The interaction among pinning effect, PSN mechanism and strain-induced grain boundary migration (SIBM) nucleation mechanism imposes challenge in quantitative study of microstructure evolution. A phase field model that describes discontinuous dynamic recrystallization process considering second phase particles (PF-DDRX-SP) which mainly distribute along grain boundaries is developed. The second phase order parameter with diffusion interface is introduced to the free energy function. With considering the pinning effect of second phase particles on grain boundaries, the grain boundary energy couples the interaction between second phase and grains. The effect of second phase on dislocation density evolution is considered in the model, and the PSN and SIBM mechanism are both executed by limiting the nucleation position to the grain boundary and the second phase interface. For the validation, the initial topology consistent with initial microstructure is applied to the simulation of DRX kinetics and flow behavior of AZ80 magnesium alloy, and great reliability and accuracy of developed model is indicated. Furthermore, through the model, the microstructure evolution with different particle size and with different volume fraction of round particle during DRX process were simulated. Finally, the results show that particles with higher dispersion reduce the number of DRX nuclei, which indicate that second phase particles suppress the nucleation at grain boundaries because of the pinning effect. More importantly, the proposed model could also quantitatively predict the relationships between the parameters of second phase particles and DRX behaviors, and enable to optimize the initial second phase structure in a uniform grain structure during thermomechanical process.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Synthesis of UO 2 nanoparticles via coulometric titration: influence of electrolytes and ligands on synthesized particle properties

Here, this study investigates the controlled synthesis of uranium oxide nanoparticles (UO 2 NPs) via coulometric titration under ambient conditions, focusing on the impact of electrolyte anions, salt concentrations, and strongly complexing ligands on particle formation and properties. Using dilute solutions of mineral acids (HNO 3 , HCl, and HClO 4 ), we demonstrate that the choice of electrolyte anion affects particle size, polydispersity, and surface charge. Particle characterization using dynamic light scattering and transmission electron microscopy shows that particles synthesized in perchlorate are largest and show a high degree of polydispersity. In comparison, particles synthesized in chloride are smaller and more uniform in size. Synthesis in nitrate yields a wide variety of particle sizes, with the major fraction (>65 %) having a smaller size than that obtained in the other two electrolytes. The presence of more strongly coordinating ligands such as sulfate and acetate modulate hydrolysis and condensation reactions, with sulfate suppressing nanoparticle formation across a wide concentration range ([SO 4 2- ] > 5 mM) and acetate enabling stable colloidal suspensions at [HOAc] ≤ 0.1 M. Synthesis in concentrated electrolytes (2 M NaNO 3 , NaCl, or NaClO 4 ) accelerates reaction kinetics but introduces challenges, as particles showed increased polydispersity and aggregation, and were more prone to oxidation. Electrolyte effects on actinide oxide nanoparticle formation are discussed and a short comparison to established nanoparticle syntheses is drawn. This work underscores the importance of tailoring synthesis parameters to achieve desired nanoparticle properties, providing valuable insights for optimizing UO 2 NP production for various applications.

Actinide hydrolysis↗
DOE OSTI2024

Prospects for a high-field, compact break-even axisymmetric mirror (BEAM) and applications

This paper explores the feasibility of a break-even-class mirror referred to as BEAM (break-even axisymmetric mirror): a neutral-beam-heated simple mirror capable of thermonuclear-grade parameters and $Q\sim 1$ conditions. Compared with earlier mirror experiments in the 1980s, BEAM would have: higher-energy neutral beams, a larger and denser plasma at higher magnetic field, both an edge and a core and capabilities to address both magnetohydrodynamic and kinetic stability of the simple mirror in higher-temperature plasmas. Axisymmetry and high-field magnets make this possible at a modest scale enabling a short development time and lower capital cost. Such a $Q\sim 1$ configuration will be useful as a fusion technology development platform, in which tritium handling, materials and blankets can be tested in a real fusion environment, and as a base for development of higher- $Q$ mirrors.

Physics↗
DOE OSTI2021

Giant field-like torque by the out-of-plane magnetic spin Hall effect in a topological antiferromagnet

Abstract Spin-orbit torques (SOT) enable efficient electrical control of the magnetic state of ferromagnets, ferrimagnets and antiferromagnets. However, the conventional SOT has severe limitation that only in-plane spins accumulate near the surface, whether interpreted as a spin Hall effect (SHE) or as an Edelstein effect. Such a SOT is not suitable for controlling perpendicular magnetization, which would be more beneficial for realizing low-power-consumption memory devices. Here we report the observation of a giant magnetic-field-like SOT in a topological antiferromagnet Mn 3 Sn, whose direction and size can be tuned by changing the order parameter direction of the antiferromagnet. To understand the magnetic SHE (MSHE)- and the conventional SHE-induced SOTs on an equal footing, we formulate them as interface spin-electric-field responses and analyzed using a macroscopic symmetry analysis and a complementary microscopic quantum kinetic theory. In this framework, the large out-of-plane spin accumulation due to the MSHE has an inter-band origin and is likely to be caused by the large momentum-dependent spin splitting in Mn 3 Sn. Our work demonstrates the unique potential of antiferromagnetic Weyl semimetals in overcoming the limitations of conventional SOTs and in realizing low-power spintronics devices with new functionalities.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Ultrafast conformational dynamics of Rydberg-excited N -methyl piperidine

We have observed the ultrafast conformational dynamics of electronically excited N-methyl piperidine (NMP) using time-resolved Rydberg fingerprint spectroscopy. Optical excitation at various wavelengths ranging from 212 nm to 229 nm leads to the 3s or 3p Rydberg states and induces coherent oscillatory motions with periods of about 700 fs. Furthermore, these coherent motions survive the internal conversion from 3p to 3s but then dephase on a time scale of a few oscillations. Intramolecular vibrational energy redistribution on a picosecond time scale leads to an equilibrium between two conformeric structures that are separated in binding energy by 0.09 eV. Model calculations using the DFT-SIC method are in excellent agreement with the experiments and identify the conformers as the chair and twist structures of NMP. The analysis of the equilibrium parameters at long time delays as a function of excitation wavelength allows for the extraction of thermodynamic parameters for the conformeric transformation. We derive an enthalpy of the chair to twist reaction in the 3s excited state of 62 meV with an entropy of 19.70 J mol –1 K –1 . An activation energy of 276 meV is also obtained with a kinetic model.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Hybrid simulations of sub-cyclotron compressional and global Alfvén eigenmode stability in spherical tokamaks

A comprehensive numerical study has been conducted in order to investigate the stability of beam-driven, sub-cyclotron-frequency compressional Alfvén eigenmodes (CAEs) and global Alfvén eigenmodes (GAEs) in low-aspect-ratio plasmas for a wide range of beam parameters. The presence of CAEs and GAEs has previously been linked to anomalous electron temperature profile flattening at high beam powers in NSTX experiments, prompting a further examination of the conditions necessary for their excitation. Linear simulations have been performed with the hybrid MHD–kinetic initial value code HYM in order to capture the general Doppler-shifted cyclotron resonance that drives the modes. Additionally, three distinct types of modes were found in the simulations—co-CAEs, cntr-GAEs, and co-GAEs—with differing spectral and stability properties. The simulations revealed that unstable GAEs are more ubiquitous than unstable CAEs, which is consistent with experimental observations, as they are excited at lower beam energies and generally have larger growth rates. Local analytic theory is used to explain key features of the simulation results, including the preferential excitation of different modes based on beam injection geometry and the growth rate dependence on the beam injection velocity, critical velocity, and degree of velocity space anisotropy. The background damping rate is inferred from the simulations and analytically estimated for relevant sources absent from the simulation model, indicating that co-CAEs are closer to marginal stability than modes driven by the cyclotron resonances.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Direct Numerical Simulation of Transitional and Turbulent Flows Over Multi-Scale Surface Roughness—Part II: The Effect of Roughness on the Performance of a High-Pressure Turbine Blade

Abstract Turbine blades generally present surface roughness introduced in the manufacturing process or caused by in-service degradation, which can have a significant impact on aero-thermal performance. A better understanding of the fundamental physical mechanisms arising from the interaction between the roughness and the turbine flow at engine-relevant conditions can provide insights for the design of blades with improved efficiency and longer operational life. To this end, a high-fidelity numerical framework combining a well-validated solver for direct numerical simulation and a second-order accurate immersed boundary method is employed to predict roughness-induced aero-thermal effects on an LS89 high-pressure turbine (HPT) blade at engine-relevant conditions. Different amplitudes and distributions of surface roughness are investigated and a reference smooth-blade simulation under the same flow conditions is conducted for comparison. Roughness of increasing amplitude progressively shifts the blade suction side boundary layer transition upstream, producing larger values of the turbulent kinetic energy and higher total wake losses. The on-surface data-capturing capabilities of the numerical framework provide direct measurements of the heat flux and the skin friction coefficient, hence offering quantitative information between the surface topology and engineering-relevant performance parameters. This work may provide a benchmark for future numerical studies of turbomachinery flows with roughness.

Engineering↗
DOE OSTI2024

$\overline{TKE}$ Parameterization and $\bar{v}$ Uncertainty Analysis for CGMF

Previous work was performed on tuning CGMF parameters for 235 U, 238 U, and Plutonium isotopes. Now work is being done to tune minor uranium isotopes. However, uranium isotopes like 232 U and 236 U have almost no experimental data. We are applying cross-isotope models to extrapolate and tune CGMF on isotopes that lack experimental data. There exist several internal CGMF physics quantities that affect the output of CGMF—multi-chance fission probability, excitation energy sharing, spin-cutoff factor, spin scaling, and fragment total kinetic energy to name a few. The mean fragment total kinetic energy, $\overline{TKE}$, is particularly interesting because of its strong anti-correlation with $\bar{v}$. We are most interested in the mean fragment total kinetic energy before neutron emissions. $\overline{TKE}$ is assumed to be pre-neutron emission unless otherwise stated. Currently in CGMF, the $\overline{TKE}$ model for 233,234,235,238 U are tuned independently to reproduce ν for the associated isotopes. In this report, we will tune a cross-isotope $\overline{TKE}$ model to experimental $\overline{TKE}$ data for 232,233,234,235,236,238 U. Because of the unreliable and sparse nature of $\overline{TKE}$ experimental data, future work will use more reliable experimental $\bar{v}$ data to infer the $\overline{TKE}$ model (and likely other internal CGMF parameters) for uranium isotopes. Such work has been performed previously using a sensitivity analysis and Kalman filter methods.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2025

Revealing the Influence of Binding Motifs on Electron Transfer and Recombination Kinetics for CdSe Quantum Dots Functionalized with a Modified Viologen

Anchoring of molecules to the surfaces of semiconductor nanocrystals (NCs) presents an opportunity to leverage the precise synthetic tunability of molecular function and the remarkable light harvesting properties of NCs to drive photochemical reactions. However, charge transfer between the two species depends not only on the energy level alignments but also on the details of their binding interactions, which are difficult to probe. Here, in this work, we characterize the binding between CdSe quantum dots (QDs) and a new phosphonated derivative of the electron acceptor methyl viologen, designed to attach to the QD surface via the phosphonate group. We use isothermal titration calorimetry to probe the thermodynamics of the QD–molecule interaction and use the parameters determined therein to analyze transient absorption spectroscopy measurements of forward and back electron transfer from QDs to the viologen. We find that the ligand-like phosphonate binding leads to an electron-transfer rate constant that is 3 orders of magnitude smaller than that for the face-on binding of the bipyridine ring of methyl viologen. Back electron transfer is also significantly slower in the derivative. Interestingly, a minor fraction of the phosphonated derivative also binds in the face-on configuration, with similar forward and back electron transfer kinetics as methyl viologen. Numerical simulations show that the ligand-like binding will lead to significantly improved quantum yields of photocatalysis over a wide range of reaction rates. By independently characterizing binding thermodynamics and charge transfer kinetics, this work reveals how the complexities underlying electron transfer at the NC–molecule interface determine photocatalytic outcomes. This work also represents a step toward controlling forward and back electron transfer kinetics via rational molecular design.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗