Solar Energy Conversion by Nanostructured TiO2View this Special Issue
Research Article | Open Access
Abou El-Maaty Aly, Ashraf Nasr, "Theoretical Study of One-Intermediate Band Quantum Dot Solar Cell", International Journal of Photoenergy, vol. 2014, Article ID 904104, 10 pages, 2014. https://doi.org/10.1155/2014/904104
Theoretical Study of One-Intermediate Band Quantum Dot Solar Cell
The intermediate bands (IBs) between the valence and conduction bands play an important role in solar cells. Because the smaller energy photons than the bandgap energy can be used to promote charge carriers transfer to the conduction band and thereby the total output current increases while maintaining a large open circuit voltage. In this paper, the influence of the new band on the power conversion efficiency for the structure of the quantum dots intermediate band solar cell (QDIBSC) is theoretically investigated and studied. The time-independent Schrödinger equation is used to determine the optimum width and location of the intermediate band. Accordingly, achievement of maximum efficiency by changing the width of quantum dots and barrier distances is studied. Theoretical determination of the power conversion efficiency under the two different ranges of QD width is presented. From the obtained results, the maximum power conversion efficiency is about 70.42% for simple cubic quantum dot crystal under full concentration light. It is strongly dependent on the width of quantum dots and barrier distances.
The intermediate band solar cells (IBSCs) have attracted great attention due to the possibility of exceeding the Shockley-Queisser (SQ) limit [1–4]. From the analysis by Luque et al. [5, 6], the IBSC’s concept yields a maximum theoretical efficiency of 63.2%, surpassing the limit of 40.7% of single gap solar cells under maximum light concentration (the sun being assumed as a blackbody at 6000 K) . Since its introduction in 1997, there have been many theoretical and experimental efforts to explore this idea . The use of quantum dot (QD) technology is proposed as a near term proof of concept of the operating principles of an intermediate band solar cell (IBSC). This intermediate band allows the extra generation of electron-hole pairs through the two-step absorption of subbandgap photons. In the first step, an electron is pumped from the valence band (VB) to the intermediate band (IB), while in the second step, another electron is launched from the IB to the CB . Quantum dot heterojunctions may implement an IBSC because of their ability to provide the three necessary bands . In comparison to conventional quantum well superlattices or multiple quantum well structures, quantum dot superlattice (QDS) that consists of multiple arrays of quantum dots has many advantages due to its modified density of electronic states and optical selection rules. For example, due to relaxed intraband optical selection rules in QDS, they are capable of absorbing normally incident radiation, for example, Indium gallium nitride alloys feature a bandgap ranging from the near infrared (~0.7 eV) to ultraviolet (~3.42 eV); this range corresponds very closely to the solar spectrum, making alloys a promising material for future solar cells. alloys solar cells have been fabricated with different indium contents and the results are encouraging . However, it is not possible in quantum well superlattices. The latter makes QDS a good candidate for infrared photodetector applications . As considered in , 3D-ordered QDS with the closely spaced quantum dots and high quality interfaces allow a strong wave function overlap and the formation of minibands. In such structures, the quantum dots play a similar role to that of atoms in real crystals. To distinguish such nanostructures from the disordered multiple arrays of quantum dots, we refer to them as quantum dot supracrystals (QDC). Other assumptions that include the rules for solving the Schrödinger equation for determining the IBs are considered in the following sections. The remainder of the paper is organized as follows. In Section 2, the basic assumptions and a superlattice model description are presented. The intermediate band energy and wavevectors of a charge carrier are calculated in Section 3. The induced photocurrent density and power conversion energy are described in Section 4. The numerical results and discussions are summarized in Section 5. Finally, a conclusion of the results is outlined in Section 5.
2. Basic Assumptions and Superlattice Model Description
In this paper, the investigations are devoted to one intermediate band solar cell. We explain the various behaviors of this model depending on the QD solar cell parameters, such as quantum dot size, interdot distance, and type of composition alloy. Meanwhile, in the next step, the multi-intermediate bands are investigated. Therefore, the following are the basic assumptions that are used in the QD one intermediate band solar cell model [11–14]. (a)The solar cell model is thick enough; carrier mobility will be high enough to ensure the full absorption of photons. All photons with energy greater than the lowest energy gap in the QDIBS model are absorbed. (b)The quasi-Fermi energy levels which are equivalent to infinite carrier mobility are constant throughout the model. (c)The transitions occurring between the bands are only radiative recombination. (d)The solar cell absorbs blackbody radiation at a sun temperature, 6000 K, and the temperature of the solar cell, 300 K, and emits blackbody radiation at 300 K only. (e)No carriers can be extracted from the intermediate band; the net pumping of electrons from the VB to the IB must equal the net pumping of electrons from the IB to the CB. (f)The shape of the QDs is cubic and they must be arranged in a periodic lattice in order to establish well-placed intermediate band boundaries. For simplicity, the energy corresponding to the top of the valence band is the same both in the barrier and the QD material; therefore there is no valence band offset and only confining potential occurs at the conduction band offset. (g)The intermediate band should be approximately half-filled with electrons in order to receive electrons from the valence energy band and pump electrons to the conduction energy band.
When charge carriers in semiconductors can be confined by potential barriers in three dimensions, it is called quantum dots. QDs periodic arrays from semiconductor with a smaller bandgap are sandwiched between two layers of another second semiconductor having a larger bandgap (n or p type). This configuration creates the potential barriers. Two potential wells are formed in this structure; one is for conduction band electrons and the other for valence band holes. The well (QDs layer) depth for electrons is called the conduction band offset, which is the difference between the conduction band edges of the well and barrier semiconductors. The well depth for holes is called the valence band offset. If the offset for either the conduction or valence bands is zero, then only one carrier will be confined in a well. In this structure, if the barrier thickness between adjacent wells prevents significant electronic coupling between the wells, then each well is electronically isolated. On the other hand, if the barrier thickness is sufficiently thin to allow electronic coupling between wells, then the electronic charge distribution can become delocalized along the direction normal to the well layers, therefore producing new minibands (see Figure 1).
The electronic coupling rapidly increases with decreasing the barrier thickness and miniband formation is very strong below 2 nm . Superlattice structures yield efficient charge transport normal to the layers because the charge carriers can move through the minibands. As a result the barrier will be narrower and the miniband and the carrier mobility will be wider and higher, respectively.
3. Wavevectors of a Charge Carrier
The wavevector of a charge carrier (single electron or hole) can be described by the time-independent Schrödinger equation, which has the following form [11, 16]: where is the Plank’s constant, is the effective mass, is a second order differential operator, is the potential energy, is the total energy of charge carrier, and is the wavevector. The Kronig-Penney model solved this equation by the one-dimensional periodic potential shown also in Figure 1.
This model assumes that the wave travelling of charge carrier is in the positive direction for one-dimensional only. Therefore, the mathematical form of the repeating unit of the potential is Here, is the conduction band offset, is the barrier width, and is the quantum dot width. The period of the considered potential is equal to . The Schrödinger equation for this model is [10, 17, 18].Consider According to the Kronig-Penney model, the solution of (3a) and (3b) can be expressed as For simplicity, one can assume the following symbols for internal terms in (4): where , are effective mass of electron in barrier region and effective mass of electron in quantum dots region, respectively.
Therefore, from (4), the factor into the first term plays an important role for investigating this proposed QDIBSC model. It can be expressed as follows: Furthermore, the other arguments in (4) for hyperbolic and sinusoidal functions can be defined as Substituting these definitions into (4), it will become asThe left-hand side of (8a), (8b), and (8c) can be represented by for all values of the ratio of total energy of electrons over conduction band offset: (). Consider Equation (9) cannot be solved analytically, but it can be solved graphically. Figure 2 shows the left-hand side of (9), , against at fixed values of and . The concerned QDIBSC model depends on the considered alloys in previous studies for comparing the results; the quantum dot width: (InAs0.9N0.1, ), barrier width: nm (GaAs0.98Sb0.02, .066 ), and eV . But in further studies, other alloys will be processed. The left-hand side is not constrained to and is a function of energy only. The right-hand side is constrained to a range of and is a function of only. The limits of the right-hand side occur at to . The two horizontal red lines in Figure 2 represent the two extreme values of . The only allowed values of are those for which lies between the two horizontal red lines; see Table 1.
Figure 2 and Table 1 show the allowed energy ranges or bands, , which fall into continuous regions separated by gaps. This distribution of allowed energies illustrates the band structure of crystalline solids. According to Figure 2 and Table 1, when the values of are from 0 to 0.2039 eV, is greater than one and so there is no real value of that satisfies (9). At .2039 eV the wavevector will be equal to zero. As varies from 0.2039 eV to 0.2369 eV, varies from +1 to −1 and varies from 0 to . As varies from 0.2369 eV to 0.8140 eV, there is no real value of that satisfies (9). At .8140 eV, the argument of the right-hand side, , will be equal to . At eV, the argument of the right-hand side, , will be equal to . Other bandgaps will periodically appear with the same behavior.
Figure 3 shows the allowed energy bands and the wavevector states in one-dimensional crystal with the same data as in Figure 2. In this case, the wavevector is changed from 0 to , and the allowed energy bands into the electron energy is assigned. Also, the curve shows that the slope is zero at the boundaries (i.e., 0 and ). Thus, the velocity of the electrons approaches zero at the boundaries. This means that the electron trajectory or momentum is confined to stay within the allowable . Figure 4 denotes more explanation for the considered energy allowed bands as in Figure 3. It shows the allowed - states with the bandgaps that appear when . The dashed line represents - states for free electrons energy. It can be obtained by letting in (4). Therefore, the free energy of electron is . After the essential equations and assumptions are demonstrated, the following section is concerned with the determination of the effect of intermediate band, alloy construction, and dot and barrier width in each of induced photocurrent density and corresponding power efficiency.
4. Induced Current Density and Power Conversion Efficiency
The photon generated current density in QDIBSC with one intermediate band is derived in this section. Then, the sensitivity of the power conversion efficiency as a function of the intermediate band energy level will be investigated. According to the assumptions in Section 2, the only radiative transitions occur between the bands therefore the generation and recombination events are represented by photon absorption and emission. Figure 5 illustrates the construction of an energy band diagram for a heterostructure in the case of one intermediate band solar cell. An electron in the valence band can be excited to either the intermediate or conduction band. Also, an electron in the intermediate band can be excited to the conduction band. Therefore, there are three upward energy transitions: , , and . represents valence to intermediate band, represents intermediate to conduction band, and represents the conventional bandgap between the valence and conduction band. The two intermediate transitions and are independent of each other, while the bandgap transition is a function of the two intermediate ones: .
The photon flux density, , is the number of photons per second per unit area per unit wavelength and behaves like a blackbody flux density and according to the Roosbroeck-Shockley equation [19–22] is given by where is the temperature, is geometric factor, is Plank constant, is Boltzmann constant, is speed of light, and is chemical potential. To simplify the analysis, assume that . Therefore the photons with energy greater than and less than are absorbed and electrons transfer from valence band to intermediate band and leave holes in the valence band. Any excess energy greater than and less than will be lost due to thermalization and carriers will relax to the band edges before another radiative occurs. One can notice that this thermalization value is very small in comparison with its value in the case of bulk based semiconductor solar cells. This means that an absorbed photon with energy greater than and less than has the same effect as an absorbed photon with energy equal to .
Photons with energy greater than and less than are absorbed and an electron transfers from intermediate band to conduction band. The excess energy behavior has the same effect as considered in the previous case regarding the difference in the transition band process. Photons with energy greater than are absorbed and an electron transfers from valence band to conduction band and creates a hole in the valence band. The absorbed photon with energy greater than has the same effect as an absorbed photon with energy equal to taking into account the excess energy process. The net photon flux is equal to the number of charge carrier flux collected at the contact. When the charge carrier flux is multiplied by the electric charge of electron, , the current density of the QDIBSC for one intermediate band is [11, 13, 14]: where is concentration factor, is temperature of sun (6000 K), is ambient temperature (300 K), is quasi-Fermi energy, and is chemical potential between conduction and intermediate bands. The terms in the first bracket represent the current density generated when the electrons transfer from the valence band to the conduction band as typically for conventional solar cell. While the terms in the second bracket represent the current density generated when the electrons transfer from the intermediate band to the conduction band. In both bracketed terms, the QDIBSC absorbs radiation from the sun at the temperature and , respectively, while it emits radiation at the temperature and a corresponding chemical potential. The current density of the QDIBSC is formulated according to the proper operation of the QDIBSC which requires that there is no current extracted from the intermediate band; that is, the current entering the intermediate band must equal the current leaving the intermediate band. Therefore, the second term in (11) can be rewritten as [14, 23, 24] The output voltage can be described as the difference of the chemical potentials between CB and VB; that is, .
In this work, the light intensity on QDIBSC is calculated by the number of suns, where 1 sun (or concentration factor ) means the standard intensity at the surface of the Earth’s atmosphere. Therefore, at the surface of the Earth’s atmosphere the power density falling on a QDIBSC is w/m2, where is Stefan’s constant and is temperature of sun (6000 K). Theoretically, the full concentration would be achieved when . The power conversion efficiency, , of the QDIBSC is dependent on , so that it varies with the level concentration of . We concentrate our study on the QDIBSC efficiencies with unconcentrated light and also compared it with full concentration light . The power conversion efficiency equation of the QDIBSC is where is open circuit voltage, is short circuit current density, is fill factor, is maximum voltage of the QDIBSC, and is maximum current density of the QDIBSC. After mathematical concepts and assumptions are defined in previous sections for assigned composition of material, the following section reports some features of the QDIBSCs performance for one intermediate band case.
5. Numerical Results and Discussions
The discussion firstly starts by manifesting the relation between QDIBSC parameters and the distribution of IB energies in the gap between valence and conduction bands. When the Schrödinger equation (1) is solved, many solutions are obtained. Some of them can be satisfied but others cannot. From the satisfied solutions, there are multi-intermediate bands as pointed from Figures 2, 3, and 4. They are essentially dependent on QDIBSC parameters such as QD, barrier widths. Here, we are concerned only with the effect of the first intermediate band, IB, into the behavior of the proposed model. From the obtained results, we found two ranges of quantum dot widths that vary the positions of intermediate bands. Other ranges of quantum dot widths will denote unachievable behavior. In the following discussions, the energy distributions, - characteristics, and corresponding power conversion efficiency are studied for each range. Comparisons between the two ranges outcomes are processed. A summary of all the related values describing the model is considered in Table 2. For the first range (QD width changes from 4–7 nm), the dependence of first intermediate band width, , on QD width at different barrier widths is depicted in Figure 6. One can recognize that the width of IB, , is decreased with an increase in the QD width. This behavior is consistent with the previous investigations that demonstrate that the decreasing of IB will assist in the fall down of the power conversion efficiency, as will be described later [11, 14]. Also the same trend will be found in the second range of QD and barrier widths, which also denotes the first IB. Also, Figure 6 illustrates that when the barrier width increases in the allowable range (2-3 nm), IB somewhat increases. At the same time, the increase of barrier width will somewhat inversely decrease the energy gap between IB width and each of VB and CB, as this behavior is shown in Figures 7(a) and 7(b). It is expected that when IB enlarges, the energy gap difference between it and the considered VB and CB will be decreased, as the total distance is constant [11, 14]. The first range power efficiency and - characteristics will be compared with the second range in the last part of these discussions.
For the second allowable range (QD width changes from 8.1–9.1 nm) and at the same barrier width range (2-3 nm), the intermediate band width, , with QD width at different values of barrier width is considered in Figure 8. When comparing this figure with the previous Figure 6, they gave the same trend. On the other hand, the varying in the first range covers a longer range of energy in comparison with the second one. Moreover, the allowable barrier width range in this second range is shifted left to be from 1.98 nm to 2.7 nm. The energy gap between the IB and each of VB and CB is illustrated in Figures 9(a) and 9(b). As a result of changing the values of , the corresponding values of energy gaps will be changed for the composition of the model, QD width (InAs0.9N0.1) and barrier width (GaAs0.98Sb0.02). Also the curves in this case tend to be straight in comparison with their equivalent in Figures 7(a) and 7(b). One can notice that the second range investigations were not considered before in previous studies. Although it gives an enhancement of power conversion efficiency as shown in the next results.
An example of the J-V characteristics for the second range of the proposed model is mentioned in Figure 10. It is important to note that the short circuit current density is related to the quantum dot width directly; meanwhile the open circuit voltage is approximately constant. Additionally, the wide range of the J-V characteristics with short circuit current and open circuit voltage values is confirmed in Table 2. This will contribute to the enhancement of obtained power conversion efficiency as shown in Figure 11. One can notice that the behavior of efficiency is in full agreement with the obtained results in . But the new contribution is that the efficiency in this case, full concentration, reached up to (70.4%) (see also Table 2). The observed enhancement of the efficiency returns to discovering second range. Comparison between J-V characteristics for fully and unconcentrated cases are described in Figure 12. From this figure and Table 2, one can recognize that the second range denotes higher values of and a small amount of increase of in comparison to the first range. But in the two considered ranges, the values of in the full concentration case are higher than that in the unconcentrated case. From a solar cell point of view, this behavior normally happens because the allowable large width of quantum dots will acquire high photons and it then excites a large number of electrons: high induced current density. For the same cases considered before, the comparison between the power conversion efficiency is depicted in Figure 13. Also, the corresponding power conversion efficiency, , , , , , and fill factor, , for each combination of QD width () and barrier width () for fully and unconcentration are considered in Table 2. The fill factor pointed in the table is determined from (13). It also indicates an enhancement of power conversion efficiency, where is greater than 93% and 90% in the considered cases, respectively. The main target of these demonstrations is to compare the obtained power conversion efficiency from the first and second ranges. When a full concentration case is taken in the theoretical calculation, higher values of each of the open circuit voltage and short circuit current density are obtained that correspondingly gave higher power conversion efficiency. When one held comparison between our obtained results and others in the fully and unconcentrated cases: rather than 63.2% [11, 14, 26, 27], and rather than of 46.8%, respectively, . The numerical results in Table 2 can be utilized also for experimental achievements. In further investigations, the effect of two intermediate bands for different alloys into the power conversion efficiency will be studied.
In this work, a theoretical study for the QDIBSC that would contain one intermediate band between the valence and conduction bands was investigated. This intermediate band exploits the low energy photons to enhance the transferring process of charge carriers from valence to conduction bands. Otherwise, it would normally be lost as heat or not absorbed by the conventional solar cell. As a sequence of utilizing the proposed model, the incoming photons would be well-matched with the energies among bands; thereby the thermalization is reduced. From the obtained results, the open circuit voltage is approximately constant while the short circuit current density is increased. As a result, the fill factor will be enhanced by more than 90%. The power conversion efficiency was calculated under two different ranges of QD width, specified by three different values of barrier width; that is, from 4–7 nm, at values 2, 2.5, and 3 nm and from 7–11 nm, at values 1.98, 2.34, and 2.7 nm, respectively. The power conversion efficiency is affected by changing the QD and barrier widths. It reaches the maximum value, 63.2%, for full concentration light case in the first range when QD width is 4 nm and barrier width is 2 nm, which agrees with previous studies. The new facts that have not been addressed before are (a) for the second range, power conversion efficiency reaches 70.4% when QD width is 8.1 nm and the barrier width is 1.98 nm in the case of full concentration and (b) the maximum efficiency for unconcentrated light case is 51.5% at the first range and 57.5% for another range. One can recognize from the obtained results an enhancement for efficiency in each of the unconcentrated cases into two ranges and the case of full concentration in the second range. These results of efficiencies create more motivation for further research in this type of structure; therefore future work will concentrate on two intermediate bands: QD solar cell and other alloys with different compositions.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
- W. Shockley and H. J. Queisser, “Detailed balance limit of efficiency of p-n junction solar cells,” Journal of Applied Physics, vol. 32, no. 3, pp. 510–519, 1961.
- A. Martí, E. Antolín, E. Cánovas, N. López, and A. Luque, “Progress in quantum-dot intermediate band solar cell research,” in Proceedings of the 21st European Photovoltaic Solar Energy Conference, Dresden, Germany, 2006.
- A. Nasr and A. Aly, “Theoretical investigation of some parameters into the behavior of quantum dot solar cells,” under publication.
- A. Nasr, “Theoretical study of the photocurrent performance into quantum dot solar cells,” Optics & Laser Technology, vol. 48, pp. 135–140, 2013.
- A. Luque and A. Martí, “Increasing the efficiency of ideal solar cells by photon induced transitions at intermediate levels,” Physical Review Letters, vol. 78, no. 26, pp. 5014–5017, 1997.
- A. Luque, P. G. Linares, E. Antolín et al., “Understanding the operation of quantum dot intermediate band solar cells,” Journal of Applied Physics, vol. 111, Article ID 044502, 2012.
- P. G. Linaresa, C. D. Farmerb, E. Antolína et al., “ As quatemary alloys for quantum dot intermediate band solar cells,” Energy Procedia, vol. 2, no. 1, pp. 133–141, 2009.
- C. - Lin, W.-L. Liu, and C.-Y. Shih, “Detailed balance model for intermediate band solar cells with photon conservation,” in Proceedings of the 37th IEEE Photovoltaic Specialists Conference (PVSC '11), 2011.
- L. Cuadra, A. Martí, and A. Luque, “Present status of intermediate band solar cell research,” Thin Solid Films, vol. 451-452, pp. 593–599, 2004.
- M. Y. Levy, C. Honsberg, A. Marti, and A. Luque, “Quantum dot intermediate band solar cell material systems with negligable valence band offsets,” in Proceedings of the 31st IEEE Photovoltaic Specialists Conference, pp. 90–93, January 2005.
- Q. Deng, X. Wang, C. Yang et al., “Theoretical study on N/GaN quantum dots solar cell,” Physica B: Condensed Matter, vol. 406, no. 1, pp. 73–76, 2011.
- O. L. Lazarenkova and A. A. Balandin, “Miniband formation in a quantum dot crystal,” Journal of Applied Physics, vol. 89, no. 10, pp. 5509–5515, 2001.
- Q. Shao, A. A. Balandin, A. I. Fedoseyev, and M. Turowski, “Intermediate-band solar cells based on quantum dot supracrystals,” Applied Physics Letters, vol. 91, no. 16, Article ID 163503, 2007.
- E. J. Steven, Quantum dot intermediate band solar cells: design criteria and optimal materials [Ph.D. thesis], Drexel University, Philadelphia, Pa, USA, 2012.
- T. Soga, Nanostructured Materials for Solar Energy Conversion, 1st edition, 2006.
- S. Birner, Modeling of Semiconductor Nanostructures and Semiconductor-Electrolyte Interfaces, 2011.
- S. P. Day, H. Zhou, and D. L. Pulfrey, “The Kronig-Penney approximation: may it live on,” IEEE Transactions on Education, vol. 33, no. 4, pp. 355–358, 1990.
- R. Aguinaldo, Modeling solutions and simulations for advanced III-V photovoltaics based on nanostructures [M.S. thesis in materials science & engineering], College of Science, Rochester Institute of Technology, Rochester, NY, USA, 2008.
- W. van Roosbroeck and W. Shockley, “Photon-radiative recombination of electrons and holes in germanium,” Physical Review, vol. 94, no. 6, pp. 1558–1560, 1954.
- A. Martí, L. Cuadra, and A. Luque, “Quasi-drift diffusion model for the quantum dot intermediate band solar cell,” IEEE Transactions on Electron Devices, vol. 49, no. 9, pp. 1632–1639, 2002.
- L. A. Kosyachenko, Solar Cells—New Aspects and Solutions, 2011.
- S. P. Bremner and C. B. Honsberg, “Intermediate band solar cell with non-ideal band structure under AMl.5 spectrum,” in Proceedings of the 38th IEEE Photovoltaic Specialists Conference (PVSC '12), 2012.
- A. Martí, L. Cuadra, and A. Luque, “Partial filling of a quantum dot intermediate band for solar cells,” IEEE Transactions on Electron Devices, vol. 48, no. 10, pp. 2394–2399, 2001.
- J. Ojajärvi, Tetrahedral chalcopyrite quantum dots in solar-cell applications [M.S. thesis], Department of Physics, University of Jyväskylä, 2010.
- A. Luque, A. Martí, and L. Cuadra, “Impact-ionization-assisted intermediate band solar cell,” IEEE Transactions on Electron Devices, vol. 50, no. 2, pp. 447–454, 2003.
- Q.-W. Deng, X.-L. Wang, C.-B. Yang et al., “Computational investigation of N/InN quantum-dot intermediate-band solar cell,” Chinese Physics Letters, vol. 28, no. 1, Article ID 018401, 2011.
- A. Nasr, “Theoretical model for observation of the conversion efficiency into quantum dot solar cells,” under publication.
Copyright © 2014 Abou El-Maaty Aly and Ashraf Nasr. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.