Si24: An Efficient Solar Cell Material - ACS Publications - American

Jul 7, 2017 - Centre for Advanced 2D Materials, National University of Singapore, 6 Science Drive 2, Singapore 117546, Singapore. •S Supporting ...
0 downloads 0 Views 2MB Size
Subscriber access provided by University of South Florida

Article 24

Si : An Efficient Solar Cell Material Jiajun Linghu, Lei Shen, Ming Yang, Shuyan Xu, and Yuan Ping Feng J. Phys. Chem. C, Just Accepted Manuscript • DOI: 10.1021/acs.jpcc.7b04032 • Publication Date (Web): 07 Jul 2017 Downloaded from http://pubs.acs.org on July 10, 2017

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

The Journal of Physical Chemistry C is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 23

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Figure 1. (a) The unit cell of Si24 with the three nonequivalent atomic positions labeled as Site 1~3. (b) PDOS of Si atoms at the three nonequivalent atomic positions. (c) Band structure of Si24 calculated with the hybrid functional. The band gap value is given in the unit of eV. 73x133mm (300 x 300 DPI)

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 2 Band structures (left panel) and PDOS (right panel) of (a) pure Si24 and (b) B-doped Si24, (c) Aldoped Si24, (d) Ga-doped Si24, (e) P-doped Si24, and (f) As-doped Si24. The Fermi level is indicated by the dashed pink line. The band gap value is given in the unit of eV. 136x457mm (300 x 300 DPI)

ACS Paragon Plus Environment

Page 2 of 23

Page 3 of 23

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Figure 3. Partial charge from the VBM or the CBM of B-doped Si24(a), and P-doped Si24(b). The Si atoms are shown in blue balls, and the dopant is shown in pink ball. 14x5mm (300 x 300 DPI)

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 4 of 23

Si24: An Efficient Solar Cell Material Jiajun Linghu, † Lei Shen, ‡ Ming Yang, ¶ Shuyan Xu, § and Yuan Ping Feng †,∥,* †

Department of Physics, 2 Science Drive 3, National University of Singapore, 117551, Singapore



Department of Mechanical Engineering, 9 Engineering Drive 1, National University of

Singapore, 117575, Singapore ¶

Institute of Materials Research and Engineering, A*-STAR, 2 Fusionopolis Way, 138634,

Singapore §

National Institute of Education, Nanyang Technological University, 637616, Singapore



Centre for Advanced 2D Materials, 6 Science Drive 2, National University of Singapore,

117546, Singapore

ACS Paragon Plus Environment

1

Page 5 of 23

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

ABSTRACT

Si24, a recently synthesized allotrope of silicon, has received much attention due to its quasi-direct band gap around 1.3 eV. To explore its potential application in solar cell device, we investigated the doping effect on the electronic properties of Si24 using first-principles calculations. It is found that Si24 can be easily doped as both p- and n-type, respectively, and the dopants are readily ionized. Furthermore, the incorporation of these dopants only reduces the band gap of Si24 slightly, which remains in the ideal region for solar cells. Boron and phosphorus are identified as the most promising elements for the p-type and n-type doping in Si24, respectively, due to their low formation energies, small ionization energies, and small reductions in the band gap. These properties suggest great potential in constructing novel Si24 based p-n junction which is highly desired for future industrial application of photovoltaic devices.

ACS Paragon Plus Environment

2

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 6 of 23

INTRODUCTION Silicon-based photovoltaic devices have been dominating the solar cell industry for several decades due to their elementary and technical advantages such as abundance, high stability, environmental-friendliness, and the relatively high photoelectric transfer efficiency. Even though the widely used crystal silicon in the diamond structure (d-Si) suffers from a number of shortcomings and other semiconductors such as GaAs show better efficiency, the economic advantages of silicon will keep it as the dominant material for the solar cell industry in the near future. The challenge is to improve the efficiency of silicon-based devices. D-Silicon is actually not an ideal material for photovoltaic applications due to two major factors. First, d-Si has low absorption efficiency in visible light region because its band gap (1.12 eV)1 is smaller than the optimal gap (1.4 eV) for achieving the maximum absorption efficiency for the AM1.5 solar spectrum.2 Second, d-Si is an indirect band gap semiconductor, in which its valence band maximum (VBM) and its conduction band minimum (CBM) appear at different crystal momentum (k-point), hence, its light absorption needs the assistance of phonon, leading to a much lower light absorption efficiency compared to direct band gap semiconductors whose VBM and CBM appears at the same k-point. Therefore, to improve the efficiency of siliconbased solar cell devices, it is desirable to have silicon allotropes with larger and direct energy gap, which has been the goal of numerous efforts in engineering silicon structures.3-7 Driven by this desire, various phases of silicon have been synthesized experimentally,8-10 among which the BC8 and R8 phases are found stable at the ambient condition.11-13 However, they are not suitable for solar cell applications because the BC8 phase is semimetallic,14 and the band gap of the R8 phase is too small (0.24 eV).15 Along with the experimental efforts, extensive

ACS Paragon Plus Environment

3

Page 7 of 23

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

computational studies have been carried out to design new phases of Si with desired properties. For example, Botti et al.16 performed a structural search using the minima hopping method and predicted several metastable silicon phases, including M-10, C2221, imma, cmcm, P-1, and P21/c structure, with band gaps ranging from 0.8 to 1.5 eV. Xiang and his colleagues17 proposed a cubic Si20 phase with a quasi-direct gap of 1.55 eV. Similarly, using particle swarm optimization algorithm, Wang et al.18 predicted six novel silicon structures, with band gaps varying between 0.81 and 1.25 eV. Lee et al.19 and Amsler et al.20 also generated a number of silicon phases with direct or quasi-direct band gaps in the range of 0.29 to 1.8 eV. Most of these theoretical predictions, however, have not been verified experimentally. More recently, Kim et al.21 synthesized a silicon phase with the cmcm symmetry. The 24atom unit cell of Si24 can be obtained by removing sodium atoms from the Na4Si24 precursor. This Si24 structure is found to have an indirect band gap of 1.29 eV and a slightly large direct band gap of 1.39 eV, which are close to the optimal value (1.4 eV) for solar cell applications. In addition, its light absorption efficiency is estimated to be higher than that of d-Si in the visible light range. These suggest that Si24 is a promising material for solar cell applications. Obviously, further studies are necessary to evaluate Si24 for solar energy application. For example, an important issue for solar cell application is whether the material can be easily doped to p- and n-type and form p-n junction, and how the doping affects its electronic structure. In the present study, we perform first-principles calculations to investigate the doping properties of Si24. Our calculations predict that both p- and n-type Si24 can be easily achieved by doping with group III and group V elements, respectively, particularly boron and phosphorous. Furthermore, Si24 is found to have similar doping properties as d-Si. Based on these properties, we can expect that p-n junction of Si24 can be easily fabricated and Si24 can be an efficient solar cell material.

ACS Paragon Plus Environment

4

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 8 of 23

COMPUTATIONAL METHODS All calculations are carried out using the density functional theory (DFT) as implemented in the Vienna ab initio simulation package (VASP), with the projector-augmented-wave (PAW)22 method for the electron-ion interaction and the Perdew-Burke-Ernzerhof (PBE)23 exchangecorrelation functional. The cutoff energy for the plane-wave expansion of electron wave function is set at 400 eV, and the force on each atom is relaxed until it is smaller than 0.01 eV/Å in all calculations. It is well known that standard DFT exchange-correlation (XC) functionals within local density approximation (LDA) or generalized gradient approximation (GGA) seriously underestimate band gaps of semiconductors by as much as 50% for some materials, due to the spurious electron self-interaction and lack of the derivative discontinuities of the XC functional24. Various approaches such as GGA+U,25 hybrid functional,26-28 GW,29 etc. have been developed to address this issue. The hybrid functional method, which is performed by mixing a fraction of nonlocal exact exchange energy with PBE exchange, is accepted as a practical approach due to its accuracy and relatively low demand for computational resources and is widely used to calculate the electronic structures of semiconductors. In this work, the HSE06 hybrid functional30-31, with 34% exact exchange energy, is applied in band structure calculation to obtain a more reliable band gap. Doping is modeled using a 3×1×1 supercell containing 72 atoms which is sufficient to minimize the potential Coulomb interaction among dopants. The Monkhorst-Pack scheme with 9×3×3 and 3×3×3 k-point meshes are used for the unit cell and supercell calculations, respectively. The formation energy of a doping defect X in charge state q is given by:32  (X  ) = tot (X  ) − tot (bulk) − ∑   + ( +  + ∆)

(1)

ACS Paragon Plus Environment

5

Page 9 of 23

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

where tot (X  ) is the total energy of the supercell with the defect, X, tot (bulk) is the total energy of bulk (undoped) Si24 supercell.  is the number of atoms of type  (host atoms or dopant) that have been added to ( >0) or removed from ( < 0) the supercell to construct the defect, and  are the corresponding chemical potential. In this work, Si24 and pure metal of aluminum, gallium, and arsenic are used as the source of the Si, Al, Ga, and As element to obtain the chemical potential, respectively. For the boron element, it has many allotropes such as αrhombohedral, α-tetragonal, β-rhombohedral, β-tetragonal, and γ-boron phase, among which the β-rhombohedral phase is most stable and selected as the origin of boron in our calculation. In addition, phosphorus owns four common allotropes, white, red, violet, and black, among which the least reactive and thermodynamically stable black phosphorous is used as the source of phosphorus element.  is the Fermi energy that is varied within the band gap with respect to the VBM ( = 0 at the VBM), hence, is affected by the reference of the VBM. However, the existence of defect would give rise to a constant shift in the electrostatic potential of VBM between defect supercell and the bulk. To align the VBM, a method was proposed by Van de Walle and his colleagues that the energy of VBM, i.e. the term  in the formula, is firstly calculated in bulk, and subsequently, the potential in the imperfect supercell is measured at where far from the defect and compared with that in the bulk, achieving the potential shift ∆. Therefore, the VBM in defect supercell can be aligned to that in bulk through adding ∆ to  .32-33 Based on the defect formation energy, the thermodynamic transition level ( /! ) is defined as the Fermi level position where the formation energies of a defect with charge states  and ! are equal:32

( /! ) =

"# ($ %& ;"( )*) + "# ($ %, ;"( )*) , +&

(2)

ACS Paragon Plus Environment

6

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 10 of 23

One thing worth to note is that, in the formation energy calculation of charged defect, the Fermi level changes from the VBM to the CBM in the system. Meanwhile, the position of transition level is decided by the energy differences between different charge states. Hence, the accuracy of the transition level is also affected by the band gap problem caused by general DFT functional. Therefore, the HSE06 method with 34% exact exchange energy is also used in all the calculations related to the formation energy of charged defect, as well as the transition levels.

RESULTS AND DISCUSSION The atomistic structure of Si24 is shown in Figure 1a. The unit cell of the orthorhombic structure contains 24 Si atoms occupying three non-equivalent sites (labeled as Site-1~3). The calculated projected density of states (PDOS), shown in Figure 1b, indicate that the Si atoms at these different sites contribute almost equally to the total density of states (DOS). As seen clearly in the calculated HSE06 band structure shown in Figure 1c, Si24 is a semiconductor with an indirect band gap of 1.29 eV which is in excellent agreement with the measured value.21 The VBM is located at the Γ point but the CBM is near the mid-point along the Γ-X line (Details of the reciprocal space and high symmetry k-points are given in Figure S1 in the supporting information). It is also noted that the direct band gap at Γ is 1.35 eV, only 0.06 eV larger than the indirect band gap. This suggests simultaneous optical transitions corresponding to the indirect and the direct band gap, therefore a high solar adsorption efficiency, and potential for further band structure engineering of Si24-based devices. In contrast, d-Si is an indirect band gap semiconductor, but the energy difference between its indirect and direct band gap is as large as 2.28 eV.34

ACS Paragon Plus Environment

7

Page 11 of 23

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Figure 1. (a) The unit cell of Si24 with the three nonequivalent atomic positions labeled as Site 1~3. (b) PDOS of Si atoms at the three nonequivalent atomic positions. (c) Band structure of Si24 calculated with the hybrid functional. The band gap value is given in the unit of eV. The basic component of a solar cell device is the p-n junction. To produce a p-n junction, it is essential to dope the semiconductor to p-type and n-type which are traditionally done by substituting Si with group III and group V elements, respectively. Here, we carry out first principles calculations to investigate the feasibility of realizing p-type and n-type Si24. Group III elements B, Al, and Ga are considered for p-type doping, while group V elements P and As are investigated for n-type doping. Doped Si24 is modeled by introducing one dopant into the 3×1×1 Si24 supercell, which corresponds to a doping concentration of 1.39%. Both the interstitial and

ACS Paragon Plus Environment

8

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 12 of 23

substitutional doping are considered. However, the calculated defect formation energies (Table 1) indicate that all dopants prefer substituting a Si atom rather than the interstitial sites. The defect formation energy for substitutional doping ranges from -0.03 eV to 0.81 eV which are much lower compared to those at the interstitial sites (9.03~11.40 eV). Due to the high formation energies, doping at the interstitial sites can essentially be ruled out. In the following, we focus our discussion on substitution doping only. Furthermore, a dopant can substitute a Si atom at any of the three inequivalent sites. The calculated formation energies differ by less than 0.1 eV and their electronic structures are also very similar. In Table 1, we list the lowest formation energy in each case. It is also noted that the defect formation energies of various dopants considered are quite low and the highest value is 0.72 eV, corresponding to substitution of Al for Si. At the same time, the intrinsic defects have much higher formation energies (Table S1 in the supporting information). This suggests that both p- and n-type doping can be easily realized in Si24 using these elements as the dopants. Results of our calculations also indicate that it is relatively easy to achieve n-type doping in Si24. In particular, the negative formation energy for PSi suggests that this doping process is an exothermic reaction,35 and can occur spontaneously.36 This could be due to the similar atomic radius of Si and P and the small structural changes induced by the substitution of P for Si atom. The P-Si bond in P-doped Si24 differs by 0.5% compared to the SiSi bond in the undoped Si24. For other dopants, the bond length varies by 2.5% to 6%.

ACS Paragon Plus Environment

9

Page 13 of 23

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Table 1. Formation energy ( E f ) for substitutional and interstitial defects, calculated ionization energy (ε ) of various dopants in Si24. The calculated ( ε d − Si Cal.) and experimental ionization energy ( ε d − Si Exp.) of the same defect in diamond silicon are included for comparison. B

Al

Ga

P

As

E f (substitutional) (eV)

0.60

0.72

0.49

-0.03

0.24

E f (interstitial) (eV)

11.40

9.56

9.03

9.99

10.04

0.08~0.12

0.10~0.14

0.09~0.13

0.09~0.13

0.11~0.15

ε d − Si Cal. (eV)

0.06

0.09

0.10

0.06

0.08

ε d − Si Exp. (eV)a

0.045

0.067

0.072

0.045

0.054

ε

a

(eV)

Reference37

For n- and p-type doping, another crucial parameter is the ionization energy which should be as low as possible so that the excess charges introduced by the dopants can break away from them and become conducting charges.37 For the dopants considered here, the ionization energy is given by the thermodynamic transition level between the neutral (0) and the +1 (for donors) or -1 (for acceptors) charged states of the defect. The thermodynamic transition level also determines the electronic behavior and is often used as the basis for experimental detection or identification of the defect.33 However, it was reported that transition level derived from defect formation energies is reliable for deep transition levels,38 but questionable for shallow levels.39-40 Here, we benchmark our calculations for Si24 against the experimental ionization energy and theoretical value, calculated using the same method, of d-Si whose doping properties are well known. As can be seen in Table 1, the calculated ionization energies of d-Si are close but generally larger than the experimental values. More importantly, the calculated ionization energies show the

ACS Paragon Plus Environment

10

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 14 of 23

correct chemical trend. By comparing the calculated ionization energy for Si24 with the value of d-Si, dopant ionization in Si24 can be assessed relative to that in d-Si. A couple of points are worth noting in the calculations. First, the dipole correction on the total energy of the system with charged defect should be scaled by the dielectric constant of the system. In the absence of the experimental dielectric constant for Si24, we used the calculated value which is usually larger than the experimental one (the difference is about 9% for d-Si, see Table S2 in supporting information). Therefore, the calculation would provide an upper bound for the ionization energy. Second, due to the anisotropic structure of Si24, the dielectric tensor has three different diagonal elements which lead to different values of the calculated ionization energy. In Table 1, the range of ionization energy calculated using these dielectric constants is listed for each dopant. On the whole, the ionization energies of various dopants in Si24 are slightly higher than but close to those in d-Si, indicating similar p- and n-type conductivity in doped Si24 and d-Si. Our calculations also show that B and P are the most favorable dopant to produce p- and n-type conductivity, respectively. Meanwhile, the ionization energy of Ga is slightly lower than that of Al, in contrast to the trend in d-Si. The effective doping of Si24 by the elements being studied is confirmed by the calculated electronic structures and PDOS, as shown in Figure 2. It is clear that doping with group III element, B, Al, and Ga, introduces an acceptor level which mixes with the valance band of Si24 and pins the Fermi level at the valence band, resulting in p-type semiconductors. On the contrary, doping with group V element, P and As, produces a donor level which is combined with the conduction band, resulting in n-type semiconductors. Orbital contributions to the states near the Fermi level can be seen in the calculated PDOS, where for p-type doping, the VBM mainly consists of p-orbitals of the Si atom and the dopant, while for n-type doping, the CBM is a

ACS Paragon Plus Environment

11

Page 15 of 23

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

hybridization of s- and p-orbitals of Si and the dopant. Compared to other dopants, the p-orbital of Ga has a higher intensity near the VBM, indicating the dominant role of the Ga-p-orbital near the Fermi level.

Figure 2 Band structures (left panel) and PDOS (right panel) of (a) pure Si24 and (b) B-doped Si24, (c) Aldoped Si24, (d) Ga-doped Si24, (e) P-doped Si24, and (f) As-doped Si24. The Fermi level is indicated by the dashed pink line. The band gap value is given in the unit of eV.

ACS Paragon Plus Environment

12

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 16 of 23

It is worth noting that doping with these elements does not introduce any obvious defect levels in the band gap of Si24. Due to the low ionization energy, the donor and acceptor levels are mixed with the CBM or VBM, respectively, resulting in a small reduction of the band gap. The amount of the band gap narrowing is about 7%, 24%, 12%, 8% and 0.8% for B, Al, Ga, P, and As doping, respectively. The shapes of the CBM and the VBM are preserved and not much affected by the doping. In addition, the extra electron and hole introduced by the dopants are delocalized, as shown in Figure 3 for P- and B-doped Si24. Similar charge delocalization is found in other n-type and p-type Si24 and the details are given in Figure S2 in the supporting information. The delocalized charge distributes along the close-packed atomic chains in the [100] direction (x-axis). Thus, good conductivity can be expected along the [100] direction for both n-type and p-type Si24.

Figure 3. Partial charge from the VBM or the CBM of B-doped Si24(a), and P-doped Si24(b). The Si atoms are shown in blue balls, and the dopant is shown in pink ball.

ACS Paragon Plus Environment

13

Page 17 of 23

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

CONCLUSION In conclusion, based on detailed first-principles calculations, we found that the formation energies of group III and group V element doped Si24 are very low, p- and n-type Si24 and their junction can be easily produced. The dopants have low ionization energies and therefore can be easily ionized. The donor level below the conduction band minimum is mixed with the conduction band and the acceptor level above the valence band maximum hybridizes with the valence band, leading to a small reduction of the band gap, without affecting the shapes of the CBM and VBM. These doping properties are very similar to that of d-Si. On the other hand, even though the band gap is reduced slightly by doping, the band gaps of B-, Ga-, P- and As-doped Si24 remain in the ideal range (1.1 ~ 1.7 eV) for solar cell applications. In addition, the existence of a direct energy gap with a gap value close to that of the indirect gap effectively makes Si24 a quasi-direct energy gap material which will lead to much more efficient light absorption, compared to d-Si. Even though overall efficiency of a solar cell device depends on other dynamic factors such as rate of electron and hole recombination, one can expect that the combined effect of the optimal band gap and the efficient light absorption will lead to higher efficiency of Si24 based devices than the current solar cell devices based on d-Si. For practical use, we suggest P and B as the best dopant for n-type and p-type doping, respectively, because of their relatively low formation energy, low ionization energy, and band gap closer to the ideal value.

ACS Paragon Plus Environment

14

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 18 of 23

ASSOCIATED CONTENT Supporting Information The reciprocal space with high-symmetry k-point path of Si24. Formation energy of intrinsic defects in Si24. Dielectric constant of d-Si and Si24. (PDF)

AUTHOR INFORMATION Corresponding Author *E-mail: [email protected]

ACKNOWLEDGMENT The authors are indebted to Prof. Fan Haiming, Dr. Zhou Jun, and Mr. Yang Tong for helpful discussions. The computational resources were provided by Centre for Advanced 2D Materials, National University of Singapore

ACS Paragon Plus Environment

15

Page 19 of 23

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

REFERENCES (1) Bludau, W.; Onton, A.; Heinke, W. Temperature Dependence of the Band Gap of Silicon. Journal of Applied Physics 1974, 45, 1846-1848. (2) Shockley, W.; Queisser, H. J. Detailed Balance Limit of Efficiency of P‐N Junction Solar Cells. Journal of Applied Physics 1961, 32, 510-519. (3) Yoshitaka, F.; Takashi, K.; Susumu, S.; Takashi, M.; Atsushi, O. A New Crystalline Phase of Four-Fold Coordinated Silicon and Germanium. New Journal of Physics 2008, 10, 083001. (4) Malone, B. D.; Cohen, M. L. Prediction of a Metastable Phase of Silicon in the Ibam Structure. Physical Review B 2012, 85, 024116. (5) Jinhui, Z.; Dongli, Y.; Kun, L.; Qianqian, W.; Zhisheng, Z.; Julong, H.; Yongjun, T. Metastable C-Centered Orthorhombic Si 8 and Ge 8. Journal of Physics: Condensed Matter 2012, 24, 405803. (6) Cohen, M. L.; Malone, B. D. Wave Function Engineering: Other Phases of Si for Photovoltaic Applications. Journal of Applied Physics 2011, 109, 102402. (7) Needs, R. J.; Mujica, A. First-Principles Pseudopotential Study of the Structural Phases of Silicon. Physical Review B 1995, 51, 9652-9660. (8) Duclos, S. J.; Vohra, Y. K.; Ruoff, A. L. Experimental Study of the Crystal Stability and Equation of State of Si to 248 Gpa. Physical Review B 1990, 41, 12021-12028. (9) McMahon, M. I.; Nelmes, R. J. New High-Pressure Phase of Si. Physical Review B 1993, 47, 8337-8340. (10) McMahon, M. I.; Nelmes, R. J.; Wright, N. G.; Allan, D. R. Pressure Dependence of the Imma Phase of Silicon. Physical Review B 1994, 50, 739-743. (11) Minomura, S.; Drickamer, H. G. Pressure Induced Phase Transitions in Silicon, Germanium and Some Ш-V Compounds. Journal of Physics and Chemistry of Solids 1962, 23, 451-456. (12) Crain, J.; Ackland, G. J.; Maclean, J. R.; Piltz, R. O.; Hatton, P. D.; Pawley, G. S. Reversible Pressure-Induced Structural Transitions between Metastable Phases of Silicon. Physical Review B 1994, 50, 13043-13046. (13) Piltz, R. O.; Maclean, J. R.; Clark, S. J.; Ackland, G. J.; Hatton, P. D.; Crain, J. Structure and Properties of Silicon Xii: A Complex Tetrahedrally Bonded Phase. Physical Review B 1995, 52, 4072-4085. (14) Besson, J. M.; Mokhtari, E. H.; Gonzalez, J.; Weill, G. Electrical Properties of Semimetallic Silicon Ш and Semiconductive Silicon Iv at Ambient Pressure. Physical Review Letters 1987, 59, 473-476. (15) Malone, B. D.; Sau, J. D.; Cohen, M. L. Ab Initio Study of the Optical Properties of SiXii. Physical Review B 2008, 78, 161202. (16) Botti, S.; Flores-Livas, J. A.; Amsler, M.; Goedecker, S.; Marques, M. A. L. Low-Energy Silicon Allotropes with Strong Absorption in the Visible for Photovoltaic Applications. Physical Review B 2012, 86, 121204. (17) Xiang, H. J.; Huang, B.; Kan, E.; Wei, S.-H.; Gong, X. G. Towards Direct-Gap Silicon Phases by the Inverse Band Structure Design Approach. Physical Review Letters 2013, 110, 118702.

ACS Paragon Plus Environment

16

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 20 of 23

(18) Wang, Q.; Xu, B.; Sun, J.; Liu, H.; Zhao, Z.; Yu, D.; Fan, C.; He, J. Direct Band Gap Silicon Allotropes. Journal of the American Chemical Society 2014, 136, 9826-9829. (19) Lee, I.-H.; Lee, J.; Oh, Y. J.; Kim, S.; Chang, K. J. Computational Search for Direct Band Gap Silicon Crystals. Physical Review B 2014, 90, 115209. (20) Amsler, M.; Botti, S.; Marques, M. A. L.; Lenosky, T. J.; Goedecker, S. Low-Density Silicon Allotropes for Photovoltaic Applications. Physical Review B 2015, 92. (21) Kim, D. Y.; Stefanoski, S.; Kurakevych, O. O.; Strobel, T. A. Synthesis of an OpenFramework Allotrope of Silicon. Nature Materials 2015, 14, 169-173. (22) Blochl, P. E. Projector Augmented-Wave Method. Physical Review B 1994, 50, 1795317979. (23) Perdew, J. P.; Burke, K.; Ernzerhof, M. Generalized Gradient Approximation Made Simple. Physical Review Letters 1996, 77, 3865-3868. (24) Perdew, J. P. Density Functional Theory and the Band Gap Problem. International Journal of Quantum Chemistry 1985, 28, 497-523. (25) Anisimov, V. I.; Solovyev, I. V.; Korotin, M. A.; Czyzyk, M. T.; Sawatzky, G. A. Density-Functional Theory and Nio Photoemission Spectra. Physical Review B 1993, 48, 1692916934. (26) Seidl, A.; Gorling, A.; Vogl, P.; Majewski, J. A.; Levy, M. Generalized Kohn-Sham Schemes and the Band-Gap Problem. Physical Review B 1996, 53, 3764-3774. (27) Perry, J. K.; Tahir-Kheli, J.; Goddard, W. A. Antiferromagnetic Band Structure of La2cuo4: Becke-3-Lee-Yang-Parr Calculations. Physical Review B 2001, 63. (28) Muscat, J.; Wander, A.; Harrison, N. M. On the Prediction of Band Gaps from Hybrid Functional Theory. Chemical Physics Letters 2001, 342, 397-401. (29) Aryasetiawan, F.; Gunnarsson, O. The Gw Method. Reports on Progress in Physics 1998, 61, 237. (30) Heyd, J.; Scuseria, G. E. Efficient Hybrid Density Functional Calculations in Solids: Assessment of the Heyd-Scuseria-Ernzerhof Screened Coulomb Hybrid Functional. Journal of Chemical Physics 2004, 121, 1187-1192. (31) Krukau, A. V.; Vydrov, O. A.; Izmaylov, A. F.; Scuseria, G. E. Influence of the Exchange Screening Parameter on the Performance of Screened Hybrid Functionals. The Journal of Chemical Physics 2006, 125, 224106. (32) Van de Walle, C. G.; Neugebauer, J. First-Principles Calculations for Defects and Impurities: Applications to Ш-Nitrides. Journal of Applied Physics 2004, 95, 3851-3879. (33) Freysoldt, C.; Grabowski, B.; Hickel, T.; Neugebauer, J.; Kresse, G.; Janotti, A.; Van de Walle, C. G. First-Principles Calculations for Point Defects in Solids. Reviews of Modern Physics 2014, 86, 253-305. (34) Hybertsen, M. S.; Louie, S. G. First-Principles Theory of Quasiparticles: Calculation of Band Gaps in Semiconductors and Insulators. Physical Review Letters 1985, 55, 1418-1421. (35) dos Santos, L. S.; Schmidt, W. G.; Rauls, E. Group-Vii Point Defects in Znse. Physical Review B 2011, 84, 115201. (36) Yamauchi, J.; Aoki, N. Deactivation Mechanism of in Atoms Doped in a Si Crystal and Reactivation Due to Codoping with B and C. Physical Review B 2005, 71. (37) Kittel, C. Introduction to Solid State Physics. John Wiley & Sons, New York, 2005, 209212. (38) Schultz, P. A. Theory of Defect Levels and the ``Band Gap Problem'' in Silicon. Physical Review Letters 2006, 96, 246401.

ACS Paragon Plus Environment

17

Page 21 of 23

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

(39) Zhang, G.; Canning, A.; Gronbech-Jensen, N.; Derenzo, S.; Wang, L.-W. Shallow Impurity Level Calculations in Semiconductors Using Ab Initio Methods. Physical Review Letters 2013, 110. (40) Wang, L.-W. Density Functional Calculations of Shallow Acceptor Levels in Si. Journal of Applied Physics 2009, 105.

ACS Paragon Plus Environment

18

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 22 of 23

TOC Graphic

ACS Paragon Plus Environment

19

Page 23 of 23

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

TOC 47x30mm (300 x 300 DPI)

ACS Paragon Plus Environment