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] Fabrizio Clarelli, Roberto Natalini . A pressure model of immune response to mycobacterium tuberculosis infection in several space dimensions. Mathematical Biosciences and Engineering, 2010, 7(2): 277-300. doi: 10.3934/mbe.2010.7.277
    [2] Chang Gong, Jennifer J. Linderman, Denise Kirschner . A population model capturing dynamics of tuberculosis granulomas predicts host infection outcomes. Mathematical Biosciences and Engineering, 2015, 12(3): 625-642. doi: 10.3934/mbe.2015.12.625
    [3] Eduardo Ibargüen-Mondragón, Lourdes Esteva, Edith Mariela Burbano-Rosero . Mathematical model for the growth of Mycobacterium tuberculosis in the granuloma. Mathematical Biosciences and Engineering, 2018, 15(2): 407-428. doi: 10.3934/mbe.2018018
    [4] Gesham Magombedze, Winston Garira, Eddie Mwenje . Modelling the human immune response mechanisms to mycobacterium tuberculosis infection in the lungs. Mathematical Biosciences and Engineering, 2006, 3(4): 661-682. doi: 10.3934/mbe.2006.3.661
    [5] Eduardo Ibargüen-Mondragón, M. Victoria Otero-Espinar, Miller Cerón Gómez . A within-host model on the interaction dynamics between innate immune cells and Mycobacterium tuberculosis. Mathematical Biosciences and Engineering, 2025, 22(3): 511-527. doi: 10.3934/mbe.2025019
    [6] Eduardo Ibarguen-Mondragon, Lourdes Esteva, Leslie Chávez-Galán . A mathematical model for cellular immunology of tuberculosis. Mathematical Biosciences and Engineering, 2011, 8(4): 973-986. doi: 10.3934/mbe.2011.8.973
    [7] Xu Zhang, Dongdong Chen, Wenmin Yang, JianhongWu . Identifying candidate diagnostic markers for tuberculosis: A critical role of co-expression and pathway analysis. Mathematical Biosciences and Engineering, 2019, 16(2): 541-552. doi: 10.3934/mbe.2019026
    [8] Benjamin H. Singer, Denise E. Kirschner . Influence of backward bifurcation on interpretation of $R_0$ in a model of epidemic tuberculosis with reinfection. Mathematical Biosciences and Engineering, 2004, 1(1): 81-93. doi: 10.3934/mbe.2004.1.81
    [9] Maoxing Liu, Jie Zhang, Zhengguang Li, Yongzheng Sun . Modeling epidemic in metapopulation networks with heterogeneous diffusion rates. Mathematical Biosciences and Engineering, 2019, 16(6): 7085-7097. doi: 10.3934/mbe.2019355
    [10] Abba B. Gumel, Baojun Song . Existence of multiple-stable equilibria for a multi-drug-resistant model of mycobacterium tuberculosis. Mathematical Biosciences and Engineering, 2008, 5(3): 437-455. doi: 10.3934/mbe.2008.5.437
  • 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.


    In December 2019, a novel coronavirus (COVID-19) was discovered [1]. COVID-19 is a respiratory infectious disease caused by the severe acute respiratory syndrome coronavirus, transmitted by contact, aerosols and inhalation of virus-infected droplets [2,3]. Unfortunately, there were no effective drugs to treat the disease in 2021 [1]. Therefore, the control protocols were mainly physical isolation, such as quarantine, contact tracing and social lockdown [4,5,6,7,8]. Therefore, it is crucial to study the effectiveness of control measures for this disease.

    The basic reproduction number $ \mathcal R_0 $ is a critical parameter in the analysis of infectious diseases. It measures the average number of secondary infections caused by a typical infectious individual in a fully susceptible population. The disease will not cause an epidemic if $ \mathcal R_{0} < 1 $ [11,12,13]. As the epidemic progresses, or when control measures such as social distancing and vaccination are implemented, susceptible individuals are depleted and some of the contacts of the infectious individuals are made to already infected individuals, which do not cause transmission. In this case, the threshold for disease spread is measured by the effective reproduction number $ \mathcal R_t $, which measures the number of secondary infections caused by a typical infectious individual in the current population (with some already infected individuals). Its value is usually calculated as the product of $ \mathcal R_0 $ and the average population susceptibility [2,4,11]. In the initial stage of disease when the number of infected individuals is only a tiny fraction of the total population, $ \mathcal R_t\approx \mathcal R_0 $. Wallinga and Teunis [9] proposed a method to approximate the effective reproduction number, which studied the average number of reproduction for patients who are infected on a given day. Note that this is not the same as the number of infections caused by a patient on a given day. In addition, this method requires the time of infection of patients, which is difficult to trace. Cori et al. [14] presented a simple method for estimating the effective reproduction number, which is based on the time series of disease occurrence. However, this method's effective reproduction number $ \mathcal R_t $ is delayed significantly. Das [4] presented a method for estimating an approximate $ \mathcal R_t $ by taking into account both the mean generation interval and the instantaneous exponential growth rate. However, the instantaneous exponential growth rate can only be estimated from a long enough time period, preventing this method from detecting sudden changes in $ \mathcal R_t $. In this paper, we propose a new method to directly calculate the real-time reproduction number through confirmed cases and evaluate the effectiveness of control measures.

    Based on reported cases of COVID-19, we use the back-calculating method [15] to obtain the number of incidences and the number of infectious patients on each day and then we use them to estimate the change in reproduction number. Moreover, we derive the impact of control measures from the change in the reproduction numbers.

    We establish the model in Section 2 and verify it using simulations in Section 3. In Section 4, we show how we applied the model to British Columbia (BC), Canada to obtain their reproduction number and the impact of control measures on the reproduction number; the results are summarized and future work is discussed in Section 5.

    In this section, we consider a discrete-time stochastic seir model in a randomly mixed population. Let $ S_t $, $ E_t $ and $ I_t $ be the number of susceptible, latent and infectious individuals on day $ t $. Because of the random-mixing assumption, the expected number of new infections on day $ t $ is

    $ Zt=βStIt.
    $
    (1)

    A newly infected patient goes through a latent period $ L $ and becomes infectious. Here $ L $ is a discrete random variable with a probability mass function $ \{p_i\}_{i = 0}^\infty $, i.e., the probability that the latent period has a length of $ i $ days is $ p_i $. This patient then goes through an infectious period $ X $ and is diagnosed. Let the probability mass function of $ X $ be $ \{q_i\}_{i = 0}^\infty $. Let $ Q_i $ be the cumulative probability function of $ X $, that is $ Q_i = \sum\limits_{j = 0}^i q_j $, $ 0 \leq i \leq \infty $. We assume that, once diagnosed, the patient is fully isolated and stops being infectious. The course of disease $ Y $ is the sum of the latent and infectious periods, i.e., $ Y = L+X $. Let $ d_i $ be the probability mass function of $ Y $. Then

    $ di=ij=0pjqij.
    $
    (2)

    By definition, the basic reproduction number, i.e., the average number of secondary infections caused by a typical infectious individual during the infectious period in a fully susceptible population, is

    $ \mathcal R_{0} = \beta E[X] , $

    where $ E[X] = \sum\limits_{i = 0}^\infty iq_i $ is the mean infectious period.

    Thus, the effective reproduction number on day $ t $ is

    $ Rt=R0St=βStE[X]=ZtItE[X].
    $
    (3)

    Note that the last step is from (1).

    Given the mean infectious period $ E[X] $, we need to estimate the number of new infections $ Z_{t} $ and the number of infectious individuals $ I_t $ so that we can estimate the effective reproduction number $ \mathcal R_t $.

    Note that $ Z_t $ and $ m_t $ have the following relationship

    $ mt=i=0Ztidi.
    $
    (4)

    With the $ m_t $ given, we need to solve $ Z_t $. Unfortunately, this is a deconvolution problem, and it is difficult to solve (see, e.g., [16]). Instead of solving it, we use the following method to approximate $ Z_t $. Suppose that the number of diagnosed cases on day $ t $, namely $ m_t $, is observed for days $ t = 0, 1, \dots, T $. A patient who is diagnosed on day $ t+i $, $ 0 \leq i \leq T-t $ was infected on day $ t $ if and only if the serial interval $ Y = i $, i.e.,

    $ Zt=Tti=0mt+idi=Tti=0ij=0mt+ipjqij.
    $
    (5)

    Note that (5) does not solve the deconvolution problem given by (4). However, we will show that this can give a good approximation for $ Z_t $, especially if $ m_t $ is approximately exponentially growing or decaying and the change in the exponential growth rate is slow (measured on the time-scale of the mean serial interval). To see this, assume that $ Z_t = Z_0\mu^t $ for a constant $ \mu > 0 $ (the exponential growth rate is thus $ \log\mu $); then, (4) becomes

    $ mt=i=0Z0μtidi=Z0μtg(1/μ),
    $
    (6)

    where $ g(x) = \sum_{i = 0}^\infty x^id_i $ is the probability generating function of the serial interval distribution $ d_i $. Substitute this into the right hand side of (5) and assume $ T\gg1 $; then,

    $ i=0mt+idi=Z0μt+ig(1/μ)di=Z0μtg(μ)g(1/μ)=Ztg(μ)g(1/μ).
    $
    (7)

    Thus, if the change in $ Z_t $ is slow, i.e., $ \mu\approx 1 $, then $ g(\mu)g(1/\mu)\approx 1 $.

    Using a similar approximation, a patient who is diagnosed on day $ t+i $ was infectious on day $ t $ because the infectious period $ X\geq i $. That is,

    $ It=Tti=0mt+iProb{Xi}=Tti=0mt+i(1Qi1).
    $
    (8)

    Thus, from (3),

    $ Rt=Tti=0ij=0mt+ipjqijTti=0mt+i(1Qi1)E[X].
    $
    (9)

    The symbols used in this article are described in detail in Table 1.

    Table 1.  Explanations of nouns appearing in the paper.
    Symbols Significance
    $ S_t $ susceptible at time $ t $
    $ E_t $ latent at time $ t $
    $ I_t $ infectious at time $ t $
    $ C_t $ the number of individuals who are infected on day $ t $
    $ L $ latent period
    $ X $ infectious period
    $ Y $ the sum of the latent and infectious periods
    $ q_i $ the probability that the infectious period has a length of $ i $ days
    $ p_i $ the probability that the latent period has a length of $ i $ days
    $ d_i $ the probability mass function of the disease course $ Y $
    $ m_t $ the number of cases reported on day $ t $

     | Show Table
    DownLoad: CSV

    Given the mean infectious period E[X], in order to find the $ 95% $ confidence interval of $ \mathcal R_t $, we will get a random sample of $ Z_{t} $ and $ I_{t} $ by using the Monte Carlo method. For the patients who are diagnosed on day $ t $, let $ \tilde Z_{t, t-i} $ be the number of those who are infected on day $ t-i $ for $ i = 0, 1, \dots $. Then, the approximation given by (5) is equivalent to the following two steps:

    $ {\rm{a}}) $ Assume that $ \tilde Z_{t, t-i} $ is multinomially distributed according to

    $ ˜Zt,timultinomial(mt,di).
    $
    (10)

    $ {\rm{b}}) $ The number of individuals who are infected on day $ t $ is

    $ Zt=Tti=0˜Zt+i,t.
    $
    (11)

    Note that the mean of (11) is given by (5).

    Similarly, for patients who show symptoms on day $ t $, the number of those who were infected on day $ t-i $ is

    $ ˜Ct,timultinomial(mt,qi).
    $
    (12)

    Thus, the number of individuals who are infected on day $ t $ is

    $ Ct=Tti=0˜Ct+i,t.
    $
    (13)

    Note that, here, $ C_t $ is a random variable, the mean of which is given by

    $ E[Ct]=Tii=0mt+iqi,
    $
    (14)

    which uses a similar approximation as (5).

    The expected number of individuals $ I_t $ who are infectious on day $ t $ is the total number of people who have become infectious but have not been removed from transmission (via recovery or isolation). Thus, $ I_t $ can be calculated as

    $ It=ti=0Citi=0mi,
    $
    (15)

    where the first term on the right-hand side is the number of patients who have become infectious before (or on) day $ t $, while the second term is the number of patients who have been diagnosed and isolated before (or on) day $ t $.

    To generate one sample of $ Z_{t} $ and $ I_t $, for each $ m_t $, $ t = 0, 1, \dots, T $, we use (12) to generate a sample for $ \tilde C_{t, t-i} $, and then use (13) to calculate $ C_{t} $. We then use the calculated $ C_{t} $ to generate a sample of $ Z_t $ using (10) and (11), and we use the calculated $ C_{t} $ to generate a sample of $ I_t $ using (10), (11) and (15). We can then use (3) to calculate a sample of the curve $ \mathcal R_t $.

    We generate $ 10^5 $ samples for $ \mathcal R_t $. For each $ t = 0, 1, \dots, T $, we use these samples to estimate the $ 95% $ confidence interval.

    To verify that our model can correctly estimate the reproduction number, we apply (9) to a dataset generated from stochastic simulations by using the method in Section 2.2 to estimate the $ 95% $ confidence interval of the reproduction number.

    Non-pharmacological intervention (NPI) measures reduce the transmission rate [17]. The latent period and infectious period are specific to the disease, and are not affected by NPI measures. Thus, we consider the following two cases.

    Case 1 Seasonal variation in the transmission rate is sometimes approximated by a sinusoidal function [10,18,19]. Here we assume $ \beta $ to be sinusoidal to verify that our method can detect continuous change in $ \beta(t) $. Specifically,

    $ β=0.2[cos(2πt365)+1].
    $
    (16)

    In this case, the infectious period is assumed to be gamma-distributed with a shape parameter of $ 3 $ and a rate parameter of $ 0.2 $; the latent period is assumed to be gamma-distributed with a shape parameter of $ 3 $ and a rate parameter of $ 0.3 $. Note that these choices only serve as a numerical example and are (no comma) not tied to a specific disease.

    Case 2 The transmission rate is assumed to be a step function to simulate a sequence of control measures that cause sudden change in $ \beta $, that is,

    $ β={0.5,t[0,40),0.4,t[40,60),0.3,t[60,80),0.2,t[80,100),0.3,t[100,).
    $
    (17)

    In this case, the infectious period is assumed to be gamma-distributed with a shape parameter of $ 3 $ and a rate parameter of $ 0.25 $; the latent period is assumed to be gamma-distributed with a shape parameter of $ 3 $ and a rate parameter of $ 0.2 $.

    In both cases, the simulated population size is $ 10^6 $. During the simulated time periods, the number of infected individuals is only a small fraction of the population size; thus, $ S(t)\approx 1 $ and the effective reproduction number is approximately the basic reproduction number.

    Figure 1 shows the comparison of the estimated reproduction number $ \mathcal R_t $ as a function of time with the true value for Case 1. Figure 2 shows the comparison for Case 2. In both cases, our method can correctly estimate the reproduction number. In addition, these figures also show that the confidence interval narrows with a larger case count.

    Figure 1.  Model validation for Case 1. The top panel shows the simulated epidemic curve (in $ \log_{10} $), while the lower panel shows the comparison of the estimated reproduction number and its confidence interval with the true value in blue. Note that, when the case count is low, the confidence interval becomes very large, causing the upper and lower bounds of the confidence interval to disappear on some days.
    Figure 2.  Model validation for Case 2. The top panel shows the simulated epidemic curve (in $ \log_{10} $), while the lower panel shows the comparison of the estimated reproduction number and its confidence interval with the true value in blue.

    Now that we have validated our method, in this section, we apply the method to study the change of the reproduction number as a function of time for the COVID-19 outbreak in BC, Canada in 2020.

    In BC, we consider the following policy changes:

    ● Provincial state of emergency was declared on March 17;

    ● Businesses reopened on May 19;

    ● Provincial state of emergency was declared on July 7;

    ● Provincial state of emergency was declared on August 5;

    ● Public K-12 schools reopened on September 10;

    ● Provincial state of emergency was declared on October 28;

    ● Provincial state of emergency was declared on December 23.

    We used the daily reported case count data for the period of March 1 to December 31, 2020 that were released from the BC Centre for Disease Control (BCCDC) as a spreadsheet. This spreadsheet has been taken offline. However, the data can still be accessed via the COVID-19 dashboard on the BCCDC website [20].

    To apply our method, we need to know the latent period distribution and the infectious period distribution. We use the latent period distribution estimated by [21], which is gamma-distributed with a mean of 5.48 days and a standard deviation of 2.72 days.

    On the other hand, regional differences in testing policy and human behavior in voluntary testing may affect when a patient is diagnosed and isolated, and, in turn, affect the infectious period. In Subsection 4.1, we estimate the infectious period distribution in BC from the daily number of diagnosed cases and symptom onsets.

    We assume that the patients will be isolated after being diagnosed. Therefore, the end of their infectious period is marked by diagnosis, not recovery. We assume that the infectious period is gamma-distributed [22,23,24,25,26], with a shape parameter $ \alpha $ and a scale parameter $ \varepsilon $.

    We digitized and tabulated the daily number of symptom onsets for the period of January 15 to June 7, 2020 from the British Columbia COVID-19 Daily Situation Report released on June 9 [27].

    Using an approximation similar to (5), the expected number of patients showing symptoms on day $ t $ can be calculated from the diagnosed cases on day $ t+i $ ($ m_{t+i} $) as

    $ λt=Tti=0mt+iqi.
    $
    (18)

    We assume that $ C_t $ is the observed symptom onset count on day $ t $; it follows a Poisson distribution with the mean $ \lambda_t $, i.e.,

    $ CtPoisson(λt)
    $
    (19)

    We use the Markov chain Monte Carlo method via the R package "R2jag" to estimate the distribution parameters $ \alpha $ and $ \varepsilon $. The prior distributions of the parameters are chosen to follow a uniform distribution with wide intervals:

    $ αU(0,10),εU(0,5)
    $
    (20)

    The results are given in Table 2. Figure 3 shows the point estimate and $ 95% $ confidence interval of the estimated density function of the infectious period distribution.

    Table 2.  The estimated parameters for the gamma-distributed COVID-19 infectious period in BC, Canada in 2020.
    Parameter Mean Sd $ 95 $% confidence interval
    $ \alpha $ 4.791 1.074 (3.133, 7.358)
    $ \varepsilon $ 1.815 0.373 (1.143, 2.610)

     | Show Table
    DownLoad: CSV
    Figure 3.  The point estimates and $ 95% $ confidence interval of the probability density function of the infectious period distribution in BC, Canada.

    Using the point estimate of the infectious period distribution in Subsection 4.1, we estimate the reproduction number as a function of time in BC, Canada.

    Figure 4 shows both the epidemic curve (reported cases) and the estimated reproduction number. This figure shows that, since the provincial state of emergency on March 17, the reproduction number was controlled to below 1 until the relaxation (business reopening) announced on May 19. The reproduction number then increased gradually after the relaxation to 1.76 in June, being largely maintained until August 1st, at which point it was about 1.60. The strengthening of control measures on August 5 reduced the reproduction number and eventually controlled it to around 1 on September 10. It then increased again to a peak value of 1.67 on October 10. It was then brought back to about unity beginning on November 13.

    Figure 4.  The top panel shows the daily reported cases (in $ \log_{10} $). The lower panel shows the estimated reproduction number and the confidence interval on each day. The vertical lines show the dates of the implementation of epidemic control measures, where the orange lines show the declaration of a state of emergency, and the green lines show the dates of relaxation.

    We have developed a novel method to estimate the change of the reproduction number with time. Using simulated data, we have shown that our method can estimate the change in the reproduction number due to either seasonal forcing or control measures. This means that our method is widely applicable to understand the change of the transmission rate.

    Applying our method to the COVID-19 outbreak in BC, Canada in 2020 shows that the strengthening of control measures such as social distancing, restricting gathering and closing schools from March 20 to the end of May successfully reduced the reproduction number to below 1, except for a period in early April (may be due to clustered cases in long-term care facilities [28], or the gathering activities of the Easter holiday). However, the reproduction number gradually increased to above 1 after business reopening in May, even though the case counts did not exhibit an immediate increase. This shows that our method is very sensitive as a tool to detect the changes in reproduction number. The same increase also happened after the school reopening in September, which eventually triggered the fast increase of cases in October and early November. Note that, during this time, the variants of concern had not shown up yet, that is, no variants of concern appeared in 2020 [29,30]. Thus, the increase of cases is mostly likely due to the relaxation of control measures.

    Not surprisingly, our estimation yields a narrower confidence interval with a larger case count. Our method also relies on reliable estimation of the latent and infectious periods, which may be difficult to estimate during the early stage of a disease outbreak. However, we have also demonstrated that our method can be adapted to estimate the infectious period from the daily counts of symptom onset and diagnosed cases.

    Another limitation of our method is that it ignores asymptomatic and pre-symptomatic transmissions, which may be an important factor driving the COVID-19 transmission. However, this may not significantly affect our method if the ratio of asymptomatic cases to all cases remains roughly constant, as the proportional fact is canceled in our formulation.

    Our method provides a new tool for analyzing host immunity resulting from the effective vaccinations, with and without NPI measures. At the onset of disease spread, the effective vaccination rate $ v $ is the primary factor influencing the number of susceptible individuals. In this case, the effective reproduction number can be expressed as $ R_t = \beta (1-v)S_tE(X) $, which is similar to (3). If the effective vaccination rate $ v $ is known, our method can estimate the temporal changes in $ \beta $. However, in the absence of information about $ v $, it is only possible to estimate the value of $ \beta(1-v) $, not the individual parameters $ \beta $ and $ v $.

    Our method can be used to study other infectious diseases as well. For instance, it can be used to investigate the influence of seasonality on the transmission of seasonal influenza, or to examine the effect of control measures on historical outbreaks, such as pandemic influenza, SARS and Ebola. Furthermore, our method can be applied to the study of vector-borne diseases, including those transmitted by mosquitoes, by extending our model to consider the disease transmission from person to person, with mosquitoes as the vectors. However, obtaining the specific changes in $ \beta $ is challenging, as the infection rate through the vector depends on the change in infected mosquitoes, resulting in a more complex dependence of $ \beta $ on mosquitoes than the simple $ SEIR $ model. Therefore, further research is needed to address this complexity. Additionally, the same generalization can be applied to sexually transmitted infections.

    The authors declare that they have not used artificial intelligence tools in the creation of this article.

    This research was supported by the National Natural Science Foundation of China (No. 12271088) (ML), the Natural Sciences Foundation of Shanghai (No. 21ZR1401000) (ML) and a discovery grant from the Natural Sciences and Engineering Research Council Canada (JM), as well as two NSERC EIDM grants (OMNI and MfPH) (JM).

    The authors declare that there is no conflict of interest.

    [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.
  • This article has been cited by:

    1. Boxiang Yu, Slim Ibrahim, Junling Ma, Bin Sun, Meili Li, Disease threshold of the effective degree SIS model on a random network with large degrees, 2025, 426, 00220396, 206, 10.1016/j.jde.2025.01.077
  • 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(6708) PDF downloads(1054) Cited by(4)

Figures and Tables

Figures(4)  /  Tables(4)

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog