Relativistische vier und zwei Spinor Finite Elemente Berechnungen zweiatomiger Moleküle

Der Vielelektronen Aspekt wird in einteilchenartigen Formulierungen berücksichtigt, entweder in Hartree-Fock Näherung oder unter dem Einschluß der Elektron-Elektron Korrelationen durch die Dichtefunktional Theorie. Da die Physik elektronischer Systeme (Atome, Moleküle, Cluster, Kondensierte Materie, Plasmen) relativistisch ist, habe ich von Anfang an die relativistische 4 Spinor Dirac Theorie eingesetzt, in jüngster Zeit aber, und das wird der hauptfortschritt in den relativistischen Beschreibung durch meine Promotionsarbeit werden, eine ebenfalls voll relativistische, auf dem sogenannten Minimax Prinzip beruhende 2-Spinor Theorie umgesetzt. Im folgenden ist eine kurze Beschreibung meiner Dissertation: Ein wesentlicher Effizienzgewinn in der relativistischen 4-Spinor Dirac Rechnungen konnte durch neuartige singuläre Koordinatentransformationen erreicht werden, so daß sich auch noch für das superschwere Th2 179+ hächste Lösungsgenauigkeiten mit moderatem Computer Aufwand ergaben, und zu zwei weiteren interessanten Veröffentlichungen führten (Publikationsliste). Trotz der damit bereits ermöglichten sehr viel effizienteren relativistischen Berechnung von Molekülen und Clustern blieben diese Rechnungen Größenordnungen aufwendiger als entsprechende nicht-relativistische. Diese behandeln das tatsächliche (relativitische) Verhalten elektronischer Systeme nur näherungsweise richtig, um so besser jedoch, je leichter die beteiligten Atome sind (kleine Kernladungszahl Z). Deshalb habe ich nach einem neuen Formalismus gesucht, der dem möglichst gut Rechnung trägt und trotzdem die Physik richtig relativistisch beschreibt. Dies gelingt durch ein 2-Spinor basierendes Minimax Prinzip: Systeme mit leichten Atomen sind voll relativistisch nunmehr nahezu ähnlich effizient beschrieben wie nicht-relativistisch, was natürlich große Hoffnungen für genaue (d.h. relativistische) Berechnungen weckt. Es ergab sich eine erste grundlegende Veröffentlichung (Publikationsliste). Die Genauigkeit in stark relativistischen Systemen wie Th2 179+ ist ähnlich oder leicht besser als in 4-Spinor Dirac-Formulierung. Die Vorteile der neuen Formulierung gehen aber entscheidend weiter: A. Die neue Minimax Formulierung der Dirac-Gl. ist frei von spuriosen Zuständen und hat keine positronischen Kontaminationen. B. Der Aufwand ist weit reduziert, da nur ein 1/3 der Matrix Elemente gegenüber 4-Spinor noch zu berechnen ist, und alle Matrixdimensionen Faktor 2 kleiner sind. C. Numerisch verhält sich die neue Formulierung ähnlilch gut wie die nichtrelativistische Schrödinger Gleichung (Obwohl es eine exakte Formulierung und keine Näherung der Dirac-Gl. ist), und hat damit bessere Konvergenzeigenschaften als 4-Spinor. Insbesondere die Fehlerwichtung (singulärer und glatter Anteil) ist in 2-Spinor anders, und diese zeigt die guten Extrapolationseigenschaften wie bei der nichtrelativistischen Schrödinger Gleichung. Die Ausweitung des Anwendungsbereichs von (relativistischen) 2-Spinor ist bereits in FEM Dirac-Fock-Slater, mit zwei Beispielen CO und N2, erfolgreich gemacht. Weitere Erweiterungen sind nahezu möglich. Siehe Minmax LCAO Nährung.

A computer-algebraic approach to the study of the symmetry properties of molecules and clusters

This thesis work is dedicated to use the computer-algebraic approach for dealing with the group symmetries and studying the symmetry properties of molecules and clusters. The Maple package Bethe, created to extract and manipulate the group-theoretical data and to simplify some of the symmetry applications, is introduced. First of all the advantages of using Bethe to generate the group theoretical data are demonstrated. In the current version, the data of 72 frequently applied point groups can be used, together with the data for all of the corresponding double groups. The emphasize of this work is placed to the applications of this package in physics of molecules and clusters. Apart from the analysis of the spectral activity of molecules with point-group symmetry, it is demonstrated how Bethe can be used to understand the field splitting in crystals or to construct the corresponding wave functions. Several examples are worked out to display (some of) the present features of the Bethe program. While we cannot show all the details explicitly, these examples certainly demonstrate the great potential in applying computer algebraic techniques to study the symmetry properties of molecules and clusters. A special attention is placed in this thesis work on the flexibility of the Bethe package, which makes it possible to implement another applications of symmetry. This implementation is very reasonable, because some of the most complicated steps of the possible future applications are already realized within the Bethe. For instance, the vibrational coordinates in terms of the internal displacement vectors for the Wilson's method and the same coordinates in terms of cartesian displacement vectors as well as the Clebsch-Gordan coefficients for the Jahn-Teller problem are generated in the present version of the program. For the Jahn-Teller problem, moreover, use of the computer-algebraic tool seems to be even inevitable, because this problem demands an analytical access to the adiabatic potential and, therefore, can not be realized by the numerical algorithm. However, the ability of the Bethe package is not exhausted by applications, mentioned in this thesis work. There are various directions in which the Bethe program could be developed in the future. Apart from (i) studying of the magnetic properties of materials and (ii) optical transitions, interest can be pointed out for (iii) the vibronic spectroscopy, and many others. Implementation of these applications into the package can make Bethe a much more powerful tool.

Theoretical study of the structural dependence of nuclear quadrupole frequencies in high-T_c superconductors

The remarkable difference between the nuclear quadrupole frequencies v_Q of Cu(1) and Cu(2) in YBa_2Cu_3O_6 and YBa_2Cu_3O_7 is analyzed. We calculate the ionic contribution to the electric field gradients and estimate, by using experimental results for Cu_2O and La_2CuO_4, the contribution of the d valence electrons. Thus, we determine v_Q1, v_Q2, and the asymmetry parameter η for YBa_2Cu_3O_6 and YBa_2Cu_3O_7. The number of holes in dthe Cu-O planes and chains is found to be important for the different behavior of v_Q1 and v_Q2.

Delocalization of a hole in van der Waals clusters: Ionization potential of rare-gas and small Hg_n clusters

The size dependence of the ionization potential I_p(n) of van der Waals (vdW) bound clusters has been calculated by using a model Hamiltonian, which includes electron hopping, vdW interactions, and charge-dipole interactions. The charge-density and dipole-density distributions for both neutral and ionized n-atom clusters are determined self-consistently. The competition between the polarization energy of the neutral atoms surrounding a partially localized hole and the tendency toward hole delocalization in the ionized clusters is found to dominate the size dependence of I_p(n). To test our theory, we culculate I_p(Xe_n) and I_p(Kr_n) for n \le 300. Good quantitative agreement with experiment is obtained. The theory is also applied to calculate I_p(Hg_n). Comparison with experiments suggests that in Hg_n^+ clusters with n \le 20 the positive charge is mainly distributed within a trimer which is situated at the center of the cluster and which polarizes the n - 3 surrounding neutral atoms.

Theory for Optical Absorption in Small Clusters: Dependence on Atomic Structure and Cluster Size

We present a theory which permits for the first time a detailed analysis of the dependence of the absorption spectrum on atomic structure and cluster size. Thus, we determine the development of the collective excitations in small clusters and show that their broadening depends sensitively on the tomic structure, in particular at the surface. Results for Hg_n^+ clusters show that the plasmon energy is close to its jellium value in the case of spherical-like structures, but is in general between w_p/ \wurzel{3} and w_p/ \wurzel{2} for compact clusters. A particular success of our theory is the identification of the excitations contributing to the absorption peaks.

Unrestricted hartree-fock calculation of the ionization potential of small Hg_n clusters

The ionization potential of small Hg_n clusters has been calculated. For the first time good agreement with experimental results has been obtained. It is shown that interatomic Coulomb interactions are important. The energy of Hg_n^+ is calculated using the unrestricted inhomogeneous Hartree-Fock approximation. As a consequence of a change in the charge distribution in Hg_n^+ , we obtain an abrupt change in the slope of the ionization potential at the critical cluster size n_cr ~ 14. The presented results are expected to be valid for covalent clusters in between ionized van der Waals clusters and metallic clusters.

Calculation of the electronic properties of neutral and ionized divalent-metal clusters

The electronic properties of neutral and ionized divalent-metal clusters have been studied using a microscopic theory, which takes into account the interplay between van der Waals (vdW) and covalent bonding in the neutral clusters, and the competition between hole delocalization and polarization energy in the ionized clusters. By calculating the ground-state energies of neutral and ionized. Hg_n clusters, we determine the size dependence of the bond character and the ionization potential I_p(n). For neutral Hg_n clusters we obtain a transition from van del Waals to covalent behaviour at the critical size n_c ~ 10-20 atoms. Results for I_p(Hg_n) with n \le 20 are in good agreement with experiments, and suggest that small Hg_n^+ clusters can be viewed as consisting of a positive trimer core Hg_3^+ surrounded by n - 3 polarized neutral atoms.

On the localization of electrons and holes in Hg_n and rare-gas clusters

We have used a microscopic theory to study the size dependence of the degree of localization of the valence electrons and holes in neutral an ionized rare-gas and Hg_n clusters. We discuss under which circumstances localization of the electrons and holes is favoured. We have calculated the ionization potential of Xe_n, Kr_n and small Hg_n clusters. Good agreement with experiments is obtained. We have also determined the dependence of the ionization potential on cluster structure.

On the transition from localized to delocalized electronic states in divalent-metal clusters

The transition from van der Waals to covalent bonding, which is expected to occur in divalent-metal clusters with increasing cluster size, is discussed. We propose a model which takes into account, within the same electronic theory, the three main competing contributions, namely the kinetic energy of the electrons, the Coulomb interactions between electrons, and the s \gdw p intraatomic transitions responsible for van der Waals like bonding. The model is solved by taking into account electron correlations using a generalized Gutzwiller approximation (slave boson method). The occurrence of electron localization is studied as a function of the interaction parameters and cluster size.

Theory for the optical properties of small Hg_n clusters

The static and dynamical polarizabilities of the Hg-dimer are calculated by using a Hubbard Hamiltonian to describe the electronic structure. The Hamiltonian is diagonalized exactly within a subspace of second-quantized electronic states from which only multiply ionized atomic configurations have been excluded. With this approximation we can describe the most important electronic transitions including the effect of charge fluctuations. We analyze the polarizability as a function of the intraatomic Coulomb interaction which represents the repulsion between electrons. We obtain that this interaction results in strong electronic correlations in the excited states and increases the first excitation energy of the dimer by 0.8 eV in comparison to a calculation which neglects correlations, resulting in a better agreement with the experiment.

Theory for the change of bond character in divalent-metal clusters

To determine the size dependence of the bonding in divalent-metal clusters we use a many-electron Hamiltonian describing the interplay between van der Waals (vdW) and covalent interactions. Using a saddle-point slave-boson method and taking into account the size-dependent screening of charge fluctuations, we obtain for Hg_n a sharp transition from vdW to covalent bonding for increasing n. We show also, by solving the model Hamiltonian exactly, that for divalent metals vdW and covalent bonding coexist already in the dimers.

Slave-boson approach to the size-dependent transition from van der Waals to covalent bonding in Hg_n clusters

We use a microscopic theory to describe the dynamics of the valence electrons in divalent-metal clusters. The theory is based on a many-body model Harniltonian H which takes into account, on the same electronic level, the van der Waals and the covalent bonding. In order to study the ground-state properties of H we have developed an extended slave-boson method. We have studied the bonding character and the degree of electronic delocalization in Hg_n clusters as a function of cluster size. Results show that, for increasing cluster size, an abrupt change occurs in the bond character from van der Waals to covalent bonding at a critical cluster size n_c ~ 10-20. This change also involves a transition from localized to delocalized valence electrons, as a consequence of the competition between both bonding mechanisms.

Electronic structure calculations of small AI_n (n = 2 - 8) clusters

The electronic states of small AI_n (n = 2 - 8) clusters have been calculated with a relativistic ab-initio MOLCAO Dirac-Fock-Slater method using numerical atomic DFS wave-functions. The excitation energies were obtained from a ground state calculation of neutral clusters, and in addition from negative clusters charged by half an electron in order to account for part of the relaxation. These energies are compared with experimental photoelectron spectra.

Nuclear charge radii from X-ray transitions in muonic atoms of carbon, nitrogen and oxygen

Energies of muonic X-rays of the K-series of carbon, nitrogen and oxygen have been measured with an accuracy of about 15 eV. Root mean square radii of the nuclear charge distributions were deduced. The results 2.49±0.05 fm for carbon, 2.55 ±0.03 fm for nitrogen and 2.71 ±0.02 fm for oxygen are in good agreement at comparable accuracy with recent electron scattering data.

Elastic nuclear scattering at intermediate and relativistic energies

The classical scattering cross section of two colliding nuclei at intermediate and relativistic energies is reevaluated. The influence of retardation and magnetic field effects is taken into account. Corrections due to electron screening as well as due to attractive nuclear forces are discussed. This paper represents an addendum to [l].