The Correlated Electronic States of a Few Polycyclic Aromatic

6 May 2019 - In recent years Polycyclic Aromatic Hydrocarbons (PAHs) have been studied for their electronic properties as they are viewed as nanodots ...
0 downloads 0 Views 402KB Size
Subscriber access provided by AUT Library

A: Spectroscopy, Molecular Structure, and Quantum Chemistry

The Correlated Electronic States of a Few Polycyclic Aromatic Hydrocarbons: A Computational Study Geetanjali Giri, Anusooya Yadthare Pati, and Sooryanarayana Ramasesha J. Phys. Chem. A, Just Accepted Manuscript • DOI: 10.1021/acs.jpca.9b02196 • Publication Date (Web): 06 May 2019 Downloaded from http://pubs.acs.org on May 7, 2019

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 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 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.

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 73 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 Correlated Electronic States of a few Polycyclic Aromatic Hydrocarbons: A Computational Study Geetanjali Giri, Y. Anusooya Pati, and S. Ramasesha∗ Solid State and Structural Chemistry Unit, Indian Institute of Science, Bangalore 560 012, India E-mail: [email protected]

∗ To

whom correspondence should be addressed

ACS Paragon Plus Environment

1

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 In recent years Polycyclic Aromatic Hydrocarbons (PAHs) have been studied for their electronic properties as they are viewed as nanodots of graphene. They have also been of interest as functional molecules for applications such as light emitting diodes and solar cells. Since last few years varying structural and chemical properties corresponding to the size and geometry of these molecules have been studied both theoretically and experimentally. In this paper, we carry out a systematic study of the electronic states of several PAHs using the Pariser-Parr-Pople model which incorporates long-range electron correlations. In all the molecules studied by us, we find that the 2A state is below the 1B state and hence none of them will be fluorescent in the gaseous phase. The singlet-triplet gap is more than one-half of the singlet-singlet gap in all cases and hence none of these PAHs can be candidates for improved solar cell efficiencies in a singlet fission. We discuss in detail the properties of the electronic states which include bond orders and spin densities (in triplets) of these systems.

Introduction Electronic structure of π− conjugated carbon systems is of enduring interest, although there have been paradigm shifts in the underlying reasons. Early interest revolved around aromaticity 1–3 and later with the advent of the era of electronic polymers, the interest has shifted to the study of relative energy level ordering of the low energy excited states. With recent focus on graphene, there has been a resurgence of interest in condensed aromatic ring systems. Thus, polycyclic aromatic hydrocarbons (PAHs) which were studied mainly for environmental impact are now of interest from the electronic structure view point as these molecules serve as model nanodots of graphene. 4–6 The effect of shape on the optical properties of nanodots of graphene, have also been studied. 7 The role of size, shape and hetero atom doping on the magnetic and electronic properties of graphene quantum dots have been studied theoretically using density functional methods. 8–10 Indeed, it has also been reported that alkali metal doped ACS Paragon Plus Environment

2

Page 2 of 73

Page 3 of 73 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

picene exhibits superconductivity below 18K 11 and K-doped 1,2:8,9-dibenzopentacene shows superconductivity above 30 K, 12 Sm-doped phenacenes show very low Tc ≈ 5K, 13,14 which adds further dimension to the study of PAHs. The spectra of PAHs have also been studied extensively, since they have been detected in the interstellar medium. 15 Experimental studies have mainly focused on the vibrational spectra of PAH molecules. Allamandola et al have assigned the unidentified infrared bands (UIBs) of interstellar particles to Raman spectra of partially hydrogenated, positively charged PAHs. 16 They have studied the infrared fluorescence spectrum of chrysene molecule. The observation of diffused interstellar absorption bands (DIBs) is attributed to the presence of larger PAHs molecule. 17 Salama et al have obtained the spectra of several neutral and ionic PAH molecules under the astrophysical laboratory conditions and compared this with the astronomical spectra. 18 These astronomical bands can be explained with electronic absorptions and/or dynamical methods such as electronic relaxations and intramolecular vibrational relaxations. Palewska et al have obtained the high resolution fluorescence and phosphorescence spectra of tetrahelicene and hexahelicene in nheptane matrix at room temperature. 19 They find the singlet-singlet transition to be at 3.1 eV whereas the singlet-triplet transition to be at 2.5 eV. The electronic absorption spectrum of chrysene molecule in boric acid glass matrix is reported by Hussain. 20 They compared their excitation energies obtained from experiments with those obtained from semiempirical AM1 method. They found a reasonably good agreement in the calculated bond lengths and bond angles with those obtained from the X-ray diffraction. Electronic absorption spectra in complex PAH molecules were obtained within semiempirical approach where restricted configuration interactions were used. 21 Ping Du et al used semi empirical quantum mechanical intermediate neglect of differential overlap (INDO) method to compare optical spectra of naphthalene and its ions with experimental data. 22 Infrared spectra of neutral and cationic PAHs are obtained within the DFT method and the neutral molecules showed a good agreement with the experimental spectra. 23 The absolute absorption crossACS Paragon Plus Environment

3

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

section of several PAH molecules and their cations were obtained within time-dependent density functional theory. 24 They find good agreement with the experimental excitation energies for neutral PAH molecules. This is not surprising as TDDFT method with adiabatic local density approximation is equivalent to single CI method 25 which has been shown to give good optical gaps due to cancellation of errors. However, with ions, this cancellation does not occur. Using time dependent DFT, Hammonds et al have obtained the electronic spectra of protonated and partially hydrogenated PAHs and found similarity in the absorption peaks of diffuse interstellar bands. 26 Size and geometry dependence of larger PAH molecules on optical spectra has been studied by Cocchi et al within ZINDO model, using only single excitation configuration interactions (SCI). 27 The role of aliphatic side groups on the vibrational spectra of PAHs molecule is studied using the DFT method. 28 A feature common to all conjugated carbon systems is that the π−electrons experience strong correlations. The role of electron correlation on the ordering of the excited state energy levels in PAHs has been studied by Mazumdar et al using Pariser-Parr-Pople model hamiltonian. 29 They find that the effect of electron-electron interactions is weaker in PAHs than in polyenes, but stronger than in the quasi-1D systems like poly-(paraphenylene) and poly(paraphenylenevinylene). Basak et al have carried out a detailed study of the linear absorption spectra of diamond shaped graphene quantum dots (where pyrene is the smallest unit), with quadruple CI and MRSDCI basis using the correlated PPP model Hamiltonian. They have compared their excitation gaps with screened and standard PPP parameters for carbon. 30 They have also studied the two-photon absorption and photoinduced absorption of these molecules. They find that the effect of electron-electron correlation gets enhanced with the increase in the size of graphene quantum dots. 31 Motivated by this large number of experimental and theoretical data on PAHs, the current study focuses on the effect of strong electron-electron correlations on one-photon and two-photon absorptions in PAHs. So far, a thorough study of this class of molecules inclusive of strong long range electron-electron correlations has not been reported to ACS Paragon Plus Environment

4

Page 4 of 73

Page 5 of 73 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 best of our knowledge. In this paper, we have compared electronic spectra of various PAHs (see Fig. 1) such as Fluoranthene, Pyrene (both 16 conjugated carbon atoms), Benzanthracene, Chrysene, Helicene, Triphenylene (18 conjugated carbon atoms), Benzopyrene (20 conjugated carbon atoms) and Picene (22 conjugated carbon atoms) using numerically exact diagrammatic valence bond (DVB) method as well as Density matrix renormalization group method (DMRG) for larger molecules. We have analyzed the bond orders and spin densities in triplet states based on the ring orientations of these systems. We find that none of the molecules belonging to this class, studied by us, fulfill the energy criterion for singlet fission. Details of methods used in this study are discussed in next section followed by a description of singlet and triplet states of all molecules considered, with their corresponding properties. We have summarized our findings in the last section.

Methodology We have considered all molecules to be geometrically planar. Ignoring hydrogen atoms in a molecule leaves us with the carbon skeleton. Carbon atoms in sp2 hybridization in aromatic hydrocarbons provides three sigma bonds with neighboring atoms resulting in a rigid framework of the system whereas 2pz orbital gives conjugated network perpendicular to the plane of the molecule. This conjugated framework in aromatic systems gives rise to many interesting properties owing to the delocalized nature of orbitals. In our studies we focus on π-electron network to explore various properties of the PAHs as this network is energetically well separated from the underlying σ−network. This approximation is the σ − π separability approximation. 32 In the present study, we have used Pariser-Parr-Pople (PPP) model Hamiltonian to describe the π−electrons. 33,34 PPP model is obtained from many-body Hamiltonian considering only πelectrons and invoking zero differential overlap (ZDO) approximation. The PPP Hamiltonian is ACS Paragon Plus Environment

5

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

[1]

[2]

[4]

Page 6 of 73

[3]

[5]

[6]

[8] [7]

Figure 1: Schematic representation of PAH molecules that have been studied. 1) Fluoranthene, 2) Pyrene, 3) Benzanthracene, 4) Chrysene, 5) Helicene, 6) Triphenylene, 7) Benzopyrene and 8) Picene.

ACS Paragon Plus Environment

6

Page 7 of 73 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

given below:

Hˆ PPP =





+∑ i

t0 (cˆ†i,σ cˆ j,σ + H.C.) + ∑ εi nˆ i i

U nˆ i (nˆ i − 1) + ∑ Vi j (nˆ i − zi )(nˆ j − z j ) 2 i> j

(1)

Here first term is the noninteracting part of the PPP Hamiltonian in which t0 is nearest-neighbour hopping term between two bonded orbitals ‘i’ and ‘j’. cˆ†i,σ (cˆi,σ ) creates (annihilates) an electron of spin σ in 2pz orbital of carbon atom at site i. εi is the site energy of ‘i’-th carbon atom. The second term corresponds to interacting part of the Hamiltonian which consists of on-site Hubbard repulsion term for two electrons occupying the same orbital. In the third term Vi j s are electron-electron repulsion integrals between electrons in two separate orbitals i and j. Vi j s are e2 interpolated between U and as ri j → ∞ using the Ohno interpolation scheme, 35 r

" Vi j = 14.397

28.794 (Ui +U j )

2

#− 1 2 2

+ ri j

(2)

The operator nˆi is the number operator, the C-C bond length is fixed at 1.4 A˚ for nearest neighbour carbon atoms. We use the standard PPP model parameters for carbon, namely, t0 = −2.4 eV, U = 11.26 eV and site energies εi = 0.0 eV for unsubstituted carbon site. zi is the number of electrons at site ‘i’ which leaves the site neutral and correspond to local chemical potential. For carbon atoms in conjugation, zi values are set to 1. The PPP model Hamiltonian conserves total spin. So, to obtain the eigenstates of PAHs within PPP model we have used the diagrammatic valence bond (DVB) method because it generates spin-adapted basis functions contrary to slater basis functions which conserves only the z-component of total spin. To reduce the computational cost further, we have exploited the combination of C2 and electron-hole (e-h) symmetry and factorized the Hilbert space into symACS Paragon Plus Environment

7

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

metry adapted bases. 36,37 Since C2 and eh symmetry operators commute, we have an Abelian group with 4 one-dimensional representation. Except benzanthracene, all molecules studied by us have C2 symmetry - either along X-axis or along Y-axis. Since these molecules are assumed to be planar, inversion symmetry and C2 along Z-axis are equivalent. So, we have also used the term C2 for inversion symmetry. The space with character +1 under C2 and +1 under e-h symmetry corresponds to the A+ space and contains covalent VB diagrams, i.e. VB diagram in which every site is neutral. The ground state of the system is found in this space. The optically connected excited state (to the ground state) is found in the B− space, where the character under both C2 and e-h symmetry is -1. This space does not contain any covalent VB diagram and hence is also called the ionic space. The e-h symmetry exists only in bipartite lattices where the transfer term is nonzero only between atoms in different sublattices. Among the PAHs we have studied only Fluoranthene does not belong to this class. We use the symmetry adapted VB basis to set-up the Hamiltonian matrix. 38 A few low-lying states in each subspace are obtained from the Rettrup algorithm. 39 With the computational resources we have, we can do exact calculations for the singlets when both C2 and eh symmetries exist, for a system with 20 carbon atoms and 20 electrons (neutral system). The dimensionalities of the subspaces for all the PAH systems we have studied is given in Table 1. Wherever we are unable to carry out exact calculations, we have resorted to the DMRG calculations. Thus triplets of benzopyrene and both singlets and triplets of picene were obtained within the DMRG approach. After obtaining the low-lying eigenstates, to fully characterize the state, we have computed the transition dipole between ground and excited states in appropriate spaces. We have also computed bond orders to understand the nature of equilibrium geometry in all eigenstates. In triplet states we have computed spin densities which can be verified from ESR studies.

ACS Paragon Plus Environment

8

Page 8 of 73

Page 9 of 73 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: Dimensionalities of the symmetry adapted spaces for PAHs studied in this paper. The total number of bases in singlet and triplet spaces for neutral systems with N = 16 are 34,763,300 and 66,745,536; with N = 18 are 449,141,836 and 901,995,588; with N = 20 are 5,924,217,936 and 12,342,120,700; with N = 22 are 79,483,257,303 and 170,724,392,916, respectively. Molecule Pyrene Benzanthracene Chrysene Helicene Triphenylene BenzoPyrene

Ne = N 16 18 18 18 18 20

Singlets

Triplets

Acov

Bionic

Bcov

8,698,485 224,573,349 112,313,738 224,573,349 112,313,738 1,481,162,738

8,695,320 224,568,487 112,303,369 224,568,487 112,303,369 1,481,122,588

16,692,888 451,003,761 225,524,745 451,003,761 225,524,745 -

Aionic 16,689,604 450,991,827 225,510,840 450,991,827 225,510,840 -

Results and Discussion All the molecules studied, with the exception of fluoranthene, have electron-hole symmetry. Thus the states in these molecules (except fluoranthene) can be classified as ’covalent’ (’+’) or ’ionic’ (’-’). The ground state is in the covalent subspace, so are the 2-photon states. Because of the non-alternancy, only fluoranthene will have non-zero dipole moment (4.22 Debye) in the ground state. The excitations to low-lying states in all molecules are shown in Fig. 2. In cases where the energy gap between levels in a molecule are small, we have ensured the linear independency of eigen states by computing the overlap of nearly degenerate eigen states. One common feature we note is that in all the molecules, the two-photon excitation is below the onephoton excitation, except in the case of fluoranthene, where the lowest excited singlet is in the B space and thus has a nonzero transition dipole. The lowest two-photon state in fluoranthene is about 0.2eV above the lowest one-photon state. In all other molecules, the lowest two-photon state is about 0.7 eV below the lowest one-photon state. Hence except fluoranthene, none of the other molecules will be fluorescent by Kasha rule. 40 It is also reasonable to assume that the (2A−1B) gap is a measure of the strength of electron correlation. We find it lowest in fluoranthene, implying that the effective electron correlations are ACS Paragon Plus Environment

9

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

3A

(3.4, 0.0)

2B 2A 1B

(0.0, 1.0)

(0.0, 5.3)

4.24 3.78 3.66 3.45

-

2.87

3B2 + 4A 2B3 + 3A1B2

2.34

2A

3

2B 3 3A 3 2A 1

(0.0, 0.3)

3.96 3.76 3.57

3

A (0.0, 0.4)

1 3B

(0.2, 0.0)

(3.2, 0.0)

(5.1, 0.0)

Page 10 of 73

(0.0, 2.8) 5.41 4.94 4.51 3.81 3.75 3.07

+

43 B2 33 B3 23 B2 3 1 B3 5A 4A 3 A2 A-

(0.0, 3.5)

5.68 5.59 5.30 4.77 4.10 3.88 3.28 3.12

(2.4, 0.0) (0.0, 5.3)

+

1B

1.86 (2.3, 0.0)

1A

1A

+

+

1A

1A

(4.2, 0.0)

Fluoranthene

3S 2S + 3S 1S +

2S

(5.9, 2.8) (5.6, 0.5)

4.85 4.53

3.84 (1.4, 3.1) 3.77 3.07

3 T2T 1T 3T 2T 1T

1S

+

1S

PYRENE

5.12 4.97

3B 2B

4.09

3A 1B

-

+

+

3.26 +

(0.8, 4.2) 5.23 5.02 (8.7, 1.9)

+

2A

2.78

+

4.21 3.93 (2.5, 3.3) 3.22

3 3A 3 2A 3 1A

3

+

3

+

3B 3 + 2B 1B

(0.2, 0.2)

5.98

5.21 4.99 (1.4, 1.0) 3.64 3.33 (0.3, 0.1) 2.34

1.96

+

+

1A

BENZANTHRACENE

+

1A CHRYSENE

Figure 2: Low-lying excitation gaps for PAH molecules. Molecules are labelled in each box. The states shown to the left are singlets and those to the right are triplets. In the case of Fluoranthene, the symmetry labels were assigned from analyzing the eigenstates and confirmed from the polarization of the transition dipole. The transition dipoles are in Debye and energies are in eV. The spacing between labels is scaled to the magnitude of the gaps. The X and Y components of the transition dipoles are given in parenthesis and the thickness of the arrow indicates the strength of the transition. State labels with ’+’ and ’-’ superscripts correspond to ’covalent’ and ’ionic’ subspaces.

ACS Paragon Plus Environment

10

Page 11 of 73 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

3B 2B 1B+ 3A

(3.7, 0.1)

4.99 (6.9, 0.0) 4.64 4.04 3.90

+

2A

3.16

3A 2A 1A +

3B + 2B + 1B

(4.0, 0.0)

3 3A

5.46 (0.0, 0.4)

3B 2B 1B+ 3A

4.98 4.14

(0.3, 0.0) 3.32 2.90 2.46

(0.0, 4.3) (0.0, 6.4)

4.22 4.10

+

2A

6.22 (2.4, 0.0)

5.37 5.11

3.38

2A 3 1A 3 + 3B 3 + 2B 3 + 1B 3

5.16 (1.1, 0.0)

4.36 3.57 3.45 2.83

(0.0, 0.8)

+

+

1A

+

1A

+

1A

1A

HELICENE

TRIPHENYLENE

3 3A

3B 2B 1B

(0.0, 2.6) (1.5, 0.0)

5.52 4.60 3.97 3.91 3.23

+

2A

3 2A 3 1A 3

(0.0, 0.7)

3B 2B 1B

5.84 5.19 4.51

+

2A 1A

(1.2, 0.2) (0.8, 0.0)

4.70

3

4.27 3.67

+

3.17

3

+

3

+

3B

2B 3

(1.4, 0.1)

4.74 4.19 3.71

3.19 (0.3, 0.3)

+

1B

+

(0.0, 4.3)

5.50 (0.1, 0.1)

3

2A

(0.0, 0.1)

3

+

3

7.12

2B

1B

1A

3A

(0.0, 2.5)

2.24

2.09

(1.3, 0.1)

+

1A

+

+

1A BENZOPYRENE

1A PICENE

Figure 2 (contd.): Low-lying excitation gaps for the PAH molecules.

ACS Paragon Plus Environment

11

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

weakest in this system. Of the alternant systems, the weakest effective correlation is in picene and all the other systems show comparable effective correlation. The strongest optical absorption to 1Bu state is in benzopyrene and the weakest in chrysene. In many of these systems, we do observe strong excitations to higher singlet states. In Table 2, we give a comparison of experimental gaps with the theoretical gaps. 41,42 We find that our theoretical gaps compare very well with the experimental gaps. For pyrene molecule, we also find good agreement with the theoretical result of Basak et al obtained using the standard PPP parameters, taking quadruple CI basis. 30 Even the lowest two-photon gap compares well with their results from the standard parameters. 31 Triphenylene molecule having C3 symmetry shows two two-photon levels below B− state. Here the excitation to the first B− state is very weak and the second level is ≈ 0.5eV higher than the experimental value. 43,44 In other cases also the experimental gaps are nearly 0.5 eV lower than the computed values, which can be attributed to solvent effects. In case of benzopyrene, the difference is slightly larger because the gap is obtained from the truncated CI basis within DMRG method. 45 Furthermore, since the molecule is not a quasi-one dimensional system (in contrast to picene), the excitation gap is slightly higher than the experimental gap. We have given one-photon, two-photon and spin gaps of these molecules in Fig. 3. We find that in all the cases except pyrene and picene, the lowest triplet state is well above half of the singlet-singlet gap and hence none of the molecules are suitable for improving the solar cell efficiency via singlet fission. However, for pyrene and picene, since [ES1 − 2 × ET 1 ] value is almost zero, these molecules can be tuned to have a favourable conditions for singlet fission, by substitution. In fact, B-N substituted pyrenes are shown to be possible candidate materials for singlet fission by Zeng et al. 46 In fluoranthene, pyrene, benzanthracene, benzopyrene and picene, the (2A − T1 ) gap is ≥ 1.0 eV. However, triphenylene has the lowest 2A − T1 gap of 0.5 eV, while the gap is 0.7 eV in helicene, and is 0.9 and 1.0 eV in chrysene and benzopyrene, respectively. We notice that the optical gap falls off rapidly from anthracene to pentacene (3.7 eV to 2.92 eV, see Table S1 in supporting information), whereas for phenacenes, it is slower ACS Paragon Plus Environment

12

Page 12 of 73

Page 13 of 73 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

(4.3 eV for phenanthrene to 3.82 eV for picene). In the case of oligoacenes, the spin gap almost reaches zero, whereas for phenacenes, it is quite large, even in case of picene (1.95 eV). Table 2: Electronic excitation gap and singlet-triplet gaps in eV for PAH molecules. Numbers in the bracket are experimental gaps obtained from references shown as superscript. Except benzopyrene and picene other gaps are numerically exact. For benzopyrene and picene, the singlet-singlet and singlet-triplet gaps are obtained from DMRG studies. Molecule Fluoranthene Pyrene Benzanthracene Chrysene Helicene Triphenylene Benzopyrene Picene

Optical gap singlet(Expt gap) triplet gap 3.45 (3.46) 41 2.34 41,42 3.75 (3.34) 1.86 41 3.77 (3.23) 1.96 3.95 (3.45) 20 2.34 19 4.04 (3.49) 2.47 43,44 5.11 (4.36) 2.83 3.91 (3.38) 45 2.28 41 3.82 (3.30) 1.95

Singlet State Properties Bond orders We have computed the bond orders of the PAH molecules in various electronic states. The bond  1 order between sites ’i’ and ’j’ is given by the expectation value of pˆi j = − Eˆi j + Eˆ ji in the 2 † state of interest, where Eˆi j = ∑σ ciσ c jσ . The larger the bond order, the shorter is the bond length and vice versa. Comparison of bond orders in a given excited state gives the propensity of a bond to shorten or lengthen with respect to the ground state bonds. At equilibrium, the ratio of bond length to bond order is expected to be the same for all bonds.

In Fig. 4, we have given the bond orders in the singlet ground state as well as the lowest onephoton excited state for all the PAHs we have studied. For comparison, we have given the bond ACS Paragon Plus Environment

13

The Journal of Physical Chemistry

5

One-photon state Two-photon state Lowest triplet state

4

Energy gaps in eV

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

1

0

e e ra en yren h h t t n e P zacen ran n o u Be Fl

ne

se ry h C

e e e e en en ren l en c y i c y l i p n P o He he nz ip e r B T

Figure 3: Energy gaps in eV, for one-photon state (black filled bar), two-photon state (blue hatched bar) and triplet states (red dotted bar) for PAHs molecules. orders and spin densities of naphthalene, anthracene, phenanthrene and biphenyl molecules in supporting information (Fig. S1 and S2).

Molecules with 16 π-electrons (Fluoranthene and Pyrene) In the ground state of fluoranthene, we note that the bond order of the bonds connecting the naphthalene unit and benzene unit are very small. This indicates that the two units are only weakly connected by the π−conjugation. The bond orders in the six membered rings are close to those of benzene in the H¨uckel model ∼0.66, although bonds which are shared between rings have slightly lower bond order. In the optically excited state the bond connecting the benzene and naphthalene units become much stronger, indicating a contraction of the distance between the two units. The other bonds also show slightly smaller but nearly uniform bond ACS Paragon Plus Environment

14

Page 14 of 73

Page 15 of 73

0.6 5

5

0.6 0.58 5 0.4 0.5

0.57

0.55 0.4 9 9 0.4

3

3 0.5 0.48

0.48

0.6

4

0

0.7

2

6

9

0.6

0.49

3 0.6

6

5 0

0.4

0.5 0.47

7

0.60

8 0 0.4 .52

7

7

0.6

[5]

6 0.5

9 2

8 0.6

9 0.5

5

0.5

4

8

0.5 9 0.4

3 0.5

0.5

1 0.5

8

0.6

0.3 0.5

2 0.6 7

7

0.5

9

0.5

7

0.4

0.72 0.60

0.41 0.33

2

9

9

6

0.4

5

0.49

0.40

0.5

0.5

0.6

0.54

0.5

1

5

0.42

0.65

0.6

0.5

0.52

0.65 0.48

2

3

0.69

0.6

0.55

0.25 0.46

0

0.4

0.51 0.43

0.5

7

0.5

0.65

6

0.6

5

7

8

0.63

0.58 0.64

7

0.5 5

0.6 7

.48

0.39

0.7

5

0.55 0.6

0.5

1

0.5

5

0.5

0.47

1

0.6

0.4

0.61

0.5

0.3 7

1

0.6

6 0.4

7

0.4 4

1

0.55

0.55

0.5

0.62

0.6

3

0.6

6

5 0.6

8 0.5 0.6

3

0.6

0.61 0.57

0.4

0

[4]

[3]

[6]

65

0.4 0 .44 0 1 0.42 0. 0.7 0. 55 0 60 .60 0.58 0.55 0.49 5 4 7 5 . . 0 0.7 0 0.4 0 7

0 7 .50 .46 0.6 0 1 .57 0 0.57 0 0.5 0.62 0.46 0.47 0.53 2 54 .64 0.5 0 0.6 0. 0.5 .53 0 2 4

0.43

8 0.7

0.7

0.70

3 0.5 0.

0.58

0.71

0.4

0.67

0

0.5

[3]

1 0.6 0. 51

0.5

0.37 0.46

0.5

0.68

6

0.5 6

.62

6

8

0.58

.63

.50

0.5

0.57 0.43

0.51

0.51

0

[2]

[1]

0.55

0.4

0

0.55

0.6

4

6 0.5 0

2 0.7 0. 5

0.5

3

0.5

5

0.6

0.64

0.5

0

0.5

0

0.5

2

2

0.5

0

0.64

8 0.6 0 .60 0.57

0.77 0.64

0.49

0.7

1

8

0.7

2 0.6 5 0.6 0.59 0.63 0.45 0.57 0.5 5 2 5 0 0.4 0. 0.28

0.58 3 0.45 0.5

3

0.77

0.6

0.5

0.54 0.52

0.6

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

0.52

0.63

0.75

[8]

[7]

Figure 4: Bond orders in singlet ground state (in red) and in dipole allowed state (in blue). The numbering for molecules is as given in Figure 1.

orders, indicating ring expansion in the excited state. This also shows that the Stokes shift in the spectral lines will be large in this system. This is in confirmity with large Stokes shift

ACS Paragon Plus Environment

15

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

observed experimentally, by Gusten and Heinrich. 47 In pyrene, the bond orders in the ground state resemble those of biphenyl connected by two ethylenic units on either side, with alternate single and double bonds, rather than the bond orders of naphthalene connected by two propylenic units at the top and bottom. This is in accordance with the most possible resonance structures for pyrene. Compared to biphenyl, the central bond becomes stronger in pyrene (see Fig. S1). In the optically excited state, all the bond orders become more uniform. This shows that at equilibrium geometry, the molecule is not distorted much from the uniform bond lengths assumed in the calculations. We find a uniform expansion of rings in the optically excited state. We should expect the Stokes shift due to relaxation of the geometry in the excited state to be small as opposed to what is expected in fluoranthene.

Molecules with 18 π-electrons (Benzanthracene, Chrysene, Helicine and Triphenylene) Benzanthracene molecule can be viewed as a phenyl ring attached either to an anthracene or to a phenanthrene. This is well reproduced by the bond order pattern in the ground state, where one of the two terminal rings have anthracenic like and the other two have phenanthrenic like bond orders (see SI Figure 1). In benzanthracene, the middle bond of armchair of the phenanthrenic moiety is the weakest while the bond opposite to that is the strongest. The ground state bond orders of both chrysene and helicene, again, resemble those of phenanthrene with a marginal change in bond orders. The middle bond of armchair of the phenanthrenic unit is the weakest and the bond opposite to that is the strongest, in accordance with the possible resonance structures. This is true even for helicene, where there are two weaker bonds at the arm chair, connecting the two naphthalenic unit. In both molecules, there is a slight enlargement of the rings in the excited state. The Stokes shift in these system are not expected to be large. The ground state bond orders in triphenylene molecule, correspond to three phenyl rings connected by a very weak bond (bond order is 0.37). The bond order pattern is reversed upon excitation, with the central ring having a more uniform and slightly stronger bonds compared to the ground ACS Paragon Plus Environment

16

Page 16 of 73

Page 17 of 73 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

state geometry. Even in this molecule, the excited state geometry has more uniform bond orders. Once again the change in geometry from the ground state to the excited state is small, which leads us to expect small Stokes shift. Notice that the effect of non-conjugation in triphenylene is reflected in the excitation gap. Among all the PAHs we have studied, triphenylene has the largest optical gap and the singlet-triplet gap (see Table 2). Among 18 π− electronic molecules, benzanthracene has the lowest gap while chrysene and helicene have almost similar gap (3.95 eV and 4.04 eV, see Table 2) owing to their similarity in conjugation pattern.

Molecules with 20 (Benzopyrene) and 22 π-electrons (Picene) Benzopyrene molecule has similar bond orders as those of triphenylene with an extra ethylenic unit connected at lower end of the molecule (see [7] in Fig. 4). Picene being a oligomer of phenanthrene, shows the ground state bond order similar to phenanthrene with minimal change in the bond orders of interior bonds. The change in equilibrium geometry from the ground to the excited state is not much, although the bond lengths in the excited states are slightly larger. Thus, even in picene, we should not expect large Stokes shift in the fluorescence.

Triplet State Properties Bond-Orders Fig. 5 gives bond orders of the PAH molecules in the lowest triplet and the lowest two photon states. In fluoranthene (1) the bond common to the two benzene rings becomes very weak and elongated while the peripheral bonds, except the one that connects the phenyl ring to the five membered ring, all have similar bond lengths at equilibrium. However the triplet bond orders show more alternation and only the bond connecting the benzene and naphthalene units stays ACS Paragon Plus Environment

17

The Journal of Physical Chemistry

0.6

4

0.6

0.32 0.61

0.5 0

0.5

6

0.5

6

[2]

[1]

1 0.7 0

3 0.6 0 .59 0.53 0.61 6 4 0 6 6 5 5 5 . . . . 0 0 0. 0 0. 0 0.5 52 45 3 0.54 0.38 0.36 0.49 3 5 4 5 6 5 . . . 0 0. 0 0 0.5 0.6 55 3 5

0.45

0 5 .57 .44 2 0.4 0. 0 0 0.6 0.5 55 .53 4 0.56 0.45 0.43 2 3 0.7 0.6 0.5 0. 47 3

7 0.6

0.6

7 0.5

0.5

4

4 0.6

2

4 0.6

0 2

0.4

0.4 8

1 0.5

3

1

6

0.6

6 0.4

2

5 0.6

3

0.5

3 0.5 0.5

6

0.29 0.44

4

0

0.40 0.60

0.47

0.5

1

0.5

0.53

0.66

0.4

7

0.53

0.37 0.52

8 0.5

0.39

9 0.4

0.6

0.4

1

3 0 0.5 .53

4 0.5

0.57

8

0.39 0.31

0.5

1 0.6

1

9 0.5

5

4

0.5 0.52

0.43 0.39

0.65 0.5

0.4

0.47 0.41

9

0.46

0.6

3

0.5

0.62 0.58

0.4

2

3 0.5

0.5

1

5 0.48

[5]

0

1

0.5

6 0.50 0.4 2

0.6

0.4

7

0.5

7

0.6

3

0.53

0.6

0.4

0.6

5

9

0.4

0.48

0.63

[6]

0.64

0.4

5

0.5

3

3

0.61

0.54

0.56

0.54

0 0.6

6

0.5

0.6

0.5 8 0.5

0.45 0.49

9

0.59

0.5

4 0.6

0.5

4

0.6

[4]

0.58 0.56

0.53

3 0.5 0. 59 0.43

.48

[3]

0.47 0.43

[3]

3

0.53

0 0.5 0.54

0.54

0

0.63 6 9 0.5 0. 0.4 0. 6 47 54 1 0.5 0.54 0.40 0.38 0.48 7 3 6 6 5 5 . . . 0 0 0.5 0 0.6 0.5 8 3 4

9 0.6 0.5

0.5

4

7

0.5

0.63 0.60 7 .6

0.56 0.68

0.52 0.63 0. 42 2 0.4

0.68

4 0.6 0. 60 0.54

0.6

2 0.5 8

0.5

0.35 0.59

.51

0.56

0.6 4

0.56 0.55

0.4

5

0.43 0.5 1 9 0.5

0.6

9

4

0.4

0.4

3

0.49 0.28

0.5

4

0.43 9 0.5 0

0.72

0.5

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 18 of 73

0.55

0.67

0.57

[8]

[7]

Figure 5: Bond orders in lowest triplet state (in magenta) and in two photon allowed state (in green). The numbering for molecules is as given in Figure 1.

elongated. In pyrene (2), bond orders in the triplet and two-photon allowed state show an opposite pattern. i.e. the bond orders are all nearly uniform in the naphthalenic unit in the triplet state and in biphenyl unit, they are alternating. On the other hand, in the two photon state, bonds are uniform in the biphenylenic unit whereas alternating in the naphthalenic unit. In benzanACS Paragon Plus Environment

18

Page 19 of 73 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

thracene (3), both the triplet and the two photon states have similar bond order patterns, with nearly uniform peripheral bonds and weaker inter-ring bonds. In chrysene as well as in helicene, both the triplet and the two photon states show bond alternation, although it is more pronounced in the two photon state. In helicene the bond lying at the C2 axis is the weakest bond both in the triplet state and two-photon allowed state (2A). In triphenylene (6), bond orders in both 2A and triplet states are similar and nearly uniform along the benzenic unit with a weaker central bond connecting the three phenyl rings. In benzopyrene (7), the interior bonds through which the C2 axis passes through, are weaker. In picene (8), the outer bonds have a marginal bond alternation in both the states, except the middle bond of the arm-chair of phenanthrenic units (0.35 ± 0.05). In the triplet state, the bonds connecting the two terminal naphthalene units to the central benzene ring are very weak (0.37 and 0.29), but only slightly weaker in the case of the 2A state (0.52 and 0.44) compared to other bonds.

Spin Densities Spin density at a site ”i” is defined as the difference in charge density of up and down spins. i.e. ρi =< ψk |a†iα aiα − a†iβ aiβ |ψk >= 2 < ψk |Sz |ψk >. In Fig. 6 are given the spin densities in the lowest energy triplet state of the PAH molecules. We find large positive spin densities at the α sites of the naphthalenic units of fluoranthene. There is slightly negative spin density at the sites common to the two benzene rings in the naphthalene unit. The total spin density is almost localised on the napthalenic unit and spin density pattern resembles that of naphthalene (see supporting information Fig. S2). On the other hand, while isolated benzene has a significant spin density at each site (0.17), benzenic unit of fluoranthene has negligible spin density (< 0.02 ± 0.01). In pyrene, the largest positive spin densities are found at the α sites (0.18) of biphenylinic unit followed by the β carbon atoms of naphthalenic unit (0.09). There is negative spin density at the two peripheral carbon atoms of biphenylinic unit (−0.05) and all other carbon ACS Paragon Plus Environment

19

The Journal of Physical Chemistry 1 2 3 4 5 6 −0.05 7 0.18 0.18 8 9 0.03 0.01 0.01 0.00 0.07 0.04 10 0.05 0.06 0.09 0.09 11 0.03 0.17 0.11 0.03 0.04 0.20 12 0.05 0.02 0.09 0.09 0.01 13 0.00 0.04 −0.03 −0.02 0.08 0.06 0.03 0.01 0.01 14 0.01 0.02 0.05 0.06 0.06 15 0.20 0.18 0.18 0.09 0.17 0.03 16 0.03 −0.05 0.03 0.10 0.07 0.04 17 18 [1] [2] [3] 19 20 21 22 0.03 0.03 23 0.05 0.08 0.09 0.09 0.04 0.04 24 0.06 0.04 0.02 25 −0.00 −0.00 0.10 0.09 0.09 0.19 −0.02 0.13 26 0.00 0.01 0.01 0.00 27 −0.02 0.19 0.10 28 0.04 0.02 0.12 0.12 0.06 29 0.04 0.05 0.05 30 0.05 0.08 31 [5] [4] 32 33 34 35 0.04 0.04 36 0.05 0.05 37 0.09 0.09 38 0.07 0.07 0.03 0.03 39 0.15 −0.00 0.15 −0.00 0.05 0.05 0.05 0.05 40 0.07 0.07 0.08 0.08 41 −0.02 −0.00 0.07 0.07 0.04 −0.00 −0.02 0.05 0.05 0.04 0.00 0.00 42 0.02 0.02 43 0.05 0.07 0.07 0.06 0.06 −0.02 0.05 0.05 0.05 −0.02 0.00 0.00 0.14 44 0.14 0.16 0.02 0.02 0.16 45 0.05 0.07 0.07 0.05 0.02 0.02 46 47 [6] [7] [8] 48 49 50 Figure 6: Spin densities in lowest triplet state are given in red color. The numbering of 51 molecules are as given in Figure 1. 52 53 54 atoms have negligible spin densities. 55 56 Unlike the bond order pattern of benzanthracene, spin density pattern neither resembles that 57 58 59 ACS Paragon Plus Environment 60

20

Page 20 of 73

Page 21 of 73 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

of anthracene nor phenanthrenic unit. A uniform spin density (0.05 ± 0.01) is seen except at the α− carbon atoms of anthracenic moiety (0.11). No negative spin density is observed in this molecule. Chrysene shows an opposite spin density pattern compared to other molecules. Here, the spin density is largest (0.19) at the two carbon atoms where two naphthelinic units are fused. It also shows negative spin densities at one of the adjacent sites followed by a large positive spin density (0.10 and 0.08) at the α sites of naphthalenic units. On the other hand, in helicene almost all the α sites have large positive spin densities (∼ 0.1) and β sites have smaller spin densities. Triphenylene has a uniform spin density of 0.05 and 0.07 at alternate sites, with the central benzenic unit having a uniform spin density of 0.05, owing to its three-fold symmetry. Benzopyrene has similar spin density pattern as that of triphenylene, with a pronounced alternation in spin densities at adjacent sites. The highest positive spin density of 0.14 is seen at the carbon atoms adjacent to the vinylenic bridge in benzopyrene. Picene molecule, on the other hand, has alternate positive and negative spin densities at the peripheral phenyl rings. There is also large positive spin densities (0.15) at four alpha carbon atoms of the interior naphthalenic units. The carbon atoms of the interior arm-chair have spin density of 0.07.

Summary In this paper we have studied the correlated electronic states of eight well known polycyclic aromatic hydrocarbons. From energetics, we find that in all these molecules, the one-photon state is above the two photon state. Hence none of these molecules are expected to be fluorescent. The optical gap obtained from our study agree well with the experimental gaps. The singlet-triplet gap in all the molecules is also quite large and does not satisfy the requirements for singlet fission. The bond orders show that the equilibrium geometries of all the states are nearly the same as assumed in the calculations. Thus, we expect small Stokes shifts in fluorescence. The ACS Paragon Plus Environment

21

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

spin densities in the triplet state are mostly confined to the α sites of the naphthalenic units in the molecules. Most molecules show very small negative spin densities at a few sites and often all the sites have only positive spin densities.

Supporting Information Available We have given the following table and figures: i) A Table for electronic excitation gap and singlet-triplet gaps for PAH molecules - naphthalene, anthracene and phenanthrene. ii) A figure for bond orders for the lowest singlet and triplet states for naphthalene, anthracene, phenanthrene and biphenyl. iii) A figure containing spin densities of the lowest triplet states for naphthalene, anthracene, phenanthrene and biphenyl. This material is available free of charge via the Internet at http: //pubs.acs.org/.

Acknowledgement Y. A. P. thanks Disha programme for women in Science (DST0128), Department of Science and Technology, Government of India, for financial support.

References (1) Salem, L.; Molecular Orbital Theory of Conjugated Systems, W. A. Benjamin: New York, (1966); Chapter 4. (2) Kuwajima, S.; Soos, Z. G. Enhanced Ring Current of Charged Annulenes, J. Am. Chem. Soc., 1986 108, 1707-1708.

ACS Paragon Plus Environment

22

Page 22 of 73

Page 23 of 73 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

(3) Gomes, J. A. N. F.; Mallion, R. B. Aromaticity and Ring Currents, Chem. Rev., 2001, 101, 1349-1384. (4) Zhang, Y.; Sheng, W.; Li, Y. Dark Excitons and Tunable Optical Gap in Graphene Nanodots, Phys. Chem. Chem. Phys., 2017, 19, 23131-23137. (5) Hai, X.; Feng, J.; Chen, X.; Wang, J. Tuning the Optical Properties of Graphene Quantum Dots for Biosensing and Bioimaging, J. Mater. Chem. B, 2018, 6, 3219-3234. (6) Chen, S.; Ullah, N.; Wang, T.; Zhang, R. Tuning the Optical Properties of Graphene Quantum Dots by Selective Oxidation: a Theoretical Perspective, J. Mater. Chem. C, 2018, 6, 6875-6883. (7) Wettstein, C.M.; Bonafe, F. P.; Oviedo, M. B.; Sanchez, C. G. Optical Properties of Graphene Nanoflakes: Shape Matters, J. Chem. Phys., 2016, 144, 224305. (8) Sharma, SRKC Y.; Bandyopadhyay, A.; Pati, S. K. Structural Stability, Electronic, Magnetic, and Optical Properties of Rectangular Graphene and Boron Nitride Quantum Dots: Effects of Size, Substitution, and Electric Field, J. Phys. Chem. C, 2013, 117, 23295-23304. (9) Alam, Sk. M.; Ananthanarayanan, A.; Huang, L.; Lim, K. H.; Chen, P. Revealing the Tunable Photoluminescence Properties of Graphene Quantum Dots, J. Mater. Chem. C, 2014, 2, 6954-6960. (10) Geethalakshmi, K. R.; Yong Ng Y.; Crespo-Otero, R. Tunable Optical Properties of OHFunctionalised Graphene Quantum Dots J. Mater. Chem. C, 2016, 4, 8429-8438. (11) Mitsuhashi, R.; Suzuki, Y.; Yamanari, Y.; Mitamura, H.; Kambe, T.; Ikeda, N.; Okamoto, H.; Fujiwara, A.; Yamaji, M.; Kawasaki, N.; Maniwa, Y.; Kubozono, Y. Superconductivity in Alkali Metal Doped Picene, Nature, 2010, 464, 76-79.

ACS Paragon Plus Environment

23

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

(12) Xue, M.; Cao, T.; Wang, D.; Wu, Y.; Yang, H.; Dong, X.; He, J.; Li, F.; Chen, G. F. Superconductivity Above 30 K in Alkali-Metal-Doped Hydrocarbon, Sci. Rep., 2012, 2, 389. (13) Artioli, G. A.; Hammerath, F.; Mozzati, M. C.; Carretta, P.; Corana, F.; Mannucci, B.; Margadonna, S.; Malavasi, L. Superconductivity in Sm-Doped [n] Phenacenes (n =3,4,5), Chem.Commun., 2015, 51, 1092-1095. (14) Nakagawa, T.; Yuan, Z.; Zhang, J.; Yusenko, K. V.; Drathen, C.; Liu, Q.; Margadonna, S.; Jin, C. Structure and Magnetic Property of Potassium Intercalated Pentacene: Observation of Superconducting Phase in KxC22 H14 , J. Phys.: Condens. Matter, 2016, 28, 484001. (15) Allamandola, L. J.; Sandford, S. A.; Wopenka, B. Interstellar Polycyclic Aromatic Hydrocarbons and Carbon in Interplanetary Dust Particles and Meteorites, Science, 1987, 237, no. 4810, 56-59. (16) Allamandola, L. J.; Tielens, A. G. G. M.; Barker, J. R., Polycyclic Aromatic Hydrocarbons and the Unidentified Infrared Emission Bands - Auto Exhaust Along the Milky Way, Astrophys. J. Part II, 1985, 290, L25-L28. (17) d’Hendecourt L.; Ehrenfreund P. Spectroscopic Properties of Polycyclic Aromatic Hydrocarbons (PAHs) and Astrophysical Implications, Adv. Space Res., 1997, 19, 1023-1032. (18) Salama, F.; Galazutdinov, G. A.; Krelowski, J.; Biennier, L.; Beletsky, Y.; Song, I-O. Polycyclic Aromatic Hydrocarbons and the Diffuse Inetestellar Bands: A Survey, Astrophys. J., 2011, 728, 154. (19) Palewska, K.; Ruziewicz, Z.; Chojnacki H.; Meister, E. C. High-resolution Electronic Spectra of Tetrahelicene and Hexahelicene in Low-Temperature Polycrystalline Matrices, Chem. Phys., 1992, 161, 437-445.

ACS Paragon Plus Environment

24

Page 24 of 73

Page 25 of 73 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

(20) Hussain, M. M. Measurement and Theoretical Characterization of Electronic Absorption Spectrum of Neutral Chrysene and Its Positive Ion in H3 BO3 Matrix, Spectrochim. Acta A, 2007, 68, 156-164. (21) Canuto, S.; Zerner, M. C.; Diercksen, G. H. F. Theoretical Studies of the Absorption Spectra of Polycyclic Aromatic Hydrocarbons, Astrophys. J. Part I, 1991, 377, 150-157. (22) Ping Du, Salama, F.; Loew, G. H. Theoretical Study of the Electronic Spectra of a Polycyclic Aromatic Hydrocarbon, Naphthalene, and its Derivatives, Chem. Phys., 1993, 173, 421-427. (23) Langhoff, S. R. Theoretical Infrared Spectra for Polycyclic Aromatic Hydrocarbon Neutrals, Cations, and Anions, J. Phys. Chem., 1996, 100, 2819-2841. (24) Malloci, G.; Mulas, G.; Joblin, C. Electronic Absorption Spectra of PAHs up to Vacuum UV Towards a Detailed Model of Interstellar PAH Photophysics, Astron. Astrophys., 2004, 426, 105117. (25) Filipp, F.; Reinhart, A. Adiabatic Time-Dependent Density Functional Methods for Excited State Properties, J. Chem. Phys., 2002, 117, 7433-7447. (26) Hammonds, M.; Pathak, A.; Sarre, P. J. TD-DFT Calculations of Electronic Spectra of Hydrogenated Protonated Polycyclic Aromatic Hydrocarbon (PAH) Molecules: Implications for the Origin of the Diffuse Interstellar Bands, Phys. Chem. Chem. Phys., 2009, 11, 44584464. (27) Cocchi, C.; Prezzi, D.; Ruini, A.; Caldas, M. J.; Molinari, E. Anisotropy and Size Effects on the Optical Spectra of Polycyclic Aromatic Hydrocarbons, J.Phys.Chem.A, 2014, 118, 6507-6513.

ACS Paragon Plus Environment

25

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

(28) SeyedAbdolreza, S.; Yong Z.; Sun, K. A Theoretical Study on the Vibrational Spectra of Polycyclic Aromatic Hydrocarbon Molecules with Aliphatic Side Groups, Astrophys. J., 2015, 801, 34. (29) Aryanpour, K.; Roberts, A.; Sandhu, A.; Rathore, R.; Shukla A.; Mazumdar, S. Subgap Two-Photon States in Polycyclic Aromatic Hydrocarbons: Evidence for Strong Electron Correlations, J. Phys. Chem. C, 2014, 118, 3331-3339. (30) Basak, T.; Chakraborty, H.; Shukla, A. Theory of Linear Optical Absorption in DiamondShaped Graphene Quantum Dots, Phys. Rev. B., 2015, 92, 205404. (31) Basak, T.; Basak T.; Shukla, A. Electron Correlation Effects and Two-Photon Absorption in Diamond-Shaped Graphene Quantum Dots, Phys. Rev. B., 2018, 98, 035401. (32) Lykos, P. G. and Parr, R. G. On the π Electron Approximation and Its Possible Refinement, J. Chem. Phys., 1956, 24, 1166-1173. (33) Pariser, R.; Parr, R. G. A Semi-Empirical Theory of the Electronic Spectra and Electronic Structure of Complex Unsaturated Molecules. II, J. Chem. Phys., 1953, 21, 767-776. (34) Pople, J. A., The Electronic Spectra of Aromatic Molecules II: A Theoretical Treatment of Excited States of Alternant Hydrocarbon Molecules Based on Self-Consistent Molecular Orbitals, Proc. Phys. Soc. A, 1955, 68, 81-89. (35) Ohno, K. Some Remarks on the Pariser-Parr-Pople Method, Theor. Chim. Acta, 1964, 2, 219-227. (36) Soos, Z. G.; Ramasesha, S. Valence Bond Approach to Exact Nonlinear Optical Properties of Conjugated Systems, J. Chem. Phys., 1989, 90, 1067-1076. (37) Ramasesha, S.; Soos, Z. G. Diagrammatic Valence-Bond Theory for Finite Model Hamiltonians, Int. J. Quantum Chem., 1984, 25, 1003-1021. ACS Paragon Plus Environment

26

Page 26 of 73

Page 27 of 73 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

(38) Ramasesha, S.; Soos, Z. G. Symmetry Adaptation of Correlated States in the Valence Bond Basis, J. Chem. Phys., 1993, 98, 4015-4022. (39) Rettrup, S. An iterative Method for Calculating Several of the Extreme Eigensolutions of Large Real Non-Symmetric Matrices, J. Comp. Phys., 1982, 45, 100-107. (40) Kasha, M.; Characterization of Electronic Transitions in Complex Molecules, Discuss. Faraday Soc., 1950, 9, 14-19. (41) Clarence, K. Jr. The Longest Wavelength Band in the Electronic Spectra of Polycyclic Aromatic Hydrocarbons for Analytical Use, Appl. Spectr., 1959, 13, 15-25. (42) Sun, Y. P. Excitation Wavelength Dependence of Pyrene Fluorescence in Supercritical Carbon Dioxide. Evidence for a Supercritical Solvent-Assisted Solute-Solute Clustering Mechanism J. Am. Chem. Soc., 1993, 115, 3340-3341. (43) Hall, R. D.; Valeur, B.; Weber, G. Polarization of the Fluorescence of Triphenylene: a Planar Molecule with Threefold Symmetry, Chem. Phys. Lett., 1985, 116, 202-205. (44) Tanikawa, T.; Saito, M.; Guo, J. D.; Nagase, S. Synthesis, Structures and Optical Properties of Trisilasumanene and its Related Compounds, Org. Biomol. Chem., 2011, 9, 17311735. (45) Schwarz F. P.; Wasik, S. P. Fluorescence Measurements of Benzene, Naphthalene, Anthracene, Pyrene, Fluoranthene, and Benzo[e]pyrene in Water, Analytical Chem., 1976, 48, 525-527. (46) Zeng, T; Mellerup, S. K.; Yang, D.; Wang, X.; Wang, S.; Stamplecoskie, K. Identifying (BN)2 pyrenes as a New Class of Singlet Fission Chromophores: Significance of Azaborine Substitution, J.Phys.Chem.Lett., 2018, 9, 2919-2927.

ACS Paragon Plus Environment

27

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

(47) Gusten, H.; Heinrich, G. Photophysical Properties of Fluoranthene and its Benzo Analogues, J. Photochem., 1982, 18, 9-17.

ACS Paragon Plus Environment

28

Page 28 of 73

Page 29 of 73

6 5 ENERGY GAPS IN eV

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

4 3 2 1 0

11 00 00 11 00 11 00 11 00 11 00 11 00 11 00 11 00 11 00 11 00 11 00 11 00 11 00 11 00 11 00 11 00 11 00 11 00 11 00 11 00 11 00 11 00 11 00 11 00 11 00 11 00 11 00 11 00 11 00 11 00 11 00 11

11 00 00 11 00 11 11 00 11 00 00 11 00 11 00 11 11 00 00 11 00 11 00 11

00 11 00 11 00 11 00 11 00 11 00 11 00 11 00 11 00 11 00 11 00 11 00 11 00 11 00 11 00 11 00 11 00 11 00 11 00 11 00 11 00 11 00 11 00 11 00 11 00 11 00 11 00 11 00 11 00 11 00 11

11 00 11 00 00 11 00 11 00 11 11 00 00 11 00 11 00 11 00 11

0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1

11 00 00 11 00 11 00 11 11 00 00 11 00 11 00 11 11 00 11 00

0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1

11 00 00 11 00 11 00 11 00 11 00 11 00 11 00 11 11 00 00 11 00 11 00 11

0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1

11 00 00 11 00 11 00 11 00 11 00 11 00 11 00 11 11 00 00 11 00 11 00 11 11 00

0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1

1 0 0 1 1 0 1 0 0 1 0 1 0 1 1 0 1 0 0 1 1 0 1 0 1 0 0 1 0 1

0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1

1 0 0 1 0 1 1 0 1 0 0 1 0 1 0 1 1 0 0 1 0 1 0 1

0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1

0 1 1 0 0 1 0 1 1 0 1 0 0 1 0 1 0 1 1 0 1 0

One−photon state (BLACK BAR), Two−photon state (BLUE BAR), Lowest triplet state (RED BAR)

ACS Paragon Plus Environment

29

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

[3]

ACS Paragon Plus Environment

Page 30 of 73

Page 31 of 73 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

The Journal of Physical Chemistry

[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

[4]

ACS Paragon Plus Environment

Page 32 of 73

Page 33 of 73 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

The Journal of Physical Chemistry

[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

[5]

ACS Paragon Plus Environment

Page 34 of 73

Page 35 of 73 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

The Journal of Physical Chemistry

[8] 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

[2]

ACS Paragon Plus Environment

Page 36 of 73

Page 37 of 73 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

The Journal of Physical Chemistry

[6]

ACS Paragon Plus Environment

The Journal of Physical Chemistry

3S 2S + 3S 1S

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

+

2S

1S

(5.9, 2.8) (5.6, 0.5)

4.85 4.53

3.84 (1.4, 3.1) 3.77 3.07

3 T2T 1T 3T 2T 1T

+

1S

Page 38 of 73

5.12 4.97

4.09 +

3.26 +

2.78

+

1.96

+

ACS Paragon Plus Environment BENZANTHRACENE

Page 39 of 73

The Journal of Physical Chemistry

3A 3

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

3B 2B 1B

2A 3 1A

(0.0, 2.6) (1.5, 0.0)

5.52 4.60 3.97 3.91 3.23

2A

3

(0.0, 0.7)

5.84 5.19 4.51

+

2B

(0.0, 0.1)

3

+

1B (0.0, 4.3)

1A

7.12

3

+

+

(0.0, 2.5)

+

1A

ACS Paragon Plus Environment

BENZOPYRENE

2.24

The Journal of Physical Chemistry 3

3B 2B

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

+

3A 1B

+

2A

(0.8, 4.2) 5.23 5.02 (8.7, 1.9)

3A 3 2A 3 1A

4.21 3.93 (2.5, 3.3) 3.22

3B 3 + 2B

3

+

3

+

1B

+

1A

+

1A

ACS Paragon Plus Environment

CHRYSENE

(0.2, 0.2)

Page 40 of 73 5.98

5.21 4.99 (1.4, 1.0) 3.64 3.33 (0.3, 0.1) 2.34

Page 41 of 73

3A

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

2B 2A 1B

The Journal of Physical Chemistry

(3.4, 0.0) (0.0, 1.0)

4.24 3.78 3.66 3.45

3

2B 3 3A 3 2A 1

3

A

1 3B

(0.0, 0.4) (0.2, 0.0)

(3.2, 0.0)

1A

1A

(0.0, 0.3)

3.96 3.76 3.57

(4.2, 0.0) ACS Paragon Plus Environment

Fluoranthene

2.87 2.34

The Journal of Physical Chemistry

3B 2B 1B+ 3A

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

(3.7, 0.1)

4.99 (6.9, 0.0) 4.64 4.04 3.90

+

2A

3.16

3A 2A 1A

+

3B + 2B + 1B

(0.0, 0.8)

+

1A

+

1A

ACS Paragon Plus Environment

HELICENE

(4.0, 0.0)

Page 42 of 73

5.46 (0.0, 0.4)

4.98 4.14

(0.3, 0.0) 3.32 2.90 2.46

Page 43 of 73

The Journal of Physical Chemistry

3A 3

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

3B 2B 1B

(0.1, 0.1)

3

(1.2, 0.2) (0.8, 0.0)

4.70

3

4.27 3.67

+

2A

2A 1A

5.50

3.17

3

+

3

+

3B

2B 3

(1.4, 0.1)

4.74 4.19 3.71

3.19 (0.3, 0.3)

+

1B

2.09

(1.3, 0.1)

+

+

1A

1A

ACS Paragon Plus Environment

PICENE

(0.0, 5.3) 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

-

3B2 + 4A 2B3 + 3A1B2 +

2A

(5.1, 0.0)

(0.0, 2.8)

The Journal of Physical Chemistry

5.41 4.94 4.51 3.81 3.75 3.07

43 B2 33 B3 23 B2 3 1 B3 5A 4A 3 A2 A-

(0.0, 3.5)

Page 44 of 73 5.68 5.59 5.30

4.77 4.10 3.88 3.28 3.12 (2.4, 0.0)

(0.0, 5.3)

+

1B

1.86 (2.3, 0.0)

+

+

1A

1A

ACS Paragon Plus Environment

PYRENE

Page 45 of 73 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

The Journal of Physical Chemistry

3 3A

3B 2B 1B+ 3A +

2A

+

1A

(0.0, 4.3)

6.22 (2.4, 0.0)

5.37 5.11

(0.0, 6.4) 4.22 4.10 3.38

3 2A 3 1A 3 + 3B 3 + 2B 3 + 1B

+

1A

ACS Paragon Plus Environment

TRIPHENYLENE

5.16 (1.1, 0.0)

4.36 3.57 3.45 2.83

The Journal of Physical Chemistry

Energy gaps in eV

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

Page 46 of 73

One-photon state Two-photon state Lowest triplet state

5

4

3

2

1

0

h

e en

u Fl

o

t n a r

e n re

Py

ar th e n zacen n e

B

ys

e en

lic

e en

e n e l

e

y y p r n e ACS Paragon zo H he Ch Plus Environment n p i Be Tr

ren

e

P

n e c i

Page 47 of 73 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

0.5

The Journal of Physical Chemistry

0.55

.60

0

0.5

7

0.62

0.6

2

4 5 . 0

0 5 . 0 0. 5

0

0.46

0.47

.46

7

0.5

2 5 . 0 4

0.5

3

[3]

ACS Paragon Plus Environment

9 6 . 0

6

0.71

0.4

7

1 5 . 0 0.53

4 6 . 0

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

8 6 . 0

Page 48 of 73

0.6

0

0.57

.72

0

0.5

0

0.55

0.7

2

0 5 . 0

6 5 . 0 0 .6

0 4 . 0

0.51

0.51

3

0.5

6

2 6 . 0

0.68

6 5 . 0

0.5

6

0.4

3

[3]

ACS Paragon Plus Environment

0.43

8 7 . 0

0.5 6

6

0 .6

0.58 0.64

0 .6

1 0.5

8

0.5

3

0.5

ACS Paragon Plus Environment

5

2 0.6

2

0

0.5

6

0.4

0.72 0.60

[7]

0.6

6

0.5

7

0.55

0.25 0.46

0.52

0.6

0.3

5

0.4

0.51 0.43

7

1 0.6 0.61

0.6

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

The Journal of Physical Chemistry

3

Page 49 of 73

0.49

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

1 6 . 0 0

.51

0.58

Page 50 of 73

3 5 . 0 0

.65

0.49

0 . 4 0 2 1 4 7 4 . 4 . . 0 0 0. 0 0 0 .60 55 .60 0.58 0.55 0.49 5 4 7 5 . . 0 0.7 0 0.4 0 7

[4]

ACS Paragon Plus Environment

0.62

1 6 . 0

0.57

0.5 0

0.7

0.5

5

0

0.5

3

0.6

0.64

0 2 . 6 65 . 0 0.59 0.63 0.45 0.57 0.5 5 2 5 4 . 0 0. 0 0.28

8

2

0.54 0.52

0.5 2

The Journal of Physical Chemistry

0.6

0.5

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

1

Page 51 of 73

[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

.58

0

0.70

Page 52 of 73

0.7

3 6 . 0

0.55

0.55

0

0.6

1 0.63

0.4

0.5

4

5 5 . 0

0.47

5 0.48 0.

47

0.60

0.47

0.39

0.7

6

0 8 . 52 4 . 0

[5]

ACS Paragon Plus Environment

7 5 . 0

7 6 . 0

0.65

9

0.63

0.75

[8]

ACS Paragon Plus Environment

0.5

0.5

2

8

0.7

0

4 0.5 9 0.4

9 0.5

0.5

0.40

2 0.5

7

0.5

7

0.4

0.41 0.33

9 0.5

1 0.6

0.42

9 0.5

5 0.5 0.54

0.65 0.48

The Journal of Physical Chemistry

0.6

0.69

9 0.6

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

8

Page 53 of 73

0.52

0.6

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

5 6 . 0

Page 54 of 73

5

0.58 5 4 . 0 0.5

0.58 0 . 45 3 5 . 0 3 0.77 0.77 0.49 0.64 0.64 0.55 0.4 9 0.5 9 3 4 5 3 . . 0 0 0.48 0.48

0.6

4

[2]

4 6 . 0

ACS Paragon Plus Environment

Page 55 of 73

0.5 5

0.4

7

5

1

0.65

0.67

[6]

ACS Paragon Plus Environment

0.5

7

8

0.5

0.6

0.5

5

7

0.37 0.46

0.6

0.5

0.6

0.55

3

0.3 7

0.57 0.43

6 0.4

0.58

5 0.6

0.6 7

0.61 0.57

8 0.5

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

The Journal of Physical Chemistry

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

4 6 0.

Page 56 of 73

0.6

0.54

9 6 . 0 0.

0 0.63

6 9 5 4 . . 51 0 0.54 0 0.47 6 5 . 0 0.54 0.40 0.38 0.48 7 3 6 6 5 5 . 0. 0.5 0. 0.5 0 0.6 8 3 4

[3]

ACS Paragon Plus Environment

Page 57 of 73 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

The Journal of Physical Chemistry

3 6 . 0

0.5

9

0.53

6 6 . 0 0.5 3 0.54 4 5 . 0 0.6 5

0.61

4 5 . 0 0

0 5 . 0 0

6 5 . 0

0.38

0.36

0.49

.52

0.5

3

.45

3 5 . 0 0

.55

[3]

ACS Paragon Plus Environment

5 6 . 0

0.62 0.58

2 0.5

1

0.5

1 0

0.6

0.40 0.60

[7]

ACS Paragon Plus Environment

2 0.6

5

0.5 7

0 .6

0.4 6

9

0.5

8

0.4

0.39 0.31

0.57

0.5

0.6

5 0.47 0.41

0 0.4

7

0.5

2

0.46

Page 58 of 73

0 .5

5 0.5 0.48

0.6

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

7

The Journal of Physical Chemistry

0.53

Page 59 of 73 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

The Journal of Physical Chemistry

1 7 . 0

0.4

8

0.45

3 5 0. 0.43

0.5 5 4 7 2 4 4 6 . . 0 0 0 0. 0.5 0.5 .53 4 5 0.56 0.45 0.43 2 3 7 . 5 0 0.6 0. 0. 47 3

[4]

ACS Paragon Plus Environment

0.5

9 0.54

0 6 . 0

4

0.4

The Journal of Physical Chemistry

0.4

0.5

8

7

0.5

4

0.5

0.53

Page 60 of 73

0.6

4

0.56 0.55

0.4

5

3

0.49 0.28

0.32 0.61

0.5

2

0.5

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

9

0.72

5 0.5 3 . 0 9

[1]

ACS Paragon Plus Environment

0.63

0

0.60 .67

Page 61 of 73 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

The Journal of Physical Chemistry

4 6 .

0

0.59

0.5

5 0.48

0.5

4 6 . 0

0.56

0.54

9

0.5

6 0.64

0.4 .51

0

6 0.50 0.

42

0.43 0.39

0.6

1

0 3 . 53 5 . 0

[5]

ACS Paragon Plus Environment

4 5 . 0 0.52

1 6 . 0

[8]

ACS Paragon Plus Environment

0.6

4 0.67

0.57

0.5

0.4

8

0.4 4

0.5

1

0.5

0.6 3

3 0.5

0.47

2

4 0.6

0

9 0.4

7 0.6

0.53

0.29 0.44

8 0.5

0.39

Page 62 of 73

6 0.4

0.6

0.66

4 0.5

9 0.5

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

0.37 0.52

The Journal of Physical Chemistry

0.55

4

0.65

Page 63 of 73 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

0.6

The Journal of Physical Chemistry

4 6 . 0

4

0.43 0 1 .59 5 . 0 .51 0.56 0.56 0.52 0.63 0. 0.68 0.68 42 0.5 2 0 4 5 0 . . 0 0 0.54 0.54 0.43 9 5 . 0 0

0.5

6

6 5 . 0

[2]

ACS Paragon Plus Environment

0.58 0.56

0.5 3

3 0.4

3

6

1 0.61

0.63

[6]

ACS Paragon Plus Environment

0.5

0.4 5

3

8

0.5

0.5

0.5

0.45 0.49

3

0.53

0.6

7

0.4

0.6

0.47 0.43

9 0.4

0.53

Page 64 of 73

3 0.6

3 0.5

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

0.6 1

The Journal of Physical Chemistry

Page 65 of 73 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

The Journal of Physical Chemistry

0.04

0.11

0.03

0.06

0.05

0.03

0.05

0.04

0.06

0.08

0.06

0.05

0.06

0.03

0.10

[3]

ACS Paragon Plus Environment

0.03

0.03 0.09

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

0.04

0.04

0.09

0.09

0.03

0.03 0.07 0.00 0.14

Page 66 of 73

0.08 0.08 0.02

0.02

0.00 0.00 0.02

0.02

[7]

ACS Paragon Plus Environment

0.07 0.00 0.14

Page 67 of 73 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

0.00

The Journal of Physical Chemistry

0.08

0.04 0.10

0.19

−0.02 0.19

0.02 0.04 0.05

0.06

0.02

−0.02

0.06

0.05

0.08

[4]

ACS Paragon Plus Environment

0.10

0.00

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

0.04

0.20 −0.03

Page 68 of 73

0.07 0.17 −0.02

0.03 0.01 0.01

0.20 0.04

0.17 0.07

[1]

ACS Paragon Plus Environment

0.03

0.02 0.02

Page 69 of 73 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

The Journal of Physical Chemistry

0.03

0.03 0.09

0.04 −0.00

0.09

0.09

0.04

−0.00

0.13

0.01

0.09 0.01

0.04

0.12 0.05

0.05

[5]

ACS Paragon Plus Environment

0.12

0.05 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

−0.02

−0.00 0.04

0.06 −0.02 0.16

0.15

0.15

−0.00

The Journal of Physical Chemistry

−0.00

−0.00

0.07

0.04

0.07 −0.02

0.02

0.02

[8]

ACS Paragon Plus Environment

0.16

0.05

Page 70 of 73

−0.02 0.06

−0.05

Page 71 of 73 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

The Journal of Physical Chemistry

0.18

0.18

0.01

0.01

0.00

0.09

0.09

0.09

0.09 0.01

0.00

0.18

0.01 0.18

−0.05

[2]

ACS Paragon Plus Environment

0.05

0.05

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

0.07

0.07 0.05

0.05 0.07

0.05 0.05 0.05

0.05 0.05

Page 72 of 73

0.07 0.05

0.05 0.07

0.07

[6]

ACS Paragon Plus Environment

0.05

Page 73 ofThe 73 Journal of Physical Chemistry

6

ENERGY GAPS IN eV

5

1 2 3 10 3 11 00 11001 11 00 00 11 11 10 00 010 00 11 00 11 24 10 00 11 1010 00 11 11 00 11 00 10 00 11 11 00 101 11 00 11 00 11 00 1010 11 10 00 00 11 0 11 00 00 11 00 Paragon 11 11 1 00 0 ACS Plus Environment 00 11 1 0 11 00 00 11 5 00 11 1 11 10 00 11 00 1 0 00 11 11 00 1 11 00 1001 00 11 0 1 00 11 1 0 11 00 00 11 00 11 010 11 1 00 0 11 00 101 00 11 11 00 11 00 6 11 00 101 00 11 0 11 00 00 11 00 11 10 11 00 0 11 00 1 0 00 11 11 00 11 00 101 00 11 0 1 00 11 1 0 11 00 11 00 07 11 00 0 00 11 10 11 00 10 One−photon state (BLACK BAR), Two−photon state (BLUE BAR), Lowest triplet state (RED BAR) 4

00 11 00 11 00 11 00 11 00 11 00 11 00 11 00 11 00 11 00 11 00 11 00 11 00 11 00 11 00 11 00 11 00 11 00 11 00 11 00 11 00 11 00 11 00 11 00 11 00 11 00 11 00 11 00 11 00 11 00 11 00 11 00 11

00 11 00 11 00 11 00 11 00 11 00 11 00 11 00 11 00 11 00 11 00 11 00 11 00 11 00 11 00 11 00 11 00 11 00 11 00 11 00 11 00 11 00 11 00 11 00 11 00 11 00 11 00 11 00 11 00 11 00 11

0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1

0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1

0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1

0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1

0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1

0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1