**3. Helical insertion of GNRs**

materials may not only be designed to be highly dynamic, displaying adaptive and self‐heal‐ ing properties, but could also help gain an understanding of the rules that govern biomolecule assembly processes [4]. In the self‐assembly process, the hollow interior of SWCNT can serve as nanometre‐sized moulds and templates [5] to control the configuration of other materi‐ als, or as a protective layer [6] to prevent the filler from oxidation and shape fragmentation. The physical and chemistry properties of the heterostructure are expected to be modified due to the interaction between SWCNT and exotic materials [7]. Thus, filling SWCNTs with chosen fillers can produce one‐dimensional nanostructures with exciting new applications [8]. A wide array of fillers, including various metal atoms [9], halides [10], C60 [11, 12] and polyacetylene [13], were found to be filled into the SWCNT, with novel configurations and properties. Strano et al. [14] found that water molecules displayed different phase transition temperatures when confined in the nanotube. C60 [15] or other spherical metallofullerenes [16] were distributed evenly along the axial direction after self‐assembling into the SWCNT just as beans distributed in pod, while those linear fillers such as DNA [17] and polyacetylene [13] were easy to insert into or wrap the SWCNT spontaneously with an interesting helical

Recently, our studies [18, 19] indicated that GNRs, the material with unique optical, mag‐ netic and electrical properties [20–23], have been succeeded in spontaneously and spirally wrapping and inserting into the SWCNT. GNRs possess intriguing electronic structures ranging from semiconducting to half‐metallic, depending on their geometry and dimen‐ sions. Thus, the helical composite structures formed by two novel carbon materials are expected to have excellent properties, such as high carrier mobility and high mechanical strength, to be used in microcircuit interconnectors, nanoelectronic devices and nanosen‐ sors [24–26]. And a convincing model on the interaction between the GNR and SWCNT may inspire great efforts in theoretical studies, synthesis and chemical modifications which focus on their electronic, biological, chemical and even magnetic properties. In this chap‐ ter, the self‐assembly behaviour of GNRs, especially multiple GNRs, induced by SWCNT is studied in detail, while the mechanism and influencing factors of self‐assembly are also described, providing an opportunity for a comprehensive and satisfactory understanding of how to control the composite structure of this SWCNT‐GNR nanomaterial hybrid. This discovery is of great significance in the exploration of the properties of the GNR‐SWCNT system and may expand the applications of GNR and SWCNT in extensive fields involving medicine, chemistry, biology and even fuel cells. The perfect GNR helix in SWCNT may be used to control the band gap of various GNT‐based nanodevices, which will pave the way

Molecular dynamics method is used to study the self‐assembly behaviour of GNRs induced by SWCNTs. The force field of condensed‐phase optimized molecular potentials for atomistic simulation studies (COMPASS) [27] is applied to model the atomic interaction. COMPASS is an ab initio force field which has been parameterized and validated using

configuration.

232 Graphene Materials - Structure, Properties and Modifications

for the progress of nanoelectronics.

**2. Simulation method**

Direct simulations in **Figure 1(a)** show the self‐assembly of six GNRs into the fixed SWCNT. The length of SWCNTs is 98.38Å and six GNRs have the same size of 5.681Å in width and 147.600Å in length. When one end of the GNR is captured by the SWCNT, the GNR adheres onto the tube wall parallel and moves towards the tube. At the beginning, six GNRs fill into the cavity of SWCNT synchronously and spontaneously under the action of van der Waals coupling. When the simulation time is up to 2.3 ns, a clear and symmetrical multiple helical configuration with uniform pitch is formed in SWCNT, like an upgrade of the double helix of DNA [34] or the triple helix of collagen [35] in biology, suggesting nature's preference for helical structure. While the processes of encapsulating and spiralling are terminated, no more axial movement of the helical GNRs has been observed due to the van der Waals potential well, but rotation is continued in their circumferential direction, resulting the GNR tails out‐ side SWCNT keeps a new subtle spiral morphology by movement inertia. In addition, for the helical GNRs in SWCNT, the essential corrugation [36] is suppressed by the enhanced interaction between GNRs and SWCNT, inspiring an efficient way to straighten the twisted and corrugated GNRs. The filling snapshots of six GNRs into the unfixed SWCNT are estab‐ lished in **Figure 1(b)** where a slight deformation of carbon nanotube appears followed by the formation of a perfect single helix of GNRs at 0.95 ns. The unfixed SWCNT‐GNR system takes shorter time to attain equilibrium, demonstrating that the deformation and movement of carbon nanotubes do not hinder the self‐assembly of GNRs, but improve their filling speed.

To further research the helical configuration of the multiple GNRs self‐assembled into SWCNT, different numbers of annular and multilayer GNRs filling into SWCNT are shown in **Figure 2**. For convenience, a new and canonical definition for (*a*<sup>1</sup> , *a*<sup>2</sup> …*a*<sup>i</sup> …*a*N) and [*a*<sup>1</sup> , *a*<sup>2</sup> …*a*<sup>i</sup> …*a*N] is set, where *a*<sup>i</sup> means the number of GNRs in *i*‐layer and parentheses and square brackets refer to the simulated model and result, respectively. **Figure 2** illustrates that different numbers of GNRs can self‐assemble into SWCNTs to form different layers of helical structures. As shown in **Figure 2(a)**, seven GNRs self‐assemble and form a monolayer helical configuration in SWCNT, while for the system of eight GNRs, a GNR deviates from the first helix into the second one in **Figure 2(b)**, sug‐ gesting that there is a quantity threshold for every GNR helix. The helical angle of inner GNRs is easily affected by the outer one to ensure the maximum overlap between GNRs in different heli‐ ces when few GNRs in the inner helix are located, as the royal blue GNR in **Figure 2(b)**. But when several GNRs locate on the inner helix (**Figure 2(c)**), their helical angle does not easily affect each other. This phenomenon might be useful when adjusting or controlling the self‐assembly mor‐ phology of GNRs in SWCNT. Another noteworthy phenomenon is that the outer GNRs can enter the inner layer spontaneously in **Figure 2(d)** and **(e)**, which is closely related to the self‐assembly mechanism of GNRs and is explained in detail in Section 4. When encapsulating the (7, 4, 2) GNRs into SWCNT in **Figure 2(f)**, the GNRs arrange densely to obtain a perfect multiple spiral structure, and no GNR transfers into the inner helix. When setting more GNRs in SWCNT like in **Figure 2(g)** and **(h)**, the self‐assembly process also exists but the helical configuration becomes imperfect due to their undistinguished layer. In addition, the well‐organized spiral configuration of (7, 4, 2) consists of an array of helicoids formed concentrically (shown in **Figure 2(j)**), similar to the stranded wire in cable (shown in **Figure 2(k)**), so the GNR@SWCNT composite structure could be expected to have excellent performance in mechanics and electrics like stranded wire.

**Figure 1.** The insertion and helix‐forming process of six GNRs into SWCNT: (a) snapshots of GNRs encapsulated in the fixed SWCNT; (b) snapshots of GNRs encapsulated in the unfixed SWCNT.

The rolled GNR, which is obtained by connecting the two ends of straight GNR, can also self‐assemble into SWCNT to form a double‐helical configuration. **Figure 3(a)** provides the representative snapshots of a rolled GNR inserted helically into the SWCNT. At the begin‐ ning, a rolled GNR is placed near the SWCNT separated by 5Å. As time continues, the rolled GNR tends to be collapsed. Due to the Van der Waals force [37] between the rolled GNR and SWCNT, the rolled GNR stretches its cross section and is close to the SWCNT at 0.1 ns. Then, the portion of the rolled GNR near the entrance is captured by the SWCNT hollow space and starts to move forward along the inner wall of the SWCNT. After an initial correlation time, the rolled GNR begins to curl in the SWCNT (at *t* = 3.5 ns). When the simulation time reaches 4.8 ns, the rolled GNR shows a clear helical conformation with a large pitch. Then, the helix becomes denser to occupy the entirety of the tubes. Eventually, a perfect DNA‐like double‐helix configuration with remarkably constant pitches is formed in the SWCNT. The result of two rolled GNRs self‐assembling into the SWCNT is shown in **Figure 3(b)**. The two rolled GNRs with the same length are placed at the two ends of the SWCNT. When the simula‐ tion begins, the two rolled GNRs are captured by the inner hollow of the SWCNT simultane‐ ously. After they contact with each other, one of them is pushed down to a two‐layer graphitic nanoribbon, and then spirals like an unrolled GNR. As the simulation time reaches 2.25 ns, helices arise at both of the rolled GNRs. At the simulation time of 5.0 ns, the two rolled GNRs display a clear helical conformation.

self‐assemble into SWCNTs to form different layers of helical structures. As shown in **Figure 2(a)**, seven GNRs self‐assemble and form a monolayer helical configuration in SWCNT, while for the system of eight GNRs, a GNR deviates from the first helix into the second one in **Figure 2(b)**, sug‐ gesting that there is a quantity threshold for every GNR helix. The helical angle of inner GNRs is easily affected by the outer one to ensure the maximum overlap between GNRs in different heli‐ ces when few GNRs in the inner helix are located, as the royal blue GNR in **Figure 2(b)**. But when several GNRs locate on the inner helix (**Figure 2(c)**), their helical angle does not easily affect each other. This phenomenon might be useful when adjusting or controlling the self‐assembly mor‐ phology of GNRs in SWCNT. Another noteworthy phenomenon is that the outer GNRs can enter the inner layer spontaneously in **Figure 2(d)** and **(e)**, which is closely related to the self‐assembly mechanism of GNRs and is explained in detail in Section 4. When encapsulating the (7, 4, 2) GNRs into SWCNT in **Figure 2(f)**, the GNRs arrange densely to obtain a perfect multiple spiral structure, and no GNR transfers into the inner helix. When setting more GNRs in SWCNT like in **Figure 2(g)** and **(h)**, the self‐assembly process also exists but the helical configuration becomes imperfect due to their undistinguished layer. In addition, the well‐organized spiral configuration of (7, 4, 2) consists of an array of helicoids formed concentrically (shown in **Figure 2(j)**), similar to the stranded wire in cable (shown in **Figure 2(k)**), so the GNR@SWCNT composite structure could be expected to have excellent performance in mechanics and electrics like stranded wire.

234 Graphene Materials - Structure, Properties and Modifications

The rolled GNR, which is obtained by connecting the two ends of straight GNR, can also self‐assemble into SWCNT to form a double‐helical configuration. **Figure 3(a)** provides the representative snapshots of a rolled GNR inserted helically into the SWCNT. At the begin‐ ning, a rolled GNR is placed near the SWCNT separated by 5Å. As time continues, the rolled GNR tends to be collapsed. Due to the Van der Waals force [37] between the rolled GNR and SWCNT, the rolled GNR stretches its cross section and is close to the SWCNT at 0.1 ns. Then, the portion of the rolled GNR near the entrance is captured by the SWCNT hollow space and starts to move forward along the inner wall of the SWCNT. After an initial correlation time, the rolled GNR begins to curl in the SWCNT (at *t* = 3.5 ns). When the simulation time reaches 4.8 ns, the rolled GNR shows a clear helical conformation with a large pitch. Then, the helix becomes denser to occupy the entirety of the tubes. Eventually, a perfect DNA‐like

**Figure 1.** The insertion and helix‐forming process of six GNRs into SWCNT: (a) snapshots of GNRs encapsulated in the

fixed SWCNT; (b) snapshots of GNRs encapsulated in the unfixed SWCNT.

**Figure 2.** The initial arrangement and equilibrium configuration of the annular and multilayer GNRs in SWCNT. From (a) to (h), the number of GNRs is 7, 8, 11, 12, 12, 13, 14 and 15, respectively; (j) Corey‐Pauling‐Koltun configuration of the (7, 4, 2) GNRs; and (k) configuration of the stranded conductor in cable.

**Figure 3.** Representative snapshots of (a) a rolled GNR and (b) two rolled GNRs self‐assembling into the SWCNT to form helical configurations. The length of the GNR is 46.40Å.

#### **4. Interaction mechanism**

Why the GNR could self‐assemble into the SWCNT and develop a helical manner? To further study its interaction mechanism, the final geometric configuration and energy evolution of this self‐assembly system are displayed in this section.

**Figure 4** shows the concentration distribution profile of the multiple‐layered GNR@SWCNT system to help characterize the parameters of the multiple helical configurations. Although the peak of concentration distribution profiles is split near the centre of SWCNT due to the rigidity of C‐C bond and the larger curvature in the inner helix, it can be measured that the separation of the adjacent layers is 3.5Å, similar to the basal plane separation in graphite or the wall thickness of the multi‐walled carbon nanotubes (3.4Å). That is to say, GNRs are arranged very closely in the limited space of SWCNTs. In the crowded environment, long molecular chains frequently adopt ordered, helical conformations, which has been confirmed by Snir and Kamien [38]. In their simulation system, a solid, impenetrable but flexible tube was immersed into a solution of hard spheres. It is amazing that the flexible tube will take on a regular helical configuration which takes up the least space. The enlarged environment entropy is considered to be the driving force of helix formation.

**Figure 4.** Concentration distribution profiles of the GNR and SWCNT in the (7, 4, 2) system in the X‐direction.

On the other hand, the spontaneous filling and intertangling of GNRs are also driven by the increasing π‐π stacking area between GNRs and SWCNT, because the system energy is released when carbon six‐membered rings are stacked. In order to confirm this view, a series of simulations about different SWCNTs are studied. **Figure 5(a)** shows that similar multiple helical configurations are taken in cylindrical SWCNT shell to maximize their π‐π stacking area, while the flat GNRs are obtained in the special cubical fixed SWCNT because this config‐ uration can guarantee the largest stacking area with the straight and stationary tube wall. The partial enlarged view of the GNR@SWCNT helical configuration is displayed in **Figure 5(b)**, demonstrating that the GNR ring is offset from SWCNT ring slightly and their centroid‐cen‐ troid distance is 3.858Å, consistent with the typical aromatic‐aromatic offset π‐π stacking [39].

**4. Interaction mechanism**

236 Graphene Materials - Structure, Properties and Modifications

this self‐assembly system are displayed in this section.

entropy is considered to be the driving force of helix formation.

Why the GNR could self‐assemble into the SWCNT and develop a helical manner? To further study its interaction mechanism, the final geometric configuration and energy evolution of

**Figure 4** shows the concentration distribution profile of the multiple‐layered GNR@SWCNT system to help characterize the parameters of the multiple helical configurations. Although the peak of concentration distribution profiles is split near the centre of SWCNT due to the rigidity of C‐C bond and the larger curvature in the inner helix, it can be measured that the separation of the adjacent layers is 3.5Å, similar to the basal plane separation in graphite or the wall thickness of the multi‐walled carbon nanotubes (3.4Å). That is to say, GNRs are arranged very closely in the limited space of SWCNTs. In the crowded environment, long molecular chains frequently adopt ordered, helical conformations, which has been confirmed by Snir and Kamien [38]. In their simulation system, a solid, impenetrable but flexible tube was immersed into a solution of hard spheres. It is amazing that the flexible tube will take on a regular helical configuration which takes up the least space. The enlarged environment

On the other hand, the spontaneous filling and intertangling of GNRs are also driven by the increasing π‐π stacking area between GNRs and SWCNT, because the system energy is released when carbon six‐membered rings are stacked. In order to confirm this view, a series of simulations about different SWCNTs are studied. **Figure 5(a)** shows that similar multiple helical configurations are taken in cylindrical SWCNT shell to maximize their π‐π stacking area, while the flat GNRs are obtained in the special cubical fixed SWCNT because this config‐ uration can guarantee the largest stacking area with the straight and stationary tube wall. The partial enlarged view of the GNR@SWCNT helical configuration is displayed in **Figure 5(b)**, demonstrating that the GNR ring is offset from SWCNT ring slightly and their centroid‐cen‐ troid distance is 3.858Å, consistent with the typical aromatic‐aromatic offset π‐π stacking [39].

**Figure 4.** Concentration distribution profiles of the GNR and SWCNT in the (7, 4, 2) system in the X‐direction.

**Figure 5.** (a) The effect of the SWCNT configurations on the self‐assembly process and (b) partial enlarged detail of the GNR@SWCNT helical configuration where SWCNT is grey and GNR is blue.

**Figure 6** illustrates the potential energy (*E*p), van der Waals interaction energy (*E*vdw) and π‐π stacking area for the GNR@SWCNT system as functions of simulation time in the self‐ assembly process. The evolution of π‐π stacking area and system energy is identical, dem‐ onstrating that the energy decreasing is caused by the π‐π stacking area expanding and the van der Waals interaction between the GNRs and SWCNT plays a dominant role in driving the formation of the helical configuration. The decline of potential energy indicates that the self‐assembly course is spontaneous. Taking the self‐assembly of (7, 4, 1) GNRs as an example, a clear change of energy gradient is shown in **Figure 6(a)**. The initial large energy gradient is contributed by the rapid expansion of π‐π stacking areas as multiple GNRs filling into SWCNT. Then, the energy gradient decreases when the GNRs reach the other end of the car‐ bon nanotube, because the tail of GNRs inserts into the SWCNT slowly with GNRs spiralling and configuration optimizing. As for the (7, 4, 1) GNRs in **Figure 6(b)**, both energy curves are divided into two stages, corresponding to two processes of structural changing. In the first stage from 0 to 1.8 ns, multiple GNRs encapsulate into carbon nanotube along their original directions and remain stable for some time. Then in the second stage (2.3–2.6 ns), two GNRs transfer into inner helical structure individually from the first and the second helix to obtain the largest π‐π stacking area for the system, so the system configuration is changed from the (7, 4, 1) to (6, 4, 2) with a rapid decrease of potential energy, suggesting that this process can improve their stability. In this process, the van der Waals energy drops 0.7 Mcal/mol and the potential energy drops just 0.5 Mcal/mol, indicating that the disappeared energy is stored efficiently by the bending of new inner GNRs.

In principle, the helical self‐assembly is determined by the competition between the van der Waals energy and the bending strain energy of GNRs [40]. The former provides an attrac‐ tive force to spiral GNRs, and the latter usually tends to maintain the GNRs in their original form [41]. During the self‐assembly process, the van der Waals interaction helps the GNRs overcome their energy barrier caused by the bending strain energy, and drives the GNRs to be trapped in the SWCNT and undergo self‐scrolling. Throughout the helix‐forming course, the van der Waals energy has partially converted into the internal energy of GNRs for their mechanical deformation and partially converted into kinetic energy.

**Figure 6.** The potential energy (*E*p) of GNR@SWCNT system and the van der Waals interaction energy (*E*vdw) between the GNR and SWCNT as functions of simulation time in the helical encapsulation of (a)13 GNRs and (b)12 GNRs, and the inserted graph is the corresponding evolution of π‐π stacking area.

#### **5. Dependence of size, function group and initial arrangement of GNRs**

In order to further exploit the effect of system size on self‐assembly, a series of simulations and parameter analyses have been completed. **Figure 7** reveals that all GNRs of different widths in the SWCNT of different diameters are able to form the perfect helical structures, with varying screw pitches. The helical angle increases with GNR width but have no clear monotone change tendency with SWCNT diameter. In the self‐assembly process, more time is needed to form a large helical angle because of the sluggish motion and coordination of multiple GNRs. Besides, the relation between system size and threshold value is also shown in **Figure 7**. It can be seen that the simulated threshold of GNRs in the first perfect helix is inverse with their width due to the limited space inside SWCNT but has a positive correlation with SWCNT diameter.

A simple mathematical model (**Figure 8**) is built to calculate the quantity threshold for the helical GNRs self‐assembled in SWCNT. Three formulas are given as follows.

The threshold of GNRs in *i*‐layer is

$$V\_i = \left\lfloor \frac{(D - 2iA)\pi}{A' + W} \right\rfloor \tag{1}$$

where *D* is the SWCNT diameter, *i* is the sequence of the GNR helix, *A* is the fixed distance between GNRs and SWCNT, *A'* is the fixed distance between adjacent turns of GNRs and *W* is the GNR width. The specific square brackets mean that the calculated value should be rounded down to the nearest whole unit. All distances are in units of Angstroms.

Self-Assembly of Graphene Nanoribbons Induced by the Carbon Nanotube http://dx.doi.org/10.5772/67413 239

overcome their energy barrier caused by the bending strain energy, and drives the GNRs to be trapped in the SWCNT and undergo self‐scrolling. Throughout the helix‐forming course, the van der Waals energy has partially converted into the internal energy of GNRs for their

**5. Dependence of size, function group and initial arrangement of GNRs**

**Figure 6.** The potential energy (*E*p) of GNR@SWCNT system and the van der Waals interaction energy (*E*vdw) between the GNR and SWCNT as functions of simulation time in the helical encapsulation of (a)13 GNRs and (b)12 GNRs, and the

with SWCNT diameter.

The threshold of GNRs in *i*‐layer is

*Vi* = ⌊

inserted graph is the corresponding evolution of π‐π stacking area.

In order to further exploit the effect of system size on self‐assembly, a series of simulations and parameter analyses have been completed. **Figure 7** reveals that all GNRs of different widths in the SWCNT of different diameters are able to form the perfect helical structures, with varying screw pitches. The helical angle increases with GNR width but have no clear monotone change tendency with SWCNT diameter. In the self‐assembly process, more time is needed to form a large helical angle because of the sluggish motion and coordination of multiple GNRs. Besides, the relation between system size and threshold value is also shown in **Figure 7**. It can be seen that the simulated threshold of GNRs in the first perfect helix is inverse with their width due to the limited space inside SWCNT but has a positive correlation

A simple mathematical model (**Figure 8**) is built to calculate the quantity threshold for the

\_\_\_\_\_\_\_\_ (*D* − 2*iA*)*π*

where *D* is the SWCNT diameter, *i* is the sequence of the GNR helix, *A* is the fixed distance between GNRs and SWCNT, *A'* is the fixed distance between adjacent turns of GNRs and *W* is the GNR width. The specific square brackets mean that the calculated value should be

*<sup>A</sup>*′ <sup>+</sup> *<sup>W</sup>* ⌋ (1)

helical GNRs self‐assembled in SWCNT. Three formulas are given as follows.

rounded down to the nearest whole unit. All distances are in units of Angstroms.

mechanical deformation and partially converted into kinetic energy.

238 Graphene Materials - Structure, Properties and Modifications

**Figure 7.** The dependence of the diameter of SWCNT and the width of GNRs in the evolution: (a) helical encapsulation of GNRs with different width into the SWCNT (20, 20); (b) helical encapsulation of GNR into the SWCNT with different diameter.

**Figure 8.** The ideal model of the cross section of GNR@SWCNT helical structure. SWCNT and GNR are referred to as grey circle and blue arc, respectively.

The whole threshold of the self‐assembly system is

$$V\_{\text{total}} = \sum\_{i=1}^{N} V\_i \tag{2}$$

where *N* is the number of GNR helix in SWCNT, which is described by

$$N = \begin{bmatrix} \frac{D}{2A} \end{bmatrix} \tag{3}$$

The above formulas are used to calculate the theoretical threshold for 5.681Å GNR in (20, 20) SWCNT, and the result is consistent with the acquired thresholds of the first, second GNR helix and whole system very well. So the formulas are proved applicable to estimate the quantity threshold for the helical GNR in SWCNT.

Then, the effect of system size on the self‐assembly of rolled GNRs is further studied. **Figure 9(a)** shows the diameter threshold of SWCNT which could be successfully encapsulated by the rolled GNRs with certain widths. It can be observed that the SWCNT diameter should be larger than the threshold to ensure the collapse and insertion of the rolled GNR. When the SWCNT diameter is too small, the van der Waals interaction cannot overcome the large energy barrier in system and the distance between GNR and SWCNT cannot maintain 3.5Å due to the confinement. Consequently, the SWCNT size should be chosen appropriately according to the GNR width to facilitate the self‐assembly. In addition, **Figure 9(b)** demonstrates that a perfect helix inside the nanotube can be obtained when the space inside the nanotube is sufficient, which is independent of the diameter of the SWCNT. But what would happen if the length of the SWCNT is too short? **Figure 9(c)** shows the snapshots of the rolled GNR self‐assembling into the short SWCNT. The rolled GNR inserts into the SWCNT successfully with the trapped portion to form a helix. After the helical configuration in SWCNT is closely packed, the front helical GNR is pushed out and the GNR tail continuously enters the SWCNT in a helical man‐ ner. Finally, the rolled GNR attempts to be symmetric to the tube, with the middle GNR in the short SWCNT forming a perfect helical configuration.

**Figure 9.** (a) The critical diameters of the SWCNTs to guarantee the insertion of the rolled GNR with different widths successfully: GNR(7.38Å)@SWCNT (13,13), GNR(9.84Å)@SWCNT (14,14), GNR(12.30Å)@SWCNT (17,17), GNR(14.76Å)@SWCNT (18,18), GNR(17.22Å)@SWCNT (19,19), respectively. (b) The helical insertion of GNR with the width of 9.84Å into the SWCNT with different diameters. (c) Representative snapshots of a rolled GNR inserting into the SWCNT (25, 25) with the length of 24.60Å to form a helical configuration.

The self‐assembly of GNRs into SWCNT can be utilized to deliver substances into the nanoscale‐ confined space without any other external force. In order to explore its possible applications, dif‐ ferent functional groups are modified on the edge of the rolled GNR. As shown in **Figure 10(a)**, the rolled GNR modified with 133 –OH can self‐assemble into the SWCNT spontaneously to form a helix, and the interval between neighbouring segments is proved to be influenced by the chemical groups –OH. Then, a series of simulations are performed on the rolled GNR that is modified with different number of –OH, as shown in **Figure 10(b)**. The results demonstrate that the rolled GNR can form a helix in the SWCNT when one end of the GNR is modified with 50,100 and 133 –OH. By contrast, the rolled GNR cannot form a helix when fully modified with 200 –OH at one end, and even cannot be trapped by the hollow interior of the SWCNT when modified with 400 –OH. Another functional group, the –H, is also chemically attached to the rolled GNR uniformly. When the dangling σ‐orbitals on the GNR are all saturated by hydrogen atoms, the rolled GNR can still self‐assemble into the SWCNT and form a helix, which is quite different from the result of the GNR modified with –OH, indicating that the functional groups on the GNRs have a significant influence on their self‐assembly. In addition, the spiral of the GNR without any dangling σ‐orbitals also reveals that the van der Waals interaction is the main driving force of the helix self‐assembly, rather than the dangling σ‐orbitals at the edge of GNRs [42]. At the same time, the stability of the final configuration in different systems has been veri‐ fied. **Figure 10(c)** shows that the potential energy of the 200 –OH system decreases to a certain value and then maintains at this level for 7 ns, indicating that the non‐helical configuration is comparably stable. The potential energy in the system of GNR modified with 133 –OH is very different (**Figure 10(d)**) as it is divided into several periods: two 'rapid drop' periods, a 'platform' period and an equilibrium state. So it may be suspected that the modification of the GNR with different functional groups can influence the final configuration of the GNR@SWCNT system.

The whole threshold of the self‐assembly system is

240 Graphene Materials - Structure, Properties and Modifications

*V*total = ∑*<sup>i</sup>*=1

*N* = ⌊

short SWCNT forming a perfect helical configuration.

SWCNT (25, 25) with the length of 24.60Å to form a helical configuration.

threshold for the helical GNR in SWCNT.

where *N* is the number of GNR helix in SWCNT, which is described by

\_\_\_ *D*

The above formulas are used to calculate the theoretical threshold for 5.681Å GNR in (20, 20) SWCNT, and the result is consistent with the acquired thresholds of the first, second GNR helix and whole system very well. So the formulas are proved applicable to estimate the quantity

Then, the effect of system size on the self‐assembly of rolled GNRs is further studied. **Figure 9(a)** shows the diameter threshold of SWCNT which could be successfully encapsulated by the rolled GNRs with certain widths. It can be observed that the SWCNT diameter should be larger than the threshold to ensure the collapse and insertion of the rolled GNR. When the SWCNT diameter is too small, the van der Waals interaction cannot overcome the large energy barrier in system and the distance between GNR and SWCNT cannot maintain 3.5Å due to the confinement. Consequently, the SWCNT size should be chosen appropriately according to the GNR width to facilitate the self‐assembly. In addition, **Figure 9(b)** demonstrates that a perfect helix inside the nanotube can be obtained when the space inside the nanotube is sufficient, which is independent of the diameter of the SWCNT. But what would happen if the length of the SWCNT is too short? **Figure 9(c)** shows the snapshots of the rolled GNR self‐assembling into the short SWCNT. The rolled GNR inserts into the SWCNT successfully with the trapped portion to form a helix. After the helical configuration in SWCNT is closely packed, the front helical GNR is pushed out and the GNR tail continuously enters the SWCNT in a helical man‐ ner. Finally, the rolled GNR attempts to be symmetric to the tube, with the middle GNR in the

**Figure 9.** (a) The critical diameters of the SWCNTs to guarantee the insertion of the rolled GNR with different widths successfully: GNR(7.38Å)@SWCNT (13,13), GNR(9.84Å)@SWCNT (14,14), GNR(12.30Å)@SWCNT (17,17), GNR(14.76Å)@SWCNT (18,18), GNR(17.22Å)@SWCNT (19,19), respectively. (b) The helical insertion of GNR with the width of 9.84Å into the SWCNT with different diameters. (c) Representative snapshots of a rolled GNR inserting into the

*<sup>N</sup> Vi* (2)

<sup>2</sup>*A*⌋ (3)

In addition, we further simulate the influence of the temperature on the final structure of the GNR@SWCNT system. It is surprising that the final composite structure of the GNR@SWCNT system can maintain stability even when the temperature is 1000 K.

**Figure 10.** (a) Self‐assembly snapshots of the GNR modified with –OH entering with the SWCNT (30, 30). The length of the GNR is 135.6Å, the width of the GNR is 14.76Å and the number of the –OH is 133. (b) Final configuration of this GNR modified with functional groups encapsulated into the SWCNT (30, 30). Total potential energy (*E*P) of the GNR‐SWCNT systems as functions of time: (c) the GNR modified with 200 –OH and (d) the GNR modified with 133 –OH.

#### **6. Helical wrapping of GNRs**

The above studies mainly focus on the insertion of GNRs into the hollow SWCNTs, but there is another form of self‐assembly for the GNR/SWCNT system. This section inves‐ tigates the collapse and wrapping characteristic of the rolled GNRs around the SWCNT. **Figure 11(a)** shows the typical snapshots of the rolled GNR spontaneously wrapping on the SWCNT (15, 15). Initially, the circular rolled GNR is positioned in the centre of the SWCNT, with the axis parallel to the SWCNT and the offset distance of about 5Å above the SWCNT. As shown in **Figure 11(a)**, the rolled GNR approaches the SWCNT rapidly due to the strong attractive force between them. During the approaching process, the rolled GNR stretches its cross section from a circle to an oval gradually because the carbon atoms close to the SWCNT endure stronger van der Waals force than the uppers. When the lower atoms in rolled GNR come into contact with the SWCNT, the upper atoms continue approaching the SWCNT, inducing the oval cross‐section collapses to a triangle one [43] that could cover the upper surface of the SWCNT completely. Then the GNR rotates and its cross section changes into a bow shape in 5.5 ns. During the above process, the rolled GNR collapses rather slowly because the van der Waals interaction on the top semi‐circular tube is so weak. But when th'e simulation reaches 9 ns, the collapse of the entire rolled GNR is further accelerated due to the increasing contact area [40] between SWCNT and GNR. The sequen‐ tial collapse of these rings manufactures a linked double‐layered GNR and then results in a double‐walled helix wrapping around the SWCNT with a large helical pitch. Regardless of the angle between the two axes of the rolled GNR and SWCNT, the GNR can wrap around the SWCNT helically. The GNR helix out of the SWCNT wall is similar to the tendril of a morning glory climbing around a tree trunk where growth pattern helps the plant to maxi‐ mize the sunlight absorption in the minimum space and grow strong. **Figure 11(b)** and **(c)** shows the concentration distribution profiles of the final configuration in the X‐direction and the total potential energy of the GNR/SWCNT system in the wrapping process. The former figure reveals that the distance between two layers of the CNR outside the tube is also about 3.5Å, and the latter energy curve has a decreasing tendency with the simulation time and finally reaches the minimum at equilibrium, similar to the insertion of GNRs into SWCNT, meaning that the wrapping self‐assembly of GNR/SWCNT has the same mecha‐ nism as the filling self‐assembly.

The different final configurations of the rolled GNRs adhering on the SWCNTs are summarized in **Figure 12**. The length of all SWCNTs (15, 15) is 73.79Å, and the width of all rolled GNRs is 14.76Å. All rolled GNRs are circular initially. As illustrated in **Figure 12**, the rolled GNRs with a length larger than the threshold of 38.33Å are collapsed. However, the rolled GNRs with diameters between 157.61 and 230.03Å cannot fully collapse, but simply deform from cylindri‐ cal symmetry to bow shapes, which is because the van der Waals interaction in these system cannot overcome the energy required for the fully mechanical deformation of the round GNRs, leading to the mechanically bistable configuration [41]. When the GNR length continues to increase more than 234.30Å, the perfect helix outside of SWCNT can be guaranteed.

**Figure 11.** (a) Representative snapshots of a rolled GNR wrapping on the SWCNT (15,15) to form a spiral configuration, which resembles the tendril climbing around a trunk (hand drawn sketches by Wei Chen). The width of the GNR is 14.76Å and the length of the SWCNT is 73.79Å. (b) Concentration distribution profiles of the GNR and SWCNT in the system in the X‐direction. (c) Total potential energy (*E*P) of the GNR‐SWCNT systems as a function of time.


**6. Helical wrapping of GNRs**

242 Graphene Materials - Structure, Properties and Modifications

nism as the filling self‐assembly.

The above studies mainly focus on the insertion of GNRs into the hollow SWCNTs, but there is another form of self‐assembly for the GNR/SWCNT system. This section inves‐ tigates the collapse and wrapping characteristic of the rolled GNRs around the SWCNT. **Figure 11(a)** shows the typical snapshots of the rolled GNR spontaneously wrapping on the SWCNT (15, 15). Initially, the circular rolled GNR is positioned in the centre of the SWCNT, with the axis parallel to the SWCNT and the offset distance of about 5Å above the SWCNT. As shown in **Figure 11(a)**, the rolled GNR approaches the SWCNT rapidly due to the strong attractive force between them. During the approaching process, the rolled GNR stretches its cross section from a circle to an oval gradually because the carbon atoms close to the SWCNT endure stronger van der Waals force than the uppers. When the lower atoms in rolled GNR come into contact with the SWCNT, the upper atoms continue approaching the SWCNT, inducing the oval cross‐section collapses to a triangle one [43] that could cover the upper surface of the SWCNT completely. Then the GNR rotates and its cross section changes into a bow shape in 5.5 ns. During the above process, the rolled GNR collapses rather slowly because the van der Waals interaction on the top semi‐circular tube is so weak. But when th'e simulation reaches 9 ns, the collapse of the entire rolled GNR is further accelerated due to the increasing contact area [40] between SWCNT and GNR. The sequen‐ tial collapse of these rings manufactures a linked double‐layered GNR and then results in a double‐walled helix wrapping around the SWCNT with a large helical pitch. Regardless of the angle between the two axes of the rolled GNR and SWCNT, the GNR can wrap around the SWCNT helically. The GNR helix out of the SWCNT wall is similar to the tendril of a morning glory climbing around a tree trunk where growth pattern helps the plant to maxi‐ mize the sunlight absorption in the minimum space and grow strong. **Figure 11(b)** and **(c)** shows the concentration distribution profiles of the final configuration in the X‐direction and the total potential energy of the GNR/SWCNT system in the wrapping process. The former figure reveals that the distance between two layers of the CNR outside the tube is also about 3.5Å, and the latter energy curve has a decreasing tendency with the simulation time and finally reaches the minimum at equilibrium, similar to the insertion of GNRs into SWCNT, meaning that the wrapping self‐assembly of GNR/SWCNT has the same mecha‐

The different final configurations of the rolled GNRs adhering on the SWCNTs are summarized in **Figure 12**. The length of all SWCNTs (15, 15) is 73.79Å, and the width of all rolled GNRs is 14.76Å. All rolled GNRs are circular initially. As illustrated in **Figure 12**, the rolled GNRs with a length larger than the threshold of 38.33Å are collapsed. However, the rolled GNRs with diameters between 157.61 and 230.03Å cannot fully collapse, but simply deform from cylindri‐ cal symmetry to bow shapes, which is because the van der Waals interaction in these system cannot overcome the energy required for the fully mechanical deformation of the round GNRs, leading to the mechanically bistable configuration [41]. When the GNR length continues to

increase more than 234.30Å, the perfect helix outside of SWCNT can be guaranteed.

**Figure 12.** The final configuration of rolled GNRs with different lengths adhering on the SWCNT (15, 15). The width of the GNRs is 14.76Å and the length of SWCNTs is 73.79Å.

#### **7. Conclusion**

In summary, this work proves that the multiple and rolled GNRs can self‐assemble into or wrap around the cylindrical SWCNT to form a well‐organized multiple helical configuration which takes the least amount of energy and takes up the least space. The self‐assembly pro‐ cess is driven by the van der Waals interaction between the GNRs and the SWCNT, because this energy is released when carbon six‐membered rings are stacked. The final self‐assembly configuration of GNRs is responsible for the combined effect of system size and functional groups. The sizes of the SWCNTs should exceed a certain value to ensure the insertion of GNRs with certain widths, while they have negligible effects on the formation of the per‐ fect helix after GNR encapsulation. The quantity threshold of GNRs inside SWCNT increases with the SWCNT diameter and decreases with the GNR width, and it can be estimated by the empirical formulas shown in this paper. The modification on GNR edge can influence the final configuration in the self‐assembly process. In addition, the results of this study also reveal that the size of GNRs and SWCNTs should meet some required conditions to guarantee the wrapping self‐assembly outside of the SWCNT in a helical form.

The helical configuration of the GNRs is similar to the ordered, spiral structure of long molec‐ ular chains such as DNA and protein in the confined cell. In fact, the spiral phenomenon is common but interesting in the nature, such as the spiral arrangement of florets in sunflower, the spiral nucleation of SiC crystal in a limited space and the spiral air flow in the centre of a tropical cyclone. What's more, Chinese old Tai Chi totem also emphasizes the philosophy of helical configuration, which implies that the universe derives from a series of helix. Thus, the findings in this paper are of great importance towards the better understanding of the helix, and of significance in predicting and controlling the self‐assembly configuration. Moreover, the unique phenomenon of self‐assembly in the GNRs/SWCNT system has potential advan‐ tages in the application as nanocontainer for drug delivery and molecular transportation. And this composite structure can also be used to fabricate functional nanodevices, such as sensor, nanoelectronic components, integrated circuits and optoelectronic devices.

### **Acknowledgements**

The authors would like to acknowledge the support from the National Natural Science Foundation of China (Grant No.51671114). This work is also supported by the Special Funding in the Project of the Taishan Scholar Construction Engineering and National Key Research Program of China (Grant No. 2016YFB0300501).

### **Author details**

Hui Li\*, Yifan Li and Wei Chen

\*Address all correspondence to: lihuilmy@hotmail.com

Key Laboratory for Liquid‐Solid Structural Evolution and Processing of Materials, Ministry of Education, Shandong University, Jinan, People's Republic of China

## **References**

**7. Conclusion**

244 Graphene Materials - Structure, Properties and Modifications

**Acknowledgements**

**Author details**

of China (Grant No. 2016YFB0300501).

Hui Li\*, Yifan Li and Wei Chen

\*Address all correspondence to: lihuilmy@hotmail.com

Education, Shandong University, Jinan, People's Republic of China

In summary, this work proves that the multiple and rolled GNRs can self‐assemble into or wrap around the cylindrical SWCNT to form a well‐organized multiple helical configuration which takes the least amount of energy and takes up the least space. The self‐assembly pro‐ cess is driven by the van der Waals interaction between the GNRs and the SWCNT, because this energy is released when carbon six‐membered rings are stacked. The final self‐assembly configuration of GNRs is responsible for the combined effect of system size and functional groups. The sizes of the SWCNTs should exceed a certain value to ensure the insertion of GNRs with certain widths, while they have negligible effects on the formation of the per‐ fect helix after GNR encapsulation. The quantity threshold of GNRs inside SWCNT increases with the SWCNT diameter and decreases with the GNR width, and it can be estimated by the empirical formulas shown in this paper. The modification on GNR edge can influence the final configuration in the self‐assembly process. In addition, the results of this study also reveal that the size of GNRs and SWCNTs should meet some required conditions to guarantee

The helical configuration of the GNRs is similar to the ordered, spiral structure of long molec‐ ular chains such as DNA and protein in the confined cell. In fact, the spiral phenomenon is common but interesting in the nature, such as the spiral arrangement of florets in sunflower, the spiral nucleation of SiC crystal in a limited space and the spiral air flow in the centre of a tropical cyclone. What's more, Chinese old Tai Chi totem also emphasizes the philosophy of helical configuration, which implies that the universe derives from a series of helix. Thus, the findings in this paper are of great importance towards the better understanding of the helix, and of significance in predicting and controlling the self‐assembly configuration. Moreover, the unique phenomenon of self‐assembly in the GNRs/SWCNT system has potential advan‐ tages in the application as nanocontainer for drug delivery and molecular transportation. And this composite structure can also be used to fabricate functional nanodevices, such as sensor,

The authors would like to acknowledge the support from the National Natural Science Foundation of China (Grant No.51671114). This work is also supported by the Special Funding in the Project of the Taishan Scholar Construction Engineering and National Key Research Program

Key Laboratory for Liquid‐Solid Structural Evolution and Processing of Materials, Ministry of

the wrapping self‐assembly outside of the SWCNT in a helical form.

nanoelectronic components, integrated circuits and optoelectronic devices.


[29] Wang Q, Duan WH, Liew KM, He XQ. Inelastic buckling of carbon nanotubes. Applied Physics Letters. 2007;**90**(3):033110. doi:10.1063/1.2432235

[15] Smith BW, Luzzi DE. Formation mechanism of fullerene peapods and coaxial tubes: A path to large scale synthesis. Chemical Physics Letters. 2000;**321**(1–2):169–74. doi:10.1016/

[16] Sato Y, Suenaga K, Okubo S, Okazaki T, Iijima S. Structures of D5d‐C80 and Ih‐Er3N@ C80 fullerenes and their rotation inside carbon nanotubes demonstrated by aberration‐ corrected electron microscopy. Nano Letters. 2007;**7**(12):3704–8. doi:10.1021/nl0720152

[17] Caoduro C, Hervouet E, Girard‐Thernier C, Gharbi T, Boulahdour H, Delage‐Mourroux R, et al. Carbon nanotubes as gene carriers: Focus on internalization pathways related to function‐ alization and properties. Acta Biomaterialia. 2017;**49**:36‐44. doi:10.1016/j.actbio.2016.11.013

[18] Li Y, Chen W, Ren H, Zhou X, Li H. Multiple helical configuration and quantity thresh‐ old of graphene nanoribbons inside a single‐walled carbon nanotube. Scientific Reports.

[19] Chen W, Li H. How does carbon nanoring deform to spiral induced by carbon nano‐

[20] Dutta S, Pati SK. Novel properties of graphene nanoribbons: A review. Journal of

[21] Bai J, Cheng R, Xiu F, Liao L, Wang M, Shailos A, et al. Very large magnetoresistance in gra‐ phene nanoribbons. Nature Nanotechnology. 2010;**5**(9):655–9. doi:10.1038/nnano.2010.154

[22] Yang L, Cohen ML, Louie SG. Excitonic effects in the optical spectra of graphene nanorib‐

[23] Ni Y, Yao K, Fu H, Gao G, Zhu S, Wang S. Spin seebeck effect and thermal colossal mag‐ netoresistance in graphene nanoribbon heterojunction. Scientific Reports. 2013;**3**:1380.

[24] Xia D, Xue Q, Xie J, Chen H, Lv C, Besenbacher F, et al. Fabrication of carbon nanoscrolls from monolayer graphene. Small. 2010;**6**(18):2010–9. doi:10.1002/smll.201000646

[25] Wang L, Zhang HW, Zhang ZQ, Zheng YG, Wang JB. Buckling behaviors of single‐walled carbon nanotubes filled with metal atoms. Applied Physics Letters. 2007;**91**(5):051122.

[26] Cicero G, Grossman JC, Schwegler E, Gygi F, Galli G. Water confined in nanotubes and between graphene sheets: A first principle study. Journal of the American Chemical

[27] Sun H. COMPASS: An ab initio force‐field optimized for condensed‐phase applica‐ tions overview with details on alkane and benzene compounds. The Journal of Physical

[28] Bunte SW, Sun H. Molecular modeling of energetic materials: The parameterization and validation of nitrate esters in the COMPASS force field. The Journal of Physical

tube? Scientific Reports. 2014;**4**:3865. doi:10.1038/srep03865

bons. Nano Letters. 2007;**7**(10):3112–5. doi:10.1021/nl0716404

Society. 2008;**130**(6):1871–8. doi:10.1021/ja074418+

Chemistry B. 1998;**102**(38):7338–64. doi:10.1021/jp980939v

Chemistry B. 2000;**104**(11):2477–89. doi:10.1021/jp991786u

Materials Chemistry. 2010;**20**(38):8207–23. doi:10.1039/C0JM00261E

S0009‐2614(00)00307‐9

246 Graphene Materials - Structure, Properties and Modifications

doi:10.1038/srep01380

doi:10.1063/1.2767249

2015;**5**:13741. doi:10.1038/srep13741

