BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Kazuhiro Yabana

Kazuhiro Yabana appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
nucl-thcond-mat.mtrl-sciphysics.opticsnucl-excond-mat.mes-hallphysics.atm-clusphysics.chem-phastro-ph.SRcond-mat.otherphysics.comp-phquant-ph

Catalog footprint

What is connected

35works
11topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

35 published item(s)

preprint2022arXiv

Electron spill-out effect on third-order optical nonlinearity of metallic nanostructure

Over the last three decades, plasmonics using metallic nanostructures has become central to nanophotonics research. Recently, its targets have been extended to nonlinear optical phenomena. In a nonlinear regime, quantum mechanical effects, as exemplified by electron spill-out on the surface of nanostructures, significantly influence the optical response. Quantum hydrodynamic theory (QHT) is a promising basis for the analysis of nonlinear optical responses involving such quantum mechanical effects. Herein, QHT was applied to calculate the third-order optical nonlinearity of a spherical metallic nanostructure. The results demonstrate how electron spill-out strongly affects plasmon resonance and third-order optical nonlinearity.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

First-principles method for nonlinear light propagation at oblique incidence

We have developed a computational method to describe the nonlinear light propagation of an intense and ultrashort pulse at oblique incidence on a flat surface. In the method, coupled equations of macroscopic light propagation and microscopic electron dynamics are simultaneously solved using a multiscale modeling. The microscopic electronic motion is described by first-principles time-dependent density functional theory. The macroscopic Maxwell equations that describe oblique light propagation are transformed into one-dimensional wave equations. As an illustration of the method, light propagation at oblique incidence on a silicon thin film is presented.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

High order harmonic generation in semiconductors driven at near- and mid-IR wavelengths

We study high order harmonics generation (HHG) in crystalline silicon and diamond subjected to near and mid-infrared laser pulses. We employ time-dependent density functional theory and solve the time-dependent Kohn-Sham equation in the single-cell geometry. We demonstrate that clear and clean HHG spectra can be generated with careful selection of the pulse duration. In addition, we simulate dephasing effects in a large silicon super-cell through displacement of atomic positions prepared by a molecular dynamics simulation. We compare our results with the previous calculations by Floss et al. [arXiv:1705.10707] [Phys. Rev. A 97, 011401(R) (2018)] on Diamond at 800 nm and by Tancogne-Dejean et al. [arXiv:1609.09298] [Phys. Rev. Lett. 118, 087403 (2017)] on Si at 3000 nm.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Infrared-Shielding of Plasmonic Random Metasurface Constructed by Cesium-Doped Tungsten Bronze

The heat-shielding properties of random metasurface, composed of spherical or spheroidal nanoparticles with random displacements and/or random deformation, were theoretically investigated using the finite difference time domain method. The effective coverage was defined using the total area of nanoparticles in the metasurface, and the robustness of the near-infrared light reflection against randomness was investigated. When the effective coverage was high, the near-infrared light reflection was reduced by at least 20% in both nanoparticle arrangement and shape randomness compared to the hexagonal close-packed perfect metasurface. In contrast, when effective coverage was low, the randomness of the nanoparticle arrangement had almost no effect on the near-infrared light reflection. Furthermore, the near-infrared light reflection performance was improved by the randomness of the nanoparticle shape.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Propagation effects in high-harmonic generation from dielectric thin films

Theoretical investigation is conducted of high-order harmonic generation (HHG) in silicon thin films to elucidate the effect of light propagation in reflected and transmitted waves. The first-principles simulations are performed of the process in which an intense pulsed light irradiates silicon thin films up to 3 $μ$m thickness. Our simulations are carried within the time-dependent density functional theory (TDDFT) with the account of coupled dynamics of the electromagnetic fields and the electronic motion. It was found that the intensity of transmission HHG gradually decreases with the thickness, while the reflection HHG becomes constant from a certain thickness. Detailed analyses show that transmission HHG have two origins: the HHG generated near the front edge and propagating to the back surface, and that generated near the back edge and emitted directly. The dominating mechanism of the transmission HHG is found to depend on the thickness of the thin film and the frequency of the HHG. At the film thickness of 1 $μ$m, the transmission HHG with the frequency below 20 eV is generated near the back edge, while that with the frequency above 20 eV is generated near the front edge and propagates from there to the back surface.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

First-principles study of ultrafast and nonlinear optical properties of graphite thin films

We theoretically investigate ultrafast and nonlinear optical properties of graphite thin films based on first-principles time-dependent density functional theory. We first calculate electron dynamics in a unit cell of graphite under a strong pulsed electric field and explore the transient optical properties of graphite. It is shown that the optical response of graphite shows a sudden change from conducting to insulating phase at a certain intensity range of the applied electric field. It also appears as a saturable absorption, the saturation in the energy transfer from the electric field to electrons. We next investigate a light propagation in graphite thin films by solving coupled dynamics of the electrons and the electromagnetic fields simultaneously. It is observed that the saturable absorption manifests in the propagation with small attenuation in the spatial region where the electric field amplitude is about $4 \sim 7 \times 10^{-2}$ V/Angstrom.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Nonlinear dynamics of electromagnetic field and valley polarization in WSe$_{2}$ monolayer

Linear and nonlinear optical response of WSe$_{2}$ monolayer is investigated by two-dimensional Maxwell plus time-dependent density functional theory with spin-orbit interaction. By applying the chiral resonant pulses, the electron dynamics along with high harmonic generation are examined at weak and strong laser fields. WSe$_{2}$ monolayer shows the linear optical response at the intensity I = 10$^{10}$~W/cm$^{2}$ while a complex nonlinear behavior is observed at I = 10$^{12}$~W/cm$^{2}$. The nonlinear response of WSe$_{2}$ monolayer in terms of saturable absorption is observed at strong laser field. By changing the chirality of the resonant light, a strong circular dichroic effect is observed in the excited state population. A relatively weak laser field shows effective valley polarization while strong field induces spin-polarized carrier peak between $K$($K'$) and $\mathitΓ$-point via nonlinear process. On the other hand, the strong laser field shows high harmonics up to the 11th order. Our results demonstrate that circularly polarized resonant pulse generate high harmonics in WSe$_{2}$ monolayer of order 3n$\pm1$.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Valley polarization control in WSe2 monolayer by a single-cycle laser pulse

Abstract The valley degree of freedom in two-dimensional materials provides an opportunity to extend the functionalities of valleytronics devices. Very short valley lifetimes demand the ultrafast control of valley pseudospin. Here, we theoretically demonstrate the control of valley pseudospin in WSe2 monolayer by single-cycle linearly polarized laser pulse. We use the asymmetric electric field controlled by the carrier-envelope phase (CEP) to make the valley polarization between K and K'-point in the Brillouin zone (BZ). Time-dependent density functional theory with spin-orbit interaction reveals that no valley asymmetry and its CEP dependence is observed within the linear-optical limit. In the nonlinear-optical regime, linearly polarized pulse induces a high degree of valley polarization and this polarization is robust against the field strength. Valley polarization strongly depends and oscillates as a function of CEP. The carrier density distribution forms nodes as the laser intensity increases, our results indicate that the position of the carrier density in the BZ can be controlled by the laser intensity. From the analysis by the massive Dirac Hamiltonian model, the nodes of the carrier density can be attributed to the Landau-Zener-Stückelberg interference of wave packets of the electron wave function.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Modulation of probe signal in coherent phonon detection revisited: Analytical and first-principles computational analyses

Modulation of probe signal in pump-probe measurements of coherent phonons in dielectrics, with and without spectral resolution, are investigated theoretically taking diamond as an example. Analytical investigation as well as first-principles calculations based on time-dependent density functional theory is utilized to clarify the mechanism of the modulation of the probe signals. Boundary and bulk effects are investigated systematically, putting emphasis on the phase relation between the modulation and the atomic motion of the coherent phonon. They are summarized as follows: Modulation by the boundary effect is in phase with the coherent phonon amplitude, while that by the bulk effect shows $π/2$ phase difference. Strong frequency dependence appears in the modulation by the bulk effect, while no frequency dependence by the boundary effect. First-principles calculations support the reliability of the analytical result.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Symmetry properties of attosecond transient absorption spectroscopy in crystalline dielectrics

We theoretically investigate a relation between the crystalline symmetry and the transient modulation of optical properties of crystalline dielectrics in pump-probe measurements using intense pump and attosecond probe fields. When the photon energy of the pump field is much below the bandgap energy, the modulation of the optical conductivity is caused by the intraband electronic motion, that is, the dynamical Franz-Keldysh effect. We analytically investigate symmetry properties of the modulated optical conductivity utilizing the Houston function, and derive a formula that relates the temporal oscillation in the absorption with the transformation properties of the modulated optical conductivity. To verify the validity of the formula, we perform real-time first-principles calculations based on the time-dependent density functional theory for a pump-probe process taking 4H-SiC crystal as an example.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Time-dependent Hartree-Fock calculations for multinucleon transfer and quasifission processes in the $^{64}$Ni+$^{238}$U reaction

Background: Multinucleon transfer (MNT) and quasifission (QF) processes are dominant processes in low-energy collisions of two heavy nuclei. They are expected to be useful to produce neutron-rich unstable nuclei. Nuclear dynamics leading to these processes depends sensitively on nuclear properties such as deformation and shell structure. Purpose: We elucidate reaction mechanisms of MNT and QF processes involving heavy deformed nuclei, making detailed comparisons between microscopic time-dependent Hartree-Fock (TDHF) calculations and measurements for the $^{64}$Ni+$^{238}$U reaction. Methods: Three-dimensional Skyrme-TDHF calculations are performed. Particle-number projection method is used to evaluate MNT cross sections from the TDHF wave function after collision. Results: Fragment masses, total kinetic energy (TKE), scattering angle, contact time, and MNT cross sections are investigated for the $^{64}$Ni+$^{238}$U reaction. They show reasonable agreements with measurements. At small impact parameters, collision dynamics depends sensitively on the orientation of deformed $^{238}$U. In tip (side) collisions, we find a larger (smaller) TKE and a shorter (longer) contact time. In tip collisions, we find a strong influence of quantum shells around $^{208}$Pb. Conclusions: It is confirmed that the TDHF calculations reasonably describe both MNT and QF processes in the $^{64}$Ni+$^{238}$U reaction. Analyses of this system indicates the significance of the nuclear structure effects such as deformation and quantum shells in nuclear reaction dynamics at low energies.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Controlling ultrafast currents by the non-linear photogalvanic effect

We theoretically investigate the effect of broken inversion symmetry on the generation and control of ultrafast currents in a transparent dielectric (SiO2) by strong femto-second optical laser pulses. Ab-initio simulations based on time-dependent density functional theory predict ultrafast DC currents that can be viewed as a non-linear photogalvanic effect. Most surprisingly, the direction of the current undergoes a sudden reversal above a critical threshold value of laser intensity I_c ~ 3.8*10^13 W/cm2. We trace this switching to the transition from non-linear polarization currents to the tunneling excitation regime. We demonstrate control of the ultrafast currents by the time delay between two laser pulses. We find the ultrafast current control by the non-linear photogalvanic effect to be remarkably robust and insensitive to laser-pulse shape and carrier-envelope phase.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Multinucleon transfer reaction in time-dependent Hartree-Fock theory

Time-dependent Hartree-Fock (TDHF) theory has achieved a remarkable success in describing and understanding nuclear many-body dynamics from nucleons' degrees of freedom. We here report our investigation of multinucleon transfer (MNT) processes employing the TDHF theory. To calculate transfer probabilities for channels specified by the number of protons and neutrons included in reaction products, a particle-number projection (PNP) method has been developed. The PNP method is also used to calculate excitation energies of reaction products. Combined use of the PNP method with a statistical model, we can evaluate MNT cross sections taking account of effects of particle evaporation. Using these methods, we evaluate MNT cross sections for $^{40,48}$Ca+$^{124}$Sn, $^{40}$Ca+$^{208}$Pb, and $^{58}$Ni+$^{208}$Pb reactions. From systematic analyses, we find that cross sections for channels with a large reaction probability are in good agreement with experimental data. However, the agreement becomes less accurate as the number of transferred nucleons increases. Possible directions to improve the description are discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Nonlinear electronic excitations in crystalline solids using meta-generalized gradient approximation and hybrid functional in time-dependent density functional theory

We develop numerical methods to calculate electron dynamics in crystalline solids in real-time time-dependent density functional theory employing exchange-correlation potentials which reproduce band gap energies of dielectrics; a meta generalized gradient approximation (meta-GGA) proposed by Tran and Blaha [Phys. Rev. Lett. 102, 226401 (2009)] (TBm-BJ) and a hybrid functional proposed by Heyd, Scuseria, and Ernzerhof [J. Chem. Phys. 118, 8207 (2003)] (HSE). In time evolution calculations employing the TB-mBJ potential, we have found it necessary to adopt a predictor-corrector step for stable time-evolution. Since energy functional is not known for the TB-mBJ potential, we propose a method to evaluate electronic excitation energy without referring to the energy functional. Calculations using the HSE hybrid functional is computationally expensive due to the nonlocal Fock-like term. We develop a computational method for the operation of the Fock-like term in Fourier space, for which we employ massively parallel computers equipped with graphic processing units. To demonstrate significances of utilizing potentials providing correct band gap energies, we compare electronic excitations induced by femtosecond laser pulses using the TB-mBJ, HSE, and a simple local density approximation (LDA). At low laser intensities, electronic excitations are found to be sensitive to the band gap energy: results using TB-mBJ and HSE are close to each other, while the excitation of the LDA calculation is more intensive than the others. At high laser intensities close to a damage threshold, we have found that electronic excitation energies are similar among the three cases.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Observation of molecular dipole excitations by attosecond self-streaking

We propose a protocol to probe the ultrafast evolution and dephasing of coherent electronic excitation in molecules in the time domain by the intrinsic streaking field generated by the molecule itself. Coherent electronic motion in the endohedral fullerene \Necsixty~is initiated by a moderately intense femtosecond UV-VIS pulse leading to coherent oscillations of the molecular dipole moment that persist after the end of the laser pulse. The resulting time-dependent molecular near-field is probed through the momentum modulation of photoemission from the central neon atom by a time-delayed attosecond XUV pulse. Our ab-initio time-dependent density functional theory and classical trajectory simulations predict that this self-streaking signal accurately traces the molecular dipole oscillations in real time. We discuss the underlying processes and give an analytical model that captures the essence of our ab-initio simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Efficient basis expansion for describing linear and nonlinear electron dynamics in crystalline solids

We propose an efficient basis expansion for electron orbitals to describe real-time electron dynamics in crystalline solids. Although a conventional grid representation in the three-dimensional Cartesian coordinates works robustly, it requires a large amount of computational resources. To reduce computational costs, we consider an expansion using basis functions with a truncation. A simple choice employing eigenstates of the ground state Hamiltonian with a truncation turned out to be useless. We have found that adding occupied eigenstates of nearby $k$-points to the truncated basis functions composed of eigenstates of the original $k$-point is crucially important. We demonstrate the usefulness of the method for linear and nonlinear electron dynamics calculations in crystalline SiO$_2$.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Particle-number projection method in time-dependent Hartree-Fock theory: Properties of reaction products

Background: The time-dependent Hartree-Fock (TDHF) theory has been successful in describing low-energy heavy ion collisions. Recently, we have shown that multinucleon transfer processes can be reasonably described in the TDHF theory combined with the particle-number projection technique. Purpose: In this work, we propose a theoretical framework to analyze properties of reaction products in TDHF calculations. Methods: TDHF calculation in three-dimensional Cartesian grid representation combined with particle number projection method. Results: We develop a theoretical framework to calculate expectation values of operators in the TDHF wave function after collision with the particle-number projection. To show how our method works in practice, the method is applied to $^{24}$O+$^{16}$O collisions for two quantities, angular momentum and excitation energy. The analyses revealed following features of the reaction: The nucleon removal proceeds gently, leaving small values of angular momentum and excitation energy in nucleon removed nuclei. Contrarily, nuclei receiving nucleons show expectation values of angular momentum and excitation energy which increase as the incident energy increases. Conclusions: We have developed a formalism to analyze properties of fragment nuclei in the TDHF theory combined with the particle-number projection technique. The method will be useful for microscopic investigations of reaction mechanisms in low-energy heavy ion collisions as well as for evaluating effects of particle evaporation on multinucleon transfer cross sections.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Strong Orientation Dependence of Multinucleon Transfer Processes in $^{238}$U+$^{124}$Sn Reaction

We theoretically investigate multinucleon transfer (MNT) processes in $^{238}$U+$^{124}$Sn reaction at $E_\mathrm{lab}=5.7$ MeV/$A$ using the time-dependent Hartree-Fock (TDHF) theory. For this reaction, measurements of MNT processes have been reported, showing substantial MNT cross sections accompanying more than ten protons. From the calculation, we find that the amount of transferred nucleons depends much on the relative orientation between the deformation axis of $^{238}$U and the relative vector connecting centers of $^{238}$U and $^{124}$Sn nuclei. We find a formation of thick neck when the $^{238}$U collides from its tip with $^{124}$Sn. However, the neck formation is substantially suppressed when $^{238}$U collides from its side. We have found that a large number of protons are transferred in the tip collision. This is caused by the breaking of the neck and subsequent absorption of nucleons in the neck region. We thus conclude that the measured MNT processes involving about ten protons originate from the neck breaking dynamics in the tip collisions of a deformed $^{238}$U nucleus.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Time-dependent Hartree-Fock calculations for multi-nucleon transfer processes: Effects of particle evaporation on production cross sections

We present a microscopic calculation of multi-nucleon transfer reactions employing the time-dependent Hartree-Fock (TDHF) theory. In our previous publication [Phys. Rev. C 88, 014614 (2013)], we reported our analysis for the multi-nucleon transfer processes for several systems. Here we discuss effects of particle evaporation processes on the production cross sections. Since particle evaporation processes may not be described adequately by the TDHF calculations, we evaluate them using a statistical model. As an input of the statistical model, excitation energies of the final fragments are necessary. We evaluate them from the TDHF wave function after collisions, extending the particle number projection technique. From the calculation, the particle evaporation effects are found to improve descriptions of the production cross sections. However, the production cross sections are still underestimated for processes where a number of protons are transferred. Possible origins of the discrepancy are discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

First-principles simulation of the optical response of bulk and thin-film α-quartz irradiated with an ultrashort intense laser pulse

A computational method based on a first-principles multiscale simulation has been used for calculating the optical response and the ablation threshold of an optical material irradiated with an ultrashort intense laser pulse. The method employs Maxwell's equations to describe laser pulse propagation and time-dependent density functional theory to describe the generation of conduction band electrons in an optical medium. Optical properties, such as reflectance and absorption, were investigated for laser intensities in the range $10^{10} \, \mathrm{W/cm^{2}}$ to $2 \times 10^{15} \, \mathrm{W/cm^{2}}$ based on the theory of generation and spatial distribution of the conduction band electrons. The method was applied to investigate the changes in the optical reflectance of $α$-quartz bulk, half-wavelength thin-film and quarter-wavelength thin-film and to estimate their ablation thresholds. Despite the adiabatic local density approximation used in calculating the exchange--correlation potential, the reflectance and the ablation threshold obtained from our method agree well with the previous theoretical and experimental results. The method can be applied to estimate the ablation thresholds for optical materials in general. The ablation threshold data can be used to design ultra-broadband high-damage-threshold coating structures.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Nuclear "pasta" structures in low-density nuclear matter and neutron star crust

In neutron star crust, non-uniform structure of nuclear matter is expected, which is called the "pasta" structure. From the recent studies of giant flares in magnetars, these structures might be related to some observables and physical quantities of the neutron star crust. To investigate the above quantities, we numerically explore the pasta structures with a fully threedimensional geometry and study the properties of low-density nuclear matter, based on the relativistic mean-field model and the Thomas-Fermi approximation. We observe typical pasta structures for fixed proton number-fraction and two of them for cold catalyzed matter. We also discuss the crystalline configuration of "pasta".

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Time-dependent Hartree-Fock calculations for multinucleon transfer processes in $^{40,48}$Ca+$^{124}$Sn, $^{40}$Ca+$^{208}$Pb, and $^{58}$Ni+$^{208}$Pb reactions

Multinucleon transfer processes in heavy-ion reactions at energies slightly above the Coulomb barrier are investigated in a fully microscopic framework of the time-dependent Hartree-Fock (TDHF) theory. Transfer probabilities are calculated from the TDHF wave function after collision using the projection operator method which has recently been proposed by Simenel (C. Simenel, Phys. Rev. Lett. 105, 192701 (2010)). We show results of the TDHF calculations for transfer cross sections of the reactions of $^{40,48}$Ca+$^{124}$Sn at $E_{lab}=$ 170, 174 MeV, $^{40}$Ca+$^{208}$Pb at $E_{lab}=$ 235, 249 MeV, and $^{58}$Ni+$^{208}$Pb at $E_{lab}=$ 328.4 MeV, for which measurements are available. We find the transfer processes show different behaviors depending on the $N/Z$ ratios of the projectile and the target, and the product of the charge numbers, $Z_P Z_T$. When the projectile and the target have different $N/Z$ ratios, fast transfer processes of a few nucleons towards the charge equilibrium of the initial system occur in reactions at large impact parameters. As the impact parameter decreases, a neck formation is responsible for the transfer. A number of nucleons are transferred by the neck breaking when two nuclei dissociate, leading to transfers of protons and neutrons in the same direction. Comparing cross sections by theory and measurements, we find the TDHF theory describes the transfer cross sections of a few nucleons reasonably. As the number of transferred nucleons increases, the agreement becomes less accurate. The TDHF calculation overestimates transfer cross sections accompanying a large number of neutrons when more than one proton are transferred. Comparing our results with those by other theories, we find the TDHF calculations give qualitatively similar results to those of direct reaction models such as GRAZING and Complex WKB.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Imaginary-time method for radiative capture reaction rate

We propose a new computational method for astrophysical reaction rate of radiative capture process. In the method, an evolution of a wave function is calculated along the imaginary-time axis which is identified as the inverse temperature. It enables direct evaluation of reaction rate as a function of temperature without solving any scattering problem. The method is tested for two-body radiative capture reaction, ${^{16}{\rm O}}(α,γ){^{20}{\rm Ne}}$, showing that it gives identical results to that calculated by the ordinary procedure. The new method will be suited for calculation of triple-alpha radiative capture rate for which an explicit construction of the scattering solution is difficult.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Three-dimensional calculation of inhomogeneous structure in low-density nuclear matter

In low-density nuclear matter which is relevant to the crust region of neutron stars and collapsing stage of supernovae, non-uniform structures called "nuclear pasta" are expected. So far, most works on nuclear pasta have used the Wigner-Seitz cell approximation with anzats about the geometrical structures like droplet, rod, slab and so on. We perform fully three-dimensional calculation of non-uniform nuclear matter for some cases with fixed proton ratios and in beta-equilibrium based on the relativistic mean-field model and the Thomas-Fermi approximation. In our calculation typical pasta structures are observed. However, there appears some difference in the density region of each pasta structure.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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).

0 citations0 reviewsNo code linkedOpen access