Source author record

P. -G. Reinhard

P. -G. Reinhard appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

98works
11topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

98 published item(s)

preprint2022arXiv

Electric and Magnetic Moments and Transition Probabilities in $^{208}$Pb $\pm 1$ Nuclei

We present moments and transition probabilities in the neighboring odd-mass nuclei of $^{208}$Pb calculated fully self-consistently from the s.p. properties of $^{208}$Pb with polarization corrections from its excitations, both given from previous Skyrme-Hartree-Fock and RPA calculations. The electric results agree nicely with the data with two very interesting exceptions. In the magnetic case we obtain similar results. We discuss also polarization contributions to the $l$-forbidden $M1$ transitions, which are, however, much too small compared to the data. With a modified external field operator which accounts effectively for mesonic and many-body effects the description of the data can be substantially improved.

preprint2021arXiv

Fine Structure of the Isovector Giant Dipole Resonance in $^{142-150}$Nd and $^{152}$Sm

Background: Inelastic proton scattering at energies of a few hundred MeV and very-forward angles including $0^\circ$ has been established as a tool to study electric-dipole strength distributions in nuclei. The present work reports a systematic investigation of the chain of stable even-mass Nd isotopes representing a transition from spherical to quadrupole-deformed nuclei. Purpose: Extraction of the equivalent photo-absorption cross sections and analysis of their fine structure in the energy region of the IsoVector Giant Dipole Resonance (IVGDR). Method: Proton inelastic scattering reactions of 200 MeV protons were measured at iThemba LABS in Cape Town, South Africa. The scattering products were momentum-analysed by the K600 magnetic spectrometer positioned at $θ_{\mathrm{Lab}}=0^\circ$. Using dispersion-matching techniques, energy resolutions of $ΔE \approx 40 - 50$ keV were obtained. After subtraction of background and contributions from other multipoles, the spectra were converted to photo-absorption cross sections using the equivalent virtual-photon method. Results: Wavelet-analysis techniques are used to extract characteristic energy scales of the fine structure of the IVGDR from the experimental data. Comparisons with the Quasiparticle-Phonon Model (QPM) and Skyrme Separable Random Phase Approximation (SSRPA) predictions provide insight into the role of different giant resonance damping mechanisms. Conclusions: Fine structure is observed even for the most deformed nuclei studied. Fragmentation of the one particle-one hole ($1p1h$) strength seems to be the main source of fine structure in both spherical and deformed nuclei. Some impact of the spreading due to coupling of the two particle-two hole ($2p2h$) states to the $1p1h$ doorway states is seen in the spherical/transitional nuclei, where calculations beyond the $1p1h$ level are available.

preprint2021arXiv

Nuclear Charge Radii of the Nickel Isotopes $^{58-68,70}$Ni

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$_{\rm sat}$, which allows an assessment of their accuracy. We find agreement with experiment in differential radii $δ\left\langle r_\mathrm{c}^2 \right\rangle$ for all employed ab initio methods and interactions, while the absolute radii are consistent with data only for NNLO$_{\rm sat}$. Within nuclear density functional theory, the Skyrme functional SV-min matches experiment more closely than the Fayans functional Fy($Δr$,HFB).

preprint2021arXiv

Self-consistent description of high-spin states in doubly magic $^{208}$Pb

We analyze recent data on a long series of high-spin states in $^{208}$Pb with a self-consistent phonon-coupling model for nuclear excitations based on the Skyrme functionals. The model is the renormalized time-blocking approximation (RenTBA) which takes the coherent one-particle-one-hole (1p1h) states of the random-phase approximation (RPA) as starting point and develops from that more complex configurations beyond RPA. To the best of our knowledge, this is the first investigation of high spin states in $^{208}$Pb using self-consistent nuclear models. The interesting point here is that complex configurations are compulsory to describe the upper end of the long spin series at all. The data thus provide an ideal testing ground for phonon-coupling models as they give direct access to complex configurations. We find that standard Skyrme functionals which perform well in ground state properties and giant resonance excitations deliver at once an agreeable description of these high spin states.

preprint2020arXiv

M1 resonance in $^{208}$Pb within the self-consistent phonon-coupling model

The main goal of the paper is to investigate theoretically the experimentally observed fragmentation of the isovector $M1$ resonance in $^{208}$Pb within a self-consistent model based on an energy-density functional (EDF) of the Skyrme type. This fragmentation (spread of the $M1$ strength) is not reproduced in a conventional one-particle--one-hole ($1p1h$) random-phase approximation (RPA) and thus has to be investigated in the framework of more complicated models. However, previously applied models of this type were not self-consistent. In the present work, we use a recently developed renormalized version of the self-consistent time blocking approximation (RenTBA) in which the $1p1h\otimes$phonon configurations are included on top of the RPA $1p1h$ configurations. We have determined several sets of the parameters of the modified Skyrme EDF fitted within the RenTBA and RPA and have found the necessary condition of producing the fragmentation of the $M1$ resonance in $^{208}$Pb in our model. We present also the results of the RenTBA and RPA calculations for the first excited states of the natural parity modes in $^{208}$Pb obtained with these modified parametrizations.

preprint2020arXiv

Measurement and microscopic description of odd-even staggering of charge radii of exotic copper isotopes

The mesoscopic nature of the atomic nucleus gives rise to a wide array of macroscopic and microscopic phenomena. The size of the nucleus is a window into this duality: while the charge radii globally scale as $A^{1/3}$, their evolution across isotopic chains reveals unanticipated structural phenomena [1-3]. The most ubiquitous of these is perhaps the Odd-Even Staggering (OES) [4]: isotopes with an odd number of neutrons are usually smaller in size than the trend of their even-neutron neighbours suggests. This OES effect varies with the number of protons and neutrons and poses a significant challenge for nuclear theory [5-7]. Here, we examine this problem with new measurements of the charge radii of short-lived copper isotopes up to the very exotic $^{78}$Cu $(Z=29, N=49)$, produced at only 20 ions/s, using the highly-sensitive Collinear Resonance Ionisation Spectroscopy (CRIS) method at ISOLDE-CERN. Due to the presence of a single proton outside of the closed Z=28 shell, these measurements provide crucial insights into the single-particle proton structure and how this affects the charge radii. We observe an unexpected reduction in the OES for isotopes approaching the $N=50$ shell gap. To describe the data, we applied models based on nuclear Density Functional Theory [2,8] (DFT) and ab-initio Valence-Space In-Medium Similarity Renormalization Group (VS-IMSRG) theory [9,10]. Through these comparisons, we demonstrate a relation between the global behavior of charge radii and the saturation density of nuclear matter, and show that the local charge radii variations, which reflect the many-body polarization effects due to the odd neutron, naturally emerge from the VS-IMSRG calculations.

preprint2020arXiv

Skyrme RPA for nuclear resonances: trouble with magnetic modes

We discuss major differences between electric and magnetic excitations in nuclei appearing in self-consistent calculation based on Skyrme energy-density functionals. Tools of analysis are Landau-Migdal parameters for bulk properties and RPA for resonance modes of $^{208}$Pb as representative of finite nuclei. We show that the relation between the effective mass and the effective particle-hole interaction, well known in the Landau-Migdal theory, explains the success of self-consistent calculations of electric transitions in such approaches. This effect, however, does not automatically exist in the magnetic case. This calls for further developments of the Skyrme functional in the spin channel.

preprint2019arXiv

A Survey of Nuclear Pasta in the Intermediate Density Regime: Structure Functions for Neutrino Scattering

Background: Nuclear pasta matter, emerging due to the competition between the long-range Coulomb force and the short-range strong force, is believed to be present in astrophysical scenarios, such as neutron stars and core-collapse supernovae. Its structure can have a high impact e.g. on neutrino transport or the tidal deformability of neutron stars. Purpose: We investigate the impact of nuclear pasta on neutrino interactions and compare the results to uniform matter. Method: We calculate the elastic and inelastic static structure factors for nuclear pasta matter using density functional theory (DFT), which contain the main nuclear input for neutrino scattering. Results: Each pasta structure leaves a unique imprint in the elastic structure factor and it is largely enhanced. The inelastic structure factors are very similar for all configurations. Conclusion: Nuclear pasta has a noticeable impact on neutrino neutral-current scattering opacities. While for inelastic reactions the cross section is reduced, the elastic coherent scattering increases dramatically. The effect can be of importance for the cooling of neutron stars as well as for core-collapse supernova models.

preprint2019arXiv

Beyond the charge radius: the information content of the fourth radial moment

Measurements of atomic transitions in different isotopes offer key information on the nuclear charge radius. The anticipated high-precision experimental techniques, augmented by atomic calculations, will soon enable extraction of the higher-order radial moments of the charge density distribution. To assess the value of such measurements for nuclear structure research, we study the information content of the fourth radial moment $\rf$ by means of nuclear density functional theory and a multiple correlation analysis. We show that $\rf$ can be directly related to the surface thickness of nuclear density, a fundamental property of the atomic nucleus that is difficult to obtain for radioactive systems. Precise knowledge of these radial moments is essential to establish reliable constraints on the existence of new forces from precision isotope shift measurements .

preprint2019arXiv

Systematics of toroidal dipole modes in Ca, Ni, Zr, and Sn isotopes

We analyze the relation between isoscalar toroidal modes and so-called pygmy dipole resonance (PDR) which both appear in the same region of low-energy dipole excitations. To this end, we use a theoretical description within the fully self-consistent Skyrme quasiparticle random-phase approximation (QRPA). Test cases are spherical nuclei $^{40,48}$Ca, $^{58,72}$Ni, $^{90,100}$Zr, and $^{100,120,132}$Sn which cover four different elements and for each element at least two isotopes with different neutron excess, one small and another large. The structure of the modes is investigated in terms of strength functions, transition densities (TD) and transition currents (TC). For all considered nuclei, we see that, independently on whether PDR strength exists or not, the flow pattern in the lower part of the "PDR energy region" is basically the isoscalar vortical toroidal motion with a minor irrotational fraction. A one-to-one correspondence between calculated TD and TC can be established. The toroidal flow appears already in the uncoupled two-quasiparticle (2qp) excitations and becomes definitively strong for the QRPA modes. Altogether, we find that low-lying dipole strength often denoted as isoscalar PDR is actually an oversimplified imitation of the basically toroidal motion in nuclei with a sufficient neutron excess.

preprint2018arXiv

Individual low-energy E1 toroidal and compression states in light nuclei: deformation effect, spectroscopy and interpretation

The existence of individual low-energy E1 toroidal and compression states (TS and CS) in $^{24}$Mg was predicted recently in the framework of quasiparticle random-phase-approximation (QRPA) model with Skyrme forces. It was shown that the strong axial deformation of $^{24}$Mg is crucial to downshift the toroidal strength to the low-energy region and thus make the TS the lowest E1(K=1) dipole state. In this study, we explain this result by simple mean-field arguments. Comparing TS in two strongly axial nuclei, $^{24}$Mg and $^{20}$Ne, we show that the lowest TS is not not a universal phenomenon but rather a peculiarity of $^{24}$Mg. The spectroscopy of TS and CS is analyzed and some additional interpretation of these states is suggested.

preprint2016arXiv

Application of an Extended Random Phase Approximation on Giant Resonances in Light, Medium and Heavy Mass Nuclei

We present results of the time blocking approximation (TBA) on giant resonances in light, medium and heavy mass nuclei. The TBA is an extension of the widely used random-phase approximation (RPA) adding complex configurations by coupling to phonon excitations. A new method for handling the single-particle continuum is developed and applied in the present calculations. We investigate in detail the dependence of the numerical results on the size of the single particle space and the number of phonons as well as on nuclear matter properties. Our approach is self-consistent, based on an energy-density functional of Skyrme type where we used seven different parameter sets. The numerical results are compared with experimental data.

preprint2016arXiv

Carbon-oxygen-neon mass nuclei in superstrong magnetic fields

The properties of $\isotope[12]{C}$, $\isotope[16]{O}$, and $\isotope[20]{Ne}$ nuclei in strong magnetic fields $B\simeq 10^{17}\,$G are studied in the context of strongly magnetized neutron stars and white dwarfs. The SKY3D code is extended to incorporate the interaction of nucleons with the magnetic field and is utilized to solve the time-independent Hartree-Fock equations with a Skyrme interaction on a Cartesian three-dimensional grid. The numerical solutions demonstrate a number of phenomena, which include a splitting of the energy levels of spin-up and -down nucleons, spontaneous rearrangment of energy levels in $\isotope[16]{O}$ at a critical field, which leads to jump-like increases of magnetization and proton current in this nucleus, and evolution of the intrinsically deformed $\isotope[20]{Ne}$ nucleus towards a more spherical shape under increasing field strength. Many of the numerical features can be understood within a simple analytical model based on the occupation by the nucleons of the lowest states of the harmonic oscillator in a magnetic field.

preprint2016arXiv

Magnetic dipole excitations of $^{50}$Cr

The low-lying $M1$-strength of the open-shell nucleus $^{50}$Cr has been studied with the method of nuclear resonance fluorescence up to 9.7 MeV, using bremsstrahlung at the superconducting Darmstadt linear electron accelerator S-DALINAC and Compton backscattered photons at the High Intensity $γ$-ray Source (HI$γ$S) facility between 6 and 9.7 MeV of the initial photon energy. Fifteen $1^{+}$ states have been observed between 3.6 and 9.7 MeV. Following our analysis, the lowest $1^{+}$ state at 3.6 MeV can be considered as an isovector orbital mode with some spin admixture. The obtained results generally match the estimations and trends typical for the scissors-like mode. Detailed calculations within the Skyrme Quasiparticle Random-Phase-Approximation method and the Large-Scale Shell Model justify our conclusions. The calculated distributions of the orbital current for the lowest $1^{+}$-state suggest the schematic view of Lipparini and Stringari (isovector rotation-like oscillations inside the rigid surface) rather than the scissors-like picture of Lo Iudice and Palumbo. The spin M1 resonance is shown to be mainly generated by spin-flip transitions between the orbitals of the $fp$-shell.

preprint2016arXiv

Nuclear charge and neutron radii and nuclear matter: trend analysis

Radii of charge and neutron distributions are fundamental nuclear properties. They depend on both nuclear interaction parameters related to the equation of state of infinite nuclear matter and on quantal shell effects, which are strongly impacted by the presence of nuclear surface. In this work, by studying the dependence of charge and neutron radii, and neutron skin, on nuclear matter parameters, we assess different mechanisms that drive nuclear sizes. We apply nuclear density functional theory using a family of Skyrme functionals obtained by means of different optimization protocols targeting specific nuclear properties. By performing the Monte-Carlo sampling of reasonable functionals around the optimal parametrization, we study correlations between nuclear matter paramaters and observables characterizing charge and neutron distributions. We demonstrate the existence of the strong converse relation between the nuclear charge radii and the saturation density of symmetric nuclear matter and also between the neutron skins and the slope of the symmetry energy. For functionals optimized to experimental binding energies only, proton and neutron radii are weakly correlated due to canceling trends from different nuclear matter parameters. We show that by requiring that the nuclear functional reproduces the empirical saturation point of symmetric nuclear matter practically fixes the charge (or proton) radii, and vice versa. The neutron skin uncertainty primarily depends on the slope of the symmetry energy. Consequently, imposing a constraint on both $ρ_0$ and $L$ practically determines the nuclear size, modulo small variations due to shell effects.

preprint2016arXiv

Sensitivity of the fusion cross section to the density dependence of the symmetry energy

It is the aim of this paper to discuss the impact of nuclear fusion on the EOS. This is a timely subject given the expected availability of increasingly exotic beams at rare isotope facilities\,\cite{balantekin2014}. In practice, we focus on $^{48}$Ca+$^{48}$Ca fusion. We employ three different approaches to calculate fusion cross-sections for a set of energy density functionals with systematically varying nuclear matter properties. Fusion calculations are performed using frozen densities, using a dynamic microscopic method based on density-constrained time-dependent Hartree-Fock (DC-TDHF) approach, as well as direct TDHF study of above barrier cross-sections. For these studies, we employ a family of Skyrme parametrizations with systematically varied nuclear matter properties. We find a slight preference for forces which deliver a slope of symmetry energy of $L\approx 50$\,MeV that corresponds to a neutron-skin thickness of $^{48}$Ca of $R_\mathrm{skin}\!=\!(0.180\!-\!0.210)$\,fm.

preprint2016arXiv

Time-dependent density functional theory with twist-averaged boundary conditions

Time-dependent density functional theory is widely used to describe excitations of many-fermion systems. In its many applications, 3D coordinate-space representation is used, and infinite-domain calculations are limited to a finite volume represented by a box. For finite quantum systems (atoms, molecules, nuclei), the commonly used periodic or reflecting boundary conditions introduce spurious quantization of the continuum states and artificial reflections from boundary; hence, an incorrect treatment of evaporated particles. These artifacts can be practically cured by introducing absorbing boundary conditions (ABC) through an absorbing potential in a certain boundary region sufficiently far from the described system. But also the calculations of infinite matter (crystal electrons, quantum fluids, neutron star crust) suffer artifacts from a finite computational box. In this regime, twist- averaged boundary conditions (TABC) have been used successfully to diminish the finite-volume effects. In this work, we extend TABC to time-dependent framework and apply it to resolve the box artifacts for finite quantum systems using as test case small- and large-amplitude nuclear vibrations. We demonstrate that by using such a method, one can reduce finite volume effects drastically without adding any additional parameters. While they are almost equivalent in the linear regime, TABC and ABC differ in the nonlinear regime in their treatment of evaporated particles.

preprint2015arXiv

Deformation effects in Giant Monopole Resonance

The isoscalar giant monopole resonance (GMR) in Samarium isotopes (from spherical $^{144}$Sm to deformed $^{148-154}$Sm) is investigated within the Skyrme random-phase-approximation (RPA) for a variety of Skyrme forces. The exact RPA and its separable version (SRPA) are used for spherical and deformed nuclei, respectively. The quadrupole deformation is shown to yield two effects: the GMR broadens and attains a two-peak structure due to the coupling with the quadrupole giant resonance.

preprint2015arXiv

Deformation-induced splitting of the monopole giant resonance in 24Mg

The strong deformation splitting of the isoscalar giant monopole resonance (ISGMR), recently observed in ($α,α'$) reaction in prolate $^{24}$Mg, is analyzed in the framework of the Skyrme quasiparticle random-phase-approximation (QRPA) approach with the Skyrme forces SkM*, SVbas and SkP$^δ$. The calculations with these forces give close results and confirm that the low-energy E0-peak is caused by the deformation-induced coupling of ISGMR with the K=0 branch of the isoscalar giant quadrupole resonance.

preprint2015arXiv

Dipole polarizability of 120Sn and nuclear energy density functionals

The electric dipole strength distribution in 120Sn between 5 and 22 MeV has been determined at RCNP Osaka from a polarization transfer analysis of proton inelastic scattering at E_0 = 295 MeV and forward angles including 0°. Combined with photoabsorption data an electric dipole polarizability α_D(120Sn) = 8.93(36) fm^3 is extracted. The dipole polarizability as isovector observable par excellence carries direct information on the nuclear symmetry energy and its density dependence. The correlation of the new value with the well established α_D(208Pb) serves as a test of its prediction by nuclear energy density functionals (EDFs). Models based on modern Skyrme interactions describe the data fairly well while most calculations based on relativistic Hamiltonians cannot.

preprint2015arXiv

Estimating the relevance of predictions from nuclear mean-field models

This contribution reviews the present status of the Skyrme-Hartree-Fock (SHF) approach as one of the leading self-consistent mean-field models in the physics of atomic nuclei. It starts with a brief summary of the formalism and strategy for proper calibration of the SHF functional. The main emphasis lies on an exploration of the reliability of predictions, particularly in the regime of extrapolations. Various strategies are discussed to explore the statistical and systematic errors of SHF. The strategies are illustrated on examples from actual applications. Variations of model and fit data are used to get an idea about systematic errors. The statistical error is evaluated in straightforward manner by statistical analysis based on $χ^2$ fits. This also allows also to evaluate the correlations (covariances) between observables which provides useful insights into the structure of the model and of the fitting strategy.

preprint2015arXiv

Fission barriers and probabilities of spontaneous fission for elements with Z$\geq$100

This is a short review of methods and results of calculations of fission barriers and fission half-lives of even-even superheavy nuclei. An approvable agreement of the following approaches is shown and discussed: The macroscopic-microscopic approach based on the stratagem of the shell correction to the liquid drop model and a vantage point of microscopic energy density functionals of Skyrme and Gogny type selfconsistently calculated within Hartree-Fock-Bogoliubov method. Mass parameters are calculated in the Hartree-Fock-Bogoliubov cranking approximation. A short part of the paper is devoted to the nuclear fission dynamics. We also discuss the predictive power of Skyrme functionals applied to key properties of the fission path of $^{266}$Hs. It applies the standard techniques of error estimates in the framework of a $χ^2$ analysis.

preprint2015arXiv

Monopole giant resonance in $^{100-132}$Sn, $^{144}$Sm and $^{208}$Pb

The isoscalar giant monopole resonance (GMR) in spherical nuclei $^{100-132}$Sn, $^{144}$Sm, and $^{208}$Pb is investigated within the Skyrme random-phase-approximation (RPA) for a variety of Skyrme forces and different pairing options. The calculated GMR strength functions are directly compared to the available experimental distributions. It is shown that, in accordance to results of other groups, description of GMR in Sn and heavier Sm/Pb nuclei needs different values of the nuclear incompressibilty, $K \approx$ 200 or 230 MeV, respectively. Thus none from the used Skyrme forces is able to describe GMR in these nuclei simultaneously. The GMR peak energy in open-shell $^{120}$Sn is found to depend on the isoscalar effective mass, which might be partly used for a solution of the above problem. Some important aspects of the problem (discrepancies of available experimental data, proper treatment of the volume and surface compression in finite nuclei, etc) are briefly discussed.

preprint2015arXiv

On the analysis of photo-electron spectra

We analyze Photo-Electron Spectra (PES) for a variety of excitation mechanisms from a simple mono-frequency laser pulse to involved combination of pulses as used, e.g., in attosecond experiments. In the case of simple pulses, the peaks in PES re- flect the occupied single-particle levels in combination with the given laser frequency. This usual, simple rule may badly fail in the case of excitation pulses with mixed frequencies and if resonant modes of the system are significantly excited. We thus develop an extension of the usual rule to cover all possible excitation scenarios, including mixed frequencies in the attosecond regime. We find that the spectral dis- tributions of dipole, monopole and quadrupole power for the given excitation taken together and properly shifted by the single-particle energies provide a pertinent picture of the PES in all situations. This leads to the derivation of a generalized relation allowing to understand photo-electron yields even in complex experimental setups.

preprint2015arXiv

Swelling of nuclei embedded in neutron-gas and consequences for fusion

Fusion of very neutron rich nuclei may be important to determine the composition and heating of the crust of accreting neutron stars. We present an exploratory study of the effect of the neutron-gas environment on the structure of nuclei and the consequences for pycnonuclear fusion cross-sections in the neutron drip region. We studied the formation and properties of Oxygen and Calcium isotopes embedded in varying neutron-gas densities. We observe that the formed isotope is the drip-line nucleus for the given effective interaction. Increasing the neutron-gas density leads to the swelling of the nuclear density. We have used these densities to study the effect of this swelling on the fusion cross-sections using the São-Paulo potential. At high neutron-gas densities the cross-section is substantially increased but at lower densities the modification is minimal.

preprint2014arXiv

A quantum relaxation-time approximation for finite fermion systems

We propose a relaxation time approximation for the description of the dynamics of strongly excited fermion systems. Our approach is based on time-dependent density functional theory at the level of the local density approximation. This mean-field picture is augmented by collisional correlations handled in relaxation time approximation which is inspired from the corresponding semi-classical picture. The method involves the estimate of microscopic relaxation rates/times which is presently taken from the well established semi-classical experience. The relaxation time approximation implies evaluation of the instantaneous equilibrium state towards which the dynamical state is progressively driven at the pace of the microscopic relaxation time. As test case, we consider Na clusters of various sizes excited either by a swift ion projectile or by a short and intense laser pulse, driven in various dynamical regimes ranging from linear to strongly non-linear reactions. We observe a strong effect of dissipation on sensitive observables such as net ionization and angular distributions of emitted electrons. The effect is especially large for moderate excitations where typical relaxation/dissipation time scales efficiently compete with ionization for dissipating the available excitation energy. Technical details on the actual procedure to implement a working recipe of such a quantum relaxation approximation are given in appendices for completeness.

preprint2014arXiv

Appearance of the Single Gyroid Network Phase in Nuclear Pasta Matter

Nuclear matter under the conditions of a supernova explosion unfolds into a rich variety of spatially structured phases, called nuclear pasta. We investigate the role of periodic network-like structures with negatively curved interfaces in nuclear pasta structures, by static and dynamic Hartree-Fock simulations in periodic lattices. As the most prominent result, we identify for the first time the {\it single gyroid} network structure of cubic chiral $I4_123$ symmetry, a well known configuration in nanostructured soft-matter systems, both as a dynamical state and as a cooled static solution. Single gyroid structures form spontaneously in the course of the dynamical simulations. Most of them are isomeric states. The very small energy differences to the ground state indicate its relevance for structures in nuclear pasta.

preprint2014arXiv

Electrons as probes of dynamics in molecules and clusters : a contribution from Time Dependent Density Functional Theory

Various ways to analyze the dynamical response of clusters and molecules to electromagnetic perturbations exist. Particularly rich information can be obtained from measuring the properties of electrons emitted in the course of the excitation dynamics. Such an analysis of electron signals covers total ionization, Photo-Electron Spectra, Photoelectron Angular Distributions, and ideally combined PES/PAD, with a long history in molecular physics, also increasingly used in cluster physics. Recent progress in the design of new light sources (high intensity and/or frequency, ultra short pulses) opens new possibilities for measurements and thus has renewed the interest on the analysis of dynamical scenarios through these observables, well beyond a simple access to a density of states. This, in turn, has motivated many theoretical investigations of the dynamics of electronic emission for molecules and clusters. A theoretical tool of choice is here Time-Dependent Density Functional Theory (TDDFT) propagated in real time and on a spatial grid, and augmented by a Self-Interaction Correction. This provides a pertinent, robust, and efficient description of electronic emission including the detailed pattern of PES and PAD. A direct comparison between experiments and well founded elaborate microscopic theories is thus readily possible, at variance with more demanding observables such as for example fragmentation or dissociation cross sections. The aim of this paper is to review the available experimental results motivating such studies, describe the theoretical tools developed on the basis of real-time and real-space TDDFT to address in a realistic manner the analysis of electronic emission following irradiation of clusters and molecules by various laser pulses, discuss representative results, and finally give some future directions of investigations.

preprint2014arXiv

Error estimates for the Skyrme-Hartree-Fock model

There are many complementing strategies to estimate the extrapolation errors of a model which was calibrated in least-squares fits. We consider the Skyrme-Hartree-Fock model for nuclear structure and dynamics and exemplify the following five strategies: uncertainties from statistical analysis, covariances between observables, trends of residuals, variation of fit data, dedicated variation of model parameters. This gives useful insight into the impact of the key fit data as they are: binding energies, charge r.m.s. radii, and charge formfactor. Amongst others, we check in particular the predictive value for observables in the stable nucleus $^{208}$Pb, the super-heavy element $^{266}$Hs, $r$-process nuclei, and neutron stars.

preprint2014arXiv

Error Estimates of Theoretical Models: a Guide

This guide offers suggestions/insights on uncertainty quantification of nuclear structure models. We discuss a simple approach to statistical error estimates, strategies to assess systematic errors, and show how to uncover inter-dependencies by correlation analysis. The basic concepts are illustrated through simple examples. By providing theoretical error bars on predicted quantities and using statistical methods to study correlations between observables, theory can significantly enhance the feedback between experiment and nuclear modeling.

preprint2014arXiv

Landau-Migdal vs. Skyrme

The magnitude and density-dependence of the non-spin dependent Landau-Migdal parameters are derived from Skyrme energy functionals and compared with the phenomenological ones. We perform RPA calculations with various approximations for the Landau-Migdal particle-hole interaction and compare them with the results obtained with the full Skyrme interaction. For the first time the next to leading order in the Landau-Migdal approach is considered in nuclear structure calculations.

preprint2014arXiv

Multi-plasmon excitations in electron spectra of small systems irradiated by swift charged projectiles

We investigate the kinetic-energy spectrum of electrons emitted from an excited many-electron system, often called photo-electron spectrum (PES). We are particularly interested on the impact of resonant modes of the system on PES. To this end, we consider three systems with strong resonances, a Mg atom, the small alkaline cluster ${{\rm K}_9}^+$, and the small carbon chain C$_3$. To avoid dominant frequencies in the excitation process, we consider a collision with a fast ion which is realized by an instantaneous boost of the valence electrons, a process which excites all frequencies with equal weight. The electron dynamics is investigated from a theoretical perspective using time-dependent density-functional theory augmented by an average-density self-interaction correction. We observe patterns which are similar to PES usually obtained after irradiation by a laser pulse, in particular the appearance of clear peaks. We show that these patterns are driven by strong resonance modes of the system. Resonances are thus found to be another source of peaks in the PES, besides photons (when present) with definite frequencies.

preprint2014arXiv

Nuclear energy density optimization: Shell structure

Nuclear density functional theory is the only microscopical theory that can be applied throughout the entire nuclear landscape. Its key ingredient is the energy density functional. In this work, we propose a new parameterization UNEDF2 of the Skyrme energy density functional. The functional optimization is carried out using the POUNDerS optimization algorithm within the framework of the Skyrme Hartree-Fock-Bogoliubov theory. Compared to the previous parameterization UNEDF1, restrictions on the tensor term of the energy density have been lifted, yielding a very general form of the energy density functional up to second order in derivatives of the one-body density matrix. In order to impose constraints on all the parameters of the functional, selected data on single-particle splittings in spherical doubly-magic nuclei have been included into the experimental dataset. The agreement with both bulk and spectroscopic nuclear properties achieved by the resulting UNEDF2 parameterization is comparable with UNEDF1. While there is a small improvement on single-particle spectra and binding energies of closed shell nuclei, the reproduction of fission barriers and fission isomer excitation energies has degraded. As compared to previous UNEDF parameterizations, the parameter confidence interval for UNEDF2 is narrower. In particular, our results overlap well with those obtained in previous systematic studies of the spin-orbit and tensor terms. UNEDF2 can be viewed as an all-around Skyrme EDF that performs reasonably well for both global nuclear properties and shell structure. However, after adding new data aiming to better constrain the nuclear functional, its quality has improved only marginally. These results suggest that the standard Skyrme energy density has reached its limits and significant changes to the form of the functional are needed.

preprint2014arXiv

Nuclear Energy Density Optimization: UNEDF2

The parameters of the UNEDF2 nuclear energy density functional (EDF) model were obtained in an optimization to experimental data consisting of nuclear binding energies, proton radii, odd-even mass staggering data, fission-isomer excitation energies, and single particle energies. In addition to parameter optimization, sensitivity analysis was done to obtain parameter uncertainties and correlations. The resulting UNEDF2 is an all-around EDF. However, the sensitivity analysis also demonstrated that the limits of current Skyrme-like EDFs have been reached and that novel approaches are called for.

preprint2014arXiv

Nuclear Pasta Matter for Different Proton Fractions

Nuclear matter under astrophysical conditions is explored with time-dependent and static Hartree-Fock calculations. The focus is in a regime of densities where matter segregates into liquid and gaseous phases unfolding a rich scenario of geometries, often called nuclear pasta shapes (e.g. spaghetti, lasagna). Particularly the appearance of the different phases depending on the proton fraction and the transition to uniform matter are investigated. In this context the neutron background density is of special interest, because it plays a crucial role for the type of pasta shape which is built. The study is performed in two dynamical ranges, one for hot matter and one at temperature zero to investigate the effect of cooling.

preprint2014arXiv

On the dynamics of photo-electrons from C$_{60}$

We explore photo-electron spectra (PES) and photo-electron angular distributions (PAD) of C$_{60}$ with time-dependent density functional theory (TDDFT) in real time. To simulate experiments in gas phase, we consider isotropic ensembles of cluster orientations and perform orientation averaging of the TDDFT calculations. First, we investigate ionization properties of C$_{60}$ by one-photon processes in the range of VUV energies. The PES map the energies of the occupied single-particle states, while the weights of the peaks in PES are given by the depletion of the corresponding level. The different influences can be disentangled by looking at PES from slightly different photon frequencies. PAD in the one-photon regime can be characterized by one parameter, the anisotropy. This single parameter unfolds worthwhile information when investigating the frequency and state dependences. We also discuss the case of multi-photon ionization induced by strong infrared laser pulses in C$_{60}$. In agreement with measurements, we find that the PES show a regular comb of peaks separated by the photon energy. Our calculations reveal that this happens because only very few occupied states of C$_{60}$ near the ionization threshold contribute to emission and that these few states happen to cooperate filling the same peaks. The PAD show a steady increase of anisotropy with increasing photon order.

preprint2014arXiv

Relation of E1 pygmy and toroidal resonances

A possible relation of the low-lying E1 (pygmy resonance) and toroidal strengths is analyzed by using Skyrme-RPA results for the strength functions, transition densities and current fields in $^{208}$Pb. It is shown that the irrotational pygmy motion can appear as a local manifestation of the collective vortical toroidal dipole resonance (TDR) at the nuclear surface. The RPA results are compared to unperturbed (1ph) ones.

preprint2013arXiv

Deformation effects in toroidal and compression dipole excitations of $^{170}$Yb: Skyrme-RPA analysis

The effect of nuclear deformation on the isoscalar toroidal and compression dipole modes in prolate $^{170}$Yb is studied in the framework of the random-phase-approximation method with a representative set of Skyrme forces (SV-bas, SLy6, SkM$^*$ and SkI3). It is shown that the deformation crucially redistributes the strength of both modes. The compression mode has the same sequence of $μ$=0 and 1 branches as the isovector giant dipole resonance where for prolate nuclei the $μ=0$ mode is lower in energy ($μ$ being the projection of the axial momentum of the mode). Instead, the toroidal mode exhibits an anomalous (opposite) sequence where the $μ$=1 branch precedes the $μ$=0 one.

preprint2013arXiv

Information content of the weak-charge form factor

Parity-violating electron scattering provides a model-independent determination of the nuclear weak-charge form factor that has widespread implications across such diverse areas as fundamental symmetries, nuclear structure, heavy-ion collisions, and neutron-star structure. We assess the impact of precise measurements of the weak-charge form factor of ${}^{48}$Ca and ${}^{208}$Pb on a variety of nuclear observables, such as the neutron skin and the electric-dipole polarizability. We use the nuclear Density Functional Theory with several accurately calibrated non-relativistic and relativistic energy density functionals. To assess the degree of correlation between nuclear observables and to explore systematic and statistical uncertainties on theoretical predictions, we employ the chi-square statistical covariance technique. We find a strong correlation between the weak-charge form factor and the neutron radius, that allows for an accurate determination of the neutron skin of neutron-rich nuclei. We determine the optimal range of the momentum transfer $q$ that maximizes the information content of the measured weak-charge form factor and quantify the uncertainties associated with the strange quark contribution. Moreover, we confirm the role of the electric-dipole polarizability as a strong isovector indicator. Accurate measurements of the weak-charge form factor of ${}^{48}$Ca and ${}^{208}$Pb will have a profound impact on many aspects of nuclear theory and hadronic measurements of neutron skins of exotic nuclei at radioactive-beam facilities.

preprint2013arXiv

Nuclear vorticity in isoscalar E1 modes: Skyrme-RPA analysis

Two basic concepts of nuclear vorticity, hydrodynamical (HD) and Rawenthall-Wambach (RW), are critically inspected. As a test case, we consider the interplay of irrotational and vortical motion in isoscalar electric dipole E1(T=0) modes in $^{208}$Pb, namely the toroidal and compression modes. The modes are described in a self-consistent random-phase-approximation (RPA) with the Skyrme force SLy6. They are examined in terms of strength functions, transition densities, current fields, and formfactors. It is shown that the RW conception (suggesting the upper component of the nuclear current as the vorticity indicator) is not robust. The HD vorticity is not easily applicable either because the definition of a velocity field is too involved in nuclear systems. Instead, the vorticity is better characterized by the toroidal strength which closely corresponds to HD treatment and is approximately decoupled from the continuity equation.

preprint2013arXiv

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.

preprint2013arXiv

Symmetry energy in nuclear density functional theory

The nuclear symmetry energy represents a response to the neutron-proton asymmetry. In this survey we discuss various aspects of symmetry energy in the framework of nuclear density functional theory, considering both non-relativistic and relativistic self-consistent mean-field realizations side-by-side. Key observables pertaining to bulk nucleonic matter and finite nuclei are reviewed. Constraints on the symmetry energy and correlations between observables and symmetry-energy parameters, using statistical covariance analysis, are investigated. Perspectives for future work are outlined in the context of ongoing experimental efforts.

preprint2013arXiv

The TDHF Code Sky3D

The nuclear mean-field model based on Skyrme forces or related density functionals has found wide-spread application to the description of nuclear ground states, collective vibrational excitations, and heavy-ion collisions. The code Sky3D solves the static or dynamic equations in a three-dimensional Cartesian mesh with isolated or periodic boundary conditions and no further symmetry assumptions. Pairing can be included in the BCS approximation. The code is implemented with a view to allow easy modifications for including additional physics or special analysis of the results.

preprint2013arXiv

Toroidal, compression, and vortical dipole strengths in $^{144-154}$Sm: Skyrme-RPA exploration of deformation effect

A comparative analysis of toroidal, compressional and vortical dipole strengths in the spherical $^{144}$Sm and the deformed $^{154}$Sm is performed within the random-phase-approximation using a set of different Skyrme forces. Isoscalar (T=0), isovector (T=1), and electromagnetic excitation channels are considered. The role of the nuclear convection $j_{\text{con}}$ and magnetization $j_{\text{mag}}$ currents is inspected. It is shown that the deformation leads to an appreciable redistribution of the strengths and causes a spectacular deformation splitting (exceeding 5 MeV) of the isoscalar compressional mode. In $^{154}$Sm, the $μ$=0 and $μ$=1 branches of the mode form well separated resonances. When stepping from $^{144}$Sm to $^{154}$Sm, we observe an increase of the toroidal, compression and vortical contributions in the low-energy region (often called pygmy resonance). The strength in this region seems to be an overlap of various excitation modes. The energy centroids of the strengths depend significantly on the isoscalar effective mass $m_0$. Skyrme forces with a large $m_0$ (typically $m_0/m \approx 0.8 - 1$) seem to be more suitable for description of experimental data for the isoscalar giant dipole resonance.

preprint2012arXiv

A critical analysis of the theoretical scheme to evaluate photoelectron spectra

We discuss in depth the validity and limitations of a theoretical scheme to evaluate photo-electron spectra (PES) through collecting the phase oscillations at a given measuring point. Problems appear if the laser pulse is still active when the first bunches of outgoing flow reach the measuring point. This limits the simple scheme for evaluation of PES to low and moderate laser intensities. Using a model of free particle plus dipole field, we develop a generalized scheme which is shown to considerably improve the results for high intensities.

preprint2012arXiv

Dynamic microscopic study of pre-equilibrium giant resonance excitation and fusion in the reactions $^{132}$Sn+$^{48}$Ca and $^{124}$Sn+$^{40}$Ca

We study pre-equilibrium giant dipole resonance excitation and fusion in the neutron-rich system $^{132}$Sn+$^{48}$Ca at energies near the Coulomb barrier, and we compare photon yields and total fusion cross sections to those of the stable system $^{124}$Sn+$^{40}$Ca. The dynamic microscopic calculations are carried out on a three-dimensional lattice using both the Time-Dependent Hartree-Fock method and the Density Constrained TDHF method. We demonstrate that the peak of the GDR excitation spectrum occurs at a substantially lower energy than expected for an equilibrated system, thus reflecting the very large prolate elongation of the dinuclear complex during the early stages of fusion. Our theoretical fusion cross-sections for both systems agree reasonably well with recent data measured at HRIBF.

preprint2012arXiv

Electric dipole polarizability and the neutron skin

The recent high-resolution measurement of the electric dipole (E1) polarizability (alphad) in 208Pb [Phys. Rev. Lett. 107, 062502 (2011)] provides a unique constraint on the neutron-skin thickness of this nucleus. The neutron-skin thickness (rskin) of 208Pb is a quantity of critical importance for our understanding of a variety of nuclear and astrophysical phenomena. To assess the model dependence of the correlation between alphad and rskin, we carry out systematic calculations for 208Pb, 132Sn, and 48Ca based on the nuclear density functional theory (DFT) using both non-relativistic and relativistic energy density functionals (EDFs). Our analysis indicates that whereas individual models exhibit a linear dependence between alphad and rskin, this correlation is not universal when one combines predictions from a host of different models. By averaging over these model predictions, we provide estimates with associated systematic errors for rskin and alphad for the nuclei under consideration. We conclude that precise measurements of rskin in both 48Ca and 208Pb---combined with the recent measurement of alphad---should significantly constrain the isovector sector of the nuclear energy density functional.

preprint2012arXiv

Energy density functional for nuclei and neutron stars

We aim to develop a nuclear energy density functional that can be simultaneously applied to finite nuclei and neutron stars. We use the self-consistent nuclear density functional theory (DFT) with Skyrme energy density functionals and covariance analysis to assess correlations between observables for finite nuclei and neutron stars. In a first step two energy functionals -- a high density energy functional giving reasonable neutron properties, and a low density functional fitted to nuclear properties -- are matched. In a second step, we optimize a new functional using exactly the same protocol as in earlier studies pertaining to nuclei but now including neutron star data. This allows direct comparisons of performance of the new functional relative to the standard one. The new functional TOV-min yields results for nuclear bulk properties (energy, r.m.s. radius, diffraction radius, surface thickness) that are of the same quality as those obtained with the established Skyrme functionals, including SV-min. When comparing SV-min and TOV-min, isoscalar nuclear matter indicators vary slightly while isovector properties are changed considerably. We discuss neutron skins, dipole polarizability, separation energies of the heaviest elements, and proton and neutron drip lines. We confirm a correlation between the neutron skin of $^{208}$Pb and the neutron star radius. We demonstrate that standard energy density functionals optimized to nuclear data do not carry information on the expected maximum neutron star mass, and that predictions can only be made within an extremely broad uncertainty band. For atomic nuclei, the new functional TOV-min performs at least as well as the standard nuclear functionals, but it also reproduces expected neutron star data within assumed error bands.

preprint2012arXiv

Existence of exotic torus configuration in high-spin excited states of $^{40}$Ca

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 $\hbar$ and an excitation energy of about 170 MeV in all calculations using various Skyrme interactions. The total angular momentum J=60 $\hbar$ consists of aligned 12 nucleons with the orbital angular momenta $Λ=+4$, +5, and +6 for spin up-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.

preprint2012arXiv

Have Superheavy Elements been Produced in Nature?

We discuss the possibility whether superheavy elements can be produced in Nature by the astrophysical rapid neutron capture process. To this end we have performed fully dynamical network r-process calculations assuming an environment with neutron-to-seed ratio large enough to produce superheavy nuclei. Our calculations include two sets of nuclear masses and fission barriers and include all possible fission channels and the associated fission yield distributions. Our calculations produce superheavy nuclei with A ~ 300 that however decay on timescales of days.

preprint2012arXiv

Information content of the low-energy electric dipole strength: correlation analysis

We study the information content carried by the electric dipole strength with respect to isovector and isoscalar indicators characterizing bulk nuclear matter and finite nuclei. To separate isoscalar and isovector modes, and low-energy strength and giant resonances, we analyze the E1 strength as a function of excitation energy $E$ and momentum transfer $q$. We use the self-consistent nuclear density functional theory with Skyrme energy density functionals, augmented by the random phase, to compute the E1 strength, and covariance analysis to assess correlations between observables. Calculations are performed for spherical, doubly-magic nuclei $^{208}$Pb and $^{132}$Sn. We demonstrate that E1 transition densities in the low-energy region below the giant dipole resonance exhibit appreciable state dependence and multi-nodal structures, which are fingerprints of weak collectivity. The correlation between the accumulated low-energy strength and symmetry energy is weak, and dramatically depends on the energy cutoff assumed. On the other hand, a strong correlation is predicted between isovector indicators and the accumulated isovector strength at $E$ around 20 MeV and momentum transfer $q\sim 0.65$ fm$^{-1}$. Momentum- and coordinate-space pattern of the low-energy dipole modes indicate a strong fragmentation into individual particle-hole excitations. The global measure of low-energy dipole strength poorly correlates with the nuclear symmetry energy and other isovector characteristics. Consequently, our results do not support the suggestion that there exists a collective "pygmy dipole resonance," which is a strong indicator of nuclear isovector properties.

preprint2012arXiv

Nuclear energy density optimization: Large deformations

A new Skyrme-like energy density suitable for studies of strongly elongated nuclei has been determined in the framework of the Hartree-Fock-Bogoliubov theory using the recently developed model-based, derivative-free optimization algorithm POUNDerS. A sensitivity analysis at the optimal solution has revealed the importance of states at large deformations in driving the parameterization of the functional. The good agreement with experimental data on masses and separation energies, achieved with the previous parameterization UNEDF0, is largely preserved. In addition, the new energy density UNEDF1 gives a much improved description of the fission barriers in 240Pu and neighboring nuclei.

preprint2012arXiv

Optimization of relativistic mean field model for finite nuclei to neutron star matter

We have optimized the parameters of extended relativistic mean-field model using a selected set of global observables which includes binding energies and charge radii for nuclei along several isotopic and isotonic chains and the iso-scalar giant monopole resonance energies for the $^{90}$Zr and $^{208}$Pb nuclei. The model parameters are further constrained by the available informations on the energy per neutron for the dilute neutron matter and bounds on the equations of state of the symmetric and asymmetric nuclear matter at supra-nuclear densities. Two new parameter sets BSP and IUFSU* are obtained, later one being the variant of recently proposed IUFSU parameter set. The BSP parametrization uses the contributions from the quartic order cross-coupling between $ω$ and $σ$ mesons to model the high density behaviour of the equation of state instead of the $ω$ meson self-coupling as in the case of IUFSU* or IUFSU. Our parameter sets yield appreciable improvements in the binding energy systematics and the equation of state for the dilute neutron matter. The importance of the quartic order $ω-σ$ cross coupling term of the extended RMF model, as often ignored, is realized.

preprint2012arXiv

Phase space structure and dynamics within the time-dependent Hartree-Fock approach

We study the equilibration and relaxation processes within the time-dependent Hartree-Fock approach using the Wigner distribution function. On the technical side we present a geometrically unrestricted framework which allows us to calculate the full six-dimensional Wigner distribution function. With the removal of geometrical constraints, we are now able to extend our previous phase-space analysis of heavy-ion collisions in the reaction plane to unrestricted mean-field simulations of nuclear matter on a three-dimensional Cartesian lattice. From the physical point of view we provide a quantitative analysis on the stopping power in TDHF. This is linked to the effect of transparency. For the medium-heavy $^{40}$Ca+$^{40}$Ca system we examine the impact of different parametrizations of the Skyrme force, energy-dependence, and the significance of extra time-odd terms in the Skyrme functional. For the first time, transparency in TDHF is observed for a heavy system, $^{24}$Mg+$^{208}$Pb.

preprint2012arXiv

Self-consistent calculations of the electric giant dipole resonances in light and heavy mass 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 $^{16}$O, $^{40}$Ca, and $^{208}$Pb are calculated using two new Skyrme interactions.

preprint2012arXiv

Sensitivity of the electric dipole polarizability to the neutron skin thickness in ${}^{208}$Pb

The static dipole polarizability, $α_{\rm D}$, in ${}^{208}$Pb has been recently measured with high-resolution via proton inelastic scattering at the Research Center for Nuclear Physics (RCNP). This observable is thought to be intimately connected with the neutron skin thickness, $r_{\rm skin}$, of the same nucleus and, more fundamentally, it is believed to be associated with the density dependence of the nuclear symmetry energy. The impact of $r_{\rm skin}$ on $α_{\rm D}$ in ${}^{208}$Pb is investigated and discussed on the basis of a large and representative set of relativistic and non-relativistic nuclear energy density functionals (EDF).

preprint2012arXiv

Single-particle dissipation in TDHF studied from a phase-space perspective

We study dissipation and relaxation processes within the time-dependent Hartree-Fock approach using the Wigner distribution function. On the technical side we present a geometrically unrestricted framework which allows us to calculate the full six-dimensional Wigner distribution function. With the removal of geometrical constraints, we are now able to extend our previous phase-space analysis of heavy-ion collisions in the reaction plane to unrestricted mean-field simulations of nuclear matter on a three-dimensional Cartesian lattice. From the physical point of view we provide a quantitative analysis on the stopping power in TDHF. This is linked to the effect of transparency. For the medium-heavy $^{40}$Ca+$^{40}$Ca system we examine the impact of different parametrizations of the Skyrme force, energy-dependence, and the significance of extra time-odd terms in the Skyrme functional.

preprint2012arXiv

Toroidal, compression, and vortical dipole strengths in 124Sn

The toroidal, compression and vortical dipole strength functions in semi-magic $^{124}$Sn (and partly in doubly-magic $^{100,132}$Sn) are analyzed within the random-phase-approximation method with the SkT6, SkI3, SLy6, SV-bas, and SkM* Skyrme forces. The isoscalar (T=0), isovector (T=1), and electromagnetic ('elm') channels are considered. Both convection $j_c$ and magnetization $j_m$ nuclear currents are taken into account. The calculations basically confirm the previous results obtained for $^{208}$Pb with the force SLy6. In particular, it is shown that the vortical and toroidal strengths are dominated by $j_c$ in T=0 channel and by $j_m$ in T=1 and 'elm' channels. The compression strength is always determined by $j_c$. It is also shown that the 'elm' strength (relevant for (e,e') reaction) is very similar to T=1 one. The toroidal mode resides in the region of the pygmy resonance. So, perhaps, this region embraces both irrotational (pygmy) and vortical (toroidal) flows.

preprint2011arXiv

Equilibration in the time-dependent Hartree-Fock approach probed with the Wigner distribution function

Calculating the Wigner distribution function in the reaction plane, we are able to probe the phase-space behavior in time-dependent Hartree-Fock during a heavy-ion collision. We compare the Wigner distribution function with the smoothed Husimi distribution function. Observables are defined to give a quantitative measure for local and global equilibration. We present different reaction scenarios by analyzing central and non-central $^{16}O+$$^{16}O$ and $^{96}Zr+$$^{132}Zn$ collisions. It is shown that the initial phase-space volumes of the fragments barely merge. The mean values of the observables are conserved in fusion reactions and indicate a "memory effect" in time-dependent Hartree-Fock. We observe strong dissipation but no evidence for complete equilibration.

preprint2011arXiv

Fission properties for r-process nuclei

We present a systematics of fission barriers and fission lifetimes for the whole landscape of super-heavy elements (SHE), i.e. nuclei with Z>100. The fission lifetimes are also compared with the alpha-decay half-lives. The survey is based on a self-consistent description in terms of the Skyrme-Hartree-Fock (SHF) approach. Results for various different SHF parameterizations are compared to explore the robustness of the predictions. The fission path is computed by quadrupole constrained SHF. The computation of fission lifetimes takes care of the crucial ingredients of the large-amplitude collective dynamics along the fission path, as self-consistent collective mass and proper quantum corrections. We discuss the different topologies of fission landscapes which occur in the realm of SHE (symmetric versus asymmetric fission, regions of triaxial fission, bi-modal fission, and the impact of asymmetric ground states). The explored region is extended deep into the regime of very neutron-rich isotopes as they are expected to be produced in the astrophysical r process.

preprint2011arXiv

General Treatment of Vortical, Toroidal, and Compression Modes

The multipole vortical, toroidal, and compression modes are analyzed. Following the vorticity concept of Ravenhall and Wambach, the vortical operator is derived and related in a simple way to the toroidal and compression operators. The strength functions and velocity fields of the modes are analyzed in $^{208}$Pb within the random-phase-approximation using the Skyrme force SLy6. Both convection and magnetization nuclear currents are taken into account. It is shown that the isoscalar (isovector) vortical and toroidal modes are dominated by the convection (magnetization) nuclear current while the compression mode is fully convective. The relation between the above concept of the vorticity to the hydrodynamical vorticity is briefly discussed.

preprint2011arXiv

Localization in light nuclei

We investigate the presence of spatial localization in nuclei using a method that maps the nucleon same-spin pair probability and is based on the density-matrix. The method is used to study spatial localization of light nuclei within the Hartree-Fock approximation. We show that the method provides an alternative tool for studying spatial localization in comparison to the localization observed from maxima in the nuclear mass density.

preprint2010arXiv

Angular distributions of electrons emitted from free and deposited Na$_8$ clusters

We explore from a theoretical perspective angular distributions of electrons emitted from a Na$_8$ cluster after excitation by a short laser pulse. The tool of the study is time-dependent density-functional theory (TDDFT) at the level of the local-density approximation (LDA) augmented by a self-interaction correction (SIC) to put emission properties in order. We consider free Na$_8$ and Na$_8$ deposited on the surfaces MgO(001) or Ar(001). For the case of free Na$_8$, we distinguish between a hypothetical situation of known cluster orientation and a more realistic ensemble of orientations. We also consider the angular distributions for emission from separate single-electron levels.

preprint2010arXiv

Deposition of Na Clusters on MgO(001)

We investigate the dynamics of deposition of small Na clusters on MgO(001) surface. A hierarchical modeling is used combining Quantum Mechanical with Molecular Mechanical (QM/MM) description. Full time-dependent density-functional theory is used for the cluster electrons while the substrate atoms are treated at a classical level. We consider Na$_6$ and Na$_8$ at various impact energies. We analyze the dependence on cluster geometry, trends with impact energy, and energy balance. We compare the results with deposit on the much softer Ar(001) surface.

preprint2010arXiv

Energy systematics of heavy nuclei -- mean field models in comparison

We compare the systematics of binding energies computed within the standard and extended versions of the relativistic mean-field (RMF) model and the Skyrme Hartree-Fock (SHF) model. The general trends for the binding energies for super-heavy nuclei are significantly different for these models. The SHF models tend to underbind the superheavy nuclei, while, RMF models show just the opposite trend. The extended RMF model seems to provide remarkable improvements over the results obtained for the standard RMF model.

preprint2010arXiv

Entrance Channel Dynamics of Hot and Cold Fusion Reactions Leading to Superheavy Elements

We investigate the entrance channel dynamics for the reactions $\mathrm{^{70}Zn}+\mathrm{^{208}Pb}$ and $\mathrm{^{48}Ca}+\mathrm{^{238}U}$ using the fully microscopic time-dependent Hartree-Fock (TDHF) theory coupled with a density constraint. We calculate excitation energies and capture cross-sections relevant for the study of superheavy formations. We discuss the deformation dependence of the ion-ion potential for the $\mathrm{^{48}Ca}+\mathrm{^{238}U}$ system and perform an alignment angle averaging for the calculation of the capture cross-section. The results show that this parameter-free approach can generate results in good agreement with experiment and other theories.

preprint2010arXiv

Exploration of a modified density dependence in the Skyrme functional

A variant of the basic Skyrme-Hartree-Fock (SHF) functional is considered dealing with a new form of density dependence. It employs only integer powers and thus will allow a more sound basis for projection schemes (particle number, angular momentum). We optimize the new functional with exactly the same adjustment strategy as used in an earlier study with a standard Skyrme functional. This allows direct comparisons of the performance of the new functional relative to the standard one. We discuss various observables: bulk properties of finite nuclei, nuclear matter, giant resonances, super-heavy elements, and energy systematics. The new functional performs at least as well as the standard one, but offers a wider range of applicability (e.g. for projection) and more flexibility in the regime of high densities.

preprint2010arXiv

Exploration of dynamical regimes of irradiated small protonated water clusters

We explore from a theoretical perspective the dynamical response of small water clusters, (H$_2$O)$_n$H$_3$O$^+$ with $n=1,2,3$, to a short laser pulse for various frequencies, from infrared (IR) to ultra-violet (UV) and intensities (from $6\times10^{13}$ W/cm$^2$ to $5\times10^{14}$ W/cm$^2$). To that end, we use time-dependent local-density approximation for the electrons, coupled to molecular dynamics for the atomic cores (TDLDA-MD). The local-density approximation is augmented by a self-interaction correction (SIC) to allow for a correct description of electron emission. For IR frequencies, we see a direct coupling of the laser field to the very light H$^+$ ions in the clusters. Resonant coupling (in the UV) and/or higher intensities lead to fast ionization with subsequent Coulomb explosion. The stability against Coulomb pressure increases with system size. Excitation to lower ionization stages induced strong ionic vibrations. These maintain rather harmonic pattern in spite of the sizeable amplitudes (often 10% of the bond length).

preprint2010arXiv

Fission of super-heavy nuclei explored with Skyrme forces

We present a large scale survey of life-times for spontaneous fission in the regime of super-heavy elements (SHE), i.e. nuclei with Z=104-122. This is done on the basis of the Skyrme-Hartree-Fock model. The axially symmetric fission path is computed using a quadrupole constraint. Self-consistent cranking is used for the collective masses and associated quantum corrections. The actual tunneling probability is estimated by the WKB approximation. Three typical Skyrme forces are used to explore the sensitivity of the results. Benchmarks in the regime Z=104-108 show an acceptable agreement. The general systematics reflects nicely the islands of shell stabilization and the crossover from $α$-decay to fission for the decay chains from the region of Z/N=118/176.

preprint2010arXiv

Microscopic Description of Nuclear Fission Dynamics

We discuss possible avenues to study fission dynamics starting from a time-dependent mean-field approach. Previous attempts to study fission dynamics using the time-dependent Hartree-Fock (TDHF) theory are analyzed. We argue that different initial conditions may be needed to describe fission dynamics depending on the specifics of the fission phenomenon and propose various approaches towards this goal. In particular, we provide preliminary calculations for studying fission following a heavy-ion reaction using TDHF with a density contraint. Regarding prompt muon-induced fission, we also suggest a new approach for combining the time-evolution of the muonic wave function with a microscopic treatment of fission dynamics via TDHF.

preprint2010arXiv

Microscopic study of the $^{132,124}$Sn+$^{96}$Zr reactions: dynamic excitation energy, energy-dependent heavy-ion potential, and capture cross section

We study reactions between neutron-rich $^{132}$Sn nucleus and $^{96}$Zr within a dynamic microscopic theory at energies in the vicinity of the ion-ion potential barrier peak, and we compare the properties to those of the stable system $^{124}$Sn+$^{96}$Zr. The calculations are carried out on a three-dimensional lattice using the density-constrained Time-Dependent Hartree-Fock method. In particular, we calculate the dynamic excitation energy $E^*(t)$ and the quadrupole moment of the dinuclear system, $Q_{20}(t)$, during the initial stages of the heavy-ion collision. Capture cross sections for the two reactions are analyzed in terms of dynamic effects and a comparison with recently measured data is given.

preprint2010arXiv

Misfits in Skyrme-Hartree-Fock

We address very briefly five critical points in the context of the Skyrme-Hartree-Fock (SHF) scheme: 1) the impossibility to consider it as an interaction, 2) a possible inconsistency of correlation corrections as, e.g., the center-of-mass correction, 3) problems to describe the giant dipole resonance (GDR) simultaneously in light and heavy nuclei, 4) deficiencies in the extrapolation of binding energies to super-heavy elements (SHE), and 5) a yet inappropriate trend in fission life-times when going to the heaviest SHE. While the first two points have more a formal bias, the other three points have practical implications and wait for solution.

preprint2010arXiv

Properties of odd nuclei and the impact of time-odd mean fields: A systematic Skyrme-Hartree-Fock analysis

We present a systematic analysis of the description of odd nuclei by the Skyrme-Hartree-Fock approach augmented with pairing in BCS approximation and blocking of the odd nucleon. Current and spin densities in the Skyrme functional produce time-odd mean fields (TOMF) for odd nuclei. Their effect on basic properties (binding energies, odd-even staggering, separation energies and spectra) is investigated for the three Skyrme parameterizations SkI3, SLy6, and SV-bas. About 1300 spherical and axially-deformed odd nuclei with 16 < Z < 92 are considered. The calculations demonstrate that the TOMF effect is generally small, although not fully negligible. The influence of the Skyrme parameterization and the consistency of the calculations are much more important. With a proper choice of the parameterization, a good description of binding energies and their differences is obtained, comparable to that for even nuclei. The description of low-energy excitation spectra of odd nuclei is of varying quality depending on the nucleus.

preprint2010arXiv

Skyrme-Hartree-Fock Description of the Dipole Strength in Neutron-Rich Tin Isotopes

Low-energy E1 strength in neutron-rich $^{132-164}$Sn isotopes is analyzed in the framework of the Skyrme random phase approximation (RPA) with different Skyrme forces. A double folding procedure is applied to take into account the energy-dependent width effects beyond RPA. All the considered Skyrme forces indicate a soft prolate deformation in the open shell isotopes $^{142-164}$Sn. The integrated E1 strength in the energy region of the pygmy resonance grows with the neutron number. The influence of deformation on the integrated strength near the particle emission thresholds (which is of a keen astrophysical interest) is strictly suppressed by the mutual compensation effect for the branches of the giant dipole resonance. The results obtained are in a good agreement with the previous findings of the relativistic mean field model.

preprint2010arXiv

Spin-flip M1 giant resonance as a challenge for Skyrme forces

Despite a great success of the Skyrme mean-field approach in exploration of nuclear dynamics, it seems to fail in description of the spin-flip M1 giant resonance. The results for different Skyrme parameterizations are contradictory and poorly agree with experiment. In particular, there is no parameterization which simultaneously describes the one-peak gross structure of M1 strength in doubly magic nuclei and two-peak structure in heavy deformed nuclei. The reason of this mismatch could lie in an unsatisfactory treatment of spin correlations and spin-orbit interaction. We discuss the present status of the problem and possible ways of its solution. In particular, we inspect i) the interplay of the collective shift and spin-orbit splitting, ii) the isovector M1 response versus isospin-mixed responses, and iii) the role of tensor and isovector spin-orbit interaction.

preprint2010arXiv

The Generalized SIC-OEP formalism and the Generalized SIC-Slater approximation (stationary and time-dependent cases)

We present a generalized formulation of the Optimized Effective Potential (OEP) approach to the Self Interaction Correction (SIC) problem in Time Dependent (TD) Density Functional Theory (DFT). The formulation relies on the introduction of a double set of single electron orbitals. It allows the derivation of a generalized Slater approximation to the full OEP formulation, which extends the domain of validity of the standard Slater approximation. We discuss both formal aspects and practical applications of the new formalism and give illustrations in cluster and molecules. The new formalism provides a valuable ansatz to more elaborate (and computationally very demanding) full TD OEP and full TD SIC calculations especially in the linear domain.

preprint2010arXiv

The information content of a new observable: the case of the nuclear neutron skin

We address two questions pertaining to the uniqueness and usefulness of a new observable: (i) Considering the current theoretical knowledge, what novel information does new measurement bring in? (ii) How can new data reduce uncertainties of current theoretical models? We illustrate these points by studying the radius of the neutron distribution of a heavy nucleus, a quantity related to the equation of state for neutron matter that determines properties of nuclei and neutron stars. By systematically varying parameters of two theoretical models and studying the resulting confidence ellipsoid, we quantify the relationships between the neutron skin and various properties of finite nuclei and infinite nuclear matter. Using the covariance analysis, we identify observables and pseudo-observables that correlate, and do not correlate, with the neutron skin. By adding the information on the neutron radius to the pool of observables determining the energy functional, we show how precise experimental determination of the neutron radius in $^{208}$Pb would reduce theoretical uncertainties on the neutron matter equation of state.

preprint2009arXiv

Dynamics of clusters and molecules in contact with an environment

We present recent theoretical investigations on the dynamics of metal clusters in contact with an environment, deposited of embedded. This concerns soft deposition as well as irradiation of the deposited/embedded clusters by intense laser pulses. We discuss examples of applications for two typical test cases, Na clusters deposited on MgO(001) surface and Na clusters in/on Ar substrate. Both environments are insulators with sizeable polarizability. They differ in their geometrical and mechanical properties.

preprint2009arXiv

Laser-driven nonlinear cluster dynamics

Laser excitation of nanometer-sized atomic and molecular clusters offers various opportunities to explore and control ultrafast many-particle dynamics. Whereas weak laser fields allow the analysis of photoionization, excited-state relaxation, and structural modifications on these finite quantum systems, large-amplitude collective electron motion and Coulomb explosion can be induced with intense laser pulses. This review provides an overview of key phenomena arising from laser-cluster interactions with focus on nonlinear optical excitations and discusses the underlying processes according to the current understanding. A brief general survey covers basic cluster properties and excitation mechanisms relevant for laser-driven cluster dynamics. Then, after an excursion in theoretical and experimental methods, results for single- and multiphoton excitations are reviewed with emphasis on signatures from time- and angular resolved photoemission. A key issue of this review is the broad spectrum of phenomena arising from clusters exposed to strong fields, where the interaction with the laser pulse creates short-lived and dense nanoplasmas. The implications for technical developments include the controlled generation of ion, electron, and radiation pulses, as will be addressed along with corresponding examples. Finally, future prospects of laser-cluster research as well as experimental and theoretical challenges are discussed.

preprint2009arXiv

Polarizibilities as a test of localized approximations to the self-interaction correction

We present applications of the recently introduced ``Generalized SIC-Slater'' scheme which provides a simple Self-Interaction Correction approximation in the framework of the Optimized Effective Potential. We focus on the computation of static polarizabilities which are known to constitute stringent tests for Density Functional Theory. We apply the new method to model H chains, but also to more realistic systems such as C4 (organic) chains, and less symmetrical systems such as a Na5 (metallic) cluster. Comparison is made with other SIC schemes, especially with the standard SIC-Slater one.

preprint2009arXiv

Skyrme-Random-Phase-Approximation description of spin-flip and orbital giant resonances

The self-consistent separable random-phase approximation (SRPA) model with Skyrme forces is extended to the case of magnetic excitations and applied to the description of spin-flip and orbital M1 giant resonances in the isotopic chain $^{142-152}$Nd. The Skyrme forces SkT6, SkM*, SLy6 and SkI3 are used. The calculations show the onset of the scissors mode with increasing deformation. A specific three-peak structure of the spin-flip response is found and explained by particular neutron and proton spin-flip transitions. Although the employed forces provide an acceptable qualitative description, the Skyrme functional still needs further improvement to reproduce quantitatively the experimental data.

preprint2006arXiv

Exploration of electronic quadrupole states in atomic clusters by two-photon processes

We analyze particular two-photon processes as possible means to explore electronic quadrupole states in free small deformed atomic clusters. The analysis is done in the time-dependent local density approximation (TDLDA). It is shown that the direct two-photon population (DTP) and off-resonant stimulated Raman (ORSR) scattering can be effectively used for excitation of the quadrupole states in high-frequency (quadrupole plasmon) and low-frequency (infrared) regions, respectively. In ORSR, isolated dipole particle-hole states as well as the tail of the dipole plasmon can serve as an intermediate state. A simultaneous study of low- and high-frequency quadrupoles, combining DTP and ORSR, is most effective. Femtosecond pulses with intensities $I = 2\cdot 10^{10} - 2\cdot 10^{11} W/cm^2$ and pulse durations $T = 200 - 500$ fs are found to be optimal. Since the low-lying quadrupole states are dominated by one single electron-hole pair, their energies, being combined with the photoelectron data for hole states, allow to get the electron spectrum above the Fermi level and thus greatly extend our knowledge on the single particle spectra of clusters. Besides, the developed schemes allow to estimate the lifetime of the quadrupole states.

preprint2005arXiv

Electronic excitations in atomic clusters: beyond dipole plasmon

Multipole electron modes beyond the Mie plasmon in atomic clusters are investigated within the time-dependent local density approximation theory (TDLDA). We consider the origin of the modes, their connection with basic cluster properties and possible routes of experimental observation. Particular attention is paid to infrared magnetic orbital modes, scissors and twist, and electric quadrupole mode. The scissors and twist modes determine orbital magnetism of clusters while the electric quadrupole mode provides direct access to the single electron spectra of the cluster. We examine two-photon processes (Raman scattering, stimulated emission pumping and stimulated adiabatic Raman passage) as the most promising tools for experimental investigation of the modes.

preprint2005arXiv

Practicable factorized TDLDA for arbitrary density- and current-dependent functionals

We propose a practicable method for describing linear dynamics of different finite Fermi systems. The method is based on a general self-consistent procedure for factorization of the two-body residual interaction. It is relevant for diverse density- and current-dependent functionals and, in fact, represents the self-consistent separable random-phase approximation (RPA), hence the name SRPA. SRPA allows to avoid diagonalization of high-rank RPA matrices and thus dwarfs the calculation expense. Besides, SRPA expressions have a transparent analytical form and so the method is very convenient for the analysis and treatment of the obtained results. SRPA demonstrates high numerical accuracy. It is very general and can be applied to diverse systems. Two very different cases, the Kohn-Sham functional for atomic clusters and Skyrme functional for atomic nuclei, are considered in detail as particular examples. SRPA treats both time-even and time-odd dynamical variables and, in this connection, we discuss the origin and properties of time-odd currents and densities in initial functionals. Finally, SRPA is compared with other self-consistent approaches for the excited states, including the coupled-cluster method.

preprint2005arXiv

Self-Consistent Separable Rpa Approach for Skyrme Forces: Axial Nuclei

The self-consistent separable RPA (random phase approximation) method is formulated for Skyrme forces with pairing. The method is based on a general self-consistent procedure for factorization of the two-body interaction. It is relevant for various density- and current-dependent functionals. The contributions of the time-even and time-odd Skyrme terms as well as of the Coulomb and pairing terms to the residual interaction are taken self-consistently into account. Most of the expression have a transparent analytical form, which makes the method convenient for the treatment and analysis. The separable character of the residual interaction allows to avoid diagonalization of high-rank RPA matrices and thus to minimize the calculation effort. The previous studies have demonstrated high numerical accuracy and efficiency of the method for spherical nuclei. In this contribution, the method is specified for axial nuclei. We provide systematic and detailed presentation of formalism and discuss different aspects of the model.

preprint2005arXiv

Two-Photon Excitation of Low-Lying Electronic Quadrupole States in Atomic Clusters

A simple scheme of population and detection of low-lying electronic quadrupole modes in free small deformed metal clusters is proposed. The scheme is analyzed in terms of the TDLDA (time-dependent local density approximation) calculations. As test case, the deformed cluster $Na^+_{11}$ is considered. Long-living quadrupole oscillations are generated via resonant two-photon (two-dipole) excitation and then detected through the appearance of satellites in the photoelectron spectra generated by a probe pulse. Femtosecond pump and probe pulses with intensities $I = 2\cdot 10^{10} - 2\cdot 10^{11} W/cm^2$ and pulse duration $T = 200 - 500$ fs are found to be optimal. The modes of interest are dominated by a single electron-hole pair and so their energies, being combined with the photoelectron data for hole states, allow to gather new information about mean-field spectra of valence electrons in the HOMO-LUMO region. Besides, the scheme allows to estimate the lifetime of electron-hole pairs and hence the relaxation time of electronic energy into ionic heat.

preprint2004arXiv

Systematics of Fission Barriers in Superheavy Elements

We investigate the systematics of fission barriers in superheavy elements in the range Z = 108-120 and N = 166-182. Results from two self-consistent models for nuclear structure, the relativistic mean-field (RMF) model as well as the non-relativistic Skyrme-Hartree-Fock approach are compared and discussed. We restrict ourselves to axially symmetric shapes, which provides an upper bound on static fission barriers. We benchmark the predictive power of the models examining the barriers and fission isomers of selected heavy actinide nuclei for which data are available. For both actinides and superheavy nuclei, the RMF model systematically predicts lower barriers than most Skyrme interactions. In particular the fission isomers are predicted too low by the RMF, which casts some doubt on recent predictions about superdeformed ground states of some superheavy nuclei. For the superheavy nuclei under investigation, fission barriers drop to small values around Z = 110, N = 180 and increase again for heavier systems. For most of the forces, there is no fission isomer for superheavy nuclei, as superdeformed states are in most cases found to be unstable with respect to octupole distortions.

preprint2003arXiv

Generator Coordinate Truncations

We investigate the accuracy of several schemes to calculate ground-state correlation energies using the generator coordinate technique. Our test-bed for the study is the $sd$ interacting boson model, equivalent to a 6-level Lipkin-type model. We find that the simplified projection of a triaxial generator coordinate state using the $S_3$ subgroup of the rotation group is not very accurate in the parameter space of the Hamiltonian of interest. On the other hand, a full rotational projection of an axial generator coordinate state gives remarkable accuracy. We also discuss the validity of the simplified treatment using the extended Gaussian overlap approximation (top-GOA), and show that it works reasonably well when the number of boson is four or larger.

preprint2003arXiv

The nonrelativistic limit of the relativistic point coupling model

We relate the relativistic finite range mean-field model (RMF-FR) to the point-coupling variant and compare the nonlinear density dependence. From this, the effective Hamiltonian of the nonlinear point-coupling model in the nonrelativistic limit is derived. Different from the nonrelativistic models, the nonlinearity in the relativistic models automatically yields contributions in the form of a weak density dependence not only in the central potential but also in the spin-orbit potential. The central potential affects the bulk and surface properties while the spin-orbit potential is crucial for the shell structure of finite nuclei. A modification in the Skyrme-Hartree-Fock model with a density-dependent spin-orbit potential inspired by the point-coupling model is suggested.

preprint2001arXiv

Pairing Interaction and Self-Consistent Densities in Neutron-Rich Nuclei

Particle and pairing densities in spherical even-even neutron-rich nuclei are studied within the Skyrme-Hartree-Fock-Bogoliubov approach with the density-dependent pairing interaction. The influence of the density dependence of the pairing interaction on asymptotic properties of nucleonic distributions are analyzed. It is demonstrated that the size of the neutron halo dramatically depends on the behavior of the pairing interaction at low density.

preprint2000arXiv

Pairing gaps from nuclear mean-field models

We discuss the pairing gap, a measure for nuclear pairing correlations, in chains of spherical, semi-magic nuclei in the framework of self-consistent nuclear mean-field models. The equations for the conventional BCS model and the approximate projection-before-variation Lipkin-Nogami method are formulated in terms of local density functionals for the effective interaction. We calculate the Lipkin-Nogami corrections of both the mean-field energy and the pairing energy. Various definitions of the pairing gap are discussed as three-point, four-point and five-point mass-difference formulae, averaged matrix elements of the pairing potential, and single-quasiparticle energies. Experimental values for the pairing gap are compared with calculations employing both a delta pairing force and a density-dependent delta interaction in the BCS and Lipkin-Nogami model. Odd-mass nuclei are calculated in the spherical blocking approximation which neglects part of the the core polarization in the odd nucleus. We find that the five-point mass difference formula gives a very robust description of the odd-even staggering, other approximations for the gap may differ from that up to 30% for certain nuclei.

preprint1999arXiv

Pairing Gap and Polarisation Effects

The phenomenological adjustment of the nuclear pairing strength is usually performed with respect to the odd-even staggering of the binding energies. We find that the results strongly depend on the way in which the ground states of the odd nuclei are computed. A thorough calculation including all time-even and time-odd polarisation effects induced by the odd nucleon produces about 30% reduced odd-even staggering as compared to the standard spherical calculations in the relativistic mean-field model. The pairing strength must be enhanced by about 20% to compensate for that effect. The enhanced strength has dramatic consequences for the predicted deformation properties of the underlying mean-field models, possibly implying that new adjustments of their parameters become necessary as well.