Research article Special Issues

Optimizing MDA and antimalarial treatment in the presence of drug resistance for effective malaria control

  • Received: 17 January 2025 Revised: 27 May 2025 Accepted: 05 June 2025 Published: 18 June 2025
  • Antimalarial drugs are critical for controlling malaria, but the emergence of drug resistance poses a significant challenge to global eradication efforts. This study explores strategies to minimize resistance prevalence and improve malaria control, particularly through the use of mass drug administration (MDA) in combination with antimalarial drugs. We develop a compartmental mathematical model that incorporates asymptomatic, paucisymptomatic, and clinical states of infection and evaluates the impact of resistance mutations on transmission dynamics. The model includes both treated and untreated states among infected and recovered individuals, with a focus on optimizing control strategies through MDA and antimalarial treatment. A global sensitivity analysis identifies the critical factors that influence malaria dynamics, including MDA coverage, treatment access for different infection states, the probability of mutation from treated sensitive human infections, to treated resistant human infections and the initial prevalence of resistance. The model is extended to include optimal control strategies that provide time-dependent control interventions for treatment and MDA. Intuitively, when the mutation rate is relatively low, the optimal strategy combines the use of antimalarial drugs and MDA, with a gradual decrease in antimalarial drug use over time, ensuring sustainable malaria control. In contrast, at higher mutation rates, the strategy prioritizes broader deployment of MDA while significantly reducing reliance on antimalarial to minimize the risk of resistance developing. Numerical simulations of the optimal control problem reinforce the importance of strategic intervention in mitigating drug resistance. This study contributes to understanding the role of MDA and treatment strategies in the control of malaria, with implications for optimizing malaria control programs in endemic regions.

    Citation: Manuela M. Nimpa, Hyacinthe N. Teytsa, Joseph Mbang, Charles S. Wondji, Ramsès Djidjou-Demasse. Optimizing MDA and antimalarial treatment in the presence of drug resistance for effective malaria control[J]. Mathematical Biosciences and Engineering, 2025, 22(8): 1898-1930. doi: 10.3934/mbe.2025069

    Related Papers:

    [1] Uichi Akiba, Jun-ichi Anzai . Cyclodextrin-containing layer-by-layer films and microcapsules: Synthesis and applications. AIMS Materials Science, 2017, 4(4): 832-846. doi: 10.3934/matersci.2017.4.832
    [2] Said Karim Shah, Jahangeer Khan, Irfan Ullah, Yaqoob Khan . Optimization of active-layer thickness, top electrode and annealing temperature for polymeric solar cells. AIMS Materials Science, 2017, 4(3): 789-799. doi: 10.3934/matersci.2017.3.789
    [3] Rucheng Zhu, Yota Mabuchi, Riteshkumar Vishwakarma, Balaram Paudel Jaisi, Haibin Li, Masami Naito, Masayoshi Umeno, Tetsuo Soga . Enhancing in-plane uniformity of graphene nanowalls using a rotating platform for solid-state lithium-ion battery. AIMS Materials Science, 2024, 11(4): 760-773. doi: 10.3934/matersci.2024037
    [4] Yuguo Tao, Vijaykumar Upadhyaya, Keenan Jones, Ajeet Rohatgi . Tunnel oxide passivated rear contact for large area n-type front junction silicon solar cells providing excellent carrier selectivity. AIMS Materials Science, 2016, 3(1): 180-189. doi: 10.3934/matersci.2016.1.180
    [5] Stefan Wagner, Astrid Pundt . Hydrogen as a probe for defects in materials: Isotherms and related microstructures of palladium-hydrogen thin films. AIMS Materials Science, 2020, 7(4): 399-419. doi: 10.3934/matersci.2020.4.399
    [6] Silvia Colodrero . Conjugated polymers as functional hole selective layers in efficient metal halide perovskite solar cells. AIMS Materials Science, 2017, 4(4): 956-969. doi: 10.3934/matersci.2017.4.956
    [7] Hmoud. Al-Dmour . Capacitance response of solar cells based on amorphous Titanium dioxide (A-TiO2) semiconducting heterojunctions. AIMS Materials Science, 2021, 8(2): 261-270. doi: 10.3934/matersci.2021017
    [8] Ouassim Hamdi, Frej Mighri, Denis Rodrigue . Piezoelectric cellular polymer films: Fabrication, properties and applications. AIMS Materials Science, 2018, 5(5): 845-869. doi: 10.3934/matersci.2018.5.845
    [9] Thomas Debesay, Sam-Shajing Sun, Messaoud Bahoura . Multi-functional organic field effect transistor based on a dual doped P3HT. AIMS Materials Science, 2021, 8(5): 823-835. doi: 10.3934/matersci.2021050
    [10] Nur Jassriatul Aida binti Jamaludin, Shanmugan Subramani . Thermal performance of LED fixed on CVD processed ZnO thin film on Al substrates at various O2 gas flow rates. AIMS Materials Science, 2018, 5(2): 246-256. doi: 10.3934/matersci.2018.2.246
  • Antimalarial drugs are critical for controlling malaria, but the emergence of drug resistance poses a significant challenge to global eradication efforts. This study explores strategies to minimize resistance prevalence and improve malaria control, particularly through the use of mass drug administration (MDA) in combination with antimalarial drugs. We develop a compartmental mathematical model that incorporates asymptomatic, paucisymptomatic, and clinical states of infection and evaluates the impact of resistance mutations on transmission dynamics. The model includes both treated and untreated states among infected and recovered individuals, with a focus on optimizing control strategies through MDA and antimalarial treatment. A global sensitivity analysis identifies the critical factors that influence malaria dynamics, including MDA coverage, treatment access for different infection states, the probability of mutation from treated sensitive human infections, to treated resistant human infections and the initial prevalence of resistance. The model is extended to include optimal control strategies that provide time-dependent control interventions for treatment and MDA. Intuitively, when the mutation rate is relatively low, the optimal strategy combines the use of antimalarial drugs and MDA, with a gradual decrease in antimalarial drug use over time, ensuring sustainable malaria control. In contrast, at higher mutation rates, the strategy prioritizes broader deployment of MDA while significantly reducing reliance on antimalarial to minimize the risk of resistance developing. Numerical simulations of the optimal control problem reinforce the importance of strategic intervention in mitigating drug resistance. This study contributes to understanding the role of MDA and treatment strategies in the control of malaria, with implications for optimizing malaria control programs in endemic regions.



    1. Introduction

    For more than five decades, quasi-two-dimensional structures in the form of quantum semiconductor films (quantum wells) have been intensively studied both theoretically and experimentally [1,2,3,4,5,6,7,8,9,10]. They very quickly found a variety of instrumental applications and occupy their solid place in modern solid-state electronics [11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28]. Such a wide application of these quasi-two-dimensional quantum structures is mainly due to two factors, the presence of which is fundamentally impossible in a three-dimensional sample. This is the phenomenon of the size quantization of charge carriers in one of the directions of motion [1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28] and corresponding to this phenomenon the specificity of the quasi-two-dimensional Fermi distribution of electrons and holes in the film [1,2,29,30,31,32,33,34,35,36,37,38]. In reality, being a three-dimensional system, a quantum semiconductor film exhibits a number of properties characteristic of 2D systems in the field of multiparticle phenomena (in fact, the two-dimensional nature of the Coulomb coupling in the sample, the two-dimensional dispersion law and the two-dimensional density of states for carriers, etc.).

    The object of our consideration in this paper will be precisely questions of the statistics of charge carriers in a quantum semiconductor film. The general questions of the properties of the Fermi-Dirac distribution in a semiconductor film, when dimensional quantization of carriers in one direction takes place are considered in Refs. [1,2,29,30,31,36,37]. In these papers, the dependence of the chemical potential of the electronic subsystem on the thickness of the film is also considered. The papers [32,33,34,35] are devoted to the determination of the chemical potential of the electronic subsystem in quasi-two-dimensional systems, and in particular, in quantum semiconductor films.

    Chemical potential and internal energy of a noninteracting Fermi gas at low temperature are evaluated using the Sommerfeld method in the fractional-and quasi-two-dimensional spaces in Ref. [34]. The effects of low dimensionality on the thermodynamics of a Fermi gas trapped by isotropic power law potentials are analyzed in Ref. [38].

    Simultaneously it is well known that along with size quantization, one of the powerful external modulating factors is an external homogeneous electrostatic field. In the case of a semiconductor film, the situation where the field is directed along the quantization direction is of greatest physical interest. In this case, the states of the discrete spectrum of charge carriers undergo an additional, field modulation. This leads to a change in the kinetic, optical, and a number of other properties of the semiconductor film [39,40,41,42,43,44]. Naturally, the external electrostatic field exerts its significant influence also on the statistical properties of a quantum semiconductor film.

    In the present paper, the influence of a strong external lateral homogeneous electrostatic field on the statistical properties of the electronic subsystem in a quantum narrow-band semiconductor film is theoretically considered, when a regime of "strong quantization" is realized for charge carriers in a film.


    2. Single-electron states in the film in the presence of strong electrostatic field

    As already noted, in this paper we investigate the case in which a regime of "strong quantization" is realized in a semiconductor film for charge carriers. For this, as is known, at least the following two important physical conditions must be satisfied [1,7,43]:

    (a) The film thickness (L) must be much less than the Bohr radius of the bulk exciton of the film material (aex)

    L2a2ex (1)

    (b) The distance between the levels of the size quantization (Δεn,n+1) must significantly exceed the energy of the thermal motion of charge carriers (kBT):

    Δεn,n+1kBT (2)

    Here kB is Boltzmann constant and T is the absolute temperature of the sample.

    The conditions (1) and (2) dictate the choice of the material under investigation, respectively.

    In a number of narrow-gap semiconductors, owing to the large static permittivity and the small effective mass of charge carriers, the Bohr radius of the bulk exciton has a value of the order of several tens of nanometers: PbS (Eg = 0.32 eV) me/mem0m0=0.04, aex = 20 nm, PbSe (Eg = 0.27 eV) me/mem0m0=0.08, aex = 47 nm, PbTe (Eg = 0.41 eV) me/mem0m0=0.04, aex = 150 nm, HgSe (Eg = 0.24 eV) me/mem0m0=0.016, aex = 52 nm, HgS (Eg = 0.5 eV) me/mem0m0=0.036, aex = 50 nm, InAs (Eg = 0.354 eV) me/mem0m0=0.023, aex = 35 nm, InSb (Eg = 0.18 eV) me/mem0m0=0.014, aex = 60 nm. Here me is the effective mass of electron, m0 is the mass of a free electron, Eg is the width of the band gap of a bulk sample (data are taken from Refs. [45,46,47,48]).

    At a film thickness L~10 nm, the condition (1) for such semiconductors begins with sufficient accuracy even at relatively high temperatures, when the temperature of the system T varies, for example, from 100 to 300 K. The energy of the thermal motion of the electrons varies in this case, respectively, in the range from 10 to 25 meV.

    We now consider single-particle states in a quantizing film in the presence of a strong transverse electrostatic field. Тhe quantizing film is approximated by the model of an infinitely deep rectangular potential well along the direction of quantization (z) [1,7,42]. Such a choice of the form of the model potential of a film is justified by the following physical factors. Firstly, a strong external field creates a deep triangular trap for the particle in the lower part of the well. Secondly, we are interested in the lower levels of quantization. But those levels in this case will be located precisely in the triangular part of the quantum well.

    Then, we can write for the potential energy of the electron in the presence of an external electric field F=F(0;0;F) in the film:

    U(z,F)={;z0,zL|e|Fz;0zL (3)

    Here e is the charge of an electron, L is the thickness of a film.

    In the absence of an external field the energy spectrum of one-electron states in the approximation of an isotropic effective mass, is given by the following expression [1,42]:

    εn(px,py)=p2x+p2y2me+π22n22meL2ε(px,py)+ε1n2;n=1,2,... (4)

    Here px,py are the components of the quasi-momentum of the particle in the plane of the film.

    We regard the external field as strong. Physically, this implies that the energy received by the particle from the field within the film (|e|FL) is much larger than the particle's size quantization energy (ε1):

    |e|FL>>ε1=π22/π222meL22meL2 (5)

    It is shown in Ref. [42], that the strong-field condition with sufficient accuracy can be written as follows

    |e|FL10ε1 (6)

    In the case of an InAs film L = 10 nm, we have for ε1:ε10.1635 eV. Accordingly, the condition (6) for the InAs film will be satisfied if the potential difference between the planes of the film is as follows: Δφ1.5 V.

    When the condition (5) is satisfied, the energy spectrum of a particle in a film in the presence of potential Eq 3 can be determined by a variation method [42]. Taking into account that it is the lower sub-bands that are important for the manifestation of the quantum size effect, we give the expressions obtained for the first two electronic levels in the film in the presence of a strong electrostatic field:

    ε1,e(F)2.345(2e2F22me)1/133=1.093ε1,e(|e|FLε1,e)2/2331.093ε1,eη2/233 (7)
    ε2,e(F)4.109(2e2F22me)1/133=1.931ε1,e(|e|FLε1,e)2/2331.931ε1,eη2/233 (8)

    Note that the conditions (5) and (6) are, in fact, also the quasi-classical conditions for the potential Eq 3. In view of this circumstance, the energy spectrum of the particle motion in the field (Eq 3) under the conditions (5) and (6) can also be determined from the Bohr-Sommerfeld quantization conditions [49]:

    π(l+34)=2mz00εl(F)|e|Fzdz;(z0=εl(F)/εl(F)|e|F|e|F); (9)

    Hence we obtain:

    εl(F)=[32π|e|F2m(l+34)]2/233;l=0,1,2,... (10)

    It is easy to see that both approaches in this case give very similar results. We obtain from the Eq 9:

    εl=0(F)=2.327(2e2F22m)1/133;εl=1(F)=4.084(2e2F22m)1/133 (11)

    We do not write out the shape of the envelope wave functions, since in the future they will not be needed.

    Let us now consider the statistical characteristics of the electron-hole subsystem in a film in the presence of an external field.

    All calculations and estimates will be made using the example of an InAs (meme=0.023m0,mlh=mlh=0.026m0,mhh=mhh=0.41m0,ε(0)g=0.354eV) quantum film at a film thickness L = 1 nm. In such a sample, the conditions (1) and (2) are satisfied with great accuracy.

    The melting temperature of the InAs crystal is equal to T0=1215 K, so that over the entire temperature range of this crystal (T[0,T0]) this semiconductor is non-degenerate: even when T~1000 K, kBT=0.1032eV<Eg=0.354eV.

    We also note that the properties of InAs quantum films are being intensively studied at the present time (for example, [49,50]) and the study proposed here is of real interest from this point of view.


    3. The statistical characteristics of charge carriers in a semiconductor film in the presence of a strong electrostatic field

    In the presence of an external field with a strength F for the energy spectrum of an electron in general form, we can now write:

    εn(p,F)=p22me+εn(F)=ε(p)+εn(F)ε;p=p(px,py,0) (12)

    Here the energy εn(F) is given by the Eqs 7–10.

    The general form of Fermi-Dirac distribution function for charge carriers in this case is as follows:

    f(εn(p,F))=[1+expεn(p,F)μkBT]1 (13)

    Here μ is the chemical potential of charge carriers' sub-system.

    The effect of size quantization, as already noted, manifests itself most effectively in the conditions when the fewest possible sub-bands are filled. Therefore, we confine ourselves here to considering the filling of only the first two size quantization sub-bands in the presence of an external field.

    In the case of the dispersion law from Eq 12, the normalization condition for the distribution function in Eq 13 for electrons and holes is represented in the following forms, respectively:

    Ne=p,nef(εne(p,F));ρe=NeV=meπ2Lneεne(F)dεeΘ(εeεne(F))1+exp[(εeμ)/(εeμ)kBTkBT] (14)
    Nh=p,nhf(εnh(p,F));ρh=NhV=mhπ2Lnhεnh(F)dεhΘ(εhεnh(F))1+exp[(μεh)/(μεh)kBTkBT] (15)

    Here and below, the index e corresponds to the electronic states, and h to states of the holes, Ne(h) is the number of particles, ρe(h) is the carrier concentration, and V=SL is the volume of the system for area S and thickness L, and Θ(x) is the Heaviside step function.

    Let us now turn to the calculation of the statistical characteristics of the quasi-two-dimensional sample in question.


    3.1. Chemical potential and concentration of charge carriers in the film in the presence of external field

    We will consider an intrinsic semiconductor with a standard parabolic dispersion law for carriers of the form Eq 12 and confine ourselves to the framework of a simple two-band model.In this case, the splitting of the valence band into bands of heavy (hh) and light (lh) holes will be taken into account. It takes place in the semiconductors cited above, including in InAs.


    3.1.1. Only the first sub-band of size quantization is filled

    As the energy reference point we take the bottom of the first electronic sub-band of size quantization. Then, if only the first sub-band is filled, for the corresponding distribution functions fe,h we have:

    fe(εe)=[1+expεeμkBT]1;εe=p22me (16)
    fh(εh)=[1+expμεhkBT]1;εh=ε1,e(F)ε(0)gε1,h(F)p22mh (17)

    Here ε(0)g is the width of band gap of a bulk sample.

    We substitute the distribution functions Eqs 16 and 17 into Eqs 14 and 15, and use after that the electroneutrality condition of the sample:

    ρe=ρlh+ρhh;(ρe=mekBTπ2Lln[1+eμ1/μ1kBTkBT];ρv=mvkBTπ2Lln[1+exp(μ1+ε(0)g+ε1,e(F)+ε1,v(F)kBT)]) (18)

    Here v=lh;hh, and the index 1 for a chemical potential means that only the first sub-band of size quantization is filled. To determine the chemical potential of charge carriers of a film when only the first sub-band is filled (μ1(T,F)), we obtain the following equation:

    meln[1+expμ1kBT]=mlhln[1+exp(μ1+ε1,e(F)+ε1,lh(F)+ε(0)gkBT)]+mhhln[1+exp(μ1+ε1,e(F)+ε1,hh(F)+ε(0)gkBT)] (19)

    Table 1 lists the values of the chemical potential (in dimensionless units y1=μ1(T,F)/μ1(T,F)ε(0)gε(0)g) of a narrow-gap InAs semiconductor film for different values of the sample temperature and external field at a film thickness L = 10 nm.

    Table 1. The values of the chemical potential μ1(T,F) of the InAs (L = 10 nm) film at different values of the sample temperature and the external field strength, when only the first size-quantized sub-band is filled.
    T (K)
    0 50 100 300 600 1000
    y1=μ1/μ1ε(0)gε(0)g η=eFL/eFLε1,eε1,e=0 0.7479 0.7271 0.7095 0.6399 0.5341 0.3819
    η=eFL/eFLε1,eε1,e=10 2.1298 2.1033 2.0859 2.0167 1.9122 1.7714
    η=eFL/eFLε1,eε1,e=30 3.8288 3.8146 3.7996 3.7398 3.6633 3.5226
    η=eFL/eFLε1,eε1,e=40 4.5912 4.5677 4.5503 4.4807 4.3762 4.2356
    η=eFL/eFLε1,eε1,e=50 5.4951 5.4658 5.4346 5.3104 5.1418 4.8909
     | Show Table
    DownLoad: CSV

    The results are obtained by numerical solution of Eq 19.

    As we see, the results of solving Eq 19 clearly show the following:

    (a) The chemical potential of the film, when only the first sub-band is filled, in the absence of the field takes its largest absolute negative value at absolute zero.With the increase in temperature μ1 decreasing in modulus, approaches linearly to the zero energy reference level, i.e., to the bottom of the first film sub-band of the conduction band.

    (b) Under the influence of an external field the maximum modular value of the chemical potential μ1(T,F) increases. This is due to the effective broadening of band gap of the sample under the influence of an external field.

    (с) At a fixed value of the temperature, the value of the chemical potential is the greater in magnitude, the greater is the intensity of the external field.

    (d) In the presence of a field, the chemical potential also increases linearly with increasing the temperature, again approaching the bottom of the conduction band.

    Figure 1 shows the graphical dependence y1(T)=μ1(T,F)/μ1(T,F)ε(0)gε(0)g for different values of the external field η=eFL/eFLε1,cε1,c.

    Figure 1. The dependence y1(T)=μ1(T,F)/μ1(T,F)ε(0)gε(0)g for values of the external field η=eFL/eFLε1,cε1,c=0;10;30;50;(x=kBT/kBTε(0)gε(0)g).

    The same results can also be obtained analytically.For this it is necessary to take into account that over the entire temperature range of the existence of the InAs crystal it can be regarded as a nondegenerate semiconductor with great accuracy. Using the approximate formula

    ln(1+x)|x1x (20)

    from Eq 19 we arrive to the following analytical expression for the chemical potential:

    μ1(T,F)=ε(0)g+ε1,e(F)+ε1,hh(F)2+kBT2{lnmhhme+mlhmhhexp[ε1,lh(F)ε1,hh(F)kBT]}ε(0)g+ε1,e(F)+ε1,hh(F)2+kBT2lnmhhme (21)

    It is easy to see that for the film thickness L = 10 nm, the results in Table 1 and, accordingly, the graphs of Figure 1 are also obtained from Eq 21. In the intrinsic semiconductor, as is known, the chemical potential of electrons and holes is not a function of the concentration (density) of charge carriers.

    However, knowing the chemical potential of the system, it is possible to determine the dependence of the carrier concentration (density) on temperature and in this case, also the dependence on the magnitude of the external field.

    Indeed, substituting the value of μ1(T,F) from Eq 21 into Eq 18, we obtain for the concentration of charge carriers:

    ρ1,e=mekBTπ2Lln{1+(mhhme)1/122exp[ε(0)g+ε1,e(F)+ε1,hh(F)2kBT]}ρ1,lh=mlhkBTπ2Lln{1+(memhh)1/122exp[ε(0)g+ε1,e(F)+ε1,hh(F)2kBTε1,lh(F)ε1,hh(F)kBT]}ρ1,hh=mhhkBTπ2Lln{1+(memhh)1/122exp[ε(0)g+ε1,e(F)+ε1,hh(F)2kBT]}} (22)

    As we see, under the considered conditions and in the approximation Eq 20, the contribution of light holes to the total concentration of positive charge carriers is close to zero, when only the first sub-band of size quantization is filled.

    So, if special accuracy is not required, then in calculations one can put:

    ρ1,lh0,ρ1,e=ρ1,hhmemhhπ2LkBTexp[ε(0)g+ε1,e(F)+ε1,hh(F)2kBT] (23)

    3.1.2. The first and second sub-bands are filled

    To determine the chemical potential (μ2) of the film, when the first and second sub-bands are filled, we now arrive to the following equation (by using Eqs 14, 15 and 18):

    meln{[1+expμ2kBT][1+expμ2ε2,e(F)+ε1,e(F)kBT]}=mlhln{[1+exp(μ2+ε1,e(F)+ε1,lh(F)+ε(0)gkBT)][1+exp(μ2+ε1,e(F)+ε2,lh(F)+ε(0)gkBT)]}+mhhln{[1+exp(μ2+ε1,e(F)+ε1,hh(F)+ε(0)gkBT)][1+exp(μ2+ε1,e(F)+ε2,hh(F)+ε(0)gkBT)]} (24)

    The numerical solution of this equation leads for to the μ2(T,F) results that practically coincide with those from Table 1. If we again take advantage of the fact that the semiconductor under consideration is nondegenerate, then we obtain the following equation for μ2(T,F):

    μ2(T,F)=ε(0)g+ε1,e(F)+ε1,hh(F)2+kBT2lnmhhme+kBT2ln[1+exp(ε2,hh(F)ε1,hh(F)kBT)+mlhmhhexp(ε1,lh(F)ε1,hh(F)kBT)] (25)

    The third term in square brackets can simply be neglected because of its apparent smallness with respect to unity, even in the absence of an external field.

    The remaining exponent in the last term makes sense only when F = 0. In the presence of a strong field, this term also tends to zero, i.e., we arrive at the following result:

    μ2(T,F)μ1(T,F) (26)

    For carriers' concentrations in the film, when the first and second sub-bands are filled we obtain, respectively:

    ρtot,e=ρ1,e+ρ2,e;ρ2,e=mekBTπ2Lln{1+mhhmeexp[εg,1(F)2kBTε2,c(F)ε1,c(F)kBT]}ρtot,lh=ρ1,lh+ρ2,lh;ρ2,lh=mlhkBTπ2Lln{1+memhhexp[εg,1(F)2kBTε2,lh(F)ε1,hh(F)kBT]}ρtot,hh=ρ1,hh+ρ2,hh;ρ2,hh=mhhkBTπ2Lln{1+memhhexp[εg,1(F)2kBTε2,hh(F)ε1,hh(F)kBT]}εg,1(F)=ε(0)g+ε1,e(F)+ε1,hh(F)} (27)

    Comparison of the expressions for ρ1 and ρ2 shows that the concentration of each type of carrier in the second sub-band of size quantization is less than the corresponding concentration in the first sub-band.


    3.2. Energy and heat capacity of charge carriers in a film in the presence of an external field

    For the energy of the system of charge carriers we can write in general form:

    E=p,nεn(p,F)f(εn(p,F)) (28)

    Here, the particle energy εn(p,F) and the distribution function f(εn(p,F)) are given by Eqs 12 and 13, respectively.

    Let us proceed to a concrete calculation of the energy and heat capacity of charge carriers in a film in the presence of an external field.


    3.2.1. Only the first sub-band is filled

    It was shown above Eq 22 that the concentration of light holes is negligible compared to the concentrations of electrons and heavy holes. Accordingly, we will confine ourselves here to calculating the energy and the heat capacity only for electrons and heavy holes. From Eq 28, taking into account Eqs 7, 8, and 14–17, for the energy (E) and heat capacity (C) of these carriers, we obtain, respectively:

    E1E1,e+E1,hh;E1,e=(mhhme)1/122(kBT)2Sπ2exp[εg,1(F)2kBT]E1,hh=(mhhme)1/122(kBT)Sπ2[εg,1(F)+kBT]exp[εg,1(F)2kBT]} (29)
    C=(ET)V;C1C1,e+C1,hh;C1,e=(mhhme)1/122kBSπ2[2kBT+εg,1(F)2]exp[εg,1(F)2kBT]C1,hh=(mhhme)1/122kBSπ2{εg,1(F)+2kBT+εg,1(F)[εg,1(F)+kBT]2kBT}exp[εg,1(F)2kBT]} (30)

    3.2.2. The first and second sub-bands are filled

    Using the value of μ2(T,F) from the Eqs 25 and 26, for the energy and heat capacity of the charge carriers from the second sub-band of size quantization, we get:

    E2,e=(mhhme)1/122(kBT)Sπ2[ε2,e(F)ε1,e(F)+kBT]exp[εg,1(F)2kBTε2,e(F)ε1,e(F)kBT]E2,hh=(mhhme)1/122(kBT)Sπ2[εg,1(F)+ε2,hh(F)ε1,hh(F)+kBT]×exp[εg,1(F)2kBTε2,hh(F)ε1,hh(F)kBT]} (31)
    C2,e=(mhhme)1/122SkBπ2{[ε2,e(F)ε1,e(F)+kBT][εg,1(F)+2ε2,e(F)2ε1,e(F)]2kBT+[ε2,e(F)ε1,e(F)+2kBT]}×exp[εg,1(F)2kBTε2,e(F)ε1,e(F)kBT]C2,hh=(mhhme)1/122SkBπ2{[εg,1(F)+Δ12,hh(F)+kBT][εg,1(F)+2Δ12,hh(F)]2kBT+[εg,1(F)+Δ12,hh(F)+2kBT]}×exp[εg,1(F)2kBTΔ12,hh(F)kBT];Δ12,hh(F)=ε2,hh(F)ε1,hh(F)} (32)

    The total energy and total heat capacity of electrons and holes in the sample will be expressed by the sum of Eqs 29–32, respectively. Let us discuss now the results obtained in this paper.


    4. Results and discussions

    In the quantized film, as it is well known, due to the quantum size effect, an effective broadening of band gap of the sample takes place:

    ε(0)gεg(L)=ε(0)g+ε1,e(L)+ε1,v(L);Δg(L)=εg(L)ε(0)g=ε1,e(L)+ε1,v(L) (33)

    For a fixed film thickness, under the influence of the external field, there is an additional broadening of band gap of the sample with respect to its previous broadening.

    As in the semiconductor InAs, mhhmlh, from the valence band only heavy holes will contribute to this broadening:

    εg(L,F)εg(L,F=0)Δg(F)=[ε1,c(F)ε1,c(L)]+[ε1,hh(F)ε1,hh(F)] (34)

    According to Eqs 7, 8, 10 and 11 with the increase in the external field, this broadening increases according to the law у~x2/3:

    Δg(F)F2/233 (35)

    As we can see, the action of the external field actually substitutes the process of reducing the thickness of the film. This process naturally has a reasonable technological limit. By superimposing an external field, it becomes possible for a given thickness of the film to vary the energy of the quantized motion of the charge carriers in a controlled manner and, accordingly, to vary the width of the band gap of the sample in a controlled manner.

    The chemical potential of a system of electrons and holes depends explicitly on the external field.

    According to the Eqs 21, 25 and 26 remaining negative in sign, with the increase in the field it increases in absolute value (F2/233).

    Both in the absence and in the presence of an external field, the chemical potential of the electrons and holes decrease by absolute value with the increase in the temperature. In the system under consideration, the contribution of light holes in the statistical characteristics of electrons and holes is negligible.

    This is explained by the fact that in the system under consideration, between the first levels of heavy and light holes, there is a significant energy gap (Δlh,hh(F=0)=0.1359eV) and, in the presence of size quantization, light holes do not participate in the formation of a new forbidden band. Physically, this is quite understandable: the particles tend to occupy the lowest possible energy states, and such are the levels of the size quantization of heavy holes in this case.

    We note that in the case of a bulk sample, where the energies of the sub-bands of heavy and light holes in the center of the Brillouin zone coincide, the secondary role of light holes is not so strongly pronounced.

    For example, for the chemical potential of a bulk InAs sample, we obtain:

    μ3Dε(0)g2+3kBT4lnmhhme+kBT2ln[1+(mlhmhh)3/322] (36)

    Hence it is clear that in the 3D case the contribution of light holes is not as small as in the quasi-two-dimensional case.

    It follows from Eqs 22 and 27 that the concentration of free charge carriers increases with the increase in temperature. At fixed values of the temperature and the external field, the carriers' concentration is maximal for the first sub-band. As the sub-band number increases, the concentration decreases. At a fixed temperature, an increase in the external field leads to a decrease in the concentration of free carriers. This is due to an increase in the width of the band gap under the action of the external field. This makes it possible, by changing the external field, to regulate the concentration of free carriers in the quasi-two-dimensional system under consideration. I.e., with the aid of an external field, we are able to control the degree of population of the sub-bands of the quantized motion of free charge carriers. Since the energy and heat capacity of a system are determined by the number of particles, the concentration behavior mentioned above also determines the behavior of the energy and heat capacity of the system of charge carriers in the film in the presence of an external field.

    The contribution of the second sub-band to the energy and heat capacity of electrons and holes is much less (respectively, exp[ε2,e(F)ε1,e(F)2kBT] and exp[ε2,hh(F)ε1,hh(F)2kBT] times) than the contribution of the first sub-band.

    As we see, the result of the influence of an external strong electric field on the statistical characteristics of a quantum semiconductor film is determined, mainly, by the influence of this field on the first film sub-band.

    In particular, as can be seen from Eqs 29–32, the decisive factor for determining the values of the energy and heat capacity of an electron-hole subsystem in a film is the factor exp[εg(F,L)/εg(F,L)2kBT2kBT]. Here εg(F,L) is the value of the band gap of the sample, taking into account its effective broadening under the action of size quantization and the external field (see Eqs 7, 8, 33 and 34). It is clear that with increasing of εg(F,L) the values of the energy and heat capacity at a given fixed temperature will decrease.

    The dependence of the energy of the electron-hole subsystem for the bulk InAs semiconductor (Figure 2a) and for the InAs film (L = 10 nm) in the temperature range from 0 to 600 K for the different values of the external field is shown (Figures 2b2d).

    Figure 2. Curve of the function (a) g1(x)=4.4x3/322(1+3x)exp(1/12x2x); (bulk sample: x=kBT/kBTε(0)gε(0)g;x[0,0.1458]). (b) g2(x)=1.98x(1+2x)exp(1/12x2x); (InAs film, L = 10 nm, x=kBT/kBTεg,1εg,1;x[0,0.097], the external field is absent). (c) g3(x)=14.67x(1+2x)exp(1/12x2x); (InAs film, L = 10 nm, the external field is present: η=10;x=kBT/kBTεg,1εg,1(F);x[0,0.0344]). (d) g4(x)=48.87x(1+2x)exp(1/12x2x); (InAs film, L = 10 nm, the external field is present: η=30;x=kBT/kBTεg,1εg,1(F);x[0,0.019]).

    A similar dependence also holds for the heat capacity of the electron-hole sub-system in a semiconductor film.

    In the end, we note that all the calculations and conclusions that have been made are also applicable to semiconductor quasi-two-dimensional structures noted at the beginning of this article.


    5. Conclusions

    In terms of the effect on single-electron states, the presence of an external field makes it possible, for a fixed film thickness and a certain temperature, to obtain for the charge carriers energies unattainable by reducing the film thickness. Simultaneously, by using an external field, the position of the size quantization levels can be moved in a controlled manner.

    In terms of the effect on the thermodynamic properties of the electron-hole subsystem, the presence of an external field makes it possible to vary the width of the forbidden band of the sample in a controlled manner. That is, an external strong electrostatic field makes it possible to control the values of such fundamental quantities of a system of free charge carriers in a film as their concentration, energy and heat capacity.


    Acknowledgments

    The author wants to express their gratitude to David Hayrapetyan and Vram Mughnetsyan for helping in numerical calculations.


    Conflict of interest

    Тhere is no conflict of interest




    [1] A. Hamilton, F. Haghpanah, M. Hasso-Agopsowicz, I. Frost, G. Lin, E. Schueller, et al., Modeling of malaria vaccine effectiveness on disease burden and drug resistance in 42 African countries, Commun. Med., 3 (2023), 144. https://doi.org/10.1038/s43856-023-00373-y doi: 10.1038/s43856-023-00373-y
    [2] J. Muwanguzi, G. Henriques, P. Sawa, T. Bousema, C. J. Sutherland, K. B. Beshir, Lack of K13 mutations in Plasmodium falciparum persisting after artemisinin combination therapy treatment of Kenyan children, Malar. J., 15 (2016), 1–6. https://doi.org/10.1186/s12936-016-1095-y doi: 10.1186/s12936-016-1095-y
    [3] P. Patel, P. K. Bharti, D. Bansal, N. A. Ali, R. K. Raman, P. K. Mohapatra, et al., Prevalence of mutations linked to antimalarial resistance in Plasmodium falciparum from Chhattisgarh, Central India: A malaria elimination point of view, Sci. Rep., 7 (2017), 16690. https://doi.org/10.1038/s41598-017-16866-5 doi: 10.1038/s41598-017-16866-5
    [4] C. W. Tarama, H. Soré, M. Siribié, S. Débé, R. Kinda, A. Ganou, et al., Plasmodium falciparum drug resistance-associated mutations in isolates from children living in endemic areas of Burkina Faso, Malar. J., 22 (2023), 213. https://doi.org/10.1186/s12936-023-04645-9 doi: 10.1186/s12936-023-04645-9
    [5] A. G. B. Tuedom, E. M. Sarah-Matio, C. E. E. Moukoko, B. L. Feufack-Donfack, C. N. Maffo, A. N. Bayibeki, et al., Antimalarial drug resistance in the Central and Adamawa regions of Cameroon: Prevalence of mutations in P. falciparum crt, Pfmdr1, Pfdhfr and Pfdhps genes, PLoS One, 16 (2021), e0256343. https://doi.org/10.1371/journal.pone.0256343 doi: 10.1371/journal.pone.0256343
    [6] M. I. Veiga, S. K. Dhingra, P. P. Henrich, J. Straimer, N. Gnädig, A. C. Uhlemann, et al., Globally prevalent PfMDR1 mutations modulate Plasmodium falciparum susceptibility to artemisinin-based combination therapies, Nat. Commun., 7 (2016), 11553. https://doi.org/10.1038/ncomms11553 doi: 10.1038/ncomms11553
    [7] World Health Organization, Malaria Microscopy Quality Assurance Manual, World Health Organization, Geneva, 2nd edition, 2016.
    [8] T. D. Nguyen, T. N.-A. Tran, D. M. Parker, N. J. White, M. F. Boni, Antimalarial mass drug administration in large populations and the evolution of drug resistance, PLoS Global Public Health, 3 (2023), e0002200. https://doi.org/10.1371/journal.pgph.0002200 doi: 10.1371/journal.pgph.0002200
    [9] J. C. Koella, R. Antia. Epidemiological models for the spread of anti-malarial resistance, Malar. J., 2 (2003), 1–11. https://doi.org/10.1186/1475-2875-2-3 doi: 10.1186/1475-2875-2-3
    [10] P. Pemberton-Ross, N. Chitnis, E. Pothin, T. A. Smith, A stochastic model for the probability of malaria extinction by mass drug administration, Malar. J., 16 (2017), 1–9. https://doi.org/10.1186/s12936-017-2010-x doi: 10.1186/s12936-017-2010-x
    [11] J. M. Tchuenche, C. Chiyaka, D. Chan, A. Matthews, G. Mayer, A mathematical model for antimalarial drug resistance, Math. Med. Biol. J. IMA, 28 (2011), 335–355. https://doi.org/10.1093/imammb/dqq017 doi: 10.1093/imammb/dqq017
    [12] D. D. Laishram, P. L. Sutton, N. Nanda, V. L. Sharma, R. C. Sobti, J. M. Carlton, et al., The complexities of malaria disease manifestations with a focus on asymptomatic malaria, Malar. J., 11 (2012), 1–15. https://doi.org/10.1186/1475-2875-11-29 doi: 10.1186/1475-2875-11-29
    [13] R. Telfils, A. Y. Dossou, A. Djènontin, E. Adimi, R. Akoho, J. Bailly, et al., Dynamics of submicroscopic and microscopic asymptomatic malaria infection and associated factors: A longitudinal study in South Benin, PLoS One, 19 (2024), e0311217. https://doi.org/10.1371/journal.pone.0311217 doi: 10.1371/journal.pone.0311217
    [14] L. S. Pontryagin, V. G. Boltyanskii, Mathematical Theory of Optimal Processes, John Wiley & Sons, 1st edition, 1962. https://doi.org/10.1201/9780203749319
    [15] N. Chitnis, J. M. Hyman, J. M. Cushing, Determining important parameters in the spread of malaria through the sensitivity analysis of a mathematical model, Bull. Math. Biol., 70 (2008), 1272–1296. https://doi.org/10.1007/s11538-008-9299-0 doi: 10.1007/s11538-008-9299-0
    [16] A. Mahamar, M. J. Smit, K. Sanogo, Y. Sinaba, S. M. Niambele, A. Sacko, et al., Artemether-lumefantrine with or without single-dose primaquine and sulfadoxine-pyrimethamine plus amodiaquine with or without single-dose tafenoquine to reduce Plasmodium falciparum transmission: A phase 2, single-blind, randomised clinical trial in Ouelessebougou, Mali, Lancet Microbe, 5 (2024), 633–644. https://doi.org/10.1016/S2666-5247(24)00023-5 doi: 10.1016/S2666-5247(24)00023-5
    [17] L. Cirera, B. Galatas, S. Alonso, K. Paaijmans, M. Mamuquele, H. Martí-Soler, et al., Moving towards malaria elimination in southern Mozambique: Cost and cost-effectiveness of mass drug administration combined with intensified malaria control, PLoS One, 15 (2020), e0235631. https://doi.org/10.1371/journal.pone.0235631 doi: 10.1371/journal.pone.0235631
    [18] S. S. Kyaw, G. Delmas, T. L. Drake, O. Celhay, W. Pan-Ngum, S. Pukrittayakamee, et al., Estimating the programmatic cost of targeted mass drug administration for malaria in Myanmar, BMC Public Health, 21 (2021), 826. https://doi.org/10.1186/s12889-021-10842-5 doi: 10.1186/s12889-021-10842-5
    [19] J. O. Yukich, C. Scott, K. Silumbe, B. A. Larson, A. Bennett, T. P. Finn, et al., Cost-effectiveness of focal mass drug administration and mass drug administration with dihydroartemisinin–piperaquine for malaria prevention in Southern Province, Zambia: Results of a Community-Randomized Controlled Trial, Am. J. Trop. Med. Hyg., 103(2_Suppl) (2020), 46–53. https://doi.org/10.4269/ajtmh.19-0661 doi: 10.4269/ajtmh.19-0661
    [20] J. K. Hale, J. P. LaSalle, M. Slemrod, Theory of a general class of dissipative processes, J. Math. Anal. Appl., 39 (1972), 177–191.
    [21] O. A. Ladyzhenskaya, On the determination of minimal global attractors for the Navier-Stokes and other partial differential equations, Russ. Math. Surv., 42 (1987), 25–60. https://doi.org/10.1070/RM1987v042n06ABEH001503 doi: 10.1070/RM1987v042n06ABEH001503
    [22] H. R. Thieme, Spectral bound and reproduction number for infinite-dimensional population structure and time heterogeneity, SIAM J. Appl. Math., 70 (2009), 188–211. https://doi.org/10.1137/080732870 doi: 10.1137/080732870
    [23] P. Magal, Compact attractors for time-periodic age-structured population models, Electron. J. Differ. Equations (EJDE), 2001 (2001), 1–35.
    [24] A. Pazy, Semigroups of Linear Operators and Applications to Partial Differential Equations, Springer-Verlag Berlin and Heidelberg GmbH & Co. K, New York Berlin Heidelberg, 1983.
    [25] W. Fleming, R. Rishel, Deterministic and Stochastic Optimal Control, Springer New York, New York, NY, 1975.
    [26] S. Lenhart, J. T. Workman, Optimal Control Applied to Biological Models, Chapman and Hall/CRC, New York, May 2007. https://doi.org/10.1201/9781420011418
    [27] J. Macki, A. Strauss, Introduction to Optimal Control Theory, Springer, New York Berlin Heidelberg, corr print edition edition, 1981.
    [28] Q. Richard, M. Choisy, T. Lefèvre, R. Djidjou-Demasse, Human-vector malaria transmission model structured by age, time since infection and waning immunity, Nonlinear Anal. Real World Appl., 63 (2022), 103393. https://doi.org/10.1016/j.nonrwa.2021.103393 doi: 10.1016/j.nonrwa.2021.103393
    [29] Q. Richard, M. Choisy, T. Lefèvre, R. Djidjou-Demasse, On the necessity of accounting for age structure in human malaria transmission modeling, Math. Biosci., 378 (2024), 109319. https://doi.org/10.1016/j.mbs.2024.109319 doi: 10.1016/j.mbs.2024.109319
    [30] B. M. Somé, E. Guissou, D. F. Da, Q. Richard, M. Choisy, K. B. Yameogo, et al., Mosquito ageing modulates the development, virulence and transmission potential of pathogens, Proc. R. Soc. B Biol. Sci., 291 (2024), 20232097. https://doi.org/10.1098/rspb.2023.2097 doi: 10.1098/rspb.2023.2097
    [31] T. Cazenave, A. Haraux, An Introduction to Semilinear Evolution Equations, Oxford University Press, 1998.
    [32] J. K. Hale, Asympstotic Behavior of Dissipative Systems, American Mathematical Society, 1988.
    [33] G. F. Webb, Compactness of bounded trajectories of dynamical systems in infinite dimensional spaces, Proc. R. Soc. Edinburgh Sec. A: Math., 84 (1979), 19–33. https://doi.org/10.1017/S0308210500016930 doi: 10.1017/S0308210500016930
    [34] P. Magal, S. Ruan, Theory and Applications of Abstract Semilinear Cauchy Problems, Manhattan, NY, USA, Springer International Publishing, 2018. https://doi.org/10.1007/978-3-030-01506-0
  • Reader Comments
  • © 2025 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(256) PDF downloads(33) Cited by(0)

Article outline

Figures and Tables

Figures(9)  /  Tables(2)

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog