BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Koji Uryu

Koji Uryu contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
gr-qcastro-ph.HEastro-ph

Trust snapshot

Quick read

Trust 21 - Emerging
14works
0followers
3topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

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

14 published item(s)

preprint2015arXiv

New code for quasiequilibrium initial data of binary neutron stars: Corotating, irrotational and slowly spinning systems

We present the extension of our \cocal~- Compact Object CALculator - code to compute general-relativistic initial data for binary compact-star systems. In particular, we construct quasiequilibrium initial data for equal-mass binaries with spins that are either aligned or antialigned with the orbital angular momentum. The Isenberg-Wilson-Mathews formalism is adopted and the constraint equations are solved using the representation formula with a suitable choice of a Green's function. We validate the new code with solutions for equal-mass binaries and explore its capabilities for a wide range of compactnesses, from a white dwarf binary with compactness $\sim 10^{-4}$, up to a highly relativistic neutron-star binary with compactness $\sim 0.22$. We also present a comparison with corotating and irrotational quasiequilibrium sequences from the spectral code \lorene [Taniguchi and Gourgoulhon, Phys. Rev. D {\bf 66}, 104019 (2002)] and with different compactness, showing that the results from the two codes agree to a precision of the order of $0.05\%$. Finally, we present equilibria for spinning configurations with a nuclear-physics equation of state in a piecewise polytropic representation.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Equilibrium solutions of relativistic rotating stars with mixed poloidal and toroidal magnetic fields

Stationary and axisymmetric solutions of relativistic rotating stars with strong mixed poloidal and toroidal magnetic fields are obtained numerically. Because of the mixed components of the magnetic field, the underlying stationary and axisymmetric spacetimes are no longer circular. These configurations are computed from the full set of the Einstein-Maxwell equations, Maxwell's equations and from first integrals and integrability conditions of the magnetohydrodynamic equilibrium equations. After a brief introduction of the formulation of the problem, we present the first results for highly deformed magnetized rotating compact stars.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Binary black hole circular orbits computed with COCAL

In this work we present our first results of binary black hole circular orbits using {\sc cocal}, the Compact Object CALculator. Using the 3+1 decomposition five equations are being solved under the assumptions of conformal flatness and maximal slicing. Excision is used and the appropriate apparent horizon boundary conditions are applied. The orbital velocity is determined by imposing a Schwarzschild behaviour at infinity. A sequence of equal mass black holes is obtained and its main physical characteristics are calculated.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

New code for equilibriums and quasiequilibrium initial data of compact objects. II. Convergence tests and comparisons of binary black hole initial data

COCAL is a code for computing equilibriums or quasiequilibrium initial data of single or binary compact objects based on finite difference methods. We present the results of supplementary convergence tests of COCAL code using time symmetric binary black hole data (Brill-Lindquist solution). Then, we compare the initial data of binary black holes on the conformally flat spatial slice obtained from COCAL and KADATH, where KADATH is a library for solving a wide class of problems in theoretical physics including relativistic compact objects with spectral methods. Data calculated from the two codes converge nicely towards each other, for close as well as largely separated circular orbits of binary black holes. Finally, as an example, a sequence of equal mass binary black hole initial data with corotating spins is calculated and compared with data in the literature.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

A new code for equilibriums and quasiequilibrium initial data of compact objects

We present a new code, named COCAL - Compact Object CALculator, for the computation of equilibriums and quasiequilibrium initial data sets of single or binary compact objects of all kinds. In the cocal code, those solutions are calculated on one or multiple spherical coordinate patches covering the initial hypersurface up to the asymptotic region. The numerical method used to solve field equations written in elliptic form is an adaptation of self-consistent field iterations in which Green's integral formula is computed using multipole expansions and standard finite difference schemes. We extended the method so that it can be used on a computational domain with excised regions for a black hole and a binary companion. Green's functions are constructed for various types of boundary conditions imposed at the surface of the excised regions for black holes. The numerical methods used in cocal are chosen to make the code simpler than any other recent initial data codes, accepting the second order accuracy for the finite difference schemes. We perform convergence tests for time symmetric single black hole data on a single coordinate patch, and binary black hole data on multiple patches. Then, we apply the code to obtain spatially conformally flat binary black hole initial data using boundary conditions including the one based on the existence of equilibrium apparent horizons.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Magnetohydrodynamics in stationary and axisymmetric spacetimes: a fully covariant approach

A fully geometrical treatment of general relativistic magnetohydrodynamics (GRMHD) is developed under the hypotheses of perfect conductivity, stationarity and axisymmetry. The spacetime is not assumed to be circular, which allows for greater generality than the Kerr-type spacetimes usually considered in GRMHD. Expressing the electromagnetic field tensor solely in terms of three scalar fields related to the spacetime symmetries, we generalize previously obtained results in various directions. In particular, we present the first relativistic version of the Soloviev transfield equation, subcases of which lead to fully covariant versions of the Grad-Shafranov equation and of the Stokes equation in the hydrodynamical limit. We have also derived, as another subcase of the relativistic Soloviev equation, the equation governing magnetohydrodynamical equilibria with purely toroidal magnetic fields in stationary and axisymmetric spacetimes.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Neutron star equation of state via gravitational wave observations

Gravitational wave observations can potentially measure properties of neutron star equations of state by measuring departures from the point-particle limit of the gravitational waveform produced in the late inspiral of a neutron star binary. Numerical simulations of inspiraling neutron star binaries computed for equations of state with varying stiffness are compared. As the stars approach their final plunge and merger, the gravitational wave phase accumulates more rapidly if the neutron stars are more compact. This suggests that gravitational wave observations at frequencies around 1 kHz will be able to measure a compactness parameter and place stringent bounds on possible neutron star equations of state. Advanced laser interferometric gravitational wave observatories will be able to tune their frequency band to optimize sensitivity in the required frequency range to make sensitive measures of the late-inspiral phase of the coalescence.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Quasi-equilibrium models of magnetized compact objects

We report work towards a relativistic formulation for modeling strongly magnetized neutron stars, rotating or in a close circular orbit around another neutron star or black hole, under the approximations of helical symmetry and ideal MHD. The quasi-stationary evolution is governed by the first law of thermodynamics for helically symmetric systems, which is generalized to include magnetic fields. The formulation involves an iterative scheme for solving the Einstein-Maxwell and relativistic MHD-Euler equations numerically. The resulting configurations for binary systems could be used as self-consistent initial data for studying their inspiral and merger.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Thermodynamics of magnetized binary compact objects

Binary systems of compact objects with electromagnetic field are modeled by helically symmetric Einstein-Maxwell spacetimes with charged and magnetized perfect fluids. Previously derived thermodynamic laws for helically-symmetric perfect-fluid spacetimes are extended to include the electromagnetic fields, and electric currents and charges; the first law is written as a relation between the change in the asymptotic Noether charge $\dl Q$ and the changes in the area and electric charge of black holes, and in the vorticity, baryon rest mass, entropy, charge and magnetic flux of the magnetized fluid. Using the conservation laws of the circulation of magnetized flow found by Bekenstein and Oron for the ideal magnetohydrodynamic (MHD) fluid, and also for the flow with zero conducting current, we show that, for nearby equilibria that conserve the quantities mentioned above, the relation $\dl Q=0$ is satisfied. We also discuss a formulation for computing numerical solutions of magnetized binary compact objects in equilibrium with emphasis on a first integral of the ideal MHD-Euler equation.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Measuring the neutron star equation of state with gravitational wave observations

We report the results of a first study that uses numerical simulations to estimate the accuracy with which one can use gravitational wave observations of double neutron star inspiral to measure parameters of the neutron-star equation of state. The simulations use the evolution and initial-data codes of Shibata and Uryu to compute the last several orbits and the merger of neutron stars, with matter described by a parametrized equation of state. Previous work suggested the use of an effective cutoff frequency to place constraints on the equation of state. We find, however, that greater accuracy is obtained by measuring departures from the point-particle limit of the gravitational waveform produced during the late inspiral. As the stars approach their final plunge and merger, the gravitational wave phase accumulates more rapidly for smaller values of the neutron star compactness (the ratio of the mass of the neutron star to its radius). We estimate that realistic equations of state will lead to gravitational waveforms that are distinguishable from point particle inspirals at an effective distance (the distance to an optimally oriented and located system that would produce an equivalent waveform amplitude) of 100 Mpc or less. As Lattimer and Prakash observed, neutron-star radius is closely tied to the pressure at density not far above nuclear. Our results suggest that broadband gravitational wave observations at frequencies between 500 and 1000 Hz will constrain this pressure, and we estimate the accuracy with which it can be measured. Related first estimates of radius measurability show that the radius can be determined to an accuracy of ~1 km at 100 Mpc.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Non-conformally flat initial data for binary compact objects

A new method is described for constructing initial data for a binary neutron-star (BNS) system in quasi-equilibrium circular orbit. Two formulations for non-conformally flat data, waveless (WL) and near-zone helically symmetric (NHS), are introduced; in each formulation, the Einstein-Euler system, written in 3+1 form on an asymptotically flat spacelike hypersurface, is exactly solved for all metric components, including the spatially non-conformally flat potentials, and for irrotational flow. A numerical method applicable to both formulations is explained with an emphasis on the imposition of a spatial gauge condition. Results are shown for solution sequences of irrotational BNS with matter approximated by parametrized equations of state that use a few segments of polytropic equations of state. The binding energy and total angular momentum of solution sequences computed within the conformally flat -- Isenberg-Wilson-Mathews (IWM) -- formulation are closer to those of the third post-Newtonian (3PN) two point particles up to the closest orbits, for the more compact stars, whereas sequences resulting from the WL/NHS formulations deviate from the 3PN curve even more for the sequences with larger compactness. We think it likely that this correction reflects an overestimation in the IWM formulation as well as in the 3PN formula, by $\sim 1$ cycle in the gravitational wave phase during the last several orbits. The work suggests that imposing spatial conformal flatness results in an underestimate of the quadrupole deformation of the components of binary neutron-star systems in the last few orbits prior to merger.

0 citations0 reviewsNo code linkedOpen access
preprint2004arXiv

Deriving formulations for numerical computation of binary neutron stars in quasicircular orbits

Two relations, the virial relation $M_{\rm ADM}=M_{\rm K}$ and the first law in the form $δM_{\rm ADM}=ΩδJ$, should be satisfied by a solution and a sequence of solutions describing binary compact objects in quasiequilibrium circular orbits. Here, $M_{\rm ADM}$, $M_{\rm K}$, $J$, and $Ω$ are the ADM mass, Komar mass, angular momentum, and orbital angular velocity, respectively. $δ$ denotes an Eulerian variation. These two conditions restrict the allowed formulations that we may adopt. First, we derive relations between $M_{\rm ADM}$ and $M_{\rm K}$ and between $δM_{\rm ADM}$ and $ΩδJ$ for general asymptotically flat spacetimes. Then, to obtain solutions that satisfy the virial relation and sequences of solutions that satisfy the first law at least approximately, we propose a formulation for computation of quasiequilibrium binary neutron stars in general relativity. In contrast to previous approaches in which a part of the Einstein equation is solved, in the new formulation, the full Einstein equation is solved with maximal slicing and in a transverse gauge for the conformal three-metric. Helical symmetry is imposed in the near zone, while in the distant zone, a waveless condition is assumed. We expect the solutions obtained in this formulation to be excellent quasiequilibria as well as initial data for numerical simulations of binary neutron star mergers.

0 citations0 reviewsNo code linkedOpen access
preprint2003arXiv

Merger of binary neutron stars of unequal mass in full general relativity

We present results of three dimensional numerical simulations of the merger of unequal-mass binary neutron stars in full general relativity. A $Γ$-law equation of state $P=(Γ-1)ρε$ is adopted, where $P$, $ρ$, $\varep$, and $Γ$ are the pressure, rest mass density, specific internal energy, and the adiabatic constant, respectively. We take $Γ=2$ and the baryon rest-mass ratio $Q_M$ to be in the range 0.85--1. The typical grid size is $(633,633,317)$ for $(x,y,z)$ . We improve several implementations since the latest work. In the present code, the radiation reaction of gravitational waves is taken into account with a good accuracy. This fact enables us to follow the coalescence all the way from the late inspiral phase through the merger phase for which the transition is triggered by the radiation reaction. It is found that if the total rest-mass of the system is more than $\sim 1.7$ times of the maximum allowed rest-mass of spherical neutron stars, a black hole is formed after the merger irrespective of the mass ratios. The gravitational waveforms and outcomes in the merger of unequal-mass binaries are compared with those in equal-mass binaries. It is found that the disk mass around the so formed black holes increases with decreasing rest-mass ratios and decreases with increasing compactness of neutron stars. The merger process and the gravitational waveforms also depend strongly on the rest-mass ratios even for the range $Q_M= 0.85$--1.

0 citations0 reviewsNo code linkedOpen access
preprint2001arXiv

Thermodynamics of binary black holes and neutron stars

We consider compact binary systems, modeled in general relativity as vacuum or perfect-fluid spacetimes with a helical Killing vector k^α, heuristically, the generator of time-translations in a corotating frame. Systems that are stationary in this sense are not asymptotically flat, but have asymptotic behavior corresponding to equal amounts of ingoing and outgoing radiation. For black-hole binaries, a rigidity theorem implies that the Killing vector lies along the horizon's generators, and from this one can deduce the zeroth law (constant surface gravity of the horizon). Remarkably, although the mass and angular momentum of such a system are not defined, there is an exact first law, relating the change in the asymptotic Noether charge to the changes in the vorticity, baryon mass, and entropy of the fluid, and in the area of black holes. Binary systems with MΩsmall have an approximate asymptopia in which one can write the first law in terms of the asymptotic mass and angular momentum. Asymptotic flatness is precise in two classes of solutions used to model binary systems: spacetimes satisfying the post-Newtonian equations, and solutions to a modified set of field equations that have a spatially conformally flat metric. (The spatial conformal flatness formalism with helical symmetry, however, is consistent with maximal slicing only if replaces the extrinsic curvature in the field equations by an artificially tracefree expression in terms of the shift vector.) For these spacetimes, nearby equilibria whose stars have the same vorticity obey the relation δM = ΩδJ, from which one can obtain a turning point criterion that governs the stability of orbits.

0 citations0 reviewsNo code linkedOpen access