Enhancing Intramolecular Ferromagnetic Coupling in Tetrathiafulvalene-Nitronyl Nitroxide-Based Compounds through Spin Polarization Mechanism.

Spin-polarized donor radicals based on tetrathiafulvalene (TTF) derivatives and nitronyl nitroxide (NN) radicals in which one-electron oxidation involves the HOMO instead of the SOMO are well known for exhibiting magnetoresistance. In particular, BTBN consists of one dibromo-TTF and one NN radical, which are linked by a phenyl coupler group. One of the key factors driving magnetoresistance is the presence of intramolecular ferromagnetic (FM) coupling between the oxidized π-donor (TTF+⋅, D unit) and NN (R unit). Here, a theoretical study is carried out to assess suitable candidates with enhanced FM coupling with respect BTBN, which is thus used as a reference. The study is conducted via in silico chemical modification of the substituents of the BTBN basic functional units (D and R radicals, C coupler) to benefit from the spin polarization mechanism to boost the intramolecular FM coupling, aiming to distort the BTBN radical arrangement within the molecular crystal as little as possible, in the event the material can be synthesized. NICSiso(1) and Wiberg's Bond Order are analyzed to further assist in identifying promising potential candidates, since the decrease in aromaticity is expected to enhance the diradical character and give rise to a larger magnetic coupling value. The most favorable diradical building block to replace the BTBN moiety results from using a hydroxyl-ethylene (-(H)C=C(OH)-) as a coupler preserving BTBN original radicals, namely, NN and TTF+⋅ units. This study aims at illustrating the feasibility of improving the intramolecular FM interaction between radical moieties, which is fully realized, as a first step towards the synthesis of new materials with (possibly) enhanced magnetoresistance properties.

Optimal Oriented External Electric Fields to Trigger a Barrierless Oxaphosphetane Ring Opening Step of the Wittig Reaction.

The Wittig reaction is one of the most important processes in organic chemistry for the asymmetric synthesis of olefinic compounds. In view of the increasingly acknowledged potentiality of the electric fields in promoting reactions, here we will consider the effect of the oriented external electric field (OEEF) on the second step of Wittig reaction (i. e. the ring opening oxaphosphetane) in a model system for non-stabilized ylides. In particular, we have determined the optimal direction and strength of the electric field that should be applied to annihilate the reaction barrier of the ring opening through the polarizable molecular electric dipole (PMED) model that we have recently developed. We conclude that the application of the optimal external electric field for the oxaphosphetane ring opening favours a Bestmann-like mechanism.

An algorithm to find the optimal oriented external electrostatic field for annihilating a reaction barrier in a polarizable molecular system.

The use of oriented external electric fields (OEEFs) to promote and control chemical reactivity has motivated many theoretical and computational studies in the last decade to model the action of OEEFs on a molecular system and its effects on chemical processes. Given a reaction, a central goal in this research area is to predict the optimal OEEF (oOEEF) required to annihilate the reaction energy barrier with the smallest possible field strength. Here, we present a model rooted in catastrophe and optimum control theories that allows us to find the oOEEF for a given reaction valley in the potential energy surface (PES). In this model, the effective (or perturbed) PES of a polarizable molecular system is constructed by adding to the original, non-perturbed, PES a term accounting for the interaction of the OEEF with the intrinsic electric dipole and polarizability of the molecular system, so called the polarizable molecular electric dipole (PMED) model. We demonstrate that the oOEEF can be established by locating a point in the original PES with unique topological properties: the optimal barrier breakdown or bond-breaking point (oBBP). The essential feature of the oBBP structure is the fact that this point maintains its topological properties for all the applied OEEFs, also for the unperturbed PES, thus becoming much more relevant than the commonly used minima and transition state structures. The PMED model proposed here has been implemented in an open access package and is shown to successfully predict the oOEEF for two processes: an isomerization reaction of a cumulene derivative and the Huisgen cycloaddition reaction.

Fine-tuning of the spin-crossover properties of Fe(III) complexes via ligand design.

Exploring the chemical space of a given ligand aiming to modulate its ligand field strength is a versatile strategy for the fine-tuning of physical properties such as the transition temperature (T1/2) of spin-crossover (SCO) complexes. The computational study presented herein aims at systematically exploring the extent to which the ligand substituent effects can modulate T1/2 in two families of Fe(III) SCO systems with a N4O2 coordination environment and at identifying the best descriptors for fast and accurate prediction of changes in T1/2 upon ligand functionalization. B3LYP* calculations show that the attachment of substituents to ß-ketoiminato fragments (L1) leads to drastic changes in T1/2, while functionalization of phenolato moieties (L2) allows for a finer degree of control over T1/2. Natural Bond Orbital (NBO) charges of the donor atoms, Hammett parameters for both para and meta-functionalization of L2, and Swain-Lupton parameters for L1 and para-functionalization of L2 have been found to be the suitable descriptors for predicting the changes in T1/2. Further analysis of the ligand-field splitting in such systems rationalizes the observed trends and shows that ligand substituents modify both the σ and π bonds between the Fe(III) center and the ligands. Thus, we provide simple yet reliable guidelines for the rational design of new SCO systems with specific values of T1/2 based on their ligand design.

Emergent Spin Frustration in Neutral Mixed-Valence 2D Conjugated Polymers: A Potential Quantum Materials Platform.

Two-dimensional conjugated polymers (2DCPs)─organic 2D materials composed of arrays of carbon sp2 centers connected by π-conjugated linkers─are attracting increasing attention due to their potential applications in device technologies. This interest stems from the ability of 2DCPs to host a range of correlated electronic and magnetic states (e.g., Mott insulators). Substitution of all carbon sp2 centers in 2DCPs by nitrogen or boron results in diamagnetic insulating states. Partial substitution of C sp2 centers by B or N atoms has not yet been considered for extended 2DCPs but has been extensively studied in the analogous neutral mixed-valence molecular systems. Here, we employ accurate first-principles calculations to predict the electronic and magnetic properties of a new class of hexagonally connected neutral mixed-valence 2DCPs in which every other C sp2 nodal center is substituted by either a N or B atom. We show that these neutral mixed-valence 2DCPs significantly energetically favor a state with emergent superexchange-mediated antiferromagnetic (AFM) interactions between C-based spin-1/2 centers on a triangular sublattice. These AFM interactions are surprisingly strong and comparable to those in the parent compounds of cuprate superconductors. The rigid and covalently linked symmetric triangular AFM lattice in these materials thus provides a highly promising and robust basis for 2D spin frustration. As such, extended mixed-valence 2DCPs are a highly attractive platform for the future bottom-up realization of a new class of all-organic quantum materials, which could host exotic correlated electronic states (e.g., unusual magnetic ordering, quantum spin liquids).

The template effect of a SiF62- guest drives the formation of a heteroleptic Fe(II) coordination helicate.

The anion SiF62- exerts a strong template effect, driving the exclusive assembly of two different bispyridylpyrazolyl ligands into a triple stranded Fe(II) dinuclear heteroleptic helicate, engendering a new class within the large family of coordination helicates.

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.

Controlling Chemical Reactivity with Optimally Oriented Electric Fields: A Generalization of the Newton Trajectory Method.

The use of oriented external electric fields (OEEF) as a tool to accelerate chemical reactions has recently attracted much interest. A new model to calculate the optimal OEEF of the least intensity to induce a barrierless chemical reaction path is presented. A suitable ansatz is provided by defining an effective potential energy surface (PES), which considers the unperturbed or original PES of the molecular reactive system and the action of a constant OEEF on the overall dipole moment of system. Based on a generalization of the Newton Trajectories (NT) method, it is demonstrated that the optimal OEEF can be determined upon locating a special point of the potential energy surface (PES), the so-called "optimal bond-breaking point" (optimal BBP), for which two different algorithms are proposed. At this point, the gradient of the original or unperturbed PES is an eigenvector of zero eigenvalue of the Hessian matrix of the effective PES. A thorough discussion of the geometrical aspects of the optimal BBP and the optimal OEEF is provided using a two-dimensional model, and numerical calculations of the optimal OEEF for a SN2 reaction and the 1,3-dipolar retrocycloaddition of isoxazole to fulminic acid plus acetylene reaction serve as a proof of concept. The knowledge of the orientation of optimal OEEF provides a practical way to reduce the effective barrier of a given chemical process.

Accurate calculation of spin-state energy gaps in Fe(III) spin-crossover systems using density functional methods.

Fe(III) complexes are receiving ever-increasing attention as spin crossover (SCO) systems because they are usually air stable, as opposed to Fe(II) complexes, which are prone to oxidation. Here, we present the first systematic study exclusively devoted to assess the accuracy of several exchange-correlation functionals when it comes to predicting the energy gap between the high-spin (S = 5/2) and the low-spin (S = 1/2) states of Fe(III) complexes. Using a dataset of 24 different Fe(III) hexacoordinated complexes, it is demonstrated that the B3LYP* functional is an excellent choice not only for predicting spin-state energy gaps for Fe(III) complexes undergoing spin-transitions but also for discriminating Fe(III) complexes that are either low- or high-spin in the whole range of temperatures. Our benchmark study has led to the identification of a very versatile Fe(III) compound whose SCO properties can be engineered upon changing a single axial ligand. Overall, this work demonstrates that B3LYP* is a reliable functional for screening new spin-crossover systems with tailored properties.

Low temperature structures and magnetic interactions in the organic-based ferromagnetic and metamagnetic polymorphs of decamethylferrocenium 7,7,8,8-tetracyano-p-quinodimethanide, [FeCp*2]˙+[TCNQ]˙.

To identify the genesis of the differing magnetic behaviors for the ferro- (FO) and metamagnetic (MM) polymorphs of [FeCp*2][TCNQ] (Cp* = pentamethylcyclopentadienide; TCNQ = 7,7,8,8-tetracyano-p-quinodimethane) the low temperature (18 ± 1 K) structures of each polymorph were determined from high-resolution synchrotron powder diffraction data. Each polymorph possesses chains of alternating S = 1/2 [FeCp*2]˙+ cations and S = 1/2 [TCNQ]˙+, but with differing relative orientations. These as well as an additional paramagnetic polymorph do not thermally interconvert. In addition, the room and low (<70 ± 10 K) temperature structures of the MM polymorph, MMRT and MMLT, respectively, differ from that previously reported at 167 K (-106 °C) MM structure, and no evidence of either phase transition was previously noted even from the magnetic data. This transition temperature and enthalpy of this phase transition for MMRT⇌MM was determined to be 226.5 ± 0.4 K (-46.7 ± 0.4 °C) and 0.68 ± 0.04 kJ mol-1 upon warming, respectively, from differential calorimetry studies (DSC). All three MM phases are triclinic (P1[combining macron]) with the room temperature phase having a doubled unit cell relative to the other two. The lower temperature phase transition involves a small rearrangement of the molecular ions and shift in lattice parameters. These three MM and FO polymorphs have been characterized and form extended 1-D chains with alternating S = 1/2 [FeCp*2]˙+ cations, and S = 1/2 [TCNQ]˙- anions, whereas the fifth, paramagnetic (P) polymorph possesses S = 0 π-[TCNQ]22- dimers. At 18 ± 1 K the intrachain FeFe separations are 10.738(2) and 10.439(3) Å for the FO and MMLT polymorphs, respectively. The key structural differences between FO and MMLT at 18 ± 1 K are the 10% shorter interchain NN and the 2.8% shorter intrachain FeFe separation present for MMLT. Computational analysis of all nearest-neighbor spin couplings for the 18 K structures of FO and MMLT indicates that the intrachain [FeCp*2]˙+[TCNQ]˙- spin couplings (H = -2Si·Sj) are the strongest (4.95 and 6.5 cm-1 for FO and MMLT, respectively), as previously hypothesized, and are ferromagnetic due to their S = 1/2 spins residing in orthogonal orbitals. The change in relative [TCNQ]˙-[TCNQ]˙- orientations leads to a computed change from the ferromagnetic interaction (0.2 cm-1) for FO to an antiferromagnetic interaction (-0.1 cm-1) for MMLT in accord with its observed antiferromagnetic ground state. Hence, the magnetic ground state cannot be solely described by the dominant magnetic interactions.

Controlling pairing of π-conjugated electrons in 2D covalent organic radical frameworks via in-plane strain.

Controlling the electronic states of molecules is a fundamental challenge for future sub-nanoscale device technologies. π-conjugated bi-radicals are very attractive systems in this respect as they possess two energetically close, but optically and magnetically distinct, electronic states: the open-shell antiferromagnetic/paramagnetic and the closed-shell quinoidal diamagnetic states. While it has been shown that it is possible to statically induce one electronic ground state or the other by chemical design, the external dynamical control of these states in a rapid and reproducible manner still awaits experimental realization. Here, via quantum chemical calculations, we demonstrate that in-plane uniaxial strain of 2D covalently linked arrays of radical units leads to smooth and reversible conformational changes at the molecular scale that, in turn, induce robust transitions between the two kinds of electronic distributions. Our results pave a general route towards the external control, and thus technological exploitation, of molecular-scale electronic states in organic 2D materials.

Twistable dipolar aryl rings as electric field actuated conformational molecular switches.

The ability to control the chemical conformation of a system via external stimuli is a promising route for developing molecular switches. For eventual deployment as viable sub-nanoscale components that are compatible with current electronic device technology, conformational switching should be controllable by a local electric field (i.e. E-field gateable) and accompanied by a rapid and significant change in conductivity. In organic chemical systems the degree of π-conjugation is linked to the degree of electronic delocalisation, and thus largely determines the conductivity. Here, by means of accurate first principles calculations, we study the prototypical biphenyl based molecular system in which the dihedral angle between the two rings determines the degree of conjugation. In order to make this an E-field gateable system we create a net dipole by asymmetrically functionalising one ring with: (i) electron withdrawing (F, Br and CN), (ii) electron donating (NH2), and (iii) mixed (NH2/NO2) substituents. In this way, the application of an E-field interacts with the dipolar system to influence the dihedral angle, thus controlling the conjugation. For all considered substituents we consider a range of E-fields, and in each case extract conformational energy profiles. Using this data we obtain the minimum E-field required to induce a barrierless switching event for each system. We further extract the estimated switching speeds, the conformational probabilities at finite temperatures, and the effect of applied E-field on electronic structure. Consideration of these data allow us to assess which factors are most important in the design of efficient gateable electrical molecular switches.

Barnes Update Applied in the Gauss-Newton Method: An Improved Algorithm to Locate Bond Breaking Points.

A mechanochemical reaction is a reaction induced by mechanical energy. A general accepted model for this type of reaction consists of a first-order perturbation on the associated potential energy surface (PES) of the unperturbed molecular system due to mechanical stress or pulling force. Within this theoretical framework, the so-called optimal barrier breakdown points or optimal bond breaking points (BBPs) are critical points of the unperturbed PES where the Hessian matrix has a zero eigenvector that coincides with the gradient vector. Optimal BBPs are "catastrophe points" that are particularly important because their associated gradient indicates how to optimally harness tensile forces to induce reactions by transforming a chemical reaction into a barrierless process. Building on a previous method based on a nonlinear least-squares minimization to locate BBPs (Bofill et al., J. Chem. Phys. 2017, 147, 152710-10), we propose a new algorithm to locate BBPs of any molecular system based on the Gauss-Newton method combined with the Barnes update for a nonsymmetric Jacobian matrix, which is shown to be more appropriate than the Broyden update. The efficiency of the new method is demonstrated for a multidimensional model PES and two medium size molecular systems of interest in enzymatic catalysis and mechanochemistry.

Refined metadynamics through canonical sampling using time-invariant bias potential: A study of polyalcohol dehydration in hot acidic solutions.

We propose a canonical sampling method to refine metadynamics simulations a posteriori, where the hills obtained from metadynamics are used as a time-invariant bias potential. In this way, the statistical error in the computed reaction barriers is reduced by an efficient sampling of the collective variable space at the free energy level of interest. This simple approach could be useful particularly when two or more free energy barriers are to be compared among chemical reactions in different or competing conditions. The method was then applied to study the acid dependence of polyalcohol dehydration reactions in high-temperature aqueous solutions. It was found that the reaction proceeds consistently via an SN 2 mechanism, whereby the free energy of protonation of the hydroxyl group created as an intermediate is affected significantly by the acidic species. Although demonstration is shown for a specific problem, the computational method suggested herein could be generally used for simulations of complex reactions in the condensed phase.

Thermal spin crossover in Fe(ii) and Fe(iii). Accurate spin state energetics at the solid state.

The thermal spin crossover (SCO) phenomenon refers to an entropy-driven spin transition in some materials based on d6-d9 transition metal complexes. While its molecular origin is well known, intricate SCO behaviours are increasingly common, in which the spin transition occurs concomitantly to e.g. phase transformations, solvent absorption/desorption, or order-disorder processes. The computational modelling of such cases is challenging, as it requires accurate spin state energies in the solid state. Density Functional Theory (DFT) is the best framework, but most DFT functionals are unable to balance the spin state energies. While a few hybrid functionals perform better, they are still too expensive for solid-state minima searches in moderate-size systems. The best alternative is to dress cheap local (LDA) or semi-local (GGA) DFT functionals with a Hubbard-type correction (DFT+U). However, the parametrization of U is not straightforward due to the lack of reference values, and because ab initio parametrization methods perform poorly. Moreover, SCO complexes undergo notable structural changes upon transition, so intra- and inter-molecular interactions might play an important role in stabilizing either spin state. As a consequence, the U parameter depends strongly on the dispersion correction scheme that is used. In this paper, we parametrize U for nine reported SCO compounds (five based on FeII, 1-5 and four based on FeIII, 6-9) when using the D3 and D3-BJ dispersion corrections. We analyze the impact of the dispersion correction treatments on the SCO energetics, structure, and the unit cell dimensions. The average U values are different for each type of metal ion (FeIIvs. FeIII), and dispersion correction scheme (D3 vs. D3-BJ) but they all show excellent transferability, with mean absolute errors (MAE) below chemical accuracy (i.e. MAE <4 kJ mol-1). This enables a better description of SCO processes and, more generally, of spin state energetics, in materials containing FeII and FeIII ions.

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.

Selective Nanomechanics of Aromatic versus Aliphatic Thiolates on Gold Surfaces.

Thiolated gold nanointerfaces play a key role in numerous fields of science, technology, as well as modern medicine to coat, functionalize, and protect. Our computational study reveals that the mechanical vs thermal stabilities of aliphatic thiolates on gold surfaces are strikingly different from those of aromatic thiolates. The aliphatic thiolates feature, at the same time, a higher thermal desorption energy but a lower mechanical rupture force than thiophenolates. Our analysis discloses that this most counterintuitive property is due to different mechanochemical detachment mechanisms. Electronic structure analyses along the detachment pathways trace this back to the distinct electronic properties of the S─Au bond in stretched nanojunctions. The discoveries that it is a higher thermal stability that entails a lower mechanical stability and that mechanical loads generate different local nanostructures depending on the nature of the thiolate are highly relevant for the rational design of improved thiol-gold nanocontacts.

Understanding Competition of Polyalcohol Dehydration Reactions in Hot Water.

Dehydration of biomass-derived polyalcohols has recently drawn attention in green chemistry as a prototype of selective reactions controllable in hot water or hot carbonated water, without any use of organic solvents or metal catalysts. Here, we report a free-energy analysis based on first-principles metadynamics and blue-moon ensemble simulations to understand the mechanism of competing intramolecular dehydration reactions of 1,2,5-pentanetriol in hot acidic water. The simulations consistently predict that the most dominant mechanism is the proton-assisted SN2 process, where the protonation of the hydroxyl group by water and the C-O bond breaking and formation occur in a single step. However the free-energy barriers are different between the reaction paths: those leading to five-membered ether products, tetrahydrofurfuryl alcohol (THFA), are few kcal/mol lower than those leading to six-membered ether products, 3-hydroxytetrahydropyran (3-HTHP). A slight difference is seen in the timing of the protonation of the hydroxyl group of THFA and 3-HTHP on their reaction pathways. The detailed mechanism found from the simulations shows how the reaction paths are selective in hot water and why the reaction rates are accelerated in acidic environments, thus giving a clear explanation of experimental findings for a broad class of competing dehydration processes of polyalcohols.

The magnetic fingerprint of dithiazolyl-based molecule magnets.

Magnetic bistability in organic-radical based materials has attracted significant interest due to its potential application in electronic devices. The first-principles bottom-up study herein presented aims at elucidating the key factors behind the different magnetic response of the low and high temperature phases of four different switchable dithiazolyl (DTA)-based compounds. The drastic change in the magnetic response upon spin transition is always due to the changes in the JAB magnetic interactions between adjacent radicals along the π-stacks of the crystal, which in turn are driven mostly by the changes in the interplanar distance and degree of lateral slippage, according to the interpretation of a series of magneto-structural correlation maps. Furthermore, specific geometrical dispositions have been recognized as a ferromagnetic fingerprint in such correlations. Our results thus show that an appropriate substitution of the chemical skeleton attached to the DTA ring could give rise to new organic materials with dominant ferromagnetic interactions.

An algorithm to locate optimal bond breaking points on a potential energy surface for applications in mechanochemistry and catalysis.

The reaction path of a mechanically induced chemical transformation changes under stress. It is well established that the force-induced structural changes of minima and saddle points, i.e., the movement of the stationary points on the original or stress-free potential energy surface, can be described by a Newton Trajectory (NT). Given a reactive molecular system, a well-fitted pulling direction, and a sufficiently large value of the force, the minimum configuration of the reactant and the saddle point configuration of a transition state collapse at a point on the corresponding NT trajectory. This point is called barrier breakdown point or bond breaking point (BBP). The Hessian matrix at the BBP has a zero eigenvector which coincides with the gradient. It indicates which force (both in magnitude and direction) should be applied to the system to induce the reaction in a barrierless process. Within the manifold of BBPs, there exist optimal BBPs which indicate what is the optimal pulling direction and what is the minimal magnitude of the force to be applied for a given mechanochemical transformation. Since these special points are very important in the context of mechanochemistry and catalysis, it is crucial to develop efficient algorithms for their location. Here, we propose a Gauss-Newton algorithm that is based on the minimization of a positively defined function (the so-called σ-function). The behavior and efficiency of the new algorithm are shown for 2D test functions and for a real chemical example.