Researcher profile

Takashi Nakatsukasa

Takashi Nakatsukasa contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
33works
0followers
4topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

33 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.

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

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

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.

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.

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

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 &#39;coexisting spherical excited 0+ state&#39; based on the naive inversion picture of the spherical and deformed configurations does not hold.

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 &#39;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.