Post-Newtonian SPH calculations of binary neutron star coalescence. I. Method and first results

We present the first results from our Post-Newtonian (PN) Smoothed Particle Hydrodynamics (SPH) code, which has been used to study the coalescence of binary neutron star (NS) systems. The Lagrangian particle-based code incorporates consistently all lowest-order (1PN) relativistic effects, as well as gravitational radiation reaction, the lowest-order dissipative term in general relativity. We test our code on sequences of single NS models of varying compactness, and we discuss ways to make PN simulations more relevant to realistic NS models. We also present a PN SPH relaxation procedure for constructing equilibrium models of synchronized binaries, and we use these equilibrium models as initial conditions for our dynamical calculations of binary coalescence. Though unphysical, since tidal synchronization is not expected in NS binaries, these initial conditions allow us to compare our PN work with previous Newtonian results. We compare calculations with and without 1PN effects, for NS with stiff equations of state, modeled as polytropes with $\Gamma=3$. We find that 1PN effects can play a major role in the coalescence, accelerating the final inspiral and causing a significant misalignment in the binary just prior to final merging. In addition, the character of the gravitational wave signal is altered dramatically, showing strong modulation of the exponentially decaying waveform near the end of the merger. We also discuss briefly the implications of our results for models of gamma-ray bursts at cosmological distances.Comment: RevTeX, 37 pages, 17 figures, to appear in Phys. Rev. D, minor corrections onl

Graviton-induced Bremsstrahlung

We discuss photon Bremsstrahlung induced by virtual graviton exchange in proton-proton interactions at hadronic colliders, resulting from the exchange of Kaluza-Klein excitations of the graviton. The relevant subprocesses, gg to G to e^+e^-gamma and q barq to e^+e^-gamma are discussed in both the ADD and the RS scenarios. Although two-body final states (or real graviton emission) would presumably be the main discovery channels, a search for three-body final states could be worthwhile since such events have characteristic features that could provide additional confirmation. In particular, the k_perp-distribution of the photon is in both scenarios harder than that of the Standard-model background.Comment: 24 pages, including figures. v2: Including initial-state Bremsstrahlung and photon k_perp spectra. Version to appear in PR

Post-Newtonian SPH calculations of binary neutron star coalescence. II. Binary mass ratio, equation of state, and spin dependence

Using our new Post-Newtonian SPH (smoothed particle hydrodynamics) code, we study the final coalescence and merging of neutron star (NS) binaries. We vary the stiffness of the equation of state (EOS) as well as the initial binary mass ratio and stellar spins. Results are compared to those of Newtonian calculations, with and without the inclusion of the gravitational radiation reaction. We find a much steeper decrease in the gravity wave peak strain and luminosity with decreasing mass ratio than would be predicted by simple point-mass formulae. For NS with softer EOS (which we model as simple $\Gamma=2$ polytropes) we find a stronger gravity wave emission, with a different morphology than for stiffer EOS (modeled as $\Gamma=3$ polytropes as in our previous work). We also calculate the coalescence of NS binaries with an irrotational initial condition, and find that the gravity wave signal is relatively suppressed compared to the synchronized case, but shows a very significant second peak of emission. Mass shedding is also greatly reduced, and occurs via a different mechanism than in the synchronized case. We discuss the implications of our results for gravity wave astronomy with laser interferometers such as LIGO, and for theoretical models of gamma-ray bursts (GRBs) based on NS mergers.Comment: RevTeX, 38 pages, 24 figures, Minor Corrections, to appear in Phys. Rev.

Irrotational binary neutron stars in quasiequilibrium

We report on numerical results from an independent formalism to describe the quasi-equilibrium structure of nonsynchronous binary neutron stars in general relativity. This is an important independent test of controversial numerical hydrodynamic simulations which suggested that nonsynchronous neutron stars in a close binary can experience compression prior to the last stable circular orbit. We show that, for compact enough stars the interior density increases slightly as irrotational binary neutron stars approach their last orbits. The magnitude of the effect, however, is much smaller than that reported in previous hydrodynamic simulations.Comment: 4 pages, 2 figures, revtex, accepted for publication in Phys. Rev.

Next-to-Leading Order Cross Sections for Tagged Reactions

We extend the phase space slicing method of Giele, Glover and Kosower for performing next-to-leading order jet cross section calculations in two important ways: we show how to include fragmentation functions and how to include massive particles. These extensions allow the application of this method to not just jet cross sections but also to cross sections in which a particular final state particle, including a $D$ or $B$-meson, is tagged.Comment: 36 pages, Latex Small corrections to text. To appear in Phys. Rev.

Towards a Realistic Neutron Star Binary Inspiral: Initial Data and Multiple Orbit Evolution in Full General Relativity

This paper reports on our effort in modeling realistic astrophysical neutron star binaries in general relativity. We analyze under what conditions the conformally flat quasiequilibrium (CFQE) approach can generate ``astrophysically relevant'' initial data, by developing an analysis that determines the violation of the CFQE approximation in the evolution of the binary described by the full Einstein theory. We show that the CFQE assumptions significantly violate the Einstein field equations for corotating neutron stars at orbital separations nearly double that of the innermost stable circular orbit (ISCO) separation, thus calling into question the astrophysical relevance of the ISCO determined in the CFQE approach. With the need to start numerical simulations at large orbital separation in mind, we push for stable and long term integrations of the full Einstein equations for the binary neutron star system. We demonstrate the stability of our numerical treatment and analyze the stringent requirements on resolution and size of the computational domain for an accurate simulation of the system.Comment: 22 pages, 18 figures, accepted to Phys. Rev.

Hydrologic Modeling and Delineation of Calumpang River Watershed using GIS and Hydrologic Model System

The Calumpang Bridge in CALABARZON Region, Southern Luzon Philippines did not withstand the river current brought by Super Typhoon Glenda (Rammasun) on July 2014, causing the destruction of one-third of the bridge, thus resulting to traffic congestion and economic distress for the residents and businessmen. This paper aims to create a calibrated hydrological model. Specifically, it focuses on delineating watershed, simulation based on observed flow data and validation that will be the basis for flood risk map for the river. The watershed has been delineated and it is calibrated, the basin was modeled using Hydrologic Model System and GIS Applications to determine hydrologic parameters which have spatial characteristic and to compute the peak Discharge and loss infiltration. The model delineated by the use of GIS software has been hydrologically corrected. The model describes the correlation of the rainfall, losses, time of concentration, Storage factor and amount of discharge. It shows that the amount of Rainfall at a given time dictates the amount of discharge the watershed generates. The higher the rainfall amount, the higher the discharge. The model was validated using Percentage Error and Rational Method, and due to the accuracy of the validated model, it provides a promising approach to Bridge Structural Design and Flood Risk Map Generation

Close encounters of a rotating star with planets in parabolic orbits of varying inclination and the formation of Hot Jupiters

(abbreviated) We extend the theory of close encounters of a planet on a parabolic orbit with a star to include the effects of tides induced on the central rotating star. Orbits with arbitrary inclination to the stellar rotation axis are considered. We obtain results both from an analytic treatment and numerical one that are in satisfactory agreement. These results are applied to the initial phase of the tidal circularisation problem. We find that both tides induced in the star and planet can lead to a significant decrease of the orbital semi-major axis for orbits having periastron distances smaller than 5-6 stellar radii (corresponding to periods $\sim 4-5$ days after the circularisation has been completed) with tides in the star being much stronger for retrograde orbits compared to prograde orbits. We use the simple Skumanich law for the stellar rotation with its rotational period equal to one month at the age of 5Gyr. The strength of tidal interactions is characterised by circularisation time scale, $t_{ev}$ defined as a time scale of evolution of the planet's semi-major axis due to tides considered as a function of orbital period $P_{obs}$ after the process of tidal circularisation has been completed. We find that the ratio of the initial circularisation time scales corresponding to prograde and retrograde orbits is of order 1.5-2 for a planet of one Jupiter mass and $P_{obs}\sim$ four days. It grows with the mass of the planet, being of order five for a five Jupiter mass planet with the same $P_{orb}$. Thus, the effect of stellar rotation may provide a bias in the formation of planetary systems having planets on close orbits around their host stars, as a consequence of planet-planet scattering, favouring systems with retrograde orbits. The results may also be applied to the problem of tidal capture of stars in young stellar clusters.Comment: to be published in Celestial Mechanics and Dynamical Astronom

Associated Higgs production with top quarks at the Large Hadron Collider: NLO QCD corrections

We present in detail the calculation of the O(alpha_s^3) inclusive total cross section for the process pp -> t-tbar-h, in the Standard Model, at the CERN Large Hadron Collider with center-of-mass energy sqrt(s_H)=14 TeV. The calculation is based on the complete set of virtual and real O(alpha_s) corrections to the parton level processes q-qbar -> t-tbar-h and gg -> t-tbar-h, as well as the tree level processes (q,qbar)g -> t-tbar-h-(q,qbar). The virtual corrections involve the computation of pentagon diagrams with several internal and external massive particles, first encountered in this process. The real corrections are computed using both the single and the two cutoff phase space slicing method. The next-to-leading order QCD corrections significantly reduce the renormalization and factorization scale dependence of the Born cross section and moderately increase the Born cross section for values of the renormalization and factorization scales above m_t.Comment: 70 pages, 12 figures, RevTeX4: one word changed in the abstract, one sentence reworded in the introduction. To appear in Phys. Rev.