Nuclear magnetic resonance chemical shifts with the statistical average of orbital-dependent model potentials in Kohn-Sham density functional theory

The performance of the SAOP potential for the calculation of NMR chemical shifts was evaluated. SAOP results show considerable improvement with respect to previous potentials, like VWN or BP86, at least for the carbon, nitrogen, oxygen, and fluorine chemical shifts. Furthermore, a few NMR calculations carried out on third period atoms (S, P, and Cl) improved when using the SAOP potential

Modelling of Mercury’s surface composition and remote detection from the orbit with the BepiColombo Mercury Planetary Orbiter

It can be assumed that the composition of Mercury’s thin gas envelope (exosphere) is related to the composition of the planets crustal materials. If this relationship is true, then inferences regarding the bulk chemistry of the planet might be made from a thorough exospheric study. The most vexing of all unsolved problems is the uncertainty in the source of each component. Historically, it has been believed that H and He come primarily from the solar wind, while Na and K originate from volatilized materials partitioned between Mercury’s crust and meteoritic impactors. The processes that eject atoms and molecules into the exosphere of Mercury are generally considered to be thermal vaporization, photonstimulated desorption (PSD), impact vaporization, and ion sputtering. Each of these processes has its own temporal and spatial dependence. The exosphere is strongly influenced by Mercury’s highly elliptical orbit and rapid orbital speed. As a consequence the surface undergoes large fluctuations in temperature and experiences differences of insolation with longitude. We will discuss these processes but focus more on the expected surface composition and solar wind particle sputtering which releases material like Ca and other elements from the surface minerals and discuss the relevance of composition modelling

Selected configuration interaction with truncation energy error and application to the Ne atom

Selected configuration interaction (SCI) for atomic and molecular electronic structure calculations is reformulated in a general framework encompassing all CI methods. The linked cluster expansion is used as an intermediate device to approximate CI coefficients BK of disconnected configurations (those that can be expressed as products of combinations of singly and doubly excited ones) in terms of CI coefficients of lower-excited configurations where each K is a linear combination of configuration-state-functions (CSFs) over all degenerate elements of K. Disconnected configurations up to sextuply excited ones are selected by Brown's energy formula, ΔEK=(E-HKK)BK2/(1-BK2), with BK determined from coefficients of singly and doubly excited configurations. The truncation energy error from disconnected configurations, Δdis, is approximated by the sum of ΔEKS of all discarded Ks. The remaining (connected) configurations are selected by thresholds based on natural orbital concepts. Given a model CI space M, a usual upper bound ES is computed by CI in a selected space S, and EM=E S+ΔEdis+δE, where δE is a residual error which can be calculated by well-defined sensitivity analyses. An SCI calculation on Ne ground state featuring 1077 orbitals is presented. Convergence to within near spectroscopic accuracy (0.5 cm-1) is achieved in a model space M of 1.4× 109 CSFs (1.1 × 1012 determinants) containing up to quadruply excited CSFs. Accurate energy contributions of quintuples and sextuples in a model space of 6.5 × 1012 CSFs are obtained. The impact of SCI on various orbital methods is discussed. Since ΔEdis can readily be calculated for very large basis sets without the need of a CI calculation, it can be used to estimate the orbital basis incompleteness error. A method for precise and efficient evaluation of ES is taken up in a companion paper

How does basis set superposition error change the potential surfaces for hydrogen-bonded dimers?

We describe a simple method to automate the geometric optimization of molecular orbital calculations of supermolecules on potential surfaces that are corrected for basis set superposition error using the counterpoise (CP) method. This method is applied to the H-bonding complexes HF/HCN, HF/H2O, and HCCH/H2O using the 6-31G(d,p) and D95 + + (d,p) basis sets at both the Hartree-Fock and second-order Møller-Plesset levels. We report the interaction energies, geometries, and vibrational frequencies of these complexes on the CP-optimized surfaces; and compare them with similar values calculated using traditional methods, including the (more traditional) single point CP correction. Upon optimization on the CP-corrected surface, the interaction energies become more negative (before vibrational corrections) and the H-bonding stretching vibrations decrease in all cases. The extent of the effects vary from extremely small to quite large depending on the complex and the calculational method. The relative magnitudes of the vibrational corrections cannot be predicted from the H-bond stretching frequencies alone

Families of periodic horseshoe orbits in the restricted three-body problem

We compute families of symmetric periodic horseshoe orbits in the restricted three-body problem. Both the planar and three-dimensional cases are considered and several families are found.We describe how these families are organized as well as the behavior along and among the families of parameters such as the Jacobi constant or the eccentricity. We also determine the stability properties of individual orbits along the families. Interestingly, we find stable horseshoe-shaped orbit up to the quite high inclination of 17◦

Rhodium(I) catalyzed [2+2+2] cycloaddition reactions: experimental and theoretical studies

The [2+2+2] cycloaddition reaction involves the formation of three carbon-carbon bonds in one single step using alkynes, alkenes, nitriles, carbonyls and other unsaturated reagents as reactants. This is one of the most elegant methods for the construction of polycyclic aromatic compounds and heteroaromatic, which have important academic and industrial uses. The thesis is divided into ten chapters including six related publications. The first study based on the Wilkinson’s catalyst, RhCl(PPh3)3, compares the reaction mechanism of the [2+2+2] cycloaddition process of acetylene with the cycloaddition obtained for the model of the complex, RhCl(PH3)3. In an attempt to reduce computational costs in DFT studies, this research project aimed to substitute PPh3 ligands for PH3, despite the electronic and steric effects produced by PPh3 ligands being significantly different to those created by PH3 ones. In this first study, detailed theoretical calculations were performed to determine the reaction mechanism of the two complexes. Despite some differences being detected, it was found that modelling PPh3 by PH3 in the catalyst helps to reduce the computational cost significantly while at the same time providing qualitatively acceptable results. Taking into account the results obtained in this earlier study, the model of the Wilkinson’s catalyst, RhCl(PH3)3, was applied to study different [2+2+2] cycloaddition reactions with unsaturated systems conducted in the laboratory. Our research group found that in the case of totally closed systems, specifically 15- and 25-membered azamacrocycles can afford benzenic compounds, except in the case of 20-membered azamacrocycle (20-MAA) which was inactive with the Wilkinson’s catalyst. In this study, theoretical calculations allowed to determine the origin of the different reactivity of the 20-MAA, where it was found that the activation barrier of the oxidative addition of two alkynes is higher than those obtained for the 15- and 25-membered macrocycles. This barrier was attributed primarily to the interaction energy, which corresponds to the energy that is released when the two deformed reagents interact in the transition state. The main factor that helped to provide an explanation to the different reactivity observed was that the 20-MAA had a more stable and delocalized HOMO orbital in the oxidative addition step. Moreover, we observed that the formation of a strained ten-membered ring during the cycloaddition of 20-MAA presents significant steric hindrance. Furthermore, in Chapter 5, an electrochemical study is presented in collaboration with Prof. Anny Jutand from Paris. This work allowed studying the main steps of the catalytic cycle of the [2+2+2] cycloaddition reaction between diynes with a monoalkyne. First kinetic data were obtained of the [2+2+2] cycloaddition process catalyzed by the Wilkinson’s catalyst, where it was observed that the rate-determining step of the reaction can change depending on the structure of the starting reagents. In the case of the [2+2+2] cycloaddition reaction involving two alkynes and one alkene in the same molecule (enediynes), it is well known that the oxidative coupling may occur between two alkynes giving the corresponding metallacyclopentadiene, or between one alkyne and the alkene affording the metallacyclopentene complex. Wilkinson’s model was used in DFT calculations to analyze the different factors that may influence in the reaction mechanism. Here it was observed that the cyclic enediynes always prefer the oxidative coupling between two alkynes moieties, while the acyclic cases have different preferences depending on the linker and the substituents used in the alkynes. Moreover, the Wilkinson’s model was used to explain the experimental results achieved in Chapter 7 where the [2+2+2] cycloaddition reaction of enediynes is studied varying the position of the double bond in the starting reagent. It was observed that enediynes type yne-ene-yne preferred the standard [2+2+2] cycloaddition reaction, while enediynes type yne-yne-ene suffered β-hydride elimination followed a reductive elimination of Wilkinson’s catalyst giving cyclohexadiene compounds, which are isomers from those that would be obtained through standard [2+2+2] cycloaddition reactions. Finally, the last chapter of this thesis is based on the use of DFT calculations to determine the reaction mechanism when the macrocycles are treated with transition metals that are inactive to the [2+2+2] cycloaddition reaction, but which are thermally active leading to new polycyclic compounds. Thus, a domino process was described combining an ene reaction and a Diels-Alder cycloaddition.

Synthesis, structure and reactivity of novel CuI, CuII and CuIII complexes containing triaza and hexaaza macrocyclic ligands

El treball de tesi s'emmarca dins del camp de la bioinorgànica, disciplina que estudia les propietats estructurals i de reactivitat dels centres actius dels enzims, servint-se de models síntètics de baix pes molecular per tal d'intentar reproduïr la reactivitat presentada per l'enzim i conèixer els mecanismes de reacció a nivell molecular que tenen lloc en els processos biològics.1 Més concretament el treball posa especial èmfasi en els processos d'activació d'oxigen molecular que tenen lloc en les metaloproteïnes de Coure del Tipus 3, com són l'hemocianina i la tirosinasa, ambdues presentant un complex dinuclear de Cu(I)) en el centre actiu de la forma reduïda, capaç d'activar l'O2 cap a espècies de tipus peròxid.2 Un altre camp d'interès ha estat l'estudi dels processos d'activació d'enllaços C-H no activats en hidrocarburs, tant per la seva importàcia a nivell industrial com per comprendre els mecanismes intrínsecs d'aquesta activació a través de metalls de trancisió.3,4 Durant el treball de tesi presentat s'ha desenvolupat la síntesi de nous complexes de Coure(I), Coure(II) y Cu(III) utilitzant lligands macrocíclics de tipus triaza i hexaaza, i s'han estudiat la seves propietats estructurals així com la seva reactivitat. La reacció dels lligands triazacíclics H32m, H2Me33m i H33m amb sals de coure(II) dóna lloc a una reacció de desproporció de Cu(II) per obtenir-se en quantitats equimolars un complex organometàl·lic de Cu(III) i un complex de Cu(I). La caracterizació estructural exhaustiva dels complexes del tipus aryl-Cu(III) evidencia la formació d'un enllaç organometàl·lic entre l'àtom de Cu(III) i el carboni més próxim de l'anell aromàtic del lligand. Aquesta reacció, a més de representar una nova forma de desproporció en la química del Cu, suposa l'activació d'un enllaç C-H aromàtic a temperatura ambient que, mitjançant l'estudi cinètic d'aquesta desproporció per espectroscòpia UV-Vis, dels càlcul de l'efecte cinètic isotòpic utilitzant el lligand deuterat en el C-H de l'anell, juntament amb el recolzament teòrics dels càlculs DFT per a la optimització de geometries d'intermedis de reacció, ens permeten proposar un mecanisme de reacció pel nostre sistema, on l'activació de l'enllaç C-H aromàtic transcorre per la formació d'un enllaç de tipus agòstic C-H ? Cu(II),5 seguit de la desprotonació del C-H aromàtic per acció d'una base i posterior transferència electrònica per obtenir el complex organometàlic de Cu(III) i el complex de de Cu(I). En quant a la reactivitat d'aquests complexes organometàl·lics aryl-Cu(III) s'ha observat que una base en medi aquós causa la inestabilitat d'aquests compostos, evolucionant cap a la inserció d'un àtom d'oxigen sobre la posició activada de l'anell aromàtic, per a donar lloc a un complex dinuclear de Cu(II) amb dos grups fenoxo actuant de pont entre els àtoms metàl·lics. La reacció transcorre per un intermedi colorejat, caracteritzat com el complex ayl-Cu(III) monodesprotonat en una de les seves amines benzíliques, els quals s'observen igualment en la reacció dels correponents complexos de Cu(I) amb oxigen molecular (O2). És en els nostres sistemes en els quals es descriu per primera vegada la participació d'intermedis organometàl·lics Cu(III)-C en processos d'hidroxilació aromàtica, tals com el desenvolupat per l'enzim tirosinasa o per alguns dels seus models químics de síntesi.6,7,8 S'han estudiat les propietats magnètiques dels quatre bis(fenoxo)complexes de Cu(II) descrits, obtenint-se uns acoplaments de tipus antiferromagnètic o ferromagnètic de diversa magnitud, depenent del solapament orbitalari a l'enllaç Cu-O, a través del qual es produeix el superintercanvi. Nous complexos de Cu(I) sintetitzats amb lligands hexaazamacrocíclics han estat estudiats, i posant especial èmfasi a la seva reactivitat respecta a l'activació d'oxigen molecular (O2). S'ha observat una reactivitat diferenciada segons la concentració de complex de Cu(I) utilitzada, de manera que a altes concentracions s'obté un carbonato complex tetranuclear de Cu(II) per fixació de CO2 atmosfèric, mentre que a baixes concentracions s'observa la hidroxilació aromàtica intramolecular d'un dels anells benzílics del lligand, reacció que presumiblement transcorre per atac electrofílic d'un peroxo complex intermedi sobre el sistema ? de l'anell.6 Els resultats obtinguts en aquest treball ens mostren la facilitat per activar enllaços C-H aromàtics per metalls de transició de la primera sèrie (Cu, Ni) quan aquests estan suficientment pròxims a l'enllaç C-H, en unes condicions de reacció molt suaus (1atm., temperatura ambient). Els nous complexos organometàl·lics Aryl-Cu(III) són el producte d'una nova reacció de desproporció de Cu(II), així com un posició aromàtica activada que podria ser el punt de partida per l'estudi de funcionalització selectiva d'aquests grups aromàtics.