Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ionization”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
NASA NTRS1994

Mass-loss rates, ionization fractions, shock velocities, and magnetic fields of stellar jets

In this paper we calculate emission-line ratios from a series of planar radiative shock models that cover a wide range of shock velocities, preshock densities, and magnetic fields. The models cover the initial conditions relevant to stellar jets, and we show how to estimate the ionization fractions and shock velocities in jets directly from observations of the strong emission lines in these flows. The ionization fractions in the HH 34, HH 47, and HH 111 jets are approximately 2%, considerably smaller than previous estimates, and the shock velocities are approximately 30 km/s. For each jet the ionization fractions were found from five different line ratios, and the estimates agree to within a factor of approximately 2. The scatter in the estimates of the shock velocities is also small (+/- 4 km/s). The low ionization fractions of stellar jets imply that the observed electron densities are much lower than the total densities, so the mass-loss rates in these flows are correspondingly higher (approximately greater than 2 x 10(exp -7) solar mass/yr). The mass-loss rates in jets are a significant fraction (1%-10%) of the disk accretion rates onto young stellar objects that drive the outflows. The momentum and energy supplied by the visible portion of a typical stellar jet are sufficient to drive a weak molecular outflow. Magnetic fields in stellar jets are difficult to measure because the line ratios from a radiative shock with a magnetic field resemble those of a lower velocity shock without a field. The observed line fluxes can in principle indicate the strength of the field if the geometry of the shocks in the jet is well known.

Hartigan, Patrick↗
NASA NTRS1994

Theoretical investigation of wave-vector-dependent analytical and numerical formulations of the interband impact-ionization transition rate for electrons in bulk silicon and GaAs

The electron interband impact-ionization rate for both silicon and gallium arsenide is calculated using an ensemble Monte Carlo simulation with the expressed purpose of comparing different formulations of the interband ionization transition rate. Specifically, three different treatments of the transition rate are examined: the traditional Keldysh formula, a new k-dependent analytical formulation first derived by W. Quade, E. Scholl, and M. Rudan (1993), and a more exact, numerical method of Y. Wang and K. F. Brennan (1994). Although the completely numerical formulation contains no adjustable parameters and as such provides a very reliable result, it is highly computationally intensive. Alternatively, the Keldysh formular, although inherently simple and computationally efficient, fails to include the k dependence as well as the details of the energy band structure. The k-dependent analytical formulation of Quade and co-workers overcomes the limitations of both of these models but at the expense of some new parameterization. It is found that the k-dependent analytical method of Quade and co-workers produces very similar results to those obtained with the completely numerical model for some quantities. Specifically, both models predict that the effective threshold for impact ionization in GaAs and silicon is quite soft, that the majority of ionization events originate from the second conduction band in both materials, and that the transition rate is k dependent. Therefore, it is concluded that the k-dependent analytical model can qualitatively reproduce results similar to those obtained with the numerical model yet with far greater computational efficiency. Nevertheless, there exist some important drawbacks to the k-dependent analytical model of Quade and co-workers: These are that it does not accurately reproduce the quantum yield data for bulk silicon, it requires determination of a new parameter, related physically to the overlap intergrals of the Bloch state which can only be adjusted by comparison to experiment, and fails to account for any wave-vector dependence of the overlap integrals. As such the transition rate may be overestimated at those points for which 'near vertical,' small change in k, transitions occur.

Kolnik, Jan↗
NASA NTRS1995

High-resolution measurements of the K-alpha spectra of low-ionizationm species of iron: A new spectral signature of nonequilibrium ionization conditions in young supernova remnants

We present the first systematic laboratory measurements of high-resolution K-alpha spectra of intermediate ions of iron, Fe X-XVII. These lines are not produced in collisional equilibrium plasmas because of the relevant charge states cannot exist at the high electron temperatures required for appreciable excitation of the K-alpha transitions. However, they can provide excellent spectral diagnostics for nonequilibrium ionization conditions, such the ionizing plasmas of young supernova remnants. To facilitate the line identifications, we compare our spectra with theoretical atomic calculations performed using multiconfiguration parametric potential and Dirac-Fock atomic codes. Our measurements also allow direct comparison with time-dependent ionization balance calculations for ionizing plasmas, and good agreement is found.

Decaux, V.↗
NASA NTRS1994

Numerical study of the wave-vector dependence of the electron interband impact ionization rate in bulk GaAs

Ensemble Monte Carlo calculations of the electron interband impact ionization rate in bulk GaAs are presented using a wave-vector (k)-dependent formulation of the ionization transition rate. The transition rate is evaluated through the use of numerically generated wavefunctions determined via a k-p calculation within the first two conduction bands at numerous points within a finely spaced three-dimensional grid in k space. The transition rate is determined to be greatest for states within the second conduction band. Is is found that the interband impact ionization transition rate in bulk GaAs is best characterized as having an exceedingly soft threshold energy. As a consequence, the dead space, defined as the distance over which the ionization probability for a given carrier is assumed to be zero, is estimated to be much larger than that estimated using a harder threshold. These results have importance in the design of the multiquantum-well avalanche photodiodes.

Wang, Yang↗
NASA NTRS1994

A discharge flow-photoionization mass spectrometric study of the FO(X 2 Pi i) radical. Photoionization efficiency spectrum and ionization energy

Photoionization efficiency spectra of FO were measured over the wavelength range 80.0-100.0 nm and in the ionization threshold region, 94.0-100.0 nm, using a discharge flow-photoionization mass spectrometer apparatus coupled to a synchrotron radiation source. FO was generated by the reaction of F2P atoms with NO3 and via a F2O2 discharge. A value of 12.78 +/- 0.03 eV was obtained for the adiabatic ionization energy of FO from photoion thresholds which corresponds to FO(+)(X 3 Sigma -) from FO(X 2 Pi i). These results, which are the first to be obtained by direct Photo-ionization mass spectrometry (PIMS) measurements, corroborate those of a photoelectron spectroscopy (PES) study; however, the ionization energy determined here is free from interferences due to other species which complicated the PES measurement. A value of 109.5 +/- 8.0 kJ/mol for Delta f H 0 298(FO) is computed from the present value of IE(FO) and a previous appearance energy measurement, and a value for the proton affinity of FO is calculated to be 511.5 +/- 10.0 kJ/mol.

Zhang, Zhengyu↗
NASA NTRS1995

Ensemble Monte Carlo calculation of the hole initiated impact ionization rate in bulk GaAs and silicon using a k-dependent, numerical transition rate formulation

The hole initiated impact ionization rate in bulk silicon and GaAs is calculated using a numerical formulation of the impact ionization transition rate incorporated into an ensemble Monte Carlo simulation. The transition rate is calculated from Fermi's golden rule using a two-body screened Coulomb interaction including a wavevector dependent dielectric function. It is found that the effective threshold for hole initiated ionization is relatively soft in both materials, that the split-off band dominates the ionization process in GaAs. and that no clear dominance by any one band is observed in silicon, though the rate out of the light hole band is greatest.

Oguzman, Ismail H.↗
NASA NTRS1996

Cosmic Metal Production and the Contribution of QSO Absorption Systems to the Ionizing Background

The recent discovery by Cowie et al. (1995) and Tytler et al. (1995) of metals in the Ly alpha clouds shows that the intergalactic medium (IGM) at high redshift is contaminated by the products of stars and suggests that ionizing photons from massive star formation may be a significant contributor to the UV background radiation at early epochs. We assess the validity of the stellar photoionization hypothesis. Based on recent computations of metal yields and 0-star Lyman continuum (Lyc) fluxes, we find that 0.2 percent of the rest-mass energy of the metals produced is radiated as Lyc. By modeling the transfer of ionizing radiation through the IGM and the rate of chemical enrichment, we demonstrate that the background intensity of photons at 1 ryd that accompanies the production of metals in the Ly alpha forest clouds may be significant, approaching 0.5 x 10(exp -21) ergs cm squared s(-1) Hz(-1) sr(-1) at z approximately equals 3 if the Lyc escape fraction is greater than of equal to 0.25. Together with quasars, massive stars could then, in principle, provide the hydrogen and helium Lyc photons required to ionize the universe at high redshifts. We propose that observations of the He2 Gunn-Peterson effect and of the metal ionization states of the Ly alpha forest and Lyman-limit absorbers should show the signature of a stellar spectrum. We also note that the stellar photoionization model fails if a large fraction of the UV radiation emitted from stars cannot escape into the IGM, as suggested by the recent Hopkins Ultraviolet Telescope observations by Leitherer et al. (1995) of low-redshift starburst galaxies, or if most of the metals observed at z is approximately 3 were produced at much earlier epochs.

Madau, Piero↗
NASA NTRS1996

Spectra of High-Ionization Seyfert 1 Galaxies: Implications for the Narrow-Line Region

We present line profiles and profile parameters for the Narrow-Line Regions (NLRs) of six Seyfert 1 galaxies with high-ionization lines: MCG 8-11-11, Mrk 79, Mrk 704, Mrk 841, NGC 4151, and NGC 5548. The sample was chosen primarily with the goal of obtaining high-quality [Fe VII] lambda6087 and, when possible, [Fe X] lambda6374 profiles to determine if these lines are more likely formed in a physically distinct 'coronal line region' or are formed throughout the NLR along with lines of lower critical density (n(sub cr)) and/or Ionization Potential (IP). We discuss correlations of velocity shift and width with n(sub cr) and IP. In some objects, lines of high IP and/or n(sub cr) are systematically broader than those of low IP/n(sub cr). Of particular interest, however, are objects that show no correlations of line width with either IP or n(sub cr). In these objects, lines of high and low IP/n(sub cr), are remarkably similar, which is difficult to reconcile with the classical picture of the NLR, in which lines of high and low IP/n(sub cr) are formed in physically distinct regions. We argue for similar spatial extents for the flux in lines with similar profiles. Here, as well as in a modeling-oriented companion paper, we develop further an idea suggested by Moore & Cohen that objects that do and do not show line width correlations with IP/n(sub cr) can both be explained in terms of a single NLR model with only a small difference in the cloud column density distinguishing the two types of object. Overall, our objects do not show correlations between the Full Width at Half-Maximum (FWHM) and IP and/or n(sub cr). The width must be defined by a parameter that is sensitive to extended profile wings in order for the correlations to result. We present models in which FWHM correlations with IP and/or n(sub cr) result only after simulating the lower spectral resolution used in previous observational studies. The models that simulate the higher spectral resolution of our observational study produce line width correlations only if the width is defined by a parameter that is more sensitive to extended profile wings than is the FWHM. Our sample of six objects is in effect augmented by incorporating the larger sample (16 objects) of Veilleux into some of our discussion. This paper focuses on new interpretations of NLR emission-line spectra and line profiles that stem directly from the observations. Paper 2 focuses on modeling and complements this paper by illustrating explicitly the effects that spatial variations in electron density, ionization parameter, and column density have on model profiles. By comparing model profiles with the observed profiles presented here, as well as with those presented by Veilleux, Paper 2 yields insight into how the electron density, ionization parameter, and column density likely vary throughout the NLR.

Moore, David↗
NASA NTRS1996

Three-Dimensional Electromagnetic Monte Carlo Particle-in-Cell Simulations of Critical Ionization Velocity Experiments in Space

Although the existence of the critical ionization velocity (CIV) is known from laboratory experiments, no agreement has been reached as to whether CIV exists in the natural space environment. In this paper we move towards more realistic models of CIV and present the first fully three-dimensional, electromagnetic particle-in-cell Monte-Carlo collision (PIC-MCC) simulations of typical space-based CIV experiments. In our model, the released neutral gas is taken to be a spherical cloud traveling across a magnetized ambient plasma. Simulations are performed for neutral clouds with various sizes and densities. The effects of the cloud parameters on ionization yield, wave energy growth, electron heating, momentum coupling, and the three-dimensional structure of the newly ionized plasma are discussed. The simulations suggest that the quantitative characteristics of momentum transfers among the ion beam, neutral cloud, and plasma waves is the key indicator of whether CIV can occur in space. The missing factors in space-based CIV experiments may be the conditions necessary for a continuous enhancement of the beam ion momentum. For a typical shaped charge release experiment, favorable CIV conditions may exist only in a very narrow, intermediate spatial region some distance from the release point due to the effects of the cloud density and size. When CIV does occur, the newly ionized plasma from the cloud forms a very complex structure due to the combined forces from the geomagnetic field, the motion induced emf, and the polarization. Hence the detection of CIV also critically depends on the sensor location.

Wang, J.↗
NASA NTRS1994

Semiclassical Study of the Wave Vector Dependence of the Interband Impact Ionization Rate in Bulk Silicon

We present calculations of the interband impact ionization rate calculated using a wave vector dependent (k-dependent) semiclassical formulation of the transition rate. The transition rate is determined using Fermi's golden rule from a two-body screened Coulomb interaction assuming energy and momentum conservation. The transition rate is calculated for the first two conduction bands of silicon by numerically integrating over the full Brillouin zone. The overlap integrals in the expression for the transition rate are determined numerically using a 15 band k-p calculation. It is found that the transition rate depends strongly on the initiating electron wave vector (k vector) and that the transition rate is greatest for electrons originating within the second conduction band than the first conduction band. An ensemble Monte Carlo simulation, which includes the numerically determined ionization transition rate as well as the full details of the first two conduction bands, is used to calculate the total impact ionization rate in bulk silicon. Good agreement with the experimentally determined electron ionization rate data is obtained.

Wang, Yang↗
NASA NTRS1994

Numerical Study of the Wave-Vector Dependence of the Electron Interband Impact Ionization Rate in Bulk GaAs

Ensemble Monte Carlo calculations of the electron interband impact ionization rate in bulk GaAs are presented using a wave-vector (k)-dependent formulation of the ionization transition rate. The transition rate is evaluated through use of numerically generated wavefunctions determined via a k-p calculation within the first two conduction bonds at numerous points within a finely spaced three-dimensional grid in k space. The transition rate is determined to be greatest for states within the second conduction band. It is found that the interband impact ionization transition rate in bulk GaAs is best characterized as having an exceedingly "soft" threshold energy. As a consequence, the dead space, defined as the distance over which the ionization probability for a given carrier is assumed to be zero, is estimated to be-much larger than that estimated using a "harder" threshold. These results have importance in the design of multiquantum-well avalanche photodiodes.

Wang, Yang↗
NASA NTRS1994

Theoretical Investigation of Wave-Vector-Dependent Analytical and Numerical Formulations of the Interband Impact-Ionization Transition Rate for Electron in Bulk Silicon and GaAs

The electron interband impact-ionization rate for both silicon and gallium arsenide is calculated using an ensemble Monte Carlo simulation with the expressed purpose of comparing different formulations of the interband ionization transition rate. Specifically, three different treatments of the transition rate are examined: the traditional Keldysh formula, a new k-dependent analytical formulation first derived by W. Quade, E Scholl, and M. Rudan, and a more exact, numerical method of Y. Wang and K. F. Brennan. Although the completely numerical formulation contains no adjustable parameters and as such provides a very reliable result, it is highly computationally intensive. Alternatively, the Keldysh formula, although inherently simple and computationally efficient, fails to include the k dependence as well as the details of the energy band structure. The k-dependent analytical formulation of Quade and co-workers overcomes the limitations of both of these models but at the expense of some new parameterization. It is found that the k-dependent analytical method of Quade and co-workers produces very similar results to those obtained with (he completely numerical model for some quantities. Specifically, both models predict that the effective threshold for impact ionization in GaAs and silicon is quite soft, that the majority of ionization events originate from the second conduction band in both materials, and that the transition rate is k dependent. Therefore, it is concluded that the k-dependent analytical model can qualitatively reproduce results similar to those obtained with the numerical model yet with far greater computational efficiency. Nevertheless, there exist some important drawbacks to the k-dependent analytical model of Quade and co-workers: These are that it does not accurately reproduce the quantum yield data for bulk silicon, it requires determination of a new parameter, related physically to (he overlap integrals of the Bloch state which can only be adjusted by comparison to experiment, and fails to account for any wave-vector dependence of the overlap integrals. As such [he transition rate may be overestimated at those points for which "near vertical," small change in k, transitions occur.

Kolnik, Jan↗
NASA NTRS2001

Viscosity Coefficient Curve Fits for Ionized Gas Species Grant Palmer

Viscosity coefficient curve fits for neutral gas species are available from many sources. Many do a good job of reproducing experimental and computational chemistry data. The curve fits are usually expressed as a function of temperature only. This is consistent with the governing equations used to derive an expression for the neutral species viscosity coefficient. Ionized species pose a more complicated problem. They are subject to electrostatic as well as intermolecular forces. The electrostatic forces are affected by a shielding phenomenon where electrons shield the electrostatic forces of positively charged ions beyond a certain distance. The viscosity coefficient for an ionized gas species is a function of both temperature and local electron number density. Currently available curve fits for ionized gas species, such as those presented by Gupta/Yos, are a function of temperature only. What they did was to assume an electron number density. The problem is that the electron number density they assumed was unrealistically high. The purpose of this paper is two-fold. First, the proper expression for determining the viscosity coefficient of an ionized species as a function of both temperature and electron number density will be presented. Then curve fit coefficients will be developed using the more realistic assumption of an equilibrium electron number density. The results will be compared against previous curve fits and against highly accurate computational chemistry data.

Palmer, Grant↗
NASA NTRS2002

Hydrodynamic Models of Line-Driven Accretion Disk Winds III: Local Ionization Equilibrium

We present time-dependent numerical hydrodynamic models of line-driven accretion disk winds in cataclysmic variable systems and calculate wind mass-loss rates and terminal velocities. The models are 2.5-dimensional, include an energy balance condition with radiative heating and cooling processes, and includes local ionization equilibrium introducing time dependence and spatial dependence on the line radiation force parameters. The radiation field is assumed to originate in an optically thick accretion disk. Wind ion populations are calculated under the assumption that local ionization equilibrium is determined by photoionization and radiative recombination, similar to a photoionized nebula. We find a steady wind flowing from the accretion disk. Radiative heating tends to maintain the temperature in the higher density wind regions near the disk surface, rather than cooling adiabatically. For a disk luminosity L (sub disk) = solar luminosity, white dwarf mass M(sub wd) = 0.6 solar mass, and white dwarf radii R(sub wd) = 0.01 solar radius, we obtain a wind mass-loss rate of M(sub wind) = 4 x 10(exp -12) solar mass yr(exp -1) and a terminal velocity of approximately 3000 km per second. These results confirm the general velocity and density structures found in our earlier constant ionization equilibrium adiabatic CV wind models. Further we establish here 2.5D numerical models that can be extended to QSO/AGN winds where the local ionization equilibrium will play a crucial role in the overall dynamics.

Pereyra, Nicolas Antonio↗
NASA NTRS2002

Dissociative Ionization of Benzene by Electron Impact

We report a theoretical study of the dissociative ionization (DI) of benzene from the low-lying ionization channels. Our approach makes use of the fact that electron motion is much faster than nuclear motion and DI is treated as a two-step process. The first step is electron-impact ionization resulting in an ion with the same nuclear geometry as the neutral molecule. In the second step the nuclei relax from the initial geometry and undergo unimolecular dissociation. For the ionization process we use the improved binary-encounter dipole (iBED) model. For the unimolecular dissociation step, we study the steepest descent reaction path to the minimum of the ion potential energy surface. The path is used to analyze the probability of unimolecular dissociation and to determine the product distributions. Our analysis of the dissociation products and the thresholds of the productions are compared with the result dissociative photoionization measurements of Feng et al. The partial oscillator strengths from Feng et al. are then used in the iBED cross section calculations.

Huo, Winifred↗
NASA NTRS2003

The Ionized Gas and Nuclear Environment in NGC 3783: Variability and Modeling of the 900 ks CHANDRA Spectrum - IV

We present a detailed analysis of the 900 ks spectrum of NGC3783 obtained by Chandra in 2000-2001 (Kaspi et al. 2002). We split the data in various ways to look for time dependent and luminosity dependent spectral variations. This analysis, the measured equivalent widths of a large number of X-ray lines, and our photoionization calculations, lead us to the following conclusions: 1) NGC 3783 fluctuated in luminosity, by a factor N 1.5, during individual 170 ks observations. The fluctuations were not associated with significant spectral variations. 2) On a longer time scale, of 20-120 days, we discovered two very different spectral shapes that are noted the high state and the low state spectra. The observed changes between the two can be described as the appearance and disappearance of a soft continuum component. The spectral variations are not related, in a simple way, to the brightening or the fading of the short wavelength continuum, as observed in other objects. NGC3783 seems to be the first AGN to show this unusual behavior. 3) The appearance of the soft continuum component is consistent with beeing the only spectral variation and there is no need to invoke changes in the absorber s opacity. In particular, all absorption lines with reliable measurements show the same equivalent width, within the observational uncertainties, during high and low states. 4) Photoionization model calculations show that a combination of three ionization components, each split into two kinematic components, explain very well the intensity of almost all absorption lines and the bound-free absorption. The components span a large range of ionization and a total column of about 3 x 10(exp 22) per square centimeter Moreover, all components are thermally stable and are situated on the vertical branch of the stability curve.. This means that they are in pressure equilibrium and perhaps occupy the same volume of space. This is the first detection of such a multi-component equilibrium gas in AGN. 5) The only real discrepancy between the model and the observations is the wavelength of the iron M-shell UTA feature. This is most likely due to an underestimation of the dielectronic recombination O VI and discuss its possible origin. 6) The lower limit on the distance of the absorbing gas in NGC3783 is between 0.2 and 3.2 pc, depending of the specific ionization component. The constant pressure assumption imposes an upper limit of about 25 pc on the distance of the least ionized gas from the central sourec.

Netzer, Hagai↗
NASA NTRS2003

Dissociative Ionization and Product Distributions of Benzene and Pyridine by Electron Impact

We report a theoretical study of the dissociative ionization (DI) and product distributions of benzene (C6H6) and pyridine (C5H5N) from their low-lying ionization channels. Our approach makes use of the fact that electronic motion is much faster than nuclear motion allowing DI to be treated as a two-step process. The first step is the electron-impact ionization resulting in an ion with the same nuclear geometry as the neutral molecule. In the second step, the nuclei relax from the initial geometry and undergo unimolecular dissociation. For the ionization process we use the improved binary-encounter dipole (iBED) model [W.M. Huo, Phys. Rev. A64,042719-I (2001)]. For the unimolecular dissociation, we use multiconfigurational self-consistent field (MCSCF) methods to determine the steepest descent pathways to the possible product channels. More accurate methods are then used to obtain better energetics of the paths which are used to determine unimolecular dissociation probabilities and product distributions. Our analysis of the dissociation products and the thresholds of their productions for benzene are compared with the recent dissociative photoionization meausurements of benzene by Feng et al. [R. Feng, G. Cooper, C.E. Brion, J. Electron Spectrosc. Relat. Phenom. 123,211 (2002)] and the dissociative photoionization measurements of pyridine by Tixier et al. [S. Tixier, G. Cooper, R. Feng, C.E. Brion, J. Electron Spectrosc. Relat. Phenom. 123,185 (2002)] using dipole (e,e+ion) coincidence spectroscopy.

Dateo, Christopher E.↗
NASA NTRS2003

Dissociative Ionization of Aromatic and Heterocyclic Molecules

Space radiation poses a major health hazard to humans in space flight. The high-energy charged particles in space radiation ranging from protons to high atomic number, high-energy (HZE) particles, and the secondary species they produce, attack DNA, cells, and tissues. Of the potential hazards, long-term health effects such as carcinogenesis are likely linked to the DNA lesions caused by secondary electrons in the 1 - 30 eV range. Dissociative ionization (DI) is one of the electron collision processes that can damage the DNA, either directly by causing a DNA lesion, or indirectly by producing radicals and cations that attack the DNA. To understand this process, we have developed a theoretical model for DI. Our model makes use of the fact that electron motion is much faster than nuclear motion and assumes DI proceeds through a two-step process. The first step is electron-impact ionization resulting in a particular state of the molecular ion in the geometry of the neutral molecule. In the second step the ion undergoes unimolecular dissociation. Thus the DI cross section sigma(sup DI)(sub a) for channel a is given by sigma(sup DI)(sub a) = sigma(sup I)(sub a) P(sub D) with sigma(sup I)(sub a) the ionization cross section of channel a and P(sub D) the dissociation probability. This model has been applied to study the DI of H2O, NH3, and CH4, with results in good agreement with experiment. The ionization cross section sigma(sup I)(sub a) was calculated using the improved binary encounter-dipole model and the unimolecular dissociation probability P(sub D) obtained by following the minimum energy path determined by the gradients and Hessians of the electronic energy with respect to the nuclear coordinates of the ion. This model is used to study the DI from the low-lying channels of benzene and pyridine to understand the different product formation in aromatic and heterocyclic molecules. DI study of the DNA base thymine is underway. Solvent effects will also be discussed.

Huo, Winifred M.↗