The growth and fluorescence of phthalocyanine monolayers, thin films and multilayers on hexagonal boron nitride

The organisation of organic materials on a surface is important in realising their potential as electronic and optoelectronic devices, and is a major focus of supramolecular chemistry and materials science. In this study we investigate the influence of surface organisation on the fluorescence of free-base phthalocyanine (H2Pc; see inset in Fig. 1(b)). Phthalocyanines, both the free-base and metallated variants, have been widely exploited in organic electronics and gas-sensing, and also provide a prototypical example of an organic semiconductor on which many fundamental studies have been undertaken. Specifically we have investigated the adsorption of H2Pc on hexagonal boron nitride (hBN) and find that well-ordered monolayers can be formed in which the molecules lie flat with lattice vectors which are aligned with a high symmetry axis of the substrate, but additional layers grow in a three-dimensional morphology which is much less ordered. While the monolayer H2Pc gives rise to a narrow intense fluorescence peak the emission from the multilayers is much broader and significantly red-shifted. H2Pc powder was purchased from Sigma-Aldrich and sublimed, under vacuum (base pressure B 10 7 mbar), onto hBN flakes which had been exfoliated, transferred to a supporting Si/SiO2 substrate and cleaned following the procedures summarised in ESI† and our recent work. The deposition rate was B0.5 monolayers per min. hBN is chosen as a substrate since it is compatible with the formation of 2D supramolecular layers, and the acquisition of images with molecular resolution using atomic force microscopy (AFM) under ambient conditions. Since hBN is an insulator the fluorescence of the adsorbed molecular layers may also be measured, facilitating investigations of the correlation between morphology and optical properties. Fig. 1(a) shows the surface after 0.7 monolayer (ML) of H2Pc was deposited on hBN while the substrate was held at room temperature. Rectangular islands are observed with typical lengths ranging from 0.3 to 2 mm, and widths from 0.2 to 0.4 mm. The profile in Fig. 1(b) shows that the height of the islands is 0.34 0.02 nm, consistent with the formation of a monolayer of H2Pc with the molecular plane lying parallel to the surface. The elongated island facets are oriented along specific directions – see histogram in Fig. 1(c) which has three peaks separated by 601. This suggests that the major axes of the islands are aligned with a high symmetry axis of hBN. If the substrate temperature is increased to 100 1C, or higher, the growth self-terminates after the formation of a monolayer. Fig. 1 (a) AC mode AFM image of 0.7 ML of H2Pc on hBN, (b) the height profile of an island along the red line marked in panel (a). (c) Histogram of island orientations.

The growth and fluorescence of phthalocyanine monolayers, thin films and multilayers on hexagonal boron nitride † Manal Alkhamisi, a Vladimir V. Korolkov Free-base phthalocyanine forms distinct interfacial phases and thin films on hexagonal boron nitride including a monolayer arrangement as determined using high resolution atomic force microscopy.
The phases reveal significant differences in photoluminescence with an intense peak for monolayer coverages of flat-lying molecules which is red-shifted in agreement with theoretical models.
The organisation of organic materials on a surface is important in realising their potential as electronic and optoelectronic devices, and is a major focus of supramolecular chemistry and materials science. 1,2 In this study we investigate the influence of surface organisation on the fluorescence of free-base phthalocyanine (H 2 Pc; see inset in Fig. 1(b)). Phthalocyanines, both the free-base and metallated variants, have been widely exploited in organic electronics and gas-sensing, and also provide a prototypical example of an organic semiconductor on which many fundamental studies have been undertaken. [3][4][5][6][7][8][9][10][11][12] Specifically we have investigated the adsorption of H 2 Pc on hexagonal boron nitride (hBN) and find that well-ordered monolayers can be formed in which the molecules lie flat with lattice vectors which are aligned with a high symmetry axis of the substrate, but additional layers grow in a three-dimensional morphology which is much less ordered. While the monolayer H 2 Pc gives rise to a narrow intense fluorescence peak the emission from the multilayers is much broader and significantly red-shifted. H 2 Pc powder was purchased from Sigma-Aldrich and sublimed, under vacuum (base pressure B 10 À7 mbar), onto hBN flakes which had been exfoliated, transferred to a supporting Si/SiO 2 substrate and cleaned following the procedures summarised in ESI † and our recent work. 13 The deposition rate was B0.5 monolayers per min. hBN is chosen as a substrate since it is compatible with the formation of 2D supramolecular layers, and the acquisition of images with molecular resolution using atomic force microscopy (AFM) under ambient conditions. Since hBN is an insulator the fluorescence of the adsorbed molecular layers may also be measured, facilitating investigations of the correlation between morphology and optical properties. 14,15 Fig. 1(a) shows the surface after 0.7 monolayer (ML) of H 2 Pc was deposited on hBN while the substrate was held at room temperature. Rectangular islands are observed with typical lengths ranging from 0.3 to 2 mm, and widths from 0.2 to 0.4 mm. The profile in Fig. 1(b) shows that the height of the islands is 0.34 AE 0.02 nm, consistent with the formation of a monolayer of H 2 Pc with the molecular plane lying parallel to the surface. The elongated island facets are oriented along specific directions -see histogram in Fig. 1(c) which has three peaks separated by 601. This suggests that the major axes of the islands are aligned with a high symmetry axis of hBN.
If the substrate temperature is increased to 100 1C, or higher, the growth self-terminates after the formation of a monolayer. The resulting islands are large (micron scale) and it is possible to acquire images with molecular scale resolution -see Fig. 2(a, b and d). The lattice vectors are marked in Fig. 2(b); we find that a Pc = b Pc = 2.18 AE 0.05 nm and the angle y = 861. The H 2 Pc molecules appear as topographically bright features with a 'crosslike' shape similar to the molecular structure ( Fig. 1(b) inset). The orientation of the lattice vectors (red arrows) of the H 2 Pc monolayer and the hBN (blue arrows) can be determined by comparing AFM images of the H 2 Pc and the hBN lattice (see Fig. 2(c and d)). The lattice vectors of hBN may be resolved in Fig. 2(c) (parallel to the blue arrows) and are overlaid on the H 2 Pc molecules in Fig. 2(d). From these images we determine that the lattice vector b Pc is aligned at angle of 301 to both hBN lattice vectors, one of the high symmetry directions of the hBN substrate, consistent with the histogram in Fig. 1(c). The axis of the H 2 Pc molecules (parallel to the C-H bonds within the central ring of the macrocycle) is misaligned from the hBN lattice vector by 4 AE 21.
A schematic of the H 2 Pc layer overlaid on hBN is shown Fig. 2(e). The molecules appear alternately bright (at the corners of the unit cell in Fig. 2(b and e)) and dark (at the centre of the unit cell). We suggest the reason for the differing AFM contrast of molecules is due to a variation in adsorption site; in our previous studies of molecular overlayers on hBN we have reported the appearance of moiré patterns 13 in AFM images arising from a mismatch in lattice constants, and thus variations of local adsorption site. The H 2 Pc unit cell is similar to that previously reported on other substrates. [16][17][18][19] We have calculated the adsorption energy and geometry of H 2 Pc using density functional theory (DFT). The structure of H 2 Pc on hBN is optimised using the range-separated hybrid oB97X-D functional, 20 which includes an empirical dispersion correction, in combination with the correlation-consistent cc-pVDZ basis set. 21 The atomic positions of a model hBN monolayer flake consisting of 65 boron atoms and 65 nitrogen atoms were initially optimized and frozen in subsequent calculations. The optimized structure for a single H 2 Pc molecule on hBN is shown in Fig. 3 and confirms a planar adsorption geometry with a molecule-surface separation, d = 0.33 nm and adsorption energy, E A = 3.35 eV. In common with our images the calculations show the molecular axis misaligned with the hBN lattice vector, but by a larger angle, 131. We attribute this difference to a combination of a low barrier for rotation and the role of intermolecular interactions. Full details of the calculations, which are performed within the QCHEM package, 22 are included in ESI. † As the H 2 Pc coverage is increased we observe the formation of clusters of large needle-like crystallites with a rather disordered morphology. Fig. 4(a) shows a large area AFM image of an hBN flake in which several large islands (identified by arrows) are formed following the deposition of H 2 Pc (thickness 1 nm; 1 ML is approximately equivalent to 0.33 nm). These islands have a height of 20-30 nm (see Fig. 4(c)) and typical length of 10-20 mm, and can be resolved using optical microscopy; Fig. 4(b) shows the same hBN flake. Interestingly the islands have a direction which is locally correlated and oriented in one of three directions subtended by 601. Following the deposition of thicker layers (Fig. 4(d)) the needle-like islands merge to cover a greater proportion of the surface, but are highly disordered. We are not able to obtain images with molecular-scale resolution in these regions. Overall our results show an abrupt transition to the growth of three-dimensional crystallites following the completion of a highly ordered monolayer. A likely explanation for this is that molecules switch from the flat-lying adsorbed state in the monolayer to a geometry in which the molecular plane is at an angle to the hBN substrate in higher layers, thus forming one-dimensional crystallites in which molecules are co-planar, maximising p-p interactions. Such a switch from a flat-lying to upright morphology has been observed previously in studies of adsorbed H 2 Pc. 4 To measure the optical properties of a H 2 Pc thin film, a Horiba LabRAM system was used in conjunction with a  660 nm laser (photon energy 1.88 eV) with a spot size B 1 mm 2 as an excitation source (see ESI † for further details). In Fig. 5 the photoluminescence spectra of H 2 Pc films with varying thickness are shown. For a thickness of 0.7 ML, in which all molecules adsorb in a flat-lying configuration, we observe a narrow intense peak centred at 1.75 eV, which we attribute to the zero phonon (0-0) peak, together with broader, less intense vibronic peaks at 1.67 and 1.58 eV. These peaks are also present for a sample with H 2 Pc thickness of 1 nm (as in Fig. 2), but undergo a small shift; the 0-0 peak is centred on 1.76 eV. There are small differences between the spectra acquired when the laser spot is incident either on or off a multilayer H 2 Pc island. In particular, an additional peak centred at 1.38 eV is slightly more intense on these islands. This peak is much more intense for a thicker sample (32 nm), for which the peaks associated with the flat-lying monolayer may no longer be clearly resolved.
We have also measured the spectrum of H 2 Pc powder on a glass slide which showed a broad peak at 1.42 eV (light blue line in Fig. 5) and we attribute the peak observed at 1.38 eV for the 32 nm thick H 2 Pc layer to the existence of a bulk-like H 2 Pc phase on the surface. The 0-0 peak for H 2 Pc is red-shifted by DE = 0.13 eV as compared with a molecule in the gas phase 23 (0-0 peak at 1.876 eV). As we, 14,15 and others, [24][25][26] have recently discussed this red-shift can be understood from a combination of effects arising from interactions with the substrate and with nearest neighbours. A non-resonant contribution to the shift due to the substrate can be calculated using time-dependent density functional theory (TD-DFT). The molecule is optimized in the excited (S 1 ) state both on (adsorbed) and off (gas phase) the hBN surface and the transition to the ground state (S 0 ) is calculated in each case (see ESI † for more details). The calculated transition energies for a molecule on the hBN surface, and in the gas phase are, respectively, 1.76 eV and 1.87 eV. These values are close to the measured peaks in fluorescence and their difference implies a non-resonant red-shift DE nr = 0.11 eV due to adsorption on hBN.
There are also contributions to the red-shift from resonant interactions with the substrate and neighbouring molecules. We have recently argued 15 that these effects depend on the refractive index, n, of the substrate and include a 'gas-surface' red-shift given by DE subs = f (d)(n 2 À 1)m 2 /4pe 0 (n 2 + 1)d 3 , where m is the transition dipole moment and d is the molecule-substrate separation, which from our DFT calculations (see ESI †) are, respectively, 8.0 Debye and 0.33 nm (for the anisotropic hBN substrate we use a geometric average 15 n = 1.87; the reduction factor f (d) is related to the molecular geometry as discussed below). This shift arises from a coupling of the transition dipole moment of the adsorbed molecule with the dielectric hBN substrate and may be described semi-classically as an interaction with an image dipole (see Kerfoot et al. 15 ).
In the point dipole approximation f (d) = 1; this is satisfied if d is much greater than the molecular dimensions. However, this condition is violated for H 2 Pc (see dimensions in Fig. 3); instead we treat the dipole moment as arising from an oscillating charge density, 27 the transition density, T(r), equal to the product of the wave functions of the lowest unoccupied molecular orbital (LUMO) and the highest occupied molecular orbital (HOMO). Previously 15 we have used the extended dipole model to determine an approximate expression for f (d), however this is expected to be valid only for linear molecules. For the cross-shaped H 2 Pc molecule we instead calculate the reduction factor due to the interaction of the charge density described by T(r) with its image, T 0 (r 0 ), given by  f ðdÞ ¼ A Ð Ð drdr 0 TðrÞTðr 0 Þ=4pe 0 jr À r 0 þ dẑj, where the constant A = m 2 /4pe 0 d 3 and T 0 (x 0 , y 0 , z 0 ) = T(x 0 , y 0 , Àz 0 ). The isosurfaces of the HOMO, LUMO and the transition density are shown in Fig. 3. Evaluating numerically this double integral using the transition density calculated using DFT gives, for d = 0.33 nm, f (d) = 0.106 and DE subs = 0.03 eV.
We might also expect a red-shift due to a screened resonant interaction with neigbouring molecules; 24 the characteristic energy for a pairwise interaction is DE nn = m 2 /2(n 2 + 1)pe 0 r nn 3 , where r nn = 1.5 nm, the nearest neighbour separation. We find DE nn = 5 meV; this red-shift is expected for molecules which have aligned transition dipole moments. From our DFT calculations we find that the transition dipole moment for H 2 Pc is oriented close to perpendicular to the C-H bonds within the central ring of the macrocycle. From our images we are unable to determine whether the molecules are aligned (i.e. whether the C-H bonds within all molecules have a common alignment). However, the contrast of the islands, which may be resolved in optical microscopy, shows no dependence on polarisation indicating that the transition dipole moments are not ordered. Assuming a random alignment, the total red-shift due to N nearest neighbours is Bw 2 NDE nn /2 where w 2 is the Franck Condon factor. 24 Since w 2 r 1 we estimate that this contribution to the gas-surface red-shift is B10 meV. Combining the various calculated contributions to the red-shift gives DE tot = DE nr + DE subs + DE nn in the range 0.14-0.15 eV in excellent agreement with the experimental value, 0.13 eV.
Overall our results demonstrate that the optoelectronic properties of molecules on surfaces depends strongly on their adsorption geometry and environment as determined, with molecular resolution, using atomic force microscopy. The absorption of H 2 Pc in a flat lying geometry results in a narrow intense emission peak, but evolves to a broad emission feature for multilayer samples.

Conflicts of interest
There are no conflicts of interest to declare.