Investigation of the electron energy distribution function in hollow-cathode glow discharges in nitrogen and oxygen

The mechanism for the formation of the inverse electron distribution function is proposed and realized experimentally in a nitrogen plasma of a hollow-cathode glow discharge. It is shown theoretically and experimentally that, for a broad range of the parameters of an N₂ discharge, it is possible to form a significant dip in the profile of the electron distribution function in the energy range ε=2–4 eV and, accordingly, to produce the inverse distribution with df(ε)/d?>0. The formation of a dip is associated with both the vibrational excitation of N₂ molecules and the characteristic features of a hollow-cathode glow discharge. In such a discharge, the applied voltage drops preferentially across a narrow cathode sheath. In the main discharge region, the electric field E is weak (E<0.1 V/cm at a pressure of about p～0.1 torr) and does not heat the discharge plasma. The gas is ionized and the ionization-produced electrons are heated by a beam of fast electrons (with an energy of about 400 eV) emitted from the cathode. A high-energy electron beam plays an important role in the formation of a dip in the profile of the electron distribution function in the energy range in which the cross section for the vibrational excitation of nitrogen molecules is maximum. A plasma with an inverted electron distribution function can be used to create a population inversion in which more impurity molecules and atoms will exist in electronically excited states.     

Electrostatic electron vortex structure in a plasma in an external magnetic field

A previously developed method for describing vortex structures is used to construct electrostatic vortices in a plasma in an external magnetic field. An equation for the radial electric field that gives rise to azimuthal electron drift in crossed electric (E _r ) and magnetic (B _z ) fields is derived without allowance for the magnetic field of the electron currents. Two types of the resulting electrostatic vortex structures with a positive and a negative electric potential at the axis are analyzed. The results obtained are compared with experimental data on vortex structures.     

Optical characteristics and parameters of the plasma of a barrier discharge excited in a mixture of mercury dibromide vapor with nitrogen and helium

Results are presented from experimental and theoretical studies of the optical characteristics and parameters of the plasma of an atmospheric-pressure barrier discharge excited in a HgBr₂: N₂: He mixture, which was used as the working medium of a small-size (with a radiation area of 8 cm²) exciplex gas-discharge radiation source. The mean radiation power of 87 mW was achieved at the radiation wavelength λ_max = 502 nm. The electron energy distribution function, the transport characteristics, the specific energy lost in the processes involving electrons, the electron temperature and density, and the rate constants of elastic and inelastic electron scattering by the components of the working mixture were calculated as functions of the reduced field E/N. The plasma of a discharge excited in a HgBr₂: N₂: He mixture can be used as the working medium of a small-size blue-green radiation source. Such a source can find application in biotechnology, photonics, and medicine and can also be used to manufacture gas-discharge display panels.     

Effect of the nonlocal nature of the electron energy spectrum on the dissociation of oxygen molecules in a discharge

The effect of the nonlocal nature of the electron distribution function on the dissociation rate of oxygen molecules in a dc glow discharge is studied. The concentration of oxygen atoms and the probability of their loss at the discharge tube wall are measured as functions of the discharge parameters by means of the timeresolved actinometric method involving argon atoms. An analysis of the measurement data in terms of both a discharge model in which the effect of the nonlocal nature of the electron energy spectrum is taken into account and a model in which this effect is ignored makes it possible to thoroughly investigate the balance of oxygen atoms in the discharge. The production rate of O(³P) atoms and their concentration in the plasma are calculated with allowance for the nonlocal nature of the electron energy distribution function. The calculated values agree well with the experimental data and differ substantially from those obtained using a spatially homogeneous distribution function.     

Method of simplified kinetic equation for simulating the electron kinetics in dense gases in strong electric fields

A method for calculating the electron kinetics in dense gases in strong electric fields is developed. The method differs from the forward-backward approximation proposed by Ra?zer and Shne?der for “high-energy” electrons in that it introduces the effective cosines of the scattering angles with respect to the electric force, μ⁺(?, E) and μ^?(?, E), which are different from +1 and ?1, as in the forward-backward approximation. The method was implemented numerically for atmospheric-pressure helium and molecular nitrogen for fields in the ranges 10–200 and 50–800 kV/cm, respectively. The cosines μ⁺(?, E) and μ^?(?, E), the frequency of “fatal” collisions making high-energy electrons to pass from an acceleration regime to a deceleration one, and the rate at which the electrons leave the low-energy reservoir with energies of ≤15 eV for nitrogen and of ≤20 eV for helium are calculated.     

Influence of the plasma density and heating power on the intensity of electron cyclotron emission in the L-2M stellarator

Results are presented from experiments on studying the plasma behavior in the L-2M stellarator in regimes with a high power deposition in electrons during electron cyclotron heating at the second harmonic of the electron gyrofrequency (X mode) at heating powers of P _in=120–400 kW and average plasma densities from n _e≤3×10¹⁹ to 0.3×10¹⁹ m^?3. It is shown that, as the plasma density decreases and the heating power increases, the electron cyclotron emission spectrum is modified; this may be attributed to a deviation of the electron energy distribution from a Maxwellian and the generation of suprathermal electrons. At low plasma densities, the emission intensity at the second harmonic of the electron gyrofrequency increases, whereas the plasma energy measured by diamagnetic diagnostics does not increase. This poses the question of the correctness of determining the plasma electron temperature by electron cyclotron emission diagnostics under these conditions.     

Comparison between an event-by-event Monte Carlo code,NOREC, and ETRAN for electron scaled point kernels between 20 keV and 1 MeV

An event-by-event Monte Carlo code called NOREC, a substantially improved version of the Oak Ridge electron transport code (OREC), was released in 2003, after a number of modifications to OREC. In spite of some earlier work, the characteristics of the code have not been clearly shown so far, especially for a wide range of electron energies. Therefore, NOREC was used in this study to generate one of the popular dosimetric quantities, the scaled point kernel, for a number of electron energies between 0.02 and 1.0 MeV. Calculated kernels were compared with the most well-known published kernels based on a condensed history Monte Carlo code, ETRAN, to show not only general agreement between the codes for the electron energy range considered but also possible differences between an event-by-event code and a condensed history code. There was general agreement between the kernels within about 5% up to 0.7 r/r ₀ for 100 keV and 1 MeV electrons. Note that r/r ₀ denotes the scaled distance, where r is the radial distance from the source to the dose point and r ₀ is the continuous slowing down approximation (CSDA) range of a mono-energetic electron. For the same range of scaled distances, the discrepancies for 20 and 500 keV electrons were up to 6 and 12%, respectively. Especially, there was more pronounced disagreement for 500 keV electrons than for 20 keV electrons. The degree of disagreement for 500 keV electrons decreased when NOREC results were compared with published EGS4/PRESTA results, producing similar agreement to other electron energies.     

Diffuse and constricted modes of a dc discharge in neon: Simulation of the hysteresis transition

Results are presented from theoretical studies of high-pressure (～100 Torr) dc discharges in neon. The diffuse and constricted discharge modes are studied using a model including the equation of balance for charged and excited particles, heat conduction equations for the neutral gas and plasma electrons, and Poisson’s equation for the radial electric field at a fixed total discharge current. A specific feature of the constricted mode in the investigated range of low fields and high degrees of ionization is that the excitation and ionization rates in the center of the discharge tube and at the periphery differ by several orders of magnitude. This implies that, in the constricted mode, the region where the electron energy distribution function is Maxwellian due to electron-electron collisions may adjoin the region (beyond the constriction zone) where the high-energy part of the distribution function is depleted. The hysteresis transition between the diffuse and constricted modes is analyzed. A transition from the constricted to the diffuse mode can be regarded as a manifestation of the nonlocal character of the formation of the electron distribution function, specifically, the diffusion of high-energy electrons capable of producing gas ionization from the central (constricted) region toward the periphery. The nonlocal formation of the distribution function is described by a nonlocal kinetic equation accounting for electron-electron collisions and electron transport along the radius of the discharge tube. Since only high-energy electrons produce gas ionization, the effect of the nonlocal formation of the electron distribution function is taken into account by introducing the effective temperature of the high-energy part of the distribution function and solving the equation for the radial profile of the high-energy part of the distribution function. This approach allows one to approximately take into account the nonlocal character of the electron distribution function without substantial expenditure of computer resources. The nonlocal model makes it possible to numerically simulate the hysteresis transition between the diffuse and constricted modes, which is impossible in the local approximation.     

Measurements and calculations of the electron distribution function in the electronegative plasma of a hollow-cathode discharge in N<Subscript>2</Subscript>: SF<Subscript>6</Subscript> mixtures

Results are presented from measurements and numerical calculations of the electron energy distribution function in the plasma of a hollow-cathode glow discharge in N₂: SF₆ mixtures. It is shown that, when a small amount of SF₆ is added to nitrogen, the number of electrons in the inverse region of the distribution function (2–6 eV) increases by about one order of magnitude. As the electric field in plasma increases to ≈0.5 V/cm, the dip in the distribution function disappears and the inversion region vanishes.     

Condensed-history Monte Carlo simulation of the dosimetric distribution of electron microbeam

To evaluate the dose distributions of an electron microbeam and to help optimization of its design, the condensed-history (CH) Monte Carlo simulation algorithm implemented in the Geant4 toolkit was selected as an alternative to the conventionally used analog algorithm. Compared to the analog algorithm, the CH algorithm is cheaper and less limited by the lack of cross-sections. And, with a properly chosen production cut for secondaries, its accuracy for the problems of microdosimetry is satisfactory. In this work, calculations of the single-event (imparted energy ɛ) size distribution f ₁(ɛ) is described, for compartments in the Orlando electron micro beam. The results agree well with those obtained by use of the analog algorithm and reported in the literature. It is shown that substituting tissue with water in HeLa cells, and replacing Mylar with water of the same mass stopping power in the substrate, makes little difference. Additionally, the neighbor-to-target ratio of average event size R _NT has been calculated as a function of the incident energy of the electrons. Comparison with analog results reported in the literature suggests that the average event size in neighbors, and hence the neighbor-to-target ratio, is sensitive to the selection of the energy threshold. Finally, the effect of finite beam radius on the event size distribution and the neighbor-to-target ratio has also been studied. All results presented suggest the condensed-history technique to provide an efficient and valuable alternative tool in the design of electron microbeams.     

Steady-state vortex structure in a plasma with a strong magnetic field

A study is made of nonquasineutral vortex structures in a plasma with a magnetic field B _z in which the charges separate on a spatial scale equal to the magnetic Debye radius r _B=B _z/4πen _e. The electric field arising due to charge separation leads to radial expansion of the ions, thereby destroying the initial electron vortex. It is shown that the ion pressure gradient stops ion expansion in a nonquasineutral electron vortex and gives rise to a steady structure with a characteristic scale on the order of r _B. With the electron inertia taken into account in the hydrodynamic approximation, the magnetic vortex structure in a hot plas mamanifests itself in the appearance of a “hole” in the plasma density.     

Electron leakage mechanism in a plasma-filled magnetically insulated transmission line

It is shown that relativistic electron current can propagate across the magnetic field B ₀ over a distance d much larger than the electron gyroradius, r ₀ ? m _e v _z c/(eB ₀) ? d. This current is driven by the Hall electric field, which is generated on a spatial scale equal to the magnetic Debye radius r _B = B ₀/(4πen _e) and causes the electrons to drift in crossed electric and magnetic fields. For a plane equilibrium current configuration, analytic profiles of the electron velocity and electron density are calculated and the electric and magnetic fields are determined. The results obtained are used to explain electron leakages in magnetically insulated transmission lines filled with a plasma expanding from the electrodes. Equations describing an equilibrium configuration of the ions and electrons that drift simultaneously across a strong magnetic field are derived.     

Effect of the radial electric field on the results of measurements of the poloidal rotation of a tokamak plasma by charge exchange recombination spectroscopy

Abtract The effect of the radial electric field E _r on the results of measurements of the poloidal rotation of a tokamak plasma by charge exchange recombination spectroscopy is considered. It is shown that the emission line shift arising from the finite lifetime of the excited state of the ions is proportional to E _r. For helium ions, the maximum shift corresponds to the poloidal rotation velocity, which is about one-third of the drift velocity in the crossed radial electric (E _r) and toroidal magnetic (B _t) fields. __________ Translated from Fizika Plazmy, Vol. 27, No. 11, 2001, pp. 1050–1052. Original Russian ? 2001 by Romannikov, Chernobai.     

Application of Polyelectrolyte Theory to the Study of the B-Z Transition in DNA (1)

Abstract

We have used the polyelectrolyte theory to study the ionic strength dependence of the B-Z equilibrium in DNA. A DNA molecule is molded as an infinitely long continuously charged cylinder of radius a with reduced linear charge density q. The parameters a and q for the B and Z forms were taken from X-ray data: a _B = 1nm, q _B = 4.2, a _z = 0.9 nm and q _z = 3.9. A simple theory shows that at low ionic strengths (when Debye screening length r _D>>a) the electrostatic free energy difference F ^el _Bz = F ^el _Z - F ^el _B increases with increasing ionic strength since q _B>q_z. At high ionic strengths (when r _D<<a) the F ^el _BZ would go on growing with increasing ionic strength if the inequality q _B/^a _B<q_z/a _z were valid. In the converse case when q _z/q _B<a_z/a _B the F ^el _BZ value decreases with increasing salt concentration at high ionic strength. Since X-ray data correspond to the latter case, theory predicts that the F ^el _BZ value reaches a maximum at an intermediate ionic strength of about 0.1 M (where r _D～a). We also performed rigorous calculations based on the Poisson-Boltzmann equation. These calculations have confirmed the above criterion of nonmonotonous behaviour of the F ^el _BZ value as a function of ionic strength. Different theoretical predictions for the B-Z transition in linear and superhelical molecules are discussed. Theory predicts specifically that at a very low ionic strength the Z form may prove to be more stable than the B form. Thus, one can observe the Z-B-Z transition with increasing ionic strength. In the light of our theoretical findings we discuss numerous experimental data on the B-Z transition in linear and superhelical DNA.     

Applications of the ETS-NOCV method in descriptions of chemical reactions

The present study characterizes changes in the electronic structure of reactants during chemical reactions based on the combined charge and energy decomposition scheme, ETS-NOCV (extended transition state–natural orbitals for chemical valence). Decomposition of the activation barrier, ΔE ^#, into stabilizing (orbital interaction, ΔE _orb, and electrostatic, ΔE _elstat) and destabilizing (Pauli repulsion, ΔE _Pauli, and geometry distortion energy, ΔE _dist) factors is discussed in detail for the following reactions: (I) hydrogen cyanide to hydrogen isocyanide, HCN → CNH isomerization; (II) Diels-Alder cycloaddition of ethene to 1,3-butadiene; and two catalytic processes, i.e., (III) insertion of ethylene into the metal-alkyl bond using half-titanocene with phenyl-phenoxy ligand catalyst; and (IV) B–H bond activation catalyzed by an Ir-containing catalyst. Various reference states for fragments were applied in ETS-NOCV analysis. We found that NOCV-based deformation densities (Δρ _i) and the corresponding energies ΔE _orb(i) obtained from the ETS-NOCV scheme provide a very useful picture, both qualitatively and quantitatively, of electronic density reorganization along the considered reaction pathways. Decomposition of the barrier ΔE^# into stabilizing and destabilizing contributions allowed us to conclude that the main factor responsible for the existence of positive values of ΔE ^# for all processes (I, II, III and IV) is Pauli interaction, which is the origin of steric repulsion. In addition, in the case of reactions II, III and IV, a significant degree of structural deformation of the reactants, as measured by the geometry distortion energy, plays an important role. Depending on the reaction type, stabilization of the transition state (relatively to the reactants) originating either from the orbital interaction term or from electrostatic attraction can be of vital importance. Finally, use of the ETS-NOCV method to describe catalytic reactions allows extraction of information on the role of catalysts in determination of ΔE ^#.     

Estimation of the electron density and radiative energy losses in a calcium plasma source based on an electron cyclotron resonance discharge

The parameters of a calcium plasma source based on an electron cyclotron resonance (ECR) discharge were calculated. The analysis was performed as applied to an ion cyclotron resonance system designed for separation of calcium isotopes. The plasma electrons in the source were heated by gyrotron microwave radiation in the zone of the inhomogeneous magnetic field. It was assumed that, in such a combined trap, the energy of the extraordinary microwave propagating from the high-field side was initially transferred to a small group of resonance electrons. As a result, two electron components with different transverse temperatures—the hot resonance component and the cold nonresonance component—were created in the plasma. The longitudinal temperatures of both components were assumed to be equal. The entire discharge space was divided into a narrow ECR zone, where resonance electrons acquired transverse energy, and the region of the discharge itself, where the gas was ionized. The transverse energy of resonance electrons was calculated by solving the equations for electron motion in an inhomogeneous magnetic field. Using the law of energy conservation and the balance condition for the number of hot electrons entering the discharge zone and cooled due to ionization and elastic collisions, the density of hot electrons was estimated and the dependence of the longitudinal temperature T _e∥ of the main (cold) electron component on the energy fraction β lost for radiation was obtained.     

Nonlocal equation for the symmetric part of the electron distribution function in an inhomogeneous plasma

A nonlocal kinetic equation is derived for the symmetric part of the distribution function of suprathermal electrons. It is shown that Albritton equations are merely local approximations to the total kinetic equation. Even in the simplest situation, the local approximations of the nonlocal effects are impossible to construct because of the interdependence of the variables. A self-similar solution to the equations under study is proposed.     

Modelling and off line optimization of batch gluconic acid fermentation

Summary The role of mathematical modelling and off line optimization for a batch fermentation process is described. The fermentation of gluconic acid by Acetobacter suboxydans ATCC 621 was studied. The model is based on a series of batch experiments in which the temperature was the only variable. The differential equations of the models were derived from these experiments to give the kinetic parameters and the parametric models varying with the temperature. The fermentation was optimized using Pontryagin's maximum principle. This gave the temperature profile of fermentation.Abbreviations x, g, l, S, c, P The concentration of cell mass, glucose, lactone, gluconic acid, 5-ketogluconic acid and total acidic products respectively - r₁, r₂ E₁ or E₂ enzyme in complex/total E₁ or E₂ enzyme content - a, b E₁ and E₂ enzyme content of unit quantity of biomass - k_i and K_j Rate constants - _max Maximum specific growth rate - y_x z₁=r₁x z₂=r₂x Yield coefficient of biomass with respect to growth on glucose - z₁ r₁x - z₂ r₂x