Recent Advances in Organic Radicals and Their Magnetism

The review presents an overview of the organic radicals that have been designed and synthesized recently, and their magnetic properties are discussed. The π-conjugated organic radicals such as phenalenyl systems, functionalized nitronylnitroxides, benzotriazinyl, bisthiazolyl, aminyl-based radicals and polyradicals, and Tetrathiafulvalene (TTF)-based H-bonded radicals have been considered. The examples show that weak supramolecular interactions play a major role in modulating the ferromagnetic and antiferromagnetic properties. The new emerging direction of zethrenes, organic polyradicals, and macrocyclic polyradicals with their attractive and discrete architectures has been deliberated. The magnetic studies delineate the singlet-triplet transitions and their corresponding energies in these organic radicals. We have also made an attempt to collate the major organic neutral radicals, radical ions and radical zwitterions that have emerged over the last century.


Introduction
Organic molecules we frequently come across in our daily lives are comprised of closed-shell electronic structures with equal numbers of electrons with up-and down-spins.Therefore, most organic compounds are diamagnetic [1].The paradigm-shifting discovery of the triphenylmethyl radical [2] by Mosses Gomberg in 1900 was a major landmark in the 20th century, which set the stage for studies of molecules with unpaired electrons, i.e., systems with open-shell electronic structures.After this discovery, several classes of organic radical systems having one or more unpaired electrons have been synthesized [3].These radical systems can be polyradicals, biradicals and diradicals.In the biradicals, the distance between the two unpaired electrons in the molecule is too large and the electron exchange interaction (J) between them is negligible or nearly negligible [4,5].On the other hand, when, in a molecule comprised of two unpaired electrons, the magnitude of the dipole-dipole interaction is large enough to produce two spin states, i.e., singlet and triplet, it is referred to as a diradical [4,6].These novel organic radical systems have initiated a wide range of highly interesting applications in redox-catalysis [7,8], living radical polymerization [9][10][11] and particularly magnetochemistry [12,13], the topic which pertains to this review.In addition, radicals find new-age applications as components in batteries [14,15], molecular spintronics [16], and the kondo effect [17], as well as in chemical biology and medicine as spin-labeling agents, in electron paramagnetic resonance (EPR) imaging and also in the diverse areas of biochemistry [18,19].
The design of pure all-organic molecular materials with magnetic properties has been one of the challenges of the last few decades.In this context, the indispensable work of Cambi and Szego on spin-crossover magnetic materials in 1930 [20] has been one of the essential reference points for chemists, physicists and materials scientists towards the generation of pure organic magnetic materials.In the pure organic compounds, as a result of the constituent light-weight elements, only the isotropic magnetic exchange interactions govern the magnetic behavior at temperatures above 0.1 K.The other magnetic interactions in organic radical systems, for example hyperfine or spin-orbit interactions, sources of magnetic anisotropies, are believed to be insignificant above this temperature.As a consequence, magnetic organic systems are described at zero applied magnetic field by the effective spin Hamiltonian approach by H = −2 ∑ J ij S i S j (1), where S ij represents the effective exchange interaction parameter for the magnetic centers i and j, having total quantum spin numbers S i and S j , respectively, with the summation operating over the adjacent pair of centers.From eqn. (1), when S ij is positive the magnetic interaction is ferromagnetic (FM), and if S ij is negative the interaction is antiferromagnetic (AFM).In general, organic magnetic materials exhibit paramagnetic properties at high temperatures.The three aspects central to designing organic molecular materials with bulk magnetic properties which should be diligently considered are: (i) the permanence of the spin-containing moieties; (ii) the magnetic interaction mechanism or the coupling routes between the neighboring spin-possessing moieties; and (iii) the transmission of the magnetic interactions along the material [21].Based on the different possibilities of the coupling routes, a FM interaction may originate if an orthogonal organization with a zero overall overlap integral between the two singly-occupied molecular orbitals (SOMOs) of the two interacting moieties is produced.In contrast, an AFM interaction results if such an integral is non-zero for non-degenerate orbitals.Magnetic interactions can take place in an intramolecular fashion between two (or more) spin-containing units within the same molecule, or in an intermolecular fashion between organic radicals governed by the isotropic exchange interactions between the unpaired electrons of the nearest molecules positioned on the respective SOMOs.Therefore, the creation of organic magnetic materials would require control over supramolecular interactions within the spin-containing sub-units, in order to promote the correct isotropic exchange interactions between the adjacent organic radicals.Supramolecular interactions such as hydrogen bonding, π-π stacking, halogen bonding [22] and bridging of radical ions through the counter-ions are the majorly applied non-covalent intermolecular interactions to assemble such molecular building blocks.
In one of the earliest reports on an organic ferromagnet, Tamura et al. in 1991 reported the existence of FM intermolecular interactions in the crystals of p-nitrophenyl-nitronylnitroxide [23], with a ferromagnetically ordered state below 0.65 K.After this pioneering report, diverse design strategies were applied towards the goal of realizing ferromagnetically ordered organic free radicals at the highest possible temperature.The result that has remained unsurpassed to date is a sulfur-nitrogen-based organic radical with a magnetic ordering at 36 K [24].
We have collated a list of diverse organic neutral radicals , radical anions , radical cations  and zwitterionic radicals [132][133][134][135][136] in Charts 1-3 according to the year of their report.It is to be noted that the list is not exhaustive.The goal has been to acquaint the reader with different classes of radicals that have emerged over time and to illustrate the growing interest towards creation of new organic radicals.Several of these radicals have been isolated and structural characterization has been carried out.Innovative design principles have led to radicals with unprecedented ambient stabilization and also a handful of radicals that can withstand even the acidic chromatography matrix during their purification.Also of tremendous interest to the synthetic radical chemistry community is to improve their opto-electronic properties, having the potential to exhibit multistability.Several of the radicals and radical ions outlined in Charts 1-3 have found tremendous applications as organic magnetic molecular materials.After presenting the collection of organic radicals, we have considered the recent (2011 until now) advances made in the design and magnetic studies of π-conjugated organic radicals such as phenalenyl systems, nitronylnitroxides, benzotriazinyl, bisthiazolyl, aminyl-based radicals and polyradicals, and H-bonded Tetrathiafulvalene (TTF)-based radicals (Sections 2 -7).The examples clearly show that weak noncovalent interactions play a major role in modulating Chart 3. Molecular structures of organic radical cations and zwitterionic radicals.

Phenalenyl-Based π-Radical
The phenalenyl (PLY) π-radical 6 is a stable hydrocarbon radical known for the last several decades [31].The stability of 6 arises from the π-electron delocalization of the unpaired electron; however, 6 could not be isolated due to its σ-dimerization.Various PLY-based organic molecular materials have been designed and synthesized, prompting experimental and theoretical studies due Chart 3. Molecular structures of organic radical cations and zwitterionic radicals.

Phenalenyl-Based π-Radical
The phenalenyl (PLY) π-radical 6 is a stable hydrocarbon radical known for the last several decades [31].The stability of 6 arises from the π-electron delocalization of the unpaired electron; however, 6 could not be isolated due to its σ-dimerization.Various PLY-based organic molecular materials have been designed and synthesized, prompting experimental and theoretical studies due to their intriguing opto-electronic and magnetic properties.The fascinating properties emerge from the ability of phenalenyl and its derivatives to form π-dimers in addition to σ-dimers.Huang et al. have reported a paramagnetic cyclo-biphenalenyl biradicaloid with chair (97) and boat conformations (98) (Figure 1a) [137].The paramagnetism emerges from the triplet of the π-dimer as the singlet states are magnetically silent.
Magnetochemistry 2016, 2, 42 6 of 35 to their intriguing opto-electronic and magnetic properties.The fascinating properties emerge from the ability of phenalenyl and its derivatives to form π-dimers in addition to σ-dimers.Huang et al.
have reported a paramagnetic cyclo-biphenalenyl biradicaloid with chair (97) and boat conformations (98) (Figure 1a) [137].The paramagnetism emerges from the triplet of the π-dimer as the singlet states are magnetically silent.Magnetic susceptibility measurements of radicals 99 and 100 determined the paramagnetic components of the magnetic susceptibility.For both radicals, the χ P follows the Curie-Weiss behavior, χ P = C/ (T − θ) and Curie constants C corresponding to 0.291 and 0.366 emu/mol, respectively, were determined (Figure 2a).The negative Weiss constants θ of −4.0 and −8.0 K of radicals 99 and 100, respectively, and the reduced spin counts realized at low temperatures pointed towards the AFM interactions concomitant to the dimeric character of the radicals (Figure 2b). to their intriguing opto-electronic and magnetic properties.The fascinating properties emerge from the ability of phenalenyl and its derivatives to form π-dimers in addition to σ-dimers.Huang et al.

Nitronylnitroxides
A diverse range of organic nitronylnitroxide-based molecules has been synthesized and their magnetic properties have been well studied by several groups.We will discuss some recent examples, which have been designed based on subtle functional group modifications, for example aminophenylnitronylnitroxides, hydroxyphenylnitronylnitroxides and biphenyl nitronylnitroxides.3a).A maximum at 40-50 K was revealed from the χ versus T plot, which is characteristic of the AFM exchange (Figure 3b).

Nitronylnitroxides
A diverse range of organic nitronylnitroxide-based molecules has been synthesized and their magnetic properties have been well studied by several groups.We will discuss some recent examples, which have been designed based on subtle functional group modifications, for example aminophenylnitronylnitroxides, hydroxyphenylnitronylnitroxides and biphenyl nitronylnitroxides.The AFM exchange interaction could be revealed as the χ versus T data for 101 was positioned below the paramagnetic Curie curve (Figure 3c).The higher temperature region of the 1/χ versus T plot demonstrated a Curie constant of 0.363 emu•K/Oe•mol, and an extrapolated Weiss constant θ = −3.4K having an alteration of the slope at about 5 K (Figure 3d).The AFM interactions at the higher temperature region were also established in the χT versus T plot, which decreased considerably in the range 50-3 K (Figure 3e).However, at 2-3 K the decrease levels out, and therefore, the susceptibility behavior has been investigated in the range of 0.3-3.6K (Figure 3f).Interestingly, competing exchange pathways were revealed in 101 as χT versus T increased after minimizing at about 2 K.A small positive Weiss constant of +0.26 K with a Curie constant of 0.2 emu•K/Oe•mol.could be determined from the 1/χ versus T plot of the 0.6-3.6K data.The AFM exchange interaction could be revealed as the χ versus T data for 101 was positioned below the paramagnetic Curie curve (Figure 3c).The higher temperature region of the 1/χ versus T plot demonstrated a Curie constant of 0.363 emu•K/Oe•mol, and an extrapolated Weiss constant θ = −3.4K having an alteration of the slope at about 5 K (Figure 3d).The AFM interactions at the higher temperature region were also established in the χT versus T plot, which decreased considerably in the range 50-3 K (Figure 3e).However, at 2-3 K the decrease levels out, and therefore, the susceptibility behavior has been investigated in the range of 0.3-3.6K (Figure 3f).Interestingly, competing exchange pathways were revealed in 101 as χT versus T increased after minimizing at about 2 K.A small positive Weiss constant of +0.26 K with a Curie constant of 0.2 emu•K/Oe•mol.could be determined from the 1/χ versus T plot of the 0.6-3.6K data.
The plot of the molar magnetization (M/Ms) versus field (H) (Figure 5a) for 103 showed a normalization to an expected field-saturated magnetization of 5585 emu/mol for S = 1/2 spin carriers.The M/Ms plot rests below the Brillouin curve as anticipated for isolated S = 1/2 spins, indicating significant intermolecular AFM exchange interactions.M/Ms vs. H is fairly linear up to 40,000 Oe at
The plot of the molar magnetization (M/Ms) versus field (H) (Figure 5a) for 103 showed a normalization to an expected field-saturated magnetization of 5585 emu/mol for S = 1/2 spin carriers.
The M/Ms plot rests below the Brillouin curve as anticipated for isolated S = 1/2 spins, indicating significant intermolecular AFM exchange interactions.M/Ms vs. H is fairly linear up to 40,000 Oe at 2.0 K with a very subtle inflection point at ~7000 Oe.In contrast, at 0.55 K the M/Ms vs. H exhibited a concave down-region from 1000 to 7000 Oe, and a concave up-region from 7000 to 40,000 Oe, as clearly seen from Figure 5a.This state is followed by a concave down-region above 40,000 Oe and finally saturation around ~50,000 Oe.
Magnetochemistry 2016, 2, 42 9 of 35 2.0 K with a very subtle inflection point at ~7000 Oe.In contrast, at 0.55 K the M/Ms vs. H exhibited a concave down-region from 1000 to 7000 Oe, and a concave up-region from 7000 to 40,000 Oe, as clearly seen from Figure 5a.This state is followed by a concave down-region above 40,000 Oe and finally saturation around ~50,000 Oe.The χdc versus T curve persists below the Curie curve, again indicating AFM exchange (Figure 5b).Furthermore, AFM exchange was also revealed from the Curie-Weiss plot of 1/χdc versus T, which produced a Weiss constant of θ = −1.4K (Figure 5c).The plot rises below 4 K monotonically and levels out below 1.5 K to give an ill-defined maximum at about 0.9 K.A zero field ac susceptibility plot over 0.35-10 K exhibited χac versus T to rise below 4 K, as is also shown in Figure 5d.There is    The χdc versus T curve persists below the Curie curve, again indicating AFM exchange (Figure 5b).Furthermore, AFM exchange was also revealed from the Curie-Weiss plot of 1/χdc versus T, which produced a Weiss constant of θ = −1.4K (Figure 5c).The plot rises below 4 K monotonically and levels out below 1.5 K to give an ill-defined maximum at about 0.9 K.A zero field ac susceptibility plot over 0.35-10 K exhibited χac versus T to rise below 4 K, as is also shown in Figure 5d.There is The χdc versus T curve persists below the Curie curve, again indicating AFM exchange (Figure 5b).Furthermore, AFM exchange was also revealed from the Curie-Weiss plot of 1/χdc versus T, which produced a Weiss constant of θ = −1.4K (Figure 5c).The plot rises below 4 K monotonically and levels out below 1.5 K to give an ill-defined maximum at about 0.9 K.A zero field ac susceptibility plot over 0.35-10 K exhibited χac versus T to rise below 4 K, as is also shown in Figure 5d.There is no indication that 103 forms a magnetically ordered state above 0.35 K since no maximum in the zero field χac versus T plot could be delineated.
have described the strong AFM interactions in biphenyl-4,4'bis(nitronylnitroxide) (104) (Chart 5) in solution as well as in the solid state [141].The strength of the AFM intra-dimer coupling constant was found to be J/K B = −14.0± 0.9 K.In addition, a magnetic field-induced transition for the field distinctly above 20 T was suggested from the inter-dimer coupling J'/k B at 1 K, which is contingent on the topology of the magnetic couplings.no indication that 103 forms a magnetically ordered state above 0.35 K since no maximum in the zero field χac versus T plot could be delineated.

Biphenyl-4,4'-bis(nitronylnitroxide)
Baumgarten et al. have described the strong AFM interactions in biphenyl-4,4'bis(nitronylnitroxide) (104) (Chart 5) in solution as well as in the solid state [141].The strength of the AFM intra-dimer coupling constant was found to be J/KB 14.0 0.9 K.In addition, a magnetic field-induced transition for the field distinctly above 20 T was suggested from the inter-dimer coupling J'/kB at 1 K, which is contingent on the topology of the magnetic couplings.

Benzotriazinyl Radicals
The  The molar susceptibility (χ) for the radicals 105 and 106 increased gradually upon cooling from 300 K and reached a broad maximum at Tmax = (50 ± 2) and (50 ± 4) K, respectively, which indicated strong antiferromagnetic interactions (Figure 6a,b).Further, χ decreased gradually down to 15 K for 105 and to 20 K for 106 before it increased subsequently at lower temperatures.The uncoupled radicals at defect sites in the lattice were thought to be responsible for this increase in χ at low temperatures.On the other hand, 107 displayed different a magnetic behavior as it followed Curie-Weiss behavior between 5 and 300 K with C = 0.338 emu•K•mol −1 close to the value expected for S = 1/2 spin (Figure 6c).The positive Weiss constant of 1.83 K indicated the presence of local FM interactions between the radicals.Upon lowering the temperature from 300 K, χT increased and attained a maximum at Tmax = (30 ± 4) K. Below this temperature, χT fell sharply, suggesting the presence of AFM interactions at low temperatures.The molar susceptibility (χ) for the radicals 105 and 106 increased gradually upon cooling from 300 K and reached a broad maximum at T max = (50 ± 2) and (50 ± 4) K, respectively, which indicated strong antiferromagnetic interactions (Figure 6a,b).Further, χ decreased gradually down to 15 K for 105 and to 20 K for 106 before it increased subsequently at lower temperatures.The uncoupled radicals at defect sites in the lattice were thought to be responsible for this increase in χ at low temperatures.On the other hand, 107 displayed different a magnetic behavior as it followed Curie-Weiss behavior between 5 and 300 K with C = 0.338 emu•K•mol −1 close to the value expected for S = 1/2 spin (Figure 6c).The positive Weiss constant of 1.83 K indicated the presence of local FM interactions between the radicals.Upon lowering the temperature from 300 K, χT increased and attained a maximum at T max = (30 ± 4) K. Below this temperature, χT fell sharply, suggesting the presence of AFM interactions at low temperatures.This group further reported the synthesis of π-stacked 1,3-diphenyl-7-(thien-2-yl)-1,4-dihydro-1,2,4-benzotriazin-4-yl (110) and the FM interactions of 1D linear chains [144].Compound 110 adopted a glided π-stacked structure of centrosymmetric dimers with alternate short and long interplanar distances of 3.48 Å and 3.52 Å, respectively.Magnetic susceptibility measurements revealed that 110 followed Curie-Weiss behavior in the 5-300 K region with C = 0.378 emu•K•mol −1 and θ = +4.72K.This result is consistent with the FM interactions between S = 1/2 radicals.
In a significant development, Koutentis et al.
reported the synthesis of 1-phenyl-3trifluoromethyl-1,4-dihydrobenzo[e]-[1,2,4]triazin-4-yl (114) [147], and it is the first example of a hydrazyl radical that demonstrates a reversible sharp spin transition between a paramagnetic and a diamagnetic phase at 58(2) K. Notably, this transition is completed within a slender temperature range of 5(1) K.The first-order transition proceeds via small changes in the intra-and inter-stack interactions between comparable structural phases.X-ray crystallography demonstrated clear changes in the intermolecular interactions at 4 K and 75 K (Figure 8).This offers significant future potential to use different stimuli to drive such phase changes.range of 5(1) K.The first-order transition proceeds via small changes in the intra-and inter-stack interactions between comparable structural phases.X-ray crystallography demonstrated clear changes in the intermolecular interactions at 4 K and 75 K (Figure 8).This offers significant future potential to use different stimuli to drive such phase changes.Casu et al. reported the synthesis of radical 116 [148].During the synthesis they considered two possible constitutional isomers, 116 and 117, due to the two different modes of fusing pyrene and Blatter radicals: 1,3-diphenyl-1,4-dihydro-1,2,4-triazin-4-yl.The DFT calculation showed significant out-of-plane twisting in the pyrene moiety and high energy in the constitutional isomer 117; thus, they selected 116.The EPR spectrum of 116 in benzene at room temperature showed hyperfine coupling.

Bis-Dithiazolyl Radicals
The first thiazolyl radical was reported in 1985 by R. Mayer [56,57].Thiazyl-based neutral radicals are highly interesting since the sulfur heteroatom augments intermolecular magnetic and electronic interactions.Among them, the radicals based on the pyridine-bridged bis-dithiazolyl framework are remarkable molecular constructs as their solid-state structures as well as their transport properties can be modulated by varying the peripheral ligands or by the replacement of sulfur by selenium.Oakley et al. have also reported the bis-diselenazolyl radicals by simple replacement of S by its heavier Se analogue.They have prepared a series of isostructural S/Secontaining variants (Chart 7: 118-121) [149].

Bis-Dithiazolyl Radicals
The first thiazolyl radical was reported in 1985 by R. Mayer [56,57].Thiazyl-based neutral radicals are highly interesting since the sulfur heteroatom augments intermolecular magnetic and electronic interactions.Among them, the radicals based on the pyridine-bridged bis-dithiazolyl framework are remarkable molecular constructs as their solid-state structures as well as their transport properties can be modulated by varying the peripheral ligands or by the replacement of sulfur by selenium.Interestingly, all the radicals were found to crystallize in non-centric space groups (Figure 9).The structures comprised of pinwheel-like clusters of radicals arranged about four centers, with each of the four radicals within the pinwheel offering the foundation for a glided, π-stacked arrangement running parallel to the c-axis.They compared the magnetic properties of the bis-diselenazolyl radical (118) with its isostructural 119, 120 and 121 radicals using the ferromagnetic absorption method.The experimental data of bis-diselenazolyl radical 118 confirmed the commencement of the long-range FM order below T C , as well as a remarkably large coercive field, H C , which was revealed by this organic ferromagnet.Interestingly, all the radicals were found to crystallize in non-centric space groups (Figure 9).The structures comprised of pinwheel-like clusters of radicals arranged about four centers, with each of the four radicals within the pinwheel offering the foundation for a glided, π-stacked arrangement running parallel to the c-axis.They compared the magnetic properties of the bis-diselenazolyl radical (118) with its isostructural 119, 120 and 121 radicals using the ferromagnetic absorption method.The experimental data of bis-diselenazolyl radical 118 confirmed the commencement of the long-range FM order below TC, as well as a remarkably large coercive field, HC, which was revealed by this organic ferromagnet.A strong AFM-coupled paramagnetic behavior could be established from the cooling curve plots of χ and χT versus T for solvated 122 (Figure 10a,b), which was measured using an external field of H = 1 kOe [150].Remarkably, a large negative θ value of −61.7 K was realized from the Curie-Weiss fit over the range 100-300 K.However, there was no substantiation for a structural or magnetic phase transition upon cooling the sample to 2 K.The high-temperature cooling curve data (with H = 1 kOe) for the unsolvated material 122 (Figure 10c,d) also suggested an AFM-coupled Curie-Weiss paramagnet with θ = −27.0K. On the other hand, upon cooling below 10 K, this magnetic material demonstrated a sharp increase in χ as well as χT.A strong AFM-coupled paramagnetic behavior could be established from the cooling curve plots of χ and χT versus T for solvated 122 (Figure 10a,b), which was measured using an external field of H = 1 kOe [150].Remarkably, a large negative θ value of −61.7 K was realized from the Curie-Weiss fit over the range 100-300 K.However, there was no substantiation for a structural or magnetic phase transition upon cooling the sample to 2 K.The high-temperature cooling curve data (with H = 1 kOe) for the unsolvated material 122 (Figure 10c,d) also suggested an AFM-coupled Curie-Weiss paramagnet with θ = −27.0K. On the other hand, upon cooling below 10 K, this magnetic material demonstrated a sharp increase in χ as well as χT.
H = 1 kOe [150].Remarkably, a large negative θ value of −61.7 K was realized from the Curie-Weiss fit over the range 100-300 K.However, there was no substantiation for a structural or magnetic phase transition upon cooling the sample to 2 K.The high-temperature cooling curve data (with H = 1 kOe) for the unsolvated material 122 (Figure 10c,d) also suggested an AFM-coupled Curie-Weiss paramagnet with θ = −27.0K. On the other hand, upon cooling below 10 K, this magnetic material demonstrated a sharp increase in χ as well as χT.Further, Oakley et al. have also reported the Br-and I-derivatives of bis-dithiazolyl radicals 124 and 125 [152].The field-cooled variable-temperature (2-300 K) χT versus T measurements for 125 at H = 1 kOe (Figure 12b) revealed strong FM exchange interactions.This was indicated by a θ value of +23.3 K obtained from a Curie-Weiss fit to the 100-300 K data.Upon cooling below 90 K, the subsequent decrease in χT signaled the onset of weaker AFM interactions leading to spin-canted antiferromagnetic ordering below T N = 35 K, as revealed by bifurcation in χT for the ZFC and FC sweeps.

Aminyl Diradicals
The nitrogen-centered aminyl-based radicals are typically short-lived, and only a handful of aminyl monoradicals are stable under ambient conditions.Recently, Rajca et al. have reported the synthesis and magnetic characterization of sterically shielded aminyl diradical 126, which possessed a triplet ground state with a ∆ comparable to that of m-xylylene 127 (Chart 8) [67].However, 126 was found to be persistent in solution at room temperature on the time scale of minutes.The triplet ground state for 126 was unambiguously confirmed by SQUID magnetometry as verified by the S = 1 paramagnetic behavior for both M as a function of H as well as the plot of χ versus T (Figure 13).A small AFM coupling due to weak intermolecular interactions in solutions of 126 was suggested based on the mean-field parameter θ ≈ −0.2 K. Notably; this is the first example of an organic diradical that is persistent in solution at room temperature up to minutes with a singlet-triplet energy gap of 10 kcal/mol.

Aminyl Diradicals
The nitrogen-centered aminyl-based radicals are typically short-lived, and only a handful of aminyl monoradicals are stable under ambient conditions.Recently, Rajca et al. have reported the synthesis and magnetic characterization of sterically shielded aminyl diradical 126, which possessed a triplet ground state with a ∆E S−T comparable to that of m-xylylene 127 (Chart 8) [67].However, 126 was found to be persistent in solution at room temperature on the time scale of minutes.The triplet ground state for 126 was unambiguously confirmed by SQUID magnetometry as verified by the S = 1 paramagnetic behavior for both M as a function of H as well as the plot of χ versus T (Figure 13).A small AFM coupling due to weak intermolecular interactions in solutions of 126 was suggested based on the mean-field parameter θ ≈ −0.2 K. Notably; this is the first example of an organic diradical that is persistent in solution at room temperature up to minutes with a singlet-triplet energy gap of 10 kcal/mol.T at H = 10 Oe.Reproduced from Reference [152].Copyright 2015, American Chemical Society.

Aminyl Diradicals
The nitrogen-centered aminyl-based radicals are typically short-lived, and only a handful of aminyl monoradicals are stable under ambient conditions.Recently, Rajca et al. have reported the synthesis and magnetic characterization of sterically shielded aminyl diradical 126, which possessed a triplet ground state with a ∆ comparable to that of m-xylylene 127 (Chart 8) [67].However, 126 was found to be persistent in solution at room temperature on the time scale of minutes.The triplet ground state for 126 was unambiguously confirmed by SQUID magnetometry as verified by the S = 1 paramagnetic behavior for both M as a function of H as well as the plot of χ versus T (Figure 13).A small AFM coupling due to weak intermolecular interactions in solutions of 126 was suggested based on the mean-field parameter θ ≈ −0.2 K. Notably; this is the first example of an organic diradical that is persistent in solution at room temperature up to minutes with a singlet-triplet energy gap of 10 kcal/mol.

Aminyl Tetraradicals
Rajca et al. also reported a series of aminyltetraradicals, 128-130, and their magnetic properties (Figure 14) [153].In 125-127, an effective 2pπ-2pπ overlap within the π-system of aminyl radicals and m-phenylenes was possible as a result of the planarity of the tetraazanonacene backbone.A quintet (S = 2) ground state could be established separated from the low spin-excited states by an energy gap of ΔE ~5 kcal•mol −1 .Interestingly, the carbon and nitrogen atoms with predominant spin densities are sterically shielded, which enhanced the stability of the tetraradicals without disturbing the planarity of the radical system.The tetraradical 129 possessed five 4-tert-alkylphenyl steric groups, and in 2-methyltetrahydrofuran at room temperature it possessed a half-life of 1 h.

Aminyl Tetraradicals
Rajca et al. also reported a series of aminyltetraradicals, 128-130, and their magnetic properties (Figure 14) [153].In 125-127, an effective 2pπ-2pπ overlap within the π-system of aminyl radicals and m-phenylenes was possible as a result of the planarity of the tetraazanonacene backbone.A quintet (S = 2) ground state could be established separated from the low spin-excited states by an energy gap of ∆E ~5 kcal•mol −1 .Interestingly, the carbon and nitrogen atoms with predominant spin densities are sterically shielded, which enhanced the stability of the tetraradicals without disturbing the planarity of the radical system.The tetraradical 129 possessed five 4-tert-alkylphenyl steric groups, and in 2-methyltetrahydrofuran at room temperature it possessed a half-life of 1 h.

Aminyl Tetraradicals
Rajca et al. also reported a series of aminyltetraradicals, 128-130, and their magnetic properties (Figure 14) [153].In 125-127, an effective 2pπ-2pπ overlap within the π-system of aminyl radicals and m-phenylenes was possible as a result of the planarity of the tetraazanonacene backbone.A quintet (S = 2) ground state could be established separated from the low spin-excited states by an energy gap of ΔE ~5 kcal•mol −1 .Interestingly, the carbon and nitrogen atoms with predominant spin densities are sterically shielded, which enhanced the stability of the tetraradicals without disturbing the planarity of the radical system.The tetraradical 129 possessed five 4-tert-alkylphenyl steric groups, and in 2-methyltetrahydrofuran at room temperature it possessed a half-life of 1 h.The SQUID magnetometry studies of 128 and 129 in 2-MeTHF verified their quintet ground states, as demonstrated by the M versus H and the χ versus T plots.The M/M sat versus H/(T − θ) plots for 128 and 129 at 1.8-5 K followed the S = 2 Brillouin function.The curvatures of the plots were independent of the radical concentration, which indicated that the measured S ≈ 2 is the ground state.Weak intermolecular AFM interactions between the tetraradicals and intra-dimer AFM exchange coupling were implied from the mean field parameter, with θ ≈ −0.5 K.

H-Bonded TTF-Based Radicals
Mori et al. reported an exceptional H-bond dynamics-based switching of electrical conductivity and magnetism in a deuterated congener, i.e., κ-D 3 catechol-fused ethylenedithiotetrathiafulvalene (Cat-EDTTTF) 2 ( 132) of an H-bonded organic conductor (Figure 15a) [154].They found κ-H and κ-D to be isostructural paramagnetic semiconductors having a dimer-Mott-type electronic structure at room temperature.Notably, only κ-D underwent a phase transition at 185 K, and it transforms to a nonmagnetic insulator with a charge-ordered electronic structure.versus T, and the solid line represents the numerical fit to the Heisenberg Hamiltonian for the dimer of S = 2 tetraradicals.Reproduced from Reference [153].Copyright 2013, American Chemical Society.
The SQUID magnetometry studies of 128 and 129 in 2-MeTHF verified their quintet ground states, as demonstrated by the M versus H and the χ versus T plots.The M/Msat versus H/(T − θ) plots for 128 and 129 at 1.8-5 K followed the S = 2 Brillouin function.The curvatures of the plots were independent of the radical concentration, which indicated that the measured S ≈ 2 is the ground state.Weak intermolecular AFM interactions between the tetraradicals and intra-dimer AFM exchange coupling were implied from the mean field parameter, with θ ≈ −0.5 K.

H-Bonded TTF-Based Radicals
Mori et al. reported an exceptional H-bond dynamics-based switching of electrical conductivity and magnetism in a deuterated congener, i.e., κ-D3 catechol-fused ethylenedithiotetrathiafulvalene (Cat-EDTTTF)2 ( 132) of an H-bonded organic conductor (Figure 15a) [154].They found κ-H and κ-D to be isostructural paramagnetic semiconductors having a dimer-Mott-type electronic structure at room temperature.Notably, only κ-D underwent a phase transition at 185 K, and it transforms to a nonmagnetic insulator with a charge-ordered electronic structure.

X-ray crystallographic analysis illustrated that the remarkable switching of the electro-physical properties initiate from D transfer or displacement within the [O•••D•••O]
−1 H-bond.This process was shown to be further complemented by electron transfer between the H-bonded Cat-EDT-TTF πsystems.This corroborated because the H-bonded D dynamics and the conducting TTF π-electron are synergistically coupled.In a pioneering experiment, single crystals of κ-D were prepared from H2Cat-EDT-TTF applying a constant current to the H-donor through a simultaneous process of H/D exchange, hydroxyl group exchange, oxidation of the TTFs, and H-bond formation.Importantly, the κ-D crystals were quantitatively deuterated and could be stored under moisture-free conditions for a long time.However, the deuteration ratio was found to be gradually lowered in a few weeks under ambient conditions, possibly as a result of the D/H exchange by atmospheric H2O.
They found that the low-temperature (LT) crystal structures of κ-D and κ-H were significantly different from each other.It was in the LT phase (~180 K) that the magnetic susceptibility (Figure 15b) and TC were observed, suggesting that the phase transition occurred in association with this crystal structure change.The crystal structure of the LT phase of κ-D at 50 K is shown in Figure 16.

X-ray crystallographic analysis illustrated that the remarkable switching of the electro-physical properties initiate from D transfer or displacement within the [O•••D•••O]
−1 H-bond.This process was shown to be further complemented by electron transfer between the H-bonded Cat-EDT-TTF π-systems.This corroborated because the H-bonded D dynamics and the conducting TTF π-electron are synergistically coupled.In a pioneering experiment, single crystals of κ-D were prepared from H 2 Cat-EDT-TTF applying a constant current to the H-donor through a simultaneous process of H/D exchange, hydroxyl group exchange, oxidation of the TTFs, and H-bond formation.Importantly, the κ-D crystals were quantitatively deuterated and could be stored under moisture-free conditions for a long time.However, the deuteration ratio was found to be gradually lowered in a few weeks under ambient conditions, possibly as a result of the D/H exchange by atmospheric H 2 O.
They found that the low-temperature (LT) crystal structures of κ-D and κ-H were significantly different from each other.It was in the LT phase (~180 K) that the magnetic susceptibility (Figure 15b) and T C were observed, suggesting that the phase transition occurred in association with this crystal structure change.The crystal structure of the LT phase of κ-D at 50 K is shown in Figure 16.

Polycyclic Aromatic Hydrocarbon Radicals
The polycyclic aromatic hydrocarbons (PAHs) with an open-shell configuration denote a unique type of hydrocarbon radical comprised of one or more π-electrons which are not firmly paired into the bonding molecular orbital in the ground state.The class of PAHs with two weakly bound πelectrons is referred to as biradicals.The two unpaired electrons can be either in a low-spin singlet or in a high-spin triplet state.A thermally excited transition from the singlet to triplet state can result in a paramagnetic property, although the molecule is a singlet in the ground state.This is only possible if the singlet-triplet energy gap (∆ ) is adequately small.Several types of stable open-shell PAHs have been synthesized: (1) quinodimethane derivatives; (2) phenalenyls; (3) teranthene; and (4) zethrenes.Molecules sustaining a characteristic resonance structure between closed shell and open shell, known as Chichibabin's hydrocarbon, have also been studied in detail.
Wu et al reported derivatives of Chichibabin's hydrocarbon, 133-CS/133-OS and 134-CS/134-OS (Figure 17) [155].The time-dependent UV-vis-NIR absorption spectra of the 133-OS biradical showed that the absorption spectrum of the biradical diminished with time and the final spectrum after the decay entirely matched with the stable quinoidal form 133-CS (Figure 18a).Further, 1 HNMR spectroscopy also nicely delineated this slow relaxation process and established a transition from an unstable 133-OS to a stable 133-CS.On the other hand, 134-OS exhibited a significantly red-shifted

Polycyclic Aromatic Hydrocarbon Radicals
The polycyclic aromatic hydrocarbons (PAHs) with an open-shell configuration denote a unique type of hydrocarbon radical comprised of one or more π-electrons which are not firmly paired into the bonding molecular orbital in the ground state.The class of PAHs with two weakly bound π-electrons is referred to as biradicals.The two unpaired electrons can be either in a low-spin singlet or in a high-spin triplet state.A thermally excited transition from the singlet to triplet state can result in a paramagnetic property, although the molecule is a singlet in the ground state.This is only possible if the singlet-triplet energy gap (∆E S−T ) is adequately small.Several types of stable open-shell PAHs have been synthesized: (1) quinodimethane derivatives; (2) phenalenyls; (3) teranthene; and (4) zethrenes.Molecules sustaining a characteristic resonance structure between closed shell and open shell, known as Chichibabin's hydrocarbon, have also been studied in detail.
Wu et al reported derivatives of Chichibabin's hydrocarbon, 133-CS/133-OS and 134-CS/134-OS (Figure 17) [155].The time-dependent UV-vis-NIR absorption spectra of the 133-OS biradical showed that the absorption spectrum of the biradical diminished with time and the final spectrum after the decay entirely matched with the stable quinoidal form 133-CS (Figure 18a).Further, 1 HNMR spectroscopy also nicely delineated this slow relaxation process and established a transition from an unstable 133-OS to a stable 133-CS.On the other hand, 134-OS exhibited a significantly red-shifted absorption in the near-infraredNIR region, which was attributed to a small highest occupied molecule orbital and lowest un-occupied molecular orbital HOMO-LUMO energy gap of 1.28 eV (Figure 18b).Notably, even after cooling to −100 • C, 134-OS did not show any NMR signals at room temperature, suggesting the presence of considerable paramagnetic species.Next, the SQUID magnetometry studies were performed for 134-OS in powder form at 5-380 K (Figure 18c).The plot of χT versus temperature and the magnetic susceptibility could be nicely fitted with the Bleaney-Bowers equation.Notably, the singlet-triplet gap was found to be 2J/kB = 166 K (1.4 kJ/mol), which suggested that the ground state for 134-OS is a triplet.Importantly, 2-OS both in solution and solid states displayed very high stability under ambient air and light conditions for months.The origin of this extraordinary stabilization was ascribed to the fluorenyl and the anthracene units, accounting for the thermodynamic and kinetic stabilization, respectively.

Zethrenes
Zethrenes with six, seven and eight rings (136, 137, and 138) have been synthesized and studied for their dynamic open-shell to close-shell properties (Chart 9) [156,157] , and heptazethrenebis(dicarboximide) (140) [156].Compound 139 was found to be a closed-shell Kekule hydrocarbon as it exhibited sharp 1 H-NMR signals at room temperature due to a large singlet-triplet energy gap and only a minute quantity of triplet species.In contrast, 140 revealed a biradical character as the solution in deuterated DCM showed no NMR signals for the heptazethrene core in the aromatic region.On the other hand, sharp peaks for the heptazethrene backbone started to emerge upon cooling at −40 °C.The absorption spectrum of 140 in chloroform exhibited absorption bands from the far-red to the NIR region with maxima at 827, 747, 701, and 641 nm, further confirming its biradical character (Figure 19a).Next, the SQUID magnetometry studies were performed for 134-OS in powder form at 5-380 K (Figure 18c).The plot of χT versus temperature and the magnetic susceptibility could be nicely fitted with the Bleaney-Bowers equation.Notably, the singlet-triplet gap was found to be 2J/kB = 166 K (1.4 kJ/mol), which suggested that the ground state for 134-OS is a triplet.Importantly, 2-OS both in solution and solid states displayed very high stability under ambient air and light conditions for months.The origin of this extraordinary stabilization was ascribed to the fluorenyl and the anthracene units, accounting for the thermodynamic and kinetic stabilization, respectively.

Zethrenes
Zethrenes with six, seven and eight rings (136, 137, and 138) have been synthesized and studied for their dynamic open-shell to close-shell properties (Chart 9) [156,157] (140) [156].Compound 139 was found to be a closed-shell Kekule hydrocarbon as it exhibited sharp 1 H-NMR signals at room temperature due to a large singlet-triplet energy gap and only a minute quantity of triplet species.In contrast, 140 revealed a biradical character as the solution in deuterated DCM showed no NMR signals for the heptazethrene core in the aromatic region.On the other hand, sharp peaks for the heptazethrene backbone started to emerge upon cooling at −40 °C.The absorption spectrum of 140 in chloroform exhibited absorption bands from the far-red to the NIR region with maxima at 827, 747, 701, and 641 nm, further confirming its biradical character (Figure 19a).Next, the SQUID magnetometry studies were performed for 134-OS in powder form at 5-380 K (Figure 18c).The plot of χT versus temperature and the magnetic susceptibility could be nicely fitted with the Bleaney-Bowers equation.Notably, the singlet-triplet gap was found to be 2J/kB = 166 K (1.4 kJ/mol), which suggested that the ground state for 134-OS is a triplet.Importantly, 2-OS both in solution and solid states displayed very high stability under ambient air and light conditions for months.The origin of this extraordinary stabilization was ascribed to the fluorenyl and the anthracene units, accounting for the thermodynamic and kinetic stabilization, respectively.

Zethrenes
Zethrenes with six, seven and eight rings (136, 137, and 138) have been synthesized and studied for their dynamic open-shell to close-shell properties (Chart 9) [156,157].Recently, Wu et al. reported a series of derivatives of zethrenes: zethrenebis(dicarboximide) (139), and heptazethrenebis(dicarboximide) (140) [156].Compound 139 was found to be a closed-shell Kekule hydrocarbon as it exhibited sharp 1 H-NMR signals at room temperature due to a large singlet-triplet energy gap and only a minute quantity of triplet species.In contrast, 140 revealed a biradical character as the solution in deuterated DCM showed no NMR signals for the heptazethrene core in the aromatic region.On the other hand, sharp peaks for the heptazethrene backbone started to emerge upon cooling at −40 • C. The absorption spectrum of 140 in chloroform exhibited absorption bands from the far-red to the NIR region with maxima at 827, 747, 701, and 641 nm, further confirming its biradical character (Figure 19a).The singlet-triplet energy gap (∆ ) of 142 was determined by SQUID measurements at 5-380 K.The magnetic measurements revealed an increasing susceptibility above 220 K (Figure 19b), and a value of 2J/kB = −1946 K (3.87 kcal/mol) was obtained by fitting the data using the Bleaney-Bowers equation.A small singlet-triplet energy gap (∆ ≈ 3.87 kcal/mol) demonstrated that the singlet open-shell ground state of 142 can be thermally excited to its triplet excited state at room temperature.The percentage of the triplet species at room temperature was assessed to be ~0.14%.
Wu et al. have subsequently reported the two stable dibenzoheptazethrene derivatives 143 and 144 (Chart 9) and their biradical characters [158].The biradical characters were examined by X-ray crystallographic analysis, spectroscopic measurements, and theoretical studies.The studies afforded the understanding of how Clar′s aromatic sextet rule can be applied to the singlet biradicaloids.A distinct ESR signal, the absence of an NMR signal at room temperature and the temperaturedependent magnetic susceptibility behavior of 144 within the range 5-380 K established a singlet biradical ground state.The fitting of the data using the Bleaney-Bowers equation revealed an exchange interaction energy (2J/kB, i.e., ∆ ) of −1859.6K (−3.7 kcal/mol) (Figure 19c).

P-Quinodimethane (p-QDM)
p-QDM is known to have a large biradical character in the ground state as a result of the retrieval of the aromaticity of the benzenoid rings, which also makes them reactive.There have been focused attempts to design and synthesize stable p-QDM derivatives and its extended congeners.Wu et al. have reported a diverse range of stable oligo(N-annulated perylene)quinodimethanes nPer-CN (n = 1-6) (145)(146)(147)(148)(149)(150), respectively, comprising up to 12 para-linked benzenoid rings (Figure 21a) [159].

P-Quinodimethane (p-QDM)
p-QDM is known to have a large biradical character in the ground state as a result of the retrieval of the aromaticity of the benzenoid rings, which also makes them reactive.There have been focused attempts to design and synthesize stable p-QDM derivatives and its extended congeners.Wu et al. have reported a diverse range of stable oligo(N-annulated perylene)quinodimethanes nPer-CN (n = 1-6) (145)(146)(147)(148)(149)(150), respectively, comprising up to 12 para-linked benzenoid rings (Figure 21a) [159].From the SQUID measurements of 146 and 147 in the powder form at 5-300 K, the singlettriplet energy gap ∆ (i.e., −2J/kB) was found to be 0.342 and 0.107 kcal/mol for 146 and 147, respectively (Figure 21b).This further confirmed that both 146 and 147 have a singlet biradical ground state.A very weak coupling between the two radicals could be judged from the small ∆ (−0.064 kcal/mol) for 148.Due to the inevitable error during the data fitting for this borderline molecule, the ground state of 148 was described as a singlet biradical with a very large biradical character based on the variable-temperature ESR experiments and computational studies.Further moving to 149 and 150, the ∆ became −0.56 and −0.88 kcal/mol, respectively, signifying a triplet biradical ground state.

Tetraradicaloids
Wu et al. further reported a series of tetraradicaloid molecules (151-162) (Chart 10) [160].The SQUID measurement of 159 showed that χT increased as the temperature was increased beyond 100 K, which implied a singlet ground state.The ∆ was calculated to be −852.6K (−1.69 kcal/mol) for 159 by fitting of the data using the Bleaney-Bowers equation.
In a very recent report, Wu et al. described the first synthesis of a relatively stable nonazethrene derivative, 163 (Figure 22a) [161].They found that 163 has an open-shell ground state since the broad 1 H-NMR spectrum at room temperature got progressively sharpened by lowering the temperature.The ESR signal was found to be strong and broad with a g value of 2.0027 for both the solid and solution states, confirming its delocalized singlet diradicaloid nature.
SQUID measurements of 163 in the temperature range of 2-380 K showed that χT increased with the increase of the temperature after 250 K (Figure 22b), which correlates to a thermal population from a singlet to a paramagnetic triplet state.The singlet-triplet energy gap (∆ ) was estimated to be −5.2 kcal/mol by fitting the data using the Bleaney-Bowers equation, which was close to the calculated value of −5.4 kcal/mol.From the SQUID measurements of 146 and 147 in the powder form at 5-300 K, the singlet-triplet energy gap ∆E S−T (i.e., −2J/k B ) was found to be 0.342 and 0.107 kcal/mol for 146 and 147, respectively (Figure 21b).This further confirmed that both 146 and 147 have a singlet biradical ground state.A very weak coupling between the two radicals could be judged from the small ∆E S−T (−0.064 kcal/mol) for 148.Due to the inevitable error during the data fitting for this borderline molecule, the ground state of 148 was described as a singlet biradical with a very large biradical character based on the variable-temperature ESR experiments and computational studies.Further moving to 149 and 150, the ∆E S−T became −0.56 and −0.88 kcal/mol, respectively, signifying a triplet biradical ground state.

Tetraradicaloids
Wu et al. further reported a series of tetraradicaloid molecules (151-162) (Chart 10) [160].The SQUID measurement of 159 showed that χT increased as the temperature was increased beyond 100 K, which implied a singlet ground state.The ∆E S−T was calculated to be −852.6K (−1.69 kcal/mol) for 159 by fitting of the data using the Bleaney-Bowers equation.
In a very recent report, Wu et al. described the first synthesis of a relatively stable nonazethrene derivative, 163 (Figure 22a) [161].They found that 163 has an open-shell ground state since the broad 1 H-NMR spectrum at room temperature got progressively sharpened by lowering the temperature.
The ESR signal was found to be strong and broad with a g value of 2.0027 for both the solid and solution states, confirming its delocalized singlet diradicaloid nature.
SQUID measurements of 163 in the temperature range of 2-380 K showed that χT increased with the increase of the temperature after 250 K (Figure 22b), which correlates to a thermal population from a singlet to a paramagnetic triplet state.The singlet-triplet energy gap (∆E S−T ) was estimated to be −5.2 kcal/mol by fitting the data using the Bleaney-Bowers equation, which was close to the calculated value of −5.4 kcal/mol.Copyright 2016, American Chemical Society.

Macrocyclic Polyradicals
In a highly innovative design, Wu et al. reported two carbazole-based macrocycles, 165 and 166, constituting four and six alternatingly arranged quinoidal and aromatic carbazole units, respectively (Figure 23) [162].The quinoidal carbazole moiety has a strong propensity to recuperate two aromatic sextet rings in the diradical form, as it can be considered to be an analogue of a pro-aromatic Tschitschibabin's hydrocarbon.Therefore, in 165 and 166, two aromatic sextets can be accomplished at each transition from the closed-shell to the open-shell diradical form, to the tetraradical form, and ultimately to the hexaradical form.

Macrocyclic Polyradicals
In a highly innovative design, Wu et al. reported two carbazole-based macrocycles, 165 and 166, constituting four and six alternatingly arranged quinoidal and aromatic carbazole units, respectively (Figure 23) [162].The quinoidal carbazole moiety has a strong propensity to recuperate two aromatic sextet rings in the diradical form, as it can be considered to be an analogue of a pro-aromatic Tschitschibabin's hydrocarbon.Therefore, in 165 and 166, two aromatic sextets can be accomplished at each transition from the closed-shell to the open-shell diradical form, to the tetraradical form, and ultimately to the hexaradical form.The magnetic properties of 165 and 166 in the powder sample were studied in the temperature range of 2-380 K, where 166 showed an increase in the χT above 35 K, which suggested a singlet ground state.In addition, another clear increase of χT around 250 K was observed, which is possibly related with excitation into other spin states or a phase transition.The magnetic properties of 165 and 166 in the powder sample were studied in the temperature range of 2-380 K, where 166 showed an increase in the χT above 35 K, which suggested a singlet ground state.In addition, another clear increase of χT around 250 K was observed, which is possibly related with excitation into other spin states or a phase transition.
The magnetization measurements for 166 established a weak magnetization originating around 200 K, reaching a maximum at 250 K, and then decreasing with the increasing temperature and finally becoming weak at 380 K (Figure 24a,b).This behavior resembles a long-range ordered FM coupling.However, the narrow hysteresis loop suggests that the magnetization should be described as super-paramagnetism.In contrast, 165 showed an increase in magnetic susceptibility above 4 K and relatively weaker magnetization at 300 K (Figure 24c,d).The singlet-triplet energy gap (∆E S−T , i.e., 2J/k B ) was found to be −0.30kcal/mol for 166 and −0.25 kcal/mol for 165.The magnetization measurements for 166 established a weak magnetization originating around 200 K, reaching a maximum at 250 K, and then decreasing with the increasing temperature and finally becoming weak at 380 K (Figure 24a,b).This behavior resembles a long-range ordered FM coupling.However, the narrow hysteresis loop suggests that the magnetization should be described as superparamagnetism.In contrast, 165 showed an increase in magnetic susceptibility above 4 K and relatively weaker magnetization at 300 K (Figure 24 c,d).The singlet-triplet energy gap (∆ , i.e., 2J/kB) was found to be −0.30kcal/mol for 166 and −0.25 kcal/mol for 165.Reproduced from Reference [162].Copyright 2016, American Chemical Society.

Conclusions
In summary, we have reviewed the recent advances made in the design and magnetic studies of organic radicals and polyradicals.The π-conjugated organic radicals such as phenalenyl systems, nitronylnitroxides, benzotriazinyl, bisthiazolyl and aminyl-based radicals and polyradicals have been considered.The examples clearly show that weak noncovalent interactions play a major role in modulating ferromagnetic and antiferromagnetic properties.We have also deliberated on the emerging direction of zethrenes, organic polyradicals and macrocyclic polyradicals.The review also presents a collation of the major organic neutral radicals, radical ions and radical zwitterions that have emerged over the period of last century.
The organic radicals and radical ions have tremendous future potential in multiple fields pertaining to purely academics as well as applications and technology in the areas of energy storage [163,164], spintronics [16], quantum information processing [165] and sensing [166][167][168].For example, the realization of a flexible radical battery from organic radicals is a major step forward Reproduced from Reference [162].Copyright 2016, American Chemical Society.

Conclusions
In summary, we have reviewed the recent advances made in the design and magnetic studies of organic radicals and polyradicals.The π-conjugated organic radicals such as phenalenyl systems, nitronylnitroxides, benzotriazinyl, bisthiazolyl and aminyl-based radicals and polyradicals have been considered.The examples clearly show that weak noncovalent interactions play a major role in modulating ferromagnetic and antiferromagnetic properties.We have also deliberated on the emerging direction of zethrenes, organic polyradicals and macrocyclic polyradicals.The review also presents a collation of the major organic neutral radicals, radical ions and radical zwitterions that have emerged over the period of last century.
The organic radicals and radical ions have tremendous future potential in multiple fields pertaining to purely academics as well as applications and technology in the areas of energy storage [163,164], spintronics [16], quantum information processing [165] and sensing [166][167][168].For example, the realization of a flexible radical battery from organic radicals is a major step forward towards the storage of solar or other renewable sources of energy.Electron and nuclear spins as quantum bits (qubits) have been a focus to gain insight and contribute in the areas of quantum information science and technology.Furthermore, the attractive optical characteristics of organic radicals with multi-channel absorption bands in the UV-vis-NIR region, their conductivity and their paramagnetism provide a unique opportunity to probe and detect toxic chemicals in an effective way in comparison to the traditional ways of chemical sensing.

Chart 2 . 2 .
Molecular structures of organic radical ions.Chart Molecular structures of organic radical ions.

Figure 1 .
Figure 1.Left column depicts the valence tautomerizations between π-dimers and right column shows the valence tautomerizations between the σ-dimers of the chair 97 and boat conformations 98 of the cyclo-biphenalenyl biradicaloid molecule.The 2e/mc (multi-centered) π-π bondings are represented by the red/green/blue dashed lines/curves.The hydrogen atoms have been removed for clarity.Reproduced from Reference [137].Copyright 2011 American Chemical Society.And the molecular structure of 99 and 100.Haddon et al. have reported synthetically intricate PLYs based on the spiro-bis(3,4,6,7tetrachalcogenide-substituted phenalenyl) boron salts and the corresponding tetrathiomethyl (99)and tetrathioethyl (100)-substituted spiro-bis(phenalenyl)boron radicals having the PLY nucleus substituted [138].The radicals 99 and 100 existed as weak π-dimers in the solid state due to steric encumbrance of the thioalkyl groups in the overlaid PLY units.Magnetic susceptibility measurements of radicals 99 and 100 determined the paramagnetic components of the magnetic susceptibility.For both radicals, the χ follows the Curie-Weiss behavior, χ / θ and Curie constants C corresponding to 0.291 and 0.366 emu/mol, respectively, were determined (Figure2a).The negative Weiss constants θ of −4.0 and −8.0 K of radicals 99 and 100, respectively, and the reduced spin counts realized at low temperatures pointed towards the AFM interactions concomitant to the dimeric character of the radicals (Figure2b).

Figure 1 .
Figure 1.Left column depicts the valence tautomerizations between π-dimers and right column shows the valence tautomerizations between the σ-dimers of the chair 97 and boat conformations 98 of the cyclo-biphenalenyl biradicaloid molecule.The 2e/mc (multi-centered) π-π bondings are represented by the red/green/blue dashed lines/curves.The hydrogen atoms have been removed for clarity.Reproduced from Reference [137].Copyright 2011 American Chemical Society.And the molecular structure of 99 and 100.

Figure 1 .
Figure 1.Left column depicts the valence tautomerizations between π-dimers and right column shows the valence tautomerizations between the σ-dimers of the chair 97 and boat conformations 98 of the cyclo-biphenalenyl biradicaloid molecule.The 2e/mc (multi-centered) π-π bondings are represented by the red/green/blue dashed lines/curves.The hydrogen atoms have been removed for clarity.Reproduced from Reference [137].Copyright 2011 American Chemical Society.And the molecular structure of 99 and 100.Haddon et al. have reported synthetically intricate PLYs based on the spiro-bis(3,4,6,7tetrachalcogenide-substituted phenalenyl) boron salts and the corresponding tetrathiomethyl (99)and tetrathioethyl (100)-substituted spiro-bis(phenalenyl)boron radicals having the PLY nucleus substituted [138].The radicals 99 and 100 existed as weak π-dimers in the solid state due to steric encumbrance of the thioalkyl groups in the overlaid PLY units.Magnetic susceptibility measurements of radicals 99 and 100 determined the paramagnetic components of the magnetic susceptibility.For both radicals, the χ follows the Curie-Weiss behavior, χ / θ and Curie constants C corresponding to 0.291 and 0.366 emu/mol, respectively, were determined (Figure2a).The negative Weiss constants θ of −4.0 and −8.0 K of radicals 99 and 100, respectively, and the reduced spin counts realized at low temperatures pointed towards the AFM interactions concomitant to the dimeric character of the radicals (Figure2b).

Figure 2 .
Figure 2. (a) Plots of the magnetic susceptibilities of 99 (red) and 100 (blue).Note the black line representing the Curie-Weiss fit of the magnetic data for 100; (b) Fraction of Curie spins obtained from χT for 99 (red) and 100 (blue).Reproduced from Reference [138].Copyright 2013 American Chemical Society.

Figure 2 .
Figure 2. (a) Plots of the magnetic susceptibilities of 99 (red) and 100 (blue).Note the black line representing the Curie-Weiss fit of the magnetic data for 100; (b) Fraction of Curie spins obtained from χT for 99 (red) and 100 (blue).Reproduced from Reference [138].Copyright 2013 American Chemical Society.

Figure 3 .
Figure 3. Plots for the paramagnetic susceptibility data of 102: (a) 1/χ vs. T plot linearly fitted and (b) χ vs. T fitted with without meanfield (broken line) and including mean field (solid line the data was acquired using the 1000 Oe external field.Paramagnetic susceptibility data for 101: (c) χ vs. T plot; with noninteracting spin curie line(d) 1/χ vs. T plot depicting linear fits to data at higher (solid line) and lower (broken line) temperature regions: (e) χT vs. T plot; and (f) low temperature 1/χ vs. T plot depicting the linear fit to data.The data obtained at T > 1.8 K at 1000 Oe are shown in triangles and the data collected by ac-susceptibility with zero field at 0.6 K < T < 4 K are shown with circles.Reproduced from Reference [139].Copyright 2011, American Chemical Society.

Figure 3 .
Figure 3. Plots for the paramagnetic susceptibility data of 102: (a) 1/χ vs. T plot linearly fitted and (b) χ vs. T fitted with without meanfield (broken line) and including mean field (solid line the data was acquired using the 1000 Oe external field.Paramagnetic susceptibility data for 101: (c) χ vs. T plot; with noninteracting spin curie line(d) 1/χ vs. T plot depicting linear fits to data at higher (solid line) and lower (broken line) temperature regions: (e) χT vs. T plot; and (f) low temperature 1/χ vs. T plot depicting the linear fit to data.The data obtained at T > 1.8 K at 1000 Oe are shown in triangles and the data collected by ac-susceptibility with zero field at 0.6 K < T < 4 K are shown with circles.Reproduced from Reference [139].Copyright 2011, American Chemical Society.

Figure 4 .
Figure 4. Structure of 103 and supramolecular interactions involving the SyrNN nitronylnitroxide units in 103: 2D network formation through hydrogen-bonded OH.O-N chains, inter-radical N-O•••C chain contacts, and radical methyl CH•••ON interactions.H atoms have been removed for clarity.Reproduced from Reference [140].Copyright 2013, American Chemical Society.

Figure 5 .
Figure 5. (a) Plot of molar magnetization normalized to saturation magnetization (M/Ms) versus field at 0.55 K (upper, black data) and 2.0 K (lower, red data).The Brillouin curve for isolated S = 1/2 spins at 0.55 K and magnetic susceptibility versus temperature is depicted by the dashed lines; (b) Plot of χdc versus T at 1000 Oe external field; (c) Plot of 1/χac versus T data in zero external field (modulation 10 Oe, frequency 155 Hz); and (d) shows χac versus T data in zero external field (modulation 10 Oe, frequency 155 Hz).The solid red lines exhibits fits to an AFM 2D square planar model without mean field correction.Reproduced from Reference [140].Copyright 2013, American Chemical Society.

Figure 4 .
Figure 4. Structure of 103 and supramolecular interactions involving the SyrNN nitronylnitroxide units in 103: 2D network formation through hydrogen-bonded OH.O-N chains, inter-radical N-O•••C chain contacts, and radical methyl CH•••ON interactions.H atoms have been removed for clarity.Reproduced from Reference [140].Copyright 2013, American Chemical Society.
with a very subtle inflection point at ~7000 Oe.In contrast, at 0.55 K the M/Ms vs. H exhibited a concave down-region from 1000 to 7000 Oe, and a concave up-region from 7000 to 40,000 Oe, as clearly seen from Figure5a.This state is followed by a concave down-region above 40,000 Oe and finally saturation around ~50,000 Oe.

Figure 4 .
Figure 4. Structure of 103 and supramolecular interactions involving the SyrNN nitronylnitroxide units in 103: 2D network formation through hydrogen-bonded OH.O-N chains, inter-radical N-O•••C chain contacts, and radical methyl CH•••ON interactions.H atoms have been removed for clarity.Reproduced from Reference [140].Copyright 2013, American Chemical Society.

Figure 5 .
Figure 5. (a) Plot of molar magnetization normalized to saturation magnetization (M/Ms) versus field at 0.55 K (upper, black data) and 2.0 K (lower, red data).The Brillouin curve for isolated S = 1/2 spins at 0.55 K and magnetic susceptibility versus temperature is depicted by the dashed lines; (b) Plot of χdc versus T at 1000 Oe external field; (c) Plot of 1/χac versus T data in zero external field (modulation 10 Oe, frequency 155 Hz); and (d) shows χac versus T data in zero external field (modulation 10 Oe, frequency 155 Hz).The solid red lines exhibits fits to an AFM 2D square planar model without mean field correction.Reproduced from Reference [140].Copyright 2013, American Chemical Society.

Figure 5 .
Figure 5. (a) Plot of molar magnetization normalized to saturation magnetization (M/Ms) versus field at 0.55 K (upper, black data) and 2.0 K (lower, red data).The Brillouin curve for isolated S = 1/2 spins at 0.55 K and magnetic susceptibility versus temperature is depicted by the dashed lines; (b) Plot of χdc versus T at 1000 Oe external field; (c) Plot of 1/χac versus T data in zero external field (modulation 10 Oe, frequency 155 Hz); and (d) shows χac versus T data in zero external field (modulation 10 Oe, frequency 155 Hz).The solid red lines exhibits fits to an AFM 2D square planar model without mean field correction.Reproduced from Reference [140].Copyright 2013, American Chemical Society.

Figure 8 .
Figure 8.A view of the unit cell of 114 down the b-axis (left) and π-stacks of 114 at 4 K (middle) and depiction of the supramolecular contacts and the interplanar distances in pairs I-II and II-III at 75 K (right).H atoms have been removed for clarity.Reproduced from Reference [147].Copyright 2014, American Chemical Society.

Figure 8 .
Figure 8.A view of the unit cell of 114 down the b-axis (left) and π-stacks of 114 at 4 K (middle) and depiction of the supramolecular contacts and the interplanar distances in pairs I-II and II-III at 75 K (right).H atoms have been removed for clarity.Reproduced from Reference [147].Copyright 2014, American Chemical Society.
Oakley et al. have also reported the bis-diselenazolyl radicals by simple replacement of S by its heavier Se analogue.They have prepared a series of isostructural S/Se-containing variants (Chart 7: 118-121) [149].

Chart 7 .
electronic interactions.Among them, the radicals based on the pyridine-bridged bis-dithiazolyl framework are remarkable molecular constructs as their solid-state structures as well as their transport properties can be modulated by varying the peripheral ligands or by the replacement of sulfur by selenium.Oakley et al. have also reported the bis-diselenazolyl radicals by simple replacement of S by its heavier Se analogue.They have prepared a series of isostructural S/Secontaining variants (Chart 7: 118-121)[149].Molecular structures of bis-dithiazolyl radical molecules.

Figure 9 .
Figure 9.A view of the crystal packing of 118: (a) Parallel and (b) perpendicular to the π-stacking direction.Radicals 119, 120, and 121 were found to be isostructural to 118.Reproduced from Reference [149].Copyright 2011, American Chemical Society.

Figure 9 .
Figure 9.A view of the crystal packing of 118: (a) Parallel and (b) perpendicular to the π-stacking direction.Radicals 119, 120, and 121 were found to be isostructural to 118.Reproduced from Reference [149].Copyright 2011, American Chemical Society.

Figure 11 .
Figure 11.(a) A view of the unit cell of 123, viewed perpendicular to the yz plane, with intermolecular S3•••O' supramolecular contacts depicted with dashed lines; (b,c) Views of the ABABAB π-stacks in 120, with S1•••S3 contacts d1 and d2.Reproduced with permission from Reference [151].Copyright 2011, the Royal Society of Chemistry.Further, Oakley et al. have also reported the Br-and I-derivatives of bis-dithiazolyl radicals 124 and 125 [152].The field-cooled variable-temperature (2-300 K) χT versus T measurements for 125 at H = 1 kOe (Figure 12b) revealed strong FM exchange interactions.This was indicated by a θ value of +23.3 K obtained from a Curie-Weiss fit to the 100-300 K data.Upon cooling below 90 K, the subsequent decrease in χT signaled the onset of weaker AFM interactions leading to spin-canted antiferromagnetic ordering below TN = 35 K, as revealed by bifurcation in χT for the ZFC and FC sweeps.

Figure 11 .
Figure 11.(a) A view of the unit cell of 123, viewed perpendicular to the yz plane, with intermolecular S3•••O' supramolecular contacts depicted with dashed lines; (b,c) Views of the ABABAB π-stacks in 120, with S1•••S3 contacts d 1 and d 2 .Reproduced with permission from Reference [151].Copyright 2011, the Royal Society of Chemistry.

Figure 15 .
Figure 15.(a) Molecule 131 and 132; (b) Magnetic susceptibility data of κ-D.Cooling and heating processes are represented by blue and red circles, respectively.Black circles represent the cooling process for κ-H.The fitting curve for κ-D by the singlet-triplet dimer model with an AFM coupling of 2J/kB ≈ −600 K is represented by the orange line.Reproduced from Reference [154].Copyright 2014, American Chemical Society.

Figure 15 .
Figure 15.(a) Molecule 131 and 132; (b) Magnetic susceptibility data of κ-D.Cooling and heating processes are represented by blue and red circles, respectively.Black circles represent the cooling process for κ-H.The fitting curve for κ-D by the singlet-triplet dimer model with an AFM coupling of 2J/k B ≈ −600 K is represented by the orange line.Reproduced from Reference [154].Copyright 2014, American Chemical Society.
near-infraredNIR region, which was attributed to a small highest occupied molecule orbital and lowest un-occupied molecular orbital HOMO-LUMO energy gap of 1.28 eV (Figure18b).Notably, even after cooling to −100 °C, 134-OS did not show any NMR signals at room temperature, suggesting the presence of considerable paramagnetic species.

Figure 18 .
Figure 18.(a) UV-vis absorption spectra of the freshly generated 133-OS biradical and its timedependent changes; (b) UV-vis-NIR absorption spectrum of 134-OS and compound 135 in DCM; (c) A plot of χT versus T of 134-OS in the SQUID measurements and fitting via the Bleaney-Bowers equation.Reproduced from Reference [155].Copyright 2012, American Chemical Society.

Figure 18 .
Figure 18.(a) UV-vis absorption spectra of the freshly generated 133-OS biradical and its timedependent changes; (b) UV-vis-NIR absorption spectrum of 134-OS and compound 135 in DCM; (c) A plot of χT versus T of 134-OS in the SQUID measurements and fitting via the Bleaney-Bowers equation.Reproduced from Reference [155].Copyright 2012, American Chemical Society.

Figure 18 .
Figure 18.(a) UV-vis absorption spectra of the freshly generated 133-OS biradical and its time-dependent changes; (b) UV-vis-NIR absorption spectrum of 134-OS and compound 135 in DCM; (c) A plot of χT versus T of 134-OS in the SQUID measurements and fitting via the Bleaney-Bowers equation.Reproduced from Reference [155].Copyright 2012, American Chemical Society.

Chart 9 .
Molecular structures of the zethrene molecules.

Figure 19 .
Figure 19.(a) Absorption spectra of 139 and 140 in chloroform; (b) Plot of χT versus T for the solid 142.Open circles represent the measured data.Fitting curve was drawn using the Bleaney-Bowers equation with g = 2.00; (c) Plot of χT versus T for the solid 144.The measured data were plotted as open circles, and the fitting curve was obtained using the Bleaney-Bowers equation with g = 2.00.Reproduced from Reference [156-158].Copyright 2011, 2012, 2013, American Chemical Society.Wu et al., in a significant modification in their molecular design, reported the kinetically blocked heptazethrene derivative (141) and octazethrene derivative (142) [157].Compound 141 was found to

Figure 19 .
Figure 19.(a) Absorption spectra of 139 and 140 in chloroform; (b) Plot of χT versus T for the solid 142.Open circles represent the measured data.Fitting curve was drawn using the Bleaney-Bowers equation with g = 2.00; (c) Plot of χT versus T for the solid 144.The measured data were plotted as open circles, and the fitting curve was obtained using the Bleaney-Bowers equation with g = 2.00.Reproduced from Reference [156-158].Copyright 2011, 2012, 2013, American Chemical Society.Wu et al., in a significant modification in their molecular design, reported the kinetically blocked heptazethrene derivative (141) and octazethrene derivative (142) [157].Compound 141 was found to

Figure 19 .
Figure 19.(a) Absorption spectra of 139 and 140 in chloroform; (b) Plot of χT versus T for the solid 142.Open circles represent the measured data.Fitting curve was drawn using the Bleaney-Bowers equation with g = 2.00; (c) Plot of χT versus T for the solid 144.The measured data were plotted as open circles, and the fitting curve was obtained using the Bleaney-Bowers equation with g = 2.00.Reproduced from Reference [156-158].Copyright 2011, 2012, 2013, American Chemical Society.

Figure 21 .
Figure 21.(a) Molecular structures of 145-150 and (b) plot of χT versus T curves of nPer-CN (n = 2-6) (146-150).The solid lines represent the fitting curves obtained from the Bleaney-Bowers equation; g = 2.00.Reproduced from Reference [159].Copyright 2013, American Chemical Society.The smallest derivative, 145, exhibited well-resolved 1 H and 13 C-NMR signals in CDCl3 at 298 K, suggesting that 145 has a closed-shell quinoidal structure in the ground state.On the other hand, the extended derivatives 146-150 were found to be NMR-silent even at low temperature (−100 °C in CD2Cl2), suggesting that the extended derivatives have a significant biradical character in the ground state.Derivatives 146-150 exhibited broad ESR signals in the solid as well as solution state with g = 2.0031.From the SQUID measurements of 146 and 147 in the powder form at 5-300 K, the singlettriplet energy gap ∆ (i.e., −2J/kB) was found to be 0.342 and 0.107 kcal/mol for 146 and 147, respectively (Figure21b).This further confirmed that both 146 and 147 have a singlet biradical ground state.A very weak coupling between the two radicals could be judged from the small ∆ (−0.064 kcal/mol) for 148.Due to the inevitable error during the data fitting for this borderline molecule, the ground state of 148 was described as a singlet biradical with a very large biradical character based on the variable-temperature ESR experiments and computational studies.Further moving to 149 and 150, the ∆ became −0.56 and −0.88 kcal/mol, respectively, signifying a triplet biradical ground state.

Figure 22 .
Figure 22.(a) Molecular structures of 163 and 164; (b) Plot of χT versus T plot for the solid sample of 163 in SQUID magnetometry studies.Open circles represent the measured data and the red curve denotes the best fit using the Bleaney-Bowers equation with ge = 2.00.Reproduced from Reference [161].Copyright 2016, American Chemical Society.

Figure 22 .
Figure 22.(a) Molecular structures of 163 and 164; (b) Plot of χT versus T plot for the solid sample of 163 in SQUID magnetometry studies.Open circles represent the measured data and the red curve denotes the best fit using the Bleaney-Bowers equation with ge = 2.00.Reproduced from Reference [161].Copyright 2016, American Chemical Society.

Figure 22 .
Figure 22.(a) Molecular structures of 163 and 164; (b) Plot of χT versus T plot for the solid sample of 163 in SQUID magnetometry studies.Open circles represent the measured data and the red curve denotes the best fit using the Bleaney-Bowers equation with g e = 2.00.Reproduced from Reference [161].Copyright 2016, American Chemical Society.

Figure 23 .
Figure 23.Molecular structures of the closed-shell and open-shell canonical forms of 165 and 166.Mes = mesityl groups.Clar's aromatic sextets are depicted in blue hexagons.Reproduced from Reference [162].Copyright 2016, American Chemical Society.

Figure 23 .
Figure 23.Molecular structures of the closed-shell and open-shell canonical forms of 165 and 166.Mes = mesityl groups.Clar's aromatic sextets are depicted in blue hexagons.Reproduced from Reference [162].Copyright 2016, American Chemical Society.

Figure 24 .
Figure 24.(a) The χT-T plot for powder 166.The measured data were plotted as open circles, and the fitting curve was drawn using the Bleaney-Bowers equation with ge = 2.00; (b) Magnetization curves at different temperatures for the powder of 166; (c) χT-T plot at different temperatures for the powder of 165.The measured data were plotted as open circles, and the fitting curve was drawn using the Bleaney-Bowers equation with ge = 2.00; (d) Magnetization curve at 300 K for the powder of 165.Reproduced from Reference [162].Copyright 2016, American Chemical Society.

Figure 24 .
Figure 24.(a) The χT-T plot for powder 166.The measured data were plotted as open circles, and the fitting curve was drawn using the Bleaney-Bowers equation with g e = 2.00; (b) Magnetization curves at different temperatures for the powder of 166; (c) χT-T plot at different temperatures for the powder of 165.The measured data were plotted as open circles, and the fitting curve was drawn using the Bleaney-Bowers equation with g e = 2.00; (d) Magnetization curve at 300 K for the powder of 165.Reproduced from Reference [162].Copyright 2016, American Chemical Society.