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.

Computation of relativistic symmetry orbitals for finite double point groups

A program is presented for the construction of relativistic symmetry-adapted molecular basis functions. It is applicable to 36 finite double point groups. The algorithm, based on the projection operator method, automatically generates linearly independent basis sets. Time reversal invariance is included in the program, leading to additional selection rules in the non-relativistic limit.

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.

A clear atomic example for the surface sensitivity of penning ionization

Using a crossed-beam apparatus with a double hemispherical electron spectrometer, we have studied the spectrum of electrons released in thermal energy ionizing collisions of metastable He^*(2^3S) atoms with ground state Yb(4f^14 6s^2 ^1S_0) atoms, thereby providing the first Penning electron spectrum of an atomic target with-4f-electrons. In contrast to the HeI (58.4nm) and NeI (73.6/74.4nm) photoelectron spectra of Yb, which show mainly 4f- and 6s-electron emission in about a 5:1 ratio, the He^*(2^3S) Penning electron spectrum is dominated by 6s-ionization, acoompnnied by some correlation- induced 6p-emission (8% Yb+( 4f^14 6p^2P) formation) and very little 4f-ionization (<_~ 2.5%). This astounding result is attributed to the electron exchange mechanism for He^*(2^3S) ionization and reflects the poor overlap of the target 4f-electron wavefunction with the 1s-hole of He^*(2^3S), as discussed on thc basis of Dirac-Fock wave functions for the Yb orbitals and through calculations of the partial ionization cross sections involving semiempirical complex potentiale. The presented case may be regarded as the elearest atomic example for the surface sensitivity of He^*(2^3S) Penning ionization observed so far.

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] 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.

Reduced L_1 level width and Coster-Kronig yields by relaxation and continuum interactions in atomic zinc

Calculations of the level width \gamma( L_1) and the f_12 and f_13 Coster-Kronig yields for atomic zinc have been performed with Dirac-Fock wave functions. For \gamma(L_1), a large deviation between theory and evaluated data exists. We include the incomplete orthogonality of the electron orbitals as well as the interchannel interaction of the decaying states. Orbital relaxation reduces the total rates in all groups of the electron-emission spectrum by about 10-20 %. Different, however, is the effect of the continuum interaction. The L_1-L_23X Coster-Kronig part of the spectrum is definitely reduced in its intensity, whereas the MM and MN spectra are slightly enhanced. This results in a reduction of Coster-Kronig yields, where for medium and heavy elements considerable discrepancies have been found in comparison to relativistic theory. Briefly, we discuss the consequences of our calculations for heavier elements.

Relativistic Dirac-Fock-Slater program to calculate potential-energy curves for diatomic molecules

A LCAO-MO (linear combination of atomic orbitals - molecular orbitals) relativistic Dirac-Fock-Slater program is presented, which allows one to calculate accurate total energies for diatomic molecules. Numerical atomic Dirac-Fock-Slater wave functions are used as basis functions. All integrations as well as the solution of the Poisson equation are done fully numerical, with a relative accuracy of 10{^-5} - 10{^-6}. The details of the method as well as first results are presented here.

Calculations of the polycentric linear molecule H^2+_3 with the finite element method

A fully numerical two-dimensional solution of the Schrödinger equation is presented for the linear polyatomic molecule H^2+_3 using the finite element method (FEM). The Coulomb singularities at the nuclei are rectified by using both a condensed element distribution around the singularities and special elements. The accuracy of the results for the 1\sigma and 2\sigma orbitals is of the order of 10^-7 au.

Ground state correlation energy of the Be-sequence for Z = 4-20 in MCDF approximation

The ground state (J = 0) electronic correlation energy of the 4-electron Be-sequence is calculated in the Multi-Configuration Dirac-Fock approximation for Z = 4-20. The 4 electrons were distributed over the configurations arising from the 1s, 2s, 2p, 3s, 3p and 3d orbitals. Theoretical values obtained here are in good agreement with experimental correlation energies.

Atomic and ionic radii of superheavy elements

Atomic and ionic radii are presented for the elements E104-E120 and E156-E172. It is shown that a number of effects correlated with the large relativistic contraction of orbitals with low angular momentum leads to smaller atoms for higher atomic numbers. It is expected that Cs is the largest atom in nature.

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

Results of relativistic (Dirac-Slater and Dirac-Fock) and nonrelativistic (Hartree-Fock-Slater) atomic and molecular calculations have been compared for the group 5 elements Nb, Ta, and Ha and their compounds MCl_5, to elucidate the influence of relativistic effects on their properties especially in going from the 5d element Ta to the 6d element Ha. The analysis of the radial distribution of the valence electrons of the metals for electronic configurations obtained as a result of the molecular calculations and their overlap with ligands show opposite trends in behavior for ns_1/2, np_l/2, and (n -1 )d_5/2 orbitals for Ta and Ha in the relativistic and nonrelativistic cases. Relativistic contraction and energetic stabilization of the ns_1/2 and np_l/2 wave functions and expansion and destabilization of the (n-1)d_5/2 orbitals make hahnium pentahalide more covalent than tantalum pentahalide and increase the bond strength. The nonrelativistic treatment of the wave functions results in an increase in ionicity of the MCl_5 molecules in going from Nb to Ha making element Ha an analog of V. Different trends for the relativistic and nonrelativistic cases are also found for ionization potentials, electronic affinities, and energies of charge-transfer transitions as well as the stability of the maximum oxidation state.

Calculation of isomer shift in Mössbauer spectroscopy

The approximations normally used in the calculation of the isomer shift are compared with the exact expressions using Dirac-Slater orbitals and a three-parameter Fermi-type nuclear charge distribution. The nonuniformity of the electronic density over the nuclear volume affects the results. Different choices of the nuclear surface thickness t and the radius c in the protonic density P_N (\gamma) also affects the isomer shift differently even though the values are chosen to yield a given value of \delta . The change in the electronic charge density which is caused by the alteration of P_N (\gamma) in the ground state and excited state of the nucleus is discussed using two extrememodels and the possible influence on the observable isomer shift is estimated.

Effect of nuclear motion on the energy eigenvalues in muonic atoms

I have investigated the effect of the nuclear motion on the energy eigenvalues in muonic atoms. In addition to the usually used reduced-mass correction, I have calculated the relativistic influences including the magnetic and retardation interaction between the nucleus and the muon for the inner orbitals of muonic atoms.

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.