Excited-state dynamics of [Mn(im)(CO)3(phen)]+: PhotoCORM, catalyst, luminescent probe?

Mn(I) α-diimine carbonyl complexes have shown promise in the development of luminescent CO release materials (photoCORMs) for diagnostic and medical applications due to their ability to balance the energy of the low-lying metal-to-ligand charge transfer (MLCT) and metal-centered (MC) states. In this work, the excited state dynamics of [Mn(im)(CO)3(phen)]+ (im = imidazole; phen = 1,10-phenanthroline) is investigated by means of wavepacket propagation on the potential energy surfaces associated with the 11 low-lying Sn singlet excited states within a vibronic coupling model in a (quasi)-diabatic representation including 16 nuclear degrees of freedom. The results show that the early time photophysics (<400 fs) is controlled by the interaction between two MC dissociative states, namely, S5 and S11, with the lowest S1-S3 MLCT bound states. In particular, the presence of S1/S5 and S2/S11 crossings within the diabatic picture along the Mn-COaxial dissociative coordinate (qMn-COaxial) favors a two-stepwise population of the dissociative states, at about 60-70 fs (S11) and 160-180 fs (S5), which reaches about 10% within 200 fs. The one-dimensional reduced densities associated with the dissociative states along qMn-COaxial as a function of time clearly point to concurrent primary processes, namely, CO release vs entrapping into the S1 and S2 potential wells of the lowest luminescent MLCT states within 400 fs, characteristics of luminescent photoCORM.

Optical absorption properties of metal-organic frameworks: solid state versus molecular perspective.

The vast chemical space of metal and ligand combinations in Transition Metal Complexes (TMCs) gives rise to a rich variety of electronic excited states with local and non-local character such as intra-ligand (IL), metal-centered (MC), metal-to-ligand (MLCT) or ligand-to-metal charge-transfer (LMCT) states. Those features are equally found in metal organic frameworks (MOFs), defined as modular materials built from metal-nodes connected through organic-ligands. Because of the electronic and structural complexity of MOFs, the computational description of their excited states is a formidable challenge for which two different approaches have been usually followed: the solid state and the molecular perspective. The first consists in analysing the frontier electronic bands and crystal orbitals of the electronic ground state (GS) in periodic boundary conditions, while the latter points to an accurate computation of the excited states in representative clusters at the molecular level. Herein, we apply both approaches to evaluate the optical absorption properties of three experimentally reported Ti(iv) mononuclear MOFs with in silico metal substitutions with Zn(ii), Cd(ii), Fe(ii), Ru(ii) and Zr(iv) ions, thus covering d10, d6 and d0 electronic configurations of 1st and 2nd row TMCs in MOFs. Our analysis captures the main electronic features attributed to these systems while we discuss the main advantages and drawbacks of both approximations.

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.

Computational Assessment of MLCT versus MC Stabilities in First-to-Third-Row d6 Pseudo-Octahedral Transition Metal Complexes.

An accurate modeling of metal-to-ligand-charge-transfer (MLCT) and metal-centered (MC) excited state energies is key to predict the photoinduced response in transition metal complexes (TMCs). Herein, the importance of the ground state and excited state reference geometries is addressed for three-prototype d6 pseudo-octahedral TMCs, each displaying a different potential energy landscape of MLCT versus MC relative stabilities. Several functionals are used within the time-dependent density functional theory (TDDFT), as well as multireference wave-function theory (MS-CASPT2), applied to [Mn(im)(CO)3 (phen)]+ , [Ru(im)2 (bpy)2 ]2+ , and [Re(im)(CO)3 (phen)]+ , (im: imidazole, bpy: bypiridine, phen: phenantroline). The results revel that TDDFT is robust except when using B3LYP functional for first-row d6 TMCs. In contrast, MS-CASPT2 calculations are strongly biased in those cases with competitive MLCT/MC states. The results reinforce the reliability of B3LYP to describe the excited states in d6 TMCs, but question the validity of assessing the density functional theory (DFT)/TDDFT performance via direct comparison with MS-CASPT2 performed at the same DFT reference geometry as a standard strategy. © 2019 Wiley Periodicals, Inc.

Excited-State Reactivity of [Mn(im)(CO)3 (phen)]+ : A Structural Exploration.

The electronic excited state reactivity of [Mn(im)(CO)3 (phen)]+ (phen = 1,10-phenanthroline; im = imidazole) ranging between 420 and 330 nm have been analyzed by means of relativistic spin-orbit time-dependent density functional theory and wavefunction approaches (state-average-complete-active-space self-consistent-field/multistate CAS second-order perturbation theory). Minimum energy conical intersection (MECI) structures and connecting pathways were explored using the artificial force induced reaction (AFIR) method. MECIs between the first and second singlet excited states (S1 /S2 -MECIs) were searched by the single-component AFIR (SC-AFIR) algorithm combined with the gradient projection type optimizer. The structural, electronic, and excited states properties of [Mn(im)(CO)3 (phen)]+ are compared to those of the Re(I) analogue [Re(im)(CO)3 (phen)]+ . The high density of excited states and the presence of low-lying metal-centered states that characterize the Mn complex add complexity to the photophysics and open various dissociative channels for both the CO and imidazole ligands. © 2018 Wiley Periodicals, Inc.

Charge-Transfer versus Charge-Separated Triplet Excited States of [ReI (dmp)(CO)3 (His124)(Trp122)]+ in Water and in Modified Pseudomonas aeruginosa Azurin Protein.

A computational investigation of the triplet excited states of a rhenium complex electronically coupled with a tryptophan side chain and bound to an azurin protein is presented. In particular, by using high-level molecular modeling, evidence is provided for how the electronic properties of the excited-state manifolds strongly depend on coupling with the environment. Indeed, only upon explicitly taking into account the protein environment can two stable triplet states of metal-to-ligand charge transfer or charge-separated nature be recovered. In addition, it is also demonstrated how the rhenium complex plus tryptophan system in an aqueous environment experiences too much flexibility, which prevents the two chromophores from being electronically coupled. This occurrence disables the formation of a charge-separated state. The successful strategy requires a multiscale approach of combining molecular dynamics and quantum chemistry. In this context, the strategy used to parameterize the force fields for the electronic triplet states of the metal complex is also presented.

Intriguing Effects of Halogen Substitution on the Photophysical Properties of 2,9-(Bis)halo-Substituted Phenanthrolinecopper(I) Complexes.

Three new copper(I) complexes [Cu(LX)2]+(PF6-) (where LX stands for 2,9-dihalo-1,10-phenanthroline and X = Cl, Br, and I) have been synthesized in order to study the impact of halogen substituents tethered in the α position of the chelating nitrogen atoms on their physical properties. The photophysical properties of these new complexes (hereafter named Cu-X) were characterized in both their ground and excited states. Femtosecond ultrafast spectroscopy revealed that early photoinduced processes are faster for Cu-I than for Cu-Cl or Cu-Br, both showing similar behaviors. Their electronic absorption and electrochemical properties are comparable to benchmark [Cu(dmp)2]+ (where dmp stands for 2,9-dimethyl-1,10-phenanthroline); furthermore, their optical features were fully reproduced by time-dependent density functional theory and ab initio molecular dynamics calculations. All three complexes are luminescent at room temperature, showing that halogen atoms bound to positions 2 and 9 of phenanthroline are sufficiently bulky to prevent strong interactions between the excited Cu complexes and solvent molecules in the coordination sphere. Their behavior in the excited state, more specifically the extent of the photoluminescence efficiency and its dependence on the temperature, is, however, strongly dependent on the nature of the halogen. A combination of ultrafast transient absorption spectroscopy, temperature-dependent steady-state fluorescence spectroscopy, and computational chemistry allows one to gain a deeper understanding of the behavior of all three complexes in their excited state.

Absorption Spectroscopy and Photophysics of a ReI -dppz Probe for DNA-Mediated Charge Transport.

Optical properties of [Re(CO)3 (dppz)(py)]+ (dppz=dipyrido[3,2-a:2',3'-c]phenazine; py=pyridine) in acetonitrile, water and DNA have been investigated based on DFT, time-dependent-DFT (TD-DFT)/ conductor-like screening model, with and without explicit solvent molecules, and molecular dynamics. Whereas implicit solvent model is not appropriate to model optical properties of dppz-substituted metal complexes, adding explicit solvent molecules in interaction with dppz stabilizes the metal-to-ligand-charge-transfer (MLCT) transitions. Classical molecular dynamics simulations point to an important conformational flexibility, as evidenced by the coexistence of two conformers A and B. When considering the conformational sampling, the lowest band of the absorption spectrum is red-shifted and broadened up to 500 nm in agreement with the experimental spectra supporting important dynamical effects. The absorption spectra of [Re(CO)3 (dppz)(py-R)]+/ GC-DNA and [Re(CO)3 (dppz)(py-R)]+ /AT-DNA (R=CH2 -CH2 -COO- ) intercalated in both major or minor grooves exhibit a lowest energy charge separated (CS) band at about 600 nm and 500 nm, respectively, corresponding mainly to excitations from guanine and adenine to dppz. These states may play a central role into DNA-mediated charge transport processes. The over stabilization of the lowest 3 ILdppz state of [Re(CO)3 (dppz)(py)]+ in water as compared to acetonitrile could be responsible for the quenching of emission in water.

Interstate vibronic coupling constants between electronic excited states for complex molecules.

In the construction of diabatic vibronic Hamiltonians for quantum dynamics in the excited-state manifold of molecules, the coupling constants are often extracted solely from information on the excited-state energies. Here, a new protocol is applied to get access to the interstate vibronic coupling constants at the time-dependent density functional theory level through the overlap integrals between excited-state adiabatic auxiliary wavefunctions. We discuss the advantages of such method and its potential for future applications to address complex systems, in particular, those where multiple electronic states are energetically closely lying and interact. We apply the protocol to the study of prototype rhenium carbonyl complexes [Re(CO)3(N,N)(L)]n+ for which non-adiabatic quantum dynamics within the linear vibronic coupling model and including spin-orbit coupling have been reported recently.

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 π 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.

Lattice-Solvent Effects in the Spin-Crossover of an Fe(II)-Based Material. The Key Role of Intermolecular Interactions between Solvent Molecules.

The spin transition of Fe(II) complexes is the subject of intensive synthetic and computational efforts. In this manuscript, we analyze the spin crossover (SCO) of [Fe(E-dpsp)2]2+ (1), which features a spin transition depending on the cocrystallizing solvent molecules. Whereas the use of acetone results in a hysteretic spin transition at ∼170 K, the use of propylene carbonate (PC) results in a permanent diamagnetic signal up to 300 K. By means of DFT+U+D2 calculations in the solid state of the material, we unravel the reasons for such different behavior. Our results allow us to ascribe the relatively low transition temperature of 1(BF4)2·acetone to the distorted arrangement of the SCO molecules in the low-spin state of the material. In turn, intermolecular interactions play the primary role in the case of 1(BF4)2·2PC. In particular, we found that solvent-solvent interactions actively promote the stability of the low-spin state due to the formation of PC dimers. These dimers would appear at larger distances in the high-spin phase, with the subsequent loss of phase stability. This is yet another proof of how subtle is the spin transition phenomenon in Fe(II)-based architectures.

Understanding room-temperature π-dimerisation of radical ions: intramolecular π-[TTF]22+ in functionalised calix[4]arenes.

Long, multicentre π-dimers of radical ions are weakly bound and can only be observed in solution at low temperature. However, recent supramolecular approaches induce the extra stabilisation required to preserve them at room temperature, by different means depending on the approach. In particular, π-[TTF]22+ dimers (TTF = tetrathiafulvalene) were detected upon oxidation of a TTF-based calix[4]arene in acetonitrile solution at room temperature, manifesting intramolecular [R-TTF]˙+[R-TTF]˙+ interactions (Chem. Commun. 2006, 2, 2233). In this work, the reasons behind the remarkable formation of these π-dimers in the calix[4]arene, [calix], molecule are unravelled by means of DFT calculations. We first demonstrate that the properties of the π-[R-TTF]22+ dimers are preserved in the [calix]2+. Most importantly, our results show that the π-dimerised and non-dimerised forms of the [calix]2+ are isoenergetic at room temperature, and that the activation energy for this process is ca. 9.5 kcal mol-1. Hence, both forms coexist in equilibrium at 298 K, as the intramolecular nature of the interaction ensures a high reaction rate. The role of the Na+ cation in preventing the π-[R-TTF]22+ dimerisation of the [calix]2+ receptor is also examined, unveiling that this effect is mostly due to the electrostatic repulsion induced by the cation. Finally, we provide a revision on room-temperature stable supramolecular long, multicentre π-dimers of radical ions, a class of systems with great potential as molecular switches.

Excited-states of a rhenium carbonyl diimine complex: solvation models, spin-orbit coupling, and vibrational sampling effects.

We present a quantum-chemical investigation of the excited states of the complex [Re(CO)3(Im)(Phen)]+ (Im = imidazole; Phen = 1,10-phenanthroline) in solution including spin-orbit couplings and vibrational sampling. To this aim, we implemented electrostatic embedding quantum mechanics/molecular mechanics (QM/MM) in the Amsterdam Density Functional program suite, suitable for time-dependent density functional calculations including spin-orbit couplings. The new implementation is employed to simulate the absorption spectrum of the complex, which is compared to the results of implicit continuum solvation and frozen-density embedding. Molecular dynamics simulations are used to sample the ground state conformations in solution. The results demonstrate that any study of the excited states of [Re(CO)3(Im)(Phen)]+ in solution and their dynamics should include extensive sampling of vibrational motion and spin-orbit couplings.

Description of excited states in [Re(Imidazole)(CO)3 (Phen)](+) including solvent and spin-orbit coupling effects: Density functional theory versus multiconfigurational wavefunction approach.

The low-lying electronic excited states of [Re(imidazole)(CO)3 (phen)](+) (phen = 1,10-phenanthroline) ranging between 420 nm and 330 nm have been calculated by means of relativistic spin-orbit time-dependent density functional theory (TD-DFT) and wavefunction approaches (state-average-CASSCF/CASPT2). A direct comparison between the theoretical absorption spectra obtained with different methods including SOC and solvent corrections for water points to the difficulties at describing on the same footing the bands generated by metal-to-ligand charge transfer (MLCT), intraligand (IL) transition, and ligand-to-Ligand- charge transfer (LLCT). While TD-DFT and three-roots-state-average CASSCF (10,10) reproduce rather well the lowest broad MLCT band observed in the experimental spectrum between 420 nm and 330 nm, more flexible wavefunctions enlarged either by the number of roots or by the number of active orbitals and electrons destabilize the MLCT states by introducing IL and LLCT character in the lowest part of the absorption spectrum. © 2016 Wiley Periodicals, Inc.

On the zeroth-order hamiltonian for CASPT2 calculations of spin crossover compounds.

Complete active space self-consistent field theory (CASSCF) calculations and subsequent second-order perturbation theory treatment (CASPT2) are discussed in the evaluation of the spin-states energy difference (ΔH(elec)) of a series of seven spin crossover (SCO) compounds. The reference values have been extracted from a combination of experimental measurements and DFT + U calculations, as discussed in a recent article (Vela et al., Phys Chem Chem Phys 2015, 17, 16306). It is definitely proven that the critical IPEA parameter used in CASPT2 calculations of ΔH(elec), a key parameter in the design of SCO compounds, should be modified with respect to its default value of 0.25 a.u. and increased up to 0.50 a.u. The satisfactory agreement observed previously in the literature might result from an error cancellation originated in the default IPEA, which overestimates the stability of the HS state, and the erroneous atomic orbital basis set contraction of carbon atoms, which stabilizes the LS states.

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

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

Electronic and Photophysical Properties of [Re (L)(CO)3(phen)](+) and [Ru(L)2(bpy)2](2+) (L = imidazole), Building Units for Long-Range Electron Transfer in Modified Blue Copper Proteins.

The electronic, optical, and photophysical properties of [Re(im)(CO)3(phen)](+) and [Ru(bpy)2(im)2](2+) (im = imidazole; phen = 1,10-phenanthroline; bpy = 2,2'-bipyridine) in water, including spin-orbit coupling (SOC) effects, were studied by means of density functional theory (DFT) and time-dependent DFT. The main features of the visible experimental absorption spectra of both molecules are well-reproduced. Whereas the theoretical spectrum of the Re(I) complex is characterized by one metal-to-ligand charge transfer (MLCTphen) state of low intensity at 394 nm and a strongly absorbing MLCTphen state calculated at 370 nm, the spectrum of the Ru(II) complex presents a high density of singlet MLCTbpy excited states with significant oscillator strengths that contribute to the two broad bands centered at 490 and 340 nm. The absorption spectrum of [Re(im) (CO)3(phen)](+) is perturbed by SOC with non-negligible mixing between the low-lying triplet and singlet absorbing states, while SOC has no effect on the absorption spectrum of [Ru(bpy)2(im)2](2+). A detailed structural investigation of the two lowest singlet and four lowest triplet excited states of the Re(I) complex point to MLCTphen (S1, S2, T1, T2) and intra-ligand ILphen (T3) localized spin-densities characterized by small contractions from both Re-N and phen CC central bonds in the MLCT states and nearly no deformation in the IL state. A mechanism of luminescent decay of [Re(im) (CO)3(phen)](+) is proposed on the basis of the calculated energy minima and wavelengths of emission for the interpretation of the three frequency/time-scale signals put in evidence by ultrafast experiments. The long-lived emissive properties of [Ru(bpy)2(im)2](2+) are analyzed on the basis of the relative energies of the two lowest (3)MLCTbpy and metal-centered (3)MC excited states. The minimum corresponding to the (3)MC spin density shows a significant structural rearrangement with an increase of the Ru-N bond distance of 0.33 Å and a closure of the N-Ru-N bond angle of 20° inducing a large distortion of the octahedral motif. The spin-density associated with the lowest (3)MLCTbpy localized on one bpy ligand suggests the presence of a second degenerate (3)MLCTbpy minimum. The luminescence of the Ru(II) complex calculated at 669 nm is partially quenched by the presence of the low (3)MC nonradiative state at 1064 nm. When interacting with modified metal-based proteins the two complexes will behave differently because of these distinctive photophysical properties.

Unravelling the Key Driving Forces of the Spin Transition in π-Dimers of Spiro-biphenalenyl-Based Radicals.

Spiro-biphenalenyl (SBP) boron radicals constitute an important family of molecules for the preparation of functional organic materials. The building blocks of several SBP-based crystals are π-dimers of these radicals, in which two phenalenyl (PLY) rings face each other and the other two PLYs point away from the superimposed PLYs. The dimers of ethyl-SBP and butyl-SBP undergo a spin transition between a diamagnetic and a paramagnetic state upon heating, while other dimers exhibit paramagnetism at all temperatures. Here, we present a computational study aimed at establishing the driving forces of the spin-transition undergone by ethyl-SBP at ∼140 K. The ground state of the π-dimers below 140 K is a singlet state in which the SBP unpaired electrons are partially localized in the superimposed PLYs. Above 140 K, the unpaired electrons are localized in the nonsuperimposed PLYs. These high-temperature structures are exclusively governed by the ground triplet state because the open-shell singlet with the unpaired electrons localized in the nonsuperimposed PLYs does not feature any minimum in the potential energy surface of the system. Furthermore, we show that the electrostatic component of the interaction energy between SBP radicals in the π-dimers is more attractive in the triplet than in the singlet, thereby partially counteracting the bonding and dispersion components, which favor the singlet. This electrostatic stabilization of the triplet state is a key driving force of the spin transition of ethyl-SBP and a key factor explaining the paramagnetic response of the π-dimers of other SBP-based crystals.

Towards an accurate and computationally-efficient modelling of Fe(II)-based spin crossover materials.

The DFT + U methodology is regarded as one of the most-promising strategies to treat the solid state of molecular materials, as it may provide good energetic accuracy at a moderate computational cost. However, a careful parametrization of the U-term is mandatory since the results may be dramatically affected by the selected value. Herein, we benchmarked the Hubbard-like U-term for seven Fe(ii)N6-based pseudo-octahedral spin crossover (SCO) compounds, using as a reference an estimation of the electronic enthalpy difference (ΔHelec) extracted from experimental data (T1/2, ΔS and ΔH). The parametrized U-value obtained for each of those seven compounds ranges from 2.37 eV to 2.97 eV, with an average value of U = 2.65 eV. Interestingly, we have found that this average value can be taken as a good starting point since it leads to an unprecedented mean absolute error (MAE) of only 4.3 kJ mol(-1) in the evaluation of ΔHelec for the studied compounds. Moreover, by comparing our results on the solid state and the gas phase of the materials, we quantify the influence of the intermolecular interactions on the relative stability of the HS and LS states, with an average effect of ca. 5 kJ mol(-1), whose sign cannot be generalized. Overall, the findings reported in this manuscript pave the way for future studies devoted to understand the crystalline phase of SCO compounds, or the adsorption of individual molecules on organic or metallic surfaces, in which the rational incorporation of the U-term within DFT + U yields the required energetic accuracy that is dramatically missing when using bare-DFT functionals.

Keys for the existence of stable dimers of bis-tetrathiafulvalene (bis-TTF)-functionalized molecular clips presenting [TTF](•+)···[TTF](•+) long, multicenter bonds at room temperature.

A systematic theoretical and computational investigation is performed to determine the keys governing the existence, in acetonitrile solutions, of dimers of bis-tetrathiafulvalene (bis-TTF)-functionalized diphenylglycoluril molecular clips (clip2(n+)) that are stable at room temperature for n ≤ 4. Although the experimental structure of these dimers in solution is unknown, electronic absorption studies suggest that they have [TTF](l+)···[TTF](m+) interactions that are preserved at room temperature (note that when l = m = 1 these interactions become long, multicenter bonds). In good agreement with the interpretation of the experimental spectroscopic data, all clip2(n+) dimers whose charge is ≤4 present an optimum geometry that, in all cases, has three short interfragment [TTF](l+)···[TTF](m+) interactions. The computed ΔG(298 K) for these optimum structures matches the available experimental data on the stability of these dimers. Such optimum geometry, combined with the zwitterionic character of the electron distribution in monomers and dimers (most of the net positive charge is equally distributed among the TTF groups, while a 1- au charge is located in the central fused five-membered rings) allows the formation of a maximum of two long, multicenter [TTF](•+)···[TTF](•+) bonds when all TTF groups host a 1+ au of charge, as in clip2(4+). However, these long, multicenter bonds alone do not account for the stability of clip2(n+) dimers at room temperature. Instead, the studies carried out here trace the origin of their stability to (1) the zwitterionic character of their charge distribution, (2) the proper geometrical shape of the interacting monomers, which allows the intercalation of their arms, thus making possible the simultaneous formation of two short contacts, both involving the positively charged TTF group of one monomer and the negatively charged central ring of the other, (3) the simultaneous presence of three short contacts among the TTF groups in the optimum geometry of the clip2(n+) dimers, which become two long, multicenter bonds and one van der Waals interaction when the four TTF groups host a 1+ charge, and (4) the net stabilizing effect of the solvent.