key: cord-0939601-rlnqubwb authors: Zhang, Chi; Tan, Junyang; Pan, Yikun; Cai, Xingke; Zou, Xiaolong; Cheng, Hui-Ming; Liu, Bilu title: Mass Production of Two-Dimensional Materials by Intermediate-Assisted Grinding Exfoliation date: 2020-07-14 journal: Natl Sci Rev DOI: 10.1093/nsr/nwz156 sha: 6885ddf1086b66065e78d7d883e90a2b2eb27480 doc_id: 939601 cord_uid: rlnqubwb The scalable and high-efficiency production of two-dimensional (2D) materials is a prerequisite to their commercial use. Currently, only graphene and graphene oxide can be produced on a ton scale, and the inability to produce other 2D materials on such a large scale hinders their technological applications. Here we report a grinding exfoliation method that uses micro-particles as force intermediates to resolve applied compressive forces into a multitude of small shear forces, inducing the highly-efficient exfoliation of layer materials. The method, referred to as intermediate-assisted grinding exfoliation (iMAGE), can be used for the large-scale production of many 2D materials. As an example, we have exfoliated bulk h-BN into 2D h-BN with large flake sizes, high quality and structural integrity, with a high exfoliation yield of 67%, a high production rate of 0.3 g h-1 and a low energy consumption of 3.01x10^6 J g-1. The production rate and energy consumption are one to two orders of magnitude better than previous results. Besides h-BN, this iMAGE technology has been used to exfoliate various layer materials such as graphite, black phosphorus, transition metal dichalcogenides, and metal oxides, proving its universality. Molybdenite concentrate, a natural low-cost and abundant mineral, was used as a demo for the large-scale exfoliation production of 2D MoS2 flakes. Our work indicates the huge potential of the iMAGE method to produce large amounts of various 2D materials, which paves the way for their commercial application. The large number of two-dimensional (2D) materials, including graphene, hexagonal boron nitride (h-BN), transition metal dichalcogenides (TMDCs) like MoS2 and WSe2, metal oxides (MxOy), black phosphorene (b-P), etc, provide a wide range of properties and numerous potential applications, [1] [2] [3] [4] [5] [6] [7] but achieving these requires scalable production. Bottom-up strategies like chemical vapor deposition (CVD) and chemical synthesis have been extensively explored but only small quantities have been produced so far. [8, 9] For example, the highest production of CVD-grown graphene is on the order of one million square meters (1,000,000 m 2 ) per year, which is equal to ~0.77 kg, while for other CVD-grown 2D materials it is much less. [10, 11] Moreover, the high growth temperatures (450-1100 °C) and time-consuming transfer processes usually needed for CVD-grown 2D materials result in high cost and a low production rate. [11] As for the chemical synthesis strategy, solvothermal treatment has been used to synthesize different kinds of 2D materials including not only TMDCs [12] but also even hard-for-CVD-synthesised MxOy [13] and b-P [14] , however, the quantity is still small. Another strategy for the preparation of 2D materials is top-down exfoliation, which involves exfoliating bulk layer materials to a monolayer or few-layer flakes. [3, 7, [15] [16] [17] [18] Micromechanical exfoliation using the Scotch ® tape can produce 2D materials with the highest quality, but in negligible amounts. [19] Although there are other techniques to exfoliate layer materials which can be scaled-up, they are only suitable for specific materials and/or require harsh conditions. [2, 3, 18, 20] So far, only graphene and graphene oxide can be prepared at the ton level and the difficulties in preparing other 2D materials in such a large quantity has significantly hindered their commercialization. [11] Among the current exfoliation techniques, ball milling [17, 21] and liquid phase exfoliation [22, 23] which use shear forces to induce slipping between adjacent layers, are promising for the scalable production of 2D materials other than graphene, such as h-BN and MoS2. However, ball milling usually needs a long time (e.g., 20 hours) to achieve high-yield (85%) exfoliation and the 2D materials obtained are small (~100 nm in lateral size) due to the large impact force of the irregularly moving hard balls [17] On the other hand, liquid phase exfoliation can produce 2D materials of high quality in a large scale, [22, 24, 25] but has a very low yield (< 3%). In addition, materials produced have low concentrations (<0.1 mg mL -1 ) [24, 25] in organic/surfactant solvents, and may need additional condensation and/or separation to concentrate them for use. It is obvious 4 that both ball milling and liquid phase exfoliation still face many challenges for the scalable production of good quality 2D materials. An ideal exfoliation method for the commercial scalable production of 2D materials should meet the following requirements. First, the exfoliated products should have a large lateral size and high quality with good structural integrity, because size and quality are the two main parameters that influence the material's properties. Second, the method should be efficient, with a high exfoliation yield and a high production rate. Third, the process should be cheap, green and require a low energy input, in particular when considering large scale production. Fourth, it would be ideal if the method could be universal for the exfoliation of different layer materials. Unfortunately, such a method is still not available. Note that layer materials like graphite and MoS2 possess low coefficients of friction and have been used as solid lubricants in industry for more than half a century. [26] They can buffer a large compressive force (Fc) between rolls by releasing a frictional force (Ff), transferred from Fc, through slip between adjacent layers in their structure. [27] By increasing Fc, Ff increases accordingly and eventually leads to the failure of the lubricant, which means irreversible sliding of the layers, i.e., exfoliation of these layer lubricants into thin 2D flakes. In principle, such a process would not have a big impact on the quality of the exfoliated material, because Ff is parallel to planes of the material and no additional chemicals are involved during the process. Note that a large force is not needed to exfoliate layer materials because of the weak van der Waals interactions between their layers, as revealed by experimental measurements [28] and theoretical calculations [29] . For example, for a graphite flake with a lateral size of 10 μm, slip between layers happens at a small shear force of 3.9 μN [28] In comparison, the shear force provided by milling rolls is on the order of tens to hundreds of Newtons, over 6-8 orders of magnitude higher than the force needed to exfoliate graphite. Therefore, the critical factor in increasing the exfoliation efficiency is how to effectively produce a frictional force (Ff) on a layer material. Based on the above analysis, we have developed an exfoliation technology which we call intermediate-assisted grinding exfoliation (iMAGE) that essentially meets these requirements. The key to our method is to use intermediate materials that increase the coefficient of friction of the mixture and effectively apply sliding frictional forces to the layer material, resulting in a dramatically increased exfoliation efficiency. Considering the case of 2D h-BN, the production rate and energy consumption can reach 0.3 g h -1 and 3.01×10 6 J g -1 , respectively, both of which are one to two orders of magnitude better than previous results. The resulting exfoliated 2D h-BN flakes have an average thickness of 4 nm and an average lateral size of 1.2 μm. This iMAGE method has been extended to a series of layer materials with different properties, including graphite, b-P, TMDCs, and metal oxides. To further demonstrate this scaling up, molybdenite concentrate, a naturally existing cheap and earth abundant mineral, has been used for the large-scale exfoliation production of 2D MoS2, making the commercial production and application of 2D MoS2 feasible. (Fig. 1b) . When the shear friction force ffi > bEe, where b and Ee are the width and exfoliation energy of the layer material, the layers will slip and become exfoliated (Fig. S1 ), analogous to the failure of the solid lubricants discussed above. [28] The effectiveness of the iMAGE method was evaluated by exfoliating h-BN because its exfoliation into a few layers (thickness < 5 nm) with a high yield and high quality is challenging. [17, 20, 21, 30, 31] In a typical experiment, bulk h-BN was ground with silicon carbide water, and after standing for 8 hours a green sediment containing SiC and un-exfoliated h-BN was seen, and the supernatant was milky white (Fig. S3) . The observation of a clear Tyndall effect in the supernatant indicates the colloidal state of the exfoliated 2D h-BN in DI water (Fig. 1c) . UV-vis-NIR 6 optical absorption measurements of this supernatant showed an optical bandgap of 5.8 eV based on the plot of (ahv)^2 versus photon energy hv, indicating the dispersed material is high-quality h-BN with its original structure, which is confirmed by its white color (Figs. 1c and S3). [32] To verify the production of exfoliated 2D h-BN, SiC particles and bulk h-BN were independently added to DI water, and both precipitated completely within 5 minutes (Fig. S4 ) with a clear supernatant. These results confirm the separation of 2D h-BN from the mixture by simply dispersing it in water and standing. Previous studies show that nanomaterials can be dispersed in solvents when they have a surface energy that matches that of the solvent. The surface energy of 2D h-BN is in the range 44-66 mJ m -2 [16] and that of DI Water is 72 mJ m -2 so this is not the best match, but it is the most commonly used moderate solvent and has an acceptable solvability for 2D h-BN. In contrast, SiC and h-BN particles cannot be stabilized in DI water due to a mismatch of surface energy, and quickly precipitate because they have densities larger than the solvent, i.e., 3.2 g cm -3 for SiC and 2. (Table S1 ). In addition, the results from high resolution transmission electron microscopy (HRTEM) observations and fast Fourier transformation (FFT) analysis show high quality 2D h-BN sheets (Fig. 1g ) without noticeable defects in the planes and along the edges. The crystalline quality of the 2D h-BN is comparable with that of samples prepared by liquid phase exfoliation [22] and is superior to that of the samples obtained from intercalation [33] and ball milling [21] , which have many structural defects and functional groups. The clear lattice fringes on the edge indicate negligible damage to the layer structure when transforming bulk h-BN into thin 2D sheets (Figs. 1h and S7). The large size and high quality of the obtained 2D h-BN must result from the unique exfoliation mode in the iMAGE process. During exfoliation, the SiC particles contact and interact with the top layers of bulk h-BN, causing them to slip with negligible damage to their inplane structure. Fig. 1h shows a ten-layer h-BN step (region ii) slipping from twenty-five-layer h-BN 7 (region iii), similar to previous in-situ TEM studies of the exfoliation of MoS2, where the sliding of 2D sheets from the mother MoS2 has been clearly identified under a shear force [34] , consistent with the exfoliation mechanism proposed above. Various spectroscopic characterization techniques confirm that the iMAGE-produced 2D h-BN has good quality with no noticeable functional groups (Figs. 1i-j; Supplementary Information, Section 1.1, Figs. S7 and S8). For example, X-ray photoemission spectroscopy (XPS) of 2D h-BN show two major peaks at 190.2 and 398.2 eV for B1s and N1s, the same as those for bulk h-BN, without any functional group-related peaks (Fig. 1i) . Fourier transformation infrared (FTIR) analysis directly proves the absence of functional groups on the iMAGE-produced 2D h-BN (Fig. 1j) . We note that the shapes of FTIR peaks for exfoliated 2D h-BN and bulk h-BN are different. This could be attributed to the reduced size and thickness after exfoliation. Compared to the lateral size of the exfoliated h-BN (1.2 μm), the size of the bulk h-BN (30 μm) is more close to the wavelength of infrared light in FTIR (0.2 μm ~ 20 μm), which would lead to a much stronger light scattering and cause the baseline drift as well as the peak broadening in bulk h-BN ( Fig. 1j and Fig. S8 ). In addition, after exfoliation, due to the significantly reduced thickness, the B-N stretching and bending could be excited more easily, which will sharpen the FTIR peaks of the exfoliated h-BN. [35] Moreover, the Raman peak position of the B-N vibration mode (E2g) for the exfoliated 2D h-BN is 1366 cm -1 , the same as that for bulk h-BN, also indicating negligible change to the intrinsic bonding (Fig. S7) . These results are expected because iMAGE is a purely mechanical process in which no chemical reactions take place. Based on the principle of friction force-induced sliding and exfoliation in the iMAGE technique, exfoliation occurs when ffi > bEe. Therefore, the parameters that influence ffi, b and Ee will affect the exfoliation efficiency. It is clear that the size of the force intermediate will affect the exfoliation yield. For the bulk h-BN with a 30 μm particle size used in this study, the frictional force ffi needs to be larger than 13.8 μN assuming the exfoliation energy Ee for h-BN and graphite are identical, i.e., 0.46 J m -2 . [28, 29] We calculated the frictional force on h-BN applied by each SiC particle (Supplementary Information, Section 1.4 and Table S2 ), demonstrating that only particles larger than 73 μm can exert an average frictional force greater than 13.8 μN. This result shows that the sizes of the raw bulk materials and the intermediate particles should be matched to obtain a high exfoliation yield. An exfoliation yield of 67% is confirmed by both direct measurements of the weight of 2D h-BN and 8 UV-vis absorption measurements ( Supplementary Information, Section 1.2, Section 1.3, Fig. S9 and S10) which is one of the highest reported values for the exfoliation of pristine 2D h-BN (Supplementary Table S1 ), and it was achieved in 270 minutes. Such fast high-yield exfoliation can be understood when compared to the ball milling method (Supplementary Table S3 ). In the iMAGE process, the weight of the intermediates was almost an order of magnitude lower than that of the balls used in ball milling, [31] however, the number of particles is four orders of magnitude higher and the surface area is double, both of which facilitate the effective transfer of frictional force to the layer material, leading to a shorter exfoliation time with a much higher yield. In addition to the exfoliation yield, production rate and energy consumption are the other two important parameters to be considered for the commercialization of the iMAGE method. Typically, when 2 g of bulk h-BN was used, the production rate of 2D h-BN could reach 0.3 g h -1 , which is around ten times that of a typical ball milling method for preparing amino-functionalized small 2D h-BN sheets [17] This rate is the highest value reported for h-BN production to date (Supplementary Table S1 ). The method achieves both a high yield and a high production rate (Fig. 2a) , and this is the first time a 2D material other than graphene has been produced with a yield of more than 50% and a production rate of over 0.1g h -1 . An annual production capability of 2D h-BN is expected to be exceeding 10 tons by the iMAGE technology ( Supplementary Information, Section 1.4) . For comparison, electrochemical intercalation [18] and liquid phase exfoliation [22] can achieve only a high yield or a high production rate, while ball milling may produce functionalized 2D materials with a high yield but a medium production rate. In addition, our method requires an energy consumption of only 3.01 × 10 6 J g -1 to exfoliate bulk h-BN, which is at least a tenth of that of other methods ( Fig. 2b and Supplementary Table S1 ). According to theoretical calculations, the minimum energy required to exfoliate 1 g h-BN is 279 J (Supplementary Table S4 ), three orders of magnitude lower than the value for our method, demonstrating that there is still much room for improvement. In short, the iMAGE method requires a low energy to produce 2D materials because of its efficient transfer of force and energy to the layer materials, while the production rate is among the highest, making it a practical method for the commercial scalable production of 2D h-BN. Besides SiC, several other hard micro-particles, such as Mo2C and cheap sea sand, were also used as intermediate materials to transmit the sliding fractional force onto bulk h-BN, and lead to the 9 efficient exfoliation of h-BN, demonstrating that the economical production of 2D materials could be realized by using cheap intermediates. The process would also be valuable for the one-step fabrication of 2D material-based composites in which an intermediate is a functional component in the composite material. In addition to h-BN, we have extended the method to a series of layer materials with a wide range of properties. Not only air-stable metallic graphite, n-type semiconducting Bi2Te3 and MoS2, and insulating mica, but also air-sensitive p-type semiconducting b-P were exfoliated into few layers (Figs. 3a, S12, S13 and S14). UV-vis-NIR absorption (Fig. 3b) shows the wavelength-independent absorbance for the above materials in the visible light range. For graphene, the curve is similar to the Scotch® tape-exfoliated [36] and CVD-grown materials [37] , but different from graphene oxide [4] . For as-exfoliated 2D TiOx, mica and Bi2Te3, the characteristic absorption bands are almost at the same positions as those of the raw bulk materials (Supplementary Fig. S12 ) and previous reports [38, 39] In addition, for Bi2Te3 the appearance of a prominent absorption peak at 265 nm reveals the production of few layer 2D Bi2Te3. [40] These results confirm that the iMAGE-produced 2D materials have good structural integrity rather than having been functionalized. Regarding the size and thickness of the exfoliated 2D materials, taking MoS2 as an example (Fig. S13 ), AFM measurements show that the exfoliated flakes have an average thickness of 10 nm and an average lateral size of 600 nm. The size of bulk raw MoS2 used for this exfoliation was less than 2 μm, which is small. By using MoS2 crystals with a larger size and optimizing the iMAGE process we could, in principle, obtain larger flakes. Raman analysis shows a redshift in the A1g phonon mode of 2D MoS2 compared to bulk MoS2, which is attributed to the lower interlayer van der Waals force after exfoliation [16] . The high quality of these exfoliated 2D MoS2 is confirmed by HRTEM examination and UV-vis-NIR absorption with characteristic absorption peaks at 672 and 610 nm (Fig. S13 ). [41, 42] The other layer materials (Figs. 3c-d and S12) were exfoliated into a few layers with an average thicknesses less than 10 nm and lateral sizes of 400-1000 nm under non-optimized conditions. Air-sensitive layer materials, such as b-P, have also been exfoliated by the iMAGE process in an argon-or nitrogen-filled glove box. The thickness of exfoliated 2D b-P is around 2-4 nm, as revealed by AFM characterization (Fig. 3d) . UV-vis-NIR absorption spectrum shows a broad absorption peak around 520 nm (Fig. S14) as the wavelength range used has a much higher energy than the bandgap of b-P, which is similar to the previous reported results. [43] Raman spectroscopy shows that the Ag 1 /Ag 2 ratios of exfoliated b-P and the bulk material were 0.88 and 0.68 (Fig. S14) , respectively, which are both larger than the oxidation criterion value (0.6) and indicates that the 2D b-P flakes obtained were not oxidized after iMAGE exfoliation. [44] As a step further, we used iMAGE to exfoliate abundant and cheap natural layer minerals to produce large quantities of 2D materials, for example MoS2. The industrial production of MoS2 ore follows the following procedure. 1) Finding a Mo-rich mine where the Mo concentration is higher than 0.1 wt% (Fig. 4a) . 2) Drilling and blasting to obtain an ore with a size of a few meters (Fig. 4b) . 3) Crushing the ore to a size of 10-20 μm. 4) Floating the ore to obtain a Mo-enriched molybdenite concentrate where Mo > 45 w% (Fig. 4c) . This concentrate can then be used for the iMAGE exfoliation to prepare 2D MoS2. Note that in the current molybdenum industry, the molybdenite concentrate is the flotation product of a natural ore and the precursor to industrial grade MoS2. Because of its earlier production in the manufacturing processes, a molybdenite concentrate has a much lower price (~$2.6×10 2 /t) than other molybdenum sources, e.g., 2.36% of the price of industrialgrade MoS2 (~$1.1×10 4 /t) and only 0.04% of the high-purity MoS2 chemical from several vendors (~$7.0×10 5 /t). The natural layer structure and low price suggest that molybdenite concentrate could be an ideal choice starting material for the scalable exfoliation of 2D MoS2. SEM inspection show that a bulk molybdenite concentrate does not have a well-defined layer structure because the main part is cubic particles with some layer-structured material on its surface (Fig. S15) . By using the iMAGE method we prepared 60 liters of few-layer molybdenite dispersed in water (0.25 mg mL -1 ) in less than 270 minutes (Fig. 4d) . The concentrations of 2D materials could be further increased, for example, to be >5 mg mL -1 , by using suitable solvents or adding surfactants in water. AFM and TEM shows that the lateral size of the material was 200-300 nm (Fig. 4e) . UV-vis-NIR absorption and Raman spectroscopy results confirmed the high quality of the MoS2 exfoliated from a molybdenite concentrate (Figs. 4e and S15), and were the same as those obtained from MoS2 exfoliated from MoS2 chemicals (Figs. 3 and S13) . These results show the scaling-up potential of the iMAGE method and suggest that it could be used to exfoliate cheap and abundant raw materials or natural minerals to produce high-quality 2D materials. We have developed a powerful iMAGE method to exfoliate various layer materials into 2D materials, which is suitable for commercial-scale production. The method converts applied macroscopic compressive forces into microscopic frictional forces with the assistance of an intermediate. These forces cause slipping of the layers in the layer materials resulting in highlyefficient exfoliation. Taking h-BN as an example, the iMAGE prepared 2D h-BN has a good quality, an average thickness of 4 nm (~12 layers), and average lateral size of 1.2 μm. Moreover, this method has one of the highest yields, the highest production rate, the lowest energy consumption, and the best scalable production ability of all currently-available techniques. The generality of this method has been proven by the effective exfoliation of several layer materials including graphite, Bi2Te3, b-P, MoS2, TiOx, h-BN, and mica, covering 2D metals, semiconductors with different bandgaps, and insulators. The scaling-up capability of the iMAGE technique has been demonstrated by exfoliating the cheap and abundant mineral, molybdenite, to produce 2D MoS2 in large quantities. The iMAGE method overcomes one of the main challenges, scalable production, in the 2D materials field, and is expected to significantly accelerate their commercialization for a wide range of applications. bulk TiOx crystal is based on a previous work [45] with some modifications, and the details will be published later. In short, the preparation process could be divided into two stages, including (i) alkali metal ion intercalation and (ii) protonation of the titanates. In the stage (i), the layered K0.8Ti1.73Li0.27O4 was prepared by grinding and heating a stoichiometric mixture of TiO2, K2CO3, and Li2CO3 at 1000 o C for two times. After the first calcination for 5 hours, the intermediate products need to be taken out and grinded, followed by the second calcination for 20 hours. In the stage (ii), the stage (i) obtained K0.8Ti1.73Li0.27O4 (1g) was added into HCl solution (1 M, 200 mL) and stirred for 4 days continuously. After stopping stirring, the product H1.07Ti1.73O4 was obtained by washing the sediment with plenty of water and drying at room temperature. Then, the H1.07Ti1.73O4 bulk underwent the iMAGE exfoliation process to prepare TiOx flakes. Bulk crystals of b-P were grown by mineralizer-assisted short-distance transport (SDT) reactions. Different liquids, including DI water (Millipore Milli-Q water purification system, 18.2 MΩ) and dimethylformamide (DMF, purity >99 wt%, Shanghai Macklin Biochemical Co. Ltd., China) were used without further purification. In a typical process, raw h-BN particles (2 g) and force intermediates (SiC, 8 g) were added to a fast rotation apparatus (RM 200, Germany), and ground for 270 min under a compressive force of around 100 Newtons. To exfoliate air-stable layer materials like h-BN, graphite, Bi2Te3, MoS2, TiOx and mica, the iMAGE process was performed in ambient atmosphere. To exfoliate air-sensitive layer materials, such as b-P, the iMAGE was performed in a glove box (500 × 600 × 600 mm 3 ), which was evacuated and then filled with argon or nitrogen gas twice before the exfoliation process. The as-ground mixture was added to liquids (DMF for b-P and water for the others). After standing for 8 hours, a sediment was clearly seen at the bottom of the container, and the exfoliated 2D materials were stably dispersed in the supernatant. To collect the exfoliated 2D materials from the supernant, the dispersion was centrifuged at 1500 rpm for 30 minutes, and then vacuum filtered onto an alumina membrane with a pore size of 20 nm. A UV-vis-NIR absorption spectrometer (Agilent Cary 5000, USA) was used to characterize the exfoliated 2D materials in dispersion. Optical microscopy (Carl Zeiss Axio Imager 2, Germany), atomic force microscopy (AFM, tapping mode, Oxford Instruments, USA), scanning electron microscopy (SEM, Hitachi SU8010, 15 kV, Japan), transmission electron microscopy (TEM, FEI Tecnai G2 F30, 300 kV, USA) and a laser particle size analyser (Malvern Zetasizer Nano-ZS90, UK) 13 were used to characterize the morphology and structure of the exfoliated 2D materials, including their lateral size, thickness, crystal quality, etc. Raman spectroscopy (514 nm laser with a power density of 1 mW cm -2 , Horiba LabRAM ER, Japan) was used to examine the quality of the bulk and exfoliated 2D materials. X-ray photoelectron spectroscopy (XPS, monochromatic Al Kα X-rays, 1486.6 eV, PHI VersaProbe II, Japan) and Fourier transform infrared spectroscopy (FTIR, Thermo Scientific Nicolet iS 50, USA) were used to examine the chemistry of the bulk and exfoliated 2D h-BN. In addition, XPS and powder X-ray diffraction (XRD, with monochromatic Cu Kα radiation λ = 0.15418 nm, Bruker D8 Advance, Germany) were used to examine the purity and crystallinity of the exfoliated 2D materials collected from the supernatant. Supplementary data are available at NSR online. Preparation of 2D material dispersions and their applications 2D metal carbides and nitrides (MXenes) for energy storage Carbon-based supercapacitors produced by activation of graphene Processable aqueous dispersions of graphene nanosheets Mass production and industrial applications of graphene materials Production of phosphorene nanoribbons Two-dimensional materials from high-throughput computational exfoliation of experimentally known compounds Fast growth of inch-sized single-crystalline graphene from a controlled single nucleus on Cu-Ni alloys Surface modification-induced phase transformation of hexagonal close-packed gold square sheets Precision synthesis versus bulk-scale fabrication of graphenes The global growth of graphene Morphology and surface chemistry engineering toward pH-universal catalysts for hydrogen evolution at high current density Generalized self-assembly of scalable two-dimensional transition metal oxide nanosheets Facile bottom-up synthesis of partially oxidized black phosphorus nanosheets as metal-free photocatalyst for hydrogen evolution Chemically derived graphene oxide: towards large-area thin-film electronics and optoelectronics Two-dimensional nanosheets produced by liquid exfoliation of layered materials Boron nitride colloidal solutions, ultralight aerogels and freestanding membranes through one-step exfoliation and functionalization An effective method for the fabrication of few-layer-thick inorganic nanosheets Electric field effect in atomically thin carbon films Controlled gas exfoliation of boron nitride into few-layered nanosheets Simultaneous Production and Functionalization of Boron Nitride Nanosheets by Sugar-Assisted Mechanochemical Exfoliation Scalable production of large quantities of defect-free few-layer graphene by shear exfoliation in liquids Graphene-based liquid crystal device Turbulence-assisted shear exfoliation of graphene using household detergent and a kitchen blender 2D-crystal-based functional inks Solid lubricants: a review Tribology of tungsten disulphide (WS2): characterization of wear-induced transfer films Measurement of the cleavage energy of graphite van der Waals bonding in layered compounds from advanced density-functional first-principles calculations Exfoliation of hexagonal boron nitride by molten hydroxides Scalable exfoliation process for highly soluble boron nitride nanoplatelets by hydroxide-assisted ball milling Tuning of the optical, electronic, and magnetic properties of boron nitride nanosheets with oxygen doping and functionalization Phase engineering of transition metal dichalcogenides Nanomechanical cleavage of molybdenum disulphide atomic layers A facile chemical exfoliation method to obtain large size boron nitride nanosheets Fine structure constant defines visual transparency of graphene Large-scale pattern growth of graphene films for stretchable transparent electrodes Controlled doping of semiconducting titania nanosheets for tailored spinelectronic materials Transforming ground mica into high-performance biomimetic polymeric mica film Low temperature reduction route to synthesise bismuth telluride (Bi2Te3) nanoparticles: structural and optical studies Lattice dynamics in mono-and few-layer sheets of WS2 and WSe2 Edge and confinement effects allow in situ measurement of size and thickness of liquid-exfoliated nanosheets Phosphorene: synthesis, scale-up, and quantitative optical spectroscopy Liquid exfoliation of solvent-stabilized few-layer black phosphorus for applications beyond electronics A mixed alkali metal titanate with the lepidocrocite-like layered structure Preparation, crystal structure, protonic form, and acid-base intercalation properties We thank Mingqiang Liu for the growth of bulk b-P crystals by short-distance transport (SDT) method,