
The density dependent term in Skyrme forces is essential, which simulates
threebody and manybody correlations beyond the lowmomentum twobody
interaction. We speculate that a single density term may be insufficient and a
higherorder density dependent term is added. The present work investigates the
influences of higherorder 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 centerofmass
corrections and LipkinNogami 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 shapecoexisting nucleus $^{40}$Mg. The calculations were
done in a fully selfconsistent continuum finiteamplitude 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 longsought corehalo 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
HartreeFockBogoliubov equation, we extend the secondorder ThomasFermi
approximation to generalized superfluid Fermi systems by including the
densitydependent effective mass and the spinorbit 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 ThomasFermi method are intended for various
inhomogeneous superfluid Fermi systems.

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

Weaklybound deformed nuclei have been studied by the Skyrme
HartreeFockBogoliubov (HFB) approach in large coordinatespace boxes. In
particular, the boxsize dependence of the HFB calculations of weaklybound
deformed nuclei are investigated, including the particle density and pairing
density distributions at nuclear surfaces, the nearthreshold 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 largebox calculations of weaklybound nuclei are important to
precisely describe exotic phenomena such as deformed halos and peninsulas of
stability beyond drip lines.

The selfconsistent HartreeFockBogoliubov problem in large boxes can be
solved accurately in the coordinate space with the recently developed solvers
HFBAX (2D) and MADNESSHFB (3D). This is essential for the description of
superfluid Fermi systems with complicated topologies and significant spatial
extend, such as fissioning nuclei, weaklybound nuclei, nuclear matter in the
neutron star rust, and ultracold Fermi atoms in elongated traps. The HFBAX
solver based on Bspline techniques uses a hybrid MPI and OpenMP programming
model for parallel computation for distributed parallel computation, within a
node multithreaded LAPACK and BLAS libraries are used to further enable
parallel calculations of large eigensystems. The MADNESSHFB solver uses a
novel multiresolution analysis based adaptive pseudospectral techniques to
enable fully parallel 3D calculations of very large systems. In this work we
present benchmark results for HFBAX and MADNESSHFB 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
SkyrmeHartreeFock model and the tiltedaxiscranking method. Restricting the
rotational axis to one of the principal axes  as done in previous cranking
calculations  two welldefined 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 higherenergy minimum is shown to be a saddle point. This
resolves the longstanding 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 quasiparticle energy spectrum of the HartreeFockBogoliubov (HFB)
equations contains discrete bound states, resonances, and nonresonant
continuum states. We study the structure of the unbound quasiparticle 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 ThomasFermi approximation to the nonresonant continuum has been shown
to be very effective, especially for coordinatespace HFB calculations in large
boxes that involve huge amounts of discretized quasiparticle 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 LarkinOvchinnikov (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 Bogoliubovde Gennes treatments.

A systematic study of fissionbarrier dependence on excitation energy has
been performed using the selfconsistent finitetemperature HartreeFock+BCS
(FTHF+BCS) formalism with the SkM* Skyrme energy density functional. The
calculations have been carried out for eveneven superheavy nuclei with Z
ranging between 110 and 124. For an accurate description of fission pathways,
the effects of triaxial and reflectionasymmetric 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 groundstate
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
heavyion fusion reactions. In this work, we investigate the isentropic fission
barriers by means of the selfconsistent 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 groundstate and saddlepoint temperatures. The effect of the particle gas
is found to be negligible in the range of temperatures studied.

Properties of strongly interacting, twocomponent finite Fermi systems are
discussed within the recently developed coordinatespace
HartreeFockBogoliubov (HFB) code {\hfbax}. Two illustrative examples are
presented: (i) weakly bound deformed Mg isotopes, and (ii) spinpolarized
atomic condensates in a strongly deformed harmonic trap.

The coordinate space formulation of the HartreeFockBogoliubov (HFB) method
enables selfconsistent treatment of meanfield and pairing in weakly bound
systems whose properties are affected by the particle continuum space. Of
particular interest are neutronrich, deformed dripline nuclei which can
exhibit novel properties associated with neutron skin. To describe such systems
theoretically, we developed an accurate 2D lattice SkyrmeHFB solver {\hfbax}
based on Bsplines. 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 selfconsistent 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 Bspline
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 coordinatespace and
configurationspace 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 meanfieldtype 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 parameterfree 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 selfconsistent
meanfield model, with the folding factor determined using the quantization
conditions of the quasibound cluster state. We have investigated the
$\alpha$particle and heaviercluster decays of trans$^{100}$Sn and
trans$^{208}$Pb nuclei, and the observed diproton emission from the proton
dripline nucleus $^{16}$Ne, showing the overall reasonable descriptions of
cluster radioactivities with calculated halflives 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 SkyrmeHartreeFock 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
groundstate deformations, nuclear radii, neutron skin thicknesses and
$\alpha$decay energies were also calculated. Density distributions were
discussed with the calculations of singleparticle 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.