Research article Special Issues

Quantitative characterization and modeling of sub-bandgap absorption features in thin oxide films from spectroscopic ellipsometry data

  • Analytic representations of the complex dielectric function, which describe various types of materials, are needed for the analysis of optical measurements, in particularly, ellipsometric data. Here, we examine an improved multi-oscillator Tauc-Lorentz (TL) model with a constraint on the band-gap parameter Eg, which forces it to be common for all TL oscillators, and possibility to represent reasonably weak absorption features below the bandgap by inclusion of additional unbounded Lorentz and/or Gaussian oscillators with transition energies located below Eg. We conclude that the proposed model is the most appropriate for the characterization of various materials with sub-band absorption features and provides meaningful value for the energy bandgap. A few examples to illustrate the use of modified model have been provided.

    Citation: Dmitriy V. Likhachev, Natalia Malkova, Leonid Poslavsky. Quantitative characterization and modeling of sub-bandgap absorption features in thin oxide films from spectroscopic ellipsometry data[J]. AIMS Materials Science, 2015, 2(4): 356-368. doi: 10.3934/matersci.2015.4.356

    Related Papers:

    [1] Bitian Wu . Day ahead scheduling model of wind power system based on fuzzy stochastic chance constraints—considering source-load dual-side uncertainty case. AIMS Energy, 2025, 13(3): 471-492. doi: 10.3934/energy.2025018
    [2] Sameer Thakare, Neeraj Dhanraj Bokde, Andrés E. Feijóo-Lorenzo . Forecasting different dimensions of liquidity in the intraday electricity markets: A review. AIMS Energy, 2023, 11(5): 918-959. doi: 10.3934/energy.2023044
    [3] Gerardo Guerra, Juan A. Martinez Velasco . A virtual power plant model for time-driven power flow calculations. AIMS Energy, 2017, 5(6): 887-911. doi: 10.3934/energy.2017.6.887
    [4] Vishal Chaudhary, Hari Mohan Dubey, Manjaree Pandit, Surender Reddy Salkuti . A chaotic Jaya algorithm for environmental economic dispatch incorporating wind and solar power. AIMS Energy, 2024, 12(1): 1-30. doi: 10.3934/energy.2024001
    [5] Yan Li, Yaheng Su, Qixin Zhao, Bala Wuda, Kaibo Qu, Lei Tang . An electricity price optimization model considering time-of-use and active distribution network efficiency improvements. AIMS Energy, 2025, 13(1): 13-34. doi: 10.3934/energy.2025002
    [6] Andrew A. Chien, Fan Yang, Chaojie Zhang . Characterizing curtailed and uneconomic renewable power in the mid-continent independent system operator. AIMS Energy, 2018, 6(2): 376-401. doi: 10.3934/energy.2018.2.376
    [7] Gerardo Guerra, Juan A. Martinez-Velasco . Optimal sizing and operation of energy storage systems considering long term assessment. AIMS Energy, 2018, 6(1): 70-96. doi: 10.3934/energy.2018.1.70
    [8] V. V. Thang, Thanhtung Ha . Optimal siting and sizing of renewable sources in distribution system planning based on life cycle cost and considering uncertainties. AIMS Energy, 2019, 7(2): 211-226. doi: 10.3934/energy.2019.2.211
    [9] Chiara D’Alpaos, Michele Moretto . Do smart grid innovations affect real estate market values?. AIMS Energy, 2019, 7(2): 141-150. doi: 10.3934/energy.2019.2.141
    [10] Dalong Guo, Chi Zhou . Potential performance analysis and future trend prediction of electric vehicle with V2G/V2H/V2B capability. AIMS Energy, 2016, 4(2): 331-346. doi: 10.3934/energy.2016.2.331
  • Analytic representations of the complex dielectric function, which describe various types of materials, are needed for the analysis of optical measurements, in particularly, ellipsometric data. Here, we examine an improved multi-oscillator Tauc-Lorentz (TL) model with a constraint on the band-gap parameter Eg, which forces it to be common for all TL oscillators, and possibility to represent reasonably weak absorption features below the bandgap by inclusion of additional unbounded Lorentz and/or Gaussian oscillators with transition energies located below Eg. We conclude that the proposed model is the most appropriate for the characterization of various materials with sub-band absorption features and provides meaningful value for the energy bandgap. A few examples to illustrate the use of modified model have been provided.


    1. Introduction

    There are many isolated islands in Japan which are not connected to the main grid. The power systems in these islands mainly depend on diesel generators (DGs), which is not economically viable. However, these remote islands are rich in natural resources and have good geographical locations that make them suitable for renewable energy generation. Therefore, renewable energy systems (RESs) such as photovoltaic (PV) and wind generator (WG) systems have been introduced [1,2,3,4,5,6,7,8]. However, generated power by RESs varies according to the time and season as well as the weather conditions such as solar radiatoin and wind speed [9,10,11]. In case where PV and WG are introduced into small power systems, the power systems are affected because the rate of RESs to load demand is relatively large. Also, there are cases wherein the operational plan created on day-ahead cannot work because of forecasting error. In such scenario, the operational plan should be revisited based on new weather information. As a measure against the generated power fluctuation of RESs, battery energy storage system (BESSs) are often utilized. BESSs are very effective system because they allow the generated power fluctuation to be suppressed by charging and discharging [12,13,14,15,16].

    Moreover, demand response (DR) programs are attracting much attention these days [17,18,19,20,21,22,23,24,25]. DR programs enable Power Companies to control the load consumption of customers, and they are divided into time-based and incentive-based programs. Several literatures have developed DR models [19,20,21,22]. Furthermore, the effectiveness of DR has been investigated and the benefit it makes to the supply-side by using demand-side resources [22,23,24,25]. However, this paper utilizes real-time pricing (RTP) as a time-based programs to correspond to the fluctuation of powers generated by RESs. RTP encourages customer to increase or reduce the load consumption by varying electricity prices.

    With the power system reform in Japan, the power markets transactions are also increasing. In particular, spot (day-ahead) market, intraday hour-ahead market and forward market have been established by the Japan electric power exchange (JEPX). Day-ahead market is wherein power transaction is conducted based on the power generation and demand plan of next day, and hour-ahead market is wherein you adjusts the plan in intraday. In the forward market, the power transaction is performed on fixed price for the relatively far future [26]. These markets make it possible to deal with uncertainty related to the generated power fluctuation [27,28,29,30,31].

    This paper presents the optimal operation method coping with the forecasting errors caused by the generated power of the RESs considering the reforecasting and re-planning, the RTP, and two kinds of power markets in multi-area small power systems. The purpose of this research is to indicate the cooperative optimal opration method for a group of isolated small power systems which have RESs and BESS to relieve a economical burden. Furthermore, the process of scheduling the plan that contains RTP and two kinds of power markets is presented, as a measure to cope with uncertainty of renewable energies. The ultimate goal of this study is to reduce the total operational cost for a group of the power systems.Two kinds of power markets are considered for power transactions in three small power systems. The power transactions for the next day are performed in the day-ahead market, and the power transactions in the hour-ahead market are conducted every 3 hours in focused day. However, reforecasing and re-planning are carried out at the same time to reduce forecasting error. Furthermore, the RTP is applied to level the load profiles by varying the electricity prices every hour. The proposed method has achieved to reduce the total operational cost of the multi-area small power systems by using the surplus power obtained through two kinds of power markets effectively even if the forecasting errors of the power generated by RESs and the load demand occur.

    The rest of the paper is organized as follows. Section 2 describes the multi-area small power systems model and scenario-based method. The transactions in two mentioned markets and the RTP are discussed in Sections 3 and 4, respectively. In Section 5, the proposed operation method is explained in detail. Section 6 indicates the simulation results and discussion. Finally, conclusion and future work are given in Section 7.


    2. Description of Power System Model

    The multi-area small power systems model assumed in this paper is illustrated in Fig. 1. It can be seen that each power system performs a power trading via power markets. Table 1 shows the installed capacity in each power system (PS). Furthermore, Fig. 2 shows the configuration of resources in the power system 1 (PS1) which includes two types of DG, PV, and BESS. The power generated by PV and load profile in PS1 are indicated by solid lines in Figs. 3(b) and 3(c). It is assumed that the hourly prediction data of demand, wind, and irradiance are obtained, and updated every 3 hours. However, uncertainties are considered by using scenario-based method. Fig. 3 shows normal distributions and the load demand and the power generated by PV. The forecasting error increases in proportion to time, and thus, the normal distributions for making scenario are assumed as shown in Fig. 3(a). In Figs. 3(b) and 3(c), the dotted lines represent the scenario examples of the load demand and the power generated by PV.

    Figure 1. Multi-area small power systems model.
    Table 1. Installed capacity of each PS.
    PS1 PS2 PS3
    Max Load(MW) 1.0 0.85 0.85
    DG1(0.5MW/unit) 1 1 1
    DG2(0.35MW/unit) 2 1 1
    PV(MW) 0.5 0.5 0.5
    WG(MW) 0.5 0.75
    BESS((0.5MW/5.0MWh)/unit) 1 1 1
     | Show Table
    DownLoad: CSV
    Figure 2. Configuration of Power System 1 (PS1).
    Figure 3. Making scenario.

    3. Power Markets

    As mentioned earlier, the day-ahead and hour-ahead markets are considered. This paper determines the bidding price and bidding power by using fuzzy inference. Furthermore, the membership function used in bidding is shown on the following equations.

    $ SF_{se}=k_{se}P_{dl}~~(P_{dl} \leq 0) $ (3.1)
    $ CRh_{se}=CRh_{se}+SF_{se} $ (3.2)
    $ CRl_{se}=CRl_{se}+SF_{se} $ (3.3)
    $ MS_{se}=\frac{(PR_{ma}-CRl_{se})}{CRh_{se}-CRl_{se}}~~~(CRl_{se} < PR_{ma} \leq CRh_{se}) $ (3.4)

    where, $SF_{se}$ is shift function for selling; $k_{se}$ is coefficient; $P_{dl}$ is the difference between load demand and generation output (expected surplus power); $CRh_{se}$ is maximum criterion price; $CRl_{se}$ is minimum criterion price; $MS_{se}$ is membership function; $PR_{ma}$ is the trading price. First, the validities obtained from membership function are evaluated as shown in Figs. 4(a) and 4(b) and then, the bidding price and bidding power are determined in Fig. 4(c) by using center-of-gravity method. As a decision criterion for the bidding price and bidding power, the minimum difference between the bidding powers in each PS is employed. Moreover, these equations are:

    $ min: \sum\limits_{PS=1}^{3}{P_{bid}^{PS}(t)} $ (3.5)
    Figure 4. Membership functions.
    $ PPSdl(t)=kl×(PPSload(t)Ni=1PPSDGi(t))+kBE×MSPSBE×PPSBE
    $
    (3.6)
    $ P_{bid}^{PS}(t)=MS_{se}^{'PS} \times P_{dl}^{PS}(t) $ (3.7)

    where $P_{bid}^{PS}(t)$ is the bidding power in $PS$ th power system at time $t$; $P_{dl}^{PS}(t)$ is the difference between load demand and generation output (expected surplus power) in $PS$ th power system at time $t$; $k_{l}$ is the weight for load demand; $P_{load}^{'PS}$ is the difference between load demand and RES output in $PS$ th power system at time $t$; $N$ is the number of DG introduced into power systems; $P_{DGi}^{PS}(t)$ is the output power of the $i$ th DG in $PS$ th power system at time $t$; $k_{BE}$ is the weight for the state of charge (SOC); $MS_{BE}^{'PS}$ is the validity of SOC in $PS$ th power system; $P_{BE}^{PS}$ is the charging and discharging power of BESS in $PS$ th power system at time $t$; $MS_{se}^{'PS}$ is the validity of the expected surplus power.


    4. Real-time Pricing

    The sigmoid function is utilized as a demand response model of customer to the electricity price. The load demands are levelized by the RTP shown in Fig. 5 and Eq. (4.1). The equation of demand amount to electricity price is:

    $ h_{rtp}(t)=\frac{1}{1+\exp (PR_{rtp}(t)-SF_{1})}+SF_{2} $ (4.1)
    $ \Delta P_{Load}(t)=h_{rtp}(t) \times \left(P'_{Load}(t)-\frac{\sum\limits_{t=1}^T{P'_{Load}(t)}}{T}\right) \label{response_amount_kWh} $ (4.2)
    $ P_{RTPLoad}(t)=P_{Load}(t)+\Delta P_{Load}(t) $ (4.3)
    Figure 5. Sigmoid function.

    where, $h_{rtp}(t)$ is the variation rate of load demand at time $t$; $PR_{rtp}(t)$ is the set price at time $t$; $SF_{1}$ is the left-right shift coefficient; $SF_{2}$ is the top-bottom shift coefficient; $\Delta P_{Load}(t)$ is the response amount at time $t$; $P'_{Load}(t)$ is the difference between load demand and RES output at time $t$; $P_{RTPLoad}(t)$ is the difference between load demand and RES output at time $t$ after RTP. The response amounts are determined based on this sigmoid function.


    5. Proposed Operation Method

    The flow of operational method proposed in this paper is shown in Fig. 6. Planning the operational plan and bidding to the day-ahead market are performed by using the forecasting load in the preceding day. In the forecast day, the reforecasting and re-planning, RTP, and transactions in the hour-ahead market are performed every 3 hour by calculating the shortage and surplus power caused by the forecasting errors.

    Figure 6. Operation scheduling flow.

    This section describes the objective function, constraints and optimization method to reduce the total operational cost of multi-power systems.


    5.1. Objective Function

    The objective function is to minimize the total operational cost of multi-power systems in this paper. The start/stop and output plan for the DG, the charging and discharging plan for the BESS, the demand response plan, and interchanging power are determined such that the objective function is minimized. Furthermore, the objective function and the total cost of multi-power systems are given by

    $ min~:~TC=\sum\limits_{s=1}^{100}p_s{(TC_{PS1}+TC_{PS2}+TC_{PS3})} $ (5.1)
    $ TC_{PS}=\sum\limits_{t=1}^{24}[\sum\limits_{i=1}^{N}{\Bigl(F_{i}(t)(P_{DGi}(t))+SU_{i}(t)+SD_{i}(t)\Bigr)}+C_{BESS}+C_{day}(t)+C_{hour}(t)] $ (5.2)

    where, $p_{s}$ is the probability of error; $s$ is the scinario number; $TC$ is the total cost; $F_{i}(t)(P_{DGi}(t))$ is the fuel cost [Yen] of the $i$ th DG at time $t$; $P_{DGi}(t)$ is the output power [kW] of the $i$ th DG at time $t$; $SU_{i}(t)$ and $SD_{i}(t)$ are the starting and stopping cost [Yen] of the $i$ th DG at time $t$; $C_{BESS}$ is the energy cost [Yen] of BESS used in one day; $C_{day}$ and $C_{hour}$ are trading prices in the day-ahead and hour-ahead markets, respectively.

    Furthermore, fuel cost is decided by:

    $ F_{i}(P_{DGi})=\frac{Q \times P_{DGi} \times PR_{c}}{Q_{c} \times \eta_{i}} $ (5.3)

    where, $Q$ is the heat value, it has been considered in this paper by $9~970$ [kJ/kW]; $PR_c$ is the fossil fuel cost which is considered as $90~000$ [Yen/kL]; $Q_c$ is the heat value of fossil fuel (heat value of fossil fuel has taken by $41.9 \times 10^6$ [kJ/kL]); $\eta_{i}(t)$ is the generating efficiency of the $i$ th DG at time $t$.


    5.2. Constraints


    5.2.1. System Power Balance Limit

    The total amount of outputs of each DG, the charging and discharging power of BESS and the interchanging power should be equal to load demand at every time as follows:

    $ \sum\limits_{i=1}^{N}{P_{DGi}(t)+P_{BESS}(t)+P_{Re}(t)+P_{day}(t)+P_{hour}(t)}=P_{RTPLoad}(t) $ (5.4)

    where, $P_{Re}(t)$ is the output [kW] of RESs at time $t$; $P_{BESS}(t)$ is the charging and discharging power [kW] of the BESS at time $t$; $P_{day}(t)$ and $P_{hour}(t)$ are the trading power [kW] in the day-ahead and hour-ahead markets, respectively. $P_{Load}(t)$ is the load demand [kW] at time $t$.


    5.2.2. Constraints Related to RTP

    The controllable load $P_{RTPLoad}$ after RTP must be within the limit of $P_{Load}^{'min}$ and $P_{Load}^{'max}$ every time as follows:

    $ P_{Load}^{'min}(t)\leq P_{RTPLoad}(t)\leq P_{Load}^{'max}(t) $ (5.5)

    where, $P_{RTPLoad}$ is the modified load demand by RTP [kW]. In addition, $P_{Load}^{'max}$ and $P_{Load}^{'min}$ represent the maximum and minimum load demand [kW], respectively. $P_{Load}^{'min}$ and $P_{Load}^{'max}$ are:

    $ P_{Load}^{'max}(t)=P'_{Load}(t)+\frac{P_{rtp}}{P_{RTPLoad}(t)} $ (5.6)
    $ P_{Load}^{'min}(t)=P'_{Load}(t)-\frac{P_{rtp}}{P_{RTPLoad}(t)} $ (5.7)

    where, $P_{rtp}$ is the response power [kW]; $P'_{Load}$ is the load demand before RTP [kW]. Furthermore, the load consumption after RTP should be equal to the load consumption before RTP as follows:

    $ W_{re}(T)=\left|\int_0^TP_{in}(t)dt-\int_0^TP_{de}(t)dt\right|=0 $ (5.8)

    where, $W_{re}(T)$ is the controlled load consumption [kW]; $P_{in}(t)$ is the increased load power [kWh]; $P_{de}(t)$ is the decreased load power [kW] as follows:

    $ P_{in}(t)=P_{RTPLoad}(t)-P'_{Load}, ~~~(P_{RTPLoad}\geq{P'_{Load}}) $ (5.9)
    $ P_{de}(t)=P'_{Load}(t)-P_{RTPLoad}, ~~~(P'_{Load}\geq{P_{RTPLoad}}) $ (5.10)

    5.2.3. Constraints Related to DG

    The generation output of each DG must be limited within a specified range as follows:

    $ P_{DGi}(t) \leq P_{DGi}^{max} $ (5.11)

    where, $P_{DGi}^{max}$ and the maximum output power [MW] of $DG_{i}$. Furthermore, the spinning reserve must be installed in each DG to cope with load demand fluctuation as follows:

    $ P_{DGi}(t)+P_{re i} \leq P_{DGi}^{max} $ (5.12)

    where, $P_{re i}$ is the spinning resereve [kW] of the $i$ th DG which is considered as 10[%].


    5.2.4. Constraints Related to BESS

    The charging and discharging power of BESS need to be within $P_{BESS}^{min}$ and $P_{BESS}^{max}$ as follows:

    $ P_{BESS}^{min}\leq P_{BESS}(t) \leq P_{BESS}^{max} $ (5.13)

    where, $P_{BESS}^{min}$ and $P_{BESS}^{max}$ are the maximum charging power [kW] and maximum discharging power of the BESS. Furthermore, the BESSs can be charged and discharged within the limit ranges of $SOC_{min}$ and $SOC_{max}$ as follows:

    $ SOC_{min} \leq SOC(t) \leq SOC_{max} $ (5.14)

    where, $SOC_{min}$ and $SOC_{max}$ are the minimum and maximum SOC [%] for BESS at time $t$, it is considered by 20% and 80% to cope with the output fluctuated of RESs.


    5.3. Optimization Method

    In this paper, the unit commitment for the DG, the charging and discharging plan for the BESS, and the trading power are determined by using this optimization method, Tabu-search. The optimization flowchart is shown in Fig. 7.

    Figure 7. Flow chart of proposed method.

    STEP1: The forecasting data of the load demand and the power generated by PV and WG are inputted.

    STEP2: The unit commitment plans of DGs are determined by the forecasting data in each power system, and then determine the bidding power to the bidding price in day-ahead market by calculating the shortage and surplus power from the determined plan.

    STEP3: The contract processing is conducted from the bidding price and power in day-ahead market, and the electricity is traded by determining the contract price and power. In addition, the operational plan is created considering the day-ahead market.

    STEP4: The reforecasting and re-planning are conducted, and calculate the shortage and surplus power caused by the forecasting error of the power supply and demand. Furthermore, the bidding power to the bidding price in hour-ahead market is determined by calculating the shortage and surplus power from the determined plan.

    STEP5: The operational plan is determined considering the hour-ahead market. Furthermore, if forecasting data is updated, the algorithm proceeds to STEP3, or it ends.


    6. Results and Discussion

    To confirm the effectiveness of the proposed method, this Section makes a comparison between three cases: the operational plan is determined without considering two kinds of power markets and the RTP in Case 1. Case 2 adresses the day-ahead market and RTP, and Case 3 looked at proposed method, the operational plan is determined by considering the reforecasting and re-planning, the day-ahead and hour-ahead markets, and the RTP. In this paper, the operational plan of one day is determined, and it is assumed that forecasting data is updated every 3 hour. Furthermore, the uncertainty regarding the power supply-demand is considered by applying scenario-based method and using 100 kinds of scenarios with respect to the load demand and the output of RESs. Three cases are operated in each scenario, and the method for confirming the effectiveness is to calculate the expectation value of the operational cost.

    The simulation results of three cases are illustrated in Figs. 8, 9 and 10, respectively. It can be seen from Fig. 8 that the number of starts and stops of the DGs in PS2 are many at morning, while that of in PS3 are many at night. This is because the DGs are required to start and stop to cope with forecasting error of the load demand and power generated by RESs. In particular, the SOC reaches 80% after 16:00 in the PS3 because the surplus power occurred in large amount caused by the forecasting error. Therefore, the DGs start and stop frequently because the BESS cannot charge. From Figs. 9 and 11(a), PS1 is obtained profit by selling electricity in day-ahead market. In PS2 and PS3, they meet the load demands in the forenoon by buying electricity in day-ahead market and charging and discharging of the BESSs. The results show that day-ahead market brings about reduction in the number of starts and stops of the DGs. At night, DG in PS3 unavoidably stopped because BESS cannot charge. Furthermore, in PS2 and PS3 in case 3, the forecasting errors of the load demand and the renewable energies are compensated by interchanging electricity through two kinds of power markets and RTP as shown in Figs. 9 and 11(b). It is clear from Figs. 8, 9 and 10 that the number of starts and stops of the DGs in proposed method is less than that of the DGs in Case 1 and Case 2. Fig. 12 shows the simulation results of RTP. The load demand is decreased when the electricity price is expensive, while it is increased when the electricity price is low as shown in Fig. 12. Thereby, it is accomplished to level the load demands in each PS. The expected values of the Fuel and operational costs in each case are listed in Table 2. It can be seen that the value of the proposed method has the lowest cost of three cases. Thus, the high-efficiency operation has been achieved, even when the forecasting error occur.

    Figure 8. Simulation results of Case 1.
    Figure 9. Simulation results of Case 2.
    Figure 10. Simulation results of Case 3.
    Figure 11. Simulation results for power markets.
    Figure 12. Simulation results of RTP.
    Table 2. Fossil fuel and operational costs.
    Case 1 Case 2 Case 3
    RTP without with with
    Day-ahead market without with with
    Hour-ahead market without without with
    Fuel cost [Yen] 2,810,100 2,168,100 1,738,900
    (100%) (77.2%) (61.9%)
    Operational cost [Yen] 2,937,047 2,294,047 1,865,847
    (100%) (78.1%) (63.5%)
     | Show Table
    DownLoad: CSV

    7. Conclusion

    In recent year, RES and BESS have been introducing into the power systems in isolated small islands in Japan because of its low economic cost as compared to DGs. However, the problem of uncertainty of renewable energies has a bad effect on the power systems, especially in the case of the small systems. therefore, as countermeasures towards the above problem, the process of creating the effective operational plan for RTP, power trading at the markets, and the unit commitment of DGs has been explained in this paper. Furthermore, the uncertainty of renewable energies is faithfully modeled in simulation by taking into account that forecast errors are increased in proportion to the time.

    This paper focuses on determining the optimal operation plan capable of coping with uncertainty regarding the power supply-demand by considering the reforecasting and re-planning, the day-ahead and hour-ahead markets, and RTP. By applying RTP, it was able to create the operational plan which can cope with the output variation of RESs by using the DR ability of customers, and to effectively use energies. The power markets enable to transact the shortage and surplus power caused by the output variations of RESs with neighborhood islands. Furthermore, the optimal amount of trading power were determined by fuzzy inference in this paper. The proposed method utilized renewable energies effectively by suppressing the surplus and shortage powers caused by a forecasting error. Fortunately the proposed method was able to reduce the total operational cost of the multi-area small power systems even if the forecasting errors relating the power supply-demand occurred.

    In this work, the optimum installed capacity of DGs, RESs and BESSs are not considered. However, optimum installed capacity is required to be considered for coping with uncertainty regarding the power supply-demand with increase in RESs. Thus, this study can be extended into the optimal operation method considering the optimum installed capacity of each system for the future research.


    Conflict of Interest

    All authors declare on conflicts of interst in this paper.


    [1] Takeuchi H, Ha D, King T-J (2004) Observation of bulk HfO2 defects by spectroscopic ellipsometry. J Vac Sci Technol A 22: 1337-1341. doi: 10.1116/1.1705593
    [2] Nguyen NV, Davydov AV, Chandler-Horowitz D, et al. (2005) Sub-bandgap defect states in polycrystalline hafnium oxide and their suppression by admixture of silicon. Appl Phys Lett 87: 192903. doi: 10.1063/1.2126136
    [3] Price J, Lysaght PS, Song SC, et al. (2007) Identification of sub-band-gap absorption features at the HfO2/Si(100) interface via spectroscopic ellipsometry. Appl Phys Lett 91: 061925. doi: 10.1063/1.2769389
    [4] Ferrieu F, Dabertrand K, Lhostis S, et al. (2007) Observation of HfO2 thin films by deep UV spectroscopic ellipsometry. J Non-Cryst Solids 353: 658-662. doi: 10.1016/j.jnoncrysol.2006.10.032
    [5] Price J, Lysaght PS, Song SC, et al. (2008) Observation of interfacial electrostatic field-induced changes in the silicon dielectric function using spectroscopic ellipsometry. Phys. Status Solidi A 205: 918-921. doi: 10.1002/pssa.200777836
    [6] Price J, Bersuker G, Lysaght PS (2009) Identification of interfacial defects in high-κ gate stack films by spectroscopic ellipsometry. J Vac Sci Technol B 27: 310-312. doi: 10.1116/1.3021045
    [7] Price J, Bersuker G, Lysaght PS (2012) Identification of electrically active defects in thin dielectric films by spectroscopic ellipsometry. J Appl Phys 111: 043507. doi: 10.1063/1.3684605
    [8] Vasić R, Consiglio S, Clark RD, et al. (2013) Multi-technique x-ray and optical characterization of crystalline phase, texture, and electronic structure of atomic layer deposited Hf1-xZrxO2 gate dielectrics deposited by a cyclical deposition and annealing scheme. J Appl Phys 113: 234101. doi: 10.1063/1.4811446
    [9] Fan X, Liu H, Zhang X (2014) Identification of optimal ALD process conditions of Nd2O3 on Si by spectroscopic ellipsometry. Appl Phys A 114: 545-550.
    [10] Fan X, Liu H, Zhang X, et al. (2015) Optical characteristics of H2O-based and O3-based HfO2 films deposited by ALD using spectroscopy ellipsometry. Appl Phys A 119: 957-963. doi: 10.1007/s00339-015-9048-9
    [11] Collins RW, Ferlauto AS (2005) Optical physics of materials, In: Tompkins HG, Irene EA (Eds.), Handbook of Ellipsometry, Norwich: William Andrew Publishing/ Noyes, 93-235.
    [12] Jellison GE Jr. (2005) Data analysis for spectroscopic ellipsometry, In: Tompkins HG, Irene EA (Eds.), Handbook of Ellipsometry, Norwich: William Andrew Publishing/ Noyes, 237-296.
    [13] Petrik P (2014) Parameterization of the dielectric function of semiconductor nanocrystals. Physica B 453: 2-7. doi: 10.1016/j.physb.2014.03.065
    [14] Cho YJ, Nguyen NV, Richter CA, et al. (2002) Spectroscopic ellipsometry characterization of high-κ dielectric HfO2 thin films and the high-temperature annealing effects on their optical properties. Appl Phys Lett 80: 1249-1251. doi: 10.1063/1.1448384
    [15] Sancho-Parramon J, Modreanu M, Bosch S, et al. (2008) Optical characterization of HfO2 by spectroscopic ellipsometry: Dispersion models and direct data inversion. Thin Solid Films 516: 7990-7995. doi: 10.1016/j.tsf.2008.04.007
    [16] Jellison, Jr GE, Modine FA (1996) Parameterization of the optical functions of amorphous materials in the interband region. Appl Phys Lett 69: 371-373; Erratum: ‘‘Parameterization of the optical functions of amorphous materials in the interband region'' [Appl. Phys. Lett. 69, 371 (1996)], idid. 69: 2137. doi: 10.1063/1.118064
    [17] Ferlauto AS, Ferreira GM, Pearce JM, et al. (2002) Analytical model for the optical functions of amorphous semiconductors from the near-infrared to ultraviolet: Applications in thin film photovoltaics. J Appl Phys 92: 2424-2436. doi: 10.1063/1.1497462
    [18] Ferlauto AS, Ferreira GM, Pearce JM, et al. (2004) Analytical model for the optical functions of amorphous semiconductors and its applications for thin film solar cells. Thin Solid Films 455-456: 388-392. doi: 10.1016/j.tsf.2003.11.234
    [19] Foldyna M, Postava K, Bouchala J, et al. (2004) Model dielectric functional of amorphous materials including Urbach tail, In: Pistora J, Postava K, Hrabovsky M, et al. (Eds.), Microwave and Optical Technology 2003, Ostrava, Czech Republic, August 11-15, 2003, SPIE Proc. 5445: 301-305.
    [20] Falahatgar SS, Ghodsi FE (2013) A developed model for the determination of the dielectric function for some absorbing thin films using pseudo-Urbach tail. Physica B 412: 4-11. doi: 10.1016/j.physb.2012.12.011
    [21] Price J, Hung PY, Rhoad T, et al. (2004) Spectroscopic ellipsometry characterization of HfxSiyOz films using the Cody-Lorentz parameterized model. Appl Phys Lett 85: 1701-1703. doi: 10.1063/1.1784889
    [22] Kamineni VK, Hilfiker JN, Freeouf JL, et al. (2011) Extension of far UV spectroscopic ellipsometry studies of high-κ dielectric films to 130 nm. Thin Solid Films 519: 2894-2898. doi: 10.1016/j.tsf.2010.12.080
    [23] Mei JJ, Chen H, Shen WZ, et al. (2006) Optical properties and local bonding configurations of hydrogenated amorphous silicon nitride thin films. J Appl Phys 100: 073516. doi: 10.1063/1.2356915
    [24] Eiamchai P, Chindaudom P, Pokaipisit A, et al. (2009) A spectroscopic ellipsometry study of TiO2 thin films prepared by ion-assisted electron-beam evaporation. Curr Appl Phys 9: 707-712. doi: 10.1016/j.cap.2008.06.011
    [25] Avci N, Smet PF, Poelman H, et al. (2009) Characterization of TiO2 powders and thin films prepared by non-aqueous sol-gel techniques. J Sol-Gel Sci Technol 52: 424-431. doi: 10.1007/s10971-009-2028-9
    [26] Peiponen K-E, Vartiainen EM (1991) Kramers-Kronig relations in optical data inversion. Phys Rev B 44: 8301-8303. doi: 10.1103/PhysRevB.44.8301
    [27] De Sousa Meneses D, Malki M, Echegut P (2006) Structure and lattice dynamics of binary lead silicate glasses investigated by infrared spectroscopy. J Non-Cryst Solids 352: 769-776. doi: 10.1016/j.jnoncrysol.2006.02.004
    [28] Jellison, Jr. GE, Modine FA (1983) Optical functions of silicon between 1.7 and 4.7 eV at elevated temperatures. Phys Rev B 27: 7466-7472.
    [29] Lautenschlager P, Garriga M, Viña L, et al. (1987) Temperature dependence of the dielectric function and interband critical points in silicon. Phys Rev B 36: 4821-4830. doi: 10.1103/PhysRevB.36.4821
    [30] Vineis CJ (2005) Complex dielectric function of biaxial tensile strained silicon by spectroscopic ellipsometry. Phys Rev B 71: 245205. doi: 10.1103/PhysRevB.71.245205
    [31] Awazu K, Kawazoe H, Saito Y, et al. (1991) Structural imperfections in silicon dioxide films identified with vacuum ultraviolet optical absorption measurements. Appl Phys Lett 59: 528-530. doi: 10.1063/1.105428
    [32] Awazu K, Kawazoe H, Muta K-i (1991) Optical properties of oxygen-deficient centers in silica glasses fabricated in H2 or vacuum ambient. J Appl Phys 70: 69-74. doi: 10.1063/1.350245
    [33] Terada N, Haga T, Miyata N, et al. (1992) Optical absorption in ultrathin silicon oxide films near the SiO2/Si interface. Phys Rev B 46: 2312-2318. doi: 10.1103/PhysRevB.46.2312
    [34] Aarik J, Mändar H, Kirm M, et al. (2004) Optical characterization of HfO2 thin films grown by atomic layer deposition. Thin Solid Films 466: 41-47. doi: 10.1016/j.tsf.2004.01.110
    [35] Lucovsky G, Zhang Y, Luning J, et al. (2005) Intrinsic band edge traps in nano-crystalline HfO2 gate dielectrics. Microelectron Eng 80: 110-113. doi: 10.1016/j.mee.2005.04.052
    [36] Hoppe EE, Sorbello RS, Aita CR (2007) Near-edge optical absorption behavior of sputter deposited hafnium dioxide. J Appl Phys 101: 123534. doi: 10.1063/1.2750406
    [37] Ferrieu F, Dabertrand K, Lhostis S, et al. (2007) Observation of HfO2 thin films by deep UV spectroscopic ellipsometry. J Non-Cryst Solids 353: 658-662. doi: 10.1016/j.jnoncrysol.2006.10.032
    [38] Martínez FL, Toledano-Luque M, Gandía JJ, et al. (2007) Optical properties and structure of HfO2 thin films grown by high pressure reactive sputtering. J Phys D: Appl Phys 40: 5256-5265. doi: 10.1088/0022-3727/40/17/037
    [39] Hill DH, Bartynski RA, Nguyen NV, et al. (2008) The relationship between local order, long range order, and sub-band-gap defects in hafnium oxide and hafnium silicate films. J Appl Phys 103: 093712. doi: 10.1063/1.2909442
    [40] Park J-W, Lee D-K, Lim D, et al. (2008) Optical properties of thermally annealed hafnium oxide and their correlation with structural change. J Appl Phys 104: 033521. doi: 10.1063/1.2961326
    [41] Bersch E, Di M, Consiglio S, et al. (2010) Complete band offset characterization of the HfO2/SiO2/Si stack using charge corrected x-ray photoelectron spectroscopy. J Appl Phys 107: 043702. doi: 10.1063/1.3284961
    [42] Xu K, Sio H, Kirillov OA, et al. (2013) Band offset determination of atomic-layer-deposited Al2O3 and HfO2 on InP by internal photoemission and spectroscopic ellipsometry. J Appl Phys 113: 024504. doi: 10.1063/1.4774038
    [43] Di M, Bersch E, Diebold AC, et al. (2011) Comparison of methods to determine bandgaps of ultrathin HfO2 films using spectroscopic ellipsometry. J Vac Sci Technol A 29: 041001.
    [44] Franta D, Ohlídal I, Nečas D, et al. (2011) Optical characterization of HfO2 thin films. Thin Solid Films 519: 6085-6091. doi: 10.1016/j.tsf.2011.03.128
    [45] Nguyen NV, Han J-P, Kim JY, et al. (2003) Optical properties of jet-vapor-deposited TiAlO and HfAlO determined by vacuum utraviolet spectroscopic ellipsometry, In: Seiler DG, Diebold AC, Shaffner TJ, et al. (Eds.), Characterization and Metrology for ULSI Technology 2003, Austin, TX, U.S.A., March 24-28, 2003, AIP Conf. Proc. 683: 181-185.
    [46] Nguyen NV, Sayan S, Levin I, et al. (2005) Optical band gaps and composition dependence of hafnium-aluminate thin films grown by atomic layer chemical vapor deposition. J Vac Sci Technol A 23: 1706-1713.
    [47] Wang XF, Li Q, Egerton RF, et al. (2007) Effect of Al addition on the microstructure and electronic structure of HfO2 film. J Appl Phys 101: 013514. doi: 10.1063/1.2405741
    [48] Park TJ, Kim JH, Jang JH, et al. (2010) Reduction of electrical defects in atomic layer deposited HfO2 films by Al doping. Chem Mater 22: 4175-4184. doi: 10.1021/cm100620x
    [49] Sandberg RL, Allred DD, Lunt S, et al. (2004) Optical properties and application of uranium-based thin films for the extreme ultraviolet and soft x-ray region, In: Soufli R, Seely JF (Eds.), Optical Constants of Materials for UV to X-Ray Wavelengths, Denver, CO, U.S.A., August 2, 2004, SPIE Proc. 5538: 107-118.
    [50] Meek TT, von Roedern B (2008) Semiconductor devices fabricated from actinide oxides. Vacuum 83: 226-228. doi: 10.1016/j.vacuum.2008.04.005
    [51] Kruschwitz CA, Mukhopadhyay S, Schwellenbach D, et al. (2014) Semiconductor neutron detectors using depleted uranium oxide, In: Burger A, Franks L, James RB, et al. (Eds.), Hard X-Ray, Gamma-Ray, and Neutron Detector Physics XVI, San Diego, CA, U.S.A., August 17, 2014, SPIE Proc. 9213: 92130C-1-92130C-9.
    [52] Chen Q, Lai X, Bai B, et al. (2010) Structural characterization and optical properties of UO2 thin films by magnetron sputtering. Appl Surf Sci 256: 3047-3050. doi: 10.1016/j.apsusc.2009.11.071
    [53] He H, Andersson DA, Allred DD, et al. (2013) Determination of the insulation gap of uranium oxides by spectroscopic ellipsometry and density functional theory. J Phys Chem C 117: 16540-16551. doi: 10.1021/jp401149m
    [54] Khilla MA, Rofail NH (1986) Optical absorption edge of uranium trioxide phases: Part I. Radiochim Acta 40: 155-158.
  • Reader Comments
  • © 2015 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(7166) PDF downloads(1056) Cited by(4)

Article outline

Figures and Tables

Figures(4)  /  Tables(4)

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog