As a new two-dimensional layered material, black phosphorus (BP) is a promising material for nanoelectronics and nano-optoelectronics. We use Raman spectroscopy and first-principles theory to report our findings related to low-frequency (LF) interlayer breathing modes (<100 cm -1 ) in few-layer BP for the first time. The breathing modes are assigned to Ag symmetry by the laser polarization dependence study and group theory analysis. Compared to the high-frequency (HF) Raman modes, the LF breathing modes are much more sensitive to interlayer coupling and thus their frequencies show much stronger dependence on the number of layers. Hence, they could be used as effective means to probe both the crystalline orientation and thickness for fewlayer BP. Furthermore, the temperature dependence study shows that the breathing modes have a harmonic behavior, in contrast to HF Raman modes which are known to exhibit anharmonicity. 2 Department of Physics, Applied Physics, and Astronomy, Rensselaer Polytechnic Institute, Troy, New York 12180, USA.3 Center for Nanophase Materials Sciences, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, USA.
2Orthorhombic black phosphorus (BP) is the most stable allotrope of phosphorus. It features a layered structure with puckered monolayers stacked by van der Waals (vdW) force 1 . Few-or single-layer BP can be mechanically exfoliated from bulk BP [2][3][4][5] . Due to BP's intrinsic thicknessdependent direct bandgap (ranging from 0.3 eV to 2.0 eV) and relatively high carrier mobility (up to ~1,000 cm 2 V -1 s -1 at room temperature) 2,4,[6][7][8][9] , it is expected to have promising applications in nanoelectronic devices [2][3][4]10 , and near and mid-infrared photodetector [11][12][13][14][15][16][17][18][19] . Recently, high performance thermoelectric devices were also predicted based on BP thin films [20][21][22] . With the surge of interest in two-dimensional (2D) materials (such as graphene and transition metal dichalcogenides (TMDs)) 23,24 , BP has become a new attraction since 2014 because it bridges the gap between graphene and TMDs, and offers the best trade-off between mobility and on-off ratio 3 . Moreover, the unique anisotropic puckered honeycomb lattice of BP leads to many novel in-plane anisotropic properties, which could lead to more applications based on BP 3,9,25,26 .Phonon behaviors play an important role in the diverse properties of materials 27 , which has been intensively studied in vdW layered materials, such as graphene and TMDs [28][29][30][31][32][33][34][35] . Raman spectroscopy is a powerful non-destructive tool to investigate the phonons and their coupling to electrons, and has been successfully applied to vdW layered materials [36][37][38][39][40] . Due to the lattice dynamics of vdW layered compounds, the phonon modes in these materials can be classified as high-frequency (HF) intralayer modes and low-frequency (LF) interlayer modes 27 . Intralayer modes involve vibrations from the intralayer chemical bonds (Fig. 1c), and the associated frequencies ...