1 Introduction
Nuclear batteries have attracted the interest of researchers as far back as the 1900s [1]. In recent decades, betavoltaic batteries have been intensively demonstrated. The structure mainly contains a beta source and an energy converter such as a p-n junction, a p-i-n junction, or a Schottky barrier diode [1-20]. In these energy converters, beta particles interact with matter and the ionizing radiation generates electron–hole pairs. The electron–hole pairs, which are generated in or near the depletion region, can be separated and eventually swept out of this region by the built-in potential barrier. Then they are collected for the radiation-induced current formation, thereby converting the decay energy into electrical energy directly [1, 10, 13]. The performance of a betavoltaic cell is determined by the radioisotope, radiation transport properties of beta particles in matter, and the energy converter. Compared with micro-chemical cells and micro-solar cells, the advantages of betavoltaic batteries mainly include easy integration on a small scale, stable performance, long service lifetime, high energy density, self-sustainability, and little maintenance [2-20]. Therefore, such batteries have become a promising micro-power source for micro-electronic mechanical systems (MEMSs), medical applications such as pacemaker power, and applications in execrable ambient such as space exploration [2-20].
Previous research has shown that SiC material and a 63Ni source are suitable for p-n-junction-based betavoltaic batteries. For example, in 2006, a micro-battery with a 1 mCi 63Ni source on an n+p SiC structure was reported. An open-circuit voltage of 0.72 V and a short-circuit current density of 16.8 nA/cm2 were measured [3]. In 2015, the Monte Carlo method was used to investigate a 4H-SiC cell with a 63Ni source. The maximum device conversion efficiency was 25.04% when the thickness of the source was 6 μm [9]. In 2016, another 4H-SiC-63Ni cell was studied by using the Monte Carlo method and a Synopsys Medici device simulator. The simulation results showed that the total conversion efficiency decreased from 13.65% to 1.52% when the thickness of the source increased from 0.1 to 5 μm [10].
In this research, we analyze the main factors that affect the energy conversion efficiencies of 4H-SiC-63Ni p-n-junction-based batteries. Based on the simulation results, theoretical maximum electrical properties and energy conversion efficiencies are calculated. Finally, the light p-type region is selected as the emitter layer and the high n-type region is used as the base layer. Then the p-n junction is designed to match the depletion region with the main energy deposition region of the 63Ni source in 4H-SiC material. An optimal design that has a maximum output power of 0.36 μW/cm2 and a device conversion efficiency of 23.5% is then proposed. These results are helpful for designing practical cells. Although 4H-SiC material and a 63Ni source are used in this study, the simulated and calculated model can be extended to explore other semiconductor materials and beta sources for betavoltaic batteries.
2 Selection of materials and sources
4H-SiC is a wide-band-gap semiconductor material. The band gap and the intrinsic carrier concentration are ∼3.25 eV and ∼9.5×10-9/cm3 at 300 K, respectively [1]. The wide band gap and the low intrinsic carrier concentration contribute to the high energy conversion efficiency and high radiation resistance of nuclear batteries [1, 10]. In the interaction of beta particles with the target material, electron–hole pairs are generated as well as heat. The fraction of energy that goes into electron–hole pair formation is calculated as the band gap energy Eg divided by the mean ionization energy Eehp. This ratio represents the ultimate conversion efficiency ≤ult of the electron-hole pairs production efficiency [1]. The mean ionization energy is confirmed as Eehp = 2.8Eg + 0.5 in Ref. [9]. As shown in Fig. 1, the ultimate conversion efficiency increases and then nearly reaches its saturation value with increasing band gap. For 4H-SiC material, the ultimate conversion efficiency is 33.8%. Thus, 4H-SiC material is fit for betavoltaic cells.
-201811/1001-8042-29-11-017/alternativeImage/1001-8042-29-11-017-F001.jpg)
The sources for betavoltaic batteries are listed in Table 1 [1, 2, 8]. In general, high incident energy results in high electrical power in a certain volume. The average decay energy Eave of a 3H source is low for high-energy devices. A 35S source is not suitable for long-term applications owing to its short half-life T1/2. A 63Ni source is the main candidate for betavoltaic cells because of its pure beta radiation, long half-life, and the moderate average decay energy. Moreover, the maximum decay energy βmax of a 63Ni source is below the threshold energy of radiation damage in most semiconductor materials [4, 5, 8]. For a 85Kr source and a 147Pm source, their maximum decay energies are high and they are both gamma emitters (85Kr: 0.4% for γ emission and 147Pm: 0.003% for γ emission) [1, 8]. The high-energy beta particles and the gamma rays may badly damage the lattice structure of semiconductor materials. A 90Y source is also a beta emitter and it is produced from the decay of a 90Sr source. Therefore, a 90Sr source is a mixture of 90Sr and 90Y sources, and the average decay energy is 1.1 MeV [8]. This decay energy is excessively high for 4H-SiC material. Therefore, a pure 63Ni source is suitable for 4H-SiC betavoltaic batteries. In addition, a 63Ni source is easier to handle because of its solid-metal form.
Isotope | 3H | 35S | 63Ni | 85Kr | 90Sr | 90Y | 147Pm |
---|---|---|---|---|---|---|---|
T1/2 (years) | 12.3 | 0.24 | 100.2 | 10.7 | 28.8 | 0.007 | 2.6 |
βmax (keV) | 18.6 | 168 | 66.7 | 670 | 546 | 2280 | 225 |
Eave (keV) | 5.7 | 48.8 | 17.4 | 251 | 195.8 | 945 | 62 |
3 Simulation results and discussion
The Monte Carlo method is used to simulate the self-absorption effect of the 63Ni source, the backscattering process, and the transport of beta particles in 4H-SiC material. Then we analyze the main factors that affect the energy conversion efficiencies of 4H-SiC-63Ni cells.
3.1 Self-absorption effect of source
Beta particles dissipate their energy within the source itself and this phenomenon is referred as self-absorption. The self-absorption loss ≤sel is defined as the energy deposited in the source divided by the total energy produced in the source. In the simulated model, the energy spectrum of the 63Ni isotope is shown in Fig. 2[19]. Its geometrical structure is cuboid and its cross-sectional area is 1×1 cm2.
-201811/1001-8042-29-11-017/alternativeImage/1001-8042-29-11-017-F002.jpg)
In Fig. 3, it can be seen that, when the thickness of the 63Ni layer increases from 2×10-3 to 10 μm, the apparent power density increases from 1×10-2 to 3.24 μW/cm2 and the self-absorption loss increases from 0.036% to 93.7%. Moreover, when the thickness of the source is 2 μm, the apparent power density is 3.1 μW/cm2 and the self-absorption loss is 70%. For 0.004-, 0.1-, 0.3-, 1-, 1.6-, 2-, 4-, 6-, and 10-μm-thick 63Ni sources, the emitted energy spectra are shown in Fig. 4. It can be seen that the peaks of these emitted energy spectra move toward the higher energy region and then mostly keep fixed as the thickness of the source increases. Finally, as shown in Fig. 5, when the thickness of the 63Ni layer increases from 2×10-3 to 10 μm, the average energy of the emitted energy spectrum increases from 17.5 keV to a maximum value of 21.24 keV at∼2 μm. These results can be explained as follows: When the source is thin, low-energy beta particles are easily absorbed while high-energy beta particles have a chance to escape from the source. However, if the source is thick enough, even high-energy beta particles in deep positions can barely emit from the beta source.
-201811/1001-8042-29-11-017/alternativeImage/1001-8042-29-11-017-F003.jpg)
-201811/1001-8042-29-11-017/alternativeImage/1001-8042-29-11-017-F004.jpg)
-201811/1001-8042-29-11-017/alternativeImage/1001-8042-29-11-017-F005.jpg)
3.2 Directional loss
The typical prototype of a p-n-junction-based betavoltaic cell is illustrated in Fig. 6[6, 12, 14]. In such a structure, less than half of the emitted beta particles can reach the surface of the energy converter. Then the directional loss ≤dir is defined as the number of beta particles that do not reach the surface of the energy converter divided by the total number of beta particles emitted from the source [9]. In a practical cell, the same energy converter is usually set on the other side of the source. Therefore, a stacked structure is the main method used to decrease directional loss [7].
-201811/1001-8042-29-11-017/alternativeImage/1001-8042-29-11-017-F006.jpg)
3.3 External interaction loss
The external interaction loss is described as the energy loss between the source and the energy converter [9]. One part of the external interaction loss is the energy deposited in the dead layer, which mainly contains the metal electrode layer, the narrow air layer, and the passivation layer on the surface of cell. Generally, the scale of the dead layer is tens of nanometers and the energy deposited in this region can be neglected [1]. Another part of the external interaction loss is the energy backscattering loss on the surface, which is a significant energy loss in betavoltaic cells. The energy backscattering coefficient ≤sca is defined as the reflected energy divided by the total energy entering the energy converter. A theoretical analysis shows that the low atomic number of 4H-SiC material and the high energy of the 63Ni source contribute to the small energy backscattering coefficient. In this study, the energy converter consists only of 4H-SiC material. Its geometrical structure is cuboid and the junction area is 1×1 cm2. Moreover, the beta source is placed close to the surface of the energy converter and the whole device is set in vacuum. In Fig. 7, it can be seen that, when the thickness of the 63Ni layer increases from 2×10-3 to 10 μm, the energy backscattering coefficient on the surface of the 4H-SiC material device decreases from 28.77% to its minimum value of 20.15% at ~ 2 μm.
-201811/1001-8042-29-11-017/alternativeImage/1001-8042-29-11-017-F007.jpg)
3.4 Energy distribution and internal interaction loss
The energy distribution of beta sources in matter can be obtained by using the Monte Carlo method. Here, for 0.1-, 0.3-, 1-, 1.6-, 2-, 4-, 6-, and 10-μm-thick 63Ni sources, the calculated results in Fig. 8 show that the power density deposition per unit length nearly exponentially decreases with the penetration depth increasing in 4H-SiC material. The results indicate that the maximum energy deposition per unit length takes place at the surface and that the energy of the beta source is quickly absorbed in matter. Further calculation shows that the width of the main energy deposition region, which absorbs >99% of the total energy deposited in 4H-SiC material, approaches 8 μm and is nearly identical for 63Ni layers of different thicknesses. This can be explained by the fact that the emitted energy spectra are similar for various source thicknesses. Therefore, a thin 4H-SiC membrane close to the source is needed to absorb the energy of the 63Ni source.
-201811/1001-8042-29-11-017/alternativeImage/1001-8042-29-11-017-F008.jpg)
Based on the beta source energy distribution law in matter, it is concluded that the p-n junction, the p-i-n junction, and the Schottky barrier diode are appropriate for betavoltaic cells, because an internal electrical field exists and their depletion regions can be designed to match the main region of ionizing radiation. Furthermore, the simulation results show that the maximum penetration depth of the 63Ni source in 4H-SiC material is nearly 23 μm and is independent of the source thickness, because the maximum penetration depth of beta particles in matter is determined by the maximum decay energy of the beta source.
In the p-n-junction-based cell, the electron–hole pairs generated in the depletion region are eventually separated and then collected for the radiation-induced current formation. The charge collection efficiency (CCE) in the depletion region is 100% [1]. However, the electron–hole pairs generated outside the depletion region, which must diffuse into the depletion region, can be separated and collected. Therefore, the depletion region and its nearby region, which is about one diffusion length of minority carrier away from the edge of depletion region, are referred to as the active region in a cell. The charge collection efficiency outside the active region can be neglected [1, 12, 15]. The energy transport efficiency ≤dpl can then be defined as the energy deposited in the active region (Pdep) divided by the apparent energy of the source (Psur) entering the energy converter. Moreover, the energy deposition outside the active region is described as the internal interaction loss in the cell.
The total energy deposited in the energy converter represents the maximum available energy for electron–hole pairs generation. Based on the simulation results, when the thickness of the 63Ni layer increases from 2×10-3 to 10 μm, the total energy deposited in 4H-SiC material increases from 3.71×10-3 to 1.53 μW/cm2. Moreover, when the thickness of the source is 2 μm, the total power density deposited in 4H-SiC material is 1.46 μW/cm2. These results are shown in Fig. 9.
-201811/1001-8042-29-11-017/alternativeImage/1001-8042-29-11-017-F009.jpg)
Further analysis shows that the energy deposited in 4H-SiC material is greater than the difference between the incident energy and the backscattering energy on the surface. This result may be attributed to multiple backscattering of beta particles on the source and device surfaces. Finally, if the energy deposited in 4H-SiC material is totally in the active region of the cell, the energy transport efficiency increases from 72% to 94% and it nearly reaches a saturation value at 1 μm when the thickness of the 63Ni layer increases from 2×10-3 to 10 μm. The result is plotted in Fig. 9.
4 Optimal cell design
4.1 Theoretical electrical properties
In a betavoltaic cell, the short-circuit current density Jsc is given by [1]
where CCE is the electron–hole pair collection efficiency, q is the unit electron charge, Pdep is the energy deposited in the active region, S is the junction area, and Eehp is the mean ionization energy. The open-circuit voltage Voc is given by [9]
where kB is Boltzmann’s constant, T is the absolute temperature, and J0 is the minimum ideal reverse saturation current density, which is confirmed as J0 = 1.5×105 exp(-Eg/kBT) in A/cm2 in Ref. [9]. During the calculations, we assume that the total energy deposited in 4H-SiC material is in the active region and that the charge collection efficiency is 100%. As shown in Fig. 10, when the thickness of the 63Ni layer increases from 2×10-3 to 10 μm, the short-circuit current density increases from 3.87×10-4 to 0.16 μA/cm2 and the open-circuit voltage increases from 2.38 to 2.53 V at 300 K.
-201811/1001-8042-29-11-017/alternativeImage/1001-8042-29-11-017-F010.jpg)
In Fig. 10, compared with the short-circuit current density, it can be seen that the open-circuit voltage is insensitive to the thickness of the 63Ni layer. The main reason for this is that the minimum ideal reverse saturation current is much smaller than the short-circuit current. Experimental tests have shown that the open-circuit voltage is mainly affected by the reverse saturation current [8]. In a p-n-junction-based energy converter, the reverse saturation current is mainly determined by the physical parameters of the material and the semiconductor manufacturing technology [1]. Moreover, the more radiation damage that exists, the greater the number of internal defects is created and then the reverse saturation current increases. Finally, the ratio, which is the open-circuit voltage to the band gap, is defined as the driving potential efficiency ≤p of the cell [1]. The calculated result shows that, when the thickness of the 63Ni layer increases from 2×10-3 to 10 μm, the theoretical maximum driving potential efficiency increases from 73.2% to 77.8% and it nearly reaches its maximum value at 2 μm in the 4H-SiC-63Ni cells.
The net current density J of the cell, which is in the reverse-biased direction, is given by [22]
where V is the voltage across the resistive load. For 0.1-, 0.3-, 1-, 1.6-, 2-, 4-, 6-, and 10-μm-thick 63Ni sources, the current density–voltage characteristics of 4H-SiC cells are shown in Fig. 11.
-201811/1001-8042-29-11-017/alternativeImage/1001-8042-29-11-017-F011.jpg)
The power density P delivered to the resistive load is given by
Based on Eq. (4), the maximum output power Pmax of the cell is given by [22]
where Vm is the voltage across the resistive load when P = Pmax. Then the fill factor FF is given by
As can be seen in Fig. 12, when the thickness of the 63Ni layer increases from 2×10-3 to 10 μm, the theoretical upper limit of the maximum output power density increases from 8.7×10-4 to 0.38 μW/cm2 and the fill factor increases from 94.04% to 94.32%. Moreover, the maximum output power density nearly reaches its saturation value when the thickness of the source is 2 μm and the fill factor nearly reaches its maximum value at 1 μm.
-201811/1001-8042-29-11-017/alternativeImage/1001-8042-29-11-017-F012.jpg)
4.2 Theoretical energy conversion efficiencies
There are usually two types of energy conversion efficiencies in betavoltaic batteries. One is the device conversion efficiency ≤dev, which is calculated by using [1]
where ≤dpl is the energy transport efficiency of the source, ≤ult is the ultimate conversion efficiency of the electron-hole pairs production efficiency, ≤p is the driving potential efficiency of the cell, and FF is the fill factor. Another type is the total conversion efficiency ≤total, which can be separated into the utilization efficiency of the source and the device conversion efficiency of the cell. The utilization efficiency of the source is mainly determined by the self-absorption loss itself and the directional loss of the cell. To minimize the directional loss, the same energy converters are usually placed on both the upper and base sides of the source in a practical cell. Therefore, the total conversion efficiency in this study is calculated by using the maximum output power of one energy converter divided by half the total power produced in the source[1, 8]:
where ≤sel is the self-absorption loss of the source.
Now, the energy conversion efficiencies of 4H-SiC-63Ni batteries can be calculated. It is shown that, when the thickness of the 63Ni layer increases from 2×10-3 to 10 μm, the theoretical maximum device conversion efficiency increases from 16.77% to the saturation value of 23.51% at ~2 μm and the total conversion efficiency decreases from 16.73% to 1.48%. The results are plotted in Fig. 13. The theoretical energy conversion efficiencies in Fig. 13 are nearly in accordance with the results reported in Refs. [9,10]. Although the energy conversion efficiencies measured in the fabricated cells are much lower than the theoretical calculated values [2,8,13,16], the results achieved in this study are helpful for optimizing the design of 4H-SiC-63Ni cells. Moreover, the semiconductor manufacturing technology and the isotope loading technique need further development. Meanwhile, radiation management is also an important issue.
-201811/1001-8042-29-11-017/alternativeImage/1001-8042-29-11-017-F013.jpg)
4.3 An optimal design of doping concentrations
In p-n-junction-based batteries, the p-type region is set as the base layer in solar cells. However, the p-type region is selected as the emitter layer and the n-type region is used as the base layer in betavoltaic cells. There are two reasons for this. The first reason is that the main energy deposition region of a beta source in matter is close to the surface of the emitter layer and it is different from that of the solar cells. The second reason is that the diffusion length of the electron is much greater than that of a hole under the same excess carrier density. In general, the diffusion length of an electron is twice or more greater than that of a hole under the same excess carrier density [22]. Moreover, it is known that the depletion region is usually limited to a few micrometers, which is smaller than the main energy deposition region in a fabricated betavoltaic cell [1]. Therefore, the energy deposited outside the depletion region is important for radiation-induced current formation. The design of the p-type region as an emitter layer contributes to improve the charge collection efficiency outside the depletion region in the cell.
For optimal performance of betavoltaic cells, the most important issue is that the depletion region of the energy converter should match the main region of the ionizing radiation well. According to the Shockley equation, for an ideal p-n homojunction, the built-in potential barrier Vbi at zero bias voltage is given by [1, 5]
where NA is the acceptor doping concentration in the p-type region, ND is the donor doping concentration in the p-type region, and ni is the intrinsic carrier concentration. The built-in potential barrier is regarded as the theoretical maximum open-circuit voltage in nuclear battery systems [18]. The width of the depletion region at zero bias voltage is given by [1, 5]
where ε0 is the vacuum dielectric constant and εs is the relative dielectric constant (which is 9.7 for 4H-SiC material) [1]. In the 4H-SiC-based junction, when both NA and ND vary between 1×1013/cm3 and 1×1018/cm3, the built-in potential barrier increases from 2.5 to 3.1 V and the width of the depletion region decreases from 23.2 to 0.08 μm at 300 K. The results are shown in Figs. 14 and 15, respectively.
-201811/1001-8042-29-11-017/alternativeImage/1001-8042-29-11-017-F014.jpg)
-201811/1001-8042-29-11-017/alternativeImage/1001-8042-29-11-017-F015.jpg)
The figures show that the lower doping concentration produces a lower built-in potential barrier and a larger width of the depletion region. In addition, the width of the depletion region is more sensitive than the built-in potential barrier to the doping concentrations. Next, the width of the depletion region, which extends into the p-type region (Xp) and the n-doped region (Xn) from the metallurgical junction, is given by [22]
and
respectively. Moreover, the relationship between Xp and Xn is expressed as Xp/Xn = ND/NA. Therefore, the design of doping concentrations is the primary method to obtain the appropriate energy converter in betavoltaic cells.
In this research, the depletion region of the energy converter should be ≥8 μm and it is better that it extends into the p-type region. From Figs. 14 and 15, we can obtain that when the light doping concentration of the p-type region is 4.7×1013/cm3, the doping concentration of the n-doped region varies between 5×1014/cm3 and 1×1018/cm3, and the width of depletion region is ∼8 μm. The result is plotted in Fig. 16. The corresponding built-in potential barrier increases from 2.64 to 2.84 V at 300 K, and the values of built-in potential barrier are larger than the theoretical maximum open-circuit voltage in 4H-SiC-63Ni betavoltaic cells. Therefore, the design of doping concentrations mentioned above is reasonable. Furthermore, in Fig. 16, it is shown that, if the light doping concentration of the p-type region is 4.7×1013/cm3, the doping concentration of the n-doped region varies between 3×1016/cm3 and 1×1018/cm3, and the depletion region almost extends into the p-type region.
-201811/1001-8042-29-11-017/alternativeImage/1001-8042-29-11-017-F016.jpg)
Finally, a feasible design of a p-n-junction-based 4H-SiC-63Ni cell is presented. Based on the typical cell structure illustrated in Fig. 6, when the thickness of the 63Ni layer is 2 μm, the thickness of the p-type region is nearly 8 μm, the doping concentration of the p-type region is 4.7×1013/cm3, and the doping concentration of the n-doped region varies between 3×1016/cm3 and 1×1018/cm3, the theoretical calculated cell performance nearly reaches its maximum. The results are listed in Table 2. In addition, the whole thickness of the energy conversion unit should be >23 μm to prevent beta particles from escaping the device for 63Ni layers with different thicknesses.
Parameter | Value |
---|---|
Jsc (μA/cm2) | 0.152 |
Voc (V) | 2.53 |
Pmax (μW/cm2) | 0.36 |
FF (%) | 94.3 |
ηdev (%) | 23.5 |
ηtotal (%) | 11.8 |
5 Conclusion
In summary, 4H-SiC p-n-junction-based betavoltaic batteries with a 63Ni source are investigated. First, the main factors that affect the energy conversion efficiencies of the cells are analyzed. Next, based on the simulation results, theoretical maximum electrical properties and energy conversion efficiencies of the cells are calculated. When the thickness of the 63Ni source increases from 2×10-3 to 10 μm, the theoretical maximum device conversion efficiency increases from 16.77% to 23.51% and the total conversion efficiency decreases from 16.73% to 1.48%. Furthermore, a feasible design with a maximum output power density of 0.36 μW/cm2 and an optimal device conversion efficiency of 23.5% is presented. Finally, the results in our work are helpful for realizing an optimal design in the practical fabrication of 4H-SiC-63Ni cells. Moreover, the simulated and calculated model can be used in p-n-junction-based batteries with other semiconductor materials and beta sources.
A review of nuclear batteries
. Pro. Nucl. Energ. 75, 117-148 (2014). doi: 10.1016/j.pnucene.2014.04.007.Gallium Nitride Schottky betavoltaic nuclear batteries
. Energ. Conv. Manage. 52, 1955-1958 (2011). doi: 10.1016/j.enconman.2010.10.048.63Ni schottky barrier nuclear battery of 4H-SiC
. Radioanal. Nucl. Chem, 287, 173-176 (2011). doi: 10.1007/s10967-010-0746-7.Optimization design and analysis of Si-63Ni betavoltaic battery
. Sci. China. Tech. Sci. 55, 990-996 (2012). doi: 10.1007/s11431-012-4752-6.Optimization design of GaN betavoltaic microbattery
. Sci. China. Tech. Sci. 55, 659-664 (2012). doi: 10.1007/s11431-011-4739-8.GaN PIN betavoltaic nuclear batteries
. Sci. China. Tech. Sci. 57, 25-28 (2014). doi: 10.1007/s11431-013-5422-z.Optimization and temperature effects on sandwich betavoltaic microbattery
. Sci. China. Tech. Sci. 57, 14-18 (2014). doi: 10.1007/s11431-013-5413-0.Development of nuclear microbattery prototype based on Schottky barrier diamond diodes
. Phys. Status. Solidi A. 11, 2539-2547 (2015). doi: 10.1002/pssa.201532214.Influences of planar source thickness on betavoltaics with different semiconductors
. Radioanal. Nucl. Chem. 304, 517-525 (2015). doi: 10.1007/s10967-014-3879-2.Prediction of 4H-SiC betavoltaic microbattery characteristics based on practical Ni-63 sources
. Appl. Radiat. Isotopes. 107, 272-277 (2016). doi: 10.1016/j.apradiso.2015.11.001.Comparative study of different metals for Schottky barrier diamond betavoltaic power converter by EBIC technique
. Phys. Status. Solidi A. 9, 2492-2497 (2016). doi: 10.1002/pssa.201533060.Optimization of Silicon parameters as a betavoltaic battery: Comparison of Si p-n and Ni/Si Schottky barrier
. Radiat. Phys. Chem. 125, 205-212 (2016). doi: 10.1016/j.radphyschem.2016.04.012.Single crystal CVD diamond membranes for betavoltaic cells
. Appl. Phys. Lett. 108, 252105:1-4 (2016). doi: 10.1063/1.4954013.Quantitative modeling of betavoltaic microbattery performance
. Sens. Actuators A. 240, 131-137 (2016). doi: 10.1016/j.sna.2016.01.028.Design and simulation of a semispherical semiconductor to construct a beta-voltaic battery using c-Si and a-Si: H materials with different doping concentration
. J. Comput. Electron. 15, 1577-1592 (2016). doi: 10.1007/s10825-016-0896-0.High efficiency 4H-SiC betavoltaic power sources using tritium radioisotopes
. Appl. Phys. Lett. 108, 013505:1-4 (2016). doi: 10.1063/1.4939203.Temperature dependence of an AlInP 63Ni betavoltaic cell
. J. Appl. Phys. 120, 144501:1-5 (2016). doi: 10.1063/1.4964504.Development of betavoltaic cell technology production based on microchannel silicon and its electrical parameters evaluation
. Appl. Radiat. Isotopes. 121, 71-75 (2017). doi: 10.1016/j.apradiso.2016.12.019.Dissertation
.Betavoltaic study of a GaN p-i-n structure grown by metal-organic vapour phase epitaxy with a Ni-63 source
. Thin Solid Films. 636, 746-750 (2017). doi: 10.1016/j.tsf.2017.07.033Vacuum degree effects on betavoltaics irradiated by 63Ni with differently apparent activity densities
. Sci. China. Tech. Sci. 60, 282-288 (2017). doi: 10.1007/s11431-016-0505-x.