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.     

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.     

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.     

A DFT study on the magnetostructural property of ferromagnetic heteroverdazyl diradicals with phenylene coupler

We have theoretically designed five different m-phenylene coupled high-spin bis-heteroverdazyl diradicals and their analogous p-phenylene coupled low-spin positional isomers. The geometry-based aromaticity index, harmonic oscillator model of aromaticity (HOMA) values for both the couplers (local HOMA), and the whole diradicals (global HOMA) have been calculated for all the diradicals. We also qualitatively relate these HOMA values with the intramolecular magnetic exchange coupling constants (J), calculated using a broken symmetry approach within unrestricted density functional theory framework. Structural aromaticity index HOMA of linkage specific benzene rings in our designed diradical systems shows that the aromatic character depends on the planarity of the molecule and it controls the sign and magnitude of J. The predicted J values are explained on the basis of spin polarization maps, average dihedral angles, and magnetic orbitals. The effect of the spin leakage phenomenon on magnetic exchange coupling constant and that on HOMA values of certain phosphaverdazyl systems has been explicitly discussed. In addition, a similar comparison is made between the calculated exchange coupling constants and corresponding HOMA values. The main novelty of this work stands on the consideration of the aromatic behavior by means of the geometrical index HOMA. We also estimate another aromaticity index, nucleus independent chemical shift (NICS) values for the phenylene coupler in each diradical to measure aromaticity and compare its value with that of HOMA. The ground state stabilities of these diradicals have also been compared.     

Experimental and theoretical studies of magnetic exchange in silole-bridged diradicals

Five bis(tert-butylnitroxide) diradicals connected by a silole (7 a-d) or a thiophene (12) ring as a coupler were studied. Compound 12 crystallizes in the orthorhombic space group Pna2(1) with a = 20.752(5), b = 5.826(5), and c = 34.309(5) A. X-ray crystal structure determination, electronic spectroscopy, variable-temperature EPR spectroscopy, SQUID measurements and DFT computations (UB3LYP/6-31+G*) were used to study the molecular conformations and electronic spin coupling in this series of molecules. Whereas compounds 7 b, 7 c, and 7 d are quite stable both in solution and in the solid state, 7 a and 12 undergo a partial electronic rearrangement to both a diamagnetic quinonoidal form and a monoradical species owing to the fact that they correspond to the open form of a pi-conjugated Kekulé structure. In the solid state, magnetic measurements indicate that the diradicals are all antiferromagnetically coupled, as expected from their topology. These interactions are best reproduced by means of a "Bleaney-Bowers" model that gives values of J = -142.0 cm(-1) for 7 a, -1.8 cm(-1) for 7 b, -1.3 cm(-1) for 7 c, -4.2 cm(-1) for 7 d, and -248.0 cm(-1) for 12. The temperature dependence of the EPR half-field transition in frozen dichloromethane solutions is consistent with singlet ground states and thermally accessible triplet states for diradicals 7 b, 7 c, and 7 d with DeltaE(T-S) values of 3.48, 2.09, and 8 cm(-1), respectively. No evidence of a populated triplet state was found for diradicals 7 a and 12. Similarities between the DeltaE(T-S) and J values (DeltaE(T-S) = -2 J) clearly show the intramolecular origin of the observed antiferromagnetic interaction. Analyses of the data with a "Karplus-Conroy"-type equation enabled us to establish that the silole ring, as a whole, allows a more efficient magnetic coupling of the two nitroxide radicals attached to its 2,5-positions than the thiophene ring. This superiority probably originates from the nonaromaticity of the silole which thus permits a better magnetic interaction through it. DFT calculations also support the experimental results, indicating that the magnetic exchange pathway preferentially involves the carbon pi system of the silole.     

Synthesis and Magnetic Properties of (Pyrrolidin‐1‐oxyl)–(Nitronyl Nitroxide)/(Iminonitroxide)‐Dyads

Unlike extensively studied diradicals linked by π‐conjugated systems, only a few studies have investigated weakly coupled diradicals linked by an sp³ carbon atom. Herein, we prepared pyrrolidin‐1‐oxyl–(nitronyl nitroxide)‐dyad 5 and pyrrolidin‐1‐oxyl–iminonitroxide‐dyad 6 . From the observed temperature dependence of the magnetic susceptibility, 5 and 6 were determined to be in singlet ground states with 2J_intra/k_B=?35.2 K and ?13.6 K, respectively. From these results and theoretical calculations of related diradicals, the spin‐polarization model counting the small spin density of the sp³ carbon atom could be used as a spin‐prediction model.     

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.     

Reversible Self-Assembling of Boryl Radical Anions to Their Diradicals with Tunable Singlet Ground States

Two novel boron-centered diradicals based on dimesitylpyridine borane ( 1 ) were synthesized by the self-assembling of the corresponding radical sodium and potassium salts, respectively. The sodium diradical was obtained by re-dissolving the crystals of the radical salt 1Na in toluene, while the potassium diradical was directly obtained by the reduction of 1 with potassium in THF. The diradicals could be converted back to their radical anions in THF solution, forming a reversible process. EPR spectroscopy and SQUID measurements, together with theoretical calculations, show that the diradicals have singlet ground states with excited triplet states. Their singlet–triplet energy gaps are tunable with metals.     

Interaction of the substrate radical and the 5'-deoxyadenosine-5'-methyl group in vitamin B(12) coenzyme-dependent ethanolamine deaminase

The distance and relative orientation of the C5' methyl group of 5'-deoxyadenosine and the substrate radical in vitamin B(12) coenzyme-dependent ethanolamine deaminase from Salmonella typhimurium have been characterized by using X-band two-pulse electron spin-echo envelope modulation (ESEEM) spectroscopy in the disordered solid state. The (S)-2-aminopropanol-generated substrate radical catalytic intermediate was prepared by cryotrapping steady-state mixtures of enzyme in which catalytically exchangeable hydrogen sites in the active site had been labeled by previous turnover on (2)H(4)-ethanolamine. Simulation of the time- and frequency-domain ESEEM requires two types of coupled (2)H. The strongly coupled (2)H has an effective dipole distance (r(eff)) of 2.2 A, and isotropic coupling constant (A(iso)) of -0.35 MHz. The weakly coupled (2)H has r(eff) = 3.8 A and A(iso) = 0 MHz. The best (2)H ESEEM time- and frequency-domain simulations are achieved with a model in which the hyperfine couplings arise from one strongly coupled hydrogen site and two equivalent weakly coupled hydrogen sites located on the C5' methyl group of 5'-deoxyadenosine. This model indicates that the unpaired electron on C1 of the substrate radical and C5' are separated by 3.2 A and are thus at closest contact. The close proximity of C1 and C5' indicates that C5' of the 5'-deoxyadenosyl moiety directly mediates radical migration between cobalt in cobalamin and the substrate/product site over a distance of 5-7 A in the active site of ethanolamine deaminase.     

Spin-orbit coupling in oxygen containing diradicals

Intermediate diradicals which occur in the Paterno-Büchi photocycloaddition and in the Norrish type I photoreactions have been calculated taking into account the spin-orbit coupling (SOC) between the singlet (S) and triplet (T) states. Reaction paths for the photocycloaddition of formaldehyde to ethene and the diradical products of the -cleavage of cyclohexanone have been optimized by the MNDO CI method for a number of different singlet and triplet states. SOC integrals are calculated by an effective one-electron approximation. Intermediate diradicals in the Paterno-Büchi reaction and the SOC effects are also studied ab initio with CAS SCF geometry optimization in a TZV basis set. Both methods predict a large SOC matrix element between the S and T states in the course of the C-C attack, while the SOC integral is two orders of magnitude smaller for the diradical produced in the C-O attack. In the Norrish type I photoreaction the oxygen atom also produces some nonzero contribution to the SOC integral which governs intersystem crossing in a ·C-C· diradical. For the diradicals produced by the -cleavage of cyclohexanone a vibronic interaction is responsible for the SOC mixing between the lowest S and T states. The importance of one-center versus two-center SOC contributions in diradicals is briefly discussed.     

Recent Advances of Stable Phenoxyl Diradicals

Organic radicals with unpaired electrons have shown great semiconducting properties with potential applications in the field of organic photodetectors, organic light-emitting diodes and organic spintronics. A major problem for limiting radicals from laboratory research to practical applications is the relatively low chemical and physical stabilities. Therefore, the right selection of radical core is the key to meaningful scientific research. Phenoxyl radical is one of the few stable radicals with spin distribution properties. Moreover, phenoxyl diradicals provide extract stability due to multiple resonance structures. Due to the long-distance spin distribution, which makes phenoxyl diradicals show interesting electronic and magnetic properties. In this review, we summarize the progress of phenoxyl diradicals in recent years in terms of syntheses, properties and future perspective.     

Stretch Effects Induced by Molecular Strain on Weakening σ-Bonds: Molecular Design of Long-Lived Diradicals (Biradicals)

Stretch effects induced by two types of molecular strain were examined by quantum chemical calculations at the B3LYP/6-31G(d), B3LYP/cc-PVDZ, CCSD/6-31G(d), and CASPT2/ANO-RCC-VTZP level of theory, to design persistent multiradicals such as localized diradicals and oxyallyls. The cooperative molecular strain (Type-1) induced by the spiro[5.5]undecane and bicyclo[2.1.0]cyclopentane structures was found to significantly destabilize in energy the ring-closed compounds of the diradicals, leading to small energy differences between the diradicals and the σ-bonded compounds. Another stretch effect (Type-2) induced by macrocyclic systems was also found to energetically destabilize the corresponding ring-closed structures of the 1,3-diradicals. The computational studies predict that the two types of stretch effects are quite effective in lowering the energy barriers of the bond-breaking reaction of the ring-closed compounds and in generating long-lived localized diradicals and oxyallyl derivatives.     

Reduction of diesters of 1,2-diols. Regioselective C-O bond cleavage of the anionic forms

The electrochemical reduction of benzoate diesters of glycols has been studied in acetonitrile and N,N-dimethylformamide as solvents. The reductions occur in two closely spaced one-electron steps, and it was found that the dianion diradicals decompose by one of two routes, depending on the substituents on the ethylene moiety: cleavage of two benzoates to produce alkene or formation of benzil by way of a postulated cyclic intermediate to produce also the dianion of the diol. These correspond to cleavage of the R-OC(O)Ar bonds and the RO-C(O)Ar bonds, respectively. When the radical formed by the former cleavage is a primary or secondary radical, the reaction is too slow to compete with the latter reaction that produces benzil. However, when that radical is either tertiary or benzylic, the former cleavage reaction is fast and no benzil is detected. The dianions of p-cyano- and p-nitrobenzoate esters are rather stable on the voltammetric time scale. However, the addition of lithium ions results in detectable formation of 4,4'-dicyanobenzil from four different p-cyanobenzoate diesters.     

Spin–orbit coupling in oxygen containing diradicals

Intermediate diradicals which occur in the Paterno–Büchi photocycloaddition and in the Norrish type I photoreactions have been calculated taking into account the spin–orbit coupling (SOC) between the singlet (S) and triplet (T) states. Reaction paths for the photocycloaddition of formaldehyde to ethene and the diradical products of the α-cleavage of cyclohexanone have been optimized by the MNDO CI method for a number of different singlet and triplet states. SOC integrals are calculated by an effective one-electron approximation. Intermediate diradicals in the Paterno–Büchi reaction and the SOC effects are also studied ab initio with CAS SCF geometry optimization in a TZV basis set. Both methods predict a large SOC matrix element between the S and T states in the course of the C–C attack, while the SOC integral is two orders of magnitude smaller for the diradical produced in the C–O attack. In the Norrish type I photoreaction the oxygen atom also produces some nonzero contribution to the SOC integral which governs intersystem crossing in a ·C–C· diradical. For the diradicals produced by the α-cleavage of cyclohexanone a vibronic interaction is responsible for the SOC mixing between the lowest S and T states. The importance of one-center versus two-center SOC contributions in diradicals is briefly discussed.     

Triradical cation of p-phenylenediamine having two nitroxide radical groups: spin alignment mediated by delocalized spin

The cationic state of a p-phenylenediamine (PDA) molecule having two nitroxide radical groups was prepared and characterized using electrochemical, electron spin resonance (ESR) spectroscopic, and absorption spectroscopic methods. The delocalized intervalence state of the p-phenylendiamine (PDA) moiety was detected in the cationic state. From the pulsed ESR measurements, it was confirmed that the delocalized spin induces parallel spin alignment between the localized two nitroxide groups which are magnetically weakly coupled in the neutral state. It was found that the resulting high-spin alignment does not seriously affect the delocalized intervalence state of the PDA radical cation.     

Ab initio STUDIES ON NITROXYL RADICAL MODELS OF ORGANIC FERROMAGNETS

The magnetic coupling between the unpaired spin of nitroxyl group (>N--O) and itsneighbor sp~2 carbon atom has been studied in detail by UHF/6--31G and MP4/6-31G methods.and a number of conjugate nitroxyl diradicals have also been studied by ab initio UHFmethod. A simple topological rule about the ferromagnetic coupling between two nitroxylgroups has been revealed by our calculation. Based on this rule as well as ab initio cal-culation, some nitroxyl radical models of organic ferromagnets (OFM) have been proposed,and the possible way of synthesis has also been suggested. Moreover, the ab initio UHF-crystal orbital (UHF--CO) method has been used to characterize the electronic band struc-ture of one of the models. The ab initio UHF--CO results show that this model is a verypromising candidate of OFM,     

The D parameter of vinyl-substituted 1,3-cyclopentanediyl triplet diradicals as a sensitive tool To determine electronic substituent effects in allylic radicals

The zero-field D parameter of the localized vinyl-substituted 1, 3-cyclopentanediyl triplet diradicals V was determined at 77 K in a 2-methyltetrahydrofuran (2-MTHF) glass matrix. Good linear correlations were obtained with the reported alpha-hyperfine coupling constants (r(2) = 0.991, n = 7) and with the semiempirically calculated (PM3) spin densities (r(2) = 0.989, n = 16) of the corresponding allylic monoradicals A. The observed substituent effects are generally larger and, thus, more accurately measured compared to the previously examined aryl-substituted 1, 3-cyclopentanediyl triplet diradicals P. The vinyl-substituted triplet diradicals reflect accurately the delocalizing propensity of substituents, either through hyperconjugative, mesomeric, or inductive effects. For the methyl group, the small but significant stabilization of adjacent radical sites has been clearly demonstrated. In the case of the halogen set, the small but definitive heavy-atom effect has been determined for the bromo (0. 0010 cm(-)(1)) and the iodo (0.0024 cm(-)(1)) substituents. The organometallic substituents SiMe(3) and SnMe(3) are shown to be weak spin acceptors, while spin delocalization for the sulfur series follows the increasing order MeS > MeSO(2) > MeSO.     

Mechanism of the forbidden [3s,5s]-sigmatropic shift: orbital symmetry influences stepwise mechanisms involving diradical intermediates

The mechanisms of [3s,5s]-sigmatropic shifts of octa-1,3,7-triene and 7-methylenenona-1,3,8-triene have been elaborated using B3LYP and BPW91 density functional theory and CASPT2 methods. These orbital symmetry forbidden rearrangements are stepwise, involving diradical intermediates. A comparison with several [3,3]-sigmatropic shifts of substituted hexadienes and of [5,5]-sigmatropic shifts that are allowed, but nevertheless follow stepwise paths, shows that the activation barrier for the disallowed [3,5] shift is significantly larger than that for the stepwise reactions that are orbital symmetry allowed. Cyclic diradicals that have an aromatic circuit of electrons including the two radical centers and conjugated pi or sigma bonds are stabilized as compared to cyclic diradicals with an antiaromatic circuit of electrons. This applies to the transition states leading to and from the diradicals and influences the activation energies of stepwise sigmatropic shifts. The magnitudes of these effects are small but will have a significant influence on the rates of competing processes. This series of calculations has been used to assess the relative capabilities of the two functionals. We find that BPW91 underestimates the endothermicity of diradical formation and the barrier to diradical formation whereas B3LYP overestimates these quantities.     

A novel nonanuclear CuII carboxylate-bridged cluster aggregate with an S= 7/2 ground spin state

The Cu(II) aggregate in [Cu(9)(cpida)(6)(MeOH)(6)].6(MeOH)[H(3)cpida = 2-(carboxyphenyl)iminodiacetic acid] is made up of two weakly ferromagnetically coupled carboxylate-bridged Cu(4) units that are antiferromagnetically linked through a central Cu(II) to give a Cu(9) core with an S= 7/2 spin ground state.     

A recipe for geometry optimization of diradicalar singlet States from broken-symmetry calculations

The equilibrium geometries of the singlet and triplet states of diradicals may be somewhat different, which may have an influence on their magnetic properties. The single-determinantal methods, such as Hartree-Fock or Kohn-Sham density functional theory, in general rely on broken-symmetry solutions to approach the singlet-state energy and geometry. An approximate spin decontamination is rather easy for the energy of this state but is rarely performed for its geometry optimization. We suggest simple procedures to estimate the optimized geometry and energy of a spin-decontaminated singlet, the accuracies of which are tested on a few organic diradicals. This technique can be generalized to interactions between higher-spin units or to multispin systems.