Next Article in Journal
Maximum Acceptable Tilt Angle for Point Autofocus Microscopy
Previous Article in Journal
NOBEL-BOX: Development of a Low-Cost Ship-Based Instrument for Ocean Monitoring
Previous Article in Special Issue
Infrared HOT Photodetectors: Status and Outlook
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A 4H-SiC CMOS Oscillator-Based Temperature Sensor Operating from 298 K up to 573 K

1
Department of Industrial Engineering, University of Salerno, Via Giovanni Paolo II, 132, 84084 Fisciano, SA, Italy
2
Fraunhofer Institute for Integrated Systems and Device Technology (IISB), Schottkystraße 10, 91058 Erlangen, Germany
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Sensors 2023, 23(24), 9653; https://doi.org/10.3390/s23249653
Submission received: 11 October 2023 / Revised: 29 November 2023 / Accepted: 3 December 2023 / Published: 6 December 2023

Abstract

:
In this paper, we propose a temperature sensor based on a 4H-SiC CMOS oscillator circuit and that is able to operate in the temperature range between 298 K and 573 K. The circuit is developed on Fraunhofer IISB’s 2 μm 4H-SiC CMOS technology and is designed for a bias voltage of 20 V and an oscillation frequency of 90 kHz at room temperature. The possibility to relate the absolute temperature with the oscillation frequency is due to the temperature dependency of the threshold voltage and of the channel mobility of the transistors. An analytical model of the frequency-temperature dependency has been developed and is used as a starting point for the design of the circuit. Once the circuit has been designed, numerical simulations are performed with the Verilog-A BSIM4SiC model, which has been opportunely tuned on Fraunhofer IISB’s 2 μm 4H-SiC CMOS technology, and their results showed almost linear frequency-temperature characteristics with a coefficient of determination that was higher than 0.9681 for all of the bias conditions, whose maximum is 0.9992 at a V D D = 12.5 V. Moreover, we considered the effects of the fabrication process through a Monte Carlo analysis, where we varied the threshold voltage and the channel mobility with different values of the Gaussian distribution variance. For example, at V D D = 20 V, a deviation of 17.4 % from the nominal characteristic is obtained for a Gaussian distribution variance of 20 % . Finally, we applied the one-point calibration procedure, and temperature errors of +8.8 K and −5.8 K were observed at V D D = 15 V.

1. Introduction

4H-polytype silicon carbide semiconductor material is widely used for high temperature applications [1,2,3], and the possibility to fabricate integrated circuits (ICs) allows for an extension of its application fields. Up to now, both unipolar and bipolar 4H-SiC IC technologies have been developed: on bipolar technology, a Bipolar Junction Transistor based on multi-epitaxial stacks is used in analog [4,5,6] and digital [7] ICs, and its performance has been demonstrated up to 873 K. The 4H-SiC Complementary Metal Oxide Semiconductor Field Effect Transistor, CMOS technology proposed by Raytheon, has been developed, and analog and digital building blocks have been fabricated [8], such as, for example, a Positive-To-Absolute-Temperature (PTAT) circuit in the range between 298 K and 573 K, and with a maximum deviation from the ideal linear curve of 33 % [9]. Recently, Fraunhofer IISB provided a 4H-SiC 2 μ m-CMOS technology [10] and several ICs have been proposed, like CMOS Complementary-To-Absolute-Temperature (CTAT), a sensor in the range between 298 K and 438 K and with a sensitivity of 7.5 mV/K [11], or a temperature sensor based on a p-n diode from 297 K and 873 K, with an R 2 = 0.9998 of the voltage-temperature characteristic. Moreover, such technology is also compatible with other device structures that are useful for sensing temperature and ultraviolet radiation [12].
All the 4H-SiC-based proposed temperature sensors transduce temperature in electrical quantities, either through the difference of Gate Source Voltages ( Δ V G S ) between two MOSFETs [11], or through the difference of diode forward voltages [2,3,9,13,14]. However, 4H-SiC diodes with good performances have vertical structures and they are incompatible with VLSI circuits, whereas, although MOSFETs can be used, they need an integrated circuit in order to read out the voltage-temperature signal. Hence, the study of 4H-SiC CMOS circuits is a relevant topic, both for investigating their potentiality as read-out circuits for sensors and to propose new types of sensors. Among these last, temperature sensors based on time conversion can be a valid, fully compatible, 4H-SiC CMOS alternative, but they have been reported only in Si CMOS technology [15,16]. It is worth noting that the fabrication process quality is still poor because the 4H-SiC CMOS technology is at an early stage, and the variation in the fabrication process parameter on circuit performances is a mandatory analysis.
In this paper, a first all-CMOS temperature sensor based on an oscillator circuit is proposed. It converts the operating temperature into an oscillation frequency of triangular or square waveform voltages. The circuit is fully compatible with the 4H-SiC 2 μ m-CMOS technology and its performances are extracted from numerical simulation up to an operating temperature of 573 K. Considering that the oscillation frequency is related to the charge and discharge of an integrated capacitance through the device currents, we can obtain a temperature-frequency conversion because those currents have a temperature dependency. The paper is organized as follows: in Section 2, the topology, the operating principle, and the design of the sensor are reported; in Section 3, numerical simulation results are shown, focusing both on the differences from the analytical design due to second-order device phenomena and on the effects of the fabrication process variations; Section 4 and Section 5 are, respectively, the layout and the conclusions.

2. The Topology of the Temperature-Sensing Oscillator

The circuit of Figure 1 is the proposed sensor and is based on an oscillator whose frequency is uniquely related to the temperature. In the following, we report the single sub-circuits and the design of the circuit.
The circuit is an astable multivibrator and it is composed of a CMOS Schmitt trigger, three CMOS inverters, a n-type MOS capacitor, and a common drain NMOSFET amplifier. The circuit can generate a square waveform voltage, O U T 3 , thanks to the CMOS inverter I N V 3 , which is opportunely sized to load capacitance C 2 , and a triangular waveform voltage at the O U T 1 . The waveforms of O U T 1 and O U T 3 are related to each other thanks to the integrator composed of I N V 2 and NMOS Cap. The sensing circuit is biased to positive voltage, V D D , and a negative voltage, V S , which needs to bias the Class-A amplifier.

2.1. CMOS Voltage Schmitt Trigger

The 4H-SiC CMOS voltage Schmitt trigger schematic and transcharacteristic are reported, respectively, in Figure 2a,b [17,18]. It has three pairs of CMOS transistors, which are M N 1 3 and M P 1 3 , and the input and output signals are, respectively, I N and O U T 0 .
Until V I N is lower than V T H N 1 , which is the threshold voltage of the NMOSFET M N 1 , the stacked PMOS transistors M P 1 and M P 2 are in linear operation region, and the source follower M N 3 is in saturation operation mode, whereas M P 3 , M N 1 , and M N 2 are turned off. The drain potential of the M N 1 is V D D V T H N 3 until V I N remains lower than V T H N 1 . When V I N V T H N 1 , M N 1 is turned on and biased in the saturation region. When V I N is higher than V D S N 1 + V T H N 2 , the transistor M N 2 turns on so that the positive feedback loop, composed of M N 2 and M N 3 , quickly drops the O U T 0 voltage to ground [19]. Meanwhile, the source follower M P 3 turns on and brings the source potential of M P 2 to a low value so that this last one turns off. The high threshold voltage, V T + , can be calculated as follows [17]:
V I N = V T + = V D D | V T H P | + K N K P V T H N 1 + K N K P + K N K P ( V D D V T H N ) ( 1 + W N 1 L N 3 W N 3 L N 1 ) ( 1 + K N K P )
where K N = K N 1 / / K N 2 and K P = K P 1 / / K P 2 , K N i = K N W N i / L N i , K P i = K P W P i / L P i . Also, K N = μ N C O X and K P = μ P C O X are, respectively, the transconductance coefficients of the NMOSFET and PMOSFET, μ N ( P ) is the carrier channel mobility, C O X is the gate oxide capacitance, and W i and L i are, respectively, the width and the length of the i-transistor. Moreover, at V I N = V T + and before the positive feedback loop intervenes, it is reasonable that V O U T 0 = V D D and that the following relation is valid [19]:
K N 1 K N 3 = V D D V T + V T + V T H N 2
In the reverse direction, i.e., where V I N reduces from V D D to ground, the symmetrical process takes place. At V I N = V D D , the transistors M N 1 , M N 2 , and M P 3 conduct, while M N 3 , M P 1 , and M P 2 are turned off. Then, at V I N V D D | V T H P | , the transistor M P 1 conducts, and when V I N becomes lower than V D D V S D P 1 | V T H P 2 | , the transistor M P 2 turns on and the positive feedback loop, composed by M P 2 and M P 3 , is closed and brings the output voltage to V D D . Simultaneously, the source follower M N 3 turns on and increases the source potential of M N 2 , inducing its turn-off. The low threshold voltage, V T , can be calculated using the following equation [17]:
V I N = V T = V D D | V T H P | + K N K P V T H N 1 + K N K P V D D | V T H P | ( 1 + W P 1 L P 3 W P 3 L P 1 ) ( 1 + K N K P )
Moreover, at V I N = V T and before the positive feedback loop intervenes, it is reasonable that V O U T 0 = 0 and the following relation is valid [19]:
K P 1 K P 3 = V T V D D V T | V T H P | 2
By defining the values of V D D , V T + , and V T , one can design the channel sizes of M N 1 and M N 3 from (2), and of M P 1 and M P 3 from (4); instead, M N 2 and M P 2 are obtained, respectively, from (1) and (3).
In the proposed circuit of Figure 1, the O U T 0 is followed by a CMOS inverter, I N V 1 , composed of M N 4 - M P 4 , in order to invert the transcharacteristic and to make it compatible with the second part of the circuit. Moreover, I N V 1 drives the loads connected to the Schmitt trigger so that the dynamic behavior is preserved, as shown in the Section 2.2, Section 2.3 and Section 2.4.

2.2. Integrator and Output Stages

The NMOS capacitor, M N , C A P , of Figure 1, is an NMOS-based capacitor and it has been preferred to MIM-caps due to the lack of a Verilog-A model of this last in our actual 4H-SiC CMOS technology. Its gate is controlled by the output of the Schmitt trigger and is driven by the CMOS inverter I N V 2 , composed by M N 5 - M P 5 , with a constant current so that a triangular waveform appears at the input of the Schmitt trigger, i.e., V O U T 1 .
The ramp time is defined by the charging and discharging of M N C A P : when V G N , C A P increases toward V D D , starting from V T , the time constant is τ P = R P C L and the charging stops at V G N , C A P = V T + because the Schmitt trigger changes state; when V G N , C A P decreases toward ground starting from V T + , the time constant is τ n = R N C L and the discharging stops at V G N , C A P = V T . The capacitance C L I is the total load capacitance at the O U T 1 terminal and the equivalent resistances of the time constants are [20]:
R N = L N μ N C O X W N ( V D D V T H N )
R P = L P μ P C O X W P ( V D D | V T H P | )
Moreover, to avoid frequency and waveform distortions, we added two output stages needed to charge the probe capacitances C 1 and C 2 , which are of the order of tens of p F . Observing the circuit of the Figure 1, one is a Class-A amplifier based on NMOS M N 7 for the triangular waveform output, V O U T 2 , and the other is a CMOS inverter I N V 3 , composed of M N 6 - M P 6 , for the square waveform output, V O U T 3 . The V O U T 2 is the same waveform of V O U T 1 with a level-shift voltage of V G S N 7 , and hence, the maximum and minimum voltages of V O U T 2 are V + = V T + V G S + and V = V T V G S 7 + , respectively, where V G S 7 + and V G S 7 are the V G S N 7 when V G 7 is equal, respectively, to the V T + and V T values. Moreover, upon observing Figure 1, the Class-A amplifier is biased by a resistor, R S and a bias voltage, V S , which are externally applied to the O U T 2 terminal.

2.3. Evaluation of the Oscillation Frequency

The oscillation frequency, f O S C , can be calculated considering the charge and discharge times of M N C A P , t C A P , together with the time propagation delays of the Schmitt trigger and of I N V 1 . Indeed, one obtains:
f O S C = 1 2 t C A P + t p H L t r i g g e r + t p L H t r i g g e r + t p H L I N V 1 + t p L H I N V 1
where t p H L ( L H ) t r i g g e r and t p H L ( L H ) I N V 1 are the time propagation delays high-low (low-high) of the Schmitt trigger and of I N V 1 , respectively, and t C A P is the mean time of the charge and discharge time of M N C A P . The latter has the following relation in order to design the circuit:
t C A P = t c h a r g e = t d i s c h a r g e
Therefore, we obtain the same time to charge and discharge M N , C A P , which means a saw-tooth waveform signal for V O U T 1 as well as for V O U T 2 , and a duty cycle of 50 % for the square waveform of V O U T 3 . To calculate f O S C , the single terms of (6) are evaluated as follows [20]:
t c h a r g e = R P 5 C L I l n V D D V T V D D V T +
t d i s c h a r g e = R N 5 C L I l n V T + V T
t p H L t r i g g e r = 0.69 ( R N 1 + R N 2 ) C L I I
t p L H t r i g g e r = 0.69 ( R P 1 + R P 2 ) C L I I
t p H L I N V 1 = 0.69 R N 4 C L I I I
t p L H I N V 1 = 0.69 R P 4 C L I I I
where the capacitances are as follows:
C L I = C N C A P + 2 C G D N 5 + 2 C G D P 5 + C D B N 5 + C D B P 5 + C G N 1 + C G N 2 + C G P 1 + C G P 2 + C G N 7
C L I I = C G N 4 + C G P 4 + 2 C G D N 2 + 2 C G D P 2 + C D B N 2 + C D B P 2 + C G N 3 + C G P 3
C L I I I = C G N 5 + C G P 5 + C G N 6 + C G P 6 + 2 C G D N 4 + 2 C G D P 4 + C D B N 4 + C D B P 4
where C G , C G D , and C D B are, respectively, the Gate, the Gate-Drain, and the Drain-Body capacitance.
The design can be simplified, assuming the following conditions:
  • for C N C A P that is higher than other capacitances, one has C L I C N C A P ;
  • for t C A P that is higher than the propagation delays, one obtains f O S C 0.5 t C A P 1 .
About the O U T 3 signal, I N V 3 loads the probe capacitance C 2 and the time propagation delay is evaluated as follows:
t p H L I N V 3 = 0.69 R N 6 C 2
t p L H I N V 3 = 0.69 R P 6 C 2
Therefore, assuming that all the channel lengths are equal to the minimum one, and that the time propagation delays are negligible for t p H L ( L H ) = 10 3 t C A P , the design is completed in this way:
  • fixing f O S C , one has t C A P = t c h a r g e = t d i s c h a r g e = 0.5 f O S C 1 ;
  • fixing C N C A P , one obtains W P 5 from (8a) and W N 5 from (8b);
  • fixing t p I N V 3 and C 2 , the channel widths of I N V 3 are calculated from (12);
  • fixing t p I N V 1 , the channel widths of I N V 1 are calculated from (10);
  • fixing t p t r i g g e r , (9) complete the equations used together with (1)–(4) to design the Schmitt trigger.
Finally, because the M N 7 class-A output stage loads the probe capacitance C 1 in order that the triangular waveform is undistorted, its current has to respect the following relation (see Figure 1):
I C 3 = C 1 Δ V Δ t = C 1 V + V t C A P
Also, considering that M N 7 is in saturation operation mode and for a selected value of V S , it is possible to write:
{ 2 V + V S R S I C 3 μ N C o x W 7 L 7 + V T H N V G S 7 + = 0 2 V V S R S I C 3 μ N C o x W 7 L 7 + V T H N V G S 7 = 0
Hence, the values of W N 7 , L N 7 , and R s can be found from (13)–(14).

2.4. Design of the Circuit

The equations from (1) to (14) are used to design our circuit with the project specifications of Table 1. The values are from a trade-off among the expected performances of the 4H-SiC CMOS technology and the possibility to make our sensor compatible with the current electronic and the occupied wafer area. The device physical parameters used in (1)–(14) are the threshold voltages of the MOSFETs, V T H N = 5.8 V and V T H P = 8 V, the channel mobility, μ N = 17.14 cm 2 V 1 s 1 and μ P = 3.52 cm 2 V 1 s 1 , and the oxide capacitance, C O X = 62.78 nFcm 2 . In the Appendix A, the device physical parameters, characteristics, and extraction procedure are reported.
By using the results of the previous subsections, the design of the proposed circuit is reported in Table 2. The proposed procedure is based on an analytical approach and allows for the definition of a first-order design of the circuit, which can be used as a starting point for the numerical simulations. Indeed, a tuning procedure is required to achieve the project specifications, as has been shown in the following section, due to second-order effects of the transistors, like 4 H S i C / S i O 2 interface defects or fabrication process non-uniformity.

3. Numerical Simulation Results and Process Variability

Once the design has been completed, numerical simulations have been performed in the Cadence Virtuoso environment [21] by using a Verilog-A BSIM MOSFETs model [22] whose parameters are opportunely tuned to fit the experimental curves of the 4H-SiC MOSFETs in the temperature range from 298 K to 573 K. The model has been developed by Fraunhofer IISB, and a more detailed description of the 4H-SiC MOSFETs is reported in Appendix A.
Numerical results report some inconsistencies with the project specifications; for example, the oscillation frequency is 85.94 kHz instead of 90 kHz, or the trigger threshold voltages are V T + = 8.14 V and V T = 4.4 V compared to 10 V and 5 V, respectively, as reported in Table 1. Keeping the same L N and L P of Table 2 but varying W N and W P , we obtained the design reported in Table 3, and the circuit results show f O S C = 93.8 kHz, V T + = 10.66 V, and V T = 5.37 V. Indeed, in Figure 3, the comparison between the two designs shows how the asymmetry of the trans-characteristics obtained using the analytical design approach, disappears in the one resulting from the tuning of the transistor sizes. Such differences can be explained by second-order effects, which are neglected in the simplified current model used to extract (1)–(14): for example, they are the saturation of the carrier velocity, the channel length modulation, the bias effect of the body, and the effects of defects at the 4 H S i C / S i O 2 interface on the MOSFET electric behavior [23]. For example, observing Figure A3b, the high density of interface defects modifies the NMOSFET channel mobility dependency on V G S compared to the typical shape in Silicon technology: indeed, in Silicon technology, we expect a step-like curve at V G S V T H N , and then, a slight decay of μ N for V G S > V T H N [24]; instead, in 4H-SiC technology, the mobility has a continuous increase with V G S , and for V G S > 15 V, it remains constant.
As our temperature sensor is an oscillator, the total harmonic distortions, THDs, of V O U T 2 and V O U T 3 for Table 3 have been evaluated. Indeed, our sensor can be followed by a read-out circuit, like a frequency counter or microcontroller, whose aim is to measure the frequency; such a measurement is as accurate as the waveform is not distorted. These respectively show triangular and square waveforms in Figure 4 at T = 298 K, with a THD, respectively, of 9.89 % and of 33.18 % , which are slightly lower than the pure symmetric waveforms, i.e., T H D = 12.1 % for the triangular and T H D = 48.3 % for the squared waveforms [25].
Moreover, an estimation of the total average power dissipation gives a value of 2.45 mW at 298 K and the bias-frequency sensitivity is 9.28 kHz/V, as shown in Figure 5.
In the following, the effects on the oscillation frequency of the temperature and of the fabrication process variations are analyzed and investigated for the design reported in Table 3.

3.1. Oscillation Frequency Dependency on the Temperature

The oscillation frequency dependency on the temperature, f O S C -T, can be shown, beginning with the following relation:
f O S C 1 t c h a r g e + t d i s c h a r g e = 1 C L I R P 5 l n V D D V T V D D V T + + R N 5 l n V T + V T
where it is assumed that t C A P ( t C M O S I N V I + t t r i g g e r ) in (6). The dependence of the MOSFETs on the temperature is described through the channel mobilities and the threshold voltages, which are explicitly reported in the Appendix A. However, for a first analysis of (15), only the channel mobility is considered because V T H N and V T H P (i) appear in V T and V T + , which are divided by V D D and they are also the arguments of a logarithm function; and (ii) are in R P 5 and R N 5 , which are divided by V D D . Hence, one obtains:
f O S C = 1 C L I 1 l n V D D V T V D D + V T + μ P ( T 0 ) C o x W P 5 L P 5 V D D 1 | V T P | V D D + l n V T + V T μ N ( T 0 ) C o x W N 5 L N 5 V D D 1 V T N V D D T T 0 α = A T T 0 α
where A includes all the parameters related to the fabrication process and is independent from the temperature, in a first approximation. In Figure 6a, the resulting f O S C -T curve at V D D = 20 V is reported and the best fitting of (16) gives A = 98.66 kHz and α = 1.22 with an R 2 = 0.9656 . We also investigated the effects of V D D on the linearity of the f O S C -T curve, varying from 12.5 V to 20 V, and an improvement is obtained at V D D = 12.5 V with an R 2 = 0.9992 , as shown in Figure 6b. The f O S C -T curve at V D D = 12.5 V is reported in Figure 6a and it is compared with the model having A = 23.88 kHz and α = 1.883 , where a better fitting is evident; however, as is expected from (16), there is a reduction of the value of f O S C from 93.8 kHz to 23.88 kHz, evaluated at T = 298 K. For completeness in Figure 6a, we report the f O S C -T characteristic for the bias voltage of 15 V and we obtain R 2 = 0.9928 , A = 45.92 kHz, and α = 1.64 .

3.2. Effects of Process Parameters Variation

Fabrication process variations are expected, and for example, in 4H-SiC CMOS technology, they can be related to variations on the activation process of the Aluminum p-type doping atoms [26], on the uniformity of the doping concentration and of the oxide thickness, or on the quality of the contact resistance of doped regions [27]. All these fabrication process variations can be modeled, in a first analysis, on variations in the channel mobility and of the threshold voltage, and then used for Monte Carlo analysis to assess the sensitivity of the circuit. In particular, the analysis is a 1000 point process Monte Carlo and consists of evaluating the f O S C -T curves, considering a Gaussian distribution for μ N ( P ) and V T H , N ( P ) , which is defined by a standard deviation ( σ ) and a mean value ( μ ) through the following distribution:
p ( x ) = 1 σ 2 π e ( x μ ) 2 2 σ 2
During the analysis, we varied the supply voltage from 12.5 V to 20 V and the ratio σ / μ of the Gaussian distribution from ± 10 % to ± 20 % , either for both parameters or singularly. To compare the cases, we use the oscillation frequency variation, f O S C , v a r , defined as follows:
f O S C , v a r ( T ) = f O S C , σ ( T ) f O S C ( T ) f O S C ( T ) 100
In Figure 7a–c, the f O S C -T curves for a σ / μ = ± 0.1 and at different values of V D D are reported. Although the curve of Figure 7a for the case V D D = 12.5 V shows a better linearity, a higher f O S C , v a r appears: indeed, it is almost 23 % , and in Figure 7d, the f O S C , v a r -T curve is shown. Instead, observing Figure 7c for V D D = 20 V, a maximum f O S C , v a r of 8 % is achieved, but a worse linearity is obtained, as shown in Figure 6. Hence, a supply voltage of 15 V allows for a good trade-off between the process parameter variation and the linearity, as is clearly shown in Figure 7b,d.
To understand the effects of the single device parameter on the performance of the circuit, we performed Monte Carlo analysis by singularly varying either V T H , N ( P ) or μ N ( P ) . In Figure 8, a σ / μ of ± 10 % for V T H , N ( P ) with a constant value of μ N ( P ) at the nominal value reported in Table 1 shows a stronger variation in the oscillation frequency as a function of the temperature than the case of a σ / μ of ± 10 % for μ N ( P ) with a constant value of V T H , N ( P ) , as reported in Figure 9. Indeed, observing Figure 8d, the maximum f O S C , v a r related to the variations in V T H , N ( P ) are 22.64 % , 11.76 % , and 7.24 % for V D D , equal to 12.5 V, 15 V, and 20 V, respectively, whereas for the variations in μ N ( P ) , they are around 5.8 % for all bias conditions (see Figure 9d). That is also confirmed, considering that the maximum variation for the cases of Figure 7 is almost similar to that of Figure 8, and in Table 4, they are reported for ease of reading. However, the best trade-off in terms of the linearity and maximum variation is still at V D D = 15 V for both cases.
To stress the effect of process variations on the circuit performance, we increase σ / μ to ± 15 % and ± 20 % for V D D = 20 V, this being the bias condition defined during the design specification (see Table 1). In Figure 10, the results for a σ / μ of both parameters indicate an expected increase in the divergence from the nominal value, and observing Figure 10c, the percentage of variation of f O S C from the nominal value decreases with the temperature, but it has a peak at around T = 400 K and a further increase over 500 K. To understand the behavior, we separately analyzed the σ / μ of the parameters, and the results are reported in Figure 11 and in Figure 12, respectively, for μ N ( P ) and V T H , N ( P ) . It is clear that the variation reduces for the μ N ( P ) -case with the increase in the temperature, whereas the V T H , N ( P ) -case is almost constant, except for a maximum at T = 400 K, which is related to that of Figure 10c. Moreover, both cases have an increase for temperatures that are higher than 500 K. In Table 5, we summarized the maximum variation for f O S C , and the greater effect of the process variations for the V T H , N ( P ) -case compared to the μ N ( P ) -case is clear.
Considering that the best trade-off between linearity and circuit specification is at V D D = 15 V, we also performed Monte Carlo analysis for this case by increasing the σ / μ , and in Figure 13, the results are reported. Observing the curves, the 15 V-case has an increased variation compared with V D D = 20 V, which is more evident for T > 400 K. It can be summarized through the f O S C , v a r -T curve of Figure 13c, with a maximum variation of 27.11 % at 425 K for a σ / μ = 0.2 . Instead, the deviation reduces in the high temperature range, contrary to the 20 V-case. Moreover, the separation of the effects of the process variation of μ N ( P ) and V T H , N ( P ) have been analyzed and shown, respectively, in Figure 14 and in Figure 15: the prominent effect of V T H , N ( P ) is evident with respect to μ N ( P ) , which has at least double the value and defines the behavior of Figure 13c. It is interesting to note that the increase of the f O S C , v a r at T > 500 K for V D D = 20 V disappears for the case of V D D = 15 V. In Table 5, we also report the maximum values of f O S C , v a r at V D D = 15 V, and higher values are shown than the 20V-case, so that although the 15V-case has a higher linearity, it is more greatly affected by process variations.

3.3. Results after Sensor Calibrations

If, on one hand, the bias voltage of V D D = 20 V reduces the dependency on the fabrication process of the oscillation frequency, on the other hand, the f O S C -T characteristic has a scarce linearity. To overcome it, (16) suggests that it is possible to apply the one-point calibration procedure [16], where f O S C -T curve is normalized by a f O S C , 1 p o i n t selected at a fixed temperature, T 1 p o i n t , as follows:
f O S C , n o r m ( T ) = f O S C ( T ) f O S C ( T 1 p o i n t ) = f O S C ( T ) f O S C , 1 p o i n t
In this way, the negative effect of the process variation can be partially eliminated. Indeed, we applied it in the case of σ / μ = ± 0.1 , both for V T H , N ( P ) and for μ N ( P ) , and in Figure 16, the results are shown for different bias conditions and T 1 p o i n t . For V D D = 12.5 V of Figure 16a, a better correction is for T 1 p o i n t = 473 K compared to the room temperature case, as well as for V D D = 15 V of Figure 16b, where T 1 p o i n t = 423 K; instead, for V D D = 20 V, a T 1 p o i n t = 298 K can be used (see Figure 16c). This last result makes the bias supply of V D D = 20 V advantageous because an easier calibration procedure is applicable. Indeed, in term of the calibration easiness, the normalization at T 1 p o i n t = 298 K is of great advantage, but in the case of V D D = 12.5 V, it has a maximum variation of almost 8.65 % from the nominal f O S C , n o r m ; instead, the same case evaluated at T 1 p o i n t = 423 K has a maximum variation of 5.93 % .
Finally, once the calibration has been done, one can extract temperature from (16) using the measured frequency as the input variable. Hence, for the case of V D D = 15 V, one can use the parameters A = 1.02 and α = 1.64 , and the calibration at T 1 p o i n t = 298 K reported in Figure 16b. In these conditions, the error between the extracted and effective temperatures, T E R R , is shown in Figure 17a, and its maximum absolute value is 8.8 K across the whole temperature range between 298 K and 573 K. Furthermore, the process variations of σ / μ = ± 0.1 , singularly either for V T H , N ( P ) or for μ N ( P ) , give errors of 8.8 K and 6.58 K, respectively. For the 20 V-case of Figure 17c, a maximum T E R R of 11.25 K has been found, both for the nominal case and for a process variation of σ / μ = ± 0.1 . The error T E R R can be reduced if a third-order curve is used as model [16], and after a similar analysis, we obtained the results reported in Figure 17b for V D D = 15 V and in Figure 17d for V D D = 20 V, where T E R R reduces, respectively, at 5 K/4.78 K and 8.15 K/4.49 K.

3.4. Comparisons with the State-of-the-Art

An exhaustive comparison between the state-of-the-art and our proposal is difficult due to the limited availability of 4H-SiC temperature-frequency converter sensors. However, in Table 6, we report sensors based on similar operating principles and on other technology in order to understand how our proposal improves the state-of-the-art. The circuit of [16] is a Silicon CMOS technology circuit based on delay-locked loops and it shows an error of between + 4 K and 4 K, which is slightly lower than ours, but in a narrower temperature range, i.e., [ 273.15 ; 373.15 ] K. Indeed, as an example from a comparison with [28], our proposal has a greater temperature error, i.e., around 10 % , but it can operate within a wider temperature range, even to 200 K, thanks to the higher performance of 4H-SiC technology compared to the Silicon one. Then, the proposals of [29,30], based on an Si 65nm-CMOS technology, use an area that is lower than ours, and in particular, [29] have 83 % power dissipation and 66 % temperature error, whereas [30] has, respectively, 44 % and 83 % compared to ours. Anyway, the dynamic power can be justified with the reduced channel length; instead, the error is lower because both the temperature variations of the Si transistor parameters are reduced due to the better S i / S i O 2 interface quality, and the temperature range is only 100 K compared to the value of 275 K for our sensor. In order to stress the comparison, we reduced the temperature range between 298 K and 443 K, extracting a T E R R of 3.15 K/4.58 K, which reduces to 3 K/3.84 K if a third-order model is used.
In Table 6, we report time-based temperature sensors fabricated using Silicon-On-Insulator technology. Ref. [31] is a Full-Depleted SOI at 28 nm CMOS and it has a very low T E R R as well as power consumption, and similarly, the SOI 32nm CMOS temperature sensors of [32] have a T E R R = ± 1.95 K. However, in both cases, although SOI-CMOS is much more mature technology than the 4H-SiC CMOS one, the temperature range is up to 385 K, which is 188 K lower than our proposal.

4. Layout

The design results of Table 3 are used to draw the final layout of the circuit reported in Figure 18 with the Cadence Virtuoso 6.1.8 Layout software. The 4H-SiC 2 μ m-CMOS process has 14 masks and two metal layers, whereas the active area of the sensor of Figure 1 is 0.163 mm 2 .

5. Conclusions

In this paper, the design of a 4H-SiC CMOS temperature sensor based on an analog oscillator is presented, as well as an analysis of its performances in terms of the process fabrication variations and the bias voltage. Unlike Silicon technology, our circuit showed a reduction in the propagation delay with an increase in the temperature, because the channel mobility improves at high temperature, and consequently, the frequency increases. The relation between the oscillation frequency and the temperature is almost linear, with an R 2 = 0.9992 at V D D = 12.5 V, but with a smaller influence of the process fabrication variation, i.e., 17.4 % for σ / μ = ± 20 % , is at V D D = 20 V, with an R 2 = 0.9681 . On the other hand, performing a one-point calibration, we obtain a temperature error of + 8.8 K and 5.8 K when V D D = 15 V.
The use of a simple model for the MOSFET current results are useful for a first-order approximation of the circuit design, but the effects of a high defects density at the S i O 2 / 4 H S i C interface should be considered for a better description of the circuit. It means that a more accurate model has to be developed and the μ N ( P ) -dependency on the temperature should be focused on.

Author Contributions

Conceptualization, N.R. and L.D.B.; methodology, N.R. and L.D.B.; software, N.R., C.R. and A.M.; validation, N.R., R.L., L.D.B., C.R. and A.M.; formal analysis, L.D.B.; investigation, N.R., R.L. and L.D.B.; resources, L.D.B. and M.R.; data curation, R.L., L.D.B. and M.R.; writing—original draft preparation, N.R., R.L., A.M., M.R., G.D.L. and L.D.B.; writing—review and editing, N.R., L.D.B., R.L., A.M., C.R. and M.R.; visualization, N.R., L.D.B. and R.L.; supervision, L.D.B., M.R. and G.D.L.; project administration, A.R. and M.R.; funding acquisition, A.R. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

A numerical simulation model has been developed and provided by Fraunhofer IISB. Please contact [email protected] for further information.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

4H-SiC 2 μ m-CMOS technology has been developed by Fraunhofer IISB (Germany), and the main physical parameters and the electrical parameters evaluated at room temperature are reported in Table A1.
Table A1. Physical and electrical parameters of the 4H-SiC CMOS technology at r.t.
Table A1. Physical and electrical parameters of the 4H-SiC CMOS technology at r.t.
ParameterUnitValue
t O X [nm]55
ϵ O X [fFcm 1 ]345.31
ϵ S i C [fFcm 1 ]855.30
C O X [nFcm 2 ]62.78
V T H N [V]5.8
V T H P [V]−8
μ N [cm 2 V 1 s 1 ]17.14
μ P [cm 2 V 1 s 1 ]3.52
The experimental transcharacteristics of Figure A1 and the output characteristics of Figure A2 are for PMOSFET and NMOSFET, with a form factor of 100 μ m/10 μ m and 10 μ m/10 μ m, respectively.
Figure A1. Comparison between the experimental and Verilog-A numerical model transcharacteristics of lateral 4H-SiC. (a) PMOSFET with W P / L P = 100 μ m/10 μ m and (b) NMOSFET with W N / L N = 10 μ m/10 μ m. All the curves are at r.t., | V D S | from 0.1 V to 20 V, and | V B S | = 0 V.
Figure A1. Comparison between the experimental and Verilog-A numerical model transcharacteristics of lateral 4H-SiC. (a) PMOSFET with W P / L P = 100 μ m/10 μ m and (b) NMOSFET with W N / L N = 10 μ m/10 μ m. All the curves are at r.t., | V D S | from 0.1 V to 20 V, and | V B S | = 0 V.
Sensors 23 09653 g0a1
Figure A2. Comparison between experimental and Verilog-A numerical model output characteristics of lateral 4H-SiC. (a) PMOSFET with W P / L P = 100 μ m/10 μ m and (b) NMOSFET with W P / L P = 10 μ m/10 μ m. All the curves are at r.t. and | V B S | = 0 V.
Figure A2. Comparison between experimental and Verilog-A numerical model output characteristics of lateral 4H-SiC. (a) PMOSFET with W P / L P = 100 μ m/10 μ m and (b) NMOSFET with W P / L P = 10 μ m/10 μ m. All the curves are at r.t. and | V B S | = 0 V.
Sensors 23 09653 g0a2
The curves of Figure A1 at | V D S | = 100 mV are used to extract both the threshold voltage, V T H , and the mobility, μ N ( P ) . By using the second derivative method [33], we extracted the threshold voltages reported in Table A1, whereas the channel mobility dependency on the V G S is obtained using the following equation:
μ = I D S V G S | V D S = 100 mV W L C o x V D S
and are reported in Figure A3a,b for PMOSFET and NMOSFET, respectively. Their mean value is the channel mobility reported in Table A1 and it is used to design the proposed circuit in Section 2. In Figure A3b, the unusual behavior of the NMOSFET channel mobility as a function of V G S is shown, and it has been described in Section 3 in order to justify the inconsistency between the expected specification and the numerical simulation results of the design in Table 2. Indeed, only PMOSFET channel mobility (see Figure A3a) shows typical step-like behavior, as well as a decrease with the increase in V G S ; instead, μ N has a continuous increase until it reaches a maximum at around V G S = 15 V.
Figure A3. Extracted mobility from the transcharacteristics of Figure A1 by using (A1) for 4H-SiC. (a) PMOSFET with W P / L P = 100 μ m/10 μ m and (b) NMOSFET with W P / L P = 10 μ m/10 μ m at | V D S | = 0.1 V, | V B S | = 0 V, and at r.t.
Figure A3. Extracted mobility from the transcharacteristics of Figure A1 by using (A1) for 4H-SiC. (a) PMOSFET with W P / L P = 100 μ m/10 μ m and (b) NMOSFET with W P / L P = 10 μ m/10 μ m at | V D S | = 0.1 V, | V B S | = 0 V, and at r.t.
Sensors 23 09653 g0a3
The temperature effect on the transistor currents are reported in Figure A4 for a | V D S | = 20 V, and it is clear that the increase in the current with the temperature occurs both for PMOSFET and NMOSFET. Such a behavior is due to a combination of phenomena: on one hand, the channel mobility increases up to 400 K and then decreases, as reported in [34] for a 4H-SiC lateral NMOS with a bulk doped N A = 5 × 10 17 cm 3 ; on the other hand, the threshold voltage continuously decreases with the temperature [35]. In order to take into account such behavior, a temperature dependency of V T H can be modeled with the following linear function [24]:
V T H ( T 0 ) V T H ( T 0 ) [ 1 + β ( T T 0 ) ]
where β is a negative constant and we can neglect it from our analysis in a first approximation, because it is of the order of mV/K; on the other side, the temperature dependency of the channel mobility is modeled as follows [24,34]:
μ μ ( T 0 ) T T 0 α
where α is a positive fitting parameter and is higher than 1.
It is worth noting that the high density of defects at the 4 H S i C / S i O 2 interface is the reason for the improvement in the MOSFETs’ performance with the increase in temperature [36], whereas 4H-SiC bipolar devices showed worst behaviors under high-temperature operations; in particular, a thermal runaway of the current is observed [37,38], limiting their applications under harsh environment conditions.
Figure A4. Comparison between experimental and Verilog-A numerical model transcharacteristics of lateral 4H-SiC. (a) PMOSFET with W P / L P = 100 μ m/10 μ m and (b) NMOSFET with W N / L N = 10 μ m/10 μ m at | V D S | = 20 V, | V B S | = 0 V, and for temperature, from 298 K to 573 K.
Figure A4. Comparison between experimental and Verilog-A numerical model transcharacteristics of lateral 4H-SiC. (a) PMOSFET with W P / L P = 100 μ m/10 μ m and (b) NMOSFET with W N / L N = 10 μ m/10 μ m at | V D S | = 20 V, | V B S | = 0 V, and for temperature, from 298 K to 573 K.
Sensors 23 09653 g0a4
The numerical simulations performed in Cadence Virtuoso are based on the compact BSIM4SiC model [22] and they have been opportunely developed by Fraunhofer IISB for the IISB’s 2 μ m 4H-SiC CMOS technology. It is a extension of the widely used BSIM4 compact model, in which the typical physical phenomena of the 4H-SiC lateral n-channel and p-channel MOSFETs are considered. Indeed, the lower quality of the 4 H S i C / S i O 2 interface induces a high density of interface trapped charges, differing significantly, in terms of the transfer and output characteristics, from the Silicon MOSFETs ones.
The effects include mobility degradation, a lower subthreshold slope, soft saturation, and an increase in the flat band and threshold voltages [22]. For an accurate description of the IISB’s 2 μ m 4H-SiC CMOS technology, the parameters of the BSIM4SiC model have been determined through a parameter extraction procedure based on the measurement data of n- and p-type transistors. The developed compact model fits the experimental curves quite well, as shown in Figure A1 and Figure A2, with a higher accuracy for the NMOSFET case with respect to the PMOSFET case, especially for the output characteristics displayed in Figure A2a. The modeling of 4H-SiC PMOSFET is generally more challenging, due to the more limited knowledge about the defects at the 4 H S i C / S i O 2 interface in the n-type body region with respect to 4H-SiC NMOSFET, which has had much more interest thanks to the power electronics applications. Then, the model also includes channel geometry scaling, which makes the parameter extraction and the achievement of satisfactory results even more complicated for all of the considered geometry. Moreover, in Figure A4, its validity is shown over a temperature range from 298 K to 573 K.

References

  1. Kimoto, T.; Cooper, J.A. Fundamentals of Silicon Carbide Technology: Growth, Characterization, Devices and Applications; John Wiley & Sons Singapore Pte. Ltd.: Singapore, 2014. [Google Scholar]
  2. Zhang, N.; Lin, C.M.; Senesky, D.G.; Pisano, A.P. Temperature sensor based on 4H-silicon carbide pn diode operational from 20 °C to 600 °C. Appl. Phys. Lett. 2014, 104, 073504. [Google Scholar] [CrossRef]
  3. Matthus, C.D.; Di Benedetto, L.; Kocher, M.; Bauer, A.J.; Licciardo, G.D.; Rubino, A.; Erlbacher, T. Feasibility of 4H-SiC pin diode for sensitive temperature measurements between 20.5 K and 802 K. IEEE Sens. J. 2019, 19, 2871–2878. [Google Scholar] [CrossRef]
  4. Hussain, M.W.; Elahipanah, H.; Zumbro, J.E.; Rodriguez, S.; Malm, B.G.; Mantooth, H.A.; Rusu, A. A SiC BJT-based negative resistance oscillator for high-temperature applications. IEEE J. Electron Devices Soc. 2019, 7, 191–195. [Google Scholar] [CrossRef]
  5. Hedayati, R.; Lanni, L.; Rusu, A.; Zetterling, C.M. Wide temperature range integrated bandgap voltage references in 4H–SiC. IEEE Electron Device Lett. 2015, 37, 146–149. [Google Scholar] [CrossRef]
  6. Hedayati, R.; Lanni, L.; Rodriguez, S.; Malm, B.G.; Rusu, A.; Zetterling, C.M. A monolithic, 500 C operational amplifier in 4H-SiC bipolar technology. IEEE Electron Device Lett. 2014, 35, 693–695. [Google Scholar]
  7. Lanni, L.; Malm, B.G.; Östling, M.; Zetterling, C.M. ECL-based SiC logic circuits for extreme temperatures. Mater. Sci. Forum. 2015, 821–823, 910–913. [Google Scholar] [CrossRef]
  8. Rahman, A.; Shepherd, P.D.; Bhuyan, S.A.; Ahmed, S.; Akula, S.K.; Caley, L.; Mantooth, H.A.; Di, J.; Francis, A.M.; Holmes, J.A. A family of CMOS analog and mixed signal circuits in SiC for high temperature electronics. In Proceedings of the 2015 IEEE Aerospace Conference, Big Sky, MT, USA, 7–14 March 2015; IEEE: Piscataway, NJ, USA, 2015; pp. 1–10. [Google Scholar]
  9. Rahman, A.; Francis, A.M.; Ahmed, S.; Akula, S.K.; Holmes, J.; Mantooth, A. High-temperature voltage and current references in silicon carbide CMOS. IEEE Trans. Electron Devices 2016, 63, 2455–2461. [Google Scholar] [CrossRef]
  10. Fraunhofer Institute for Integrated Systems and Device Technology (IISB). Available online: https://www.iisb.fraunhofer.de/ (accessed on 30 September 2023).
  11. Romijn, J.; Middelburg, L.M.; Vollebregt, S.; El Mansouri, B.; van Zeijl, H.W.; May, A.; Erlbacher, T.; Zhang, G.; Sarro, P.M. Resistive and CTAT temperature sensors in a silicon carbide CMOS technology. In Proceedings of the 2021 IEEE Sensors, Sydney, Australia, 31 October–3 November 2021; IEEE: Piscataway, NJ, USA, 2021; pp. 1–4. [Google Scholar]
  12. Di Benedetto, L.; Licciardo, G.; Erlbacher, T.; Bauer, A.; Rubino, A. A 4H-SiC UV phototransistor with excellent optical gain based on controlled potential barrier. IEEE Trans. Electron Devices 2019, 67, 154–159. [Google Scholar] [CrossRef]
  13. Mo, J.; Li, J.; Zhang, Y.; Romijn, J.; May, A.; Erlbacher, T.; Zhang, G.; Vollebregt, S. A Highly Linear Temperature Sensor Operating up to 600 °C in a 4H-SiC CMOS Technology. IEEE Electron Device Lett. 2023, 44, 995–998. [Google Scholar] [CrossRef]
  14. Di Benedetto, L.; Licciardo, G.D.; Rao, S.; Pangallo, G.; Della Corte, F.G.; Alfredo, R. V2O5/4H-SiC Schottky diode temperature sensor: Experiments and model. IEEE Trans. Electron Devices 2018, 65, 687–694. [Google Scholar] [CrossRef]
  15. Tolić, I.P.; Mikulić, J.; Schatzberger, G.; Barić, A. Design of CMOS temperature sensors based on ring oscillators in 180-nm and 110-nm technology. In Proceedings of the 2020 43rd International Convention on Information, Communication and Electronic Technology (MIPRO), Opatija, Croatia, 28 September–2 October 2020; IEEE: Piscataway, NJ, USA, 2020; pp. 104–108. [Google Scholar]
  16. Ha, D.; Woo, K.; Meninger, S.; Xanthopoulos, T.; Crain, E.; Ham, D. Time-domain CMOS temperature sensors with dual delay-locked loops for microprocessor thermal monitoring. IEEE Trans. Very Large Scale Integr. (VLSI) Syst. 2011, 20, 1590–1601. [Google Scholar] [CrossRef]
  17. Dokic, B. CMOS Schmitt triggers. In Proceedings of the IEE Proceedings G. Electronic Circuits and Systems; 1984; Volume 131, pp. 197–202. Available online: https://digital-library.theiet.org/content/journals/10.1049/ip-g-1.1984.0037 (accessed on 5 September 2023).
  18. Jorgensen, J.M. CMOS Schmitt Trigger. U.S. Patent 3,984,703, 5 October 1976. [Google Scholar]
  19. Filanovsky, I.; Baltes, H. CMOS Schmitt trigger design. IEEE Trans. Circuits Syst. I Fundam. Theory Appl. 1994, 41, 46–49. [Google Scholar] [CrossRef]
  20. Rabaey, J.M.; Chandrakasan, A.P.; Nikolic, B. Digital Integrated Circuits; Prentice Hall: Englewood Cliffs, NJ, USA, 2002; Volume 2. [Google Scholar]
  21. Cadence Design Systems. Cadence Virtuoso. Available online: https://www.cadence.com/en_US/home.html (accessed on 30 September 2023).
  22. Ahmed, S.; Rashid, A.U.; Hossain, M.M.; Vrotsos, T.; Francis, A.M.; Mantooth, H.A. DC modeling and geometry scaling of SiC low-voltage MOSFETs for integrated circuit design. IEEE J. Emerg. Sel. Top. Power Electron. 2019, 7, 1574–1583. [Google Scholar] [CrossRef]
  23. Mudholkar, M.; Mantooth, H.A. Characterization and modeling of 4H-SiC lateral MOSFETs for integrated circuit design. IEEE Trans. Electron Devices 2013, 60, 1923–1930. [Google Scholar] [CrossRef]
  24. Arora, N. Mosfet Modeling for VLSI Simulation: Theory and Practice; World Scientific: Singapore, 2007. [Google Scholar]
  25. Blagouchine, I.V.; Moreau, E. Analytic method for the computation of the total harmonic distortion by the Cauchy method of residues. IEEE Trans. Commun. 2011, 59, 2478–2491. [Google Scholar] [CrossRef]
  26. Nipoti, R.; Di Benedetto, L.; Albonetti, C.; Bellone, S. Al+ implanted anode for 4H-SiC pin diodes. ECS Trans. 2013, 50, 391. [Google Scholar] [CrossRef]
  27. May, A.; Rommel, M.; Beuer, S.; Erlbacher, T. Via size-dependent properties of TiAl ohmic contacts on 4H-SiC. Mater. Sci. Forum 2022, 1062, 185–189. [Google Scholar] [CrossRef]
  28. Makinwa, K.A.; Snoeij, M.F. A CMOS Temperature-to-Frequency Converter with an Inaccuracy of Less than 0.5 °C (3σ) from −40 °C to 105 °C. IEee J. Solid-State Circuits 2006, 41, 2992–2997. [Google Scholar] [CrossRef]
  29. Anand, T.; Makinwa, K.A.; Hanumolu, P.K. A VCO based highly digital temperature sensor with 0.034 °C/mV supply sensitivity. IEEE J. Solid-State Circuits 2016, 51, 2651–2663. [Google Scholar] [CrossRef]
  30. Hwang, S.; Koo, J.; Kim, K.; Lee, H.; Kim, C. A 0.008 mm2 500 μW 469 kS/s Frequency-to-Digital Converter Based CMOS Temperature Sensor with Process Variation Compensation. IEEE Trans. Circuits Syst. I Regul. Pap. 2013, 60, 2241–2248. [Google Scholar] [CrossRef]
  31. Cochet, M.; Keller, B.; Clerc, S.; Abouzeid, F.; Cathelin, A.; Autran, J.L.; Roche, P.; Nikolić, B. A 225 μm2 Probe Single-Point Calibration Digital Temperature Sensor Using Body-Bias Adjustment in 28 nm FD-SOI CMOS. IEEE-Solid-State Circuits Lett. 2018, 1, 14–17. [Google Scholar] [CrossRef]
  32. Chowdhury, G.; Hassibi, A. An on-chip temperature sensor with a self-discharging diode in 32-nm SOI CMOS. IEEE Trans. Circuits Syst. II Express Briefs 2012, 59, 568–572. [Google Scholar] [CrossRef]
  33. Ortiz-Conde, A.; Sánchez, F.G.; Liou, J.J.; Cerdeira, A.; Estrada, M.; Yue, Y. A review of recent MOSFET threshold voltage extraction methods. Microelectron. Reliab. 2002, 42, 583–596. [Google Scholar] [CrossRef]
  34. Uhnevionak, V.; Burenkov, A.; Strenger, C.; Ortiz, G.; Bedel-Pereira, E.; Mortet, V.; Cristiano, F.; Bauer, A.J.; Pichler, P. Comprehensive study of the electron scattering mechanisms in 4H-SiC MOSFETs. IEEE Trans. Electron Devices 2015, 62, 2562–2570. [Google Scholar] [CrossRef]
  35. Rumyantsev, S.; Shur, M.; Levinshtein, M.; Ivanov, P.; Palmour, J.; Agarwal, A.; Hull, B.; Ryu, S.H. Channel mobility and on-resistance of vertical double implanted 4H-SiC MOSFETs at elevated temperatures. Semicond. Sci. Technol. 2009, 24, 075011. [Google Scholar] [CrossRef]
  36. Le-Huu, M.; Schmitt, H.; Noll, S.; Grieb, M.; Schrey, F.F.; Bauer, A.J.; Frey, L.; Ryssel, H. Investigation of the reliability of 4H–SiC MOS devices for high temperature applications. Microelectron. Reliab. 2011, 51, 1346–1350. [Google Scholar] [CrossRef]
  37. Di Benedetto, L.; Licciardo, G.D.; Nipoti, R.; Bellone, S. On the crossing-point of 4H-SiC power diodes characteristics. IEEE Electron Device Lett. 2013, 35, 244–246. [Google Scholar] [CrossRef]
  38. Li, X.; Luo, Y.; Fursin, L.; Zhao, J.; Pan, M.; Alexandrov, P.; Weiner, M. On the temperature coefficient of 4H-SiC BJT current gain. Solid-State Electron. 2003, 47, 233–239. [Google Scholar] [CrossRef]
Figure 1. The topology of the proposed temperature sensor, based on an oscillator circuit.
Figure 1. The topology of the proposed temperature sensor, based on an oscillator circuit.
Sensors 23 09653 g001
Figure 2. (a) CMOS voltage Schmitt trigger schematic and (b) its transcharacteristic.
Figure 2. (a) CMOS voltage Schmitt trigger schematic and (b) its transcharacteristic.
Sensors 23 09653 g002
Figure 3. Numerically simulated trans-characteristics of the Schmitt trigger from analytical and tuning designs at T = 298 K.
Figure 3. Numerically simulated trans-characteristics of the Schmitt trigger from analytical and tuning designs at T = 298 K.
Sensors 23 09653 g003
Figure 4. Waveform signals of (a) V O U T 2 and (b) of V O U T 3 at T = 298 K.
Figure 4. Waveform signals of (a) V O U T 2 and (b) of V O U T 3 at T = 298 K.
Sensors 23 09653 g004
Figure 5. Oscillation frequency, f O S C , dependency on power supply voltage, V D D , at T = 298 K, showing a sensitivity of 9.28 kHz/V.
Figure 5. Oscillation frequency, f O S C , dependency on power supply voltage, V D D , at T = 298 K, showing a sensitivity of 9.28 kHz/V.
Sensors 23 09653 g005
Figure 6. (a) Comparisons of the f O S C -T curves between numerical simulation results, obtained from the Verilog-A BSIM MOSFETs model and the analytical model of (16) at V D D = 12.5 V and 20 V. (b) R 2 V D D curve obtained from the f O S C -T curves.
Figure 6. (a) Comparisons of the f O S C -T curves between numerical simulation results, obtained from the Verilog-A BSIM MOSFETs model and the analytical model of (16) at V D D = 12.5 V and 20 V. (b) R 2 V D D curve obtained from the f O S C -T curves.
Sensors 23 09653 g006
Figure 7. Results of a 1000 points process Monte Carlo analysis for (a) V D D = 12.5 V, (b) V D D = 15 V, and (c) V D D = 20 V, with σ / μ = ± 0.1 for V T H , N ( P ) and for μ N ( P ) . (d) f O S C , v a r as function of the temperature in terms of V D D resulting from (ac).
Figure 7. Results of a 1000 points process Monte Carlo analysis for (a) V D D = 12.5 V, (b) V D D = 15 V, and (c) V D D = 20 V, with σ / μ = ± 0.1 for V T H , N ( P ) and for μ N ( P ) . (d) f O S C , v a r as function of the temperature in terms of V D D resulting from (ac).
Sensors 23 09653 g007
Figure 8. Results of a 1000 points process Monte Carlo analysis for (a) V D D = 12.5 V, (b) V D D = 15 V, and (c) V D D = 20 V, with σ / μ = ± 0.1 for only V T H , N ( P ) , whereas μ N ( P ) is constant at nominal value of Table 1. (d) f O S C , v a r as function of the temperature in terms of V D D resulting from (ac).
Figure 8. Results of a 1000 points process Monte Carlo analysis for (a) V D D = 12.5 V, (b) V D D = 15 V, and (c) V D D = 20 V, with σ / μ = ± 0.1 for only V T H , N ( P ) , whereas μ N ( P ) is constant at nominal value of Table 1. (d) f O S C , v a r as function of the temperature in terms of V D D resulting from (ac).
Sensors 23 09653 g008
Figure 9. Results of a 1000 points process Monte Carlo analysis for (a) V D D = 12.5 V, (b) V D D = 15 V, and (c) V D D = 20 V, with σ / μ = ± 0.1 for only μ N ( P ) , whereas V T H , N ( P ) is constant at nominal value of Table 1. (d) f O S C , v a r as function of the temperature in terms of V D D resulting from (ac).
Figure 9. Results of a 1000 points process Monte Carlo analysis for (a) V D D = 12.5 V, (b) V D D = 15 V, and (c) V D D = 20 V, with σ / μ = ± 0.1 for only μ N ( P ) , whereas V T H , N ( P ) is constant at nominal value of Table 1. (d) f O S C , v a r as function of the temperature in terms of V D D resulting from (ac).
Sensors 23 09653 g009
Figure 10. Results of a 1000 points process Monte Carlo analysis for (a) σ / μ = ± 0.15 , (b) σ / μ = ± 0.2 for μ N ( P ) and V T H , N ( P ) at V D D = 20 V. (c) f O S C , v a r as function of the temperature for different values of σ of μ N ( P ) and V T H , N ( P ) at V D D = 20 V.
Figure 10. Results of a 1000 points process Monte Carlo analysis for (a) σ / μ = ± 0.15 , (b) σ / μ = ± 0.2 for μ N ( P ) and V T H , N ( P ) at V D D = 20 V. (c) f O S C , v a r as function of the temperature for different values of σ of μ N ( P ) and V T H , N ( P ) at V D D = 20 V.
Sensors 23 09653 g010
Figure 11. Results of a 1000 points process Monte Carlo analysis for (a) σ / μ = ± 0.15 and (b) σ / μ = ± 0.2 , for only μ N ( P ) , whereas V T H , N ( P ) is constant at nominal value of Table 1 at V D D = 20 V. (c) f O S C , v a r as function of the temperature for different values of σ of μ N ( P ) with V T H , N ( P ) at nominal value and at V D D = 20 V.
Figure 11. Results of a 1000 points process Monte Carlo analysis for (a) σ / μ = ± 0.15 and (b) σ / μ = ± 0.2 , for only μ N ( P ) , whereas V T H , N ( P ) is constant at nominal value of Table 1 at V D D = 20 V. (c) f O S C , v a r as function of the temperature for different values of σ of μ N ( P ) with V T H , N ( P ) at nominal value and at V D D = 20 V.
Sensors 23 09653 g011
Figure 12. Results of a 1000 points process Monte Carlo analysis for (a) σ / μ = ± 0.15 and (b) σ / μ = ± 0.2 for only V T H , N ( P ) , whereas μ N ( P ) is constant at nominal value of Table 1 at V D D = 20 V. (c) f O S C , v a r as function of the temperature for different values of σ of V T H , N ( P ) with μ N ( P ) at nominal value and at V D D = 20 V.
Figure 12. Results of a 1000 points process Monte Carlo analysis for (a) σ / μ = ± 0.15 and (b) σ / μ = ± 0.2 for only V T H , N ( P ) , whereas μ N ( P ) is constant at nominal value of Table 1 at V D D = 20 V. (c) f O S C , v a r as function of the temperature for different values of σ of V T H , N ( P ) with μ N ( P ) at nominal value and at V D D = 20 V.
Sensors 23 09653 g012
Figure 13. Results of a 1000 points process Monte Carlo analysis for (a) σ / μ = ± 0.15 and (b) σ / μ = ± 0.2 for μ N ( P ) and V T H , N ( P ) at V D D = 15 V. (c) f O S C , v a r as function of the temperature for different values of σ of μ N ( P ) and V T H , N ( P ) at V D D = 15 V.
Figure 13. Results of a 1000 points process Monte Carlo analysis for (a) σ / μ = ± 0.15 and (b) σ / μ = ± 0.2 for μ N ( P ) and V T H , N ( P ) at V D D = 15 V. (c) f O S C , v a r as function of the temperature for different values of σ of μ N ( P ) and V T H , N ( P ) at V D D = 15 V.
Sensors 23 09653 g013
Figure 14. Results of a 1000 points process Monte Carlo analysis for (a) σ / μ = ± 0.15 and (b) σ / μ = ± 0.2 for only μ N ( P ) , whereas V T H , N ( P ) is constant at nominal value of Table 1 at V D D = 15 V. (c) f O S C , v a r as function of the temperature for different values of σ of μ N ( P ) with V T H , N ( P ) at nominal value and at V D D = 15 V.
Figure 14. Results of a 1000 points process Monte Carlo analysis for (a) σ / μ = ± 0.15 and (b) σ / μ = ± 0.2 for only μ N ( P ) , whereas V T H , N ( P ) is constant at nominal value of Table 1 at V D D = 15 V. (c) f O S C , v a r as function of the temperature for different values of σ of μ N ( P ) with V T H , N ( P ) at nominal value and at V D D = 15 V.
Sensors 23 09653 g014
Figure 15. Results of a 1000 points process Monte Carlo analysis for (a) σ / μ = ± 0.15 and (b) σ / μ = ± 0.2 for only V T H , N ( P ) , whereas μ N ( P ) is constant at nominal value of Table 1 at V D D = 15 V. (c) f O S C , v a r as function of the temperature for different values of σ of V T H , N ( P ) with μ N ( P ) at nominal value and at V D D = 15 V.
Figure 15. Results of a 1000 points process Monte Carlo analysis for (a) σ / μ = ± 0.15 and (b) σ / μ = ± 0.2 for only V T H , N ( P ) , whereas μ N ( P ) is constant at nominal value of Table 1 at V D D = 15 V. (c) f O S C , v a r as function of the temperature for different values of σ of V T H , N ( P ) with μ N ( P ) at nominal value and at V D D = 15 V.
Sensors 23 09653 g015
Figure 16. f O S C , n o r m -T curves from the results of 1000 points Monte Carlo analysis with σ / μ = ± 0.1 for V T H , N ( P ) and μ N ( P ) . The nominal temperature and bias voltages values are, respectively, (a) V D D = 12.5 V, T = 298.15 K, and 473.15 K; (b) V D D = 15 V, T = 298.15 K, and 423.15 K; and (c) V D D = 20 V, T = 298.15 K.
Figure 16. f O S C , n o r m -T curves from the results of 1000 points Monte Carlo analysis with σ / μ = ± 0.1 for V T H , N ( P ) and μ N ( P ) . The nominal temperature and bias voltages values are, respectively, (a) V D D = 12.5 V, T = 298.15 K, and 473.15 K; (b) V D D = 15 V, T = 298.15 K, and 423.15 K; and (c) V D D = 20 V, T = 298.15 K.
Sensors 23 09653 g016
Figure 17. Temperature error as function of the temperature of the f O S C , n o r m in Figure 16 with σ / μ = ± 0.1 for V T H , N ( P ) and μ N ( P ) . The error is with respect to (16) at V D D (a) 15 V and (c) 20 V, and to a third-order curve model at V D D (b) 15 V and (d) 20 V.
Figure 17. Temperature error as function of the temperature of the f O S C , n o r m in Figure 16 with σ / μ = ± 0.1 for V T H , N ( P ) and μ N ( P ) . The error is with respect to (16) at V D D (a) 15 V and (c) 20 V, and to a third-order curve model at V D D (b) 15 V and (d) 20 V.
Sensors 23 09653 g017
Figure 18. Layout of the proposed sensor.
Figure 18. Layout of the proposed sensor.
Sensors 23 09653 g018
Table 1. Circuit specifications.
Table 1. Circuit specifications.
ParameterUnitValue
V D D [V]20
V S [V]−8
f O S C [kHz]90
W N C A P L N C A P μ m μ m 270 270
C 1 [pF]8
C 2 [pF]8
L m i n [ μ m]6
V T + [V]10
V T [V]5
Table 2. Results of the analytical design of the proposed circuit using (1)–(14).
Table 2. Results of the analytical design of the proposed circuit using (1)–(14).
ParameterDeviceValueUnit
W L M N 1 36 / 6 μ m μ m
M N 2 6 / 6
M N 3 6 / 6
M P 1 6 / 6
M P 2 6 / 6
M P 3 12 / 6
M N 4 6 / 6
M P 4 6 / 6
M N 5 6 / 6
M P 5 12 / 6
M N 6 12 / 6
M P 6 42 / 6
M N 7 18 / 6
R S 94[ K Ω ]
Table 3. Results of the design after the tuning using numerical simulations. The changed parameters compared to Table 2 are highlighted in red-bold.
Table 3. Results of the design after the tuning using numerical simulations. The changed parameters compared to Table 2 are highlighted in red-bold.
ParameterDeviceValueUnit
W L M N 1 6/6 μ m μ m
M N 2 6 / 6
M N 3 6 / 6
M P 1 6 / 6
M P 2 6 / 6
M P 3 12 / 6
M N 4 6 / 6
M P 4 6 / 6
M N 5 12/6
M P 5 32/6
M N 6 24/6
M P 6 64/6
M N 7 24/6
R S 60[ K Ω ]
Table 4. Maximum f O S C , v a r for σ / μ = ± 0.1 and at various V D D .
Table 4. Maximum f O S C , v a r for σ / μ = ± 0.1 and at various V D D .
V DD f OSC , var
σ / μ | V TH , N ( P ) σ / μ | μ N ( P ) σ / μ | V TH , N ( P ) μ N ( P )
12.5  V 22.64 % 5.86 % 23.33 %
15 V 11.76 % 5.72 % 12.41 %
20 V 7.24 % 5.98 % 8.41 %
Table 5. f O S C , v a r maximum value for different σ / μ and at V D D = 20 V.
Table 5. f O S C , v a r maximum value for different σ / μ and at V D D = 20 V.
V TH , N ( P ) μ N ( P ) V TH , N ( P ) μ N ( P )
σ / μ ± 10 % ± 15 % ± 20 % ± 10 % ± 15 % ± 20 % ± 10 % ± 15 % ± 20 %
V D D = 15 V 11.82 % 18.04 % 26.64 % 5.72 % 8.73 % 11.99 % 12.42 % 19.86 % 27.11 %
V D D = 20 V 7.24 % 10.92 % 14.66 % 5.98 % 9.03 % 12.36 % 8.41 % 12.77 % 17.4 %
Table 6. Comparison between our proposal and the state-of-the-art.
Table 6. Comparison between our proposal and the state-of-the-art.
RangeErrorBias Volt. V DD PowerAreaTech.
[K][K][V][mW][mm2]
This work (Num.Sim.)298/573−5.8/8.8150.890.1634H-SiC 2  μ m-CMOS
[16]273/373−4/41.21.20.12Si 130 nm-CMOS
[28]233/378± 0.5 ( 3 σ ) 52.52.3Si 0.7 μ m-CMOS
[29]273/373−3/310.1540.004Si 65 nm-CMOS
[30]273/383−1.5/1.510.50.008Si 65 nm-CMOS
[31]272/385−2.5/1.2 ( V D D  = 0.6 V)
−1.4/1.3 ( V D D  = 1.2 V)
0.6 V/1.2 V0.0560.001FD-SOI 28 nm-CMOS
[32]273/373±1.95 (3 σ ) 1.65V0.10.001SOI 32 nm-CMOS
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Rinaldi, N.; Liguori, R.; May, A.; Rossi, C.; Rommel, M.; Rubino, A.; Licciardo, G.D.; Di Benedetto, L. A 4H-SiC CMOS Oscillator-Based Temperature Sensor Operating from 298 K up to 573 K. Sensors 2023, 23, 9653. https://doi.org/10.3390/s23249653

AMA Style

Rinaldi N, Liguori R, May A, Rossi C, Rommel M, Rubino A, Licciardo GD, Di Benedetto L. A 4H-SiC CMOS Oscillator-Based Temperature Sensor Operating from 298 K up to 573 K. Sensors. 2023; 23(24):9653. https://doi.org/10.3390/s23249653

Chicago/Turabian Style

Rinaldi, Nicola, Rosalba Liguori, Alexander May, Chiara Rossi, Mathias Rommel, Alfredo Rubino, Gian Domenico Licciardo, and Luigi Di Benedetto. 2023. "A 4H-SiC CMOS Oscillator-Based Temperature Sensor Operating from 298 K up to 573 K" Sensors 23, no. 24: 9653. https://doi.org/10.3390/s23249653

APA Style

Rinaldi, N., Liguori, R., May, A., Rossi, C., Rommel, M., Rubino, A., Licciardo, G. D., & Di Benedetto, L. (2023). A 4H-SiC CMOS Oscillator-Based Temperature Sensor Operating from 298 K up to 573 K. Sensors, 23(24), 9653. https://doi.org/10.3390/s23249653

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