Next Article in Journal
The Efficiency Comparison of Hydro Turbines for Micro Power Plant from Free Vortex
Next Article in Special Issue
Gate Current in p-GaN Gate HEMTs as a Channel Temperature Sensitive Parameter: A Comparative Study between Schottky- and Ohmic-Gate GaN HEMTs
Previous Article in Journal
Optimization and Field Test of a Chelating Acid System for Scaled Gas Wells in the Hechuan Gas Field
Previous Article in Special Issue
Cascaded Smart Gate Drivers for Modular Multilevel Converters Control: A Decentralized Voltage Balancing Algorithm
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

VDS and VGS Depolarization Effect on SiC MOSFET Short-Circuit Withstand Capability Considering Partial Safe Failure-Mode

1
IRT Saint-Exupery, CS34436, 3 Rue Tarfaya, 31400 Toulouse, France
2
Laplace, University of Toulouse, CNRS, INPT, UPS, 2 Rue Camichel BP7122, 31071 Toulouse, France
*
Author to whom correspondence should be addressed.
Energies 2021, 14(23), 7960; https://doi.org/10.3390/en14237960
Submission received: 20 September 2021 / Revised: 9 November 2021 / Accepted: 23 November 2021 / Published: 29 November 2021
(This article belongs to the Special Issue Safety Design and Management of Power Devices including Gate-Drivers)

Abstract

:
This paper presents a detailed analysis of 1200 V Silicon Carbide (SiC) power MOSFET exhibiting different short-circuit failure mechanisms and improvement in reliability by VDS and VGS depolarization. The device robustness has undergone an incremental pulse under different density decreasing; either drain-source voltage or gate-driver voltage. Unlike silicon device, the SiC MOSFET failure mechanism firstly displays specific gradual gate-cracks mechanism and progressive gate-damage accumulations greater than 4 µs/9 J·cm−2. Secondly, a classical drain-source thermal runaway appears, as for silicon devices, in a time greater than 9 µs. Correlations with short-circuit energy measurements and temperature simulations are investigated. It is shown that the first mechanism is an incremental soft gate-failure-mode which can be easily used to detect and protect the device by a direct feedback on the gate-driver. Furthermore, it is highlighted that this new mechanism can be sufficiently consolidated to avoid the second drain-source mechanism which is a hard-failure-mode. For this purpose, it is proposed to sufficiently depolarize the on-state gate-drive voltage to reduce the chip heating-rate and thus to decouple the failure modes. The device is much more robust with a short-circuit withstand time higher than 10 µs, as in silicon, no risk of thermal runaway and with an acceptable penalty on RDS-ON.

1. Introduction

The 4H-Silicon Carbide (4H-SiC) material has much better physical characteristics than silicon (Si). With a threefold higher band-gap energy, a tenfold higher breakdown field, and an extremely low intrinsic free carrier concentration, 4H-SiC is a leading material for high-voltage power devices. Due to the difficulties of deep and wide P+ doping implantation, bipolar devices such as BJT (Bipolar Junction Transistor) and IGBT (Insulated Gate Bipolar Transistor) are not as easy to process as for silicon. A unipolar SiC device such as MOSFET (Metal Oxide Silicon Field Effect Transistor) is preferred. However, non-native oxide based SiO2-SiC structures still suffer from imperfections at the SiC interface leading to reduced mobility in the channel with gate-threshold voltage drift. However, SiC MOSFETs offer several advantages compared to silicon device (low on-state and switching losses, higher switching frequencies, and higher temperature stability), which makes them an alternative solution of IGBT or silicon Super-Junction MOSFETs, up to 30% to 50% losses reduction can be achieved depending on switching frequency and voltage rating [1]. SiC MOSFETs are taking over different markets, where they are found in different applications, especially on-board chargers or power trains in automotive for medium power and railways for high power. However, in extreme operation, short-circuit (SC) robustness is a key application requirement for semiconductor power devices. The short-circuit (SC) withstanding time (TSCWT) of standard and commercial SiC MOSFETs is much lower than that of silicon devices [2]: from 10 µs, which is now practically the IGBT standard, to a few µs for recent commercial MOSFET device. This difference is due to the continuous reduction of the chip area and thus the increase of the current density in the extreme short-circuit operation. Indeed, manufacturers optimize SiC MOSFET from the point of view of reducing switching times and energies by reducing input–cross–output structural capacitors (Ciess, Crss, Coss) of the chip, with the same apparent drain-source on-state resistance (RDSON). This implies faster gate-driver protection, which must also remain robust against switching oscillations at turn-on. In order to compensate the low short-circuit robustness of SiC MOSFET, specific designs have been analyzed to find the best trade-off between TSCWT and RDSON·cm2 (normalized on-state drain-source resistance). A first approach consists in optimizing the gate geometry while remaining in usual planar structure. In [3], authors show that a thinning of the gate from 55 nm under VGS = 20 V to 27 nm under VGS = 10 V allows to increase TSCWT by 47% at the cost of only a 12% increase in TSCWT. In contrast, in [4], the authors show that a thickening of the gate from 50 nm to 85 nm allows the TSCWT to be increased by 50% but at the cost of a more penalizing increase of 30% in RDSON·cm2. Source region (N+) around the gate is also an interesting parameter. In [5], authors show that adjusting the doping downwards is the most efficient parameter to both increase the TSCWT from 4 µs to 10 µs with a +20% penalty on RDSON·cm2. Another approach is to redesign the gate shape. In [6], authors show that an adapted trench gate structure allows the maximum field amplitude and temperature level to be reduced to the same short circuit time as a planar gate. The time for which failure will be visible on such a device will therefore be higher, resulting in a more robust component. In [7], trench structure can also be completed by a so-called super-junction structure which amplifies the above properties with an increase of 10% to 50% of the TSCWT at the cost of a very acceptable increase of the RDSON·cm2. All these approaches are interesting in principle, to reach the best trade-off between TSCWT and RDSON·cm2, but they require a more expensive design and manufacturing effort than a standard design. So, in this paper, authors propose solutions to improve robustness based on standard and commercial components. In addition, a particular effort is made to finely analyze the device’s degradation modes and specifically those that preserve the application’s operating safety by avoiding a drain-source short-circuited state.
The SC performances and failure mechanisms of standard and commercial SiC MOSFET and Si IGBT have been discussed in different studies [8,9,10]. The SiC MOSFET mainly presents two failure mechanisms, under SC operations, both known as drain-source fail-to-open-circuit (FTO, or soft failure mode) and drain-source fail-to-short-circuit (FTS, or hard failure mode) [9,11,12]. Clearly, these signatures depend on the power density applied and then the temperature rise within the device. The FTS mode is derived from thermal runaway by the parasitic BJT p-well region triggering the loss of the turn-off and the failure of the chip. On the other hand, the specific FTO mode; which is by far the preferred mode; results at VGS FAILED < VGSTH leading to transistor shutting off, an intrinsic fault self-containment capability with enabling the system safe operation without complex or costly auxiliary components. This failure mode can be due mainly to the progressive ILD (Internal Layer Dielectric) cracks, fusion, and diffusion of the Al top metal, leading to the partial shorting of the gate and source terminals. This is related to extreme thermo-mechanical stresses and heat in SC between the polysilicon gate and the Al source.
In this context, a direct approach to make commercial SiC MOSFET more robust and safer is to reduce the power density applied to the chip, either by the drain voltage depolarization or by the gate voltage depolarization. In [13], authors present an analysis of the main parameters VDS bias, VGS bias, and RG (external gate resistance) on short-circuit electrical stresses with an estimation of the transient chip temperature based on a TCAD model. Parameter VDS bias is the most important one to reduce the temperature, but the authors do not present any estimate or measurement of the corresponding TSCWT and no consequence on the failure modes. In [14], authors propose a fairly comprehensive exploration of the reduction of VGS and VDS from 800 V to 400 V in planar and trench structures. At 800 V VDS bias, authors show that the FTS mode cannot be avoided and that a VGS depolarization from 15 V bias to 12 V bias increases TSCWT by 40%. At 600 V VDS bias, TSCWT increases by 53% and the FTS mode can be neutralized to leave the FTO mode visible, provided that the channel turn-off is reinforced with a VGS bias of −10 V instead of −5 V. Finally, authors also show that at 400 V VDS bias, both structures exhibit a FTO mode with a 250% higher TSCW. However, although interesting, this article does not detail the degradation mechanisms of the gate.
This paper aims to present different failure mechanisms of a recent, second generation, ST-Microelectronics™ device 1200 V, 75 mΩ, 33 A @Tcase = 25 °C planar-gate type device, in an automotive-grade 3-lead TO247 case [15], under incremental-pulse short-circuit behavior. The DC testing voltage was at nominal and reduced 600 V operation instead of standard 800 V used in extreme short-circuit assessment. The impact of gate-source decreased from 18 V to 10 V is also analyzed on the TSCWT and RDSON trade-off.
In contrast to the state-of-the-art literature, we will pay particular attention to the mechanisms and failure modes specific to the SiC MOSFET under reduced or depolarized voltage conditions, in particular the properties that preserve the electrical safety of the device and the absence of drain-source short-circuit.
The paper structure is divided by two main approaches: the drain-source voltage VDS depolarization and gate-driver voltage VDrv depolarization. Each section will discuss both the drain-source voltage and the gate-driver voltage depolarization. The first section is focused on the short-circuit electrical behavior under incremental pulse width, whereas the second section is focused on the electrical effect after the first failure. The paper points out the effect of progressive gate-damage accumulation, by discussing observations over thermal and energy dissipation. It will show that this specific SiC MOSFET mechanism can be easily detected to self-protect the device with a feedback on the gate-driver. Moreover, the authors studied the depolarize on-state gate-drive voltage sufficiently to reduce the chip heating-rate and thus increasing the robustness of the short-circuit withstand time above 10 µs, as in silicon, without the risk of thermal runaway [16].

2. Short-Circuit Stress Analysis

2.1. Applied Stress and Power Test Bench

In order to electrically identify and characterize the device short-circuit withstand time, short-circuit energy, and failure mechanism, the device was placed under extreme short-circuit mode; which is known as Hard Switch Fault HSF [16]. The DUT (Device Under Test) was placed on the low side of a phase-leg power test bench presented in Figure 1. On the high side of the power test bench, a connection was made between the drain and source electrode of the upper device of the leg by using a copper line in order to simulate the HSF short-circuit behavior. The power test bench presented in Figure 1 was supported by a fast high current electronic breaker and safety fuses, to limit the final damage of the device from a possible physical explosion. The power test bench also consisted of an adapted gate driver with on-state and off-state adjustable voltage, controlled by an FPGA (field programmable gate arrays) and a touch screen, in order to ease the short-circuit test manipulation.
One should note, there is a straightforward relationship between the power density in short-circuit operation and the maximum temperature (1) [17]. This relationship can help to easily identify different features about the DUT. This relation assumes a homogeneous surface density on the top of the chip which is well suited to the case of the SiC MOSFET. Indeed, its epitaxy is thin and highly doped compared to silicon device, with a strong gradient of narrow field, very close to the chip surface. Moreover, ref. [9] presents that both failure modes depend on the temperature increase and the growth rate within the device. In order to avoid a fast critical FTS failure, the temperature increase within the device should be controlled or limited. One can use Equation (1) to decrease the junction temperature by decreasing the power density. The power density can be decreased either by decreasing the drain voltage VDS, or by decreasing the density of saturation current by depolarizing the gate-drive voltage VDrv. This paper presents both methods. The first approach was decreasing the drain voltage by using 600 V DC bus voltage instead of the nominal 800 V used in application. The second approach was to further decrease the power density by VGS depolarization.
Δ T jmax ( t ) = 2 × V DS × J sat ρ c p λ π t
where r is the density of SiC, cp is the specific heat of SiC and l is the thermal conductivity of SiC. VDS is the drain-source voltage. JSat is the saturation channel current. t is the transient time, proportional to short-circuit time in our case.

2.2. DC Volatage Depolarization Method

2.2.1. Failure Mode Sequences for a Given Device

Sequence 0:
The applied test is a successive incremental short-circuit pulse with an initial pulse of 0.5 µs and a step of 0.25 µs, until the device reaches its physical destruction. Different DUTs have been tested. For this study, seven components have been tested. The test is under a nominal VDrv (−5 V/18 V) driving with a gate resistor RG = 47 Ω and a reduced VBus = 600 V. The DUT was designed for 1.2 kV compared to VStandard = 2/3 × 1200 kV = 800 V, which is the specified range of the maximum transient bus voltage typically in reverse operation of automotive power train. The RG gate resistor is deliberately chosen to be high in order to avoid turn-off power surges and voltage spikes instead of using a dedicated soft turn-off circuit. Moreover, in this study this gate resistor is used to ease local measurement of dynamic, temporal gate leakage IGS(t). Most of the figures are oscilloscope views illustrating VDrv, VGS, VDS, ID, and IG. The gate current IG is given by the difference between VDrv and VGS divided by the gate resistor.
Figure 2a presents preliminary oscilloscope view at tPulse = 1 µs of the device number #3. This figure is introduced in order to get familiarized with a single pulse waveform. Figure 2b illustrates the incremental short-circuit test. The gate current illustrated is a cropped view of the whole signal, the positive and negative current peaks are not illustrated; only to enhance and amplify the current leakage signal if illustrated. Each component is led to its final destruction. During the test run, different failure modes and degradations have been noticed.
Sequence 1:
After several incremental pulses, the first degradation that has been noticed is that VGS starts to slightly decrease (point I. in Figure 3); (green waveform, Figure 3 at 3.5 µs for DUT #3). This behavior is conditioned by electric field stress in the oxide combined with high temperature in the gate-oxide caused by extreme power density applied, in the channel and the JFET regions [18]. This behavior is known by the gate current-leakage, (point II. In Figure 3). The physical origin of the gate leakage-current is still under study by the authors, but it is strongly assumed to be related to indirect tunneling effect at medium electric field and high temperature (Schottky emission). In Figure 3, the VGS drop is directly translated to an ohmic drop which is a high gate leakage current through the external gate resistor 47 Ω, as illustrated in the red waveform (IG). This non-permanent and transient gate leakage current is a great marker in order to detect the presence of a short-circuit, as also suggested in [19,20].
Sequence 2:
Afterwards, the pulse duration increases, therefore an increase in the temperature, eventually leads to a significant drain current leakage increase (point III. In Figure 3); known as the current tail [21]. This drain leakage is visible in Figure 3 at 4.25 µs.
As the pulse length gets longer, VGS drop becomes significant, which helps reduce the dissipated power and brings the device to a less thermal stable point. This can result in a delayed thermal runaway. Moreover, the high temperature during SC events exerts thermo-mechanical stresses to the oxide due to a significant misalignment between thermal expansion coefficients αAL and αSiO2 (ratio ≈ 40) [22]. In this case, with lower VDS, the DUT goes into a partial and delayed gate damage as shown in Figure 4, which is an irreversible phenomenon, (point IV. & V. in Figure 4).
The delayed crack [22] as presented in Figure 4 is probably due to a mechanical time constant plus the penetration duration or the diffusion of the metal in the cracks.
As presented in Figure 4 when the pulse duration reaches 5.25 µs, the gate-source voltage shows a depolarization from −5 V to −4 V (point V. in Figure 4), after 5 µs delay, which is translated by an increase in the gate current leakage. Furthermore, The DUT presents an appearance of a first failure mode, known as a soft gate damage.
The gate failure seems to be a reaction of a transient electro-thermal-mechanical process. In order to study this first failure mode in depth, a cell-level 1D thermo-metallurgical simulation is presented in the fourth section. This failure behavior can be the outcome of the increase of the temperature resulting into aluminum melting and stress on SiO2 which exceeds it mechanical stress [23].
Sequence 3:
After the appearance of the first partial gate damage, the next pulse; which is 250 ns longer than the previous one (5.25 µs); clearly indicates the confirmation of a permanent partial gate damage (point V. in Figure 5). The waveforms show that VGS is softly depolarized from 18 V to 15 V and from −5 V to −4 V. Moreover, IDMax is significantly reduced (point VII. in Figure 5), as well as the power density. These reductions allow the device to support few more short-circuit pulses, and allowing the temperature and the power density within the device to be furthermore increased. On the other hand, the gate current leakage is increased, and will remain very high for future pulses.
However, such a high level of permanent gate-current leakage or gate-voltage collapse can be easily detected by the driver. This property can be used to locally self-turn-off the device to prevent it from getting too close to its critical failure mode, FTS [19,24].
Sequence 4:
Afterwards, the device was exposed to longer and longer pulses duration, Figure 5. In the process the device showed more and more drain leakage (point VII. in Figure 5), again with an increase on the temperature and the power density, which brought the device to a second gate damage (point VI. in Figure 5). This resulted in a reduction of the drain leakage tail (point VIII. in Figure 5), and further depolarization of VGS (point VI. in Figure 5). This pushed the device to a further extreme operational point, and continued until the device reached a final critical failure, FTS @7.5 µs (point IX. In Figure 5).
These different post cracks were a clear marker of a permeant physical degradation of the device. Different studies show the presence of the gate damage failure mode [25,26], and the successive degradations [27]. However, some authors have shown that it is also possible to partially restore the gate-integrity by a dedicated original gate-bias procedure, in off-line operation, provided that the gate is not totally cracked [28].

2.2.2. DC Voltage Depolarization Method under Several Devices

The sequences in the section above and the study of different devices helped to illustrate the timeline presented in Figure 6. The timeline presents, the appearance of gate and drain leakage as well as the gate damage failure and the final failure FTS respectively.
The incremental pulse duration test is applied on different devices separately. After the gate-damage, each device shows its unique crack behavior and depolarization amplitude; regarding the decrease of the drain current ID and the amplitude depolarization of VGS. Moreover, the number of post-cracks differs from a device to another, where it can be one crack or up to four cracks, as illustrated in Figure 5 and Figure 7.
For a large scope of study, different tests have been applied on different devices:
  • Different SC initial pulse duration:
Oscilloscope waveforms of a different device are presented in Figure 7. The applied test on this device was only different by its initial pulse. This initial pulse duration time was 3 µs instead of 0.5 µs. The purpose of this applied test was to see the effect of the gate-damage time by reducing the early stress on the device. One should note that the degradation of the device was not noticed when a pulse longer than the average gate-damage pulse (5 µs) was applied to a non-stressed device. Figure 7 presents the behavior of the DUT #4 after gate damage which occurs at 5 µs. Similar behavior has been noticed for those devices with a different initial pulse. As the one can see; after the presence of gate damage (point V.); the presence of point VI., VII., VIII., and IX. are visible with different amplitude.
  • Long SC pulse:
In the case where the pulse width equaled 110% of the gate-damage pulse or longer, the device went directly to a critical failure mode, FTS. The device had no time to consolidate the gate-damage behavior, therefore there were no visible cracks and no SC cycling capability. It could be due to the absence of the non-visible cracks, or fatigue, and no partial degradation before reaching certain energy amount.
Different key parameters of several devices are presented in Table 1. The appearance of the gate damage was reproducible around 4.75 µs and 5.25 µs, contrary to the FTS failure mode. Moreover, the case temperature has no clear influence on the gate-damage failure time.

2.3. Gate Voltage Depolarization Method

This section is dedicated to the second depolarization method. This method is based on the direct decrease of VDrv voltage. Figure 8 shows the results of short circuit tests carried out on the DUT #8 under (VDS = 600 V, VDrv = 10 V) using a direct depolarization of VGS; bringing the maximum short-circuit current from 250 A; observed for DUT #4; to less than 115 A; for DUT #8 (i.e., 50% reduction with a more softened form).
An increase in the Tsc with a theoretical ratio close to four (×4) is to be expected, according to Relation (1) and considering an invariant failure temperature threshold between the two bias modes of the gate voltage, (quadratic effect between Tsc and power density).
Indeed, it can be observed that direct depolarization of VGS method showed promising results exhibiting a TSC equal to 16 µs; compared to 5 µs for DUT #4. This confirmed the previous analysis while being slightly less advantageous. However, the method allowed exceeding 10μs (standard for IGBT device), with the appearance of FTO mode under a VDS close to the nominal value.
Figure 8. presents the first failure mode (gate damage or partial FTO), which occurred around 11 µs, and the full FTO failure, which appeared at 16 µs (point X. in Figure 8). The impact of a direct VGS depolarization on RDS-ON is illustrated in Section 3.2.

3. Electrical and Structural Analysis

3.1. Energy Study under Short-Circuit

The gate damage and the degradation that followed on the device can be seen and interpreted from the energy dissipated in the device. The degradations were translated by several post-cracks. The energy dissipated (EDis) inside the DUT was calculated from the VDS and ID waveforms captured during the SC phase by integrating VDS × ID over the duration of the SC pulse TSC-Pulse (2).
E Dis = 0 t SC Pulse V Ds × T D   d t
Figure 9a,b present the short-circuit energy of the DUT #3 and DUT #4 respectively. The figures perfectly exhibit the depolarization dips of the device and the presence of the successive gate damage including the amplitude of the cracks. The DUT #3 shows at the 1st appearance of the gate damage that the dissipated energy was reduced from 556 mJ to 534 mJ allowing the device to withstand more temperature increase. The same reduction could be noticed for the cracks that followed. Similarities to this behavior are also reproducible for different DUT, as presented in Figure 9c for the DUT #6.
For all DUT tested the gate-damage occurred between 530 mJ and 550 mJ. On the other hand, the FTS occurred between 581 mJ and 606 mJ. The robustness margin was small between the soft gate-damage and the hard failure mode. This small margin was visible when 110% of the gate-damage pulse was applied on an unstressed device; exhibiting directly an FTS.
From Figure 10, one can estimate that improving the gate, i.e., by increasing the mechanical strength of the ILD using a material other than SiO2, might slightly increase the apparent withstand time from 5.25 towards 5.75 µs, emphasizing a critical failure “FTS” under a nominal VGS bias. This improvement will hide the gate-depolarization phenomena therefore avoiding the safe failure mode. In contrast to the safe damage behavior at 5.25 µs, the device still has switching capability (where VGS is de-polarized, and the device only softly damaged and not destroyed).

3.2. RDSON Analysis

The occurrence of the first crack led to a depolarization of the gate voltage, which is reduced by few volts at most. Figure 11 shows that this depolarization affects the RDS-ON increase by no more than 10%. In general cases where conduction and switching losses are equal, this penalty on the RDS-ON will only impact the total losses by 5% which is very reasonable. However, this result has to be considered on a case-by-case scenario depending on the design of the component and the shape of its ID (VDS) characteristic.
If this first crack is detected by the driver and the fault is inhibited, the same component can be reused without having to be replaced. However, the gate-driver supply must feed a permanent current overload of several tens of mA. An excess of power to be transferred to the gate which can reach +1 W in the driver power supply. This over-supply is significant and should be taken into account at the design stage. Furthermore, this over-supply can also be used as a means of detecting the occurrence of a hard crack to manage the shutdown of the converter and then the replacement of the component at the appropriate time regarding the mission profile.
The gate depolarization described in Section 2.2.1 obviously also involves an increase in RDS-ON. However, this penalty is more acceptable for the SiC MOSFET because this component has a channel resistance that decreases at intermediate temperatures (up to 300 °C) compared to the other series resistance components RJFET and RDRIFT. This effect is related to Coulomb-like mobility through scattering of ionized impurities. As the channel resistance becomes of majority value in depolarization mode, this property gives an overall resistance that is more temperature stable than a silicon MOSFET. For the DUT studied in this paper, RDS-ON (VGS) is drawn in Figure 12 and points out an increase of 23%. In the general case, conduction and switching losses are equal, this penalty on the RDS-ON will only impact the total losses by half (11.5%) which is acceptable.

4. Thermal 1D Simulation Analysis

In order to support the gate-damage analysis, thermal modelling and simulation of the chip temperature were carried out, as detailed in [29]. The model takes into account the temperature dependence of the conductivity and the specific heat of SiC and aluminum, as well as the solid–liquid phase transition of the Al top layer. Using structural computer aided design tools (Comsol™ was used in this case), a cell-level 1Dx/Ex—electric field, Jx—current density model was carried-out. The software receives as input the current waveform IDS(t) of the DUT #4 from the experimental bench at a constant voltage VDS = 600 V and VGS = 18 V/−5 V.
Figure 13 shows the temperature rise at Tpulse_G_D = 5 µs of all the main layers of the chip during the short-circuit. After only 3 µs, it is visible that the Al layer was above its melting point with a high thermal energy injected into its liquid phase. Moreover, according to [29], at the Al melting temperature, the thick oxide was probably already cracked because the mechanical stress exceeded its mechanical strength which is confirmed previously by the gate-crack electrical observation on Figure 6. The final and maximum temperature in the chip was evaluated by simulation at 1280 °C. This value remained lower than the polysilicon level fusion (1400 °C), SiO2 level fusion (1700 °C) and nitride passivation layer Si3N4 (1900 °C). Then, all these material regions should not be severely damaged up to 5 µs short-circuit operation.
In order to explain the physical origin of FTO failure mode; Figure 14 illustrates the electro-thermal simulation of DUT #8, a device that had undergone a gate depolarization method; VDrv = 10 V. The simulation shows that the temperature on the top layer of the Al reached the start point of the metallurgical phase transition (solid–liquid) around 8 µs, i.e., more than twice the delay time of DUT #4. The temperature of the top Al layer starts rising again and follows the thermal dynamics of the junction after about 8 µs of stable phase.
The final and maximum temperature in the chip is simulated at 1242 °C, which is only 3% lower than in case DUT #4, while the short-circuit pulse times are in a ratio of 3.2 (5 µs vs. 16 µs). It seems therefore that this threshold temperature could be a salient indicator of the onset of the gate degradation mechanism.
After the end of the pulse, it is clearly seen that the temperature at the bottom Al layer became lower than the temperature at the top Al layer: the metal layer restituted thermal energy in a liquid state. When the temperature reached the phase transition zone, the heat sink reactivated and considerably reduced the thermal dynamics of the bottom Al layer and the junction. Afterwards, the accumulated energy was released during the following microsecond, leading to the solidification of the whole Al layer. Through a thermo-mechanical stress point of view, it is obvious that the estimated temperature of the metal top was above 600 K, which led to the mechanical rupture of the ILD—SiO2 to 1.4 GPa, by compression to the Al layer on the top of SiO2 [29].

5. Discussion

In order to fully understand the effect of both depolarization methods (VGS and VDS) a comparative study is presented in Table 2, including results of the same rating device under nominal conditions at VDS = 800 V [15]. As one can see, under nominal conditions the device withstands around TSCWT = 1.62 µs; with a critical FTS mode. On the other hand, using VDS depolarization method, the device withstands around TSCWT = 4.9 µs, with a safe failure mode (gate damage). Finally, using VGS depolarization method, the device withstands around TSCWT = 9.6 µs, which is close to what the IGBT standard TSCWT (10 µs), with the same safe failure mode and only a 23% RDSON increase. Overall, we can note that the damage or failure temperature was decreased (1) from one method to another; (respectively: 1320 °C, 1207 °C, 1056 °C); which allowed the device to withstand more the SC; (respectively: 1.62 µs, 4.9 µs, 9.6 µs). It is also remarkable to note that the damage mode corresponding to the strong degradation of the gate (2) is well associated with comparable temperature values for VDS and VGS depolarization (respectively: 1207 °C and 1241 °C (2)). The gate leakage on the other hand is not visible under VGS depolarization, ((3), Table 2).
The experimental observations on the gate voltage in short-circuit operation and the thermal simulations make it possible to establish relevant correlations: A threshold temperature between 1200 °C and 1300 °C seems to correspond to a limit of use of the component. This is an interesting result, although it is to be expected that the duration of the pulse (total and after the melt time) will have an influence on the level of gate-damage its repetitive short-circuit cycles capability, but this last point is a research prospect for single pulse test and repetitive pulses. However, even if this threshold is exceeded, the gate-damage mode remains soft and progressive: it can be easily detected and managed by monitoring the gate voltage at the buffer output or indirectly by oversupply. The depolarization of the gate voltage from 18 V to 10 V allows the short-circuit withstand time to reach the usual 10 µs which is a standard for IGBT technology which is still the most common power device in medium and high power. This property can be inexpensive on losses, depending on the case, it also allows considerably the reduction of the level of electric field stress on the gate oxide and improves reliability on long term operation.

6. Conclusions

Unlike silicon power MOSFETs, SiC devices showed different failure mechanism at reduced drain—source voltage 600 V compared to nominal 800 V for a 1.2 kV rating with nominal gate drive voltage. At 600 V the power density was decreased showing a succession of gate degradations as a non-catastrophic failure-mode. For the 75 mOhms studied device the gate damage occurs between 530 mJ and 550 mJ around 5 µs which is perfectly reproducible. The gate damage produces a soft or hard depolarization on VGS which keeps the thermal runaway further away. After the first crack, the device was operational with a small penalty on RDS-ON. At the end, the device could withstand further longer short-circuit with successive cracks and gate depolarization which further reducesd the power density to be dissipated. This SiC device also showed further robustness on short-circuit withstand time using a gate-drive voltage depolarization. This method also showed a first failure mode gate damage, and a final failure as a complete fail-to-open mode giving VGS < VGSTH.
These results complement and clarify many previous results obtained by the authors or their collaborators in SiC planar gate-MOSFET technology. This technology is the most widely used in SiC. The trench-gate structures are not considered in this paper because they are more complex to realize with more singular properties and with a more sensitive variation of the threshold gate voltage according to the gate polarization. It would require a more specific study in line with the perspective works of the authors.

Author Contributions

Conceptualization, Y.B., F.R. and W.J.; methodology, Y.B., F.R. and W.J.; software, Y.B.; validation, Y.B., F.R. and W.J.; supervision, F.R. and J.-M.R. All authors have read and agreed to the published version of the manuscript.

Funding

This work was made possible in large part by direct or indirect funding from the French National Research Agency (ANR). A large part of the results was obtained within the framework of the internal SiCRET “Silicon Carbide Reliability Evaluation for Transport” project Ref: CDP-E-001-068-V0 led by the Saint-Exupéry Institute of Technological Research in Toulouse (IRT SE).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Acknowledgments

We would like to thank Safran Tech™ division from Safran™ company and the Power Electronics team, in particular its leader Stéphane Azzopardi, for their scientific and financial support in the development of this topic and of part of the results obtained (part. Section 2.2.2).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Baliga, B.J. Silicon Carbide Power Devices; World Scientific Publishing Co. Pte. Ltd.: Singapore, 2005. [Google Scholar]
  2. Wang, Z.; Shi, X.; Xue, Y.; Tolbert, L.; Wang, F.; Blalock, B.J. Design and Performance Evaluation of Overcurrent Protection Schemes for Silicon Carbide (SiC) Power MOSFETs. IEEE Trans. Ind. Electron. 2014, 61, 5570–5581. [Google Scholar] [CrossRef]
  3. Agarwal, A.; Kanale, A.; Han, K.; Baliga, B.J. Switching and Short-Circuit Performance of 27 nm Gate Oxide, 650 V SiC Planar-Gate MOSFETs with 10 to 15 V Gate Drive Voltage. In Proceedings of the 2020 32nd International Symposium on Power Semiconductor Devices and ICs (ISPSD), Vienna, Austria, 17–21 May 2020; pp. 250–253. [Google Scholar] [CrossRef]
  4. An, J.; Namai, M.; Yano, H.; Iwamuro, N.; Kobayashi, Y.; Harada, S. Methodology for enhanced short-circuit capability of SiC MOSFETs. In Proceedings of the 2018 IEEE 30th International Symposium on Power Semiconductor Devices and ICs (ISPSD), Chicago, IL, USA, 13–17 May 2018; pp. 391–394. [Google Scholar] [CrossRef]
  5. Bolotnikov, A.; Losee, P.A.; Ghandi, R.; Halverson, A.; Stevanovic, L. Optimization of 1700 V SiC MOSFET for Short Circuit Ruggedness. Mater. Sci. Forum 2019, 963, 801–804. [Google Scholar] [CrossRef]
  6. Yang, T.; Wang, Y.; Yue, R. SiC Trench MOSFET With Reduced Switching Loss and Increased Short-Circuit Capability. IEEE Trans. Electron Devices 2020, 67, 3685–3690. [Google Scholar] [CrossRef]
  7. Okada, M.; Kyogoku, S.; Kumazawa, T.; Saito, J.; Morimoto, T.; Takei, M.; Harada, S. Superior Short-Circuit Performance of SiC Superjunction MOSFET. In Proceedings of the 2020 32nd International Symposium on Power Semiconductor Devices and ICs (ISPSD), Vienna, Austria, 17–21 May 2020; pp. 70–73. [Google Scholar] [CrossRef]
  8. Boige, F.; Richardeau, F.; Lefebvre, S.; Blaquière, J.-M.; Guibaud, G.; Bourennane, A. Ensure an original and safe “fail-to-open” mode in planar and trench power SiC MOSFET devices in extreme short-circuit operation. Microelectron. Reliab. 2018, 88-90, 598–603. [Google Scholar] [CrossRef] [Green Version]
  9. Romano, G.; Fayyaz, A.; Riccio, M.; Maresca, L.; Breglio, G.; Castellazzi, A.; Irace, A. A Comprehensive Study of Short-Circuit Ruggedness of Silicon Carbide Power MOSFETs. IEEE J. Emerg. Sel. Top. Power Electron. 2016, 4, 978–987. [Google Scholar] [CrossRef]
  10. Chen, C.; Labrousse, D.; Lefebvre, S.; Petit, M.; Buttay, C.; Morel, H. Study of short-circuit robustness of SiC MOSFETs, analysis of the failure modes and comparison with BJTs. Microelectron. Reliab. 2015, 55, 1708–1713. [Google Scholar] [CrossRef] [Green Version]
  11. Othman, D.; Lefebvre, S.; Berkani, M.; Khatir, Z.; Ibrahim, A.; Bouzourene, A. Robustness of 1.2 kV SiC MOSFET devices. Microelectron. Reliab. 2013, 53, 1735–1738. [Google Scholar] [CrossRef]
  12. Richardeau, F.; Boige, F. Circuit-Type modelling of SiC power MOSFET in short-circuit operation including selective fail-to-open and fail-to-short modes competition. In Proceedings of the 30th ESREF, Toulouse, France, 23–26 September 2019. [Google Scholar]
  13. Li, X.; Deng, X.; Zhu, H.; Wen, Y.; Li, X.; Zhang, Y.; Li, Z.; Zhang, B. Analysis of Short-Circuit Behavior for SiC MOSFETs with Various Circuit Characteristics. In Proceedings of the 2020 17th China International Forum on Solid State Lighting & 2020 International Forum on Wide Bandgap Semiconductors China (SSLChina: IFWS), Shenzhen, China, 23–25 November 2020; pp. 54–57. [Google Scholar] [CrossRef]
  14. Namai, M.; An, J.; Yano, H.; Iwamuro, N. Investigation of short-circuit failure mechanisms of SiC MOSFETs by varying DC bus voltage. Jpn. J. Appl. Phys. 2018, 57, 074102. [Google Scholar] [CrossRef]
  15. Pulvirenti, M.; Cavallaro, D.; Bentivegna, N.; Cascino, S.; Zanetti, E.; Saggio, M. Thermal Model Development for SiC MOSFETs Robustness Analysis under Repetitive Short Circuit Tests. In Proceedings of the 2020 22nd European Conference on Power Electronics and Applications (EPE’20 ECCE Europe), Lyon, France, 7–10 September 2020; pp. P.1–P.10. [Google Scholar] [CrossRef]
  16. Qin, H.; Dong, Y.; Xu, K.; Xu, H.; Fu, D.; Wang, S.; Zhao, C. A comprehensive study of the short-circuit characteristics of SiC MOSFETs. In Proceedings of the 2017 12th IEEE Conference on Industrial Electronics and Applications (ICIEA), Siem Reap, Cambodia, 18–20 June 2017; pp. 332–336. [Google Scholar] [CrossRef]
  17. Clemente, S. Transient thermal response of power semiconductors to short power pulses. IEEE Trans. Power Electron. 1993, 8, 337–341. [Google Scholar] [CrossRef]
  18. Boige, F.; Tremouilles, D.; Richardeau, F. Physical Origin of the Gate Current Surge During Short-Circuit Operation of SiC MOSFET. IEEE Electron Device Lett. 2019, 40, 666–669. [Google Scholar] [CrossRef]
  19. Richardeau, F.; Boige, F.; Lefebvre, S. Gate Leakage-Current, Damaged Gate and Open-Circuit Failure-Mode of Recent SiC Power Mosfet; IEEE: Nottingham, UK, 2018. [Google Scholar]
  20. Unger, C.; Pfost, M. Investigation of Gate and Drain Leakage Currents During the Short Circuit of SiC-MOSFETs. In Proceedings of the 2020 32nd (ISPSD), Vienna, Austria, 17–21 May 2020; pp. 238–241. [Google Scholar]
  21. Boige, F.; Richardeau, F.; Trémouilles, D.; Lefebvre, S.; Guibaud, G. Investigation on Damaged Planar-Oxide of 1200 V SiC Power MOSFETs in Non-Destructive Short-Circuit Operation; Elsevier: Amsterdam, The Netherlands, 2017. [Google Scholar]
  22. Liu, J.; Zhang, G.; Wang, B.; Li, W.; Wang, J. Gate Failure Physics of SiC MOSFETs Under Short-Circuit Stress. IEEE Electron Device Lett. 2020, 41, 103–106. [Google Scholar] [CrossRef]
  23. Boige, F. Caractérisation et Modélisation électrothermique Compacte étendue du MOSFET SiC en Régime Extrême de Fonctionnement Incluant ses Modes de Défaillance: Application à la Conception d’une Protection Intégrée au Plus Proche du Circuit de Commande. Ph.D. Thesis, Institut National Polytechnique de Toulouse, Toulouse, France, 2019. [Google Scholar]
  24. Cavallaro, D.; Pulvirenti, M.; Zanetti, E.; Saggio, M.G. Capability of SiC MOSFETs under Short-Circuit Tests and Development of a Thermal Model by Finite Element Analysis. MSF 2019, 963, 788–791. [Google Scholar] [CrossRef]
  25. Reigosa, P.D.; Iannuzzo, F.; Ceccarelli, L. Effect of short-circuit stress on the degradation of the SiO2 dielectric in SiC power MOSFETs. Microelectron. Reliab. 2018, 88–90, 577–583. [Google Scholar] [CrossRef]
  26. Chen, C.; Nguyen, T.; Labrousse, D.; Lefebvre, S.; Buttay, C.; Morel, H. New definition of critical energy for SiC MOSFET robustness under short circuit operations: The repetitive critical energy. Microelectron. Reliab. 2020, 114, 113839. [Google Scholar] [CrossRef]
  27. Castellazzi, A.; Richardeau, F.; Borghese, A.; Boige, F.; Fayyaz, A.; Irace, A.; Guibaud, G.; Chazal, V. Gate-damage accumulation and off-line recovery in SiC power MOSFETs with soft short-circuit failure mode. Microelectron. Reliab. 2020, 114, 113943. [Google Scholar] [CrossRef]
  28. Richardeau, F.; Boige, F.; Castellazzi, A.; Chazal, V.; Fayyaz, A.; Borghese, A.; Irace, A.; Guibaud, G. SiC MOSFETs soft and hard failure modes: Functional analysis and structural characterization. In Proceedings of the 2020 32nd International Symposium on Power Semiconductor Devices and ICs (ISPSD), Vienna, Austria, 13–18 September 2020; pp. 170–173. [Google Scholar] [CrossRef]
  29. Boige, F.; Richardeau, F.; Lefebvre, S.; Cousineau, M. SiC Power MOSFET in Short-Circuit Operation: Electro-Thermal Macro-Modelling Combining Physical and Numerical Approaches with Circuit-Type Implementation; Matcom: Amsterdam, The Netherlands, 2019; pp. 375–386. [Google Scholar]
Figure 1. (a) Schematic view of the power test bench used for short-circuit test. (b) V1-2018 Laplace Laboratory test-bench: Max 600 V, 400 A, 100 µs, 125 °C.
Figure 1. (a) Schematic view of the power test bench used for short-circuit test. (b) V1-2018 Laplace Laboratory test-bench: Max 600 V, 400 A, 100 µs, 125 °C.
Energies 14 07960 g001
Figure 2. (a) Oscilloscope waveforms at 1 µs of the DUT #3 under HSF, @VDS = 600 V, RG = 47 Ω, VDrv = −5 V/18 V. (b) Applied control SC signal test.
Figure 2. (a) Oscilloscope waveforms at 1 µs of the DUT #3 under HSF, @VDS = 600 V, RG = 47 Ω, VDrv = −5 V/18 V. (b) Applied control SC signal test.
Energies 14 07960 g002
Figure 3. Stacked incremental pulse duration oscilloscope waveforms from 0.5 µs to 4.5 µs with a step of 250 ns. @VDS = 600 V, RG = 47 Ω, VDrv = −5 V/18 V, DUT #3.
Figure 3. Stacked incremental pulse duration oscilloscope waveforms from 0.5 µs to 4.5 µs with a step of 250 ns. @VDS = 600 V, RG = 47 Ω, VDrv = −5 V/18 V, DUT #3.
Energies 14 07960 g003
Figure 4. Continuity of the stacked incremental pulse duration oscilloscope waveforms from 4.75 µs to 5.25 µs with a step of 250 ns. @VDS = 600 V, RG = 47 Ω, VDrv = −5 V/18 V, DUT #3.
Figure 4. Continuity of the stacked incremental pulse duration oscilloscope waveforms from 4.75 µs to 5.25 µs with a step of 250 ns. @VDS = 600 V, RG = 47 Ω, VDrv = −5 V/18 V, DUT #3.
Energies 14 07960 g004
Figure 5. Continuity of the stacked incremental pulse duration oscilloscope waveforms from 5.5 µs to 7.5 µs with a step of 250 ns. @VDS = 600 V, RG = 47 Ω, VDrv = −5 V/18 V, DUT #3, 1st SC test. (NB: gate damage at 5.25 µs).
Figure 5. Continuity of the stacked incremental pulse duration oscilloscope waveforms from 5.5 µs to 7.5 µs with a step of 250 ns. @VDS = 600 V, RG = 47 Ω, VDrv = −5 V/18 V, DUT #3, 1st SC test. (NB: gate damage at 5.25 µs).
Energies 14 07960 g005
Figure 6. Timeline illustration exhibiting the main pulses: G_L, gate_leakage, D_L, drain_leakage, G_D, gate-damage; FTS, pulse_final.
Figure 6. Timeline illustration exhibiting the main pulses: G_L, gate_leakage, D_L, drain_leakage, G_D, gate-damage; FTS, pulse_final.
Energies 14 07960 g006
Figure 7. Stacked incremental pulse duration oscilloscope waveforms from 5.25 µs to 17.75 µs with a step of 250 ns. @VDS = 600 V, RG = 47 Ω, VDrv = −5 V/18 V, DUT #4, 2nd SC test. (NB: gate damage at 5 µs).
Figure 7. Stacked incremental pulse duration oscilloscope waveforms from 5.25 µs to 17.75 µs with a step of 250 ns. @VDS = 600 V, RG = 47 Ω, VDrv = −5 V/18 V, DUT #4, 2nd SC test. (NB: gate damage at 5 µs).
Energies 14 07960 g007
Figure 8. Stacked incremental pulse duration oscilloscope waveforms from 1 µs to 16 µs with a step of 1 µs. @VDS = 600 V, RG = 47 Ω, VDrv = −5 V/10 V, DUT #8.
Figure 8. Stacked incremental pulse duration oscilloscope waveforms from 1 µs to 16 µs with a step of 1 µs. @VDS = 600 V, RG = 47 Ω, VDrv = −5 V/10 V, DUT #8.
Energies 14 07960 g008
Figure 9. Energy dissipated for each incremental pulse test, exhibiting the failure energy behavior. (a) DUT #3, (b) DUT #4, (c) DUT #6, @VDS = 600 V, RG = 47 Ω, VDrv = − 5 V/18 V.
Figure 9. Energy dissipated for each incremental pulse test, exhibiting the failure energy behavior. (a) DUT #3, (b) DUT #4, (c) DUT #6, @VDS = 600 V, RG = 47 Ω, VDrv = − 5 V/18 V.
Energies 14 07960 g009
Figure 10. Estimated behavior from the energy extrapolation.
Figure 10. Estimated behavior from the energy extrapolation.
Energies 14 07960 g010
Figure 11. Gate-crack impact on RDS-ON increase using VDS depolarization.
Figure 11. Gate-crack impact on RDS-ON increase using VDS depolarization.
Energies 14 07960 g011
Figure 12. RDSon (VGS) under VDrv depolarization @ID = 20 A, TCASE = 200 °C from Data-Sheet device SCTW40N120G2VAG; https://www.st.com (accessed on 25 November 2021).
Figure 12. RDSon (VGS) under VDrv depolarization @ID = 20 A, TCASE = 200 °C from Data-Sheet device SCTW40N120G2VAG; https://www.st.com (accessed on 25 November 2021).
Energies 14 07960 g012
Figure 13. 1D chip temperature simulation in short-circuit operation for VDS = 600 V, VDrv = −5 V/18 V, RG = 47 Ω, TPulse = 5 µs, Tcase = 25 °C.
Figure 13. 1D chip temperature simulation in short-circuit operation for VDS = 600 V, VDrv = −5 V/18 V, RG = 47 Ω, TPulse = 5 µs, Tcase = 25 °C.
Energies 14 07960 g013
Figure 14. Estimated temperature of the main chip layers under Comsol for the DUT #8 at VDS = 600 V, VDrv = −5 V/10 V and Tcase = 25 °C at TSC = 16 µs.
Figure 14. Estimated temperature of the main chip layers under Comsol for the DUT #8 at VDS = 600 V, VDrv = −5 V/10 V and Tcase = 25 °C at TSC = 16 µs.
Energies 14 07960 g014
Table 1. Different key parameters of several devices.
Table 1. Different key parameters of several devices.
DUT#2 @
25 °C
#3 @
25 °C
#4 @
25 °C
#6 @
25 °C
#16 @
125 °C
#19 @
125 °C
Test1st Test1st Test2nd Test1st Test1st Test1st Test
VGS_th [V]3.853.313.272.893.062.99
TPulse_G_L [µs]3.53.53.533.753.5
TPulse_D_L [µs]4.254.253.753.754.754.25
TPulse_G_D [µs]4.755.2554.754.755
TPulse_Final [µs]8.757.517.756.5NA5
ID_Max [A]260285285320245270
EFinal [mJ]NA606592581NA504
EFinal·cm−2 [J·cm−2]NA10.19.879.68NA8.40
Crack Delay [µs]2.4 µs5.5 µs5 µs4 µs>4.25 µs6.25 µs
Pulses Nb. before Damage18209181819
Failure TypeGD + FTSGD + FTSGD + FTSGD + FTSGD + FTSGD followed FTS
Table 2. Average key parameters between VDS and VGS depolarization compared to nominal operating conditions [15].
Table 2. Average key parameters between VDS and VGS depolarization compared to nominal operating conditions [15].
[15]VDS Depolarization MethodVGS Depolarization Method
ConditionsVBus = 800 V, RG = 4.7 Ω, VGS = −5 V/+18 VVBus = 600 V, RG = 47 Ω, VGS = −5 V/+18 VVBus = 600 V, RG = 47 Ω, VGS = −5 V/+10 V
Time [µs]T° [C]Time [µs]T° [C]Time [µs]T° [C]
TPulse_G_L--3.25917N.V (3)N.V (3)
TPulse_D_L--4.2510828.8977
TPulse_G_D--4.91207 (1),(2)9.61056 (1)
TPulse_Final1.621320 (1)9Destroyed161241 (2)
EG_D [mJ]-540-
EFinal [mJ]-600-
Failure modeFTSGD + FTSGD + FTO
N.V: none value, (1) first light gate-damage or failure temperature, (2) hard gate-damage temperature, (3) N.V: none value.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Barazi, Y.; Richardeau, F.; Jouha, W.; Reynes, J.-M. VDS and VGS Depolarization Effect on SiC MOSFET Short-Circuit Withstand Capability Considering Partial Safe Failure-Mode. Energies 2021, 14, 7960. https://doi.org/10.3390/en14237960

AMA Style

Barazi Y, Richardeau F, Jouha W, Reynes J-M. VDS and VGS Depolarization Effect on SiC MOSFET Short-Circuit Withstand Capability Considering Partial Safe Failure-Mode. Energies. 2021; 14(23):7960. https://doi.org/10.3390/en14237960

Chicago/Turabian Style

Barazi, Yazan, Frédéric Richardeau, Wadia Jouha, and Jean-Michel Reynes. 2021. "VDS and VGS Depolarization Effect on SiC MOSFET Short-Circuit Withstand Capability Considering Partial Safe Failure-Mode" Energies 14, no. 23: 7960. https://doi.org/10.3390/en14237960

APA Style

Barazi, Y., Richardeau, F., Jouha, W., & Reynes, J.-M. (2021). VDS and VGS Depolarization Effect on SiC MOSFET Short-Circuit Withstand Capability Considering Partial Safe Failure-Mode. Energies, 14(23), 7960. https://doi.org/10.3390/en14237960

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