Search

We present an overview of the Hartree-Fock-Bogoliubov (HFB) theory of nucleonic superfluidity for finite nuclei. After introducing basic concepts related to pairing correlations, we show how the correlated pairs are incorporated into the HFB wave function. Thereafter, we present derivation and structure of the HFB equations within the superfluid nuclear density functional formalism and discuss several aspects of the theory, including the unitarity of the Bogoliubov transformation in truncated single-particle and quasiparticle spaces, form of the pairing functional, structure of the HFB continuum, regularization and renormalization of pairing fields, and treatment of pairing in systems with odd particle numbers., Comment: Contributed chapter in "50 Years of Nuclear BCS", edited by R. A. Broglia and V. Zelevinsky, references corrected

Background: Recent experiments on beta-delayed fission in the mercury-lead region and the discovery of asym- metric fission in 180 Hg [1] have stimulated theoretical interest in the mechanism of fission in heavy nuclei. Purpose: We study fission modes and fusion valleys in 180 Hg and 198 Hg to reveal the role of shell effects in pre-scission region and explain the experimentally observed fragment mass asymmetry and its variation with A. Methods: We use the self-consistent nuclear density functional theory employing Skyrme and Gogny energy density functionals. Results: The potential energy surfaces in multi-dimensional space of collective coordinates, including elongation, triaxiality, reflection-asymmetry, and necking, are calculated for 180 Hg and 198 Hg. The asymmetric fission valleys - well separated from fusion valleys associated with nearly spherical fragments - are found in in both cases. The density distributions at scission configurations are studied and related to the experimentally observed mass splits. Conclusions: The energy density functionals SkM\ast and D1S give a very consistent description of the fission process in 180 Hg and 198 Hg. We predict a transition from asymmetric fission in 180 Hg towards more symmetric distribution of fission fragments in 198 Hg. For 180 Hg, both models yield 100 Ru/80 Kr as the most probable split. For 198 Hg, the most likely split is 108 Ru/90 Kr in HFB-D1S and 110 Ru/88 Kr in HFB-SkM\ast., Comment: 6 pages, 5 figures, to be published in Physical Review C

The demands of cutting-edge science are driving the need for larger and faster computing resources. With the rapidly growing scale of computing systems and the prospect of technologically disruptive architectures to meet these needs, scientists face the challenge of effectively using complex computational resources to advance scientific discovery. Multidisciplinary collaborating networks of researchers with diverse scientific backgrounds are needed to address these complex challenges. The UNEDF SciDAC collaboration of nuclear theorists, applied mathematicians, and computer scientists is developing a comprehensive description of nuclei and their reactions that delivers maximum predictive power with quantified uncertainties. This paper describes UNEDF and identifies attributes that classify it as a successful computational collaboration. We illustrate significant milestones accomplished by UNEDF through integrative solutions using the most reliable theoretical approaches, most advanced algorithms, and leadership-class computational resources., Comment: 12 pages, 15 figures, Proceedings for Conference on Computational Physics (CCP2011)

The self-consistent Hartree-Fock-Bogoliubov problem in large boxes can be solved accurately in the coordinate space with the recently developed solvers HFB-AX (2D) and MADNESS-HFB (3D). This is essential for the description of superfluid Fermi systems with complicated topologies and significant spatial extend, such as fissioning nuclei, weakly-bound nuclei, nuclear matter in the neutron star rust, and ultracold Fermi atoms in elongated traps. The HFB-AX solver based on B-spline techniques uses a hybrid MPI and OpenMP programming model for parallel computation for distributed parallel computation, within a node multi-threaded LAPACK and BLAS libraries are used to further enable parallel calculations of large eigensystems. The MADNESS-HFB solver uses a novel multi-resolution analysis based adaptive pseudo-spectral techniques to enable fully parallel 3D calculations of very large systems. In this work we present benchmark results for HFB-AX and MADNESS-HFB on ultracold trapped fermions., Comment: Conference on Computational Physics (CCP 2011) Proceeding

Background: An asymptotic normalization coefficient (ANC) characterizes the asymptotic form of a one-nucleon overlap integral required for description of nucleon-removal reactions. Purpose: We investigate the impact of the particle continuum on proton and neutron ANCs for mirror systems from $p$- and $sd$-shell regions. Method: We use the real-energy and complex-energy continuum shell model approaches. Results: We studied the general structure of the single-particle ANCs as a function of the binding energy and orbital angular momentum. We computed ANCs in mirror nuclei for different physical situations, including capture reactions to weakly-bound and unbound states. Conclusions: We demonstrated that the single-particle ANCs exhibit generic behavior that is different for charged and neutral particles. We verified the previously proposed relation [N.K. Timofeyuk, R.C. Johnson, and A.M. Mukhamedzhanov, Phys. Rev. Lett. 91, 232501 (2003); Phys. Rev. Lett. 97, 069904(E) (2006); N. K. Timofeyuk and P. Descouvemont, Phys. Rev. C 72, 064324 (2005)] between proton and neutron mirror ANCs. We find minor modifications if the spectroscopic strength is either localized in a single state or broadly distributed. For cases when several states couple strongly to the decay channel, these modifications may reach 30%., Comment: 18 pages, 17 figures, 9 tables

Clustering is one of the most complex phenomena known to the structure of atomic nuclei. A comprehensive description of this ubiquitous phenomenon goes beyond standard shell model and cluster model frameworks. We argue that clustering is a consequence of an openness of the nuclear many-body system. To illustrate this point, we study the near-threshold behavior of exceptional points in ^16 Ne and ^24 S., Comment: Proceedings "Frontier Issues in Physics of Exotic Nuclei" (YKIS2011/DCEN2011), 14 pages, 9 figures

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., Comment: Manuscript contains 5 pages, 2 figures, and 1 table. Submitted to Physical Review Letters

Stimulated by recent experimental discoveries, triaxial strongly deformed (TSD) states in $^{158}$Er at ultrahigh spins have been studied by means of the Skyrme-Hartree-Fock model and the tilted-axis-cranking method. Restricting the rotational axis to one of the principal axes -- as done in previous cranking calculations -- two well-defined TSD minima in the total Routhian surface are found for a given configuration: one with positive and another with negative triaxial deformation $\gamma$. By allowing the rotational axis to change direction, the higher-energy minimum is shown to be a saddle point. This resolves the long-standing question of the physical interpretation of the two triaxial minima at a very similar quadrupole shape obtained in the principal axis cranking approach. Several TSD configurations have been predicted, including a highly deformed band expected to cross lesser elongated TSD bands at the highest spins. Its transitional quadrupole moment $Q_t \approx 10.5$\,eb is close to the measured value of $\sim$11\,eb; hence, it is a candidate for the structure observed in experiment., Comment: 5 pages, 5 figures

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., Comment: 16 pages, 11 figures, accepted for publication in Phys. Rev. C

At the heart of many scientific applications is the solution of algebraic systems, such as linear systems of equations, eigenvalue problems, and optimization problems, to name a few. TOPS, which stands for Towards Optimal Petascale Simulations, is a SciDAC applied math center focused on the development of solvers for tackling these algebraic systems, as well as the deployment of such technologies in large-scale scientific applications of interest to the U.S. Department of Energy. In this paper, we highlight some of the solver technologies we have developed in optimization and matrix computations. We also describe some accomplishments achieved using these technologies in UNEDF, a SciDAC application project on nuclear physics., Comment: SciDAC 2011 Conference, July 10-14, 2011, Denver, CO; 5 pages, 2 tables, 2 figures

Within the complex-energy configuration interaction framework, we study correlations of valence neutrons to explain the behavior of charge radii in the neutron halo nuclei $^{6,8}$He. We find that the experimentally observed decrease of the charge radius between $^6$He and $^8$He is caused by a subtle interplay between three effects: dineutron correlations, a spin-orbit contribution to the charge radius, and a core swelling effect. We demonstrate that two-neutron angular correlations in the $2^+_1$ resonance of $^6$He differ markedly from the ground-state correlations in $^{6,8}$He., Comment: 5 pages, 5 figures

Reliable calculations of the structure of heavy elements are crucial to address fundamental science questions such as the origin of the elements in the universe. Applications relevant for energy production, medicine, or national security also rely on theoretical predictions of basic properties of atomic nuclei. Heavy elements are best described within the nuclear density functional theory (DFT) and its various extensions. While relatively mature, DFT has never been implemented in its full power, as it relies on a very large number (~ 10^9-10^12) of expensive calculations (~ day). The advent of leadership-class computers, as well as dedicated large-scale collaborative efforts such as the SciDAC 2 UNEDF project, have dramatically changed the field. This article gives an overview of the various computational challenges related to the nuclear DFT, as well as some of the recent achievements., Comment: Proceeding of the Invited Talk given at the SciDAC 2011 conference, Jul. 10-15, 2011, Denver, CO

The UNEDF SciDAC collaboration of nuclear theorists, applied mathematicians, and computer scientists is developing a comprehensive description of nuclei and their reactions that delivers maximum predictive power with quantified uncertainties. This paper illustrates significant milestones accomplished by UNEDF through integration of the theoretical approaches, advanced numerical algorithms, and leadership class computational resources., Comment: 5 pages, 8 figures

We present an efficient method for calculating strength functions using the finite amplitude method (FAM) for deformed superfluid heavy nuclei within the framework of the nuclear density functional theory. We demonstrate that FAM reproduces strength functions obtained with the fully self-consistent quasi-particle random-phase approximation (QRPA) at a fraction of computational cost. As a demonstration, we compute the isoscalar and isovector monopole strength for strongly deformed configurations in $^{240}$Pu by considering huge quasi-particle QRPA spaces. Our approach to FAM, based on Broyden's iterative procedure, opens the possibility for large-scale calculations of strength distributions in well-bound and weakly bound nuclei across the nuclear landscape., Comment: 5 pages, 3 figures

The quasi-particle energy spectrum of the Hartree-Fock-Bogoliubov (HFB) equations contains discrete bound states, resonances, and non-resonant continuum states. We study the structure of the unbound quasi-particle spectrum of weakly bound nuclei within several methods that do not rely on imposing scattering or outgoing boundary conditions. Various approximations are examined to estimate resonance widths. It is shown that the stabilization method works well for all HFB resonances except for very narrow ones. The Thomas-Fermi approximation to the non-resonant continuum has been shown to be very effective, especially for coordinate-space HFB calculations in large boxes that involve huge amounts of discretized quasi-particle continuum states., Comment: 12 pages,11 figures,submitted to PRC

The superallowed beta-decay rates that provide stringent constraints on physics beyond the Standard Model of particle physics are affected by nuclear structure effects through isospin-breaking corrections. The self-consistent isospin- and angular-momentum-projected nuclear density functional theory is used for the first time to compute those corrections for a number of Fermi transitions in nuclei from A=10 to A=74. The resulting leading element of the CKM matrix, |V_{ud}|= 0.97447(23), agrees well with the recent result by Towner and Hardy [Phys. Rev. C {\bf 77}, 025501 (2008)]., Comment: 4 pages, 2 figures; replaced with accepted version

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., Comment: 16 pages, submitted to Phys. Rev. C

We present the isospin- and angular-momentum-projected nuclear density functional theory (DFT) and its applications to the isospin-breaking corrections to the superallowed beta-decay rates in the vicinity of the N=Z line. A preliminary value obtained for the Cabbibo-Kobayashi-Maskawa matrix element, |V_{ud}|=0.97463(24), agrees well with the recent estimate by Towner and Hardy [Phys. Rev. C{\bf 77}, 025501 (2008)]. We also discuss new opportunities to study the symmetry energy by using the isospin-projected DFT., Comment: Presented at the 17th Nuclear Physics Conference "Maria & Pierre Curie" on "Symmetry and symmetry breaking in nuclear physics", Kazimierz Dolny, 22-25 September, 2010, Poland

Theoretical approaches that use one-body densities as dynamical variables, such as Hartree-Fock or the density functional theory (DFT), break isospin symmetry both explicitly, by virtue of charge-dependent interactions, and spontaneously. To restore the spontaneously broken isospin symmetry, we implemented the isospin-projection scheme on top of the Skyrme-DFT approach. This development allows for consistent treatment of isospin mixing in both ground and exited nuclear states. In this study, we apply this method to evaluate the isospin impurities in ground states of even-even and odd-odd N~Z nuclei. By including simultaneous isospin and angular-momentum projection, we compute the isospin-breaking corrections to the 0+ --> 0+ superallowed beta-decay., Comment: Presented at the International Conference on Nuclear Structure "Nuclear Landscape at the Limits", Zakopane, Poland, 2010

Collective mass tensor derived from the cranking approximation to the adiabatic time-dependent Hartree-Fock-Bogoliubov (ATDHFB) approach is compared with that obtained in the Gaussian Overlap Approximation (GOA) to the generator coordinate method. Illustrative calculations are carried out for one-dimensional quadrupole fission pathways in 256Fm. It is shown that the collective mass exhibits strong variations with the quadrupole collective coordinate. These variations are related to the changes in the intrinsic shell structure. The differences between collective inertia obtained in cranking and perturbative cranking approximations to ATDHFB, and within GOA, are discussed., Comment: 9 pages, RevTeX, 4 figures

The augmented Lagrangiam method (ALM), widely used in quantum chemistry constrained optimization problems, is applied in the context of the nuclear Density Functional Theory (DFT) in the self-consistent constrained Skyrme Hartree-Fock-Bogoliubov (CHFB) variant. The ALM allows precise calculations of multidimensional energy surfaces in the space of collective coordinates that are needed to, e.g., determine fission pathways and saddle points; it improves accuracy of computed derivatives with respect to collective variables that are used to determine collective inertia; and is well adapted to supercomputer applications., Comment: 6 pages, 3 figures; to appear in Eur. Phys. J. A

We carry out state-of-the-art optimization of a nuclear energy density of Skyrme type in the framework of the Hartree-Fock-Bogoliubov (HFB) theory. The particle-hole and particle-particle channels are optimized simultaneously, and the experimental data set includes both spherical and deformed nuclei. The new model-based, derivative-free optimization algorithm used in this work has been found to be significantly better than standard optimization methods in terms of reliability, speed, accuracy, and precision. The resulting parameter set UNEDFpre results in good agreement with experimental masses, radii, and deformations and seems to be free of finite-size instabilities. An estimate of the reliability of the obtained parameterization is given, based on standard statistical methods. We discuss new physics insights offered by the advanced covariance analysis., Comment: 19 pages, 11 figures

By applying the recently proposed antisymmetric superfluid local density approximation (ASLDA) to strongly interacting polarized atomic gases at unitarity in very elongated traps, we find families of Larkin-Ovchinnikov (LO) type of solutions with prominent transversal oscillation of pairing potential. These LO states coexist with a superfluid state having a smooth pairing potential. We suggest that the LO phase could be accessible experimentally by increasing adiabatically the trap aspect ratio. We show that the local asymmetry effects contained in ASLDA do not support a deformed superfluid core predicted by previous Bogoliubov-de Gennes treatments., Comment: 4 pages, 4 figures; revision; accepted for publication in Phys.Rev.A (Rapid Communication)

Naive dimensional analysis based on chiral effective theory, when adapted to nuclear energy density functionals, prescribes natural units and a hierarchy of contributions that could be used to constrain fits of generalized functionals. By applying these units, a large sample of Skyrme parametrizations is examined for naturalness, which is signaled by dimensionless coupling constants of order one. The bulk of the parameters are found to be natural, with an underlying scale consistent with other determinations. Significant deviations from unity are associated with deficiencies in the corresponding terms of particular functionals or with an incomplete optimization procedure., Comment: 5 pages, 2 figures, accepted for publication in Phys. Rev. C

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., Comment: 5 pages, 3 figures

Is there a connection between the branch point singularity at the particle emission threshold and the appearance of cluster states which reveal the structure of a corresponding reaction channel? Which nuclear states are most impacted by the coupling to the scattering continuum? What should be the most important steps in developing the theory that will truly unify nuclear structure and nuclear reactions? The common denominator of these questions is the continuum shell-model approach to bound and unbound nuclear states, nuclear decays, and reactions., Comment: Contribution to the "Focus issue on Open Problems in Nuclear Structure", Journal of Physics G

Spatial symmetries of the densities appearing in the nuclear Density Functional Theory are discussed. General forms of the local densities are derived by using methods of construction of isotropic tensor fields. The spherical and axial cases are considered., Comment: 12 LaTeX pages, submitted to International Journal of Modern Physics E

Isospin symmetry of atomic nuclei is explicitly broken by the charge-dependent interactions, primarily the Coulomb force. Within the nuclear density functional theory, isospin is also broken spontaneously. We propose a projection scheme rooted in a mean field theory, that allows the consistent treatment of isospin breaking in both ground and exited nuclear states. We demonstrate that this scheme is essentially free from spurious divergences plaguing particle-number and angular-momentum restoration approaches. Applications of the new technique include excited high-spin states in medium-mass N=Z nuclei, such as superdeformed bands and many-particle-many-hole terminating states., Comment: 13 pages, 7 figures, submitted to Phys. Rev. C

The isospin breaking effects due to the Coulomb interaction in weakly-bound nuclei are studied using the Gamow Shell Model, a complex-energy configuration interaction approach which simultaneously takes into account many-body correlations between valence nucleons and continuum effects. We investigate the near-threshold behavior of one-nucleon spectroscopic factors and the structure of wave functions along an isomultiplet. Illustrative calculations are carried out for the T=1 isobaric triplet. By using a shell-model Hamiltonian consisting of an isoscalar nuclear interaction and the Coulomb term, we demonstrate that for weakly bound or unbound systems the structure of isobaric analog states varies within the isotriplet and impacts the energy dependence of spectroscopic factors. We discuss the partial dynamical isospin symmetry present in isospin-stretched systems, in spite of the Coulomb interaction that gives rise to large mirror symmetry breaking effects., Comment: 10 pages, 5 figures ; published version

Symmetry properties of densities and mean fields appearing in the nuclear Density Functional Theory with pairing are studied. We consider energy functionals that depend only on local densities and their derivatives. The most important self-consistent symmetries are discussed: spherical, axial, space-inversion, and mirror symmetries. In each case, the consequences of breaking or conserving the time-reversal and/or proton-neutron symmetries are discussed and summarized in a tabulated form, useful in practical applications., Comment: 26 RevTex pages, 1 eps figure, 9 tables, submitted to Physical Review C

We study one-quasiproton excitations in the rare-earth region in the framework of the nuclear Density Functional Theory in the Skyrme-Hartree-Fock-Bogoliubov variant. The blocking prescription is implemented exactly, with the time-odd mean field fully taken into account. The equal filling approximation is compared with the exact blocking procedure. We show that both procedures are strictly equivalent when the time-odd channel is neglected, and discuss how nuclear alignment properties affect the time-odd fields. The impact of time-odd fields on calculated one-quasiproton bandhead energies is found to be rather small, of the order of 100-200 keV; hence, the equal filling approximation is sufficiently precise for most practical applications. The triaxial polarization of the core induced by the odd particle is studied. We also briefly discuss the occurrence of finite-size spin instabilities that are present in calculations for odd-mass nuclei when certain Skyrme functionals are employed., Comment: Accepted for publication in Phys. Rev. C after third report from the referee. One reference fixed

Our understanding of nuclear fission, a fundamental nuclear decay, is still incomplete due to the complexity of the process. In this paper, we describe a study of spontaneous fission using the symmetry-unrestricted nuclear density functional theory. Our results show that the observed bimodal fission can be explained in terms of pathways in multidimensional collective space corresponding to different geometries of fission products. We also predict a new phenomenon of trimodal spontaneous fission for some rutherfordium, seaborgium, and hassium isotopes., Comment: 6 pages, 6 figures; accepted for publication in Phys. Rev. C

A systematic study of fission-barrier dependence on excitation energy has been performed using the self-consistent finite-temperature Hartree-Fock+BCS (FT-HF+BCS) formalism with the SkM* Skyrme energy density functional. The calculations have been carried out for even-even superheavy nuclei with Z ranging between 110 and 124. For an accurate description of fission pathways, the effects of triaxial and reflection-asymmetric degrees of freedom have been fully incorporated. Our survey demonstrates that the dependence of isentropic fission barriers on excitation energy changes rapidly with particle number, pointing to the importance of shell effects even at large excitation energies characteristic of compound nuclei. The fastest decrease of fission barriers with excitation energy is predicted for deformed nuclei around N=164 and spherical nuclei around N=184 that are strongly stabilized by ground-state shell effects. For nuclei 240Pu and 256Fm, which exhibit asymmetric spontaneous fission, our calculations predict a transition to symmetric fission at high excitation energies due to the thermal quenching of static reflection asymmetric deformations., Comment: 5 pages, 4 figures

We present the self-consistent, non-perturbative analysis of isospin mixing using the nuclear density functional approach and the rediagonalization of the Coulomb interaction in the good-isospin basis. The largest isospin-breaking effects are predicted for N = Z nuclei and they quickly fall with the neutron excess. The unphysical isospin violation on the mean-field level, caused by the neutron excess, is eliminated by the proposed method. We find a significant dependence of the magnitude of isospin breaking on the parametrization of the nuclear interaction term. A rough correlation has been found between the isospin mixing parameter and the difference of proton and neutron rms radii. The theoretical framework described in this study is well suited to describe a variety of phenomena associated with isospin violation in nuclei, in particular the isospin symmetry-breaking corrections to superallowed Fermi beta decays., Comment: submitted to Physical Review Letters

The dependence of fission barriers on the excitation energy of the compound nucleus impacts the survival probability of superheavy nuclei synthesized in heavy-ion fusion reactions. In this work, we investigate the isentropic fission barriers by means of the self-consistent nuclear density functional theory. The relationship between isothermal and isentropic descriptions is demonstrated. Calculations have been carried out for $^{264}$Fm, $^{272}$Ds, $^{278}$112, $^{292}$114, and $^{312}$124. For nuclei around $^{278}$112 produced in "cold fusion" reactions, we predict a more rapid decrease of fission barriers with excitation energy as compared to the nuclei around $^{292}$114 synthesized in "hot fusion" experiments. This is explained in terms of the difference between the ground-state and saddle-point temperatures. The effect of the particle gas is found to be negligible in the range of temperatures studied., Comment: 4 pages, 5 figures(revised according to referee's comments)

Rotational motion of heated 72-Ge is studied within the microscopic Shell Model Monte Carlo approach. We investigate the the angular momentum alignment and nuclear pairing correlations associated with J-pi Cooper pairs as a function of the rotational frequency and temperature. The reentrance of pairing correlations with temperature is predicted at high rotational frequencies. It manifests itself through the anomalous behavior of specific heat and level density., Comment: 4 pages; 4 figures

Properties of strongly interacting, two-component finite Fermi systems are discussed within the recently developed coordinate-space Hartree-Fock-Bogoliubov (HFB) code {\hfbax}. Two illustrative examples are presented: (i) weakly bound deformed Mg isotopes, and (ii) spin-polarized atomic condensates in a strongly deformed harmonic trap., Comment: 4 pages, 2 figures, ENAM 2008 conference proceedings (EPJA)

We survey odd-even nuclear binding energy staggering using density functional theory with several treatments of the pairing interaction including the BCS, Hartree-Fock-Bogoliubov, and the Hartree-Fock-Bogoliubov with the Lipkin-Nogami approximation. We calculate the second difference of binding energies and compare with 443 measured neutron energy differences in isotope chains and 418 measured proton energy differences in isotone chains. The particle-hole part of the energy functional is taken as the SLy4 Skyrme parametrization and the pairing part of the functional is based on a contact interaction with possible density dependence. An important feature of the data, reproduced by the theory, is the sharp gap quenching at magic numbers. With the strength of the interaction as a free parameter, the theory can reproduce the data to an rms accuracy of about 0.25 MeV. This is slightly better than a single-parameter phenomenological description but slightly poorer than the usual two-parameter phenomenological form C/A^alpha . The following conclusions can be made about the performance of common parametrization of the pairing interaction: (i) there is a weak preference for a surface-peaked neutron-neutron pairing, which might be attributable to many-body effects; (ii) a larger strength is required in the proton pairing channel than in the neutron pairing channel; (iii) pairing strengths adjusted to the well-known spherical isotope chains are too weak to give a good overall fit to the mass differences., Comment: 13 pages, 9 figures

This work reviews foundations and applications of the complex-energy continuum shell model that provides a consistent many-body description of bound states, resonances, and scattering states. The model can be considered a quasi-stationary open quantum system extension of the standard configuration interaction approach for well-bound (closed) systems., Comment: Topical Review, J. Phys. G, Nucl. Part. Phys, in press (2008)

We have extended the density matrix renormalization group (DMRG) approach to two-fluid open many-fermion systems governed by complex-symmetric Hamiltonians. The applications are carried out for three- and four-nucleon (proton-neutron) systems within the Gamow Shell Model (GSM) in the complex-energy plane. We study necessary and sufficient conditions for the GSM+DMRG method to yield the correct ground state eigenvalue and discuss different truncation schemes within DMRG. The proposed approach will enable configuration interaction studies of weakly-bound and unbound strongly interacting complex systems which, because of a prohibitively large size of Fock space, cannot be treated by means of the direct diagonalization., Comment: 13 pages, 15 figures

The coordinate space formulation of the Hartree-Fock-Bogoliubov (HFB) method enables self-consistent treatment of mean-field and pairing in weakly bound systems whose properties are affected by the particle continuum space. Of particular interest are neutron-rich, deformed drip-line nuclei which can exhibit novel properties associated with neutron skin. To describe such systems theoretically, we developed an accurate 2D lattice Skyrme-HFB solver {\hfbax} based on B-splines. Compared to previous implementations, we made a number of improvements aimed at boosting the solver's performance. These include: explicit imposition of axiality and space inversion, use of the modified Broyden's method to solve self-consistent equations, and a partial parallelization of the code. {\hfbax} has been benchmarked against other HFB solvers, both spherical and deformed, and the accuracy of the B-spline expansion was tested by employing the multiresolution wavelet method. Illustrative calculations are carried out for stable and weakly bound nuclei at spherical and very deformed shapes, including constrained fission pathways. In addition to providing new physics insights, {\hfbax} can serve as a useful tool to assess the reliability and applicability of coordinate-space and configuration-space HFB solvers, both existing and in development., Comment: 12 pages,7 figs, submitted to Phys. Rev. C

In the framework of the Density Functional Theory for superconductors, we study the restoration of the particle number symmetry by means of the projection technique. Conceptual problems are outlined and numerical difficulties are discussed. Both are related to the fact that neither the many-body Hamiltonian nor the wave function of the system appear explicitly in the Density Functional Theory. Similar obstacles are encountered in self-consistent theories utilizing density-dependent effective interactions., Comment: 18 RevTex pages, 12 figures, submitted to Physical Review C

The additivity principle of the extreme shell model stipulates that an average value of a one-body operator be equal to the sum of the core contribution and effective contributions of valence (particle or hole) nucleons. For quadrupole moment and angular momentum operators, we test this principle for highly and superdeformed rotational bands in the A~130 nuclei. Calculations are done in the self-consistent cranked non-relativistic Hartree-Fock and relativistic Hartree mean-field approaches. Results indicate that the additivity principle is a valid concept that justifies the use of an extreme single-particle model in an unpaired regime typical of high angular momenta., Comment: 14 pages, 6 figures, submitted to Physics Review C

The presence of a particle continuum, both of a resonant and non-resonant character, can significantly impact spectroscopic properties of weakly bound nuclei and excited nuclear states close to, and above, the particle emission threshold. In the framework of the continuum shell model in the complex momentum-plane, the so-called Gamow Shell Model, we discuss salient effects of the continuum coupling on the one-neutron overlap integrals and the associated spectroscopic factors in neutron-rich helium and oxygen nuclei. In particular, we demonstrate a characteristic near-threshold energy dependence of the spectroscopic factors for different l-waves. We show also that the realistic radial overlap functions, which are needed for the description of transfer reactions, can be generated by single-particle wave functions of the appropriately chosen complex potential., Comment: 9 figures; 23 pages; corrected version; accepted in Nuclear Physics A

In the threshold region, the cross section and the associated overlap integral obey the Wigner threshold law that results in the Wigner-cusp phenomenon. Due to flux conservation, a cusp anomaly in one channel manifests itself in other open channels, even if their respective thresholds appear at a different energy. The shape of a threshold cusp depends on the orbital angular momentum of a scattered particle; hence, studies of Wigner anomalies in weakly bound nuclei with several low-lying thresholds can provide valuable spectroscopic information. In this work, we investigate the threshold behavior of spectroscopic factors in neutron-rich drip-line nuclei using the Gamow Shell Model, which takes into account many-body correlations and the continuum effects. The presence of threshold anomalies is demonstrated and the implications for spectroscopic factors are discussed., Comment: Accepted in Physical Review C Figure corrected

Theoretical predictions and experimental discoveries for neutron-rich, short-lived nuclei far from stability indicate that the familiar concept of nucleonic shell structure should be considered as less robust than previously thought. The notion of single-particle motion in exotic nuclei is reviewed with a particular focus on three aspects: (i) variations of nuclear mean field with neutron excess due to tensor interactions; (ii) importance of many-body correlations; and (iii) influence of open channels on properties of weakly bound and unbound nuclear states., Comment: 14 pages, 7 figures, submitted to Progress in Particle and Nuclear Physics, Proc. of the International School of Nuclear Physics 28th Course, Radioactive Beams, Nuclear Dynamics and Astrophysics, Erice-Sicily: 16 - 24 September 2006

We study nuclear high-spin states undergoing the transition to the fully stretched configuration with maximum angular momentum I_max within the space of valence nucleons. To this end, we perform a systematic theoretical analysis of non-fully-stretched I_max-2 and I_max-1 f_{7/2}^n seniority isomers and d_{3/2}^{-1} f_{7/2}^{n+1} intruder states in the A~44 nuclei from the lower-fp shell. We employ two theoretical approaches: (i) the density functional theory based on the cranked self-consistent Skyrme-Hartree-Fock method, and (ii) the nuclear shell model in the full sdfp configuration space allowing for 1p-1h cross-shell excitations. We emphasize the importance of restoration of broken angular momentum symmetry inherently obscuring the mean-field treatment of high-spin states. Overall good agreement with experimental data is obtained., Comment: submitted to Phys. Rev. C