Defect-assisted reversible phase transition in mono- and few-layer ReS2

Defect-assisted reversible phase transition in mono- and few-layer ReS2

Introduction

Over the past two decades, following the seminal work on graphene1, there has been a substantial growth of interest in atomically thin 2D materials. Among these, transition metal dichalcogenides (TMDs) have emerged as a particularly intriguing family due to their remarkable performance in various nanophotonic2,3, nanoelectronic4, and other nanotechnology applications5. One highly appealing aspect of TMDs is their versatile nature, which allows for their tuning during device fabrication via various control parameters.

TMDs form a diverse material platform comprising transition metal (e.g., Mo, W, Re, Nb, Ta, etc.) and chalcogenide elements (e.g., S, Se, Te). The flexibility to choose these constituent elements during synthesis enables the tunability of TMDs’ optical6 and electronic properties7. Furthermore, stacking atomically thin 2D materials in van der Waals heterostructures8 offers unprecedented control over their properties and can lead to unconventional quantum phenomena, such as superconductivity in graphene-based structures9,10, Moiré excitons and correlated states in TMDs11,12,13,14,15,16, topological excitons17, advanced excitonic devices18, and exciton polaritons19,20,21. Following the fabrication of large TMD flakes, achieved through methods like exfoliation22 or epitaxial single crystal growth23,24, further modifications of material properties can be realized through nanopatterning, e.g., via electron beam lithography (EBL) and chemical etching. Nanopatterning TMDs has enabled the production of high-quality factor (high-Q) nanophotonic structures25,26,27,28, exploration of higher order modes and anapole states29,30,31,32, and lasing from indirect bandgap semiconductor WS2 microdisks33. Expanding control over TMD properties beyond conventional fabrication methods will enable more precise engineering of TMD device performance and deeper insights into the fundamental physics of TMD-related phenomena34,35.

The ability to actively and reversibly tune TMD properties post-fabrication is of high interest within the 2D material community36. Realizing this control to its full potential could expand the range of applications for TMDs, including neuromorphic computing37,38 as well as optically rewritable electronic circuits and memories39. One approach for achieving active control over TMD electronic properties involves photo-induced modifications of the concentration of substitutional and interstitial point defects40. Another approach leverages the substantial changes in properties associated with material polymorphism or phase transitions. For instance, entire TMD flakes can undergo phase transitions driven by electrostatic doping41 or mechanical strain42. While these approaches show promise, they typically lack temporal stability and/or precise spatial control at the nanometer scale. Laser-induced phase patterning offers more precise spatial control but is often irreversible due to ablation43,44, similar to thermal heating during annealing in 2D materials45. Nevertheless, promising routes for photo-induced phase switching exist46,47, such as the reversible transition of TaS2 into a hidden stable quantum state, enabling drastic changes in conductivity48,49. Therefore, exploring new mechanisms for phase transitions in TMDs and understanding their associated property changes are crucial for achieving improved active control.

Rhenium disulfide (ReS2) is a promising TMD material for optical, electronic, and gas-sensing applications50,51,52,53,54,55. Its unique features include being a direct bandgap semiconductor in a bulk van der Waals crystal56 and its high optical anisotropy, making it ideal for directionally selective photonic applications6,57,58,59,60. Previous investigations have focused on ReS2’s atomic structure and phases using density functional theory (DFT)56 and scanning transmission electron microscopy (STEM)61. Additionally, second-harmonic generation (SHG) studies have linked its activity to stacking and thickness variations62. Notably, ReS2 monolayers (ML) and bilayers (BL) have been demonstrated to exhibit a reversible phase transition under near-infrared femtosecond (fs) laser irradiation, characterized by a significant change in SHG activity and polarization-resolved SHG63. It was hypothesized from the available data that ReS2 transitions from the (1{rm{T}}^{prime}) phase into the 1H phase63. Despite this promising initial work, questions remain regarding the involved transition phases and the mechanism enabling the transition.

Here, we investigate the phase transition in ReS2 using a combination of electronic structure calculations, optical SHG microscopy, STEM imaging and diffraction, and electronic transport measurements. Our study results in several key observations: (i) femtosecond laser irradiation induces a phase transition in ReS2 which converts back to the initial phase with time, (ii) the initial and the “switched” phases are identified as ({rm{T}}^{primeprime}) and ({rm{H}}^{prime}), respectively, with the latter being metastable (as opposed to its unstable H counterpart), and (iii) defects, particularly sulfur vacancies, play a crucial role in facilitating the (1{rm{T}}^{primeprime})-1({rm{H}}^{prime}) transition. Our investigation provides a deeper understanding of the possible phases of ML and BL ReS2 and reveals the importance of vacancy defect formation for the mechanisms of phase switching.

Results and discussion

System under study

The main premise of our work is depicted in Fig. 1. Under intense fs-laser illumination or influence of high-energy electrons, (1{rm{T}}^{primeprime}) ML and (2{rm{T}}^{primeprime}) BL flakes undergo a phase transition into what we hypothesize to be the (1{rm{H}}^{prime}) phase, a detailed description of which is provided in the Electronic band structure calculation section below. This phase transition is facilitated by the generation of a sufficient amount of sulfur vacancies that allow the transition to occur by providing the necessary geometrical freedom. The ({rm{H}}^{prime}) phase is metastable and is protected by an activation energy barrier (Ea), illustrated schematically in Fig. 1a. Successful realization of this phase transition is visually evident in a significantly amplified SHG signal that strongly depends on the crystalline symmetry of the atomic structure, as shown in Fig. 1b for ML ReS2 using DFT calculations. The atomic rearrangement during the phase transition from (1{rm{T}}^{primeprime}) to (1{rm{H}}^{prime}) is illustrated in Fig. 1c.

Fig. 1: Overview.
Defect-assisted reversible phase transition in mono- and few-layer ReS2

a Schematic of the physical process underpinning the laser-induced phase transition. Femtosecond laser pulses in the near-IR range induce a phase transition from the ({rm{T}}^{primeprime}) ground state of ReS2 flakes into the ({rm{H}}^{prime}) transient state manifested by a structural change. b SHG polarization-resolved intensity of (1{rm{T}}^{primeprime}) and (1{rm{H}}^{prime}) monolayer flakes calculated using DFT. c Structural changes in the lattice symmetry upon the phase transition between (1{rm{T}}^{primeprime}) and (1{rm{H}}^{prime}) monolayer flakes of ReS2.

Full size image

Our experimental investigation starts with the fabrication of ML and BL ReS2 flakes on two different substrates via mechanical exfoliation and dry transfer from a host crystal. Using this approach, we deposited flakes onto: (i) a silicon/silicon dioxide (Si/SiO2) substrate for optical and field-effect transistor (FET) transport studies, and (ii) a silicon nitride (SiN) membrane for STEM and electron diffraction studies. More details about the fabrication process are provided in the “Methods” section.

The thickness of the fabricated flakes was identified in three steps: optical imaging, SHG intensity, and atomic-force microscopy (AFM). Optical contrast in reflection microscopy provides an accurate estimation of the thickness of the flake64,65. Owing to different symmetries upon stacking, the ReS2 flakes demonstrate SHG activity depending on the number of layers62. Finally, AFM reveals the exact thickness of the most promising flakes (see Supplementary Information (SI), Supplementary Fig. 10).

Evidencing the phase transition in mono- and bilayer ReS2 by electron microscopy

To directly characterize the unperturbed starting phase of ML and BL ReS2, as well as the evolution of their atomic structure during defect generation, we conducted experimental and simulated atomic resolution annular dark-field scanning transmission electron microscopy (ADF STEM) and position averaged convergent beam electron diffraction (PACBED) investigations. ADF STEM imaging enables one to directly visualize the positions of the Re atoms (see “Methods” section for more details about the selection of STEM imaging modes and resulting image contrast) and, therefore, the atomic structure evolution. PACBED patterns provide a wealth of structural information about crystalline samples66,67, such as the phase, by recording momentum-resolved measurements of the forward scattered electrons. PACBED enables the collection of this information from extremely small areas of the sample (down to a single unit cell or even atom68), which can be controlled by selecting the sample area where the STEM probe is scanned. PACBED is particularly important to this work because it is sensitive to the electron scattering from both the Re and S atoms, and thus provides a S-sensitive phase characterization tool, in contrast to ADF STEM imaging. To enable these investigations, exfoliated ML and BL flakes were transferred to ~20 nm thick SiN membrane TEM grids (see “Methods” section for further details).

We identify the structure of unperturbed ML ReS2 using ADF STEM imaging and PACBED while considering 5 possible ML ReS2 phases: 1H, 1T, (1{rm{T}}^{prime}), (1{rm{T}}^{primeprime}), and (1{rm{H}}^{prime}). The experimental atomic resolution ADF STEM image of unperturbed ML ReS2 in Fig. 2a and the associated template-averaged image in Fig. 2a,i reveals the positions of the Re atoms and the presence of a double-distorted structure that is indicated by the red diamonds. A simulated ADF STEM image of the (1{rm{T}}^{primeprime}) phase is shown in the lower right inset image in Fig. 2a,ii. The similarity between the experimental (Fig. 2a,i) and simulated (Fig. 2a,ii) ADF STEM images, as well as the dissimilarity between the experimental ADF STEM image and simulated images of the other possible ML ReS2 phases (see Supplementary Fig. 3), provide evidence that unperturbed ML ReS2 has the (1{rm{T}}^{primeprime}) phase. In addition, the experimental PACBED pattern from ML ReS2 (Fig. 2b,i) is in best agreement with the simulated PACBED pattern from the (1{rm{T}}^{primeprime}) phase (Fig. 2b,ii and Supplementary Fig. 3) providing further evidence that unperturbed ML ReS2 has the (1{rm{T}}^{primeprime}) phase. Our findings are in agreement with previous studies69, further supporting that monolayer ReS2 adopts the (1{rm{T}}^{primeprime}) configuration.

Fig. 2: ML and BL ReS2 STEM observations.
figure 2

a Experimental ADF STEM image of undamaged ML ReS2. The scale bar is 5 nm. The lower left inset image is a template-matched average image produced from the entire image in (a). The lower right inset image is a simulated ADF STEM image of the (1{rm{T}}^{primeprime}) phase. The scale bar is 1 nm in both inset images. b(i) Experimental PACBED pattern of undamaged ML ReS2. b(ii) Simulated PACBED pattern of the (1{rm{T}}^{primeprime}) phase. The scale bar is 0.6 inverse Angstroms in both b(i) and b(ii). c Experimental ADF STEM images from different regions of electron-irradiated ML ReS2. The scale bar is 1 nm. c(i) is from an undamaged (1{rm{T}}^{primeprime}) region. c(ii–iv) are from different damaged regions that exhibit a loss of the (1{rm{T}}^{primeprime}) phase but maintain a dominant distorted axis that rotates by 60° around the out-of-plane c-axis to produce three possible in-plane orientations. The green markers represent the rearrangement of Re atoms upon phase transition. d Experimental ADF STEM image of undamaged BL ReS2. The lower left inset image is a template-matched average image produced from the entire image in (d). The lower right inset image is a simulated ADF STEM image of the (2{rm{T}}^{primeprime})e phase. The scale bar is 1 nm in both inset images. e(i) Experimental PACBED pattern of undamaged BL ReS2. e(ii) Simulated PACBED pattern of the (2{rm{T}}^{primeprime})e phase. The scale bar is 0.6 inverse Angstroms in both e(i) and e(ii). f Experimental ADF STEM images extracted from an image series showing the atomic structure evolution of BL ReS2 due to electron beam irradiation. The scale bar is 1 nm. f(i) shows the initial undamaged structure that is consistent with the (2{rm{T}}^{primeprime})e phase. f(ii–iv) shows three distinct structures that are observed in sequence during electron beam irradiation, which are followed by complete destruction of the 2D crystalline lattice. The green markers represent the rearrangement of Re atoms upon phase transition.

Full size image

Electron beam irradiation is known to introduce damage to TMD materials, primarily through a knock-on sputtering process that introduces chalcogenide vacancies70. For the unperturbed phase measurements, special experimental design and analysis were conducted to ensure there were negligible effects from the electron beam irradiation on the phase identification (see “Methods” section for more details). However, we also leveraged the chalcogenide vacancies production mechanism to design ADF STEM imaging experiments that resolve the atomic structure evolution of ML and BL ReS2 during increasing S-vacancy concentrations. This is achieved by acquiring time-resolved ADF STEM image series with carefully selected electron beam energy, electron dose, and STEM image scan parameters (see “Methods” section).

During time-resolved ADF STEM imaging experiments on ML ReS2 (see Supplementary Movies 1 and 2 for example image series), electron irradiation damage occurs very rapidly, partially because experiments are limited to using a 200 keV electron beam due to sample constraints (see “Methods” section for more details). Similar behavior was observed in ML MoTe2 during lower voltage TEM experiments71. The image series shown in Supplementary Movies 1 and 2 reveal that ML ReS2 undergoes a complete destruction of its 2D crystallinity after a few tens of image frames (see “Methods” section for dose values). After a complete loss of the 2D crystalline structure, it appears as though Re nanoclusters remain, providing evidence that the damage occurs through a sulfur sputtering and vacancy generation mechanism. Interestingly, just prior to the loss of the 2D crystallinity, ML ReS2 undergoes a rearrangement of its atomic structure from a double-distorted (1{rm{T}}^{primeprime}) phase to what appears like nano-sized regions of a single distorted nonstoichiometric phase that can have three possible in-plane orientations. Figure 2c,i–iv shows extracted images from areas adjacent to and after Supplementary Movie 2. Figure 2c,i shows the starting (1{rm{T}}^{primeprime}) phase, while Fig. 2c,ii–iv show the three possible orientations of the resulting single distorted structure that are rotated by 60° from each other. Supplementary Fig. 2 provides details on the image locations of Fig. 2c,i–iv and a Fast Fourier transform analysis of these nano-sized regions that confirms a clear change in structure and a 60° rotation between possible single distorted phase regions. Because the ADF STEM images are only sensitive to the Re atoms and not the S atoms, we can not unambiguously identify the resulting single distorted phase using this measurement. However, this investigation reveals ML ReS2 transitions from the (1{rm{T}}^{primeprime}) phase to a single distorted nonstoichiometric phase that is consistent with (1{rm{H}}^{prime}) (Supplementary Fig. 3). Furthermore, the (1{rm{H}}^{prime}) phase is consistent with our SHG data and DFT calculations, as we show below.

We now use ADF STEM imaging and PACBED to identify the structure of unperturbed BL ReS2. For this analysis, we consider 15 possible BL ReS2 phases, including 2H, 2T, (2{rm{T}}^{prime}), and twelve different (2{rm{T}}^{prime}) phases (named (2{rm{T}}^{primeprime})c-n) that have previously been explored as possible stable phases72. The experimental atomic resolution ADF STEM image of unperturbed BL ReS2 in Fig. 2d and the associated template-averaged image in Fig. 2d,i reveals the presence of a stacked double-distorted structure that is indicated by the red and blue diamonds. A simulated ADF STEM image of the (2{rm{T}}^{primeprime})e phase is shown in Fig. 2d,ii. The similarity between the experimental, Fig. 2d,i, and simulated, Fig. 2d,ii ADF STEM images, as well as the dissimilarity between the experimental ADF STEM image and simulated images of the other possible BL ReS2 phases (see Supplementary Figs. 7 and 8), provides evidence that unperturbed BL ReS2 has the (2{rm{T}}^{primeprime})e phase. In addition, the experimental PACBED pattern from BL ReS2, Fig. 2e,i, is in best agreement with the simulated PACBED pattern from the (2{rm{T}}^{primeprime})e phase (Fig. 2e,ii, Supplementary Fig. 6, Supplementary Fig. 9) providing further evidence that unperturbed BL ReS2 has the (2{rm{T}}^{primeprime})e phase. After considering all 12 possible phases, we can conclude unperturbed BL ReS2 has the (2{rm{T}}^{primeprime}) phase and most likely the (2{rm{T}}^{primeprime})e variant.

Time-resolved ADF STEM imaging experiments also provide evidence for electron beam-induced phase transitions in BL ReS2. Supplementary Movies 7 and 8 (see “Methods” section for data processing details and Supplementary Movies 5-6 for unprocessed versions) are two consecutive image series of the same area that show BL ReS2 exhibits obvious structural changes during electron beam irradiation before eventually losing its 2D crystallinity. This is similar to what was observed for ML ReS2, except BL ReS2 requires a substantially higher electron dose for a complete loss of 2D crystallinity (see “Methods” section for dose values). Additionally, it undergoes what appears like three distinct atomic structure transitions before it loses its 2D crystallinity. Fig. 2f,i–iv are extracted images from Supplementary Movies 7 and 8 at four different times (see “Methods” section). Fig. 2f,i shows the initial (2{rm{T}}^{primeprime})e phase characterized by the two stacked red and blue diamonds. After electron irradiation and the generation of S-vacancies, the structure transitions through 3 distinct phases shown in Fig. 2f,ii–iv. Fig. 2f,ii and iii show the atomic structure transitions through what appears like two different single distorted phases. Following this, the final phase is characterized by the loss of the double and single distortion axes of the Re atoms and the presence of a hexagonal structure (see green hexagon in Fig. 2f,iv). Based on the atomic column contrast in Fig. 2f,iv, each column appears to contain two Re atoms stacked in an A-A stacking configuration. Because the S atoms are invisible in these ADF STEM images and they are critical to identifying an H or T phase, it is not possible to unambiguously identify the transitional phases from this data (the T phase is, however, inconsistent with our SHG data, discussed below). Despite this, we can conclude the electron irradiation and the subsequent generation of S-vacancies causes BL ReS2 to evolve from a starting double-distorted (2{rm{T}}^{primeprime}) phase into nano-sized domains of two different nonstoichiometric single distorted phases before losing its Re distortion axes completely and transitioning into an A-A stacked hexagonal nonstoichiometric structure with a subsequent loss of its 2D crystallinity, leaving behind Re nanoclusters. Although we can not extract the S-vacancy concentration throughout this evolution, we expect the S-vacancy concentration to continuously increase during this process with the total electron dose, highlighting the critical role of defects in enabling the ({rm{T}}^{primeprime})({rm{H}}^{prime}) transition.

Light-induced phase transition

We now turn our attention to light-induced phase transition in ReS2. In these experiments, we use a Si/SiO2 substrate with a 285-nm-thick SiO2 layer. We start our analysis by studying the time-dependent formation of the SHG-active (1{rm{H}}^{prime}) phase. In this experiment, the fs-laser beam was focused on a specific location on the ML ReS2 flake, with a new spot selected for each excitation power. The SHG intensity was continuously monitored as a function of time during the exposure. Fig. 3a depicts the time evolution of the SHG intensity signal from different areas of the ReS2 ML flake. The colors of the curves from blue to red demonstrate the increase of the incident 1040 nm wavelength fs-laser power. One can see similar features in every curve for different powers: (i) some time delay, or, as we introduce, the lag time τlag before the SHG intensity begins to substantially increase (the rate is the time in the power of three), (ii) flattening or saturation of the signal after reaching a certain level of counts. The analysis of the lag time τlag or the inverted lag time ({tau }_{{rm{lag}}}^{-1}) as a function of incident power will give us information on how soon the phase transition begins. Fig. 3b depicts exactly this with a clear third-order dependence of the incident power, meaning doubling the power will result in an 8-fold quicker initiation of the phase transition of the ReS2 ML on a 285 nm Si/SiO2 substrate. This finding is of interest for patterning extended ReS2 areas, as shown in comparison in the insets of Fig. 3b before exposure (upper) and after exposure (lower), since it allows us to precisely determine the necessary and sufficient conditions for the phase transition to occur. The reason behind such two-phase behavior and the dependence of the lag time on fs-laser power might arise from two key processes – multiple photon absorption of the flake for transition into the (1{rm{H}}^{prime}) phase and the phase transition nucleation and growth process, which requires defect formation. As we discuss in detail below, the latter is in line with the increase in the SHG intensity signal of the patterned ReS2 bilayer with the number of cycles of fs-laser patterning that the particular sample underwent, as shown in Fig. 5a – crosses. Additionally, it has been reported that ReS2 MLs tend to host point defects more likely than, for example, MoS2 as the covalent bond is rather soft in the (1{rm{T}}^{primeprime}) double-distorted phase73,74. Further evidence for defect formation comes from the FET data (see Supplementary Fig. 13). Notably, the current in the patterned structures drops by nearly two orders of magnitude, indicating defect formation. This observation is in line with previous studies that report degradation in transport properties due to the formation of sulfur vacancies in e.g. MoS275.

Fig. 3: Defect-induced phase transition in ReS2 flakes.
figure 3

a Time evolution of SHG counts for ReS2 monolayer at 1040 nm fs-laser-induced point (static position of the laser beam) patterning process at different incident powers. τlag is the specific time of the SHG counts dramatic increase. The dashed line corresponds to the third-power slope. The inset depicts the optical image of the flake, with the black square indicating the area of interest. The scale bar is 5 μm. b The rate of the phase transition (inverted τlag) as a function of incident power of the fs-laser at 1040 nm pump wavelength. The color of the experimental points corresponds to the same color as in the time evolution curve in (a). The dashed line corresponds to the third-power slope. Insets depict the SHG map before (upper) and after (lower) laser point patterning. c and d are polarization-resolved SHG polar plots obtained from respective areas of (e) – red cross (e.ii – pristine) and red diamond (e.iii – exposed to moderate UV light). The flake after the formation of fs-laser-induced defects studied in (a, b) is subsequently exposed to UV light (e.iii). e also shows an optical image of the flake (e.i).

Full size image

Once the ML is heavily patterned by fs-laser (bottom inset Fig. 3b) and, therefore, hosts enough defects, one can induce a phase transition from (1{rm{T}}^{primeprime}) into (1{rm{H}}^{prime}) phase using an ultraviolet (UV) source at significantly smaller power compared to the fs-laser (see “Methods” section). Notably, without such pre-patterning by the fs-laser, UV light alone can not trigger the phase transition (see Supplementary Fig. 12), suggesting that it is insufficient to generate defects. This also suggests that switching and defect formation processes can be effectively decoupled, allowing UV-induced switching to occur without the need for further defect generation. In addition, we stress that the UV experiment was performed across the entire flake using an unfocused light, affecting both areas directly exposed to the fs-laser pre-patterning and those farther away. Both regions successfully transitioned to the new phase, supporting the nucleation and growth hypothesis underlying this phase transition. The polarization-resolved SHG plots in Fig. 3c, d before and after patterning manifest significant changes in the structure upon UV exposure. The pronounced 3 pairs of 2-fold symmetric SHG profiles with 60° shifts between each other in the patterned case indicate the formation of a different phase. The three 2-fold symmetries are in line with the STEM observations in Fig. 2c and data provided in SI, which indicate that upon phase transition ML ReS2 forms distinct domains with the b-axes in each domain oriented 60° with respect to each other. Individually, each domain yields an SHG spectrum with 2-fold symmetry, as calculated in Fig. 1b, but there is no guarantee that in a given SHG measurement the surface area of the three mutually rotated (1{rm{H}}^{prime}) domains will be equal. Hence, in general, the SHG signal is expected to be characterized by three 2-fold symmetric profiles. Moreover, the intensity of the SHG from the ML in Fig. 3e confirms the phase transition. The AFM data (see Supplementary Fig. 10) of the studied ReS2 flake confirms its monolayer nature.

We validate the defect-assisted nature of the light-induced phase transition in ReS2 through various experimental techniques, including fs-laser patterning, SHG polarization-resolved spectroscopy, UV-induced phase switching, and photoluminescence (PL) (the latter is discussed in SI, Supplementary Fig. 11). Furthermore, using the Si/SiO2 substrate allows for conducting both optical and FET experiments and therefore gaining a deeper understanding of the electronic conductivity properties (see Supplementary Fig. 13). We note that the final phase of the ReS2 can not be unambiguously identified solely from the SHG data, however, we conclude that ReS2 switches into either H or ({rm{H}}^{prime}) due to SHG intensity and polarization-resolved data pointing to the broken inversion symmetry of the lattice. Together with ADF STEM presented earlier and the subsequent discussion on DFT, we propose the hypothesis of ({rm{T}}^{primeprime})({rm{H}}^{prime}) phase transitions, as this scenario aligns consistently with all the experimental and theoretical methods used in our study as well as with observations from previously published work63. Specifically, our observations indicate that a certain level of defects is necessary to facilitate the ({rm{T}}^{primeprime})({rm{H}}^{prime}) phase transition. These defects are induced either by fs-laser irradiation or exposure to energetic electrons, presenting an excited state pathway for the ({rm{T}}^{primeprime})({rm{H}}^{prime}) transition. The (1{rm{H}}^{prime}) phase is metastable due to the presence of an activation energy barrier, and once defects are introduced, the transition from the (1{rm{T}}^{primeprime}) phase to the (1{rm{H}}^{prime}) phase can be initiated by various external stimuli, such as UV light, fs-laser pulses, or fast electrons. The formation of defects is irreversible, as evidenced by the irreversibility of the PL and FET measurements, which do not revert to their initial performance (see SI for details). However, the phase change in terms of SHG activity can be almost completely reversed at elevated temperatures, as discussed below.

Electronic structure calculations

Our STEM and SHG data indicate a phase transition to a state with broken inversion symmetry. A reasonable first assumption is that the transition involves the standard undistorted H phase, similar to that of monolayers MoS2, WS2, and others63. However, electronic structure calculations reveal that this phase is unstable in monolayer ReS252,69, see Fig. 4. This suggests that another phase, exhibiting broken inversion symmetry, likely underlies the observed SHG and STEM behaviors. STEM data further reveal that the “switched” phase is distorted (see Fig. 2c), confirming its distinction from the 1H phase. Hence, we propose that this phase corresponds to (1{rm{H}}^{prime}) (the distorted 1H phase; space group Amm2, see Fig. 1c). Our numerical simulations, incorporating phononic and electronic band structure diagrams (see Fig. 4), indicate the stability of the (1{rm{H}}^{prime}) phase and its potential accessibility from (1{rm{T}}^{primeprime}) via high-energy excitations utilizing photons (both multi- and single-photon processes are possible) or electrons (see Fig. S2). These calculations confirm that the 1H phase is unstable and therefore unlikely to be responsible for the observed behavior. Instead, ReS2 tends to switch into a metastable (1{rm{H}}^{prime}) semimetal state, which aligns with our experimental observations in ADF STEM and SHG. To our knowledge, this is the first demonstration of the existence and stability of the (1{rm{H}}^{prime}) phase in ReS2.

Fig. 4: DFT analysis of monolayer ReS2 phases.
figure 4

Left: stable (1{rm{T}}^{primeprime}) semiconducting phase. Middle: unstable 1H phase. Right: metastable (1{rm{H}}^{prime}) semimetal phase.

Full size image

Our calculations suggest the involvement of sulfur vacancies in facilitating the ({rm{T}}^{primeprime})({rm{H}}^{prime}) phase transition. Furthermore, they also suggest the reduced energy required to form sulfur vacancies in the (1{rm{H}}^{prime}) phase compared to (1{rm{T}}^{primeprime}). The role of the sulfur defects consists in providing the necessary space for the transition to occur. Importantly, both the direct ({rm{T}}^{primeprime}) to ({rm{H}}^{prime}) and reverse ({rm{H}}^{prime}) to ({rm{T}}^{primeprime}) phase transitions require synchronized motion of large numbers of sulfur atoms, which can only occur when there is sufficient geometrical freedom to perform such motion. Additionally, molecular dynamics simulations demonstrate that the back-transition from (1{rm{H}}^{prime}) to (1{rm{T}}^{primeprime}) requires S-vacancies and represents a first-order phase transition involving nucleation and growth of one (switched) phase within the continuum of another (unswitched) phase (data not shown), confirming the role of sulfur vacancies. These predictions align with both light- and electron-induced (1{rm{T}}^{primeprime})(1{rm{H}}^{prime}) phase transitions observed in SHG, STEM, and FET measurements, as all of them require and indicate defect formation prior to the transition (for instance, see UV-induced phase change data, Fig. 3c–e, and Supplementary Figs. 12 and 13). Finally, it is important to note that the reverse transition from (1{rm{H}}^{prime}) to (1{rm{T}}^{primeprime}) can not be induced by light due to the significant differences in their electronic band structures and corresponding absorption spectra (see Fig. 4).

The SHG pattern of the proposed (1{rm{H}}^{prime}) phase can be simulated by using the χ(2) tensor obtained from density functional theory (DFT) and yields a polarization-dependent profile with two distinct lobes (Fig. 1b). Our experimental SHG data suggest the existence of three possible orientation domains within the (1{rm{H}}^{prime}) phase (rotated with respect to each other by 60°). Their stochastic combination accounts for a 6-fold-like symmetric pattern (Fig. 3d), where SHG lobes vary in length, reflecting the prevalence of specific orientations in the final orientation domain composite (see STEM data in Fig. 2c).

Although numerical electronic structure calculations offer substantial support for the existence and relative stability of the (1{rm{H}}^{prime}) phase, alongside information about sulfur vacancies and second-order optical nonlinearities, they lack precision in estimating the activation energy barrier for the reverse ({rm{H}}^{prime}) to ({rm{T}}^{primeprime}) transition, which is crucial to characterize the stability of the ({rm{H}}^{prime}) phase. Thus, in the subsequent section, we shift our focus to experimental observations of this back-switching process. The metastability of the ({rm{H}}^{prime}) phase implies its eventual inevitable return to the stable ({rm{T}}^{primeprime}) phase, with the switching rate following a temperature-dependent Arrhenius-like behavior, as we demonstrate below.

Estimation of the activation energy barrier in bilayer ReS2

We use a pristine BL flake of ReS2 to study the activation energy barrier height for the reverse transition from the ({rm{H}}^{prime}) to the ({rm{T}}^{primeprime}) phase. Unlike an ML (see Fig. 3), the BL is highly SHG-active62,63, which is reflected in the upper inset with a cyan frame of Fig. 5a. Upon fs-laser patterning of the area with specific step size (50 nm), exposure time (400 ms), and power (8 mW), the flake significantly reduces the counts of the SHG signal, indicating the phase transition from the ({rm{T}}^{primeprime}) into ({rm{H}}^{prime}) phase (see Fig. 5a, lower inset with a magenta frame). After some time, the sample recovers its initial SHG intensity counts, signaling the reverse transition to the pristine state. By applying the exposure conditions multiple times after each recovery, we repeated the patterning process 7 times for the bilayer to prove the reversibility of the phase transition, as shown in Fig. 5a. The circles depict the recovered (or pristine for the first cyan circle) SHG-active ({rm{T}}^{primeprime}) bilayer states, whereas the crosses mark the patterned SHG-inactive ({rm{H}}^{prime}) ones. Notably, we observe a clear trend that the difference in SHG signal between the circles and crosses is reducing with the number of cycles. This may be attributed to incomplete phase switching caused by irreversible degradation or permanent damage of the material triggered by multiple exposures to the fs-laser.

Fig. 5: Fs-laser-induced patterning mechanism.
figure 5

a The average SHG intensity over the flake upon several patterning-recovery cycles of a bilayer ReS2 flake. Circles correspond to the counts of SHG-active (2{rm{T}}^{primeprime}) of a bilayer, while crosses correspond to SHG counts of a (2{rm{H}}^{prime}) phase. The right side panels (highlighted in cyan and magenta) demonstrate the SHG maps of a bilayer before and after the 1st patterning, respectively. All the patterning cycles were carried out at the same conditions of the fs-laser – 8 mW average power, 400 ms exposure, and a 50 nm spatial step. The recovery conditions are specified above each recovered point. The colored regions correspond to the conditions of the recovery of the same color in (b) and (c). b Recovery FOM(t,T) as a function of time for different temperature values that stimulates the recovery process. The lines indicate a single exponential fit for the recovery process. c Arrhenius plot showing logarithm of recovery rate (1/τ) vs. inverted temperature at which the recovery process was measured. The linear fitting reflects the energy activation barrier, Ea = 1.10 ± 0.53 eV, between the phases. Colors and symbols correspond to the same colors and symbols as in (b). Inset is the sketch of the energy diagram corresponding to the thermally-induced ({rm{H}}^{prime})({rm{T}}^{primeprime}) recovery process.

Full size image

Another important parameter is the time the BL ReS2 flake resides in the metastable ({rm{H}}^{prime}) phase. Importantly, the recovery time is strongly temperature-dependent (see Fig. 5a). At room temperature, the sample takes several days to return to its initial state, whereas at elevated temperatures, the relaxation process occurs within minutes. By measuring the duration of this recovery process at different temperatures, we aim to conduct an Arrhenius-like analysis to estimate the activation energy barrier (Ea) between the ground ({rm{T}}^{primeprime}) phase and the metastable ({rm{H}}^{prime}) phase. With that in mind, we introduce the phase transition recovery figure-of-merit (FOM), which relies on the SHG intensity as a measure of the phase transition. If the Irec is the SHG intensity of the recovered bilayer flake, then FOM(t,T) = (frac{{I}_{{rm{T}}}(t)}{{I}_{{rm{rec}}}}), where T is the temperature of the flake, t is the time after the patterning took place, and IT(t) – the intensity of the SHG of the flake at a particular moment t and at specific temperature T. For a range of temperatures from 60 °C to 80 °C with a 5 °C step size, we performed a series of the measurements of the FOM(t,T) as shown in Fig. 5b, for the same sample from Fig. 5a. FOM = 0 right at the moment of a freshly patterned flake ((2{rm{H}}^{prime})), whereas FOM = 1 means that the sample recovered to its initial SHG intensity ((2{rm{T}}^{primeprime})). The sample recovery curves at different temperatures can be fitted with a single exponential growth function – FOM(t,T) = 1 − et/τ, with a characteristic time scale, τ. The rate of the reverse phase transition is, therefore, τ−1. The normalized SHG intensities as a function of recovery time (symbols), along with the corresponding fits (solid lines), are shown in Fig. 5b for various temperatures. Then, for every particular temperature, we plot the logarithm of τ−1 multiplied by the value of the Boltzmann constant. Furthermore, we fit the ({rm{-ln}}(tau )) versus 1/T using a linear function and obtain Ea = 1.10 ± 0.53 eV with 95% confidence (see Fig. 5c). The experimental point corresponding to 70 °C was excluded from this analysis due to an insufficient number of measurements, resulting in large errors in this case. The fitting procedure includes the error bars of the extracted values of τ for every temperature condition. Based on these experimental findings, we deduce that the stabilization energy barrier for the (2{rm{H}}^{prime}) phase is of the order of 1 eV. The phase transition diagram is schematically depicted in Fig. 5c. Initially, the ({rm{H}}^{prime}) phase can not form due to the insufficient number of defects. However, upon fs-laser-induced multiphoton absorption, combined with the simultaneous creation of S-vacancies, a new excited state relaxation pathway opens up, enabling the formation of the metastable ({rm{H}}^{prime}) phase through the first-order phase transition, which involves nucleation and growth.

In this study, we present evidence of a reversible phase transition occurring in mono-, bi-, and trilayer ReS2. The transition occurs from the ({rm{T}}^{primeprime}) to the ({rm{H}}^{prime}) phase at room temperature and ambient conditions. Through a combination of several experimental and theoretical methods, including electron microscopy, second-order nonlinear spectroscopy, photoluminescence, electronic transport, and DFT calculations, we identify that the transition is consistent with the formation of the metastable ({rm{H}}^{prime}) (distorted H) phase. Notably, our data suggest a crucial role of sulfur vacancies in enabling this transition. Our analysis also confirms the ground state phase of the ReS2 ML and BL flakes to be (1{rm{T}}^{primeprime}) and (2{rm{T}}^{primeprime}), respectively. Furthermore, we experimentally determine the activation energy barrier for the reverse ({rm{H}}^{prime}) to ({rm{T}}^{primeprime}) transition in BL ReS2 to be approximately 1.10 ± 0.53 eV. From a broader perspective, these results expand the knowledge of light- and electron-induced phase-switching materials and offer insights into the potential for optically programmable properties of TMD-based devices.

Methods

Sample fabrication

The ReS2 flakes were mechanically exfoliated from the high-quality bulk crystal (HQ Graphene) using scotch tape with a subsequent transfer on polydimethylsiloxane stamps (PDMS) in a cleanroom (Nanofabrication laboratory at Chalmers within MyFab). For the STEM experiments, ReS2 flakes were transferred to 20 nm thick SiN membrane TEM grids (simpore.com) using an all-dry transfer method76. For the laser patterning and all the optical measurements, the flakes were transferred from the PDMS stamps to the glass slides. For the electrical transport measurements, ReS2 flakes were transferred to the 285-nm thick Si/SiO2 highly n-doped substrates (Graphene Supermarket, USA) with pre-fabricated 20 nm Cr/200 nm Au back contact. To produce source and drain top contacts for ReS2 flakes, a ~300 nm layer of 950 PMMA A4 resist (MicroChem, USA) was applied by spin-coating. PMMA was patterned by e-beam lithography employing a Raith EBPG 5200 (Germany) system operating at 100 kV accelerating voltage and 30 nA current. A dose of 1300 μC/cm2 was applied, and 1:3 MIBK:IPA mixture was used as a developer. Then, the 5 nm Cr/200 nm Au layer was evaporated by Lesker PVD 225 e-beam evaporation system (Kurt J. Lesker Company, Germany). Overnight lift-off in acetone was used to finalize the contact fabrication.

Optical measurements

Second-harmonic generation (SHG) measurements were carried out using (690–1040 nm) Ti:sapphire femtosecond laser (MaiTai HP-Newport Spectra-Physics) with a ~100 fs pulse duration and 80 MHz repetition rate as an excitation source. The light was polarized linearly with a broadband linear polarizer and focused through an objective (Nikon, 40×, NA = 0.95). The positioning of the sample was controlled both in the plane (xy) and the focus (z) by the piezoelectric motorized stage. The signal was collected with an avalanche photo-diode for mapping (APD, IDQ, ID100 Visible Single-Photon Detector) or spectrometer (Andor 500i, equipped with Newton 920 CCD camera) through an optical fiber. The polarization-resolved SHG measurements were performed on the same setup using λ/2 plate to rotate the linear polarization of the excitation light. For SHG imaging, we utilized a tightly focused 1040 nm laser beam with 2 mW power, employing a 250 nm scanning step and a 5 ms exposure time per spot. Under these imaging conditions, the ReS2 phase transition was not triggered.

Photoluminescence (PL) measurements were carried out on the inverted optical microscopy setup using a high-NA oil immersion objective (Nikon, 60×, NA = 1.49). The excitation 447 nm continuous-wave laser was coupled to an optical fiber and focused on the bottom side of the sample through the objective. A 488 nm long-pass filter is placed in the detection path to remove the laser excitation from the PL signal. The setup allows for rotating the incident linear polarization and offers the possibility for an analyzer in the signal collection channel.

Light-induced phase transition conditions

In the case of fs-laser-induced transitions, defect formation, and the light-induced phase transition occur at the same time. For patterning the ReS2 flake, we utilized a tightly focused (close to diffraction limit) 1040 nm laser beam with 8 mW power (average intensity of ~106 W/cm2, peak intensity of ~1.25 × 1011 W/cm2), 80 MHz repetition rate, and a 100 fs pulse duration. The sample patterning conditions were consistent across all fs-laser experiments, employing a 50 nm scanning step and a 400 ms exposure time per spot. These conditions resulted in significantly higher laser exposure compared to the SHG-imaging conditions specified above and were sufficient to induce the phase transition.

For phase transition using UV illumination, ReS2 flakes were exposed to an unfocused quasi-CW 237 nm nanosecond laser with an average intensity of 0.01 W/cm2 (peak intensity of 1182 W/cm2) for one hour.

STEM measurements and calculations

STEM imaging and PACBED experiments were carried out on a JEOL Mono NEO ARM 200F microscope. ADF STEM imaging was conducted using an accelerating voltage of 200 kV, a probe convergence half-angle of 26.7 mrad, an ADF detector collection angle range of 65–146 mrad, and a probe current of 44 picoampere (pA). In order to improve the quality of the ADF STEM images for the atomic structure analyses, a combination of non-rigid registration of image series77 and template matching was used depending on the specific data set. PACBED experiments66 were conducted using an accelerating voltage of 80 kV, a probe convergence half-angle of 6.35 mrad, and a probe current of <3 pA. To resolve the ML and BL ReS PACBED patterns more clearly, careful corrections of the detector dark signal and SiN support signal were manually performed. ADF STEM and PACBED simulations were calculated using the frozen phonon multislice algorithm implemented in the GPU-accelerated software Prismatic 2.078.

The ADF STEM collection angle range of 65–146 mrad produces Z-contrast images79,80,81 where primarily the Re atom positions are visible while the S positions are invisible. Other collection angles more suitable for light-element imaging, such as annular bright field STEM imaging82, yielded poor results due to the signal from the SiN support.

The accelerating voltage of 200 kV was selected because lower voltages made it more difficult to achieve sufficient atomic contrast of the Re atoms above the SiN support background intensity. The accelerating voltage of 200 kV is above the known knock-on threshold energy for sulfur vacancy generation in TMD materials70. However, the 44 picoamperes (pA) probe current and 4–32 μs/pixel dwell times were sufficient to image ReS2 before significant structural damage occurred. These beam conditions enabled the observation of both the initial undamaged atomic structure and the dynamic evolution of the atomic structure as sulfur vacancies were generated.

For PACBED experiments, an accelerating voltage of 80 kV was beneficial for two primary reasons. First, it minimized damage to the sample and enabled long PACBED exposure times, which were required to reveal the weak ML and BL diffraction signals above the noise level. 80 kV is slightly above the knock-on damage threshold for most TMD materials70, but subsequent PACBED acquisitions of the same area showed no signs of a changing structure and damage. Second, lowering the accelerating voltage increases the Bragg scattering angles of the ReS2 diffraction signals. This effectively separates the diffraction disks, causes less disc overlap, and makes it easier to identify PACBED patterns from different phases. PACBED patterns were acquired by scanning the STEM beam over a 10–20 nm2 area for 300 s while continuously acquiring the diffraction pattern. This acquisition technique enables sufficient sampling and position-averaging of probe positions over the structure to produce high-quality PACBED patterns, and it distributes the electron dose over a larger area to minimize sample damage. To see the ML and BL ReS2 PACBED patterns more clearly, careful corrections of the dark signal and SiN support signal were manually performed. The dark signal was corrected by subtracting a reference dark pattern that was acquired directly after the ReS2 PACBED pattern. The SiN support adds a diffuse amorphous signal to the PACBED pattern. This SiN support signal was removed by subtracting a reference pattern that was acquired from only the SiN support near the flake. The dark and SiN reference patterns were acquired with the same acquisition parameters as the ReS2 patterns. The resulting PACBED patterns were binned by 4 × 4 pixels to enhance the signal-to-noise ratio and gain more clarity in the diffraction signals.

The ADF STEM image shown in Fig. 2a is a single-pass STEM image acquired with a size of 1024 × 1024 pixels, a dwell time of 32 μs/pixel, a pixel size of 0.039 nm/pixel, a flyback time of 500 μs, and a total electron dose of 6 × 104 e Å−2. All of the doses reported in this publication include the dose associated with the line flyback time. The template-matched average image shown in the lower left inset image of Fig. 2a is the average of 3195 template images from within Fig. 2a. The template images that were selected for the averaging were determined using an in-house written IGOR Pro script that locates the maxima in the normalized cross-correlation between a user-selected input template image and all other possible template images (see Supplementary Fig. 1 for more details).

The ADF STEM images shown in Fig. 2c are sub-regions extracted from a single-pass STEM image that was acquired directly after a higher magnification image series which damaged the ML. The single-pass image is shown in Supplementary Fig. 2 and was acquired with a size of 1024 × 1024 pixels, a dwell time of 64 μs/pixel, a pixel size of 0.048 nm/pixel, a flyback time of 500 μs, and a total electron dose of 7.6 × 104 e Å−2. Supplementary Fig. 2 also shows the locations of the cropped images displayed in Fig. 2c and their associated FFTs, which further show the changes in atomic structure and orientation. The higher magnification image series that damaged the ML is shown in Supplementary Movie 2 and was acquired at the center of the single-pass image shown in Supplementary Fig. 2. The image series was acquired with 100 frames, a size of 512 × 512 pixels, a dwell time of 16 μs/pixel, a pixel size of 0.02 nm/pixel, a flyback time of 500 μs, and a total electron dose of 1.2 × 107 e Å−2. An additional ADF STEM image series that reveals the dynamics of the structural changes caused by electron irradiation is shown in Supplementary Movie 1. It was acquired with 100 frames, a size of 1024 × 1024 pixels, a dwell time of 4 μs/pixel, a pixel size of 0.016 nm/pixel, a flyback time of 500 μs, and a total electron dose of 4.8 × 106 e Å−2. The series displayed in Supplementary Movie 1 was spatially binned by 4 × 4 pixels and smoothed with a 3 × 3 pixels Gaussian to improve the signal-to-noise ratio and enhance the visibility of the atomic structure.

The ADF STEM image shown in Fig. 2d is an average of an image time series that was non-rigidly registered77 to correct the effects of the microscope, sample, and environmental instabilities. The image series was acquired with 50 frames, a size of 512 × 512 pixels, a dwell time of 8 μs/pixel, a pixel size of 0.04 nm/pixel, a flyback time of 500 μs, and a total electron dose of 7.9 × 105 e Å−2. The non-rigidly registered and a temporally binned version are shown in Supplementary Movies 3 and 4, respectively. These movies, in combination with Supplementary Fig. 4, show that negligible structural damage was caused by the electron irradiation. The template-matched average image shown in the lower left inset image of Fig. 2d is the average of 843 template images from within Fig. 2d that were determined using the same process discussed above for Fig. 2a (see Supplementary Fig. 5 for more details).

The ADF STEM images shown in Fig. 2f are sub-regions extracted from frames of an image time series that damaged the bilayer. The total image series is composed of two separate consecutively acquired image series of the same sample area, each of which was acquired with 50 frames, a size of 512 × 512 pixels, a dwell time of 16 μs/pixel, a pixel size of 0.0395 nm/pixel, a flyback time of 500 μs, and a total electron dose of 1.5 × 106 e Å−2. Each image series was non-rigidly registered77 to correct the effects of the microscope, sample, and environmental instabilities. Both series were temporally binned by 5 frames to improve the signal-to-noise ratio and the visibility of the structural changes, producing an image series with 10 frames. The images shown in Fig. 2f (main text) are extracted regions from (i) frame 1 of 10 in series 1, (ii) frame 10 of 10 in series 1, (iii) frame 4 of 10 in series 2, and (iv) frame 10 of 10 in series 2. These 4 images have been smoothed with a 3 × 3 pixel Gaussian to improve the signal-to-noise ratio and enhance the visibility of the atomic structure. The non-rigidly registered image series are provided in Supplementary Movies 5 and 6 and the temporally binned image series are provided in Supplementary Movies 7 and 8.

ADF STEM and PACBED simulations were calculated using the frozen phonon multislice algorithm implemented in the GPU-accelerated software Prismatic 2.078. DFT relaxed atomic models of the different ML and BL phases were used as input structures for the calculations. All calculations utilized root-mean-square displacements of 0.045 Å and 0.077 Å for the Re and S atoms, respectively. ML ADF STEM image simulations were performed with an accelerating voltage of 200 keV, a STEM probe convergence half-angle of 26.69 mrad, 10 phonon configurations, atomic models with a size of around 150 Å × 150 Å, a potential space sampling of 0.04 Å, an ADF detector range of 50–155 mrad, and a probe step size of 0.1 Å. ML PACBED simulations were performed with an accelerating voltage of 80 kV, a STEM probe convergence half-angle of 6.35 mrad, 10 phonon configurations, atomic models with a size of around 150 Å × 150 Å, a potential space sampling of 0.05 Å, and a probe step size of around 1.6 Å, which is sufficient for Nyquist sampling. Bilayer ADF STEM image simulations were performed with an accelerating voltage of 200 kV, a STEM probe convergence half-angle of 26.69 mrad, 10 phonon configurations, atomic models with a size of around 50 Å × 50 Å, a potential space sampling of 0.03 Å, an ADF detector range of 50–155 mrad, and a probe step size of 0.05 Å. Bilayer PACBED simulations were performed with an accelerating voltage of 80 kV, a STEM probe convergence half-angle of 6.35 mrad, 30 phonon configurations, atomic models with a size of around 100 Å × 100 Å, a potential space sampling of 0.1 Å, and a probe step size of around 1.6 Å, which is sufficient for Nyquist sampling.

Electrical transport measurements

FET measurements were conducted using two source channel -measurement units (B2901A, Keysight) connected through a General Purpose Interface Bus (GPIB) cable with a computer. The typical FET measurements were performed at room temperature. The back-gate voltage was systematically adjusted within the range of −10 V to 10 V, while simultaneously tuning the source-drain voltage between −1 V and 1 V.

Electronic structure calculations

The energetics of different phases, their dynamic stability, defect formation energies, and electronic band structures were obtained via DFT calculations that were carried out using the projector augmented-wave method83 as implemented in the Vienna ab initio simulation package84,85. The exchange-correlation potential was represented using the van der Waals density functional with consistent exchange method86,87. All calculations were carried out using a plane-wave energy cutoff of 520 eV and for the calculation of the forces, a finer support grid was employed to improve their numerical accuracy. The Brillouin zone was sampled using Γ-centered grid with a linear (overrightarrow{k})-point spacing of about 0.25 Å−1 and Gaussian smearing with a width of 0.1 eV. This setup has been previously employed and tested for other layered materials88,89.

For the ideal 1A and Amm2 structures the electronic band structures were furthermore computed using the G0W0 method90 on a Γ-centered 8 × 8 × 1 grid. These calculations employed 282 empty bands, thereby including states up to at least 28 eV above the valence band edge (1 A) or Fermi level (Amm2), and used 50 imaginary frequency and time grid points.

Phonon dispersions were computed via the PHONOPY package91 in combination with the HIPHIVE package92 using force constants from DFT calculations.

The χ(2) polarizability tensor was calculated using the method and tools outlined by Taghizadeh et al.93. Briefly, we calculated the nonlinear optical response using the DFT open-source GPAW package94,95. We used the PBE exchange-correlation functional96, the Kohn–Sham orbitals were expanded using a plane-wave basis set with an energy cutoff of 500 eV, we used a 21 × 21 × 3 Monkhorst–Pack grid for the k-mesh and a grid spacing of less than 0.2 Å. The number of empty bands included in the sum over bands was set to three times the number of occupied bands. We used Fermi–Dirac occupation number smearing with a factor of 20 meV and line-shape broadening of 50 meV. Time-reversal symmetry was imposed in order to reduce the k-integrals to half.

Related Articles

Solution-processable 2D materials for monolithic 3D memory-sensing-computing platforms: opportunities and challenges

Solution-processable 2D materials (2DMs) are gaining attention for applications in logic, memory, and sensing devices. This review surveys recent advancements in memristors, transistors, and sensors using 2DMs, focusing on their charge transport mechanisms and integration into silicon CMOS platforms. We highlight key challenges posed by the material’s nanosheet morphology and defect dynamics and discuss future potential for monolithic 3D integration with CMOS technology.

Photovoltaic bioelectronics merging biology with new generation semiconductors and light in biophotovoltaics photobiomodulation and biosensing

This review covers advancements in biosensing, biophotovoltaics, and photobiomodulation, focusing on the synergistic use of light, biomaterials, cells or tissues, interfaced with photosensitive dye-sensitized, perovskite, and conjugated polymer organic semiconductors or nanoparticles. Integration of semiconductor and biological systems, using non-invasive light-probes or -stimuli for both sensing and controlling biological behavior, has led to groundbreaking applications like artificial retinas. From fusion of photovoltaics and biology, a new research field emerges: photovoltaic bioelectronics.

Memristors based on two-dimensional h-BN materials: synthesis, mechanism, optimization and application

Memristors offer vast application opportunities in storage, logic devices, and computation due to their nonvolatility, low power consumption, and fast operational speeds. Two-dimensional materials, characterized by their novel mechanisms, ultra-thin channels, high mechanical flexibility, and superior electrical properties, demonstrate immense potential in the domain of high-density, fast, and energy-efficient memristors. Hexagonal boron nitride (h-BN), as a new two-dimensional material, has the characteristics of high thermal conductivity, flexibility, and low power consumption, and has a significant application prospect in the field of memristor. In this paper, the recent research progress of the h-BN memristor is reviewed from the aspects of device fabrication, resistance mechanism, and application prospect.

First-principles and machine-learning approaches for interpreting and predicting the properties of MXenes

MXenes are a versatile family of 2D inorganic materials with applications in energy storage, shielding, sensing, and catalysis. This review highlights computational studies using density functional theory and machine-learning approaches to explore their structure (stacking, functionalization, doping), properties (electronic, mechanical, magnetic), and application potential. Key advances and challenges are critically examined, offering insights into applying computational research to transition these materials from the lab to practical use.

3D printing of micro-nano devices and their applications

In recent years, the utilization of 3D printing technology in micro and nano device manufacturing has garnered significant attention. Advancements in 3D printing have enabled achieving sub-micron level precision. Unlike conventional micro-machining techniques, 3D printing offers versatility in material selection, such as polymers. 3D printing technology has been gradually applied to the general field of microelectronic devices such as sensors, actuators and flexible electronics due to its adaptability and efficacy in microgeometric design and manufacturing processes. Furthermore, 3D printing technology has also been instrumental in the fabrication of microfluidic devices, both through direct and indirect processes. This paper provides an overview of the evolving landscape of 3D printing technology, delineating the essential materials and processes involved in fabricating microelectronic and microfluidic devices in recent times. Additionally, it synthesizes the diverse applications of these technologies across different domains.

Responses

Your email address will not be published. Required fields are marked *