Nobutsugu Hamamotoa,
Hiromitsu Sonodaa,
Michinori Sumimotob,
Kenji Horib and
Hitoshi Fujimoto*a
aDepartment of Chemistry, Faculty of Science, Kumamoto University, Kurokami 2-39-1, Kumamoto 860-8555, Japan. E-mail: fuji@aster.sci.kumamoto-u.ac.jp
bGraduate School of Science and Engineering, Yamaguchi University, Tokiwadai, Ube 755-8611, Japan
First published on 27th January 2017
Lead(II) phthalocyanine (PbPc) crystallises into two modifications of monoclinic and triclinic systems. In order to investigate the solid state electronic structures, we calculated three model dimers and two model trimers of PbPc using the M06 method of density functional theory (DFT). The geometries of these models for the crystal polymorphs were optimised, and the stability of each model was discussed along with the Gibbs free energy when comparing with the case of the triclinic system of tin(II) phthalocyanine. The optimised geometries of these models simulated well the molecular structures and orientations in the crystalline solids of PbPc. The time-dependent DFT (TD-DFT) calculation has been applied to these model dimers and trimers in order to investigate the excited states. The results also reproduced well the observed UV-vis spectra of PbPc in the solid state.
Among many MPc compounds, lead(II) phthalocyanine (PbPc) has a peculiar structure with C4v symmetry like a shuttlecock. It has been reported that there exists two polymorphs in PbPc crystals, the monoclinic system with space group P21/b5 and the triclinic one with P.6 Thin films with both crystalline polymorphs showed different sensitivity against air, nitrogen oxide, and water vapour.7 PbPc has been also used as non-linear optical materials,8 organic solar cells,9 and so on. Tin(II) phthalocyanine (SnPc) also has a C4v shuttlecock-shaped structure; however, the triclinic polymorph has been only reported for SnPc10,11 up to now to our knowledge.
Density functional theory (DFT) methods have been applied to investigate various physical properties of several MPc compounds.12–19 The time-dependent DFT (TD-DFT) calculations have been especially used to discuss the UV-vis spectra of them. These works showed that DFT methods are efficient, reliable and effective to investigate the molecular and electronic structures of MPcs for the present. In the previous report, we successfully investigated physical properties for the triclinic solid of SnPc with these methods using the dimer models.20
In this study, we extended the same approach as SnPc to investigating the electronic structures of two polymorphic PbPc solids. Three types of dimers were extracted from the two crystal systems; Type 1 from the monoclinic polymorph, Types 2 and 3 from the triclinic one (see Scheme 1). The geometries of these dimers were optimised using the DFT method, and the results were compared to the observed molecular structures in solids. We conducted an investigation into electronic structures of both ground and excited states for crystals with using these model dimers. The results of the TD-DFT calculations were examined carefully with comparing to the observed electronic absorption spectra for thin solid films. In order to check the propriety of the dimer model, we also compared the same calculation results for the model trimers of Types 4 and 5 (see Scheme 1) for monoclinic and triclinic polymorphs, respectively, as those for the dimers. Through these investigations, we will give an insight into the origin of polymorphs in SnPc and PbPc.
![]() | ||
Scheme 1 Structures of the model dimers and trimers of PbPc. The distances between Pb atoms are referred from ref. 5 (Type 1) and 6 (Types 2 and 3) for monoclinic and triclinic crystals, respectively. |
In order to give insight into the electronic structures of PbPc polymorphic solids, we optimised the geometries of the model dimers (Types 1–3) and trimers (Types 4 and 5) in Scheme 1 under the following conditions: (1) the structural parameters of each PbPc moiety were fixed as those of the PbPc monomer, which were optimised under constraint to keep C4v symmetry using the parameter shown in Fig. 1. (2) The distance between the PbPc moieties was optimised with the dimer of Type 1 to keep C4v symmetry, and the distance and the tilt angle between them were optimised simultaneously with those of Types 2 and 3 under constraint to keep C2h symmetries. (3) We calculated the trimers under constraint to keep C4v and C1 symmetries for Types 4 and 5, respectively. The two distances between the PbPc moieties were optimised independently with the trimer of Type 4, and the distances and tilt angles between them were optimised with the trimer of Type 5.
All these models had vibration modes with small negative frequency. Our purpose is to reproduce molecular orientations in PbPc solids for investigating their electronic structures; therefore, we proceeded with this study using these geometries. The differences of the total energy between the models and the corresponding components are listed in Table S1 of the ESI.†
Here, we would like to compare the optimised structures of PbPc and SnPc with the M06 and M06-2X methods. The Pc rings of both compounds have a saucer-shaped structure, where the central metal atom sticks out of the molecule on the opposite side of its curvature. The calculation results showed that the Mulliken charge of the Pb atom in PbPc is estimated to be +1.1 e and +1.2 e by the M06 and M06-2X methods, respectively, which are similar to the value of +0.9 e and +1.0 e for the Sn atom of SnPc by the M06 and M06-2X methods, respectively. It should be noted, as pointed out in the previous work,12 that there exist little differences of the bond lengths and bond angles in the region of the Pc rings between SnPc and PbPc, where the curvature of the Pc ring differs just a little. Contrary to this similarity in the atomic charges and molecular structures, the dipole moments were quite different between PbPc and SnPc, where both methods gave the same values of 1.0 and −0.6 Debye for PbPc and SnPc, respectively, along the four-fold rotational axis. This might be originated from the small difference of about 0.2 e in the atomic charges of the central metal atoms and in the atomic distances between the central metal and its neighbour nitrogen atoms. The distances were estimated to be 2.340 and 2.334 Å for PbPc by the M06 and M06-2X methods, respectively, where those of SnPc were to be 2.272 and 2.270 Å by these two methods, respectively. These values were similar to the value of 2.279 Å obtained by the previous work.12
![]() | ||
Fig. 2 Optimised geometries of the model dimers of PbPc. The relative Gibbs free energies (kcal mol−1) to the sum of the monomers are shown in parentheses. |
There is no report on the SnPc dimer similar to Type 1; therefore, we additionally calculated for the Type 1 dimer of SnPc using the same method with this study. The relative Gibbs free energy was interestingly estimated to be +2.2 kcal mol−1. This value shows clearly that the model dimer of the monoclinic crystal is unstable as compared, at least, to those of the triclinic crystal for SnPc.
The distances between two Pb atoms were evaluated to be 3.656, 7.976 and 6.195 Å for the optimised geometries of Types 1–3, respectively. These values were consistent with the corresponding experimental values of 3.73 Å for the monoclinic crystal5 and of 8.26 and 6.52 Å for the triclinic one.6
The M06-2X method gave nearly similar results to the M06 one in stabilities of the dimers. The differences between the Gibbs free energies of the dimers and two monomers were estimated to be −5.9, −8.3 and −12.9 kcal mol−1 for Types 1–3, respectively. The distances between Pb atoms were evaluated to be 3.748, 7.815 and 6.278 Å for Types 1–3, respectively, for the optimised geometries with the M06-2X method. These values are similar to those obtained by the M06 method; therefore, we chose the M06 method in the following discussion for the simple reason to compare the present results with those reported previously for SnPc.12
Fig. 3 illustrates the optimised structures of the two model trimers, Types 4 and 5. The relative Gibbs free energies of these trimers were much lower than the sum total of three monomers by 15.2 and 17.8 kcal mol−1 for the trimers of Types 4 and 5, respectively. It should be noted that the relative energies of these two trimers become closer to each other as compared to those of the dimers. It was also shown that the Type 5 trimer is less stable about 5 kcal mol−1 than the corresponding trimer of SnPc, for which the relative Gibbs free energy was estimated to be −23.4 kcal mol−1.20 In addition, we also performed the calculations for the similar SnPc trimer to the Type 4 one, and the relative Gibbs free energy was estimated to be −0.9 kcal mol−1. These results showed clearly that the PbPc molecules are much stable in a columnar structure of the monoclinic crystal, contrary to SnPc. This is consistent with the experimental facts that PbPc crystallises in both triclinic and monoclinic crystal systems, but the triclinic crystal system was only reported on SnPc up to now. This remarkable difference between SnPc and PbPc might be originated from the small differences in the curvature of the molecules and the projection of the metal atoms. It might also be thought that a charge distribution on the molecule will affect the stable arrangement for each pair of the PbPc molecules in solids.
![]() | ||
Fig. 3 Optimised geometries of the model trimers of PbPc. The relative Gibbs free energies (kcal mol−1) to the sum of the monomers are shown in parentheses. |
The distances between two Pb atoms of the optimised geometries were evaluated to be 3.656 and 3.654 Å for the Type 4 trimer and 8.034 and 6.272 Å for the Type 5 one. It would be mentioned that these values for the Type 5 trimer become much closer to the distances of 8.26 and 6.52 Å in the triclinic crystal6 as compared to those calculated for the dimers. In contrast with this, the distances between the Pb atoms for the Type 4 trimer were calculated to be nearly the same value to that of the dimer and the value of 3.73 Å observed in the monoclinic crystal.5 These facts, in conjunction with the relative Gibbs free energy, might indicate that the dimer is sufficient to investigate the molecular and electronic structures in the monoclinic system theoretically. On the other hand, it would be concluded properly from the crystal structure that the triclinic crystal of PbPc requires at least the trimer for investigation. In passing, the differences between the observed and calculated values may be due to the effect of the crystal packing.
Fig. 4 indicates the energy diagrams of the monomer, the three dimers and the two trimers in the region of the highest occupied (HOMO) and the lowest unoccupied (LUMO) molecular orbitals. The orbital symmetries are labelled under the C4v symmetry group for the monomer, the Type 1 dimer and the Type 4 trimer, under the C2h symmetry group for the dimers of Types 2 and 3, and under the C1 symmetry group for the Type 5 trimer. The energy diagrams of the dimers basically consisted of two orbitals in pairs, which were constructed of interactions between two moieties in bonding and anti-bonding ways. In the case of the trimers, there existed the non-bonding type molecular orbital in addition to each pair of bonding and anti-bonding type orbitals. The energy values of these molecular orbitals are compiled in Table S3 of the ESI,† where the corresponding molecular orbitals are indicated in the same colour. The atomic contributions to the several molecular orbitals near the frontier orbitals are shown schematically in Fig. S3–S7 of the ESI† for the model dimers and trimers of Types 1–5, respectively.
The HOMO of the PbPc monomer has a large contribution from carbon atoms of the inner C8N8 ring and the next occupied orbital to the HOMO (HOMO − 1) has a strong contribution from s-type atomic orbital of a Pb atom as shown in Fig. 4. It would be reasonable to think that the energy splitting between the bonding and anti-bonding type orbitals depends on the strength of the interaction between PbPc moieties of the dimers. From the calculation results, the energy differences between the HOMO and HOMO − 1 were estimated to be 0.61, 0.09 and 0.34 eV for the dimers of Types 1–3, respectively. It has been reported that the values for SnPc are estimated to be 0.30 and 0.34 eV for the Type 2 (concave) and Type 3 (convex) dimers, respectively.20 In the case of the trimers, these energy differences between the bonding and anti-bonding type orbitals were evaluated to be 0.92 and 0.36 eV for the trimers of Types 4 and 5, respectively, which were a little larger than those of the corresponding dimers. The value of the Type 5 PbPc trimer was a little smaller than that of the SnPc trimer, which has been reported as to be 0.43 eV.20
We also calculated the ionisation potentials and the electron affinities by the ΔSCF method, and the results are compiled in Table 1. For comparison, the values for SnPc are also shown in parenthesis, where the ionisation potentials are based on the previous work20 for the monomer, the dimers of Types 2 and 3, and the trimer of Type 5. It has been reported35 that the bilayer film of PbPc shows two bands split with about 0.35 eV in the top region of the valence band, while the monolayer film has a single peak in the same energy region of the photoelectron spectra. The energy difference of 0.34 eV between the HOMO and the HOMO − 1 of the Type 3 dimer is fairly consistent to their conclusions. The calculation results would give an interesting indication that the ionisation threshold of the PbPc solids will be affected strongly by their polymorphic systems. Even in the thin films of PbPc such as the bilayer one, a molecular arrangement is crucial to photoelectric properties of the film. The molecular arrangement of Type 1 will give the most small ionisation energy, while those of Types 2 and 3 will affect slightly on the first ionisation energy but strongly on the second ionisation one. In contrast with this, the photoelectric properties for the thin films of SnPc would be insensitive to the molecular arrangement.
Ionisation potentiala | Electron affinitya | ||
---|---|---|---|
a Vertical ionisation potential and electron affinity by the ΔSCF.b Ref. 20. | |||
Monomer | 6.19 (6.26)b | 2.02 (2.12) | |
Dimer | Type 1 | 5.66 (5.76) | 2.39 (2.42) |
Type 2 | 5.91 (5.87)b | 2.15 (2.29) | |
Type 3 | 5.82 (5.84)b | 2.36 (2.45) | |
Trimer | Type 4 | 5.37 (5.49) | — |
Type 5 | 5.69 (5.68)b | — |
State | Main configuration (|C| ≧ 0.30) | Ea | fb | pc | ||
---|---|---|---|---|---|---|
a Excitation energy in eV.b Oscillator strength.c Transition moment direction. | ||||||
Type 1 | C4v | 3E | ±0.47(44a1 → 67e)[45%] − 0.40(26a2 → 68e)[32%] | 1.90 | 0.219 | x + y |
4E | +0.46(44a1 → 67e)[42%] ± 0.32(26a2 → 68e)[20%] | 2.05 | 0.256 | x + y | ||
Type 2 | C2h | 2Bu | +0.54(61au → 62bg)[58%] + 0.43(61bg → 62au)[36%] | 1.93 | 0.568 | x + y |
2Au | +0.58(61au → 74ag)[66%] − 0.37(61bg → 74bu)[28%] | 1.94 | 0.492 | z | ||
3Au | +0.43(73ag → 62au)[37%] + 0.52(73bu → 62bg)[55%] | 2.94 | 0.049 | z | ||
3Bu | −0.46(73ag → 74bu)[43%] + 0.50(73bu → 74ag)[50%] | 2.95 | 0.057 | x + y | ||
Type 3 | C2h | 1Bu | +0.70(64bg → 65au)[98%] | 1.60 | 0.481 | x + y |
1Au | +0.67(64bg → 71bu)[90%] | 1.68 | 0.116 | z | ||
2Au | +0.66(64au → 71ag)[87%] | 2.06 | 0.441 | z | ||
2Bu | +0.69(64au → 65bg)[95%] | 2.11 | 0.287 | x + y | ||
3Au | +0.58(70ag → 65au)[66%] − 0.37(70bu → 65bg)[27%] | 2.83 | 0.050 | z | ||
3Bu | +0.43(70ag → 71bu)[37%] + 0.53(70bu → 71ag)[56%] | 2.83 | 0.094 | x + y |
The electronic ground states of the monomer, and the model dimers and trimers under investigation are all totally symmetric singlet states, which are labelled as 1A1, 1Ag and 1A under the C4v, C2h and C1 point groups, respectively; therefore, the optical allowed excited states belong to the 1A1 or 1E, 1Au or 1Bu, and 1A (all states) representations, respectively. The results for the monomer were similar to those of the previous work12 in spite of the different functional used. The calculated results showed clearly that the solid state effects would be observed obviously in the Q-band region as broadening with some features, but not in the B-band region.
In the case of the monoclinic models (Types 1 and 4), there would exist low energy excitations in the region below 1.7 eV with small oscillator strengths. These states consisted of one-electron transitions between the molecular orbitals originated from the HOMO and LUMO of the monomer, which could be assigned as the π–π* transitions. Contrary to these small oscillator strengths of these states, the excited states around 2 eV had relatively strong oscillator strengths. These states contained one-electron transitions from the molecular orbitals with large contribution from the s-type atomic orbital of the central metal atom in addition to the π–π* transitions; therefore these excited states would be assigned to the mixed states with the π–π* and the metal to ligand charge-transfer (MLCT) transitions. A series of the excited states would contribute together to a strong absorption profile over 3.5 eV, which generates the B-band. It would be concluded from these results that the Q-band observed for the monoclinic system has a strong contribution of the MLCT. This result differs remarkably from those discussed for the free molecules.12,18
The triclinic models (Types 2, 3 and 5) showed a different excitation profile as compared to the monoclinic ones. The excited states below 2.1 eV consisted of one-electron transitions between the molecular orbitals originated from the HOMO and LUMO of the monomer. These states could be assigned to the π–π* transition states. On the contrary, the excited states with contributions of the MLCT transition showed small oscillator strengths, and their excitation energies were estimated to be about 2.9 eV. It would be concluded for the triclinic system that the Q-band has the π–π* transition states as a main component, and the MLCT band be observed at the region between the Q- and B-band with a small intensity as compared to the Q-band. These excitation features for the triclinic solid are similar to those of the free molecule except the band broadening due to interactions between molecules in solids.
Fig. 5 shows the observed absorption spectra of the thin PbPc film before and after annealing by broken and solid curves, respectively, in the bottom panel. The absorption spectrum for a chloroform solution is also shown in the figure by a green curve. As expected by the calculations, the absorption features of the thin film differ remarkably from those of a solution especially in the Q-band region. Before annealing, the film showed an absorption band at 1.7 eV with a shoulder at about 1.9 eV in the Q-band region, and the broad absorption band around 3.6 eV as the B-band. The weak feature around 2.8 eV was also observed in the region between the Q- and B-bands. After annealing, the thin film exhibited three absorption features at 1.40, 1.74 and 1.85 eV in the Q-band region, where the feature at 1.4 eV had the strongest intensity among them. The position of the B-band was independent of annealing. It should be noted that the feature at 2.8 eV was intensified by annealing as compared to the others. These experimental results were the same as those reported previously.36–38 Particularly, Miyamoto et al. have also reported that the monoclinic crystalline film of PbPc changes its crystalline form to the triclinic system by annealing through the detailed investigation.36 The present calculation results with the model dimers represented well all of the experimental results. The present work would also show that the electronic structure in a solid state can be investigated sufficiently with using the relatively small models chosen appropriately.
The stabilities of the model dimers and trimers were estimated from the differences in the Gibbs free energies between them and the corresponding components. The results revealed clearly that PbPc would show a possibility to take both molecular orientations of the monoclinic and triclinic crystal systems. Contrary to this, SnPc may prefer the triclinic molecular orientation strongly to the monoclinic one.
The excitation energies of the three dimers and two trimers were calculated by the TD-DFT method. The theoretical absorption spectra could explain well the experimental results observed for the thin film including its evolution along with polymorphic change. The calculated results would solid state effects were obviously observed especially in the Q-band region. It would be interestingly concluded from the calculation results that the π–π* and MLCT transitions contribute differently to the excited states of the monoclinic and triclinic solids of PbPc. Particularly, the main contribution to the Q-band could be assigned as the π–π* and MLCT transitions for the triclinic and monoclinic systems, respectively.
Finally, this work has been shown that the electronic structures in the solid state can also be investigated theoretically with using the small models properly selected.
Footnote |
† Electronic supplementary information (ESI) available: Cartesian coordinates, several molecular orbitals and all calculated excitation energies for all species, and the results of elemental analysis for sublimed PbPc samples and X-ray powder pattern recorded with the evaporated thin film of PbPc. See DOI: 10.1039/c6ra27269j |
This journal is © The Royal Society of Chemistry 2017 |