Tuning the HOMO–LUMO Energy Gap of Small Diamondoids Using

Feb 20, 2017 - All calculations were performed within the Kohn–Sham framework using the Gaussian09 software package.(40) During the inverse design s...
5 downloads 13 Views 7MB Size
Subscriber access provided by University of Newcastle, Australia

Article

Tuning the HOMO-LUMO energy gap of small diamondoids using inverse molecular design Jos L. Teunissen, Frank De Proft, and Freija De Vleeschouwer J. Chem. Theory Comput., Just Accepted Manuscript • DOI: 10.1021/acs.jctc.6b01074 • Publication Date (Web): 20 Feb 2017 Downloaded from http://pubs.acs.org on February 23, 2017

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.

Journal of Chemical Theory and Computation 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 25

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

Journal of Chemical Theory and Computation

Tuning the HOMO-LUMO energy gap of small diamondoids using inverse molecular design Jos L. Teunissen, Frank De Proft, Freija De Vleeschouwer* [email protected] Research Group of General Chemistry, Vrije Universiteit Brussel (VUB) Pleinlaan 2, 1050 Brussels (Belgium)

Abstract Functionalized diamondoids show great potential as building blocks for various new optoelectronic applications. However, until now, only simple mono- and double substitutions were investigated. In this work, we considered up to 10 and 6 sites for functionalization of the two smallest diamondoids, adamantane and diamantane, respectively, in search for diamondoid derivatives with a minimal and a maximal HOMO-LUMO energy gap. To this end, the energy gap was optimized systematically using an inverse molecular design methodology based on the best-first-search algorithm combined with a Monte Carlo component to escape local optima. Adamantane derivatives were found with HOMO-LUMO gaps ranging from 2.42 eV to 10.63 eV, 9.45 eV being the energy gap of pure adamantane. For diamantane, similar values were obtained. The structures with the lowest HOMO-LUMO gaps showed apparent push-pull character. The push character is mainly formed by sulfur or nitrogen dopants and thiol groups, whereas the pull character is predominantly determined by the presence of electron-withdrawing nitro or carbonyl groups assisted by amino and hydroxyl groups via the formation of intramolecular hydrogen bonds. In contrast, maximal HOMO-LUMO gaps were obtained by introducing numerous electronegative groups.

Introduction Nanoparticles play an increasingly important role in modern technology.1,2 Especially carbon nanoparticles like fullerenes and nanotubes are nowadays widely used. Diamondoids, however, are still only scarcely assessed. Diamondoids or nanodiamonds are a class of carbonnanoparticles that consist of fused diamond cages and are fully sp3-hybridized, very stable3, and perfectly size-selectable. The first members of the series of diamondoids, having respectively one, two and three cages, are adamantane, diamantane and triamantane. Since small diamondoids have been extracted from natural oil only recently,4 there is a growing research interest in using diamondoids. Among others, diamondoids find possible applications in the field of optoelectronics. Functionalized diamondoids can for example be used for bioimaging1,5, when the optical gap would be reduced to the visible region. On the other hand, diamondoids have their value in semi-conductors or other electronic devices, when the electrical conductivity is increased by functionalization.6 Moreover, diamondoids can form monolayers on gold surfaces7 or self-assemble into more complex structures.8 Pure diamondoids have, however, in accordance with pure diamond, a high hardness9 and a very low conductivity.2 Henceforth, for

1 ACS Paragon Plus Environment

Journal of Chemical Theory and Computation

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 2 of 25

useful behavior, one tries to manipulate pure diamondoids to exhibit enriched functionality, for example by introducing other chemical substituents into diamondoids. Several types of functionalization have been applied: diamondoids with fluoroalkyl groups show negative electron affinities and electron-phonon scattering,10 amine-functionalized diamondoids display low enhanced conductivity11, and thionated diamondoids have low optical gaps.1 All these properties are largely determined by the position of the frontier molecular orbitals. Hence, in functionalization studies generally the effect of doping and substitution on electronic properties such as the ionization potential, the electron affinity, and the HOMO-LUMO energy gap is examined,12 in addition to their function as binding sites to make usable nanostructures.2 In particular, some initial studies showed that push-pull doping13,14 is the most effective approach to decrease HOMO-LUMO gaps and that the most electron-withdrawing functional group dominates this reduction.15 Fokin and Schreiner also found that group effects are not additive, meaning that the final property value is not simply an summation of individual contributions from substituents or dopants; therefore, double functionalization sometimes shows rather little effect. Moreover, these authors concluded that the nature of the functional groups seems to be significantly more crucial than the size of the diamondoids itself.15 Until now, most studies focused only on mono- and double functionalization of diamondoids. To explore the world of diamondoids to a larger extend, we thoroughly investigated to what degree the electronic properties of adamantane and diamantane can be tuned by multiple functionalization. Since the number of possible chemical structures rises vastly when more sites for functionalization are taken into account, an inverse molecular design approach was employed to search systematically for structures with ideal properties. We probed the boundaries of the HOMO-LUMO energy gap by searching for their minimal and maximal values in a given chemical space of diamondoid derivatives. In addition, a detailed evaluation of the functionalization effect on the different properties has been carried out.

Computational Methodology In the last decade growing attention has been paid to the design and discovery of valuable molecular structures, using quantum chemistry. A fundamental challenge in molecular design and discovery is the large number of possible structures that are accessible through the systematic variation of the composition of the molecular system.16-18 Direct approaches in molecular design are e.g. the systematic modification of structural motives with known favorable chemical properties, quantitative structure–activity relationship methods, and virtual high-throughput-screening of existing molecular libraries.19-23 In this work, however, an inverse molecular design methodology24-32 is applied to design diamondoid derivatives displaying optimal HOMO-LUMO energy gaps. In this methodology, one searches for an optimal external potential, i.e. the potential due to the nuclei of the system, corresponding to a molecular system with the associated target property. One of the most straightforward search methods is the discrete best-first-search (BFS) algorithm.33 The algorithm evaluates the property value of systematically chemically modified structures and then adopts the most optimal modification. The best first search (BFS) algorithm is known to be quite robust and efficient for single property optimizations, when the search space is of moderate size. In those cases, it generally performs better than the genetic algorithm (GA) method and the dead-end elimination (DEE) method.30 Balamurugan et al.30 have shown that the same performance of GA can be achieved with BFS using 20 times lower computational cost for a system with 12 available sites and a substituent library of 13, which results in a chemical space size of approximately 2.3 1013.

2 ACS Paragon Plus Environment

Page 3 of 25

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

Journal of Chemical Theory and Computation

In this study, the imposed structural changes are the placements of various functional groups (substitution or doping) on certain atomic positions, so-called sites, within an initially determined molecular framework. On each site, the different functional groups are ranked according to their corresponding property value and, subsequently, the most optimally modified structure is used as a new starting point for the next site to visit. Convergence is reached when the current optimum cannot be improved anymore by modifying a single site. This procedure only converges to a global optimum when the independent site approximation more or less holds. The independent site approximation assumes that the sites can be optimized individually, meaning that the global property value can be written as a sum of individual contributions from functional groups on each of the sites. More details about the algorithm can be found in literature.30,34,35 However, in many cases the independent site approximation is not perfectly valid and, therefore, one could end up in local optima. The simplest solution is to run the design procedure multiple times, using a different starting structure and/or site order.36 However, this way, costly computer time could be wasted on structures that are far from optimal. In addition, variations on the BFS procedure have been proposed that allow to move away from the local optimum by introducing randomization in the early stages of the optimization process, with the BFS converged solution as the new initial structure.30,37 Hence, based on the work of Hu, Beratan and Yang,38,39 we have adapted their gradientdirected Monte Carlo approach to the discrete best first search algorithm. This approach works as follows. When a complete BFS procedure has been performed and converged, we have available, for a given optimal structure, also all the structures with only one site changed. That is, for a converged optimum, 𝑋opt 𝑏! , 𝑏! , … , 𝑏! , … , 𝑏! , with bi the functional group on site i, all the property values of structures with only one site changed, 𝑋!! 𝑏! , 𝑏! , … , 𝑥! , … , 𝑏! , are known. When the independent site approximation would be valid, one can calculate the property value P of every new structure, 𝑋 ! 𝑥! , 𝑥! , … , 𝑥! , … , 𝑥! with possibly N sites modified with respect to the converged optimum, by the property value for the Xopt structure and a difference term: 𝑃 𝑋 ! = 𝑃(𝑋opt ) + ∆𝑃(𝑋 ! )

(1)

This difference term is the sum of the difference in property value between the converged optimum and the structures with only one site of the optimal configuration changed to the functional group present in XN: ∆𝑃(𝑋 ! ) =

! !!!

𝑃 𝑋!! 𝑏! , 𝑏! , … , 𝑥! , … , 𝑏!

− 𝑃 𝑋opt 𝑏! , 𝑏! , … , 𝑏! , … , 𝑏!

(2)

To escape from possible local traps, after each complete BFS convergence, a random structure XN is generated and accepted with a probability, 𝑝 ∆𝑃 𝑋

!

=𝑒

!

∆! !! !

= 𝑒 !!∆!

(3)

using the Metropolis algorithm. kB is the Boltzmann constant and T is the temperature. The temperature controls the acceptance probability and is increased slowly when no structures are accepted. All calculations were performed within the Kohn-Sham framework using the Gaussian09 software package.40 During the inverse design search, all geometries were optimized at the B3PW91/6-31G(d,p)41,42 level of theory and the same computational methodology was used to obtain the HOMO and LUMO energies, in accordance with the work of Fokin and Schreiner.10 It

3 ACS Paragon Plus Environment

Journal of Chemical Theory and Computation

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 4 of 25

was shown that the B3PW9141 density functional outperforms B3LYP in describing polycyclic fully sp3-hybridized structures,43 while B3PW91 and B3LYP show similar good agreement in describing hydrogen bonds.44-46 In addition, the B3PW91 functional outperforms all other density functionals in computing the band gap for inorganic semiconductors.47 In all calculations, 631G(d,p)42 was used as basis set. No diffuse functions were added, despite the spread-out character of the lowest unoccupied molecular orbital, being delocalized around the surface of the diamondoids. Diamondoids, however, often possess negative electron affinities48-50 and for these cases it has been shown that the LUMO energy shows a much poorer agreement with experimentally established trends when diffuse functions are included in the basis set.51,52 To validate the computational approach, the adiabatic ionization potential was calculated as an energy difference between the neutral diamondoid and the cation. B3PW91/6-31G(d,p) reproduces the experimentally observed adiabatic ionization potential reasonably well: 9.23 eV vs. 8.96 eV for adamantane and 8.80 eV vs. 8.43 eV for diamantane.53 Furthermore, we have tested the B3PW91 functional for an additional 15 substituted small diamondoids. The results are listed in the Supporting Information. On average, B3PW91/6-31G(d,p) computed IP values underestimate the experimental values by 0.39 eV. Using the same functional but a larger basis set, B3PW91/6-311++G(d,p), the IP difference decreases to 0.24 eV. Correlation diagrams between experimental and computed values result in correlation coefficients R2 of 0.87 for the smaller basis set and 0.93 for the larger basis set.

Results and discussion To start the inverse design search, one has to establish the properties of interest and the molecular frameworks, with indication of the sites and functional groups that will be used. In this study, the target property is the HOMO-LUMO energy gap. For the molecular framework, the two smallest possible diamondoids, adamantane and diamantane, are chosen. In adamantane the position of each carbon atom is treated as an adjustable site, with four tertiary and six secondary carbon positions, whereas for diamantane only six out of fourteen carbon atoms are handled as a site for functionalization to keep our chemical compound space and computational cost manageable. Diamantane has three symmetrically different carbon positions: two apical tertiary atoms, six medial tertiary atoms, and six secondary atoms. Of each type, two of them are used as a site for functionalization. The scaffolds with indication of the tunable sites are shown in Figure 1. The set of functional groups consist of three groups for internal doping: (N), (O), (S), and twelve groups for external doping or substitution: (C)-COOH, (C)-NO2, (C)-CN, (C)-OH, (C)SH, (C)-NH2, (C)-CH3 (C)-CF3, (C)-Cl, (C)-F, (C)-H, (C)=O. This means that there are twelve functional groups (11 substituents and 1 dopant) available for tertiary carbon positions in the molecular framework and an additional three for the secondary positions (substituent (C)=O and dopants (O) and (S)), by removing also the other hydrogen atom on that site. All of these functional groups have been already used in the synthesis of mono-substituted diamondoids.54-56 Remark that in principle one could introduce two functional groups on secondary carbon atoms, but they are more difficult to synthesize and only sparse examples are found in literature. In this work, we focused on neutral functional groups, but one could also obtain low HOMO-LUMO gap structures by introducing charged groups.15 Based on the number of sites and functional groups as listed above, the investigated chemical compound spaces of adamantane and diamantane derivatives are respectively 2.4*1011 and 4.6*106 (neglecting symmetrically identical molecules).



4 ACS Paragon Plus Environment

Page 5 of 25

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

Journal of Chemical Theory and Computation

1 5

2

6

5

7 1 8

3

4

10 9

(a)

3 4 6 2

(b)

Figure 1. Molecular framework of (a) adamantane and (b) diamantane, with indication of the adjustable sites considered in our work. In adamantane, there are 4 tertiary (1, 2, 3, 4) and 6 secondary (5, 6, 7, 8, 9, 10) positions, whereas in diamantane, 2 apical tertiary (1, 2), 2 medial tertiary (3, 4), and 2 secondary (5, 6) sites are chosen.

Table 1. The HOMO and LUMO of pure adamantane and diamantane, with indication of their corresponding orbital energies, are depicted (isovalue: 0.02). In addition, the HOMO-LUMO energy gap (HLg) and the vertical and adiabatic ionization potential (IPvert, IPad) are given in eV.

Diamondoid

Structure 2

6 1 8

3

LUMO

HLg

IPvert

IPad

9.45 eV

9.29 eV

8.96 eV

9.00 eV

8.92 eV

8.43 eV

5

7

Adamantane

HOMO

4

10 9

εHOMO =-7.43 eV

εLUMO = 2.02 eV

1 5 3

Diamantane

4 6 2

εHOMO = -7.02 eV

εLUMO = 1.98 eV

A risk of using automated algorithms is to lose the insight of the underlying basic principles. Therefore, in order to gain initial understanding of complicated diamondoid derivatives, we start by considering the frontier orbitals of the unsubstituted and mono-substituted adamantane and diamantane structures. The HOMO and LUMO of pure adamantane and diamantane are delocalized over the whole molecule, with the LUMO showing Rydberg state character.57 The HOMO is mainly associated with the carbon-carbon bonds, whereas the LUMO shows anti-bonding character on the carbon-hydrogen bonds and is more diffuse and delocalized over the surface of the whole diamondoid. The orbitals are depicted in Table 1. The LUMO energies for both diamondoids are positive, indicating that these molecules probably have negative electron affinities. Both the HOMO-LUMO energy gap and the ionization potential are around 0.5 eV smaller for the larger diamantane, compared to adamantane. This is in accordance

5 ACS Paragon Plus Environment

Journal of Chemical Theory and Computation

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 6 of 25

with the trend that the energy gap of increasing diamondoid size converges towards the bandgap value of approximately 5.0 eV for diamond. Functionalization of diamondoids by doping and substitution changes the electronic properties. In most cases this is attributed to the placement of localized electronic states into the HOMO-LUMO gap, either occupied or unoccupied.14 We have computed the HOMO-LUMO gap for all possible mono-substituted adamantane and diamantane structures, with functionalization on the different types of carbon positions, as shown in Figure 2. Sulfur doping (thia) or substitution with a carbonyl, nitro, thiol, or carboxyl group introduces both a new HOMO and a new LUMO level into the HOMO-LUMO gap of the pure diamondoids. These groups therefore result in the most significant decrease of the HOMO-LUMO gap. The largest reduction of 3.5 eV by a single functionalization occurs by incorporating a carbonyl or nitro group. Nitrogen doping (aza) and amino groups only introduce a new HOMO level, whereas a new LUMO level only appears when chloro or cyano groups are included. An example of each case is depicted in Figure 3. When visualizing the orbitals, it is observed that the HOMOs of the amine, alcohol, and thiol compounds are the most strongly localized on the lone pair of the functional group. Aza/thia doping, on the other hand, gives rise to HOMO levels that are largely located on the dopant but also still partly delocalized in the cage. Fluor and trifluoromethyl groups, on the contrary, stabilize the pure diamondoid’s HOMO and thus increase the HOMOLUMO gap.

Figure 2. HOMO-LUMO energy gap, LUMO and HOMO energy level values (in eV) for monosubstituted diamondoids. The values for adamantane (with distinction for secondary and tertiary sites) are given in green whereas for diamantane (with distinction between secondary, apical tertiary, and medial tertiary sites) are blue.

Mono-substituted adamantanes and diamantanes show very similar behavior with the only difference that, when no new HOMO levels are introduced in the gap, the HOMO levels of diamantane are higher in energy.

6 ACS Paragon Plus Environment

Page 7 of 25

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

Journal of Chemical Theory and Computation

Effects of substitution on different positions in the molecule are very small. Most notable is that nitrogen doping on secondary positions results in a lower-lying HOMO level than doping on a tertiary position. This complies with the fact that upon electron removal the positive charge on tertiary nitrogen atoms is more stabilized by hyperconjugation than on secondary nitrogen atoms.

Figure 3. MO-diagram showing the introduction of new HOMO and LUMO levels inside the HOMOLUMO gap of pure adamantane. 2-Aminoadamantane introduces a new nN-HOMO level. 2-Cyanoadamantane introduces two new, almost degenerate π*C ≡N anti-bonding π-type LUMO levels. Adamantanone introduces both a new HOMO and a new LUMO level.

Recently, a few studies on multiple substitutions have shown that the most significantly decreased HOMO-LUMO gap is obtained when the HOMO and LUMO are located on different functional groups.14,15 Remark, however, that group effects are not additive. For example, placing one nitro group on adamantane reduces the HOMO-LUMO gap from 9.44 to 6.15 eV. Introducing another nitro group further reduces the gap to 6.01 eV, while 10 nitro groups finally decrease the HOMO-LUMO gap to 5.64 eV, merely an additional 5% lower compared to the reduction obtained with a single nitro substitution. Hence, to obtain a larger closing of HOMOLUMO gaps, one has to establish diamondoids with different groups that cooperate together. Because these optimal multiple functionalized diamondoids are much less intuitive to find, we turn to the Inverse Molecular Design (IMD) approach.

Inverse molecular design: minimizing the HOMO-LUMO energy gap We have explored the limits of the HOMO-LUMO energy gap of adamantane and diamantane derivatives by a minimization and maximization via the best first search algorithm, using the substituents and sites as mentioned in the previous section.

7 ACS Paragon Plus Environment

Journal of Chemical Theory and Computation

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 25

We first discuss the results for the HOMO-LUMO gap minimization of adamantane. Two searches have been performed, both converging to a minimum within 4 global iterations (over all sites and substituents). This means that, without considering additional Monte Carlo steps, less than 500 structures were visited in each search, which is much smaller than the size of the investigated chemical space (2.4*1011). A first BFS optimization started from a random start configuration having a HOMO-LUMO gap of 4.7 eV, already well below the energy gap for the unsubstituted adamantane. This can be attributed to the fact that the configuration already contained a carbonyl and a nitro group, as well as a nitrogen dopant (see optimization pathways in Supporting Information). Therefore, from the beginning push-pull character is present. The first BFS procedure resulted in an energy gap minimum of 2.85 eV. Subsequent BFS runs, connected via two Monte Carlo steps, further decreased the HOMO-LUMO gap to 2.75 and 2.67 eV, respectively. The resulting structures are depicted in Table 2, whereas the different chemical modifications along the entire optimization process and plots of the frontier orbitals of the structures in Table 2 can be found in the Supporting Information. In the three optima, the LUMO is almost completely localized on the nitro group, assisted by the neighboring hydroxyl and amino groups through intramolecular hydrogen bonding. In fact, a chain of weaker to stronger hydrogen bonds is present in all optima, and this is confirmed by plotting the non-covalent interactions in the structures using NCIplot58,59 (cf. section 8 in the Supporting Information. The HOMO, on the other hand, is mainly located on the sulfur dopant and to a lesser degree on the nearby lying amino groups (on sites 7 and/or 9) in optimal configuration 1 and 2, whereas in the final configuration, ada_min_gap1, the HOMO is located both on the sulfur dopant on site 8 and the thiol group on site 6. Note that in all cases the HOMO and LUMO are spatially separated. All in all, the three minima look very alike, hence the similar energy gap values. We have also depicted the convergence behavior (Figure 4(a)). In the first global iteration, the gap reduces continuously to approximately 3.2 eV. In the three subsequent iterations, the gap decreases further though very slowly to 2.85 eV. Next, a Monte Carlo step is performed. Within one global iteration, already a slightly better energy gap value is obtained, that is further reduced in the next iterations. After another MC step, it takes another three global iterations to find ada_min_gap1. Note that in 8 of the 10 global iteration sections, the spread in HOMO- LUMO gap is very large when site 5 is changed. Site 5, containing the nitro group in the minimal gap configuration, is by far the most important site in this structure, with energy gap values ranging over almost 4 eV. This is in accordance with the results for the mono-substituted adamantane derivatives. Table 2. The optimal chemical structures obtained from three successive best first search runs, connected via Monte Carlo steps (MC1 and MC2), of adamantane with a minimal HOMO-LUMO gap (HLg). More details can be found in the Supporting Information.

Numbered structure

Optimal configuration 1 OH

2

6

5

7 1 8

3

4

10 9

NH 2

H2N

OH

NH 2

CH 3 OH

H3 C S H3 C

CF3

O N

HO O

CH 3

H2 N

OH

HLg = 2.85 eV

NH 2

HLg = 2.75 eV

O N

HS O

S H3 C NH 2



Final configuration (after MC1) (ada_min_gap1)

O N

O

Optimal configuration 2 (after MC1)

O

F

CH 3

H2 N

OH S H2 N OH

HLg = 2.67 eV 8

ACS Paragon Plus Environment

Page 9 of 25

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

Journal of Chemical Theory and Computation

(a)

(b)

Figure 4. Chronological order of the BFS for minimal HOMO-LUMO gap adamantane derivatives: path to ada_min_gap1 (a) and ada_min_gap2 (b). The black lines separate the global iterations. Per separate color the lower limit represents the current structure with the minimal HOMO-LUMO gap. The site order varies in every global iteration. In (a), two MC steps were performed, dividing the full optimization pathway in 3 BFS convergence sections with different number of global iterations (cf. numbering in plot).

Since the initial structure of the above-described BFS optimization already had a quite low energy gap, we performed a second BFS optimization. This time, the random start configuration, as shown in Table 3, had a HOMO-LUMO gap of 6.5 eV, and did not contain nitro, carbonyl or thiol substituents. However, in the first four sites that are visited immediately two nitro groups (sites 7 and 4), a thiol group, and a sulfur dopant were placed (site 6 and 8, respectively). These modifications result in a largely decreased HOMO-LUMO gap of 3.81 eV. At this point, pushpull character is again present: the LUMO is localized on the nitro group on site 4, while the thiol and sulfur yield a high HOMO level. In the following two global iterations, the HOMO-LUMO gap further reduces to 2.42 eV (cf. Figure 4(b)). The different chemical modifications along the complete optimization process are given in the Supporting Information. In this particular case, the use of Monte Carlo steps after the first BFS convergence did not result in a lower gap configuration. As is clear from Figure 4, the energy gap variations within a particular site are far less outspoken for the optimization to ada_min_gap2 than for that to ada_min_gap1, even though both minima hold a nitro group. In fact, 8 out of 10 functionalizations are the same in both structures, albeit on different positions. It seems that in ada_min_gap2 more sites are significantly contributing to the minimization of the energy gap.



9 ACS Paragon Plus Environment

Journal of Chemical Theory and Computation

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 25

Table 3. Three relevant chemical structures obtained from a best first search for an adamantane derivative with a minimal HOMO-LUMO gap (HLg).

Numbered structure

Random start configuration

Optimum after 1st global iteration O

4 10

9

NH 2

HO

5 3

6

H3 C

7

CF3

O

8

NH 2

HO H2 N

SH CF3

HS S H2 N

CN CN

O N

HN

1 2

O

N

Cl

O

HLg = 6.51 eV

HO HO

O

HLg = 3.62 eV

O

S

H2 N

NH 2 H2 N

N

Cl

Final configuration (ada_min_gap2)

SH

O NH 2

HLg = 2.42 eV

The final configuration, as can be seen in Table 3, has no symmetry. Nevertheless, some local symmetry can be found around the thiol group as the two adjacent sites (1 and 2) are substituted with amino groups and around the nitro group as two of the neighboring sites hold hydroxyl groups. Notably, in this optimum the nitro group is placed on a tertiary site, whereas for ada_min_gap1 NO2 is positioned on a secondary site. The frontier orbitals are shown in Table 4. The HOMO, with an energy level of -5.51 eV, is mainly formed by the non-bonding electron pairs of the thiol group and to a lower extent of the amino group on site 7. These thiol and amino groups are electron-donating and hence responsible for the pushing character by increasing the electron density of the structure. Also the amino groups on sites 1 and 2, seemingly not participating in forming the HOMO, have a considerable impact on increasing the HOMO’s energy level. The LUMO of ada_min_gap2 (-3.09 eV) is a π*-MO located on the nitro group. The nitro group is electron-withdrawing and hence removes electron density from the cage. The π*-MO LUMO of nitroadamantane lies at -1.68 eV. By introducing the neighboring hydroxyl groups, the LUMO energy is further decreased to -2.43 eV. However, the LUMO of ada_min_gap is not located on those neighboring groups. The hydroxyl groups seem to assist the nitro group through intramolecular hydrogen bonding, as we also observed for the ada_min_gap1 structure. Again, we have computed the non-covalent interactions using NCIplot.58,59 In Figure 5, we clearly see two strong hydrogen bonds between the nitro group and the adjacent hydroxyl groups, as well as a weaker attractive interaction between the sulfur dopant and the nitro group. Many other weak attractive interactions are present, as is displayed in an additional NCI plot in the Supporting Information. We have also looked into the site vs. functionalization dependency of the HOMO-LUMO energy gap w.r.t. the optimal structure. One can appreciate the importance of each site or the functionalization on that site by looking at how the property value changes when only that particular site is changed, while the rest of the optimum is kept fixed. The last cycle in Figure 4 shows the variation in energy gap values for the last global iteration, in which no more chemical modifications to the optimum configuration take place. In the Supporting information more detailed plots are included. We will briefly discuss for ada_min_gap2. Site 4, containing the nitro group, is the most important site in this structure, because modifying this group while keeping the other sites fixed induces a minimal increase in HOMO-LUMO gaps of 1.5 eV. Most other sites display more or less the same range in energy gaps, though with a different functionalization preference. Only site 3 seems to be less influential, with many modifications hardly affecting the

10 ACS Paragon Plus Environment

Page 11 of 25

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

Journal of Chemical Theory and Computation

HOMO-LUMO gap. Note that, if the independent site approximation (ISA) would be valid, similar sites would show a similar property vs. substituent trend. In this case, however, the behavior of all the functionalizations per site is very different. This deviation from ISA can be attributed to the fact that in ada_min_gap2 the frontier orbitals become quite localized on the dopants and/or functional groups. Table 4. Optimal structures for the HOMO-LUMO gaps of adamantane and diamantane derivatives.

Framework

Structure

HOMO

LUMO

HLg (eV)

ada_min_gap2 O

O

N

4 10

9 5

1

HO HO

S

2.42

3 8

2 6

H2 N

NH 2

7

O

H2 N

NH 2

SH

εHOMO = -5.51 eV

εLUMO = -3.08 eV

dia_min_gap H

1

O

O

5 N O

3

2.94

H

4

O N

6 2 N

H3 C

εHOMO = 5.57 eV

εLUMO = -2.64 eV

ada_max_gap F F

2

6 7 1 8

3

F

5

10.63

4

10

F

9

F

O O

O F

εHOMO = -8.70 eV

εLUMO = 1.93 eV

dia_max_gap CF3

1

F

5 3

10.08

CF3

CF3

4 6 2

F CF3



εHOMO = -8.45 eV

εLUMO = 1.63 eV

11 ACS Paragon Plus Environment

Journal of Chemical Theory and Computation

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 25

Figure 5. Non-covalent interactions (with the electron density ρ ranging from 0.020 to 0.050 a.u.; the reduced density gradient s cutoff = 0.5 a.u. and NCI color scale = -0.03 < ρ < 0.03 a.u.) for the ada_min_gap2 structure using NCIplot.59 Blue surfaces indicate strong attractive interactions such as hydrogen bonds.

Next, we turn to the energy gap minimization of diamantane. Our lowest HOMO-LUMOgap diamantane derivative was found after 6 global iterations with one Monte Carlo step in between, the summary of the procedure being depicted in Figure 6. Subsequently, two other Monte Carlo steps and BFSs were performed, one yielding a minimum of 3.0 eV and the other again our lowest minimum structure. The minimal HOMO-LUMO gap structure (denoted as dia_min_gap) is depicted in Table 4 and has an energy gap of 2.94 eV. This is considerably lower than the 9.00 eV found for the unsubstituted diamantane and closer to the minimal gap of the adamantane derivative (2.42 eV), with the former having only 6 functional sites. In fact, the dia_min_gap structure looks quite similar to the adamantane case, with a nitro group assisted by hydrogen-bonded hydroxyl groups on the one side and electron-rich dopants on the other side. Again, the different optimal structures along the BFS optimization are given in the Supporting Information. During the optimization process, the variation in HOMO-LUMO gap values over all possible functionalizations becomes larger for each of the sites. However, for most sites there exist several functionalization possibilities to obtain a similar result.

Figure 6. Chronological order of the BFS for a minimal HOMO-LUMO gap diamantane derivative; the black lines separate the global iterations. Per separate color the lower limit represents the current structure with the minimal HOMO-LUMO gap. After two global iterations, the first BFS converged to a value of 4.09 eV. Subsequently, a Monte Carlo step was undertaken. The second BFS converged in four global iterations to the final value of 2.94 eV.



12 ACS Paragon Plus Environment

Page 13 of 25

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

Journal of Chemical Theory and Computation

Again, we observe a spatial separation between the frontier molecular orbitals and a similar trend in HOMO and LUMO. The high-lying HOMO level is localized in the cage, mainly on the free electron pairs of the two nitrogen atoms. The presence of the methyl group in between can be explained by the fact that it lifts the HOMO energy level by stabilization of possible positive charge on the nitrogen atoms by hyperconjugation. As in the case of adamantane, a low-lying LUMO is observed, located mainly outside the cage, on the nitro group, assisted through strong intramolecular hydrogen bonding by the two neighboring hydroxyl groups. The partially positively charged hydrogen atoms from the two hydroxyl groups stabilize partial negative charge on the oxygen atoms of the nitro group. The confirming NCI plot is given in the Supporting information. The strength of the hydrogen bonds was determined by the energy difference of the conformations with the hydrogens of the hydroxyl groups pointing towards and away from the nitro group. The individual hydrogen bond strength was calculated to be approximately 5 kcal mol-1. Because the result is so largely influenced by the presence of the nitro group, the final result will depend on which site is changed first in the optimization process. The probability of the first chemical modification being the introduction of a nitro group is rather high, because including this group gives by far the lowest-lying LUMO (cf. Figure 2). Therefore, searches were performed using the three symmetrically distinctive sites (1,3,5) as the starting site. This indeed resulted in three different optima with the nitro group present on the starting site. The structures are depicted in Table 5. The minimal gaps are 3.13, 3.51, and 2.94 eV, respectively; the last one corresponds to dia_min_gap presented above. The resulting minima look very similar. The HOMO part is identical in the three cases. It can be seen that the nitro group always participates in intramolecular hydrogen bonding with its neighbor(s). When the nitro group is put on a tertiary carbon atom, the hydrogen bonding is evidently limited to one neighbor, due to our choice of adjustable sites. This, however, confirms the importance of intramolecular hydrogen bonding in lowering the LUMO energy. Table 5. Three structures resulting from a BFS run to a diamantane derivative with the lowest HOMOLUMO gap with starting site 1, 3 and 5, respectively.

Apical site O

Medial site

Secondary site H

O H

N

N

O

H N

O

H

O N O

O N

N

O N

O

N H3 C

H3 C

H

N

N H3 C

N

N

HOMO-LUMO gap = 3.13 eV HOMO-LUMO gap = 3.51 eV HOMO-LUMO gap = 2.94 eV Since the nitro-group is so significantly determining the results, it is interesting to examine which minima can be found when the nitro group is removed from the functionalization library. The start and converged configuration for each BFS procedure, followed by a Monte Carlo step, are given in the Supporting Information. In total, we have performed six subsequent Monte Carlo steps, in order to escape possible local minima. It turns out that the carbonyl group takes over the role as electron-withdrawing group on which the LUMO mainly resides, in accordance with our findings for the mono-substituted diamantanes. Even though the carbonyl group can only be 13 ACS Paragon Plus Environment

Journal of Chemical Theory and Computation

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 14 of 25

placed at a secondary site, the BFS results still depend a lot on the start structure. Almost every consecutive Monte Carlo step led to a more suitable structure, ultimately lowering the energy gap of the initial minimum by 0.3 eV. This lowest-energy-gap structure was the result of 3 out of 7 Monte Carlo directed BFS searches, indicating that this might be the global minimum in our search space. All minima have the LUMO corresponding to a σ*-orbital located on the carbonyl group. The carbonyl has in all cases a hydroxyl group as a neighbor forming an intramolecular hydrogen bond. The structural part on the other end of the carbonyl group is different each time, containing either a sulfur atom or, one or two nitrogen atoms. Only the lowest minimum structure shows similar ending as dia_min_gap in Table 4.

Inverse molecular design: maximizing the HOMO-LUMO energy gap Although fewer applications are foreseen for large-gap diamondoid derivatives, we are interested in how much the HOMO-LUMO energy gap can be stretched. Based on the monosubstituted diamondoid results in Figure 2, a competition between hydrogen, fluorine, trifluoromethyl, and methyl groups is expected. In addition, we anticipate only a small energy gap increase compared to the pristine diamondoid structures. After five global iterations an adamantane structure is found with a HOMO-LUMO gap of 10.63 eV, an increase of merely 1.18 eV compared to the unsubstituted case. This structure, depicted in Table 4, has C3v symmetry and contains methyl, trifluoromethyl, fluoro groups, and oxygen dopants. For diamantane, also after five global iterations, a structure is found containing only fluoro and fluoromethyl groups, having a HOMO-LUMO gap of 10.08 eV. The two ends of the structure are identical in this case. For both maxima, no new HOMO or LUMO levels are introduced. As illustrated in Figure 7, the increase in energy gap can be attributed to a lowering of the HOMO level. Functionalization with fluorine and trifluoromethyl already lowers the HOMO energy level in the mono-substituted cases. The LUMO is still a surface state of the C-H anti-bonding sigma-bonds. Mono-substitution of a diamondoid site with methyl, alcohol, fluorine, or trifluoromethyl, or doping with oxygen has no significant effect on the LUMO energy level (cf. Figure 2). Hence, the dia_max_gap structure is easy to understand. The ada_max_gap structure contains electron-withdrawing fluorine and trifluoromethyl groups, but also three oxygen atoms. Oxa-adamantane, however, has a HOMO level of -6.3 eV, around 1 eV more than for pure adamantane. So surprisingly, the behavior of the oxygen dopant is changed when positioned between fluorine groups. More detailed information about the maximization of HOMO-LUMO gaps of adamantane and diamantane are gathered in the Supporting Information. The general trend of energy gap vs. substituent is in these cases much more consistent among all sites. Effects of interchanging methyl, fluoromethyl, fluoro, and oxygen are rather small. All give HOMO-LUMO gaps around 10 eV. Because the HOMO and LUMO are now both delocalized in or over the cage and not reside on a single functional group, all sites show similar behavior. Hence, the independent site approximation is largely valid in this subclass of diamondoid chemical space. Note that, like for the mono-substituted diamondoids, carbonyl is predicted to have a larger impact on reducing the HOMO-LUMO energy gap than the nitro group.



14 ACS Paragon Plus Environment

Page 15 of 25

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

Journal of Chemical Theory and Computation

Figure 7. The frontier molecular orbital levels of pure adamantane and diamantane and the minimal and maximal HOMO-LUMO energy gap structures.

The HOMO-LUMO gap: HOMO versus LUMO A HOMO-LUMO energy gap optimization is in essence a linear scalarization of a multipleobjective problem. The actual optimization property is a linear combination of individual optimization parameters. In this case, the HOMO-LUMO energy gap is a linear combination of two optimization parameters, the HOMO and LUMO energy with coefficients -1 and 1, respectively. A more thorough understanding of the HOMO-LUMO energy gap optimization could thus be obtained by understanding how the individual parameters behave under optimization. The HOMO and LUMO energy levels could be subjected to individual maximization and minimization. Because we are mainly interested in minimal HOMO-LUMO gaps, only the HOMO maximization and the LUMO minimization were performed. The results are presented in Table 6. Note that via Koopmans’ theorem the HOMO energy level can be related to the ionization potential (IP). Therefore, HOMO maximization can be conceptually understood as a minimization of the IP. Dopants or functional groups with free electron pairs, such as nitrogen and sulfur, and amino, thiol, and carbonyl groups, having high-lying HOMOs, can introduce a new HOMO level into the HOMO-LUMO gap of the pure diamondoids and as such assist in maximizing the HOMO. Accordingly, the resulting structures of the HOMO maximization contain mainly amino groups and nitrogen substituents. It can be observed, however, that in dia_max_homo, possessing Ci symmetry, the lone pairs of the amino groups do not contribute to the HOMO, although they might influence the energy level. On the other hand, functional groups with free electron pairs, such as nitro, carbonyl, and carboxyl groups, having low-lying LUMOs, can introduce a new LUMO level into the energy gap of the pure diamondoids, which helps to minimize the LUMO. The structures of the LUMO minimization indeed contain many nitro and carbonyl groups. In ada_min_lumo, with Cs symmetry, the LUMO

15 ACS Paragon Plus Environment

Journal of Chemical Theory and Computation

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 16 of 25

is localized both over the entire cage and on the substituents’ lone pairs. This is not the case for dia_min_lumo, for which the LUMO resides outside the cage, on three out of four nitro groups. The cyano groups help to stabilize the LUMO, but not through their lone pairs. Also note that the LUMO of ada_min_lumo is even lower than the HOMO of ada_max_homo. When we compare these orbital energy values with the corresponding values in the minimal HOMO-LUMO gap structures, we observe a further increase in the HOMO level of around 0.80.9 eV for ada_max_homo and dia_max_homo, and a lowering in the LUMO level of approximately 0.5 eV for dia_min_lumo and a significant 1.8 (!) eV for ada_min_lumo. Nevertheless, the same functionalization patterns appear in both sets of structures. This confirms that, instead of selecting substituents or dopants that perform well on average by both lowering the LUMO and increasing the HOMO via stabilization/destabilization effects or the introduction of new relevant orbital levels, push-pull doping is the most effective. Incorporating electronwithdrawing functional groups on the one end and electron-donating groups on the other end results in the largest energy gap reduction, spatially separating the two frontier molecular orbitals. Table 6. Diamondoid derivatives found by searching for maximizations of the HOMO energy level of adamantane and diamantane and for minimizations of the LUMO energy level of adamantane and diamantane

ada_max_homo

dia_max_homo

NH 2

ada_min_lumo O

NH 2

N

N

N

O

O

H2 N H2 N

N

N NH H2 N

-4.68 eV

O

O

O

N

N

N

O

N

N

O O

O

O

N O

-4.65 eV

O O

O

H2 N

N

O

O

N

N

dia_min_lumo

O

O O

-4.90 eV

N N O

O

-3.16 eV

The optimal structures from the diamantane derivatives were rather symmetrical and intuitive, while the ada_min_gap structure is much more complex. We therefore propose some adamantane derivatives with fewer substituted sites that still have significantly lowered HOMOLUMO gap. These structures should all have some push-pull character. The push component is formed primarily by sulfur (on secondary site) or nitrogen dopants. Nitrogen dopants give higher HOMO levels when placed on tertiary rather than on secondary sites. Multiple nitrogen dopants can be included on different tertiary sites that are next nearest neighbors, as is the case in the dia_min_gap structure. Furthermore, the pulling component should consist of an electronwithdrawing group, a nitro or a carbonyl group. Placing neighboring electronegative groups lowers the LUMO energy levels even further, as is done here with two hydroxyl groups that form hydrogen bonds with the nitro group. Doubly substituted adamantane with a nitro-group and a 16 ACS Paragon Plus Environment

Page 17 of 25

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

Journal of Chemical Theory and Computation

secondary sulfur or tertiary nitrogen dopant lowers the gap to 4.2 eV or 4.1 eV, respectively (Table 7). Triply and quaternary substitution of adamantane by adding hydroxyl groups besides the nitro group lowers the gap to 3.9 and 3.5 eV, respectively. Subsequent addition of two amine groups on the next nearest neighbors of the sulfur atom results in a symmetrical structure with a HOMO-LUMO energy gap of 3.0 eV. Interestingly, when adamantane is substituted in a similar pattern as the dia_min_gap structure, the HOMO-LUMO gap is also 3.0 eV. Table 7. Two proposed pathways with increasing number of functionalizations leading to structures with small HOMO-LUMO energy gaps, namely ada_min_gap_red1 and ada_min_gap_red2. The upper part follows a path towards structures similar to ada_min_gap while the lower part follows the path of increasing functionalizations towards the same substitution pattern of the dia_min_gap structure. OH NO 2

OH

NO 2

H2 N

NO 2

NO 2

OH

OH

NO 2 S

S

4.23 eV

S

OH S

3.91 eV

NH 2

3.54 eV

3.02 eV OH

OH

9.45 eV

NO 2

NO 2

6.14 eV

NO 2

NO 2

N

N

N

N

4.00 eV

3.88 eV

N

3.05 eV

N

OH H3 C

N

OH

2.99 eV

The suggested structures were based on ada_min_gap1 (and earlier minima in that BFS) and dia_min_gap. We also considered reducing the number of substituents of ada_min_gap2, the structure with the lowest energy gap. However, it turned out that nearly all functionalizations in that structure have a non-negligible impact on the energy gap reduction. To support this idea, we built up the final structure, starting from pure adamantane and adding step by step the functionalization having the largest impact at this particular step. Since a carbonyl group in a mono-substituted adamantane results in a larger energy gap decrease than a nitro group, two different pathways were considered: one starting with the introduction of a carbonyl group and one introducing firstly a nitro group. We will briefly discuss the NO2-pathway, as shown in Figure 8. You can consult the Supporting Information for more details. The largest energy gap reduction happens in the first two steps, with NO2 causing a significant lowering of the LUMO level and S causing a large increase in the HOMO level, inducing already push-pull character. The next three steps further decrease the LUMO level, especially via the inclusion of hydrogenbonded hydroxyl groups. The functionalizations added in steps 6 to 10 are then used to further push up the HOMO level. In this pathway, it is clear that lowering the LUMO energy level is more beneficial than increasing the HOMO level. However, in order to obtain a similar energy gap value as the proposed reduced structures in Table 7, we would have to include 8 of the 10 functionalizations. This clearly shows the synergetic effect of substituents and dopants in ada_min_gap2 to attain a large HOMO-LUMO gap reduction. What is very interesting to see is that throughout the build-up, the HOMO-LUMO spatial separation is manifesting more distinctly, with the LUMO residing on the nitro-group and the HOMO shifting from the sulfur dopant, over CO and S, then to NH2 on site 1, to finally being located on SH. The structure is thus being divided in a part that helps to lower the LUMO and a part that aids in increasing the HOMO.

17 ACS Paragon Plus Environment

Journal of Chemical Theory and Computation

4

location of HOMO

6.16

HOMO energy

NO 2

NH 2

H 2N

4.50

4.02

3.83

3.70

3.33

3.07

2.68

2.42

CNHH (1)

CSH (6)

CNHH (2)

CNHH (7)

CNHH (3)

NH 2

CO (8)

COH (10)

S (9)

CNOO (4)

-1

4.21

O

COH (5)

H 2N SH

1

energy gap

S

HO

2

0

LUMO energy

HO

3

reduction [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

Page 18 of 25

-2 Figure 8. The reduction in HOMO and LUMO energy level, and the resulting energy gap, when adding the functionalization with the largest impact (the site number is included in parenthesis). The numbers above the bars indicate the energy gap value at that stage. In addition, the increasing HOMO-LUMO spatial separation is shown by the dashed lines.

Even lower HOMO-LUMO gaps might be obtained by using e.g. charged functional groups, diamondoids with more cages like triamantane, tetramantane etc., or by using other functionalizations such as boron or phosphor dopants, or large conjugated functional groups as phenyl and sulfonyl. To conclude, we would like to shed some light on how the HOMO-LUMO gap using the Kohn-Sham (KS) frontier orbitals from DFT can be interpreted. For molecular systems, the KS gap in fact becomes a good approximation to the molecular optical gap.60 In addition, we have performed extra calculations to more accurately describe the molecular “band gap” for a selection of diamondoids. In bulk materials, the band gap is a well-defined quantity. The equivalent in molecular systems is the so-called fundamental gap, being the energy difference between the vertical electron affinity (EA) and the vertical ionization potential (IP). It is important to note, however, that the band gap for solids is typically considerably smaller in energy than the fundamental gap in molecules.61 There is thus a clear distinction between the optical gap and the fundamental gap with the difference being a measure of the electron–hole pair binding energy.61 For a detailed discussion on the difference between HOMO-LUMO energy gap and IP-EA in DFT, we refer to Baerends et al.,60 Tozer et al.,62 and Teale et al. 63. It is therefore no surprise that there exist large differences between the two sets of quantities. In Table 8, we have collected the IP-EA results for pure adamantane and diamantane, the reduced structures ada_min_gap_red1 and ada_min_gap_red2 (Table 7), and dia_min_gap using the following levels of theory: B3PW91/6-31G(d,p), MP2/cc-pVDZ, MP4(SDQ)/cc-pVDZ, and CCSD(T)/cc-pVDZ. The results for other computational methods are reported in the Supporting Information. Moreover, to obtain more accurate EA and IP values, extrapolations were performed toward the CCSD(T) level of theory in the limit of an asymptotically complete basis set 18 ACS Paragon Plus Environment

Page 19 of 25

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

Journal of Chemical Theory and Computation

(CCSD(T)/cc-pV∞Z) using the principles of a focal point analysis (FPA)64-71. In such an approach, the faster convergence of the higher order correlation corrections to the calculated energy differences is exploited in well-suited extrapolations of results obtained using CCSD(T) theory. More specifically, the so-called FPA-QZ method 72-75 was used, already proven to be able to reproduce IPs and positive EAs much better than 0.1 eV, and negative EAs better than 0.2 eV accuracy. Remark that for the structures with a positive EA, namely ada_min_gap_red1, ada_min_gap_red2, and dia_min_gap, the DFT functional B3PW91 with basis set 6-31G(d,p) is deviating less than 0.3 eV from our benchmark data. Table 8. Computed IP-EA values (in eV) with IP the vertical ionization potential and EA the vertical electron affinity, using different levels of theory. CCSD(T)/cc-pV∞Z using the FPA-QZ method serves as benchmark.

B3PW91 MP2 MP4(SDQ) CCSD(T) FPA-QZ 6-31(d,p) cc-pVDZ cc-pVDZ cc-pVDZ adamantane 12.77 12.57 13.36 13.02 11.43 diamantane 12.27 12.14 12.92 12.30 10.78 ada_min_gap_red1 7.39 8.34 8.19 8.13 7.64 ada_min_gap_red2 7.44 8.61 8.17 8.25 7.72 dia_min_gap 7.22 8.15 7.96 7.93 7.42 structure

Conclusions We used an Inverse Molecular Design approach to search for diamondoid derivatives with minimal and maximal HOMO-LUMO gaps. Using Monte Carlo steps within the design algorithm effectively helped to escape local optima towards more optimal structures and also decreased the number of calculations that are far from optimal. These combined techniques resulted in the discovery of several adamantane and diamantane derivatives with HOMO-LUMO gaps even smaller than 3.0 eV, and without incorporating charged functional groups. The minimal gap structures clearly show push-pull character with the HOMO and LUMO located on respectively electron-donating (n-doping) and electron-withdrawing (p-doping) groups. The structural element showing the lowest-lying LUMO was identified to be an electron withdrawing nitro group that is assisted by neighboring hydroxyl or amino groups forming intramolecular hydrogen bonds. A high-lying HOMO was mostly obtained by electron-donating dopants (S, N) in the diamondoid cages. Sulfur and nitrogen dopants were assigned to respectively secondary and tertiary sites. Groups that stabilize positive charge were placed close to these dopants, mostly stabilized by hyperconjugation. In all cases, we observed a spatial separation of the frontier molecular orbitals. Maximization of the HOMO-LUMO gap led to diamondoid derivatives with only electronegative groups as trifluoromethyl and fluorine groups, sometimes accompanied by oxygen dopants. The HOMO and LUMO were, as in pure diamondoids, located on or over the whole cage of the molecule and only the HOMO energy level was decreased. No push-pull doping was present. Optionally, larger diamondoids or diamondoids with charged groups could bring down the HOMO-LUMO gap even further, but similar trends and cooperation of functional groups are expected.



19 ACS Paragon Plus Environment

Journal of Chemical Theory and Computation

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 20 of 25

Acknowledgments F.D.V. acknowledges the Research Foundation-Flanders (FWO) for a postdoctoral fellowship (1227014N). F.D.P. acknowledges the Fund for Scientific Research-Flanders (FWO) and the Free University of Brussels (VUB) for continuous support to his research group, in particular the VUB for awarding a Strategic Research Program (SRP) to the ALGC research group started on January 1, 2013. J.L.T also wishes to acknowledge the Free University of Brussels (VUB) for support via the SRP.

Supporting Information The Supporting Information contains computed adiabatic ionization potential values and their comparison with experimental values, the HOMO-LUMO gap data for the mono-substituted diamondoids, the optimization pathways for all BFS runs, plots of the BFS convergence and the substituent vs. energy gap correlation for each site, plotted HOMO and LUMO of the first adamantane energy gap minimization, a table showing the Monte Carlo steps of the diamantane gap minimization without NO2 in the functionalization library, Noncovalent Interaction index plots for all minima, data on the build-up of structure ada_min_gap2, and energy data used for the molecular energy gap calculations.

References (1) Vörös, M.; Demjén, T.; Szilvási, T.; Gali, A. Tuning the Optical Gap of Nanometer-Size Diamond Cages by Sulfurization: A Time-Dependent Density Functional Study. Phys. Rev. Lett. 2012, 108, 267401. (2) McIntosh, G. C.; Yoon, M.; Berber, S.; Tománek, D. Diamond Fragments as Building Blocks of Functional Nanostructures. Phys. Rev. B 2004, 70, 45401. (3) Ashfold, M. N. R.; May, P. W.; Rego, C. A.; Everitt, N. M. Thin Film Diamond by Chemical Vapour Deposition Methods. Chem Soc Rev 1994, 23, 21–30. (4) Dahl, J. E.; Liu, S. G.; Carlson, R. M. K. Isolation and Structure of Higher Diamondoids, Nanometer-Sized Diamond Molecules. Science 2003, 299, 96–99. (5) Gonçalves, M. S. T. Fluorescent Labeling of Biomolecules with Organic Probes. Chem. Rev. 2009, 109, 190–212. (6) Wang, Y.; Kioupakis, E.; Lu, X.; Wegner, D.; Yamachika, R.; Dahl, J. E.; Carlson, R. M. K.; Louie, S. G.; Crommie, M. F. Spatially Resolved Electronic and Vibronic Properties of Single Diamondoid Molecules. Nat. Mater. 2008, 7, 38–42. (7) Sakai, Y.; Nguyen, G. D.; Capaz, R. B.; Coh, S.; Pechenezhskiy, I. V.; Hong, X.; Wang, F.; Crommie, M. F.; Saito, S.; Louie, S. G.; Cohen, M. L. Intermolecular Interactions and Substrate Effects for an Adamantane Monolayer on the Au(111) Surface. Phys. Rev. B 2013, 88, 235407. (8) Xue, Y.; Mansoori, G. A. Self-Assembly of Diamondoid Molecules and Derivatives (MD Simulations and DFT Calculations). Int. J. Mol. Sci. 2010, 11, 288–303. (9) Clay, W. A.; Dahl, J. E. P.; Carlson, R. M. K.; Melosh, N. A.; Shen, Z.-X. Physical Properties of Materials Derived from Diamondoid Molecules. Rep. Prog. Phys. 2015, 78, 16501. (10) Clay, W. A.; Liu, Z.; Yang, W.; Fabbri, J. D.; Dahl, J. E.; Carlson, R. M. K.; Sun, Y.; Schreiner, P. R.; Fokin, A. A.; Tkachenko, B. A.; Fokina, N. A.; Pianetta, P. A.; Melosh, N.; Shen, Z. X.

20 ACS Paragon Plus Environment

Page 21 of 25

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

Journal of Chemical Theory and Computation

Origin of the Monochromatic Photoemission Peak in Diamondoid Monolayers. Nano Lett. 2009, 9, 57–61. (11) Xue, Y.; Mansoori, G. A. Quantum Conductance and Electronic Properties of Lower Diamondoid Molecules and Derivatives. Int. J. Nanosci. 2008, 7, 63–72. (12) Rander, T.; Staiger, M.; Richter, R.; Zimmermann, T.; Landt, L.; Wolter, D.; Dahl, J. E.; Carlson, R. M. K.; Tkachenko, B. A.; Fokina, N. A.; Schreiner, P. R.; Möller, T.; Bostedt, C. Electronic Structure Tuning of Diamondoids through Functionalization. J. Chem. Phys. 2013, 138, 024310. (13) Landt, L.; Staiger, M.; Wolter, D.; Klünder, K.; Zimmermann, P.; Willey, T. M.; van Buuren, T.; Brehmer, D.; Schreiner, P. R.; Tkachenko, B. A.; Fokin, A. A.; Möller, T.; Bostedt, C. The Influence of a Single Thiol Group on the Electronic and Optical Properties of the Smallest Diamondoid Adamantane. J. Chem. Phys. 2010, 132, 024710. (14) Adhikari, B.; Fyta, M. Towards Double-Functionalized Small Diamondoids: Selective Electronic Band-Gap Tuning. Nanotechnology 2015, 26, 35701. (15) Fokin, A. A.; Schreiner, P. R. Band Gap Tuning in Nanodiamonds: First Principle Computational Studies. Mol. Phys. 2009, 107, 823–830. (16) Ertl, P. Cheminformatics Analysis of Organic Substituents: Identification of the Most Common Substituents, Calculation of Substituent Properties, and Automatic Identification of Drug-like Bioisosteric Groups. J. Chem. Inf. Comput. Sci. 2003, 43, 374– 380. (17) Lipinski, C.; Hopkins, A. Navigating Chemical Space for Biology and Medicine. Nature 2004, 432, 855–861. (18) Dobson, C. M. Chemical Space and Biology. Nature 2004, 432, 824–828. (19) Subramaniam, S.; Mehrotra, M.; Gupta, D. Virtual High Throughput Screening (vHTS)-a Perspective. Bioinformation 2008, 3, 14–17. (20) Curtarolo, S.; Hart, G. L. W.; Nardelli, M. B.; Mingo, N.; Sanvito, S.; Levy, O. The HighThroughput Highway to Computational Materials Design. Nat. Mater. 2013, 12, 191– 201. (21) Pyzer-Knapp, E. O.; Suh, C.; Gómez-Bombarelli, R.; Aguilera-Iparraguirre, J.; AspuruGuzik, A. What Is High-Throughput Virtual Screening? A Perspective from Organic Materials Discovery. Annu. Rev. Mater. Sci. 2015, 45, 195-216. (22) Er, S.; Suh, C.; P. Marshak, M.; Aspuru-Guzik, A. Computational Design of Molecules for an All-Quinone Redox Flow Battery. Chem. Sci. 2015, 6, 885–893. (23) Gómez-Bombarelli, R.; Aguilera-Iparraguirre, J.; Hirzel, T. D.; Duvenaud, D.; Maclaurin, D.; Blood-Forsythe, M. A.; Chae, H. S.; Einzinger, M.; Ha, D.-G.; Wu, T.; Markopoulos, G.; Jeon, S.; Kang, H.; Miyazaki, H.; Numata, M.; Kim, S.; Huang, W.; Ik Hong, S.; Baldo, M.; Adams, R. P.; Aspuru-Guzik, A. Design of Efficient Molecular Organic Light-Emitting Diodes by a High-Throughput Virtual Screening and Experimental Approach. Nat. Mater. 2016, 15, 1120–1127. (24) Kuhn, C.; Beratan, D. N. Inverse Strategies for Molecular Design. J. Phys. Chem. 1996, 100, 10595–10599. (25) Franceschetti, A.; Zunger, A. The Inverse Band-Structure Problem of Finding an Atomic Configuration with given Electronic Properties. Nature 1999, 402, 60–63. (26) Jóhannesson, G. H.; Bligaard, T.; Ruban, A. V.; Skriver, H. L.; Jacobsen, K. W.; Nørskov, J. K. Combined Electronic Structure and Evolutionary Search Approach to Materials Design. Phys. Rev. Lett. 2002, 88, 255506.

21 ACS Paragon Plus Environment

Journal of Chemical Theory and Computation

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 25

(27) von Lilienfeld, O. A.; Lins, R. D.; Rothlisberger, U. Variational Particle Number Approach for Rational Compound Design. Phys. Rev. Lett. 2005, 95, 153002. (28) Kim, K.; Graf, P. A.; Jones, W. B. A Genetic Algorithm Based Inverse Band Structure Method for Semiconductor Alloys. J. Comput. Phys. 2005, 208, 735–760. (29) Wang, M.; Hu, X.; Beratan, D. N.; Yang, W. Designing Molecules by Optimizing Potentials. J. Am. Chem. Soc. 2006, 128, 3228–3232. (30) Balamurugan, D.; Yang, W.; Beratan, D. N. Exploring Chemical Space with Discrete, Gradient, and Hybrid Optimization Methods. J. Chem. Phys. 2008, 129, 174105. (31) Kanal, I. Y.; Owens, S. G.; Bechtel, J. S.; Hutchison, G. R. Efficient Computational Screening of Organic Polymer Photovoltaics. J. Phys. Chem. Lett. 2013, 4, 1613–1623. (32) Ng, L. Y.; Chong, F. K.; Chemmangattuvalappil, N. G. Challenges and Opportunities in Computer-Aided Molecular Design. Comput. Chem. Eng. 2015, 81, 115–129. (33) Pearl; J.; Korf, R. E. Search Techniques. Annu. Rev. Comput. Sci. 1987, 2, 451–467. (34) De Vleeschouwer, F.; Yang, W.; Beratan, D. N.; Geerlings, P.; De Proft, F. Inverse Design of Molecules with Optimal Reactivity Properties: Acidity of 2-Naphthol Derivatives. Phys. Chem. Chem. Phys. 2012, 14, 16002. (35) De Vleeschouwer, F.; Geerlings, P.; De Proft, F. Molecular Property Optimizations with Boundary Conditions through the Best First Search Scheme. Chemphyschem 2016, 17, 1414–1424. (36) Elward, J. M.; Rinderspacher, B. C. Smooth Heuristic Optimization on a Complex Chemical Subspace. Phys. Chem. Chem. Phys. 2015, 17, 24322–24335. (37) Xiao, D.; Yang, W.; Beratan, D. N. Inverse Molecular Design in a Tight-Binding Framework. J. Chem. Phys. 2008, 129, 44106. (38) Hu, X.; Beratan, D. N.; Yang, W. A Gradient-Directed Monte Carlo Approach to Molecular Design. J. Chem. Phys. 2008, 129, 64102. (39) Hu, X.; Beratan, D. N.; Yang, W. A Gradient-Directed Monte Carlo Method for Global Optimization in a Discrete Space: Application to Protein Sequence Design and Folding. J. Chem. Phys. 2009, 131, 154117. (40) Frisch, M. J.; Trucks, G. W.; Schlegel, H. B.; Scuseria, G. E.; Robb, M. A.; Cheeseman, J. R.; Scalmani, G.; Barone, V.; Petersson, G. A.; Nakatsuji, H.; Li, X.; Caricato, M.; Marenich, A.; Bloino, J.; Janesko, B. G.; Gomperts, R.; Mennucci, B.; Hratchian, H. P.; Ortiz, J. V.; Izmaylov, A. F.; Sonnenberg, J. L.; Williams-Young, D.; Ding, F.; Lipparini, F.; Egidi, F.; Goings, J.; Peng, B.; Petrone, A.; Henderson, T.; Ranasinghe, D.; Zakrzewski, V. G.; Gao, J.; Rega, N.; Zheng, G.; Liang, W.; Hada, M.; Ehara, M.; Toyota, K.; Fukuda, R.; Hasegawa, J.; Ishida, M.; Nakajima, T.; Honda, Y.; Kitao, O.; Nakai, H.; Vreven, T.; Throssell, K.; Montgomery, Jr., J. A.; Peralta, J. E.; Ogliaro, F.; Bearpark, M.; Heyd, J. J.; Brothers, E.; Kudin, K. N.; Staroverov, V. N.; Keith, T.; Kobayashi, R.; Normand, J.; Raghavachari, K.; Rendell, A.; Burant, J. C.; Iyengar, S. S.; Tomasi, J.; Cossi, M.; Millam, J. M.; Klene, M.; Adamo, C.; Cammi, R.; Ochterski, J. W.; Martin, R. L.; Morokuma, K.; Farkas, O.; Foresman, J. B. and Fox, D. J. Gaussian 09, Revision ; Gaussian, Inc.: Wallingford, CT, USA, 2009. (41) Perdew, J. P.; Wang, Y. Accurate and Simple Analytic Representation of the ElectronGas Correlation Energy. Phys. Rev. B 1992, 45, 13244–13249. (42) Hariharan, P. C.; Pople, J. A. The Influence of Polarization Functions on Molecular Orbital Hydrogenation Energies. Theor. Chim. Acta 1973, 28, 213–222.



22 ACS Paragon Plus Environment

Page 23 of 25

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

Journal of Chemical Theory and Computation

(43) Schreiner, P. R.; Fokin, A. A.; Pascal, Robert A.; de Meijere, A. Many Density Functional Theory Approaches Fail To Give Reliable Large Hydrocarbon Isomer Energy Differences. Org. Lett. 2006, 8, 3635–3638. (44) Zhang, C.; Hu, J.; Chen, S.; Ji, F. Theoretical Study of Hydrogen Bonding Interactions on MDI-Based Polyurethane. J. Mol. Model. 2010, 16, 1391–1399. (45) Chandra, A. K.; Nguyen, M. T. Quantum Chemical Study of the Hydrogen-Bonded C4H2···HCl Complex. J. Chem. Res. (S) 1997, 216–217. (46) Staroverov, V. N.; Scuseria, G. E.; Tao, J.; Perdew, J. P. Comparative Assessment of a New Nonempirical Density Functional: Molecules and Hydrogen-Bonded Complexes. J. Chem. Phys. 2003, 119, 12129–12137. (47) Xiao, H.; Tahir-Kheli, J.; Goddard, W. A. Accurate Band Gaps for Semiconductors from Density Functional Theory. J. Phys. Chem. Lett. 2011, 2, 212–217. (48) Yang, W. L.; Fabbri, J. D.; Willey, T. M.; Lee, J. R. I.; Dahl, J. E.; Carlson, R. M. K.; Schreiner, P. R.; Fokin, A. A.; Tkachenko, B. A.; Fokina, N. A.; Meevasana, W.; Mannella, N.; Tanaka, K.; Zhou, X. J.; van Buuren, T.; Kelly, M. A.; Hussain, Z.; Melosh, N. A.; Shen, Z.-X. Monochromatic Electron Photoemission from Diamondoid Monolayers. Science 2007, 316, 1460–1462. (49) Drummond, N. D.; Williamson, A. J.; Needs, R. J.; Galli, G. Electron Emission from Diamondoids: A Diffusion Quantum Monte Carlo Study. Phys. Rev. Lett. 2005, 95, 96801. (50) Lu, A. J.; Pan, B. C.; Han, J. G. Electronic and Vibrational Properties of Diamondlike Hydrocarbons. Phys. Rev. B 2005, 72, 35447. (51) Heinrich, N.; Koch, W.; Frenking, G. On the Use of Koopmans’ Theorem to Estimate Negative Electron Affinities. Chem. Phys. Lett. 1986, 124, 20–25. (52) Guerra, M. On the Use of Diffuse Functions for Estimating Negative Electron Affinities with LCAO Methods. Chem. Phys. Lett. 1990, 167, 315–319. (53) Lenzke, K.; Landt, L.; Hoener, M.; Thomas, H.; Dahl, J. E.; Liu, S. G.; Carlson, R. M. K.; Möller, T.; Bostedt, C. Experimental Determination of the Ionization Potentials of the First Five Members of the Nanodiamond Series. J. Chem. Phys. 2007, 127, 84320. (54) Fokin, A. A.; Tkachenko, B. A.; Gunchenko, P. A.; Gusev, D. V.; Schreiner, P. R. Functionalized Nanodiamonds Part I. An Experimental Assessment of Diamantane and Computational Predictions for Higher Diamondoids. Chem. Eur. J. 2005, 11, 7091– 7101. (55) Fokina, N. A.; Tkachenko, B. A.; Merz, A.; Serafin, M.; Dahl, J. E. P.; Carlson, R. M. K.; Fokin, A. A.; Schreiner, P. R. Hydroxy Derivatives of Diamantane, Triamantane, and [121]Tetramantane: Selective Preparation of Bis-Apical Derivatives. Eur. J. Org. Chem. 2007, 28, 4738–4745. (56) Schreiner, P. R.; Fokina, N. A.; Tkachenko, B. A.; Hausmann, H.; Serafin, M.; Dahl, J. E. P.; Liu, S.; Carlson, R. M. K.; Fokin, A. A. Functionalized Nanodiamonds: Triamantane and [121]Tetramantane. J. Org. Chem. 2006, 71, 6709–6720. (57) Demján, T.; Vörös, M.; Palummo, M.; Gali, A. Electronic and Optical Properties of Pure and Modified Diamondoids Studied by Many-Body Perturbation Theory and TimeDependent Density Functional Theory. J. Chem. Phys. 2014, 141, 64308. (58) Johnson, E. R.; Keinan, S.; Mori-Sánchez, P.; Contreras-García, J.; Cohen, A. J.; Yang, W. Revealing Noncovalent Interactions. J. Am. Chem. Soc. 2010, 132, 6498–6506.



23 ACS Paragon Plus Environment

Journal of Chemical Theory and Computation

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 24 of 25

(59) Contreras-García, J.; Johnson, E. R.; Keinan, S.; Chaudret, R.; Piquemal, J.-P.; Beratan, D. N.; Yang, W. NCIPLOT: A Program for Plotting Noncovalent Interaction Regions. J. Chem. Theory Comput. 2011, 7, 625–632. (60) Baerends, E. J.; Gritsenko, O. V.; van Meer, R. The Kohn-Sham Gap, the Fundamental Gap and the Optical Gap: the Physical Meaning of Occupied and Virtual Kohn-Sham Orbital Energies. Phys. Chem. Chem. Phys. 2013, 15, 16408-16425. (61) Bredas, J.-L. Mind the Gap! Mater. Horiz. 2014, 1, 17–19. (62) Tozer, D. J. and De Proft, F. Computation of the Hardness and the Problem of Negative Electron Affinities in Density Functional Theory J. Phys. Chem. A 2005, 109, 89238929. (63) Teale, A. M.; De Proft, F.; Tozer, D. J. Orbital Energies and Negative Electron Affinities from Density Functional Theory: Insight from the Integer Discontinuity. J. Chem. Phys. 2008, 129, 044110. (64) Allinger, N. L.; Fermann, J. T.; Allen, W. D. and Schaefer III, H. F. The Torsional Conformations of Butane: Definitive Energetics from Ab Initio Methods. J. Chem. Phys. 1997, 106, 5143-5150. (65) Salam, A. and Deleuze, M. S. High-level theoretical study of the conformational equilibrium of n-pentane. J. Chem. Phys. 2002, 116, 1296-1302. (66) Kwasniewski, S. P.; Claes, L.; François, J.-P. and Deleuze, M. S. High level theoretical study of the structure and rotational barriers of trans-stilbene. J. Chem. Phys. 2003, 118, 7823-7836. (67) East, A. L. L. and Allen, W. D. The Heat of Formation of NCO. J. Chem. Phys. 1993, 99, 4638-4650. (68) Császár, A. G.; Allen, W. D. and Schaefer III, H. F. In Pursuit of the Ab Initio Limit for Conformational Energy Prototypes. J. Chem. Phys. 1998, 108, 9751-9764. (69) Huang, Y. R.; Knippenberg, S.; Hajgató, B.; François, J.-P.; Deng, J. K.; Deleuze M. S. Imaging momentum orbital densities of conformationally versatile molecules: A benchmark theoretical study of the molecular and electronic structures of dimethoxymethane. J. Phys. Chem. A 2007, 111, 5879-5897. (70) Knippenberg, S.; Huang, Y. R.; Hajgató, B.; François, J-P.; Deng, J. K.; Deleuze M. S. Probing molecular conformations in momentum space: The case of n-pentane. J. Chem. Phys. 2007, 127, 174306. (71) Monten, R.; Hajgató, B.; Deleuze, M. S. Many-body calculations of molecular electric polarizabilities in asymptotically complete basis sets. Mol. Phys. 2011, 109, 23172339. (72) Hajgató, B.; Deleuze, M.S.; Tozer, D.J.; De Proft, F. A benchmark theoretical study of the electron affinities of benzene and linear acenes. J. Chem. Phys. 2008, 129, 084308. (73) Hajgató, B.; Szieberth, D.; Geerlings, P.; De Proft, F.; Deleuze, M.S. A benchmark theoretical study of the electronic ground state and of the singlet-triplet split of benzene and linear acenes. J. Chem. Phys. 2009, 131, 224321. (74) B. Hajgató, B.; Huzak, M.; Deleuze, M.S. Focal Point Analysis of the Singlet-Triplet Energy Gap of Octacene and Larger Acenes. J. Phys. Chem. A. 2011, 115, 9282-9293. (75) Huzak, M.; Hajgató, B.; Deleuze, M. S. Benchmark theoretical study of the ionization energies, electron affinities and singlet-triplet energy gaps of azulene, phenanthrene, pyrene, chrysene and perylene. Chem. Phys. 2012, 406, 55-64.



24 ACS Paragon Plus Environment

Page 25 of 25

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

Journal of Chemical Theory and Computation

Table of Contents / Abstract graphic: CF3

F CF3

F3C F CF3

MAX O

H

PURE

O N+ O

N

O– H

N

MIN



25 ACS Paragon Plus Environment