Combinatorial nuclear level density by a Monte Carlo method

We present a new combinatorial method for the calculation of the nuclear level density. It is based on a Monte Carlo technique, in order to avoid a direct counting procedure which is generally impracticable for high-A nuclei. The Monte Carlo simulation, making use of the Metropolis sampling scheme, allows a computationally fast estimate of the level density for many fermion systems in large shell model spaces. We emphasize the advantages of this Monte Carlo approach, particularly concerning the prediction of the spin and parity distributions of the excited states, and compare our results with those derived from a traditional combinatorial or a statistical method. Such a Monte Carlo technique seems very promising to determine accurate level densities in a large energy range for nuclear reaction calculations.Comment: 30 pages, LaTex, 7 figures (6 Postscript figures included). Fig. 6 upon request to the autho

Recent breakthroughs in Skyrme-Hartree-Fock-Bogoliubov mass formulas

We review our recent achievements in the construction of microscopic mass tables based on the Hartree-Fock-Bogoliubov method with Skyrme effective interactions. In the latest of our series of HFB-mass models, we have obtained our best fit ever to essentially all the available mass data, by treating the pairing more realistically than in any of our earlier models. The rms deviation on the 2149 measured masses of nuclei with N and Z>8 has been reduced for the first time in a mean field approach to 0.581 MeV. With the additional constraint on the neutron-matter equation of state, this new force is thus very well-suited for the study of neutron-rich nuclei and for the description of astrophysical environments like supernova cores and neutron-star crusts.Comment: Proceedings of the Fifth International Conference on Exotic Nuclei and Atomic Masses, September 7-13 2008, Ryn (Poland). To appear in the European Physical Journal

Shell structure of superheavy nuclei in self-consistent mean-field models

We study the extrapolation of nuclear shell structure to the region of superheavy nuclei in self-consistent mean-field models -- the Skyrme-Hartree-Fock approach and the relativistic mean-field model -- using a large number of parameterizations. Results obtained with the Folded-Yukawa potential are shown for comparison. We focus on differences in the isospin dependence of the spin-orbit interaction and the effective mass between the models and their influence on single-particle spectra. While all relativistic models give a reasonable description of spin-orbit splittings, all non-relativistic models show a wrong trend with mass number. The spin-orbit splitting of heavy nuclei might be overestimated by 40%-80%. Spherical doubly-magic superheavy nuclei are found at (Z=114,N=184), (Z=120,N=172) or (Z=126,N=184) depending on the parameterization. The Z=114 proton shell closure, which is related to a large spin-orbit splitting of proton 2f states, is predicted only by forces which by far overestimate the proton spin-orbit splitting in Pb208. The Z=120 and N=172 shell closures predicted by the relativistic models and some Skyrme interactions are found to be related to a central depression of the nuclear density distribution. This effect cannot appear in macroscopic-microscopic models which have a limited freedom for the density distribution only. In summary, our findings give a strong argument for (Z=120,N=172) to be the next spherical doubly-magic superheavy nucleus.Comment: 22 pages REVTeX, 16 eps figures, accepted for publication in Phys. Rev.

Nuclear Ground State Observables and QCD Scaling in a Refined Relativistic Point Coupling Model

We present results obtained in the calculation of nuclear ground state properties in relativistic Hartree approximation using a Lagrangian whose QCD-scaled coupling constants are all natural (dimensionless and of order 1). Our model consists of four-, six-, and eight-fermion point couplings (contact interactions) together with derivative terms representing, respectively, two-, three-, and four-body forces and the finite ranges of the corresponding mesonic interactions. The coupling constants have been determined in a self-consistent procedure that solves the model equations for representative nuclei simultaneously in a generalized nonlinear least-squares adjustment algorithm. The extracted coupling constants allow us to predict ground state properties of a much larger set of even-even nuclei to good accuracy. The fact that the extracted coupling constants are all natural leads to the conclusion that QCD scaling and chiral symmetry apply to finite nuclei.Comment: 44 pages, 13 figures, 9 tables, REVTEX, accepted for publication in Phys. Rev.

Theoretical Aspects of Science with Radioactive Nuclear Beams

Physics of radioactive nuclear beams is one of the main frontiers of nuclear science today. Experimentally, thanks to technological developments, we are on the verge of invading the territory of extreme N/Z ratios in an unprecedented way. Theoretically, nuclear exotica represent a formidable challenge for the nuclear many-body theories and their power to predict nuclear properties in nuclear terra incognita. It is important to remember that the lesson learned by going to the limits of the nuclear binding is also important for normal nuclei from the neighborhood of the beta stability valley. And, of course, radioactive nuclei are crucial astrophysically; they pave the highway along which the nuclear material is transported up in the proton and neutron numbers during the complicated synthesis process in stars.Comment: 26 ReVTeX pages, 11 Postscript figures, uses epsf.sty, to be published in: Theme Issue on Science with Beams of Radioactive Nuclei, Philosophical Transactions, ed. by W. Gelletl

Axially symmetric Hartree-Fock-Bogoliubov Calculations for Nuclei Near the Drip-Lines

Nuclei far from stability are studied by solving the Hartree-Fock-Bogoliubov (HFB) equations, which describe the self-consistent mean field theory with pairing interaction. Calculations for even-even nuclei are carried out on two-dimensional axially symmetric lattice, in coordinate space. The quasiparticle continuum wavefunctions are considered for energies up to 60 MeV. Nuclei near the drip lines have a strong coupling between weakly bound states and the particle continuum. This method gives a proper description of the ground state properties of such nuclei. High accuracy is achieved by representing the operators and wavefunctions using the technique of basis-splines. The detailed representation of the HFB equations in cylindrical coordinates is discussed. Calculations of observables for nuclei near the neutron drip line are presented to demonstrate the reliability of the method.Comment: 13 pages, 4 figures. Submitted to Physical Review C on 05/08/02. Revised on Dec/0

Pairing Properties In Relativistic Mean Field Models Obtained From Effective Field Theory

We apply recently developed effective field theory nuclear models in mean field approximation (parameter sets G1 and G2) to describe ground-state properties of nuclei from the valley of $\beta$-stability up to the drip lines. For faster calculations of open-shell nuclei we employ a modified BCS approach which takes into account quasi-bound levels owing to their centrifugal barrier, with a constant pairing strength. We test this simple prescription by comparing with available Hartree-plus-Bogoliubov results. Using the new effective parameter sets we then compute separation energies, density distributions and spin--orbit potentials in isotopic (isotonic) chains of nuclei with magic neutron (proton) numbers. The new forces describe the experimental systematics similarly to conventional non-linear $\sigma-\omega$ relativistic force parameters like NL3.Comment: 29 pages, 17 figures, accepted for publication in PR

Shell Corrections of Superheavy Nuclei in Self-Consistent Calculations

Shell corrections to the nuclear binding energy as a measure of shell effects in superheavy nuclei are studied within the self-consistent Skyrme-Hartree-Fock and Relativistic Mean-Field theories. Due to the presence of low-lying proton continuum resulting in a free particle gas, special attention is paid to the treatment of single-particle level density. To cure the pathological behavior of shell correction around the particle threshold, the method based on the Green's function approach has been adopted. It is demonstrated that for the vast majority of Skyrme interactions commonly employed in nuclear structure calculations, the strongest shell stabilization appears for Z=124, and 126, and for N=184. On the other hand, in the relativistic approaches the strongest spherical shell effect appears systematically for Z=120 and N=172. This difference has probably its roots in the spin-orbit potential. We have also shown that, in contrast to shell corrections which are fairly independent on the force, macroscopic energies extracted from self-consistent calculations strongly depend on the actual force parametrisation used. That is, the A and Z dependence of mass surface when extrapolating to unknown superheavy nuclei is prone to significant theoretical uncertainties.Comment: 14 pages REVTeX, 8 eps figures, submitted to Phys. Rev.

Cohomological tautness for Riemannian foliations

In this paper we present some new results on the tautness of Riemannian foliations in their historical context. The first part of the paper gives a short history of the problem. For a closed manifold, the tautness of a Riemannian foliation can be characterized cohomologically. We extend this cohomological characterization to a class of foliations which includes the foliated strata of any singular Riemannian foliation of a closed manifold