Abstract
Thermally induced structural transformation of 2D materials opens unique avenues for generating other 2D materials by physical methods. Imaging these transitions in real time provides insight into synthesis routes and property tuning. We have used in situ transmission electron microscopy (TEM) to follow thermally induced structural transformations in layered CoSe2. Three transformation processes are observed: orthorhombic to cubic-CoSe2, cubic-CoSe2 to hexagonal-CoSe, and hexagonal to tetragonal-CoSe. In particular, the unit-cell-thick orthorhombic structure of CoSe2 transforms into cubic-CoSe2 via rearrangement of lattice atoms. Cubic-CoSe2 transforms to hexagonal-CoSe at elevated temperatures through the removal of chalcogen atoms. All nanosheets transform to basal-plane-oriented hexagonal 2D CoSe. Finally, the hexagonal to tetragonal transformation in CoSe is a rapid process wherein the layered morphology of hexagonal-CoSe is broken and islands of tetragonal-CoSe are formed. Our results provide nanoscopic insights into the transformation processes of 2D CoSe2 which can be used to generate these intriguing 2D materials and to tune their properties by modifying their structures for electro-catalytic and electronic applications.
Similar content being viewed by others
Introduction
Two-dimensional transition metal dichalcogenides (2D TMDs) have attracted intense interest recently due to their layer dependent physical and chemical properties and their potential applications in various fields1,2,3,4. Many TMDs can adopt different crystal structures exhibiting diverse electronic and magnetic properties. Tuning the structures of these TMDs through various modes is an interesting topic in 2D materials5,6,7. As an example, single-layered MoS2 has two different symmetries depending on the position of sulfur atoms, semiconducting trigonal prismatic (2H) and metallic octahedral (1T). The structural transition between these two phases has been observed through in situ microscopy8. Various other materials can be interconverted by numerous stimuli, such as substitutional doping9,10, applying pressure11, electron irradiation8,12, etc.
Layered metal chalcogenides show several stable phases with different stoichiometry. Layered CoSe2 also possesses multiple stable phases with distinct properties and stoichiometries. Orthorhombic marcasite (o-CoSe2) is one of these structures, and is catalytically active for the hydrogen evolution reaction (HER) in the context of water splitting13,14. CoSe2 can also be found in the cubic phase (c-CoSe2) which belongs to the pyrite-type structures with Se22− dumbbells and with Co2+ in an octahedral coordination15. Recently, c-CoSe2 has been investigated as a promising candidate for HER as well14,16,17,18. The transformation from orthorhombic to cubic in CoSe2 has been reported also by other studies15,19,20. Ex situ heating experiments confirmed the transformation of orthorhombic to the cubic structure of CoSe221. Cobalt selenide in the different stoichiometry of monochalcogenide CoSe also exhibits two different structures with distinct properties: the hexagonal phase and the tetragonal phase22,23. The hexagonal form exhibits a non-layered structure23 while the tetragonal form exhibits a layered structure24,25. Hexagonal-CoSe nanoparticles are explored as promising catalysts in various fields26,27.
In this work, we studied the thermal evolution by in situ heating in a transmission electron microscope (TEM), which permits real-time imaging of transforming nanostructures at atomic resolution28. Our findings provide understanding of the structural transformation mechanisms in 2D CoSe2 and 2D CoSe. This is useful for tuning the structure–property relationships of the material by controlling process conditions, in the present case, by means of heating.
Results and discussion
Exfoliation and characterization of CoSe2 nanosheets
Figure 1 shows plan-view TEM and high-resolution TEM (HR-TEM) images of few-layers-thick o-CoSe2 nanoflakes (Fig. 1a, b), exfoliated, and deposited onto a SiN support membrane. The HR-TEM image in Fig. 1b displays the high crystalline quality of o-CoSe2 in thin areas and confirms the orthorhombic symmetry of Co atomic columns in [010] direction. The corresponding fast Fourier transform (FFT) pattern of the HR-TEM image shown in the inset of Fig. 1b indicates the orthorhombic crystal structure of CoSe2 (see Supplementary Fig. 1 for more structural details). The filtered HR-TEM image and atomic model of o-CoSe2 (Fig. 1c, d) confirm the orthorhombic configuration. Energy-dispersive X-ray spectrometry (EDS) elemental mapping images (Fig. 1g, h) confirm the homogeneous distribution of Co (indicated by red) and Se (indicated by green) in the flakes. The EDS spectrum (Supplementary Fig. 2) reveals that Co and Se are the only the elements in the prepared sample with a 1:2 ratio for Co/Se. Atomic force microscopy (AFM) imaging of a monolayer of o-CoSe2 (Fig. 1j) confirms the single unit cell thickness of ~0.6 nm. Vertically oriented o-CoSe2 nanoflakes on carbon fiber were used to exfoliate ultrathin nanosheets as shown in the scanning electron microscopy image in Fig. 1k. Few-layered (4–5 layers) thick o-CoSe2 nanosheets were used as starting materials to investigate thermally induced structural transformations. Experimental details are provided in the “Methods” section.
Heating-induced transformations monitored by in situ TEM
We performed many in situ experiments, which all showed a similar sequence of structural transformations. Here we provide a detailed overview of the successive transformations and show several highlights of transformation events, while a selection of additional in situ observations and corresponding analysis are shown in Supplementary Figs. 3–19 of the Supplementary Information (SI). The transformations took place everywhere, also in areas not previously exposed to the electron beam, as discussed also in the Supplementary Text in the SI.
The first transformation is the orthorhombic to cubic transformation of CoSe2. Figure 2 shows a single area where the thermal evolution could be monitored over multiple transformations. Crystallographic details and lattice parameters of the CoSe2 and CoSe phases can be found in Supplementary Tables 2 and 3 of the SI. The TEM image of Fig. 2a shows the starting point of o-CoSe2, in a [100] zone axis (ZA) projection. As is also clear from the FFT (Fig. 2a′), the structure is orthorhombic and this was the best orientation to investigate any heating-induced structural changes. The o-CoSe2 nanosheets were heated in the TEM from room temperature to 1100 °C in steps of 20 °C. The orthorhombic structure of CoSe2 remains stable during the increase in temperature up to 580 °C without showing any significant structural changes. This can be observed from the HR-TEM image in Fig. 2a which was recorded at 580 °C. The corresponding FFT pattern is shown in Fig. 2a′ and the spots indicate the crystal structure of the orthorhombic phase projected along the [100] axis. Monitoring the change to the cubic phase is not trivial as during this change the stoichiometry and approximate mass density of the material is retained, so that the phases can only be distinguished by inspection of atomic lattice fringes. Movies capturing various stages of the evolution were recorded while increasing the temperature. The transition from the orthorhombic to cubic phase goes through one intermediate stage where both phases can be distinguished in Fig. 2b and in the FFT pattern shown in Fig. 2b′. This transition stage suggests a rearrangement of Co and Se lattice atoms while retaining the CoSe2 stoichiometry. The metal dichalcogenides exist in orthorhombic marcasite-type or cubic pyrite-type structures in which the metal atoms are octahedrally bonded to adjacent chalcogen atoms. The structural difference is mainly related to the octahedral linkage as the octahedra are edge-shared in marcasite and corner-shared in pyrite20,29,30. The transformations from marcasite to pyrite structures involve rotation of metal-chalcogen chains in marcasite that produces disordered pyrite structures20. This explains how the orthorhombic (marcasite type) structure converts to cubic (pyrite type) while retaining the same stoichiometry. The o-CoSe2 and c-CoSe2 phases have a similarity in structure with a small lattice mismatch, and there are reports of easy epitaxial growth between them30,31.
The full conversion of the orthorhombic phase to the cubic structure was observed after increasing the temperature to 740 °C. At this temperature the entire nanocrystal of o-CoSe2 had completely converted into cubic-CoSe2, as can be seen from the TEM image in Fig. 2c. The corresponding FFT pattern in Fig. 2c′ confirms the cubic structure where a few new spots can be seen that were absent for o-CoSe2. Increasing the temperature further makes c-CoSe2 unstable and undergo a subsequent structural transition ending in the hexagonal-CoSe phase (h-CoSe) at 860 °C. This second structural transformation took place via an intermediate transition phase where the presence of both cubic and hexagonal domains of CoSe2 and CoSe were seen (Fig. 2d–f). Details on both of these transition processes are discussed in the following sections.
The first transition, from orthorhombic to cubic-CoSe2, involves an intermediate phase during which signatures of both phases are found. Supplementary Figure 3a and a′ in the Supplementary Information shows an o-CoSe2 structure with (020), (004), and (024) lattice fringes at 660 °C. This suggests thermal stability of the orthorhombic phase at this high temperature, but it lasts only for a short period of time, as at this temperature the transition process starts. Supplementary Fig. 3a–d shows stills of the movie at 0, 100, 150, and 200 s, respectively. The appearance of the intermediate transition phase can be seen from the FFT in Supplementary Fig. 3d′, where new spots in the FFT appeared and correspond to the cubic phase. The transitionary phase occured in the temperature window of 660–740 °C, after which a pure cubic phase was observed at 740 °C, as shown in Supplementary Fig. 3e,f. We inspected other areas on the heating chip to confirm this transition at 740 °C and found that the transformation process is uniform over the entire region on the heating chip. This confirms a very low or negligible effect of the electron beam on the process at these high temperatures. To investigate the change in shape and size of the nanosheet under observation, we plotted time-labeled contours of the nanosheet area (Supplementary Fig. 3g), which indicates a minor change in the shape and size of that area. The plot shows a change in area of ~1 nm2, which is negligible and suggests no loss of either Co or Se atoms. To investigate the stoichiometry of the transformed cubic phase we employed EDS elemental mapping as shown in Supplementary Fig. 4 and found the results to be in agreement with a 1:2 ratio for Co/Se.
As shown in Supplementary Video 1, we observed this structural transformation at identical temperature conditions (740 °C) on another TEM microscope as well, where a lower acceleration voltage of 200 kV was used (see “Methods” for details). It shows a smooth transition from the orthorhombic to the cubic phase in a slightly thicker, lath-shaped nanosheet of CoSe2. Supplementary Figure 5 shows stills from the recorded movie (Supplementary Video 1) of the transformation with initial, midway, and final structures. The initial structure is orthorhombic oriented along the [010] ZA with (200) and (002) lattice fringes (Supplementary Fig. 6a′). The transformed cubic structure oriented along \(\left[ {\bar 111} \right]\) ZA is shown in Supplementary Fig. 6b with corresponding FFT pattern in Supplementary Fig. 6b′. According to mechanical alloying theory32, thermal displacement of Co and/or Se atoms to the nearest neighboring sites produces disorder and instability in the orthorhombic structure. This instability results in the partial conversion to a cubic structure and the appearance of co-existence of both structures, followed by a further diminishing of the orthorhombic structure and the transformation of all crystals to the cubic structure.
The second transition from c-CoSe2 to h-CoSe also involves one intermediate phase that can be seen in Supplementary Fig. 7. This transition process was observed in the temperature range of 740–860 °C, where the cubic structure continuously went through small changes. The transition process starts at 760 °C at which the cubic structure starts to convert into an unknown intermediate phase. This pathway is likely initiated by the removal of Se atoms from the CoSe2 lattice after increasing the temperature to 760 °C. For other TMDs it was previously observed that chalcogen atoms can be removed due to thermal and electron irradiation effects33. This removal of Se atoms changes the c-CoSe2 configuration with a Co:Se stoichiometry of 1:2, to CoSe with stoichiometry of 1:1. EDS elemental mapping of the transformed structure shown in Supplementary Fig. 8 shows a 1:1 ratio of Co and Se in the h-CoSe structure. Local thinning of the area during the transition process can be observed from the TEM images in Supplementary Fig. 7a–f, and could be the result of the removal of Se atoms.
Changes in shape and reduction in size is another noticeable modification during this transition. This can be seen from the TEM images in Supplementary Fig. 7a–f. Time-labeled contours of this area are shown in Supplementary Fig. 9g. The change during this transformation is much more pronounced as the decrease in area against time (Supplementary Fig. 9h) is nearly 10 times greater than in the case of the orthorhombic to cubic transition. This could be the result of thinning of the flakes and change of shape due to the removal of Se atoms along with the rearrangement of Co and Se atoms into a non-layered hexagonal-CoSe lattice structure. At 860 °C, the h-CoSe phase has a small area compared to the initial area of the cubic phase of CoSe2 at 740 °C, and becomes thinner as well, as can be seen in Supplementary Fig. 7a, f. As shown in Supplementary Fig. 9a–f, the intermediate phase starts to convert and stabilizes to h-CoSe at a temperature of 800 °C. The time series of TEM images shows a long process of 24 min during which the temperature was kept at 800 °C. The detailed time series is shown in Supplementary Figs. 11 and 12. At 840 °C h-CoSe (Supplementary Fig. 10) undergoes a relaxation process, which ends with the stable h-CoSe structure. We have observed movement and rearrangement of Co atomic columns during this transition process, as can also be seen from the movie recordings (Supplementary Video 2). According to mechanical alloying theory32, an increase in the concentration of defects in the cubic structure (in our case, removal of chalcogen atoms (Se atoms) from the unit cell thick nanosheets) increases the Gibbs free energy, which favors movement of Co and/or Se atoms to the nearest neighboring sites and increases the instability. Further increase in the concentration of defects leads to a critical value of the Gibbs free energy, after which the cubic structure collapses and transforms to the stable hexagonal-CoSe structure32. Further increasing the temperature (above 900 °C) results in another transition process in which the hexagonal structure of CoSe transforms into a tetragonal structure. The layers of hexagonal-CoSe break into islands of tetragonal-CoSe. During this transition, hexagonal layers were found to etch away slowly and form islands of tetragonal-CoSe; these islands were found to be stable at 900–1000 °C.
The detailed analysis of the structures before and after transformation, from orthorhombic to cubic-CoSe2, to hexagonal-CoSe, and to tetragonal-CoSe, is shown in Fig. 3. We started with the orthorhombic structure of CoSe2, and the analysis of its structural details is described in Fig. 3a–d. A part of a nanoflake with an orthorhombic structure is shown in Fig. 3a. This nanocrystal undergoes a structural transformation upon heating and obtains a stable cubic structure at the respective temperature (740 °C). Figure 3e shows the HR-TEM image of a transformed structure and agrees well with the cubic structure for CoSe2. The experimental image shows a square pattern as the contrast from the Co atoms (that are in a square pattern in this projection, see the atomic model shown in Fig. 3h) is much stronger than the contrast from the Se atoms, as was also confirmed by HR-TEM imaging simulations provided in Supplementary Fig. 19 of the SI. The FFT analysis in Fig. 3f shows (220), \(\left( {2\bar 20} \right)\), and (400) atomic lattice fringes, corresponding to interplanar distances of 0.207 and 0.146 nm, respectively, and indexed along the [001] ZA. As the correct identification of the transformation from orthorhombic to cubic is critical, measurements of interplanar distances played an important role in distinguishing the structures. In the next transition process, cubic-CoSe2 transforms into hexagonal-CoSe. As CoSe has two different structures, non-layered hexagonal and layered tetragonal, it is of importance to verify the structure of the transformation product. The features of the formed CoSe agree well with the hexagonal structure as shown in Fig. 3i, where the basal-plane oriented structure can be seen. The FFT analysis (Fig. 3j) of this structure clearly indicates that the cubic-CoSe2 completely transformed into hexagonal-CoSe. The filtered image (Fig. 3k) shows the hexagonal nature of the flake which matches well with the schematic atomic model (Fig. 3l) for h-CoSe (CIF IDs for all the structures are given in Supplementary Table 3 in SI). The hexagonal-CoSe structure has a hexagonal unit cell with P63/mmc space group where each Co atom is surrounded by six Se atoms (Fig. 3l). The unmasked FFTs for all the structures along with the HR-TEM images are shown in Supplementary Fig. 15.
After increasing the temperature above 900 °C, yet another phase, tetragonal-CoSe, was formed. This transformation was observed to take place rapidly where the nanosheets of hexagonal-CoSe break up into islands of tetragonal-CoSe (Fig. 3m) no longer displaying a layered structure. Analyis of the FFT pattern (Fig. 3n) confirms the aforementioned observations with atomic lattice fringes (110), \(\left( {1\bar 10} \right)\), and (020) with the crystal oriented along the [001] ZA. It can be seen from the Fourier-filtered image (Fig. 3o) and the schematic atomic model (Fig. 3p) that the observed structure is tetragonal-CoSe. Tetragonal-CoSe has a tetragonal unit cell in P4/nmm space group in which each Co atom is surrounded by four Se atoms (Fig. 3p). Structural details of the four phases can be found in Supplementary Tables 2 and 3.
Among the two different structures of CoSe, the hexagonal structure is of great scientific importance because of its applications in particular in catalysis26,27 and it is of interest to investigate its observed structure through lattice resolved TEM images. Figure 4 shows a detailed description of hexagonal-CoSe with the help of HR-TEM images, interplanar distance measurements, and a schematic atomic model. Figure 4a shows a layered structure of the basal-plane-oriented hexagonal-CoSe. Though hexagonal-CoSe has a non-layered structure23, the transformation process results in a layered structure with a unit cell thickness of hexagonal-CoSe. The atomic model in Fig. 4c mimics the filtered image in Fig. 4b. It shows how the hexagonal structure of CoSe consists of Co and Se atoms forming a hexagon. The interplanar distances found in the HR-TEM image (Fig. 4a, b) of 0.35 and 0.18 nm agree very well with the interplanar distances of the (100) and (110) planes, respectively, of the hexagonal-CoSe structure. Vertically oriented layers of hexagonal-CoSe are shown and described in Supplementary Fig. 13 (Supplementary Fig. 14 shows more details on t-CoSe). While mainly observing basal-plane-oriented layers of hexagonal-CoSe, some of the edges rolled up and can be seen as vertically oriented layers. This helped to perform interlayer measurements, where the layers were found to be separated by approximately 0.5 nm from each other (Supplementary Fig. 13d, e). The ball and stick model used to mimick the HR-TEM image of edges through [010] ZA agrees with the observed results and confirms the layered structure of hexagonal-CoSe.
Figure 5 shows yet another illustrative example of subsequent transformations that were monitored in situ in the same area of the same heating chip, confirming the generally observed transformation sequence. Furthermore, the structural transformations are observed at the different windows of the heating chip in a single experiment which suggests that these transitions occurred everywhere and not only at those areas previously illuminated by the electron beam (see Supplementary Figs. 16–18 for more details). TEM imaging simulations were carried out based on the structures refined by means of the density functional theory (DFT) calculations, and are shown in Supplementary Fig. 19, agreeing well with the experimental results.
Density functional theory calculations
In order to gain insight into the energetics of the transformations, DFT calculations were conducted using the plane-wave Vienna Ab initio Simulation Package34,35 employing the projector augmented wave approach36,37 and the generalized gradient approximation (GGA) with the exchange-correlation functional by Perdew, Burke and Ernzerhof (PBE)38. The structures were also evaluated using a van der Waals density functional (optB88-vdW)39,40. Computational details are provided in the “Methods” section.
Figure 6a provides an overview of the crystal structures of all phases involved in the transformations. Transformation (1) where the orthorhombic-CoSe2 phase transforms into the cubic-CoSe2 phase was already described above, and can be explained from rotations of metal-chalcogen chains20. The subsequent transformation (2), from c-CoSe2 to h-CoSe, and transformation (3), from h-CoSe to t-CoSe, are experimentally observed in the present work. The intermediate phase that was observed during transformation (2) may correspond to the Co3Se4 structure41 and we have included this phase in the computations. The sequence ends with the tetragonal t-CoSe phase which clearly consists of 2D atomic planes that are stacked along the c-axis. The experimental value of its c lattice parameter is 5.25 Å (Supplementary Table 2). This value was strongly overestimated when using the standard GGA-PBE functional because of the neglect of van der Waals interactions, while better agreement was found when using the optB88-vdW functional (Supplementary Table 4). For this reason, we show in Fig. 6b the relative stability of all phases calculated using opB88-vdW. The formation energies of all phases are listed in Supplementary Table 4 and the DFT-optimized crystal structures are provided as Supplementary Data as well.
The formation energies displayed in Fig. 6b were calculated as \(E^{\mathrm {form}}\left( {{\mathrm {Co}}_x{\mathrm {Se}}_y} \right) = E\left( {{\mathrm {Co}}_{x}{\mathrm {Se}}_y} \right) - {{x}}\,E\left( {{\mathrm {Co}}} \right) - {{y}}\,E\left( {{\mathrm {Se}}} \right),\)where all energies are expressed in eV/atom; E(CoxSey) is the energy of the particular compound, and E(Co) and E(Se) are the energies of the elemental phases of Co and Se, respectively. All formation energies displayed in Fig. 6b are thus defined with respect to the elemental phases, and are plotted as a function of their composition. The black line is the common tangent line connecting the most favorable phases. It is clear that the transformation sequence displayed in Fig. 6a is feasible from the point of view of the total energies: The o-CoSe2 and c-CoSe2 formation energies are on or close to the tangent line and differ only slightly in magnitude, implying that they are of similar thermodynamic stability. Both c-CoSe2 and hexagonal-CoSe are quite close to the tangent line, rendering them both relatively stable at their respective compositions. The monoclinic Co3Se4 phase is compositionally and energetically positioned in between c-CoSe2 and h-CoSe, and is also close to the tangent line. The Co3Se4 phase exhibits both cubic and hexagonal structural features as visualized in Supplementary Fig. 20. Finally, the last transformation from h-CoSe to tetragonal t-CoSe is clearly energetically unfavorable according to the phase stability diagram of Fig. 6b. It is very likely, however, that the planar-like structure of t-CoSe can well accommodate phonons, rendering it a more stable phase at elevated temperatures (the DFT calculations are valid only for a temperature of 0 K and phonon energies are not included).
Figure 6c shows the strong similarly between the \(\left( {\bar 1100} \right)\) atomic planes of the h-CoSe phase and the (001) atomic planes of the t-CoSe phase. Exactly these planes are experimentally observed on-edge at the left-hand side of Fig. 4a, where \(\left( {\bar 1100} \right)\) layers are peeling off from the hexagonal-CoSe to form the (001) layers of tetragonal-CoSe. As is clear from the perspective views of these {CoSe} layers at the bottom of panel c in Figure 6, the Co sublattice in both layers is simple square (exactly square in the c-plane of t-CoSe, and approximately square in the \(\left( {\bar 1100} \right)\) plane of h-CoSe). The main difference is that in the h-CoSe structure, the up/down alternation of the Se atoms is only in one direction (along the c-axis), while the up/down alternation of the Se atoms in the t-CoSe structure is in both directions. However, a change in buckling of the Se atoms around the Co plane can very likely be induced by thermal excitations, resulting in these layers to be peeled off from the h-CoSe phase as observed in Fig. 4a, and leading to the formation of the t-CoSe phase.
In conclusion, we used in situ electron microscopy to investigate the heat-induced structural transformation in 2D CoSe2 nanosheets. Transformations from the orthorhombic to the cubic structure of CoSe2, from cubic-CoSe2 to hexagonal-CoSe, and from hexagonal-CoSe to tetragonal-CoSe, were observed in real time in the TEM. Different pathways for distinct transformations produced layered structures of cubic-CoSe2 and hexagonal-CoSe. The initial transformation process to cubic-CoSe2 is driven by heat-induced transitions with rearrangement of atoms in the lattice. The subsequent transformation is driven by a loss of chalcogen atom due to heating at elevated temperatures, and leads to a mix of di- and monochalcogenides, followed by the complete transformation to hexagonal-CoSe. By using and intentionally applying these heating-induced transformations to orthorhombic-CoSe2 nanosheets, high-quality ultrathin layers of the other phases (cubic-CoSe2, hexagonal-CoSe, and tetragonal-CoSe) can be prepared. Furthermore, the formation of well-defined CoSe2/CoSe heterostructures may be explored by precise tuning of the processing conditions. The transformation mechanisms elucidated here help understanding and controlling the conversion of layered cobalt chalcogenides and other layered materials exhibiting different stable structures with various stoichiometries, where transformations are driven by thermal effects and removal of chalcogen atoms.
Methods
Synthesis of o-CoSe2 nanosheets
The o-CoSe2 nanosheets were prepared by the hydrothermal method. Lamellar hybrid CoSe2/DETA (diethylenetriamine) intermediates were synthesized as described in the existing literature42. As prepared samples were further exfoliated into fairly thin nanosheets using a probe sonicator. In this procedure, 10 mg of CoSe2 was dispersed in 20 ml of ethanol and then ultrasonicated for 6 h. After ultrasonic treatment, the resultant dispersions were centrifuged at 3000 r.p.m. for 10 min to remove the unexfoliated component. This procedure yielded very thin nanosheets of CoSe2 as shown in Fig. 1. If the ultrasonication was performed for a shorter period of time, slightly thicker nanosheets such as the lath-shaped nanosheets shown in Supplementary Fig. 5 were found in the dispersion as well.
TEM details and settings
The in situ TEM study of orthorhombic-CoSe2 nanosheets converting to hexagonal-CoSe nanosheets was performed on a Cs-corrected TEM (FEI Titan 80-300) operating at an acceleration voltage of 300 kV. Where mentioned, experiments were conducted on a FEI Talos F200X TEM operating at 200 kV, as well. All EDS chemical mapping experiments were performed on the Talos F200X TEM equipped with a Chemi-STEM elemental analysis setup. CoSe2 samples for in situ heating experiments were prepared by drop casting it directly onto MEMS chips with SiN membranes, which were mounted on a DENSsolutions heating holder. The temperature was increased from room temperature to 1100 °C in steps of 20 °C. A Direct Electron camera on the FEI Titan TEM was used to record movies (frame rate 25 frames/s) and images for observing dynamical structural transformations in HR-TEM mode. The total screen current (fluscreen) was typically 5 nA and the beam exposure on the DE camera was limited to 4000 electrons per frame. Scanning TEM (STEM) imaging on the Talos F200X was conducted using a probe current of 30 pA and a dwell time per pixel of 4.0 μs. For all the filtered images in this article, the contrast was improved between the material and amorphous region in the background by applying a mask on the amorphous region followed by the inverse fast Fourier transformation (IFFT) for better display purposes. The TEM imaging simulations (Supplementary Fig. 19) were performed using QSTEM. The simulated images were generated with the following settings: accelerating voltage: 200 kV, objective aperture: 15 mrad, convergence angle: 1 mrad, focal spread: 2 nm, defocus: 71, 10, 26, and 20 nm for orthorhombic, cubic, hexagonal, and tetragonal structures, respectively.
Atomic force microscopy
AFM measurements were performed on a JPK Nanowizard II with a Nikon Eclipse Ti-U inverted microscope. The machine was placed on top of an active vibration isolation table together with a Halcyonics Active Vibration Isolation Unit. The whole setup was placed inside a large acoustic vibration box to further isolate the setup from external vibrations. We have used Bruker OTESPA-R3 AFM tips.
DFT calculations
The energy cutoff and the density of the k-meshes were tested to ensure that the energies converged to well within 0.5 meV/atom. The cutoff energies of the wave functions and augmentation functions were 500 and 700 eV, respectively, and the wave functions were sampled in k-space using the method by Monkhorst and Pack. The cubic c-CoSe2 structure was sampled on a 12 × 12 × 12 k-mesh corresponding to a linear reciprocal spacing of 0.014 Å−1. All other structures were evaluated with k-meshes that were at least equally dense. Both non-magnetic and magnetic (ferromagnetic and antiferromagnetic) input configurations were considered. All magnetic solutions that were found, were not energetically more favorable than the non-magnetic solutions within the accuracy of the method, and here we report only the non-magnetic solutions. Full relaxation with respect to cell dimensions and atomic coordinates was applied to find lower-energy configurations. Zero-point vibration contributions are neglected and the calculations are valid for a temperature of 0 K and a pressure of 0 Pa. The optB88-vdW optimized structure files (POSCAR files) are included as Supplementary Data.
Data availability
Most of the data generated and analyzed during this study is included in this Article and in the Supplementary Information file. DFT-calculated structure files are provided as Supplementary Data. Raw experimental data and calculation data will be stored in the Data Repository of Utrecht University, and are available from the authors upon request.
References
Cheng, L. et al. Ultrathin WS2 nanoflakes as a high‐performance electrocatalyst for the hydrogen evolution reaction. Angew. Chem. Int. Ed. 53, 7860–7863 (2014).
Cheng, R. et al. Few-layer molybdenum disulfide transistors and circuits for high-speed flexible electronics. Nat. Commun. 5, 5143 (2014).
Steinleitner, P. et al. Direct observation of ultrafast exciton formation in a monolayer of WSe2. Nano Lett. 17, 1455–1460 (2017).
Zhou, H. et al. Large area growth and electrical properties of p-type WSe2 atomic layers. Nano Lett. 15, 709–713 (2014).
Chhowalla, M. et al. The chemistry of two-dimensional layered transition metal dichalcogenide nanosheets. Nat. Chem. 5, 263 (2013).
Empante, T. A. et al. Chemical vapor deposition growth of few-layer MoTe2 in the 2H, 1T′, and 1T phases: tunable properties of MoTe2 films. ACS Nano 11, 900–905 (2017).
Zhou, Y. et al. Pressure-induced Td to 1T′ structural phase transition in WTe2. AIP Adv. 6, 075008 (2016).
Lin, Y.-C., Dumcenco, D. O., Huang, Y.-S. & Suenaga, K. Atomic mechanism of the semiconducting-to-metallic phase transition in single-layered MoS2. Nat. Nanotechnol. 9, 391 (2014).
Enyashin, A. N. et al. New route for stabilization of 1T-WS2 and MoS2 phases. J. Phys. Chem. C 115, 24586–24591 (2011).
Kappera, R. et al. Phase-engineered low-resistance contacts for ultrathin MoS2 transistors. Nat. Mater. 13, 1128 (2014).
Nayak, A. P. et al. Pressure-induced semiconducting to metallic transition in multilayered molybdenum disulphide. Nat. Commun. 5, 3731 (2014).
Wang, L., Xu, Z., Wang, W. & Bai, X. Atomic mechanism of dynamic electrochemical lithiation processes of MoS2 nanosheets. J. Am. Chem. Soc. 136, 6693–6697 (2014).
Basu, M. et al. Heterostructure of Si and CoSe2: a promising photocathode based on a non‐noble metal catalyst for photoelectrochemical hydrogen evolution. Angew. Chem. Int. Ed. 54, 6211–6216 (2015).
Chen, P. et al. Phase‐transformation engineering in cobalt diselenide realizing enhanced catalytic activity for hydrogen evolution in an alkaline medium. Adv. Mater. 28, 7527–7532 (2016).
Brostigen, G. & Kjekshus, A. Relationships between the structure types pyrite, marcasite, and arsenopyrite. Acta Chem. Scand. 24, 2983–2992 (1970).
Gao, M.-R. et al. An efficient molybdenum disulfide/cobalt diselenide hybrid catalyst for electrochemical hydrogen generation. Nat. Commun. 6, 5982 (2015).
Xu, Y. F., Gao, M. R., Zheng, Y. R., Jiang, J. & Yu, S. H. Nickel/nickel (II) oxide nanoparticles anchored onto cobalt (IV) diselenide nanobelts for the electrochemical production of hydrogen. Angew. Chem. Int. Ed. 52, 8546–8550 (2013).
Zheng, Y.-R. et al. Cobalt diselenide nanobelts grafted on carbon fiber felt: an efficient and robust 3D cathode for hydrogen production. Chem. Sci. 6, 4594–4598 (2015).
Dodony, I., Posfal, M. & Buseck, P. R. Structural relationship between pyrite and marcasite. Am. Mineral. 81, 119–125 (1996).
Fleet, M. Structural aspects of the marcasite-pyrite transformation. Canad. Mineral. 10, 225–231 (1970).
Zhang, H. et al. Polymorphic CoSe2 with mixed orthorhombic and cubic phases for highly efficient hydrogen evolution reaction. ACS Appl. Mater. Int. 7, 1772–1779 (2015).
Campos, C. et al. Hexagonal CoSe formation in mechanical alloyed Co75Se25 mixture. Sol. Stat. Commun. 131, 265–270 (2004).
Umeyama, N. et al. Synthesis and magnetic properties of NiSe, NiTe, CoSe, and CoTe. Jap. J. Appl. Phys. 51, 053001 (2012).
Wilfong, B. et al. Frustrated magnetism in the tetragonal CoSe analog of superconducting FeSe. Phys. Rev. B 97, 104408 (2018).
Zhou, X. et al. Metastable layered cobalt chalcogenides from topochemical deintercalation. J. Am. Chem. Soc. 138, 16432–16442 (2016).
Nekooi, P., Akbari, M. & Amini, M. K. CoSe nanoparticles prepared by the microwave-assisted polyol method as an alcohol and formic acid tolerant oxygen reduction catalyst. Int. J. Hydrogen Energy 35, 6392–6398 (2010).
Zhang, Y. et al. Nitrogen-doped yolk–shell-structured CoSe/C dodecahedra for high-performance sodium ion batteries. ACS Appl. Mater. Int. 9, 3624–3633 (2017).
van Huis, M. A. et al. Atomic imaging of phase transitions and morphology transformations in nanocrystals. Adv. Mater. 21, 4992–4995 (2009).
Cabán-Acevedo, M., Faber, M. S., Tan, Y., Hamers, R. J. & Jin, S. Synthesis and properties of semiconducting iron pyrite (FeS2) nanowires. Nano Lett. 12, 1977–1982 (2012).
Sun, R., Chan, M. & Ceder, G. First-principles electronic structure and relative stability of pyrite and marcasite: Implications for photovoltaic performance. Phys. Rev. B 83, 235311 (2011).
Kong, D., Wang, H., Lu, Z. & Cui, Y. CoSe2 nanoparticles grown on carbon fiber paper: an efficient and stable electrocatalyst for hydrogen evolution reaction. J. Am. Chem. Soc. 136, 4897–4900 (2014).
Campos, C. et al. Structural studies of cobalt selenides prepared by mechanical alloying. Phys. B Condens. Matter 324, 409–418 (2002).
Sutter, E. et al. Electron-beam induced transformations of layered tin dichalcogenides. Nano Lett. 16, 4410–4416 (2016).
Kresse, G. & Furthmüller, J. Efficiency of ab-initio total energy calculations for metals and semiconductors using a plane-wave basis set. Comput. Mater. Sci. 6, 15–50 (1996).
Kresse, G. & Furthmüller, J. Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set. Phys. Rev. B 54, 11169 (1996).
Blöchl, P. E. Projector augmented-wave method. Phys. Rev. B 50, 17953 (1994).
Kresse, G. & Joubert, D. From ultrasoft pseudopotentials to the projector augmented-wave method. Phys. Rev. B 59, 1758 (1999).
Perdew, J. P., Burke, K. & Ernzerhof, M. Generalized gradient approximation made simple. Phys. Rev. Lett. 77, 3865 (1996).
Dion, M., Rydberg, H., Schröder, E., Langreth, D. C. & Lundqvist, B. I. Van der Waals density functional for general geometries. Phys. Rev. Lett. 92, 246401 (2004).
Klimeš, J., Bowler, D. R. & Michaelides, A. Chemical accuracy for the van der Waals density functional. J. Phys. Condens. Matter 22, 022201 (2009).
Garcı́a-Garcı́a, F. J., Larsson, A.-K., Norèn, L. & Withers, R. L. The crystal structures of Co3Se4 and Co7Se8. Sol. Stat. Sci. 6, 725–733 (2004).
Liang, L. et al. Metallic single‐unit‐cell orthorhombic cobalt diselenide atomic layers: robust water‐electrolysis catalysts. Angew. Chem. Int. Ed. 54, 12004–12008 (2015).
Acknowledgements
This project was financially supported by the European Research Council through an ERC Consolidator Grant NANO-INSITU (No. 683076).
Author information
Authors and Affiliations
Contributions
D.S.G. conceived the study and performed the in situ TEM experiments. B.T. and G.L. prepared the samples and M.A.M. supervised the work. C.P. conducted the AFM measurements and D.S.G. analyzed the data. H.v.G. performed and analyzed the DFT calculations. D.G. and M.A.v.H. wrote the manuscript. All authors commented on the manuscript and approved the final version.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary information
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Gavhane, D.S., van Gog, H., Thombare, B. et al. In situ electron microscopy study of structural transformations in 2D CoSe2. npj 2D Mater Appl 5, 24 (2021). https://doi.org/10.1038/s41699-021-00206-3
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41699-021-00206-3
This article is cited by
-
Time-resolved transmission electron microscopy for nanoscale chemical dynamics
Nature Reviews Chemistry (2023)