In this Letter, we demonstrate a Pt/BaTiO3/Al0.58Ga0.42N lateral heterojunction diode with enhanced breakdown characteristics. By using BaTiO3, a high-k material, as a dielectric material between the anode and semiconductor, the peak electric field at the anode edge near the cathode was significantly reduced and an average breakdown field exceeding 8 MV/cm was achieved for devices with an anode to cathode spacing of <0.2 μm. In contrast, Pt/Al0.58Ga0.42N control Schottky diodes displayed an average breakdown field of ∼4 MV/cm for devices with similar dimensions. The use of a high-k dielectric can more effectively utilize the high breakdown fields in ultra-wide bandgap materials by proper management of the electric field. This demonstration thus provides a framework to realize ultra-scaled lateral devices with improved breakdown characteristics.

The AlxGa1−xN material system has garnered significant research interest due to its attractive material properties such as an extremely high critical breakdown field, predicted to be above 9 MV/cm for x ≥ 0.5, good low-field transport properties, and high electron saturation velocity comparable to GaN.1–4 As a result of these properties, AlxGa1−xN has a higher predicted Johnson Figure of Merit,

(JFOM=vsatFBR2π),

compared to GaN (vsat is the saturated electron velocity and FBR is the critical breakdown field). This can enable ultra-scaled devices, which can surpass the state of the art GaN-based amplifiers in terms of current density, breakdown voltage, and output power density at mm-wave and THz frequencies.2–4 In addition, the large critical breakdown field of AlGaN could also enable more favorable device characteristics for normally off devices used in power switching applications.5 

While various figures of merit suggest the excellent potential of these materials for future device technology, approaching the breakdown limit in an actual device remains a standing challenge. Material breakdown limits for wide-band gap semiconductors have been reached previously only with PN junctions.6 Achieving the theoretical material breakdown field in lateral diodes and field effect transistors is challenging due to non-uniform electric field distribution in the depletion region (which causes electric field peaking) and is further limited by tunneling breakdown at the Schottky gate or anode electrode. These effects are more significant in the case of ultra-wide band gap semiconductors, which have average breakdown fields estimated to be higher than 7 MV/cm. AlGaN devices that employ metal–semiconductor based Schottky anodes or gates are limited by the breakdown strength of the Schottky diodes, which typically have barrier heights of less than 2 eV. Breakdown fields up to 3.9 MV/cm have been demonstrated for Al0.70Ga0.30N/Al0.50Ga0.50N high electron mobility transistors (HEMTs) and 2.86 MV/cm for Schottky gate electrodes for Al0.70Ga0.30N metal semiconductor field effect transistors (MESFETs).7,8 However, these numbers are still significantly lower than the estimated material breakdown field (11 MV/cm) at this composition. Conventional dielectrics such as SiOX and Al2O3 inserted between the anode or gate metal and semiconductor do not offer a significantly better performance since they are limited by the breakdown of the metal–dielectric junction.

A solution to this problem, proposed recently, is to use extreme permittivity oxides inserted between the metal and semiconductor layers and between the gate and drain region.9,10 Extreme permittivity dielectrics enable improved breakdown due to two reasons. First, a high dielectric constant layer greatly reduces the peaking or non-uniformity in the lateral electric field between the gate to drain for lateral transistors or anode to cathode for lateral diodes and results in a higher breakdown voltage and higher average breakdown field due to a more uniform electric field distribution. This is illustrated in Fig. 1, which shows a comparison of the total electric field and potential distribution profiles between the anode and cathode of a heterojunction diode in the semiconductor (channel) near the semiconductor–dielectric interface (inset of Fig. 1) with a 25 nm low-k dielectric (εr = 10) and a high k dielectric (εr = 100) for an applied anode to cathode voltage of 150 V and an anode to cathode spacing of 0.2 μm. The channel material assumed here is AlGaN with a sheet charge density, ns, of 1.5 × 1013 cm−2. The addition of the high permittivity dielectric expands the depletion region due to the negative polarization charge induced by the field. Therefore, the voltage drop occurs over a larger gate-drain depletion width, thereby decreasing the field in the gate-drain region. In addition to reducing peak fields due to depletion, the high permittivity material below the gate or anode metal helps to reduce leakage, which suppresses gate or anode leakage related breakdown. This can be understood by considering the continuity of the displacement field across the dielectric/semiconductor junction, which suppresses the electric field within the high dielectric constant material.

FIG. 1.

Comparison of the total electric field and potential distribution between the anode and the cathode at the dielectric-channel interface for a lateral heterojunction diode (inset) with a 25 nm low-k dielectric (εr = 10) and high-k dielectric (εr = 100) for an applied voltage, VAC = −150 V. The dashed red line in the inset shows the cutline where the electric field and potential distribution have been simulated.

FIG. 1.

Comparison of the total electric field and potential distribution between the anode and the cathode at the dielectric-channel interface for a lateral heterojunction diode (inset) with a 25 nm low-k dielectric (εr = 10) and high-k dielectric (εr = 100) for an applied voltage, VAC = −150 V. The dashed red line in the inset shows the cutline where the electric field and potential distribution have been simulated.

Close modal

Barium titanate, BaTiO3, a perovskite oxide displaying a high dielectric constant, is a suitable candidate for such applications. In this report, we demonstrate a >8 MV/cm average breakdown field for BaTiO3/Al0.58Ga0.42N (predicted critical breakdown field of Al0.58Ga0.42N ∼9.7 MV/cm) heterojunction diodes, which represents an important step toward achieving a true material breakdown field in lateral devices. In comparison, control Schottky diodes displayed an average breakdown field of ∼4 MV/cm, which is significantly lower than the material breakdown field limit and representative of typical observations for Schottky-based devices.7,8

Low-pressure (LP) Metal-Organic Chemical Vapor Deposition (MOCVD) was used to grow the active layers for the reported device structures on 3 μm thick high quality AlN-(0001) sapphire templates. These template layers were also deposited by LPMOCVD at growth temperatures close to 1250 °C. Their RMS surface roughness was measured using an atomic-force microscope (AFM) to be 1.44 nm, as shown in Fig. 2(a). The epilayer consisted of a 500 nm thick undoped i-Al0.58Ga0.42N buffer layer grown pseudomorphically over the entire 2-in. diameter AlN template at a growth temperature close to 1100 °C. This was followed by the growth of a 60 nm thick [Si+]-doped n-type Al0.58Ga0.42N layer, with a doping concentration of 4 × 1018 cm−3 on top of the buffer layer as reported elsewhere.11 High-resolution x-ray diffraction (XRD) spectra (BEDE D1 High Resolution Triple Axis X-ray Diffraction System) [Fig. 2(b)] were used to determine the composition.

FIG. 2.

(a) AFM scan of the surface of the sample showed an RMS surface roughness of 1.44 nm and (b) XRD was used to confirm the composition of the MOCVD grown film.

FIG. 2.

(a) AFM scan of the surface of the sample showed an RMS surface roughness of 1.44 nm and (b) XRD was used to confirm the composition of the MOCVD grown film.

Close modal

Selective area regrowth by molecular beam epitaxy (MBE) was used to achieve Ohmic contact to the MOCVD grown epilayer. 500 nm of silicon dioxide (SiO2) was deposited and patterned on the MOCVD-grown channel layer to act as a hard mask for MBE contact regrowth. Heavily [Si+]-doped Al0.58Ga0.42N (50 nm) was then regrown in the selectively formed pits using MBE to completely fill the etched regions. This was capped with a heavily [Si+]-doped reverse Al-composition graded AlGaN layer of thickness of 50 nm and a doping of 1 × 1020 cm−3. The purpose of the reverse Al-composition graded n++ AlGaN layer is to minimize abrupt conduction band offsets and facilitate low resistance non-alloyed Ohmic contact formation.7,12

Following MBE contact regrowth, the SiO2 regrowth mask was removed using a diluted buffered oxide etch (BOE) solution, and a metal stack of Ti/Al/Ni/Au (20/120/30/100 nm) was deposited on the regrown contact regions via an e-beam evaporator to form Ohmic contacts. This was followed by device isolation by an inductively coupled plasma-reactive ion etching (ICP-RIE) system with an RIE power of 30 W and a pressure of 5 mTorr. BaTiO3 was deposited via RF sputtering at 630 °C in an oxygen ambient. The dielectric constant of the BaTiO3 films was estimated to be around 60. The dielectric constant of BaTiO3 is a function of the grain size of the deposited film. For the deposition conditions used, the deposited film is expected to be amorphous or nano-crystalline.13 The observed dielectric constant of the deposited films is consistent with previous reports.14 For the fabrication of the control samples, BaTiO3 was etched away completely using SF6 plasma ICP-RIE etching. For the heterojunction diodes, BaTiO3 was etched away in the contact regions only. Anodes were formed by the deposition of a Pt/Au (60/100 nm) metal stack via e-beam evaporation for both the control (Schottky) and the heterojunction diodes. The schematic of the fabricated heterojunction diode is shown in Fig. 3(a). Figure 3(b) shows the equilibrium energy-band diagrams calculated using a self-consistent one-dimensional Poisson–Schrödinger solver, under the anode, for the BaTiO3 sample.15 The conduction band offset and Schottky barrier height were estimated from the known electron affinity and work function values for the materials.16,17 The schematic and energy band-diagram of the control (Schottky) device, under the anode, are shown in Figs. 3(c) and 3(d), respectively.

FIG. 3.

(a) Schematic of the lateral Pt/BaTiO3/Al0.58Ga0.42N heterojunction diode (εr,BaTiO3 = 100), (b) band-diagram under the anode for the heterojunction diode, (c) schematic of the fabricated lateral Pt/Al0.58Ga0.42N Schottky diode, and (d) band-diagram under the anode for the Al0.58Ga0.42N Schottky diode.

FIG. 3.

(a) Schematic of the lateral Pt/BaTiO3/Al0.58Ga0.42N heterojunction diode (εr,BaTiO3 = 100), (b) band-diagram under the anode for the heterojunction diode, (c) schematic of the fabricated lateral Pt/Al0.58Ga0.42N Schottky diode, and (d) band-diagram under the anode for the Al0.58Ga0.42N Schottky diode.

Close modal

Hall measurements were performed on four-terminal ungated van der Pauw (VDP) structures without BaTiO3 and indicated a sheet resistance of 6.2 kΩ/□, a Hall mobility of 65 cm2/V s, and a sheet charge density of 1.55 × 1013 cm−2. Electrical characterization of the devices was performed using an Agilent B1500 parameter analyzer. The inset of Fig. 4 shows the capacitance voltage characteristics measured on a 16 μm × 100 μm heterojunction diode. The integrated charge density, ns, was found to be 1.06 × 1013 cm−2. The measured charge density and the simulated charge density are shown in Fig. 4. The extracted doping density is found to be 4 × 1018 cm−3, as expected, near the top of the structure, and the lower density at larger depletion depths is likely due to errors introduced by capacitor loss.

FIG. 4.

Comparison of the measured charge profile to the simulated charge profile of the BaTiO3/Al0.58Ga0.42N diode; (inset) capacitance–voltage measurement results from BaTiO3/Al0.58Ga0.42N diodes.

FIG. 4.

Comparison of the measured charge profile to the simulated charge profile of the BaTiO3/Al0.58Ga0.42N diode; (inset) capacitance–voltage measurement results from BaTiO3/Al0.58Ga0.42N diodes.

Close modal

Figure 5(a) shows the forward characteristics of a Schottky diode and a heterojunction diode with an anode to cathode spacing of 700 nm. The turn-on voltages of the Schottky diode and the heterojunction diode are 0.9 V and 1.5 V, respectively. Both these devices display a significantly lower turn-on voltage compared to what is typically observed for AlGaN PN/PIN diodes. This suggests that the BaTiO3/AlGaN interface presents relatively small barriers to electron transport, in contrast to conventional metal-oxide-semiconductor junctions where a large voltage would be needed to achieve turn-on of the current. The differential on-resistance (RON), calculated from the forward IV characteristics, was 31 mΩ cm2 for the Schottky diode and 56 mΩ cm2 for the heterojunction diode. The increase in RON of the heterojunction device is likely due to the presence of the BaTiO3 layer, which is nominally undoped. Breakdown characteristics were evaluated for several randomly selected Schottky and heterojunction diodes. Figure 5(b) shows the breakdown characteristics for a Schottky diode and a heterojunction diode, with anode to cathode spacings of 0.19 μm and 0.18 μm, respectively, as measured using a scanning electron microscope (SEM). While the Schottky diode breaks down at an anode-to-cathode voltage of 42 V, the heterojunction diode breaks down at a much higher voltage of 155 V, despite both devices having a similar anode to cathode spacing. The improvement in the average breakdown field is attributed to a reduced peak field and the suppression of anode leakage current due to the BaTiO3 layer. Assuming a total depleted sheet charge density of 1.55 × 1013 cm−2, the peak vertical field (see the supplementary material for details) was estimated to be 3.15 MV/cm. Assuming a pinch-off voltage of −12 V, the average lateral field of the Schottky diode was estimated to be 1.58 MV/cm (see the supplementary material), which gives a total average field of 3.5 MV/cm (see the supplementary material). The peak vertical field and the average lateral field of the heterojunction diode were estimated to be 3.15 MV/cm and 7.94 MV/cm, respectively, leading to an average total breakdown electric field of 8.5 MV/cm. The vertical field should be the same in both devices since the total number of ionized donors in the depletion region should be the same for both devices. Multiple randomly chosen devices were measured for both the Schottky diodes and the heterojunction diodes, and the results are shown in Fig. 6(a). In general, for the heterojunction diodes, for an anode to cathode spacing in the range of 0.2–0.3 μm, the average breakdown fields observed were in the range of 6–8 MV/cm. Figure 6(b) shows the highest reported average breakdown fields for various semiconductor materials as a function of the bandgap. The value reported in this work, 8.5 MV/cm, is the highest reported experimental average breakdown field for any semiconductor material.

FIG. 5.

Comparison between the Schottky barrier diode and the BaTiO3/Al0.58Ga0.42N heterojunction diode: (a) forward IV characteristics and (b) breakdown performance; (inset) anode to cathode spacing for the heterojunction diode with VBR = −155 V, as measured using a scanning electron microscope (SEM).

FIG. 5.

Comparison between the Schottky barrier diode and the BaTiO3/Al0.58Ga0.42N heterojunction diode: (a) forward IV characteristics and (b) breakdown performance; (inset) anode to cathode spacing for the heterojunction diode with VBR = −155 V, as measured using a scanning electron microscope (SEM).

Close modal
FIG. 6.

(a) Breakdown performance for lateral Pt/BaTiO3/Al0.58Ga0.42N heterojunction diodes and AlGaN Schottky barrier devices for LAC < 500 nm and (b) comparison of experimental average breakdown fields in various semiconductor devices achieved to date6,8,9,18–23 (pink: GaN reports, green: Ga2O3 reports, red: AlGaN reports, and brown: diamond reports). The dashed orange line shows the phenomenological fit of the critical breakdown field (FBR) as a function of bandgap (EG) for direct bandgap semiconductors.1 The solid black circles represent the predicted critical breakdown fields for various materials.20,22,24–26

FIG. 6.

(a) Breakdown performance for lateral Pt/BaTiO3/Al0.58Ga0.42N heterojunction diodes and AlGaN Schottky barrier devices for LAC < 500 nm and (b) comparison of experimental average breakdown fields in various semiconductor devices achieved to date6,8,9,18–23 (pink: GaN reports, green: Ga2O3 reports, red: AlGaN reports, and brown: diamond reports). The dashed orange line shows the phenomenological fit of the critical breakdown field (FBR) as a function of bandgap (EG) for direct bandgap semiconductors.1 The solid black circles represent the predicted critical breakdown fields for various materials.20,22,24–26

Close modal

The breakdown voltage achieved for the diode with an average breakdown field of 8.5 MV/cm was 155 V, which is lower than what is expected from the theoretically estimated maximum breakdown field of 9.7 MV/cm and a flat electric field distribution. Therefore, there exists room for optimization of this BaTiO3/AlGaN heterostructure design including optimization of the BaTiO3 thickness and growth conditions together with the exploration of potential integration of field plate structures.

One important benefit of using BaTiO3 as a gate or anode dielectric is that it enables the fabrication of ultra-scaled devices while meeting the breakdown voltage requirements due to the much higher breakdown fields achievable in these structures. This is particularly beneficial for low mobility materials like AlGaN and Ga2O3, which require scaled devices for high performance RF electronics to compensate for the higher sheet resistance compared to GaN. The presence of the high-k dielectric also leads to scaling of the gate to channel distance, which can result in lower output conductance and higher transconductance.10 However, it should be noted that in such a heterojunction with a high-k dielectric, the drain depletion region is longer than the case with a low-k dielectric and a correspondingly higher gate to channel capacitance due to the high-permittivity region between the gate and drain. As a result, the peak cutoff frequency is expected to be lower than a device with a low-k barrier.10 It is therefore important to properly optimize the thickness of the high-permittivity region to manage the trade-off between field management and gate–drain capacitance. Scaling, enabled by the high-k heterostructure devices, is also expected to play an important role in next generation power electronics, whereby the on-resistance of the device can be significantly reduced due to shortened device dimensions while at the same time delivering higher output power more efficiently.

In conclusion, we have demonstrated Pt/BaTiO3/Al0.58Ga0.42N lateral heterojunction diodes with enhanced breakdown characteristics. By using BaTiO3 as an anode and an anode to cathode region dielectric, the peak field at the anode edge was significantly reduced and an average breakdown field exceeding 8 MV/cm was achieved for devices with an anode to cathode spacing of <0.2 μm. In contrast, Pt/Al0.58Ga0.42N control Schottky diodes displayed an average breakdown field of ∼4 MV/cm for devices with similar dimensions. The use of a high-k dielectric can more effectively utilize the high breakdown fields in ultra-wide bandgap materials by proper electric field management. This demonstration thus provides a framework to realize ultra-scaled lateral devices with improved breakdown characteristics.

See the supplementary material for a description of the estimation method of the average breakdown electric field.

The authors acknowledge funding from the Air Force Office of Scientific Research (AFOSR Grant No. FA9550-17-1-0227, Program Manager Kenneth Goretta) and the DARPA DREaM program (No. ONR N00014-18-1-2033, Program Manager Dr. Young-Kai Chen, monitored by the Office of Naval Research, Program Manager Dr. Paul Maki).

1.
J. L.
Hudgins
,
G. S.
Simin
,
E.
Santi
, and
M. A.
Khan
,
IEEE Trans. Power Electron.
18
(
3
),
907
(
2003
).
2.
A. F. M.
Anwar
,
S.
Wu
, and
R. T.
Webster
,
IEEE Trans. Electron Devices
48
(
3
),
567
(
2001
).
3.
M.
Farahmand
,
C.
Garetto
,
E.
Bellotti
,
K. F.
Brennan
,
M.
Goano
,
E.
Ghillino
,
G.
Ghione
,
J. D.
Albrecht
, and
P. P.
Ruden
,
IEEE Trans. Electron Devices
48
(
3
),
535
(
2001
).
4.
T.
Razzak
,
H.
Xue
,
Z.
Xia
,
S.
Hwang
,
A.
Khan
,
W.
Lu
, and
S.
Rajan
, paper presented at the 2018 IEEE MTT-S International Microwave Workshop Series on Advanced Materials and Processes for RF and THz Applications (IMWS-AMP),
2018
.
5.
S.
Bajaj
,
T.-H.
Hung
,
F.
Akyol
,
D.
Nath
, and
S.
Rajan
,
Appl. Phys. Lett.
105
(
26
),
263503
(
2014
).
6.
M.
Qi
,
K.
Nomoto
,
M.
Zhu
,
Z.
Hu
,
Y.
Zhao
,
V.
Protasenko
,
B.
Song
,
X.
Yan
,
G.
Li
,
J.
Verma
,
S.
Bader
,
P.
Fay
,
H. G.
Xing
, and
D.
Jena
,
Appl. Phys. Lett.
107
(
23
),
232101
(
2015
).
7.
T.
Razzak
,
S.
Hwang
,
A.
Coleman
,
H.
Xue
,
S. H.
Sohel
,
S.
Bajaj
,
Y.
Zhang
,
W.
Lu
,
A.
Khan
, and
S.
Rajan
,
Appl. Phys. Lett.
115
(
4
),
043502
(
2019
).
8.
H.
Xue
,
S.
Hwang
,
T.
Razzak
,
C.
Lee
,
G. C.
Ortiz
,
Z.
Xia
,
S.
Hasan Sohel
,
J.
Hwang
,
S.
Rajan
,
A.
Khan
, and
W.
Lu
,
Solid-State Electron.
164
,
107696
(
2020
).
9.
Z.
Xia
,
C.
Wang
,
H.
Chandrasekar
,
W.
Moore
,
A.
Lee
,
N. K.
Kalarickal
,
F.
Yang
, and
S.
Rajan
, paper presented at the 2019 76th Device Research Conference (DRC),
2018
.
10.
Z.
Xia
,
C.
Wang
,
N. K.
Kalarickal
,
S.
Stemmer
, and
S.
Rajan
,
IEEE Trans. Electron Devices
66
(
2
),
896
(
2019
).
11.
S.
Muhtadi
,
S.
Hwang
,
A.
Coleman
,
F.
Asif
,
A.
Lunev
,
M. V. S.
Chandrashekhar
, and
A.
Khan
,
Appl. Phys. Lett.
110
(
19
),
193501
(
2017
).
12.
S.
Bajaj
,
F.
Akyol
,
S.
Krishnamoorthy
,
Y.
Zhang
, and
S.
Rajan
,
Appl. Phys. Lett.
109
(
13
),
133508
(
2016
).
13.
N.-Y.
Lee
,
T.
Sekine
,
Y.
Ito
, and
K.
Uchino
,
Jpn. J. Appl. Phys., Part 1
33
(
3R
),
1484
(
1994
).
14.
Z. Q.
Shi
,
Q. X.
Jia
, and
W. A.
Anderson
,
J. Electron. Mater.
20
(
11
),
939
(
1991
).
15.
M.
Grundmann
, BandEng: Poisson–Schrödinger Solver Software 2004, http://my.ece.ucsb.edu/mgrundmann/bandeng/.
16.
H. S.
Nalwa
,
Handbook of Thin Film Materials: Ferroelectric and Dielectric Thin Films
(
Academic Press
,
2002
).
17.
J.
Li
,
N.
Li
,
C.
Ge
,
H.
Huang
,
Y.
Sun
,
P.
Gao
,
M.
He
,
C.
Wang
,
G.
Yang
, and
K.
Jin
,
iScience
16
,
368
(
2019
).
18.
A. A.
Allerman
,
A. M.
Armstrong
,
A. J.
Fischer
,
J. R.
Dickerson
,
M. H.
Crawford
,
M. P.
King
,
M. W.
Moseley
,
J. J.
Wierer
, and
R. J.
Kaplar
,
Electron. Lett.
52
(
15
),
1319
(
2016
).
19.
S.
Bajaj
,
A.
Allerman
,
A.
Armstrong
,
T.
Razzak
,
V.
Talesara
,
W.
Sun
,
S. H.
Sohel
,
Y.
Zhang
,
W.
Lu
, and
A. R.
Arehart
,
IEEE Electron Device Lett.
39
(
2
),
256
(
2018
).
20.
A. J.
Green
,
K. D.
Chabak
,
E. R.
Heller
,
R. C.
Fitch
,
M.
Baldini
,
A.
Fiedler
,
K.
Irmscher
,
G.
Wagner
,
Z.
Galazka
, and
S. E.
Tetlak
,
IEEE Electron Device Lett.
37
(
7
),
902
(
2016
).
21.
A.
Nishikawa
,
K.
Kumakura
,
T.
Akasaka
, and
T.
Makimoto
,
Appl. Phys. Lett.
88
(
17
),
173508
(
2006
).
22.
A.
Nishikawa
,
K.
Kumakura
, and
T.
Makimoto
,
Jpn. J. Appl. Phys., Part 1
46
(
4S
),
2316
(
2007
).
23.
P.-N.
Volpe
,
P.
Muret
,
J.
Pernot
,
F.
Omnès
,
T.
Teraji
,
Y.
Koide
,
F.
Jomard
,
D.
Planson
,
P.
Brosselard
, and
N.
Dheilly
,
Appl. Phys. Lett.
97
(
22
),
223501
(
2010
).
24.
C. F.
Zhe
,
Silicon Carbide: Materials, Processing & Devices
(
CRC Press
,
2003
).
25.
R.
Kaplar
,
A. A.
Allerman
,
A.
Armstrong
,
J. J.
Wierer
,
A. J.
Fischer
,
A. G.
Baca
, and
J. C.
Neely
, No. SAND2015-5644C,
Sandia National Lab.
,
Albuquerque
, NM.
26.
K.
Driche
,
H.
Umezawa
,
N.
Rouger
,
G.
Chicot
, and
E.
Gheeraert
,
Jpn. J. Appl. Phys., Part 1
56
(
4S
),
04CR12
(
2017
).

Supplementary Material