THE JOURNAL OF CHEMICAL PHYSICS 138, 014306 (2013)
DFT study of vibronic properties of d8 (Ni-, Pd-, and Pt-) phthalocyanines D. Pouladsaz,1,a) M. Schreiber,2 and T. G. Gopakumar3 1
Max Planck Institute for the Physics of Complex Systems, D-01187 Dresden, Germany Institute of Physics, Chemnitz University of Technology, D-09126 Chemnitz, Germany 3 Institute of Experimental and Applied Physics, University of Kiel, D-24098 Kiel, Germany 2
(Received 6 July 2012; accepted 6 December 2012; published online 3 January 2013) By means of density functional theory, we have studied the electronic structure and vibronic properties of single neutral NiPc, PdPc, and PtPc molecules and their singly and doubly ionized cations and anions. In particular, the vibronic couplings and reorganization energies of all systems are compared. Partitioning of the reorganization energy, corresponding to the photoelectron spectra of the first and second ionizations of studied molecules, into normal mode contributions shows that the major contributions are due to several vibrational modes with a1g symmetry and energies lower than 1600 cm−1 . The results reveal that the reorganization energy due to the singly positive ionization in the studied molecules is up to about one order of magnitude less than other reorganization energies. This makes these metal phthalocyanines, from the perspective of intramolecular reorganization energies, attractive as electron donor for intramolecular electron transfer in electron acceptor-donor systems. © 2013 American Institute of Physics. [http://dx.doi.org/10.1063/1.4773014] I. INTRODUCTION
Metal phthalocyanines (Fig. 1), a class of weakly semiconducting organic dye molecules, have attracted considerable attention due to their electronic, optical, and structural properties. The wide range of their application includes chemical sensors,1–3 devices using photoelectric conversion,4 hetero-junction organic solar cells,5, 6 optoelectronic devices,7–9 and data storage.10 Especially, their application in photoconductive and photovoltaic response is promising.11–13 Particularly, the photoinduced electron transfer in a system composed of fullerenes as electron-acceptor molecules covalently linked to phthalocyanines as photoactive donors has been observed.14 Better understanding of the influence of molecular structures of the donors and acceptors on photoinduced electron transfer reaction rates helps to improve the capture and storage of solar energy. In π -conjugated systems, the electronic transport properties are manipulated by the strong coupling between geometric and electronic structure.15–17 The electronic transport is efficient when the electronic coupling between adjacent molecules is maximized and the reorganization energy is minimized.18–21 Small reorganization energies in electron transfer reactions result in acceleration of the charge separation process and deceleration of the charge recombination process. We consider only the intramolecular reorganization energy which originates from the change in optimized geometry of the molecule due to the gain or loss of electrons.22–24 Density functional theory (DFT) has already been employed to study the decomposition of the reorganization energy in several organic molecules, including indolocarbazoles and several triarylamine derivatives, among the normal modes, i.e., in the internal coordinate space, in order to identify the molecular segment which has most cona) Electronic mail:
[email protected].
0021-9606/2013/138(1)/014306/6/$30.00
tributions to the reorganization energy.25 In case of metalfree phthalocyanines, the charge transport in their liquid crystalline phase has been investigated calculating the reorganization energy.26 By means of DFT, we study the electronic structure and vibronic properties of single neutral NiPc, PdPc, and PtPc molecules and their singly and doubly ionized cations and anions. The results point to a collective interaction of several vibrational modes with the molecular electronic structure rather than one particular mode. The gas-phase experiments of the reactions that remove/ donate electrons from/to the molecules which are singly or even multiply charged, e.g., buckminsterfullerene,27 provide a perspective to the fundamental process of electron transfer. The double-electron transfer reaction has been previously observed between triply charged buckminsterfullerene, and the polycyclic aromatic hydrocarbon molecules anthracene, corannulene, benzo[rst]pentaphene, and pyrene.28, 29 Besides, the stability of doubly ionized cations of metal phthalocyanines have been already reported.30 Therefore, we have also considered the double ionization in our study. II. COMPUTATIONAL METHODS
At the DFT level of theory, we previously31 performed geometry optimization for D4h symmetric single neutral molecules of NiPc, PdPc, and PtPc, along with singly and doubly charged cations and anions of all molecules. We showed that ionization shrinks the HOMO-LUMO gap of the studied molecules. In the current study, those DFT calculations, which utilized the hybrid functional B3LYP32–34 and def-TZVP basis set,35, 36 are extended to the determination of vibrational frequencies and normal modes. The calculations were performed with TURBOMOLE 5.10.37, 38 Moreover, we performed similar calculations with def2-TZVP39 basis set and def2-ecp40 effective core potential which includes relativistic effects for heavy transition metals (Pd and Pt). These
138, 014306-1
© 2013 American Institute of Physics
Downloaded 06 Aug 2013 to 193.175.8.59. This article is copyrighted as indicated in the abstract. Reuse of AIP content is subject to the terms at: http://jcp.aip.org/about/rights_and_permissions
014306-2
Pouladsaz, Schreiber, and Gopakumar
J. Chem. Phys. 138, 014306 (2013) TABLE II. HOMO-LUMO gap energy (eV) of neutral and ionized NiPc, PdPc, and PtPc. Doubly negative ionization
FIG. 1. Structure of a metal phthalocyanine molecule. M is the metal atom (e.g., Ni, Pd, or Pt).
calculations of the electronic structure give very similar results. In Tables I and II, we compare the energies of the frontier molecular orbitals and the HOMO-LUMO gap energies for both basis sets. To calculate the reorganization energy, we consider two potential energy surfaces for electronic states 1 and 2 which correspond to the neutral state and the ionized state of the molecule. The intramolecular reorganization energy λ consists of two terms called relaxation energies21, 41, 42 λ = λ(1) + λ(2) ,
(1)
where λ(1) is the difference of energies of state 1 in geometries 1 and 2, in which the geometry i denotes the molecular geometry optimized in state i. Similarly, λ(2) is the difference of energies of state 2 in geometries 1 and 2 (see Fig. 2). Expanding the potential energies of the neutral and the ionized states in a power series of the normal coordinates indicates how each vibrational mode contributes to the reorganization energy.42 The Huang-Rhys factors are the contribution coefficients which equal the square of so-called vibronic coupling constants (g2 ). Similar to our previous study,43 the HuangRhys factors corresponding to reorganization energies λ(1)
Singly negative ionization
Neutral molecule
NiPc PdPc PtPc
1.056 1.065 1.072
def2-TZVP + def2-ecp 1.907 2.209 0.541 2.277 0.543 2.334
NiPc PdPc PtPc
1.056 1.062 1.068
1.736 0.540 0.542
def-TZVP 2.232 2.283 2.332
Singly positive ionization
Doubly positive ionization
1.104 1.101 1.105
2.170 1.917 1.682
1.107 1.100 1.102
2.027 1.893 1.689
and λ(2) were calculated by projecting the deformation pattern onto the vibrational modes of the relevant neutral and charged molecules, respectively. The normal modes of the neutral and ionized states are generally related by a transformation consisting of a translation and a rotation. The rotation matrix is the so-called Duschinsky matrix44 which describes the mixing of the normal modes. We calculated the Duschinsky matrices by the DUSHIN program.41 III. RESULTS AND DISCUSSIONS
The results of projecting the deformation pattern onto the vibrational modes based on DFT/B3LYP calculations with def2-TZVP and def2-ecp effective core potential for NiPc are shown in Figs. 3–6 (for the comparison of the results of both basis sets for all studied molecules, see the supplementary material47 ), where the relaxation energies λ(1) and λ(2) are partitioned into contributions from the vibrational modes. In these figures, the panels on the left show the Huang-Rhys factors related to λ(1) where the deformation pattern is projected onto the vibrational modes of the relevant neutral molecule, whereas the panels on the right show the Huang-Rhys factors related to λ(2) where the deformation pattern is projected onto the vibrational modes of the relevant ionized molecule. The total relaxation energies, obtained from the normal mode analysis, are shown in Table III. Our calculations show that all vibrational modes contributing to the reorganization energy have a1g symmetry. Moreover, the vibronic features of
TABLE I. Energy (eV) of calculated frontier molecular orbitals for NiPc, PdPc, and PtPc. HOMO-1/HOMO-2 (eg )
HOMO (a1u )
LUMO/LUMO+1 (eg )
NiPc PdPc PtPc
−6.772 −6.586 −6.360
def2-TZVP + def2-ecp −5.117 −5.200 −5.236
NiPc PdPc PtPc
−6.657 −6.626 −6.435
−5.176 −5.255 −5.292
HOMO-LUMO gap (eV) Theora
Exptlb
−2.908 −2.923 −2.902
2.209 2.277 2.334
2.647 2.7 2.777
−2.944 −2.972 −2.960
2.232 2.283 2.332
2.647 2.7 2.777
def-TZVP
a b
B3LYP. Reference 31.
Downloaded 06 Aug 2013 to 193.175.8.59. This article is copyrighted as indicated in the abstract. Reuse of AIP content is subject to the terms at: http://jcp.aip.org/about/rights_and_permissions
014306-3
Pouladsaz, Schreiber, and Gopakumar
J. Chem. Phys. 138, 014306 (2013)
FIG. 2. Schematic view of potential energy surfaces for two electronic states as a function of molecular geometry Q, showing the vertical transitions (dashed lines), geometry deformation Q, and relaxation energies λ(1) and λ(2) .
λ(1) and λ(2) , corresponding to the molecules NiPc, PdPc, and PtPc with the same ionization, are roughly similar regardless of the central metal atom. This can be understood by noting that the influence of the central metal atom on the vibrational modes, contributing to the reorganization energy, is negligible. It can also be noted that the DFT calculations of reorganization energies and vibronic couplings with both basis sets give very similar results. Except for singly negative ionizations, the reorganization energy λ decreases when the atomic number of the central metal atom increases. In positive ionizations, the difference λ of the relaxation energies is less than 0.3 meV, which shows that the curvature of the PES of the molecule has been changed only slightly due to the ionization. This value is larger for negative ionizations (Table IV). Moreover, Table IV shows that, in each of the studied molecules, the reorganization energy due to the geometry deformation in singly charged cations is up to about one order of magnitude less than the reorganization energy due to other degrees of ionization, which clarifies why metal phthalocya-
nines have the potential of playing a role as electron donor for intramolecular electron transfer in acceptor-donor systems (e.g., see Ref. 14). Besides, the partitioning of the relaxation energy into normal mode contributions shows that the major contributions are due to several vibrational modes with wave numbers lower than 1600 cm−1 , except two CH stretching modes in the range of 3100–3400 cm−1 which have negligible contribution to the relaxation energy. In positive ionizations, the main contribution to the relaxation energy comes from highenergy vibrations and the contribution from low-energy vibrations to λ is negligible. This high-energy contribution is in fact distributed over several vibrational modes with wave numbers in the range of 1000–1600 cm−1 . For NiPc−1 , there is a very elongated mode in the low-energy region which makes the contribution of other vibrational modes to the reorganization energy negligible. For other negative ionizations, although there are very elongated modes in the low-energy region (between 200 and 350 cm−1 ), the contribution of other vibrational modes above 350 cm−1 is still significant. The difference in the relaxation of negatively and positively ionized molecules can be discussed on the basis of the electronic structure of molecular orbitals through which the vibrational relaxations occur. Positively charged molecules are formed by emptying the HOMO. As the HOMO possesses π -electron character, a positively charged molecule is relaxed via the vibrational excitations involving π -electrons (1000–1600 cm−1 , related to stretching modes of C=C and C=N). As the amount of charge increases (doubly ionized states), more π -electrons are involved in the relaxation process. This is seen as the increased contribution of frequencies in the range of 1000–1600 cm−1 (related to stretching modes of C=C and C=N). When an electron is added to the molecule, its degenerate LUMO and LUMO+1, which possess a certain dxz /dyz (14%) character, are occupied. Thus, the metal N in-plane bond bending in the range of 200–350 cm−1 may be involved in the relaxation of negatively ionized states as seen in the Huang-Rhys factors. It has been shown before for purely aromatic molecules that high frequency stretching modes of C=C are primarily involved in the relaxation during hole/electron transport.45, 46 Adding a heteroatom like S into the carbon chain of conjugated systems seems to strongly influence the low frequency modes of C–C bonds.45, 46
0.05
vibronic coupling g2
vibronic coupling g2
0.05 0.04 0.03 0.02 0.01 0.00
0
1000
2000
3000
wavenumber / cm-1
0.04 0.03 0.02 0.01 0.00
0
1000
2000
3000
wavenumber / cm-1
FIG. 3. The Huang-Rhys factors of the vibrational modes of neutral NiPc (left) and NiPc+1 (right) contributing to the relaxation energies λ(1) and λ(2) , due to the geometry deformation in singly positive ionization.
Downloaded 06 Aug 2013 to 193.175.8.59. This article is copyrighted as indicated in the abstract. Reuse of AIP content is subject to the terms at: http://jcp.aip.org/about/rights_and_permissions
014306-4
Pouladsaz, Schreiber, and Gopakumar
J. Chem. Phys. 138, 014306 (2013)
TABLE III. The relaxation energies λ(1) and λ(2) in meV, due to the geometry deformation in singly and doubly charged cations and anions of the molecules NiPc, PdPc, and PtPc. Doubly negative ionization
Singly negative ionization
Singly positive ionization
λ(2)
Doubly positive ionization
λ(1)
λ(2)
λ(1)
λ(1)
λ(2)
λ(1)
λ(2)
NiPc PdPc PtPc
121 130 137
118 128 134
def2-TZVP + def2-ecp 207 193 19 29 29 18 31 31 18
19 18 18
66 63 63
66 63 63
NiPc PdPc PtPc
121 130 135
119 128 132
181 29 30
def-TZVP 172 29 30
19 18 18
66 64 63
66 63 63
5.0
vibronic coupling g2
vibronic coupling g2
5.0
19 18 18
4.0 3.0 2.0 1.0 0.0 0
1000
2000
3000
4.0 3.0 2.0 1.0 0.0 0
wavenumber / cm-1
1000
2000
3000
wavenumber / cm-1
FIG. 4. As Fig. 3, but for NiPc (left) and NiPc−1 (right), singly negative ionization. 0.20
vibronic coupling g2
vibronic coupling g2
0.20
0.15
0.10
0.05
0.00 0
1000
2000
0.15
0.10
0.05
0.00
3000
0
wavenumber / cm-1
1000
2000
3000
wavenumber / cm-1
FIG. 5. As Fig. 3, but for NiPc (left) and NiPc+2 (right), doubly positive ionization. 0.6
vibronic coupling g2
vibronic coupling g2
0.6
0.4
0.2
0.0 0
1000
2000
3000
wavenumber / cm-1
0.4
0.2
0.0 0
1000
2000
3000
wavenumber / cm-1
FIG. 6. As Fig. 3, but for NiPc (left) and NiPc−2 (right), doubly negative ionization.
Downloaded 06 Aug 2013 to 193.175.8.59. This article is copyrighted as indicated in the abstract. Reuse of AIP content is subject to the terms at: http://jcp.aip.org/about/rights_and_permissions
014306-5
Pouladsaz, Schreiber, and Gopakumar
J. Chem. Phys. 138, 014306 (2013)
TABLE IV. The reorganization energy λ and absolute difference of the relaxation energies λ in meV, due to the geometry deformation in singly and doubly charged cations and anions of the molecules NiPc, PdPc, and PtPc. Doubly negative ionization λ
λ
NiPc PdPc PtPc
239 258 271
3 3 3
NiPc PdPc PtPc
240 258 267
3 3 2
Singly negative ionization
Singly positive ionization
Doubly positive ionization
λ
λ
λ
λ
λ
39 37 36
0 0 0
133 127 126
0 0 0
38 36 36
0 0 0
132 127 126
0 0 0
λ
def2-TZVP + def2-ecp 400 13 58 0 62 0 def-TZVP 353 9 58 0 60 0
The efficiency of these results can be checked by measuring the vibrational structure of the first and second ionization peaks of the ultraviolet photoelectron spectroscopy spectra for NiPc, PdPc, and PtPc. In order to improve our results, we computed the Duschinsky rotation matrices of normal modes. The results point to the absence of mixing of ag modes among the vibrations for cations, anions, and dianions of the studied molecules. There is only mixing of one pair of ag modes in NiPc+2 (around 1140–1210 cm−1 ), whereas two pairs of ag modes (around 1140–1210 cm−1 and 1410–1470 cm−1 ) of PdPc+2 and PtPc+2 are mixed. In all cases, the mixing of the modes is less than 20% (detailed in the supplementary material47 ). Therefore, neglecting the mixing of a few vibrational modes is still reasonable.
IV. CONCLUSION
We have presented a DFT study of the reorganization energy for the first and second ionizations of NiPc, PdPc, and PtPc molecules. Since, in each of these molecules, the reorganization energy due to singly positive ionization is less than other counterparts, from the perspective of intramolecular reorganization energies, the results indicate that d8 metal phthalocyanines, can better play the role of an electron donor, rather than an electron acceptor, in intramolecular electron acceptor-donor systems, although for a complete conclusion, the intermolecular reorganization energies, which are beyond the aim of our study, should also be taken into account. The partitioning of the reorganization energy into normal mode contributions shows that the major contributions are due to several vibrational modes lower than 1600 cm−1 . The only exception is NiPc−1 , for which there is a very elongated mode in the low-energy region which makes the contribution of other vibrational modes in the reorganization energy negligible. Moreover, there is only mixing of a few ag modes in the region of 1000–1600 cm−1 among the vibrations contributing to the reorganization energy for anions and dications of the studied molecules, excluding NiPc−1 , due to the Duschinsky rotation. The DFT calculations with effective core potential which considers relativistic effects give very similar results to those without effective core potential.
ACKNOWLEDGMENTS
The first author wishes to thank A. Abbasi, L. Gisslén, R. Scholz, A. Langari, F. Petraki, and T. Chassé for helpful discussions, and S. Wüster for reading the first draft of this article and giving valuable comments. 1 H.
Mockert, D. Schmeißer, and W. Göpel, Sens. Actuators, B 19, 159 (1989). 2 R. Zhou, F. Josse, W. Göpel, Z. Z. Öztürk, and Ö. Bekaroglu, Appl. Organomet. Chem. 10, 557 (1996). 3 S. A. Wring, J. P. Hart, L. Bracey, and B. Birch, Anal. Chim. Acta 231, 203 (1990). 4 M. A. Loi, P. Denk, H. Hoppe, H. Neugebauer, C. Winder, D. Meissner, C. Brabec, N. S. Sariciftci, A. Gouloumis, P. Vázquez, and T. Torres, J. Mater. Chem. 13, 700 (2000). 5 K. Walzer, B. Maennig, M. Pfeiffer, and K. Leo, Chem. Rev. 107, 1233 (2007). 6 S. Siebentritt, S. Günster, and D. Meissner, Synth. Met. 41, 1173 (1991). 7 S. Wang, Y. Liu, X. Huang, G. Yu, and D. Zhu, J. Phys. Chem. B 107, 12639 (2003). 8 K. Xiao, Y. Liu, X. Huang, Y. Xu, G. Yu, and D. Zhu, J. Phys. Chem. B 107, 9226 (2003). 9 Q. Yong, G. Yudi, W. Peng, and W. Liduo, Appl. Phys. Lett. 80, 2628 (2002). 10 L. Ma, Q. Xu, and Y. Yanga, Appl. Phys. Lett. 84, 4908 (2004). 11 J. B. Whitlock, P. Panayotatos, G. D. Sharma, M. D. Cox, R. R. Sauers, and G. R. Bird, Opt. Eng. 32, 1921 (1993). 12 G. McHale, M. I. Newton, P. D. Hooper, and M. R. Willis, Opt. Eng. 6, 89 (1996). 13 K.-Y. Law, Chem. Rev. 93, 449 (1993). 14 H. Neugebauer, M. A. Loi, C. Winder, N. S. Sariciftci, G. Cerullo, A. Gouloumis, P. Vázquez, and T. Torres, Sol. Energy Mater. Sol. Cells 83, 201 (2004). 15 W. P. Su, J. R. Schrieffer, and A. J. Heeger, Phys. Rev. Lett. 42, 1698 (1979). 16 J. L. Bredas and G. B. Street, Acc. Chem. Res. 18, 309 (1985). 17 W. R. Salaneck, R. H. Friend, and J. L. Bredas, Phys. Rep. 319, 231 (1999). 18 R. A. Marcus, Discuss. Faraday Soc. 29, 21 (1960). 19 R. A. Marcus and N. Sutin, Comments Inorg. Chem. 5, 119 (1986). 20 N. S. Hush, Coord. Chem. Rev. 64, 135 (1985). 21 P. F. Barbara, T. J. Meyer, and M. A. Ratner, J. Phys. Chem. 100, 13148 (1996). 22 A. Devos and M. Lannoo, Phys. Rev. B 58, 8236 (1998). 23 T. Kato and T. J. Yamabe, J. Chem. Phys. 115, 8592 (2001). 24 V. Coropceanu, J. M. Andre, M. Malagoli, and J. L. Bredas, Theor. Chem. Acc. 110, 59 (2003). 25 H. Geng, Y. Niu, Q. Peng, Zh. Shuai, V. Coropceanu, and J.-L. Brédas, J. Chem. Phys. 135, 104703 (2011). 26 J. Tant, Y. H. Geerts, M. Lehmann, V. De Cupere, G. Zucchi, B. W. Laursen, Th. Bjørnholm, V. Lemaur, V. Marcq, A. Burquel, E. Hennebicq,
Downloaded 06 Aug 2013 to 193.175.8.59. This article is copyrighted as indicated in the abstract. Reuse of AIP content is subject to the terms at: http://jcp.aip.org/about/rights_and_permissions
014306-6
Pouladsaz, Schreiber, and Gopakumar
F. Gardebien, P. Viville, D. Beljonne, R. Lazzaroni, and J. Cornil, J. Phys. Chem. B 109, 20315 (2005). 27 D. K. Bohme, Int. Rev. Phys. Chem. 13, 163 (1994). 28 Gh. Jahavery, H. Becker, S. Petrie, P.-Ch. Cheng, H. Schwarz, L. T. Scott, and D. K. Bohme, Org. Mass Spectrosc. 28, 1005 (1993). 29 D. Schröder, J. Loos, H. Schwarz, R. Thissen, D. V. Preda, L. T. Scott, D. Caraiman, M. V. Frach, and D. K. Böhme, Helv. Chim. Acta 84, 1625 (2001). 30 S. M. Schildcrout, J. Am. Chem. Soc. 105, 3852 (1983). 31 T. G. Gopakumar, J. Meiss, D. Pouladsaz, and M. Hietschold, J. Phys. Chem. C 112, 2529 (2008). 32 A. D. Becke, Phys. Rev. A 38, 3098 (1988). 33 A. D. Becke, J. Chem. Phys. 98, 5648 (1993). 34 C. Lee, W. Yang, and R. G. Parr, Phys. Rev. B 37, 785 (1988). 35 A. Schaefer, C. Huber, and R. Ahlrichs, J. Chem. Phys. 100, 5829 (1994). 36 K. Eichkorn, F. Weigend, O. Treutler, and R. Ahlrichs, Theor. Chem. Acc. 97, 119 (1997).
J. Chem. Phys. 138, 014306 (2013) 37 R.
Ahlrichs, M. Bär, H.-P. Baron et al., TURBOMOLE 5.10, Universität Karlsruhe, Germany, 2008, see http://www.turbomole.com. 38 R. Ahlrichs, M. Bär, M. Häser, H. Horn, and C. Kölmel, Chem. Phys. Lett. 162, 165 (1989). 39 F. Weigend and R. Ahlrichs, Phys. Chem. Chem. Phys. 7, 3297 (2005). 40 D. Andrae, U. Häussermann, M. Dolg, H. Stoll, and H. Preuss, Theor. Chim. Acta 77, 123 (1990). 41 J. R. Reimers, J. Chem. Phys. 115, 9103 (2001). 42 M. Malagoli, V. Coropceanu, D. A. da Silva Filho, and J. L. Brédas, J. Chem. Phys. 120, 7490 (2004). 43 A. Staicu, G. Rouillé, T. Henning, F. Huisken, D. Pouladsaz, and R. Scholz, J. Chem. Phys. 129, 074302 (2008). 44 F. Duschinsky, Acta Physicochim. URSS 7, 551 (1937). 45 J. L. Bredas, D. Beljonne, V. Coropceanu, and J. Cornil, Chem. Rev. 104, 4971 (2004). 46 V. Coropceanu, J. Cornil, D. A. da Silva Filho, Y. Olivier, R. Silbey, and J. L. Bredas, Chem. Rev. 107, 926 (2007). 47 See the supplementary material at http://dx.doi.org/10.1063/1.4773014 for the comparison of the results of both basis sets for all studied molecules.
Downloaded 06 Aug 2013 to 193.175.8.59. This article is copyrighted as indicated in the abstract. Reuse of AIP content is subject to the terms at: http://jcp.aip.org/about/rights_and_permissions