Next Article in Journal
Plasmon-Induced Transparency for Tunable Atom Trapping in a Chiral Metamaterial Structure
Next Article in Special Issue
Electrical and Optical Characterization of CsPbCl3 Films around the High-Temperature Phase Transitions
Previous Article in Journal
Pharmacokinetics of PEGylated Gold Nanoparticles: In Vitro—In Vivo Correlation
Previous Article in Special Issue
The Low-Temperature Photocurrent Spectrum of Monolayer MoSe2: Excitonic Features and Gate Voltage Dependence
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Stability and Bandgap Engineering of In1−xGaxSe Monolayer

by
Mattia Salomone
1,*,
Federico Raffone
2,*,
Michele Re Fiorentin
2,
Francesca Risplendi
1 and
Giancarlo Cicero
1
1
Dipartimento di Scienza Applicata e Tecnologia, Politecnico di Torino, 10129 Torino, Italy
2
Center for Sustainable Future Technologies, Istituto Italiano di Tecnologia, 10144 Torino, Italy
*
Authors to whom correspondence should be addressed.
Nanomaterials 2022, 12(3), 515; https://doi.org/10.3390/nano12030515
Submission received: 31 December 2021 / Revised: 26 January 2022 / Accepted: 28 January 2022 / Published: 1 February 2022
(This article belongs to the Special Issue Optoelectronic Properties and Applications of Nanomaterials)

Abstract

:
Bandgap engineering of semiconductor materials represents a crucial step for their employment in optoelectronics and photonics. It offers the opportunity to tailor their electronic and optical properties, increasing the degree of freedom in designing new devices and widening the range of their possible applications. Here, we report the bandgap engineering of a layered InSe monolayer, a superior electronic and optical material, by substituting In atoms with Ga atoms. We developed a theoretical understanding of I n 1 x G a x S e stability and electronic properties in its whole compositional range ( x = 0 1 ) through first-principles density functional theory calculations, the cluster expansion method, and kinetic Monte Carlo simulations. Our findings highlight the possibility of modulating the InGaSe bandgap by ≈0.41 eV and reveal that this compound is an excellent candidate to be employed in many optoelectronic and photonic devices.

1. Introduction

Photonic technologies represent a rapidly emerging field, which is considered one of the most promising directions for industrial and scientific innovation, thanks to the wide variety of possible applications, ranging from quantum cryptography and quantum computing to sensors for biomedical imaging and light-emitting diodes [1,2,3,4,5,6,7]. Compared to traditional bulk photonic materials such as GaAs and Si, two dimensional (2D) materials own naturally passivated surfaces without any dangling bond, allowing for their easy integration with photonic structures [8,9,10,11]. They offer the opportunity of creating vertical heterostructures without the “lattice mismatch” issue, because of weak interactions between different monolayers (van der Waals forces). In addition, most 2D structures interact strongly with light and can be used to cover a wide electromagnetic spectrum thanks to their diverse electronic properties, ranging from insulating hexagonal boron nitride (hBN) [8,12] and semiconducting transition metal dichalcogenides (TMDs) [13,14] to semi-metallic graphene [15].
In this regard, InSe is a layered III-IV chalcogenides material that in recent years has attracted interest due to its optical and electronic properties [16,17]. In its monolayer phase, it shares the promising characteristics listed so far for 2D materials, and, in addition, it shows an extremely high carrier mobility (≈10,000 cm 2 V 1 s 1 ) [18]. For photonic applications, it is required to develop 2D materials with continuous tunable bandgaps, and doping is one of the best techniques that can be exploited for this purpose [19]. Substitutional defects in a pristine monolayer may lead to a change in the lattice structure of the host material, resulting in a wide range bandgap tunability [20,21]. In our recent work [22], we defected 2D InSe with different substitutional impurities to generate well-defined and sharp defect states within the bandgap of the material, and we showed the possibility of using these new energy levels to achieve efficient emissions of single photons [23,24]. Here, considering the excellent optoelectronic properties and potential photonic applications shared by InSe and GaSe [25,26,27], and taking into account that they are almost lattice-matched and possess thickness-dependent tunable bandgaps [28,29] and very similar band structures (Figure 1b), we report a study of InGaSe compounds. In particular, we demonstrate the possibility of modulating the bandgap of an InSe monolayer by substituting In with Ga atoms. In the first part of the work, we study the thermodynamic stability and lowest energy structures of I n 1 x G a x S e compounds by means of the cluster expansion (CE) method coupled with density functional theory (DFT) simulations. We then verify I n 1 x G a x S e kinetic stability and present the more energetically favored lattice configurations for different indium to gallium fractions with a kinetic Monte Carlo (KMC) approach. In the second part of the work, we study the electronic properties of the more stable structures for different Ga concentrations in terms of band diagrams and density of states (DOS).

2. Methods

The CE is a method used to calculate the total energy of a compound existing in different configurations. As such, it is often exploited to predict stability and concentration-dependent phase diagrams of alloys/compounds. The lattice of a binary compound can be thought as composed of sites, which can be occupied either by atom A or atom B. This aspect is modeled in CE by means of an occupational variable σ I which can assume the values ± 1 , depending on a the atomic species occupying the I-th site. The total energy of a certain compound configuration is then given by the summation of the energy contribution associated with clusters (or figures) of these occupational variables (singlet, pairs, triplets... up to the M-body cluster, where M is the total number of sites in the lattice). The energy contribution of each cluster is referred to as the specific energetic contribution (ECI), indicated by J. The general Hamiltonian of a compound existing in several possible configurations can be written as [30]
H ( σ ) = f m f J f I f σ I .
Here, m f is the figure multiplicity which takes into account degenerate clusters, and σ is the configuration vector, which describes all the possible arrangements of atom A and B. The summation is over all clusters, which can involve one site (for singlets), two sites (for pairs), and so forth up to M. Because of the extremely high number of possible clusters, the summation is truncated. It was demonstrated that even with a finite number of terms the configurational energy can be represented with good precision [31,32].
DFT simulations were used to extract the ECIs necessary to find the Hamiltonian of the system by calculating the energy of a small set of structures with known lattice configurations [33]. Equation (1) is then reversed and filled with the data from DFT so that the ECIs are derived. The CE Hamiltonian, obtained by the calculated ECIs, is then used to quickly predict the energy of all configurations in a concentration range, allowing also for the study of structures that are too computationally demanding for DFT. The degree of accuracy of the expansion is monitored by a cross validation score, χ V L , which evaluates the difference between predicted and known energies (calculated by DFT) for different sets of data. We assume that when its value is 0.008 eV, which is well below the one used in the literature [34], the expansion guarantees reliable results. CE was performed by means of the Alloy Theoretic Automated Toolkit (ATAT) software package [33] and DFT calculations with the Quantum Espresso package [35,36].
For DFT simulations, we exploited ultrasoft pseudopotentials [37] to describe the electron-ion interaction and the gradient corrected Perdew Burke Ernzerhof (PBE) functional [38] to describe the exchange-correlation effects. The Kohn-Sham wavefunctions were expanded in plane waves with a kinetic energy cutoff of 80 Ry, to ensure converged stress values during lattice parameters optimization. The Brillouin Zone of the primitive cell was sampled by employing a 6 × 6 × 1 Monkhorst-Pack grid [39] for monolayer unit cells and reduced accordingly as the supercell dimensions increased to ensure the same k-point density. All calculations were performed without considering the effect of spin-orbit coupling (SOC), since its inclusion does not strongly affect InSe and GaSe band structures [22,40,41].
The KMC approach was used to study the kinetic stability and the more favored lattice configurations of InGaSe for different gallium concentrations. The process was carried out by simulating a sequence of diffusive events, each of which was selected on the basis of its transition rate [31,42,43]. The latter were affected by the energetics of the system calculated by means of the ECIs obtained previously with the CE. The simulations were performed exploiting the Zacros software [44,45]. We used this code to model the diffusion of gallium defects through the compound structure via positional exchanges between Ga and In. The lattice was composed of two stacked layers (Se atoms were not active participants in the simulated processes, because Ga substitutes only In; therefore, they were not considered).

3. Results and Discussion

Figure 2 shows the phase diagram obtained through CE. On the x axis appears the Ga concentration in I n 1 x G a x S e (x = 0 corresponds to pure InSe, while x = 1 to pure GaSe), on the y axis, instead, are reported the formation energies calculated as
E F O R M = E I n 1 x G a x S e [ ( 1 x ) μ I n + x μ G a + μ S e ]
where x is the Ga concentration, E I n 1 x G a x S e is the total energy of the corresponding compound structure and μ I n , μ S e , and μ G a are the chemical potentials of In, Se, and Ga in their stable bulk phases (I4/mmm tetragonal indium, P 3 1 21 trigonal selenium, and Cmce gallium). Blue dots indicate energies of structures calculated with DFT, while gray crosses represent the energies of the structures that were predicted through the cluster expansion. Black dots indicate, instead, the energies of the ground states associated with the two pure phases, and the black dashed line that connects them marks the ground-state energy for any In/Ga fraction.
Analyzing the phase diagram, we are able to identify, for each concentration of Ga, the more stable structure, i.e., the one with a formation energy closer to the predicted ground state for a specific compound composition.
We found that a cluster expansion of 18 figures was sufficient to achieve a cross-validation score χ V L of ≈0.008 eV, which satisfied our precision requirements. In the expansion, we considered the energetic contribution of the background (i.e., the one associated with the pristine InSe monolayer), the singlet interaction (related to the interaction of a defect with itself), and pair interactions (nearest neighbors’ figures, shortened to NN for simplicity) up to a diameter of 16.3 Å, labeled from NN2 to NN16. All the interactions discussed in the following involve atoms belonging to the same atomic “plane” (see Figure 1a as a reference), while interactions between interplane pairs were associated with smaller ECIs; therefore, their contribution to the total energy of the system was minor.
The compound total energy was increased (lowered) by positive (negative) ECIs, if we are considering interactions between atoms of the same species. Therefore, it is more probable to observe atoms of the same type forming patterns reproducing the figures associated with negative ECIs. In Figure 3, we report the least (Figure 3a) and the most favored (Figure 3b) patterns associated with pairs of the same species. It is apparent that figures NN2, NN4, and NN6, which represent close neighbor couples, were observed with small probabilities, while figures NN12, NN14, and NN16 were more likely to occur. In the supporting information we report all the figures considered in the expansion (Figure S1) along with their respective ECIs (Table S1).
Some of the predicted energies (Figure 2, gray crosses) were within the cross-validation score sensitivity (their energy difference with respect to the predicted ground state for that specific Ga concentration was smaller than χ V L ). Considering also that the CE was performed using as reference the two pure phases of InSe and GaSe, the interface energy between the two material was not considered explicitly; instead, we studied the kinetic stability of the I n 1 x G a x S e structure by KMC simulations. We verified the possibility of creating mixed phases by studying the kinetic stability of the I n 1 x G a x S e compounds. The obtained results provided information about the equilibrium In/Ga arrangement within the lattice. Since we are interested in understanding whether the two pure phases tend to separate when mixed together, we used lattices containing already formed clusters as the initial structures. Figure 4a,c reports the results for two significant cases: a GaSe quantum dot ( x = 0.14 , supercell area of ≈26 nm 2 ) and a GaSe stripe ( x = 0.50 , supercell area of ≈65 nm 2 ).
Observing Figure 4b,d, respectively, it can be noticed that InSe and GaSe (well separated at the beginning of the simulations) tended to mix, forming a homogeneous compound. The initial structures presented an interface energy, which had a dominant effect, and led to cluster dissolution. This findings were supported by the analysis of ECIs obtained from the cluster expansion. Indeed, the formation of clusters was hindered because it would increase the total energy of the compound, according to the ECIs shown in Figure 3a, for the whole gallium concentration range. Nevertheless, one can distinguish a recurring pattern (Figure S3), which was more often observed for low Ga concentration. This arrangement almost reflects the periodic repetition of the most stable figures obtained by the CE, reported in Figure 3b.
Our findings highlight that when indium atoms were substituted with gallium atoms, the two pure phases tended to mix, forming homogeneous stable lattice configurations. This result shows that, for this specific system, the interface energy between the two pure phases (not considered by the CE) was relevant and cannot be neglected to uniquely identify the equilibrium structure, hence, a KMC approach is necessary.
Starting from the results provided by the CE, we chose the lowest energy InGaSe structures in the whole concentration range and performed an in depth analysis of the electronic properties of the material. DFT calculations showed that the bandgap changed monotonically with the gallium concentration: the higher the percentage of Ga, the wider the bandgap (Figure 5a). In particular, E G d (direct, calculated in Γ ) ranged from 1.47 eV (pure InSe) to 1.88 eV (pristine GaSe), while E G i (indirect, calculated between the valence band maximum, VBM, and the conduction band minimum, CBm) varied from 1.39 eV to 1.78 eV (see the supporting information for more details). The CB and VB of InSe at the Γ point showed a strong contribution of In s-orbitals and Se p z -orbitals, respectively, (in agreement with other works [46,47,48]). However, the topmost levels of the valence band were also influenced by p orbitals of indium. These results were confirmed for high concentrations of In, such as, the I n 0.83 G a 0.17 S e compound, which was the structure with the lowest Ga percentage between those studied to monitor the bandgap behavior. In Figure 6a–c, we plotted the k-resolved projected DOS, to highlight the contribution of each orbital to the band diagram. Very similar bandgap behaviors were obtained in other studies for other compounds containing In and Ga, in particular for I n 1 x G a x P and I n 1 x G a x A s [49].
Finally, we aligned the I n 1 x G a x S e band diagrams (using the vacuum level as common reference) to assess the possibility of creating heterostructures with useful band alignment and studied the optical properties of the compound (this analysis is reported in the supporting information). Heterostructures are crucial components in many electronic, photonic, and optoelectronic devices such as field effect transistors, solar cells, and sensors; hence, their possible implementation using InGaSe would certainly widen the range of possible applications for this compounds. In Figure 5b, in particular, we reported the comparison considering direct bandgaps. We observed a sequence of type I heterostructures between InSe and the other compounds. The valence band maxima of I n 0.83 G a 0.17 S e , I n 0.50 G a 0.50 S e and I n 0.17 G a 0.83 S e were very close in energy; therefore, heterojunctions involving only these materials were not considered. Knowing that k B T at room temperature was ≈0.025 eV, it was apparent that most of the highlighted bandgap discontinuities were greater than this value, ranging from ≈2 k B T ( Δ E V InSe/ I n 0.17 G a 0.83 S e ) to 22 k B T ( Δ E C InSe/GaSe). This fact implies that, in principle, it is possible to make heterojunctions featuring efficient charge separation. As an example, one may build a multijunction solar cell by stacking InGaSe layers at increasing Ga concentrations and obtain a structure able to efficiently harvest a good portion of the solar spectrum.

4. Conclusions

In conclusion, we predicted I n 1 x G a x S e structure and thermodynamic stability exploiting the CE method. In addition, we overcame the limitation of neglecting the interface energy between the two pure phases with a KMC approach, and we showed the possibility of creating kinetically stable compounds, which did not exhibit phase separation and were characterized by concentration dependent bandgaps. The tuning was very effective and, considering the two pure phases, the variations of the energy gap values were 0.41 eV (direct) and 0.39 eV (indirect). In addition, heterostructures made of InSe/GaSe and a generic I n 1 x G a x S e compound exhibited high energy band discontinuities and can be used to actuate either charge separation or carrier confinement. These findings point out that InSe is an excellent candidate to be used in optoelectronics and photonics, offering many degrees of freedom in designing new devices for a wide range of possible applications.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/nano12030515/s1, Figure S1: pair interactions obtained from CE. Figure S2: mapping of the 3D crystal structure to the CE model. Figure S3: figures associated to the most negative ECIs. Figure S4: DFT-obtained bandgaps and their interpolations with a linear function and a quadratic form with bowing parameter C. Figure S5: optical absorption spectra of pristine InSe, GaSe and three representative InGaSe compounds, computed at the IP-RPA level, along the zig-zag and armchair directions.; Table S1: Effective Cluster Interactions of the figures obtained from the CE method.

Author Contributions

Conceptualization, G.C.; Calculations, M.S.; Data analysis and discussion, M.S., F.R. (Federico Raffone), M.R.F., F.R. (Francesca Risplendi) and G.C. All authors have contributed to writing and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The authors acknowledge CINECA for the availability of high performance computing resources under the Iscra-B initiative, as well as the computational support provided by HPC@POLITO (http://www.hpc.polito.it).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Hu, H.; Tang, B.; Wan, H.; Sun, H.; Zhou, S.; Dai, J.; Chen, C.; Liu, S.; Guo, L.J. Boosted ultraviolet electroluminescence of InGaN/AlGaN quantum structures grown on high-index contrast patterned sapphire with silica array. Nano Energy 2020, 69, 104427. [Google Scholar] [CrossRef]
  2. Zhao, X.; Tang, B.; Gong, L.; Bai, J.; Ping, J.; Zhou, S. Rational construction of staggered InGaN quantum wells for efficient yellow light-emitting diodes. Appl. Phys. Lett. 2021, 118, 182102. [Google Scholar] [CrossRef]
  3. Zhou, S.; Liu, X.; Yan, H.; Chen, Z.; Liu, Y.; Liu, S. Highly efficient GaN-based high-power flipchip light-emitting diodes. Opt. Express 2019, 27, A669–A692. [Google Scholar] [CrossRef] [PubMed]
  4. Chitraleema, C.; Nick, V.; Dirk, E. Advances in quantum light emission from 2D materials. Nanophotonics 2019, 8, 2017–2032. [Google Scholar] [CrossRef]
  5. Lodahl, P.; Mahmoodian, S.; Stobbe, S. Interfacing single photons and single quantum dots with photonic nanostructures. Rev. Mod. Phys. 2015, 87, 347–400. [Google Scholar] [CrossRef] [Green Version]
  6. Aspuru-Guzik, A.; Walther, P. Photonic quantum simulators. Nat. Phys. 2012, 8, 285–291. [Google Scholar] [CrossRef] [Green Version]
  7. Lo, H.; Curty, M.; Tamaki, K. Secure quantum key distribution. Nat. Photonics 2014, 8, 595–604. [Google Scholar] [CrossRef] [Green Version]
  8. Xia, F.; Wang, H.; Xiao, D.; Dubey, M.; Ramasubramaniam, A. Two-dimensional material nanophotonics. Nat. Photonics 2014, 8, 899–907. [Google Scholar] [CrossRef]
  9. Gan, X.; Shiue, R.J.; Gao, Y.; Meric, I.; Heinz, T.F.; Shepard, K.; Hone, J.; Assefa, S.; Englund, D. Chip-integrated ultrafast graphene photodetector with high responsivity. Nat. Photonics 2013, 7, 883–887. [Google Scholar] [CrossRef]
  10. Furchi, M.; Urich, A.; Pospischil, A.; Lilley, G.; Unterrainer, K.; Detz, H.; Klang, P.; Andrews, A.M.; Schrenk, W.; Strasser, G.; et al. Microcavity-integrated graphene photodetector. Nano Lett. 2012, 12, 2773–2777. [Google Scholar] [CrossRef]
  11. Gan, X.; Mak, K.F.; Gao, Y.; You, Y.; Hatami, F.; Hone, J.; Heinz, T.F.; Englund, D. Strong enhancement of light-matter interaction in graphene coupled to a photonic crystal nanocavity. Nano Lett. 2012, 12, 5626–5631. [Google Scholar] [CrossRef]
  12. Dean, C.R.; Young, A.F.; Meric, I.; Lee, C.; Wang, L.; Sorgenfrei, S.; Watanabe, K.; Taniguchi, T.; Kim, P.; Shepard, K.L.; et al. Boron nitride substrates for high-quality graphene electronics. Nat. Nanotechnol. 2010, 5, 722–726. [Google Scholar] [CrossRef]
  13. Mak, K.F.; Lee, C.; Hone, J.; Shan, J.; Heinz, T.F. Atomically thin MoS2: A New Direct-Gap Semiconductor. Phys. Rev. Lett. 2010, 105, 136805. [Google Scholar] [CrossRef] [Green Version]
  14. Eda, G.; Maier, S.A. Two-dimensional crystals: Managing light for optoelectronics. ACS Nano 2013, 7, 5660–5665. [Google Scholar] [CrossRef] [PubMed]
  15. Gul, H.Z.; Sakong, W.; Ji, H.; Torres, J.; Yi, H.; Ghimire, M.K.; Yoon, J.H.; Yun, M.H.; Hwang, H.R.; Lee, Y.H.; et al. Semimetallic Graphene for Infrared Sensing. ACS Appl. Mater. Interfaces 2019, 11, 19565–19571. [Google Scholar] [CrossRef]
  16. Yang, Z.; Jie, W.; Mak, C.H.; Lin, S.; Lin, H.; Yang, X.; Yan, F.; Lau, S.P.; Hao, J. Wafer-Scale Synthesis of High-Quality Semiconducting Two-Dimensional Layered InSe with Broadband Photoresponse. ACS Nano 2017, 11, 4225–4236. [Google Scholar] [CrossRef]
  17. Lei, S.; Ge, L.; Najmaei, S.; George, A.; Kappera, R.; Lou, J.; Chhowalla, M.; Yamaguchi, H.; Gupta, G.; Vajtai, R.; et al. Evolution of the Electronic Band Structure and Efficient Photo-Detection in Atomic Layers of InSe. ACS Nano 2014, 8, 1263–1272. [Google Scholar] [CrossRef] [PubMed]
  18. Bandurin, D.A.; Tyurnina, A.V.; Geliang, L.Y.; Mishchenko, A.; Zólyomi, V.; Morozov, S.V.; Kumar, R.K.; Gorbachev, R.V.; Kudrynskyi, Z.R.; Pezzini, S.; et al. High electron mobility, quantum Hall effect and anomalous optical response in atomically thin InSe. Nat. Nanotechnol. 2016, 12, 223–227. [Google Scholar] [CrossRef]
  19. Susarla, S.; Kutana, A.; Hachtel, J.A.; Kochat, V.; Apte, A.; Vajtai, R.; Idrobo, J.C.; Yakobson, B.I.; Tiwary, C.S.; Ajayan, P.M. Quaternary 2D transition metal dichalcogenides (TMDS) with tunable bandgap. Adv. Mater. 2017, 29, 1702457. [Google Scholar] [CrossRef] [PubMed]
  20. Raffone, F.; Ataca, C.; Grossman, J.C.; Cicero, G. MoS2 Enhanced T-Phase Stabilization and Tunability through Alloying. J. Phys. Chem. Lett. 2016, 7, 2304–2309. [Google Scholar] [CrossRef]
  21. Verma, A.K.; Raffone, F.; Cicero, G. Prediction of the structural and electronic properties of MoxTi1-xS2 monolayers via first principle simulations. Nanomater. Nanotechnol. 2020, 10, 1–6. [Google Scholar] [CrossRef]
  22. Salomone, M.; Fiorentin, M.R.; Cicero, G.; Risplendi, F. Point Defects in Two-Dimensional Indium Selenide as Tunable Single-Photon Sources. J. Phys. Chem. Lett. 2021, 12, 10947–10952. [Google Scholar] [CrossRef]
  23. Eisaman, M.D.; Fan, J.; Migdall, A.; Polyakov, S.V. Invited Review Article: Single-photon sources and detectors. Rev. Sci. Instrum. 2011, 82, 071101. [Google Scholar] [CrossRef] [Green Version]
  24. Darquie, B. Controlled Single-Photon Emission from a Single Trapped Two-Level Atom. Science 2005, 309, 454–456. [Google Scholar] [CrossRef] [Green Version]
  25. Tonndorf, P.; Schwarz, S.; Kern, J.; Niehues, I.; Pozo-Zamudio, O.D.; Dmitriev, A.I.; Bakhtinov, A.P.; Borisenko, D.N.; Kolesnikov, N.N.; Tartakovskii, A.I.; et al. Single-photon emitters in GaSe. 2D Mater. 2017, 4, 021010. [Google Scholar] [CrossRef]
  26. Zhou, J.; Shi, J.; Zeng, Q.; Chen, Y.; Niu, L.; Liu, F.; Yu, T.; Suenaga, K.; Liu, X.; Lin, J.; et al. InSe monolayer: Synthesis, structure and ultra-high second-harmonic generation. 2D Mater. 2018, 5, 025019. [Google Scholar] [CrossRef]
  27. Zhou, X.; Cheng, J.; Zhou, Y.; Cao, T.; Hong, H.; Liao, Z.; Wu, S.; Peng, H.; Liu, K.; Yu, D. Strong Second-Harmonic Generation in Atomic Layered GaSe. J. Am. Chem. Soc. 2015, 137, 7994–7997. [Google Scholar] [CrossRef]
  28. Balakrishnan, N.; Kudrynskyi, Z.R.; Smith, E.F.; Fay, M.W.; Makarovsky, O.; Kovalyuk, Z.D.; Eaves, L.; Beton, P.H.; Patanè, A. Engineering p – n junctions and bandgap tuning of InSe nanolayers by controlled oxidation. 2D Mater. 2017, 4, 025043. [Google Scholar] [CrossRef]
  29. Li, Y.; Wang, T.; Wu, M.; Cao, T.; Chen, Y.; Sankar, R.; Ulaganathan, R.K.; Chou, F.; Wetzel, C.; Xu, C.Y.; et al. Ultrasensitive tunability of the direct bandgap of 2D InSe flakes via strain engineering. 2D Mater. 2018, 5, 021002. [Google Scholar] [CrossRef] [Green Version]
  30. Wu, Q.; He, B.; Song, T.; Gao, J.; Shi, S. Cluster expansion method and its application in computational materials science. Comput. Mater. Sci. 2016, 125, 243–254. [Google Scholar] [CrossRef]
  31. Raffone, F.; Savazzi, F.; Cicero, G. Controlled Pore Generation in Single-Layer Graphene Oxide for Membrane Desalination. J. Phys. Chem. Lett. 2019, 10, 7492–7497. [Google Scholar] [CrossRef]
  32. Connolly, J.W.D.; Williams, A.R. Density-functional theory applied to phase transformations in transition-metal alloys. Phys. Rev. B Condens. Matter Mater. Phys. 1983, 27, 5169. [Google Scholar] [CrossRef]
  33. van de Walle, A.; Ceder, G. Automating first-principles phase diagram calculations. J. Phase Equilib. 2002, 23, 348. [Google Scholar] [CrossRef] [Green Version]
  34. van de Walle, A.; Asta, M.; Ceder, G. The alloy theoretic automated toolkit: A user guide. Calphad 2002, 26, 539–553. [Google Scholar] [CrossRef] [Green Version]
  35. Giannozzi, P.; Baroni, S.; Bonini, N.; Calandra, M.; Car, R.; Cavazzoni, C.; Ceresoli, D.; Chiarotti, G.L.; Cococcioni, M.; Dabo, I.; et al. QUANTUM ESPRESSO: A modular and open-source software project for quantum simulations of materials. J. Phys.-Condens. Matter 2009, 21, 395502. [Google Scholar] [CrossRef]
  36. Giannozzi, P.; Andreussi, O.; Brumme, T.; Bunau, O.; Nardelli, M.B.; Calandra, M.; Car, R.; Cavazzoni, C.; Ceresoli, D.; Cococcioni, M.; et al. Advanced capabilities for materials modelling with quantum espresso. J. Phys.-Condens. Matter 2017, 29, 465901. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Schlipf, M.; Gygi, F. Optimization algorithm for the generation of ONCV pseudopotentials. Comput. Phys. Commun. 2015, 196, 36–44. [Google Scholar] [CrossRef] [Green Version]
  38. Perdew, J.; Burke, K.; Wang, Y. Generalized gradient approximation for the echange-correlation hole of many electron system. Phys. Rev. B 1996, 54, 16533–16539. [Google Scholar] [CrossRef] [Green Version]
  39. Monkhorst, H.J.; Pack, J.D. Special points for brillouin-zone integrations. Phys. Rev. B 1976, 13, 5188. [Google Scholar] [CrossRef]
  40. Magorrian, S.J.; Zolyomi, V.; Fal’ko, V.I. Electronic and optical properties of two-dimensional InSe from a DFT-parametrized tight-binding model. Phys. Rev. B 2016, 94, 245431. [Google Scholar] [CrossRef] [Green Version]
  41. Ghalouci, L.; Benbahi, B.; Hiadsi, S.; Abidri, B.; Vergoten, G.; Ghalouci, F. First principle investigation into hexagonal and cubic structures of Gallium Selenide. Comput. Mater. Sci. 2013, 67, 73–82. [Google Scholar] [CrossRef]
  42. Jansen, A.P.J. An Introduction to Kinetic Monte Carlo Simulations of Surface Reactions; Springer: Berlin, Germany, 2012; p. 856. [Google Scholar]
  43. Raffone, F.; Cicero, G. Unveiling the fundamental role of temperature in RRAM switching mechanism by multiscale simulations. ACS Appl. Mater. Interfaces 2018, 10, 7512–7519. [Google Scholar] [CrossRef] [PubMed]
  44. Nielsen, J.; D’Avezac, M.; Hetherington, J.; Stamatakis, M. Parallel kinetic Monte Carlo simulation framework incorporating accurate models of adsorbate lateral interactions. J. Chem. Phys. 2013, 224706, 139. [Google Scholar] [CrossRef] [Green Version]
  45. Stamatakis, M.; Vlachos, D.G. A graph-theoretical kinetic Monte Carlo framework for on-lattice chemical kinetics. J. Chem. Phys. 2011, 134, 214115. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Li, H.; Han, X.; Pan, D.; Yan, X.; Wang, H.W.; Wu, C.; Cheng, G.; Zhang, H.; Yang, S.; Li, B.; et al. Bandgap engineering of InSe Single Crystals through S substitution. Cryst. Growth Des. 2018, 18, 2899–2904. [Google Scholar] [CrossRef]
  47. Ma, Y.; Dai, Y.; Yu, L.; Liu, C.; Huang, B. Engineering a topological phase transition in β-InSe via strain. New J. Phys. 2013, 15, 073008. [Google Scholar] [CrossRef]
  48. Zolyomi, V.; Drummond, N.D.; Falko, V.I. Electrons and phonons in single layers of hexagonal indium chalcogenides from ab initio calculations. Phys. Rev. B Condens. Matter Mater. Phys. 2014, 89, 205416. [Google Scholar] [CrossRef] [Green Version]
  49. Nahory, R.E.; Pollack, M.A.; Johnston, W.D., Jr.; Barns, R.L. Band gap versus composition and demonstration of Vegard’s law for In1−xGaxAsyP1−y lattice matched to InP. Appl. Phys. Lett. 2008, 33, 659–661. [Google Scholar] [CrossRef]
Figure 1. Top and side views (a) of the geometry of a pristine InSe (GaSe) monolayer. Selenium is represented in blue, while indium (gallium) is in orange. The monolayer is formed by four stacked atomic planes; therefore, two different levels in the In/Ga atoms arrangement can be distinguished (highlighted with the two horizontal dashed lines). This lattice arrangement is also shared by InGaSe compounds. (b) InSe band structure (black lines) and GaSe band structure (dashed red lines) calculated at the DFT-PBE level. In both cases, the vacuum level ( E 0 ) is set to 0 eV.
Figure 1. Top and side views (a) of the geometry of a pristine InSe (GaSe) monolayer. Selenium is represented in blue, while indium (gallium) is in orange. The monolayer is formed by four stacked atomic planes; therefore, two different levels in the In/Ga atoms arrangement can be distinguished (highlighted with the two horizontal dashed lines). This lattice arrangement is also shared by InGaSe compounds. (b) InSe band structure (black lines) and GaSe band structure (dashed red lines) calculated at the DFT-PBE level. In both cases, the vacuum level ( E 0 ) is set to 0 eV.
Nanomaterials 12 00515 g001
Figure 2. Phase diagram of the I n 1 x G a x S e : gray crosses represent the energies of the structures predicted through the CE, blue dots indicate the energies found with DFT calculations performed on generated structures, black dots represent, instead, the energies of the two ground states associated with the two pure phases of InSe and GaSe, and the black dashed line marks the ground-state energy for all In/Ga fraction.
Figure 2. Phase diagram of the I n 1 x G a x S e : gray crosses represent the energies of the structures predicted through the CE, blue dots indicate the energies found with DFT calculations performed on generated structures, black dots represent, instead, the energies of the two ground states associated with the two pure phases of InSe and GaSe, and the black dashed line marks the ground-state energy for all In/Ga fraction.
Nanomaterials 12 00515 g002
Figure 3. Representation of the less (a) and the more (b) stable pair figures obtained with the CE. ECIs related to red interactions (b) are negative; therefore, these configurations lower the total energy of the structure. Blue figures (a), instead, are associated with positive J coefficients, which implies that their formation is not energetically favored. Indium atoms are represented in black, selenium atoms in gray, and gallium atoms in blue (a) and red (b).
Figure 3. Representation of the less (a) and the more (b) stable pair figures obtained with the CE. ECIs related to red interactions (b) are negative; therefore, these configurations lower the total energy of the structure. Blue figures (a), instead, are associated with positive J coefficients, which implies that their formation is not energetically favored. Indium atoms are represented in black, selenium atoms in gray, and gallium atoms in blue (a) and red (b).
Nanomaterials 12 00515 g003
Figure 4. Panels (a,c) show the two initial KMC structures considered: a quantum dot and a stripe, respectively. Atoms belonging to the "bottom" plane are transparent, while their color represents the atomic species, gallium (red) and indium (light gray). Panels (b,d) report the equilibrium structures obtained from the two corresponding initial configurations shown in (a,c). It is apparent that in both cases the two pure phases tend to mix.
Figure 4. Panels (a,c) show the two initial KMC structures considered: a quantum dot and a stripe, respectively. Atoms belonging to the "bottom" plane are transparent, while their color represents the atomic species, gallium (red) and indium (light gray). Panels (b,d) report the equilibrium structures obtained from the two corresponding initial configurations shown in (a,c). It is apparent that in both cases the two pure phases tend to mix.
Nanomaterials 12 00515 g004
Figure 5. (a) I n 1 x G a x S e bandgap (direct blue, indirect gray) at varying Ga concentrations. (b) E G d ( x ) alignment with respect to the vacuum level for some of the InGaSe structures.
Figure 5. (a) I n 1 x G a x S e bandgap (direct blue, indirect gray) at varying Ga concentrations. (b) E G d ( x ) alignment with respect to the vacuum level for some of the InGaSe structures.
Nanomaterials 12 00515 g005
Figure 6. Band diagrams and k-resolved projected density of states, for I n 0.83 G a 0.17 S e (ac) and I n 0.17 G a 0.83 S e (d). Selenium p orbitals provide a very similar contribution for both the considered compounds. S and p orbitals give a strong contribution both to conduction and valence bands, when the Ga percentage is low (panels b,c). When the Ga concentration is high, Ga s (not shown here) and p orbitals give a larger contribution.
Figure 6. Band diagrams and k-resolved projected density of states, for I n 0.83 G a 0.17 S e (ac) and I n 0.17 G a 0.83 S e (d). Selenium p orbitals provide a very similar contribution for both the considered compounds. S and p orbitals give a strong contribution both to conduction and valence bands, when the Ga percentage is low (panels b,c). When the Ga concentration is high, Ga s (not shown here) and p orbitals give a larger contribution.
Nanomaterials 12 00515 g006
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Salomone, M.; Raffone, F.; Re Fiorentin, M.; Risplendi, F.; Cicero, G. Stability and Bandgap Engineering of In1−xGaxSe Monolayer. Nanomaterials 2022, 12, 515. https://doi.org/10.3390/nano12030515

AMA Style

Salomone M, Raffone F, Re Fiorentin M, Risplendi F, Cicero G. Stability and Bandgap Engineering of In1−xGaxSe Monolayer. Nanomaterials. 2022; 12(3):515. https://doi.org/10.3390/nano12030515

Chicago/Turabian Style

Salomone, Mattia, Federico Raffone, Michele Re Fiorentin, Francesca Risplendi, and Giancarlo Cicero. 2022. "Stability and Bandgap Engineering of In1−xGaxSe Monolayer" Nanomaterials 12, no. 3: 515. https://doi.org/10.3390/nano12030515

APA Style

Salomone, M., Raffone, F., Re Fiorentin, M., Risplendi, F., & Cicero, G. (2022). Stability and Bandgap Engineering of In1−xGaxSe Monolayer. Nanomaterials, 12(3), 515. https://doi.org/10.3390/nano12030515

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop