From Magnets to Metals: The Response of Tetragonal Bisdiselenazolyl Radicals to Pressure

Magnetic coupling and spin ordering in bisdithiazolyl, thiaselenazolyl, and bisdiselenazolyl molecular materials

The use of purely organic materials is a promising approach for the miniaturization of devices due to their interesting optical, electronic and magnetic properties. Bisdithiazolyl-based bisDTA compounds have emerged as promising candidates for radical-based single component conductors exhibiting simultaneously magnetic properties. Our computational work focuses on the intriguing magnetism of 4 isostructural pyridine-bridged bisDTA-multifunctional materials triggered by their magnetic and conducting properties being strongly dependent on the different S/Se ratios in the neutral radical skeleton: specifically, bisdithiazolyl (S,S) displays no magnetic order at low temperatures, thiaselenazolyl (Se,S) exhibits spin-canted antiferromagnetism (AFM), and both (S,Se) and bisdiselenazolyl (Se,Se) behave as bulk ferromagnets (FM). Our results reveal that (1) the magnetic response depends on the existence of an intricate network of both AFM and FM spin exchange J_AB couplings between neighbouring radicals; and (2) the structural arrangement of π-stacked pairs of radicals sits on a point in the configurational space that is very close to a crossover region where J_AB switches from AFM to FM. Indeed, for bulk FM, the experimental response is only accounted for when considering an ab initio optimised crystal structure able to portray adequately the electronic structure of bisDTAs in the region close to the temperature at which magnetic ordering emerges. Magneto-structural correlation maps show the large sensitivity of J_AB to very small structural changes with temperature along the π-stacks that lead to drastic changes in the magnetic properties. Clearly, the understanding of magnetism in the title bisDTA compounds is decisive to rationally tailor the properties of multifunctional materials by subtle structural modifications of their crystal packing.

Electronic structure and magnetic coupling in selenium substituted pyridine-bridged bisdithiazolyl multifunctional molecular materials

Bisdithiazolyl radicals have furnished in recent years multiple examples of molecular materials with promising conductive and magnetic properties. The electronic band structure and magnetic ordering in four different isostructural pyridine-bridged bisdithiazolyl and Selenium substituted compounds have been studied by means of hybrid DFT based methods as implemented in the CRYSTAL code. The full rationalization of the properties of these multifunctional magnetic molecular materials requires a careful description of their complex open-shell electronic structure. The results describe the systems as narrow band (0.2-0.3 eV dispersion) open-shell semiconductors with a gap of 1.15-1.40 eV between the valence and conducting bands. The bands defining the insulating gap are dominated by orbital contributions arising from the heteroatoms sitting in the outer rings. A low energy closed-shell metallic solution is found at 0.25-0.35 eV above the magnetic solutions thus suggesting a complex mechanism for electric conduction with band and hopping contributions. The observed trend of the conductivity is in line with the variation of the insulating gap but more rigorous modelling is required to take into account the details of the band structure of the systems. For all the systems the spin density is well localised on the molecular units and is independent of the magnetic solution. Thus the system can be described as an ensemble of well-defined S = 1/2 magnetic centres using a two-body Heisenberg-Dirac-van Vleck spin Hamiltonian. The lowest energy electronic solutions are in line with the observed magnetic behaviour at low temperature. The set of competing magnetic exchange interactions that emerges from using a suitable mapping to consistently describe the low energy magnetic solutions explains the variety of magnetic responses (absence of long-range magnetic order, antiferromagnetism or ferromagnetism) of the four studied compounds at low temperatures.

A first-order phase transition in Blatter's radical at high pressure

The crystal structure of Blatter's radical (1,3-diphenyl-1,4-dihydrobenzo[e][1,2,4]triazin-4-yl) has been investigated between ambient pressure and 6.07 GPa. The sample remains in a compressed form of the ambient-pressure phase up to 5.34 GPa, the largest direction of strain being parallel to the direction of π-stacking interactions. The bulk modulus is 7.4 (6) GPa, with a pressure derivative equal to 9.33 (11). As pressure increases, the phenyl groups attached to the N1 and C3 positions of the triazinyl moieties of neighbouring pairs of molecules approach each other, causing the former to begin to rotate between 3.42 to 5.34 GPa. The onset of this phenyl rotation may be interpreted as a second-order phase transition which introduces a new mode for accommodating pressure. It is premonitory to a first-order isosymmetric phase transition which occurs on increasing pressure from 5.34 to 5.54 GPa. Although the phase transition is driven by volume minimization, rather than relief of unfavourable contacts, it is accompanied by a sharp jump in the orientation of the rotation angle of the phenyl group. DFT calculations suggest that the adoption of a more planar conformation by the triazinyl moiety at the phase transition can be attributed to relief of intramolecular H...H contacts at the transition. Although no dimerization of the radicals occurs, the π-stacking interactions are compressed by 0.341 (3) Å between ambient pressure and 6.07 GPa.

Low temperature insights into the crystal and magnetic structure of a neutral radical ferromagnet

The crystal structure of the radical ferromagnet 1a at 2 K reveals a contraction in the unit cell c constant which, at the molecular level, translates into a decrease in slippage of the radical π-stacks and an increase in ferromagnetic exchange interactions along the stacking axis. The results of BS-DFT calculations using long-range corrected functionals are consistent with an overall ferromagnetic topology.

Heisenberg Spin Chains via Chalcogen Bonding: Noncovalent S···O Contacts Enable Long-Range Magnetic Order

The new radical ligand 5,8-dimethyl-1,4-dioxonaphtho[2,3-d][1,2,3]dithiazolyl (1) is reported. Two crystal polymorphs, 1α and 1β, differing in their pancake-bonded dimerization motif and S···O contact network, are identified. The self-assembly of Mn(II) metal ions with 1 leads to the formation of [Mn(hfac)₂]₃(1)₂ that exhibits a Mn(II)-radical-Mn(II)-radical-Mn(II) linear arrangement of three Mn(hfac)₂ units bridged by two radical ligands (hfac = 1,1,1,5,5,5-hexafluoroacetylacetonato-). Characterization by single-crystal X-ray diffraction of this Mn(II) complex packing structure reveals close noncovalent S···O contacts between the [Mn(hfac)₂]₃(1)₂ units in one dimension along the b-c direction. The magnetic properties of the coordination complex are characterized by dc and ac susceptibility measurements on a microcrystalline solid. The magnetic data down to 4.8 K indicate the presence of effective ferromagnetic interactions (J/k_B = +0.16 K) between the molecular S_T = 13/2 units along the supramolecular chain involving noncovalent S···O contacts. Below 2.9 K, a non-zero out-of-phase component appears in the ac susceptibility, indicating the presence of a three-dimensional magnetic phase transition.

Crystal structure and metallization mechanism of the π-radical metal

Radical electrons tend to localize on individual molecules, resulting in an insulating (Mott-Hubbard) bandgap in the solid state. Herein, we report the crystal structure and intrinsic electronic properties of the first single crystal of a π-radical metal, tetrathiafulvalene-extended dicarboxylate (TED). The electrical conductivity is up to 30 000 S cm-1 at 2 K and 2300 S cm-1 at room temperature. Temperature dependence of resistivity obeys a T 3 power-law above T > 100 K, indicating a new type of metal. X-ray crystallographic analysis clarifies the planar TED molecule, with a symmetric intramolecular hydrogen bond, is stacked along longitudinal (the a-axis) and transverse (the b-axis) directions. The π-orbitals are distributed to avoid strong local interactions. First-principles electronic calculations reveal the origin of the metallization giving rise to a wide bandwidth exceeding 1 eV near the Fermi level. TED demonstrates the effect of two-dimensional stacking of π-orbitals on electron delocalization, where a high carrier mobility of 31.6 cm2 V-1 s-1 (113 K) is achieved.

A chemical perspective on high pressure crystal structures and properties

The general availability of third generation synchrotron sources has ushered in a new era of high pressure research. The crystal structure of materials under compression can now be determined by X-ray diffraction using powder samples and, more recently, from multi-nano single crystal diffraction. Concurrently, these experimental advancements are accompanied by a rapid increase in computational capacity and capability, enabling the application of sophisticated quantum calculations to explore a variety of material properties. One of the early surprises is the finding that simple metallic elements do not conform to the general expectation of adopting 3D close-pack structures at high pressure. Instead, many novel open structures have been identified with no known analogues at ambient pressure. The occurrence of these structural types appears to be random with no rules governing their formation. The adoption of an open structure at high pressure suggested the presence of directional bonds. Therefore, a localized atomic hybrid orbital description of the chemical bonding may be appropriate. Here, the theoretical foundation and experimental evidence supporting this approach to the elucidation of the high pressure crystal structures of group I and II elements and polyhydrides are reviewed. It is desirable and advantageous to extend and apply established chemical principles to the study of the chemistry and chemical bonding of materials at high pressure.

Effect of Structure on the Spin Switching and Magnetic Bistability of Solid-State Aryl Dicyanomethyl Monoradicals and Diradicals

Stable organic radicals with switchable spin states have applications in medicine, biology, and material science. An emerging class of such spin-switchable radicals is based on dicyanomethyl radicals, which are typically thermally and air-stable species that form weakly bonded (closed-shell singlet) dimers at a lower temperature that rupture into electron paramagnetic resonance-active diradicals at a higher temperature. However, thus far, the study of these dicyanomethyl radicals has focused on their solution-phase behavior. An understanding of how chemical structure affects the solid-state spin switching behavior for these radicals is unknown. Here, we examine the solid-state spin crossover behavior of 6 monoradicals and 10 tethered diradicals and demonstrate that these species also undergo spin switching in the solid state. We find that the susceptibility for solid-state spin switching for the intermolecular dimers is weakly correlated to the solution-phase Gibbs free energies of dimerization, but no apparent correlations are seen between the solution-state free energies for the intramolecular dimerization and the solid-state behavior. Furthermore, intramolecular diradical dimers have greatly enhanced temperature-responsive behavior compared to their intermolecular counterparts. Crystalline and amorphous powders of the same radicals feature similar spin switching behavior, but the crystalline materials have slower bond-rupture kinetics at higher temperatures, suggesting that solid-state packing effects are an important kinetic consideration. An interesting feature of these systems is that, upon cooling down to room temperature after heating, some radicals remain trapped in the solids, indicating magnetic bistability, while others partially or fully return to the diamagnetic dimers. This work provides insights into how chemical structure affects spin crossover in the solid state for this new class of air-stable radicals, the knowledge of importance for the construction of dynamically responsive solid-state materials, and organic spin crossover polymers.

The Importance of Electronic Dimensionality in Multiorbital Radical Conductors

The exceptional performance of oxobenzene-bridged bis-1,2,3-dithiazolyls 6 as single-component neutral radical conductors arises from the presence of a low-lying π-lowest unoccupied molecular orbital, which reduces the potential barrier to charge transport and increases the kinetic stabilization energy of the metallic state. As part of ongoing efforts to modify the solid-state structures and transport properties of these so-called multiorbital materials, we report the preparation and characterization of the acetoxy, methoxy, and thiomethyl derivatives 6 (R = OAc, OMe, SMe). The crystal structures are based on ribbonlike arrays of radicals laced together by S···N' and S···O' secondary bonding interactions. The steric and electronic effects of the exocyclic ligands varies, affording one-dimensional (1D) π-stacked radicals for R = OAc, 1D cofacial dimer π-stacks for R = SMe, and a pseudo two-dimensional (2D) brick-wall arrangement for R = OMe. Variable-temperature magnetic and conductivity measurements reveal strong antiferromagnetic interactions and Mott insulating behavior for the two radical-based structures (R = OAc, OMe), with lower room-temperature conductivities (σ_RT ≈ 1 × 10^-4 and ∼1 × 10^-3 S cm^-1, respectively) and higher thermal activation energies ( E_act = 0.24 and 0.21 eV, respectively) than found for the ideal 2D brick-wall structure of 6 (R = F), where σ_RT ≈ 1 × 10^-2 S cm^-1 and E_act = 0.10 eV. The performance of R = OMe, OAc relative to that of R = F, is consistent with the results of density functional theory band electronic structure calculations, which indicate a lower kinetic stabilization energy of the putative metallic state arising from their reduced electronic dimensionality.

Pancakes under Pressure: A Case Study on Isostructural Dithia- and Diselenadiazolyl Radical Dimers

The isostructural dimers of the 1,4-phenylene-bridged bis-1,2,3,5-dithia- and bis-1,2,3,5-diselenadiazolyl diradicals 1,4-S/Se are small band gap semiconductors. The response of their molecular and solid state electronic structures to pressure has been explored over the range 0-10 GPa. The crystal structures, which consist of cofacially aligned (pancake) π-dimers packed into herringbone arrays, experience a continuous, near-isotropic compression. While the intramolecular covalent E-E (E = S/Se) bonds remain relatively unchanged with pressurization, the intradimer E···E separations are significantly shortened. Molecular and band electronic structure calculations using density functional theory methods indicate that compression of the π-dimers leads to a widening of the gap Δ E between the highest occupied and lowest unoccupied molecular orbitals of the dimer, an effect that offsets the expected decrease in the valence-to-conduction band gap E_g occasioned by pressure-induced spreading of the valence and conduction bands. Consistent with the predicted consequences of this competition between intra- and interdimer interactions, variable temperature high pressure conductivity measurements reveal at best an order-of-magnitude increase in conductivity with pressure for the two compounds over the pressure range 0-10 GPa. While a small reduction in the thermal activation energy E_act with increasing pressure is observed, extrapolation of the rate of decrease suggests a projected onset of metallization ( E_act ≈ 0) in excess of 20 GPa.

Three-Dimensional Magnetic Exchange Networks in Trigonal Bisdithiazolyl Radicals

The N-methyl-4-phenyl-pyridine-bridged bisdithiazolyl radical PhBPMe is polymorphic, crystallizing from cold acetonitrile in a trigonal α-phase, space group P3₁21, and from hot dichloroethane in an orthorhombic β-phase, space group Pca2₁. The crystal structures of both phases consist of slipped π-stacks of undimerized radicals aligned laterally into herringbone arrays. In the β-phase, there are two independent radicals in the asymmetric unit, and the resulting π-stacks form corrugated layers interspersed by methyl and phenyl groups which block the approach of neighboring radicals. In the α-phase, the methyl/phenyl groups and the radical π-stacks separately form spirals about 3₁ axes, the latter giving rise to a 3D network of close radical/radical contacts. Variable temperature magnetic susceptibility measurements on the β-phase indicate strong antiferromagnetic coupling. Weaker but predominantly antiferromagnetic interactions (θ = -20.7 K) are observed in the α-phase. A high temperature series expansion analysis of the magnetic data for the α-phase affords antiferromagnetic exchange energies for the one- and two-step radical/radical interactions about the 3₁ spirals ( J₁ = -1.2 K, J₂ = -10.9 K, respectively), with weak ferromagnetic interactions along the π-stacks ( J_π = +1.8 K). Despite the presence of a 3D network based on the dominant J₂ interactions, which affords two independent bipartite sublattices, no evidence of bulk antiferromagnetic order has been observed above T = 2 K. The magnetic results are discussed in light of exchange energies calculated using density functional theory broken symmetry methods.

Ferromagnetic interactions in a 1D Heisenberg linear chain of 1-phenyl-3,7-bis(trifluoromethyl)-1,4-dihydro-1,2,4-benzotriazin-4-yls

1-Phenyl-3,7-bis(trifluoromethyl)-1,4-dihydro-1,2,4-benzotriazin-4-yl (2), was characterized by single crystal X-ray diffractometry and variable temperature SQUID magnetometry to investigate its structure-magnetism correlation.

Benzoquinone-Bridged Heterocyclic Zwitterions as Building Blocks for Molecular Semiconductors and Metals

In pursuit of closed-shell building blocks for single-component organic semiconductors and metals, we have prepared benzoquino-bis-1,2,3-thiaselenazole QS, a heterocyclic selenium-based zwitterion with a small gap (λ_max = 729 nm) between its highest occupied and lowest unoccupied molecular orbitals. In the solid state, QS exists in two crystalline phases and one nanocrystalline phase. The structures of the crystalline phases (space groups R3 c and P2₁/ c) have been determined by high-resolution powder X-ray diffraction methods at ambient and elevated pressures (0-15 GPa), and their crystal packing patterns have been compared with that of the related all-sulfur zwitterion benzoquino-bis-1,2,3-dithiazole QT (space group Cmc2₁). Structural differences between the S- and Se-based materials are interpreted in terms of local intermolecular S/Se···N'/O' secondary bonding interactions, the strength of which varies with the nature of the chalcogen (S vs Se). While the perfectly two-dimensional "brick-wall" packing pattern associated with the Cmc2₁ phase of QT is not found for QS, all three phases of QS are nonetheless small band gap semiconductors, with σ_RT ranging from 10^-5 S cm^-1 for the P2₁/ c phase to 10^-3 S cm^-1 for the R3 c phase. The bandwidths of the valence and conduction bands increase with applied pressure, leading to an increase in conductivity and a decrease in thermal activation energy E_act. For the R3 c phase, band gap closure to yield an organic molecular metal with a σ_RT of ∼10² S cm^-1 occurs at 6 GPa. Band gaps estimated from density functional theory band structure calculations on the ambient- and high-pressure crystal structures of QT and QS correlate well with those obtained experimentally.

Synthesis of new hybrid 1,4-thiazinyl-1,2,3-dithiazolyl radicals via Smiles rearrangement

The condensation reaction of 2-aminobenzenethiols and 3-aminopyrazinethiols with 2-amino-6-fluoro-N-methylpyridinium triflate afforded thioether derivatives that were found to undergo Smiles rearrangement and cyclocondensation with sulphur monochloride to yield new hybrid 1,4-thiazine-1,2,3-dithiazolylium cations. The synthesized cations were readily reduced to the corresponding stable neutral radicals with spin densities delocalized over both 1,4-thiazinyl and 1,2,3-dithiazolyl moieties.

Singlet fission in pancake-bonded systems

Open-shell aggregates with pancake bonding are found to cause highly efficient singlet fission and large charge transport simultaneously.

The quest for rationalizing the magnetism in purely organic semiquinone-bridged bisdithiazolyl molecular magnets

Semiquinone-bridged bisdithiazolyl-based radicals (XBBO) are appealing purely organic magnetic building blocks for the synthesis of new functional materials. Remarkably, for the phenyl-derivative PhBBO, the rationalization of its magnetism becomes a proof of concept that DFT can dramatically fail to evaluate JAB magnetic interactions between purely organic radical pairs. Instead, wavefunction-based methods are required. Once JAB's are fully characterized, the magnetic topology of PhBBO is disclosed to consist of ferromagnetic FM π-stacks that are very weakly coupled (by FM and AFM JAB interactions). The magnetic susceptibility χT(T) and magnetization M(H) of PhBBO are then calculated using a first-principles bottom-up approach. The study of the unit cell contraction upon cooling from room temperature to zero-Kelvin is relevant to propose a suitable model for the phase transition that occurs at 4.5 K. A simplistic picture tells us that the antiparallel-aligned 1D-FM-chains convert into domains of weakly either FM- or AFM-coupled 1D-FM-chains. Accordingly, the presence of these domains may introduce geometrical spin frustration below 4.5 K.

Halide Influence on Molecular and Supramolecular Arrangements of Iron Complexes with a 3,5-Bis(2-Pyridyl)-1,2,4,6-Thiatriazine Ligand

A series of iron centered complexes, namely, [Fe(Py2TTA)Cl2] (1), [Fe(Py2TTA)Br2] (2), and [Fe(μ-F)(Py2TTAO)F]∞ (3), were isolated via complexation of 3,5-bis(2-pyridyl)-1,2,4,6-thiatriazine (Py2TTAH) with various ferric halides (e.g., FeF3, FeCl3, and FeBr3). Comparison of the optical and electrochemical spectroscopy, structural analysis, and magnetic studies reveal numerous similarities between the chlorido (1) and bromido (2) derivatives, which crystallize as discrete five-coordinate iron centered complexes with coordination geometries that are intermediate between trigonal bipyramidal and square base pyramid. Conversely, the fluorido derivative (3) results in a completely different structure due to oxidation of the ligand and the formation of a one-dimensional coordination polymer held together through a bridging fluoride ion. Consequently, the spectroscopic and magnetic behavior of 3 differs significantly compared with 1 and 2. Complexes 1 and 2 exhibit paramagnetic properties typical for a mononuclear S = 5/2 system with weak intermolecular antiferromagnetic interactions at low temperatures, whereas complex 3 demonstrates significant exchange couplings within the chain and weak antiferromagnetic interchain interactions, which stabilize a canted antiferromagnetic state below 4.2 K.

Pushing TC to 27.5 K in a heavy atom radical ferromagnet

In the solid state the iodo-substituted bisdiselenazolyl radical 1c orders as a bulk ferromagnet with T_C = 10.5 K. With the application of pressure T_C rises rapidly, reaching a value of 27.5 K at 2.4 GPa.

The metallic state in neutral radical conductors: dimensionality, pressure and multiple orbital effects

Pressure-induced changes in the solid-state structures and transport properties of three oxobenzene-bridged bisdithiazolyl radicals 2 (R = H, F, Ph) over the range 0-15 GPa are described. All three materials experience compression of their π-stacked architecture, be it (i) 1D ABABAB π-stack (R = Ph), (ii) quasi-1D slipped π-stack (R = H), or (iii) 2D brick-wall π-stack (R = F). While R = H undergoes two structural phase transitions, neither of R = F, Ph display any phase change. All three radicals order as spin-canted antiferromagnets, but spin-canted ordering is lost at pressures <1.5 GPa. At room temperature, their electrical conductivity increases rapidly with pressure, and the thermal activation energy for conduction Eact is eliminated at pressures ranging from ∼3 GPa for R = F to ∼12 GPa for R = Ph, heralding formation of a highly correlated (or bad) metallic state. For R = F, H the pressure-induced Mott insulator to metal conversion has been tracked by measurements of optical conductivity at ambient temperature and electrical resistivity at low temperature. For R = F compression to 6.2 GPa leads to a quasiquadratic temperature dependence of the resistivity over the range 5-300 K, consistent with formation of a 2D Fermi liquid state. DFT band structure calculations suggest that the ease of metallization of these radicals can be ascribed to their multiorbital character. Mixing and overlap of SOMO- and LUMO-based bands affords an increased kinetic energy stabilization of the metallic state relative to a single SOMO-based band system.

Ferromagnetic ordering in the organic radical cation salt BBDTA·Au(CN)2at 8.2 K

An organic radical cation salt, BBDTA·Au(CN)₂, with a slipped π-stacking columnar structure and intercolumnar short contacts, shows ferromagnetic ordering at 8.2 K, the highest reported temperature among the BBDTA⁺cation salts.

Electrochemically informed synthesis: oxidation versus coordination of 5,6-bis(phenylchalcogeno)acenaphthenes

Chalcogen dications: Facile synthesis of E--E bonded dications can be readily achieved. Radical cations are identified as the intermediates.

Synthesis of tetrachalcogenide-substituted phenalenyl derivatives: preparation and solid-state characterization of bis(3,4,6,7-tetrathioalkyl-phenalenyl)boron radicals

We report the synthesis and properties of a series of spiro-bis(3,4,6,7-tetrachalcogenide-substituted-phenalenyl)boron salts and two of the corresponding tetrathioalkyl-substituted spiro-bis(phenalenyl)boron radicals [tetrathiomethyl (10) and tetrathioethyl (11)] in which all of the active positions of the phenalenyl (PLY) nucleus are functionalized. In the solid state, radicals 10 and 11 exist as a weak π-dimers due to the steric congestion of the thioalkyl groups in the superimposed PLY units. As a result, the spins are localized in the isolated (nonsuperimposed) PLY rings, and the structure, magnetic susceptibility measurements, and band structure calculations confirm that these PLY units are unable to undergo strong intermolecular interaction as a result of the orientation of the thioalkyl groups.

Experimental and computational study on the structure and properties of Herz cations and radicals: 1,2,3-benzodithiazolium, 1,2,3-benzodithiazolyl, and their Se congeners

Salts of 1,2,3-benzodithiazolium (1), 2,1,3-benzothiaselenazolium (3), and 1,2,3-benzodiselenazolium (4) (Herz cations), namely, [1][BF4], [1][SbCl6], [3][BF4], [3][GaCl4], [3][SbCl6], and [4][GaCl4], were prepared from the corresponding chlorides and NaBF4, GaCl3, or SbCl5. It was found that [1][SbCl6] and [3][SbCl6] spontaneously transform in MeCN solution to [1]3[SbCl6]2[Cl] and [3]3[SbCl6]2[Cl], respectively. [1][BF4], [1]3[SbCl6]2[Cl], [3][BF4], [3]3[SbCl6]2[Cl], and [4][GaCl4] were structurally characterized by X-ray diffraction (XRD). In solution, these [BF4](-) and [GaCl4](-) salts as well as [1][GaCl4], [2][GaCl4], [3][GaCl4], [3][Cl], and [4][Cl] were characterized by multinuclear nuclear magnetic resonance (NMR). The corresponding Herz radicals 1(•)-4(•) were obtained in toluene and DCM solutions by the reduction of the appropriate salts with Ph3Sb and characterized by EPR. Cations 1-4 and radicals 1(•)-4(•) were investigated computationally at the density functional theory (DFT) and second-order Møller-Plesset (MP2) levels of theory. The B1B95/cc-pVTZ method was found to satisfactorily reproduce the experimental geometries of 1-4; an increase in the basis set size to cc-pVQZ results in only minor changes. For both 1-4 and 1(•)-4(•), the Hirshfeld charges and bond orders, as well as the Hirshfeld spin densities for the radicals, were calculated using the B1B95/cc-pVQZ method. It was found for both the cations and the radicals that replacing S atoms with Se atoms leads to considerable changes in the atomic charges, bond lengths, and bond orders only at the involved and the neighboring sites. According to the calculations, 60% of the positive charge in the cations and 80% of the spin density in the radicals is localized on the heterocycles, with the spin density distributions being very similar for all radicals 1(•)-4(•). For the cations 1-4, the NICS values (B3LYP/cc-pVTZ for B1B95/cc-pVTZ geometries) lie in the narrow range from -5.5 ppm to -6.6 ppm for the carbocycles, and from -14.4 ppm to -15.5 ppm for heterocycles, clearly indicating the aromaticity of the cations. Calculations on radical dimers [1(•)]2-[4(•)]2 revealed, with only one exception, positive dimerization energies, i.e., the dimers are inherently unstable in the gas phase.