Open Access
15 April 2022 Measuring the soft x-ray quantum efficiency of a hybrid CMOS detector
Author Affiliations +
Abstract

Next-generation x-ray observatories, such as the Lynx X-ray Observatory Mission Concept or other similar concepts in the coming decade, will require detectors with high quantum efficiency (QE) across the soft x-ray band to observe the faint objects that drive their mission science objectives. Hybrid CMOS detectors (HCDs), a form of active-pixel sensor, are promising candidates for use on these missions because of their fast read-out, low power consumption, and intrinsic radiation hardness. We present QE measurements of a Teledyne H2RG HCD, performed using a gas-flow proportional counter as a reference detector. We find that this detector achieves high QE across the soft x-ray band, with an effective QE of 94.6  ±  1.1  %   at the Mn Kα  /  Kβ energies (5.90/6.49 keV), 98.3  ±  1.9  %   at the Al Kα energy (1.49 keV), 85.6  ±  2.8  %   at the O Kα energy (0.52 keV), and 61.3  ±  1.1  %   at the C Kα energy (0.28 keV). These values are in good agreement with our model, based on the absorption of detector layers. We find similar results in a more restrictive analysis considering only high-quality events, with only somewhat reduced QE at lower energies.

1.

Introduction

Next-generation x-ray observatories, potentially with capabilities similar to those of the proposed Lynx X-Ray Observatory,1,2 will provide unprecedented observations of the early universe. These missions will require large effective areas to enable detections of distant, faint sources. Detectors with capabilities surpassing those of the charge-coupled devices (CCDs) used in the focal planes of current missions will be required for future telescopes. The detectors used in these observatories will need high quantum efficiency (QE), the fraction of incident x-rays individually detected, to fully utilize the large collecting areas provided by the mirror assemblies and maximize sensitivity to faint sources. Future missions also require detectors with fast read-out to avoid pile-up, an effect in which multiple x-rays land in the same pixel during a single frame, when observing bright sources. Without a major increase in the detector frame rate, pile-up would be exacerbated by the greater effective areas of these missions.

Hybrid CMOS detectors (HCDs) are promising candidates for use in these future x-ray observatories. HCDs are composed of a silicon absorber layer bonded to a complementary metal–oxide–semiconductor (CMOS) read-out integrated circuit (ROIC). The absorber layer produces photoelectrons when absorbing an x-ray and the ROIC reads the amount of charge generated. Each pixel in an HCD contains its own ROIC, removing the need to transfer charge through other pixels. This provides significantly faster readout, intrinsic radiation hardness, and lower power demands compared with the CCDs on current missions. This hybridized structure, shown in Fig. 1, allows for independent optimization of the ROIC and absorber layers during fabrication. HCDs provide high QE across the soft x-ray band, due to their thin inactive surface layers and large depletion depths (100  μm in current models) enabled by the use of high-resistivity silicon in the absorber layer.

Fig. 1

Diagram detailing the structure of silicon HCDs. In these detectors, x-rays absorbed in the absorber layer create a charge cloud which is collected and transferred to the CMOS ROIC to be measured and read out. The energy of an x-ray will impact its absorption depth and the size of its charge cloud. Adapted from Ref. 3.

JATIS_8_2_026003_f001.png

While the Penn State High-Energy Astrophysics Detector and Instrumentation Laboratory has previously measured a high QE on an H1RG HCD at the Mn Kα/Kβ lines,3 further measurements are necessary to characterize the soft x-ray QE on other detectors and over a range of energies. Characterization of the QE at low x-ray energies (0.2 to 1.5 keV) is crucial to detector development, as good low-energy performance is vital to the science goals of future missions. For example, high QE at low energies would allow a mission like Lynx to more efficiently observe high-redshift quasars, whose power-law spectra have significant contribution in this band, and allow imaging of the halos of nearby galaxies, with low-energy thermal emission peaks.1 In this work, we present measurements of the QE of an H2RG HCD at moderate and low energies. In particular, we report on the QE of this H2RG HCD at the Mn Kα/Kβ (5.90/6.49 keV), Al Kα (1.49 keV), O Kα (0.52 keV), and C Kα (0.28 keV) characteristic x-ray emission lines. Much of this work was presented at the 2021 SPIE Optics and Photonics meeting and is adapted from the corresponding proceedings,4 with further expansions relative to that work.

2.

Quantum Efficiency Model

We use a one-dimensional (1D) slab absorption model to estimate the QE of HCDs. In this model, the QE is given by the fraction of x-rays transmitted through inactive surface layers and absorbed by the silicon absorber layer. This can be expressed in the following form:

Eq. (1)

QE(E)=(1eσSi(E)nSitabs)iNeσi(E)niti,
where tabs is the thickness of the absorbing layer (with the cross section and number density for silicon) and the product is over the atomic species in the inactive surface layer(s) with atomic cross sections σi, number densities ni, and corresponding layer thicknesses ti. The H2RG has a 100-μm Si absorber layer and a single inactive surface layer composed of SiO2 with a nominal thickness of 1028 Å5. In this work, we calculate transmission values using the absorption coefficients provided by Ref. 6. The resulting H2RG QE model is shown in Fig. 2.

Fig. 2

Estimate of the H2RG QE based on a 1D slab absorption model, using the nominal thicknesses of a 1028-Å SiO2 oxide layer and a 100-μm absorber layer. The QE decreases at lower energies due to the increasing fraction of x-rays absorbed by the oxide layer and decreases again at high energies due to an increasing fraction of x-rays which are transmitted through the entire absorbing layer.

JATIS_8_2_026003_f002.png

At energies below 5 keV, few x-rays are transmitted through the Si absorbing layer, and the QE is primarily determined by the transmittance of the inactive surface layer. Above 5 keV, almost all x-rays are transmitted through the inactive surface layer, and the QE is primarily determined by the fraction of x-rays absorbed by the Si absorbing layer. While the simple slab absorption model has proven accurate in past HCD QE measurements,3,7,8 effects ignored by this simple 1D model, such as charge traps near the Si/SiO2 interface, may cause deviations from model values. This model does not consider the detector response to absorbed x-rays, which will impact the QE achievable by a detector. Charge diffusion and excessive read noise may cause a suboptimal detector response to absorbed events, preventing their energies from being properly identified.9 These various effects can combine to complicate the response of pixelated silicon detectors, making measurements essential to truly understand their behaviour. QE measurements over a range of energies are therefore necessary to fully characterize the performance of detectors, to determine their suitability for future missions, and to guide the development of future detectors.

3.

Experimental Methods

QE measurements require a separate, well-characterized reference detector to determine the incident x-ray flux. We used a gas-flow proportional counter (PC) as our reference detector because of its known QE and high response across the soft x-ray band. During the experiment, both the H2RG and PC were simultaneously exposed at a similar distance to the same source. The flux measured by the H2RG was compared with the flux measured by the reference detector to determine the QE. The QE measurements were conducted in the 47-m x-ray beamline operated by the Penn State High-Energy Astrophysics Detector and Instrumentation Laboratory. The mounting structure containing both detectors can be seen in and outside of the beamline in Fig. 3. The large distance between the source and detectors as well as the simultaneous data collection minimized the difference in flux incident upon both detectors. The beamline is maintained at high vacuum (<105 mbar) to minimize attenuation of x-rays and avoid the condensation of water vapor onto the HCD (which is cooled during operation).

Fig. 3

(a) The front of the detector mounting structure shown outside the vacuum chamber, with the PC (left) and the H2RG (right), and (b) the back of the same structure shown in place for the QE measurements inside the beamline. A Si PIN diode detector was also mounted between these detectors but was not used during QE measurements.

JATIS_8_2_026003_f003.png

3.1.

H2RG

We conducted our QE measurements on a modified x-ray H2RG fabricated by Teledyne Imaging Sensors. This detector uses a standard H2RG ROIC, an array of 2048×2048 pixels with an 18  μm pitch, bonded to an array of 1024×1024 36μm pitch Si absorber cells, with only one out of every four ROIC pixels connected to the absorber layer. This design was selected to increase the distance between pixels to reduce the effect of interpixel capacitance, a form of crosstalk which decreased the energy resolution of older-generation H1RG HCDs.10 While extensive characterization measurements have been taken for this detector,11 and it even has flight heritage on the Water Recovery X-ray Rocket mission,12 no QE measurements had been made on the x-ray H2RG prior to this work.

X-ray HCDs must be cooled during operation to minimize dark current and instrumental background noise. The detector was cooled using a copper strap thermally connecting the device to an internal liquid nitrogen dewar and the temperature was regulated using a heater. When possible, the detector was maintained at 160 K during measurements. This temperature is slightly higher than the optimal operating temperature of 150  K, leading to increased noise and reduced energy resolution; however, we found that operating below 160 K allowed for deposition of an ice layer on the detector surface, which reduced the measured QE. This effect was observed as a gradual degradation of the QE with time. We installed a residual gas analyzer (RGA) to measure the partial pressure of water vapor in the chamber and confirmed that it was lower than the deposition pressure at 160 K during our QE measurements.13 During the lowest-energy measurement (C Kα), the detector was cooled to 140 K to suppress the instrument noise peak sufficiently to obtain satisfactory low-energy data. The deposition pressure of water vapor at this temperature is less than the partial pressure which was measured by the RGA during these measurements, allowing a layer of ice to form on the detector surface.

3.2.

Proportional Counter

The PC and H2RG were mounted side-by-side and equidistant from the midline of the chamber, as shown in Fig. 3. This arrangement provides roughly equal flux incident on both detectors. P-10 gas (90% argon and 10% methane, with a tolerance of 2%) was flowed through the PC at a rate of 200  cc/min. The temperature of the gas was measured to be 290±5  K. The pressure of the gas was recorded for each measurement (as it was affected by variations in atmospheric pressure), with a typical value of 1030 mbar and an estimated uncertainty of 2.7 mbar. The gas is separated from vacuum by a thin window installed in the aperture of the detector. The window was manufactured by Luxel Corporation for high transmittance of low-energy x-rays and was constructed from aluminized polyimide, supported by a stainless steel mesh.

For the PC to act as a reference detector and accurately measure the absolute flux produced by the source, its QE must be well-characterized. PCs are commonly used as reference detectors because their QE can be reliably determined through measurements of gas properties and window transmittance. We calculated the QE using a model similar to the HCD model described in Sec. 2, where the PC QE is given by the fraction of x-rays transmitted through the window which are then absorbed by the P-10 gas. The gas absorption was calculated based on the pressure, temperature, and composition of the P-10 gas, while the transmittance of the window in the PC aperture was measured to ensure accuracy.

Measurements of the window transmittance were necessary to obtain an accurate absolute calibration of the PC, as slight variations in the thickness of the window layers can cause significant changes in its transmittance. In the low-energy band (0.2 to 1.5 keV), almost all x-rays transmitted by the window are absorbed by the P-10 gas, so the QE of the detector is primarily a function of the window transmittance. Calibration measurements are also important at high energies, as the stainless steel mesh used to support the window acts as a binary mask, blocking x-rays at all energies throughout the soft x-ray band.

3.2.1.

Proportional counter calibration

We measured the transmittance of the PC window using a secondary window which was also manufactured by Luxel with the same specifications as the primary window. The secondary window was installed in the aperture of the PC during the calibration measurements. We then mounted the primary window directly in front of the secondary window and compared the count rates when the primary window was present and when it was removed. The ratio of the observed count rates in these measurements gives the transmittance of the window at a given energy.

We used radioactive sources because these provide a known flux, decaying at a known rate. An Fe-55 source was used for a measurement of the effective window transmittance at the Mn Kα/Kβ energies. We made transmittance measurements at lower energies using targets fluoresced by Polonium-210 alpha-particle sources. Measurements were made at the Al Kα, F Kα (0.68 keV), and C Kα lines using fluoresced aluminum, Teflon, and polyethylene targets, respectively. These measurements were directly used in the calculation of PC QE for the Mn Kα/Kβ, Al Kα, and C Kα H2RG QE calculations.

We did not make a measurement at the O Kα energy; however, we were able to fit a model to calculate the transmittance at this energy. In our model, the window layer thicknesses were fixed at the values provided by Luxel (listed in Table 1), and only the mesh transmittance was allowed to vary. A mesh transmittance of 0.780±0.005 provided the best fit to our measurements. We found that this fitted model agrees quite well with our measurements and used it to determine the window transmittance at the O Kα energy. This fit is shown in Fig. 4, along with the measured transmittance values and the modeled PC QE.

Table 1

Properties of the various layers of the PC window. Nominal thicknesses and estimates of the uncertainty for each layer were provided by Luxel.14 These values were used to model the PC window transmittance and calculate the PC QE at the O Kα energy. The densities of the aluminum and carbon layers are assumed to be 90% of their bulk densities because of the method of fabrication.

Window layerThickness (nm)Uncertainty (nm)
Carbon18.00.9
Polyimide440.022.0
Aluminum19.00.95
Aluminum oxide3.01.0

Fig. 4

Measurements of PC window transmittance at four energies shown on top of the fitted transmittance model (solid black line). The modeled PC QE is shown by the dashed gray line. Note that the QE and window transmittance models only diverge at higher energies, where a significant fraction of x-rays are transmitted through the gas undetected.

JATIS_8_2_026003_f004.png

3.3.

X-Ray Sources

We used two different sources of x-rays for our QE measurements. A radioactive Fe-55 source produced x-rays at the Mn Kα and Kβ energies, and a Manson electron-impact x-ray source produced lower-energy x-rays at the Al, O, and C Kα lines. These sources required careful calibration to ensure that both the energies and spatial distribution of the x-rays were well known. Both sources and their calibrations are described in this section.

3.3.1.

Fe-55 source

We used an Fe-55 radioactive source to produce emission at 5.9 and 6.5 keV. We mounted the Fe-55 source on the detector end of the beamline to produce a sufficient count rate. Unlike the Manson source, the Fe-55 source emits isotropically, so the incident flux is only a function of distance to the source. The Fe-55 source is not exactly equidistant from both detectors (2.450±0.010 m from the H2RG and 2.455±0.010 m from the PC), leading the flux incident on the H2RG to be 0.57±0.04% greater than the flux incident on the PC. The source was placed on an electronic linear stage so that it could be moved behind a shutter when not in use, allowing for dark or Manson measurements without opening the chamber.

The PC does not have sufficient energy resolution to distinguish between the Kα and Kβ peaks. Because the QE of the PC differs significantly at these energies (as shown in Fig. 4), knowing the fraction of incident x-rays at each energy is essential to accurately determining the fraction of x-rays absorbed by the P-10 gas. The intrinsic ratio of Mn Kα to Kβ emission is 0.1195;15 however, the emitted line ratio is increased by the greater transmission of Kβ x-rays through this Fe-55 source and its casing. We calculated the emitted Kβ to Kα x-ray ratio to be 0.132±0.002 based on the thicknesses of the Fe-55 source layer and Ni casing provided by the manufacturer.

3.3.2.

Manson source

A Manson ultrasoft x-ray source, a type of electron-impact source, was used to produce x-rays for measurements at lower energies. We used an aluminum anode in the Manson source for QE measurements at the Al Kα line (1.49 keV), a heavily oxidized magnesium anode for measurements at the O Kα line (0.52 keV), and a graphite anode for measurements at the C Kα line (0.28 keV). In addition to characteristic lines, the Manson source also emits continuum x-rays. These x-rays must be suppressed to prevent x-rays at different energies from being included in the characteristic x-ray fluxes used in the QE measurements. We installed carefully designed filters in the source end of the beamline to block continuum x-rays. These filters are also required to block optical light, produced by the Manson source, predominantly from the hot filament. If not blocked, this light would be registered by the detector, adding noise to our measurements and potentially damaging the CMOS ROIC. We used an aluminum filter, constructed from a piece of heavy-duty aluminum foil (25  μm), for the Al Kα measurements. The lower-energy measurements required more carefully designed filters made from thin foils and films, manufactured by Lebow Company, allowing for adequate transmittance of low-energy x-rays. The filter for the O Kα measurement was constructed from a 1-μm Cr foil supported by a 0.5-μm Mylar film, and the filter for the C Kα measurement was constructed from a 0.1-μm Al foil supported by a 5-μm parylene N film. Models of the transmittance for each filter are shown in Fig. 5. The filters do not completely eliminate the continuum emission; however, they provide sufficient separation between the characteristic x-ray and continuum peaks so both can be modeled and the counts can be separated.

Fig. 5

Model of filter transmittance as a function of energy for the filters used with the Manson source. The absorption edges make these filters extremely efficient at blocking continuum x-rays just above the characteristic energy.

JATIS_8_2_026003_f005.png

While the large distance between the Manson source and the detectors (46±0.1  m) allowed for fairly uniform flux incident on the detectors, we found that it was neither perfectly uniform nor perfectly symmetric. The detectors were positioned on the mounting structure such that the distances from each detector to the source would be identical to a tolerance of 1 cm, much too small to explain the observed differences in incident flux. These spatial variations are instead thought to result from inherent properties of the source setup (e.g., variations in the filter thicknesses). We obtained a calibration for each Manson source and filter configuration using a detector mounted on a linear stage. We measured the count rate at the position of both the H2RG and PC. These measurements were repeated to account for variability in the absolute flux produced by the Manson source. We found that the percentage of the flux incident upon the H2RG relative to that incident on the PC was 99.1±1.0% in the Al Kα configuration, 100.0±0.2% in the O Kα configuration, and 99.3±0.4% in the C Kα configuration.

4.

Analysis and Results

During data acquisition, both detectors measured the incident flux produced by a given source simultaneously. At each of the energies measured, we recorded data with the detectors exposed to the source then subsequently recorded background measurements with the source blocked. Each measurement lasted for a duration of ∼10 min (with exact live times for the detectors differing slightly). In this section, we describe the methods used to analyze this data, discuss the QE calculation, and present the results of our measurements.

4.1.

H2RG Data Reduction

During these measurements, we operated the H2RG in “up-the-ramp” mode, in which many frames were taken successively without clearing charge using nondestructive readout. The charge accumulated in a single frame was found by taking the difference in pixel values between consecutive frames. Row and channel noise were then removed from the image using a boxcar smoothing algorithm. To identify x-ray events, we found pixels which had charge depositions that were local maxima and above a set “event” threshold. If the charge in a surrounding pixel exceeded a “split” threshold, it was included in the total charge of the x-ray event.

We also used the split threshold to assign a grade to each event, based on the distribution of pixels included in the events. This allows us to characterize the degree of charge spreading and enables discrimination between events corresponding to x-rays and those caused by cosmic rays or noise. The selection of allowed grades has a significant effect on the measured QE. Including all events, regardless of grade, will give an idealized QE for a detector, while a flight-like analysis which restricts the allowed event grades may give a more realistic QE. In this work, we investigate both approaches, measuring the QE when considering all event grades and when only considering events with high-quality grades. For our more restrictive analysis, we select events with grades 0 to 12 in the Swift XRT photon-counting mode grading scheme,16 allowing for the inclusion of 1 to 4 pixel events with shapes expected for x-ray events.

We use different event and split thresholds for each energy to optimally exclude noise and include charge in each analysis. These different values are shown in Table 2 in terms of the measured detector noise (σn). We use high thresholds in the Mn Kα/Kβ and Al Kα analyses to clearly separate pixels with charge from x-ray events from those with spurious values due to noise. The lower-energy O Kα and C Kα analyses require lower event thresholds, as higher values can miss a significant fraction of events and reduce the measured QE. We also reduce the split thresholds to more accurately sum charge for low-energy events. Given these 2σn split thresholds, we expect a considerable number of x-ray events to include split pixels with high values due to random noise and not charge spreading, artificially skewing the energy distribution and giving some events low-quality grades.

Table 2

Thresholds used for x-ray event identification and evaluation at different energies, given in terms of the detector noise (σn). Note that the noise value for the C Kα measurement is lower due to the colder detector operating temperature used during this measurement.

Emission lineσn (e−)Event threshold (σn)Split threshold (σn)
Mn Kα/Kβ12105
Al Kα1284
O Kα1232
C Kα642

Our H2RG is an engineering-grade detector, with a greater number of defects in the silicon lattice of the absorber layer compared to a science-grade sensor. It has many charge traps which create areas with significant dark current and limited dynamic range. We chose to exclude x-ray events from pixels in regions associated with these charge traps, as these pixels may have reduced sensitivity to incident x-rays. A science-grade detector would have fewer defects, so this measurement better reflects the expected QE of the detectors which would be used on a mission. We also excluded events from pixels too close to the edge of the detector, where edge effects reduce the detector sensitivity. After masking off these pixels, we are left with 85.7% of the total detector area as active. The area of the detector was scaled by this value in calculations of the measured flux and the QE.

4.2.

Mn Kα/Kβ Measurement

Figure 6 shows the H2RG and PC spectra from the Fe-55 QE measurement. Both the H2RG analysis including all event grades and that including only grades 0 to 12 are shown in this plot, revealing only a small reduction in the number of events due to the grades cut. Due to the poor energy resolution of the PC and H2RG (when operated at 160 K), the Mn Kα and Kβ peaks are blended together in the spectra from both detectors. For this reason, we make an effective QE measurement at the combined Mn Kα and Kβ energies, giving the fraction of all incident Mn K x-rays detected.

Fig. 6

H2RG and PC spectra from the Fe-55 QE measurements. Both spectra include combined Mn Kα/Kβ peaks with escape peaks at lower energies. The (scaled) background spectra from the dark measurement are shown by the dashed lines. The bounds used to sum the counts and determine the fluxes are shown by the dotted red lines. The H2RG spectra with all event grades are shown in black, while the spectra with grades 0 to 12 are shown in blue. Applying this grades cut did not significantly change the Fe-55 or background spectra, as shown by the nearly overlapping lines. The location of the Mn Kα peak from this spectrum was used to calculate the gain used in all of the plotted H2RG spectra in this work (PC gain varied between measurements and was fit separately in each measurement).

JATIS_8_2_026003_f006.png

We determine the count rates for both detectors by summing all counts between specified bounds. The sums of the counts in the background spectra within these bounds are subtracted from these count rates, removing the instrumental backgrounds and leaving only the count rates from real x-ray events. This is an appropriate method for measuring the count rate of Mn K x-rays, as no other sources of characteristic or continuum emission were present during this measurement. Spectra from both detectors include escape peaks, as the Mn Kα/Kβ x-rays are more energetic than the binding energies of the Ar and Si K shells. The escape peaks were included within the bounds of the sums because they correspond to real Mn K x-rays registered by the detectors whose energies can be accurately determined given a proper detector response matrix. The H2RG spectrum also includes substantial high- and low-energy tails, likely due to thresholding effects; however, as most counts from these features are included in our bounds, they should have little effect on the measured QE.

4.3.

Al Kα, O Kα, and C Kα Measurements

The Al Kα, O Kα, and C Kα spectra contain continuum emission from the Manson source in addition to the characteristic x-ray peaks. This continuum emission must be separated out to obtain accurate QE measurements. For this reason, we chose to model and fit these spectra, allowing us to use only the model components corresponding to the characteristic emission in our QE analysis. We removed instrument background counts from each of the spectra before fitting by subtracting the binned counts of the dark spectra (scaled based on live time) from those of the measurement spectra. The models, described below, were fit using a least-squares algorithm. The total number of counts was found by integrating the area under the fitted curves of the model components corresponding to characteristic x-ray peaks, with lower integration bounds of zero and upper bounds set individually to fully assess the flux for each measurement.

4.3.1.

PC Models

The PC spectra, fitted models, and model components are shown in Fig. 7. The primary features of these spectra are skewed characteristic peaks with low-energy shelves along with a peak from continuum emission at higher energies. The low-energy shelves are the results of events with partial charge collection, where a fraction of the charge recombines before amplification and collection. We modeled the characteristic peaks and low-energy shelves using the distributions described in Ref. 17. The characteristic peaks were modeled as Prescott functions:18

Eq. (2)

P(x)=P0x01/4x3/44πQexp[(xx0)2/Q],
where x0 is a parameter which defines the location of the peak, Q defines its width, and P0 defines its amplitude. The low-energy shelves were modeled as the following:

Eq. (3)

S(x)=(a+bx)(1erf[xx022x0Q]),
where x0 and Q are parameters from the Prescott function and a and b are parameters defining the amplitude and slope of the shelf. In the Al Kα spectrum, the slope of the shelf (b) is set to zero because this parameter is not needed to achieve a good fit. The higher-energy continuum peaks are modeled as skew-normal curves in the Al and O Kα spectra and as a Gaussian curve in the C Kα spectrum. The high-energy peak in the O Kα spectrum also includes Mg Kα counts; however, a single skew-normal curve was sufficient to model both the continuum and Mg Kα components due to the limited spectral resolution of the PC.

Fig. 7

PC spectra for the Al, O, and C Kα QE measurements along with the spectral models (red dash-dotted lines) and residuals plotted below. The model components included in the PC characteristic flux (Prescott function and low-energy shelf) are shown by the blue dashed lines while the excluded continuum flux is shown by the gray dotted lines.

JATIS_8_2_026003_f007.png

4.3.2.

H2RG Models

We fit two different models for each of the H2RG spectra, one for the spectrum including all events and another for those with only grades 0 to 12. We run two separate analyses for the different grades filters, fitting a model separately for each spectrum. The Al Kα spectra for both grades filters are shown in Fig. 8 along with the fitted models and model components. In both analyses, the characteristic peak is modeled as a sum of two Gaussians, and the high-energy continuum peak is modeled as a skew-normal curve. We find that the skewness of the main peak at low energies, necessitating the addition of a second Gaussian, is present regardless of the split threshold or grade selection. This leads us to believe that it is a real feature of the spectrum and not simply a feature of the detector response. We believe this is due to the presence of continuum emission which is transmitted through the filter between energies of 1 and 1.5 keV (see Fig. 5). Even if this is the case, the inclusion of this continuum should not significantly bias our measurement, as the filter transmission drops sharply at lower energies and the H2RG QE in this energy range is fairly uniform (see Fig. 2).

Fig. 8

The H2RG spectra from the Al Kα measurement along with the models and residuals. The full model is shown as the red dash-dotted line, while the characteristic model components are shown by the blue dashed lines and the excluded continuum model is shown by the gray dotted line. Both grades filters result in similar spectra, with slightly better differentiation between the characteristic and continuum peaks in the graded spectrum.

JATIS_8_2_026003_f008.png

Figure 9 shows the spectra from the O Kα measurement. The characteristic oxygen peaks are modeled as simple Gaussians. For this measurement, we used an oxidized magnesium anode in the Manson source, so the high-energy peak contains a Mg Kα line (1.25 keV) in addition to continuum emission. We modeled this peak as a combination of two Gaussians.

Fig. 9

H2RG spectra from the O Kα measurement along with the models (red dash-dotted lines) and residuals. Each of the model components is shown, including the Gaussian curve for the characteristic peak (blue dashed line) and the two Gaussians for the continuum/Mg Kα peak (gray dotted lines).

JATIS_8_2_026003_f009.png

Figure 10 shows the spectra from the C Kα measurement. In the all-events model, the C Kα peak was modeled as a mixture of a skew-normal and normal curve. This combination was selected to sum all C Kα events, matching the skewness of the peak caused by the interaction of noise and thresholding effects. After applying the grades cut, we find that the skewness of the C Kα peak is reduced while the high-energy tail is still present. We ignore the tail in the grades 0 to 12 analysis, using a single Gaussian to model the C Kα peak, to focus on only high-quality events whose energy could be accurately determined from a simple analysis. Including events in the tail would lead to a 0.4% increase in the count rate and QE in the high-quality event analysis. In both analyses, the continuum peak is modeled as a Gaussian.

Fig. 10

H2RG spectra from the C Kα measurement along with the models and residuals. The full models are shown by the red dash-dotted lines, with the included model components shown by blue dashed lines and the excluded components shown by gray dotted lines. Note that we used different models for the C Kα peak in the all event grades and grades 0 to 12 analyses.

JATIS_8_2_026003_f010.png

4.4.

QE Calculation

The QE is given by the ratio of the flux measured by the H2RG to the flux incident on the detector. We determine the absolute flux incident on the PC using the flux measured by the PC and the modeled PC QE (based on window transmission measurements and gas absorption). We calculate the flux incident on the H2RG by scaling the incident PC flux by an incident flux ratio (k), determined using the Manson beam calibration measurements (described in Sec. 3.3.2) or the distances of the detectors from the Fe-55 source. These factors are combined to obtain the following expression for the H2RG QE:

Eq. (4)

QEH2RG=QEPCFH2RGkFPC.
The values used in this calculation are shown in Table 3 for each of the different energies.

Table 3

Values used for the calculation of the H2RG QE, including the fluxes measured by both detectors, the modeled PC QE, and the ratio of the flux incident on the H2RG to that incident on the PC. Fluxes are given in units of counts s−1  cm−2. The incident flux ratio for the C Kα measurement accounts for the estimated ice attenuation in addition to the measured nonuniformity of the beam.

Emission linePC fluxPC QEIncident flux ratioH2RG flux
Grades 0 to 12All event grades
Mn Kα/Kβ13.41±0.0769.3±0.7%1.0057±0.000417.91±0.0518.40±0.06
Al Kα39.42±0.2474.9±1.1%0.991±0.01050.46±0.1451.30±0.15
O Kα19.18±0.0938.5±1.2%1.000±0.00235.61±0.1642.70±0.22
C Kα86.14±0.3653.3±0.8%0.950±0.00879.71±0.5694.10±0.39

The C Kα analysis is complicated by the presence of an ice layer on the detector during this measurement, due to the lower operating temperature used. We correct for the resulting attenuation by measuring the thickness of the ice layer through observations of the decrease in the QE over the course of multiple hours. We find that the thickness of the ice layer increased linearly with time, allowing us to infer the thickness of the layer by determining the accumulation rate and the amount of time below the deposition temperature. Our measurements are shown along with the resulting fit in Fig. 11. We determine that the thickness of the ice layer was 83±14 nm at the time of our measurement, leading to a transmission of 95.6±0.7% for C Kα x-rays through the ice layer.

Fig. 11

Measurement of the change in ice thickness as a function of time after the first measurement. The data points come from the observed decrease in QE relative to this initial measurement. The dashed line shows the linear fit used to determine the rate of ice accumulation with the uncertainty indicated by the gray shaded region. This rate was used to calculate the thickness of ice on the detector surface during the C Kα measurement and correct for attenuation due to this layer.

JATIS_8_2_026003_f011.png

4.5.

Error Analysis

We compute the error on our QE measurements using a Monte Carlo analysis, in which we repeatedly sample input variables from their respective distributions to sample the QE error distribution. These input variables include the observed count rates, PC dimensions and gas properties, and incident flux and PC window calibration measurements. We use 105 simulations, a large enough sample to provide accurate error estimates. We test the resulting distributions for normality using a Shapiro–Wilk test19 and find that all our QE measurement errors are consistent with Gaussian distributions with P-values greater than 0.05. All uncertainty values provided in this work are the 1σ values taken from these distributions.

4.6.

QE Results

Our QE analysis including all event grades yields results that are in good agreement with our simple slab-absorption model. These results are shown along with the model values (using the nominal values of the absorber and oxide layer thicknesses) in Table 4 and in Fig. 12. At the Mn Kα/Kβ energies, we measured the effective QE to be 94.6±1.1%. This value is in good agreement with our modeled effective QE value of 96.2%, computed by taking the weighted average of the modeled QE at the Mn Kα and Kβ energies (using the calculated Mn Kα/Kβ ratio specific to our source). It is also consistent with the previous measurement made on an H1RG HCD, with a value of 97±5%.3 This H2RG is expected to have similar QE to H1RGs at moderate energies due to their identical absorber layer thicknesses (100μm), even as the H1RGs have thinner (250 Å) oxide layers and some devices have deposited aluminum optical blocking layers.

Table 4

Comparison of H2RG QE measurements from the analysis with high-quality grades (grades 0 to 12) and with all event grades, shown along the model QE values, which give the fraction of all incident x-ray events expected to be registered by the detector.

Emission lineEnergy (keV)Measured QE (grades 0 to 12)Measured QE (all event grades)Model QE (all x-ray events)
Mn Kα/Kβ5.90/6.4992.1±1.1%94.6±1.1%96.2%
Al Kα1.4996.7±1.9%98.3±1.9%97.5%
O Kα0.5271.3±2.3%85.6±2.8%89.8%
C Kα0.2851.9±1.0%61.3±1.1%62.4%

Fig. 12

Measured QE values (including all event grades) overplotted with the slab-absorption H2RG QE model (using nominal thicknesses of detector layers). The Fe-55 measurement is plotted at the weighted mean of the Mn Kα and Kβ energies.

JATIS_8_2_026003_f012.png

We measured the QE at the Al Kα energy to be 98.3±1.9%, also showing excellent agreement with our model and with previous H1RG measurements, which demonstrated QE90% for these detectors between 1 and 3 keV.7 As previously noted, a small fraction of continuum x-ray events with energies between 1 and 1.5 keV may be included in this analysis. This should not significantly affect the analysis, as the QE at these energies should not differ significantly from that at the Al Kα energy. We find that the QE is still high at low energies and remains consistent with our model, with measured values of 85.6±2.8% at the O Kα energy and 61.3±1.1% at the C Kα energy. As the absorption length for these energies is quite small, these measurements indicate that the charge collection efficiency is still high for x-rays absorbed near the surface of the absorber layer.

Our analysis of high-quality events with grades 0 to 12 still yields high QE, showing only a slight reduction due to the exclusion of events with low-quality grades. In this analysis, we find an effective QE of 92.1±1.1% at the Mn Kα/Kβ energies. At the Al Kα energy, we find a QE of 96.7±1.9%. At these higher energies, the small difference in QE from the two analyses likely results from events excluded due to excess charge spreading. The combination of a 15 V substrate voltage, large pixel pitch (36  μm), and high split thresholds makes this a minor effect. A larger fraction of events are excluded by the grades cut in the lower-energy analyses, as more pixels with high values due to noise are able to trigger the lower split threshold, resulting in events with low-quality grades. We find a QE of 71.3±2.3% at the O Kα energy and a value of 51.9±1.0% at the C Kα energy. The next generation of HCDs is expected to have lower read noise and should be able to achieve QE values which are impacted less by grade cuts, as fewer split pixels would be triggered by noise. These future devices should also have a more Gaussian response at low energies, as thresholds can be better optimized to discriminate pixels to which charge has spread from those with noise.

These measurements represent high QE, comparable to state-of-the-art back-illuminated x-ray CCDs, such as those used in Chandra ACIS,20 XMM-Newton EPIC,21,22 or eROSITA.23,24 Our H2RG absorber layer was fabricated with a design targeted for optical detectors and was designed with a thicker-than-necessary oxide layer (acting as an optical antireflection coating). Other x-ray HCDs with thinner oxide layers should theoretically be capable of even higher low-energy QE than that measured on our H2RG. We anticipate future x-ray HCD models, with lower read noise and absorber layers fabricated explicitly for x-ray applications, to achieve a substantially better low-energy response. The Lynx High-Definition X-ray Imager25 would need high QE, with requirements of QE>85% between 0.5 and 10 keV and QE>10% between 0.2 and 0.5 keV.1 Our measurements indicate that HCDs are capable of meeting these requirements; however, increasing the thickness of the absorber layer may be necessary to achieve the required QE at the highest energies and decreasing the read noise may be required to achieve higher QE at lower energies when only considering high-quality events.

5.

Rotation Measurement

As further verification of our low-energy QE measurements and our slab-absorption model, we conducted a rotation experiment to measure the relative QE at the C Kα energy as a function of incident angle. Changing the incident angle changes the thickness of the oxide layer a photon must traverse to be detected, allowing us to verify the thickness of the oxide layer and probe the detector response to x-rays absorbed at various depths in the detector.

For the experiment, we mounted the H2RG HCD, internal dewar, and electronics on a rotatable gimbal, shown in Fig. 13. We positioned the detector such that its face was centered on the gimbal’s axis of rotation, allowing only its angle relative to the source, and not its position, to change between measurements. The PC was mounted beside the gimbal and did not rotate with the H2RG assembly. It was used as a reference detector to correct for variations in the flux produced by the source. We used the same source configuration as was used for the C Kα QE measurement, with the Manson source (with the graphite anode) and Al/parylene N filter.

Fig. 13

Setup of the rotation experiment inside the beamline. The detector can be seen mounted on the gimbal in front of its readout electronics and liquid nitrogen dewar, rotated to an angle of 60  deg off axis. The PC is also visible on the right side of the image.

JATIS_8_2_026003_f013.png

During the measurement, we recorded the count rate of the H2RG as a function of incident angle relative to the source. Our analysis for this measurement was similar to that for the all-events carbon analysis, except we fixed parameters of the higher-energy normal component of the C Kα model to scale with those of the main skew-normal component, allowing for better fits to the higher-angle measurements with fewer counts. We also corrected for attenuation due to the ice layer, in a similar manner to the correction described in Sec. 4.4.

Figure 14 shows the results of the rotation measurement, giving the count rate at each angle relative to an initial face-on measurement, corrected for variation in the source intensity and attenuation by the ice layer. Our model accounts for the reduced solid area visible to the source (by a factor of cos(i), where i is the angle of incidence) and greater attenuation as x-rays travel through a larger distance in the oxide layer (d1/cos(i)). We fit this model to our rotation measurements, allowing the thickness of the oxide layer to vary, and found a best-fit value of 1151±103 Å, in reasonable agreement with the nominal value of 1028 Å SiO2. This best-fit model is shown along with the nominal model in Fig. 14. This rotation measurement provides further validation of the slab-absorption model, indicating that it accurately describes charge collection in the H2RG, and provides verification of our low-energy QE results.

Fig. 14

Rotation experiment measurements showing the count rate relative to a prior measurement (face-on at 0 deg), corrected for variation in the Manson flux and ice attenuation. The gray dashed line shows the expected decrease in count rate due to the decrease in solid area exposed to the source. The model for the expected decrease in count rate due to both the decrease in solid area and the increase in path length through the oxide layer is shown for both the nominal (red dash-dotted) and best-fit (solid black) SiO2 layer thicknesses. Residuals are shown for the best-fit model.

JATIS_8_2_026003_f014.png

6.

Conclusion

We have measured the QE of a Teledyne Imaging Sensors H2RG HCD at multiple energies using a gas-flow PC as a reference detector. Our results show that this detector is capable of achieving high QE across the soft x-ray band, with measured values of 94.6±1.1% at the Mn Kα/Kβ energies (5.90/6.49 keV), 98.3±1.9% at the Al Kα energy (1.49 keV), 85.6±2.8% at the O Kα energy (0.52 keV), and 61.3±1.1% at the C Kα energy (0.28 keV) when including all x-ray events grades. These values are consistent with our model, which gives the QE as the fraction of events transmitted by the surface oxide layer that are then absorbed by the absorber layer. Our measurements of the relative QE at the C Kα energy as a function of angle are also in agreement with this model. These measurements demonstrate that the detector maintains high charge collection efficiency across the absorber layer and that the QE is not significantly degraded by other factors unaccounted for in the model.

We find the detector still maintains high QE when considering only events with high-quality grades (i.e., grade 0 to 12 x-ray events, those preferred for studies requiring higher spectral resolution). In this more restrictive analysis, we measured the QE to be 92.1±1.1% at the Mn Kα/Kβ energies, 96.7±1.9% at the Al Kα energy, 71.3±2.3% at the O Kα energy, and 51.9±1.0% at the C Kα energy. The reduction of QE due to the grades cut is larger at the lower O Kα and C Kα energies, due to noise artificially triggering split-pixel events. Future devices with lower read noise should achieve a better low-energy response. In addition to reducing the instrumental background, this would allow fewer low-energy events to include split pixels triggered by noise, allowing the QE for high-quality events to approach the higher values for all events presented in this work.

HCDs are capable of providing high QE at soft x-ray energies, helping enable future missions to attain the large effective areas necessary to achieve their ambitious science goals. The faster frame rates enabled by HCDs and other active-pixel sensors will allow thinner optical and UV-blocking filters to sufficiently suppress these backgrounds, thus enabling larger effective areas. We also anticipate future models of HCDs with lower noise and thinner surface oxide layers. These factors should further improve the low-energy QE achievable by HCDs, making them excellent candidates for missions which require good response across the soft x-ray band.

Acknowledgments

This work was supported by a NASA Space Technology Graduate Research Opportunity (Grant No. 80NSSC20K1210), as well as NASA Grant No. 80NSSC20K0778. We would also like to acknowledge useful discussions and advice from Yibin Bai at Teledyne Imaging Sensors and Ben Zeiger at Luxel. We would also like to thank David Palmer for suggesting the rotation measurement as a method of characterizing detector QE. In addition, we would like to thank the anonymous reviewers for their helpful comments and suggestions. The authors have identified no financial interests or other potential conflicts of interest to disclose.

References

1. 

The Lynx Team, “The Lynx Mission concept study final report,” (2019). Google Scholar

2. 

J. A. Gaskin et al., “Lynx x-ray observatory: an overview,” J. Astron. Telesc. Instrum. Syst., 5 021001 (2019). https://doi.org/10.1117/1.JATIS.5.2.021001 Google Scholar

3. 

S. D. Bongiorno et al., “Measuring the x-ray quantum efficiency of a hybrid CMOS detector with 55Fe,” Nucl. Instrum. Methods Phys. Res., 785 87 –92 (2015). https://doi.org/10.1016/j.nima.2015.02.043 NIMRD9 0167-5087 Google Scholar

4. 

J. M. Colosimo et al., “Measuring the quantum efficiency of x-ray hybrid CMOS detectors,” Proc. SPIE, 11821 118210T (2021). https://doi.org/10.1117/12.2594040 PSISDG 0277-786X Google Scholar

5. 

Y. Bai, (2021). Google Scholar

6. 

B. L. Henke, E. M. Gullikson and J. C. Davis, “X-ray interactions: photoabsorption, scattering, transmission, and reflection at E = 50-30,000 eV, Z = 1-92,” At. Data Nucl. Data Tables, 54 181 –342 (1993). https://doi.org/10.1006/adnd.1993.1013 ADNDAT 0092-640X Google Scholar

7. 

Z. R. Prieskorn et al., “Soft x-ray quantum efficiency of silicon hybrid CMOS detectors,” Proc. SPIE, 9154 915410 (2014). https://doi.org/10.1117/12.2056566 PSISDG 0277-786X Google Scholar

8. 

A. T. Kenter et al., “Rockwell CMOS hybrid imager as a soft x-ray imaging spectrometer,” Proc. SPIE, 5898 58981P (2005). https://doi.org/10.1117/12.632560 PSISDG 0277-786X Google Scholar

9. 

E. D. Miller et al., “The effects of charge diffusion on soft x-ray response for future high-resolution imagers,” Proc. SPIE, 10699 106995R (2018). https://doi.org/10.1117/12.2313565 PSISDG 0277-786X Google Scholar

10. 

A. D. Falcone et al., “Recent progress on developments and characterization of hybrid CMOS X-ray detectors,” Proc. SPIE, 8453 84530E (2012). https://doi.org/10.1117/12.925518 PSISDG 0277-786X Google Scholar

11. 

Z. Prieskorn et al., “Characterization of Si hybrid CMOS detectors for use in the soft x-ray band,” Nucl. Instrum. Methods Phys. Res., 717 83 –93 (2013). https://doi.org/10.1016/j.nima.2013.03.057 NIMRD9 0167-5087 Google Scholar

12. 

M. Wages et al., “Flight camera package design, calibration, and performance for the water recovery x-ray rocket mission,” Proc. SPIE, 11118 111180D (2019). https://doi.org/10.1117/12.2529361 PSISDG 0277-786X Google Scholar

13. 

W. Wagner et al., “New equations for the sublimation pressure and melting pressure of H2O Ice Ih,” J. Phys. Chem. Ref. Data, 40 043103 (2011). https://doi.org/10.1063/1.3657937 JPCRBU 0047-2689 Google Scholar

14. 

B. Zeiger, (2019). Google Scholar

15. 

J. H. Scofield, “Relativistic Hartree–Slater values for K and L X-ray emission rates,” At. Data Nucl. Data Tables, 14 121 (1974). https://doi.org/10.1016/S0092-640X(74)80019-7 ADNDAT 0092-640X Google Scholar

16. 

D. N. Burrows et al., “The swift x-ray telescope,” Space Sci. Rev., 120 (3-4), 165 –195 (2005). https://doi.org/10.1007/s11214-005-5097-2 SPSRA4 0038-6308 Google Scholar

17. 

J. M. Auerhammer et al., “High-accuracy calibration of the HXDS flow-proportional counter for AXAF at the PTB laboratory at BESSY,” Proc. SPIE, 3444 19 –29 (1998). https://doi.org/10.1117/12.331242 PSISDG 0277-786X Google Scholar

18. 

J. R. Prescott, “Photomultiplier single-electron statistics and the shape of the ideal scintillation line,” Nucl. Instr. Meth., 22 256 –268 (1963). https://doi.org/10.1016/0029-554X(63)90253-2 Google Scholar

19. 

S. S. Shapiro and M. B. Wilk, “An analysis of variance test for normality (complete samples),” Biometrika, 52 591 –611 (1965). https://doi.org/10.1093/biomet/52.3-4.591 BIOKAX 0006-3444 Google Scholar

20. 

G. P. Garmire et al., “Advanced CCD imaging spectrometer (ACIS) instrument on the Chandra X-ray Observatory,” Proc. SPIE, 4851 28 –44 (2003). https://doi.org/10.1117/12.461599 PSISDG 0277-786X Google Scholar

21. 

R. Hartmann et al., “Quantum efficiency of the XMM pn-CCD camera,” Proc. SPIE, 3765 703 –713 (1999). https://doi.org/10.1117/12.366552 PSISDG 0277-786X Google Scholar

22. 

L. Strüder et al., “The European photon imaging camera on XMM-Newton: the pn-CCD camera,” Astron. Astrophys., 365 L18 –L26 (2001). https://doi.org/10.1051/0004-6361:20000066 AAEJAF 0004-6361 Google Scholar

23. 

S. Ebermayer et al., “Quantum efficiency measurements of eROSITA pnCCDs,” Proc. SPIE, 7742 77420U (2010). https://doi.org/10.1117/12.856939 PSISDG 0277-786X Google Scholar

24. 

N. Meidinger et al., “eROSITA camera array on the SRG satellite,” J. Astron. Telesc. Instrum. Syst., 7 025004 (2021). https://doi.org/10.1117/1.JATIS.7.2.025004 Google Scholar

25. 

A. D. Falcone et al., “Overview of the high-definition x-ray imager instrument on the Lynx X-ray surveyor,” J. Astron. Telesc. Instrum. Syst., 5 021019 (2019). https://doi.org/10.1117/1.JATIS.5.2.021019 Google Scholar

Biography

Joseph M. Colosimo is a graduate student at the Pennsylvania State University, currently working toward his PhD in astronomy and astrophysics. He received his BS degree in physics from Georgia Institute of Technology in 2018. His research interests include the characterization of x-ray detectors and their applications to high-energy astrophysics.

Abraham D. Falcone is a research professor of astronomy and astrophysics at Pennsylvania State University. He is the author of more than 200 refereed journal publications. He leads research in fields ranging from x-ray and gamma-ray astrophysics instrumentation to the study of particle acceleration at high energy astrophysical sites such as active galactic nuclei, gamma ray bursts, and x-ray binaries. He is a member of the Swift team that was awarded the Bruno Rossi prize in 2007, as well as a VERITAS collaboration member, and a collaborator on multiple future high energy astronomy mission development efforts. He is the principal investigator for the NASA BlackCAT CubeSat mission, as well as several astronomical detector development programs.

Mitchell Wages is a research and development engineer at the Pennsylvania State University. He received his BS degree in astronomy and astrophysics from the Pennsylvania State University in 2016. His research interests include x-ray instrumentation, vacuum technologies, CubeSats, and mechanical design.

Samuel V. Hull is a postdoctoral program fellow at NASA’s Goddard Space Flight Center. He received his BA degree in physics and mathematics from Washington University in St. Louis in 2015, and his PhD in astronomy and astrophysics from Pennsylvania State University in 2021. His primary research interest is in the development of astronomical x-ray detectors for future missions, including both solid-state detectors and transition-edge sensors.

Daniel M. LaRocca is a postdoctoral scholar at the Pennsylvania State University. He received his PhD in physics from the University of Iowa in 2020 and his BA degree in physics from Illinois Wesleyan University in 2013. He was a principal member of the HaloSat team and now is a member of the BlackCAT team. His research interests include x-ray detectors, x-ray astronomy using small satellites or CubeSats, diffuse x-ray emission, and multi-messenger astrophysics.

Biographies of the other authors are not available.

CC BY: © The Authors. Published by SPIE under a Creative Commons Attribution 4.0 Unported License. Distribution or reproduction of this work in whole or in part requires full attribution of the original publication, including its DOI.
Joseph M. Colosimo, Abraham D. Falcone, Mitchell J. Wages, Samuel V. Hull, Daniel M. LaRocca, David N. Burrows, Cole R. Armstrong, Gooderham McCormick, Mitchell Range, and Fredric Hancock "Measuring the soft x-ray quantum efficiency of a hybrid CMOS detector," Journal of Astronomical Telescopes, Instruments, and Systems 8(2), 026003 (15 April 2022). https://doi.org/10.1117/1.JATIS.8.2.026003
Received: 16 November 2021; Accepted: 25 March 2022; Published: 15 April 2022
Lens.org Logo
CITATIONS
Cited by 2 scholarly publications.
Advertisement
Advertisement
KEYWORDS
Quantum efficiency

X-rays

Sensors

X-ray detectors

Aluminum

CMOS sensors

Manganese

RELATED CONTENT


Back to Top