Search

Supramolecular gels form a class of soft materials that has been heavily explored by the chemical community in the past 20 years. While a multitude of experimental techniques has demonstrated its usefulness when characterizing these materials, the potential value of computational techniques has received much less attention. This review aims to provide a complete overview of studies that employ computational tools to obtain a better fundamental understanding of the self-assembly behavior of supramolecular gels or to accelerate their development by means of prediction. As such, we hope to stimulate researchers to consider using computational tools when investigating these intriguing materials. In the concluding remarks, we address future challenges faced by the field and formulate our vision on how computational methods could help overcoming them. ispartof: Gels vol:7 issue:3 pages:1-18 ispartof: location:Switzerland status: Published online

Non-covalent interactions involving aromatic rings play a central role in many areas of modern chemistry. In medicinal and bioorganic chemistry, the intermolecular interactions between the aromatic side chains of amino acids, such as phenylalanine and tyrosine, are of great interest. To enhance the affinity between such aromatic side chains, halogenation is a promising modification strategy. In the current work, the nature and strength of halogenated p-p stacked phenylalanine (Phe) dimers have been investigated using density functional theory, energy decomposition analyses and the non-covalent interaction (NCI) method. Our analysis shows that increasing the degree of halogenation enhances the strength of the stacking interactions and, moreover, the heavier halides (Cl, Br and I) lead to stronger interactions compared to the lighter F. This effect was traced back to local secondary interactions of the halide with the aliphatic C-H bonds of the phenylalanine side chain. Based on the computational findings, a set of peptide hydrogelators was synthesized, and the resulting hydrogel properties were further investigated via dynamic rheometry. Experimental observations can be correlated to the trends found in the theoretical analysis, suggesting that local interactions indeed play a noticeable role in enhancing peptide-based hydrogel strength. This journal is

Metabolomics is recognized as a crucial scientific domain, promising to advance our understanding of cell biology, physiology and medicine. Tandem mass spectrometry (MS/MS) has strong potential to elucidate the metabolites in complex biological samples and to become a standard tool complementing established techniques. However, despite its potential for answering many key questions, a major challenge in the use of tandem mass spectrometry for characterizing metabolites lies in a lack of computational tools for accurate annotation and structure identification allowing us to turn complex data into molecular knowledge. Chemo-informatics and related machine-learning in silico fragmentation tools have already been established and used for different classes of metabolites. For the classes of metabolites where existing chemo-informatics approaches produce insufficiently accurate predictions a supervised machine learning based strategy can be used to predict possible molecular structures from “unassigned” experimental tandem MS data. Here, we propose a new innovative in silico approach employing quantum mechanical (QM) methods in order to predict ion formation and subsequent fragmentation patterns of arbitrary small molecules and validate putative annotations of tandem mass spectrometry (MS) data. The focus is on the evaluation of a new conceptual density functional theory (CDFT) nuclear reactivity descriptor of the nuclear Fukui function type, that characterizes the forces that the atomic nuclei experience due to proton attachment and captures the onset of the change in the nuclear positions induced by it. A series of test compounds for which high quality experimental data exist and that were investigated before in a more approximate theoretical framework have been examined. The output of these calculations provides a list of the most probable molecular structures predicted to match the experimental tandem MS spectrum (“de novo metabolite identification”).

Judiciously applied oriented external electric fields (OEEFs) exert catalytic effects on the kinetics and improve the thermodynamics of chemical reactions. Herein, we examine the ability of OEEFs to assist catalysts and show that the rate of oxidative addition between palladium catalysts and alkyl/aryl electrophiles can be controlled by an OEEF applied along the direction of electron reorganization (the "reaction axis"). The concerted mechanism of oxidative addition proceeds through a transition state with moderate charge transfer character. We demonstrate that OEEFs along the reaction axis can control this charge transfer and impart electrostatic catalysis. When the applied field exceeds a certain critical value (∼0.15 V/Å), we observed a mechanistic crossover from the concerted to a dissociative CSNAr type of reactivity for aryl electrophiles. To our surprise, alkyl electrophiles follow a hitherto unexplored SN2 pathway for the reaction with large transition state stabilization at relatively low OEEFs. A valence-bond state correlation diagram (VBSCD) is employed to comprehend the results. Finally, although the catalytic effect of salt additives in oxidative addition is known, its mechanism is still under debate. Our findings further show evidence that salt additives exert electric-field effects on the rate of cross-coupling reactions, and their cocatalytic effects can be judiciously reproduced by applied external electric fields. As such, we propose that the use of additives (anionic or cationic) is an experimentally viable strategy to implement external electric-field effects in routinely used oxidative addition catalysis.

With the idea of using alchemical derivatives to explore in an efficient, computer- and cost-effective way Chemical Space was launched several years ago. In the context of Conceptual DFT response functions, these energies vs nuclear charge derivatives permit the estimatation of the energy of transmutants of a given starting or reference molecule showing different nuclear compositions. After an explorative study on small and planar molecules (Balawender et al. J. Chem. Theory Comput. 2013, 9, 5327) by the present authors of this paper, the present study fully exploits the computational advantages of the alchemical derivatives in larger three-dimensional systems. Starting from a single reference calculation on C60, the complete BN substitution pattern, from single substituted C58BN via the belt (C20(BN)20 and the ball C12(BN)24 structures to the fully substituted (BN)30, is explored. Successive and simultaneous substitution strategies are followed and compared, indicating that both techniques yield identical results up to 13 substitutions but that for higher substitutions the simultaneous approach needs to be taken. Due to the cost-efficiency of the algorithm this path can indeed be followed as opposed to earlier work in the literature where for each step a full SCF calculation was at stake leading to prohibitively large computational demands for adopting the simultaneous approach. Previously formulated rules governing the substitution pattern by Kar and co-workers are scrutinized in this context and reformulated giving chemical insight in the gradual substitution process and the relative energies of the isomers. In its present form the method offers an interesting venue to study BN substitution patterns in higher fullerenes and graphene and in general paves the way for more efficient exploration of the Chemical Space.

The most relevant manifestations of ligand noninnocence of quinone and bipyridine derivatives are thoroughly scrutinized and discussed through an extensive and systematic set of octahedral ruthenium complexes, [(en)2RuL]z, in four oxidation states (z = +3, +2, +1, and 0). The characteristic structural deformation of ligands upon coordination/noninnocence is put into context with the underlying electronic structure of the complexes and its change upon reduction. In addition, by means of decomposing the corresponding reductions into electron transfer and structural relaxation subprocesses, the energetic contribution of these structural deformations to the redox energetics is revealed. The change of molecular electron density upon metal- and ligand-centered reductions is also visualized and shown to provide novel insights into the corresponding redox processes. Moreover, the charge distribution of the π-subspace is straightforwardly examined and used as indicator of ligand noninnocence in the distinct oxidation states of the complexes. The aromatization/dearomatization processes of ligand backbones are also monitored using magnetic (NICS) and electronic (PDI) indicators of aromaticity, and the consequences to noninnocent behavior are discussed. Finally, the recently developed effective oxidation state (EOS) analysis is utilized, on the one hand, to test its applicability for complexes containing noninnocent ligands, and, on the other hand, to provide new insights into the magnitude of state mixings in the investigated complexes. The effect of ligand substitution, nature of donor atom, ligand frame modification on these manifestations, and measures is discussed in an intuitive and pedagogical manner M.G. and P.S. are thankful for financial support from the FEDER grant UNGI10-4E-801 and MINECO grant CTQ2014-59212-P/BQU

Pauling long bond orders have previously been used to obtain qualitative insight in the transmission through nanographenes. Here we show that this long bond order, the atom–atom polarizability (a measure for delocalization) and the transmission probability are intimately linked and that their relationships are valid for all alternant hydrocarbons. These relationships allow a simple rationalization of the transport properties of a variety of molecules considered in molecular electronics. As an example, some molecular wires such as oligo( p -phenylene) are studied, leading to a simple explanation for the experimentally observed exponential decay of the transmission probability with the number of phenyl units.

The nature and origin of ion-π and ion-σ interactions has been systematically investigated using dispersion-corrected density functional theory and the recently developed noncovalent interaction (NCI) method. A detailed analysis of these interactions is performed with the aim to identify the requirements that have to be fulfilled by the molecular system for strong ion-ligand interactions. Interestingly, our results indicate that aliphatic systems, such as cyclohexane, can interact as strong as aromatic ones with both cations and anions, despite of having a negligible quadrupole moment. In fact, cyclohexane binds anions stronger than benzene itself but slightly weaker that hexafluorobenzene. The NCI method reveals that the interaction between the ions and three C–H bonds of the saturated fragment are responsible for the surprisingly strong ion-σ interaction. A weakening of the ion-σ interactions is observed in the order: Li + > F − > Na + > Cl − > Br − ≈ K + . In addition, a complete Ziegler–Rauk type energy decomposition analysis has been carried out in order to reveal the origins of the thermodynamic driving force for complex formations. The electron density deformation upon complex formation has been scrutinized with a complementary NOCV analysis allowing the identification of molecular orbital interaction contributions to the stabilization. Based on these analysis, it is shown that the formally anion-π interaction is rather an anion-σ ∗ interaction.

Metal and ligand-based reductions have been modeled in octahedral ruthenium complexes revealing metal–ligand interactions as the profound driving force for the redox-active behaviour of orthoquinoid-type ligands. Through an extensive investigation of redox-active ligands we revealed the most critical factors that facilitate or suppress redox-activity of ligands in metal complexes, from which basic rules for designing non-innocent/redox-active ligands can be put forward. These rules also allow rational redox-leveling, i.e. the moderation of redox potentials of ligand-centred electron transfer processes, potentially leading to catalysts with low overpotential in multielectron activation processes.

The bonding in heavier Group 14 zero-valent complexes of a general formula L2 E (E=Si-Pb; L=phosphine, N-heterocyclic and acyclic carbene, cyclic tetrylene and carbon monoxide) is probed by combining valence bond (VB) theory and maximum probability domain (MPD) approaches. All studied complexes are initially evaluated on the basis of the structural parameters and the shape of frontier orbitals revealing a bent structural motif and the presence of two lone pairs at the central E atom. For the VB calculations three resonance structures are suggested, representing the "ylidone", "ylidene" and "bent allene" structures, respectively. The influence of both ligands and central atoms on the bonding situation is clearly expressed in different weights of the resonance structures for the particular complexes. In general, the bonding in the studied E0 compounds, the tetrylones, is best described as a resonating combination of "ylidone" and "ylidene" structures with a minor contribution of the "bent allene" structure. Moreover, the VB calculations allow for a straightforward assessment of the π-backbonding (E→L) stabilization energy. The validity of the suggested resonance model is further confirmed by the complementary MPD calculations focusing on the E lone pair region as well as the E-L bonding region. Likewise, the MPD method reveals a strong influence of the σ-donating and π-accepting properties of the ligand. In particular, either one single domain or two symmetrical domains are found in the lone pair region of the central atom, supporting the predominance of either the "ylidene" or "ylidone" structures having one or two lone pairs at the central atom, respectively. Furthermore, the calculated average populations in the lone pair MPDs correlate very well with the natural bond orbital (NBO) populations, and can be related to the average number of electrons that is backdonated to the ligands.

In this study, an intuitive concept is derived, which explains the characteristic dependence of spin-state energetics on the exact exchange admixture of DFT functionals in the case of octahedral transition metal complexes. The change in electron density distributions upon varying the admixture, c 3, in the B3LYP functional is analyzed for archetype ionic and covalent systems as well as for the Fe 2+ ion in an ideal octahedral field. An understanding of how the DFT description of the electronic structure of octahedral complexes changes as a function of c 3 is sought. A systematic spin-state energy analysis of 50 octahedral complexes of various metals and ligands with consistent experimental data is presented, allowing the derivation, in theory, of an optimal c 3 value for each system. The notion that the admixture dependence of DFT spin-state energetics stems from the treatment of nondynamic electrons arising from the mixing of (M–L z2) 0(d z2) 2 and (M–L x2-y2 0(d( x2-y2)) 2 configurations into the dominant (M–L x2-y2) 2(d x2-y2) 0 and (M–L x2-y2) 2(d x2-y2) 0 ones in the low(er) spin states is put forward. That is, in the effort to mimic such electron–electron interactions, E x LDA overestimates, whereas exact exchange downplays the contribution of this type of electron correlation to the stability of low(er) spin states, leading to the widespread practical observation that the higher the exact exchange admixture, the more stable the high-spin-state configuration.

The chalcogen bond has been acknowledged as an influential noncovalent interaction (NCI) between an electron-deficient chalcogen (donor) and a Lewis base (acceptor). This work explores the main features of chalcogen bonding through a large-scale computational study on a series of donors and acceptors spanning a wide range in strength and character of this type of bond: (benzo)chalcogenadiazoles (with Ch = Te/Se/S) versus halides and neutral Lewis bases with O, N, and C as donor atoms. We start from Pearson's hard and soft acids and bases (HSAB) principle, where the hard nature of the chalcogen bond is quantified through the molecular electrostatic potential and the soft nature through the Fukui function. The σ-holes are more pronounced when going down in the periodic table and their directionality matches the structural orientation of donors and acceptors in the complexes. The Fukui functions point toward an n→σ*-type interaction. The initial conjectures are further scrutinized using quantum mechanical methods, mostly relating to the systems' electron density. A Ziegler–Rauk energy decomposition analysis shows that electrostatics plays a distinctly larger role for the soft halides than for the hard, uncharged acceptors, associated with the softness matching within the HSAB principle. The natural orbital for chemical valence analysis confirms the n→σ* electron donation mechanism. Finally, the electron density and local density energy at the bond critical point in the quantum theory of atoms in molecules study and the position of the spikes in the reduced density gradient versus density plot in the NCI theory situate the chalcogen bond in the same range as strong hydrogen bonds. © 2017 Wiley Periodicals, Inc.

Even though the concept of aromaticity and antiaromaticity is extremely important and widely used, there still exist lots of controversies in the literature, which are believed to be originated from the fact that there are so many aromatic types discovered and at the same time there are many aromaticity indexes proposed. In this work, using seven series of substituted fulvene derivatives as an example and with the information-theoretic approach in density functional reactivity theory, we examine these concepts from a different perspective. We investigate the changing patterns of Shannon entropy, Fisher information, Ghosh–Berkowitz–Parr entropy, information gain, Onicescu information energy, and relative Renyi entropy on the ring carbon atoms of these systems. Meanwhile, we also consider variation trends of four representative kinds of aromaticity indexes such as FLU, HOMA, ASE and NICS. Statistical analyses among these quantities show that with the same ring structure of the derivatives, both information-theoretic quantities and aromaticity indexes obey the same changing pattern, which are valid across all seven systems studied. However, cross correlations between these two sets of quantities yield two completely opposite patterns. These ring-structure dependent correlations are in good agreement with Hückel's 4n + 2 rule of aromaticity and 4n rule of antiaromaticity. Our results should provide a novel and complementary viewpoint on how aromaticity and antiaromaticity should be appreciated and categorized. More studies are in progress to further our understanding about the matter.

Functionalized diamondoids show great potential as building blocks for various new optoelectronic applications. However, until now, only simple mono and double substitutions were investigated. In this work, we considered up to 10 and 6 sites for functionalization of the two smallest diamondoids, adamantane and diamantane, respectively, in search for diamondoid derivatives with a minimal and maximal HOMO-LUMO energy gap. To this end, the energy gap was optimized systematically using an inverse molecular design methodology based on the best-first search algorithm combined with a Monte Carlo component to escape local optima. Adamantane derivatives were found with HOMO-LUMO gaps ranging from 2.42 to 10.63 eV, with 9.45 eV being the energy gap of pure adamantane. For diamantane, similar values were obtained. The structures with the lowest HOMO-LUMO gaps showed apparent push-pull character. The push character is mainly formed by sulfur or nitrogen dopants and thiol groups, whereas the pull character is predominantly determined by the presence of electron-withdrawing nitro or carbonyl groups assisted by amino and hydroxyl groups via the formation of intramolecular hydrogen bonds. In contrast, maximal HOMO-LUMO gaps were obtained by introducing numerous electronegative groups.

We consider the problem of defining an appropriate local descriptor corresponding to an arbitrary global descriptor. Although it does not seem easy to rigorously and uniquely define local analogues of derived global descriptors (e.g., the electrophilicity) or the fundamental global descriptors associated with the canonical ensemble (e.g., the hardness), the local response of these global descriptors can be defined unambiguously. We look at the local response of the global electrophilicity and compare it to the conventional, ad hoc, definition of the local electrophilicity. The local response of global nucleofugality and electrofugality is also discussed.

To rationalize the driving factors in the competition of halogen bonding and hydrogen bonding, the complexes of the combined halogen-/hydrogen-bond donor difluoroiodomethane with the Lewis bases trimethylphosphine, dimethyl sulfide and chloromethane are studied. For all Lewis bases,ab initiocalculations lead to halogen- and hydrogen-bonded complexes. Fourier transform–IR experiments involving solutions of mixtures of difluoroiodomethane with trimethylphosphine(-d9) or dimethyl sulfide(-d6) in liquid krypton confirm the coexistence of a halogen-bonded and hydrogen-bonded complex. Also for solutions containing chloromethane, evidence of the formation of binary associations is found, but no definitive assignment of the multiple complex bands could be made. Using van't Hoff plots, the experimental complexation enthalpies for the halogen- and hydrogen-bonded complex of difluoroiodomethane with trimethylphosphine are determined to be −15.4 (4) and −10.5 (3) kJ mol−1, respectively, while for the halogen- and hydrogen-bonded complexes with dimethyl sulfide, the values are −11.3 (5) and −7.7 (6) kJ mol−1, respectively. The experimental observation that for both trimethylphospine and dimethyl sulfide the halogen-bonded complex is more stable than the hydrogen-bonded complex supports the finding that softer Lewis bases tend to favor iodine halogen bonding over hydrogen bonding.

The role of the solvent and the influence of dynamics on the kinetics and mechanism of the SNAr reaction of several halonitrobenzenes in liquid ammonia, using both static calculations and dynamic ab initio molecular dynamics simulations, are investigated. A combination of metadynamics and committor analysis methods reveals how this reaction can change from a concerted, one-step mechanism in gas phase to a stepwise pathway, involving a metastable Meisenheimer complex, in liquid ammonia. This clearly establishes, among others, the important role of the solvent and highlights the fact that accurately treating solvation is of crucial importance to correctly unravel the reaction mechanism. It is indeed shown that H-bond formation of the reacting NH3 with the solvent drastically reduces the barrier of NH3 addition. The halide elimination step, however, is greatly facilitated by proton transfer from the reacting NH3 to the solvent. Furthermore, the free energy surface strongly depends on the halide substituent and the number of electron-withdrawing nitro substituents.

In this theoretical and experimental study, the ability of carbonyl fluoride (COF2) and carbonyl chloride fluoride (COFCl) to form noncovalent interactions with the Lewis base dimethyl ether (DME) is assessed. From ab initio calculations, two stable complexes are found for COF2·DME, both formed through a lone pair···π interaction. FTIR measurements on liquefied noble gas solutions, supported by ab initio calculations, statistical thermodynamical calculations and Monte Carle Free Energy Perturbation calculations, show that a 1:1 lone pair···π bonded complex is found in solution, with an experimental complexation enthalpy of -14.5(3) kJ mol-1. For COFCl·DME three lone pair···π complexes, as well as a Cl···O halogen bonded complex, are found from ab initio calculations. Experimentally, clear complex bands for 1:1 lone pair···π complexes are observed, with an experimental complexation enthalpy of -11.4(2) kJ mol-1. Furthermore, indications of the presence of a small amount of the halogen bonded complex are also observed.

International audience; In recent years, expanded porphyrins have emerged as a promising class of π-conjugated molecules that display unique electronic, optical and conformational properties. Several expanded porphyrins can switch between planar and twisted conformations, which have different photophysical properties. Such a change of topology involves a Hückel–Möbius aromaticity switch in a single molecule and it can be induced by solvent, pH and metallation. These features make expanded porphyrins suitable for the development of a novel type of molecular switches for molecular electronic devices. Octaphyrins consisting of eight pyrrole rings, exhibit twisted-Hückel, Möbius and Hückel π-conjugation topologies depending on the oxidation and protonation state, with distinct electronic structures and aromaticity. Our working hypothesis is that a significant change in the conductance of expanded porphyrins will be observed after the topology switching. Despite the potential of Hückel–Möbius systems as conductance switches, the relationship between the conductance and the molecular topology is not yet understood. We have explored the performance of local descriptors of conductivity in simple molecules, as well as the relationship with conductance. Since these indexes provide a qualitative measure of delocalization and conductance in the probe molecules, we have carried out a local analysis of electrical conductance changes as a function of the π-conjugation in two examples. In one of them, the locality of the electronic changes ensures the ability of these indexes to describe the conductance as local. Moreover, it enables to identify which conformational switch would be more efficient from an electronic device perspective. However, we also show that it is not always possible to reduce conductance changes to one bond, and in those molecules where a deep rearrangement occurs far from the structural perturbation, local measures show a limited efficiency. This is a first step for the description of the connection between the molecular structure and conductance in molecular switches.

In this study, the Young’s and flexural moduli of single- and double-layer graphene have been theoretically investigated using periodic boundary condition (PBC) density functional theory (DFT) with the PBE, HSE06H, and M06L functionals in conjunction with the 6-31G* and the 3-21G basis sets. The unit-cell size and shape dependence as well as the directional dependencies of the mechanical properties have also been investigated. Additionally, the calculated stretching and flexural strain–stress curves are reported. Finally, finite-element simulations have been performed so as to find a homogeneous equivalent isotropic transverse material for single-layer graphene, in order to reproduce mechanical behavior in both tensile and bending sollicitations.

The fundamental principles of pericyclic reactions are governed by the Woodward-Hoffmann rules, which state that these reactions can only take place if the symmetries of the reactants' molecular orbitals and the products' molecular orbitals are the same. As such, these rules rely on the nodal structure of either the wave function or the frontier molecular orbitals, so it is unclear how these rules can be recovered in the density functional reactivity theory (or "conceptual DFT"), where the basic quantity is the strictly positive electron density. A third, nonsymmetry based approach to predict the outcome of pericyclic reactions is due to Zimmerman which uses the concept of the aromatic transition states: allowed reactions possess aromatic transition states, while forbidden reactions possess antiaromatic transition states. Based on our recent work on cycloadditions, we investigate the initial response of the chemical hardness, a central DFT based reactivity index, along the reaction profiles of a series of electrocyclizations. For a number of cases, we also compute complete initial reaction coordinate (IRC) paths and hardness profiles. We find that the hardness response is always higher for the allowed modes than for the forbidden modes. This suggests that the initial hardness response along the IRC is the key for casting the Woodward-Hoffmann rules into conceptual DFT.

The difference in reactivity between the activated 2-bromomethyl-1-tosylaziridine and the nonactivated 1-benzyl-2-(bromomethyl)aziridine with respect to sodium methoxide was analyzed by means of DFT calculations within the supermolecule approach, taking into account explicit solvent molecules. In addition, the reactivity of epibromohydrin with regard to sodium methoxide was assessed as well. The barriers for direct displacement of bromide by methoxide in methanol are comparable for all three heterocyclic species under study. However, ring opening was found to be only feasible for the epoxide and the activated aziridine, and not for the nonactivated aziridine. According to these computational analyses, the synthesis of chiral 2-substituted 1-tosylaziridines can take place with inversion (through ring opening/ring closure) or retention (through direct bromide displacement) of configuration upon treatment of the corresponding 2-(bromomethyl)aziridines with 1 equiv of a nucleophile, whereas chiral 2-substituted 1-benzylaziridines are selectively obtained with retention of configuration (via direct bromide displacement). Furthermore, the computational results showed that explicit accounting for solvent molecules is required to describe the free energy profile correctly. To verify the computational findings experimentally, chiral 1-benzyl-2-(bromomethyl)aziridines and 2-bromomethyl-1-tosylaziridines were treated with sodium methoxide in methanol. The presented work concerning the reactivity of 2-bromomethyl-1-tosylaziridine stands in contrast to the behavior of the corresponding 1-tosyl-2-(tosyloxymethyl)aziridine with respect to nucleophiles, which undergoes a clean ring-opening/ring-closure process with inversion of configuration at the asymmetric aziridine carbon atom.

Density function theory calculations were carried out on 14T clusters wherein one of the 14 Si atoms was substituted by either B, Al, or Ga. The clusters were neutralized by one alkaline cation. Two accessible oxygen atoms, O1 and O4, bonded to the substituted atom are considered for their basic properties. The computed proton affinities and NO+ binding energies show the same trends as experimentally observed. These data point out that the O4 oxygen is the most basic for the B substitution and O1 is the most basic for the Al and Ga substitutions. In view of the hard character of the oxygen atom, hard descriptors are computed. Among them, the charge neither reproduces correctly the effect of the alkaline cation nor the effect of the isomorphous substitution. The local hardness reproduces well the isomorphous substitution but fails to predict the most basic site or the effect of the alkaline cation. The molecular electrostatic potential appears to be the best hard descriptor. It predicts well the effect of ...

Radicals can be regarded as electrophilic/nucleophilic, depending on their tendency to attack sites of relatively higher/lower electron density. In this paper, an electrophilicity scale, global as well as local, and a nucleophilicity scale for 35 radicals is reported. The global electrophilicity scale correlates well with the nucleophilicity scale, suggesting that these concepts are inversely related.

In this paper, we verify the usefulness of the alchemical derivatives in the prediction of chemical properties. We concentrate on the stability of the transmutation products, where the term "transmutation" means the change of the nuclear charge at an atomic site at constant number of electrons. As illustrative transmutations showing the potential of the method in exploring chemical space, we present some examples of increasing complexity starting with the deprotonation, continuing with the transmutation of the nitrogen molecule, and ending with the substitution of isoelectronic B-N units for C-C units and N units for C-H units in carbocyclic systems. The basis set influence on the qualitative and quantitative accuracies of the alchemical predictions was investigated. The alchemical deprotonation energy (from the second order Taylor expansion) correlates well with the vertical deprotonation energy and can be used as a preliminary indicator for the experimental deprotonation energy. The results of calculations for the BN derivatives of benzene and pyrene show that this method has great potential for efficient and accurate scanning of chemical space.

The heretofore unpredictable behavior of [26] and [28]hexaphyrins upon metalation has been elucidated through quantum chemical calculations. It is demonstrated that the molecular topology of Group 10 and Group 11 metal complexes of hexaphyrins depends on sensitive interplay between the intrinsic ligand strain and the metal–ligand interaction strength. As such, the aromaticity of the ligand and effective charge of the metal are revealed as key factors determining the binding mode and the preference for Möbius or Hückel structures. These findings offer a new perspective to rationalize experimental observations for metalated hexaphyrins. More importantly, the proposed guidelines could be useful for designing novel complexes of hexaphyrins, such as a hitherto unknown Möbius [26]hexaphyrin complex.

International audience; The electron density changes from reactants towards the transition state of a chemical reaction is expressed as a linear combination of the state-specific dual descriptors (SSDD) of the corresponding reactant complexes. Consequently, the SSDD can be expected to bear important resemblance to the so-called natural orbitals for chemical valence (NOCV), introduced as the orbitals that diagonalize the deformation density matrix of interacting molecules. This agreement is shown for three case studies: the complexation of a Lewis acid with a Lewis base, a S(N)2 nucleophilic substitution reaction and a Diels-Alder cycloaddition reaction. As such, the SSDD computed for reactant complexes are shown to provide important information about charge transfer interactions during a chemical reaction.

I. Introduction: Conceptual vs Fundamental andComputational Aspects of DFT1793II. Fundamental and Computational Aspects of DFT 1795A. The Basics of DFT: The Hohenberg−KohnTheorems1795B. DFT as a Tool for Calculating Atomic andMolecular Properties: The Kohn−ShamEquations1796C. Electronic Chemical Potential andElectronegativity: Bridging Computational andConceptual DFT1797III. DFT-Based Concepts and Principles 1798A. General Scheme: Nalewajski’s ChargeSensitivity Analysis1798B. Concepts and Their Calculation 18001. Electronegativity and the ElectronicChemical Potential18002. Global Hardness and Softness 18023. The Electronic Fukui Function, LocalSoftness, and Softness Kernel18074. Local Hardness and Hardness Kernel 18135. The Molecular Shape FunctionsSimilarity 18146. The Nuclear Fukui Function and ItsDerivatives18167. Spin-Polarized Generalizations 18198. Solvent Effects 18209. Time Evolution of Reactivity Indices 1821C. Principles 18221. Sanderson’s Electronegativity EqualizationPrinciple18222. Pearson’s Hard and Soft Acids andBases Principle18253. The Maximum Hardness Principle 1829IV. Applications 1833A. Atoms and Functional Groups 1833B. Molecular Properties 18381. Dipole Moment, Hardness, Softness, andRelated Properties18382. Conformation 18403. Aromaticity 1840C. Reactivity 18421. Introduction 18422. Comparison of Intramolecular ReactivitySequences18443. Comparison of Intermolecular ReactivitySequences18494. Excited States 1857D. Clusters and Catalysis 1858V. Conclusions 1860VI. Glossary of Most Important Symbols andAcronyms1860VII. Acknowledgments 1861VIII. Note Added in Proof 1862IX. References 1865

The acid zeolite catalyzed ethylation and isopropylation of benzene with ethene and propene is investigated at the B3LYP/6-31G* level of calculation with a T 4 -cluster representing the Bronsted acid site of the zeolite. After geometry optimization of reactants, transition structures, and products, a kinetic and a sensitivity analysis were used to model the reactions. These reactions proceed in one step, without formation of alkoxy species or charged intermediates. Both reaction rate constants and local HSAB properties indicate that isopropylation is favored over ethylation, in agreement with experimental observations.

Density Functional Theory is situated within the evolution of Quantum Chemistry as a facilitator of computations and a provider of new, chemical insights. The importance of the latter branch of DFT, conceptual DFT is highlighted following Parr's dictum "to calculate a molecule is not to understand it". An overview is given of the most important reactivity descriptors and the principles they are couched in. Examples are given on the evolution of the structure-property-wave function triangle which can be considered as the central paradigm of molecular quantum chemistry to (for many purposes) a structure-property-density triangle. Both kinetic as well as thermodynamic aspects can be included when further linking reactivity to the property vertex. In the field of organic chemistry, the ab initio calculation of functional group properties and their use in studies on acidity and basicity is discussed together with the use of DFT descriptors to study the kinetics of SN2 reactions and the regioselectivity in Diels Alder reactions. Similarity in reactivity is illustrated via a study on peptide isosteres. In the field of inorganic chemistry non empirical studies of adsorption of small molecules in zeolite cages are discussed providing Henry constants and separation constants, the latter in remarkable good agreement with experiments. Possible refinements in a conceptual DFT context are presented. Finally an example from biochemistry is discussed : the influence of point mutations on the catalytic activity of subtilisin.

Protein thiol/sulfenic acid oxidation potentials provide a tool to select specific oxidation agents, but are experimentally difficult to obtain. Here, insights into the thiol sulfenylation thermodynamics are obtained from model calculations on small systems and from a quantum mechanics/molecular mechanics (QM/MM) analysis on human 2-Cys peroxiredoxin thioredoxin peroxidase B (Tpx-B). To study thiol sulfenylation in Tpx-B, our recently developed computational method to determine reduction potentials relatively compared to a reference system and based on reaction energies reduction potential from electronic energies is updated. Tpx-B forms a sulfenic acid (R-SO−) on one of its active site cysteines during reactive oxygen scavenging. The observed effect of the conserved active site residues is consistent with the observed hydrogen bond interactions in the QM/MM optimized Tpx-B structures and with free energy calculations on small model systems. The ligand effect could be linked to the complexation energies of ligand L with CH3S− and CH3SO−. Compared to QM only calculations on Tpx-B’s active site, the QM/MM calculations give an improved understanding of sulfenylation thermodynamics by showing that other residues from the protein environment other than the active site residues can play an important role.

The reactivity and complexation properties of dicoordinated Sn(II) and Sn(0) compounds are reviewed. The (dominant) electrophilicity of the stannylenes was confirmed and quantified through density functional theory (DFT) based reactivity indices. For these compounds, combining theoretical DFT calculations and experimental nuclear magnetic resonance (NMR) spectroscopic results has evidenced their potential to undergo π-complexation from aromatic π clouds in addition to significantly stronger σ-complexation. Moreover, their potential as Lewis bases was scrutinized in their interactions and reactions with iron and tungsten carbonyl Lewis acids. Finally, a prospective comparison of the reactivity of divalent stannylenes and stannylones, with a 0 oxidation state at the Sn atom, is presented.

Noncovalent interactions involving aromatic rings, such as π-stacking and CH/π interactions, are central to many areas of modern chemistry. However, recent studies proved that aromaticity is not required for stacking interactions, since similar interaction energies were computed for several aromatic and aliphatic dimers. Herein, the nature and origin of π/π, σ/σ, and σ/π dispersion interactions has been investigated by using dispersion-corrected density functional theory, energy decomposition analysis, and the recently developed noncovalent interaction (NCI) method. Our analysis shows that π/π and σ/σ stacking interactions are equally important for the benzene and cyclohexane dimers, explaining why both compounds have similar boiling points. Also, similar dispersion forces are found in the benzene⋅⋅⋅methane and cyclohexane⋅⋅⋅methane complexes. However, for systems larger than naphthalene, there are enhanced stacking interactions in the aromatic dimers adopting a parallel-displaced configuration compared to the analogous saturated systems. Although dispersion plays a decisive role in stabilizing all the complexes, the origin of the π/π, σ/σ, and σ/π interactions is different. The NCI method reveals that the dispersion interactions between the hydrogen atoms are responsible for the surprisingly strong aliphatic interactions. Moreover, whereas σ/σ and σ/π interactions are local, the π/π stacking are inherently delocalized, which give rise to a non-additive effect. These new types of dispersion interactions between saturated groups can be exploited in the rational design of novel carbon materials.

The mechanism of the decomposition of cyclopropenones and cyclopropenethiones giving acetylenes (R1–CC–R2 with R1 = R2 = H, F) plus carbon monoxide and carbon monothioxide, respectively (CX, X = O, S) has been probed using DFT (B3LYP/6-311G(d,p)) and MO (CCSD(T)/6-31G(d)) calculations. It turns out that the decomposition is stepwise, forming first an intermediate which has the properties of a semi-carbene, semi-zwitterion, whose structure is best described in both cases as a resonance hybrid between a carbene and a zwitterion. Using the PCM model, study of the solvent effect on the reactions in three continua (water, acetonitrile and benzene) shows that solvents do not affect the shape of the potential energy surfaces but tend to stabilize all the isomers. Estimation of the first vertical excited states by CIS and TD-B3LYP methods suggests that the photochemical reactions are likely to take place in the ground state rather than in an excited state. Hardness and polarizability profiles along the reaction paths show that there is a maximum in the polarizability profile besides an inverse relationship between hardness and polarizability. Fluorine substitution tends to stabilize the intermediate making it a genuinely detectable moiety.

Halogen bonds between the trifluoromethyl halides CF3Cl, CF3Br and CF3I, and dimethyl ether, dimethyl sulfide, trimethylamine and trimethyl phosphine were investigated using Pearson's hard and soft acids and bases (HSAB) concept with conceptual DFT reactivity indices, the Ziegler-Rauk-type energy-decomposition analysis, the natural orbital for chemical valence (NOCV) framework and the non-covalent interaction (NCI) index. It is found that the relative importance of electrostatic and orbital (charge transfer) interactions varies as a function of both the donor and acceptor molecules. Hard and soft interactions were distinguished and characterised by atomic charges, electrophilicity and local softness indices. Dual-descriptor plots indicate an orbital sigma hole on the halogen similar to the electrostatic s hole manifested in the molecular electrostatic potential. The predicted high halogen-bond-acceptor affinity of N-heterocyclic carbenes was evidenced in the highest complexation energy for the hitherto unknown CF3I center dot NHC complex. The dominant NOCV orbital represents an electron-density deformation according to a n ->sigma*-type interaction. The characteristic signal found in the reduced density gradient versus electron-density diagram corresponds to the non-covalent interaction between contact atoms in the NCI plots, which is the manifestation of halogen bonding within the NCI theory. The unexpected C-X bond strengthening observed in several cases was rationalised within the molecular orbital framework.

Tuning the band gap of graphene nanoribbons by chemical edge functionalisation is a promising approach towards future electronic devices based on graphene. The band gap is closely related to the aromaticity distribution and therefore tailoring the aromaticity patterns is a rational way for controlling the band gap. In the present work, it is shown how the three distinct classes of aromaticity patterns already found for armchair graphene nanoribbons can be rationally tuned by chemical edge functionalisation to modify their electronic arrangement and band gap. The electronic structure and the aromaticity distribution are studied using DFT calculations and through a series of delocalisation and geometry analysis methods, like the six-centre index (SCI) and the mean bond length (MBL) geometry descriptor. Novel aromaticity patterns are found for fluorine and nitrogen functionalisation characterised as inverted incomplete-Clar, and broken-Kekulé classes, while oxygen and nitrogen functionalisation is found to cut and extend the aromatic system, respectively. All these different arrangements of aromatic rings along the structure of graphene nanoribbons are explained using Clar's sextet theory, and a mesomeric effect mechanism for fluorine and nitrogen. In all cases, the changes in the aromaticity patterns are related to changes in the band gap. The energy and stability of the different edge functionalised graphene nanoribbons are also studied. An overall picture of edge effects, aromaticity patterns, and band gap tuning is provided.

A density functional theory (DFT) study is performed to determine the stability of the complexes formed between either the N or O site of 3-methyl-4-pyrimidone and 1-methyl-2-pyrimidone molecules and different ligands. The studied ligands are boron and alkali Lewis acids, namely, B(CH(3))(3), HB(CH(3))(2), H(2)B(CH(3)), BH(3), H(2)BF, HBF(2), BF(3), Li(+), Na(+), and K(+). The acids are divided into two groups according to their hardness. The reactivity predictions, according to the molecular electrostatic potential (MEP) map and the natural bond orbital (NBO) analysis, are in agreement with the calculated relative stabilities. Our findings reveal a strong regioselectivity with borane and its derivatives preferring the nitrogen site in both pyrimidone isomers, while a preference for oxygen is observed for the alkali acids in the 3-methyl-4-pyrimidone molecule. The complexation of 1-methyl-2-pyrimidone with these hard alkali acids does not show any discrimination between the two sites due to the presence of a continuous delocalized density region between the nitrogen and the oxygen atoms. The preference of boron Lewis acids toward the N site is due to the stronger B-N bond as compared to the B-O bond. The influence of fluorine or methyl substitution on the boron atom is discussed through natural orbital analysis (NBO) concentrating on the overlap of the boron empty p-orbital with the F lone pairs and methyl hyperconjugation, respectively. The electrophilicity of the boron acids gives a good overall picture of the interaction capabilities with the Lewis base.

Recently, hexaphyrins have emerged as a promising class of π-conjugated molecules that display a range of interesting electronic, optical, and conformational properties, including the formation of stable Mobius aromatic systems. Besides the Mobius topology, hexaphyrins can adopt a variety of conformations with Huckel and twisted Huckel topologies, which can be interconverted under certain conditions. To determine the optimum conditions for viable Mobius topologies, the conformational preferences of [26]- and [28]hexaphyrins and the dynamic interconversion between the Mobius and Huckel topologies were investigated by density functional calculations. In the absence of meso substituents, [26]hexaphyrin prefers a planar dumbbell conformation, strongly aromatic and relatively strain free. The Mobius topology is highly improbable: the most stable tautomer is 33 kcal mol(-1) higher in energy than the global minimum. On the other hand, the Mobius conformer of [28]hexaphyrin is only 6.5 kcal mol(-1) higher in energy than the most stable dumbbell conformation. This marked difference is due to aromatic stabilization in the Mobius 4n electron macrocycle as opposed to antiaromatic destabilization in the 4n+2 electron system, as revealed by several energetic, magnetic, structural, and reactivity indices of aromaticity. For [28]hexaphyrins, the computed activation barrier for interconversion between the Mobius aromatic and Huckel antiaromatic conformers ranges from 7.2 to 10.2 kcal mol(-1), in very good agreement with the available experimental data. The conformation of the hexaphyrin macrocycle is strongly dependent on oxidation state and solvent, and this feature creates a promising platform for the development of molecular switches.

We have investigated aqueous Al-dimer complexes using density functional theory methods (e.g. the B3LYP exchange-correlation functional and the 6-311++G(d,p) basis set). In these calculations interactions between the Al(3+) cations and the H(2)O or OH(-) coordinating ligands are considered explicitly while the second hydration shell and remaining solvent are treated as a continuum under the IEF-PCM formalism. The Al-dimer chemical reactivity is discussed by analysis of changes in geometry, electronic structure and Gibbs free energy of formation, relative to two independent Al(H(2)O) monomers, as a function of water and hydroxide coordination. Our results indicate that the mechanism of cooperativity, i.e. decreased Al-water bond stability with increasing OH(-) coordination and increased water ligand hydrolysis as complex CN decreases, is operating on the dimer species and that, therefore, a wide variety of dimer species are available. While the stability of these species is observed to be dependent on the number of water and hydroxide ligands, the hydroxide bridging structure (singly, doubly and triply bridged species are considered) does not appear to correlate with dimer stability. Interestingly, intra-molecular H-bonds (in the form of the well known H(3)O bridge as well as two bridging structures, H(4)O(2) and H(2)O, that have not, to our knowledge, been previously considered) are observed to influence dimer stability. The evaluation of the equilibrium mole fraction of the dimer species in equilibrium with the aqueous Al(3+) monomer species of our previous study displays the qualitatively correct trend of solution composition as pH increases, namely monomeric aqueous Al(3+) and Al(OH) complexes dominate at low and high pH, respectively, and all remaining monomer and dimer species exist at intermediate pH. Further refinement of our data set by eliminating dimer complexes with OH/Al ratios greater than 2.6 brings our predicted equilibrium mole fraction distributions into excellent agreement with experimental observations. The triply bridged dimer is observed in low amounts while the singly and doubly bridged dimers dominate our model system at pH = ∼4-7.

We analyse the electronic structure and aromaticity of graphene nanoribbons and carbon nanotubes through a series of delocalisation and geometry analysis methods. In particular, the six-centre index (SCI) is found to be in good agreement with the mean bond length (MBL) and ring bond dispersion (RBD) geometry descriptors. Based on DFT periodic calculations, three distinct classes of aromaticity patterns have been found for armchair graphene nanoribbons, appearing periodically as the width of the ribbon is increased. The periodicity in the band gap is found to be related to these aromaticity patterns. Also, the appearance of such distinct aromaticity distribution is explained within the framework of the Clar’s sextet theory. Both delocalisation and geometry analysis methods are shown to be very fast and reliable tools for easily analysing the aromaticity in carbon nanosystems.

An experimental and theoretical DFT study was carried out on the solution behavior in [D(7)]DMF for bis-chelate complex [Pd(L)(2)](BF(4))(2)center dot 2CH(3)CN (L = 4-phenyl-1-(2-picolyl)-1,2,3-triazole). In structure of [Pd(L)(2)](2+), the central square-planar palladium(II) cation is trans-chelated by two L substrates, each through the pyridine and the triazole N2 nitrogen atoms, forming two six-membered metallacycles. These can adopt boat-like conformations anti-trans-[Pd(L)(2)](2+) and syn-trans-[Pd(L)(2)](2+) in which the picolyl methylene carbons are anti or syn, respectively, relative to the palladium coordination plane. In solution, the boat-to-boat inversion at both metallacycles takes place. The conformers are in a dynamic equilibrium, which was monitored by variable-temperature (VT) (1)H NMR spectroscopy in the temperature range of 223-353 K. The equilibrium lies on the side of the anti-trans-[Pd(L)(2)](2+) conformer and the corresponding reaction enthalpy and entropy is estimated to be 0.6 +/- 0.5kcal mol(-1) and 0.8 +/- 1 cal mol(-1) K(-1), respectively. From the full-line-shape analysis of resonances in the VT (1)H NMR spectra, the activation enthalpy and activation entropy was determined to be 13.0 +/- 0.4 kcal mol(-1) and 2.7 +/- 1.6 cal mol(-1) K(-1), respectively. The activation entropy close to zero suggests a nondissociative mechanism for the isomerisation. DFT investigation revealed that the isomerisation proceeds through a one step mechanism with a barrier of 11.40 kcal mol(-1). The structures of the syn and anti conformers as well as that of the transition state were characterized. Energy decomposition analysis was carried out in order to explore the origins of the stability difference between the syn and anti isomers. (C) 2011 Elsevier Ltd. All rights reserved.

The concept of bond ionicity, obtained via a valence bond analysis, is invoked in the interpretation of the catalytic activity of supported vanadium oxides, in analogy with previous work conducted within the framework of conceptual DFT. For a set of model clusters representing the vanadium oxide supported on SiO(2), Al(2)O(3), TiO(2), ZrO(2), the ionic character of the vanadium-oxygen bond, involved in the dissociative adsorption of methanol on the catalyst, wasquantified. Detailed scrutiny shows that this ionicity increases from the Al through the Zr support, in agreement with the increasing catalytic activity through this series; the case of the Si supported oxide is found to be an exception however, giving rise to the most ionic V-O bond of the different compounds studied. This finding is confirmed by calculations on smaller clusters focusing on detail in the pi back bonding. (C) 2011 Elsevier B. V. All rights reserved.

Density functional theory calculations at the B3LYP/SVP and B3LYP/6-311G(d) levels were carried out for a series of XH(3)B(80) complexes with X = {N, P, As, B, Al}. To probe the regioselectivity of B(80), the electronic Fukui function, the molecular electrostatic potential (MEP), and the natural bond orbital (NBO) were determined. These indices were shown to provide reliable guides to predict the relative reactivities of the boron buckyball sites. Thermodynamic stabilities of the complexes formed by the reaction of B(80) with nucleophiles (NH(3), PH(3), AsH(3)) and electrophiles (BH(3), AlH(3)) are in good agreement with the prediction of regioselectivity indicated on the basis of Fukui and MEP indices. The qualitative results suggest the boron buckyball to be an amphoteric and hard molecule. It has two distinct reactive sites localized on caps and frame, which act as acids and bases, respectively. Most of the complexes are stable with formation energies comparable to that of the analogous complexes of the borane molecule, BH(3)BH(3), BH(3)NH(3), and BH(3)AlH(3). The B-H-B bond characteristics of diborane are recovered in B(80)BH(3). Exohedral complexes are more stable than endohedral complexes. The most stable complexes are those with NH(3) on the caps and BH(3) on the pentagonal ring of B(80).

The S E Ar reaction was scrutinized using a quantitative, fragment-based interaction energy decomposition analysis (EDA) as implemented in the ADF program. The interaction energy between monosubstituted benzene derivatives and a model electrophile at the onset of the reaction was studied and decomposed into Pauli repulsion, electrostatic interaction and orbital interaction terms for a plane parallel to the molecular plane. As such experimentally observed selectivity patterns arise in the orbital interaction, stressing again the role of both electrophile and benzene derivative in determining the regioselectivity. Using the HSAB principle and MO theory, the orientation mechanism of the SEAr is further explored.

Three intrinsic group properties, the group electronegativity, hardness, and softness, were obtained for 30 organic groups in a nonempirical way starting from the definitions within the framework of density functional theory and applying limited configuration interaction (CISD). The calculated group electronegativities correlatcd well with most scales introduced for this quantity. The values were compared with the field factor F and with the 13 C 1 J CC (ipso-ortho) coupling constant in monosubstituted benzenes, the most important experimental quantity known to reflect substituent electronegativities

This letter is concerned with the calculation and interpretation of the atom-condensed linear response kernel as defined within conceptual density functional theory. An easily applicable and accurate methodology following from second order perturbation theory is used to analyse a series of cyclohexane and benzene derivatives. We show that the variation of the condensed linear response matrix elements with internuclear distance provides a means to differentiate between inductive, resonance and hyperconjugation effects and conclude that this quantity measures the extent of electron delocalisation.

Experimentalists and quantum chemists are living in a different world. A wealth of theoretical enzymology-related publications is hardly known by experimentalists, and vice versa. Our aim is to bring both worlds together and to show the powerful possibilities of a multidisciplinary approach to study subtle details of complicated enzymatic processes to a broad readership. MD simulations and QM/MM approaches often focus on the calculation of reaction paths based on activation energies, which is a time-consuming task. A valuable alternative is the reactivity descriptors founded in conceptual DFT like softness, electrophilicity, and the Fukui function, which describe the kinetic aspects of a reaction in terms of the response to perturbations in N and/or upsilon(r), typical for a chemical reaction, of the reagents in the ground state. As such, the relative energies at the beginning of the reaction predict a sequence of activation energies only based on the properties of the reactants (Figure 5 ). In 2003, Geerlings et al. published a key review giving a detailed description of the principles and concepts of conceptual DFT and highlighting its success to study generalized acid/base reactions including addition, substitution, and elimination reactions. Since the time that this review appeared, conceptual DFT has proven its strength in literally hundreds of papers with application to organic and inorganic reactions. Its role in unravelling enzymatic reaction mechanisms, in handling experimentally difficult accessible biochemical problems, and in the interpretation of biochemical experimental observations is emerging and very promising.

Although extensive theoretical and experimental research has been conducted on fluorinated fullerenes, little detailed information exists on their solubility in different solvents. However, this solubility is crucial for their processability and possible application. In this work, we predict the solubility of fluorinated C(60) in various polar and non-polar solvents, based on a correlation between experimentally measured solubilities for C(60) from the literature and theoretically predicted solubilisation energies. These solubilisation energies are predicted using the polarizable continuum model (PCM) at the ab initio Hartree-Fock 6-31+G* level of theory. In particular, the solubilities are predicted for C(60)F(2)(n) (n = 1-10) isomers, part of the addition route to saturnene C(60)F(20). An increasing solubility is found for more polar solvents with higher degree of fluorination. With these results, we can determine the minimal fluorination degree necessary for possible solubilisation in a given solvent, and offer new perspectives to separate and purify species with different degrees of fluorination.