Next Article in Journal
Estimation of Axle Torque for an Agricultural Tractor Using an Artificial Neural Network
Next Article in Special Issue
3D Dosimetry Based on LiMgPO4 OSL Silicone Foils: Facilitating the Verification of Eye-Ball Cancer Proton Radiotherapy
Previous Article in Journal
EEG-Based Sleep Staging Analysis with Functional Connectivity
Previous Article in Special Issue
Performance of Optical Coupling Materials in Scintillation Detectors Post Temperature Exposure
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Fully Depleted Monolithic Active Microstrip Sensors: TCAD Simulation Study of an Innovative Design Concept

by
Lorenzo De Cilladi
1,2,*,
Thomas Corradino
3,4,
Gian-Franco Dalla Betta
3,4,
Coralie Neubüser
4 and
Lucio Pancheri
3,4 on behalf of the ARCADIA Collaboration
1
Dipartimento di Fisica, Università degli Studi di Torino, 10125 Torino, Italy
2
Istituto Nazionale di Fisica Nucleare (INFN), Sezione di Torino, 10125 Torino, Italy
3
Dipartimento di Ingegneria Industriale, Università degli Studi di Trento, 38123 Trento, Italy
4
Trento Institute for Fundamental Physics and Applications—Istituto Nazionale di Fisica Nucleare (TIFPA-INFN), 38123 Trento, Italy
*
Author to whom correspondence should be addressed.
Sensors 2021, 21(6), 1990; https://doi.org/10.3390/s21061990
Submission received: 18 January 2021 / Revised: 24 February 2021 / Accepted: 4 March 2021 / Published: 11 March 2021
(This article belongs to the Special Issue Ionizing Radiation Sensor and Detector)

Abstract

:
The paper presents the simulation studies of 10 μ m pitch microstrips on a fully depleted monolithic active CMOS technology and describes their potential to provide a new and cost-effective solution for particle tracking and timing applications. The Fully Depleted Monolithic Active Microstrip Sensors (FD-MAMS) described in this work, which are developed within the framework of the ARCADIA project, are compliant with commercial CMOS fabrication processes. A set of Technology Computer-Aided Design (TCAD) parametric simulations was performed in the perspective of an upcoming engineering production run with the aim of designing FD-MAMS, studying their electrical characteristics, and optimizing the sensor layout for enhanced performance in terms of low capacitance, fast charge collection, and low-power operation. A fine pitch of 10 μ m was chosen to provide high spatial resolution. This small pitch still allows readout electronics to be monolithically integrated in the inter-strip regions, enabling the segmentation of long strips and the implementation of distributed readout architectures. The effects of surface radiation damage expected for total ionizing doses of the order of 10 to 10 5 krad were also modeled in the simulations. The results of the simulations exhibit promising performance in terms of timing and low power consumption and motivate R&D efforts to further develop FD-MAMS; the results will be experimentally verified through measurements on the test structures that will be available from mid-2021.

1. Introduction

Charged particle tracking and timing are fundamental tools for both physics research and for numerous applications. Although a number of detection techniques are available, silicon detectors have become largely employed due to their versatility and to the parallel strong developments of the semiconductor industry. Various types of silicon sensors have been developed to meet the specific requirements of different experiments and applications, such as high spatial resolution, fast charge collection, low power consumption, high radiation tolerance, and low cost per unit area.
Silicon detectors are divided in two categories, namely hybrid and monolithic detectors. The former are made of two separate silicon elements, the sensor and the chip, which are interconnected through external bump or wire bonding. While the sensor hosts the sensing volume only, the chip integrates the front-end readout electronics. On the contrary, monolithic sensors, which are emerging as a valid alternative to hybrid detectors, embed the front-end electronics in the same silicon substrate which hosts the sensing volume, with benefits in terms of material budget, production yield, and fabrication cost, as they are produced with commercial microelectronics processes [1,2,3].
Due to their characteristics, monolithic sensors have recently raised a wide interest in different research fields; studies, proposals, and developments have been made for applications in high energy physics (HEP) [4,5,6], X-ray imaging [7,8], medical particle imaging [9], and space experiments [10].
The state of the art includes three main types of monolithic sensors. The first type, called Depleted Field Effect Transistors (DEPFETs), is capable of low noise operation, thanks to low sensor input capacitance [11]. DEPFET detectors have been developed and used for HEP applications [12], for X-ray imaging in space [13] and for free electron laser experiments [14]. The main limitation of DEPFETs is the need to reset their internal gate which can be quickly saturated by the leakage [15], thus making this technology not suitable for environments with high levels of non-ionizing radiation.
A second approach consists in the SOI (Silicon-On-Insulator) monolithic sensors. SOI sensors embed a buried-oxide layer separating a thin low-resistivity silicon layer, which hosts the integrated readout circuitry, from a thicker high-resistivity substrate, which serves as the sensitive detection region [16,17]. This technology allows a low capacitance to be obtained [17]; however, SOI sensors suffer from back-gate effect and have a reduced radiation hardness, due to accumulation of positive holes charges in the buried oxide layer after irradiation [18]. Strategies have been found to overcome these limitations and to recover from the Total Ionizing Dose (TID) [19], but, as a consequence, the fabrication process of SOI sensors have become highly specialized and not compliant with standard microelectronics production processes. This results in increased cost per unit area, which is a critical issue for large-area detector applications.
A third category of monolithic sensors is represented by CMOS (Complementary Metal-Oxide Semiconductor) sensors [20]. CMOS sensors were already in use for light detection when they were first proposed for charged particle tracking at the beginning of the 2000s [21]. Over the last few years, important advancements in CMOS sensors allowed them to be employed in many applications, eventually leading to large scale productions for particle trackers at collider experiments. The STAR pixel detector, which took data at the Relativistic Heavy Ion Collider (RHIC) from 2014 to 2016, was the first large area monolithic pixel tracker ever built, for a total of 0.16 m 2 [22]. These dimensions have been exceeded by the newly-constructed Inner Tracking System of the ALICE experiment at CERN, in which a total detector surface of about 10 m 2 is covered by ALPIDE CMOS monolithic active (A monolithic sensors is called “active” if it integrates a signal amplifier inside each pixel or strip.) pixel sensors (MAPS) [4].
These achievements demonstrate the level of maturity and reliability that CMOS sensors have recently reached. However, there is still room for further improvements, especially in terms of charge collection speed and radiation hardness, and possibility to push previous limits in terms of low power density, high spatial resolution, and signal to noise ratio (SNR).
Pixel detectors are the first choice for small scale applications and for vertex trackers at collider experiments [2] as they have an intrinsic capability of providing a two-dimensional position information [23]. On the other hand, microstrip sensors [24] are widely used as particle detectors for space applications and are a competitive option for particle trackers due to their high spatial resolution, simpler readout and much lower power density (i.e., power consumption per unit area) compared to pixel detectors. Particle collider experiments have employed silicon hybrid strip sensors in the past and are still developing and assembling new trackers based on this technology, as in the case of the Phase-2 Upgrades of the CMS Outer Tracker [25] and of the ATLAS Strip Inner Tracker [26]. Recent space experiments equipped with silicon hybrid microstrip trackers include FERMI-LAT [27], DAMPE [28], PAMELA [29], and AMS-02 [30]. Strip-like sensors integrated in a monolithic technology have been proposed by combining the outputs of 55 μ m × 55 μ m [31] or 40 μ m × 600 μ m [32] pixels in each column or row of a pixel matrix.
Spatial resolution of 1.25–1.3 μ m was achieved using hybrid silicon microstrip sensors with 25 μ m pitch [33]. However, it has recently been demonstrated with fully depleted double-SOI monolithic pixel sensors that the 1 μ m limit can be exceeded by semiconductor detectors [34]. The keys to a high spatial resolution with analogue readout are a fine microstrip pitch, a low sensor thickness to reduce Coulomb scattering and delta-ray emission, and an increased SNR, which can be achieved by reducing the leakage current and the sensor input capacitance to the readout electronics, but which is ultimately limited by the noise of the front-end electronics [24,33,35].
This paper presents the first investigation, design and simulation studies of CMOS Fully Depleted Monolithic Active Microstrip Sensors (FD-MAMS) with 10 μ m pitch for charged particle detection. Properly optimized sensor layouts may allow sub-micron resolution, improved radiation hardness, and fast timing performance, thanks to full depletion [6,36] in a power-saving and cost-effective commercial technology. Moreover, a further advantage of monolithic microstrips is the potential complexity reduction of the detector assembly compared to hybrid microstrip detectors. In fact, since many readout functions can be monolithically integrated on the same chip which hosts the sensing volume, 1-by-1 strip bonding to the external readout electronics would not be needed anymore. Hence, we studied and designed the FD-MAMS within the framework of the INFN ARCADIA (Advanced Readout CMOS Architectures with Depleted Integrated sensor Arrays) project in order to provide an innovative solution for satellite-based space trackers and for particle detectors at future collider experiments.
The results of the Technology Computer-Aided Design (TCAD) simulations (The TCAD simulations were produced using the Synopsys® Sentaurus (Version O-2018.06-SP2) software.) which allowed different MAMS design options to be compared in terms of sensor capacitance, reference voltage values, leakage current, and charge collection time and efficiency are presented; the effects of the inclusion of a silicon dioxide (SiO 2 ) layer on top of the sensor and of surface radiation damage on the sensor operating parameters are explored; the study of charge sharing between groups of adjacent strips when particles with different Linear Energy Transfer (LET) traverse the sensor is reported. A selection of the MAMS presented in this paper is going to be implemented in test structures which were submitted in November 2020 for an engineering production run.
The paper is organized as follows: Section 2 presents the sensor concept for the ARCADIA fully depleted CMOS monolithic microstrip sensors and illustrates the set of parametric simulations that were performed for the sensor design; Section 3 describes and discusses the results of the simulations; Section 4 presents the conclusions, the future perspectives, and the planned tests for the ARCADIA monolithic microstrip sensors.

2. The ARCADIA Sensor Concept

The ARCADIA project and its precursor, SEED (Sensor with Embedded Electronics Development), designed an innovative sensor concept [37,38] based on a modified 110 nm CMOS process developed in collaboration with LFoundry and compatible with their standard 110 nm CMOS process. Up to 6 metal layers can be stacked on top of the sensor, for a total metal and insulator thickness of about 4–5 μ m. The ARCADIA collaboration is developing a scalable event-driven readout architecture to cover detection surfaces of O(cm 2 ) while maintaining ultra-low power consumption. The target for pixel sensors is 10 mW/cm 2 to 20 mW/cm 2 at rates in the order of 100 MHz, but, for less dense particle environments (e.g., in space applications), a dedicated low-power operation mode implements a cyclic pulling of the data packets from each section of the pixel matrix and disables most of the serializers and data transceivers, further reducing the total power consumption of the chip.
In our project, an n-on-n sensor concept enabling full substrate depletion over tens or hundreds of microns and allowing full CMOS electronics to be implemented was employed. A simplified view of the sensor cross section is visible in Figure 1. The process allows to achieve sensor thicknesses from 50 to 400 μ m. A high resistivity n-type substrate was used and constitutes the active volume. The sensing n-well node, located on top of the sensor, collects the electrons produced by ionization due to particles traversing the active detection volume.
N-doped and p-doped wells intended to host pMOSFETs and nMOSFETs, respectively, are shielded by a deep p-well, which allows the integration of full CMOS electronics and, hence, more complex digital functions, when necessary. In fact, the deep p-well prevents the n-wells hosting pMOSFETs from competing with the n-doped sensing node in the collection of the charge, thus avoiding loss of charge collection efficiency.
A p+ boron-doped region sits at the backside of the n-substrate, thus forming a pn-junction; when a negative bias voltage V b a c k is applied to the backside p+ contact, sensor depletion starts from the pn-junction at the bottom of the sensor and eventually extends to the whole sensor, if the backside voltage is sufficiently large. Since the high voltage needed for sensor depletion is applied at the backside, it is possible to maintain the voltage V n w e l l applied to the front n-well electrode below 1 V and to use low-voltage integrated electronics (1.2 V transistors), which is more radiation-resistant and has lower noise. Full sensor depletion allows fast charge collection by drift (beneficial to enhance the timing performance), higher charge collection efficiency, deeper collection depth, and larger SNR; it also leads to improved radiation tolerance, as charge losses by trapping are reduced [5]. Since thicker sensors need higher backside bias voltage to reach full depletion, termination structures composed of multiple floating guard rings are used to avoid early breakdown at the edges of the pn-junction.
An additional n-type epitaxial layer, with lower resistivity than the substrate, is integrated between the n-type substrate and the deep p-wells. Its aim is to better control the potential barrier below the deep p-well, in order to delay the onset of the punch through current described in details in Section 2.2.2.
The feasibility of this sensor concept and approach to Fully Depleted monolithic CMOS sensors was proven in the framework of the SEED project [37,38], and the design activities described in this article have as a starting point the experimental and simulation work performed in that project. The main constraint on the design was to reduce to the minimum the modifications to the foundry’s standard fabrication process, in order to guarantee easy portability to different commercial fabrication processes. In the top side of the sensor, the only custom implant is the deep p-well. An additional constraint was the request to keep a low voltage on the sensor top surface to deal with possible embedded 1.2 V transistors.
The upcoming ARCADIA engineering run will include different design options of FD-CMOS monolithic sensors, both pixelated and strip-like. Large-area ( 1.3 × 1.3 cm 2 ) pixel demonstrators with embedded CMOS electronics and pixel test structures ( 0.5 × 0.5 and 1.5 × 1.5 mm 2 ) without integrated readout circuitry [39] are foreseen, with pitches ranging from 10 to 50 μ m. The test structures will include, as well, the innovative MAMS and will allow a detailed characterization of these sensors. The 3D TCAD simulations performed to design the first FD-MAMS will be presented and discussed in the following.

2.1. TCAD Simulations

Three-dimensional TCAD simulations were employed as a tool to design the sensors and study their performance. The use of 3D simulations is necessary to have a more realistic domain and results which are more accurate and less affected by boundary conditions. Furthermore, we were also interested in studying the charge collection dynamics after a particle crosses the sensor, and this is more straightforward with 3D simulations. A fine pitch of 10 μ m was chosen for the microstrips in order to explore the characteristics and performance of a sensor layout which pushes the requirements on both spatial and timing resolution. Different sensor thicknesses foreseen for the production runs were simulated. Variations in the sensor layout and operating parameters were tested to study and optimize the sensor response. The simulated sensor options take into account the limitations imposed by the foundry’s fabrication process, especially for the n-well and p-well sizes. The strip simulations investigated sensor layouts which pushed the design to the limits of the process requirements.
All the TCAD simulations were performed at a temperature of 300 K. A standard simulation domain including three 50 μ m long, 50 μ m thick, 10 μ m pitch microstrips is shown as an example in Figure 2a. The n-doped substrate is shown in green, the epitaxial layer in yellow, the microstrip sensing n-wells in red, the p-wells in blue, and the less doped deep p-well in light blue. The default value for V n w e l l is 0.8 V. The p-wells, instead, are kept at a voltage V p w e l l = 0 V.
One of the simulated sensor options has been specifically designed to allow for CMOS digital library cells to be integrated along the strips and is shown in Figure 2b. This sensor variant would allow the deployment of complex CMOS digital functions along the strip for distributed signal processing. We observed that the n-wells dedicated to the implementation of PMOS transistors and shielded by the deep p-well do not significantly influence the electrical characteristics of the detector in the TCAD simulation results. Therefore, we did not include them in the simulations.
The deep p-well can be removed in the test structures that will be used to characterize the sensor (see Figure 2c), and the necessary CMOS front-end electronics can be deployed at the end of the strips in the chip periphery. Sensors without the deep p-well were simulated, as well.
Different n-well, p-well, and deep p-well sizes were considered to find the optimal layout in terms of sensor performance. Simulations were also employed to predict the effects that possible production uncertainties can have on the sensor operating parameters and electrical characteristics. For instance, the thickness and resistivity of the epitaxial layer may vary within a confidence range around their typical specified value (see Appendix A). Three-dimensional simulations for the different cases were run and compared. Some simulation parameters were fine-tuned using characterization results from a previous set of test structures, produced in the framework of the SEED project [37].

2.2. Electrical and Transient Simulation

In this section, the simulations performed to extract the sensor electrical characteristics and to study the charge collection dynamics are briefly illustrated. Shared definitions and conventions on simulation setups and operating parameters were agreed for all the ARCADIA sensor simulations and are also described by Neubüser et al. [39]. The strip length in the upcoming production run will be 1.2 cm. However, MAMS with lengths of 50 μ m were simulated in order to run the TCAD simulations in a reasonable computational time. The results were then scaled to the desired length.

2.2.1. Depletion Voltage

Sensor depletion starts at the backside, where the pn-junction between the n-type substrate and the p+ contact is located. If no negative bias voltage is applied to the backside contact, the sensor is not fully depleted, and the collection n-wells are not isolated. This means that a resistive path exists between the n-type sensing nodes (see Figure 3a). Therefore, if a voltage difference is applied between two adjacent n-wells, a current will flow between them.
As the negative voltage applied to the backside contact increases, the space charge region enlarges through the high resistivity substrate, eventually merging with the depletion volume which surrounds the pn-junctions formed between the n-type substrate or epitaxial layer and the deep p-wells. At this point, the sensor is fully depleted, the resistive path between the sensing nodes is closed, and the collection n-wells are isolated; this is shown in Figure 3b. In this condition, no current (except for the leakage current) will flow among adjacent n-wells even when different voltages are applied to them.
This behavior can be observed in the orange example IV (current-voltage) curve in Figure 4 ( I n w e l l , u n b a l a n c e d ). The simulated domain shown in Figure 2a was used. In this simulation, a voltage unbalance of 10 mV was applied between adjacent strips: the first n-well was biased at 0.79 V, the central one at 0.8 V, and the third one at 0.81 V. The curve shows the current measured at the sensing node of the central strip as a function of | V b a c k | . A current of about 1 nA is measured at V b a c k = 0 V. As the backside voltages increases and the space charge region enlarges, the current starts decreasing, and eventually reaches a plateau at a current of about 10 5 nA. This baseline corresponds to the leakage current (green IV curve, I n w e l l , l e a k a g e ). The backside voltage at which the single microstrips become isolated and the plateau is reached is the sensor depletion voltage V d p l ; this voltage is evaluated as the intersection point between the exponential decay fitting of the IV curve decreasing segment and the baseline.
Figure 5 shows the simulated electrostatic potential and electric field maps at V b a c k = V d p l in a cross section of a 3-strip domain with all the n-wells at V n w e l l = 0.8 V. Electric field lines are plotted on top of both the electrostatic potential and the electric field maps.

2.2.2. Punch-Through

If V b a c k exceeds a certain value, a hole current flowing between the shallow p-doped backside region and the (deep) p-well exponentially increases. This condition is known as punch-through, and the hole current is the punch-through current [40]. We define the voltage corresponding to the onset of the punch-through as V p t . The onset of the punch-through currents can be observed from the blue IV curve in Figure 4 ( I p w e l l ), which shows the absolute value of the current measured at the top p-well contacts as a function of | V b a c k | . The dip in the curve, corresponding to the point of sign inversion of the current, was defined as V p t . The simulation domain includes three 50 μ m long, 50 μ m thick, 10 μ m pitch microstrips. In this case, the n-wells are all biased at V n w e l l = 0.8 V, which is the default value.
Sensor operation in low punch-through regime can be tolerated, whereas a punch-through current density exceeding a value of the order of 100 μ A/cm 2 ought to be avoided, as it determines a substantial increase in the power consumption of the whole detector. For this reason, we chose V b a c k = V p t as a safe reference sensor operating voltage; this is the operating point for all the results shown in the following, if not stated differently. The sensor power density can be defined as p d = V b a c k × ( I p w e l l + I n w e l l ) A , where I n w e l l and I p w e l l are the currents flowing at the sensing node and at the top p-well contacts, respectively, and A is the top surface area of the simulated microstrip domain. In order to quantify the maximum acceptable backside bias voltage that limits the absorbed power density, the value V p d at which p d = 0.1 mW/cm 2 was extracted from the simulated IV curves (see Figure 4).
Figure 6 shows the hole current density at two different | V b a c k | > | V p t | in the simulation domain used to extract the I p w e l l curve of Figure 4. In Figure 6a, a backside voltage exceeding V p t by 1 V was chosen, while, in Figure 6b, V b a c k was set to V p d . An increase in the hole current density of several orders of magnitude can be observed below the deep p-wells and in the substrate.
Care had to be taken to ensure that | V d p l | < | V p t | in the designed sensors. In this way, full depletion is reached before the onset of the punch-through. Moreover, the voltage operating range between V d p l and V p t , defined as Δ V o p = | V p t V d p l | , should be large enough to ensure safe operation in full depletion before the onset of the punch-through even if deviations from the simulated design occur in the sensor fabrication process.

2.2.3. Leakage Current

The same sensor domain and n-well voltage configuration used for the extraction of V p t was also used to evaluate the sensor leakage current I l e a k . The leakage current is defined as the current flowing at the collection nodes in full depletion and in absence of external stimuli, such as particles or radiation. The leakage current as a function of the backside bias voltage is shown in Figure 4 as a green curve ( I l e a k ). In the example shown in Figure 4, a value of 10 fA was extracted for I l e a k at V b a c k = V p t .

2.2.4. Sensor Capacitance

The sensor CV (capacitance-voltage) curve was simulated through AC simulations with a frequency of 10 kHz using the same sensor domain employed for V p t and I l e a k evaluation, with V n w e l l = 0.8 V. The major contribution to the sensor capacitance C s e n s , which is the input capacitance seen by the DC-coupled front-end electronics, originates from the lateral capacitance between the collection n-well and the surrounding p-wells. Thus, it is important to minimize this contribution by a careful selection of the distance between the edge of the collection n-well and the p-wells; we call this distance “gap” (see Figure 1). An example CV curve is shown in red in Figure 4, with the capacitance per unit length considered. In the example of Figure 4, a value of about 0.33 fF/ μ m was obtained at V b a c k = V p t .
It has to be mentioned that in these sensors the depletion voltage does not necessarily correspond to the voltage of minimum capacitance. The reason for this is the presence of the epitaxial layer, which is located far from the backside pn-junction and has a lower resistivity than the substrate. Therefore, the depletion of the epitaxial layer begins after the depletion of the substrate and progresses more slowly with voltage. Full depletion of the whole sensor, including the epitaxial layer, and minimum capacitance are only reached at | V b a c k | > | V d p l | . From this point, both capacitance and leakage current values will be intended at V b a c k = V p t .
A central focus of the design was the minimization of the sensor capacitance. In fact, low input capacitance to the DC-coupled CMOS readout electronics allows for low-noise readout, low analog power [5], and, in particular, SNR maximization. Large input capacitance worsens the noise levels and the speed of the front-end electronics [6].

2.2.5. Surface Radiation Damage

In the set of simulations performed to study the properties of MAMS, a silicon dioxide (SiO 2 ) layer was added on the top-side of the sensor. In addition to this, surface damage was modeled to evaluate the effects of Total Ionizing Dose (TID) on the sensor electrical properties.
The impact of surface radiation damage was modeled following the AIDA-2020-D7.4 report [41]. The model introduces fixed positive oxide charges and band-gap acceptor/donor defect levels (trap states) at the Si-SiO 2 interface. The concentrations of oxide charges and defect levels start from a fixed value before irradiation (i.e., with the only inclusion of the SiO 2 surface layer, at d o s e = 0 ) and increase with the dose provided to the sensors. The dependence of the oxide charge density Q o x (charges × cm 2 ), of the acceptor integrated interface trap state density N i n t a c c (cm 2 ), and of the donor integrated interface trap state density N i n t d o n (cm 2 ) on the dose is shown in Figure 7. Pre-irradiation values, shown as dotted horizontal lines in Figure 7, are Q o x = 6.5 × 10 10 charges × cm 2 , N i n t a c c = 2.0 × 10 9 cm 2 , and N i n t d o n = 2.0 × 10 9 cm 2 .
In the simulations, the effects of the inclusion of the SiO 2 layer and of the radiation damage on the leakage current, sensor capacitance, depletion voltage, and punch-through voltage were investigated and will be discussed in Section 3.

2.2.6. Transient Simulations

TCAD transient simulations were run to study the sensor charge collection process in response to particles traversing the simulated microstrip domain. These simulations also let us identify the most relevant layout parameters to be optimized for improving the sensor performance in terms of fast and uniform charge collection irrespective of the particle incidence position. The transient simulations employ the Synopsys® Sentaurus TCAD HeavyIon model, described in Sentaurus Device User Guide [42]. The HeavyIon model gives an analytical description of the amount of charge generated within a 3D cylindrical distribution along the incident particle track. Two main parameters have to be passed to the HeavyIon model: the Linear Energy Transfer (LET), defined as the average deposited charge per unit length, and the transverse size of the charge deposition volume generated around the particle trajectory. We chose the charge transverse distribution profile to be gaussian around the particle track.
Two extreme cases in terms of particle impact position were studied to evaluate the uniformity of charge collection time and charge collection efficiency. Particle trajectories perpendicular to the sensor surface were considered. In the best-case scenario, the particle impact point corresponds to the centre of a microstrip, which is the centre of a collection n-well. On the contrary, in the worst-case scenario, the particle traverses the sensor at the edge between two adjacent microstrips, i.e., in the middle of a p-well. In Figure 8, the two cases and the corresponding numbering of the strips are illustrated. This conventional strip nomenclature will be used in the following when referring to transient simulations.
In order to save computational time, a reduced TCAD simulation domain that employs the symmetries was used. This reduced domain corresponds to a quarter of the full domain, with the particle incident in the corner of the domain instead of in the centre. An example for the best-case scenario is shown in Figure 9. The collected charge and current signals were then scaled to reproduce the full domain case, which includes nine or ten 100 μ m long microstrips in the best-case and worst-case scenario, respectively (Figure 8). These numbers and size of strips guarantee that the amount of deposited charge reaching the borders of the simulation domain is negligible. The correctness of this strategy was verified and confirmed by comparing the results of a simulation with a quarter domain and of a simulation with full domain.
An example of current signals I n w e l l ( t ) measured at the microstrip sensing nodes when a particle crosses the microstrip domain is shown in Figure 10a. We defined as charge collection efficiency for the i-th strip (CCE i ) the integral of the current signal I n w e l l , i ( t ) extracted from the i-th strip and normalized at the total charge Q t o t deposited in the sensor by the particle, according to the formula:
C C E i ( t ) = 0 t I n w e l l , i ( t ) d t Q t o t = 0 t I n w e l l , i ( t ) d t L E T × d S i ,
where d S i is the sensor thickness. The total charge collection efficiency CCE for the whole simulated domain is defined as:
C C E ( t ) = i = 1 N s t r i p s C C E i ( t ) ,
where N s t r i p s is the total number of strips in the simulated domain. The total CCE at the end of the charge collection process (i.e., at t = t e n d = 30 ns, which was observed to be enough for complete charge collection) has to be equal to 100% in the absence of recombination:
C C E ( t = t m a x ) = 100 % .
The CCE i as a function of time is shown in Figure 10b, for strip number 1. The times needed for collecting the 95% and 99% of the total deposited charge were evaluated and referred to as t 95 and t 99 , respectively. These values were compared for different design options and used to select the fastest sensor layouts.
The spatial mesh of the transient simulations was forced to be finer around the particle trajectory to more accurately simulate the charge deposition and the drift of electrons and holes from their generation points along the particle track towards the electrodes. Additionally, the time step of the transient simulations was fine tuned to guarantee the necessary accuracy while keeping the computational time requirement economical. We observed that these adjustments prevented the simulations from giving unphysical results.

2.3. Determination of the LET for Heavy Nuclei

Since MAMS are an interesting candidate for tracking detectors in space applications, the charge collection was studied not only for minimum ionizing particles (MIPs) but also for heavy nuclei of interest for in-orbit astroparticle experiments. The LET values of carbon and oxygen ions were studied in Geant4 (version 10.6 patch 01) simulations. The simulation setup consists of a 50 μ m thick silicon layer immersed in air, with a transverse size of 1 × 1 cm 2 . The particle gun was positioned 15 cm in front of the centre of the silicon layer. The G4EmPenelopePhysics physics list was used to model the electromagnetic processes and the necessary precision on the energy deposited within the silicon was achieved with a maximum step size of 1 μ m [43]. This setup has been validated with simulations of 300 MeV muons, where the known LET for MIPs of about 80 electron-hole (e-h) pairs per μ m [44] (or 1.28 × 10 5 pC per μ m) in silicon could be reproduced. The LETs for carbon (C 12 + ) and oxygen (O 16 + ) ions at their minimum ionization were computed from their most probable energy loss (i.e., the most probable value of the straggling or Landau functions [45,46]). Figure 11 shows the LET as a function of the particle energy obtained for C and O ions traversing 50 μ m of silicon. The energies E m i n at which C and O ions are at the minimum of ionization were found to be 35 GeV and 60 GeV, respectively. The corresponding LETs are 45.6 × 10 5 pC/ μ m and 83.0 × 10 5 pC/ μ m, which results in 36 and 65 times the MIP value. This is consistent with the expected scaling from the Bethe-Bloch formula.
We were especially interested in studying the charge sharing among the microstrips surrounding the particle impact point and the charge collection time at different LETs. This will be reported and discussed in Section 3.

3. Results and Discussions

In this section, the results of the TCAD simulations will be presented in connection with the design objectives. As mentioned in Section 1, the main targets of the FD-MAMS design were the following.
  • To enhance the spatial resolution. The microstrip pitch of 10 μ m means an intrinsic spatial resolution of p i t c h 12 2.9 μ m with a digital readout, which can be improved, thanks to charge sharing and with an analog readout.
  • To minimize the sensor capacitance C s e n s at V b a c k = V p t .
  • To obtain fast and uniform charge collection, irrespective of the particle incidence position. This will enhance the sensor timing capabilities and will reduce the dead-time between successive particle detections.
For reasons of space available for MAMS in the first ARCADIA engineering run, only a few sensor layout options could be included. Hence, a set of parametric simulations was needed to identify the best performing sensor layouts. The deep p-well, when present, was kept the same size as the p-well. The expression “p-well and deep p-well” will be contracted and referred to as “(deep) p-well”. In the legends of the figures, the abbreviation “dpw” will be used for deep p-well.

3.1. SiO 2 Layer and Surface Damage

A first group of TCAD simulation studies was aimed at investigating the effects of the SiO 2 layer and of surface TID damage on the FD-MAMS characteristics. The model we employed was presented in Section 2.2.5. As can be seen from Figure 12, for one of the selected 50 μ m thick microstrip layouts, the inclusion of the SiO 2 layer with a minimum concentration of traps and oxide charges ( d o s e = 0 ) determines a small increase of about 5% in the leakage current I l e a k from 20.8 fA to 22.0 fA. The sensor capacitance C s e n s is strongly affected by the inclusion of the SiO 2 layer, as it increases by 31% from 0.26 fF/ μ m to 0.34 fF/ μ m. Both I l e a k and C s e n s are found to rise with increasing dose. The minimum dose that we considered is 50 krad, as the model is not validated for lower doses [41]. Figure 13, instead, shows the effect of the SiO 2 layer and of the TID on V d p l and on V p t . The effect of the dose on these two values is smaller than in the case of I l e a k and C s e n s . Furthermore, V d p l and V p t are influenced by the dose in opposite directions, which results in a slight increase in the operating range Δ V o p with increasing dose.
The reason for the capacitance increase even after the simple inclusion of the SiO 2 layer was found to be the introduction of positive oxide charges at the Si-SiO 2 interface [39]. In fact, the model that we adopted foresees a significant positive oxide charge concentration Q o x = 6.5 × 10 10 charges/cm 2 already at dose = 0. The positive oxide charges attract free electrons from the n-type silicon epitaxial layer towards the Si-SiO 2 interface in the gap. This determines an increase in the electron concentration in the gap, as illustrated in Figure 14 and Figure 15, which behaves as an extension of the collection n-well.

3.2. Capacitance Minimization

The sizes of both the collection n-well and of the gap were found to contribute to C s e n s and were adjusted to reach minimum capacitance. The SiO 2 layer influences C s e n s in different ways for different gap sizes. Thus, the sensor capacitance was evaluated with and without the SiO 2 layer to quantify its effect. Figure 16 shows the trend of C s e n s as a function of the gap size for 50 μ m thick microstrips. The different sensor thicknesses considered (50, 100, and 300 μ m) were found not to influence the sensor capacitance. Both the case with fixed minimum-size n-well and variable (deep) p-well (blue curves) and the case with fixed minimum-size (deep) p-well and variable n-well (orange curve) were studied. The dash-dotted lines refer to simulations without the surface SiO 2 layer, whereas solid lines to the case with SiO 2 layer included with minimal oxide charge and trap concentration.
The reason for which smaller gaps with fixed n-wells could not be investigated is referred to as channel choking, a condition that inhibits sensor operation; this condition is explained in Section 3.3. The vertical grey band in Figure 16 and in the following ones corresponds to the forbidden region due to the constraints on n-well and (deep) p-well minimum sizes imposed by the fabrication process. The leftmost limit of the grey band is still permitted.
Variations of n-well and of (deep) p-well size do not lead to the same C s e n s for the same gap size. A fixed n-well size with SiO 2 layer included shows a trend that is not monotonic, but has a minimum at slightly less than 0.34 fF/ μ m. This effect is caused by the electron accumulation in the gap at the Si-SiO 2 interface. However, the difference in C s e n s between the minimum-capacitance option and the sensor layout at the edge of the forbidden region is lower than 2%. There was, as expected, no benefit from having large n-wells. The sensor capacitance increases with the n-well size, as can be seen from the blue curve in Figure 16. Therefore, the best layout for minimum C s e n s has the smallest possible n-well size and sufficiently small (deep) p-well.
After fixing the n-well to the minimum possible size, we compared the sensor capacitance for layouts with deep p-well (Figure 17, orange curve) and without deep p-well (Figure 17, green curve). It was observed that removing the deep p-well helps in reducing the sensor capacitance.

3.3. Reference and Operating Voltages

The influence of the n-well size on the operating voltages was negligible compared to the effect of the (deep) p-well size. Figure 18 presents the effect of the (deep) p-well size effect on V d p l and on V p t for 50 μ m thick sensors with minimum-sized n-well. Both the cases with (orange curves) and without deep p-well (green curves) were considered. The voltage values are reported for the case of dose = 0.
In all the layouts considered in Figure 18, the onset of the punch through happens at voltages sufficiently larger than the depletion voltage. Outside of the forbidden region (grey band), the operating range Δ V o p is always between 4.2 V and 6.2 V, or between the 23% and the 41% of V d p l . This is sufficient for safe sensor operation, even in the hypothesis of possible doping inhomogeneities among adjacent microstrips or slight deviations from the doping design values. Similar observations on Δ V o p have been made for 100 μ m thick and 300 μ m thick sensors.
It can be observed in Figure 18 that smaller (deep) p-wells result in increased V d p l and V p t . This can be interpreted as follows. Large p-doped surfaces below the (deep) p-wells create wider pn-junctions with the n-doped epitaxial layer, thus facilitating the depletion of the underlying epitaxial layer at lower voltages. On the other hand, large (deep) p-wells also lower the potential barrier that prevents the direct flow of holes towards the substrate. This results in the earlier onset of the punch through hole current between the (deep) p-wells and the backside p+ region. Sensors without the deep p-well showed higher reference voltages. In fact, the presence of a deep p-well reduces the epitaxial layer thickness below the p-wells, thus requiring a lower voltage to achieve both full depletion and the onset of punch-through currents.
For sensors with too large deep p-well, a phenomenon that we defined as channel choking was observed. It consists in the closure of the conductive channel below the collection n-well due to the lateral merging of the closely adjacent depletion regions formed at the junctions between the deep p-wells and the n-epitaxial layer. In this condition, even though the space charge region of the backside junction has not reached the surface yet, the n-wells are already isolated from one another at Vback = 0 V. Hence, the process of charge collection is inhibited by the strong potential barrier present below the n-wells. No channel choking was observed for sensor layouts without the deep p-well.
With a minimum-sized collection n-well and within the foundry’s process requirements, a maximum gap of about 2.1 μ m is allowed to have enough space for monolithically integrated CMOS electronics in the inter-strip region. This specific layout option is not affected by the channel choking and is, hence, operational. If, however, NMOS-only electronics is integrated along the strips, the smallest p-well size considered in our study can be used.
For completeness, Figure 19a illustrates the dependence of V p t and V d p l on the sensor thickness for the sensor layout with minimum sizes for the n-well and for the (deep) p-well. The trend is linear over a wide range of thicknesses, both with and without the deep p-well. In addition, the operating voltage Δ V o p = V p t V d p l linearly increases with the sensor thickness, as shown in Figure 19b. The sensor thickness investigated was extended down to 20 μ m, well below the smallest thickness (i.e., 50 μ m) of the sensors that will be produced in the first ARCADIA engineering run, as the study of thin sensors was functional for improving the sensor timing performance in the simulations.
The voltage V p d at which the power density is 0.1 mW/cm 2 was found to be about 4–5 V above V p t for 50 μ m thick microstrips, 7–8 V for 100 μ m thick microstrips, and 18–20 V for 300 μ m thick microstrips when the deep p-well was included.
The V n w e l l voltage was varied with the aim of finding possible improvements in the sensor performances. The results for V d p l and V p t are shown in Figure 20a, where the vertical red line indicates the default value of 0.8 V. A minimum V n w e l l of about 0.5 V is necessary to satisfy the condition | V p t | > | V d p l | . Moreover, an increase in V n w e l l has several interesting effects. First of all, it allows the sensor full depletion to be reached at lower (in absolute value) backside voltages. Secondly, it also shifts the onset of the punch through towards larger | V b a c k | , thus increasing the operating range Δ V o p . Finally, as shown in Figure 20b, higher V n w e l l implies lower sensor capacitance. At V n w e l l = 3 V, the capacitance can be reduced to about 0.2 fF/ μ m.

3.4. Charge Collection Studies

As described in Section 2.2.6, TCAD transient simulations were used to study the charge collection dynamics. In order to select the layouts with the optimal performance in terms of fast and uniform charge collection, the effect of the (deep) p-well size on the charge collection time at V b a c k = V p t was evaluated. The time t 95 needed to collect 95% of the total charge deposited in the simulated sensor domain is plotted in Figure 21 for 50 μ m thick sensors and LET = 1.28 × 10 5 pC/ μ m (1 MIP) as a function of the gap size and with fixed minimum n-well size.
Microstrips with large gaps, and hence small (deep) p-wells, are to be preferred for fast charge collection. The reason for this is a higher | V p t | , which enables sensor operation at a larger | V b a c k | . The consequent stronger electric field in the sensor results in higher charge velocity in the silicon substrate. For the same reason, microstrip sensors without deep p-well revealed a significantly faster charge collection in both the best-case and the worst-case scenario. The options with small (deep) p-wells also show uniform charge collection for different particle incidence positions. The difference in t 95 for the best-case and worst-case scenarios is below 0.1 ns for the fastest permitted options. This result is also achieved, thanks to the fine microstrip pitch of 10 μ m. The channel choking, as described in Section 3.3, limits the deep p-well size as the potential barrier below the sensing node slows down the electron collection. This problem, as shown in Figure 21, can be avoided by removing the deep p-well.
A gap size of about 2.1 μ m, which guarantees enough space for embedded CMOS electronics in the inter-strip region according to the foundry’s layout rules, is not in the channel choking region. Hence, it preserves the sensor functionality without a dramatic slow down of the charge collection process.
Figure 22 demonstrates that the proposed MAMS guarantee fast sensor response also under heavily ionizing particles. The charge collection time is only weakly proportional to the charge deposited by the incident particle within an LET range of 1.28 × 10 5 pC/ μ m to 128 × 10 5 pC/ μ m. A 50 μ m thick sensor was considered in Figure 22, and the LET values corresponding to 1 MIP, carbon (C) ion, and oxygen (O) ion at their minimum of ionization are highlighted as vertical green lines. t 99 is added to show that the time needed for complete charge collection is only slightly larger than t 95 , due to a small fraction of charge collected by the strips adjacent to the central one. t 95 and t 99 were never found to exceed 2 ns and 3 ns, respectively, in 50 μ m thick sensors.
As we discussed in Section 3.4, the first strategy for improving the timing performance of the proposed microstrip sensors is to remove the deep p-well in order to obtain higher | V p t | . However, we also investigated other ways to increase | V p t | and to speed up the charge collection. In particular, as shown in Section 3.3, a larger V n w e l l shifts the onset of the punch-through current towards higher | V b a c k | . Hence, we explored the effects of V n w e l l on the charge collection time.
In a strip readout system, timing information can be retrieved only from the strips collecting most of the charge (i.e., strip(s) number 1, following the nomenclature of Figure 8), as they provide a signal with sufficiently large SNR. Therefore, in order to study the sensor timing performance and after verifying through t 95 that the total deposited charge is quickly collected in the whole simulation domain, we considered the time t 95 c e n t r a l needed to collect 95% of the charge in the central strip(s).
Figure 23 shows the dependence of t 95 c e n t r a l at V b a c k = V p t and with LET = 1.28 × 10 5 pC/ μ m on the voltage applied to the sensing node. A 50 μ m thick sensor with a layout optimized for fast charge collection was considered. A significant improvement could be reached at larger V n w e l l . For the option without deep p-well and at V n w e l l = 3 V, t 95 c e n t r a l is 0.84 ns in the best-case and 0.94 ns in the worst-case scenario. If we assume an electron drift saturation velocity of 1 × 10 7 cm/s in silicon at a temperature of 300 K [47], the minimum drift time for electrons that have to cover a 50 μ m distance is 0.5 ns. This explains the saturation observed in Figure 23 and demonstrates the fast charge collection and the promising timing capabilities of the proposed MAMS.
A way to further reduce the collection time is to explore thinner sensors. Figure 24 demonstrates that the charge collection time t 95 c e n t r a l is proportional to the sensor thickness. For these simulations, V n w e l l was set to 0.8 V and a 1 MIP LET was considered. Even at thicknesses as large as 300 μ m, t 95 c e n t r a l does not exceed 6 ns. In the best-case scenario without the deep p-well, reducing the sensor thickness from 50 μ m to 40 μ m, 30 μ m, and 20 μ m results in a decrease in t 95 c e n t r a l of 15%, 33%, and 50%, respectively. Analogous proportionality was observed for t 95 . Therefore, for future production runs, thinner sensors could be considered for the enhancement of the timing performance.
A set of TCAD simulations was dedicated to study the charge sharing among adjacent microstrips when particles with different LETs traverse the sensor. Charge sharing is relevant for improving the spatial resolution, especially with analog readout, and is enhanced by fine microstrip pitches and large sensor thicknesses. On the contrary, it is reduced at higher V b a c k for a fixed sensor thickness.
In Figure 25, the case of a 300 μ m thick sensor at V b a c k = V p t is presented for the best-case scenario. The total charge collected by each strip (identified using the nomenclature of Figure 8) is plotted versus the LET. The black horizontal line indicates a possible charge threshold corresponding to 10% of a MIP at the single strip level. A comparison with the sensors that will be produced in the first ARCADIA engineering run will allow deeper investigation on the charge sharing, a fine tuning of the simulations and studies aimed at evaluating the spatial resolution of 10 μ m pitch MAMS.

4. Conclusions

In this work, we presented detailed TCAD simulations of CMOS-based FD-MAMS, which may find use for tracking and timing in particle and nuclear physics, space and medical applications. The results of the TCAD simulations, performed to design the 10 μ m pitch FD-MAMS, demonstrate their fast and uniform charge collection, which encourages their practicality for various applications, even under heavily ionizing particles. The effect of surface ionizing radiation damage was investigated, and the layout parameters were varied to achieve a minimum capacitance, beneficial for electronic noise reduction. The possibility to operate the sensor in full depletion and at low-power density (i.e., before the onset of the punch through current) was verified in the simulations. A preference for small collection diodes and small (deep) p-wells emerged for obtaining lower capacitance and faster sensor response. Additionally, even with a pitch of only 10 μ m, it was found that there is enough space for hosting transistors in the inter strip region according to the foundry layout rules without slowing down the charge collection process. The first FD-MAMS samples will be produced in the upcoming ARCADIA engineering production run in mid 2021 and will allow the simulation results to be compared with experimental data from electrical characterization, laser and beam irradiation tests. The promising results of the FD-MAMS simulations will translate into further R&D activities to enhance the sensor performance in terms of low capacitance and high timing and spatial resolution.

Author Contributions

Data curation, L.D.C.; Formal analysis, L.D.C.; Investigation, L.D.C.; Supervision, C.N. and L.P. Writing – original draft, L.D.C.; Writing—review & editing, T.C., G.-F.D.B., C.N. and L.P. All authors have read and agreed to the published version of the manuscript.

Funding

The research activity presented in this article has been carried out in the framework of the ARCADIA experiment funded by the Istituto Nazionale di Fisica Nucleare (INFN), CSN5. The activity has also been supported by the project “Dipartimento di Eccellenza”, Physics Department of the University of Torino (Dipartimento di Fisica—Università degli Studi di Torino), Italy, funded by MUR (Ministero dell’Università e della Ricerca).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Expected Effects from Epitaxial Layer Thicknesses

Possible variations in the epitaxial layer thickness of (−15%; +30%) communicated by the foundry with respect to the reference value induced us to investigate their effect on the operating parameters. While the sensor capacitance was observed not to be influenced, both V d p l and V p t showed a linear dependence on the epitaxial layer thickness. This behavior is presented in Figure A1.
Figure A1. V d p l and V p t as a function of the epitaxial layer thickness, expressed as percentage variation with respect to the reference thickness.
Figure A1. V d p l and V p t as a function of the epitaxial layer thickness, expressed as percentage variation with respect to the reference thickness.
Sensors 21 01990 g0a1

References

  1. Garcia-Sciveres, M.; Wermes, N. A review of advances in pixel detectors for experiments with high rate and radiation. Rep. Prog. Phys. 2018, 81, 66101. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Wermes, N. Pixel detectors for charged particles. Nucl. Instrum. Methods Phys. Res. Sect. A 2009, 604, 370–379. [Google Scholar] [CrossRef] [Green Version]
  3. Wermes, N. Pixel detectors... where do we stand? Nucl. Instrum. Methods Phys. Res. Sect. A 2019, 924, 44–50. [Google Scholar] [CrossRef] [Green Version]
  4. Mager, M.; The ALICE Collaboration. ALPIDE, the Monolithic Active Pixel Sensor for the ALICE ITS upgrade. Nucl. Instrum. Methods Phys. Res. Sect. A 2016, 824, 434–438. [Google Scholar] [CrossRef] [Green Version]
  5. Pernegger, H.; Bates, R.; Buttar, C.; Dalla, M.; van Hoorne, J.W.; Kugathasan, T.; Maneuski, D.; Musa, L.; Riedler, P.; Riegel, C.; et al. First tests of a novel radiation hard CMOS sensor process for Depleted Monolithic Active Pixel Sensors. J. Instrum. 2017, 12, P06008. [Google Scholar] [CrossRef]
  6. Wang, T.; Rymaszewski, P.; Barbero, M.; Degerli, Y.; Godiot, S.; Guilloux, F.; Hemperek, T.; Hirono, T.; Krüger, H.; Liu, J.; et al. Development of a depleted monolithic CMOS sensor in a 150 nm CMOS technology for the ATLAS inner tracker upgrade. J. Instrum. 2017, 12, C01039. [Google Scholar] [CrossRef]
  7. Wunderer, C.B.; Marras, A.; Bayer, M.; Glaser, L.; Göttlicher, P.; Lange, S.; Pithan, F.; Scholz, F.; Seltmann, J.; Shevyakov, I.; et al. The PERCIVAL soft X-ray imager. J. Instrum. 2014, 9, C03056. [Google Scholar] [CrossRef]
  8. Hatsui, T.; Omodani, M.; Kudo, T.; Kobayashi, K.; Imamura, T.; Ohmoto, T.; Iwata, A.; Ono, S.; Kirihara, Y.; Kameshima, T.; et al. A direct-detection X-ray CMOS image sensor with 500 μm thick high resistivity silicon. In Proceedings of the 2013 International Image Sensor Workshop, Snowbird, UT, USA, 12–16 June 2013; Volume 3, p. 4. Available online: http://www.imagesensors.org/Past%20Workshops/2013%20Workshop/2013%20Papers/03-5_058_hatsui_paper.pdf (accessed on 15 December 2020).
  9. Mattiazzo, S.; Baruffaldi, F.; Bisello, D.; Di Ruzza, B.; Giubilato, P.; Iuppa, R.; La Tessa, C.; Pantano, D.; Pozzobon, N.; Ricci, E.; et al. iMPACT: An innovative tracker and calorimeter for proton computed tomography. IEEE Tran. Rad. Plasma Med. Sci. 2017, 2, 345–352. [Google Scholar] [CrossRef]
  10. Scotti, V.; Osteria, G.; The CSES-LIMADOU Collaboration. The High Energy Particle Detector onboard CSES-02 satellite. In Proceedings of the 36th International Cosmic Ray Conference (ICRC), Madison, WI, USA, 24 July–1 August 2019; Available online: https://pos.sissa.it/358/135/pdf (accessed on 4 December 2020).
  11. Rummel, S.; Andricek, L.; Moser, H.-G.; Richter, R. Intrinsic properties of DEPFET active pixel sensors. J. Instrum. 2009, 4, P03003. [Google Scholar] [CrossRef] [Green Version]
  12. Mariñas, C.; Vos, M. The Belle-II DEPFET pixel detector: A step forward in vertexing in the superKEKB flavour factory. Nucl. Instrum. Methods Phys. Res. Sect. A 2011, 650, 59–63. [Google Scholar] [CrossRef]
  13. Treis, J.; Andricek, L.; Aschauer, F.; Heinzinger, K.; Herrmann, S.; Hilchenbach, M.; Lauf, T.; Lechner, P.; Lutz, G.; Majewski, P.; et al. MIXS on BepiColombo and its DEPFET based focal plane instrumentation. Nucl. Instrum. Methods Phys. Res. Sect. A 2010, 624, 540–547. [Google Scholar] [CrossRef]
  14. Lutz, G.; Lechner, P.; Porro, M.; Strüder, L.; De Vita, G. DEPFET sensor with intrinsic signal compression developed for use at the XFEL free electron laser radiation source. Nucl. Instrum. Methods Phys. Res. Sect. A 2010, 624, 528–532. [Google Scholar] [CrossRef]
  15. Lutz, G.; Porro, M.; Aschauer, S.; Wölfel, S.; Strüder, L. The DEPFET sensor-amplifier structure: A method to beat 1/f noise and reach sub-electron noise in pixel detectors. Sensors 2016, 16, 608. [Google Scholar] [CrossRef] [Green Version]
  16. Kucewicz, W.; Bulgheroni, A.; Caccia, M.; Grabiec, P.; Marczewski, J.; Niemiec, H. Development of monolithic active pixel detector in SOI technology. Nucl. Instrum. Methods Phys. Res. Sect. A 2005, 541, 172–177. [Google Scholar] [CrossRef]
  17. Lan, H.; Wang, Y.; Chen, J.; Cao, F.; Hu, H.F. SoI Monolithic Active Pixel Sensors for Radiation Detection Applications: A Review. IEEE Sens. J. 2015, 15, 2732–2746. [Google Scholar] [CrossRef]
  18. Hagino, K.; Yarita, K.; Negishi, K.; Oono, K.; Hayashida, M.; Kitajima, M.; Kohmura, T.; Tsuru, T.G.; Tanaka, T.; Uchida, H.; et al. Radiation damage effects on double-SOI pixel sensors for X-ray astronomy. Nucl. Instrum. Methods Phys. Res. Sect. A 2020, 978, 164435. [Google Scholar] [CrossRef]
  19. Asano, M.; Sekigawa, D.; Hara, K.; Aoyagi, W.; Honda, S.; Tobita, N.; Arai, Y.; Miyoshi, T.; Kurachi, I.; Tsuboyama, T.; et al. Characteristics of non-irradiated and irradiated double SOI integration type pixel sensor. Nucl. Instrum. Methods Phys. Res. Sect. A 2016, 831, 315–321. [Google Scholar] [CrossRef] [Green Version]
  20. Turchetta, R.; French, M.; Manolopoulos, S.; Tyndel, M.; Allport, P.; Bates, R.; O’Shea, V.; Hall, G.; Raymond, M. Monolithic active pixel sensors (MAPS) in a VLSI CMOS technology. Nucl. Instrum. Methods Phys. Res. Sect. A 2003, 501, 251–259. [Google Scholar] [CrossRef]
  21. Turchetta, R.; Berst, J.; Casadei, B.; Claus, G.; Colledani, C.; Dulinski, W.; Hu, Y.; Husson, D.; Le Normand, J.; Riester, J.; et al. A monolithic active pixel sensor for charged particle tracking and imaging using standard VLSI CMOS technology. Nucl. Instrum. Methods Phys. Res. Sect. A 2001, 458, 677–689. [Google Scholar] [CrossRef]
  22. Contin, G.; Greiner, L.; Schambach, J.; Szelezniak, M.; Anderssen, E.; Bell, J.; Cepeda, M.; Johnson, T.; Qiu, H.; Ritter, H.-G.; et al. The STAR MAPS-based PiXeL detector. Nucl. Instrum. Methods Phys. Res. Sect. A 2018, 907, 60–80. [Google Scholar] [CrossRef] [Green Version]
  23. Rossi, L.; Fischer, P.; Rohe, T.; Wermes, N. Pixel Detectors: From Fundamentals to Applications; Springer Science & Business Media: Berlin/Heidelberg, Germany, 2006; ISBN 9783540283331. [Google Scholar]
  24. Peisert, A. Silicon microstrip detectors. In: Instrumentation in High Energy Physics. World Sci. Adv. Ser. Direct. High Energy Phys. 1992, 9, 1–79. [Google Scholar] [CrossRef]
  25. Chowdhury, S.R.; The CMS collaboration. The Phase-2 Upgrade of the CMS Outer Tracker. Nucl. Instrum. Methods Phys. Res. Sect. A 2020, 979, 164432. [Google Scholar] [CrossRef]
  26. David, C. A new strips tracker for the upgraded ATLAS ITk detector. J. Instrum. 2018, 13, C01003. [Google Scholar] [CrossRef] [Green Version]
  27. Atwood, W.; Bagagli, R.; Baldini, L.; Bellazzini, R.; Barbiellini, G.; Belli, F.; Borden, T.; Brez, A.; Brigida, M.; Caliandro, G.; et al. Design and initial tests of the Tracker-converter of the Gamma-ray Large Area Space Telescope. Astropart. Phys. 2007, 28, 422–434. [Google Scholar] [CrossRef] [Green Version]
  28. Azzarello, P.; Ambrosi, G.; Asfandiyarov, R.; Bernardini, P.; Bertucci, B.; Bolognini, A.; Cadoux, F.; Caprai, M.; De Mitri, I.; Domenjoz, M.; et al. The DAMPE silicon–tungsten tracker. Nucl. Instrum. Methods Phys. Res. Sect. A 2016, 831, 378–384. [Google Scholar] [CrossRef] [Green Version]
  29. Adriani, O.; Bonechi, L.; Bongi, M.; Castellini, G.; D’Alessandro, R.; Gabbanini, A.; Grandi, M.; Papini, P.; Ricciarini, S.B.; Spillantini, P.; et al. The magnetic spectrometer of the PAMELA satellite experiment. Nucl. Instrum. Methods Phys. Res. Sect. A 2003, 511, 72–75. [Google Scholar] [CrossRef]
  30. Lübelsmeyer, K.; von Dratzig, A.S.; Wlochal, M.; Ambrosi, G.; Azzarello, P.; Battiston, R.; Becker, R.; Becker, U.; Bertucci, B.; Bollweg, K.; et al. Upgrade of the Alpha Magnetic Spectrometer (AMS-02) for long term operation on the International Space Station (ISS). Nucl. Instrum. Methods Phys. Res. Sect. A 2011, 654, 639–648. [Google Scholar] [CrossRef] [Green Version]
  31. Benhammadi, S.; Dopke, J.; Guerrini, N.; Phillips, P.; Sedgwick, I.; Villani, G.; Wilson, F.; Zhang, Z.; Warren, M.; Allport, P.P.; et al. DECAL: A Reconfigurable Monolithic Active Pixel Sensor for use in Calorimetry and Tracking. In Topical Workshop on Electronics for Particle Physics; 2019; Volume 2, p. 6. Available online: https://s3.cern.ch/inspire-prod-files-4/4f56e107d64ecc5b165344c1695f60c1 (accessed on 7 December 2020).
  32. Han, Y.; Zhu, H.; Affolder, A.; Arndt, K.; Bates, R.; Benoit, M.; Di Bello, F.; Blue, A.; Bortoletto, D.; Buckland, M.; et al. Study of CMOS strip sensor for future silicon tracker. Nucl. Instrum. Methods Phys. Res. Sect. A 2020, 981, 164520. [Google Scholar] [CrossRef]
  33. Straver, J.; Toker, O.; Weilhammer, P.; Colledani, C.; Dulinski, W.; Turchetta, R.; Bosisio, L. One micron spatial resolution with silicon strip detectors. Nucl. Instrum. Methods Phys. Res. Sect. A 1994, 348, 485–490. [Google Scholar] [CrossRef] [Green Version]
  34. Sekigawa, D.; Endo, S.; Aoyagi, W.; Hara, K.; Honda, S.; Tsuboyama, T.; Yamada, M.; Ono, S.; Togawa, M.; Ikegmi, Y.; et al. Fine-pixel detector FPIX realizing sub-micron spatial resolution developed based on FD-SOI technology. In Proceedings of the International Conference on Technology and Instrumentation in Particle Physics, Beijing, China, 22–26 May 2017; Springer: Singapore, 2017; pp. 331–338. Available online: https://link.springer.com/chapter/10.1007/978-981-13-1316-5_62 (accessed on 21 November 2020).
  35. Turchetta, R. Spatial resolution of silicon microstrip detectors. Nucl. Instrum. Methods Phys. Res. Sect. A 1993, 335, 44–58. [Google Scholar] [CrossRef]
  36. Snoeys, W.; Rinella, G.A.; Hillemanns, H.; Kugathasan, T.; Mager, M.; Musa, L.; Riedler, P.; Reidt, F.; Van Hoorne, J.; Fenigstein, A.; et al. A process modification for CMOS monolithic active pixel sensors for enhanced depletion, timing performance and radiation tolerance. Nucl. Instrum. Methods Phys. Res. Sect. A 2017, 871, 90–96. [Google Scholar] [CrossRef]
  37. Pancheri, L.; Olave, J.; Panati, S.; Rivetti, A.; Cossio, F.; Rolo, M.; DeMaria, N.; Giubilato, P.; Pantano, D.; Mattiazzo, S. A 110 nm CMOS process for fully-depleted pixel sensors. J. Instrum. 2019, 14, C06016. [Google Scholar] [CrossRef]
  38. Pancheri, L.; Giampaolo, R.A.; Di Salvo, A.; Mattiazzo, S.; Corradino, T.; Giubilato, P.; Santoro, R.; Caccia, M.; Margutti, G.; Olave, J.E.; et al. Fully Depleted MAPS in 110-nm CMOS Process with 100–300 μm Active Substrate. IEEE Trans. Electron. Devices 2020, 67, 2393–2399. [Google Scholar] [CrossRef]
  39. Neubüser, C.; Corradino, T.; Dalla Betta, G.-F.; De Cilladi, L.; Pancheri, L. Sensor design optimization of innovative low-power, large area MAPS for HEP and applied science. arXiv 2020, arXiv:2011.09723. [Google Scholar]
  40. Chu, J.L.; Persky, G.; Sze, S.M. Thermionic injection and space-charge-limited current in reach-through p+np+ structures. J. Appl. Phys. 1972, 43, 3510–3515. [Google Scholar] [CrossRef]
  41. Passeri, D.; Morozzi, A. TCAD Radiation Damage Model. AIDA-2020-D7.4 Deliverable Report. 2019. Available online: http://cds.cern.ch/search?p=AIDA-2020-D7.4 (accessed on 21 November 2020).
  42. Sentaurus Device User Guide, Version O-2018.06; Synopsys Inc.: Mountain View, CA, USA, 2018; pp. 681–683.
  43. Centis Vignali, M. Available online: https://github.com/mcentis/muonOnSilicon (accessed on 10 July 2020).
  44. Olive, K.A.; Particle Data Group. Particle Detectors at Accelerators. Chin. Phys. C 2014, 38, 090001. [Google Scholar] [CrossRef]
  45. Tanabashi, M.; Particle Data Group. Passage of Particles Through Matter. Phys. Rev. D 2018, 98, 030001, (2019 update). [Google Scholar] [CrossRef] [Green Version]
  46. Bichsel, H. A method to improve tracking and particle identification in TPCs and silicon detectors. Nucl. Instrum. Methods Phys. Res. Sect. A 2006, 562, 154–197. [Google Scholar] [CrossRef]
  47. Canali, C.; The CMS Collaboration. Electron and hole drift velocity measurements in silicon and their empirical relation to electric field and temperature. IEEE Trans. Electron. Devices 1975, 22, 1045–1047. [Google Scholar] [CrossRef]
Figure 1. ARCADIA monolithic sensor concept. The dotted arrows indicate the drift path of electrons (e) and holes (h) generated by a particle crossing the sensor. The voltages V n w e l l and V b a c k applied to the sensor contacts are shown in green. Not to scale.
Figure 1. ARCADIA monolithic sensor concept. The dotted arrows indicate the drift path of electrons (e) and holes (h) generated by a particle crossing the sensor. The voltages V n w e l l and V b a c k applied to the sensor contacts are shown in green. Not to scale.
Sensors 21 01990 g001
Figure 2. Example Technology Computer-Aided Design (TCAD) 3D sensor domains for ARCADIA microstrips (top row) and corresponding cross sections (bottom row). (a) Standard simulation domain for sensors with the deep p-well. (b) Addition of n-wells above the deep p-wells. (c) Simulated ARCADIA microstrips without deep p-wells.
Figure 2. Example Technology Computer-Aided Design (TCAD) 3D sensor domains for ARCADIA microstrips (top row) and corresponding cross sections (bottom row). (a) Standard simulation domain for sensors with the deep p-well. (b) Addition of n-wells above the deep p-wells. (c) Simulated ARCADIA microstrips without deep p-wells.
Sensors 21 01990 g002
Figure 3. Depletion process in ARCADIA microstrips. (a): cross section of a sensor before full depletion is reached. (b): cross section of a fully depleted Monolithic Active Microstrip Sensors (MAMS). The orange lines indicate the edge of the depletion region.
Figure 3. Depletion process in ARCADIA microstrips. (a): cross section of a sensor before full depletion is reached. (b): cross section of a fully depleted Monolithic Active Microstrip Sensors (MAMS). The orange lines indicate the edge of the depletion region.
Sensors 21 01990 g003
Figure 4. Example characteristic IV and CV curves extracted from TCAD simulations of ARCADIA monolithic sensors. The red vertical axis refer to the sensor capacitance ( C s e n s ) CV curve.
Figure 4. Example characteristic IV and CV curves extracted from TCAD simulations of ARCADIA monolithic sensors. The red vertical axis refer to the sensor capacitance ( C s e n s ) CV curve.
Sensors 21 01990 g004
Figure 5. Electrostatic potential map (a) and electric field map (b) for a group of three ARCADIA microstrip sensors at V b a c k = V d p l . The electric field lines are plotted on top of both maps.
Figure 5. Electrostatic potential map (a) and electric field map (b) for a group of three ARCADIA microstrip sensors at V b a c k = V d p l . The electric field lines are plotted on top of both maps.
Sensors 21 01990 g005
Figure 6. Hole current density in a simulated sensor domain including three microstrips in punch-through condition at two different V b a c k .
Figure 6. Hole current density in a simulated sensor domain including three microstrips in punch-through condition at two different V b a c k .
Sensors 21 01990 g006
Figure 7. Dependence of the oxide charge density Q o x , acceptor integrated interface trap state density N i n t a c c , and donor integrated interface trap state density N i n t d o n on the dose for the surface radiation damage model described by Passeri et al. [41]. Pre-irradiation values are shown as horizontal dotted lines.
Figure 7. Dependence of the oxide charge density Q o x , acceptor integrated interface trap state density N i n t a c c , and donor integrated interface trap state density N i n t d o n on the dose for the surface radiation damage model described by Passeri et al. [41]. Pre-irradiation values are shown as horizontal dotted lines.
Sensors 21 01990 g007
Figure 8. Best-case and worst-case scenarios considered in the TCAD transient simulations. The microstrips are represented as adjacent grey blocks, and the particle traversing the domain is shown as an orange cylinder. The nomenclature used to identify the microstrips (from 1 to 5) is illustrated.
Figure 8. Best-case and worst-case scenarios considered in the TCAD transient simulations. The microstrips are represented as adjacent grey blocks, and the particle traversing the domain is shown as an orange cylinder. The nomenclature used to identify the microstrips (from 1 to 5) is illustrated.
Sensors 21 01990 g008
Figure 9. Example reduced TCAD domain used in transient simulations (best-case scenario). The microstrips are labeled following the nomenclature illustrated in Figure 8. A crossing particle is represented as an orange cylinder hitting the corner of the simulated reduced domain.
Figure 9. Example reduced TCAD domain used in transient simulations (best-case scenario). The microstrips are labeled following the nomenclature illustrated in Figure 8. A crossing particle is represented as an orange cylinder hitting the corner of the simulated reduced domain.
Sensors 21 01990 g009
Figure 10. Simulated current signals (a) and corresponding charge collection efficiency CCE i (b) in the best-case and worst-case scenarios for strip number 1 in an example 50 μ m thick microstrip domain. A particle track with an Linear Energy Transfer (LET) of 1.28 × 10 5 pC/ μ m was simulated.
Figure 10. Simulated current signals (a) and corresponding charge collection efficiency CCE i (b) in the best-case and worst-case scenarios for strip number 1 in an example 50 μ m thick microstrip domain. A particle track with an Linear Energy Transfer (LET) of 1.28 × 10 5 pC/ μ m was simulated.
Sensors 21 01990 g010
Figure 11. Dependence of the LET on the energy of carbon ions (C 12 + , blue) and oxygen ions (O 16 + , orange) incident on 50 μ m thick silicon. The LET values were evaluated through Geant4 simulations. The red vertical lines indicate the minimum ionization energies for the two particle species.
Figure 11. Dependence of the LET on the energy of carbon ions (C 12 + , blue) and oxygen ions (O 16 + , orange) incident on 50 μ m thick silicon. The LET values were evaluated through Geant4 simulations. The red vertical lines indicate the minimum ionization energies for the two particle species.
Sensors 21 01990 g011
Figure 12. Leakage current I l e a k (green) and sensor capacitance C s e n s (red) as a function of the total ionizing dose for a 50 μ m thick microstrip sensor. The values obtained in simulations with and without the SiO 2 layer in the absence of irradiation are shown as horizontal lines and referred to as “dose = 0” and “no SiO 2 layer”, respectively.
Figure 12. Leakage current I l e a k (green) and sensor capacitance C s e n s (red) as a function of the total ionizing dose for a 50 μ m thick microstrip sensor. The values obtained in simulations with and without the SiO 2 layer in the absence of irradiation are shown as horizontal lines and referred to as “dose = 0” and “no SiO 2 layer”, respectively.
Sensors 21 01990 g012
Figure 13. Depletion voltage V d p l (blue) and punch-through voltage V p t (orange) as a function of the total ionizing dose for a 50 μ m thick microstrip sensor. The values obtained in simulations with and without the silicon dioxide layer in the absence of irradiation are shown as horizontal lines.
Figure 13. Depletion voltage V d p l (blue) and punch-through voltage V p t (orange) as a function of the total ionizing dose for a 50 μ m thick microstrip sensor. The values obtained in simulations with and without the silicon dioxide layer in the absence of irradiation are shown as horizontal lines.
Sensors 21 01990 g013
Figure 14. Schematic illustration of the electron accumulation in the gap between the collection n-well and the surrounding p-wells due to the positive oxide charges introduced at the Si-silicon dioxide (SiO 2 ) interface.
Figure 14. Schematic illustration of the electron accumulation in the gap between the collection n-well and the surrounding p-wells due to the positive oxide charges introduced at the Si-silicon dioxide (SiO 2 ) interface.
Sensors 21 01990 g014
Figure 15. Electron density in an example microstrip simulation domain without (a) and with (b) the SiO 2 layer on top of the sensors.
Figure 15. Electron density in an example microstrip simulation domain without (a) and with (b) the SiO 2 layer on top of the sensors.
Sensors 21 01990 g015
Figure 16. C s e n s as a function of the gap size for different sensor layout configuration. The vertical grey band is the forbidden region due to fabrication constraints; its leftmost limit is still permitted.
Figure 16. C s e n s as a function of the gap size for different sensor layout configuration. The vertical grey band is the forbidden region due to fabrication constraints; its leftmost limit is still permitted.
Sensors 21 01990 g016
Figure 17. Sensor capacitance C s e n s as a function of the gap size for different sensor layout configurations with and without the deep p-well.
Figure 17. Sensor capacitance C s e n s as a function of the gap size for different sensor layout configurations with and without the deep p-well.
Sensors 21 01990 g017
Figure 18. Sensor depletion voltage V d p l and punch-through voltage V p t as a function of the gap size for different sensor layout configurations. The orange region indicates the forbidden region due to the observed channel choking.
Figure 18. Sensor depletion voltage V d p l and punch-through voltage V p t as a function of the gap size for different sensor layout configurations. The orange region indicates the forbidden region due to the observed channel choking.
Sensors 21 01990 g018
Figure 19. Dependence of V d p l and V p t (a) and of the operating voltage range Δ V o p (b) on the sensor thickness.
Figure 19. Dependence of V d p l and V p t (a) and of the operating voltage range Δ V o p (b) on the sensor thickness.
Sensors 21 01990 g019
Figure 20. V d p l and V p t (a) and sensor capacitance (b) as a function of V n w e l l . The vertical red line indicates the default value of V n w e l l = 0.8 V.
Figure 20. V d p l and V p t (a) and sensor capacitance (b) as a function of V n w e l l . The vertical red line indicates the default value of V n w e l l = 0.8 V.
Sensors 21 01990 g020
Figure 21. t 95 as a function of the gap size for best-case and worst-case scenarios.
Figure 21. t 95 as a function of the gap size for best-case and worst-case scenarios.
Sensors 21 01990 g021
Figure 22. t 95 (blue) and t 99 (red) as a function of the LET for best-case and worst-case scenarios.
Figure 22. t 95 (blue) and t 99 (red) as a function of the LET for best-case and worst-case scenarios.
Sensors 21 01990 g022
Figure 23. t 95 c e n t r a l as a function of the voltage V n w e l l applied to the sensing node for best-case and worst-case scenarios. The vertical red line indicates the default value of V n w e l l = 0.8 V.
Figure 23. t 95 c e n t r a l as a function of the voltage V n w e l l applied to the sensing node for best-case and worst-case scenarios. The vertical red line indicates the default value of V n w e l l = 0.8 V.
Sensors 21 01990 g023
Figure 24. t 95 c e n t r a l as a function of the sensor thickness for best-case and worst-case scenarios.
Figure 24. t 95 c e n t r a l as a function of the sensor thickness for best-case and worst-case scenarios.
Sensors 21 01990 g024
Figure 25. Charge sharing among adjacent microstrips. The total charge collected by strips 1 to 5 (following the nomenclature illustrated in Figure 8) is shown as a function of the LET.
Figure 25. Charge sharing among adjacent microstrips. The total charge collected by strips 1 to 5 (following the nomenclature illustrated in Figure 8) is shown as a function of the LET.
Sensors 21 01990 g025
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

De Cilladi, L.; Corradino, T.; Dalla Betta, G.-F.; Neubüser, C.; Pancheri, L., on behalf of the ARCADIA Collaboration. Fully Depleted Monolithic Active Microstrip Sensors: TCAD Simulation Study of an Innovative Design Concept. Sensors 2021, 21, 1990. https://doi.org/10.3390/s21061990

AMA Style

De Cilladi L, Corradino T, Dalla Betta G-F, Neubüser C, Pancheri L on behalf of the ARCADIA Collaboration. Fully Depleted Monolithic Active Microstrip Sensors: TCAD Simulation Study of an Innovative Design Concept. Sensors. 2021; 21(6):1990. https://doi.org/10.3390/s21061990

Chicago/Turabian Style

De Cilladi, Lorenzo, Thomas Corradino, Gian-Franco Dalla Betta, Coralie Neubüser, and Lucio Pancheri on behalf of the ARCADIA Collaboration. 2021. "Fully Depleted Monolithic Active Microstrip Sensors: TCAD Simulation Study of an Innovative Design Concept" Sensors 21, no. 6: 1990. https://doi.org/10.3390/s21061990

APA Style

De Cilladi, L., Corradino, T., Dalla Betta, G. -F., Neubüser, C., & Pancheri, L., on behalf of the ARCADIA Collaboration. (2021). Fully Depleted Monolithic Active Microstrip Sensors: TCAD Simulation Study of an Innovative Design Concept. Sensors, 21(6), 1990. https://doi.org/10.3390/s21061990

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