Efficient numerical methods for solving Poisson equation constraint optimal control problems with random coefficient are discussed in this paper. By applying the finite element method and the Monte Carlo approximation, the original optimal control problem is discretized and transformed into an optimization problem. Taking advantage of the separable structures, Algorithm 1 is proposed for solving the problem, where an alternating direction method of multiplier is used. Both computational and storage costs of this algorithm are very high. In order to reduce the computational cost, Algorithm 2 is proposed, where the multi-modes expansion is introduced and applied. Further, for reducing the storage cost, we propose Algorithm 3 based on Algorithm 2. The main idea is that the random term is shifted to the objective functional, which could be computed in advance. Therefore, we only need to solve a deterministic optimization problem, which could reduce all the costs significantly. Moreover, the convergence analyses of the proposed algorithms are established, and numerical simulations are carried out to test the performances of them.
Citation: Xiaowei Pang, Haiming Song, Xiaoshen Wang, Jiachuan Zhang. Efficient numerical methods for elliptic optimal control problems with random coefficient[J]. Electronic Research Archive, 2020, 28(2): 1001-1022. doi: 10.3934/era.2020053
[1] | Iryna Markevich, Tetyana Stara, Larysa Khomenkova, Volodymyr Kushnirenko, Lyudmyla Borkovska . Photoluminescence engineering in polycrystalline ZnO and ZnO-based compounds. AIMS Materials Science, 2016, 3(2): 508-524. doi: 10.3934/matersci.2016.2.508 |
[2] | R.A. Silva, C.O. Soares, R. Afonso, M.D. Carvalho, A.C. Tavares, M.E. Melo Jorge, A. Gomes, M.I. da Silva Pereira, C.M. Rangel . Synthesis and electrocatalytic properties of La0.8Sr0.2FeO3−δ perovskite oxide for oxygen reactions. AIMS Materials Science, 2017, 4(4): 991-1009. doi: 10.3934/matersci.2017.4.991 |
[3] | Imaddin A. Al-Omari, Muna D. Al-Mamari, D.J. Sellmyer . Tuning the giant Magnetocaloric Effect and refrigerant capacity in Gd1–xYxCrO3 (0.0 ≤ x ≤ 0.9) perovskites nanoparticles. AIMS Materials Science, 2022, 9(2): 297-310. doi: 10.3934/matersci.2022018 |
[4] | Larysa Khomenkova, Mykola Baran, Jedrzej Jedrzejewski, Caroline Bonafos, Vincent Paillard, Yevgen Venger, Isaac Balberg, Nadiia Korsunska . Silicon nanocrystals embedded in oxide films grown by magnetron sputtering. AIMS Materials Science, 2016, 3(2): 538-561. doi: 10.3934/matersci.2016.2.538 |
[5] | Yana Fajar Prakasa, Sumari Sumari, Aman Santoso, Muhammad Roy Asrori, Ririn Cahyanti . The performance of radar absorption of MnxFe3–xO4/rGO nanocomposites prepared from iron sand beach and coconut shell waste. AIMS Materials Science, 2023, 10(2): 227-248. doi: 10.3934/matersci.2023013 |
[6] | Etsana Kiros Ashebir, Berhe Tadese Abay, Taame Abraha Berhe . Sustainable A2BⅠBⅢX6 based lead free perovskite solar cells: The challenges and research roadmap for power conversion efficiency improvement. AIMS Materials Science, 2024, 11(4): 712-759. doi: 10.3934/matersci.2024036 |
[7] | Alexandre Lavrov, Elvia Anabela Chavez Panduro, Kamila Gawel, Malin Torsæter . Electrophoresis-induced structural changes at cement-steel interface. AIMS Materials Science, 2018, 5(3): 414-421. doi: 10.3934/matersci.2018.3.414 |
[8] | Jamal Alnofiay, Ahmed Al-Shahrie, Elsayed Shalaan . Green synthesis of high-performance gallium oxide supercapacitor: A path to outstanding energy density. AIMS Materials Science, 2024, 11(6): 1065-1082. doi: 10.3934/matersci.2024051 |
[9] | Finn Zahari, Mirko Hansen, Thomas Mussenbrock, Martin Ziegler, Hermann Kohlstedt . Pattern recognition with TiOx-based memristive devices. AIMS Materials Science, 2015, 2(3): 203-216. doi: 10.3934/matersci.2015.3.203 |
[10] | M. F. Achoi, S. Kato, N. Kishi, T. Soga . Improved photovoltaic properties of ((CH3NH3)1-xCsx)3Bi2I9: (x = 0-1.0) hybrid perovskite solar cells via a hot immersion method. AIMS Materials Science, 2024, 11(4): 605-619. doi: 10.3934/matersci.2024031 |
Efficient numerical methods for solving Poisson equation constraint optimal control problems with random coefficient are discussed in this paper. By applying the finite element method and the Monte Carlo approximation, the original optimal control problem is discretized and transformed into an optimization problem. Taking advantage of the separable structures, Algorithm 1 is proposed for solving the problem, where an alternating direction method of multiplier is used. Both computational and storage costs of this algorithm are very high. In order to reduce the computational cost, Algorithm 2 is proposed, where the multi-modes expansion is introduced and applied. Further, for reducing the storage cost, we propose Algorithm 3 based on Algorithm 2. The main idea is that the random term is shifted to the objective functional, which could be computed in advance. Therefore, we only need to solve a deterministic optimization problem, which could reduce all the costs significantly. Moreover, the convergence analyses of the proposed algorithms are established, and numerical simulations are carried out to test the performances of them.
Antimonide based quantum wells and nanostructures continue to attract considerable interest due to their potential for accessing technologically important applications in the spectral range beyond 1.5 µm. In particular GaInSb/GaSb offers a route towards the development of sources and detectors for use in remote gas sensing instrumentation, medical diagnostics and night vision technologies [1,2]. GaSb-based diode lasers, LEDs and photodetectors have been developed with excellent performance characteristics operating at many of the key wavelengths including methane [3], CO2 [4] and CO. However, the lack of semi-insulating GaSb substrates precludes the development of monolithic focal plane arrays and more complex devices, while the relatively high cost of GaSb wafers makes large volume applications cost-sensitive. Consequently, it is worthwhile investigating alternative growth methods for GaInSb based device structures on inexpensive GaAs substrates, which are readily available as semi-insulating wafers and for which there exists a mature processing technology. Recently, the growth of GaSb on GaAs has been successfully demonstrated using the interface misfit (IMF) [5,6,7] technique, which accommodates the lattice mismatch by forming lateral misfit dislocations along the GaSb-GaAs interface rather than threading dislocations in the GaSb bulk. However, samples grown using this molecular beam epitaxy (MBE) technique have shown reduced optical quality compared with homoepitaxial material since not all the dislocations are accommodated at the heterointerface and the remaining threading dislocations can still be problematic. In the present work we report on the MBE growth of AlGaSb as a ternary IMF metamorphic buffer layer on GaAs and study the properties of GaInSb quantum wells grown on this novel relaxed quasi substrate. AlGaSb conveniently provides substantial offsets in both the conduction and valence bands and results in improved photoluminescence efficiency from AlGaSb/InGaSb MQWs.
The Ga1-xInxSb MQW structures were grown on a (100) orientated n-GaSb substrate, or in the case of the IMF samples a (100) n-GaAs substrate, by solid source MBE using a VG-V80H reactor fitted with a Sb valved cracker cell. The GaSb MQW sample (QA216) in Figure1(a) was grown at a substrate temperature of 510 ℃ using growth rates of 0.65 ML/s for the GaSb and 1 ML/s for the Ga1-xInxSb QW respectively. Once a buffer layer of GaSb was grown, GaInSb QWs with 50 nm GaSb barriers were repeated 5 times. The GaSb IMF (QA276) in Figure 1(b) was grown using these conditions with the GaSb IMF growth based on the method used previously for GaSb IMF on GaAs [8]. The ternary AlGaSb IMF sample (QJ410), in Figure1(c) was grown in a similar manner with a 1μm AlGaSb buffer layer. For the AlGaSb MQW sample a growth rate of 0.46 ML/s with a substrate temperature of 515 ℃ was used. After the IMF growth a 1 μm buffer had been grown, 5x InGaSb QWs with 50 nm Al0.53Ga0.47Sb barriers were grown before finally capping with 100 nm of AlGaSb. Details of the resulting structures are given in Table 1.
QA216 | QA276 | QJ410 | ||||
Buffer layer | GaSb | GaSb IMF | AlGaSb IMF | |||
QW thickness (Lz)/ nm | 4±1 | 4±1 | 7±4 | |||
QW In content (%) | 34 | 32 | 28 | |||
Strain (%) | 2.1 | 2.0 | 1.4 | |||
PL peak | Wavelength (μm) | 1.77 | 1.77 | 1.70 | ||
Energy (eV) | 0.699 | 0.703 | 0.728 | |||
Relative PL intensity (10 K) | 1 | 0.025 | 0.66 | |||
FWHM (meV)at 4 K | 15.4 | 16.7 | 20.5 | |||
Confinement energy (meV) | e- | 49 | 48 | 506 | ||
h+ | 98 | 98 | 345 | |||
Z parameter | 4K | 1.9 | 1.2 | 1.8 | ||
225K | 2.2 | - | 2.3 |
The structural properties of the resulting samples were studied using high resolution x-ray diffraction (XRD) and transmission electron microscopy (TEM), to obtain information about dislocations, In content, thickness and perfection of the MQWs. In order to probe the optical properties of the samples, PL spectroscopy (4-300 K) was carried out using a continuous flow helium cryostat (Oxford Instruments Optistat). Excitation at 514 nm was provided by an Ar+ ion laser giving an adjustable power density of 0.4 to 10 Wcm−2 at the sample surface. The PL emission was detected and analysed using a Bentham M300 monochromator having a grating blazed at 3.5 µm, with a 77 K InSb photodiode detector and lock-in amplifier.
The high resolution XRD spectra obtained from each of the samples is shown in Figure 2. Sample (QA216) which was grown homoepitaxially on GaSb, has clearly defined satellite peaks with a full width half maximum (FWHM) of 78 arcsecs and a GaSb peak with FWHM of 27 arcsecs, indicating high crystalline quality. For sample (QA276) the satellite peaks are less distinct, being reduced in intensity with a FWHM of ~200 arcsecs, while the peak from the GaSb grown IMF layer is broadened with a FWHM of ~220 arcsecs; although this behaviour is consistent with the broadening of quasi substrate peaks seen in other IMF samples [9].
The XRD spectrum of QJ410 with the ternary AlGaSb IMF metamorphic buffer structure shown in Figure 1(c) has no observable satellite peaks and a rather broad AlGaSb “quasi-substrate” broad band with a FWHM of 473 arcsecs. The lack of satellite peaks is normally indicative of a mixed composition layer, or a structure which contains several different lattice parameters. To help elucidate this, TEM images of this sample are shown at different magnifications in Figure 3. The ternary IMF dislocation array is clearly evident at the AlGaSb-GaAs interface in Figure 3(a) and the MQW are shown in the higher magnification image of Figure3(b). Figure 3(c) shows that the overall quality of the metamorphic buffer using the AlGaSb is high with relatively few threading dislocations. Abrupt thickness fluctuations in the QWs are clearly observable under high magnification in Figure 3(b). The target QW thickness was 6 nm but the actual QWs vary in a step-like fashion with thickness from 7 nm to 12 nm across the sample, which is sufficient to prevent the observation of the satellite peaks in the XRD spectrum. The appearance of these abrupt thickness variations is reminiscent of the early stages of quantum dot formation as the strain limit is approached but may also be related to the growth temperature, since 515 ℃ was used throughout to avoid the growth interrupt normally employed to grow InGaSb and AlGaSb layers at their respective optimum temperatures, (480 ℃ and 515 ℃ respectively).
A comparison of the 4 K PL spectra from each of the MQW samples is shown in Figure 4. The relative intensities, peak energy and FWHM are given in Table 1. The peak energies are in reasonable agreement with the calculated values obtained using a Schrödinger solver for a finite quantum well within the effective mass approximation. The corresponding confinement energies for electrons and holes are given in Table 1. The calculations are based on the indium compositions and QW thickness values obtained from the growth parameters, XRD measurements and TEM images. Although the MQW in samples QA216 and QA276 have nominally the same structure, the QA276 GaSb IMF sample has a slightly higher peak energy of 0.703 eV compared with 0.699 eV due to a decrease in Indium content of 1.7% which originates from a change in growth rate when growing on the IMF buffer. Meanwhile, in sample QJ410 the higher band offset provided by the Al0.53Ga0.47Sb produces higher confinement and consequently the peak energy is further increased to 0.729eV.
The FWHM of the samples at 4 K (Table 1) compare favourably with samples grown previously both on native GaSb substrates and on Si [10]. The temperature dependent PL of QJ410 exhibits emission up to room temperature due to the stronger carrier confinement provided by the AlGaSb barriers. The FWHM increases to 39.5 meV as shown in Figure4(b), and an exponential tail is observed which extends up to higher energy which is associated with band filling and thermal broadening in the QWs.
The power dependence of the PL emission from each of the samples was also measured to determine the dominant recombination process based on the PL intensity (L) vs pump power (I) expression; where the exponent Z= 1 corresponds to Shockley-Read-Hall (SRH), Z = 2 corresponds to radiative recombination and Z = 3 corresponds to non-radiative Auger recombination respectively [11]. The corresponding Z parameter are given in Table 1 from which we observe that recombination is predominantly radiative in these samples except for QA276 where SRH recombination is more significant at 10 K. These values are also plotted in Figure 5 showing the much lower value for QA276 which is attributed to defects in the lattice structure whilst showing a nominally constant radiative-dominated process in both QA216 and QJ410. Due to the limitations of this method, differences in Z paramter of ~0.1 can be considered negligible due to the assumptions in the model including the need for a constant number of carriers.
The temperature dependence of the integrated PL emission intensity is shown in Figure 6. The high initial intensity of QA216 at 10 K is attributed to better crystalline perfection. The recombination is predominantly radiative and the thermal quenching is determined by limited carrier (electron) confinement in the MQWs. Sample QA276 quenches at about the same rate as QJ410 but has the lowest 4 K PL intensity which is as expected from the increased threading dislocation density originating from the heteroepitaxial growth mediated by the IMF buffer layer. By comparison the 4 K PL intensity of the AlGaSb ternary IMF sample (QJ410) is about one order higher over the entire temperature range. The increased confinement arising from the AlGaSb barriers helps to maintain PL emission up to room temperature compared to the other two samples which are quenched above ~250 K. This is also supported by the nominally equal recombination processes in QA216 and QJ410 but weaker emission from the smaller barriers of QA216. Both IMF based samples (QA276 and QJ410) are also more susceptible to carrier recombination in the barriers due to the higher levels of threading dislocations and structure defects compared with the homoepitaxial sample QA216.
In summary, we have reported the MBE growth of strained InGaSb/(Al)GaSb MQW structures containing ~30% In on GaSb and GaAs substrates using IMF buffers of both GaSb and ternary AlGaSb. The structural properties of the different samples and the effect on the PL emission efficiency and thermal quenching were compared. The structural perfection and crystallinity of the homoepitaxial MQW was found to be superior and the thermal quenching in this case is determined mainly by the electron-hole confinement in the QW. Transferring this structure directly onto GaAs using a GaSb IMF approach resulted in inferior PL emission intensity at 4 K and thermal quenching of the emission above 250 K, which we attributed to non-radiative SRH recombination within the GaSb barriers. By comparison, it has been shown that when a ternary AlGaSb IMF is used instead, the 4 K PL intensity is considerably recovered and the PL emission persists up to room temperature, due to the improved crystallinity of the ternary IMF and the additional carrier confinement for both electrons and holes arising from the AlGaSb barriers. In each case the PL emission peaks were reconciled with the calculated values. Although further optimisation is required to improve the thickness uniformity of the InGaSb MQW, we have shown that Al0.53Ga0.47Sb/Ga0.7In0.3Sb MQWs grown on GaAs using a ternary Al0.53Ga0.47Sb IMF buffer layer can provide a viable alternative to growth on GaSb substrates to access longer wavelength applications since it facilitates the use of both inexpensive and semi-insulating substrates.
This work was supported by The Centre for Global Eco-Innovation, Oxley Developments Ltd & European Regional Development Fund 2007-13. Thanks are also given to Dr Richard Beanland of Warwick University for the TEM images.
All Authors declare no conflicts of interest in this paper
[1] |
A stochastic collocation method for elliptic partial differential equations with random input data. SIAM Rev (2010) 52: 317-355. ![]() |
[2] |
A stochastic collocation method for elliptic partial differential equations with random input data. SIAM J. Numer. Anal. (2007) 45: 1005-1034. ![]() |
[3] |
Solving elliptic boundary value problems with uncertain coefficients by the finite element method: The stochastic formulation. Comput. Methods. Appl. Mech. Engrg. (2005) 194: 1251-1294. ![]() |
[4] |
Multi-level Monte Carlo finite element method for elliptic PDEs with stochastic coefficients. Numer. Math. (2011) 119: 123-161. ![]() |
[5] |
R. E. Caflisch, Monte-Carlo and quasi-Monte Carlo methods, in Acta Numerica, Acta Numer., 7, Cambridge University Press, Cambridge, 1998, 1–49. doi: 10.1017/S0962492900002804
![]() |
[6] |
complexity analysis of the alternating direction method of multipliers. Sci. China Math. (2019) 62: 795-808. ![]() |
[7] |
An efficient Monte Carlo method for optimal control problems with uncertainty. Comput. Optim. Appl. (2003) 26: 219-230. ![]() |
[8] |
Taylor approximation and variance reduction for PDE-constrained optimal control under uncertainty. J. Comput. Phys. (2019) 385: 163-186. ![]() |
[9] | P. G. Ciarlet, The Finite Element Method for Elliptic Problems, Studies in Mathematics and its Applications, 4, North-Holland Publishing Co., Amsterdam-New York-Oxford, 1978. |
[10] |
An efficient numerical method for acoustic wave scattering in random media. SIAM/ASA J. Uncertain. Quantif. (2015) 3: 790-822. ![]() |
[11] |
A multimodes Monte Carlo finite element method for elliptic partial differential equations with random coefficients. Int. J. Uncertain. Quantif. (2016) 6: 429-443. ![]() |
[12] |
D. Gilbarg and N. S. Trudinger, Elliptic Partial Differential Equations of Second Order, Classics in Mathematics, Springer-Verlag, Berlin, 2001. doi: 10.1007/978-3-642-61798-0
![]() |
[13] |
Error estimates of stochastic optimal Neumann boundary control problems. SIAM J. Numer. Anal. (2011) 49: 1532-1552. ![]() |
[14] |
An alternating direction method of multipliers for the optimization problem constrained with a stationary Maxwell system. Commun. Comput. Phys. (2018) 24: 1435-1454. ![]() |
[15] |
M. Hinze, R. Pinnau, M. Ulbrich and S. Ulbrich, Optimization with PDE Constraints, Mathematical Modelling: Theory and Applications, 23, Springer, New York, 2009. doi: 10.1007/978-1-4020-8839-1
![]() |
[16] |
Finite element approximations of stochastic optimal control problems constrained by stochastic elliptic PDEs. J. Math. Anal. Appl. (2011) 384: 87-103. ![]() |
[17] |
D. P. Kouri, M. Heinkenschloss, D. Ridzal and B. G. van Bloemen Waanders, A trust-region algorithm with adaptive stochastic collocation for PDE optimization under uncertainty, SIAM J. Sci. Comput., 35 (2013), A1847–A1879. doi: 10.1137/120892362
![]() |
[18] |
Analytic regularity and GPC approximation for control problems constrained by linear parametric elliptic and parabolic PDEs. SIAM J. Control Optim. (2013) 51: 2442-2471. ![]() |
[19] |
An efficient alternating direction method of multipliers for optimal control problems constrained by random Helmholtz equations. Numer. Algorithms (2018) 78: 161-191. ![]() |
[20] |
On the Schrödinger equation and the eigenvalue problem. Comm. Math. Phys. (1983) 88: 309-318. ![]() |
[21] |
J. C. De los Reyes, Numerical PDE-Constrained Optimization, SpringerBriefs in Optimization, Springer, Cham, 2015. doi: 10.1007/978-3-319-13395-9
![]() |
[22] |
Stochastic Galerkin method for elliptic SPDEs: A white noise approach. Discrete Contin. Dyn. Syst. Ser. B (2006) 6: 941-955. ![]() |
[23] | Stochastic Galerkin method for constrained optimal control problem governed by an elliptic integro-differential PDE with stochastic coefficients. Int. J. Numer. Anal. Model. (2015) 12: 593-616. |
[24] |
Lower bounds for higher eigenvalues by finite difference methods. Pacific J. Math. (1958) 8: 339-368. ![]() |
[25] |
Fast numerical methods for robust optimal design. Eng. Optim. (2008) 40: 489-504. ![]() |
[26] |
An alternating direction method of multipliers for elliptic equation constrained optimization problem. Sci. China Math. (2017) 60: 361-378. ![]() |
1. | L. Chenini, A. Aissat, M. Halbwax, J.P. Vilcot, Performance simulation of an InGaSb/GaSb based quantum well structure for laser diode applications, 2023, 03759601, 128711, 10.1016/j.physleta.2023.128711 | |
2. | L. S. Lunin, M. L. Lunina, A. S. Pashchenko, A. V. Donskaya, Effect of Concentration Supercooling on the Structure and Properties of GaInAsSbP/GaP Heterostructures, 2025, 19, 1027-4510, 21, 10.1134/S1027451025700041 |
QA216 | QA276 | QJ410 | ||||
Buffer layer | GaSb | GaSb IMF | AlGaSb IMF | |||
QW thickness (Lz)/ nm | 4±1 | 4±1 | 7±4 | |||
QW In content (%) | 34 | 32 | 28 | |||
Strain (%) | 2.1 | 2.0 | 1.4 | |||
PL peak | Wavelength (μm) | 1.77 | 1.77 | 1.70 | ||
Energy (eV) | 0.699 | 0.703 | 0.728 | |||
Relative PL intensity (10 K) | 1 | 0.025 | 0.66 | |||
FWHM (meV)at 4 K | 15.4 | 16.7 | 20.5 | |||
Confinement energy (meV) | e- | 49 | 48 | 506 | ||
h+ | 98 | 98 | 345 | |||
Z parameter | 4K | 1.9 | 1.2 | 1.8 | ||
225K | 2.2 | - | 2.3 |
QA216 | QA276 | QJ410 | ||||
Buffer layer | GaSb | GaSb IMF | AlGaSb IMF | |||
QW thickness (Lz)/ nm | 4±1 | 4±1 | 7±4 | |||
QW In content (%) | 34 | 32 | 28 | |||
Strain (%) | 2.1 | 2.0 | 1.4 | |||
PL peak | Wavelength (μm) | 1.77 | 1.77 | 1.70 | ||
Energy (eV) | 0.699 | 0.703 | 0.728 | |||
Relative PL intensity (10 K) | 1 | 0.025 | 0.66 | |||
FWHM (meV)at 4 K | 15.4 | 16.7 | 20.5 | |||
Confinement energy (meV) | e- | 49 | 48 | 506 | ||
h+ | 98 | 98 | 345 | |||
Z parameter | 4K | 1.9 | 1.2 | 1.8 | ||
225K | 2.2 | - | 2.3 |