Enhancing the Intermolecular Interaction by Cyano Substitution in CuPc

†Department of Physics and Center for Energy Science, Indian Institute of Science Education and. Research (IISER) ... and scanning tunneling spectro...
1 downloads 8 Views 1MB Size
Subscriber access provided by ECU Libraries

Article

Enhancing the Intermolecular Interaction by Cyano Substitution in CuPc Rejaul Sk, Srilatha Arra, Barun Dhara, Joel S. Miller, Mukul Kabir, and Aparna Deshpande J. Phys. Chem. C, Just Accepted Manuscript • DOI: 10.1021/acs.jpcc.7b09225 • Publication Date (Web): 15 Dec 2017 Downloaded from http://pubs.acs.org on December 16, 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.

The Journal of Physical Chemistry C 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 31 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

Enhancing the Intermolecular Interaction by Cyano Substitution in CuPc Rejaul Sk,# Srilatha Arra,$ Barun Dhara,$ Joel S. Miller,§ Mukul Kabir,† Aparna Deshpande†#*

#

Scanning Tunneling Microscopy and Atom Manipulation Laboratory, Department of Physics, h-

cross, Indian Institute of Science Education and Research (IISER) Pune 411008, India $

Department of Chemistry, Indian Institute of Science Education and Research (IISER) Pune

411008, India †

Department of Physics and Center for Energy Science, Indian Institute of Science Education and

Research (IISER) Pune 411008, India §

Department of Chemistry, University of Utah, Salt Lake City, UT 84112-0850, USA

Author Information *[email protected]

ACS Paragon Plus Environment

1

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 2 of 31

ABSTRACT: On-surface molecular self-assembly is one of the key paradigms for understanding intermolecular interactions and molecule-substrate interactions at

the atomic scale.

Phthalocyanines are planar π conjugated systems capable of self-assembly and can act as versatile, robust, and tunable templates for surface functionalization. One of the ways to tailor the properties of phthalocyanines is by pendant group substitution. How such a scheme brings about changes in the properties of the phthalocyanines at the nanoscale has not been greatly explored. Here we present an atomic scale picture of the self-assembly of copper phthalocyanine, CuPc, and compare it with its cyano analogue, CuPc(CN)8,on Au(111) using scanning tunneling microscopy (STM) and scanning tunneling spectroscopy (STS) in ultra-high vacuum (UHV) at 77 K. STM imaging reveals a tetramer unit cell to be the hallmark of each assembly. The periodicity of herringbone reconstruction of Au(111) is unchanged upon CuPc(CN)8 adsorption whereas for CuPc adsorption this periodicity changes. STS measurements show an increment in the HOMO-LUMO gap from CuPc to CuPc(CN)8. Extensive ab initio calculations within the density functional theory (DFT) match well with the experimental observations. The STM imaging shows adsorption induced organizational chirality for both assemblies. For CuPc(CN)8 at LUMO energy the individual molecule exhibits an orbital energy dependent chirality on top of the existing organizational chirality. It remains achiral at HOMO energy and within the HOMO-LUMO gap. No such peculiarity is seen in the CuPc assembly. This energy-selective chiral picture of CuPc(CN)8 is ascribed to the cyano groups which participate in the antiparallel dipolar coupling thereby enhancing the intermolecular interaction in the CuPc(CN)8 assembly. Thus our atomically resolved topographic and spectroscopic studies supplemented by the DFT calculations demonstrate that pendant group substitution is an effective strategy for tweaking intermolecular interactions and for surface functionalization.

ACS Paragon Plus Environment

2

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

1. INTRODUCTION

Self-assembly is a ubiquitous phenomenon and has its hallmark in various systems at all length scales ranging from carbon atoms in interstellar medium1, organic molecules in the interstellar ices2, proteins in cephalopods3 to three dimensional photonic crystals4. Molecular self-assembly in solution as well as in solid phase relies on collective strength of non-covalent interactions5. For polymerization induced by visible light, nanoparticles in different morphologies of self-assembly can be realized in solution6, whereas for substrate-supported self-assembly of diacetylene molecules, polymerization using ultraviolet (UV) light has been demonstrated to template the surface7. Conjugated polymers with their propensity for supramolecular self-assembly have proved to be promising for applications in optoelectronic devices8-11. Surface self-assembly entails a rich mix of physisorption and chemisorption giving rise to a co-ordinated network of atoms and molecules on surfaces. It indulges in covalent as well as non-covalent bonds, thereby making this supramolecular

chemistry

indispensable

for

interdisciplinary

advances

towards

new

functionalities12. In case of organic molecular self-assembly, phthalocyanines (Pcs) are one of the prime members of organic functional materials with applications in molecular photovoltaics. Their potential to enhance the efficiency of solar cells by using different anchoring groups or ligands has been well-demonstrated13. Surface-supported Pcs like (RuPc)2 dimers are promising entities for spintronics and nanoscale devices14. The introduction of functional groups to the periphery of the Pc core opens up the possibility of more in-plane intermolecular interactions between adsorbed molecules on a surface15. Adding an electron withdrawing group to phthalocyanine tunes it remarkably for high performance electronic device applications16. It is an effective way for tailoring different types of self-assemblies and molecular architectures17 with desirable properties for targeted applications in organic electronic and optoelectronic devices18-19. Many molecular

ACS Paragon Plus Environment

3

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 4 of 31

assemblies have directional intermolecular interactions guiding a wide range of one dimensional (1D) and two dimensional (2D) supramolecular self-assemblies20-21 while in other assemblies a balance of the interaction strengths between the intermolecular interaction and molecule-substrate interaction22 forms a bimolecular structure23 making the self-assembly relevant for spintronics24. Molecular self-assembly on 2D materials is also pertinent from the viewpoint of surface functionalization of the 2D substrates. This functionalization, either covalent or non-covalent, has the potential to exploit the inherent properties of the 2D materials and to enhance them for dedicated objectives like flex electronics25. Chiral molecules carry properties of mirror symmetry and handedness. The adsorption and assembly of chiral molecules on a surface leads to symmetry breaking leaving chiral signatures that include molecule-surface interactions26. If both, the molecule and the surface, are achiral the local chirality at a single molecule upon adsorption can lead to supramolecular chirality at the assembly level through the weak, non-covalent bonding interactions27. These interactions are substantially weaker than the covalent bonding within a molecule. These non-covalent interactions can lead to chirality in different forms, e.g. surface induced chirality in a single molecule28 or in a 2D supramolecular assembly29-30. Achiral molecules with a specific geometry can form the self-assembly through intermolecular interactions like hydrogen bonding. Such a self-assembly shows organizational chirality31-33. Also the molecules can encounter conformational changes upon adsorption and exhibit chirality that arises from a change in the conformation of the molecules34-35. Apart from these chiral manifestations, there is one more important kind of chirality which is completely electronic in nature and it has recently been observed in copper phthalocyanine (CuPc) on Ag(100)36 and ZnPc on Cu(100)37. This type of chirality is independent of the geometrical structure of the molecule, but it is energy dependent. It occurs at specific energies, which means it represents the chiral nature of the molecule at a

ACS Paragon Plus Environment

4

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

specific orbital energy. How the pendant group substitution of a phthalocyanine can enhance its intermolecular interactions and introduce energy specific chirality in its self-assembly forms the premise of our study. Herein, using STM and STS at 77 K in UHV, we explore the self-assembly of CuPc(CN)8 on Au(111). We then compare it with the self-assembly of CuPc on Au(111). Both these molecules assemble forming a tetramer unit cell on the surface. They exhibit adsorption induced organizational chirality which arises due to the characteristic arrangement of molecules in the assembly. The herringbone reconstruction of Au(111) remains unperturbed upon CuPc(CN)8 assembly but for the CuPc assembly it increases, which is indicative of a stronger moleculesubstrate interaction than intermolecular interaction for CuPc as compared to CuPc(CN)8. STS measurements pin the HOMO and LUMO energy levels for the molecule in its self-assembled mode. A remarkable observation from STS is that as the chemical structure of the molecule is changed from CuPc to CuPc(CN)8 the HOMO-LUMO gap increases from 3.0 eV for CuPc to 4.46 eV. To verify and to get more insight into these results we have done extensive calculations using density functional theory (DFT). Both our experimental observations are well supported by these calculations. Apart from the organizational chirality, the voltage dependent STM images reveal the chiral nature of individual CuPc(CN)8 molecules on Au(111) only at LUMO energy whereas such chirality is not present in the assembly of CuPc on Au(111). This voltage dependent chirality in the 2D assembly of CuPc(CN)8 on Au(111) is characterized by the orbital energy dependence, thereby making it completely electronic in nature. At LUMO energy, the handedness of each chiral molecule is decided by the direction of the dipolar coupling interaction in the tetramer unit cell. Based on the comparison of CuPc on Au(111) and CuPc(CN)8 on Au(111) we claim that the antiparallel dipolar coupling interaction present in the CuPc(CN)8 tetramer, a consequence of the

ACS Paragon Plus Environment

5

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 6 of 31

cyano substitution of CuPc, strengthens the intermolecular interaction and plays a pivotal role in the chiral manifestation in the CuPc(CN)8 molecular self-assembly at LUMO energy. Thus our work, experimental findings from STM and STS, and supporting DFT calculations, gives us a holistic picture of the system under study, and installs a footprint arising from the cyano substitution of CuPc.

2. METHODS 2.1 Sample Preparation and STM Measurements. The experiments were performed using an Omicron ultra-high vacuum low temperature scanning tunneling microscope (UHV-LTSTM) at 77 K with a base pressure of 5x10-11 mbar. The Au(111) surface was prepared by subsequent Ar+ ion sputtering annealing cycles with 1 keV beam voltage and annealing at 750 K. The surface quality was checked with STM at 77 K. The images revealed herringbone reconstruction and atomic resolution that confirmed an atomically clean surface. CuPc(CN)8 was synthesized using a method in the literature38 and was evaporated from a resistively heated quartz crucible at 650 K onto a clean Au(111) surface held at room temperature. In another experimental run CuPc, bought from Sigma Aldrich and used as is, was deposited at 670 K onto a clean Au(111) surface. See the supplementary information (Figure S1) for the chemical structure of CuPc and CuPc(CN)8. STM imaging and spectroscopy were carried out at 77 K. An electrochemically etched tungsten tip was used for imaging. The images were processed using the image analysis software (SPIP 6.0.9, Image Metrology, Denmark). The local density of states (LDOS) was determined by scanning tunneling spectroscopy (STS). Using standard lock-in technique with a modulation voltage of 20 mV at 652 Hz frequency, the tunneling voltage (V) was varied within ±2.5 V and

ACS Paragon Plus Environment

6

Page 7 of 31 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

the corresponding variation in the current (I) was obtained as a dI/dV spectrum, thus giving the tunneling conductance for the CuPc(CN)8/Au(111) system and CuPc/Au(111) system. 2.2 Computational Details. The calculations were carried out using the spin-polarized density functional theory within the projector augmented wave formalism39 as implemented in VASP40-41. The wave functions were expanded in plane wave basis with 500 eV cutoff. The exchange-correlation energy were described with Perdew-Burke-Ernzerhof (PBE) form of generalized gradient approximation42 which very well describes the structural and electronic properties of bulk-Au and CuPc and CuPc(CN)8 molecules in the gas phase. To better model the interaction between the molecules and the Au(111) substrate, the van der Waals interaction was incorporated using DFT-D3 formalism with zero damping43. The Au substrate included four Au layers with 10×9 lateral supercell, which were separated by 18 Å vacuum. While the first two layers were fixed to the bulk, the top layers and the molecules were fully relaxed until the forces were smaller than 0.01 eV/Å. The Brillouin zone was sampled using 4×4×1 Monkhorst-Pack 𝑘grid. A 30×30×12 Å3 supercell was used to study the molecules in the gas phase, and the reciprocal space integration was carried using the Γ point. The molecular adsorption energy was calculated as 𝐸ads = 𝐸molecule/Au(111) − 𝐸Au(111) − 𝐸molecule for CuPc/Au(111) and CuPc(CN)8/Au(111).

3. RESULTS AND DISCUSSION CuPc(CN)8 molecules38 were deposited on the Au(111) surface. These molecules assemble in the form of hundred nanometer-sized rectangular clusters scattered over the Au(111) surface and in some cases the assembly spills over the step edges. Figure 1a shows an STM image of one cluster of CuPc(CN)8 where these square-planar molecules are assembled in a rectangular array. Through the array and in the background of the array the bright lines are the lines of the herringbone

ACS Paragon Plus Environment

7

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 8 of 31

reconstruction. The white trace on top of the cluster shows the line profile of the reconstruction on the surface and over the array. This reconstruction is attributed to the atoms in the topmost layer being compressed by 4.5% in the [11̅0] direction thereby placing 23 atoms over the 22 atoms in the underneath layer. To remain thermodynamically stable, this compression leads to the formation of a periodic array of cubic close packing face centered cubic (fcc) and hexagonal close packing (hcp) regions. The interface between the hcp and fcc regions is marked by the dotted line, Figure 1a. These regions are separated by a region where atoms go through a bridge-site transition that appears as a bright line in the STM image called “soliton wall”44. A pair of soliton walls is called “herringbone” and their average separation is 6.34 nm on Au(111). This separation, also referred to as herringbone reconstruction periodicity, is a sensitive tool to measure the interaction strength of an absorbed molecule upon a substrate. For a strong interaction between the adsorbed molecules on Au(111) the herringbone separation increases from 6.34 nm and disappears when all the 4.5% Au atoms are removed depending on the strength of the interaction45-48. For the CuPc(CN)8 assembly the herringbone spacing is 6.30 nm (Figure 1a) which is almost the same as its free space value. Hence, the self-assembly does not remove the reconstruction, indicating a weak van der Waals bonding between the CuPc(CN)8 and the Au(111) surface.

ACS Paragon Plus Environment

8

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

Figure 1. (a) STM image of CuPc(CN)8 clusters on Au(111) substrate with herringbone reconstruction in background (b) An enlarged cluster showing molecular self-assembly, referred to as configuration α, with the inset showing a ball-and-stick model of CuPc(CN)8 superimposed on a single molecule. The unit vectors 𝑎⃗ and 𝑏⃗⃗ are shown to specify the unit cell in the assembly. (c) Enlarged cluster with molecular assembly referred to as configuration β, with the unit vectors 𝑎⃗ and 𝑏⃗⃗ showing the unit cell in the assembly. (d) STM image of CuPc showing a monolayer coverage on Au(111) with herringbone reconstruction in background. All STM images are taken at 0.2 nA and 1.0 V. (e) and (f) STM images represent configurations α and β for CuPc. (g) and (h) are the schematic of the tetrameric unit cell for the α, and β configurations for CuPc and CuPc(CN)8 respectively and δ is the angle between the unit cell vector of the assembly and the molecular symmetry axis M1. (a) and (d) scale bar = 10 nm and (b), (c), (e), (f) scale bar = 2 nm. The 4-fold symmetry of CuPc(CN)8 is observed in the self-assembly as a 2D square lattice. Figures 1b and 1c are the close-up STM images of the assembly where the 2D lattice has unit cell parameters 𝑎⃗ = 𝑏⃗⃗ = 1.56 𝑛𝑚 and the molecules have a high packing density of 0.41

ACS Paragon Plus Environment

9

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 10 of 31

molecule/nm2. Figure 1b corresponds to the α configuration and Figure 1c shows the β configuration. Both these are mirror symmetric to each other. This assembly of CuPc(CN)8 over Au(111) shows the adsorption induced organizational chirality. Figure 1g is the schematic for the tetrameric unit cell of CuPc(CN)8. The configuration is referred to as α when the angle δ between the unit cell vector of the assembly and the molecular symmetry axis M1 of CuPc(CN)8 is 32±2° and β when it is 56±2°. In both configurations the unit vectors are 𝑎⃗ = 𝑏⃗⃗ = 1.56 𝑛𝑚. The configurations α and β are unique to each molecular cluster, no mixed configuration was seen in any cluster. The black dotted lines in Figure 1g denote the probable intermolecular H…NC hydrogen bonding between CuPc(CN)8 molecules that mediates the intermolecular attraction and gives rise to a 2D close-packed assembly on Au(111), and guides the pattern formation. This is similar to the self-assembly of ZnPcCl8 where the 2D hydrogen bonded network is formed due to the H•••Cl interaction46. In this case, it is the highly electronegative nitrogen belonging to the cyano group that interacts with two hydrogen atoms of the nearest neighboring molecule. In addition to the hydrogen bonding there is an intrinsic dipole due to each cyano group that also interacts with each other in the molecular assembly. This interaction is the anti-parallel dipolar coupling. This interaction was also observed in the self-assembly of cyano-functionalized porphyrins20 and triarylamines49 on Au(111). Both the hydrogen bonding and dipolar coupling intermolecular interactions stabilize the self-assembly giving rise to submonolayer domains of different sizes and shapes. No single CuPc(CN)8 molecule adsorption was observed on the Au(111) surface even at low coverage as the substrate was held at room temperature during molecule deposition but observing single molecule at very low temperature like 5K is plausible. Only 2D regions with tetramers were observed. Also conspicuous by its absence was the formation

ACS Paragon Plus Environment

10

Page 11 of 31 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

of 1D chain. This could be due to the high affinity for hydrogen bonding in all the four in-plane directions. The self-assembly of CuPc is shown in Figure 1d. These molecules form a monolayer coverage (unlike clusters for CuPc(CN)8) going over the step edges of Au(111). The assembly is planar and each molecule has a square geometry as expected. The herringbone reconstruction can be seen through the monolayer. The herringbone reconstruction periodicity for this assembly is 7.2±0.3 nm, larger than the unperturbed 6.3 nm that is seen for CuPc(CN)8 as discussed above. This indicates that, unlike CuPc(CN)8, there is a stronger bonding of CuPc with Au(111). This enhanced periodicity is a direct pointer to the stronger molecule-substrate interaction in the CuPc/Au(111) assembly. There are 2 different configurations of the molecule within the assembly, α shown in the STM image in Figure 1e, and β in the STM image in Figure 1f. These configurations are enantiomers and show their chiral nature on the surface only upon adsorption as neither Au(111) nor CuPc is chiral. Hence the chirality in the self-assembly is induced upon adherence to the surface. The schematic of the tetrameric unit cell for CuPc self-assembly is shown in Figure 1g for configuration α, and in Figure 1h for configuration β. The CuPc configuration is referred to as α (Figure 1e) when the angle δ between the unit cell vector of the assembly and the molecular symmetry axis M1 of CuPc is 32±2° and β (Figure 1f) when it is 59±2°. In both configurations the unit vectors are 𝑎⃗ = 𝑏⃗⃗ = 1.53 ± 0.1 𝑛𝑚 and the packing density is 0.43 molecule/nm2. The black dotted lines for the CuPc unit cell denote the probable intermolecular van der Waals interactions between the CuPc molecules that bind the assembly in addition to the stronger molecule-substrate interaction evidenced by the 7.2 ± 0.3 nm herringbone reconstruction periodicity that makes the CuPc/Au(111) self-assembly stable.

ACS Paragon Plus Environment

11

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 12 of 31

During STS measurements, the tunneling of electrons through the HOMO and LUMO gives rise to peaks in dI/dV. The dI/dV spectra were taken at the center of the molecule, for both, CuPc(CN)8 and CuPc, and were measured for several molecules located in different positions on the substrate. The average of dI/dV spectra taken over CuPc(CN)8 molecule (orange curve), CuPc molecule (green curve), and taken on the bare Au(111) substrate(black curve) are shown in Figure 2. As seen from Figure 2, for CuPc(CN)8 the LUMO occurs at 2.0 V and the HOMO occurs at -2.3 V with reference to the Fermi energy at V=0. For CuPc the LUMO occurs at 1.5 V and the HOMO occurs at -1.5 V. The key result of these spectra is that the HOMO-LUMO gap increases from 3.0 V to 4.3 V as we go from CuPc to CuPc(CN)8. For CuPc(CN)8 the STS curve measured over different molecules in position of the cluster but no significant change observed with respect to the STS position (see supplementary Figure S4).Since the spectra were recorded at different sites on Au(111) any substrate site-dependent orbital energy fluctuation, as seen in other systems, can be ruled out.50 Spin density functional theory (SDFT) calculations have shown that the band gap increases from CuPc (1.1 eV) to F16CuPc (1.4 eV).51 For our system, we speculate that the increase in the HOMO-LUMO gap from CuPc to CuPc(CN)8 is a consequence of the role played by the cyano groups. The HOMO-LUMO gap is also known to be affected by the molecule-substrate interaction. As the molecule-substrate interaction decreases the HOMO-LUMO gap increases. This is due to the reduced screening of the substrate that enhances the Coulomb repulsion for the tunneling electrons52-53. Our observation of the increase of HOMO-LUMO gap for the CuPc(CN)8/Au(111) system, where the molecule-substrate interaction is weaker compared to the CuPc/Au(111) system, as evidenced by the intact herringbone periodicity and by the presence of surface state54 in the STS data for the CuPc(CN)8/Au(111) system, is consistent and in line with the expected behavior of such molecule-substrate systems.

ACS Paragon Plus Environment

12

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

Figure 2. dI/dV spectra over the Cu atom for CuPc(CN)8 (orange) and CuPc (green). The HOMO and LUMO peaks for CuPc(CN)8 at -2.3 V and 2.0 V and the same for CuPc are at -1.5 V and 1.5 V. The Au(111) surface state at -0.5 V is indicated by the black trace in the dI/dV spectra over Au(111)(black). The blue dotted oval highlights the absence of the Au(111) surface state in CuPc gap, the presence and a slight shift of the Au(111) surface state towards the Fermi level in the CuPc(CN)8 gap and the surface state of bare Au(111) surface. The self-assembly of CuPc is stronger than that of CuPc(CN)8 as seen from the change in the periodicity of the herringbone reconstruction of Au(111) through the assembly. This is also reflected in the dI/dV spectrum of CuPc (green trace in Figure 2) where the Au(111) surface state is not seen. However, the surface state can be seen in the dI/dV spectrum of CuPc(CN) 8, (orange trace in Figure 2). Thus our inference of the molecule-substrate interaction for the CuPc/Au(111) system being stronger than the CuPc(CN)8/Au(111) system is validated by the change in the herringbone reconstruction periodicity and the quenching of the surface state for the CuPc/Au(111) system and the trend of the HOMO-LUMO gap variation. The slight shift of the Au(111) surface

ACS Paragon Plus Environment

13

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 14 of 31

state towards the Fermi level in case of CuPc(CN)8 can be ascribed to the phenomenon of molecular adsorption upon Au(111) albeit with a weak van der Waals interaction55. To corroborate the complex structural modulation of Au(111) substrate and the electronic interaction of CuPc and CuPc(CN)8 molecules with the substrate, we have performed extensive ab initio calculations within density functional theory implemented in VASP40-41. An unreconstructed Au(111) slab model is a good approximation for studying the electronic properties single CuPc and CuPc(CN)8 molecules along with the structural changes of the substrate that may occur due molecular adsorption. There are four high symmetry absorption sites on Au(111) surface namely the top, hcp-hollow, fcc-hollow, and bridge sites. However, the orientation of the absorbed molecules is complex due to their rotational degree of freedom with respect to the crystalline direction of the substrate. We find that both CuPc and CuPc(CN)8 molecules are adsorbed at the hcp-hollow site. The most stable azimuthal orientation of the Cu-ligand axis for both molecules is rotated by ~30O from the [11̅0] direction of the Au(111) surface [Figures 3a and 3b]. Such azimuthal rotation increases the binding energy by 0.25 and 0.06 eV for CuPc and CuPc(CN)8 respectively, from the corresponding hcp-hollow site absorption. The CN-group in R-CuPc(CN)8 strongly influences the interaction with the substrate and binds strongly (-5.86 eV) compared to the R-CuPc, which is -4.66 eV. The trend in the binding energy can be explained in terms of the interaction between the sp2 orbitals localized on aza-N atoms with the Au atoms below36. The average vertical distance between aza-N atoms for R-CuPc(CN)8 with the Au(111) substrate is found to be much shorter 3.11 Å compared to that for the R-CuPc, which is 3.23 Å. Thus, RCuPc(CN)8 binds strongly through the increased interaction between aza-N atoms and Ausubstrate. Consequently, in their optimized configurations a remarkable net charge transfer of 0.55 and 1.22 electrons into the CuPc and CuPc(CN)8 molecules, respectively, obtained from the Bader

ACS Paragon Plus Environment

14

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

charge analysis56 explain the trend in adsorption energy. The experimental trend observed in the herringbone reconstruction can be qualitatively explained by investigating the structural details of the top Au-layer without and with absorbed molecule. Within the PBE exchange-correlation functional56 the distance between the successive Au(111) layers for bulk-Au is 𝑎⁄√3 ~ 2.40 Å, where 𝑎 is the lattice parameter 4.16 Å. Due to the lack of periodicity, the top layer of the bare Au(111) surface is reconstructed such that the vertical distance of the top layer with the layer below decreases to 2.37 Å, whereas the layers below are unperturbed. While the adsorption of RCuPc(CN)8 does not perturb the picture above, the R-CuPc molecule pulls up the Au-atoms below by an amount of 0.01 Å. This explains the observation of modified herringbone reconstruction periodicity in STM images for CuPc molecular assemblies. Thus, while the periodicity of the herringbone reconstruction is unperturbed for R-CuPc(CN)8 absorption, remaining intact at 6.3 Å, for R-CuPc this periodicity increases to 7.2 ± 0.3 Å. The calculated spin-polarized density of states (DOS) indicates that the magnetic moment of 1 and 0.92 𝜇𝐵 for CuPc and CuPc(CN)8, respectively, arises from the singly occupied localized Cu𝑑𝑥 2 −𝑦 2 orbital, which lies at -0.5 and -0.85 eV with respect to the Fermi level, Figure 3c. This is in contrast to their gas phase counterparts, where Cu-𝑑𝑥 2 −𝑦 2 lies at ~ -0.45 eV for both molecules, and thus the cyano group does not affect the Cu-𝑑 orbitals in the gas phase. Note that 𝑑𝑥 2 −𝑦 2 orbital does not contribute to the tunneling current, and the DOS indicates that the HOMO is of Cu-𝑑𝑧 2 and exposed to the STM-tip. Similar to 𝑑𝑥 2 −𝑦 2 orbital, in the gas phase the cyano group does not affect the energy position of the 𝑑𝑧 2 , which lies at -2.51 and -2.48 eV for CuPc and CuPc(CN)8 respectively. However, the picture changes drastically while the molecules are adsorbed on the Au(111) substrate. Owing to the presence of the cyano group, the Cu-𝑑𝑧 2 HOMO is pushed deeper in energy at -2.85 eV due to stronger aza-N and Au interaction for

ACS Paragon Plus Environment

15

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 16 of 31

CuPc(CN)8/Au(111). However, this is unperturbed for R-CuPc/Au(111). This picture qualitatively explains the experimental trend in HOMO. The degenerate 𝑑𝑥𝑧 and 𝑑𝑦𝑧 for both molecules lie deeper in energy compared to the corresponding 𝑑𝑧 2 orbital.

Figure 3. (a) Configuration of the adsorbed molecules and corresponding density of states. Both (a) CuPc and (b) CuPc(CN)8 are adsorbed at the high-symmetry hcp-hollow site on Au(111) substrate with a ~300 azimuthal rotation of the Cu-ligand about the [11̅0] direction of the substrate. (c) The spin-polarized density of states projected on Cu, calculated with 0.1 eV smearing and compared with the corresponding molecule in the gas phase. The Cu-𝑑𝑥 2 −𝑦 2 and Cu-𝑑𝑧 2 are shown. While the singly occupied planer 𝑑𝑥 2 −𝑦 2 orbital contributes to the magnetic moment, the 𝑑𝑧 2 orbital interacts with the STM-tip and generates the tunneling current. In accordance with the

ACS Paragon Plus Environment

16

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

experimental results, the cyano group and concurrent adsorption perturbs the 𝑑𝑧 2 HOMO and pushes it away from the Fermi level 𝐸𝐹 . With HOMO and LUMO energy levels as reference points voltage dependent imaging was carried out for both CuPc(CN)8 and CuPc assemblies to see if the conformation of the selfassembled molecules undergoes any change as a function of the voltage. Figure 4 shows selective images related to this study. Note that the adsorption induced organizational chirality is present at all the voltages. For CuPc(CN)8 on Au(111) no chiral contrast is seen in the HOMO-LUMO energy gap and at HOMO energy. This is clear from the STM images at -2.2 V (Figures 4a and 4b) where the conformation remains unchanged. Apart from the organizational chirality an unexpected chiral conformational change is observed at 2.2 V, the LUMO energy, when the individual molecule in the tetramer conformation becomes chiral, Figure 4c. As chirality of individual molecule is dependent on the specific energy so it is electronic nature. The molecular assembly undergoes a corresponding change that is enhanced with increasing positive bias voltage. This electronic chirality emerges at the onset of LUMO energies from 1.8 V and remains chiral up to applied bias of 2.5 V (Figure S2). This energy range matches with the onset of the LUMO peak in the dI/dV spectra, Figure 2. The underlying organizational chirality remains intact and this energy specific electronic chirality in individual molecule emerges on top of the organizational chirality. In this assembly individual molecule at LUMO energy becomes chiral so the chirality is propagating from the organizational chirality to electronic chirality in individual molecule. This is evident from Figure 4 and is discussed in detail further ahead. For CuPc on Au(111), at HOMO energy, there is a slight change in the appearance of the individual molecules in the assembly, Figure 4d. Each lobe of the 4-lobed molecule looks squareish as opposed to the oval shape of the lobes. The oval shape is commonly seen at any energy

ACS Paragon Plus Environment

17

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 18 of 31

in the gap, Figure 4e, and at LUMO energy, Figure 4f. Thus no change in the shape of the lobe occurs at any energy in the gap or at LUMO energy. The slight deviation of the lobe shape at HOMO (-1.5 V) (which starts at -1.2 V, see Figure S3), could be due to the interaction of the CuPc HOMO orbital with the Au(111) surface as the CuPc/Au(111) interaction is strong. A bias dependent change in the molecular structure has been reported for the BTDA-TCNQ/Au(111) supramolecular assembly at LUMO energy due to the unoccupied states at LUMO energy57. The emergence of supramolecular chirality in CuPc(CN)8 at LUMO energy is a striking observation and the absence of chirality in CuPc assembly is a pointer to the role played by the 8 cyano groups in contributing to this orbital dependent electronic chirality that shows up only at the specific LUMO energy. Such a chiral behavior has been reported for CuPc/Ag(100) for single CuPc molecules and in the CuPc self-assembly as well. The asymmetric charge transfer between the molecules and the substrate is reported to induce this chirality36.

Figure 4. (a-c) are the voltage dependent STM images of CuPc(CN)8 on Au(111) at voltages corresponding to HOMO , in the gap, and LUMO for the same area taken with constant current 0.2 nA. Among them (c) at 2.2 V shows chirality. (d-f) are the voltage dependent STM images of

ACS Paragon Plus Environment

18

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

CuPc on Au(111) at voltages corresponding to HOMO, in the gap and LUMO for the same area taken with constant current 0.2 nA. No chirality is seen for CuPc assembly at any voltage. Scale bar = 2 nm for all images.

As the supramolecular self-assembly of CuPc(CN)8 on Au(111) becomes chiral at LUMO energy, the handedness emerges. The enantiomeric configurations α and β exhibit left and right handed chirality at LUMO energy. Shown in Figures 5a and 5b are the STM images taken with bias voltage 2.2 V, the LUMO energy. The schematic tetrameric unit cells are shown in Figures 5c and 5d. The organizational chirality, which is seen independent of bias voltages, as discussed after Figure 1, is present in the images in Figures 5a and 5b too. But the emergence of handedness at LUMO energy is an unprecedented observation. This is explained using the ball-and-stick model of CuPc(CN)8 shown in Figures 5c and 5d. For the α configuration, the green oval in the ball and stick model of Figure 5c is a schematic representation of the antiparallel dipolar coupling between the nearest neighbor cyano groups. This coupling is shown by the white dotted oval in the STM image, Figure 5a, of the α configuration. This explains the left handed chirality seen in the STM image. Similarly, for β configuration, the green oval in the ball-and-stick model of Figure 5d representing the antiparallel dipolar coupling between the nearest neighbor –CN groups is reflected in the white dotted oval in the STM image, Figure 5b, of the β configuration. This explains the right handed chirality as observed in the STM image.

ACS Paragon Plus Environment

19

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 20 of 31

Figure 5. STM images of CuPc(CN)8 on Au(111). Images a and b (10 nm X 10 nm) are taken at the LUMO energy of 2.2 V, I = 0.2 nA. They are labeled ‘LEFT’ and ‘RIGHT’ due to their chirality; (c) and (d) illustrate the tetrameric assembly. The chiral directions are indicated by the red solid arrows. Black dotted lines denote hydrogen bonding. Green solid ovals represent the antiparallel dipolar coupling that gives higher contrast in STM images at LUMO energy. In a and b the white dotted ovals represent the corresponding green oval in illustrations (c) and (d). The STM images, Figures 5a and 5b reveal that the region between the two lobes involved in the antiparallel dipolar coupling becomes bright. In both configurations each molecule has four such antiparallel dipolar coupling interactions involved with each lobe. For configuration α, this occurs on the right of the molecular axis (red dotted arrow in Figure 5c) leading to left handedness; and for configuration β, the four antiparallel dipolar coupling interactions involved with each lobe occur on the left of the molecular axis (red dotted arrow in Figure 5d) giving right handedness for the β configuration. The organizational chirality uniquely defines the direction of the antiparallel

ACS Paragon Plus Environment

20

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

dipolar coupling. This coupling provides the handedness in the electronic chirality in the individual molecule in each of the configurations α and β at LUMO energy. Thus the electronic and organizational chirality are coupled and are irreversible. Thus, we are able to comment from the sub-molecular resolution of the STM images of the tetrameric unit cell of the CuPc(CN) 8 selfassembly on Au(111), that the antiparallel dipolar coupling seems to be the driving force for the observed chirality at LUMO energy in the CuPc(CN)8-Au(111) system. We emphasize that the chirality is “turned on” only at the LUMO energy and supports our observation of the weak interaction of the CuPc(CN)8 with the Au(111) substrate thereby underscoring the lack of contribution from the substrate to this unique observation of left handed and right handed chiral assemblies. The intermolecular interaction in the CuPc(CN)8 assembly is augmented by the cyano groups, thus making the cyano substitution of CuPc a clear contributor to this chiral behavior.

4. CONCLUSION In summary, we have presented an atomic scale investigation of the influence of the cyano substitution of CuPc on its molecular self-assembly. STM images of CuPc(CN)8 on Au(111) showed a self-assembly with a tetrameric unit cell that exhibited adsorption induced structural chirality and kept the herringbone reconstruction periodicity unperturbed. Monolayer coverage of CuPc on Au(111) too showed tetramers in two mirror symmetric configurations. But the herringbone reconstruction periodicity was modified. This is indicative of a stronger moleculesubstrate interaction between the CuPc/Au(111) assembly than the CuPc(CN)8/Au(111) assembly. From the STS measurements the HOMO-LUMO gap showed a drastic increment from 3.0 eV for CuPc/Au(111) system to 4.3 eV for the CuPc(CN)8/Au(111) system. We attribute these results to the partial cyanation of CuPc. These experimental findings of herringbone periodicity modification

ACS Paragon Plus Environment

21

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 22 of 31

and the variation of the HOMO-LUMO gap were validated by ab initio DFT calculations. Voltage dependent STM images unveil the orbital energy specific chirality of the CuPc(CN)8 molecules and superimpose right-handedness and left-handedness on the existing structural chirality within the assembly. However, for CuPc no such energy dependent chiral behavior was witnessed at any orbital energy. Our work thus presents a strategy to functionalize molecules for a preferential interaction among them or with the underlying substrate, thereby enabling surface functionalization and templating for applications in catalysis and gas sensing. ASSOCIATED CONTENT Supporting Information: This material is available free of charge via the Internet at http://pubs.acs.org. A ball-and-stick model of CuPc and CuPc(CN)8 molecules highlighting that the -CN group replaces 8 outer H atoms in CuPc, a series of voltage dependent STM images of CuPc(CN)8 on Au(111) to show the onset of chirality at the LUMO energy of CuPc(CN)8, a series of voltage dependent STM images of CuPc on Au(111) for comparison to show that there is no voltage dependent chiral behaviour in this case, STS Curves measured at different positions of the cluster of CuPc(CN)8 molecules on Au(111).

AUTHOR INFORMATION Corresponding Author *E-mail: [email protected]. Phone: +91(0)20 25908063 Notes The authors declare no competing financial interests.

ACS Paragon Plus Environment

22

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

ACKNOWLEDGMENTS AD thanks Dr. Nirmalya Ballav for insightful discussions. We acknowledge the supercomputing facilities at the Centre for Development of Advanced Computing, Pune; Inter University Accelerator Centre, Delhi; and at the Center for Computational Materials Science, Institute of Materials Research, Tohoku University. Funding sources: This work was financially supported by the Indian Institute of Science Education

and

Research

(IISER),

Pune

and DAE-BRNS

(India,

Project

No.

2011/20/37C/17/BRNS). M. K. acknowledges the funding from the Department of Science and Technology, Government of India under Ramanujan Fellowship, and Nano Mission project SR/NM/TP-13/2016. SA and BD is grateful to CSIR Government of India for financial support in the form of research fellowship.

REFERENCES 1.

Patra, N.; Král, P.; Sadeghpour, H. R. Nucleation and Stabilization of Carbon-Rich

Structures in Interstellar Media. The Astrophysical Journal 2014, 785, 6. 2.

Dworkin, J. P.; Deamer, D. W.; Sandford, S. A.; Allamandola, L. J. Self-Assembling

Amphiphilic Molecules: Synthesis in Simulated Interstellar/Precometary Ices. Proc. Natl. Acad. Scie. U.S.A. 2001, 98, 815-819. 3.

Naughton, K. L.; Phan, L.; Leung, E. M.; Kautz, R.; Lin, Q.; Van Dyke, Y. Marmiroli, B.;

Sartori, B.; Arvai, A.; Li, S.; et al. Self-Assembly of the Cephalopod Protein Reflectin. Adv. Mater. 2016, 28, 8405-8412.

ACS Paragon Plus Environment

23

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

4.

Page 24 of 31

Hynninen, A.-P.; Thijssen, J. H. J.; Vermolen, E. C. M.; Dijkstra, M.; van Blaaderen, A.

Self-Assembly Route for Photonic Crystals with a Bandgap in the Visible Region. Nat. Mater. 2007, 6, 202-205. 5.

Evans, R. C. Harnessing. Self-Assembly Strategies for the Rational Design of Conjugated

Polymer Based Materials. J. Mater. Chem. C 2013, 1, 4190-4200. 6.

Yeow, J.; Sugita, O. R.; Boyer, C. Visible Light-Mediated Polymerization-Induced Self-

Assembly in the Absence of External Catalyst or Initiator. ACS Macro Lett. 2016, 5, 558-564. 7.

Deshpande, A.; Sham, C.-H.; Alaboson, J. M. P.; Mullin, J. M.; Schatz, G. C.; Hersam, M.

C. Self-Assembly and Photopolymerization of Sub-2 Nm One-Dimensional Organic Nanostructures on Graphene. J. Am. Chem. Soc. 2012, 134, 16759-16764. 8.

Lin, J.; Yu, Z.; Zhu, W.; Xing, G.; Lin, Z.; Yang, S.; Xie, L.; Niu, C.; Huang, W. A π-

Conjugated Polymer Gelator from Polyfluorene-Based Poly(Tertiary Alcohol) via the HydrogenBonded Supramolecular Functionalization. Polym. Chem. 2013, 4, 477-483. 9.

Lin, J. Y.; Liu, B.; Yu, M. N.; Ou, C. J.; Lei, Z. F.; Liu, F.; Wang, X. H.; Xie, L.; Zhu, W.

S.; Ling, H.; et al. Understanding the Molecular Gelation Processes of Heteroatomic Conjugated Polymers for Stable Blue Polymer Light-Emitting Diodes. J. Mater. Chem. C 2017, 5, 6762-6770. 10.

Yu, M.; Liu, B.; Lin, J.; Bai, L.; Ling, H.; Zhu, W.; Xie, L.; Yi, M.; Wang, J.; Huang, W.

One-Step Preparation of Conjugated Homopolymer Sub-Microspheres via a Controllable Supramolecular Approach toward Optoelectronic Applications. RSC Adv. 2017, 7, 14688-14693. 11.

Lin, J.-Y.; Wong, J.; Xie, L.-H.; Dong, X.-C.; Yang, H. Y.; Huang, W. Hydrogen-Bonded

Supramolecular Conjugated Polymer Nanoparticles for White Light-Emitting Devices. Macromol. Rapid Commun. 2014, 35, 895-900.

ACS Paragon Plus Environment

24

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

12.

Mali, K. S.; Pearce, N.; De Feyter, S.; Champness, N. R. Frontiers of Supramolecular

Chemistry at Solid Surfaces. Chem. Soc. Rev. 2017, 46, 2520-2542. 13.

Ragoussi, M.-E.; Cid, J.-J.; Yum, J.-H.; de la Torre, G.; Di Censo, D.; Grätzel, M.;

Nazeeruddin, M. K.; Torres, T. Carboxyethynyl Anchoring Ligands: A Means to Improving the Efficiency of Phthalocyanine-Sensitized Solar Cells. Angew. Chem. Int. Ed. 2012, 51, 4375-4378. 14.

Alippi, P.; Lanzilotto, V.; Paoletti, A. M.; Mattioli, G.; Zanotti, G.; Pennesi, G.; Filippone,

F.; Cossaro, A.; Verdini, A.; Morgante, A. et al. A Ru-Ru Pair Housed in Ruthenium Phthalocyanine: The Role of a "Cage" Architecture in the Molecule Coupling with the Ag(111) Surface. Phys. Chem. Chem. Phys. 2017, 19, 1449-1457. 15.

Jiang, J. Functional Phthalocyanine Molecular Materials; Springer 2010; Vol. 135.

16.

Bao, Z.; Lovinger, A. J.; Brown, J. New Air-Stable N-Channel Organic Thin Film

Transistors. J. Am. Chem. Soc. 1998, 120, 207-208. 17.

Barrena, E.; de Oteyza, D. G.; Dosch, H.; Wakayama, Y. 2d Supramolecular Self-

Assembly of Binary Organic Monolayers. ChemPhysChem 2007, 8, 1915-1918. 18.

Babudri, F.; Farinola, G. M.; Naso, F.; Ragni, R. Fluorinated Organic Materials for

Electronic and Optoelectronic Applications: The Role of the Fluorine Atom. Chem. Commun. 2007, 1003-1022. 19.

Brinkmann, H.; Kelting, C.; Makarov, S.; Tsaryova, O.; Schnurpfeil, G.; Wöhrle, D.;

Schlettwein, D. Fluorinated Phthalocyanines as Molecular Semiconductor Thin Films. phys. stat. sol. (a) 2008, 205, 409-420. 20.

Yokoyama, T.; Yokoyama, S.; Kamikado, T.; Okuno, Y.; Mashiko, S. Selective Assembly

on a Surface of Supramolecular Aggregates with Controlled Size and Shape. Nature 2001, 413, 619-621.

ACS Paragon Plus Environment

25

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

21.

Page 26 of 31

Müller , K.; Moreno-López, J. C.; Gottardi, S.; Meinhardt, U.; Yildirim, H.; Kara, A.;

Kivala, M.; Stöhr, M. Cyano-Functionalized Triarylamines on Coinage Metal Surfaces: Interplay of Intermolecular and Molecule–Substrate Interactions. Chem. – Eur. J. 2016, 22, 581-589. 22.

de Oteyza, D. G.; Silanes, I.; Ruiz-Osés, M.; Barrena, E.; Doyle, B. P.; Arnau, A.; Dosch,

H.; Wakayama, Y.; Ortega, J. E. Balancing Intermolecular and Molecule–Substrate Interactions in Supramolecular Assemblies. Adv. Funct. Mater. 2009, 19, 259-264. 23.

Hipps, K. W.; Scudiero, L.; Barlow, D. E.; Cooke, M. P. A Self-Organized 2-Dimensional

Bifunctional Structure Formed by Supramolecular Design. J. Am. Chem. Soc. 2002, 124, 21262127. 24.

Wäckerlin, C.; Nowakowski, J. ; Liu, S. X.; Jaggi, M.; Siewert, D.; Girovsky, J.; Shchyrba,

A.; Hählen, T.; Kleibert, A.; Oppeneer, P. M.; et al. Two-Dimensional Supramolecular Electron Spin Arrays. Adv. Mater. 2013, 25, 2404-2408. 25.

Avijit, K.; Kaustuv, B.; Peter, L. Molecular Assembly on Two-Dimensional Materials.

Nanotechnology 2017, 28, 082001. 26.

Ernst, K.-H. Molecular Chirality at Surfaces. physica status solidi (b) 2012, 249, 2057-

2088. 27.

Raval, R. Chiral Expression from Molecular Assemblies at Metal Surfaces: Insights from

Surface Science Techniques. Chem. Soc. Rev. 2009, 38, 707-721. 28.

Lopinski, G. P.; Moffatt, D. J.; Wayner, D. D. M.; Wolkow, R. A. Determination of the

Absolute Chirality of Individual Adsorbed Molecules Using the Scanning Tunnelling Microscope. Nature 1998, 392, 909-911. 29.

Ortega Lorenzo, M.; Baddeley, C. J.; Muryn, C.; Raval, R. Extended Surface Chirality

from Supramolecular Assemblies of Adsorbed Chiral Molecules. Nature 2000, 404, 376-379.

ACS Paragon Plus Environment

26

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

30.

Charra, F.; Cousty, J. Surface-Induced Chirality in a Self-Assembled Monolayer of

Discotic Liquid Crystal. Phys. Rev. Lett. 1998, 80, 1682-1685. 31.

Liriano, M. L.; Carrasco, J.; Lewis, E. A.; Murphy, C. J.; Lawton, T. J.; Marcinkowski, M.

D.; Therrien, A. J.; Michaelides, A.; Sykes, E. C. H. The Interplay of Covalency, Hydrogen Bonding, and Dispersion Leads to a Long Range Chiral Network: The Example of 2-Butanol. J. Chem. Phys. 2016, 144, 094703. 32.

Böhringer, M.; Morgenstern, K.; Schneider, W.-D.; Berndt, R.; Mauri, F.; De Vita, A.; Car,

R. Two-Dimensional Self-Assembly of Supramolecular Clusters and Chains. Phys. Rev. Lett. 1999, 83, 324-327. 33.

Weckesser, J.; De Vita, A.; Barth, J.; Cai, C.; Kern, K. Mesoscopic Correlation of

Supramolecular Chirality in One-Dimensional Hydrogen-Bonded Assemblies. Phys. Rev. Lett. 2001, 87, 096101. 34.

Weigelt, S.; Busse, C.; Petersen, L.; Rauls, E.; Hammer, B.; Gothelf, K. V.; Besenbacher,

F.; Linderoth, T. R. Chiral Switching by Spontaneous Conformational Change in Adsorbed Organic Molecules. Nat. Mater. 2006, 5, 112-117. 35.

Booth, N. A.; Woodruff, D. P.; Schaff, O.; Gießel, T.; Lindsay, R.; Baumgärtel, P.;

Bradshaw, A. M. Determination of the Local Structure of Glycine Adsorbed on Cu(110). Surf. Sci. 1998, 397, 258-269. 36.

Mugarza, A.; Lorente, N.; Ordejón, P.; Krull, C.; Stepanow, S.; Bocquet, M.-L.; Fraxedas,

J.; Ceballos, G.; Gambardella, P. Orbital Specific Chirality and Homochiral Self-Assembly of Achiral Molecules Induced by Charge Transfer and Spontaneous Symmetry Breaking. Phys. Rev. Lett. 2010, 105, 115702.

ACS Paragon Plus Environment

27

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

37.

Page 28 of 31

Chen, F.; Chen, X.; Liu, L.; Song, X.; Liu, S.; Liu, J.; Ouyang, H.; Cai, Y.; Liu, X.; Pan,

H. Chiral Recognition of Zinc Phthalocyanine on Cu (100) Surface. Appl. Phys.Lett. 2012, 100, 081602. 38.

Wöhrle, D.; Meyer, G.; Wahl, B. Polymere Phthalocyanine Und Ihre Vorstufen, 1.

Reaktive Oktafunktionelle Phthalocyanine Aus 1,2,4,5-Tetracyanbenzol. Makromol.Chem. 1980, 181, 2127-2135. 39.

Blöchl, P. E. Projector Augmented-Wave Method. Phys. Rev. B 1994, 50, 17953-17979.

40.

Kresse, G.; Hafner, J. Ab Initio. Phys. Rev. B 1993, 47, 558-561.

41.

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

Perdew, J. P.; Burke, K.; Ernzerhof, M. Generalized Gradient Approximation Made

Simple. Phys. Rev. Lett. 1996, 77, 3865-3868. 43.

Grimme, S.; Antony, J.; Ehrlich, S.; Krieg, H. A Consistent and Accurate Ab Initio

Parametrization of Density Functional Dispersion Correction (Dft-D) for the 94 Elements H-Pu. J. Chem. Phys. 2010, 132, 154104. 44.

Barth, J.; Brune, H.; Ertl, G.; Behm, R. Scanning Tunneling Microscopy Observations on

the Reconstructed Au (111) Surface: Atomic Structure, Long-Range Superstructure, Rotational Domains, and Surface Defects. Phys. Rev. B 1990, 42, 9307. 45.

Jewell, A. D.; Kyran, S. J.; Rabinovich, D.; Sykes, E. C. H. Effect of Head-Group

Chemistry on Surface-Mediated Molecular Self-Assembly. Chem. – Eur. J. 2012, 18, 7169-7178. 46.

Abel, M.; Oison, V.; Koudia, M.; Maurel, C.; Katan, C.; Porte, L. Designing a New Two‐

Dimensional Molecular Layout by Hydrogen Bonding. ChemPhysChem 2006, 7, 82-85.

ACS Paragon Plus Environment

28

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

47.

Bellisario, D. O.; Jewell, A. D.; Tierney, H. L.; Baber, A. E.; Sykes, E. C. H. Adsorption,

Assembly, and Dynamics of Dibutyl Sulfide on Au{111}. J. Phys. Chem. C 2010, 114, 1458314589. 48.

Jewell, A. D.; Tierney, H. L.; Sykes, E. C. H. Gently Lifting Gold’s Herringbone

Reconstruction: Trimethylphosphine on Au (111). Phys. Rev. B 2010, 82, 205401. 49.

Gottardi, S.; Müller, K.; Moreno-López, J. C.; Yildirim, H.; Meinhardt, U.; Kivala, M.;

Kara, A.; Stöhr, M. Cyano-Functionalized Triarylamines on Au(111): Competing Intermolecular Versus Molecule/Substrate Interactions. Adv. Mater. Inter. 2014, 1, 1300025-n/a. 50.

Järvinen, P.; Hämäläinen, S. K.; Ijäs, M.; Harju, A.; Liljeroth, P. Self-Assembly and Orbital

Imaging of Metal Phthalocyanines on a Graphene Model Surface. J. Phys. Chem. C 2014, 118, 13320-13325. 51.

Arillo-Flores, O.; Fadlallah, M.; Schuster, C.; Eckern, U.; Romero, A. Magnetic,

Electronic, and Vibrational Properties of Metal and Fluorinated Metal Phthalocyanines. Phys. Rev. B 2013, 87, 165115. 52.

Hesper, R.; Tjeng, L. H.; Sawatzky, G. A. Strongly Reduced Band Gap in a Correlated

Insulator in Close Proximity to a Metal. EPL (Europhysics Letters) 1997, 40, 177. 53.

Franke, K. J.; Schulze, G.; Henningsen, N.; Fernández-Torrente, I.; Pascual, J. I.; Zarwell,

S.; Rück-Braun, K.; Cobian, M.; Lorente, N. Reducing the Molecule-Substrate Coupling in ${\Mathrm{C}}_{60}$-Based Nanostructures by Molecular Interactions. Phys. Rev. Lett. 2008, 100, 036807. 54.

Chen, W.; Madhavan, V.; Jamneala, T.; Crommie, M. Scanning Tunneling Microscopy

Observation of an Electronic Superlattice at the Surface of Clean Gold. Phys. Rev. Lett. 1998, 80, 1469.

ACS Paragon Plus Environment

29

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 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

55.

Page 30 of 31

Ruiz-Osés, M.; de Oteyza, D. G.; Fernández-Torrente, I.; Gonzalez-Lakunza, N.; Schmidt-

Weber, P. M.; Kampen, T.; Horn, K.; Gourdon, A.; Arnau, A.; Ortega, J. E. Non-Covalent Interactions in Supramolecular Assemblies Investigated with Electron Spectroscopies. ChemPhysChem 2009, 10, 896-900. 56.

Bader, R. F. W.; Henneker, W. H.; Cade, P. E. Molecular Charge Distributions and

Chemical Binding. J. Chem. Phys. 1967, 46, 3341-3363. 57.

Kim, J.-H.; Ribierre, J.-C.; Yang, Y. S.; Adachi, C.; Kawai, M.; Jung, J.; Fukushima, T.;

Kim, Y. Seamless Growth of a Supramolecular Carpet. Nat. Commun. 2016, 7, 10653.

ACS Paragon Plus Environment

30

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

(TOC Graphic)

ACS Paragon Plus Environment

31