Third post-Newtonian angular momentum flux and the secular evolution of orbital elements for inspiralling compact binaries in quasi-elliptical orbits

The angular-momentum flux from an inspiralling binary system of compact objects moving in quasi-elliptical orbits is computed at the third post-Newtonian (3PN) order using the multipolar post-Minkowskian wave generation formalism. The 3PN angular-momentum flux involves the instantaneous, tail, and tail-of-tails contributions as for the 3PN energy flux, and in addition a contribution due to nonlinear memory. We average the angular-momentum flux over the binary's orbit using the 3PN quasi-Keplerian representation of elliptical orbits. The averaged angular-momentum flux provides the final input needed for gravitational-wave phasing of binaries moving in quasi-elliptical orbits. We obtain the evolution of orbital elements under 3PN gravitational radiation reaction in the quasi-elliptic case. For small eccentricities, we give simpler limiting expressions relevant for phasing up to order e(2). This work is important for the construction of templates for quasi-eccentric binaries, and for the comparison of post-Newtonian results with the numerical relativity simulations of the plunge and merger of eccentric binaries.

The third post-Newtonian gravitational wave polarizations and associated spherical harmonic modes for inspiralling compact binaries in quasi-circular orbits

The gravitational waveform (GWF) generated by inspiralling compact binaries moving in quasi-circular orbits is computed at the third post-Newtonian (3PN) approximation to general relativity. Our motivation is two-fold: (i) to provide accurate templates for the data analysis of gravitational wave inspiral signals in laser interferometric detectors; (ii) to provide the associated spin-weighted spherical harmonic decomposition to facilitate comparison and match of the high post-Newtonian prediction for the inspiral waveform to the numerically-generated waveforms for the merger and ringdown. This extension of the GWF by half a PN order (with respect to previous work at 2.5PN order) is based on the algorithm of the multipolar post-Minkowskian formalism, and mandates the computation of the relations between the radiative, canonical and source multipole moments for general sources at 3PN order. We also obtain the 3PN extension of the source multipole moments in the case of compact binaries, and compute the contributions of hereditary terms (tails, tails-of-tails and memory integrals) up to 3PN order. The end results are given for both the complete plus and cross polarizations and the separate spin-weighted spherical harmonic modes.

Possible use of self-calibration to reduce systematic uncertainties in determining distance-redshift relation via gravitational radiation from merging binaries

By observing mergers of compact objects, future gravity wave experiments would measure the luminosity distance to a large number of sources to a high precision but not their redshifts. Given the directional sensitivity of an experiment, a fraction of such sources (gold plated) can be identified optically as single objects in the direction of the source. We show that if an approximate distance-redshift relation is known then it is possible to statistically resolve those sources that have multiple galaxies in the beam. We study the feasibility of using gold plated sources to iteratively resolve the unresolved sources, obtain the self-calibrated best possible distance-redshift relation and provide an analytical expression for the accuracy achievable. We derive the lower limit on the total number of sources that is needed to achieve this accuracy through self-calibration. We show that this limit depends exponentially on the beam width and give estimates for various experimental parameters representative of future gravitational wave experiments DECIGO and BBO.

2.5PN linear momentum flux from inspiralling compact binaries in quasicircular orbits and associated recoil: Nonspinning case

Anisotropic emission of gravitational waves (GWs) from inspiralling compact binaries leads to the loss of linear momentum and hence gravitational recoil of the system. The loss rate of linear momentum in the far-zone of the source (a nonspinning binary system of black holes in quasicircular orbit) is investigated at the 2.5 post-Newtonian (PN) order and used to provide an analytical expression in harmonic coordinates for the 2.5PN accurate recoil velocity of the binary accumulated in the inspiral phase. The maximum recoil velocity of the binary system at the end of its inspiral phase (i.e at the innermost stable circular orbit (ISCO)) estimated by the 2.5PN formula is of the order of 4 km s(-1) which is smaller than the 2PN estimate of 22 km s(-1). Going beyond inspiral, we also provide an estimate of the more important contribution to the recoil velocity from the plunge phase. The maximum recoil velocity at the end of the plunge, involving contributions both from inspiral and plunge phase, for a binary with symmetric mass ratio nu = 0.2 is of the order of 182 km s(-1).

A Suzaku view of IGR J16393-4643

The pulsar IGR J16393-4643 belongs to a class of highly absorbed supergiant high-mass X-ray binaries (HMXBs), characterized by a very high column density of absorbing matter. We present the results of simultaneous broad-band pulsation and spectrum analysis from a 44-ks Suzaku observation of the source. The orbital intensity profile created with the Swift Burst Alert Telescope (Swift-BAT) light curve shows an indication of IGR J16393-4643 being an eclipsing system with a short eclipse semi-angle theta(E) similar to 17 degrees. For a supergiant companion star with a 20-R-circle dot radius, this implies an inclination of the orbital plane in the range 39 degrees-57 degrees, whereas for a main-sequence B star as the companion with a 10-R-circle dot radius, the inclination of the orbital plane is in the range 60 degrees-77 degrees. Pulse profiles created for different energy bands have complex morphology, which shows some energy dependence and increases in pulse fraction with energy. We have also investigated broad-band spectral characteristics, phase-averaged spectra and resolving the pulse phase into peak and trough phases. The phase-averaged spectrum has a very high N-H(similar to 3 x 10(23) cm(-2)) and is described by a power law (Gamma similar to 0.9) with a high-energy cut-off above 20 keV. We find a change in the spectral index in the peak and trough phases, implying an underlying change in the source spectrum.

Massive black-hole binary inspirals: results from the LISA parameter estimation taskforce

The LISA Parameter Estimation Taskforce was formed in September 2007 to provide the LISA Project with vetted codes, source distribution models and results related to parameter estimation. The Taskforce's goal is to be able to quickly calculate the impact of any mission design changes on LISA's science capabilities, based on reasonable estimates of the distribution of astrophysical sources in the universe. This paper describes our Taskforce's work on massive black-hole binaries (MBHBs). Given present uncertainties in the formation history of MBHBs, we adopt four different population models, based on (i) whether the initial black-hole seeds are small or large and (ii) whether accretion is efficient or inefficient at spinning up the holes. We compare four largely independent codes for calculating LISA's parameter-estimation capabilities. All codes are based on the Fisher-matrix approximation, but in the past they used somewhat different signal models, source parametrizations and noise curves. We show that once these differences are removed, the four codes give results in extremely close agreement with each other. Using a code that includes both spin precession and higher harmonics in the gravitational-wave signal, we carry out Monte Carlo simulations and determine the number of events that can be detected and accurately localized in our four population models.

"Central Atoms" models for ternary silicate and alumino-silicate melts

The �Central Atoms� model presented by the authors in an earlier paper is extended to ternary silicate and alumino-silicate melts. The model is applied to the CaO-FeO-SiO2 and the CaO-Al2O3-SiO2 system. Use is made of the parameters from the relevant binaries only. The agreement between experimental and calculated isoactivity curves is good in all cases.

Description of accretion induced outflows from ultra-luminous sources to under-luminous AGNs

We study the energetics of the accretion-induced outflow and then plausible jet around black holes/compact objects using a newly developed disc-outflow coupled model. Inter-connecting dynamics of outflow and accretion essentially upholds the conservation laws. The energetics depend strongly on the viscosity parameter α and the cooling factor f which exhibit several interesting features. The bolometric luminosities of ultra-luminous X-ray binaries (e.g. SS433) and family of highly luminous AGNs and quasars can be reproduced by the model under the super-Eddington accretion flows. Under appropriate conditions, low-luminous AGNs (e.g. Sagittarius A*) also fit reasonably well with the luminosity corresponding to a sub-Eddington accretion flow with f→1.

Transition from radiatively inefficient to cooling dominated phase in two temperature accretion disks around black holes

We investigate the transition of a radiatively inefficient phase of a viscous two temperature accreting flow to a cooling dominated phase and vice versa around black holes. Based on a global sub-Keplerian accretion disk model in steady state, including explicit cooling processes self-consistently, we show that general advective accretion flow passes through various phases during its infall towards a black hole. Bremsstrahlung, synchrotron and inverse Comptonization of soft photons are considered as possible cooling mechanisms. Hence the flow governs a much lower electron temperature similar to 10(8) - 10(9.5) K compared to the hot protons of temperature similar to 10(10.2) - 10(11.8) K in the range of the accretion rate in Eddington units 0.01 less than or simiar to (M) over dot less than or similar to 100. Therefore, the solutions may potentially explain the hard X-rays and the gamma-rays emitted from AGNs and X-ray binaries. We finally compare the solutions for two different regimes of viscosity and conclude that a weakly viscous flow is expected to be cooling dominated compared to its highly viscous counterpart which is radiatively inefficient. The flow is successfully able to reproduce the observed minosities of the under-fed AGNs and quasars (e.g. Sgr A*), ultra-luminous X-ray sources (e.g. SS433), as well as the highly luminous AGNs and ultra-luminous quasars (e.g. PKS 0743-67) at different combinations of the mass accretion rate and ratio of specific heats.

Effect of pulse profile variations on measurement of eccentricity in orbits of Cen X-3 and SMC X-1

It has long been argued that better timing precision allowed by satellites like Rossi X-ray Timing Explorer (RXTE) will allow us to measure the orbital eccentricity and the angle of periastron of some of the bright persistent high-mass X-ray binaries (HMXBs) and hence a possible measurement of apsidal motion in these system. Measuring the rate of apsidal motion allows one to estimate the apsidal motion constant of the mass losing companion star and hence allows for the direct testing of the stellar structure models for these giant stars present in the HMXBs. In the present paper, we use the archival RXTE data of two bright persistent sources, namely Cen X-3 and SMC X-1, to measure the very small orbital eccentricity and the angle of periastron. We find that the small variations in the pulse profiles of these sources, rather than the intrinsic time resolution provided by RXTE, limit the accuracy with which we can measure arrival time of the pulses from these sources. This influences the accuracy with which one can measure the orbital parameters, especially the very small eccentricity and the angle of periastron in these sources. The observations of SMC X-1 in the year 2000 were taken during the high-flux state of the source and we could determine the orbital eccentricity and omega using this data set.

Two-temperature accretion around rotating black holes: a description of the general advective flow paradigm in the presence of various cooling processes to explain low to high luminous sources

We investigate viscous two-temperature accretion disc flows around rotating black holes. We describe the global solution of accretion flows with a sub-Keplerian angular momentum profile, by solving the underlying conservation equations including explicit cooling processes self-consistently. Bremsstrahlung, synchrotron and inverse Comptonization of soft photons are considered as possible cooling mechanisms. We focus on the set of solutions for sub-Eddington, Eddington and super-Eddington mass accretion rates around Schwarzschild and Kerr black holes with a Kerr parameter of 0.998. It is found that the flow, during its infall from the Keplerian to sub-Kepleria transition region to the black hole event horizon, passes through various phases of advection: the general advective paradigm to the radiatively inefficient phase, and vice versa. Hence, the flow governs a much lower electron temperature similar to 10(8)-10(9.5) K, in the range of accretion rate in Eddington units 0.01 less than or similar to (M) over dot less than or similar to 100, compared to the hot protons of temperature similar to 10(10.2)-10(11.8) K. Therefore, the solution may potentially explain the hard X-rays and gamma-rays emitted from active galactic nuclei (AGNs) and X-ray binaries. We then compare the solutions for two different regimes of viscosity. We conclude that a weakly viscous flow is expected to be cooling dominated, particularly at the inner region of the disc, compared to its highly viscous counterpart, which is radiatively inefficient. With all the solutions in hand, we finally reproduce the observed luminosities of the underfed AGNs and quasars (e. g. Sgr A*) to ultraluminous X-ray sources (e. g. SS433), at different combinations of input parameters, such as the mass accretion rate and the ratio of specific heats. The set of solutions also predicts appropriately the luminosity observed in highly luminous AGNs and ultraluminous quasars (e. g. PKS 0743-67).

Microarchitecture sensitive empirical models for compiler optimizations.

This paper proposes the use of empirical modeling techniques for building microarchitecture sensitive models for compiler optimizations. The models we build relate program performance to settings of compiler optimization flags, associated heuristics and key microarchitectural parameters. Unlike traditional analytical modeling methods, this relationship is learned entirely from data obtained by measuring performance at a small number of carefully selected compiler/microarchitecture configurations. We evaluate three different learning techniques in this context viz. linear regression, adaptive regression splines and radial basis function networks. We use the generated models to a) predict program performance at arbitrary compiler/microarchitecture configurations, b) quantify the significance of complex interactions between optimizations and the microarchitecture, and c) efficiently search for'optimal' settings of optimization flags and heuristics for any given microarchitectural configuration. Our evaluation using benchmarks from the SPEC CPU2000 suits suggests that accurate models (< 5% average error in prediction) can be generated using a reasonable number of simulations. We also find that using compiler settings prescribed by a model-based search can improve program performance by as much as 19% (with an average of 9.5%) over highly optimized binaries.

Orbital evolution studies of two HMXB pulsars

High mass X-ray binary (HMXB) pulsars are of two types, persistent and transient. 4U1538-52 is a persistent HMXB whose orbit was previously measured to be circular but the RXTE observations revealed an eccentric orbit. We observed this system with RXTE-PCA in August 2003 and our timing analysis supports the eccentric orbit of the system. However, we do not find any evidence for orbital evolution. Rotational and tidal interactions between the stars of a closed binary system result in apsidal motion which can be measured in systems with eccentric orbit. 4U0115+63 is a Be-transient HMXB whose eccentric orbit was well-determined during its 1978 outburst. We report preliminary results from analysis of data obtained during the 1999 outburst of this source with the RXTE-PCA.

Apsidal motion in 4U 0115+63 and orbital parameters of 2S 1417-624 and V0332+53

Be/X-ray binary pulsars have wide eccentric orbits and hence the angle of periastron of the orbit is very well defined in these sources. The presence of an X-ray pulsar allows for accurate measurements of orbital elements. A Be star usually is a rapidly rotating star and hence will deviate from spherical geometry. The tidal interaction between the neutron star and the Be star will add to the distortion of the Be star and alter its mass distribution. Thus a measurable rate of apsidal motion is expected from these systems. In this paper, we present the first conclusive detection of apsidal motion of the binary 4U 0115+63. We also present new and accurate orbital parameters of the Be/X-ray binaries V0332+53 and 2S 1417-624.

Computer estimation of the Cd---Hg---Te phase diagram

The compounds CdHgTe and its constituent binaries CdTe, HgTe, and CdHg are semiconductors which are used in thermal, infrared, nuclear, thermoelectric and other photo sensitive devices. The compound CdHgTe has a Sphaleritic structure of possible type A1IIB1IIC6VI. The TERCP program of Kaufman is used to estimate the stable regions of the ternary phase diagram using available thermodynamic data. It was found that there was little variation in stochiometry with temperature. The compositions were calculated for temperatures ranging from 325K to 100K and the compositional limits were Cd13−20Hg12−01Te75−79, Hg varying most. By comparison with a similar compound, Cd In2Te4 of forbidden band width. 88 to .90 e.V., similar properties are postulated for Cd1Hg1Te6 with applications in the infra red region of the spectrum at 300K where this composition is given by TERCP at the limit of stability.