Nuclear Charge Radii of the Nickel Isotopes

Collinear laser spectroscopy is performed on the nickel isotopes ^{58-68,70}Ni, using a time-resolved photon counting system. From the measured isotope shifts, nuclear charge radii R_{c} are extracted and compared to theoretical results. Three ab initio approaches all employ, among others, the chiral interaction NNLO_{sat}, which allows an assessment of their accuracy. We find agreement with experiment in differential radii δ⟨r_{c}^{2}⟩ for all employed ab initio methods and interactions, while the absolute radii are consistent with data only for NNLO_{sat}. Within nuclear density functional theory, the Skyrme functional SV-min matches experiment more closely than the Fayans functional Fy(Δr,HFB).

Laser Spectroscopy of Neutron-Rich Tin Isotopes: A Discontinuity in Charge Radii across the N=82 Shell Closure.

The change in mean-square nuclear charge radii δ⟨r^{2}⟩ along the even-A tin isotopic chain ^{108-134}Sn has been investigated by means of collinear laser spectroscopy at ISOLDE/CERN using the atomic transitions 5p^{2} ^{1}S_{0}→5p6 s^{1}P_{1} and 5p^{2} ^{3}P_{0}→5p6s ^{3}P_{1}. With the determination of the charge radius of ^{134}Sn and corrected values for some of the neutron-rich isotopes, the evolution of the charge radii across the N=82 shell closure is established. A clear kink at the doubly magic ^{132}Sn is revealed, similar to what has been observed at N=82 in other isotopic chains with larger proton numbers, and at the N=126 shell closure in doubly magic ^{208}Pb. While most standard nuclear density functional calculations struggle with a consistent explanation of these discontinuities, we demonstrate that a recently developed Fayans energy density functional provides a coherent description of the kinks at both doubly magic nuclei, ^{132}Sn and ^{208}Pb, without sacrificing the overall performance. A multiple correlation analysis leads to the conclusion that both kinks are related to pairing and surface effects.

From Calcium to Cadmium: Testing the Pairing Functional through Charge Radii Measurements of

Differences in mean-square nuclear charge radii of ^{100-130}Cd are extracted from high-resolution collinear laser spectroscopy of the 5s ^{2}S_{1/2}→5p ^{2}P_{3/2} transition of the ion and from the 5s5p ^{3}P_{2}→5s6s ^{3}S_{1} transition in atomic Cd. The radii show a smooth parabolic behavior on top of a linear trend and a regular odd-even staggering across the almost complete sdgh shell. They serve as a first test for a recently established new Fayans functional and show a remarkably good agreement in the trend as well as in the total nuclear charge radius.

Individual Low-Energy Toroidal Dipole State in

The low-energy dipole excitations in ^{24}Mg are investigated within the Skyrme quasiparticle random phase approximation for axial nuclei. The calculations with the force SLy6 reveal a remarkable feature: the lowest I^{π}K=1^{-}1 excitation (E=7.92 MeV) in ^{24}Mg is a vortical toroidal state (TS) representing a specific vortex-antivortex realization of the well-known spherical Hill's vortex in a strongly deformed axial confinement. This is a striking example of an individual TS which can be much more easily discriminated in experiment than the toroidal dipole resonance embracing many states. The TS acquires the lowest energy due to the huge prolate axial deformation in ^{24}Mg. The result persists for different Skyrme parametrizations (SLy6, SVbas, SkM*). We analyze spectroscopic properties of the TS and its relation with the cluster structure of ^{24}Mg. Similar TSs could exist in other highly prolate light nuclei. They could serve as promising tests for various reactions to probe a vortical (toroidal) nuclear flow.

Towards the analysis of attosecond dynamics in complex systems.

We study from a theoretical perspective the ionization of molecules and clusters induced by irradiation of a combined two-color laser field consisting of a train of attosecond XUV pulses in the presence of an IR field. We use time-dependent density-functional theory (TDDFT) in real time and real space as a theoretical tool. The calculated results are compared to experimental data when available. We also compare TDDFT with results obtained using a time-dependent Schrödinger equation (TDSE), which is well suited to simple systems while TDDFT allows dealing with more complex molecules and clusters. As a key observable, we study ionization versus delay time of the XUV pulses with respect to the IR background pulse. Experiments in simple atoms (He and Ar) show a regular modulation of this signal with half the IR period. This feature is recovered by TDDFT as well as by the TDSE (although total ionization differs by an order of magnitude). As more complex systems, we consider a C3 chain molecule and Na clusters. Here we encounter a different picture as the ionization signal develops a more involved pattern with several peaks per half IR period and as the TDSE produces a different pattern to TDDFT. Both effects could be related to the appearance of strong resonance modes in these more complex systems.

Charge Radii of Neutron Deficient

Bunched-beam collinear laser spectroscopy is performed on neutron deficient ^{52,53}Fe prepared through in-flight separation followed by a gas stopping. This novel scheme is a major step to reach nuclides far from the stability line in laser spectroscopy. Differential mean-square charge radii δ⟨r^{2}⟩ of ^{52,53}Fe are determined relative to stable ^{56}Fe as δ⟨r^{2}⟩^{56,52}=-0.034(13) fm^{2} and δ⟨r^{2}⟩^{56,53}=-0.218(13) fm^{2}, respectively, from the isotope shift of atomic hyperfine structures. The multiconfiguration Dirac-Fock method is used to calculate atomic factors to deduce δ⟨r^{2}⟩. The values of δ⟨r^{2}⟩ exhibit a minimum at the N=28 neutron shell closure. The nuclear density functional theory with Fayans and Skyrme energy density functionals is used to interpret the data. The trend of δ⟨r^{2}⟩ along the Fe isotopic chain results from an interplay between single-particle shell structure, pairing, and polarization effects and provides important data for understanding the intricate trend in the δ⟨r^{2}⟩ of closed-shell Ca isotopes.

Finite Nuclei in the Quark-Meson Coupling Model.

We report the first use of the effective quark-meson coupling (QMC) energy density functional (EDF), derived from a quark model of hadron structure, to study a broad range of ground state properties of even-even nuclei across the periodic table in the nonrelativistic Hartree-Fock+BCS framework. The novelty of the QMC model is that the nuclear medium effects are treated through modification of the internal structure of the nucleon. The density dependence is microscopically derived and the spin-orbit term arises naturally. The QMC EDF depends on a single set of four adjustable parameters having a clear physics basis. When applied to diverse ground state data the QMC EDF already produces, in its present simple form, overall agreement with experiment of a quality comparable to a representative Skyrme EDF. There exist, however, multiple Skyrme parameter sets, frequently tailored to describe selected nuclear phenomena. The QMC EDF set of fewer parameters, derived in this work, is not open to such variation, chosen set being applied, without adjustment, to both the properties of finite nuclei and nuclear matter.

Towards time-dependent current-density-functional theory in the non-linear regime.

Time-Dependent Density-Functional Theory (TDDFT) is a well-established theoretical approach to describe and understand irradiation processes in clusters and molecules. However, within the so-called adiabatic local density approximation (ALDA) to the exchange-correlation (xc) potential, TDDFT can show insufficiencies, particularly in violently dynamical processes. This is because within ALDA the xc potential is instantaneous and is a local functional of the density, which means that this approximation neglects memory effects and long-range effects. A way to go beyond ALDA is to use Time-Dependent Current-Density-Functional Theory (TDCDFT), in which the basic quantity is the current density rather than the density as in TDDFT. This has been shown to offer an adequate account of dissipation in the linear domain when the Vignale-Kohn (VK) functional is used. Here, we go beyond the linear regime and we explore this formulation in the time domain. In this case, the equations become very involved making the computation out of reach; we hence propose an approximation to the VK functional which allows us to calculate the dynamics in real time and at the same time to keep most of the physics described by the VK functional. We apply this formulation to the calculation of the time-dependent dipole moment of Ca, Mg and Na2. Our results show trends similar to what was previously observed in model systems or within linear response. In the non-linear domain, our results show that relaxation times do not decrease with increasing deposited excitation energy, which sets some limitations to the practical use of TDCDFT in such a domain of excitations.

Probing time-dependent molecular dipoles on the attosecond time scale.

Photoinduced molecular processes start with the interaction of the instantaneous electric field of the incident light with the electronic degrees of freedom. This early attosecond electronic motion impacts the fate of the photoinduced reactions. We report the first observation of attosecond time scale electron dynamics in a series of small- and medium-sized neutral molecules (N(2), CO(2), and C(2)H(4)), monitoring time-dependent variations of the parent molecular ion yield in the ionization by an attosecond pulse, and thereby probing the time-dependent dipole induced by a moderately strong near-infrared laser field. This approach can be generalized to other molecular species and may be regarded as a first example of molecular attosecond Stark spectroscopy.

On transition rates in surface hopping.

Trajectory surface hopping (TSH) is one of the most widely used quantum-classical algorithms for nonadiabatic molecular dynamics. Despite its empirical effectiveness and popularity, a rigorous derivation of TSH as the classical limit of a combined quantum electron-nuclear dynamics is still missing. In this work, we aim to elucidate the theoretical basis for the widely used hopping rules. Naturally, we concentrate thereby on the formal aspects of the TSH. Using a Gaussian wave packet limit, we derive the transition rates governing the hopping process at a simple avoided level crossing. In this derivation, which gives insight into the physics underlying the hopping process, some essential features of the standard TSH algorithm are retrieved, namely (i) non-zero electronic transition rate ("hopping probability") at avoided crossings; (ii) rescaling of the nuclear velocities to conserve total energy; (iii) electronic transition rates linear in the nonadiabatic coupling vectors. The well-known Landau-Zener model is then used for illustration.

Self-consistent calculations of the electric giant dipole resonances in light and heavy nuclei.

While bulk properties of stable nuclei are successfully reproduced by mean-field theories employing effective interactions, the dependence of the centroid energy of the electric giant dipole resonance on the nucleon number A is not. This problem is cured by considering many-particle correlations beyond mean-field theory, which we do within the quasiparticle time blocking approximation. The electric giant dipole resonances in 16O, 40Ca, and 208Pb are calculated using two new Skyrme interactions.

Existence of an exotic torus configuration in high-spin excited states of 40Ca.

We investigate the possibility of the existence of the exotic torus configuration in the high-spin excited states of (40)Ca. We here consider the spin alignments about the symmetry axis. To this end, we use a three-dimensional cranked Skyrme Hartree-Fock method and search for stable single-particle configurations. We find one stable state with the torus configuration at the total angular momentum J=60 h and an excitation energy of about 170 MeV in all calculations using various Skyrme interactions. The total angular momentum J=60 h consists of aligned 12 nucleons with the orbital angular momenta Λ=+4, +5, and +6 for spin-up or -down neutrons and protons. The obtained results strongly suggest that a macroscopic amount of circulating current breaking the time-reversal symmetry emerges in the high-spin excited state of (40)Ca.

Time-dependent density-functional theory with a self-interaction correction.

We discuss an implementation of the self-interaction correction for the local-density approximation to time-dependent density-functional theory. A variational formulation is given, taking care of the necessary constraints. A manageable and transparent propagation scheme using two sets of wave functions is proposed and applied to laser excitation with subsequent ionization of a dimer molecule.

Wave packet simulation of dense hydrogen.

Dense hydrogen is studied in the framework of wave packet molecular dynamics. In this semiquantal many-body simulation method the electrons are represented by wave packets which are suitably parametrized. The equilibrium properties and time evolution of the system are obtained with the help of a variational principle. At room temperature the results for the isotherms are in good agreement with anvil experiments. At higher densities beyond the range of the experimental data a transition from a molecular to a metallic state is predicted. The wave packets become delocalized and the electrical conductivity increases sharply. The phase diagram is calculated in a wide range of the pressure-density-temperature space. The observed transition from the molecular to metallic state is accompanied by an increase in density in agreement with recent reverberating shock wave experiments.

From self-consistent covariant effective field theories to their galilean-invariant counterparts.

We discuss how to obtain the nonrelativistic limit of a self-consistent relativistic effective field theory for dynamic problems. It is shown that the standard v/c expansion yields Galilean invariance only to first order in v/c, whereas second order is required to obtain important contributions such as the spin-orbit force. We propose a modified procedure which is a mapping rather than a strict v/c expansion.

Role of boundary conditions in dynamic studies of nuclear giant resonances and collisions.

Absorbing boundary conditions are often employed in time-dependent mean-field calculations to cope with the problem of emitted particles which would otherwise return back onto the system and falsify the dynamical evolution. We scrutinize two widely used methods, imaginary potentials and gradual attenuation by a mask function. To that end, we consider breathing oscillations of a nucleus computed on a radial one-dimensional grid in coordinate space. The most critical test case is the computation of resonance spectra in the (linear) domain of small amplitude motion. Absorbing bounds turn out to provide a reliable alternative to fully fledged continuum random phase approximation (RPA) calculations, although rather large absorbing bounds are required to simulate reliably well continuum conditions. We also investigate the computation of observables in the nonlinear domain. This regime turns out to be less demanding. Smaller absorbing margin suffice to achieve the wanted absorption effect.

Simple DFT model of clusters embedded in rare gas matrix: trapping sites and spectroscopic properties of Na embedded in Ar.

We present a theoretical model to study the dynamics of metallic clusters embedded in a rare gas matrix. We describe the active electrons of the embedded cluster using time dependent density functional theory, while the surrounding matrix is described in terms of classical molecular dynamics of polarizable atoms. The coupling between the cluster and the rare gas atoms is deduced from the work of Gross and Spiegelmann [J. Chem. Phys. 108, 4148 (1998)] and reformulated explicitly in a simple and efficient density functional form. The electron rare gas interaction takes the form of an averaged dipole fluctuation term, which retains the van der Waals long range interaction, and a short range repulsive pseudopotential, which accounts for the Pauli repulsion of the electron by the rare gas atom. We applied our model to Na clusters embedded in Ar matrix. For the latter we developed an efficient local pseudopotential, which allows studying systems containing more than 10(3) Ar atoms. We show that large systems are indeed necessary to account properly for long range polarization of the matrix, that competes with the matrix confinement effect. We focus our study on Na(2), Na(4), and Na(8). For each system, we have determined the geometry of the most favorable trapping site by means of damped molecular dynamics. We present the effect of matrix embedding on the optical absorption spectrum. For Na(2), the trapping site can be unambiguously identified by comparison of the absorption spectrum with experiment. For Na(4) the spectrum of the embedded cluster is significantly different from the free cluster spectrum, while for Na(8) differences are less pronounced.

Crossed beam pump and probe dynamics in metal clusters.

We investigate theoretically pump and probe dynamics in metal clusters with crossed-polarization laser beams. We explore the relation between the pronounced Mie plasmon resonance and the laser frequency. The resonance moves with the cluster radius and splits according to the actual deformation. We demonstrate that probe pulses with different (linear) polarization axes allow one to resolve the global shape oscillations of the cluster in monopole and quadrupole degrees of freedom.

Twist mode in spherical alkali metal clusters.

A remarkable orbital quadrupole magnetic resonance, so-called twist mode, is predicted in alkali metal clusters where it is represented by Ipi = 2(-) low-energy excitations of valence electrons with strong M2 transitions to the ground state. We treat the twist by both macroscopic and microscopic ways. In the latter case, the shell structure of clusters is fully exploited, which is crucial for the considered size region ( 8