Local character of magnetic coupling in ionic solids

Magnetic interactions in ionic solids are studied using parameter-free methods designed to provide accurate energy differences associated with quantum states defining the Heisenberg constant J. For a series of ionic solids including KNiF3, K2NiF4, KCuF3, K2CuF4, and high- Tc parent compound La2CuO4, the J experimental value is quantitatively reproduced. This result has fundamental implications because J values have been calculated from a finite cluster model whereas experiments refer to infinite solids. The present study permits us to firmly establish that in these wide-gap insulators, J is determined from strongly local electronic interactions involving two magnetic centers only thus providing an ab initio support to commonly used model Hamiltonians.

MgO/Ag(001) interface structure and STM images from first principles

The interface of MgO/Ag(001) has been studied with density functional theory applied to slabs. We have found that regular MgO films show a small adhesion to the silver substrate, the binding can be increased in off-stoichiometric regimes, either by the presence of O vacancies at the oxide film or by a small excess of O atoms at the interface between the ceramic to the metal. By means of theoretical methods, the scanning tunneling microscopy signatures of these films is also analyzed in some detail. For defect free deposits containing 1 or 2 ML and at low voltages, tunnelling takes place from the surface Ag substrate, and at large positive voltages Mg atoms are imaged. If defects, oxygen vacancies, are present on the surface of the oxide they introduce much easier channels for tunnelling resulting in big protrusions and controlling the shape of the image, the extra O stored at the interface can also be detected for very thin films.

Shifts of electronic KX-rays in muonic atoms

Energies of electronic K X-rays in muonic atoms were calculated for muons in various outer orbitals and for different numbers of electrons. Energy shifts were obtained with respect to characteristic X-rays belonging to nuclear charge (Z - 1) and their possible observation is discussed. The shifts in muonic Sn as an example amount to 19, 37, and 59 eV for the muon in 5g, 6h, and 7i states respectively. However, shifts due to the number of electrons present and the electron vacancy distribution in the L-shell are significantly larger. Accurate measurements of the K X-ray energies would therefore enable us to learn more about the electronic structure during the muonic cascade.

Total differential scattering cross section of {Ar^+}-Ar at 15 to 400 keV

We report on the measurement of the total differential scattering cross section of {Ar^+}-Ar at laboratory energies between 15 and 400 keV. Using an ab initio relativistic molecular program which calculates the interatomic potential energy curve with high accuracy, we are able to reproduce the detailed structure found in the experiment.

Volume isotope shifts in low lying transitions of Au I

Electronic factors in the volume isotope shift have been calculated in an ab initio way with the relativistic Dirac-Fock method for a number of different optical single/and two-photon transitions in Au I. The agreement with a semi-empirical method is within 10% for the resonance transition. For this one and a few other transitions the effect of core excitation has been analyzed with the Multi-configuration Dirac-Fock method as well, and it was found to reduce the electronic factor in the order of 5 %.

Femtosecond real-time probing of reactions

Femtosecond reaction dynamics of OClO in a supersonic molecular beam are reported. The system is excited to the A^2A_2 state with a femtosecond pulse, covering a range of excitation in the symmetric stretch between v_1 = 17 to v_1 = 11 (308-352 nm). A time-delayed femtosecond probe pulse ionizes the OClO, and OClO^+ is detected. This ion has not been observed in previous experiments because of its ultrafast fragmentation. Transients are reported for the mass of the parent OClO as well as the mass of the ClO. Apparent biexponential decays are observed and related to the fragmentation dynamics: OClO+hv \rightarrow (OClO)^{(++)*} \rightarrow ClO+O \rightarrow Cl+O_2. Clusters of OClO with water (OClO)_n (H_2 0)_m with n from 1 to 3 and m from 0 to 3 are also observed. The dynamics of the fragmentation reveal the nuclear motions and the electronic coupling between surfaces. The time scale for bond breakage is in the range of 300-500 fs, depending on v_1; surface crossing to form new intermediates is a pathway for the two channels of fragmentation: ClO+O (primary) and Cl+O_2 (minor). Comparisons with results of ab initio calculations are made.

Theoretical study of magnetism, structure and chemical order in transition-metal alloy clusters

Research on transition-metal nanoalloy clusters composed of a few atoms is fascinating by their unusual properties due to the interplay among the structure, chemical order and magnetism. Such nanoalloy clusters, can be used to construct nanometer devices for technological applications by manipulating their remarkable magnetic, chemical and optical properties. Determining the nanoscopic features exhibited by the magnetic alloy clusters signifies the need for a systematic global and local exploration of their potential-energy surface in order to identify all the relevant energetically low-lying magnetic isomers. In this thesis the sampling of the potential-energy surface has been performed by employing the state-of-the-art spin-polarized density-functional theory in combination with graph theory and the basin-hopping global optimization techniques. This combination is vital for a quantitative analysis of the quantum mechanical energetics. The first approach, i.e., spin-polarized density-functional theory together with the graph theory method, is applied to study the Fe$_m$Rh$_n$ and Co$_m$Pd$_n$ clusters having $N = m+n \leq 8$ atoms. We carried out a thorough and systematic sampling of the potential-energy surface by taking into account all possible initial cluster topologies, all different distributions of the two kinds of atoms within the cluster, the entire concentration range between the pure limits, and different initial magnetic configurations such as ferro- and anti-ferromagnetic coupling. The remarkable magnetic properties shown by FeRh and CoPd nanoclusters are attributed to the extremely reduced coordination number together with the charge transfer from 3$d$ to 4$d$ elements. The second approach, i.e., spin-polarized density-functional theory together with the basin-hopping method is applied to study the small Fe$_6$, Fe$_3$Rh$_3$ and Rh$_6$ and the larger Fe$_{13}$, Fe$_6$Rh$_7$ and Rh$_{13}$ clusters as illustrative benchmark systems. This method is able to identify the true ground-state structures of Fe$_6$ and Fe$_3$Rh$_3$ which were not obtained by using the first approach. However, both approaches predict a similar cluster for the ground-state of Rh$_6$. Moreover, the computational time taken by this approach is found to be significantly lower than the first approach. The ground-state structure of Fe$_{13}$ cluster is found to be an icosahedral structure, whereas Rh$_{13}$ and Fe$_6$Rh$_7$ isomers relax into cage-like and layered-like structures, respectively. All the clusters display a remarkable variety of structural and magnetic behaviors. It is observed that the isomers having similar shape with small distortion with respect to each other can exhibit quite different magnetic moments. This has been interpreted as a probable artifact of spin-rotational symmetry breaking introduced by the spin-polarized GGA. The possibility of combining the spin-polarized density-functional theory with some other global optimization techniques such as minima-hopping method could be the next step in this direction. This combination is expected to be an ideal sampling approach having the advantage of avoiding efficiently the search over irrelevant regions of the potential energy surface.

Many-body theory of laser-induced ultrafast demagnetization and angular momentum transfer in ferromagnetic transition metals

An electronic theory is developed, which describes the ultrafast demagnetization in itinerant ferromagnets following the absorption of a femtosecond laser pulse. The present work intends to elucidate the microscopic physics of this ultrafast phenomenon by identifying its fundamental mechanisms. In particular, it aims to reveal the nature of the involved spin excitations and angular-momentum transfer between spin and lattice, which are still subjects of intensive debate. In the first preliminary part of the thesis the initial stage of the laser-induced demagnetization process is considered. In this stage the electronic system is highly excited by spin-conserving elementary excitations involved in the laser-pulse absorption, while the spin or magnon degrees of freedom remain very weakly excited. The role of electron-hole excitations on the stability of the magnetic order of one- and two-dimensional 3d transition metals (TMs) is investigated by using ab initio density-functional theory. The results show that the local magnetic moments are remarkably stable even at very high levels of local energy density and, therefore, indicate that these moments preserve their identity throughout the entire demagnetization process. In the second main part of the thesis a many-body theory is proposed, which takes into account these local magnetic moments and the local character of the involved spin excitations such as spin fluctuations from the very beginning. In this approach the relevant valence 3d and 4p electrons are described in terms of a multiband model Hamiltonian which includes Coulomb interactions, interatomic hybridizations, spin-orbit interactions, as well as the coupling to the time-dependent laser field on the same footing. An exact numerical time evolution is performed for small ferromagnetic TM clusters. The dynamical simulations show that after ultra-short laser pulse absorption the magnetization of these clusters decreases on a time scale of hundred femtoseconds. In particular, the results reproduce the experimentally observed laser-induced demagnetization in ferromagnets and demonstrate that this effect can be explained in terms of the following purely electronic non-adiabatic mechanism: First, on a time scale of 10–100 fs after laser excitation the spin-orbit coupling yields local angular-momentum transfer between the spins and the electron orbits, while subsequently the orbital angular momentum is very rapidly quenched in the lattice on the time scale of one femtosecond due to interatomic electron hoppings. In combination, these two processes result in a demagnetization within hundred or a few hundred femtoseconds after laser-pulse absorption.

C-H···O H-bonded complexes: how does basis set superposition error change their potential-energy surfaces?

Geometries, vibrational frequencies, and interaction energies of the CNH⋯O3 and HCCH⋯O3 complexes are calculated in a counterpoise-corrected (CP-corrected) potential-energy surface (PES) that corrects for the basis set superposition error (BSSE). Ab initio calculations are performed at the Hartree-Fock (HF) and second-order Møller-Plesset (MP2) levels, using the 6-31G(d,p) and D95++(d,p) basis sets. Interaction energies are presented including corrections for zero-point vibrational energy (ZPVE) and thermal correction to enthalpy at 298 K. The CP-corrected and conventional PES are compared; the unconnected PES obtained using the larger basis set including diffuse functions exhibits a double well shape, whereas use of the 6-31G(d,p) basis set leads to a flat single-well profile. The CP-corrected PES has always a multiple-well shape. In particular, it is shown that the CP-corrected PES using the smaller basis set is qualitatively analogous to that obtained with the larger basis sets, so the CP method becomes useful to correctly describe large systems, where the use of small basis sets may be necessary

Analysis of chemical bonding and aromaticity from electronic delocalization descriptors

Interactions between electrons determine the structure and properties of matter from molecules to solids. Therefore, the understanding of the electronic structure of molecules will enable us to extract relevant chemical information. In the first part of this thesis, we focus our attention on the analysis of chemical bonding by means of the Electron Localization Function (ELF) and the Domain-Averaged Fermi Hole analysis (DAFH). In the second part, we assess the performance of some indicators of aromaticity by analyzing their advantages and drawbacks. We propose a series of tests based on well-known aromaticity trends that can be applied to evaluate the aromaticity of current and future indicators of aromaticity in both organic and inorganic species. Moreover, we investigate the nature of electron delocalization in both aromatic and antiaromatic systems in the light of Hückel’s (4n + 2) rule. Finally, we analyze the phenomenon of multiple aromaticity in all-metal clusters.

A potential energy surface for the ground state of acetylene, C2H2

A potential energy function has been derived for the ground state surface of C2H2 as a many-body expansion. The 2- and 3-body terms have been obtained by preliminary investigation of the ground state surfaces of CH2( 3B1) and C2H( 2Σ+). A 4-body term has been derived which reproduces the energy, geometry and harmonic force field of C2H2. The potential has a secondary minimum corresponding to the vinylidene structure and the geometry and energy of this are in close agreement with predictions from ab initio calculations. The saddle point for the HCCH-H2CC rearrangement is predicted to lie 2•530 eV above the acetylene minimum.

A potential energy surface for the ground state of formaldehyde, H2CO(1A1)

A model potential energy function for the ground state of H2CO has been derived which covers the whole space of the six internal coordinates. This potential reproduces the experimental energy, geometry and quadratic force field of formaldehyde, and dissociates correctly to all possible atom, diatom and triatom fragments. Thus there are good reasons for believing it to be close to the true potential energy surface except in regions where both hydrogen atoms are close to the oxygen. It leads to the prediction that there should be a metastable singlet hydroxycarbene HCOH which has a planar trans structure and an energy of 2•31 eV above that of equilibrium formaldehyde. The reaction path for dissociation into H2 + CO is predicted to pass through a low symmetry transition state with an activation energy of 4•8 eV. Both of these predictions are in good agreement with recently published ab initio calculations.

Analytical functions for the ground state potential surfaces of C3 (X 1 Σ g+) and HCN (X 1 Σ +)

Analytic functions have been obtained to represent the potential energy surfaces of C3 and HCN in their ground electronic states. These functions closely reproduce the available data on the energy, geometry, and force constants in all stable conformations, as well as data on the various dissociation products, and ab initio calculations of the energy at other conformations. The form of the resulting surfaces are portrayed in various ways and discussed briefly.

The general harmonic force field of methyl fluoride

The complete general harmonic force field of methyl flouride was recalculated using the most recent literature frequency, Coriolis ζ, and centrifugal distortion data for 12CH3F, 13CH3F, 12CD3F, 12CHD2F and 12CH2DF. The anharmonic corrections applied to the observed frequency data and the adopted molecular geometry are considered to be more realistic than those used hitherto. There is excellent overall agreement between the fitted force constants and the highest quality ab initio force field currently available.