Full Paper
allow the desired lateral polymerization leading to 2DP, we
noted that for a structurally related monomer it was recently
shown in a rather comprehensive crystallization study[7b] that
etf can be rendered into ftf, if the appropriate solvent and the
optimum crystallization conditions are chosen. While this
paper was under review, crystallization of compound 1 into
layered ftf-stacked crystals was also successful.[15]
neglect attractive long-range dispersion effects. They can be
efficiently included through a dispersion correction for atom
pairs at sufficiently large distance. Because of the interactions
between the three anthracene moieties, we expected that dis-
persion interactions play an important role and thus employed
Grimme’s dispersion corrections D and D3.[16] Detailed results
obtained for different computational set-ups are reported in
the Supporting Information. We found that 1B is energetically
favored over 1A by about ꢀ3 to ꢀ4 kcalmolꢀ1, if we take dis-
persion interaction energies into account. By contrast, without
dispersion interactions, 1B is disfavored over 1A by around
1 kcalmolꢀ1. Attractive dispersion interactions also affect the
molecular structures. Thus, in 1B, the dispersion-affected ftf-
stacking blades are by 1.7 ꢂ closer together than in the struc-
ture for which dispersion was not taken into account in the
structure optimization. Comparing dispersion corrected, as well
as dispersion-free results, structures 1A and B are more similar
in energy (DEel(1A!1B) ranges from ꢀ10 to 1 kcalmolꢀ1) than
As was stated before, a computational analysis was initiated
in the beginning of our synthetic efforts to predict the confor-
mational stability of double-decker 1. The intrinsic bent shape
of its oxygen bond can, in principle, lead to various conformers
of 1. We focused on three (Figure 3a): the C3h symmetric rotor-
structures 1A and
C
(DEel(1A!1C) ranges from 14 to
19 kcalmolꢀ1). Because structure 1C has a significantly higher
electronic energy than structure 1A, it is not likely that 1C will
be formed from 1A.
To draw conclusions about transition probabilities from one
minimum to the other, it is necessary to analyze transition bar-
riers between the minima. The calculations revealed that for
the transformation from structure 1A to 1B, at least
19 kcalmolꢀ1 have to be provided to reach the first transition
state, TS-1, and after a shallow energy minimum (1AB),
a second transition state, TS-2, of about the same electronic
energy as TS-1, was found, which connects 1AB with structure
1B. In Figure 3b, we show the rearrangement energy path cal-
culated with BP86, BP86-D3 and TPSS-D3.[17] The electronic-
energy profiles are very similar, and the addition of the D3 dis-
persion correction to the BP86 functional does not change the
electronic energy profile significantly. Due to the double-
decker nature of 1, it is reasonable that the transition from 1A
to 1B occurs in two steps and requires two transition states.
This implies that the C-O-C angles change one after the other
with an energetically shallow intermediate, in which the C-O-C
angles are of the same size, but the upper and the lower part
of the blades are twisted in different directions (1AB). The elec-
tronic-energy difference between the minimum 1AB and the
two transition states is very small (about 1 kcalmolꢀ1). By care-
fully probing our results on technical artifacts (such as optimi-
zation thresholds or numerical grid sizes), we confirmed the
existence of a small dip in the potential energy path from 1A
to 1B (see the Supporting Information).
Figure 3. Results of quantum chemical calculations. a) Quantum chemically
optimized conformations of compound 1. b) Electronic-energy transition
path from structure 1A to 1B calculated with density functionals BP86,
BP86-D3, and TPSS-D3 (all structures were fully optimized with a given func-
tional; see the Supporting Information) and with the def2-TZVP basis set.
The BP86/def2-TZVP-optimized structures are depicted. For TS-1 and TS-2,
the mode with imaginary frequency is shown. Color code: C gray; H white;
N blue; O red.
shaped form (1A), a non-C3 symmetric form (1B), which is ac-
cessible from 1A through a flipping motion of one of the
three blades, hereby twisting the two overlaying oxygen link-
ers to the opposite side, and another C3 symmetric screw-
shaped form (1C), in which the three oxygen bonds connect-
ing one triazine core face in opposite direction compared to
the other ones. To investigate the stability of these conformers,
as well as potential rearrangement pathways, density function-
al theory (DFT) calculations were performed. A detailed assess-
ment of our computation methodology is provided in the Sup-
porting Information.
The calculated energy barrier theoretically permits transfor-
mations from isolated 1A to isolated 1B at elevated tempera-
tures. However, under experimental conditions, also solvent–
solute interactions have to be taken into account. The solubility
of 1 was not known in advance. DMSO is known for its general
strong solvation and was therefore chosen as an example sol-
vent for BOMD simulations of explicitly solvated structures 1A
and 1B to predict whether these interactions would play a role.
The analysis of BOMD snapshot structures of the solvated
conformers indicates that conformer 1A is slightly more stable
First the structures 1A, B, and C were optimized. It is impor-
tant to note that standard DFT calculations usually completely
Chem. Eur. J. 2014, 20, 6934 – 6938
6937
ꢁ 2014 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim