We report electronic structure calculations in chemically functionalized linear cubane-based chains. The effects of covalent chemical attachments on chain transport properties are examined with nonorthogonal tight-binding model (NTBM) considering Landauer-Büttiker formalism. The covalent bonding of even a single-type functional group is shown to considerably alter the conductance of the chain. For similar radical doping density, electronic characteristics are found to range from insulator to narrow-gap semiconductor depending on the nature of the covalent bonding. Therefore it has become possible to tune electronic properties of the cubane-based one-dimensional oligomers by the functionalization for nanoelectronic applications.
1. Introduction
Hydrocarbon cubane C8H8 synthesized in 1964 [1] remains the focus of modern research because of its unusual geometry, energy, and electronic properties. In cubane molecule the valence angles between the C–C covalent bonds are 90°, distinct from the value 109.5°, typical of sp3-hybridized carbon compounds. Nevertheless, despite the strained carbon cage (strain energy is about 670 kJ/mol [2]) cubane molecules demonstrate surprisingly high kinetic stability [3] and can form various derivatives, in which one or more hydrogen atoms are replaced by functional groups or radicals. For example, cubane-based systems with methyl groups [4], propyl groups [4], carboxyl groups [5], and nitrogroups [6–8] are successfully synthesized. Moreover, dehydrogenated cubanes can form oligomers through covalent bonding. The cubane dimer (cubylcubane) [9], linear chains constructed from elementary cubyl fragments [5], and cubane-based solids [10, 11] are intensively studied as well. Although a number of cubane-based oligomers remain hypothetical and are not obtained experimentally, they are examined theoretically, including numerical simulation. Both ab initio [12, 13] and semiempirical [14, 15] calculations predict thermodynamic and kinetic stabilities of one- and two-dimensional nanostructures built from cubylene units (chains, rings, and networks). These systems are of particular interest because of their remarkable electronic properties [16–19]. For example, polycubane chains linked by various small molecules (“molecular bridges”) demonstrate a variety of electronic characteristics: from strong insulators to metals [17]. Thus, it becomes possible to design novel advanced materials with adjustable electronic structure. Recently reported charge transfer between donor and acceptor groups in cubane-based one-dimensional systems makes these nanostructures perspective materials to be used in pull-push devices [18].
Linear chains built from cubylene units [5] are the simplest members of the cubane-based structures family. Ab initio calculations [13] predict that the sufficiently long chain is a strong insulator with the energy gap equal to 4.59 eV, but this value can be reduced through doping the initial system by donor/acceptor functional groups [13, 17]. Nevertheless, the mechanism of the effect of chemical functionalization on the energy gap of cubane-based nanostructures is still not clear.
In this paper we present electronic structure calculations for the linear chains built from substituted cubylene units. In each unit one or more hydrogen atoms are replaced by functional groups or radicals. We restrict ourselves to H–C–N–O chemical groups with various electronegativities, for which corresponding substituted cubanes are already synthesized or can be easily prepared in the future. The main purpose of this work is to demonstrate the electronic properties of substituted cubane-based chains dependence of the chemical groups attached and to calculate their conductance.
2. Materials and Methods
In our study we use the nonorthogonal tight-binding model originally developed for H–C–N–O systems [20, 21] as well as high-level density functional theory (DFT). Density functional theory with B3LYP [22, 23] and PBE [24] functionals with the electron basis sets of 6-311G/6-311G(d,p) and 6-31+G(d,p) [25] was used. All DFT calculations were performed using the GAMESS program package [26]. The parameters’ fitting of tight-binding model is based on the criterion of the best correspondence between the computed and experimental (not derived from ab initio calculations) values of binding energies, bond lengths, and valence angles of several selected small HkClNmOn molecules (for parameterization details see [20, 27]). To validate the model developed, we computed the binding energies, bond lengths, and valence angles for various H–C–N–O compounds not used in the fitting procedure and compared them with the experimental data [28]. The resulting tight-binding potential has a good transferability and is well suited to modeling various HkClNmOn molecules [27] from small clusters to large systems as graphene, diamond, different peptides, and so forth. Such tight-binding approach is commonly used for electronic structure calculations [29–32], including HOMO-LUMO gaps, density of states, and transmission coefficients.
The samples studied are constructed from eleven functionalized cubylene fragments. Each structure is obtained from pristine chain (see Figure 1(a)) by replacing a hydrogen atom with chemical functional group or radical: CH3, OH, NH2, CN, COOH (Figure 1(b)), NO2, and NO (Figure 1(c)). The choice of these functional groups will be clear from the analysis below. So, the total number of oligomers investigated is equal to eight.
Fragments of the chemically functionalized cubane-based chains studied: (a) pristine chain; (b) COOH-functionalized chain; (c) NO-functionalized chain.
To obtain the equilibrium structures of the cubane-based oligomers we applied the method of structural relaxation by using the same nonorthogonal tight-binding model (NTBM) [20, 21, 27]. The corresponding initial configuration relaxed to a state with the local or global energy minimum under the influence of intramolecular forces only. First of all, the forces acting on all atoms were calculated using the Hellmann-Feynman theorem. Next, atoms were shifted in the direction of the forces obtained proportional to the corresponding forces. Then the relaxation step was repeated. The whole structure is relaxed up to residual atomic forces smaller than 10-4 eV/Å. Note that for all cubane-based chains studied we also calculated the frequency spectra in the framework of the same level of theory after the geometry optimization. For each metastable configuration, the presence of local minima of energy was confirmed by the reality of all frequencies. We have verified on the smaller oligomers that optimizing the structures in this way leads to minor changes in the energy and electronic characteristics (e.g., binding energies and HOMO-LUMO gaps) as compared to high-level DFT relaxation. Results obtained for HOMO-LUMO gaps for cubane-based oligomers constructed from two and four fragments (dimers and tetramers) are presented in Table 1. From Table 1 one can see that the accuracy of the NTBM is in good agreement with the DFT approach (except, maybe, CN-functionalization).
HOMO-LUMO gaps ΔHL for various cubane-based oligomers constructed from two and four fragments.
Radical
ΔHL (eV)
Dimers
Tetramers
NTBM
DFT/PBE/6-31+G(d,p)
NTBM
DFT/PBE/6-31+G(d,p)
H (pristine chain)
5.17
5.31
4.91
5.02
CH3
5.05
5.06
4.80
4.81
OH
4.59
4.34
4.38
4.13
NH2
4.49
3.92
4.27
3.73
CN
3.40
5.21
3.33
4.97
COOH
3.01
3.64
2.98
3.62
NO2
2.27
3.20
2.25
3.06
NO
0.60
1.17
0.60
1.08
To calculate the electronic density of states and the conductance of infinite pristine and functionalized cubane-based chains we have applied the Landauer-Büttiker formalism [33–36]. The general illustration of our scheme is presented in Figure 2.
General scheme of the considered samples comprising three parts: semi-infinite left and right leads and device area, consisting of three substituted cubylene fragments, sandwiched between them.
The device region (D) consisting of three substituted cubylene units is connected to left (L) and right (R) semi-infinite leads. The device area is described by the Hamiltonian HD and corresponding overlap matrix SD. The conductance at zero temperature and in low bias limit can be obtained from Landauer’s theory for transport(1)GE=G0TrΓLGΓRG†,where G is the conductance expressed within Fisher-Lee relationship [34], while G0=2e2/h is the conductance quantum, E is the energy of the incident charge carrier, and Tr denotes the trace of corresponding operator; G is the retarded Green’s function of the functionalized linear cubane-based chain and matrices ΓL,R=iΣL,R-ΣR,L† consider the coupling to the right or left leads through their self-energies ΣL,R. The self-energies of the device-leads system can be described as(2)ΣL=E+iηSLD-HLD†·E+iηSL0-HL0+E+iηSL01-HL01†T-L-1·E+iηSLD-HLD;ΣR=E+iηSDR-HDR·E+iηSR0-HR0+E+iηSR01-HR01†TR-1·E+iηSDR-HDR†,where η is an arbitrary small quantity; HLD, HDR, SLD, and SDR are the coupling Hamiltonians and corresponding overlap matrices describing connections device left lead and device right lead; HL0, HR0, HL01, HR01, SL0, SR0, SL01, and SR01 describe the left/right lead unit cell and coupling between the neighboring unit cells of the lead; T-L and TR are the lead transfer matrices that can be easily computed from HL0, HL01, SL0, SL01 and HR0, HR01, SR0, SR01 via an iterative procedure (see [37–40] for details). Note that in our case left and right leads as well as device region have the same molecular structure. Therefore HL0≡HR0≡HD, SL0≡SR0≡SD, and HL01≡HR01≡HLD≡HDR, SL01≡SR01≡SLD≡SDR. The retarded Green’s function is given by the expression(3)GE=E+iηSD-HD-ΣL-ΣR-1.The electronic density of states (DOS) is given by(4)DOSE=-12πTrG-G†SD.
The electronic Hamiltonians are calculated from overlap matrix elements using extended Hückel approximation [41] with Anderson distant dependence for Wolfsberg-Helmholtz parameter [42]. The overlap integrals are calculated using standard Slater-Koster-Roothaan procedure [43, 44]. Therefore, the analytical expressions for the overlap integrals greatly simplify the calculation of Hamiltonian. The abovementioned fact as well as good accuracy of NTBM model make this approach very efficient for recursive Green’s function method as described in [39].
3. Results and Discussion
At first we focus on the case of finite linear cubane-based chains (oligomers) constructed from eleven substituted cubylene fragments. In each fragment only one hydrogen atom is replaced by the following functional groups or radicals: OH, CN, NO, NH2, NO2, CH3, COOH, N, O, CO, and NH. Geometry optimization data demonstrate the instability of cubane-based chains with N, O, CO, and NH groups. Only univalent radicals conserve the cubic identity of the elementary units. Thus, polyvalent radicals mentioned above will not be further considered here. We calculated binding energies Eb and HOMO-LUMO gaps ΔHL of the remaining chains. The binding energy Eb of the chain per atom is determined by the equation(5)EbeV/atom=1NatkEH+lEC+mEN+nEO-Etotchain,where Nat=k+l+m+n is the total number of atoms in the chain, Etotchain is the total chain (oligomer) energy, E(H), E(C), E(N), and E(O) are the energies of the isolated hydrogen, carbon, nitrogen, and oxygen atoms, respectively. As evident from Table 2ΔHL is very sensitive to functional group attached. The narrowest gap is observed for NO-doping (Figure 1(c)); it is equal to 0.60 eV. Note that HOMO-LUMO gap is not directly dependent on the binding energy (see Table 2). Therefore, it seems possible to adjust the electronic structure in a wide range by radical substitution of hydrogen from cubylene unit without loss of thermodynamic stability of the samples.
Binding energies Eb and HOMO-LUMO gaps ΔHL for various cubane-based oligomers obtained within NTBM model.
Radical
ΔHL (eV)
Eb (eV/atom)
H (pristine chain)
4.82
4.69
CH3
4.72
4.59
OH
4.30
4.65
NH2
4.19
4.54
CN
3.30
4.96
COOH
3.04
4.77
NO2
2.25
4.68
NO
0.60
4.69
To study the doping effect in more detail we replace the second hydrogen atom in each elementary unit of NO-doped chain by another NO-group (see Figure 3). As a result ΔHL slightly decreases from 0.60 to 0.51 eV. Note that the subsequent increase of the number of NO-groups in the chain leads to the stability loss of cubic framework. Probably it may be concerned with the adjacency of functional groups attached. Earlier we found that their “effective repulsion” makes the synthesis of isolated cubanes with adjacent nitrogen-containing groups quite difficult [20]. Thus, we were unable to make the HOMO-LUMO gap disappear and convert the cubane-based chain to metallic state. However, we assume the existence of such combination of radicals (not only H–C–N–O) which can close the HOMO-LUMO gap in electronic structure of cubane-based oligomers.
Fragment of the double NO-functionalized cubane-based chain.
We further explore the conduction regime in chemically functionalized infinite cubane-based chains by calculating the conductance for chains with varying radicals attached. The dependence of the functionalized chain conductance as well as DOS as a function of the radical attached is illustrated in Figure 4 for a set of typical configurations. It is shown that the functionalization alters the electronic properties of the chain since the dopants states are located near the Fermi level. So, the pristine cubane-based chain is an insulator, while NO-group substituted chain is a narrow-gap semiconductor.
Density of electronic states (blue) and conductance (red, in units of G0=2e2/h) of the infinite cubane-based chains for the pristine case (a) and for the COOH- and NO-functionalized cases ((b) and (c), resp.).
Note that energy gaps for infinite substituted cubane-based chains are in good agreement with the HOMO-LUMO gaps of corresponding eleven-fragment samples. Our numerical calculations show that in pristine cubane-based chain HOMO-LUMO excitation was characterized by electronic density transfer to the ends of chain. By contrast, in narrow-gap substituted cubane-based chains (e.g., NO-functionalized chain) electronic density rearranges to the functional groups. So, the fine adjustment of HOMO-LUMO gap has its origins in electronic density transfer to the corresponding radicals. Therefore the HOMO-LUMO gap depends on the functional groups selected. The important consequence is that the functionalization of the cubane-based chain leads to the opening of new conduction channels (see Figure 4): for functionalized chains the energy gap becomes narrower since the dopants states are located near the Fermi level. Therefore, they exhibit semiconducting properties, and conductivity becomes possible at nonzero temperatures. This fact indicates that the functionalized cubane-based chains can be potentially used as effective conductors in nanodevices.
4. Conclusions
We have studied electronic properties of the cubane-based compounds containing different chemical groups with various electronegativities. It is shown that such type of functionalization alters the electronic properties of the chain since the dopants states are located near the Fermi level. Some functional groups or radicals can significantly reduce the energy gap by opening the additional conduction channels in the molecular structure. Moreover, the increase of the number of chemical groups (e.g., NO-groups) leads to further decreasing the gap, but on the other hand it may lead to the stability loss and decomposition of the chain. So, we were unable to convert the cubane-based chain to metallic state by removing the hydrogen and adding the functional group instead of it. Maybe it is possible to make the gap disappear by implementing the additional “molecular bridges” between the substituted cubylene fragments of the cubane-based compound. Anyway it has become possible to tune electronic properties of the cubane-based chains by the chemical functionalization that will be useful for nanoelectronic applications.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
This work was supported by the Russian Foundation for Basic Research, Project no. 14-02-31416-mol_a. The useful comments of anonymous reviewers are gratefully acknowledged.
EatonP. E.ColeT. W.Jr.The cubane system196486596296410.1021/ja01059a0722-s2.0-0000483172AgapitoF.SantosR. C.Borges dos SantosR. M.Martinho SimõesJ. A.The thermochemistry of cubane 50 years after its synthesis: a high-level theoretical study of cubane and its derivatives2015119122998300710.1021/jp511756vMaslovM. M.LobanovD. A.PodlivaevA. I.OpenovL. A.Thermal stability of cubane C8H8200951364564810.1134/s10637834090303422-s2.0-62249124665EatonP. E.LiJ.UpadhyayaS. P.Synthesis of methylcubane and cyclopropylcubane. The cubane-1,4-diyl route199560496696810.1021/jo00109a0302-s2.0-0012336161EatonP. E.PramodK.EmrickT.GilardiR.Building with cubane-1,4-diyl. Synthesis of aryl-substituted cubanes, p-[n]cubyls, and cubane-separated bis(arenes)1999121174111412310.1021/ja983441fEatonP. E.GilardiR. L.ZhangM.-X.Polynitrocubanes: advanced high-density, high-energy materials20001215114311482-s2.0-0034246156ZhangM. X.EatonP. E.GilardiR. L.Hepta- and octanitrocubanes2000392401404EatonP. E.ZhangM.-X.GilardiR. L.GelberN.IyerS.SurapaneniR.Octanitrocubane: a new nitrocarbon20022711610.1002/1521-4087(200203)27:1lt;1::aid-prep1>3.0.co;2-62-s2.0-0036502463GilardiR.MagginiM.EatonP. E.X-ray structures of cubylcubane and 2-tert-butylcubylcubane: short cage-cage bonds1988110217232723410.1021/ja00229a058YildirimT.Structure and dynamics from combined neutron scattering and first-principles studies20002611-220521610.1016/s0301-0104(00)00231-7YildirimT.CiraciS.KiliçÇ.BuldumA.First-principles investigation of structural and electronic properties of solid cubane and its doped derivatives200062117625763310.1103/physrevb.62.76252-s2.0-0034664480ValenciaF.RomeroA. H.KiwiM.RamírezR.Toro-LabbeA.Ab initio study of cubyl chains and networks2004121189172917710.1063/1.18054922-s2.0-9744236583HerreraB.ValenciaF.RomeroA. H.KiwiM.RamírezR.Toro-LabbéA.Cubane oligomers: a density functional theory study20067691–318318710.1016/j.theochem.2006.04.0532-s2.0-33748599255MaslovM. M.Theoretical modeling of the cubane-based chains, networks and bulks200812901204510.1088/1742-6596/129/1/012045MaslovM. M.KatinK. P.AvkhadievaA. I.PodlivaevA. I.Dynamic modeling of cyclotetracubyl thermal decomposition20148215215710.1134/S19907931140201952-s2.0-84901007040KonstantinovaE.CamiloA.Jr.BaroneP. M. V. B.DantasS. O.GalvãoD. S.Some electronic properties of saturated and unsaturated cubane oligomers using DFT-based calculations20088681–3374110.1016/j.theochem.2008.07.0342-s2.0-53249121013ValenciaF.RomeroA. H.KiwiM.RamírezR.Toro-LabbeA.Polycubanes linked with C2, N2, NO, and NS: from insulating to metallic behavior200571303341010.1103/physrevb.71.0334102-s2.0-15444379021KonstantinovaE.BaroneP. M. V. B.DantasS. O.On the dependence of electronic properties on geometric structure for cubane derivatives and oligomers with donor and acceptor groups20099091–311612110.1016/j.theochem.2009.06.0012-s2.0-67650707939PalP. P.PatiR.Charge transport in strongly coupled molecular junctions: ‘in-phase’ and ‘out-of-Phase’ contribution to electron tunneling201111535175641757310.1021/jp204711c2-s2.0-80052311093KatinK. P.MaslovM. M.Thermal stability of nitro derivatives of hydrocarbon cubane20115577077910.1134/S1990793111090181Nonorthogonal Tight-Binding Model version number 1, http://ntbm.info/LeeC.YangW.ParrR. G.Development of the Colle-Salvetti correlation-energy formula into a functional of the electron density198837278578910.1103/physrevb.37.785BeckeA. D.Density-functional thermochemistry. III. The role of exact exchange1993987, article 564810.1063/1.464913PerdewJ. P.BurkeK.ErnzerhofM.Generalized gradient approximation made simple199677183865386810.1103/PhysRevLett.77.38652-s2.0-4243943295KrishnanR.BinkleyJ. S.SeegerR.PopleJ. A.Self-consistent molecular orbital methods. XX. A basis set for correlated wave functions198072165065410.1063/1.4389552-s2.0-26844534384SchmidtM. W.BaldridgeK. K.BoatzJ. A.ElbertS. T.GordonM. S.JensenJ. H.KosekiS.MatsunagaN.NguyenK. A.SuS.WindusT. L.DupuisM.MontgomeryJ. A.General atomic and molecular electronic structure system199314111347136310.1002/jcc.5401411122-s2.0-84893169025MaslovM. M.PodlivaevA. I.KatinK. P.Nonorthogonal tight-binding model with H–C–N–O parameterisation201510.1080/08927022.2015.1044453NIST Standard Reference Database 101, Computational Chemistry Comparison and Benchmark DataBase, http://cccbdb.nist.gov/JavanM. B.Electronic transport properties of linear nC20 (n≤5) oligomers: theoretical investigation20156713514210.1016/j.physe.2014.11.0082-s2.0-84918811482GhobadiN.AbdiY.Device characteristics and tight binding based modeling of bilayer graphene field-effect transistor20131361082108910.1016/j.cap.2013.02.0162-s2.0-84877581380RabaniH.MardaaniM.ShahrakiA. M.Analytical study of the ballistic transport of ladder-like graphene nanoribbons within the tight-binding approach20135910611410.1016/j.spmi.2013.03.0252-s2.0-84877355175SukkabotW.Structural properties of SiC zinc-blende and wurtzite nanostructures: atomistic tight-binding theory20154011712210.1016/j.mssp.2015.06.079LandauerR.Electrical resistance of disordered one-dimensional lattices19702117286386710.1080/14786437008238472FisherD. S.LeeP. A.Relation between conductivity and transmission matrix19812312, article 685110.1103/physrevb.23.68512-s2.0-4243971222BüttikerM.ImryY.LandauerR.PinhasS.Generalized many-channel conductance formula with application to small rings198531106207621510.1103/physrevb.31.62072-s2.0-8744270531BüttikerM.Four-terminal phase-coherent conductance198657141761176410.1103/physrevlett.57.1761SanchoM. P. L.SanchoJ. M. L.RubioJ.Quick iterative scheme for the calculation of transfer matrices: application to Mo (100)19841451205121510.1088/0305-4608/14/5/0162-s2.0-0000821265SanchoM. P. L.SanchoJ. M. L.RubioJ.Highly convergent schemes for the calculation of bulk and surface Green functions1985154, article 85110.1088/0305-4608/15/4/0092-s2.0-0001124898NardelliM. B.Electronic transport in extended systems: application to carbon nanotubes199960117828783310.1103/physrevb.60.7828MarkussenT.RuraliR.BrandbygeM.JauhoA.-P.Electronic transport through Si nanowires: role of bulk and surface disorder2006742424531310.1103/physrevb.74.2453132-s2.0-33845519735HoffmannR.An extended Hückel theory. I. Hydrocarbons1963396139710.1063/1.1734456AndersonA. B.Derivation of the extended Hückel method with corrections: one electron molecular orbital theory for energy level and structure determinations19756231187118810.1063/1.430562SlaterJ. C.KosterG. F.Simplified LCAO method for the periodic potential problem19549461498152410.1103/physrev.94.14982-s2.0-2442522754RoothaanC. C. J.A study of two-center integrals useful in calculations on molecular structure. I1951191445145810.1063/1.1748100