Wurtzite materials in alloys of rock salt compounds

Abstract

Materials with crystal structures containing tetrahedral motifs are preferable for optoelectronic applications because they often have direct band gaps and low electron effective masses. However, crystal structures of manganese chalcogenides typically contain octahedral motifs, such as in rock salt (RS) MnS and MnSe materials. Here, we experimentally show that MnS_1−xSe_x alloys with tetrahedrally bonded wurtzite (WZ) structure can form between MnSe and MnS parent compounds with octahedral RS structures, at S-rich compositions (x < 0.4) and low synthesis temperatures (∼300 °C). The calculated mixing enthalpies of MnS_1−xSe_x alloys in RS and WZ structures cannot explain this experimental observation, so we hypothesize that WZ stabilization may be related to smaller structure density and lower surface energy compared with RS. The resulting WZ MnS_1−xSe_x alloys have 3.0–3.2 eV optical absorption onset and lower electrical conductivity (<0.0001 S/cm) than the parent RS compounds. These experimental measurement results are consistent with computationally predicted band gaps and effective masses.

Introduction

Binary metal-chalcogenide semiconductors are w idely used in optoelectronic devices [1]. Structures containing tetrahedrally coordinated (T _d) bonding motifs are often preferable for these applications because they have wide direct band gaps and low electron effective masses, and because they are structurally compatible with other chalcogenide materials. For example, CdTe and CdSe used as absorber layers in thin-film solar cells and ZnS, ZnSe, and CdS used as transport layers in electronics all assume either zinc blende (ZB) or wurtzite (WZ) structures [2, 3]. Manganese chalcogenides are also compelling materials because they are chemically compatible with other binary chalcogenides and have excellent semiconducting properties. For example, MnO, MnS, and MnSe are wide band gap materials, MnTe exhibits high p-type conductivity, and MnCh (Ch = S, Se, Te) are important for instilling magnetism in semiconductors [4]. However, MnCh ground states contain octahedrally coordinated (O _h) structural motifs, which are not generally compatible with other commonly used T _d structures such as WZ and ZB: MnO, MnS, and MnSe crystallize in the rock salt (RS) structure, and MnTe crystalizes in the nickeline (NC) structure. Besides, the MnCh with O _h structures show indirect band gaps [5], further hindering their optoelectronic applications. Therefore, stabilizing T _d polymorphs of MnCh may enable their future deployment as active optical materials in optoelectronic devices.

There have been several prior reports about changing O _h structures (like RS and NC) to T _d structures (like ZB and WZ) in MnCh by heterostructural alloying with T _d ZnCh leading to novel optoelectronic properties. In the oxides, the insulating RS phase of MnO was changed to the semiconducting WZ phase in the Mn_1−xZn_xO alloy at Zn compositions above x = 0.3 and applied as water splitting photoelectrodes [6]. Similarly, in the selenides, the RS structure of MnSe was converted to the WZ structure in the Zn_1−xMn_xSe alloys [7]. In the tellurides, by alloying NC MnTe with trace amount of ZB ZnTe, the ground state switched from the narrow band gap NC structure to the wide band gap WZ structure [8, 9]. Despite the success of alloying for stabilization of WZ MnO, MnSe, and MnTe, there is only one report on using alloying to stabilize WZ MnS [4]. However, such stabilization of WZ MnS is likely energetically feasible because WZ MnS occurs naturally [10], and has been synthesized via chemical bath deposition [11], heteroepitaxial growth on GaAs [12], and other methods [13].

Studies demonstrating that structures containing T _d motifs can form from alloys containing MnCh end-member compounds with O _h motifs have been reported. For example, the WZ structure was observed at intermediate compositions in related MnSe_1−xTe_x alloys formed by alloying RS MnSe and NC MnTe [14]. The proposed stabilization mechanism in that case was that the less dense metastable phases (T _d polymorphs) can form in alloys of more dense end-member compounds (O _h structures) due to the lower competition for space in the less dense phases [14]. Depending on the mixing enthalpies and end-member formation energies, this stabilization should also be feasible in iso-structural alloy systems [15]. A computational study was also reported which predicts that T _d phases can form by alloying end-member compounds with the same O _h structures [16]; however, this theory has not been demonstrated experimentally yet. We also note that alloy stabilization of T _d phases by anion substitution (e.g., Se–Te) compared with alloy stabilization by cation substitution (e.g., Mn–Zn) may be more beneficial in optoelectronic devices. This is because in binary compound semiconductors, the conduction band minimum is primarily composed of cation rather than anion states [5], so anion alloying may lead to less perturbation of the low electron effective mass in T _d structures.

Here, we report the stabilization of a structure with T _d bonding motifs in alloys containing two end-member compounds with O _h motifs. Specifically, we find that the WZ structure forms at intermediate compositions in MnS_1−xSe_x alloys, despite end-member MnS and MnSe compounds exhibiting RS structures. Structural and optoelectronic properties of MnS_1−xSe_x (0 ≤ x ≤ 1) alloy thin films prepared by combinatorial sputtering are reported as a function of chemical compositions and synthesis temperatures. Although films over the entire MnS_1−xSe_x composition range crystalize in the RS structure at higher synthesis temperature (420 °C), at lower temperatures (280–340 °C) the WZ phase appears at S-rich compositions (x < 0.4). The optical absorption onset of WZ MnS_1−xSe_x alloys is 3.0–3.2 eV, whereas the conductivity is <0.0001 S/cm, indicating its potential for use as a wide band gap semiconductor in various optoelectronic applications. Theoretical mixing enthalpies, electron/hole effective masses, and band gaps are also reported to help better understand the properties of MnS_1−xSe_x alloys, and possible mechanisms of the WZ phase stabilization are discussed.

Results and discussion

Structural properties

Figure 1(a) shows the synthesis phase map of MnS_1−xSe_x alloys as a function of synthesis temperature and chemical composition, based on the X-ray diffraction (XRD) and X-ray fluorescence (XRF) data. As expected, the MnSe and MnS end-member compounds deposited as reference samples exhibit diffraction patterns corresponding to the ground state RS structure at 280–420 °C. However, their crystal structure changes when they form MnS_1−xSe_x alloys at lower temperatures within this range. At 420 °C, the MnS_1−xSe_x alloys have the RS structure across the entire composition range, the same as the end-member compound, and similar to other isostructural alloys (e.g., WZ Al_xGa_1−xN [17] or ZB CdTe_1−xSe_x) [18]. At 340 °C, the MnS_1−xSe_x alloys still exhibit peaks indicative of the RS structure through the whole composition range, but an additional signal from a WZ phase is observed at S-rich compositions (x < 0.3). At 280 °C, a similar combination of RS + WZ is observed for 0.3 < x < 0.4, but for x < 0.3, the RS signal disappears completely, and all XRD peaks can be attributed to the WZ-phase. Figure 1(b) shows the XRD heat map for MnS_1−xSe_x alloys synthesized at 280 °C, and similar figures for films grown at 420 °C and 340 °C that were used to populate points in the phase map can be found in Fig. S1 in the Supplementary material (SM). In the past, similar metastable WZ-MnSe_1−xTe_x [14] alloys have been shown to decompose into stable parent compounds close to 400 °C. Overall, these results indicate that the WZ phase can be obtained under nonequilibrium growth conditions at lower temperatures, where the phase transformation into the RS ground state is kinetically hindered.

Figure 1: (a) Phase map of MnS_1−xSe_x alloys as a function of chemical composition and synthesis temperature. (b) XRD heat map of MnS_1−xSe_x alloys at 280 °C, showing only WZ peaks appearing at x < 0.3 composition.

To support the WZ phase assignment, a sample with MnS_0.92Se_0.08 composition synthesized at 280 °C was chosen as the representative for more detailed characterization. Figure 2(a) shows the XRD pattern for this sample as a function of 2Θ angle, with χ angle integration from −126° to −54° (detector at −90°), and a fixed φ angle. The WZ phase shows strong (002) preferential orientation, with only one prominent peak at 2Θ ∼ 30°, similar to WZ MnS thin films reported in the literature [19]. Comparison of XRD patterns for WZ-MnS_0.92Se_0.08 with RS-MnS and RS-MnSe is shown in Fig. S2(a) in SM. To further investigate the phase assignment, an XRD pole figure was collected on a different instrument as a function of φ (0–360°) and χ (0–90°) for the measured 2Θ peak position at ∼30°. Figure 2(b) shows that the peak intensity is strongest at χ = 90°, where the incident and scattered X-ray plane is normal to the substrate surface, further confirming the c-axis preferential orientation of WZ thin films. Similar full-range XRD measurements were conducted for possible 2Θ positions of the RS phase, but none of the other peaks were detected, ruling out the presence of the RS phase in MnS_0.92Se_0.08 synthesized at 280 °C.

Figure 2: (a) XRD pattern as a function of 2Θ for MnS_0.92Se_0.08 synthesized at 280 °C, showing strong preferential orientation of WZ phase. (b) Full-range XRD scan as a function of φ and χ at 2Θ ∼ 30° WZ peak for MnS_0.92Se_0.08 deposited at 280 °C.

The existence of the WZ phase in MnS_1−xSe_x alloys is also supported by selected area electron diffraction (SAED) measurements. Figure 3(a) shows the measured SAED pattern for MnS_0.92Se_0.08 synthesized at 280 °C and Fig. 3(b) is the simulated reference pattern based on the [100] zone axis of the WZ and RS phase. The measured SAED pattern matches well with the simulated WZ pattern but is different from the RS pattern, and unambiguously confirms the formation of the WZ phase at this composition in MnS_1−xSe_x alloys. However, minor splitting appears for the (0112) peak family, which may be due additional WZ grains with different orientations, or due to minor RS impurity phase. A transmission electron microscopy (TEM) image of this sample point is also provided as Fig. S3 in SM, to show a cross-sectional view of the WZ phase with ∼100-nm grain size that may be possible to increase by annealing. The estimate film thickness from the cross-sectional TEM image is 150 nm, which is similar to 140 nm determined from the XRF model.

Figure 3: (a) SAED patterns of MnS_0.92Se_0.08 grown at 280 °C. (b) Simulated WZ [100] and RS [100] reference pattern from a single crystal.

Computational results

In an attempt to understand the origin of the appearance of the WZ phase between RS end-member compounds, 0 K mixing enthalpies of MnS_1−xSe_x alloys as function of x (∆H _m(x)) were calculated from first principles, and are shown in Fig. 4(a). Generally, the mixing enthalpy of the RS phase is lower than that of the WZ phase along the entire composition range, predicting that the alloys should assume the RS phase under equilibrium conditions at 0 K. For the typical synthesis temperature (280–420 °C), the difference in entropic contribution to the free energy between two solid phases (WZ and RS) is likely to be smaller than the calculated difference in enthalpies of mixing (15–45 meV/at). These computational results are consistent with the measured RS crystal structure of MnS_1−xSe_x alloys synthesized at the higher temperature of 420 °C [Fig. 1(a)]. However, at 280 °C, the WZ phase is experimentally observed (Fig. 1) at sulfur-rich compositions. At this lower temperature, formation of the stable RS structure may be kinetically hindered, facilitating formation of metastable WZ structures of the MnS_1−xSe_x alloys. It is particularly surprising that the WZ phase was observed at S-rich compositions [Fig. 1(b)], where the calculated WZ-to-RS polymorph energy difference is the highest [Fig. 4(a)].

Figure 4: (a) Calculated mixing enthalpy as a function of x in MnS_1−xSe_x. (b) Calculated volume per formula unit as a function of x in MnS_1−xSe_x.

Synthesis of the higher energy WZ polymorph in MnS_1−xSe_x alloys by co-sputtering from RS-MnS and RS-MnSe targets (Fig. 1) is quite surprising in light of our theoretical calculations [Fig. 4(a)], but not uncommon in a broader sense. In general, preferential growth of WZ structure has been observed across different material chemistries (I–VII, II–VI, and II–V) synthesized by a wide range of techniques (solid, liquid, and vapor growths). For example, higher energy WZ polymorphs of I–VII materials (LiBr [20] and InCl [21]) have been synthesized by crystallization from amorphous precursors deposited by thermal evaporation on cooled substrates. As another example, preference toward WZ polymorph stabilization at low crystal sizes has been observed in pressure-dependent studies of WZ-to-RS transformations of II–VI (CdSe) nanocrystals prepared by colloidal synthesis techniques [22, 23]. Similar preference of the WZ over the ZB polymorph is also known in vapor–liquid–solid (VLS) growth of III–V (GaAs [24] and InAs [25]) nanowires. This broad range of materials and methods suggests a universal mechanism of WZ polymorph preference over other possible structures.

We speculate that the formation of the WZ phase at S-rich compositions of co-sputtered MnS_1−xSe_x alloys reported here may be related to Ostwald's “step rule” [26]. This rule is based on empirical observations that during the crystallization process metastable structures often form before stable phases. More recently, this rule has been rationalized by measurements of lower surface energies of the metastable polymorphs compared to stable structures [27]. In the case of MnS_1−xSe_x alloys reported here, nucleation of the T _d-bonded WZ metastable polymorph may be related to its lower surface energy relative to the O _h-bonded polymorph. Such lower surface energy is likely because of the lower number of broken bonds at the surface of a growing particle, but explicit theoretical calculations would be needed to evaluate this hypothesis. This relation to surface energy is also suggested by the observed strong preferential orientation of the MnS_1−xSe_x alloys along (002) WZ and (111) RS growth directions (Fig. 2), which are structurally related and highly polar. More calculations would be required to rule out potential other hypothesis, such as defect-induced stabilization of the WZ polymorph observed here.

Optoelectronic properties

Optical properties of MnS_1−xSe_x were also investigated, both experimentally and theoretically. Figure 5 shows the measured absorption heat map of MnS_1−xSe_x as a function of x where an absorption spectrum at each x is displayed as a false color map. The theoretically predicted band gaps of WZ (direct) and RS (indirect) phases are shown as markers for comparison. At 420 °C, the alloys assume the RS phase across the entire composition range [Fig. 5(a)]. Both the predicted and experimental band gaps decrease nearly linearly as x increases, indicating the weak interaction between elemental constituents in RS MnS_1−xSe_x alloys. This linear change may help to fine-tune optical properties of RS MnS_1−xSe_x alloys, in contrast to the band gap bowing phenomenon commonly observed in other semiconductors due to the alloying interaction effect [16]. The measured absorption onset (defined here as an energy at which the absorption coefficient is 10⁵ cm⁻¹) is comparable between WZ samples [Fig. 5(b)] and RS samples [Fig. 5(a)] at the same compositions, both in the 3.0–3.2 eV range. The experimental and computational optical absorption onset trends appear cross, as the composition is changed from MnS to MnSe, but quantification of these differences is difficult, given the present data. Comparison of optical absorption spectra for WZ–MnS_0.92Se_0.08 with RS–MnS and RS–MnSe is shown in Fig. S2(b) in SM. Despite the similar absorption onset, WZ alloys have direct optical band gaps, whereas RS alloys have indirect gaps. The direct band gap of the WZ phase enables its potential application in optoelectronic devices.

Figure 5: Optical absorption spectra as a function of x in MnS_1−xSe_x alloys synthesized at (a) 420 °C and (b) 280 °C. The discrete markers are the predicted band gaps of MnS_1−xSe_x alloys in the WZ and RS structures.

Electrical conductivities of the MnS_1−xSe_x alloys synthesized at 420 and 280 °C are summarized in Fig. 6(a). The conductivities of the MnS_1−xSe_x alloys with a WZ structure are lower than the measurement limit of the 4-point probe instrument; thus, it is estimated to be below 0.0001 S/cm. For the MnS_1−xSe_x alloys with a RS structure, the conductivities increase from 0.0001 S/cm to 1 S/cm with increasing x, suggesting that MnSe has shallower defects than MnS. Computationally predicted effective masses of electrons and holes in MnS_1−xSe_x alloys are listed in Fig. 6(b). In MnS_1−xSe_x alloys with a RS structure, the effective masses of electrons and holes both decrease with x, which partly explains the increase of conductivities in Fig. 6(a). For the MnS_1−xSe_x alloys with a WZ structure, both effective masses are relatively insensitive to chemical composition x, and hole effective mass is quite scattered. This scatter may be related to finite size effect in supercell-based alloy calculations (Fig. S4), which may be amplified in WZ valence band due to p–d coupling effects [28]. It is also noteworthy that the effective mass of electrons in the WZ structure (0.3m ₀ at x = 0.2) with a direct band gap is much lower than that of electrons in the RS structure (1.5m ₀ at x = 0.2) with an indirect band gap, suggesting potential optoelectronic applications of WZ alloys if they could be doped n-type. Such differences in effective masses between tetrahedrally coordinated (WZ and ZB) and octahedrally coordinated (RS and NC) polymorphs have been reported in other binary chalcogenides [8, 14] and in ternary nitrides [29, 30, 31].

Figure 6: (a) Measured electrical conductivity of MnS_1−xSe_x alloys with the RS structure as a function of x. (b) Computationally predicted density of states effective masses (m*) of electrons and holes in MnS_1−xSe_x alloys with RS and WZ structures.

Summary

The structural and optoelectronic properties of MnS_1−xSe_x alloys (0 ≤ x ≤ 1) as a function of chemical composition and synthesis temperature are reported. Despite the fact that the end-member compounds MnS and MnSe assume RS as the experimentally observed and calculated ground state structure, a metastable WZ structure appears at S-rich compositions grown at a lower temperature of 280 °C. This is in contrast to higher temperature growths at 420 °C where the alloys remain in the RS over the entire composition range. The stabilization of the WZ-phase at lower deposition temperatures is surprising in light of theoretical calculations that predict the mixing enthalpy for WZ structures to be higher than for RS structures at all compositions. The stabilization of WZ polymorph may be attributed to lower surface energy of the WZ structure compared with the RS structure, as observed in other materials. In addition, the higher density of RS structure may account for the preference of the formation of WZ phase at S-rich compositions. Computationally predicted band gaps of MnS_1−xSe_x alloys in RS and WZ structures correlate well with the experimentally measured optical absorption onsets. In addition, the RS MnS_1−xSe_x alloys are measured to be more electrically conductive for the Se-rich composition, and calculated electron and hole effective masses are consistent with this observation.

Methods

Experiments

MnS_1−xSe_x (0 ≤ x ≤ 1) alloys were prepared by combinatorial RF sputtering (AJA International, Orion 8). Two sputter targets with a diameter of 50.8 mm, one of MnSe and one of MnS, were placed opposite to each other in confocal geometry to generate a composition gradient covering the entire composition range in 3–4 depositions. The sputtering power was 8–50 W. For end-member samples, only a single source was used. Fused silica was used as substrates and partially mounted on the sample platen by silver paint. This partial contact induces a thermal gradient across the substrate orthogonal to the composition gradient. The sample platen and substrates were heated by a SiC heater to achieve the maximum substrate surface temperatures of 420 °C, 340 °C, and 280 °C. The chamber was pumped down to a base pressure of <1 × 10⁻⁶ torr before thin films were deposited in high-purity Ar at a working pressure of 1.3 × 10⁻³ torr and an Ar flow rate of 16 sccm. The target–substrate distance was 134 mm. More detailed synthesis information can be found in our previous papers [32, 33].

After synthesis, the MnS_1−xSe_x thin-film libraries were characterized as a function of position at 44 locations on a rectangular 4 × 11 grid. Thin-film properties obtained from these mapping techniques included: composition and thickness measurement by X-ray fluorescence (XRF; Fischerscope XDV-SDD), crystal structure measured by X-ray diffraction (XRD; Bruker D8) with Cu K_α radiation and a 2D detector in an approximate Bragg-Brentano geometry, optical transmittance and reflectance measured by a custom UV-Vis instrument using ocean optics spectrometers covering wavelengths of 300–1100 nm, and sheet resistance collected by a custom 4-point probe instrument [34]. The XRF measurements of MnS_1−xSe_x alloy composition were verified by measuring MnS and MnSe parent compounds with known composition, and the XRF-derived thickness was verified by a cross-sectional TEM measurement. The absorption spectra and electrical conductivities were determined from these measurements based on thickness data obtained from XRF. The large amounts of data that resulted from spatially resolved characterization were analyzed using the COMBIgor software package (https://www.combigor.com/) [35] and will be made publicly available at HTEM DB (https://htem.nrel.gov/) [36].

The X-ray pole figure of the WZ (002) peak was collected on a Rigaku Smartlab diffractometer using Cu-K_α radiation and at 2Θ = 30°. TEM was performed on a FEI Tecnai ST30 at 300 kV. A Pt capping layer was deposited on the surface of the samples to protect the thin-film structure. TEM lamella were prepared by the standard focused ion beam (FIB) lift-out technique and followed by Ga ion milling to reduce the final thickness to less than 100 nm. Then, a Fischione Nanomill was used to remove preparation damage on the TEM foils using low energy cleaning at 500 eV under vacuum of 10⁻⁷ torr at −170 °C for ±10°.

Computations

Periodic boundary condition density functional theory calculations were performed within the Vienna Ab Initio Simulation Package [37] (VASP) using projector augmented-wave (PAW) pseudopotentials [38]. To obtain accurate polymorph energies of MnS and MnSe polymorphs, we used the random phase approximation based on the adiabatic-connection fluctuation–dissipation theorem (RPA-ACFDT) [39]. For RPA calculations of MnS, the energy cutoff for the plane-wave basis was 360 eV while the energy cutoff for the response function was 240 eV with a 4 atom unit cell. DFT + U with U = 3 eV for Mn-d orbitals was utilized as the input wave function for obtaining the Hartree Fock energy and the ACFDT-RPA correlation energy. The energy difference with respect to the magnetic ground state was determined at the DFT + U level and subtracted from the RPA energy. Polymorph energies for MnSe were calculated in a previous publication in a similar way [14].

For calculations of MnS_xSe_1−x alloy energetics, we used the Perdew–Burke–Ernzerhof (PBE) generalized gradient approximation exchange-correlation functional with a Hubbard correction of U = 3 eV for Mn-d orbitals [40]. Mixing enthalpies were calculated using a 96 atom supercell with random cation configurations generated using the special quasirandom structures method [41]. Calculations were conducted at 17 compositions across the phase space for each crystal structure. Calculations utilized a Gamma centered 2 × 2 × 2 k-point mesh with an energy cutoff of 400 eV. Alloy interaction parameters (Ω) for each polymorph were determined by fitting ΔH _m(x) using a second-order polynomial.

The theoretical band gaps and DOS effective masses of the alloys were calculated through a linear weighting of the DFT-predicted values to the end-member calculated within the GW approximation [42]. For density of states effective mass calculations, a denser k-mesh with 1000 k-points per reciprocal atom was used, and the results obtained from an energy window corresponding to effective temperatures of 500–3000 K. Correcting the GW self-energy with an attractive on-site potential for the d orbitals has been shown to result in accurate band gaps for transition metal compounds which are generally underestimated by DFT, as described in detail in prior publications [43]. However, the GW method is computationally prohibitive for the large supercells required for the alloys studied here. Therefore, the GW-computed values of the end-members were used to correct the DFT values for the alloys across the composition space. The computational results of the GW approximation are available from the NREL materials database (http://materials.nrel.gov) [43, 44]. The “effective GW” band gaps and effective masses for alloys of MnS_1−xSe_x are given by a simple linear weighting: $E_x^{{\rm{GW}}} = E_x^{{\rm{DFT}}} \left( {x{{E_{x = 1}^{{\rm{GW}}} } \over {E_{x = 1}^{{\rm{DFT}}} }} + {{\left( {1 - x} \right)E_{x = 0}^{{\rm{GW}}} } \over {E_{x = 0}^{{\rm{DFT}}} }}} \right)$, as described in previous publications [45]. This approach is reasonable because the band gaps and effective masses are underestimated by DFT, but the changes in the electronic structure across the composition range are generally correct in these materials.

Supplementary material

To view the supplementary material for this article, visit https://doi.org/10.1557/jmr.2019.402.