Density functional theory has been used in combination with the correlation consistent and polarization consistent basis sets to investigate the structures and energetics for a series of first-row closed shell and several second-row molecules of potential importance in atmospheric chemistry. The impact of basis set choice upon molecular description has been examined, and irregular convergence of molecular properties with respect to increasing basis set size for several functionals and molecules has been observed. The possible reasons and solutions for this unexpected behavior including the effect of contraction and uncontraction, of the basis set diffuse sp basis functions, basis set superposition error (BSSE) and core-valence sets also have been examined.

Employing correlation consistent basis sets coupled with electronic structure methods has enabled accurate predictions of chemical properties for second- and third-row main group and transition metal molecular species. For third-row (Ga-Kr) molecules, the performance of the correlation consistent basis sets (cc-pVnZ, n=D, T, Q, 5) for computing energetic (e.g., atomization energies, ionization energies, electron and proton affinities) and structural properties using the ab initio coupled cluster method including single, double, and quasiperturbative triple excitations [CCSD(T)] and the B3LYP density functional method was examined. The impact of relativistic corrections on these molecular properties was determined utilizing the Douglas-Kroll (cc-pVnZ-DK) and pseudopotential (cc-pVnZ-PP) forms of the correlation consistent basis sets. This work was extended to the characterization of molecular properties of novel chemically bonded krypton species, including HKrCl, FKrCF3, FKrSiF3, FKrGeF3, FKrCCF, and FKrCCKrF, and provided the first evidence of krypton bonding to germanium and the first di-krypton system. For second-row (Al-Ar) species, the construction of the core-valence correlation consistent basis sets, cc-pCVnZ was reexamined, and a revised series, cc-pCV(n+d)Z, was developed as a complement to the augmented tight-d valence series, cc-pV(n+d)Z. Benchmark calculations were performed to show the utility of these new sets for second-row species. Finally, the correlation consistent basis sets were used to study the structural and spectroscopic properties of Au(CO)Cl, providing conclusive evidence that luminescence in the solid-state can be attributed to oligomeric species rather than to the monomer.

The development of the semi-empirical atomistic potential called the embedded atom method (EAM) has allowed for the efficient modeling of solid-state environments, at a lower computational cost than afforded by density functional theory (DFT). This offers the capability of EAM to model the energetics of solid-state phases of varying coordination, including defects, such as vacancies and self-interstitials. This dissertation highlights the development and application of two EAMs: a Ti potential constructed with the multi-state modified embedded atom method (MS-MEAM), and a Ni potential constructed with the fragment Hamiltonian (FH) method. Both potentials exhibit flexibility in the description of different solid-states phases and applications. This dissertation also outlines two applications of DFT. First, a study of structure and stability for solid-state forms of NixCy (in which x and y are integers) is investigated using plane-wave DFT. A ground state phase for Ni2C is elucidated and compared to known and hypothesized forms of NixCy. Also, a set of correlation consistent basis sets, previously constructed using the B3LYP and BLYP density functionals, are studied. They are compared to the well-known to the correlation consistent basis sets that were constructed with higher-level ab initio methodologies through computations of enthalpies of formation and combustion enthalpies. The computational accuracy with regard to experiment is reported.

In this thesis, the implementation of the incremental scheme for open-shell systems with unrestricted Hartree-Fock reference wave functions is described. The implemented scheme is tested within robustness and performance with respect to the accuracy in the energy and the computation times. New approaches are discussed to implement a fully automated incremental scheme in combination with the domain-specific basis set approximation. The alpha Domain Partitioning and Template Equalization are presented to handle unrestricted wave functions for the local correlation treatment. Both orbital schemes are analyzed with a test set of structures and reactions. As a further goal, the DSBSenv orbital basis sets and auxiliary basis sets are optimized to be used as environmental basis in the domain-specific basis set approach. The performance with respect to the accuracy and computation times is analyzed with a test set of structures and reactions. In another project, a scheme for the optimization of auxiliary basis sets for uranium is presented. This scheme was used to optimize the MP2Fit auxiliary basis sets for uranium. These auxiliary basis enable density fitting in quantum chemical methods and the application of the incremental scheme for systems containing uranium. Another project was the systematical analysis of the binding energies of four water dodecamers. The incremental scheme in combination with the CCSD(T) and CCSD(T)(F12*) method were used to calculate benchmark energies for these large clusters.

The development of the correlation consistent basis sets, cc-pVnZ (where n = D, T, Q, etc.) have allowed for the systematic elucidation of the intrinsic accuracy of ab initio quantum chemical methods. In density functional theory (DFT), where the cc-pVnZ basis sets are not necessarily optimal in their current form, the elucidation of the intrinsic accuracy of DFT methods cannot always be accomplished. This dissertation outlines investigations into the basis set requirements for DFT and how the intrinsic accuracy of DFT methods may be determined with a prescription involving recontraction of the cc-pVnZ basis sets for specific density functionals. Next, the development and benchmarks of a set of cc-pVnZ basis sets designed for the s-block atoms lithium, beryllium, sodium, and magnesium are presented. Computed atomic and molecular properties agree well with reliable experimental data, demonstrating the accuracy of these new s-block basis sets. In addition to the development of cc-pVnZ basis sets, the development of a new, efficient formulism of the correlation consistent Composite Approach (ccCA) using the resolution of the identity (RI) approximation is employed. The new formulism, denoted 'RI-ccCA,' has marked efficiency in terms of computational time and storage, compared with the ccCA formulism, without the introduction of significant error. Finally, this dissertation reports three separate investigations of the properties of FOOF-like, germanium arsenide, and silicon hydride/halide molecules using high accuracy ab initio methods and the cc-pVnZ basis sets.

In this thesis, the implementation of the incremental scheme for open-shell systems with unrestricted Hartree-Fock reference wave functions is described. The implemented scheme is tested within robustness and performance with respect to the accuracy in the energy and the computation times. New approaches are discussed to implement a fully automated incremental scheme in combination with the domain-specific basis set approximation. The alpha Domain Partitioning and Template Equalization are presented to handle unrestricted wave functions for the local correlation treatment. Both orbital schemes are analyzed with a test set of structures and reactions. As a further goal, the DSBSenv orbital basis sets and auxiliary basis sets are optimized to be used as environmental basis in the domain-specific basis set approach. The performance with respect to the accuracy and computation times is analyzed with a test set of structures and reactions. In another project, a scheme for the optimization of auxiliary basis sets for uranium is presented. This scheme was used to optimize the MP2Fit auxiliary basis sets for uranium. These auxiliary basis enable density fitting in quantum chemical methods and the application of the incremental scheme for systems containing uranium. Another project was the systematical analysis of the binding energies of four water dodecamers. The incremental scheme in combination with the CCSD(T) and CCSD(T)(F12*) method were used to calculate benchmark energies for these large clusters.