The density dependent term in Skyrme forces is essential, which simulates three-body and many-body correlations beyond the low-momentum two-body interaction. We speculate that a single density term may be insufficient and a higher-order density dependent term is added. The present work investigates the influences of higher-order density dependencies based on extended UNEDF0 and SkM* forces. The global descriptions of nuclear masses and charge radii have been presented. Consequently the extended UNEDF0 force gives a global rms error on binding energies of 1.29 MeV. The influences on fission barriers and equation of state have also been investigated. The perspectives to improve Skyrme forces have also been discussed, including global center-of-mass corrections and Lipkin-Nogami pairing corrections.

In order to explore the nature of collective modes in weakly bound nuclei, we have investigated deformation effects and surface flow patterns of isovector dipole modes in a shape-coexisting nucleus $^{40}$Mg. The calculations were done in a fully self-consistent continuum finite-amplitude Quasiparticle Random Phase Approximation (QRPA) in a large deformed spatial mesh. An unexpected result of pygmy and giant dipole modes having disproportionate deformation splittings in strength functions was obtained. Furthermore, the transition current densities demonstrate that the long-sought core-halo oscillation in pygmy resonances is collective and compressional, corresponding to the lowest excitation energy and the simplest quantum flow topology. Our calculations show that surface flow patterns become more complicated as excitation energies increase.

Using the $\hbar$-expansion of the Green's function of the Hartree-Fock-Bogoliubov equation, we extend the second-order Thomas-Fermi approximation to generalized superfluid Fermi systems by including the density-dependent effective mass and the spin-orbit potential. We first implement and examine the full correction terms over different energy intervals of the quasiparticle spectra in calculations of finite nuclei. Final applications of this generalized Thomas-Fermi method are intended for various inhomogeneous superfluid Fermi systems.

Based on the Hartree-Fock-Bogoliubov solutions in large deformed coordinate spaces, the finite amplitude method for quasiparticle random phase approximation (FAM-QRPA) has been implemented, providing a suitable approach to probe collective excitations of weakly-bound nuclei embedded in the continuum. The monopole excitation modes in Magnesium isotopes up to the neutron drip line have been studied with the FAM-QRPA framework on both the coordinate-space and harmonic oscillator basis methods. Enhanced soft monopole strengths and collectivity as a result of weak-binding effects have been unambiguously demonstrated.

Weakly-bound deformed nuclei have been studied by the Skyrme Hartree-Fock-Bogoliubov (HFB) approach in large coordinate-space boxes. In particular, the box-size dependence of the HFB calculations of weakly-bound deformed nuclei are investigated, including the particle density and pairing density distributions at nuclear surfaces, the near-threshold resonant and continuum quasiparticle spectra, and energetic properties. The box size may have larger influences in pairing properties than in other bulk properties. We demonstrate that large-box calculations of weakly-bound nuclei are important to precisely describe exotic phenomena such as deformed halos and peninsulas of stability beyond drip lines.

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.

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.

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.

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.

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.

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.

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.

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.

The $\alpha$ particle and $^6$He emissions from possible molecular states in beryllium and carbon isotopes have been studied using a mean-field-type cluster potential. Calculations can reproduce well the $\alpha$-decay widths of excited states in $^{8}$Be, $^{12}$C and $^{20}$Ne. For the nucleus $^{10}$Be, we discussed the $\alpha$-decay widths with different shapes or decay modes, in order to understand the very different decay widths of two excited states. The widths of $^{6}$He decay from $^{12}$Be and $\alpha$ decays from $^{13,14}$C are predicted, which could be useful for future experiments.

We present parameter-free tunneling calculations for various cluster radioactivities including the diproton decays of atomic nuclei. An uniform folded cluster potential has been suggested that is based on a self-consistent mean-field model, with the folding factor determined using the quantization conditions of the quasibound cluster state. We have investigated the $\alpha$-particle and heavier-cluster decays of trans-$^{100}$Sn and trans-$^{208}$Pb nuclei, and the observed diproton emission from the proton drip-line nucleus $^{16}$Ne, showing the overall reasonable descriptions of cluster radioactivities with calculated half-lives agreeing well with experimental data. We have also predicted the properties of yet unobserved cluster decays of the exotic nuclei $^{112,114}$Ba, $^{104}$Te and $^{38}$Ti.

We employed the Skyrme-Hartree-Fock model to investigate the density distributions and their dependence on nuclear shapes and isospins in the superheavy mass region. Different Skyrme forces were used for the calculations with a special comparison to the experimental data in $^{208}$Pb. The ground-state deformations, nuclear radii, neutron skin thicknesses and $\alpha$-decay energies were also calculated. Density distributions were discussed with the calculations of single-particle wavefunctions and shell fillings. Calculations show that deformations have considerable effects on the density distributions, with a detailed discussion on the $^{292}$120 nucleus. Earlier predictions of remarkably low central density are not supported when deformation is allowed for.