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

Hirotada Okawa

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

ResearcherUnclaimed source record
gr-qcastro-ph.HEhep-phhep-thastro-ph.COastro-ph.IMastro-ph.SRmath.NANumerical Analysis

Catalog footprint

What is connected

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

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

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

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

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

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

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Hyperbolic Self-Gravity Solver for Large Scale Hydrodynamical Simulations

A new computationally efficient method has been introduced to treat self-gravity in mesh based hydrodynamical simulations. It is applied simply by slightly modifying the Poisson equation into an inhomogeneous wave equation. This roughly corresponds to the weak field limit of the Einstein equations in general relativity, and as long as the gravitation propagation speed is taken to be larger than the hydrodynamical characteristic speed, the results agree with solutions for the Poisson equation. The solutions almost perfectly agree if the domain is taken large enough, or appropriate boundary conditions are given. Our new method can not only significantly reduce the computational time compared with existent methods, but is also fully compatible with massive parallel computation, nested grids and adaptive mesh refinement techniques, all of which can accelerate the progress in computational astrophysics and cosmology.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Interaction between bosonic dark matter and stars

We provide a detailed analysis of how bosonic dark matter "condensates" interact with compact stars, extending significantly the results of a recent Letter. We focus on bosonic fields with mass $m_B$, such as axions, axion-like candidates and hidden photons. Self-gravitating bosonic fields generically form "breathing" configurations, where both the spacetime geometry and the field oscillate, and can interact and cluster at the center of stars. We construct stellar configurations formed by a perfect fluid and a bosonic condensate, and which may describe the late stages of dark-matter accretion onto stars, in dark matter-rich environments. These composite stars oscillate at a frequency which is a multiple of $f=2.5\times 10^{14}\,\left(m_{B}c^2/eV\right)\,{\rm Hz}$. Using perturbative analysis and Numerical Relativity techniques, we show that these stars are generically stable, and we provide criteria for instability. Our results also indicate that the growth of the dark matter core is halted close to the Chandrasekhar limit. We thus dispel a myth concerning dark matter accretion by stars: dark mater accretion does not necessarily lead to the destruction of the star, nor to collapse to a black hole. Finally, we argue that stars with long-lived bosonic cores may also develop in other theories with effective mass couplings, such as (massless) scalar-tensor theories.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Accretion of dark matter by stars

Searches for dark matter imprints are one of the most active areas of current research. We focus here on light fields with mass $m_B$, such as axions and axion-like candidates. Using perturbative techniques and full-blown nonlinear Numerical Relativity methods, we show that (i) dark matter can pile up in the center of stars, leading to configurations and geometries oscillating with frequency which is a multiple of f=$2.5 10^{14}$ $m_B c^2$/eV Hz. These configurations are stable throughout most of the parameter space, and arise out of credible mechanisms for dark-matter capture. Stars with bosonic cores may also develop in other theories with effective mass couplings, such as (massless) scalar-tensor theories. We also show that (ii) collapse of the host star to a black hole is avoided by efficient gravitational cooling mechanisms.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Black hole-neutron star binary merger: Dependence on black hole spin orientation and equations of state

We systematically performed numerical-relativity simulations for black hole (BH) - neutron star (NS) binary mergers with a variety of the BH spin orientation and equations of state (EOS) of the NS. The initial misalignment angles of the BH spin are chosen in the range of i_tilt,0 = 30--90[deg.]. We employed four models of NS EOS with which the compactness of the NS is in the range of C = M_NS/R_NS = 0.138--0.180, where M_NS and R_NS are the mass and the radius of the NS, respectively. The mass ratio of the BH to the NS, Q = M_BH/M_NS, and the dimensionless spin parameter of the BH, chi, are chosen to be Q = 5 and chi = 0.75, together with M_NS = 1.35 M_sun. We obtain the following results: (i) The inclination angle of i_tilt,0 < 70[deg.] and i_tilt,0 < 50[deg.] are required for the formation of a remnant disk with its mass larger than 0.1 M_sun for the case C = 0.140 and C = 0.160, respectively, while the disk mass is always smaller than 0.1M_sun for C = 0.175. The ejecta with its mass larger than 0.01 M_sun is obtained for i_tilt,0 < 85[deg.] with C = 0.140, for i_tilt,0 < 65[deg.] with C = 0.160, and for i_tilt,0 < 30[deg.] with C = 0.175. (ii) The rotational axis of the dense part of the remnant disk is approximately aligned with the remnant BH spin for i_tilt,0 = 30[deg.]. On the other hand, the disk axis is misaligned initially with ~ 30[deg.] for i_tilt,0 = 60[deg.], and the alignment with the remnant BH spin is achieved at ~ 50--60 ms after the onset of merger. The accretion time scale of the remnant disk is typically ~ 100 ms and depends only weakly on the misalignment angle and the EOS. (iii) The ejecta velocity is typically ~ 0.2--0.3c and depends only weakly on i_tilt,0 and the EOS of the NS, while the morphology of the ejecta depends on its mass. (iv) The gravitational-wave spectra contains the information of the NS compactness in the cutoff frequency for i_tilt,0 < 60[deg.].

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Collapse of massive fields in anti-de Sitter spacetime

Gravitational collapse in asymptotically anti-de Sitter spacetime has a rich but poorly-understood structure. There are strong indications that some families of initial data form "bound" states, which are regular everywhere, while other families seem to always collapse to black holes. Here, we investigate the collapse of massive scalar fields in anti-de Sitter, with enlarged freedom in the initial data setup, such as several distinct wavepackets, gravitationally interacting with each other. Our results are fully consistent with previous findings in the literature: massive fields, which have a fully resonant spectra, collapse at (arbitrarily?) small amplitude for some classes of initial data, and form oscillating stars for others. We find evidence that initial data consisting on several wavepackets may allow efficient exchange of energy between them, and delay the collapse substantially, or avoid it altogether. When the AdS boundary is artificially changed so that the spectrum is no longer resonant, cascading to higher frequencies may still be present. Finally, we comment on the asymptotically flat counterparts.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Dynamical mass ejection from black hole-neutron star binaries

We investigate properties of material ejected dynamically in the merger of black hole-neutron star binaries by numerical-relativity simulations. We systematically study the dependence of ejecta properties on the mass ratio of the binary, spin of the black hole, and equation of state of the neutron-star matter. Dynamical mass ejection is driven primarily by tidal torque, and the ejecta is much more anisotropic than that from binary neutron star mergers. In particular, the dynamical ejecta is concentrated around the orbital plane with a half opening angle of 10--20deg and often sweeps out only a half of the plane. The ejecta mass can be as large as ~0.1M_sun, and the velocity is subrelativistic with ~0.2--0.3c for typical cases. The ratio of the ejecta mass to the bound mass (disk and fallback components) is larger, and the ejecta velocity is larger, for larger values of the binary mass ratio, i.e., for larger values of the black-hole mass. The remnant black hole-disk system receives a kick velocity of O(100)km/s due to the ejecta linear momentum, and this easily dominates the kick velocity due to gravitational radiation. Structures of postmerger material, velocity distribution of the dynamical ejecta, fallback rates, and gravitational waves are also investigated. We also discuss the effect of ejecta anisotropy on electromagnetic counterparts, specifically a macronova/kilonova and synchrotron radio emission, developing analytic models.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Exploring tidal effects of coalescing binary neutron stars in numerical relativity II: Longterm simulations

We perform new longterm (15-16 orbits) simulations of coalescing binary neutron stars in numerical relativity using an updated Einstein's equation solver, employing low-eccentricity initial data, and modeling the neutron stars by a piecewise polytropic equation of state. A convergence study shows that our new results converge more rapidly than the third order and using the determined convergence order, we construct an extrapolated waveform for which the estimated total phase error should be less than 1 radian. We then compare the extrapolated waveforms with those calculated by the latest effective-one-body (EOB) formalism in which the so-called tidal deformability, higher post-Newtonian corrections, and gravitational self-force effects are taken into account. We show that for a binary of compact neutron stars with their radius 11.1 km, the waveform by the EOB formalism agrees quite well with the numerical waveform so that the total phase error is smaller than 1 radian for the total phase of $\sim 200$ radian up to the merger. By contrast, for a binary of less compact neutron stars with their radius 13.6 km, the EOB and numerical waveforms disagree with each other in the last few wave cycles, resulting in the total phase error of $\sim 3$ radian.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Testing General Relativity with Present and Future Astrophysical Observations

One century after its formulation, Einstein's general relativity has made remarkable predictions and turned out to be compatible with all experimental tests. Most of these tests probe the theory in the weak-field regime, and there are theoretical and experimental reasons to believe that general relativity should be modified when gravitational fields are strong and spacetime curvature is large. The best astrophysical laboratories to probe strong-field gravity are black holes and neutron stars, whether isolated or in binary systems. We review the motivations to consider extensions of general relativity. We present a (necessarily incomplete) catalog of modified theories of gravity for which strong-field predictions have been computed and contrasted to Einstein's theory, and we summarize our current understanding of the structure and dynamics of compact objects in these theories. We discuss current bounds on modified gravity from binary pulsar and cosmological observations, and we highlight the potential of future gravitational wave measurements to inform us on the behavior of gravity in the strong-field regime.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Black Hole Universe with $Λ$

Time evolution of a black hole lattice universe with a positive cosmological constant $Λ$ is simulated. The vacuum Einstein equations are numerically solved in a cubic box with a black hole in the center. Periodic boundary conditions on all pairs of opposite faces are imposed. Configurations of marginally trapped surfaces are analyzed. We describe the time evolution of not only black hole horizons, but also cosmological horizons. Defining the effective scale factor by using the area of a surface of the cubic box, we compare it with that in the spatially flat dust dominated FLRW universe with the same value of $Λ$. It is found that the behaviour of the effective scale factor is well approximated by that in the FLRW universe. Our result suggests that local inhomogeneities do not significantly affect the global expansion law of the universe irrespective of the value of $Λ$.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Black holes and fundamental fields in Numerical Relativity: initial data construction and evolution of bound states

Fundamental fields are a natural outcome in cosmology and particle physics and might therefore serve as a proxy for more complex interactions. The equivalence principle implies that all forms of matter gravitate, and one therefore expects relevant, universal imprints of new physics in strong field gravity, such as that encountered close to black holes. Fundamental fields in the vicinities of supermassive black holes give rise to extremely long-lived, or even unstable, configurations which slowly extract angular momentum from the black hole or simply evolve non-linearly over long timescales, with important implications for particle physics and gravitational-wave physics. Here, we perform a fully non-linear study of scalar-field condensates around rotating black holes. We provide novel ways to specify initial data for the Einstein-Klein-Gordon system, with potential applications in a variety of scenarios. Our numerical results confirm the existence of long-lived bar-modes which act as lighthouses for gravitational wave emission: the scalar field condenses outside the black hole geometry and acts as a constant frequency gravitational-wave source for very long timescales. This effect could turn out to be a potential signature of beyond standard model physics and also a promising source of gravitational waves for future gravitational wave detectors.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Black holes and fundamental fields: hair, kicks and a gravitational "Magnus" effect

Scalar fields pervade theoretical physics and are a fundamental ingredient to solve the dark matter problem, to realize the Peccei-Quinn mechanism in QCD or the string-axiverse scenario. They are also a useful proxy for more complex matter interactions, such as accretion disks or matter in extreme conditions. Here, we study the collision between scalar "clouds" and rotating black holes. For the first time we are able to compare analytic estimates and strong field, nonlinear numerical calculations for this problem. As the black hole pierces through the cloud it accretes according to the Bondi-Hoyle prediction, but is deflected through a purely kinematic gravitational "anti-Magnus" effect, which we predict to be present also during the interaction of black holes with accretion disks. After the interaction is over, we find large recoil velocities in the transverse direction. The end-state of the process belongs to the vacuum Kerr family if the scalar is massless, but can be a hairy black hole when the fundamental scalar is massive.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Collapse of self-interacting fields in asymptotically flat spacetimes: do self-interactions render Minkowski spacetime unstable?

The nonlinear instability of anti-de Sitter spacetime has recently been established with the striking result that generic initial data collapses to form black holes. This outcome suggests that confined matter generically collapses, and that collapse can only be halted -- at most -- by nonlinear bound states. Here we provide evidence that such mechanism can operate even in asymptotically flat spacetimes, by studying the evolution of the Einstein-Klein-Gordon system for a self-interacting scalar field. We show that (i) configurations which do not collapse promptly can do so after successive reflections off the potential barrier, but (ii) that at intermediate amplitudes and Compton wavelengths, collapse to black holes is replaced by the appearance of oscillating soliton stars. Finally, (iii) for very small initial amplitudes, the field disperses away in a manner consistent with power-law tails of massive fields. Minkowski is stable against gravitational collapse. Our results provide one further piece to the rich phenomenology of gravitational collapse and show the important interplay between bound states, blueshift, dissipation and confinement effects.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Higher dimensional Numerical Relativity: code comparison

The nonlinear behavior of higher dimensional black hole spacetimes is of interest in several contexts, ranging from an understanding of cosmic censorship to black hole production in high-energy collisions. However, nonlinear numerical evolutions of higher dimensional black hole spacetimes are tremendously complex, involving different diagnostic tools and "dimensional reduction methods". In this work we compare two different successful codes to evolve Einstein's equations in higher dimensions, and show that the results of such different procedures agree to numerical precision, when applied to the collision from rest of two equal-mass black holes. We calculate the total radiated energy to be E/M=9x10^{-4} in five dimensions and E/M=8.1x10^{-4} in six dimensions.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Light rings as observational evidence for event horizons: long-lived modes, ergoregions and nonlinear instabilities of ultracompact objects

Ultracompact objects are self-gravitating systems with a light ring. It was recently suggested that fluctuations in the background of these objects are extremely long-lived and might turn unstable at the nonlinear level, if the object is not endowed with a horizon. If correct, this result has important consequences: objects with a light ring are black holes. In other words, the nonlinear instability of ultracompact stars would provide a strong argument in favor of the "black hole hypothesis," once electromagnetic or gravitational-wave observations confirm the existence of light rings. Here we explore in some depth the mode structure of ultracompact stars, in particular constant-density stars and gravastars. We show that the existence of very long-lived modes -- localized near a second, stable null geodesic -- is a generic feature of gravitational perturbations of such configurations. Already at the linear level, such modes become unstable if the object rotates sufficiently fast to develop an ergoregion. Finally, we conjecture that the long-lived modes become unstable under fragmentation via a Dyson-Chandrasekhar-Fermi mechanism at the nonlinear level. Depending on the structure of the star, it is also possible that nonlinearities lead to the formation of small black holes close to the stable light ring. Our results suggest that the mere observation of a light ring is a strong evidence for the existence of black holes.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

On the nonlinear instability of confined geometries

The discovery of a "weakly-turbulent" instability of anti-de Sitter spacetime supports the idea that confined fluctuations eventually collapse to black holes and suggests that similar phenomena might be possible in asymptotically-flat spacetime, for example in the context of spherically symmetric oscillations of stars or nonradial pulsations of ultracompact objects. Here we present a detailed study of the evolution of the Einstein-Klein-Gordon system in a cavity, with different types of deformations of the spectrum, including a mass term for the scalar and Neumann conditions at the boundary. We provide numerical evidence that gravitational collapse always occurs, at least for amplitudes that are three orders of magnitude smaller than Choptuik's critical value and corresponding to more than $10^5$ reflections before collapse. The collapse time scales as the inverse square of the initial amplitude in the small-amplitude limit. In addition, we find that fields with nonresonant spectrum collapse earlier than in the fully-resonant case, a result that is at odds with the current understanding of the process. Energy is transferred through a direct cascade to high frequencies when the spectrum is resonant, but we observe both direct- and inverse-cascade effects for nonresonant spectra. Our results indicate that a fully-resonant spectrum might not be a crucial ingredient of the conjectured turbulent instability and that other mechanisms might be relevant. We discuss how a definitive answer to this problem is essentially impossible within the present framework.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

An Alternative Numerical Method for the Stationary Pulsar Magnetosphere

Stationary pulsar magnetospheres in the force-free system are governed by the pulsar equation. In 1999, Contopoulos, Kazanas, and Fendt (hereafter CKF) numerically solved the pulsar equation and obtained a pulsar magnetosphere model called the CKF solution that has both closed and open magnetic field lines. The CKF solution is a successful solution, but it contains a poloidal current sheet that flows along the last open field line. This current sheet is artificially added to make the current system closed. In this paper, we suggest an alternative method to solve the pulsar equation and construct pulsar magnetosphere models without a current sheet. In our method, the pulsar equation is decomposed into Ampere's law and the force-free condition. We numerically solve these equations simultaneously with a fixed poloidal current. As a result, we obtain a pulsar magnetosphere model without a current sheet, which is similar to the CKF solution near the neutron star and has a jet-like structure at a distance along the pole. In addition, we discuss physical properties of the model and find that the force-free condition breaks down in a vicinity of the light cylinder due to dissipation that is included implicitly in the numerical method.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Black Hole Universe: Time Evolution

Time evolution of a black hole lattice universe is simulated. The vacuum Einstein equations in a cubic box with a black hole at the origin are numerically solved with periodic boundary conditions on all pairs of opposite faces. Defining effective scale factors by using the area of a surface and the length of an edge of the cubic box, we compare them with that in the Einstein-deSitter universe. It is found that the behaviour of the effective scale factors is well approximated by that in the Einstein-deSitter universe. Our result suggests that local inhomogeneities do not significantly affect the global expansion law of the universe even if the inhomogeneity is extremely nonlinear.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Coalescence of binary neutron stars in a scalar-tensor theory of gravity

We carry out numerical-relativity simulations of coalescing binary neutron stars in a scalar-tensor theory that admits spontaneous scalarization. We model neutron stars with realistic equations of state. We choose the free parameters of the theory taking into account the constraints imposed by the latest observations of neutron-star-- white-dwarf binaries with pulsar timing. We show that even within those severe constraints, scalarization can still affect the evolution of the binary neutron stars not only during the late inspiral, but also during the merger stage. We also confirm that even when both neutron stars have quite small scalar charge at large separations, they can be strongly scalarized dynamically during the final stages of the inspiral. In particular, we identify the binary parameters for which scalarization occurs either during the late inspiral or only after the onset of the merger when a remnant, supramassive or hypermassive neutron star is formed. We also discuss how those results can impact the extraction of physical information on gravitational waves once they are detected.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Initial Conditions for Numerical Relativity -- Introduction to numerical methods for solving elliptic PDEs

Numerical relativity became a powerful tool to investigate the dynamics of binary problems with black holes or neutron stars as well as the very structure of General Relativity. Although public numerical relativity codes are available to evolve such systems, a proper understanding of the methods involved is quite important. Here we focus on the numerical solution of elliptic partial differential equations. Such equations arise when preparing initial data for numerical relativity, but also for monitoring the evolution of black holes. Because such elliptic equations play an important role in many branches of physics, we give an overview of the topic, and show how to numerically solve them with simple examples and sample codes written in C++ and Fortran90 for beginners in numerical relativity or other fields requiring numerical expertise.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

NR/HEP: roadmap for the future

Physics in curved spacetime describes a multitude of phenomena, ranging from astrophysics to high energy physics. The last few years have witnessed further progress on several fronts, including the accurate numerical evolution of the gravitational field equations, which now allows highly nonlinear phenomena to be tamed. Numerical relativity simulations, originally developed to understand strong field astrophysical processes, could prove extremely useful to understand high-energy physics processes like trans-Planckian scattering and gauge-gravity dualities. We present a concise and comprehensive overview of the state-of-the-art and important open problems in the field(s), along with guidelines for the next years. This writeup is a summary of the "NR/HEP Workshop" held in Madeira, Portugal from August 31st to September 3rd 2011.

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

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

0 citations0 reviewsNo code linkedOpen access