Toward Molecular Magnets of Organic Origin via Anion−π Interaction

Oct 20, 2016 - A rigorous description of molecular electronic structure within a molecular orbital (MO) framework can be obtained using the full confi...
0 downloads 8 Views 3MB Size
Article pubs.acs.org/JPCA

Toward Molecular Magnets of Organic Origin via Anion−π Interaction Involving m‑Aminyl Diradical: A Theoretical Study Published as part of The Journal of Physical Chemistry virtual special issue “Mark S. Gordon Festschrift”. Debojit Bhattacharya,† Suranjan Shil,‡ Anirban Misra,§ Laimutis Bytautas,*,∥ and Douglas J. Klein† †

Department of Marine Sciences, Texas A&M University at Galveston, Galveston, Texas 77553, United States Center for Atomic Scale Materials Design, Department of Physics, Technical University of Denmark, Kgs. Lyngby 2800, Denmark § Department of Chemistry, University of North Bengal, Darjeeling, PIN. 734013, West Bengal, India ∥ Department of Chemistry, Galveston College, 4015 Avenue Q, Galveston, Texas 77550, United States ‡

S Supporting Information *

ABSTRACT: Here we study a set of novel magnetic organic molecular species with different halide ions (fluoride, chloride, bromide) absorbed ∼2 Å above or below the center of an aromatic π-ring in an m-aminyl diradical. Focus is on the nature of anion−π interaction and its impact on magnetic properties, specifically on magnetic anisotropy and on intramolecular magnetic exchange coupling. In the development of single molecule magnets, magnetic anisotropy is considered to be the most influential factor. A new insight regarding the magnetic anisotropy that determines the barrier height for relaxation of magnetization of m-aminyl diradical-derived anionic complexes is obtained from calculations of the axial zero-field-splitting (ZFS) parameter D. The noncovalent anion−π interaction strongly influences magnetic anisotropy in m-aminyl−halide diradical complexes. In particular, the change of D values from positive (for the m-aminyl diradical, m-aminyl diradical/fluoride, and maminyl diradical/chloride complexes) to negative D-values in m-aminyl diradical complexes containing bromide signals a change from oblate to prolate type of spin-density distribution. Furthermore, the noncovalent halide−π interactions lead to large values of intramolecular magnetic exchange coupling coefficients J exhibiting a ferromagnetic sign. The magnitude of J steadily increases going from anionic complexes containing fluoride to chloride and then to bromide. Relations are sought between the magnetic exchange coupling coefficients J and aromaticity, namely structural HOMA (harmonic oscillator model of aromaticity) and magnetic NICS (nucleus independent chemical shift) aromaticity indices, in particular, the NICSzz(+1) component. Finally, possible numerical checks on the conditions relating to validity of the well-known Yamaguchi’s formula for calculating the exchange coupling coefficient J in diradical systems are discussed.

1. INTRODUCTION Investigations of novel organic molecules having high-spin ground states have intensified in the last few decades yielding a significant number of exciting candidates for magnetic materials.1 The ability of such materials to display tunable macroscopic magnetic properties makes them of high interest for practical applications in spintronics and as memory or sensing materials,2−4 superconductivity,5,6 photomagnetic devices,7−9 and in other areas. Recent literature10−24 contains examples of diradical systems containing π-conjugated fragments as couplers between unpaired spin centers (SC). It is quite common to use (five- or six-membered) aromatic rings as couplers due to their ability to maintain near-planar geometries. Although the use of homoatomic (i.e., carbon-based) aromatic rings is well documented,10−15 the use of heteroatomic aromatic rings as couplers, in addition to carbon atoms also containing nitrogen, boron, sulfur, phosphorus, and other atoms, has been explored as well.16−24 It is of interest to note that, in general, in isomer classes of organic compounds the fraction of diradicals relative to nonradicals (closed-shell © XXXX American Chemical Society

systems) increases dramatically as the number of atoms increases (see, e.g., ref 25 for a case of acyclic conjugated polyenoids). Theoretical models26−28 that use low-level approximations, like Hartree−Fock (HF) theory, in general, are quite useful in describing large molecules at a low cost. The best performance of HF approximation is expected for systems of closed-shell type using restricted HF theory (RHF). For open-shell systems, a better description is often obtained from unrestricted or spinpolarized UHF theory where α and β electrons are allowed to occupy different orbitals in a single UHF determinant. More accurate energies (and densities) can be obtained using unrestricted (spin-polarized) density functional theory (DFT). Although the Kohn−Sham (KS) formulation29 of the density functional theory is conceptually and computationally very similar to the Hartree−Fock Theory, unrestricted KS-DFT Received: September 23, 2016 Revised: October 20, 2016 Published: October 20, 2016 A

DOI: 10.1021/acs.jpca.6b09666 J. Phys. Chem. A XXXX, XXX, XXX−XXX

Article

The Journal of Physical Chemistry A

coupling coefficient J. For a diradical, the high-spin state is triplet and the low-spin state is singlet, and J may be obtained as E(S = 1) − E(S = 0) = −2J. Thus, J > 0 indicates ferromagnetic character whereas J < 0 reflects antiferromagnetic character of the ground state. For π-conjugated planar bipartite systems the sign of J can be predicted using starred/unstarred neighboringsite alternation rule as pioneered by Ovchinnikov48 (see also refs 10, 11, and 33−41). On the one hand, in Kekulé-type molecules antiferromagnetic coupling with a large HOMO− LUMO gap is quite common with the HOMO doubly occupied (antiparallel-spins) and a singlet ground state (see, e.g., ref 49). On the other hand, non-Kekulé molecules, i.e., fully conjugated systems that cannot be represented by neighbor-paired Kekulé structures, can contain several π-nonbonded atoms.50−53 NonKekulé diradicals have two nearly degenerate singly occupied frontier MOs (SOMOs).44−46 According to Borden and Davidson,54 the two SOMOs can be either nondisjoint (spatially coincident on some of the atoms) or disjoint (with all amplitudes on different atoms). Borden and Davidson argue that for nondisjoint nonbonding SOMOs, exchange coupling is of ferromagnetic (Hund-rule) type, whereas for disjoint SOMOs antiferromagnetic exchange coupling often occurs (particularly when the two SOMO amplitudes are on adjacent atoms). Though such singlet ground states with disjoint SOMOs is observed in many systems,44,54 there are exceptions (see, e.g., refs 19 and 22), where a refined rule10,11 predicts antiferromagnetic or ferromagnetic character depending on the (shortest) lengths of coupling pathways between disjoint orbital amplitudes (thereby bringing the predictions for alternant systems into coincidence with Ovchinnikov’s rule48). Thus, the disjoint/nondisjoint SOMO’s rules should be used with care, especially, in nonalternant or heteroatomic systems.55 A remarkable feature observed in π-conjugated diradicals is correlation between the aromaticity of a central fragment (the coupler) and the ground-state multiplicity of the overall structure where two radical centers are separated by this coupler. Indeed, in many studies56−67 it has been observed that the increase in a degree of aromaticity for an aromatic coupler linking two radical centers is associated with the enhancement of high-spin (ferromagnetic) character in diradicals.18,21−23,68−72 In 2015, Ito and Nakano23 on the basis of the heteroacene models involving nitrogen atoms demonstrated that the increase in the aromaticity of the central ring(s) favors the diradical form over the zwitterionic form of resonance structures in their heteroacene systems. This evidence23 suggests that the aromaticity of central rings enhances the diradical character of nitrogen-containing heteroacenes where the unpaired spin densities are localized at nitrogen atoms. In the present study we explore the link between the exchange coupling constant J and aromaticity in terms of several aromaticity indices. We emphasize that aromaticity is, in general, a multidimensional or partially ordered characteristic73−78 (see also a comprehensive discussion by Schleyer and co-workers on this subject in ref 64 and also in recent reviews79−81). In general, aromaticity criteria can be structural, energetic, magnetic, electric, or reactivity-oriented. It is natural to expect that the exchange coupling coefficient J would correlate closely with magneto-electric descriptors, such as nuclear independent chemical shifts (NICS)64,82 or ring currents.83,84 Of the various NICS indices, the NICSzz,66,67,82 a zz-component of the NICS tensor has been found to correlate quite well with the ring currents.83,84 Also, our recent study18 found an excellent correlation between the NICSzz values and

is often more successful in describing open-shell systems, like diradicals. Individual KS-DFT implementations (or methods) vary from each other on the basis of specif ic choices for a functional form corresponding to the exchange−correlation energy.26,27 The frequent success of such DFT functionals notwithstanding, more sophisticated electronic structure theories, like multiconfiguration pair-density functional theory (MC-PDFT), sometimes are required to describe large molecules at a satisfactory level, still seeking to keep the computational cost at a minimum to make computations tractable.30 A rigorous description of molecular electronic structure within a molecular orbital (MO) framework can be obtained using the full configuration interaction (FCI) approach where a wave function of the system is represented by a superposition of all possible Slater determinants (configurations). For example, at the zeroth order the electronic structure of closedshell molecules at their equilibrium geometries is ordinarily well described by a wave function that is represented by a superposition of Slater determinants where all MOs are doubly occupied.31,32 If an ab initio wave function is dominated by a single closed-shell determinant, then the natural orbital occupation numbers are close to either 2 (doubly occupied) or 0 (unoccupied), indicating a singlet ground state. However, in more complex situations, e.g., in molecular dissociations, when a wave function has many important configurations, the NO occupation numbers can deviate drastically from either 2 or 0. A typical illustration is dissociation of the ground-state H2 molecule using a minimal basis set28 when the wave function can be well represented by two doubly occupied configurations. At the equilibrium distance, the natural orbital occupation numbers are 2 (bonding-MO) and 0 (antibonding-MO), whereas at large internuclear distances, the NO occupation numbers are 1 (bonding-MO) and 1 (antibonding-MO) representing two singly occupied hydrogen atomic orbitals as expected. Here, taking plus and minus combinations of bonding and antibonding molecular orbitals yields localized, atomic-like orbitals on each atom. Radical character can also arise on graphene nanoribbons with unpaired electron density localized near edges.33−41 In wave function-based approaches, the degree of polyradical character of such systems can be well represented by the deviation of the natural orbital occupation numbers from 2 (doubly occupied) and 0 (unoccupied), as shown by Lischka and co-workers.39−41 In contrast to closed-shell singlet-state (S = 0) molecules, the unpaired spin alignment42−44 in high-spin molecules can lead to a net magnetic moment (e.g., for a triplet state with S = 1). Here, S represents the total spin quantum number, and the spin multiplicity can be calculated as 2S + 1, yielding, for example, 3 for triplet (S = 1) and 1 for singlet (S = 0) electronic states. Diradicals45,46 having triplet ground states are popular potential candidates for the theoretical design of molecular magnets. Recent experimental research efforts by Rajca and coworkers44,47 have focused on organic molecules possessing high-spin ground states that are significantly lower (at least 2−3 kcal/mol) in energy compared to their excited low-spin states. In the design of high-spin organic molecules, it may be desirable44 to connect radical fragments or radical centers in such a way that the unpaired spins would couple ferromagnetically producing high-spin ground states (S ≥ 1). If the ferromagnetic spin coupling is strong, then the high-spin ground state is well separated from the low-spin excited states as reflected by a large (positive) magnitude of the exchangeB

DOI: 10.1021/acs.jpca.6b09666 J. Phys. Chem. A XXXX, XXX, XXX−XXX

Article

The Journal of Physical Chemistry A the exchange coupling constant J for a number of diradicals with heteroatomic couplers (containing boron and nitrogen atoms). In general, negative values of NICS indicate aromaticity whereas positive values of NICS reflect antiaromatic character. One should keep in mind, however, that in some special cases, for example, involving the d-orbitals of transition metal atoms, the correlation between structural aromaticity index like HOMA and its magnetic counterpart NICS might exhibit opposing trends.85,86 Also, some findings87,88 seem to suggest that negative values of NICS may not reflect aromatic character in some specialized cases, most notably involving transition metals. Nevertheless, with a possible exception of such special cases, the overall usefulness of the NICS index in describing aromaticity and its correlation with the J coefficient for organic systems has been well documented.18,64,82 Furthermore, as has been shown recently in ab initio studies by Trinquier and Malrieu,89 aromaticity can prevent electron pairing and impose polyradical character in certain conjugated alternant hydrocarbons. Indeed, such is the (common) case of the higher polyacenes, which, though having a singlet ground state (according to Ovchinnikov’s rule48), demonstrate90 a strong radicaloid character in central rings with a countervailing aromatic character in noncentral rings. A similar effect accounts36 for unpaired π-electron density near suitable (e.g., “zig-zag”) edges of graphene. A promising avenue to enhance ferromagnetic coupling between unpaired spin centers connected by a π-conjugated coupler is to moderate it via noncovalent ion−π interaction. Ali and Oppeneer reported68 that at equilibrium geometries the interactions between an anion (F−, Cl−, and Br−) and coupler fragment in 1,3-phenylene-based bis(aminoxyl) diradical yield a significant increase in exchange coupling constant J, whereas the interaction between a cation (Li+, Na+, K+) and coupler fragments in 1,3-phenylene-based bis(aminoxyl) diradical yields a slight decrease in J values, making these systems less ferromagnetic. Thus, this study68 suggests that the use of interaction between anion and π-conjugated coupler can be an effective way to enhance ferromagnetic interaction between radical centers. Inspired by the result of Ali and Oppeneer,68 in this work we pursue the anion/coupler noncovalent interaction strategy in exploring the effect of anions not only in enhancing ferromagnetism but also, most importantly, in influencing magnetic anisotropy of organic diradical complexes. Although Ali and Oppeneer68 discussed the influence of noncovalent cation/anion−π interactions on the magnetic exchange phenomenon in the case of 1,3-phenylene-based bis (aminoxyl) diradical, these authors have not explored magnetic anisotropy. Magnetic anisotropy is considered to be the most influential factor in the design of single molecule magnets because it controls the barrier height for relaxation of magnetization.91−106 Due to the critical importance of magnetic anisotropy in development of molecular magnets, in this study we focus on the magnetic anisotropy in terms of the axial zero-field-splitting (ZFS) parameter D. To this end, we choose m-aminyl diradical as our parent system (Figure 1), which has been recently isolated by Rajca and co-workers.12 Starting with the parent structure of the m-aminyl diradical, we explore magnetic properties in its original form as well as under the influence of noncovalent interactions due to the presence of halides F−, Cl−, and Br−. In the m-aminyl diradical complexes the halide ions F−, Cl−, and Br− are located directly above the center of the six-membered ring (A) (Figure 2) of the parent

Figure 1. m-Aminyl diradical synthesized by Rajca and co-workers12 taken as a parent diradical in this study. Ring “A” is the aromatic sextet above (and below) which the halides (F−, Cl−, and Br−) are placed and n-Pr represents the n-propyl group.

Figure 2. Optimized distances (d and d*) of the anions from the center of the π-ring surface. The arrow shows the direction of charge transfer from halide ions to the sextet ring (A). Here, X is the representative of halide ions. (a) depicts anion−π interaction when a halide is placed above the aromatic π-ring. (b) represents a situation when two halides of the same type are placed on both sides of the aromatic π-ring.

diradical (Figure 1) at a distance of about 2 Å. In addition, we also consider a Br−/m-aminyl diradical/Br− complex with the aminyl diradical being sandwiched by two Br− ions. It should be mentioned that, in general, noncovalent interaction between an electron deficient π-system and an anion is of importance in supramolecular chemistry,107 playing a major role not only in solid state but also in solution.108 The common classical charge transfer of anion−π systems offers an exciting avenue to investigate the magnetic properties of various molecular complexes, and X-ray crystallography may be used to demonstrate anion−π interactions.109,110 The judicial design of efficient anion receptors is an important task in investigations of environmental, biological and medical significance.111 Indeed, various ion interactions with neutral closed-shell systems as well as diradicals have been a subject of intense research for many years. For example, pH levels are highly relevant to the stability of certain organomercury compounds112 and can also influence high-to-low-spin ground-state switchover in some organic compounds.113,114 Thus, the current investigation of the noncovalent interactions between various halide anions and the central πring in the m-aminyl diradical12 should be of interest, keeping in mind its potential use as strong organic magnets. C

DOI: 10.1021/acs.jpca.6b09666 J. Phys. Chem. A XXXX, XXX, XXX−XXX

Article

The Journal of Physical Chemistry A

2. COMPUTATIONAL DETAILS The halide anions (F−, Cl−, and Br−) are placed at about 2 Å above (and below in the case of two Br− ions) the ring center of the parent diradical to produce four halide/m-aminyl diradical complexes (Figures 1 and 2 and Table 1). The geometries of

core do not span the same space. Moreover, the a and b orbitals do not span the same space (for the nT and BS cases), which is weaker than requiring that the canonical SCF-like orbtials be the same in these two cases, it being understood that the canonicalness condition of most computations enable mixings of these two near-degenerate SOMOs. There are different possible checks that the Yamaguchi conditions are close to being satisfied. For instance, the expectation ⟨nT|(S⃗)·(S⃗)|nT⟩ = (S2̂ )nT should be close to that for a (pure) triplet, namely, S(S + 1) = 2. And further, ⟨BS|(S⃗)· (S⃗)|BS⟩ = (Ŝ2)BS should be close to that for a state that is half-and-half triplet and singlet, namely 1 1 1(1 + 1) + 2 0(0 + 1) = 1. Another check involves consid2 eration of the fractions wS of spin-S contributions to expectation values that are made by the different wave functions (to be used in the Yamaguchi formula). At least for the nominal triplet |nT⟩, one can bound its excess-spin weight (i.e., its net weights exceeding S = 1) w>1 = w2 + w3 + ..., if one has available the expectation ⟨Ŝ2⟩nT. Indeed,

Table 1. Optimized Distances (Å) from the Center of the Ring “A” (Figure 1)a diradical 2 3 4 5

(parent−F−) (parent−Cl−) (parent−Br−) (Br−−parent−Br−)

UB3LYP

UM06-2X

2.47625 3.38430 3.34714 3.66740/3.93132

2.58004 3.04066 3.17501 3.34411/3.43508

Basis sets are 6-31G(d,p). Here, “parent” structure implies m-aminyl diradical.

a

these five diradicals (parent, parent−F−, parent−Cl−, parent− Br−, Br−−parent−Br−; here “parent” structure implies m-aminyl diradical) are fully optimized at the UB3LYP/6-31G(d,p)115−118 level of theory. Further, the geometries of these five diradicals have been fully optimized using the M06-2X functional.119,120 The magnetic exchange interaction between two spin centers is expressed by the Heisenberg effective Hamiltonian, Ĥ = −2JS1̂ ·S2̂

2

⟨S ̂ ⟩nT =

S≥1

1 2̂ 1 ⟨S − 2⟩nT = 4 4

(1)

(E BS − ET ) 2 2 (⟨S ̂ ⟩T − ⟨S ̂ ⟩BS )

=

(3)

Note that when the ⟨Ŝ2⟩T value is very close to 2, eq 1 leads to ΔES−T ≈ 2J. The Yamaguchi formula of eq 2 is based on assumptions, that there are nominal triplet (nT) and broken-symmetry (BS) states of a form and

∑ wS(nT)· S(S + 1) − 2 4

(6)

1 2̂ 1 ⟨S ⟩nT − 4 2

(7)

Many other (partial) checks are also conceivable, but the S2̂ expectations and w>1 values are ordinarily available in most of quantum chemical packages in use today. The degeneracy of the triplet (S = 1) microstates MS = +1, 0, −1 at the nonrelativistic level of theory can be lifted upon the inclusion of dipolar spin−spin couplings.129 This effect is known as the zero-field-splitting (ZFS). In the DFT framework, the ZFS tensor may be obtained using the Pederson and Khanna method.130 The majority of calculations performed in this study are done by the Gaussian’09W131 quantum chemical package, whereas the axial ZFS parameter D is calculated by the ORCA132 package of programs using UB3LYP functional and 6-31G(d,p) basis set (unrestricted Kohn−Sham DFT framework). The ZFS value arising from the spin−spin (SS) interactions can be estimated through an effective spin Hamiltonian Ĥ ZFS = Σ μ,ν D μν Ŝ μ Ŝ ν , where D μν are the components of the ZFS tensor, Ŝμ is the μth Cartesian component of the total electron spin operator.133 The ZFS129,133−139 axial parameter D is determined from the tensor 1 components as D = Dzz − 2 (Dxx + Dyy) The method of Pederson and Khanna is known to yield a correct sign of the ZFS parameter D.133−137 The zero-field splitting arises from anisotropic electron spin distribution, with the axial parameter D reflecting the energetic separation between MS levels 0 and ±1. If two unpaired spins are localized in parts of a molecule that are separated by a large distance, the magnitude of D is expected to be small (because of its dipolar distance

where EBS and ET are energies of the broken-symmetry singlet (BS) and the unrestricted (nominally) triplet states. Here, ⟨S2̂ ⟩T and ⟨Ŝ2⟩BS are the corresponding average spin-square values. One can estimate the energy separation between the spinprojected singlet and the unrestricted DFT triplet (pure singlet−triplet splitting ΔES−T = ES − ET), by applying Ginsberg’s128 formula:

|nT⟩ = |core × aα × bα⟩

S≥1

But as all the coefficients in this last sum are ≥1, it is seen that this last sum is ≥w>1, for our nT state. That is,

(2)

2

(5)

∑ {S(S + 1) − 2}·wS(nT)

S≥2

w> 1(nT) ≤

ΔES − T = ⟨S ̂ ⟩T × J

S≥2

and

where Ŝ1 and S2̂ are local spin operators on the two component radicals and the singlet−triplet energy gap is given by 2J. The positive and negative sign of J reflects ferromagnetic and antiferromagnetic interactions, respectively.121−125 The J values were calculated using the well-known Yamaguchi formula126,127 J=

∑ S(S + 1)·wS = 2w1 + ∑ wS·S(S + 1)

|BS⟩ = |core × aα × bβ⟩ (4)

where the cores are completely singlet and are the same for each of these two states. As a result, |nT⟩ is a pure triplet and |BS⟩ is an equal mixture of singlet and triplet, and yet further the Yamaguchi formula then gives the splitting between the singlet and triplet. The UHF wave functions, |nT⟩ and |BS⟩ are not precisely of this form in eq 4, because the core is different in the two cases, the core is not a pure singlet, and the spaces spanned by the (spin-free) a and b orbitals are different in the two cases. Indeed, in computational practice the core is in each case a single Slater determinant, and the spin−σ orbitals in the D

DOI: 10.1021/acs.jpca.6b09666 J. Phys. Chem. A XXXX, XXX, XXX−XXX

Article

The Journal of Physical Chemistry A

coefficient (spin-coupling constant) J is one of the most frequently used quantities for describing magnetic properties of molecular systems. On the one hand, when the J value is large and positive, it means that triplet ground state is well below the lowest singlet state. On the other hand, when the value of J is negative, it implies that the singlet state is the ground state. The data for exchange coupling coefficients J at the UB3LYP/631G(d,p) level of theory are displayed in Table 2. Also, in

dependence). The value of D is expected to be positive for an oblate spin distribution where a flattening in one direction is observed, whereas D should be negative for a prolate spindensity distribution. For a positive D value the components MS = ±1 are higher in energy than the MS = 0 level, whereas for negative D the components MS = ±1 are lower than the MS = 0 component. To quantify the degree of aromaticity of the central ring in our halide/m-aminyl diradical complexes, we have used several aromaticity indices, namely, the harmonic oscillator model of aromaticity (HOMA) and different nuclear independent chemical shifts (NICS). We have calculated the harmonic oscillator model of aromaticity as an effective, structure-based aromaticity index. The HOMA index140,141 is defined as α HOMA = 1 − n

Table 2. Optimized UB3LYP Level Absolute Energies (E, hartree), ⟨S2⟩, and Intramolecular Magnetic Exchange Coupling Constants (J, cm−1) Using the 6-31G(d,p) Basis Set for All the Diradicalsa UB3LYP/6-31G(d,p)

n

diradical

2

∑ (R 0 − R i ) i=1

1 (parent)

(8)

where α is an empirical constant for a particular bond in a π system, n is the number of bonds taken in the summation, Ri is the calculated length of a specific bond in a given ring, and R0 is the optimal aromatic bond length taken from literature. Here, we use α = 98.89 Å−2 and R0 = 1.397 Å to calculate the HOMA index of a central ring in m-aminyl diradical-complexes. We have also calculated the NICS index which represents a magnetic index of aromaticity, as proposed by Schleyer and coworkers.63 The aromaticity index NICS represents a negatively signed absolute shielding measured at any desired point of a given system. A more negative value of NICS implies a higher degree of aromaticity. We used gauge independent atomic orbital (GIAO) methodology142 to calculate our NICS data for the ring labeled as “A” (Figure 1). The NICS(0) symbol implies that the probe is placed at the center of the ring, the NICS(+1) symbol means that the probe is 1 Å above the ring center on the same side as the halide ion, and the NICS(−1) symbol implies that the probe is on the opposite side (at −1 Å from the ring center) of the halide ion. The reason behind the calculation of the NICS index at 1 Å above the ring is that the effect of the ring current due to π-electrons is larger and also the effect of σ aromaticity is negligible at this distance. We have also calculated the zz-component of NICS, i.e., NICSzz, namely NICSzz(+1) and NICSzz(−1). The value of NICSzz(+1) represents the zz-component of NICS tensor calculated at 1 Å above the ring center of m-aminyl diradical complex, and the value of NICSzz(−1) is that 1 Å below the ring center. In our earlier studies18 we found that the NICSzz(+1) index has exhibited a strong correlation with the exchange coupling coefficient J as well as with aromaticity index HOMA whereas no such correlations have been observed using the total NICS(+1) index. The advantage of the zz-component NICSzz compared to the total NICS index in describing aromaticity of π-delocalized systems has also been emphasized by Schleyer and co-workers.64 Thus, in this study we shall focus on the NICSzz(+1) component of the NICS index rather than total indices (NICS(0), NICS(−1), and NICS(+1)) to examine a possible correlation of with the exchange coupling constant J in our five diradical complexes.

2 (parent−F−) 3 (parent−Cl−) 4 (parent−Br−) 5 (Br−−parent−Br−) a

E ⟨S2⟩ E ⟨S2⟩ E ⟨S2⟩ E ⟨S2⟩ E ⟨S2⟩

BS

triplet

J, cm−1

−1356.64244 0.978 −1456.49149 1.024 −1816.91163 1.001 −3928.42952 1.013 −6500.14403 0.982

−1356.65351 2.068 −1456.50339 2.059 −1816.92437 2.064 −3928.44261 2.062 −6500.15755 2.066

2229 2524 2630 2739 2727

Here, “parent” structure implies the m-aminyl diradical.

Table 2 are listed the total energy values (in hartree) for brokenspin (BS) state and triplet electronic states as well as their corresponding ⟨S2̂ ⟩ values. It is of interest to note that for the singlet states we have ⟨Ŝ2⟩ = 0 whereas for triplet states one has ⟨Ŝ2⟩ = 2. As one can see from Table 2, for UHF triplet states we have ⟨S2̂ ⟩ values that are very close to 2 whereas for UHF broken-symmetry states the values are very close to 1, i.e., halfway between a singlet and a triplet. This then provides a good check for the validity of the Yamaguchi formula, and in addition, from formula 7 we see that our (nominal) triplet state has a weight fraction of nontriplet contribution of no more than ≈0.8% for our various molecular species. The last column in Table 2 lists the J values (cm−1) for each diradical, as calculated using eq 2. The spin-projected singlet and triplet energy splittings for each of our five diradicals can be estimated using eq 3. Although the UB3LYP/6-31G(d,p) level of theory has been used extensively in describing open-shell systems, it is of interest to explore the sensitivity of our calculations with respect to different DFT functionals and/or basis sets. To this end, we selected UM06-2X functional119,120 that performed quite well in our recent studies of cis−trans conformers of coupled bis-imino nitroxide diradicals.9 The UM06-2X calculations have been performed using 6-31G(d,p) and 6311++G(d,p) basis sets. The relevant data are displayed in Table 3. When comparing the J data for diradicals 1−5 in Tables 2 and 3, one can see that the neutral structure of maminyl diradical has J values that are 10−20% smaller compared to those of anionic complexes. Furthermore, using the 631G(d,p) basis set, UB3LYP level of theory yields J values that exhibit a very slight but increasing trend going from system 2 to 4, whereas the UM06-2X level of theory yields J values that are practically constant and identical for these complexes. Finally, the results for the UM06-2X level of theory and larger basis set, i.e., 6-311++G(d,p), show (Table 3) that all anionic species exhibit J values that are about 700 cm−1 larger (i.e., more

3. RESULTS AND DISCUSSION 3.1. Ground-State Multiplicity and Exchange Coupling Coefficients. Electronic structure theory based on DFT functionals has been used extensively for predicting magnetic properties of organic diradicals. The exchange coupling E

DOI: 10.1021/acs.jpca.6b09666 J. Phys. Chem. A XXXX, XXX, XXX−XXX

Article

The Journal of Physical Chemistry A

Table 3. Total Energies (E, hartree), ⟨S2⟩, and Intramolecular Magnetic Exchange Coupling Constants (J, cm−1) Using the UM06-2X Level of Theory and 6-31G(d,p) and 6-311++G(d,p) Basis Sets for the Diradical Complexesa UM06-2X/6-31G(d,p) diradical 1 (parent) 2 (parent−F−) 3 (parent−Cl−) 4 (parent−Br−) 5 (Br−−parent−Br−) a

E ⟨S2⟩ E ⟨S2⟩ E ⟨S2⟩ E ⟨S2⟩ E ⟨S2⟩

UM06-2X/6-311++G(d,p)

BS

triplet

J

BS

triplet

J

−1356.00862 0.983 −1455.82961 1.007 −1816.25014 0.996 −3927.87599 1.003 −6499.66857 1.022

−1356.02089 2.070 −1455.84420 2.072 −1816.26469 2.069 −3927.89074 2.068 −6499.68551 2.071

2477

−1356.30877 0.977 −1456.19085 0.999 −1816.59888 0.987 −3930.58122 0.993 −6504.77994 0.980

−1356.32010 2.073 −1456.20510 2.076 −1816.61322 2.072 −3930.59573 2.072 −6504.79521 2.074

2269

3007 2976 3040 3544

2904 2901 2951 3063

Here, “parent” structure implies the m-aminyl diradical.

Figure 3. Spin density plots of m-aminyl diradical systems investigated in the present work at the UB3LYP/6-31G(d,p) level of theory. The hydrogen atoms are removed due to picture clarity. The mauve and the sea green colors represent (α) and (β) spin densities, respectively.

ferromagnetic) compared to the original m-aminyl diradical. Thus, these findings suggest that in the quest for diradical systems having large singlet−triplet gap and high stability at room temperature, the anionic m-aminyl diradical complexes should be of high practical interest, especially keeping in mind that these diradicals can be excellent candidates to serve as building blocks for high-spin polyradicals12,44 in the development of robust magnetic polymers.143−146

To summarize, all the results obtained using UB3LYP and UM06-2X functionals and basis sets described above clearly demonstrate that the halides F−, Cl−, and Br− significantly enhance the ferromagnetic character of the m-aminyl diradical via noncovalent interaction. For the sake of simplicity, in the present study we primarily focus on results obtained using the UB3LYP/6-31G(d,p) level of theory to discuss magnetic properties of our diradical systems. F

DOI: 10.1021/acs.jpca.6b09666 J. Phys. Chem. A XXXX, XXX, XXX−XXX

Article

The Journal of Physical Chemistry A 3.2. Spin-Density Distribution. Tables 2 and 3 display the data for the broken-symmetry state (BS) and nominal triplet for each of the five m-aminyl diradical-derived systems. One can estimate the energy separation between the spin-projected singlet and the unrestricted DFT triplet (representing pure singlet−triplet splitting ΔES−T = ES − ET), using the formula in eq 3 after the spin-coupling coefficient J is obtained from Yamaguchi’s formula in eq 2. Because J coefficients are positive for all five complexes, the positive value of ΔES−T indicates that triplets lie below singlet states, thus yielding a triplet as the ground state for all five systems. A further indication of the dominance of a triplet ground state in these systems can be seen from spin-density plots displayed in Figure 3. The spindensity plots in Figure 3 are determined as differences between (α) and (β) atomic spin density contributions. From the graphical representation of spin-density plots shown in Figure 3 one can observe an alternating spin-density distribution pattern in all five diradical complexes, a pattern characteristic of spinpolarization and triplet ground states.10,11,33,34,147 In terms of frontier molecular orbitals, different singly occupied molecular orbitals (SOMOs) may arise if the energetic separation between two SOMOs, having orbital energies ESMO2 and ESMO1 is sufficiently small. In particular, on the basis of arguments put forward by Hund that lead to the establishment of the well-known Hund’s rule,148 when the ΔESMO‑SMO = (ESMO2 − ESMO1) value is sufficiently small (i.e., when there is near-degeneracy), the ground state is expected to be a triplet. Note that all orbitals discussed in this study represent Kohn− Sham type29 molecular orbitals. SOMO1 and SOMO2, as well as LUMOs for all five diradical systems are shown in Figure 4. In 1968, Hoffmann, Zeiss, and Van Dine149 using extended Hückel calculations to provide an empirical criterion stating that if ΔESMO‑SMO < 1.5 eV, then two nonbonding electrons will favor occupying different molecular orbitals with parallel-spin orientations resulting in a triplet ground state. More recently, Constantinides and co-workers,147 based on calculations using B3LYP level of theory for a set of about 14 polyheteroacenes concluded that systems with ΔESMO‑SMO > 1.3 eV are clear singlets. Molecules with small ΔESMO‑SMO values have a tendency to have large positive ΔES−T = ES − ET (singlet− triplet splittings, here ES = singlet energy and ET = triplet energy) and vice versa, molecules with large ΔESOMO−SOMO values have small or negative ΔES−T’s. In the present study we find that SOMO−SOMO energy gaps are very small and for any of our diradical structures do not exceed 0.8 eV (SOMO1, SOMO2, and LUMO total orbital energies are listed in Figure 4) and, according to the criteria presented above, are consistent with the fact that these diradical structures have triplet ground states. 3.3. Aromaticity and Magnetism. Both aromaticity and magnetism are of fundamental significance in discussions of physical and chemical properties of molecular systems. Unlike magnetism, however, aromaticity is best described as a partially ordered quantity having a multidimensional character.73−78 It has been observed18,21−23,68−72 that there is a close correlation between the aromaticity of coupler fragments that link the unpaired-spin centers and exchange coupling coefficients J in diradicals. The majority of the evidence suggests that an increase in aromaticity of such couplers favors a ferromagnetic trend. Although various aromaticity descriptors exist in the literature, our recent study of borazine couplers18 suggests that NICSzz and HOMA indices should be well suited to describe a correlation between magnetic properties in terms of exchange

Figure 4. Molecular orbitals (SOMO1, SOMO2, and LUMO) of the diradicals (iso-value 0.02). Below each orbital its orbital energy is listed in electronvolts. The results correspond to the UB3LYP level of theory using the 6-31G(d,p) basis set. The maroon and the green color represent different signs (phases) of the orbital coefficients.

coupling constants J and aromaticity of the π-conjugated couplers in our diradical complexes. The most relevant among the NICS indices was found to be the NICSzz(+1) component. The relevant data are listed in Tables 2, 3, and 4, and displayed graphically in Figure 5, where panel A represents the (UB3LYP/6-31G(d,p)) data and panel B represents the (UM06-2X/6-31G(d,p)) results. It is clear that the (UB3LYP/6-31G(d,p)) data reflect a higher degree of correlation between J and NICSzz(+1) than the (UM06-2X/ 6-31G(d,p)) results. However, on the average, both DFT functionals exhibit the same trend, namely, that as NICSzz(+1) values become more negative, indicating an increasing degree of aromaticity of a coupler, the J values become larger, reflecting the increase in ferromagnetic character of a diradical complex. G

DOI: 10.1021/acs.jpca.6b09666 J. Phys. Chem. A XXXX, XXX, XXX−XXX

Article

The Journal of Physical Chemistry A

Table 4. Various Nuclear Independent Chemical Shift (NICS) Values (ppm) for the Different Prototype Systems Considered in This Study Evaluated at Optimized Geometriesb at the UB3LYP/6-31G(d,p) Level of Theorya UB3LYP diradical 1 2 3 4 5

(parent) (parent−F−) (parent−Cl−) (parent−Br−) (Br−−parent−Br−)

UM06-2X

NICS(0)

NICS(+1)/ NICS(−1)

NICSzz(+1)/ NICSzz(−1)

NICS(0)

NICS(+1)/NICS(−1)

NICSzz(+1)/ NICSzz(−1)

−0.0817 0.6721 0.1489 0.2647 0.0236

−3.4775 −2.7190/−2.4734 −2.9489/−3.2563 −2.7853/−3.1840 −2.7549/−3.1779

−2.5392 −5.3613/−1.6173 −7.6200/−2.7172 −9.5351/−2.8847 −9.7402/−7.8475

−0.7747 1.3703 1.2237 1.4416 1.3660

−3.7064 −3.2721/−2.1103 −2.9965/−2.9069 −2.7370/−3.0000 −2.4849/−2.7454

−3.0239 −8.5134/−2.6352 −11.9708/−1.6903 −12.8532/−3.2519 −12.4035/−10.5995

Here, “parent” structure implies m-aminyl diradical. bThe NICS(0), NICS(+1), NICS(−1), NICSzz(+1) and NICSzz(−1) are the representative of the NICS values at the center, 1 Å below and 1 Å above of the sextet ring ‘A’ of the parent aminyl diradical. Note that 1 Å above of the ring “A” means that the dummy atom is placed in between the sextet “A” and the halide anion. NICS(+1) and NICSzz (+1) are the values of NICS when a probe is placed 1 Å above the plane of the ring toward the anion, but NICS(−1) and NICSzz(−1) are the values of NICS when a probe is placed 1 Å below the plane of the ring in opposite to the anion.

a

Figure 5. (A) Plot of exchange-coupling constants J (cm−1) versus NICSzz(+1) (units of ppm) at the UB3LYP/6-31G(d,p) level of theory. (B) Plot of exchange-coupling constants J (cm−1) versus NICSzz(+1) (units of ppm) at the UM06-2X/6-31G(d,p) level of theory. Here, NICSzz(+1) represents NICSzz values when a probe is placed 1 Å above the center of the hexagon ring “A” toward the anion (see text).

correlation between HOMA and the J constant (6-31G(d,p) basis set). Namely, as HOMA values come closer to 1 (highest degree of aromaticity), the values of J become larger, implying greater ferromagnetic character. Of course, as discussed in our earlier study,18 it is somewhat expected that the correlation between the exchange coupling coefficient J and various aromaticity indices, i.e., HOMA and NICSzz(+1), may not be perfect for complex heteroatomic systems, such as considered in this study. 3.4. Zero-Field Splitting. It is well recognized91−97 that magnetic anisotropy is the most influential factor in design of single-molecule magnets because it controls the barrier height for (slow) relaxation of magnetization. The sign and magnitude of the axial zero-field splitting parameter (D) plays a significant role101 in this context. As mentioned earlier, for the m-aminyl diradical Rajca and co-workers12 reported the experimental ZFS parameter D of 1.8 × 10−2 cm−1. The calculation of the axial parameter D using B3LYP/EPR-II level of theory yields a slightly larger value (relative to the experiment) of 3.42 × 10−2 cm−1. Despite this small deviation, the positive sign for D (>0) clearly indicates an oblate-like shape for the overall spin density (i.e., f lattening in one direction). Rajca et al.12 also noted that, by contrast, a similar system, m-phenylene−nitroxide diradical, exhibits a prolate-like overall spin density (i.e., elongated shape in

This observation indicates the same trend that has been observed in our earlier study18 for borazine couplers. The correlation between another aromaticity index, HOMA, and the exchange coupling constant J also indicates a reasonable correlation between these quantities for the set of our five diradical systems (see the data in Table 5 and Figure 6). The data obtained by both DFT functionals, i.e., UB3LYP and UM06-2X, exhibit a similar behavior in describing the Table 5. Harmonic Oscillator Model of Aromaticity (HOMA) Values Using UB3LYP/6-31G(d,p) and UM062X/6-31G(d,p) Levels of Theorya HOMA diradical 1 2 3 4 5

(parent) (parent−F−) (parent−Cl−) (parent−Br−) (Br−−parent−Br−)

UB3LYP/6-31G(d,p)

UM06-2X/6-31G(d,p)

0.8566 0.8750 0.8806 0.8790 0.9158

0.8767 0.9070 0.8897 0.8897 0.9145

a

Results correspond to the six membered ring above which the anions are placed. HOMA model constants are α = 98.89 Å−2 and R0 = 1.397 Å (see text). Here, “parent” structure implies m-aminyl diradical. H

DOI: 10.1021/acs.jpca.6b09666 J. Phys. Chem. A XXXX, XXX, XXX−XXX

Article

The Journal of Physical Chemistry A

Figure 6. (A) Plot of exchange-coupling constants J (cm−1) versus HOMA aromaticity indices at UB3LYP level of theory. (B) Plot of exchangecoupling constants J (cm−1) versus HOMA aromaticity indices at UM06-2X level of theory. All calculations use 6-31G(d,p) basis sets.

energy than the MS = 0 level.100,102 In general, especially for very high-spin S states, a negative value of the ZFS-parameter D assures that the maximum spin projection (MS = ±S) along the z-axis will be the lowest in energy, whereas smaller magnitude MS projections will be higher in energy. As a consequence, a classical rotation of the projection from MS = +S to MS = −S will require an excitation over a significant energy barrier.101−103 Thus, it is generally accepted that, in the absence of quantum tunneling, molecules having D < 0 values are of critical importance in the design of single-molecule or singlechain magnets because the ground state for such systems exhibits magnetic bistability.101 In particular, an energy barrier (also known as the magnetic anisotropy energy barrier) that separates MS = +S and MS = −S sublevels is equal102 to |S2·D| (for integer S). Thus, for the triplet state (S = 1), the energy barrier between sublevels MS = +1 and MS = −1 is equal to |2D|. For the parent m-aminyl diradical our calculated ZFS axial parameter D of 3.42 × 10−2 cm−1 practically coincides with the value reported by Rajca et al.12 The calculated values for the axial ZFS parameter D in m-aminyl diradical complexes 2 and 3 containing F− and Cl− ions, respectively, exhibit positive values of a similar magnitude of about 3.2 × 10−2 cm−1. On the contrary, for m-aminyl diradical/Br− and Br−/m-aminyl diradical/Br− complexes, the calculations yield negative values of −0.60 × 10−2 and −0.93 × 10−2 cm−1, respectively (Table 6). The change of axial parameter D values going from positive (for structures 1, 2, and 3) to negative in Br−-containing maminyl diradical complexes 4 and 5, indicates a transition from oblate-like overall spin-density to prolate-like overall spin-density distribution (Figure 3). Such oblate-to-prolate spin transition (with reorientation of the ZFS tensor) has been often discussed in the context of building conjugated diradical-based systems from monomers to polymer chains.101−106 Furthermore, the magnitude of the axial parameter D for Br−-containing maminyl diradical complexes 4 and 5 is signif icantly smaller compared to the values for structures representing parent maminyl diradical, F−-containing m-aminyl diradical, and Cl−containing m-aminyl diradical complexes, respectively. One can reason that the spin−spin interaction between unpaired spins makes use of bromide’s more diffused valence atomic orbitals

one direction) with a negative value of parameter D (D < 0). In fact, there are many examples of molecular systems exhibiting negative D values (see, e.g., refs 98−103). In this study, we perform ZFS calculations using the UB3LYP/6-31G(d,p) level of theory, reproducing the same value of axial parameter D for the parent structure, i.e., maminyl diradical, just as reported by Rajca and co-workers12 (Table 6). In addition, we calculated axial parameter D values Table 6. Static Zero-Field Splitting (ZFS) Axial Parameter D Values (cm−1) Obtained at the UB3LYP/6-31G(d,p) Level of Theorya diradical 1 2 3 4 5 a

(parent) (parent−F−) (parent−Cl−) (parent−Br−) (Br−−parent−Br−)

D (cm−1) 3.42 3.11 3.20 −0.60 −0.93

× × × × ×

10−2 10−2 10−2 10−2 10−2

Here, “parent” structure implies m-aminyl diradical.

for anion-containing complexes 2, 3, 4, and 5 to investigate the extent of triplet-state delocalization and shape of overall spin density in these, as yet unexplored, systems. The ZFS interaction being approximated by a point-dipole model is useful in providing an estimate of the average interelectron distance between unpaired spins in the triplet state and, consequently, the extent of triplet-state delocalization. On the one hand, the magnitude of the axial ZFS parameter D is indicative of how strongly the unpaired spins interact in a triplet state. On the other hand, the sign of D indicates the shape of the overall spin density (oblate shape when D > 0 and prolate shape when D < 0). For the triplet (S = 1), the three magnetic sublevels MS = −1, 0, +1, split into 0 and ±1 levels due to the inclusion of spin−spin interaction between unpaired spins. If the zero-field-splitting axial parameter D has a negative sign, then the degeneracy in magnetic sublevels MS is removed by placing higher magnitude MS levels lower in energy, whereas the situation is reversed for the positive values of D. Thus, if D > 0, then MS = ±1 sublevels will be higher in energy compared to MS = 0; and if D < 0, the MS = ±1 sublevels will be lower in I

DOI: 10.1021/acs.jpca.6b09666 J. Phys. Chem. A XXXX, XXX, XXX−XXX

Article

The Journal of Physical Chemistry A

design of single-molecule or single-chain organic magnets giving rise to a magnetically bistable ground state.101 The Br−/ m-aminyl diradical/Br− complex exhibits the negative D value of the largest magnitude. This result would suggest that the development of organic single molecule magnets can make use of atoms/ions having larger (more diffused) atomic orbitals in their electronic valence shells, e.g., Br− ions. An even more promising candidate to yield negative D values of large magnitude is the m-aminyl diradical−halide complex with two iodide (I− ions) sandwiching the m-aminyl “parent” structure. Furthermore, all five m-aminyl diradical-derived complexes display large and positive exchange coupling coefficient J values. An increase of J values is observed going from the neutral maminyl diradical (“parent”) structure to anionic m-aminyl diradical complexes involving the halides F−, Cl−, and Br−. A further reflection of high-spin (i.e., triplet) ground states of our diradical complexes was found in the alternating total spindensity distribution pattern, which is characteristic of spinpolarization and triplet ground states. A reasonable correlation has been observed between the exchange coupling coefficient J and aromaticity indices NICSzz(+1) and HOMA, which is consistent with our earlier study.18 Although studies of magnetic properties of organic aromatic molecules have a very long history,150 it is very satisfying to observe the manifestation of a relationship between magnetism and aromaticity in such diverse complexes as investigated in this study. This study reveals a somewhat different trend for J and D quantities when we compare our five diradicals derived from the m-aminyl “parent” structure. In particular, although the magnitude of the exchange coefficient J for the neutral m-aminyl diradical is much smaller compared to the anionic m-aminyl− halide complexes, by contrast, the calculated axial ZFS parameter D values display a very different trend. For the neutral m-aminyl diradical as well as the anionic m-aminyl diradical−halide complexes involving F− and Cl− ions, the ZFS parameter D values are almost identical and positive. However, the anionic m-aminyl diradical−halide complexes involving Br− ions exhibit negative and smaller D values. The extent of tripletstate spin-density distribution pattern changes considerably going from neutral m-aminyl diradical as well as m-aminyl diradical−halide complexes with F− and Cl− ions, to m-aminyl diradical/Br− and Br−/m-aminyl diradical/Br− complexes, highlighting the oblate-to-prolate spin-density transition with a reorientation of the ZFS tensor. Overall, the results presented in this study are expected to be useful in the rational design of high-spin organic molecular magnets.

reflecting the increase in an average distance between unpaired electron spins. As a result of the increase of an average distance between unpaired spins, the strength of the spin−spin interaction is reduced, which, in turn, narrows the splitting between MS = 0 and MS = ±1 components of the triplet state, as is reflected by the diminished absolute value of D in complexes containing bromide (Table 6 and Figure 3). The practical implications for design of single molecule magnets that involve the use of halide anions favor the use of atoms/ions having more diffused atomic orbitals in their electronic valence shells, namely, bromide ions. The Br−/m-aminyl diradical/Br− complex where the Br− ions are sandwiching the m-aminyl “parent” structure exhibits a negative D value of the largest magnitude. Because molecules having a negative D value of large magnitude are essential101 in the design of single-molecule or single-chain organic magnets, a very promising candidate is the m-aminyl diradical−halide complex with two iodide (I− ions) sandwiching the m-aminyl “parent” structure. Thus, it is beneficial to pursue a new strategy of building single molecule organic magnets having negative (and large magnitude) D values, and therefore large barrier heights for (slow) relaxation of magnetization by using heavier halide ions, for example, Br− or I− ions. Such external ionic species could be stabilized by inserting them in molecular crystals. It is worthwhile to discuss the overall spin-density distribution displayed in Figure 3. In particular, although large spin-density contributions can be seen on anionic species of Cl− and Br− in 3, 4, and 5, the spin-density contribution on F− in the complex 2 is rather small. The axial ZFS parameter D data (positive versus negative values) in Table 6 seem to reflect the fact that the spin−spin interactions are more similar in structures 1, 2, and 3 than in structures 4 and 5. In Figure 3 we observe that in m-aminyl diradical complexes with F−, Cl−, and Br−, the unpaired spin density on anion increases dramatically going from F− to Cl− and to Br−. We also observe that the unpaired spin-density type on halides is of the same type as on unpaired spin centers localized on the nitrogen atoms represented by blue-colored spheres in Figure 3. These results suggest that the presence of an external ion near aromatic coupler significantly affects the strength of the magnetic exchange interaction as can be seen in a steady increase of the magnitude of the magnetic exchange interaction constant J as we go from 1 to 5 structures, as shown in Table 2 (UB3LYP/6-31G(d,p)). This sequence correlates well with the steady increase of the unpaired-spin density on halide ions going from structures 1 to 5, as seen in Figure 3 (UB3LYP/6-31G(d,p)).

4. CONCLUSIONS In this investigation we explored the nature of noncovalent anion−π interaction and its impact on magnetic properties, specifically on magnetic anisotropy and on intramolecular magnetic exchange coupling. We find that noncovalent anion−π interaction strongly influences magnetic anisotropy in m-aminyl−halide diradical complexes, which is a significant new result, keeping in mind that for single molecule magnets magnetic anisotropy is the most influential factor. This study demonstrates that the sign and magnitude of the axial zerofield-splitting (ZFS) parameter D depends on type and geometrical arrangement of halide ions positioned above (or above and below) the central ring of the m-aminyl diradical. Only the anionic m-aminyl diradical−halide complexes involving Br− ions exhibit negative D values. This fact is important because molecules having D < 0 are essential in the



ASSOCIATED CONTENT

* Supporting Information S

The Supporting Information is available free of charge on the ACS Publications website at DOI: 10.1021/acs.jpca.6b09666. Mulliken charges in the halide ions (in electron charge units, |e|) for selected DFT functionals, optimized XYZ coordinates and the respective optimized geometries of all the magnetic species in the gas phase, Mulliken charge distribution structures, complete reference to the Gaussian’09 W package of quantum chemistry programs (PDF) J

DOI: 10.1021/acs.jpca.6b09666 J. Phys. Chem. A XXXX, XXX, XXX−XXX

Article

The Journal of Physical Chemistry A



(17) Ko, K. C.; Cho, D.; Lee, J. Y. Scaling Approach for Intramolecular Magnetic Coupling Constants of Organic Diradicals. J. Phys. Chem. A 2013, 117, 3561−3568. (18) Bhattacharya, D.; Shil, S.; Misra, A.; Bytautas, L.; Klein, D. J. Borazine: Spin Blocker or Not ? Phys. Chem. Chem. Phys. 2015, 17, 14223−14237. (19) Shil, S.; Roy, M.; Misra, A. Role of the Coupler to Design Organic Magnetic Molecules: LUMO Plays an Important Role in Magnetic Exchange. RSC Adv. 2015, 5, 105574−105582. (20) Cho, D.; Ko, K. C.; Ikabata, Y.; Wakayama, K.; Yoshikawa, T.; Nakai, H.; Lee, J. Y. Effect of Hartree-Fock Exact Exchange on Intramolecular Magnetic Coupling Constants of Organic Radicals. J. Chem. Phys. 2015, 142, 024318. (21) Zhang, F.; Song, X.; Bu, Y. Redox-Modulated Magnetic Transformations between Ferro- and Antiferromagnetism in Organic Systems: Rational Design of Magnetic Organic Molecular Switches. J. Phys. Chem. C 2015, 119, 27930−27937. (22) Ali, Md. E.; Datta, S. N. Broken-Symmetry Density Functional Theory Investigation on Bis-Nitronyl Nitroxide Diradicals: Influence of Length and Aromaticity of Couplers. J. Phys. Chem. A 2006, 110, 2776−2784. (23) Ito, S.; Nakano, M. Theoretical Molecular Design of Heteroacenes for Singlet Fission: Tuning the Diradical Character by Modifying π-Conjugation Length and Aromaticity. J. Phys. Chem. C 2015, 119, 148−157. (24) Cho, D.; Ko, K. C.; Lee, J. Y. Organic Magnetic Diradicals (Radical-Coupler-Radical): Standardization of Couplers for Strong Ferromagnetism. J. Phys. Chem. A 2014, 118, 5112−5121. (25) Bytautas, L.; Klein, D. J. Isomer Combinatorics for Acyclic Conjugated Polyenes: Enumeration and Beyond. Theor. Chem. Acc. 1999, 101, 371−387. (26) Jensen, F. Introduction to Computational Chemistry; John Wiley & Sons, Ltd.: Chichester, England, 2007. (27) Cramer, C. J. Essentials of Computational Chemistry: Theories and Models; John Wiley & Sons Ltd.: Chichester, England, 2004. (28) Szabo, A.; Ostlund, N. S. Modern Quantum Chemistry: Introduction to Advanced Electronic Structure Theory; Mc-Graw Hill, Inc.: New York, 1989. (29) Kohn, W.; Sham, L. J. Self-Consistent Equations Including Exchange and Correlation Effects. Phys. Rev. 1965, 140, A1133− A1138. (30) Odoh, S. O.; Li Manni, G.; Carlson, R. K.; Truhlar, D. G.; Gagliardi, L. Separated-Pair Approximation and Separated-Pair PairDensity Functional Theory. Chem. Sci. 2016, 7, 2399−2413. (31) Bytautas, L.; Henderson, T. M.; Jiménez-Hoyos, C. A.; Ellis, J. K.; Scuseria, G. E. Seniority and Orbital Symmetry as Tools for Establishing a Full Configuration Interaction Hierarchy. J. Chem. Phys. 2011, 135, 044119. (32) Bytautas, L.; Scuseria, G. E.; Ruedenberg, K. Seniority Number Description of Potential Energy Surfaces: Symmetric Dissociation of Water, N2, C2 and Be2. J. Chem. Phys. 2015, 143, 094105. (33) Klein, D. J.; Nelin, C.; Alexander, S.; Matsen, F. A. High-Spin Hydrocarbons. J. Chem. Phys. 1982, 77, 3101−3108. (34) Klein, D. J.; Alexander, S. A. In Organic Polyradicals, High-Spin Hydrocarbons, and Organic Ferromagnets, Chemical Applications of Topology and Graph Theory; King, R. B., Ed.; Elsevier: Amsterdam, 1987; pp 404−419. (35) Lieb, E. H. Two Theorems on the Hubbard Model. Phys. Rev. Lett. 1989, 62, 1201−1204. (36) Klein, D. J.; Bytautas, L. Graphitic Edges and Unpaired πElectron Spins. J. Phys. Chem. A 1999, 103, 5196−5210. (37) Ivanciuc, O.; Klein, D. J.; Bytautas, L. Unpaired π-Spin Density in Defected Graphite. Carbon 2002, 40, 2063−2083. (38) Ivanciuc, O.; Bytautas, L.; Klein, D. J. Mean-Field ResonatingValence-Bond Theory for Unpaired π-Electrons in Benzenoid Carbon Species. J. Chem. Phys. 2002, 116, 4735−4748. (39) Plasser, F.; Pašalić, H.; Gerzabek, M. H.; Libisch, F.; Reiter, R.; Burgdörfer, J.; Müller, T.; Shepard, R.; Lischka, H. The Multiradical

AUTHOR INFORMATION

Corresponding Author

*L. Bytautas. E-mail: [email protected]. Phone number: 1-409944-1273. Notes

The authors declare no competing financial interest.



ACKNOWLEDGMENTS The authors thank Md. Ehesan Ali for useful discussions and interest in this work. D.B. and D.J.K. acknowledge the support (via grant BD-0894) of the Welch Foundation, of Houston, Texas, USA. A.M. is grateful to DST, India, for financial support.



REFERENCES

(1) Miller, J. S. Magnetically Ordered Molecule-Based Materials. Chem. Soc. Rev. 2011, 40, 3266−3296. (2) Wolf, S. A.; Awschalom, D. D.; Buhraman, R. A.; Daughton, J. M.; von Molnár, S.; Roukes, M. L.; Chtchelkanova, A. Y.; Treger, D. M. Spintronics: A Spin-Based Electronics Vision for the Future. Science 2001, 294, 1488−1495. (3) Sugawara, T.; Matsushita, M. M. Spintronics in Organic πElectronic Systems. J. Mater. Chem. 2009, 19, 1738−1753. (4) Ratera, I.; Veciana, J. Playing with Organic Radicals as Building Blocks for Functional Molecular Materials. Chem. Soc. Rev. 2012, 41, 303−349. (5) Kobayashi, H.; Kobayashi, A.; Cassoux, P. BETS as a Source of Molecular Magnetic Superconductors (BETS = bis(ethylenedithio)tetraselenafulvalene). Chem. Soc. Rev. 2000, 29, 325−333. (6) Uji, S.; Shingawa, H.; Terashima, T.; Yakabe, T.; Terai, Y.; Tokumoto, M.; Kobayashi, A.; Tanaka, H.; Kobayashi, H. MagneticField-Induced Superconductivity in a Two-Dimensional Organic Conductor. Nature 2001, 410, 908−910. (7) Matsuda, K.; Matsuo, M.; Irie, M. Photoswitching of Intramolecular Magnetic Interaction Using Diarylethene with Oligothiophene π-Conjugated Chain. J. Org. Chem. 2001, 66, 8799−8803. (8) Tanifuji, N.; Irie, M.; Matsuda, K. New Photoswitching Unit for Magnetic Interaction: Diarylethene with 2,5-Bis(arylethynyl)-3-thienyl Group. J. Am. Chem. Soc. 2005, 127, 13344−13353. (9) Bhattacharya, D.; Shil, S.; Misra, A.; Bytautas, L.; Klein, D. J. Photomagnetic and Nonlinear Optical Properties in Cis-Trans Green Fluoroprotein Chromophore Coupled Bis-Imino Nitroxide Diradicals. Int. J. Quantum Chem. 2015, 115, 1561−1572. (10) Datta, S. N.; Trindle, C. D.; Illas, F. Theoretical and Computational Aspects of Magnetic Organic Molecules; Imperial College Press: London, 2014. (11) Trindle, C.; Datta, S. N. Molecular Orbital Studies on the Spin States of Nitroxide Species: Bis- and Trisnitroxymetaphenylene, 1,1Bisnitroxyphenylethylene, and 4,6-Dimethoxy-1,3-Dialkylnitroxy-Benzenes. Int. J. Quantum Chem. 1996, 57, 781−799. (12) Rajca, A.; Olankitwanit, A.; Rajca, S. Triplet Ground State Derivative of Aza-m-xylylene Diradical with Large Singlet-Triplet Energy Gap. J. Am. Chem. Soc. 2011, 133, 4750−4753. (13) Boratyński, P. J.; Pink, M.; Rajca, S.; Rajca, A. Isolation of the Triplet Ground State Aminyl Diradical. Angew. Chem., Int. Ed. 2010, 49, 5459−5462. (14) Olankitwanit, A.; Pink, M.; Rajca, S.; Rajca, A. Synthesis of Azam-Xylylene Diradicals with Large Singlet-Triplet Energy Gap and Statistical Analyses of Their Spectra. J. Am. Chem. Soc. 2014, 136, 14277−14288. (15) Ko, K. C.; Park, Y. G.; Cho, D.; Lee, J. Y. Simple but Useful Scheme Toward Understanding of Intramolecular Magnetic Interactions: Benzene-Bridged Oxoverdazyl Diradicals. J. Phys. Chem. A 2014, 118, 9596−9606. (16) Ali, Md. E.; Staemmler, V.; Illas, F.; Oppeneer, P. M. Designing the Redox-Driven Switching of Ferro- to Antiferromagnetic Couplings in Organic Diradicals. J. Chem. Theory Comput. 2013, 9, 5216−5220. K

DOI: 10.1021/acs.jpca.6b09666 J. Phys. Chem. A XXXX, XXX, XXX−XXX

Article

The Journal of Physical Chemistry A Character of One- and Two-Dimensional Graphene Nanoribbons. Angew. Chem., Int. Ed. 2013, 52, 2581−2584. (40) Cui, Z.-H.; Lischka, H.; Mueller, T.; Plasser, F.; Kertesz, M. Study of the Diradicaloid Character in a Prototypical Pancake-Bonded Dimer: The Stacked Tetracyanoethylene (TCNE) Anion Dimer and the Neutral K2 TCNE2 Complex. ChemPhysChem 2014, 15, 165−176. (41) Machado, F. B. C.; Aquino, A. J. A.; Lischka, H. The Electronic States of a Double Carbon Vacancy Defect in Pyrene: a Model Study for Graphene. Phys. Chem. Chem. Phys. 2015, 17, 12778−12785. (42) Cho, D.; Ko, K. C.; Park, H.; Lee, J. Y. Ferromagnetic Graphene Nanoribbons: Edge Termination with Organic Radicals. J. Phys. Chem. C 2015, 119, 10109−10115. (43) Ma, H.; Liu, C.; Zhang, C.; Jiang, Y. Theoretical Study of Very High Spin Organic π-Conjugated Polyradicals. J. Phys. Chem. A 2007, 111, 9471−9478. (44) Gallagher, N. M.; Olankitwanit, A.; Rajca, A. High-Spin Organic Molecules. J. Org. Chem. 2015, 80, 1291−1298. (45) Salem, L.; Rowland, C. The Electronic Properties of Diradicals. Angew. Chem., Int. Ed. Engl. 1972, 11, 92−111. (46) Abe, M. Diradicals. Chem. Rev. 2013, 113, 7011−7088. (47) Rajca, A. Organic Diradicals and Polyradicals: From Spin Coupling to Magnetism? Chem. Rev. 1994, 94, 871−893. (48) Ovchinnikov, A. A. Multiplicity of the Ground State of Large Alternant Organic Molecules with Conjugated Bonds. Theor. Chim. Acta 1978, 47, 297−304. (49) Klein, D. J.; Yang, Y.; Ye, D. HOMO-LUMO Gaps for SubGraphenic and Sub-Buckytubic Species. Proc. R. Soc. London, Ser. A 2015, 471, 20150183. (50) Berson, J. A. A New Class of Non-Kekulé Molecules with Tunable Singlet−Triplet Energy Spacings. Acc. Chem. Res. 1997, 30, 238−244. (51) Cramer, C. J. Paul Dowd and Diradicals. J. Chem. Soc., Perkin Trans. 2 1998, 2, 1007−1014. (52) Lineberger, W. C.; Borden, T. W. The Synergy Between Qualitative Theory, Quantitative Calculations, and Direct Experiments in Understanding, Calculating, and Measuring the Energy Differences Between the Lowest Singlet and Triplet States of Organic Diradicals. Phys. Chem. Chem. Phys. 2011, 13, 11792−11813. (53) Barborini, M.; Coccia, E. Investigating Disjoint Non-Kekulé Diradicals with Quantum Monte Carlo: The Tetramethyleneethane Molecule through the Jastrow Antisymmetrized Geminal Power Wave Function. J. Chem. Theory Comput. 2015, 11, 5696−5704. (54) Borden, W. T.; Davidson, E. R. Theoretical Studies of Diradicals Containing Four.pi. Electrons. Acc. Chem. Res. 1981, 14, 69−76. (55) Gryn’ova, G.; Coote, M. L. Theory and Practice of Uncommon Molecular Electronic Configurations. WIREs Comput. Mol. Sci. 2015, 5, 440−459. (56) Ehrenfest, P. Opmerkingen over het Diamagnetisme van vast Bismuth. Physica 1925, 5, 388. (57) Pauling, L. The Diamagnetic Anisotropy of Aromatic Molecules. J. Chem. Phys. 1936, 4, 673−677. (58) London, F. Théorie Quantique des Courants Interatomiques dans les Combinaisons Aromatiques. J. Phys. Radium 1937, 8, 397− 409. (59) Pople, J. A. Molecular Orbital Theory of Aromatic Ring Currents. Mol. Phys. 1958, 1, 175−180. (60) Flygare, W. H. Magnetic Interactions in Molecules and an Analysis of Molecular Electronic Charge Distribution from Magnetic Parameters. Chem. Rev. 1974, 74, 653−687. (61) Schaad, L. J.; Hess, B. A., Jr. Dewar Resonance Energy. Chem. Rev. 2001, 101, 1465−1476. (62) Krygowski, T. M.; Cyrański, M. K. Structural Aspects of Aromaticity. Chem. Rev. 2001, 101, 1385−1420. (63) Schleyer, P. von R.; Maerker, C.; Dransfeld, A.; Jiao, H.; van Eikema Hommes, N. J. R. Nucleus-Independent Chemical Shifts: A Simple and Efficient Aromaticity Probe. J. Am. Chem. Soc. 1996, 118, 6317−6318.

(64) Chen, Z.; Wannere, C. S.; Corminboeuf, C.; Puchta, R.; Schleyer, P. v. R. Nucleus-Independent Chemical Shifts (NICS) as an Aromaticity Criterion. Chem. Rev. 2005, 105, 3842−3888. (65) Schleyer, P. v. R.; Manoharan, M.; Jiao, H.; Stahl, F. The Acenes: Is There a Relationship between Aromatic Stabilization and Reactivity? Org. Lett. 2001, 3, 3643−3646. (66) Fallah-Bagher-Shaidaei, H.; Wannere, C. S.; Corminboeuf, C.; Puchta, R.; Schleyer, P. v R. Which NICS Aromaticity Index for Planar π Rings Is Best? Org. Lett. 2006, 8, 863−866. (67) Corminboeuf, C.; Heine, T.; Seifert, G.; Schleyer, P. v. R.; Weber, J. Induced Magnetic Fields in Aromatic [n]-Annulenes Interpretation of NICS Tensor Components. Phys. Chem. Chem. Phys. 2004, 6, 273−276. (68) Ali, M. E.; Oppeneer, P. Influence of Noncovalent Cation/ Anion−π Interactions on the Magnetic Exchange Phenomenon. J. Phys. Chem. Lett. 2011, 2, 939−943. (69) Güell, M.; Poater, J.; Luis, J.; Mó, O.; Yáñez, M.; Solà, M. Aromaticity Analysis of Lithium Cation/ π Complexes of Aromatic Systems. ChemPhysChem 2005, 6, 2552−2561. (70) Barone, V.; Cacelli, I.; Ferretti, A. Magnetic Coupling in BisNitronylnitroxide Radicals: The Role of Aromatic Bridges. J. Chem. Phys. 2009, 130, 094306. (71) Bhattacharya, D.; Misra, A. Density Functional Theory Based Study of Magnetic Interaction in Bis-Oxoverdazyl Diradicals Connected by Different Aromatic Couplers. J. Phys. Chem. A 2009, 113, 5470−5475. (72) Latif, I.; Panda, A.; Datta, S. N. Very Strongly Ferromagnetically Coupled Diradicals from Mixed Radical Centers: Nitronyl Nitroxide Coupled to Oxoverdazyl via Polyene Spacers. J. Phys. Chem. A 2009, 113, 1595−1600. (73) Klein, D. J. Aromaticity via Kekule Structures and Conjugated Circuits. J. Chem. Educ. 1992, 69, 691−694. (74) Klein, D. J.; Babić, D. Partial Orderings in Chemistry. J. Chem. Inf. Comp. Sci. 1997, 37, 656−671. (75) Jug, K.; Köster, A. M. Aromaticity as a Multi-Dimensional Phenomenon. J. Phys. Org. Chem. 1991, 4, 163. (76) Bultinck, P.; Fias, S.; Ponec, R. Local Aromaticity in Polycyclic Aromatic Hydrocarbons: Electron Delocalization versus Magnetic Indices. Chem. - Eur. J. 2006, 12, 8813−8818. (77) Katritzky, A.; Barczynski, P.; Musumarra, G.; Pisano, D.; Szafran, M. Aromaticity as a Quantitative Concept. 1. A Statistical Demonstration of the Orthogonality of Classical and Magnetic Aromaticity in Five- and Six-Membered Heterocycles. J. Am. Chem. Soc. 1989, 111, 7−15. (78) Klein, D. J. Prolegomenon on Partial Orderings in Chemistry. Comm. Math. Computer Chemistry (MATCH) 2000, 42, 7−21. (79) Zeng, Z.; Shi, X.; Chi, C.; Navarrete, J. T. L.; Casado, J.; Wu, J. Pro-Aromatic and Anti-Aromatic π-Conjugated Molecules: an Irresistible Wish to be Diradicals. Chem. Soc. Rev. 2015, 44, 6578− 6596. (80) Gershoni-Poranne, R.; Stanger, A. Magnetic Criteria of Aromaticity. Chem. Soc. Rev. 2015, 44, 6597−6615. (81) Cocq, K.; Lepetit, C.; Maraval, V.; Chauvin, R. “Carboaromaticity” and novel carbo-aromatic compounds. Chem. Soc. Rev. 2015, 44, 6535−6559. (82) Tsipis, A. C. Efficiency of the NICSzz-scan curves to probe the antiaromaticity of organic and inorganic rings/cages. Phys. Chem. Chem. Phys. 2009, 11, 8244−8261. (83) Steiner, E.; Fowler, P. W. Patterns of Ring Currents in Conjugated Molecules: A Few-Electron Model Based on Orbital Contributions. J. Phys. Chem. A 2001, 105, 9553−9562. (84) Fowler, P. W.; Lilington, M.; Olson, L. P. Aromaticity, πelectron delocalization, and ring currents. Pure Appl. Chem. 2007, 79, 969−979. (85) Radenković, S.; Kojić, J.; Petronijević, J.; Antić, M. Effect of Benzo-Annelation on Local Aromaticity in Heterocyclic Conjugated Compounds. J. Phys. Chem. A 2014, 118, 11591−11601. (86) Radenković, S.; Tošović, J.; Nikolić, J. D. Local Aromaticity in Naphtho-Annelated Fluoranthenes: Can the Five-Membered Rings be L

DOI: 10.1021/acs.jpca.6b09666 J. Phys. Chem. A XXXX, XXX, XXX−XXX

Article

The Journal of Physical Chemistry A More Aromatic than the Six-Membered Rings? J. Phys. Chem. A 2015, 119, 4972−4982. (87) Islas, R.; Martínez-Guajardo, G.; Jiménez-Halla, J. O. C.; Solà, M.; Merino, G. Not All That Has a Negative NICS Is Aromatic: The Case of the H-Bonded Cyclic Trimer of HF. J. Chem. Theory Comput. 2010, 6, 1131−1135. (88) Foroutan-Nejad, C. Is NICS a Reliable Aromaticity Index for Transition Metal Clusters? Theor. Chem. Acc. 2015, 134, 8. (89) Trinquier, G.; Malrieu, J.-P. Kekulé versus Lewis: When Aromaticity Prevents Electron Pairing and Imposes Polyradical Character. Chem. - Eur. J. 2015, 21, 814−828. (90) Bhattacharya, D.; Panda, A.; Misra, A.; Klein, D. J. Clar Theory Extended for Polyacenes and Beyond. J. Phys. Chem. A 2014, 118, 4325−4338. (91) Singh, S. K.; Gupta, T.; Rajaraman, G. Magnetic Anisotropy and Mechanism of Magnetic Relaxation in Er(III) Single-Ion Magnets. Inorg. Chem. 2014, 53, 10835−10845. (92) Singh, S. K.; Rajaraman, G. Deciphering the Origin of Giant Magnetic Anisotropy and Fast Quantum Tunnelling in Rhenium(IV) Single-Molecule Magnets. Nat. Commun. 2016, 7, 10669. (93) Neese, F.; Pantazis, D. A. What is not Required to Make a Single Molecule Magnet. Faraday Discuss. 2011, 148, 229−238. (94) Ruiz, E.; Cirera, J.; Cano, J.; Alvarez, S.; Loose, C.; Kortus, J. Can Large Magnetic Anisotropy and High Spin Really Coexist? Chem. Commun. 2008, 52−54. (95) Singh, S. K.; Rajaraman, G. Probing the Origin of Magnetic Anisotropy in a Dinuclear {MnIIICuII} Single-Molecule Magnet: The Role of Exchange Anisotropy. Chem. - Eur. J. 2014, 20, 5214−5218. (96) Waldmann, O. A. A Criterion for the Anisotropy Barrier in Single-Molecule Magnets. Inorg. Chem. 2007, 46, 10035−10037. (97) Cirera, J.; Ruiz, E.; Alvarez, S.; Neese, F.; Kortus, J. How to Build Molecules with Large Magnetic Anisotropy. Chem. - Eur. J. 2009, 15, 4078−4087. (98) Thomson, A. J.; Greenwood, C.; Gadsby, P. M. A.; Noote, N. “Nature of Ferryl, Fe(IV) = 0, Haem in Peroxidases. In Cytochrome Systems: Molecular Biology and Bioenergetics;Papa, S., Chance, B., Ernster, L., Eds.; Plenum Press: New York, 1987; p 349. (99) Long, J. R. Molecular Cluster Magnets. In The Chemistry of Nanostructured Materials; Yang, P., Ed.; World Scientific Publishing Co.: Singapore, 2003; p 291. (100) Sandberg, K. A. Ph.D. Thesis, Department of Chemistry, NC State University, 1998. (101) Feng, X.; Liu, J.; Harris, T. D.; Hill, S.; Long, J. R. Slow Magnetic Relaxation Induced by a Large Transverse Zero-Field Splitting in a MnIIReIV(CN)2 Single-Chain Magnet. J. Am. Chem. Soc. 2012, 134, 7521−7529. (102) Pedersen, K. S.; Bendix, J.; Clérac, R. Single-Molecule Magnet Engineering: Building-Block Approaches. Chem. Commun. 2014, 50, 4396−4415. (103) Tait, C. E.; Neuhaus, P.; Anderson, H. L.; Timmel, C. R. Triplet State Delocalization in a Conjugated Porphyrin Dimer Probed by Transient Electron Paramagnetic Resonance Techniques. J. Am. Chem. Soc. 2015, 137, 6670−6679. (104) Angiolillo, P. J.; Uyeda, H. T.; Duncan, T. V.; Therien, M. J. Impact of Electronic Asymmetry on Photoexcited Triplet-State Spin Distributions in Conjugated Porphyrin Oligomers Probed via EPR Spectroscopy. J. Phys. Chem. B 2004, 108, 11893−11903. (105) Bennati, M.; Grupp, A.; Mehring, M.; Bäuerle, P. Pulsed EPR Spectroscopy of the Photoexcited Triplet States of Thiophene Oligomers in Frozen Solution. J. Phys. Chem. 1996, 100, 2849−2853. (106) Riplinger, C.; Kao, J. P. Y.; Rosen, G. M.; Kathirvelu, V.; Eaton, G. R.; Eaton, S. S.; Kutateladze, A.; Neese, F. Interaction of Radical Pairs Through-Bond and Through-Space: Scope and Limitations of the Point−Dipole Approximation in Electron Paramagnetic Resonance Spectroscopy. J. Am. Chem. Soc. 2009, 131, 10092−10106. (107) Beer, P. D.; Gale, P. A. Anion Recognition and Sensing: The State of the Art and Future Perspectives. Angew. Chem., Int. Ed. 2001, 40, 486−516.

(108) Garau, C.; Frontera, A.; Quiñonero, D.; Ballester, P.; Costa, A.; Deyá, P. M. Cation-π vs Anion-π Interactions: a Complete π-Orbital Analysis. Chem. Phys. Lett. 2004, 399, 220−225. (109) Frański, R.; Gierczyk, B.; Schroeder, G. Anion−π InteractionsInteractions Between Benzo-Crown Ether Metal Cation Complexes and Counter Ions. J. Am. Soc. Mass Spectrom. 2009, 20, 257−262. (110) Gamez, P.; Mooibroek, T. J.; Teat, S. J.; Reedijk, J. Anion Binding Involving π-Acidic Heteroaromatic Rings. Acc. Chem. Res. 2007, 40, 435−444. (111) Schottel, B. L.; Chifotides, H. T.; Dunber, K. R. Anion-π Interactions. Chem. Soc. Rev. 2008, 37, 68−83. (112) Bytautas, L. Stability of Dimethylmercury and Related Mercury-Containing Compounds with Respect to Selected Chemical Species Found in Aqueous Environment. Croat. Chem. Acta 2013, 86, 453−462. (113) Ishiguro, K.; Ozaki, M.; Sekine, N.; Sawaki, Y. Acid−Base Equilibrium between Phenoxyl−Nitronyl Nitroxide Biradical and Closed-Shell Cation. A Magnetic pH Sensor. J. Am. Chem. Soc. 1997, 119, 3625−3626. (114) Gryn’ova, G.; Marshall, D. L.; Blanksby, S. J.; Coote, M. L. Switching Radical Stability by pH-Induced Orbital Conversion. Nat. Chem. 2013, 5, 474−481. (115) Lee, C.; Yang, W.; Parr, R. G. Development of the ColleSalvetti Correlation-Energy Formula into a Functional of the Electron Density. Phys. Rev. B: Condens. Matter Mater. Phys. 1988, 37, 785−789. (116) Becke, A. D. Density-Functional Thermochemistry. III. The Role of Exact Exchange. J. Chem. Phys. 1993, 98, 5648−5652. (117) Hehre, W. J.; Ditchfield, R.; Pople, J. A. SelfConsistent Molecular Orbital Methods. XII. Further Extensions of Gaussian Type Basis Sets for Use in Molecular Orbital Studies of Organic Molecules. J. Chem. Phys. 1972, 56, 2257−2261. (118) Hariharan, P. C.; Pople, J. A. The Influence of Polarization Functions on Molecular Orbital Hydrogenation Energies. Theor. Chim. Acta 1973, 28, 213−222. (119) Zhao, Y.; Truhlar, D. G. The M06 Suite of Density Functionals for Main Group Thermochemistry, Thermochemical Kinetics, Noncovalent Interactions, Excited States, and Transition Elements: Two New Functionals and Systematic Testing of Four M06-Class Functionals and 12 Other Functionals. Theor. Chem. Acc. 2008, 120, 215−241. (120) Zhao, Y.; Truhlar, D. G. Applications and Validations of the Minnesota Density Functionals. Chem. Phys. Lett. 2011, 502, 1−13. (121) Polo, V.; Alberola, A.; Andres, J.; Anthony, J.; Pilkington, M. Towards Understanding of Magnetic Interactions Within a Series of Tetrathiafulvalene−π Conjugated-Verdazyl Diradical Cation System: a Density Functional Theory Study. Phys. Chem. Chem. Phys. 2008, 10, 857−864. (122) Bhattacharya, D.; Shil, S.; Misra, A.; Klein, D. J. Intramolecular Ferromagnetic Coupling in Bis-Oxoverdazyl and Bis-Thioxoverdazyl Diradicals with Polyacene Spacers. Theor. Chem. Acc. 2010, 127, 57− 67. (123) Bhattacharya, D.; Shil, S.; Panda, A.; Misra, A. A DFT Study on the Magnetostructural Property of Ferromagnetic Heteroverdazyl Diradicals with Phenylene Coupler. J. Phys. Chem. A 2010, 114, 11833−11841. (124) Bhattacharya, D.; Shil, S.; Misra, A. Photoresponsive Magnetization Reversal in Green Fluorescent Protein Chromophore Based Diradicals. J. Photochem. Photobiol., A 2011, 217, 402−410. (125) Bhattacharya, D.; Panda, A.; Shil, S.; Goswami, T.; Misra, A. A Theoretical Study on Photomagnetic Fluorescent Protein Chromophore Coupled Diradicals and their Possible Applications. Phys. Chem. Chem. Phys. 2012, 14, 6905−6913. (126) Yamaguchi, K.; Takahara, Y.; Fueno, T.; Nasu, K. Ab Initio MO Calculations of Effective Exchange Integrals between TransitionMetal Ions via Oxygen Dianions: Nature of the Copper-Oxygen Bonds and Superconductivity. Jpn. J. Appl. Phys. 1987, 26, L1362−L1364. (127) Yamaguchi, K.; Jensen, F.; Dorigo, A.; Houk, K. N. A Spin Correction Procedure for Unrestricted Hartree-Fock and MøllerM

DOI: 10.1021/acs.jpca.6b09666 J. Phys. Chem. A XXXX, XXX, XXX−XXX

Article

The Journal of Physical Chemistry A Plesset Wavefunctions for Singlet Diradicals and Polyradicals. Chem. Phys. Lett. 1988, 149, 537−542. (128) Ginsberg, A. P. Magnetic Exchange in Transition Metal Complexes. 12. Calculation of Cluster Exchange Coupling Constants with the X.alpha.-Scattered Wave Method. J. Am. Chem. Soc. 1980, 102, 111−117. (129) Sinnecker, S.; Neese, F. Spin−Spin Contributions to the ZeroField Splitting Tensor in Organic Triplets, Carbenes and Biradicals−A Density Functional and Ab Initio Study. J. Phys. Chem. A 2006, 110, 12267−12275. (130) Pederson, M. R.; Khanna, S. N. Magnetic Anisotropy Barrier for Spin Tunneling in Mn12O12 Molecules. Phys. Rev. B: Condens. Matter Mater. Phys. 1999, 60, 9566−9572. (131) Frisch, M. J.; Trucks, G. W.; Schlegel, H. B.; Scuseria, G. E.; Robb, M. A.; Cheeseman, J. R.; Scalmani, G.; Barone, V.; Mennucci, B.; Petersson, G. A.; Nakatsuji, H.; et al. Gaussian’09W, Revision A.02; Gaussian, Inc.: Wallingford, CT, 2009. (132) Neese, F. ORCA − An Ab-initio Density Functional and Semiempirical Program Package, Version 2.8−20; Max-Planck Institute for Bioinorganic Chemistry: Mulheim an der Ruhr, Germany, 2010. (133) Aquino, F.; Rodriguez, J. H. Accurate Calculation of Zero-Field Splittings of (Bio)inorganic Complexes: Application to an {FeNO}7 (S = 3/2) Compound. J. Phys. Chem. A 2009, 113, 9150−9156. (134) Takeda, R.; Mitsuo, S.; Yamanaka, S.; Yamaguchi, K. Density functional study of zero-field splitting. Polyhedron 2005, 24, 2238− 2241. (135) Kortus, J.; Pederson, M. R.; Baruah, T.; Bernstein, N.; Hellberg, C. S. Density Functional Studies of Single Molecule Magnets. Polyhedron 2003, 22, 1871−1876. (136) Park, K.; Pederson, M. R.; Richardson, S. L.; Aliaga-Alcalde, N.; Christou, G. Density-Functional Theory Calculation of the Intermolecular Exchange Interaction in the Magnetic Mn4 Dimer. Phys. Rev. B: Condens. Matter Mater. Phys. 2003, 68, 020405. (137) Baruah, T.; Pederson, M. R. Density Functional Study of the Conformers of Co4-Based Single-Molecule Magnet. Int. J. Quantum Chem. 2003, 93, 324−331. (138) Loboda, O.; Minaev, B.; Vahtras, O.; Schimmelpfennig, B.; Ǻ gren, H.; Ruud, K.; Jonsson, D. Ab Initio Calculations of Zero-Field Splitting Parameters in Linear Polyacenes. Chem. Phys. 2003, 286, 127−137. (139) Boča, R. Theoretical Foundations on Molecular Magnetism; Elsevier: Amsterdam, 1999. (140) Kruszewski, J.; Krygowski, T. M. Definition of Aromaticity Basing on the Harmonic Oscillator Model. Tetrahedron Lett. 1972, 13, 3839−3842. (141) Krygowski, T. M. Crystallographic Studies of Inter- and Intramolecular Interactions Reflected in Aromatic Character of.pi.Electron Systems. J. Chem. Inf. Comput. Sci. 1993, 33, 70−78. (142) Wolinski, K.; Hinton, J. F.; Pulay, P. Efficient Implementation of the Gauge-Independent Atomic Orbital Method for NMR Chemical Shift Calculations. J. Am. Chem. Soc. 1990, 112, 8251−8260. (143) Rajca, A.; Wongsriratanakul, J.; Rajca, S. Magnetic Ordering in an Organic Polymer. Science 2001, 294, 1503−1505. (144) Rajca, A. The Physical Organic Chemistry of Very High-Spin Polyradicals. Adv. Phys. Org. Chem. 2005, 40, 153−199. (145) Herrmann, C.; Solomon, G. C.; Ratner, M. A. Organic Radicals As Spin Filters. J. Am. Chem. Soc. 2010, 132, 3682−3684. (146) Sarbadhikary, P.; Shil, S.; Panda, A.; Misra, A. A Perspective on Designing Chiral Organic Magnetic Molecules with Unusual Behavior in Magnetic Exchange Coupling. J. Org. Chem. 2016, 81, 5623−5630. (147) Constantinides, C. P.; Koutentis, P. A.; Schatz, J. A DFT Study of the Ground State Multiplicities of Linear vs Angular Polyheteroacenes. J. Am. Chem. Soc. 2004, 126, 16232−16241. (148) Hund, F. Zur Deutung der Molekelspektren. IV. Eur. Phys. J. A 1928, 51, 759−795. (149) Hoffmann, R.; Zeiss, G. D.; Van Dine, G. W. The Electronic Structure of Methylenes. J. Am. Chem. Soc. 1968, 90, 1485−1499. (150) Herges, R. Magnetic Properties of Aromatic Compounds and Aromatic Transition States, in The Chemical Bond: Chemical Bonding

Across the Periodic Table; Frenking, G., Shaik, S., Ed.; Wiley-VCH Verlag GmbH & Co. KGaA: Weinheim, Germany, 2014; doi 10.1002/ 9783527664658.ch13.

N

DOI: 10.1021/acs.jpca.6b09666 J. Phys. Chem. A XXXX, XXX, XXX−XXX