Luiz Felipe C. Pereira

Modelling heat conduction in polycrystalline hexagonal boron-nitride films

We conducted extensive molecular dynamics simulations to investigate the thermal conductivity of polycrystalline hexagonal boron-nitride (h-BN) films. To this aim, we constructed large atomistic models of polycrystalline h-BN sheets with random and uniform grain configuration. By performing equilibrium molecular dynamics (EMD) simulations, we investigated the influence of the average grain size on the thermal conductivity of polycrystalline h-BN films at various temperatures. Using the EMD results, we constructed finite element models of polycrystalline h-BN sheets to probe the thermal conductivity of samples with larger grain sizes. Our multiscale investigations not only provide a general viewpoint regarding the heat conduction in h-BN films but also propose that polycrystalline h-BN sheets present high thermal conductivity comparable to monocrystalline sheets.

Anisotropic thermal conductivity and mechanical properties of phagraphene: a molecular dynamics study

Phagraphene is a novel 2D carbon allotrope with interesting electronic properties which has been recently theoretically proposed. Phagraphene is similar to a defective graphene structure with an arrangement of pentagonal, heptagonal and hexagonal rings. In this study we investigate the thermal conductivity and mechanical properties of phagraphene using molecular dynamics simulations. Using the non-equilibrium molecular dynamics method, we found the thermal conductivity of phagraphene to be anisotropic, with room temperature values of 218 ± 20 W m−1 K−1 along the armchair direction and 285 ± 29 W m−1 K−1 along the zigzag direction. Both values are one order of magnitude smaller than for pristine graphene. The analysis of the phonon group velocities also shows a significant reduction in this quantity for phagraphene in comparison to graphene. By performing uniaxial tensile simulations, we studied the deformation process and mechanical response of phagraphene. We found that phagraphene exhibits a remarkable high tensile strength around 85 ± 2 GPa, whereas its elastic modulus is also anisotropic along the in-plane directions, with values of 870 ± 15 GPa and 800 ± 14 GPa for the armchair and zigzag directions, respectively. The lower thermal conductivity of phagraphene along with its predicted electronic properties suggests that it could be a better candidate than graphene in future carbon-based thermoelectric devices.

Thermal conductivity decomposition in two-dimensional materials: Application to graphene

Two-dimensional materials have unusual phonon spectra due to the presence of flexural (out-of-plane) modes. Although molecular dynamics simulations have been extensively used to study heat transport in such materials, conventional formalisms treat the phonon dynamics isotropically. Here, we decompose the microscopic heat current in atomistic simulations into in-plane and out-of-plane components, corresponding to in-plane and out-of-plane phonon dynamics, respectively. This decomposition allows for direct computation of the corresponding thermal conductivity components in two-dimensional materials. We apply this decomposition to study heat transport in suspended graphene, using both equilibrium and nonequilibrium molecular dynamics simulations. We show that the flexural component is responsible for about two-thirds of the total thermal conductivity in unstrained graphene, and the acoustic flexural component is responsible for the logarithmic divergence of the conductivity when a sufficiently large tensile strain is applied.

Bimodal grain-size scaling of thermal transport in polycrystalline graphene from large-scale molecular dynamics simulations

Grain boundaries in graphene are inherent in wafer-scale samples prepared by chemical vapor deposition. They can strongly influence the mechanical properties and electronic and heat transport in graphene. In this work, we employ extensive molecular dynamics simulations to study thermal transport in large suspended polycrystalline graphene samples. Samples of different controlled grain sizes are prepared by a recently developed efficient multiscale approach based on the phase field crystal model. In contrast to previous works, our results show that the scaling of the thermal conductivity with the grain size implies bimodal behavior with two effective Kapitza lengths. The scaling is dominated by the out-of-plane (flexural) phonons with a Kapitza length that is an order of magnitude larger than that of the in-plane phonons. We also show that, to get quantitative agreement with the most recent experiments, quantum corrections need to be applied to both the Kapitza conductance of grain boundaries and the thermal conductivity of pristine graphene, and the corresponding Kapitza lengths must be renormalized accordingly.

Controlling resonant tunneling in graphene via Fermi velocity engineering

We investigate the resonant tunneling in a single layer graphene superlattice with modulated energy gap and Fermi velocity via an effective Dirac-like Hamiltonian. We calculate the transmission coefficient with the transfer matrix method and analyze the effect of a Fermi velocity modulation on the electronic transmission, in the case of normal and oblique incidence. We find it is possible to manipulate the electronic transmission in graphene by Fermi velocity engineering, and show that it is possible to tune the transmitivity from 0 to 1. We also analyze how a Fermi velocity modulation influences the total conductance and the Fano factor. Our results are relevant for the development of novel graphene-based electronic devices.

Thermal Conductivity of Graphene-hBN Superlattice Ribbons

Superlattices are ideal model systems for the realization and understanding of coherent (wave-like) and incoherent (particle-like) phonon thermal transport. Single layer heterostructures of graphene and hexagonal boron nitride have been produced recently with sharp edges and controlled domain sizes. In this study we employ nonequilibrium molecular dynamics simulations to investigate the thermal conductivity of superlattice nanoribbons with equal-sized domains of graphene and hexagonal boron nitride. We analyze the dependence of the conductivity with the domain sizes, and with the total length of the ribbons. We determine that the thermal conductivity reaches a minimum value of 89 W m−1K−1 for ribbons with a superlattice period of 3.43 nm. The effective phonon mean free path is also determined and shows a minimum value of 32 nm for the same superlattice period. Our results also reveal that a crossover from coherent to incoherent phonon transport is present at room temperature for BNC nanoribbons, as the superlattice period becomes comparable to the phonon coherence length. Analyzing phonon populations relative to the smallest superlattice period, we attribute the minimum thermal conductivity to a reduction in the population of flexural phonons when the superlattice period equals 3.43 nm. The ability to manipulate thermal conductivity using superlattice-based two-dimensional materials, such as graphene-hBN nanoribbons, opens up opportunities for application in future nanostructured thermoelectric devices.

BxCyNz hybrid graphenylene: stability and electronic properties

Interest in hybrid monolayers with arrangements that differ from that of the honeycomb lattice has been growing. However, systematic investigations on the properties of these structures are still lacking. In this work, we combined density functional theory (DFT) and molecular dynamics (MD) simulations to study the stability and electronic properties of nanosheets composed of B, C, and N atoms arranged in the pattern of the carbon allotrope graphenylene. We considered twenty structures with varied atomic arrangements and stoichiometries, which we call BxCyNz hybrid graphenylenes. We calculated the formation energy for each arrangement, and found that it decreases as the number of B–C and N–C bonds decreases. We also found that the structure with minimum energy has stoichiometry B2CN and an atomic arrangement with BN and C stripes connected along the zigzag direction. Regarding the electronic properties, we found that all investigated structures are semiconductors, with band gaps ranging from 0.14 to 1.65 eV. Finally, we found that the optimized hybrid lattices presented pores of varied sizes and shapes. This diversity in pore geometry suggests that these structures might be particularly suited for molecular sieve applications.

Tuning the Fano factor of graphene via Fermi velocity modulation

Abstract In this work we investigate the influence of a Fermi velocity modulation on the Fano factor of periodic and quasi-periodic graphene superlattices. We consider the continuum model and use the transfer matrix method to solve the Dirac-like equation for graphene where the electrostatic potential, energy gap and Fermi velocity are piecewise constant functions of the position x . We found that in the presence of an energy gap, it is possible to tune the energy of the Fano factor peak and consequently the location of the Dirac point, by a modulation in the Fermi velocity. Hence, the peak of the Fano factor can be used experimentally to identify the Dirac point. We show that for higher values of the Fermi velocity the Fano factor goes below 1/3 at the Dirac point. Furthermore, we show that in periodic superlattices the location of Fano factor peaks is symmetric when the Fermi velocity v A and v B is exchanged, however by introducing quasi-periodicity the symmetry is lost. The Fano factor usually holds a universal value for a specific transport regime, which reveals that the possibility of controlling it in graphene is a notable result.