C O M M U N I C A T I O N S
These conclusions on stability of the prochiral or chiral dimers
of 1 have been further confirmed by MD simulations including
solvent molecules. For this, the solvent was introduced in the model
by adding two layers of solvent molecules on the top of the surface
and by using the dielectric constant of 1-phenyloctane. The results
show that both dimers are stabilized by a few kcal/mol due to the
presence of the solvent molecules. However, since the stabilization
is comparable for both dimers, the PP dimer remains significantly
more stable than the PM dimer. For the latter, the kinetics of the
dissociation slows down, probably due to a cage effect of the solvent
molecules. Finally, an additional indication about the structure of
the dimer building block of the monolayer arises from the fact that
with PM dimers it is not possible to build a monolayer matching
the structure observed with STM. In contrast, the monolayer of PP
dimers as modeled by force-field simulations is in good agreement
with the experiment.
In conclusion, we have shown that chiral domains are formed at
the liquid/solid interface by achiral diarylethenes because of surface-
induced atropoisomerism. On the surface of HOPG, these molecules
with a locked helicity form dimers by shape complementarity.
Molecular dynamics show that these dimers of diarylethenes are
chiral dimers (PP and MM), rather than prochiral dimers (PM and
MP). These results add to our understanding of how chirality can
emerge in 2D systems.
Figure 3. Molecular dynamics trajectory showing the molecule-molecule
binding energy for two types of dimers (heterodimers PM and homodimers
PP) on HOPG, at room temperature, extracted from the total energy of the
system. Blue line: PM dimer. Red line: PP dimer. The curve shows that,
after 20 ps, the heterodimer dissociates, whereas the homodimer remains
stable.
Table 1. Average Molecule-Molecule Binding Energy and
Standard Deviation for the Two Dimers on Graphite, Based on the
MD Simulations Shown in Figure 3
dimer
binding energy (kcal/mol)
STD (kcal/mol)
PM
PP
-1.85
-9.27
2.21
1.45
conformer is unlikely to allow for stabilizing molecule-molecule
or molecule-substrate interactions. Moreover, as all experiments
are carried out in the dark, we exclude that 1 adsorbs in a
photoisomerized form. We therefore assume that 1 exclusively
physisorbs in its antiparallel conformation on HOPG.
Acknowledgment. This work has been supported by The
Netherlands Organization for Scientific Research (NWO-CW)
through a VENI grant (N.K.), the Marie Curie RTN CHEXTAN
(MRTN-CT-2004-512161), FNRS, the Fund of Scientific Research-
Flanders (FWO), and the Belgian Federal Science Policy Office
through IAP-6/27.
Dimer formation, as observed experimentally, suggests that an
energetically favorable interaction occurs between either prochiral
(PM/MP) or chiral (PP and MM) pairs of adsorbed molecules. In
order to understand the stereochemistry of these dimeric building
blocks, molecular dynamics (MD) have been performed for PP and
for PM dimers on the HOPG surface, at rt and in a dry situation,
that is, without solvent molecules.13 The calculations show a
dissociation of the prochiral dimer on the surface, while the chiral
dimer remains intact, due to favorable interactions between the
molecules and the surface; for example, the alkyl chains can be
fully adsorbed on graphite and among the molecules themselves.
In order to evaluate the strength of the interactions between the
molecules, the molecule-molecule binding energy for homodimers
and heterodimers on HOPG (Figure 3) has been extracted from
the total energy of the system (including the surface binding energy),
over a MD trajectory of 50 ps. From these curves, it appears that
on graphite the PP dimer is clearly more stable than its PM
counterpart, which dissociates after 20 ps. High-temperature simula-
tions at 400 K confirmed this trend. Table 1 summarizes these
results by showing the average binding energy and the standard
deviation for the two types of dimers. A detailed analysis of the
components of the binding energies of the dimers shows that the
main difference in their stability arises from weak van der Waals
interactions acting between the two helices. Therefore, we propose
that the difference of stability between the dimers is due to a better
locking between the conjugated heads of the molecules in the PP
(and MM) dimers rather than from a better adsorption of the alkyl
chains on the surface of HOPG. This means that both the molecular
locking, in terms of relative position of the helices (which determine
the shape of the dimer), and the binding energy depend on the
details of the molecular geometries. The fact that PP and MM
dimers are the most stable dimers added to the fact that the two
domains of Figure 2a are mirror images and form an angle which
is not 60° with respect to each other suggest that these domains
are formed by either PP or MM dimers.
Supporting Information Available: Synthesis and characterization
of 1, schematic representation of its parallel conformation, and
description of STM experiments. This material is available free of
References
(1) Kitzerow, H.-S.; Bahr, C. Chirality in Liquid Crystals; Springer: Berlin,
2001.
(2) (a) Eelkema, R.; Pollard, M. M.; Vicario, J.; Katsonis, N.; Serrano Ramon,
B.; Bastiaansen, C. W. M.; Broer, D. J.; Feringa, B. L. Nature 2006, 440,
163. (b) Engelkamp, H.; Middelbeek, S.; Nolte, R. J. M. Science 1999,
284, 785-788.
(3) (a) Perez-Garcia, L.; Amabilino, D. B. Chem. Soc. ReV. 2007, 36, 941-
967. (b) Ernst, K. H. in Supramolecular Chirality; Crego-Calama, M.,
Reinhoudt, D., Eds.; Topics of Current Chemistry 269; Springer-Verlag:
Berlin, 2006; p 209 and references cited therein.
(4) (a) Charra, F.; Cousty, J. Phys. ReV. Lett. 1998, 80, 1682. (b) Weckesser,
J.; De Vita, A.; Barth, J. V.; Cai, C.; Kern, K. Phys. ReV. Lett. 2001, 87,
96101. (c) Messina, P.; Dmitriev, A.; Lin, N.; Spillmann, H.; Abel, M.;
Barth, J. V.; Kern, K. J. Am. Chem. Soc. 2002, 124, 14000-14001.
(5) De Feyter, S.; De Schryver, F. C. In Scanning Probe Microscopies Beyond
Imaging, Manipulation of Molecules and Nanostructures; Samori, P., Ed.;
Wiley-VCH: New York, 2006; pp 3-35.
(6) Fasel, R.; Parschau, M.; Ernst, K. H. Nature 2006, 439, 449-452.
(7) Blum, M. C.; Cavar, E.; Pivetta, M.; Patthey, F.; Schneider, W. D. Angew.
Chem., Int. Ed. 2005, 44, 5334-5337.
(8) For stereochemical definitions, see: Eliel, E. L.; Wilen, S. H. Stereo-
chemistry of Organic Compounds; Wiley: New York, 1994.
(9) A schematic representation of the parallel conformation of diarylethene
1 is shown in Supporting Information.
(10) (a) de Jong, J. J. D.; Lucas, L. N.; Kellogg, R. M.; van Esch, J. H.; Feringa,
B. L. Science 2004, 304, 278-281. (b) de Jong, J. J. D.; Tiemersma-
Wegman, T. D.; van Esch, J. H.; Feringa, B. L. J. Am. Chem. Soc. 2005,
127, 13804-13805.
(11) (a) Katsonis, N.; Marchenko, A.; Fichou, D. J. Am. Chem. Soc. 2003,
125, 13682-13683. (b) Lacaze, E.; Alba, M.; Goldman, M.; Michel, J.
P.; Rieutord, F. Eur. Phys. J. B 2004, 39, 261-272.
(12) Irie, M. Chem. ReV. 2000, 100, 1685-1716.
(13) Molecular dynamics simulations have been performed using the free
molecular modelling package TINKER4.2 and the MM3(2000) force field.
JA075917D
9
J. AM. CHEM. SOC. VOL. 130, NO. 2, 2008 387