We construct compact and high accuracy Reduced Basis (RB) representations of single and multiple quasinormal modes (QNMs). The RB method determines a hierarchical and relatively small set of the most relevant waveforms. We find that the exponential convergence of the method allows for a dramatic compression of template banks used for ringdown searches. Compressing a catalog with a minimal match $\MMm=0.99$, we find that the selected RB waveforms are able to represent {\em any} QNM, including those not in the original bank, with extremely high accuracy, typically less than $10^{-13}$. We then extend our studies to two-mode QNMs. Inclusion of a second mode is expected to help with detection, and might make it possible to infer details of the progenitor of the final black hole. We find that the number of RB waveforms needed to represent any two-mode ringdown waveform with the above high accuracy is {\em smaller} than the number of metric-based, one-mode templates with $\MMm=0.99$. For unconstrained two-modes, which would allow for consistency tests of General Relativity, our high accuracy RB has around $10^4$ {\em fewer} waveforms than the number of metric-based templates for $\MMm=0.99$. The number of RB elements grows only linearly with the number of multipole modes versus exponentially with the standard approach, resulting in very compact representations even for many multiple modes. The results of this paper open the possibility of searches of multi-mode ringdown gravitational waves.

We introduce a reduced basis approach as a new paradigm for modeling, representing and searching for gravitational waves. We construct waveform catalogs for non-spinning compact binary coalescences, and we find that for accuracies of 99% and 99.999% the method generates a factor of about $10-10^5$ fewer templates than standard placement methods. The continuum of gravitational waves can be represented by a finite and comparatively compact basis. The method is robust under variations in the noise of detectors, implying that only a single catalog needs to be generated.

Disks around young stars are known to evolve from optically thick, gas-dominated protoplanetary disks to optically thin, almost gas-free debris disks. It is thought that the primordial gas is largely removed at ages of ~10 Myr, but it is difficult to discern the true gas densities from gas observations. This suggests using observations of dust: it has been argued that gas, if present with higher densities, would lead to flatter radial profiles of the dust density and surface brightness than those actually observed. However, here we show that these profiles are surprisingly insensitive to variation of the parameters of a central star, location of the dust-producing planetesimal belt, dustiness of the disk and - most importantly - the parameters of the ambient gas. This result holds for a wide range of gas densities (three orders of magnitude), for different radial distributions of the gas temperature, and different gas compositions. The brightness profile slopes of -3...-4 we find are the same that were theoretically found for gas-free debris disks, and they are the same as actually retrieved from observations of many debris disks. Our specific results for three young (10-30 Myr old), spatially resolved, edge-on debris disks (beta Pic, HD 32297, and AU Mic) show that the observed radial profiles of the surface brightness do not pose any stringent constraints on the gas component of the disk. We cannot exclude that outer parts of the systems may have retained substantial amounts of primordial gas which is not evident in the gas observations (e.g. as much as 50 Earth masses for beta Pic). However, the possibility that gas, most likely secondary, is only present in little to moderate amounts, as deduced from gas detections (e.g. ~0.05 Earth masses in the beta Pic disk), remains open, too.

We report on early results of a numerical and statistical study of binary black hole inspirals. The two black holes are evolved using post-Newtonian approximations starting with initially randomly distributed spin vectors. We characterize certain aspects of the distribution shortly before merger. In particular we note the uniform distribution of black hole spin vector dot products shortly before merger and a high correlation between the initial and final black hole spin vector dot products in the equal-mass, maximally spinning case. These simulations were performed on Graphics Processing Units, and we demonstrate a speed-up of a factor 50 over a more conventional CPU implementation.

The spin of the final black hole in the coalescence of nonspinning black holes is determined by the ``residual'' orbital angular momentum of the binary. This residual momentum consists of the orbital angular momentum that the binary is not able to shed in the process of merging. We study the angular momentum radiated, the spin of the final black hole and the gravitational bursts in a series of orbits ranging from almost direct infall to numerous orbits before infall that exhibit multiple bursts of radiation in the merger process. We show that the final black hole gets a maximum spin parameter $a/M_h \le 0.78$, and this maximum occurs for initial orbital angular momentum $L \approx M^2_h$.

Generic inspirals and mergers of binary black holes produce beamed emission of gravitational radiation that can lead to a gravitational recoil or kick of the final black hole. The kick velocity depends on the mass ratio and spins of the binary as well as on the dynamics of the binary configuration. Studies have focused so far on the most astrophysically relevant configuration of quasi-circular inspirals, for which kicks as large as 3,300 km/s have been found. We present the first study of gravitational recoil in hyperbolic encounters. Contrary to quasi-circular configurations, in which the beamed radiation tends to average during the inspiral, radiation from hyperbolic encounters is plunge dominated, resulting in an enhancement of preferential beaming. As a consequence, it is possible to achieve kick velocities as large as 10,000 km/s.

We present results from numerical relativity simulations of equal mass, non-spinning binary black hole inspirals and mergers with initial eccentricities e <= 0.8 and coordinate separations D >= 12 M of up to 9 orbits (18 gravitational wave cycles). We extract the mass M_f and spin a_f of the final black hole and find, for eccentricities e < 0.4, that a_f/M_f = 0.69 and M_f/M_adm = 0.96 are independent of the initial eccentricity, suggesting that the binary has circularized by the merger time. For e > 0.5, the black holes plunge rather than orbit, and we obtain a maximum spin parameter a_f/M_f = 0.72 around e = 0.5.

Although known for almost a century, the photophoretic force has only recently been considered in astrophysical context for the first time. In our work, we have examined the effect of photophoresis, acting together with stellar gravity, radiation pressure, and gas drag, on the evolution of solids in transitional circumstellar disks. We have applied our calculations to four different systems: the disks of HR 4796A and HD 141569A, which are several Myr old AB-type stars, and two hypothetical systems that correspond to the solar nebula after disk dispersal has progressed sufficiently for the disk to become optically thin. Our results suggest that solid objects migrate inward or outward, until they reach a certain size-dependent stability distance from the star. The larger the bodies, the closer to the star they tend to accumulate. Photophoresis increases the stability radii, moving objects to larger distances. What is more, photophoresis may cause formation of a belt of objects, but only in a certain range of sizes and only around low-luminosity stars. The effects of photophoresis are noticeable in the size range from several micrometers to several centimeters (for older transitional disks) or even several meters (for younger, more gaseous, ones). We argue that due to gas damping, rotation does not substantially inhibit photophoresis.

We present a study of spinning black hole binaries focusing on the spin dynamics of the individual black holes as well as on the gravitational recoil acquired by the black hole produced by the merger. We consider two series of initial spin orientations away from the binary orbital plane. In one of the series, the spins are anti-aligned; for the second series, one of the spins points away from the binary along the line separating the black holes. We find a remarkable agreement between the spin dynamics predicted at 2nd post-Newtonian order and those from numerical relativity. For each configuration, we compute the kick of the final black hole. We use the kick estimates from the series with anti-aligned spins to fit the parameters in the \KKF{,} and verify that the recoil along the direction of the orbital angular momentum is $\propto \sin\theta$ and on the orbital plane $\propto \cos\theta$, with $\theta$ the angle between the spin directions and the orbital angular momentum. We also find that the black hole spins can be well estimated by evaluating the isolated horizon spin on spheres of constant coordinate radius.

We present results from fully nonlinear simulations of unequal mass binary black holes plunging from close separations well inside the innermost stable circular orbit with mass ratios q = M_1/M_2 = {1,0.85,0.78,0.55,0.32}, or equivalently, with reduced mass parameters $\eta=M_1M_2/(M_1+M_2)^2 = {0.25, 0.248, 0.246, 0.229, 0.183}$. For each case, the initial binary orbital parameters are chosen from the Cook-Baumgarte equal-mass ISCO configuration. We show waveforms of the dominant l=2,3 modes and compute estimates of energy and angular momentum radiated. For the plunges from the close separations considered, we measure kick velocities from gravitational radiation recoil in the range 25-82 km/s. Due to the initial close separations our kick velocity estimates should be understood as a lower bound. The close configurations considered are also likely to contain significant eccentricities influencing the recoil velocity.

Tremendous progress has been made towards the solution of the binary-black-hole problem in numerical relativity. The waveforms produced by numerical relativity will play a role in gravitational wave detection as either test-beds for analytic template banks or as template banks themselves. As the parameter space explored by numerical relativity expands, the importance of quantifying the effect that each parameter has on first the detection of gravitational waves and then the parameter estimation of their sources increases. In light of this, we present a study of equal-mass, spinning binary-black-hole evolutions through matched filtering techniques commonly used in data analysis. We study how the match between two numerical waveforms varies with numerical resolution, initial angular momentum of the black holes and the inclination angle between the source and the detector. This study is limited by the fact that the spinning black-hole-binaries are oriented axially and the waveforms only contain approximately two and a half orbits before merger. We find that for detection purposes, spinning black holes require the inclusion of the higher harmonics in addition to the dominant mode, a condition that becomes more important as the black-hole-spins increase. In addition, we conduct a preliminary investigation of how well a template of fixed spin and inclination angle can detect target templates of arbitrary spin and inclination for the axial case considered here.

The inspiral and merger of binary black holes will likely involve black holes with both unequal masses and arbitrary spins. The gravitational radiation emitted by these binaries will carry angular as well as linear momentum. A net flux of emitted linear momentum implies that the black hole produced by the merger will experience a recoil or kick. Previous studies have focused on the recoil velocity from unequal mass, non-spinning binaries. We present results from simulations of equal mass but spinning black hole binaries and show how a significant gravitational recoil can also be obtained in these situations. We consider the case of black holes with opposite spins of magnitude $a$ aligned/anti-aligned with the orbital angular momentum, with $a$ the dimensionless spin parameters of the individual holes. For the initial setups under consideration, we find a recoil velocity of $V = 475 \KMS a$. Supermassive black hole mergers producing kicks of this magnitude could result in the ejection from the cores of dwarf galaxies of the final hole produced by the collision.

We have mapped the inner 360'' regions of M51 in the 158micron [CII] line at 55'' spatial resolution using the Far-infrared Imaging Fabry-Perot Interferometer (FIFI) on the Kuiper Airborne Observatory (KAO). The emission is peaked at the nucleus, but is detectable over the entire region mapped, which covers much of the optical disk of the galaxy. There are also two strong secondary peaks at ~43% to 70% of the nuclear value located roughly 120'' to the north-east, and south-west of the nucleus. These secondary peaks are at the same distance from the nucleus as the corotation radius of the density wave pattern. The density wave also terminates at this location, and the outlying spiral structure is attributed to material clumping due to the interaction between M51 and NGC5195. This orbit crowding results in cloud-cloud collisions, stimulating star formation, that we see as enhanced [CII] line emission. The [CII] emission at the peaks originates mainly from photodissociation regions (PDRs) formed on the surfaces of molecular clouds that are exposed to OB starlight, so that these [CII] peaks trace star formation peaks in M51. The total mass of [CII] emitting photodissociated gas is ~2.6x10^{8} M_{sun}, or about 2% of the molecular gas as estimated from its CO(1-0) line emission. At the peak [CII] positions, the PDR gas mass to total gas mass fraction is somewhat higher, 3-17%, and at the secondary peaks the mass fraction of the [CII] emitting photodissociated gas can be as high as 72% of the molecular mass.... (continued)

We present observations of NGC 4038/39 in the [CII] 158 micron fine structure line taken with the MPE/UCB Far-infrared Imaging Fabry-Perot Interferometer (FIFI) on the KAO with 55" resolution. The total [CII] luminosity of the Antennae is 3.7x10^8 L(sun), which is about 1% of the FIR luminosity. The main part of the [CII] emission probably arises from photodissociation regions (PDRs). Up to 1/3 of the observed [CII] emission may originate from the CNM for higher temperatures and densities. From PDR models we derive densities and FUV intensities for the nuclei and the interaction zone. The minimum masses in the [CII] emitting regions in the interaction zone and the nuclei are a few x 10^(7) M(sun). We find that the [CII] to CO intensity ratio at the interaction zone is a factor of 2.6 lower than usually observed in starburst galaxies, but still a factor of about 1.3 to 1.4 higher than at the nuclei of NGC 4038/39, indicating no global starburst is taking place in the Antennae. [CII] emission arising partly from confined starburst regions and partly from surrounding quiescent clouds could explain the observed [CII] radiation at the interaction zone and the nuclei. There are small confined regions with high star formation activity in the interaction zone and with a lower star formation activity in the nuclei. This supports the high density and high FUV intensity for the PDRs in the interaction zone and the nuclei.