Study of Li Adsorption on Graphdiyne Using Hybrid DFT Calculations

May 3, 2018 - adsorption energy comparable to that of a periodic system. Our ... the atom-center spheres with the radius of 5.5 Å determined by the...
546 downloads 3 Views 2MB Size
Forum Article Cite This: ACS Appl. Mater. Interfaces XXXX, XXX, XXX−XXX

www.acsami.org

Study of Li Adsorption on Graphdiyne Using Hybrid DFT Calculations Jaewook Kim, Sungwoo Kang, Jaechang Lim, and Woo Youn Kim* Department of Chemistry, Korea Advanced Institute of Science and Technology (KAIST), 291 Daehak-ro, Yuseong-gu, Daejeon 34141, Republic of Korea ABSTRACT: Promising applications of graphdiyne have often been initiated by theoretical predictions especially using DFT known as the most powerful first-principles electronic structure calculation method. However, there is no systematic study on the reliability of DFT for the prediction of the electronic properties of the graphdiyne. Here, we performed a study of Li adsorption on the graphdiyne using hybrid DFT with LC-ωPBE and compared the results with those of PBE, because accurate prediction of the Li adsorption is important for performance as a Li storage that was first theoretically suggested and then experimentally realized. Our results show that PBE overestimates the adsorption energy inside a pore and the barrier height at the transition state of in-plane diffusion compared to the those of LCωPBE. In particular, LC-ωPBE predicted almost barrier-less in-plane diffusion of Li on the graphdiyne because of the presence of both in-plane and out-of-plane π orbitals. Also, LC-ωPBE favors a high spin state due to the exact exchange energy when several Li atoms are adsorbed on the graphdiyne, whereas PBE favors a low spin state. Thus, the use of the hybrid DFT is critical for reliable predictions on the electronic properties of the graphdiyne. KEYWORDS: graphdiyne, lithium storage, energy storage, density functional theory, hybrid functional

1. INTRODUCTION Graphdiyne (GDY) has attracted much attention as a new 2D carbon allotrope consisting of sp- and sp2-hybridized carbons, since synthesis of large-scale GDYs became possible.1,2 It shows excellent conductivity and mechanical properties like graphene, while at the same time exhibiting semiconducting properties due to large pores formed with 6-alkyne bonds in contrast to the semimetallic graphene.3−9 For these reasons, the GDY has been actively studied for various applications such as energy storage,10−16 a catalyst,17−20 an optoelectronic device,21−23 and a permeable membrane.24−26 Along the recent trend in new materials research, application of the GDY has often been initiated by computational predictions and then realized in experiments. For example, the active study of the GDY for optoelectronic devices has been stimulated by the prediction of the natural band gap and high hole mobility of the GDY.27,28 Also, the feasibility of the GDY as a Li storage material has been proposed by density functional theory (DFT) calculations in 2012 and 2013 and proven by experiments in 2015.10−13 Therefore, accurate prediction of the GDY properties using valid computational methods is vital for successful applications. The GDY properties can be predicted by DFT which has been most widely used for electronic structure calculations because of its reasonable reliability and efficiency.29 Most previous DFT studies on the GDY under periodic boundary conditions adopted the generalized gradient approximation (GGA) such as PBE30 for the exchange-correlation functional.5,8,10,11,16,24,27,28,31−34 However, there is no systematic study on the reliability of the GGA method. In general, hybrid DFT methods made by including a fractional amount of exact exchange energy in GGA show better performance in many © XXXX American Chemical Society

aspects than GGA. A benchmark study on thermodynamic properties shows that the hybrid DFT is far more accurate than PBE; for example, the mean absolute errors (MAEs) of LCωPBE35 and B3LYP36 for the atomization energies of the G2 test set37 are 5.6 and 7.5 kcal/mol, respectively, whereas that of PBE is 17.9 kcal/mol. For the activation energies in the HTBH38/0438 and NHTBH38/0439 test sets, LC-ωPBE shows much smaller MAE (2.4 kcal/mol) than that of PBE (MAE = 10.7 kcal/mol).40,41 Given that a reaction rate varies exponentially with the activation energy, use of the hybrid DFT for activation energies is more critical than for adsorption energies. In particular, magnetic properties and band gaps are even more sensitive to the choice of appropriate hybrid functionals whose accuracy can be tuned by the amount of exact exchange.42 In fact, a few previous works used hybrid DFT methods for adsorption energy calculations, but they employed a small flake of GDY.19,20,43−47 That is because a computational load of the hybrid DFT is extremely high under periodic boundary conditions. Alternatively, namely the GGA+U48 method devised for strong on-site interactions can be utilized for periodic systems including transition metals, since it shows similar results with those of hybrid functionals.49,50 However, it Special Issue: Graphdiyne Materials: Preparation, Structure, and Function Received: March 1, 2018 Accepted: May 3, 2018

A

DOI: 10.1021/acsami.8b03482 ACS Appl. Mater. Interfaces XXXX, XXX, XXX−XXX

Forum Article

ACS Applied Materials & Interfaces

where EGDY, ELi, atom, and Etotal denote the energies of the GDY flake, a Li atom, and the adsorbed system, respectively. To make sure that Eads of the flake is quantitatively comparable to that of the corresponding bulk system, we investigated the convergence of Eads for a single Li atom as increasing the flake size from a single triangular unit (GDY-1tri) to a seven triangular units (GDY-7tri) as depicted in Figure 1. Eads of PBE

would not be appropriate for GDY systems with no transition metals. In this work, we investigated the effect of the hybrid DFT on the prediction of the electronic properties of GDY. Specifically, we studied the adsorption energy of Li on the GDY as a candidate for a Li storage, which is one of the most promising applications of the GDY. We compared the results of LC-ωPBE with those of PBE. LC-ωPBE as a well-known long-range corrected hybrid functional is known to provide more accurate thermodynamic and activation energies than those of B3LYP as a de facto standard hybrid method in chemistry. In addition, accelerated hybrid calculations are feasible with LC-ωPBE by employing an approximation known as the 2-gau method.41 We proposed a new approach of the hybrid DFT based on a local exact exchange potential method and showed that LC-ωPBE modified by our new approach together with the 2-gau method is just about twice slower than PBE in a grid-based DFT code, enabling us to apply it to large systems.51 Thus, we were able to employ a sufficiently large GDY flake to give converged adsorption energy comparable to that of a periodic system. Our study shows that LC-ωPBE and PBE give considerably different adsorption properties of Li on the GDY, implying that the use of the hybrid DFT is critical for reliable predictions on the electronic properties of the GDY.

2. COMPUTATIONAL METHODS To determine an appropriate size of the GDY flake, we checked out the Li adsorption energy convergence with respect to the flake size. The flake geometry should be the same as that of the periodic system for comparison to the bulk value. The bulk geometries of the GDY with and without Li were obtained from periodic calculations with PBE augmented by the D3 dispersion correction52 with the BeckeJohnson damping53 as implemented in VASP 5.4;54 as already known, structural change by the D3 correction was negligible.10 We used the energy cutoff of 400 eV and the k-points of 7 × 7 × 1. In both GDY and Li-adsorbed GDY (GDY+Li), the unit cell was set to a = b = 9.465 Å, c = 20 Å, α = β = 90°, and γ = 60°. Then, the optimized GDY structure was cut into a flake, and the resulting dangling bonds were capped with hydrogen atoms. The flake size was increased by repeating the unit cell geometry. The actual adsorption energy on the GDY flake was calculated using our in-house grid-based code, ACE-Molecule,51,55−58 for direct comparison with the results of VASP, because ACE-Molecule also employs the projector augmented wave (PAW)59 method for core electrons like VASP. In addition, it supports an accelerated LC-ωPBE method for the efficient hybrid DFT calculation of large GDY flakes. To achieve the accuracy of 0.1 eV (∼3 kcal/mol) for the adsorption energy, the grid spacing was set to 0.45 Bohr, which corresponds to the kinetic energy cutoff of ∼663 eV in a plane-wave code like VASP. The simulation box was made by superimposing all of the atom-center spheres with the radius of 5.5 Å determined by the total energy convergence test. For diffusion study, the transition state of the in-plane diffusion of Li was obtained by using the TS optimization method as implemented in GAUSSIAN0960 with PBE/631+G(d). Then, the activation energy on the given geometry was calculated by ACE-Molecule using the same simulation setting with that of the adsorption energy. All orbital figures were drawn with the isovalue of 0.02.

Figure 1. Structure of the Li-adsorbed GDY flake with various sizes and the corresponding adsorption energies with PBE/LC-ωPBE. Each structure corresponds to (a) single triangular unit, (b) two triangular units, (c) four triangular units, and (d) seven triangular units.

almost converged at the size of GDY-4tri, whereas LC-ωPBE converged at the size of GDY-7tri. In the case of PBE, we also performed a periodic calculation on the gamma point for a single Li adsorbed supercell of 3 × 3 × 1 with the plane-wave cutoff of 400 eV. For comparison between the bulk and the flake, we imposed the doublet spin state for the periodic system. The Li adsorption energy of the bulk system was 2.24 eV, which is somewhat higher than the converged value of the flake (1.82 eV). We noted that the adsorption energy of GDY2tri from GAUSSIAN09 with PBE/cc-pVTZ (aug-cc-pVTZ for Li) (1.68 eV) was even lower than that of ACE-Molecule. These differences are probably due to the difference between the methods used for core electrons in each package; GAUSSIAN09 treated all electrons explicitly, ACE-Molecule used the PAW setup files provided by GPAW,61 while VASP adopted its own files. We note that for all the flake sizes, the Eads of LC-ωPBE was significantly lower than that of PBE. To understand such a tendency, we examined the electronic mechanism of Li adsorption to the GDY flake. We analyzed the orbitals of the Li atom, GDY-1tri, and GDY-1tri+Li. The 2s orbital energy of Li was 2 eV higher than the highest occupied molecular orbital (HOMO) energy of GDY-1tri for both PBE and LC-ωPBE. In contrast, the lowest unoccupied molecular orbital (LUMO) energy of GDY-1tri was very close to the 2s orbital energy of Li. Therefore, the 2s orbital of Li forms a bond with the LUMO of GDY by sharing the single valence electron of Li, resulting in the electron transfer from Li to GDY. The LUMO of GDY is

3. RESULTS AND DISCUSSION 3.1. Li Adsorption Energy. We first compared the adsorption energies of Li on the GDY flake obtained with PBE and LC-ωPBE. The adsorption energy (Eads) was calculated as follows. Eads = −(Etotal − EGDY − E Li,atom) B

DOI: 10.1021/acsami.8b03482 ACS Appl. Mater. Interfaces XXXX, XXX, XXX−XXX

Forum Article

ACS Applied Materials & Interfaces

Figure 2. Orbital energy diagram of the Li-adsorbed GDY flake and the corresponding orbitals. All results were obtained with PBE.

made of π orbitals perpendicular to the GDY plane, and hence it is delocalized across the GDY flake. The resulting singly occupied molecular orbital (SOMO) of the GDY+Li maintained the original shape of the LUMO of the free GDY flake. Unlike graphene, the GDY has additional π-bonding orbitals parallel to the GDY plane which comes from acetylenic chains. Li after giving its electron to the LUMO of the GDY attracts electrons in the in-plane π-bonding orbitals, leading to substantial distortion of those orbitals as shown in Figure 2. As a result, the energy levels of the corresponding orbitals dropped down by a few eV (Figure 2). As discussed above, the Li adsorption caused the electron transfer from the localized 2s orbital of Li to the delocalized out-of-plane π orbital of the GDY flake. This explains the reason why PBE overestimates the adsorption energy compared to that of LC-ωPBE. In general, PBE favors a delocalized state of the electron density due to the so-called self-interaction error. Such a spurious repulsion becomes weaker as the electron moves from the localized 2s to the delocalized π orbital, resulting in the overestimation of Eads. In contrast, the exact exchange energy included in hybrid functionals like LCωPBE alleviates the spurious self-interaction error, giving rise to the smaller adsorption energy than that of PBE. Our analysis provides an orbital-based insight about the most preferred Li adsorption site; the Li cation is strongly stabilized by the electron-rich in-plane π orbitals of two adjacent acetylenic chains simultaneously, whereas Li on the out-ofplane is stabilized by only one π orbital. The previous calculation with PBE under a periodic boundary condition showed that the in-plane position is the most preferred adsorption site of Li.10,11 Unfortunately, one cannot perform periodic calculations with LC-ωPBE at present. Instead, we expect that PBE and LC-ωPBE will give similar electronic structures including the π orbitals except for energy gaps, considering the fact that various calculations using LDA, PBE, HSE,62 and G0W063 for a pristine GDY gave similar electronic structures but with different band gaps.3,6,9 Therefore, the adsorption mechanism due to the π orbitals and thus the energetically preferred adsorption site would be the same for PBE and LC-ωPBE. 3.2. Effect of Defects on the Li Adsorption. For the practical use of GDY as a Li storage, it is important to consider the effect of various defects. In fact, the GDY is known to have a low concentration of defects.2 A previous report shows that the electronic structure of the GDY does not change significantly due to carbon vacancies except for the band gap change from 0.43 to 0.37 eV.8 Since the Li adsorption is

affected by the LUMO energy level, it is expected that the carbon vacancies may increase the adsorption energy a little bit due to the smaller band gap. The in-plane diffusion barrier will not change much, because both the in-plane and out-of-plane π orbitals remain intact in the presence of the vacancies. Therefore, overall behavior would be similar in the presence of a small number of the vacancies. However, an aging process of the GDY naturally increases oxygenic and nitrogenic defects.32 The chemisorption of such atoms needs sp2 or sp3 hybridized carbon atoms, giving rise to breaking up the sp hybridization of the GDY frontier orbitals. As a result, it is expected that the adsorption mechanism of Li would be different around such defect sites. 3.3. Transition State Energy of in-Plane Diffusion of Li. To investigate the effect of the hybrid DFT on the in-plane diffusion rate of Li for energy storage application, we compared the energy barrier of LC-ωPBE at the transition state along an in-plane diffusion path across an acetylenic chain with that of PBE. The GDY-2tri was used as a model system. To mimic a periodic calculation using the nudged elastic band method in the literature,10 we changed only the angle between Li and the in-plane vector while maintaining the other geometric factors as shown in Figure 3(a). Then, we investigated the change in the Li adsorption energy according to the angle (θ). The barrier height of PBE (0.53 eV) was close to that of periodic

Figure 3. (a) Structure of the GDY-2tri+Li for in-plane diffusion study. (b) Out-of-plane and in-plane π MOs of the free GDY flake. (c) Change in the out-of-plane and in-plane π MOs of the GDY-2tri+Li with respect to the Li position. C

DOI: 10.1021/acsami.8b03482 ACS Appl. Mater. Interfaces XXXX, XXX, XXX−XXX

Forum Article

ACS Applied Materials & Interfaces calculations with PBE (0.75 eV)10 and with B8864+LYP65 (0.52 eV).11 However, LC-ωPBE gave considerably lower barrier height (0.33 eV). Like in the stable adsorption case, the electron in the 2s orbital of Li was transferred to the delocalized out-of-plane π orbital of the GDY, resulting in the SOMO having a similar shape with the LUMO of the free GDY flake (Figure 3c). The lower barrier height of LC-ωPBE can be explained for the same reason as the lower adsorption energy of LC-ωPBE. To further elucidate the origin of such a low energy barrier, we examined the change in the shapes of the SOMO and the π orbital of the GDY-Li complex with respect to the Li position. On the minimum energy position (θ = 0°), the Li ion made by the electron transfer was stabilized by the distortion of the inplane π orbital (the bottom figure in the first column in Figure 3c). On the transition state (θ = 0°), the Li ion is on top of the triple bond of two carbon atoms and thus interacts mainly with the out-of-plane π orbital. In this case, there would be no stabilization effect by the in-plane π orbital (the bottom figure in the rightmost column in Figure 3c). Instead, the out-of-plane π orbital was distorted by the attraction from the Li ion to stabilize the complex. On the position in-between the minimum energy and transition state (θ = 45°), the Li ion interacts with both the in-plane and out-of-plane π orbitals in part, and thus the energy difference from the other two states becomes small. Consequently, the small energy barrier is due to the presence of both the in-plane and out-of-plane π orbitals forming the triple bond between carbons. Therefore, Li can be stabilized by either in-plane or out-of-plane π orbitals at any angle around the triple bond. To examine the local structural relaxation effect on the barrier height at the transition state, we performed the TS optimization while maintaining the GDY framework with fixed positions of all side atoms. Figure 4 shows the result. We used

structural relaxation during the diffusion gave rise to the barrierless in-plane diffusion of Li, which is a unique feature of GDY due to the presence of sp-hybridized orbitals not being in other 2D carbon allotropes such as graphene. 3.4. Spin State. We investigated the effect of the hybrid DFT on the spin state of the Li-adsorbed GDY flake. To consider the fully charged state that has been computationally discovered, we supposed that three Li atoms occupied the 18carbon cavity of GDY. For the sake of computational efficiency, we used the GDY-2tri with six Li atoms. The GDY-2tri+6Li complex can have four possible spin states: singlet, triplet, quintet, and septet. Table 1 shows the relative energy of the Table 1. Relative Energy of the GDY-2tri+6Li Complex at Various Spin States with Respect to the Lowest Energy State (units in eV) singlet triplet quintet septet

PBE

LC-ωPBE

0.26 0.00 0.25 1.18

0.58 0.73 0.00 0.37

complex for each state with respect to the lowest energy state. PBE favors the triplet state, while LC-ωPBE has the quintet as the lowest energy state. It is well-known that hybrid functionals favor high spin states because of the Pauli Exclusion principle implied in the exact exchange energy, whereas pure functionals having no exact exchange term favor low spin states. The results in Table 1 also follow the aforementioned trend. 3.5. Practical aspects. To investigate the effect of the hybrid DFT on GDY-based Li-ion battery performance, we roughly estimated the open circuit voltage (OCV), storage capacity, and Li ion mobility using the results of PBE and LCωPBE. The OCV can be deduced directly from the adsorption energy of Li with the following formula. VOC = (Eads − E Li,coh)/e

where Eads is the adsorption energy per adsorbed Li atom, ELi,coh is the cohesive energy of the bulk Li with the body-centered cubic unit cell, and e is the electric charge of an electron. We used the converged adsorption energies of Li given in Figure 1 and the experimental value of 1.63 eV for ELi,coh.66 The PBE result predicted a positive OCV value (0.19 V), which is consistent to that of the previous periodic calculations (0.1−0.5 V for LixC6 case with x from 0.43 to 11.1816). The reason for the significantly lower OCV value in our calculation is due to the use of the experimental cohesive energy (ELi,coh of PBE = 1.53 eV67) and the GDY flake. Not so surprisingly, LC-ωPBE gave a negative OCV value (−0.86 V) due to the substantially lower adsorption energy than that of PBE. It is known that the OCV should be in the range of 0−2 V to be qualified as an anode material when the bulk Li is used as a reference electrode.16 However, the negative OCV value of LC-ωPBE does not mean that the GDY is not qualified as an anode. The negative value comes from the crude approximation used here. For more accurate estimation of the OCV, one has to consider the entropy contribution of Li to the adsorption free energy. The entropy of Li will be increased from that of the bulk by the Li adsorption. Another important factor is to consider a multilayer GDY. It is expected that the adsorption energy of Li will be larger in a multilayer than in a single layer. In fact, the experiment of Huang et al. showed that in addition

Figure 4. Transition state geometry for the Li in-plane diffusion on the GDY flake: (a) top view, (b) side view along the arrow in panel a, and (c) tilted-side view.

the same geometry at the transition state obtained from PBE for LC-ωPBE calculation. At the transition state, the Li atom pushed underlying carbon atoms in the acetylenic chain. Thus, the angle between the three carbon atoms was changed from 180° to 155.2° as shown in Figure 4c. As a result, the 2p orbitals of the carbon atoms forming the out-of-plane π bond directed along the Li atom. Electrons in those orbitals were attracted by the Li ion left after the electron transfer, giving rise to the considerable stabilization of the transition state. The resulting energy barrier of PBE at the transition state was 0.13 eV, while that of LC-ωPBE was very close to zero. Thus, the D

DOI: 10.1021/acsami.8b03482 ACS Appl. Mater. Interfaces XXXX, XXX, XXX−XXX

ACS Applied Materials & Interfaces



to the surface contribution, the interlayer insertion and extraction of Li in the multilayer GDY are important to make a high capacity Li storage.13 Therefore, the negative OCV value predicted by LC-ωPBE can be likely compensated by both the larger adsorption energy in the multilayer GDY and the positive entropy contribution, which may result in an overall positive value in the range of 0−2 V. Apparently, the OCV value of PBE will further increase by the same reason. The lower adsorption and OCV value of LC-ωPBE will lead to a lower Li storage capacity than predicted by PBE. In fact, Sun et al. estimated the reversible Li capacity as 744 mAh/g through periodic calculations with PBE on a monolayer GDY.10 However, experimental values using GDY films of 10.9 and 22.1 μm thickness were significantly lower as 520 and 285 mAh/g, respectively,13 than the predicted value. Considering the intercalation effect and the entropy contribution, it is suspicious that PBE strongly overestimates the Li storage capacity. Quantitative comparison between theoretical and experimental values may not be possible due to the lack of information on the intercalation effect and entropy contribution. However, the trend of LC-ωPBE that predicts a lower OCV value and a lower storage capacity than those of PBE supports the experimental results. Similarly, the same trend of LC-ωPBE on the barrier height affects the Li ion mobility and diffusion dynamics which are critical for fast charging and discharging rates. The almost barrierless in-plane diffusion of Li predicted by LC-ωPBE implies high battery rate performance as was reported in experiments.13 It should be noted that these arguments are not based on quantitative analysis. Therefore, further works should be followed for more concrete study using hybrid functionals.

Forum Article

AUTHOR INFORMATION

Corresponding Author

*E-mail: [email protected]. ORCID

Jaewook Kim: 0000-0002-7209-1939 Woo Youn Kim: 0000-0001-7152-2111 Notes

The authors declare no competing financial interest.



ACKNOWLEDGMENTS This work was supported by Basic Science Research Programs (NRF-2015R1A1A1A05001480) funded by the Korea government [MSIP] and by KISTI supercomputing center through the strategic support program (KSC-2017-C3-0034).



REFERENCES

(1) Haley, M. M.; Brand, S. C.; Pak, J. J. Carbon Networks Based on Dehydrobenzoannulenes: Synthesis of Graphdiyne Substructures. Angew. Chem., Int. Ed. Engl. 1997, 36, 836−838. (2) Li, G.; Li, Y.; Liu, H.; Guo, Y.; Li, Y.; Zhu, D. Architecture of Graphdiyne Nanoscale Films. Chem. Commun. 2010, 46, 3256−3258. (3) Luo, G.; Qian, X.; Liu, H.; Qin, R.; Zhou, J.; Li, L.; Gao, Z.; Wang, E.; Mei, W.-N.; Lu, J.; Li, Y.; Nagase, S. Quasiparticle Energies and Excitonic Effects of the Two-Dimensional Carbon Allotrope Graphdiyne: Theory and Experiment. Phys. Rev. B: Condens. Matter Mater. Phys. 2011, 84, 075439. (4) Cranford, S. W.; Brommer, D. B.; Buehler, M. J. Extended Graphynes: Simple Scaling Laws for Stiffness, Strength and Fracture. Nanoscale 2012, 4, 7797−7809. (5) Narita, N.; Nagai, S.; Suzuki, S.; Nakao, K. Optimized Geometries and Electronic Structures of Graphyne and Its Family. Phys. Rev. B: Condens. Matter Mater. Phys. 1998, 58, 11009−11014. (6) Kim, H.; Kim, Y.; Kim, J.; Kim, W. Y. Computational Searching for New Stable Graphyne Structures and Their Electronic Properties. Carbon 2016, 98, 404−410. (7) Kang, B.; Moon, J. H.; Lee, J. Y. Size Dependent Electronic Band Structures of β- and γ-Graphyne Nanotubes. RSC Adv. 2015, 5, 80118−80121. (8) Ghorbanzadeh Ahangari, M. Effect of Defect and Temperature on the Mechanical and Electronic Properties of Graphdiyne: A Theoretical Study. Phys. E 2015, 66, 140−147. (9) Pari, S.; Cuéllar, A.; Wong, B. M. Structural and Electronic Properties of Graphdiyne Carbon Nanotubes from Large-Scale DFT Calculations. J. Phys. Chem. C 2016, 120, 18871−18877. (10) Sun, C.; Searles, D. J. Lithium Storage on Graphdiyne Predicted by DFT Calculations. J. Phys. Chem. C 2012, 116, 26222−26226. (11) Zhang, H.; Xia, Y.; Bu, H.; Wang, X.; Zhang, M.; Luo, Y.; Zhao, M. Graphdiyne: A Promising Anode Material for Lithium Ion Batteries with High Capacity and Rate Capability. J. Appl. Phys. 2013, 113, 044309. (12) Zhang, S.; Liu, H.; Huang, C.; Cui, G.; Li, Y. Bulk Graphdiyne Powder Applied for Highly Efficient Lithium Storage. Chem. Commun. 2015, 51, 1834−1837. (13) Huang, C.; Zhang, S.; Liu, H.; Li, Y.; Cui, G.; Li, Y. Graphdiyne for high capacity and long-life lithium storage. Nano Energy 2015, 11, 481−489. (14) Krishnamoorthy, K.; Thangavel, S.; Chelora Veetil, J.; Raju, N.; Venugopal, G.; Kim, S. J. Graphdiyne Nanostructures as a New Electrode Material for Electrochemical Supercapacitors. Int. J. Hydrogen Energy 2016, 41, 1672−1678. (15) Zhang, S.; Du, H.; He, J.; Huang, C.; Liu, H.; Cui, G.; Li, Y. Nitrogen-Doped Graphdiyne Applied for Lithium-Ion Storage. ACS Appl. Mater. Interfaces 2016, 8, 8467−8473. (16) Jang, B.; Koo, J.; Park, M.; Lee, H.; Nam, J.; Kwon, Y.; Lee, H. Graphdiyne as a high-capacity lithium ion battery anode material. Appl. Phys. Lett. 2013, 103, 263904.

4. CONCLUSION We investigated the effect of hybrid DFT on the Li adsorption and diffusion on graphdiyne (GDY) through a comparative study with PBE known as a standard approach to material research. The results show that LC-ωPBE predicts substantially smaller adsorption energy and almost zero barrier height for inplane diffusion. The barrier-less diffusion was due to the presence of both in-plane and out-of-plane π orbitals forming the triple bond between carbons in an acetylenic chain. The Li cation on the GDY can be stabilized by strong interaction with either in-plane or out-of-plane π orbitals at any position along the diffusion path, resulting in little energy change. That is a unique feature of 2D carbon allotropes with acetylenic chains. When several Li atoms occupied the GDY pores, LC-ωPBE favored a higher spin state because the exact exchange interaction in the hybrid functional stabilizes parallel spins. Our study shows that the use of hybrid DFT is vital for reliable predictions on the electronic properties of new materials such as GDY. Nonetheless, it is seldom used in the computational study of bulk materials, since hybrid DFT calculations are very demanding for periodic systems. There have been continuous efforts for the development of faster hybrid methods.68−70 We have also developed a new hybrid approach based on a local exact exchange potential for that purpose, but it does not work for periodic systems at present.51 We are planning to extend it for periodic calculations in the near future. We believe that such a fast hybrid method will further promote successful applications of GDY through more reliable predictions. E

DOI: 10.1021/acsami.8b03482 ACS Appl. Mater. Interfaces XXXX, XXX, XXX−XXX

Forum Article

ACS Applied Materials & Interfaces

Spectra Using Density Functional Force Fields. J. Phys. Chem. 1994, 98, 11623−11627. (37) Curtiss, L. A.; Raghavachari, K.; Redfern, P. C.; Pople, J. A. Assessment of Gaussian-2 and Density Functional Theories for the Computation of Enthalpies of Formation. J. Chem. Phys. 1997, 106, 1063−1079. (38) Zhao, Y.; Lynch, B. J.; Truhlar, D. G. Development and Assessment of a New Hybrid Density Functional Model for Thermochemical Kinetics. J. Phys. Chem. A 2004, 108, 2715−2719. (39) Zhao, Y.; González-García, N.; Truhlar, D. G. Benchmark Database of Barrier Heights for Heavy Atom Transfer, Nucleophilic Substitution, Association, and Unimolecular Reactions and Its Use to Test Theoretical Methods. J. Phys. Chem. A 2005, 109, 2012−2018. (40) Song, J.-W.; Yamashita, K.; Hirao, K. Gaussian Attenuation Hybrid Scheme Applied to the Ernzerhof-Perdew Exchange Hole Model (Gau-PBEh). J. Chem. Phys. 2012, 137, 244105. (41) Song, J.-W.; Hirao, K. Long-Range Corrected Density Functional Theory with Accelerated Hartree-Fock Exchange Integration Using a Two-Gaussian Operator [LC-ωPBE(2Gau)]. J. Chem. Phys. 2015, 143, 144112. (42) Marsman, M.; Paier, J.; Stroppa, A.; Kresse, G. Hybrid Functionals Applied to Extended Systems. J. Phys.: Condens. Matter 2008, 20, 064201. (43) Liu, R.; Liu, H.; Li, Y.; Yi, Y.; Shang, X.; Zhang, S.; Yu, X.; Zhang, S.; Cao, H.; Zhang, G. Nitrogen-Doped Graphdiyne as a MetalFree Catalyst for High-Performance Oxygen Reduction Reactions. Nanoscale 2014, 6, 11336−11343. (44) Desroches, M.; Courtemanche, M.-A.; Rioux, G.; Morin, J.-F. Synthesis and Properties of Rhomboidal Macrocyclic Subunits of Graphdiyne-Like Nanoribbons. J. Org. Chem. 2015, 80, 10634−10642. (45) Zhang, S.; He, J.; Zheng, J.; Huang, C.; Lv, Q.; Wang, K.; Wang, N.; Lan, Z. Porous Graphdiyne Applied for Sodium Ion Storage. J. Mater. Chem. A 2017, 5, 2045−2051. (46) Chandra Shekar, S.; Swathi, R. S. Stability of Nucleobases and Base Pairs Adsorbed on Graphyne and Graphdiyne. J. Phys. Chem. C 2014, 118, 4516−4528. (47) Chandra Shekar, S.; Swathi, R. S. Molecular Switching on Graphyne and Graphdiyne: Realizing Functional Carbon Networks in Synergy with Graphene. Carbon 2018, 126, 489−499. (48) Liechtenstein, A. I.; Anisimov, V. I.; Zaanen, J. DensityFunctional Theory and Strong Interactions: Orbital Ordering in MottHubbard Insulators. Phys. Rev. B: Condens. Matter Mater. Phys. 1995, 52, R5467−R5470. (49) Wehling, T. O.; Lichtenstein, A. I.; Katsnelson, M. I. TransitionMetal Adatoms on Graphene: Influence of Local Coulomb Interactions on Chemical Bonding and Magnetic Moments. Phys. Rev. B: Condens. Matter Mater. Phys. 2011, 84, 235110. (50) He, J.; Ma, S. Y.; Zhou, P.; Zhang, C. X.; He, C.; Sun, L. Z. Magnetic Properties of Single Transition-Metal Atom Absorbed Graphdiyne and Graphyne Sheet from DFT+U Calculations. J. Phys. Chem. C 2012, 116, 26313−26321. (51) Kim, J.; Kang, S.; Lim, J.; Hwang, S.-Y.; Kim, W. Y. Practical Approach for Fast Hybrid Density Functional Calculations in RealSpace Numerical Grid Methods. Comput. Phys. Commun. 2018, DOI: 10.1016/j.cpc.2018.04.010. (52) Grimme, S.; Antony, J.; Ehrlich, S.; Krieg, H. A Consistent and Accurate Ab Initio Parametrization of Density Functional Dispersion Correction (DFT-D) for the 94 Elements H-Pu. J. Chem. Phys. 2010, 132, 154104. (53) Grimme, S.; Ehrlich, S.; Goerigk, L. Effect of the Damping Function in Dispersion Corrected Density Functional Theory. J. Comput. Chem. 2011, 32, 1456−1465. (54) Kresse, G.; Furthmüller, J. Efficient Iterative Schemes for Ab Initio Total-Energy Calculations Using a Plane-Wave Basis Set. Phys. Rev. B: Condens. Matter Mater. Phys. 1996, 54, 11169−11186. (55) Choi, S.; Hong, K.; Kim, J.; Kim, W. Y. Accuracy of LagrangeSinc Functions as a Basis Set for Electronic Structure Calculations of Atoms and Molecules. J. Chem. Phys. 2015, 142, 094116.

(17) Zhang, X.; Zhu, M.; Chen, P.; Li, Y.; Liu, H.; Li, Y.; Liu, M. Pristine Graphdiyne-Hybridized Photocatalysts Using Graphene Oxide as a Dual-Functional Coupling Reagent. Phys. Chem. Chem. Phys. 2015, 17, 1217−1225. (18) Li, Y.; Guo, C.; Li, J.; Liao, W.; Li, Z.; Zhang, J.; Chen, C. Pyrolysis-Induced Synthesis of Iron and Nitrogen-Containing Carbon Nanolayers Modified Graphdiyne Nanostructure as a Promising CoreShell Electrocatalyst for Oxygen Reduction Reaction. Carbon 2017, 119, 201−210. (19) Zhang, S.; Cai, Y.; He, H.; Zhang, Y.; Liu, R.; Cao, H.; Wang, M.; Liu, J.; Zhang, G.; Li, Y.; Liu, H.; Li, B. Heteroatom Doped Graphdiyne as Efficient Metal-Free Electrocatalyst for Oxygen Reduction Reaction in Alkaline Medium. J. Mater. Chem. A 2016, 4, 4738−4744. (20) Xue, Y.; Guo, Y.; Yi, Y.; Li, Y.; Liu, H.; Li, D.; Yang, W.; Li, Y. Self-Catalyzed Growth of Cu@graphdiyne Core−shell Nanowires Array for High Efficient Hydrogen Evolution Cathode. Nano Energy 2016, 30, 858−866. (21) Jin, Z.; Zhou, Q.; Chen, Y.; Mao, P.; Li, H.; Liu, H.; Wang, J.; Li, Y. Graphdiyne:ZnO Nanocomposites for High-Performance UV Photodetectors. Adv. Mater. 2016, 28, 3697−3702. (22) Li, J.; Gao, X.; Liu, B.; Feng, Q.; Li, X.-B.; Huang, M.-Y.; Liu, Z.; Zhang, J.; Tung, C.-H.; Wu, L.-Z. Graphdiyne: A Metal-Free Material as Hole Transfer Layer To Fabricate Quantum Dot-Sensitized Photocathodes for Hydrogen Production. J. Am. Chem. Soc. 2016, 138, 3954−3957. (23) Xiao, J.; Shi, J.; Liu, H.; Xu, Y.; Lv, S.; Luo, Y.; Li, D.; Meng, Q.; Li, Y. Efficient CH3NH3PbI3 Perovskite Solar Cells Based on Graphdiyne (GD)-Modified P3HT Hole-Transporting Material. Adv. Energy Mater. 2015, 5, 1401943. (24) Jiao, Y.; Du, A.; Hankel, M.; Zhu, Z.; Rudolph, V.; Smith, S. C. Graphdiyne: A Versatile Nanomaterial for Electronics and Hydrogen Purification. Chem. Commun. 2011, 47, 11843−11845. (25) Meng, Z.; Zhang, X.; Zhang, Y.; Gao, H.; Wang, Y.; Shi, Q.; Rao, D.; Liu, Y.; Deng, K.; Lu, R. Graphdiyne as a High-Efficiency Membrane for Separating Oxygen from Harmful Gases: A FirstPrinciples Study. ACS Appl. Mater. Interfaces 2016, 8, 28166−28170. (26) Bartolomei, M.; Carmona-Novillo, E.; Hernández, M. I.; Campos-Martínez, J.; Pirani, F.; Giorgi, G. Graphdiyne Pores: “Ad Hoc” Openings for Helium Separation Applications. J. Phys. Chem. C 2014, 118, 29966−29972. (27) Long, M.; Tang, L.; Wang, D.; Li, Y.; Shuai, Z. Electronic Structure and Carrier Mobility in Graphdiyne Sheet and Nanoribbons: Theoretical Predictions. ACS Nano 2011, 5, 2593−2600. (28) Jalili, S.; Houshmand, F.; Schofield, J. Study of Carrier Mobility of Tubular and Planar Graphdiyne. Appl. Phys. A: Mater. Sci. Process. 2015, 119, 571−579. (29) Burke, K. Perspective on Density Functional Theory. J. Chem. Phys. 2012, 136, 150901. (30) Perdew, J. P.; Burke, K.; Ernzerhof, M. Generalized Gradient Approximation Made Simple. Phys. Rev. Lett. 1996, 77, 3865−3868. (31) Seif, A.; López, M. J.; Granja-DelRío, A.; Azizi, K.; Alonso, J. A. Adsorption and Growth of Palladium Clusters on Graphdiyne. Phys. Chem. Chem. Phys. 2017, 19, 19094−19102. (32) Li, Y.; Xu, L.; Liu, H.; Li, Y. Graphdiyne and Graphyne: From Theoretical Predictions to Practical Construction. Chem. Soc. Rev. 2014, 43, 2572−2586. (33) Bu, H.; Zhao, M.; Zhang, H.; Wang, X.; Xi, Y.; Wang, Z. Isoelectronic Doping of Graphdiyne with Boron and Nitrogen: Stable Configurations and Band Gap Modification. J. Phys. Chem. A 2012, 116, 3934−3939. (34) Yang, N.; Liu, Y.; Wen, H.; Tang, Z.; Zhao, H.; Li, Y.; Wang, D. Photocatalytic Properties of Graphdiyne and Graphene Modified TiO2: From Theory to Experiment. ACS Nano 2013, 7, 1504−1512. (35) Vydrov, O. A.; Scuseria, G. E. Assessment of a Long-Range Corrected Hybrid Functional. J. Chem. Phys. 2006, 125, 234109. (36) Stephens, P. J.; Devlin, F. J.; Chabalowski, C. F.; Frisch, M. J. Ab Initio Calculation of Vibrational Absorption and Circular Dichroism F

DOI: 10.1021/acsami.8b03482 ACS Appl. Mater. Interfaces XXXX, XXX, XXX−XXX

Forum Article

ACS Applied Materials & Interfaces (56) Kang, S.; Ryu, S.; Choi, S.; Kim, J.; Hong, K.; Kim, W. Y. Update to ACE-Molecule: Projector Augmented Wave Method on Lagrange-Sinc Basis Set. Int. J. Quantum Chem. 2016, 116, 644−650. (57) Kim, J.; Hong, K.; Hwang, S.-Y.; Ryu, S.; Choi, S.; Kim, W. Y. Effects of the Locality of a Potential Derived from Hybrid Density Functionals on Kohn−Sham Orbitals and Excited States. Phys. Chem. Chem. Phys. 2017, 19, 10177−10186. (58) Lim, J.; Choi, S.; Kang, S.; Kim, J.; Hong, K.; Kim, W. Y. Improvement of Initial Guess via Grid-Cutting for Efficient Grid-Based Density Functional Calculations. Int. J. Quantum Chem. 2016, 116, 1397−1403. (59) Blöchl, P. E. Projector Augmented-Wave Method. Phys. Rev. B: Condens. Matter Mater. Phys. 1994, 50, 17953−17979. (60) Frisch, M. J.; Trucks, G. W.; Schlegel, H. B.; Scuseria, G. E.; Robb, M. A.; Cheeseman, J. R.; Scalmani, G.; Barone, V.; Mennucci, B.; Petersson, G. A.; Nakatsuji, H.; Caricato, M.; Li, X.; Hratchian, H. P.; Izmaylov, A. F.; Bloino, J.; Zheng, G.; Sonnenberg, J. L.; Hada, M.; Ehara, M.; Toyota, K.; Fukuda, R.; Hasegawa, J.; Ishida, M.; Nakajima, T.; Honda, Y.; Kitao, O.; Nakai, H.; Vreven, T.; Montgomery, J. A., Peralta, J. E.; Ogliaro, F.; Bearpark, M.; Heyd, J. J.; Brothers, E.; Kudin, K. N.; Staroverov, V. N.; Keith, T.; Kobayashi, R.; Normand, J.; Raghavachari, K.; Rendell, A.; Burant, J. C.; Iyengar, S. S.; Tomasi, J.; Cossi, M.; Rega, N.; Millam, M. J.; Klene, M.; Knox, J. E.; Cross, J. B.; Bakken, V.; Adamo, C.; Jaramillo, J.; Gomperts, R.; Stratmann, R. E.; Yazyev, O.; Austin, A. J.; Cammi, R.; Pomelli, C.; Ochterski, J. W.; Martin, R. L.; Morokuma, K.; Zakrzewski, V. G.; Voth, G. A.; Salvador, P.; Dannenberg, J. J.; Dapprich, S.; Daniels, A. D.; Farkas, Ö .; Foresman, J. B.; Ortiz, J. V.; Cioslowski, J.; Fox, D. J. Gaussian 09, Revision D.01; Gaussian, Inc.: Wallingford, CT, 2016. (61) Enkovaara, J.; Rostgaard, C.; Mortensen, J. J.; Chen, J.; Dułak, M.; Ferrighi, L.; Gavnholt, J.; Glinsvad, C.; Haikola, V.; Hansen, H. A.; Kristoffersen, H. H.; Kuisma, M.; Larsen, A. H.; Lehtovaara, L.; Ljungberg, M.; Lopez-Acevedo, O.; Moses, P. G.; Ojanen, J.; Olsen, T.; Petzold, V.; Romero, N. A.; Stausholm-Møller, J.; Strange, M.; Tritsaris, G. A.; Vanin, M.; Walter, M.; Hammer, B.; Häkkinen, H.; Madsen, G. K. H.; Nieminen, R. M.; Nørskov, J. K.; Puska, M.; Rantala, T. T.; Schiøtz, J.; Thygesen, K. S.; Jacobsen, K. W. Electronic Structure Calculations with GPAW: A Real-Space Implementation of the Projector Augmented-Wave Method. J. Phys.: Condens. Matter 2010, 22, 253202. (62) Heyd, J.; Scuseria, G. E. Efficient hybrid density functional calculations in solids: Assessment of the Heyd−Scuseria−Ernzerhof screened Coulomb hybrid functional. J. Chem. Phys. 2004, 121, 1187− 1192. (63) Onida, G.; Reining, L.; Rubio, A. Electronic excitations: densityfunctional versus many-body Green’s-function approaches. Rev. Mod. Phys. 2002, 74, 601−659. (64) Becke, A. D. Density-Functional Exchange-Energy Approximation with Correct Asymptotic Behavior. Phys. Rev. A: At., Mol., Opt. Phys. 1988, 38, 3098−3100. (65) Lee, C.; Yang, W.; Parr, R. G. Development of the Colle-Salvetti Correlation-Energy Formula into a Functional of the Electron Density. Phys. Rev. B: Condens. Matter Mater. Phys. 1988, 37, 785−789. (66) Kittel, C. In Introduction to Solid State Physics, 8th ed; John Wiley & Sons: Hoboken, NJ, 2005; pp 50−50. (67) Faglioni, F.; Merinov, B. V.; Goddard, W. A. Structural and Electronic Properties of Graphdiyne Carbon Nanotubes from LargeScale DFT Calculations. J. Phys. Chem. C 2016, 120, 27104−27108. (68) Hu, W.; Lin, L.; Banerjee, A. S.; Vecharynski, E.; Yang, C. Adaptively Compressed Exchange Operator for Large-Scale Hybrid Density Functional Calculations with Applications to the Adsorption of Water on Silicene. J. Chem. Theory Comput. 2017, 13, 1188−1198. (69) Marsili, M.; Umari, P. Method for the Fast Evaluation of Fock Exchange for Nonlocalized Wave Functions. Phys. Rev. B: Condens. Matter Mater. Phys. 2013, 87, 205110. (70) Duchemin, I.; Gygi, F. A Scalable and Accurate Algorithm for the Computation of Hartree−Fock Exchange. Comput. Phys. Commun. 2010, 181, 855−860.

G

DOI: 10.1021/acsami.8b03482 ACS Appl. Mater. Interfaces XXXX, XXX, XXX−XXX