Theoretical studies of pi-pi interactions on several dimers of curved polycyclic aromatic systems have been carried out. In the first part, dispersion corrected density functional theory methods (DC-DFT) were used to evaluate the basis set superposition errors (BSSE) in dispersion interactions of the corannulene dimer, and the accuracy of the calculations using DC-DFT methods was compared with high level benchmark calculations. In these calculations, Grimme’s B97D DC-DFT method provided reasonably accurate results with the benchmark calculations. In addition, BSSE obtained with the B97D method along with cc-pVQZ basis set was negligible. Furthermore, a series of calculations were carried out to obtain the pi-pi interaction energy and most stable conformation for the sumanene dimer system. In these calculations, Grimme’s B97D method was used. The potential energy minimum of the sumanene dimer was determined as the concave-convex stacked arrangement with one monomer unit rotated to 60°. The binding energy of the dimer was found to be 19.34 kcal/mol with a 3.72 angstrom distance between two monomer units. Dimers of three different heterosumanenes along with the parent sumanene were also studied. In this set of calculations, two different concave-convex dimer motifs were chosen, eclipsed and staggered (60° rotated). For all the heterosumanenes, as well as the parent sumanene, the staggered conformation is the most stable geometry. The parent sumanene had the highest binding energy. The –NH substituted sumanene produced the second highest binding energy, while the –O analog was the weakest bonded dimer. Finally, dispersion calculations were carried out for the planar aromatic compound of triphenylene. The pi-system of the dimer was distorted by rotating one monomer unit around the principle axis and parallel displacing one monomer unit relative to the other one. Among the rotational dimers, the 39° rotated dimer was the minimum energy conformation. Interaction energy of that dimer was 14.42 kcal/mol with 3.40 angstrom separation between monomers at the B97D/cc-pVQZ level. The parallel displaced minimum energy dimer has a binding energy about 1.0 kcal/mol smaller than the rotational minimum energy geometry.

Les polyphénols naturels forment des complexes non-covalents dans lesquels le π-stacking et les liaisons hydrogène jouent un rôle clé dans la stabilisation. Les calculs DFT incluant la dispersion (DFT-D), la description des processus d'agrégation non-covalente de produits naturels devient fiable. Dans ce travail, les méthodes DFT-D sont appliquées à i) la compréhension de la biosynthèse stéréo- et régio-sélective des oligostilbenoïdes, ii) la prédiction de l'agrégation des antioxydants naturels au sein de la membrane bicouche lipidique, qui pourrait rationaliser la synergie de la vitamine E, la vitamine C et polyphénols dans leur action antioxydante, et iii) la modulation des propriétés optiques de dérivés de chalcones. / Natural polyphenols form non-covalent complexes in which π-stacking and H-bonding play a key stabilizing role. The dispersion-corrected DFT calculations have paved the way towards reliable description of aggregation processes of natural products. In this work, these methods are applied at i) understanding of stereo- and regio-selective oligostilbenoids biosynthesis; ii) predicting natural antioxidant aggregation within lipid bilayer membrane, which may allow rationalizing the synergism of vitamin E, vitamin C and polyphenols in their antioxidant action; and iii) modulating optical properties of chalcone derivatives.

Dans le cadre de la recherche de vecteurs d'énergie “propres”, le borazane et ses dérivés amine-boranes sont devenus des candidats intéressants en tant que matériaux de stockage de l'hydrogène en raison de leur pourcentage massique relativement élevé en hydrogène (19,6% pour borazane) et de la réversibilité potentielle de la réaction de déshydrogénation. Pour des applications réelles, le contrôle des réactions se produisant à la température ambiante est fondamental. Dans ce contexte, la compréhension du processus de la déshydrogénation/déshydrocouplage catalytique de l'amine-borane apparaît comme un élément important. Dans cette thèse, les catalyseurs de types métallocènes du Groupe IV (Cp2M, M = Ti, Zr et Hf) sont étudiés en détail. Le déshydrocouplage de HMe2N·BH3 catalysé par le titanocène a été étudié à la fois expérimentalement et théoriquement mais aucun accord n'avait été atteint auparavant. Dans ce travail, les caractérisations systématiques des interactions 3-centre 2-électron M···H-B impliquées dans les intermédiaires réactionnels ont été réalisées avec des approches topologiques QTAIM et ELF. Par la suite, des mécanismes réactionnels détaillés ont été étudiés. Les résultats théoriques ont démontré que la méthode DFT corrigée avec la dispersion (DFT-D) étaient nécessaire et suffisantes pour une description énergétique correcte des chemins réactionnels. Mon travail a également permis l'identification d'un complexe de van der Waals jouant un rôle clé dans le mécanisme réactionnel en accord avec les observations expérimentales. / With the increasing demand of clean energy carriers, ammonia borane and its related amine-borane compounds have emerged as attractive candidates for hydrogen storage materials due to their relatively high weight percentage of available hydrogen (19.6% for ammonia borane) as well as the potential reversibility for the hydrogen release reactions. Actual applications would benefit from controlled reactions occurring close to room-temperature. In this context, catalytic dehydrogenation/dehydrocoupling of amine-borane appears as a promising solution. In this thesis the Group IV metallocene (Cp2M, M = Ti, Zr and Hf) are mainly discussed. The dehydrocoupling of HMe2N·BH3 catalyzed by titanocene was investigated both experimentally and theoretically but no agreement were reached. In this work, systematic characterization of M···H-B 3-center 2-electron interactions involved in reaction intermediates were carried out with QTAIM and ELF topological approaches. Afterwards, detailed mechanisms were further studied. Computational results have demonstrated that the dispersion corrected DFT (DFT-D) method was indispensable for a correct enegetic prediction for reaction pathways. The identification of a van der Waals complexe also plays a central role for a reaction mechanism with good agreement with experimental observations.

On account of the high temperatures and pressures experienced by energetic materials during deflagration and detonation, it is important to know not only the physical properties of these materials at ambient temperatures and pressures, but also to understand how their structure and properties are affected by extreme conditions. Combined computational and experimental investigations of the effects of high pressures on the structure and properties of several energetic materials are described herein. A comparison of the performances of different pseudopotentials and density functional theory (DFT) dispersion correction schemes in calculating crystal geometries and vibrational frequencies of crystalline ammonium perchlorate at high pressure is described. The results highlight the fact that care must be taken when choosing pseudopotentials for high-pressure studies. A comprehensive comparison of calculated vibrational modes (including symmetry) with experiment has been performed, with the frequencies of all internal modes predicted to lie within 5% of experimental values. This study established that no significant improvements in the calculation of crystal geometries of ammonium perchlorate are obtained by employing DFT-D corrections. The enthalpy of fusion (ΔHfus) of the highly metastable β-form of RDX (cyclotrimethylenetrinitramine) was determined to be 12.63 ± 0.28 kJ mol-1. DFT-D calculations of the lattice energies of the α- and β-forms of RDX are described. Furthermore, the response of the lattice parameters and unit-cell volumes to pressure for the α-, γ- and ε-forms of RDX calculated using DFT-D are in very good agreement with experimental data. Phonon calculations provide good agreement with vibrational frequencies obtained from Raman spectroscopy, and a predicted inelastic neutron scattering (INS) spectrum of α-RDX shows excellent agreement with experimental INS data recorded as part of this study. The results of the high-pressure phonon calculations have been used to show that the heat capacities of the α-, γ- and ε- forms of RDX are only weakly affected by pressure. DFT-D calculations have been utilised to describe accurately the structure and properties of both β-HMX (Cyclotetramethylenetetranitramine) and α-FOX-7 (1,1-Diamino-2,2-dinitroethylene) as a function of pressure. This work presents data for the experimental hydrostatic compression of both deuterated β-HMX and α-FOX-7 performed using neutron powder diffraction at the ISIS Neutron and Muon facility, in addition to experimental determinations of the INS spectra of both β-HMX and α-FOX-7. The DFT-D hydrostatic compression studies for both materials reproduce the experimental compression trends. Furthermore, the calculated vibrational properties as a function of pressure were in very good agreement with available experimental data. The results of the phonon calculations were then used to predict the effect of pressure on the heat capacities of β-HMX and α-FOX-7. These predictions suggest a very weak pressure dependence of heat capacities (approximately -1 J K-1 mol-1 GPa-1) for these materials. This work demonstrates that the DFT-D model performs extremely well over a range of conditions, and is able to describe accurately intramolecular and intermolecular interactions, and thus the structure and properties of organic molecular nitramine crystals. The computational model was therefore used to predict the high-pressure hydrostatic compression behaviour of a related nitramine, CL-20 (2,4,6,8,10,12-Hexanitro-2,4,6,8,10,12-hexaazaisowurtzitane), the results of which highlighted possible discrepancies in the experimental high-pressure X-ray diffraction data recorded for ε-CL-20. This prompted a high-pressure neutron powder diffraction study, which showed good agreement with the computational results, thereby highlighting radiation damage in the X-ray experiments.