Thermal Rectification in Ultra-narrow Hydrogen Functionalized Graphene: A Non- Equilibrium Molecular Dynamics Study

doi:10.21203/rs.3.rs-1743864/v1

Download PDF

Research Article

Thermal Rectification in Ultra-narrow Hydrogen Functionalized Graphene: A Non- Equilibrium Molecular Dynamics Study

https://doi.org/10.21203/rs.3.rs-1743864/v1

This work is licensed under a CC BY 4.0 License

You are reading this latest preprint version

In this study, the non-equilibrium molecular dynamics simulation (NEMD) has been used to evaluate the thermal properties, especially the rectification of ultra-narrow edge-functionalized graphene with hydrogen atoms. The system's small width equals 4.91 Å (equivalent to two hexagonal rings). The dependence of the thermal rectification on the mean temperature, hydrogen concentration, and temperature difference between the two baths was investigated. Results reveal that the thermal rectification increases up to 100% at 550 K by increasing the mean temperature. Also, it is disclosed that hydrogen concentration plays a vibrant role in thermal rectification. Furthermore, a thoroughgoing rectification is obtained in the half fully hydrogenated system due to maximum phonon scattering at the interface. Furthermore, the effects of temperature difference of baths ΔT on thermal rectification has been calculated. As a result, the thermal rectification decreases even though the current heat increases with ΔT. Finally, the thermal resistance at the interface using a mismatching factor between the two-phonon density of states (DOS) on both sides of the interface has been explained.

Thermal rectification

ultra-narrow graphene

edge-functionalized

molecular dynamics

Graphene [1], a two-dimensional (2D) monolayer sheet with honeycomb form, has drawn more attention in the past decade due to its outstanding properties. Graphene shows excellent optical and electronic [2,3], mechanical [4], and thermal properties [5]. Besides graphene, novel 2D materials have been discovered or synthesized functionalized graphene for different purposes, such as carbon nitride [6], molybdenum disulfide [7], and borophene [8], phosphorene [9], graphene [10], and so on. Functional atoms or groups can affect physical properties. For example, graphane was theoretically predicted [10] and synthesized [11,12] by hydrogenating both sides of the graphene in an alternating manner. In this situation, the sp² bond type is transformed to sp³ and affects all properties. In addition, the number of hydrogen atoms on graphane can be changed to control properties.

The hybrid of 2D monolayers such as graphene and graphane (or functionalized graphene with an arbitrary number of hydrogen atoms) is also enjoyable. Thermal rectification in such systems was investigated [13,14] recently. In thermal rectifier (thermal diode) systems, the heat current across the system depends significantly on the thermal gradient direction. Considering the vital role of electrical diodes in electronics, thermal diodes have substantial implications for diverse thermal management issues, ranging from nanoscale calorimeters to microelectronic processors to macroscopic refrigerators and energy-saving buildings [15]. Therefore, the heat current prefers one direction to another one is observable. Rajabpour et al. [13] have studied the graphene and graphane junction to estimate thermal rectification. They obtained thermal rectification of about 20% through NEMD simulation. Moreover, they also examined thermal properties at the interface, which is responsible for phonon scattering, leading to thermal resistance (known as Kapitza resistance). They also showed that the Kapitza resistance depends on the temperature gradient, which causes the heat current to prefer one direction to the opposite.

Shavikloo et al. [14] studied thermal rectification in graphene and graphane with grain boundary interfaces between two sides. They examined thermal rectification by randomly removing hydrogen atoms from graphane. They found that the thermal rectification is less than 10% at room temperature. Furthermore, they investigated the effect of the hydrogen concentration on thermal rectification in the mentioned hybrid system. Removing hydrogen help us to tune the thermal properties so that they have succeeded in increasing thermal rectification 17% by removing 40% of hydrogen atoms than fully hydrogenated graphane. In the two above studies, the system's width was about 10 nm, and periodic boundary conditions were applied to in-plane directions. The width and boundary condition can be influenced the thermal properties. Functionalization with hydrogen has been used in other systems such as single-walled carbon nanotube [16], nanoporous graphene [17], and pristine graphene [18] with particular patterns. Moreover, thermal rectification has also been seen in other systems. In graphene-C₃N junction [19,20], graphene-C₃B [21], graphene-pillared [22], radial rectification [23], telescopic silicon nanowire [24], and so on. Thermal rectification in ultra-narrow graphene with functionalized edge via NEMD simulation has been investigated. The dependency of thermal rectification on temperature and the temperature difference between baths and hydrogen concentration was studied. Additionally, the phonon density of states was obtained to explain thermal resistance at the interface.

In semiconductors, the contribution of electrons in the thermal transport is minor than the phonon's contribution. Therefore, we ignore the electrons and compute phonon thermal current in our ultra-narrow system. This study implemented all NEMD simulations using a Large-scale Atomic/Molecular Massively Parallel Simulator (LAMMPS) [25] as a molecular dynamics program. The desired structure is indicated in Fig. 1 with a length of 20 nm and a width of 4.91 Å. As seen, both edges in the left half of the structure are fully functionalized with hydrogen atoms, while the right side has no hydrogen.

The hydrogen concentration can be different in some simulations due to investigating the effect of that on thermal rectification. Periodic boundary conditions along the x-direction, but free boundary conditions were assumed in other directions have been applied. Moreover, to model the interactions in this system between C-C and C-H, we employed the Adaptive Intermolecular Reactive Empirical Bond Order (AIREBO) [26] potential with a cutoff 3.0 (σ scale factor). This potential is a three-body type and is appropriate for modeling breaking and creating bonds and the description of hydrocarbon systems. Velocity Verlet algorithm is used to integrate the equation of motion of atoms with a time step of 0.2 fs. The time step was considered small because we have hydrogen atoms in this system that can oscillate with high frequency.

First of all, to relax the system and eliminate all extra stress, the NPT ensemble with Nose-Hoover [27] barostat and thermostat was applied to the system along x-direction so that the pressure tensor element P_xx becomes zero at room temperature. In the mentioned situation, the system can expand or shrink if potential desired. After relaxing, as represented in Fig. 1, two slabs at two ends of the system were considered fixed atoms so that the particles inside of them could not move in during simulation time. Moreover, two other slabs near two fixed layers were selected as thermal baths (hot and cold sources). We used an NVT ensemble with a Nose-Hoover thermostat to keep the temperature fixed in thermal sources. The NVE ensemble was applied elsewhere.

Here, to estimate the thermal rectification, we use the equation below [17],

$TR \left(\%\right)= \frac{{Q}_{L\to R}-{Q}_{R\to L}}{{Q}_{R\to L}}\times 100$	(1)
where ${Q}_{L\to R}$ (${Q}_{R\to L}$) is the heat current from left to right (right to the left). We assume that in Eq. 1, ${Q}_{L\to R} > {Q}_{R\to L}$.

Moreover, the phonon density of states (DOS) on both sides of the interface was calculated using the Fourier transform of the velocity autocorrelation function. Therefore, we choose two slabs on both sides of the interface and use the equation below [28,29],

$DOS\left(\omega \right)= \sum _{s}\frac{{m}_{s}}{{k}_{B}T}{\int }_{0}^{\infty }{e}^{-i\omega t}<{v}\left(0\right).{v}\left(t\right){>}_{s}dt$

(2)

Where s is the atom type and symbol < > indicates the velocity autocorrelation function. Parameters of m, k_B and T are the mass, Boltzmann constant, and system temperature. The $\omega$ is the angular frequency.

To obtain heat current across the system (along x-direction), the accumulative energies added to and removed from the heat baths were calculated (see Fig. 2). Indeed, the slope of the linear fit to energy profile $dE/dt$ is the heat current. As indicated in Fig. 2, the input and output energies from the baths and their slope are equal. It shows that the heat current reached to steady state.

The heat current depends significantly on the temperature gradient sign because there is a mass asymmetry across the system. It causes the phonons to prefer to move more in one direction than the opposite one, usually from thick to thin.

Moreover, we present the sensibility of the heat current and the thermal rectification of the mean temperature of the system. Therefore, we enhance the temperature from 200 to 650 K.

As shown in Fig. 3(a), the heat current decreases with increasing the temperature. The behavior of the heat current in both temperature gradient biases arises from two phenomena. The phonon-phonon scattering will increase by increasing temperature, which is a negative factor for heat current. On the other hand, by increasing temperature, the high frequency (high energy) phonons will excite, leading to increased heat current (positive factor). Competition between these two factors specifies the behavior of the heat current. Here, the simulation results indicate that the phonon-phonon scattering has been dominant. Furthermore, the heat current along the positive direction is greater than the negative one. It is the origin of thermal rectification.

In Fig. 3(b), the thermal rectification increases with increasing temperature. It should be noted that the decreasing rate of negative heat current than positive leads to increase thermal rectification. The magnitude of thermal rectification is large and greater than 30% in this case. The estimated trend in Fig. 3(b) was obtained via two fitted lines in Fig. 3(a). As long as the heat flux behavior is linear, the estimated trend remains valid. In analogy with other studies by Rajabpour [13], Shavikloo [14], Yousefi [17], Cartoixà [24], we get a higher thermal rectification. Therefore, we introduce a promising thermal rectifier.

Interestingly, thermal rectification under various conditions explored, for example, different ΔT. This parameter can be changed under application circumstances. Therefore, we check its effect on heat current and thermal rectification.

According to Fourier's law ($J=-\kappa \nabla T$), the heat flux is proportional to ΔT. Therefore, we expect that $J$ increases with increasing ΔT. The NEMD result confirms this proportionality, and the heat current increases with increasing ΔT as shown in Fig. 4(a). For the positive direction, the slope of the fitted line is greater than the negative due to thermal conductivity depending on the temperature gradient sign. Although the difference between two heat currents gets larger, the heat current for the negative direction gets more significant too. Therefore, we obtain a decreasing behavior of thermal rectification as indicated in Fig. 4(b). Thermal rectifications are in the range of (~ 20 %, ~ 80%) forthe range of ΔT (20 K, 50 K) that we studied.

One of the crucial factors that help us tune the thermal rectification is hydrogen concentration in the half of the system. According to Figs. 5(a) and (b), we explored the effect of removed hydrogen atoms in the system in 5% − 95%. For any amount of removed hydrogen, the heat current in the positive direction is more significant than the negative. For 95%, due to the dilution of hydrogen atoms on the left side, the heat currents increase in both directions and close to each other because the interface is vanishing. On the other hand, the difference between the heat currents arises from thermal rectification. Figure 5(b) shows that the NEMD result indicates that the thermal rectification suppresses in the mentioned range by removing hydrogen. With a slight decrease in hydrogen (5%), the thermal rectification drops from 40–18%. With more removing the hydrogen atoms up to 50%, thermal rectification increases slightly. This behavior was reported already elsewhere [14]. This shows that the hydrogens are responsible for phonon scattering and thermal rectification. Thermal rectification beyond 65% decreases because the phonon scattering decreases at the interface in a negative direction, leading to an increase in heat current.

To understand the underlying mechanism of the thermal rectification, we should investigate the origin of the thermal resistance at the interface (known as Kapitza) [30,31]. Kapitza originates from phonon scattering at the interface. Therefore, a formalism called interface conductance modal analysis (ICMA) developed by Gordiz and Henry [32,33] applies to knowing the number of phonons is conductance or thermal resistance at the interface. The phonon density of states (DOS) on both sides of the interface in positive and negative temperature gradient signs specify mismatching (see Fig. 6). The mismatching indicates that some phonon modes are present on only one side. This is attributed to the Kapitza resistance at the interface, which depends on thermal gradient bias. Moreover, the mismatching factor between two DOSs on both sides of the interface can be calculated through the equation below [17],

$\psi = \frac{{\int }_{0}^{\infty }{D}_{1}\left(\omega \right){D}_{2}\left(\omega \right)d\omega }{{\int }_{0}^{\infty }{D}_{1}\left(\omega \right)d\omega {\int }_{0}^{\infty }{D}_{2}\left(\omega \right)d\omega }$

(3)

Where ${D}_{1}$ and ${D}_{2}$ are the DOS on the left and the right side, respectively. For positive and negative temperature gradients, the mismatching factor is 0.01253 and 0.01332, respectively, for $T=300 K$ and$\varDelta T=40 K.$

Using a series of NEMD simulations, we reported the thermal rectification in ultra-narrow hydrogen functionalized graphene. Thermal rectification is significantly high and even increases up to 100% in temperature 550 K. The dependency of thermal rectification on mean temperature, ΔT, and hydrogen atoms concentration were investigated. We found that the thermal rectification in the narrow structure in this work increases with increasing mean temperature while then decreases with increasing ΔT. Furthermore, we can tune the thermal rectification between 10%-30% by removing some percent of hydrogen atoms. Finally, by calculating the phonon density of states on both sides of the interface, the underlying mechanism of the phonon scattering at the interface was shown.

Funding: N/A

Conflicts of interest/Competing interests: N/A

Availability of data and material: N/A

Authors' contributions: All authors have the same contribution.

Ethics approval: On behalf of all authors the corresponding author approves.

Consent to participate: On behalf of all authors the corresponding author consents.

Consent for publication: On behalf of all authors the corresponding author consents published under the subscription model.

1. A.K. Geim, K.S. Novoselov, The rise of graphene, Nat. Mater. 6 (2007) 183–191. doi:10.1038/nmat1849.

2. A.H. Castro Neto, F. Guinea, N.M.R. Peres, K.S. Novoselov, A.K. Geim, The electronic properties of graphene, Rev. Mod. Phys. 81 (2009) 109–162. doi:10.1103/RevModPhys.81.109.

3. W. Hu, Z. Li, J. Yang, W. Hu, Z. Li, J. Yang, Electronic and optical properties of graphene and graphitic ZnO nanocomposite Electronic and optical properties of graphene and graphitic ZnO nanocomposite structures, J. Chem. Phys. 124706 (2015) 124706. doi:10.1063/1.4796602.

4. F. Rezaee, F. Yousefi, F. Khoeini, Heat transfer in strained twin graphene: A non-equilibrium molecular dynamics simulation, Phys. A Stat. Mech. Its Appl. 564 (2021) 125542. doi:10.1016/j.physa.2020.125542.

5. A. Kumar, K. Sharma, A.R. Dixit, A review on the mechanical and thermal properties of graphene and graphene-based polymer nanocomposites: understanding of modelling and MD simulation, Https://Doi.Org/10.1080/08927022.2019.1680844. 46 (2019) 136–154. doi:10.1080/08927022.2019.1680844.

6. L.-B. Shi, Y.-Y. Zhang, X.-M. Xiu, H.-K. Dong, Structural characteristics and strain behavior of two-dimensional C3N : First principles calculations, Carbon N. Y. 134 (2018) 103–111. doi:10.1016/j.carbon.2018.03.076.

7. D. Kong, H. Wang, J.J. Cha, M. Pasta, K.J. Koski, J. Yao, Y. Cui, Synthesis of MoS 2 and MoSe 2 Films with Vertically Aligned Layers, Nano Lett. 13 (2013) 1341–1347. doi:10.1021/nl400258t.

8. A.J. Mannix, X.-F. Zhou, B. Kiraly, J.D. Wood, D. Alducin, B.D. Myers, X. Liu, B.L. Fisher, U. Santiago, J.R. Guest, M.J. Yacaman, A. Ponce, A.R. Oganov, M.C. Hersam, N.P. Guisinger, Synthesis of borophenes: Anisotropic, two-dimensional boron polymorphs, Science (80-. ). 350 (2015) 1513–1516. doi:10.1126/science.aad1080.

9. S. Das, M. Demarteau, A. Roelofs, Ambipolar Phosphorene Field Effect Transistor, ACS Nano. 8 (2014) 11730–11738. doi:10.1021/nn505868h.

10. J.O. Sofo, A.S. Chaudhari, G.D. Barber, Graphane: A two-dimensional hydrocarbon, Phys. Rev. B - Condens. Matter Mater. Phys. 75 (2007). doi:10.1103/PhysRevB.75.153401.

11. D.C. Elias, R.R. Nair, T.M.G. Mohiuddin, S. V. Morozov, P. Blake, M.P. Halsall, A.C. Ferrari, D.W. Boukhvalov, M.I. Katsnelson, A.K. Geim, K.S. Novoselov, Control of Graphene's Properties by Reversible Hydrogenation: Evidence for Graphane, Science (80-. ). 323 (2009) 610–613. doi:10.1126/science.1167130.

12. A.K. Singh, B.I. Yakobson, Electronics and Magnetism of Patterned Graphene Nanoroads, Nano Lett. 9 (2009) 1540–1543. doi:10.1021/nl803622c.

13. A. Rajabpour, S.M. Vaez Allaei, F. Kowsary, Interface thermal resistance and thermal rectification in hybrid graphene-graphane nanoribbons: A non-equilibrium molecular dynamics study, Appl. Phys. Lett. 99 (2011) 051917. doi:10.1063/1.3622480.

14. M. Shavikloo, S. Kimiagar, Thermal rectification in partially hydrogenated graphene with grain boundary, a non-equilibrium molecular dynamics study, Comput. Mater. Sci. 139 (2017) 330–334. doi:10.1016/j.commatsci.2017.08.024.

15. C. W. Chang, D. Okawa, A. Majumdar, A. Zettl, Solid-State Thermal Rectifier, Science (80-. ). 314 (2006) 1121. http://www.sciencemag.org/content/314/5802/1121.full.pdf.

16. K. Gordiz, S. Mehdi Vaez Allaei, Thermal rectification in pristine-hydrogenated carbon nanotube junction: A molecular dynamics study, J. Appl. Phys. 115 (2014) 163512. doi:10.1063/1.4873124.

17. F. Yousefi, F. Khoeini, A. Rajabpour, Thermal conductivity and thermal rectification of nanoporous graphene: A molecular dynamics simulation, Int. J. Heat Mass Transf. 146 (2020) 118884. doi:10.1016/j.ijheatmasstransfer.2019.118884.

18. S.K. Chien, Y.T. Yang, C.K. Chen, Influence of hydrogen functionalization on thermal conductivity of graphene: Non-equilibrium molecular dynamics simulations, Appl. Phys. Lett. (2011). doi:10.1063/1.3543622.

19. O. Farzadian, C. Spitas, K. V Kostas, Graphene-carbon nitride interface-geometry effects on thermal rectification: a molecular dynamics simulation, Nanotechnology. 32 (2021) 215403. doi:10.1088/1361-6528/abe786.

20. O. Farzadian, A. Razeghiyadaki, C. Spitas, K. V Kostas, Phonon thermal rectification in hybrid graphene-$\mathrm{C_{3}N}$: a molecular dynamics simulation, Nanotechnology. 31 (2020) 485401. doi:10.1088/1361-6528/abb04b.

21. L. Kiani, J. Hasanzadeh, F. Yousefi, P.A. Anaraki, Phonon modes contribution in thermal rectification in graphene-C3B junction: A molecular dynamics study, Phys. E Low-Dimensional Syst. Nanostructures. (2021) 114724. doi:10.1016/j.physe.2021.114724.

22. F. Yousefi, F. Khoeini, A. Rajabpour, Thermal rectification and interfacial thermal resistance in hybrid pillared-graphene and graphene: a molecular dynamics and continuum approach, Nanotechnology. 31 (2020) 285707. doi:10.1088/1361-6528/ab8420.

23. F. Yousefi, M. Shavikloo, M. Mohammadi, Non-equilibrium molecular dynamics study on radial thermal conductivity and thermal rectification of graphene, Mol. Simul. 45 (2019) 646–651. doi:10.1080/08927022.2019.1578354.

24. X. Cartoixà, L. Colombo, R. Rurali, Thermal Rectification by Design in Telescopic Si Nanowires, Nano Lett. 15 (2015) 8255–8259. doi:10.1021/acs.nanolett.5b03781.

25. S. Plimpton, Fast Parallel Algorithms for Short-Range Molecular Dynamics, J. Comput. Phys. 117 (1995) 1–19. doi:10.1006/jcph.1995.1039.

26. S.J. Stuart, A.B. Tutein, J.A. Harrison, A reactive potential for hydrocarbons with intermolecular interactions, J. Chem. Phys. 112 (2000) 6472–6486. doi:10.1063/1.481208.

27. S. Melchionna, G. Ciccotti, B. Lee Holian, Hoover NPT dynamics for systems varying in shape and size, Mol. Phys. 78 (1993) 533–544. doi:10.1080/00268979300100371.

28. M. Khalkhali, A. Rajabpour, F. Khoeini, Thermal transport across grain boundaries in polycrystalline silicene: A multiscale modeling, Sci. Rep. 9 (2019) 5684. doi:10.1038/s41598-019-42187-w.

29. B. Li, J. Lan, L. Wang, Interface thermal resistance between dissimilar anharmonic lattices, Phys. Rev. Lett. 95 (2005) 1–5. doi:10.1103/PhysRevLett.95.104302.

30. H. Budd, J. Vannimenus, Thermal Boundary Resistance, Phys. Rev. Lett. 26 (1971) 1637–1640. doi:10.1103/PhysRevLett.26.1637.

31. A. Rajabpour, S. Volz, Thermal boundary resistance from mode energy relaxation times: Case study of argon-like crystals by molecular dynamics, J. Appl. Phys. 108 (2010) 094324. doi:10.1063/1.3500526.

32. K. Gordiz, A. Henry, A formalism for calculating the modal contributions to thermal interface conductance, New J. Phys. 17 (2015). doi:10.1088/1367-2630/17/10/103002.

33. K. Gordiz, A. Henry, Phonon transport at interfaces: Determining the correct modes of vibration, J. Appl. Phys. 119 (2016). doi:10.1063/1.4939207.

No competing interests reported.

Download PDF

Editorial decision: Major revision
12 Aug, 2022
Reviews received at journal
11 Aug, 2022
Reviews received at journal
26 Jul, 2022
Reviewers agreed at journal
17 Jul, 2022
Reviewers invited by journal
29 Jun, 2022
Editor assigned by journal
17 Jun, 2022
Submission checks completed at journal
17 Jun, 2022
First submitted to journal
10 Jun, 2022

You are reading this latest preprint version

Thermal Rectification in Ultra-narrow Hydrogen Functionalized Graphene: A Non- Equilibrium Molecular Dynamics Study

Status:

Version 1

Abstract

Figures

Introduction

Computational Method

Results And Discussion

Conclusion

Declarations

References

Additional Declarations

Status:

Version 1

\(TR \left(\%\right)= \frac{{Q}_{L\to R}-{Q}_{R\to L}}{{Q}_{R\to L}}\times 100\)	(1)
where \({Q}_{L\to R}\) (\({Q}_{R\to L}\)) is the heat current from left to right (right to the left). We assume that in Eq. 1, \({Q}_{L\to R} > {Q}_{R\to L}\).