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.

London or dispersion interactions are weak van der Waals (vdW) interactions. They are important in determining the structure and properties of many chemical and biochemical systems. In this thesis, an optimizer using the nonempirical generalized gradient approximation (GGA) functional PW86+PBE+XDM, to capture van der Waals interactions, is presented. The work in this thesis covers the assessment of a variety of basis sets for their ability to reproduce accurate GGA repulsive and binding energies. Selected basis sets were then used to compute binding energies of 65 vdW complexes at equilibrium. This functional was also tested for binding energies of two sets of vdW complexes at distorted geometries. The last part deals with forces to investigate their accuracy using PW86+PBE+XDM in order to build an optimizer for vdW complexes using a nonempirical DFT method. Eventually, after confirming a high reproducibility of the optimizer on the geometries and binding energies, it was used in two biologically relevant applications. This optimizer is a unique tool to compute deformation energies with a nonempirical DFT method. The second part of this thesis covers a biologically relevant application where a conventional DFT is used. This application is related to the carrier of the genetic codes in living cells, DNA. DNA undergoes harmful mutations under external perturbations such as applied external electric fields. In this study, DNA base pairs were first mimicked by a simpler model, namely, the formic acid dimer. The effect of applied external electric fields on the geometries of the formic acid dimer is studied. The effect of these applied fields on the potential energy surface, the barrier height and the frequency of the double proton transfer in the formic acid dimer are also investigated. The study was then repeated on DNA base pairs to study the effect of an external applied electric field on the tunneling corrected rate constants of the double proton transfer reactions in AT and GC.

Engineering new solid propellant materials requires optimization of several factors, to include energy density, burn rate, sensitivity, and environmental impact. Equally important is the need for materials that will maintain their mechanical properties and thermal stability during long periods of storage. The nanoscale materials considered in this dissertation are proposed metal additives that may enhance energy density and improve combustion in a composite rocket motor. Density Functional Theory methods are used to determine cluster geometries, bond strengths, and energy densities. The ground-state geometries and electron affinities (EAs) for MnxO?: x = 3, 4, y = 1, 2 clusters were calculated with GGA, and estimates for the vertical detachment energies compare well with experimental results. It was found that the presence of oxygen influences the overall cluster moment and spin configuration, stabilizing ferrimagnetic and antiferromagnetic isomers. The calculated EAs range from 1.29-1.84 eV, which is considerably lower than the 3.0-5.0 eV EAs characteristic of current propellant oxidizers. Their use as solid propellant additives is limited. The structures and bonding of a range of Al-cyclopentadienyl cluster compounds were studied with multilayer quantum mechanics/molecular mechanics (QM:MM) methods. The organometallic Al-ligand bonds are generally 55-85 kcal/mol and are much stronger than Al-Al interactions. This suggests that thermal decomposition in these clusters will proceed via the loss of surface metal-ligand units. The energy density of the large clusters is calculated to be nearly 60% that of pure aluminum. These organometallic cluster systems may provide a route to extremely rapid Al combustion in solid rocket motors. Lastly, the properties of COOH-terminated passivating agents were modeled with the GPW method. It is confirmed that fluorinated polymers bind to both Al(111) and Al(100) at two Al surface sites. The oligomers HCOOH, CH3CH2COOH, and CF3CF2COOH chemisorb onto Al(111) with adsorption energies of 10-45 kcal/mol. The preferred contact angle for the organic chains is 65-85 degrees, and adsorption energy weakens slightly with increasing chain length. Despite their relatively weak adsorption energies, fluorinated polymers have elevated melting temperatures, making them good passivation materials for micron-scale Al fuel particles.

The present work describes new aspects of organic and supramolecular chemistry. The scientific contribution consists of two parts, which focus on the development of receptors for the sulfate anion and quantitative assessment of stacking interactions between an anthracene dye and nucleobases in an aqueous solution. In Chapter 1, basic concepts concerning supramolecular chemistry and recognition of cations and anions are discussed, as well as modern methods for the determination of binding constants. Particular attention is paid to fluorescence sensing of ions and underlying mechanisms of binding-induced fluorescence responses. Chapter 2 is dedicated to the design and synthesis of new fluorescent sulfate receptors functioning in aqueous solution. After a short review of the most effective sulfate receptors/probes created so far, a new design of PET probes for sulfate sensing is presented. The syntheses and anion binding properties of new compounds are described. The experimental data obtained for the receptors are discussed in detail to reveal the origin of high selectivity towards sulfate. Chapter 3 explores the importance of nucleobase–arene stacking interactions in recognition of nucleotides by synthetic receptors. Various experimental and theoretical approaches are presented to assess dispersion interactions between aromatic rings and nucleobases in the receptor–nucleotide complexes.

Most density functionals lack to correctly account for long-range London dispersion interactions, and numerous a posteriori correction schemes have been proposed in recent years. In van der Waals structures, the interlayer distance controls the proximity effect on the electronic structure, and the interlayer interaction energy indicates the possibility to mechanically exfoliate a layered material. For upcoming twisted van der Waals heterostructures, a reliable but efficient and scalable theoretical scheme to correctly predict the interlayer distance is required. Therefore, the performance of a series of popular London dispersion corrections combined with computationally affordable density functionals is validated. As reference data, the experimental interlayer distance of layered bulk materials is used, and corresponding interlayer interaction energies are calculated using the random phase approximation. We demonstrate that the SCAN-rVV10 and PBE-rVV10L functionals predict interlayer interaction energies and interlayer distances of the studied layered systems within the range of the defined error limits of 10 meV per atom and 0.12 Å, respectively. Semi-empirical and empirical dispersion-corrected functionals show significantly larger error bars, with PBE+dDsC performing best with comparable quality of geometries, but with higher interlayer interaction energy error limits of ≈20 meV per atom.