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
[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.
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.
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)
L2≪a2ex | (1) |
(b) The distance between the levels of the size quantization
Δεn,n+1≫kBT | (2) |
Here
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)
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
U(z,F)={∞;z≤0,z≥L|e|Fz;0≤z≤L | (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+p2y2m∗e+π2ℏ2n22m∗eL2≡ε(px,py)+ε1n2;n=1,2,... | (4) |
Here
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>>ε1=π2ℏ2/π2ℏ22m∗eL22m∗eL2 | (5) |
It is shown in Ref. [42], that the strong-field condition with sufficient accuracy can be written as follows
|e|FL≥10ε1 | (6) |
In the case of an InAs film L = 10 nm, we have for
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(ℏ2e2F22m∗e)1/133=1.093ε1,e(|e|FLε1,e)2/233≡1.093ε1,eη2/233 | (7) |
ε2,e(F)≃4.109(ℏ2e2F22m∗e)1/133=1.931ε1,e(|e|FLε1,e)2/233≡1.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)=√2m∗z0∫0√εl(F)−|e|Fzdz;(z0=εl(F)/εl(F)|e|F|e|F); | (9) |
Hence we obtain:
εl(F)=[32πℏ|e|F√2m∗(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
The melting temperature of the InAs crystal is equal to
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.
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)=→p22m∗e+εn(F)=ε(→p)+εn(F)≡ε;→p=→p(px,py,0) | (12) |
Here the energy
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
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=m∗eπℏ2L∑ne∞∫εne(F)dεeΘ(εe−εne(F))1+exp[(εe−μ)/(εe−μ)kBTkBT] | (14) |
Nh=∑→p,nhf(εnh(→p,F));ρh=NhV=m∗hπℏ2L∑nh∞∫ε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,
Let us now turn to the calculation of the statistical characteristics of the quasi-two-dimensional sample in question.
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.
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(εe)=[1+expεe−μkBT]−1;εe=p22m∗e | (16) |
fh(εh)=[1+expμ−εhkBT]−1;εh=−ε1,e(F)−ε(0)g−ε1,h(F)−p22m∗h | (17) |
Here
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
m∗eln[1+expμ1kBT]=m∗lhln[1+exp(−μ1+ε1,e(F)+ε1,lh(F)+ε(0)gkBT)]+m∗hhln[1+exp(−μ1+ε1,e(F)+ε1,hh(F)+ε(0)gkBT)] | (19) |
Table 1 lists the values of the chemical potential (in dimensionless units
T (K) | |||||||
0 | 50 | 100 | 300 | 600 | 1000 | ||
0.7479 | 0.7271 | 0.7095 | 0.6399 | 0.5341 | 0.3819 | ||
2.1298 | 2.1033 | 2.0859 | 2.0167 | 1.9122 | 1.7714 | ||
3.8288 | 3.8146 | 3.7996 | 3.7398 | 3.6633 | 3.5226 | ||
4.5912 | 4.5677 | 4.5503 | 4.4807 | 4.3762 | 4.2356 | ||
5.4951 | 5.4658 | 5.4346 | 5.3104 | 5.1418 | 4.8909 |
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
(b) Under the influence of an external field the maximum modular value of the chemical potential
(с) 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
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)|x≪1∼x | (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,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,lh≈0,ρ1,e=ρ1,hh≈√memhhπℏ2LkBTexp[−ε(0)g+ε1,e(F)+ε1,hh(F)2kBT] | (23) |
To determine the chemical potential
m∗eln{[1+expμ2kBT][1+expμ2−ε2,e(F)+ε1,e(F)kBT]}=m∗lhln{[1+exp(−μ2+ε1,e(F)+ε1,lh(F)+ε(0)gkBT)][1+exp(−μ2+ε1,e(F)+ε2,lh(F)+ε(0)gkBT)]}+m∗hhln{[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)=−ε(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
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
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.
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:
E1∼E1,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=(∂E∂T)V;C1∼C1,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) |
Using the value of
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.
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,
ε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 (
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 (
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,
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
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 2b–2d).
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.
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.
The author wants to express their gratitude to David Hayrapetyan and Vram Mughnetsyan for helping in numerical calculations.
Т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 |
T (K) | |||||||
0 | 50 | 100 | 300 | 600 | 1000 | ||
0.7479 | 0.7271 | 0.7095 | 0.6399 | 0.5341 | 0.3819 | ||
2.1298 | 2.1033 | 2.0859 | 2.0167 | 1.9122 | 1.7714 | ||
3.8288 | 3.8146 | 3.7996 | 3.7398 | 3.6633 | 3.5226 | ||
4.5912 | 4.5677 | 4.5503 | 4.4807 | 4.3762 | 4.2356 | ||
5.4951 | 5.4658 | 5.4346 | 5.3104 | 5.1418 | 4.8909 |
T (K) | |||||||
0 | 50 | 100 | 300 | 600 | 1000 | ||
0.7479 | 0.7271 | 0.7095 | 0.6399 | 0.5341 | 0.3819 | ||
2.1298 | 2.1033 | 2.0859 | 2.0167 | 1.9122 | 1.7714 | ||
3.8288 | 3.8146 | 3.7996 | 3.7398 | 3.6633 | 3.5226 | ||
4.5912 | 4.5677 | 4.5503 | 4.4807 | 4.3762 | 4.2356 | ||
5.4951 | 5.4658 | 5.4346 | 5.3104 | 5.1418 | 4.8909 |