Researcher profile

Hirotada Okawa

Hirotada Okawa contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
9works
0followers
7topics
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

9 published item(s)

preprint2022arXiv

A novel formulation for the evolution of relativistic rotating stars

We present a new formulation to construct numerically equilibrium configurations of rotating stars in general relativity. Having in mind the application to their quasi static evolutions, we adopt a Lagrangian formulation of our own devising, in which we solve force balance equations to seek for the positions of fluid elements assigned to the grid points, instead of the ordinary Eulerian formulation. Unlike previous works in the literature, we do not employ the first integral of the Euler equation, which is not obtained by an analytic integration in general. We assign a mass, specific angular momentum and entropy to each fluid element in contrast to the previous methods, in which the spatial distribution of the angular velocity or angular momentum is specified. Those distributions are determined after the positions of all fluid elements (or grid points) are derived in our formulation. We solve the large system of algebraic nonlinear equations that are obtained by discretizing the time-independent Euler and Einstein equations in the finite-elements method by using our new multi-dimensional root-finding scheme, named the W4 method. To demonstrate the capability of our new formulation, we construct some rotational configurations both barotropic and baroclinic. We also solve three evolutionary sequences that mimic the cooling, mass-loss, and mass-accretion as simple toy models.

preprint2022arXiv

Primordial black hole formation from massless scalar isocurvature

We numerically study the primordial black hole (PBH) formation by an isocurvature perturbation of a massless scalar field on super Hubble scales in the radiation-dominated universe. As a first step we perform simulations of spherically symmetric configurations. For the initial condition, we employ the spatial gradient expansion and provide the general form of the growing mode solutions valid up through the second order in this expansion. The initial scalar field profile is assumed to be Gaussian with a characteristic comoving wavenumber $k$; $\sim\exp(-k^2R^2)$, where $R$ is the radial coordinate. We find that a PBH is formed for a sufficiently large amplitude of the scalar field profile. Nevertheless, we find that the late time behavior of the gravitational collapse is dominated by the dynamics of the fluid but not by the scalar field, which is analogous to the PBH formation from an adiabatic perturbation in the radiation-dominated universe.

preprint2022arXiv

Principal-Axis Analysis of the Eddington Tensor for the Early Post-Bounce Phase of Rotational Core-Collapse Supernovae

Using full Boltzmann neutrino transport, we performed two-dimensional (2D) core-collapse supernova simulations in axisymmetry for two progenitor models with 11.2M and 15.0M both rotational and non-rotational. We employed the results obtained in the early post-bounce phase (t < 20 ms) to assess performance under rapid rotation of some closure relations commonly employed in the truncated moment method. We first made a comparison in 1D under spherical symmetry, though, of the Eddington factor p defined in the fluid rest frame (FR). We confirmed that the maximum entropy closure for the Fermionic distribution (MEFD) performs better than others near the proto-neutron star surface, where p < 1/3 occurs, but does not work well even in 1D when the phase space occupancy satisfies e < 0.5 together with p < 1/3, the condition known to be not represented by MEFD. For the 2D models with the rapid rotation, we employed the principal-axis analysis of the Eddington tensor. We paid particular attention to the direction of the longest principal axis. We observed in FR that it is aligned neither with the radial direction nor with the neutrino flux in 2D, particularly so in convective and/or rapidly rotating regions, the fact not accommodated in the moment method. We repeated the same analysis in the laboratory frame (LB) and found again that the direction of the longest principal axis is not well reproduced by MEFD because the interpolation between the optically thick and thin limits is not very accurate in this frame.

preprint2022arXiv

Singularity-Avoiding Multi-Dimensional Root-Finder

We proposed in this paper a new method, which we named the W4 method, to solve nonlinear equation systems. It may be regarded as an extension of the Newton-Raphson~(NR) method to be used when the method fails. Indeed our method can be applied not only to ordinary problems with non-singular Jacobian matrices but also to problems with singular Jacobians, which essentially all previous methods that employ the inversion of the Jacobian matrix have failed to solve. In this article, we demonstrate that (i) our new scheme can define a non-singular iteration map even for those problems by utilizing the singular value decomposition, (ii) a series of vectors in the new iteration map converges to the right solution under a certain condition, (iii) the standard two-dimensional problems in the literature that no single method proposed so far has been able to solve completely are all solved by our new method.

preprint2021arXiv

ECO-spotting: looking for extremely compact objects with bosonic fields

Black holes are thought to describe the geometry of massive, dark compact objects in the universe. To further support and quantify this long-held belief requires knowledge of possible, if exotic alternatives. Here, we wish to understand how compact can self-gravitating solutions be. We discuss theories with a well-posed initial value problem, consisting in either a single self-interacting scalar, vector or both. We focus on spherically symmetric solutions, investigating the influence of self-interacting potentials into the compactness of the solutions, in particular those that allow for flat-spacetime solutions. We are able to connect such stars to hairy black hole solutions, which emerge as a zero-mass black hole. We show that such stars can have light rings, but their compactness is never parametrically close to that of black holes. The challenge of finding black hole mimickers to investigate full numerical-relativity binary setups remains open.

preprint2020arXiv

Fast collective neutrino oscillations inside the neutrino sphere in core-collapse supernovae

Neutrinos are believed to have a key role in the explosion mechanism of core-collapse supernovae as they carry most of the energy released by the gravitational collapse of a massive star. If their flavor is converted fast inside the neutrino sphere, the supernova explosion may be influenced. This paper is reporting the results of the extended work of our previous paper. We perform a thorough survey of the ELN crossing in one of our self-consistent, realistic Boltzmann simulations in two spatial dimensions under axisymmetry for the existence of the crossings between $ν_e$ and $\barν_e$ angular distributions, or the electron lepton number (ELN) crossing. We report for the first time the positive detections deep inside the core of the massive star in the vicinity of neutrino sphere at $r$ $\approx$ 16 - 21 km. We find that low values of the electron fraction $Y_e$ produced by convective motions together with the appearance of light elements are critically important to give rise to the ELN crossing by enhancing the chemical potential difference between proton and neutron, and hence by mitigating the Fermi-degeneracy of $ν_e$. Since the region of positive detection are sustained and, in fact, expanding with time, it may have an impact on the explosion of core-collapse supernovae, observational neutrino astronomy and nucleosynthesis of heavy nuclei.

preprint2020arXiv

Gravitational Waves from Hierarchical Triple Systems with Kozai-Lidov Oscillation

We study gravitational waves from a hierarchical three-body system up to first-order postNewtonian approximation. Under certain conditions, the existence of a nearby third body can cause periodic exchange between eccentricity of an inner binary and relative inclination, known as Kozai-Lidov oscillations. We analyze features of the waveform from the inner binary system undergoing such oscillations. We find that variation caused due to the tertiary companion can be observed in the gravitational waveforms and energy spectra, which should be compared with those from isolated binaries and coplanar three-body system. The detections from future space interferometers will make possible the investigation of the gravitational wave spectrum in mHz range and may fetch signals by sources addressed.

preprint2011arXiv

Binary Neutron Star Mergers: Dependence on the Nuclear Equation of State

We perform a numerical-relativity simulation for the merger of binary neutron stars with 6 nuclear-theory-based equations of state (EOSs) described by piecewise polytropes. Our purpose is to explore the dependence of the dynamical behavior of the binary neutron star merger and resulting gravitational waveforms on the EOS of the supernuclear-density matter. The numerical results show that the merger process and the first outcome are classified into three types; (i) a black hole is promptly formed, (ii) a short-lived hypermassive neutron star (HMNS) is formed, (iii) a long-lived HMNS is formed. The type of the merger depends strongly on the EOS and on the total mass of the binaries. For the EOS with which the maximum mass is larger than 2Msun, the lifetime of the HMNS is longer than 10 ms for a total mass m_0=2.7Msun. A recent radio observation suggests that the maximum mass of spherical neutron stars is M_max \geq 1.97\pm 0.04Msun in one σlevel. This fact and our results support the possible existence of a HMNS soon after the onset of the merger for a typical binary neutron star with m_0=2.7Msun. We also show that the torus mass surrounding the remnant black hole is correlated with the type of the merger process; the torus mass could be large, \geq 0.1Msun, in the case that a long-lived HMNS is formed. We also show that gravitational waves carry information of the merger process, the remnant, and the torus mass surrounding a black hole.

preprint2011arXiv

Is super-Planckian physics visible? -- Scattering of black holes in 5 dimensions

It may be widely believed that probing short-distance physics is limited by the presence of the Planck energy scale above which scale any information is cloaked behind a horizon. If this hypothesis is correct, we could observe quantum behavior of gravity only through a black hole of Planck mass. We numerically show that in a scattering of two black holes in the 5-dimensional spacetime, a visible domain, whose curvature radius is much shorter than the Planck length, can be formed. Our result indicates that super-Planckian phenomena may be observed without an obstruction by horizon formation in particle accelerators.