Scalable Exfoliation of Bulk MoS2 to Single- and Few-Layers Using Toroidal Taylor Vortices

The production of a large amount of high-quality transition metal dichalcogenides is critical for their use in industrial applications. Here, we demonstrate the scalable exfoliation of bulk molybdenum disulfide (MoS2) powders into single- or few-layer nanosheets using the Taylor-Couette flow. The toroidal Taylor vortices generated in the Taylor-Couette flow provide efficient mixing and high shear stresses on the surfaces of materials, resulting in a more efficient exfoliation of the layered materials. The bulk MoS2 powders dispersed in N-methyl-2-pyrrolidone (NMP) were exfoliated with the Taylor-Couette flow by varying the process parameters, including the initial concentration of MoS2 in the NMP, rotation speed of the reactor, reaction time, and temperature. With a batch process at an optimal condition, half of the exfoliated MoS2 nanosheets were thinner than ~3 nm, corresponding to single to ~4 layers. The spectroscopic and microscopic analysis revealed that the exfoliated MoS2 nanosheets contained the same quality as the bulk powders without any contamination or modification. Furthermore, the continuous exfoliation of MoS2 was demonstrated by the Taylor-Couette flow reactor, which produced an exfoliated MoS2 solution with a concentration of ~0.102 mg/mL. This technique is a promising way for the scalable production of single- or few-layer MoS2 nanosheets without using hazardous intercalation materials.

The large-scale production of TMDs is a crucial and essential step for their realization in industrial applications. Liquid-phase exfoliation is a promising means for producing high-quality ultra-thin nanosheets with a high throughput [12]. Lithium-based intercalation has been a classical exfoliation method for MoS 2 [13]. During the liquid-phase exfoliation, two steps are involved; (i) initially, the foreign species intercalate between the adjacent layers, resulting in an increase of the interlayer spacing and a decrease of the interlayer interaction, and (ii) each layer is exfoliated by subsequent treatments using external forces such as sonication. However, the lithium-based exfoliation has several issues such as a long intercalation time, flammability of lithium in air, and tedious purifying steps, which are undesirable for an environmentally friendly and low-cost production in the industry [14]. Furthermore, the lithium intercalation induces the phase transformation of MoS 2 from semiconducting 2H-MoS 2 to metallic 1T-MoS 2 , requiring additional treatments to recover its semiconducting structure [15].
As an alternative to the lithium-based exfoliation, Coleman et al. have reported simple liquid-phase exfoliation of bulk MoS 2 powders using sonication with common solvents [16]. The exfoliated MoS 2 nanosheets had thicknesses of~3 nm to~12 nm corresponding to~4 to~17 layers, as the thickness of the monolayer MoS 2 is~0.7 nm. Although the sonication treatment is efficient to apply high energy, it may induce defects like tears and pinholes in exfoliated layers. In addition, it was difficult to obtain the monolayer MoS 2 when they used the sonication treatment. Therefore, Coleman and his colleagues utilized strong shear forces to efficiently exfoliate the MoS 2 after dispersing the MoS 2 powders in aqueous surfactant solutions [17]. Even though the work demonstrated a potential of shear exfoliation on a large scale, the exfoliation is highly dependent on the surfactant concentration and the surfactant may need to be removed for further applications.
In order to avoid the use of intercalating agents and to improve the exfoliation efficiency, we demonstrate a facile exfoliation method using toroidal Taylor vortices in the Taylor-Couette flow. The Taylor-Couette flow is generated in a small gap between two concentric cylinders; the inner cylinder rotates while the outer one is fixed. If the shaft rotates at a higher rate than a critical point, the centrifugal force exceeds the viscous force and toroidal Taylor vortices are formed. The formation of these vortices in the annulus along with the cylinders generates high shear forces with homogeneous turbulent mixing [18]. Therefore, the Taylor-Couette flow has been used for enhancing the dispersion of particles and chemical reactions [19][20][21]. Recently, it has also been utilized in exfoliating graphite into graphene or graphene oxide nanosheets [22][23][24]. In this work, we applied the Taylor-Couette flow to exfoliate bulk MoS 2 powders into single or few layers. Instead of using an organic intercalating agent, a common organic solvent was used for helping with the exfoliation and dispersion of the MoS 2 nanosheets. Furthermore, the Taylor-Couette flow reaction was adapted for a continuous exfoliation process of MoS 2 nanosheets, presenting a potential for scalable production.

Exfoliation of MoS 2 Layers with Taylor-Couette Flow
Bulk MoS 2 powders (CAS-804169, Sigma Aldrich, St. Louis, MO, USA) with lateral sizes ranging from sub micrometer to a few micrometers were initially grinded by a mortar and pestle, prior to exfoliation in the Taylor-Couette flow. A small amount of N-methyl-2-pyrrolidone (NMP, Dae-Jung Chemicals and Metals) was added into the powders and grinded for about 3 h. The shapes of the raw and grinded MoS 2 powders are shown in Figure S1 (Supplementary Materials).
The Taylor-Couette flow reactor (Laminar Co. Ltd., Seongnam, South Korea) consisted of two concentric cylinders with a small gap distance of 1.125 mm (Figure 1). The inner cylinder rotates while the outer cylinder is stationary. The gap between the inner shaft and inner wall of the outer jacket had a volume capacity of 12 mL and was surrounded with circulating hot water to control the reaction temperature. The grinded MoS 2 powders mixed with NMP were fed into the gap of the reactor and the inner shaft was rotated for a given reaction time. After completing the exfoliation process, the MoS 2 dispersed in the NMP was obtained and centrifuged at 1500 rpm for 30 min to remove the unexfoliated MoS 2 flakes. The upper half of the centrifuged solution was decanted for collecting the exfoliated MoS 2 nanosheets.

Characterization
The morphology of the exfoliated MoS2 nanosheets was observed by scanning electron microscopy (SEM, Jeol JSM-7600) after vacuum filtration with an anodic aluminum oxide (AAO) membrane. High-resolution transmission electron microscopy (TEM, Jeol ARF 200F, Tokyo, Japan) with CEOS Cs aberration correctors was used to investigate the atomic structures of the exfoliated nanosheets. The thicknesses of the exfoliated MoS2 nanosheets were characterized by atomic force microscopy (AFM, Seiko Instruments Inc., Chiba, Japan, SII-SPA-300HV) in a tapping mode after placing MoS2 nanosheets on a mechanically cleaved mica surface. The chemical structures of the exfoliated MoS2 were characterized with X-ray photoelectron microscopy (XPS, Thermo-Scientific, Waltham, MA, USA, ESCALAB-250 with a monochromated Al Kα radiation), Raman spectroscopy (Nanobase, with a 532 nm wavelength laser), and a UV-VIS spectrophotometer (Varian, Cary-5000 UV-vis-NIR, Palo Alto, CA, USA). For the XPS analysis, a survey scan was performed to identify the peaks of all of the elements, and calibration was done with respect to the C 1s peak at 285.0 eV. The peak deconvolution was done and all of the peaks were fit with Gaussian-Lorentzian functions after the subtraction of Shirley-type baselines. To estimate the concentration of the exfoliated MoS2 in a solution, the mass of the exfoliated nanosheets was measured with a microbalance, after vacuumfiltering the solution on an AAO membrane and drying at 60 °C for 5 h in a vacuum.

Exfoliation of Bulk MoS2 in NMP
As an exfoliation and dispersion solvent, NMP was selected for this work. For an efficient exfoliation of MoS2, the surface energy of the solvent needs to be matched with that of MoS2 in a thermodynamic point of view [16,25,26]. NMP has a surface tension (σ) of ~40 mJ/m 2 and its surface energy (γ) is governed by γ = σ + TSS, where T is the temperature, SS is the surface entropy (~0.1 mJ/K·m 2 ), and TSS = ~29 mJ/m 2 for liquids at room temperature [26]. Thus, the surface energy of the NMP closely matches that of the MoS2 (~75 mJ/m 2 ) [26]. However, considering only the surface energy (or surface tension) value is not sufficient for effective exfoliation, as it has been observed that many solvents having similar surface energies exhibited a worse exfoliation [16]. The surface tension can be divided into two components, polar and dispersive interactions. It has been found that the ratio of polar (σ p ) to dispersive (σ d ) components of the surface tension is more important for selecting a suitable solvent when exfoliating 2D materials [27,28]. The polar to dispersive component ratios (σ p /σ d ) for NMP and MoS2 are 0.396 and 0.449, respectively [29]. Thus, NMP is a proper solvent for the exfoliation process of MoS2.
Prior to the exfoliation process in the Taylor-Couette flow, the bulk MoS2 powders were grinded with a small amount of NMP. The grinding process provided thinner MoS2 flakes with reduced lateral sizes, helping the subsequent exfoliation of MoS2 in the Talyor-Couette flow ( Figure S1 in Supplementary Materials).
In the Taylor-Couette flow, the toroidal Taylor vortices are formed when the Taylor number (Ta) exceeds the critical value of ~1700. The Taylor number is calculated by the following [24,30]:

Characterization
The morphology of the exfoliated MoS 2 nanosheets was observed by scanning electron microscopy (SEM, Jeol JSM-7600) after vacuum filtration with an anodic aluminum oxide (AAO) membrane. High-resolution transmission electron microscopy (TEM, Jeol ARF 200F, Tokyo, Japan) with CEOS Cs aberration correctors was used to investigate the atomic structures of the exfoliated nanosheets. The thicknesses of the exfoliated MoS 2 nanosheets were characterized by atomic force microscopy (AFM, Seiko Instruments Inc., Chiba, Japan, SII-SPA-300HV) in a tapping mode after placing MoS 2 nanosheets on a mechanically cleaved mica surface. The chemical structures of the exfoliated MoS 2 were characterized with X-ray photoelectron microscopy (XPS, Thermo-Scientific, Waltham, MA, USA, ESCALAB-250 with a monochromated Al K α radiation), Raman spectroscopy (Nanobase, with a 532 nm wavelength laser), and a UV-VIS spectrophotometer (Varian, Cary-5000 UV-vis-NIR, Palo Alto, CA, USA). For the XPS analysis, a survey scan was performed to identify the peaks of all of the elements, and calibration was done with respect to the C 1s peak at 285.0 eV. The peak deconvolution was done and all of the peaks were fit with Gaussian-Lorentzian functions after the subtraction of Shirley-type baselines. To estimate the concentration of the exfoliated MoS 2 in a solution, the mass of the exfoliated nanosheets was measured with a microbalance, after vacuum-filtering the solution on an AAO membrane and drying at 60 • C for 5 h in a vacuum.

Exfoliation of Bulk MoS 2 in NMP
As an exfoliation and dispersion solvent, NMP was selected for this work. For an efficient exfoliation of MoS 2 , the surface energy of the solvent needs to be matched with that of MoS 2 in a thermodynamic point of view [16,25,26]. NMP has a surface tension (σ) of~40 mJ/m 2 and its surface energy (γ) is governed by γ = σ + TS S , where T is the temperature, S S is the surface entropy (~0.1 mJ/K·m 2 ), and TS S =~29 mJ/m 2 for liquids at room temperature [26]. Thus, the surface energy of the NMP closely matches that of the MoS 2 (~75 mJ/m 2 ) [26]. However, considering only the surface energy (or surface tension) value is not sufficient for effective exfoliation, as it has been observed that many solvents having similar surface energies exhibited a worse exfoliation [16]. The surface tension can be divided into two components, polar and dispersive interactions. It has been found that the ratio of polar (σ p ) to dispersive (σ d ) components of the surface tension is more important for selecting a suitable solvent when exfoliating 2D materials [27,28]. The polar to dispersive component ratios (σ p /σ d ) for NMP and MoS 2 are 0.396 and 0.449, respectively [29]. Thus, NMP is a proper solvent for the exfoliation process of MoS 2 .
Prior to the exfoliation process in the Taylor-Couette flow, the bulk MoS 2 powders were grinded with a small amount of NMP. The grinding process provided thinner MoS 2 flakes with reduced lateral sizes, helping the subsequent exfoliation of MoS 2 in the Talyor-Couette flow ( Figure S1 in Supplementary Materials).
In the Taylor-Couette flow, the toroidal Taylor vortices are formed when the Taylor number (T a ) exceeds the critical value of~1700. The Taylor number is calculated by the following [24,30]: where Ω is the angular velocity, ν is the kinematic viscosity, and R o and R i are the radii of the outer and inner cylinders, respectively. Based on the equation, the Taylor number is above the critical value even at a low rotation speed in this work, T a =~1.25 × 10 4 for 500 rpm, which confirms the formation of the Taylor vortices over the whole experiments. In addition, it has been known that the average wall shear rate in the Taylor-Couette flow increases with the increasing the angular velocity [31]. Thus, increasing the rotational speed of the inner cylinder in the Taylor-Couette flow reactor is expected to enhance the exfoliation efficiency. The wall shear stress in the Taylor-Couette flow can be estimated by using the computational fluid dynamics (CFD) simulation. The previous work reported that the maximum shear stress of the Taylor-Couette flow in NMP increased from 3.98 Pa at a rotation speed of 500 rpm to 26.91 Pa at 3000 rpm, where the reactor had a gap size of 2.5 mm and an inner shaft diameter of 52 mm [24]. Even though the reactor dimensions were different from those in our work, the CFD results showed that the shear stress in the flow is highly dependent on the rotation speed.

Parametric Study for Optimizing the Exfoliation Process
Several factors can affect the exfoliation process with the Taylor-Couette flow. To find an optimal exfoliation condition for a higher concentration of exfoliated MoS 2 , four parameters were investigated, namely: the initial concentration of MoS 2 in a solution, the reaction time, the rotation speed of the inner cylinder, and the reaction temperature. Firstly, the initial MoS 2 concentration in NMP varied from 1 mg/mL to 50 mg/mL at a rotation speed of 3000 rpm, a reaction temperature of 90 • C, and a reaction time of 120 min. It is found that the maximum concentration of the exfoliated MoS 2 was 0.0419 mg/mL at the initial MoS 2 concentration of 10 mg/mL ( Figure 2a). However, further increasing the initial MoS 2 concentration leads an adverse effect on the exfoliated quantity. It might be due to the increased viscosity of the MoS 2 solution with more MoS 2 , resulting in lower shear forces in the Taylor-Couette flow. Thus, the optimum condition for the initial MoS 2 solution concentration was 10 mg/mL, and we kept this for the other exfoliation experiments.
As discussed above, the local shear forces are highly dependent on the rotation speed in the Taylor-Couette flow. We varied the rotation speed of the inner cylinder from 500 rpm to 3000 rpm at a fixed reaction time of 120 min and a reaction temperature of 90 • C. It can be clearly noticed from Figure 2b that the concentration of the exfoliated MoS 2 increases when increasing the rotation speed. The results confirm the presence of a higher shear rate and wall shear stress at a higher rotation speed in the Taylor-Couette flow, which is critical for the efficient exfoliation of MoS 2 . Therefore, we have chosen the rotation speed of 3000 rpm (the maximum speed of the reactor) for all of the experiments in this study.
Another important experimental parameter is the reaction time, because other existing exfoliation methods consume too much time to exfoliate the layers of their bulk materials. Thus, a method that can save a process time with efficient exfoliation results would be a promising choice for the scalable production of few-layer MoS 2 . Figure 2c shows the concentration change of the exfoliated MoS 2 by varying the reaction time from 15 min to 180 min, while fixing the rotation speed (3000 rpm) and the reaction temperature (90 • C). It is observed that the concentration of the exfoliated MoS 2 increases with the increasing the reaction time.
reaction time of 120 min. It is found that the maximum concentration of the exfoliated MoS2 was ~0.0419 mg/mL at the initial MoS2 concentration of 10 mg/mL (Figure 2a). However, further increasing the initial MoS2 concentration leads an adverse effect on the exfoliated quantity. It might be due to the increased viscosity of the MoS2 solution with more MoS2, resulting in lower shear forces in the Taylor-Couette flow. Thus, the optimum condition for the initial MoS2 solution concentration was 10 mg/mL, and we kept this for the other exfoliation experiments.  Also, the reaction time was varied from 30 • C to 90 • C by circulating water through the outside cylinder, while fixing the rotation speed (3000 rpm) and the reaction time (120 min). We observed that the concentration of the exfoliated MoS 2 nanosheets increased when increasing the reaction temperature, as shown in Figure 2d. Generally, the viscosity of the liquid decreases with an increase in temperature. The decrease in the kinematic viscosity indicates the increase of the Taylor number, meaning an increased shear rate. Therefore, the elevated temperature helped the exfoliation of MoS 2 in the Taylor-Couette flow. Furthermore, this may be partially due to the auto-oxidation of NMP at an elevated temperature [29]. It is known that NMP can be automatically oxidized at an elevated temperature with the presence of oxygen and water from the atmosphere. This leads to the formation of active intermediate species, such as radicals and hydroperoxide, which may oxidize the MoS 2 edges and help the exfoliation of MoS 2 [29]. Figure 3a shows a SEM image of the MoS 2 nanosheets exfoliated at the optimized reaction condition. The exfoliated nanosheets were vacuum-filtered and placed on an AAO filter membrane. It can be clearly seen that the nanosheets are relatively thin compared to original MoS 2 powders ( Figure S1 in Supplementary Materials). The nanosheets are thin enough to see through the pore structures of the AAO membrane. The atomic structures of the exfoliated MoS 2 were analyzed by TEM, as shown in Figure 3b,c. The nanosheets that appeared to be well exfoliated had single or few layers. Figure 3b shows a monolayer MoS 2 nanosheet with a typical hexagonal pattern in a selected area electron diffraction (SAED). The high-resolution TEM images show the high-quality crystalline structures of MoS 2 (Figure 3c), confirming that the Taylor-Couette flow reaction did not damage or destroy the intrinsic hexagonal structure of MoS 2 . structures of the AAO membrane. The atomic structures of the exfoliated MoS2 were analyzed by TEM, as shown in Figure 3b,c. The nanosheets that appeared to be well exfoliated had single or few layers. Figure 3b shows a monolayer MoS2 nanosheet with a typical hexagonal pattern in a selected area electron diffraction (SAED). The high-resolution TEM images show the high-quality crystalline structures of MoS2 (Figure 3c), confirming that the Taylor-Couette flow reaction did not damage or destroy the intrinsic hexagonal structure of MoS2.  The size and thickness of the exfoliated MoS 2 nanosheets were investigated by AFM measurements, as shown in Figure 4. The mechanically cleaved mica was used as a substrate for the AFM samples, because of its atomically flat surface. Figure 4a shows a morphology image of the deposited MoS 2 nanosheets; a scan in a smaller area is shown in Figure S2 (Supplementary Materials). Figure 4b shows the histograms for the corresponding thickness and lateral size of the sample. More than 130 nanosheets were measured and about 50% of the measured nanosheets were thinner than~3 nm, indicating that the exfoliated MoS 2 nanosheets were predominately thinner than~4 layers, by assuming that monolayer thickness is~0.7 nm. However, the exfoliated nanosheets show smaller lateral sizes compared to the electrochemically exfoliated MoS 2 nanosheets [32]. The size and thickness of the exfoliated MoS2 nanosheets were investigated by AFM measurements, as shown in Figure 4. The mechanically cleaved mica was used as a substrate for the AFM samples, because of its atomically flat surface. Figure 4a shows a morphology image of the deposited MoS2 nanosheets; a scan in a smaller area is shown in Figure S2 (Supplementary Materials). Figure 4b shows the histograms for the corresponding thickness and lateral size of the sample. More than 130 nanosheets were measured and about 50% of the measured nanosheets were thinner than ~3 nm, indicating that the exfoliated MoS2 nanosheets were predominately thinner than ~4 layers, by assuming that monolayer thickness is ~0.7 nm. However, the exfoliated nanosheets show smaller lateral sizes compared to the electrochemically exfoliated MoS2 nanosheets [32].  Figure 5a shows the typical Raman spectra of the bulk and exfoliated MoS2 nanosheets. Two characteristic peaks are found at 384 (E 1 2g) and 410 cm −1 (A1g) for the bulk MoS2. The E 1 2g and A1g are known to be related to in-plane opposite vibrations of sulfur/molybdenum atoms and out-of-plane vibrations of sulfur atoms, respectively [32,33]. As the bulk MoS2 was exfoliated to nanosheets, the E 1 2g and A1g peaks shifted to 382 and 406 cm −1 , and the frequency difference between the two peaks decreased from 26 cm −1 to 24 cm −1 . The change of positions and the frequency difference of the E 1 2g and A1g peaks are attributed to the change of the MoS2 layer numbers [33]. The frequency difference of 24 cm −1 of the exfoliated MoS2 nanosheets corresponds to roughly 4-5 layers [33].  and A 1g are known to be related to in-plane opposite vibrations of sulfur/molybdenum atoms and out-of-plane vibrations of sulfur atoms, respectively [32,33]. As the bulk MoS 2 was exfoliated to nanosheets, the E 1 2g and A 1g peaks shifted to 382 and 406 cm −1 , and the frequency difference between the two peaks decreased from 26 cm −1 to 24 cm −1 . The change of positions and the frequency difference of the E 1 2g and A 1g peaks are attributed to the change of the MoS 2 layer numbers [33]. The frequency difference of 24 cm −1 of the exfoliated MoS 2 nanosheets corresponds to roughly 4-5 layers [33]. Figure 5a shows the typical Raman spectra of the bulk and exfoliated MoS2 nanosheets. Two characteristic peaks are found at 384 (E 1 2g) and 410 cm −1 (A1g) for the bulk MoS2. The E 1 2g and A1g are known to be related to in-plane opposite vibrations of sulfur/molybdenum atoms and out-of-plane vibrations of sulfur atoms, respectively [32,33]. As the bulk MoS2 was exfoliated to nanosheets, the E 1 2g and A1g peaks shifted to 382 and 406 cm −1 , and the frequency difference between the two peaks decreased from 26 cm −1 to 24 cm −1 . The change of positions and the frequency difference of the E 1 2g and A1g peaks are attributed to the change of the MoS2 layer numbers [33]. The frequency difference of 24 cm −1 of the exfoliated MoS2 nanosheets corresponds to roughly 4-5 layers [33].   [34]. Therefore, the reaction in the Talyor-Couette flow exfoliated the bulk MoS 2 and decreased its layer numbers.

Characterization of Exfoliated MoS 2 Nanosheets
The chemical nature of the exfoliated MoS 2 nanosheets was investigated using X-ray photoelectron spectroscopy (XPS). Figure 6 shows the XPS spectra of the experimental and fitted data for the bulk and exfoliated MoS 2 samples. For the exfoliated MoS 2 nanosheets, there are peaks for Mo 3d 5/2 and Mo 3d 3/2 at 229.4 and 232.6 eV, respectively, and additional peaks for S 2p 3/2 and S 2p 1/2 at 162.2 and 163.4 eV, respectively [32]. For the bulk MoS 2 , the peaks for Mo 3d 5/2 , Mo 3d 3/2 , S 2p 3/2 , and S 2p 1/2 are positioned at 229.1, 232.3, 161.9, and 163.0 eV, respectively. The peaks for Mo and S in the XPS spectra of the exfoliated nanosheets are very similar to those of the bulk powders, indicating that the exfoliation process in the Taylor-Couette flow did not alter the original chemical structures of the bulk MoS 2 . In addition, the peak positions and shapes for the Mo and S of the exfoliated nanosheets are the same as those of 2H-MoS 2 nanosheets, which were obtained by chemical exfoliation and subsequent thermal annealing [35]. This confirms that the exfoliated nanosheets kept semiconducting the 2H-MoS 2 phase after the exfoliation process.

Continuous Production of Exfoliated MoS 2
In order to achieve a high-throughput production of the few-layer MoS 2 nanosheets, the Taylor-Couette flow reactor was adapted for a continuous flow process. Figure 7a shows the schematic illustration of the continuous flow method. In this method, a MoS 2 solution was continuously fed into the cylinder gap through the bottom of the reactor using a peristaltic pump with a constant flow rate. The inner cylinder was continuously rotating during the reaction and the excess of the solution was collected in the flask and again recirculated inside the reactor. In the whole process, the MoS 2 dispersion was continuously sheared between the rotating inner and the stationary outer shaft of the reactor. The optimal conditions (initial MoS 2 concentration = 10 mg/mL, rotation speed = 3000 rpm, and reaction temperature = 90 • C) were used for this continuous exfoliation.
XPS spectra of the exfoliated nanosheets are very similar to those of the bulk powders, indicating that the exfoliation process in the Taylor-Couette flow did not alter the original chemical structures of the bulk MoS2. In addition, the peak positions and shapes for the Mo and S of the exfoliated nanosheets are the same as those of 2H-MoS2 nanosheets, which were obtained by chemical exfoliation and subsequent thermal annealing [35]. This confirms that the exfoliated nanosheets kept semiconducting the 2H-MoS2 phase after the exfoliation process.

Continuous Production of Exfoliated MoS2
In order to achieve a high-throughput production of the few-layer MoS2 nanosheets, the Taylor-Couette flow reactor was adapted for a continuous flow process. Figure 7a shows the schematic illustration of the continuous flow method. In this method, a MoS2 solution was continuously fed into the cylinder gap through the bottom of the reactor using a peristaltic pump with a constant flow rate. The inner cylinder was continuously rotating during the reaction and the excess of the solution was collected in the flask and again recirculated inside the reactor. In the whole process, the MoS2 dispersion was continuously sheared between the rotating inner and the stationary outer shaft of the reactor. The optimal conditions (initial MoS2 concentration = 10 mg/mL, rotation speed = 3000 rpm, and reaction temperature = 90 °C) were used for this continuous exfoliation. Figure 7b,c shows the change of the concentration of the exfoliated MoS2 nanosheets after the reactions for 2 h. The slow feeding of the MoS2 solution generated a higher concentration of the exfoliated MoS2 solution. Interestingly, the maximum concentration of the exfoliated nanosheets was 0.102 mg/mL at the lowest feeding rate of 0.09 mL/min, which was higher compared to the nanosheets exfoliated in a batch process (0.0419 mg/mL). In this work, the flow reactor is vertically positioned. Thus, the bulk MoS2 powders could sink on the bottom of the reactor by gravity when they were exfoliated in a batch process. For the continuous exfoliation, the liquid could push up the bulk MoS2 powders and might induce a better dispersion and more shear interactions during the process, resulting in more exfoliation. However, with increasing the feeding rate, the MoS2 particles might pass through the Taylor vortices with reduced shear interactions, providing a low exfoliation efficiency.
There have been several reports on the high-yield exfoliation of MoS2. For example, lithium intercalation and the subsequent sonication treatments generated thin MoS2 layers in a 92% yield [36]. However, the use of lithium may not be desirable for the environmentally friendly and low-cost production of MoS2 nanosheets because of several issues related to lithium, such as long intercalation time, flammability in air, and purification. Thus, a simple method using shear exfoliation in common solvents has been demonstrated for graphene with a concentration of ~0.3 mg/mL in a below 0.1% yield [37]. Our work demonstrated MoS2 exfoliation into single and few layers in reasonably high yield compared to the shear exfoliation of graphite in NMP. Thus, this continuous method using the Taylor-Couette flow leads to a more productive and stable exfoliation process for MoS2 and other 2D materials.

Conclusions
In this study, a facile and efficient method was presented for the exfoliation of bulk MoS2 layered crystals into a few layers. The toroidal Taylor vortices in the Taylor-Couette flow generated strong shear stress on MoS2, resulting in an efficient exfoliation in NMP. We found that half of the exfoliated nanosheets were thinner than ~3 nm (about ~4 layers) and the atomic and chemical characteristics were not altered by the exfoliation process. In addition, a continuous production system using the Taylor-Couette flow was demonstrated for exfoliating MoS2. The maximum concentration of the exfoliated nanosheets was ~0.102 mg/mL for the 2-h process. This method has potential for scalable  Interestingly, the maximum concentration of the exfoliated nanosheets was 0.102 mg/mL at the lowest feeding rate of 0.09 mL/min, which was higher compared to the nanosheets exfoliated in a batch process (0.0419 mg/mL). In this work, the flow reactor is vertically positioned. Thus, the bulk MoS 2 powders could sink on the bottom of the reactor by gravity when they were exfoliated in a batch process. For the continuous exfoliation, the liquid could push up the bulk MoS 2 powders and might induce a better dispersion and more shear interactions during the process, resulting in more exfoliation. However, with increasing the feeding rate, the MoS 2 particles might pass through the Taylor vortices with reduced shear interactions, providing a low exfoliation efficiency.
There have been several reports on the high-yield exfoliation of MoS 2 . For example, lithium intercalation and the subsequent sonication treatments generated thin MoS 2 layers in a 92% yield [36]. However, the use of lithium may not be desirable for the environmentally friendly and low-cost production of MoS 2 nanosheets because of several issues related to lithium, such as long intercalation time, flammability in air, and purification. Thus, a simple method using shear exfoliation in common solvents has been demonstrated for graphene with a concentration of~0.3 mg/mL in a below 0.1% yield [37]. Our work demonstrated MoS 2 exfoliation into single and few layers in reasonably high yield compared to the shear exfoliation of graphite in NMP. Thus, this continuous method using the Taylor-Couette flow leads to a more productive and stable exfoliation process for MoS 2 and other 2D materials.

Conclusions
In this study, a facile and efficient method was presented for the exfoliation of bulk MoS 2 layered crystals into a few layers. The toroidal Taylor vortices in the Taylor-Couette flow generated strong shear stress on MoS 2 , resulting in an efficient exfoliation in NMP. We found that half of the exfoliated nanosheets were thinner than~3 nm (about~4 layers) and the atomic and chemical characteristics were not altered by the exfoliation process. In addition, a continuous production system using the Taylor-Couette flow was demonstrated for exfoliating MoS 2 . The maximum concentration of the exfoliated nanosheets was~0.102 mg/mL for the 2-h process. This method has potential for scalable production of few-layer MoS 2 nanosheets, which is critical for various device applications in the future.

Conflicts of Interest:
The authors declare no conflict of interest.