Abstract

Nanocomposites of silicon nanoparticles (Si NPs) dispersed in between graphene layers emerge as potential anode materials of high-charge capacity for lithium-ion batteries. A key design requirement is to keep Si NPs dispersed without aggregation. Experimental design of the Si NP dispersion in graphene layers has remained largely empirical. Through extensive molecular dynamics simulations, we determine a critical NP dispersion distance as the function of NP size, below which Si NPs in between graphene layers evolve to bundle together. These results offer crucial and quantitative guidance for designing NP-graphene nanocomposite anode materials with high charge capacity.


Rechargeable lithium-ion batteries with high energy capacity and long cycle life are in great demand for applications such as portable devices and electric vehicles [13]. Silicon has the highest known theoretical charge capacity, more than ten times higher than the existing graphite anodes [4]. However, the huge volume changes (up to 400%) of silicon upon charging/discharging cause its pulverization and capacity fading, posing a significant challenge to use silicon as anodes [4, 5]. Nanocomposites of silicon nanoparticles (Si NPs) dispersed in between graphene layers are emerging as a novel and promising solution to high-performance Si-based anode materials [6, 7]. Nanoscale silicon materials (e.g., nanoparticles and nanowires) have been shown to be able to sustain much larger volume change without pulverization than their bulk counterpart [6, 8]. The graphene layers of high mechanical flexibility and large surface area serve as a percolating scaffold network to help accommodate the volume change of the Si NPs during charging/discharging and ensure facile electron transport. A key design requirement of such nanocomposite anode materials is to keep the Si NPs evenly dispersed in between the graphene layers [6, 9]. Given their huge surface area over volume ratio, Si NPs that are too close to each other tend to bundle together to reduce surface energy. The aggregations of the Si NPs become more prone to pulverization than dispersed Si NPs. On the other hand, the charging capacity of the nanocomposite anode material scales with the total volume of the Si NPs that can be dispersed in between the graphene layers. Therefore, the optimal design of Si NP-graphene-based anode materials hinges upon determining the maximum dispersion density of the Si NPs in the graphene layers without Si NP aggregation. Such a critical dispersion of the Si NPs in the graphene layers remains elusive so far, as most experimental studies in this regard are largely empirical [6, 7]. To address such a crucial but large unexplored issue, in this paper we conduct extensive molecular dynamics simulations to determine a critical dispersion distance of the Si NPs as the function of NP size, below which Si NPs in between two graphene layers evolve to bundle together.

The dispersion of the Si NPs in between graphene layers is governed by the energetic interplay between the graphene layers and the Si NPs [10]. The typical size of the Si NPs dispersed in between graphene layers is a few nanometers in diameter, much larger than the equilibrium graphene-graphene interlayer distance (~0.34 nm). As a result, the graphene layers are locally separated by the Si NPs and assume corrugated morphology to wrap around the Si NPs [6]. Such corrugated morphology of the graphene layers results from the following energetic interplay. On the one hand, the graphene-graphene interlayer interaction is mainly van der Waals type. Therefore, the graphene layers tend to assume their equilibrium distance to reduce the interaction energy. On the other hand, the corrugation of graphene layers due to Si NP intercalation causes the bending and stretching of the graphene, and thus leads to an increase of the graphene strain energy. The tighter the graphene layers wrap around the Si NPs, the higher the resulting graphene strain energy, and the lower the graphene-graphene interaction energy. The equilibrium morphology of the graphene-NP structure is dictated by minimizing the total free energy of the structure, which includes two dominant contributions, that is, the graphene strain energy and the graphene-graphene interaction energy. Given the large in-plane dimension of the graphene (on the order of microns or higher), the weak interaction between the graphene and Si NPs (also van der Waals type) is of secondary significance in the total free energy. So is the strain energy of the Si NPs, which is negligible due to the weak graphene-Si interaction and out-of-plane flexibility of the graphene.

We first conduct molecular dynamics (MD) simulations to determine the equilibrium morphology of the graphene layers with a single Si NP intercalated in between, which corresponds to the limiting case when the Si NPs are widely dispersed in between graphene layers, thus the interaction among Si NPs is negligible. In consideration of the computational cost, we model a graphene bilayer with a Si NP intercalated in between, corresponding to the two closest neighbor graphene layers that wrap around the Si NP in a real material. Such a simplified model can capture the dominant energetic interplay, given that the graphene-graphene interlayer interaction decays quickly as their distance increases, and thus the contribution from the next neighbor graphene layers to the interaction energy is of secondary significance.

In the MD simulations, each graphene layer has a size of 50 nm by 70 nm. Single-crystal Si NPs of various diameters (i.e., 2~6 nm) are used to study the effect of NP size on the corrugated graphene morphology. The carbon–carbon (C–C) covalent bonds in the graphene are described by the adaptive intermolecular reactive bond order (AIREBO) potential [11]. The nonbonded graphene-graphene interlayer interaction and the graphene-Si interaction are described by two Lennard-Jones pair potentials [12], 𝑉 C C ( 𝑟 ) = 4 𝜀 C C ( 𝜎 1 2 C C / 𝑟 1 2 𝜎 1 2 C C / 𝑟 6 ) and 𝑉 C S i ( 𝑟 ) = 4 𝜀 C S i ( 𝜎 1 2 C S i / 𝑟 1 2 𝜎 1 2 C S i / 𝑟 6 ) , respectively, where 𝜀 C C = 0 . 0 0 2 8 4 e V , 𝜎 C C = 0 . 3 4 n m , 𝜀 C S i = 0 . 0 0 2 1 3 e V , and 𝜎 C S i = 0 . 1 5 n m . To reduce the computation cost, the Si NP is assumed to be rigid (i.e., whose bonding energy remains as a constant). In each MD simulation case, the graphene bilayer, with free boundary conditions for both layers, is prescribed with an initial morphology that, near the Si NP, it bulges out into a conical dome in each layer to house the Si NP inside and far away from the Si NP, it remains flat with an interlayer distance of 0.6 nm. The Si NP is first fixed, and the graphene bilayer is equilibrated for 30 ps to minimize the dependence of the final results on the initial prescribed morphology, then the Si NP is set free to evolve with the graphene bilayer until an equilibrium is reached. The MD simulations are carried out using LAMMPS with Canonical Ensemble at temperature 300 K and with time step 1 fs.

Figure 1 shows the simulated morphology of graphene with a single Si NP intercalated in between of various diameters, 𝑑 N P = 2 n m , 3 n m , 4 n m , 5 n m , and 6 n m , respectively. For visual clarity, the top panel in Figure 1 shows the top graphene layer and the bottom one shows the Si nanoparticle and the bottom graphene layer. If the Si NP is small (e.g., 𝑑 N P = 2 n m ), each graphene layer bulges out into a dome to wrap around the Si NP. As the size of the Si NP increases, the graphene layers start to form ridged morphology to wrap around the Si NP. The locations of such ridges are approximately complementary in the top and bottom graphene layers. The formation of ridged morphology in the graphene bilayer can be understood as follows. As the size of the Si NP increases, forming a smooth conical dome in each graphene layer to wrap around the Si NP requires accommodating the out-of-plane deflection of the graphene bilayer through the in-plane stretching in radial direction and compressing in hoop direction of the graphene bilayer. Such in-plane deformation of the graphene is energetically unfavorable, given graphene’s ultrahigh in-plane stiffness. By contrast, given the much amenable out-of-plane flexibility of graphene, accommodating the out-of-plane deflection of the graphene bilayer by bending each graphene layer to form ridges corresponds to a lower total free energy state, thus is more favorable.

We next consider how the structural morphology of the Si NP-graphene nanocomposites evolves as the dispersion distance between the Si NPs varies. To capture the dominant underlying physics of such a morphologic evolution within reasonable computation cost, we model a graphene bilayer with two neighboring Si NPs intercalated in between. Similar materials parameters and numerical strategy as those described above are used in the MD simulations. For a given Si NP size, MD simulations are conducted over a range of Si NP dispersion distance. Two modes of morphologic evolution emerge from the simulations, as illustrated in Figure 2. If the dispersion distance between Si NPs, denoted by 𝑆 , is sufficiently large (e.g., 𝑆 = 37 nm for 𝑑 N P = 5 n m ), the graphene bilayer in between two neighboring Si NPs is first corrugated to form a long tunnel (Figure 2(b)). The tunnel formation, however, results in increased strain energy in the graphene, therefore is energetically unfavorable. As a result, the long tunnel in the graphene bilayer evolves to disappear, preventing the two Si NPs to migrate toward each other (Figures 2(c) and 2(d)). At equilibrium, the Si NPs remain dispersed and are individually wrapped by the graphene (Figure 2(e)). The corrugated morphology of the graphene near each Si NP is similar to that depicted in Figure 1. If the dispersion distance between the Si NPs is small (e.g., 𝑆 = 25 nm for 𝑑 N P = 5 n m ), the corrugated graphene bilayer between the two neighboring Si NPs can form a short tunnel, which facilitates the migration of the two Si NPs toward each other. As the two Si NPs evolve to come closer, the length of the tunnel becomes shorter (e.g., Figures 2(h) and 2(i)), resulting in a decrease of the strain energy of the graphene, which is thus energetically favorable. Eventually, the two Si NPs evolve to form a bundle, which is wrapped together by the graphene bilayer (Figure 2(j)).

Also emerging from extensive MD simulations is a critical dispersion distance of Si NPs, 𝑆 c r , above which the Si NPs remain dispersed and below which the Si NPs evolve to bundle together. Figure 3 plots 𝑆 c r as a function of the diameter of the Si NP, 𝑑 N P . There is an approximately linear dependence of 𝑆 c r on 𝑑 N P , that is, 𝑆 c r 1 0 𝑑 N P 1 6 . 8 . Results in Figure 3 can serve as guidance for the material and structural design of Si NP-graphene nanocomposites as anode materials for lithium-ion batteries. For example, for a given design criteria of charging capacity, the maximum volume expansion of the Si NPs upon lithiation can be estimated. The corresponding enlarged size of the Si NPs, instead of the size of pristine Si NPs, should be used to determine a critical dispersion distance to prevent the aggregation of the Si NPs.

In summary, we conduct extensive MD simulations to investigate the morphologic evolution of Si NPs intercalated in between graphene layers, from which a critical dispersion distance of the Si NPs can be determined below which neighboring Si NPs evolve to form a bundle. A roughly linear dependence of the critical distance of the Si NPs on its size is revealed. The ridged morphology of graphene has been observed in recent experiments of substrate-supported graphene with NPs intercalated between the graphene and the substrate [13]. It is further shown that tunneling ridges form between NPs below a critical dispersion distance of the NPs, though the rigid substrate to some extent inhibits the neighboring NPs evolving to bundle together. In our MD simulations, we consider the interaction between two closest neighbor NPs of same size intercalated between graphene layers. In real materials, the variation of the NP sizes and the randomness of the NP distribution render more complex morphologic features that pose significant challenge to be characterized solely by modeling. Nonetheless, results from the present study captures the dominant underlying energetics of NP intercalated in between graphene layers. Furthermore, since the graphene-NP interaction is of secondary importance in determining the corrugated graphene morphology, results from the present study are generally applicable to graphene layers with other anode materials in NP form (e.g., SnO2 and TiO2) intercalated in between [1416]. While the mechanistic understanding and quantitative determination of the critical dispersion distance of Si NPs intercalated between graphene layers shed crucial light on optimizing the design of Si NP-graphene nanocomposites as the high-performance anode materials in lithium-ion batteries, systematic experimental explorations are desirable to validate and further leverage the full potential of the present study.

Acknowledgments

This paper is supported by National Science Foundation (Grants no. 1069076 and no. 1129826) and a UMD GRB Summer Research Award. The authors acknowledge the constructive comments from anonymous reviewers.

Supplementary Materials

Two videos showing the morphologic evolution of the Si NP-graphene system corresponding to the dispersed and bundled cases as illustrated in Figure 2(a)–(e) and (f)–(j), respectively.

  1. Supplementary Video 1
  2. Supplementary Video 2