Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “thermal plasma”

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 955 records · Page 53
All sourcesNASA NTRSDOE OSTI
NASA NTRS1986

The effect of beam-driven return current instability on solar hard X-ray bursts

The problem of electrostatic wave generation by a return current driven by a small area electron beam during solar hard X-ray bursts is discussed. The marginal stability method is used to solve numerically the electron and ion heating equations for a prescribed beam current evolution. When ion-acoustic waves are considered, the method appears satisfactory and, following an initial phase of Coulomb resistivity in which T sub e/T sub i rise, predicts a rapid heating of substantial plasma volumes by anomalous ohmic dissipation. This hot plasma emits so much thermal bremsstrahlung that, contrary to previous expectations, the unstable beam-plasma system actually emits more hard X-rays than does the beam in the purely collisional thick target regime relevant to larger injection areas. Inclusion of ion-cyclotron waves results in ion-acoustic wave onset at lower T sub e/T sub i and a marginal stability treatment yields unphysical results.

Cromwell, D.↗
NASA NTRS2023

Data-constrained Solar Modeling with GX Simulator

To facilitate the study of solarflares and active regions, we have created a modeling framework, the freelydistributed GX Simulator IDL package, that combines 3D magnetic and plasma structures with thermal andnonthermal models of the chromosphere, transition region, and corona. Its object-based modular architecture,which runs on Windows, Mac, and Unix/Linux platforms, offers the ability to either import 3D density andtemperature distribution models, or to assign numerically defined coronal or chromospheric temperatures anddensities, or their distributions, to each individual voxel. GX Simulator can apply parametric heating modelsinvolving average properties of the magneticfield lines crossing a given voxel, as well as compute and investigatethe spatial and spectral properties of radio,(sub)millimeter, EUV, and X-ray emissions calculated from the model,and quantitatively compare them with observations. The package includes a fully automatic model productionpipeline that, based on minimal users input, downloads the required SDO/HMI vector magneticfield data,performs potential or nonlinear force-freefield extrapolations, populates the magneticfield skeleton withparameterized heated plasma coronal models that assume either steady-state or impulsive plasma heating, andgenerates non-LTE density and temperature distribution models of the chromosphere that are constrained byphotospheric measurements. The standardized models produced by this pipeline may be further customizedthrough specialized IDL scripts, or a set of interactive tools provided by the graphical user interface. Here, wedescribe the GX Simulator framework and its applications.

Solar active regions↗
NASA NTRS1979

Two-dimensional potential double layers and discrete auroras

This paper is concerned with the formation of the acceleration region for electrons which produce the visible auroral arc and with the formation of the inverted V precipitation region. The former is embedded in the latter, and both are associated with field-aligned current sheets carried by plasma sheet electrons. It is shown that an electron current sheet driven from the plasma sheet into the ionosphere leads to the formation of a two-dimensional potential double layer. For a current sheet of a thickness less than the proton gyrodiameter solutions are obtained in which the field-aligned potential drop is distributed over a length much greater than the Debye length. For a current sheet of a thickness much greater than the proton gyrodiameter solutions are obtained in which the potential drop is confined to a distance on the order of the Debye length. The electric field in the two-dimensional double-layer model is the zeroth-order field inherent to the current sheet configuration, in contrast to those models in which the electric field is attributed to the first-order field due to current instabilities or turbulences. The maximum potential in the two-dimensional double-layer models is on the order of the thermal energy of plasma sheet protons, which ranges from 1 to 10 keV.

Kan, J. R.↗
NASA NTRS2022

Energy Partition at Collisionless Supercritical Quasiperpendicular Shocks

Collisionless shocks in astrophysical plasmas are important thermalizers, converting some of the incident flow energy into thermal energy, and non-thermalizers, partitioning that energy in unequal ways to different particle species, sub-populations thereof, and field components. This partition problem, or equivalently the shock equation of state, lies at the heart of shock physics. Here we employ systematically a framework to capture all the incident and downstream energy fluxes at two example traversals of the Earth’s bow shock by the Magnetospheric Multiscale Mission. Here and traditionally such data has to be augmented by information from other spacecraft, e.g., to provide more accurate measurements of the cold solar wind beam. With some care and fortuitous choices, the energy fluxes are constant, including instantaneous measurements through the shock layer. The dominant incident proton ram energy is converted primarily into downstream proton enthalpy flux, the majority of which is actually carried by a small fraction of suprathermal protons. Fluctuations include both real and instrumental effects. Separating these, resolving the solar wind beam, and other considerations point the way to a dedicated mission to solve this energy partition problem across a full range of plasma and shock conditions.

Steven J Schwartz↗
NASA NTRS1971

Active cleaning technique for removing contamination from optical surfaces in space

Plasma cleaning experiments were completed on hydrocarbon contaminant films, experiments were initiated to determine a satisfactory technique for depositing contaminant films, and an experiment was conducted to determine whether specimens are being thermally cleaned rather than plasma cleaned. Results of plasma cleaning experiments on hydrocarbon contaminant films showed that the optical properties of mirrors and gratings could be satisfactorily restored. Results on fused silica optical flats were inconclusive because of the insensitivity of measurement techniques to the contaminant films. White thermal control surfaces, degraded by the hydrocarbon contaminant film, could not be restored by oxygen plasma exposure. The reflectance of silvered FEP Teflon thermal control surfaces could be restored by plasma cleaning. Experiments with a silicone contaminant indicated that it could not be easily polymerized onto surfaces with ultraviolet radiation. Results of the thermal cleaning experiment showed that the polymerized hydrocarbon contaminant film could not be removed by heating in vacuum to a temperature in excess of that expected during plasma cleaning.

Gillette, R. B.↗
DOE OSTI2022

3D ion gyro-orbit heat load predictions for NSTX-U

Abstract High power tokamaks operate with divertor heat loads capable of destroying the plasma facing components (PFCs). High fidelity heat load predictions are necessary to ascertain the PFC state for design and during operation. Typical heat flux calculations are 2D, time invariant, and assume that power flows directly along the magnetic field lines (the optical approximation). These assumptions neglect the complex 3D geometries employed to protect the PFCs, the time varying nature of the plasma and PFC thermal state, and the helical trajectories of ions with finite Larmor radii (the gyro-orbit approximation). An integrated software framework, the heat flux engineering analysis toolkit (HEAT), was developed to generate time varying optical heat loads applied to real engineering computer aided design (CAD) (Looby et al 2022 Fusion Sci. Technol. 78 10–27). Recently, an ion-gyro orbit module has been added to HEAT. This module calculates the helical trajectories of ions as they gyrate about the magnetic field lines using kinetic theory macro-particles to accelerate the calculation. First, the new gyro-orbit module will be presented. Next, a comparison to existing research is performed. Finally, an analysis of the gyro-orbit heat loads for NSTX-U is presented for diverted discharges using the engineering CAD models utilized for PFC fabrication. Including these gyro-orbit effects can enhance the PFC performance by ‘smearing’ out the magnetic shadows associated with the castellated fish-scaled geometry. Simultaneously, the helical trajectories can degrade performance when they load narrow regions on edges and corners with high heat fluxes. Analysis of the trade-offs between these competing effects is included, and regions for further investigation are identified.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Pressure Response Optimization of an Eddy Current-Driven Flyer Plate Valve for the ITER Shattered Pellet Injection System

One technique for mitigating disruptions in a tokamak is shattered pellet injection (SPI). SPI is a process in which a large solid pellet consisting of deuterium, neon, or argon is desublimated in a pipe gun barsrel and launched downstream. Pellets are shattered just before entering the plasma by an impact with an angled tube. Injection of these materials into the plasma radiates stored thermal energy, limits current decay rates, suppresses the generation of runaway electrons, and dissipates runaway electrons if necessary. A critical element of the SPI system is a fast-acting valve that releases high-pressure gas to dislodge and accelerate pellets directly, or indirectly via a mechanical punch. A prototype valve sized for the ITER SPI system has been designed and fabricated. A pulsed high-voltage power supply energizes the valve’s internal magnetic coil, which induces eddy currents in the adjacent flyer plate resulting in a repulsive force between the flyer plate and the coil. The flyer plate action lifts a valve seat, allowing high-pressure gas to flow from the valve plenum to the downstream (breech) location of the pellet or mechanical punch. All of the valve’s internal components are designed to operate in ITER-level static background magnetic fields. Here, a study was conducted to optimize the downstream pressure response for a range of valve sizes and operating pressures. In particular, the study analyzes the breech pressure response associated with varying plenum pressures as well as varying breech volumes. A computational fluid dynamics simulation was built in STAR-CCM+ and validated against data from laboratory experiments. The resulting simulation outputs, in the form of downstream responses for a variety of initial plenum pressures and breech volumes, will be used as a complement to experimental data to ensure the pressure pulse is suitable for pellet survivability. These data, combined with studies on pellet shear strength and shock response, will be applied to optimization of overall operating parameters of the ITER SPI system.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Thermal stability of oxidized ultrathin Li films on TZM for plasma facing components

For improved understanding of the behavior of Li plasma facing components (PFCs) in the presence of oxygen impurities, we report on the thermal stability of ultrathin (up to 1.0 nm) LiO x films on a titanium-zirconium-molybdenum (TZM) alloy substrate over the temperature range of 310–1400 K. LiO x films were prepared by either post-oxidation of Li films or Li films deposited on an oxygen-precovered TZM surface. Thermal stability of these films was studied using a combination of temperature programmed desorption (TPD), X-ray photoelectron spectroscopy (XPS), and low-energy ion scattering (LEIS). For post-oxidized Li films, TPD showed that no Li desorption occurred until temperatures above 620 K, and then Li desorbed from the surface via at least three desorption peaks. The O 1s XPS spectra on the post-oxidized films at 310 K indicate the formation of lithium oxide (Li2O) and peroxide/hydroxide (Li 2 O 2 /LiOH). The peroxide/hydroxide converts to oxide after heating to 680 K with no desorption of Li or O 2 , and then this film decomposes to liberate Li into the gas phase while leaving oxygen at the TZM surface. Heating the LiO x films to 1070 K led to a solid-state reaction with surface Mo atoms of the substrate to form a condensed binary lithium molybdenum oxide (LixMoOy) phase, which upon further heating decomposed to liberate gas phase binary lithium molybdenum oxide species. Li deposition on an oxygen-precovered TZM surface formed a complex, highly stable, oxygen rich, Li-O-Mo interfacial oxide. Li desorption from multilayer Li films on oxygen-precovered TZM surfaces occurred in a metallic Li multilayer peak and three other oxide-derived peaks. These results support the improved thermal stability of oxidized lithium as compared to metallic Li and have implications for Mo impurity transport into the plasma at high temperatures.

36 MATERIALS SCIENCE↗
DOE OSTI2025

On the role of ion temperature gradient turbulence in driving ion thermal transport in neutral beam injection-heated L-mode plasmas in a superconducting tokamak

In this paper, we report a detailed experimental study of the role of ion-temperature-gradient driven turbulence in driving ion thermal transport in neutral beam injection (NBI)-heated L-mode plasmas in Experimental Advanced Superconducting Tokamak (EAST) (Wan et al 2000 Nucl. Fusion 40 1057) coupled with linear and nonlinear gyrokinetic simulations. Significant ion-scale turbulence k ⟂ < cm −1 (k ⟂ < 1.5, where k ⟂ where is the perpendicular wavenumber and ρ s is the ion gyroradius calculated using local electron temperature T e ), measured with a microwave reflectometer, is observed in the plasma core at for r/a ≈ 0.265 and 0.5 (where r is half the diameter of the closed flux surface at a given radial position, and a is half the diameter of the last closed flux surface). Local linear stability analysis with experimental equilibrium quantities at these two radial locations using the GS2 gyrokinetic code shows that the most unstable ion-scale micro-instability is the ion temperature gradient (ITG) mode. Since the computed maximum ion-scale linear growth rates are larger than the local Waltz-Miller E × B shearing rate (Waltz and Miller 1999 Phys. Plasmas 6 4265), the E × B shear is unable to suppress ion-scale turbulence, consistent with the experimental observation of ion-scale turbulence. Experimental ion and electron thermal transports, calculated with power balance analysis, are both anomalous with ion thermal transport being the dominant channel. Nonlinear gyrokinetic simulations using the electrostatic global particle-in-cell gyrokinetic tokamak simulation (GTS) code (Weixing Wang et al 2010 Phys. Plasmas 17 072511) show good quantitative agreement between predicted ion thermal transport and that from the power balance analysis. These results demonstrate for the first time the important role of ITG turbulence in driving thermal transport in NBI-heated L-mode plasmas in EAST, showing the effectiveness of GTS code in simulating electrostatic turbulence in EAST and the need of suppressing ITG turbulence to improve energy confinement in EAST.

Sun, P. J. [Chinese Academy of Sciences (CAS), Hef↗
DOE OSTI2023

K-shell radiation and bright spot characteristics of high-energy-density Fe-Cr-Ni plasmas influenced by X-pinch load geometry

K-shell radiation is diagnostically advantageous for the study of high-energy-density plasmas. Here, the current work investigates the influence of load geometry on K-shell x-rays in Z-pinch plasmas produced on the 1-MA Zebra generator. Stainless steel (Fe: 69%, Cr: 20%, Ni: 9%) X-pinch wire loads are fielded with an interwire angle of 31° (small-angle) or 62.5° (large-angle), respectively and studied using various x-ray diode, imaging, and spectroscopic diagnostics. The large-angle geometry produced large, individual x-ray-emitting bright spots (> 3 keV) at the wire cross-point averaging areas ≥ 1 mm 2 , with noteworthy soft x-ray (> 0.75 keV) bursts and radiation yields ≤ 14.6 kJ. The small-angle geometry produced multiple x-ray-emitting bright spots extending the pinch axis with smaller average sizes (≤ 0.5 mm 2 ), consistent with stronger current, increased radiation yield (≥ 15.5 kJ), and notable hard x-ray (> 9 keV) bursts, which peak with largest cumulative bright spot size. Spectroscopic analysis is performed with non-LTE collisional-radiative modeling, indicating a cold, nonthermal plasma region radiating from neutral to Ne-like ions, a hot, thermal region radiating from Fe and Cr He- and Li-like ions, and an intermediate region radiating Kα satellites from Li-like to O-like ions. Modeling of the nonthermal and satellite features imply a fractional hot electron abundance ~0.1% and ~0.5% for large- and small-angle geometries, respectively. Relative intensity analysis is performed on Heα and Kα lines, revealing intensity ratios that deviate from the original wire composition. Optical depth estimates allude to optically thick, thermal K-shell plasmas with densities consistent with non-LTE modeling.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2021

Correlation between ferroelectricity and ferroelectric orthorhombic phase of Hf x Zr 1-x O 2 thin films using synchrotron x-ray analysis

The change in the interplanar spacing ( d -spacing) including the ferroelectric orthorhombic (O) phase in the low-temperature fabricated Hf x Zr 1-x O 2 (HZO) films was studied using synchrotron grazing-incidence wide-angle x-ray scattering analysis. The 10-nm-thick HZO films were fabricated by thermal and plasma-enhanced atomic layer deposition (TH- and PE-ALD) methods using H 2 O gas and O 2 plasma as oxidants, respectively, and a post-metallization annealing (PMA) was performed at 300–400 °C. The d -spacing of the mixture of (111)-, (101)-, and (111)-planes of O, tetragonal (T), and cubic (C) phases, respectively, for the TH- and PE-ALD HZO films increased up to 2.99 Å with an increase in PMA temperature, while the d -spacing estimated by conventional x-ray diffraction was 2.92 Å regardless of the PMA temperature. The remanent polarization (2 P r = P r + - P r - ) of the HZO films increased as the PMA temperature increased. It is clear that the 2 P r value satisfied a linear relationship as a function of the d -spacing of O(111)/T(101)/C(111) phases. Furthermore, the wake-up effect was found to depend on the ferroelectric O phase formation. The wake-up effect was significantly reduced in both the TH- and PE-ALD HZO films after the PMA at 400 °C due to the increase in the ferroelectric O phase formation. The leakage current density ( J )–electric field properties of the PE-ALD HZO film with the lowest d -spacing were divided into three steps, such as low, middle, and large J values, in the wake-up (10 3 cycles), pristine (10 0 cycle), and fatigue (10 7 cycles) states, respectively. Therefore, an analysis of the ferroelectric O phase is very important for understanding the ferroelectricity including endurance.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Diamond surface functionalization via visible light–driven C–H activation for nanoscale quantum sensing

Nitrogen-vacancy (NV) centers in diamond are a promising platform for nanoscale NMR sensing. Despite significant progress toward using NV centers to detect and localize nuclear spins down to the single spin level, NV-based spectroscopy of individual, intact, arbitrary target molecules remains elusive. Such sensing requires that target molecules are immobilized within nanometers of NV centers with long spin coherence. The inert nature of diamond typically requires harsh functionalization techniques such as thermal annealing or plasma processing, limiting the scope of functional groups that can be attached to the surface. Solution-phase chemical methods can be readily generalized to install diverse functional groups, but they have not been widely explored for single-crystal diamond surfaces. Moreover, realizing shallow NV centers with long spin coherence times requires highly ordered single-crystal surfaces, and solution-phase functionalization has not yet been shown with such demanding conditions. In this work, we report a versatile strategy to directly functionalize C–H bonds on single-crystal diamond surfaces under ambient conditions using visible light, forming C–F, C–Cl, C–S, and C–N bonds at the surface. This method is compatible with NV centers within 10 nm of the surface with spin coherence times comparable to the state of the art. As a proof-of-principle demonstration, we use shallow ensembles of NV centers to detect nuclear spins from surface-bound functional groups. Our approach to surface functionalization opens the door to deploying NV centers as a tool for chemical sensing and single-molecule spectroscopy.

Science & Technology - Other Topics↗
DOE OSTI2021

Chirping ion cyclotron emission (ICE) on NSTX-U

In this work, we report on the discovery on NSTX-U of a qualitatively new type of non-thermal emission from plasmas in the ion-cyclotron frequency range (ICE). This new type of ICE is weakly damped, γ damp /ω ICE ≈ 0.14%, low order eigenmodes, excited by a fast ion population radially located near an internal transport barrier. The ICE appears as repetitive bursts with durations of ≈1 ms, and the bursts exhibit frequency chirps upwards and/or downwards with maximum δf/f ≈ 1%. We refer to this chirping ICE as ch-ICE. The bursts are longer than the typical ICE bursts of ≈100 μs reported previously for NSTX and NSTX-U, for simplicity, NSTX(-U). The chirping resembles that seen in modeling of weakly unstable fast particle driven instabilities (cf, Berk et al 1997 Phys. Lett. A 234 213). These data also provide some preliminary indications of non-linear coupling between the ch-ICE harmonics. Despite extensive experimental and theoretical studies of ICE, theoretical understanding of ICE remains qualitative at best. The detailed experimental observations of ICE from NSTX(-U) find a wide variation in the ICE characteristics. The new observations of ch-ICE presented here, combined with the characteristics of ICE on NSTX(-U) previously reported (Fredrickson et al 2019 Phys. Plasmas 26 032111), could guide the development of more complete theoretical models of ICE.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Dependence of ion-cyclotron range of frequencies instabilities on species mix and fast-ion distribution: II. Time evolution

Dedicated DIII-D experiments studied instabilities below the ion cyclotron frequency f ci and at cyclotron harmonics by injecting thirteen different neutral beam populations into plasmas with different thermal compositions of H and D at five different values of toroidal field. The resulting database of 2529 conditions is described in Part 1: Stability. This paper focuses on the temporal evolution. For some beam geometries, the ion cyclotron emission (ICE) is much stronger at the beginning of injection, suggesting the importance of ‘bump-on-tail’ energy drive. Sawtooth crashes cause the emission to cease for some beam types and emission bands, but the opposite effect occurs for others. When unstable, for virtually all conditions and frequency bands, the emission bursts in a cycle reminiscent of predator-prey oscillations. Evidence for nonlinear interaction between ICE harmonics is observed. Strongly unstable modes typically appear within 1 ms from the start of beam injection and often decay within 0.1 ms when the beam turns off, but both the initial rise and final decay of the mode power exhibit a wide variety of behavior.

electromagnetic ion cyclotron instabilities↗
DOE OSTI2021

Search for GeV neutrino emission during intense gamma-ray solar flares with the IceCube Neutrino Observatory

Solar flares convert magnetic energy into thermal and nonthermal plasma energy, the latter implying particle acceleration of charged particles such as protons. Additionally, protons are injected out of the coronal acceleration region and can interact with dense plasma in the lower solar atmosphere, producing mesons that subsequently decay into gamma rays and neutrinos at O(MeV-GeV) energies. We present the results of the first search for GeV neutrinos emitted during solar flares carried out with the IceCube Neutrino Observatory. While the experiment was originally designed to detect neutrinos with energies between 10 GeV and a few PeV, a new approach allowing for a O(GeV) energy threshold will be presented. The resulting limits allow us to constrain some of the theoretical estimates of the expected neutrino flux.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Simulating of magnetic reconnection in solar flares

Magnetic reconnection is an astrophysical process where neighboring magnetic field lines running anti-parallel are reconfigured. Observations of solar flares are what pushed for serious research regarding magnetic reconnection as it seemed to be the underlying mechanism. The reconfiguration of the field lines results in an explosive release of energy as magnetic field energy is converted into plasma kinetic and thermal energies. Our work begins with Athena++, an astrophysical magnetohydrodynamic (MHD) simulation code, and an existing magnetic reconnection setup.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Theoretical studies of chemical reactions related to the formation and growth of polycyclic aromatic hydrocarbons (PAH) and molecular properties of their key intermediates (Final Progress Report)

The formation mechanisms of polycyclic aromatic hydrocarbons, (PAHs) – organic molecules carrying fused benzene rings – are of great interest to scientists and engineers due to their importance in combustion chemistry and astrochemistry. On Earth, PAHs are largely produced in incomplete combustion of fossil fuel and are considered as critical precursors to unwanted soot particles leading to combustion inefficiency and causing air pollution along with detrimental health effects. Simple PAH molecules initially formed in the gas phase, are further involved in a build-up process in combustion flames leading to larger PAH, bowl-shaped nanostructures, fullerenes, and solid-phase species including carbonaceous dust, graphene particles, and soot. In deep space, PAH and their derivatives are potential key intermediates and nucleation sites leading eventually to carbonaceous nanoparticles (“interstellar grains”). Therefore, the understanding of the key processes in the synthesis of PAHs along with their precursors and their degradation mechanisms in combustion systems and in interstellar, circumstellar, and planetary atmospheric environments will provide critical insights into how complex aromatic structures, carbonaceous nanoparticles, and fullerenes are formed and destroyed. Achieving this understanding is an important step in the development of the efficient combustion processes and of the ecofriendly devices with reduced environmental pollution as well as technological strategies for the production of hydrogen and solid carbon through thermal or plasma-assisted pyrolysis of natural gas and biomass. Also, the understanding of the key processes of PAH and soot growth will help in our comprehension of chemical evolution in the universe. Detailed information on the mechanisms and reliable rate constants of the key elementary chemical reactions involved in PAH formation and destruction processes and in inception of soot particles is often missing, with the main deficiencies being the absence of temperature- and pressure-dependent rate constants for the broad range of conditions occurring in various terrestrial and interstellar processes and the lack of data on the reaction products and their branching ratios. Complementary to experimental studies, these gaps in knowledge can be filled by using quantum chemical calculations of reaction potential energy surfaces providing us with accurate energies of reaction products, intermediates, and transition states, revealing the reaction mechanism, and giving the molecular properties required to compute rate constants for relevant reaction steps and product branching ratios using the RRKM-Master Equation (ME) method. Molecular dynamics (MD) simulations can be used in cases when a reaction rate cannot be properly described by statistical theories. During the terminal renewal project period we employed these ab initio/RRKM-ME and MD approaches to complete our studies on several key reactions relevant to the formation/growth of PAH and inception of soot particles including (1) the reaction mechanism and kinetics of the resonance stabilized fulvenallenyl radical with propargyl and C 3 H 4 isomers; (2) the reaction mechanism and kinetics for the C + indene and C 2 + styrene reactions producing naphthyl or azulenyl radicals in low-temperature environments; (3) the MD study of non-equilibrium dimerization of acepyrene and coronene and its radical. The information derived from our theoretical calculations contributed to a better fundamental understanding of the reaction mechanisms and provide missing critical kinetic data to improve combustion models of hydrocarbon fuels and astrochemical models of the growth of carbonaceous molecules and particles in cold molecular clouds, circumstellar envelopes, and planetary atmospheres.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗