Next Article in Journal
Memory Effect on a LDH/zeolite A Composite: An XRD In Situ Study
Next Article in Special Issue
DABCO Derived Nitrogen-Doped Carbon Nanotubes for Oxygen Reduction Reaction (ORR) and Removal of Hexavalent Chromium from Contaminated Water
Previous Article in Journal
Effect of Pre-Stretch on the Precipitation Behavior and the Mechanical Properties of 2219 Al Alloy
Previous Article in Special Issue
Effects of Boron Carbide on Coking Behavior and Chemical Structure of High Volatile Coking Coal during Carbonization
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Crumpled Graphene-Storage Media for Hydrogen and Metal Nanoclusters

by
Liliya R. Safina
1,*,
Karina A. Krylova
2,3,
Ramil T. Murzaev
2,
Julia A. Baimova
2,3 and
Radik R. Mulyukov
1,2
1
Ufa State Petroleum Technological University, Kosmonavtov Str. 1, 450062 Ufa, Russia
2
Institute for Metals Superplasticity Problems of the Russian Academy of Sciences, Khalturina 39, 450001 Ufa, Russia
3
Bashkir State University, Validy Str. 32, 450076 Ufa, Russia
*
Author to whom correspondence should be addressed.
Materials 2021, 14(9), 2098; https://doi.org/10.3390/ma14092098
Submission received: 17 March 2021 / Revised: 5 April 2021 / Accepted: 14 April 2021 / Published: 21 April 2021
(This article belongs to the Special Issue Carbon-Based Materials)

Abstract

:
Understanding the structural behavior of graphene flake, which is the structural unit of bulk crumpled graphene, is of high importance, especially when it is in contact with the other types of atoms. In the present work, crumpled graphene is considered as storage media for two types of nanoclusters—nickel and hydrogen. Crumpled graphene consists of crumpled graphene flakes bonded by weak van der Waals forces and can be considered an excellent container for different atoms. Molecular dynamics simulation is used to study the behavior of the graphene flake filled with the nickel nanocluster or hydrogen molecules. The simulation results reveal that graphene flake can be considered a perfect container for metal nanocluster since graphene can easily cover it. Hydrogen molecules can be stored on graphene flake at 77 K, however, the amount of hydrogen is low. Thus, additional treatment is required to increase the amount of stored hydrogen. Remarkably, the size dependence of the structural behavior of the graphene flake filled with both nickel and hydrogen atoms is found. The size of the filling cluster should be chosen in comparison with the specific surface area of graphene flake.

1. Introduction

The transformation of nanoscale structural elements into three-dimensional (3D) complex architecture is currently an important task of materials science. Since the discovery of fullerenes in 1985, a lot of new carbon structures have been proposed. Carbon polymorphs can be used to obtain nanostructures with unique mechanical and physical properties applicable in nanoelectronics, energy storage devices, sensors, supercapacitors, Li-ion batteries, etc. [1,2,3,4,5,6,7].
Graphene is capable of enhancing the performance, functionality as well as durability of many applications. The one-atom thin structure can serve as the platform for other materials especially since the graphene layer can be bent or crumpled. Extensive studies have been carried out in recent decades to investigate the crumpling behavior of thin sheets like graphene, by both theoretical and experimental methods [8,9,10,11,12,13,14,15,16,17,18,19]. It was shown that crumpled structures can have excellent compression and aggregation-resistant properties [5,6,9,10,20,21,22]. Crumpled graphene can be used for the production of new supercapacitors [19]. However, one of their new and important applications is that of a natural container for other atoms.
It is known that coupling metal nanoclusters with carbon materials can efficiently promote catalysis and electrocatalytic activities [23,24] or prevent the corrosion of metals [25]. redAmong metals, some can be easily attached to its surface, while others even repulse from the surface. The bigger the carbon solubility of metal, the more amounts of carbon can be attracted to the metal surface. Metals such as Ni, Cu, Pt, or Au are the most preferable as catalysts for the growth of graphene [26,27,28]. However, nickel is one of the most broadly used metals for a range of graphene applications, including the growth of carbon nanotubes (CNTs) and graphene [29,30,31]. Lattice mismatch is one of the important factors for choosing the metal for metal–graphene interface and Ni(111) surface is the closest matched interface with respect to the graphene of all transition metals [26]. In [29], a simple and scalable method for the synthesis of hollow graphene balls using a Ni nanocluster template was developed. The interface between graphene and Ni attracts considerable interest due to the possibility of the synthesis of large-area graphene on metal substrate [32,33,34,35,36]. Since the fabrication of graphene layers on the metal substrate was very successful, the idea of the fabrication of metal–graphene layered composites was raised by [37,38,39,40]. Graphene added into the copper and nickel matrix by chemical vapor deposition can considerably improve the strength of metal [38]. The overview of graphene–nickel composites can be found in [39]. Both the theoretical and experimental works on Ni-graphene layered composites showed that the mechanical properties of metal-matrix composite reinforced with the graphene layer considerably improved. In the present work, the other approach to obtain the Ni-graphene composite was applied: introducing metal nanoparticles to the crumpled graphene matrix. In [41,42,43,44] it was shown that such composites can be obtained through compression combined with high temperatures.
The other important advantages of such structures are a large specific surface area (SSA) and, respectively, a high rate of gas adsorption, which makes it possible to predict their use in hydrogen storage [45,46,47,48,49,50,51,52,53]. Various structural parameters can affect the degree of hydrogen accumulation [45]. For example, the amount of adsorbed hydrogen increases with an increase in the CNT diameter, since this increases the surface area on which hydrogen can be adsorbed. An increase in the distance between graphene sheets in the structure also leads to an increase in the gravimetric hydrogen absorption density. Another promising carbon structure is two layers of graphene connected by short CNTs with graphene cones added on top of the structure. The idea of such a structure appeared when it was shown that the hydrogenation of graphene can lead to the buckling of the graphene sheet [54]: the amount of accumulated hydrogen increased from only 3 wt% for the simple graphene plane to 20 wt% for graphene plane with graphene cone. Carbon nanostructures doped with alkali metals (Li and K) can adsorb even more hydrogen [48,49,50,51,52,53]: 14 wt% (Li) and 20 wt% (K) of hydrogen at moderate conditions, in contradiction with the lower values reported later [50]. The Li-doped activated carbon [55] can store from 2.1 to 2.6 wt% of H 2 at 77 K and at 2 MPa, which shows that at given pressure–temperature conditions the amount of stored hydrogen can be increased. It should be mentioned that the required amount of stored hydrogen for carbon nanostructures was still not achieved. Thus, an active search for new materials and structures for hydrogen storage and transportation is of high importance.
In the present work, two types of nanofillers for crumpled graphene flake were considered—Ni and hydrogen nanoclusters. Thus, the idea of using crumpled graphene flake for the storage of different atoms using Ni and hydrogen clusters as the example was realized. It is important to understand the dynamics of interaction between single graphene flake and a metal or non-metal nanocluster before the study of storage in three-dimensional crumpled graphene. The behavior of the hybrid structure is studied by atomistic simulation at different temperatures. For the carbon–nickel system, room temperature (300 K) and temperature close to the melting point of Ni (1000 K) were chosen. However, for the carbon–hydrogen system, the temperature of liquid nitrogen (77 K) and room temperature were chosen. The dynamics of the structure were studied during exposure for 20 ps.

2. Materials and Methods

To study the interaction between graphene flake (GF) and atomic nanocluster, an initial structure composed of the GF of one size ( N C = 252) and atomic clusters of different sizes (N = 21; 38; 47; 66; 78) inside the flake was considered (see Figure 1a). One graphene flake was a small CNT(11,11) 1.3 nm-long with two atomic rows deleted along the axis of the nanotube. Since the size of this small sample of graphene along three dimensions is very small (1.3 nm along the z axis, 3.6 nm along the y axis and one-atom-thick along the x axis) it is called “flake” rather than “nanoribbon”. This GF was filled with nanoclusters of two types: (i) nickel nanocluster and (ii) hydrogen nanocluster.
The maximum size of the nanocluster was chosen to almost completely fill the cavity of the graphene flake. The size of the structural elements is shown in Table 1. The distance between the edge of the nanocluster and the side of the flake can be defined as a = ( D d ) / 2 . However, these parameters affect the structural transformation just in the case of metal nanocluster since it is a solid particle. a hydrogen cluster is gas and small H atoms can easily spread inside GF.
It should be mentioned, that nanocluster Ni 78 was considered here as nanocluster of critical size for the chosen size of the GF. However, for the metal cluster, this size was too big and the sides of GF were closer than it should have been in a real system. To eliminate the negative effect of the overly sized nanocluster, nanocluster Ni 66 was considered the characteristic for the case of metal–graphene interaction. While for the case of hydrogen–graphene interaction, an even bigger size of the nanocluster could be considered since the size of the hydrogen atoms was much smaller than the size of Ni atoms.
The simulation was conducted using a large-scale atomic/molecular massively parallel simulator (LAMMPS). Equations of motion for the atoms were integrated numerically using the fourth-order Verlet method with the time step of 0.2 fs. The Nose–Hoover thermostat was used to control the system temperature. The periodic boundary conditions were applied in all directions, however, the simulation box is much bigger than the size of GF filled with nanocluster. The adaptive intermolecular reactive empirical bond order potential (AIREBO) [56] was used to describe the interatomic interactions between carbon atoms, including both covalent bonds in the basal plane of graphene and van der Waals interactions between GF and nanocluster. The simulation configurations were visualized by Visual Molecular Dynamics (VMD) Software [57].

2.1. Graphene with Nickel Nanocluster

Graphene flake with the nanocluster was exposed at 300 K and 1000 K for 20 ps to study the dynamics of nanocluster inside graphene flake. Previously [41], it was shown that the melting temperature of the Ni nanocluster was about 1300 K. Thus, in the present work, the highest considered temperature was 1000 K. At this temperature, Ni nanocluster was close to melting but not melted yet, even for the smallest nanocluster, which consisted of 21 atoms. Initially, Ni atoms were packed into the face-centered cubic lattice.
To describe the interatomic interaction between Ni–Ni and Ni–C, Morse interatomic potential was used with the parameters D e = 0.4205 eV, R e = 2.78 Å and β = 1.4199 1/Å for Ni–Ni [58]; and D e = 0.433 eV, R e = 2.316 Å, β = 3.244 1/Å for Ni–C. The parameters of the Morse potential for describing the interaction of nickel and carbon atoms were recently proposed using ab-initio simulation [59,60].

2.2. Graphene with Hydrogen Nanocluster

Graphene flake with hydrogen nanocluster was exposed at 77 K and 300 K for 20 ps. The temperature of 77 K was chosen since it was previously shown that better sorption of hydrogen can be found at 77 K [61]. The temperature of 300 K was chosen to study the dynamics at room temperature, however, it is known that this temperature is too high for hydrogen storage [46].
To describe the interatomic interaction between C and H, AIREBO interatomic potential was used.
It should be mentioned that the initially obtained structure contains atomic hydrogen, however, H atoms transform to H 2 molecules and just several H atoms remain single (see Figure 1b). As it was shown, the lowest binding energy between graphene and H 2 was observed when the distance between C and H 2 in the range from 2.9 to 3.2 Å [61]. Hydrogen molecules can be bonded by van der Waals interaction when they move close enough to the side of GF. Single hydrogen atoms can be chemically adsorbed on graphene by covalent bonding.

3. Results

3.1. Graphene with Ni Nanocluster

At first, the evolution of the potential energy of the system during the crumpling process for graphene flake filled within spherical the Ni nanocluster was analyzed. In Figure 2, potential energy as the function of exposure time at 300 K (a) and 1000 K (b) for five types of nanoclusters was shown. It was found that the total potential energy of the system was saturated to a practically constant value at the end of the equilibration process, indicating that the system reached equilibrium and a stable state. All the changes in the energy curves correspond to some structural changes.
The longest time of stabilization at 300 K was found for Ni 21 (5 ps). Curves for nanoclusters of close diameter are almost the same. Graphene flake with Ni 38 and Ni 47 reach the equilibrium state at about 4 ps of exposure, while the GF with Ni 66 and Ni 78 reach equilibrium state at 3.2 ps. The bigger the diameter of the nanocluster, the less the equilibration time. This can be explained by the mutual arrangement of the nanocluster and GF. For a small nanocluster, the distance a is two times higher than for the biggest one which means that the time required to attach the nanocluster by the graphene flake is longer. For nanocluster Ni 78 , 3 ps is enough for GF to fully cover the nanocluster, while for Ni 21 , not only coverage took place, but also the further crumpling of the flake with the changing of the round shape of nanocluster.
At higher temperature (1000 K), again a strong correlation between equilibration time and the size of the nanocluster was found. For GF with Ni 21 nanocluster, the transformation was fast since the temperature was close to the melting point and the nanocluster can be easily destroyed. Temperature fluctuations facilitate the crumpling process and the appearance of new bonds between the edges of GF. For bigger nanoclusters, the temperature slightly affects the time of equilibration. Again, the same values of equilibration time were obtained for Ni 38 and Ni 47 (about 4 ps) and for Ni 66 and Ni 78 (about 3.2 ps). However, temperature decreases the total potential energy of the system.
In Figure 3, the process of crumpling for nanoclusters Ni 21 (a), Ni 47 (b) and Ni 66 (c) was shown in details after exposure at 300 K (blue line) and 1000 K (red line). At the beginning of the crumpling process, GF startED to change its round shape, and the fcc crystalline order of Ni nanocluster WAs destroyed. Metal atoms are attracted by the graphene surface and tend to occupy equilibrium positions above the center of the carbon hexagon, which was also mentioned in [41]. Ni atoms WEre interacting with graphene hexagons by van der Waals forces. The edges of GF can be bonded during exposure.
At first, consider GF with Ni 21 (see Figure 3a, lower line of snapshots). At 300 K, GF lost its round shape at about 0.5 ps, edges of GF move towards each other, and at t = 1.5 ps one covalent bond appeared between two edges. Nanocluster disturbed by the temperature and several atoms attach to the graphene plane (state II). This leads to the spreading of the Ni atoms over the graphene plane (state III). Simultaneously, such a small nanocluster allowed the graphene flake to easily bend and the structural unit transforms to the bi-layer graphene with Ni atoms spread inside (state IV). Several more covalent bonds appeared on the side edges of GF. At 300 K, not many covalent bonds between the edges of GF can appear, since the temperature is far from the melting temperature of graphene [62,63].
Increase in the exposure temperature to 1000 K facilitates the transformation. At 300 K, equilibration time is equal to 5 ps, while at 1000 K, it is equal to 3.5 ps. Subsquently, the structure came to a stable state and no further changes were observed despite slight thermal fluctuations. At 1000 K, more carbon atoms on the edges of GF found neighbors since the structure is disturbed by thermal fluctuations.
For structures with Ni 47 and Ni 66 , the behavior is qualitatively close. An almost round Ni nanocluster inside the graphene flake was observed at the initial state at both 300 K and 1000 K (as can be seen in Figure 3b,c). The edges of GF can attach to each other with the formation of new covalent bonds. During exposure, graphene flake transforms into a capsule containing nickel nanocluster. The nickel nanocluster almost completely fills the graphene flake, in comparison with Ni 21 . As a result, it is difficult to deform such a structure. However, at 1000 K, the nanocluster is more disturbed and the stable state (state III, the top line of snapshots) reached at a lower equilibration time. The nickel nanocluster became more planar since GF is rigid and can bend the soft metal nanocluster. In the case of Ni 47 , GF transforms the “bag” for nickel, while Ni 66 nanocluster is too big and edges on both sides cannot be attached. In the case of Ni 66 , GF just covers the Ni nanocluster as much is possible.
A graphene flake with Ni 38 behaves in a similar way to Ni 47 , while GF with Ni 78 behaves similarly to Ni 66 . Note that GF always tends to wrap the metal cluster.
One of the important applications of such structural elements is the fabrication of composites. In Figure 4, the initial structure (a, a’) was composed of graphene flake filled with Ni nanocluster Ni 21 and Ni 78 correspondingly. The initial structure is compressed at 1000 K to obtain graphene–nickel composite material (c, c’). Composite can be fabricated under hydrostatic compression at high temperatures [41,42,43,44].

3.2. Graphene with Hydrogen Nanocluster

In Figure 5, the potential energy of GF filled with hydrogen cluster during exposure at 77 K (a) and 300 K (b) is presented as the function of equilibration time. Similar to GF with the Ni nanocluster, five structural units were divided into three groups—GF with 21 hydrogen atoms; GF with 38 and 47 H atoms; and GF with 66 and 78 H atoms. Despite that the hydrogen atoms transform into hydrogen molecules, the initial number of atoms was also used, since the number of H 2 molecules and single hydrogen atoms can change from one simulation run to another which depends on thermal fluctuations.
At 77 K (Figure 5a), the equilibrium state is the state when hydrogen atoms found their sits and GF change the cylinder shape to the one with minimal energy. At 300 K (Figure 5a), the equilibrium state is the state when all hydrogen atoms disrobed from the side of GF or even fully leave the cavity of GF. This would be discussed later together with the description of corresponding snapshots. The biggest drop of the energy value took place during the first picosecond and connected with the change of initial nonequilibrium configuration of the cluster and slightly with the change of shape of the GF.
At 77 K, the time of equilibration was the longest (3 ps) for hydrogen clusters consist of 77, 66, and 38 atoms, and the shortest (1.5 ps) is for 21 and 47 atoms. At that temperature, the size of the cluster plays quite an important role which is connected with the number of sites for hydrogen on the side of GF. All hydrogen molecules and atoms can easily find sites on graphene for an initial number of atoms less than 47, while for bigger clusters, the number of sites is not enough and some molecules and atoms will move outside GF or seek better sits near GF.
At 300 K, the main point is the process of dehydrogenation which is quite quick at such a high temperature, and for clusters with 38–78 atoms, about 1.5 ps is enough for all molecules to detach the side of GF and move outside. The long time of equilibration for cluster H 21 can be explained simply: there are about nine molecules and three atoms and all the time spent for the slow movement of this hydrogen. A large number of atoms in the cluster pushes the sides of GF and opens it much faster.
To understand the dynamics of the hydrogen sorption/desorption, snapshots of the structural units altogether with their potential energy are presented in Figure 6 for three groups. If the hydrogen cluster is quite small (21 atoms), there are a lot of vacant places on the graphene plane to attach hydrogen. Only several atoms or molecules can move outside GF at a low temperature equal to 77 K. Most of the hydrogen molecules attached by van der Waals force to graphene and slightly moving along graphene flake.
At 300 K, hydrogen molecules and atoms have no chance to form even a weak bond to graphene. At 300 K, all hydrogen atoms move outside GF which is quite understandable. Numerous theoretical and experimental works confirm the sorption of molecular hydrogen on carbon nanostructures [50,64,65,66] in a specific temperature range (50–200 K). Here, such a high temperature was used to analyze how hydrogen will leave the cavity of GF. As it can be seen from Figure 6a, this process is very fast and at 4 ps, only two molecules stay inside GF (stage IV, upper line).
A large number of hydrogen atoms almost completely fills the cavity of the GF. At 77 K, hydrogen atoms are easily converted into hydrogen molecules and under the action of chemical attraction or van der Waals forces, attach to the walls of the graphene flake. Mainly, hydrogen, which is located in the center of the flake under the influence of temperature, tends to leave the graphene cavity. Therefore, in structures with 66 and 78 hydrogen atoms, the stabilization of the structure takes longer (3 ps), in contrast to small hydrogen clusters, such as 21 and 47 (1.5 ps).
For a bigger nanocluster, energy curves almost coincide for two temperatures, since hydrogen molecules at 77 K can find the places for sorption just after the first steps. In this case, hydrogen atoms are placed near graphene and link it as far as the van der Waals radius of the molecule reaches. Thus, the hydrogen adsorption capacities depend considerably on the initial size of the hydrogen cluster. The SSA of the GF should be big enough for a chosen number of H atoms and H 2 molecules. Here, SSA is equal to 1153.9 m/g 2 , which is enough to settle down 21–38 H atoms, but not enough for a bigger number of H atoms.
In Figure 7, snapshots of GF with 78 hydrogen atoms during exposure at 77 K for 20 ps are presented. As it can be seen, some hydrogen molecules attach the opposite side of GF since there were no sites inside the cavity. During exposure, even at 77 K, GF opens and then closes again. If the exposure time would be increased to even 100 ps, GF will move like the wings of a flying bird with adsorbed hydrogen atoms. This state is equilibrium and can be preserved for a long time.
As well as for a metal nanocluster, the three-dimensional structure of crumpled graphene filled with hydrogen is presented in Figure 8. The corresponding structure was considered in [46,47], where the effect of hydrostatic compression on the possibility of hydrogen storage was studied. In Figure 8a, the initial structure of crumpled graphene was presented. As it can be seen, initially there are too many channels for hydrogen to move out of the structure. Thus, additional treatment was required to save hydrogen inside the pores of crumple graphene. In Figure 8b, the structure after 40% of hydrostatic compression was presented. In such a compressed structure, hydrogen can be stored much more effectively than in undeformed crumpled graphene [46,47]. It can be concluded that crumpled graphene is an effective storage media for hydrogen. However, a search for the improvement of the quantity of stored hydrogen should be found, for instance, in the introduction of metal atoms. From this point of view, understanding the interaction between hydrogen and metal nanoclusters and GF is of high importance.

4. Conclusions

Molecular dynamics simulation is used to study the dynamics of graphene flakes filled with the nanoclusters of two different types: metal and non-metal. The obtained results can shed the light on understanding the possibility of using crumpled graphene as a storage media from the point of single flake behavior.
It was found that cavities of crumpled graphene (the structure consists of crumpled graphene flakes) can be used as containers for metal nanoclusters—for instance, for Ni. A nanocluster consisting of 21 to 78 Ni atoms was considered. A graphene flake can easily cover nanocluster, however, the dynamics of the interaction strongly depend on the nanocluster size. Small nanoclusters can be easily bent by rigid graphene flake, while the biggest conserve the shape. Such a structure, composed of Ni nanoclusters and graphene flakes, can be further used to obtain composite material with improved mechanical properties.
The problem of hydrogen storage has been of high importance for decades and the application of carbon nanostructures from this point of view also looks promising. Crumpled graphene has a high specific surface area, light weight and a lot of pores and cavities which can be filled with hydrogen. It was shown that at 77 K, hydrogen molecules and atoms can be absorbed by both sides of graphene flakes. However, a single flake cannot store enough hydrogen for practical application. However, when graphene flakes composed in another structure and with additional treatment like hydrostatic compression [46,47], it can be successfully used for hydrogen storage and transportation.

Author Contributions

Software and initial simulation setup were conducted by L.R.S. and R.T.M. The molecular dynamics study and drafting the manuscript were conducted by J.A.B., K.A.K. and R.R.M. All authors have read and agreed to the published version of the manuscript.

Funding

Part of this research was supported by the grant of Russian Science Foundation (No. 20-72-10112). Work of R.R.M. and R.T.M. supported by the program of fundamental researches of Government Academy of Sciences of IMSP RAS.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
MDMolecular Dynamics
CNTCarbon Nanotubes
GFGraphene Flake
SSASpecific Surface Area

References

  1. Iwan, A.; Malinowski, M.; Pasciak, G. Polymer fuel cell components modified by graphene: Electrodes, electrolytes and bipolar plates. Renew. Sustain. Energy Rev. 2015, 49, 954–967. [Google Scholar] [CrossRef]
  2. Baimova, J.A.; Rysaeva, L.K.; Liu, B.; Dmitriev, S.V.; Zhou, K. From flat graphene to bulk carbon nanostructures. Phys. Status Solidi (b) 2015, 252, 1502–1507. [Google Scholar] [CrossRef]
  3. Xiao, J.; Mei, D.; Li, X.; Xu, W.; Wang, D.; Graff, G.L.; Bennett, W.D.; Nie, Z.; Saraf, L.V.; Aksay, I.A.; et al. Hierarchically Porous Graphene as a Lithium–Air Battery Electrode. Nano Lett. 2011, 11, 5071–5078. [Google Scholar] [CrossRef] [PubMed]
  4. Liu, T.; Zhang, L.; Cheng, B.; Hu, X.; Yu, J. Holey Graphene for Electrochemical Energy Storage. Cell Rep. Phys. Sci. 2020, 1, 100215. [Google Scholar] [CrossRef]
  5. Baimova, Y.A.; Murzaev, R.T.; Dmitriev, S.V. Mechanical properties of bulk carbon nanomaterials. Phys. Solid State 2014, 56, 2010–2016. [Google Scholar] [CrossRef]
  6. Niu, L.; Xie, J.; Chen, P.; Li, G.; Zhang, X. Quasi-static compression properties of graphene aerogel. Diam. Relat. Mater. 2021, 111, 108225. [Google Scholar] [CrossRef]
  7. Yola, M.L. Development of Novel Nanocomposites Based on Graphene/Graphene Oxide and Electrochemical Sensor Applications. Curr. Anal. Chem. 2019, 15, 159–165. [Google Scholar] [CrossRef]
  8. Meyer, J.C.; Geim, A.K.; Katsnelson, M.I.; Novoselov, K.S.; Booth, T.J.; Roth, S. The structure of suspended graphene sheets. Nature 2007, 446, 60–63. [Google Scholar] [CrossRef]
  9. Tallinen, T.; Åström, J.A.; Timonen, J. The effect of plasticity in crumpling of thin sheets. Nat. Mater. 2008, 8, 25–29. [Google Scholar] [CrossRef]
  10. Matan, K.; Williams, R.B.; Witten, T.A.; Nagel, S.R. Crumpling a Thin Sheet. Phys. Rev. Lett. 2002, 88. [Google Scholar] [CrossRef] [Green Version]
  11. Zhang, L.; Zhang, F.; Yang, X.; Long, G.; Wu, Y.; Zhang, T.; Leng, K.; Huang, Y.; Ma, Y.; Yu, A.; et al. Porous 3D graphene-based bulk materials with exceptional high surface area and excellent conductivity for supercapacitors. Sci. Rep. 2013, 3. [Google Scholar] [CrossRef] [Green Version]
  12. Balankin, A.S.; Ochoa, D.S.; León, E.P.; de Oca, R.C.M.; Rangel, A.H.; Cruz, M.Á.M. Power law scaling of lateral deformations with universal Poisson’s index for randomly folded thin sheets. Phys. Rev. B 2008, 77. [Google Scholar] [CrossRef]
  13. Baimova, J.A.; Liu, B.; Zhou, K. Folding and crumpling of graphene under biaxial compression. Lett. Mater. 2014, 4, 96–99. [Google Scholar] [CrossRef] [Green Version]
  14. Krainyukova, N.V.; Zubarev, E.N. Carbon Honeycomb High Capacity Storage for Gaseous and Liquid Species. Phys. Rev. Lett. 2016, 116. [Google Scholar] [CrossRef]
  15. Semat, J.M.; Rashmi, W.; Mahesh, V.; Khalid, M.; Priyanka, J. Synthesis of crumpled graphene by fast cooling method. In Proceedings of the International Engineering Research Conference-12th Eureca 2019; AIP Publishing: Melville, NY, USA, 2019. [Google Scholar] [CrossRef]
  16. Savin, A.; Savina, O.I. The Effect of Layers Interaction on the Stiffness of Bending Deformations of Multilayered Carbon Nanoribbons. Phys. Solid State 2019, 61, 686–692. [Google Scholar] [CrossRef]
  17. van Bruggen, E.; van der Linden, E.; Habibi, M. Tailoring relaxation dynamics and mechanical memory of crumpled materials by friction and ductility. Soft Matter 2019, 15, 1633–1639. [Google Scholar] [CrossRef]
  18. Liao, Y.; Li, Z.; Xia, W. Size-dependent structural behaviors of crumpled graphene sheets. Carbon 2021, 174, 148–157. [Google Scholar] [CrossRef]
  19. Liu, S.; Wang, A.; Li, Q.; Wu, J.; Chiou, K.; Huang, J.; Luo, J. Crumpled Graphene Balls Stabilized Dendrite-free Lithium Metal Anodes. Joule 2018, 2, 184–193. [Google Scholar] [CrossRef] [Green Version]
  20. Luo, J.; Jang, H.D.; Sun, T.; Xiao, L.; He, Z.; Katsoulidis, A.P.; Kanatzidis, M.G.; Gibson, J.M.; Huang, J. Compression and Aggregation-Resistant Particles of Crumpled Soft Sheets. ACS Nano 2011, 5, 8943–8949. [Google Scholar] [CrossRef]
  21. Baimova, J.A.; Liu, B.; Dmitriev, S.V.; Zhou, K. Mechanical properties of crumpled graphene under hydrostatic and uniaxial compression. J. Phys. D Appl. Phys. 2015, 48, 095302. [Google Scholar] [CrossRef]
  22. Khan, M.B.; Wang, C.; Wang, S.; Fang, D.; Chen, S. The mechanical property and microscopic deformation mechanism of nanoparticle-contained graphene foam materials under uniaxial compression. Nanotechnology 2020, 32, 115701. [Google Scholar] [CrossRef] [PubMed]
  23. Cui, X.; Ren, P.; Deng, D.; Deng, J.; Bao, X. Single layer graphene encapsulating non-precious metals as high-performance electrocatalysts for water oxidation. Energy Environ. Sci. 2016, 9, 123–129. [Google Scholar] [CrossRef]
  24. Wang, X.X.; Tan, Z.H.; Zeng, M.; Wang, J.N. Carbon nanocages: A new support material for Pt catalyst with remarkably high durability. Sci. Rep. 2014, 4. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Deng, J.; Ren, P.; Deng, D.; Bao, X. Enhanced Electron Penetration through an Ultrathin Graphene Layer for Highly Efficient Catalysis of the Hydrogen Evolution Reaction. Angew. Chem. Int. Ed. 2015, 54, 2100–2104. [Google Scholar] [CrossRef] [PubMed]
  26. Dahal, A.; Batzill, M. Graphene–nickel interfaces: A review. Nanoscale 2014, 6, 2548. [Google Scholar] [CrossRef] [PubMed]
  27. Sutter, P.; Sadowski, J.T.; Sutter, E. Graphene on Pt(111): Growth and substrate interaction. Phys. Rev. B 2009, 80, 245411. [Google Scholar] [CrossRef] [Green Version]
  28. Gao, L.; Guest, J.R.; Guisinger, N.P. Epitaxial Graphene on Cu(111). Nano Lett. 2010, 10, 3512–3516. [Google Scholar] [CrossRef]
  29. Yoon, S.M.; Choi, W.M.; Baik, H.; Shin, H.J.; Song, I.; Kwon, M.S.; Bae, J.J.; Kim, H.; Lee, Y.H.; Choi, J.Y. Synthesis of Multilayer Graphene Balls by Carbon Segregation from Nickel Nanoparticles. ACS Nano 2012, 6, 6803–6811. [Google Scholar] [CrossRef]
  30. Ji, Z.; Wang, Y.; Yu, Q.; Shen, X.; Li, N.; Ma, H.; Yang, J.; Wang, J. One-step thermal synthesis of nickel nanoparticles modified graphene sheets for enzymeless glucose detection. J. Colloid Interface Sci. 2017, 506, 678–684. [Google Scholar] [CrossRef]
  31. Ai, L.; Tian, T.; Jiang, J. Ultrathin Graphene Layers Encapsulating Nickel Nanoparticles Derived Metal–Organic Frameworks for Highly Efficient Electrocatalytic Hydrogen and Oxygen Evolution Reactions. ACS Sustain. Chem. Eng. 2017, 5, 4771–4777. [Google Scholar] [CrossRef]
  32. Batzill, M. The surface science of graphene: Metal interfaces, CVD synthesis, nanoribbons, chemical modifications, and defects. Surf. Sci. Rep. 2012, 67, 83–115. [Google Scholar] [CrossRef]
  33. Omukai, A.; Yoshimura, A.; Watanabe, F.; Tachibana, M. Graphene films synthesized by chemical vapor deposition with ethanol. Trans. Mater. Res. Soc. Jpn. 2011, 36, 359–362. [Google Scholar] [CrossRef] [Green Version]
  34. Li, X.; Cai, W.; Colombo, L.; Ruoff, R.S. Evolution of Graphene Growth on Ni and Cu by Carbon Isotope Labeling. Nano Lett. 2009, 9, 4268–4272. [Google Scholar] [CrossRef] [Green Version]
  35. Chae, S.J.; Güneş, F.; Kim, K.K.; Kim, E.S.; Han, G.H.; Kim, S.M.; Shin, H.J.; Yoon, S.M.; Choi, J.Y.; Park, M.H.; et al. Synthesis of Large-Area Graphene Layers on Poly-Nickel Substrate by Chemical Vapor Deposition: Wrinkle Formation. Adv. Mater. 2009, 21, 2328–2333. [Google Scholar] [CrossRef]
  36. Reina, A.; Thiele, S.; Jia, X.; Bhaviripudi, S.; Dresselhaus, M.S.; Schaefer, J.A.; Kong, J. Growth of large-area single- and Bi-layer graphene by controlled carbon precipitation on polycrystalline Ni surfaces. Nano Res. 2009, 2, 509–516. [Google Scholar] [CrossRef] [Green Version]
  37. Chang, S.W.; Nair, A.K.; Buehler, M.J. Nanoindentation study of size effects in nickel–graphene nanocomposites. Philos. Mag. Lett. 2013, 93, 196–203. [Google Scholar] [CrossRef]
  38. Kim, Y.; Lee, J.; Yeom, M.S.; Shin, J.W.; Kim, H.; Cui, Y.; Kysar, J.W.; Hone, J.; Jung, Y.; Jeon, S.; et al. Strengthening effect of single-atomic-layer graphene in metal–graphene nanolayered composites. Nat. Commun. 2013, 4. [Google Scholar] [CrossRef]
  39. Kuang, D.; Xu, L.; Liu, L.; Hu, W.; Wu, Y. Graphene–nickel composites. Appl. Surf. Sci. 2013, 273, 484–490. [Google Scholar] [CrossRef]
  40. Kim, D.J.; Truong, Q.T.; Kim, J.I.; Suh, Y.; Moon, J.; Lee, S.E.; Hong, B.H.; Woo, Y.S. Ultrahigh-strength multi-layer graphene-coated Ni film with interface-induced hardening. Carbon 2021, 178, 497–505. [Google Scholar] [CrossRef]
  41. Safina, L.R.; Baimova, J.A.; Mulyukov, R.R. Nickel nanoparticles inside carbon nanostructures: Atomistic simulation. Mech. Adv. Mater. Mod. Process. 2019, 5. [Google Scholar] [CrossRef]
  42. Safina, L.R.; Krylova, K.A. Effect of particle size on the formation of the composite structure in Ni-graphene system: Atomistic simulation. J. Phys. Conf. Ser. 2020, 1435, 012067. [Google Scholar] [CrossRef]
  43. Safina, L.L.; Baimova, J.A. Molecular dynamics simulation of fabrication of Ni-graphene composite: Temperature effect. Micro Nano Lett. 2020, 15, 176–180. [Google Scholar] [CrossRef]
  44. Safina, L.; Baimova, J.; Krylova, K.; Murzaev, R.; Mulyukov, R. Simulation of metal-graphene composites by molecular dynamics: A review. Lett. Mater. 2020, 10, 351–360. [Google Scholar] [CrossRef]
  45. Wu, C.D.; Fang, T.H.; Lo, J.Y. Effects of pressure, temperature, and geometric structure of pillared graphene on hydrogen storage capacity. Int. J. Hydrog. Energy 2012, 37, 14211–14216. [Google Scholar] [CrossRef]
  46. Krylova, K.; Baimova, J.; Mulyukov, R. Effect of deformation on dehydrogenation mechanisms of crumpled graphene: Molecular dynamics simulation. Lett. Mater. 2019, 9, 81–85. [Google Scholar] [CrossRef]
  47. Krylova, K.A.; Baimova, J.A.; Lobzenko, I.P.; Rudskoy, A.I. Crumpled graphene as a hydrogen storage media: Atomistic simulation. Phys. B Condens. Matter 2020, 583, 412020. [Google Scholar] [CrossRef]
  48. Dimitrakakis, G.K.; Tylianakis, E.; Froudakis, G.E. Pillared Graphene: A New 3-D Network Nanostructure for Enhanced Hydrogen Storage. Nano Lett. 2008, 8, 3166–3170. [Google Scholar] [CrossRef] [PubMed]
  49. Ranjbar, A.; Khazaei, M.; Venkataramanan, N.S.; Lee, H.; Kawazoe, Y. Chemical engineering of adamantane by lithium functionalization: A first-principles density functional theory study. Phys. Rev. B 2011, 83. [Google Scholar] [CrossRef] [Green Version]
  50. Yang, R.T. Hydrogen storage by alkali-doped carbon nanotubes–revisited. Carbon 2000, 38, 623–626. [Google Scholar] [CrossRef]
  51. Khantha, M.; Cordero, N.A.; Molina, L.M.; Alonso, J.A.; Girifalco, L.A. Interaction of lithium with graphene: Anab initiostudy. Phys. Rev. B 2004, 70. [Google Scholar] [CrossRef]
  52. Medeiros, P.V.C.; de Brito Mota, F.; Mascarenhas, A.J.S.; de Castilho, C.M.C. Adsorption of monovalent metal atoms on graphene: A theoretical approach. Nanotechnology 2010, 21, 115701. [Google Scholar] [CrossRef]
  53. Yoshida, A.; Okuyama, T.; Terada, T.; Naito, S. Reversible hydrogen storage/release phenomena on lithium fulleride (LinC60) and their mechanistic investigation by solid-state NMR spectroscopy. J. Mater. Chem. 2011, 21, 9480. [Google Scholar] [CrossRef]
  54. Waqar, Z. Hydrogen accumulation in graphite and etching of graphite on hydrogen desorption. J. Mater. Chem. 2007, 42, 1169–1176. [Google Scholar] [CrossRef]
  55. oś, S.; Duclaux, L.; Letellier, M.; Azaïs, P. Study of Adsorption Properties on Lithium Doped Activated Carbon Materials. Acta Phys. Pol. A 2005, 108, 371–377. [Google Scholar] [CrossRef]
  56. Stuart, S.J.; Tutein, A.B.; Harrison, J.A. A reactive potential for hydrocarbons with intermolecular interactions. J. Chem. Phys. 2000, 112, 6472–6486. [Google Scholar] [CrossRef] [Green Version]
  57. Humphrey, W.; Dalke, A.; Schulten, K. VMD: Visual molecular dynamics. J. Mol. Gragh. 1996, 14, 33–38. [Google Scholar] [CrossRef]
  58. Girifalco, L.A.; Weizer, V.G. Application of the Morse Potential Function to Cubic Metals. Phys. Rev. 1959, 114, 687–690. [Google Scholar] [CrossRef]
  59. Katin, K.P.; Prudkovskiy, V.S.; Maslov, M.M. Molecular dynamics simulation of nickel-coated graphene bending. Micro Nano Lett. 2018, 13, 160–164. [Google Scholar] [CrossRef]
  60. Galashev, A.Y.; Katin, K.P.; Maslov, M.M. Morse parameters for the interaction of metals with graphene and silicene. Phys. Lett. A 2019, 383, 252–258. [Google Scholar] [CrossRef]
  61. Petucci, J.; LeBlond, C.; Karimi, M.; Vidali, G. Diffusion, adsorption, and desorption of molecular hydrogen on graphene and in graphite. J. Chem. Phys. 2013, 139, 044706. [Google Scholar] [CrossRef] [Green Version]
  62. Ganz, E.; Ganz, A.B.; Yang, L.M.; Dornfeld, M. The initial stages of melting of graphene between 4000 K and 6000 K. Phys. Chem. Chem. Phys. 2017, 19, 3756–3762. [Google Scholar] [CrossRef] [PubMed]
  63. Los, J.H.; Zakharchenko, K.V.; Katsnelson, M.I.; Fasolino, A. Melting temperature of graphene. Phys. Rev. B 2015, 91, 045415. [Google Scholar] [CrossRef] [Green Version]
  64. Ye, Y.; Ahn, C.C.; Witham, C.; Fultz, B.; Liu, J.; Rinzler, A.G.; Colbert, D.; Smith, K.A.; Smalley, R.E. Hydrogen adsorption and cohesive energy of single-walled carbon nanotubes. Appl. Phys. Lett. 1999, 74, 2307–2309. [Google Scholar] [CrossRef] [Green Version]
  65. Pinkerton, F.E.; Wicke, B.G.; Olk, C.H.; Tibbetts, G.G.; Meisner, G.P.; Meyer, M.S.; Herbst, J.F. Thermogravimetric Measurement of Hydrogen Absorption in Alkali-Modified Carbon Materials. J. Phys. Chem. B 2000, 104, 9460–9467. [Google Scholar] [CrossRef]
  66. Schimmel, H.G.; Kearley, G.J.; Nijkamp, M.G.; Visser, C.T.; de Jong, K.P.; Mulder, F.M. Hydrogen Adsorption in Carbon Nanostructures: Comparison of Nanotubes, Fibers, and Coals. Chem. A Eur. J. 2003, 9, 4764–4770. [Google Scholar] [CrossRef] [Green Version]
Figure 1. (a) Initial structure: graphene flake with the diameter D filled with nanocluster with the diameter d; (b) atomic and molecular hydrogen inside the graphene flake.
Figure 1. (a) Initial structure: graphene flake with the diameter D filled with nanocluster with the diameter d; (b) atomic and molecular hydrogen inside the graphene flake.
Materials 14 02098 g001
Figure 2. Potential energy as the function of exposure time at 300 K (a) and 1000 K (b) for five types of nanoclusters.
Figure 2. Potential energy as the function of exposure time at 300 K (a) and 1000 K (b) for five types of nanoclusters.
Materials 14 02098 g002
Figure 3. Potential energy as the function of exposure time for Ni 21 (a), Ni 47 (b) and Ni 66 (c). Corresponding snapshots of the structure of crumpled graphene flake filled with nanocluster during exposure at 300 K (bottom line of snapshots) and 1000 K (upper line of snapshots). Carbon atoms are shown by violet and nickel atoms are shown by green.
Figure 3. Potential energy as the function of exposure time for Ni 21 (a), Ni 47 (b) and Ni 66 (c). Corresponding snapshots of the structure of crumpled graphene flake filled with nanocluster during exposure at 300 K (bottom line of snapshots) and 1000 K (upper line of snapshots). Carbon atoms are shown by violet and nickel atoms are shown by green.
Materials 14 02098 g003
Figure 4. Snapshots of the structure of crumpled graphene filled with (ac) Ni 21 and (a’c’) Ni 78 nanocluster. Initial structure was obtained by annealing at 300 K, while composite was obtained by hydrostatic compression at 1000 K. Ni atoms are shown by purple and C atoms are shown by green color.
Figure 4. Snapshots of the structure of crumpled graphene filled with (ac) Ni 21 and (a’c’) Ni 78 nanocluster. Initial structure was obtained by annealing at 300 K, while composite was obtained by hydrostatic compression at 1000 K. Ni atoms are shown by purple and C atoms are shown by green color.
Materials 14 02098 g004
Figure 5. Potential energy as the function of exposure time at 77 K (a) and 300 K (b) for five types of hydrogen nanoclusters.
Figure 5. Potential energy as the function of exposure time at 77 K (a) and 300 K (b) for five types of hydrogen nanoclusters.
Materials 14 02098 g005
Figure 6. Potential energy as the function of exposure time for three types of hydrogen nanoclusters: (a) 21 atoms; (b) 47 atoms; and (c) 78 atoms. Corresponding snapshots of graphene flake filled with nanocluster during exposure at 77 K (bottom line of snapshots) and 300 K (upper line of snapshots). Carbon atoms are shown by violet and hydrogen atoms are shown by blue.
Figure 6. Potential energy as the function of exposure time for three types of hydrogen nanoclusters: (a) 21 atoms; (b) 47 atoms; and (c) 78 atoms. Corresponding snapshots of graphene flake filled with nanocluster during exposure at 77 K (bottom line of snapshots) and 300 K (upper line of snapshots). Carbon atoms are shown by violet and hydrogen atoms are shown by blue.
Materials 14 02098 g006
Figure 7. Snapshots of graphene flake with 78 hydrogen atoms during exposure at 77 K for 20 ps. Colors are as in Figure 6.
Figure 7. Snapshots of graphene flake with 78 hydrogen atoms during exposure at 77 K for 20 ps. Colors are as in Figure 6.
Materials 14 02098 g007
Figure 8. Snapshots of crumpled graphene filled with hydrogen atoms before (a) and after (b) hydrostatic compression. Colors are as in Figure 6.
Figure 8. Snapshots of crumpled graphene filled with hydrogen atoms before (a) and after (b) hydrostatic compression. Colors are as in Figure 6.
Materials 14 02098 g008
Table 1. The size of the structural elements.
Table 1. The size of the structural elements.
D, ÅNd, Åa, Å
216.24.25
388.13.30
14.7479.02.85
669.42.65
7810.42.15
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Safina, L.R.; Krylova, K.A.; Murzaev, R.T.; Baimova, J.A.; Mulyukov, R.R. Crumpled Graphene-Storage Media for Hydrogen and Metal Nanoclusters. Materials 2021, 14, 2098. https://doi.org/10.3390/ma14092098

AMA Style

Safina LR, Krylova KA, Murzaev RT, Baimova JA, Mulyukov RR. Crumpled Graphene-Storage Media for Hydrogen and Metal Nanoclusters. Materials. 2021; 14(9):2098. https://doi.org/10.3390/ma14092098

Chicago/Turabian Style

Safina, Liliya R., Karina A. Krylova, Ramil T. Murzaev, Julia A. Baimova, and Radik R. Mulyukov. 2021. "Crumpled Graphene-Storage Media for Hydrogen and Metal Nanoclusters" Materials 14, no. 9: 2098. https://doi.org/10.3390/ma14092098

APA Style

Safina, L. R., Krylova, K. A., Murzaev, R. T., Baimova, J. A., & Mulyukov, R. R. (2021). Crumpled Graphene-Storage Media for Hydrogen and Metal Nanoclusters. Materials, 14(9), 2098. https://doi.org/10.3390/ma14092098

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop