Polyacene spacers in intramolecular magnetic coupling

We predict the intramolecular magnetic exchange coupling constant (J) for eleven nitronyl nitroxide diradicals (NN) with different linear and angular polyacene couplers from broken-symmetry density functional treatment. For the linear acene couplers, J initially decreases with increase in the number of fused rings. But from anthracene coupler onward, the J value increases with the number of benzenoid rings due to an increasing diradical character of the coupler moiety. The J value for the diradical with a fused bent coupler is always found to be smaller than that for a diradical with a linear coupler of the same size. The nuclear independent chemical shift (NICS) is calculated, and it is observed that the average of the NICS values per benzenoid ring in the diradical is less than that in the normal polyacene molecule. An empirical formula for the magnetic exchange coupling constant of a NN diradical with an aromatic spacer is obtained by combining the Wiberg bond order (BO), the angle of twist (phi) of the monoradical (NN) plane from the plane of the coupler, and the NICS values. A comparison of the formula with the computed values reveals that, from tetracene onward, the diradical nature of the linear acene couplers becomes prominent thereby leading to an increase in the ferromagnetic coupling constant. Isotropic hyperfine coupling constants are calculated by using a polarized continuum model for the diradicals in different solvents and in vacuum.     

Systematic approach to design organic magnetic molecules: strongly coupled diradicals with ethylene coupler

The intramolecular magnetic coupling constant (J) values of diradical systems linked with two monoradicals through a coupler (para-substituted phenyl acetylene (Model I), meta-substituted phenyl acetylene (Model II), ethylene (Model III)) were investigated by unrestricted density functional theory calculations. We divided eight monoradicals into α-group and β-group according to Mulliken spin density values of the connected atoms. The overall trends in the strength of magnetic interactions of diradicals were found to be identical in three different model systems. The diradicals with para-substituted phenyl acetylene coupler resulted in almost twice stronger intramolecular magnetic coupling interactions of the corresponding diradicals as compared to the meta-substituted one with opposite magnetism. NN-Ethylene-PO (nitronyl nitroxide radical coupled to phenoxyl radical via ethylene coupler) was calculated to have the strongest magnetic coupling constant with ferromagnetism, and to be even stronger (more than twice) than NN-ethylene-NN (nitronyl nitroxide diradical with ethylene coupler), which was reported to have strong antiferromagnetic interactions in a previous experiment. It was found that the spin density values of the connected atoms are closely related to the determination of magnetic interactions and J values. The spin states of the ground state in diradical systems were explained by means of the spin alternation rule.     

Broken-symmetry density functional theory investigation on bis-nitronyl nitroxide diradicals: influence of length and aromaticity of couplers

A series of nitronyl nitroxide (NN) diradicals with linear conjugated couplers and another series with aromatic couplers have been investigated by the broken-symmetry (BS) DFT approach. The overlap integral between the magnetically active orbitals in the BS state has been explicitly computed and used for the evaluation of the magnetic exchange coupling constant (J). The calculated J values are in very good agreement with the observed values in the literature. The magnitude of J depends on the length of the coupler as well as the conformation of the radical units. The aromaticity of the spacer decreases the strength of the exchange coupling constant. The SOMO-SOMO energy splitting analysis, where SOMO stands for the singly occupied molecular orbital, and the calculation of electron paramagnetic resonance (EPR) parameters have also been carried out. The computed hyperfine coupling constants support the intramolecular magnetic interactions. The nature of magnetic exchange coupling constant can also be predicted from the shape of the SOMOs as well as the spin alternation rule in the unrestricted Hartree-Fock (UHF) treatment. It is found that pi-conjugation along with the spin-polarization plays the major role in controlling the magnitude and sign of the coupling constant.     

High magnetic exchange coupling constants: a density functional theory based study of substituted schlenk diradicals

The Schlenk diradical has been known since 1915. After a detailed experimental work by Rajca, its magnetic nature has remained more or less unexplored. We have investigated by quantum chemical calculations the nature of magnetic coupling in 11 substituted Schlenk diradicals. Substitution has been considered at the fifth carbon atom of the meta-phenylene moiety. The UB3LYP method has been used to study 12 diradicals including the original one. The 6-311G(d,p) basis set has been employed for optimization of molecular geometry in both singlet and triplet states for each species. The singlet optimization has led to the optimization of the broken-symmetry structure for 10 species including the unsubstituted one. This development makes it possible to carry out further broken symmetry calculations in two ways. The triplet calculation has been done using 6-311++G(d,p) basis set and the optimized triplet geometry in both procedures. The broken symmetry calculations have used the optimized geometries of either the triplet states or the broken symmetry solutions. The first method leads to the prediction of electron paramagnetic resonance (EPR) compatible magnetic exchange coupling constant (J) in the range 517-617 cm(-1). A direct optimization of the broken symmetry geometry gives rise to a lower estimate of J, in the range of 411-525 cm(-1) and compatible with macroscopic Curie studies. The calculated J for the unsubstituted Schlenk diradical is 512 cm(-1) that can be compared with 455 cm(-1) estimated by Rajca. In both cases, introduction of groups with +M and +I effects (Ingold's notation) decreases the J value from that for the unsubstituted Schlenk diradical while -I and -M groups at the same position increases J. These trends have been explained in terms of Hammett constants, atomic spin densities, and dihedral angles.     

Towards understanding of magnetic interactions within a series of tetrathiafulvalene-pi conjugated-verdazyl diradical cation system: a density functional theory study

The intramolecular magnetic exchange coupling constants (J) for a series of tetrathiafulvalene (TTF) and verdazyl diradical cations connected by a range of pi conjugated linkers have been investigated by means of methodology based on unrestricted density functional theory. The magnetic interaction between radicals is transmitted via pi-electron conjugation for all considered compounds. The calculation of J yields strong or medium ferromagnetic coupling interactions (in the range of 56 and 300 K) for diradical cations connected by linkers with an even number of carbon atoms that are able to provide a spin polarization pathway, while antiferromagnetic coupling is predicted when linkers with an odd number of carbon atoms are employed. The topological analysis of spin density distributions have been used to reveal the effects of the spin polarization on both linkers and spin carriers. The absence of heteroatoms that impede the spin polarization pathway, and the existence of a unique spin polarization path instead of several possible competitive routes are factors which contribute to large positive J values favoring ferromagnetic interactions between the two terminal pi-radicals. The magnitude of J depends strongly on the planarity of the molecular structure of the diradical cation since a more effective orbital overlap between the two pi-systems can be achieved. Hence, the dependence of J on the torsion angle (theta) of each spin carrier has been analyzed. In this respect, our findings show that this geometrical distortion reduces largely the calculated J values for ferromagnetic couplings, leading to weak antiferromagnetic interactions for a torsion angle of 90 degrees .     

Molecular tailoring and prediction of strongly ferromagnetically coupled trimethylenemethane-based nitroxide diradicals

We have investigated the magnetic properties of four recently synthesized stable TMM-type nitroxide diradicals. Four new diradicals are proposed by tailoring one of the species in such a way that both conjugation and planarity increase. As a remarkable consequence, the intramolecular ferromagnetic exchange interaction was found to be quite high in the proposed radicals. The calculated coupling constants were in the range of +102 to +140 cm-1. The MO and spin density analysis are provided to interpret the exchange interactions. We observed the existence of intramolecular pi-pi-interactions for the species 2. This slightly increased the J value by shortening the length of the spacer between the two spin sources.     

Ab initio quantum chemical investigation of intramolecular magnetic interaction in some diradical derivatives of imino nitroxide and nitronyl nitroxide

The magnetic properties of the monoradicals 2-(4-phenyl acetylene)-4,4,5,5-tetramethyl-4,5-dihydro-1H-imidozolyl-oxyl (1) and 2-(4-phenyl acetylene)-4,4,5,5-tetramethyl-4,5-dihydro-1H-imidazole-1-oxyl-3-oxide (2) and the diradicals 2,2'-(1,2-ethynediyldi-4,1-phenylene)bis[4,4,5,5-tetramethyl-4,5-dihydro-1H-imidozolyl-oxyl] (3), 2,2'-(1,2-ethynediyldi-4,1 3,1-phenylene)bis[4,4,5,5-tetramethyl-4,5-dihydro-1H-imidozolyl-oxyl] (4), and 2,2'-(1,2-ethynediyldi-4,1 3,1-phenylene)bis[4,4,5,5-tetramethyl-4,5-dihydro-1H-imidazole-1-oxyl-3-oxide] (5) are investigated by ab initio quantum chemical methods. The rule of spin alternation in the unrestricted Hartree-Fock (UHF) method clearly shows that the radical sites are antiferromagnetically coupled in 3 and ferromagnetically coupled in 4 and 5, which is consistent with a previous experiment. The molecular geometries are optimized at Hartree-Fock levels. This is followed by single-point calculations using the density functional (UB3LYP) treatment and the multiconfigurational complete active space self-consistent field (CASSCF) methodology. Magnetic exchange coupling constants are determined from the broken-symmetry approach. The calculated J values, -3.60 cm(-1) for 3, 0.16 cm(-1) for 4, and 0.67 cm(-1) for 5, are in excellent agreement with the observed values. Because of the very large size of the diradicals 3-5, the CASSCF (10,10) calculations cannot yield realistic J values. Nevertheless, the CASSCF calculations support the antiferromagnetic nature of the magnetic coupling in 3 and the ferromagnetic nature of the coupling in 4 and 5. The existence of an intramolecular magnetic coupling in 3-5 is also confirmed through computations of the isotropic hyperfine coupling constants for monoradicals 1 and 2 as well as diradicals 3-5.     

Photomagnetic and nonlinear optical properties in cis‐trans green fluoroprotein chromophore coupled Bis‐imino nitroxide diradicals

The current study extends an earlier investigation (Bhattacharya, et al., Phys. Chem. Chem. Phys. 2012, 14, 6905) to further explore various photomagnetic and optical properties of bis‐imino nitroxide, that is, (IN)₂‐based green fluorescent protein (GFP) chromophore coupled diradicals revealing new significant features. The conversion mechanisms of selected trans‐isomers into their corresponding cis‐conformers are discussed in detailed using a number of recently‐developed density functional theory (DFT) functionals based on the Minnesota suite of DFT‐models as well as using some other DFT functionals developed earlier. To provide a more in‐depth analysis of variations in magnetic properties as trans‐conformers (singlet ground‐state) convert into their cis‐analogues (triplet ground‐state), the changes in exchange magnetic coupling constants J are compared with the variation of the selected aromaticity indices. The aromaticity indices include the nuclear independent chemical shift [NICS(0)] values calculated at the center of ring structures and the harmonic oscillator model of aromaticity. Furthermore, the investigation of static nonlinear optical response properties in the (IN)₂‐based GFP chromophore coupled diradicals reveal unusually large static first hyperpolarizabilities for these systems which is highly significant for practical applications in optics and optoelectronics. © 2015 Wiley Periodicals, Inc.     

Intramolecular ferromagnetic coupling in bis-oxoverdazyl and bis-thioxoverdazyl diradicals with polyacene spacers

We predict the intramolecular exchange coupling constant (J) for 10 different oxo- and thioxo-verdazyl-based hi-spin ground-state diradicals with linear polyacene couplers of varying length using the broken symmetry approach in an unrestricted DFT framework. The magnetic characteristics of these systems are explained using the spin-density distribution, and an analysis is made by “magnetic” orbitals. The nuclear independent chemical shift (NICS) values have been calculated for the diradicals. The average NICS(1) (1 Å above the ring surface) value per benzenoid ring increases as the size of the coupler increases. So-called ΔNICS(1) values [the difference among average NICS(1) per benzenoid ring in the coupler and the NICS(1) of the linear polyacene molecule] are correlated with J values. Bond orders and hyperfine coupling constants have also been evaluated and analyzed for the diradicals.     

Conformationally constrained, stable, triplet ground state (S = 1) nitroxide diradicals. Antiferromagnetic chains of S = 1 diradicals

Nitroxide diradicals, in which nitroxides are annelated to m-phenylene forming tricyclic benzobisoxazine-like structures, have been synthesized and characterized by X-ray crystallography, magnetic resonance (EPR and 1H NMR) spectroscopy, as well as magnetic studies in solution and in solid state. For the octamethyl derivative of benzobisoxazine nitroxide diradical, the conformationally constrained nitroxide moieties are coplanar with the m-phenylene, leading to large values of 2J (2J/k > 200 K in solution and 2J/k > 300 K in the solid state). For the diradical, in which all ortho and para positions of the m-phenylene are sterically shielded, distortion of the nitroxide moieties from coplanarity is moderate, such that the singlet-triplet gaps remain large in both solution (2J/k > 200 K) and the solid state (2J/k approximately 400-800 K), though an onset of thermal depopulation of the triplet ground state is detectable near room temperature. These diradicals have robust triplet ground states with strong ferromagnetic coupling and good stability at ambient conditions. Magnetic behavior of the nitroxide diradicals at low temperature is best fit to the model of one-dimensional S = 1 Heisenberg chains with intrachain antiferromagnetic coupling. The antiferromagnetic coupling between the S = 1 diradicals may be associated with the methyl nitroxide C-H- - -O contacts, including nonclassical hydrogen bonds. These unprecedented organic S = 1 antiferromagnetic chains are highly isotropic, compared to those of the extensively studied Ni(II)-based chains.     

The D Parameter (EPR Zero-Field Splitting) of Localized 1,3-Cyclopentanediyl Triplet Diradicals as a Measure of Electronic Substituent Effects on the Spin Densities in Para-Substituted Benzyl-Type Radicals

The zero-field splitting parameters D of the symmetrically disubstituted and unsymmetrically monosubstituted 1,3-diaryl-1,3-cyclopentanediyl triplet diradicals 1, 2 (X = p-MeO, p-Me, p-Cl, p-NH(2), p-CO(2)Me, p-CN, p-NO(2)), and 5 were determined in 2-methyltetrahydrofuran glass at 77 K. The linear plot (m = 0.558, r(2) = 0.993) of the experimental D values for the symmetrically disubstituted derivatives versus the corresponding monosubstituted ones reveals that the electronic substituent effects are additive and implies (except for the magnetic dipolar interaction) that each benzyl-type radical site acts independently in the localized diradicals. This additivity permits us to view these triplet diradicals as a composite of the two separate monoradical components and allows us to assess valuable electronic properties of benzyl-type monoradicals from the D parameter of the triplet diradical species. A theoretical analysis shows that the D parameter is a measure of the spin density rho at the benzylic positions and the inter-radical distance d in localized diradicals. A good correlation exists between the D parameter of these triplet diradicals (constant inter-radical distance d) and the EPR hyperfine coupling constants of the corresponding benzyl-type monoradicals, which establishes that the observed electronic substituent effects reflect changes in the spin densities at the radical sites. The novel DeltaD scale allows us to quantify spectroscopically the para substituent effect on the spin delocalization at the benzylic position.     

On the photomagnetism of nitronyl nitroxide, imino nitroxide, and verdazyl-substituted azobenzene

The cis- and trans-azobenzenes are known as photochromic isomers with the trans- converting into the cis-form and vice versa upon irradiation with specific wavelengths. We have quantum chemically investigated the cis- and trans-forms of substituted azobenzene diradicals, with two nitronyl nitroxides, imino nitoxides, or verdazyls at para positions and serving as monoradical centers, to determine whether they can exhibit a photoassisted magnetic crossover. Geometries of both substituted and unsubstituted molecules have been optimized by density functional (DF) method UB3LYP using the 6-311G(d,p) basis set. Optimization of the geometry of the cis isomers has required special care. Single point singlet, triplet, and broken symmetry calculations have been done using 6-311++G(3df, 3pd) basis set. The magnetic exchange coupling constants have been estimated from the broken symmetry calculations. Absorption wavelengths have been estimated for both substituted and unsubstituted species from time-dependent DF treatment using restricted spin-polarized methodology RB3LYP and 6-311++G(3df, 3pd) basis set. From the similarity in the calculated absorption wavelengths for the unsubstituted and substituted azobenzenes, and the increased oscillator strengths (f) for the substituted species, we predict that the diradical isomers would be strongly photochromic. From our triplet state and broken symmetry calculations, we predict that both the cis- and the trans-diradicals are antiferromagnetically coupled. This prediction is consistent with the spin alternation rule, and the possibility of a magnetic crossover is nonexistent for these species.     

Magnetic Modulation in Heteroallene and Heterocumulene Based tert-butyl Nitroxide Diradicals: Spin Delocalization and Conformation Play Crucial Roles

We have theoretically investigated the magnetic properties of heteroallene (>C=C=X−) and heterocumulene (>C=C=C=X−) based tert-butyl nitroxide diradicals (X is P/As). Calculation of magnetic exchange coupling constant (J) shows ferromagnetic interaction in heteroallene based diradicals. Whereas, in heterocumulene based diradicals, tuning of J value from antiferro- to ferro-magnetic state is observed from Z- to E- isomer. Delocalization of spin density from radical site to the coupler (in planar arrangement) is observed in spin distribution analysis which is also advocated by molecular orbital analysis. The typical feature of tert-butyl nitroxide radical creates spin delocalization along with spin polarization within the coupler. The J values of all the diradicals strongly depend on the dihedral angle between radical center and coupler. Magneto-structural correlation shows that the change in dihedral angle tunes the magnetic property for both the Z- and E- isomers of heterocumulenes depending on the spin accumulation on two nearby magnetic centers. The extent of spin delocalization and conformation of spin centers on the molecular axis are important for the different J values observed in our designed systems.     

Theoretical study of the magnetic behavior of ferric wheels.

Calculations of exchange coupling constants (J) based on density functional theory for eight complete, nonmodeled ferric wheels have been performed, and a comparison with the values obtained from magnetic susceptibility data is presented. The calculated J values obtained with a generalized-gradient approximation (GGA) functional are in good agreement with the experiment probably because the inclusion of pseudopotentials partially compensates the overestimation of the spin delocalization. The magnetostructural correlation obtained shows a strong dependence of the exchange coupling on both the Fe-O-Fe bond angle and the Fe-O bond distance. This correlation holds for both the ferric wheels and the alkoxo-bridged Fe(III) dinuclear complexes reported in the literature.     

Singlet�Ctriplet energy gap for trimethylenemethane, oxyallyl diradical, and related species: single- and multireference computational results

We have investigated the through-bond exchange interactions in three non-Kekulé hydrocarbon diradicals on the basis of single- and multireference coupled cluster and related broken-symmetry (BS) methods. The singlet?Ctriplet energy gap (S-T gap) and diradical characters for these species are evaluated. It is found that the spin contamination involved in the BS solutions is non-negligible and the approximate spin-projection method greatly improves the usual BS solutions. As for Mukherjee??s state-specific multireference coupled cluster (MkMRCC) computations, the size-consistent correction with the UHF localized natural orbitals (ULO) is useful to obtain the qualitatively correct 2J values.     

Density functional theoretical investigation of the aromatic nature of BN substituted benzene and four ring polyaromatic hydrocarbons

We have studied the topological and local aromaticity of BN-substituted benzene, pyrene, chrysene, triphenylene and tetracene molecules. The nucleus-independent chemical shielding (NICS), harmonic oscillator model of aromaticity (HOMA), para-delocalization index (PDI) and aromatic fluctuation index (FLU) have been calculated to quantify aromaticity in terms of magnetic and structural criteria. We find that charge separations due to the introduction of heteroatoms largely affect both the local and topological aromaticity of these molecules. Our studies show that the presence of any kind of heteroatom in the ring not only reduces the local delocalization in the six membered ring, but also affects strongly the topological aromaticity. In fact, the relative orders of the topological and local aromaticity depend strongly on the position of the heteroatoms in the structure. In general, more ring shared BN containing molecules are less aromatic than the less ring shared BN molecules. In addition our results provide evidence that the structural stability of the molecule is dominated by the σ bond rather than the π bond.     

Spin delocalization by triple-bonded functionalities in propargyl and heteropropargyl radicals,assessed from the EPR-spectral D parameter of 1,3-cyclopentanediyl triplet diradicals

The cyclopentane-1,3-diyl triplet diradicals T and T' with the triplet-bonded acetylene, cyano, and isocyano functionalities at one of the radical sites are readily prepared from the corresponding azoalkanes by photodenitrogenation in a 2-methyltetrahydrofuran (2-MTHF) matrix at 77 K. The EPR-spectral D values of these triplet diradicals show that the spin delocalizing ability of the triple-bonded pi substituent follows the order -Ctbd1;CH > -NC approximately -CN. Good correlations of the D values have been obtained with the hyperfine coupling constants (a(H)) and with the calculated spin densities (PM3/AUHF-CI method) of the corresponding monoradicals M. The propargyl-type mesomeric structure is favored over the allenyl-type contributor for all three triple-bonded functionalities; spin delocalization is less pronounced in the heteropropargyl derivatives due to the electronegativity effect of the nitrogen atom.     

Steric influence of E- or Z-monosubstituted and terminally disubstituted vinyl groups in matrix-isolated cyclopentane-1,3-diyl triplet diradicals as probed by the zero-field EPR D parameter

The zero-field D parameter of the localized E, Z, and disubstituted vinyl cyclopentane-1,3-diyl triplet diradicals V was determined at 77 K in a 2-methyltetrahydrofuran (2-MTHF) matrix. Good linear correlations were obtained with the semiempirically (PM3; r(2) = 0. 991, n = 19) and DFT (B3LYP; r(2) = 0.998, n = 7) calculated spin densities of the triplet diradicals. The D values for the disubstituted triplet diradicals V are generally larger than the corresponding monosubstituted ones and, thus, the former are less well-delocalized and thereby more poorly stabilized. For the E- and Z-diastereomeric pairs V, only marginal changes in the theoretical assessed spin densities as well as in the D values have been found. Steric effects operate and distort the conformation of the vinyl substituent in the triplet diradical V. This is adequately reproduced by theoretical calculations. For the diphenyl-substituted triplet diradical Vl, for example, they show a more or less planar alignment of the E-phenyl group and the allylic pi system (torsion angle 12 degrees ) and thus optimal delocalization of spin, whereas the Z-phenyl group is twisted about 78 degrees out of plane and therefore is not involved in the delocalization and stabilization of spin. This results in a slightly higher D value (0.0368 cm(-)(1)), and the spin is delocalized less than in the monosubstituted E-Vm one (0.0357 cm(-)(1)).     

Experimental and theoretical studies on ferromagnetically coupled metal complexes with imino nitroxides

Copper(II), zinc(II), and nickel(II) complexes with tridentate imino nitroxyl diradicals, [CuCl(bisimpy)(MeOH)](PF(6)) (1), [ZnCl(2)(bisimpy)] (2), and [NiCl(bisimpy)(H(2)O)(2)]Cl x 2H(2)O (3) (bisimpy = 2,6-bis(1'-oxyl-4',4',5',5'-tetramethyl-4',5'-dihydro-1'H-imidazol-2'-yl)pyridine), were prepared, and their magnetic properties were studied. In 1, the Cu(II) ion has a square pyramidal coordination geometry, of which the equatorial coordination sites are occupied by three nitrogen atoms from the bisimpy and a chloride ion. The coordination geometry of the Zn(II) ion in 2 can be described as a trigonal bipyramid, with two chloride ions and a bisimpy. In 3, the Ni(II) ion has a distorted octahedral coordination geometry, of which four coordination sites are coordinated by the bisimpy and chloride ion, and two water molecules occupy the remaining cis positions. Magnetic susceptibility and EPR measurements revealed that in 1 and 3 the Cu(II) and Ni(II) ions with imino nitroxyl diradicals were ferromagnetically coupled, with the coupling constants J (H = -2J(ij) summation operator S(i)S(j)) of +165(1) and 109(2) cm(-1), respectively, and the intraligand ferromagnetic interactions in 1-3 were very weak. DFT molecular orbital calculations were performed on the diradical ligand, 1, and 2 to study the spin density distribution before and after coordination to the metal ions.     

Size dependences of the diradical character and the second hyperpolarizabilities in dicyclopenta-fused acenes: relationships with their aromaticity/antiaromaticity

Using long-range corrected density functional theory, the relationships between the electronic, magnetic, and nonlinear optical properties are drawn for two families of organic compounds, the dicyclopenta-fused acenes (DPAs) and the polyacenes (PAs), containing up to N = 12 fused rings. First, the longitudinal second hyperpolarizability (γ) of singlet DPAs is significantly enhanced with increasing system size, in comparison to PAs. This behavior is associated with an increase in the longitudinal spin polarization between the terminal five-membered rings of DPAs and is consistent with previous studies where γ is maximized for intermediate diradical character. The size dependence of the diradical character is also found to cause a hump in the γ/N evolution for singlet DPAs around N = 8. In fact, in the case of singlet PAs, the diradical characters y(0) and y(1), the various magnetic properties and the γ/N values vary monotonically with N, whereas for singlet DPAs, the shielding, the magnetizability, and the γ/N values exhibit extrema near N = 8 due to the appearance of transversal spin polarization in the middle six-membered rings in addition to the longitudinal spin polarization between the terminal five-membered rings. Moreover, it is shown that for singlet DPAs the longitudinal spin polarization (characterized by y(0)) is associated with the antiaromaticity (N ≤ 3) and the slight- or non-aromaticity (N ≥ 4) of the terminal five-membered rings, whereas the appearance of transversal spin-polarization (characterized by y(1)) is associated with the decrease in the aromaticity in the inner six-membered rings as shown for large PAs. Therefore, the exceptional behaviors in singlet DPAs for small N (N < 9) are caused by the increase in diradical character y(0) correlated with the anti-aromaticity or the slight-/non-aromaticity of terminal rings and the corresponding emergence of a global aromatic character. Such a relationship between the aromaticity/antiaromaticity and the diradical character is useful for designing real open-shell NLO molecules through the control of their diradical characters.