Yufeng
Guo
* and
Wanlin
Guo
State Key Laboratory of Mechanics and Control of Mechanical Structures and MOE Key Laboratory for Intelligent Nano Materials and Devices, Institute of Nanoscience, Nanjing University of Aeronautics and Astronautics, Nanjing, 210016, China. E-mail: yfguo@nuaa.edu.cn; wlguo@nuaa.edu.cn; Fax: +86 25 84895827; Tel: +86 25 84890513
First published on 2nd November 2012
We studied the thermophoretic motion of wrinkles formed in substrate-supported graphene sheets by nonequilibrium molecular dynamics simulations. We found that a single wrinkle moves along applied temperature gradient with a constant acceleration that is linearly proportional to temperature deviation between the heating and cooling sides of the graphene sheet. Like a solitary wave, the atoms of the single wrinkle drift upwards and downwards, which prompts the wrinkle to move forwards. The driving force for such thermophoretic movement can be mainly attributed to a lower free energy of the wrinkle back root when it is transformed from the front root. We establish a motion equation to describe the soliton-like thermophoresis of a single graphene wrinkle based on the Korteweg–de Vries equation. Similar motions are also observed for wrinkles formed in a Cu-supported graphene sheet. These findings provide an energy conversion mechanism by using graphene wrinkle thermophoresis.
How to manipulate and control wrinkle morphology is a key issue for utilizing wrinkled graphene characteristics to design functional devices. Previous experiments have shown that the geometry of pre-existing wrinkles or ridges in graphene sheets can be modified using atomic force microscopy (AFM) tips.12,15 Moreover, thermally induced compressive strain creates periodic graphene wrinkles with controllable wavelength and amplitude.19,20 During the working process of graphene-based devices, temperature might vary in different parts of the graphene. Wrinkles are actually hard to avoid in the graphene plane while the physical properties of graphene are closely related to its structure. However, the effects of temperature change on wrinkle structure and dynamic behaviors, to our best knowledge, are seldom considered in previous studies. On the other hand, thermal load, such as thermophoresis caused by temperature gradient, is a commonly used way to explore thermally actuated phenomena at the nanoscale, and also useful for converting thermal energy into desired mechanical energy.21–25 Further investigations on the response of wrinkles to thermal influence are therefore necessary for actual applications of wrinkled graphene.
In this study, we show by nonequilibrium molecular dynamics (MD) simulations that a single wrinkle formed in a substrate-supported graphene sheet will move from hot to cold regions when subjected to a temperature gradient. Thermophoresis of the wrinkle is a constant acceleration motion, and its acceleration increases linearly with increasing the temperature deviation between the heating and cooling sides of the sheet. The atoms of the single wrinkle drift upwards and downwards, moving the wrinkle forwards. This kind of motion is similar to the dynamic behaviors of a solitary wave or single-soliton described by the Korteweg–de Vries (KdV) equation, where the solitary wave is localized and travels with unchanged shape.26,27 The driving force for such thermophoretic movement can be attributed to the fact that the front root of the wrinkle transforming into the back root will have a lower free energy. Based on the KdV equation, we establish a motion equation to describe the soliton-like thermophoretic movement of a single graphene wrinkle. Furthermore, similar accelerated motion has been observed for two wrinkles formed in a Cu-supported graphene sheet, which move together from the hot to cold regions, but the wrinkle in the hotter region has a higher acceleration than that of the wrinkle in the cooler region. Our results unveil a novel energy conversion through thermophoresis of a graphene wrinkle.
![]() | ||
Fig. 1 (a) A wrinkle is created by applying a −6.1% uniaxial compression in the x direction of the graphene sheet supported by the graphite substrate (top), and a snapshot of the thermophoretic movement of the wrinkle under a temperature deviation of 100 K between the heating and cooling regions (bottom). (b) The temperature distribution along the graphene sheet averaged from 20 to 100 ps before the wrinkle starts to move. Here l is the length of the graphene sheet between the heating and cooling regions. The averaged temperature along the graphene sheet is about 300 K. |
![]() | ||
Fig. 2 (a) Translational displacement of the wrinkle of ε = −6.1% as a function of time under different temperature deviations ΔT. Here we set the time to zero when the wrinkle begins to move. The solid lines are the quadratic polynomial fits to the MD data. The inset shows the structure between the back root rb and front root rf of the wrinkle. (b)The variation of the translational acceleration aw with temperature gradient ΔT/l. |
Temperature also varies along the wrinkle, and the atoms at higher temperature possess higher potential energies. The translational movement of the wrinkle is actually activated by its carbon atoms detaching from substrate (ascending in the z direction) and attaching (descending) with substrate under the temperature gradient. Fig. 3(a) and (b) show some snapshots of the motion trajectory of the chosen atoms (cyan color) at the wrinkle front root under ΔT = 100 K (ΔT/l = 3.84 K nm−1). The atoms at the front root firstly rise up and transform into the wrinkle top, then they fall down to the substrate and become the back root. This kind of motion is similar to the behavior of a solitary wave. The displacements of the chosen atoms in the z and x directions with time are shown in Fig. 3(c). Besides the vibration in the z direction, the atoms are also moving forwards along the x direction. The translational displacement in the x direction for the atoms completing the transition from the front root to the back root is just the value of the wrinkle length subtracting its width, while the wrinkle is moving forward with a distance equal to the wrinkle length. The corresponding free energies that are the sum of kinetic and potential energies for the chosen atoms are shown in Fig. 3(d). The energy increases as the atoms at the front root rise up, as denoted by the red arrow in Fig. 3(d), which means that the ascending process is endothermic. On the contrary, the energy decreases when the wrinkle top becomes the back root, and the descending process is exothermic. It is clearly shown in Fig. 3(d) that the energy change ΔEb between the wrinkle top and back root is larger than that of ΔEf between the top and front root, here ΔEb = Etop − Erb and ΔEf = Etop − Erf. The atoms at the wrinkle front root move forward to a colder region when they transform into the back root, as shown by the x displacements in Fig. 3(c), so that Erb is lower than the initial Erf. As a result, the thermophoresis of the wrinkle can occur only if ΔEb > ΔEf. The energy difference between Erb and Erf depends on the magnitude of the temperature gradient, which is the main driving force for the thermophoretic motion of the graphene wrinkle.
![]() | ||
Fig. 3 The chosen carbon atoms highlighted by a cyan circle (a) ascend and (b) descend when the wrinkle is moving forward under ΔT/l = 3.84 K nm−1. The black arrows denote the locations of these atoms, and the horizontal arrows show the translational direction of the wrinkle. (c) The displacements in the z and x directions, and (d) the free energy with time for the chosen atoms shown in (a) and (b). The shaded region represents the descending of the chosen atoms. Erf, Erb and Etop are the free energies for the chosen atoms located at the front root, back root and top of the wrinkle, respectively. Here we set the time to zero when the atom begins to move. |
With increasing compression on the graphene sheet, the wrinkle height and length increase but the width decreases, and the shape of the wrinkle is changed as well, as shown by the insets in Fig. 4. Applying the thermal difference of ΔT = 200 K (ΔT/l = 7.69 K nm−1) can still motivate the wrinkle to move from the hot site to the cold site. Similarly, all the MD motion data can be well fitted by a quadratic polynomial (Fig. 4), but the translation acceleration decreases with the strain increase. This is because the change of wrinkle shape increases its effective mass and extends the time for the front root transforming into the back root. However, the thermally driven motion will be halted when the highly deformed wrinkle is formed. Under ε = −21.2%, the front and back bottoms of the wrinkle bind together by non-bonding van der Waals interactions. This wrinkle cannot move in the presence of ΔT = 200 K or even higher, as the interlayer relative friction between the wrinkle bottoms is large enough to lock its movement.
![]() | ||
Fig. 4 Translational displacements of the wrinkles of ε = −10.6% and −19.7% as a function of time under ΔT/l = 7.69 K nm−1. The solid lines are quadratic polynomial fits of ![]() |
The thermophoretic movement and structural shape of the graphene wrinkles are very similar to a solitary wave and its motion characteristics. The motion equation of a single-soliton is usually described by the KdV equation:26,39
One solution to this equation is: , here c is considered to be the velocity of the soliton. According to our MD simulation results and fitted quadratic polynomial, we assume that the wrinkle at different positions has a structural form: uw = 3cwsec h2ξ, here
. cw, a and b are constants, and cw can be adjusted to make uw consistent with the shape of the wrinkle. Then we can demonstrate that the uw satisfies the motion equation below:
As the shape of the wrinkle is unchanged during traveling, the wrinkle motion can be considered as a kind of single-soliton. Compared with the KdV equation, there is an additional item . This is the driving force from thermal gradient for the accelerated motion of the wrinkles. The Hamiltonian of this soliton derived from the motion equation is:
To understand the effect of substrate, we have studied thermophoretic motion of a wrinkle formed in a Cu-supported graphene sheet. We chose a graphene sheet with a length of 58.8 nm placed on a (100) Cu substrate. The C–C and Cu–Cu interactions are described by the AIREBO potential and embedded atom method (EAM) potential,40 respectively. Here we use the 6–12 L-J potential to describe the C–Cu interaction with parameters σc–cu = 3.225 Å and εc–cu = 0.019996 eV, which have been proved to be valid for non-bonding C–Cu interaction.41,42 The interlayer interaction between the graphene and the Cu substrate is stronger than that between two graphene sheets. A wrinkle is created by applying a uniaxial compression of −5.4%. Being the same as previous case, the wrinkle starts to move after a temperature gradient is established. Fig. 5 shows the variation of the translational displacement of the wrinkle with time under ΔT/l = 7.05 K nm−1. All the motion data can be well fitted by the quadratic polynomial: . Therefore, soliton-like thermophoresis is also available for a single graphene wrinkle on Cu substrate.
![]() | ||
Fig. 5 Translational displacement of the wrinkle of ε = −5.4% formed in Cu-supported graphene sheet as a function of time under ΔT/l = 7.05 K nm−1. The solid lines are quadratic polynomial fits to the MD data. The inset shows the structure of the wrinkle. |
Furthermore, we have studied thermophoretic motion of two graphene wrinkles on Cu substrate. The two wrinkles are created by applying a uniaxial compression of −10.88% with imposing a constraint in the middle part of the sheet. Then the whole system is relaxed without this constraint under 300 K for 200 ps. The two wrinkles are quite stable in the Cu-supported graphene, as shown in Fig. 6(a). The same simulation setting (Fig. 1) is applied for this system and a temperature gradient of ΔT/l = 3.53 K nm−1 is formed along the sheet before the wrinkles move. Fig. 6(b) shows the displacements of the two wrinkles with time. The wrinkle (W1) in the hotter region comes into motion first, and the second wrinkle (W2) in the cooler region begins to move after an approximate 120 ps delay. Both the W1 and W2 motion data are well fitted by the quadratic polynomial: . The corresponding accelerations for the W1 and W2 are 122 and 112 nm ns−2, respectively. This is because that the size and mass of wrinkle 1 are smaller than that of wrinkle 2. Due to a larger acceleration, the W1 moves faster than the W2 so that they crash with each other after 260 ps. The strong interaction between the crashed wrinkles leads to a structure shift, so the movement of the stacked wrinkles under thermal gradient is different from that of a single wrinkle, and finally they stay at the cooler site. These results indicate that the thermal gradient can drive graphene wrinkles moving to the cooler region and forming a tightly stacking wrinkle structure [the top inset in Fig. 6(b)] on a metal substrate. When the same system is put on a flat graphite substrate, owing to the weaker interlayer interaction it is difficult to stabilize two such wrinkles in the graphene plane, which will transform into one big wrinkle after structural relaxation at 300 K.
![]() | ||
Fig. 6 (a) Two wrinkles are created by applying a uniaxial compression of ε = −10.88% in the x direction of the graphene sheet supported by Cu substrate. (b) Displacements of the two wrinkles W1 and W2 in the x direction as a function of time under ΔT/l = 3.53 K nm−1. The solid lines are the quadratic polynomial fits of ![]() |
This journal is © The Royal Society of Chemistry 2013 |