Thomas G.
Dane
a,
Julia E.
Bartenstein
a,
Beatrice
Sironi
a,
Benjamin M.
Mills
a,
O.
Alexander Bell
a,
J.
Emyr Macdonald
b,
Thomas
Arnold
c,
Charl F. J.
Faul
*a and
Wuge H.
Briscoe
*a
aSchool of Chemistry, University of Bristol, Cantock's Close, Bristol BS8 1TS, UK. E-mail: charl.faul@bris.ac.uk; wuge.briscoe@bris.ac.uk
bSchool of Physics and Astronomy, Cardiff University, Cardiff CF24 3AA, UK
cDiamond Light Source Ltd., Diamond House, Harwell Science and Innovation Campus, Didcot, Oxfordshire, UK
First published on 12th August 2016
The influence of processing conditions on the thin film microstructure is a fundamental question that must be understood to improve the performance of solution-processed organic electronic materials. Using grazing-incidence X-ray diffraction, we have studied the structure of thin films of a tetra(aniline)-surfactant complex prepared by drop-casting from five solvents (hexane, chloroform, tetrahydrofuran, dichloromethane and ethanol), selected to cover a range of polarities. We found that the structure, level of order and degree of orientation relative to the substrate were extremely sensitive to the solvent used. We have attempted to correlate such solvent sensitivity with a variety of solvent physical parameters. Of particular significance is the observation of a sharp structural transition in the thin films cast from more polar solvents; such films presented significantly greater crystallinity as measured by the coherence length and paracrystalline disorder parameter. We attribute this higher structural order to enhanced dissociation of the acid surfactant in the more polar solvents, which in turn promotes complex formation. Furthermore, the more polar solvents provide more effective screening of (i) the attractive ionic interaction between oppositely charged molecules, providing greater opportunity for dynamic reorganisation of the supramolecular aggregates into more perfect structures; and (ii) the repulsive interaction between the positively charged blocks permitting a solvophobic-driven aggregation of the aromatic surfaces during solvent evaporation.
One important factor in the self-assembly behaviour of such organic molecules is the processing solvent. The cost, toxicity and properties of the solvent are also relevant commercial and environmental considerations. The processing solvent has been shown to influence the microstructure and charge carrier mobility in devices based on the well-studied semiconducting polymer poly(3-hexylthiophene) (P3HT).5–11 Some investigations, for example, suggest that films spin-cast from high boiling point solvents demonstrate a higher degree of structural order and orientation than from lower boiling point solvents.8,9 This observation was ascribed to the longer evaporation time, which permitted the formation of thermodynamically stable structures. It was also found that the P3HT aromatic planes were parallel to the substrate when drop-cast from chloroform, but perpendicular when cast from dichloromethane (DCM); such induced differences would have obvious consequences on the direction of charge transport through the film.10
The bulk-phase supramolecular structure of poly(aniline)—a conducting polymer offering environmental stability and switchable conductivity and optical transitions—is strongly dependent on the processing conditions;12–19 however, there have been very few studies on such dependencies for thin films. The use of oligomers of poly(aniline) offers a means to overcome the polydispersity and structural disorder inherent to such polymers,20 and can induce high levels of supramolecular order via an ionic self-assembly (ISA) route.21,22 Such model compounds therefore provide opportunities to understand physical mechanisms and control subtle factors that govern structure20 and functionality.21
We have previously investigated the thin self-assembly behaviour of aniline oligomers when deposited by drop casting. This process, of relevance to widely used thin film deposition processes, proved to be a facile route to highly ordered thin films. We found that the oligomer thin film self-assembly behaviour depends strongly on the molecular architecture and thermal treatment.25,26
Of particular interest is the TANI(BEHP)0.5 oligomer–surfactant complex, shown in Fig. 1(a), which is formed from phenyl/phenyl end-capped tetra(aniline) (Ph/Ph TANI) and two equivalents of the acid surfactant dopant bis(ethyl hexyl) hydrogen phosphate (BEHP). This system forms a highly ordered bilayer-type structure as represented in Fig. 1(b). In the bulk phase, TANI(BEHP)0.5 forms a rhombohedral 2D unit cell with lattice parameters a = 2.35 nm, b = 2.7 nm and γ = 80°.21 Thin films drop-cast from tetrahydrofuran (THF) adopt a hexagonal lattice (a = b = 2.53 nm, γ = 120°), whereby the bilayers preferentially orient parallel to the substrate.25,26
Grazing-incidence X-ray diffraction (GIXD) is a valuable tool in ascertaining structural information of soft matter thin films.27 In combination with a 2D detector, the structure normal to the surface (out-of-plane, qz) and within the surface plane (in-plane, qxy) can be simultaneously probed (as represented pictorially in Fig. 1(c)). Here we investigate the effect of solvent used for drop-casting on the resulting structure and degree of order in TANI(BEHP)0.5 thin films using GIXD. We determine the level of order within the films by considering the broadening of diffraction peaks and relate these structural characteristics to physical properties of the solvents (such as boiling point, vapour pressure and polarity). We have considered five solvents (including the previously studied THF as a control), and have summarised their relevant physical properties in Table 1. It is expected that the results from this initial study will provide fundamental insight into the role of the solvent, and thus stimulate further investigations into the influence of solvent properties on structure formation, self-assembly and function.
Solvent | ε 23 | Hansen solubility parameters (Pa1/2)24 | b.p.23 (°C) | v.p.23 (kPa) | d 23 (g mL−1) | |||
---|---|---|---|---|---|---|---|---|
δ Total | δ D | δ P | δ H | |||||
Hexane | 1.89 | 14.9 | 14.9 | 0.0 | 0.0 | 68.7 | 20.2 | 0.661 |
Chloroform | 4.81 | 19.0 | 17.8 | 3.1 | 5.7 | 61.2 | 26.2 | 1.479 |
THF | 7.52 | 19.4 | 16.8 | 5.7 | 8.0 | 65.0 | 21.6 | 0.883 |
DCM | 8.93 | 20.3 | 18.2 | 6.3 | 6.1 | 40.0 | 58.2 | 1.327 |
Ethanol | 25.3 | 26.5 | 15.8 | 8.8 | 19.4 | 78.3 | 7.87 | 0.789 |
Films were produced by drop-casting as follows: a droplet of solution (40 μL) was placed on the substrate and allowed to evaporate for one hour in a closed atmosphere saturated with vapour from the solvent in question (by exposure to an open vial of the solvent (∼5 mL)). The films were then examined by GIXD within 24 hours of film preparation.
GIXD measurements were performed on beamline I07, Diamond Light Source, UK. The X-ray beam energy was E = 10 keV (wavelength λ = 1.24 Å) and the beam size was ∼300 × 300 μm (FWHM) with an approximately Gaussian intensity profile. An incident angle of αi = 0.36° (∼2 × αc, the critical angle) was used, as this allowed for complete illumination of the film with minimal substrate-reflected beam which could complicate analysis (see Section S3 in the ESI† for further discussion). Data were collected on a PILATUS 2M detector (DECTRIS) at a distance of 340.02 mm from the sample, which was calibrated by silicon powder in transmission mode. Samples were enclosed in a helium-filled chamber,29 mounted on a hexapod attached to a 2+3 circle Huber diffractometer.
The GIXD scattering geometry is shown in Fig. 1(c). GIXD data were reduced using a custom-made software package written in IgorPro (Wavemetrics®).26 The diffraction patterns were transformed from the raw images into reciprocal space maps, i.e., scattered intensity as a function of the in-plane (qxy) and out-of-plane (qz) scattering vectors. Diffraction features along the out-of-plane direction reveals information about structuring normal to the surface, whilst in-plane features are related to the ordering within the plane of the film. Out-of-plane line profiles (intensity vs. qz) were extracted by integrating a sector of data at χ = 90° (where χ is the azimuthal angle relative to the surface plane, cf.Fig. 1), with an angular width of Δχ = 25°. In-plane line profiles (intensity vs. qz) were extracted similarly at χ = 5°, with an angular width of Δχ = 5°. The larger angular integration width for the out-of-plane profiles was chosen to account for the splitting of q-space at higher qz values. Azimuthal line profiles (intensity vs. χ) were assessed to determine the degree of crystallite orientation relative to the surface and were extracted at a fixed q ∼ 0.55 Å−1 with a radial width of Δq = 0.05 Å−1. The integration regions used to extract out-of-plane, in-plane and azimuthal line profiles were the same for each film and are shown in Fig. S2 in the ESI.†
There are a number of features in these data to be highlighted: firstly, for each of the films examined, five orders of the (h0) family of lamellar reflections were present in the GIXD patterns (and out-of-plane line profiles, which are indexed with purple labels), suggesting the formation of the TANI(BEHP)0.5 lamellar structure with the (10) planes preferentially oriented parallel to the substrate. For all systems, this lattice spacing was 2.27–2.28 nm, which is in good agreement with one of the three polymorphs found in the previous study (2.29 nm),26 except for the hexane-cast film, which presented d(10) = 2.21 nm. There is a diffuse isotropic ring in the GIXD patterns (and present as very broad peaks in the line profiles) at q ∼ 1.4 Å−1, which arises from the disordered alkyl chains of the BEHP dopant. The scattering intensity of this feature relative to that of the lamellar reflections varies between samples. Lastly, an important feature in the diffraction patterns is the reflection arising from the π-stacking interaction of the TANI units. When present, this feature occurs at q ∼ 1.75 Å−1 (lattice spacing d ∼ 0.36 nm) and is most apparent in the in-plane line profiles, which indicates that these planes preferentially orient normal to the surface. This observation is consistent with the (10) preferential orientation relative to the surface. These features in Fig. 2 are labelled with orange arrows on the in-plane line profiles. Fig. 2(e) shows the GIXD data from the THF-cast thin film, which is consistent with our previous studies.26 Despite the presence of common features, there were a number of crucial differences between the data for each of the five solvents used for film casting.
In contrast to the other four films, the GIXD pattern of the hexane-cast film showed isotropic Debye–Scherrer diffraction rings with little azimuthal intensity variation (cf.Fig. 2(a)), indicating that the crystallites were randomly oriented relative to the substrate plane. A reflection at d = 0.36 nm, present only in the in-plane line profile, indicates that some π-stacking was present in the film, with these domains preferentially oriented perpendicular to the substrate. Furthermore, the (11), (21) and (31) reflections, absent from the films cast from other solvents, were observed here for the hexane-cast film.
The chloroform- and THF-cast films presented very similar features in their GIXD patterns (cf.Fig. 2(c) and (e)), and line profiles ((d) and (f)). Both films displayed hexagonal symmetry, evident from the presence of the (02) reflections at 60° to the (20) reflection (cf. azimuthal profile in Fig. 3). The main difference between these two films was the presence of a π-stacking reflection (d = 0.37 nm) for the THF-cast film, which was not present for the chloroform-cast film.
![]() | ||
Fig. 3 Azimuthal integration about the (20) reflection for each of the TANI(BEHP)0.5 films. Dotted lines represent Lorentzian fits to the intensity maxima at χ = 90°. |
Both the DCM- and ethanol-cast films also displayed the typical lamellar reflections with a relatively high degree of preferential orientation of crystallites parallel to the substrate. The GIXD patterns for both films showed many additional reflections, evident from the peaks in the range 1.25 Å−1 < qxy < 2.5 Å−1 in the in-plane line profiles (highlighted with orange stars in Fig. 2(h) and (j)). These sharp reflections extended to the limit of observable reciprocal space, suggesting highly crystalline structures. This observation is consistent with the AFM images showing well-defined crystallites (Fig. S1(i) and (j) in the ESI†). It has not yet been possible to interpret the detailed structure of this phase, though our investigations are ongoing. No hexagonal symmetry was observed in either the DCM- or ethanol-cast films. For both samples, we observed π-stacking reflections of lattice spacing d = 0.35 nm.
Note that for the DCM- and ethanol-cast films, the reflections were split into two sets of closely overlapping peaks, visible in both the 2D patterns and the line profiles, with expanded views shown in Fig. S3 in the ESI.† We attribute it to two possible explanations.‡ First, polymorphism could be present (two structures with differing lattice parameters, well-known for organic thin films). Alternatively, the large footprint of the X-ray beam caused smearing of the scattering pattern from isolated thicker regions of the film. These two out-of-plane reflection series are referred to below with the notation (I) and (II).
Solvent | Lamellar spacing d (nm) | Coherence length La (nm) | Disorder parameter g (%) | Azimuthal FWHM Δχ(20) (°) |
---|---|---|---|---|
Hexane | 2.21 | 22.4 ± 0.3 | 2.4 ± 0.3 | — |
Chloroform | 2.27 | 23.4 ± 0.7 | 2.4 ± 0.3 | 4.1 ± 0.3 |
THF | 2.28 | 23.7 ± 0.3 | 2.6 ± 0.4 | 3.5 ± 0.2 |
DCM | ||||
(I) | 2.30 | 69.3 ± 2.6 | 2.0 ± 0.6 | — |
(II) | 2.26 | 66.5 ± 1.7 | 1.6 ± 0.6 | — |
Average | 2.28 | 67.9 ± 2.2 | 1.8 ± 0.6 | 18.8 ± 0.4 |
Ethanol | ||||
(I) | 2.29 | 66.6 ± 1.1 | 1.7 ± 0.4 | — |
(II) | 2.25 | 86.6 ± 1.6 | 1.7 ± 0.3 | — |
Average | 2.27 | 76.6 ± 1.4 | 1.7 ± 0.4 | 19.4 ± 0.2 |
The FWHM Δq and d-spacing values for each of the TANI(BEHP)0.5 samples cast from different solvents were determined from fitting the out-of-plane line profiles (Fig. 2) using the MultiPeak fitting routines in IgorPro (Wavemetrics) (fitted data are shown in Fig. S5 in the ESI†). The paracrystalline disorder plots are shown in Fig. 4, with the data points shown as open symbols and the fits by dashed lines.¶ Note that the data for the second phases (II) of the DCM- and ethanol-cast films are shown as filled symbols and the fits as solid lines. These data are also plotted separately for clarity in Fig. S6 in the ESI.† The full numerical fitting results are summarised in Table 2.
The paracrystalline disorder plots reveal a clear difference between the level of order in the films cast from more polar solvents (DCM and ethanol) and those cast from less polar solvents (hexane, chloroform and THF). The fitting for the polar solvents gave a lower intercept on the y-axis and a smaller gradient, indicating both larger domain sizes (∼70 nm vs. ∼23 nm) and less paracrystalline disorder (∼1.7% vs. ∼2.4%) in these films compared with those of the less polar solvents. We now proceed to relate the origin of these structural differences with the physical properties of the solvents.
It is desirable to correlate structural characteristics with certain solvent physical parameters, which would offer clear guidance for solvent choice and processing conditions. Previously, the solvent boiling point, vapour pressure and dielectric constant have been used to explain observed variations in the thin film structure.5–11 Hildebrand and Scott proposed a solubility parameter δ defined as the square root of the cohesive energy density given by33
![]() | (1) |
![]() | (2) |
There appears to be no correlation between the level of structural order and the boiling point of the solvent (cf. Fig. S7(a) in the ESI†). Whilst previous studies on P3HT have reported correlations between the boiling point and thin film structure,8,9 the similarity between the diffraction patterns for the DCM- and ethanol-cast films (boiling points 40 °C and 77 °C, respectively) suggests that, as a parameter reflecting the total intermolecular attractions between the solvent molecules, the solvent boiling point cannot explain the structural differences as a result of TANI(BEHP)0.5 self-assembly in different solvents.
It is well known that the residual patterns from the evaporative drying process of a sessile drop on a surface depend on the evaporation rate.35,36 However, the vapour pressure (which is a measure of the evaporation rate of the solvents) does not show any correlation with the thin film structure (cf. Fig. S7(b) in the ESI†). This is most apparent in the DCM and ethanol cast films, which have significantly different vapour pressures (58.2 and 7.87 kPa, respectively) yet similar thin film structures. The microscopic morphology and topography of the residual patterns depend on the capillary and Marangoni flows which are steered by evaporation and thus depend on the evaporation rate. Indeed, our AFM imaging (Fig. S1 in the ESI†) has revealed different topographies of the films drop cast from different solvents. However, the GIXS results examine the structure on a different length scale (<5 nm) reflecting the intermolecular packing. Thus, the nanoscopic structural order we observe is not affected by the different evaporation rates, whist the residual patterns and topography on a larger length scale are.
The dispersion and hydrogen-bonding HSPs also do not correlate with the film structure (cf. Fig. S9 in the ESI†). However, the polarity of the solvent measured both by the dielectric constant (cf. Fig. S8(a) in the ESI†) and the polar HSP correlate more closely with the resulting film order as shown Fig. 5. Note that the average values of La and g have been plotted for the (I) and (II) phases of the DCM- and ethanol-cast films. The structural order within the film also correlates with the total solubility parameter (cf. Fig. S8(b), though not δ–δP as shown in Fig. S10, ESI†). However, given the dispersion and hydrogen-bonding HSPs do not seem to play a significant role, it can be inferred that the polarity of the solvent is here the dominating factor.
The plot of coherence length and paracrystalline disorder parameter as a function of the polar HSP shown in Fig. 5 indicates that the ionic self-assembly of TANI(BEHP)0.5 is promoted in a more polar environment and that there exists a critical polarity, above which the resulting level of structural order in the films is significantly greater than below this polarity. This critical polarity occurs at around 6 Pa1/2, and we should note that this value is specific to the chemical system investigated.
A variety of dynamic processes and interactions are present in solutions of TANI and BEHP. These can be summarized as follows: (a) formation of the charged species via dissociation of the acid surfactant (and protonation of the basic TANI molecule): TANI + 2BEHP ⇌ TANI2+ + 2BEHP−; (b) ionic interactions between the charged species: TANI2+ + 2BEHP− ⇌ TANI(BEHP)2; (c) solubility governed by solvent–molecule interactions for the three equilibrium states (uncharged TANI and BEHP, charged TANI2+ and BEHP− and charge-neutralized TANI(BEHP)2); (d) aggregation driven by van der Waals, π-stacking and hydrogen bonding interactions; (e) entropy of aggregation. The relative importance of these processes (in addition to being dependent on the solvent properties) is also dependent on concentration and are thus temporally dynamic during the evaporation process.
Whilst we have only considered five solvents, we hypothesize that the clear correlation between the degree of structural order and solvent polarity (and lack of correlation with any other solvent properties) is due to an increased driving force for self-assembly in more polar solvents, which can be explained as follows. In more polar solvents (higher dielectric constant), the degree of dissociation of the acid surfactant is greater. Whilst the strength of the ionic interaction between the charged species is weaker in more polar solvents (due to greater charge-screening of this interaction), the increased population of the charged species promotes the formation of the TANI(BEHP)0.5 complex. In the more polar solvents, the weaker ionic interactions between the ionic species means that there is a greater opportunity for dynamic reorganisation of the charged species into more thermodynamically stable structures during aggregation in the final stages of evaporation, leading to an increase in the overall structural order.
Furthermore, the charge screening also reduces repulsive interactions between charged nitrogen atoms on the TANI units, which facilitates aggregation through π-stacking. This solvophobic interaction is particularly favourable in more polar solvents so as to minimise contact of the hydrophobic aromatic planes with the polar medium, and may serve to direct further self-assembly. This hypothesis is supported by the observation of both enhanced crystallinity in the π-stacking direction and larger bilayer stacking domains for films cast from the more polar solvents, DCM and ethanol.
Footnotes |
† Electronic supplementary information (ESI) available. See DOI: 10.1039/c6cp05221e |
‡ We rule out the possibility of the incoming X-ray beam reflecting from the underlying substrate and acting as a second illumination source (see ESI† for full discussion). |
§ The “missing” data can be obtained by setting the incident angle such that the specular condition coincides with the Bragg reflection,7 which was not possible due to limitations in the experimental configuration. |
¶ Note that we plot only the orders h = 2 to 5 because at the first order reflection, the I and II phases for the DCM- and ethanol-cast films are too closely overlapping to separately resolve. |
This journal is © the Owner Societies 2016 |