Source author record

M. Bender

M. Bender 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

47works
7topics
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

47 published item(s)

preprint2026arXiv

Skyrme-Hartree-Fock-Bogoliubov mass models on a 3D mesh: V. The N2LO extension of the Skyrme EDF

We present BSkG5, the latest entry in the Brussels-Skyrme-on-a-Grid (BSkG) series and the first large-scale nuclear structure model based on next-to-next-to-leading order (N2LO) Skyrme energy density functional (EDF). By extending the traditional Skyrme EDF ansatz with central terms containing up to four gradients, we are able to combine an excellent global description of nuclear ground state properties with a stiff equation of state for pure neutron matter that is consistent with all astronomical observations of neutron stars. More precisely, the new model matches the accuracy of earlier BSkG models but with two parameters less: we achieve root-mean-square deviations of 0.649 MeV for 2457 atomic masses, 0.0267 fm for 810 charge radii, and 0.43 MeV for 45 primary fission barriers of actinide nuclei. We demonstrate that the complexities of N2LO EDFs are not insurmountable, even for demanding many-body calculations.

preprint2021arXiv

$τ$ lepton mass measurement at Belle II

The reconstruction of tau-pair production, $e^{+}e^{-} \to τ^{+}τ^{-}$, from the subsequent 3-prong ($τ^{+} \rightarrow π^{+} π^{-} π^{+} \barν_τ$) and 1-prong ($τ^{-} \to \ell^{-} \barν_{\ell} ν_τ$, $τ^{-} \to h^{-} ν_τ$ or $τ^{-} \to π^{-} π^0 ν_τ$) decays, is presented using 8.8 fb$^{-1}$ of $e^{+}e^{-}$ collision data of Belle II at the center-of-mass energy $\sqrt{s} = m_{Υ(4S)}$. The pseudomass technique developed by the ARGUS experiment is used to measure the $τ$-lepton mass $m_τ$ in the 3-prong $τ^{+} \to π^{+} π^{-} π^{+} \barν_τ $ decay, resulting in $m_τ = 1777.28 \pm 0.75~{\rm (stat.)} \pm 0.33 ~{\rm (sys.)}~{\rm{MeV}/\rm{c}^2}$.

preprint2020arXiv

A calibration of the Belle II hadronic tag-side reconstruction algorithm with $B \rightarrow X\ell ν$ decays

Tag-side reconstruction is an important method for reconstructing $B$ meson decays with missing energy. The Belle II tag-side reconstruction algorithm, Full Event Interpretation, relies on a hierarchical reconstruction of $B$ meson decays with multivariate classification employed at each stage of reconstruction. Given the large numbers of classifiers employed and decay chains reconstructed, the performance of the algorithm on data and simulation differs significantly. Here, calibration factors are derived for hadronic tag-side $B$ decays by measuring a signal side decay, $B \rightarrow X\ell ν$, in $34.6$ fb$^{-1}$ of Belle II data. For a very loose selection on the tag-side $B$ multivariate classifier, the calibration factors are $0.65 \pm 0.02$ and $0.83 \pm 0.03$ for tag-side $B^{+}$ and $B^{0}$ mesons, respectively.

preprint2020arXiv

Charmless $B$ decay reconstruction in 2019 Belle II data

We report on the reconstruction of various charmless $B$ decays from electron-positron collisions at the energy corresponding to the $Υ(4S)$ resonance collected with the Belle II detector at the SuperKEKB collider. We use simulation to devise optimized event selections and apply them to the full data set collected in 2019, corresponding to 8.7\,fb$^{-1}$ of integrated luminosity. We fit the difference between half of the collision energy and the $B$ candidate energy (in the $Υ(4S)$ frame) for events restricted to a signal-rich range in beam-energy-constrained mass to search for charmless signals. Signal yields of approximately 80, 15, 20, 30, 90, and 160 decays are reconstructed for the channels $B^0 \to K^+π^-$, $B^0 \to π^+π^-$, $B^+ \to K^0_S(\to π^+π^-)π^+$, $B^+ \to K^+π^0(\to γγ)$, $B^+ \to K^+K^-K^+$, and $B^+ \to K^+π^-π^+$, respectively. Yields and background contaminations are compatible with those expected from simulation and comparable with those obtained by the Belle experiment. The results show a good understanding of the detector performance and offer a reliable basis to assess projections for future reach.

preprint2020arXiv

Exclusive $B^0 \to π^- \ell^+ ν_\ell$ Decays with Hadronic Full Event Interpretation Tagging in 34.6 fb$^{-1}$ of Belle II Data

We present the results of the re-discovery of the decay $B^0 \to π^- \ell^+ ν_\ell$ in 34.6 fb$^{-1}$ of Belle II data using hadronic $B$-tagging via the Full Event Interpretation algorithm. We observe 21 signal events on a background of 155 in a fit to the distribution of the square of the missing mass, $M_{\mathrm{miss}}^2$, with a significance of 5.69$σ$, and determine a total branching fraction of (1.58 $\pm$ 0.43$_{\mathrm{stat}}$ $\pm$ 0.07$_{\mathrm{sys}}$) $\times 10^{-4}$.

preprint2020arXiv

First flavor tagging calibration using 2019 Belle II data

We report on the first calibration of the standard Belle II $B$-flavor tagger using the full data set collected at the $Υ(4{\rm S})$ resonance in 2019 with the Belle II detector at the SuperKEKB collider, corresponding to 8.7 fb$^{-1}$ of integrated luminosity. The calibration is performed by reconstructing various hadronic charmed $B$-meson decays with flavor-specific final states. We use simulation to optimize our event selection criteria and to train the flavor tagging algorithm. We determine the tagging efficiency and the fraction of wrongly identified tag-side $B$~candidates from a measurement of the time-integrated $B^0-\overline{B}^0$ mixing probability. The total effective efficiency is measured to be $\varepsilon_{\rm eff} = \big(33.8 \pm 3.6(\text{stat}) \pm 1.6(\text{sys})\big)\%$, which is in good agreement with the predictions from simulation and comparable with the best one obtained by the Belle experiment. The results show a good understanding of the detector performance and offer a basis for future calibrations.

preprint2020arXiv

First glimpse of the $N=82$ shell closure below $Z=50$ from masses of neutron-rich cadmium isotopes and isomers

We probe the $N=82$ nuclear shell closure by mass measurements of neutron-rich cadmium isotopes with the ISOLTRAP spectrometer at ISOLDE-CERN. The new mass of $^{132}$Cd offers the first value of the $N=82$, two-neutron shell gap below $Z=50$ and confirms the phenomenon of mutually enhanced magicity at $^{132}$Sn. Using the recently implemented phase-imaging ion-cyclotron-resonance method, the ordering of the low-lying isomers in $^{129}$Cd and their energies are determined. The new experimental findings are used to test large-scale shell-model, mean-field and beyond-mean-field calculations, as well as the ab initio valence-space in-medium similarity renormalization group.

preprint2020arXiv

In-beam gamma-ray and electron spectroscopy of $^{249,251}$Md

The odd-$Z$ $^{251}$Md nucleus was studied using combined $γ$-ray and conversion-electron in-beam spectroscopy. Besides the previously observed rotational band based on the $[521]1/2^-$ configuration, another rotational structure has been identified using $γ$-$γ$ coincidences. The use of electron spectroscopy allowed the rotational bands to be observed over a larger rotational frequency range. Using the transition intensities that depend on the gyromagnetic factor, a $[514]7/2^-$ single-particle configuration has been inferred for this band, i.e., the ground-state band. A physical background that dominates the electron spectrum with an intensity of $\simeq$ 60% was well reproduced by simulating a set of unresolved excited bands. Moreover, a detailed analysis of the intensity profile as a function of the angular momentum provided a method for deriving the orbital gyromagnetic factor, namely $g_K = 0.69^{+0.19}_{-0.16}$ for the ground-state band. The odd-$Z$ $^{249}$Md was studied using $γ$-ray in-beam spectroscopy. Evidence for octupole correlations resulting from the mixing of the $Δl = Δj = 3$ $[521]3/2^-$ and $[633]7/2^+$ Nilsson orbitals were found in both $^{249,251}$Md. A surprising similarity of the $^{251}$Md ground-state band transition energies with those of the excited band of $^{255}$Lr has been discussed in terms of identical bands. Skyrme-Hartree-Fock-Bogoliubov calculations were performed to investigate the origin of the similarities between these bands.

preprint2020arXiv

Measurement of Hadronic Mass Moments $\langle M_X^n \rangle $ in $B \rightarrow X_c \ell ν$ Decays at Belle II

We present measurements of the first six hadronic mass moments in semileptonic $B \rightarrow X_c \ell ν$ decays. The hadronic mass moments, together with other observables of inclusive $B$ decays, can be used to determine the CKM matrix element $|{V_{cb}}|$ and mass of the $b$-quark $m_b$ in the context of Heavy Quark Expansions of QCD. The Belle~II data recorded at the $Υ(4S)$ resonance in 2019 and 2020 (March-July), corresponding to an integrated luminosity of $34.6\;\mathrm{fb}^{-1}$, is used for this measurement. The decay $Υ(4S) \rightarrow B \overline{B}$ is reconstructed by applying the hadronic tagging algorithm provided by the Full Event Interpretation to fully reconstruct one $B$ meson. The second $B$ meson is reconstructed inclusively by selecting a high-momentum lepton. The $X_c$ system is identified by the remaining reconstructed tracks and clusters in the electromagnetic calorimeter. We report preliminary results for the hadronic mass moments $\langle M_X^n \rangle $ with $n=1,\dots,6$, measured as a function of a lower cut on the lepton momentum in the signal $B$ rest frame.

preprint2020arXiv

Measurement of the $B^0$ lifetime using fully reconstructed hadronic decays in the 2019 Belle II dataset

This document presents the measurement of $B^0$ meson lifetimes using the 2019 Belle II dataset that corresponds to an integrated luminosity of $8.7 \pm 0.2$ fb$^{-1}$. Each candidate is fully reconstructed with hadronic decay final states on the signal side, while the rest-of-event technique allows to infer the decay vertex position on the other (tag) side. $B^0$ lifetime is extracted from an unbinned maximum likelihood fit to the distribution of the difference between the signal side $B^0$ candidate and the tag side decay times. The measured lifetime is $τ_{B^0} = 1.48 \pm 0.28 \pm 0.06 $ ps, where the first uncertainty is statistical and the second is systematic.

preprint2020arXiv

Measurement of the branching fraction ${\cal B}(\bar B^0\to D^{*+}\ell^-\barν_\ell)$ with early Belle II data

We measure the branching fractions for the decays $\bar{B}^0 \to D^{*+} e^{-} \barν_e$ and $\bar{B}^0 \to D^{*+} μ^{-} \barν_μ$ using $8.70 \pm 0.09~{\rm fb}^{-1}$ of data collected by the Belle II experiment at the SuperKEKB asymmetric-energy $e^+ e^-$ collider. Candidate signal decays are reconstructed with the subsequent decays $D^{*+}\to D^0 π^+$ and $D^0\to K^-π^+$. We obtain the results ${\cal B}(\bar{B}^0 \to D^{*+} e^{-} \barν_e)$ = $(4.55\pm0.14(\mathrm{stat})\pm0.35 (\mathrm{syst})) \%$ and ${\cal B}(\bar{B}^0 \to D^{*+} μ^{-} \barν_μ)$ = $(4.84\pm 0.13(\mathrm{stat})\pm0.37(\mathrm{syst})) \%$, in agreement with the world averages. The measurements serve to validate the full chain of detector operation and calibration, data collection and processing, and production of physics results in the case of semileptonic $B$-meson decays.

preprint2020arXiv

Measurement of the semileptonic $\bar{B}^0 \to D^{*+} \ell^{-} ν_{\ell}$ branching fraction with fully reconstructed $B$ meson decays and 34.6 fb$^{-1}$of Belle II data

We present a first measurement of the $\bar{B^{0}} \rightarrow D^{*+} \ell^{-} ν_{\ell}$ branching fraction using fully reconstructed $B$ meson decays employing the Full Event Interpretation algorithm. Collision events corresponding to an integrated luminosity of \lumi are analyzed, which were recorded by the Belle~II detector operated at the SuperKEKB accelerator complex. We measure$\cal{B}(\bar{B^{0}} \rightarrow D^{*+} \ell^{-} ν_{\ell}) =4.51 \pm 0.41_{stat}\pm0.27_{syst} \pm0.45_{π_s}$, with the first and second error denoting the statistical and systematic uncertainty, respectively, and the third dominant uncertainty is from the slow pion reconstruction efficiency.

preprint2020arXiv

Measurements of branching fractions and CP-violating charge asymmetries in charmless $B$ decays reconstructed in 2019--2020 Belle~II data

We report on first measurements of branching fractions~($\mathcal{B}$) and CP-violating charge asymmetries~($\mathcal{A}$) in charmless $B$ decays at Belle~II. We use a sample of electron-positron collisions collected in 2019 and 2020 at the $Υ(4S)$ resonance and corresponding to $34.6$\,fb$^{-1}$ of integrated luminosity. We use simulation to determine optimized event selections. The $ΔE$ distributions of the resulting samples, restricted in $M_{\rm bc}$, are fit to determine signal yields. Signal yields are corrected for efficiencies determined from simulation and control data samples to obtain branching fractions and CP-violating asymmetries for flavour-specific channels. These are the first measurements in charmless decays reported by Belle~II. Results are compatible with known determinations and show detector performance comparable with the best Belle results offering a reliable basis to assess projections for future reach.

preprint2020arXiv

Rediscovery of $B \to ϕK^{(*)}$ decays and measurement of the longitudinal polarization fraction $f_L$ in $B \to ϕK^{*}$ decays using the Summer 2020 Belle II dataset

We utilize a sample of 34.6 fb$^{-1}$, collected by the Belle II experiment at the SuperKEKB asymmetric energy $e^+e^-$ collider, to search for the $B^+ \to ϕK^+$, $B^+ \to ϕK^{*+}$, $B^0 \to ϕK^0_S$, and $B^0 \to ϕK^{*0}$ decays. Charmless hadronic $B$ decays represent an important part of the Belle II physics program, and are an ideal benchmark to test the detector capabilities in terms of tracking efficiency, charged particle identification, vertexing, and advanced analysis techniques. Each channel is observed with a significance that exceeds 5 standard deviations, and we obtain measurements of their branching ratios that are in good agreement with the world averages. For the $B \to ϕK^*$ modes, we also perform a measurement of the longitudinal polarization fraction $f_L$.

preprint2020arXiv

Studies of the semileptonic $\bar B^0\to D^{*+}\ell^-\barν_\ell$ and $B^-\to D^{0}\ell^-\barν_\ell$ decay processes with 34.6 fb$^{-1}$ of Belle II data

We report measurements of the $\bar{B}^0 \to D^{*+} \ell^{-} \barν_l$ and $B^- \to D^{0} \ell^{-} \barν_l$ processes using 34.6 fb$^{-1}$ of collision events recorded by the Belle II experiment at the SuperKEKB asymmetric-energy $e^+ e^-$ collider. For the $B^-\to D^{0}\ell^-\barν_\ell$ channel, we present first studies that isolate this decay from other semileptonic processes and backgrounds. We report a measurement of the $\bar{B}^0 \to D^{*+} \ell^{-} \barν_l$ branching fraction and obtain ${\cal B}(\bar{B}^0 \to D^{*+} \ell^{-} \barν_l) = \left(4.60 \pm 0.05_{\mathrm{stat}}\pm0.17_{\mathrm{syst}} \pm 0.45_{π_s}\right) \%$, in agreement with the world average. Here, the uncertainties are statistical, systematic, and related to slow pion reconstruction, respectively. The systematic uncertainties are limited by the statistics of auxiliary measurements and will improve in the future. We also report differential branching fractions in five bins of the hadronic recoil parameter $w$ for $\bar{B}^0 \to D^{*+} \ell^{-} \barν_l$, unfolded to account for resolution and efficiency effects.

preprint2016arXiv

Shell-model calculation of isospin-symmetry breaking correction to superallowed Fermi beta-decay

We investigate the radial-overlap part of the isospin-symmetry breaking correction to superallowed $0^+\to 0^+$-decay using the shell-model approach similar to that of Refs. [1, 2]. The 8 sd-shell emitters with masses between $A=22$ and $A=38$ have been re-examined. The Fermi matrix element is evaluated with realistic spherical single-particle wave functions, obtained from spherical Woods-Saxon (WS) or Hartree-Fock (HF) potentials, fine-tuned to reproduce the experimental data on charge radii and separation energies for nuclei of interest. The elaborated adjustment procedure removes any sensitivity of the correction to a specific parametrisation of the WS potential or to various versions of the Skyrme interaction. The present results are generally in good agreement with those reported in Refs. [3, 4]. At the same time, we find that the calculations with HF wave functions result in systematically lower values of the correction.

preprint2016arXiv

Symmetry-unrestricted Skyrme mean-field study of heavy nuclei

In the light of recent experimental developments, increasing attention is devoted to nuclear phenomena related to rotational excitations of exotic intrinsic nuclear configurations that often lack symmetries present in the majority of nuclei. Examples include configurations with a non-vanishing octupole moment. In order to describe this kind of states, we have developed a new computer code to solve the self-consistent mean-field equations, able to use most of today's effective Skyrme interactions and working in coordinate-space. We report on the development of MOCCa, a code based on the same principles as EV8, but offering the user individual control on many symmetry assumptions. In addition, the HF+BCS pairing treatment of EV8 has been generalised to the full machinery of Hartree-Fock-Bogoliubov transformations. We discuss as example the static fission barrier of $^{226}$Ra, prefacing extended studies in the region, using the recent series of Skyrme parameterizations SLy5s1 through SLy5s8.

preprint2015arXiv

Ab initio-driven nuclear energy density functional method. A proposal for safe/correlated/improvable parametrizations of the off-diagonal EDF kernels

This programmatic paper lays down the possibility to reconcile the necessity to resum many-body correlations into the energy kernel with the fact that safe multi-reference energy density functional (EDF) calculations cannot be achieved whenever the Pauli principle is not strictly enforced, as is for example the case when many-body correlations are parametrized under the form of empirical density dependencies. Our proposal is to exploit a newly developed ab initio many-body formalism to guide the construction of safe, explicitly correlated and systematically improvable parametrizations of the {\it off-diagonal} energy and norm kernels that lie at the heart of the nuclear EDF method. The many-body formalism of interest relies on the concepts of symmetry breaking {\it and} restoration that have made the fortune of the nuclear EDF method and is, as such, amenable to this guidance. After elaborating on our proposal, we briefly outline the project we plan to execute in the years to come.

preprint2015arXiv

Beyond mean-field study of elastic and inelastic electron scattering off nuclei

Electron scattering provides a powerful tool to determine charge distributions and transition densities of nuclei. This tool will soon be available for short-lived neutron-rich nuclei. [Purpose] Beyond mean-field methods have been successfully applied to the study of excitation spectra of nuclei in the whole nuclear chart. These methods permit to determine energies and transition probabilities starting from an effective in-medium nucleon-nucleon interaction but without other phenomenological ingredients. Such a method has recently been extended to calculate the charge density of nuclei deformed at the mean-field level of approximation [J. M. Yao et al., Phys. Rev. C86, 014310 (2012)]. The aim of this work is to further extend the method to the determination of transition densities between low-lying excited states. [Method] The starting point of our method is a set of Hartree-Fock-Bogoliubov wave functions generated with a constraint on the axial quadrupole moment and using a Skyrme energy density functional. Correlations beyond the mean field are introduced by projecting mean-field wave functions on angular-momentum and particle number and by mixing the symmetry restored wave functions.[Results] We give in this paper detailed formulae derived for the calculation of densities and form factors. These formulae are rather easy to obtain when both initial and final states are $0^+$ states but are far from being trivial when one of the states has a finite $J$-value. Illustrative applications to $^{24}$Mg and to the even-mass $^{58-68}$Ni have permitted to analyse the main features of our method, in particular the effect of deformation on densities and form factors. An illustration calculation of both elastic and inelastic scattering form factors is presented....

preprint2015arXiv

Collectivity in the light radon nuclei measured directly via Coulomb excitation

Background: Shape coexistence in heavy nuclei poses a strong challenge to state-of-the-art nuclear models, where several competing shape minima are found close to the ground state. A classic region for investigating this phenomenon is in the region around $Z=82$ and the neutron mid-shell at $N=104$. Purpose: Evidence for shape coexistence has been inferred from $α$-decay measurements, laser spectroscopy and in-beam measurements. While the latter allow the pattern of excited states and rotational band structures to be mapped out, a detailed understanding of shape coexistence can only come from measurements of electromagnetic matrix elements. Method: Secondary, radioactive ion beams of $^{202}$Rn and $^{204}$Rn were studied by means of low-energy Coulomb excitation at the REX-ISOLDE facility in CERN. Results: The electric-quadrupole ($E2$) matrix element connecting the ground state and first-excited $2^{+}_{1}$ state was extracted for both $^{202}$Rn and $^{204}$Rn, corresponding to ${B(E2;2^{+}_{1} \to 2^{+}_{1})=29^{+8}_{-8}}$ W.u. and $43^{+17}_{-12}$ W.u., respectively. Additionally, $E2$ matrix elements connecting the $2^{+}_{1}$ state with the $4^{+}_{1}$ and $2^{+}_{2}$ states were determined in $^{202}$Rn. No excited $0^{+}$ states were observed in the current data set, possibly due to a limited population of second-order processes at the currently-available beam energies. Conclusions: The results are discussed in terms of collectivity and the deformation of both nuclei studied is deduced to be weak, as expected from the low-lying level-energy schemes. Comparisons are also made to state-of-the-art beyond-mean-field model calculations and the magnitude of the transitional quadrupole moments are well reproduced.

preprint2015arXiv

Numerical accuracy of mean-field calculations in coordinate space

Background: Mean-field methods based on an energy density functional (EDF) are powerful tools used to describe many properties of nuclei in the entirety of the nuclear chart. The accuracy required on energies for nuclear physics and astrophysics applications is of the order of 500 keV and much effort is undertaken to build EDFs that meet this requirement. Purpose: The mean-field calculations have to be accurate enough in order to preserve the accuracy of the EDF. We study this numerical accuracy in detail for a specific numerical choice of representation for the mean-field equations that can accommodate any kind of symmetry breaking. Method: The method that we use is a particular implementation of 3-dimensional mesh calculations. Its numerical accuracy is governed by three main factors: the size of the box in which the nucleus is confined, the way numerical derivatives are calculated and the distance between the points on the mesh. Results: We have examined the dependence of the results on these three factors for spherical doubly-magic nuclei, neutron-rich $^{34}$Ne, the fission barrier of $^{240}$Pu and isotopic chains around Z = 50. Conclusions: Mesh calculations offer the user extensive control over the numerical accuracy of the solution scheme. By making appropriate choices for the numerical scheme the achievable accuracy is well below the model uncertainties of mean-field methods.

preprint2015arXiv

Properties of nuclei in the nobelium region studied within the covariant, Skyrme, and Gogny energy density functionals

We calculate properties of the ground and excited states of nuclei in the nobelium region for proton and neutron numbers of 92 <= Z <= 104 and 144 <= N <= 156, respectively. We use three different energy-density-functional (EDF) approaches, based on covariant, Skyrme, and Gogny functionals, each within two different parameter sets. A comparative analysis of the results obtained for odd-even mass staggerings, quasiparticle spectra, and moments of inertia allows us to identify single-particle and shell effects that are characteristic to these different models and to illustrate possible systematic uncertainties related to using the EDF modelling

preprint2014arXiv

Beyond Mean-Field Calculations for Odd-A Nuclei

Beyond mean-field methods are very successful tools for the description of large-amplitude collective motion for even-even atomic nuclei. The state-of-the-art framework of these methods consists in a Generator Coordinate Method based on angular-momentum and particle-number projected triaxially deformed Hatree-Fock-Bogoliubov (HFB) states. The extension of this scheme to odd-mass nuclei is a long-standing challenge. We present for the first time such an extension, where the Generator Coordinate space is built from self-consistently blocked one-quasiparticle HFB states. One of the key points for this success is that the same Skyrme interaction is used for the mean-field and the pairing channels, thus avoiding problems related to the violation of the Pauli principle. An application to 25Mg illustrates the power of our method, as agreement with experiment is obtained for the spectrum, electromagnetic moments, and transition strengths, for both positive and negative parity states and without the necessity for effective charges or effective moments. Although the effective interaction still requires improvement, our study opens the way to systematically describe odd-A nuclei throughout the nuclear chart.

preprint2014arXiv

Local formation of nitrogen-vacancy centers in diamond by swift heavy ions

We exposed nitrogen-implanted diamonds to beams of swift uranium and gold ions (~1 GeV) and find that these irradiations lead directly to the formation of nitrogen vacancy (NV) centers, without thermal annealing. We compare the photoluminescence intensities of swift heavy ion activated NV- centers to those formed by irradiation with low-energy electrons and by thermal annealing. NV- yields from irradiations with swift heavy ions are 0.1 of yields from low energy electrons and 0.02 of yields from thermal annealing. We discuss possible mechanisms of NV-center formation by swift heavy ions such as electronic excitations and thermal spikes. While forming NV centers with low efficiency, swift heavy ions enable the formation of three dimensional NV- assemblies over relatively large distances of tens of micrometers. Further, our results show that NV-center formation is a local probe of (partial) lattice damage relaxation induced by electronic excitations from swift heavy ions in diamond.

preprint2014arXiv

Solution of the Skyrme HF+BCS equation on a 3D mesh. II. A new version of the Ev8 code

We describe a new version of the EV8 code that solves the nuclear Skyrme-Hartree-Fock+BCS problem using a 3-dimensional cartesian mesh. Several new features have been implemented with respect to the earlier version published in 2005. In particular, the numerical accuracy has been improved for a given mesh size by (i) implementing a new solver to determine the Coulomb potential for protons (ii) implementing a more precise method to calculate the derivatives on a mesh that had already been implemented earlier in our beyond-mean-field codes. The code has been made very flexible to enable the use of a large variety of Skyrme energy density functionals that have been introduced in the last years. Finally, the treatment of the constraints that can be introduced in the mean-field equations has been improved. The code Ev8 is today the tool of choice to study the variation of the energy of a nucleus from its ground state to very elongated or triaxial deformations with a well-controlled accuracy.

preprint2013arXiv

Radiation Hardness of Graphene and MoS2 Field Effect Devices Against Swift Heavy Ion Irradiation

We have investigated the deterioration of field effect transistors based on twodimensional materials due to irradiation with swift heavy ions. Devices were prepared with exfoliated single layers of MoS2 and graphene, respectively. They were characterized before and after irradiation with 1.14 GeV U228+2 ions using three different fluences. By electrical characterization, atomic force microscopy and Raman spectroscopy we show that the irradiation leads to significant changes of structural and electrical properties. At the highest fluence of 4 x 102^11 ions/cm^2, the MoS2 transistor is destroyed, while the graphene based device remains operational, albeit with an inferior performance.

preprint2013arXiv

Skyrme functional from a three-body pseudo-potential of second-order in gradients. Formalism for central terms

In one way or the other, all modern parametrizations of the nuclear energy density functional (EDF) do not respect the exchange symmetry associated with Pauli's principle. It has been recently shown that this practice jeopardizes multi-reference (MR) EDF calculations by contaminating the energy with spurious self-interactions that, for example, lead to finite steps or even divergences when plotting it as a function of collective coordinates. As of today, the only viable option to bypass these pathologies is to rely on EDF kernels that enforce Pauli's principle from the outset by strictly and exactly deriving from a genuine, i.e. density-independent, Hamilton operator. We wish to develop the most general Skyrme-like EDF parametrization containing linear, bilinear and trilinear terms in the density matrices with up to two gradients, under the key constraint that it derives strictly from an effective Hamilton operator. The most general three-body Skyrme-like pseudo-potential containing up to two gradient operators is constructed to generate the trilinear part. The present study is limited to central terms. Spin-orbit and tensor will be addressed in a forthcoming paper. (See paper for full abstract)

preprint2013arXiv

Spurious finite-size instabilities in nuclear energy density functionals

It is known that some well-established parametrizations of the EDF do not always provide converged results for nuclei and a qualitative link between this finding and the appearance of finite-size instabilities of SNM near saturation density when computed within the RPA has been pointed out. We seek for a quantitative and systematic connection between the impossibility to converge self-consistent calculations of nuclei and the occurrence of finite-size instabilities in SNM for the example of scalar-isovector (S=0, T=1) instabilities of the standard Skyrme EDF. We aim to establish a stability criterion based on computationally-friendly RPA calculations of SNM that is independent on the functional form of the EDF and that can be utilized during the adjustment of its coupling constants. Tuning the coupling constant $C^{ρΔρ}_{1}$ of the gradient term that triggers scalar-isovector instabilities of the standard Skyrme EDF, we find that the occurrence of instabilities in finite nuclei depends strongly on the numerical scheme used to solve the self-consistent mean-field equations. The link to instabilities of SNM is made by extracting the lowest density $ρ_{\text{crit}}$ at which a pole appears at zero energy in the RPA response function when employing the critical value of the coupling constant $C^{ρΔρ}_{1}$ extracted in nuclei. Our analysis suggests a two-fold stability criterion to avoid scalar-isovector instabilities: (i) The density $ρ_{\text{min}}$ corresponding to the lowest pole in the RPA response function should be larger than about 1.2 times the saturation density; (ii) one needs to verify that $ρ_{p}(q_{\text{pq}})$ exhibits a distinct global minimum and is not a decreasing function for large transferred momenta.

preprint2012arXiv

New parameterization of Skyrme's interaction for regularized multi-reference energy density functional calculations

[Background] Symmetry restoration and configuration mixing in the spirit of the generator coordinate method based on energy density functionals have become widely used techniques in low-energy nuclear structure physics. Recently, it has been pointed out that these techniques are ill-defined for standard Skyrme functionals, and a regularization procedure has been proposed to remove the resulting spuriosities from such calculations. This procedure imposes an integer power of the density for the density dependent terms of the functional. At present, only dated parameterizations of the Skyrme interaction fulfill this condition. [Purpose] To construct a set of parameterizations of the Skyrme energy density functional for multi-reference energy density functional calculations with regularization using the state-of-the-art fitting protocols. [Method] The parameterizations were adjusted to reproduce ground state properties of a selected set of doubly magic nuclei and properties of nuclear matter. Subsequently, these parameter sets were validated against properties of spherical and deformed nuclei. [Results] Our parameter sets successfully reproduce the experimental binding energies and charge radii for a wide range of singly-magic nuclei. Compared to the widely used SLy5 and to the SIII parameterization that has integer powers of the density, a significant improvement of the reproduction of the data is observed. Similarly, a good description of the deformation properties at $A\sim 80$ was obtained. [Conclusions] We have constructed new Skyrme parameterizations with integer powers of the density and validated them against a broad set of experimental data for spherical and deformed nuclei. These parameterizations are tailor-made for regularized multi-reference energy density functional calculations and can be used to study correlations beyond the mean-field in atomic nuclei.

preprint2012arXiv

Skyrme pseudo-potential-based EDF parametrisation for spuriousity-free MR-EDF calculations

First exploratory steps towards a pseudo-potential-based Skyrme energy density functional for spuriousity-free multi-reference calculations are presented. A qualitatively acceptable fit can be accomplished by adding simple three- and four-body contact terms to the standard central plus spin-orbit two-body terms. To achieve quantitative predictive power, higher-order terms, e.g. velocity-dependent three-body terms, will be required.

preprint2012arXiv

Systematics of low-lying states of even-even nuclei in the neutron-deficient lead region from a beyond-mean-field calculation

Nuclei located in the neutron-deficient Pb region have a complex structure, rapidly evolving as a function of neutron and proton numbers. The most famous example is $^{186}$Pb where the three lowest levels are $0^+$ states, the two excited $0^+$ ones being located at low excitation energy around 600 keV. Coexisting structures with different properties are found in the neighboring nuclei. Many experiments have been performed over the last few years in which in-band and out-of-band $γ$ transition probabilities have been measured. A detailed interpretation of experimental data requires the use of a method going beyond a mean-field approach that permits to determine spectra and transition probabilities. Such methods have already been applied to selected isotopes in this mass region. Our aim is to provide a systematic investigation of this mass region in order to determine how well experimental data can be understood using a state-of-the-art method for nuclear structure. The starting point of our method is a set of mean-field wave functions generated with a constraint on the axial quadrupole moment and using a Skyrme energy density functional. Correlations beyond the mean field are introduced by projecting mean-field wave functions on angular-momentum and particle number and by mixing the symmetry restored wave functions as a function of the axial quadrupole moment. A detailed comparison with the available data is performed for energies, charge radii, spectroscopic quadrupole moments, $E0$ and $E2$ transition probabilities for the isotopic chains of neutron deficient Hg, Pb, Po and Rn. The connection between our results and the underlying mean field is also analyzed. Qualitative agreement with the data is obtained although our results indicate that the actual energy density functionals have to be improved further to achieve a quantitative agreement.

preprint2011arXiv

Symmetry restoration for odd-mass nuclei with a Skyrme energy density functional

In these proceedings, we report first results for particle-number and angular-momentum projection of self-consistently blocked triaxial one-quasiparticle HFB states for the description of odd-A nuclei in the context of regularized multi-reference energy density functionals, using the entire model space of occupied single-particle states. The SIII parameterization of the Skyrme energy functional and a volume-type pairing interaction are used.

preprint2011arXiv

The tensor part of the Skyrme energy density functional. III. Time-odd terms at high spin

This article extends previous studies on the effect of tensor terms in the Skyrme energy density functional by breaking of time-reversal invariance. We have systematically probed the impact of tensor terms on properties of superdeformed rotational bands calculated within the cranked Hartree-Fock-Bogoliubov approach for different parameterizations covering a wide range of values for the isoscalar and isovector tensor coupling constants. We analyze in detail the contribution of the tensor terms to the energies and dynamical moments of inertia and study their impact on quasi-particle spectra. Special attention is devoted to the time-odd tensor terms, the effect of variations of their coupling constants and finite-size instabilities.

preprint2010arXiv

Early onset of ground-state deformation in the neutron-deficient polonium isotopes

In-source resonant ionization laser spectroscopy of the even-$A$ polonium isotopes $^{192-210,216,218}$Po has been performed using the $6p^37s$ $^5S_2$ to $6p^37p$ $^5P_2$ ($λ=843.38$ nm) transition in the polonium atom (Po-I) at the CERN ISOLDE facility. The comparison of the measured isotope shifts in $^{200-210}$Po with a previous data set allows to test for the first time recent large-scale atomic calculations that are essential to extract the changes in the mean-square charge radius of the atomic nucleus. When going to lighter masses, a surprisingly large and early departure from sphericity is observed, which is only partly reproduced by Beyond Mean Field calculations.

preprint2010arXiv

Regularization of Multi-Reference Energy Density Functional Calculations

We report the first application of a recently proposed regularization procedure for multi-reference energy density functionals, which removes spurious divergent or non-continuous contributions to the binding energy, to a general configuration mixing. As an example, we present a calculation that corresponds to the particle-number and angular momentum projection of axially symmetric time-reversal invariant quasiparticle vacua of different quadrupole deformation for the nucleus O18. The SIII parameterization of the Skyrme energy functional is used.

preprint2010arXiv

Shape coexistence at the proton drip-line: First identification of excited states in 180Pb

Excited states in the extremely neutron-deficient nucleus, 180Pb, have been identified for the first time using the JUROGAM II array in conjunction with the RITU recoil separator at the Accelerator Laboratory of the University of Jyvaskyla. This study lies at the limit of what is presently achievable with in-beam spectroscopy, with an estimated cross-section of only 10 nb for the 92Mo(90Zr,2n)180Pb reaction. A continuation of the trend observed in 182Pb and 184Pb is seen, where the prolate minimum continues to rise beyond the N=104 mid-shell with respect to the spherical ground state. Beyond mean-field calculations are in reasonable correspondence with the trends deduced from experiment.

preprint2010arXiv

Surface Symmetry Energy of Nuclear Energy Density Functionals

We study the bulk deformation properties of the Skyrme nuclear energy density functionals. Following simple arguments based on the leptodermous expansion and liquid drop model, we apply the nuclear density functional theory to assess the role of the surface symmetry energy in nuclei. To this end, we validate the commonly used functional parametrizations against the data on excitation energies of superdeformed band-heads in Hg and Pb isotopes, and fission isomers in actinide nuclei. After subtracting shell effects, the results of our self-consistent calculations are consistent with macroscopic arguments and indicate that experimental data on strongly deformed configurations in neutron-rich nuclei are essential for optimizing future nuclear energy density functionals. The resulting survey provides a useful benchmark for further theoretical improvements. Unlike in nuclei close to the stability valley, whose macroscopic deformability hangs on the balance of surface and Coulomb terms, the deformability of neutron-rich nuclei strongly depends on the surface-symmetry energy; hence, its proper determination is crucial for the stability of deformed phases of the neutron- rich matter and description of fission rates for r-process nucleosynthesis.

preprint2009arXiv

Configuration Mixing within the Energy Density Functional Formalism: Removing Spurious Contributions from Non-Diagonal Energy Kernels

Multi-reference calculations along the lines of the Generator Coordinate Method or the restoration of broken symmetries within the nuclear Energy Density Functional (EDF) framework are becoming a standard tool in nuclear structure physics. These calculations rely on the extension of a single-reference energy functional, of the Gogny or the Skyrme types, to non-diagonal energy kernels. There is no rigorous constructive framework for this extension so far. The commonly accepted way proceeds by formal analogy with the expressions obtained when applying the generalized Wick theorem to the non-diagonal matrix element of a Hamilton operator between two product states. It is pointed out that this procedure is ill-defined when extended to EDF calculations as the generalized Wick theorem is taken outside of its range of applicability. In particular, such a procedure is responsible for the appearance of spurious divergences and steps in multi-reference EDF energies, as was recently observed in calculations restoring particle number or angular momentum. In the present work, we give a formal analysis of the origin of this problem for calculations with and without pairing, i.e. constructing the density matrices from either Slater determinants or quasi-particle vacua. We propose a correction to energy kernels that removes the divergences and steps, and which is applicable to calculations based on any symmetry restoration or generator coordinate. The method is formally illustrated for particle number restoration and is specified to configuration mixing calculations based on Slater determinants.

preprint2009arXiv

Particle-Number Restoration within the Energy Density Functional Formalism

We give a detailed analysis of the origin of spurious divergences and finite steps that have been recently identified in particle-number restoration calculations within the nuclear energy density functional framework. We isolate two distinct levels of spurious contributions to the energy. The first one is encoded in the definition of the basic energy density functional itself whereas the second one relates to the canonical procedure followed to extend the use of the energy density functional to multi-reference calculations. The first level of spuriosity relates to the long-known self-interaction problem and to the newly discussed self-pairing interaction process which might appear when describing paired systems with energy functional methods using auxiliary reference states of Bogoliubov or BCS type. A minimal correction to the second level of spuriosity to the multi-reference nuclear energy density functional proposed in [D. Lacroix, T. Duguet, M. Bender, arXiv:0809.2041] is shown to remove completely the anomalies encountered in particle-number restored calculations. In particular, it restores sum-rules over (positive) particle numbers that are to be fulfilled by the particle-number-restored formalism. The correction is found to be on the order of several hundreds of keVs up to about 1 MeV in realistic calculations, which is small compared to the total binding energy, but often accounts for a substantial percentage of the energy gain from particle-number restoration and is on the same energy scale as the excitations one addresses with multi-reference energy density functional methods.

preprint2009arXiv

Particle-Number Restoration within the Energy Density Functional formalism: Nonviability of terms depending on noninteger powers of the density matrices

We discuss the origin of pathological behaviors that have been recently identified in particle-number-restoration calculations performed within the nuclear energy density functional framework. A regularization method that removes the problematic terms from the multi-reference energy density functional and which applies (i) to any symmetry restoration- and/or generator-coordinate-method-based configuration mixing calculation and (ii) to energy density functionals depending only on integer powers of the density matrices, was proposed in [D. Lacroix, T. Duguet, M. Bender, arXiv:0809.2041] and implemented for particle-number restoration calculations in [M. Bender, T. Duguet, D. Lacroix, arXiv:0809.2045]. In the present paper, we address the viability of non-integer powers of the density matrices in the nuclear energy density functional. Our discussion builds upon the analysis already carried out in [J. Dobaczewski \emph{et al.}, Phys. Rev. C \textbf{76}, 054315 (2007)]. First, we propose to reduce the pathological nature of terms depending on a non-integer power of the density matrices by regularizing the fraction that relates to the integer part of the exponent using the method proposed in [D. Lacroix, T. Duguet, M. Bender, arXiv:0809.2041]. Then, we discuss the spurious features brought about by the remaining fractional power. Finally, we conclude that non-integer powers of the density matrices are not viable and should be avoided in the first place when constructing nuclear energy density functionals that are eventually meant to be used in multi-reference calculations.

preprint2007arXiv

The tensor part of the Skyrme energy density functional. I. Spherical nuclei

We perform a systematic study of the impact of the J^2 tensor term in the Skyrme energy functional on properties of spherical nuclei. In the Skyrme energy functional, the tensor terms originate both from zero-range central and tensor forces. We build a set of 36 parameterizations, which covers a wide range of the parameter space of the isoscalar and isovector tensor term coupling constants, with a fit protocol very similar to that of the successful SLy parameterizations. We analyze the impact of the tensor terms on a large variety of observables in spherical mean-field calculations, such as the spin-orbit splittings and single-particle spectra of doubly-magic nuclei, the evolution of spin-orbit splittings along chains of semi-magic nuclei, mass residuals of spherical nuclei, and known anomalies of charge radii. Our main conclusion is that the currently used central and spin-orbit parts of the Skyrme energy density functional are not flexible enough to allow for the presence of large tensor terms.

preprint2004arXiv

Self-Consistent Description of Multipole Strength in Exotic Nuclei I: Method

We use the canonical Hartree-Fock-Bogoliubov basis to implement a completely self-consistent quasiparticle-random-phase approximation with arbitrary Skyrme energy density functionals and density-dependent pairing functionals. The point of the approach is to accurately describe multipole strength functions in spherical even-even nuclei, including weakly-bound drip-line systems. We describe the method and carefully test its accuracy, particularly in handling spurious modes. To illustrate our approach, we calculate isoscalar and isovector monopole, dipole, and quadrupole strength functions in several Sn isotopes, both in the stable region and at the drip lines.

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.

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.