Research Progress of the Liquid-Phase Exfoliation and Stable Dispersion Mechanism and Method of Graphene

Graphene is a honeycomb hexagonal 2D crystal nanomaterial with a thickness of only 0.334 nm. It has been widely used and studied because of its ultra-thin 2D nano-characteristics and excellent electrical, thermal, optical and mechanical properties. With the continuous in-depth study of graphene, the liquid phase dispersion of graphene has also achieved breakthroughs. This review summarises the research progress of the liquid phase exfoliation mechanism, exfoliation method, stable dispersion mechanism and dispersion method of graphene in recent years. The research situations of Hamaker constant theory in exfoliation mechanism and Hansen solubility coefficient theory in stable dispersion mechanism are mainly discussed. The shortcomings of the research are summarised, and analyses the graphene liquid phase dispersed in the important challenge in the future. The stable dispersion method of graphene is also summarised. In the future, the π-π interaction is the most potential method for studying graphene stabilization.


INTRODUCTION
Materials are the basis and forerunner for the development of human society (Du, 2007;Wegst et al., 2015), and the birth of new materials is an important milestone of social progress. New materials technology of modern industries supports human civilization. Carbon materials are one of the richest resources on earth and play an important role in the history of human development (Cheng, 1998). After the discovery of fullerenes in 1985 and carbon nanotubes (CNTs) in 1991, the existence of materials with zero-dimensional (Fuller sphere), one-dimensional (1D) (CNTs), and three-dimensional (3D) (graphite) carbon has been confirmed by experiments (Cheng et al., 2006;Zhang, 2018). In 2004, Geim and Novoselov of Manchester University adopted the tape stripping method to prepare high-quality single-layer graphene for the first time (Novoselov et al., 2004), overturning the recognition that the finite temperature, classical thermal theory does not allow the perception of the existence of two-dimensional (2D) crystals (Xu et al., 2009). Graphene is a 2D nanomaterial composed of a single layer of carbon atoms and is closely mixed with sp 2 (Si and Samulski, 2008;Zhu et al., 2013). Its surface has a large area of π-π conjugated bond and has excellent electrical properties (at room temperature, electron mobility: 15,000 cm 2 [V·s] ) (Allen et al., 2010;Green and Hersam, 2010;Brownson, 2013), thermal properties (thermal conductivity: 5,300 W [m·K] , which is 10 times higher than that of copper and aluminum) (Sadasivuni et al., 2014;Xu et al., 2016), mechanical properties (tensile strength: 130 GPa, which is 100 times that of steel) (Sadasivuni et al., 2014;Long et al., 2016), optical properties (single-layer transmittance: 97%, 2.3% for each additional layer absorbance) (Zhang et al., 2009;Cao and Guo, 2016), and a large specific surface area (2,630 m 2 /g) (Kim et al., 2010;Zou et al., 2014). Given graphene's excellent 2D nano-features, the modification of SnO 2 with highly conductive 2D naphthalene diimide graphene can increase surface hydrophobicity and form van der Waals interaction between the surfactant and the organic-inorganic hybrid lead halide perovskite compounds. As a result, highly efficient perovskite solar cells with power conversion efficiency of 20.2% can be achieved with an improved fill factor of 82%, which could be mainly attributed to the augmented charge extraction and transport (Zhao X. J. et al., 2018). In anti-corrosion coating applications, graphene can block oxygen molecules and corrosion factors from penetrating into the coating substrate and improve the salt spray resistance of the coating by 20% and the antibacterial property by 50% (Gu et al., 2015(Gu et al., , 2016Wang et al., 2016;Zhao et al., 2017). As lithium-ion battery electrode materials, graphene and lithium iron phosphate are mixed, which increases the reversible specific capacity of the battery by 1.4 times Uchaker and Cao, 2014). As a lubricating oil additive, the 2D layered structure of graphene has a friction coefficient reduced by up to 30% compared with ordinary lubricating oil through the interlayer slip effect Liu et al., 2019). However, the van der Waals force between graphene sheets is prone to agglomeration (Li and Shi, 2014;Xie et al., 2014). As a result, it is neither lipophilic nor hydrophilic (Li et al., 2008), and it is difficult to uniformly disperse in various liquid media in practical applications (Vallés et al., 2008;Wan et al., 2013;Cui and Zhou, 2018). Obviously, the liquid-phase dispersion of graphene is the core problem that must be overcome in the experimental study of nanomaterials Coletti et al., 2013;Wang and Chen, 2015;Sahoo and Ramaprabhu, 2017;Yu et al., 2018) to promote its industrial production (Yang K. et al., 2010;Secor et al., 2013;Wang and Chen, 2015;Fadavi et al., 2016;Huan-Yan et al., 2018).
Graphene liquid-phase dispersion can be divided into two parts: graphene exfoliation (because the degree of exfoliation of graphene determines its dispersion grade) and graphene stable dispersion. The exfoliation theory is the major discipline problem involving molecular inter-atomic forces, and its final goal is to overcome the graphene sheet interlayer van der Waals force and graphene layer interaction between functional groups. Blocking the adjacent layers of graphene to each other achieves monolayer of graphene. To quantify the graphene interlayer van der Waals force, researchers established the Hamaker constant (Diba et al., 2016) theory on the basis of intermolecular forces. The stable dispersion theory major discipline involved in the problem is the interface between the molecules of the solution theory. That is, an appropriate dispersant-solvent system that interacts with the graphene surface must be selected to improve the liquid medium interface properties and realize the uniform and stable dispersion of graphene. To cope with graphene matching an appropriate dispersant-solvent system, researchers studied the potential energy of the graphene surface and its relationship with the potential energy of the dispersant-solvent system, eventually believing the stability of the graphene dispersion and dispersant-solvent system of Hansen solubility coefficient (HSP) (Lim et al., 2014).
In this paper, the liquid-phase dispersion theory of graphene is reviewed from the aspects of liquid-phase exfoliation and stable dispersion. Firstly, aiming at overcoming the van der Waals force between graphene sheets, the progress of graphene liquidphase exfoliation theory in recent years was summarized. The research results of the Hamaker constant theory were reviewed in each research group, and the advantages and disadvantages of the main exfoliation methods (e.g., ultrasonic degradation, mechanical exfoliation, and electrochemical exfoliation) were reviewed and analyzed. Then, in the mechanism of stable dispersion of graphene, the theory of the HSP was discussed, which was the basis for selecting a suitable dispersant-solvent system to interact with graphene. Finally, the research progress of stable dispersion methods of graphene was summarized, especially the π-π interaction in the stable dispersion methods of graphene.

Mechanism of Graphene Liquid-Phase Exfoliation
Currently, compared with the low-yield solid phase (Novoselov et al., 2004;Gao et al., 2018;Melios et al., 2018) (e.g., micromechanical exfoliation and epitaxial growth on SiC) method and the high-energy-consuming gas phase method (e.g., chemical vapor deposition) Mao et al., 2018;Xia et al., 2018;Cheng et al., 2019), the liquid-phase exfoliation method has the advantages of higher processability, lower cost (Artur and Paolo, 2016), shorter preparation period, and largerscale preparation (Artur and Paolo, 2014;Paton et al., 2014;Yang et al., 2015;Zhang et al., 2017). Liquid-phase exfoliation of graphene is a process in which graphite (or multilayer graphene) is converted into a single layer (or a few layers) of graphene by overcoming the van der Waals force between sheets (Allen et al., 2010). The van der Waals force between graphene sheets is mainly composed of Keesom interaction, Debye interaction, and dispersive force (Alecrim et al., 2015), where the dispersive force is the most important interaction force that constitutes the van der Waals force between atoms and molecules, and the intensity is generally less than the Coulomb force or hydrogen bonding (Israelachvili, 2011). With different exfoliation methods to overcome the van der Waals force between layers, the dispersant-solvent system enters the graphene sheet layer and interacts with its surface to achieve stable dispersion of the liquid phase (Fu and Yang, 2013). Therefore, the degree of dispersion determines the dispersion level. In 2007, Huang (2007) quantified the van der Waals forces  between nanoparticles for the first time by using the formula F A = AD 24H 2 , where A is the Hamaker constant (Diba et al., 2016) of the nanoparticle in vacuum, H is the particle spacing, and D is the particle diameter. Coleman (2009) found that CNTs attract each other because of the van der Waals force interaction between adjacent carbon atoms, in which the potential energy of these paired interactions meets V r = − C r 6 (where r is the atomic spacing and C is a constant: A = π 2 ρ 2 C, where ρ is the number of atoms in the unit volume in the cylinder) and A is the Hamaker constant, which depends on the total interaction force between adjacent CNTs and also the sum of the paired interaction forces of all carbon atoms. Coleman (2012) found that the theory of graphene exfoliation was similar to that of CNT exfoliation. These potential energies are larger than the paired interacting atoms and generally interact with molecules. The dispersantsolvent/graphene system involved in the liquid-phase dispersion of olefins also has an above potential energy, and the theoretical basis for achieving the successful exfoliation of graphene is to offset the above attraction. Coleman (2012) believed that the essence of graphene exfoliation is that the Hamaker constant was calculated to accurately quantify the van der Waals forces between graphene sheets. Although the researchers first used the developed parametric model and force-distance relationship to calculate the Hamaker constant using atomic force microscopy (AFM) analysis (Chu et al., 2017;Chiou et al., 2018), Gou et al. (2016) believed that the theoretical Hamaker constant of colloid should be determined by dynamic light scattering. In addition, for the ionic graphene dispersion system, when calculating the Hamaker constant by dynamic light scattering technique, the Hofmeister effect must be considered to obtain the accurate Hamaker constant, because different types of salt ions determine different Hofmeister effects Luo et al., 2018). It has certain advantages, but a standard method for determining the Hamaker constant is currently lacking. Therefore, the Hamaker constant theory needs to be further studied in the future.

Graphene Liquid-Phase Exfoliation Method
The main liquid-phase exfoliation methods include ultrasonic degradation (Narayan et al., 2017;, mechanical exfoliation (Su et al., 2011;Jung et al., 2015;Ping et al., 2017), and electrochemical exfoliation (Parvez et al., 2014;Jung et al., 2015;Shinde et al., 2016). In the study of the liquidphase exfoliation of graphene, the most widely used method is to use natural flake graphite as raw material. In the dispersantsolvent system, graphite is first intercalated to obtain expanded graphite, and then different methods are employed to overcome the van der Waals force between graphite layers according to the Hamaker constant theory to achieve the exfoliation and dispersion of graphene (Dappe et al., 2006;Boström and Sernelius, 2012;Davidson et al., 2014;Lee et al., 2015). The dispersion mechanism and process are shown in Figure 1.

Ultrasonic Degradation
Ultrasonic degradation has the advantages of low cost and simple operation , but long-term usage of ultrasonic graphene causes defects, and layer number renders the lamella radial size difficult to control (Khan et al., 2010) and enlarges the scale and the ultrasonic generator (Tian et al., 2017). The traditional ultrasonic degradation process is no longer adopted in present experimental studies and industrial production, and many studies aimed for improvement and innovation (Moorthy et al., 2015;Abo-Zahhad et al., 2016;Narayan et al., 2017;Gai et al., 2018). Gao et al. (2017) combined an ultrasonic generator with a supercritical reactor (Figure 2) by using ultrasonic degradation and supercritical CO 2 (scCO 2 ) excellent penetration of air to natural graphite in the scCO 2 , H 2 O, and ethanol system. High-quality graphene (>50%) was obtained in the mixed system. The method for industrial production provides a green, fast, and scalable way of graphene dispersion. Phiri et al. (2017) improved the traditional process of graphene preparation by ultrasonically exfoliating graphite in water medium. In the process of ultrasonic exfoliation, graphene was exfoliated and dispersed by means of a high shear colloidal mixer for 2 h to obtain a graphene dispersion with a concentration of up to 1.1 mg/ml. The graphene fractions prepared by ultrasonic peeling graphite for 100 h were separated by traditional process. Graphene dispersions of the same quality can be prepared in 0.5 h with the improved process. The graphene prepared by the improved process is analyzed by AFM and Raman. The dispersions contain plenty of defectfree and unmodified graphene layers. Hadi et al. (2018) mixed magnetic nanoparticles (Fe 3 O 4 ) and graphite evenly during the preparation of graphene dispersion through the liquid-phase peeling of graphite. During ultrasonic exfoliation, the intense shear force between graphite nanosheets promotes the exfoliation of graphene, greatly shortens the ultrasonic time, and avoids long-term ultrasound bringing defects to graphene surface.

Mechanical Exfoliation
Because of the different parameters of the ultrasonic exfoliation process (e.g., time, frequency, power, and shape of the inner cavity), the experimental results are greatly affected (Río et al., 2017). At present, in the mechanical peeling process (including shearing and ball milling) (Deng et al., 2016;Liu Y. J. et al., 2018), the combination of shear and ball milling exfoliation process has been identified as the trend of replacing the ultrasonic peeling process in the future (Tour, 2014;Bonaccorso et al., 2016). Novoselov et al. (2004) first prepared graphene by the scotch tape stripping method and won the 2010 Nobel Prize in Physics. But the efficiency of this method is extremely limited to laboratory research (Yi and Shen, 2015). The advantage of mechanical exfoliation is that high-quality graphene can be obtained, but the mechanical exfoliation theory is more complicated than other stripping methods, and the stripping process is affected by many uncontrollable factors (e.g., the thickness of the graphene, and horizontal size). Chen et al. (2012) used the mechanical peeling method in the experimental study. In the N-methyl pyrrolidone (NMP) system, the shear force produced by high-speed mechanical exfoliation interacts with van der Waals force between graphite sheets in an inclined 45 • test tube (Figure 3a), resulting in relative slippage between graphite sheets, good peeling effect of the dispersion system assisted by NMP and overall dispersion. The process does not introduce defects in graphene (Figures 3b-d). Paton et al.

FIGURE 2 | (A)
The scCO 2 , aqueous solution, and graphite particle state are distributed before and after the ultrasonic wave is applied to the system. (B) Pilot plant for the large-scale production of graphene in scCO 2 /polyvinylpyrrolidone (PVP)/H 2 O systems by ultrasonic-assisted stripping of natural graphite (where UB is the ultrasonic bath, TC is the temperature controller, PI is the pressure indicator, and PC is pressure controller; BPR is the back-pressure regulator and CV is the check valve; x is the raw material supply system, y is the CO 2 feeding system, z is the heating/mixing system, { is the ultrasonic coupling supercritical reactor, | is the product collector, and } is the exhaust gas collector).
(2014) extended the theory of Chen et al. (2012) and proposed the theory of shear rate: where γ min is the shear rate, E S,G and E S,L are the surface energies of graphene and solvent, the surface energy of NMP is 69 mJ m 2 , η is the viscosity of the liquid, and L is the length of the graphite sheet. L = 300-800 nm of graphene was measured by transmission electron microscopy (TEM), and the minimum shear rate γ min ≈ 10 4 s −1 was calculated. E S,G was calculated by Formula (1): E S,G = 70.5-71 mJ m 2 , which is very close to the surface potential energy (E S,G = 69 mJ m 2 ) for the automatic exfoliation of graphene liquid phase. To verify the surface energy theory of graphene liquid-phase auto-exfoliation, many researchers have obtained the best surface energy through experiments, but the numerical values are different, and a unified theory is currently lacking (Kozbial et al., 2014). Single surface energy theory cannot completely explain the mechanism of graphene liquid-phase auto-exfoliation (Coleman, 2012) and thus has been replaced by the Hamaker constant theory.

Electrochemical Exfoliation
Electrochemical exfoliation method has become one of the most researched potentials of graphene liquid-phase exfoliation theory (Parvez et al., 2014) because of its advantages of being fast and has high efficiency, environmental protection, and low cost, etc. (Parvez et al., 2014;Artur and Paolo, 2016;Li et al., 2016;Ping et al., 2017). Su et al. (2011) exfoliated graphite by electrochemical stripping to prepare graphene dispersion, characterized by AFM, TEM, scanning tunneling microscopy (STM), Raman and attenuated total reflection-Fourier-transform infrared spectroscopy (ATR-FTIR). More than 60% of graphene in the dispersion was AB stack. Bilayer graphene has a lateral dimension of 30 µm on the surface, and its electrical properties were superior to those of reduced graphene oxide (rGO). Shinde et al. (2016) used hydrodynamics to prepare graphene by the shear-assisted electrochemical exfoliation method for the first time in the study of the electrochemical exfoliation of graphene. The graphene prepared at a high potential was thick and contained cracks. At low voltage (<1 V) potential, shear force had a more obvious assistant effect on electrochemical exfoliation (the exfoliation mechanism is shown in Figure 4). The monolayer, bilayer, and trilayer in the dispersion graphene occupied a high proportion. The Raman characterization showed only a few defects on the surface of graphene (I D /I G = 0.21-0.32), and these defects were distributed at the edge of graphene. Shi et al. (2018) demonstrated a non-electrified electrochemical exfoliation method to produce high-quality graphene (NEEG) sheets. By direct electrochemical reaction between graphite powders and metallic Li in 1 M of LiPF 6 /propylene carbonate electrolyte, continuous graphite exfoliation is achieved with a high yield of excess 80% without any consumption of electric energy. The as-prepared NEEG has high quality with few defects (I D /I G = 0.45), a high C/O ratio of 27.74, and excellent electronic conductivity of 102.5 S/cm. As conductive additives in Li 4 Ti 5 O 12 electrode for lithium ion batteries, NEEG contributes to higher capacity and initial coulombic efficiency than common thermally rGO. Therefore, electrochemical exfoliation method is promising for large-scale production of high-quality graphene nanosheets .

Liquid-Phase Stable Dispersion Mechanism of Graphene
The ultimate aim of the study on the stable dispersion of graphene in the liquid phase is to find the most suitable FIGURE 4 | Schematic of the experimental principle of an electrochemical microreactor. Graphite crystals were both reactor walls and working electrodes that apply an electric field at the same time as high shear. H, Q, and γ are the height between the two electrodes, the electrolyte flow rate, and the shear rate, respectively. dispersant-solvent system to interact with graphene in many dispersant-solvent systems in order to improve the interfacial properties of graphene in liquid medium and realize the uniform and stable dispersion of graphene. Vallés et al. (2008) experimentally showed that the surface tension of the dispersantsolvent system suitable for the liquid-phase dispersion of graphene should be equal to the surface tension of graphene. Coleman (2009) first proposed the dispersion theory of surface energy in CNT dispersion, extended this theory to the theory of liquid-phase dispersion of graphene, and successfully measured the surface tension of graphene to be 70 mJ/m 2 . Then, Catheline et al. (2011) used KC 8 intercalated graphite to prepare uniform graphene dispersion by liquid-phase exfoliation in NMP. Although the process of KC 8 dissolution could not explained, the Hansen dissolution coefficient theory was used in the NMP dispersion system. The reason for the uniform and stable dispersion of graphene was researched; that is, the free energy of the graphene layer intercalated by KC 8 was much less than the free energy between the original graphene sheets. Coleman (2012) continued to study the surface energy theory. In a series of solvents, graphite powder was ultrasonically exfoliated and dispersed at low power. The surface energy theory was not a strong basis for selecting the best dispersant-solvent system. In considering that not all dispersant-solvent systems have a surface energy close to 70 mJ/m 2 , the graphene in the liquidphase medium was automatically dispersed. Coleman proposed the HSP theory and expressed the content of HSP whose sum of squares is equal to the sum of the squares of the Hildebrand dissolution parameter δ, which satisfies where δ D , δ P , and δ H are the dispersibility, polarity, and hydrogen bond solubility parameters, respectively. Formula (2) shows that the value of HSP must be accurately calculated and a suitable dispersant-solvent system must be selected to obtain stable dispersed graphene dispersions. This process allows graphene to be dispersed stably in liquid medium in the form of single or few layers. Based on the study of Coleman, Yi et al. (2013) extended the study of HSP. In liquid medium, the distance between solvent molecule 1 and solute molecule 2 was defined as Ra = 4 (δ D1 − δ D2 ) 2 + (δ P1 − δ P2 ) 2 + (δ H1 + δ H2 ) 2 [the parameters in this formula were the same as those in Formula (2)]. When Ra was closer, the solubility of HSP was the highest. The theory was used to analyze the acetone/water mixture system. When the mass fraction of acetone was 75%, the best HSP was obtained. Song et al. (2013) found that the non-covalent interaction between graphene and 1-pyridinebutyric acid (PBA) has the lowest degree of damage to the structure of graphene. The results of experiments on various dispersant-solvent systems showed that graphene exhibits the best dispersion in the PBAacetone system. Song et al. (2013) considered that the interaction of δ H and δ H in HSP was the result of the calculation of δ P + δ H = 17.0, which satisfies the predicted value of HSP (generally δ P + δ H = 13-29). At the same time, inspired by the experimental conclusion of Park et al. (2009), graphene bound by non-covalent bonds can also be well-dispersed in the mixture of methanol and water (δ P + δ H > 29). Song et al. (2013) considered that COOH groups on PBA and liquid-phase systems can also be well-dispersed in the mixture of methanol and water (δ P + δ H > 29). Interaction occurs; and graphene, which was non-covalent bonded by PBA, was well-dispersed in liquid medium. At present, the HSP theory is the most breakthrough theory of liquid-phase stable dispersion of graphene. Its advantage is to accurately calculate the dispersibility and solubility of graphene in solventdispersant system. However, the theoretical calculation is complicated and requires high experimental analysis. It has not been widely used in practical experimental research and industrial production.
The concentration of graphene is a key parameter in the liquid-phase dispersion theory of graphene. However, many researchers failed to discuss the initial concentration of graphene, the volume of solvent, and the key parameters (e.g., time, rotating speed, power, and temperature) of liquid-phase exfoliation in the dispersion liquid system and to compare them with the graphene concentrations reported by other groups in the experimental discussion. For example, Yang and Yang (2018) used an electrochemical exfoliation method in the NMP system to strip a graphite rod intercalated with sodium sulfate to prepare a graphene dispersion and only reported the concentration of graphene in the dispersion (11.47 mg/ml). Obviously, it is meaningless to ignore the key parameters of the experiment to discuss graphene concentration. Fortunately, Artur and Paolo (2016) studied the graphene concentration in the graphene liquid-phase dispersion problem and discussed three important ratios: (1) mass ratio of starting graphite to graphene (Y, %); (2) ratio of singlelayer graphene to multilayer graphene (W, %); and (3) ratio of single-layer graphene to the graphene dispersion system (V, %).

Liquid-Phase Stable Dispersion Method π-π Interaction
In a review of the functionalization of non-covalent bonds of graphene, Georgakilas et al. (2016) reported the existence of π-π interactions: (1) π bonds exist between two molecules; (2) the interaction with two the geometry of the molecule is related; and (3) two molecules containing the π structure overlap, and the flatness of the surface of the nano-thin film of the two conjugated structures influences the π-π interaction. Wang et al. (2014) used the density functional theory of verified dispersion correction in combination with the zeroorder symmetric adaptive perturbation theory to study the ππ interaction between graphene and aromatic molecules and to verify the dispersion-modified density functional theory. The zero-order symmetric adaptive perturbation theory of spin component scaling analyses the π-π interaction between graphene and aromatic molecules, which is believed to depend largely on the dispersion force and electrostatic interaction and is found to be more abundant. When a negative electron fluorine atom is substituted to an aromatic group, the π-π interaction becomes stronger and the liquid-phase dispersion property of graphene improves. Unfortunately, Wang et al. (2014) have only made a simple speculation on the π-π interaction and have not been able to further study it. In the field of liquid-phase dispersion of graphene, the nature of the interaction between π-π remains controversial to date (Ding et al., 2018). Wang X. D. et al. (2018) used strong non-covalent bonding (e.g., π-π interaction, hydrogen bonding, and GH-C interaction) between molecules to prevent the agglomeration of rGO lamellae and obtained uniformly dispersed rGO/poly(vinyl alcohol) (PVA) nanocomposites, as shown in Figure 5. The mechanical properties (including Young's modulus and tensile strength) and thermal properties of the composites can be significantly improved by adding a small amount of rGO into PVA.
Only some researchers speculated a reason for the formation of electron-rich exchange between the π-conjugated structure of graphene surface and aromatic molecules. In addition, the π-π interaction between aromatic molecules and graphene not only improves the liquid-phase dispersion but also adjusts the band gap to precisely control its electronic properties (Nduwimana and Wang, 2009;Zhang et al., 2011;Deka and Chowdhury, 2016). The π-π interaction method includes pyrene derivatives, biomolecules, CNT, and GO.

Pyrene Derivatives
Pyrene derivatives were first used in the liquid-phase dispersion of CNTs. Most researchers used NMP as solvent and pyrene derivatives as dispersant to prepare CNTs (Ehli et al., 2006;Meuer et al., 2009)   in a water system. The graphene dispersed in the pyrene derivative system had a dispersion stability of several months. Among different pyrene derivative systems, the effect of the sodium pyrene sulfonate (Py-SASS) dispersion system was the best. As shown in Figure 6A, the concentration of graphene in liquid medium was as high as 0.8-1 mg/ml, as shown in Table 1. Parviz et al. (2012) found that dispersants were automatically adsorbed onto the surface of graphene when the polarity of solvents and dispersants was significantly different. A study on the role of pyrene derivative functional groups in liquid-phase dispersion found that functional groups with high electronegativity effectively promote the adsorption of dispersants onto the surface of graphene ( Figure 6B). Zhang (2016) synthesized a new dispersant, sodium pyrene-butane-1-sulfonate (Py-C 4 -SASS), for the stripping and dispersion of graphene in aqueous media. Parviz et al. (2012) found through a comparative analysis that the new dispersant and Py-SASS, a bismuth-based graphene dispersant, could stably diffuse graphene in an aqueous medium. The exfoliation effect of the two dispersant systems was characterized by Raman and AFM spectroscopic methods (Figure 7A). The ratio of single layer and less layer (≤4 layers) was extremely high, and the radial dimension (0.5 µm) of Py-C 4 -SASS-assisted exfoliated graphene is larger than the radial dimension (0.2 µm) of Py-SASS-assisted exfoliation. The UV-Vis spectroscopy was also performed to test the two dispersion systems. The concentration of graphene in the Py-C 4 -SASS system was twice that of the Py-SASS system ( Figure 7B). Unfortunately, Zhang failed to demonstrate the underlying cause of the difference in the ability to strip and disperse graphene between Py-SASS and Py-C 4 -SASS. Chen C. et al. (2017) used poly(2-butylaniline) (P2BA) as a dispersant to peel off a few layers of graphene in organic solvents. The strong π-π interaction between P2BA with π conjugate structure and graphene solves the problem of graphene agglomeration in organic solvents.

Biomolecules
Among the various reagents that interact with graphene, the advantages of biomolecules are low cost, environmental friendliness, and biocompatibility . Biomolecules that can interact with graphene include proteins and nucleotides (Ehli et al., 2006;Zhang et al., 2011;Deka and Chowdhury, 2016) (DNA/RNA), polysaccharides, and bile salts (Chang et al., 2010;Paredes and Villar-Rodil, 2016). The non-covalent bonding of protein and graphene in an aqueous medium produces a strong steric hindrance effect, which sharply reduces the free energy of graphene in the liquid medium, preventing the mutual graphene layers from each other. The rate of proximity is important to achieve uniform and stable dispersion of the graphene liquid phase (Bourlinos et al., 2009;Laaksonen et al., 2010). Flavin mononucleotides can also be strongly adsorbed onto the surface of rGO (Ayán-Varela et al., 2015;Yoon et al., 2015;Munuera et al., 2016). For example, Vashist and Luong (2015) studied target proteins and found that nucleotides DNA/RNA were highly adsorbed onto the surface of electrochemically and thermochemically reduced rGO. Compared with the surfactant dispersion system, although the concentration of graphene was much lower in the biomolecular system than in the surface active system (Seo et al., 2011;Notley, 2012;Zhu et al., 2014), the advantages of the nucleotide (DNA/RNA) dispersion system,  such as biocompatibility and low cost, are not comparable with those of other dispersion systems (Georgakilas et al., 2012;Ayan-Varela et al., 2017).

Carbon Nanotubes
As a 1D nanomaterial, CNTs form a six-membered ring structure from carbon atoms and contain plenty of conjugated large π bonds (Ehli et al., 2006;. It can interact with graphene through π-π interaction to form a steric hindrance effect (Zhang and Huang, 2000;Zhou and Zhang, 2015;Cao et al., 2017) to achieve the stable liquid-phase dispersion of graphene. Yue et al. (2014) prepared epoxy resin fillers by ultrasonic degradation. A study on the dispersion effects of graphene nanoplates (GNPs), CNTs, and GNP/CNT in liquid medium by the volume sedimentation method found that the GNP/CNT (2:8 mixture) system exhibits the best effect of graphene stable dispersion. As shown in Figure 8a, scanning electron microscopy (SEM) characterization showed that GNP and CNTs can prevent each other from agglomerating. Pure GNP system has higher dispersion than the GNP/CNT system, as shown in Figures 8b-g. Ji et al. (2016) made use of sulfonated CNT (SCNT) and successfully prepared a high concentration of graphene. SCNT played a tremendous role in not only suppressing the agglomeration of graphene but also bridging graphene sheets to improve the interlamellar conductivity. Compared with other dispersant systems in the liquid-phase dispersion of graphene, CNTs and graphene are excellent carbon nanomaterials, and their mixed dispersion system has a wide range of applications in energy storage and conduction (Cui and Li, 2002;Dimitrios et al., 2006;Tang et al., 2014;Zhao et al., 2014;Ji et al., 2016). In the process of preparing PS/PMMA composites, Chen J. W. et al. (2017) developed an effective thermodynamic method to balance the π-π interaction between PS and multi-walled CNTs and the PMMA between multi-walled CNTs. The dipole-dipole interaction with the carboxyl group enables uniform dispersion of the multi-walled CNTs in the composite. Fortunately, the mechanism of the π-π interaction can also be analyzed.

Graphene Oxide
GO is amphiphilic (Geng et al., 2017) and reduces the solidliquid interfacial energy in liquid media (Tian et al., 2011;Shao et al., 2014). Tung et al. (2016) used GO as a dispersant to prepare graphene by stripping graphite in liquid phase. GO as a surfactant not only assists graphite stripping but also renders it uniformly and stably dispersed in liquid phase through the π-π interaction. UV-Vis, SEM, and TEM analyses showed that the stripping yield of single-layer and less-layer graphene in the dispersion is high (40% of the initial graphite mass), and the quality of graphene is excellent (C/O = 21.5, I D I G = 0.12, conductivity 6.2 × 10 4 m s ).

Interaction With Polymers
The non-covalent bond with polymers can achieve both the dispersion of the graphene inorganic solvent and the dispersion of the organic solvent. For example, Choi et al. (2010) first prepared rGO dispersions in pure water by the chemical reduction of GO. Amino groups on polystyrene capped by polymer amine were successfully grafted with carboxylic acid groups remaining on the surface of rGO, forming non-covalent bonds with graphene and transforming graphene from the aqueous phase to the organic phase by ultrasonic degradation. After FTIR and Raman were analyzed, the protonated amino group of polystyrene interacts with the carboxylic acid group on the surface of rGO, realizing non-covalent bonding with graphene and improving the dispersion of graphene in organic media. Tchernook et al. (2014) used polyethylene nanoparticles to add a surfactant to the aqueous system to polymerize, and its dispersion process is shown in Figure 9. The test proves that the graphene is not modified in the aqueous polyethylene/graphene mixed system and the polyethylene nanoparticles disperse the graphene. It has a steric hindrance effect (Li et al., 2012) to achieve the uniform and stable dispersion of graphene.  prepared water-borne graphene dispersions from graphite by electrochemical exfoliation with alkaline solution of hydrolyzed styrene-maleic anhydride copolymer as electrolyte. X-ray powder diffraction (XRD), AFM, and TEM analyses showed that monolayer graphene exists in graphene in the dispersions. FTIR and thermogravimetric analysis demonstrated that the hydrolyzed styrene-maleic anhydride polymer is combined with the graphene surface by a non-covalent bond to obtain a high concentration (>1 mg/ml) and stable aqueous graphene dispersion.

Surfactant Interaction
Surfactants can adjust the surface energy of graphene and different media to promote its stable dispersion. Surfactants have been widely studied as effective agents for the liquidphase stable dispersion of graphene (Coleman, 2012;. Wei et al. (2011) compared the dispersing ability of different surfactants to graphene in the experiment of preparing graphene dispersion solution by exfoliation graphite in liquid phase. Results show that non-ionic surfactant PVP exerts greater dispersing ability, better biocompatibility, environmental friendliness, and lower cost than do other dispersants. Vadukumpully et al. (2011) applied the cationic surfactant hexadecyl trimethyl ammonium bromide (CTAB) as a dispersing agent to ultrasonically exfoliate graphite in a glacial acetic acid system, and they found that a large number of singlelayer and less-layer graphene are stably dispersed in the N, Ndimethylformamide-CTAB system. Tkalya et al. (2012) reviewed the liquid-phase dispersion theory of graphene by surfactants and summarized the surface active agents, including CTAB, PVP, sodium dodecyl sulfate, sodium dodecylbenzenesulfonate (SDBS), lithium dodecyl sulfate, tetradecyl trimethyl ammonium  bromide, sodium cholate, sodium deoxycholate, taurocholic acid, Tween 20, Tween 80, and Triton X-100, among others. Wang et al. (2017) used spray coating to prepare graphene films on glass and n-Si substrates through ultrasonic degradation and utilized anionic surfactant SDBS as a dispersant to prepare graphene dispersion in pure water. The liquid was subjected to UV-Vis analysis to obtain a stable graphene dispersion with a concentration of 15% SDBS. The surface morphology of the graphene has a blade-like edge structure. At the same time, the transmittance of visible light is more than 82% (about seven layers) (Selvakumar et al., 2018). This result proves that SDBS has strong exfoliation and stable dispersing ability to graphene. The surfactant system exerts the best dispersion effect on graphene liquid-phase dispersion, but it cannot avoid introducing other side effects to the dispersion liquid system because of its own defects while solving graphene liquid-phase stable dispersion Bepete et al., 2017;Shabafrooz et al., 2018). For example, Tkalya et al. (2012) found that nonionic surfactants have long hydrophilic end and large volume characteristics. While preventing graphene agglomeration, the large volume effect of non-ionic surfactants can also lead to large distance between adjacent graphene sheets, which affects the conductivity of graphene and the permeability of graphene membranes in the dispersion system. Bepete et al. (2017) applied mechanical exfoliation theory to peel graphite after intercalation of the compound in water to obtain graphenide (negatively charged graphene) solutions. By using the advantages of outstanding volatility of tetrahydrofuran, graphene, and tetrahydrofuran can be mixed in water at room temperature to obtain a high concentration (0.16 g/L) and uniformly dispersed aqueous graphene dispersion. Bellani et al. (2019) prepared graphene dispersion with high concentration and excellent dispersion stability in NMP by wet-jet grinding exfoliation and solvent-exchange process. Compared with traditional exfoliation and dispersion technologies (e.g., ultrasonic exfoliation and mechanical peeling), the above process extremely improved the exfoliated efficiency of graphene and satisfied the industrial requirement of large-scale preparation graphene. Although the dispersion of graphene is not prepared by using any surfactant, the large scale makes use of tetrahydrofuran, and NMP in the industry poses a certain safety hazard to human beings (Rodriguez-Laguna et al., 2018;Shabafrooz et al., 2018). Therefore, the future theory of choosing green dispersants or avoiding using dispersants is the focus of the study on graphene liquid dispersion. For example, Zhang et al. (2018) exfoliated and dispersed graphene with nanocrystalline cellulose (NCC). The addition of NCC considerably reduces the deposition of graphene in the dispersion system. The NCC dispersion method is an efficient, green, and scalable method for preparing graphene fractions. The method of dispersing liquid has a wide potential for the industrial application of graphene in the future. Wang S. F. et al. (2018) used Li + and OH − between graphene sheets. The positive charge of Li + moves to the adjacent C atoms, and the adjacent OH − interacts with it, forming a hydroxyl group added to the surface of the graphene and increasing the hydrophilicity of the graphene (Figure 10A). The graphene dispersion in an aqueous medium was improved without using a dispersant, and a graphene dispersion having a concentration of 0.09 mg/ml was prepared (Figures 10B,C).

Hydrogen Bond Interaction
Hydrogen bonding is a strong intermolecular force (bonding energy is 5-30 kJ/mol). The formation of intermolecular hydrogen bonds is profitable to the dispersion and dissolution between substances . Generally, the certain quantitative oxygen-containing groups will remain on the surface of graphene prepared by oxidation-reduction method (Aradhana et al., 2018) and ultrasonic graphite in water for a long time . The oxygen-containing groups can produce strong hydrogen bond with dispersant-solvent system to improve the interfacial properties of graphene and achieve stable dispersion of graphene in liquid phase Qin et al., 2018). For example,  prepared a PVA/graphene mixture by dispersing rGO uniformly and stably in PVA. FTIR analysis was carried out because the residual oxygen-containing functional groups (such as carboxyl and hydroxyl) on the surface of rGO were combined with hydroxyl groups in PVA through hydrogen bonding, and the liquid-phase dispersion of graphene was well realized. León et al. (2013) used melamine as a dispersant to peel graphite in different liquid media (including organic solvents and aqueous solutions) by ball milling to obtain graphene dispersions with good dispersion stability. In the same solvent system, the exfoliation and dispersion ability of graphite with electron-rich benzene derivatives and triazine derivatives were tested, respectively. The present equivalent benzene derivatives do not exhibit the ability of peeling and dispersing graphite with melamine in the same solvent system. However, the aminotrigine system has a strong peeling and dispersing ability for graphene. The analysis and calculation showed that the aminotrigine derivatives could form a wide range of hydrogen bond 2D structures on the surface of graphene and has a good peeling and dispersing ability for graphene. The formation of hydrogen bonding makes it possible for the multipoint interaction of graphene surface in liquid medium. The stable combination of dispersant and oxygen-containing groups on graphene surface (as indicated by the green dotted line in Figure 5) can improve the dispersion stability of graphene in liquid medium.

SUMMARY AND OUTLOOK
In recent years, the most breakthrough theories in the study of liquid-phase stripping and stable dispersion of graphene are the Hamaker constant and the HSP. The two theories reveal the essence of graphene liquid-phase exfoliation and stable dispersion from the intermolecular force and the interfacial force between solute and solvent, respectively. It is also one of the important approaches to solve the problem of graphene agglomeration and realize the uniform and stable dispersion of graphene liquid phase. Coleman (2012) proposed the Hamaker constant theory to quantify the van der Waals force between graphene lamellae, and the HSP theory was used to analyze the relationship between different dispersant-solvent systems and the interface of graphene. However, the Hamaker constant theory still cannot accurately calculate the van der Waals force between graphene sheets. The theoretical calculation of the Hansen dissolution coefficient is complex, and the level of experimental analysis is high. Therefore, certain challenges remain in the theoretical study of liquid-phase exfoliation and dispersion of graphene. Among the many methods of stable dispersion of graphene, the π-π interaction provides a good protection for the original properties of graphene and an ideal stable dispersion effect. However, the essential mechanism of the π-π interaction is unclear. Only some researchers speculated that the π-conjugated structure on graphene surface is caused by the electronrich exchange phenomenon produced by aromatic molecules. The real reason for the π-π interaction remains unclear. In summary, the current research on the liquid-phase stripping and stable dispersion theory of graphene is not comprehensive enough. Future studies should focus on solving the problem of graphene dispersion in nanomaterial research and practical industrial applications.

AUTHOR CONTRIBUTIONS
LLi, LJ, MZ, and YM jointly determined the research theme. LLi, SY, and LJ organized the database. LLi, MZ, XL, and LLiu conducted statistics and analysis. LLi, LJ, and MZ wrote the first draft. HZ, ZM, and ZL wrote part of the first draft. All authors participated in the revision of the manuscript, read, and approved the submitted version.