Next Article in Journal
Exploration of the Interaction Strength at the Interface of Anionic Chalcogen Anchors and Gold (111)-Based Nanomaterials
Next Article in Special Issue
The External Electric Field-Induced Tunability of the Schottky Barrier Height in Graphene/AlN Interface: A Study by First-Principles
Previous Article in Journal
Comparative Study of Different Acidic Surface Structures in Solid Catalysts Applied for the Isobutene Dimerization Reaction
Previous Article in Special Issue
Insights on the Excitation Spectrum of Graphene Contacted with a Pt Skin
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

TAO-DFT Study on the Electronic Properties of Diamond-Shaped Graphene Nanoflakes

by
Hong-Jui Huang
1,†,
Sonai Seenithurai
1,† and
Jeng-Da Chai
1,2,*
1
Department of Physics, National Taiwan University, Taipei 10617, Taiwan
2
Center for Theoretical Physics and Center for Quantum Science and Engineering, National Taiwan University, Taipei 10617, Taiwan
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Nanomaterials 2020, 10(6), 1236; https://doi.org/10.3390/nano10061236
Submission received: 18 May 2020 / Revised: 12 June 2020 / Accepted: 23 June 2020 / Published: 25 June 2020
(This article belongs to the Special Issue Science and Technology of Graphene)

Abstract

:
At the nanoscale, it has been rather troublesome to properly explore the properties associated with electronic systems exhibiting a radical nature using traditional electronic structure methods. Graphene nanoflakes, which are graphene nanostructures of different shapes and sizes, are typical examples. Recently, TAO-DFT (i.e., thermally-assisted-occupation density functional theory) has been formulated to tackle such challenging problems. As a result, we adopt TAO-DFT to explore the electronic properties associated with diamond-shaped graphene nanoflakes with n = 2–15 benzenoid rings fused together at each side, designated as n-pyrenes (as they could be expanded from pyrene). For all the n values considered, n-pyrenes are ground-state singlets. With increasing the size of n-pyrene, the singlet-triplet energy gap, vertical ionization potential, and fundamental gap monotonically decrease, while the vertical electron affinity and symmetrized von Neumann entropy (which is a quantitative measure of radical nature) monotonically increase. When n increases, there is a smooth transition from the nonradical character of the smaller n-pyrenes to the increasing polyradical nature of the larger n-pyrenes. Furthermore, the latter is shown to be related to the increasing concentration of active orbitals on the zigzag edges of the larger n-pyrenes.

1. Introduction

Graphene, a wonder material consisting of a hexagonal arrangement of carbon atoms in a plane, has revolutionized science and technology over the past few decades [1,2]. The covalent bonds formed from the s p 2 -hybridized carbon atoms are responsible for the structure of graphene. Besides, graphene exhibits the Dirac cone that is suitable for electronics applications. Graphene also has excellent electronic, mechanical, thermal, and optical properties. However, graphene has been found to be a zero-gap semiconductor or semimetal. Therefore, in order to operate graphene as a semiconductor material in electronic devices, it is essential to adopt some scenarios to introduce a band gap in graphene, such as defect formation, doping, functionalization, and many others.
Among them, one scenario to tune the properties of graphene is to cut an infinite graphene sheet into nano-sized graphene fragments, such as graphene nanoflakes (GNFs). GNFs can be classified based on their shapes (e.g., rectangular [3,4,5,6,7], disk-like [8], hexagon-shaped [9,10,11,12,13], triangle-shaped [11,14,15], bowtie-shaped [16], and many other shaped [17,18,19,20,21] GNFs) and their edges (e.g., zigzag and/or armchair [11] GNFs). In particular, zigzag GNFs could possess room-temperature magnetism as well as adjustable energy gaps [13]. Also, the optical and electronic properties associated with graphene nanoribbons (GNRs) (i.e., rectangular GNFs) could be tuned by using uniaxial strain [7]. GNRs (i.e., rectangular GNFs) have large hydrogen storage capacity [4] and interesting magnetic properties [3]. Besides, triangle-shaped GNFs linked by carbon chains can have intrinsic ferromagnetism [19]. In addition, topological frustration in GNFs can yield large net spin [17]. Moreover, the edge states and termination of edges using different atoms lead to interesting properties.
The structural versatility of GNFs gives rise to an interesting avenue for the exploration of the properties of GNFs with different shapes, edges, sizes, lengths/breadths, edge-termination, pores, and so forth [11,14,21,22]. The interesting electronic and magnetic properties and tunability of properties make these GNFs potential candidates for electronics, spintronics [16], nanoscale molecular logic gates [23], and so on. Besides, GNFs can be used as the building blocks of novel nanomaterials [19]. Large-spin GNFs can be used in nanoscale spintronics devices [17]. In addition, graphene quantum dots are potential candidates for spin memory, transistors, and solid-state qubits in quantum computer applications [24,25].
Despite flourishing activities in studying the properties and applications of different GNFs, there are only scarce studies on diamond-shaped GNFs (i.e., rhombic GNFs). A diamond-shaped GNF can be viewed as two interconnected triangle-shaped GNFs [15]. It is also worth mentioning that triangulene, which can be viewed as a small triangle-shaped GNF, has been recently synthesized by Pavliček et al. [26]. Since the larger triangle-shaped GNFs has been found to possess pronounced polyradical character [15], it may be of interest to investigate whether the polyradical character can persist for the closely related diamond-shaped GNFs. Therefore, in this work, we explore the electronic properties associated with diamond-shaped GNFs with n benzenoid rings fused together at each side (designated as n-pyrenes), which could be expanded from pyrene (see Figure 1). Note that n-pyrene (with the chemical formula C 2 n 2 + 4 n H 4 n + 2 ), which contains N = 12 n 2 + 28 n + 2 electrons, belongs to the class of PAHs (i.e., polycyclic aromatic hydrocarbons). Particularly, 2-pyrene (i.e., pyrene), a naturally occurring molecule in interstellar space [27], has been synthesized by Förster and Kasper [28]. However, owing to the difficulties in synthesizing the larger n-pyrenes, there are few experimental studies of the larger n-pyrenes.
Consequently, theoretical studies can be complementary, and could offer useful information on the properties of n-pyrenes (i.e., diamond-shaped GNFs of different sizes). Because of the pronounced radical character found in the closely related triangle-shaped GNFs [15] and because of their low dimensionality [29], n-pyrenes are expected to possess radical character, and hence they can be troublesome electronic systems for common computational approaches. For systems possessing radical character, KS-DFT (i.e., Kohn-Sham density functional theory) [30] with the commonly employed XC (i.e., exchange-correlation) energy functionals is prone to predict unreliable results [31,32]. In order to explore the properties associated with such electronic systems, one generally resorts to multi-reference (MR) electronic structure methods, including the complete-active-space self-consistent-field (CASSCF) and related methods [33,34,35,36,37,38,39]. While these MR electronic structure methods are reliably accurate, they are computationally expensive, and hence intractable for large electronic systems. Therefore, efficient methods are necessary for performing the electronic structure calculations on the larger n-pyrenes.
Recently, TAO-DFT (i.e., thermally-assisted-occupation density functional theory) [40] has been developed for studying the electronic properties associated with nanosystems exhibiting radical character. In strong contrast to KS-DFT, we emphasize that TAO-DFT is a DFT (i.e., density functional theory) employing fractional orbital occupation numbers (given by the Fermi-Dirac distribution function with the fictitious temperature θ ). In TAO-DFT, an entropy contribution component, which depends on the θ value as well as the orbital occupation numbers, is capable of providing an adequate description for strong static correlation even when the very simple local density approximation (LDA) XC energy functional is employed. Relatively complicated semilocal [41], common hybrid [42] or long-range corrected hybrid [42,43] XC energy functionals could be used for performing the TAO-DFT calculations as well. Note that TAO-DFT, which is similar to KS-DFT in computational complexity, reduces to KS-DFT for electronic systems exhibiting nonradical character. Besides, to improve the overall accuracy of TAO-DFT for diverse applications, an approach that determines the self-consistent value of θ in TAO-DFT has been recently formulated [44]. Also, to comment on its applicability, TAO-DFT has been widely adopted for electronic structure calculations [5,12,15,45,46,47,48,49,50,51,52,53], hydrogen storage applications [46,48,49], and vibrational analysis [54], especially for nanosystems possessing radical character. Furthermore, in several recent investigations [5,37,39,40,42,45], the occupation numbers of orbitals from TAO-DFT have been shown to be close to the occupation numbers of natural orbitals obtained with the variational two-electron reduced-density-matrix-driven CASSCF (v2RDM-CASSCF) method, which is an accurate MR electronic structure method, leading to a qualitatively similar tendency in describing the radical nature associated with various PAHs.
In particular, TAO-DFT has been employed to explore the electronic properties associated with zigzag GNFs of different shapes and sizes, including zigzag GNRs (i.e., rectangular GNFs) [5], hexagon-shaped GNFs [12], and triangle-shaped GNFs [15] in recent years. According to the findings, the electronic properties associated with these GNFs are very different. Therefore, here we continue employing TAO-DFT to explore the electronic properties associated with n-pyrenes with n = 2–15 (i.e., diamond-shaped GNFs of different sizes), which could pave the way for their potential applications in electronics and optoelectronics.

2. Computational Details

Geometry optimizations and single-point energy calculations are all performed using TAO-LDA [40], which is TAO-DFT employing the LDA θ -dependent energy functional and XC energy functional, where the recommended value of θ = 7 × 10 3 (hartree) is employed [40]. All computational results are obtained with Q-Chem 4.4 [55], adopting the 6-31G(d) basis set (i.e., a valence double-zeta polarized basis set) as well as a numerical quadrature which contains 75 grid points in the Euler-Maclaurin radial grid and 302 grid points in the Lebedev angular grid.
All calculations are performed under isolated boundary conditions (i.e., well suited for studying molecules). A convergence tolerance of 10 8 (hartree) is set for single-point energy calculations. For the convergence criteria of geometry optimizations, we adopt the default tolerances of Q-Chem (i.e., geometry optimizations are considered converged, when a convergence tolerance of 3 × 10 4 (hartree/bohr) on the maximum gradient component is reached, and either a convergence tolerance of 1.2 × 10 3 (bohr) on the maximum atomic displacement or a convergence tolerance of 10 6 (hartree) on the energy change of successive optimization cycles is reached).

3. Results and Discussion

3.1. Singlet-Triplet Energy Gap

To determine the ground state, we first optimize the structures of n-pyrene associated with the lowest singlet state and lowest triplet state, by spin-unrestricted TAO-LDA, and thereafter calculate the difference between the lowest singlet energy ( E singlet ) and lowest triplet energy ( E triplet ) associated with n-pyrene using [5,12,15,40,41,42,45,46,47,48,49,50,51]
E ST = E triplet E singlet ,
where E ST is the singlet-triplet energy gap of n-pyrene.
As can be seen in Figure 2, the E ST value of n-pyrene remains positive for each n value (also see Table S1 in Supplementary Information (SI)). Therefore, all the n-pyrenes studied are ground-state singlets. As the size of n-pyrene increases, the E ST value monotonically decreases at a faster rate for smaller n = 2–4, followed by a slower rate of monotonic decrease for larger values of n. Understanding the E ST values is crucially important for certain applications, such as the thermal energy conversion where the singlet-fission phenomenon is utilized [56]. Besides, exploring systems with extremely small (or even vanishing) values of E ST is also demanding for the thermally activated delayed fluorescence (TADF) applications [57,58,59]. On the basis of the E ST values reported, the larger n-pyrenes can be useful for the singlet-fission phenomenon and TADF applications.
KS-DFT employing the widely adopted XC energy functionals can suffer from the unphysical symmetry-breaking problems in the respective spin-unrestricted calculations [31,32], especially for systems with radical character. For an exact theory (where there is no unphysical symmetry-breaking problem) [36,40,60], the spin-unrestricted and spin-restricted calculations must yield identical energy. To investigate whether spin-unrestricted TAO-LDA resolves unphysical symmetry-breaking problems, we additionally perform geometry optimizations for the lowest singlet states of n-pyrenes by spin-restricted TAO-LDA. The requirement of symmetry can indeed be satisfied for n = 2–15 (i.e., for all the n values considered), as the TAO-LDA spin-unrestricted and spin-restricted singlet energies of n-pyrenes are essentially the same.

3.2. Vertical Ionization Potential/Electron Affinity as Well as Fundamental Gap

Now, we explore if n-pyrenes (containing N electrons) could be useful for applications in optoelectronics. At the ground-state geometry, we perform spin-unrestricted TAO-LDA calculations to compute the vertical ionization potential of ground-state n-pyrene using [5,12,15,41,42,46,47,48,49,50,51]
IP v = E N 1 E N ,
the vertical electron affinity of ground-state n-pyrene using
EA v = E N E N + 1 ,
and the fundamental gap (i.e., IP v EA v ) of ground-state n-pyrene using
E g = E N + 1 + E N 1 2 E N ,
where E N , E N 1 , and E N + 1 are the energies associated with neutral n-pyrene, cationic n-pyrene, and anionic n-pyrene, respectively.
The IP v (Figure 3), EA v (Figure 4), and E g (Figure 5) values of ground-state n-pyrene are presented as functions of n (also see Table S2 in SI). With the increase of the size of n-pyrene, the IP v value decreases monotonically, and the EA v value increases monotonically, which results in a monotonic decrease in the E g value (i.e., the difference between the IP v value and the EA v value). Interestingly, n-pyrenes (with n = 5–15) possess the E g values lying between 1 eV and 3 eV, showing promise for their applications in nanophotonics.

3.3. Symmetrized von Neumann Entropy

For a quantitative measure of the radical character of ground-state n-pyrene, we compute the symmetrized von Neumann entropy [5,12,15,41,42,46,47,48,49,50,51,60]
S vN = 1 2 σ = , i = 1 f i , σ ln ( f i , σ ) + ( 1 f i , σ ) ln ( 1 f i , σ ) ,
by spin-unrestricted TAO-LDA. Note that f i , σ , which ranges from zero to one, is the occupation number of the i th σ -spin (i.e., up-spin or down-spin) orbital from spin-unrestricted TAO-LDA, and is closely related to the occupation number of the i th natural orbital of σ -spin [40,41,42]. For an electronic system with nonradical character, the occupation numbers of all spin-orbitals are in the vicinity of either 0 or 1, making only a vanishingly small contribution to the respective value of S vN . Nevertheless, for an electronic system possessing pronounced radical character, the occupation numbers of active spin-orbitals (which are the spin-orbitals with significant fractional occupations) can be rather different from 0 and 1 (e.g., in the range of 0.1 to 0.9); hence, the respective value of S vN can be greatly increased when the number of active spin-orbitals increases and/or the occupation numbers associated with active spin-orbitals become closer to 0.5.
As presented in Figure 6, the S vN value of ground-state n-pyrene increases monotonically with increasing system size, which implies that the larger ground-state n-pyrenes could exhibit increasing polyradical nature (also see Table S2 in SI).

3.4. Active Orbital Occupation Numbers

In order to unlock the mystery of increasing S vN value with the size of n-pyrene, we present the active orbital occupation numbers of ground-state n-pyrene (with N electrons), computed using spin-restricted TAO-LDA (see Figure 7). Here, we define the HOMO (i.e., highest occupied molecular orbital) as the ( N / 2 ) th orbital, the LUMO (i.e., lowest unoccupied molecular orbital) as the ( N / 2 + 1 ) th orbital, and so on [5,12,15,40,42,45,47,50,51]. Besides, the active orbitals are regarded as the orbitals possessing an occupation number in the range of 0.2 to 1.8.
For smaller values of n (up to n = 3), the occupation numbers of all orbitals are in the vicinity of either 0 or 2. Accordingly, 2-pyrene and 3-pyrene are expected to possess nonradical character. However, when n increases, the number of active orbitals increases and/or the occupation numbers associated with active orbitals become closer to 1, indicating that the larger ground-state n-pyrenes could exhibit increasing polyradical nature (also see Table S3 in SI).
From the occupation number analysis, the smaller ground-state n-pyrenes (e.g., n 3 ) exhibit nonradical nature, which is consistent with the analyses of the other electronic properties of these stable n-pyrenes (e.g., the larger values of E ST , the larger values of E g , and the smaller values of S vN ). On the other hand, the larger ground-state n-pyrenes (e.g., n > 3 ) exhibit increasing polyradical nature, which is also consistent with the analyses of the other electronic properties of these relatively unstable n-pyrenes (e.g., the smaller values of E ST , the smaller values of E g , and the larger values of S vN ).

3.5. Real-Space Representation of Active Orbitals

We explore the real-space representation of active orbitals for the ground states of some representative n-pyrenes, such as 2-pyrene (see Figure 8), 4-pyrene (see Figure 9), 6-pyrene (see Figure 10), 8-pyrene (see Figure 11), and 10-pyrene (see Figure 12), computed using spin-restricted TAO-LDA. As shown, the active orbitals of the smaller n-pyrenes (e.g., n 3 ) appear delocalized over the whole n-pyrenes. Nonetheless, when n increases, more and more active orbitals become concentrated on the zigzag edges of the larger n-pyrenes, which is expected to be intimately related to the increasing polyradical nature of the larger n-pyrenes. This implies that when available, these edge states should be energetically more preferred than other states. Since more edge states become available for the larger n-pyrenes, more unpaired electrons become concentrated on the zigzag edges of the larger n-pyrenes. Nevertheless, the reason of this correlation (i.e., the connection between edge states and radical formation) remains unclear, and hence, it can be of great interest to develop a simple model to qualitatively explain this correlation in the future.

4. Conclusions

In conclusion, the radical character associated with GNFs has posed a serious challenge to common computational approaches. The challenge has been greatly resolved by the recently formulated TAO-DFT, which is a highly efficient computational method for exploring the electronic properties associated with nanosystems exhibiting radical nature. Here, we have adopted TAO-DFT to explore the electronic properties associated with n-pyrenes with n = 2–15. The larger ground-state n-pyrenes have been found to be polyradicals, playing a significant role in the determination of their electronic properties. Because of the polyradical character associated with the larger n-pyrenes, studying their electronic properties using KS-DFT with the commonly employed XC energy functionals may not be reliable, and studying their electronic properties using accurate MR electronic structure methods may not be computationally feasible. Accordingly, studying the electronic properties of n-pyrenes using TAO-DFT in the present work is definitely justified.
From our TAO-DFT findings, n-pyrenes are ground-state singlets for n = 2–15 (i.e., for all the n values considered). With increasing the size of n-pyrene, the value of E ST , the value of IP v , and the value of E g monotonically decrease, while the value of EA v and the value of S vN monotonically increase. When n increases, there is a smooth transition from the nonradical nature of the smaller n-pyrenes (e.g., n 3 ) to the increasing polyradical nature of the larger n-pyrenes (e.g., n > 3 ). Furthermore, the latter is expected to be intimately correlated with the increasing concentration of active orbitals on the zigzag edges of the larger n-pyrenes.

Supplementary Materials

The following are available online at https://www.mdpi.com/2079-4991/10/6/1236/s1, Table S1: Singlet-triplet energy gap E ST (in kcal/mol) of n-pyrene, calculated by spin-unrestricted TAO-LDA, Table S2: Vertical ionization potential IP v (in eV), vertical electron affinity EA v (in eV), fundamental gap E g (in eV), and symmetrized von Neumann entropy S vN of ground-state n-pyrene, calculated by spin-unrestricted TAO-LDA, Table S3: Active orbital occupation numbers (HOMO−9, HOMO−8, …, HOMO, LUMO, …, LUMO+8, and LUMO+9) of ground-state n-pyrene, calculated by spin-restricted TAO-LDA. For brevity, HOMO/LUMO is denoted as H/L.

Author Contributions

Conceptualization, H.-J.H., S.S. and J.-D.C.; Data curation, H.-J.H. and S.S.; Formal analysis, H.-J.H. and S.S.; Funding acquisition, J.-D.C.; Investigation, H.-J.H. and S.S.; Methodology, H.-J.H., S.S. and J.-D.C.; Project administration, J.-D.C.; Resources, J.-D.C.; Software, J.-D.C.; Supervision, J.-D.C.; Validation, H.-J.H., S.S. and J.-D.C.; Visualization, H.-J.H. and S.S.; Writing—original draft, S.S.; Writing—review & editing, J.-D.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Ministry of Science and Technology of Taiwan, grant number MOST107-2628-M-002-005-MY3.

Acknowledgments

We thank the Ministry of Science and Technology of Taiwan for the financial support (MOST107-2628-M-002-005-MY3), National Taiwan University for the financial support (NTU-CDP-105R7818), and the NCTS of Taiwan.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Geim, A.K.; Novoselov, K.S. The rise of graphene. Nat. Mater. 2007, 6, 183–191. [Google Scholar] [CrossRef]
  2. Novoselov, K.S.; Fal’ko, V.I.; Colombo, L.; Gellert, P.R.; Schwab, M.G.; Kim, K. A roadmap for graphene. Nature 2012, 490, 192–200. [Google Scholar] [CrossRef]
  3. Owens, F.J. Electronic and magnetic properties of armchair and zigzag graphene nanoribbons. J. Chem. Phys. 2008, 128, 194701. [Google Scholar] [CrossRef]
  4. Lee, H.; Ihm, J.; Cohen, M.L.; Louie, S.G. Calcium-decorated graphene-based nanostructures for hydrogen storage. Nano Lett. 2010, 10, 793–798. [Google Scholar] [CrossRef] [Green Version]
  5. Wu, C.-S.; Chai, J.-D. Electronic properties of zigzag graphene nanoribbons studied by TAO-DFT. J. Chem. Theory Comput. 2015, 11, 2003–2011. [Google Scholar] [CrossRef]
  6. Yeh, C.-N.; Lee, P.-Y.; Chai, J.-D. Electronic and optical properties of the narrowest armchair graphene nanoribbons studied by density functional methods. Aust. J. Chem. 2016, 69, 960–968. [Google Scholar] [CrossRef] [Green Version]
  7. Mokkath, J.H. Optical response tuning in graphene nanoflakes: A computational study. Chem. Phys. Lett. 2018, 692, 102–105. [Google Scholar] [CrossRef]
  8. Georgakilas, V.; Vrettos, K.; Katomeri, K.; Kouloumpis, A.; Dimos, K.; Gournis, D.; Zboril, R. Highly dispersible disk-like graphene nanoflakes. Nanoscale 2015, 7, 15059–15064. [Google Scholar] [CrossRef]
  9. Nagai, H.; Nakano, M.; Yoneda, K.; Fukui, H.; Minami, T.; Bonness, S.; Kishi, R.; Takahashi, H.; Kubo, T.; Kamada, K.; et al. Theoretical study on third-order nonlinear optical properties in hexagonal graphene nanoflakes: Edge shape effect. Chem. Phys. Lett. 2009, 477, 355–359. [Google Scholar] [CrossRef]
  10. Hu, W.; Lin, L.; Yang, C.; Yang, J. Electronic structure and aromaticity of large-scale hexagonal graphene nanoflakes. J. Chem. Phys. 2014, 141, 214704. [Google Scholar] [CrossRef] [Green Version]
  11. Mansilla Wettstein, C.; Bonafé, F.P.; Oviedo, M.B.; Sánchez, C.G. Optical properties of graphene nanoflakes: Shape matters. J. Chem. Phys. 2016, 144, 224305. [Google Scholar] [CrossRef] [Green Version]
  12. Yeh, C.-N.; Wu, C.; Su, H.; Chai, J.-D. Electronic properties of the coronene series from thermally-assisted-occupation density functional theory. RSC Adv. 2018, 8, 34350–34358. [Google Scholar] [CrossRef] [Green Version]
  13. Hu, W.; Huang, Y.; Qin, X.; Lin, L.; Kan, E.; Li, X.; Yang, C.; Yang, J. Room-temperature magnetism and tunable energy gaps in edge-passivated zigzag graphene quantum dots. npj 2D Mater. Appl. 2019, 3, 17. [Google Scholar] [CrossRef]
  14. Shi, H.; Barnard, A.S.; Snook, I.K. Modelling the role of size, edge structure and terminations on the electronic properties of trigonal graphene nanoflakes. Nanotechnology 2012, 23, 065707. [Google Scholar] [CrossRef]
  15. Deng, Q.; Chai, J.-D. Electronic properties of triangle-shaped graphene nanoflakes from TAO-DFT. ACS Omega 2019, 4, 14202–14210. [Google Scholar] [CrossRef] [Green Version]
  16. Ge, Y.; Ji, J.; Shen, Z.; Zhang, Q.; Jian, A.; Duan, Q.; Wang, C.; Jiang, J.; Zhang, W.; Sang, S. First principles study of magnetism induced by topological frustration of bowtie-shaped graphene nanoflake. Carbon 2018, 127, 432–436. [Google Scholar] [CrossRef] [Green Version]
  17. Wang, W.L.; Meng, S.; Kaxiras, E. Graphene nanoflakes with large spin. Nano Lett. 2008, 8, 241–245. [Google Scholar] [CrossRef]
  18. Kuc, A.; Heine, T.; Seifert, G. Structural and electronic properties of graphene nanoflakes. Phys. Rev. B 2010, 81, 085430. [Google Scholar] [CrossRef] [Green Version]
  19. Zhou, J.; Wang, Q.; Sun, Q.; Jena, P. Intrinsic ferromagnetism in two-dimensional carbon structures: Triangular graphene nanoflakes linked by carbon chains. Phys. Rev. B 2011, 84, 081402(R). [Google Scholar] [CrossRef] [Green Version]
  20. Singh, S.K.; Neek-Amal, M.; Peeters, F.M. Electronic properties of graphene nano-flakes: Energy gap, permanent dipole, termination effect, and Raman spectroscopy. J. Chem. Phys. 2014, 140, 074304. [Google Scholar] [CrossRef] [Green Version]
  21. Wohner, N.; Lam, P.; Sattler, K. Energetic stability of graphene nanoflakes and nanocones. Carbon 2014, 67, 721–735. [Google Scholar] [CrossRef]
  22. Fernandez, M.; Shi, H.; Barnard, A.S. Geometrical features can predict electronic properties of graphene nanoflakes. Carbon 2016, 103, 142–150. [Google Scholar] [CrossRef]
  23. Wang, W.L.; Yazyev, O.V.; Meng, S.; Kaxiras, E. Topological frustration in graphene nanoflakes: Magnetic order and spin logic devices. Phys. Rev. Lett. 2009, 102, 157201. [Google Scholar] [CrossRef] [PubMed]
  24. Trauzettel, B.; Bulaev, D.V.; Loss, D.; Burkard, G. Spin qubits in graphene quantum dots. Nat. Phys. 2007, 3, 192–196. [Google Scholar] [CrossRef] [Green Version]
  25. Droth, M.; Burkard, G. Spintronics with graphene quantum dots. Phys. Status Solidi RRL 2016, 10, 75–90. [Google Scholar] [CrossRef]
  26. Pavliček, N.; Mistry, A.; Majzik, Z.; Moll, N.; Meyer, G.; Fox, D.J.; Gross, L. Synthesis and characterization of triangulene. Nat. Nanotechnol. 2017, 12, 308–311. [Google Scholar] [CrossRef]
  27. Zhao, L.; Kaiser, R.I.; Xu, B.; Ablikim, U.; Ahmed, M.; Joshi, D.; Veber, G.; Fischer, F.R.; Mebel, A.M. Pyrene synthesis in circumstellar envelopes and its role in the formation of 2D nanostructures. Nat. Astron. 2018, 2, 413–419. [Google Scholar] [CrossRef]
  28. Förster, T.; Kasper, K. Ein Konzentrationsumschlag der Fluoreszenz. Z. Physik. Chem. (Frankfurt) 1954, 1, 275–277. [Google Scholar] [CrossRef]
  29. Brus, L. Size, dimensionality, and strong electron correlation in nanoscience. Acc. Chem. Res. 2014, 47, 2951–2959. [Google Scholar] [CrossRef]
  30. Kohn, W.; Sham, L.J. Self-consistent equations including exchange and correlation effects. Phys. Rev. 1965, 140, A1133–A1138. [Google Scholar] [CrossRef] [Green Version]
  31. Cohen, A.J.; Mori-Sánchez, P.; Yang, W. Insights into current limitations of density functional theory. Science 2008, 321, 792–794. [Google Scholar] [CrossRef] [Green Version]
  32. Cohen, A.J.; Mori-Sánchez, P.; Yang, W. Challenges for density functional theory. Chem. Rev. 2012, 112, 289–320. [Google Scholar] [CrossRef] [PubMed]
  33. Andersson, K.; Malmqvist, P.-Å.; Roos, B.O. Second-order perturbation theory with a complete active space self-consistent field reference function. J. Chem. Phys. 1992, 96, 1218. [Google Scholar] [CrossRef]
  34. Hachmann, J.; Dorando, J.J.; Aviles, M.; Chan, G.K.L. The radical character of the acenes: A density matrix renormalization group study. J. Chem. Phys. 2007, 127, 134309. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Gidofalvi, G.; Mazziotti, D.A. Active-space two-electron reduced-density-matrix method: Complete active-space calculations without diagonalization of the N-electron hamiltonian. J. Chem. Phys. 2008, 129, 134108. [Google Scholar] [CrossRef]
  36. Gryn’ova, G.; Coote, M.L.; Corminboeuf, C. Theory and practice of uncommon molecular electronic configurations. WIREs Comput. Mol. Sci. 2015, 5, 440–459. [Google Scholar] [CrossRef] [Green Version]
  37. Fosso-Tande, J.; Nguyen, T.-S.; Gidofalvi, G.; DePrince, A.E., III. Large-scale variational two-electron reduced-density-matrix-driven complete active space self-consistent field methods. J. Chem. Theory Comput. 2016, 12, 2260–2271. [Google Scholar] [CrossRef]
  38. Battaglia, S.; Faginas-Lago, N.; Andrae, D.; Evangelisti, S.; Leininger, T. Increasing radical character of large [n]cyclacenes unveiled by wave function theory. J. Phys. Chem. A 2017, 121, 3746–3756. [Google Scholar] [CrossRef]
  39. Mullinax, J.W.; Maradzike, E.; Koulias, L.N.; Mostafanejad, M.; Epifanovsky, E.; Gidofalvi, G.; DePrince, A.E., III. Heterogeneous CPU + GPU algorithm for variational two-electron reduced-density matrix-driven complete active-space self-consistent field theory. J. Chem. Theory Comput. 2019, 15, 6164–6178. [Google Scholar] [CrossRef]
  40. Chai, J.-D. Density functional theory with fractional orbital occupations. J. Chem. Phys. 2012, 136, 154104. [Google Scholar] [CrossRef]
  41. Chai, J.-D. Thermally-assisted-occupation density functional theory with generalized-gradient approximations. J. Chem. Phys. 2014, 140, 18A521. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Chai, J.-D. Role of exact exchange in thermally-assisted-occupation density functional theory: A proposal of new hybrid schemes. J. Chem. Phys. 2017, 146, 044102. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Xuan, F.; Chai, J.-D.; Su, H. Local density approximation for the short-range exchange free energy functional. ACS Omega 2019, 4, 7675–7683. [Google Scholar] [CrossRef] [PubMed]
  44. Lin, C.-Y.; Hui, K.; Chung, J.-H.; Chai, J.-D. Self-consistent determination of the fictitious temperature in thermally-assisted-occupation density functional theory. RSC Adv. 2017, 7, 50496–50507. [Google Scholar] [CrossRef] [Green Version]
  45. Yeh, C.-N.; Chai, J.-D. Role of Kekulé and non-Kekulé structures in the radical character of alternant polycyclic aromatic hydrocarbons: A TAO-DFT study. Sci. Rep. 2016, 6, 30562. [Google Scholar] [CrossRef]
  46. Seenithurai, S.; Chai, J.-D. Effect of Li adsorption on the electronic and hydrogen storage properties of acenes: A dispersion-corrected TAO-DFT study. Sci. Rep. 2016, 6, 33081. [Google Scholar] [CrossRef] [Green Version]
  47. Wu, C.-S.; Lee, P.-Y.; Chai, J.-D. Electronic properties of cyclacenes from TAO-DFT. Sci. Rep. 2016, 6, 37249. [Google Scholar] [CrossRef]
  48. Seenithurai, S.; Chai, J.-D. Effect of Li termination on the electronic and hydrogen storage properties of linear carbon chains: A TAO-DFT study. Sci. Rep. 2017, 7, 4966. [Google Scholar] [CrossRef] [Green Version]
  49. Seenithurai, S.; Chai, J.-D. Electronic and hydrogen storage properties of Li-terminated linear boron chains studied by TAO-DFT. Sci. Rep. 2018, 8, 13538. [Google Scholar] [CrossRef]
  50. Chung, J.-H.; Chai, J.-D. Electronic properties of Möbius cyclacenes studied by thermally-assisted-occupation density functional theory. Sci. Rep. 2019, 9, 2907. [Google Scholar] [CrossRef] [Green Version]
  51. Seenithurai, S.; Chai, J.-D. Electronic properties of linear and cyclic boron nanoribbons from thermally-assisted-occupation density functional theory. Sci. Rep. 2019, 9, 12139. [Google Scholar] [CrossRef] [PubMed]
  52. Hanson-Heine, M.W.D.; Rogers, D.M.; Woodward, S.; Hirst, J.D. Dewar benzenoids discovered in carbon nanobelts. J. Phys. Chem. Lett. 2020, 11, 3769–3772. [Google Scholar] [CrossRef] [PubMed]
  53. Manassir, M.; Pakiari, A.H. Total non-Lewis structures: An application to predict the stability and reactivity of linear and angular polyacenes. J. Mol. Graph. Model. 2020. [Google Scholar] [CrossRef]
  54. Hanson-Heine, M.W.D. Static correlation in vibrational frequencies studied using thermally-assisted-occupation density functional theory. Chem. Phys. Lett. 2020, 739, 137012. [Google Scholar] [CrossRef]
  55. Shao, Y.; Gan, Z.; Epifanovsky, E.; Gilbert, A.T.; Wormit, M.; Kussmann, J.; Lange, A.W.; Behn, A.; Deng, J.; Feng, X.; et al. Advances in molecular quantum chemistry contained in the Q-Chem 4 program package. Mol. Phys. 2015, 113, 184–215. [Google Scholar] [CrossRef] [Green Version]
  56. Su, Y.; Wang, X.; Wang, L.; Zhang, Z.; Wang, X.; Song, Y.; Power, P.P. Thermally controlling the singlet-triplet energy gap of a diradical in the solid state. Chem. Sci. 2016, 7, 6514–6518. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  57. Yu, L.; Wu, Z.; Xie, G.; Zhong, C.; Zhu, Z.; Cong, H.; Ma, D.; Yang, C. Achieving a balance between small singlet-triplet energy splitting and high fluorescence radiative rate in a quinoxaline-based orange-red thermally activated delayed fluorescence emitter. Chem. Commun. 2016, 52, 11012–11015. [Google Scholar] [CrossRef] [PubMed]
  58. Rajamalli, P.; Senthilkumar, N.; Gandeepan, P.; Ren-Wu, C.-C.; Lin, H.W.; Cheng, C.H. A method for reducing the singlet-triplet energy gaps of TADF materials for improving the blue OLED efficiency. ACS Appl. Mater. Interfaces 2016, 8, 27026–27034. [Google Scholar] [CrossRef]
  59. Rajamalli, P.; Senthilkumar, N.; Gandeepan, P.; Huang, P.-Y.; Huang, M.-J.; Ren-Wu, C.-Z.; Yang, C.-Y.; Chiu, M.-J.; Chu, L.-K.; Lin, H.-W.; et al. A new molecular design based on thermally activated delayed fluorescence for highly efficient organic light emitting diodes. J. Am. Chem. Soc. 2016, 138, 628–634. [Google Scholar] [CrossRef]
  60. Rivero, P.; Jiménez-Hoyos, C.A.; Scuseria, G.E. Entanglement and polyradical character of polycyclic aromatic hydrocarbons predicted by projected Hartree-Fock theory. J. Phys. Chem. B 2013, 117, 12750–12758. [Google Scholar] [CrossRef]
Figure 1. Structure of n-pyrene, which contains n benzenoid rings fused together at each side.
Figure 1. Structure of n-pyrene, which contains n benzenoid rings fused together at each side.
Nanomaterials 10 01236 g001
Figure 2. Singlet-triplet energy gap of n-pyrene ((a) n = 2–15 and (b) n = 4–15), calculated by spin-unrestricted TAO-LDA.
Figure 2. Singlet-triplet energy gap of n-pyrene ((a) n = 2–15 and (b) n = 4–15), calculated by spin-unrestricted TAO-LDA.
Nanomaterials 10 01236 g002
Figure 3. Vertical ionization potential of ground-state n-pyrene, calculated by spin-unrestricted TAO-LDA.
Figure 3. Vertical ionization potential of ground-state n-pyrene, calculated by spin-unrestricted TAO-LDA.
Nanomaterials 10 01236 g003
Figure 4. Vertical electron affinity of ground-state n-pyrene, calculated by spin-unrestricted TAO-LDA.
Figure 4. Vertical electron affinity of ground-state n-pyrene, calculated by spin-unrestricted TAO-LDA.
Nanomaterials 10 01236 g004
Figure 5. Fundamental gap of ground-state n-pyrene, calculated by spin-unrestricted TAO-LDA.
Figure 5. Fundamental gap of ground-state n-pyrene, calculated by spin-unrestricted TAO-LDA.
Nanomaterials 10 01236 g005
Figure 6. Symmetrized von Neumann entropy of ground-state n-pyrene, calculated by spin-unrestricted TAO-LDA.
Figure 6. Symmetrized von Neumann entropy of ground-state n-pyrene, calculated by spin-unrestricted TAO-LDA.
Nanomaterials 10 01236 g006
Figure 7. Active orbital occupation numbers (H−9, H−8, …, H, L, …, L+8, and L+9) of ground-state n-pyrene, calculated by spin-restricted TAO-LDA. HOMO/LUMO is denoted as H/L for brevity.
Figure 7. Active orbital occupation numbers (H−9, H−8, …, H, L, …, L+8, and L+9) of ground-state n-pyrene, calculated by spin-restricted TAO-LDA. HOMO/LUMO is denoted as H/L for brevity.
Nanomaterials 10 01236 g007
Figure 8. Real-space representation of H−2 (2.000), H−1 (2.000), H (1.998), L (0.002), L+1 (0.000), and L+2 (0.000) of ground-state 2-pyrene, calculated by spin-restricted TAO-LDA. Here the isovalue is 0.02 e/Å 3 . The occupation numbers of active orbitals are given in parentheses, and HOMO/LUMO is denoted as H/L for brevity.
Figure 8. Real-space representation of H−2 (2.000), H−1 (2.000), H (1.998), L (0.002), L+1 (0.000), and L+2 (0.000) of ground-state 2-pyrene, calculated by spin-restricted TAO-LDA. Here the isovalue is 0.02 e/Å 3 . The occupation numbers of active orbitals are given in parentheses, and HOMO/LUMO is denoted as H/L for brevity.
Nanomaterials 10 01236 g008
Figure 9. Real-space representation of H−2 (2.000), H−1 (1.976), H (1.669), L (0.334), L+1 (0.021), and L+2 (0.000) of ground-state 4-pyrene, calculated by spin-restricted TAO-LDA. Here the isovalue is 0.02 e/Å 3 . The occupation numbers of active orbitals are given in parentheses, and HOMO/LUMO is denoted as H/L for brevity.
Figure 9. Real-space representation of H−2 (2.000), H−1 (1.976), H (1.669), L (0.334), L+1 (0.021), and L+2 (0.000) of ground-state 4-pyrene, calculated by spin-restricted TAO-LDA. Here the isovalue is 0.02 e/Å 3 . The occupation numbers of active orbitals are given in parentheses, and HOMO/LUMO is denoted as H/L for brevity.
Nanomaterials 10 01236 g009
Figure 10. Real-space representation of H−2 (1.905), H−1 (1.524), H (1.161), L (0.852), L+1 (0.478), and L+2 (0.080) of ground-state 6-pyrene, calculated by spin-restricted TAO-LDA. Here the isovalue is 0.02 e/Å 3 . The occupation numbers of active orbitals are given in parentheses, and HOMO/LUMO is denoted as H/L for brevity.
Figure 10. Real-space representation of H−2 (1.905), H−1 (1.524), H (1.161), L (0.852), L+1 (0.478), and L+2 (0.080) of ground-state 6-pyrene, calculated by spin-restricted TAO-LDA. Here the isovalue is 0.02 e/Å 3 . The occupation numbers of active orbitals are given in parentheses, and HOMO/LUMO is denoted as H/L for brevity.
Nanomaterials 10 01236 g010
Figure 11. Real-space representation of H−2 (1.404), H−1 (1.141), H (1.028), L (0.992), L+1 (0.895), and L+2 (0.576) of ground-state 8-pyrene, calculated by spin-restricted TAO-LDA. Here the isovalue is 0.02 e/Å 3 . The occupation numbers of active orbitals are given in parentheses, and HOMO/LUMO is denoted as H/L for brevity.
Figure 11. Real-space representation of H−2 (1.404), H−1 (1.141), H (1.028), L (0.992), L+1 (0.895), and L+2 (0.576) of ground-state 8-pyrene, calculated by spin-restricted TAO-LDA. Here the isovalue is 0.02 e/Å 3 . The occupation numbers of active orbitals are given in parentheses, and HOMO/LUMO is denoted as H/L for brevity.
Nanomaterials 10 01236 g011
Figure 12. Real-space representation of H−2 (1.133), H−1 (1.040), H (1.038), L (1.011), L+1 (0.981), and L+2 (0.872) of ground-state 10-pyrene, calculated by spin-restricted TAO-LDA. Here the isovalue is 0.02 e/Å 3 . The occupation numbers of active orbitals are given in parentheses, and HOMO/LUMO is denoted as H/L for brevity.
Figure 12. Real-space representation of H−2 (1.133), H−1 (1.040), H (1.038), L (1.011), L+1 (0.981), and L+2 (0.872) of ground-state 10-pyrene, calculated by spin-restricted TAO-LDA. Here the isovalue is 0.02 e/Å 3 . The occupation numbers of active orbitals are given in parentheses, and HOMO/LUMO is denoted as H/L for brevity.
Nanomaterials 10 01236 g012

Share and Cite

MDPI and ACS Style

Huang, H.-J.; Seenithurai, S.; Chai, J.-D. TAO-DFT Study on the Electronic Properties of Diamond-Shaped Graphene Nanoflakes. Nanomaterials 2020, 10, 1236. https://doi.org/10.3390/nano10061236

AMA Style

Huang H-J, Seenithurai S, Chai J-D. TAO-DFT Study on the Electronic Properties of Diamond-Shaped Graphene Nanoflakes. Nanomaterials. 2020; 10(6):1236. https://doi.org/10.3390/nano10061236

Chicago/Turabian Style

Huang, Hong-Jui, Sonai Seenithurai, and Jeng-Da Chai. 2020. "TAO-DFT Study on the Electronic Properties of Diamond-Shaped Graphene Nanoflakes" Nanomaterials 10, no. 6: 1236. https://doi.org/10.3390/nano10061236

APA Style

Huang, H. -J., Seenithurai, S., & Chai, J. -D. (2020). TAO-DFT Study on the Electronic Properties of Diamond-Shaped Graphene Nanoflakes. Nanomaterials, 10(6), 1236. https://doi.org/10.3390/nano10061236

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