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.

An ab-initio calculation of the Coulomb explosion of N_2 after heavy-ion bombardement

Self-consistent-field calculations for the total potential energy of highly ionized N_2 molecules are presented. We compare these calculations to the experimentally observed energy released in the Coulomb explosion of ionized N_2 molecules created after collision with fast heavy ions. The most important electronic states of the fragment ions are determined.

Non-relativistic and relativistic molecular calculations for the chalcogen hexafluorides: SF_6, SeF_6, TeF_6, PoF_6

Non-relativistic Hartree-Fock-Slater and relativistic Dirac-Slater self-consistent orbital models are applied for the analysis of the electronic structure of the chalcogen hexafluorides: SF_6, SeF_6, TeF_6 and PoF_6. The molecular eigenfunctions and eigenvalues are generated using the discrete variational method (DVM) with numerical basis functions. The results obtained for SF_6 are compared with other ab initio calculations. Information about relativistic level shifts and spin-orbit splitting has been obtained by comparison between the non-relativistic and relativistic results.

Ab initio Molecular Dynamics Simulations of the Structural Response of Solids to Ultrashort Laser and XUV Pulses

The theoretical model and underlying physics described in this thesis are about the interaction of femtosecond-laser and XUV pulses with solids. The key to understand the basics of such interaction is to study the structural response of the materials after laser interaction. Depending on the laser characteristics, laser-solid interaction can result in a wide range of structural responses such as solid-solid phase transitions, vacuum phonon squeezing, ultrafast melting, generation of coherent phonons, etc. During my research work, I have modeled the systems irradiated by low-, medium- and high-laser intensities, and studied different types of structural dynamics of solids at various laser fluences.

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.

Threshold behaviour of L x-ray excitation in Xe - Ag collisions

In the collision system Xe - Ag, the thresholds for excitation of quasimolecular L radiation and characteristic Ag L radiation have been found to lie at about 5 MeV and 1 MeV, respectively. These results are discussed on the basis of ab initio calculations of the screened interaction potential and the electron-correlation diagram.

Quadrupole moments of radium isotopes from the 7p^2P_3/2 hyperfine structure in Ra II

The hyperfine structure and isotope shift of ^{221- 226}Ra and ^{212, 214}Ra have been measured in the ionic (Ra 11) transition 7s^2 S_{1/2} - 7p ^2 P_{3/2} (\lamda = 381.4 nm). The method of on-line collinear fast-beam laser spectroscopy has been applied using frequency-doubling of cw dye laser radiation in an external ring cavity. The magnetic hyperfine fields are compared with semi-empirical and ab initio calculations. The analysis of the quadrupole splitting by the same method yields the following, improved values of spectroscopic quadrupole moments: Q_s(^221 Ra)= 1.978(7)b, Q_s (^223 Ra)= 1.254(3)b and the reanalyzed values Q_s(^209 Ra) = 0.40(2)b, Q_s(^211 Ra) = 0.48(2)b, Q_s(^227 Ra)= 1.58(3)b, Q_s (^229 Ra) = 3.09(4)b with an additional scaling uncertainty of ±5%. Furthermore, the J-dependence of the isotope shift is analyzed in both Ra II transitions connecting the 7s^2 S_{1/2} ground state with the first excited doublet 7p^ P_{1/2} and 7p^ P_{3/2}.

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.

Atomistic-continuum modeling of ultrafast laser-induced melting of silicon targets

In this work, we present an atomistic-continuum model for simulations of ultrafast laser-induced melting processes in semiconductors on the example of silicon. The kinetics of transient non-equilibrium phase transition mechanisms is addressed with MD method on the atomic level, whereas the laser light absorption, strong generated electron-phonon nonequilibrium, fast heat conduction, and photo-excited free carrier diffusion are accounted for with a continuum TTM-like model (called nTTM). First, we independently consider the applications of nTTM and MD for the description of silicon, and then construct the combined MD-nTTM model. Its development and thorough testing is followed by a comprehensive computational study of fast nonequilibrium processes induced in silicon by an ultrashort laser irradiation. The new model allowed to investigate the effect of laser-induced pressure and temperature of the lattice on the melting kinetics. Two competing melting mechanisms, heterogeneous and homogeneous, were identified in our big-scale simulations. Apart from the classical heterogeneous melting mechanism, the nucleation of the liquid phase homogeneously inside the material significantly contributes to the melting process. The simulations showed, that due to the open diamond structure of the crystal, the laser-generated internal compressive stresses reduce the crystal stability against the homogeneous melting. Consequently, the latter can take a massive character within several picoseconds upon the laser heating. Due to the large negative volume of melting of silicon, the material contracts upon the phase transition, relaxes the compressive stresses, and the subsequent melting proceeds heterogeneously until the excess of thermal energy is consumed. A series of simulations for a range of absorbed fluences allowed us to find the threshold fluence value at which homogeneous liquid nucleation starts contributing to the classical heterogeneous propagation of the solid-liquid interface. A series of simulations for a range of the material thicknesses showed that the sample width we chosen in our simulations (800 nm) corresponds to a thick sample. Additionally, in order to support the main conclusions, the results were verified for a different interatomic potential. Possible improvements of the model to account for nonthermal effects are discussed and certain restrictions on the suitable interatomic potentials are found. As a first step towards the inclusion of these effects into MD-nTTM, we performed nanometer-scale MD simulations with a new interatomic potential, designed to reproduce ab initio calculations at the laser-induced electronic temperature of 18946 K. The simulations demonstrated that, similarly to thermal melting, nonthermal phase transition occurs through nucleation. A series of simulations showed that higher (lower) initial pressure reinforces (hinders) the creation and the growth of nonthermal liquid nuclei. For the example of Si, the laser melting kinetics of semiconductors was found to be noticeably different from that of metals with a face-centered cubic crystal structure. The results of this study, therefore, have important implications for interpretation of experimental data on the kinetics of melting process of semiconductors.

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.

New calculations of P(b) curves for 1s_\sigma excitation in low-Z (Z\le10) ion-atom scattering

Ab initio fully relativistic SCF molecular calculations of energy eigenvalues as well as coupling-matrix elements are used to calculate the 1s_\sigma excitation differential cross section for Ne-Ne and Ne-O in ion-atom collisions. A relativistic perturbation treatment which allows a direct comparison with analogous non-relativistic calculations is also performed.

Inclusive probability calculations for the K-vacancy transfer in collisions of S{^15+} on Ar

Using the single-particle amplitudes from a 20-level coupled-channel calculation with ab initio relativistic self consistent LCAO-MO Dirac-Fock-Slater energy eigenvalues and matrix elements we calculate within the frame of the inclusive probability formalism impact-parameter-dependent K-hole transfer probabilities. As an example we show results for the heavy asymmetric collision system S{^15+} on Ar for impact energies from 4.7 to 16 MeV. The inclusive probability formalism which reinstates the many-particle aspect of the collision system permits a qualitative and quantitative agreement with the experiment which is not achieved by the single-particle picture.

{2p_\pi}-{2p_\sigma} crossing in heavy symmetric ion-atom collisions: I. Level structure

Ab initio self-consistent DFS calculations are performed for five different symmetric atomic systems from Ar-Ar to Pb-Pb. The level structure for the {2p_\pi}-{2p_\sigma} crossing as function of the united atomic charge Z_u is studied and interpreted. Manybody effects, spin-orbit splitting, direct relativistic effects as well as indirect relativistic effects are differently important for different Z_u. For the I-I system a comparison with other calculations is given.

Influence of occupation number of single particle levels on K - K charge transfer in collisions of 90 keV - Ne{^9+} on Ne

The influence of the occupation of the single particle levels on the impact parameter dependent K - K charge transfer occuring in collisions of 90 keV Ne{^9+} on Ne was studied using coupled channel calculations. The energy eigenvalues and matrixelements for the single particle levels were taken from ab initio self consistent MO-LCAO-DIRAC-FOCK-SLATER calculations with occupation numbers corresponding to the single particle amplitudes given by the coupled channel calculations.