Computational prediction of optimal metal ions to induce coordinated polymerization of muscle-like [c2]daisy chains

Recently, a muscle-like organometallic polymer has been successfully synthesized using Fe 2+ as a linker atom. The polymer exhibits acid–base controllable muscle-like expansion and contraction on the micrometer scale. Further development could be facilitated by revealing the polymerization mechanism and by searching for optimal linker atoms. In this work, we have examined possible equilibrium and intermediate polymer structures, which consist of [c2]daisy chains linked by divalent transition metal ions (Sc 2+ , Ti 2+ , Fe 2+ , Co 2+ , Ni 2+ or Zn 2+ ) with various hexa-coordination arrangements, based on calculations using density functional theory. We find that the metal linkers in polymers are weaker in acid than in base due to excess positive charges on the polymer, leading to their thermodynamical instability or even decomposition. This can explain the experimental difficulty in improving the degree of polymerization for metal-linked polymers. We also find that the polymers with either Fe 2+ or Co 2+ are the most favorable, with the latter extending 1.4% longer than with the former. Since Fe 2+ has been confirmed experimentally to be a successful linker, Co 2+ would function equally well and thus could be used as an alternative choice for polymerization.


Introduction
For decades, scientists have synthesized and examined muscle-like functional materials endowed with mechanical expansion and contraction triggered by chemical, [1][2][3][4][5][6][7] electrochemical, 1,[8][9][10][11][12][13][14][15] or photochemical 10,[16][17][18][19][20][21] stimuli under human control.In 2012, Giuseppone et al. 22 successfully synthesized a polymer containing [c2]daisy chains and Fe 2+ ions, which is acid-base controllable and micrometer-scale responsive while having muscle-like capability.As revealed by measurements of light scattering and small-angle neutron scattering, 22 the unit lengths of the extended and contracted polymer units were found to be 5.4 and 3.2 nm, respectively.The 2.2 nm length difference in each unit is amplified in a linear polymer with about 3000 units to cause a stretch to 15 900 from 9400 nm, thus realizing muscle-like motion in a quasi-macroscopic regime.This achievement to some extent overcame the challenge of a low degree of polymerization and significantly advanced the technological prospect of practical muscle-like functional materials. 23 [c2]daisy chain is a kind of doubly-threaded rotaxane.Each [c2]daisy chain supramolecule consists of two identical rodand-ring molecules that are threaded to each other; they have the capability of telescoping as well as folding, which enables expansion and contraction (shown in Fig. 1) of each monomer.By adding the linking metal ion Fe 2+ , polymerization can be induced.As schematically shown in Fig. 1, the polymer unit is composed of one [c2]daisy chain and one Fe 2+ ion, which are linked by organometallic hexa-coordination.The number of units (n) in the experimentally synthesized polymer was found to be about 3000 by Giuseppone's group. 22Each half of the [c2]daisy chain is threaded through the other half's dibenzo [24]crown-8 (DB24C8) ether ring, while the bulky stopper (terpyridine group together with its adjacent ethinyl and phenyl species) prevents dethreading and also serves as a polymerizing linkage to the neighboring [c2]daisy chain by means of a bridging Fe 2+ .Each half of the [c2]daisy chain contains two recognition sites located in the thread, which serve as ''docking stations'' for the contracted vs. extended forms of the polymer: one is a triazolium ion (C 2 N 3 + ) and the other is a neutral amino group (NH) at high pH, which changes to a secondary dialkyl ammonium ion (NH 2  This journal is © the Owner Societies 2016 environment.The length of each unit varies dynamically as the binding affinity of DB24C8 for the two docking stations changes with pH: in acidic solution (e.g., trifluoroacetic acid), the NH 2 + preferentially attaches to the crown ether, resulting in the extended form; in basic solution (e.g., sodium hydroxide), the deprotonated NH has lower affinity for the crown ether, which therefore binds with the triazolium ion, favoring the contracted form.The contraction/extension cycle of the above described single [c2]daisy chain 24 supramolecule has been studied in our previous computational simulations 25 using the monomer of this muscle-like polymer. 22,23ur systematic quantum chemical and molecular dynamics calculations have revealed that the polymer length change is not only due to the linear telescoping between the two docking stations but also due to the folding/unfolding of the two bulky stoppers, 25 which is a striking difference as illustrated in Fig. 1 (top vs. bottom).
Following a bottom-up strategy, the experimental quasimacroscopic muscle-like behavior originates from the uniform change of many identical organometallic units.Thus, extending our prior study of the behavior of an individual [c2]daisy chain, investigations of the characteristics of the polymer containing different transition metal ions carried out in this work can clarify important aspects of this system and shed light on the mechanism of polymerization induced by the metal ions, which in turn can help in the design of optimal linkers.

Models and computational details
In this work, we choose Sc 2+ , Ti 2+ , Fe 2+ , Co 2+ , Ni 2+ or Zn 2+ to evaluate their effects on the polymer structures as they belong to the same period of elements with 3d 1 , 3d 2 , 3d 6 , 3d 7 , 3d 8 or 3d 10 valence electrons, respectively.It is known that the Lewis acidity 26,27 of Fe 2+ , Co 2+ , Ni 2+ and Zn 2+ ranges between soft and hard, giving some flexibility to coordinate with soft or hard Lewis bases.Nonetheless, their chemical characteristics may differentiate their role in inducing coordinated polymerization in our case and thus allow optimizing their use in applications.Although Sc 2+ and Ti 2+ are not the most frequent oxidation states of these elements, they are chosen to track the dependence of polymer energetics on the number of 3d electrons.Since Fe 2+ or Zn 2+ was used in experiment 22 to induce polymerization, we use them as the comparative reference in this work.
This work uses the same approach and terminology as our previous study of the monomeric [c2]daisy chain. 25 ), which are then left out of the calculation.The metal ions are, before further optimization, inserted at equal distances (B2 Å) from nearby nitrogen atoms of the terpyridines so as to hexa-coordinate with two [c2]daisy chains, which connect head-to-tail through a metal ion due to a periodic boundary condition.In the dimensions perpendicular to the chain, there are large intermolecular distances to avoid discernible interaction between parallel polymers.
In the experiment by Giuseppone et al. 22 the cationic [c2]daisy chains and Fe 2+ ions were polymerized in an ionic liquid, where stoichiometric PF 6 À and Cl À counterions accompany the polymers to neutralize the charge in the environment of organic solvents (acetonitrile CH 3 CN and chloroform CHCl 3 in a 1 : 1 mixture solution À and Cl À counterions in our calculations as they usually interact only indirectly with the polymer and float around the cationic polymer in complexes with solvent molecules.We have previously 25 presented the stable complexes PF 6 À (CH 3 CN) n or PF 6 À (CHCl 3 ) n (n = 1-8) and verified that these solvents provide a soft negative electrostatic potential environment for [c2]daisy chains, which thus avoids strong electrostatic interactions between the [c2]daisy chains and the surrounding counter-ions while keeping the [c2]daisy chain in a functioning state during the muscle-like motion.To simplify the calculation model, we thus feel justified to study the cationic polymer under vacuum.
Due to the fact that [c2]daisy chains are of large size and contain weak interactions such as hydrogen bonds and/or p-p stacking, the self-consistent charge density-functional tightbinding dispersion (SCC-DFTB-D) method [28][29][30][31][32][33][34][35] is first utilized to optimize the molecular and unit cell structure.Dispersion contributions are included via Lennard-Jones potentials, which take the parameters from the universal force field.The Slater-Koster (SK) potentials of metals with C, H, O, N are taken from the parameter sets trans3d 36,37 and znorg, [38][39][40] which have been well tested to predict comparably accurate geometries using DFT.Applying periodic boundary conditions, the lattice component force is calculated together with the interatomic force during each SCC cycle.If the structure reaches the force optimization criterion but the lattice does not, the unit cell is automatically modified to other cell lengths until the lattice force converges.A global search is carried out to find the most favorable equilibrated and intermediate conformations of the polymer units corresponding to lowest total energies.
To verify the reliability of DFTB in predicting geometric structures in our work, we firstly study small complexes with hexa-coordination bonds formed between two terpyridines and one metal ion by SCC-DFTB-D and oB97XD/6-31G* optimization, as listed in Table S1 (ESI †).The results confirm that SCC-DFTB-D can successfully predict similar hexa-coordination geometries as DFT and can deal with periodic and complicated polymer structures in this study.To obtain reliable molecular orbitals and energetic data, single-point calculations using the geometries derived at the SCC-DFTB-D level of theory are carried out using oB97XD/ 6-31G*, a long-range hybrid density functional theory with dispersion corrections, available in the Gaussian 09 package. 41ote that the oB97XD 42 functional yields satisfactory accuracy for thermochemistry, kinetics, and non-covalent interactions.

Polymer units
The combined optimization of geometry and cell parameters yields polymer structures containing transition metal ions (M 2+ : Sc 2+ , Ti 2+ , Fe 2+ , Co 2+ , Ni 2+ , or Zn 2+ ).Fig. 2 [proton-cont 4+ -Ti 2+ ] n are increased to 2.055 and 2.057 Å, indicating bond breaking.In [proton-cont 4+ -Zn 2+ ] n , two DB24C8 rings are opened considerably and the hydrocarbon rods are no longer kept linear, counteracting the lengthening step.Thus, it appears that the conformations under the acidic conditions become unstable and may even decompose.The N-C bond cleavage phenomenon is worth noticing by the experimentalists as it probably impedes the smooth stretching of the polymer in its practical application.
To quantify the hexa-coordination behavior, we check the distances between the M 2+ and the six N atoms to which the M 2+ ion is coordinated in the terpyridine groups of [c2]daisy chains, as listed in Table 1.In most conformations, two shorter and four longer M 2+ -N bond lengths reflect reasonable metal coordination as suggested by the d 2 -sp 3 orbital hybrid.However, one particularly long distance (42.5 Å) appears in the systems  This suggests that the conformations are becoming more distorted and weaker in acid than in base.Ni 2+ is not a good linker ion for this muscle-like polymer as it cannot hexa-coordinate well with [c2]daisy chains, as shown by one unusual Ni 2+ -N distance (B3 Å) in all its configurations.The more uneven hexa-coordination bonds via Ni 2+ will probably make the polymer more easily disordered than other ions so as to fail in undergoing multiple muscle-like cycles.Fe 2+ -N bond lengths are in the range 1.9-2.1 Å (a distance larger than 2.5 Å is obtained in the unstable intermediate and not considered here), basically consistent with the experimental value of 2.06 Å, as reported in a similar complex. 43This again reflects that SCC-DFTB-D can ensure comparable coordination bond lengths with the experiment.
To examine the impact of the metal linker on the polymer length, we define the unit length (L) as the distance between two successive M 2+ ions, as labeled in Fig. 2 and Fig. S1 (ESI †) and then compared in Fig. 3.The unit length data are all listed in Table S2 (ESI †).Firstly, it should be noted that the theoretical L value of [ext 4+ -Fe 2+ ] n is found to be 5.38 nm, in good agreement with the corresponding experimental value of 5.4 nm.Furthermore, the muscle-like motion is clearly exhibited by the unit length changes among these four conformations.Starting from the equilibrated state in acid, as the solution is changed from acidic to basic, the [ext 4+ -M 2+ ] n state (leftmost solid red bars in Fig. 3) converts to the [deproton-ext 2+ -M 2+ ] n state (next bars toward the right) with the lengths slightly decreasing by 0.3-0.4nm.This change is induced by the slight folding of the bulky stoppers as displayed in Fig. 2 and Fig. S1 (ESI †).Next, the contraction process from the [deproton-ext 2+ -M 2+ ] n intermediate to the [cont 2+ -M 2+ ] n equilibrium (next blue bars toward the right) causes the unit length to greatly decrease by 1.3-1.5 nm.As the solution is changed from basic to acidic, [cont 2+ -M 2+ ] n becomes [proton-cont 4+ -M 2+ ] n (rightmost bars) and the unfolding of the bulky stoppers is clearly manifested as the unit length dramatically increases by 0.9-1.0nm.The cycle is completed when the extension process from the [proton-cont 4+ -M 2+ ] n intermediate to the [ext 4+ -M 2+ ] n equilibrium (leftmost solid red bars) shows unit lengths increasing by 0.6-0.8nm.This indicates that the pathway of the contraction and stretching in the polymer is basically consistent with our analysis for the individual [c2]daisy chain supramolecule and confirms that the folding of the two bulky stoppers occurs in the last phase of contraction while unfolding takes place in the first phase of stretching. 25The response of the polymer to the acid is obvious at first, namely, the initial unfolding process significantly increasing the polymer length.Afterwards, the linear extension in the following process is hampered by a higher barrier and becomes difficult as revealed in our previous study. 25By contrast, the structural response of the polymer to the base is not as apparent as to the acid at first, while the subsequent contraction is smoother than stretching with a lower energy barrier.Notably, the unit length difference between [ext 4+ -M 2+ ] n and [cont 2+ -M 2+ ] n equilibria plays a key role in determining the relative polymer contraction under acid-base actuation.As shown in Fig. 3, the highest relative contraction, 31.9%, is found in the case of Co 2+ , while the second and third highest, are 31.6%and 31.2% for Ni 2+ and Fe 2+ .Ni 2+ is excluded here due to unsuccessful [proton-cont 4+ -Ni 2+ ] n configuration.The main advantage of Co 2+ is that the [ext 4+ -Co 2+ ] n is perfectly extended 1.4% longer than [ext 4+ -Fe 2+ ] n , while [cont 2+ -Co 2+ ] n remains relatively as short as [cont 2+ -Fe 2+ ] n .Considering that Fe 2+ has been successfully used in experiment, we thus conjecture that the polymer containing Co 2+ might function equally well as that containing Fe 2+ .

Hexa-coordination of M 2+ with [c2]daisy chains
Integrated muscle-like polymerization was realized in experiment by using metal ions to link daisy chains together.A single M 2+ connects two [c2]daisy chains by forming the bis(terpyridine)-M 2+ coordination complex with coordination of the metal ion to 6 N atoms in 6 pyridine rings.Increasing the degree of polymerization to lengthen the polymer is a key for enabling macroscopic expansion and contraction for large-scale mechanical applications.Hence, it is important to study the coordination of M 2+ with rotaxanes at the molecular level.
It should be noted here that, compared to the [c2]daisy chain monomers, the polymerization reaction reorients 4 of the 6 pyridine rings in the two adjacent [c2]daisy chains such that 6 N atoms can bond to the linking metal ion: this converts the two adjacent terpyridines from trans to cis configuration.The trans form is a natural and stable state for [c2]daisy chains before polymerization and needs to be activated to the cis form and then to coordinate with M 2+ for creating a polymer.This trans to cis conversion deforms the terpyridines moderately.For example, the N-C-C-N dihedral angle highlighted by enlarged balls in the trans and cis forms is shown in Fig. S2 (ESI †) for the four conformers of the monomer.While terpyridine is close to planar in the monomeric trans form, it is tilted by about 201 between two pyridine parts in the monomeric cis form, except for cont 2+ which is planar, as shown in Fig. S3 (ESI †).The inevitable deformation will, to some extent, impact the coordination with metal ions in a symmetric coordination field.Nonetheless, the metal ions may try to flatten the terpyridine as they react with the [c2]daisy chains.The energy differences in the trans to cis conversions for the four states range from 151 to 156 kcal mol À1 at the oB97XD/6-31G* level of theory, cf.Fig. S3 (ESI †).This means that cis-terpyridine will be automatically converted to the trans form once the M 2+ -N coordination bond is broken and that more input energy is required to recover the polymerization reaction.
The formation of hexa-coordination bonds in an octahedral environment would require the use of the unoccupied 3d, 4s, 4p orbitals of the metal ions to interact with orbitals occupied by twelve electrons of the [c2]daisy chain.This is facilitated by two terpyridines each consisting of three nitrogen atoms each of which has a s donor electron pair.Fig. S4 (ESI †) shows the twelve electrons localized at two terpyridine groups, confirming that the polymerization reaction occurs at the two ends.The six 3d electrons of Fe 2+ fill the 3d xy , 3d xz and 3d yz orbitals, leaving behind six empty 3d z 2, 3d x 2 Àdy 2, 4s, 4p x , 4p y , 4p z orbitals pointing at the N ligands.Six fully filled bonding orbitals and six empty anti-bonding orbitals will be formed in the Fe-based hexa-coordination compound.For Co 2+ , Ni 2+ and Zn 2+ , one, two, or four more 3d electrons should fill the ''anti-bonding'' (e g *) orbitals, destabilizing more the hexa-coordination compound than Fe 2+ .For Sc 2+ and Ti 2+ , the only one or two 3d electrons should also fill the ''non-bonding'' (t 2g ) orbitals but their relatively large ion radii may be not as favorable as Fe 2+ for forming stable hexa-coordination bonds.Thus, Fe 2+ is the best choice.On the other hand, the Sc 2+ or Co 2+ has an unpaired electron as it hexa-coordinates with [c2]daisy chains.Thus, the spin signal from the Sc-based and Co-based polymers may also help in probing some changes from the polar or non-polar environment provided by the ionic liquid.
We have discussed above how the hexa-coordination is impacted by the conformational change.To examine the linking strength, i.e., the thermal stability of the organometallic coordinated polymers, we have calculated the hexa-coordination reaction heats (DE hexa ) between two [c2]daisy chains and M 2+ according to the equation: Here, all terpyridine groups in all hexa-coordination compounds are in cis form, and [c2]daisy chains can be cis-ext 4+ , cis-deproton-ext 2+ , ciscont 2+ , and cis-proton-cont 4+ .The structures containing Co 2+ are shown in Fig. S5 (ESI †) as an example, where the Co 2+ and six coordinated N are highlighted by enlarged balls.The polymerization reaction heats (DE hexa ) for M 2+ linking two [c2]daisy chains in four configurations were calculated at the oB97XD/ 6-31G* level of theory based on the SCC-DFTB-D geometries, as shown in Fig. 4. The single-point energy corrections at the higher level of DFT based on the structures predicted by SCC-DFTB-D can guarantee a correct trend of the hexa-coordination thermal stability (see Fig. S6, ESI †).Among the metal ions that we studied, the linking strength is in the order: Fe 2+ 4 Co 2+ 4 Ni 2+ 4 Zn 2+ E Ti 2+ 4 Sc 2+ .The structure with Fe 2+ presents the highest stability whereas the one with Sc 2+ is the least stable.This leaves Co 2+ as the theoretically next favorable choice for being applied in this muscle-like polymer.Our results in fact support the experimental observation 23 that, when Zn 2+ instead of Fe 2+ was used as the transition metal, the partial depolymerization of its extended polymer occurred while the contracted polymer was successfully obtained.
Moreover, Fig. 4 indicates that the polymerization of [c2]daisy chains with M 2+ occurs in the order: deproton-ext 2+ 4 cont 2+ 4 ext 4+ 4 proton-cont 4+ .The contraction in base is energetically more favorable than the stretching in acid, with the reaction heats for the contracted case being about twice as large as for the extended case.A smaller reaction heat would easily lead to thermal decomposition of the extended polymer, especially for the protonated intermediates, and thus stops expansion during a This journal is © the Owner Societies 2016 cyclic motion.It is actually difficult to polymerize cationic monomers, as more positive charges on ext 4+ or proton-cont 4+ may hinder the polymerization reaction due to a larger electrostatic repulsion between the tetra-cationic [c2]daisy chains and the di-cationic M 2+ .We have verified this by calculating the reaction heat for two neutral cis-terpyridine molecules coordinating with one M 2+ , which represents the linking segment in all polymer structures.As shown in Fig. 4, their reaction heats are all increased by the same amount, much larger than for any polymerization.Clearly, in the experiment there are PF 6 À and Cl À counterions in solution to compensate charged monomers and to help alleviate the Coulomb repulsion from its positive charges and ensure that polymerization occurs.Otherwise the hexa-coordination of M 2+ with the [c2]daisy chains would be weakened considerably.However, from a different perspective, counterions are actually counter-productive because this polymer relies on two positive recognition sites to perform contraction and expansion.We believe that this is the main reason why in experiment it has been difficult to further improve the degree of polymerization.We note that Fig. 4 merely shows a general trend of the polymerization ability for a charged [c2]daisy chain monomer as the computational model does not include the accompanying counterions.Nevertheless, it tells us that further improvement of the thermal stability of the polymer by selecting an optimal metal ion such as Fe 2+ and Co 2+ or moving to neutral [c2]daisy chain supramolecules may help in increasing the degree of polymerization.
To further reflect the interaction strength between the metal ion and the [c2]daisy chain, we additionally calculate the Mayer bond order [44][45][46] (MBO) in each hexa-coordination compound, as shown in Fig. 5.It is surprisingly found that the MBOs from Co 2+ are strongest except for coordinating with cis-deprotonext 2+ , while those from Zn 2+ are weakest.The MBO order is: Co

Conclusions and perspectives
The are both optimal metal ions for polymerization by theoretical predictions.Furthermore, the extended polymer linked by Co 2+ is 1.4% longer than that by Fe 2+ , producing a 0.7% higher relative contraction than with Fe 2+ .As Fe 2+ has been shown in experiment to be a successful linker, Co 2+ ions would function  equally well and thus could be used as an alternative choice for polymerization.

Fig. 1
Fig. 1 Schematic representations of the polymerized [c2]daisy chain linked with Fe 2+ ions (yellow balls).Each half of the [c2]daisy chain consists of a hydrocarbon rod (straight segments) with two docking stations (blue and purple balls), a crown ether ring DB24C8 (red ring), and a terpyridine end group (green arc).Two identical halves are mutually threaded to form a [c2]daisy chain unit, docked (top panel) at the triazolium + station (in high pH) in the contracted state, or (bottom panel) at the NH 2 + station (in low pH) in the extended state.
(a) and (b) and Fig. S1(a)-(d) (cf.ESI †) display two equilibrium structures and two intermediate polymer structures of three units of an infinite polymer chain for each of the six metal ions.Hereafter, the extended equilibrium and protonated intermediate polymers under acidic conditions are called [ext 4+ -M 2+ ] n and [proton-cont 4+ -M 2+ ] n , respectively, while the contracted equilibrium and deprotonated intermediates in base are denoted as [cont 2+ -M 2+ ] n and [deprotonext 2+ -M 2+ ] n , respectively.Here n stands for infinite repetition in the polymers.The muscle-like motion cycles among four conformations follow the generic sequence shown in Scheme 1.Interestingly, our calculations show that not all of the metal ions can induce stable muscle-like polymers.While polymerization is accomplished when using Sc 2+ , Fe 2+ or Co 2+ , the use of Ti 2+ , Ni 2+ , or Zn 2+ induces instabilities in the protonated intermediate polymer, as shown in Fig. S1(b)-(d) (ESI †).In [proton-cont 4+ -Ti 2+ ] n and [proton-cont 4+ -Ni 2+ ] n , the N-C bonds near NH 2 + break up, fragmenting the [c2]daisy chains and interrupting further extension.For example, the distances between N (of NH 2 + ) and the neighboring C (of CH 2 ) in and [deproton-ext 2+ -Ni 2+ ] n .

Fig. 2 Scheme 1
Fig. 2 The equilibrium and intermediate polymer structures in three linked units containing (a) Fe 2+ or (b) Co 2+ , from SCC-DFTB-D calculations under vacuum.Gray, white, red, and blue atom colors represent C, H, O, and N, respectively.(Note that the terpyridine and benzene groups in the bulky stoppers are connected by -CC-groups, which in this representation look like abnormally long single bonds.)

2+ 4
Ti 2+ 4 Fe 2+ 4 Ni 2+ 4 Sc 2+ 4 Zn 2+ .It is seen that Co 2+ is indeed an option for being applied to the muscle-like [c2]daisy chain polymer.Lastly, an MBO analysis can help reflect the N-C bond weakening in the C-N(H/H 2 + )-C sites of all hexa-coordination compounds, as shown in Fig. 5.Under the basic conditions, the MBO values of the NH group with two neighboring alkyl C are 11-20% higher than those of NH 2 + under the acidic conditions, especially compared to the protonated intermediates.This indicates that the N-C bond strength is significantly reduced in the acidic environment, analogous to the hexa-coordination thermal stability.Since structural deformation in all hexa-coordination compounds is not as strong as in the realistic periodic polymers, we do not see the breakage of the N-C bond in compounds.However, it is expected from bond order results that the N-C breakage in the polymers will most likely occur in the most heavily distorted [proton-cont 4+ -M 2+ ] n , which is supported by the N-C cleavages in [proton-cont 4+ -Ti 2+ ] n and [proton-cont 4+ -Ni 2+ ] n .The N-C bond order change has no direct relationship with the metal type, see the four almost flat curves in Fig. 5, and will be significantly affected by the acid-base environment change.
study of equilibrium and intermediate polymer configurations in acid or base formed by linking muscle-like [c2]daisy chains with various transition metal ions reveals stable polymers linked by Sc 2+ , Fe 2+ and Co 2+ instead of Ti 2+ , Ni 2+ and Zn 2+ .The conformations in acid suffer strong distortions quite easily and can even decompose, especially for the protonated contracted polymer.The bond between N (of NH 2 + ) and its neighbor C (of CH 2 ) is much weakened due to acid catalysis.The positive excess charges on the [c2]daisy chains greatly reduce the reaction heat for polymerization and weaken the hexa-coordination links with metal ions, showing the difficulty to improve the degree of polymerization to macroscale as experienced in experiment.In hexa-coordination compounds, the thermal stability vary with the order: Fe 2+ 4 Co 2+ 4 Ni 2+ 4 Zn 2+ E Ti 2+ 4 Sc 2+ , while the Mayer bond order has the order: Co 2+ 4 Ti 2+ 4 Fe 2+ 4 Ni 2+ 4 Sc 2+ 4 Zn 2+ .Thus, Fe 2+ and Co 2+

Fig. 5
Fig. 5 Calculated Mayer bond orders of the hexa-coordination bonds (in columns) and the C-N bonds at the C-N(H/H 2 + )-C sites (in lines) at the level of oB97XD/6-31G*.
).The neutral polymer unit in acid or base is actually composed of C 116 H 134 O 16 N 14 Fe(PF 6 ) 4 Cl 2 or C 116 H 132 O 16 N 14 Fe(PF 6 ) 2 Cl 2 .However, we do not include the PF 6

Table 1 M
2+-N distances in polymers calculated by SCC-DFTB-D; the six values in each situation correspond to the hexa-coordination of the ion to six N atoms; italic values are abnormally long, indicating instability of those conformations (unit: Å) This journal is © the Owner Societies 2016 Phys.Chem.Chem.Phys., 2016, 18, 7419--7426 | 7423