scispace - formally typeset
Search or ask a question
Journal ArticleDOI

Origin of Bistability in the Butyl-Substituted Spirobiphenalenyl-Based Neutral Radical Material

07 Jun 2017-Chemistry: A European Journal (John Wiley & Sons, Ltd)-Vol. 23, Iss: 32, pp 7772-7784
TL;DR: The 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.
Abstract: One of the most remarkable bistable materials reported so far is made of π 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 π dimers coupled with an order-disorder transition involving the butyl chains linked to the nitrogen atoms of the superimposed phenalenyl rings of the π 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 π 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.

Summary (2 min read)

Introduction

  • Therefore, the studies aimed at elucidating the origin of these barriers and at establishing the role of cooperative effects have the potential to offer most valuable hints on how to devise new bistable materials with improved properties.
  • Below the spin transition temperature, the structures of the πdimers are governed by the potential energy surface (PES) of the ground singlet state (1Ag state), whose minimum structure features a partial localization of the unpaired electrons of each SBP radical in the superimposed phenalenyl (sup-PLY) rings, that is, on the phenalenyl (PLY) units directly involved in the π-dimer .

1) Phase transition of butyl-SBP: a spin transition coupled with a conformational

  • The optimized structures of the LT and HT polymorphs were obtained by means of variable-cell geometry relaxations, in which the atomic positions and the lattice parameters are optimized simultaneously.
  • A constant number of plane waves imply no Pulay stress but a decreasing precision of the calculation as the volume of the cell increases.
  • The large cutoff employed in these calculations ensures that the artifacts arising from this change of precision are negligible.
  • The starting atomic positions and initial lattice parameters for the relaxation of the LS(gau) and HS (anti) polymorphs were taken from the X-ray resolved structures of the LT and HT phases of 2 at 100 and 360 K, respectively.
  • The optimizations of the isolated π-dimers of 2 (carried out with the goal of evaluating the gas-phase ΔEadia values) were also done with plane wave pseudopotential calculations using Vanderbilt ultrasoft pseudopotentials.

3) Driving forces of the phase transition of butyl-SBP. Order-disorder transition

  • The vibrational entropy of the different polymorphs and the isolated π-dimers was evaluated after computing the vibrational frequencies of these systems in the harmonic approximation.
  • The vibrational frequencies in the condensed phase were calculated by means of a finite-difference normal-mode analysis of the optimized structures.
  • It thus follows that the computational strategy has been properly set up.
  • All dynamic simulations were performed in the canonical ensemble using   23   the Nosé-Hoover chain thermostats100 in order to control the kinetic energy of the nuclei and the fictitious kinetic energy of the orbitals.
  • In the simulations of the LT phase, the electronic structure of the π-dimers was that corresponding to their singlet ground state.

4) Origin of the hysteresis and the coupling between the spin transition and the

  • Conformational change in the phase transition of butyl-SBP For the LS(gau) polymorph, in turn, five calculations were performed along the same rotation coordinate between the -107° and 0° values of the θ dihedral angle.
  • To study the elementary steps of the LS(gau) à LS(anti) phase transition, the intermediate states connecting these two polymorphs (3gauche-1anti, 2gauche-2anti and 1gauche-3anti) were obtained by means of variable-cell geometry relaxations, in which the atomic positions and the lattice parameters are optimized simultaneously.
  • Cell parameters of the reported LT-340 and HT-340 X-ray crystal structures of 2. Table S3.
  • Cell parameters of the LT-0 minimum energy structure of butyl-SBP and of the LT crystallographic structures resolved at 100 and 340 K.

Acknowledgments

  • 83 The ΔEadia value between both gauche-IN LS and HS minima is ca. 1 kcal/mol larger than the ΔEadia value between anti polymorphs .
  • Therefore, the large increase in vibrational entropy associated with the conformational switch explains why the system undergo first the conformational switch even if the energy gap that is cleared in this transition is larger than that of the spin switch.

Did you find this useful? Give us your feedback

Figures (11)

Content maybe subject to copyright    Report

!
1!
The origin of bistability in the butyl-substituted spiro-
biphenalenyl-based neutral radical material
!
Maria Fumanal
, Juan J. Novoa, Jordi Ribas-Arino*
Departament de Química Física and IQTCUB, Facultat de Química, Universitat
de Barcelona, Av. Diagonal 645, 08028-Barcelona (Spain)
* jordi.ribas.jr@gmail.com, j.ribas@ub.edu
Present address: Laboratoire de Chimie Quantique, Institut de Chimie UMR7177
CNRS-Université de Strasbourg, 1 Rue Blaise Pascal BP 296/R8, F-67007 Strasbourg,
France
!
!

!
2!
!
Abstract
One of the most remarkable bistable materials so far reported is made of π-dimers of a
butyl-substituted spiro-biphenalenyl 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 and is centered at 335 K. Here, we
present a computational study aimed at deciphering the origin of this hysteresis. We
show that the phase transition of butyl-SBP consists of a spin transition of their
constituent π-dimers coupled with an order-disorder transition involving the butyl
chains linked to the N atoms of the superimposed phenalenyl rings of the π-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 N 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 π-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 the 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.
!

!
3!
!
Introduction
Bistability is an intriguing phenomenon exhibited by a few materials that present two
stable phases that can both exist within a given range of temperatures. Molecule-
based bistable materials have been the subject of intense research during the last
years because they hold great promise for application in sensors, displays and
switching devices.
1,2,3 ,4 ,5
The numerous examples of molecular bistable materials
include: materials based on transition metal complexes undergoing spin
transitions
6,7,8,9,10,11,12,13,14
, organic spin-transition materials
15,16,17,18,19,20,21,22
, compounds
whose phase transition is induced by a charge transfer between an electron-donor and
an electron-acceptor
23, 24,25 ,26 ,27
, compounds featuring charge-transfer-induced spin
transitions
28 , 29 , 30
, inorganic-organic hybrid frameworks undergoing phase
transitions
31,32
, molecular crystals whose phase transitions are triggered by changes in
the orientation of molecules
33
. The transition temperature and the hysteresis loop width
of a bistable material are crucial parameters in determining whether its bistability can
be harnessed in technological applications. These two parameters, in turn, depend on
the intermolecular interactions within the crystal and on the energy barriers associated
with the lattice reorganization upon phase transition. For most of the bistable
compounds reported so far, very little is known about either the origin of the energy
barriers associated with their phase transitions (i.e, whether the energy barrier of the
overall phase transition is dominated by a single molecular rearrangement or whether
the barrier is the result of the contributions of different reorganization events) nor the
role of structural cooperativity in promoting such phase transitions. It is clear that the
lack of this sort of information poses a major obstacle for the rational design of new
derivatives of a given bistable parent compound with the goal of fine tuning its
transition temperature and its hysteresis loop width. Therefore, the studies aimed at
elucidating the origin of these barriers and at establishing the role of cooperative
effects have the potential to offer most valuable hints on how to devise new bistable
materials with improved properties. Here, on the basis of a computational study, we
disclose the origin of the hysteretic phase transition of a phenalenyl-based butyl-
substituted neutral radical, which is one of the most prominent compounds within the
family of bistable materials.
Phenalenyl (PLY) is an odd-alternant hydrocarbon neutral radical arising from a
triangular fusion of three benzene rings. This open-shell molecule has emerged in the
past years as one of the most versatile building blocks for functional molecular devices

!
4!
and materials.
16,34,35,36,37,38,39,40
The numerous spiro-biphenalenyl (SBP) boron radicals
reported by Haddon and coworkers constitute a very important class of PLY
derivatives.
41,42,43,44,45,46,47,48,49,50,51,52,53,54,55,56,57
SBPs present two nearly perpendicular
phenalenyl units connected through a boron spiro-linkage. The N- and O-functionalized
SBPs (ie, SBPs in which each phenalenyl unit is bonded to the central boron atom via
an oxygen and a nitrogen atom) exhibit diverse packing motifs in the solid state, and
hence different physical properties, depending on the substituents attached to the
nitrogen atom. Ethyl (1) and butyl-substituted (2) SBPs (see Figure 1) present a crystal
structure containing π -dimers as the basic building block (see Figure 2 and Figure S1).
These two compounds undergo a phase transition that is accompanied by a change in
their optical, conductive and magnetic properties.
16,42
The phase transition of ethyl-SBP
is reversible and occurs at about 140 K, while that of butyl-SBP occurs with an
hysteretic loop 25-K wide and is centered at a much higher temperature (~ 335 K). At
this point, it is worth mentioning that butyl-SBP is one of the few multifunctional bistable
materials that switch the response in multiple physical channels upon phase
transition.
25,30,32,26
Besides, the volume of the crystals of butyl-SBP significantly change
upon phase transition; specifically, a notable expansion of the crystal is observed when
the system switches from its low-temperature (LT) phase to its high-temperature (HT)
phase.
58
This volume change in response to external stimuli is currently a sought-after
phenomenon in the context of new functional materials due to its potential applicability
to microscale or nanoscale actuators.
33
The experimental
58,59
and theoretical studies
60,61,62,63,64
conducted over the last years on
ethyl- and butyl-SBP have culminated in a clear understanding of their electronic
structure and the different magnetic and conducting properties of their phases. Upon
phase transition in the heating mode, the constituent π-dimers of these materials
undergo a spin transition from a closed-shell diamagnetic singlet state to an open-shell
paramagnetic state. Below the spin transition temperature, the structures of the π-
dimers are governed by the potential energy surface (PES) of the ground singlet state
(
1
A
g
state), whose minimum structure features a partial localization of the unpaired
electrons of each SBP radical in the superimposed phenalenyl (sup-PLY) rings, that is,
on the phenalenyl (PLY) units directly involved in the π-dimer (see Figure 2a). The
strong coupling between the SBP unpaired electrons in this configuration leads to a
magnetically silent state, and, thus, to a diamagnetic LT phase. Above the spin
transition temperature, the π-dimers adopt a configuration characterized by a
localization of the SBP unpaired electrons in the nonsuperimposed phenalenyl (non-
PLY) units, that is, on the PLYs not directly involved in the π-dimer (see Figure 2b),

!
5!
which leads to a paramagnetic phase. This configuration is exclusively governed by the
PES of the ground triplet state (
3
A
u
state) because the corresponding open-shell singlet
does not feature any minimum in that region of the PES even if it lies slightly below in
energy than the triplet state. In a recent article
64
, we have shown that the high-spin
(HS) state is energetically competitive with the low-spin (LS) state because the
electrostatic component of the interaction energy between SBP radicals in the π-
dimers is more attractive in the high-temperature
3
A
u
state than in the low-temperature
1
A
g
state. This electrostatic stabilization of the high-temperature
3
A
u
state was ascribed
to the zwitterionic nature of the SBP moieties, in particular, to the interaction between
the positively-charged superimposed PLYs in the triplet state (Figure 2b) and the
negatively-charged spiro-linkages with the central boron atom. These electrostatic
interactions also explain why the unpaired electrons prefer to localize on the
nonsuperimposed PLYs in the high-temperature triplet state.
64
Despite the current good understanding of the electronic structure of the π-dimers of
ethyl- and butyl-SBP and several theoretical studies on other phenalenyl-based
systems
65,66,67,68,69,70,71,72,73,74,75,76,77,78
, there are two crucial questions concerning the
phase transitions of ethyl- and butyl-SBP that remain unsettled, namely: i) why is the
transition temperature of butyl-SBP so much higher than that of ethyl-SBP?, and ii) why
does butyl-SBP display an hysteretic phase transition, in contrast with ethyl-SBP,
which features a smooth phase transition? A meticulous study carried out by Haddon
and coworkers in Ref. 58 on numerous crystal structures of butyl-SBP at different
temperatures led to the suggestion that the HT phase is the thermodynamically stable
phase within the bistability region, while the existence of the LT phase within the
hysteretic loop was rationalized on the basis of the large energy barrier that the system
needs to overcome when switching from LT to HT. Even if this barrier was estimated to
be larger than 24 kcal/mol, the specific molecular rearrangements responsible for that
barrier were not identified. In the computational study herein presented, not only do we
provide a rationale for the higher spin-transition temperature of butyl-SBP but also
disclose the hitherto elusive origin of its hysteresis loop. In particular, our study reveals
that the bistability arises from a very simple molecular rearrangement, namely, a
conformational rearrangement of the butyl groups attached to the SBP radicals.

Citations
More filters
Journal ArticleDOI
TL;DR: In this paper, the authors describe the fabrication of bench-stable radical ions under ambient conditions, which is of utmost importance from the perspective of materials and structural (solid-state) chemistry.
Abstract: Fabrication of bench-stable radical ions under ambient conditions is of utmost significance from the perspective of materials and structural (solid-state) chemistry. Two exceptionally stable benzot...

8 citations

Journal ArticleDOI
TL;DR: In this article, the authors reveal that regular π-stacks of planar dithiazolyl radicals can be rendered stable, in some range of temperatures, via two different mechanisms.
Abstract: Materials based on regular π-stacks of planar organic radicals are intensively pursued by virtue of their technologically relevant properties. Yet, these π-stacks are commonly unstable against π-dimerization. In this computational study, we reveal that regular π-stacks of planar dithiazolyl radicals can be rendered stable, in some range of temperatures, via two different mechanisms. When the radicals of a π-stack are both longitudinally and latitudinally slipped with respect to each other, the corresponding regular π-stacked configuration is associated with a locally stable minimum in the potential energy surface of the system. Conversely, those regular π-stacks in which radicals are latitudinally slipped with respect to each other are stable as a result of a dynamic interconversion between two degenerate dimerized configurations. The existence of two stabilization mechanisms, which can be traced back to the bonding properties of isolated π-dimers, translates into two different ways of exploiting spin-Peierls-like transitions in switchable dithiazolyl-based materials.

5 citations

Journal ArticleDOI
TL;DR: In this paper, an electrochromic phenalenyl based neutral radical small molecule was proposed to switch from transmissive to black states in both the visible and short-wave IR range.
Abstract: Applications for energy saving smart windows require materials which can switch from transmissive to black states in both the visible and short-wave IR range. We introduce an electrochromic phenalenyl based neutral radical small molecule and design devices capable of modulating light in both of these ranges.

4 citations

Journal ArticleDOI
TL;DR: substitutional doping was applied to [Bu]2 B by crystallizing solid solutions of bistable [Bu]:2 B and its non-radical-containing counterpart [Bu]; the abrupt transformation of the crystal structure to that of the pure [Bu?]2 Be crystal packing was observed, accompanied by a complete collapse of the hysteresis loop.
Abstract: The ability to tune the physical properties of bistable organic functional materials by means of chemistry can facilitate their development for molecular electronic switching components. The butylamine-containing biphenalenyl boron neutral radical, [Bu]2 B, crystalline compound has recently attracted significant attention by displaying a hysteretic phase transition accompanied by simultaneous bistability in magnetic, electrical, and optical properties close to room temperature. In this report, substitutional doping was applied to [Bu]2 B by crystallizing solid solutions of bistable [Bu]2 B and its non-radical-containing counterpart [Bu]2 Be. With increasing doping degree, the hysteretic phase transition is gradually suppressed in terms of reducing the height, but conserves the width of the hysteresis loop as observed through magnetic susceptibility and electrical conductivity measurements. At the critical doping level of about 6 %, the abrupt transformation of the crystal structure to that of the pure [Bu]2 Be crystal packing was observed, accompanied by a complete collapse of the hysteresis loop. Further study of the structure-properties relationships of bistable neutral radical conductors based on the [Bu]2 B host can be conducted utilizing a variety of biphenalenyl-based molecular conductors.

3 citations

Journal ArticleDOI
TL;DR: In this paper, the synthesis, structure and characterization of an Fe(II) complex, FeL2 (L = N′-(pyridin-2-ylmethylene)-4-(undec-10-en-1-yloxy)benzohydrazide), having long aliphatic chains with terminal CC bonds was reported.
Abstract: Spin-crossover (SCO) molecules have been proposed for applications in molecular spintronic devices due to their bistability which can be reversibly switched by an external stimulus. Here, we report the synthesis, structure and characterization of an Fe(II) complex, FeL2 (L = N′-(pyridin-2-ylmethylene)-4-(undec-10-en-1-yloxy)benzohydrazide), having long aliphatic chains with terminal CC bonds. Variable-temperature Single Crystal X-ray Diffraction (SCXRD) and magnetic susceptibility measurements reveal that the SCO behavior with a thermal hysteresis loop originates from the combination of the change of the Fe(II) center coordination environment and the thermal motion of long aliphatic chains. Differential Scanning Calorimetry (DSC), variable-temperature Powder X-ray Diffraction (PXRD) and Polarized Optical Microscopy (POM) measurements suggest that the latter is also associated with a crystal–liquid crystal (LC) phase transition at high temperatures. Together, these results provide the prerequisites for potential applications of bistable SCO molecules in spintronic devices and liquid crystal displays.

2 citations

References
More filters
Journal ArticleDOI
TL;DR: A simple derivation of a simple GGA is presented, in which all parameters (other than those in LSD) are fundamental constants, and only general features of the detailed construction underlying the Perdew-Wang 1991 (PW91) GGA are invoked.
Abstract: Generalized gradient approximations (GGA’s) for the exchange-correlation energy improve upon the local spin density (LSD) description of atoms, molecules, and solids. We present a simple derivation of a simple GGA, in which all parameters (other than those in LSD) are fundamental constants. Only general features of the detailed construction underlying the Perdew-Wang 1991 (PW91) GGA are invoked. Improvements over PW91 include an accurate description of the linear response of the uniform electron gas, correct behavior under uniform scaling, and a smoother potential. [S0031-9007(96)01479-2] PACS numbers: 71.15.Mb, 71.45.Gm Kohn-Sham density functional theory [1,2] is widely used for self-consistent-field electronic structure calculations of the ground-state properties of atoms, molecules, and solids. In this theory, only the exchange-correlation energy EXC › EX 1 EC as a functional of the electron spin densities n"srd and n#srd must be approximated. The most popular functionals have a form appropriate for slowly varying densities: the local spin density (LSD) approximation Z d 3 rn e unif

146,533 citations

Journal ArticleDOI
TL;DR: A new density functional of the generalized gradient approximation (GGA) type for general chemistry applications termed B97‐D is proposed, based on Becke's power‐series ansatz from 1997, and is explicitly parameterized by including damped atom‐pairwise dispersion corrections of the form C6 · R−6.
Abstract: A new density functional (DF) of the generalized gradient approximation (GGA) type for general chemistry applications termed B97-D is proposed. It is based on Becke's power-series ansatz from 1997 and is explicitly parameterized by including damped atom-pairwise dispersion corrections of the form C(6) x R(-6). A general computational scheme for the parameters used in this correction has been established and parameters for elements up to xenon and a scaling factor for the dispersion part for several common density functionals (BLYP, PBE, TPSS, B3LYP) are reported. The new functional is tested in comparison with other GGAs and the B3LYP hybrid functional on standard thermochemical benchmark sets, for 40 noncovalently bound complexes, including large stacked aromatic molecules and group II element clusters, and for the computation of molecular geometries. Further cross-validation tests were performed for organometallic reactions and other difficult problems for standard functionals. In summary, it is found that B97-D belongs to one of the most accurate general purpose GGAs, reaching, for example for the G97/2 set of heat of formations, a mean absolute deviation of only 3.8 kcal mol(-1). The performance for noncovalently bound systems including many pure van der Waals complexes is exceptionally good, reaching on the average CCSD(T) accuracy. The basic strategy in the development to restrict the density functional description to shorter electron correlation lengths scales and to describe situations with medium to large interatomic distances by damped C(6) x R(-6) terms seems to be very successful, as demonstrated for some notoriously difficult reactions. As an example, for the isomerization of larger branched to linear alkanes, B97-D is the only DF available that yields the right sign for the energy difference. From a practical point of view, the new functional seems to be quite robust and it is thus suggested as an efficient and accurate quantum chemical method for large systems where dispersion forces are of general importance.

23,058 citations

Journal ArticleDOI
TL;DR: QUANTUM ESPRESSO as discussed by the authors is an integrated suite of computer codes for electronic-structure calculations and materials modeling, based on density functional theory, plane waves, and pseudopotentials (norm-conserving, ultrasoft, and projector-augmented wave).
Abstract: QUANTUM ESPRESSO is an integrated suite of computer codes for electronic-structure calculations and materials modeling, based on density-functional theory, plane waves, and pseudopotentials (norm-conserving, ultrasoft, and projector-augmented wave). The acronym ESPRESSO stands for opEn Source Package for Research in Electronic Structure, Simulation, and Optimization. It is freely available to researchers around the world under the terms of the GNU General Public License. QUANTUM ESPRESSO builds upon newly-restructured electronic-structure codes that have been developed and tested by some of the original authors of novel electronic-structure algorithms and applied in the last twenty years by some of the leading materials modeling groups worldwide. Innovation and efficiency are still its main focus, with special attention paid to massively parallel architectures, and a great effort being devoted to user friendliness. QUANTUM ESPRESSO is evolving towards a distribution of independent and interoperable codes in the spirit of an open-source project, where researchers active in the field of electronic-structure calculations are encouraged to participate in the project by contributing their own codes or by implementing their own ideas into existing codes.

19,985 citations

Journal ArticleDOI
TL;DR: Novel features are that the pseudopotential itself becomes charge-state dependent, the usual norm-conservation constraint does not apply, and a generalized eigenproblem is introduced.
Abstract: A new approach to the construction of first-principles pseudopotentials is described. The method allows transferability to be improved systematically while holding the cutoff radius fixed, even for large cutoff radii. Novel features are that the pseudopotential itself becomes charge-state dependent, the usual norm-conservation constraint does not apply, and a generalized eigenproblem is introduced. The potentials have a separable form well suited for plane-wave solid-state calculations, and show promise for application to first-row and transition-metal systems.

18,782 citations

Journal ArticleDOI
TL;DR: Quantum ESPRESSO as discussed by the authors is an integrated suite of computer codes for electronic-structure calculations and materials modeling, based on density functional theory, plane waves, and pseudopotentials (norm-conserving, ultrasoft, and projector-augmented wave).
Abstract: Quantum ESPRESSO is an integrated suite of computer codes for electronic-structure calculations and materials modeling, based on density-functional theory, plane waves, and pseudopotentials (norm-conserving, ultrasoft, and projector-augmented wave). Quantum ESPRESSO stands for "opEn Source Package for Research in Electronic Structure, Simulation, and Optimization". It is freely available to researchers around the world under the terms of the GNU General Public License. Quantum ESPRESSO builds upon newly-restructured electronic-structure codes that have been developed and tested by some of the original authors of novel electronic-structure algorithms and applied in the last twenty years by some of the leading materials modeling groups worldwide. Innovation and efficiency are still its main focus, with special attention paid to massively-parallel architectures, and a great effort being devoted to user friendliness. Quantum ESPRESSO is evolving towards a distribution of independent and inter-operable codes in the spirit of an open-source project, where researchers active in the field of electronic-structure calculations are encouraged to participate in the project by contributing their own codes or by implementing their own ideas into existing codes.

13,052 citations