SEARCH · Search NASA

Results for “kinetics and energetics”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 217 records · Page 12

All sources NASA NTRS DOE OSTI

NASA NTRS2003

The Mass of a Solar Quiescent Prominence

This paper follows up on our recent paper on the role of prominence mass in the storage of magnetic energy for driving a coronal mass ejection (CME). The previous paper erroneously rejected a set of sheet- prominence solutions, the recovery of which allows for a simple theoretical estimate of the mass of a quiescent prominence. For coronal fields of 5-10 G, these hydromagnetic solutions suggest that a prominence mass of (1-26) x 10(exp 6) g is needed to hold detached magnetic fields of intensity comparable to the coronal fields in an unbounded atmosphere such that the global magnetic field is energetically able to spontaneously open up and still have enough energy to account for the kinetic and gravitational potential energies carried away in a CME. This simple result is discussed in relation to observed prominence magnetic field intensities, densities, and masses, pointing to the relevance of such observations to the question of magnetic energy storage in the solar corona.

Low, B. C.↗

NASA NTRS1993

Dispersed ion structures at the poleward edge of the auroral oval: Low-altitude observations and numerical modeling

We have compared the AUREOL 3 (A3) observations of auroral ion precipitation, particularly ion beams, with the results from the global kinetic model of magnetotail plasma of Ashour-Abdalla et al. (1993). We have identified 101 energetic keV H(+) velocity dispersed precipitating ion structures (VDIS) with fluxes above 10(exp -3) ergs./sq cm./s in the A3 record between the end of 1981 and mid-1984. These beams display a systematic increase in energy with increasing latitude and were observed in a narrow region within less than 1 deg in latitude of the polar cap boundary. The VDIS are the most distinctive feature in the auroral zone of the plasma sheet boundary layer. We report first on a statistical analysis of the possible ralationships between magnetic activity or substorm phase and the VDIS properties. Our particle simulations of the precipitating ions have been extended by using a series of modified versions of the Tsyganenko (1989) magnetic field model and by varying the cross-magnetosphere electric field. In the simulations, plasma from a mantle source is subject to strong nonlinear acceleration, forming beams which flow along the PSBL. Only 3 to 4% of these beams precipitate into the ionosphere to form the VDIS while the majority return to the equatorial plane after mirroring and form the thermalized central plasma sheet. The final energy and the dispersion of the beams in the model depend on the amplitude of the cross-tail electric field. Two unsual observations of low-energy (less than 5 keV) O(+) VDIS, shifted by 4 deg 5 deg in invariant latitude equatorward of H(+) VDIS are analyzed in detail. The sparsity of such O(+) events and the absence of the changes in the flux and frequency of occurrence indicate a solar wind origin for the plasma. Finally, large-scale kinetic modeling, even with its simplifications and assumptions (e.g., static magnetic field, solar wind source), reproduces low-altitude auroral ion features fairly well; it may therefore be presented as an appropriate framework into which data on energization and transport of the hot plasma, obtained in the equatorial plane, could be inserted in the near future.

Bosqued, J. M.↗

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 OSTI2024

TATB thermal decomposition: Expanding the molecular profile with cryo-focused pyrolysis GC-MS

Understanding the molecular composition of high explosives during thermal decomposition is vital for predicting the sensitivity, safety, and performance of explosive materials. The thermal decomposition of 1,3,5-triamino-2,4,6-trinitrobenzene (TATB) has been linked to the formation of furazans through a series of dehydration reactions of the NO2 and NH2 groups on the phenyl ring, along with breakdown into small molecules (≤120 amu). Molecular identification of compounds formed in this transformation of the furazans to light gases has been lacking. To address this, we have applied a pseudo-confined sampling system in a cryo-focused pyrolysis gas chromatography-mass spectrometry (pyGC-MS) system to molecularly identify these intermediates. By design, sublimation of TATB, which has complicated MS analyses of thermal degradation, was significantly reduced and additional compounds were identified with potential structural information. In addition to the known furazan compounds, one of these compounds forms from the loss of oxygen from benzo-trifurazan (F3) and produces an open ring structure that may be the first step in the formation of lower molecular weight furazan breakdown products. The loss of a nitro group from benzo-monofurazan (F1) was also discovered and implicates the formation of oxidizing NO2 gas in the thermal decomposition mechanism. So these findings are vital for understanding the proper heat flow from energetic materials on a molecular level, necessary when measuring enthalpy and developing decomposition models based on kinetic parameters.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

The theory of kinetic effects on resistive wall mode stability in tokamaks

Tokamak fusion plasmas benefit from high pressures but are then susceptible to modes of instability. These magnetohydrodynamic (MHD) modes are macroscopic distortions of the plasma, but certain collective motions of individual particles can provide stabilizing effects opposing them. The presence of a resistive wall slows the mode growth, converting a kink to a resistive wall mode (RWM). A kinetic MHD model includes Maxwell's equations, ideal MHD constraints, and kinetic effects included through the pressure tensor, calculated with the perturbed drift-kinetic distribution function of the particles. The kinetic stabilizing effects on the RWM arise through resonances between the plasma rotation and particle drift motions: precession, bounce, and transit. A match between particle motions and the mode allows efficient transfer of energy that would otherwise drive the growth of the mode, thus damping the growth. The first approach to calculating RWM stability is to write a set of equations for the complex mode frequency in terms of known quantities and then to solve the system. The “energy principle” approach, which has the advantage of clarity in distinguishing the various stabilizing and destabilizing effects, is to change the force balance equation into an equation in terms of changes of kinetic and potential energies, and then to write a dispersion relation for the mode frequency in terms of those quantities. These methods have been used in various benchmarked codes to calculate kinetic effects on RWM stability. Importantly, the theory has illuminated the important roles of plasma rotation, energetic particles, and collisions in RWM stability.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

Fast nonthermal processes in pulsed laser deposition

Pulsed laser deposition (PLD) is typically used to grow epitaxial thin films of quantum materials. In this research, we use in situ x-ray scattering to study homoepitaxy of SrTiO 3 by energetic deposition (e-PLD) versus PLD thermalized by a He background gas (th-PLD). Energetic PLD suppresses the lateral growth of two-dimensional islands, which suggests that particles with kinetic energies of ~ 100 eV break up smaller islands. Fast interlayer transport occurs for th-PLD as well as e-PLD, implying a process operating on submicrosecond time scales for incident particles with kinetic energies below 10 eV.

36 MATERIALS SCIENCE↗

NASA NTRS2000

High-LET radiation-induced aberrations in prematurely condensed G2 chromosomes of human fibroblasts

PURPOSE: To determine the number of initial chromatid breaks induced by low- or high-LET irradiations, and to compare the kinetics of chromatid break rejoining for radiations of different quality. MATERIAL AND METHODS: Exponentially growing human fibroblast cells AG1522 were irradiated with gamma-rays, energetic carbon (290MeV/u), silicon (490MeV/u) and iron (200 and 600 MeV/u). Chromosomes were prematurely condensed using calyculin A. Chromatid breaks and exchanges in G2 cells were scored. PCC were collected after several post-irradiation incubation times, ranging from 5 to 600 min. RESULTS: The kinetics of chromatid break rejoining following low- or high-LET irradiation consisted of two exponential components representing a rapid and a slow time constant. Chromatid breaks decreased rapidly during the first 10min after exposure, then continued to decrease at a slower rate. The rejoining kinetics were similar for exposure to each type of radiation. Chromatid exchanges were also formed quickly. Compared to low-LET radiation, isochromatid breaks were produced more frequently and the proportion of unrejoined breaks was higher for high-LET radiation. CONCLUSIONS: Compared with gamma-rays, isochromatid breaks were observed more frequently in high-LET irradiated samples, suggesting that an increase in isochromatid breaks is a signature of high-LET radiation exposure.

Non-NASA Center↗

DOE OSTI2026

The Effect of Large-scale Magnetic Fluctuations on Energetic Particle Profiles across Interplanetary Shocks

In situ observations of energetic particles associated with interplanetary (IP) shocks, known as energetic storm particle (ESP) events, often show time-intensity profiles quite different from predictions of classical diffusive shock acceleration theory. We use numerical simulations, including test-particle simulations and hybrid simulations (with fluid electrons and kinetic protons) for shocks propagating into a turbulent magnetic field to study the ESP intensity-time profiles across a strong IP shock. We find that several types of energetic particle intensity profiles similar to in situ observations can be produced in our simulations. The peak of energetic particle count is often near, but can be shifted from, the locations of the shock front, with the peak usually being observed downstream of the shock. These findings may help understand particle acceleration at both traveling IP shocks and shocks in other heliospheric and astrophysical environments.

Lo, Ming-Yuan [New Mexico Consortium (NMC), Los Al↗

DOE OSTI2021

Hydrodechlorination of 1,2-Dichloroethane on Platinum Catalysts: Insights from Reaction Kinetics Experiments, Density Functional Theory, and Microkinetic Modeling

Catalytic hydrodechlorination is a promising strategy for treating industrial 1,2-dichloroethane wastes, for which Pt and Pt-based alloy catalysts are widely used. Here, we performed a detailed mechanistic study for 1,2-dichloroethane hydrodechlorination on Pt using a synergistic approach combining density functional theory (DFT) calculations, reaction kinetics experiments, and microkinetic modeling. Using planewave DFT calculations, we evaluated the reaction energy and activation energy barrier of each elementary step involved in the reaction network on Pt(111). The calculated energetics were then incorporated into a comprehensive mean-field microkinetic model accounting for a total of 65 elementary steps. The model-predicted reaction rates were compared with the results from our reaction kinetics experiments on SiO2-supported Pt catalysts. Our results indicated that the hydrodechlorination of 1,2-dichloroethane on Pt(111) starts with a H-removal step; then, it proceeds through a sequence of alternating dechlorination and dehydrogenation steps until vinylidene (CH 2 C*) is formed; finally, CH 2 C* is hydrogenated to the final product, ethane, sequentially via vinyl (CH 2 CH*), ethylene, and ethyl (CH 3 CH 2 *) intermediates. After model parameter adjustments, we achieved good agreement between our theoretical model and experimental results; the adjustments to the calculated parameters are consistent with the typically anticipated coverage effects. Furthermore, our study offers valuable mechanistic insights, which are useful for improving catalysts for this chemistry.

1,2-dichloroethane↗

DOE OSTI2023

Machine learning-based ethylene and carbon monoxide estimation, real-time optimization, and multivariable feedback control of an experimental electrochemical reactor

Electrochemical reduction of CO 2 gas is a novel CO 2 utilization technique that has the potential to mitigate the global climate crisis caused by anthropogenic CO 2 emissions, and enable the large-scale storage of energy generated from renewable sources in the form of carbon-based chemicals and fuels. However, due to the complexity of the electrochemical reactions, the explicit first-principles models for CO2 reduction are not available yet, and there has been a limited effort to develop process modeling, optimization and control of CO 2 electrochemical reactors. To this end, a rotating cylinder electrode (RCE) reactor has been constructed at UCLA to understand the mass transfer and reaction kinetics effects separately on the productivity. In the RCE reactor, the applied potential strongly influences the reaction energetics and the electrode rotation speed affects the hydrodynamic boundary layer and modifies the film mass transfer coefficient, which involves convective and diffusive transport. Further, the present work aims to develop a multi-input multi-output (MIMO) control scheme for the RCE reactor that integrates techniques from artificial and recurrent neural network modeling, nonlinear optimization, and process controller design. Specifically, production rates of two products from the experimental reactor, ethylene and carbon monoxide, are controlled by manipulating two inputs, applied potential and catalyst rotation speed. Process dynamics and controllability are analyzed, a feedback control strategy is designed and the controllers are tuned accordingly. The experimental electrochemical cell is employed to gather data for process modeling and implement the multivariable control system. Finally, the experimental results are presented which demonstrate excellent closed-loop performance by the control system and regulation of the outputs at three different set-points including an economically-optimal set-point.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS2013

Galactic Cosmic Ray Event-Based Risk Model (GERM) Code

This software describes the transport and energy deposition of the passage of galactic cosmic rays in astronaut tissues during space travel, or heavy ion beams in patients in cancer therapy. Space radiation risk is a probability distribution, and time-dependent biological events must be accounted for physical description of space radiation transport in tissues and cells. A stochastic model can calculate the probability density directly without unverified assumptions about shape of probability density function. The prior art of transport codes calculates the average flux and dose of particles behind spacecraft and tissue shielding. Because of the signaling times for activation and relaxation in the cell and tissue, transport code must describe temporal and microspatial density of functions to correlate DNA and oxidative damage with non-targeted effects of signals, bystander, etc. These are absolutely ignored or impossible in the prior art. The GERM code provides scientists data interpretation of experiments; modeling of beam line, shielding of target samples, and sample holders; and estimation of basic physical and biological outputs of their experiments. For mono-energetic ion beams, basic physical and biological properties are calculated for a selected ion type, such as kinetic energy, mass, charge number, absorbed dose, or fluence. Evaluated quantities are linear energy transfer (LET), range (R), absorption and fragmentation cross-sections, and the probability of nuclear interactions after 1 or 5 cm of water equivalent material. In addition, a set of biophysical properties is evaluated, such as the Poisson distribution for a specified cellular area, cell survival curves, and DNA damage yields per cell. Also, the GERM code calculates the radiation transport of the beam line for either a fixed number of user-specified depths or at multiple positions along the Bragg curve of the particle in a selected material. The GERM code makes the numerical estimates of basic physical and biophysical quantities of high-energy protons and heavy ions that have been studied at the NASA Space Radiation Laboratory (NSRL) for the purpose of simulating space radiation biological effects. In the first option, properties of monoenergetic beams are treated. In the second option, the transport of beams in different materials is treated. Similar biophysical properties as in the first option are evaluated for the primary ion and its secondary particles. Additional properties related to the nuclear fragmentation of the beam are evaluated. The GERM code is a computationally efficient Monte-Carlo heavy-ion-beam model. It includes accurate models of LET, range, residual energy, and straggling, and the quantum multiple scattering fragmentation (QMSGRG) nuclear database.

Cucinotta, Francis A.↗

NASA NTRS2020

Nuclear Fusion Reactions in Deuterated Metals

Nuclear fusion reactions of D-D are examined in an environment comprised of high density cold fuel embedded in metal lattices in which a small fuel portion is activated by hot neutrons. Such an environment provides for enhanced screening of the Coulomb barrier due to conduction and shell electrons of the metal lattice, or by plasma induced by ionizing radiation (γ quanta). We show that neutrons are far more efficient than energetic charged particles, such as light particles (e−, e+) or heavy particles (p, d, α) in transferring kinetic energy to fuel nuclei (D) to initiate fusion processes. It is well-known that screening increases the probability of tunneling through the Coulomb barrier. Electron screening also significantly increases the probability of large- versus small-angle Coulomb scattering of the reacting nuclei to enable subsequent nuclear reactions via tunneling. This probability is incorporated into the astrophysical factor S(E). Aspects of screening effects to enable calculation of nuclear reaction rates are also evaluated, including Coulomb scattering and localized heating of the cold fuel, primary D-D reactions, and subsequent reactions with both the fuel and the lattice nuclei. The effect of screening for enhancement of the total nuclear reaction rate is a function of multiple parameters including fuel temperature and the relative scattering probability between the fuel and lattice metal nuclei. Screening also significantly increases the probability of interaction between hot fuel and lattice nuclei increasing the likelilhood of Oppenheimer-Phillips processes opening a potential route to reaction multiplication. We demonstrate that the screened Coulomb potential of the target ion is determined by the nonlinear Vlasov potential and not by the Debye potential. In general, the effect of screening becomes important at low kinetic energy of the projectile. We examine the range of applicability of both the analytical and asymptotic expressions for the well-known electron screening lattice potential energy Ue, which is valid only for E >> Ue (E is the energy in the center of mass reference frame). We demonstrate that for E ≤ Ue, a direct calculation of Gamow factor for screened Coulomb potential is required to avoid unreasonably high values of the enhancement factor f (E) by the analytical—and more so by the asymptotic—formulas.

Vladimir Pines↗

NASA NTRS2008

MESSENGER: Exploring Mercury's Magnetosphere

The MESSENGER mission to Mercury offers our first opportunity to explore this planet's miniature magnetosphere since Mariner 10's brief fly-bys in 1974-5. Mercury's magnetosphere is unique in many respects. The magnetosphere of Mercury is the smallest in the solar system with its magnetic field typically standing off the solar wind only - 1000 to 2000 km above the surface. For this reason there are no closed dri-fi paths for energetic particles and, hence, no radiation belts; the characteristic time scales for wave propagation and convective transport are short possibly coupling kinetic and fluid modes; magnetic reconnection at the dayside magnetopause may erode the subsolar magnetosphere allowing solar wind ions to directly impact the dayside regolith; inductive currents in Mercury's interior should act to modify the solar In addition, Mercury's magnetosphere is the only one with its defining magnetic flux tubes rooted in a planetary regolith as opposed to an atmosphere with a conductive ionosphere. This lack of an ionosphere is thought to be the underlying reason for the brevity of the very intense, but short lived, approx. 1-2 min, substorm-like energetic particle events observed by Mariner 10 in Mercury's magnetic tail. In this seminar, we review what we think we know about Mercury's magnetosphere and describe the MESSENGER science team's strategy for obtaining answers to the outstanding science questions surrounding the interaction of the solar wind with Mercury and its small, but dynamic magnetosphere.

Slavin, James A.↗

NASA NTRS2005

MESSENGER: Exploring Mercury's Magnetosphere

The MESSENGER mission to Mercury offers our first opportunity to explore this planet s miniature magnetosphere since the brief flybys of Mariner 10. Mercury s magnetosphere is unique in many respects. The magnetosphere of Mercury is among the smallest in the solar system; its magnetic field typically stands off the solar wind only - 1000 to 2000 km above the surface. For this reason there are no closed drift paths for energetic particles and, hence, no radiation belts. The characteristic time scales for wave propagation and convective transport are short and kinetic and fluid modes may be coupled. Magnetic reconnection at the dayside magnetopause may erode the subsolar magnetosphere allowing solar wind ions to impact directly the regolith. Inductive currents in Mercury s interior may act to modify the solar wind interaction by resisting changes due to solar wind pressure variations. Indeed, observations of these induction effects may be an important source of information on the state of Mercury s interior. In addition, Mercury s magnetosphere is the only one with its defining magnetic flux tubes rooted in a planetary regolith as opposed to an atmosphere with a conductive ionospheric layer. This lack of an ionosphere is probably the underlying reason for the brevity of the very intense, but short-lived, - 1-2 min, substorm-like energetic particle events observed by Mariner 10 during its first traversal of Mercury s magnetic tail. Because of Mercury s proximity to the sun, 0.3 - 0.5 AU, this magnetosphere experiences the most extreme driving forces in the solar system. All of these factors are expected to produce complicated interactions involving the exchange and re-cycling of neutrals and ions between the solar wind, magnetosphere, and regolith. The electrodynamics of Mercury s magnetosphere are expected to be equally complex, with strong forcing by the solar wind, magnetic reconnection at the magnetopause and in the tail, and the pick-up of planetary ions all driving field-aligned electric currents. However, these field-aligned currents do not close in an ionosphere, but in some other manner. In addition to the insights- into magnetospheric physics offered by study of the solar wind - Mercury system, quantitative specification of the "external" magnetic field generated by magnetospheric currents is necessary for accurate determination of the strength and multi-polar decomposition of Mercury s intrinsic magnetic field. MESSENGER S highly capable instrumentation and broad orbital coverage will greatly advance our understanding of both the origin of Mercury s magnetic field and the acceleration of charged particles in small magnetospheres. In. this article, we review what is known about Mercury s magnetosphere and describe the MESSENGER science team s strategy for obtaining answers to the outstanding science questions surrounding the interaction of the solar wind with Mercury and its small, but dynamic, magnetosphere.

Slavin, James A.↗

DOE OSTI2026

Concerted Electron-Ion Transport by Polyacrylonitrile Elucidated with Reactive Deep Learning Potentials

Charge transport in polymers, such as polyacrylonitrile (PAN), is crucial for electronics and energy storage. For instance, PAN can transport cations e.g., Li + , by facilitating dynamic cation-nitrile coordination in batteries. However, little is known regarding the underlying role of complex reactive polymer configurations. Herein, we develop a deep-learning potential, trained on ab initio energies and forces of nonequilibrium reactive PAN configurations, to unravel the kinetics of PAN cyclization initiated by a nucleophile (OH – dissociated from LiOH) attacking the terminal nitrile carbon. We find, based on the reaction free-energetics, rates, and charge analysis, that the nucleophile attack producing the first ring is the rate-limiting step, which subsequently triggers Li + -coupled electron transfer along the PAN backbone, causing ∼10 4 times faster sequential ring-formation of the remaining nitriles. PAN’s extended configurations, where dipolar and H-bonding interactions are minimal, enable such rapid kinetics. By validating our computational findings with IR and NMR experiments, we establish a pathway for designing reactive polymers with enhanced charge transport for energy applications.

Chahal-Crockett, Rajni [Oak Ridge National Laborat↗

DOE OSTI2021

Emission of electromagnetic waves as a stopping mechanism for nonlinear collisionless ionization waves in a high-β regime

A high energy density plasma embedded in a neutral gas is able to launch an outward-propagating nonlinear electrostatic ionization wave that traps energetic electrons. The trapping maintains a strong sheath electric field, enabling rapid and long-lasting wave propagation aided by field ionization. Here, using 1D3V kinetic simulations, we examine the propagation of the ionization wave in the presence of a transverse MG-level magnetic field with the objective to identify qualitative changes in a regime where the initial thermal pressure of the plasma exceeds the pressure of the magnetic field (Beta > 1). Our key finding is that the magnetic field stops the propagation by causing the energetic electrons sustaining the wave to lose their energy by emitting an electromagnetic wave. The emission is accompanied by the magnetic field expulsion from the plasma and an increased electron loss from the trapping wave structure. The described effect provides a mechanism mitigating rapid plasma expansion for those applications that involve an embedded plasma, such as high-flux neutron production from laser-irradiated deuterium gas jets.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Leaching Model of Radionuclides in Metal-Organic Framework Particles

Metal-organic frameworks (MOFs) have been used to sequester radionuclides and seal them inside of porous scaffolds using postsynthetic modification procedures. Experiments show that certain Zr-MOF with different capping linkers significantly affects the radionuclide release kinetics. In this work, we developed a leaching model of radionuclides in Zr-MOF particles. The model assumes that uranyl species occupy two energetically favored sites: the metal node and the MOF pores. For a given overall concentration of uranyl species, the partitions of uranyl species at the metal nodes and within the pores are determined by their chemical potentials. The model also considers the effect of particle surface and concentration on chemical potentials and diffusivity. The effects of spatial and structural dependent chemical potentials and diffusivity as well as particle sizes on leaching kinetics are investigated with the model. Predicted and measured uranyl leaching kinetics in Zr-MOF particles under batch experiments are compared. The results demonstrate the model’s capability for exploring the mechanisms of leaching and provide guidance for material design.

Metal-Organic Framework, Diffusion, uranyl, Leachi↗

DOE OSTI2023

Distinguishing photo-induced oxygen attack on alkyl chain versus conjugated backbone for alkylthienyl-benzodithiophene (BDTT)-based push–pull polymers

Synthetic design has enabled increasing power conversion efficiency advances in organic photovoltaics (OPV). One continuing knowledge gap is detailed understanding of single-material chemical (photo)stability. Many considerations are based on prior OPV-relevant donor homopolymer systems rather than the next-generation push–pull architectures. Generally, energetic offsets between the lowest occupied molecular orbital of the donor and molecular oxygen are assumed to dictate kinetics of photo-induced charge transfer to a super oxide radical. Herein we determine the ambient-induced photo-degradation pathways, presented as proposed site-specific arrow-pushing mechanisms, for five donor polymers all containing the same push unit – alkylthienyl-substituted-benzodithiophene (BDTT) – but with chemically distinct pull units. The donor-only polymer films were subject to controlled photobleaching in air as an accelerated degradation approach, coupled with simultaneous monitoring of absorptance. Sample subsets were periodically removed and analyzed with X-ray photoelectron spectroscopy, to evaluate near-surface chemical composition and oxygen additions to hetero reporter atoms on either the BDTT push unit or distinct pull unit (i.e., sulphur and nitrogen). Here, this methodology allows us to distinguish between different mechanisms of bond cleavage and formation. Overall, we find that neither polymer redox properties nor individual push or pull unit stability are sufficient to predict photo-oxidative degradation of these polymers. Rather, there is a greater dependence on the susceptibility of unique structural groups within a polymeric system. Alkyl chain oxygen addition is generally the first attack site and direct sulfur oxidation on the conjugated backbone occurs after saturation of the alkyl chain initiation sites. This work provides a standard that could be used to evaluate relative photo-oxidative (in)stability for new OPV materials quickly – prior to time-consuming device optimization – and demonstrates an effective methodology for correlating optical degradation with chemical structure alterations via spectroscopic signatures to guide synthetic design.

14 SOLAR ENERGY↗