In one embodiment, an apparatus comprises: a detection pixel including a first material that is responsive to incident photons to produce Compton electrons and an electric or magnetic response in the first material, wherein the incident photons and the Compton electrons each have a characteristic travel distance in the first material; and wherein the detection pixel has at least one dimension that is selected according to at least one of the characteristic travel distance of the incident photons and the characteristic travel distance of the Compton electrons.
In another embodiment, a system comprises an array of detection pixels, each detection pixel in the array of detection pixels including a first material that is responsive to incident photons to produce Compton electrons and an electric or magnetic response in the first material, wherein the incident photons and the Compton electrons each have a characteristic travel distance in the first material; and wherein each detection pixel in the array of detection pixels has at least one dimension that is selected according to at least one of the characteristic travel distance of the incident photons and the characteristic travel distance of the Compton electrons.
In another embodiment, an x-ray imaging system comprises: a source of x-rays; an array of detection pixels, each detection pixel in the array of detection pixels including a first material that is responsive to the x-rays to produce Compton electrons; and detection circuitry operably connected to the array of detection pixels and configured to measure a change in the first material produced by the Compton electrons.
In another embodiment, a method corresponding to a detection pixel including a first material having an electric or magnetic response to electromagnetic energy comprises: receiving information corresponding to a resonance of a detection circuit that includes the detection pixel, wherein the resonance is characterized by a resonant frequency and a linewidth; and determining a change in the resonant frequency of the detection circuit to determine the presence of electromagnetic energy in the detection pixel.
The foregoing summary is illustrative only and is not intended to be in any way limiting. In addition to the illustrative aspects, embodiments, and features described above, further aspects, embodiments, and features will become apparent by reference to the drawings and the following detailed description.
The patent or application file contains at least one drawing executed in color. Copies of this patent or patent application publication with color drawing(s) will be provided by the Office upon request and payment of the necessary fee.
In the following detailed description, reference is made to the accompanying drawings, which form a part hereof. In the drawings, similar symbols typically identify similar components, unless context dictates otherwise. The illustrative embodiments described in the detailed description, drawings, and claims are not meant to be limiting. Other embodiments may be utilized, and other changes may be made, without departing from the spirit or scope of the subject matter presented here.
In some embodiments the material of the pixel 102 includes a ferroelectric material, a piezoelectric material, or a multiferroic material. These materials may have values of electronic and magnetic permeability that change when exposed to incident photons in the energy range of interest. Examples of ferroelectric materials include, but are not limited to, titanates of lead, bismuth, barium, zirconium, and their alloys (a prototype being lead-zirconate-titanate, i.e., PZT) which would be appropriate for the higher energy range of photons. An example of a material that exhibits both ferroelectricity and piezoelectricity is a non-linear optical material such as lithium niobate, which would be appropriate for the lower energy range of the photons, and where the detection may be through optical methods. Examples of multiferroic materials include perovskites such as bismuth manganite, bismuth ferrite, and the related family of transition metal and rare earth manganites and ferrites. Other multiferroic materials include the non-perovskite and non-oxide compounds such as spinel compounds.
In one exemplary embodiment, the pixel 102 includes a ferroelectric material, and the dimensions and material of the pixel 102 are selected according to the following considerations. The detector can be configured such that the volume of interaction is comparable to the volume of an individual pixel, to yield a measurable response. The material is chosen so that the Compton scattering process dominates over the photoelectric scattering process for the incident high energy photons used in the specific application. The incident photons and the Compton electrons each have a characteristic travel distance in the material, which is dependent on the energy of the incident photons and the average atomic density of the detector material. In the scattering process, the incident photon travels a distance L1 before a scattering event, leading to a Compton electron. As the fraction of energy that is transferred to the Compton electron is large, the Compton electrons 108 will then have a high kinetic energy and a travel length L2 that is significant. As these electrons 108 travel through the material, they disturb the lattice structure sufficiently to reduce the intrinsic dielectric (and/or magnetic permeability) constant and yield an effective volume of interaction. Detailed calculations demonstrating these concepts are included in Appendix A, which is attached hereto. Calculations show that optimal pixel size in the range of L1 to L2, however actual selection of pixel size will depend on considerations specific to an application and may vary from this optimal pixel size.
In some embodiments, such as applications involving higher energy photons, the pixel 102 includes a material with relatively large atomic numbers to provide relatively higher stopping power. The higher atomic density results in relatively shorter lengths L1 and L2. In the calculations, for a typical material such as PZT, L1 is approximately 10 micrometers and L2 is approximately 2 micrometers, defining the range of dimensions for the pixel. In another example such as lithium niobate, the lower atomic number would make the material appropriate for lower energy photons, and yield larger travel lengths L1 and L2 while maintain a sufficient ratio. For each case, the detector pixel dimensions would be chosen accordingly.
Although the pixel 102 shown in
In other detection configurations, such as embodiments in which the pixel material is a multiferroic, the pixel may be configured as an inductor in a resonance circuit, operating on the same principals as a ferroelectric capacitor in a resonance circuit.
In some embodiments the array 400 of pixels 102 forms an imaging system. The dimensions of the pixels 102 may be selected according to the considerations as previously described herein. The pixels are each operably connected to the detection circuitry 402, and the detection circuitry may include a multichannel analyzer configured to identify which pixels have received photons.
In some embodiments the pixel materials in the array are not all the same. For example, pixels having different materials may be incorporated in an array where the materials have different ranges of electromagnetic energy they are capable of detecting, thus providing for a detector with a broader range of detection capability than a detector whose pixels are all the same.
The system includes detection circuitry (not shown) as previously described herein, where the detection circuitry is operably connected to the array of detection pixels and configured to measure a change in the first material produced by the Compton electrons, again as previously described herein. In some embodiments the detection circuitry is configured to compute a tomographic image of a subject, i.e., a CT scan. The source of x-rays and the array of detection pixels may be configured on a stage or other type of mount that is movable relative to a subject in order to scan the subject. The system may in some embodiments be configured to determine one or more properties corresponding to the incoming photons, such as an energy and/or an energy range, and/or an intensity.
Disclosed below are calculations corresponding to the high efficiency photon detection approach described herein. It is comprised of four sections describing: 1) dominance of the Compton process in comparison to the photoelectric process for the appropriate choice of material and x-ray energies, 2) demonstration of significant distortion of a prototype lattice in the presence of Compton electrons, 3) calculation of the number and travel length of the electrons for a specified x-ray energy to determine the volume of interaction and thus the size of the detector pixel, and 4) determination of single-photon sensitivity of the detector pixel.
1. Calculation of the Compton Cross-Section
High energy photons interact with matter through several possible paths, determined by the energy of the photons and the atomic number of the constituent materials in which they travel. Classically, the sum of all these paths results in an overall intensity that drops exponentially with travel length as the x-rays are transformed into other forms of energy. The length scale over which the intensity decreases is classically characterized by a mass attenuation parameter. The PbTixZr1-xO3 (PZT) compound mass attenuation is calculated according to
where
and wi are the mass attenuation constants of the constituent atoms and the corresponding atomic weight fractions, respectively. The resulting effective attenuation length calculated from the tabulated data is shown in
X-ray imaging for security and medicine have energies ranging from tens to hundreds of keV. In this energy range, there are two possible interaction mechanisms: Compton and photoelectric scattering. Compton scattering results in an ejected outer shell electron and a significant portion of the x-ray energy is transferred into the kinetic energy of the scattered electron. In contrast, the photoelectric process involves scattering of inner-shell electrons, inducing an inner shell electronic transition where part of the energy is transferred to secondary x-ray photons. As the travel length of electrons is significantly shorter than photons in matter, the Compton process results in a larger quantity of energy being dissipated in a smaller volume. This causes sufficient perturbation of the lattice structure to alter its intrinsic dielectric value for a piezoelectric or ferroelectric material (and permeability in the case of a multiferroic). The result is a dramatic reduction in the large dielectric constant (permeability), yielding a measurable electronic response due to a change in the electronic impedance of the detector circuit.
The detector material should be of the appropriate atomic density so that the Compton process dominates the scattering at the operating energy of the x-ray source. The Compton process dominates at higher x-ray energies and lower (Z) detector density. Higher energy applications, such as computed tomography and security screening, would use higher (averaged) atomic number detector materials such as PZT, whereas lower energy applications would use less dense materials such as LiNbO3. Calculation of the Compton cross-section also yields the average number of Compton electrons produced in the detector by a single photon, and thus the magnitude of the response (change in capacitance or inductance).
The relative probabilities of Compton scattering and the photoelectric effect can be directly assessed from the ratio between the corresponding scattering cross-sections as functions of the incident photon energy and the detector (Z) density. For Compton scattering, the well-established Klein-Nishina formalism is used:
where r0=2.817×10−15 m is the classical electron radius and
The value of k determines whether scattering occurs within the classical or the relativistic frame of reference associated with the moving electron. For the highest incoming photon energies considered below (150 keV), k=0.29, and thus relativistic effects do not play a significant role in these considerations. The effective scattering cross-section of the photoelectric effect in the non-relativistic treatment is:
where α=1/137 is the fine structure constant. In
is plotted for the constituent atoms, as well as for PZT as a whole, as a function of incident photon energy.
The data in
2. Demonstration of Large Lattice Distortion in the Presence of Compton Electrons
The dielectric constant modification by local charging, as well as by the lattice distortions caused by electric charging may be obtained by quasi ab initio density functional theory (DFT) calculations. The numerical calculation presented here is a concurrent cell and structure optimization performed at an externally applied isotropic pressure of 1.0 atm. The simulated system corresponds to x=1 (PbTiO3). The simulation is implemented with use of the Perdew-Burke-Ernzerhof (PBE) exchange correlation function with the appropriately selected basis functions. The simulation utilizes the Gaussian plane wave (GPW) approach to representing the electron wave-functions within the CP2K simulator. The basic structure of the calculation is the energy-optimized electrically neutral structure shown in
Positive charge: The electrically charged system (total charge +1e, or +e/8 per unit cell due to a leaving Compton electron, corresponding to the inner cone in
Negative charge: Extreme anisotropic strain is observed along the X-axis. The resulting lattice periodicity is: X: 5.06 Å (27% stretching), Y: 3.86 Å (3.26% compression), Z: 3.86 Å (3.26% compression).
The anisotropy of strain indicates a significant modification of the dielectric constant from direct theoretical calculations at high charge densities. The strain values can be used for an empirically based quantitative estimation of the dielectric constant modification. Existing experimental data or the measured PZT voltage/field-to-strain constants can be used for estimation.
3. Calculation of Travel Length and Volume of Interaction
As the x-ray fluence within the material decreases exponentially with the penetration distance, the probability of Compton scattering quickly decreases with decreasing incident photon energy as the single photon continues to propagate within the detector. Thus Compton scattering only dominates the earlier stages of scattering for the selected material leading eventually to a photoelectric absorption resulting in complete absorption of the photon. To quantitatively illustrate the discussion above, a Monte-Carlo simulation of the scattering process is implemented. The Klein-Nishina description of each scattering event is used:
where λi and λf are the photon wavelengths before and after a given scattering event, respectively; m, Eb, and θ are the electron's rest mass, the atomic binding energy, and the photon's scattering angle relative to its initial direction, respectively. The prefactor term does not play a significant role during the relatively early stages of the scattering process (Eb<<hc/λi), which is where Compton scattering dominates, as discussed above.
The stochastic portion of the Monte Carlo simulation is set by the random scattering angle φ of the recoiling Compton electron (relative to the photon's direction prior to each scattering event). The probability distribution between φ=0 and φ=2π can be assumed to be uniform. From conservation of momentum, at every scattering event
The results are shown in
4. Calculation of Response to a Single X-Ray Photon
The kinetic energy of the high-energy electrons is eventually converted to lattice vibrations. The highly transient effect of the initial lattice shockwaves generated by the electrons is of interest, as it contributes to disrupting the local polarizability of the material.
The electrons undergoing collisional scattering in the crystal will produce spherical lattice shockwaves, effectively forming a conical wave-front along their paths within the detector material. Because the electrons' velocities (and energies) decrease exponentially in time, the “final” travel distance of the electron is asymptotic. The time it takes the electron to reach the asymptote is calculated to be 10 ps, which is the effective onset time for disruption. The speed of sound in PZT is vs=5000 m/s, and yields a radius of the wave-front cone of 50 nm. With a travel length of a few micrometers, the effective total volume of the resulting highly perturbed lattice can be calculated.
Because the orientation of the wave-fronts is spherically symmetric, a further simplification is possible as long as the effective total volume of the perturbed region is preserved (right pane of
Assuming the geometry and the rectangular mesh as shown in the right pane of
where N is the total number of perturbed regions (equal to the number of Compton electrons) and A and d are the cross-section and the height of the detector, respectively. Along the horizontal (x) axis, within each “slice” the mesh elements are effectively connected in parallel, while the “slices” are in series. Within each slice, the effective capacitance is therefore
where Δx and ΔA=Δx2 are the element height and the cross-sectional area of each rectangular element, respectively; ε1 and ε2 are the effective dielectric constants of bulk PZT and the highly perturbed region, respectively. The sums are over the elements contained within each slice. With use of Eq. (6), Eq. (7) yields
The overall effective capacitance is obtained from all slices
resulting in a simple expression:
where Ω is the effective volume of all perturbed regions and Ωt is the volume of the entire detector. The result is successfully reduced to phase-separate volumetric mixing. Here ε2<<ε1, as the transient volume distortion yields a large reduction in the dielectric constant and Eq. (10) reduces to
Using the above values, each single photon will result in
or 0.1%. As the intrinsic dielectric constant of the material is high, the change is readily measurable through electronics based on resonance circuitry, or other forms.
While various aspects and embodiments have been disclosed herein, other aspects and embodiments will be apparent to those skilled in the art. The various aspects and embodiments disclosed herein are for purposes of illustration and are not intended to be limiting, with the true scope and spirit being indicated by the following claims.
If an Application Data Sheet (ADS) has been filed on the filing date of this application, it is incorporated by reference herein. Any applications claimed on the ADS for priority under 35 U.S.C. §§119, 120, 121, or 365(c), and any and all parent, grandparent, great-grandparent, etc. applications of such applications, are also incorporated by reference, including any priority claims made in those applications and any material incorporated by reference, to the extent such subject matter is not inconsistent herewith. The present application claims the benefit of the earliest available effective filing date(s) from the following listed application(s) (the “Priority Applications”), if any, listed below (e.g., claims earliest available priority dates for other than provisional patent applications or claims benefits under 35 USC §119(e) for provisional patent applications, for any and all parent, grandparent, great-grandparent, etc. applications of the Priority Application(s)). This application claims priority under 35 U.S.C. §119(e) to Provisional Patent App. No. 62/267,677, filed on Dec. 15, 2015, titled “High Efficiency Photon Detection,” which application is hereby incorporated by reference in its entirety. If the listings of applications provided above are inconsistent with the listings provided via an ADS, it is the intent of the Applicant to claim priority to each application that appears in the Domestic Benefit/National Stage Information section of the ADS and to each application that appears in the Priority Applications section of this application. All subject matter of the Priority Applications and of any and all applications related to the Priority Applications by priority claims (directly or indirectly), including any priority claims made and subject matter incorporated by reference therein as of the filing date of the instant application, is incorporated herein by reference to the extent such subject matter is not inconsistent herewith.
Number | Date | Country | |
---|---|---|---|
62267677 | Dec 2015 | US |