A multivariate variational objective analysis-assimilation method. Part 2: Case study results with and without satellite data

The variational multivariate assimilation method described in a companion paper by Achtemeier and Ochs is applied to conventional and conventional plus satellite data. Ground-based and space-based meteorological data are weighted according to the respective measurement errors and blended into a data set that is a solution of numerical forms of the two nonlinear horizontal momentum equations, the hydrostatic equation, and an integrated continuity equation for a dry atmosphere. The analyses serve first, to evaluate the accuracy of the model, and second to contrast the analyses with and without satellite data. Evaluation criteria measure the extent to which: (1) the assimilated fields satisfy the dynamical constraints, (2) the assimilated fields depart from the observations, and (3) the assimilated fields are judged to be realistic through pattern analysis. The last criterion requires that the signs, magnitudes, and patterns of the hypersensitive vertical velocity and local tendencies of the horizontal velocity components be physically consistent with respect to the larger scale weather systems

Toward stable 3D numerical evolutions of black-hole spacetimes

Three dimensional (3D) numerical evolutions of static black holes with excision are presented. These evolutions extend to about 8000M, where M is the mass of the black hole. This degree of stability is achieved by using growth-rate estimates to guide the fine tuning of the parameters in a multi-parameter family of symmetric hyperbolic representations of the Einstein evolution equations. These evolutions were performed using a fixed gauge in order to separate the intrinsic stability of the evolution equations from the effects of stability-enhancing gauge choices.Comment: 4 pages, 5 figures. To appear in Phys. Rev. D. Minor additions to text for clarification. Added short paragraph about inner boundary dependenc

Spin effects in gravitational radiation backreaction II. Finite mass effects

A convenient formalism for averaging the losses produced by gravitational radiation backreaction over one orbital period was developed in an earlier paper. In the present paper we generalize this formalism to include the case of a closed system composed from two bodies of comparable masses, one of them having the spin S. We employ the equations of motion given by Barker and O'Connell, where terms up to linear order in the spin (the spin-orbit interaction terms) are kept. To obtain the radiative losses up to terms linear in the spin, the equations of motion are taken to the same order. Then the magnitude L of the angular momentum L, the angle kappa subtended by S and L and the energy E are conserved. The analysis of the radial motion leads to a new parametrization of the orbit. From the instantaneous gravitational radiation losses computed by Kidder the leading terms and the spin-orbit terms are taken. Following Apostolatos, Cutler, Sussman and Thorne, the evolution of the vectors S and L in the momentary plane spanned by these vectors is separated from the evolution of the plane in space. The radiation-induced change in the spin is smaller than the leading-order spin terms in the momentary angular momentum loss. This enables us to compute the averaged losses in the constants of motion E, L and L_S=L cos kappa. In the latter, the radiative spin loss terms average to zero. An alternative description using the orbital elements a,e and kappa is given. The finite mass effects contribute terms, comparable in magnitude, to the basic, test-particle spin terms in the averaged losses.Comment: 12 pages, 1 figure, Phys.Rev.D15, March, 199

Variational objective analyses for cyclone studies

The basic analysis equations, i.e., the two horizontal momentum equations, the hydrostatic equation, and the integrated continuity equation were derived for the nonlinear vertical coordinate, nondimensionalized, and expressed in finite differences on a staggered grid. Special care was taken to transform the hydrostatic equation and the pressure gradient terms of the horizontal momentum equations to nearly eliminate truncation error over steeply sloping terrain. This formulation also eliminated explicit reference to orographically induced variations in the thermodynamic variables so that the variational adjustments are on the scale of the meteorological perturbations. The analysis equations were subjected to the Euler-Lagrange operations as expressed for finite differences and an additional set of five partial differential equations was derived, bringing to nine the number of equations in Model I. Higher order terms, terms containing observed quantities, and terms containing none of the variables to be adjusted were grouped into forcing functions and the equations were solved for the zero order terms. Zero order variables were eliminated between these equations and there resulted two diagnostic equations which take the form of general linear second order partial differential equations with nonconstant coefficients

Introductory lectures on the Effective One Body formalism

The Effective One Body (EOB) formalism is an analytical approach which aims at providing an accurate description of the motion and radiation of coalescing binary black holes. We present a brief review of the basic elements of this approach.Comment: 22 pages, 3 figures, lectures given at the Second ICRANet Stueckelberg Workshop on Relativistic Field Theories (Pescara, Italy, September 3-8, 2007); to be published in the International Journal of Modern Physics

Numerical simulations of neutron star-black hole binaries in the near-equal-mass regime

Simulations of neutron star-black hole (NSBH) binaries generally consider black holes with masses in the range $(5-10)M_\odot$, where we expect to find most stellar mass black holes. The existence of lower mass black holes, however, cannot be theoretically ruled out. Low-mass black holes in binary systems with a neutron star companion could mimic neutron star-neutron (NSNS) binaries, as they power similar gravitational wave (GW) and electromagnetic (EM) signals. To understand the differences and similarities between NSNS mergers and low-mass NSBH mergers, numerical simulations are required. Here, we perform a set of simulations of low-mass NSBH mergers, including systems compatible with GW170817. Our simulations use a composition and temperature dependent equation of state (DD2) and approximate neutrino transport, but no magnetic fields. We find that low-mass NSBH mergers produce remnant disks significantly less massive than previously expected, and consistent with the post-merger outflow mass inferred from GW170817 for moderately asymmetric mass ratio. The dynamical ejecta produced by systems compatible with GW170817 is negligible except if the mass ratio and black hole spin are at the edge of the allowed parameter space. That dynamical ejecta is cold, neutron-rich, and surprisingly slow for ejecta produced during the tidal disruption of a neutron star : $v\sim (0.1-0.15)c$. We also find that the final mass of the remnant black hole is consistent with existing analytical predictions, while the final spin of that black hole is noticeably larger than expected -- up to $\chi_{\rm BH}=0.84$ for our equal mass case

Gravitational radiation reaction in compact binary systems: Contribution of the quadrupole-monopole interaction

The radiation reaction in compact spinning binaries on eccentric orbits due to the quadrupole-monopole interaction is studied. This contribution is of second post-Newtonian order. As result of the precession of spins the magnitude $L$ of the orbital angular momentum is not conserved. Therefore a proper characterization of the perturbed radial motion is provided by the energy $E$ and angular average $\bar{L}$. As powerful computing tools, the generalized true and eccentric anomaly parametrizations are introduced. Then the secular losses in energy and magnitude of orbital angular momentum together with the secular evolution of the relative orientations of the orbital angular momentum and spins are found for eccentric orbits by use of the residue theorem. The circular orbit limit of the energy loss agrees with Poisson's earlier result.Comment: accepted for publication in Phys. Rev.

Spin-spin effects in radiating compact binaries

The dynamics of a binary system with two spinning components on an eccentric orbit is studied, with the inclusion of the spin-spin interaction terms appearing at the second post-Newtonian order. A generalized true anomaly parametrization properly describes the radial component of the motion. The average over one radial period of the magnitude of the orbital angular momentum $\bar{L}$ is found to have no nonradiative secular change. All spin-spin terms in the secular radiative loss of the energy and magnitude of orbital angular momentum are given in terms of $\bar{L}$ and other constants of the motion. Among them, self-interaction spin effects are found, representing the second post-Newtonian correction to the 3/2 post-Newtonian order Lense-Thirring approximation.Comment: 12 pages, to appear in Phys. Rev.

Solving the Darwin problem in the first post-Newtonian approximation of general relativity

We analytically calculate the equilibrium sequence of the corotating binary stars of incompressible fluid in the first post-Newtonian(PN) approximation of general relativity. By calculating the total energy and total angular momentum of the system as a function of the orbital separation, we investigate the innermost stable circular orbit for corotating binary(we call it ISCCO). It is found that by the first PN effect, the orbital separation of the binary at the ISCCO becomes small with increase of the compactness of each star, and as a result, the orbital angular velocity at the ISCCO increases. These behaviors agree with previous numerical works.Comment: 33 pages, revtex, 4 figures(eps), accepted for publication in Phys. Rev.