Source author record

Takashi Nakatsukasa

Takashi Nakatsukasa 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

50works
9topics
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

50 published item(s)

preprint2021arXiv

Microscopic collective inertial masses for nuclear reaction in the presence of nucleonic effective mass

Collective inertial mass coefficients with respect to translational, relative, and rotational motions are microscopically calculated, along the collective reaction path self-consistently determined, based on the adiabatic self-consistent collective coordinate (ASCC) method. The impact of the time-odd component of the mean-field potential on the inertial masses are investigated. The results are compared with those calculated with the cranking formulae. The inertial masses based on the ASCC method reproduce the exact total nuclear mass for the translational motion as well as the exact reduced masses as the asymptotic values for the relative and rotational motions. In contrast, the cranking formulae fail to do so. This is due to the fact that the (local) Galilean invariance is properly restored in the ASCC method, but violated in the cranking formulae. A model Hamiltonian for low-energy nuclear reaction is constructed with the microscopically derived potentials and inertial masses. The astrophysical S-factors are calculated, which indicates the importance of microscopic calculation of proper inertial masses.

preprint2020arXiv

Collective inertial masses in nuclear reactions

Towards the microscopic theoretical description for large amplitude collective dynamics, we calculate the coefficients of inertial masses for low-energy nuclear reactions. Under the scheme of energy density functional, we apply the adiabatic self-consistent collective coordinate (ASCC) method, as well as the Inglis' cranking formula to calculate the inertias for the translational and the relative motions, in addition to those for the rotational motion. Taking the scattering between two $α$ particles as an example, we investigate the impact of the time-odd components of the mean-field potential on the collective inertial masses. The ASCC method asymptotically reproduces the exact masses for both the relative and translational motions. On the other hand, the cranking formula fails to do so when the time-odd components exist.

preprint2020arXiv

Coordinate-space solver for finite-temperature Hartree-Fock-Bogoliubov calculation using the shifted Krylov method

In order to study structure of proto-neutron stars and those in subsequent cooling stages, it is of great interest to calculate inhomogeneous hot and cold nuclear matter in a variety of phases. The finite-temperature Hartree-Fock-Bogoliubov (FT-HFB) theory is a primary choice for this purpose, however, its numerical calculation for superfluid (superconducting) many-fermion systems in three dimensions requires enormous computational costs. To study a variety of phases in the crust of hot and cold neutron stars, we propose an efficient method to perform the FT-HFB calculation with the three-dimensional (3D) coordinate-space representation. Recently, an efficient method based on the contour integral of Green's function with the shifted conjugate-orthogonal conjugate-gradient method has been proposed [Phys. Rev. C 95, 044302 (2017)]. We extend the method to the finite temperature, using the shifted conjugate-orthogonal conjugate-residual method. We benchmark the 3D coordinate-space solver of the FT-HFB calculation for hot isolated nuclei and fcc phase in the inner crust of neutron stars at finite temperature. The computational performance of the present method is demonstrated. Different critical temperatures of the quadrupole and the octupole deformations are confirmed for $^{146}$Ba. The robustness of the shape coexistence feature in $^{184}$Hg is examined. For the neutron-star crust, the deformed neutron-rich Se nuclei embedded in the sea of superfluid low-density neutrons appear in the fcc phase at the nucleon density of 0.045 fm$^{-3}$ and the temperature of $k_B T=200$ keV. The efficiency of the developed solver is demonstrated for nuclei and inhomogeneous nuclear matter at finite temperature. It may provide a standard tool for nuclear physics, especially for the structure of the hot and cold neutron-star matters.

preprint2020arXiv

Finite-amplitude method for collective inertia in spontaneous fission

Background: Microscopic description of spontaneous fission is one of the most challenging subjects in nuclear physics. It is necessary to evaluate the collective potential and the collective inertia along a fission path for a description of quantum tunneling in spontaneous or low-energy fission. In past studies of the fission dynamics based on nuclear energy density functional (EDF) theory, the collective inertia has been evaluated with the cranking approximation, which neglects dynamical residual effects. Purpose: The purpose is to provide a reliable and efficient method to include dynamical residual effects in the collective inertia for fission dynamics. Methods: We use the local quasiparticle random-phase approximation (LQRPA) to evaluate the collective inertia along a fission path obtained by the constrained Hartree-Fock-Bogoliubov method with the Skyrme EDF. The finite-amplitude method (FAM) with a contour integration technique enables us to efficiently compute the collective inertia in a large model space. Results: We evaluate the FAM-QRPA collective inertia along a symmetric fission path in $^{240}$Pu and $^{256}$Fm. The FAM-QRPA inertia is significantly larger than the one of the cranking approximation, and shows pronounced peaks around the ground state and the fission isomer. This is due to dynamical residual effects. Conclusions: To describe the spontaneous or low-energy fission, we provide a reliable and efficient method to construct the collective inertia with dynamical residual effects that have been neglected in most of EDF-based works in the past. We show the importance of dynamical residual effects to the collective inertia. This work will be a starting point for a systematic study of fission dynamics in heavy and superheavy nuclei to microscopically describe the nuclear large-amplitude collective motions.

preprint2016arXiv

Microscopic derivation of the Bohr-Mottelson collective Hamiltonian and its application to quadrupole shape dynamics

We discuss the nature of the low-frequency quadrupole vibrations from small-amplitude to large-amplitude regimes. We consider full five-dimensional quadrupole dynamics including three-dimensional rotations restoring the broken symmetries as well as axially symmetric and asymmetric shape fluctuations. Assuming that the time-evolution of the self-consistent mean field is determined by five pairs of collective coordinates and collective momenta, we microscopically derive the collective Hamiltonian of Bohr and Mottelson, which describes low-frequency quadrupole dynamics. We show that the five-dimensional collective Schrödinger equation is capable of describing large-amplitude quadrupole shape dynamics seen as shape coexistence/mixing phenomena. We summarize the modern concepts of microscopic theory of large-amplitude collective motion, which is underlying the microscopic derivation of the Bohr-Mottelson collective Hamiltonian.

preprint2016arXiv

Multiple period states of the superfluid Fermi gas in an optical lattice

We study multiple period states of a two-component unpolarized superfluid Fermi gas in an optical lattice along the Bardeen-Cooper-Schrieffer (BCS) to Bose-Einstein condensate (BEC) crossover. The existence of states whose period is a multiple of the lattice spacing is a direct consequence of the non-linear behavior of the gas, which is due to the presence of the order parameter associated with superfluidity. By solving Bogoliubov-de Gennes equations for a superfluid flow with finite quasimomentum, we find that, in the BCS side of the crossover, the multiple period states can be energetically favorable compared to the normal Bloch states and their survival time against dynamical instability drastically increases, suggesting that these states can be accessible in current experiments, in sharp contrast to the situation in BECs.

preprint2016arXiv

Self-consistent collective coordinate for reaction path and inertial mass

We propose a numerical method to determine the optimal collective reaction path for the nucleus-nucleus collision, based on the adiabatic self-consistent collective coordinate (ASCC) method. We use an iterative method combining the imaginary-time evolution and the finite amplitude method, for the solution of the ASCC coupled equations. It is applied to the simplest case, the $α-α$ scattering. We determine the collective path, the potential, and the inertial mass. The results are compared with other methods, such as the constrained Hartree-Fock method, the Inglis's cranking formula, and the adiabatic time-dependent Hartree-Fock (ATDHF) method.

preprint2016arXiv

Time-dependent density-functional description of nuclear dynamics

We present the basic concepts and recent developments in the time-dependent density functional theory (TDDFT) for describing nuclear dynamics at low energy. The symmetry breaking is inherent in nuclear energy density functionals (EDFs), which provides a practical description of important correlations at the ground state. Properties of elementary modes of excitation are strongly influenced by the symmetry breaking and can be studied with TDDFT. In particular, a number of recent developments in the linear response calculation have demonstrated their usefulness in description of collective modes of excitation in nuclei. Unrestricted real-time calculations have also become available in recent years, with new developments for quantitative description of nuclear collision phenomena. There are, however, limitations in the real-time approach; for instance, it cannot describe the many-body quantum tunneling. Thus, we treat the quantum fluctuations associated with slow collective motions assuming that time evolution of densities are determined by a few collective coordinates and momenta. The concept of collective submanifold is introduced in the phase space associated with the TDDFT and used to quantize the collective dynamics. Selected applications are presented to demonstrate the usefulness and quality of the new approaches. Finally, conceptual differences between nuclear and electronic TDDFT are discussed, with some recent applications to studies of electron dynamics in the linear response and under a strong laser field.

preprint2015arXiv

Microscopic derivation of the quadrupole collective Hamiltonian for shape coexistence/mixing dynamics

Assuming that the time-evolution of the self-consistent mean field is determined by five pairs of collective coordinate and collective momentum, we microscopically derive the collective Hamiltonian for low-frequency quadrupole modes of excitation. We show that the five-dimensional collective Schrödinger equation is capable of describing large-amplitude quadrupole shape dynamics seen as shape coexistence/mixing phenomena. We focus on basic ideas and recent advances of the approaches based on the time-dependent mean-field theory, but relations to other time-independent approaches are also briefly discussed.

preprint2015arXiv

Time-dependent density functional studies of nuclear quantum dynamics in large amplitudes

The time-dependent density functional theory (TDDFT) provides a unified description of the structure and reaction. The linear approximation leads to the random-phase approximation (RPA) which is capable of describing a variety of collective motion in a harmonic regime. Beyond the linear regime, we present applications of the TDDFT to nuclear fusion and fission reaction. In particular, the extraction of the internuclear potential and the inertial mass parameter is performed using two different methods. A fusion hindrance mechanism for heavy systems is investigated from the microscopic point of view. The canonical collective variables are determined by the adiabatic self-consistent collective coordinate method. Preliminary results of the spontaneous fission path, the potential, and the collective mass parameter are shown for 8Be --> alpha+alpha.

preprint2014arXiv

Finite amplitude method in linear response TDDFT calculations

The finite amplitude method is a feasible and efficient method for the linear response calculation based on the time-dependent density functional theory. It was originally proposed as a method to calculate the strength functions. Recently, new techniques have been developed for computation of normal modes (eigenmodes) of the quasiparticle-random-phase approximation. Recent advances associated with the finite amplitude method are reviewed.

preprint2014arXiv

Mean-Field Calculation Based on Proton-Neutron Mixed Energy Density Functionals

We have performed calculations based on the Skyrme energy density functional (EDF) that includes arbitrary mixing between protons and neutrons. In this framework, single-particle states are generalized as mixtures of proton and neutron components. The model assumes that the Skyrme EDF is invariant under the rotation in isospin space and the Coulomb force is the only source of the isospin symmetry breaking. To control the isospin of the system, we employ the isocranking method, which is analogous to the standard cranking approach used for describing high-spin states. Here, we present results of the isocranking calculations performed for the isobaric analog states in $A = 40$ and $A = 54$ nuclei.

preprint2014arXiv

Repulsive aspects of pairing correlation in nuclear fusion reaction

Numerical simulation on nuclear collisions are performed using the canonical-basis time-dependent Hartree-Fock-Bogoliubov theory (Cb-TDHFB) in the three-dimensional coordinate space. Comparing results of the Cb-TDHFB and the conventional time-dependent Hartree-Fock (TDHF) calculations, we study effects of the pairing correlation on fusion reaction of $^{22}$O+$^{22}$O, $^{52}$Ca+$^{52}$Ca, and $^{22}$O+$^{52}$Ca, using the Skyrme SkM$^*$ functional and a contact-type pairing energy functional. Although current results are yet preliminary, they may suggest that the pairing correlation could hinder the fusion probability at energies in the vicinity of the Coulomb barrier height. We also perform a calculation for heavier nuclei, $^{96}$Zn+$^{124}$Sn, which seems to suggest a similar hindrance effect.

preprint2014arXiv

Systematic investigation of low-lying dipole modes using the canonical-basis time-dependent Hartree-Fock-Bogoliubov theory

Systematic investigations of the electric dipole (E1) modes of excitation are performed using the canonical-basis time-dependent Hartree-Fock-Bogoliubov (Cb-TDHFB) theory. The Cb-TDHFB is able to describe dynamical pairing correlations in excited states of nuclear systems. We apply the method to the real-time calculation of linear response in even-even nuclei with Skyrme functionals. Effects of shell structure, neutron skin, deformation, and neutron chemical potential (separation energy) are studied in a systematic way. This reveals a number of characteristic features of the low-energy E1 modes. We also find a universal behavior in the low-energy E1 modes for heavy neutron-rich isotopes, which suggests the emergence of decoupled E1 peaks beyond N = 82.

preprint2013arXiv

Efficient calculation for the quasiparticle random-phase approximation matrix

We present an efficient numerical technique to evaluate the matrix of the (quasiparticle)-random-phase approximation, using the finite amplitude method (FAM). The method is tested in calculation of monopole excitations in 120Sn, compared with result obtained with the former iterative FAM. The neutron-pair-transfer modes are calculated with the present method and their character change in neutron-rich Pb isotopes is discussed. Computational aspects of different FAM approaches are also discussed for future applications to a large-scale computation.

preprint2013arXiv

Energy-density-functional calculations including the proton-neutron mixing

We present results of calculations based on the Skyrme energy density functional including the arbitrary mixing between protons and neutrons. In this framework, single-particle states are superpositions of proton and neutron components and the energy density functional is fully invariant with respect to three-dimensional rotations in the isospin space. The isospin of the system is controlled by means of the isocranking method, which carries over the standard cranking approach to the isospin space. We show numerical results of the isocranking calculations performed for isobaric analogue states in the A=14 and $A=40-56$ nuclei. We also present such results obtained for high-isospin states in $^{48}$Cr, with constraints on the isospin implemented by using the augmented Lagrange method.

preprint2013arXiv

Feasibility of the finite amplitude method in covariant density functional theory

Self-consistent relativistic random-phase approximation (RPA) in the radial coordinate representation is established by using the finite amplitude method (FAM). Taking the isoscalar giant monopole resonance in spherical nuclei as example, the feasibility of the FAM for the covariant density functionals is demonstrated, and the newly developed methods are verified by the conventional RPA calculations. In the present relativistic RPA calculations, the effects of the Dirac sea can be automatically taken into account in the coordinate-space representation. The rearrangement terms due to the density-dependent couplings can be implicitly calculated without extra computational costs in both iterative and matrix FAM schemes.

preprint2013arXiv

Finite-amplitude method: An extension to the covariant density functionals

The finite-amplitude method (FAM) is one of the most promising methods for optimizing the computational performance of the random-phase approximation (RPA) calculations in deformed nuclei. In this report, we will mainly focus on our recent progress in the self-consistent relativistic RPA established by using the FAM. It is found that the effects of Dirac sea can be taken into account implicitly in the coordinate-space representation and the rearrangement terms due to the density-dependent couplings can be treated without extra computational costs.

preprint2013arXiv

Low-energy E1 strength in select nuclei: Possible constraints on the neutron skins and the symmetry energy

Correlations between low-lying electric dipole (E1) strength and neutron skin thickness are systematically investigated with the fully self-consistent random-phase approximation using the Skyrme energy functionals. The presence of strong correlation among these quantities is currently under dispute. We find that the strong correlation is present in properly selected nuclei, namely in spherical neutron-rich nuclei in the region where the neutron Fermi levels are located at orbits with low orbital angular momenta. The significant correlation between the fraction of the energy-weighted sum value and the slope of the symmetry energy is also observed. The deformation in the ground state seems to weaken the correlation.

preprint2013arXiv

Pairing Effects in Nuclear Fusion Reaction

We simulate a heavy-ion collision using the canonical-basis time-dependent Hartree-Fock-Bogoliubov theory (Cb-TDHFB) treating pairing correlation in the three-dimensional coordinate space. We apply the Cb-TDHFB to 22O+22O collision with a contact-type pairing energy functional, and compare results of Cb-TDHFB and TDHF to investigate the effects of pairing correlations in nuclear fusion. Our results seem to indicate that pairing effects do not increase the fusion cross section in this system.

preprint2013arXiv

Shape evolution of giant resonances in Nd and Sm isotopes

Giant multipole resonances in Nd and Sm isotopes are studied by employing the quasiparticle-random-phase approximation on the basis of the Skyrme energy-density-functional method. Deformation effects on giant resonances are investigated in these isotopes which manifest a typical nuclear shape change from spherical to prolate shapes. The peak energy, the broadening, and the deformation splitting of the isoscalar giant monopole (ISGMR) and quadrupole (ISGQR) resonances agree well with measurements. The magnitude of the peak splitting and the fraction of the energy-weighted strength in the lower peak of the ISGMR reflect the nuclear deformation. The experimental data on ISGMR, ISGDR, and ISGQR are consistent with the nuclear-matter incompressibility $K \simeq 210-230$ MeV and the effective mass $m^*_0/m \simeq 0.8-0.9$. However, the high-energy octupole resonance (HEOR) in $^{144}$Sm seems to indicates a smaller effective mass, $m^*_0/m \simeq 0.7-0.8$. A further precise measurement of HEOR is desired to determine the effective mass.

preprint2013arXiv

Systematic investigation of E1 strength for the isotopes from Z = 28 to 50

We carry out a systematic study of electric dipole mode (E1) for neutron-rich isotopes from nickel (Z=28) to tin (Z=50) using a time-dependent mean-field theory. Our time-dependent scheme is the canonical-basis time-dependent Hartree-Fock-Bogoliubov theory which can self-consistently describe nuclear dynamics with pairing correlation. We focus our discussion on the pygmy dipole resonance (PDR) and E1 polarizability. The correlation between neutron-skin thickness and PDR strongly depends on the neutron number, but the correlation between the skin thickness and E1 polarizability is much more stable.

preprint2012arXiv

Density functional approaches to collective phenomena in nuclei: Time-dependent density-functional theory for perturbative and non-perturbative nuclear dynamics

We present the basic concepts and our recent developments in the density functional approaches with the Skyrme functionals for describing nuclear dynamics at low energy. The time-dependent density-functional theory (TDDFT) is utilized for the exact linear response with an external perturbation. For description of collective dynamics beyond the perturbative regime, we present a theory of a decoupled collective submanifold to describe for a slow motion based on the TDDFT. Selected applications are shown to demonstrate the quality of their performance and feasibility. Advantages and disadvantages in the numerical aspects are also discussed.

preprint2012arXiv

Linear-response calculation in the time-dependent density functional theory

Linear response calculations based on the time-dependent density-functional theory are presented. Especially, we report results of the finite amplitude method which we have recently proposed as an alternative and feasible approach to the (quasiparticle-)random-phase approximation. Calculated properties of the giant resonances and low-energy E1 modes are discussed. We found a universal linear correlation between the low-energy E1 strength and the neutron skin thickness.

preprint2012arXiv

Real-time calculations of many-body dynamics in quantum systems

Real-time computation of time-dependent quantum mechanical problems are presented for nuclear many-body problems. Quantum tunneling in nuclear fusion at low energy is described using a time-dependent wave packet. A real-time method of calculating strength functions using the time-dependent Schroedinger equation is utilized to properly treat the continuum boundary condition. To go beyond the few-body models,we resort to the density-functional theory. The nuclear mean-field models are briefly reviewed to illustrate its foundation and necessity of state dependence in effective interactions. This state dependence is successfully taken into account by the density dependence, leading to the energy density functional. Photoabsorption cross sections in 238U are calculated with the real-time method for the time-dependent density-functional theory.

preprint2012arXiv

Shape transition and fluctuation in neutron-rich Cr isotopes around N = 40

The spherical-to-prolate shape transition in neutron-rich Cr isotopes from N = 34 to 42 is studied by solving the collective Schr\"odinger equation for the five-dimensional quadrupole collective Hamiltonian. The collective potential and inertial functions are microscopically derived with use of the constrained Hartree-Fock-Bogoliubov plus local quasiparticle random-phase approximation method. Nature of the quadrupole collectivity of low-lying states is discussed by evaluating excitation spectra and electric quadrupole moments and transition strengths. The result of calculation indicates that Cr isotopes around 64Cr are prolately deformed but still possess transitional character; large-amplitude shape fluctuations dominate in their low-lying states.

preprint2012arXiv

Systematic study of low-lying E1 strength using the time-dependent mean field theory

We carry out systematic investigation of electric dipole (E1) mode from light to heavy nuclei, using a new time-dependent mean field theory: the Canonical-basis Time-Dependent Hartree-Fock-Bogoliubov (Cb-TDHFB) theory. The Cb-TDHFB in the three-dimensional coordinate space representation can deal with pairing correlation and any kind of deformation in the time-dependent framework. We report the neutron-number dependence of the low-energy E1 mode for light (A < 40) and heavy isotopes (A > 100) around N = 82.

preprint2011arXiv

Emergence of pygmy dipole resonances: Magic numbers and neutron skins

The pygmy dipole resonances (PDR) for even-even nuclei in 8=<Z=<40 are studied performing a systematic calculation of the random-phase approximation with the Skyrme functional of SkM*. The calculation is fully self-consistent and does not assume any symmetry in the nuclear shape of the ground state. In every isotopic chain, the PDR emerges by showing a peak of the E1 strength at energies less than 10 MeV. The E1 strength of the PDR strongly depends on the position of the Fermi level and shows a clear correlation with the occupation of the orbits with the orbital angular momenta less than 3\hbar (l =< 2). We also found a strong correlation between the isotopic dependence of the neutron skin thickness and the pygmy dipole strength. The fraction of the energy weighted strength exhausted by the PDR and the neutron skin thickness show a linear correlation with the universal rate of about 0.2/fm.

preprint2011arXiv

Finite amplitude method for the quasi-particle-random-phase approximation

We present the finite amplitude method (FAM) for superfluid systems. A Hartree-Fock-Bogoliubov code may be transformed into a code of the quasi-particle-random-phase approximation (QRPA) with simple modifications. This technique has advantages over the conventional QRPA calculations, such as coding feasibility and computational cost. We perform the fully self-consistent linear-response calculation for a spherical neutron-rich nucleus 174Sn, modifying the HFBRAD code, to demonstrate the accuracy, feasibility, and usefulness of the FAM.

preprint2011arXiv

Fragmentation of electric dipole strength in N=82 isotones

Fragmentation of the dipole strength in the N=82 isotones 140Ce, 142Nd and 144Sm is calculated using the second random-phase approximation (SRPA). In comparison with the result of the random-phase approximation (RPA), the SRPA provides the additional damping of the giant dipole resonance and the redistribution of the low-energy dipole strength. Properties of the low-energy dipole states are significantly changed by the coupling to two-particle-two-hole (2p2h) states, which are also sensitive to the correlation among the 2p2h states. Comparison with available experimental data shows a reasonable agreement for the low-energy E1 strength distribution.

preprint2011arXiv

Microscopic approach to large-amplitude deformation dynamics with local QRPA inertial masses

We have developed a new method for determining microscopically the fivedimensional quadrupole collective Hamiltonian, on the basis of the adiabatic self-consistent collective coordinate method. This method consists of the constrained Hartree-Fock-Bogoliubov (HFB) equation and the local QRPA (LQRPA) equations, which are an extension of the usual QRPA (quasiparticle random phase approximation) to non-HFB-equilibrium points, on top of the CHFB states. One of the advantages of our method is that the inertial functions calculated with this method contain the contributions of the time-odd components of the mean field, which are ignored in the widely-used cranking formula. We illustrate usefulness of our method by applying to oblate-prolate shape coexistence in 72Kr and shape phase transition in neutron-rich Cr isotopes around N=40.

preprint2011arXiv

Microscopic description of large-amplitude shape-mixing dynamics with inertial functions derived in local quasiparticle random-phase approximation

On the basis of the adiabatic self-consistent collective coordinate method, we develop an efficient microscopic method of deriving the five-dimensional quadrupole collective Hamiltonian and illustrate its usefulness by applying it to the oblate-prolate shape coexistence/mixing phenomena in proton-rich 68,70,72Se. In this method, the vibrational and rotational collective masses (inertial functions) are determined by local normal modes built on constrained Hartree-Fock-Bogoliubov states. Numerical calculations are carried out using the pairing-plus-quadrupole Hamiltonian including the quadrupole-pairing interaction. It is shown that the time-odd components of the moving mean-field significantly increase the vibrational and rotational collective masses in comparison with the Inglis-Belyaev cranking masses. Solving the collective Schroedinger equation, we evaluate excitation spectra, quadrupole transitions and moments. Results of the numerical calculation are in excellent agreement with recent experimental data and indicate that the low-lying states of these nuclei are characterized as an intermediate situation between the oblate-prolate shape coexistence and the so-called gamma unstable situation where large-amplitude triaxial-shape fluctuations play a dominant role.

preprint2011arXiv

Microscopic description of large-amplitude shape-mixing dynamics with local QRPA inertial functions

We introduce a microscopic approach to derive all the inertial functions in the five-dimensional quadrupole collective Hamiltonian. Local normal modes are evaluated on the constrained mean field in the quasiparticle random-phase approximation in order to derive the inertial functions. The collective Hamiltonians for neutron-rich Mg isotopes are determined with use of this approach, and the shape coexistence/mixing around the N = 20 region is analyzed.

preprint2011arXiv

Self-consistent description of nuclear photoabsorption cross sections

Several approaches to photonuclear reactions, based on the time-dependent density-functional theory, have been developed recently. The standard linearization leads to the random-phase approximation (RPA) or the quasiparticle-random-phase approximation (QRPA). We have developed a parallelized QRPA computer program for axially deformed nuclei. We also present a feasible approach to the (Q)RPA calculation, that is the finite amplitude method (FAM). We show results of photoabsorption cross sections for deformed nuclei using the QRPA and FAM calculations. Finally, the canonical-basis approach to the time-dependent Hartree-Fock-Bogoliubov method is presented, to demonstrate its feasibility and usefulness.

preprint2011arXiv

Shape fluctuations in the ground and excited 0+ states of 30Mg and 32Mg

Large-amplitude collective dynamics of shape phase transition in the low-lying states of 30-36Mg is investigated by solving the five-dimensional (5D) quadrupole collective Schroedinger equation. The collective masses and potentials of the 5D collective Hamiltonian are microscopically derived with use of the constrained Hartree-Fock-Bogoliubov plus local quasiparticle RPA method. Good agreement with the recent experimental data is obtained for the excited 0+ states as well as the ground bands. For 30Mg, the shape coexistence picture that the deformed excited 0+ state coexists with the spherical ground state approximately holds. On the other hand, large-amplitude quadrupole-shape fluctuations dominate in both the ground and the excited 0+ states in 32Mg, so that the interpretation of 'coexisting spherical excited 0+ state' based on the naive inversion picture of the spherical and deformed configurations does not hold.

preprint2011arXiv

Study of pygmy dipole resonance with a new time-dependent mean field theory

We carry out a systematic investigation on the low-energy electric dipole strength, which is often called pygmy dipole resonances (PDR), using the canonical-basis time-dependent Hartree-Fock-Bogoliubov (Cb-TDHFB) method. The Cb-TDHFB is a new method which is derived from TDHFB with an approximation analogous to the BCS theory that the pair potential is assumed to be diagonal in the time-dependent canonical basis. We apply the method to linear-response calculation for even-even nuclei. We report the neutron-number dependence of PDR in light (A < 70) and heavy isotopes (A > 100) around N = 82.

preprint2010arXiv

Canonical-basis time-dependent Hartree-Fock-Bogoliubov theory and linear-response calculations

We present simple equations for a canonical-basis formulation of the time-dependent Hartree-Fock-Bogoliubov (TDHFB) theory. The equations are obtained from the TDHFB theory with an approximation that the pair potential is assumed to be diagonal in the canonical basis. The canonical-basis formulation significantly reduces the computational cost. We apply the method to linear-response calculations for even-even light nuclei and demonstrate its capability and accuracy by comparing our results with recent calculations of the quasi-particle random-phase approximation with Skyrme functionals. We show systematic studies of E1 strength distributions for Ne and Mg isotopes. The evolution of the low-lying pygmy strength seems to be determined by the interplay of several factors, including the neutron excess, separation energy, neutron shell effects, deformation, and pairing.

preprint2010arXiv

Dipole responses in Nd and Sm isotopes with shape transitions

Photoabsorption cross sections of Nd and Sm isotopes from spherical to deformed even nuclei are systematically investigated by means of the quasiparticle-random-phase approximation based on the Hartree-Fock-Bogoliubov ground states (HFB+QRPA) using the Skyrme energy density functional. The gradual onset of deformation in the ground states as increasing the neutron number leads to characteristic features of the shape phase transition. The calculation well reproduce the isotopic dependence of broadening and emergence of a double-peak structure in the cross sections without any adjustable parameter. We also find that the deformation plays a significant role for low-energy dipole strengths. The $E1$ strengths are fragmented and considerably lowered in energy. The summed $E1$ strength up to 10 MeV is enhanced by a factor of five or more.

preprint2010arXiv

Linear response calculation using the canonical-basis TDHFB with a schematic pairing functional

A canonical-basis formulation of the time-dependent Hartree-Fock-Bogoliubov (TDHFB) theory is obtained with an approximation that the pair potential is assumed to be diagonal in the time-dependent canonical basis. The canonical-basis formulation significantly reduces the computational cost. We apply the method to linear-response calculations for even-even nuclei. E1 strength distributions for proton-rich Mg isotopes are systematically calculated. The calculation suggests strong Landau damping of giant dipole resonance for drip-line nuclei.

preprint2009arXiv

A Model Analysis of Triaxial Deformation Dynamics in Oblate-Prolate Shape Coexistence Phenomena

From a viewpoint of oblate-prolate symmetry and its breaking, we adopt the quadrupole collective Hamiltonian to study dynamics of triaxial deformation in shape coexistence phenomena. It accommodates the axially symmetric rotor model, the $γ$-unstable model, the rigid triaxial rotor model and an ideal situation for the oblate-prolate shape coexistence as particular cases. Numerical solutions of this model yield a number of interesting suggestions. (1) The relative energy of the excited 0+ state can be a signature of the potential shape along the $γ$ direction. (2) Specific E2 transition probabilities are sensitive to the breaking of the oblate-prolate symmetry. (3) Nuclear rotation may induce the localization of collective wave functions in the ($β, γ$) deformation space.

preprint2009arXiv

N-Body Nuclear Forces at Short Distances in Holographic QCD

We provide a calculation of N-body (N>2) nucleon interactions at short distances in holographic QCD. In the Sakai-Sugimoto model of large N_c massless QCD, N baryons are described by N Yang-Mills instantons in 5 spacetime dimensions. We compute a classical short distance interaction hamiltonian for N 'tHooft instantons. This corresponds to N baryons sharing identical classical spins and isospins. We find that genuine N-body nuclear forces turn out to vanish for N>2, at the leading order. This suggests that classical N-body forces are always suppressed compared with 2-body forces.

preprint2009arXiv

Self-consistent calculation of nuclear photoabsorption cross section: Finite amplitude method with Skyrme functionals in the three-dimensional real space

The finite amplitude method (FAM), which we have recently proposed (T. Nakatsukasa, T. Inakura, and K. Yabana, Phys. Rev. C 76, 024318 (2007)), simplifies significantly the fully self-consistent RPA calculation. Employing the FAM, we are conducting systematic, fully self-consistent response calculations for a wide mass region. This paper is intended to present a computational scheme to be used in the systematic investigation and to show the performance of the FAM for a realistic Skyrme energy functional. We implemented the method in the mixed representation in which the forward and backward RPA amplitudes are represented by indices of single-particle orbitals for occupied states and the spatial grid points for unoccupied states. We solve the linear response equation for a given frequency. The equation is a linear algebraic problem with a sparse non-hermitian matrix, which is solved with an iterative method. We show results of the dipole response for selected spherical and deformed nuclei. The peak energies of the giant dipole resonance agree well with measurements for heavy nuclei, while they are systematically underestimated for light nuclei. We also discuss the width of the giant dipole resonance in the fully self-consistent RPA calculation.

preprint2007arXiv

Finite amplitude method for the RPA solution

We propose a practical method to solve the random-phase approximation (RPA) in the self-consistent Hartree-Fock (HF) and density-functional theory. The method is based on numerical evaluation of the residual interactions utilizing finite amplitude of single-particle wave functions. The method only requires calculations of the single-particle Hamiltonian constructed with independent bra and ket states. Using the present method, the RPA calculation becomes possible with a little extension of a numerical code of the static HF calculation. We demonstrate usefulness and accuracy of the present method performing test calculations for isoscalar responses in deformed 20Ne.

preprint2007arXiv

Gauge-Invariant Formulation of Adiabatic Self-Consistent Collective Coordinate Method

The adiabatic self-consistent collective coordinate (ASCC) method is a practical microscopic theory of large-amplitude collective motions in nuclei with superfluidity. We show that its basic equations are invariant against transformations involving the gauge angle in the particle-number space. By virtue of this invariance, a clean separation between the large-amplitude collective motion and the pairing rotational motion can be achieved, enabling us to restore the particle-number symmetry broken by the Hartree-Fock-Bogoliubov (HFB) approximation. We formulate the ASCC method explicitly in a gauge-invariant form. In solving the ASCC equations, it is necessary to fix the gauge. Applying this new formulation to the multi-O(4) model, we compare different gauge-fixing procedures and demonstrate that calculations using different gauges indeed yield the same results for gauge-invariant quantities, such as the collective path and quantum spectra. We suggest a gauge-fixing prescription that seems most convenient in realistic calculations.

preprint2006arXiv

Suppressed fusion cross section for neutron halo nuclei

Fusion reactions of neutron-halo nuclei are investigated theoretically with a three-body model. The time-dependent wave-packet method is used to solve the three-body Schrodinger equation. The halo neutron behaves as a spectator during the Coulomb dissociation process of the projectile. The fusion cross sections of 11Be-209Bi and 6He-238U are calculated and are compared with measurements. Our calculation indicates that the fusion cross section is slightly hindered by the presence of weakly bound neutrons.

preprint2005arXiv

Effects of Time-Odd Components in Mean Field on Large Amplitude Collective Dynamics

We apply the adiabatic self-consistent collective coordinate (ASCC) method to the multi-O(4) model and study collective mass (inertia function) of the many-body tunneling motion. Comparing results with those of the exact diagonalization, we show that the ASCC method succeeds in describing gradual change of excitation spectra from an anharmonic vibration about the spherical shape to a doublet pattern associated with a deformed double-well potential possessing the oblate-prolate symmetry. The collective mass is significantly increased by the quadrupole-pairing contribution to time-odd components of the moving mean field. In contrast, the cranking (Inglis-Belyaev) mass based on the constrained mean field, which ignores the time-odd components, is smaller than the ASCC mass and fails to reproduce the exact spectra.

preprint2000arXiv

Photoabsorption spectra in the continuum of molecules and atomic clusters

We present linear response theories in the continuum capable of describing photoionization spectra and dynamic polarizabilities of finite systems with no spatial symmetry. Our formulations are based on the time-dependent local density approximation with uniform grid representation in the three-dimensional Cartesian coordinate. Effects of the continuum are taken into account either with a Green's function method or with a complex absorbing potential in a real-time method. The two methods are applied to a negatively charged cluster in the spherical jellium model and to some small molecules (silane, acetylene and ethylene).