First principles based multiparadigm modeling of electronic structures and dynamics

Electronic structures and dynamics are the key to linking the material composition and structure to functionality and performance.

An essential issue in developing semiconductor devices for photovoltaics is to design materials with optimal band gaps and relative positioning of band levels. Approximate DFT methods have been justified to predict band gaps from KS/GKS eigenvalues, but the accuracy is decisively dependent on the choice of XC functionals. We show here for CuInSe₂ and CuGaSe₂, the parent compounds of the promising CIGS solar cells, conventional LDA and GGA obtain gaps of 0.0-0.01 and 0.02-0.24 eV (versus experimental values of 1.04 and 1.67 eV), while the historically first global hybrid functional, B3PW91, is surprisingly the best, with band gaps of 1.07 and 1.58 eV. Furthermore, we show that for 27 related binary and ternary semiconductors, B3PW91 predicts gaps with a MAD of only 0.09 eV, which is substantially better than all modern hybrid functionals, including B3LYP (MAD of 0.19 eV) and screened hybrid functional HSE06 (MAD of 0.18 eV).

The laboratory performance of CIGS solar cells (> 20% efficiency) makes them promising candidate photovoltaic devices. However, there remains little understanding of how defects at the CIGS/CdS interface affect the band offsets and interfacial energies, and hence the performance of manufactured devices. To determine these relationships, we use the B3PW91 hybrid functional of DFT with the AEP method that we validate to provide very accurate descriptions of both band gaps and band offsets. This confirms the weak dependence of band offsets on surface orientation observed experimentally. We predict that the CBO of perfect CuInSe₂/CdS interface is large, 0.79 eV, which would dramatically degrade performance. Moreover we show that band gap widening induced by Ga adjusts only the VBO, and we find that Cd impurities do not significantly affect the CBO. Thus we show that Cu vacancies at the interface play the key role in enabling the tunability of CBO. We predict that Na further improves the CBO through electrostatically elevating the valence levels to decrease the CBO, explaining the observed essential role of Na for high performance. Moreover we find that K leads to a dramatic decrease in the CBO to 0.05 eV, much better than Na. We suggest that the efficiency of CIGS devices might be improved substantially by tuning the ratio of Na to K, with the improved phase stability of Na balancing phase instability from K. All these defects reduce interfacial stability slightly, but not significantly.

A number of exotic structures have been formed through high pressure chemistry, but applications have been hindered by difficulties in recovering the high pressure phase to ambient conditions (i.e., one atmosphere and room temperature). Here we use dispersion-corrected DFT (PBE-ulg flavor) to predict that above 60 GPa the most stable form of N₂O (the laughing gas in its molecular form) is a 1D polymer with an all-nitrogen backbone analogous to cis-polyacetylene in which alternate N are bonded (ionic covalent) to O. The analogous trans-polymer is only 0.03-0.10 eV/molecular unit less stable. Upon relaxation to ambient conditions both polymers relax below 14 GPa to the same stable non-planar trans-polymer, accompanied by possible electronic structure transitions. The predicted phonon spectrum and dissociation kinetics validate the stability of this trans-poly-NNO at ambient conditions, which has potential applications as a new type of conducting polymer with all-nitrogen chains and as a high-energy oxidizer for rocket propulsion. This work illustrates in silico materials discovery particularly in the realm of extreme conditions.

Modeling non-adiabatic electron dynamics has been a long-standing challenge for computational chemistry and materials science, and the eFF method presents a cost-efficient alternative. However, due to the deficiency of FSG representation, eFF is limited to low-Z elements with electrons of predominant s-character. To overcome this, we introduce a formal set of ECP extensions that enable accurate description of p-block elements. The extensions consist of a model representing the core electrons with the nucleus as a single pseudo particle represented by FSG, interacting with valence electrons through ECPs. We demonstrate and validate the ECP extensions for complex bonding structures, geometries, and energetics of systems with p-block character (C, O, Al, Si) and apply them to study materials under extreme mechanical loading conditions.

Despite its success, the eFF framework has some limitations, originated from both the design of Pauli potentials and the FSG representation. To overcome these, we develop a new framework of two-level hierarchy that is a more rigorous and accurate successor to the eFF method. The fundamental level, GHA-QM, is based on a new set of Pauli potentials that renders exact QM level of accuracy for any FSG represented electron systems. To achieve this, we start with using exactly derived energy expressions for the same spin electron pair, and fitting a simple functional form, inspired by DFT, against open singlet electron pair curves (H₂ systems). Symmetric and asymmetric scaling factors are then introduced at this level to recover the QM total energies of multiple electron pair systems from the sum of local interactions. To complement the imperfect FSG representation, the AMPERE extension is implemented, and aims at embedding the interactions associated with both the cusp condition and explicit nodal structures. The whole GHA-QM+AMPERE framework is tested on H element, and the preliminary results are promising.

Studies of the representation of interatomic distances in molecules and crystal as sum of atomic radii

The electron diffraction investigation of the following compounds has been carried out: sulfur, sulfur nitride, realgar, arsenic trisulfide, spiropentane, dimethyltrisulfide, cis and trans lewisite, methylal, and ethylene glycol.

The crystal structures of the following salts have been determined by x-ray diffraction: silver molybdateand hydrazinium dichloride.

Suggested revisions of the covalent radii for B, Si, P, Ge, As, Sn, Sb, and Pb have been made, and values for the covalent radii of Al, Ga, In, Ti, and Bi have been proposed.

The Schomaker-Stevenson revision of the additivity rule for single covalent bond distances has been used in conjunction with the revised radii. Agreement with experiment is in general better with the revised radii than with the former radii and additivity.

The principle of ionic bond character in addition to that present in a normal covalent bond has been applied to the observed structures of numerous molecules. It leads to a method of interpretation which is at least as consistent as the theory of multiple bond formation.

The revision of the additivity rule has been extended to double bonds. An encouraging beginning along these lines has been made, but additional experimental data are needed for clarification.

Design, Synthesis, and Study of Novel Platforms for Iron-N2 Chemistry and Photoinduced, Copper-mediated C-N Bond Formation

Several new ligand platforms designed to support iron dinitrogen chemistry have been developed. First, we report Fe complexes of a tris(phosphino)alkyl (CP^iPr₃) ligand featuring an axial carbon donor intended to conceptually model the interstitial carbide atom of the nitrogenase iron-molybdenum cofactor (FeMoco). It is established that in this scaffold, the iron center binds dinitrogen trans to the C_alkyl anchor in three structurally characterized oxidation states. Fe-C_alkyl lengthening is observed upon reduction, reflective of significant ionic character in the Fe-C_alkyl interaction. The anionic (CP^iPr₃)FeN₂^- species can be functionalized by a silyl electrophile to generate (CP^iPr₃)Fe-N₂SiR₃. This species also functions as a modest catalyst for the reduction of N₂ to NH₃. Next, we introduce a new binucleating ligand scaffold that supports an Fe(μ-SAr)Fe diiron subunit that coordinates dinitrogen (N₂-Fe(μ-SAr)Fe-N₂) across at least three oxidation states (Fe^IIFe^II, Fe^IIFe^I, and Fe^IFe^I). Despite the sulfur-rich coordination environment of iron in FeMoco, synthetic examples of transition metal model complexes that bind N₂ and also feature sulfur donor ligands remain scarce; these complexes thus represent an unusual series of low-valent diiron complexes featuring thiolate and dinitrogen ligands. The (N₂-Fe(μ-SAr)Fe-N₂) system undergoes reduction of the bound N₂ to produce NH₃ (~50% yield) and can efficiently catalyze the disproportionation of N₂H₄ to NH₃ and N₂. The present scaffold also supports dinitrogen binding concomitant with hydride as a co-ligand. Next, inspired by the importance of secondary-sphere interactions in many metalloenzymes, we present complexes of iron in two new ligand scaffolds ([SiP^NMe₃] and [SiP^iPr₂P^NMe]) that incorporate hydrogen-bond acceptors (tertiary amines) which engage in interactions with nitrogenous substrates bound to the iron center (NH₃ and N₂H₄). Cation binding is also facilitated in anionic Fe(0)-N₂ complexes. While Fe-N₂ complexes of a related ligand ([SiP^iPr₃]) lacking hydrogen-bond acceptors produce a substantial amount of ammonia when treated with acid and reductant, the presence of the pendant amines instead facilitates the formation of metal hydride species.

Additionally, we present the development and mechanistic study of copper-mediated and copper-catalyzed photoinduced C-N bond forming reactions. Irradiation of a copper-amido complex, ((m-tol)₃P)₂Cu(carbazolide), in the presence of aryl halides furnishes N-phenylcarbazole under mild conditions. The mechanism likely proceeds via single-electron transfer from an excited state of the copper complex to the aryl halide, generating an aryl radical. An array of experimental data are consistent with a radical intermediate, including a cyclization/stereochemical investigation and a reactivity study, providing the first substantial experimental support for the viability of a radical pathway for Ullmann C-N bond formation. The copper complex can also be used as a precatalyst for Ullmann C-N couplings. We also disclose further study of catalytic C_alkyl-N couplings using a CuI precatalyst, and discuss the likely role of [Cu(carbazolide)₂]^- and [Cu(carbazolide)₃]^- species as intermediates in these reactions.

Finally, we report a series of four-coordinate, pseudotetrahedral P₃Fe^II-X complexes supported by tris(phosphine)borate ([PhBP₃Fe^R]^-) and phosphiniminato X-type ligands (-N=PR'₃) that in combination tune the spin-crossover behavior of the system. Low-coordinate transition metal complexes such as these that undergo reversible spin-crossover remain rare, and the spin equilibria of these systems have been studied in detail by a suite of spectroscopic techniques.

A model for energy and morphology of crystalline grain boundaries with arbitrary geometric character

It has been well-established that interfaces in crystalline materials are key players in the mechanics of a variety of mesoscopic processes such as solidification, recrystallization, grain boundary migration, and severe plastic deformation. In particular, interfaces with complex morphologies have been observed to play a crucial role in many micromechanical phenomena such as grain boundary migration, stability, and twinning. Interfaces are a unique type of material defect in that they demonstrate a breadth of behavior and characteristics eluding simplified descriptions. Indeed, modeling the complex and diverse behavior of interfaces is still an active area of research, and to the author's knowledge there are as yet no predictive models for the energy and morphology of interfaces with arbitrary character. The aim of this thesis is to develop a novel model for interface energy and morphology that i) provides accurate results (especially regarding "energy cusp" locations) for interfaces with arbitrary character, ii) depends on a small set of material parameters, and iii) is fast enough to incorporate into large scale simulations.

In the first half of the work, a model for planar, immiscible grain boundary is formulated. By building on the assumption that anisotropic grain boundary energetics are dominated by geometry and crystallography, a construction on lattice density functions (referred to as "covariance") is introduced that provides a geometric measure of the order of an interface. Covariance forms the basis for a fully general model of the energy of a planar interface, and it is demonstrated by comparison with a wide selection of molecular dynamics energy data for FCC and BCC tilt and twist boundaries that the model accurately reproduces the energy landscape using only three material parameters. It is observed that the planar constraint on the model is, in some cases, over-restrictive; this motivates an extension of the model.

In the second half of the work, the theory of faceting in interfaces is developed and applied to the planar interface model for grain boundaries. Building on previous work in mathematics and materials science, an algorithm is formulated that returns the minimal possible energy attainable by relaxation and the corresponding relaxed morphology for a given planar energy model. It is shown that the relaxation significantly improves the energy results of the planar covariance model for FCC and BCC tilt and twist boundaries. The ability of the model to accurately predict faceting patterns is demonstrated by comparison to molecular dynamics energy data and experimental morphological observation for asymmetric tilt grain boundaries. It is also demonstrated that by varying the temperature in the planar covariance model, it is possible to reproduce a priori the experimentally observed effects of temperature on facet formation.

Finally, the range and scope of the covariance and relaxation models, having been demonstrated by means of extensive MD and experimental comparison, future applications and implementations of the model are explored.