skip to main content


Title: The role of negative hydroxyl ions in the electron generation and breakdown during plasma formation in liquid water
The role of negative hydroxyl ions in liquid-phase plasma discharge formation is investigated using an inhouse modeling framework. Two tunneling sources for electrons are considered—tunneling ionization of water molecules and tunneling detachment of negative hydroxyl ions together with additional reaction steps. The simulations are conducted for a needle-like powered electrode with two different nanosecond rise time voltage profiles—a linear and an exponential rise. Both the profiles have a maximum voltage of 15 kV. The predictions show that the electron detachment, which has a much lower threshold energy requirement, provides a stream of electrons at low applied voltage during the initial rise time. The electrical forces from the electron detachment process generate stronger compression but a weaker expansion regime in the liquid resulting in ∼40% increase in the density and only ∼1% decrease. The electron detachment tunneling process is found to be not limited by the electric field, but rather by the availability of negative hydroxyl ions in the system and ceases when these ions are depleted. The tunnel ionization of water molecules forms the electron wave at a higher applied voltage, but the resulting peak electron number density is typically six orders of magnitude larger than the detachment tunneling. The higher electron number density allows the recycling of depleted negative hydroxyl ions in the system and can reestablish tunneling detachment. In addition, the system experiences a larger variation in density; specifically, a decrease in density due to tunnel ionization. The prediction also shows that irrespective of the initial electron sources (i.e. tunnel ionization or tunnel detachment) the reduced electric field is not sufficient enough to allow electron impact ionization to be active and make a significant contribution. Path flux analysis is conducted to determine the kinetics responsible for the recycling of the negative hydroxyl ions.  more » « less
Award ID(s):
1707282
NSF-PAR ID:
10273200
Author(s) / Creator(s):
;
Date Published:
Journal Name:
Plasma sources science technology
Volume:
30
ISSN:
1361-6595
Page Range / eLocation ID:
065025
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract Self-Organized Patterns (SOPs) at plasma-liquid interface in atmospheric pressure plasma discharges refer to the formation of intricate and puzzling structures due to the interplay of electrodynamic and hydrodynamic processes. Studies conducted to date have shown that this phenomenon results in the formation of distinctive patterns such as circular ring, star, gear, dots, spikes, etc., and primarily depends on working gas, electrolyte type, gap distance, current, conductivity, etc. However, an adequate understanding of how these patterns change from one type to another is still not available. This study aims to elucidate the influence of initial liquid conductance ( σ i ) on the temporal evolution of SOPs in liquid-anode discharges. The discharge was generated in a pin-to-liquid anode configuration at a constant helium (He) flow rate of 500 sccm and DC applied voltage of 6 kV at a gap distance of 12 mm. Through the gradual increment of σ i from 1.8 μ S to 4820 μ S, we observe that the trend in the evolution of SOPS takes place as solid discs, spikes, dots, rings, double rings, and stars. The continuous formation of reactive species onto the liquid anode in all conductive solutions results in a decrease in pH, an increase in bulk liquid temperature, and an increase in total dissolved solutes, and these have been confirmed through experimental measurements. Observations using optical emission spectroscopy show that the electrons at the plasma-liquid interface participate in the reduction of cations followed by their excitation & ionization due to which electron density as well as emissions from excited species (mainly hydroxyl radicals & excited nitrogen) decrease with time. Our investigation provides experimental evidence on the presence of cations at the plasma-liquid interface required for SOP formation. 
    more » « less
  2. Abstract

    This paper presents a computational model to study ion and electron transportation and current-voltage characteristics inside a methane-oxygen flame. A commercial software is used to develop the model by splitting the simulation into the combustion and electrochemical transportation parts. A laboratory experiment is used to compare the results from the model. The initial and boundary conditions represented in the model are similar to the experimental conditions in the laboratory experiment.

    In the combustion part, the general GRI3.0 mechanism plus three additional ionization reactions are applied and results are then used as input into the electrochemical transportation part. A particular inspection line is created to analyze the results of the electrochemical transportation part. Ion, electron number density, and current density are studied along the interval from −40V to 40V electric potential. The ions are heavier and more difficult to move than electrons. The results show that at both torch and work surfaces charged sheaths are formed and cause three different regions of current-voltage relations.

     
    more » « less
  3. Resonant tunneling diodes (RTDs) have come full-circle in the past 10 years after their demonstration in the early 1990s as the fastest room-temperature semiconductor oscillator, displaying experimental results up to 712 GHz and fmax values exceeding 1.0 THz [1]. Now the RTD is once again the preeminent electronic oscillator above 1.0 THz and is being implemented as a coherent source [2] and a self-oscillating mixer [3], amongst other applications. This paper concerns RTD electroluminescence – an effect that has been studied very little in the past 30+ years of RTD development, and not at room temperature. We present experiments and modeling of an n-type In0.53Ga0.47As/AlAs double-barrier RTD operating as a cross-gap light emitter at ~300K. The MBE-growth stack is shown in Fig. 1(a). A 15-μm-diam-mesa device was defined by standard planar processing including a top annular ohmic contact with a 5-μm-diam pinhole in the center to couple out enough of the internal emission for accurate free-space power measurements [4]. The emission spectra have the behavior displayed in Fig. 1(b), parameterized by bias voltage (VB). The long wavelength emission edge is at  = 1684 nm - close to the In0.53Ga0.47As bandgap energy of Ug ≈ 0.75 eV at 300 K. The spectral peaks for VB = 2.8 and 3.0 V both occur around  = 1550 nm (h = 0.75 eV), so blue-shifted relative to the peak of the “ideal”, bulk InGaAs emission spectrum shown in Fig. 1(b) [5]. These results are consistent with the model displayed in Fig. 1(c), whereby the broad emission peak is attributed to the radiative recombination between electrons accumulated on the emitter side, and holes generated on the emitter side by interband tunneling with current density Jinter. The blue-shifted main peak is attributed to the quantum-size effect on the emitter side, which creates a radiative recombination rate RN,2 comparable to the band-edge cross-gap rate RN,1. Further support for this model is provided by the shorter wavelength and weaker emission peak shown in Fig. 1(b) around = 1148 nm. Our quantum mechanical calculations attribute this to radiative recombination RR,3 in the RTD quantum well between the electron ground-state level E1,e, and the hole level E1,h. To further test the model and estimate quantum efficiencies, we conducted optical power measurements using a large-area Ge photodiode located ≈3 mm away from the RTD pinhole, and having spectral response between 800 and 1800 nm with a peak responsivity of ≈0.85 A/W at  =1550 nm. Simultaneous I-V and L-V plots were obtained and are plotted in Fig. 2(a) with positive bias on the top contact (emitter on the bottom). The I-V curve displays a pronounced NDR region having a current peak-to-valley current ratio of 10.7 (typical for In0.53Ga0.47As RTDs). The external quantum efficiency (EQE) was calculated from EQE = e∙IP/(∙IE∙h) where IP is the photodiode dc current and IE the RTD current. The plot of EQE is shown in Fig. 2(b) where we see a very rapid rise with VB, but a maximum value (at VB= 3.0 V) of only ≈2×10-5. To extract the internal quantum efficiency (IQE), we use the expression EQE= c ∙i ∙r ≡ c∙IQE where ci, and r are the optical-coupling, electrical-injection, and radiative recombination efficiencies, respectively [6]. Our separate optical calculations yield c≈3.4×10-4 (limited primarily by the small pinhole) from which we obtain the curve of IQE plotted in Fig. 2(b) (right-hand scale). The maximum value of IQE (again at VB = 3.0 V) is 6.0%. From the implicit definition of IQE in terms of i and r given above, and the fact that the recombination efficiency in In0.53Ga0.47As is likely limited by Auger scattering, this result for IQE suggests that i might be significantly high. To estimate i, we have used the experimental total current of Fig. 2(a), the Kane two-band model of interband tunneling [7] computed in conjunction with a solution to Poisson’s equation across the entire structure, and a rate-equation model of Auger recombination on the emitter side [6] assuming a free-electron density of 2×1018 cm3. We focus on the high-bias regime above VB = 2.5 V of Fig. 2(a) where most of the interband tunneling should occur in the depletion region on the collector side [Jinter,2 in Fig. 1(c)]. And because of the high-quality of the InGaAs/AlAs heterostructure (very few traps or deep levels), most of the holes should reach the emitter side by some combination of drift, diffusion, and tunneling through the valence-band double barriers (Type-I offset) between InGaAs and AlAs. The computed interband current density Jinter is shown in Fig. 3(a) along with the total current density Jtot. At the maximum Jinter (at VB=3.0 V) of 7.4×102 A/cm2, we get i = Jinter/Jtot = 0.18, which is surprisingly high considering there is no p-type doping in the device. When combined with the Auger-limited r of 0.41 and c ≈ 3.4×10-4, we find a model value of IQE = 7.4% in good agreement with experiment. This leads to the model values for EQE plotted in Fig. 2(b) - also in good agreement with experiment. Finally, we address the high Jinter and consider a possible universal nature of the light-emission mechanism. Fig. 3(b) shows the tunneling probability T according to the Kane two-band model in the three materials, In0.53Ga0.47As, GaAs, and GaN, following our observation of a similar electroluminescence mechanism in GaN/AlN RTDs (due to strong polarization field of wurtzite structures) [8]. The expression is Tinter = (2/9)∙exp[(-2 ∙Ug 2 ∙me)/(2h∙P∙E)], where Ug is the bandgap energy, P is the valence-to-conduction-band momentum matrix element, and E is the electric field. Values for the highest calculated internal E fields for the InGaAs and GaN are also shown, indicating that Tinter in those structures approaches values of ~10-5. As shown, a GaAs RTD would require an internal field of ~6×105 V/cm, which is rarely realized in standard GaAs RTDs, perhaps explaining why there have been few if any reports of room-temperature electroluminescence in the GaAs devices. [1] E.R. Brown,et al., Appl. Phys. Lett., vol. 58, 2291, 1991. [5] S. Sze, Physics of Semiconductor Devices, 2nd Ed. 12.2.1 (Wiley, 1981). [2] M. Feiginov et al., Appl. Phys. Lett., 99, 233506, 2011. [6] L. Coldren, Diode Lasers and Photonic Integrated Circuits, (Wiley, 1995). [3] Y. Nishida et al., Nature Sci. Reports, 9, 18125, 2019. [7] E.O. Kane, J. of Appl. Phy 32, 83 (1961). [4] P. Fakhimi, et al., 2019 DRC Conference Digest. [8] T. Growden, et al., Nature Light: Science & Applications 7, 17150 (2018). [5] S. Sze, Physics of Semiconductor Devices, 2nd Ed. 12.2.1 (Wiley, 1981). [6] L. Coldren, Diode Lasers and Photonic Integrated Circuits, (Wiley, 1995). [7] E.O. Kane, J. of Appl. Phy 32, 83 (1961). [8] T. Growden, et al., Nature Light: Science & Applications 7, 17150 (2018). 
    more » « less
  4. Recent advancement in the switching of perpendicular magnetic tunnel junctions with an electric field has been a milestone for realizing ultra-low energy memory and computing devices. To integrate with current spin-transfer torque-magnetic tunnel junction and spin–orbit torque-magnetic tunnel junction devices, the typical linear fJ/V m range voltage controlled magnetic anisotropy (VCMA) needs to be significantly enhanced with approaches that include new materials or stack engineering. A possible bidirectional and 1.1 pJ/V m VCMA effect has been predicted by using heavily electron-depleted Fe/MgO interfaces. To improve upon existing VCMA technology, we have proposed inserting high work function materials underneath the magnetic layer. This will deplete electrons from the magnetic layer biasing the gating window into the electron-depleted regime, where the pJ/V m and bidirectional VCMA effect was predicted. We have demonstrated tunable control of the Ta/Pd(x)/Ta underlayer's work function. By varying the Pd thickness (x) from 0 to 10 nm, we have observed a tunable change in the Ta layer's work function from 4.32 to 4.90 eV. To investigate the extent of the electron depletion as a function of the Pd thickness in the underlayer, we have performed DFT calculations on supercells of Ta/Pd(x)/Ta/CoFe/MgO, which demonstrate that electron depletion will not be fully screened at the CoFe/MgO interface. Gated pillar devices with Hall cross geometries were fabricated and tested to extract the anisotropy change as a function of applied gate voltage for samples with various Pd thicknesses. The electron-depleted Pd samples show three to six times VCMA improvement compared to the electron accumulated Ta control sample. 
    more » « less
  5. Abstract Non-equilibrium plasmas derive their low temperature reactivity from producing and driving energetic electrons and active species under large electric fields. Therefore, the impact of reactants on the plasma properties including electron number density, electric field, and electron temperature is critical for applications such as plasma methane (CH 4 ) reforming. Due to experimental complexity, electron properties and the electric field are rarely measured together in the same discharge. In this work, we combine time-resolved Thomson scattering and electric field induced second harmonic generation to probe electron temperature, electron density, and electric field strength in a 60 Torr CH 4 /Ar nanosecond-pulsed dielectric barrier discharge while varying the CH 4 mole fraction from 0% to 8%. These measurements are compared to a 1D numerical model to benchmark its predictions and identify areas of uncertainty. Nonlinear coupling between CH 4 addition, electron temperature, electron density, and the electric field was directly observed. Contrary to previous measurements in He, the electron temperature increased with CH 4 mole fraction. This rise in electron temperature is identified as electron heating by residual electric fields that increased with larger CH 4 mole fraction. Moreover, the electron number density has been found to decrease rapidly with the increase of methane mole fraction. Comparison of these measurements with the model yielded better agreement at higher CH 4 mole fractions and with the usage of ab initio calculated Ar electron-impact cross-sections from the B-spline R-matrix database. Furthermore, the calculated plasma properties are shown to be sensitive to the residual surface charge implanted on the quartz dielectric surfaces. Without considering surface charge in the simulations, the calculated electric field profiles agreed well with the measurements, but the electron properties were underpredicted by more than a factor of three. Therefore, measurements of either the electric field or electron properties measurements alone are insufficient to fully validate modeling predictions. 
    more » « less