Isaac
Alcón
^{a},
Daniel
Reta
^{a},
Iberio de P. R.
Moreira
^{a} and
Stefan T.
Bromley
*^{ab}
^{a}Institut de Química Teòrica i Computacional de la Universitat de Barcelona (IQTC-UB), Departament de Ciència de Materiales i Química Física de la Universitat de Barcelona, C/Martí I Franqués 1, 08028 Barcelona, Spain. E-mail: s.bromley@ub.edu
^{b}Institució Catalana de Recerca i Estudis Avançats (ICREA), 08010 Barcelona, Spain
First published on 31st August 2016
Triarylmethyls (TAMs) are prominent highly attractive open shell organic molecular building blocks for materials science, having been used in breakthrough syntheses of organic magnetic polymers and metal organic frameworks. With their radical π-conjugated nature and a proven capacity to possess high stability via suitable chemical design, TAMs display a variety of desirable characteristics which can be exploited for a wide range of applications. Due to their particular molecular and electronic structure, the spin localization in TAMs almost entirely depends on the dihedral angles of their three aryl rings with respect to the central methyl carbon atom plane, which opens up the possibility of controlling their fundamental properties by twisting the three aryl rings. Aryl ring twist angles can be tuned to a single value by specific chemical functionalisation but controlling them by external means in organic materials or devices represents a challenging task which has not yet been experimentally achieved. Herein, through rational chemical design we propose two 2D covalent organic frameworks (2D-COFs) based on specific TAM building blocks. By employing ab initio computational modeling we demonstrate that it is possible to externally manipulate the aryl ring twist angles in these 2D-linked TAM frameworks by external mechanical means. Furthermore, we show this structural manipulation allows for finely tuning the most important characteristics of these materials such as spin localization, optical electronic transitions and magnetic interactions. Due to the enormous technological potential offered by this new class of material and the fact that our work is guided by real advances in organic materials synthesis, we believe that our predictions will inspire the experimental realization of radical-2D-COFs with externally controllable characteristics.
Fig. 1 Generic molecular structure of TAMs. φ_{1}, φ_{2} and φ_{3} represent the aryl ring twist angles mainly determining the spin density distribution associated with the unpaired electron (green) in TAMs.^{11} |
Recently, we used ab initio computational modelling to show that, due to the π-conjugated nature of the unpaired electron in TAMs, the degree of spin delocalization in this type of molecule almost entirely depends on the dihedral angle of each aryl ring with respect to the central methyl carbon atom plane (φ_{1}, φ_{2}, φ_{3} in Fig. 1).^{11} Based on π-overlap arguments,^{20–23} we found that a linear relationship exists between the spin density on αC and the average of the cosine squared of the three dihedrals (〈cos^{2}φ〉 = (cos^{2}φ_{1} + cos^{2}φ_{2} + cos^{2}φ_{3})/3). Moreover, we found that this relationship is independent of both the chemical functionalization of the considered TAM derivative and the system's temperature. In TAMs, the presence of the unpaired electron is responsible for their most interesting properties. Hence, due to the spin-localization structural relationship in these radicals, the manipulation of aryl ring twist angles represents a potentially powerful tool for controlling these properties in a precise manner.
Aryl ring twist angles in π-conjugated molecules are mainly determined by their chemical functionalisation through the correspondingly induced steric hindrance.^{24,25} Some works studying bi-phenyl molecular junctions have pointed to the possibility of using electric fields^{26} or electro-chemical potentials^{27} to externally manipulate dihedral angles within molecular devices. However, chemical design is, by far, the most utilized approach to finely tune this structural parameter to a certain value in order to determine the associated electrical,^{24,28} magnetic,^{29} optical^{30} and electro-chemical properties.^{25} The challenging task of externally tuning aryl ring twist angles in single-molecule devices comes from the extreme difficulties associated with: (i) controlling the orientation of the molecular units within the device^{15} and, (ii) the need to direct the external stimulus uniquely to twist aryl rings without causing other physico-chemical phenomena which may also interfere with the output characteristics of the system.^{10}
Recently, it has been experimentally demonstrated that applying a uniaxial strain onto graphene mono- and multi-layers greatly influences the electrical properties of the 2D material due to the perturbation of its honey-comb like structure.^{31–37} In a very similar way, externally applied strain onto a suitably designed TAM-based material could induce an aryl ring twist within every TAM unit composing the material's skeleton which, in turn, should open the possibility of externally controlling its most fundamental characteristics.^{11} External manipulation of aryl ring twist angles in π-conjugated materials though, as far as we know, has not yet been experimentally achieved, highlighting the challenging task that it represents. To accomplish this objective there are some important criteria regarding the structural characteristics of the TAM-based material to be designed:
(1) The main skeleton of the material must be composed of TAM building blocks in such a way that an external perturbation of the material's structure will induce a significant change of the conformation of the TAM building units composing it.
(2) The bonding that sustains the basic structure of the material must be strong enough to withstand the mechanical loads required to overcome the sterical hindrance between aryl rings within TAM units.
(3) The material's skeleton should be ordered as to ensure a homogeneous structural response when applying the external mechanical stimulus, leading to correspondingly smooth, sensitive and reproducible observable physico-chemical changes.
Among the different TAM-based functional materials reported to date there are two examples that come close to fulfilling our three criteria. The first is the organic polymer magnet based on a triphenylmethyl derivative (i.e. fully hydrogenated TAMs herein abbreviated as TPMs) reported by Rajca et al. in 2001 (Fig. 2a).^{13} In this case, the material's skeleton is composed of TPM units which are covalently bonded (criteria 1 and 2) but, as being a 1D organic polymer prepared in solution, its structure consists of a disordered organic network. As a consequence, the application of any type of external mechanical load would tend to generate different types of uncontrolled inelastic structural changes, preventing thus the use of aryl ring twist angles to uniquely determine the characteristics of the stretched material. The second interesting example is a metal–organic framework (MOF) based on a perchloro-triarylmethyl (i.e. fully chlorinated TAMs, herein abbreviated as PTMs) reported by Maspoch et al. in 2003 (Fig. 2b).^{18} In this case the PTM units also make up the main structure of the material which, in turn, possesses crystalline ordering (criteria 1 and 3). However, the bonding interactions between building blocks are of a weak electrostatic nature and, thus, the strain release would probably change the relative positions of PTM units within the framework before perturbing their inner molecular conformation (i.e. aryl ring twist angles), thus also precluding the use of this interesting material for our goal.
Fig. 2 (a) Organic polymer magnet based on the triphenylmethyl (TPM) reported in 2001.^{13} (b) Metal–organic framework based on a perchloro-triarylmethyl (PTM) reported on 2003.^{18} (c) STM image of a 2D-covalent organic framework (2D-COF) formed by the on surface self-reaction of 1,3,5-tri-iodo-benzene.^{77} (d) STM image of a tri-carboxylic PTM monolayer on copper.^{89} (e) TAM-based 2D-covalent organic framework proposed in this study. |
Over the last few years a new type of organic material which could realize the proposed three requirements has emerged: so called 2D-covalent organic frameworks^{38–41} (or 2D-COFs). 2D-COFs are organic bi-dimensional layers where initially discrete molecular building blocks occupy ordered positions within a covalent framework (Fig. 2c). The final properties of 2D-COFs entirely depend on the constituent molecular units, their connectivity (i.e. how building blocks are covalently bonded) and the utilized methodology to prepare them (preparation can be performed either in solution^{42,43} or on surface^{44,45}). 2D-COFs may be understood as a class of material which bridges between organic polymers and MOFs (2a and b), and which possesses the appropriate characteristics for our purpose: i.e. covalent bonding and crystalline ordering. A few computational studies have demonstrated the unique characteristics of polyradical 2D-COFs^{46–52} but, due to the inherent instability of open shell molecules, as far as we know, no 2D-COF whose structure is composed of radical building blocks has yet been experimentally reported. However, the demonstrated chemical persistence of TAMs both in solution^{53} and on surfaces^{54,55} (Fig. 2d), the use of TAMs as building blocks for other types of materials^{13,18} (Fig. 2a and b) and the impressive recent developments in the field of 2D-COFs^{40,56} represent strong support in favour of the possible preparation of TAM-based 2D-COFs in the near future (Fig. 2e).
In this study we have rationally designed two different 2D-COFs based on the triphenylmethyl (TPM 2D-COF) and the perchloro-triarylmethyl (PTM 2D-COF), which represent the most important families of TAM derivatives for materials science to date. By means of density functional theory (DFT) calculations we have optimized the corresponding 2D periodic model structures and studied their resulting physico-chemical characteristics. Our results show that both proposed 2D-COFs are stable planar nanoporous materials consisting of ordered arrays of linked open-shell centres. Application of a uniaxial strain to both networks gives rise to a smooth perturbation of the molecular structure within every TAM unit composing the material. In agreement with our previous study on single molecule TAMs,^{11} the strain-induced aryl ring twist gives rise to a monotonic and homogeneous change in the spin localization within every radical centre. Accordingly, the electronic and magnetic properties of the materials are also found to change, demonstrating the validity of our approach to design 2D multi-functional materials with externally controllable characteristics. Further, we have carried out ab initio molecular dynamics (AIMD) simulations to test the effect of finite temperatures on the studied properties of the 2D networks. Here we observe a random small deviation in the studied variables produced by thermally activated bond vibrations. However, as we observed for isolated TAM monomers, the average value of each variable over time is clearly determined by the corresponding aryl ring twist angles. By comparing the different behaviours at finite temperatures of the differently functionalized 2D-COFs (TPM and PTM based) we also demonstrate the essential role of chemical functionality in determining the controllability of the material's properties by external manipulation of aryl ring twist angles.
To mimic the application of a uniaxial strain on our 2D materials, we performed a series of restricted optimizations (following the cascade methodology indicated above), systematically increasing one of the in-plane cell parameters (by 0–5 Å) allowing the other in-plane cell parameter and internal atomic positions to relax. We note that due to the symmetry of the considered 2D-COF, increasing either the a or b cell parameter is equivalent. Biaxial isotropic strain was also investigated by applying the same degree of distortion to both a and b cell parameters simultaneously and optimizing all the internal independent coordinates. This kind of isotropic distortion was found to directly affect the sigma-bonded skeleton and the energy variation was significantly larger than the separate a (or b) uniaxial strain, even for small distortions. Results for biaxial strain are given in the ESI.† Herein we will focus our study on the uniaxial distorted structures.
The thermally-activated response of the materials was modelled via ab initio molecular dynamics (AIMD) simulations. Calculations were run for 5 ps (1 ps of equilibration plus 4 ps of production) at 300 K with fixed lattice parameters using the Bussi–Donadio–Parrinello^{60} thermostat and the hybrid PBE0 functional.
All the above described DFT-based calculations were performed using the FHI-AIMS code,^{61,62} employing the ferromagnetic solution of the system. All reported atomically partitioned spin populations were calculated using the Hirshfeld method.^{63}
For the calculation of the relevant (i.e.: non zero) magnetic coupling constants, we used the CRYSTAL09 code^{64,65} and the B3LYP^{66} hybrid density functional based method with a standard 6-21G* basis set for all atoms.^{67,68} ITOL values have been fixed to 7,7,7,7,14 to force stringent numerical convergence of energy and gradients and a shrinking factor of 6 within the Monkhorst–Pack scheme to define the reciprocal space mesh in the 2D irreducible Brillouin zone leading to a set of 20 k-points.
The description of the magnetic properties of a system of interacting localised magnetic moments is based on the Heisenberg–Dirac–Van Vleck (HDVV) spin Hamiltonian^{69}
(1) |
(2) |
Based on our designed structure (3a) we constructed two different networks with R = H (TPM 2D-COF) and R = Cl (PTM 2D-COF) whose atom coordinates and cell parameters were fully optimized (see Methodology section). We note that both our networks were constructed with all TAM units having the same helicity of their propeller-like structure; i.e. with all aryl rings within the unit cell twisted in the same direction. Experimentally, it is possible that networks may be obtained with some TAM centres possessing different helicities (i.e. plus and minus). However, for the properties studied in this work we do not expect this small conformal variability to significantly affect our predictions. As shown in Fig. 3b and c, after optimization, both the TPM and PTM 2D-COFs become, effectively, planar networks with stable radical centres placed in ordered positions within the 2D nano-porous array (pore radius: 0.7 nm). Benzene rings and αC radical centres (purple and red coloured in Fig. 3a, respectively) lay parallel to the network plane, whereas the corresponding aryl rings are partially twisted by 35° and 53° for the TPM and PTM 2D-COFs, respectively (see y view in Fig. 3b and c top). These angle values are slightly higher than those ones for the corresponding single molecules (34° and 49°, respectively), which may be attributed to the slightly higher steric congestion within the 2D materials (especially for the PTM case).
Examining the spin distribution (see Fig. 3b and c bottom), it can be seen that all TAM units within both 2D-COFs present spin density (alpha = green, beta = purple) which is mainly located on every αC with very little, or no spin population on the connecting benzene rings. Moreover, the TAM units are, in both cases, ferromagnetically coupled (i.e. all unpaired electrons present the same spin alignment). These two features are characteristic of non-Kekulé structures which, besides avoiding electron–electron pairing (i.e. quinoidization), induce FM spin alignments.^{78} The FM solutions lay 16.2 and 1.6 meV (per unit cell) below the first excited AFM solutions for the TPM and PTM 2D-COFs, respectively. The lower FM stability for the PTM 2D-COF may be associated with a weaker interaction between its unpaired electrons due to the much more pronounced twisted conformation of its aryl rings (〈φ_{i}〉 = 53°) and the subsequent suppression of spin density on the connecting benzene units. We also note that the closed-shell solutions for the TPM and PTM 2D-COFs are 1.5 and 2.2 eV higher (per unit cell) than the corresponding open shell FM solutions.
The obtained bandgap values of 2.17 and 2.86 eV for the TPM and PTM 2D-COFs, respectively, are indicative of a localization of the electronic states near the Fermi energy level. This suggests that unpaired electrons are not delocalized through the whole π-conjugated electronic system, but rather confined within each corresponding TAM unit. Therefore, our designed TPM and PTM 2D-COFs are not in-plane electrical conductors as other 2D covalent π-conjugated organic materials (such as graphene, where a local perturbation has an effect on the whole π electronic system) but, more likely, covalent frameworks composed of electronically isolated TAM radical centres. However, due to the electro-active nature of TAMs already experimentally exploited for the preparation of non-volatile memory devices,^{8} these TAM-based 2D-COFs may be potentially used as flexible data storage materials where every TAM centre within the ordered array may behave as a logic bit with two or more redox states and different output signals.
Fig. 4 Fig. 4 shows the z and y views of the relaxed structures (top) and most stretched (+22%) structures (bottom) for the TPM (4a) and PTM (4b) 2D-COFs. As it can be observed, uniaxially stretching the networks effectively generates a twist on some of the aryl rings of the constituent TAMs. Specifically, for both cases, two of the three aryl rings within each TAM unit twist towards more orthogonal positions (1 and 2 in 4a and b) with respect to the 2D network's plane, whereas the third one becomes more flattened respectively (3 in 4a and b, see twist angle values below each conformer). This twisting effect arises from the compression of the structure of the material in the in-plane direction perpendicular to the applied strain and the subsequent increase in steric hindrance between adjacent aryl rings. In Fig. 4c it can be seen that the energetic cost associated with uniaxially stretching the 2D-COFs by more than +20% of their initial size is rather modest (ca. 0.1 eV per atom). Graphene monolayers, for comparison, suffer irreversible ruptures when being uniaxially stretched by +3%.^{32,34,37} Moreover, our calculated Young's modulus for both networks (58 and 55 GPa for TPM and PTM 2D-COFs, respectively, as calculated for in-plane strains covering 0–15%) is around 20 times smaller than that of graphene (1024 GPa),^{84} which corroborates the relatively high elasticity of our designed covalent networks. This property may be understood by their nanoporous structure which offers a number of mechanical relaxation modes (i.e. aryl ring bending, twisting, bond lengthening, etc.) upon stretching as compared to the perfectly flat and densely tessellated graphene. We also note that, with uniaxial strain, the nanopores of the 2D-COFs become elongated parallel to the stretched direction (1 nm) and narrower with respect to the perpendicular direction (0.3 nm), as compared to the more radially symmetric pores (0.7 nm diameter) of the unstrained structures.
In Fig. 4d we plot the average twist angle of all aryl rings within the unit cell with respect to the corresponding central αC plane (〈φ_{i}〉) against uniaxial strain. Here it can be seen that, although not all aryl rings rotate in the same direction (see angle values in 4a and b), their average twist angle smoothly increases during the stretching process for both networks. Moreover, as already expected, the TPM 2D-COF, with a lower steric hindrance imposed by its phenyl rings, presents a wider variation of 〈φ_{i}〉 (blue curve in 4d) as compared to the more sterically hindered PTM 2D-COF (green curve in 4d).
As it can be seen in Fig. 5c, within both 2D-COFs the average αC partitioned spin population is linearly dependent on 〈cos^{2}φ_{i}〉 for the differently stretched structures. Hence, the spin-localization structural dependence for single molecule TAMs is applicable to our mechanically stretched TAM-based 2D-COFs, confirming the power of aryl ring twist angles to control the localization of unpaired electrons within our designed materials. It is also worth noting the different responses of the two differently functionalized networks. The TPM 2D-COF, which can accommodate wider twisting angles (see Fig. 4d) presents an accordingly larger spin localization variation (blue points in 5c) as compared to the more structurally constrained chlorinated PTM 2D-COF (green points in 5c).
As it can be seen in Fig. 6a, in the relaxed conformation the π-conjugated double-bond system (πO) lays at low energies as a set of bands with non-negligible electronic dispersion, which is in accordance with its intrinsically delocalized nature. Above them and energetically well separated, the singly occupied levels (SO) corresponding to the three unpaired electrons within the unit cell appear as flat electronic bands, suggesting the lack of significant delocalization of these states. The singly unoccupied counterparts (SU) appear at higher energies showing, again, little band dispersion. The flattening of the SO and SU levels are in full accordance with the mostly localized unpaired electrons on αC positions within our networks, as measured by the corresponding spin density distribution (Fig. 5). Above the SO and SU levels, the anti-bonding unoccupied bands of the π-conjugated double-bond system (πU) appear with significant dispersion.
Upon stretching (Fig. 6b), the π-conjugated double-bond bands (πO and πU) become substantially flattened, which indicates that aryl ring twisting cuts the delocalization of the π-conjugated double bond electrons. This effect is in agreement with the experimentally demonstrated electrical conductivity dependence on twist angles in π-conjugated single-molecule devices.^{24} The singly occupied/unoccupied levels (SO and SU) do not present large changes with respect to the shape of their bands, but their relative energy difference (ΔE_{SO–SU}) increases upon stretching, as indicated by black arrows (Fig. 6a and b). Such type of energy gap variation has also been observed in differently functionalized bi-phenyl compounds where a linear correlation was obtained with the cos^{2} of the corresponding twist angles based, also, on the π-conjugation dependence on the structural parameter.^{30} Therefore, due to the effective localization of unpaired electrons by twisting aryl rings (Fig. 5), we should expect a similar correlation between ΔE_{SO–SU} and the average of aryl ring twist angles (〈cos^{2}φ_{i}〉) during the uniaxial strain application. In Fig. 6c we plot ΔE_{SO–SU} against 〈cos^{2}φ_{i}〉 (where 〈cos^{2}φ_{i}〉 = (∑cos^{2}φ_{i})/9 and φ_{i} is the twist angle of each aryl ring with respect to the corresponding αC plane) for the differently stretched conformations for both the TPM and PTM 2D-COFs (blue and green, respectively).
As it can be seen in Fig. 6c there exists an excellent correlation between ΔE_{SO–SU} and 〈cos^{2}φ_{i}〉 for the differently stretched TPM (blue) and PTM (green) 2D-COFs. In this case the link between the two parameters seems to be independent on the chemical functionalization of the networks, since for both 2D-COFs ΔE_{SO–SU} follows the same correlation line with 〈cos^{2}φ_{i}〉. We found the same phenomenon when systematically twisting aryl rings in single TAM molecules by constrained optimizations (see ESI†), confirming the generality of this powerful structural–electronic relationship for TAM-based systems. Therefore, the intimate link between the delocalization of electrons and their energies^{30} allows tuning of the electronic band structure of the material by manipulation of aryl ring twist angles (see Fig. 6a and b). In this way, it is possible for instance, to elastically tune ΔE_{SO–SU} by up to 0.7 eV within the TPM 2D-COF by simply stretching the network by a 22% of its initial size (blue points in Fig. 6c). We note that these singly occupied/unoccupied levels (SO and SU) participate in a number of important phenomena in TAM-based materials, such as colour absorption bands,^{8} fluorescence^{6} or enhanced single molecule electrical conductivities.^{7} Therefore, the results of Fig. 6 demonstrate that it would be possible to finely tune optical and electrical (e.g. sheet to substrate) properties of the material by stretching it in one direction.
(cos^{2}φ_{1}·cos^{2}θ_{1}) × (cos^{2}φ_{2}·cos^{2}θ_{2}) | (3) |
(〈cos^{2}φ_{i}〉·〈cos^{2}θ_{i}〉)^{2} | (4) |
In Fig. 7c we plot ΔE_{FM–AFM} against (〈cos^{2}φ_{i}〉·〈cos^{2}θ_{i}〉)^{2} for the uniaxially stretched TPM and PTM 2D-COFs (blue and green points, respectively). As we can see, ΔE_{FM–AFM} linearly varies with (〈cos^{2}φ_{i}〉·〈cos^{2}θ_{i}〉)^{2} for both TAM 2D-COFs due to the fundamental link between the aryl ring twist angles and the localization of unpaired electrons (Fig. 5). More specifically, as it can be seen in the ESI,† the magnetic coupling constants between first neighbours (J_{1}) split into two different new magnetic coupling constants upon stretching (J_{1b} and J_{1c}). This is a direct consequence of the uniaxial applied strain and the associated reduction in the symmetry of the triangle relating the three nearest radical centres (see ESI†), which goes from equilateral (relaxed) to isosceles (22% stretched).
More generally, as shown in Fig. 7c, for the most stretched conformations (22% strained) of both materials the ferro- and anti-ferromagnetic solutions become nearly degenerate. This corresponds to a situation that has been mechanically induced where unpaired electrons do not have any energetically preferred spin alignment. Thus the ground ferromagnetic state of both TAM 2D-COFs can be tailored by applying increasing strains, due to the corresponding twisting of aryl rings. Reports of a synthesized tri-TPM derivative, which may be understood as the analogue oligomer of our designed TPM 2D-COF (see ESI† for molecular structure), showed this derivative presents a quartet state (i.e. ferromagnetic interactions) at 93 K. Due to the very small obtained J coupling constants for both materials though, any spin alignment would be expected only below 10 K (as previously reported for TAM-based magnetic polymers^{13}) This would thus limit the potential of our currently proposed materials as organic ferromagnets. However, our results demonstrate the existence of a general, simple and powerful way to externally tune magnetic interactions in multi-radical π-conjugated 2D systems and, hence, we believe these ideas may be applied to other organic ferromagnets, or may inspire the rational design and preparation of similar organic materials with mechanically switchable magnetic properties.
In the top panels of both Fig. 8a (TPM 2D-COF) and 8b (PTM 2D-COF) we can see how the average aryl ring twist angle (i.e. 〈φ_{i}〉) vary versus time at 300 K. Dashed lines represent the corresponding values at 0 K. Thus, it can be seen that the effect of finite temperatures does not importantly vary the average aryl ring twist angle from the corresponding value at 0 K, but rather introduces a random oscillation around the most stable conformation. Moving to the spin localization (Fig. 8a and b, middle panels) and ΔE_{SO–SU} variations (Fig. 8a and b bottom panel) versus time we can see that very similar results are obtained. For both quantities the predicted values at 0 K (dashed lines) are well maintained at 300 K and thermal energy simply introduces fast random fluctuations due to thermal vibrations. By comparing the three graphs within each network we can see how the rotational oscillations of aryl rings are approximately followed by the spin localization and the ΔE_{SO–SU} fluctuations. This demonstrates that the intimate link between the structural parameter and the physico-chemical properties within both 2D-COFs is well maintained at room temperature. This behaviour was also observed for single molecule TAMs.^{11}
However, the most relevant conclusions from Fig. 8 are extracted by comparing the results from the two differently functionalized networks. In the TPM 2D-COF (Fig. 8a) thermal vibrations do not prevent distinguishing the different signals from the relaxed and most stretched conformations, as the perturbation produced by the uniaxial strain is substantially higher than the thermal noise for all studied properties. For the PTM 2D-COF (Fig. 8b) though, the mechanically induced changes on both the spin localization and ΔE_{SO–SU} are so small that they are similar in magnitude to the background noise produced by thermal fluctuations (Fig. 8b), hence making their detection over time much more difficult or even not possible. Thus, these results unambiguously demonstrate the key role of chemical functionalization in order to exploit the power of aryl ring twist angles to externally control fundamental characteristics of the materials. On the basis of these results (Fig. 8), the TPM 2D-COF appears to be a better candidate for such a purpose. However, as previously explained, it is more feasible to prepare a PTM 2D-COF in the lab due to the higher stability of the fully chlorinated TAM derivatives. Therefore, an optimal candidate to design a TAM-based 2D-COF in the lab with externally controllable characteristics would be an intermediate derivative (between the TPM and PTM) with the most appropriate balance between chemical stability and structural flexibility. Such TAMs may be found by use of appropriately substituted PTM derivatives^{86,87} or, more generally, in the extensive library of synthesized TAM derivatives.^{2}
By mechanically stretching the networks in one of the in-plane directions, we show that most important properties of the designed materials can be finely and elastically tuned. Properties such as spin localization, ferromagnetic interactions and electronic energy levels can all be smoothly controlled by applying a uniaxial strain. This behaviour relies on the manipulation of TAM's aryl rings twist angles upon stretching the organic networks, which is confirmed by the linear dependence that all the studied properties show versus the averaged cosines squared of the corresponding dihedral angles.
By performing AIMD simulations at 300 K for the relaxed and most stretched conformations of the two 2D networks, we show that steric hindrance, determined by the chemical functionalization of the networks, plays a crucial role in determining the ease by which aryl rings can be twisted, and thus the external control of materials' properties. PTMs are the most promising building blocks to be used for preparing a TAM-based 2D-COF. However, for better control of properties via strain-induced aryl ring twist angles, a TPM-based material would be preferential. To avoid the technical difficulties of preparing a stable material based on TPM derivatives, a specifically designed TAM derivative presenting radical stability and aryl ring structural freedom (i.e. between PTM and TPM) would probably be the optimal practical building block for realizing TAM-based 2D-COFs with controllable characteristics by mechanical means. Fortunately, the wide range of existing and ever increasing experimentally synthesised TAM derivatives greatly facilitates finding such a molecular candidate.
In summary, we believe that this work represents the first instance where a material has been specifically and rationally designed to enable control of its most fundamental properties by external manipulation of its aryl ring twist angles. The power of this structural parameter in π-conjugated single molecule systems has been theoretical and experimentally demonstrated.^{24,25,28–30} Our work shows how this fundamental molecular level link between structure and properties can be entirely transferred to extended materials, opening thus its potential use which may give rise to unprecedented technological applications. We hope this work becomes the starting line for the exploration of a new class of functional materials where aryl ring twist angles can be internally tuned (e.g. by chemical design) or externally varied (e.g. by mechanical strain, electric/magnetic fields, light irradiation) to yield molecular scale control of material's properties. Although of a predictive nature, our work is strongly grounded in the reported synthetic chemistry of both persistent radicals and 2D-COFs. As such we strongly believe in the chemical feasibility of our presented proposal which may inspire the experimental realization of similar extended 2D systems in the close future.
Footnote |
† Electronic supplementary information (ESI) available: (1) plot of twist angle vs. energy profiles for TPM and PTM molecules, (2) plot of bi-axial vs. uni-axial strain, (3) plots of structure and spin density for PTM 2D-COF upon strain, (4) plots of band structure of PTM 2D-COF upon strain, (5) plot of SOMO–SUMO energy difference for TPM and PTM molecules, (6) full details of magnetic coupling calculations. See DOI: 10.1039/c6sc01412g |
This journal is © The Royal Society of Chemistry 2017 |