J. Mol. Struct. 1137, 186-192 (2017)

DOI: 10.1016/j.molstruc.2017.02.032

The structure of N-arylindazoles and their aza-derivatives in the solid state: A systematic analysis of the Cambridge Structural Database coupled with DFT calculations

A search in the Cambridge Structural Database for N-aryl indazoles and their aza derivatives affords 227 structures (183 1-aryl and 44 2-aryl). To discuss their structures, DFT calculations on 20 model compounds were carried out. The geometry of the five-membered ring (the pyrazole) and the conformation if the N-aryl substituent were analyzed.

Struct. Chem. 28, 345-355 (2017)

DOI: 10.1007/s11224-016-0882-y

An insight on the aromatic changes in closed shell icosagen, tetrel, and pnictogen phenalenyl derivatives

A computational study of the aromatic and antiaromatic characteristics of closed shell charged phenalenyl (PLY+1 and PLY−1) upon replacement of the central carbon atom by icosagen (B, Al and Ga), tetrel (Si and Ge) and pnictogen (N, P and As) atoms comprising systems in which the icosagen and pnictogen derivatives considered are neutral while the tetrel ones are anions or cations, has been carried out at the B3LYP/6–311++G(d,p) computational level. By substitution, two different kinds of structures have been obtained, one planar (N and B) and another one bowl-shaped depending on the size of the central atom. In terms of aromaticity, the substitution of the central C atom causes a loss of the aromatic character in all cases as indicated by nucleus-independent chemical shifts (NICS) profiles and NICS values on the 0.001 au isosurface. Regarding the charge, PLY+1 presents larger electron delocalisation than PLY−1, phenomenon associated with aromaticity. Furthermore, the current density maps for those planar systems corroborate NICS findings, showing anticlockwise currents in PLY+1 (like in benzene) but clockwise in PLY-N0 and PLY-B0, indicating aromatic and antiaromatic behaviour, respectively.

J. Org. Chem., 81, 7612–7625 (2016)

DOI: 10.1021/acs.joc.6b01330

15N NMR Spectroscopy, X-ray and Neutron Diffraction, Quantum-Chemical Calculations, and UV/vis-Spectrophotometric Titrations as Complementary Techniques for the Analysis of Pyridine-Supported Bicyclic Guanidine Superbases

Pyridine substituted with one and two bicyclic guanidine groups has been studied as a potential source of superbases. 2-{hpp}C5H4N (I) and 2,6-{hpp}2C5H3N (II) (hppH = 1,3,4,6,7,8-hexahydro-2H-pyrimido[1,2-a]pyrimidine) were protonated using [HNEt3][BPh4] to afford [I-H][BPh4] (1a), [II-H][BPh4] (2), and [II-H2][BPh4]2 (3). Solution-state 1H and 15N NMR spectroscopy shows a symmetrical cation in 2, indicating a facile proton-exchange process in solution. Solid-state 15N NMR data differentiates between the two groups, indicating a mixed guanidine/guanidinium. X-ray diffraction data are consistent with protonation at the imine nitrogen, confirmed for 1a by single-crystal neutron diffraction. The crystal structure of 1a shows association of two [I-H]+ cations within a cage of [BPh4]− anions. Computational analysis performed in the gas phase and in MeCN solution shows that the free energy barrier to transfer a proton between imino centers in [II-H]+ is 1 order of magnitude lower in MeCN than in the gas phase. The results provide evidence that linking hpp groups with the pyridyl group stabilizes the protonation center, thereby increasing the intrinsic basicity in the gas phase, while the bulk prevents efficient cation solvation, resulting in diminished pKa(MeCN) values. Spectrophotometrically measured pKa values are in excellent agreement with calculated values and confirm that I and II are superbases in solution.

Chem. Commun., 2016, 52, 9656-9659.

DOI: 10.1039/c6cc04350j

On the existence of intramolecular one-electron Be–Be bonds

Although the Be–Be bond is extremely weak in Be2 dimers, we have shown that rather stable Be–Be one-electron sigma bonds are formed upon electron attachment to 1,8-diBeX-naphthalene derivatives. Wavefunction analyses corroborate the formation of Be–Be covalent linkages in which the extra electron is accommodated between the Be atoms as reflected in the dramatic shortening of the Be–Be distance with respect to the corresponding neutral molecule.

Graphical abstract: On the existence of intramolecular one-electron Be–Be bonds

Bull. Chem. Soc. Japan, 89, 762-769 (2016)

DOI: 10.1246/bcsj.20160026

The Thermodynamic Stability of Adamantylideneadamantane and Its Proton- and Electron-Exchanges. Comparison with Simple Alkenes

We report herein the results of an experimental and computational study of adamantylideneadamantane (1) and a variety of substituted ethylenic hydrocarbons. The standard enthalpy of formation in the gas phase as well as the gas-phase basicity (GA) of 1were experimentally determined for the first time, respectively by calorimetric techniques and FT-ICR spectrometry. In parallel, computational studies at the MP2/6-311+G(d,p), G3(MP2), and G3 levels were performed on the neutral (1) and protonated (1H+). The agreement with experimental results was very good. The structures of 1 and 1H+ were subject to treatment by “Atoms in Molecules” in order to assess the characteristics of the closest H···H interactions involving both adamantane moieties. Also, the second-order perturbation analysis within the Natural Bond Orbital Theory methodology shows four degenerate charge-transfer interactions between the σ C–H bond of one of the adamantyl subunits towards the σ* C–H of the other adamantyl subunit. The standard enthalpies of formation of new adamantyl compounds were obtained using our experimental data. The computational study of a variety of ethylenic compounds including cyclohexylidenecyclohexane and several alkyl-substituted ethylenes using isodesmic and homodesmotic reactions was carried out. This study was extended to their proton affinities and gas-phase basicities.

Chem. Eur. J. 22, 10651–10660 (2016)

DOI: 10.1002/chem.201600932

Saddle-Shaped Cyclic Indole Tetramers: 3D Electroactive Molecules

We present a joint theoretical and experimental study of a series of cyclic indole tetramers aimed at understanding the fundamental electronic properties of this 3D platform and evaluating its potential in the construction of new semiconductors. To this end, we combined absorption and Raman spectroscopy, cyclic voltammetry, and spectroelectrochemistry with DFT calculations. Our results suggest that this platform can be easily and reversibly oxidized. Additionally, it has a HOMO that matches very well with the workfunction of gold, therefore charge injection from a gold electrode is expected to occur without significant barriers. Interestingly, the cyclic tetraindoles allow for good electron delocalization in spite of their saddle-shaped structures. The steric constraints introduced by N-substitution significantly inhibits ring inversion of the central cyclooctatetraene unit, whereas it only barely affects the optical and electrochemical properties (a slightly higher oxidation potential and a blueshifted absorption upon alkylation are observed).

Tetrahedron 72, 4690-4699 (2016)

DOI: 10.1016/j.tet.2016.06.050

Aromatic changes in isoelectronic derivatives of phenalenyl radicals by central carbon replacement

A computational study of the aromatic characteristics of phenalenyl (PLY) upon replacement of the central carbon atom by Si, Ge, N, P, As, B, Al and Ga atoms comprising isoelectronic series with different charges (0, +1, −1) has been performed at the B3LYP/6-311++G(d,p) level. Two different geometries have been obtained, one planar and another one bowl-shaped. PLY derivatives exhibit a loss of the aromatic character in all cases indicated by NICS (Nucleus-Independent Chemical Shifts) values. These aromatic features are also in agreement with the pyramidalisation undergone by the central atom. The current density maps for those planar systems corroborate NICS findings. These variations in the aromaticity have been rationalised in terms of charge localisation. Harmonic oscillator model of aromaticity (HOMA) and para-delocalisation indices (PDI) have been used to evaluate the electron density delocalisation. HOMA values do not correlate with NICS but PDI indices follow the same trend than NICS.

 Image for unlabelled figure

Chem. Eur. J. 22, 9804-9811 (2016)

DOI: 10.1002/chem.201601201

The Conformational Map of Volatile Anesthetics: Enflurane Revisited

Previous ambiguities in the conformational and structural landscape of the volatile anesthetic enflurane have been solved combining microwave spectroscopy in a jet expansion and ab initio calculations. The broadband (2–18 GHz) rotational spectra identified three different rotamers, sharing a common trans ether skeleton but differing in the ±gauche/trans position of the terminal chlorine atom. For each chlorine conformation two different gauche orientations were predicted for the opposite difluoromethyl group, but only one is experimentally observable due to collisional relaxation in the jet. The experimental dataset comprised nine different isotopologues (35Cl, 37Cl, 13C) and a large number (>6500) of rotational transitions. The inertial data provided structural information using the substitution and effective procedures. The structural preferences were rationalized with additional ab initio, natural-bond-orbital and non-covalent-interaction analysis, which suggest that plausible anomeric effects at the difluoromethyl group could be overridden by other intramolecular effects. The difluoromethyl orientation thus reflects a minimization of inter-fluorine repulsions while maximizing F⋅⋅⋅H attractive interactions. A comparison with previous electron diffraction and spectroscopic data in the gas and condensed phases finally resulted in a comprehensive description of this ether, completing a rotational description of the most common multi-halogenated anesthetics.

Tet. Lett. 57, 1838-1842 (2016)

DOI: 10.1016/j.tetlet.2016.03.044

Theoretical studies on the conformation of large carbocyclic rings. I. 5,6,11,12,17,18-Hexahydrotribenzo[a,e,i]cyclododecane (1,2;5,6;9,10-tribenzododeca-1,5,8-triene)

Tribenzododecatriene is a simple and symmetric molecule that constitutes an excellent model for comparing structural properties both measured and computed. B3LYP/6-311++G(d,p) calculations perfectly reproduced the geometry and NMR properties of this compound in its C2 geometry.

Image for unlabelled figure

Struc. Chem. 27, 705-730 (2016)

DOI: 10.1007/s11224-015-0704-7

The structure of β-diketones related to curcumin determined by X-ray crystallography, NMR (solution and solid state) and theoretical calculations

Structural data are reported on sixteen ketoenols of β-diketones: solution NMR, solid-state NMR (CPMAS and MAS) and X-ray crystallography (four compounds, where three are new). The emphasis is on the tautomerism between both ketoenols, in solution and in the solid state. GIAO/B3LYP/6-311++G(d,p) and Quantum ESPRESSO (QE) calculations were used and compared. For average values, the GIAO/DMSO-PCM is enough, but splittings can only be approached by using QE. A case of rotational disorder has been analyzed. Some anomalies related to C–F bonds and to the C–CF3 group have been detected.

Struct. Chem. 27, 91-99 (2016)

Download from the journal web page, DOI: 10.1007/s11224-015-0721-6 

Aromaticity of benzenoid hydrocarbons with inserted –B=B– and –BH–BH– groups: a comparison

Structures of selected polycyclic conjugated hydrocarbons with –B=B– and –BH–BH– moieties inserted in different places were calculated at the B3LYP/6-311++G** level and their aromatic properties evaluated. HOMA, NICS(0), NICS(1)zz, Λ and PDI indices were used for studying their aromatic properties. Both optimized planar (as in parent hydrocarbons) and non-planar structures were taken into account. It is shown that insertion of both types of boron groups disturbs and decreases the aromaticity of the corresponding hydrocarbons. The decreasing effect of the –BH–BH– group is much stronger. What is quite intriguing is that it appears that non-planar structures of the studied compounds have a little higher aromaticity than the strictly planar ones. Mutual correlations between results obtained by different aromaticity indices are calculated and thoroughly discussed.

Phys. Chem. Chem. Phys., 17, 14961-14971 (2015)

DOI: 10.1039/C5CP00876J

Influence of fluoro and cyano substituents in the aromatic and antiaromatic characteristics of cyclooctatetraene

An exhaustive and systematic study of the structural and electronic properties of cyclooctatetraene (COT) upon substitution of hydrogen atoms by fluoro and cyano groups has been carried out in order to analyse the influence of both substituents on the aromaticity. We found that C–C distances decrease with fluoro substitution while in cyano derivatives the opposite happens. All the compounds retain their original structural type, with the exception of the cyano derivatives; thus, compounds 25CN6T, 27CN6T and 30CN8T show boat-like structure, whereas compounds 20CN5T,26CN6T, and 29CN7T present twisted structures. Regarding the relative energies of those compounds with the same number of substitutions, it was found that compounds where the X groups were more separated among them were the most stable ones. Inversion barriers (ΔETS) were found to increase with the number of substitutions; in the case of fluoro derivatives these barriers have a two-fold, increase compared to the parent compound while in the cyano ones a three-fold increase was observed. The aromatic character based on the NICS values, was found to increase in the ground singlet states and in the transition states of both fluoro and cyano derivatives. For triplet states, a decrease of the aromatic behaviour was found upon substitution. NICS profiles and 3D NICS isosurfaces confirm such findings. Finally, HOMA indexes corroborate the aromatic changes described by the NICS values, although, no good correlations between both quantities were found.

Graphical abstract: Influence of fluoro and cyano substituents in the aromatic and antiaromatic characteristics of cyclooctatetraene

ChemPhysChem, 16, 465-469 (2015)

DOI: 10.1002/cphc.201402711

Linear Free-Energy Relationships between a Single Gas-Phase Ab Initio Equilibrium Bond Length and Experimental pKa Values in Aqueous Solution

Remarkably simple yet effective linear free energy relationships were discovered between a single ab initio computed bond length in the gas phase and experimental pKa values in aqueous solution. The formation of these relationships is driven by chemical features such as functional groups, meta/para substitution and tautomerism. The high structural content of the ab initio bond length makes a given data set essentially divide itself into high correlation subsets (HCSs). Surprisingly, all molecules in a given high correlation subset share the same conformation in the gas phase. Here we show that accurate pKa values can be predicted from such HCSs. This is achieved within an accuracy of 0.2 pKa units for 5 drug molecules.

ChemPhysChem 15, 3493-3501 (2014)

DOI: 10.1002/cphc.201402442

A Computational Study on 3-Azonia-, 3-Phosphonia-, and 3-Arsoniaspiro[2.2]pentanes and Related Three-Membered Heterocycles

A theoretical study at the ab initio MP2/6-311++G(d,p) level of theory is carried out to characterize several heterocyclic spiro[2.2]pentane cations with N, P, and As as spiro atoms. The strain and relative stability of the spiropentanes are obtained through isodesmic reactions. Nucleus-independent chemical shifts (NICS) and 3D NICS isosurfaces show σ-aromatic characteristics, similar to those found in cyclopropane. The interaction with the Cl anion, which results in four different stationary structures, is studied and characterized by means of the atoms in molecules methodology, and Cl⋅⋅⋅pnicogen, Cl⋅⋅⋅H, and Cl⋅⋅⋅C interactions are found. The most stable structure in all cases corresponds to opening of one of the three-membered rings, due to the attack of the Cl atom, and C[BOND]Cl bond formation. Furthermore, the reaction with the 3-boranuidaspiro[2.2]pentane anion results in the formation of a new compound through cleavage of one ring of both reactants.

Eur. J. Org. Chem. 2014, 5370-5377 (2014)

DOI: 10.1002/ejoc.201402509

Aromaticity and Magnetic Properties of 1- and 2-Indenones and Their Aza Derivatives

The almost unknown series of 1- and 2-indenones and their aza-derivatives have been analyzed with the aim of understanding the higher stability and aromaticity of the 1-one compared to the 2-one series. The effect of tetrafluorination of the six-membered ring on their relative stability and aromaticity and the corresponding 19F–19F spin–spin coupling constants (SSCC) have been analyzed for the compounds of these two series. Magnetic NICS and electronic PDI, FLU, and MCI aromaticity criteria, complemented by an energy decomposition analysis, allow the higher stability and aromaticity of 1-indenones compared with 2-indenones to be attributed to the π-component of the orbital interaction term. The reduced aromatic character of the 2-one series also provides an explanation for the 19F–19F SSCC observed in the tetrafluorinated 2-one compounds, which differ significantly from those of 1,2,3,4-tetrafluorobenzene.

Helv. Chim. Acta 97, 471-490 (2014)

DOI: 10.1002/hlca.201300395

A Spectroscopic Study of Colchicine in the Solid State and in Solution by Multinuclear Magnetic Resonance and Vibrational Circular Dichroism

Although almost 200-years-old, several unknown aspects remain to be explored of colchicine, the unique available drug for acute flares of gout. In this article, we report density-functional theory (DFT) studies of geometry, energy, and NMR; 1H-, 13C-, and 15N-NMR chemical shifts and some spin-spin coupling constants, including the complete analysis of the saturated part (ring B); the assignment of both enantiomers by NMR using a chiral solvating agent; solid-state NMR experiments of the different forms of natural and racemic colchicine, and IR and vibrational circular dichroism (VCD) studies of these same forms.

Organometallics, 33, 616–619 (2014)

DOI: 10.1021/om401215m

Arene–Ruthenium Chemistry and Brønsted Acid Catalysis of a Chiral Phosphane-Hydroxyl Ligand

The phosphane-hydroxyl ligand (SC1,RC2)-Ph2PC(Ph)HC(OH)HCH2OMe (POH) displays κ1P, κ2P,O, κ3P,O,O′, and anionic deprotonated κ3P,O,O′ coordination modes toward the metallic (η6-p-MeC6H4iPr)Ru moiety. The hydroxyl group in [(η6-p-MeC6H4iPr)Ru(κ3P,O,O′-POH)][SbF6]2(3) has strong acidic properties. This complex catalyzes the Diels–Alder reaction between cyclopentadiene and trans-β-nitrostyrene or methacrolein as well as the Friedel–Crafts reaction between indole and methyl 3,3,3,-trifluoropyruvate or trans-β-nitrostyrene. Solution studies strongly indicate that complex 3 acts as a metallic Lewis acid assisted Brønsted-acid catalyst.

ARKIVOC, ii, 319-332 (2014)

Download from the journal page

Preparation of 2,2-dimethylchroman-4-ones from 5-alkyl-substituted resorcinols: microwave-assisted synthesis and theoretical calculations

The influence of different 5-alkyl-substituted resorcinols on the formation of 2,2-dimethyl-chroman-4-ones is examined experimentally and theoretically. Structures are fully assigned by means of experimental and theoretical 13C and 1H NMR chemical shifts. Based on experimental and theoretical calculations of Friedel-Crafts acylation, it is possible to explain the formation of 2,2-dimethyl-5-hydroxychroman-4-ones and/or 2,2-dimethyl-7-hydroxychroman-4-ones. Evaluation of their biological activity as cannabinoid receptor ligands is also reported. 

J. Phys. Org. Chem. 2013, 26 791–796

DOI: 10.1002/poc.3159

Relationship between experimental pKa values in aqueous solution and a gas phase bond length in bicyclo[2.2.2]octane and cubane carboxylic acids

Linear correlations were established between the calculated bond lengths and the pKa or σI values for a series of 4-substituted bicyclo[2.2.2]octane-1-carboxylic acid and 4-cubane-1-carboxylic acid derivatives. The bond lengths have been calculated at a modest computational level, HF/6-31G(d), both in the gas phase and with the continuum solvation model, polarisable continuum model (PCM). In general, the best correlations are obtained when the PCM model is taken into account, especially when neutral and charged molecules are considered together. The best models in each case show square correlation coefficients (R2) larger than 0.9 and indicate that they can be used as predictive tools. These results expand previous results that indicate the possibility of a relationship between gas phase bond length and pKa values in aqueous solution and indicate that such relationships are more general than hitherto expected. 

Cent. Eur. J. Chem. 11, 1711-1722 (2013)

DOI: 10.2478/s11532-013-0311-7

A theoretical study of the limits of the acidity of carbon acids in phase transfer catalysis in water and in liquid ammonia

The acidities of a large number of carbon acids have been theoretically calculated for the gas-phase and for DMSO solution. The gas-phase values, both ΔH and ΔG, are very well correlated with the available experimental data. From the calculated ΔG values in DMSO and the pKas in the same solvent, a homogeneous set of pKa (DMSO) values was devised that was used to generate pKa(water). These last pKas were used to establish the limits of the acidity of carbon acids for reactions under PTC conditions both alkylations and H/D exchange. A step further led to the pK as in liquid ammonia and from them to the virtual use of PTC using liquid ammonia instead of water.

Phys. Chem. Chem. Phys. 15, 14026-14036 (2013)

DOI: 10.1039/C3CP51963E

Modeling the mechanism of glycosylation reactions between ethanol, 1,2-ethanediol and methoxymethanol

The mechanism of the SN2 model glycosylation reaction between ethanol, 1,2-ethanediol and methoxymethanol has been studied theoretically at the B3LYP/6-311+G(d,p) computational level. Three different types of reactions have been explored: (i) the exchange of hydroxyl groups between these model systems; (ii) the basic catalysis reactions by combination of the substrates as glycosyl donors (neutral species) and acceptors (enolate species); and (iii) the effect on the reaction profile of an explicit H2O molecule in the reactions considered in (ii). The reaction force, the electronic chemical potential and the reaction electronic flux have been characterized for the reaction path in each case. Energy calculations show that methoxymethanol is the worst glycosyl donor model among the ones studied here, while 1,2-ethanediol is the best, having the lowest activation barrier of 74.7 kJ mol−1 for the reaction between this one and the ethanolate as the glycosyl acceptor model. In general, the presence of direct interactions between the atoms involved in the penta-coordinated TS increases the activation energies of the processes.

Struct. Chem. 24, 542-548 (2013)


HOMA parameters for the boron–boron bond: How the introduction of a BB bond influences the aromaticity of selected hydrocarbons

An extension of the harmonic oscillator model of aromaticity (HOMA) model for systems with boron–boron bonds is presented. For the first time, the parameters of the HOMA model are estimated using only theoretically calculated bond lengths. The HOMA parameters obtained make geometric aromaticity studies possible for a large number of compounds containing the boron–boron bond. The derived HOMA parameters have been used to investigate how the introduction of the boron–boron moiety in the structure of selected hydrocarbons modifies their aromaticity. The conclusion is that the insertion of a boron–boron bond usually strongly decreases the aromaticity of the boron-containing compounds in comparison to their parent hydrocarbons.

Chem. Eur. J. 18, 13379–13387 (2012)


Regioselectivity in the intramolecular heck reaction of a series of cyclic sulfonamides: an experimental and computational study.

Regioselectivity in the intramolecular Heck reaction of a series of N-sulfonyl-2,5-dihydro-3-substituted pyrroles was studied. These substrates are unbiased in terms of the formed ring size of the new heterocycle. Results indicate that high levels of regioselectivity are observed under a range of conditions, and that there is an underlying propensity for carbon-carbon bond formation at the most hindered end of the alkene. For two examples (3-Me and 3-tBu), DFT calculations were performed and indicate that in both cases, the modelled transition state for carbopalladation is energetically lower for the experimentally preferred isomer.

J. Mol. Struct., 1015,138-146 (2012)


The tautomeric structures of 3(5),3′(5′)-azopyrazole [(E)-1,2-di(1H-pyrazol-3(5)-yl)diazene)]: The combined use of NMR and electronic spectroscopies with DFT calculations

An azo derivative of 1H-pyrazole has been chosen to study the information obtained from different techniques for determining the structure of a non-crystalline compound in the solid-state and in solution. Syn–antiisomerism of the azo group, prototropic tautomerism of the 1H-pyrazole and rotation about the pyrazole–azo group resulted in 20 structures that were analyzed. Energy calculations, 13C and 15N chemical shifts, 1H–1H coupling constants and electronic spectra reduced the 20 possible structures to only one, the 3,3′-Z,Z-anti-azopyrazole.

Struct. Chem. 23, 595-600 (2012)

DOI: 10.1007/s11224-011-9907-8

Calculation of the HOMA model parameters for the carbon-boron bond

An extension of the harmonic oscillator model of aromaticity (HOMA) model to systems with carbon–boron bonds is presented. Model parameters were estimated using experimental and theoretical bond lengths. It is shown that both approaches produce very similar HOMA models. In the second part of the article, the aromaticity levels of several model compounds containing carbon–boron bonds are calculated using the previously obtained parameters. The results of these calculations are compared with those provided by other aromaticity indices. The aromaticity of boron-containing compounds is also compared with the aromaticity of analogous compounds containing carbon and nitrogen.

Struct. Chem. 23, 1245-1252 (2012)

DOI: 10.1007/s11224-012-9991-4

A theoretical study of six-membered rings containing the -N=S-S=N- motif

The sulfur electronic structure allows for the existence of hypervalent molecules that were discussed as early as 1939 by Schomaker and Pauling. In the present paper, we report a theoretical study of three hypothetical possible six-membered rings that contain the –N=S–S=N– motif carried out at the MP2/aug-cc-pVDZ level. Three minima have been found, of which two of them are six-membered rings with all covalent bonds. Since they do not have a planar structure, we have calculated the aromaticity index NICS through the symmetry axis perpendicular to the mid-plane and at various points in space. Benzene, cyclohexane, and 1,2,5-thiadiazole have served as points of comparison.

Phys. Chem. Chem. Phys. 14, 11468-11477 (2012)

DOI: 10.1039/C2CP41230F

Strong interactions between copper halides and unsaturated systems: new metallocycles? Or the importance of deformation

The complexes formed by CuF with CC double and triple bonds have been studied at the MP2 and CCSD(T) computational levels. The interaction of CuF with acetylene, ethylene and their fluoro derivatives is very strong, with interaction energies close to those of conventional covalent bonds. Hence, these complexes could be actually viewed as a new kind of metallocycles, with significantly strong Cu–C linkages. All electronic indexes analyzed by means of the AIM, ELF and NBO formalisms, indicate that the strength of the interaction should increase with the number of fluorine substituents in both series of compounds. Surprisingly, however, although both series of compounds exhibit the same bonding arrangements, they follow opposite stability trends and the expected increase of the interaction energies with the number of fluorine substituents is only observed in the acetylene series. The reason for this unexpected behavior is once more associated with the effects triggered by the distortion of the interacting subunits. Deformation not only has a direct energetic cost but dramatically affects the intrinsic properties of the interacting systems.

J. Phys. Chem. A 116, 8250-8259 (2012)

DOI: 10.1021/jp304495f

Mechanisms of Formation of Hemiacetals: Intrinsic Reactivity Analysis

The reaction mechanism of the hemiacetal formation from formaldehyde and methanol has been studied theoretically at the B3LYP/6-311++G(d,p) level. In addition to the study of the reaction between the isolated reactants, three different kinds of catalysis have been explored. The first one examines the use of assistants, especially bridging water molecules, in the proton transfer process. The second one attempts to increase the local electrophilicity of the carbon atom in formaldehyde with the presence of a Brønsted acid (H+ or H3O+). The last one considers the combined effect of both catalytic strategies. The reaction force, the electronic chemical potential, and the reaction electronic flux have been characterized for the reaction path in each case. In general, it has been found that structural rearrangements represent an important energetic penalty during the activation process. The barriers for the reactions catalyzed by Brønsted acids show a high percentage of electronic reorganization contribution. The catalytic effects for the reactions assisted by water molecules are due to a reduction of the strain in the transition state structures. The reaction that includes both acid catalysis and proton assistance transfer shows the lowest energy barrier (25.0 kJ mol–1).

Tetrahedron 67, 7316-7320 (2011)


Isomerization barriers in bis(4H-thiopyran) and in bothioxanthenes

The inversion and rotation mechanisms for the isomerization of Feringa’s bithioxanthenes existing in two conformations, up/up and up/down, have been calculated at the B3LYP/6-31G(d) and B3LYP/6-311++G(d,p) levels. The inversion mechanism that maintains the double bond nature of the central bond is a classical one but the rotation mechanisms that break the double bond to form a biradical needs to explore the singlet and triplet states. To do this we have removed the four fused phenyl rings of bithioxanthene and calculated at the CASSCF and CASPT2 levels bis(4H-thiopyran) proving that B3LYP calculations yield reasonable results for the rotation barriers.

J. Phys. Org. Chem. 24, 744-750 (2011)

DOI: 10.1002/poc.1813

Thermodynamic and kinetic effects of Lewis acid complexation on a Schiff base present in two tautomeric forms

The effect of complexation with Lewis acids on the tautomeric equilibrium of a derivative of pyridoxal has been explored using density functional theory. Three complexation sites (pyridine nitrogen, aromatic ring, and carbonyl group) have been considered. The tautomeric equilibrium can be modulated with the formation of complexes at the different sites. The changes observed in the geometric, energetic and electronic properties of the tautomers and transition states have been analyzed. Relationships have been found between those parameters, including a relationship between energetic and geometric parameters in agreement with the Hammond postulate.

Struct. Chem. 22, 707-715 (2011)

DOI: 10.1007/s11224-011-9749-4

Modeling the allosteric effect: modification of the tautomerism by intermolecular interactions and extension to molecular wires

B3LYP/6−31G(d) and B3LYP/6−311++G(d,p) calculations were carried out on quinolone and its four azaderivatives (15 tautomers), five wire models (10 tautomers) for proton transfer formed of anthracenol and azaacridines and the corresponding five crown ethers (ten tautomers) in the anthracenol part. The wires are formed by four carbon atoms existing either as cumulenes (pentaene) or polyynes (diyne). On these structures the effect of hydrogen bonds with HF as hydrogen bond donor, protonation on the aza N atoms and coordination with Li+ on the same positions were studied. The resulting energies were analyzed taking into account proximity effects.

Comp. Theor. Chem. 966, 334-339 (2011)

DOI: 10.1016/j.comptc.2011.03.034

A theoretical study of the conformation, basicity and NMR properties of 2,2'-, 3,3'- and 4,4'-bipyridines and their conjugated acids

A theoretical study at the B3LYP/6-311++G(d,p) level has been carried out on 2,2′-, 3,3′- and 4,4′-bipyridines, as well as on their monoprotonated and diprotonated forms. The geometries, torsion angles, and the energies of the minima and transition states have been calculated with good agreement with previous calculations and with most of the scarce available data. The absolute shieldings were calculated on these geometries using the GIAO approximation and transformed into chemical shifts using previously established relationships.

Inorg. Chem. 50, 4212-4220 (2011)

DOI: 10.1021/ic102314r

Unprecedented 1,3-Diaza[3]ferrocenophane Scaffold as Molecular Probe for Anions

The guanidine unit in the guise of 2-aminoimidazole in the new structural motif 2-arylamino-1,3-diaza[3]ferrocenophane 4 acts as a binding site for anions. The electrochemical behavior of this compound has been studied by cyclic voltammetry (CV) and differential pulse voltammetry (DPV) and was found to exhibit a quasi reversible oxidation peak, associated to the Fe(II)/Fe(III) redox couple (Ep = 440 mV), and a non-reversible oxidation wave (Ep = 817 mV), probably associated to the oxidation of the C═N unit present in the guanidine bridge. Recognition of AcO, PhCO2, F, Cl, and Br anions by the free receptor and the less basic anions Br, Cl, and NO3 by its monoprotonated form takes place by unusual redox-ratiometric measurements and spectroscopic (1H NMR and UV−vis) changes.

Chem. Phys. Lett. 515, 210-213 (2011)

DOI: 10.1016/j.cplett.2011.09.037

Multiple 3c-2e bonding of methane with metal cations

Methane is a greenhouse gas but at the same time a possible starting material in chemical industry. Here, we demonstrate by computational methods that methane could form in gas phase stable adducts with alkaline and earth alkaline metal cations. The interaction can be described as multiple 3csingle bond2e bonding between a cation and Csingle bondH bonds. In the alkaline metal adducts the open mode prevails, while for the earth alkaline moieties it is difficult to define, whether there is either close or open mode prevailing.

Comput. Theor. Chem. 974, 37-42 (2011)


A theoretical study of the conformation of 2,2’-bifuran,2,2’-bithiophene, 2,2’-bitellurophene and mixed derivates: chalcogen-chalcogen interactions or dipole-dipole effects.

Using as models 2,2′-biheterocycles (furan, thiophene, selenophene, tellurophene) the possibility of stabilizing effects due to the proximity of two chalcogen atoms (O, S, Se, Te) was explored. The geometries of the ground states were compared satisfactorily with the available experimental results. The torsional barriers through planar and orthogonal transition states were calculated and discussed in function of the nature of the chalcogen atoms. Special attention was devoted to dipole moments that were calculated by a vector summation with good results

Chem. Eur. J. 16, 11897-11905 (2010)

DOI: 10.1002/chem.201001254

New Insights into Factors Influencing B N Bonding in X:BH3-nFn and X:BH3-nCln for X=N2, HCN, LiCN, H2CNH, NF3, NH3 and n=0-3: The Importance of Deformation

Understanding the bonding in complexes X:BH3−nFn and X:BH3−nCln, for X=N2, HCN, LiCN, H2CNH, NF3, NH3 with n=0–3, is a challenging task. The trends in calculated binding energies cannot be explained in terms of any of the usual indexes, including π donation from the halogen lone pairs to the p(π) empty orbital on B, deformation energies, charge capacities, or LUMO energies, which are normally invoked to explain the higher Lewis acidity of BCl3 relative to BF3. The results of the high-level G3B3 ab initio calculations reported in this study suggest that the interaction energies of these complexes are determined by a combination of at least three factors. These include the decrease in the electron-accepting ability of B as a result of π donation by the halogen atom, the increase in the electron-acceptor capacity of B due to deformation of the acid, and the large increase in the deformation energy of the acid with increasing halogen substitution. The dominant effects are those derived from the electronic effects of acid deformation. Deformation not only has direct energetic consequences, which are reflected in the large differences between dissociation (D0) and interaction (Eint) energies, but also leads to an enhancement of the intrinsic acidities of BH3−nFn and BH3−nCln moieties by lowering the LUMO energies to very different extents, consistent with the frontier orbital model of chemical reactivity. Although this lowering depends on both the number and the nature of the halogen substituents, binding energies do not systematically increase or decrease as the number of halogen atoms increases.

J. Am. Chem. Soc. 132, 7862-7863 (2010)

DOI: 10.1021/ja1025062

Crystallization of an Achiral Cyclohexanone Ethylene Ketal in Enantiomorphs and Determination of the Absolute Structure

The achiral 4-methoxy-4-(p-methoxyphenyl)-cyclohexanone ethylene ketal (1) resolves spontaneously. The crystal structure is solved in chiral spatial group P2(1). Because compound 1 is composed of only light atoms (C, H, O) it is not possible to determine its absolute structure configuration. An efficient procedure for the absolute structure configuration determination of flexible molecules containing only light atoms is proposed, based on the combination of X-ray diffraction, solid-state VCD, and DFT calculations.

J. Phys. Org. Chem. 23, 526-535 (2010)

DOI: 10.1002/poc.1636

Double proton transfer in crystals of 1,3,4,6,7,8-hexahydro-2H-pyrimido[1,2-a] pyrimidine (hppH): 13C and 15N CPMAS NMR study of (hppH)2

In this paper, we report an example of intermolecular solid-state proton transfer in the bicyclic guanidine, hppH. A combination of X-ray crystallography, CPMAS NMR (13C and 15N) and theoretical calculations allows us to determine that a double proton transfer takes place in the (hppH)2 dimer with an activation energy of about 50 kJ mol−1. According to the B3LYP/6-311++G(d,p) calculations, the double proton transfer occurs non-symmetrically through a zwitterion.

Phys. Chem. Chem. Phys. 12, 5101-5108 (2010)

DOI: 10.1039/b924322d

On the electronic structure and stability of icosahedral r-X2Z10H12 and Z12H122- clusters; r = {ortho, meta, para}, X = {C, Si}, Z = {B, Al}

We report on the electronic structure of the 12-vertexicosahedralclustersr-X2Z10H12and Z12H122–, where X = {C,Si} and Z = {B,Al}. The least stablecluster—with the lowest HOMO–LUMO gap (Eg)—corresponds to theortho-X2Z10H12isomers for all values of X = {C,Si} and Z = {B,Al}. The well-known energetic orderE(para)<E(meta)<E(ortho)forr-carboranes is also valid for all compounds exceptr-C2Al10H12. Substitution of two atoms ofcarbonorsiliconinto theicosahedralcage B12H122−enhances considerably the stability of the system as analyzed fromEggaps, as opposite to Al12H122−, where similar gaps are found upon doublecarbonorsiliconsubstitution regardless of the positions in the cage. In order to highlight similarities and differences in the titleclusters, topological analysis of the electron density was performed, together with analysis of the deviation from polyhedronicosahedralform with (i) volumes, skewness and kurtosis calculations; and (ii) continuous shape measures.

J. Mol. Stuct. 969, 106-110 (2010)

DOI: 10.1016/j.molstruc.2010.01.049

The solid-state structure of primary fatty amines: True amines or ammonium amides?

A computational and experimental NMR study of octadecylamine and hexadecylamine proves that these compounds behave in the solid-state like all other primary amines and not as ammonium amide inner salts.

J. Chem. Inf. Model. 50, 87-96 (2010)

DOI: 10.1021/ci900396k

Prediction of the Basicities of Pyridines in the Gas Phase and in Aqueous Solution

The basicities of 125 pyridine derivatives in the gas phase and in water have been correlated with the electron density properties within the framework of quantum topological molecular similarity (QTMS). We used the theory of quantum chemical topology (QCT) to provide ab initio descriptors that are able to predict pKb values. Partial least squares (PLS) and the machine-learning technique Kriging generated validated models. Properties were considered for systems in their neutral and protonated forms. The compounds were divided into a training set, used to develop the models, and a test set, for which the predicted values of the different models were compared with the experimental ones. The results were found to be good for those compounds with substituents in the meta and para positions, whereas the use of Kriging was required to obtain reasonable results when ortho derivatives were included. The basicity was found to be better described in the gas phase than in water. Special attention was paid to external validation.

J. Chem. Theory Comput. 5, 2763-2771 (2009)

DOI: 10.1021/ct900364y

Beryllium Bonds, Do They Exist?

The complexes between BeX2 (X = H, F, Cl, OH) with different Lewis bases have been investigated through the use of B3LYP, MP2, and CCSD(T) approaches. This theoretical survey showed that these complexes are stabilized through the interaction between the Be atom and the basic center of the base, which are characterized by electron densities at the corresponding bond critical points larger than those found in conventional hydrogen bonds (HBs). Actually, all bonding indices indicate that, although these interactions that we named “beryllium bonds” are in general significantly stronger than HBs, they share many common features. Both interactions have a dominant electrostatic character but also some covalent contributions associated with a non-negligible electron transfer between the interacting subunits. This electron transfer, which in HBs takes place from the HB acceptor lone-pairs toward the σYH* antibonding orbital of the HB donor, in beryllium bonds goes from the lone pairs of the Lewis base toward the empty p orbital of Be and the σBeX* antibonding orbital. Accordingly, a significant distortion of the BeX2 subunit, which in the complex becomes nonlinear, takes place. Concomitantly, a significant red-shifting of the X−Be−X antisymmetric stretching frequencies and a significant lengthening of the X−Be bonds occur. The presence of the beryllium bond results in a significant blue-shifting of the X−Be−X symmetric stretch.

Croat. Chem. Acta 82, 173-183 (2009)


Barriers about Double Carbon-Nitrogen Bond in Imine Derivatives (Aldimines, Oximes, Hydrazones, Azines)

The paper presents the results referring to the inversion mechanism of imines and their derivatives (hydrazones, oximes, azines). The calculated barriers [B3LYP/6-311++G(d,p) and G3B3)] are in good agreement with the scarce existing data. The transition states correspond in all cases to a pure nitrogen inversion except in the case of azines where they have some rotation character. The electron properties of the minima and the transition states have been characterized, allowing explanation of the geometrical changes observed in the process.

Theochem, 896, 92-95 (2009)

DOI: 10.1016/j.theochem.2008.11.005

Does 7-norcaranone (bicyclo[4.1.0]heptan-7-one) exist?

Theoretical calculations have been carried on four compounds, bicyclo[4.1.0]heptan-7-one cis and trans, 2-cyclohepten-1-one and cycloheptenyl cation to determine their UV, IR and NMR spectra (chemical shifts). In the case of the known 2-cyclohepten-1-one the agreement is satisfactory. Predictive properties have been calculated for the elusive bicyclo[4.1.0]heptan-7-one cis molecule.


Tetrahedron 64, 8089-8094 (2008)

Doi: 10.1016/j.tet.2008.06.065

A study of the tautomerism of β-dicarbonyl compounds with special emphasis on curcuminoids

Six β-diketones related to curcumin and curcumin itself have been studied by 13C NMR spectroscopy in chloroform in order to determine the equilibrium constant between the two keto/enol tautomers. In order to do this GIAO/B3LYP/6-31G∗∗ calculations of absolute shieldings (σ, ppm) were carried out. To establish relationships between σ and experimental chemical shifts (δ, ppm), three simple β-diketones (acetylacetone, dibenzoylmethane and benzoylacetone) have been studied. The preference for different groups to be conjugated with the Cdouble bond; length as m-dashO has been determined.

J. Mol. Struct. 891, 325-328 (2008)

Doi: 10.1016/j.molstruc.2008.04.002

Desmotropy in reduced plumbagins: alpha- and beta-Dihydroplumbagins

A report of 1933 about the desmotropy of dihydroplumbagin was examined using both DFT calculations and high field NMR experiments. It was concluded that the report was essentially correct because both tautomers were isolated and characterized.

Struct. Chem. 19, 339-352 (2008)

Doi: 10.1007/s11224-008-9290-2

Substitution effects on neutral and protonated pyridine derivatives along the periodic table

A theoretical study of the monosubstitution effects of all atoms of the second and third rows of the periodic table on the α, β and γ positions of neutral and protonated pyridine has been carried out by means of B3LYP/6-31 + G(d,p) DFT calculations. The geometric and electronic properties, calculated using the Atoms in Molecules methodology, and the electrostatic potential have been analysed. Concurrently, three separate aromaticity indexes (NICS(0), NICS(1) and HOMA) have been evaluated and compared to the above results. Furthermore, the effect of protonation on these parameters has been investigated. A comparison with analogous results for benzene derivatives has also been carried out.

Russ. J. Gen. Chem. 78, 784-792 (2008)

Doi: 10.1134/S1070363208040385

Acid Properties of Porphyrins and Related Systems

Structural, energetic, and electronic aspects of a series of porphyrins and related compounds in their neutral, anion, and dianion electronic states have been studied by means of DFT calculations [B3LYP/6-31+G(d,p)]. Analysis of the electron density puts forward a number of different kinds of interactions between the inner atoms. Intramolecular correlations similar to those found in intermolecular N-H···N interactions have been found.

J. Phys. Org. Chem. 21, 381-386 (2008)

Doi: 10.1002/poc.1350

Theoretical study of racemization in chiral alkenylidene truxenes

DFT based methods have been used to study the racemization process of chiral truxene derivatives. Four minima and the transition state connecting them have been characterized. It has been found that the energetic results for the transition states are correlated for each compound and also among the different systems considered. Thus, a unique equation has been found to fit all the values obtained. Model systems have shown that the calculated barriers are not inherent to the truxene structure. The AIM analysis has shown the presence of a large number of intramolecular closed-shell interactions while electron density and Laplacian at the bond critical points are exponentially correlated with the interatomic distance. 

Arkivoc vii, 48-56 (2008)


Computational studies of the structure of aldazines and ketazines. Part 1. Simple compounds

The molecular structure of seven simple aliphatic and aromatic aldazines and ketazines has been calculated at the B3LYP/6-311++G** computational level. Geometries,  E/Z- isomerism and conformation of the compounds have been compared with the available data (electron diffraction, X-ray crystallography). IR and Raman spectra have been calculated and compared with experimental ones. 

Theochem, 847, 25-31 (2007)

Doi: 10.1016/j.theochem.2007.08.026

Computational studies of the structure of aldazines and ketazines: Part 2. Halogen and alpha,beta-unsaturated derivatives

Nine azines bearing F, Cl, Br, CF3, CF2CF3, CH(CF3)2, CHsingle bondCHdouble bond; length as m-dashCH2, CHsingle bondCHdouble bond; length as m-dashCH(CH3) and CHsingle bondCHdouble bond; length as m-dashCH(C6H5) substituents have been theoretically studied at the B3LYP/6-311++G(d,p) level. Concerning the E/Z isomerism, in all cases, the large groups prefer to be “outside” (in general, E), i.e. for conjugated azines the extended conformation (EEEE) is always preferred. In the case of the CH(CF3)2 substituent, the imine–enamine tautomerism has been studied, the imine one being clearly more stable. The differences in energy between the different isomers and conformers have been statistically analyzed.

Int. J. Mass Spectr., 267, 302-307 (2007)

Doi: 10.1016/j.ijms.2007.02.056

The intrinsic (gas-phase) acidities of bridgehead alcohols. An experimental (FT-ICR) and computational study

The gas-phase acidities of 1-adamantanol and perfluoro1-adamantanol were determined by means of Fourier transform ion cyclotron resonance spectrometry (FT-ICR). The acidity of perfluoro1-adamantanol seems to be the highest ever reported for an alcohol. A computational study of these species and their anions at both the MP2/6-311 + G(d,p) and B3LYP/6-311 + G(d,p) levels was performed. Also studied were the tertiary alcohols (including their perfluorinated forms) derived from norbornane, bicyclo[2.2.2]octane and cubane. It was found that: (i) the intrinsic acidity of non-fluorinated bridgehead alcohols increases with the strain of the hydrocarbon framework and, (ii) perfluorination of these compounds strongly increases their acidity and, likely, significantly modifies their internal strain.

Theor. Chem. Acc. 118, 533-539 (2007)

Doi: 10.1007/s00214-007-0340-4

Theoretical study of neutral and protonated triple bonded molecules formed between C, N, Si, P, B and Al

A theoretical study of the possible protonation sites of simple molecules formed by C, N, Si, P, B and Al that present a triple bond between those atoms has been carried out. The calculations performed include MP2 and CCSD(T) methods with the aug-cc-pVTZ basis set. The nature of the protonated species has been analyzed with the Atoms In Molecules methodology.

Phys. Chem. Chem. Phys. 9, 3970-3977 (2007)

Doi: 10.1039/b702480k

Unusual substituent effects on the bonding of iminoboranes

Substituent effects on iminoboranes XBNH, HBNX and XBNX (X = H, CH3, NH2, OH, F) have been analyzed in the framework of the NBO, AIM and ELF approaches, using B3LYP/6-311++G(d,p) optimized geometries and electron densities. Boron-substituted derivatives, XBNH, are more stable than the corresponding nitrogen-substituted isomers HBNX, with the energy difference increasing as the electron withdrawing character of the substituent increases. The BN linkage is not much affected by N-substitution, but it is significantly altered when the substituent is attached to the boron atom in both XBNH and XBNX series of compounds. Moreover, substituent effects on the structures of iminoboranes are opposite those observed for the corresponding isoelectronic acetylene derivatives. The ELF analysis indicates that electron-withdrawing substituents enhance the localization of electrons in a torus around the CC or the BN axis. As a result, although electron density is depleted at the bcp, the bond does not necessarily become weaker, since density increases around the periphery, a phenomenon named the “hole” effect. The dissimilarities between acetylene and iminoborane derivatives are primarily a consequence of the significant distortion of this torus in the latter, due to the large difference between the electronegativities of B and N, which leads to a large contribution of the Image ID:b702480k-t1.gifresonance structure in some cases. The “hole” effect is reflected in a reasonable correlation between the Laplacian of the electron density at the bcp and the BN bond length.

Dalton Trans, 3995-399 (2007)

Doi: 10.1039/b705703b

A theoretical and experimental study of the fluxional behaviour of molybdenum dihydrobis- and hydrotris-pyrazolylborates

The fluxional barrier of (dicarbonyl)[dihydrobis(3,5-dimethylpyrazol-1-yl)borato][η-(1,2,3)-2-methylpropen-1-yl]-molybdenum (1) has been measured and a complete assignment of its 1H,13C and 15N NMR signals has been carried out. Theoretical calculations at the B3LYP/LANL2DZ level including GIAO absolute shieldings (σ) have allowed to analyze the molecular contributions to the barrier as well as to assign some signals involved in the fluxional process.

Organomet. 25, 5638-5647 (2006)

Doi: 10.1021/om0608197

On the Existence of a-Agostic Bonds: Bonding Analyses of Titanium Alkyl Complexes

The main characteristics of the so-called “agostic bonding” in titanium complexes were computationally investigated in terms of diverse analyses, including the quantum theory of atoms in molecules (QTAIM) and the electron localization function (ELF). Computations on a set of titanium-based molecular models that were presumably able to present α- and β-like bonding showed a clear distinction between α- and β-bonding schemes. In view of the geometric, energetic, and electronic data of these molecules, we concluded that the geometries presenting a Ti···H approximation similar to α bonding are not the result of a bonding attraction. On the contrary, their origin arises from a short-range repulsion between the metal core and the lone pairs of the alkylidene group, the latter pivoting in its own plane around carbon and thereby allowing simultaneously a closer approach between the Ti and the C atoms and indirectly resulting in a short Ti···H distance. Additionally, we found that the electronic and geometric distortion of the C−H bond present in the agostic bonding are not univocally linked to this kind of bonding; instead, this originates in the close distance of the C−H bond and any metal center, independently of whether it is due to an agostic bond or not. Therefore, the lengthening of the C−H bond and the reduction of its electron density at the bond critical point should not be considered as indication of agostic bonding but as a side effect.

Chem. Phys. Lett. 429, 23-26 (2006)

Doi: 10.1016/j.cplett.2006.08.012

The structure of protonated HCP: A classical or non-classical ion?

The bond dissociation enthalpies in n-silanes and n-germanes have been studied with the G2MP2 and G3B3 computational methods. The results have been compared with those of the corresponding n-alkanes. The differences observed in the three families of compounds have been rationalized using the calculated hyperpolarization effect in each family. In addition, the conformational profiles of X2H6 molecules, X = C, Si, and Ge, have been compared.

Chem. Phys. Lett. 429, 58-61 (2006)

Doi: 10.1016/j.cplett.2006.07.098

Theoretical study of the bond energy in n-silanes and n-germanes: Comparison with n-alkanes

The bond dissociation enthalpies in n-silanes and n-germanes have been studied with the G2MP2 and G3B3 computational methods. The results have been compared with those of the corresponding n-alkanes. The differences observed in the three families of compounds have been rationalized using the calculated hyperpolarization effect in each family. In addition, the conformational profiles of X2H6 molecules, X = C, Si, and Ge, have been compared.

Chem. Phys. Lett. 425, 221-224 (2006)

Doi: 10.1016/j.cplett.2006.05.050

The carbon-carbon bond dissociation energy as a function of the chain length

All the C–C fragmentations of the linear n-alkanes, CnH2n+2, for n = 2 to n = 10 and the final C–H bond breaking have been calculated with the composite methods G2, G2MP2, G3, G3MP2B3 and G3B3. The results have been compared with the available experimental data.

Org. Biomol. Chem., 3 , 3905-3910 (2005)

Doi: 10.1039/b510535h

3-(2-Pyridyl)-[1,2,3]triazolo[1,5-a]pyridines. An experimental and theoretical (DFT) study of the ring-chain isomerization 

An experimental (1H NMR) and theoretical (DFT) study of the ring–chain–ring isomerization of 3-(2-pyridyl)-[1,2,3]triazolo[1,5-a]pyrid-7-yl derivatives (A) into 6-{[1,2,3]triazolo[1,5-a]pyrid-3-yl}-2-pyridyl derivatives (B) has been carried out. Based on the calculations, a mechanism of several steps will be proposed. The experimental results as well as the calculations lead to the conclusion that the A–B ratio depends on the electronic properties of the substituents.

J. Phys. Chem. A, 109, 4359-4365 (2005)

Doi: 10.1021/jp050857o

Periodic Trends in Bond Dissociation Energies. A Theoretical Study

Bond dissociation energies (BDEs) of all possible A−X single bonds involving the first- and second-row atoms, from Li to Cl, where the free valences are saturated by hydrogens, have been estimated through the use of the G3-theory and at the B3LYP/6-311+G(3df,2pd)//B3LYP/6-31G(2df,p) DFT level of theory. BDEs exhibit a periodical behavior. The A−X (A = Li, Be, B, Na, Mg, Al, and Si) BDEs show a steady increase along the first and the second row of the periodic table as a function of the atomic number Z(X). For A−X bonds involving electronegative atoms (A = C, N, O, F, P, S, and Cl) the bond energies achieve a maximum around Z(X) = 5. The same behavior is observed when BDEs are plotted against the electronegativity χ(X) of the atom X. Thus, for A−X bonds (A = Li, Be, B, Na, Mg, Al, Si), the BDEs for a fixed A increases, grosso modo, as the electronegativity differences between X and A increase, with some exceptions, which reflect the differences in the relaxation energies of the radicals produced upon the bond cleavage. A similar trend, albeit less pronounced, is found for single A−X bonds, where A = C, N, O, F, P, S, and Cl. However, there is an additional feature embodied in the enhancement of the strength of the A−boron bonds due to the ability of boron to act as a strong electron acceptor. The trends in bond lengths and charge densities at the bond critical points are in line with the aforementioned behavior.

Arkivoc, ix, 305-320 (2005).


Influence of protonation on the properties derived from electron density

The effect of protonation on the electron density of a series of small molecules has been studied by means of the quantum theory of the Atoms in Molecules (QTAIM) methodology. The nonuniform redistribution of the electron density upon protonation of the systems produces an increase of the atomic energy of the protonated heteroatom and a reduction of the atomic charge for the protonated oxygen, while an increase of the charge on the protonated nitrogens is observed. The total volume of the protonated systems is smaller than that of the corresponding neutral ones. Thus, all the groups that form the molecule suffer a reduction of the volume when the systems become protonated. The transferability of the methyl group properties between the molecules studied has been analyzed. 

The Chemistry of Cyclobutanes, pg. 177-211

Acidity and basicity of cyclobutanes

Book Chapter: The Chemistry of Cyclobutanes, pg. 177-211, Ed. Z. Rappoport and J.F. Liebaman, Patai Series: The Chemistry of Functional Groups

Struct. Chem. 16, 77-79 (2005)

Doi: 10.1007/s11224-005-1089-9

Polyynes vs. Cumulenes: Their Possible Use as Molecular Wires

Polyynes and cumulenes from 2–12 atoms have been calculated at the B3LYP/6-311++G** level and their energies compared using an isodesmic reaction. The chain length has been modeled empirically affording an equation that predicts substantial variation for long chains.


Tet. Lett. 45, 4127-4129 (2004)

Doi: 10.1016/j.tetlet.2004.03.141

Classical versus redox tautomerism: substituent effects on the keto/enol and sulfoxide/sulfenic acid equilibria

MP2/6-311+G** ab initio calculations have been carried out for a classical example of tautomerism, the keto/enol [RCOCH3/RC(OH)CH2] and an example of redox tautomerism, the sulfoxide/sulfenic acid [RS(O)H/RSOH]. Eleven R substituents have been examined. Both equilibria show proportional energies and similar dependence on the Swain–Lupton Fand Rparameters.

Struct. Chem. 15, 599-604 (2004)

Doi: 10.1007/s11224-004-0735-y

The Position of the Lone Pair in Amines: A Comparison Between Bader's AIM Approach and Pure Geometrical Considerations

The position of the lone pair of 17 amines has been located with the help of Bader’s AIM methodology. This position has been compared with those estimated with simple geometrical models. It appears that these models are convenient only for some symmetrical or undistorted cases.

Collect. Czech. Chem. Commun. 69, 2134-2146 (2004)

Doi: 10.1135/cccc20042134

Protic acidity of some aliphatic and alicyclic hydrocarbons in the gas phase and in solution. An empirical and computational link

Here we report the results of the study of a set of thirteen hydrocarbons R-H not leading to extensively charge-delocalized anions R- upon ionization, R-H → R- + H+. It essentially involves the following: (i) The computation at the G2 level of the changes in thermodynamic state functions for this process in the gas phase. The reliability of the computational method was further assessed by comparison with CCSD(T)/Aug-cc-pVTZ results. (ii) A direct comparison of the experimentally available thermodynamic (gas phase) and kinetic (solution) data pertaining to this reaction. (iii) A careful re-examination of the experimental data sets and the Brønsted-type relationships derived therefrom by using the computed thermodynamic data. This treatment suggests the existence of mechanistic features affecting the experimental data and indicates the need for further experimental and computational work.

Chem. Phys. Lett. 398, 560-563 (2004)

Doi: 10.1016/j.cplett.2004.09.127

The enormous apparent gas-phase acidity of cubylamine

The high acidity of cubylamine (1NH2) seems to originate in the release of strain energy attending the breaking of some C–C bonds in 1NH. This process is greatly facilitated by the strong stereoelectronic interactions in 1NH. The anionic species thus formed are less strained, and their corresponding conjugate acids seem unable to ‘borrow strength’ from the residual strain, at least within the time-scale of the FT ICR experiments.

Angew. Chem. Int. Ed. 42, 2281-2284 (2003)

Doi: 10.1002/anie.200350993

Stereoelectronic, Strain, and Medium Effects on the Protonation of Cubylamine, a Janus-like Base

The cubic structure of cubylamine (1NH2) is maintained on protonation in aqueous solution, but in the gas phase 1NH2 is irreversibly protonated on the carbon framework, leading to the formation of iminium ions through ring opening (see scheme).