Solution of the time dependent Dirac-Fock-Slater equation for many-electron collisions systems using a time window method

We present a new scheme to solve the time dependent Dirac-Fock-Slater equation (TDDFS) for heavy many electron ion-atom collision systems. Up to now time independent self consistent molecular orbitals have been used to expand the time dependent wavefunction and rather complicated potential coupling matrix elements have been neglected. Our idea is to minimize the potential coupling by using the time dependent electronic density to generate molecular basis functions. We present the first results for 16 MeV S{^16+} on Ar.

Theory of laser-induced ultrafast structural changes in solids

The present thesis is a contribution to the study of laser-solid interaction. Despite the numerous applications resulting from the recent use of laser technology, there is still a lack of satisfactory answers to theoretical questions regarding the mechanism leading to the structural changes induced by femtosecond lasers in materials. We provide here theoretical approaches for the description of the structural response of different solids (cerium, samarium sulfide, bismuth and germanium) to femtosecond laser excitation. Particular interest is given to the description of the effects of the laser pulse on the electronic systems and changes of the potential energy surface for the ions. Although the general approach of laser-excited solids remains the same, the potential energy surface which drives the structural changes is calculated with different theoretical models for each material. This is due to the difference of the electronic properties of the studied systems. We use the Falicov model combined with an hydrodynamic method to study photoinduced phase changes in cerium. The local density approximation (LDA) together with the Hubbard-type Hamiltonian (LDA+U) in the framework of density functional theory (DFT) is used to describe the structural properties of samarium sulfide. We parametrize the time-dependent potential energy surface (calculated using DFT+ LDA) of bismuth on which we perform quantum dynamical simulations to study the experimentally observed amplitude collapse and revival of coherent $A_{1g}$ phonons. On the basis of a time-dependent potential energy surface calculated from a non-orthogonal tight binding Hamiltonian, we perform molecular dynamics simulation to analyze the time evolution (coherent phonons, ultrafast nonthermal melting) of germanium under laser excitation. The thermodynamic equilibrium properties of germanium are also reported. With the obtained results we are able to give many clarifications and interpretations of experimental results and also make predictions.

Many-electron relativistic calculation and interpretation of atomic processes in time dependent heavy-ion scattering

The time dependence of a heavy-ion-atom collision system is solved via a set of coupled channel equations using energy eigenvalues and matrix elements from a self-consistent field relativistic molecular many-electron Dirac-Fock-Slater calculation. Within this independent particle model we give a full many-particle interpretation by performing a small number of single-particle calculations. First results for the P(b) curves for the Ne K-hole excitation for the systems F{^8+} - Ne and F{^6+} - Ne as examples are discussed.

Electron-Nuclear Correlation in Laser-Driven Diatomic Molecules

The interaction of short intense laser pulses with atoms/molecules produces a multitude of highly nonlinear processes requiring a non-perturbative treatment. Detailed study of these highly nonlinear processes by numerically solving the time-dependent Schrodinger equation becomes a daunting task when the number of degrees of freedom is large. Also the coupling between the electronic and nuclear degrees of freedom further aggravates the computational problems. In the present work we show that the time-dependent Hartree (TDH) approximation, which neglects the correlation effects, gives unreliable description of the system dynamics both in the absence and presence of an external field. A theoretical framework is required that treats the electrons and nuclei on equal footing and fully quantum mechanically. To address this issue we discuss two approaches, namely the multicomponent density functional theory (MCDFT) and the multiconfiguration time-dependent Hartree (MCTDH) method, that go beyond the TDH approximation and describe the correlated electron-nuclear dynamics accurately. In the MCDFT framework, where the time-dependent electronic and nuclear densities are the basic variables, we discuss an algorithm to calculate the exact Kohn-Sham (KS) potentials for small model systems. By simulating the photodissociation process in a model hydrogen molecular ion, we show that the exact KS potentials contain all the many-body effects and give an insight into the system dynamics. In the MCTDH approach, the wave function is expanded as a sum of products of single-particle functions (SPFs). The MCTDH method is able to describe the electron-nuclear correlation effects as the SPFs and the expansion coefficients evolve in time and give an accurate description of the system dynamics. We show that the MCTDH method is suitable to study a variety of processes such as the fragmentation of molecules, high-order harmonic generation, the two-center interference effect, and the lochfrass effect. We discuss these phenomena in a model hydrogen molecular ion and a model hydrogen molecule. Inclusion of absorbing boundaries in the mean-field approximation and its consequences are discussed using the model hydrogen molecular ion. To this end, two types of calculations are considered: (i) a variational approach with a complex absorbing potential included in the full many-particle Hamiltonian and (ii) an approach in the spirit of time-dependent density functional theory (TDDFT), including complex absorbing potentials in the single-particle equations. It is elucidated that for small grids the TDDFT approach is superior to the variational approach.

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.

Variationelle Zeitdiskretisierungen höherer Ordnung der Konvektions-Diffusions-Reaktions-Gleichungen in zeitabhängigen Gebieten

Wir betrachten zeitabhängige Konvektions-Diffusions-Reaktions-Gleichungen in zeitabhängi- gen Gebieten, wobei die Bewegung des Gebietsrandes bekannt ist. Die zeitliche Entwicklung des Gebietes wird durch die ALE-Formulierung behandelt, die die Nachteile der klassischen Euler- und Lagrange-Betrachtungsweisen behebt. Die Position des Randes und seine Geschwindigkeit werden dabei so in das Gebietsinnere fortgesetzt, dass starke Gitterdeformationen verhindert werden. Als Zeitdiskretisierungen höherer Ordnung werden stetige Galerkin-Petrov-Verfahren (cGP) und unstetige Galerkin-Verfahren (dG) auf Probleme in zeitabhängigen Gebieten angewendet. Weiterhin werden das C 1 -stetige Galerkin-Petrov-Verfahren und das C 0 -stetige Galerkin- Verfahren vorgestellt. Deren Lösungen lassen sich auch in zeitabhängigen Gebieten durch ein einfaches einheitliches Postprocessing aus der Lösung des cGP-Problems bzw. dG-Problems erhalten. Für Problemstellungen in festen Gebieten und mit zeitlich konstanten Konvektions- und Reaktionstermen werden Stabilitätsresultate sowie optimale Fehlerabschätzungen für die nachbereiteten Lösungen der cGP-Verfahren und der dG-Verfahren angegeben. Für zeitabhängige Konvektions-Diffusions-Reaktions-Gleichungen in zeitabhängigen Gebieten präsentieren wir konservative und nicht-konservative Formulierungen, wobei eine besondere Aufmerksamkeit der Behandlung der Zeitableitung und der Gittergeschwindigkeit gilt. Stabilität und optimale Fehlerschätzungen für die in der Zeit semi-diskretisierten konservativen und nicht-konservativen Formulierungen werden vorgestellt. Abschließend wird das volldiskretisierte Problem betrachtet, wobei eine Finite-Elemente-Methode zur Ortsdiskretisierung der Konvektions-Diffusions-Reaktions-Gleichungen in zeitabhängigen Gebieten im ALE-Rahmen einbezogen wurde. Darüber hinaus wird eine lokale Projektionsstabilisierung (LPS) eingesetzt, um der Konvektionsdominanz Rechnung zu tragen. Weiterhin wird numerisch untersucht, wie sich die Approximation der Gebietsgeschwindigkeit auf die Genauigkeit der Zeitdiskretisierungsverfahren auswirkt.

Femtosecond time-resolved molecular multiphoton ionization and fragmentation of Na_2

We present a comparison between experimental and theoretical results for pump/probe multiphoton ionizing transitions of the sodium dimer, initiated by femtosecond laser pulses. It is shown that the motion of vibrational wavepackets in two electronic states is probed simultaneously and their dynamics is reflected in the total Na^+_2 ion signal which is recorded as a function of the time delay between pump and probe pulse. The time dependent quantum calculations demonstrate that two ionization pathways leading to the same final states of the molecularion exist: one gives an oscillating contribution to the ion signal, the other yields a constant background. From additional measurements of the Na^+ -transient photofragmentation spectrum it is deduced that another ionization process leading to different final ionic states exists. The process includes the excitation of a doubly excitedbound Rydberg state. This conclusion is supported by the theoretical simulation.

Physik der Elektronen in Nanostrukturen und niedrigdimensionalen Systemen - Studie zweier Beispiele

Als Beispiele für die vielfältigen Phänomene der Physik der Elektronen in niedrigdimensionalen Systemen wurden in dieser Arbeit das Cu(110)(2x1)O-Adsorbatsystem und die violette Li0.9Mo6O17-Bronze untersucht. Das Adsorbatsystem bildet selbstorganisierte quasi-eindimensionale Nanostrukturen auf einer Kupferoberfläche. Die Li-Bronze ist ein Material, das aufgrund seiner Kristallstruktur quasi-eindimensionale elektronische Eigenschaften im Volumen aufweist. Auf der Cu(110)(2x1)O-Oberfläche kann durch Variation der Sauerstoffbedeckung die Größe der streifenartigen CuO-Domänen geändert werden und damit der Übergang von zwei Dimensionen auf eine Dimension untersucht werden. Der Einfluss der Dimensionalität wurde anhand eines unbesetzten elektronischen Oberflächenzustandes studiert. Dessen Energieposition (untere Bandkante) verschiebt mit zunehmender Einschränkung (schmalere CuO-Streifen) zu größeren Energien hin. Dies ist ein bekannter quantenmechanischer Effekt und relativ gut verstanden. Zusätzlich wurde die Lebensdauer des Zustandes auf der voll bedeckten Oberfläche (zwei Dimensionen) ermittelt und deren Veränderung mit der Breite der CuO-Streifen untersucht. Es zeigt sich, dass die Lebensdauer auf schmaleren CuO-Streifen drastisch abnimmt. Dieses Ergebnis ist neu. Es kann im Rahmen eines Fabry-Perot-Modells als Streuung in Zustände außerhalb der CuO-Streifen verstanden werden. Außer den gerade beschriebenen Effekten war es möglich die Ladungsdichte des diskutierten Zustandes orts- und energieabhängig auf den CuO-Streifen zu studieren. Die Li0.9Mo6O17-Bronze wurde im Hinblick auf das Verhalten der elektronischen Zustandsdichte an der Fermikante untersucht. Diese Fragestellung ist besonders wegen der Quasieindimensionalität des Materials interessant. Die Messungen von STS-Spektren in der Nähe der Fermienergie zeigen, dass die Elektronen in der Li0.9Mo6O17-Bronze eine sogenannte Luttingerflüssigkeit ausbilden, die anstatt einer Fermiflüssigkeit in eindimensionalen elektronischen Systemen erwartet wird. Bisher wurde Luttingerflüssigkeitsverhalten erst bei wenigen Materialien und Systemen experimentell nachgewiesen, obschon die theoretischen Voraussagen mehr als 30 Jahre zurückliegen. Ein Charakteristikum einer Luttingerflüssigkeit ist die Abnahme der Zustandsdichte an der Fermienergie mit einem Potenzgesetz. Dieses Verhalten wurde in STS-Spektren dieser Arbeit beobachtet und quantitativ im Rahmen eines Luttingerflüssigkeitsmodells beschrieben. Auch die Temperaturabhängigkeit des Phänomens im Bereich von 5K bis 55K ist konsistent mit der Beschreibung durch eine Luttingerflüssigkeit. Generell zeigen die Untersuchungen dieser Arbeit, dass die Dimensionalität, insbesondere deren Einschränkung, einen deutlichen Einfluss auf die elektronischen Eigenschaften von Systemen und Materialien haben kann.

Intrinsic channel closing in strong-field single ionization of H2

Die Ionisation von H2 in intensiven Laserpulsen wird mit Hilfe der numerischen Integration der zeitabhängigen Schrödingergleichung für ein Einelektronenmodell untersucht, das die Vibrationsbewegung berücksichtigt. Die Spektren der kinetischen Elektronenenergie hängen stark von der Vibrationsquantenzahl des erzeugten H2+ Ions ab. Für bestimmte Vibrationszustände ist die Ausbeute der Elektronen in der Mitte des Plateaus stark erhöht. Der Effekt wird "channel closings" zugeschrieben, die in Atomen durch Variation der Laserintensität beobachtet wurden. The ionization of H2 in intense laser pulses is studied by numerical integration of the time-dependent Schrödinger equation for a single-active-electron model including the vibrational motion. The electron kinetic energy spectra in high-order above-threshold ionization are strongly dependent on the vibrational quantum number of the created H2+ ion. For certain vibrational states, the electron yield in the mid-plateau region is strongly enhanced. The effect is attributed to channel closings, which were previously observed in atoms by varying the laser intensity.

Relativistic ab initio calculations for ion-atom collisions

Within the independent particle model we solve the time-dependent single-particle equation using ab initio SCF-DIRAC-FOCK-SLATER wavefunctions as a basis. To reinstate the many-particle aspect of the collision system we use the inclusive probability formalism to answer experimental questions. As an example we show an application to the case of S{^15+} on Ar where experimental data on the K-K charge transfer are available for a wide range of impact energies from 4.7 to 90 MeV. Our molecular adiabatic calculations and the evaluation using the inclusive probability formalism show good results in the low energy range from 4.7 to 16 MeV impact energy.

Inclusive probabilities for the scattering system 16 MeV-S{^16+} on Ar

We performed ab initio calculations of many particle inclusive probabilities for the scattering system 16 MeV-S{^16+} on Ar. The solution of the time-dependent DIRAC-FOCK-SLATER-equation is achieved via a set of coupled-channel equations with energy eigenvalues and matrix elements which are given by static SCF molecular many electron calculations.

Multiple vacancy inclusive probabilities for the scattering system 16 MeV S{^16+} on Ar

To evaluate single and double K-shell inclusive charge transfer probabilities in ion-atom collisions we solve the time-dependent Dirac equation. By expanding the timedependent wavefunction in a set of molecular basis states the time-dependent equation reduces to a set of coupled-channel equations. The energy eigenvalues and matrix elements are taken from self-consistent relativistic molecular many-electron Dirac-Fock-Slater calculations. We present many-electron inclusive probabilities for different final configurations as a function of impact parameter for single and double K-shell vacancy production in collisions of bare S on Ar.

Tomographic imaging of molecular orbitals in length and velocity form

Recently Itatani et al. [Nature 432, 876 (2004)] introduced the new concept of molecular orbital tomography, where high harmonic generation (HHG) is used to image electronic wave functions. We describe an alternative reconstruction form, using momentum instead of dipole matrix elements for the electron recombination step in HHG. We show that using this velocity-form reconstruction, one obtains better results than using the original length-form reconstruction. We provide numerical evidence for our claim that one has to resort to extremely short pulses to perform the reconstruction for an orbital with arbitrary symmetry. The numerical evidence is based on the exact solution of the time-dependent Schrödinger equation for 2D model systems to simulate the experiment. Furthermore we show that in the case of cylindrically symmetric orbitals, such as the N2 orbital that was reconstructed in the original work, one can obtain the full 3D wave function and not only a 2D projection of it. Vor kurzem führten Itatani et al. [Nature 432, 876 (2004)] das Konzept der Molelkülorbital-Tomographie ein. Hierbei wird die Erzeugung hoher Harmonischer verwendet, um Bilder von elektronischen Wellenfunktionen zu gewinnen. Wir beschreiben eine alternative Form der Rekonstruktion, die auf Impuls- statt Dipol-Matrixelementen für den Rekombinationsschritt bei der Erzeugung der Harmonischen basiert. Wir zeigen, dass diese "Geschwindigkeitsform" der Rekonstruktion bessere Ergebnisse als die ursprüngliche "Längenform" liefert. Wir zeigen numerische Beweise für unsere Behauptung, dass man zu extrem kurzen Laserpulsen gehen muss, um Orbitale mit beliebiger Symmetrie zu rekonstruieren. Diese Ergebnisse basieren auf der exakten Lösung der zeitabhängigen Schrödingergleichung für 2D-Modellsysteme. Wir zeigen ferner, dass für zylindersymmetrische Orbitale wie das N2-Orbital, welches in der oben zitierten Arbeit rekonstruiert wurde, das volle 3D-Orbital rekonstruiert werden kann, nicht nur seine 2D-Projektion.

Realistic many electron coupled channel calculations in heavy ion scattering

The time dependent Dirac equation which describes a heavy ion-atom collision system is solved via a set of coupled channel equations with energy eigenvalues and matrix elements which are given by a selfconsistent field many electron calculation. After a brief discussion of the theoretical approximations and the connection of the many particle with the one particle interpretation we discuss first results for the systems F{^8+} - Ne and F{^6+} - Ne. The resulting P(b) curves for the creation of a Ne K-hole are in good agreement with the experimental results.

High laser field effects in multiphoton ionization of Na_2

Femtosecond pump/probe multiphoton ionization experiments on Na_2 molecules are performed. The dependence of the total Na^+_2 ion signal on the delay time and the intensity of the femtosecond laser pulses is studied in detail. It is observed that molecular vibrational wavepacket motion in different electronic states dominates the time dependence of the ion signal. For higher laser intensities the relative contributions from the A ^1 \summe^+_u and the 2 ^1 \produkt__g states change dramatically, indicating the increasing importance of a two-electron versus a one-electron process. For even stronger fields (10 ^12 W/ cm²) a vibrational wavepacket in the electronic ground state X ^1 \summe^+_g is formed and its dynamics is also observed in the transient Na^+_2 signal. Time-dependent quantum calculations are presented. The theoretical results agree well with the experiment.