Next Article in Journal
Real-Time Plant Leaf Counting Using Deep Object Detection Networks
Previous Article in Journal
Identification of Bridge Key Performance Indicators Using Survival Analysis for Future Network-Wide Structural Health Monitoring
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Toward the Super Temporal Resolution Image Sensor with a Germanium Photodiode for Visible Light

1
School of Science and Engineering, Ritsumeikan University, 1-1-1 Noji-Higashi, Kusatsu, Shiga 525-8577, Japan
2
School of Electronics and Telecommunications, Hanoi University of Science and Technology, 1 Dai Co Viet, Bach Khoa, Hai Ba Trung, Hanoi 100803, Vietnam
3
IMEC, 3001 Leuven, Belgium
4
Faculty of Information Science and Technology, Osaka Institute of Technology, Hirakata Campus, 1-79-1 Kitayama, Hirakata City, Osaka 573-0196, Japan
5
Link Research Corporation, 291-4, Kuno, Odawara-shi, Kanagawa 250-0055, Japan
6
Graduate School of Engineering, Osaka University, 1-1 Yamada-oka, Suita, Osaka 565-0871, Japan
7
Department of Civil and Environmental Engineering, School of Science and Engineering, Kindai University, 3-4-1 Kowakae, Higashi-Osaka, Osaka 577-8502, Japan
8
Advanced Quantum Architecture Laboratory, EPFL, Rue de la Maladière 71b, Case Postale 526, CH-2002 Neuchâtel, Switzerland
*
Author to whom correspondence should be addressed.
Submission received: 30 October 2020 / Revised: 19 November 2020 / Accepted: 24 November 2020 / Published: 2 December 2020
(This article belongs to the Special Issue Advances in High-Speed CMOS Image Sensor and Related Technologies)

Abstract

:
The theoretical temporal resolution limit t T of a silicon photodiode (Si PD) is 11.1 ps. We call “super temporal resolution” the temporal resolution that is shorter than that limit. To achieve this resolution, Germanium is selected as a candidate material for the photodiode (Ge PD) for visible light since the absorption coefficient of Ge for the wavelength is several tens of times higher than that of Si, allowing a very thin PD. On the other hand, the saturation drift velocity of electrons in Ge is about 2/3 of that in Si. The ratio suggests an ultra-short propagation time of electrons in the Ge PD. However, the diffusion coefficient of electrons in Ge is four times higher than that of Si. Therefore, Monte Carlo simulations were applied to analyze the temporal resolution of the Ge PD. The estimated theoretical temporal resolution limit is 0.26 ps, while the practical limit is 1.41 ps. To achieve a super temporal resolution better than 11.1 ps, the driver circuit must operate at least 100 GHz. It is thus proposed to develop, at first, a short-wavelength infrared (SWIR) ultra-high-speed image sensor with a thicker and wider Ge PD, and then gradually decrease the size along with the progress of the driver circuits.

1. Introduction

The theoretical temporal resolution limit of a silicon photodiode (Si PD) is 11.1 ps (90.1 Gfps) [1]. However, users in advanced sciences require much faster multiframing than this limit. For example, a pump-probe method is a common technique to measure phenomena changing in several femtoseconds to nanoseconds. If the change is captured with one-shot imaging, the measurement performance is significantly improved. Furthermore, the pump-probe method can be applied to only reversible phenomena. This paper is a pilot study toward such an ultrafast continuous imaging. We propose the definition of “super temporal resolution (STR)” as a temporal resolution less than the theoretical limit of the Si PD since silicon image sensors are commonly used for high-speed imaging.
The triangles in Figure 1 show imaging devices having achieved STR [2,3,4,5,6,7,8,9,10,11,12,13,14]. Some of which are presented with motion pictures on the page for Invited Speakers on the webpage of ICHSIP31 [15].
In 1991, the frame rate of the fastest image sensor in the world was 4500 fps [16], have exponentially increased [17,18,19,20,21], and, in 2018, it reached 100 Mfps (equivalent to the frame interval of 10 ns) [22,23,24,25,26,27]. For X-ray imaging, a temporal resolution of 3 ns (333 Mfps) has been achieved using a silicon image sensor [23,27]. The large number of electrons generated by a single X-ray photon allows an increase of the frame rate without the suppression of the noise associated with the conversion of electron packets to voltage signals at an extremely high operating rate. For visible light, the current highest frame rate is about 100 Mfps.
Numerous devices have been developed for high-speed imaging. However, once an image sensor has achieved a high frame rate, other devices developed for imaging for a similar frame rate were immediately replaced by a camera with the image sensor. Except silicon image sensors, there still remain the imaging devices which achieved the STR. The frame rate of about 50 ps may be achieved by an image sensor with existing technologies [25]. This is close to the theoretical limit. Now, it is worth exploring image sensor technologies to go beyond the theoretical limit of the Si PD.
Germanium is a promising material for a high-speed transistor and for a photodiode in an SWIR range [28,29,30,31,32,33,34]. Thermal imaging together with motion imaging is very important in ultra-high-speed imaging. For example, in the automobile industry, high-speed imaging of heat generation provides useful information for improving the combustion efficiency of engines, tire wheels at sudden braking, a blast of an airbag, and laser and plasma processes in machining. However, the highest frame rate of existing IR video cameras is only 10 kfps for a reasonable pixel count. The nominal highest frame rates in the catalogues are about 100 kfps. However, the pixel count for the frame rate is impractically small.
STR imaging with a Ge PD for visible light arose in our discussion on the development of a high-speed SWIR image sensor. For a Si PD, the temporal resolution limit is approximately proportional to the ratio of the penetration depth δ of light to the saturation drift velocity v c [1]. The penetration depth of visible light to Ge is 1/50 to 1/90 of that to Si with a large variation depending on the source literature, and the saturation drift velocity v c in Ge is about 2/3 of that in Si. Therefore, it is expected that the temporal resolution limit of the Ge PD for the visible light is less than 1/30 ( 1/50 × 3/2) of the theoretical limit of the Si PD.
This article is organized as follows. Definitions on the temporal resolution are summarized. Tools for the analyses are explained, including the exact formulation based on the drift-diffusion equation and the flux equation (the DDF model), the approximate explicit expression derived by us [1] with their asymptotic expressions for the infinitesimal thickness of the photodiodes, and the Monte Carlo simulation codes. Models and the parameters for the analyses are summarized. Then, the temporal resolutions of the Si PD and the Ge PD for visible light are analyzed. The Si PD has been analyzed by our inhouse simulation code for Si semiconductors. For the analysis of the Ge PD, a Sentaurus MC simulation code by Synopsys is introduced. However, standard deviations of the propagation times of charges between electrode pairs are not included in the output list, which are the basis to evaluate the temporal resolution. Therefore, a method to estimate the temporal resolution with the Sentaurus code was proposed. To keep consistency of the analyses, both simulation codes are applied to a very thin Si PD to confirm agreement of the results. Then, the temporal resolution of the Ge PD for visible light are evaluated with the proposed method.

2. Definitions of Temporal Resolution and Tools for Analyses

2.1. Criteria for Temporal Resolution

The spread of the arrival time of signal electrons to a detection plane is represented by the standard deviation σ . As shown in Figure 2, the sum of the ordinates of two Gaussian distributions has no separable peak for μ 2   μ 1 < 2 σ , where ( μ 2 > μ 1 ) , and μ 2   and     μ 1 are the averages of the Gaussian distributions. Physically, the time difference of a pair of incident light pulses to the backside of an image sensor is kept as the difference of the averages of the arrival times as shown in Figure 3. Therefore, we defined the theoretical temporal resolution limit as t T = 2 σ . However, 16% of the distribution function overlaps the next one. Therefore, we also propose another criterion to define the practical temporal resolution t P = t 95 t 5 where t 95   and   t 5 are, respectively, the times when 95% and 5% of the signal electrons arrive at the detection plane. For the Gaussian distribution, t P = 3.30 σ . Empirically, this relationship fits well even for a skewed distribution of the arrival time with a long tail.
To evaluate the temporal resolution, the thickness of the photodiode needs to be specified. For the theoretical analysis, we employed a thickness W T = δ , where δ is the average penetration depth, 1/(absorption rate α ), of the incident light. For this condition, 63.2% of incident photons are converted to electron-hole pairs, and 36.8% are passing through the photodiode. Ultra-high-speed imaging requires the highest possible sensitivity. Therefore, we also use another thickness, W P = 3 δ , for practical evaluation for which 95.0% of the incident photons undergo photoelectron conversion.
Therefore, in this article, the temporal resolution is evaluated for a combination of the thickness of the photodiode and the standard deviation of the arrival time; W T ,   t T = δ ,   2 σ and W P ,   t P = 3 δ ,   3.3 σ for the theoretical and the practical evaluations.

2.2. Temporal Resolution Based on DDF Models

The standard deviation of the arrival time can be exactly formulated based on the drift-diffusion equation and the flux equation. Hereafter, the model expressed by the set of the equations is called the DDF model. For a simple model as shown in Figure 4, the standard deviation is expressed with a set of equations shown in Appendix A. The expression requires numerical integration to obtain the standard deviation.
Etoh et al. derived the approximate explicit expression which accurately agrees with the numerical result of the DDF model except for a very thin photodiode [1]. The approximate expression is shown in Appendix B.
The thickness of the Ge PD for visible light is extremely thin due to much higher absorption coefficient of Ge than that of Si. Therefore, we also derived the asymptotic expressions for the infinitesimal thickness W of the photodiode. For the approximate expression, the asymptotic expressions can be derived as shown in Appendix C. Analytical derivation of the limit may be possible for the infinitesimal W . For the DDF model, we found the asymptotic expression through numerical integration of the exact formulation for various parameter values. The standard deviation always converges to 5 D c / v c 2 , which is physically incorrect. It may be due to the expression formulated with continuous variables, while the stepwise motion of the electrons becomes significant in such a thin layer.

2.3. Monte Carlo Simulation Codes

The Monte Carlo (MC) simulation code for Si used in this paper is an inhouse one developed by Kunikiyo, Kamakura et al. in Osaka University for homogeneous intrinsic silicon [35]. The model employs realistic physical models for (1) the energy band structure, (2) the phonon dispersion relations, (3) the electron-phonon scattering rate, (4) the impact ionization rate, and (5) the conservation of both energy and momentum in each scattering process to properly represent the final state of the scattering electron. The model requires no adjustable parameters and the results quantitively agree to the experimental data.
We do not have an MC simulation code for Ge. Therefore, we employed the Sentaurus MC simulation code marketed by Synopsys. The code is applicable to a Si1−xGex compound, where x denotes the ratio of Ge. However, it is not directly applicable for our analysis. For example, the standard deviation to estimate the temporal resolution is not included in the output list, and the source of electrons cannot be placed in the bulk. The program can be easily modified, but the source program is not open. Therefore, we developed a method to utilize the code for our purpose without modification of the code (Appendix D).

3. Numerical Analyses

3.1. Model Parameters

Figure 5 shows the dependence of the drift velocity and the longitudinal diffusion coefficient on the electric field for Si and Ge [36,37,38]. The temperature is mostly 300 K. We defined the critical field as the field at which the drift velocity is 95% of the saturation value, since the saturation value keeps elongated at a constant value for a wide range of fields above their critical values (truncated in Figure 5). Table 1 shows for Si and Ge, (1) critical field, (2) drift velocity and diffusion coefficient for the critical field, and (3) penetration depth of light of wavelengths of 550 nm and 1 μ m [39,40].
Average penetration depths (absorption coefficients) largely vary, depending on the source in the literature. For example, the penetration depth of light at 550 nm in Si can be 1.08, 1.44, and 1.73 μ m , depending on the references [39,40,41]. We employed the value 1.73 μ m from the latest reference [39], providing the largest standard deviation and the safety-side estimation of the temporal resolution limit.
For Ge, some sources are available for the absorption coefficients in the SWIR range. However, we found only one data set for visible light in [40] (pp. 465–478) published in 1985. The absorption coefficient based on the data set is easily extracted by the code supplied by [42]. Figure 6 seems to be generated from the same source and copied several times by many people. If we calculate the penetration depth of Ge for 550 nm directly from [40] or [42], the value is 20.0 nm, which is employed in this paper. However, from Figure 6, it is 26.7 nm, which may be due to rounding during multiple copying.
Considering the distributions of data on the penetration depths in the literature, the ratio of the penetration depth of Ge to Si varies from about 1/50 to 1/90. Therefore, the values presented in the conclusions of this paper are the tentative values. Accurate data for the penetration depths are expected to be specified by a professional society.
From Table 1, the ratio of the penetration depth for Ge and Si for light at 550 nm is 1/86.5 (=0.020/1.73). The critical field of Ge, 4.25 kV/cm, is about 17% of that of Si, 25 kV/cm. The saturation velocities and the dielectric constants of Ge are, respectively, about 2/3 and 3/2 of those of Si. The diffusion coefficient for Ge is about four times larger than that of Si.

3.2. Analysis for Si PD

In Figure 7, the green line σ d d f shows the standard deviation of the arrival time calculated by numerical integration of the expression for the DDF model shown in Appendix A with the parameters shown in Table 1. The orange dotted line σ a p p shows the approximate expression of the standard deviation shown in Appendix B with the mixing component σ m i x , that is due to the distribution of the penetration depth of incident light and the drift velocity, and the diffusion component σ d i f f . The figure suggests that:
  • there are three ranges: the mixing (drift) effect is dominant for 0.14   μ m < W < 100   μ m , which covers the whole practically meaningful range, and the diffusion effect is dominant for the very thin range W < 0.14   μ m and for the very thick range W > 120   μ m   (not drawn);
  • for W > 0.14   μ m , the approximate expressions almost perfectly agree with the exact solution from the DDF model; for W < 0.14   μ m , the discrepancy increases and the DDF model asymptotically converges to a constant value, σ d d f L = 5 D c / v c 2 ;
  • the range used in practical applications, δ < W < 3   δ , is included in the drift-dominant range as concluded in our previous paper [1];
  • The theoretical temporal resolution limit for W T ,   t T is δ ,   2 σ = 1.73   μ m ,   11.1   ps , and the practical limit for W P ,   t P is 3 δ ,   3.3 σ =   5.19   μ m ,   45.2   ps .
Since the penetration depth of visible light is more than 1 μ m , except for light close to the UV region, our approximate expression can be used for estimation of the temporal resolution limit for the thickness range. However, the penetration depth of 550 nm light to Ge is 20.0 nm. For W < 1 μ m , these values significantly deviate from each other. Therefore, we evaluated the standard deviation of the arrival time by using Monte Carlo simulations. For W   > 1 μ m , the MC simulation provides the standard deviation σ M C very close to σ d d f from the DDF model and to σ a p p from our approximate expression [1].
Previously, we have used the inhouse MC simulation code for Si. However, we do not have the code for Ge. Therefore, the Sentaurus MC simulation code by Synopsys [44] is applied to the analysis on the Ge PD as shown in the next subsection.
The penetration depth δ of 550 nm light in a Ge PD is 20 nm, and 3 δ is 60 nm. Therefore, the temporal resolution of Ge PD is analyzed later for the thickness range including these values. The consistency of the analyses in this range by using these two different simulation codes are confirmed for a Si PD. The standard deviations σ M C and σ M C S of the arrival time for 10 nm < W < 100 nm are evaluated and compared, where σ M C and σ M C S are the standard deviations calculated with our inhouse code and by Sentaurus (with post-processing described in Appendix D). Figure 8 shows the comparison result.
The standard deviations σ M C and σ M C S agree well around W = 100   nm . However, they show a difference of 28% at W = 10   nm .
Some other interesting facts are also observed as follows:
  • the Monte Carlo simulation shows that σ M C monotonically decreases when the thickness of the photodiode decreases, while the numerical solution σ d d f of the DDF model converges to a constant value for an infinitesimal thickness;
  • for 0.1   μ m < W < 0.5   μ m , σ M C departs from σ d d f , and goes along the approximate expression σ a p p ;
  • for W < 0.1   μ m , σ M C   resulting from our inhouse MC simulation code departs from σ a p p , and, for W < 0.03   μ m , it becomes parallel to the mixing (drift) component σ m i x L = 1 / 12   W / v c with the slope proportional to W , suggesting that the motion of the signal electrons converges to a ballistic motion (see Appendix C).

3.3. Super Temporal Resolution Limit for a Ge PD

We used Sentaurus for the analysis of the Ge PD. A method is developed to estimate the standard deviation of the arrival time from standard outputs from the code (Appendix D).
Figure 9 shows that the theoretical and the practical temporal resolution limits for a Ge PD receiving incident light of 550 nm are, respectively, 0.26 and 1.41 ps for ( W T , σ ) = (20 nm, 0.13 ps) and ( W P ,   3.3 σ ) = (60 nm, 0.43 ps). The accuracy may be within 30%, considering the result for Si.
The Sentaurus MC simulation of σ M C S appears to approach σ m i x L   = 1 / 12   W / v c , representing the mixing effect due to the drift motion and the distribution of the penetration depth. However, in Figure 8 for Si, the standard deviations calculated with our code and the Sentaurus code seem parallel at W = 10   nm , showing 28% difference. Further research will make it clear which is correct, asymptotic or parallel approach to σ m i x L . Based on the result, it may be possible to propose an approximate expression of the standard deviation for W < 1   μ m which is smoothly connected to our approximate expression σ a p p for W > 1   μ m .

3.4. Toward Super Temporal Resolution through SWIR Imaging

At this moment, it is not practical to try to directly achieve the theoretical temporal resolution limit of Ge PDs. The highest f T of transistors is higher than 100 GHz (equivalent to 10 ps). Germanium is a popular material for a photodiode for SWIR imaging up to 1.5 μ m light. With some strain to the crystal, the sensitivity can be increased to 1.55 μ m at which an optical fiber for communication takes the minimum absorption coefficient. Therefore, we will propose to develop an ultra-high-speed SWIR image sensor of the temporal resolution about 100 ps at first, targeting the STR of 10 ps by gradually thinning the thickness of the Ge PD following the progress of semiconductor technology. It may take a decade as shown in Figure 1. Therefore, the temporal resolution of the Ge PD for the longer wavelength is also analyzed.
Figure 10a,b shows standard deviation results for the Ge PD when incident light with the wavelengths λ of 1 and 1.5 μ m. The practical temporal resolutions for the wavelengths are 61.5 ps for 1.58 μ m and 180.8 ps for 6.12 μ m.
The result suggests development of a SWIR ultra-high-speed image sensor with the thickness of the Ge PD of 1 5   μ m at first, which has many important practical applications.
In the SWIR range, the drift motion effect gradually increases. For λ = 1   μ m in Figure 10a, the diffusion and the drift components are comparable, which makes a slight hump in the shape of the temporal resolution. For λ = 1.5   μ m in Figure 10b, the drift component is dominant between δ < W < 3 δ . In these ranges, our approximate expression accurately follows the DDF model, even including the humps. In the range, MC simulations are not applied, since the results may fit our approximate expression as experienced in the simulations for a Si PD.
For the thickness of the photodiode W > 1   μ m , the temporal resolution limit is accurately calculated by our simple approximate expression. However, for W < 1   μ m , it is expected for someone to derive a new expression of the standard deviation based on the physics including the ballistic motion, and other important physical components.

4. Dark Current

Figure 11 shows a pixel of a backside-illuminated burst image sensors with a Si charge collection pyramid that we previously proposed [45], and with a thin Ge PD. The upper layer of the former structure is a Si PD of about 10   μ m thick and the lower layer is a silicon diffusion circuit layer of about 1   μ m thick or less. The pyramid efficiently works to collect signal electrons in a very short time with a 100% fill factor. In Figure 11b, the Si PD is replaced by a Ge PD of about 0.1   μ m thick. To suppress spread of the arrival time for STR imaging, the lateral size of the backside surface of the Ge PD is about 200 nm, which requires an advanced semiconductor process and an on-chip microlens array.
A Ge PD generates a huge dark current in SWIR imaging at a video rate of 30 fps. However, the dark current during a frame interval of 10 ps for STR imaging is 1/33,333,333,333 of that for video-rate SWIR imaging. Thickness of a Ge PD for STR imaging is about 1/100 of that for SWIR imaging, further reducing the dark current to 1/100. For the same thickness, the dark current is proportional to the square of the lateral size. When the pixel size is about 1   μ m , the dark current further reduces to 1/25. Therefore, deep cooling may be unnecessary for the STR image sensor with a Ge PD. Instead, the most important design issue is to efficiently suppress migration of the dark current to the Si circuit layer during transfer of voltage signals converted from signal electron packets in floating diffusions to an analogue or digital memory unit in a functional chip stacked to the sensor chip.

5. Concluding Remarks

5.1. Conclusions

Super temporal resolution (STR) is defined as a resolution better than the theoretical temporal limit 11.1 ps for silicon. This paper reports a pilot study to achieve STR with an image sensor. Germanium is employed as a candidate photodiode material in the visible light range to show the fundamental feasibility.
Previously, we defined the theoretical temporal resolution limit as t T = 2 σ   for   W T = δ , where σ is the standard deviation of the arrival time of signal electrons,   δ is the average penetration depth of incident light, and W T is the thickness of the photodiode [1]. However, the condition results in 16% of temporal crosstalk and 37% loss of incident light passing through the photodiode. Therefore, the practical temporal resolution limit is defined in this paper as t P = 3.3 σ for W P = 3 δ . For this condition, both the temporal crosstalk and passing-through electrons are both 5%. Then, for silicon, t T = 11.1   ps and t P = 45.2   ps .
Monte Carlo simulations are applied to an intrinsic Ge PD to estimate these values. The theoretical and the practical temporal resolution limits of the Ge PD are estimated as 0.26 and 1.41 ps, respectively.
For W > 1   μ m , our simple approximate expression of the standard deviation of the arrival time agrees almost perfectly to the MC simulation results and to the expression based on exact formulation, though, requiring numerical integration. For W < 1   μ m , a simple approximate expression is expected to be derived.
The dark current is a major issue for SWIR imaging with a Ge PD. However, it is not a serious problem during an image capturing operation with a STR image sensor. Instead, it is a major design issue to efficiently suppress migration of the dark current to the circuit layer after the image capturing operation.

5.2. Remarks

The proposed strategy to estimate the temporal resolution limits may be acceptable as a reasonable approximation. However, the values of the temporal resolution limits should be reexamined based on the accurate measurement of the average penetration depths of light which widely vary depending in the literature. We asked Synopsys to improve the Sentaurus MC simulation code for directly calculating standard deviations of the arrival times, which is not supported by the current version. The new version will appear in the near future.
The germanium process has not been fully developed for practical applications. The driver circuit at this moment cannot operate at such a high-speed as the STR image sensor. It is proposed to develop a SWIR ultra-high-speed image sensor with a wider Ge PD of 1 μ m thick for NIR light or 5 μ m thick for 1.5 μ m light as shown in Figure 10, and decrease the size toward the STR along with progress of relating technologies. The SWIR ultra-high-speed image sensor itself has many important practical applications.
It may take a decade to achieve the STR image sensor as predicted in the trend line in Figure 1. Further research with other materials and structures is encouraged.

Author Contributions

Contributions of the authors are as follows: N.H.N.: Numerical and theoretical analyses, and writing; A.Q.N.: Device and Monte Carlo simulations for silicon; F.M.B.: Advice on the Sentaurus MC simulation code; Y.K.: Development of the silicon MC simulation code; H.M.: Development of the device simulator; P.M.: Communication to IMEC; T.H., T.S. and H.W.: Supply of information on a Ge photodiode; K.S.: A project organizer; K.T.: A project organizer; E.C.: Driver circuit and Ge process; T.G.E.: Proposal of the concept, theoretical analyses and writing. All authors have read and agreed to the published version of the manuscript.

Funding

This research was financially supported by the Japanese government; the Grant-in-Aids for Scientific Research: 18H01548(B) (Rep.: K. Takehara), and 19H02204(B) (Rep.: K. Shimonomura). The full-band Monte Carlo simulation code for silicon was developed and provided by T. Kunikiyo and Y. Kamakura [35].

Conflicts of Interest

The authors declare no conflict of interest.

List of Symbols

SymbolDescription
α Absorption rate of incident light
δ Penetration depth of incident light
v c Drift velocity at saturation
D c Diffusion coefficient at the drift velocity saturated
μ Mean
t T Theoretical temporal resolution limit
t P Practical temporal resolution
σ Standard deviation of the arrival time of signal electrons
σ d d f Standard deviation of the drift-diffusion-flux model
σ d d f L Limit of the standard deviation of the drift-diffusion-flux model
σ a p p Standard deviation of the approximation model
σ m i x Standard deviation of the mixing component
σ d i f f Standard deviation of the diffusion component
σ a p p L σ a p p for the infinitesimal thickness of the photodiode
σ m i x L σ m i x for the infinitesimal thickness of the photodiode
σ d i f f L σ d d f for the infinitesimal thickness of the photodiode
σ M C Standard deviation calculated with the inhouse MC simulation code
σ M C S Standard deviation calculated with the Sentaurus MC simulation code
W Thickness of the sensor
W T W for the theoretical temporal resolution limit
W P W for the practical temporal resolution limit

Appendix A. Exact Formulation of the Temporal Resolution

The expression of the temporal resolution is exactly formulated by the drift-diffusion equation and the flux equation. Numerical integration is necessary.
Table A1. Exact formulated expressions of the temporal resolution.
Table A1. Exact formulated expressions of the temporal resolution.
Strict Expression Asymptotic   Expressions   for   W     0
The Gaussian drift-diffusion equation for a single pulse g z ,   t   = 1 / 4 π D t × exp z v t 2 / 4 D t
The flux passing a detection plane h z , t   = v g z , t D g z , t / z = 1 / 2 v + z / t g z , t
The penetration depth distribution k s   =   1 / δ e s / δ
Convolution of k s and h z , t f t   = 0 W k s h W s ,   t d s
The 0th moment (Absorption rate) p = 0 f t d t W = W / δ
The 1st moment E t   = 1 / p × 0 t f t   d t D / 2 = D / v 2
The 2nd moment E t 2   = 1 / p × 0 t 2 f t   d t 3 / 2 D 2 = 6 D / v 2 2
Variance σ d d f 2 = E t 2 E t 2 σ d d f L   2 = 5 / 4 D 2 = 5 D / v 2 2
Temporal resolution Δ t d d f = 2 σ d d f 5 D = 2 5 D / v 2
Normalized parameters W = W / δ ,   t = δ / v ,   D = 2 D / v 2

Appendix B. An Explicit Approximate Expression for the Temporal Resolution

Table A2 shows an approximate expression derived by Etoh et al. [1].
Table A2. Approximate expressions of the temporal resolution.
Table A2. Approximate expressions of the temporal resolution.
Approximation FormulaSolutionsAsymptotic Expressions
for   W     0
The penetration depth distribution k s   =   1 / δ e s / δ
Absorption rate p = 0 W k s   d s 1 e W W = W / δ
Average arrival time t r =   W s / v
E t r 1 p 0 t r   k s   d s t W / p 1
E t r 2 1 p 0 t r 2   k s   d s t 2 W 2 2 W / p + 2
Mixing component of variance σ m i x 2 = E t r 2 E t r 2 t 2 1 W 2 × 1 p / p 2 W t 2 / 12 = W 2 / 12 v 2
Diffusion component of variance σ d i f f 2 = 2 D / v 2   E t r D t W p / p D t W / 2 = D W / v 3
Variance σ a p p 2 = σ m i x 2 + σ d i f f 2 D t W / 2 = D W / v 3
Approximation for a whole thickness σ a p p S 2 = σ a p p 2 + σ d d f L   2 σ a p p S 2 perfectly fits σ d d f 2
( σ d d f 2 and σ d d f L   2 are in Appendix A)
Temporal resolution Δ t a p p = 2 σ a p p
Normalized parameters W = W / δ ,   t = δ / v ,   D = 2 D / v 2

Appendix C. Derivations of the Asymptotic Expressions

The Taylor series of the absorption rate for the infinitesimal W :
p = 1 e W = 1 n = 0 f n 0 n ! W 0 n = W W 2 2 + W 3 6
The limit of the variance of the mixing component σ m 2 .
The Taylor series of p is expanded to the polynomial of the degree 3:
σ m i x 2 = t 2 1 W 2 × 1 p p 2 = t 2 1 W 2 × 1 W + W 2 / 2 W 3 / 6 W W 2 / 2 + W 3 / 6   2
= t 2 W 2 12 × 1 + W 2 / 3 1 W / 2 + W 2 / 6   2
lim W 0 σ m i x 2 t 2 W 2 12 × lim W 0 1 + W 2 / 3 1 W / 2 + W 2 6 2 = t 2 W 2 12 = 1 12 ( W v 2 )
The limit of the variance of diffusion component σ d 2 .
Similarly, the Taylor series of p is expanded to the polynomial of the degree 2:
σ d i f f 2 = D t W p p = D t × W W + W 2 / 2 W W 2 / 2 = D t W 2 × 1 1 W / 2
lim W 0 σ d i f f 2 D t W 2 × lim W 0 1 1 W / 2 = D t W 2

Appendix D. A Method to Estimate the Standard Deviation by Using the Sentaurus MC Simulation Code

There are some free software programs for MC simulations of Ge. However, as user’s guides are not available or insufficient for these programs, we chose Sentaurus. However, the code does not output even a standard deviation of the arrival time. Therefore, we developed a method to approximately estimate the standard deviation of the arrival time from standard outputs from the code.
In Sentaurus MC simulation, when a generated particle reaches an electrode, the next particle is generated at one of the electrodes one after another. The propagation time of each particle randomly distributes. The generation probability of each electrode is proportional to the product of the length of the electrode and the electron density at the adjacent element to the electrode. Therefore, the electrode from which the next electron is generated is also randomly selected. The propagation time of each particle is accumulated until the total time reaches the prespecified time interval t M for each calculation step of the averages.
Bufler proposed a special treatment for a small value of t M for our application. The currents between electrode pairs are calculated based on the particle counting. The current at an electrode is calculated from the product of the unit charge q and N k , the net number of the particles during the k-th calculation step. Although N k is not output from the simulation code, it can be calculated from the value of the current. The current I D (A/ μ m ) is calculated by the following formula:
I D = ( q × N k × w ) / t M
where q = 1.6 × 10−19 C is the unit electron charge, N k n is the number of electrons that have propagated to the front side, and w ( μ m 2 / cm 3 ) is a weight calculated by integration of the electron density over the device. Therefore, the average propagation time t k for the calculation step k is:
t k = t M / N k = I D / q × w
There are two problems:
(1)
N k is the number of all generated electrons, including electrons taking paths which are not our target, especially, the number of electrons absorbed at the source electrode and injected from another electrode,
(2)
N k is a random number.
To solve the first problem (1), we introduced a technique that a thin p-layer is placed in front of each electrode to reduce the electron density to suppress the generation probability of electrons, except the source electrode.
To handle the second problem (2), Etoh proposed an approximate formula to estimate the original standard deviation σ 0 from a sequence of averages t k , each calculated from N k samples, where k = 1,..., K , and N k randomly distributes.
Hereafter, t k is replaced by x k as a standard notation in statistics. Let μ denote the population mean (original average). The standard deviation σ k of the k-th group with N k samples is σ 0 / N k . Therefore, the standardized value X k   =   x k μ / σ k = x k μ × N k / σ 0 . The mean and the variance of a standardized value X k , respectively, distribute around 0 and 1. Therefore, it is expected for a sufficiently large K that the sum of the squares of X k approaches K , then,
k = 1 K X k 2 = k = 1 K x k μ 2   N k / σ 0 2 = K .
Therefore,
σ 0 2 = 1 / K × k = 1 K x k μ 2   N k
Figure A1. The model to estimate the original standard deviation σ 0 by using the Sentaurus MC simulation code (W = 100, 50, and 10 nm; the front-side doping for the Ge PD is 1015 cm−3).
Figure A1. The model to estimate the original standard deviation σ 0 by using the Sentaurus MC simulation code (W = 100, 50, and 10 nm; the front-side doping for the Ge PD is 1015 cm−3).
Sensors 20 06895 g0a1
The voltages of the backside, front-side, and source electrodes are adjusted to make the field approximately equal to the saturation field shown in Table 1. For a given s, the number of the electrons taking the path A in the figure is calculated from the current by assuming the propagation times of the electrons taking E, F, and G are negligibly small, and the generation probabilities of electrons taking C and D are suppressed, respectively, by the strong backside p-doping, and the adverse field and the front-side p-doping. The electrons taking B can be an error source, but the generation probability may not be so large, especially when s increases. Then, the average propagation time for s is calculated by dividing the time for the calculation step by the number of the electrons taking A. The simulation is repeated, and the variance is estimated by using Equation (A4) for a value of s.
Then, the variances are calculated for s = 0%, 10%, 20%, …, 100% of W, and averaged by weighting the probability of the penetration depth of light to estimate the variance of the propagation time considering the penetration depth distribution. The standard deviation is calculated as the square root of the variance to estimate the temporal resolution.

References

  1. Etoh, T.G.; Nguyen, A.Q.; Kamakura, Y.; Shimonomura, K.; Le, T.Y.; Mori, N. The theoretical highest frame rate of silicon image sensors. Sensors 2017, 17, 483. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Shiraga, H.; Miyanaga, N.; Heya, M.; Nakasuji, M.; Aoki, Y.; Azechi, H.; Yamanaka, T.; Mima, K. Ultrafast two-dimensional x-ray imaging with x-ray streak cameras for laser fusion research. Rev. Sci. Instrum. 1997, 68, 745–749. [Google Scholar] [CrossRef]
  3. Kubota, T.; Awatsuji, Y. Observation of light propagation by holography with a picosecond pulsed laser. Opt. Lett. 2002, 27, 815–817. [Google Scholar] [CrossRef] [PubMed]
  4. Kubota, T.; Komai, K.; Yamagiwa, M.; Awatsuji, Y. Moving picture recording and observation of three-dimensional image of femtosecond light pulse propagation. Opt. Express 2007, 15, 14348–14354. [Google Scholar] [CrossRef]
  5. Kakue, T.; Tosa, K.; Yuasa, J.; Tahara, T.; Awatsuji, Y.; Nishio, K.; Ura, S.; Kubota, T. Digital Light-in-Flight Recording by Holography by Use of a Femtosecond Pulsed Laser. IEEE J. Sel. Top. Quantum Electron. 2012, 18, 479–485. [Google Scholar] [CrossRef]
  6. Nagel, S.R.; Hilsabeck, T.J.; Bell, P.M.; Bradley, D.K.; Ayers, M.J.; Barrios, M.A.; Felker, B.; Smith, R.F.; Collins, G.W.; Jones, O.S.; et al. Dilation X-ray imager a new/faster gated X-ray imager for the NIF. Rev. Sci. Instrum. 2012, 83, 10E116. [Google Scholar] [CrossRef]
  7. Nakagawa, K.; Iwasaki, A.; Oishi, Y.; Horisaki, R.; Tsukamoto, A.; Nakamura, A.; Hirosawa, K.; Liao, H.; Ushida, T.; Goda, K.; et al. Sequentially timed all-optical mapping photography (STAMP). Nat. Photonics 2014, 8, 695–700. [Google Scholar] [CrossRef]
  8. Gao, L.; Liang, J.; Li, C.; Wang, L.V. Single-shot compressed ultrafast photography at one hundred billion frames per second. Nature 2014, 516, 74–77. [Google Scholar] [CrossRef]
  9. Cai, H.; Zhao, X.; Liu, J.; Xie, W.; Bai, Y.; Lei, Y.; Liao, Y.; Niu, H. Dilation framing camera with 4 ps resolution. Apl Photonics 2016, 1, 16101. [Google Scholar] [CrossRef] [Green Version]
  10. Suzuki, T.; Hida, R.; Yamaguchi, Y.; Nakagawa, K.; Saiki, T.; Kannari, F. Single-shot 25-frame burst imaging of ultrafast phase transition of Ge2Sb2Te5 with a sub-picosecond resolution. Appl. Phys. Express 2017, 10, 92502. [Google Scholar] [CrossRef]
  11. Liang, J.; Ma, C.; Zhu, L.; Chen, Y.; Gao, L.; Wang, L.V. Single-shot real-time video recording of a photonic Mach cone induced by a scattered light pulse. Sci. Adv. 2017, 3, e1601814. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Inoue, T.; Matsunaka, A.; Funahashi, A.; Okuda, T.; Nishio, K.; Awatsuji, Y. Spatiotemporal observations of light propagation in multiple polarization states. Opt. Lett. 2019, 44, 2069–2072. [Google Scholar] [CrossRef] [PubMed]
  13. Kim, T.; Liang, J.; Zhu, L.; Wang, L.V. Picosecond-resolution phase-sensitive imaging of transparent objects in a single shot. Sci. Adv. 2020, 6, eaay6200. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Wang, P.; Liang, J.; Wang, L.V. Single-shot ultrafast imaging attaining 70 trillion frames per second. Nat. Commun. 2020, 11, 2091. [Google Scholar] [CrossRef]
  15. 31st International Congress of High-Speed Imaging and Photonics. Available online: https://www.ile.osaka-u.ac.jp/research/fps/ichsip-31/invited/index.html (accessed on 23 November 2020).
  16. Etoh, T. A high-speed video camera operating at 4500 fps. J. Inst. Telev. Eng. Jpn. 1992, 46, 543–545. (In Japanese) [Google Scholar]
  17. Kosonocky, W.; Yang, G.; Ye, C.; Kabra, R.; Xie, L.; Lawrence, J.; Mastrocolla, V.; Shallcross, F.; Patel, V. 360 × 360-element very-high-frame-rate burst image sensor. In Proceedings of the Digest of Technical Papers, 1996 IEEE International Solid-State Circuits Conference (ISSCC), San Fransisco, CA, USA, 10 February 1996; pp. 182–183. [Google Scholar]
  18. Etoh, T.G.; Poggemann, D.; Ruckelshausen, A.; Theuwissen, A.; Kreider, G.; Folkerts, H.-O.; Mutoh, H.; Kondo, Y.; Maruno, H.; Takubo, K.; et al. A CCD image sensor of 1 Mframes/s for continuous image capturing 103 frames. In Proceedings of the 2002 IEEE International Solid-State Circuits Conference, Digest of Technical Papers (Cat. No.02CH37315), San Francisco, CA, USA, 7 February 2002; pp. 46–443. [Google Scholar]
  19. Etoh, T.G.; Nguyen, D.H.; Dao, S.V.T.; Vo, C.L.; Tanaka, M.; Takehara, K.; Okinaka, T.; van Kuijk, H.; Klaassens, W.; Bosiers, J.; et al. A 16 Mfps 165kpixel backside-illuminated CCD. In Proceedings of the 2011 IEEE International Solid-State Circuits Conference, San Francisco, CA, USA, 20–24 February 2011; pp. 406–408. [Google Scholar]
  20. Crooks, J.; Marsh, B.; Turchetta, R.; Taylor, K.; Chan, W.; Lahav, A.; Fenigstein, A. Kirana: A solid-state megapixel uCMOS image sensor for ultrahigh speed imaging. In Proceedings of the Sensors, Cameras, and Systems for Industrial and Scientific Applications XIV, Burlingame, CA, USA, 19 February 2013; Volume 8659, pp. 36–49. [Google Scholar]
  21. Tochigi, Y.; Hanzawa, K.; Kato, Y.; Kuroda, R.; Mutoh, H.; Hirose, R.; Tominaga, H.; Takubo, K.; Kondo, Y.; Sugawa, S. A Global-Shutter CMOS Image Sensor with Readout Speed of 1-Tpixel/s Burst and 780-Mpixel/s Continuous. IEEE J. Solid-State Circuits 2013, 48, 329–338. [Google Scholar] [CrossRef]
  22. Mochizuki, F.; Kagawa, K.; Okihara, S.; Seo, M.; Zhang, B.; Takasawa, T.; Yasutomi, K.; Kawahito, S. Single-shot 200 Mfps 5×3-aperture compressive CMOS imager. In Proceedings of the 2015 IEEE International Solid-State Circuits Conference (ISSCC), San Francisco, CA, USA, 22–26 February 2015; Volume 58, pp. 116–117. [Google Scholar]
  23. Claus, L.; Fang, L.; Kay, R.; Kimmel, M.; Long, J.; Robertson, G.; Sanchez, M.; Stahoviak, J.; Trotter, D.; Porter, J.L. An overview of the Ultrafast X-ray Imager (UXI) program at Sandia Labs. In Proceedings of the Target Diagnostics Physics and Engineering for Inertial Confinement Fusion IV, San Diego, CA, USA, 31 August 2015; Volume 9591, p. 95910. [Google Scholar]
  24. Dao, V.T.S.; Ngo, N.; Nguyen, A.Q.; Morimoto, K.; Shimonomura, K.; Goetschalckx, P.; Haspeslagh, L.; De Moor, P.; Takehara, K.; Etoh, T.G. An Image Signal Accumulation Multi-Collection-Gate Image Sensor Operating at 25 Mfps with 32 × 32 Pixels and 1220 In-Pixel Frame Memory. Sensors 2018, 18, 3112. [Google Scholar] [CrossRef] [Green Version]
  25. Etoh, T.G.; Okinaka, T.; Takano, Y.; Takehara, K.; Nakano, H.; Shimonomura, K.; Ando, T.; Ngo, N.; Kamakura, Y.; Dao, V.T.S.; et al. Light-in-flight imaging by a silicon image sensor: Toward the theoretical highest frame rate. Sensors 2019, 19, 2247. [Google Scholar] [CrossRef] [Green Version]
  26. Kuroda, R.; Suzuki, M.; Sugawa, S. Over 100 million frames per second high speed global shutter CMOS image sensor. In Proceedings of the 32nd International Congress on High-Speed Imaging and Photonics, Enschede, The Netherlands, 28 January 2019; Volume 11051, pp. 57–62. [Google Scholar]
  27. Hart, P.A.; Carpenter, A.; Claus, L.; Damiani, D.; Dayton, M.; Decker, F.-J.; Gleason, A.; Heimann, P.; Hurd, E.; McBride, E.; et al. First X-ray test of the Icarus nanosecond-gated camera. In Proceedings of the X-ray Free-Electron Lasers: Advances in Source Development and Instrumentation V, Prague, Czech Republic, 24 April 2019; Volume 110380, p. 110380Q. [Google Scholar]
  28. Ogawa, S.; Asahara, R.; Minoura, Y.; Sako, H.; Kawasaki, N.; Yamada, I.; Miyamoto, T.; Hosoi, T.; Shimura, T.; Watanabe, H. Insights into thermal diffusion of germanium and oxygen atoms in HfO2/GeO2/Ge gate stacks and their suppressed reaction with atomically thin AlOx interlayers. J. Appl. Phys. 2015, 118, 235704. [Google Scholar] [CrossRef]
  29. Shimura, T.; Matsue, M.; Tominaga, K.; Kajimura, K.; Amamoto, T.; Hosoi, T.; Watanabe, H. Enhancement of photoluminescence from n-type tensile-strained GeSn wires on an insulator fabricated by lateral liquid-phase epitaxy. Appl. Phys. Lett. 2015, 107, 221109. [Google Scholar] [CrossRef]
  30. Oka, H.; Tomita, T.; Hosoi, T.; Shimura, T.; Watanabe, H. Lightly doped n-type tensile-strained single-crystalline {GeSn}-on-insulator structures formed by lateral liquid-phase crystallization. Appl. Phys. Express 2017, 11, 1011304. [Google Scholar] [CrossRef]
  31. Wada, Y.; Inoue, K.; Hosoi, T.; Shimura, T.; Watanabe, H. Demonstration of mm long nearly intrinsic {GeSn} single-crystalline wires on quartz substrate fabricated by nucleation-controlled liquid-phase crystallization. Jpn. J. Appl. Phys. 2019, 58, SBBK01. [Google Scholar] [CrossRef]
  32. Gariepy, G.; Krstajić, N.; Henderson, R.; Li, C.; Thomson, R.R.; Buller, G.S.; Heshmat, B.; Raskar, R.; Leach, J.; Faccio, D. Single-photon sensitive light-in-fight imaging. Nat. Commun. 2015, 6, 6021. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Vines, P.; Kuzmenko, K.; Kirdoda, J.; Dumas, D.C.S.; Mirza, M.M.; Millar, R.W.; Paul, D.J.; Buller, G.S. High performance planar germanium-on-silicon single-photon avalanche diode detectors. Nat. Commun. 2019, 10, 1086. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Martinez, N.J.D.; Derose, C.T.; Brock, R.W.; Starbuck, A.L.; Pomerene, A.T.; Lentine, A.L.; Trotter, D.C.; Davids, P.S. High performance waveguide-coupled Ge-on-Si linear mode avalanche photodiodes. Opt. Express 2016, 24, 19072–19081. [Google Scholar] [CrossRef]
  35. Kunikiyo, T.; Takenaka, M.; Kamakura, Y.; Yamaji, M.; Mizuno, H.; Morifuji, M.; Taniguchi, K.; Hamaguchi, C. A Monte Carlo simulation of anisotropic electron transport in silicon including full band structure and anisotropic impact-ionization model. J. Appl. Phys. 1994, 75, 297–312. [Google Scholar] [CrossRef]
  36. Jacoboni, C.; Canali, C.; Ottaviani, G.; Quaranta, A. A review of some charge transport properties of silicon. Solid-State Electron. 1977, 20, 77–89. [Google Scholar] [CrossRef]
  37. Jacoboni, C.; Nava, F.; Canali, C.; Ottaviani, G. Electron drift velocity and diffusivity in germanium. Phys. Rev. B 1981, 24, 1014–1026. [Google Scholar] [CrossRef] [Green Version]
  38. Omar, M.A.; Reggiani, L. Drift velocity and diffusivity of hot carriers in germanium: Model calculations. Solid-State Electron. 1987, 30, 1351–1354. [Google Scholar] [CrossRef]
  39. Jones, M.H.; Jones, S.H. Optical Properties of Silicon. Available online: https://www.univie.ac.at/photovoltaik/vorlesung/ss2014/unit4/optical_properties_silicon.pdf (accessed on 23 November 2020).
  40. Palik, E.D. Handbook of Optical Constants of Solids; Academic Press: New York, NY, USA, 1985; Volume 1, ISBN 978-0-08-054721-3. [Google Scholar]
  41. Green, M.A.; Keevers, M.J. Optical properties of intrinsic silicon at 300 K. Prog. Photovolt. Res. Appl. 1995, 3, 189–192. [Google Scholar] [CrossRef]
  42. Optical Absorption Coefficient Calculator. Available online: https://cleanroom.byu.edu/opticalcalc (accessed on 23 November 2020).
  43. Nemecek, A.; Zach, G.; Swoboda, R.; Oberhauser, K.; Zimmermann, H. Integrated BiCMOS pin photodetectors with high bandwidth and high responsivity. IEEE J. Sel. Top. Quantum Electron. 2006, 12, 1469–1475. [Google Scholar] [CrossRef]
  44. Sentaurus Device Monte Carlo User Guide, O-2018.06; Synopsys: Mountain View, CA, USA, 2018.
  45. Nguyen, A.Q.; Dao, V.T.S.; Shimonomura, K.; Takehara, K.; Etoh, T.G. Toward the Ultimate High-speed Image Sensor: From 10 ns to 50 ps. Sensors 2018, 18, 2407. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. Evolution of ultra-high-speed image sensors and imaging devices having achieved the super temporal resolution. The theoretical temporal resolution limit for a Ge PD is presented later in Section 3.3 of this paper. The slanted dashed line assumes that a chip with in-situ driver circuits is stacked to the sensor chip.
Figure 1. Evolution of ultra-high-speed image sensors and imaging devices having achieved the super temporal resolution. The theoretical temporal resolution limit for a Ge PD is presented later in Section 3.3 of this paper. The slanted dashed line assumes that a chip with in-situ driver circuits is stacked to the sensor chip.
Sensors 20 06895 g001
Figure 2. Definitions of the temporal resolution. (a) The theoretical temporal resolution limit, t T = 2 σ , (b) a practical temporal resolution, t P = t 95 t 05 3.3 σ , (c) for reference.
Figure 2. Definitions of the temporal resolution. (a) The theoretical temporal resolution limit, t T = 2 σ , (b) a practical temporal resolution, t P = t 95 t 05 3.3 σ , (c) for reference.
Sensors 20 06895 g002
Figure 3. The theoretical temporal resolution limit: Δ t = t T   = 2 σ .
Figure 3. The theoretical temporal resolution limit: Δ t = t T   = 2 σ .
Sensors 20 06895 g003
Figure 4. A model to analyze the standard deviation σ of the arrival time of electrons to the detection plane. W: thickness of the photodiode, s: the generation site of an electron distributing in the depth direction due to the exponential distribution of the penetration depth of light. The bulk is made of an intrinsic silicon. From Figure 5a, the electric fields between the backside and the detection plane are set at 2.5 V / μ m for simulations in Si and 0.42 V / μ m for simulations in Ge at which the drift velocities are 95% of the fully saturated values (see also Table 1).
Figure 4. A model to analyze the standard deviation σ of the arrival time of electrons to the detection plane. W: thickness of the photodiode, s: the generation site of an electron distributing in the depth direction due to the exponential distribution of the penetration depth of light. The bulk is made of an intrinsic silicon. From Figure 5a, the electric fields between the backside and the detection plane are set at 2.5 V / μ m for simulations in Si and 0.42 V / μ m for simulations in Ge at which the drift velocities are 95% of the fully saturated values (see also Table 1).
Sensors 20 06895 g004
Figure 5. The drift velocity of electrons and the longitudinal diffusion coefficient in Si and Ge: triangles in (a) are the experimental data after [36,37], and triangles and circles in (b) are the experimental data after [36,38]. Lines are calculated using Monte Carlo simulation codes.
Figure 5. The drift velocity of electrons and the longitudinal diffusion coefficient in Si and Ge: triangles in (a) are the experimental data after [36,37], and triangles and circles in (b) are the experimental data after [36,38]. Lines are calculated using Monte Carlo simulation codes.
Sensors 20 06895 g005
Figure 6. Absorption coefficients for Si, Ge, GaAs and InP [43]. The original date is assembled in [40].
Figure 6. Absorption coefficients for Si, Ge, GaAs and InP [43]. The original date is assembled in [40].
Sensors 20 06895 g006
Figure 7. The temporal resolution of the Si PD with incident light of 550 nm along the <111> direction. Solid line σ d d f : exact formulation of the DDF model; dashed line σ a p p : our approximate expression; σ d d f L : asymptotic value of σ d d f for an infinitesimal W ; σ d i f f and σ m i x : diffusion and mixing components of σ a p p where mixing is due to a combined effect of the drift velocity and the distribution of the penetration depth.
Figure 7. The temporal resolution of the Si PD with incident light of 550 nm along the <111> direction. Solid line σ d d f : exact formulation of the DDF model; dashed line σ a p p : our approximate expression; σ d d f L : asymptotic value of σ d d f for an infinitesimal W ; σ d i f f and σ m i x : diffusion and mixing components of σ a p p where mixing is due to a combined effect of the drift velocity and the distribution of the penetration depth.
Sensors 20 06895 g007
Figure 8. Standard deviation of various relevant parameters for a Si PD for 0.01 < W < 1   μ m . Wavelength of light λ = 550   nm ; σ M C (blue line) and σ M C S (red line): standard deviations of the arrival times calculated by our inhouse MC simulation code and by the Sentaurus MC simulation code; σ d i f f L and σ m i x L ; asymptotic expression of σ d i f f and σ m i x ; for an infinitesimal W; σ d i f f is a straight line proportional to W 1 , extending σ d i f f L (black line).
Figure 8. Standard deviation of various relevant parameters for a Si PD for 0.01 < W < 1   μ m . Wavelength of light λ = 550   nm ; σ M C (blue line) and σ M C S (red line): standard deviations of the arrival times calculated by our inhouse MC simulation code and by the Sentaurus MC simulation code; σ d i f f L and σ m i x L ; asymptotic expression of σ d i f f and σ m i x ; for an infinitesimal W; σ d i f f is a straight line proportional to W 1 , extending σ d i f f L (black line).
Sensors 20 06895 g008
Figure 9. Standard deviation of arrival time and temporal resolutions of the Ge PD for 550 nm illumination with the average penetration depth of 20 nm. σ M C S is calculated by a newly proposed method based on Sentaurus MC simulations. The other results in the plots are calculated with equations in Appendix A and Appendix B for the values in Table 1. The diffusion component σ d i f f is not drawn since it completely fits σ a p p for W < 0.1   μ m for Ge PD.
Figure 9. Standard deviation of arrival time and temporal resolutions of the Ge PD for 550 nm illumination with the average penetration depth of 20 nm. σ M C S is calculated by a newly proposed method based on Sentaurus MC simulations. The other results in the plots are calculated with equations in Appendix A and Appendix B for the values in Table 1. The diffusion component σ d i f f is not drawn since it completely fits σ a p p for W < 0.1   μ m for Ge PD.
Sensors 20 06895 g009
Figure 10. Standard deviations for Ge PDs: (a) wavelength of incident light λ = 1   μ m and the average penetration depth δ = 0.526   μ m , (b) λ = 1.5   μ m ,   δ = 2.04   μ m .
Figure 10. Standard deviations for Ge PDs: (a) wavelength of incident light λ = 1   μ m and the average penetration depth δ = 0.526   μ m , (b) λ = 1.5   μ m ,   δ = 2.04   μ m .
Sensors 20 06895 g010
Figure 11. A Si charge collection pyramid (a) vs. a Ge PD (b).
Figure 11. A Si charge collection pyramid (a) vs. a Ge PD (b).
Sensors 20 06895 g011
Table 1. Parameters used in the analyses of the standard deviation of the arrival time.
Table 1. Parameters used in the analyses of the standard deviation of the arrival time.
Si
E < 111 >
Ge
E < 100 >
Wavelength0.55 μ m0.55 μ m1 μ m
Penetration depth1.73 μ m20.0 nm526 nm
Critical E-field25 kV/cm4.2 kV/cm
Drift velocity9.19 × 106 cm/s5.8 × 106 cm/s
Diffusion coefficient10.8 cm2/s42.5 cm2/s
(1) Si and Ge are intrinsic. (2) The critical field is defined as the field at which the drift velocity is 95% of the saturated value. (3) After figures in [36,38], we measured that, for Si, the diffusion coefficients at the critical fields are 13.1, 12.1, 10.8, and 10.8 cm2/s for 77, 160, 200, and 300 K, respectively. They take close values for 200 and 300 K. For Ge, the diffusion coefficients at critical fields are 48.5, and 42.5 cm2/s at 77 and 190 K. Therefore, the diffusion coefficient for Ge at 300 K is assumed to be 42.5 cm2/s.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Ngo, N.H.; Nguyen, A.Q.; Bufler, F.M.; Kamakura, Y.; Mutoh, H.; Shimura, T.; Hosoi, T.; Watanabe, H.; Matagne, P.; Shimonomura, K.; et al. Toward the Super Temporal Resolution Image Sensor with a Germanium Photodiode for Visible Light. Sensors 2020, 20, 6895. https://doi.org/10.3390/s20236895

AMA Style

Ngo NH, Nguyen AQ, Bufler FM, Kamakura Y, Mutoh H, Shimura T, Hosoi T, Watanabe H, Matagne P, Shimonomura K, et al. Toward the Super Temporal Resolution Image Sensor with a Germanium Photodiode for Visible Light. Sensors. 2020; 20(23):6895. https://doi.org/10.3390/s20236895

Chicago/Turabian Style

Ngo, Nguyen Hoai, Anh Quang Nguyen, Fabian M. Bufler, Yoshinari Kamakura, Hideki Mutoh, Takayoshi Shimura, Takuji Hosoi, Heiji Watanabe, Philippe Matagne, Kazuhiro Shimonomura, and et al. 2020. "Toward the Super Temporal Resolution Image Sensor with a Germanium Photodiode for Visible Light" Sensors 20, no. 23: 6895. https://doi.org/10.3390/s20236895

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop