Carbon Nano-onions as Photosensitizers: Stacking-Induced Red-Shift

Jan 4, 2018 - We show that the likely reason for this immobility of the outer layers is the superadditivity of the interaction energies between indivi...
0 downloads 11 Views 13MB Size
Subscriber access provided by UNIV OF TASMANIA

Article

Carbon Nano-Onions as Photo-Sensitizers: Stacking Induced Red-Shift Muhammad Ali Hashmi, and Matthias Lein J. Phys. Chem. C, Just Accepted Manuscript • DOI: 10.1021/acs.jpcc.7b11421 • Publication Date (Web): 04 Jan 2018 Downloaded from http://pubs.acs.org on January 4, 2018

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

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

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

The Journal of Physical Chemistry

Carbon Nano-Onions as Photo-Sensitizers: Stacking Induced Red-Shift Muhammad Ali Hashmi† and Matthias Lein∗,†,‡ †School of Chemical and Physical Sciences (SCPS), Victoria University of Wellington, New Zealand ‡Centre for Theoretical Chemistry and Physics (CTCP), New Zealand Institute for Advanced Study, Massey University Auckland E-mail: [email protected]

1

ACS Paragon Plus Environment

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

Abstract We have investigated the double-layered carbon nano-onions (CNOs) C20 @C60 , C20 @C80 , C60 @C180 , C60 @C240 , C80 @C240 , and C240 @C540 as well as their triplelayered analogues C60 @C240 @C540 , and C80 @C240 @C540 with high-level electronic structure calculations. We were able to show that earlier work, which showed the free rotation of the inner Fullerene in a CNO, doesn’t extend to multi-layered CNOs. We show that the likely reason for this immobility of the outer layers is the superadditivity of the interaction energies between individual shells; i.e. the total interaction energy of all shells is larger than the sum of all individual interaction energies between pairs of shells. We also show how the electronic states of individual Fullerenes are polarised, but essentially preserved in CNOs, and how charge-transfer excitations between layers arise, that are significantly red-shifted to lower absorption energies compared to the free Fullerenes. This allows CNOs to be utilised as photo-sensitizers in nanotechological applications.

1

Introduction

Carbon nano-onions (CNOs) are nanoparticles of pure carbon made from stacked fullerene shells and are variously known in the literature as carbon nano-onions, onion-like carbon, onion-like Fullerenes multishell Fullerenes, hyperfullerenes and nested Fullerenes. These nanoparticles may be spherical or polyhedral in shape and are often smaller than 10 nm in diameter. 1 CNOs were first discovered by Iijima, 2 several years before Kroto and Smalley’s seminal synthesis of Buckminster fullerene 3 (see Figure 1), and later described in detail by Ugarte. 4 Since their discovery, CNOs have been prepared in a number of experiments 5–10 and their structure, 11 morphology, 12 and other properties have been investigated systematically. 13–17 CNOs have been suggested for use in applications like energy storage, 18,19 high-performance wear resistance materials, 20,21 super-conductive materials, 22–24 and biomaterials. 20,25–27 Their synthesis, chemistry, and applications have been reviewed in detail. 28 2

ACS Paragon Plus Environment

Page 2 of 30

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

The Journal of Physical Chemistry

Figure 1: Structures of fullerenes optimized at the RI-PBE0-D3/def2-SVP level of theory. (a) C20 , (b) C60 , (c) C80 , (d) C180 , (e) C240 , (f) C540 . Even though an economical method that allows the synthesis of CNOs from nanodiamond particles with high yields exists, 29 it is still experimentally diffcult to synthesise large quantities of double and triple layered CNOs selectively. Hence, these are usually explored in more detail using theoretical studies. High level theoretical studies on CNOs are also limited due to their larger number of constituent atoms. Literature witnesses little density functional theory (DFT) studies on CNOs. The first report on the total energy of C60 @C240 @C540 using DFT calculations was published in 1998 by Terrones et al. 30 Later on T¨ urker studied the bucky onion of C20 @C60 at semi empirical calculations and showed that it is a highly endothermic system. 31 He further studied the formation and stability of C60 @C180 system on AM1 level and showed that the structure is stable in terms of total energy but its formation is endothermic. 32 Zheng et al. studied C20 @C60 at DFT-B3LYP/3-21G* level of theory and reported the same results that the system is stable but highly exothermic. 33 The reason for that is the shorter inter-shell distance than the observed average inter-layer distance in C60 @C240 which is 3.4 ˚ A. 34 Zope and Dunlap developed an analytical variant of DFT called grid-free ADFT which can study larger fullerenes at all-electron level. 35 Zope studied the electronic structure and static dipole polarizability of C60 @C240 using ADFT/6-311G** level 3

ACS Paragon Plus Environment

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

of theory and determined that the outer C240 almost completely shields the inner C60 from the applied static electric field. 36 Another study on nested fullerenes described the movement of C20 inside C240 and C60 inside C540 in potential wells in C20 @C240 and C60 @C540 , respectively, and concluded that C20 is neither necessarily positioned in the center of outer C240 shell nor attached to its walls. It has been shown that C20 moves with a regularity inside the outer shell of C240 even at low temperatures. 37 Recently, Salcedo et al. studied two complexes of nested fullerenes, i.e. C20 @C116 and C60 @C180 using m05-2X/6-31G** level of theory. It was concluded that the former system has a lot of interactions between the two cages but the latter was thermodynamically not feasible and high energy. It has also been observed/proposed that electrons can move between the two cages of the nested fullerene in the form of free radicals. Very recently, Zope et al. studied the site specific atomic polarizabilities of CNOs including C60 @C180 and C60 @C240 and concluded that the encapsulation of C60 inside these fullerenes reduces its polarizability by 83% and 75%, respectively. It has been proposed that the outer atoms in the fullerene shell have major contribution towards the fullerene polarizability. 38 In the present work, we studied the structural properties of C20 @C60 , C20 @C80 , C60 @C180 , C60 @C240 , C80 @C240 , C240 @C540 , C60 @C240 @C540 , and C80 @C240 @C540 CNOs using state of the art DFT calculations and the energetics of the gyroscopic motion of the inner fullerene cage. The stability of these systems alongside their interaction energies and formation energies per carbon atom have also been studied and discussed. Single point energy calculations were also performed on the optimized geometries to shed light on frontier orbitals of these systems.

2

Computational Details

Computations have been carried out using Adamo and Barrone’s hybrid formulation of Perdew, Burke, and Ernzerhof’s GGA density functional 39–41 (PBE0) in conjunction with

4

ACS Paragon Plus Environment

Page 4 of 30

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

The Journal of Physical Chemistry

Grimmes empirical dispersion correction (D3) 42 for structural optimisations. Second order Møller-Plesset perturbation theory (MP2) calculations were carried out for selected singlepoint energy calculations. 43–46 The resolution of the identity (RI) approximation 47–49 within the chain of spheres (COS) approximation 50–52 was employed with corresponding auxiliary basis set 53 for computational efficiency. Alrichs double-ζ basis set 54–59 (def2-SVP) was used for structural optimisation throughout. Frequency calculations were carried out for the smaller systems (i.e. those with fewer than 200 atoms, cf. discussion) at the same (density functional) level of theory to confirm the nature of the stationary points. All examined structures were identified as minima on the potential energy hypersurface by the absence of imaginary frequencies. We abbreviate these levels of theory as RI-PBE0-D3/def2-SVP and RI-MP2/def2-SVP respectively. Single point energy calculations were carried out with a larger basis set of triple-ζ quality 54–59 (def2-TZVP) using the previously obtained structures at the RI-PBE0-D3/def2TZVP or RI-MP2/def2-TZVP level of theory with a very fine grid for the numerical integration and very tight SCF convergence criteria in order to achieve greater numerical accuracy. Excited state spectra were calculated with the ZINDO/S method 60–63 at the previously obtained geometries. These calculations were performed using the quantum chemical methods described above as implemented in ORCA 3.0.3. 64 Frontier molecular orbital (FMO) visualisations were produced with GaussView 5.0.9 65 and are based on single point PBE0-D3/def2-SVP calculations using Gaussian 09 (Revision D.01) 66 with the previously obtained geometries. Molecular structures were visualised using CylView. 67

5

ACS Paragon Plus Environment

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

3

Results and discussion

Of all possible carbon nano-onions, only a few have been characterised and studied experimentally. Among those that we were able to identify in the literature are the doubly-layered C60 @C240 and C240 @C540 , as well as the triply-layered C80 @C240 @C540 and C60 @C240 @C540 . 68,69 Hence, we base our theoretical investigation on these experimentally confirmed systems as well as a few others that are structurally similar. CNOs with many more than three layers have been observed with high-resolution transmission electron microscopy (HRTEM) and we aim to shed light on the interaction of these larger systems through accurate computation of smaller analogues and subsequent extrapolation. The main difficulty from a computational perspective is the relatively large number of electrons in the systems. Currently, the only quantum chemical method that allows the routine calculation of systems of this size is density functional theory in the Kohn-Sham formulation (KS-DFT). This particular approach of density functional approximations (DFA) has long been known to suffer from a poor description of dispersion interactions of the type that can be expected to be of significance in the kind of stacked or nested π-systems that we are investigating here. However, us and others were able to demonstrate that the combination of a hybrid-DFA (like PBE0) with an empirical correction for dispersion interactions (like D3) can successfully describe the intra and intermolecular forces even with a moderately sized basis-set (like def2-SVP). 70–73 In order to explore the systematic error that is associated with the RI-PBE0-D3/def2SVP and RI-PBE0-D3/def2-TZVP methods that we use in our discussion, we also calculated the interaction energies of some smaller CNOs with a highly accurate ab initio method, RI-MP2/def2-TZVP. This method includes a much more accurate treatment of the intermolecular dispersion forces between the different shells of the CNO as an implicit part of its treatment of dynamic electron correlation through many body perturbation theory. It has been already described in the literature that CNOs stabilize the inner fullerene cage which plays the role of a guest in these cages. 74 The electronic behaviour, structural changes during the formation of onions, the gyroscopic movement of the guest cages, and 6

ACS Paragon Plus Environment

Page 6 of 30

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

The Journal of Physical Chemistry

the volume of these CNOs have been discussed in detail below.

3.1

Structural Properties of Carbon Nano Onions

Initial structures for the C20 , C60 , C80 , C180 , C240 and C540 fullerenes were obtained from the online database of the Michigan State University Computational Nanotechnology Lab. 75 While all, apart from C80 (which has D5h symmetry), fullerene starting geometries possess high (icosahedral) symmetry, they were optimised without any symmetry constraints in either the SCF or the geometry optimisation cycles. Nevertheless, all structures that were obtained after optimisation had retained the same point group symmetry (see Figure 1).

Figure 2: Structures of double-layered CNOs optimized at the RI-PBE0-D3/def2-SVP level of theory. (a). C20 @C60 , (b). C20 @C80 , (c). C60 @C180 , (d). C60 @C240 , (e). C80 @C240 , (f). C240 @C540 . In the next step, the doubly layered structures for C20 @C60 , C20 @C80 , C60 @C180 , C60 @C240 , C80 @C240 and C240 @C540 were produced by stacking the optimised structures of the individual fullerenes and subsequent structural optimisation. In this procedure, the two constituent fullerenes were aligned along their respective principal axes of inertia in all three possible orientations. In this work, we report only the lowest energy structure that was obtained for each CNO (see Figure 2). 7

ACS Paragon Plus Environment

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

Page 8 of 30

Figure 3: Structures of triple-layered CNOs optimized at the RI-PBE0-D3/def2-SVP level of theory. (a). C60 @C240 @C540 (b). C80 @C240 @C540 Finally, the triply layered structures C60 @C240 @C540 and C80 @C240 @C540 were produced and then optimised in the same fashion (see Figure 3). To compare all fullerenes and CNOs on the same basis, we calculated their energies of formation (Eform ) per atom (see Equation 1), with respect to the energy of a single free carbon atom (EC ).

Eform =

Etotal − nEC n

(1)

The resulting values of Eform (see Table 1) confirm that all systems considered in this study are indeed thermodynamically stable. Table 1 also shows that fullerenes with a number of carbon atoms in the range of n = 20 to 540, the formation is more exothermic for larger fullerenes than for smaller ones. The relationship between size and inter-fullerene interaction in the CNOs is more complicated, but not counter-intuitive. The sterically crowded systems C20 @C60 and C20 @C80 show values for Eform that is much closer to the value of C20 than it is to the much more stable C60 or C80 . This is an indication that the interaction between the two fullerenes in the two CNOs is repulsive, which is a likely consequence of the close approach of the two nested fullerenes. The C60 @C180 CNO shows a value of Eform that is a little higher than the arithmetic mean of the (unweighted) Eform values of the individual fullerenes. Once, the number of atoms in the two constituent 8

ACS Paragon Plus Environment

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

The Journal of Physical Chemistry

Table 1: Fullerene and CNO formation energies per atom at the RI-PBE0D3/def2-TZVP level. All energies are in kcal mol−1 . Fullerene/CNO C20 C60 C80 C180 C240 C540 C20 @C60 C20 @C80 C60 @C180 C60 @C240 C80 @C240 C240 @C540 C60 @C240 @C540 C80 @C240 @C540

Eform -126.075 -145.640 -146.708 -150.704 -151.634 -153.164 -128.817 -139.771 -148.886 -151.007 -150.878 -153.284 -152.959 -152.916

fullerenes (i.e. 60 and 180) are considered, it becomes clear that the loss incurred for the 180 carbon atoms of C180 is not compensated by the gain of the 60 carbon atoms of C60 and we can conclude that C60 @C180 is also a system that is dominated by repulsion between the two layers. The remaining doubly and triply layered CNOs all show Eform values that large enough to indicate attractive inter-layer interactions. This data is confirmed once interaction energies Eint are calculated (see Equation 2) directly from the total energies of the CNOs (ECNO ) and the total energies of the inner and outer layers (Eouter and Einner ).

Eint = ECNO − (Eouter + Einner )

(2)

The Eint values (see Table 2) show that there are three noticeable trends within the interaction energies. Firstly and most obviously, we can immediately see that there is a certain size match that is necessary for the layers of a nested fullerene to interact in a stabilising fashion. The C20 fullerene is, for example, too large to fit into the cavity of the C60 fullerene. This results in a very large repulsive interaction energy. This repulsion is an

9

ACS Paragon Plus Environment

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

Page 10 of 30

Table 2: Interaction energies between individual CNO layers. Energies were calculated at the RI-PBE0-D3/def2-TZVP//RI-PBE0-D3/def2-SVP level of theory. C60 @C240 @C540 and C80 @C240 @C540 , were calculated at the RI-PBE0D3/def2-SVP level, energies in parentheses at RI-MP2/def2-TZVP//RI-PBE0D3/def2-SVP. All energies are in kcal mol−1 and interlayer distances in ˚ A. Cage 1 C20 C20 C60 C80 C60 C240 C60 @C240 C80 @C240

Cage 2 CNO C60 C20 @C60 C80 C20 @C80 C180 C60 @C180 C240 C80 @C240 C240 C60 @C240 C540 C240 @C540 C540 C60 @C240 @C540 C540 C80 @C240 @C540

Eint 1011.68 (933.27) 325.03 (93.95) 195.17 (-260.61) -81.90 (-511.55) -143.96 (-359.13) -489.62 -658.26 -670.18

distance 2.94 3.14 3.34 3.48 3.70 3.61 3.69, 3.61 3.49, 3.62

indication that while the C20 @C60 CNO is thermodynamically stable, it is unlikely to be formed because of the kinetic barrier that would have to be overcome. A similar situation can be observed when the C60 @C240 CNO is compared to C80 @C240 . Both CNOs interact attractively and their Eform differ in less than 0.1%, but the smaller C60 fullerene fits much better inside the cavity of C240 and hence the interaction energy of C60 @C240 is 43% larger than that of C80 @C240 . Secondly, larger fullerenes interact more strongly than smaller fullerenes in size-matched CNOs. For example, the interaction between C240 and C540 in C240 @C540 is more than three times larger than the interaction of C60 and C240 in C60 @C240 . This is the result of the fact that the dispersion interactions between the CNO layers scale with the size of the surface area of the constituent fullerenes. Thirdly and most importantly, the interactions between individual layers of CNOs are superadditive. This means that the overall interactions of all layers are larger than the sum of the individual interactions between separate layers. Table 2 shows that the interaction of C60 and C240 in C60 @C240 amounts to -143.96 kcal/mol and the interaction of C240 and C540 in C240 @C540 amounts to -489.62 kcal/mol. And yet, the interaction of the C60 @C240 CNO with its C540 outer shell in the C60 @C240 @C540 CNO amounts to -658.26 kcal/mol which

10

ACS Paragon Plus Environment

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

The Journal of Physical Chemistry

is 25 kcal/mol higher than the sum of the two individual interactions. The same is true, to an even greater extent, for the other triply layered CNO, C80 @C240 @C540 , as well. We believe that this cumulative effect explains the observed thermodynamic stability of CNOs with many layers. In effect, every additional layer of a growing CNO should be at least as stable as the previous layer, which indicates that the energetics of spherical CNO growth are somewhat analogue to the energetics of crystal growth in general although the mechanism is radically different. 76 In order to validate the DFT results, we also carried out correlated wave function calculations at the MP2 level of theory. The corresponding interaction energies are shown (in parentheses) in Table 2. Although different in detail, the overall trends of the MP2 results agree with the DFT calculations. The MP2 calculations indicate, in general, even stronger binding interactions than the DFT calculations. Hence, the MP2 results suggest that C60 @C180 has a stabilising interaction energy, whereas DFT predicts that the interaction in C60 @C180 is, in fact, destabilising. Whether this is an artefact of the well known tendency of MP2 to over-bind or a real property of the system is beyond the scope of this investigation. 77 Most importantly, the MP2 results agree with the overall conclusion that the larger CNOs are well stabilised. As far as the average distances between the individual CNO layers are concerned, it is apparent that all CNOs whose formation is predicted to be kinetically inhibited, i.e. C20 @C60 , C20 @C80 and C60 @C180 , exhibit shorter distances than the experimentally observed 3.35 ˚ A for a five-layered CNO. 15 This lends further credence to our prediction of the large kinetic barrier to their formation. The larger CNOs all show longer interlayer distances than experimentally predicted. We believe that this is mostly due to the relatively high uncertainty that is inherent in HRTEM measurements for materials such as CNOs. We have also calculated the diameters across (see Tables 3 and 4) and volumes enclosed by (see Table 5) the individual constituent fullerenes in CNOs. This allows us to examine the effects of the pressure build-up as subsequent layers of the CNOs form. The relatively large

11

ACS Paragon Plus Environment

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

Page 12 of 30

interaction energies (specifically those of the smallest CNOs) could potentially be accompanied by large structural changes. Surprisingly, none of the fullerene carbon cages experience much change during the formation of the CNOs. Neither do the inner fullerenes contract by a large amount, nor do the outer fullerenes expand much. Even in the high-energy systems C20 @C60 and C20 @C80 , the diameter of the C20 fullerene changes by a maximum of 1.5% (which causes a contraction of less than 6% by volume). The outer fullerenes, C60 and C80 respectively, undergo expansion of a similar magnitude. The changes to the larger CNOs are even smaller than the ones just discussed to the extent that the volume encompassed by the C540 moiety in C80 @C240 @C540 , is just a little more than 0.5% larger than the volume encompassed by a free C540 molecule. Table 3: Diameters of fullerene cages independently and inside a double layer CNO (all values in pm. Note that 1000 pm are 1 nm). Diameter C20 C60 C80 C180 C240 C540

406 709 819 1,225 1,413 2,107

C20 @C60 406 752 -

C20 @C80 400 834 -

Diameter in CNOs C60 @C180 C60 @C240 C80 @C240 696 710 814 1,239 1,411 1,418 -

C240 @C540 1,413 2,109

Table 4: Diameters of fullerene cages independently and inside a triple layer CNO (all values in pm. Note that 1000 pm are 1 nm). Diameter C60 C80 C240 C540

709 819 1,413 2,107

Diameter in CNOs C60 @C240 @C540 C80 @C240 @C540 710 817 1,412 1,419 2,110 2,111

12

ACS Paragon Plus Environment

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

The Journal of Physical Chemistry

Table 5: Volumes of free fullerenes and fullerene moieties inside CNOs. Relative volumes are normalized against C60 . The volume shown is for the moiety printed in bold. The units are bohr3 . fullerene/CNO C20 C20 @C60 C20 @C80 C60 C20 @C60 C60 @C180 C60 @C240 C80 C20 @C80 C80 @C240 C180 C60 @C180 C240 C60 @C240 C80 @C240 C540 C240 @C540 C240 @C540 C60 @C240 @C540 C60 @C240 @C540 C60 @C240 @C540 C80 @C240 @C540 C80 @C240 @C540 C80 @C240 @C540

Volume 23.61 23.20 22.23 161.66 192.11 153.63 162.70 258.01 272.88 254.06 915.50 949.56 1418.35 1414.59 1437.26 4791.47 1419.96 4805.12 162.43 1417.23 4806.79 256.03 1439.16 4821.04

13

Relative Volume 0.146 0.144 0.138 1.000 1.188 0.950 1.006 1.596 1.688 1.572 5.663 5.587 8.774 8.750 8.891 29.639 8.784 29.724 1.005 8.767 29.734 1.584 8.902 29.822

ACS Paragon Plus Environment

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

3.2

Internal Rotation of Carbon Nano Onions

Another interesting question about CNOs is whether the individual fullerene layers can rotate feely with respect to each other like the layers of a rolling-element bearing. Because of the non-covalent nature of the inter-layer interactions, the barrier to internal rotation is expected to be low, 38 at least for systems with reasonable inter-layer distances. To shed light on this issue, we performed a scan of the potential energy changes that are associated with the rotation of the inner fullerene with respect to the outer fullerene in C60 @C240 , C80 @C240 and C240 @C540 . To this effect, a series of structures with the inner fullerene rotated around an arbitrary axis in steps of 10 degrees were created (See Figure S1). It can be seen that the cavity inside the C240 moiety provides ample space for an almost free rotation of both C60 and C80 . The potential energy curve for the C60 Buckminster fullerene is essentially flat and the corresponding curve for C80 does not exhibit a much higher barrier. The curve describing the rotation of the C240 fullerene inside a C540 moiety however, does provide evidence for a substantial barrier to rotation. To ensure that our choice of axis didnt prejudice the outcome, the calculation for the C240 @C540 CNO was repeated for two more axes, perpendicular to the first one. The results are essentially identical and are shown in the supplementary information. Even when we consider the fact that the curves obtained in our calculations provides only an upper boundary to the true value, it appears reasonable to estimate that the barrier to rotation is approximately 30 kcal/mol (see Figure S2). Although more accurate calculations that allow individual structures to relax might resolve to a lower barrier, at this stage we feel confident to predict that there will be no observable rotation between these two fullerenes. Whether this finding extends to larger fullerenes in even bigger CNOs remains to be seen. However, given that the interaction energies appear to be superadditive, we can use this as an indication that the outer layers of larger CNOs might interact strongly enough to make rotation impossible. A related, but separate question is whether any gyroscopic motion would be spectroscop14

ACS Paragon Plus Environment

Page 14 of 30

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

The Journal of Physical Chemistry

ically observable. Traditionally, Raman spectroscopy has been the method of choice for the investigation of CNOs and the focus has been on distinguishing CNOs from other carbon nano structures 78–80 The experimental difficulties in producing pure carbon nano onions (i.e. free from other carbon nano compounds) with a well-defined monodispersity are often cited as reasons for the use of Raman spectroscopy over other methods. 81 All fullerenes in this investigation (and many fullerenes in general) are centrosymmetric and hence have no permanent dipole moment, which means that single-molecule IR spectroscopy is not possible for those species. This is reproduced (to within numerical accuracy) in our calculations (see Table S3) and the double-layer CNO structures we calculate show the same behaviour. The only exception is C80 @C240 , which has a small, non-zero permanent dipole moment. Consequently, as far as the relative rotation of two objects with no (or a very small) permanent dipole moment in objects that themselves have no (or a very small) permanent dipole moment is concerned, it is a priori unclear if the breaking of the symmetry of any CNO through rotation of one fullerene with respect to the other in a CNO is enough to produce a change in the dipole moment. We have computed the permanent dipole moment of all structures used to create the potential energy scans reported in Table S4. From this data, it can be concluded that neither the rotation of C60 inside C240 , nor the rotation of C240 inside C540 change the dipole moment of the CNO to any significant extent. The rotation of C80 inside C240 however, shows variations of up to 0.44 Debye. While this value may only serve as an upper boundary to the actual changes to the permanent dipole moment upon rotation, it gives an indication that this motion might be able to be spectroscopically probed in an IR experiment.

3.3

Electronic Structure of Carbon Nano Onions

To gain some insight into the electronic structure of CNOs, we compare the frontier orbitals of the constituent fullerenes to the frontier orbitals of the resulting CNOs (see Figure 4). Inspection of the frontier orbital plots immediately reveals that while the individual orbitals 15

ACS Paragon Plus Environment

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

adjust their energy in the electrostatic potential the other CNO component creates, very little orbital mixing between the constituent parts of the CNO is happening. This is in a sense unsurprising, since no covalent interaction between individual CNO layers was expected. This also means however, that the molecular orbital scheme of CNOs consists of a sequence of states that are either localised on the inner or the outer fullerene, but not delocalised over both. Those cases where the localisation of the HOMO and LUMO alternates between the inner and outer fullerene are promising because they indicate that charge-transfer excited states might be possible within in the material as well as through the stabilisation of charge transfer states in donor-acceptor heterojunctions.

Figure 4: Frontier orbitals of fullerenes and CNOs calculated at PBE0-D3/def2-SVP. All energies are given in kcal mol−1 . Recently, Sol`a and Voityuk published an in-depth analysis of excited states in the C60 @C240 16

ACS Paragon Plus Environment

Page 16 of 30

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

The Journal of Physical Chemistry

CNO that correspond to charge-transfer between the inner and outer shells. 82 They found that the structure of CNOs stabilises charge-transfer states where one electron is promoted from the inner to the outer Fullerene forming C60+ @C240 – .

Figure 5: UV-Vis spectra of the C60 , C240 and C540 Fullerenes and the C60 @C240 , C240 @C540 and C60 @C240 @C540 CNOs. The black lines indicate the frequencies of individual excitations and their height is proportional to their oscillator stength. The red line indicates the chargetransfer excitation and the insets show the molecular orbitals involved in that transition. Note that the spectra of C240 @C540 and C60 @C240 @C540 are truncated at approximately 480 nm because of the limitation to 256 excited states in the ZINDO/S calculations. Since they were able to demonstrate the reliability of the ZINDO method for the description of CNO excitations, we have adopted the same methodology and applied it to the CNOs in this investigation. The simulated UV-Vis spectra that result from these calculations are shown in Figure 5. It can be seen that the optical absorption is more red-shifted 17

ACS Paragon Plus Environment

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

the larger the Fullerene becomes (left panels). It is also obvious that the main features of the individual spectra transfer directly to the spectra of the CNOs (right panels). For example: The strong absorptions around 450 nm and 380 nm in the C240 spectrum can be found in the spectrum of the C60 @C240 CNO as well. However, two distinct differences can be observed. Firstly, it becomes clear that there is some dispersion of the excitation energies. For example, the strong absorption in the C240 at 542 nm shifts to 473 nm in the spectrum of C60 @C240 . Secondly, many new absorptions that are present in neither the C60 nor the C240 spectrum appear. This is partly due to the fact that some formerly forbidden transitions become allowed in the CNO of reduced symmetry, and partly because new excitations from states on C60 to C240 (and vice-versa) appear in the spectrum. This aspect of the C60 @C240 spectrum has been discussed in detail by Sol`a. 82 The final two panels that show the UV-Vis spectra of C240 @C540 and C60 @C240 @C540 confirm that this is also the case for these two - much larger - carbon nano onions. It is indeed evident from our data that the lowest charge-transfer excitation (Figure 5, shown in red) is lower in energy than any transition of the individual Fullerenes that constitute the two CNOs in question.

4

Conclusions

We have investigated double-layered carbon nano-onions as well as two triple-layered analogues, C60 @C240 @C540 , and C80 @C240 @C540 , that can be obtained by stacking two doublelayered CNOs. Using high-level electronic structure calculations, we were able to show that the free rotation of the inner Fullerene inside a CNO against the outer shell depends on sizematching of the two adjacent pairs and is furthermore restricted to the innermost layer. We have found that the strong dispersion interactions between individual layers of CNOs are superadditive. This means that the total interaction energy of three shells was larger than the sum of the two individual interaction energies of the two double-layered CNOs that comprise

18

ACS Paragon Plus Environment

Page 18 of 30

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

The Journal of Physical Chemistry

the three-layered CNO. This result was found for both triple-layered CNOs in our investigation. This means that the previously reported free rotation of the inner Fullerene in CNOs is unlikely to occur in the outer shells of CNOs with many layers. We furthermore conjecture that this superadditivity is important for the experimentally observed growth of CNOs with many layers. We also show that the molecular orbitals of individual Fullerenes do not mix when CNOs are formed, but the energy levels of these molecular orbitals shift in energy according to their new environment. This is a consequence of the fact that the interaction between individual layers of CNOs is neither covalent nor ionic, but dominated by dispersion interactions. Our excited-state calculations confirm that this finding extends to the UV-vis spectra of CNOs, which comprises the spectra of individual Fullerenes energy shifted and overlaid on top of each other. The excited states that correspond to charge-transfer between different layers of a particular CNO are particularly interesting for applications. We we also able to show that the lowest such transition occurs at much lower energies in C60 @C240 and C60 @C240 @C540 than previously reported for C60 @C240 .

Acknowledgement The authors thank Julia Schacht, Yasir Altaf, Dr. Renee Goreham, Dr. Siobhan Bradley, and Dr. Khurshid Ayub for useful discussions. MAH acknowledges Victoria University of Wellington for a Victoria Doctoral Scholarship. Additional computer time was provided by the Victoria University of Wellington High Performance Computer Facility SciFacHPC and by the New Zealand eScience Infrastructure (NeSI) project NESI181 through the Pan cluster at the Centre for eResearch at the University of Auckland.

5

Supplementary Information

Optimized geometries of all structures and animations of the internal rotations are available online. 19

ACS Paragon Plus Environment

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

6

Conflict of Interest

The authors declare that they have no conflict of interest.

Notes and References (1) Suarez-Martinez, I.; Grobert, N.; Ewels, C. P. Nomenclature of sp2 Carbon Nanoforms. Carbon 2012, 50, 741 – 747. (2) Iijima, S. Direct Observation of the Tetrahedral Bonding in Graphitized Carbon Black by High Resolution Electron Microscopy. J. Cryst. Growth 1980, 50, 675 – 683. (3) Kroto, H. W.; Heath, J. R.; O’Brien, S. C.; Curl, R. F.; Smalley, R. E. C60 : Buckminsterfullerene. Nature 1985, 318, 162–163. (4) Ugarte, D. Curling and Closure of Graphitic Networks under Electron-Beam Irradiation. Nature 1992, 359, 707–709. (5) de Heer, W. A.; Ugarte, D. Carbon Onions Produced by Heat Treatment of Carbon Soot and their Relation to the 217.5 nm Interstellar Absorption Feature. Chem. Phys. Lett. 1993, 207, 480–486. (6) Golberg, D.; Bando, Y.; Bourgeois, L.; Kurashima, K. Atomic Resolution of SingleWalled Carbon Nanotubes using a Field Emission High-Resolution Transmission Electron Microscope. Carbon 1999, 37, 1858–1860. (7) Mordkovich, V. Z.; Umnov, A. G.; Inoshita, T.; Endo, M. The Observation of Multiwall Fullerenes in Thermally Treated Laser Pyrolysis Carbon Blacks. Carbon 1999, 37, 1855–1858. (8) Sano, N.; Wang, H.; Chhowalla, M.; Alexandrou, I.; Amaratunga, G. A. J. Nanotechnology: Synthesis of Carbon ’Onions’ in Water. Nature 2001, 414, 506–507. 20

ACS Paragon Plus Environment

Page 20 of 30

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

The Journal of Physical Chemistry

(9) Rettenbacher, A. S.; Elliott, B.; Hudson, J. S.; Amirkhanian, A.; Echegoyen, L. Preparation and Functionalization of Multilayer Fullerenes (Carbon Nano-Onions). Chem. Eur. J. 2006, 12, 376–387. (10) Zou, Q.; Wang, M. Z.; Li, Y. G. Onion-Like Carbon Synthesis by Annealing Nanodiamond at Lower Temperature and Vacuum. J. Exp. Nanosci. 2010, 5, 375–382. (11) Terrones, M.; Terrones, G.; Terrones, H. In Science of Crystal Structures: Highlights in Crystallography; Hargittai, I., Hargittai, B., Eds.; Springer International Publishing: Cham, 2015; pp 101–112. (12) Srolovitz, D. J.; Safran, S. A.; Homyonfer, M.; Tenne, R. Morphology of Nested Fullerenes. Phys. Rev. Lett. 1995, 74, 1779–1782. (13) Zheng, Y.; Zhu, P. Carbon Nano-Onions: Large-Scale Preparation, Functionalization and their Application as Anode Material for Rechargeable Lithium Ion batteries. RSC Adv. 2016, 6, 92285–92298. (14) Tomita, S.; Sakurai, T.; Ohta, H.; Fujii, M.; Hayashi, S. Structure and Electronic Properties of Carbon Onions. J. Chem. Phys. 2001, 114, 7477–7482. (15) Zou, Q.; Li, Y. G.; Wang, M. Z.; Zou, L. H.; Zhao, Y. C. Characterisation of Onion-Like Carbon Fabricated by Annealing Nanodiamond. Mat. Res. Innov. 2010, 14, 285–288. (16) Pichler, T.; Knupfer, M.; Golden, M. S.; Fink, J.; Cabioc’h, T. Electronic Structure and Optical Properties of Concentric-Shell Fullerenes from Electron-Energy-Loss Spectroscopy in Transmission. Phys. Rev. B 2001, 63, 155415. (17) Ryu, Z.; Zheng, J.; Wang, M.; Zhang, B. Synthesis and Characterization of Silicon Carbide Whiskers. Carbon 2001, 39, 1929–1930. (18) Zhang, C.; Li, J.; Liu, E.; He, C.; Shi, C.; Du, X.; Hauge, R. H.; Zhao, N. Synthesis

21

ACS Paragon Plus Environment

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

Page 22 of 30

of Hollow Carbon Nano-Onions and their Use for Electrochemical Hydrogen storage. Carbon 2012, 50, 3513–3521. (19) Wang, Q.; Sun, X.; He, D.; Zhang, J. Preparation and Study of Carbon Nano-Onion for Lithium Storage. Mater. Chem. Phys. 2013, 139, 333–337. (20) Camisasca, A.; Giordani, S. Carbon Nano-Onions in Biomedical Applications: Promising Theranostic Agents. Inorg. Chim. Acta 2017, 468, 67–76. (21) Matsumoto, N.; Joly-Pottuz, L.; Kinoshita, H.; Ohmae, N. Application of Onion-Like Carbon to Micro and Nanotribology. Diam. Relat. Mater. 2007, 16, 1227–1230. (22) Guo, J.; Yang, X.; Yao, Y.; Wang, X.; Liu, X.; Xu, B. Pt/Onion-Like Fullerenes as Catalyst for Direct Methanol Fuel Cell. Rare Metals 2006, 25, 305–308. (23) J¨ackel, N.; Weingarth, D.; Zeiger, M.; Aslan, M.; Grobelsek, I.; Presser, V. Comparison of Carbon Onions and Carbon Blacks as Conductive Additives for Carbon Supercapacitors in Organic Electrolytes. J. Power Sources 2014, 272, 1122–1133. (24) Pech, D.; Brunet, M.; Durou, H.; Huang, P.; Mochalin, V.; Gogotsi, Y.; Taberna, P.L.; Simon, P. Ultrahigh-Power Micrometre-Sized Supercapacitors based on Onion-Like carbon. Nat. Nano 2010, 5, 651–654. (25) Lettieri, S.; Camisasca, A.; d’Amora, M.; Diaspro, A.; Uchida, T.; Nakajima, Y.; Yanagisawa, K.; Maekawa, T.; Giordani, S. Far-Red Fluorescent Carbon Nano-Onions as a Biocompatible Platform for Cellular Imaging. RSC Adv. 2017, 7, 45676–45681. (26) Marchesano, V.;

Ambrosone, A.;

Bartelmess, J.;

Strisciante, F.;

Tino, A.;

Echegoyen, L.; Tortiglione, C.; Giordani, S. Impact of Carbon Nano-Onions on Hydra Vulgaris as a Model Organism for Nanoecotoxicology. Nanomaterials 2015, 5, 1331–1350.

22

ACS Paragon Plus Environment

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

The Journal of Physical Chemistry

(27) Palanker, D. V. Photovoltaic Restoration of Sight in Rodents with Retinal Degeneration (Conference Presentation). Proc. SPIE 2017, 10052, 100520I. (28) Mykhailiv, O.; Zubyk, H.; Plonska-Brzezinska, M. E. Carbon Nano-Onions: Unique Carbon Nanostructures with Fascinating Properties and their Potential Applications. Inorg. Chim. Acta 2017, 468, 49–66. (29) Kuznetsov, V. L.; Chuvilin, A. L.; Butenko, Y. V.; Mal’kov, I. Y.; Titov, V. M. OnionLike Carbon from Ultra-Disperse Diamond. Chem. Phys. Lett. 1994, 222, 343 – 348. (30) Heggie, M. I.; Terrones, M.; Eggen, B. R.; Jungnickel, G.; Jones, R.; Latham, C. D.; Briddon, P. R.; Terrones, H. Quantitative Density-Functional Study of Nested Fullerenes. Phys. Rev. B 1998, 57, 13339–13342. (31) T¨ urker, L. A Bucky Onion from C20 and C60 —an AM1 Treatment. J. Mol. Struct. THEOCHEM 2001, 545, 207–214. (32) T¨ urker, L.; G¨ um¨ u¸s, S. An AM1 Study on C60 @C180 system. J. Mol. Struct. THEOCHEM 2004, 674, 15–18. (33) Liu, F.; Meng, L.; Zheng, S. Density Functional Studies on a Novel Double-Shell Fullerene C20 @C60 . J. Mol. Struct. THEOCHEM 2005, 725, 17–21. (34) Mordkovich, V. Z. The Observation of Large Concentric Shell Fullerenes and Fullerenelike Nanoparticles in Laser Pyrolysis Carbon Blacks. Chem. Mater. 2000, 12, 2813– 2818. (35) Zope, R. R.; Dunlap, B. I. Slater’s Exchange Parameters α for Analytic and Variational Xα Calculations. J. Chem. Theor. Comput. 2005, 1, 1193–1200. (36) Rajendra, R. Z. Electronic Structure and Static Dipole Polarizability of C60 @C240 . J. Phys. B 2008, 41, 085101.

23

ACS Paragon Plus Environment

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

(37) Glukhova, O. E.; Kolesnikova, A. S.; Slepchenkov, M. M.; Shunaev, V. V. Moving of Fullerene between Potential Wells in the External Icosahedral Shell. J. Comp. Chem. 2014, 35, 1270–1277. (38) Zope, R. R.; Bhusal, S.; Basurto, L.; Baruah, T.; Jackson, K. Site Specific Atomic Polarizabilities in Endohedral Fullerenes and Carbon Onions. J. Chem. Phys. 2015, 143, 084306. (39) Perdew, J. P.; Burke, K.; Ernzerhof, M. Generalized Gradient Approximation Made Simple. Phys. Rev. Lett. 1996, 77, 3865–3868. (40) Perdew, J. P.; Burke, K.; Ernzerhof, M. Errata: Generalized Gradient Approximation Made Simple. Phys. Rev. Lett. 1997, 78, 1396. (41) Adamo, C.; Barone, V. Toward Reliable Density Functional Methods Without Adjustable Parameters: The PBE0 Model. J. Chem. Phys. 1999, 110, 6158–6169. (42) Grimme, S.; Antony, J.; Ehrlich, S.; Krieg, H. A Consistent and Accurate Ab Initio Parameterization of Density Functional Dispersion Correction (DFT-D) for the 94 Elements H-Pu. J. Chem. Phys. 2010, 132, 154104. (43) Pinski, P.; Riplinger, C.; Valeev, E. F.; Neese, F. Sparse Maps–A Systematic Infrastructure for Reduced-Scaling Electronic Structure Methods. I. An Efficient and Simple Linear Scaling Local MP2 Method that uses an Intermediate Basis of Pair Natural Orbitals. J. Chem. Phys. 2015, 143, 034108. (44) Riplinger, C.; Pinski, P.; Becker, U.; Valeev, E. F.; Neese, F. Sparse Maps–A Systematic Infrastructure for Reduced-Scaling Electronic Structure Methods. II. Linear Scaling Domain Based Pair Natural Orbital Coupled Cluster Theory. J. Chem. Phys. 2015, 144, 024109.

24

ACS Paragon Plus Environment

Page 24 of 30

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

The Journal of Physical Chemistry

(45) Guo, Y.; Sivalingam, K.; Valeev, E. F.; Neese, F. SparseMaps–A Systematic Infrastructure for Reduced-Scaling Electronic Structure Methods. III. Linear-Scaling Multireference Domain-Based Pair Natural Orbital N-Electron Valence Perturbation Theory. J. Chem. Phys. 2015, 144, 094111. (46) Pavoevi, F.; Pinski, P.; Riplinger, C.; Neese, F.; Valeev, E. F. SparseMaps–A Systematic Infrastructure for Reduced-Scaling Electronic Structure Methods. IV. Linear-Scaling Second-Order Explicitly Correlated Energy with Pair Natural Orbitals. J. Chem. Phys. 2015, 144, 144109. (47) Dunlap, B. I.; Connolly, J. W. D.; Sabin, J. R. On some Approximations in Applications of Xα Theory. J. Chem. Phys. 1979, 71, 3396–3402. (48) Vahtras, O.; Almlf, J.; Feyereisen, M. W. Integral Approximations for LCAO-SCF Calculations. Chem. Phys. Lett. 1993, 213, 514–518. (49) Neese, F. An Improvement of the Resolution of the Identity Approximation for the Formation of the Coulomb Matrix. J. Comp. Chem. 2003, 24, 1740–1747. (50) Izsak, R.; Neese, F. An Overlap Fitted Chain of Spheres Exchange Method. J. Chem. Phys. 2011, 135, 144105. (51) Kossmann, S.; Neese, F. Efficient Structure Optimization with Second-Order ManyBody Perturbation Theory: The RIJCOSX-MP2 Method. J. Chem. Theory Comput. 2010, 6, 2325–2338. (52) Neese, F.; Wennmohs, F.; Hansen, A.; Becker, U. Efficient, Approximate and Parallel Hartree–Fock and Hybrid DFT Calculations. A ‘Chain-of-Spheres’ Algorithm for the Hartree–Fock Exchange. Chem. Phys. 2009, 356, 98–109. (53) Weigend, F. Accurate Coulomb-fitting basis sets for H to Rn. Phys. Chem. Chem. Phys. 2006, 8, 1057–1065. 25

ACS Paragon Plus Environment

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

(54) Dunning Jr., T. H. Gaussian Basis Sets for use in Correlated Molecular Calculations. I. The Atoms Boron through Neon and Hydrogen. J. Chem. Phys. 1989, 90, 1007–1023. (55) Kendall, R. A.; Dunning Jr., T. H.; Harrison, R. J. Electron Affinities of the First-Row Atoms Revisited. Systematic Basis Sets and Wave Functions. J. Chem. Phys. 1992, 96, 6796–6806. (56) Woon, D. E.; Dunning Jr., T. H. Gaussian-Basis Sets for Use in Correlated Molecular Calculations. 3. The Atoms Aluminum through Argon. J. Chem. Phys. 1993, 98, 1358– 1371. (57) Peterson, K. A.; Woon, D. E.; Dunning Jr., T. H. Benchmark Calculations with Correlated Molecular Wave Functions. IV. The Classical Barrier Height of the H+H2 → H2 +H reaction. J. Chem. Phys. 1994, 100 . (58) Wilson, A. K.; van Mourik, T.; Dunning Jr., T. H. Gaussian Basis Sets for Use in Correlated Molecular Calculations. VI. Sextuple Zeta Correlation Consistent Basis Sets for Boron Through Neon. J. Mol. Struct. THEOCHEM 1996, 388, 339–349. (59) Davidson, E. R. Comment on ”Comment on Dunning’s Correlation-Consistent Basis Sets”. Chem. Phys. Lett. 1996, 260, 514–518. (60) Ridley, J.; Zerner, M. An Intermediate Neglect of Differential Overlap Technique for Spectroscopy: Pyrrole and the Azines. Theor. Chem. Acc. 1973, 32, 111–134. (61) Bacon, A. D.; Zerner, M. C. An Intermediate Neglect of Differential Overlap Theory for Transition Metal Complexes: Fe, Co and Cu Chlorides. Theor. Chem. Acc 1979, 53, 21–54. (62) Cory, M. G.; Zerner, M. C. Metal-Ligand Exchange Coupling in Transition-Metal Complexes. Chem. Rev. 1991, 91, 813–822.

26

ACS Paragon Plus Environment

Page 26 of 30

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

The Journal of Physical Chemistry

(63) Kotzian, M.; R¨osch, N.; Zerner, M. C. Intermediate Neglect of Differential Overlap Spectroscopic Studies on Lanthanide Complexes. Theor. Chem. Acc 1992, 81, 201– 222. (64) Neese, F. The ORCA Program System. Wiley Interdiscip. Rev.: Comput. Mol. Sci. 2012, 2, 73–78. (65) Dennington, R.; Keith, T.; Millam, J. GaussView Version 5.0.9. Semichem Inc., Shawnee Mission, KS, 2009. (66) 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. et al. Gaussian 09 Revision D.01. Gaussian Inc. Wallingford CT, 2009. (67) Legault, C. Y. CYLview. 1.0b. Universit´e de Sherbrooke, http://www.cylview.org, 2009. (68) Mordkovich, V. Z.; Umnov, A. G.; Inoshita, T. Nanostructure of Laser Pyrolysis Carbon Blacks: Observation of Multiwall Fullerenes. Int. J. Inorg. Mat. 2000, 2, 347–353. (69) Mordkovich, V.; Shiratori, Y.; Hiraoka, H.; Umnov, A.; Takeuchi, Y. A path to larger yields of multishell fullerenes. Carbon 2001, 39, 1938–1941. (70) Welsh, I.; Lein, M. Accurate Density Functional Theory Description of Binding Constants and NMR Chemical Shifts of Weakly Interacting Complexes of C60 with Corannulene-Based Molecular Bowls. J. Comp. Chem. 2014, 35, 181–191. (71) Kleingeld, R.; Wang, J.; Lein, M. Buckminster Fullerene Adhesion on Graphene Flakes: Numerical Accuracy of Dispersion Corrected DFT. Polyhedron 2016, 114, 110–117. (72) Goerigk, L.; Grimme, S. A Thorough Benchmark of Density Functional Methods for General Main Group Thermochemistry, Kinetics, and Noncovalent Interactions. Phys. Chem. Chem. Phys. 2011, 13, 6670–6688.

27

ACS Paragon Plus Environment

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

(73) Li, A.; Muddana, H. S.; Gilson, M. K. Quantum Mechanical Calculation of Noncovalent Interactions: A Large-Scale Evaluation of PMx, DFT, and SAPT Approaches. J. Chem. Theor. Comput. 2014, 10, 1563–1575. (74) Mar, N.; Sansores, L. E.; Muhl, S.; Ramos, E.; Salcedo, R. Symmetric Nested Complexes of Fullerenes. J. Mol. Model. 2015, 21, 1–7. (75) Tomanek, D.; Frederick, N. Guide through the Nanocarbon Jungle: Buckyballs, Nanotubes, Graphene, and Beyond. http://www.nanotube.msu.edu/fullerene/ fullerene-isomers.html. (76) Irle, S.; Zheng, G.; Wang, Z.; Morokuma, K. The C60 Formation Puzzle “Solved”: QM/MD Simulations Reveal the Shrinking Hot Giant Road of the Dynamic Fullerene Self-Assembly Mechanism. J. Phys. Chem. B 2006, 110, 14531–14545. (77) Preliminary SCS-MP2 calculations (which potentially mitigate against the overbinding nature of unscaled MP2 calculations) yield lower interaction energies. (78) Roy, D.; Chhowalla, M.; Wang, H.; Sano, N.; Alexandrou, I.; Clyne, T.; Amaratunga, G. Characterisation of Carbon Nano-Onions using Raman Spectroscopy. Chem. Phys. Lett. 2003, 373, 52 – 56. (79) Wang, X.; Xu, B.; Liu, X.; Jia, H.; Hideki, I. The Raman Spectrum of Nano-Structured Onion-Like Fullerenes. Physica B: Condens. Matter 2005, 357, 277 – 281. (80) Codorniu Pujals, D.; Arias de Fuentes, O.; Desd´ın Garc´ıa, L. F.; Cazzanelli, E.; Caputi, L. S. Raman Spectroscopy of Polyhedral Carbon Nano-Onions. Appl. Phys. A 2015, 120, 1339–1345. (81) Obraztsova, E.; Fujii, M.; Hayashi, S.; Kuznetsov, V.; Butenko, Y.; Chuvilin, A. Raman Identification of Onion-Like Carbon. Carbon 1998, 36, 821 – 826.

28

ACS Paragon Plus Environment

Page 28 of 30

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

The Journal of Physical Chemistry

(82) Voityuk, A. A.; Sol`a, M. Photoinduced Charge Separation in the Carbon Nano-Onion C60@C240. J. Phys. Chem. A 2016, 120, 5798–5804.

29

ACS Paragon Plus Environment

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

TOC Graphic

30

ACS Paragon Plus Environment

Page 30 of 30