Toward High-Temperature Light-Induced Spin-State Trapping in Spin-Crossover Materials: The Interplay of Collective and Molecular Effects.

Spin-crossover (SCO) materials display many fascinating behaviors including collective phase transitions and spin-state switching controlled by external stimuli, e.g., light and electrical currents. As single-molecule switches, they have been fêted for numerous practical applications, but these remain largely unrealized-partly because of the difficulty of switching these materials at high temperatures. We introduce a semiempirical microscopic model of SCO materials combining crystal field theory with elastic intermolecular interactions. For realistic parameters, this model reproduces the key experimental results including thermally induced phase transitions, light-induced spin-state trapping (LIESST), and reverse-LIESST. Notably, we reproduce and explain the experimentally observed relationship between the critical temperature of the thermal transition, T1/2, and the highest temperature for which the trapped state is stable, TLIESST, and explain why increasing the stiffness of the coordination sphere increases TLIESST. We propose strategies to design SCO materials with higher TLIESST: optimizing the spin-orbit coupling via heavier atoms (particularly in the inner coordination sphere) and minimizing the enthalpy difference between the high-spin (HS) and low-spin (LS) states. However, the most dramatic increases arise from increasing the cooperativity of the spin-state transition by increasing the rigidity of the crystal. Increased crystal rigidity can also stabilize the HS state to low temperatures on thermal cycling yet leave the LS state stable at high temperatures following, for example, reverse-LIESST. We show that such highly cooperative systems offer a realistic route to robust room-temperature switching, demonstrate this in silico, and discuss material design rationale to realize this.

Multi-Redox Responsive Behavior in a Mixed-Valence Semiconducting Framework Based on Bis-[1,2,5]-thiadiazolo-tetracyanoquinodimethane.

The two-dimensional (2-D) framework, [Cu(BTDAT)(MeOH)] {BTDAT = bis-[1,2,5]-thiadiazolo-tetracyanoquinodimethane}, possesses remarkable multi-step redox properties, with electrochemical studies revealing six quasi-stable redox states in the solid state. In situ electron paramagnetic resonance and visible-near infrared spectroelectrochemistry elucidated the mechanism for these multi-step redox processes, as well as the optical and electrochromic behavior of the BTDAT ligand and framework. In studying the structural, spectroscopic, and electronic properties of [Cu(BTDAT)(MeOH)], the as-synthesized framework was found to exist in a mixed-valence state with thermally-activated semiconducting behavior. In addition to pressed pellet conductivity measurements, single-crystal conductivity measurements using a pre-patterned polydimethylsiloxane layer on a silicon substrate provide important insights into the anisotropic conduction pathways. As an avenue to further understand the electronic state of [Cu(BTDAT)(MeOH)], computational band structure calculations predicted delocalized electronic transport in the framework. On the balance of probabilities, we propose that [Cu(BTDAT)(MeOH)] is a Mott insulator (i.e., electron correlations cause a metal-insulator transition). This implies that the conductivity is incoherent. However, we are unable to distinguish between activated transport due to Coulombically bound electron-hole pairs and a hopping mechanism. The combined electrochemical, electronic, and optical properties of [Cu(BTDAT)(MeOH)] shine a new light on the experimental and theoretical challenges for electroactive framework materials, which are implicated as the basis of advanced optoelectronic and electrochromic devices.

Guest-Induced Multistep to Single-Step Spin-Crossover Switching in a 2-D Hofmann-Like Framework with an Amide-Appended Ligand.

A detailed study of the two-dimensional (2-D) Hofmann-like framework [Fe(furpy)2Pd(CN)4]·nG (furpy: N-(pyridin-4-yl)furan-2-carboxamide, G = H2O,EtOH (A·H2O,Et), and H2O (A·H2O)) is presented, including the structural and spin-crossover (SCO) implications of subtle guest modification. This 2-D framework is characterized by undulating Hofmann layers and an array of interlayer spacing environments─this is a strategic approach that we achieve by the inclusion of a ligand with multiple host-host and host-guest interaction sites. Variable-temperature magnetic susceptibility studies reveal an asymmetric multistep SCO for A·H2O,Et and an abrupt single-step SCO for A·H2O with an upshift in transition temperature of ∼75 K. Single-crystal analyses show a primitive orthorhombic symmetry for A·H2O,Et characterized by a unique FeII center─the multistep SCO character is attributed to local ligand orientation. Counterintuitively, A·H2O shows a triclinic symmetry with two inequivalent FeII centers that undergo a cooperative single-step high-spin (HS)-to-low-spin (LS) transition. We conduct detailed structure-function analyses to understand how the guest ethanol influences the delicate balance between framework communication and, therefore, the local structure and spin-state transition mechanism.

Regulation of Multistep Spin Crossover Across Multiple Stimuli in a 2-D Framework Material.

We investigate the effects of a broad array of external stimuli on the structural, spin-crossover (SCO) properties and nature of the elastic interaction within the two-dimensional Hofmann framework material [Fe(cintrz)2Pd(CN)4]·guest (cintrz = N-cinnamalidene 4-amino-1,2,4-triazole; A·guest; guest = 3H2O, 2H2O, and Ø). This framework exhibits a delicate balance between ferro- and antiferro-elastic interaction characters; we show that manipulation of the pore contents across guests = 3H2O, 2H2O, and Ø can be exploited to regulate this balance. In A·3H2O, the dominant antiferroelastic interaction character between neighboring FeII sites sees the low-temperature persistence of the mixed spin-state species {HS-LS} for {Fe1-Fe2} (HS = high spin, LS = low spin). Elastic interaction strain is responsible for stabilizing the {HS-LS} state and can be overcome by three mechanisms: (1) partial (2H2O) or complete (Ø) guest removal, (2) irradiation via the reverse light-induced excited spin-state trapping (LIESST) effect (λ = 830 nm), and (3) the application of external hydrostatic pressure. Combining experimental data with elastic models presents a clear interpretation that while guest molecules cause a negative chemical pressure, they also have consequences for the elastic interactions between metals beyond the simple chemical pressure picture typically proposed.

Hierarchical Spin-Crossover Cooperativity in Hybrid 1D Chains of FeII -1,2,4-Triazole Trimers Linked by [Au(CN)2 ]- Bridges.

Foremost, practical applications of spin-crossover (SCO) materials require control of the nature of the spin-state coupling. In existing SCO materials, there is a single, well-defined dimensionality relevant to the switching behavior. A new material, consisting of 1,2,4-triazole-based trimers coordinated into 1D chains by [Au(CN)2 ]- and spaced by anions and exchangeable guests, underwent SCO defined by elastic coupling across multiple dimensional hierarchies. Detailed structural, vibrational, and theoretical studies conclusively confirmed that intra-trimer coupling was an order of magnitude greater than the intramolecular coupling, which was an order of magnitude greater than intermolecular coupling. As such, a clear hierarchy on the nature of elastic coupling in SCO materials was ascertained for the first time, which is a necessary step for the technological development of molecular switching materials.

Quasi-one dimensional magnetic interactions in the three-dimensional hyper-honeycomb framework [(C2H5)3NH]2Cu2(C2O4)3.

The Cu(ii) ions in [(C2H5)3NH]2Cu2(C2O4)3 form a hyperhoneycomb lattice and show no indication of long-range magnetic order down to 60 mK. It has therefore been suggested that [(C2H5)3NH]2Cu2(C2O4)3 is a three dimensional quantum spin liquid. We construct a tight-binding model of [(C2H5)3NH]2Cu2(C2O4)3 from Wannier orbital overlaps. Including interactions within the Jahn-Teller distorted Cu-centered eg Wannier orbitals leads to a highly anisotropic effective Heisenberg model. We show that this anisotropy arrises from interference between different superexchange pathways. This demonstrates that when two (or more) orbitals contribute to the localised spin superexchange can be significantly richer than in the textbook single orbital case. The hyper-honeycomb lattice contains two symmetry distinct sublattices of Cu atoms arranged in coupled chains. We show that one sublattice is strongly dimerized, the other forms isotropic antiferromagnetic chains. Integrating out the strongest (intradimer) exchange interactions leaves extremely weakly coupled Heisenberg chains.

Quantitative calculations of the non-radiative rate of phosphorescent Ir(III) complexes.

It has recently been proposed that the dominant non-radiative decay mechanism in blue Ir(iii) phosphors at room temperature is due to the low-lying non-radiative metal-centred triplet states. These are populated thermally via an activated transition from the highly radiative metal-to-ligand-charge-transfer states that are initially populated due to intersystem crossing following the radiative or electronic excitation of the phosphor. We apply transition state theory to quantitatively calculate the non-radiative decay rate of a family of Ir(iii) complexes containing N-heterocyclic carbene (NHC) ligands. We compare the, computationally inexpensive, one-dimensional theory with the, more accurate, multi-dimensional theory. Both methods find a non-radiative rate with an Arrhenius form (knr = kae-ΔE/kBT). The pre-exponential factors, ka, and activation energies, ΔE, are evaluated via density functional theory (DFT). The multi-dimensional theory shows that there is an order of magnitude variation in ka within this family of materials (between 3 × 1011 s-1 and 3 × 1012 s-1). This is not captured by the one-dimensional theory, which predicts very uniform rate constants in the middle of this range (∼1012 s-1). Nevertheless, the activated process involved, and the linear relationship between ka and knr, mean that ka plays a subtle role in determining knr. Consistent with this we find that both methods capture the trend observed experimentally in the non-radiative rates. Furthermore, the magnitude of the calculated knr is similar in both methods and in good agreement with experimental values [except for one complex with a very shallow activation barrier (<0.1 eV)]. It has previously been demonstrated that radiative decay rates can be accurately calculated from DFT. Combined with our results for the non-radiative rates, this implies that DFT methods can accurately predict the emission efficiency in Ir(iii) phosphors. Therefore, DFT calculations are both fast and accurate enough to play a significant role in the design of new deep blue Ir(iii) phosphors with high emission efficiency. Even the one-dimensional theory provides reasonable agreement with experiment. This suggests that a funneling approach - where only the best performing molecules, according to the one-dimensional theory, are studied in the more laborious multi-dimensional framework - could be a powerful strategy for designing active materials for phosphorescent organic light-emitting diodes (PHOLEDs) from first principles.

Nonradiative Decay and Stability of N-Heterocyclic Carbene Iridium(III) Complexes.

Devices based on deep-blue emitting iridium(III) complexes with N-heterocyclic carbene (NHC) ligands have recently been shown to give excellent performance as phosphorescent organic light-emitting diodes (PHOLEDs). To facilitate the design of even better deep-blue phosphorescent emitters, we carried out density functional theory (DFT) calculations of the lowest triplet (T1) potential-energy surfaces upon lengthening the iridium-ligand (Ir-C) bonds. Relativistic time dependent DFT calculations demonstrate that this changes the nature of T1 from a highly emissive metal-to-ligand charge transfer (3MLCT) state to a metal centered (3MC) state where the radiative decay rate is orders of magnitude slower than that of the 3MLCT state. We identify the elongation of an Ir-C bond on the NHC group as the pathway with the lowest energy barrier between the 3MLCT and 3MC states for all complexes studied and show that the barrier height is correlated with the experimentally measured nonradiative decay rate. This suggests that the thermal population of 3MC states is the dominant nonradiative decay mechanism at room temperature. We show that the 3MLCT → 3MC transition is reversible, in marked contrast to deep-blue phosphors containing coordinating nitrogen atoms, where the population of 3MC states breaks Ir-N bonds. This suggests that, as well as improved efficiency, blue PHOLEDs containing phosphors where the metal is only coordinated by carbon atoms will have improved device lifetimes.

Effect of n-propyl substituents on the emission properties of blue phosphorescent iridium(iii) complexes.

Ligand substitution is often used for tuning the emission color of phosphorescent iridium(iii) complexes that are used in organic light-emitting diodes. However, in addition to tuning the emission color, the substituents can also affect the radiative and non-radiative decay rates of the excited state and hence the photoluminescence quantum yield. Understanding the substituent effect is therefore important for the design of new iridium(iii) complexes with specific emission properties. Using (time dependent) density functional methods, we investigate the substituent effect of n-propyl groups on the structure, emission color, and emission efficiency of fac-tris(1-methyl-5-phenyl-[1,2,4]triazolyl)iridium(iii) based phosphorescent complexes by comparing the calculated results for structural models with and without the n-propyl substituents. We find that attachment of the n-propyl groups increases the length of three Ir-N bonds, and although the emission color does not change significantly, the radiative and non-radiative rates do, leading to a prediction of enhanced blue phosphorescence emission efficiency. Furthermore, the calculations show that the attachment of the n-propyl groups leads to a larger activation energy to degradation and the formation of dark states.

Elucidating the Spatial Arrangement of Emitter Molecules in Organic Light-Emitting Diode Films.

The effect of varying the emitter concentration on the structural properties of an archetypal phosphorescent blend consisting of 4,4'-bis(N-carbazolyl)biphenyl and tris(2-phenylpyridyl)iridium(III) has been investigated using non-equilibrium molecular dynamics (MD) simulations that mimic the process of vacuum deposition. By comparison with reflectometry measurements, we show that the simulations provide an accurate model of the average density of such films. The emitter molecules were found not to be evenly distributed throughout film, but rather they can form networks that provide charge and/or energy migration pathways, even at emitter concentrations as low as ≈5 weight percent. At slightly higher concentrations, percolated networks form that span the entire system. While such networks would give improved charge transport, they could also lead to more non-radiative pathways for the emissive state and a resultant loss of efficiency.

Bond Fission and Non-Radiative Decay in Iridium(III) Complexes.

We investigate the role of metal-ligand bond fission in the nonradiative decay of excited states in iridium(III) complexes with applications in blue organic light-emitting diodes (OLEDs). We report density functional theory (DFT) calculations of the potential energy surfaces upon lengthening an iridium-nitrogen (Ir-N) bond. In all cases we find that for bond lengths comparable to those of the ground state the lowest energy state is a triplet with significant metal-to-ligand change transfer character ((3)MLCT). But, as the Ir-N bond is lengthened there is a sudden transition to a regime where the lowest excited state is a triplet with significant metal centered character ((3)MC). Time-dependent DFT relativistic calculations including spin-orbit coupling perturbatively show that the radiative decay rate from the (3)MC state is orders of magnitude slower than that from the (3)MLCT state. The calculated barrier height between the (3)MLCT and (3)MC regimes is clearly correlated with previously measured nonradiative decay rates, suggesting that thermal population of the (3)MC state is the dominant nonradiative decay process at ambient temperature. In particular, fluorination both drives the emission of these complexes to a deeper blue color and lowers the (3)MLCT-(3)MC barrier. If the Ir-N bond is shortened in the (3)MC state another N atom is pushed away from the Ir, resulting in the breaking of this bond, suggesting that once the Ir-N bond breaks the damage to the complex is permanent-this will have important implications for the lifetimes of devices using this type of complex as the active material. The consequences of these results for the design of more efficient blue phosphors for OLED applications are discussed.

Interplay of Zero-Field Splitting and Excited State Geometry Relaxation in fac-Ir(ppy)3.

The lowest energy triplet state, T1, of organometallic complexes based on iridium(III) is of fundamental interest, as the behavior of molecules in this state determines the suitability of the complex for use in many applications, e.g., organic light-emitting diodes. Previous characterization of T1 in fac-Ir(ppy)3 suggests that the trigonal symmetry of the complex is weakly broken in the excited state. Here we report relativistic time dependent density functional calculations of the zero-field splitting (ZFS) of fac-Ir(ppy)3 in the ground state (S0) and lowest energy triplet (T1) geometries and at intermediate geometries. We show that the energy scale of the geometry relaxation in the T1 state is large compared to the ZFS. Thus, the natural analysis of the ZFS and the radiative decay rates, based on the assumption that the structural distortion is a small perturbation, fails dramatically. In contrast, our calculations of these quantities are in good agreement with experiment.

Role of semiconductivity and ion transport in the electrical conduction of melanin.

Melanins are pigmentary macromolecules found throughout the biosphere that, in the 1970s, were discovered to conduct electricity and display bistable switching. Since then, it has been widely believed that melanins are naturally occurring amorphous organic semiconductors. Here, we report electrical conductivity, muon spin relaxation, and electron paramagnetic resonance measurements of melanin as the environmental humidity is varied. We show that hydration of melanin shifts the comproportionation equilibrium so as to dope electrons and protons into the system. This equilibrium defines the relative proportions of hydroxyquinone, semiquinone, and quinone species in the macromolecule. As such, the mechanism explains why melanin at neutral pH only conducts when "wet" and suggests that both carriers play a role in the conductivity. Understanding that melanin is an electronic-ionic hybrid conductor rather than an amorphous organic semiconductor opens exciting possibilities for bioelectronic applications such as ion-to-electron transduction given its biocompatibility.

Complex relaxation of trapped spin-states in spin crossover materials.

A diverse range of relaxation dynamics of trapped spin-states are observed in spin crossover (SCO) materials, including exponential, sigmoidal, stretched exponential, multi-step, and mixed kinetics. We reproduce and explain this full range of relaxation behaviours using a semi-empirical, semi-classical model that combines crystal field theory with elastic inter-molecular interactions. We show that frustrated intermolecular interactions, which are responsible for multistep thermal transitions, also lead to multiple energetically competitive ordered phases, even in systems that contain only one crystallographically distinct SCO site. This rugged free energy landscape leads to dynamic disorder and thence the complex dynamics widely observed in SCO materials. Similar mechanisms are vital for understanding dynamics of more complex materials from proteins to quantum materials.

Co-existence of five- and six-coordinate iron(II) species captured in a geometrically strained spin-crossover Hofmann framework.

Inclusion of an angular bridging ligand, 4,2':6',4''-terpyridine (TPy), into a Hofmann-type framework produces an irregular network in which six- and five-coordinate FeII species co-exist. The octahedral sites show thermally-induced spin-crossover (SCO) and the rare five-coordinate FeII sites are high-spin.

Dual-supramolecular contacts induce extreme Hofmann framework distortion and multi-stepped spin-crossover.

An extended nitro-functionalised 1,2,4-triazole ligand has been used to induce considerable lattice distortion in a 2-D Hofmann framework material via competing supramolecular interactions. Single crystal X-ray diffraction analyses on [Fe3(N-cintrz)6(Pd(CN)4)3]·6H2O (N-cintrz: (E)-3-(2-nitrophenyl)acrylaldehyde) reveal a substantial deviation from a regular Hofmann structure, in particular as the intra- and inter-layer contacts are dominated by hydrogen-bonding interactions rather than the typical π-stacking arrays. Also, the 2-D Hofmann layers show an assortment of ligand conformations and local FeII coordination environments driven by the optimisation of competing supramolecular contacts. Temperature-dependent magnetic susceptibility measurements reveal a two-step spin crossover (SCO) transition. Variable temperature structural analyses show that the two crystallographically distinct FeII centres, which are arranged in stripes (2 : 1 ratio) within each Hofmann layer, undergo a cooperative HS ↔ HS/LS ↔ LS (HS = high spin, LS = low spin) transition without periodic spin-state ordering. The mismatch between crystallographic (2 : 1) and spin-state (1 : 1) periodicity at the HS : LS step provides key insight into the competition (frustration) between elastic interactions and crystallographically driven order.

Hydration-controlled X-band EPR spectroscopy: a tool for unravelling the complexities of the solid-state free radical in eumelanin.

Melanin, the human skin pigment, is found everywhere in nature. Recently it has gained significant attention for its potential bioelectronic properties. However, there remain significant obstacles in realizing its electronic potential, in particular, the identity of the solid-state free radical in eumelanin, which has been implicated in charge transport. We have therefore undertaken a hydration-controlled continuous-wave electron paramagnetic resonance study on solid-state eumelanin. Herein we show that the EPR signal from solid-state eumelanin arises predominantly from a carbon-centered radical but with an additional semiquinone free radical component. Furthermore, the spin densities of both of these radicals can be manipulated using water and pH. In the case of the semiquinone radical, the comproportionation reaction governs the pH- and hydration-dependent behavior. In contrast, the mechanism underlying the carbon-centered radical's pH- and hydration-dependent behavior is not clear; consequently, we have proposed a new destacking model in which the intermolecular structure of melanin is disordered due to π-π destacking, brought about by the addition of water or increased pH, which increases the proportion of semiquinone radicals via the comproportionation reaction.