# Optical properties and electrical transport of thin films of terbium(III) bis(phthalocyanine) on cobalt

Peter Robaschik, Pablo F Siles, Daniel Bülz, Peter Richter, Manuel Monecke, Michael Fronk, Svetlana Klyatskaya, Daniel Grimm, Oliver G Schmidt, Mario Ruben, et al.

► To cite this version:

Peter Robaschik, Pablo F Siles, Daniel Bülz, Peter Richter, Manuel Monecke, et al. Optical properties and electrical transport of thin films of terbium(III) bis(phthalocyanine) on cobalt. Beilstein Journal of Nanotechnology, 2014, 5, pp.2070 - 2078. 10.3762/bjnano.5.215. hal-05032642

HAL Id: hal-05032642

https://hal.science/hal-05032642v1

Submitted on 14 Apr 2025

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Optical properties and electrical transport of thin films of terbium(III) bis(phthalocyanine) on cobalt

Peter Robaschik¹, Pablo F. Siles²,³, Daniel Bülz¹, Peter Richter¹, Manuel Monecke¹, Michael Fronk¹, Svetlana Klyatskaya⁴, Daniel Grimm²,³, Oliver G. Schmidt²,³, Mario Ruben⁴,⁵, Dietrich R. T. Zahn¹ and Georgeta Salvan*¹

Full Research Paper

Open Access

Address:

¹Semiconductor Physics, Technische Universität Chemnitz, Reichenhainer Straße 70, 09107 Chemnitz, Germany, ²Material Systems for Nanoelectronics, Technische Universität Chemnitz, Reichenhainer Straße 70, 09107 Chemnitz, Germany, ³Institute for Integrative Nanosciences, IFW Dresden, Helmholtzstraße 20, 01069 Dresden, Germany, ⁴Institute of Nanotechnology, Karlsruhe Institute of Technology (KIT), 76344 Eggenstein-Leopoldshafen, Germany, and ⁵Université de Strasbourg, Institut de Physique et de Chimie des Materiaux de Strasbourg, CNRS UMP 7504, 23 Rue du Loess, 67034 Strasbourg Cedex 2, France

Email:

Georgeta Salvan* - salvan@physik.tu-chemnitz.de

Keywords:

current sensing AFM; ellipsometry; spintronics; TbPc₂; transport properties

Beilstein J. Nanotechnol. **2014**, *5*, 2070–2078.
doi:10.3762/bjnano.5.215

Received: 14 July 2014

Accepted: 15 October 2014

Published: 11 November 2014

This article is part of the Thematic Series "Molecular materials – towards quantum properties".

Associate Editor: P. Ziemann

© 2014 Robaschik et al; licensee Beilstein-Institut.

License and terms: see end of document.

## Abstract

The optical and electrical properties of terbium(III) bis(phthalocyanine) (TbPc₂) films on cobalt substrates were studied using variable angle spectroscopic ellipsometry (VASE) and current sensing atomic force microscopy (cs-AFM). Thin films of TbPc₂ with a thickness between 18 nm and 87 nm were prepared by organic molecular beam deposition onto a cobalt layer grown by electron beam evaporation. The molecular orientation of the molecules on the metallic film was estimated from the analysis of the spectroscopic ellipsometry data. A detailed analysis of the AFM topography shows that the TbPc₂ films consist of islands which increase in size with the thickness of the organic film. Furthermore, the cs-AFM technique allows local variations of the organic film topography to be correlated with electrical transport properties. Local current mapping as well as local I–V spectroscopy shows that despite the granular structure of the films, the electrical transport is uniform through the organic films on the microscale. The AFM-based electrical measurements allow the local charge carrier mobility of the TbPc₂ thin films to be quantified with nanoscale resolution.

<footnote>* Corresponding author</footnote>
Introduction

Molecular spintronic devices could bring a new era of informa-
tion technology, as the materials are inexpensive and have a
potentially higher efficiency than conventional electronic
devices [1-9]. Therefore, many studies were carried out to iden-
tify organic molecules with suitable properties for spintronics
over the past few years. Terbium(III) bis(phthalocyanine)
(TbPc2) is an excellent candidate to provide all the necessary
features for molecular spintronics, as it is both an organic semi-
conductor and a single molecule magnet (SMM). TbPc2 was
previously implemented in an organic field effect transistor
(OFET) as a hole transporting layer [10] and recently
Urdampilleta et al. reported a supramolecular spin valve made
of a carbon nanotube (CNT) covered by only a few TbPc2
molecules [5]. For electronic and spintronic devices it is crucial
to know and to control the molecular orientation on the device-
related substrates. The TbPc2/Co heterojunction was already
proposed to serve as a model system for a SMM semicon-
ducting layer on top of a ferromagnetic electrode for a future
spintronic device. The chemical and magnetic properties of this
interface were investigated by Klar et al. and it was found that
the magnetic moment of the Tb couples antiferromagnetically to
the Co substrate [11]. In this work we focus on the study of
other device-relevant aspects: the influence of the film thick-
ness, morphology, and molecular orientation on the electrical
transport in TbPc2 layers on polycrystalline cobalt films. The
TbPc2 molecule and the investigated layer stack are schemati-
cally shown in Figure 1, which includes a sketch of the molec-
ular orientation which will be discussed later. The present
manuscript reveals the optical and electrical properties of TbPc2
films with different thicknesses on Co substrates as well as the
molecular tilt angle and grain size distribution of the samples

using spectroscopic ellipsometry, AC atomic force microscopy,
and current sensing atomic force microscopy. Topographic and
electrical AFM techniques provide a reliable method to investi-
gate and correlate the structural and local electrical properties of
TbPc₂ thin films. This knowledge is crucial for the implementa-
tion and fabrication of TbPc₂-based devices.

<figure><img src="image_2.png" /><figcaption>**Figure 1:** TbPc₂ molecule (left). Investigated layer stack: TbPc₂ thin films on cobalt grown on SiO₂/Si(111).</figcaption></figure>

Results and Discussion
Spectroscopic ellipsometry

For the ellipsometric analysis, four TbPc₂ films with different thicknesses (18 ± 1 nm, 41 ± 1 nm, 58 ± 2 nm, and 87 ± 3 nm) were prepared by organic molecular beam deposition. The underlying 25 nm thick cobalt film was prepared by electron beam evaporation. Both depositions were performed without breaking the vacuum in between to avoid any oxidation of the Co films. Figure 2 shows the dielectric function, $\tilde{\epsilon} = \epsilon_1 + i\epsilon_2$, of Co taken from [12] and a 58 nm TbPc₂ film on Co. TbPc₂ films exhibit a uniaxial anisotropy of the dielectric function similar to

<figure><img src="image_3.png" /><figcaption>**Figure 2:** Dielectric function of a TbPc₂ film on cobalt. The blue lines and the red lines represent the real part (left) and the imaginary part (right), respectively, of the dielectric function. The top graphs show the dielectric function of Co and bottom graphs are obtained from 58 nm TbPc₂ on top of the Co layer.</figcaption></figure>
many planar phthalocyanines, for example, CuPc [13] or H₂Pc [14]. Consequently, the dielectric function parallel (in-plane) and perpendicular (out-of-plane) to the sample surface differs ($ε_{xx} = ε_{yy} ≠ ε_{zz}$). For a numerical analysis of the measured spectroscopic ellipsometry data a model employing eleven Gaussian oscillators was used. The energy position and the full width half maximum (FWHM) of the oscillators were fixed for all samples, while their amplitudes were allowed to vary during the Kramers–Kronig consistent fitting procedure (experimental data and model fit for one sample: Supporting Information File 1, Figure S1). From the numerical analysis we can extract the anisotropic dielectric function as well as the thickness of the different layers in the stack.

The most prominent absorption bands, namely the Q and B band, of phthalocyanines are highlighted in Figure 2. They correspond to ligand-related π–π* transitions [15].

For organic semiconductors we can consider the relative magnetic permeability to be $μ_r ≈ 1$. Thus, we can easily determine the optical constants ($\tilde{n} = n + ik$) from the diagonal elements $ε_{ii}$ of the dielectric tensor using Equation 1:

$$n_{ii}^2(\omega) = \mu_r \cdot \varepsilon_{ii}(\omega) \approx \varepsilon_{ii}(\omega), \quad (1)$$

where $i = x, y, z$.

From the anisotropy of the extinction coefficient $k$ in the Q band region (1.6–2.0 eV), it is possible to estimate the molecular orientation [13,14]. According to the model described in [16] we assume that two electronic transition dipole moments in the Q band are parallel to the ligand plane and that all molecules have the same tilt angle $α$ with respect to the substrate, but with a random azimuthal orientation in the substrate plane. It was shown that under these assumptions the average molecular tilt angle $α$ can be written as [16]:

$$\alpha = \arccos \sqrt{\frac{2 - A_{\text{out}} / A_{\text{in}}}{2 + A_{\text{out}} / A_{\text{in}}}}, \quad (2)$$

where $A_{in}$ and $A_{out}$ are the areas under the in-plane and out-of-plane component of the extinction coefficient, respectively (see Figure 2). The resulting average angles between the molecules and the substrate plane can be found in Figure 3. The Sessoli group reported TbPc₂ molecules which lay in a monolayer evaporated on polycrystalline gold and cobalt as well as standing molecules in a 200 nm thick film by exploiting synchrotron methods [17,18]. Thus, we expect that the first TbPc₂ layers will similarly lie on the polycrystalline Co films

due to their low roughness (rms: 0.4 nm, obtained by AFM measurements). At this point it should be mentioned that the ellipsometry investigations of the molecular orientation in films with smaller thickness is difficult due to the low sensitivity to the out-of-plane component of the optical constants. In particular, this is due to the presence of a metallic layer beneath the organic layer. The molecular tilt angle increases with increasing films thickness, which was also shown for H₂Pc on PTCDA by utilizing spectroscopic ellipsometry and magneto-optical Kerr effect spectroscopy (MOKE) [16]. According to Equation 2, only values between 0 and 2 are allowed for $A_{out}/A_{in}$. Nevertheless, for the thickest TbPc₂ layer, the ratio $A_{out}/A_{in}$ slightly exceeds the limit of 2. This could be related to errors in the ellipsometry fit. Therefore, we assume standing molecules for films thicker than 80 nm.

<figure><img src="image_4.png" /></figure>

<table><caption>**Figure 3:** Definition of the molecular tilt angle (top). Average tilt angle of the TbPc₂ molecules on cobalt (bottom). The thickness of the films was estimated from the ellipsometry data.</caption><thead><tr><th colspan="3">Co</th></tr><tr><th>Thickness d / nm</th><th>Ratio A<sub>out</sub> / A<sub>in</sub></th><th>Angle α / deg</th></tr></thead><tbody><tr><td>18 ± 1</td><td>1.48 ± 0.05</td><td>67 ± 2</td></tr><tr><td>41 ± 1</td><td>1.87 ± 0.05</td><td>79 ± 2</td></tr><tr><td>58 ± 2</td><td>1.97 ± 0.05</td><td>85 ± 4</td></tr><tr><td>87 ± 3</td><td>2.02 ± 0.05</td><td>90 ± 5</td></tr></tbody></table>

## Topography characteristics

To investigate the topographic characteristics of the TbPc₂ thin films on Co substrates, a detailed analysis of the grain size evolution is performed as a function of the organic film thickness. Figure 4 shows AFM images with areas of 2 × 2 µm² for TbPc₂ films as well as the Co substrates. An increase of the roughness is observed as a function of the thickness of the organic film (see inset in Figure 5a). Scan profiles in Figure 4 show the average grain heights over the sample surface. These profiles also support a clear variation of the topographic characteristics of the TbPc₂ films as the thickness is increased.

In order to perform a detailed quantification of the TbPc₂ topographic grain characteristics, a statistical analysis via
<figure><img src="image_5.png" /><figcaption>**Figure 4:** AFM topography characteristics of TbPc₂ thin films. Line scan profiles and AFM surface images for TbPc₂ films of 18, 41, 58 and 87 nm deposited on a 25 nm thick Co film.</figcaption></figure>

histograms is used to calculate the average grain diameter and height from the topography images shown in Figure 4. Figure 5a shows that the height of the grains follows a linear increase while the average grain diameter tends to reach a maximum size of about 38–40 nm with increasing organic film thickness, as expected for an unheated substrate during the deposition [19]. As an example of the statistical analysis performed, Figure 5b shows a histogram of the grain height and diameter for the case of the 87 nm TbPc₂ film. This analysis considers areas of 2 × 2 µm² shown in Figure 4, which contain approximately 2 × 10³ grains on the surface. The statistical analysis was performed on different locations of the organic films, revealing similar results. This reinforces the reliability of the data and dismisses the possibility of error during the AFM measurements on a particular location of the sample surface.

Further information on the statistical analysis can be found in Supporting Information File 1, Figure S2 and Figure S3.

<figure><img src="image_6.png" /><figcaption>**Figure 5:** AFM statistical analysis of TbPc₂ thin films. (a) Average grain diameter and height as a function of the organic film thickness. Dotted lines are guide to the eye to show the tendency of the grain height and grain diameter with the thickness. The inset shows the dependence of the roughness of the organic films thickness. (b) Diameter and height histograms for an 87 nm TbPc₂ film. Error bars (with sizes comparable to plotted dot symbols) in (a) are obtained from Gaussian fits in histograms as shown in (b).</figcaption></figure>

## Current sensing atomic force microscopy

Conductive atomic force microscopy techniques are well-established methods for local electrical characterization in organic materials [20-24]. In this work, we employ the cs-AFM technique in order to investigate the local transport properties of TbPc₂ thin films on Co substrates. Due to its high reproducibility and versatility, cs-AFM allows for local current-voltage (I-V) spectroscopy as well as current mapping of particular areas of interest on the organic film. I-V spectroscopy is realized by subsequent sweeping of the applied bias while the AFM conductive probe is located at a fixed location on the sample surface. This procedure is repeated several times over different locations to ensure the reproducibility of the electrical response. On the other hand, in current mapping experiments, a fixed voltage is applied at the probe-sample interface while the probe is scanned over a specific area of the sample surface. These high resolution current maps (512 × 512 points) allow simultaneous information to be obtained regarding
topographic and electrical current through the organic material, and establish a direct correlation of the organic topography with electrical characteristics. Also, by repeating current mapping at different locations and with different voltages, we are able to reconstruct the I–V characteristics.

Figure 6a shows a schematic diagram of the set up for local electrical measurements. A conductive AFM probe placed directly in contact with the TbPc₂ surface plays the role of the top electrode, while the Co bottom film acts as a back electrode. Samples for cs-AFM measurements were deposited on a Si substrate with a top SiO₂ layer of 1 µm in order to eliminate possible leakage current.

Figure 6b shows a 5 × 5 µm² topography image of an 80 nm thick TbPc₂ film which has the highest roughness, with respect to Figure 4 and Figure 5. It is worth pointing out that for the case of electrical AFM measurements, special conductive probes with a larger radius compared to the high resolution probes used in the topographic studies (see Experimental section for more details) are utilized. Therefore, a lower topography resolution is expected in topography images acquired during electrical measurements (e.g., Figure 6b). The electrical response of the film, for the same location indicated in Figure 6b, is explored under different applied voltages, as indicated in Figure 6c–e. We observe highly stable and uniformly transport characteristics for all TbPc₂ film thicknesses investigated. This suggests a uniform distribution of the electrical charge flow through the organic film.

To further investigate the transport mechanism in TbPc₂ thin films, a series of local I–V spectroscopy measurements on different locations along the organic material were performed. Figure 7a shows the transport response for the case of 20 and

80 nm thick TbPc₂ films. Here, the AFM probe is fixed at one single point on the surface while the voltage is swept for around 20 consecutive cycles (grey and black areas). Solid lines show the average electrical current. In order to verify the reproducibility of the I–V spectroscopy results, a series of current maps were also obtained at different applied voltages for both organic films. I–V characteristics were then reconstructed by obtaining the average current corresponding to 512 × 512 data

<figure><img src="image_7.png" /><figcaption>**Figure 6:** cs-AFM electrical measurements. (a) Electrical setup employed for local electrical measurements via cs-AFM. (b) AFM topography image of an 87 nm thick TbPc₂ film (5 × 5 µm²). Current maps for the same location indicated in (b) for the case of an applied voltage of 0.6 V (c), 1.0 V (d) and 1.5 V (e).</figcaption></figure>

<figure><img src="image_8.png" /><figcaption>**Figure 7:** Transport mechanism for TbPc₂ thin films. Red and blue solid lines indicate the average of 20 local I–V spectroscopy cycles. (a) Current–voltage characteristics for TbPc₂ thin films. Grey and black data in (a) represent the local I–V spectroscopy cycles. Dotted data represent the current value obtained via current maps. Each dot represents the average of 512 × 512 data points acquired during the AFM scanning. (b) Current density response. Solid green lines and dashed blue lines correspond to the linear and space-charge-limited transport regimes, respectively.</figcaption></figure>
points from current maps as the ones shown in Figure 6c–e. Solid dots in Figure 7a correspond to the I–V characteristics reconstructed with current maps and indicate the high reproducibility of the transport measurements for TbPc₂ organic films performed via cs-AFM.

The current values measured for the TbPc₂ thin films are close to those reported for the case of CuPc [23] and up to three orders of magnitude higher than those reported for similar organic films investigated via cs-AFM techniques such as metalloporphyrin thin films on Ni substrates [20]. In general, due to their planar structure, phthalocyanines are known to exhibit a higher hole mobility as compared to the porphyrins [25], which could lead to a higher current. When comparing the topographic characteristics of TbPc₂ samples with 20 nm and 80 nm (see Figures 4 and 5), the variation with grain size would induce more grain boundaries for electron scattering in the thinner samples when compared with the ticker samples. Thus, a different amount of scattering centers might be responsible for the fact that the electric current scales gentler than the expected laws with the film thickness (1/L).

Figure 7b shows the current density–voltage characteristics for the TbPc₂ films. In order to obtain the current density from the cs-AFM measurements, we follow a similar approach as presented by Reid et al. [26], where the contact area between probe and sample is determined assuming a tip indentation of 1 nm for the same kind of Pt-coated hemispherical probes used in this work. This results in a circular contact area with a diameter ($P_d$) of 14 nm. Here, we ensure that the probe–sample force is kept constant during the electrical measurements and no extra force, which could eventually modify the contact area, is applied.

The electrical response of the TbPc₂ thin films presents a transition from a linear ohmic-like transport regime for low voltages to a square law dependence for high voltages. These results appear to be in agreement with a space-charge-limited current process (SCLC). According to G. Horowitz et al., the linear current–voltage characteristics can exist in the SCLC model and come from electrons hopping from one insulating state to the next [27]. Hence, the transport regime of TbPc₂ thin films would follow such transition behavior:

$$J_{\text{ohmic}} = \frac{\sigma V}{L}, \quad J_{\text{sc}} = \frac{9}{8} \varepsilon \varepsilon_0 \mu \frac{V^2}{L^3}, \quad (3)$$

where $J_{\text{ohmic}}$ and $J_{\text{sc}}$ are the current densities for each regime, $\sigma$ is the low voltage conductivity, $L$ is the thickness of the organic film, $\varepsilon$ is the relative dielectric constant, $\varepsilon_0$ is the

permittivity of free space, $\mu$ is the charge carrier mobility and $V$ is the applied voltage. From the ohmic-like regime (see Figure 7b) we obtained parameters of $\sigma/L$ equal to 1.38 and 0.97 kA/cm² V for the samples of 20 and 80 nm, respectively. If the mobility is higher along the crystalline direction, where an overlap of the $\pi$ orbitals of the neighboring molecules exists, then the mobility might be higher perpendicular to the molecular plane as opposed to parallel to it. Assuming the same kind of charge for the carriers and the same charge density per unit volume in all films, a reduction in the ratio $\sigma/L$ corresponds to a reduction in the carrier mobility in the direction perpendicular to the substrate. This is in agreement with the increasingly standing molecules. The validity of the model applied for the ohmic region should be granted in both cases, since the condition is that the amount of injected carriers is low compared to that of thermally generated carriers.

Once the modelling of the current density characteristics for TbPc₂ thin films is derived from cs-AFM measurements, it is necessary to account for the electrode geometry in the AFM electrical experiments (AFM probe geometry). We apply the modelling approach described in [26], as we use the same AFM probe geometry. Here, the system is treated as an intermediate case between a situation with infinite plane–plane electrodes and a situation with point–plane electrodes. For the electrode-modified electrical characteristics, a semi-empirical expression for the current density determined via cs-AFM is given by [23,26]:

$$J_{\text{sc}} = 8.2 \varepsilon \varepsilon_0 \mu V^2 (7.8 \pm 1) \left(\frac{1}{P_d}\right)^{1.6 \pm 0.1} L^{-1.4}. \quad (4)$$

By comparing Equations 3 and 4, one can see that current density measurements via cs-AFM consider a thickness dependence from $L^{-3}$ to $L^{-1.4}$ for the case of the Pt-coated AFM probes used in this work [26]. Applying this geometry-corrected model to our current density measurements, it is possible to quantify the charge carrier mobility for the TbPc₂ thin films on a nanometer scale. The values of $L$ and $P_d$ are fixed by the film thickness and tip geometry, respectively. For the dielectric constant we consider a minimum value of $\varepsilon = 4.5$, which is the real part of the dielectric function at the lowest photon energy (1.3 eV) used in our ellipsometry experiment (see Figure 2). As a maximum value for the dielectric constant we used $\varepsilon = 13$, which was previously proposed to describe the hole transport along the phthalocyanine ligand plane in CuPc films [28]. This would be the case if all phthalocyanine molecules were perfectly aligned with their ligand plane perpendicular to the electrodes. From the SCLC fitting in Figure 7b we obtain charge carrier mobility values from $0.80 \times 10^{-4} \text{ cm}^2 \text{ V}^{-1} \text{ s}^{-1}$
(for ε = 13) to 2.31 × 10⁻⁴ cm² V⁻¹ s⁻¹ (for ε = 4.5) for the 20 nm TbPc₂ film. For the 80 nm thick film, the mobility values range from 0.35 × 10⁻⁴ cm² V⁻¹ s⁻¹ (for ε = 13) to 1.01 × 10⁻⁴ cm² V⁻¹ s⁻¹ (for ε = 4.5). The variation of the mobility values across the film surface in a scan window of 5 × 5 µm² is below 15% (see Supporting Information File 1, Figure S4). These values are comparable to the hole field effect mobility values determined from OFET measurements by Katoh et al. [10] Noteworthy, the mobility is expected to increase with increasing size of the crystalline grains in the film [29]. A decrease in the hole mobility values with increasing film thickness is therefore at a first glance surprising. However, considering the anisotropic hole transport in the phthalocyanines (with higher mobility in the direction of the π–π stacking in a molecular crystal, i.e., perpendicular to the Pc ligand), a change in the molecular orientation from more lying to more standing molecules with increasing film thickness can be responsible for the decrease in the hole mobility values with increasing film thickness. By taking into account the homogeneous current maps recorded at different voltages (see Figure 6), the excellent agreement between the local I–V characteristics and the average current values obtained via current maps (Figure 7), we conclude that the mobility values estimated above are representative for the area of the respective TbPc₂ film.

## Conclusion

In this work we present the optical, topographic and electrical properties of TbPc₂ thin films on cobalt by utilizing ellipsometry and AFM techniques. The ellipsometric studies allowed us to determine the average molecular tilt angle in the TbPc₂ films and this evaluation revealed an evolution from nearly lying molecules in the first layers to standing molecules in a thick film. The current flow through our organic layers is homogeneous within a standard deviation of about 10%, with lower values at the grain boundaries as compared to the top of the grains (see Figure S5 of Supporting Information File 1). A statistical analysis was conducted to determine the size of the grains and it was shown that the lateral expansion of the grains appears to saturate in films at a thickness higher than 58 nm. The I–V characteristics indicate that the transport through the films with thicknesses of 20 nm and 80 nm is governed by the SCLC regime. By applying a SCLC model adapted for the I–V characteristics obtained from cs-AFM measurements, we estimate the hole mobility in TbPc₂ films on cobalt substrates to be in the range from 0.35 × 10⁻⁴ cm² V⁻¹ s⁻¹ to 2.31 × 10⁻⁴ cm² V⁻¹ s⁻¹, depending on the film thickness and the dielectric constant considered. The AFM-based approach implemented here allows important transport properties such as current density homogeneity and the local charge carrier mobility to be quantified. The nanoscale resolution achieved

here for the characterization of organic systems such as TbPc₂ thin films is crucial for future molecular spintronics applications.

## Experimental

### Sample preparation

Cobalt films were grown by electron beam evaporation on a Si(111) substrate covered by an oxide layer with different thickness (VASE, AFM: 1.5 nm and cs-AFM: 1 µm) at a rate of 1.25 nm/min under UHV conditions (10⁻⁸ mbar). The substrates were cleaned in acetone and ethanol for 5 minutes each in an ultrasonic bath. On top of the cobalt, the TbPc₂ films were prepared by organic molecular beam deposition at a rate of 0.5 nm/min at a pressure below 10⁻⁷ mbar. The evaporation took place at a temperature of about 400 °C in the Knudsen cell. The samples were kept at room temperature during all depositions. The preservation of the molecules in a film was checked by UV–vis and Raman spectroscopy.

### Ellipsometry measurements

VASE measurements were performed ex situ with a Woollam T-Solar Ellipsometer in the spectral range of 0.7–5.0 eV with an energy step width of 0.02 eV. Three different angles of incidence (50°, 60° and 70°) were exploited for increased sensitivity of the optical anisotropy of the films. The initially linear polarized light becomes elliptically polarized light after reflection on the sample. The elliptical polarization state is described by the experimentally measured quantities Ψ and Δ, according to

$$ \rho = \frac{r_p}{r_s} = \tan(\Psi)e^{i\Delta}, \quad (5) $$

where ρ is the ratio between the Fresnel reflection coefficients for p- and s-polarized light. The ellipsometry data were analyzed using the WVASE 32 software from J. A. Woollam Co. [30].

### AFM measurements

AFM measurements for topography analysis and electrical characterization were performed on an Agilent 5500 AFM system. Measurements were performed under a controlled N₂ environment to preserve the integrity and avoid exposure of the organic films to ambient conditions. Topography measurements were performed in AC tapping mode, which guarantees minimal contact between the AFM probe and the organic film. Ultra sharp (4–10 nm radius) Olympus cantilevers allowed high sensitivity measurements. cs-AFM measurements were performed in contact mode using special Pt-coated Si cantilevers with a spring constant of 0.2 N/m and typical radii of about 20–25 nm. The voltage is applied directly to the
bottom Co electrode. The grounded conductive cantilever is therefore used as a top electrode for local I–V spectroscopy as well as current mapping experiments. AFM topography analysis and current maps images were analyzed using WSxM and Gwyddion software packages [31,32].

## Supporting Information

The Supporting Information shows the raw data obtained from ellipsometry and the corresponding model fit for one sample. Furthermore, a statistical analysis of the AFM data is included.

### Supporting Information File 1

Ellipsometry and AFM analysis.
[http://www.beilstein-journals.org/bjnano/content/supplementary/2190-4286-5-215-S1.pdf]

## Acknowledgements

We gratefully acknowledge the financial support from the Deutsche Forschungsgemeinschat (DFG) in the frames of the Research Unit FOR 1154 “Towards Molecular Spintronics”.

## References

1. Rocha, A. R.; Garcia-Suárez, V. M.; Bailey, S. W.; Lambert, C. J.; Ferrer, J.; Sanvito, S. *Nat. Mater.* **2005**, *4*, 335–339.
doi:10.1038/nmat1349

2. Bogani, L.; Wernsdorfer, W. *Nat. Mater.* **2008**, *7*, 179–186.
doi:10.1038/nmat2133

3. Gambardella, P.; Stepanow, S.; Dmitriev, A.; Honolka, J.; de Groot, F. M. F.; Lingenfelder, M.; Sen Gupta, S.; Sarma, D. D.; Bencok, P.; Stanescu, S.; Clair, S.; Pons, S.; Lin, N.; Seitsonen, A. P.; Brune, H.; Barth, J. V.; Kern, K. *Nat. Mater.* **2009**, *8*, 189–193.
doi:10.1038/nmat2376

4. Mannini, M.; Pineider, F.; Sainctavit, P.; Danieli, C.; Otero, E.; Sciancalepore, C.; Talarico, A. M.; Arrio, M.-A.; Cornia, A.; Gatteshi, D.; Sessoli, R. *Nat. Mater.* **2009**, *8*, 194–197. doi:10.1038/nmat2374

5. Urdampilleta, M.; Klyatskaya, S.; Cleziou, J.-P.; Ruben, M.; Wernsdorfer, W. *Nat. Mater.* **2011**, *10*, 502–506.
doi:10.1038/nmat3050

6. Candini, A.; Klyatskaya, S.; Ruben, M.; Wernsdorfer, W.; Affronte, M. *Nano Lett.* **2011**, *11*, 2634–2639. doi:10.1021/nl2006142

7. Vincent, R.; Klyatskaya, S.; Ruben, M.; Wernsdorfer, W.; Balestro, F. *Nature* **2012**, *488*, 357–360. doi:10.1038/nature11341

8. Schwöbel, J.; Fu, Y.; Brede, J.; Dilullo, A.; Hoffmann, G.; Klyatskaya, S.; Ruben, M.; Wiesendanger, R. *Nat. Commun.* **2012**, *3*, 953. doi:10.1038/ncomms1953

9. Ganzhorn, M.; Klyatskaya, S.; Ruben, M.; Wernsdorfer, W. *Nat. Nanotechnol.* **2013**, *8*, 165–169. doi:10.1038/nnano.2012.258

10. Katoh, K.; Komeda, T.; Yamashita, M. *Dalton Trans.* **2010**, *39*, 4708–4723. doi:10.1039/b926121d

11. Klar, D.; Klyatskaya, S.; Candini, A.; Krumme, B.; Kummer, K.; Ohresser, P.; Corradini, V.; de Renzi, V.; Biagi, R.; Joly, L.; Kappler, J.-P.; del Pennino, U.; Affronte, M.; Wende, H.; Ruben, M. *Beilstein J. Nanotechnol.* **2013**, *4*, 320–324. doi:10.3762/bjnano.4.36

12. Palik, E. D. *Handbook of Optical Constants of Solids II*; Academic Press, 1991.

13. Gordan, O. D.; Friedrich, M.; Zahn, D. R. T. *Org. Electron.* **2004**, *5*, 291–297. doi:10.1016/j.orgel.2004.10.001

14. Gordan, O. D.; Friedrich, M.; Zahn, D. R. T. *Thin Solid Films* **2004**, *455*, 551–556. doi:10.1016/j.tsf.2003.11.209

15. Yoshino, K.; Lee, S. B.; Sonoda, T.; Kawagishi, H.; Hidayat, R.; Nakayama, K.; Ozaki, M.; Ban, K.; Nishizawa, K.; Ohta, K.; Shirai, H. *J. Appl. Phys.* **2000**, *88*, 7137–7143. doi:10.1063/1.1316050

16. Bräuer, B.; Fronk, M.; Lehmann, D.; Zahn, D. R. T.; Salvan, G. *J. Phys. Chem. B* **2009**, *113*, 14957–14961. doi:10.1021/jp9073224

17. Margheriti, L.; Chiappe, D.; Mannini, M.; Car, P.-E.; Sainctavit, P.; Arrio, M.-A.; de Mongeot, F. B.; Cezar, J. C.; Piras, F. M.; Magnani, A.; Otero, E.; Caneschi, A.; Sessoli, R. *Adv. Mater.* **2010**, *22*, 5488–5493. doi:10.1002/adma.201003275

18. Malavolti, L.; Poggini, L.; Margheriti, L.; Chiappe, D.; Graziosi, P.; Cortigiani, B.; Lanzilotto, V.; Buatier de Mongeot, F.; Ohresser, P.; Otero, E.; Choueikani, F.; Sainctavit, Ph.; Bergenti, I.; Dediu, V. A.; Mannini, M.; Sessoli, R. *Chem. Commun.* **2013**, *49*, 11506–11508. doi:10.1039/c3cc46868b

19. Khrishnakumar, K. P.; Menon, C. S. *Mater. Lett.* **2001**, *48*, 64–73. doi:10.1016/S0167-577X(00)00281-0

20. Siles, P. F.; Bof Bufon, C. C.; Grimm, D.; Jalil, A. R.; Mende, C.; Lungwitz, F.; Salvan, G.; Zahn, D. R. T.; Lang, H.; Schmidt, O. G. *Org. Electron.* **2014**, *15*, 1432–1439. doi:10.1016/j.orgel.2014.04.004

21. Vervacke, C.; Bof Bufon, C. C.; Thurmer, D. J.; Siles, P. F.; Schmidt, O. G. *Anal. Chem.* **2012**, *84*, 8399–8406. doi:10.1021/ac301959e

22. Alexeev, A.; Loos, J.; Koetse, M. M. *Ultramicroscopy* **2006**, *106*, 191–199. doi:10.1016/j.ultramic.2005.07.003

23. MacDonald, G. A.; Veneman, P. A.; Placenica, D.; Armstrong, N. R. *ACS Nano* **2012**, *6*, 9623–9636. doi:10.1021/nn303043y

24. Hendriksen, B. L. M.; Martin, F.; Qi, Y.; Mauldin, C.; Vukmirovic, N.; Ren, J.; Wormeester, H.; Katan, A. J.; Altoe, V.; Aloni, S.; Fréchet, J. M. J.; Wang, L.-W.; Salmeron, M. *Nano Lett.* **2011**, *11*, 4107–4112. doi:10.1021/nl202720y

25. Walter, M. G.; Rudine, A. B.; Wamser, C. C. *J. Porphyrins Phthalocyanines* **2010**, *14*, 759–792. doi:10.1142/S1088424610002689

26. Reid, O. G.; Munechika, K.; Ginger, D. S. *Nano Lett.* **2008**, *8*, 1602–1609. doi:10.1021/nl080155l

27. Horowitz, G.; Fichou, D.; Peng, X.; Delannoy, P. *J. Phys. France* **1990**, *51*, 1489–1499. doi:10.1051/jphys:019900510130148900

28. Bof Bufon, C. C.; Vervacke, C.; Thurmer, D. J.; Fronk, M.; Salvan, G.; Lindner, S.; Knupfer, M.; Zahn, D. R. T.; Schmidt, O. G. *J. Phys. Chem. C* **2014**, *118*, 7272–7279. doi:10.1021/jp409617r

29. Yakuphanoglu, F.; Caglar, M.; Caglar, Y.; Ilican, S. *Synth. Met.* **2010**, *160*, 1520–1523. doi:10.1016/j.synthmet.2010.05.013

30. WVASE®, J.A. Woollam Co. Inc.: Lincoln, NE, 2014.

31. Horcas, I.; Fernández, R.; Gómez-Rodríguez, J. M.; Colchero, J.; Gómez-Herrero, J.; Baro, A. M. *Rev. Sci. Instrum.* **2007**, *78*, 013705. doi:10.1063/1.2432410

32. Gwyddion, 2.38; Klapetek, P.; Nečaš, D.; Brno, Czech Republic, 2014.
License and Terms

This is an Open Access article under the terms of the
Creative Commons Attribution License
([http://creativecommons.org/licenses/by/2.0](http://creativecommons.org/licenses/by/2.0)), which
permits unrestricted use, distribution, and reproduction in
any medium, provided the original work is properly cited.

The license is subject to the *Beilstein Journal of*
*Nanotechnology* terms and conditions:
([http://www.beilstein-journals.org/bjnano](http://www.beilstein-journals.org/bjnano))

The definitive version of this article is the electronic one
which can be found at:
[doi:10.3762/bjnano.5.215](http://doi.org/10.3762/bjnano.5.215)