Electronic Properties of Homo- and Heterobilayer Graphyne: The Idea

Nov 2, 2016 - Our predicted energy- and charge-storing capacities are greater than those predicted for nanocapacitors based on graphene and hexagonal ...
0 downloads 0 Views 3MB Size
Subscriber access provided by UNIV OF WATERLOO

Article

Electronic Properties of Homo and Hetero Bilayer Graphyne: The Idea of a Nanocapacitor Barnali Bhattacharya, Utpal Sarkar, and Nicola Seriani J. Phys. Chem. C, Just Accepted Manuscript • DOI: 10.1021/acs.jpcc.6b07092 • Publication Date (Web): 02 Nov 2016 Downloaded from http://pubs.acs.org on November 8, 2016

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

The Journal of Physical Chemistry C is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 40

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Electronic Properties of Homo and Hetero Bilayer Graphyne: The Idea of a Nanocapacitor Barnali Bhattacharya1, Utpal Sarkar 1* and Nicola Seriani 2* 1

Department of Physics, Assam University, Silchar-788011, INDIA

2

The Abdus Salam ICTP, Trieste, Italy

*

To whom correspondence should be addressed. E-mail [email protected] (Utpal

Sarkar), [email protected] (Nicola Seriani)

ACS Paragon Plus Environment

1

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 2 of 40

ABSTRACT

We have investigated the capacitive behavior of bilayer graphyne and its boron nitride derivatives by first-principles simulations based on density functional theory, including van der Waals interactions. Our predicted energy and charge storing capacity are greater than those predicted for nanocapacitors based on graphene and hexagonal boron nitride. In the most stable configuration, the two layers are stacked on top of each other, just as in bulk graphyne. The stacking arrangement has a strong effect on the electronic properties of the system: the stable stacking configurations for the graphyne systems are semiconductors with direct band gaps of 0.38 eV and 0.50 eV respectively. Substitutional boron–nitrogen doping provides a way to tune the band gap of the system. The band gap generally increases in presence of the dopants, but the value of the band gap depends on the substitution sites. This suggests that controlled boron nitride doping of graphyne could be a useful and flexible method to build nanoscale electronic and optoelectronic devices.

ACS Paragon Plus Environment

2

Page 3 of 40

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

INTRODUCTION Electronics is facing great challenges as the size of typical devices shrink down to few nanometers in size, at a level where the precise atomic structure of the device and quantum effects play a fundamental role. Therefore the search for novel materials and new approaches to nanoscale electronics is being run extensively. In this context, low-dimensional carbon-based layered materials attracted considerable attention. 1-4 Due to the flexibility of carbon in forming allotropes and the wide range of possible applications, tremendous efforts have been devoted to the discovery and study of new possible carbon allotropes during the past few decades.2-4 Newly synthesized important members of the carbon family include zero–dimensional fullerene,2 quasi– one–dimensional carbon nanotubes,3 two–dimensional graphene.4 In particular, bilayer and fewlayer graphene have already been studied at large scale, including assemblies of graphene with hexagonal boron nitride (h-BN)5-7 or with SiC8 or with Phosphorene.9 These systems exhibit some outstanding electronic, optical and mechanical properties, including capacitive behavior. 1014

It is assumed that these carbon allotropes and carbon-based molecules with their tunable

properties can potentially become building blocks in next-generation electronic and optoelectronic devices, thereby increasing the interest in new carbon allotropes. Baughman et al. in 198715 proposed one such 2D carbon allotrope built from graphene where acetylenic linkages (–C≡C–) are present in the unit cell, instead of the double bond (C=C) as in graphene. Depending on the various arrangements of acetylenic linkages and percentage of the inserted acetylenic linkages (–C≡C–), graphyne is mainly divided into three categories– α graphyne, β graphyne and γ graphyne. The α graphyne is the most symmetrical one and resembles graphene. It can be obtained from graphene by replacing every C-C bond of graphene with acetylenic linkages (–C≡C–), thus 100% acetylenic groups are present in the α graphyne. Whereas in β

ACS Paragon Plus Environment

3

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 4 of 40

graphyne, acetylenic groups are inserted into two-thirds of C-C bonds of graphene and in γ graphyne, only one-third of C-C bonds of graphene contain acetylenic groups. Due to the structural similarity, the band structure of α graphyne is very similar to graphene, containing the same type of Dirac cone, only the pseudospin in α graphyne has the opposite direction to that in graphene.16 The presence of Dirac cones are also observed in β graphyne, but it is not like graphene.16 A previous study16-22 by some of the authors shows that monolayer γ-graphyne is a direct band gap semiconductor and its band gap is tunable under different circumstances. Due to the tunable band gap property graphyne can be used in building p-n junction. In addition, graphyne has high stability in comparison to graphene. Since gamma graphyne is the lowest energy member among the three types of graphyne and semiconductor with a small band gap, we focus our study on γ- graphyne. Graphyne induce a rich variety of optical, electronic and elastic properties17-20,23-24 due to presence of triple bonds. Recently, substructures of γ- graphyne have been fabricated as low-dimensional nanostructures.25-29 Moreover, another family member of the γ-graphyne family, with di-acetylenic linkages, graphydine, has been synthesized in the form of large-area-multilayer films,30 flakes,31 nanotubes, 32 and nanowires.33 Much less attention has been devoted to the graphyne bilayer. In fact, theoretical investigations exist only for bilayers of graphydine,34 for α - and β -graphyne,35-37 showing the possibility to tune the band gap by doping with boron nitride and on external electric field. Moreover, a subtle dependence of the gap on position and concentration of the dopants is also discussed. The band gap modulation of α and β graphyne due to applied strain and electric field35-38 highlights the possibility of using them in electronic and electro-mechanical devices. Till today all the reported study of bilayer graphyne (α graphyne, β graphyne and graphdiyne) explored the electronic structure and properties under different stacking modes,37 external electric field,35 strain38 and

ACS Paragon Plus Environment

4

Page 5 of 40

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

substitutional doping.36 The electronic properties of three dimensional forms of graphyne39 and lithium decorated multilayer graphyne as a Li ion battery anode40 are reported also. It seems that bilayer systems offer a wide range of possibilities to control the gap, including acting on the stacking of the bilayer. But, to the best of our knowledge, the band gap tailoring of γ- graphyne due to combined effect of systematic BN doping and external electric field has not been explored so far. Moreover, the application of bilayer graphyne as a charge storage medium and its possible use as nanocapacitor is still unrevealed. Through this study, in particular, we are interested in the possibility to tune γ-graphyne bilayer properties in order to employ it as a nanocapacitor. To nullify the major drawbacks of the conventional capacitor such as comparatively heavy weight, slow response time due to the motion of ions in electrochemical reactions and to meet the up growing demand for lightweight, recyclable, efficient, and high-capacity energy storage capacitor the researcher are testing the capacitive behavior of nanoscale material and designed nanoscale nanocapacitor. Nanoscale capacitors have been recently shown to display very interesting functional behavior, such as faster charge storing and releasing capacity, higher amounts of charge delivering capability at higher power rates, longer life with short load,41-46 that might be relevant for practical applications. Other carbon allotropes have already been investigated for this purpose. Carbon nanotube based capacitors44,45,47 was proven to be highly efficient for energy storage in connection with photovoltaic devices. Graphene has been proposed as ideal nanocapacitor material48-51 and has recently shown high capacitance per unit area. Recent calculations of periodic hetero-layer systems seem to confirm this potential, especially for graphene–silicene bilayers13 and graphene-h-BN bilayers.10-12 In this paper, we have investigated graphyne-based bilayers by first-principles simulations to assess their properties in view of a possible employment as nanocapacitors. We have considered

ACS Paragon Plus Environment

5

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 6 of 40

the effects of different stacking modes and of doping through boron nitride. Our approach is based on the charge redistribution at the interface of bilayer graphyne because of the externally or intrinsically induced interfacial electric field that mimics the mechanism of a parallel plate capacitor. In the next section, the details of the employed computational techniques are reported. First, the results on the relation between stacking and electronic structure are reported; then, the effects of BN doping are presented, and finally the capacitive behavior of this system is investigated. A summary concludes the article. COMPUTATIONAL DETAILS We have done density functional theory (DFT) calculations with Perdew-Zunger (LDA) functional to describe exchange and correlation. Van der Waals interactions were included by using the second version of the van der Waals density functional (vdW-DF2) theory proposed by Lee et al.52 All atomic coordinates are fully relaxed using the conjugate gradient algorithm, until the forces on each atom are equal to or less than 1× 10−9 Ry/au and the energy difference between the consecutive scf steps is less than 1 × 10−9 Ry. 10 k-points have been chosen along the periodic direction of unit cell and 4 k-points are chosen along the vacuum direction. The vacuum was set to 15 Å to avoid the interactions between periodically repeated images of the bilayer system. Structure optimizations have been carried out with the Quantum-ESPRESSO 5.1.2 package,53 which employs a plane-wave basis set to describe wavefunctions. A cutoff of 65 Ry and 520 Ry have been employed, respectively for the wavefunctions and for the charge density. The smearing scheme by Marzari and Vanderbilt has been used for orbital occupations.54-55 The scalar relativistic ultrasoft pseudopotential [version ‘2.0.1’] with non-linear core correction has been employed in the calculation. Capacitor properties such as dipole

ACS Paragon Plus Environment

6

Page 7 of 40

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

moment (D), excess charge (Q), energy stored by the capacitor (E c) and capacitance (C) have been calculated using the SIESTA 3.2 code.56 Here, a further self-consistent calculation has been done for the optimized geometry with an electric field. SIESTA code works with atomic orbital (AO) basis sets while in Quantum-ESPRESSO plane wave basis set is used. For the same exchange-correlation functional, the unoccupied states in the band diagram can be reproduced better using well-converged plane-wave (PW) basis set than the AO basis sets.57 Thus, the plane wave basis set is a better choice for optimizing and electronic properties calculation of periodic system, and so we have carried out optimization and electronic properties calculation of bilayer system in Quantum-ESPRESSO suit. For capacitive properties calculation, we have used SIESTA code since it directly provides the electric dipole moment for a finite electric field. The interplanar binding energy or adhesion energy of a bilayer has been calculated as (1) (2) where

and

are the total energy of the bilayer system and of the monolayer

respectively. The formation energy of a BN doped hetero-bilayer has been calculated as-

(3)

=

where

and

are the total energies of BN doped bilayer graphyne and the perfect bilayer

graphyne. The number m is the number of CC pairs that are substituted by the BN pairs in the

ACS Paragon Plus Environment

7

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 8 of 40

unit cell. μC–C and μB–N are the energies per unit cell of graphene and h-BN monolayer respectively. RESULTS AND DISCUSSION First, the stacking of the bilayer has been investigated. There are three main types of stacking arrangements, shown in Figure 1; two of them, namely AA and AB, resemble those of graphene. In the AA stacking, the two layers are orientationally identical and lie exactly on top of each other. In the AB stacking the two layers are orientationally mismatched. In addition, one more stacking mode is expected in bilayer graphyne, namely A⍺Aβ, where the hexagonal hollow of one layer is stacked on the triangular void of the other layer.

Figure 1. Optimized geometries with unit cell of: (a) AA stacked graphyne; (b) AB stacked graphyne; (c) A⍺Aβ stacked graphyne Structural parameters and binding energies of bilayer graphynes are presented in Table 1 and Table 2. ‘Table 1’ represents pristine bilayer graphyne while ‘Table 2’ represents BN doped bilayer graphyne. All bilayers are more stable than the monolayer, with the A⍺Aβ structure having the largest binding energy. This is the same stacking found to be stable in bulk graphyne

ACS Paragon Plus Environment

8

Page 9 of 40

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

and bulk graphdiyne.39,58 For all systems the inclusion of van der Waals corrections is crucial for a correct description of the interlayer interactions (see the Supporting Information Table S1 for a comparison). In fact, the interlayer distance increases when vdW interactions are included. This counterintuitive result is in agreement with what found in other layered systems.59-62 Table 1. Calculated lattice constant (a); interlayer distance (d), binding energy (BE), and energy band gap with van der Waals (vdW) correction of bilayer graphyne with different stacking modes.

Stacking

Lattice parameter a Interlayer distance d (Å) Binding energy (eV)

Band gap (eV)

(Å) AA

6.900

3.720

-0.463

0.120

AB

6.900

3.520

-0.558

0.510

A⍺Aβ

6.900

3.510

-0.571

0.380

Figure 2. Band structure and projected density of states (PDOS) of the A⍺Aβ stacked graphyne

ACS Paragon Plus Environment

9

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 10 of 40

The electronic band diagram of bilayer graphyne is very similar to that of monolayer graphyne.16-20 The valence band maximum (VBM) and conduction band minimum (CBM) are located at the M point of the Brillouin zone, making the bilayer a small-gap semiconductor. The band structure of the stable A⍺Aβ configuration is shown in Figure 2 (a), and displays a gap of 0.38 eV. In the other configurations, the band structure is very similar (see Supporting Information Figure S1). Each band of the monolayer is split in two bands differing by a small energy, due to the interlayer interaction. The projected density of states of bilayer graphyne is depicted in Figure 2 (b). The contributions of each layer to the total DOS is exactly equal at any energy. The states near the gap, including HOMO and LUMO, are mainly due to the pz orbitals of carbon irrespective of the stacking arrangement (see Supporting Information Figure S2). The density of states (DOS) and projected density of states (PDOS) of monolayer graphyne [Supporting Information Figure S3] and most stable bilayer graphyne show that the contribution of px+py is comparatively larger in a certain energy range of VB and CB for bilayer graphyne compare to monolayer graphyne. Although near the Fermi level (depicted as dotted line in Figure S3) only the pz orbital contributed to DOS for both mono layer and bilayer graphyne and the noticeable contribution of px+py starts at -1.710 eV at VB and 3.202 eV at CB. For monolayer graphyne, throughout the considered energy range, contribution of pz orbital is larger than px+py orbital. Only in the energy range -2.262 to -2.720 eV and 4.500 to 4.700 eV, the contribution of px+py orbital is comparable with pz orbital. But for bilayer graphyne, in the energy range -1.900 to -3.100 eV and 3.400 to 4.700 eV the contribution of px+py orbital is greater than pz orbital. BN doping:

ACS Paragon Plus Environment

10

Page 11 of 40

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Co-doping by boron and nitrogen is a common method to engineer the band gap of carbon systems. Since A⍺Aβ is the stable structure of the pure bilayer, we have investigated BN doping in this structure only. We have considered nine main configurations for the doped bilayer. In three of them, there is one BN unit per layer and the two layers are structurally equivalent (homo-layer structures), while in the other six the two layers are not equivalent (hetero-layer structures). The most stable homo- and hetero-layer structures are presented in Figure 3, while the others are shown in the Supporting Information Figure S4. The structures differ from one another for the position of the BN units, which can be at the chain sites or at the ring sites. Moreover, we have considered the possibility of the full substitution of carbon with BN in one layer, leading to the formation of one layer of h-BN.

Figure 3. Optimized geometry of: (a) most stable hetero layer “pristine + graphyne like BN sheet” (b) most stable homo layer “graphyne like BN sheet + graphyne like BN sheet”

ACS Paragon Plus Environment

11

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 12 of 40

BN doping induces structural changes in the bilayer, in particular in the interlayer distance and lattice constant. BN doping leads to a decrease in the interlayer distance when BN is present in both layers while BN doping of only one layer leads to an increase in the interlayer distance. Similar to bilayer and trilayer graphdiyne

34

and nitrogen doped graphene bilayer63 interlayer

mixing of п orbitals is expected in bilayer graphyne. Our crystal orbital overlap population64 (COOP) calculation shows the existence of covalent-like п-п interaction63 between the layers. It has been predicted from the calculated integrated crystal orbital overlap population (INTCOOP) [supporting information Table S2] that, most of the п-п interactions are bonding nature which plays an important role in determining the interlayer distance.63 In addition to п-п bonding, electrostatic, dispersive and Pauli interactions are also crucial for anchoring the layers at equilibrium distance. The delicate interplay between electrostatic interaction, van der Waals interaction, Pauli repulsion and covalent-like п-п interaction between the layers, actually dictates the equilibrium interlayer distance. One of the possible explanations of increased interlayer distance of the ‘bilayer containing BN at only one layer’ may be the reduction of п-п bonding interaction between the layers compare to pristine bilayer [see Table S2 in the supporting information]. The п-п bonding interaction is negligible for the structure ‘pristine + graphyne like BN sheet’ in which interlayer separation is highest and interlayer binding energy is lowest. The significant reduction of interlayer distance for structures ‘BN at chain + BN at chain’; ‘BN at chain + BN at ring’; ‘BN at ring + BN at ring’; may be lies on the enhancement of п-п bonding interaction between the layers compare to the pristine bilayer. But, the reduction of the interlayer distance for structures ‘graphyne like BN sheet + graphyne like BN sheet’; ‘BN sheet + BN at linear chain site’; ‘BN sheet + BN at linear ring site’ cannot be explained on the basis of п-п bonding interaction only. For the bilayer structures having B-N bond in both layers, due to the

ACS Paragon Plus Environment

12

Page 13 of 40

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

polarity of B-N bond, effective charges are developed in the atomic centres, as a result of which electrostatic interaction develop. Although electrostatic interaction exists but previous studies65-67 show that interlayer binding hardly depends on electrostatic interaction; the van der Waals interaction is mainly responsible for fixing the interlayer distance apart from п-п interaction.62 However, the proper reason of variation of interlayer distance with BN doping at various positions of this type of polar bilayer material is a subject of ongoing research. Table 2. Calculated lattice constant (a), interlayer distance (d), binding energy (BE) and band gap (Eg) of homo layer and hetero layer graphyne derivatives. Systems

BN at linear chain site

Lattice parameter a Interlayer distance d

Binding Energy BE Band gap (eV)

(Å)

(eV)

(Å) 6.970

3.490

-0.558

1.440

7.000

3.460

-0.612

2.660

7.000

3.490

-0.531

3.950

6.930

3.520

-0.517

0.550

+ BN at linear chain site BN at ring site + BN at ring site Graphyne

like

BN

sheet + graphyne like BN sheet Pristine graphyne + BN at the linear chain site

ACS Paragon Plus Environment

13

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Pristine graphyne +

Page 14 of 40

6.950

3.520

-0.503

0.600

6.950

3.530

-0.463

0.590

6.980

3.490

-0.558

1.560

6.980

3.490

-0.531

1.560

7.000

3.470

-0.544

2.800

BN at ring site

Pristine graphyne + graphyne

like

BN

sheet BN at linear chain site + BN at ring site BN sheet + BN at linear chain site

BN sheet + BN at ring site

The in-plane lattice constant is calculated by plotting the total energy as a function of lattice spacing. The total energy with full relaxation of all the atoms, as a function of lattice spacing, is obtained by varying the lattice constant a. We have considered seven equispaced lattice constants. The value of ‘a’ at which the total energy is minimum for this graph is identified as inplane lattice constant. The in-plane lattice constant is minimal for the pristine bilayer and increases through BN doping. We have calculated the binding energy of the bilayer with respect to the monolayer (Table 2). With our convention, a negative value of the binding energy stands for an exothermic binding. In fact, bilayer formation is always exothermic. On the other side, the

ACS Paragon Plus Environment

14

Page 15 of 40

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

BN doping process is endothermic except two cases having negative formation energies (see Table 3). Among the homo-bilayer structures, the bilayer formed by graphyne-like BN sheets is energetically most favorable as it has the lowest formation energy (-2.068 eV) indicating that it is relatively easy to substitute BN throughout the structure of graphyne bilayer. On the other hand, among the hetero layer structures, bilayer originate due to coupling of pristine graphyne and graphyne like BN sheet is most stable as it has negative formation energy. Table 3. Formation energies of homo and hetero bilayer graphyne derivatives.

Systems

Formation energy (eV) Homo layers

BN at linear chain site + BN at linear chain site

3.497

BN at ring site + BN at ring site

3.687

Graphyne like BN sheet + graphyne like BN sheet

-2.068

Hetero layers Pristine graphyne + BN at the linear chain site

1.796

Pristine graphyne + BN at ring site

1.932

Pristine graphyne + graphyne like BN sheet

-0.939

BN at linear chain site + BN at ring site

3.619

BN sheet + BN at linear chain site

0.721

ACS Paragon Plus Environment

15

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

BN sheet + BN at ring site

Page 16 of 40

0.830

The effect of BN doping on electronic structure: The electronic band structures of selected BN-derivatives of bilayer graphyne are shown in Figure 4 and Supporting Information Figure S5. Doping increases the band gap, and the graphyne-like BN bilayer has the largest gap of 3.95 eV. Still, all systems have a direct band gap. In systems with BN doping, the band gap depends on the position of the BN pair.

Figure 4. Band diagram of structures (a) BN at the linear chain positions of one layer, (b) one layer is graphyne-like BN-sheet and other is pristine, (c) both layers is graphyne-like BN-sheet

Figure 5. DOS and PDOS diagram of (a) BN at the linear chain positions of one layer, (c) one

ACS Paragon Plus Environment

16

Page 17 of 40

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

layer is graphyne-like BN-sheet and the other is pristine, and (i) each layer is a graphyne-like BN-sheet Also in these systems, as in the pristine bilayer graphyne, the states at the top of the valence band and at the bottom of the conduction band are mainly due to the pz orbitals. Carbon contributes to both valence and conduction bands; nitrogen contributes more to the valence band and boron to the conduction band [see Figure 5 and Supporting Information Figure S6]. Actually, due to the large electronegativity difference between boron and nitrogen, charge is mainly transferred from boron to nitrogen, as confirmed by Hirshfield charge68 and Bader charge69-71 analysis [see the Supporting Information Table S4-S7]. Calculated charge analysis shows that a negative electron cloud resides around the N atom while electron deficiency is observed near the B atom. Thus all kind of B-N bonds acquires some ionic character while all C-C bonds are covalent in nature (Figure 6). The presence of B-N unit in graphyne not only transfer charge between B and N atom but also modify the charge cloud of neighbouring carbon atom. For example, Bader charge analysis shows that, in ‘graphyne with BN at chain’ monolayer of ‘pristine + graphyne with BN at chain’, the negative electron cloud around the nitrogen not only attributed by the transferred electrons from boron but some electrons are also transferred to nitrogen from the adjacent carbon atom. Similarly, a fractional amount of charge is also transferred from boron to adjacent carbon atom. Thus, in case of systems where the dopant BN sits at chain sites, the carbon atom adjacent to nitrogen also loses some charge and acts as a donor while the carbon atom adjacent to boron gains some charge and acts as an acceptor [see the Supporting Information Table S6]. The opposite effect is observed for systems with BN at ring sites [Supporting Information Table S7]. For example in ’graphyne with BN at ring’ of ‘pristine + graphyne with BN at ring’ bilayer, the

ACS Paragon Plus Environment

17

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 18 of 40

carbon atom adjacent to boron acts as a donor by transferring electron to boron and the carbon atom adjacent to nitrogen behaves as an acceptor by acquiring electron.

Figure 6. Nature of bond in all homo and hetero bilayer structures (a) covalent C-C bond (b) BN bond having ionic character Assessment of potential of graphyne bilayer as a nanocapacitor: Given the great interest in nanocapacitors, 10-13,50-51 in this section we are going to consider the possibility to employ these systems as nanocapacitors. Nanoscale capacitors store energy by accumulating opposite electric charges in two plates separated by a distance. The charge separation is achieved by applying an external electric field. It is therefore natural to consider our bilayer systems as nanocapacitors. To this aim, we have calculated dipole moment, capacitance,

ACS Paragon Plus Environment

18

Page 19 of 40

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

excess charge stored on each plate and energy stored by the capacitor. The total charge on each plate is obtained from DFT calculations by using the formula =P/d

(4)

where P is the dipole moment and d is the interlayer distance, which acts as separation thickness between two plates of the capacitor. The separation of the two charged plates is kept at the nano scale and the energy stored in the nanocapacitor can be obtained from the energy difference between the total energy under electric field and total energy without field: (5)

If m is the total mass of all atoms in unit cell, then the capacitance per unit mass is expressed by (6) In fact, in the hetero-bilayers considered here, a charge transfer takes place between the two layers, even in absence of an external electric field, resulting in a permanent dipole and in an electrostatic attraction between the two layers. The transferred charge is however quite small, upto 0.03 electronic charges per unit cell, and therefore the effect of this charge on the capacitive behaviour is minor. Still, this results in a difference of electrostatic potential of 0.2 to 0.5 V between the two layers, which might be useful if exploited in electronic devices. On the contrary, in homo-bilayer systems there is no charge transfer without an external electric field. With increasing electric field, the charge storage between two plates increases linearly [see Figure 7 (a) and 8 (a)]. Our calculation shows that for both homo bilayer and hetero bilayer, the charge stored on each plate (Qs) with same electric field is comparatively larger than graphene/hBN/graphene dielectric nanocapacitor with several h-BN layers sandwiched between two

ACS Paragon Plus Environment

19

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

graphene plates proposed by

Page 20 of 40

z elik et al.10 In addition, the calculated energy storing capacity

(Ec) [see Figure 7 (b) and 8 (b)] shows the superiority of our proposed homo bilayer and hetero bilayer over other proposed nanocapacitor.10 It indicates that more energy can be stored in our considered periodic bilayer structure in comparison with graphene/h-BN/graphene structure,10 which can be released in the external circuit. The capacitance of homo bilayer systems as a function of increasing electric field is depicted in Figure 7 (c). For electric fields up to ±0.2 V/ , dipole moment [see Supporting Information Figure S7 (a)], stored charge and energy are nearly the same for all homo-bilayer systems. At larger fields, differences start to be evident. We observe dependence of the capacitance from the electric field, with values remaining in the range between 25 and 41 F/g in the whole range of electric fields. At non-zero electric field, the capacitance of the BN doped homo-bilayers remains below that of the pristine graphyne bilayer.

Figure 7. Calculated (a) energy stored (b) charge stored (c) capacitance of homo-bilayer graphyne derivatives in dependence of an external electric field Also in the case of hetero-bilayers, the capacitance is smaller than that of the pristine bilayer graphyne and depends on the electric field. An interesting feature of these bilayers is that the capacitive behaviour is not symmetric with respect to a sign change of the electric field. This is probably due to the charge transfer between the two layers present in absence of an external electric field, which is responsible for a built-in electric field. We assume that this charge

ACS Paragon Plus Environment

20

Page 21 of 40

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

transfer is also responsible for the fluctuating behaviour of the capacitance in a region around zero electric field, giving rise to a Z shaped fluctuation region. A small electric field may invert the charge transfer, thereby inducing a large change of charge with a small change of energy. This results in a large capacitance at small electric fields, which then decreases again as the reversal of charge transfer is completed.

Figure 8. Calculated (a) energy stored (b) charge stored (c) capacitance of hetero-bilayer graphyne derivatives in dependence of an external electric field These bilayers show an interesting behaviour which can be exploited in nanocapacitors. The value of the capacitance of hetero bilayer systems is lower than the pristine bilayer graphyne [except ‘pristine + graphyne with BN at chain’ with positive electric field], however, our calculated capacitance value is higher than other systems reported in the literature, such as the graphene-h-BN-graphene systems considered in Ref. 10 except the case when only two h-BN layers are sandwiched between graphene layers. Indeed, the graphene-h-BN-graphene was shown to possess a capacitance of 5.6-24.5 F/g, while here we predict values in the range of 2040 F/g, depending on the exact structure and the magnitude of the external electric field. Moreover, the hetero-bilayers display a built-in electric field due to charge transfer between the

ACS Paragon Plus Environment

21

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 22 of 40

layers, which make the capacitive behaviour asymmetric with respect to the direction of the external electric field, and might be useful on its own for electronic applications.

CONCLUSIONS In summary, structure and electronic properties of bilayer graphyne have been calculated using density functional theory with van der Waals corrections. The highest stability for pristine bilayer systems is found for the A⍺Bβ configuration, while the AA configuration (with the two layers on top of each other) is the most unstable. The stacking arrangement plays a crucial role in determining the gap. The most stable structure, A⍺Bβ, is a semiconductor with a small band gap of 0.38 eV, which is much smaller than the band gap of monolayer graphyne. In contrast, the second most stable structure has a band gap of 0.50 eV, larger than monolayer graphyne. Interestingly, the AA stacked graphyne has a tiny band gap of 0.12 eV. Doping with boron nitride changes the electronic properties so that BN-doped bilayer graphyne becomes a wide band gap semiconductor. Boron and nitrogen doping increases the band gap of the most stable configuration and the increase depends on the BN substitution site. Such properties can make bilayer graphyne a useful and flexible material for applications in electronics, as for field effect transistors and for the creation of heterostructures, if spatially selective doping with a high degree of control can be achieved. Finally, we have investigated the capacitive behavior of the bilayers. They display a high capacitance value, and high energy storage capacity, in comparison to other 2D carbon-based materials. Moreover, the asymmetric behaviour of the capacitance with respect to the sign of the applied bias could also lead to intriguing applications in electronics.

ACS Paragon Plus Environment

22

Page 23 of 40

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Summarizing, our work predicts pristine and boron nitride doped bilayer graphyne systems to display a variety of properties that could in principle be controlled and tuned at the nanoscale, and should be useful for nanoelectronic applications. ASSOCIATED CONTENT Supporting Information The Supporting Information includes details information of structural parameters (lattice constant, interlayer distance, binding energy (BE), and energy band gap) of pristine graphyne bilayer with different stacking mode including with and without van der Waals correction. It also contains the band structure plot; DOS and PDOS plot and charge analysis of some BN doped graphyne bilayer system. This material is available free of charge via the Internet at http://pubs.acs.org. AUTHOR INFORMATION Corresponding Author Dr. Utpal Sarkar 1, E-mail address: [email protected] (Utpal Sarkar), Tel. +919401542687 Dr. Nicola Seriani 2, E-mail address: [email protected] (Nicola Seriani), Tel. +39 040 2240 279 Notes: The authors declare no competing financial interest. Author Contributions The manuscript was written through contributions of all authors. All authors have given approval to the final version of the manuscript.

ACS Paragon Plus Environment

23

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 24 of 40

ACKNOWLEDGMENT This research is supported by Assam University, Silchar, India and The Abdus Salam International Centre for Theoretical Physics, Trieste, Italy. REFERENCES 1) Novoselov, K. S.; Geim, A. K.; Morozov, S. V.; Jiang, D.; Zhang, Y.; Dubonos, S. V.; Grigorieva, I. V.; Firsov, A. A. Electric Field Effect in Atomically Thin Carbon Films. Science 2004, 306, 666-669. 2) Kroto, H. W.; Allaf, A. W.; Balm, S. P. C60: Buckminsterfullerene. Chem. Rev. 1991, 91, 1213-1235. 3) Iijima, S. Helical Microtubules of Graphitic Carbon. Nature 1991, 354, 56-58. 4) Geim, A. K.; Novoselov, K. S. The Rise of Graphene. Nat. Mater. 2007, 6, 183-191. 5) Moon, P.; Koshino, M. Electronic Properties of Graphene/Hexagonal-Boron-Nitride Moiré Superlattice. Phys. Rev. B 2014, 90, 155406 (1-12). 6) Abergel, D. S. L.; Mucha-Kruczyński, M. Infrared Absorption of Closely Aligned Heterostructures of Monolayer and Bilayer Graphene with Hexagonal Boron Nitride. Phys. Rev. B 2015, 92, 115430 (1-6). 7) Ansari, R.; Malakpour, S.; Ajori, S. Structural and Elastic Properties of Hybrid Bilayer Graphene/h-BN with Different Interlayer Distances Using DFT. Superlattices and Microstructures 2014, 72, 230–237. 8) Cavallucci, T.; Tozzini, V. Multistable Rippling of Graphene on SiC: A Density Functional Theory Study. J. Phys. Chem. C 2016, 120, 7670−7677.

ACS Paragon Plus Environment

24

Page 25 of 40

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

9) Ca, Y.; Zhang, G.; Zhang, Y. W. Electronic Properties of Phosphorene/Graphene and Phosphorene/Hexagonal Boron Nitride Heterostructures. J. Phys. Chem. C 2015, 119, 13929−13936. 10) Özçelik, V. O.; Ciraci, S. Nanoscale Dielectric Capacitors Composed of Graphene and Boron Nitride Layers: A First-Principles Study of High Capacitance at Nanoscale. J. Phys. Chem. C 2013, 117, 15327−1533. 11) Monajjemi, M. Metal-doped Graphene Layers Composed with Boron Nitride–Graphene as an Insulator: A Nano-capacitor. J Mol Model 2014, 20, 2507 (1-8). 12) Shi, G.; Hanlumyuang, Y.; Liu, Z.; Gong, Y.; Gao, W.; Li, B.; Kono, J.; Lou, J.; Vajtai, R.; Sharma, P.; Ajayan, P. M. Boron Nitride−Graphene Nanocapacitor and the Origins of Anomalous Size-Dependent Increase of Capacitance. NanoLett. 2014, 14, 1739−1744. 13) Peymanirad, F.; Neek-Amal, M.; Beheshtian, J.; Peeters, F. M. Graphene-Silicene Bilayer: A Nanocapacitor with Permanent Dipole and Piezoelectricity Effect. Phys. Rev. B 2015, 92,155113 (1-7). 14) Bhattacharya, B.; Sarkar, U. Graphyne–Graphene (Nitride) Heterostructure as Nanocapacitor. Chem. Phys. 2016, doi:10.1016/j.chemphys.2016.05.004 15) Baughman, R. H.; Eckhardt, H.; Kertesz, M. Structure‐Property Predictions for New Planar Forms of Carbon: Layered Phases Containing sp2 and sp Atoms. J Chem Phys 1987, 87, 6687-6699. 16) Kim, B. G.; Choi, H. J. Graphyne: Hexagonal Network of Carbon with Versatile Dirac Cones. Phys. Rev. B 2012, 86, 115435(1-5). 17) Singh, N. B.; Bhattacharya, B.; Sarkar, U. A First Principle Study of Pristine and BNDoped Graphyne Family. Struct. Chem. 2014, 25, 1695–1710.

ACS Paragon Plus Environment

25

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 26 of 40

18) Narita, N.; Nagai, S.; Suzuki, S.; and Nakao, K. Optimized Geometries and Electronic Structures of Graphyne and Its Family, Phys. Rev B, 1998, 58, 11009-11014. 19) Kang, J.; Li, J.; Wu, F.; Li, S. S.; Xia, J. B. Elastic, Electronic, and Optical Properties of Two-Dimensional Graphyne Sheet, J. Phys. Chem. C 2011, 115, 20466–20470. 20) Zhou, J.; Lv, K.; Wang, Q.; Chen, X. S.; Sun, Q.; Jena, P. Electronic Structures and Bonding of Graphyne Sheet and Its BN Analog. J. Chem. Phys. 2011, 134, 174701-5. 21) Sarkar, U.; Bhattacharya, B.; Seriani, N. First Principle Study of Sodium Decorated Graphyne. Chem. Phys. 2015, 461, 74–80. 22) Koo, J.; Huang, B.; Lee, H.; Kim, G.; Nam, J.; Kwon, Y.; Lee, H. Tailoring the Electronic Band Gap of Graphyne. J. Phys. Chem. C 2014, 118, 2463−2468. 23) Soni, H. R.; Jha, P. K. Vibrational and Elastic Properties of 2D Carbon Allotropes: A First Principles Study. Solid State Commun. 2014, 189, 58–62. 24) Bhattacharya, B.; Singh, N. B.; Sarkar, U. Pristine and BN Doped Graphyne Derivatives for UV Light Protection. Int J Quantum Chem 2015, 115, 820–829. 25) Kehoe, J. M.; Kiley, J. H.; English, J. J.; Johnson, C. A.; Petersen, R. C.; Haley, M. M.; Carbon Networks Based on Dehydrobenzoannulenes.3. Synthesis of Graphyne Substructures. Org. Lett. 2000, 2, 969–972. 26) Johnson, A. C.; Lu, Y.; Haley, M. M. Carbon Networks Based on Benzocyclynes. 6. Synthesis of Graphyne Substructures via Directed Alkyne Metathesis. Org. Lett. 2007, 9, 3725−3728. 27) Yoshimura, T.; Inaba, A.; Sonoda, M.; Tahara, K.; Tobe, Y.; Williams, R. V. Synthesis and Properties of Trefoil-Shaped Tris-(hexadehydrotribenzo[12]annulene) and Tris(tetradehydrotribenzo[12]annulene).Org. Lett. 2006, 8, 2933−2936.

ACS Paragon Plus Environment

26

Page 27 of 40

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

28) Haley, M. M. Synthesis and Properties of Annulenic Subunits of Graphyne and Graphdiyne Nanoarchitectures. Pure Appl. Chem. 2008, 80, 519–532. 29) Liu, H. B.; Xu, J. L.; Li, Y. J.; Li, Y. L. Aggregate Nanostructures of Organic Molecular Materials. Accounts Chem. Res., 2010, 43, 1496-1508. 30) Li, G. X.; Li, Y. L.; Liu, H. B.; Guo, Y. B.; Li, Y. J.; Zhu, D. B. Architecture of Graphdiyne Nanoscale Films. Chem. Commun. 2010, 46, 3256−3258. 31) Haley,

M.

M.;

Brand,

S.

C.;

Park,

J.

J.

Carbon

Networks

Based

on

Dehydrobenzoannulenes: Synthesis of Graphdiyne Substructures. Angew. Chem., Int. Ed. 1997, 36, 836−838. 32) Li, G.; Li, Y.; Qian, X.; Liu, H.; Lin, H.; Chen, N.; Li, Y. Construction of Tubular Molecule Aggregations of Graphdiyne for Highly Efficient Field Emission. J. Phys. Chem. C 2011, 115, 2611–2615. 33) Qian, X.; Ning, Z.; Li, Y.; Liu, H.; Ouyang, C.; Chen, Q.; Li, Y. Construction of Graphdiyne Nanowires with High-conductivity and Mobility. Dalton Transactions, 2012, 41, 730-733. 34) Zheng, Q.; Luo, G.; Liu, Q.; Quhe, R.; Zheng, J.; Tang, K.; Gao, Z.; Nagase, S.; Lu, J. Structural and Electronic Properties of Bilayer and Trilayer Graphdiyne. Nanoscale 2012, 4, 3990-6. 35) Majidi, R.; Karami, A. Electronic Properties of Bilayer and Trilayer Graphyne in the Presence of Electric Field. Struct. Chem. 2014, 25, 853-858. 36) Majidi, R.; Karami, A. Electronic Properties of BN-doped Bilayer Graphene and Graphyne in the Presence of Electric Field. Molecular Physics 2013, 111, 3194-3199.

ACS Paragon Plus Environment

27

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 28 of 40

37) Leenaerts, O.; Partoens, B.; Peeters, F. M. Tunable Double Dirac Cone Spectrum in Bilayer α-Graphyne. Appl. Phys. Lett. 2013, 103, 013105(1-4). 38) Hang, Y.; Wu, W.Z.; Yu, J.; Guo, W. L. Tuning the Energy Gap of Bilayer -Graphyne by Applying Strain and Electric Field. Chin. Phys. B. 2016, 25, 023102(1-4). 39) Narita, N.; Nagai, S.; Suzuki, S.; Nakao, K. Electronic Structure of Three-Dimensional Graphyne. Phys. Rev. B 2000, 62, 11146-11151. 40) Hwang, H. J.; Koo, J.; Park, M.; Park, N.; Kwon, Y.; Lee, H. Multilayer Graphynes for Lithium Ion Battery Anode. J. Phys. Chem. C 2013, 117, 6919−6923. 41) Stengel, M.; Spaldin, N. A. Origin of the Dielectric Dead Layer in Nanoscale Capacitors. Nature 2006, 443, 679−682 42) Stoller, M. D.; Park, S.; Zhu, Y.; An, J.; Ruoff, R. S. Graphene Based Ultracapacitors. Nano Lett. 2008, 8, 3498−3502. 43) Pandolfo, A. G.; Hollenkamp, A. F. Carbon Properties and Their Role in Supercapacitors. J. Power Sources 2006, 157, 11−27. 44) Liu, C. G.; Liu, M.; Li, F.; Cheng, H. M. Frequency Response Characteristics of SingleWalled Carbon Nanotubes as Supercapacitor Electrode Material. Appl. Phys. Lett. 2008, 92, 143108-3. 45) Niu, C. M.; Sichel, E.; Hoch, R.; Moy, D.; Tennent, H. High Power Electrochemical Capacitors Based on Carbon Nanotube Electrodes. Appl. Phys. Lett. 1997, 70, 1480-1482. 46) Tang, H.; Sodano, H. A. Ultra High Energy Density Nanocomposite Capacitors with Fast Discharge Using Ba0.2Sr0.8TiO3 Nanowires. Nano Lett. 2013, 13, 1373−1379. 47) Sorel, S.; Khan, U.; Coleman, J. N. Flexible, Transparent Dielectric Capacitors with Nanostructured Electrodes. Appl. Phys. Lett. 2012, 101, 101106.

ACS Paragon Plus Environment

28

Page 29 of 40

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

48) Vivekchand, S. R. C.; Rout, C. S.; Subrahmanyam, K. S.; Govindaraj, A.; Rao, C. N. Graphene-Based Electrochemical Supercapacitors. J. Chem. Sci. 2008, 120, 9−13. 49) Yoo, J. J.; Balakrishnan, K.; Huang, J.; Meunier, V.; Sumpter, B. G.; Srivastava, A.; Conway, M.; Reddy, A. L. M.; Yu, J.; Vajtai, R.; Ajayan, P. M. Ultrathin Planar Graphene Supercapacitors. Nano Lett. 2011, 11, 1423−1427. 50) Subramaniam, C.K.; Maiyalagan, T. Double Layer Energy Storage in Graphene - A Study. Micro and Nanosystems, 2012, 4, 180-185. 51) Ji, H.; Zhao, X.; Qiao, Z.; Jung, J.; Zhu, Y.; Lu, Y.; Zhang, L. L.; MacDonald, A. H.; Ruoff, R. S.; Capacitance of Carbon-Based Electrical Double-Layer Capacitors. Nat Commun 2014, 5, 3317(1-7). 52) Lee, K.; Murray, É. D.; Kong, L.; Lundqvist, B. I.; Langreth, D. C. Higher-Accuracy Van der Waals Density Functional. Phys. Rev. B 2010, 82, 081101(1-4). 53) Giannozzi, P. QUANTUM ESPRESSO: A Modular and Open-Source Software Project for Quantum Simulations of Materials. J. Phys. Condens. Matter 2009, 21, 395502(1-19). 54) Vanderbilt, D. Soft Self-consistent Pseudopotentials in a Generalized Eigenvalue Formalism. Phys. Rev. B 1990, 41, 7892-7895. 55) Marzari, N.; Vanderbilt, D.; Vita, A. D.; Payne, M. C. Thermal Contraction and Disordering of the Al(110) Surface. Phys. Rev. Lett. 1999, 82, 3296-3299. 56) Soler, J. M.; Artacho, E.; Gale, J. D.; García, A.; Junquera, J.; Ordejón, P.; Portal, D. S. The SIESTA Method for Ab Initio Order-N Materials Simulation. J. Phys. Condens. Matter 2002, 14, 2745–2779. 57) Agapito, L. A.; Fornari, M.; Ceresoli, D.; Ferretti, A.; Curtarolo, S.; Nardelli, M. B. Accurate Tight-Binding Hamiltonians for Two-Dimensional and Layered materials. Phys. Rev. B 2016, 93, 125137 (1-8).

ACS Paragon Plus Environment

29

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 30 of 40

58) Luo, G.; Zheng, Q.; Mei, W. N.; Lu, J.; Nagase, S. Structural, Electronic, and Optical Properties of Bulk Graphdiyne. J. Phys. Chem. C 2013, 117, 13072−13079. 59) Ukpong, A. M. First Principles Study of Van Der Waals Heterobilayers. Computational. Condens. Matter 2015, 2, 1-10. 60) Hasegawa, M.; Nishidate, K.; Iyetomi, H. Energetics of Interlayer Binding in Graphite: The Semiempirical Approach Revisited. Phys. Rev. B 2007, 76, 115424 (1-8). 61) Xiao, J.; Long, M.; Li, X.; Zhang, Q.; Xu, H.; Chan, K. S. Effects of Van Der Waals Interaction and Electric Field on the Electronic Structure of Bilayer MoS 2. J. Phys.: Condens. Matter 2014, 26, 405302 (1-8). 62) Birowska, M.; Milowska, K.; Majewski, J. A. Van Der Waals Density Functionals for Graphene Layers and Graphite. Acta Phys Pol A, 2011, 120, 845-848. 63) Tian, Y. H; Huang, J.; Sheng, X.; Sumpter, B. G.; Yoon, M.; Kertesz, M. NitrogenDoping Enables Covalent-Like п-п Bonding between Graphenes. Nano Lett. 2015, 15, 5482-5491. 64) Hoffmann, R. Solids and Surfaces: A Chemist’s View of Bonding in Extended Structures. Wiley-VCH, New York, 1988. 65) Hod, O. Graphite and Hexagonal Boron-Nitride Have the Same Interlayer Distance. Why? J. Chem. Theory Comput. 2012, 8, 1360−1369. 66) Marom, N.; Bernstein, J.; Garel, J.; Tkatchenko, A.; Joselevich, E.; Kronik, L.; Hod, O. Stacking and Registry Effects in Layered Materials: The Case of Hexagonal Boron Nitride. PRL 2010, 105, 046801 (1-4). 67) Leven, I.; Maaravi, T,; Azuri, I.; Kronik, L.; Hod, O. Interlayer Potential for Graphene/hBN Heterostructures. J. Chem. Theory Comput. 2016, 12, 2896-2905.

ACS Paragon Plus Environment

30

Page 31 of 40

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

68) Hirshfeld, F. L. Bonded-Atom Fragments for Describing Molecular Charge Densities. Theoret. Chim. Acta 1977, 44, 129-138. 69) Bader, R. F. W.; Henneker, W. H.; and Cade, P. E. Molecular Charge Distributions and Chemical Binding. J. Chem. Phys. 1967, 46, 3341-3363. 70) Bader, R. F. W. Atoms in Molecules - A Quantum Theory, Oxford University Press, Oxford, U.K., 1990. 71) Henkelman, G.; Arnaldsson, A.; Jónsson, H. A Fast and Robust Algorithm for Bader Decomposition of Charge Density. Comput. Mater. Sci. 2006, 36, 254-360.

Table of Contents (TOC)

ACS Paragon Plus Environment

31

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 1. Optimized geometries with unit cell of: (a) AA stacked graphyne; (b) AB stacked graphyne; (c) A⍺Aβ stacked graphyne 71x28mm (300 x 300 DPI)

ACS Paragon Plus Environment

Page 32 of 40

Page 33 of 40

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Figure 2. Band structure and projected density of states of the A⍺Aβ stacked graphyne 177x81mm (300 x 300 DPI)

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 3. Optimized geometry of: (a) most stable hetero layer “pristine + graphyne like BN sheet” 65x23mm (300 x 300 DPI)

ACS Paragon Plus Environment

Page 34 of 40

Page 35 of 40

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Figure 3. Optimized geometry of: (b) most stable homo layer “graphyne like BN sheet + graphyne like BN sheet” 66x24mm (300 x 300 DPI)

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 4. Band diagram of structures (a) BN at the linear chain positions of one layer, (b) one layer is graphyne-like BN-sheet and other is pristine, (c) both layers is graphyne-like BN-sheet 177x58mm (300 x 300 DPI)

ACS Paragon Plus Environment

Page 36 of 40

Page 37 of 40

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Figure 5. DOS and PDOS diagram of (a) BN at the linear chain positions of one layer, (c) one layer is graphyne-like BN-sheet and the other is pristine, and (i) each layer is a graphyne-like BN-sheet 177x61mm (300 x 300 DPI)

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 6. Nature of bond in all homo and hetero bilayer structures (a) covalent C-C bond (b) B-N bond having ionic character. 177x138mm (300 x 300 DPI)

ACS Paragon Plus Environment

Page 38 of 40

Page 39 of 40

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Figure 7. Calculated (a) energy stored (b) charge stored (c) capacitance of homo-bilayer graphyne derivatives in dependence of an external electric field 177x59mm (300 x 300 DPI)

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 8. Calculated (a) energy stored (b) charge stored (c) capacitance of hetero-bilayer graphyne derivatives in dependence of an external electric field 177x62mm (300 x 300 DPI)

ACS Paragon Plus Environment

Page 40 of 40