Citation: Małgorzata Wawrzyniak-Adamczewska, Małgorzata Wierzbowska, Juan José Meléndez. Effect of graphene substrate on the spectroscopic properties of photovoltaic molecules: role of the in-plane and out-of-plane π-bonds[J]. AIMS Materials Science, 2017, 4(1): 89-101. doi: 10.3934/matersci.2017.1.89
[1] | Andrei Korobeinikov, Aleksei Archibasov, Vladimir Sobolev . Order reduction for an RNA virus evolution model. Mathematical Biosciences and Engineering, 2015, 12(5): 1007-1016. doi: 10.3934/mbe.2015.12.1007 |
[2] | Nara Bobko, Jorge P. Zubelli . A singularly perturbed HIV model with treatment and antigenic variation. Mathematical Biosciences and Engineering, 2015, 12(1): 1-21. doi: 10.3934/mbe.2015.12.1 |
[3] | Helen Moore, Weiqing Gu . A mathematical model for treatment-resistant mutations of HIV. Mathematical Biosciences and Engineering, 2005, 2(2): 363-380. doi: 10.3934/mbe.2005.2.363 |
[4] | M. B. A. Mansour . Computation of traveling wave fronts for a nonlinear diffusion-advection model. Mathematical Biosciences and Engineering, 2009, 6(1): 83-91. doi: 10.3934/mbe.2009.6.83 |
[5] | Zhaohui Yuan, Xingfu Zou . Global threshold dynamics in an HIV virus model with nonlinear infection rate and distributed invasion and production delays. Mathematical Biosciences and Engineering, 2013, 10(2): 483-498. doi: 10.3934/mbe.2013.10.483 |
[6] | Gesham Magombedze, Winston Garira, Eddie Mwenje . Modelling the immunopathogenesis of HIV-1 infection and the effect of multidrug therapy: The role of fusion inhibitors in HAART. Mathematical Biosciences and Engineering, 2008, 5(3): 485-504. doi: 10.3934/mbe.2008.5.485 |
[7] | Thanh Nam Nguyen, Jean Clairambault, Thierry Jaffredo, Benoît Perthame, Delphine Salort . Adaptive dynamics of hematopoietic stem cells and their supporting stroma: a model and mathematical analysis. Mathematical Biosciences and Engineering, 2019, 16(5): 4818-4845. doi: 10.3934/mbe.2019243 |
[8] | Tinevimbo Shiri, Winston Garira, Senelani D. Musekwa . A two-strain HIV-1 mathematical model to assess the effects of chemotherapy on disease parameters. Mathematical Biosciences and Engineering, 2005, 2(4): 811-832. doi: 10.3934/mbe.2005.2.811 |
[9] | Avinash Shankaranarayanan, Hsiu-Chuan Wei . Mathematical modeling of SARS-nCoV-2 virus in Tamil Nadu, South India. Mathematical Biosciences and Engineering, 2022, 19(11): 11324-11344. doi: 10.3934/mbe.2022527 |
[10] | Ariel D. Weinberger, Alan S. Perelson . Persistence and emergence of X4 virus in HIV infection. Mathematical Biosciences and Engineering, 2011, 8(2): 605-626. doi: 10.3934/mbe.2011.8.605 |
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] | Rashid bin Mohd Yusoff A, Dai L, Cheng HM, et al. (2015) Graphene based energy devices. Nanoscale 7: 6881–6882. |
[2] | Ozyilmaz B, Saha S, Toh CT, et al. (2013) Photovoltaic cell with graphene-ferroelectric electrode. United States Patent Application US 14/385, 680. |
[3] |
Loh KP, Tong SW, Wu J (2016) Graphene and graphene-like molecules: prospects in solar cells. J Am Chem Soc 138: 1095–1102. doi: 10.1021/jacs.5b10917
![]() |
[4] |
Reinert F (2013) Graphene gets molecules into order. Nat Phys 9: 321–322. doi: 10.1038/nphys2643
![]() |
[5] |
Colson J W, Woll AR, Mukherjee A, et al. (2011) Oriented 2D covalent organic framework thin films on single-layer graphene. Science 332: 228–231. doi: 10.1126/science.1202747
![]() |
[6] |
Wierzbowska M, Wawrzyniak-Adamczewska M (2016) Cascade Donor-Acceptor Organic Ferroelectric Layers, between Graphene Sheets, for Solar Cell Applications. RSC Adv 6: 49988–49994. doi: 10.1039/C6RA07221F
![]() |
[7] |
Wigner E, Bardeen J (1935) Theory of the Work Functions of Monovalent Metals. Phys Rev 48: 84. doi: 10.1103/PhysRev.48.84
![]() |
[8] | Wawrzyniak-Adamczewska M,Wierzbowska M (2016) Separate-Path electron and hole transport across π-stacked ferroelectrics for photovoltaic applications. J Phys Chem C 120: 7748–7756. |
[9] | Ljungberg MP, Vänskä O, Koval P, et al. (2016) Charge transfer states at the interface of the pentacene monolayer on TiO2 and their influence on the optical spectrum. Available from: http://www.arXiv.org/abs/1610.01789. |
[10] | Masiak PM, Wierzbowska M (2017) Ferroelectric π-stacks of molecules with the energy gaps in the sunlight range. J Mater Sci [In Press]. |
[11] |
Xu L, Yu Y, Lin J, et al. (2016) On-surface synthesis of two-dimensional imine polymers with a tunable band gap: a combined STM, DFT and Monte Carlo investigation. Nanoscale 8: 8568. doi: 10.1039/C5NR07663C
![]() |
[12] |
Giannozzi P, Baroni S, Bonini N, et al. (2009) QUANTUM ESPRESSO: a modular and opensource software project for quantum simulation of materials. J Phys: Condens Matter 21: 395502. doi: 10.1088/0953-8984/21/39/395502
![]() |
[13] |
Perdew JP, Burke K, Ernzerhof M (1996) Generalized Gradient Approximation made simple. Phys Rev Lett 77: 3865–3868. doi: 10.1103/PhysRevLett.77.3865
![]() |
[14] |
Monkhorst HD, Pack JD (1976) Special points for Brillouin-zone integrations. Phys Rev B 13: 5188–5192. doi: 10.1103/PhysRevB.13.5188
![]() |
[15] | Fletcher R (1987) Practical methods of optimization, 2nd. Eds., John Wiley & sons, Chichester, UK. |
[16] |
Marini A, Hogan C, Grüning M, et al. (2009) YAMBO: an ab initio tool for excited state calculations. Comp Phys Comm 180: 1392–1403. doi: 10.1016/j.cpc.2009.02.003
![]() |
[17] |
Filippetti A, Spaldin NA (2003) Self-interaction-corrected pseudopotential scheme for magnetic and strongly-correlated systems. Phys Rev B 67: 125109. doi: 10.1103/PhysRevB.67.125109
![]() |
[18] |
Wierzbowska M, Majewski JA (2011) Forces and Atomic Relaxation in Density Functional Theory with the Pseudopotential Self-Interaction Correction. Phys Rev B 84: 245129. doi: 10.1103/PhysRevB.84.245129
![]() |
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 |
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 |