Source author record

Paul-Gerhard Reinhard

Paul-Gerhard 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

14works
10topics
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

14 published item(s)

preprint2023arXiv

Combined theoretical analysis of the parity-violating asymmetry for ${}^{48}$Ca and ${}^{208}Pb$

The recent experimental determination of the parity violating asymmetry $A_{\rm pv}$ in ${}^{48}$Ca and ${}^{208}$Pb at Jefferson Lab is important for our understanding on how neutrons and protons arrange themselves inside the atomic nucleus. To better understand the impact of these measurements, we present a rigorous theoretical investigation of $A_{\rm pv}$ in ${}^{48}$Ca and ${}^{208}$Pb and assess the associated uncertainties. We complement our study by inspecting the static electric dipole polarizability in these nuclei. The analysis is carried out within nuclear energy density functional theory with quantified input. We conclude that the simultaneous accurate description of $A_{\rm pv}$ in ${}^{48}$Ca and ${}^{208}$Pb cannot be achieved by our models that accommodate a pool of global nuclear properties, such as masses and charge radii, throughout the nuclear chart, and describe -- within one standard deviation -- the experimental dipole polarizabilities $α_{\rm D}$ in these nuclei.

preprint2022arXiv

Evidence of Two-Source King Plot Nonlinearity in Spectroscopic Search for New Boson

Optical precision spectroscopy of isotope shifts can be used to test for new forces beyond the Standard Model, and to determine basic properties of atomic nuclei. We measure isotope shifts on the highly forbidden ${}^2S_{1/2} \rightarrow {}^2F_{7/2}$ octupole transition of trapped $^{168,170,172,174,176}$Yb ions. When combined with previous measurements in Yb$^+$ and very recent measurements in Yb, the data reveal a King plot nonlinearity of up to 240$σ$. The trends exhibited by experimental data are explained by nuclear density functional theory calculations with the Fayans functional. We also find, with 4.3$σ$ confidence, that there is a second distinct source of nonlinearity, and discuss its possible origin.

preprint2022arXiv

Information content of the differences in the charge radii of mirror nuclei

Differences in the charge radii of mirror nuclei have been recently suggested to contain information on the slope of the symmetry energy L. To test this hypothesis, we perform statistical correlation analysis using quantified energy density functionals that are consistent with our previous knowledge on global nuclear observables such as binding energies and charge radii. We conclude that the difference in charge radii between a mirror pair, R_mir, is an inferior isovector indicator compared to other observables, such at the neutron skin or electric dipole polarizability. In particular, this quantity correlates poorly with both the neutron skin and L. We demonstrate that R_mir is influenced by pairing correlations in the presence of low-lying proton continuum in the proton-rich mirror-partner nucleus. Considering the large theoretical uncertainties on R_mir, we conclude that the precise data on mirror charge radii cannot provide a stringent constraint on L.

preprint2022arXiv

Statistical correlations of nuclear quadrupole deformations and charge radii

Shape deformations and charge radii, basic properties of atomic nuclei, are influenced by both the global features of the nuclear force and the nucleonic shell structure. As functions of proton and neutron number, both quantities show regular patterns and, for nuclei away from magic numbers, they change very smoothly from nucleus to nucleus. In this paper, we explain how the local shell effects are impacting the statistical correlations between quadrupole deformations and charge radii in well-deformed even-even Er, Yb, and Hf isotopes. This implies, in turn, that sudden changes in correlations can be useful indicators of underlying shell effects. Our theoretical analysis is performed in the framework of self-consistent mean-field theory using quantified energy density functionals and density-dependent pairing forces. The statistical analysis is carried out by means of the linear least-square regression. The local variations of nuclear quadrupole deformations and charge radii, explained in terms of occupations individual deformed Hartree-Fock orbits, make and imprint on statistical correlations of computed observables. While the calculated deformations or charge radii are, in some cases, correlated with those of their even-even neighbors, the correlations seem to deteriorate rapidly with particle number. The statistical correlations between nuclear deformations and charge radii of different nuclei are affected by the underlying shell structure. Even for well deformed and superfluid nuclei for which these observables change smoothly, the correlation range usually does not exceed $ΔN=4$ and $ΔZ=4$, i.e., it is rather short. This result suggests that the frequently made assumption of reduced statistical errors for the differences between smoothly-varying observables cannot be generally justified.

preprint2022arXiv

Three-dimensional Skyrme Hartree-Fock-Bogoliubov solver in coordinate-space representation

The coordinate-space representation of the Hartree-Fock-Bogoliubov theory is the method of choice to study weakly bound nuclei whose properties are affected by the quasiparticle continuum space. To describe such systems, we developed a three-dimensional Skyrme-Hartree-Fock-Bogoliubov solver HFBFFT based on the existing, highly optimized and parallelized Skyrme-Hartree-Fock code Sky3D. The code does not impose any self-consistent spatial symmetries such as mirror inversions or parity. The underlying equations are solved in HFBFFT directly in the canonical basis using the fast Fourier transform. To remedy the problems with pairing collapse, we implemented the soft energy cutoff and pairing annealing. The convergence of HFB solutions was improved by a sub-iteration method. The Hermiticity violation of differential operators brought by Fourier-transform-based differentiation has also been solved. The accuracy and performance of HFBFFT were tested by benchmarking it against other HFB codes, both spherical and deformed, for a set of nuclei, both well-bound and weakly-bound.

preprint2021arXiv

Universal trend of charge radii of even-even Ca-Zn nuclei

Radii of nuclear charge distributions carry information about the strong and electromagnetic forces acting inside the atomic nucleus. While the global behavior of nuclear charge radii is governed by the bulk properties of nuclear matter, their local trends are affected by quantum motion of proton and neutron nuclear constituents. The measured differential charge radii $δ\langle r^2_c\rangle$ between neutron numbers $N=28$ and $N=40$ exhibit a universal pattern as a function of $n=N-28$ that is independent of the atomic number. Here we analyze this remarkable behavior in even-even nuclei from calcium to zinc using two state-of-the-art theories based on quantified nuclear interactions: the ab-initio coupled cluster theory and nuclear density functional theory. Both theories reproduce the smooth rise of differential charge radii and their weak dependence on the atomic number. By considering a large set of isotopic chains, we show that this trend can be captured by just two parameters: the slope and curvature of ${δ\langle r^2_c\rangle(n)}$. We demonstrate that these parameters show appreciable model dependence, and the statistical analysis indicates that they are not correlated with any single model property, i.e., they are impacted by both bulk nuclear properties as well as shell structure.

preprint2020arXiv

Rate for Laser-Induced Nuclear Dipole Absorption

Using the Brink-Axel hypothesis we derive the rate $R$ for nuclear dipole excitation by a laser pulse carrying $N \gg 1$ photons with average energy $\hbar ω_0 \approx 5$ MeV. As expected $R \propto (\hbar ω_0)^3$. The rate is also proportional to the aperure $α$ of the laser pulse. Perhaps less expected is the fact that $R \propto N$, irrespective of the degree of coherence of the laser pulse. The expression for $R$, derived for a nearly stationary laser pulse, is valid also for short times and can, thus, be used in simulations via rate equations of multiple nuclear dipole excitations by a single pulse. The explicit dependence of $R$ on the parameters of the laser pulse and on nuclear parameters given in the paper should help to optimize experiments on laser-nucleus reactions.

preprint2020arXiv

Statistical aspects of nuclear mass models

We study the information content of nuclear masses from the perspective of global models of nuclear binding energies. To this end, we employ a number of statistical methods and diagnostic tools, including Bayesian calibration, Bayesian model averaging, chi-square correlation analysis, principal component analysis, and empirical coverage probability. Using a Bayesian framework, we investigate the structure of the 4-parameter Liquid Drop Model by considering discrepant mass domains for calibration. We then use the chi-square correlation framework to analyze the 14-parameter Skyrme energy density functional calibrated using homogeneous and heterogeneous datasets. We show that a quite dramatic parameter reduction can be achieved in both cases. The advantage of Bayesian model averaging for improving uncertainty quantification is demonstrated. The statistical approaches used are pedagogically described; in this context this work can serve as a guide for future applications.

preprint2016arXiv

Laser-Matter Interaction: Classical Regime versus Quantum Regime

Doppler backscattering of optical laser photons on a "flying mirror" of relativistic electrons promises to yield coherent photons with MeV-range energies. We compare the nuclear interaction of such a laser pulse with the standard atom-laser interaction. The mean-field description of atoms must be replaced by a rate equation and the classical field strength, far too faint in nuclei, by the dipole transition rate. Significant nuclear excitation occurs for photon numbers much smaller than typical for atoms. That drastically reduces the requirements on the experimental realization of a "flying mirror".

preprint2016arXiv

Stochastic TDHF in an exactly solvable model

We apply in a schematic model a theory beyond mean-field, namely Stochastic Time-Dependent Hartree-Fock (STDHF), which includes dynamical electron-electron collisions on top of an incoherent ensemble of mean-field states by occasional 2-particle-2-hole ($2p2h$) jumps. The model considered here is inspired by a Lipkin-Meshkov-Glick model of $Ω$ particles distributed into two bands of energy and coupled by a two-body interaction. Such a model can be exactly solved (numerically though) for small $Ω$. It therefore allows a direct comparison of STDHF and the exact propagation. The systematic impact of the model parameters as the density of states, the excitation energy and the bandwidth is presented and discussed. The time evolution of the STDHF compares fairly well with the exact entropy, as soon as the excitation energy is sufficiently large to allow $2p2h$ transitions. Limitations concerning low energy excitations and memory effects are also discussed.

preprint2016arXiv

Strong-field effects in the photo-emission spectrum of the C$_{60}$ fullerene

Considering C$_{60}$ as a model system for describing field emission from the extremity of a carbon nanotip, we explore electron emission from this fullerene excited by an intense, near-infrared, few-cycle laser pulse ($10^{13}$-$10^{14}~{\rm W/cm}^2$, 912 nm, 8-cycle). To this end, we use time-dependent density functional theory augmented by a self-interaction correction. The ionic background of C$_{60}$ is described by a soft jellium model. Particular attention is paid to the high energy electrons. Comparing the spectra at different emission angles, we find that, as a major result of this study, the photoelectrons are strongly peaked along the laser polarization axis forming a highly collimated electron beam in the forward direction, especially for the high energy electrons. Moreover, the high-energy plateau cut-off found in the simulations agrees well with estimates from the classical three-step model. We also investigate the build-up of the high-energy part of a photoelectron spectrum by a time-resolved analysis. In particular, the modulation on the plateau can be interpreted as contributions from intracycle and intercycle interferences.

preprint2014arXiv

Optimizing relativistic energy density functionals: covariance analysis

The stability of model parameters for a class of relativistic energy density functionals, characterized by contact (point-coupling) effective inter-nucleon interactions and density-dependent coupling parameters, is analyzed using methods of statistical analysis. A set of pseudo-observables in infinite and semi-infinite nuclear matter is used to define a quality measure $χ^2$ for subsequent analysis. We calculate uncertainties of model parameters and correlation coefficients between parameters, and determine the eigenvectors and eigenvalues of the matrix of second derivatives of $χ^2$ at the minimum. This allows to examine the stability of the density functional in nuclear matter, and to deduce weakly and strongly constrained combinations of parameters. In addition, we also compute uncertainties of observables that are not included in the calculation of $χ^2$: binding energy of asymmetric nuclear matter, surface thickness of semi-infinite nuclear matter, binding energies and charge radii of finite nuclei.

preprint2013arXiv

Koopmans' condition in self-interaction corrected density functional theory

We investigate from a practitioner's point of view the computation of the ionization potential (IP) within density functional theory (DFT). DFT with (semi-)local energy-density functionals is plagued by a self-interaction error which hampers the computation of IP from the single-particle energy of the highest occupied molecular orbital (HOMO). The problem may be cured by a self interaction correction (SIC) for which there exist various approximate treatments. We compare the performance of the SIC proposed by Perdew and Zunger with the very simple average-density SIC (ADSIC) for a large variety of atoms and molecules up to larger systems as carbon rings and chains. Both approaches to SIC provide a large improvement to the quality of the IP if calculated from the HOMO level. The surprising result is that the simple ADSIC performs even better than the original Perdew-Zunger SIC (PZSIC) in the majority of the studied cases.

preprint2013arXiv

Time-Dependent Hartree-Fock Approach to Nuclear Pasta at Finite Temperature

We present simulations of neutron-rich matter at subnuclear densities, like supernova matter, with the time-dependent Hartree-Fock approximation at temperatures of several MeV. The initial state consists of $α$ particles randomly distributed in space that have a Maxwell-Boltzmann distribution in momentum space. Adding a neutron background initialized with Fermi distributed plane waves the calculations reflect a reasonable approximation of astrophysical matter. This matter evolves into spherical, rod-like, and slab-like shapes and mixtures thereof. The simulations employ a full Skyrme interaction in a periodic three-dimensional grid. By an improved morphological analysis based on Minkowski functionals, all eight pasta shapes can be uniquely identified by the sign of only two valuations, namely the Euler characteristic and the integral mean curvature. In addition, we propose the variance in the cell density distribution as a measure to distinguish pasta matter from uniform matter.