Active remote detection of radioactivity based on electromagnetic signatures

Abstract

A system for the active remote detection of radioactivity from a target of interest includes a first laser source for generating an ionizing laser beam when remotely directed on a radioactive target of interest, a second laser source for generating a laser probe beam on the radioactive target of interest, and a spectrometer configured to measure the frequency modulation of the probe beam caused by the ionization from the radioactive target of interest.

Claims

1. A method for active remote detection of radioactivity from a radioactive target of interest that emits radiation causing ions to form in air surrounding the radioactive target of interest, the method comprising: generating a laser beam using a first laser source, wherein the laser beam is configured to photo-detach the ions formed by the radioactive target of interest; generating a laser probe beam using a second laser source; directing the laser beam and the laser probe beam onto the radioactive target of interest; measuring, based on a time evolution of electron density near the radioactive target of interest, a frequency modulation of the laser probe beam; determining, based on the frequency modulation, an electromagnetic signature of the radioactive target of interest; and identifying, based on the electromagnetic signature, a type of radioactive material in the radioactive target of interest.

2. The method of claim 1, wherein the first laser source has a peak intensity of I.sub.peak=160 GW/cm.sup.2 and a pulse duration of τ.sub.laser=1 nsec.

3. The method of claim 1, wherein the laser probe beam has a frequency in the range of 90 GHz to 110 GHz.

4. The method of claim 1, wherein identifying the type of radioactive material in the radioactive target of interest further comprises: detecting a generation of a spark/air breakdown initiated by the laser beam.

5. The method of claim 1, wherein the laser beam initiates an avalanche ionization in air surrounding the radioactive target of interest.

6. The method of claim 5, wherein the avalanche ionization initiates an increase in electron density in the air surrounding the radioactive target of interest.

7. The method of claim 6, wherein the increase in electron density induces the frequency modulation.

8. The method of claim 1, wherein the measured frequency modulation corresponds to an increase in the electron density over time.

9. A system for the active remote detection of radioactivity from a radioactive target of interest, the system comprising: a first laser source configured to generate a laser beam, wherein the laser beam is configured to photo-detach ions formed by the radioactive target of interest; a second laser source configured to generate a laser probe beam directed on the radioactive target of interest; and a spectrometer configured to measure, based on a time evolution of electron density near the radioactive target of interest, a frequency modulation of the laser probe beam caused by an increase in the electron density, initiated by the laser beam, in air surrounding the radioactive target of interest.

10. The system of claim 9, wherein a radioactive composition of the radioactive target of interest is determined by: measuring a frequency modulation of the laser probe beam, determining, based on the frequency modulation, an electromagnetic signature of the radioactive target of interest, and identifying, based on the electromagnetic signature, the radioactive target of interest as having a specific type of radioactive material.

11. The system of claim 9, wherein the first laser source has a peak intensity of I.sub.peak=160 GW/cm.sup.2 and a pulse duration of τ.sub.laser=1 nsec.

12. The system of claim 9, wherein the laser probe beam has a frequency in the range of 90 GHz to 110 GHz.

13. The system of claim 9, wherein the laser beam is configured to initiate an avalanche ionization in the air surrounding the radioactive target of interest.

14. The system of claim 13, wherein the avalanche ionization initiates the increase in electron density.

15. The system of claim 14, wherein the frequency modulation of the laser probe beam correlates to a type of radioactive material.

16. The system of claim 9, wherein the spectrometer is configured to detect frequencies in the range of +/−10% of a frequency of the laser probe beam.

17. The system of claim 9, wherein the spectrometer is further configured to: detect a generation of a spark/air breakdown initiated by the laser beam.

18. A system for the active remote detection of radioactivity from a radioactive source that emits gamma rays causing O.sub.2.sup.− ions to form in air surrounding the radioactive source, the system comprising: a first laser source configured to generate a laser beam, wherein the laser beam has a peak intensity and a pulse duration sufficient to photo-detach electrons from the O.sub.2.sup.− ions and initiate avalanche ionization in air surrounding the radioactive source when the laser beam is directed on the radioactive source from a selected distance; a second laser source configured to generate a laser probe beam modulated by the avalanche ionization; and a spectrometer configured to measure, based on a time evolution of electron density near the radioactive target of interest, a frequency modulation of the laser probe beam caused by the avalanche ionization.

19. The system of claim 18, wherein a radioactive composition of the radioactive source is determined by: measuring a frequency modulation of the laser probe beam, determining, based on the frequency modulation, an electromagnetic signature of the radioactive source, and identifying, based on the electromagnetic signature, the radioactive source as having a specific type of radioactive material.

20. The system of claim 19, wherein the first laser source has a peak intensity of I.sub.peak=160 GW/cm.sup.2 and a pulse duration of τ.sub.laser=1 sec.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

(1) FIG. 1 is a schematic illustration of the active remote radioactivity detection concept according to the invention;

(2) FIG. 2 is a graph showing the radiation enhancement factor plotted as a function of distance from the radioactive source for 1 mg and 10 mg of .sup.60Co according to the invention;

(3) FIG. 3A is a graph showing electron density as a function of time in the absence of any external radioactivity, where the laser parameters are λ=1 μm, I.sub.peak=160 GW/cm.sup.2, τ.sub.laser=1 nsec and FIG. 3B is a graph showing electron density versus time in the presence of radioactivity, with the same laser parameters as in FIG. 3A, according to the invention; and

(4) FIG. 4 is a graph showing the fractional frequency shift Δω/ω.sub.o[%] versus time in the presence of radioactive material α.sub.rad=10.sup.3 at the probe interaction distance of L=10 cm, with the same laser parameters as in FIG. 3A, according to the invention; and

(5) FIG. 5 is a graph showing the fractional frequency shift Δω/ω.sub.o[%] versus time and probe interaction distance L in the presence of radioactive material (α.sub.rad=10.sup.3) according to the invention with the laser parameters the same as in FIG. 3a.

DETAILED DESCRIPTION OF THE INVENTION

Definition

(6) The term “electromagnetic signature” as used herein means, for example, the identifiable correlation between the modulated laser probe beam and the specific activity (radiation level) attributable to a particular radioactive material, as discussed below.

(7) Referring now to FIG. 1 that schematically illustrates the detection concept of the invention, detection system 100 includes an ionizing laser source 102 for generating a photo-detaching and ionizing beam 104 and a second laser source 106 for generating a laser probe beam 108 each from a selected distance onto a target of interest 110. When the target 110 is a radioactive material it emits gamma rays that ionize the surrounding air. The ionized (liberated) electrons rapidly attach to oxygen molecules forming O.sub.2.sup.− ions. The density of O.sub.2.sup.− around radioactive material can be several orders of magnitude greater than background levels. The elevated population of O.sub.2.sup.− extends several meters around the radioactive material. Electrons are easily photo-detached from O.sub.2.sup.− ions by laser radiation. The photo-detached electrons, in the presence of laser radiation, initiate avalanche ionization which results in a rapid increase in electron density. The rise in electron density induces a frequency modulation on a probe beam that correlates to and identifies the particular radioactive material based on its unique specific activity, thereby providing an electromagnetic signature.

(8) The detector: A spectrometer, capable of detecting frequencies in the range of +/−10% of the frequency of the probe laser beam.

(9) The selection of the distance of lasers 102 and 106 from the target 110 may be based on a number of factors such as prevailing atmospheric conditions.

(10) Propagation of high power short pulse lasers in the atmosphere over extended distances has been analyzed and experimentally characterized, e.g. as described in P. Sprangle, J. R. Peñano, and B. Hafizi, Phys. Rev. E 66, 046418 (2002). Since the negative ions produced by the radioactive material have a low ionization potential (0.46 eV) they can be photo-detached by laser radiation (˜0.8-1 μm). The invention's approach is based on the probe beam 108 undergoing a frequency modulation while propagating in a temporally increasing electron density. The frequency modulation on the probe beam 108 becomes a spectral signature for the presence of radioactive material.

(11) Radiation Enhancement Factor

(12) A gamma ray propagating through matter can interact through several processes, including Rayleigh scattering, photoelectric and Compton effects, pair production, and so forth. In air, photoelectric absorption dominates at low photon energies (<25 keV) while at high energies (˜25 keV-3 MeV) Compton processes dominate. As the gamma ray propagates in air it loses energy in a cascading process and its mean free path L.sub.γ decreases. A 1 MeV gamma ray has a mean free path in air of L.sub.γ≈130 m.

(13) The ionization rate due to background (ambient) radioactivity is (d N.sub.e/dt).sub.amb=Q.sub.rad. At or near ground level, the background ionization rate is typically in the range Q.sub.rad˜10-30 pairs/(cm.sup.3-sec). The gamma rays emitted by radioactive material ionize the surrounding air. In the presence of radioactive material the ionization rate (due to only radiation) can be greatly enhanced by a factor α.sub.rad>>1 and (d N.sub.e/dt).sub.radα.sub.radQ.sub.rad. For a radioactive material of mass M.sub.rad the number of disintegrations per second is v.sub.rad=M.sub.radA.sub.rad, where A.sub.rad is the specific activity associated with the material. For example, for .sup.60Co, A.sub.rad=1.1×10.sup.3 Ci/g=4.1×10.sup.13 disintegration/(g-sec). In the case of .sup.60Co each disintegration results in two gammas of energy E.sub.γ,max=1.173 MeV and E.sub.γ,max=1.332 MeV which have a range in air of ˜130 m. In air the high energy gammas generate high energy electrons, via Compton and photoelectric processes, which undergo a cascading process to sufficiently low energy to attach to O.sub.2 molecules forming O.sub.2.sup.− ions. In the cascading process the electrons lose an amount of energy ΔE≈34 eV per collision in air which results in both ionization and electronic excitation. A high energy electron with energy E.sub.e therefore generates ˜E.sub.e/ΔE low energy electrons. An electron having an energy of 1 MeV has a range in air of 4.6 m.

(14) For a small spherical source of radioactivity the steady state density of emitted gamma rays is N.sub.γ=(v.sub.radκ.sub.γ/4πc R.sup.2)exp(−R/L.sub.γ) where R is the distance from the radioactive material, L.sub.γ is the effective range (mean free path) of the gamma rays in air which is a function of the gamma ray energy, E.sub.γ, and κ.sub.γ is the number of gammas emitted per disintegration. The rate of change of electron density is a ∂N.sub.e/∂t≈(α.sub.rad+1) Q.sub.rad+air chemistry and ionization terms, where

(15) ${\alpha }_{\mathrm{rad}}\approx c.Math.{\sigma }_{\gamma -e}.Math.N_{\mathrm{air}}{N}_{\gamma }\frac{.Math.E_e.Math.}{\Delta E}\frac{1}{Q_{\mathrm{rad}}}\approx \frac{v_{\mathrm{rad}}{\kappa }_{\gamma }}{4\pi .Math.L_{\gamma -e}.Math.}\frac{.Math.E_e.Math.}{\Delta E}\frac{1}{Q_{\mathrm{rad}}}\frac{\exp \left(-R/L_{\gamma }\right)}{R^2},$
is the radiation enhancement factor, E.sub.e is the average electron energy, σ.sub.γ-e is the effective average cross section for electron generation by gammas, i.e., Compton absorption and photoelectric processes, L.sub.γ-e=(σ.sub.γ-eN.sub.air).sup.−1 is the average mean free path for electron generation by gammas and N.sub.air=2.7×10.sup.19 cm.sup.−3 is the air density at STP. In the absence of radioactive material α.sub.rad=0. In FIG. 2 the radiation enhancement factor α.sub.rad is plotted as a function of the distance from the radioactive source R. This plot is for samples containing 1 mg and 10 mg of .sup.60Co and indicates that the enhanced level is significant for ranges extending up to several meters. As an example, for E.sub.e=0.5 MeV, M.sub.rad=10 mg, v.sub.rad=M.sub.radA.sub.rad=8.2×10.sup.11 disintegrations/sec, κ.sub.γ=2, R=50 cm, L.sub.γ-e=100 m and Q.sub.rad=20 disintegrations/(cm.sup.3-sec) the radiation enhancement factor is α.sub.rad≈2×10.sup.6 which is far above the background level.

(16) Electron and Ion Density Evolution (Air Chemistry)

(17) To determine the frequency modulation on a probe pulse it is necessary to follow the time evolution of the electron and negative ion density, which are sensitive functions of air chemistry processes (see, e.g., M. Capitelli, C. M. Ferreira, B. F. Gordiets and A. I. Osipov, Plasma Kinetics in Atmospheric Gases (Springer-Verlag, NY 2010)) and electron heating by the laser radiation. The source terms for the electrons include radioactivity, detachment, photo-detachment and photo-ionization, while the loss terms include various attachment and recombination processes including aerosols. The expressions for the rate of change of electron density N.sub.E and negative ion density N_(see, e.g., R. F. Fernsler, A. W. Ali, J. R. Greig and I. M. Vitkovitsky, “The NRL CHMAIR Code: A Disturbed Sea Level Air Chemistry Code,” NRL Memorandum Report 4110 (1979); A. W. Ali, “Electron Energy Loss Rates in Air,” NRL Memorandum Report 5400 (1984); L. G. Christophorou, Atomic and Molecular Radiation Physics (Wiley-Interscience, London, UK, 1971), p. 530; P. Sprangle, J. Peñano, B. Hafizi, D. Gordon and M. Scully, Appl. Phys. Lett. 98, 211102 (2011)) are
∂N.sub.e/∂t=(1+α.sub.rad)Q.sub.rad+S.sub.e−L.sub.e, ∂N_/∂t=S_−L_,
where S.sub.e represents the various electron source terms, L.sub.e is the electron loss terms, S_ represents the ion source, L_ is the ion loss terms (see, e.g., P. Sprangle, B. Hafizi, H. Milchberg, G. Nusinovich and A. Zigler, Physics of Plasmas (to be published, 2013)).

(18) The effect of radioactivity is represented by the first term on the right hand side of the electron rate equation. The steady state electron and negative ion densities are given by
N.sub.e≈(β.sub.nN.sub.n/η)√{square root over ((1+α.sub.rad)Q.sub.rad/β.sub.+)}+(1+α.sub.rad)Q.sub.rad/η≈(β.sub.nN.sub.n/η)√{square root over ((1+α.sub.rad)Q.sub.rad/β.sub.+)}, and N_≈√{square root over ((1+α.sub.rad)Q.sub.rad/β.sub.+)},
where N.sub.n is the neutral air density (N.sub.n˜N.sub.air for low levels of ionization), β.sub.+≈2×10.sup.−6 cm.sup.3/sec is the recombination rate, η≈10.sup.8 sec.sup.−1 is the attachment rate and β.sub.n≈(5−10)×10.sup.−19 cm.sup.3/sec is the negative ion detachment rate due to collisions with neutrals. In the absence of radioactive material (α.sub.rad=0) the ambient (background) electron and negative ion density are N.sub.e10.sup.−3 cm.sup.−3 and N_≈3×10.sup.3 cm.sup.−3, where Q.sub.rad=20 cm.sup.−3 sec.sup.−1. At a distance of 4 m from a radioactive source containing 10 mg of .sup.60Co, the radioactive enhancement factor is α.sub.rad=2.2×10.sup.4 and the elevated electron and negative ion densities are N.sub.e0.2 cm.sup.−3 and N_≈4.7×10.sup.5 cm.sup.−3.

(19) The ionization potential of O.sub.2.sup.− is 0.46 eV and therefore can undergo single-photon photo-detachment with laser radiation of wavelength λ=1 μm (1.24 eV) or λ=0.8 μm (1.55 eV). The photo-detachment rate is v.sub.opt=σ.sub.optcN.sub.ph=σ.sub.optI.sub.o/ω, where c N.sub.ph=I.sub.o/ω is the incident photon flux, I.sub.o is the laser intensity and σ.sub.opt is the photo-detachment cross section. The experimental value for the single-photon ionization cross section of O.sub.2.sup.− is σ.sub.opt(λ=1 μm)≈4.5×10.sup.−19 cm.sup.2 and σ.sub.opt(λ=0.8 μm)≈7.5×10.sup.−19 cm.sup.2. The single-photon ionization rate for O.sub.2.sup.− is therefore,

(20) $v_{\mathrm{opt}}\left[{\sec }^{-1}\right]=I_o\left[W/{\mathrm{cm}}^2\right]\{\begin{array}{cc}2.3,& \lambda =1\mathrm{.Math.m}\\ 3,& \lambda =0.8\mathrm{.Math.m}\end{array},$

(21) The various source and loss terms, in particular the collisional ionization rate, are functions of the electron temperature. The electron temperature is determined by the collisional electron heating (Ohmic heating) by the laser radiation and the cooling effect resulting from excitation of vibrational modes of the air molecules. The equation for the electron temperature T.sub.e (see, e.g., Ya. B. Zel'dovich and Yu. P. Raizer, Physics of Shock Waves and High-Temperature Hydrodynamic Phenomena (Dover, Mineola, N.Y., 2002)),
(3/2)∂(N.sub.eT.sub.e)/∂t=J.Math.E+(3/2)(N.sub.e/τ.sub.cool)(T.sub.e−T.sub.eo)−U.sub.ion∂N.sub.e/∂t,
where T.sub.e is the electron temperature, J.Math.E is the Ohmic heating rate, τ.sub.cool is the electron cooling time due to inelastic collisions, T.sub.eo=0.025 eV is the ambient electron temperature and U.sub.ion is the effective ionization potential of air (˜10 eV).

(22) Frequency Modulation on a Probe Beam

(23) A probe beam propagating through a region of space in which the electron density is changing with time will undergo a frequency change. The electron density in the vicinity of the radioactive source and under the influence of the laser radiation varies in space and in time. Consequently the frequency/wavenumber of an electromagnetic probe beam propagating in such a plasma will vary in space and in time. The one-dimensional wave equation (∂.sup.2/∂z.sup.2−c.sup.−2∂.sup.2/∂t.sup.2)A(z,t)=c.sup.−2ω.sub.p.sup.2(z,t)A(z,t) can be used to determine the frequency/wavenumber shift. Here, ω.sub.p(z,t)=(4πq.sup.2N.sub.e(z,t)/m).sup.1/2 is the plasma frequency and A(z,t) is the vector potential associated with the probe. The vector potential can be expressed in terms of slowly-varying amplitude and phase, in the form A(z,t)=(½)B(z,t)exp[i(k.sub.oz−ω.sub.ot)+θ(z,t)]+c.c., where ω.sub.o is the incident probe frequency and θ(z,t) is the phase. The frequency modulation on the probe beam is given by ω.sub.probe(z,t)=ω.sub.o+Δω(z,t), where Δω(z,t)=ω.sub.po.sup.2/(2ω.sub.o)exp(v.sub.iont)(1−exp(−v.sub.ionz/c)) and v.sub.ion is the ionization frequency. The maximum fractional frequency shift occurs for z>c/v.sub.ion˜1 cm and is Δω.sub.max/ω.sub.o=(ω.sub.po.sup.2/2ω.sub.o.sup.2)exp(v.sub.iont). The effective ionization rate can vary widely but is typically v.sub.ion˜10.sup.11 sec.sup.−1.

(24) Radioactivity Detection Example

(25) The radioactive source is assumed strong enough to produce a radiation enhancement factor of α.sub.rad=10.sup.3. The radiation enhancement factor value is consistent with low quantities of radioactive material as shown in FIG. 2. We take the ionizing laser to have a peak intensity of I.sub.peak=160 GW/cm.sup.2 and pulse duration of τ.sub.laser=1 nsec. In these examples, the probe beam is taken to be a millimeter wave source of frequency f.sub.probe=94 GHz, (λ.sub.probe=3.2 mm). The critical electron density, associated with the probe frequency, ω.sub.probe=ω.sub.p,crit=5.64×10.sup.4n.sub.e,crit.sup.1/2[cm.sup.−3], is n.sub.e,crit=10.sup.14 cm.sup.−3. The background radiation is taken to be Q.sub.rad=30 disintegrations/(cm.sup.3-sec) In the absence of radioactivity, FIG. 3A, the ionizing laser intensity is just below the breakdown level, i.e., the electron density is low, and there is virtually no frequency modulation on the probe beam. FIG. 3B shows the electron density as a function of time in the presence of radioactive material (α.sub.rad=10.sup.3). The electron density at the end of the ionizing laser pulse approaches the value of n.sub.e=10.sup.13 cm.sup.−3 which is an order of magnitude less than the critical electron density.

(26) The frequency modulation on the probe millimeter wave beam is shown in FIG. 4. In the absence of radioactive material there is no frequency modulation on the probe. However, for α.sub.rad=10.sup.3 the fractional frequency modulation is significant and equal to ˜5%, which is readily detectable. The fractional frequency shift on the probe as a function of both axial interaction distance L and time is shown in FIG. 5.

(27) While the present invention has been described with respect to exemplary embodiments thereof, it will be understood by those of ordinary skill in the art that variations and modifications can be effected within the scope and spirit of the invention. For example, the probe laser can comprise a probe millimeter wave source or a microwave source. Also, the presence of an ionizing radioactive material may further be detected by the generation of a spark/air breakdown upon ionization when exposed to the ionizing laser beam. In addition, the invention may be applied to/include spectroscopic signatures from other species present in the atmosphere such as the 337 nm line of nitrogen molecules. Alternatives electromagnetic signatures include i) backscattering and frequency upshifting of radiation from energetic electrons generated by the gamma rays, ii) spectroscopic signature from other molecular constituents in the atmosphere that are excited by the gamma rays and iii) photo-detachment of electrons from atomic oxygen, O.sup.−.