Next Article in Journal
Effect of Cu on Performance of Self-Dispersing Ni-Catalyst in Production of Carbon Nanofibers from Ethylene
Next Article in Special Issue
Large-Separation Behavior of the Casimir–Polder Force from Real Graphene Sheet Deposited on a Dielectric Substrate
Previous Article in Journal
Electrical Resistance Evolution of Graphite and Talc Geological Heterostructures under Progressive Metamorphism
Previous Article in Special Issue
Casimir–Polder Force on Atoms or Nanoparticles from Gapped and Doped Graphene: Asymptotic Behavior at Large Separations
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Effective Quantum Graph Models of Some Nonequilateral Graphyne Materials

by
César R. de Oliveira
*,† and
Vinícius L. Rocha
Departamento de Matemática, Universidade Federal de São Carlos, Cx. P. 676, São Carlos 13565-905, SP, Brazil
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Submission received: 25 May 2023 / Revised: 10 July 2023 / Accepted: 3 August 2023 / Published: 8 August 2023

Abstract

:
It is shown that it is possible to adapt the quantum graph model of graphene to some types of nonequilateral graphynes considered in the literature; we also discuss the corresponding nanotubes. The proposed models are, in fact, effective models and are obtained through selected boundary conditions and an ad hoc prescription. We analytically recover some results from the literature, in particular, the presence of Dirac cones for α -, β - and ( 6 , 6 , 12 ) -graphynes; for γ -graphyne, our model presents a band gap (according to the literature), but only for a range of parameters, with a transition at a certain point with quadratic touch and then the presence of Dirac cones.
PACS:
81.05.Uw; 73.63.Bd; 73.43.Cd

1. Introduction

There are many proposals of two-dimensional carbon allotropes in the literature, with graphene being the most prominent due to its peculiar electronic properties, and it has been synthesized since 2004 (see the review [1]). Theoretical techniques to study graphene, both single and multilayer cases, include ab initio calculations (comprising the density functional theory) and tight-binding models. On the rigorous mathematical side, there are three approaches for single layer graphene in the literature: one continuous model by Fefferman and Weinstein [2] proving the “generic” presence of Dirac cones in the dispersion relations; interesting symmetry arguments by Berkolaiko and Comech [3], and this method applies to different settings, including some multilayer models; and a quantum graph model (QGM), mainly due to Amovilli, Leys and March [4] and Kuchment and Post [5]. The latter approach has two advantages: it is simpler then the previous one and permits an accurate spectral analysis including details of Dirac cones; however, this technique requires that all edges in the graph have the same length, that is, an equilateral graph. Recently, the present authors have proposed adaptations in the boundary conditions so that the AA- and AB-stacked graphene multilayer could (at least qualitatively) be modeled by QGMs [6,7]; there is a work that also covers a few layers of graphene combined with hexagonal boron nitride [8] (again with a QGM).
In this work, we discuss how to (approximately) describe some distinct structures of graphynes via QGMs, including cases with edges of different lengths, again by playing with boundary conditions, and our main interest is in the possible presence of Dirac cones; ad hoc assumptions will be employed, and our results are compatible with some known results in the literature. We also discuss graphyne nanotubes in this context. Although our approach may seem rather exploratory, the results are mathematically correct, and it was not without some surprise that such proposal has worked!
The structure of graphynes are obtained from that of graphene. Recall that all carbon atoms in the graphene equilateral hexagonal lattice (we reserve the term honeycomb to such a regular hexagonal lattice) have sp 2 hybridization with single and double bonds. Carbon atoms in graphynes have both sp 2 and sp hybridizations, and they are usually obtained by replacing selected single C ( sp 2 ) C ( sp 2 ) bonds in graphene with
C ( sp 2 ) C ( sp ) C ( sp ) C ( sp 2 )
ones [9], which then also include triple bonds. Graphynes were predicted by Baughman, Eckhardt and Kertesz [10] in 1987 and only small samples have been synthesized; however, recently [11,12], a scalable synthesis was obtained, and some authors advocate that graphyne might replace graphene [13].
Some graphynes are expected to present Dirac cones, which were theoretically found in first-principle calculations [14], tight-binding models [9] and, in a particular case for which all edges have the same length [15], through a QGM (and Dirac cones were found). Perhaps this work could be seen as a step forward for some parts of [15].
Here, we extend the QGM to four types of graphynes considered in [9,14] (see also [16]); they are shown in Figure 1 and named α -, β -, γ - and (6,6,12)-graphynes. In short, we always think of the graphene honeycomb lattice, and for each edge, we introduce a positive parameter in the boundary conditions (see (4)) at vertices, in order to take into account the intensity of the chemical bond in each edge; the stronger the bond, the larger the parameter value. This will permit us to model some graphynes whose edges have different lengths, i.e., nonequilateral ones, in terms of the honeycomb structure; it should be clear that we have just effective models, which include some heuristic arguments, since the honeycomb structure of such models may differ from the graphyne structures (and this is a strength of our approach).
After modeling with QGMs, in many instances, the calculations are rigorous and analytical, so they are without approximations (sometimes we appeal to the plots of the graphs of functions to guide us about the possible existence of zeros). We summarize our findings as follows:
  • α -graphyne has the peculiarity of edges with three double bonds, and others with two single and a triple bonds. We argue (see ahead) that, from the point of view of QGMs, they have comparable bond strengths, so our model reproduces the one of graphene [5] in this case. Hence, we do not explore α -graphene in this work (we just include some remarks for completeness with respect to [9,14]).
  • In accordance with other works [9,14], we have obtained that Dirac points are present in β -graphyne.
  • In the literature [9], a band gap was found (with no Dirac cone, of course) for the γ -graphyne. Our QGM shows a richer structure in this case; there is a transition at a certain parameter value t = t c (see ahead for details), with a gap band between valence and conduction bands for 0 < t < t c , Dirac cones for t c < t < 1 and at the transition point t = t c , there is a parabolic touch. We have an explicit description of the gap size as a function of the parameters. Such properties should be of interest to be confirmed by experimentalists, in particular because it opens a potential for practical applications of γ -graphyne.
    It is worth mentioning that this is one of the graphynes in which a scalable synthesis has been recently obtained [11].
  • Our effective QGM of (6,6,12)-graphyne also presents Dirac cones for all values of the parameters, with two different Dirac points, in agreement with [9,14]. This case is technically more involved than the others, so that we have checked that there are no other touch points by looking at the graph of the dispersion relations (instead of just analytical expressions as in the other cases). We have not found that it is “self-doped” as reported in [14] (note that this property was not reported in a tight-binding calculation [9] either).
  • For the graphynes we consider, we discuss conditions so that the Dirac cones stay present after building graphyne tubes.
  • In all cases with Dirac cones, they are explicitly located and the cone slopes are given in terms of the parameters of the model.
In Section 2, the general idea of the QGM is presented in the specific case of the honeycomb lattice (since we propose to (effectively) think of all considered cases with the structure of this lattice), with some heuristics and the proposed choices of boundary conditions. In Section 3, a general discussion about dispersion relations of such QGMs is performed, which is then, in Section 4, specialized to Dirac cones in each considered graphyne type. Nanotubes are the subject of Section 5. Conclusions appear in Section 6.

2. Honeycomb Quantum Graph Model

2.1. Graphene Graph Model: A Short Account

Now, we briefly discuss the graphene quantum graph employed by Amoville, Leys and March [4] and mainly by Kuchment and Post [5] (see these works for precise details). Let G be the honeycomb lattice, the hexagonal 2D structure shown in Figure 2a (see Section 2.1.1), and let the Hilbert space L 2 ( G ) : = e L 2 ( e ) consist of all the square integrable functions on G. The free (i.e., no potential) graphene Hamiltonian H acting in functions u L 2 ( G ) along each edge e is given by
H u e ( x ) = d 2 u e d x 2 ( x ) ,
satisfying the so-called Neumann vertex conditions at the vertices. These conditions requires the continuity of the functions and vanishing of the total flux at each vertex v, that is, u e 1 ( v ) = u e 2 ( v ) = u e 3 ( v ) and u e 1 ( v ) + u e 2 ( v ) + u e 3 ( v ) = 0 , respectively.
From the Floquet–Bloch theory, the study of the Hamiltonian H is reduced to the study of the family of Bloch Hamiltonians H ( θ ) , where θ = ( θ 1 , θ 2 ) is the quasimomentum in the Brillouin zone B : = [ π , π ] 2 . These operators H ( θ ) act the same way H does, but in functions that satisfy, in addition to the Neumann vertex conditions, the Floquet condition (see (6)). Moreover, H ( θ ) has a purely discrete spectrum and the graph of the function θ σ ( H ( θ ) ) is the dispersion relation of H. The range of this function is the spectrum of H, that is, σ ( H ) is given by the reunion of σ ( H ( θ ) ) , θ B .
Thus, in order to determine the spectrum of H and its dispersion relation, we have to solve the eigenvalue problem
H ( θ ) u = λ u , λ R ,
which is done by employing two auxiliary operators, that is, one of them with the action (1) and Dirichlet boundary condition, and the other a periodic (one-dimensional) Hill operator.
After calculations, one obtains an explicit dispersion relation depending on the function (compare with (15))
g ( θ ) = ( 1 + e i θ 1 + e i θ 2 ) ( 1 + e i θ 1 + e i θ 2 ) = 1 + 8 cos θ 1 θ 2 2 cos θ 1 2 cos θ 2 2 .
In particular, it can be shown that the dispersion relation has Dirac cones at the points ± ( 2 π / 3 , 2 π / 3 ) of the Brillouin zone.
In the following, particularly in Section 3 and Section 4, we shall present details of this method adapted and applied to some graphynes.

2.1.1. Graphyne Sheet Model

The honeycomb lattice G is generated by the union of two triangular sublattices, g A (with vertices of type A) and g B (with vertices of type B), where
g A : = Z E 1 + Z E 2 and g B : = g A + ( 1 , 0 ) ,
with E 1 = ( 3 / 2 , 3 / 2 ) and E 2 = ( 0 , 3 ) being the lattice vectors (see Figure 2a). Consider the action of the group Z 2 on G using shift by the vectors p 1 E 1 + p 2 E 2 , p = ( p 1 , p 2 ) Z 2 , with a fundamental domain (or Wigner–Seitz cell) W = { v 1 , v 2 , a 1 , a 2 , a 3 } (see Figure 2b).
We identify each edge e of G with the segment [ 0 , 1 ] , denote by E ( G ) the set of edges of G and by E v ( G ) , the set of three edges that contains the vertex v.
The proposed Hamiltonian operator H acts along each edge e E ( G ) as the free operator (see Remark 4) in suitable units,
H u e ( x ) = d 2 u e d x 2 ( x ) , x e ,
satisfying the modified Neumann vertex conditions in each vertex v of G; let E v ( G ) = { e 1 , e 2 , e 3 } :
(i)
Modified continuity condition:
u e 1 ( v ) t 1 = u e 2 ( v ) t 2 = u e 3 ( v ) t 3 .
(ii)
Modified zero total flux condition:
t 1 u e 1 ( v ) + t 2 u e 2 ( v ) + t 3 u e 3 ( v ) = 0 .
The derivatives u e j ( v ) are directed from v to the other vertex of e j and t 1 , t 2 and t 3 are the positive interaction parameters between the edge bonds. Our idea is that the larger the bond strength in an edge (so the larger the parameter value), the larger its contribution to the flux at each vertex and the larger the value of the function at the vertex (due to the division by the parameter). The values of such parameters depend on each considered graphyne.
We call the general operator H a graphyne operator. Under these conditions, the graphyne operator is self-adjoint, as discussed in Appendix A.
Remark 1. 
From (4), one may suppose that the maximum value among the parameters t 1 , t 2 , t 3 takes the value 1. For instance, in case t 1 takes the maximum value, one may write (4) in the form
t 1 u e 1 ( v ) + t ˜ 2 u e 2 ( v ) + t ˜ 3 u e 3 ( v ) = 0 ,
with 0 < t ˜ 2 = t 2 / t 1 1 and 0 < t ˜ 3 = t 3 / t 1 1 .
Remark 2. 
Note that graphene, as studied in [5], is recovered by selecting t 1 = t 2 = t 3 . This simple remark has an important consequence for this work: it tells us that in the QGM of graphene, the single and double bonds have effectively the same intensity, and we will assume this while modeling graphynes. We have checked, in some cases, that our results are essentially the same whether we slightly distinguish the intensities of single and double bonds (not shown here), but we pay the price of much more complicated expressions to deal with.
By taking into account the above remarks, let us describe some heuristics and the proposed parameter ranges in each graphyne case shown in Figure 1.
(i)
The lattice is always the honeycomb one, i.e., the equilateral hexagonal lattice in Figure 2a; we technically require that all edges have the same length (e.g., due to the representation (18));
(ii)
For each graphyne, we look at its structure (Figure 1) and select edges that will be associated with a 1 , a 2 , a 3 in the fundamental domain W of the honeycomb lattice (Figure 2b), thus identifying the values of parameters t 1 , t 2 , t 3 (recall that the larger the bond strength, the larger the corresponding parameter; note that a 1 differs from the other edges since it is the unique that connects the two vertices in W).
(iii)
By looking at Figure 1, we identify four types of edge bonds b 1 , , b 4 , whose bond strengths will be probed (and compared) by their well-known enthalpy values (http://www.wiredchemist.com/chemistry/data/bond_energies_lengths.html, accessed on 23 May 2023):
  • b 1 representing C C and known enthalpy h 1 = 346 kJ/mol;
  • b 2 representing C = C and enthalpy h 2 = 602 kJ/mol;
  • b 3 representing C = C = C = C and enthalpy h 3 = 1806 kJ/mol (this bond only occurs in α -graphyne);
  • b 4 representing C C C C and enthalpy h 4 = 1527 kJ/mol (this bond is part of all graphyne compositions).
In graphene, there occurs only b 1 and b 2 bonds, and they are considered indistinguishable from the point of view of QGMs [5]; their enthalpy difference is h 2 h 1 = 256 kJ/mol. Hence, bonds with an enthalpy difference of this order will not be distinguished in our modeling, and this is the case of b 3 and b 4 , whose enthalpy difference is 279 kJ/mol (and ( h 3 h 4 ) / ( h 2 h 1 ) 1.09 ). The other enthalpy values differ from at least 925 kJ/mol (the difference between b 2 and b 4 with ( h 4 h 2 ) / ( h 2 h 1 ) 3.6 ), and so they will be considered distinct in the modeling.
(iv)
We need a criterion to associate two vertices and three graphyne edges in its fundamental domain (Figure 1) with those in W that could depict the main graphyne characteristics. It is natural to pick the most common configuration that appears in the fundamental cell of each graphyne; for β - and γ -graphynes, this procedure works, but not for the ( 6 , 6 , 12 ) -graphyne. So, for the latter case, we have worked with the options (that includes the bond b 4 ) and selected the one that has recovered results in [9,14].
Let’s apply the above procedure to the graphynes in Figure 1; see selections that are labeled a 1 , a 2 , a 3 in that figure:
  • α -graphyne: as anticipated, all edge bonds are supposed to be indistinguishable from the QGM viewpoint, so t 1 = t 2 = t 3 = 1 ; it then coincides with the graphene QGM (which is consistent to known results).
  • β -graphyne: we have t 2 = t 3 = 1 (since a 2 ,   a 3 have b 4 bonds) and 0 < t 1 < 1 (because a 1 has a b 2 bond).
  • γ -graphyne: we have t 1 = 1 and 0 < t 2 = t 3 < 1 .
  • ( 6 , 6 , 12 ) -graphyne: we take t 1 = t 2 = 1 and 0 < t 3 < 1 .
Note that the choice associated with the a 1 edge, in all cases, point to a corner of the fundamental domain (there are two possibilities for the ( 6 , 6 , 12 ) case), but we do not have any justification to take this as a criterion to select graphyne edges to be mapped to the honeycomb fundamental domain.
Remark 3. 
With the benzene ring being a resonant structure, its bonds, represented by successive alternating single and double bonds, are rigorously equivalent. Therefore, considering them equivalent has both experimental and theoretical foundations from a quantum mechanical point of view. This way, graphynes can be reduced to two different parameters: one associated with the benzene ring and the other with the acetylene group C C C C . However, the use of enthalpy values give quantitative relations with a potential to be applied to other situations.
Remark 4. 
As done in previous works [5,6,7,15], similar qualitative results are obtained by adding an even potential to the operator (2). We have opted to keep things simpler by selecting a free operator in each edge; furthermore, this choice gives explicit expressions for some spectral quantities.

2.1.2. Graphyne Nanotube Model

A carbon nanotube is a tube made of carbon with a diameter in nanoscale. There are many types of these materials (see [17] for a more detailed discussion of nanotubes), but here, we will keep our attention to the single-walled carbon nanotubes (SWCNTs) that can be obtained by cutouts of a two-dimensional hexagonal lattice. Depending on the cutout, we have different classifications of the SWCNTs. For instance, Figure 3a shows the so-called zig-zag nanotube, since it has a zig-zag closed path around the tube. Figure 3b shows the armchair nanotube, which is encircled by a closed path in the form of an armchair. Other types of nanotubes considered here are the chiral nanotubes.
Here, we discuss the propose of the corresponding single-wall graphyne nanotube models. A more detailed discussion and classification of graphene nanotubes can be find in [17] and the references therein.
Let p R 2 { 0 } be a vector of the lattice of translation symmetries of the quantum graph G, that is, we can write p as
p = p 1 E 1 + p 2 E 2 , p 1 , p 2 Z .
Let p be the equivalence relation given by z 1 p z 2 if, and only if, z 2 z 1 = q . p , with z 1 , z 2 G and q as an integer number. The graph G p obtained as the quotient of G with respect to the equivalence relation p is proposed to model a graphyne nanotube. This graph is naturally isometrically embedded into the cylinder R 2 / p . If p = ( p 1 , p 2 ) , we denote G p = G ( p 1 , p 2 ) . There are several types of nanotubes. For instance, G ( N , 0 ) and G ( 0 , N ) are the so-called zig-zag nanotubes, G ( N , N ) are the armchair nanotubes and the chiral are nanotubes with the form G ( p 1 , p 2 ) , with p 1 p 2 , p 1 , p 2 0 .
The corresponding Hamiltonian operator will be denoted by H p , called the graphyne nanotube operator, which is defined exactly as the graphyne operator H above, with modified Neumann vertex conditions (2), and taking into account the above additional symmetry (see Section 5).

3. Dispersion Relation

Now, we derive the dispersion relation of the graphyne operator H; it is based on the Floquet–Bloch theory [18,19,20,21] combined with the idea [5] of considering spectral points outside the spectrum of the edge Dirichlet operator. We begin with the general model; then, we specialize to each graphyne case.
For each quasimomentum  θ = ( θ 1 , θ 2 ) in the Brillouin zone B, let H ( θ ) be the Bloch Hamiltonian acting on functions satisfying the conditions (3) and (4) and also the Floquet condition
u ( x + p 1 E 1 + p 2 E 2 ) = e i ( p 1 θ 1 + p 2 θ 2 ) u ( x ) ,
for any p = ( p 1 , p 2 ) Z 2 and x W . The spectra of H ( θ ) is constituted only of eigenvalues and denoted by σ ( H ( θ ) ) = { λ n ( θ ) } n , and the spectrum of H is the union of these spectra (see [18]):
σ ( H ) = θ B σ ( H ( θ ) ) .
Recall that the function θ { λ n ( θ ) } n is the dispersion relation of H. Thus, to determine σ ( H ) and the dispersion relation of H, we have to solve the following eigenvalue problem for each quasimomentum θ B :
H ( θ ) u = λ u , λ R .
In order to solve (8), consider the auxiliary operator H D , with action
H D u ( x ) = d 2 u d x 2 ( x )
and the Dirichlet boundary condition, that is, u ( 0 ) = u ( 1 ) = 0 . It is well known that H D has the purely discrete spectrum σ ( H D ) = { k 2 π 2 } k 1 . If λ σ ( H D ) , then there exist two linearly independent solutions, φ λ , 0 and φ λ , 1 , of the eigenvalue problem
φ = λ φ , λ R ,
such that
φ λ , 0 ( 0 ) = φ λ , 1 ( 1 ) = 1 ,
φ λ , 0 ( 1 ) = φ λ , 1 ( 0 ) = 0 ,
φ λ , 1 ( x ) = φ λ , 0 ( 1 x ) , x [ 0 , 1 ] .
Explicitly,
φ λ , 0 ( x ) = sin ( λ ( 1 x ) ) sin λ , φ λ , 1 ( x ) = sin ( λ x ) sin λ .
The quotient
η ( λ ) : = φ λ , 1 ( 1 ) φ λ , 1 ( 0 ) = cos λ
is well defined for λ σ ( H D ) .
Let λ σ ( H D ) and 0 < t j 1 , with j = 1 , 2 , 3 . Then, we claim that the real number λ belongs to the spectrum of H if, and only if, there exists a quasimomentum θ B such that
η ( λ ) = ± f ( θ ) T ,
with T = t 1 2 + t 2 2 + t 3 2 and
f ( θ ) = t 1 4 + t 2 4 + t 3 4 + 2 t 1 2 t 2 2 cos ( θ 1 ) + 2 t 1 2 t 3 2 cos ( θ 2 ) + 2 t 2 2 t 3 2 cos ( θ 1 θ 2 ) .
In order to conclude this, first, note that by the Floquet condition (6), we have (see Figure 2b)
u a 1 ( v 2 ) = e i θ 2 u a 2 ( v 3 ) = e i θ 2 u a 3 ( v 4 ) .
Thus, combining (16) with the modified Neumann vertex conditions, we obtain
u a 1 ( 0 ) / t 1 = u a 2 ( 0 ) / t 2 = u a 3 ( 0 ) / t 3 : = A u a 1 ( 1 ) / t 1 = e i θ 1 u a 2 ( 1 ) / t 2 = e i θ 2 u a 3 ( 1 ) / t 3 : = B t 1 u a 1 ( 0 ) + t 2 u a 2 ( 0 ) + t 3 u a 3 ( 0 ) = 0 t 1 u a 1 ( 1 ) + t 2 e i θ 1 u a 2 ( 1 ) + t 3 e i θ 2 u a 3 ( 1 ) = 0 .
Thus, for u to be an eigenfunction of H ( θ ) , it must satisfy (17).
Let λ σ ( H D ) and consider φ λ , 0 and φ λ , 1 as the linearly independent solutions of (10) that satisfy (11)–(13). We can represent u a j as
u a 1 = t 1 ( A φ λ , 0 + B φ λ , 1 ) u a 2 = t 2 ( A φ λ , 0 + e i θ 1 B φ λ , 1 ) u a 3 = t 3 ( A φ λ , 0 + e i θ 2 B φ λ , 1 ) .
It is easy to see that the function defined by (18) satisfies the modified continuity condition (3) of the modified Neumann vertex condition and solves the eigenvalue problem (10). It remains to be verified that (18) satisfies the zero total flux condition (4). Substituting (18) into the last two equations of (17), we obtain
( t 1 2 + t 2 2 + t 3 2 ) φ λ , 0 ( 0 ) A + ( t 1 2 + t 2 2 e i θ 1 + t 3 2 e i θ 2 ) φ λ , 1 ( 0 ) B = 0 ( t 1 2 + t 2 2 e i θ 1 + t 3 2 e i θ 2 ) φ λ , 0 ( 1 ) A + ( t 1 2 + t 2 2 + t 3 2 ) φ λ , 1 ( 1 ) B = 0 .
With (13), φ λ , 0 ( 0 ) = φ λ , 1 ( 1 ) and φ λ , 0 ( 1 ) = φ λ , 1 ( 0 ) . Substituting this into (19), then dividing by φ 1 ( 0 ) 0 and multiplying by 1 its second equation, we obtain
T η ( λ ) A + F ¯ ( θ ) B = 0 F ( θ ) A T η ( λ ) B = 0 ,
with F ( θ ) = t 1 2 + t 2 2 e i θ 1 + t 3 2 e i θ 2 , F ¯ the complex conjugate of F and T = t 1 2 + t 2 2 + t 3 2 . The determinant δ of this system equals
δ = T 2 η ( λ ) 2 f ( θ ) ,
with f ( θ ) = F ( θ ) F ¯ ( θ ) , which is exactly (15). Hence, if there exists a quasimomentum θ B such that η ( λ ) is one solution of δ = 0 , that is,
η ( λ ) = ± f ( θ ) T ,
it follows that the representation (18) solves the eigenvalue problem (8) and so λ σ ( H ) by (7). Therefore, the above claim is justified.
By (14), except possibly for λ σ ( H D ) , a discrete sequence of numbers, the dispersion relation of H is given by
cos ( λ ) = ± f ( θ ) T , θ B ,
with f ( θ ) given by (15). This description of the dispersion relation (23) of the graphyne operator H will allow us to study the possible presence of Dirac cones. In the next section, we present such an analysis.

4. Dirac Cones

We make use of (23) to study the possible presence of Dirac cones. We specialize in each case: β -graphyne, γ -graphyne and ( 6 , 6 , 12 ) -graphyne. Recall that a Dirac cone is a point where, in the lowest order approximation, the valence and conduction bands linearly touch each other, and the quasimomentum θ D B for which a Dirac cone occurs is called a D-point. In symbols, if θ D B is a D-point, then there is a constant γ 0 so that
λ ( θ ) λ ( θ D ) = ± γ | θ θ D | + O ( | θ θ D | 2 ) + O ( λ ( θ ) λ ( θ D ) 2 ) ,
with the “−” and “+” signs for the valence and conduction bands, respectively.
Let λ σ ( H D ) . By (23), if η ± ( λ , θ ) are the two roots of (21), given by
η ± ( λ , θ ) = ± f ( θ ) t 1 2 + t 2 2 + t 3 2 ,
in order to obtain Dirac cones, one must find D-point candidates θ D B and expand D ( λ ) = cos ( λ ) and η ± ( λ , θ ) around λ ( θ D ) and θ D , respectively. Then, if θ D is a D-point, expanding D ( λ ) around λ ( θ D ) ,
D ( λ ( θ ) ) ) = D ( λ ( θ D ) ) + D ( λ ( θ D ) ) ( λ ( θ ) λ ( θ D ) ) + O ( ( λ ( θ ) λ ( θ D ) ) 2 ) ,
which implies
D ( λ ( θ ) ) D ( λ ( θ D ) ) = c ( θ D ) ( λ ( θ ) λ ( θ D ) ) + O ( ( λ ( θ ) λ ( θ D ) ) 2 ) ,
with c ( θ D ) = λ ( θ D ) sin λ ( θ D ) / λ ( θ D ) ; finally,
λ ( θ ) = arccos η ± ( λ , θ ) + k π 2 , k Z .
Hence, if θ D is a D-point candidate, we have the explicit parameter c ( θ D ) . Therefore, it remains to be analyzed, in each type of graphyne, the possible presence of D-points.
In what follows, we will use the notations H β , H γ and H ( 6 , 6 , 12 ) to represent the β -, γ - and ( 6 , 6 , 12 ) -graphyne operators, respectively. The same indication will be employed to the roots η ± and the linear coefficients γ of (24), with   = β , γ , ( 6 , 6 , 12 ) . Recall that the α -graphyne QGM coincides with the graphene one.
  • β -graphyne
Recall that by inspecting the β -graphyne structure in Figure 1b, along with the associated bonds in the honeycomb fundamental domain in Figure 2b, we have proposed the following relations between parameters t 1 < t 2 = t 3 = 1 . Hence, the roots of (21) are given by
η ± β ( λ , θ ) = ± f ( θ ) 2 + t 1 2 ,
with
f ( θ ) = f β ( θ ) = 2 + t 1 4 + 2 t 1 2 ( cos ( θ 1 ) + cos ( θ 2 ) ) + 2 cos ( θ 1 θ 2 ) .
We have found that the roots η + β ( λ , θ ) and η β ( λ , θ ) touch each other in two points of the Brillouin zone B. Indeed, let the diagonal
B d : = { θ B : θ 1 = θ 2 }
in the Brillouin zone. Then, for θ B d , we have that η + β ( λ , θ ) = η β ( λ , θ ) = 0 if, and only if,
  f ( θ )   =   2 + t 1 4 + 2 t 1 2 cos ( θ 1 ) + cos ( θ 2 ) + 2 cos ( θ 1 θ 2 )     =   2 + t 1 4 + 2 t 1 2 cos ( θ 1 ) + cos ( θ 1 ) + 2 cos ( θ 1 + θ 1 ) (32)   =   4 cos 2 ( θ 1 ) + 4 t 1 2 cos ( θ 1 ) + t 1 4 (33)   =   ( t 1 2 + 2 cos ( θ 1 ) ) 2 = 0
which occurs in ± θ 1 D , with θ 1 D = arccos ( t 1 2 / 2 ) π . We will confirm now that these points θ 1 D are, in fact, D-points of the dispersion relation of the β -graphyne operator H β by expanding η ± β ( λ , θ ) around θ 1 D (around θ 1 D is similar). Let θ B d .
By expanding t 1 2 + 2 cos ( θ 1 ) around θ 1 D , we obtain
t 1 2 + 2 cos ( θ 1 ) = b ( θ D ) θ 1 θ 1 D + O ( θ 1 θ 1 D 2 ) ,
with b ( θ D ) = 2 sin θ 1 D . It then follows that
η ± β ( λ , θ 1 ) η ± β ( λ , θ 1 D ) = ± | b ( θ D ) | 2 + t 1 2 | θ 1 θ 1 D | + O | θ 1 θ 1 D | 2 .
Thus, combining (27) and (35),
λ ( θ ) λ ( θ D ) = ± γ β | θ θ D | + O | θ θ D | 2 + O ( λ ( θ ) λ ( θ D ) ) 2 ,
that is, (24) holds with
γ β = | b ( θ D ) | ( 2 + t 1 2 ) c ( θ D ) , θ D = ( θ 1 D , θ 1 D )
For θ D , the process is similar.
Therefore, the dispersion relation of the β -graphyne operator H β has two Dirac cones in the Brillouin zone (see Figure 4) for all 0 < t 1 < 1 , which occur at the D-points ± θ D , where θ D = ( θ 1 D , θ 1 D ) and θ 1 D = arccos ( t 1 2 / 2 ) π . Also, by analyzing the computational 3D plot of the dispersion relation, it was observed that the Dirac cones of the β -graphyne operator H β are situated at the vertices of a hexagon.
  • γ -graphyne
As already mentioned, the structure of the γ -graphyne in Figure 1c indicates the relations 1 = t 1 > t 2 = t 3 > 0 . In this case, the roots of (21) are
η ± γ ( λ , θ ) = ± f ( θ ) 1 + 2 t 2 2 ,
with
f ( θ ) = f γ ( θ ) = 1 + 2 t 2 4 + 2 t 2 2 ( cos ( θ 1 ) + cos ( θ 2 ) ) + 2 t 2 4 cos ( θ 1 θ 2 ) .
Here, we will also consider θ B d , that is, θ 2 = θ 1 . Thus, the roots (37) are equivalent to
η ± γ ( λ , θ ) = ± | 1 + 2 t 2 2 cos θ 1 | 1 + 2 t 2 2 .
Now, we analyze the possible presence of Dirac cones in the dispersion relation of the γ -graphyne operator H γ . We study (39) in three situations:
(i)
0 < t 2 < 2 / 2 . In this case, the function 1 + 2 t 2 2 cos θ 1 > 0 , for all values of θ 1 . Thus,
η ± γ ( λ , θ ) = ± 1 + 2 t 2 2 cos θ 1 1 + 2 t 2 2 .
The minimum and maximum of η + γ ( λ , θ ) and η γ ( λ , θ ) , respectively, occur at θ 1 = ± π , with values η ± γ ( λ , ± π ) = ± 1 1 + 2 t 2 2 . Note that the behavior in these points is parabolic. Then, the roots η ± γ ( λ , θ ) do not touch each other and, therefore, the dispersion relation of H γ does not have any Dirac cones in this parameter range (see Figure 5a).
(ii)
t 2 = 2 / 2 . In this case,
η ± γ ( λ , θ ) = ± 1 2 ± cos θ 1 2 ,
with parabolic touches occurring at θ 1 ± π , with value η ± γ ( λ , θ 1 ) = 0 , which proves that the dispersion relation of H γ does not have any Dirac cones (see Figure 5b).
(iii)
2 / 2 < t 2 < 1 . Differently from the cases (i) and (ii), the Dirac cones are present in this situation. In fact, we have that η + γ ( λ , θ ) = η γ ( λ , θ ) = 0 if, and only if, θ 1 = ± θ 1 D , with θ 1 D = arccos ( 1 2 t 2 2 ) . Expanding η ± γ ( λ , θ ) around ± θ 1 D , in the analogous way we have done in the β -graphyne case, we obtain
η ± γ ( λ , θ ) η ± γ ( λ , θ D ) = ± γ ˜ γ | θ θ D | O | θ θ D | 2 ,
with θ D = ( θ 1 D , θ 1 D ) and
γ ˜ γ = 4 t 2 4 1 1 + 2 t 2 2 > 0 .
Analogously, one deals with θ D . Combining (27) and (42), we obtain
λ ( θ ) λ ( θ D ) = ± γ ¯ γ | θ θ D | + O | θ θ D | 2 + O ( λ ( θ ) λ ( θ D ) ) 2 ,
with γ ¯ γ = γ ˜ γ / c ( θ D ) . Note that the γ ¯ is the linear coefficient, while the upper index “ γ ” indicates that we are considering the γ -graphyne. Therefore, for this parameter range, the dispersion relation of the γ -graphyne operator H γ has Dirac cones in its dispersion relation (see Figure 5c).
By analyzing the computational 3D plot of the dispersion relation, it was observed that the Dirac cones of the H γ is similar to the β -graphyne case, i.e., the Dirac cones are also situated at the vertices of a hexagon.
  • ( 6 , 6 , 12 ) -graphyne
In the case of ( 6 , 6 , 12 ) -graphyne (see Figure 1d), the proposed parameter relations are t 1 = t 2 = 1 and 0 < t 3 < 1 . Hence,
η ± ( 6 , 6 , 12 ) ( λ , θ ) = ± f ( θ ) 2 + t 3 2 ,
with
f ( θ ) = f ( 6 , 6 , 12 ) ( θ ) = 2 + t 3 4 + 2 cos θ 1 + 2 t 3 2 cos θ 2 + cos ( θ 1 θ 2 ) .
This case is tricker than the previous ones; in β - and γ -graphynes, we have taken the restriction to the segment B d in the Brillouin zone, since it was observed that the (possible) Dirac cones were present in the diagonal θ 2 = θ 1 . However, this does not happen for ( 6 , 6 , 12 ) -graphyne; we have found that the two Dirac points occur at the line
θ 2 = r ( θ 1 ) : = ( h ( ) / ) θ 1 ,
with = arccos t 3 4 2 2 and
h ( z ) : = z arcsin 1 t 3 2 sin z π .
Thus, let B r be the following restriction of the Brillouin zone B:
B r : = { θ B : θ 2 = r ( θ 1 ) } .
We shall check now that the points
θ D + : = ( , h ( ) ) and θ D = ( , h ( ) ) ,
in B r , are D-points. Let θ B r . Expanding f ( θ ) around , we obtain
(48) f ( θ ) =   f ( ) + f ( ) ( θ 1 ) + f ( ) 2 ( θ 1 ) 2 + O ( θ 1 ) 3 (49)     =   f ( ) 2 ( θ 1 ) 2 + O ( θ 1 ) 3 ,
since f ( ) = f ( ) = 0 , with f meaning the derivative of f with respect to θ 1 . It implies that
η ± ( 6 , 6 , 12 ) ( λ , θ ) η ± ( 6 , 6 , 12 ) ( λ , ) = ± γ ˜ ( 6 , 6 , 12 ) | θ θ D + | + O ( θ 1 ) 2 ,
with γ ˜ ( 6 , 6 , 12 ) = f ( ) / ( 2 ( 2 + t 3 2 ) ) . Combining (50) and (27), we obtain
λ ( θ ) λ ( θ D + ) = γ ( 6 , 6 , 12 ) | θ θ D + | + O | θ θ D + | 2 + O ( λ ( θ ) λ ( θ D + ) ) 2 ,
with
γ ( 6 , 6 , 12 ) = γ ˜ ( 6 , 6 , 12 ) c ( θ D + ) , θ D + = ( , h ( ) ) .
Analogously, one deals with θ D .
Therefore, the dispersion relation of the ( 6 , 6 , 12 ) -graphyne operator H ( 6 , 6 , 12 ) has Dirac cones on the Brillouin zone (Figure 6), as found in [9] via a tight-binding model and in [14] via first principles calculations, and here for all allowed parameter values.
Differently from the β - and γ -graphyne, which the Dirac cones are situated at the vertices of a hexagon, the Dirac cones of the dispersion relation of ( 6 , 6 , 12 ) -graphyne occurs at the vertices of a rhombus. In a similar way, one checks that the other two Dirac cones occuring at the D-points are ± θ ¯ D , with θ ¯ D = ( ¯ , ¯ / 2 ) and ¯ = arccos ( 2 t 3 4 ) / 2 .

5. Graphyne Nanotubes

5.1. Spectra of Graphyne Nanotubes

Let p = ( p 1 , p 2 ) Z 2 { 0 } be a vector of the lattice of translation symmetries of the graph G, denoting the corresponding nanotube by G p = G ( p 1 , p 2 ) and the graphyne nanotube operator H p = H ( p 1 , p 2 ) , as defined at the end of Section 2.
As in Section 3, we apply the Floquet–Bloch theory that provides the Bloch Hamiltonian operators H p ( θ ) for each quasimomentum θ in the Brillouin zone B and the decomposition
σ ( H p ) = θ B H p ( θ ) ,
with σ ( H p ( θ ) ) = { λ p , n ( θ ) } n , purely discrete. Since a function u on G p lifts to a p-periodic function on G, that is,
u ( x + p 1 E 1 + p 2 E 2 ) = u ( x ) ,
then, by the Floquet condition (6), it follows that
p · θ = p 1 θ 1 + p 2 θ 2 2 π Z .
Hence, we consider the restriction B p B of the Brillouin zone given by
B p = { θ B : p · θ 2 π Z } .
Therefore, (51) turns into
σ ( H p ) = θ B p H p ( θ ) ,
and the dispersion relation for H p is just the dispersion relation of H (see Section 3) restricted to B p , that is, it is given by
D ( λ ) = ± 2 f ( θ ) T , λ σ ( H D ) , θ B p ,
where D ( λ ) = cos ( λ ) (see (23)), f ( θ ) = F ( θ ) F ¯ ( θ ) , F ( θ ) = t 1 2 + t 2 2 e i θ 1 + t 3 2 e i θ 2 and T = t 1 2 + t 2 2 + t 3 2 .

5.2. Dirac Cones

Now, we analyze the possible presence of Dirac cones in the dispersion relation (55) of the graphyne nanotube operator H p . As before, we analyze separately the β , γ and ( 6 , 6 , 12 ) -graphyne.
Let θ D B be a D-point of the graphyne operator H (see (24)). Since the dispersion relation of H p is the dispersion relation of H restricted to B p , then θ D is a D-point of the nanotube graphyne operator H p if, and only if, θ D belongs to the restriction B p . Thus, given the D-point θ D in each case of graphyne, we determine p in order to have θ D B p .
  • β -graphyne nanotube
As discussed in Section 4, the dispersion relation of the β -graphyne operator H β has D-points at ± θ D , with
θ D = ( θ 1 D , θ 1 D ) , θ 1 D = arccos ( t 1 2 / 2 ) π .
Hence, the condition for θ D B p reduces to (similarly for θ D )
p · θ D = p 1 · θ 1 D p 2 · θ 1 D = arccos ( t 1 2 / 2 ) ( p 1 p 2 ) + π ( p 1 p 2 ) 2 π Z .
The above condition implies that the difference p 1 p 2 must be an even integer, that is, p 1 p 2 = 2 q , with q Z . Then,
arccos ( t 1 2 / 2 ) ( p 1 p 2 ) = 2 arccos ( t 1 2 / 2 ) q = 2 π r , r Z ,
that is,
t 1 2 = 2 cos r q π , r , q Z .
Because of the parameter t 1 ( 0 , 1 ) , we have 0 < 2 cos ( r π / q ) < 1 . Let the function
g β ( x ) : = 2 cos ( π x ) .
It is easy to see that 0 < g β ( x ) < 1 if, and only if, x C β , with
C β : = n N { 0 } ( J n + J n ) ,
J n + = 1 3 + n , 1 2 + n , n even 1 2 + n , 2 3 + n , n odd J n = 1 2 n , 1 3 n , n even 2 3 n , 1 3 n , n odd .
Hence, t 1 ( 0 , 1 ) if, and only if, r π / q C β . For instance, if p = ( p 1 , p 2 ) is such that p 1 p 2 = 14 , then if r = 1 and q = 7 , it follows that r π / q = π / 7 J 0 + , which implies that 0 < t 1 2 = 2 cos ( π / 7 ) < 1 (and so t 1 ( 0 , 1 ) ) and, thus, p · θ D 2 π Z .
Therefore, ± θ D are D-points of the dispersion relation of the β -graphyne nanotube operator H p β , for p = ( p 1 , p 2 ) , if, and only if, there exists r , q Z , such that p 1 p 2 = 2 q and r π / q C β .
  • γ -graphyne nanotube
From Section 4, we know that ± θ D , with θ D = ( θ 1 D , θ 1 D ) , θ 1 D = arccos ( 1 / 2 t 2 2 ) , are D-points of γ -graphyne operator H γ if 2 / 2 < t 2 < 1 . Hence, given p = ( p 1 , p 2 ) Z 2 , then ± θ D are D-points of H p γ if, and only if, p · ( ± θ D ) B p . By following the steps in the discussion of the β -graphyne nanotube case, we obtain that
t 2 2 = 1 2 cos r q π ,
with r , q Z and p 1 p 2 = 2 q . Let
g γ ( x ) : = 1 / 2 cos ( π x ) .
Then, it follows that 2 / 2 < g γ ( x ) < 1 if, and only if, x C γ , with C γ = n Z J n , n odd and J n = ( n 1 / 3 , n + 1 / 3 ) . For instance, if p is such that p 1 p 2 = 6 , then for r = 1 and q = 3 , it follows that r π / q = π / 3 J 1 and, then t 2 ( 2 / 2 , 1 ) .
Therefore, for t 2 ( 2 / 2 , 1 ) , ± θ D are D-points of the dispersion relation of the γ -graphyne nanotube operator H p γ , for p = ( p 1 , p 2 ) , if, and only if, there exists r , q Z , such that p 1 p 2 = 2 q and r π / q C γ .
If t 2 = 2 / 2 , we have found that the dispersion relation of H γ has parabolic touches at ± ( π , π ) . By imposing the restriction B p to these points, we obtain that p 1 p 2 must be even. Thus, if p = ( p 1 , p 2 ) is such that p 1 p 2 = 2 q , q Z , then it follows that the dispersion relation of H p γ has parabolic touches at ± ( π , π ) .
  • ( 6 , 6 , 12 ) -graphyne nanotube
We have found in Section 4 that, for t 3 ( 0 , 1 ) ,
θ D ± = ( ± , h ( ± ) ) ,
with = arccos t 3 4 2 2 and h ( z ) = z arcsin 1 t 3 2 sin z π , which are D-points for H p ( 6 , 6 , 12 ) . Due to the complexity of the function h ( x ) and the value of , we were just able to (analytically) analyze the possible presence of Dirac cones in the dispersion relation of the zig-zag nanotubes G p ( 6 , 6 , 12 ) , with p = ( N , 0 ) and N Z . By imposing that θ D ± B p and following the steps of β -graphyne case, we obtain
t 3 2 = 2 cos q N 2 π + 2 , q Z .
Hence, t 3 ( 0 , 1 ) if, and only if, 2 q π / N C ( 6 , 6 , 12 ) , with
C ( 6 , 6 , 12 ) = n Z J n , J n = n 2 1 6 , n 2 + 1 6 , n odd .
Therefore, θ D ± are D-points for the zig-zag nanotube C ( N , 0 ) ( 6 , 6 , 12 ) if, and only if, there exists q Z such that 2 q π / N C ( 6 , 6 , 12 ) .

6. Summary and Conclusions

In some situations, like graphene, QGMs are useful as models that usually permit an analytical approach with explicit calculations. For a single graphene sheet, it was mainly developed in [4,5] and has become standard, but requires that all graph edges have the same length; it was also applied to a particular equilateral graphyne in [15]. Motivated by previous works on multilayer graphene [6,7], the present authors have suggested some adaptations (1) in the boundary conditions, (2) selection of graphyne edges and (3) heuristic choices of parameters that take into account the intensity of chemical bonds, which have permitted an effective graph modeling of some graphynes.
We have proposed QGMs of the graphynes discussed in [9,14], in order to investigate Dirac cones; the main differences to graphene are that usually their hexagonal structure have edges of different lengths which carry different chemical bonds (in each graphyne, there are two or three types of bonds that repeat periodically). In principle, it is not immediate how to model such situations via quantum graphs; hence, we had to make some hypotheses which were based on heuristic observations and ad hoc procedures. We have also discussed the presence of Dirac cones for the corresponding graphyne nanotubes.
Our general approach was to model through the honeycomb structure of graphene in all cases, but introducing a positive parameter for each edge: there are three parameters, corresponding to the three edge intensities in the associated fundamental domain W of the honeycomb lattice (see Figure 2b). The idea is that the stronger the bond in an edge, the more influence it should have on the boundary condition balance (see Equations (3) and (4)). Next, we have tried to propose a way to associate edges of the graphyne to the ones ( a 1 , a 2 , a 3 ) of W; of course, the choice must be a triple that appears in the graphyne that is being considered, and the idea was to select the one that is more abundant in each graphyne fundamental domain (which are indicated in Figure 1). However, this is not the case of ( 6 , 6 , 12 ) -graphyne, as explained in previous sections. It was then recovered, via QGMs, results from tight-binding and first-principle calculations in the literature, although our results for the γ -graphyne present a transition from bands with a gap to bands with Dirac cones, as a parameter changes.
Our effective models pick the most influential three-edge structure of each graphyne and, also to simplify technical matters, some carbon bonds are considered to have similar strengths (and thus imply the same value of the corresponding parameters in (3) and (4)). For example, the QGM of graphene in the literature does not distinguish single and double carbon bonds, whose difference of enthalpies is 256 kJ/mol; thus, for each graphyne, we have associated the same parameter value to edges whose difference of enthalpies is “close” to such value. This is particularly important (and it was a guide) for the α -graphyne, whose proposed QGM then coincides with the one of graphene.
It would be interesting whether the found transition in γ -graphyne could be replicated either experimentally or by varying parameters in other theoretical calculations (e.g., tight-binding ones). We finish by mentioning that our proposal, to adapt QGMs to some graphynes, is a combination of heuristics, effective models and, after modeling, our results are mathematically rigorous.

Author Contributions

Conceptualization, C.R.d.O. and V.L.R.; methodology, C.R.d.O. and V.L.R.; validation V.L.R.; formal analysis, V.L.R.; investigation, C.R.d.O. and V.L.R.; data curation, V.L.R.; writing—original draft preparation, V.L.R.; writing—review and editing, C.R.d.O.; visualization, V.L.R. All authors have read and agreed to the published version of the manuscript.

Funding

CRdO thanks the partial support by CNPq (a Brazilian government agency, under contract 303689/2021-8). VLR thanks the financial support by CAPES (a Brazilian government agency).

Data Availability Statement

No new data were created.

Conflicts of Interest

Authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
MDPIMultidisciplinary Digital Publishing Institute
DOAJDirectory of open access journals
TLAThree letter acronym
LDLinear dichroism

Appendix A. Self-Adjointness

In this appendix, we show how to conclude that the boundary conditions (3) and (4) imply that the graphyne operator H is self-adjoint. By Theorems 1.4.4 and 1.4.11 in [22], in order to check the self-adjointness of H, it is necessary and sufficient that, for each vertex v with E v ( G ) = { e 1 , e 2 , e 3 } , there exist two 3 × 3 matrix A v and B v such that:
(i)
The 3 × 6 matrix A v B v has maximal rank;
(ii)
The matrix A v B v * is self-adjoint, where B v * is the adjoint of B v ;
(iii)
A v F ( v ) = B v F ( v ) , where the vector F ( v ) and F ( v ) are given by
F ( v ) : = u e 1 ( v ) u e 2 ( v ) u e 3 ( v )
and
F ( v ) : = u e 1 ( v ) u e 2 ( v ) u e 3 ( v ) .
Consider a vertex v. We choose the following constant 3 × 3 -matrices
A = t 2 t 1 0 0 t 3 t 2 0 0 0 and B = 0 0 0 0 0 0 t 1 t 2 t 3 .
These matrices satisfy the three conditions (i), (ii) and (iii) above for nonzero parameters t 1 , t 2 , t 3 . E.g., A B * = 0 , which is self-adjoint. Therefore, the operator H is self-adjoint. Also note that the (iii) above implies the proposed boundary conditions (3) and (4).

References

  1. Castro Neto, A.H.; Guinea, F.; Peres, N.M.R.; Novoselov, K.S.; Geim, A.K. The electronic properties of graphene. Rev. Mod. Phys. 2014, 81, 109–162. [Google Scholar] [CrossRef] [Green Version]
  2. Fefferman, C.L.; Weinstein, M.I. Honeycomb lattice potentials and Dirac cones. J. Am. Math. Soc. 2012, 25, 1169–1220. [Google Scholar] [CrossRef] [Green Version]
  3. Berkolaiko, G.; Comech, A. Symmetry and Dirac points in graphene spectrum. J. Spectr. Theory 2018, 8, 1099–1147. [Google Scholar] [CrossRef] [Green Version]
  4. Amovilli, C.; Leys, F.; March, N. Electronic energy spectrum of two-dimensional solids and a chain of C atoms from a quantum network model. J. Math. Chem. 2004, 36, 93–112. [Google Scholar] [CrossRef]
  5. Kuchment, P.; Post, O. On the spectra of carbon nano-structures. Commun. Math. Phys. 2007, 275, 805–826. [Google Scholar] [CrossRef] [Green Version]
  6. de Oliveira, C.R.; Rocha, V.L. Dirac cones for bi- and trilayer Bernal-stacked graphene in a quantum graph model. J. Phys. A Math. Theor. 2020, 53, 505201. [Google Scholar] [CrossRef]
  7. de Oliveira, C.R.; Rocha, V.L. From multilayer AA-stacked graphene sheets to graphite: Graph models and Dirac cone. Z. Naturforsch. A 2021, 76, 371–384. [Google Scholar] [CrossRef]
  8. de Oliveira, C.R.; Rocha, V.L.; Souza, O.N. Boron nitride and graphene heterostructures modeled by quantum graphs. (preprint, submitted for publication).
  9. Liu, Z.; Yu, G.; Yao, H.; Liu, L.; Jiang, L.; Zheng, Y. A simple tight-binding model for typical graphyne structures. New J. Phys. 2012, 14, 113007. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  10. Baughman, R.H.; Eckhardt, H. Structure–property predictions for new planar forms of carbon: Layered phases containing sp2 and sp atoms. J. Chem. Phys. 1987, 87, 6687–6699. [Google Scholar] [CrossRef]
  11. Desyatkin, V.G.; Martin, W.B.; Aliev, A.E.; Chapman, N.E.; Fonseca, A.F.; Galvão, D.S.; Miller, E.R.; Stone, K.H.; Wang, Z.; Zakhidov, D.; et al. Scalable synthesis and characterization of multilayer γ-graphyne, new carbon crystals with a small direct band gap. J. Am. Chem. Soc. 2022, 144, 17999–18008. [Google Scholar] [CrossRef] [PubMed]
  12. Kang, J.; Wei, Z.; Li, J. Graphyne and its family: Recent theoretical advances. Acs Appl. Mater. Interfaces 2019, 11, 2692–2706. [Google Scholar] [CrossRef] [PubMed]
  13. Rawat, S. Graphene is a Nobel Prize-Winning “Wonder Material” Graphyne Might Replace it. Big Think. The Future, 5 August 2022. Available online: https://bigthink.com/the-future/graphyne/ (accessed on 12 May 2023).
  14. Malko, D.; Neiss, C.; Viñes, F.; Görling, A. Competition for graphene: Graphynes with direction-dependent Dirac cones. Phys. Rev. Lett. 2012, 108, 086804. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Do, N.T.; Kuchment, P. Quantum graph spectra of a graphyne structure. Nanoscale Syst. Math. Model. Theory Appl. 2013, 2, 107–123. [Google Scholar] [CrossRef]
  16. Enyanshin, A.; Ivanovskii, A. Graphene alloptropes: Stability, structural and electronic properties from DF-TB calculations. Phys. Status Solidi (b) 2011, 248, 1879–1883. [Google Scholar]
  17. Harris, P. Carbon Nano-Tubes and Related Structures; Cambridge University Press: Cambridge, UK, 2002. [Google Scholar]
  18. Brown, M.B.; Eastham, M.S.P.; Schmidt, K.M. Periodic Differential Operators; Birkhäuser: Basel, Switzerland, 2013. [Google Scholar]
  19. Eastham, M.S.P. The Spectral Theory of Periodic Differential Equations; Scottish Acad. Press: Edinburgh/London, UK, 1973. [Google Scholar]
  20. Kuchment, P. Floquet Theory for Partial Differential Equations; Birkhäuser: New York, NY, USA, 1993. [Google Scholar]
  21. Reed, M.; Simon, B. Methods of Modern Mathematical Physics IV: Analysis of Operators; Academic Press: New York, NY, USA, 1978. [Google Scholar]
  22. Berkolaiko, G.; Kuchment, P. Introduction to Quantum Graphs; American Mathematical Society: Providence, RI, USA, 2012. [Google Scholar]
Figure 1. Structures of a layer of different types of graphynes (as shown). Empty and filled balls represent sp and sp 2 carbon atoms, respectively. The dashed lines delimit the fundamental domain in each case. A possible choice of edges a 1 , a 2 , a 3 (see the main text) associated with parameters t 1 , t 2 , t 3 , respectively, is indicated in each case. (a) α -graphyne; (b) β -graphyne. (c) γ -graphyne. (d) (6,6,12)-graphyne.
Figure 1. Structures of a layer of different types of graphynes (as shown). Empty and filled balls represent sp and sp 2 carbon atoms, respectively. The dashed lines delimit the fundamental domain in each case. A possible choice of edges a 1 , a 2 , a 3 (see the main text) associated with parameters t 1 , t 2 , t 3 , respectively, is indicated in each case. (a) α -graphyne; (b) β -graphyne. (c) γ -graphyne. (d) (6,6,12)-graphyne.
Carbon 09 00076 g001
Figure 2. The honeycomb lattice G. (a) The lattice with lattice vectors E 1 and E 2 ; some type-A and type-B vertices are labeled. (b) Its fundamental domain W; it contains three edges, a 1 , a 2 and a 3 , and two vertices, v 1 and v 2 .
Figure 2. The honeycomb lattice G. (a) The lattice with lattice vectors E 1 and E 2 ; some type-A and type-B vertices are labeled. (b) Its fundamental domain W; it contains three edges, a 1 , a 2 and a 3 , and two vertices, v 1 and v 2 .
Carbon 09 00076 g002
Figure 3. Zig-zag and Armchair nanotubes. The red line indicates the path that defines the type of the nanotube. (a) (10,0) Zig-zag nanotube. (b) (5,5) Armchair nanotube.
Figure 3. Zig-zag and Armchair nanotubes. The red line indicates the path that defines the type of the nanotube. (a) (10,0) Zig-zag nanotube. (b) (5,5) Armchair nanotube.
Carbon 09 00076 g003
Figure 4. Dispersion relations of β -graphyne, with t 1 = 0.55 , restricted to the diagonal B d ; valence (solid line) and conduction (dashed line) bands, as well as two Dirac cones, are shown.
Figure 4. Dispersion relations of β -graphyne, with t 1 = 0.55 , restricted to the diagonal B d ; valence (solid line) and conduction (dashed line) bands, as well as two Dirac cones, are shown.
Carbon 09 00076 g004
Figure 5. Dispersion relations of γ -graphyne restricted to the diagonal B d ; solid and dashed lines illustrate valence and conduction bands, respectively. (a) Parameter t 2 = 0.55 and θ 1 [ 0 , 2 π ] ; there is no touch. (b) t 2 = 2 / 2 and θ 1 [ 0 , 2 π ] ; there is a parabolic touch. (c) t 2 = 0.84 and θ 1 [ 0 , π ] ; a Dirac cone is shown.
Figure 5. Dispersion relations of γ -graphyne restricted to the diagonal B d ; solid and dashed lines illustrate valence and conduction bands, respectively. (a) Parameter t 2 = 0.55 and θ 1 [ 0 , 2 π ] ; there is no touch. (b) t 2 = 2 / 2 and θ 1 [ 0 , 2 π ] ; there is a parabolic touch. (c) t 2 = 0.84 and θ 1 [ 0 , π ] ; a Dirac cone is shown.
Carbon 09 00076 g005
Figure 6. Dispersion relations of ( 6 , 6 , 12 ) -graphyne, with t 3 = 0.55 , restricted to the line B r ; valence (solid line) and conduction (dashed line) bands, as well as two Dirac cones, are shown.
Figure 6. Dispersion relations of ( 6 , 6 , 12 ) -graphyne, with t 3 = 0.55 , restricted to the line B r ; valence (solid line) and conduction (dashed line) bands, as well as two Dirac cones, are shown.
Carbon 09 00076 g006
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

de Oliveira, C.R.; Rocha, V.L. Effective Quantum Graph Models of Some Nonequilateral Graphyne Materials. C 2023, 9, 76. https://doi.org/10.3390/c9030076

AMA Style

de Oliveira CR, Rocha VL. Effective Quantum Graph Models of Some Nonequilateral Graphyne Materials. C. 2023; 9(3):76. https://doi.org/10.3390/c9030076

Chicago/Turabian Style

de Oliveira, César R., and Vinícius L. Rocha. 2023. "Effective Quantum Graph Models of Some Nonequilateral Graphyne Materials" C 9, no. 3: 76. https://doi.org/10.3390/c9030076

APA Style

de Oliveira, C. R., & Rocha, V. L. (2023). Effective Quantum Graph Models of Some Nonequilateral Graphyne Materials. C, 9(3), 76. https://doi.org/10.3390/c9030076

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