- About this Journal
- Abstracting and Indexing
- Aims and Scope
- Annual Issues
- Article Processing Charges
- Articles in Press
- Author Guidelines
- Bibliographic Information
- Citations to this Journal
- Contact Information
- Editorial Board
- Editorial Workflow
- Free eTOC Alerts
- Publication Ethics
- Reviewers Acknowledgment
- Submit a Manuscript
- Subscription Information
- Table of Contents
International Journal of Photoenergy
Volume 2013 (2013), Article ID 421076, 9 pages
Numerical Analysis of Copper-Indium-Gallium-Diselenide-Based Solar Cells by SCAPS-1D
1Laboratoire des Matériaux et Environnement (LA.M.E), UFR-SEA, Université de Ouagadougou, BP 7021, Ouaga 03, Burkina Faso
2Département de Physique, Faculté des Science, Université de Yaoundé I, BP 812, Yaoundé, Cameroon
Received 8 June 2013; Revised 5 August 2013; Accepted 7 August 2013
Academic Editor: Cooper Harold Langford
Copyright © 2013 S. Ouédraogo et al. 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.
We used a one-dimensional simulation program Solar Cell Capacitance Simulator in 1 Dimension (SCAPS-1D) to investigate Copper-Indium-Gallium-Diselenide- (CIGS-) based solar cells properties. Starting with a conventional ZnO-B/i-ZnO/CdS/CIGS structure, we simulated the parameters of current-voltage characteristics and showed how the absorber layer thickness, hole density, and band gap influence the short-circuit current density (), open-circuit voltage (), fill factor (FF), and efficiency of solar cell. Our simulation results showed that all electrical parameters are greatly affected by the absorber thickness (w) below 1000 nm, due to the increase of back-contact recombination and very poor absorption. Increasing hole density (p) or absorber band gap () improves and leads to high efficiency, which equals value of 16.1% when p = 1016 cm−3 and eV. In order to reduce back-contact recombination, the effect of a very thin layer with high band gap inserted near the back contact and acting as electrons reflector, the so-called back-electron reflector (EBR), has been investigated. The performances of the solar cells are significantly improved, when ultrathin absorbers (w < 500 nm) are used; the corresponding gain of due to the EBR is 3 mA/cm2. Our results are in good agreement with those reported in the literature from experiments.
Continuously increasing demand for photovoltaic (PV) modules and the need for low-cost PV options have stretched these advantages to the limit and have exposed some inherent disadvantages of c-Si technology, such as the scarcity of feedstock material, costly processing of materials, and device fabrication steps, as well as the inability for monolithic interconnections .
Thin films solar cells, mainly CIGS in this context, are enrolled as an alternative to the silicon technology. One of the advantages of thin films solar cells based on CIGS is the reduction of production cost compared to crystalline silicon sector, related to the power used during the deposition process. The recent performance in the laboratory for CIGS solar cell is 20.3% ; close to crystalline silicon whose performance is around 25%. Increasing this performance focused researchers, and it is becoming a central topic in the field of thin films solar cells. Despite the fact that there exists an extensive literature on CIGS solar cells, which has best performances compared to other thin film photovoltaic sector, scientific knowledge of this family of solar cells are not yet exhaustive. This concerns mainly the loss mechanisms in the cell, the substitution of the toxic CdS layer by other alternative layers, and the reduction of the absorber thickness beyond 1000 nm [3, 4]. In order to improve CIGS solar cell performances, it is necessary to increase the understanding of the basic factors limiting the electrical parameters of the cell.
The purpose of this work is to study, using SCAPS-1D  simulation package, factors limiting the performance of CIGS solar cells. We examine the influence of the absorber layer thickness, band gap, and hole density, as well as the effect of the back-electron reflector (EBR), on the electrical parameters of CIGS solar cell.
2. Materials and Methods
2.1. Cell Structure
The structure of the solar cell is (Ni/Al)/MgF2/ZnO:B/i-ZnO/CdS/OVC/CIGS/Mo/substrate (Figure 1). The key parts of the cell are the CIGS absorber and the CdS buffer layer. The layer between the CdS and CIGS absorber is a thin layer named ordered vacancy compound (OVC). This layer, is formed by the interface states between the buffer layer and the absorber . The OVC layer is considered to be beneficial to the performance of CIGS cells because the electrical junction is shifted away from the high-recombination interface between the CdS and CIGS layer, and hence, the recombination rate is reduced . A ZnO intrinsic layer (i-ZnO) and boron-doped ZnO (ZnO:B) layer are deposited on the top of the buffer layer. These two layers are commonly referred to as transparent conductive oxide (TCO), because of their wide band gap which makes them transparent to most of the solar spectrum. Most solar cells based on CIGS use ZnO:Al as TCO. But in this paper, we use the boron-doped ZnO as TCO. Indeed, the ZnO:Al has absorption losses, leading to a decrease of the quantum efficiency of the solar cells in the near infrared regions. A boron doping would be more beneficial for the solar cells . The TCO is covered with an antireflection layer MgF2, which increases the absorption of photons in the absorber.
2.2. Numerical Modeling and Material Parameters
CIGS polycrystalline solar cells are complicated structures, due to the large number of layers and the fact that the effects of some particular phenomena, mechanisms or material parameters often result from intuition. Numerical simulations can be used to provide insight to interpret measurements and to assess the potential merits of a cell structure. Indeed, once multiple measurements are (more or less) quantitatively described, the simulations can be used to analyse the effect of the variation of material parameters, that is, the presence or absence of particular properties, or variation of all properties in the range of values, to obtain the optimal values for optimizing the solar cells efficiencies, and should give the manufacturers additional ideas of how to vary their production methods to improve the product performance. Several software, among which SCAPS-1D , ASA , PC-1D , AFORS-HET , and AMPS-1D , have been developed in order to simulate the functioning of thin film solar cells. SCAPS is widely used for the simulation of CIGS- and CdTe- based solar cells. The good agreements between SCAPS-1D simulation results and the existing experimental ones  motivated us to use this tool in our work. SCAPS calculates the steady-state band diagram, recombination profile, and carrier transport in one dimension, based on Poisson equation together with hole and electron continuity equations. Recombination currents are calculated with the Shockley-Read-Hall (SRH) model for bulk defects and an extension of the SRH model for interface defects. To keep the model as simple as possible, one type of single level defects is introduced in each layer. These are all compensating defects positioned at the intrinsic level that is close to midgap. To pin the Fermi level at the absorber (CIGS)/OVC and OVC/(CdS) layer interface, neutral defects were placed 0.2 eV below the conduction band. These have a small capture cross-section to separate between pinning and recombination parameters of the OVC. The optical and electrical parameters used in this paper are derived from numerical models [5, 15, 16]. The influence of the series resistance and shunt resistance are not taken into consideration. Band discontinuity at the interface of the different materials is assumed small and neglected. The solar spectrum AM.1.5 is used for this numerical simulation, and the temperature is set at 300 K. Table 1 summarizes the parameters of the different layers used in this paper.
Figure 2 shows a superposition of J-V curve of experimental data  and a simulated CIGS solar cell for our baseline, where the absorber thickness taken as the default thickness is 1800 nm. The resulting performance parameters of the open-circuit voltage (), short-circuit current density (), fill factor (FF), and efficiency are determined using J-V characteristics and are shown in Table 2.
Although the experimental values are slightly higher than the simulated ones, excellent agreement was observed in the J-V curve. The agreement between experiment and simulation is good for the J-V characteristic and validates our set of parameters as a baseline for simulating the influence of the variation of absorber parameters on the solar cell performances.
3. Results and Discussion
3.1. Effect of Absorber Layer Thickness on the Solar Cell Characteristics
The standard thickness of the Cu(In,Ga) Se2 layer in CIGS solar cells is about 3000 nm. If this thickness could be reduced, with no or only minor loss in the performance, the deposition time of CIGS layer would be reduced for thinner CIGS layers. A thinner CIGS layer would reduce the direct materials usage and thereby the materials costs. A reduction of materials usage is important for indium (In) and gallium (Ga) since the supply of these metals might become an issue if CIGS thin-film solar cells are produced in large volumes . However, the reduction of the absorber thickness is associated with a number of problems [3, 4]. With SCAPS-1D, the properties of the different layers are kept constant while varying the absorber thickness, in order to obtain qualitative information of the absorber layer thickness on the solar cell electrical parameters. The band gap of the absorber is also kept constant to 1.15 eV. Figure 3 shows the absorber layer thickness () effect on the electrical parameters of the solar cells. For all electrical parameters (, , FF, and efficiency), we can distinguish two zones of the absorber thickness () which influence strongly the electrical parameters. The first zone is nm, and the second one corresponds to nm.
All electrical parameters decrease significantly in the first zone. The short-circuit current density () is the parameter which is most affected by the decrease of the absorber thickness. This is mainly due to the recombination of photogenerated electrons at the back contact (Mo) [3, 4]. For a thinner absorber, photons of short wavelength (higher energy) penetrate deeply into the absorber and generate electron holes near the back contact, which is a zone of high recombination, thereby resulting in the decrease of the current density. The passes from 29.4 mA/cm2 for nm to 5.3 mA/cm2 for nm, that is, a loss of 24.1 mA/cm2 corresponding to the recombination current at the back contact. The quantum efficiency curve (Figure 4) shows that the absorption of the incident light greatly reduces in the first zone of thickness ( nm); mainly the wavelength is greater than 500 nm, related to the combined effect of light transmission and high recombination of electrons at the CIGS/Mo interface. The decrease of the open-circuit voltage is related to degradation of the junction in the case of ultrathin thickness. The decrease of and leads to the decrease of the efficiency of the solar cells according to (1), where FF is the fill factor.
Consider The cell efficiency decreases from 14.13% for nm to 1.35% for nm. The decrease of the efficiency is even more abrupt when the thickness of the absorber exceeds 500 nm, due to the simultaneous action of degradation of the absorption and electrons capture by the back contact. When the thickness of the absorber decreases, fewer photons are absorbed in the absorber. The quantity of electron hole is reduced, which reduces the efficiency of the solar cell.
In the second zone ( nm), all electrical parameters are almost constant. However, the short-circuit current density passes from 29.3 mA/cm2 for nm to 31.53 mA/cm2 for nm, that is, a gain of 2.23 mA/cm2, leading the cell efficiency to 16.10%. This is related to the increase in the collection of photogenerated carriers and complete absorption of the photons, as shown in Figure 4 (second zone). The quantum efficiency of the cells is high in this zone of thickness and reaches 90%. More photons are absorbed for nm, which increased the solar cell performances. Attempts are made to compare the simulation results with reported experimental data for the CIGS cells. The results are shown in Figure 3, where the solid red circles denote the experimental values extracted from . Excellent agreement was obtained between the simulated solar cell electrical parameters (, , FF and efficiency) trends and experimental values. Although the experimental values are slightly higher than the simulated ones because an intentional high defect density is introduced into the layers, we find that the main tendencies of the experimental data are reproduced by the simulation.
In sum, and decrease very significantly for thicknesses nm. These observations are also similar to those of references .
3.2. Effect of Absorber Hole Density
Figure 5 shows the influence of hole density on , , FF, and efficiency, for different thicknesses of the absorber. The open-circuit voltage () increases significantly with increasing absorber doping. This increase is particularly important when the thickness of the absorber is larger. The reaches a peak at cm−3, beyond this value, there is a saturation of the open-circuit voltage independently of the absorber thickness.
This saturation is related to the dependence of the space charge region width (SCRW) and the hole density in the absorber. If we assume that no voltage is applied to the diode and the build-in potential is equal to unity, the SCRW is given by approximation as the following equation: where is the relative permittivity of the absorber, the absolute vacuum permittivity, and the elementary charge. Increasing hole density () reduces the space charge region width, and thus, causes the saturation of the for large values of . For nm, the thickness of the absorber may become order of magnitude or smaller than the space charge region width and tends to be fully depleted . In this case, the effect of doping on the open-circuit voltage remains insignificant. The short-circuit current density follows the same trend for all thicknesses of the absorber. It decreases significantly with increasing hole density. However, the decrease of is more brutal for nm. For this range of thickness, increasing the doping may lead to a reduction of the space charge width below the thickness of the absorber, thereby reducing the collection of the carriers by the junction.
The fill factor is affected by the increase of the hole density. For all thicknesses, the voltage gain is higher than the loss of the ; thus, the overall efficiency of the device increases with doping and reaches a peak of 16.10% for cm−3. Beyond cm−3, the combined effect of the saturation of and the brutal drop of for all thicknesses makes the performance of the device almost constant independently of doping.
3.3. Effect of Ga-Grading on the Cell Performance
The band-gap of CIGS can be tuned from 1.04 eV (pure CIS) to 1.65 (pure CGSe) depending on the Ga content. The Ga dependence of the band gap follows the equation  where denotes the proportion of gallium in the absorber, that is, the ratio Ga/(Ga + In), is the optical bowing coefficient for which values between 0.11 and 0.24 have been reported .
Thereafter, we assume that the band gap of the absorber , bulk defect densities , absorption coefficients, and electron affinities  is varied according with the Ga contents in the absorber. The other parameters of the absorber as well as the properties of the other layers of the cells are assumed constant in this simulation. With a uniform band gap profile, we analyze the effect of increasing absorber band gap on the electrical parameters.
As shown in Figure 6, the open-circuit voltage () increases with the increase of the band gap, almost independently of the absorber thickness. However, this increase is not proportional to the band gap. The dependence of with the band gap is linear for eV, corresponding to , and less linear . The short-circuit current density decreases dramatically with increasing gallium concentration but strongly depends on the thickness of the absorber. For nm, the decrease of is brutal, due to the combined effect of the decrease of absorption coefficient with increasing the band gap  and the reduction of absorption in the long wavelength region of the solar spectrum in the thin absorber. These results are also confirmed in [22, 23]. The efficiency of the solar cell as well as the fill factor (FF) increases with the band gap for eV; beyond this value, the increase of the band gap by introducing the gallium has no effect and becomes detrimental to the performance of the device. These results are in good agreement with the experimental results from the literature. The best solar cell was obtained with a ratio = Ga/(In + Ga) = 0.3 , which corresponds to a band gap of 1.15 eV. Despite attempts to increase the band gap and therefore increase the , CIGS solar cell shows poor electrical characteristics when the rate of gallium exceeds 0.3%, that is, eV. Rau et al.  and Hanna et al.  have shown experimentally that the increase of the band gap by introducing gallium creates defects in the volume of the absorber. These defects become detrimental to the performance of the device when the gap exceeds 1.15 eV ( = 30%) corresponding to the minimum defect density .
3.4. Potential Improvement by Using Back-Electron Reflector
The short-circuit current density () is the most affected parameter by the variation of the thickness, especially for ultrathin absorber ( nm), where the capture of electrons by the back contact is established as the main source of reduction of the [3, 4, 26].
At the CIGS/Mo interface, it would be desirable to keep the photoelectron away from this interface, which is expected to have a relatively high recombination velocity. By using a very thin layer rich in gallium (Ga) at the CIGS/Mo interface, commonly called back-electron reflector (EBR), we can keep the high conductivity for the majority holes and at the same time reflect the minority electrons.
To better understand the beneficial effect of the EBR on the electrical parameters, all the parameters of layers are kept constant, except the band gap and the thickness of the absorber. The properties of EBR layer are summarized in Table 3. In order to obtain qualitative information on the beneficial effect of the EBR, we show in Figure 7, the quantities Δ, Δ, ΔFF, and Δefficiency, which represent the electrical parameters gained with the EBR, as a function of absorber thickness and absorber band gap.
Except the fill factor (FF), which shows a loss for eV, the other electrical parameters have a gain due to the EBR.
The gain is very important when the thickness of the absorber is very small ( nm). For thickness greater than 1000 nm, the effect of EBR on the electrical parameters is insignificant, given the fact that the thickness of the absorber is thick enough that the absorption takes place in the bulk of the absorber, away from the interface CIGS/Mo. However, when the absorber is thin ( nm), photons of long wavelengths pass through the absorber and generate carriers at the CIGS/Mo interface, which is a zone of high recombination. The presence of the EBR can repel the electrons away from this interface, avoiding their capture by the Mo, and the short-circuit current density () increases. For thin thicknesses, the gain of due to the EBR is 3 mA/cm2 on average, and it is less than 0.5 mA/cm2 for thickness greater than 1000 nm. These results agree very well with those of Kanevce , who also showed that the gain of due to the EBR for ultrathin absorber is about 3 mA/cm2.
Using SCAPS-1D package, we analyzed the variation of the absorber layer thickness, absorber holes density, the band gap, and the effect of the introduction of the EBR on the electrical parameters of a CIGS solar cells. We have shown the following.(i)Electrical parameters of CIGS solar cell are affected by reducing absorber thickness, but the most significant loss is the short-circuit current density that shows a loss of 24.1 mA/cm2 when the thickness decreases from 1000 to 100 nm, due to the increasing of recombination at the CIGS/Mo interface.(ii)Increasing hole density in the absorber can significantly increase the performance of the device. This increase is mainly due to the gain of the open-circuit voltage () when the hole density is lower than 1016 cm−3.(iii)Ga grading can increase the performance of the solar cells. However, the short-circuit current density decreases considerably with increasing Ga in the absorber but strongly depends on the absorber layer thickness. The dependence of the with the band gap is linear for eV.(iv)The use of the EBR at CIGS/Mo interface can increase the short-circuit current density, especially, when the absorber thickness is less than 1000 nm. The gain due to the EBR is estimated at 3 mA/cm2.
|AM.1.5:||Standard terrestrial solar spectrum “Air Mass 1.5”|
|SCAPS-1D:||Solar Cell Capacitance Simulator in 1 Dimension|
|CIGS:||Copper Indium Gallium Diselenide|
|SCRW:||Space Charge Region Width|
|:||Band-gap energy of the semiconductor|
|:||Short-circuit current density|
|, :||Electron and hole mobility|
|, :||Effective density of states in the conduction and valence band|
|, :||Electron and hole capture cross section|
|, :||Electron and hole thermal velocity|
|:||Absorber layer width|
The authors acknowledge the use of SCAPS-1D program developed by Marc Burgelman and colleagues at the University of Gent in all the simulations reported in this paper. The stay of S. Ouédraogo at the University of Yaoundé I was supported by a fellowships from PIMASO, a program financed by the European Union.
- T. M. Razykov, C. S. Ferekides, D. Morel, E. Stefanakos, H. S. Ullal, and H. M. Upadhyaya, “Solar photovoltaic electricity: current status and future prospects,” Solar Energy, vol. 85, no. 8, pp. 1580–1608, 2011.
- P. Jackson, D. Hariskos, E. Lotter et al., “New world record efficiency for Cu(In,Ga)Se2 thin-film solar cells beyond 20%,” Progress in Photovoltaics, vol. 19, no. 7, pp. 894–897, 2011.
- M. Gloeckler and J. R. Sites, “Potential of submicrometer thickness Cu (In,Ga) Se2 solar cells,” Journal of Applied Physics, vol. 98, no. 10, Article ID 103703, pp. 1–7, 2005.
- Z. Jehl, F. Erfurth, N. Naghavi et al., “Thinning of CIGS solar cells: part II: cell characterizations,” Thin Solid Films, vol. 519, no. 21, pp. 7212–7215, 2011.
- J. Pettersson, C. Platzer-Björkman, U. Zimmermann, and M. Edoff, “Baseline model of graded-absorber Cu(In,Ga)Se2 solar cells applied to cells with buffer layers,” Thin Solid Films, vol. 519, no. 21, pp. 7476–7480, 2011.
- S. Degrave, M. Burgelman, and P. Nollet, “Modelling of polycrystalline thin film solar cells: new features in scaps version 2.3,” in Proceddings of the 3rd World Conference on Photovoltaic Energy Conversion, pp. 487–490, May 2003.
- A. Bouloufa, K. Djessas, and A. Zegadi, “Numerical simulation of solar cells by AMPS-1D,” Thin Solid Films, vol. 515, no. 15, pp. 6285–6287, 2007.
- J. Song, S. S. Li, C. H. Huang, O. D. Crisalle, and T. J. Anderson, “Device modeling and simulation of the performance of Cu solar cells,” Solid-State Electronics, vol. 48, no. 1, pp. 73–79, 2004.
- Y. Hagiwara, T. Nakada, and A. Kunioka, “Improved Jsc in CIGS thin film solar cells using a transparent conducting ZnO:B window layer,” Solar Energy Materials and Solar Cells, vol. 67, no. 1–4, pp. 267–271, 2001.
- A. M. K. Dagamseh, B. Vet, P. Šutta, and M. Zeman, “Modelling and optimization of a-Si:H solar cells with ZnO:Al back reflector,” Solar Energy Materials and Solar Cells, vol. 94, no. 12, pp. 2119–2123, 2010.
- P. A. Basore and D. A. Clugston, “PC1D version 4 for windows: from analysis to design,” in Proceedings of the 25th IEEE Photovoltaic Specialists Conference, pp. 377–381, May 1996.
- S. Singh, S. Kumar, and N. Dwivedi, “Band gap optimization of p-i-n layers of a-Si:H by computer aided simulation for development of efficient solar cell,” Solar Energy, vol. 86, no. 5, pp. 1470–1476, 2012.
- N. Hernández-Como and A. Morales-Acevedo, “Simulation of hetero-junction silicon solar cells with AMPS-1D,” Solar Energy Materials and Solar Cells, vol. 94, no. 1, pp. 62–67, 2010.
- A. Niemegeers, S. Gillis, and M. Burgelman, “A user program for realistic simulation of polycristalline heterojunction solar cells: SCAPS-1D,” in Proceedings of the 2nd World Conference on Photovoltaic Energy Conversion, Wien, July 1998.
- R. Scheer, “Towards an electronic model for solar cells,” Thin Solid Films, vol. 519, no. 21, pp. 7472–7475, 2011.
- M. Gloeckler, A. L. Fahrenbruch, and J. R. Sites, “Numerical modeling of CIGS and CdTe solar cells: setting the baseline,” in Proceddings of the 3rd World Conference on Photovoltaic Energy Conversion, pp. 491–494, May 2003.
- O. Lundberg, M. Bodegård, J. Malmström, and L. Stolt, “Influence of the Cu(In, Ga)Se2 thickness and Ga grading on solar cell performance,” Progress in Photovoltaics, vol. 11, no. 2, pp. 77–88, 2003.
- P. Chelvanathan, M. I. Hossain, and N. Amin, “Performance analysis of copper-indium-gallium-diselenide (CIGS) solar cells with various buffer layers by SCAPS,” Current Applied Physics, vol. 10, no. 3, pp. S387–S391, 2010.
- O. Lundberg, M. Edoff, and L. Stolt, “The effect of Ga-grading in CIGS thin film solar cells,” Thin Solid Films, vol. 480-481, pp. 520–525, 2005.
- G. Hanna, A. Jasenek, U. Rau, and H. W. Schock, “Influence of the Ga-content on the bulk defect densities of Cu(In,Ga)Se2,” Thin Solid Films, vol. 387, no. 1-2, pp. 71–73, 2001.
- T. Minemoto, T. Matsui, H. Takakura et al., “Theoretical analysis of the effect of conduction band offset of window/CIS layers on performance of CIS solar cells using device simulation,” Solar Energy Materials and Solar Cells, vol. 67, no. 1–4, pp. 83–88, 2001.
- C.-H. Huang, “Effects of Ga content on Cu(In,Ga)Se2 solar cells studied by numerical modeling,” Journal of Physics and Chemistry of Solids, vol. 69, no. 2-3, pp. 330–334, 2008.
- R. Kniese, D. Hariskos, G. Voorwinden, U. Rau, and M. Powalla, “High band gap Cu(In,Ga)Se2 solar cells and modules prepared with in-line co-evaporation,” Thin Solid Films, vol. 431-432, pp. 543–547, 2003.
- N. Amin, P. Chelvanathan, M. I. Hossain, and K. Sopian, “Numerical modelling of ultra thin Cu(In,Ga)Se2 solar cells,” in Proceedings of the 6th International Conference on Materials for Advanced Technologies (ICMAT '11), pp. 291–298, July 2011.
- U. Rau, M. Schmidt, A. Jasenek, G. Hanna, and H. W. Schock, “Electrical characterization of Cu(In,Ga)Se2 thin-film solar cells and the role of defects for the device performance,” Solar Energy Materials and Solar Cells, vol. 67, no. 1–4, pp. 137–143, 2001.
- Z. J. L. Kao, N. Naghavi, F. Erfurth et al., “Towards ultrathin copper indium gallium diselenide solar cells: Proof of concept study by chemical etching and gold back contact engineering,” Progress in Photovoltaics, vol. 20, pp. 582–587, 2012.
- A. Kanevce, Anticipated performance of Cu(In, Ga)Se2 solar cells in the thin-film limit [Doctoral thesis], Colorado State University, 2007.