Simulating Monolithic Active Pixel Sensors: A Technology-Independent Approach Using Generic Doping Profiles

The optimisation of the sensitive region of CMOS sensors with complex non-uniform electric fields requires precise simulations, and this can be achieved by a combination of electrostatic field simulations and Monte Carlo methods. This paper presents the guiding principles of such simulations, using a CMOS pixel sensor with a small collection electrode and a high-resistivity epitaxial layer as an example. The full simulation workflow is described, along with possible pitfalls and how to avoid them. For commercial CMOS processes, detailed doping profiles are confidential, but the presented method provides an optimisation tool that is sufficiently accurate to investigate sensor behaviour and trade-offs of different sensor designs without knowledge of proprietary information. The workflow starts with detailed electric field finite element method simulations in TCAD, using generic doping profiles. Examples of the effect of varying different parameters of the simulated sensor are shown, as well as the creation of weighting fields, and transient pulse simulations. The fields resulting from TCAD simulations can be imported into the Allpix Squared Monte Carlo simulation framework, which enables high-statistics simulations, including modelling of stochastic fluctuations from the underlying physics processes of particle interaction. Example Monte Carlo simulation setups are presented and the different parts of a simulation chain are described. Simulation studies from small collection electrode CMOS sensors are presented, and example results are shown for both single sensors and multiple sensors in a test beam telescope configuration. The studies shown are those typically performed on sensor prototypes in test beam campaigns, and a comparison is made to test beam data, showing a maximum deviation of 4% and demonstrating that the approach is viable for generating realistic results.