Active remote detection of radioactivity based on electromagnetic signatures

Information

  • Patent Grant
  • 9835761
  • Patent Number
    9,835,761
  • Date Filed
    Friday, January 16, 2015
    9 years ago
  • Date Issued
    Tuesday, December 5, 2017
    6 years ago
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.
Description
FIELD OF THE INVENTION

The invention is directed to the detection of radioactive emission sources, and more particularly to determining a radioactive composition of the source by obtaining an electromagnetic signature associated with the specific type of radioactive emission source.


BACKGROUND OF THE INVENTION

Existing methods for detection of radioactive materials have very limited range (less than a few meters). In addition existing methods are passive resulting in limited sensitivity. The Geiger-Muller tube is the most common type of radioactivity detector. A recently proposed active radioactivity detection concept is based on a high power THz pulse inducing avalanche breakdown and spark formation in the vicinity of the radioactive material, but has very limited stand-off detection range.


It would therefore be desirable to remotely detect radioactive materials with greater precision and specificity.


BRIEF SUMMARY OF THE INVENTION

According to the invention, 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.


A laser-based radioactivity detection concept is a significant advance over existing methods and may have unique advantages depending on the stand-off distance and atmospheric conditions. The invention provides the active remote detection of radioactivity from radioactive sources based on their specific activity (radiation level), enabling stand-off detection at distances of greater than 100 m.





BRIEF DESCRIPTION OF THE DRAWINGS


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



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 60Co according to the invention;



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, Ipeak=160 GW/cm2, τ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



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



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





DETAILED DESCRIPTION OF THE INVENTION
Definition

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.


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 O2 ions. The density of O2 around radioactive material can be several orders of magnitude greater than background levels. The elevated population of O2 extends several meters around the radioactive material. Electrons are easily photo-detached from O2 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.


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


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.


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.


Radiation Enhancement Factor


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γ decreases. A 1 MeV gamma ray has a mean free path in air of Lγ≈130 m.


The ionization rate due to background (ambient) radioactivity is (d Ne/dt)amb=Qrad. At or near ground level, the background ionization rate is typically in the range Qrad˜10-30 pairs/(cm3-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 αrad>>1 and (d Ne/dt)radαradQrad. For a radioactive material of mass Mrad the number of disintegrations per second is vrad=MradArad, where Arad is the specific activity associated with the material. For example, for 60Co, Arad=1.1×103 Ci/g=4.1×1013 disintegration/(g-sec). In the case of 60Co each disintegration results in two gammas of energy Eγ,max=1.173 MeV and Eγ,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 O2 molecules forming O2 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 Ee therefore generates ˜Ee/ΔE low energy electrons. An electron having an energy of 1 MeV has a range in air of 4.6 m.


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








α
rad



c




σ

γ
-
e






N
air



N
γ






E
e




Δ





E




1

Q
rad








v
rad



κ
γ



4

π




L

γ
-
e











E
e




Δ





E




1

Q
rad





exp


(


-
R

/

L
γ


)



R
2




,





is the radiation enhancement factor, custom characterEecustom character is the average electron energy, custom characterσγ-ecustom character is the effective average cross section for electron generation by gammas, i.e., Compton absorption and photoelectric processes, custom characterLγ-ecustom character=(custom characterσγ-ecustom characterNair)−1 is the average mean free path for electron generation by gammas and Nair=2.7×1019 cm−3 is the air density at STP. In the absence of radioactive material αrad=0. In FIG. 2 the radiation enhancement factor α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 60Co and indicates that the enhanced level is significant for ranges extending up to several meters. As an example, for custom characterEecustom character=0.5 MeV, Mrad=10 mg, vrad=MradArad=8.2×1011 disintegrations/sec, κγ=2, R=50 cm, custom characterLγ-ecustom character=100 m and Qrad=20 disintegrations/(cm3-sec) the radiation enhancement factor is αrad≈2×106 which is far above the background level.


Electron and Ion Density Evolution (Air Chemistry)


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 NE 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

Ne/∂t=(1+αrad)Qrad+Se−Le, ∂N_/∂t=S_−L_,

where Se represents the various electron source terms, Le 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)).


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

Ne≈(βnNn/η)√{square root over ((1+αrad)Qrad+)}+(1+αrad)Qrad/η≈(βnNn/η)√{square root over ((1+αrad)Qrad+)}, and N_≈√{square root over ((1+αrad)Qrad+)},

where Nn is the neutral air density (Nn˜Nair for low levels of ionization), β+≈2×10−6 cm3/sec is the recombination rate, η≈108 sec−1 is the attachment rate and βn≈(5−10)×10−19 cm3/sec is the negative ion detachment rate due to collisions with neutrals. In the absence of radioactive material (αrad=0) the ambient (background) electron and negative ion density are Ne10−3 cm−3 and N_≈3×103 cm−3, where Qrad=20 cm−3 sec−1. At a distance of 4 m from a radioactive source containing 10 mg of 60Co, the radioactive enhancement factor is αrad=2.2×104 and the elevated electron and negative ion densities are Ne0.2 cm−3 and N_≈4.7×105 cm−3.


The ionization potential of O2 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 voptoptcNphoptIo/custom characterω, where c Nph=Io/custom characterω is the incident photon flux, Io is the laser intensity and σopt is the photo-detachment cross section. The experimental value for the single-photon ionization cross section of O2 is σopt(λ=1 μm)≈4.5×10−19 cm2 and σopt(λ=0.8 μm)≈7.5×10−19 cm2. The single-photon ionization rate for O2 is therefore,








v
opt



[

sec

-
1


]


=



I
o



[

W
/

cm
2


]




{





2.3
,




λ
=

1





µm







3
,




λ
=

0.8





µm





,







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 Te (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)∂(NeTe)/∂t=custom characterJ·Ecustom character+(3/2)(Necool)(Te−Teo)−Uion∂Ne/∂t,

where Te is the electron temperature, custom characterJ·Ecustom character is the Ohmic heating rate, τcool is the electron cooling time due to inelastic collisions, Teo=0.025 eV is the ambient electron temperature and Uion is the effective ionization potential of air (˜10 eV).


Frequency Modulation on a Probe Beam


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 (∂2/∂z2−c−22/∂t2)A(z,t)=c−2ωp2(z,t)A(z,t) can be used to determine the frequency/wavenumber shift. Here, ωp(z,t)=(4πq2Ne(z,t)/m)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(koz−ωot)+θ(z,t)]+c.c., where ωo is the incident probe frequency and θ(z,t) is the phase. The frequency modulation on the probe beam is given by ωprobe(z,t)=ωo+Δω(z,t), where Δω(z,t)=ωpo2/(2ωo)exp(viont)(1−exp(−vionz/c)) and vion is the ionization frequency. The maximum fractional frequency shift occurs for z>c/vion˜1 cm and is Δωmaxo=(ωpo2/2ωo2)exp(viont). The effective ionization rate can vary widely but is typically vion˜1011 sec−1.


Radioactivity Detection Example


The radioactive source is assumed strong enough to produce a radiation enhancement factor of αrad=103. 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 Ipeak=160 GW/cm2 and pulse duration of τlaser=1 nsec. In these examples, the probe beam is taken to be a millimeter wave source of frequency fprobe=94 GHz, (λprobe=3.2 mm). The critical electron density, associated with the probe frequency, ωprobep,crit=5.64×104ne,crit1/2[cm−3], is ne,crit=1014 cm−3. The background radiation is taken to be Qrad=30 disintegrations/(cm3-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 (αrad=103). The electron density at the end of the ionizing laser pulse approaches the value of ne=1013 cm−3 which is an order of magnitude less than the critical electron density.


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 αrad=103 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.


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.

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; andidentifying, 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 Ipeak=160 GW/cm2 and a pulse duration of τ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; anda 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, andidentifying, 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 Ipeak=160 GW/cm2 and a pulse duration of τ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 O2− 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 O2− 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; anda 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, andidentifying, 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 Ipeak=160 GW/cm2 and a pulse duration of τlaser=1 sec.
CROSS-REFERENCE TO RELATED APPLICATIONS

This Application claims the benefit of U.S. Provisional Application 61/935,903 filed on Feb. 5, 2014, incorporated herein by reference.

US Referenced Citations (4)
Number Name Date Kind
20020047543 Sugai Apr 2002 A1
20090001889 Miles Jan 2009 A1
20120170599 Sprangle Jul 2012 A1
20120305773 Wu Dec 2012 A1
Non-Patent Literature Citations (8)
Entry
P. Sprangle, J. R. Peñano, and B. Hafizi, Propagation of intense short laser pulses in the atmosphere, Phys. Rev. E 66, 046418 (2002).
P. Sprangle, J. Peñano, B. Hafizi, D. Gordon and M. Scully, Remotely induced atmospheric lasing, Appl. Phys. Lett. 98, 211102 (2011).
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).
V.L. Granatstein and G.S. Nusinovich, Detecting excess ionizing radiation by electromagnetic breakdown of air,J. Appl. Phys. 108, 063304 (2010).
G.S. Nusinovich, P. Sprangle, C.R. Talamas and V.L. Granatstein, Range, resolution and power of THz systems for remote detection of concealed radioactive materials,J. Appl. Phys. 109, 083303 (2011).
G.S. Nusinovich, R. Pu, T.M. Antonsen, O.V. Sinitsyn, J. Rodgers, A. Mohamed, J. Silverman, M. Al-Sheikhly, Y.S. Dimant, G.M. Milikh, M.Yu, Glyavin, A.G. Luchinin, E.A. Kopelovich and V.L. Granatstein, Development of THz-range Gyrotrons for Detection of Concealed Radioactive Materials,J. Infrared Milli. Terahz. Waves, 32 380 (2011).
Y.S. Dimant, G.S. Nusinovich, P. Sprangle, J. Peñano, C.A. Romero-Talamas and V.L. Granatstein, Propagation of gamma rays and production of free electrons in air, J. Appl. Phys. 112, 083303 (2012).
Related Publications (1)
Number Date Country
20160377761 A1 Dec 2016 US
Provisional Applications (1)
Number Date Country
61935903 Feb 2014 US