The invention relates to solid-state high energy radiation detectors, more specifically to scintillating detectors made of direct-gap semiconductors.
This patent application builds upon an earlier U.S. patent application Ser. No. 11/144,443, filed Jun. 6, 2005 and titled “Semiconductor scintillation high-energy radiation detector”. In that application, a semiconductor scintillator detector is disclosed in which ionizing radiation creates electron-hole pairs that subsequently undergo interband recombination, producing infrared light to be registered by a photo-detector. To make the semiconductor essentially transparent to its own interband radiation, a direct gap semiconductor, such as InP or GaAs, is heavily doped with shallow donors to produce the blue shift of the absorption edge, an effect known as the Burstein shift. Owing to the high Fermi level EF for electrons, in such materials the absorption edge is shifted to shorter wavelengths, compared to the interband emission spectrum. The absorption of interband radiation is suppressed by the absence of vacant electronic states under the Fermi level and the suppression factor depends exponentially on the EF/kT ratio.
The heavy doping makes the semiconductor more transparent to the scintillation photons and enable to extract said photons from deep inside the semiconductor. Limitations on the absorption length are discussed further below. In addition, the heavy doping shortens the radiative recombination time τ of holes, according to τ=1/BND, where ND is the donor concentration and B=10−10 cm3/s is the radiative recombination coefficient. For the doping level of 1019 cm−3 in InP one has τ˜10−9 s. Typically, the non-radiative time is about 10−7 s. Therefore, the device experiences little loss through the non-radiative channel of recombination and provides a fast response time of 1 ns.
The calculated absorption length λ in InP at the optical wavelength of 0.92 μm, typical for the InP interband emission spectrum is shown in
One of the main objectives of the present invention is to lift said thickness limitation by integrating multiple semiconductor slabs into a single block, each of said slabs provided with its own epitaxial photo-detector for the registration of scintillator signal, while the slab thickness is kept in the range of about 0.5 mm, to minimize the light absorption and thus maximize the light delivery to the photo-detector.
Another objective of the present invention is to disclose a new method of registration of the radiation events using registration of the radiation penetration into the depth of the detector which provides a new method of analysis of the radiation energy.
Still another main objective of the present invention is to disclose new analysis methods that become available with 2-dimensional (2D) pixellation of the photo-detector layer in each slab, thereby producing a possibility of a 3-dimensional (3D) pixellation of registered ionization events.
Finally, still another main objective of the present invention is to disclose an analysis methodology of the inventive radiation detector that enables interpreting the measured 3D-pixellated data for the determination of both the radiation energy and direction to the radiation source.
To circumvent the problem of thickness limitation, arising from the finite absorption length for the scintillator light in the semiconductor slab, a new design is disclosed, according to the present invention, wherein the detector represents a multiple stack of relatively thin semiconductor slabs, each slab being thinner than the scintillator light absorption length. Each slab is further endowed with an epitaxially grown photo-detector that has high absorption coefficient for scintillator light and a sufficient thickness to fully absorb said scintillator light. The stack of such slabs according to the present invention can accommodate a long absorption length of the high-energy gamma radiation without any loss in either the scintillator yield or the speed of response. In the preferred embodiment, the epitaxial photodetector is implemented as a p-i-n junction made of direct-gap semiconductor material with the absorption coefficient of about 104 cm−1 and of sufficient thickness that the light-to-current conversion efficiency approaches 100%. This provides an advantage over photo-multipliers used in conventional scintillators, where said efficiency is typically only 25%, as discussed in the well-known textbook by G. F. Knoll, “Radiation detection and measurement”, 3rd edition, Wiley, N.Y. (2000). Another advantage of the inventive gamma-radiation scintillator over the photo-multiplier based scintillators is owing to its achievable miniature dimensions, where the total volume can be as low as ˜1 cm−3.
The epitaxially grown photodetector material is chosen to possess a substantially the same refractive index as the bulk material of the slab, thus minimizing the optical losses from internal light reflection at the epitaxial interface between the slab and the photodetector. This feature serves to maximize the collection of scintillator light.
For InP-based detector slab, the preferred lattice-matched material is the quaternary compound InGaAsP, which can be employed to grow the photo-detector structure with an energy gap varied widely between 1.35 eV and 0.8 eV. For GaAs, the preferred epitaxial material for the photo-detector is a dilute-nitride InGaAs/N stress-compensated structure that allows growth of a photo-sensitive layer with the energy gap below that of GaAs.
In the preferred embodiment, the epitaxial photo-diode is implemented as a p-i-n junction. As is well known to those skilled in the art, the photosensitive intrinsic “i” layer must be fully depleted. In the inventive design the total thickness of the photo-sensitive layer to be depleted is only ˜2 μm, which can be easily depleted even when the background dopant density level is relatively high. For example, for background donor concentration of 1015 cm−3, one needs only about 3V of reverse biasing for full depletion.
The simplest, non-pixellated detector embodiment, according to the present invention, comprises a stack of semiconductor slabs with a parallel connection of all epitaxial p-i-n junctions in the stack. This results in an integrated radiation detector response from its full volume. This electrical response is then amplified by a single amplifier circuit, located outside of the radiation detector. It is preferable to use Si technology for electronic signal processing and attach a Si chip. This chip may include an analog-to-digital conversion circuitry that delivers amplified digital information to the recording device or a computer.
Another preferred embodiment employs one amplifier per each slab. The amplified signal from each slab is converted to a digital form and then delivered it to the signal analyzing system. One advantage of this embodiment is that the signal to be delivered is no longer a one-nanosecond analog pulse, but rather a train of standard digital pulses encoding the available information. The information available includes the index of the active slab, i.e. that slab where the signal has fired and the time when the signal occurred. The index information relates to the z-coordinate of the ionization characterization, where the z-axis is along the direction of slab stacking. This design enables tracing of the one-dimensional signal dependence on the depth within the stack. The resultant one-dimensional signal profile can be analyzed to provide an additional characterization of the energy of ionizing radiation according to its penetration depth dependence.
Still another preferred embodiment employs a two-dimensional photo-detector array on the surface of each slab in the stack. This provides a 2D pixellation that enables recording the registration of xy coordinates of an ionization event in each slab. The entire stack of detector slabs can then be viewed as a three-dimensional integration of 2D-pixellated detector arrays. This provides a 3D pixellation of the ionization events.
The contemplated typical thickness of the single scintillator slab is in the range between 0.25 and 2.5 mm, as dictated by the absorption length data, such as in
The lateral pixellation of the photosensitive layer in each slab can be arranged in a number of ways, known to those skilled in the art of photo-diode matrices. In the preferred embodiment, the lateral pixellation is accomplished by providing a set of horizontal and vertical stripes of contact layers of n+ and p+ polarity, thereby forming a pixel matrix wherein every pixel is determined by the intersection of one vertical and one horizontal line connected to the pixel p-i-n junction. All horizontal lines and all vertical lines in every slab are connected to their respective common electrodes in order to enable the application of a reverse-biased voltage to the entire pixel array. To provide pixel addressing with electronic amplification of the signals, a row and a column of amplifiers are introduced in such a way that every pixel is uniquely identified by the action of two amplifiers. Such a pixellated matrix architecture allows collection of the 3D information per each ionization event.
The information obtained in 3D-pixellated detector structure enables the implementation of novel analysis techniques, based on the kinematics of the interaction of gamma radiation with the pixel material. Interacting incident gamma quanta experience a series of Compton scattering events and ultimately, when the gamma energy becomes sufficiently low, a photoelectric absorption. Each interaction event results in a partial transfer of the gamma quantum energy to form electron and hole pairs in a particular pixel. These e-h pairs recombine to produce the scintillating photons that in the inventive device are registered pixel-wise. Typically, one expects 3 to 6 pixels to simultaneously become active upon the absorption of a single gamma quantum of exemplary energy 0.662 MeV. The measured information includes the 3D position coordinates (x, y, and z) of the active pixels and the total number (N) of scintillating photons produced in each active pixel. The number N serves as a measure of the energy deposited in said pixel.
The detector output lines 23 and 24 enter a Si electronic chip 26 to amplify the integrated signal and convert it into the digital form prior to delivering it through the line 27 to recording, analyzing and storing system outside the integrated block. The chip is attached to the substrate 28 (glass) which supports the entire device.
As discussed above, the chosen semiconductor material is optically direct to provide high efficiency of scintillating light radiation. It is also heavily doped to minimize the light losses in the slab, as well as to reduce the non-radiative component of the electron-hole recombination and decrease the emission response time. Two best candidates for the slab materials are n+ InP and n+ GaAs, since they are the most technologically advanced and mass produced for the optoelectronic industry.
The preferred detector embodiment uses standard InP or GaAs wafers as semiconductor slabs. Besides their lower cost, such an approach allows one to choose the wafer thickness below the penetration depth of ionizing radiation, while the combined thickness of all wafers in the integrated block can be made thicker than the absorption length of the ionizing radiation. For example, 0.3 mm-thick wafers of InP or GaAs are commercially available and are thinner than the characteristic absorption length for the γ-radiation energy above ˜80 KeV, see the above cited textbook by G. F. Knoll. Photons of lower energy, X-rays, will be stopped primarily in the first layer.
The top area of each slab in
Those skilled in the art will recognize that owing to the random nature of scintillating radiation approximately 50% of the scintillating light intensity will go to one surface of the slab and 50% to the other. In order to optimize the light collection efficiency detector layers can be deposited on both sides of the slab. In the preferred embodiment, however, detector layers are deposited on one side only, the other side being coated with a light reflecting mirror structure (not shown) to reflect light at the surface. Such mirror structures are well-known to those skilled in the art; exemplarily, they can be made of a thin metal film, such as Al, or thin dielectric/metal layers, such as SiO2/Al, of about 0.1 μm combined thickness.
The 2μ-thick undoped InGaAsP sensitive layer grown on n+InP wafer is followed by a ˜0.1 μm-thin p+InP contact layer. The contact lines 23 and 24 in
This embodiment enables one to trace the one-dimensional (z-coordinate) dependence of the scintillating response on the depth within the stack. This design offers the valuable possibility to measure the depth profile of radiation penetration and thus also estimate the incident radiation energy. The resolution of this method depends on the wafer thickness: smaller thickness of each sub-detector results in more precise energy measurement.
b shows the side view of the detector. Notches 42 in each p-i-n junction layer indicate etching of this layer into three columns, while 3 pixels in each column are connected with metal lines 43. Only right side of the Si chip 44 is shown. Each chip is attached to its wafer using the support glass plate 46 on top of which three metal lines 45 are deposited to connect the rows 40 to the Si chip 44. Line 47 serves as a signal output from every chip.
In this design, each pixel circuit formed at an intersection of a single horizontal and a single vertical line comprises two amplifiers. This architecture is further illustrated in the circuit diagram of
Every pixel is individually addressed, and signal registration from the electronic chips will identify the signal amplitude, the timing and the location of the ionization event within both the pixel matrix area and the detector depth. Thus, such a design allows 3-dimensional signal registration. Such a 3D detector can perform all the functions of the discussed above non-pixellated versions shown in
Finally, the same pixellated structure can be made when every wafer, shown in
a shows the preferred circuit architecture for the pixellated detector. Each pixel within a single slab is uniquely identified by the two amplifiers that carry the electrical current generated by a given interaction event. The firing pixel and the two associated amplifiers are further illustrated in
An important advantage of a 3D pixellated radiation detector is that it offers the valuable possibility for data analysis based on simultaneous signal registration by several pixels in the stack. This enables a direct measurement of the incident particle energy, said measurement being complementary to the conventional statistical spectroscopy and free from complications associated with Compton escape processes. This goes to the heart of the contemporary homeland security needs, where accurate spectral characterization of detected gamma radiation is of the essence to avoid “false alarms”. The proposed technology offers unprecedented fidelity in isotope discrimination.
cos θi=1+Ei-1−1−Ei−1 (1a)
L
i
=E
i-1
−E
i (1b)
where Ei-1 and Ei are, respectively, the gamma-particle energies before and after the i-th interaction that deposits the energy Li in the semiconductor. All energies in equation (1) above and equation (2) below are measured in the units of electron rest mass me≈511 KeV. The angle θi is the angle of scattering in the i-th interaction, illustrated in
Similarly, the angle θ1 is determined in terms of L1, L2 and θ2. The knowledge of θ1 does not yet determine the direction {circumflex over (n)}0 of the incident quantum, which is placed on the cone of angle θ1 about the known direction {circumflex over (n)}1. Intersections of such cones for two or more incident gamma quanta coming from the same point source, uniquely determines the direction to the source.
The described method of analysis based on equations (1), often referred in the literature as “Compton telescope”, is well-known, see for example, S. E. Boggs and P. Jean, “Event reconstruction in high resolution Compton telescopes”, Astron. Astrophys. Suppl. Ser. 145, 311{321 (2000). It is usually implemented by using an assembly of discrete Ge or Si diode detectors. A large number of such detectors is necessary, for otherwise the events where three detectors fire at the same time would be extremely rare. The large number of discrete detectors leads to a bulky size and considerable cost of such detector assemblies. An advantage of the inventive integrated detector is its relatively small size and low cost combined with the possibility of a very high density of pixels that guarantees a high rate of three-pixel count. Also the high speed of response of the inventive detector enables temporal separation of different incoming particles even when their rate of incidence is itself high. When fully optimized, the inventive detector is expected to distinguish between quanta arriving only 10 nanoseconds apart.
A practical use of the Compton telescope analysis method requires that for a given incident particle or gamma quantum the order of interactions is correctly identified as first, second and third. It should be noted that all three (and possibly more than three) interactions occur essentially at the same time, since the inventive detector cannot resolve them in time, which would require picosecond resolution. No known at this time detection principle can be deployed to order the interaction according to their time of occurrence. However, those skilled in the art know that interactions can be ordered with a reasonable confidence according to their intensity and the cluster geometry. This ordering is enabled by the directionality and the energy-transfer of Compton scattering, as embodied in the well-known Klein-Nishina formula, describing the anisotropic scattering cross-section, σ(θi), viz.
The ordering procedure is often referred to as tracking algorithm or event reconstruction. Current tracking algorithms claim up to 70% success in correctly ordering events produced by a 1 MeV photon.
The anisotropy of the Compton scattering cross-section can also be employed for a different type of analysis that is advantageous for the determination of direction to source and is enabled by the inventive detector. The new analysis method is illustrated in
and then draw a vector {right arrow over (ρ)}1 from {right arrow over (C)} to the first interaction.
Because of the directional (axial) symmetry of Compton interaction, the average vector {right arrow over (ρ)}1 averaged over a statistical ensemble comprising multiple realizations of the cluster, corresponding to different incoming photons, tends to a mean value {right arrow over (ρ)} which is parallel to {circumflex over (n)}0 and points in the direction to the source. It is easy to write an expression for the value of {right arrow over (ρ)}1 averaged over a sub-ensemble of Nn clusters, that comprise exactly n interactions,
where {right arrow over (δ)} is an error vector normal to {right arrow over (ρ)} and pn is the probability of correctly guessing the first interaction in a cluster of n points. So long as the guess is much better than random, pn>1/n, the procedure converges rapidly, as has been verified by a Monte Carlo simulation. We are free to analyze any sub-ensemble separately, the full ensemble comprising N=N1+N2+ . . . +Nn sub-ensembles. The final {right arrow over (ρ)}1 is a statistical average of the values {right arrow over (ρ)}1 obtained for the partial sub-ensembles.
Both the polar and azimuth angles (θ,φ) of {circumflex over (n)}0 are determined. As we examine more and more incoming γ particles, the precision improves. For N=1,000 we find {circumflex over (n)}0 to within about 2°, for N=104 to about 0.5°. The relative error associated with the finite pixel size d goes as (d/D)2, where D=D(E0) is the typical linear dimension of a cluster. For E0=660 keV and d≦1 mm, this error does not affect the above estimates of the precision. Neither does the error associated with the Doppler effect, which only makes a small contribution to the {right arrow over (δ)} in Eq. (4).