Search

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 initio calculations lead to halogen- and hydrogen-bonded complexes. Fourier transform-IR experiments involving solutions of mixtures of difluoroiodomethane with trimethylphosphine(-d 9 ) or dimethyl sulfide(-d 6 ) 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.

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.

Dehydroascorbate reductase (DHAR) catalyzes the glutathione (GSH)-dependent reduction of dehydroascorbate and plays a direct role in regenerating ascorbic acid, an essential plant antioxidant vital for defense against oxidative stress. DHAR enzymes bear close structural homology to the glutathione transferase (GST) superfamily of enzymes and contain the same active site motif, but most GSTs do not exhibit DHAR activity. The presence of a cysteine at the active site is essential for the catalytic functioning of DHAR, as mutation of this cysteine abolishes the activity. Here we present the crystal structure of DHAR2 from Arabidopsis thaliana with GSH bound to the catalytic cysteine. This structure reveals localized conformational differences around the active site which distinguishes the GSH-bound DHAR2 structure from that of DHAR1. We also unraveled the enzymatic step in which DHAR releases oxidized glutathione (GSSG). To consolidate our structural and kinetic findings, we investigated potential conformational flexibility in DHAR2 by normal mode analysis and found that subdomain mobility could be linked to GSH binding or GSSG release., Competing Interests: The authors declare no competing financial interests.

Electrophilic aromatic sulfonation of benzene with sulfur trioxide is studied with ab initio molecular dynamics simulations in gas phase, and in explicit noncomplexing (CCl 3 F) and complexing (CH 3 NO 2 ) solvent models. We investigate different possible reaction pathways, the number of SO 3 molecules participating in the reaction, and the influence of the solvent. Our simulations confirm the existence of a low-energy concerted pathway with formation of a cyclic transition state with two SO 3 molecules. Based on the simulation results, we propose a sequence of elementary reaction steps and a kinetic model compatible with experimental data. Furthermore, a new alternative reaction pathway is proposed in complexing solvent, involving two SO 3 and one CH 3 NO 2 .

The interactions of axial substituents in monosubstituted cyclohexane rings are studied in this work using an array of different computational techniques. Additionally, the anomalous axial preference for some bulky substituents is related to stabilizing dispersion interactions. We find that the C-H ax ···Y ax contacts for various substituents with distances ranging from 2 to ∼5 Å may include attractive dispersion forces that can affect the conformational equilibrium; these forces co-exist with Pauli repulsive forces effected by Y ax group due to van der Waals sphere penetration. At distances between 2 and 3 Å stabilizing electron transfer interactions were calculated and the combination of natural bond orbital and QTAIM analysis showed that, in certain cases, Y ax  =  t Bu, C ax -O or C ax  = O or S ax  = O or C ax  = S this interaction can be characterized as an improper H-bond. DFT-D3 and non-covalent interactions calculations (NCIs) in cyclohexane derivatives with Y ax  = SiOR 3 including H Yax ···H cy surfaces at distances ranging between 4 and 6 Å suggest that dispersion has a clear effect on the experimentally observed stabilization of the axial conformer. NCIs computed from the reduced density gradient help to visually identify and analyze these interactions. © 2016 Wiley Periodicals, Inc., (© 2016 Wiley Periodicals, Inc.)

The range-separation parameter in tuned, range-separated exchange-correlation functionals is investigated in two contexts. First, the system dependence of the parameter is investigated for a series of systems obtained by successively ionizing a single species, paying particular attention to the degree of linearity in the energy versus electron number curve. The parameter exhibits significant system dependence and, therefore, achieving near-linearity in one segment of the curve leads to strong nonlinearity in other segments. This provides a challenging test case for the development of new functionals designed to overcome the known problems of this class of functional. Next, the study considers whether a range-separation parameter tuned to a Koopmans energy condition is also applicable for the analogous density condition. This is tested by comparing two formulations of the Fukui function of conceptual density functional theory, for three representative systems. Both formulations yield the same general features and are not highly sensitive to the range-separation parameter. However, the agreement between the two is near-optimal when the energy-tuned parameter is used, indicating that this parameter is applicable for the analogous density condition.

In many established methods, identification of hydrogen bonds (H-bonds) is primarily based on pairwise comparison of distances between atoms. These methods often give rise to systematic errors when sulfur is involved. A more accurate method is the non-covalent interaction index, which determines the strength of the H-bonds based on the associated electron density and its gradient. We applied the NCI index on the active site of a single-cysteine peroxiredoxin. We found a different sulfur hydrogen-bonding network to that typically found by established methods, and we propose a more accurate equation for determining sulfur H-bonds based on geometrical criteria. This new algorithm will be implemented in the next release of the widely-used CHARMM program (version 41b), and will be particularly useful for analyzing water molecule-mediated H-bonds involving different atom types. Furthermore, based on the identification of the weakest sulfur-water H-bond, the location of hydrogen peroxide for the nucleophilic attack by the cysteine sulfur can be predicted. In general, current methods to determine H-bonds will need to be reevaluated, thereby leading to better understanding of the catalytic mechanisms in which sulfur chemistry is involved.

There is an increasing interest in applying quantum chemistry to rationally design novel compounds with some desired characteristics. Furthermore, many applications require more than one property to be optimal. In this Concept, several inverse design strategies, based on the discrete best first search scheme, are introduced that allow for the simultaneous optimization of multiple properties or the optimization of the most vital target property with constraints for secondary properties. A detailed assessment of the different optimization techniques is carried out, and special attention is paid to improve the cost efficacy and performance by tuning the process parameters. Our suggested protocol allows for a more successful optimization routine when additional boundary conditions are desired., (© 2016 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.)

Several expanded porphyrins switch between Hückel, Möbius and twisted-Hückel topologies, encoding different aromaticity and NLO properties. Since the topological switch can be induced by different external stimuli, expanded porphyrins represent a promising platform to develop molecular switches for molecular electronic devices. In order to determine the optimum conditions for efficient molecular switches from octaphyrins, we have carried out a comprehensive quantum chemical study focusing on the conformational preferences and aromaticity of [36]octaphyrins. Different external stimuli for triggering the topological switch have been considered in our work, such as protonation and redox reactions. Importantly, the structure-property relationships between the molecular conformation, the number of π-electrons and aromaticity in octaphyrins have been established by using energetic, magnetic, structural and reactivity descriptors. Remarkably, we found that the aromaticity of octaphyrins is highly dependent on the π-conjugation topology and the number of π-electrons and it can be modulated by protonation and redox reactions. A non-aromatic figure-eight conformation is strongly preferred by neutral [36]octaphyrins that switches to a Möbius aromatic conformation upon protonation. Such a change of topology involves an aromaticity switch in a single molecule and is accompanied by a drastic change in the NLO properties. In contrast, the twisted-Hückel topology remains the most stable one in the oxidized and reduced species, but the aromaticity is totally reversed upon redox reactions. Aromaticity is shown to be a key concept in expanded porphyrins, determining the electronic, magnetic and NLO properties of these macrocycles.

Our attempts to synthesise N→M intramolecularly coordinated diorganometallic hydrides L2MH2 [M=Si (4), Ge (5), Sn (6)] containing the CH=N imine group (in which L is C,N-chelating ligand {2-[(2,6-iPr2C6H3)N=CH]C6 H4}(-)) yielded 1,1'-bis(2,6-diisopropylphenyl)-2,2'-spriobi[benzo[c][1,2]azasilole] (7), 1,1'-bis(2,6-diisopropylphenyl)-2,2'-spriobi[benzo[c][1,2]azagermole] (8) and C,N-chelated homoleptic stannylene L2Sn (10), respectively. Compounds 7 and 8 are an outcome of a spontaneous double hydrometallation of the two CH=N imine moieties induced by N→M intramolecular coordination (M=Si, Ge) in the absence of any catalyst. In contrast, the diorganotin hydride L2SnH2 (6) is redox-unstable and the reduction of the tin centre with the elimination of H2 provided the C,N-chelated homoleptic stannylene L2Sn (10). Compounds 7 and 8 were characterised by NMR spectroscopy and X-ray diffraction analysis. Because the proposed N→M intramolecularly coordinated diorganometallic hydrides L2MH2 [M=Si (4), Ge (5), Sn (6)] revealed two different types of reduction reactions, DFT calculations were performed to gain an insight into the structures and bonding of the non-isolable diorganometallic hydrides as well as the products of their subsequent reactions. Furthermore, the thermodynamic profiles of the different reaction pathways with respect to the central metal atom were also investigated., (© 2016 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.)

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.

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.

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., (© 2015 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.)

The reduction of N,C,N-chelated bismuth chlorides [C6H3-2,6-(CH=NR)2]BiCl2 [where R = tBu (1), 2',6'-Me2C6H3 (2), or 4'-Me2NC6H4 (3)] or N,C-chelated analogues [C6H2-2-(CH=N-2',6'-iPr2C6H3)-4,6-(tBu)2]BiCl2 (4) and [C6H2-2-(CH2NEt2)-4,6-(tBu)2]BiCl2 (5) is reported. Reduction of compounds 1-3 gave monomeric N,C,N-chelated bismuthinidenes [C6H3-2,6-(CH=NR)2]Bi [where R = tBu (6), 2',6'-Me2C6H3 (7) or 4'-Me2NC6H4 (8)]. Similarly, the reduction of 4 led to the isolation of the compound [C6H2-2-(CH=N-2',6'-iPr2C6H3)-4,6-(tBu)2]Bi (9) as an unprecedented two-coordinated bismuthinidene that has been structurally characterized. In contrast, the dibismuthene {[C6H2-2-(CH2NEt2)-4,6-(tBu)2]Bi}2 (10) was obtained by the reduction of 5. Compounds 6-10 were characterized by using (1)H and (13)C NMR spectroscopy and their structures, except for 7, were determined with the help of single-crystal X-ray diffraction analysis. It is clear that the structure of the reduced products (bismuthinidene versus dibismuthene) is ligand-dependent and particularly influenced by the strength of the N→Bi intramolecular interaction(s). Therefore, a theoretical survey describing the bonding situation in the studied compounds and related bismuth(I) systems is included. Importantly, we found that the C3NBi chelating ring in the two-coordinated bismuthinidene 9 exhibits significant aromatic character by delocalization of the bismuth lone pair., (© 2015 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.)

Two ethylene-bridged methylaluminium amidinates and one aluminium amidinate containing three terminal trimethylstannyl-ethynyl groups interconnected by π-coordinated potassium ions were prepared in situ. The re-oxidation of the ethylene-bridged compound by iodine followed by further reduction using the same activation procedure demonstrated the versatility of the approach. The reactivity of an ethylene-bridged methylaluminum amidinate towards HCl was examined to demonstrate the building block concept. DFT calculations were performed to gain insight into the mechanism of the in situ activation of diphenylacetylene.

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 molecular electrostatic potential surface of unsaturated, locally electron-deficient molecules shows a positive region perpendicular to (a part of) the molecular framework. In recent years it has been shown both theoretically and experimentally that molecules are able to form noncovalent interactions with Lewis bases through this π-hole. When studying unsaturated perfluorohalogenated molecules containing a higher halogen atom, a second electropositive region is also observed near the halogen atom. This region, often denoted as a σ-hole, allows the molecules to interact with Lewis bases and form a halogen bond. To experimentally characterize the competition between both these noncovalent interactions, Fourier transform infrared and Raman spectra of liquefied noble gas solutions containing perfluorohalogenated ethylene derivatives (C2F3X; X = F, Cl, Br, or I) and trimethylamine(-d9) were investigated. Analysis of the spectra shows that in mixed solutions of trimethylamine(-d9) and C2F4 or C2F3Cl lone pair···π complex is present, while evidence for halogen-bonded complex is found in solutions containing trimethylamine(-d9) and C2F3Cl, C2F3Br, or C2F3I. For all species observed, complexation enthalpies were determined, the values varying between -4.9(1) and -24.4 kJ mol(-1).

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

The contributions of covalent and noncovalent interactions to the formation of classical adducts of bulky Lewis acids and bases and frustrated Lewis pairs (FLPs) were scrutinized by using various conceptual quantum chemical techniques. Significantly negative complexation energies were calculated for fourteen investigated Lewis pairs containing bases and acids with substituents of various sizes. A Ziegler-Rauk-type energy decomposition analysis confirmed that two types of Lewis pairs can be distinguished on the basis of the nature of the primary interactions between reactants; dative-bond formation and concomitant charge transfer from the Lewis base to the acid is the dominant and most stabilizing factor in the formation of Lewis acid-base adducts, whereas weak interactions are the main thermodynamic driving force (>50 %) for FLPs. Moreover, the ease and extent of structural deformation of the monomers appears to be a key component in the formation of the former type of Lewis pairs. A Natural Orbital for Chemical Valence (NOCV) analysis, which was used to visualize and quantify the charge transfer between the base and the acid, clearly showed the importance and lack of this type of interaction for adducts and FLPs, respectively. The Noncovalent Interaction (NCI) method revealed several kinds of weak interactions between the acid and base components, such as dispersion, π-π stacking, C-H⋅⋅⋅π interaction, weak hydrogen bonding, halogen bonding, and weak acid-base interactions, whereas the Quantum Theory of Atoms in Molecules (QTAIM) provided further conceptual insight into strong acid-base interactions., (© 2015 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.)

Atomic charges are a key concept to give more insight into the electronic structure and chemical reactivity. The Hirshfeld-I partitioning scheme applied to the model protein human 2-cysteine peroxiredoxin thioredoxin peroxidase B is used to investigate how large a protein fragment needs to be in order to achieve convergence of the atomic charge of both neutral and negatively charged residues. Convergence in atomic charges is rapidly reached for neutral residues, but not for negatively charged ones. This study pinpoints difficulties on the road toward accurate modeling of negatively charged residues of large biomolecular systems in a multiscale approach.

Inspection of the electrostatic potential of C2F3X (X = F, Cl, Br, and I) revealed a second electropositive region in the immediate vicinity of the C═C double bond apart from the σ hole of chlorine, bromine, and iodine, leading to C(sp(2))-X···Y halogen bonding, through which complexes stabilized by so-called lone pair···π interactions can be formed. Consequently, the experimental studies for the complexes of dimethyl ether with C2F3X (X = F, Cl, Br, and I) not only allowed one to experimentally characterize and rationalize the effects of hybridization on halogen bonding but, for the first time, also allowed the competition of C-X···Y halogen bonding and lone pair···π interactions to be studied at thermodynamic equilibrium. Analysis of the infrared and Raman spectra reveals that in the cryosolutions of dimethyl ether and C2F3I, solely the halogen-bonded complex is present, whereas C2F3Br and C2F3Cl give rise to a lone pair···π bonded complex as well as a halogen-bonded complex. Mixtures of dimethyl ether with C2F4 solely yield a lone pair···π bonded complex. The experimentally derived complexation enthalpies for the halogen bonded complexes are found to be -14.2(5) kJ mol(-1) for C2F3I·DME and -9.3(5) kJ mol(-1) for C2F3Br·DME. For the complexes of C2F3Cl with dimethyl ether, no experimental complexation enthalpy could be obtained, whereas the C2F4·DME complex has a complexation enthalpy of -5.5(3) kJ mol(-1). The observed trends have been rationalized with the aid of an interaction energy decomposition analysis (EDA) coupled to a Natural Orbital for Chemical Valence (NOCV) analysis and also using the noncovalent interaction index method.

The reactivity of three phosphanes (PhP(NHR)2 [R = t-Bu (1), Ph (2)] and PhP(NEt2)(NHDip) (3) (where Dip = 2,6-i-Pr2C6H3)) with n-Bu2Mg and MeMgBr is presented. In the case of 1, the reaction with n-Bu2Mg gave [PhP(NHt-Bu)(Nt-Bu)]Mg(n-Bu) (4) or [PhP(NHt-Bu)(Nt-Bu)]2Mg (5) depending on the stoichiometry. The treatment of 1 with MeMgBr led to the phosphinate [Ph(H)P(Nt-Bu)2]2Mg (7) as a result of both the NH→PH tautomeric transformation and elimination of MgBr2 from the non-isolable intermediate [PhP(NHt-Bu)(Nt-Bu)]MgBr(THF) (6). Phosphane 2 reacted with n-Bu2Mg in a 1:1 molar ratio resulting in the formation of {[PhP(NPh)2]2Mg(THF)2}2 (8), but the analogous reaction in a 2:1 molar ratio yielded phosphinate [Ph(H)P(NPh)2]2Mg(THF) (9). The heteroleptic compound [Ph(H)P(NPh)2]MgBr(THF)2 (10) was obtained by the reaction of 2 with MeMgBr. Finally, the reaction of 3 with n-Bu2Mg and MeMgBr produced compounds [PhP(NEt2)(NDip)]2Mg (11) and {[PhP(NEt2)(NDip)2]Mg(μ-Br)(THF)}2 (12), respectively. All products were characterized using (1)H, (13)C{(1)H} and (31)P NMR spectroscopy and, except for compounds 4 and 6, their molecular structures were determined using single-crystal X-ray diffraction analysis. In addition, a theoretical study on plausible isomers of 10 was performed to provide additional evidence for the presence of a syn- and anti-isomer in dynamic equilibrium in a solution of 10.

The atom-atom polarizability and the transmission probability at the Fermi level, as obtained through the source-and-sink-potential method for every possible configuration of contacts simultaneously, are compared for polycyclic aromatic compounds. This comparison leads to the conjecture that a positive atom-atom polarizability is a necessary condition for transmission to take place in alternant hydrocarbons without non-bonding orbitals and that the relative transmission probability for different configurations of the contacts can be predicted by analyzing the corresponding atom-atom polarizability. A theoretical link between the two considered properties is derived, leading to a mathematical explanation for the observed trends for transmission based on the atom-atom polarizability.

The establishment of an internally consistent scale of noble gas electronegativities is a long-standing problem. In the present study, the problem is attacked via the Mulliken definition, which in recent years gained widespread use to its natural appearance in the context of conceptual density functional theory. Basic ingredients of this scale are the electron affinity and the ionization potential. Whereas the latter can be computed routinely, the instability of the anion makes the judicious choice of computational technique for evaluating electron affinities much more tricky. We opted for Puiatti's approach, extrapolating the energy of high ε solvent stabilized anions to the ε = 1 (gas phase) case. The results give negative electron affinity values, monotonically increasing (except for helium which is an outlier in most of the story) to almost zero at eka-radon in agreement with high level calculations. The stability of the B3LYP results is successfully tested both via improving the level of theory (CCSD(T)) and expanding the basis set. Combined with the ionization energies (in good agreement with experiment), an electronegativity scale is obtained displaying (1) a monotonic decrease of χ when going down the periodic table, (2) top values not for the noble gases but for the halogens, as opposed to most (extrapolation) procedures of existing scales, invariably placing the noble gases on top, and (3) noble gases having electronegativities close to the chalcogens. In the accompanying hardness scale (hardly, if ever, discussed in the literature) the noble gases turn out to be by far the farthest the hardest elements, again with a continuous decrease with increasing Z. Combining χ value of the halogens and the noble gases the Ng(δ+)F(δ-) bond polarity emerging from ab initio calculations naturally emerges. In conclusion, the chemistry of the noble gases is for a large part determined by their extreme hardness, equivalent to a high resistance to change in its electronic population coupled to their high electronegativity.

Novel spiropseudoindoxyls were synthesized in high yields via a fully regioselective Au(III)-catalyzed cycloisomerization of easily obtainable o-nitrophenylpropiolamides, followed by an intramolecular dipolar cycloaddition as key steps. This one-pot cascade reaction resulted in new tetracyclic pseudoindoxyls, which were hydrogenated toward the title compounds as single diastereomers via N-O cleavage. The mechanism of the gold catalyzed cycloisomerization was studied by DFT and suggests a reaction path without the intermediacy of gold carbenoid species in these cases.

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.

Within the context of spin polarized conceptual density functional theory, the spin polarized linear response functions are introduced both in the [N, N(s)] and [N(α), N(β)] representations. The mathematical relations between the spin polarized linear response functions in both representations are examined and an analytical expression for the spin polarized linear response functions in the [N(α), N(β)] representation is derived. The spin polarized linear response functions were calculated for all atoms up to and including argon. To simplify the plotting of our results, we integrated χ(r, r') to a quantity χ(r, r'), circumventing the θ and ϕ dependence. This allows us to plot and to investigate the periodicity throughout the first three rows in the periodic table within the two different representations. For the first time, χ(αβ)(r, r'), χ(βα)(r, r'), and χ(SS)(r, r') plots have been calculated and discussed. By integration of the spin polarized linear response functions, different components to the polarisability, α(αα), α(αβ), α(βα), and α(ββ) have been calculated.

Hydrogen peroxide is a natural oxidant that can oxidize protein thiols (RSH) via sulfenic acid (RSOH) and sulfinic acid (RSO2H) to sulfonic acid (RSO3H). In this paper, we study the complete anionic and neutral oxidation pathway from thiol to sulfonic acid. Reaction barriers and reaction free energies for all three oxidation steps are computed, both for the isolated substrates and for the substrates in the presence of different model ligands (CH4, H2O, NH3) mimicking the enzymatic environment. We found for all three barriers that the anionic thiolate is more reactive than the neutral thiol. However, the assistance of the environment in the neutral pathway in a solvent-assisted proton-exchange (SAPE) mechanism can lower the reaction barrier noticeably. Polar ligands can decrease the reaction barriers, whereas apolar ligands do not influence the barrier heights. The same holds for the reaction energies: they decrease (become more negative) in the presence of polar ligands whereas apolar ligands do not have an influence. The consistently negative consecutive reaction energies for the oxidation in the anionic pathway when going from thiolate over sulfenic and sulfinic acid to sulfonic acid are in agreement with biological reversibility.

Most of the work done on the linear response kernel χ(r,r') has focussed on its atom-atom condensed form χAB. Our previous work [Boisdenghien et al., J. Chem. Theory Comput., 2013, 9, 1007] was the first effort to truly focus on the non-condensed form of this function for closed (sub)shell atoms in a systematic fashion. In this work, we extend our method to the open shell case. To simplify the plotting of our results, we average our results to a symmetrical quantity χ(r,r'). This allows us to plot the linear response kernel for all elements up to and including argon and to investigate the periodicity throughout the first three rows in the periodic table and in the different representations of χ(r,r'). Within the context of Spin Polarized Conceptual Density Functional Theory, the first two-dimensional plots of spin polarized linear response functions are presented and commented on for some selected cases on the basis of the atomic ground state electronic configurations. Using the relation between the linear response kernel and the polarizability we compare the values of the polarizability tensor calculated using our method to high-level values.

N-fused pentaphyrins (NFP) are the stable forms of fully meso-aryl pentaphyrins(1.1.1.1.1). In order to determine the optimum conditions for viable Möbius topologies of these porphyrinoids, the conformational preferences, Hückel-Möbius interconversion pathways and aromaticity of [22] and [24]NFP have been investigated using density functional theory calculations. The conformation of the macrocycle is shown to be strongly dependent on the oxidation state and the macrocyclic aromaticity. [22]NFP prefers a highly aromatic and relatively strain-free Hückel conformation. However, antiaromatic Hückel and weakly aromatic Möbius conformers coexist in dynamic equilibrium in [24]NFP. The Hückel-Möbius aromaticity switch requires very low activation energy barriers (Ea = 3-4 kcal mol(-1)). Interestingly, the balance between Möbius and Hückel conformations in [24]NFP can be controlled by meso-substituents. The structure-property relationship between the molecular conformation, number of π electrons and aromaticity has been established in our study using energetic, magnetic, structural, and reactivity descriptors of aromaticity. Although the Möbius topology is indeed accessible for [24]NFP, it does not exhibit a distinct macrocyclic aromaticity mainly due to the large dihedral angles around the molecular twist. Regarding the computational methodology, B3LYP and M06 show the best overall performance for describing the experimental geometries of NFP and, importantly, our computational results support the experimental evidence available for N-fused pentaphyrins.

The essential aspects of zero-temperature grand-canonical ensemble density-functional theory are reviewed in the context of spin-density-functional theory and are used to highlight the assumption of symmetry between electron addition and subtraction that underlies the corrected Koopmans approach of Tozer and De Proft (TDP) for computing electron affinities. The issue of symmetry is then investigated in a systematic study of atomic electron affinities, comparing TDP affinities with those from a conventional Koopmans evaluation and electronic energy differences. Although it cannot compete with affinities determined from energy differences, the TDP expression yields results that are a significant improvement over those from the conventional Koopmans expression. Key insight into the results from both expressions is provided by an analysis of plots of the electronic energy as a function of the number of electrons, which highlight the extent of symmetry between addition and subtraction. The accuracy of the TDP affinities is closely related to the nature of the orbitals involved in the electron addition and subtraction, being particularly poor in cases where there is a change in principal quantum number, but relatively accurate within a single manifold of orbitals. The analysis is then extended to a consideration of the ground state Mulliken electronegativity and chemical hardness. The findings further emphasize the key role of symmetry in determining the quality of the results.

Within the context of reactivity descriptors known in conceptual DFT, the linear response function (χ(r,r')) remained nearly unexploited. Although well known, in its time dependent form, in the solid state physics and time-dependent DFT communities the study of the "chemistry" present in the kernel was, until recently, relatively unexplored. The evaluation of the linear response function as such and its study in the time independent form are highlighted in the present review. On the fundamental side, the focus is on the approaches of increasing complexity to compute and represent χ(r,r'), its visualisation going from plots of the unintegrated χ(r,r') to an atom condensed matrix. The study on atoms reveals its physical significance, retrieving atomic shell structure, while the results on molecules illustrate that a variety of chemical concepts are retrieved: inductive and mesomeric effects, electron delocalisation, aromaticity and anti-aromaticity, σ and π aromaticity,…. The applications show that the chemistry of aliphatic (saturated and unsaturated) chains, saturated and aromatic/anti-aromatic rings, organic, inorganic or metallic in nature, can be retrieved via the linear response function, including the variation of the electronic structure of the reagents along a reaction path. The connection of the linear response function with the concept of nearsightedness and the alchemical derivatives is also highlighted.

Novel tricyclic tetrahydroazepinones were synthesized via an in situ Diels-Alder reaction of furan with cyclic allenamides. These reactive intermediates are the first examples of cyclic seven-membered allenamides and were prepared starting from N-(2-chloroallyl)-2-allylglycine derivatives via ring-closing metathesis followed by dehydrochlorination. The trapping of the intermediate cycloallene with furan occurred endo- and regioselectively and provided a convenient entry into new building blocks for medicinal chemistry. The diastereoselectivity of the cycloaddition was confirmed using quantum chemical computations.

Insight into the key factors driving the competition of halogen and hydrogen bonds is obtained by studying the affinity of the Lewis bases trimethylamine (TMA), dimethyl ether (DME), and methyl fluoride (MF) towards difluoroiodomethane (CHF(2) I). Analysis of the infrared and Raman spectra of solutions in liquid krypton containing mixtures of TMA and CHF(2) I and of DME and CHF(2) I reveals that for these Lewis bases hydrogen and halogen-bonded complexes appear simultaneously. In contrast, only a hydrogen-bonded complex is formed for the mixtures of CHF(2) I and MF. The complexation enthalpies for the C-H⋅⋅⋅Y hydrogen-bonded complexes with TMA, DME, and MF are determined to be -14.7(2), -10.5(5) and -5.1(6) kJ mol(-1), respectively. The values for the C-I⋅⋅⋅Y halogen-bonded isomers are -19.0(3) kJ mol(-1) for TMA and -9.9(8) kJ mol(-1) for DME. Generalization of the observed trends suggests that, at least for the bases studied here, softer Lewis bases such as TMA favor halogen bonding, whereas harder bases such as MF show a substantial preference for hydrogen bonding., (© 2014 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.)

Novel organic-inorganic hybrids of various sizes were generated by reaction of 1,8-octanediphosphonic acid (ODP) and (NH4)6Mo7O24 in aqueous solution. The formation of rodlike hybrids with variable numbers of covalently bound ODP and polyoxomolybdate (POM) units can be tuned as a function of increasing (NH4)6Mo7O24 concentration at fixed ODP concentration. The chemical structure of the ODP/POM hybrids was characterized by (1)H, (31)P, and (95)Mo NMR spectroscopy. Heteronuclear (31)P DOSY (diffusion- ordered NMR spectroscopy) and molecular mechanics (MM) calculations were applied to determine the size and shape of the nanosized hybrids generated at various ODP/POM ratios. For this purpose, the structures of ODP/POM hybrids with variable numbers of ODP and POM units were optimized by MM and then approximated as cylinder-shaped objects by using a recently described mathematical algorithm. The thus-obtained cylinder length and diameter were further used to calculate the expected diffusion coefficients of the ODP/POM hybrids. Comparison of the calculated and experimentally determined diffusion coefficients led to the most probable ODP/POM hybrid length for each sample composition. The (31)P DOSY results show that the length of the hybrids increases with increasing POM concentration and reaches a maximum corresponding to an average of 8 ODP/7 POM units per chain at a sample composition of 20 mM ODP and 14 mM POM. With excess POM, above the latter concentration, the formation of shorter-chain hybrids terminated by Mo7 clusters at one or both ends was evidenced on further increasing the POM concentration. The results demonstrate that the combination of (31)P DOSY and MM, although virtually unexplored in POM chemistry, is a powerful innovative strategy for the detailed characterization of nanosized organic-inorganic POM-based hybrids in solution., (© 2014 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.)

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., (© 2014 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.)

We present the results of a theoretical investigation focusing on the solvent structure surrounding the -1, 0 and +1 charged species of F, Cl, Br and I halogen atoms and F2, Cl2, Br2 and I2 di-halogen molecules in a methanol solvent and its influence on the electronic structure of the solute molecules. Our results show a large stabilizing effect arising from the solute-solvent interactions. Well-formed first solvation shells are observed for all species, the structure of which is strongly influenced by the charge of the solute species. Detailed analysis reveals that coordination number, CN, solvent orientation, θ, and solute-solvent distance, d, are important structural characteristics which are coupled to changes in the electronic structure of the solute. We propose that the fundamental chemistry of any solute species is generally regulated by these solvent degrees of freedom.

Our attempts to synthesize the N→Si intramolecularly coordinated organosilanes Ph2 L(1) SiH (1 a), PhL(1) SiH2 (2 a), Ph2 L(2) SiH (3 a), and PhL(2) SiH2 (4 a) containing a CH=N imine group (in which L(1) is the C,N-chelating ligand {2-[CH=N(C6 H3 -2,6-iPr2)]C6 H4}(-) and L(2) is {2-[CH=N(tBu)]C6 H4}(-)) yielded 1-[2,6-bis(diisopropyl)phenyl]-2,2-diphenyl-1-aza-silole (1), 1-[2,6-bis(diisopropyl)phenyl]-2-phenyl-2-hydrido-1-aza-silole (2), 1-tert-butyl-2,2-diphenyl-1-aza-silole (3), and 1-tert-butyl-2-phenyl-2-hydrido-1-aza-silole (4), respectively. Isolated organosilicon amides 1-4 are an outcome of the spontaneous hydrosilylation of the CH=N imine moiety induced by N→Si intramolecular coordination. Compounds 1-4 were characterized by NMR spectroscopy and X-ray diffraction analysis. The geometries of organosilanes 1 a-4 a and their corresponding hydrosilylated products 1-4 were optimized and fully characterized at the B3LYP/6-31++G(d,p) level of theory. The molecular structure determination of 1-3 suggested the presence of a Si-N double bond. Natural bond orbital (NBO) analysis, however, shows a very strong donor-acceptor interaction between the lone pair of the nitrogen atom and the formal empty p orbital on the silicon and therefore, the calculations show that the Si-N bond is highly polarized pointing to a predominantly zwitterionic Si(+) N(-) bond in 1-4. Since compounds 1-4 are hydrosilylated products of 1 a-4 a, the free energies (ΔG298), enthalpies (ΔH298), and entropies (ΔH298) were computed for the hydrosilylation reaction of 1 a-4 a with both B3LYP and B3LYP-D methods. On the basis of the very negative ΔG298 values, the hydrosilylation reaction is highly exergonic and compounds 1 a-4 a are spontaneously transformed into 1-4 in the absence of a catalyst., (Copyright © 2014 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.)

Homobimetallic metallophilic interactions between copper, silver, and gold-based [(NHC)MX]-type complexes (NHC=N-heterocyclic carbene, i.e, 1,3,4-trimethyl-4,5-dihydro-1H-1,2,4-triazol-5-ylidene; X=F, Cl, Br, I) were investigated by means of ab initio interaction energies, Ziegler-Rauk-type energy-decomposition analysis, the natural orbital for chemical valence (NOCV) framework, and the noncovalent interaction (NCI) index. It was found that the dimers of these complexes predominantly adopt a head-to-tail arrangement with typical M⋅⋅⋅M distance of 3.04-3.64 Å, in good agreement with the experimental X-ray structure determined for [{(NHC)AuCl}2 ], which has an Au⋅⋅⋅Au distance of 3.33 Å. The interaction energies between silver- and gold-based monomers are calculated to be about -25 kcal mol(-1) , whereas that for the Cu congener is significantly lower (-19.7 kcal mol(-1) ). With the inclusion of thermal and solvent contributions, both of which are destabilizing, by about 15 and 8 kcal mol(-1) , respectively, an equilibrium process is predicted for the formation of dimer complexes. Energy-decomposition analysis revealed a dominant electrostatic contribution to the interaction energy, besides significantly stabilizing dispersion and orbital interactions. This electrostatic contribution is rationalized by NHC(δ(+) )⋅⋅⋅halogen(δ(-) ) interactions between monomers, as demonstrated by electrostatic potentials and derived charges. The dominant NOCV orbital indicates weakening of the π backdonation in the monomers on dimer formation, whereas the second most dominant NOCV represents an electron-density deformation according to the formation of a very weak M⋅⋅⋅M bond. One of the characteristic signals found in the reduced density gradient versus electron density diagram corresponds to the noncovalent interactions between the metal centers of the monomers in the NCI plots, which is the manifestation of metallophilic interaction., (Copyright © 2014 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.)

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 trans effect and trans influence were investigated and rationalized in the aminolysis, a typical nucleophilic substitution reaction, of trans-TPtCl2NH3 complexes (T = NH3, PH3, CO and C2H4) using energy decomposition analysis, both along the reaction paths and on the stationary points, and Natural Orbital for Chemical Valence analysis. In order to scrutinize the underlying principles and the origin of the kinetic trans effect, plausible structural constraints were introduced in the decomposition analysis, which allowed eliminating the distance dependence of the interaction energy components. It was established that the trans effect can be rationalized with the interaction of the TPtCl2 and NH3 fragments in the reactant state and TPtCl2 and (NH3)2 fragments in the transition state. It was evinced quantitatively that the σ-donor ability of T indeed controls the stability of the reactant, whereas in the case of π-acids, backdonation stabilizes the transition state, for which conceptually two mechanisms are available: intrinsic and induced π-backdonation. In the destabilization of the reactant and also in the labilization of the leaving group (trans influence) repulsion plays a more important role than orbital sharing effects, which are the cornerstones of the widely accepted interpretations of the trans influence, such as competition for donation or limitation of the donation of the leaving group by the trans ligand T. This repulsive interaction was rationalized both in terms of donated electron density and also in the molecular orbital framework. NOCV orbitals indeed clearly show that the σ-trans effect can be envisioned as a donation from the trans ligand not only to the metal but also to the σ* orbital of the metal-leaving group bond, which manifests as a repulsion between the metal and the leaving group.

The reactivity of dicoordinated Sn(0) compounds, stannylones, is probed using density functional theory (DFT)-based reactivity indices and compared with the reactivity of dicoordinated Sn(II) compounds, stannylenes. For the former compounds, the influence of different types of electron-donating ligands, such as cyclic and acyclic carbenes, stannylenes and phosphines, on the reactivity of the central Sn atom is analyzed in detail. Sn(0) compounds are found to be relatively soft systems with a high nucleophilicity, and the plots of the Fukui function f(-) for an electrophilic attack consistently predict the highest reactivity on the Sn atom. Next, complexes of dicoordinated Sn compounds with different Lewis acids of variable hardness are computed. In a first part, the double-base character of stannylones is demonstrated in interactions with the hardest Lewis acid H(+). Both the first and second proton affinities (PAs) are high and are well correlated with the atomic charge on the Sn atom, probing its local hardness. These observations are also in line with electrostatic potential plots that demonstrate that the tin atom in Sn(0) compounds bears a higher negative charge in comparison to Sn(II) compounds. Stannylones and stannylenes can be distinguished from each other by the partial charges at Sn and by various reactivity indices. It also becomes clear that there is a smooth transition between the two classes of compounds. We furthermore demonstrate both from DFT-based reactivity indices and from energy decomposition analysis, combined with natural orbitals for chemical valence (EDA-NOCV), that the monocomplexed stannylones are still nucleophilic and as reactive towards a second Lewis acid as towards the first one. The dominating interaction is a strong σ-type interaction from the Sn atom towards the Lewis acid. The interaction energy is higher for complexes with the cation Ag(+) than with the non-charged electrophiles BH(3), BF(3), and AlCl(3)., (Copyright © 2013 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim.)

Low energy electron-attachment-induced damage in DNA, where dissociation channels may involve multiple bonds including complex bond rearrangements and significant nuclear motions, is analyzed here. Quantum mechanics/molecular mechanics (QM/MM) calculations reveal how rearrangements of electron density after vertical electron attachment modulate the position and dynamics of the atomic nuclei in DNA. The nuclear motions involve the elongation of the P-O (P-O(3') and P-O(5')) and C-C (C(3')-C(4') and C(4')-C(5')) bonds for which the acquired kinetic energy becomes high enough so that the neighboring C(3')-O(3') or C(5')-O(5') phosphodiester bond may break almost immediately. Such dynamic behavior should happen on a very short time scale, within 15-30 fs, which is of the same order of magnitude as the time scale predicted for the excess electron to localize around the nucleobases. This result indicates that the C-O phosphodiester bonds can break before electron transfer from the backbone to the base.

The surprising difference in the cationic ring-opening polymerization rate of 2-cyclopropyl-2-oxazoline versus 2- n -propyl-2-oxazoline and 2-isopropyl-2-oxazoline was investigated both experimentally and theoretically. The polymerization kinetics of all three oxazolines were experimentally measured in acetonitrile at 140 °C, and the polymerization rate constant ( k p ) was found to decrease in the order c -PropOx > n -PropOx > i -PropOx. Theoretical free energy calculations confirmed the trend for k p , and a set of DFT-based reactivity descriptors, electrostatics, and frontier molecular orbitals were studied to detect the factors controlling this peculiar behavior. Our results show that the observed reactivity is dictated by electrostatic effects. More in particular, the charge on the nitrogen atom of the monomer, used to measure its nucleophilicity, was the most negative for c -PropOx. Furthermore, the electrophilicity of the cations does not change substantially, and thus, the nucleophilicity of the monomers is the driving factor for k p .

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.

meso-Substituted pentaphyrins(1.1.1.1.1) were unexpectedly isolated as N-fused species under Rothemund-type conditions. The reaction mechanism is unknown at present, but the first example of a nonfused [22]pentaphyrin was reported in 2012. Here, the conformational preferences and N-fusion reaction of [22]- and [24]pentaphyrins have been investigated using density functional calculations, together with their aromaticity-molecular topology relationships. Two global minima are found for the unsubstituted [22]pentaphyrin corresponding to T0 and T0(4,D) Hückel structures. Möbius transition states are located in the interconversion pathways with activation barriers of 27 kcal mol(-1). Conversely, [24]pentaphyrin is able to switch between Hückel and Möbius conformers with very low activation barriers. However, nonfused [24]pentaphyrins are unstable and spontaneously undergo an N-fusion reaction driven by the strain release. On the contrary, nonfused [22]pentaphyrins could be isolated if a T0(4,D) conformation is adopted. Importantly, conformational control of pentaphyrins can be achieved by meso-substituents. Two stable conformations (T0(4,D) and T0(A,D)) are found for the nonfused [22]pentaphyrin, which are delicately balanced by the number of substituents. The T0(A,D) conformation is preferred by fully meso-aryl pentaphyrins, which is converted to the N-fused species. Interestingly, the removal of one aryl group prevents the N-fusion reaction, providing stable aromatic nonfused [22]pentaphyrins in excellent agreement with the experimental results.

The linear response kernel is used to gain insight into the aromatic behavior of the less classical metal aromatic E4(2-) and CE4(2-) (E = Al, Ga) clusters. The effect of the systematic replacement of the aluminum atoms in Al4(2-) and CAl4(2-) by germanium atoms is studied using, Al3Ge-, Al2Ge2, AlGe3+, Ge4(2+), CAl3Ge-, CAl2Ge2, CAlGe3+, and CGe4(2+). The results are compared with the values of the delocalization index (δ(1,3)) and nucleus independent chemical shifts (NICS(zz)). Unintegrated plots of the linear response, computed for the first time on molecules, are used to analyze the delocalization in these clusters. All aromaticity indices studied, the linear response, δ(1,3), and NICS(zz), predict that the systems with a central carbon are less aromatic than the systems without a central carbon atom. Also, the linear response is more pronounced in the σ-electron density than in the π-density, pointing out that the systems are mainly σ-aromatic.

The effect of non-polar and polar ligands and of monovalent cations on the one-electron reduction potential of the thiyl radical and the disulfide bond was evaluated. The reduction potentials E° for the CH3S(.)-nL/CH3S(-)-nL and CH3SSCH3-L/CH3SSCH3(.-)-L redox couples were calculated at the B3LYP, M06-2X and MP2 levels of theory, with n=1, 2 and L=CH4, C2H4, H2O, CH3OH, NH3, CH3COOH, CH3CONH2, NH4(+), Na(+), K(+) and Li(+). Non-polar ligands decrease the E° value of the thiyl radical and disulfide bond, while neutral polar ligands favour electron uptake. Charged polar ligands and cations favour electron capture by the thiyl radical while disfavouring electron uptake by the disulfide bond. Thus, the same type of ligand can have a different effect on E° depending on the redox couple. Therefore, properties of an isolated ligand cannot uniquely determine E°. The ligand effects on E° are discussed in terms of the vertical electron affinity and reorganization energy, as well as molecular orbital theory. For a given redox couple, the ligand type influences the nature of the anion formed upon electron capture and the corresponding reorganization process towards the reduced geometry., (Copyright © 2013 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.)

2-Picolyl azide reacts with cis-[PtCl2(DMSO)2] to form the diimino complex [Pt(II)Cl2{NH═C(H)Py}] with subsequent dinitrogen liberation. The formation of the latter complex is scrutinized in a combined experimental and theoretical analysis. We establish in silico that the transformation involves a highly reactive intermediate containing a Pt═N double bond formed after the extrusion of N2 from the azide functionality. The prerequisites for N2 liberation and for the stabilization of the nitrene-related intermediate are analyzed in detail.

In this study, for the first time inverse design was applied to search for the intrinsically most stable radical system in a predefined chemical space of enormous size by scanning in a rational way that entire chemical space. The focus was predominantly on thermodynamic stabilization effects, such as stabilization through resonance. Two different properties were optimized: a newly introduced descriptor called the radical delocalization value and the intrinsic stability via a previously established bond dissociation enthalpy model. The thiadiazinyl radical was chosen as case study of this new approach of inverse design in stable radical chemistry. The resulting optimal structure is found to be highly stable, intrinsically more so than other well-known stable radicals, such as verdazyls and N,N-diphenyl-N'-picrylhydrazyl, and even rivaling the intrinsic stability of nitrogen monoxide.