Ultrafast symmetry-breaking charge separation in Perylenemonoimide-embedded multichromophores: impact of regioisomerism.

Symmetry-breaking charge separation (SB-CS) has recently evolved as an emerging concept offering its potential to the latest generation of organic photovoltaics. However there are several concerns that need to be addressed to reach the state-of-the-art in SB-CS chemistry, for instance, the desirable molecular geometry, interchromophoric distance and extent of electronic coupling. To shed light on those features, it is reported herein, that ortho-functionalized perylene monoimide (PMI) constituted regioisomeric dimer and trimer derivatives with varied molecular twisting and electronic conjugation have been synthesized. In steady-state photophysical studies, all the dimers and trimer derivatives exhibit a larger bathochromic shift in the emission spectra and a significant reduction of fluorescence quantum yield in polar DMF. Among the series of multichromophores, ortho- and self-coupled dimers display the strikingly different optical feature of SB-CS with a very fast charge separation rate (τCS = 80.2 ps) upon photoexcitation in DMF, which is unveiled by femtosecond transient absorption (fs-TA) studies. The SB-CS for two dimers is well-supported by the formation of PMI˙+ and PMI˙- bands in the fs-TA spectra. Further analysis of fs-TA data revealed that, among the other multichromophores the trimer also exhibits a clear charge separation, whereas SB-CS signatures are less prominent, but can not be completely disregarded, for the meta- and para-dimers. Additionally, the charge separation dynamics of those above-mentioned PMI derivatives are devoid of a kinetically favorable excimer or triplet formation. The evidence of a profound charge transfer phenomenon in the ortho-dimer is characterized by density functional theory (DFT) calculations on excited state electronic structures. The excitonic communications in the excited state electronic arrangements unravel the key role of dihedral twisting in SB-CS. The thermodynamic feasibility of CS (ΔGCS) and activation barrier (ΔG≠) of the derivatives in DMF are established from the Rehm-Weller equation and Marcus's theory, respectively. This work is an in-depth study of the effect of mutual orientation of PMIs and regioisomerism in determining sustainable guidelines for using SB-CS.

First-principles design of heavy-atom-free singlet oxygen photosensitizers for photodynamic therapy.

In photodynamic therapy (PDT) treatment, heavy-atom-free photosensitizers (PSs) are a great source of singlet oxygen photosensitizer. Reactive oxygen species (ROS) are produced by an energy transfer from the lowest energy triplet excited state to the molecular oxygen of cancer cells. To clarify the photophysical characteristics in the excited states of a few experimentally identified thionated (>C=S) molecules and their oxygenated congeners (>C=O), a quantum chemical study is conducted. This study illustrates the properties of the excited states in oxygen congeners that render them unsuitable for PDT treatment. Concurrently, a hierarchy is presented based on the utility of the lowest-energy triplet excitons of thionated compounds. Their non-radiative decay rates are calculated for reverse-ISC and inter-system crossover (ISC) processes. In addition, the vibronic importance of C=O and C=S bonds is clarified by the computation of the Huang-Rhys factor, effective vibrational mode, and reorganization energy inside the Marcus-Levich-Jörtner system. ROS generation in thionated PSs exceeds their oxygen congeners as kf ≪ kISC, where radiative decay rate is designated as kf. As a result, the current work offers a calculated strategy for analyzing the effectiveness of thionated photosensitizers in PDT.

Pressure Induced Alteration and Stabilization of Intermolecular Stacking in Square-planar Osmium tetracarbonyl.

Osmium carbonyls are well known to form stable 18-electron complexes like Os(CO)5 , Os2 (CO)9 and Os3 (CO)12 having both bridging and terminal carbonyls. For osmium tetra-carbonyl, Os(CO)4 solid-state packing significantly alters the ground-state structure. The gas-phase stable see-saw geometry converts to a square-planar structure in solid state. Highly efficient intermolecular stacking between Os(CO)4 units assists this transformation. Each Os(CO)4 molecule is stacked in a staggered orientation with respect to each other. Pressure induces a [Xe]4f14 5d6 6s2 (S=2)→[Xe]4f14 5d8 (S=0) electronic transition in osmium stabilize a square planar osmium tetra-carbonyl. Under the influence of isotropic pressure, the molecules not only come closer to each other but their relative orientations also get significantly altered. Calculations show that at P=1 GPa and above, the eclipsed orientation for the intermolecular stacking gets preferred over the staggered form. The staggered→eclipsed intermolecular stacking orientation under pressure is shown to be controlled by London dispersion interactions.

Low Potential CO2 Reduction by Inert Fe(II)-Macrobicyclic Complex: A New Concept of Cavity Assisted CO2 Activation.

The advantage of a pre-organized π-cavity of Fe(II) complex of a newly developed macrobicycle cryptand is explored for CO2 reduction by overcoming the problem of high overpotential associated with the inert nature of the cryptate. Thus, a bipyridine-centered tritopic macrobicycle having a molecular π-cavity capable of forming Fe(II) complex as well as potential for CO2 encapsulation is synthesized. The inert Fe(II)-cryptate shows much lower potential in cyclic voltammetry than the Fe(II)-tris-dimethylbipyridine (Fe-MBP) core. Interestingly, this cryptate shows electrochemical CO2 reduction at a considerably lower potential than the Fe-MBP inert core. Therefore, this study represents that a well-structured π-cavity may generate a new series of molecular catalysts for the CO2 reduction reaction (CO2 RR), even with the inert metal complexes.

Modulation of Delayed Fluorescence Guided by Conformational Effect-Mediated Thermally Enhanced Phosphorescence in Phenothiazines-Quinoline-Cl Conjugates.

Triplet energy harvesting via thermally activated delayed fluorescence (TADF) from pure organic systems has attracted great attention in organic light-emitting diodes, sensing, and photocatalysis. However, the realization of thermally enhanced phosphorescence (TEP)-guided efficient TADF with a high rate of reverse intersystem crossing (kRISC) still needs to be discovered. Herein, we report two phenothiazine-quinoline conjugates (P2QC, P2QMC) comprising two phenothiazine donors covalently attached to the chlorine-substituted quinolinyl acceptor. Spectroscopic analysis in conjunction with quantum chemistry calculations reveals that TEP in P2QC originated due to slow internal conversion from higher-lying triplet to lowest triplet (T2' → T1') of the quasi-axial (QA) conformer and TADF (kRISC = 1.44 × 108 s-1) originated from the quasi-equatorial (QE) conformer caused by a low singlet-triplet gap (ΔES1-T1 = 0.11 eV) and triplet energy transfer from QA to QE owing to the degenerate ground state of the conformers. In contrast, TADF (kRISC = 0.74 × 108 s-1) and dual phosphorescence under ambient conditions are observed in P2QMC. This study provides a sustainable guideline for developing efficient TADF emitters via conformation effects and energy transfer mechanisms.

Autonomous self-healing organic crystals for nonlinear optics.

Non-centrosymmetric molecular crystals have a plethora of applications, such as piezoelectric transducers, energy storage and nonlinear optical materials owing to their unique structural order which is absent in other synthetic materials. As most crystals are brittle, their efficiency declines upon prolonged usage due to fatigue or catastrophic failure, limiting their utilities. Some natural substances, like bone, enamel, leaf and skin, function efficiently, last a life-time, thanks to their inherent self-healing nature. Therefore, incorporating self-healing ability in crystalline materials will greatly broaden their scope. Here, we report single crystals of a dibenzoate derivative, capable of self-healing within milliseconds via autonomous actuation. Systematic quantitative experiments reveal the limit of mechanical forces that the self-healing crystals can withstand. As a proof-of-concept, we also demonstrate that our self-healed crystals can retain their second harmonic generation (SHG) with high efficiency. Kinematic analysis of the actuation in our system also revealed its impressive performance parameters, and shows actuation response times in the millisecond range.

Redox-Triggered Activation of Heavy-Atom-Free Photosensitizer and Implications in Targeted Photodynamic Therapy.

A strategy for a redox-activatable heavy-atom-free photosensitizer (PS) based on thiolated naphthalimide has been demonstrated. The PS exhibits excellent reactive oxygen species (ROS) generation in the monomeric state. However, when encapsulated in a disulfide containing bioreducible amphiphilic triblock copolymer aggregate (polymersome), the PS exhibits aggregation in the confined hydrophobic environment, which results in a smaller exciton exchange rate between the singlet and triplet excited states (TDDFT studies), and consequently, the ROS generation ability of the PS was almost fully diminished. Such a PS (in the dormant state)-loaded redox-responsive polymersome showed excellent cellular uptake and intracellular release of the PS in its active form, which enabled cell killing upon light irradiation due to ROS generation. In a control experiment involving aggregates of a similar block copolymer, but lacking the bioreducible disulfide linkage, no intracellular reactivation of the PS was noticed, highlighting the importance of stimuli-responsive polymer assemblies in the area of targeted photodynamic therapy.

Electrocatalytic Reduction of Nitrogen to Ammonia Using Tiara-like Phenylethanethiolated Nickel Cluster.

The fixing of N2 to NH3 is challenging due to the inertness of the N≡N bond. Commercially, ammonia production depends on the energy-consuming Haber-Bosch (H-B) process, which emits CO2 while using fossil fuels as the sources of hydrogen and energy. An alternative method for NH3 production is the electrochemical nitrogen reduction reaction (NRR) process as it is powered by renewable energy sources. Here, we report a tiara-like nickel-thiolate cluster, [Ni6 (PET)12 ] (where, PET=2-phenylethanethiol)] as an efficient electro-catalyst for the electrochemical NRR at ambient conditions. Ammonia (NH3 : 16.2±0.8 µg h-1 cm-2 ) was the only nitrogenous product over the potential of -2.3 V vs. Fc + /Fc with a Faradaic efficiency of 25%±1.7. Based on theoretical calculations, NRR by [Ni6 (PET)12 ] proceeds through both the distal and alternating pathways with an onset potential of -1.84 V vs. RHE (i.e., -2.46 V vs. Fc + /Fc ) which corroborates with the experimental findings.

A highly contorted push-pull naphthalenediimide dimer and evidence of intramolecular singlet exciton fission.

Singlet fission is a process by which two molecular triplet excitons are generated subsequent to the absorption of one photon. Molecules that enable singlet fission have triplet state energy at least half of the bright singlet state energy. This stringent energy criteria have challenged chemists to device new molecular and supramolecular design principles to modulate the singlet-triplet energy gap and build singlet fission systems from a wide range of organic chromophores. Herein, we report for the first time intramolecular singlet fission in the seminal naphthalenediimide (NDI) scaffold constrained in a push-pull cyclophane architecture, while individually the NDI chromophore does not satisfy the energy criterion. The challenging synthesis of this highly contorted push-pull cyclophane is possible from the preorganized pincer-like precursor. The special architecture establishes the shortest co-facial NDI⋯NDI contacts (3.084 Å) realized to date. Using broadband femtosecond transient absorption, we find that the correlated T-T pair forms rapidly within 380 fs of photoexcitation. Electronic structure calculations at the level of state-averaged CASSCF (ne,mo)/XMCQDPT2 support the existence of the multi-excitonic T-T pair state, thereby confirming the first example of singlet exciton fission in a NDI scaffold.

Understanding the Regioselectivity of Ion-Pair-Assisted Meta-Selective C(sp2)-H Activation in Conformationally Flexible Arylammonium Salts.

The lack of directionality and the long-range nature of Coulomb interactions have been a bottleneck to achieve chemically precise C-H activation using ion-pairs. Recent report by Phipps and co-workers of the ion-pair-directed regioselective Iridium-catalyzed borylation opens a new direction toward harnessing noncovalent interactions for C-H activation. In this article, the mechanism and specific role of ion-pairing are investigated using density functional theory (DFT). Computational studies reveal that meta C-H activation is kinetically more favorable than the para analogue due to stronger electrostatic interactions between the ion-pairs in closer proximity [d(NMe3+···SO3-)TSP1m = 3.93 Å versus d(NMe3+···SO3-)TSP1p = 4.30 Å]. The electrostatic interactions overwhelm the Pauli repulsion and distortion interactions incurred in bringing the oppositely charged ions in close contact for the rate-limiting meta transition state (TSP1m). Multiple linear regression shows that the free energies of activation correlate well with descriptors like the charge densities on the meta carbon and Ir atom along with that on the cation and anion with R2 = 0.74. Tuned range-separated DFT calculations demonstrate accurately the localization of charge separation in the reactant complex and transition state for the meta selectivity.

Spin Alternation Rule in USCF for Through-Bond Magnetic Coupling─A New Look: Why and When Does It Arise and How To See It.

This work takes a new look at the spin alternation rule in unrestricted self-consistent-field (USCF) calculations in terms of structural characteristics such as periodicity, impurity location, and Coulomb exchange. For clarity, the systems considered are biradicals produced from linear conjugated hydrocarbons. Both site-parametrized Hamiltonian models for theoretical analysis and spin unrestricted density functional theory (DFT) calculations are used. Theoretical analysis leads to the following conclusions: (1) The diradical state is an excited state of a linear chain of N conjugated carbon atoms (when N is about ≤ 10). Spin alternation is a consequence of the (truncated) periodic symmetry combined with filling each closed-shell pi orbital with two electrons and each singly occupied molecular orbital (SOMO) with one electron. Spin polarization is evident in triplet (T) and broken symmetry (BS) solutions for an odd N and only in the T solution for an even N. Spin alternation is visible in the BS for an odd N and always remains muted in the calculated T. (2) For a doped chain with two radical centers, spin alternation is generally visible in the BS for an odd N. The sign of spin population on the radical centers in the BS indicates the stable spin. For radical centers separated by an odd (even) number of pz electrons, spin alternation favors T (S) state with FM (AFM) interaction. Spin oscillation remains less transparent for an even N without exchange. (3) In an unrestricted treatment with exchange, spin alternation becomes observable. Without SCF iterations, the more stable state can be identified from a clear spin oscillation in the BS. An irregular oscillation indicates a possible singlet ground state. These observations are supported by density functional calculations using the B3LYP functional and the 6-311+g(d,p) basis set on linear decapentaene diradicals with nitronyl nitroxide moieties substituted on two sets of conjugated atoms, (3,9) and (3,10). Because of the SCF procedure, one finds spin alternation in the T (BS) solution and erratic oscillation in the BS (T) solution of the 3,9 (3,10) diradical in respective equilibrium geometries. The ground state is T (S). DFT adiabatic coupling constants, SOMO energies, spin population plots, and SOMO lobe diagrams compare well with molecular electronic characteristics from theoretical analysis using Hamiltonian parameters.

Harnessing Noncovalent Interactions for a Directed Evolution of a Six-Component Molecular Crystal.

Building up on weak orthogonal interactions in supramolecular chemistry, a six-component crystal is designed. Using five distinctly different noncovalent forces, namely, hydrogen bonding, halogen bonding, cation-π, anion-π, and ion-pair interactions, three six-component crystals were designed with crown-ether (I), thiourea (II), 2,3,5,6-tetrafluoro-1,4-dibromobenzene (III), lone-pair donating anion (IV), ammonium cation (V), and electron-rich aromatic ring (VI). The M06-2X functional which is highly suitable in describing other weak interactions fails for ion-pairs. Tuned range-separated (RS)-DFT calculations are found to be capable in describing the ionic interactions in molecular solids. Molecular dynamics simulations show that the predicted multicomponent crystals are stable at room temperature and reducing the ionic charges for the ion-pairs destabilizes them. The strong electrostatic interactions between the three ion-pairs, NH4+···ClO4-, NH4+···HSO4-, and NH4+···HCO3- is the primary driving force for the stabilization of the six-component crystal. Using a hybrid of strong and weak intermolecular interactions, one may generate exotic molecular complexity like n-component crystals.

Origin of Intense Luminescence from Supramolecular 2D Molecular Crystals.

Luminescence enhancement in 2D molecular crystals (2D crystals) is promising for a variety of optical applications, yet the availability is limited because of unclear mechanism and inefficient design strategy of luminescence control. Herein, the room temperature phosphorescence from micron long molecular thin free-standing 2D crystals of a mono-cyclometalated Ir(III) complex designed at the water surface is reported. A large luminescence enhancement is observed from the 2D crystals at 300 K, which is comparable with the rigidified solution at 77 K suggesting room temperature phosphorescence origin of the luminescence. In situ synchrotron grazing incidence X-ray diffraction measurements determine the constituent centered rectangular unit cells with precise molecular conformation that promotes the formation of 2D crystals. The molecular crystal design leads to a reduced singlet-triplet energy gap (ΔEST ) and mixing of singlet-triplet states by spin-orbit coupling (SOC) for efficient intersystem crossing, which explains the phosphorescence origin at room temperature and luminescence enhancement. The supramolecular assembly process provides an elegant design strategy to realize room temperature phosphorescence from 2D crystals by rigid intermolecular interactions.

Tuning of the optoelectronic properties of peptide-appended core-substituted naphthalenediimides: the role of self-assembly of two positional isomers.

This study demonstrates how the self-assembly pattern of two different and isomeric peptide-appended core-substituted naphthalenediimides (NDIs) affects the modulation of their optoelectronic properties. Two isomeric peptide-attached NDIs were synthesized, purified and characterized. Interchanging the position of attachment of the peptide units and the alkyl chains in the NDI has altered the respective self-assembling patterns of these isomeric molecules in the aggregated states. The isomer having a peptide moiety in the core position and the alkyl chain in the imide position (compound N1) forms face to face stacking or 'H' aggregates in aliphatic solvents including n-hexane, and n-decane, whereas compound N2, in which the peptide moiety is at the imide position and the alkyl chain is attached at the core position of NDI exhibits edge to edge stacking or J aggregates under the same conditions as it is evident from their UV-vis studies. The H aggregated species (obtained from N1) show inter-connected nanofibers, whereas the J aggregated species (obtained from N2) exhibit the morphology of helical nanoribbons. FT-IR and X-ray diffraction studies are in favor of the same aggregation behavior. The individual packing patterns of these two peptide-based isomers have a direct impact on their respective electrical conductivity. Interestingly, the H aggregated species shows 100 times greater current conductivity than that of the J aggregate. Moreover, it is only the H aggregated species that exhibits a photocurrent, and no such photocurrent response is observed with the J aggregates. Computational studies also support that different types of aggregation patterns are formed by these two isomeric molecules in the same solvent system. This unique example of tuning of optoelectronic behavior holds future promise for the development of new peptide-conjugated π-functional materials.

Polymorphism Dependent 9-Phosphoanthracene Derivative Exhibiting Thermally Activated Delayed Fluorescence: A Computational Investigation.

Polymorphs of anthracene derivatives exhibit diverse photophysical properties that can help to develop efficient organic-based photovoltaic devices. 10-Anthryl-9-phosphoanthracene (10-APA) shows different photophysical behaviors for the solid state due to its variety in crystalline arrangement. Herein, we investigate the ground and excited-state properties of the monomer and two different polymorphs of 10-APA from first-principles. Calculations reveal that strong spin-orbit coupling (SOC) between first excited singlet state (S1) and triplet manifolds at their S1-optimized geometries enabling the reverse intersystem crossing (RISC). The electron-vibration coupling (Huang-Rhys factor) in the excited state is the most relevant factor here. For both ISC and RISC, a similarity in Huang-Rhys factors for the molecular vibration along the π···π stacking at low-frequency region makes the rates effective. On the other side, the nonvanishing vibronic relaxation modes provide a relatively slower RISC rate in the red crystal. However, for the red crystal, small reorganization energy (λ) and large Huang-Rhys factor toward S1 → S0 conversion reduce nonradiative decay, leading to a prompt fluorescence. As the feasibility of S1 ↔ T1 conversion increases in the yellow dimer, it allows a delay in fluorescence emission, leading to thermally activated delayed fluorescence (TADF).

Paradoxical design of a serendipitous pyrazolate bridging mode: a pragmatic strategy for inducing ineluctable ferromagnetic coupling.

In this contribution we have carried out a systematic magnetostructural investigation to establish a robust one-to-one correlation between the quasi-orthogonal bridging mode of a pyrazolate ring and ferromagnetic coupling. Generating a complex with an elusive quasi-orthogonal pyrazolate bridging is a challenging task but would ineluctably result in a ferromagnetic exchange pathway. Notwithstanding the rarity, we report herein a series of bis-pyrazolato copper complexes. We have successfully exploited a so-called hypothetical-deductive model on a particular set of ligand systems that forced the pyrazolate moiety to adopt an unusual bridging mode with the M-Npz-Npz-M torsion angles in the range from 49.7° to 72.8°. The corroborating variable temperature direct current (DC) magnetic susceptibility data unequivocally confirm the ferromagnetic coupling for the complexes with the torsion angles greater than 71.37°. Furthermore, the experimental results are in excellent agreement with theoretical calculations. Based on density functional theory (DFT) calculations, again a one-to-one correspondence is made between the ligand structure and magnetic behaviour. The diradical character (y0) of the complexes is correlated with the extent of bonding interactions between the Cu centers and hence, their ferromagnetic or antiferromagnetic nature. The broken symmetry (BS) calculations on the magnetically active molecular orbitals indicate the essential magnetic behaviour of the complexes, while the EPR g-tensor calculations confirm that dx2-y2 is the magnetic orbital.

Transition-State-like Planar Structures for Amine Inversion with Ultralong C-C Bonds in Diamino-o-carborane and Diamino-o-dodecahedron.

Umbrella-like inversion of pyramidalized amines proceed through a planar transition state (TS). Stabilization of the TS through N(lone-pair) → σ*(C-C) "negative hyperconjugation" in diamino-o-carborane (1) causes rapid inversion in the amine, which results in the observation of a planarized -NH2 from the X-ray crystal structure. This proceeds through quantum mechanical tunneling across the small and narrow barrier (low pyramidalization). Tuning this secondary orbital (donor-acceptor) interaction for various derivatives of 1 and diamino-o-dodecahedron (2) provides a rational approach toward increasing dC-C to as large as 2.001, 2.011, and 1.807 Å for 1b (amino oxide-o-carborane), 1i (di-N,N-dimethylamino-o-carborane), and 2g (di-N,N-diisopropylamino-o-dodecahedron), respectively.

Remote Functionalization through Symmetric or Asymmetric Substitutions Control the Pathway of Intermolecular Singlet Fission.

Singlet fission (SF) produces two coupled triplet excitons from a high energy singlet excitation. The mechanism of SF in a variety of phenyl (-Ph) substituted pentacene is systematically studied through both ab initio and density functional theory calculations. Two classes of substitution to pentacene are considered, namely, symmetric configuration with four Ph groups (TPP) and an asymmetric configuration with two Ph groups (DPP). The positions of the singlet and triplet states are determined by calibrating the active space through state averaged complete active space self-consistent field (SA-CASSCF) calculations. The SF rates are computed based on restricted active space with single and double spin flip wave functions (RAS-SF and RAS-2SF), which are analyzed based on different intermolecular π-stacking patterns of TPP and DPP. The contribution of charge transfer (CT) state near the multiexciton (ME) state plays a significant role for SF efficiency. The role of excimer formation is supportive for ME generation [J. Am. Chem. Soc. 2016, 138, 617], and hence it is critically studied. The ME generation in TPP is a slower process and occurs through an excimer-mediated path with a large coupling between the first singlet excited state and ME state. On the other hand, DPP exhibits a relatively faster SF rate through the formation of a ME state via low-lying CT state, especially the slip-stacked dimers. The present computation elegantly demonstrates the crucial role of functional group substitution in the structure of SF active molecules in determining the efficiency of fission dynamics.

Calculation of linear and nonlinear optical properties of azobenzene derivatives with Kohn-Sham and coupled-cluster methods.

Linear polarizabilities (α) and second hyperpolarizabilities (γ) of unsubstituted azobenzenes and 'push-pull' azobenzene derivatives are investigated using Kohn-Sham theory (KST) and coupled-cluster (CC) approaches. Various standard exchange-correlation functionals as well as a non-empirically 'tuned' long-range corrected (LC) functional with range-separated exchange are used in the KST calculations. When compared to correlated ab initio calculations and measurements, the tuned functional gives accurate low-energy excitation energies, especially for charge transfer (CT) transitions, and performs well for α. Basis set and solvent effects are also studied. In contrast to expectations, but in agreement with a prior study of π-conjugated systems that do not have low-energy CT excitations, the improvements of the CT excitation energies for the push-pull π-chromophores due to tuning do not go along with clear improvements of γ toward the CC reference data, likely due to the importance of the dynamic electron correlation for this property.

Quantum chemical investigation of meta-xylylene based one-dimensional polymer chain.

We have investigated unsubstituted and methyl substituted polyradical chains of meta-xylylene by using density functional theory-broken symmetry methodology (DFT-BS). Optimization of geometry in the high-spin and low-spin states have been done at B3LYP/6-31G(d,p) and M06-2X/6-31G(d,p) levels in unrestricted methodology. Single-point calculations on the high-spin optimized geometries have been done by using the 6-311G(d,p) basis set. Each polyradical has been found to be nonplanar with a high-spin ground state. Each has a coupling constant larger than thermal energy. For each group of polyradicals, the coupling constant has been found to exponentially decrease with increase in the number of phenylene groups. The B3LYP infinite chain limit has been estimated for both the unsubstituted and substituted polyradicals. The individual inter-radical-site coupling constants have been estimated for the triradicals and tetraradicals using HDVV Hamiltonian in ORCA 3.0.1 code. These are also generally large and positive, revealing a strong intersite ferromagnetic interaction. The intersite coupling constant too decreases with increasing distance between the radical centers. Finally, we have used CRYSTAL09 package for calculations on the infinitely long one-dimensional and periodic polyradical chains. The coupling constants estimated from the periodic calculations are quite large at about 500 cm(-1) and somewhat greater than the limiting values calculated for the polyradicals with an increasing number of phenylene groups. This happens as the individual polyradicals of increasing size progressively deviate from periodicity, and thereby lessens the strength of through-bond spin-spin coupling. The calculated band gap of ∼4.5 eV indicates that the infinitely long one-dimensional chains must be ferromagnetic and electron insulators.