Next Article in Journal
A Chemical Damage Creep Model of Rock Considering the Influence of Triaxial Stress
Previous Article in Journal
Selective Laser Melting and Mechanical Properties of Stainless Steels
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Thermal Conductance of Copper–Graphene Interface: A Molecular Simulation

1
State Key Laboratory of Nonferrous Metals and Processes, GRINM Group Co., Ltd., Beijing 100088, China
2
GRIMAT Engineering Institute Co., Ltd., Beijing 101407, China
3
General Research Institute for Nonferrous Metals, Beijing 100088, China
*
Author to whom correspondence should be addressed.
Materials 2022, 15(21), 7588; https://doi.org/10.3390/ma15217588
Submission received: 28 September 2022 / Revised: 14 October 2022 / Accepted: 19 October 2022 / Published: 28 October 2022

Abstract

:
Copper is often used as a heat-dissipating material due to its high thermal conductivity. In order to improve its heat dissipation performance, one of the feasible methods is to compound copper with appropriate reinforcing phases. With excellent thermal properties, graphene has become an ideal reinforcing phase and displays great application prospects in metal matrix composites. However, systematic theoretical research is lacking on the thermal conductivity of the copper–graphene interface and associated affecting factors. Molecular dynamics simulation was used to simulate the interfacial thermal conductivity of copper/graphene composites, and the effects of graphene layer number, atomic structure, matrix length, and graphene vacancy rate on thermal boundary conductance (TBC) were investigated. The results show that TBC decreases with an increase in graphene layers and converges when the number of graphene layers is above five. The atomic structure of the copper matrix affects the TBC, which achieves the highest value with the (011) plane at the interface. The length of the copper matrix has little effect on the TBC. As the vacancy rate is between 0 and 4%, TBC increases with the vacancy rate. Our results present insights for future thermal management optimization based on copper matrix composites.

1. Introduction

Electronic devices are increasingly becoming more integrated, multi-functional, and smaller as a result of the rapid development of the electronics industry. The increment in power density per unit volume produces a great deal of heat, which needs to be quickly dissipated to avoid negative impacts on the efficiency and longevity of equipment [1,2,3]. By incorporating a reinforcement phase with a low density and low expansion coefficient into the conventional metallic heat dissipation material, a metal matrix composite possesses high thermal conductivity and a low expansion coefficient as a typical thermal management composite material [4,5]. Metal matrix composites (MMCs) have been used for a variety of industrial applications, especially thermal management devices and aerospace materials [6,7,8,9]. Copper is a good composite matrix due to its excellent thermal conductivity, high melting temperature (1083 °C), good chemical stability, and corrosion resistance [10]. By incorporating a suitable reinforcing phase, Cu-based MMCs provide the opportunity to customize copper’s metallic characteristics to satisfy the demands of thermal management. Some carbon materials, such as diamond, graphite, and graphene, are good candidates for the reinforcing phase [5,11,12,13].
Graphene has attracted the attention of many researchers due to its excellent thermal properties [14,15]. In experiments, many studies have been conducted on composite materials prepared with copper as a matrix and graphene as a reinforcing phase. Chu et al. enabled the in-plane thermal conductivity of graphene/copper composites to reach 525 W/mK by adding 35 vol% graphene nanoplatelet [16]. Gao et al. reported that the thermal conductivity of the composite reached a peak of 397 W/mK when the graphene oxide content reached 0.3 wt.% [17]. In another work, Omayma et al. added graphene nanosheets to Cu/WC-TiC-Co [18]. When the graphene nanosheet content increases from 0 to 1wt.%, the thermal conductivity of the composite increases from 190 W/mK to 351 W/mK. Enhancing the heat transfer between adjacent materials is the key to improving the thermal conductivity of composite materials, so the thermal boundary conductance (TBC) between materials cannot be ignored in the study of material thermal properties. Although previous experiments have demonstrated that graphene as a reinforcing phase can effectively improve the thermal conductivity of copper matrix composites, data on the copper/graphene interface are relatively scarce. Measuring TBC is an important method to characterize the interfacial thermal properties of different materials, and it is impacted by defects, atomic structure, and interfacial bonding strength [19,20,21]. Therefore, it is necessary to systematically study the TBC of Cu/graphene. Molecular dynamics have been widely applied to investigate the heat transport process across solid–solid interfaces. It has been shown that interfacial phonon scattering plays an important role in heat transport at the graphene/metal interface at the temperature range from 50 to 500 K [22]. Problems regarding phonon scattering can be well explored utilizing molecular dynamics.
In this work, we have studied the heat transfer characteristics and mechanism of the Cu/Gr interface using the nonequilibrium molecular dynamics (NEMD) method. Table 1 lists the studies of different metal matrix/graphene composites using the NEMD method. The effects of graphene layer number, Cu matrix atomic structure, and graphene vacancy rate on the TBC have been systematically studied. Meanwhile, the phonon coupling factor is used to explain the underlying mechanism.

2. Model and Methods

The large-scale atomic/molecular massively parallel simulator (LAMMPS) program is applied to carry out the molecular dynamics (MD) simulations [26,27,28]. Reliable forcefields are required to achieve accurate MD simulations. However, there is not yet a single forcefield that can be used to describe a system made of Cu and C. In this work, we adopted a hybrid method including the well-known exact forcefields, which describe the Cu–Cu and C–C interactions. To explain forcefields of copper atoms and carbon atoms, respectively, the embedded atom method (EAM) [29,30] and optimized AIREBO potentials [31,32] are used in this work. Meanwhile, the Lennard-Jones (LJ) function V(rij) = 4ε[(σ/rij)12 − (σ/rij)6] with the values σCu−Gr = 0.30825 nm and εCu−Gr = 0.02578 eV describes the non-bonded van der Waals interactions between Cu and C atoms [33,34]. Based on these forcefield descriptions, the constructed models are first minimized to achieve equivalent states and then undergo the NEMD methods to calculate the TBC. Below are the details about model construction, NEMD methods, and underlying phonon spectra.

2.1. Model

Figure 1a depicts the Cu/Gr composite model, which consists of two blocks of face-centered cubic (FCC) Cu crystals sandwiching one n-layer graphene block (n = 1–7). The value n here refers to the number of layers of graphene. Lx, Ly, and Lt stand for the Cu block’s length, breadth, and thickness, respectively. To guarantee tiny lattice mismatch strain imposed in the Cu/Gr interface structure, we chose different periods while similar model lengths for the copper and graphene models (Figure 1c). Graphene is subjected to tensile strains of 0.407% and 0.085% at two corresponding directions, respectively, in the simulation box, which is set as 4.69 nm in both the x- and y-directions. Table 2 lists the model parameters of different Cu crystal planes when the number of graphene layers n = 1. The interlayer gap between graphene layers is 0.34 nm, while the lamellar distance at the neighboring Cu/Gr layers is relaxed to 0.30 nm [35]. Two layers of Cu atoms on each side in the z direction must be fixed. Meanwhile, two layers of atoms adjacent to the fixed layers serve as the heat source (hot) and heat sink (cold) thermostats, respectively.

2.2. NEMD Method

The NEMD method establishes a non-equilibrium thermal conductivity process by exerting the disturbance of the system and then calculates the thermal conductivity through Fourier’s law of thermal conductivity [36,37]. As shown in Figure 1b, the NEMD method involves placing a cold source and a heat source at both the material’s ends, adding heat to the heat source, subtracting equal amount of heat from the cold source, and, finally, calculating the heat transfer and temperature gradient as a result. Simulations were run with periodic boundary conditions carried out along the x and y (in-plane) directions and free boundary conditions in the z (out-of-plane) direction in order to obtain the TBC based on the NEMD method. Using the Langevin thermostat, we adjusted the heat and cold resources to temperatures of 400 K and 200 K, respectively [38]. The imposed temperature difference activates the process of directional heat transfer of the system. A temperature gradient was produced after 0.5 ns, with each component time-independent. Each block of the Cu model was partitioned into 26 equivalent slabs along z direction (the heat flux) to evaluate the temperature distribution (Figure 1b), and the velocity data during the last 1 ns were collected for vibrational analysis. The energy change over time (heat flux) can be calculated as follows: J = d Δ E(t)/dt; here, Δ E means the energy change. The TBC can be obtained from G = J/A Δ T; here, A indicates the cross-section area of heat transfer. Here, we calculate T by averaging the two temperature differences that occurred adjacent to the Cu/Gr interfaces’ heat source and heat sink.

2.3. Phonon Spectra

We studied the vibrational density of states (VDOS) across the interfaces to clarify the process of interfacial heat transfer. First, the velocity autocorrelation function (Vacf) [39] can be obtained from Lammps by analyzing the velocities of each atom in different directions:
Vacf ( t ) = 1 N i = 1 N v i ( 0 ) v i ( t ) ,
Then, the VDOS is calculated through the Fourier transform of the Vacf:
VDOS ( ω ) = e i ω t VACF ( t ) dt ,
The out-of-plane phonon dominates the thermal transfer process between the interfacial graphene and copper layers, and the low-frequency zone is primarily responsible for phonon coupling. S, the phonon spectrum parameter, defines the overlap as:
S = 0 P Gr ( ω ) P Cu ( ω ) d ω 0 P Gr ( ω ) d ω × 0 P Cu ( ω ) d ω ,
where PGr( ω ) and PCu( ω ) denote the phonon spectra at frequency ω of graphene and Cu atoms at the interface, respectively.

3. Results and Discussion

3.1. Effect of Graphene Layer Number

We used the Cu(001)/Gr interface model to study the effect of graphene layer number on TBC. The simulations were performed at 300 K through the NEMD approach, and Figure 1a depicts the associated schematic model. The two adjacent atom layers were set as cold (heat sink) and hot (heat source) thermostats, respectively, and both ends were fixed. As Figure 2 displays, the TBC reaches a maximum value of 640.19 MW/m2K when n = 1. With the increase in n, the TBC gradually decreases until it converges at n = 5. There is a reduction of 78% from 640.19 MW/m2K (n = 1) to 140.44 MW/m2K (n = 5).
The interfacial phonon scattering across the interface enhanced by the increment in graphene thickness results in a reduction in TBC. Similar phenomena can be observed in research on heat transport across the gra/Al, gra/Ti, and gra/Ni interfaces. To clarify the underlying mechanisms, we calculated the vibrational densities of different modes (in-plane and out-of-plane) of the interfacial copper and graphene layers and obtained the corresponding phonon coupling strengths as shown in Figure 3a, respectively. In this model, the VDOS values of graphene and copper are highly overlapped in the low-frequency region (0–10 THz). The results show that the magnitude of overlap factor S (Figure 3b) in the out-of-plane direction is four times higher than the value in the in-plane direction, so out-of-plane phonon vibration of graphene plays a dominant role in TBC. The phonon coupling strength along the out-of-plane direction monotonically decreases with the increase in the number of graphene layers, which displays a similar trend with TBC.

3.2. Effect of the Atomic Structure

We create three distinct crystal surfaces—(001), (011), and (111)—of the Cu crystals coupled with graphene at the interface to assess the effect of distinct atomic structures upon thermal transport across the Cu/Gr interface (Figure 1c). To study the influence of atomic structure on thermal transport across the Cu/Gr interface, we construct three different crystal surfaces—(001), (011), and (111)—of the Cu crystals connected with the graphene sheet at the interface (Figure 1c). The variations in the lattice structure inevitably cause the variation in the bonding scheme at the interface. In order to minimize the strain imposed by lattice mismatch, the cross-sectional area (A) was also varied when these three atomic structures were re-matched with monolayer graphene according to their dimensions in the x- and y-directions. The model parameters with distinct crystal planes of Cu can be observed in Table 2. The simulation results at room temperature of 300 K show that the TBC of Cu(011)/Gr (673.61 MW/m2K) is about 4.2% and 5.2% higher than the values of the Cu(001)/Gr (640.19 MW/m2K) and Cu(111)/Gr (644.41 MW/m2K) models, respectively. Figure 4b displays the VDOS values of carbon and copper atoms in graphene and copper layers on diverse contact interfaces. The change curve of S with TBC is depicted in Figure 4c, and the change trend of TBC is in good agreement with S, demonstrating the variety of TBC in various crystal planes. For example, the higher value of S indicates that more low-frequency phonons are activated to promote the interfacial heat transport process. Compared with the VDOS of single graphene (Figure 4a) under the same simulation conditions, the peak densities of the in-plane and out-of-plane state densities of graphene after being compounded with copper moved to the low-frequency region.
We calculated the potential energy distribution on different atomic structures by shifting the graphene layer on the ground in order to further confirm the veracity of the findings. The calculation results are shown in Figure 5. The difference values between the maximum and minimum potential energy of the (001), (011), and (111) surfaces are 0.8 × 10−3 eV, 1.0 × 10−2 eV, and 1.8 × 10−3 eV, respectively. This means that the energy required to move graphene from a stable position to an unstable position is higher on the (011) surface. Therefore, graphene has the strongest binding ability bonding with the (011) surface. Stronger interfacial bonding is good for phonon transmission, as shown in Figure 4b, which indicates that the (011) model possesses the largest S = 0.0987. Therefore, the effect of the atomic structure can be clarified through the potential energy distribution.
We further explore the TBC by elongating the length of the copper matrix, i.e., setting the length of the copper matrix as 3.6 nm, 5.4 nm, and 7.2 nm, respectively. As shown in Figure 6, the TBC is almost independent of the thickness of Cu block Lt.

3.3. Effect of Graphene Vacancy Rate

The influence of the material defects on thermal conductivity cannot be ignored. According to the synthesis approach, there can be diverse types of defects in graphene. In this work, we take the vacancy defect into consideration. The vacancy defect concentrations were sampled from 0 to 4%, and the models were established every 0.5% for simulation calculation. We produce vacancies by randomly removing a given number of carbon atoms from the graphene structure. The results show that, with the increase in the vacancy defect concentration, the interface thermal conductivity of different models fluctuates but shows an overall upward trend. The vacancy defect concentration of the model with the highest interfacial thermal conductivity is 3.5%, and the value is 977 MW/m2K, which is 52.6% higher than that of the model without vacancy defect (Figure 7a).
Figure 7b shows the VDOS values of copper and carbon atoms. Previous work has shown that the VDOS values of graphene and copper mainly coincide in the low-frequency region (less than 10 THz). With the introduction of vacancies in graphene, the out-of-plane VDOS of graphene gradually shifts towards the low-frequency region, thus increasing the coupling of the out-of-plane vibration between graphene and copper (Figure 7c). Previous studies have also demonstrated that atoms will lose sp2 bonds after removing their neighboring atoms. Consequently, their out-of-plane vibrations will be less constrained than those of pristine graphene and, therefore, possess more low-frequency patterns [40].

4. Conclusions

In summary, the interfacial thermal conductivity of the copper–graphene interface was systematically studied by means of MD simulations. The key findings are concluded as follows. First, the thermal conductivity of the interface decreases with the increase in the number of graphene layers. Compared with single-layer graphene, the thermal conductivity of n-layer (n > 5) graphene decreases by about 78%. The underlying mechanism is that the increase in graphene layers enhances phonon scattering, which weakens the phonon coupling between copper and graphene in the low-frequency region. Second, Cu(011)/Gr possesses the highest TBC (673.61 MW/m2K), which has an increment of 4.2% and 5.2%, compared with the values of the Cu(001)/Gr (640.19 MW/m2K) and Cu(111)/Gr (644.41 MW/m2K) models, respectively. The effect of atomic structures can be clarified through the potential energy distribution, which indicates that a higher energy barrier within the interface means stronger interfacial bonding and is good for phonon transmission. In addition, the interface thermal conductivity increases gradually when the vacancy defect density of graphene increases from 0 to 3.5%. The increasing trend can be attributed to variation in the vibration frequency schemes caused by the fracture between carbon atoms caused around the vacancy defects.
Compared with previous research regarding the Cu/Gr interface, we have investigated the influence of different factors on the TBC more systematically and gained a comprehensive understanding of the thermal conductance of the copper–graphene interface. The obtained findings are helpful regarding thermal transport mechanisms in future application of Cu/Gr MMCs in thermal management materials.

Author Contributions

Software, J.Z.; Writing—original draft, J.Z.; Writing—review & editing, S.H., Z.X., H.G. and H.Y. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Youth Fund Project of GRINM, grant number 66232166.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The raw data cannot be shared at this time as the data are also part of an ongoing study.

Acknowledgments

We acknowledge the support of Youth Fund Project of GRINM and State Key Laboratory of Nonferrous Metals and Processes.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Arden, L.M.; Li, S. Emerging challenges and materials for thermal management of electronics. Mater. Today 2014, 17, 163–174. [Google Scholar] [CrossRef]
  2. Zhang, J.; Shao, D.; Jiang, L.; Zhang, G.; Wu, H.; Day, R.; Jiang, W. Advanced thermal management system driven by phase change materials for power lithium-ion batteries: A review. Renew. Sustain. Energy Rev. 2022, 159, 112207. [Google Scholar] [CrossRef]
  3. Peng, Y.; Cui, Y. Advanced Textiles for Personal Thermal Management and Energy. Joule 2020, 4, 724–742. [Google Scholar] [CrossRef]
  4. Ömer, G.; Nihal, B. A short review on mechanical properties of graphene reinforced metal matrix composites. J. Mater. Res. Technol. 2020, 9, 6808–6833. [Google Scholar] [CrossRef]
  5. Andrey, M.A.; Fedor, M.S.; Andrey, I.A.; Vladimir, I.N. Mechanical properties of a diamond–copper composite with high thermal conductivity. Mater. Des. 2015, 87, 527–539. [Google Scholar] [CrossRef]
  6. Vasanthakumar, P.; Sekar, K. Processing and preparation of aerospace-grade aluminium hybrid metal matrix composite in a modified stir casting furnace integrated with mechanical supersonic vibration squeeze infiltration method. Mater. Today Commun. 2021, 26, 101732. [Google Scholar] [CrossRef]
  7. Baig, M.M.A.; Hassan, S.F.; Saheb, N.; Patel, F. Metal Matrix Composite in Heat Sink Application: Reinforcement, Processing, and Properties. Materials 2021, 14, 6257. [Google Scholar] [CrossRef]
  8. Kafil, M.R.; Eric, D.; Graham, L.W.C.; Anthony, J.R. Present and future thermal interface materials for electronic devices. Int. Mater. Rev. 2018, 63, 1–21. [Google Scholar] [CrossRef]
  9. Zulfiqar, K.; Hafiz, M.A. Air cooled heat sink geometries subjected to forced flow: A critical review. Int. J. Heat Mass Transf. 2019, 130, 141–161. [Google Scholar] [CrossRef]
  10. Tai-Fa, Y. Fabrication and thermal analysis of a copper/diamond/copper thermal spreading device. Surf. Coat. Technol. 2007, 202, 1208–1213. [Google Scholar] [CrossRef]
  11. Cao, H.J.; Tan, Z.Q.; Lu, M.H.; Ji, G.; Yan, X.J.; Di, C.; Yuan, M.Y.; Guo, Q.; Su, Y.S.; Ahmed, A.; et al. Graphene interlayer for enhanced interface thermal conductance in metal matrix composites: An approach beyond surface metallization and matrix alloying. Carbon 2019, 150, 60–68. [Google Scholar] [CrossRef]
  12. Pan, Y.; He, X.; Ren, S.; Wu, M.; Qu, X. Optimized thermal conductivity of diamond/Cu composite prepared with tungstencopper-coated diamond particles by vacuum sintering technique. Vacuum 2018, 153, 74–81. [Google Scholar] [CrossRef]
  13. Zhang, H.; Chao, M.; Zhang, H.; Tang, A.; Ren, B.; He, X. Microstructure and thermal properties of copper matrix composites reinforced by chromium-coated discontinuous graphite fibers. Appl. Therm. Eng. 2014, 73, 739–744. [Google Scholar] [CrossRef]
  14. Wei, N.; Chen, Y.; Cai, K.; Zhang, Y.Y.; Pei, Q.X.; Zheng, J.C.; Mai, Y.W.; Zhao, J.H. Unusual thermal properties of graphene origami crease: A molecular dynamics study. Green Energy Environ. 2022, 7, 86–94. [Google Scholar] [CrossRef]
  15. Rafiee, M.A.; Rafiee, J.; Wang, Z.; Song, H.H.; Yu, Z.Z.; Koratkar, N. Enhanced Mechanical Properties of Nanocomposites at Low Graphene Content. Appl. Phys. Lett. 2018, 3, 3884–3890. [Google Scholar] [CrossRef]
  16. Chu, K.; Wang, X.-H.; Li, Y.-B.; Huang, D.-J.; Geng, Z.-R.; Zhao, X.-L.; Liu, H.; Zhang, H. Thermal properties of graphene/metal composites with aligned graphene. Mater. Des. 2018, 140, 85–94. [Google Scholar] [CrossRef]
  17. Gao, X.; Yue, H.; Guo, E.; Zhang, H.; Lin, X.; Yao, L.; Wang, B. Mechanical properties and thermal conductivity of graphene reinforced copper matrix composites. Powder Technol. 2016, 301, 601–607. [Google Scholar] [CrossRef]
  18. El-Kady, O.; Yehia, H.M.; Nouh, F. Preparation and characterization of Cu/(WC-TiC-Co)/graphene nano-composites as a suitable material for heat sink by powder metallurgy method. Int. J. Refract. Met. Hard Mater. 2019, 79, 108–114. [Google Scholar] [CrossRef]
  19. Köck, T.; Brendel, A.; Bolt, H. Interface reactions between silicon carbide and interlayers in silicon carbide-copper metal-matrix composites. J. Nucl. Mater. 2007, 362, 197–201. [Google Scholar] [CrossRef]
  20. Chen, L.; Huang, Z.; Kumar, S. Impact of bonding at multi-layer graphene/metal Interfaces on thermal boundary conductance. RSC Adv. 2014, 4, 35852–35861. [Google Scholar] [CrossRef]
  21. Liu, B.; Khvesyuk, V.I.; Barinov, A.A.; Wang, M. Effect of interfacial roughness on thermal boundary conductance: An elastic wave model using the Kirchhoff approximation. Int. J. Mech. Sci. 2022, 218, 106993. [Google Scholar] [CrossRef]
  22. Hong, Y.; Li, L.; Zeng, X.C.; Zhang, J.C. Tuning thermal contact conductance at graphene–copper interface via surface nanoengineering. Nanoscale 2015, 7, 6286–6294. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Wei, N.; Zhou, C.; Li, Z.; Ou, B.; Zhao, K.; Yu, P.; Li, S.; Zhao, J. Thermal conductivity of Aluminum/Graphene metal-matrix composites: From the thermal boundary conductance to thermal regulation. Mater. Today Commun. 2022, 30, 103147. [Google Scholar] [CrossRef]
  24. Ou, B.; Yan, J.; Wang, Q.; Lu, L. Thermal Conductance of Graphene-Titanium Interface: A Molecular Simulation. Molecules 2022, 27, 905. [Google Scholar] [CrossRef]
  25. Ghosh, S.; Bao, W.; Nika, D.L.; Subrina, S.; Pokatilov, E.P.; Lau, C.N.; Balandin, A.A. Dimensional crossover of thermal transport in few-layer graphene. Nat. Mater. 2010, 9, 555–558. [Google Scholar] [CrossRef]
  26. Dumont, D.; Bougeard, D. Simulation of vibrational spectra of polydichlorophosphazene by molecular dynamics calculations: A study of the relaxed and stretched polymer. Comput. Theor. Polym. Sci. 1999, 9, 89–97. [Google Scholar] [CrossRef]
  27. Fukui, K.; Sumpter, B.G.; Barnes, M.D.; Noid, D.W. Molecular dynamics studies of the structure and properties of polymer nano-particles. Comput. Theor. Polym. Sci. 1999, 9, 245–254. [Google Scholar] [CrossRef]
  28. Plimpton, S. Fast Parallel Algorithms for Short-Range Molecular Dynamics. J. Comput. Phys. 1995, 117, 1–19. [Google Scholar] [CrossRef] [Green Version]
  29. Daw, M.S.; Baskes, M.I. Embedded-atom method: Derivation and application to impurities, surfaces, and other defects in metals. Phys. Rev. B 1984, 29, 6443–6453. [Google Scholar] [CrossRef] [Green Version]
  30. Fischer, F.; Schmitz, G.; Eich, S.M. A systematic study of grain boundary segregation and grain boundary formation energy using a new copper–nickel embedded-atom potential. Acta Mater. 2019, 176, 220–231. [Google Scholar] [CrossRef]
  31. Brenner, D.W.; Shenderova, O.A.; Harrison, J.A.; Stuart, S.J.; Ni, B.; Sinnott, S.B. A second-generation reactive empirical bond order (REBO) potential energy expression for hydrocarbons. J. Phys. Condens. Mat. 2002, 14, 783–802. [Google Scholar] [CrossRef] [Green Version]
  32. Dhaliwal, G.; Nair, P.B.; Singh, C.V. Uncertainty analysis and estimation of robust AIREBO parameters for graphene. Carbon 2019, 142, 300–310. [Google Scholar] [CrossRef]
  33. Long, X.J.; Li, B.; Wang, L.; Huang, J.Y.; Zhu, J.; Luo, S.N. Shock response of Cu/graphene nanolayered composites. Carbon 2016, 103, 457–463. [Google Scholar] [CrossRef]
  34. Weng, S.; Ning, H.; Fu, T.; Hu, N.; Zhao, Y.; Huang, C.; Peng, X. Molecular dynamics study of strengthening mechanism of nanolaminated graphene/Cu composites under compression. Sci. Rep. 2018, 8, 3089. [Google Scholar] [CrossRef] [Green Version]
  35. Cao, B.Y.; Zou, J.H.; Hu, G.J.; Cao, G.X. Enhanced thermal transport across multilayer graphene and water by interlayer functionalization. Appl. Phys. Lett. 2018, 112, 041603. [Google Scholar] [CrossRef] [Green Version]
  36. Maiti, A.; Mahan, G.D.; Pantelides, S.T. Dynamical simulations of nonequilibrium processes—Heat flow and the Kapitza resistance across grain boundaries. Solid State Commun. 1997, 102, 517–521. [Google Scholar] [CrossRef]
  37. Li, Q.; Yu, Y.; Liu, Y.; Liu, C.; Lin, L. Thermal Properties of the Mixed n-Octadecane/Cu Nanoparticle Nanofluids during Phase Transition: A Molecular Dynamics Study. Materials 2017, 10, 38. [Google Scholar] [CrossRef] [Green Version]
  38. Schneider, T.; Stoll, E. Molecular-dynamics study of a three-dimensional one-component model for distortive phase transitions. Phys. Rev. B 1978, 17, 1302–1322. [Google Scholar] [CrossRef]
  39. Mollow, B.R. Power Spectrum of Light Scattered by Two-Level Systems. Phys. Rev. 1969, 188, 1969–1975. [Google Scholar] [CrossRef]
  40. Liu, Y.; Hu, C.; Huang, J.; Sumpter, B.G.; Qiao, R. Tuning interfacial thermal conductance of graphene embedded in soft materials by vacancy defects. J. Chem. Phys. 2015, 142, 244703. [Google Scholar] [CrossRef]
Figure 1. (a) Schematic of the NEMD method. (b) Temperature distribution of the copper/graphene model. The temperature differences near the hot and cold zones are denoted by Δ Thot and Δ Tcold, respectively. (c) Atomic structures of copper (100), (110), and (111) surfaces.
Figure 1. (a) Schematic of the NEMD method. (b) Temperature distribution of the copper/graphene model. The temperature differences near the hot and cold zones are denoted by Δ Thot and Δ Tcold, respectively. (c) Atomic structures of copper (100), (110), and (111) surfaces.
Materials 15 07588 g001
Figure 2. Effect of graphene layer number on TBC.
Figure 2. Effect of graphene layer number on TBC.
Materials 15 07588 g002
Figure 3. (a) Vibrational spectra of different modes of Cu and graphene atoms at the interface; n depicts the number of graphene layers. (b) The overlap concentration of phonon spectrum S with different layers; Sout and Sin represent the corresponding values of in-plane and out-of-plane modes, respectively.
Figure 3. (a) Vibrational spectra of different modes of Cu and graphene atoms at the interface; n depicts the number of graphene layers. (b) The overlap concentration of phonon spectrum S with different layers; Sout and Sin represent the corresponding values of in-plane and out-of-plane modes, respectively.
Materials 15 07588 g003
Figure 4. (a) Phonon spectra of single graphene at 300 K. (b) Out-of-plane and in-plane vibrational spectra of different surfaces. (c) The overlap parameter S of the VDOS and corresponding TBC with distinct surfaces.
Figure 4. (a) Phonon spectra of single graphene at 300 K. (b) Out-of-plane and in-plane vibrational spectra of different surfaces. (c) The overlap parameter S of the VDOS and corresponding TBC with distinct surfaces.
Materials 15 07588 g004
Figure 5. (ac) Potential energy distribution contour of adjacent graphene layers with different atomic structures. Δ E represents the difference value between the maximum and minimum potential energy.
Figure 5. (ac) Potential energy distribution contour of adjacent graphene layers with different atomic structures. Δ E represents the difference value between the maximum and minimum potential energy.
Materials 15 07588 g005
Figure 6. Effect of Cu-bulk thickness on TBC.
Figure 6. Effect of Cu-bulk thickness on TBC.
Materials 15 07588 g006
Figure 7. (a) Variation in Z-direction interfacial thermal conductance of graphene as a function of graphene defect density. (b) Out-of-plane and in-plane VDOS values of Cu and graphene atoms at the interface with diverse density of vacancy defect. (c) The overlap of phonon spectrum S with different density of vacancy defect; Sout and Sin represent the corresponding values of in-plane and out-of-plane modes, respectively.
Figure 7. (a) Variation in Z-direction interfacial thermal conductance of graphene as a function of graphene defect density. (b) Out-of-plane and in-plane VDOS values of Cu and graphene atoms at the interface with diverse density of vacancy defect. (c) The overlap of phonon spectrum S with different density of vacancy defect; Sout and Sin represent the corresponding values of in-plane and out-of-plane modes, respectively.
Materials 15 07588 g007
Table 1. The studies of different metal matrix/graphene composites using NEMD method.
Table 1. The studies of different metal matrix/graphene composites using NEMD method.
CompositeMethodLayer NumberAtomic StructureVacancy RateOtherRef.
Cu/GrNEMD(011) > (111) > (001)Length of copper matrixThis work
Al/GrNEMD(110) > (111) > (100)×Strain[23]
Ti/GrNEMD(0001) > ( 11 2 ¯ 0 ) > ( 10 1 ¯ 0 )×Strain[24]
Ni/GrNEMD×××[25]
Table 2. The model parameters choosing different crystal planes of Cu and the corresponding mismatch strain of graphene.
Table 2. The model parameters choosing different crystal planes of Cu and the corresponding mismatch strain of graphene.
Crystal PlaneLx (nm) ε Gr−xLy (nm) ε Gr−yLt (nm)
Cu(001)4.690.407%4.690.085%5.42
Cu(011)2.430.57%4.52−3.69%4.98
Cu(111)2.430.57%5.49−0.89%6.04
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zhu, J.; Huang, S.; Xie, Z.; Guo, H.; Yang, H. Thermal Conductance of Copper–Graphene Interface: A Molecular Simulation. Materials 2022, 15, 7588. https://doi.org/10.3390/ma15217588

AMA Style

Zhu J, Huang S, Xie Z, Guo H, Yang H. Thermal Conductance of Copper–Graphene Interface: A Molecular Simulation. Materials. 2022; 15(21):7588. https://doi.org/10.3390/ma15217588

Chicago/Turabian Style

Zhu, Jiarui, Shuhui Huang, Zhongnan Xie, Hong Guo, and Hui Yang. 2022. "Thermal Conductance of Copper–Graphene Interface: A Molecular Simulation" Materials 15, no. 21: 7588. https://doi.org/10.3390/ma15217588

APA Style

Zhu, J., Huang, S., Xie, Z., Guo, H., & Yang, H. (2022). Thermal Conductance of Copper–Graphene Interface: A Molecular Simulation. Materials, 15(21), 7588. https://doi.org/10.3390/ma15217588

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