Raman and coherent phonon spectroscopies were used to investigate the thickness-dependent phononic properties of ultrathin single-crystal Bi films prepared by molecular beam epitaxy on Si(111) substrates. Both the A1g and Eg Raman peaks disappeared in the Raman spectra of a 4 nm Bi film, indicating a complete transition from the low-symmetry A7 structure to the high-symmetry A17 structure. Coherent phonon signals of the A1g mode also showed a strong dependence on the film thickness, where thin samples (≤15 nm) exhibited lower phonon frequency and shorter phonon lifetimes than the thick samples (≥30 nm). This difference is attributed to a shallower energy potential barrier caused by both a permanent phase transition, which is determined by the film thickness, and a temporary structural transition by photo-excited carriers. Our results not only provide evidence of a phase transition from the A7 to the A17 structure with the decreasing Bi film thickness but also reveal the influence of this phase transition on phonon dynamics. Understanding these material performance traits will facilitate modern application of Bi thin films in electronic devices.
Bismuth, with a Peierls distortion structure,1 has been extensively studied since the 1960s because of its unique atomic and electronic structures. The slight overlap between the conduction and valence bands makes it a semimetal with a low number of carriers, thus an exceptionally long electron mean free path at room temperature (∼2 μm).2 In recent years, with the development of a suitable growth technique, high quality bismuth films can be obtained through molecular beam epitaxy (MBE). This innovation has attracted renewed interest to investigate fundamental physics in this prototype semimetal, such as ultrafast bond softening,3 exceptional surface-state spin and valley properties,4,5 semimetal-semiconductor (SM-SC) transition,6,7 superconductivity,2 topological insulator states, 2D spintronic phenomena,8–11 transient high-symmetry phase transformation,12 and strong anharmonic scattering.13,14 All these properties, combined with a negative real and small imaginary part of the dielectric constant and the strong inter-band transition, make the MBE-grown Bi thin film promising for application in inter-band plasmonics.15
Usually, a bulk Bi crystal finds natural stability in an A7 structure with a Rhombohedral primitive cell containing two Bi atoms per cell, and an intrinsic crystal growth direction is along hexagonal (001) [Fig. 1(a)]. Some different phases were observed in the early stages of material growth by MBE, where the first several monolayers (MLs) of Bismuth were contrived to adapt to the crystal structure of the substrate. Recently, Fang et al. found a temporary phase showing periodic regions when growing Bi on the NbSe2 substrate with MBE.16 Back in 2004, Nagao et al. conducted extensive studies to show that with a Si(111) substrate, the first four MLs had a pseudocubic puckered 2D structure, similar to that of black phosphorus (BP).17,18 This puckered structure could either be the real A17 structure [Fig. 1(b)] or the A7 structure along the {012} direction [Fig. 1(c)]. Reflection high-energy electron diffraction (RHEED), which is generally used to monitor the crystal structure in situ during an MBE growth, could not distinguish A17 from A7 {012} because they produce very similar diffraction patterns due to the same out-of-plane lattice constant.18 As Bi films grow thicker, they gradually transition to the conventional A7 (001) structure. The threshold for the number of MLs for this transition to occur spans from 4 MLs to 7 MLs, the true value of which is dependent on the substrate temperature. The BP-like A17 structure possesses a cubic unit cell with only one atom per unit cell. Due to the fundamental differences between the A17 and A7 crystal structures, they may be distinguished through their phononic properties. When a crystal has two atoms per primitive cell, it has 6 phonon branches: 3 acoustic and 3 optical. For A7 Bi, the signature optical branches are the A1g and Eg phonons, which polarize perpendicular and parallel to the hexagonal (001) plane, as described in Fig. 2(a). For the A17 structure, there are only three acoustic branches because of containing one atom per primitive cell. As a result, no A1g or Eg phonons can be observed in the A17 structure.
Schematics of bismuth crystal structures in top view: (a) A7 (rhombohedral) structure along the hexagonal (001) direction; (b) A17 (pseudocubic, black-phosphorus) structure; and (c) A7 (rhombohedral) structure along the hexagonal {012} direction. The A17 and A7 {012} structures share the same lattice constant and atomic density, therefore producing a very similar diffraction pattern in the RHEED spectra.
Schematics of bismuth crystal structures in top view: (a) A7 (rhombohedral) structure along the hexagonal (001) direction; (b) A17 (pseudocubic, black-phosphorus) structure; and (c) A7 (rhombohedral) structure along the hexagonal {012} direction. The A17 and A7 {012} structures share the same lattice constant and atomic density, therefore producing a very similar diffraction pattern in the RHEED spectra.
(a) Schematic diagram of the A1g (longitudinal, orange arrow) and Eg (transverse, blue arrow) modes in an A7 rhombohedral unit cell. The coordinates are in the hexagonal indexing, i.e., z axis is hexagonal (001), corresponding to the (111) plane in rhombohedral indexing. At equilibrium, the Bi atom inside the unit cell (green atom) sits at (0.48:0.52) along the diagonal axis. The dashed circle marks the position of the center along the diagonal axis. (b) Raman spectra of Bi thin films ranging from 4 nm to 80 nm. The signature Eg and A1g modes in 80 nm and 50 nm thin films are marked by orange vertical lines as references.
(a) Schematic diagram of the A1g (longitudinal, orange arrow) and Eg (transverse, blue arrow) modes in an A7 rhombohedral unit cell. The coordinates are in the hexagonal indexing, i.e., z axis is hexagonal (001), corresponding to the (111) plane in rhombohedral indexing. At equilibrium, the Bi atom inside the unit cell (green atom) sits at (0.48:0.52) along the diagonal axis. The dashed circle marks the position of the center along the diagonal axis. (b) Raman spectra of Bi thin films ranging from 4 nm to 80 nm. The signature Eg and A1g modes in 80 nm and 50 nm thin films are marked by orange vertical lines as references.
We performed Raman spectroscopy on the MBE synthesized Bi films on Si(111) varying from 4 nm to 80 nm (see Figs. S1–S3) utilizing an inViaTM confocal Raman microscope (Renishaw) operating with 532 nm excitation. These results are shown in Fig. 2(b), presenting several noteworthy features. Raman spectra of 80 nm and 50 nm totally overlap, indicating that in this thickness range, Bi films show bulk properties. When the material thickness decreases, both A1g and Eg peaks shift to higher frequency with decreasing intensity. This is usually attributed to the strain generated at the Bi/Si interface due to the lattice mismatch. The subsequent Raman spectra of 15 nm and 8 nm samples proved to be almost the same. At 6 nm, both A1g and Eg peaks become much broader and shift a little bit back to a lower frequency. At 4 nm, both A1g and Eg peaks suddenly disappear completely. The disappearance of both the A1g and Eg peaks in a 4 nm thick film is a strong evidence of phase transition from an A7 to an A17 structure, considering the fact that the A17 structure does not have any optical phonon branches. The complex features at 6 nm might come from the transition region from A7 to A17. Using Raman Spectroscopy, we confirmed the phase transition in the extremely thin Bi film grown on the Si(111) substrate. A simple picture to describe this transition process would be the Bi atom inside the A7 primitive cell [green atom in Fig. 2(a)], which moves gradually from the off-center position along the diagonal axis to the center (dashed circle). At the same time, the Rhombohedral cell gradually shifts to a pseudocubic shape. A transition region from the pure A7 to A17 phase is expected, which is neither standard A7 nor A17. We can call this region the pseudocubic phase. Here, the spatial extension of the transition region depends on the growth conditions.
Transient reflectivity change of bismuth thin films on the silicon substrate at fluences of 0.36 mJ/cm2 (a), 1.08 mJ/cm2 (b), and 1.81 mJ/cm2 (c).
Transient reflectivity change of bismuth thin films on the silicon substrate at fluences of 0.36 mJ/cm2 (a), 1.08 mJ/cm2 (b), and 1.81 mJ/cm2 (c).
To further investigate the influence of this phase transition process on phonon dynamics, we conducted coherent phonon spectroscopic measurements. Degenerate pump-probe experiments were performed in a non-collinear reflection geometry at room temperature using a mode-locked Ti:Sapphire femtosecond laser (Spitfire ACE, Spectra Physics). Both pump and probe pulses have a central wavelength of 800 nm and a pulse width of ∼258 fs (FWHM) coming onto the sample surface at a repetition rate of 5 kHz. The absorption depth for bismuth at 800 nm is 14.7–15 nm.19,20 Pump and probe beams were focused onto the sample surface with spot sizes (diameter at the intensity) of 84 μm and 42 μm, respectively. The pump beam was modulated at 585 Hz by a mechanical chopper, and the signal from the probe beam was recorded using a lock-in amplifier (model 7265, Signal Recovery). Coherent phonons (CPs) have been used extensively to probe material phase change processes, even at ultrafast time scales.21–23 Figure 3 shows the CP signals in all samples under different pump fluences. In our experiment, we only observed A1g mode (see the FFT spectrum in Fig. S5). For all fluences, there is a sign change, from positive to negative, of the first peak when the film thickness is below 15 nm. Using the transfer matrix method (see Fig. S6),24 we have found that this sign change is a result of multi-reflection of light waves when the film thickness is comparable with an optical penetration depth (∼15 nm in Bi at 800 nm). The raw data consist of both CP signals and electronic background, and to extract the pure CP signals, we applied a digital bandpass filter to remove the slowly varying electronic background. Then, the CP signals can be fitted using damping harmonic oscillators to obtain CP properties, such as the phonon lifetime () and phonon frequency (f),
where , , , , and are the amplitude, dephasing time, initial frequency, the linear chirp rate, and initial phase, respectively. For the 4 nm sample, containing no obvious signals, neither electron nor CP oscillations were observed for the pump fluence up to 1.81 mJ/cm2 [Fig. 3(c)]. The signal at a higher fluence of 3.98 mJ/cm2 is from Si (see Fig. S7), and the Bi thin film has been damaged.
Figure 4 presents the CP lifetimes and frequencies as a function of pump fluences, F [Figs. 4(a) and 4(b)]. In Bismuth, it is well accepted that CPs are excited through the DECP process (displace excitation of coherent phonons).25 From the extracted phonon frequency and dephasing time, the excited phonon mode is the optical A1g. From the Raman spectra in Fig. 2(b), both Eg and A1g phonon modes disappear in the 4 nm sample. Consequently, there is no surprise in the disappearance of CP signals in the 4 nm sample [Fig. 3(c)]. CP measurements also confirm that there is a complete phase transition in the 4 nm sample. Figures 4(a) and 4(b) show that both τ and f decrease with pump fluence in all samples. However, the decreasing trend clearly shows two distinct groups. For samples with thicknesses of 80 nm, 50 nm, and 30 nm, the trends indicating a decrease with pump fluence (dτ/dF, df/dF) are almost the same. For samples of 15 nm, 8 nm, and 6 nm, the decreasing trends with fluence are more than 3× faster.
Fluence-dependent coherent A1g dephasing time, , (a) and initial frequency, , (b) obtained from fitting for different bismuth thicknesses ranging from 6 nm to 80 nm. The solid lines in (a) and (b) are given as guides to the eye.
Fluence-dependent coherent A1g dephasing time, , (a) and initial frequency, , (b) obtained from fitting for different bismuth thicknesses ranging from 6 nm to 80 nm. The solid lines in (a) and (b) are given as guides to the eye.
The significant difference between the two groups of samples can be attributed to the phase transition from the A7 to A17 structure. Beyond 30 nm, the pump/probe pulses mainly excite/probe the region of the pure A7 phase with an optical penetration depth of δ ∼15 nm at 800 nm. When the film thickness goes down to a value comparable to δ, the laser pulses start to excite/probe the pseudocubic phase (transition region) and the A17 phase partially. The pure A17 phase could extend up to a thickness of 4 nm, depending on the growth conditions.17,18,26 Since the pure A17 phase does not contribute to the CP signals, the differences between these two groups could come from the pseudocubic phase. First, the faster decrease in phonon frequency with F indicates stronger bond softening effects. In Fig. 2(a), the atom inside the primitive cells moves toward the center in the pseudocubic phase. During this process, the total system energy of the pseudocubic phase should be higher than that of both the A7 and A17 phases. The potential field of the pseudocubic phase has a larger anharmonicity than the A7 and A17 phases, which could be disturbed more easily through the use of an external electromagnetic field. Second, the faster decrease in τ indicates stronger scattering. In the pseudocubic phase, the phonon frequencies at different thicknesses are expected to have a continuous change when transitioning from A7 to A17 structures, as they occur gradually. Naturally, this smooth transition would bring more phonon–phonon scattering channels.
Beside the phase transition relationship driven by the film thickness, another possible contribution of the observed phenomenon is a light-driven temporary phase transition. Teitelbaum et al. reported a real-time observation of a high-symmetry phase in both 20 nm and 275 nm bismuth films using single-shot pump-probe spectroscopy.12 When more than 2% of the total valence electrons are excited, where 1% corresponds to a carrier density of 1.43 × 1021/cm3, it has been proposed that the potential barrier of the A7 structure could become shallower so much, so that Peierls distortion could be temporarily removed.12,27 A temporary phase transition from A7 to a pseudocubic structure is possible. The main evidence of phase transition in Teitelbaum's work is the fast decrease in phonon frequency with pump fluence and the eventual disappearance of coherent phonon signals. Samples used in Teitelbaum's work were prepared using a sputtering technique on sapphire and glass substrates. Hence, all their samples are polycrystalline, having a pure A7 phase where there is no permanent phase transition with the thickness. When they extrapolate the fluence-dependent phonon frequency to zero pump fluence, the results in both thick and thin samples tend to converge to the same point. It was concluded that the phonon ground states in all their samples are the same. In our results, we observed a similar phenomenon for samples of 80 nm, 50 nm, and 30 nm, signifying that these samples share the same phonon ground state. However, in the thin films (especially in 8 nm and 6 nm), the extrapolated phonon frequency at F = 0 is higher than that of the thick samples. This indicates that the phonon ground state in these thin films is different from that of the thick films. This also signifies a permanent phase transition with the thickness. According to the estimation of the average excited carrier density in our samples,19 under the fluence of 1.08 , the photoexcited carrier densities are 0.5 × 1021/cm3, 0.7 × 1021/cm3, 1.1 × 1021/cm3, 1.6 × 1021/cm3, 1.9 × 1021/cm3, and 2.0 × 1021/cm3 for 80 nm, 50 nm, 30 nm, 15 nm, 8 nm, and 6 nm, respectively. Due to the fact that we used regular pump-probe spectroscopy to excite and detect the coherent phonons, where multiple photons are absorbed, the samples were damaged before the pump fluence reached a threshold to completely remove the potential barrier of Peierls distortion. However, the much faster down-shift of phonon frequency in thinner films (15 nm, 8 nm, and 6 nm) with pump fluence indicates that the effect of a shallower potential barrier is much more obvious at these excited carrier densities. Considering the existence of a permanent phase transition region in these thin films, where the potential barrier is already flatter than the pure A7 structure, it is justifiable to recognize the fast decreasing trend of phonon frequency in relation to a relatively low carrier density.
In summary, we have demonstrated a layer-dependent phase transition from the A7 to the A17 crystal structure in MBE-grown Bi thin films via Raman spectroscopy. Both the A1g and Eg phonon modes disappear abruptly in 4 nm thick films, indicating a complete transition to the A17 phase. Ultrafast coherent phonon dynamics bifurcate into two distinct groups: thick (80 nm, 50 nm, and 30 nm) and thin (15 nm, 8 nm, and 6 nm) samples. Thin samples (≤15 nm) exhibited lower phonon frequency and shorter phonon lifetimes than the thick samples (≥30 nm). We attribute this behavior to two causes: the quasi-cubic phase in the transition region and the temporary phase transition driven by photo-excited carriers. Our results not only provided evidence of a transition to the A17 phase in the ultrathin MBE-grown Bi film on the Si(111) substrate but also reveal the influence of this phase transition on phonon dynamics. Understanding these material aspects is imperative to the facilitation of Bi thin film applications in cutting-edge electronic devices.
See the supplementary material for the sample information, Transfer Matrix Method, and FFT spectrum.
We acknowledge the help of Xianghai Meng with Raman measurements in Jung-fu Lin's laboratory. The authors acknowledge support from the National Science Foundation (NASCENT, Grant No. EEC-1160494; CAREER, Grant Nos. CBET-1351881 and CBET-1707080), the Center for Dynamics and Control of Materials (Grant No. DMR-1720595), and the Semiconductor Research Corporation and Texas Instruments Fellowships.
DATA AVAILABILITY
The data that support the findings of this study are available within this article and its supplementary material.