Abstract
Irradiation dynamics of a single graphene sheet bombarded by hydrogen atoms is studied in the incident energy range of 0.1 to 200 eV. Results for reflection, transmission, and adsorption probabilities, as well as effects of a single adsorbed atom to the electronic properties of graphene, are obtained by the quantumclassical Monte Carlo molecular dynamics within a selfconsistentchargedensity functional tight binding formalism We compare these results with those, distinctly different, obtained by the classical molecular dynamics.
PACS: 61.80.Az, 61.48.Gh, 61.80.Jh, 34.50.Dy.
Keywords:
Graphene; DFTB; Hydrogen detection; HOMOLUMO gap; Molecular dynamicsBackground
The sp^{2 }hybridized carbon allotrope, graphene, has recently shown particular promise in applications such as nanoscale electronics, hydrogen storage [1], and nanosensors. This is due to the material's outstanding thermal and electronic properties. The sensitivity of the electronic properties of a single graphene sheet to small defects in its 2D crystal structure and chemical composition indicates a possibility of its application as a fewparticle detector [24]. Graphenebased electronics in space vehicles might also be sensitive to the damages caused by cosmic radiation containing a wide spectrum of particles, a significant component of which would be light atoms from the solar wind. The significance of studies of graphene bombarded by hydrogenic atoms in understanding the damages of the CFC carbon tiles in the divertor of a fusion reactor (ITER) to the plasma irradiation has also been stressed recently [5,6]. These defects include lattice defects, with possible creation of vacancies, as well as chemical changes induced by the hydrogen sticking to the lattice [7,8]. The resultant changes in the electronic conductance due to changes in the electronic structure have also been studied [3,9]. For example, work by Deretzis et al. [2] has shown that even single vacancy deformations in graphene nanoribbons can have measurable effects on the material's conduction properties. These applications all motivate our study of energetic particle impact with graphene.
In this paper, we study the perpendicular impact of hydrogen on a single graphene sheet over more than three decades of impact energies (0.1 to 200 eV) using methods of quantumclassical Monte Carlo molecular dynamics. Our approach is described in detail in the second section entitled 'Methods'. The irradiated target was an infinite graphene sheet obtained by applying 2D periodic boundary conditions to a graphene cell of size 29.12 × 28.53 Å (336 C atoms). The graphene was prepared at a temperature of 300 K by a NoseHoover thermostat and left free during each collision event, which lasted 200 to 500 fs, depending on the impact energy. The irradiation was performed by more than 1,000 independent trajectories for each impact energy, with randomly chosen position of emission of an atom above the surface of the graphene cell. In this method, the total electronic energy of the system is solved quantummechanically at the beginning of each time step (on the order of a femtosecond), maintaining fixed atom positions; after incorporating the nucleusnucleus interaction into the total electronic energy, forces on each atom are updated, and the atoms are moved classically within the time step. The electronic structure is solved here by the selfconsistentchargedensity functional tight binding (SCCDFTB) method [1012]. To allow for the highenergy impact, we fit the original SCCDFTB parameters [13] at close distances (< 0.2 Å) to the binary ZieglerBiersackLittmark (ZBL) [14] repulsive potentials.
Results for reflection and transmission probabilities, angular distributions, and adsorption probabilities at low energies (0.1 to 1 eV) are shown and analyzed in the first part of the 'Results and discussion' section, entitled 'Irradiation dynamics and effects on electronic structure'. Additionally, changes in the molecular orbital levels close to the Fermi energy, which influence the nonequilibrium ballistic electron transport properties (i.e., the electric conductance) of the system, are calculated and characterized by the changes, ∆E_{lh}, in the difference, E_{lh}, of the (discrete) lowest unoccupied molecular orbital and highest occupied molecular orbital energies in response to the hydrogen adsorption. These changes are indicative of possible changes in the graphene sheet conductance. They are, surprisingly, on the order of 1 eV and depend on the vibrational energy of the adsorbed hydrogen. Adsorption occurs only for the lowenergy impacts (< 1 eV). This confirms some predictions in literature on the extreme sensitivity of the highest occupied molecular orbital (HOMO)lowest unoccupied molecular orbital (LUMO) gap and transport properties of graphene and SWCNT to the adsorption of hydrogen and other atoms and molecules [1519].
In the second part of the 'Results and discussion' section entitled 'Comparison with classical molecular dynamics', we perform classical molecular dynamics (CMD) calculations with two stateoftheart bond order hydrocarbon potentials, reactive empirical bondorder (REBO) [20] and adaptive intermolecular reactive empirical bond order (AIREBO) [21]. We use the corrected set of the classical potentials [22] to allow high impact energies and compare the classical MD probabilities with our quantumclassical results. Although CMD with these potentials is significantly faster than SCCDFTB, allowing for longer timescales, larger systems, and greater energy ranges to be studied, it turns out that the classical potentials are of limited applicability for the studied system and dynamics. We hope that this data motivates improvements to these potentials since their speed is very attractive for radiation damagetype problems. Our conclusions are given in the final section.
Methods
To simulate effects of irradiation on graphene, one can apply direct molecular dynamics methods in which electronic structure is treated explicitly using quantum mechanics, while the motion of the nuclei is described by the means of the classical dynamics. This allows one to accurately describe bond breaking and formation as well as the interatomic potentials. Such an approach is, however, computationally very expensive, which greatly limits the system sizes, timescales, and choice of quantum mechanicsbased methods. To mimic the dynamics observed by experiment, we apply a Monte Carlo approach to the trajectories, i.e., using a large number of trajectories, randomly varying 'impact parameters' to obtain acceptable statistics of the collision events. Even using this approach, we must use a less expensive and more approximate quantummechanical approach. Here, we use the SCCDFTB method, an approximate density functional theory (DFT) method in which only valence electron interactions are considered. Although a full DFT treatment would be ideal, this is currently too expensive computationally, even for a handful of trajectories. In SCCDFTB, the total electronic densities and energies are expressed by solution of the Schrodinger equation in the KohnSham form, using predetermined Hamiltonian and overlap integrals as well as repulsive splines fit to reference systems (socalled SlaterKoster parameters). The tight binding methods applied to the large (solidstate) systems have a long history. Here, we use a selfconsistent charge version developed by Bremen Group (Bremen, Germany) [1012]. SCC is a secondorder correction term in the DFTB total energy involving interactions between localized fluctuations of the electron density; it uses an iterative procedure to converge on the new electron density at each time step. In this SCCDFTB method, spin polarization is neglected. We employed a FermiDirac smearing with electronic temperature T_{el }= 1,000 K, which has a similar effect to averaging over many electronic states near the Fermi level.
To safely allow for highenergy bombardment simulations (in our case 200 eV), we use a refitted version of the original DFTB PBC03 [13] parameters obtained by fitting to the ZBL [14] repulsive interactions at short distances (< 0.2 Å). The PBC03 parameters used here have already shown good results for the hydrogenation of periodic graphene [23] at thermal energies. We show in Figure 1 the potential energy curves of a hydrogen atom interacting with a coronene molecule obtained by the SCCDFTB using PBC03/ZBL parameters and by DFT using a local density approximation functional [24]. At distances closer than 1.5 Å, agreement between DFT [25] and DFTB potentials is quite good. Between 1.5 and 4 Å, DFTB potentials overestimate bond strength, and wells are about 0.5 Å closer to the surface than their DFTB counterparts. Also notable is the lack of convergence of the three potentials until they approach 0 eV. Although SCCDFTB underestimates bonding at the bond center and lattice point positions, these are qualitatively similar to DFT potentials [25]. The problem of thermal atom adsorption gave rise to many experimental and theoretical papers [7,1519] and references therein. The previously reported SCCDFTB studies [26] of collisioninduced reactions in carbon materials within the same energy range considered here were in excellent agreement with experimental findings. Additional comparisons between DFT and SCCDFTB are contained in the studies of Zheng et al. and Elstner [2729].
Figure 2 compares the SCCDFTB potential energy of the hydrogengraphene and hydrogencoronene interactions as a function of zposition above the graphene/coronene plane. The coronene potentials show bonding that is roughly 1 eV weaker and a potential barrier at the hexagon center that is 1 eV higher, reflecting the changes in electronic structure between hydrogenterminated and periodic sp^{2 }carbon. Despite these differences, the forms of the Hgraphene and Hcoronene interactions are very similar. Thus, the agreement of SCCDFTB with DFT calculations of the coronene molecule in Figure 1 indicates that the PBC03/ZBL SCCDFTB parameters are as acceptable for use with graphene as the DFT approach.
Figure 2. Potential energy of the Hgraphene and Hcoronene interactions at analogous points in the lattice. As calculated by the SCCDFTB with the PBC03/ZBL parameters.
Notably, there are two bonding regions in the Hgraphene potential. For incidence directly upon a lattice carbon, the potential minimum occurs at approximately 1.1 Å, while incidence upon a CC bond center shows a shallower potential with minimum close to 1.0 Å. Indeed, there are many potential wells in the 3D multibody potential that are responsible for capturing impinging hydrogen atoms; these will later be shown to have an effect on the electronic structure of resultant Hgraphene surfaces. There are repulsive barriers at the bond center and hexagon center of heights 17 and 2.5 eV, respectively. Notably, hydrogen encounters no barrier before entering the potential well when incident directly on a lattice carbon.
About ten per decade incident kinetic energies ranging from 0.1 to 200 eV are considered for the impinging hydrogen atom. While cumulative bombardment is not investigated, 1,008 single impact simulations are performed for each incident energy; this is achieved using 1,008 processors, one for each trajectory, on the Kraken Cray XT5 supercomputer (National Institute of Computational Sciences, University of Tennessee, Knoxville, TN, USA). The target graphene surface described in the 'Background' is situated in the z = 0 plane and periodically extended in the xy coordinate plane. To simulate the bombardment in a realworld environment, the sample is thermostated (via NoseHoover scheme) to 300 K before bombardment and left to evolve freely during approximately 0.1 to 1 ps (depending on incident energy) simulation time. The impinging hydrogen atom is released from a random (x, y) position in the z = 10 Å plane, with velocity perpendicular to the graphene sheet.
Results and discussion
Irradiation dynamics and effects on electronic structure
Three outcomes of the bombardment are observed: reflection, transmission, and adsorption of the incident hydrogen atom; no sputtering of any type was observed in our quantumclassical approach. Figure 3 shows the probabilities of these processes as a function of incident Hatom energy. At 20 eV and above, transmission is the dominant process, as expected from the potentials in Figure 2. At the midrange energies of 1 to 10 eV, reflection is primarily observed, with a peak at 2 eV. At 1 eV, H still transfers enough kinetic energy to the target carbon atoms to allow its bonding in the wells of depth approximately 0.5 eV near the lattice points and bond centers. As the incident energy becomes comparable with the depth of this and smaller wells, adsorption becomes the dominant process as expected.
Figure 3. Probabilities of reflection, transmission, and adsorption as a function of incident kinetic energy.
Reflection of the incident hydrogen can occur at all points in the graphene lattice. As can be seen in Figure 2, the threshold for transmission is approximately 2.5 eV at the hexagon center. These atoms are still of insufficient energy to penetrate the barrier at the CC bond position, so those that do not impact near the center of the hexagon are reflected (see Figure 4).
Figure 4. Positions of reflection, transmission, and adsorption events for the quantumclassical calculations. In a representative graphene hexagon, using SCCDFTB. Adsorption (left) shows clustering of hydrogen atoms around the lattice carbons. Reflection (center) is most probable at the perimeter of the hexagon where interactions are strongest. Transmission (right) can occur at most points in the lattice for high energies but tends to occur at the hexagon center due to the low barrier.
By examining the position within the hexagon where incident atoms are reflected, transmitted, or adsorbed, one can infer the form of the manybody potential at nonsymmetrical parts of the lattice. Figure 4 shows the hexagonlocalized reflection, transmission, and adsorption for several energies. Lattice positions represented in Figure 4 are the turning points for reflection, closest approach positions for transmission, and final xy positions for adsorption. Adsorbed atoms are clustered around the carbon atoms, often showing some lateral vibration.
Reflection is distributed evenly around the perimeter of the hexagon, indicating that incident atoms are deflected away from the hexagon center due to the relatively low force experienced here. Also due to the weak interaction at the hexagon center, it is the most probable location for transmission to occur. Thus, atoms incident upon or deflected toward this position are both able to penetrate. These results agree with those from a previous study [30], which found that reflection occurs at all points in the hexagon, and transmission is most probable near the hexagon center.
The scattering of incident noble gas atoms has been investigated at high energies (keV), where transmitted particles were found to have very little angular deflection while leaving the graphene relatively unaffected [8]. Here, similar results are found for hydrogen at lower energies. Figure 5 shows the angular cross sections of reflection and transmission for our SCCDFTB results. Radii are normalized to unity for the purpose of comparison. Cross sections are calculated according to
Figure 5. Angular distributions of reflected (θ < 90°) and transmitted (θ > 90°) hydrogen atoms. Distributions found to fit the data are shown in black.
Here, 1/N_{max }normalizes the distribution, and the differential solid angle dΩ becomes 2π sin θ dθ due to the azimultal symmetry of the problem. N(θ ± Δθ/2) is the number of atoms scattered into a bin of width Δθ centered at polar angle θ.
Small changes in the x or ycomponents of an atom's linear momentum are much more visible for low incident energies, where these changes can be comparable to the initial momentum. In the SCCDFTB simulations, atoms with such low incident energy tend to reflect when not adsorbed, and the reflected angular distribution shows much more scattering. Transmitting hydrogen atoms in these simulations tend to have higher incident energies, so the small x or yforces don't produce a significant angular displacement of their momenta. While atoms incident at 5 and 10 eV have a wider distribution than at the higher energies, they tend to penetrate only near the center, where the HC interactions are weakest.
The dominance of adsorption in SCCDFTB simulations at impact energies below 1 eV provides enough statistical weight for an investigation of the effects of Hadsorption on the E_{lh }quantity of the affected graphene. However, roughly a third of the incident atoms are found to bond to the surface after initially being reflected at a large angle relative to their initial momenta. These 'wandering' hydrogen atoms, primarily seen at 0.5 eV incidence, generally drift above the graphene surface at a distance of about 3 Å for 2 to 5 fs before falling toward a lattice carbon and adsorbing. Roughly 10% of these 'wanderers' do not bond to a carbon within the simulation time. Therefore, while they are counted as adsorbed in Figure 3, they are ignored in the henceforth analysis to reduce uncertainties.
The graphene band gap is often computed using a band structure or density of states calculation. However, the graphene system studied here is subject to thermal motion as well as bombardment, and the impinging particle should not be included in Brillouin zone integration. As discussed earlier, we simply define a quantity E_{lh }by subtracting the energy of the highest occupied orbital from that of the lowest unoccupied orbital. The 1,000K electronic temperature used creates a 'smearing' of the orbital occupations near the Fermi level. We use occupations of 1.8 for h (analogous to the HOMO) and 0.2 for l (analogous to the LUMO). This allows us to accomplish significant statistics while accounting for the different sites of adsorption and variety of vibrational states in which atom is adsorbed. The system is a 336atom supercell, equivalent to an 18 × 18 × 1 kpoint grid.
Figure 6 shows contour plots of the two equivalent potential wells for hydrogen, corresponding to two adjacent C atoms of graphene. The depth of the wells is about 0.61 eV. Thus, when the kinetic energy of H is comparable to the well depth, excited vibrational motion is possible after adsorption; to account for this, we average the change in E_{lh }over a number of time steps at the end of the simulation. When doing this timeaveraging, it is important to avoid time steps at which some of the hydrogen atoms have not yet bound to the graphene surface. Figure 7 displays the standard deviation of the hydrogen zposition distribution averaged over 1, 4, 12, and 24 fs of simulation time. In all cases, the mean value is within a single standard deviation of 1.2 Å.
Figure 6. A contour plot of the potential energy of a Hatom. In vicinity of the two adjacent carbon bonding centers (C) in graphene, Z being the direction orthogonal to the graphene. The depths of the wells in which hydrogen bonds are equal.
Figure 7. Standard deviation of hydrogen zposition distribution as a function of averaging time.
One can see that the standard deviation, representative of the distribution's width, is higher for low averaging times. Additionally, atoms incident with higher kinetic energy are adsorbed with greater vibrational energy, so they display a wider distribution of zpositions. As shown in Figure 7, the wider distributions that come with this higher vibrational energy produce a smaller change in the E_{lh }on average.
Previous studies [31,32] have found that a tuning of the graphene band structure can be achieved by partial or full hydrogenation of nanoribbons, achieving band gaps of 0.43 to approximately 4.0 eV. The results obtained here support this, showing a sensitivity of the graphene band gap to even a single stuck hydrogen atom. At the largest averaging time considered here, the average change in E_{lh }is 171.5 meV for 0.1 eV incidence, 165.1 meV for 0.2 eV incidence, and 157.7 meV for 0.5 eV incidence (Figure 8). As discussed above, E_{lh }is not equal to the band gap, though it is correlated with it. There is a nonlinear relationship between the change in E_{lh }and zposition of hydrogen, which is the source of the difference between these results.
Figure 8. Mean change in E_{lh }as a function of averaging time for three incident kinetic energies.
Figure 9 shows the change in the E_{lh }quantity as a function of adsorbate distance for an ideal graphene plane. Since the hydrogen is directly above a lattice carbon, the minimum of the potential well is located at 1.2 Å. The average minima and maxima of lowenergy (0.1 eV incidence) and highenergy (0.5 eV incidence) oscillations are 0.1 and 0.3 Å, respectively. The change in E_{lh }decreases with increasing adsorbate distance, and the flattening observed below 1 Å causes the minimum position in oscillation to affect the gap less than the maximum position. Thus, larger oscillation amplitudes produce, on average, a smaller change in the E_{lh }quantity. These averages are shown as thick dashed lines and have a difference of 30 meV, which is on the order of the 14 meV difference between average changes induced by 0.1 and 0.5 eV bombardments.
Figure 9. Mean change in E_{lh }as a function of adsorbate hydrogen distance. Displayed are maximum, minimum, and average changes for typical large and small oscillation amplitudes resulting from 0.5 and 0.1 eV bombardments, respectively. Calculations are performed using an ideal graphene plane.
Comparison with classical molecular dynamics
In the classical molecular dynamics approach, the physical accuracy of the simulation is determined mainly by the quality of the interatomic potentials. Like its predecessor, the REBO potential, AIREBO is a member of the classical bondorder family of potentials [20,21] of the TersoffBrenner type, which provides a good description of the covalent bonds for nonpolar systems. The REBO potential is short ranged (< 2 Å) and, therefore, considerably less costly to use in computation but might not be suitable for collisions where longrange interactions are important, or for describing the coupling of adjacent graphene planes. REBO is also known for its poor treatment of conjugated couplings [20]. The AIREBO contains improved descriptions of the torsional and longrange van der Waals interactions (< 11 Å) as well as improved bonding interactions. The ability to use a classical (if reactive) molecular dynamics approach for the bombardment problem is highly desirable since these approaches are orders of magnitude computationally cheaper than even SCCDFTB.
Figure 10 shows a comparison of the refitted [22] AIREBO and REBO Hgraphene potentials with that of SCCDFTB shown in Figure 2. The CMD potentials were calculated using a 480atom graphene cluster, i.e., no periodic boundary conditions. All three potentials display a well near 1.2 Å for incidence upon a lattice carbon, and the subsequent repulsive barriers agree well as they are all fit to ZBL [22]. However, REBO and AIREBO predict 0.5 and 1.0 eV barriers, respectively, before the potential wells. As these barriers are not present in the SCCDFTB potential, it is expected that REBO and AIREBO result in different dynamics at lowenergy bombardment. The dissimilarities are even more distinct for the other positions in the lattice. AIREBO predicts a potential barrier of over 20 eV, peaking at about 1.35 Å, for incidence on the CC bond center. Neither DFTB nor REBO agree with this barrier, which is produced by the longrange LennardJones terms in AIREBO, since the AIREBO and REBO results are indistinguishable at distances less than 1 Å. Another peak of 20 eV height is found at z = 0 for AIREBO and REBO, only 2 eV higher than the corresponding SCCDFTB curve. REBO is consistently 2 to 5 eV more repulsive than DFTB, but qualitatively very similar. The most distinctive difference between these three potentials is their treatment of the graphene πorbitals. Clearly, the AIREBO LennardJones interactions coming from the six adjacent carbons produce a potential barrier at the graphene hexagon center that is more than 60 eV (525%) higher than the potential in its predecessor, which is in turn roughly 10 eV (380%) higher than DFTB. The REBO potentials clearly agree much more with the DFTB calculations than those of AIREBO, which indicates that the LennardJones interactions which produce the observed potential barriers likely overestimate the hydrogengraphene interaction.
Figure 10. Potential energy of the Hgraphene interaction at canonical points in the lattice. As calculated in DFTB (solid), AIREBO (single dash), and REBO (double dash).
Ito et al. [6,30], Nakamura et al. [5], and Saito et al. [33] have done a comprehensive study of the response of a single graphene sheet (reflection, transmission, and absorption) to the impact of hydrogen atoms and its isotopes in an energy range below 200 eV. They used classical molecular dynamics simulations with the shortrange (< 2 Å) modified Brenner (REBO) potential. Unfortunately, the distribution of barriers and wells is not clear for their modification of the potential; however, their calculation of the reflection, transmission, and adsorption probabilities upon normal impact shows good qualitative agreement with our AIREBO calculations, as illustrated in Figure 11. Classical MD AIREBO and REBO results are significantly different than the quantumclassical SCCDFTB results mainly due to the presence of the potential barriers observed in Figure 10. The AIREBO and REBO graphene potentials are more repulsive than those of SCCDFTB, which results in a 15 eV higher threshold for transmission to occur. However, all methods converge to 100% transmission at high energies, as has been observed in previous studies [8]. While REBO shows a much higher peak in adsorption probability, the presence of the aforementioned barriers in both potentials result in a dominance of reflection at low energies, inconsistent with the results of SCCDFTB presented above.
Figure 11. Probabilities of reflection, transmission, and adsorption as calculated by AIREBO and REBO [5,30]. The presence of potential barriers before potential wells (see Figure 10) results primarily in reflection at low incident energies.
Another result of the increased repulsiveness of AIREBO is the occurrence of physical carbon sputtering upon impact of hydrogen. Figure 12 shows the sputtering yield as a function of incident kinetic energy for AIREBO calculations. If E_{d }is a carbon atom displacement energy from the rapheme, then the kinetic energy of the impact atom in the headon binary collision is The known energy for displacing one atom from a pristine rapheme is 22.2 eV, which yields (H) = 78.2 eV. Consistently, the sputtering yields in Figure 12 for all sputtered species are zero at 20 eV and start rising from 50 eV impact energy (where the sputtering yield is approximately 0.002). After a peak of 0.0325 at about 200 eV, they then decrease with higher incident energy. Chemical sputtering, i.e., production of CH, which is a secondorder process (breaking of a carbon bond followed by capture of H by the carbon atom) here, is quite improbable, and its yield stays well below 0.005. We note again that no sputtering, physical or chemical, is measured in the SCCDFTB simulations in the considered range of impact H energies (< 200 eV).
Figure 12. Sputtering yields of C and CH as determined by AIREBO simulations.
Lastly, the CMD calculations result in larger angular scattering effects than SCCDFTB, as can be seen in Figure 13. The reason for the markedly different distribution is again in the potential barriers that arise from the LennardJones interactions. The stronger interaction produces a more significant change in the incident hydrogen x or ymomentum, resulting in a cos(θθo) distribution with maximum at roughly 37°. Transmitting atoms also interact with the potential barriers on the opposite side of the surface, which are responsible for deflecting these atoms and producing the much wider distribution than observed in the SCCDFTB calculations.
Figure 13. Angular distributions of reflected (θ < 90°) and transmitted (θ > 90°) hydrogen atoms. As calculated in AIREBO CMD simulations.
Conclusions
Understanding the effects of irradiation is paramount in developing graphenebased nanosensors and nanoelectronics. Thus, in this work, simulations of singlelayer graphene bombarded by hydrogen atoms for a wide range of incident energies were carried out using quantumclassical molecular dynamics based on the selfconsistentchargedensity functional tight binding method for treatment of the electron dynamics, combined with classical dynamics of the nuclei. The effects of this bombardment on the graphene sheet and the scattered particle distributions were analyzed in terms of reflection, transmission, and adsorption probabilities and angular distributions. Particularly significant effects of adsorption on the graphene E_{lh }quantity, analogous to the HOMOLUMO gap in clusters, were investigated, predicting a notable change of the graphene electrical conductivity for even one Hatom chemisorbed. Adsorption was found to be the dominant process below 1 eV, with transmission dominating above 20 eV and reflection dominating at the intermediate energies. Reflection was found to have a more significant scattering effect than transmission.
A comparison between results of the SCCDFTB simulations and classical MD simulations employing the AIREBO potential was made, showing a significant difference in the calculated probabilities and chemistry, mainly caused by differences in the multibody potentials. The AIREBO Hgraphene potential overestimates (in comparison to SCCDFTB) the interaction at the hexagon center (πorbital) and CC bond center (σorbital) lattice positions. A comparison of REBO and AIREBO showed that the overestimate is a result of the LennardJones terms in AIREBO. The effect of this added repulsiveness permeated all of the dynamics, producing wider scattering, a much smaller adsorption probability, and nonzero sputtering yields. Refitting of these terms may significantly improve the accuracy of AIREBO.
Changes in the graphene E_{lh }quantity, qualitatively associated to the HL gap and electric conductance of graphene, were found to depend on incident atom energy. Using an averaging time of 32 fs, in addition to averaging over all adsorbed trajectories, the adsorption effect on the E_{lh }differed by roughly 10 meV between incident energies. By virtue of higher vibrational energy, larger incident kinetic energies are found to have a smaller effect on the band gap, as shown in Figure 9. Further characterization of the E_{lh }changes and/or adsorbed vibrational modes could support the application of graphene in a hypersensitive slow singleparticle detector in agreement with the sensitivity to a single biomolecule being coupled to a graphene sheet [34]. This hypersensitivity of the E_{lh }quantity to hydrogen adsorption indicates that the functionality of graphenebased nanoelectronics could be adversely affected by the irradiation by light, chemically reactive species.
Competing interests
The authors declare that they have no competing interests.
Authors' contributions
RCE carried out computations, and analyzed results with PSK and JD. RCE and PSK prepared the manuscript, which was finalized together with JD, JJ, and PRCK. JJ provided help in using DFTB. All authors read and approved the final manuscript.
Acknowledgements
We acknowledge the support of the Offices of Fusion Energy Sciences (US DOE) (PSK) and the ORNL LDRD program (PSK and JD) as well as the ORISE SULI program (RCE). We acknowledge the support of the TG NSF program for the use of the NICS computer facilities (Kraken). JJ acknowledges the support by the SC/TNEPSCoR grant. Research by PRCK was supported by the Scientific User Facilities Division, U.S. Department of Energy.
References

Dimitrakakis G, Tylianakis E, Froudakis G: Pillared graphene: a new 3D innovative network nanostructure for enhanced hydrogen storage.
Nano Lett 2008, 8(10):31663170. PubMed Abstract  Publisher Full Text

Deretzis I, Fiori G, Iannaccone G, Piccitto G, La Magna A: Quantum transport modeling of defected graphene nanoribbons.
Physica E, in press.
doi: 10.1016/j.physe.2010.06.024

Gorjizadeh N, Kawazoe Y: Chemical functionalization of graphene nanoribbons.

Schedin F, Geim AK, Morozov SV, Hill EW, Blake P, Katsnelson MI, Novoselov KS: Detection of individual gas molecules adsorbed on graphene.
Nat Mater 2007, 6:652655. PubMed Abstract  Publisher Full Text

Nakamura H, Takayama A, Ito A: Molecular dynamics simulation of hydrogen isotope injection into graphene.
Contrib Plasma Phys 2008, 48:265269. Publisher Full Text

Ito A, Nakamura H: Molecular dynamics simulation of bombardment of hydrogen atoms on graphite surface.

Krasheninnikov A, Nordlund K: Ion and electron irradiationinduced effects in nanostructured materials.
J Appl Phys 2010, 107:071301. Publisher Full Text

Lehtinen O, Kotakoski J, Krasheninnikov AV, Tolvanen A, Nordlund K, Keinonen J: Effects of ion bombardment on a twodimensional target: atomistic simulations of graphene irradiation.

Wakabayashi K, Takane Y, Yamamoto M, Sigrist M: Electronic transport properties of graphene nanoribbons.
N J Phys 2009, 11:095016. Publisher Full Text

Porezag D, Frauenheim T, Kohler T, Seifert G, Kaschner R: Construction of tightbindinglike potentials on the basis of densityfunctional theory: application to carbon.
Phys Rev B 1995, 51:1294712957. Publisher Full Text

Elstner M, Porezag D, Jungnickel G, Elsner J, Haugk M, Frauenheim T, Suhai S, Seifert G: Selfconsistentcharge densityfunctional tightbinding method for simulations of complex materials properties.
Phys Rev B 1998, 58:72607268. Publisher Full Text

Oliviera AF, Seifert G, Heine T, Duarte HA: Densityfunctional based tightbinding: an approximate DFT method.
J Braz Chem Soc 2009, 20(7):11931205. Publisher Full Text

Rauls E, Elsner J, Gutierrez R, Frauenheim T: Stoichiometric and nonstoichiometric (1010) and (1120) surfaces in 2HSiC: a theoretical study.
Solid State Comm 1999, 111(8):459464. Publisher Full Text

Ziegler JF, Biersack JP, Littmark U: The Stopping and Range of Ions in Matter. New York: Pergamon; 1985.

Androitis AN, Menon M, Srivastava D, Froudakis G: Extreme hydrogen sensitivity of the transport properties of singlewall carbonnanotube capsules.

Berashevich J, Chakraborty T: Tunable band gap and magnetic ordering by adsorption of molecules on graphene.

Gao H, Wan L, Zhao J, Ding F, Lu J: Band gap tuning of hydrogenated graphene: H coverage and configuration dependence.
J Phys Chem C 2011, 115(8):32363242. Publisher Full Text

Elias DC, Nair RR, Mohiuddin TMG, Morozov SV, Blake P, Halsall MP, Ferrari AC, Boukhvalov DW, Katsnelson MI, Geim AK, Novoselov KS: Control of graphene's properties by reversible hydrogenation.
Science 2009, 323:610630. PubMed Abstract  Publisher Full Text

McKay H, Wales DJ, Jenkins SJ, Verges JA, de Andres PL: Hydrogen on graphene under stress: molecular dissociation and gap opening.

Brenner DW, Shenderova OA, Harrison JA, Stuart SJ, Ni B, Sinnott SB: A secondgeneration reactive empirical bond order (REBO) potential energy expression for hydrocarbons.
J Phys Condens Matter 2002, 14:783802. Publisher Full Text

Stuart SJ, Tutein AB, Harrison JA: A reactive potential for hydrocarbons with intermolecular interactions.
J Chem Phys 2000, 112:64726486. Publisher Full Text

Kent PRC, Dadras J, Krstic PS: Improved hydrocarbon potentials for sputtering studies.
J Nucl Mater 2011, 415:S183S186. Publisher Full Text

Yang M, Nurbawono A, Zhang C, Feng YP, Ariando : Twodimensional graphene superlattice made with partial hydrogenation.
App Phys Lett 2010, 96:193115. Publisher Full Text

Vosko SH, Wilk LH, Nussair M: Accurate spindependent electron liquid correlation energies for local spin density calculations: a critical analysis.
Can J Phys 1980, 58(8):12001211. Publisher Full Text

Jeloaica L, Sidis V: DFT investigation of the adsorption of atomic hydrogen on a clustermodel graphite surface.
Chem Phys Lett 1999, 300:157162. Publisher Full Text

Jakowski J, Irle S, Mrokuma K: Collisioninduced fusion of two C_{60 }fullerenes: quantum chemical molecular dynamics simulations.

Zheng G, Irle S, Morokuma K: Performance of the DFTB method in comparison to DFT and semiempirical methods for geometries and energies of C20C86 fullerene isomers.
Chem Phys Lett 2005, 412:210216. Publisher Full Text

Zheng G, Lundberg M, Jakowski J, Morokuma K: Implementation and benchmark tests of the DFTB method and its application to the ONIOM method.
Int J Quantum Chem 2009, 109:18411854. Publisher Full Text

Elstner M: The SCCDFTB method and its application to biological systems.

Ito A, Nakamura H, Takayama A: Molecular dynamics simulation of the chemical interaction between hydrogen atom and graphene.
J Phys Society Japan 2008, 77:114602. Publisher Full Text

Zhou J, Wu MM, Zhou X, Sun Q: Tuning electronic and magnetic properties of graphene by surface modification.
App Phys Lett 2009, 95(10):103108. Publisher Full Text

Saito S, Ito A, Nakamura H: Incident angle dependence of reactions between graphene and hydrogen atom by molecular dynamics simulation.
Annual Report of National Institute for Fusion Science 2010., 958

Nelson T, Zhang B, Prezhdo OV: Detection of nucleic acids with graphene nanopores: ab initio characterization of a novel sequencing device.
Nano Lett 2010, 10:32373242. PubMed Abstract  Publisher Full Text