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 more » 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. « less
Authors:
;
Award ID(s):
1707282
Publication Date:
NSF-PAR ID:
10273200
Journal Name:
Plasma sources science technology
Volume:
30
Page Range or eLocation-ID:
065025
ISSN:
1361-6595
Sponsoring Org:
National Science Foundation
More Like this
  1. 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.more »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).« 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.

  3. Decavanadate (V 10 O 28 6− or V10) is a paradigmatic member of the polyoxidometalate (POM) family, which has been attracting much attention within both materials/inorganic and biomedical communities due to its unique structural and electrochemical properties. In this work we explored the utility of high-resolution electrospray ionization (ESI) mass spectrometry (MS) and ion exclusion chromatography LC/MS for structural analysis of V10 species in aqueous solutions. While ESI generates abundant molecular ions representing the intact V10 species, their isotopic distributions show significant deviations from the theoretical ones. A combination of high-resolution MS measurements and hydrogen/deuterium exchange allows these deviations to be investigated and interpreted as a result of partial reduction of V10. While the redox processes are known to occur in the ESI interface and influence the oxidation state of redox-active analytes, the LC/MS measurements using ion exclusion chromatography provide unequivocal evidence that the mixed-valence V10 species exist in solution, as extracted ion chromatograms representing V10 molecular ions at different oxidation states exhibit distinct elution profiles. The spontaneous reduction of V10 in solution is seen even in the presence of hydrogen peroxide and has not been previously observed. The susceptibility to reduction of V10 is likely to be shared bymore »other redox active POMs. In addition to the molecular V10 ions, a high-abundance ionic signal for a V 10 O 26 2− anion was displayed in the negative-ion ESI mass spectra. None of the V 10 O 26 cations were detected in ESI MS, and only a low-abundance signal was observed for V 10 O 26 anions with a single negative charge, indicating that the presence of abundant V 10 O 26 2− anions in ESI MS reflects gas-phase instability of V 10 O 28 anions carrying two charges. The gas-phase origin of the V 10 O 26 2− anion was confirmed in tandem MS measurements, where mild collisional activation was applied to V10 molecular ions with an even number of hydrogen atoms (H 4 V 10 O 28 2− ), resulting in a facile loss of H 2 O molecules and giving rise to V 10 O 26 2− as the lowest-mass fragment ion. Water loss was also observed for V 10 O 28 anions carrying an odd number of hydrogen atoms ( e.g. , H 5 V 10 O 28 − ), followed by a less efficient and incomplete removal of an OH˙ radical, giving rise to both HV 10 O 26 − and V 10 O 25 − fragment ions. Importantly, at least one hydrogen atom was required for ion fragmentation in the gas phase, as no further dissociation was observed for any hydrogen-free V10 ionic species. The presented workflow allows a distinction to be readily made between the spectral features revealing the presence of non-canonical POM species in the bulk solution from those that arise due to physical and chemical processes occurring in the ESI interface and/or the gas phase.« less
  4. Pinhole-free ultrathin films of polyelectrolyte complex assembled using layer-by-layer deposition were used to evaluate electron transfer from a redox species in solution to an electrode over the distance range 1 to 9 nm. Over this thickness, the polyelectrolytes employed wet the surface and the polymer molecules flatten to less than their equilibrium size in 3-dimensions. A decay constant β for current as a function of distance of about 0.3 nm-1 placed this system in the regime expected for multistep hopping versus a one-step tunneling event. Discreet hopping sites within the films were identified as ferrocyanide ions with an equilibrium concentration of 0.032 M and an average separation of 3.7 nm. The Butler-Volmer (BV) expression for electron transfer as a function of overpotential was modified by distributing the applied voltage evenly amongst the hopping sites. This modified BV expression fit both the distance dependence and the applied potential dependence well, wherein the only freely adjustable parameter was the electron transfer coefficient. The finding that β is simply the inverse of the hopping range is consistent with previous conclusions that electrons within conjugated molecule sites are delocalized, or, for non-conjugated systems, spread over more than one repeat unit by lattice distortions
  5. Excess electrons in liquid acetonitrile are of particular interest because they exist in two different forms in equilibrium: they can be present as traditional solvated electrons in a cavity, and they can form some type of solvated molecular anion. Studies of small acetonitrile cluster anions in the gas phase show two isomers with distinct vertical detachment energies, and it is tempting to presume that the two gas-phase cluster anion isomers are precursors of the two excess electron species present in bulk solution. In this paper, we perform DFT-based ab initio molecular dynamics simulations of acetonitrile cluster anions to understand the electronic species that are present and why they have different binding energies. Using a long-range-corrected density functional that was optimally tuned to describe acetonitrile cluster anion structures, we have theoretically explored the chemistry of (CH3CN)n¯ cluster anions with sizes n=5,7 and 10. Since the temperature of the experimental cluster anions is not known, we performed two sets of simulations that investigated how the way in which the cluster anions are prepared affects the excess electron binding motif: one set of simulations simply attached excess electrons to neutral (CH3CN)n clusters, providing little opportunity for the clusters to relax in the presencemore »of the excess electron, while the other set allowed the cluster anions to thermally equilibrate near room temperature. We find that both sets of simulations show three distinct electron binding motifs: electrons can attach to the surface of the cluster (dipole-bound) or be present as either solvated monomer anions, CH3CN¯, or as solvated molecular dimer anions, (CH3CN)2¯. All three species have higher binding energies at larger cluster sizes. Thermal equilibration strongly favors the formation of the valence-bound molecular anions relative to surface-bound excess electrons, and the dimer anion becomes more stable than the monomer anion and surface-bound species as the cluster size increases. The calculated photoelectron spectra from our simulations in which there was poor thermal equilibration are in good agreement with experiment, suggesting assignment of the two experimental cluster anion isomers as the surface-bound electron and the solvated molecular dimer anion. The simulations also suggest that the shoulder seen experimentally on the low-energy isomer's detachment peak is not part of a vibronic progression but instead results from molecular monomer anions. Nowhere in the size range that we explore do we see evidence for a non-valence, cavity-bound interior-solvated electron, indicating that this species is likely only accessible at larger sizes with good thermal equilibration.« less