Influence of Size Effect on the Electronic and Elastic Properties of Graphane Nanoflakes: Quantum Chemical and Empirical Investigations
By application of empirical method it is found that graphene nanoflake (graphane) saturated by hydrogen is not elastic material. In this case, the modulus of the elastic compression of graphane depends on its size, allowing us to identify the linear parameters of graphane with maximum Young’s modulus for this material. The electronic structure of graphane nanoflakes was calculated by means of the semiempirical tight-binding method. It is found that graphane nanoflakes can be characterized as dielectric. The energy gap of these particles decreases with increasing of the length tending to a certain value. At the same time, the ionization potential of graphane also decreases. A comparative analysis of the calculated values with the same parameters of single-walled nanotubes is performed.
Graphane is a principal new and prospect material in nanoelectronics. Graphane was synthesized for the first time by the hydrogenation of graphene—a monolayer of graphite . The crystal structure of graphane is two-dimensional and hexagonal. The hydrogen atoms are attached at opposite sides of the carbon atoms plane. It is necessary to distinguish the graphane nanoflakes and nanoribbons. The sizes of nanoflakes in the different directions differ no more than three times and do not exceed 100 nanometers.
The discovery of graphane has created prerequisites for the investigations of its properties and for the search of the probable application. For example, graphane has unique optical properties . It is found that the dielectric functions of graphane nanoribbons spectra do not depend on the ribbon edge shapes and widths. Also, in the result of the graphane optical properties calculations, the presence of moderate anisotropy with respect to the type of light polarization is established. In , authors investigate thermal conduction in graphane nanoribbons by the nonequilibrium Green function method. It is demonstrated that thermal conduction can be effectively tuned by controlling the edge shape, width, and hydrogen vacancy concentration of graphane nanoribbons. Specifically, ballistic thermal conductance of graphane nanoribbons generally grows with the ribbon width. Big attention of researchers is devoted to investigation of the magnetic properties of graphane. In , it is shown that in the whole temperature range up to room temperature the hydrogenated graphene demonstrates a weak ferromagnetism. The origin of the magnetism is also determined to arise from the hydrogenated graphene itself.
Due to its unique properties, graphane has many applications. In particular, this material can find the application in the hydrogen energetics. It is established that the heating of graphane leads to release of the atomic hydrogen . Therefore, graphane can be considered as one of the effective methods of hydrogen storage. Another significant application of graphane is its application in nanoelectronics as a base for the printed circuits with conducting and nonconducting sites . One more possible field of graphane application is biosensing. In experimental work , the properties of graphane were considered for application in electrochemical oxidation of biomarkers. That can lead to finding of harmful ways for detection of biomarkers.
In present work , authors have studied theoretically the possibility of a superconducting state in two-dimensional hole-doped graphane. According to results of this work, the hole-doped graphane is a candidate for being a high temperature superconductor with a critical temperature as high as those of copper oxides.
One more perspective applied field of graphane is its use as thermoelectric materials for thermoelectric device. It is predicted that disordered armchair graphane nanoribbons with low thermal conductivity are promising candidates for constructing thermoelectric materials .
Although the properties of graphane nanoribbons have been studied extensively, the influence of size effect on the properties of graphane remained unexplored. At the same time, the variation of the geometrical dimensions of the objects could be an effective way to control the properties of graphane. The aim of this work is the theoretical investigation of the electronic and elastic properties of graphane nanoflakes by means of the quantum-chemical and empirical methods.
2. Materials and Methods
2.1. Empirical Method
To study the elastic properties and deformations of graphane nanoflakes we apply the empirical method based on the bond-order potential developed by Brenner . The total energy is described by the sum of the binding energy , the torsional energy , and the van der Waals energy :
Each pair of covalently bonded atoms interacts via potential-energy:
This is the binding energy. Here is the repulsive pair term, is the attractive pair term, and is the distance between the chemically bonded atoms with numbers and . The function is a many-body term. This term was introduced to describe the specificity of the - interaction. So, the value of the binding energy depends on the position and chemical identity of atoms.
To describe nonbonded interaction between the atoms we have attended the torsional interaction energy to the total system energy. The torsional energy is given by the following formula:
The torsional potential is given as a function of dihedral angle . The torsion angle is defined in the usual way as the angle between the plane defined by the vectors and and that defined by and . Here atoms and are not chemically bonded.
Van der Waals energy defines the interaction between nonbonded atoms:
Van der Waals interaction energy is described by the Morse potential aswhere is the average bond energy, is the repulsion nucleus energy, and coefficients nm−1 and nm−1. Calculated interatomic distance in graphite by described empirical method is equal to 0.1421 nm.
2.2. Tight-Binding Method
Within this method the total energy is calculated by the following formula:where is the bond structure energy calculated as the sum of energies of the single-particle occupied states. is the interelectron and internuclear repulsion term.
The bond structure energy is determined by the following formula:
This expression is the sum of energies of the molecular orbitals obtained by diagonalizing the Hamiltonian, is the number of the occupied orbitals, and is the energy of the single-particle orbitals.
The phenomenon energy can be expressed as a sum of two-body potentials:where , are the number of the interaction atoms; , are the Cartesian coordinates.
The overlap matrix elements are calculated by the formulas which take into account four types of interaction, , , , and , and the pair repulsive potential:where and are the orbital moments of wave function and presents the bond type ( or ). The atomic terms , are the atomic orbital energies of carbon which are located on the main diagonal of the Hamiltonian. The values of the , atomic terms, the parameters (), and the equilibrium overlap integrals (, , , and ) are given in Table 1. The coefficients , , and are the power exponents, the coefficients and are determined by the overall shape and steepness of the function, and the parameter is a cutoff distance for the hopping matrix elements and repulsive interactions.
The geometric and energy parameters for the carbon nanotubes and fullerenes obtained by the tight binding method are in good agreement with experimental data [11, 12] and ab initio data . The distances between carbon atoms and ionization potential are equal to 0.14 nm and 7.61 eV, respectively.
3. Results and Discussion
3.1. Elastic Properties of Graphane Nanoflakes
Elastic properties of graphane nanoflakes are investigated by means of the empirical method based on Brenner potential. We investigated graphane nanoflakes with lengths of 1–5 nm and width of 2.4 nm. As a result of the optimization it is established that C-C bond length is equal to 0.155 nm, and C-H bond length is equal to 0.112 nm. Obtained results are in good agreement with the results of DFT PBE calculations (0.152 and 0.111 nm for C-C and C-H, resp.) .
In our study, graphane nanoflakes are subjected to axial compression load and stretching load. The edges of the structures were rigidly fixed. Reoptimization of the structures is carried out on each stage of numerical experiment of compression. We considered cases of compression for graphane nanoflakes on 1, 2, and 3%. We observed a change of the graphane structure from planar to a wave-like at the 3% of compression (Figure 1). Thus, the behavior of graphane nanoflake at axial compression is similar to behavior of graphene under axial compression deformation . The change of the deformation energy with axial compression is presented in Figure 2. From the figure it can be seen that the deformation energy increases quadratically with compression. However, a sharp energy decrease occurs at the 2.5% compression. One can assume that this fact indicates the phase transition from elastic deformation to plastic.
We calculated pseudo-Young’s modulus for the investigated structure of graphane nanoflakes. Pseudo-Young’s modulus is calculated on the base of the structure edge length. The deforming force acts on this edge. Pseudo-Young’s modulus was chosen as a parameter of the graphane elasticity, since the graphene is a nanoplate. Therefore, such estimation of the elasticity is more justified than classical Young’s modulus. The algorithm for calculation of pseudo-Young’s modulus is as follows.(1)At the initial stage the optimization of the structure by means of the total energy minimization was performed.(2)Then the optimized structure was compressed on 1%. After that the edges of the structure were rigidly fixed, and a structure was optimized again.(3)Pseudo-Young’s modulus was calculated according to the formula where is the length of graphane nanoflake; is the edge length of graphane nanoflake; is the force, necessary for tube stretching. This force is determined according to the formula where is energy of the elastic extension and is the length of the structure deformation.
Pseudo-Young’s modulus was calculated according to the presented algorithm for graphane nanoflakes with length of 2.4 nm and width in the range of 1–5 nm. The plot of the change for the elasticity parameter is shown in Figure 3. From the plot it can be seen that pseudo-Young’s modulus does not change monotonically. It sharply increases starting with 3 nm length. This fact indicates an increase of elasticity. We can conclude that the elastic properties of the graphane nanoflakes are markedly inferior to the nanotubes on the elasticity. Pseudo-Young’s modulus of the graphane nanoflakes reduced by ~72% in comparison to graphene (0.057–0.08 TPa·nm of the graphane nanoflakes and 0.22–0.25 TPa·nm  of zigzag and armchair nanotubes). Graphane nanoflakes show behavior that is typical for the case of plastic deformation.
3.2. Electronic Properties of Graphane Nanoflakes
Electronic properties of graphane nanoflakes were studied by a semiempirical tight-binding method. At first the atomic structure for graphane nanoflake of different sizes was optimized (by means of the total energy minimization by coordinates of the atoms). Then, the ionization potential and the energy gap of the spectrum were calculated by the electronic spectrum (the energy interval between the last occupied level and the first unoccupied level).
We carried out a set of the calculations for the ionization potential and the energy gap at the consecutive increase of the atomic skeleton for graphane particle. A plot of the ionization potential of graphane nanoflake on its length is presented in Figure 4. The linear decrease of the ionization potential with the increasing of the structure size is observed. It can be assumed that the ionization potential will tend to a work function of graphane with the further increase of the length.
The energy gap of the electron spectrum is also reduced with the increasing of the atomic skeleton of graphane particles. The plot of the energy gap of graphane nanoflakes on its length is presented in Figure 5. As can be seen from the plot, the energy gap takes on a value of about 7 eV. Therefore, we can conclude that graphane nanoflakes are good dielectrics.
The analysis of publications devoted to investigation of electron structures of hydrogenated graphene showed that electron properties of just periodic graphane nanostructures were studied. In our work appearance of size effect in electron and energetic features of graphane was proved for the first time. For periodic graphane structure of armchair-type it is established that the ionization potential is equal to 4.73 eV and the energy gap of spectrum is equal to 6.05 eV . Analyzing the dependencies of the ionization potential and energy gap on graphane length presented in our paper, it can be concluded that the electron and energetic characteristics of graphane with growing of atomic skeleton will tend to values presented in .
The current paper presented results of the theoretical investigation of size effect influence on the elastic and electronic properties of graphane nanoflakes. Theoretical study of the deformation for graphane nanoflakes revealed that graphane nanoflakes are not elastic material. This is evidenced by small pseudo-Young’s modulus as compared with value of nanotubes and the appearance of the plastic deformation of the structure at already 3% compression.
Also it is found that graphane nanoflakes are dielectric with sufficiently large energy gap of the electronic spectrum. Therefore, one can obtain contact tracks with electronic conductivity of desired configuration by removing of the hydrogen atoms from graphane surface. The remainder of the surface will remain dielectric.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
The authors gratefully acknowledge funding of this work by Ministry of Education in the framework of the project of the state task (no. 3.1155.2014/K) and by Presidential scholarship 2013–2015 (Project no. SP-2302.2013.1).
T. Yatsuhashi and N. Nakashima, “Direct ionization desorption of fullerene by intense femtosecond laser fields,” The Review of Laser Engineering. Supplements, vol. 36, pp. 1000–1003, 2008.View at: Google Scholar
L. Stobinski, J. Peszke, and H.-M. Lin, “Computational studies of SWCNTs capped by hemispheres of C60 fullerene, based on semi-empirical methods,” Reviews on Advanced Materials Science, vol. 5, no. 4, pp. 363–370, 2003.View at: Google Scholar
O. Glukhova and M. Slepchenkov, “Influence of the curvature of deformed graphene nanoribbons on their electronic and adsorptive properties: theoretical investigation based on the analysis of the local stress field for an atomic grid,” Nanoscale, vol. 4, no. 11, pp. 3335–3344, 2012.View at: Publisher Site | Google Scholar
O. E. Glukhova and O. A. Terent'ev, “Theoretical study of the dependences of the young's and torsion moduli of thin single-layer carbon zigzag and armchair nanotubes on the geometric parameters,” Physics of the Solid State, vol. 48, no. 7, pp. 1411–1417, 2006.View at: Publisher Site | Google Scholar