Interlocking Molecular Gear Chains Built on Surfaces

3 days ago - analyze such behaviors in order to enable the future design and ..... we already have more inclusive data (for three-slave chains) in ...
0 downloads 4 Views 2MB Size
Subscriber access provided by Kaohsiung Medical University

Surfaces, Interfaces, and Catalysis; Physical Properties of Nanomaterials and Materials

Interlocking Molecular Gear Chains Built on Surfaces Rundong Zhao, Fei Qi, Yan-Ling Zhao, Klaus E. Hermann, Rui-Qin Zhang, and Michel A. Van Hove J. Phys. Chem. Lett., Just Accepted Manuscript • DOI: 10.1021/acs.jpclett.8b00676 • Publication Date (Web): 01 May 2018 Downloaded from http://pubs.acs.org on May 3, 2018

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a 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 19 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 Letters

Interlocking molecular gear chains built on surfaces Rundong Zhao1,†, Fei Qi1, Yan-Ling Zhao2, Klaus E. Hermann3, Rui-Qin Zhang2,4 and Michel A. Van Hove1* 1

Institute of Computational and Theoretical Studies & Department of Physics, Hong Kong Baptist University, Hong Kong SAR, China 2 Department of Physics, City University of Hong Kong, Hong Kong SAR, China 3 Inorganic Chemistry Department, Fritz-Haber-Institute der Max-Planck-Gesellschaft, Berlin, Germany 4 Beijing Computational Science Research Center, Beijing, China * Corresponding author; E-mail: [email protected] † Current address: Department of Mechanical Engineering and Materials Science, Duke University, Durham, NC 27707, USA

Abstract Periodic chains of molecular gears in which molecules couple with each other and rotate on surfaces have been previously explored by us theoretically using ab initio simulation tools. Based on the knowledge and experience gained about the interactions between neighboring molecular gears, we here explore the transmission of rotational motion and energy over larger distances, namely through a longer chain of gear-like passive “slave” molecules. Such microscopic gears exhibit quite different behaviors compared to rigid cogwheels in the macroscopic world, due to their structural flexibility affecting intermolecular interaction. Here, we investigate the capabilities of such gear chains and reveal the mechanisms of the transmission process in terms of both quantum-level density functional theory (DFT) and simple classical mechanics. We find that the transmission of rotation along gear chains depends strongly on the gear-gear distance: short distances can cause tilting of gears and even irregular “creep-then-jump” (or “stick-slip”) motion or expulsion of gears; long gear-gear distances cause weak coupling between gears, slipping and skipping. More importantly, for transmission of rotation at intermediate gear-gear distances, our modeling clearly exhibits the relative roles of several important factors: flexibility of gear arms, axles and supports, as well as resulting rotational delays, slippages, and thermal and other effects. These studies therefore allow better informed design of future molecular machine components involving motors, gears, axles, etc. Keywords: molecular machines; molecular gears; gear chains; correlated rotations; ab initio calculations; graphene

Introduction Rotational molecules, ranging from single small molecules to larger nanoscale systems, have been explored widely in the past years due to their potential as molecular motors1-3, and further, in building up larger cooperating molecular machines.4-6 On one hand, molecular motors are 1

ACS Paragon Plus Environment

The Journal of Physical Chemistry Letters 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 19

expected to transfer external energy (the input, i.e. light1, 7, electricity8, 9, chemical energy10-12, etc.) to unidirectional mechanical motion (the output), thus providing a driving force to run molecular machines such as a nanocar13-15. On the other hand, there is still a lack of fundamental understanding that is required to achieve the molecular-scale control of mechanical devices operating at the nanoscale, analogous to the rigid gears/cogwheels in a macroscopic machine system. In the past years, gear(wheel)-like molecules16, self-locked single molecules17 or interlocking molecule networks18 on the nano/molecular scale have been reported; on the micrometer scale, a solid-state SiO2 nano-gear with diameter of around 100 nm was fabricated on highly ordered pyrolytic graphite (HOPG)19. The design of molecular gears is still a largely unexplored area at the current stage. In recent years, we have started investigating this important topic. In a recent publication20 we reported theoretical explorations of basic gear-gear interactions for relatively simple but widely relevant models of two kinds of molecular gears, viz. 3-tooth PF3 molecules attached to a Cu(111) surface and multi-arm (5-tooth) molecular gears attached to a graphene sheet, such that externally powered “driver molecules” induce the rotation of passive “slave molecules”. These molecular gears were computationally simulated using periodic models with an infinite linear chain of alternating driver and slave gears, in which two driver molecules flank each slave molecule. This geometry allowed us to study in significant and useful detail how neighboring molecular gears interact and transmit rotation from one gear to the next. In the present work, we significantly extend those results to examine longer chains of passive slaves: this will illuminate the capabilities and necessary conditions to transmit rotation over longer distances through realistic assemblies of multiple slave gears. In the present work, we theoretically examine longer chains consisting of the same molecular gears that we proposed previously20. These model systems serve as generic examples that allow us to reach general conclusions valid for a much wider range of situations. We show here that such molecules can form interlocking finite one-dimensional chains or even two-dimensional networks of gears that can transmit rotational motion across multiple gears. Compared to macroscopic rigid gears, the softness (flexibility) of the molecular-scale gears, of their axles and of their supports can cause quite different behaviors, such as a “creep-then-jump” or “stick-slip” effect, as well as slipping and skipping. It is therefore important to analyze such behaviors in order to enable the future design and fabrication of effective molecular machines consisting of several linked subsystems such as molecular motors, gears, molecular drive axles, etc. In our modeling, we focus on “low-speed” non-vibrational interactions between gears. The simulations do not include “motor” molecules that convert energy input into driving forces; instead, the rotation will be externally applied to “driver” molecules that can transmit this rotation to a chain or 2-D array of “slave” molecules whose behavior is then studied; our “driver” molecules are thus generic. As such, we are not attempting to model surface systems that might be realized experimentally exactly as modeled, since we are not including the many complications of a motor molecule: the focus is instead on computationally analyzing the nature and behavior of chains of molecular gears. Thus, this work is a theoretical investigation which aims to provide a general theoretical foundation for a variety of future experimental designs, not only for one particular experimental realization.

Modeling and computational methods 2

ACS Paragon Plus Environment

Page 3 of 19 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 Letters

Our previous work20 proposed a 5-arm gear made of a carbon ring (the carbon ring is extended with appropriate chemical groups such as cyano and ethynyl groups or fluorine atoms) which can be attached to a graphene sheet with an intermediate metal atom acting as the pivot between the gear C-ring and a graphene C-ring. This design was inspired by the structures of sandwich compounds that have been studied in organometallic chemistry for several decades.21-26 According to the 18-electron principle which provides stable chemical bonds (in particular the mechanism of the formation of ferrocene21), the intermediate metal atoms in the gear systems cannot be arbitrarily chosen: manganese and its congeners in the Periodic Table have been proved to be the exclusively suitable intermediates for a 5-membered carbon ring to stand firmly but rotate freely on a graphene surface20; similarly, here we prove that chromium, molybdenum and tungsten atoms should be adopted to support a 6-membered carbon ring (see part I of the Supplementary Information). Similar sandwich structures can also be found from experimental syntheses, such as benzene(cyclopentadienyl)manganese(I)27 (C6H6MnC5H5) and (ƞ6-graphene)2Cr28-30. Using DFT calculations, we study the interlocking and rotation process between two or more such molecular gears (like those illustrated in Figure 1 later in the text) by imposing successive rotations on a “driver” molecule (at one end of the gear chain), in order to obtain resulting forces, torques and rotations on the other “slave” molecules of the chain driven by it. In the simulation, we start from an arbitrary “driver” orientation defined as 0˚ in Figure 1 or 2 (or any other similar cases): for this initial driver orientation, while the drivers are kept aligned and fixed, the slaves and the graphene sheet are fully relaxed, yielding negligible net forces and torques on all atoms and molecules. The drivers are then made to rotate stepwise in a smooth quasi-static motion, to simulate how actual motor molecules might rotate. We allow full relaxation of the slave molecules at each imposed driver rotation angle, such that the slaves can follow any appropriate energy-optimized pathway, including rotation, tilting, bending, slipping, breaking up, etc. (the internal driver structure remains rigid during its imposed rotation, so it can also represent other types of motor molecules with possibly different properties). Specifically, the simulation process involves two steps for each driver orientation. In “step 1”, we manually rotate the driver molecule rigidly clockwise by 10˚ about its rotation axis (which is close to being above the center of a carbon ring of the graphene sheet and perpendicular to the sheet), while keeping the slave molecules fixed: this generates a strained unstable situation with non-zero forces and torques acting on the slave molecules. In “step 2”, we re-optimize the slave molecules while keeping the driver rigid and fixed at its 10˚ rotation geometry, resulting in again negligible net forces and torques. This 2-step process is repeated iteratively until the driver molecule reaches the rotation angle of 80˚ (larger than the 72˚ between the arms of a 5-arm gear), which covers the symmetry reduced range of rotations according to the 5-fold symmetry of the 5-arm gear (for 6-arm gears, it is 60˚). For more detailed descriptions, please see our previous work on the mechanisms of gear-gear interactions20. The reason for choosing a rotation step of 10˚ instead of a value closely related to the 5-arm symmetry of the molecule such as 12˚ (360˚/5=72˚, 72˚/6=12˚) is twofold: first, there is an arbitrariness in favoring the 5-fold symmetry of the 5-arm molecules vs. the 6-fold symmetry of the graphene support; second – as shown later in the discussion part – the behavior of the slave molecules is actually asynchronous to the driver’s rotation, with significant angular delays and bending of the arms, so that a symmetry-linked 12˚ step value is not necessarily more meaningful. Since the single molecules considered here show only moderate distortions during rotation, 3

ACS Paragon Plus Environment

The Journal of Physical Chemistry Letters 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 19

the approximation of taking the drivers as rigid bodies is reasonable. However, for a supramolecular system that is very flexible and contains many degrees of rotation and bending, this assumption might significantly restrict the rotational motion pathways and the scenario would often be more complex.31 The finite gear chains are examined by periodic slab calculations using the Vienna Ab-initio Simulation Package (VASP)32, 33, within density functional theory (DFT), and applying the generalized gradient approximation (GGA) of Perdew, Burke, and Ernzerhof (PBE)34. Dispersion corrections due to van der Waals interaction were included according to Grimme35. To simulate the behavior of rotating multi-arm gears built on graphene sheets, the Brillouin zone was sampled by 2  1  1 k-points for the structures shown in Figures 1, 2 and 3, and by 1  1  1 k-points for the 2D networks reported in Figure 5. For all geometry optimizations, the projector-augmented wave (PAW) method36, 37 implemented in the VASP code was adopted, with a maximum 400eV cutoff for the plane-wave basis set.

Results and discussion Gear chains with one driver molecule and several slave molecules. The largest shortcoming of PF3 gears (attached to the Cu(111) surface in our prior work20) is that they have only 3 cogs/arms. This makes the interlocking between gears too flexible and loose: skipping (i.e. jumping of arms past each other) is favored by the wide 120˚ gaps between their 3 arms, by the inherent softness of arm-arm interactions, by the flexibility of molecular arms, as well as by the thermal vibrations (unlike the case of hard metallic teeth in macroscopic gears). For reference, we first consider the shortest chain possible: a driver and a slave (unlike in our earlier work20, here the slave is in contact with only one driver, not two). Compared to our earlier simulation results20 for the same gear-gear distance, we find that, if a slave is only coupled to a single driver instead of two, the weaker interlocking further increases the risk of skipping: with a single driver, the slave only manages to follow the driver until a driver orientation of about 30˚ before stopping further rotation, rather than 60˚ with two drivers (see Figure S1 in the Supplementary Information for illustration). This skipping severely limits the potential use of such 3-arm gears. Most real applications will need more tightly coupled gears, especially when a longer chain of slave gears needs to be driven. To this end, we earlier designed multi-arm gears, using 5-membered carbon rings substituted with cyano groups (-CN) and their capability in transmitting rotation to adjacent gears has been investigated before.20 Those simulations were also conducted using alternating drivers and slaves in infinite periodic linear chains (like those of the PF3 gears).

4

ACS Paragon Plus Environment

Page 5 of 19 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 Letters

Figure 1. Coupling of two adjacent molecular gears of 5-membered carbon rings with cyano (-CN) groups (driver and slave as marked). Each gear is attached to the graphene sheet with an intermediate manganese atom (blue) acting as pivot. The intermolecular distance is equal to 3a (where a is the graphene lattice constant). The CN arms of the slave molecule are labelled A-E. The curved red arrows illustrate their (counterclockwise) rotation while black arrows indicate the (clockwise) rotation of the driver molecules. All molecules have one of their nitrogen atoms highlighted in red to emphasize the rotations. The three panels correspond to three orientations of the driver molecules: 0˚, 30˚ and 60˚, respectively, see text. We again start with the case of 1 driver + 1 slave, as depicted in Figure 1, and will then build upon that simpler case to 1 driver + 1 slave + 1 slave + 1 slave, as depicted later in Figure 2. The gears are placed in close mutual distance (3 carbon rings apart on the graphene sheet, viz. 3 graphene lattice constants) in order to make the intermolecular interaction relatively strong. At closer intermolecular distances (such as 2 carbon rings apart) the graphene surface is computationally found to be torn up due to the excessively strong interaction; in that case, the molecules would not even be pushed apart to other sites, because the binding between the rotor, the metal atom and the graphene sheet is very robust.20 Figure 1 shows that this 5-membered ring by itself can successfully drive a neighboring slave gear at an intermolecular distance of 3a, i.e. 3 graphene lattice constants: the 5-arm slave gear can accomplish a rotation of 1/5 of a full cycle and thus can complete the 360˚ rotation due to its C5 symmetry. The simulation covers the driver’s rotation from 0˚ to 80˚ which guarantees that the slave has gone beyond 72° (1/5 of a full cycle); it also makes clear that a second and subsequent rotations through 72° can take place in this case (as we will show below, a larger spacing between gears may result in “skipping” and thus stopping of the slave gears’ cyclic rotation, which is manifest already in the first period before 72°). We find that the planes of the C5 rings are tilted with respect to the graphene surface, and that the CN bonds are bent (upward or downward) due to the strong intermolecular interaction at this relatively close gear-gear separation, compared with the larger gear-gear separation modeled earlier20, which does not cause any observable tilting. The angles between the CN bonds and the graphene plane range from -10˚ to +35˚: those CN bonds nearer to the driver are significantly bent upward by up to +35˚ and those on the other side are bent down toward the graphene plane by about -10˚, as illustrated in the side view in later Figure 3(A).

5

ACS Paragon Plus Environment

The Journal of Physical Chemistry Letters 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 19

Figure 2. Chain of four adjacent molecular gears of 5-membered carbon rings with cyano (-CN) groups in which the gear at one chain end acts as a driver, rotating clockwise as shown by black arrows, and the other three gears are slaves. The intermolecular distance is 3a. Each gear is attached to the graphene sheet with a central manganese atom (blue) acting as pivot. All molecules have one of their nitrogen atoms highlighted in red to emphasize the rotations. We then elongated the chain by adding more slave gears and observed their behavior. Results for the case of a chain of 3 slave gears are shown in Figure 2 and movie M1. Similar to the 1-slave case of Figure 1, the 3-slave gears tilt and their CN bonds bend by about -10˚ to +30˚ (see also the side view in Figure 3(B)). We observe that, as the driver is rotated, the last slave gear (slave #3) undergoes almost no rotation at first, and then jumps suddenly soon after the driver reaches about 70˚, so as to almost catch up with the driver: thus, slave #3 follows the driver with a “creep-then-jump” motion that is very similar to the “stick-slip” motion familiar in friction and earthquakes, for example. This type of motion also happens to a lesser extent for slave #2, while slave #1 shows a much reduced version of this motion by following the driver rotation more closely. Interestingly, we have never observed in this system that the slave gears escape rotating: it appears that the arms in this system cannot skip past each other, so any delayed rotations are caught up and compensated for later by larger forward jumps. (We also note that the “jump angle” 70˚ of slave #3 is close to 72˚, which corresponds to a fifth of a full cycle, i.e. the repeat angle of a 5-arm gear; however, we argue below that there is no reason for the first jump to occur in general at or near the driver rotation of 72˚, i.e. we believe this proximity to be a mere coincidence and actually to depend on the chain length, i.e. the number of gears in the chain.)

Figure 3. (A) and (B) give side views of the 60˚ rotation geometries in Figures 1 and 2, respectively, looking parallel to the surface and perpendicularly to the chains.

6

ACS Paragon Plus Environment

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

Figure 4. Induced rotation angles of the three slave gears #1, #2 and #3 due to the rotation of the driver gear in Figure 2. The intermolecular distance is 3a. The imposed uniformly changing orientation of the driver gear is shown as the straight black line. The rotation of each slave is defined as the average of the rotations of its five carbon atoms projected onto the graphene surface, since these 3 gears are not parallel with the graphene surface and the individual arm intervals are twisted away from 72˚. To better understand the observed creep-then-jump behavior, the calculated rotation angles are plotted in Figure 4 for each gear in this system: while the driver rotates smoothly by a constant 10˚ at each step (straight black line), the first slave gear (red line) increasingly lags behind the driver by 0 to ~15˚ until the driver reaches ~70˚, then jumps slightly ahead, thus reducing its lag again beyond ~70˚; the next slave gear (slave #2, blue line) greatly magnifies the lag of slave #1, while slave #3 (green line) accentuates it further, almost to the point of not rotating at all before the driver has reached ~70˚, but still jumping to almost catch up beyond 70˚. We see the same behavior repeating itself in the driver orientation range of ~70˚ to ~140˚, and then again after ~140˚. (The range from ~70˚ to ~140˚ again reminds us of the fifth-of-cycle angle of 72˚: in this case, we will argue that the repetition angle is indeed 72˚.) Influence of gear-gear distance. Finite linear chains with a driver molecule on one end followed by a number of slave molecules (such as MnC5(CN)5 gears on graphene) were shown in the last section above to be capable of transmitting rotational motion. However, we also find sizable tilting of the gears, which leads to considerable potential barriers in rotation (up to about 0.3 eV, depending on the degree of tilting, compared with a barrier much smaller than 0.01 eV in the parallel/non-tilting case): these barriers cause irregular “creep-then-jump” or “stick-slip” motion. The tilting is attributed to the relatively small gear-gear distance equivalent to 3 carbon rings on the graphene sheet (viz. 3 graphene lattice constants). Whereas a smaller mutual distance can break up the graphene sheets due to the excessively strong intermolecular interaction, a larger 7

ACS Paragon Plus Environment

The Journal of Physical Chemistry Letters 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 19

mutual distance will reduce the tilting due to weaker intermolecular interaction, but will also increase the risk of skipping between adjacent gears. Thus, the choice of different gear-gear distances needs to be examined carefully and designed rationally so as to reach an optimum situation where the gears can rotate as freely as possible (with small tilting angles and rotational barriers) but can also couple closely with each other. By adopting different equivalent adsorption sites on the graphene lattice (with lattice constant a), gear chains can assume only discrete gear-gear distances: for instance, between the distances 3a and 4a, only √13a (~3.6a) is a possible gear-gear distance on the hexagonal graphene lattice (see Figures 1, 2 and 5 for illustration). For the case of 1 driver and 1 slave, the system with gear-gear distance 3a has been described above (see Figure 1) and shows successful rotation but considerable tilting of the slaves. In contrast, for a gear-gear distance of 4a, our simulations show that the arrangement will lead to skipping between gears. However, for a gear-gear distance of √13a, the chain yields both continuous rotation and negligible tilting of the slaves, making this geometry the most favorable among 2-gear chains (see also Figure S2 in the Supplementary Information for detailed simulation results). An a priori surprising behavior occurs when we elongate the 2-gear √13a chain to 3 gears. While the slave in the 2-gear chain rotates smoothly with the driver (cf. Figure S2), both slaves in the 3-gear chain soon stop rotating (Figure 5(A)). Addition of the 2nd slave behind the 1st slave is sufficient to cause permanent sticking of both slaves once slave #1 has reached an almost fixed orientation such that two of its arms (C and D) are about equidistant to the driver center (cf. Figure 5(A)): Movie M2 in the Supplemental Information helps visualize this behavior. Even more surprising is that this behavior also occurs for the larger gear-gear distance 4a (Figure 5(B)). In this case both slaves stop rotating even sooner, since the additional space between gears allows the driver to slip past slave #1 more easily; also the slaves adopt almost fixed orientations, due to the weaker driver-slave interactions at the larger separation. Results for the 3a distance are not plotted as we already have more inclusive data (for 3-slave chains) in Figure 4: for the case of the 2-slave chain with 3a gear-gear distance, the difference is that the lag in rotation and the jumping of the slaves are much smaller than those exhibited in the 3-slave case.

8

ACS Paragon Plus Environment

Page 9 of 19 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 Letters

Figure 5. Rotation angles of two slave gears #1 and #2 in a 3-gear chain with two different gear-gear distances, due to smooth rotation of the driver gear (see the caption of Figure 4 for details). The gears are separated by (A) √13a and (B) 4a (a is the graphene lattice constant). The geometries at which the skipping between gears becomes most obvious are plotted in the insets: for (A) the skipping occurs by 90º, as denoted with the ellipse on the curves; for (B) the skipping occurs by 60º as denoted. It is intuitively clear that the 2-slave case causes a larger resistance to rotation than the 1-slave case (and even more so the 3-slave case), by cumulating the rotational barriers of both slave molecules: this can roughly double the total rotational barrier which the system has to overcome. Adding even more slaves may magnify this effect further, but not necessarily to unlimited levels: each slave gear can adopt a different “out-of-phase” orientation, so that barriers do not necessarily add up “in phase”; also, more complex behavior can be envisaged for longer chains, such as caterpillar- or snake-like reptation whereby only one slave at a time passes over its barrier, causing a ripple effect that propagates down the chain. Gear-gear distance (unit: a = lattice constant)

Number of slaves

Continuity of rotation

Tilting and bending effect

Rotational barrier (eV)

~0.25

~0.3

3a

1

Yes

The arms on the interacting side bend by +15º to +35º and the carbon ring is almost not tilted. See figure 3(A).

3a

2

Yes

The arms bend strongly by ~+15º to +35º on one side and by around -10º on the other

9

ACS Paragon Plus Environment

The Journal of Physical Chemistry Letters 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 19

interacting side; the carbon rings tilt by ~+5º to +30º for different orientations. See figure 3(B). 3a

3

Yes

As for 2 slaves.

~0.3

√13(≈ 3.6)

1

Yes

No

~0

√13(≈ 3.6)

2

No

No

~0

√13(≈ 3.6) 4a

3

No

No

~0

1

No

No

~0

4a

2

No

No

~0

4a

3

No

No

~0

Table 1. A summary and comparison of the behavior of various chains of MnC5(CN)5 gears on graphene, including the continuity of slave rotation (i.e. whether the slaves continue to rotate or stop rotating due to slipping/skipping), the tilting and bending of the gears (“No” indicates no significant effect), and their rotational barriers. The physical effect of the gear-gear distance can be summarized in a more general picture: (1) Assuming an infinite separation between driver and slave molecules leads to total slippage (i.e., negligible barrier to rotation, and no rotation of the slaves when the driver rotates). (2) When the molecules move closer together, the slave starts to feel the driver’s rotation, but slippage/skipping between the arms/teeth is likely to occur, due to the weak intermolecular interaction. (3) At smaller distances, the gears are coupled more tightly and slippage/skipping is excluded; however, this is also accompanied by increasing tilting and other deformations of the gears, which can increase the barrier to rotation and cause “creep-then-jump” or “stick-slip” behavior. (4) Even closer gear-gear distances may lead to different consequences, such as, in order of severity: gears chemically bond to each other, gears are expelled to other surface sites, gears are detached from the surface, gears are destroyed or the flexible surface (such as a graphene sheet) is torn up. In the process of rotation, there will be a lag due to the softness of the potential and the bending of the arms; this lag accumulates from slave to slave as the gear chain is elongated. The reason for the gears’ sudden jumps (such as those shown in Figure 4 at 70º and 140º) is that the support (e.g. graphene) creates a rotational barrier for the slave, which therefore will refuse to turn over the barrier unless a strong enough force is applied to it by the driver. For different cases explored here, a summary of the rotational behavior is shown in Table 1 for comparison.

Mechanisms of rotation from the perspectives of total energy and torque analysis. Taking the most interesting case of the gear chain shown in Figure 2 as an example, it is clear that the primary cause of the lags in rotation of slaves is that the arms of one gear do not fit snuggly between the arms of the next gear, due to the large opening angle between arms, as well as the softness of the arm-arm interaction and the distortions of the gear arms, of the gear axles and of the support. The resulting flexibility leaves “wiggle” room that causes an angular delay in the slave gears’ rotation (macroscopic metallic gears have negligible softness and flexibility but they do experience a small amount of free wiggle room before the hard teeth push each other, which still causes a small delay in rotation, usually softened by lubrication). This angular delay adds up from one slave gear to the next, but not simply linearly due to the additional influence of individual rotation barriers. A large contribution to the “sticky” part of the stick-slip behavior (“creep-then-jump” effect) is the tilting 10

ACS Paragon Plus Environment

Page 11 of 19 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 Letters

of the molecules, especially of the last slave which is relatively free to tilt: the tilting makes the arms scrape against the graphene support, which has a bumpy corrugated surface with significant barriers to rotation of the arms. Therefore, the more tilted slaves jump abruptly over those graphene-induced barriers rather than co-rotating smoothly with the driver. This behavior is akin to the jerky stick-slip motion of chalk that screeches across a blackboard or of geologic faults that slide slowly between sudden earthquake events. The stick-slip behavior will also be a major contribution to energy dissipation during rotation, as the sudden jumps can easily excite vibrations in gears and support. We see that the soft interaction, flexibility and tilting prevent transmitting the rotation rigidly and immediately through a long molecular gear chain. Instead, as the driver rotates, stress and distortion build up gradually along the chain, especially when slaves are blocked from rotating by the support due to tilting: this distortion is most clearly visible in Figure 2 at 60˚ in the form of lateral displacements of the slave gear centers from the underlying C6-ring centers, as the slave gears are pushed aside by the force from the driver gear; this also implies tilting of the rotation axle away from the surface normal, as shown in Figure 3; in addition, there is some bending of the arms of the gears, both in and out of the plane of the molecule. The stress is thus more or less evenly distributed as elastic energy along the chain length. The accumulated tension is suddenly released when it has grown large enough to finally overcome the barrier for rotating the more distant slaves, as seen from the red curve (“Step 2”) of Figure 6: the gradually increasing total energy from around 20˚ to 70˚ and from around 90˚ to 140˚ reflects the accumulation of the stress as the slaves increasingly lag behind, while its sudden drops around 70˚ and 140˚ correspond to the jumps shown in Figure 4.

Figure 6. The total torque of slave #1 and the total energy of the whole 4-gear chain system shown in Figure 2. As detailed in the text, the simulation at each rotation angle consists of two 11

ACS Paragon Plus Environment

The Journal of Physical Chemistry Letters 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 19

steps: (1) rotate the driver clockwise by 10˚ while keeping all the slaves fixed ‒ the black and blue curves indicate, respectively, the total energy of the full system and the total torque acting on slave #1 at this step; (2) re-optimize all the slaves while fixing the driver at its rotated geometry (see red curve for the total energy). As detailed in the Modeling and computational methods section, our stepwise quasi-static model imposes driving forces (which can be reflected by the torque) in the 1st step, followed in step 2 by full relaxation of the slaves. The energy difference between step 1 and step 2 can thus also indirectly reflect the strength of the driving forces: at most rotation angles (30˚-170˚), the difference in energy ranges from 0.29 to 0.85 eV (6.69 to 19.60 kcal/mol), much larger than that in the PF3 gears (0.02 to 0.09 eV, or 0.43 to 2.11 kcal/mol)20. The torque analysis scheme proposed in our previous work38 provides a very convenient and simple method to measure the rotation direction and rotation strength of rotor molecules: instead of a geometry optimization that involves many steps of self-consistent field (SCF) cycles, the Pulay forces obtained from a single-point SCF calculation can serve to estimate the rotational trend of the rotor molecule. The application to molecules of the classical-mechanics concept of torque is analogous to the common use of quantum-based forces and Newtonian mechanics in ab initio geometric optimization and in molecular dynamics38. Taking as an example the 5-arm gear which can rotate about the manganese atom (the pivot), the total torque acting on this molecule (and thus inducing rotation) can be defined simply as:

 = ∑    ,

 is the where  runs over all the atoms (10 C atoms and 5 N atoms) that constitute the rotor,  Pulay force acting on atom i, and  is the corresponding vector from the pivot (the manganese atom) to the position of atom i. In the present case, the effective torque component perpendicular to the graphene sheet is of most interest, while its other components lead to the tilting of the rotor. The first step of the 2-step process creates driving forces that can be observed only on the still-fixed slave #1: forces observed on the atoms of slave #2 and #3 at step 1 are almost zero in the present scheme, because they do not directly interact with the driver. However, such forces directly imposed on slave #1 can later (in step 2) be transformed to stresses that more evenly spread along the chain. Consequently, the strengths of such forces on slave #1 ‒ and thus the effective torque they cause ‒ are quite helpful for us to understand the strength of the direct driving force which is expected to affect the whole chain. This effective torque is plotted as the blue curve in Figure 6 and is shown to be positive at all rotation angles, implying a unidirectional counterclockwise rotation of slave #1. Its value ranges from 0.59 to 7.95 nN·nm, compared to 0.18 to 0.30 nN·nm in the PF3 case20, showing a several times stronger driving force. The shape of the torque curve also reflects the process of accumulating and releasing the stress: the torque on slave #1 grows to produce a larger stress and larger total energy, while the torque almost vanishes during the slave’s jumps. With the revealed mechanism, the longer we make the chain, the harder it will be to transmit the driver rotation to the most distant slave gear: in the particular case of our model system, four or more slave gears will thus require an input driver rotation angle considerably larger than ~70˚ to build up enough tension along the chain to make the most distant slave gears rotate; we expect the required input rotation to be roughly proportional to the chain length, as the elastic energy is 12

ACS Paragon Plus Environment

Page 13 of 19 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 Letters

spread out along the length of the chain but needs to exceed a critical value near the last gears. Nevertheless, our system has 72˚ rotational periodicity and this implies that subsequent rotational jumps (after the first one) will repeat themselves for every 72˚ of input driver rotation angle, as seen in Figure 4 for the second jump. To be more precise, we observe in Figure 4 that the initial 70˚ or so of driver rotation induces and builds up an initial lag in the slave rotation: in this phase the driver pushes against the slaves which resist the slaves’ own first rotation barrier (partly due to the more distant slaves), causing the slave arms to bend under the applied force and torque and the slave rotation to lag increasingly behind that of the driver. After the initial 60˚ of driver rotation, the slaves start to yield and then jump suddenly over the barrier around 70˚-80˚. Then, in the range 80˚-100˚, the slaves still catch up to the driver orientation, but more slowly and gradually. The range from ~40˚ to ~100˚ is then seen to repeat itself almost identically as a second cycle (110˚-170˚), indicating that already by about 40˚ the system has adopted a periodically repeating rotation mode in which the initial lag is incorporated and maintained. A deeper understanding of the strength of such intermolecular interactions needs an energy decomposition scheme39-41, by which the intermolecular interaction can be decomposed into four contributions: the direct Coulomb interaction between the two unperturbed molecules, the Pauli repulsion between the overlapping electrons, the orbital interaction, and the van der Waals interaction. Many analyses on systems very similar to those investigated here have been reported previously20, 42, 43, thus we do not perform the analogous calculations here.

Other analogous cases: 6-arm gears and gear networks formed with them. The studies illustrated in Figures 1 and 2 adopt 5-arm molecular gears; several 6-arm gears that have C6 symmetry (and thus are more in accord with the C6 symmetry of graphene) also transmit rotation well; since the 6-arm gears are found to behave very much like the 5-arm gears, we here focus on other aspects to expand our understanding. We start with the simple case of benzene rings substituted with 6 cyano (-CN) groups, modeled using periodic cells – each unit cell contains 1 driver and 1 slave and the gear chain contains alternating drivers and slaves, spaced 4a apart. We used chromium atoms as the pivots, as discussed in the Modeling and computational methods section. This simulation is a logical extension of our previous work20 that modeled an infinite chain consisting of alternating 5-arm drivers and slaves. As expected, the slaves keep pace with the rotating drivers at all rotation angles, as shown in Figure S3 in section IV of the Supplementary Information. We further modeled gears with longer arms: the behavior of the 6-carbon benzene ring substituted with 6 ethynyl (-C≡CH) groups44 and chromium pivot atoms was also studied using the same scheme. Results are shown in Figure S4. This gear has 6 arms with each arm being longer than the -CN arms in the former studies. Its performance is found to be similar to that of the C5-based gear except that the longer ethynyl arms can bend tangentially relatively more during the rotation (the slave arms are often somewhat curved and the angles between them are not exactly 72˚).

13

ACS Paragon Plus Environment

The Journal of Physical Chemistry Letters 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 19

Figure 7. Two typical models of gear networks consisting of C6(CN)6 molecules, using √13 gear-gear distances: (A) 1 driver at the center of the rectangular structure can directly drive 4 slaves. (B) 1 driver at one corner of the square structure can drive 3 slaves through two converging linkage paths. The purple atoms denote chromium atoms ‒ the pivots of the 6-arm gears. Both (A) and (B) demonstrate the configuration of the gears at the driver orientation of 0˚. See Figure S5 for the results at 20˚ and 40˚. Our main interest here is the behavior of finite non-periodic chains of such gears. Taking the C6 symmetry into account, we adopt C6 (CN)6 molecules to mesh into 2D networks consisting of 1 driver surrounded by 4 slaves (see Figure 7(A)), or 1 driver at the corner of a diamond and 3 slaves at its other corners (see Figure 7(B)). Such structures are commonly used in macroscopic gear systems existing in e.g. clockworks, gear boxes, etc. As demonstrated in Figure 7 (and Figure S5, for more details at the rotation angles of 20˚ and 40˚ of the drivers), these mesh structures function well: Figures 7(A) and S5(A) show that a single driver is able to drive another slave without skipping/slippage between their arms, and such 6-arm gears can be logically extended to finite gear chains like those shown in Figures 1 and 2; Figures 7(B) and S5(B) show that a driver at one corner can transmit its rotational motion to other gears via two paths (or even more, if needed), and this propagation can logistically be extended to larger 2D square or rectangular networks, based on the conclusions obtained from the aforementioned studies. We see that the slaves shown in Figures 7 and S5 are slightly shifted off their centers (the pivots, viz. the chromium atoms) due to the repulsion between the gears, but this does not hinder the gear-gear cooperation.

Conclusions The present study, based on and extending our earlier work20, investigates several prototypes of finite molecular gear chains. We show that gears with more arms couple better. Thus, the triple-arm gears studied in the previous work, such as PF3 or NH3 molecules attached to copper surfaces, are less effective in driving multiple slaves with only 1 driver, due to their sparse teeth. Multi-arm gears, e.g. 5-membered carbon rings substituted with 5 -CN groups or 6-membered carbon rings substituted with 6 -CN or 6 -C≡CH groups attached to graphene sheets, are better 14

ACS Paragon Plus Environment

Page 15 of 19 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 Letters

able to transmit rotation. However, we also find that transmission of rotation beyond a few slave gears exhibits a “creep-then-jump” or “stick-slip” effect between gears, resulting both from soft and elastic interactions between and within gears, and from the scraping of tilted molecules against the supporting surface. In the present work, we have modeled the cases of chains of 1 driver and 1, 2 and 3 slaves, as well as two extended 2D network models. For chains with more slaves the conclusion can be deduced logically: inasmuch as the gears are closely coupled (and thus skipping is prevented), the slave at the extremity of the chains can always follow the driver, even if with irregular rotation – as long as enough input driving force has been imposed. It is also noted that the assigned gear-gear distance can strongly influence the transmission of rotation, by affecting the strength of the intermolecular interaction and by generating the tilting of gears. We note that a quasi-classical picture of these rotations is realistic and intuitive at a qualitative level: there is no need to invoke quantum mechanics to predict the general behavior of the gear rotation and gear-gear coupling. At the quantitative level, quantum effects clearly underlie these motions, e.g. for rotational and translational barrier heights, while quantum tunneling, thermal effects and energy dissipation will also have quantitative consequences. Such studies can help designing gear systems which are indispensable in composing complicated machine systems, and may contribute to a blossoming of the development of nanoscale molecular robotics.

Acknowledgements This work was supported by the Collaborative Research Fund of Hong Kong Research Grants Council (Project No. C2014-15G). We also acknowledge the computing resources of the Tianhe2-JK cluster at the Beijing Computational Science Research Center and of the Tianhe2 cluster at the National Supercomputer Center in Guangzhou, China. K. Hermann is grateful for support by the Institute of Computational and Theoretical Studies (ICTS)at HKBU; ICTS is supported by the Institute of Creativity, which is sponsored by Hung Hin Shiu Charitable Foundation (孔憲紹慈善基金贊助).

Supporting Information Available I. Stability of the 6-membered carbon ring supported by a chromium atom or its congeners; II. limitation of the PF3 gear: skipping (slippage) between arms of different gears in a finite gear chain; III. the gear chain consisting of 1 driver and 1 slave with a gear-gear distance of √13a; IV. 6-arm gears with chromium pivots; V. 2D networks made with 6-arm gears; VI. Movies of gear rotations within molecular chains of 3 and 2 slave molecules.

References 15

ACS Paragon Plus Environment

The Journal of Physical Chemistry Letters 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.

Page 16 of 19

Koumura, N.; Zijlstra, R. W. J.; van Delden, R. A.; Harada, N.; Feringa, B. L., Light-driven

monodirectional molecular rotor. Nature 1999, 401, 152-155. 2.

van Delden, R. A.; ter Wiel, M. K. J.; Pollard, M. M.; Vicario, J.; Koumura, N.; Feringa, B. L.,

Unidirectional molecular motor on a gold surface. Nature 2005, 437, 1337-1340. 3.

Peplow, M., The tiniest lego: a tale of nanoscale motors, switches and pumps. Nature 2015, 525,

18-21. 4.

Kay, E. R.; Leigh, D. A.; Zerbetto, F., Synthetic molecular motors and mechanical machines. Angew.

Chem. Int. Ed. 2007, 46, 72-191. 5.

Izatt, R. M., Macrocyclic and supramolecular chemistry. Wiley 2016.

6.

Balzani, V.; Credi, A.; Venturi, M., Molecular devices and machines: concepts and perspectives for

the nanoworld. Wiley 2008. 7.

Balzani, V.; Clemente-León, M.; Credi, A.; Ferrer, B.; Venturi, M.; Flood, A. H.; Stoddart, J. F.,

Autonomous artificial nanomotor powered by sunlight. Proc. Natl Acad. Sci. USA 2006, 103, 1178-1183. 8.

Tierney, H. L.; Murphy, C. J.; Jewell, A. D.; Baber, A. E.; Iski, E. V.; Khodaverdian, H. Y.; McGuire, A.

F.; Klebanov, N.; Sykes, E. C. H., Experimental demonstration of a single-molecule electric motor. Nature Nanotechnology 2011, 6, 625-629. 9.

Horinek, D.; Michl, J., Molecular Dynamics Simulation of an Electric Field Driven Dipolar

Molecular Rotor Attached to a Quartz Glass Surface. J. Am. Chem. Soc. 2003, 125, 11900-11910. 10. Kelly, T. R.; Silva, R. A.; De Silva, H.; Jasmin, S.; Zhao, Y., A Rationally Designed Prototype of a Molecular Motor. J. Am. Chem. Soc. 2000, 122, 6935-6949. 11. Hernández, J. V.; Kay, E. R.; Leigh, D. A., A reversible synthetic rotary molecular motor. Science 2004, 306, 1532-1537. 12. Fletcher, S. P.; Dumur, F.; Pollard, M. M.; Feringa, B. L., A Reversible, Unidirectional Molecular Rotary Motor Driven by Chemical Energy. Science 2005, 310, 80-82. 13. Kudernac, T.; Ruangsupapichat, N.; Parschau, M.; Macia, B.; Katsonis, N.; Harutyunyan, S. R.; Ernst, K.-H.; Feringa, B. L., Electrically driven directional motion of a four-wheeled molecule on a metal surface. Nature 2011, 479, 208-211. 14. Shirai, Y.; Osgood, A. J.; Zhao, Y.; Kelly, K. F.; Tour, J. M., Directional Control in Thermally Driven Single-Molecule Nanocars. Nano Letters 2005, 5, 2330-2334. 15. Morin, J.-F.; Shirai, Y.; Tour, J. M., En Route to a Motorized Nanocar. Organic Letters 2006, 8, 1713-1716. 16. Nakamura, M.; Kishimoto, K.; Kobori, Y.; Abe, T.; Yoza, K.; Kobayashi, K., Self-Assembled Molecular Gear: A 4:1 Complex of Rh(III)Cl Tetraarylporphyrin and Tetra(p-pyridyl)cavitand. J. Am. Chem. Soc. 2016, 138, 12564-12577. 17. Štacko, P.; Kistemaker, J. C. M.; van Leeuwen, T.; Chang, M.-C.; Otten, E.; Feringa, B. L., Locked synchronous rotor motion in a molecular motor. Science 2017, 356, 964. 18. Zhang, Y.; Kersell, H.; Stefak, R.; Echeverria, J.; Iancu, V.; Perera, U. G. E.; Li, Y.; Deshpande, A.; Braun, K.-F.; Joachim, C.; Rapenne, G.; Hla, S.-W., Simultaneous and coordinated rotational switching of all molecular rotors in a network. Nature Nanotechnology 2016, 11, 706-712. 19. Yang, J.; Deng, J.; Troadec, C.; Ondarçuhu, T.; Joachim, C., Solid-state SiO2 nano-gears AFM tip manipulation on HOPG. Nanotechnology 2014, 25, 465305. 20. Zhao, R.; Zhao, Y.-L.; Qi, F.; Hermann, K.; Zhang, R.-Q.; Van Hove, M. A., Interlocking mechanism between molecular gears attached to surfaces. ACS Nano 2018, 12, 3020-3029, see also arXiv-Version: 16

ACS Paragon Plus Environment

Page 17 of 19 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 Letters

http://arxiv.org/abs/1802.01802. 21. Dunitz, J. D.; Orgel, L. E.; Rich, A., The crystal structure of ferrocene. Acta Cryst. 1956, 9, 373. 22. Hudson, R. D. A., Ferrocene polymers: current architectures, syntheses and utility. Journal of Organometallic Chemistry 2001, 637, 47-69. 23. Rayón, V. M.; Frenking, G., Bis(benzene)chromium is a δ-bonded molecule and ferrocene is a π-bonded molecule. Organometallics 2003, 22, 3304. 24. Ramos, E.; Martínez, A.; Rios, C.; Salcedo, R., Which is the best sandwich compound? Hexaphenylbenzene substituted by sandwich compounds bearing Sc, Cr, and Fe. Phys. Chem. A 2015, 119 (47), 11523. 25. Écija, D.; Auwärter, W.; Vijayaraghavan, S.; Seufert, K.; Bischoff, F.; Tashiro, K.; Barth, J. V., Assembly and Manipulation of Rotatable Cerium Porphyrinato Sandwich Complexes on a Surface. Angewandte Chemie International Edition 2011, 50, 3872-3877. 26. Komeda, T.; Isshiki, H.; Liu, J.; Zhang, Y.-F.; Lorente, N.; Katoh, K.; Breedlove, B. K.; Yamashita, M., Observation and electric current control of a local spin in a single-molecule magnet. Nature Communications 2011, 2, 217. 27. Denning, R. G.; Wentworth, R. A. D., A Mass Spectroscopic Study of Benzene(cyclopentadienyl)manganese(I). J. Am. Chem. Soc. 1966, 88, 4619-4622. 28. Sarkar, S.; Niyogi, S.; Bekyarova, E.; Haddon, R. C., Organometallic chemistry of extended periodic π-electron systems: hexahapto-chromium complexes of graphene and single-walled carbon nanotubes. Chemical Science 2011, 2, 1326. 29. Bekyarova, E.; Sarkar, S.; Wang, F.; Itkis, M. E.; Kalinina, I.; Tian, X.; Haddon, R. C., Effect of covalent chemistry on the electronic structure and properties of carbon nanotubes and graphene. Acc. Chem. Res. 2013, 46, 65. 30. Zhang, Z.; Turner, C. H., Structural and electronic properties of carbon nanotubes and graphenes functionalized with cyclopentadienyl-transition metal complexes: a DFT study. J. Phys. Chem. C 2013, 117, 8758. 31. Palma, C.-A.; Björk, J.; Rao, F.; Kühne, D.; Klappenberger, F.; Barth, J. V., Topological Dynamics in Supramolecular Rotors. Nano Letters 2014, 14, 4461-4468. 32. Kresse, G.; Furthmüller, J., Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set. Phys. Rev. B 1996, 54, 11169. 33. Kresse, G.; Furthmüller, J., Efficiency of ab-initio total energy calculations for metals and semiconductors using a plane-wave basis set. Comput. Mater. Sci. 1996, 6, 15. 34. Perdew, J. P.; Burke, K.; Ernzerhof, M., Generalized gradient approximation made simple. Phys. Rev. Lett. 1996, 77, 3865. 35. S. Grimme; Ehrlich, S.; Goerigk, L., Effect of the damping function in dispersion corrected density functional theory. J. Comput. Chem. 2011, 32, 1456-1465. 36. Blöchl, P. E., Projector augmented-wave method. Phys. Rev. B 1994, 50, 17953. 37. Kresse, G.; Joubert, D., From ultrasoft pseudopotentials to the projector augmented-wave method. Phys. Rev. B 1999, 59, 1758. 38. Zhang, R.-Q.; Zhao, Y.-L.; Qi, F.; Hermann, K.; Van Hove, M. A., Intramolecular torque, an indicator of the internal rotation direction of rotor molecules and similar systems. Phys. Chem. Chem. Phys. 2016, 18, 29665; see also arXiv-Version http://arxiv.org/abs/1607.00593. 39. Morokuma, K., Molecular orbital studies of hydrogen bonds. III. C=O···H-O hydrogen bond in H2CO···H20 and H2CO···2H20. J. Chem. Phys. 1971, 55, 1236-1244. 17

ACS Paragon Plus Environment

The Journal of Physical Chemistry Letters 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 19

40. Ziegler, T.; Rauk, A., On the calculation of bonding energies by the Hartree Fock Slater method. Theoret. Chim. Acta 1977, 46, 1-10. 41. Bagus, P. S.; Hermann, K.; Bauschlicher Jr., C. W., On the nature of the bonding of lone pair ligands to a transition metal J. Chem. Phys. 1984, 81, 1966. 42. Zhao, R.; Zhang, R. Q., A new insight into π–π stacking involving remarkable orbital interactions. Phys. Chem. Chem. Phys. 2016, 18, 25452-25457. 43. Zhao, R.; Zhang, R.-Q., Beyond the electrostatic model: the significant roles of orbital interaction and the dispersion effect in aqueous-π systems. Phys. Chem. Chem. Phys. 2017, 19, 1298. 44. Diercks, R.; Armstrong, J. C.; Boese, R.; Vollhard, K. P. C., Hexaethynylbenzene. Angew. Chem. Int. Ed. 1986, 25, 268-269.

18

ACS Paragon Plus Environment

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

TOC GRAPHIC

ACS Paragon Plus Environment