Systematic multi-configuration calculations on structures and properties of complex atoms

Die relativistische Multikonfigurations Dirac-Fock (MCDF) Methode ist gegenwärtig eines der am häufigsten benutzten Verfahren zur Berechnung der elektronischen Struktur und der Eigenschaften freier Atome. In diesem Verfahren werden die Wellenfunktionen ausgewählter atomarer Zustände als eine Linearkombination von sogenannten Konfigurationszuständen (CSF - Configuration State Functions) konstruiert, die in einem Teilraum des N-Elektronen Hilbert-Raumes eine (Vielteilchen-)Basis aufspannen. Die konkrete Konstruktion dieser Basis entscheidet letzlich über die Güte der Wellenfunktionen, die üblicherweise mit Hilfe einer Variation des Erwartungswertes zum no-pair Dirac-Coulomb Hamiltonoperators gewonnen werden. Mit Hilfe von MCDF Wellenfunktionen können die dominanten relativistischen und Korrelationseffekte in freien Atomen allgemein recht gut erfaßt und verstanden werden. Außer der instantanen Coulombabstoßung zwischen allen Elektronenpaaren werden dabei auch die relativistischen Korrekturen zur Elektron-Elektron Wechselwirkung, d.h. die magnetischen und Retardierungsbeiträge in der Wechselwirkung der Elektronen untereinander, die Ankopplung der Elektronen an das Strahlungsfeld sowie der Einfluß eines ausgedehnten Kernmodells erfaßt. Im Vergleich mit früheren MCDF Rechnungen werden in den in dieser Arbeit diskutierten Fallstudien Wellenfunktionsentwicklungen verwendet, die um 1-2 Größenordnungen aufwendiger sind und daher systematische Untersuchungen inzwischen auch an Atomen mit offenen d- und f-Schalen erlauben. Eine spontane Emission oder Absorption von Photonen kann bei freien Atomen theoretisch am einfachsten mit Hilfe von Übergangswahrscheinlichkeiten erfaßt werden. Solche Daten werden heute in vielen Forschungsbereichen benötigt, wobei neben den traditionellen Gebieten der Fusionsforschung und Astrophysik zunehmend auch neue Forschungsrichtungen (z.B. Nanostrukturforschung und Röntgenlithographie) zunehmend ins Blickfeld rücken. Um die Zuverlässigkeit unserer theoretischen Vorhersagen zu erhöhen, wurde in dieser Arbeit insbesondere die Relaxation der gebundenen Elektronendichte, die rechentechnisch einen deutlich größeren Aufwand erfordert, detailliert untersucht. Eine Berücksichtigung dieser Relaxationseffekte führt oftmals auch zu einer deutlich besseren Übereinstimmung mit experimentellen Werten, insbesondere für dn=1 Übergänge sowie für schwache und Interkombinationslinien, die innerhalb einer Hauptschale (dn=0) vorkommen. Unsere in den vergangenen Jahren verbesserten Rechnungen zu den Wellenfunktionen und Übergangswahrscheinlichkeiten zeigen deutlich den Fortschritt bei der Behandlung komplexer Atome. Gleichzeitig kann dieses neue Herangehen künftig aber auch auf (i) kompliziertere Schalensstrukturen, (ii) die Untersuchung von Zwei-Elektronen-ein-Photon (TEOP) Übergängen sowie (iii) auf eine Reihe weiterer atomarer Eigenschaften übertragen werden, die bekanntermaßen empflindlich von der Relaxation der Elektronendichte abhängen. Dies sind bspw. Augerzerfälle, die atomare Photoionisation oder auch strahlende und dielektronische Rekombinationsprozesse, die theoretisch bisher nur selten überhaupt in der Dirac-Fock Näherung betrachtet wurden.

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.

A semi-phenomenological approach to the structure and transport properties of macromolecules in solution

Motiviert durch die Lebenswissenschaften (Life sciences) haben sich Untersuchungen zur Dynamik von Makromolekülen in Lösungen in den vergangenen Jahren zu einem zukunftsweisenden Forschungsgebiet etabliert, dessen Anwendungen von der Biophysik über die physikalische Chemie bis hin zu den Materialwissenschaften reichen. Neben zahlreichen experimentellen Forschungsprogrammen zur räumlichen Struktur und den Transporteigenschaften grosser MolekÄule, wie sie heute praktisch an allen (Synchrotron-) Strahlungsquellen und den Laboren der Biophysik anzutreffen sind, werden gegenwärtig daher auch umfangreiche theoretische Anstrengungen unternommen, um das Diffusionsverhalten von Makromolekülen besser zu erklären. Um neue Wege für eine quantitative Vorhersagen des Translations- und Rotationsverhaltens grosser Moleküle zu erkunden, wurde in dieser Arbeit ein semiphänomenologischer Ansatz verfolgt. Dieser Ansatz erlaubte es, ausgehend von der Hamiltonschen Mechanik des Gesamtsystems 'Molekül + Lösung', eine Mastergleichung für die Phasenraumdichte der Makromoleküle herzuleiten, die den Einfluss der Lösung mittels effektiver Reibungstensoren erfasst. Im Rahmen dieses Ansatzes gelingt es z.B. (i) sowohl den Einfluss der Wechselwirkung zwischen den makromolekularen Gruppen (den sogenannten molekularen beads) und den Lösungsteilchen zu analysieren als auch (ii) die Diffusionseigen schaften für veschiedene thermodynamische Umgebungen zu untersuchen. Ferner gelang es auf der Basis dieser Näherung, die Rotationsbewegung von grossen Molekülen zu beschreiben, die einseitig auf einer Oberfläche festgeheftet sind. Im Vergleich zu den aufwendigen molekulardynamischen (MD) Simulationen grosser Moleküle zeichnet sich die hier dargestellte Methode vor allem durch ihren hohen `Effizienzgewinn' aus, der für komplexe Systeme leicht mehr als fünf Grössenordnungen betragen kann. Dieser Gewinn an Rechenzeit erlaubt bspw. Anwendungen, wie sie mit MD Simulationen wohl auch zukünftig nicht oder nur sehr zögerlich aufgegriffen werden können. Denkbare Anwendungsgebiete dieser Näherung betreffen dabei nicht nur dichte Lösungen, in denen auch die Wechselwirkungen der molekularen beads zu benachbarten Makromolekülen eine Rolle spielt, sondern auch Untersuchungen zu ionischen Flüssigkeiten oder zur Topologie grosser Moleküle.

Transport properties of one-dimensional, disordered two-band systems

We have studied the transport properties of disordered one-dimensional two-band systems. The model includes a narrow d band hybridised with an s band. The Landauer formula was used in the case of a very narrow d band or in the case of short chains. The results were compared with the localisation length of the wavefunctions calculated by the transfer matrix method, which allows the use of very lang chains, and an excellent agreement was obtained.

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.

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.

Chemical and physical properties of superheavy elements

A knowledge of the physical and chemical properties of superheavy elements is expected to be of great value for the detection of these elements, owing to the need for chemical separation in their isolation and identification. The methods for predicting their electronic structures, expected trends in their chemical and physical properties and the results of such predictions for the individual superheavy elements are reviewed. The periodic table is extended up to element 172.

Relativistic effects in physics and chemistry of element 105. [Part] I.

A detailed study of the electronic structure and bonding of the pentahalides of group 5 elements V, Nb, Ta, and element 105, hahnium (and Pa) has been carried out using relativistic molecular cluster Dirac-Slater discrete-variational method. A number of calculations have been performed for different geometries and molecular bond distances. The character of the bonding has been analyzed using the Mulliken population analysis of the molecular orbitals. It is shown that hahnium is a typical group 5 element. In a great number of properties it continues trends in the group. Some peculiarities in the electronic structure of HaCl_5 result from relativistic effects.

Relativistic effects in physics and chemistry of element 105. [Part] II.

Relativistic self-consistent charge Dirac-Slater discrete variational method calculations have been done for the series of molecules MBr_5, where M = Nb, Ta, Pa, and element 105, Ha. The electronic structure data show that the trends within the group 5 pentabromides resemble those for the corresponding pentaclorides with the latter being more ionic. Estimation of the volatility of group 5 bromides has been done on the basis of the molecular orbital calculations. According to the results of the theoretical interpretation HaBr_5 seems to be more volatile than NbBr_5 and TaBr_5.

Relativistic effects in physics and chemistry of element 105. [Part] III.

Electronic structures of MOCl_3 and MOBr_3 molecules, where M = V, Nb, Ta, Pa, and element 105, hahnium, have been calculated using the relativistic Dirac-Slater discrete variational method. The character of bonding has been analyzed using the Mulliken population analysis of the molecular orbitals. It was shown that hahnium oxytrihalides have similar properties to oxytrihalides of Nb and Ta and that hahnium has the highest tendency to form double bond with oxygen. Some peculiarities in the electronic structure of HaOCl_3 and HaOBr_3 result from relativistic effects. Volatilities of the oxytrihalides in comparison with the corresponding pentahalides were considered using results of the present calculations. Higher ionic character and lower covalency as well as the presence of dipole moments in MOX_3 (X = Cl, Br) molecules compared to analogous MX_5 ones are the factors contributing to their lower volatilities.

Theoretical study of the physicochemical properties of the light transactinides

A theoretical study of the physicochemical properties of elements 104, 105, and 106 and their compounds in the gas phase and aqueous solutions has been undertaken using relativistic atomic and molecular codes. Trends in properties such as bonding, ionization potentials, electron affinities, energies of electronic transitions, stabilities of oxidation states etc. have been defined within the corresponding chemical groups and within the transactinides. These trends are shown to be determined by increasing relativistic effects within the groups. The behaviour of some gas phase compounds and complexes in solutions is predicted for the gas chromatography and solvent extraction experiments. Redox potentials in aqueous solutions of these elements are estimated.

Electronic structure and properties of the group 4, 5 and 6 highest chlorides including elements 104, 105, and 106

Results of the Dirac-Slater discrete variational calculations for the group 4, 5, and 6 highest chlorides including elements 104, 105, and 106 have shown that the groups are not identical with respect to trends in the electronic structure and bonding. The charge density distribution data show that notwithstanding the basic increase in covalency within the groups this increase diminishes in going from group 4 to group 6. As a result, E106Cl_6 will be less stable toward thermal decomposition than WCl_6, which is confirmed by an estimated low E106-Cl bond energy. \delta H_form equal to -90.3 ± 6 kcal/rnol is obtained for E106Cl_6 in the gas phase, which is indicative of a very low stability of this compound. The stability of the maximum oxidation state is shown to decrease in the direction E104(+4) > E105(+5) > E106(+6).

Prediction of some thermodynamic properties of selected compounds of element 104

A set of parametrized equations has been published by Bratsch and Lagowski for calculating thermodynamic properties of the lanthanides, actinides, element 104, and certainrelated elements. Since these equations were applied to element 104, new values for the first four ionization energies and radii of the ions of charge +1, +2, +3, and +4 have been calculated for this element. The parametrized equations are used here with these new values to calculate some thermodynamic properties of element 104.

Predicted properties of the superheavy elements, II. Element 111, Eka-Gold

The chemical properties of element 111, eka-gold, are predicted through the use of the periodic table, relativistic Hartee-Fock-Slater calculations, and various qualitative theories which have established their usefulness in understanding and correlating properties of molecules. The results indicate that element 111 will be like Au(III) in its chemistry with little or no tendency to show stability in the I or II states. There is a possibility that the 111 - ion, analogous to the auride ion, will be stable.

Predicted properties of the superheavy elements. III. Element 115, Eka-Bismuth

Element 115 is expected to be in group V-a of the periodic table and have most stable oxidation states of I and III. The oxidation state of I, which plays a minor role in bismuth chemistry, should be a major factor in 115 chemistry. This change will arise because of the large relativistic splitting of the spherically symmetric 7p_l/2 shell from the 7P_3/2 shell. Element 115 will therefore have a single 7p_3/2 electron outside a 7p^2_1/2 closed shell. The magnitude of the first ionization energy and ionic radius suggest a chemistry similar to Tl^+. Similar considerations suggest that 115^3+ will have a chemistry similar to Bi^3+. Hydrolysis will therefore be easy and relatively strongly complexing anions of strong acids will be needed in general to effect studies of complexation chemistry. Some other properties of 115 predicted are as follows: ionization potentials I 5.2 eV, II 18.1 eV, III 27.4 eV, IV 48.5 eV, 0 \rightarrow 5^+ 159 eV; heat of sublimation, 34 kcal (g-atom)^-1; atomic radius, 2.0 A; ionic radius, 115^+ 1.5 A, 115^3+ 1.0 A; entropy, 16 cal deg^-1 (g-atom)^-l (25°); standard electrode potential 115^+ |115, -1.5 V; melting and boiling points are similar to element 113.