Study on Preparation and Characterization of Graphene Based on Ball Milling Method

Taking advantage of the state-of-the-art graphene preparation technologies in China and abroad, this work studied the preparation processes of graphene for road applications based on the preliminary high-speed vibration ball milling method. This work combined numerical modeling and microscopic experiments. The preparation parameters were optimized through a multifactor and multilevel test scheme. The materials, influencing factors, and parameters in the preparation process were systematically studied. A graphene preparation method was proposed that considered vibration frequency, static filling rate, material-to-ball volume ratio, and the void percentages of media. Flake graphite, aluminum powder, and 304 stainless-steel grit were used as the preparation materials. The preparation parameters and process model were established based on the uniform design method. The preparation parameters were proposed, calculated, and optimized. Microscopic analysis showed that the proposed preparation method can improve the quality of graphene. This study provides a new source of raw materials for the application of graphene in road engineering.

It is particularly difficult to peel off graphene nanosheets from graphite because of the van der Waals force between the graphite layers. In the past decade, graphene preparation technology has developed rapidly in China and abroad with the widespread attention given to carbon nanomaterials. At present, common graphene preparation methods include mechanical exfoliation, oxidation-reduction, liquid-phase exfoliation, electrochemical exfoliation, supercritical fluid exfoliation, chemical vapor deposition (CVD), and epitaxial growth [11][12][13][14][15][16]. Graphenes that are prepared by different methods have different qualities and application prospects. For example, CVD-prepared graphene sheets have a large area, good uniformity, high quality, and controllable number of layers. They are suitable for high-end fields such as optics, electronics, sensors, energy, and biology. However, they are rarely used as road materials because of the strict preparation process and high cost. Currently, graphene oxide and (nonfunctional group) graphene are used in the road field. Relevant research at home and abroad mainly focuses on the performance of graphene composite road materials, For example, Habib et al.'s [17] research suggests that graphene oxide asphalt helps improve the resistance to deformation of the road surface and reduce rutting in the road surface. According to Li's [18] research, the high-temperature performance and plasticity of the graphene-modified asphalt have been significantly improved, He's [19] research shows that the addition of graphene improves the low temperature stability of asphalt to a certain extent. However, there is no research on the preparation process of road graphene. The mechanical exfoliation method is relatively simple. Graphene materials prepared by this method have low size requirements and low cost. High-performance graphene composite road materials can be developed by combining the mechanically exfoliated graphene with road engineering materials to employ the advantages of the graphene. This study was funded by the National Natural Science Foundation of China. According to a patented graphene preparation technology, this study aimed to develop graphene with good quality and low cost based on mechanical exfoliation and ball milling. A numerical method was employed to optimize the graphene preparation parameters. This work provides a new source of raw materials for the application of graphene in road engineering and provides a reference for the development of graphene composite road materials [20][21][22][23][24][25][26][27].

Preparation of Graphene by Ball Milling
2.1. Materials and Equipment. In this study, a QM-3B high-speed vibrating ball mill (specification: 80 ml, vibration frequency: 1200 r/min, motor specification: 220 V, 180 W) produced by Nanjing University Instrument Factory was used. X-ray diffraction (XRD) was performed with a D8-ADVANCE (Bruker, Germany) X-ray diffractometer (Cu-Kα radiation source, voltage 40 kV, current 40 mA, scan speed 0.1 sec/step, wavelength 1.5418 Å). The specific surface area was measured on a Quadrasorb SI analyzer (Quantachrome Instruments). High-resolution transmission electron microscopy (HRTEM) was performed with a Tecnai G2 F20 field emission transmission electron microscope (FEI, USA).

Working
Principle of Ball Milling. In this study, graphene was prepared using a QM-3B high-speed vibrating ball mill (specification: 80 ml; vibration frequency: 1200 r/min; motor specification: 220 V, 180 W) [28]. 304 stainless-steel balls were used as grinding balls. The working principle is as follows: the high-speed rotation of the motor drives the cylinder to rotate at a high speed on the springs. This process is accompanied by simple harmonic vibration. During the rotation process, the frictional resistance between the cylinder and the steel ball causes accelerated motion of the steel ball. The kinetic energy gained by the steel ball is then transferred to the material as the ball collides with the material and other steel balls. As a result, the material is ground by the balls. At the same time, the generated vibration prevents the steel balls from centrifugal movement due to the high-speed rotation and provides more impact energy to ensure the balls have enough kinetic energy. Therefore, the kinetic energy of the high-speed vibrating ball mill is generally greater than that of the planetary ball mill. However, the cylinder volume of the high-speed vibrating ball mill is only 80 ml. Its grinding capacity is small, so this type of ball mill is mostly used for laboratory grinding studies. There are few studies on the working parameters of high-speed vibrating ball mills. The parameters are generally designed based on the theoretical calculation of tumbling ball mills and are optimized based on experience [29,30]. González-Domínguez [31] prepared graphene by planetary ball milling. The raw material ratio is only 25 mg, and the ball mill speed is only 100 r/min. In Buzaglo's [32] study, the content of graphite was only 36 mg. In summary, these two ball milling methods have low production efficiency and high cost, so they are not suitable for road materials.

Influencing Factors of Grinding Efficiency.
During the ball milling process, the rotation and vibration of the cylinder drive the steel balls and materials to rise. Then, the steel balls and materials drop into the cylinder. The steel balls and the materials crush and collide with each other to generate impact energy, shear energy, and friction energy, thereby achieving the purpose of grinding. The grinding process of the ball mill is essentially the movement process of discrete particles. At present, most of the domestic and foreign studies on the grinding process of ball mills determine the relevant parameters by discrete element methods and molecular dynamics [33][34][35]. The ball mill grinding is a complex process with many influencing factors, and these influencing factors restrict and affect each other. The grinding efficiency of the ball mill is determined by the rotation speed of the ball mill, the grinding media (material, size, static filling rate, etc.), the material-to-ball volume ratio (MBVR), the grading of the grinding media, and the grinding time.
(1) The rotation speed is defined as the ratio of the actual rotation speed n to the critical rotation speed n C , which is denoted as Ψ: (2) The static filling rate refers to the ratio of the loose volume of the material to the effective volume of the ball mill in the static state, which is denoted φ, where R is the radius of the cylinder and h is the distance 2 Journal of Nanomaterials from the center of the cylinder bottom to the top surface of the material (3) The MBVR is the ratio of the volume of the material V m to the void volume between the steel balls V t in the static state, which is denoted φ m (4) The void percentage of the media is defined as the ratio of the void volume to the bulk volume of the grinding media, which is denoted μ. V media is the total volume of the media including voids, V single is the volume of a single grinding medium particle, and N is the number of medium particles In studying the grinding process, there are many differences in the boundary conditions between the theoretical calculation and the actual process. For example, the lubricating effect of the graphite will reduce the frictional resistance between the cylinder and the steel balls; relative sliding movement exists between the steel balls and the cylinder, which lowers the proportion of useful work and the mechanical efficiency; the grading of the steel balls will change the way energy is transferred. It is difficult to apply the theoretical calculation results directly to the actual process. Therefore, the process parameters need to be determined by a combination of theoretical calculation and empirical analysis. Peng Huang, Guojun Shi, and others have performed discrete element simulation analysis on ball mills of different specifications and found that the optimal rotation speed was 80-85%. According to production practices, the grinding effect is good when the filling rate is in the range of 20-50%, and the percentage of voids between the steel balls is 0.35 [36][37][38][39][40][41][42][43][44][45][46][47].

Process Parameters of Graphene Preparation by a Ball
Mill. The raw material used to prepare graphene was flake graphite [48]. In the study of Prato and Regev [31,32], melamine and fully conjugated aromatic diluents were used as grinding aids, respectively. However, the surface energy and cohesive energy of molecular crystals are much smaller than that of aluminum; the aluminum powder is used as the grinding aid in this experiment. Graphite is a crystalline mineral of pure carbon. The crystal lattice of graphite has a hexagonal layered structure, and the layered structures can slide against each other in parallel direction. Aluminum has strong plasticity and ductility. The surface of unoxidized aluminum has a strong adsorption capacity. During the ball milling process, the kinetic energy and shear energy generated by the high-speed impact of the grinding balls can break the aluminum to form fresh aluminum surfaces. The surface of the graphite is continuously adsorbed by the freshly generated aluminum surface. The impact from the balls causes the fresh aluminum to peel off the graphite layer by layer, yielding a layered structure of graphite with graphene in the product.
According to the theoretical calculation, the filling rate was kept between 40% and 50%. In Prato's study, ten stainless steel balls (1 cm in diameter) were used for grinding, but one size of steel balls would make the energy transfer efficiency lower, so three grades of steel balls were used (8 mm particle size, 10 g, 5 pieces; 5 mm particle size, 70 g, 137 pieces; 3 mm particle size, 98 g, 889 pieces). As the rotation speed of the high-speed vibrating ball mill was not adjustable, the factory default setting of 1200 r/min was used. It is generally believed that the ball mill has good grinding efficiency when the interstitial filling fraction of the material is 0.6-1.2 [49,50]. The parameter setting of the ball mill is shown in Table 1.

Process Parameters and Selection
Basis. The process parameters were investigated based on patent CN104889411A [29]. The parameter combinations are shown in Table 2.
The densities of aluminum and graphite are similar. To ensure full contact between aluminum and graphite, the content of aluminum should not be less than that of graphite. On this basis, processes G-1 through G-4 with different combinations of aluminum particle size, aluminum content, and grinding time were obtained. It was found that graphene was obtained in G-1 and G-2 but was not obtained in G-3 or G-4. This may be because the energy of the steel balls could not be effectively transferred to the material when the particle size of aluminum powder was too small. 304 stainless-steel grit with a particle size of 0.4 mm was added to effectively disperse energy and increase the number of contact points for load transfer. Thus, process combinations G-5 through G-8 were obtained. The volume of steel grinding grit was determined based on the volumes of graphite and Journal of Nanomaterials aluminum powder. The volume of the steel grit was chosen to be equal to that of the material to ensure effective contact and was then converted to the corresponding mass.
2.6. Microscopic Analysis of Graphene. The ground products showed two types of macroscopic appearance. As shown in Figure 1, one was silver-gray powder with metallic luster, and the other was pure black powder. According to the optical properties of graphene, the black powder can be used as an indicator to preliminary determine whether graphene is successfully prepared.
When light is shone on graphene, the electrons in the valence band (the energy band occupied by valence electrons) are excited to the conduction band (the energy space formed by free electrons) by absorbing the energy of the photons. The photoconductivity of single-layer graphene depends on the fine structure constant α.
The linear transmittance of suspended single-layer graphene is [51]: where α is the fine structure constant, which is about 1/137. The light absorbance of single-layer graphene is 2.3%. In the visible light region, the reflectance of single-layer graphene is less than 0.1%, and the reflectance of 10-layer graphene is only 2%, which suggests that graphene has extremely strong light absorption. The transmittance of single-layer graphene is as high as 97.7%, and only 0.1% of the visible light will be reflected. Therefore, single-layer graphene is transparent under light. Within several layers of graphene, each additional layer increases the absorbance by 2.3%. The few-and multiple-layer graphenes appear black due to their absorbance and reflectance [52]. Graphites with different structures have different physical properties, and the 50-mesh flake graphite itself has a silver-gray metallic luster.
The products of processes G-1 and G-2 were black powders (Figure 1(b)), which preliminarily indicated that the products contained a large amount of graphene. The grinding time of G-1 was too long, and the efficiency of G-1 was low compared to that of G-2. The G-2 product was selected for the microscopic characterization and was compared with DY-1 (graphene produced by Deyang Carbonene Technology Co., Ltd.). X-ray diffraction (XRD), specific surface area measurement, and HRTEM were performed for the microscopic characterization.
2.7. XRD Experiment. The principle of X-ray crystal structure analysis is to measure the angle θ using the diffraction of Xrays with a known wavelength and then calculate the crystal plane spacing d. X-rays are electromagnetic waves with short wavelengths (approximately 20 to 0.06 Å). They can generate diffraction peaks when passing through the gaps between the graphite layers. The interlayer spacing of graphite is 3.35 Å, and the corresponding peak of the 3.35 Å crystal plane spacing is approximately 26°. The lower the peak intensity, the lower the crystal content and the higher the graphene content. XRD was performed with a D8 Advance diffractometer (Bruker, Germany). The result of the flake graphite analysis is shown in Figure 2. Figure 2 shows that the peak of the graphite crystal appeared at approximately 26°. The peak value of the flake graphite crystal reached 374,826, and the relative strength of the crystal was high, which indicated a good crystallinity of the sample. XRD was also performed on DY-1 and G-2 and is shown in Figure 3.   Journal of Nanomaterials Figures 2 and 3 show that the peak value of DY-1 was 4090, which was higher than that of G-2 (1767). This result indicated that G-2 had lower graphite crystal content, higher graphene content, and therefore higher product quality than DY-1.
2.8. Brunauer-Emmett-Teller (BET) Experiment. The BET method for specific surface area measurement is widely used in the study of particle surface adsorption. It has become the most widely used and most accurate method in contemporary research [53]. The specific surface area measurement in this study was performed with a Quadrasorb SI instrument (Quantachrome Instruments). The result is shown in Table 3. Table 3 shows that the specific surface area of G-2 was approximately 20 times that of DY-1. The adsorption capability of G-2 was much higher than that of DY-1.

HRTEM Characterization of Graphene.
A HRTEM is a high-resolution transmission electron microscope that consists of an electro-optical part, a vacuum part, and an electronic part. An accelerated and focused electron beam is projected onto a target sample, and the electrons collide with the atoms in the sample and scatter. The scattering angle of the electron beam is related to information such as the density, thickness, and spacing of the sample. Thus, the electrons after the collision carry the relevant information of the sample, and such information can be converted back to images [54]. In this study, HRTEM was performed on a Tecnai G2 F20 field-emission TEM (FEI, USA). The characterization results of the two samples are shown in Figures 4 and 5. Figures 4 and 5 show obvious layered structures in the two samples. G-2 and DY-1 show shadows in the 50 nm scanning range. Scanning the shadow in the 5 nm range, a large number of layers are observed. This is because the gra-phene obtained is not completely flat, and some of it has wrinkles due to curling. It can be seen that G-2 and DY-1 are not all single-layer graphene. It contains a small amount of single-layer graphene, but multilayer graphene [55] is the majority, so the sample of G-2 is not uniform. G-2 contained more dark spots in the 50 nm scanning range, which indicated that the impurity content in G-2 was higher than that in DY-1. Further cleaning and purification should be performed on G-2. In the application of engineering materials, graphene-modified asphalt has a very small amount of graphene. The effect of trace impurities in G-2 on engineering materials has yet to be verified.
Microscopic analysis and comparison of the two samples showed that the folds of sample DY-1 are more dense, which indicates that the G-2 prepared by the vibrating ball mill was significantly better than the DY-1 purchased on the market. However, the quality of the G-2 graphene was still not high. Therefore, further optimization of the process was needed.

Process Optimization Design for Vibration
Ball Milling 3.1. Factor Selection. According to Table 2 and the preliminary test results, it was considered that the G-1, G-2, G-7, and G-8 processes were best out of the eight preparation processes. Therefore, XRD was performed on samples prepared by these four processes, and the results are shown in Figure 6. It can be seen from Figure 6 that the peak intensity of G-1 was lower than that of G-2, but the efficiency of G-1 was low because the grinding time was too long. Therefore, the grinding time was fixed at 30 min to ensure grinding efficiency. G-8 graphite had the lowest peak value, so the optimization was performed based on the G-8 process.
Through previous experience and exploration, we found that many factors affected the quality of the graphene, such as the particle size of the material, MBVR, grinding time, and the grading and content of the steel grit. Since there are many influencing factors and these factors interact with each other, some variables were determined first to reduce the number of experiments. The particle size of the aluminum powder was selected first by the control variates method based on the G-8 process. However, a large amount of chromium impurity was found in the G-8 sample in the subsequent product purity test. It was determined that the chromium originated from the 40-mesh aluminum. Therefore, aluminum powders with four particle sizes (40 mesh, 100 mesh, 200 mesh, and 300 mesh) were reselected for testing. The process parameters are shown in Table 4. Table 4 shows that the XRD peak intensity of G-9 multilayer graphene was higher than that of G-8 multilayer graphene, which indicates that chromium was helpful for the grinding of graphene. Although the G-8 multilayer graphene had higher quality, the chromium in the product could not be removed in the posttreatment, so the produced graphene could not be purified. In addition, the chromium content could not be quantitatively analyzed, and the calibration was difficult. Therefore, process G-8 was discarded. The experimental results of processes G-9 to G-12 were  5 Journal of Nanomaterials compared. The relationship between the peak value and the aluminum particle size is shown in Figure 7.
It can be seen from Figure 7 that under the same process conditions, the XRD peak intensity of the produced multilayer graphene was the lowest when the aluminum particle size was 200 mesh. Therefore, the aluminum particle size was first set at 200 mesh, and the design of process parameters was carried out based on this condition.

Process Parameter Design.
Uniform experimental design was adopted for the process parameter design. The uniform experimental design was proposed by Chinese mathematicians Kaitai Fang and Yuan Wang in 1978. It is one of the main methods of fractional factorial design. Compared with the common orthogonal experimental design, the uniform experimental design gives the researchers more choices, and the desired results can be obtained with fewer experiments [56,57]. The percentage of voids μ between the grinding balls was measured to be 0.35. Ball mills have a good grinding efficiency when the interstitial filling fraction of the material is 0.6-1.2 [51][52][53]. This study selected 1.2 as the limit value, and the amount of the material was controlled in this range. The limit amount of flake graphite (50 mesh) and aluminum (200 mesh) were both set to 5 g based on the size of the sample cylinder. Nine levels were fractioned in the range of 1 g to 5 g. Three particle sizes (0.2 mm, 0.4 mm, and 0.8 mm) of steel grit were selected. To construct a uniform design table with equal-level factors, the steel grit was set to have a single particle size or a mixture of multiple particle sizes. The three particle sizes were configured in accordance with the configuration principle of grinding balls. To facilitate the energy dispersion and increase the contact points, the following principles were employed: equal proportions, fewer 0.4 mm balls but more 0.2 mm and 0.8 mm balls, and fewer balls with larger sizes. The total mass of steel grit was kept unchanged at 6 g to ensure that the bulk volumes of steel grits were basically the same. In this paper, three factors were selected: X 1 -graphite content, X 2 -aluminum content, and X 3 -steel grit grading, each having 9 levels. The factor-level table is shown in Table 5.
The design was assisted by the Data Processing System (DPS) analysis software. DPS is a data processing system that integrates functions such as numerical calculation, statistical analysis, model simulation, and chart and table creation. It     Journal of Nanomaterials has many experimental design functions, including uniform design and mixture uniform design [58,59]. The obtained uniform design scheme and the corresponding XRD results are shown in Table 6. A mathematical model was established based on Table 6 with X 1 , X 2 , and X 3 as independent variables and Y 1 as a response variable. The independent variables X 1 and X 2 and the dependent variable Y 1 could be directly input into the model. However, the independent variable X 3 is a single-particle-size steel grit or a triple-particle-size steel grit (both types of grit had the same mass) and could not be directly inputted into the model. Therefore, the steel grit grading was first converted into a specific surface area and then inputted into the modeling calculation. The specific surface area of steel grit was calculated as follows: ρ-steel grit density-7.93 g/cm 3 D-particle size The specific surface areas of the steel grits with particle sizes of 0.8 mm, 0.4 mm, and 0.2 mm calculated by equation (6) were 3/3172 m 2 /g, 3/1586 m 2 /g, and 3/793 m 2 /g, respectively. Then, the specific surface area corresponding to different X 3 can be obtained. At the same time, the G-13, G-14, and G-21 test groups with large Y 1 and obviously poor process capabilities were excluded from the test scheme. The results after data processing are shown in Table 7.

Model Establishment and Solution
A mathematical model was established using graphite concentration (X 1 ), aluminum concentration (X 2 ), and the specific surface area of steel grit (X 3 ) as independent variables and the XRD peak intensity (Y 1 ) as the response value. The model was optimized to obtain the preparation parameters. Graphene was prepared based on the simulated preparation parameters and was analyzed to verify the reliability and significance of the model.

Model Establishment.
Partial least squares regression analysis is a method to study the quantitative relationship between multiple explanatory variables and multiple response variables. Therefore, the partial least squares quadratic polynomial regression, partial least squares (considering interaction terms) regression, and partial least squares (considering quadratic terms) regression were first selected in the DPS software to establish the model. The model calculation parameters are shown in Table 8. Table 8 shows that when the model was established using the partial least squares (considering quadratic terms) regression, the sum of the squared errors was the smallest, the coefficient of determination R 2 was the highest (approaching 1), and the PRESS statistic was the smallest. Therefore, we chose to use the partial least squares (considering quadratic terms) regression for modeling.
First optimization (1stOpt) is a set of integrated software tools for mathematical analysis and optimization. It was independently developed by 7D-Soft High Technology Inc. (China) [60] and has strong computing ability in the fields of nonlinear regression, curve fitting, and parameter estimation for nonlinear complex engineering models. Therefore, the 1stOpt software was used for modeling, and the model was compared with that established by the DPS software. The partial least squares (considering quadratic terms) model was fitted by 1stOpt [61,62]. The model parameters are shown in Table 9.
The closer the correlation coefficient is to 1, the better the model is; the smaller the sum of squared residuals is, the better the model is. Therefore, the 1stOpt model was better than the DPS model.  The 1stOpt model is shown in equation (7): ð7Þ X 1 -graphite content (g) X 2 -aluminum content (g) X 3 -specific surface area of steel grit (m 2 /g) Y-XRD peak value 4.2. Optimal Solution of the Model. The approach to obtain the optimal solution of the model is to find the values of X 1 , X 2 , and X 3 so that Y is equal to 0. Y = 0 indicates that the crystal structure of graphite has completely disappeared, the graphene content in the product is high, the sample quality is high, and the purpose of optimization is achieved. Equa-tion (7) was solved by 1stOpt taking into consideration the XRD results of the nine tests in the uniform design. When X 1 was greater than X 2 , the XRD peak intensities of the obtained sample were all very high. To achieve a better exfoliation effect, the contents of graphite and aluminum should satisfy X 2 ≥ X 1 , and the total mass of X 1 and X 2 must be kept within a reasonable range. If the volume is too large and the optimal range of the interstitial filling fraction (0.6-1.2) is exceeded, the grinding efficiency will decrease. Therefore, the constraint conditions were set as follows: X 1 ≤ X 2 and X 1 + X 2 ≤ 7 g (interstitial filling fraction = 0:6 was set as the boundary condition). Since there are multiple solutions to the equation, the optimized solution summary is shown in Table 10. The X 3 obtained in Table 10 was the specific surface area, which had to be converted into a preparation parameter. The values of X 3 in Table 10 were compared with the values of X 3 in Table 7, but no corresponding values were found. Therefore, the particle size combination closest to each solution was first determined based on Table 7. Two specific surface areas closest to each X 3 solution were chosen from Table 7 and were then converted to a combination of the steel grit particle sizes. The optimized experimental scheme is shown in Table 11.

Microscopic Analysis.
Samples were prepared according to the six schemes in Table 10 and were characterized by XRD. The test results of samples prepared by the above six processes were compared with that of the best sample in Table 7 (G-19). The results are in Figure 8.
G-19 was the best process in the uniform design. G-22 to G-27 were processes obtained by model optimization. The XRD peak intensities of the G-22 to G-27 samples were lower than that of the G-19 sample, which indicated that the crystal contents of graphite in G-22 to G-27 were further reduced, and the qualities of multi-layer graphene in G-22 to G-27 products were improved. The effectiveness of process parameter optimization was confirmed. The peak value of the G-22 product was 134, so G-22 was the best process in the current optimization scheme.
Representative processes (G-2, G-19, and G-22) were selected for specific surface area measurement. The BET results are shown in Table 12. Table 12 shows that the specific surface area of the G-22 product was the largest. At the same time, the peak intensity of the G-22 product was the lowest by XRD. Therefore, the sample prepared by the G-22 process had the highest graphene content, and G-22 was the best process.

Conclusion
(1) Multilayer graphene with good quality can be effectively prepared by mechanical exfoliation of graphite by ball milling (2) The process parameters of multilayer graphene production by a vibrating ball mill were designed and optimized using the uniform design method and mathematical modeling. The preparation materials  Table 6: Uniform design scheme (note: the total mass of X 3 is kept at 6 g; Y 1 is the XRD result).
Test no  (3) Many factors affect the multilayer graphene production by ball mills. In this study, the vibration frequency was 1200 r/min, the static filling rate was 40%, the MBVR was 51%, and the void percentage of the media was 0.38. More influencing factors can be selected in the future for further optimization of the preparation process (4) Our findings show that for the preparation of multilayer graphene by a vibrating ball mill, the process parameters can be designed and optimized with the assistance of mathematic modeling. However, the yield of multilayer graphene produced in the laboratory was low (5) This study provides a parameter design reference for the preparation of multilayer graphene by planetary ball mills

Data Availability
The data used to support the findings of this study are included within the article.

Conflicts of Interest
The authors declare no conflict of interest.