The complete and accurate thermodynamic and kinetic description of any systemis crucialfor understanding and predicting its properties. A particular interest is in systemsthat are used for some practical applications and have to be constantly improved usingmodification of their composition and structure. This task can be quite accuratelysolved at a fundamental level by density functional theory methods. Thesemethods areapplied to two practically important systems Fe-Cr and TiC-ZrC.The elastic properties of pure iron and substitutionally disordered Fe-Cr alloy are investigatedas a function of temperature and concentration using first-principles electronicstructurecalculations by the exact muffin-tin orbitals method. The temperature effectson the elastic properties are included via the electronic, magnetic, and lattice expansioncontributions. It is shown that the degree of magnetic order in both pure iron andFe90Cr10 alloy mainly determines the dramatic change of the elastic anisotropy of thesematerials at elevated temperatures. A peculiarity in the concentration dependence ofthe elastic constants in Fe-rich alloys is demonstrated and related to a change in theFermi surface topology.A thermodynamic model for the magnetic alloys is developed from first principles andapplied to the calculation of bcc Fe-Cr phase diagram. Various contributions to the freeenergy (magnetic, electronic, and phonon) are estimated and included in the model. Inparticular, it is found that magnetic short range order effects are important just abovethe Curie temperature. The model is applied for calculating phase equilibria in disorderedbcc Fe-Cr alloys. Model calculations reproduce a feature known as a Nishizawahorn for the Fe-rich high-temperature part of the phase diagram.The investigation of the TiC-ZrC system includes a detailed study of the defect formationenergies and migration barriers of point defects and defect complexes involvedin the diffusion process. It is found, using ab initio atomistic simulations of vacancymediateddiffusion processes in TiC and ZrC, that a special self-diffusion mechanism isoperative for metal atom diffusion in sub-stoichiometric carbides. It involves a noveltype of a stable point defect, a metal vacancy ”dressed” in a shell of carbon vacancies.It is shown that this vacancy cluster is strongly bound and can propagate through thelattice without dissociating. / <p>QC 20120604</p> / HERO-M

Ab-initio results of the core-level shift and the distribution about the average for the 3d5/2 electrons of Ag, Pd and 2p3/2 of Cu are presented for the face-centered-cubic AgPd and CuPd random alloys. The complete screening model, which includes both initial and final states effects in the same scheme, has been used in the investigations. The alloys have been modeled with a supercell containing 256 atoms. Density-functional theory calculations are carried out using the locally self consistent Green's function approach. Results from the calculations clearly shows that the core-level shift distributions characteristic is Gaussian, but the components reveals a substantial difference in the FWHM (Full-Width at Half-Maximum). Comparison between the experimental and the calculated broadening shows a remarkable agreement.

The vertical electronic excitation energies for the narrow-bandgap polymers LBPF, EP37 and EP62 have been calculated using Density Functional Theory (DFT). Also the vertical excitation energies for the acceptor unit of LBPF have been calculated using the Hartree-Fock (HF), DFT and Coupled Cluster (CC) methods. The calculations cover the visible and infrared wave length region and two strong transitions are obtained, one corresponding to the pi to pi* transition and one corresponding to the pi to Acceptor transition. The excitation energies obtained from DFT are below the corresponding experimental results and attempts have therefore been made to perform bench-marking calculations using a hierarchy of CC methods.

The absorption spectra of a recently designed low-bandgap conjugated polymer has been studied using the semi-empirical method ZINDO and TDDFT/B3LYP/6-31G. The vertical excitation energies have been calculated for monomer up to hexamer. Two main absorption peaks can be seen, the one largest in wavelength corresponding to a HOMO to LUMO transition, and one involving higher order excitations. TDDFT results are red-shifted compared to the ZINDO results. Comparison with experiment yields that short conjugation lengths are dominating. This is possibly due to steric interactions between polymer chains, breaking the conjugation length. Such effects are also studied.

This thesis describes an investigation of the molecular mechanism underlying solar conversion processes that occur in Type I photosynthetic reaction centers, in which P700 plays a central role. Static Fourier transform infrared (FTIR) difference spectroscopy (DS) was used to probe the electronic and structural organization of P700 and P700+. In combination with isotope labeling and site directed mutagenesis we have investigated how protein interactions such as histidine ligation and hydrogen bonding modulate this organization. Comparison of (P700+-P700) FTIR difference spectra (DS) obtained using wild type and mutant PS I led us to suggest that the 131 keto carbonyl group of PA is essentially free from hydrogen bonding in the ground state. Upon cation formation, this hydrogen bonding becomes stronger, probably because of a cation induced reorientation of the hydroxyl group of a nearby threonine residue. We also tentatively suggested that a difference band at 1639(-)/1660(+) cm-1 in (P700+-P700) FTIR DS might be due to a C=C mode of the imidazole side chain of the ligating histidine residues. Most of this thesis is geared towards investigating the validity of this interpretation. (P700+-P700) FTIR DS obtained using mutant PS I particles in which hydrogen bonding to P700 is altered can be reconciled within the context of our new interpretation. (P700+-P700) FTIR DS obtained using uniformly 2H, 15N, and 13C labeled PS I particles also support our new interpretation, and indicate that the difference band at 1639(-)/ 1660(+) cm-1 cannot be associated with a strongly hydrogen bonded keto carbonyl group of PA. To investigate if the imidazole side-chain of ligating histidine residues could contribute to bands in (P700+-P700) FTIR DS vibrational mode frequencies and intensities for several protonation forms of 4-methylimidazole were calculated. The calculations suggest that the 1639(-)/1660(+) cm-1 band in (P700+-P700) FTIR DS may not be due to a C=C mode of the imidazole side chain of the ligating histidine residues. Thus we have produced data that suggests neither of the proposed interpretations alone can adequately explain the origin of the 1639(-)/1660(+) cm-1 difference band in (P700+-P700) FTIR DS. The origin of the 1639(-)/1660(+) cm-1 difference band in (P700+-P700) FTIR DS is therefore still an open question.

In this work the electronic and structural properties of anatase TiO2 (100) surface and gold adsorption have been investigated by using the first-principles calculations based on density functional theory (DFT). TiO2 is a wide band-gap material and to this effects it finds numerous applications in technology such as, cleaning of water, self-cleaning, coating, solar cells and so on. Primarily, the relation between the surface energy of the anatase (100)-1x1 phase and the TiO2-layers is examined. After an appropriate atomic layer has been chosen according to the stationary state of the TiO2 slab, the adsorption behavior of the Au atom and in the different combinations are searched for both the surface and the surface which is supported by a single Au atom/atoms. It has been observed that a single Au atom tends to adsorb to the surface which has an impurity of Au atom or atoms. Although, the high metal concentration on the surface have increased the strength of the adsorption, it is indicated that the system gains a metallic property which is believed to cause problems in the applications. In addition, the gold clusters of the dimer (Au2) and the trimer (Au3) have been adsorbed on the surface and their behavior on the surface is investigate. It is observed that the interaction between Au atoms in the atomic cluster each other is stronger than that of gold clusters and the surface.

The correlation between NMR chemical shifts of interstitial atoms and electronic structures of boron- and carbon-centered hexazirconium halide clusters was investigated by density functional theory (DFT) calculation. The influences of bridging halide and terminal ligand variations on electronic structure were examined respectively. Inverse proportionality was found between the chemical shifts and the calculated energy gaps between two Kohn-Sham orbitals of t1u symmetry, which arose from the bonding and antibonding interaction between the zirconium cage bonding orbitals and the interstitial 2p orbitals. Chemical shielding properties of the interstitial atoms were calculated with Gauge Including Atomic Orbital (GIAO) method. Stepwise ligand substitution of terminal chlorides on [(Zr6CCl12)Cl6]4-cluster by tri(n-butyl)-phosphine oxide (Bu3PO) was conducted with the aid of TlPF6. Composition of the reaction mixtures was analyzed by use of both 13C and 31P NMR. A preliminary scheme for synthesis and separation of [(Zr6CCl12)Cl6-x(Bu3PO)x]x-4 (x = 3 ?? 5) mixture based on solubility difference was reevaluated. Three 1,10-phenanthroline based bidentate ligands, namely, 2,9-Bis(diphenyl-phosphinyl)-1,10-phenanthroline, 2,9-Bis(diethoxyphosphoryl)-1,10-phenanthroline, and 2,9-Bis(di-n-butoxyphosphoryl)-1,10-phenantholine, were synthesized for bridge-chelating the hexazirconium clusters. Coordination chemistry of these ligands with the [Zr6BCl12] and [Zr6CCl12] clusters was subject to preliminary investigation.

Metal hydrides are promising candidates for H2 storage, but high stability and poor kinetics are the important challenges which have to be solved for vehicular applications. Most of recent experimental reports for improving thermodynamics of metal hydrides have been focused on lowering reaction enthalpies of a metal hydride by mixing other compounds. However, finding out metal hydride mixtures satisfying favorable thermodynamics among a large number of possible metal hydride mixtures is inefficient and thus a systematic approach is required for an efficient and rigorous solution. Our approaches introduced in this thesis allow a systematic screening of promising metal hydrides or their mixtures from all possible metal hydrides and their mixtures. Our approaches basically suggest two directions for improving metal hydride thermodynamics. First, our calculations for examining the relation between the particle size of simple metal hydrides and thermodynamics of their decomposition reactions provide that the relation would depend on the total surface energy difference between a metal and its hydride form. It ultimately suggests that we will be able to screen metal hydride nanoparticles having favorable thermodynamics from all possible metal hydrides by examining the total surface differences. Second, more importantly, we suggest that our thermodynamic calculations combined with the grand canonical linear programming method and updated database efficiently and rigorously screen potential promising bulk metal hydrides and their mixtures from a large collection of possible combinations. The screened promising metal hydrides and their mixtures can release H2 via single step or multi step. Our additional free energy calculations for a few selected promising single step reactions and their metastable paths show that we can identify the most stable free energy paths for any selected reactant mixtures. In this thesis, we also demonstrate that a total free energy minimization method can predict the possible evolution of impurity other than H2 for several specified mixtures. However, it is not ready to predict reaction thermodynamics from a large number of compounds.

Theoretical strength of solids is defined as the ultimate strength beyond which plastic deformation, fracture, or decohesion would occur. Understanding the microscopic origin from quantum mechanics and thermoelastic formulation is of great importance to mechanical properties and engineering design of various solids. While quite a few theory models have been made in the past century by several generations of scientists, including Frankel and Born, a general and convincing framework has not been fully established. We study this issue from three respects: (1) Unify various elastic stability criteria for solids that determine an upper bound of theoretical strength; (2) with ab initio method, we test the elastic stability conditions of crystal Au. The phenomenon of bifurcation is observed: under hydrostatic expansion, the rhombohedral modulus reaches zero first of all; while under uniaxial tensile stress, the tetragonal shear modulus first reaches zero; (3) propose a nonlinear theoretical formulation of stability criterion. As an analytic method, this scheme is quite simple, in the mean time, it saves computation resource.

Plastic deformation in crystalline materials largely depends on the properties of dislocations, in particular their mobility. While continuum description of deformation of a crystalline metal can be made reasonably well by considering the elastic properties of dislocations and neglecting the core, crystallographic aspects of dislocation motion require precise understanding of the core effects. The concept of the generalized stacking fault (GSF) energy was proposed as means to describe this. GSF energy, a fundamental property of a given material, can be determined using first principles total energy calculations. In this thesis, we use GSF to understand some of the intriguing mechanical responses recently observed in some metallic systems. First, we examine the structures and stacking fault energies in Mg-Zn-Y alloy system. This system is unique in the sense that trace additions of Zn and or Y result in long period stacking sequences such as 6l and 14l, as reported in recent literature. Further, these alloys exhibit extraordinary mechanical properties. We attempt to rationalize these experimental observations through first principles calculations of energies of periodic structures with different stacking sequences and stacking faults. For pure Mg, we find that the 6-layer structure with the ABACAB stacking is most stable after the lowest energy hcp structure with ABAB stacking. Charge density analysis shows that the 2l and 6l structures are electronically similar, which might be a cause for better stability of 6l structure over a 4l sequence or other periodic structures. Addition of 2 atomic% Y leads to stabilization of the structure to 6l sequence whereas the addition of 2 atomic% Zn makes the 6l energetically comparable to that of the hcp. Stacking fault (SF) on the basal plane of 6l structure is higher in energy than that of the hcp 2l Mg, which further increases upon Y doping and decreases significantly with Zn doping. SF energy surface for the prismatic slip indicates dissociation of dislocations in alloys with a 6l structure. Thus, in an Mg-Zn-Y alloy, Y stabilizes the long periodicity, while Zn doping provides a synergistic effect in improving the mechanical properties alongwith strengthening due to long periodic phases. Our investigation of surface properties and magnetism in Ni revealed that, the universal binding energy relation (UBER) derived earlier to describe the cohesion between two rigid atomic planes, does not accurately capture the cohesive properties when the cleavage cracked surfaces are allowed to relax through atomic displacements. We find that two characteristic length-scales are involved in the cleavage of a crystal accompanied by structural relaxation at the cleaved surface. Based on that, we suggest a modified functional form of UBER that is analytical and at the same time accurately models the properties of relaxed surfaces upon cleavage. We demonstrate the generality as well as the validity of this modified UBER through first-principles density functional theory calculations of cleavage in fcc, bcc, and hcp metals, as well as covalently bonded materials. We also found that the cohesive law (stress-displacement relation) differs significantly in the case where cracked surfaces are allowed to relax, with lower peak stresses occuring at higher displacements. We have attempted understanding these ideas through images obtained from electronic densities and eigen states. Our work should be useful in providing inputs to multi-scale simulations of fracture in materials. The third phase of the work reports the stacking fault energy and twinning in Ni with a particular emphasis on the size effect. Experimental and computational research on Nan crystalline metals (mostly on Ni) indicates unique facets of dislocation activity (extended partial dislocations) and modes of deformation (twinning). In order to capture the intrinsic scaling eject in the nano-regime, it is imperative to account for the complex electronic structure of the metal in question. The stacking fault (SF) and twinning fault (TF) energies in nano-thin elm of Ni with 7, 13, 19, and 25 layers of (111) planes were determined using rest-principles density functional theory (DFT) total energy calculations. Generalized planar fault (GPF) energy curves of the nano-thin alms show higher extreme vis-a-vis the bulk, indicating that creation of SFs in nano-Ni is relatively difficult. In contrast, the ratios of energy barriers relevant to nucleation of dislocations and twinning support the observed enhanced tendency for extended partial dislocation formation and twinning in the nano-thin films in comparison with bulk. Our results should be useful in understanding deformation behavior of nano-structured Ni-based alloys used as advanced structural materials.