{'fraction_non_alphanumeric': 0.03808362369337979, 'fraction_numerical': 0.032926829268292684, 'mean_word_length': 4.5130618517095655, 'pattern_counts': {'":': 0, '<': 3, '<?xml version=': 0, '>': 0, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 1, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Theoretical investigation of interactions between palladium and fullerene in polymer

Applying density functional theory (DFT) calculations, we predict the structural and electronic properties of different types of palladium–fullerene polymers. We examine the structures of one(1-D), two(2-D), and three-dimensional (3-D) polymers. We find that the most stable polymer is that represented by bonding via the [6,6] position of the fullerenemolecules with Pd in a distorted tetrahedral coordination. Special attention is paid to the electronic structure. We demonstrate clearly that changes in the Pd coordination geometry strongly affect the projected density-of-states picture of the 4d orbitals. The energy band gaps in the 1-D and 2-D systems obviously differ from that in the 3-D one; thus, we can directly modify the electronic properties of polymers. The results at the AM1* level of theory for the reduced 1-D polymer show that isolated polarons are the preferred electronic states.


Introduction
Many studies of fullerene-containing polymers have been reported recently, because of their potential application in many areas such as solar energy conversion, bio-and chemisensors, catalysis, adsorption and separation, energy accumulation, electroactive batteries, hydrogen storage, and electronics. [1][2][3][4] Fullerene cages can bond in many different ways to form a polymeric network. 2 Like alkenes, fullerenes can form homopolymers through [2 + 2] cycloadditions. [5][6][7][8][9] They can also be incorporated into a polymeric backbone [10][11][12] or attached to the main chain as pendant substituents. [13][14][15][16] An attractive approach to forming polymeric structures in which fullerene moieties are incorporated into the main polymeric backbone is to bond fullerene cages with transition metal atoms or their complexes. [17][18][19][20][21][22] To date, coordination polymers with fullerene cages covalently bonded to palladium and platinum atoms have been synthesised. 17,18,20,22 Structures with fullerenes coordinated through -Ir(CO) 2and -Rh(CF 3 CO 2 ) 2 -Rh-units have also been reported. 20, 21 Saito and co-workers recently reported the formation of one-dimensional (1-D) fullerene polymeric chains in which C 60 moieties are bonded to Ni[P(CH 3 ) 3 ] 2 units in h 2 -fashion. 19 Coordination polymers of fullerene and palladium, C 60 Pd x (x ¼ 1-3), have been investigated most intensively. They can be formed under chemical 17,18 and electrochemical conditions. 20,22 The chemical synthesis is carried out in benzene solution containing fullerene and a zerovalent palladium complex as precursors of the polymerisation process. Electrochemical polymerisation occurs during the reduction of Pd(II) ions in solution, in the presence of C 60 . The composition and properties of C 60 Pd x polymers depend on the concentration of the precursors in solution. By changing the polymerisation conditions, it is possible to produce materials that range from long 1-D chains to a highly cross-linked threedimensional (3-D) networks.
Very few studies have investigated the structure of C 60 Pd x and other fullerene-based coordination polymers. Electron diffraction and high-resolution electron microscopy demonstrated that C 60 Pd 3 exhibits an ordered body-centred cubic (bcc) structure of fullerene units with a slight rhombohedral distortion. 23 The fullerene cage is octahedrally coordinated to six palladium atoms. Each palladium atom is bonded to two C 60 units. X-ray photoelectron spectra of poly-Pd 3 C 60 indicated that charge is transferred from the metal atoms to the C 60 moieties during bond formation. 24 The structure of the polymer is closely related to its electronic properties.
The electrical properties of the C 60 Pd 3 polymer are particularly important in terms of potential application of this material in charge storage devices 25-28 and solar energy conversion systems. 29 The C 60 Pd 3 polymer shows low electrical conductivity in its neutral state. Reduction of the polymer results in a signicant increased conductivity 30 and creates mixed valence states along the polymeric chain. Electron exchange between the redox centres is responsible for the increase in the conductivity of the polymeric material.
Experimental investigations of the electronic properties of conducting polymers can be signicantly supported by theoretical studies. For a number of oligomeric and polymeric materials, such theoretical studies have contributed greatly to the rationalisation of their properties. [31][32][33][34] They also enable prediction of the electronic behaviour of yet unknown polymers. 35 Theoretical prediction of the electronic structure of macromolecular systems is an extremely important step towards the rational design of high-performance materials.
Up to now, there have been only a few theoretical studies on the interaction of fullerene with palladium. [36][37][38] Gal'pern and Stankevich 36,37 considered the structure and electronic conguration of such complexes with different ratios of both components within the DFT approximation. They showed the formation of the quasi-one-dimensional polymer structure C 60 Pd 3 in which the adjacent fullerenes are bound via Pd 3 clusters. More recently, El Mahdy 38 studied hydrogen adsorption and dissociation on Pd-doped C 60 . Different types of fullerene-metal complexes with PH 3 ligands were analysed, showing that the bond dissociation energies increase in the order Pd < Pt < Ni. 39 Accurate theoretical calculations of metalcontaining fullerene complexes concerning the electronic structures and optical spectra properties have been reported by Zhou and Zhao. 40 In this paper, we analyse the geometrical and electronic structure of different types of C 60 Pd polymers using rstprinciples periodic density functional theory (DFT) with some contribution from the semiempirical AM1* method. It is well accepted that expanding the dimensionality of polymers may signicantly change their physical and chemical properties. In our studies, we consider 1-D, two-dimensional (2-D), and 3-D polymer structures. For the 1-D polymer, we analyse two types of Pd binding motif on the fullerene molecules and its coordination environment. Following these results, we study the 2-D and 3-D polymers. Here, we also report the electronic structure of all the polymers, which is of signicant importance for the fundamental understanding of their transport, optical, and electrochemical properties. Finally, we present the results of calculations obtained at the AM1* level of theory for the relative stability of a polaron pair vs. a bipolaron for the reduced 1-D polymer.

Computational methods
First-principles calculations were computed using the DFT periodic approach with a plane-wave basis set. The calculations were performed using the PWSCF 41 program to solve the Kohn-Sham equations and obtain the electronic ground state. The interaction between the ions and the valence electrons is described by Vanderbilt ultraso pseudopotentials. 42 For Pd and C, the 4d electrons on the one hand and the 2s and 2p electrons on the other hand were treated as valence electrons. The remaining electrons were kept frozen. Additional test calculations for a Pd pseudopotential that also included semicore states (4s, 4p) resulted in a 0.9 kcal mol À1 smaller binding energy than that including only 4d electrons; therefore, we used the latter type of pseudopotential. A plane-wave cutoff of 30 Ry was chosen for these calculations. Increasing the cutoff energy from 30 to 35 Ry caused an increase in the binding energy of less than 0.1 kcal mol À1 . The generalised gradient approximation (GGA)type Perdew-Burke-Ernzerhof (PBE) 43 exchange-correlation functional with Grimme's correction (D2) 44 was used to describe the van der Waals interactions throughout this work. This approach (DFT-D2) consists of adding a semiempirical dispersion potential (a pairwise force eld) to the Kohn-Sham DFT energy. The total energy was minimised using the Hellman-Feynman forces. The examined systems were fully relaxed until the residual forces on the ions were below 5 Â 10 À3 eVÅ À1 .
Periodic boundary conditions and a supercell approximation were applied for all the polymers studied. The 1-D polymer was described as an innite chain oriented along the z axis. Two supercells containing 61 atoms (one C 60 Pd unit) and 122 atoms (two C 60 Pd units) were used. The vacuum space in the lateral direction was set to 13Å. For the supercell with one C 60 Pd unit, the geometry was fully optimised using a 1 Â 1 Â 1, 1 Â 1 Â 2, 1 Â 1 Â 4, or 1 Â 1 Â 6 Monkhorst-Pack k-point mesh 45 for the integration in the rst Brillouin zone. The relative total energy differences between two subsequent k-point grids were À5.5, À0.1, and 0.0 kcal mol À1 , respectively. Therefore, a 1 Â 1 Â 4 grid for the supercell with one C 60 Pd unit and a 1 Â 1 Â 2 grid for the supercell with two C 60 Pd units were used. For the 2-D and 3-D polymers, we used tetragonal and bcc supercells, respectively; there were two types of supercell, small and large, with corresponding k-point meshes of 4 Â 4 Â 4 and 2 Â 2 Â 2.
To calculate the density of states (DOS), linear tetrahedron smearing was used. For better visualisation of the DOS, the peak width for broadening was set to 0.05 eV.
All the gures were produced using the visualisation programs VESTA 50 and XcrySDen. 51 3 Results and discussion 3.1 1-D polymers 3.1.1 Geometry. The fullerene molecule with icosahedral symmetry has two types of bonds, namely, 6-6 bonds (1.39Å) and 6-5 bonds (1.45Å). 52 Many different organometallic complexes of fullerene have been prepared to date, the vast majority of which show h 2 -coordination at the [6,6] position. 53 The stability of this structure is conrmed by theoretical studies of the interaction between C 60 and transition metal atoms (Co, Rh, Ir, Ni, Pd, and Pt). 54 The second most stable coordination type involves h 2 -coordination at the [6,5] position. 54 Therefore, we considered only these two energetically most favourable positions.
First, we examined the structure of the 1-D C 60 Pd polymer represented by two congurations with different Pd binding motifs on the fullerene molecules, namely, [5,6] and [6,6] binding, which are denoted as 56p and 66p, respectively. The initial (1 Â 1 Â 1) supercell contains one C 60 Pd unit in which the Pd atom is bound to two fullerene molecules by two h 2 -type Pd-C bonds. Next, we estimated the equilibrium Pd-C 60 -Pd distance (d). By varying the distance from 10.97 to 11.29Å in 0.02Å steps during the geometry optimisation, we obtained 17 structures for both polymers. The stability of these structures was compared using their total energy differences, where the most stable structure was used as a reference (E diff ) and set to 0 kcal mol À1 . The obtained equilibrium distances for the 56p and 66p polymers are 11.17 and 11.19Å, respectively (Fig. 1). The total energy difference between the most stable congurations of the 66p and 56p polymers is À6.7 kcal mol À1 , conrming that binding via the [6,6] position is energetically favourable. For both polymers, the Pd atoms exhibit a distorted square-planar geometry involving coordination to two C 60 molecules (four Pd-C bonds). It is noteworthy that varying the distance (d) from about 10.05 to 10.29Å changes the relative total energy by less than 1 kcal mol À1 , which shows that the Pd-C bonds are very exible. The equilibrium Pd-C bond lengths and the C-Pd-C angles are 2.184Å and 39.8 and 2.174Å and 38.8 for 56p and 66p, respectively. The corresponding 5-6 and 6-6 bonds adjacent to the Pd are elongated by 0.033 and 0.044 A, respectively, compared to those in the isolated molecule.
It is very likely that a more stable structure can be obtained by changing the orientation of the C 60 moieties with respect to the palladium. Therefore, in the second step, we analysed the gradual change in the Pd coordination geometry from distorted square-planar to distorted tetrahedral. This was done by sequential rotation of one C 60 molecule around the Pd/Pd axis in a (1 Â 1 Â 2) supercell containing two C 60 Pd units, in order to determine the energetic stability. The rotation angle (R) was dened as the dihedral angle between four carbon atoms directly bound to Pd. The procedure involved changing the value of the rotation angle in 15 steps in combination with geometrical optimisation. Here, we also included a mixed polymer type, denoted as 56-66p, with an alternating two [5,6] and two [6,6] palladium-carbon bond motif. Taking into account the symmetry of the individual polymer chains, rotation between 0 and 90 for 66p and 56-66p and between 0 and 180 for 56p were considered. Altogether, 27 nonequivalent congurations were created; their energies are summarised in Table 1. Application of the above procedure to the mixed 56-66p polymer led to nonlinear polymer with two different values of the rotation angle for each C 60 . Rotation of the C 60 molecule in the supercell representing the 66p polymer gives a much lower energy minimum than that of the 56p and mixed 56-66p polymers ( Table 1). The calculations yield an optimal dihedral angle (R opt ) of 88.5 and 17.0 /39.6 for the two latter polymers, respectively, where the total energies are only 3.4 and 2.7 kcal mol À1 lower than the initial ones, respectively. Furthermore, we note here that the relative total energies in the whole range of the rotation angle are very close to each other that should lead to an easy rotation of the fullerene molecules around the longitudinal axis. In contrast, the 66p polymer is characterised by the existence of a stable global minimum (R opt ¼ 88.2 ) that is 9.4 kcal mol À1 more favourable than that for R opt ¼ 0 . Based on Table 1, it can be seen that the system easily reaches the global minimum starting from different initial dihedral angles, while the R opt of 59.9 and 74.0 correspond to local minima (7.9 and 9.0 kcal mol À1 , respectively) which are energetically close to the global minimum. The most stable congurations of the C 60 Pd polymer, namely, 56p(R À 90), 66p(R À 90), and 56-66p(R À 15), are shown in Fig. 2 with the relevant Pd-C bond lengths. For clarity, we refer to them using the initial value of the rotation angle (R). The 66p(R À 90) polymer is 7.4 and 19.4 kcal mol À1 more stable than the 56-66p(R À 15) and 56p(R À 90) ones (Table 2), respectively. For both the linear 56p(R À 90) and 66p(R À 90) polymers, the preferred coordination of Pd is distorted tetrahedral, whereas for the nonlinear 56-66p(R À 15), it is distorted square-planar. The corresponding Pd-C distances (2.105-2.303Å) predicted by the DFT calculations for all the systems are presented in Fig. 2.
To date, there are no experimental data for the 1-D C 60 Pd polymer; however, slightly shorter bond lengths (2.087(8)-2.149(8)Å) were reported for a 1-D Ni-bridged fullerene polymer, 19 where Ni is coordinated to the 6-6 position of the fullerene molecule (h 2 -type bond) as well as to two trimethylphosphine ligands, which form a four-fold coordination environment around the metal centre. However, the Ni/Ni distance within the polymer chains was found to be slightly longer (11.230(6)Å at 280 K) than the Pd/Pd distances. Recently, 1-D C 60 Ru polymeric chains were theoretically studied by Leng et al. 55 using periodic DFT calculations. According to their calculations, the most stable polymer is that with Ru atoms in a h 2(6) -h 6  coordination geometry. However, they did not apply a systematic procedure to nd the most stable polymer structure.
3.1.2 Electronic structures. The band structure and DOS for three types of the 1-D polymers are shown in Fig. 3a-c. The le side of each graph shows the band structure, whereas the right side shows the DOS. The partial contributions from Pd states, obtained by projecting the self-consistent DOS onto the corresponding atomic orbitals (projected density of states, PDOS) are also presented. The Fermi energy is set to 0 eV. The calculated highest occupied molecular orbital (HOMO)-lowest unoccupied molecular orbital (LUMO) gaps for the 56p(R À 90), 56-66p(R À 15), and 66p(R À 90) polymers are 1.11, 1.13, and 1.39 eV, respectively. For comparison, the calculated HOMO-LUMO gap of the isolated C 60 molecule is 1.64 eV, in very good agreement with other GGA calculations. 56,57 Binding via the [6,6] position increases the band gap by 0.28 eV with respect to that for binding via the [5,6] position.
The inuence of the Pd coordination environment on the electronic structure is rather small and does not seem to play an important role. All the band structures show many dispersionless bands, which suggest strong localisation of these electronic states. The band dispersion appears mainly between À2.75 and À1.8 eV and originates from the overlap between the Pd 4d and C 2p/2s orbitals (Fig. 3). The conduction band (CB) consists of three separate groups of bands located in the 1-4 eV range. The electronic density isosurface [r(r) ¼ 0.001 eÅ À3 ] of the LUMO for the most stable 66p(R À 90) polymer is presented in Fig. 4. The charge density distribution is delocalized over the whole fullerene molecules, whereas the contribution from the palladium states is small. Consequently, conduction occurs predominantly by electron hopping through the fullerene networks. The band structure for the mixed 56-66p(R À 15) polymer, which contains Pd in a distorted square-planar coordination, differs somewhat from that of others containing Pd in a distorted tetrahedral coordination. Although the DOSs of all the 1-D polymers show some differences, the general appearance of the individual peaks is very similar. A comparison of the carbon PDOSs of isolated and in-polymer fullerene (Fig. 3b) shows that the interaction with Pd modies its HOMO and HOMOÀ1 levels. They become more dispersed owing to hybridisation with the Pd 4d states, which lie in the between À2.75 eV and the Fermi level. The remaining part of the valence band (VB) is rather weakly perturbed and slightly downshied from fullerene gas phase positions. A closer inspection reveals  that for the 66p(R À 90) polymer the top of the VB and the bottom of the CB are much narrower as compared to the others. This could be validated by spectroscopy (e.g. valence photoemission spectroscopy). Upon interaction with fullerene, the 4d states of Pd split into individual components. As seen in Fig. 5, the splitting is very similar for the 56p(R À 90) and 66p(R À 90) polymers with Pd in a distorted tetrahedral coordination. In both cases, the lowestlying d zx and d zy states (À2.75 to À1.8 eV, where the d zx states are covered by the d zy ones) make a similar contribution to the hybridisation with the C states. The other two low-lying 4d states, namely, d x 2 Ày 2 and d xy , are located predominantly in the À1.0 to À0.5 eV range, whereas d z 2 is more densely populated between À0.5 eV and the Fermi level. For the 56-66p(R À 15) polymer with Pd in a distorted square-planar coordination, the splitting of the 4d states is different. The highest intensity of the 4d states falls between À1.0 and À0.5 eV and consists of d zx and d xy for Pd1 and of d x 2 Ày 2 and d xy for Pd2 (Fig. 5c). The lowestenergy states are found to be d zy , whereas the d z 2 orbitals are localised close to or at the top of the VBs. Although the changes in the Pd coordination obviously have a strong effect on the PDOS picture of the 4d orbitals, both the CB and the VB of the 1-D polymer are less sensitive towards it.
The formation of Pd-C bonds is the results of charge transfer from Pd atoms to the interstitial C atoms of C 60 . For all three polymers, the obtained Löwdin charges 58 are almost the same, with an increase of +0.05e in the charge on the C atoms bound to the Pd atoms with respect to a fullerene gas phase. The variation in the Löwdin charges on all the other carbon atoms was negligible (<0.01e). We also estimated a charge transfer of  +0.21e to each C 60 molecule as the difference in the charge of carbon atoms in the gas phase relative to that in the polymer. This result shows that the charge transfer basically involves only two C atoms directly bound to each Pd atom. Taking into account the limitations of population analysis, this nding may be interpreted as a slight polarization of the Pd atoms, meaning that they provide a small amount of electron density to the fullerenes but remain in the neutral oxidation state.

2-D and 3-D polymers
3.2.1 Geometry. As we found that linking by the [6,6] position in the 1-D C 60 Pd polymer is energetically the most favourable, we then investigated the 2-D and 3-D structures. There are no experimental data for the 2-D palladium-fullerene polymer; however, the 3-D C 60 Pd 3 polymeric structure was reportedly produced under a large excess of the palladium complex. 17 Experimental techniques showed that the 3-D polymer exhibits an ordered bcc structure of fullerene units with a slight rhombohedral distortion. 23 Each palladium atom is bonded to two C 60 units, and each fullerene cage is octahedrally coordinated by six palladium atoms. Extended X-ray absorption ne structure studies of the C 60 Pd n (n ¼ 1.0-3.2) polymer series indicated that the Pd atoms are located at the centre of the C]C bond between two hexagons of the fullerene cage. 59 Two structures were chosen to represent the 2-D polymer; the rst contains one C 60 Pd 2 unit formula (62 atoms), and the second contains four C 60 Pd 2 unit formulas (248 atoms). The second structure, which has a distorted tetrahedral coordination of Pd (Fig. 2d), is energetically more stable by 4.5 kcal mol À1 (per C 60 Pd 2 unit) than the rst one, which has a distorted square-planar coordination. Both supercells are tetragonal with an optimised distance (d) of 11.13Å (22.26Å for the large supercell). The introduction of the second dimension clearly reduces the Pd-C 60 -Pd distance relative to that of the most stable 1-D polymer (11.19Å).
To create the 3-D polymer, we took two bcc supercells including 2 C 60 Pd 3 units (126 atoms) and 16 C 60 Pd 3 units (1008 atoms). The second structure, with a distorted tetrahedral coordination of Pd (Fig. 2e), is 5.6 kcal mol À1 (per C 60 Pd 3 unit) more stable than the rst one. We obtained a lattice constant of 11.21Å in both cases (22.42Å for the large supercell), which is slightly smaller (0.8%) than the experimental value of 11.3Å. 23 The carbon atoms in neighbouring C 60 molecules interact via van der Waals forces with calculated minimum distances of about 3.24Å. For comparison, the shortest experimental C-C distance between two neighbouring C 60 molecules in the facecentred cubic (fcc) C 60 crystalline structure is about 3.4Å. 60 This shows no signicant shrinkage of the fullerene sublattice bounded by Pd atoms compared to the fcc fullerene crystal. Note that different van der Waals-corrected DFT methods give slightly different van der Waals distances. 61 3.2.2 Electronic structures. The calculated energy bands of the 2-D and 3-D polymers are shown in Fig. 3d and e. The DOSs for both polymers are not presented, as a picture similar to that for the 1-D 66p(R À 90) polymer emerges. For the 2-D polymer, the band structure (G-Z) is very close to that of the 1-D one, with an energy band gap of 1.33 eV (compared to 1.39 eV for the 1-D polymer). Dispersed bands, appearing mainly between À2.75 and À1.75 eV, originate from the overlap between Pd 4d (d zx and d zy ) and C 2p/2s, and the CB also consists of three separate groups of bands.
For the 3-D polymer, we observe strong reduction of the electronic band gap relative to that of the 2-D structure (from 1.33 to 0.61 eV). The band structure of the 3-D system, which exhibits a larger band dispersion in this case, is also different from that of the 2-D system. Because the 3-D polymer structure is represented by bcc structure, and hence is twice as dense as the simple cubic structure, which has a relatively loose network, one may expect different band structures. The results clearly indicate that the electronic properties change with the dimensionality of the polymer, and this nding should be considered during the construction of electronic devices.

Polarons and bipolarons
The n-doped C 60 Pd polymer exhibits electrochemical activity at negative potentials, which is related to reduction of the fullerene molecules. 30 This process, which is responsible for the observed conductivity, produces bipolarons that move along the polymer chain in an electrolyte solution. 30 For p-doped oligomers, theoretical DFT studies have shown that such bipolaronic states (dications) are usually delocalised over several monomer units (e.g., for polythiophene 62 and polypyrrole 63 ). Theoretical studies unambiguously showed that the description of polaronbipolaron equilibrium depends on the polymer length. 63 Therefore, to comprehensively study the relative stability of the polaron pair and bipolaron (dianion) for the n-doped C 60 Pd polymer, we took a long chain of the 1-D 66p(R À 90) polymer including 16 C 60 Pd units, but at the expense of lower accuracy using the fast semiempirical AM1* technique. 46,47 The calculations were performed with the frozen geometry taken from the DFT(PBE + D2) calculations. For the singlet state, we used spinrestricted and spin-unrestricted wave functions, where the rst corresponds to a bipolaron state, and the second corresponds to a polaron pair (with some contribution from the bipolaron state). In addition, a triplet state was also calculated. The singlet case is signicantly more stable than the triplet for separate charges. The bipolaron state did not converge, but it appears to be even more unstable than the triplet. Fig. 6 shows the singly occupied molecular orbitals (SOMOs) for the spin-unrestricted (singlet and triplet) wave function. For the singlet state, the SOMOs have coefficients at the third and fourth opposite terminal fullerene molecules, whereas for the triplet, the coef-cients include the seven fullerene molecules on both ends of the polymer and are separated by two C 60 Pd units.
To conclude, it is clear that in the reduced 1-D polymer, isolated polarons are the most stable, and, according to the experimental results, we can only speculate that bipolarons form if the charge carrier density is high. 30 However, further studies using introduced counterions and more sophisticated methods are needed to verify the current results.

Conclusions
We analysed the structural and electronic properties of different types of palladium-fullerene polymers using a periodic DFT approach. We studied 1-D, 2-D, and 3-D polymer structures. For each system, the optimal lattice constant was found.
Detailed calculations were performed for the 1-D polymer with two types of Pd-C connections via the [5,6] and [6,6] positions, including also a mixed type with an alternating two [5,6] and two [6,6] palladium-carbon bond motif. We found that the most stable polymer type is represented by bonding via the [6,6] position, whereas the least stable is that with bonding via the [5,6] position. For both polymers, the most stable coordination of Pd is distorted tetrahedral, whereas for the mixed polymer, it is distorted square-planar. The tendency of the geometry to change from distorted square-planar to distorted tetrahedral is especially evident for the most stable 1-D polymer. For all three types of 1-D polymer, a comparison of the carbon PDOS of isolated and in-polymer fullerene shows that the interaction with Pd modies its HOMO and HOMOÀ1 states. The changes in the Pd coordination strongly affect the PDOS picture of the 4d orbitals, however, both the valence and the conduction band are less sensitive towards it.
The 2-D and 3-D structures were created from the most stable 1-D polymer, where the 3-D system is described by a bcc supercell that shows good agreement with the experimental lattice constant. We found that the band structure depends on the spatial dimensionality of the polymer structure. Those of the 1-D and 2-D systems are very similar, in contrast to that of the 3-D system, where a strong reduction in the band gap is observed. Therefore, we should expect the electronic properties of the 3-D polymer to differ from those of the 1-D and 2-D polymers. Going further in this area, we showed, using the semiempirical AM1* method, that isolated polarons are the preferred electronic states for the reduced polymer. We believe our combined periodic DFT and semiempirical study will further motivate future studies on different metal-fullerene polymers or similar systems.