Processing math: 100%
Review Topical Sections

Overview of thin film deposition techniques

  • Received: 06 September 2018 Accepted: 23 January 2019 Published: 13 March 2019
  • Surface properties of the material can affect the efficiency and behavior of the material when in service. Modifying and tuning these surface properties to meet the specific demand for better performance is feasible and has been vastly employed in a different aspect of life. This can be achieved by coating the surface via deposition of the thin film. This study provides a review of the existing literature of different deposition techniques used for surface modification and coating. The two major areas of interest discussed are physical and chemical vapor deposition techniques, and the area of applications of surface coating was briefly highlighted in this report.

    Citation: Olayinka Oluwatosin Abegunde, Esther Titilayo Akinlabi, Oluseyi Philip Oladijo, Stephen Akinlabi, Albert Uchenna Ude. Overview of thin film deposition techniques[J]. AIMS Materials Science, 2019, 6(2): 174-199. doi: 10.3934/matersci.2019.2.174

    Related Papers:

    [1] Shu Xu, Yichang Wei, Abdul Hafeez Laghari, Xianming Yang, Tongchao Wang . Modelling effect of different irrigation methods on spring maize yield, water and nitrogen use efficiencies in the North China Plain. Mathematical Biosciences and Engineering, 2021, 18(6): 9651-9668. doi: 10.3934/mbe.2021472
    [2] Zixiao Xiong, Xining Li, Ming Ye, Qimin Zhang . Finite-time stability and optimal control of an impulsive stochastic reaction-diffusion vegetation-water system driven by Lˊevy process with time-varying delay. Mathematical Biosciences and Engineering, 2021, 18(6): 8462-8498. doi: 10.3934/mbe.2021419
    [3] Chenyang Zhan, Min Tian, Yang Liu, Jie Zhou, Xiang Yi . A novel greedy adaptive ant colony algorithm for shortest path of irrigation groups. Mathematical Biosciences and Engineering, 2022, 19(9): 9018-9038. doi: 10.3934/mbe.2022419
    [4] Alejandro Rincón, Fredy E. Hoyos, John E. Candelo-Becerra . Comparison, validation and improvement of empirical soil moisture models for conditions in Colombia. Mathematical Biosciences and Engineering, 2023, 20(10): 17747-17782. doi: 10.3934/mbe.2023789
    [5] Favour Adenugba, Sanjay Misra, Rytis Maskeliūnas, Robertas Damaševičius, Egidijus Kazanavičius . Smart irrigation system for environmental sustainability in Africa: An Internet of Everything (IoE) approach. Mathematical Biosciences and Engineering, 2019, 16(5): 5490-5503. doi: 10.3934/mbe.2019273
    [6] Todorka Glushkova, Stanimir Stoyanov, Lyubka Doukovska, Jordan Todorov, Ivan Stoyanov . Modeling of an irrigation system in a virtual physical space. Mathematical Biosciences and Engineering, 2021, 18(5): 6841-6856. doi: 10.3934/mbe.2021340
    [7] Weitao Liu, Qiushuang Zheng, Lifu Pang, Weimeng Dou, Xiangxi Meng . Study of roof water inrush forecasting based on EM-FAHP two-factor model. Mathematical Biosciences and Engineering, 2021, 18(5): 4987-5005. doi: 10.3934/mbe.2021254
    [8] Xiaomei Bao, Canrong Tian . Turing patterns in a networked vegetation model. Mathematical Biosciences and Engineering, 2024, 21(11): 7601-7620. doi: 10.3934/mbe.2024334
    [9] B. W. Kooi, D. Bontje, M. Liebig . Model analysis of a simple aquatic ecosystems with sublethal toxic effects. Mathematical Biosciences and Engineering, 2008, 5(4): 771-787. doi: 10.3934/mbe.2008.5.771
    [10] Vincenzo Luongo, Maria Rosaria Mattei, Luigi Frunzo, Berardino D'Acunto, Kunal Gupta, Shankararaman Chellam, Nick Cogan . A transient biological fouling model for constant flux microfiltration. Mathematical Biosciences and Engineering, 2023, 20(1): 1274-1296. doi: 10.3934/mbe.2023058
  • Surface properties of the material can affect the efficiency and behavior of the material when in service. Modifying and tuning these surface properties to meet the specific demand for better performance is feasible and has been vastly employed in a different aspect of life. This can be achieved by coating the surface via deposition of the thin film. This study provides a review of the existing literature of different deposition techniques used for surface modification and coating. The two major areas of interest discussed are physical and chemical vapor deposition techniques, and the area of applications of surface coating was briefly highlighted in this report.


    The increasing temperatures resulting from global climate change are affecting the fragile equilibrium of ecosystems worldwide [1]. Mediterranean regions are affected by heat waves and decreased rainfall, which affects the water supply and food production systems [2]. In Mediterranean and semi-arid climate areas, artificial irrigation is necessary to maintain global agriculture production during the growing season (spring to summer). The main source of irrigation water is either reservoirs that usually fill with rain or snowmelt during the fall and winter seasons. In recent years, the lack of precipitation has substantially decreased water availability for seasonal irrigation. Dry soils and low water availability reduce the plants' evapotranspiration fluxes, which generates a physiological response, called water stress, due to water scarcity [3]. When plants experience water stress, a physiological response is induced that reduces their photosynthetic rates. This negatively impacts the growth of the plant [4,5,6,7], compromising the balance between water availability in the soil and the crop. Farmers were forced to adjust irrigation practices to improve irrigation water productivity.

    Theoretically, the implementation of an irrigation strategy should focus on the quantity and timing of irrigation to meet two criteria: (Ⅰ) supplying the plants' water needs of evapotranspiration and (Ⅱ) keeping the available soil water content within a certain range [8]. However, in the field, irrigation decisions are generally empirical by observing signs of water stress in plants [9].

    There are different technical recommendations for establishing an irrigation strategy to reduce total water consumption during the growing season without negatively affecting the fruit yield and quality.

    Among the different technical recommendations, one highly recommended strategy is regulated deficit irrigation (RDI). RDI consists of irrigating only when the soil water level is minimal or until the physiological response of the plant is not affected, producing controlled levels of water stress [10]. Several field studies suggest their implementation to optimize water use without affecting the soil-plant dynamics [11,12,13,14]. However, due to the empirical nature of RDI, one needs to set-up the experimental plots before their application in the entire field, requiring time and costs.

    The construction of mathematical models facilitates our understanding of agricultural phenomena and offers advantages for simulating different scenarios of dryness by considering the complexes involved in determining how dryness affects the plant's production. These models can be used to analyze the possible effects of implementing RDI strategies, thereby providing the theoretical behavior of plant-water dynamics in the soil. Mathematical modeling offers an excellent alternative to simulate the interaction between the plant and water in the soil. Mathematical models applied to agricultural systems incorporate mass balance equations, and cause-effect relationships between variables, among others [15]. Recently, some models for the growth of crops under water deficit have been presented [4,16], considering advantages such as simulating multiple scenarios without affecting live systems and reducing costs associated with experimental field studies.

    Indicators are numerical variables that provide simple and relevant information [17], which can be environmental, physical, physiological, or biological. In agricultural systems, indicators have been used for decision-making concerning particular phenomena or conditions subject to a specific value or range of values. However, the construction of these indicators is subject to multiple variables that contemplate the dynamics of the phenomenon studied.

    For example, physiological indicators such as the stem water potential (Ψ), stomatal conductance (gs), and the Crop Water Stress Index (CWSI), among others, have been used to monitor the water status of plants and to determine when to irrigate [18,19,20,21,22]. These indicators allow the producers to verify whether or not the plants are suffering from water stress [23,24]. However, measuring these indicators in the field can be limited by their cost, time cunsumption, and representativeness, mainly because they are obtained from sentinel plants.

    Based on those mentioned above, the construction of a non-physiological indicator to decide the moment of irrigation could be an important proposal for simulation, thereby reflecting the interaction among the soil water content and the plant water behavior. Indicators are widely used in mathematical modeling, providing a theoretical basis for decision-making, considering simulated situations. The next generation matrix method developed by [25] is a mathematical method used in dynamic and biological systems to calculate the local stability of a system of ordinary differential equations. This method has been extended and developed in epidemiology by [26] to determine a threshold parameter known as the basic reproduction number. This parameter is used to establish either the persistence or absence of a disease, considering compartmentalized mathematical models [26]. In biology, this kind of model has been applied in several ways, such as calculating the source-sink dynamics in marine meta populations [27] and the evolutionary invasion analysis [28]. Additionally, the next-generation method has been applied in agricultural systems to study the behavior of plant diseases [29,30,31,32,33,34].

    This work proposes the development of a soil water indicator based on the next-generation matrix method to determine the minimum soil water amount required for optimal plant growth and fruit production. As far as the author's know, while the next-generation matrix method has been mainly used in biotic systems, its application to abiotic phenomena such as the water fluxes dynamic between soil and plants has not been explored. It is hypothesized that this method could provide an effective theoretical indicator for irrigation decisions.

    This work presents the construction of the indicator based on a theoretical mathematical model that considers two irrigation scenarios: a constant rate of irrigation without water limitations and an impulsive irrigation strategy with cycles of irrigation and water cut.

    The construction of a mathematical model that describes the soil-crop dynamics implies identifying relevant variables and parameters, such as genotype, soil water availability, energy, soil type, crop management conditions, fertilization, depth roots, and environmental factors, among other variables [4,35,36,37,38].

    To approach the mathematical model describing soil moisture-crop dynamics carried out through the relationship between the amount of water in the soil and the growth of the fruit tree, we consider the following assumptions presented in [4,39]:

    ⅰ. Crop growth dynamics are influenced by the interaction between energy, water, and vegetative growth variables.

    ⅱ. The crops respond immediately to irrigation application.

    ⅲ. Adult crop and a suitable soil for crop growth.

    ⅳ. Optimal agronomic management conditions, including pest control and fertilizer management.

    ⅴ. The system´s energy is constant.

    ⅵ. The system only considers the application of irrigation as a source of water; other possible sources, such as rainfall or groundwater, as not been considered in the model construction.

    Table 1 provides a summary of the parameters that affect the growth dynamics of plants that produce fruits, according to our model.

    Table 1.  Parameters used in the plants that produce fruits growth model.
    Parameter Meaning
    q Accumulated energy constant
    r Intrinsic growth rate of plants that produce fruits
    N Fruit tree carrying capacity
    p Rainfall contribution rate
    I Irrigation amount
    β Evapotranspiration rate
    γ Photosynthetic contribution rate
    ω Mortality rate of plants that produce fruits

     | Show Table
    DownLoad: CSV

    We denote W=W(t) as the water amount in the soil in time t and C=C(t) as the concentration of biomass in the plant that produces fruits at time t, both state variables. The variation of the amount of water in the soil concerning to time denoted W(t) is determined by the relationship between environmental and physiological variables. That is, there is a loss due to the energy (q) - water (W) interaction (Evapotranspiration) at a β rate. In addition, an alternative amount of water is used for the growth of the plant that produces fruit at a rate of r. Additionally, the variation of the biomass concentration concerning to time, denoted C(t), is positively affected by the contribution of water to the growth of the plant that produces fruits at the rate of r. An entry is presented in the variation of the biomass concentration due to the photosynthetic contribution to the growth of the plant given by the relation energy (q) - water W and vegetative growth C at a rate γ and is considered an output given by the term ωC, which represents a loss of biomass concentration due to natural death at a rate of ω.

    The growth dynamics of the trees that produce fruits are represented in Figure 1, where E = energy; W = water amount in soil; C = crop biomass amount; r = intrinsic growth rate of crops; N = crop carrying capacity; β = evapotranspiration rate; γ = photosynthetic contribution rate; ω = mortality rate of crops; p = rainfall contribution rate; and I = irrigation amount.

    Figure 1.  Dynamics flow. where E = energy; W = water amount in soil; C = crop biomass amount; r = intrinsic growth rate of crop; N = crop carrying capacity; β = evapotranspiration rate; γ = photosynthetic contribution rate; ω = mortality rate of crops; p = rainfall contribution rate and I = irrigation amount. Modified from [4].

    The system of differential equations given by (2.1) presents the theoretical mathematical model that describes the growth of crops with continuous irrigation, as mentioned before:

    {W(t)=(p+I(t))βqW(t)rC(t)(1C(t)N)W(t)C(t)=rC(t)(1C(t)N)W(t)+γqC(t)W(t)(C+1)(W+1)ωC(t) (2.1)

    where I(t) is non-negative, smooth, and bounded function. The qualitative behavior of this model with p=0 and I(t)=0 has been studied in [4], and the results show the local and global stability of the system. The results show that a scenario of water deficit can compromise the growth and development rates of the crop; if this scenario continues, the persistence of the crop is compromised. Therefore, it is important to find an indicator that identifies the point at which the amount of water in the soil does not contribute to the growth of the crops.

    Definition 2.1. The soil water indicator wd is a numerical value based on model parameters that indicate the minimum amount of water that the crop requires for their growth.

    The soil water indicator will be important to determine the water requirements in a crop and the levels at which the water resource should be managed. The indicator (wd) is calculated for system (2.1) following the next generation matrix methodology established in [26].

    Let I(t) be a smooth, monotonous, and bounded function, in the absence of crop the system (2.1) admit the solution (W,0), where W=limteβqt[t0I(s)eβqsds].

    Additionally,

    F=(rC(1CN)W+γqWC(C+1)(W+1)0),V=(ωCβqW+rC(1CN)WI(t)),

    the matrices F and V are the Jacobian matrices of F and V, respectively, and have been defined in [26].

    F(W,0)=(rW+γqWW+1000),V(W,0)=(ω0rWβq),
    V1(W,0)=1ωβq(βq0rWω).

    The next-generation matrix is described as

    FV1(W,0)=(rW+γqWW+1000)(1ω0rWωβq1βq)=(1ω(rW+pWW+1)000).

    The eigenvalues of matrix FV1 are λ1=1ω(rW+γqWW+1), λ2=0.

    max{λ1,λ2} = 1ω(rW+γqWW+1).

    Therefore, the soil water indication (wd) is:

    wd=1ω(rW+γqWW+1). (2.2)

    In what follows, the water soil indicator will be obtained for two irrigation scenarios, and based on this indicator, a criterion will be defined to determine the minimum viable amount for crop growth based on the local stability of the system solutions.

    The calculation of the soil water indicator is hypothetically illustrated in this example. It is not a common scenario in practice, as it is inefficient to consider unlimited water resources. However, despite being unrealistic, this theoretical case serves as an example to analyze the indicator and its relevant parameters. Let us assume that I(t) = I constant. The system (2.1) can be solved using this assumption.

    This example is a hypothetical illustration for the calculation of the soil water indicator. This case is not common in practice because it is not very efficient to consider unlimited water resources. However, despite being unrealistic, this theoretical case provides an example to analyze the indicator and its relevant parameters. Suppose I(t)=I constant, the system (2.1) admits the solution (W,0), where W=Im, with m=βq.

    Additionally,

    F=(rC(1CN)W+γqWC(C+1)(W+1)0),V=(ωCmW+rC(1CN)WI),

    the matrices F and V are determined as

    F(W,0)=(rW+pWW+1000),V(W,0)=(ω0rWm),
    V1(W,0)=1ωm(m0rWω),
    FV1(W,0)=(rW+pWW+1000)(1ω0rWωm1m)=(1ω(rW+pWW+1)000)

    with p=γq, the eigenvalues of matrix FV1 are λ1=1ω(rW+pWW+1), λ2=0.

    max{λ1,λ2} = 1ω(rW+pWW+1).

    Therefore, the soil water indication (wd) is:

    wd=Iω(rm+pI+m). (2.3)

    Lemma 1. If wd<1 the solution (W,0) of system (2.1) is asymptotically stable.

    Proof. The Jacobian matrix of the system (2.1) is

    J=(rW+2rWCNmrC(1CN)rW2rWCN+pW(C+1)2(W+1)ωrC(1CN)+pC(W+1)2(C+1))

    The Jacobian matrix of (2.1) evaluated in the point (W,0) is

    J(W,0)=(rWmrW+pW(W+1)ω0)

    and characteristic polynomial given by:

    P(λ)=λ2+rWλ+m(rW+pW(W+1)ω)

    The Routh-Hurwitz criterion [40] for a quadratic equation states that P(λ) has roots with negative real part if rW>0 and m(rW+pW(W+1)ω)>0, with W=Im and wd=Iω(rm+pI+m). This is rIm>0 and mwd>0 if and only if wd<1. Therefore (W,0) is asymptotically stable.

    In terms of model interpretation, the equilibrium stability implies that the crop tends to decrease until it inevitably disappears when the indicator of water in the soil is less than one. Thus, this result can be interpreted as a criterion that allows defining a threshold for the application of irrigation, and this situation should be avoided.

    The studied scenario provides relevant information on the parameters required to obtain the indicator. In order to establish what are the most relevant parameters in the variation of the value of this indicator, a sensitivity analysis was carried out (see appendix for details). However, it is necessary to consider other irrigation scenarios that reflect the application of this methodology in the field.

    A more realistic scenario for describing the growth dynamics of crops is to carry out irrigation applications at specific points in time. This scenario has been studied from numerical simulations by [41], who proposes the impulsive model described in the equation (2.4)

    {W(t)=βqW(t)rC(t)(1C(t)N)W(t)C(t)=rC(t)(1C(t)N)W(t)+γqC(t)W(t)(C(t)+1)(W(t)+1)ωC(t)}iftnT,W(t+)=W(t)+IC(t+)=C(t)}ift=nT, (2.4)

    where W=W(t),C=C(t) are nonnegative. This model is defined in space

    Ω={(W,C)R2:W0,0C<N},

    and the parameters are all positive in the space ρ={(q,r,N,I,β,γ,ω)R7+}.

    It is important to consider the previously proposed system as a realistic model of the phenomenon to be studied to prove that its solutions are bounded (i.e. there is α depending only on the initial conditions such that for all t>t0, (W(t),C(t))α).

    Lemma 2. The solutions of system (2.4) are bounded.

    Proof. We denote tn=nT, nZ. Clearly 0<t1<t2<...<tn<..., inf{tntn1:n2}0, and limntn=. We define the set σn={(t,W,C)[0,]×R2:t=tn}. At most, the integral curves of the system (2.4) meets the hypersurface σn once. If I is constant, it is possible to consider a continuous function for each choice of n. From the construction of the model, it is clear that the functions that define the system (2.4) are continuously differentiable and locally Lipschitz. Moreover, they satisfy the theorem of existence and uniqueness in each time interval [nT,(n+1)T]. If the solution (W(t),C(t)) is defined for any fixed initial condition in the interval [t0,t] with t>0 constant, from the above conditions, the result follows from Lemma 1 in [42].

    To determine the soil water indicator for (2.4), we follow the methodology outlined in [26,43].

    System (2.4) admits a periodic solution (W,0), where

    W=I1emTem(tnT);t(nT,(n+1)T). (2.5)

    Additionally,

    F=(rC(1CN)W+γqWC(C+1)(W+1)0),V=(ωCβqW+rC(1CN)W),
    (C(nT+)W(nT+))=(1001)(C(nT)W(nT))+(0I).

    According to Theorem 2.1 in [43], the soil water indicator is the solution of ρ(U(T,0,χ))=1, where ρ(U(t,s,χ)) is the evolution operator of the system

    {S(t)=(ω+rW+γqW(W+1)χ)S}iftnT,S(nT+)=S(nT)}ift=nT, (2.6)

    Then ρ(U(T,0,χ))=exp(T0(ω+rW+γqW(W+1)χ)dt). Solving the polynomial about χ,

    exp(T0(ω+rW+γqW(W+1)χ)dt)=1,
    T0(ω+rW+γqW(W+1)χ)dt=0,
    T0ωdt+T0(rW+γqW(W+1)χ)dt=0,
    T0(rW+γqW(W+1))dt=ωTχ,
    T0rWdt+T0γqW(W+1)dt=ωTχ.

    Solving each integral we have:

    Step 1:

    T0rWdt=rIm.

    Step 2:

    T0γqW(W+1)dt=γqm[ln(KemT+1(K+1))], with K=I1emT.

    Therefore, the soil water indicator for the model (2.4) is:

    wd=1ωmT[rIγqln(KemT+1K+1)]. (2.7)

    Equation (2.7) shows the relationship between the parameters considered in the dynamics. Based on this relationship, it is posible to establish criteria for the efficiency of the water-soil-crop dynamics that guarantees crop growth.

    Lemma 3. If wd<1 the periodic solution (W,0) of system (2.4) is asymptotically stable.

    Proof. The linearization of system (2.4) evaluated in (W,0) is given by

    {X(t)=(rW+pW(W+1)ω00m)X(t)}iftnT,X(nT+)=(1001)X(nT)}ift=nT, (2.8)

    the monodromy matrix of the impulsive system (2.8) is

    QΦM(T)=(1001)(ΦFV(T)00Φm(T)) = (exp(T0(rW+pW(W+1)ω)dt)00exp(T0mdt))

    Since wd=exp(T0(rW+pW(W+1)ω)dt)<1, we have ρ(QΦM(T))<1. Therefore the periodic solution (W,0) is asymptotically stable.

    After creating a theoretical mathematical model, validation is necessary to assess its robustness in characterizing the examined phenomenon. This study did not have field data for the model's validation. Thus, this process was performed by reusing published data. In this case, these data were extracted from the works of Ooi et al. [44] and Filippucci et al. [45]. The process consisted of hand-selecting data from some of the original works using the WebPlotDigitizer https://automeris.io/WebPlotDigitizer/, accessed on December 15, 2021 and March 20, 2023, respectively). After their selection, they were extracted and exported to Comma Separated Values (CSV) files for post-processing.

    Regarding the work of Ooi et al. [44], the authors presented the results of the automation of an apple tree orchard irrigation experiment using wireless sensor network technologies in Australia. In this case, for calibration and validation of the proposed models, data was extracted from Figure 2, which describes the measured soil moisture dynamics in the apple orchard, considering irrigation and water-cut cycles. The work of Filippucci et al. [45] presented the results of an irrigation experiment on irrigated tomatoes in Italy. In this case, the trend data for soil moisture was extracted from Figure 3. Please revise their original works for more details about Ooi et al. [44] and Filippucci et al. [45].

    Figure 2.  Model fit from soil moisture data (data extracted from 1 Ooi et al. [44] and 2 Filippucci et al. [45]).
    Figure 3.  Model performance evaluation.

    The soil moisture database obtained from Ooi et al. [44] consisted of 653 pairs of data. In the case of the data from Filippucci et al. [45], 258 pairs of data were record for the soil moisture content (percent). To ensure data independence for calibration and validation, these databases were split into two parts (9 and 642 pairs of data from Ooi et al. [44], and 9 and 249 pairs of data from Filippucci et al. [45].

    The parameterization of the proposed model was carried out using a non-linear least-squares curve fitting method [46], using a specific script written in Matlab©R2019a (Mathworks Inc., Natick, MA, USA). Considering the soil-moisture dynamic, the calibration allowed us to solve the system (2.1) numerically, for both irrigation (water refill) and no-irrigation events (soil-water depletion between each irrigation), adjusting the output W(t) for the soil water content, considering I(t) = 0 (non-irrigation case).

    The calibration allowed for the parameter's adjusting of Table 2 to minimize the difference between the simulated and measured data. The calibrated impulsive models were then validated using the separate database mentioned above. For this purpose, the actual data of soil moisture were compared against those simulated using the classical linear regression method suggested by Mayer and Butler [47]. In this case, the statistical deviance parameters used were the Root Mean Square Error (RMSE), Mean Absolute Error (MAE), and Mean Absolute Percent Error (MAPE). Additionally, the index of agreement (d)} [48,49] and Lin's concordance index (ρ) [50] was added to this analysis (Eq (3.1)).

    MAE=|obsest|nRMSE=(obsest)2nMAPE=100[(|obsest|/|obs|)]nd=1[(est¯obs)(obs¯obs)]2(|est¯obs|+|obs¯obs|)2ρ=2sobsests2obs+s2est+(¯obs¯est)2 (3.1)

    where obs: observed data; est: estimated data; n: the number of pairs [47]; ¯obs: mean of observed; ¯est: mean of estimate data; s2obs: variance of observed data and s2est: variance of estimated data [48,49,50].

    Table 2.  Fitted parameters of the model.
    Parameters From database 1 From database 2
    β 0.0358 0.0004
    q 0.7541 0.0081
    r 0.0002 0.0675
    γ 0.0021 0.0032
    ω 0.0569 0.0037
    where data were extracted from1 Ooi et al. [44] and 2 Filippucci et al. [45].

     | Show Table
    DownLoad: CSV

    Figure 2 shows that the trends of the dynamics of filling and emptying of water in the soil were imitated by the proposed impulsive model (blue line), compared with the real values (points), for the examples obtained from the works from Ooi et al. [44] and from Filippucci et al. [45]. (Figures 2(a) and (b), respectively) In both cases, the observed data trends showed some regularity in the irrigations (irrigation frequency), fulfilling the main assumptions necessary to apply the proposed model.

    On the other hand, the trends of the impulsive model proposed (Eq (2.4)) contrasted against the observed data and had limitations on properly reflecting the soil-moisture dynamics. Figure 2(a) shows} that the modeled values could not reach the highest soil moisture values, presenting a difference close to 1 at the maximum peaks. Regarding the lowest soil moisture values, the model imitated them well. In Figure 2(a), the trends of the soil water content observed showed values in the range of 36 to 42 . In the case of Figure 2b, there was a lag between the observed and modeled soil moisture values, particularly from days 10 to 25. The range of observed values is between 42 and 84 mm, with a difference between the minimum and maximum achievable values of 12 mm. When comparing the observed values against the modeled values using the linear regression method, it was observed that in the first database, there was a trend of the cloud of points to stay above the 1:1 line. In this case, the statistician analysis (Table 3) indicated an MAE, RMSE, MAPE, d, and ρ of 1.19, 1.57, 0.03, 0.70, and 0.46 mm, respectively. Regarding the second database, it shows a uniform cloud of points around the 1:1 line. The statistical analysis indicated a poor performance compared to the first one (MAE = 13.28 mm, RMSE = 16.56 mm, MAPE = 0.24, d = 0.52, and ρ=0.17). These contrasting results show that the performance of the proposed impulsive model will be highly dependent on the assumptions used to calibrate the proposed parameters (Table 2). According to [16,51], these kinds of models are generalist, explaining the studied phenomenon. However, their performance could be poor compared to empirical models specifically calibrated from the data. As indicated herein, compared to Filippucci et al. [45], the database from Ooi et al. [44] presented regular patterns of irrigation events, despite the original data presented in Filippucci et al. [45], which was irregular. This second scenario did not agree with the main assumption of the proposed theoretical model, which was reflected in their poorest performance after the calibration.

    Table 3.  Simulation performance.
    MAE RMSE MAPE d ρ
    Database 1 1.19 mm 1.57 mm 0.03 0.70 0.46
    Database 2 13.28 mm 16.56 mm 0.24 0.52 0.17
    where: Mean Absolute Error (MAE); Root Mean Square Error (RMSE); Mean Absolute Percent Error (MAPE); Index of agreement (d); Lin's concordance index (ρ)

     | Show Table
    DownLoad: CSV

    After carrying out the calibration and validation of the impulsive model, the values established in Table 2 were used to calculate the soil water indicator associated with the equation of the impulsive model (2.4). Initially, we calculate the indicator with the parameters set for database 1, this is:

    wd=1ωmT[rIγqln(KemT+1K+1)]=10.003(0.00360.0015(0.0233))=1.21.

    In this case, the parameters involved in the indicator came from a study where adjustment with the observed data is acceptable. Therefore, the value of the indicator provides a reliable numerical value that implies that the amount of water applied in the period of time presented in Figure 2 (a) favors the water-in-soil-crop interaction.

    The soil water indicator for the parameters associated with database 2 is given by:

    wd=1ωmT[rIγqln(KemT+1K+1)]=116683(0.008770.00002(0.00001))=146.

    In the case of database 2, the indicator is wd>1, which implies that the amount of water supplied in the time described in Figure 2(b) favors the interaction water- cultivation in the soil. However, in this case, the parameters used to calculate this indicator are the result of poor validation, which implies that this numerical value is not reliable for this case. The effectiveness of the indicator will depend on the values of the parameters involved. Thus, the more robust the model, the more accurate the indicator will be.

    To the best of the author's knowledge, various works in the literature highlight the importance of an adequate irrigation technique [8,9,44]. The aim is to establish criteria for the quantification and timing of irrigation. These criteria are based on physiological indicators such as Ψ, gs, and CWSI, among others. In our case, we propose a non-physiological indicator of water in the soil that depends on the parameters established in the water-soil-crop dynamics. The next-generation matrix technique was used to calculate the soil water indicator in this study. This mathematical technique was applied to a non-linear model of differential equations that describe the dynamic between water in soil and crop.

    This dynamic between the water content in soil and crops is presented in two different scenarios. The first scenario describes a relationship between two variables: the amount of water in the soil W(t) and trees that produce fruits biomass amount C(t) at time t; this relationship is modeled in Eq (2.1). In this model the only input of water is given by the application of constant irrigation. The calculation of the soil water indicator is presented in Eq (2.3), which shows that its value depends on the established relationship between the incorporated parameters. In addition, the indicator of water in the soil was considered for the analysis of the local stability of System (2.1). However, this scenario suggests a good example for understanding the calculation of the indicator, and may not reflect the current reality when considering a constant irrigation supply in times of water scarcity and climate crisis.

    The second presented scenario has considered the same relationship between the variables of System (2.1), but it differs in the application of irrigation. In this case, the water supply in the soil is described through the impulsive model given by Eq (2.4), where the irrigation is incorporated at regular intervals. In addition, the soil water indicator Eq (2.7) has been calculated and was used to determine the local stability of the impulsive system (2.4).

    Biological systems present complex relationships that can be difficult to represent due to the number of variables and parameters that must be considered. In that instance, the importance of mathematical modeling in agricultural systems lies in the contribution and understanding of the relationships between climatic variables, water conditions and crops. In this sense, the construction of theoretical models allows for a better understanding of a studied phenomenon. However, the validity of these models depends on the calibration of parameters adjusted to field conditions [8]. In the model that describes the} soil water-crop dynamics presented in this work and used to calculate the soil water indicator, was validated with data provided by [44], showed that the amount of water in the soil presents a trend similar to the data obtained (Figure 2(a)) from the work of [44], the errors between both works are presented in Table 3.

    In addition, as a second exercise, the proposed theoretical model has been compared with data obtained from the work of Filippucci et al. [45] applied to a tomato field. In Figure 2(b), it can be seen that the model behaves similarly to the observed data. However, there is a lag in the irrigation times, which can occur because the observed data tend to present regularity over time; however, this cannot be guaranteed as required by the model assumptions. The errors for this case are presented in Table 3. Furthermore, Lin's concordance index was calculated for both works, and the results obtained suggest a poor performance of the proposed model. This result is because the proposed theoretical mathematical model is an approximate and simplified representation of reality. The proposed model is simple. Its construction is based on restrictive hypotheses and with the minimum possible parameters to describe water-soil-crop dynamics, which do not necessarily reflect actual conditions. On the other hand, the observed data correspond to an empirical model, designed to directly represent those observations [16]. Therefore, a theoretical mathematical model does not necessarily perform well compared to an empirical one because they are based on theoretical assumptions and mathematical equations, not descriptions of observable data [16,51,52,53]. In both cases, the indicator is wd>1, and implies that the amount of water supplied favors the interaction between crop water in the soil. The soil water indicator was calculated in both Ooi et al. [44] and Filippucci et al. [45]. The results obtained suggest that the amount of water supplied is adequate to guarantee water-soil-crop dynamics. However, the results presented depend on the validity of the parameters involved, rather than the method used to calculate them.

    One of the main contributions of this model is the extension of the next-generation matrix method, widely used to establish thresholds for the control of pests in plants. In this case, this method has been used to establish the soil water indicator. The results of this work show that it is possible to describe the crop-water interaction in the soil through mathematical models in a simple way. In addition, a non-physiological indicator of water in the soil wd that provides information on the commitment or permanence of water over time was established. However, it is important to continue investigating the interaction between water in the soil and crop from other scenarios. For example, one should incorporate the irrigation supply at intervals. These possible studies broaden the understanding of the phenomenon and provide a different study approach than the one commonly used in agricultural systems.

    The authors declare they have not used Artificial Intelligence (AI) tools in the creation of this article.

    The present research was partially supported by the Ph.D. scholarship of the Universidad Católica del Maule, and the Gobierno de Chile through the Agencia Nacional de Investigación y Desarrollo (ANID) throughout the "Programa FONDECYT Iniciación en la Investigación, año 2017" (grant No. 11170323) and partially funded by UCM-IN-23202 internal grant.

    The authors declare no conflict of interest.

    Sensitivity analysis

    In the sensitivity analysis developed, we have used the technique proposed in [54,55]. The technique proposes a formula to establish the sensitivity index of all the parameters, this is defined by

    Λwd1y=wd1yxywd1

    where y represent all parameters and wd1=Iω(rβq+γqI+βq). Then

    wd1ωxωwd1=1<0wd1rxrwd1=r(I+βq)r(I+βq)+βγq20wd1γxγwd1=γβq2r(I+βq)+βγq20wd1βxβwd1=r(I+βq)2+γβ2q3(I+βq)[r(I+βq)+γβq2]0wd1qxqwd1=r(I+βq)2+γβq2I(I+βq)[r(I+βq)+γβq2]wd1IxIwd1=r(I+βq)2+γq(βq)2(I+βq)[r(I+βq)+γβq2]0

    In our case, this sensitivity analysis provides a fundamental tool to know the level of influence of each parameter on the soil water indicator. The analysis shows that by increasing the parameters r,γ,I, and considering the other parameters constant, the soil water indicator increases. While the increase of the ω,β parameters decrease the value of the soil water indicator, which can generate compromises in the growth rates of the fruit tree. The analysis for the parameter q indicates that if I>rβγ, then the contribution is negative for 2rβI4rγβI32(γβIrβ2)<q<2rβI+4rγβI32(γβIrβ2).



    [1] Stokes J (2003) Production of coated and free-standing engineering components using the HVOF (High Velocity Oxy-Fuel) process [PhD thesis]. Dublin City University.
    [2] Martin P (2010) Deposition technologies: an overview, Handbook of deposition technologies for films and coatings, 3 Eds., Oxford: Elsevier.
    [3] Halling J (1986) The tribology of surface coatings, particularly ceramics. P I Mech Eng C-J Mec 200: 31–40. doi: 10.1243/PIME_PROC_1986_200_091_02
    [4] Halling J, Nuri K (1985) The elastic contact of rough surfaces and its importance in the reduction of wear. P I Mech Eng C-J Mec 199: 139–144. doi: 10.1243/PIME_PROC_1985_199_104_02
    [5] Martin PM (2009) Handbook of deposition technologies for films and coatings: science, applications and technology, William Andrew.
    [6] Bräuer G (2014) Magnetron Sputtering, In: Hashmi S, Comprehensive Materials Processing.
    [7] Martin P (2011) Introduction to surface engineering and functionally engineered materials, John Wiley & Sons.
    [8] Bhushan B, Gupta BK (1991) Handbook of tribology: materials, coatings, and surface treatments.
    [9] Toyserkani E, Khajepour A, Corbin SF (2004) Laser cladding, CRC Press.
    [10] Mahamood RM, Akinlabi ET (2016) Laser Additive Manufacturing, In: Akinlabi ET, Mahamood RM, Akinlabi SA, Anonymous Advanced Manufacturing Techniques Using Laser Material Processing, IGI Global, 1–23.
    [11] Kashani H, Amadeh A, Ghasemi H (2007) Room and high temperature wear behaviors of nickel and cobalt base weld overlay coatings on hot forging dies. Wear 262: 800–806. doi: 10.1016/j.wear.2006.08.028
    [12] Padaki M, Isloor AM, Nagaraja K, et al. (2011) Conversion of microfiltration membrane into nanofiltration membrane by vapour phase deposition of aluminium for desalination application. Desalination 274: 177–181. doi: 10.1016/j.desal.2011.02.007
    [13] Mattox DM (2010) Handbook of physical vapor deposition (PVD) processing, William Andrew.
    [14] Park S, Ikegami T, Ebihara K, et al. (2006) Structure and properties of transparent conductive doped ZnO films by pulsed laser deposition. Appl Surf Sci 253: 1522–1527. doi: 10.1016/j.apsusc.2006.02.046
    [15] Helmersson U, Lattemann M, Bohlmark J, et al. (2006) Review Ionized physical vapor deposition (IPVD): A review of technology and applications. Thin Solid Films 513: 1–24. doi: 10.1016/j.tsf.2006.03.033
    [16] Somogyvári Z, Langer GA, Erdélyi G, et al. (2012) Sputtering yields for low-energy Ar+- and Ne+-ion bombardment. Vacuum 86: 1979–1982. doi: 10.1016/j.vacuum.2012.03.055
    [17] Roy D, Halder N, Chowdhury T, et al. (2015) Effects of Sputtering Process Parameters for PVD Based MEMS Design. IOSR J VLSI Signal Process 5: 69–77.
    [18] Carlsson J, Martin PM (2010) Chapter 7: Chemical Vapor Deposition, In: Martin PM, Handbook of Deposition Technologies for Films and Coatings: science, applications and technology, 3 Eds., Boston: William Andrew Publishing, 314–363.
    [19] Voevodin A, O'Neill J, Prasad S, et al. (1999) Nanocrystalline WC and WC/a-C composite coatings produced from intersected plasma fluxes at low deposition temperatures. J Vac Sci Technol A 17: 986–992. doi: 10.1116/1.581674
    [20] Mattox DM (2010) Handbook of Physical Vapor Deposition (PVD) processing, William Andrew. 
    [21] Trajkovska-Petkoska A, Nasov I (2014) Surface engineering of polymers: Case study: PVD coatings on polymers. Zaštita materijala 55: 3–10. doi: 10.5937/ZasMat1401003T
    [22] Pretorius R, Marais T, Theron C (1993) Thin film compound phase formation sequence: An effective heat of formation model. Mater Sci Rep 10: 1–83. doi: 10.1016/0920-2307(93)90003-W
    [23] Alami J (2005) Plasma Characterization Thin Film Growth and Analysis in Highly Ionized Magnetron Sputtering [PhD thesis]. Linköping University.
    [24] Toku H, Pessoa RS, Maciel HS, et al. (2010) Influence of Process Parameters on the Growth of Pure-Phase Anatase and Rutile TiO2 Thin Films Deposited by Low Temperature Reactive Magnetron Sputtering. Braz J Phys 40: 340–343.
    [25] Seshan K (2012) Handbook of thin film deposition, William Andrew.
    [26] Bunshah RF (1982) Deposition technologies for films and coatings: developments and applications, Noyes Publications.
    [27] Mattox DM (1998) Atomistic Film Growth and Some Growth-Related Film Properties, Handbook of Physical Vapor Deposition (PVD) Processing.
    [28] Mahan JE (2000) Physical vapor deposition of thin films, Wiley-VCH.
    [29] Elshabini A, Elshabini-Riad AA, Barlow FD (1998) Thin film technology handbook, McGraw-Hill Professional.
    [30] Holmberg K, Mathews A (1994) Coatings tribology: a concept, critical aspects and future directions. Thin Solid Films 253: 173–178. doi: 10.1016/0040-6090(94)90315-8
    [31] Kelly P, Arnell R (2000) Magnetron sputtering: a review of recent developments and applications. Vacuum 56: 159–172. doi: 10.1016/S0042-207X(99)00189-X
    [32] Kelly P, Arnell R, Ahmed W (1993) Some recent applications of materials deposited by unbalanced magnetron sputiering. Surf Eng 9: 287–292. doi: 10.1179/sur.1993.9.4.287
    [33] Ricciardi S (2012) Surface Chemical Functionalization based on Plasma Techniques, Lamberts Academics Publishing.
    [34] Cao G (2004) Nanostructures and nanomaterials: synthesis, properties and applications, 2 Eds., World Scientific.
    [35] Leydecker S (2008) Nano materials: in architecture, interior architecture and design, Walter de Gruyter.
    [36] Phillips RW, Raksha V (2003) Methods for producing enhanced interference pigments. U.S. Patent 6,524,381.
    [37] Herman MA, Sitter H (2012) Molecular beam epitaxy: fundamentals and current status, Springer Science & Business Media.
    [38] Chen Y, Bagnall D, Koh H, et al. (1998) Plasma assisted molecular beam epitaxy of ZnO on c-plane sapphire: growth and characterization. J Appl Phys 84: 3912–3918. doi: 10.1063/1.368595
    [39] Rinaldi F (2002) Basics of molecular beam epitaxy (MBE). Annual Report 2002, Optoelectronics Department, University of Ulm, 1–8.
    [40] Cho AY, Reinhart FK (1975) MBE technique for fabricating semiconductor devices having low series resistance. U.S. Patent 3,915,765.
    [41] Mergel D, Buschendorf D, Eggert S, et al. (2000) Density and refractive index of TiO2 films prepared by reactive evaporation. Thin Solid Films 371: 218–224. doi: 10.1016/S0040-6090(00)01015-4
    [42] Pulker HK, Paesold G, Ritter E (1976) Refractive indices of TiO2 films produced by reactive evaporation of various titanium–oxygen phases. Appl Optics 15: 2986–2991. doi: 10.1364/AO.15.002986
    [43] Terashima T, Iijima K, Yamamoto K, et al. (1988) Single-crystal YBa2Cu3O7−x thin films by activated reactive evaporation. Jpn J Appl Phys 27: L91. doi: 10.1143/JJAP.27.L91
    [44] Bunshah R, Raghuram A (1972) Activated reactive evaporation process for high rate deposition of compounds. J Vac Sci Technol 9: 1385–1388. doi: 10.1116/1.1317045
    [45] Bunshah R (1981) The activated reactive evaporation process: Developments and applications. Thin Solid Films 80: 255–261. doi: 10.1016/0040-6090(81)90231-5
    [46] Schultz PG, Xiang X, Goldwasser I, et al. (2006) Combinatorial synthesis and screening of non-biological polymers. U.S. Patent 7,034,091.
    [47] Seyfert U, Heisig U, Teschner G, et al. (2015) 40 Years of Industrial Magnetron Sputtering in Europe. SVC Bulletin Fall 2015: 22–26. 
    [48] Seshan K (2001) Handbook of thin-film deposition processes and Techniques, Principles, Methods, Equipment and Applications, Noyes Publications/William Andrew Publishing.
    [49] Teixeira V, Cui H, Meng L, et al. (2002) Amorphous ITO thin films prepared by DC sputtering for electrochromic applications. Thin Solid Films 420: 70–75.
    [50] Utsumi K, Iigusa H, Tokumaru R, et al. (2003) Study on In2O3–SnO2 transparent and conductive films prepared by d.c. sputtering using high density ceramic targets. Thin Solid Films 445: 229–234.
    [51] Gómez A, Galeano A, Saldarriaga W, et al. (2015) Deposition of YBaCo4O7 thin films on (001)-SrTiO3 substrates by dc sputtering. Vacuum 119: 7–14. doi: 10.1016/j.vacuum.2015.04.020
    [52] Cash Jr JH, Cunningham JA (1972) Rf sputtering method. U.S. Patent 3,677,924.
    [53] You T, Niwa O, Tomita M, et al. (2002) Characterization and electrochemical properties of highly dispersed copper oxide/hydroxide nanoparticles in graphite-like carbon films prepared by RF sputtering method. Electrochem Commun 4: 468–471. doi: 10.1016/S1388-2481(02)00340-5
    [54] Torng C, Sivertsen JM, Judy JH, et al. (1990) Structure and bonding studies of the C: N thin films produced by rf sputtering method. J Mater Res 5: 2490–2496. doi: 10.1557/JMR.1990.2490
    [55] Kelly PJ, Arnell RD (2000) Magnetron sputtering: a review of recent developments and applications. Vacuum 56: 159–172. doi: 10.1016/S0042-207X(99)00189-X
    [56] Constantin DG, Apreutesei M, Arvinte R, et al. (2011) Magnetron sputtering technique used for coatings deposition; technologies and applications. 7th International Conference on Materials Science and Engineering.
    [57] Magnetron Sputtering Technology, 2019. Available from: http://www.directvacuum.com/pdf/ what_is_sputtering.pdf.
    [58] David JC (2014) Making Magnetron Sputtering Work: Modelling Reactive Sputtering Dynamics, Part 1. Available from: https://www.svc.org/DigitalLibrary/documents/2014_Fall_DC.pdf.
    [59] Arnell RD, Kelly PJ (1999) Recent advances in magnetron sputtering. Surf Coat Tech 112: 170–176. doi: 10.1016/S0257-8972(98)00749-X
    [60] Musil J (1998) Recent advances in magnetron sputtering technology. Surf Coat Tech 100: 280–286.
    [61] Suzuki Y, Teranishi K (2009) Reactive sputtering method. U.S. Patent 7,575,661.
    [62] Ishihara M, Li S, Yumoto H, et al. (1998) Control of preferential orientation of AlN films prepared by the reactive sputtering method. Thin Solid Films 316: 152–157. doi: 10.1016/S0040-6090(98)00406-4
    [63] Chistyakov R (2006) High-power pulsed magnetron sputtering. U.S. Patent 7,147,759.
    [64] Kouznetsov V, Macák K, Schneider JM, et al. (1999) A novel pulsed magnetron sputter technique utilizing very high target power densities. Surf Coat Tech 122: 290–293. doi: 10.1016/S0257-8972(99)00292-3
    [65] Sarakinos K, Alami J, Konstantinidis S (2010) High power pulsed magnetron sputtering: A review on scientific and engineering state of the art. Surf Coat Tech 204: 1661–1684. doi: 10.1016/j.surfcoat.2009.11.013
    [66] Alami J, Bolz S, Sarakinos K (2009) High power pulsed magnetron sputtering: Fundamentals and applications. J Alloy Compd 483: 530–534. doi: 10.1016/j.jallcom.2008.08.104
    [67] Wasa K, Hayakawa S (1992) Handbook of sputter deposition technology, Noyes Publications
    [68] Wasa K, Kitabatake M, Adachi H (2004) Thin film materials technology: sputtering of control compound materials, Springer Science & Business Media.
    [69] Wasa K, Kitabatake M, Adachi H (2004) Thin Film Processes, In: Wasa K, Kitabatake M, Adachi H, Thin Film Materials Technology, Norwich, NY: William Andrew Publishing, 17–69.
    [70] Liu X, Poon RWY, Kwok SCH, et al. (2004) Plasma surface modification of titanium for hard tissue replacements. Surf Coat Tech 186: 227–233. doi: 10.1016/j.surfcoat.2004.02.045
    [71] Liu X, Chu PK, Ding C (2005) Surface modification of titanium, titanium alloys, and related materials for biomedical applications. Mat Sci Eng R 47: 49–121.
    [72] Albetran HMM (2016) Synthesis and characterisation of nanostructured TiO2 for photocatalytic applications [PhD thesis]. Curtin University.
    [73] Miller T (2010) An Investigation into the Growth and Characterisation of Thin Film Radioluminescent Phosphors for Neutron Diffraction Analysis [PhD thesis]. Nottingham Trent University.
    [74] Pessoa R, Fraga M, Santos L, et al. (2015) Nanostructured thin films based on TiO2 and/or SiC for use in photoelectrochemical cells: A review of the material characteristics, synthesis and recent applications. Mat Sci Semicon Proc 29: 56–68. doi: 10.1016/j.mssp.2014.05.053
    [75] Kommu S, Wilson GM, Khomami B (2000) A Theoretical/Experimental Study of Silicon Epitaxy in Horizontal Single‐Wafer Chemical Vapor Deposition Reactors. J Electrochem Soc 147: 1538–1550. doi: 10.1149/1.1393391
    [76] Pedersen H, Elliott SD (2014) Studying chemical vapor deposition processes with theoretical chemistry. Theor Chem Acc 133: 1476. doi: 10.1007/s00214-014-1476-7
    [77] Fuller CB, Mahoney MW, Bingel WH (2006) Friction stir weld tool and method. U.S. Patent 6,994,242.
    [78] Klamklang S (2007) Restaurant wastewater treatment by electrochemical oxidation in continuous process [PhD thesis].
    [79] Chemical Vapour Deposition, 2016. Available from: http://users.wfu.edu/ucerkb/Nan242/L09-CVD_a.pdf.
    [80] Wang DN, White JM, Law KS, et al. (1991) Thermal CVD/PECVD reactor and use for thermal chemical vapor deposition of silicon dioxide and in-situ multi-step planarized process. U.S. Patent 5,000,113.
    [81] Hirose Y, Terasawa Y (1986) Synthesis of diamond thin films by thermal CVD using organic compounds. Jpn J Appl Phys 25: L519. doi: 10.1143/JJAP.25.L519
    [82] Petzoldt F, Piglmayer K, Kräuter W, et al. (1984) Lateral growth rates in laser CVD of microstructures. Appl Phys A-Mater 35: 155–159. doi: 10.1007/BF00616969
    [83] Sousa P, Silvestre A, Popovici N, et al. (2005) Morphological and structural characterization of CrO2/Cr2O3 films grown by Laser-CVD. Appl Surf Sci 247: 423–428. doi: 10.1016/j.apsusc.2005.01.061
    [84] Matsui S, Kaito T, Fujita J, et al. (2000) Three-dimensional nanostructure fabrication by focused-ion-beam chemical vapor deposition. J Vac Sci Technol B 18: 3181–3184. doi: 10.1116/1.1319689
    [85] Inoue K, Michimori M, Okuyama M, et al. (1987) Low temperature growth of SiO2 thin film by double-excitation photo-CVD. Jpn J Appl Phys 26: 805.
    [86] Tanimoto S, Matsui M, Kamisako K, et al. (1992) Investigation on leakage current reduction of photo‐CVD tantalum oxide films accomplished by active oxygen annealing. J Electrochem Soc 139: 320–328. doi: 10.1149/1.2069193
    [87] Price J, Wu S (1987) Plasma enhanced CVD. U.S. Patent 4,692,343.
    [88] Li Y, Mann D, Rolandi M, et al. (2004) Preferential growth of semiconducting single-walled carbon nanotubes by a plasma enhanced CVD method. Nano Lett 4: 317–321. doi: 10.1021/nl035097c
    [89] Hozumi A, Takai O (1997) Preparation of ultra water-repellent films by microwave plasma-enhanced CVD. Thin Solid Films 303: 222–225. doi: 10.1016/S0040-6090(97)00076-X
    [90] Hitchman ML, Jensen KF (1993) Chemical vapor deposition: principles and applications, Elsevier.
    [91] Graniel O, Weber M, Balme S, et al. (2018) Atomic layer deposition for biosensing applications. Biosens Bioelectron 122: 147–159. doi: 10.1016/j.bios.2018.09.038
    [92] Suntola T, Antson J (1977) Method for producing compound thin films. U.S. Patent 4,058,430.
    [93] Guo HC, Ye E, Li Z, et al. (2017) Recent progress of atomic layer deposition on polymeric materials. Mat Sci Eng C-Mater 70: 1182–1191. doi: 10.1016/j.msec.2016.01.093
    [94] George SM (2009) Atomic layer deposition: an overview. Chem Rev 110: 111–131.
    [95] Miikkulainen V, Leskelä M, Ritala M, et al. (2013) Crystallinity of inorganic films grown by atomic layer deposition: Overview and general trends. J Appl Phys 113: 2.
    [96] Bohr MT, Chau RS, Ghani T, et al. (2007) The high-k solution. IEEE Spectrum 44: 29–35. doi: 10.1109/MSPEC.2007.4337663
    [97] Yan B, Li X, Bai Z, et al. (2017) A review of atomic layer deposition providing high performance lithium sulfur batteries. J Power Sources 338: 34–48. doi: 10.1016/j.jpowsour.2016.10.097
    [98] Ozer N, Lampert CM (1998) Electrochromic characterization of sol–gel deposited coatings. Sol Energ Mat Sol C 54: 147–156. doi: 10.1016/S0927-0248(98)00065-8
    [99] Ghoranneviss ZG (2016) Effects of various deposition times and RF powers on CdTe thin film growth using magnetron sputtering. J Theor Appl Phys 10: 225–231. doi: 10.1007/s40094-016-0219-7
    [100] Alfonso E, Cubillos G, Olaya J (2012) Thin film growth through sputtering technique and its applications, INTECH Open Access Publisher.
    [101] Lewis B, Anderson JC (1978) Nucleation and growth of thin films, New York: Academic Press.
    [102] Venables JA, Spiller GDT (1983) Nucleation and growth of thin films, In: Surface Mobilities on Solid Materials, Springer, 341–404.
    [103] Lewis B, Campbell D (1967) Nucleation and initial-growth behavior of thin-film deposits. J Vac Sci Technol 4: 209–218. doi: 10.1116/1.1492548
    [104] Lane G, Anderson J (1975) The nucleation and initial growth of gold films deposited onto sodium chloride by ion-beam sputtering. Thin Solid Films 26: 5–23. doi: 10.1016/0040-6090(75)90164-9
    [105] Movchan B, Demchishin A (1969) Structure and properties of thick condensates of nickel, titanium, tungsten, aluminum oxides, and zirconium dioxide in vacuum. Fiz Metal Metalloved 28: 653–660.
    [106] Thornton JA (1977) High rate thick film growth. Annu Rev Mater Sci 7: 239–260. doi: 10.1146/annurev.ms.07.080177.001323
    [107] Messier R, Giri A, Roy R (1984) Revised structure zone model for thin film physical structure. J Vac Sci Technol A 2: 500–503. doi: 10.1116/1.572604
    [108] Singh V (2004) Synthesis, structure, and tribological behavior of nanocomposite DLC based thin films [PhD thesis]. Louisiana State University.
    [109] Dresden-rossendorf I, Doctor G (2010) Growth, structure and magnetic properties of magnetron sputtered FePt thin films [PhD thesis]. Technische Universität Dresden.
    [110] Wang B, Fu X, Song S, et al. (2018) Simulation and Optimization of Film Thickness Uniformity in Physical Vapor Deposition. Coatings 8: 325. doi: 10.3390/coatings8090325
    [111] Divi S, Chatterjee A (2014) Study of Silicon Thin Film Growth at High Deposition Rates Using Parallel Replica Molecular Dynamics Simulations. Energy Procedia 54: 270–280. doi: 10.1016/j.egypro.2014.07.270
    [112] Allen MP (2004) Introduction to molecular dynamics simulation, In: Attig N, Binder K, Grubmuller H, Computational soft matter: from synthetic polymers to proteins, 23: 1–28.
    [113] Sutmann G (2002) Classical molecular dynamics and parallel computing, Germany: FZJ-ZAM-IB.
    [114] Rasoulian S, Ricardez-Sandoval LA (2015) Worst-case and distributional robustness analysis of a thin film deposition process. 9th International Symposium on Advanced Control of Chemical Processes, 7–10.
    [115] Li J, Croiset E, Ricardez-Sandoval L (2015) Carbon nanotube growth: First-principles-based kinetic Monte Carlo model. J Catal 326: 15–25. doi: 10.1016/j.jcat.2015.03.010
    [116] Family F (1986) Scaling of rough surfaces: effects of surface diffusion. J Phys A 19: 441. doi: 10.1088/0305-4470/19/8/006
    [117] Tang L, Nattermann T (1991) Kinetic roughening in molecular-beam epitaxy. Phys Rev Lett 66: 2899. doi: 10.1103/PhysRevLett.66.2899
    [118] Bhute VJ, Chatterjee A (2013) Building a kinetic Monte Carlo model with a chosen accuracy. J Chem Phys 138: 244112. doi: 10.1063/1.4812319
    [119] Bhute VJ, Chatterjee A (2013) Accuracy of a Markov state model generated by searching for basin escape pathways. J Chem Phys 138: 084103. doi: 10.1063/1.4792439
    [120] Lou Y, Christofides PD (2003) Estimation and control of surface roughness in thin film growth using kinetic Monte-Carlo models. Chem Eng Sci 58: 3115–3129. doi: 10.1016/S0009-2509(03)00166-0
    [121] Theodoropoulou A, Adomaitis RA, Zafiriou E (1998) Model reduction for optimization of rapid thermal chemical vapor deposition systems. IEEE T Semiconduct M 11: 85–98. doi: 10.1109/66.661288
    [122] Middlebrooks SA, Rawlings JB (2007) Model Predictive Control of Si1−xGex Thin Film Chemical–Vapor Deposition. IEEE T Semiconduct M 20: 114–125. doi: 10.1109/TSM.2007.895203
    [123] Rasoulian S, Ricardez-Sandoval LA (2015) Robust multivariable estimation and control in an epitaxial thin film growth process under uncertainty. J Process Contr 34: 70–81. doi: 10.1016/j.jprocont.2015.07.002
    [124] Ricardez‐Sandoval LA (2011) Current challenges in the design and control of multiscale systems. Can J Chem Eng 89: 1324–1341. doi: 10.1002/cjce.20607
    [125] Vlachos DG (2005) A review of multiscale analysis: examples from systems biology, materials engineering, and other fluid–surface interacting systems. Adv Chem Eng 30: 1–61. doi: 10.1016/S0065-2377(05)30001-9
    [126] Chaffart D, Ricardez-Sandoval LA (2018) Optimization and control of a thin film growth process: A hybrid first principles/artificial neural network based multiscale modelling approach. Comput Chem Eng 119: 465–479. doi: 10.1016/j.compchemeng.2018.08.029
    [127] Jensen KF, Rodgers ST, Venkataramani R (1998) Multiscale modeling of thin film growth. Curr Opin Solid St M 3: 562–569. doi: 10.1016/S1359-0286(98)80026-0
    [128] Baumann F, Chopp D, De La Rubia TD, et al. (2001) Multiscale modeling of thin-film deposition: applications to Si device processing. MRS Bull 26: 182–189. doi: 10.1557/mrs2001.40
    [129] Zhang P, Zheng X, Wu S, et al. (2004) Kinetic Monte Carlo simulation of Cu thin film growth. Vacuum 72: 405–410. doi: 10.1016/j.vacuum.2003.08.013
    [130] Evans RD, Ricardez-Sandoval LA (2014) Multi-scenario modelling of uncertainty in stochastic chemical systems. J Comput Phys 273: 374–392. doi: 10.1016/j.jcp.2014.05.028
    [131] Chaffart D, Ricardez-Sandoval LA (2017) Robust dynamic optimization in heterogeneous multiscale catalytic flow reactors using polynomial chaos expansion. J Process Contr 60: 128–140. doi: 10.1016/j.jprocont.2017.07.002
    [132] Rasoulian S, Ricardez-Sandoval LA (2016) Stochastic nonlinear model predictive control applied to a thin film deposition process under uncertainty. Chem Eng Sci 140: 90–103. doi: 10.1016/j.ces.2015.10.004
    [133] Rasoulian S, Ricardez-Sandoval LA (2015) A robust nonlinear model predictive controller for a multiscale thin film deposition process. Chem Eng Sci 136: 38–49. doi: 10.1016/j.ces.2015.02.002
    [134] Allgöwer F, Findeisen R, Nagy ZK (2004) Nonlinear model predictive control: From theory to application. J Chin Inst Chem Engrs 35: 299–315.
    [135] Hu G, Orkoulas G, Christofides PD (2009) Modeling and control of film porosity in thin film deposition. Chem Eng Sci 64: 3668–3682. doi: 10.1016/j.ces.2009.05.008
    [136] Li C, Song S, Gibson D, et al. (2017) Modeling and validation of uniform large-area optical coating deposition on a rotating drum using microwave plasma reactive sputtering. Appl Optics 56: C65–C70. doi: 10.1364/AO.56.000C65
    [137] Heirung TAN, Paulson JA, O'Leary J, et al. (2018) Stochastic model predictive control-how does it work? Comput Chem Eng 114: 158–170. doi: 10.1016/j.compchemeng.2017.10.026
    [138] Koronaki E, Gkinis P, Beex L, et al. (2018) Classification of states and model order reduction of large scale Chemical Vapor Deposition processes with solution multiplicity. Comput Chem Eng 121: 148–157.
    [139] Divi S, Chatterjee A (2014) Study of Silicon Thin Film Growth at High Deposition Rates Using Parallel Replica Molecular Dynamics Simulations. Energy Procedia 54: 270–280. doi: 10.1016/j.egypro.2014.07.270
  • Reader Comments
  • © 2019 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(37105) PDF downloads(17320) Cited by(190)

Figures and Tables

Figures(4)

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog