3.1. The Isothermal Orientation Processes of 79C19/GRA at 450K and 490K
We acquire a greater depth of understanding the ordered structure formation of n-alkane via MD simulations. Our simulations show that C19 molecules are perpendicular to the graphene surface at relative low temperature (T ≤ 480K) and parallel to the graphene surface at relative high temperature (T > 480K). This is similar to our previous work [46–47]. The chain fold also is found for long n-alkane chains in our simulations. Thus, we also show the two typical isothermal orientation processes for 79C19/GRA at 450K and 490K in Fig. 1(a) and (b), respectively.
At the beginning of the simulation, C19 molecules are adsorbed on graphene surface and form disordered structure. At 450K, some C19 molecules are adsorbed near the surface and are parallel to the surface, and others are still disordered at 200ps. Local orientation appears at about 1000ps. After that, more C19 molecules move into the ordered region slowly and the orientation region become larger gradually. Most of C19 molecules are perpendicular to the surface at 2000ps. The perpendicular order degree increases with the simulation proceeding. Chain fold occurs when the ordered structure forms. The chain fold leads to the height of the perpendicular ordered region decreasing gradually. This has not been found when study the orientation process of short n-alkanes on or between graphene sheet(s) [46, 47]. At the end of the simulation, most of the C19 molecules are perpendicular to the surface, and they are also parallel to each other. It can also be seen that the ends of some C19 molecules near the surface is curved in the perpendicular ordered region.
At 490K, some of C19 molecules are also adsorbed near the surface and others are still disordered before 1000ps. Local orientation appears at about 2500ps. Some of C19 molecules are curved and others are extended in the first layer near the graphene surface before 2500ps. After that, more C19 molecules enter the ordered region quickly. At 3000ps, most of C19 molecules are parallel to graphene surface and form a six-layer structure. As the simulation going, the order degree of the six-layer structure increases. The molecules arranging in the six-layer structure are extended and parallel to each other. The effect of graphene on the orientation of C19 molecules results in the orientation starting in the first layer near the graphene surface and then in the next layer. At the end of the simulation, almost all C19 molecules arrange in the six-layer structure.
In summary, the two kinds of orientation processes are considered as a three-step process (adsorption, orientation, and growth). And C19 molecules tend to be perpendicular to graphene surface at relative low temperature and form parallel multi-layer structure at relative high temperature. This is similar to the isothermal orientation of n-decane between two-layer graphene [46]. But the C19 molecules fold in the perpendicular ordered region on graphene surface at relative low temperature. Our simulations also indicate that when the chain length is long enough, chain fold only occurs in perpendicular orientation. We compare the final conformations of n-alkane on graphene surface at different temperatures, the effect of chain length on the orientation of n-alkane on graphene surface can be identified and the critical chain length of n-alkanes with perpendicular and parallel orientation on the graphene surface can also be obtained.
3.2. The Effect of Temperature on the Orientation of n-Alkane on Graphene
The n-alkane chains may be perpendicular orientation, parallel orientation or perpendicular orientation at relative low temperature and parallel orientation at relative high temperature on the graphene surface. Figure 2 displays the final conformations of C13, C19, C25 molecules on graphene surface with the three typical orientations at different temperatures.
For 115C13/GRA, C13 molecules arrange in a disordered state at 300K. Almost all C13 molecules are perpendicular to the surface at 400K, 450K and 490K. The perpendicular order degree increases with temperature increasing. At 490K, some C13 molecules are desorbed from the surface and appear below the surface due to periodic boundary conditions. At 500K, more C13 molecules are desorbed and a disordered structure forms on graphene surface. For 79C19/GRA, a disordered structure form on graphene surface at 300K. C19 molecules are perpendicular to the surface at 400K and 450K. And the height of the perpendicular ordered regions are smaller than the extended C19 molecule shown in Fig. 2(b). The length of extended C19 molecule is 22.59Å. Chain fold occurs when the perpendicular ordered structure of C19 molecules forms. At 490K, almost all C19 molecules are parallel to graphene surface and are parallel to each other. They form a parallel multi-layer structure. At 500K, C19 molecules form a disordered structure on graphene surface. For 60C25/GRA, C25 molecules are disordered at 300K, and a few molecules adsorbed near the surface are parallel to the surface. More C25 molecules are parallel to graphene surface at 400K. At 480K and 500K, parallel multi-layer structures form on graphene surface. Almost all C25 molecules arrange in the ordered structure and are parallel to each other. At 550K, C25 molecules also form a disordered structure on graphene surface.
The n-alkane can form disordered structure, perpendicular ordered structure, parallel multi-layer structure on graphene surface. The final conformations of ten kinds of n-alkane molecules on graphene surface at all the temperatures are shown in figure S1 of the Supporting Information. Ten kinds of n-alkane molecules relax slowly and form disordered structures on graphene surface at 300K. The 115C13, 107C14, 100C15 and 94C16 only form the perpendicular ordered structures on graphene surface. The 88C17, 83C18 and 79C19 are perpendicular to graphene surface at relative low temperature and form parallel multi-layer ordered structure at relative high temperature. The 75C20, 60C25 and 50C30 only form parallel multi-layer structure. When temperature is high enough, the n-alkane molecules form disordered structures on graphene surface. In addition, when n-alkane molecule is longer than C15, chain fold occurs with perpendicular ordered structure forming on the surface. Some n-alkane chains fold into two or more stems during the MD simulation. And n-alkane chains adopt extended conformations in parallel multi-layer ordered structure on the surface.
The chain length affects the orientation of n-alkanes on graphene surface. When the n-alkane molecule is shorter than that of C17, they arrange perpendicularly to the graphene surface. And when n-alkane molecule is longer than C19, they form parallel multi-layer ordered structure on graphene surface. C17, C18 and C19 molecules are perpendicular to graphene surface at relative low temperature and form parallel multi-layer structure at relative high temperature. Thus the longest n-alkane chain that can form perpendicular orientation is C16 molecule on graphene surface. The shortest n-alkane chain that forms parallel multi-layer ordered structure is C20 molecule. And when the chain length is longer than C15, chain fold only occurs in perpendicular orientation region.
Is the orientation of (C17-C19) reversible with temperature i.e., perpendicular orientation at lower temperature and parallel orientation at higher temperature when annealing procedures are performed on them? This is interesting and important in the devising and manufacturing this kind of assembling material. Therefore, annealing procedure is performed on 83C18/GRA at the cooling rate of 10K/2000ps with the same parameters as our isothermal MD simulations. As the orientation of the simulated n-alkanes is sensitive to the cooling rate, another two annealing procedures are simulated at the cooling rate of 10K/1000ps and 10K/500ps. Snapshots of the last conformation of 83C18/GRA for various temperatures during the three annealing procedures are shown in figure S2 of the Supporting Information. The C18 molecules are first disordered at high temperatures and parallel to the surface with temperature decreasing during the annealing procedure at the cooling rate of 10K/2000ps. The C18 molecules first form disordered structure at high temperatures and are perpendicular to the surface with temperature decreasing at the cooling rate of 10K/1000ps and 10K/500ps. The orientation of C18 molecules are irreversible with temperature. The temperature of parallel orientation appearing is higher than that of perpendicular orientation. The C18 molecules are favorite forming perpendicular orientation at fast cooling rate and parallel orientation at slow cooling rate. Only one orientation appears for C18 during the annealing procedure. The parallel and perpendicular orientation can both obtained by isothermal MD simulations at different temperatures. These are useful information for manufacturing this kind of assembling material.
3.3. The Orientation Parameter of n-Alkanes on Graphene
In order to describe the orientation structure of n-alkanes with different chain lengths at different temperatures quantitatively, we calculate the global orientation order parameter g-OPb and the z axis bond-orientation order parameter of the n-alkanes. The global orientation order parameter g-OPb can describe the conformation of all alkane molecules, which is defined as
where φ refers to the angle between two neighboring subbond vectors, the center of the two adjacent bonds are connected to form the subbond vector (see the schematic drawing in Fig. 3). When g-is 1.0, all of the alkane chains are extended. The z axis bond-orientation order parameter is used to describe the orientation of alkane molecules on graphene surface, which is defined as
where ϕ is the angle between the subbond vector and the z axis. When the subbond vector is completely parallel to the z axis, = 1.0; when the subbond vector is perpendicular to the z axis, = -0.5; when the subbond vector is random to the z axis, = 0.0.
Figure 4 and 5 show the global orientation order parameter g-OPb and the z axis bond-orientation order parameter for C13, C19 and C25 molecules on graphene at different temperatures, respectively. For 115C13/GRA, the is about − 0.05 at 300K, indicating that C13 molecules are laid randomly on graphene surface. The increases from 0.7 to 0.9 and the g-OPb fluctuates around 0.95 from 400K to 490K, indicating that more C13 molecules arrange in the perpendicular ordered structure and adopt extended conformation. The is near − 0.1 and the g-OPb is about 0.55 at 500 K, indicating that disorder structure forms on the surface at 500 K.
For 79C19/GRA, the increases from − 0.1 to 0.7 in the temperature range of 300K ~ 480K, meaning that more C19 molecules change from random to perpendicular orientation to the surface. The is -0.4 at 490K, indicating that the C19 molecules are parallel to the surface. The g-OPb fluctuates around 0.73 in the temperature range of 300K ~ 450K, showing that C19 molecules are folding in the ordered structure. The g-OPb increases to 0.88 at 490K. This illustrates that most of C19 molecules adopt extended conformation at 490K. The is 0.0 and the g-OPb is 0.58 at 500 K, indicating that disordered structure forms on the surface at 500 K.
For 60C25/GRA, the g-OPb increases from 0.72 to 0.89 in the temperature range of 300K ~ 500K, which indicates that most of C25 molecules adopt extended conformation, and a small number of molecules are curved. The decreases to -0.4 from 300K to 450K and fluctuates around − 0.4 in the temperature range of 450K ~ 500K, indicating that C25 molecules are parallel to the surface. And the is -0.1 and the g-OPb is 0.55 at 550K as a result of a disordered structure forms on the surface.
Simulation temperature and chain length are important in the orientation of n-alkanes on graphene surface. Figure S3 and S4 show the global orientation order parameter g-OPb and the z axis bond-orientation order parameter for the ten kinds of n-alkane molecules on graphene at different temperatures, respectively. Ten kinds of n-alkane molecules are all randomly on graphene surface at 300K. When the chains are shorter than C17, molecules are perpendicular to the surface. For long chains (> C19), molecules are parallel to the surface. The C17, C18 and C19 molecules are perpendicular to surface at relative low temperature and form parallel multi-layer structure at relative high temperature on surface. When the temperature is high enough, ten kinds of n-alkane molecules all form disordered structure on surface.
3.4. Fold Behavior of n-Alkane Chains on Graphene Surface
In the previous sections, the orientation structure of n-alkane on graphene surface has been discussed. The chain fold happens when long n-alkane form perpendicular ordered structure on graphene surface. This section will explore the chain fold in details. First, we calculate the average end-to-end distance (Red) of n-alkane at different temperatures. It refers to the line distance from one end of the molecule to the other. When the molecular is extended, Red is equal to the length of extended molecule (Rall−trans). When the molecule curved or chain folding occurs, Red is smaller than the Rall−trans.
The average end-to-end distance of the last 1000 configurations of the ten kinds of n-alkane molecules at different temperatures are calculated and provided in figure S5 of the Supporting Information. The average Red of the ten kinds of n-alkane molecules are much less than their Rall−trans at 300K, meaning that molecules are curved at 300K. The average Red of 115C13, 107C14, 100C15, 75C20, 60C25 and 50C30 are close to their Rall−trans with the simulation temperature increasing, meaning that more molecules adopt extended conformation. But the 115C13, 107C14 and 100C15 are extended in the perpendicular ordered structure and the 75C20, 60C25 and 50C30 are extended in the parallel ordered structure on graphene surface. For 94C16/GRA, 88C17/GRA, 83C18/GRA and 79C19/GRA, the average Red at relative low temperature are less than the average Red of disordered structures, meaning that more molecules adopt curved conformation. Chain fold happens during the chains perpendicular orientation on the surface. The average Red is close to the Rall−trans when the four kinds of n-alkane molecules parallel orientation at relative high temperature on graphene surface. This means that most molecules are extended and chain fold does not occur in the parallel ordered structure. The average Red of the ten kinds of n-alkane molecules is less than their Rall−trans and close to the average Red of disordered structures when temperature is high enough, meaning that molecules are curved and disordered on graphene surface. In brief, the occurrence of chain fold is related to the chain length. And chain fold happens when the chains are longer than C15 and form perpendicular ordered structure on graphene surface. This also shows the three kinds of orientations for n-alkanes on graphene surface. The ratio of average Red to Rall−trans (Red/Rall−trans) of three typical systems (115C13/GRA, 79C19/GRA and 60C25/GRA) are shown in Fig. 6 as an example.
The time evolution of the average end-to-end distances during the two orientation process can help us to understand the chain fold process in detail. In order to check the chain fold behavior of 79C19/GRA with perpendicular ordered structure at 450K and parallel multi-layer structure 490K, Fig. 7 shows the time evolution of the average end-to-end distances of 79C19/GRA on graphene surface at 450K and 490K.
At 450K, the average Red of 79C19 increases to 17.0Å before 600ps, decreases to 10.3Å from 600ps to 8000ps gradually, and fluctuates around 10.3Å after 8000ps, which are smaller than their Rall−trans (22.59Å). The ratio of average Red to the film thickness of 79C19 is about 0.72. This means that the formation of perpendicular orientation is accompanied by chain fold. The height of the perpendicular ordered regions are smaller than the extended C19 molecule. It may be speculated that the height of the perpendicular ordered regions decreases owing to the chain fold during the simulation proceeding. At 490K, the average Red of 79C19 fluctuates around 15.8Å before 2000ps, increases to 20.0Å from 2000ps to 5000ps, and fluctuates around 20.0Å after 5000ps, which are close to their Rall−trans. This shows that most of C19 molecules change from curved conformation to extended conformation gradually during the orientation process. The C19 molecules do not fold in the parallel ordered structure.
The C-C-C-C torsion angle of C19 molecules are also calculated to describe the conformation at 450K and 490K. Figure 8 gives the torsion angle distributions calculated from the last 1000ps of the simulations. There are two kinds of peaks in the distribution. The first one is the torsion angle between 30° ~ 100° or -30° ~ -100°. The peak indicates the presence of the gauche conformation. The second one is the torsion angle between 150° ~ 180° or -150° ~ -180°. The peak means that the presence of the trans conformation. The second peak is much higher than the first one. This means most of C19 molecules adopt trans conformation. For the peak between 30° ~ 100°or -30° ~ -100°, the height at 450K is higher than that at 490K. This indicates that more C19 molecules are gauche conformation at 450K. This results from the occurrence of chain fold in the perpendicular ordered structure at 450K.
There are 16 torsion angles in C19 molecule, which are denoted in turn. Figure 9 shows the final conformation of 79C19/GRA at 450K. Five typical C19 molecules are selected to study, which denoted by A (blue), B (orange), C (green), D (purple) and E (red). It can be seen that C19 A is extended, C19 B and C fold into two stems, C19 D occurs a non-integer fold, and C19 E folds into three stems. The torsion angles of C19 A, B, C, D and E at 10,000ps NVT MD simulations at 450K are listed in Table 1. The negative values of torsion angles in the Table 1 are shown as their absolute value.
Torsion angles of C19 A are all bigger than 150°, which indicates that C19 A is extended on graphene surface. There are fold and extended regions in the fold C19 molecules (B, C, D and E). The extended region is trans conformation. The torsion angles 1 ~ 6, 12 ~ 16 of C19 B, 2 ~ 7, 13 ~ 16 of C19 C, 1 ~ 8 of C19 D are about 170°. The conformation of the fold region is the gauche-gauche-trans-gauche-gauche or gauche-gauche-trans-gauche. This is similar with the structure of chain fold in the polyethylene [25, 57]. The torsion angles 7 ~ 11 of C19 B, 8 ~ 12 of C19 C, 9 ~ 12 of C19 D are in the fold regions. C19 E folds into three stems in Fig. 9. Torsion angles 2, 3, 5, 6, 8, 11, 12, 15 of C19 E are all less than 100°, indicating that C19 E twists heavily.
Liu et al. [33] found more detailed information of chain fold during n-alkane orientation from the evolution of the end-to-end distance. The change of C19 can be found from the time evolution of Red. The time evolution of the Red of C19 A, B, C, D and E at 450K are also shown in Fig. 10. The C19 A are curved at the beginning of the simulation and changes to extended conformation after 2000ps. The C19 B is bent from 1900ps to 8000ps and then folds into two stems again after 8000ps. The C19 C happens to fold before 2000ps and folds into two stems from 2000 to 10000ps. The C19 D folds into two stems at 5000ps and changes to a non-integer fold at the end of simulation. The C19 E twists at about 600 and 3200ps and fold into two stems at 4000ps. The C19 E folds into three stems finally.

When the conformation of n-alkane chain changes, the backbone torsion angles will change correspondingly. We choose three typical backbone torsion angles 1, 6, 8 of C19 B, which are at the end of chain, in the extended stem and in the fold region of the molecule after 10,000 NVT MD simulations, respectively. Figure 11 shows the evolutions of torsion angles 1, 6, 8 of C19 B on graphene surface during MD simulation at 450 K.
Three torsion angles fluctuate between 20° and 180° during C19 B twisting. The C19 B bends after 1900ps. Torsion angle 6 fluctuates between 150° and 180° from 1900ps to 4700ps, meaning that torsion angle 6 is in the extended stem. Torsion angles 1 and 8 fluctuate wildly from 1900ps to 4700ps. In the simulated time of 4700ps ~ 8000ps, torsion angles 1 and 6 fluctuate between 150° and 180°, meaning that they are in the extended stem. Torsion angle 8 fluctuates between 40° and 80°, meaning that torsion angle 8 is in the fold region. C19 B happens to fold at 8000ps and folds into two stems from 8000 to 10000ps again. Three torsion angles don’t change much after 8000ps. And torsion angles 1 and 6 are in the extended stem and torsion angle 8 is in the fold region.
When n-alkane chain fold occurs, torsion angle may move back and forth in the fold region and the extended stem. Some torsion angles in the fold region usually change to gauche conformation, and others in the extended stem usually adopt trans conformation. After the occurrence of chain fold, torsion angles fluctuate. Torsion angles in the fold region fluctuate stronger than that in the extended stem.
3.5. The Effect of the Interaction Energy on the Orientation
In this section, we will explore the effect of the interaction energy on the orientation of n-alkane with different chain lengths on graphene surface. We calculate the interaction energy between n-alkane chain (Echain), the interaction energy between molecule and graphene surface (Eint).
Eint = Etotal -(Echain + Eplane) (3)
where Etotal is total energy of the system, Eplane is energy of graphene surface. Figure 12 shows the variation of Eint, Echain and Eint-Echain for 115C13/GRA, 79C19/GRA and 60C25/GRA with increasing temperature.
For 115C13/GRA, Eint changes little as the temperature increases. Echain decreases from 300 to 450K and increases from 450 to 500K quickly. Eint is bigger than Echain at 300K, and Eint-Echain is about 2000kcal/mol. C13 molecules are disordered on the surface. This is due to n-alkane relaxing slowly at low temperatures. Echain are much lower than Eint, and Eint-Echain is about 4000kcal/mol from 370 to 490K. C13 molecules are perpendicular to graphene surface in this temperature range. Echain is as much as Eint at 500K when C13 molecules form a disordered structure.
For 79C19/GRA, Echain is lower than Eint, and Eint-Echain is about 2000kcal/mol at 300K. C19 molecules are disordered on the surface because n-alkane relaxes slowly. Echain decreases and Eint increases slowly from 300 to 480K. In addition, Echain is much lower than Eint, and Eint-Echain increases to about 3000kcal/mol gradually. C19 molecules tend to adopt perpendicular orientation on surface. When the temperature above 480K, Echain increases and Eint decreases suddenly. Echain is lower than Eint and Eint-Echain decreases to about 2000kcal/mol at 490K. C19 molecules form parallel ordered structure on the surface at 490K. Echain is close to Eint and Eint-Echain is about 500kcal/mol at 500K. C19 molecules form disordered structure on the surface.
For 60C25/GRA, Echain decreases from 300 to 480K, increases after 480K quickly. Eint changes little from 300 to 480K, decreases after 480K suddenly, and increases from 500 to 550K. A disordered structure form on the surface at 300K owing to n-alkane relaxes slowly. And Eint-Echain fluctuates around 2000kcal/mol, C25 molecules are parallel to the surface in the temperature range of 400 ~ 500K. C25 molecules form a disordered structure on the surface at 550K when Echain is close to Eint and Eint-Echain is about 500kcal/mol.
It is found that the ordered structure of n-alkanes is governed by Eint and Echain. When Eint-Echain is high enough at relative high temperature, n-alkanes can form ordered structure on graphene surface. Molecules are perpendicular to the surface when Eint-Echain of C13 molecules is about 4000kcal/mol. C19 and C25 molecules are parallel to the surface when Eint-Echain is about 2000kcal/mol. Although Eint-Echain is about 2000kcal/mol at low temperatures, C13, C19, C25 molecules are disordered on the surface. It's because n-alkane molecule relaxes slowly at low temperatures. When Eint-Echain is about 500kcal/mol at high temperatures, three kinds of n-alkanes form disordered structure. Eint-Echain of parallel ordered structure is much lower than that of perpendicular ordered structure. And Eint-Echain of disordered structure is lower than that of parallel ordered structure on the surface. This is consistent with the results of previous studies [46, 47, 58]. However, C19 molecules form perpendicular ordered structure when Eint-Echain is less than 4000kcal/mol because of chain fold.
As the simulation temperature changes, Eint and Echain change correspondingly and n-alkane molecules form different orientation structures. The difference between Eint and Echain of ten kinds of n-alkane molecules on graphene surface at different temperatures are shown in table S1 of the Supporting Information. Only when Eint-Echain is high enough, n-alkane molecules can form ordered structure on graphene surface. In addition, even if Eint-Echain is high, the ordered structure can’t form when the temperature is too low. When Eint-Echain is about 4000kcal/mol, short n-alkanes are perpendicular to the surface. When Eint-Echain is about 2000kcal/mol, n-alkanes are parallel to the surface. But when long n-alkanes form perpendicular ordered structure on the surface, Eint-Echain is less than 4000kcal/mol due to chain fold. When Eint-Echain is small, n-alkanes are disordered on the surface.