Origin of bistability in the butyl-substituted spirobiphenalenyl-based neutral radical material

One of the most remarkable bistable materials reported so far is made of pi dimers of a butyl-substituted spirobiphenalenyl boron radical (butyl-SBP). The phase transition of this material, which is accompanied by changes in its optical, conductive, and magnetic properties, occurs with a hysteretic loop 25 K wide centered at 335 K. Herein, a computational study is presented aimed at deciphering the origin of this hysteresis. The phase transition of butyl-SBP consists of a spin transition of the constituent pi dimers coupled with an order-disorder transition involving the butyl chains linked to the nitrogen atoms of the superimposed phenalenyl rings of the pi dimer. Below 335 K, the terminal methyl group of the butyl chains adopts a gauche conformation with respect to the methylene unit bonded to the nitrogen atom. Above 335 K, the methyl group is in an anti conformation and exhibits dynamic disorder. The gauche -> anti conformational rearrangement triggers the spin transition of the pi dimers and is responsible for the hysteretic behavior of butyl-SBP. Specifically, the onset of the phase transition in the heating mode, and thus, the width of the hysteresis loop, are governed by the high energy cost and strong structural cooperative effects associated with this conformational change. Our results show that coupling a spin switch with a conformational switch in a molecular crystal provides a promising strategy in the design of new bistable materials

A first-principles analysis of the magnetism of CuII polynuclear coordination complexes: the case of [Cu4(bpy)4(aspartate)2(H2O)3](ClO4)4·2.5H2O

The magnetic structure of the [Cu4(bpy)4(aspartate)2(H2O)3](ClO4)4·2.5 H2Ocrystal - using fractional coordinates determined at room-temperature ¿ has beenanalysed in detail. This analysis has been carried out by extending our first principlesbottom-up theoretical approach, which was initially designed to study through-spacemagnetic interactions, to handle through-bond magnetic interactions. The only input datarequired by this approach are the values of the computed JAB exchange parameters for allthe unique pairs of spin-containing centres. The results allow the magnetic structure ofthe crystal, which presents two types of isolated tetranuclear CuII clusters, to be definedin quantitative terms. Each of these clusters presents ferro and antiferromagneticinteractions, the former being stronger, although outnumbered by the latter. Thecomputed magnetic susceptibility curve shows the same qualitative features as theexperimental data. However, there are small differences that are presumed to beassociated with the use of room-temperature crystal coordinates

Interplay between the Gentlest Ascent Dynamics Method and Conjugate Directions to Locate Transition States

An algorithm to locate transition states on a potential energy surface (PES) is proposed and described. The technique is based on the GAD method where the gradient of the PES is projected into a given direction and also perpendicular to it. In the proposed method, named GAD-CD, the projection is not only applied to the gradient but also to the Hessian matrix. Then, the resulting Hessian matrix is block diagonal. The direction is updated according to the GAD method. Furthermore, to ensure stability and to avoid a high computational cost, a trust region technique is incorporated and the Hessian matrix is updated at each iteration. The performance of the algorithm in comparison with the standard ascent dynamics is discussed for a simple two dimensional model PES. Its efficiency for describing the reaction mechanisms involving small and medium size molecular systems is demonstrated for five molecular systems of interest

Towards the tailored design of benzotriazinyl-based organic radicals displaying a spin transition

The mechanism of the phase transition of 1-phenyl-3-trifluoromethyl- 1,4-dihydrobenzo[e][1,2,4]triazin-4-yl (1), the first reported triazinyl radical to present such a feature, is unveiled. In so doing, we identify the key ingredients that are crucial to enable the phase transition in this family of radicals, and how those can be exploited by a rational design of the spin-carrying units

Structural, ferroelectric, and optical properties of Bi3+ doped YFeO3: A first‐principles study

The orthoferrites with the general formula RFeO3 (R = Ho, Er, Lu, Sc, and Y) have recently attracted a great deal of attention because they are promising candidates for a second generation of multiferroic materials. In this computational work, the structural, ferroelectric and optical properties of the YFeO3 perovskite oxide (YFO) and a Bi‐doped YFeO3 were analyzed. Bi‐substitution in YFO leads to an increase of its lattice parameters by virtue of the larger ionic radius of Bi3+. Both compounds exhibit a G‐type antiferromagnetic ground state. The calculations disclose a significant spontaneous polarization along the [101] direction of YFO‐Bi, which originates in the asymmetric distribution of the charges around the Bi3+ ions, as a result of the Bi‐6s electrons. The electric polarizability of YFO is increased upon Bi3+‐doping and the more significant components of the real permittivity tensor of YFO‐Bi are those associated with the direction along which the maximum value of spontaneous polarization is observed. The spontaneous polarization of YFO‐Bi found in this work reveals that this compound holds the potential for the next generation of multi ferroic materials

Tracing the sources of the different magnetic behavior in the two phases of the bistable (BDTA)2[Co(mnt)2] compound

A complete computational study of the magnetic properties of the two known phases of the bistable (BDTA)2[Co(mnt)2] compound is presented. The origin of their different magnetic properties can be traced to a variation in the values of the g tensor, together with a hitherto unknown change in the JAB values and their magnetic topology

Formation of long, multicenter pi-[TCNE](2)(2-) dimers in solution: solvation and stability assessed through molecular dynamics simulations

Purely organic radical ions dimerize in solution at low temperature, forming long, multicenter bonds, despite the metastability of the isolated dimers. Here, we present the first computational study of these pi-dimers in solution, with explicit consideration of solvent molecules and finite temperature effects. By means of force-field and ab initio molecular dynamics and free energy simulations, the structure and stability of pi-[TCNE](2)(2-) (TCNE = tetracyanoethylene) dimers in dichloromethane have been evaluated. Although the dimers dissociate at room temperature, they are stable at 175 K and their structure is similar to the one in the solid state, with a cofacial arrangement of the radicals at an inter-planar separation of approximately 3.0 angstrom. The pi-[TCNE](2)(2-) dimers form dissociated ion pairs with the NBu4+ counterions, and their first solvation shell comprises approximately 20 CH2Cl2 molecules. Among them, the eight molecules distributed along the equatorial plane of the dimer play a key role in stabilizing the dimer through bridging C-H center dot center dot center dot N contacts. The calculated free energy of dimerization of TCNE center dot- in solution at 175 K is -5.5 kcal mol(-1). These results provide the first quantitative model describing the pairing of radical ions in solution, and demonstrate the key role of solvation forces on the dimerization process

Effect of La3+/Sr2+ ordering on the magnetic properties of La2/3Sr1/3MnO3 by first principles calculations

In this work, using DFT + U formalism, we investigate the effect of order-disorder in the A-site occupation byLa3+and Sr2+on the stability of the ferromagnetic order in La2/3Sr1/3MnO3with−Rc3symmetry. To date, adetailed theoretical discussion of such phenomenon, using a combination of different representations of theelectronic structure, is still missing in the Literature. We employed structural models consisting of 120 atomsupercells constructed according to the precise stoichiometry of the compound. Two configurations, describingrandomized and ordered occupation of the La3+/Sr2+ions, were evaluated. We demonstrate that the ferro-magnetic arrangement of La2/3Sr1/3MnO3with randomly distributed La3+and Sr2+ions is more stable. In suchconfiguration wefind that the Mn3+and Mn4+ions are not distinguished, favoring the double-exchange me-chanism, enhanced by the higher degree of covalence in the MneO bonds near the Fermi level between thespin-upMn-egorbitals and the O-porbitals

Insights into the magnetism and phase transitions of organic radical-based materials

Organic radicals have been consistently regarded as promising building blocks for the next generation of applied materials. Multiple radical families have been developed and characterized in the last decades, fostered by the ever-growing capabilities of organic synthesis. Thiazyl-, spiro-biphenalenyl-, 1,2,4-benzotriazinyl-, and nitroxide-based radicals have furnished striking examples of metal-free switchable materials, whose phase transitions are accompanied by changes in the magnetic, optical and/or electrical response. While similar in origin, their actual mechanism, driving force(s), and spin state stabilities often depict a different landscape. Fundamental knowledge on such aspects, as well as on the underlying network of spin exchange couplings and non-covalent interactions (including pancake bonding), are key to understand their spin transition, and the tailored modification of their properties. These complex features cannot be extracted based solely on experimental input, but rely on a computational interpretation. In this Perspective article, we discuss the insight gained from computational modelling into the magnetism and phase transitions of organic radical-based materials. We focus on the key importance of dynamic effects due to the labile nature of π-stack interactions assembling those materials, the structural distortions driven by spin changes, the coupling between electronic structure and order-disorder transitions, and the dependence of spin correlation upon temperature. All these phenomena uncovered by simulations should assist in the rational design of new dynamic organic crystals

Bistability in organic magnetic materials: a comparative study of the key differences between hysteretic and non-hysteretic spin transitions in dithiazolyl radicals

Dithiazolyl (DTA)- based radicals have furnished many examples of organic spin- transition materials, some of them occurring with hysteresis and some others without. Herein, we present a combined computational and experimental study aimed at deciphering the factors controlling the existence or absence of hysteresis by comparing the phase transitions of 4- cyanobenzo- 1,3,2- dithiazolyl and 1,3,5trithia- 2,4,6- triazapentalenyl radicals, which are prototypical examples of non- bistable and bistable spin transitions, respectively. Both materials present low- temperature diamagnetic and high- temperature paramagnetic structures, characterized by dimerized (center dot center dot center dot A-A center dot center dot center dot A-A center dot center dot center dot) n and regular (center dot center dot center dot A center dot center dot center dot A center dot center dot center dot A center dot center dot center dot A center dot center dot center dot) n pi- stacks of radicals, respectively. We show that the regular pi-stacks are not potential energy minima but average structures arising from a dynamic inter-conversion between two degenerate dimerized configurations: (center dot center dot center dot A-A center dot center dot center dot A-A center dot center dot center dot)n (-A center dot center dot center dot A-A center dot center dot center dot A-)n. The emergence of this intrastack dynamics upon heating gives rise to a second-order phase transition that is responsible for the change in the dominant magnetic interactions of the system. This suggests that the promotion of a (center dot center dot center dot A-A center dot center dot center dot A-A center dot center dot center dot)n (-A center dot center dot center dot A-A center dot center dot center dot A-) n dynamics is a general mechanism for triggering spin transitions in DTA-based materials. Yet, this intra-stack dynamics does not suffice to generate bistability, which also requires a rearrangement of the intermolecular bonds between the pstacks via a first-order phase transition