A.
Portone
a,
L.
Bellucci
a,
D.
Convertino
bd,
F.
Mezzadri
fg,
G.
Piccinini
bc,
M. A.
Giambra
e,
V.
Miseikis
bd,
F.
Rossi
f,
C.
Coletti
bd and
F.
Fabbri
*a
aNEST, Istituto Nanoscienze – CNR, Scuola Normale Superiore, Piazza San Silvestro 12, 56127 Pisa, Italy. E-mail: filippo.fabbri@nano.cnr.it
bCNI@NEST, Istituto Italiano di Tecnologia, Piazza San Silvestro 12, 56127 Pisa, Italy
cScuola Normale Superiore, Piazza San Silvestro 12, 56127 Pisa, Italy
dGraphene Labs, Istituto Italiano di Tecnologia, Via Morego 30, 16163 Genova, Italy
eCNIT, Sant'Anna, Via G. Moruzzi 1, Pisa, 56124, Italy
fIMEM-CNR, Parco Area delle Scienze 37/a, Parma, 43124, Italy
gDepartment of Chemistry, Life Sciences and Environmental Sustainability, University of Parma, Parco Area delle Scienze 11/A, 43124 Parma, Italy
First published on 2nd February 2021
The employment of two-dimensional materials, as growth substrates or buffer layers, enables the epitaxial growth of layered materials with different crystalline symmetries with a preferential crystalline orientation and the synthesis of heterostructures with a large lattice constant mismatch. In this work, we employ single crystalline graphene to modify the sulfurization dynamics of copper foil for the deterministic synthesis of large-area Cu9S5 crystals. Molecular dynamics simulations using the Reax force-field are used to mimic the sulfurization process of a series of different atomistic systems specifically built to understand the role of graphene during the sulphur atom attack over the Cu(111) surface. Cu9S5 flakes show a flat morphology with an average lateral size of hundreds of micrometers. Cu9S5 presents a direct band-gap of 2.5 eV evaluated with light absorption and light emission spectroscopies. Electrical characterization shows that the Cu9S5 crystals present high p-type doping with a hole mobility of 2 cm2 V−1 s−1.
Graphene has also been widely used as a growth substrate for the synthesis of different classes of materials and nanostructures. In particular, graphene has been employed as a substrate in the synthesis of InAs nanowires with preferential orientation of NW facets9,10 or in the production of double-heterojunction InAs/Gr/InAs nanowires.11 Similar results were obtained in the synthesis of ZnO nanorods on graphene resulting in a high vertical alignment with an epitaxial relationship between ZnO and graphene.12 Recently a new approach has been demonstrated in the growth of 3D crystalline materials using a 2D material as a buffer layer in the growth process.13–15 This process has been successfully employed in the development of low dislocation density GaN based light emitting diodes (LEDs) on graphene coated ZnO.16 Similarly, AlGaN, a deep-UV LED with promising light emission properties, has been realized on high quality AlN/Gr/AlN structures, where graphene plays a fundamental role in the synthesis of a high quality AlN substrate.17 A similar strategy has been developed for the heteroepitaxy of III–V materials as the growth on graphene-coated substrates enables spontaneous relaxation of the epilayer before the accumulated elastic energy created by the lattice mismatch induces misfit dislocations. This approach has been demonstrated on low misfit strain InGaP on GaAs (0.74% misfit strain) and on high misfit strain GaP on GaAs (3.7% misfit strain).14
In the case of the synthesis of 2D layered materials or three-dimensional (3D) crystalline layers on top of layered materials this particular method has been identified as van der Waals epitaxy (vdWE) by Koma.18
vdWE enables the epitaxial growth of layered materials with different crystalline symmetries to the substrate. In addition, this method can be used even when the lattice constant mismatch is as high as 50% and no excessive strain exists in the epitaxial layer. The main advantages of vdWE include rotational alignment with the substrate, strain-free growth, and no misfit dislocations despite the significant lattice mismatch.19
The vdWE approach includes two different configurations. The first is the growth of 2D layered materials onto 3D crystal substrates, which have saturated dangling bonds on their clean surfaces. Then, a single-crystalline film of a layered material can be heteroepitaxially grown on such a dangling-bond-terminated surface via van der Waals-like weak force. For example, layered material (NbSe2, MoSe2, GaSe, InSe) films have been grown on such saturated dangling bond surfaces as H-terminated Si(111), S or Se-terminated GaAs(111), and F-terminated CaF2(111).19 The second layout of the vdWE is the synthesis of a layered material onto another layered substrate, which has no dangling bond on its cleaved surface. This approach has been widely employed in the growth of transition metal dichalcogenides (TMDs) on top of graphene20,21 or mica substrates.22,23 For instance, the CVD synthesis of WS2 monolayers on top of epitaxial graphene or of CVD graphene has demonstrated a preferential orientation between 0° and 60° of monolayer WS2.24–27 Another example of this configuration is the synthesis of a mono- or few layer metal monochalcogenide on top of another TMD material:19 the CVD synthesis of monolayer GaSe on monolayer MoSe2, despite a large lattice misfit, results in bilayer heterojunctions with atomically sharp interfaces.28 Even in the case of two-dimensional monoelemental materials the vdWE method has been successful in the synthesis of few layer antimonene29 or tellurium30 on mica and silicene on graphite.31,32 In this work, we utilize graphene single-crystals to modify the sulfurization dynamics of copper foil for the deterministic synthesis of large-area Cu9S5 crystals. We are able to synthesize Cu9S5 flakes (also known as digenite) with a lateral size of more than 50 μm. The graphene crystals are completely etched during the growth process. Comprehensive analysis of the optical and electrical properties is carried out after dry stamp transfer of the Cu9S5 crystals. In particular, the evaluation of the optical band-gap has been carried out on the crystals revealing that the Cu9S5 has an optical direct band-gap of about 2.5 eV. The electrical characterization shows that the Cu9S5 crystals present high p-type doping with a hole mobility of 2 cm2 V−1 s−1.
Graphene single crystal arrays with a lateral size of about 100 μm are synthesized via CVD on electropolished copper (Cu) foils (Alfa Aesar, 99.8%) following the procedure described by Miseikis et al.33,34 Graphene was then synthesized at a temperature of 1060 °C inside a cold-wall CVD system (Aixtron BM) flowing methane, hydrogen, and argon.
Scanning electron microscopy (SEM) was performed to study the morphology of Cu2S flakes. The SEM analysis is carried out using an in-lens detector on a Zeiss Merlin FEG-SEM microscope, with an accelerating voltage of 5 kV and an electron current of 96 pA, in order to minimize the damage due to the electron beam irradiation.
AFM analysis was performed with an AFM+ from Anasys Instruments, operated in tapping-mode. AFM micrographs were analysed using the software Gwyddion 2.45.
Raman spectroscopy and mapping were carried out with a Renishaw InVia system, equipped with a confocal microscope, a 473 nm excitation laser and a 2400 line per mm grating (spectral resolution 0.5 cm−1). All the analyses were performed with the following parameters: excitation laser power, 0.5 mW and acquisition time for each spectrum, 10 s. The maps were obtained with the same acquisition parameters and a pixel size of 2 μm × 2 μm. PL spectroscopy was performed with the same system using the 473 nm excitation laser. All the spectra were obtained with the following parameters: excitation laser power, 0.5 mW and acquisition time for each spectrum, 30 s.
Ultraviolet-visible (UV-Vis) absorbance spectra were collected from Cu9S5 flakes transferred on a thick PDMS film by using a Lambda 950 spectrophotometer (Perkin Elmer Inc.) equipped with deuterium and tungsten halogen light sources, a double holographic grating monochromator and a high-sensitivity photomultiplier tube. Spectra were acquired in the wavelength range from 250 nm to 650 nm with a spectral resolution of 1 nm.
Energy Dispersive X-ray (EDX) spectroscopy and mapping were performed on a JEOL JEM-2200FS Transmission Electron Microscope operated at 200 kV in scanning mode (STEM, spot size 2 nm), using an Oxford XPLORE detector with drift correction. The samples were deposited by drop-casting on carbon coated nickel grids.
I–V curves and Hall measurements were carried out at room temperature with an Ecopia HMS-3000 Hall measurement system. Cu9S5 flakes are brought into contact by means of silver paste and the tip of the manipulator is positioned on the PDMS surface covered with a drop of silver paste. Due to the extreme brittleness of digenite, the possible force applied by the manipulator can break the flake. The I–V curves are obtained applying a bias from −0.5 V to 0.5 V with a step of 50 mV. In the case of Hall measurements, the probe current was selected such that the voltage signal was 0.2 V with a magnetic field of 0.37 T.
Powder X-ray diffraction data were acquired using a Rigaku Smartlab XE diffractometer. The CuKα incident beam was parallelized through a parabolic mirror in order to limit the sample roughness effects on the peak shape and position. Measurements were carried out in the Bragg–Brentano configuration with 5° soller slits both on the incident and diffracted beam. An HyPix3000 detector was operated in 1D mode. Samples were placed on a zero-background sample holder and data were collected in the 10°–80° 2θ range with 0.01° steps and 2° min−1 speed.
Visual inspections of MD trajectories and figures were performed with VMD1.9.3.41
Fig. 1d shows the SEM image of a 4 × 4 array of Cu9S5, as identified by Raman analysis. Indeed, it is clear from the micrograph that the graphene crystals assist the growth of Cu9S5, as the hexagonal shape of the graphene flake is transferred to the Cu9S5compounds. It is possible to note the presence of cracks outside the array related to the extreme brittleness of digenite, as previously reported.50 The SEM image of a single flake is reported in Fig. 1e. The high magnification SEM image of a Cu9S5 flake reveals that the areas where the graphene crystal is present are recessed with respect to the surrounding dendritic areas. The SEM analysis reveals that the flakes are composed of grains with different sizes. In order to evaluate the morphology of the grains comprising the Cu9S5 flake, atomic force microscopy (AFM) is carried out. The average root mean square (RMS) of the whole area is 700 pm meanwhile the micrometric areas, highlighted by the white square, present a lower RMS (340 pm). The Cu9S5 flake shows a flat surface in comparison with CuXS compounds, obtained with similar strategies.50 The line profile of nanometric pits is reported in Fig. S5.† The morphological analysis suggests that the presence of graphene crystals affects the growth direction of the underlying Cu9S5, resulting in the alignment of the basal plane of graphene and of Cu9S5, albeit the graphene is completely removed, by chemical etching or mechanical detachment, during the growth process. This hypothesis is supported by the extremely low RMS and by the recessed morphology of the Cu9S5 flakes. In fact, the raised dendritic morphology of the surrounding areas outside the flakes is probably due to a random orientation of the Cu9S5 and CuS phases. A similar effect was previously reported in the case of the sulfurization of copper thin films, where the increase of the copper thickness results in the synthesis of vertically aligned β-Cu2S lamellas.51
Until now the synthesis of Cu9S5 has been carried out mainly with solution processes.52–54 Recently a novel method, similar to that presented in this work, is reported for the high yield synthesis of bulk Cu9S5 by vapour phase sulfurization of copper foil.50
To characterize the formation of Cu9S5 flakes at the molecular-atomistic level, molecular dynamics (MD) simulations were performed considering different model systems composed of a copper Cu(111) surface exposed to a sulphur/argon mixture or screened from the S/Ar atoms by single layer graphene or by graphene ribbons (see Fig. 2 and S6†). For this purpose, we used the empirical reactive force-field55 (ReaxFF) as implemented in the LAMMPS software package,56,57 exploiting the recent parametrization developed for Cu/S/C/H.36 The use of such a force-field allows us to empirically model reactive events occurring during the sulfurization process through a bond-order formalism. All the systems were treated as periodic by applying a periodic boundary condition scheme (see the Experimental section). In this framework, we investigated the role of graphene on top of the Cu(111) surface during the initial stage (about 1 ns) of the sulfurization process that eventually leads to the synthesis of Cu9S5. All simulations were initialized and relaxed at constant temperature (1000 K) for 1 ns employing the same computational protocol (see methods). The relaxed structures were then optimized through a minimization protocol (see methods) and are shown in Fig. 2. The trajectories of the MD simulations were then inspected and used to grasp the response of the systems to the chemical attack of sulphur (S). The Cu(111) surface completely covered by graphene during sulphur attack is stable (Fig. 2a). The binding energy (Eb) per carbon atom, i.e. the interaction energy between graphene and the metal surface was about −40 meV, which is in agreement with the values evaluated at the density functional level of theory;40 graphene adheres to the surface and being chemically inert protects the Cu atoms from the chemical attack. Sulphur atoms in excess remain in the box simulation, occasionally forming small chains or cyclic molecules.58,59 The uncovered Cu(111) is attacked by S atoms during the MD simulation (Fig. 2b). The Cu atoms of the first layer are extruded from the surface and react with chemisorbed sulphur atoms. Reactions can also occur among sulphur atoms, leading to the formation of small S chains and S–S bonds interconnecting Cu atoms (see Movie 1†). This is likely disadvantageous for further S attack and the diffusion of S into the inner layers. The formation of these sulphur chains reminds of the structure of copper mixed monosulphide/disulphide (e.g., CuS or CuS2) that differs from digenite, which contains only isolated sulphide anions, by the presence of direct S–S bonds (Fig. S7†). The behaviour of the surfaces covered with graphene ribbons is consistent with that observed in the two previous cases. Graphene ribbons, protecting Cu from the S attack, remain attached to the metal surface (Fig. 2c). During the simulation, some sulphur atoms approach the hydrogenated edges and partially intercalate between the graphene/Cu interface. However, no sulphur atoms crossed the interface, as sulphur intercalation is hindered by the presence of graphene adhering to the surface. This means that the intercalation process is a slow process regulated by an activation energy (e.g., the energy due to the partial detachment and/or large deformation of the edges of the graphene) which prevents its observation in the time scale of standard MD simulations.60 The relaxed structures of two systems with different numbers of intercalated sulfur atoms are shown in Fig. 2d and in Fig. 2e. During the simulation large rearrangements of Cu and S atoms occur. Graphene remains attached to the surface, eventually driving the intercalation of single sulphur atoms over the Cu(111) surface. The presence of graphene prevents the intercalation of sulphur chains and shields the intercalated S atoms from the S/Ar mixture, disadvantageous for uncontrolled formation of S–S bonds or further S chains (see Movie 2†). Therefore, during the sulfurization process, graphene promotes the formation of a more ordered morphology of the surface (compare Fig. 2b with Fig. 2d and e). Remarkably, during the relaxation of the system with the graphene covered Cu(111) surface, no sulphur chains and S–S bonds form and the sulphur atoms are evenly distributed among Cu atoms of the metal surface. The bond energy for the carbon atoms decreases as the number of S increases, being −17 meV for the system shown in Fig. 2d and −19 meV for the system in Fig. 2e. The intercalation of sulphur therefore favours the detachment of graphene. However, during the MD simulations it was not possible to observe the release of graphene due to the pressure of the gas phase (i.e. the Ar/S mixture) which acts on the periodic surface of the graphene, which in this specific case prevents its detachment.
In order to analyse the Cu9S5 flakes we employ a two-step PDMS stamp process in order to transfer only the flakes. The transfer process (Fig. 3a) is based on a PDMS thick film that is applied to the surface of the substrate and gently pressed. Due to the extreme brittleness of digenite50 and the recessed morphology of the Cu9S5 flakes when the PDMS film is peeled off, the surrounding material remains attached to the film while the flakes remain on the residual copper foil substrate. With a second step, similar to the first one, it is possible to detach the flakes from the residual copper foil. Then, the flakes are analysed on the PDMS substrate or they are detached by gentle sonication in isopropanol. The optical image of the transferred flakes is reported in Fig. S8 of the ESI.† In order to evaluate the thickness of the Cu9S5 flake, we measure the thickness of the pristine copper foil, after the synthesis of graphene, and after the transfer of the sulfurized layer. The pristine copper foil is 25 μm thick, and it decreases to 8 μm due to the electrochemical polishing, necessary for the graphene synthesis. The thickness of the residual copper foil, obtained by stylus profilometry, decreases down to 5 μm after the sulfurization process, revealing that the thickness of the transferred sulfurized layer is about 3 μm. The flakes can be transferred on a nickel TEM grid in order to study the structural and compositional properties. The low magnification TEM image is reported in Fig. 3b, revealing that the flakes are damaged during the transfer process and they are not electron transparent due to the micrometric thickness. The TEM-EDX spectrum (Fig. 3c) presents the characteristic S and Cu peaks. Additional peaks are due to the nickel of the TEM grid and the Fe and Co of the TEM sample holder; no other impurity elements are detected. Elemental maps (Fig. 3d) show a continuous distribution of S and Cu, confirming that the growth process allows a good spatial uniformity of the sample.
Fig. 3 (a) Schematic of the two-step PDMS stamp transfer, (b) low magnification TEM image. (c) EDS spectrum and (d) EDS maps of copper (cyan) and sulfur (red). |
After the transfer onto the PDMS thick film (inset of Fig. 4a), we carried out optical and electrical characterization of the Cu9S5 crystals. In general, the bandgap of CuXS compounds can vary from 1.2 eV to 2.5 eV, depending on different factors, such as the crystalline phase, the stoichiometry, the nanostructure.61 In particular, bulk Cu9S5 (digenite phase) has a direct bandgap of 2.5 eV,50 however the optical bandgap of nanomaterials is deeply linked with their morphology, in fact the reported bandgap of Cu9S5 varies from 1.73 eV (film)52 to 1.93 eV (nanocube)62 and 2.1 eV for 2D triangular flakes.43
Fig. 4 (a) Tauc plot of Cu9S5 flakes on the PDMS thick film. (b) PL spectrum. (c) I–V curve of the Cu9S5 flake, the inset shows the sketch of the acquisition system. |
The absorbance spectra of Cu9S5 crystals on the PDMS thick film are presented in Fig. S9.† The optical absorption data were analysed to determine the optical bandgap values using the Tauc relation for the direct allowed electronic transition.63Fig. 4a shows the Tauc plot. It is possible to extrapolate the optical indirect band-gap of Cu9S5 flakes, with two different procedures. One is the standard intercept method following the Tauc relation63
αhν = C(hν − Eg)2 | (1) |
The light emission properties are studied by means of photoluminescence (PL) spectroscopy. The PL spectrum of Cu9S5 flakes, excited with a 473 nm laser line (Fig. 4b), presents an intense emission at 2.48 eV with a weak shoulder on the low energy side. The spectrum can be deconvolved evaluating the position of the single emissions, resulting in two main peaks set at 2.48 eV (500 nm) and 2.13 eV (557 nm). The high-yield sharp emission at 2.48 eV is in good agreement with the direct optical band-gap evaluation of 2.5 eV reported in Fig. 3a. Therefore, it is possible to attribute the broad light emissions at 2.13 nm to intra-band-gap states related to point or extended defects. Nevertheless, light emissions with similar energies in 2D-β-Cu2S flakes were attributed to the band-to-band emission of the different spurious CuXS phases.51 In our case the Raman analyses performed within the hexagonal flakes do not reveal the presence of any other crystalline structure than Cu9S5, ruling out the possible contributions of the band-to-band radiative recombination of other CuXS phases. In addition, the band-gap evaluation is important to study both the optical and the electrical properties, in fact CuXS is widely known to be a p-type semiconductor where the hole doping is due to the intrinsic concentration of the copper vacancies. It has been reported that the CuXS band-gap can increase due to the Burstein–Moss effect caused by the high concentration of charged defects.64Fig. 4c shows the I–V curves of the Cu9S5 flakes. In the inset of Fig. 4c a schematic representation of the acquisition system is illustrated. The I–V curves present a strong non-ohmic behaviour. This particular feature can be due to multiple Schottky diodes and resistors in parallel, with slightly different breakdown voltages and resistances, changing the shape of the I–V curve towards the non-ohmic behaviour. In addition, tunnelling currents, related to the high carrier concentration, can cause the non-ohmic behavior.65 The majority carrier sign, hole mobility and the net carrier concentration are determined, remaining in the linear range of current using a 0.2 V bias, by room-temperature Hall measurements. The conductivity measured using a 4-point probe on bulk flakes was very high, about 4.8 (Ω cm)−1. Cu9S5 flakes present high p-type doping with a carrier concentration of (8 ± 1) × 1019 cm−3 and a hole mobility of 2.0 ± 0.5 cm2 V−1 s−1. These data are in agreement with the doping concentration and mobility of bulk Cu9S5.50
The high doping concentration is due to the high concentration of copper vacancies,66 which makes the Cu9S5 flakes suitable as hole transport layers in perovskites52 or dye sensitized53 solar cells and in MoS2 based photodetectors.43
Footnote |
† Electronic supplementary information (ESI) available. See DOI: 10.1039/d0na00997k |
This journal is © The Royal Society of Chemistry 2021 |