• Journal of Infrared and Millimeter Waves
  • Vol. 43, Issue 2, 268 (2024)
Shuang WU1、2, Qing-Hua LIANG1, Hong-Lei CHEN1, and Rui-Jun DING1、*
Author Affiliations
  • 1Key Laboratory of Infrared Imaging Materials and Detectors,Shanghai Institute of Technical Physics,Chinese Academy of Sciences,Shanghai 200083,China
  • 2University of Chinese Academy of Sciences,Beijing 100049,China
  • show less
    DOI: 10.11972/j.issn.1001-9014.2024.02.017 Cite this Article
    Shuang WU, Qing-Hua LIANG, Hong-Lei CHEN, Rui-Jun DING. Low noise ROIC integrated with correlated double sampling with adjustable intervals for hyperspectral applications[J]. Journal of Infrared and Millimeter Waves, 2024, 43(2): 268 Copy Citation Text show less

    Abstract

    Low noise is a key requirement of readout integrated circuit (ROIC) in hyperspectral applications for its low radiation. Correlated double sampling (CDS) is commonly used to suppress noise. In this paper, CDS is improved by adjusting the time interval between the clamp and sample-and-hold (SH), which can filter low-frequency noise flexibly. A 640×512, 15 μm pixel pitch ROIC is designed and fabricated in 180 nm CMOS process. The input stage consists of low-noise capacitive trans-impedance amplifier (CTIA) and CDS with adjustable intervals (AICDS). A timing generator is proposed to extend the CDS reset time from 0 to 270 clock cycles. By extending the reset time to decrease the time interval, the noise electrons are significantly decreased from 39 e- to 18.3 e-. The SPECTRE simulation and the experimental results corroborate that the proposed structure AICDS can optimize noise performance of hyperspectral ROIC, thus can be widely used.

    Introduction

    Hyperspectral technologies have played a great role in vegetation monitoring,water resource management,geology and land cover1. It acquires continuous,narrow-band image data with a high spectral resolution as shown in Fig. 1(a). On one hand,it "captures" most of the subtle changes in the spectrum of a feature2,which can be used to identify,classify,or quantitatively analyze substances. On the other hand,the light radiation dispersed and focused on hyperspectral infrared focal plane arrays (IRFPA),as presented in Fig. 1(b),is extremely low,resulting in high requirements for low noise. Table 13-6 shows the comparison of several hyperspectral sensors launched in recent years. The readout noise of readout integrated circuit (ROIC) is one of the key requirements for hyperspectral applications.

    InstrumentCRISMPRISMAGF-5 AHSIMAJISHyspIRI
    Platform nameMROPRISMASAST3000JUICEHyspIRI
    countryUSAItalyChinaEuropeUSA
    Launched year20052018201820222023
    Spectral bands-249330-220
    Spectral range/µm0.36~3.90.4~2.50.4~2.50.4~5.70.38~2.5
    Spectral/nm6.55105~103~710
    Format640×4801 000×256512×5121 024×1 024-
    detectorsHgCdTeHgCdTeHgCdTe--
    Pixel pitch/µm27303015-
    ROIC input stageCTIACTIACTIACTIA-
    Readout noise/e-<100<150/<35060/275<170-

    Table 1. Comparison of specifications of CRISM, PRISMA, Gaofen-5(GF-5), MAJIS and HyspIRI

    Acquisition of hyperspectral data: (a) hyperspectral data cube; (b) schematic diagram of a hyperspectral imaging system

    Figure 1.Acquisition of hyperspectral data: (a) hyperspectral data cube; (b) schematic diagram of a hyperspectral imaging system

    Different ROIC input stages are designed depending on the wavelength of infrared radiation. Capacitive trans-impedance amplifier (CTIA) can achieve high sensitivity,high linearity,and high injection efficiency7-8. It’s commonly used for hyperspectral shortwave IRFPA ROIC input stage as shown in Table 1. However,its complex op-amp structure and reset switch not only occupy a large area,but also introduce thermal and KTC noise.

    To suppress the low-frequency noise (1/f noise) and KTC noise,the correlated double sampling (CDS) has been widely studied in the field of IRFPA ROIC. The principle of CDS is to store 1/f noise on one clock phase and then subtract it from subsequent clock phases. The utilization of CDS in the column pitch or common output stage of IRFPA are reported in Refs. [9-12]. In recent years,to pursue the objective of increasing the frame rate,the CDS is integrated in the pixel as presented in Refs. [13-15]. As for the noise analysis,a method to evaluate the signal and noise of imagers with CTIA-CDS ROICs is developed by Jerris F. Johnson16,and a cyclostationarity-based analytical model of charge amplification with CDS is proposed in Ref. [17]. The analytical study of CDS KTC noise is presented in Refs. [18-19]. These analyses provided designers with efficient solutions to noise assessment and ROIC optimization.

    However,to enhance the CDS noise immunity,the influence of time interval between the clamp and sample-and-hold (SH) processes of CDS needs analysis and verification.

    For this purpose,we first expound the noise mechanism of CDS with adjustable intervals (AICDS). Then,the noise power spectral density (PSD) function of combined CTIA-CDS is calculated. Next,to verify the theory,a low-noise pixel circuit is designed,which integrates low-noise CTIA and AICDS controlled by a timing generator. In section 3,the circuit performance is simulated and the noise electron number is calculated. Finally,we draw conclusions.

    1 Noise analysis

    This section introduces the noise mechanism of CDS and combined CTIA-CDS as functions of time interval between the clamp and SH.

    1.1 Noise mechanism of CDS

    The proposed CDS circuit is depicted in Fig. 2(a),which consists of the sample capacitor C0,reset switch S2,sample-hold switch S4,sample-hold capacitor Csh and its reset switchS3. The processed result is output to the column stage through the source-follower with a switch.

    Correlated double sampling: (a) structure of CDS; (b) timing diagram of CDS,equivalent circuits: (c) reset of C0 step at t2; (d) reset of Csh and clamp step at t3; (e) sample-hold step at t4

    Figure 2.Correlated double sampling: (a) structure of CDS; (b) timing diagram of CDS,equivalent circuits: (c) reset of C0 step at t2; (d) reset of Csh and clamp step at t3; (e) sample-hold step at t4

    Referring to the operating timing shown in Fig. 2(b),the equivalent circuits of different steps are shown in Figs. 2 (c),2(d),and 2(e). The input signal is first sampled on C0. Secondly,when S2 is off at t2,the right plate of C0 floats,clamping the charge stored in C0. Lastly,the redistribution of stored charge on C0 and Csh achieves when S4 turns on. As a consequence,the signals at different times are subtracted. The specific analysis of noise in a working cycle is as follows.

    We consider vn,i(t) as the noise of time t,whose RMS is the square root of the sum of the independent noise powers. At the moment t2,the noise charge stored on C0 is:

    Qn=C0vn,i(t2) .

    At the moment t4,the noise voltage input to the left pole plate of C0 becomes vn(t4). According to the principle of charge conservation,the charges at the two moments correspond to:

    C0vn,it2=C0vn,it4-vn,ot-Cshvn,ot.

    The output noise voltage is:

    vn,ot=C0C0+Cshvn,it4-vn,it2 .

    Applying the Laplace transform to Eq.(3)

    vn,os=C0C0+Cshvn,is1-e-sTd .

    Converting this to the frequency domain,the modulus of the transfer function is:

    HDω=C0C0+Csh21-cos ωTd=2C0C0+Cshsin ( ωTd/2) .

    The red curve in Fig. 3(a) is the amplitude-frequency plot of the flicker noise,represented as:

    1/f Noise Spectral Density reshaped by CDS with different Td : (a) the spectral density of 1/f noise; (b) comparison of two CDS transfer functions with different Td; (c) comparison of two output noise spectral density functions with different Td

    Figure 3.1/f Noise Spectral Density reshaped by CDS with different Td : (a) the spectral density of 1/f noise; (b) comparison of two CDS transfer functions with different Td; (c) comparison of two output noise spectral density functions with different Td

    Vn¯=KCoxWL1f,

    where K is a quantity related to the manufacturing process,Cox is the gate-oxide capacitance per unit area and WL is the channel area. Here,we assumeKCoxWL=1nV to simplify calculations.

    The transfer function of two different CDS is shown in Fig. 3(b),one is sin ( ωTd/2),the other is sin ( ω2Td/2). Here,we ignore the influence of capacitance and set Td=0.5 μs. The input noise is reshaped through CDS to the output noise presented in Fig. 3(c). By reducing Td to lengthen the transfer function's period,it can be determined that the noise power in the low-frequency section may be decreased.

    1.2 Output noise PSD of combined CTIA-CDS circuit

    The subtraction step of CDS can be considered as a delay-subtractor 20 as analyzed in the previous section. The CTIA,which is a switched-capacitor integrator,is equal to a first-order low-pass filter. The unit circuit structure model is shown in Fig. 4.

    Structure of combined CTIA-CDS ROIC: (a) the schematic structure; (b) equivalent model structure

    Figure 4.Structure of combined CTIA-CDS ROIC: (a) the schematic structure; (b) equivalent model structure

    The noise power density of first-order low-pass filter:

    S0ω=Sw1+ωτ02

    where the time constant τ0=RC.

    The noise power spectrum after the delayed subtractor is:

    SDω=S0ωHDω2 .

    It can be seen from Eqs. (5) and (8) that reducing Td can decrease the noise power in the low-frequency section.

    2 Proposed ROIC input stage

    To verify the analysis above,a new ROIC input stage composed of low noise CTIA and AICDS is designed.

    2.1 Combined CTIA-AICDS circuit

    The combined CTIA-AICDS circuit proposed is shown in Fig. 5. The detector is directly connected to CTIA,which contains a reset switch S1 and the feedback capacitor Cint. Then,the AICDS receives the integral signal. Finally,the modified signal is transmitted through SF to the column stage.

    Combined CTIA-AICDS ROIC

    Figure 5.Combined CTIA-AICDS ROIC

    The conventional design uses a differential pair as the operational amplifier,where the MOS occupies nearly twice the area. Its input reference thermal noise voltage is :

    Vn,in2¯=8kT23gm1+2gm33gm12 .

    This design uses a cascode type operational amplifier with an input reference thermal noise voltage of:

    Vn,in2¯=4kT23gm1+2gm33gm12 .

    It is close to half of the differential pair op-amp. This noise cannot be eliminated by the subsequent CDS structure 21,but will be multiplied through two sampling process,having a large impact on the circuit noise performance. Therefore,the cascode op-amp has benefits on noise reduction and area saving.

    Considering to further reduce the noise of CTIA,the following guidelines must be adhered to in the design.

    The operating current ID of the cascode must be low. The small operating current ensures a smaller power consumption,reducing the detector image interference caused by the heat of the pixel circuit. ID is determined by Vbias,which is provided through the current mirror in the bias circuit. The current mirror of the selectable channel generates ID around 100nA in this application.

    WL of input transistor M1 must be large whereas WL of load transistor M3 must be small.

    These two conditions make M1 work in the sub-threshold region,where the leakage current of M1 is:

    ID=I0expVGSζVT

    where ζ > 1 is a non-ideal factor and VT=kT/q.

    At this point,ID is exponentially related to VGS and the trans-conductance of MOS device M1 is:

    gm1=dIDdVGS=1ζVTID .

    The gain of the operational amplifier is

    AV=-gm11+gm2ro2ro1+ro2ro3 .

    It guarantees a large gain of M1 in the subthreshold region as we can conclude from Eqs. (12) and (13).

    The operating timing of the pixel circuit is shown in Fig. 6. Adjusting the additional CDS reset time Tr can regulate Td. At the beginning of each frame,S1 andS2 turn on to reset Cint and C0. After enough reset time,S1 is disconnected to start the integration process. Due to the voltage regulation characteristic of the op-amp,Vin remains unchanged and V0 is integrated as photocurrent is injected into the feedback capacitor. Subsequently,the transmission to CDS of V0 operates as described in section 2.

    Timing diagram and node voltage of the proposed ROIC pixel

    Figure 6.Timing diagram and node voltage of the proposed ROIC pixel

    2.2 Timing generator of AICDS

    The proposed timing generator can adjust the additional CDS reset time Tr relative to integration-reset time t1 as exhibited in Fig. 6. The CTIA integral reset signal intrst generates two pulse signals φ0 and φ1 to mark its rising and falling edges. As seen in Fig. 7,the division,counter,comparator,and latch make up the CDS timing control mechanism. It operates as follows.

    The proposed timing generator of AICDS

    Figure 7.The proposed timing generator of AICDS

    a. The high-frequency master clock CLK is first divided into a clock,whose cycle period is determined by the required minimum time step.

    b. The counter generates 4 bits Gray-code from the divided clock.

    c. Comparing Gray-codes with control bits TR[3:0],it generates a pulse signal φc to determine the falling edge of Tr when digital bits equal.

    d. φ1 and φc determine the operational time of the counter.

    e. φ0 and φc are injected to the latch to generate cdsrst.

    A time regulation step of 18 clock cycles is achieved by dividing the frequency with 9 DFFs. With 4 digital bits,a CDS reset time regulation ranging from 0 μs to 27 μs can be achieved for the master frequency CLK of 10 MHz,as shown in Table 2.

    TR3:0Clock cycles
    00000
    000118
    001036
    1111270

    Table 2. CDS reset time regulation

    3 Circuit performance and noise simulation

    In order to confirm the performance of circuit design,we have simulated the schematic using SPECTRE of Cadence. This circuit is designed in standard 180 nm CMOS process with 3.3 V supply voltage. Integral capacitanceCint=4 fF,CDS sample capacitor C0=200 fF and SH capacitor CSH=20 fF.

    3.1 Circuit performance simulation

    The performance of CTIA amplifier is simulated and shown in Table 3. Sweeping the injection photocurrent from 200 pA to 900 pA,the output waveform of V0 and V2 is presented in Fig. 8.

    Amplifier parameter
    gain59.6 dB
    Phase margin/°60
    Bandwidth/Hz8M
    Swing/V2.23(1-3.23)
    pixel pitch/μm15
    DC current/nA68~100

    Table 3. The performance of CTIA amplifier

    Integration and dynamic simulation

    Figure 8.Integration and dynamic simulation

    Figure 9 shows the simulation results for the CDS reset time regulation. Setting TCDSRST3:0=1001,the simulation waveforms demonstrate that Tr is 9 μs,which is consistent with intended specification.

    Time sequence transient simulation

    Figure 9.Time sequence transient simulation

    3.2 Noise simulation

    The common simulation method is to sum the noise of each separate cell. Since the CTIA reset noise and CDS of this circuit are correlated,the traditional simulation method cannot accurately reflect the impact of the interval on the noise. This design adopts transient noise simulation,which can accurately reflect the noise magnitude during the hold phase.

    We evaluate the equivalent noise charge (ENC) of input node to verify the noise performance of AICDS. The ENC μn is described as follows:

    μn=VneG

    where Vn is the noise voltage of output node and electron charge e=1.6×10-19 C. The conversion gain G is presented as:

    G=V2μe

    where V2 is the output voltage and μe is the electrons injected into the ROIC.

    In this design,μe is integrated from t2 to t4 as shown in Fig. 6.

    μe=IdetTd=IdetTint-Tr

    where Tint=34 μs.

    As for V2,we can calculate it from Eq. (3).

    V2=C0C0+CshIdetTdCint .

    Then the conversion gain G is:

    G=C0C0+Csh1Cint .

    For the case of Cdet=40 fF and photocurrent Idet=20 pA,24999 transient noise simulations are performed and sampled.

    For the results,the histogram of output voltage V2 under different settings of Tr is shown in Fig. 10. The results are finally brought into Eq. (14) to calculate the number of noise electrons.

    The histogram of V2 under different settings of Tr

    Figure 10.The histogram of V2 under different settings of Tr

    Noise electrons under different settings of Tr are shown in Fig. 11. The noise electron number of readout circuit ranges from 48.72 to 22.15,with Tr varying from 1 μs to 11 μs. After Tr=11 μs,the noise is stable at about 22 e- which is mainly decided by the post-stages.

    ENC and SNR vary with Tr

    Figure 11.ENC and SNR vary with Tr

    As Tr increases,the SNR can reach 42.98 dB at Tr=8 μs. However,after that,the signal quantity declines with decreased Td whereas the noise stays constant,resulting in attenuation of SNR. The simulation results show that extending Tr,in other words,decreasing Td,can suppress the noise of combined CTIA-CDS ROIC.

    4 Experimental results

    The full system is designed and fabricated in TSMC 0.18 μm CMOS technology. The fabricated chip is shown in Fig. 12,annotated with the floorplan. It consists of 640×512 array of 15 μm×15 μm pixel area. To meet the requirement of low noise,the pixel array is located at the center of the chip,while separated analog and digital blocks are placed at the edge.

    Layout and floorplan of the fabricated chip

    Figure 12.Layout and floorplan of the fabricated chip

    The platform designed to test the system is shown in Fig. 1322. LVDS module can generate the working time sequence and collect the output signal produced by the chip. The LabView-based test platform can control voltage,clock,and waveform. In addition,the acquired readout signal of the pixel array can be displayed as a grayscale map.

    Schematic diagram of experimental platform

    Figure 13.Schematic diagram of experimental platform

    To compare the noise performance of input stage with and without CDS,we designed a test schematic to turn off the CDS function as shown in Fig.14 (a). The timing diagram is presented in Fig.14 (b),where we can see S2 is invariably off,while S5 closes to be constantly controlled by a digital bit. As shown in Figs. 15(a) and 15(b),we get the grayscale map over 640×512 pixel array with CDS off and on. Intuitively,the uniformity with CDS on is far better than that with CDS off,which demonstrates that the noise can be suppressed by CDS. A closer look can be taken by acquiring partial 16×16 array to check the noise distribution,which is shown in Figs.15 (c) and 15(d). The noise with CDS on is far lower and much flatter than that with CDS off.

    Structure of CDS test: (a) the schematic structure;(b) timing diagram

    Figure 14.Structure of CDS test: (a) the schematic structure;(b) timing diagram

    The output voltage acquisition: (a) the grayscale map over 640×512 array when CDS-off;(b) the grayscale map over 640×512 array when CDS-on,the output noise acquisition;(c) the histogram of 16×16 segment when CDS-off;(d) the histogram of 16×16 segment when CDS-on

    Figure 15.The output voltage acquisition: (a) the grayscale map over 640×512 array when CDS-off;(b) the grayscale map over 640×512 array when CDS-on,the output noise acquisition;(c) the histogram of 16×16 segment when CDS-off;(d) the histogram of 16×16 segment when CDS-on

    The clock frequency is set at 500 kHz,which means the adjustment step is 36 μs. The additional CDS reset time Tr is adjusted and settled in condition of Tint=360 μs and Tint=640 μs. We measure the influence of Tr on the noise through the ratio of γ=Tr/Tint as presented in Fig. 16,where we can conclude that:

    The output noise voltage in condition of (a) CDS off and (b) varying ratio γ when CDS turns on

    Figure 16.The output noise voltage in condition of (a) CDS off and (b) varying ratio γ when CDS turns on

    a. The noise electron is up to 18.3 e- optimally by adjusting Tr,while the noise number is 39 in the case of turning off CDS. It shows an attenuation of 52%;

    b. With CDS on,the noise has a tendency of decreasing with the increment of Tr,whose slope is -24 µV.

    5 Conclusion

    In this paper,CDS with adjustable intervals has been proposed as a new technology to improve the noise performance of ROIC. It has been analyzed that decreasing the time interval between the clamp and SH can reduce the pre-stage noise. The feasibility of the technique is proved in a combined CTIA-AICDS controlled by a timing generator in 180 nm CMOS process. 24999 transient noise simulations demonstrate that by adjusting the interval,the ENC of input stage attenuates by 54%,and SNR is enhanced by 5dB. The experimental test results verify that the proposed structure can suppress the noise by 52% than CDS off and reduce the noise with a slope of -24 µV. It can be broadly used in Hyperspectral applications of low background radiation.

    References

    [1] J. Transon, R. d'Andrimont, A. Maugnard et al. Survey of Hyperspectral Earth Observation Applications from Space in the Sentinel-2 Context. Remote Sensing, 10, 157(2018).

    [2] S. K. Roy, S. Manna, T. Song et al. Attention-Based Adaptive SpectralSpatial Kernel ResNet for Hyperspectral Image Classification. Ieee Transactions on Geoscience and Remote Sensing, 59, 7831-7843(2021).

    [3] S. Murchie, R. Arvidson, P. Bedini et al. Compact reconnaissance Imaging Spectrometer for Mars (CRISM) on Mars Reconnaissance Orbiter (MRO). Journal of Geophysical Research-Planets, 112, 05-03(2007).

    [4] Yi-Bin BAI, William Tennant, Selmer Anglin et al. 4Kx4K Format, 10 mu m Pixel Pitch H4RG-10 Hybrid CMOS Silicon Visible Focal Plane Array for Space Astronomy, 8453(2012).

    [5] Xiao-Ning HU, Ai-Bo HUANG, Qing-Junet al LIAO. Large Format High SNR SWIR HgCdTe/Si FPA With Multiple-choice Gain for Hyperspectral Detection, 10213(2017).

    [6] Yin-Nian LIU, De-Xin SUN, Xiao-Ning HU et al. Development of visible and short-wave infrared hyperspectral imager onboard GF-5 satellite. Journal of Remote Sensing(Chinese), 24, 333-344(2020).

    [7] Peng-Yun SONG, Zhen-Hua YE, Ai-Bo HUANG et al. Theoretical investigation on input properties of DI and CTIA readout integrated circuit. Optical and Quantum Electronics, 48, 185(2016).

    [8] Mei ZOU, Ji-Qing ZHANG, Sheng-You ZHONG et al. Low-light-level CMOS imaging sensor with CTIA and digital correlated double sampling. Analog Integrated Circuits and Signal Processing, 101, 449-461(2019).

    [9] Chih-Cheng Hsieh, Chung-Yu Wu, Tai-Ping Sun. A new cryogenic CMOS readout structure for infrared focal plane array. Ieee Journal of Solid-State Circuits, 32, 1192-1199(1997).

    [10] S. W. Han, E. Yoon. Area-efficient correlated double sampling scheme with single sampling capacitor for CMOS image sensors. Electronics Letters, 42, 335-337(2006).

    [11] Xi-Qu CHEN, LYU Qiang. A versatile CMOS readout integrated circuit for microbolometric infrared focal plane arrays. Optik, 124, 4639-4641(2013).

    [12] Pan WANG, Rui-Jun DING, Zhen-Hua YE. High frequency weak signal analog chain design of short-wavelength IRFPAs. Infrared and Laser engineering, 43, 1370-1374(2014).

    [13] Y. S. Kim, D. H. Woo, Y. M. Jo et al. Low-Noise and Wide-Dynamic-Range ROIC With a Self-Selected Capacitor for SWIR Focal Plane Arrays. IEEE Sensors Journal, 17, 179-184(2017).

    [14] Zhang-Cheng HUANG, Song-Lei HUANG, Xu-Quan WANG et al. A 1024x512 ROIC with 30 mu m pixel pitch and 250Hz high frame rate for shortwave infrared detector, 10624, 2018(2018).

    [15] Yi ZHUO, Wen-Gao LU, Shan-Zhe YU et al. A Low-Noise CTIA-based pixel with CDS for SWIR Focal Plane Arrays, 258-262(2021).

    [16] J. F. Johnson, T. S. Lomheim. Focal-Plane Signal and Noise Model-CTIA ROIC. Ieee Transactions on Electron Devices, 56, 2506-2515(2009).

    [17] P. Pittet, G. N. Lu, L. Quiquerez et al. Cyclostationarity-based model for noise analysis of charge amplification with correlated double sampling. Analog Integrated Circuits and Signal Processing, 65, 407-413(2010).

    [18] S. Im, S.-G. Park. Thermal noise analysis of switched-capacitor integrators with correlated double sampling. International Journal of Circuit Theory and Applications, 44, 2101-2113(2016).

    [19] A. Caizzone, A. Boukhayma, C. Enz. An Accurate kTC Noise Analysis of CDS Circuits, 22-25(2018).

    [20] H. M. Wey, W. Guggenbuhl. Noise transfer characteristics of a correlated double sampling circuit. Ieee Transactions on Circuits and Systems, 33, 1028-1030(1986).

    [21] R. J. Kansy. Response of a correlated double sampling circuit to 1/f noise. Ieee Journal of Solid-State Circuits, 15, 373-375(1980).

    [22] Yi-Qun CEN, Jun-Ling ZHANG, Hong-Lei CHEN et al. Research on testing technology of column-level ADC in IRFPA digital readout circuits. Infrared and Laser Engineering, 49, 102-109(2020).

    Shuang WU, Qing-Hua LIANG, Hong-Lei CHEN, Rui-Jun DING. Low noise ROIC integrated with correlated double sampling with adjustable intervals for hyperspectral applications[J]. Journal of Infrared and Millimeter Waves, 2024, 43(2): 268
    Download Citation