-
In the case of ground-based telescopes equipped with adaptive optics systems,
the point spread function (PSF) is only poorly known or completely unknown.
Moreover, an accurate modeling of the PSF is in general not available.
Therefore in several imaging situations the so-called blind deconvolution
methods, aiming at estimating both the scientific target and the PSF from the
detected image, can be useful. A blind deconvolution problem is severely
ill-posed and, in order to reduce the extremely large number of possible
solutions, it is necessary to introduce sensible constraints on both the
scientific target and the PSF. In a previous paper we proposed a sound
mathematical approach based on a suitable inexact alternating minimization
strategy for minimizing the generalized Kullback-Leibler divergence, assuring
global convergence. In the framework of this method we showed that an important
constraint on the PSF is the upper bound which can be derived from the
knowledge of its Strehl ratio. The efficacy of the approach was demonstrated by
means of numerical simulations. In this paper, besides improving the previous
approach by the use of a further constraint on the unknown scientific target,
we extend it to the case of multiple images of the same target obtained with
different PSFs. The main application we have in mind is to Fizeau
interferometry. As it is known this is a special feature of the Large Binocular
Telescope (LBT). The method is applied to realistic simulations of imaging both
by single mirrors and Fizeau interferometers. Successes and failures of the
method in the imaging of stellar fields are demonstrated in simple cases. These
preliminary results look promising at least in specific situations. The IDL
code of the proposed method is available on request and will be included in the
forthcoming version of the Software Package AIRY (v.6.1).
-
In this paper we present new optimization strategies for the reconstruction
of X-ray images of solar flares by means of the data collected by the Reuven
Ramaty High Energy Solar Spectroscopic Imager (RHESSI). The imaging concept of
the satellite is based of rotating modulation collimator instruments, which
allow the use of both Fourier imaging approaches and reconstruction techniques
based on the straightforward inversion of the modulated count profiles.
Although in the last decade a greater attention has been devoted to the former
strategies due to their very limited computational cost, here we consider the
latter model and investigate the effectiveness of different accelerated
gradient methods for the solution of the corresponding constrained minimization
problem. Moreover, regularization is introduced through either an early
stopping of the iterative procedure, or a Tikhonov term added to the
discrepancy function, by means of a discrepancy principle accounting for the
Poisson nature of the noise affecting the data.
-
Compact binary systems with neutron stars or black holes are one of the most
promising sources for ground-based gravitational wave detectors. Gravitational
radiation encodes rich information about source physics; thus parameter
estimation and model selection are crucial analysis steps for any detection
candidate events. Detailed models of the anticipated waveforms enable inference
on several parameters, such as component masses, spins, sky location and
distance that are essential for new astrophysical studies of these sources.
However, accurate measurements of these parameters and discrimination of models
describing the underlying physics are complicated by artifacts in the data,
uncertainties in the waveform models and in the calibration of the detectors.
Here we report such measurements on a selection of simulated signals added
either in hardware or software to the data collected by the two LIGO
instruments and the Virgo detector during their most recent joint science run,
including a "blind injection" where the signal was not initially revealed to
the collaboration. We exemplify the ability to extract information about the
source physics on signals that cover the neutron star and black hole parameter
space over the individual mass range 1 Msun - 25 Msun and the full range of
spin parameters. The cases reported in this study provide a snap-shot of the
status of parameter estimation in preparation for the operation of advanced
detectors.
-
In this paper we propose a blind deconvolution method which applies to data
perturbed by Poisson noise. The objective function is a generalized
Kullback-Leibler divergence, depending on both the unknown object and unknown
point spread function (PSF), without the addition of regularization terms;
constrained minimization, with suitable convex constraints on both unknowns, is
considered. The problem is nonconvex and we propose to solve it by means of an
inexact alternating minimization method, whose global convergence to stationary
points of the objective function has been recently proved in a general setting.
The method is iterative and each iteration, also called outer iteration,
consists of alternating an update of the object and the PSF by means of fixed
numbers of iterations, also called inner iterations, of the scaled gradient
projection (SGP) method. The use of SGP has two advantages: first, it allows to
prove global convergence of the blind method; secondly, it allows the
introduction of different constraints on the object and the PSF. The specific
constraint on the PSF, besides non-negativity and normalization, is an upper
bound derived from the so-called Strehl ratio, which is the ratio between the
peak value of an aberrated versus a perfect wavefront. Therefore a typical
application is the imaging of modern telescopes equipped with adaptive optics
systems for partial correction of the aberrations due to atmospheric
turbulence. In the paper we describe the algorithm and we recall the results
leading to its convergence. Moreover we illustrate its effectiveness by means
of numerical experiments whose results indicate that the method, pushed to
convergence, is very promising in the reconstruction of non-dense stellar
clusters. The case of more complex astronomical targets is also considered, but
in this case regularization by early stopping of the outer iterations is
required.
-
We report a search for gravitational waves from the inspiral, merger and
ringdown of binary black holes (BBH) with total mass between 25 and 100 solar
masses, in data taken at the LIGO and Virgo observatories between July 7, 2009
and October 20, 2010. The maximum sensitive distance of the detectors over this
period for a (20,20) Msun coalescence was 300 Mpc. No gravitational wave
signals were found. We thus report upper limits on the astrophysical
coalescence rates of BBH as a function of the component masses for non-spinning
components, and also evaluate the dependence of the search sensitivity on
component spins aligned with the orbital angular momentum. We find an upper
limit at 90% confidence on the coalescence rate of BBH with non-spinning
components of mass between 19 and 28 Msun of 3.3 \times 10^-7 mergers /Mpc^3
/yr.
-
We present the results of the first search for gravitational wave bursts
associated with high energy neutrinos. Together, these messengers could reveal
new, hidden sources that are not observed by conventional photon astronomy,
particularly at high energy. Our search uses neutrinos detected by the
underwater neutrino telescope ANTARES in its 5 line configuration during the
period January - September 2007, which coincided with the fifth and first
science runs of LIGO and Virgo, respectively. The LIGO-Virgo data were analysed
for candidate gravitational-wave signals coincident in time and direction with
the neutrino events. No significant coincident events were observed. We place
limits on the density of joint high energy neutrino - gravitational wave
emission events in the local universe, and compare them with densities of
merger and core-collapse events.
-
We present the first multi-wavelength follow-up observations of two candidate
gravitational-wave (GW) transient events recorded by LIGO and Virgo in their
2009-2010 science run. The events were selected with low latency by the network
of GW detectors and their candidate sky locations were observed by the Swift
observatory. Image transient detection was used to analyze the collected
electromagnetic data, which were found to be consistent with background.
Off-line analysis of the GW data alone has also established that the selected
GW events show no evidence of an astrophysical origin; one of them is
consistent with background and the other one was a test, part of a "blind
injection challenge". With this work we demonstrate the feasibility of rapid
follow-ups of GW transients and establish the sensitivity improvement joint
electromagnetic and GW observations could bring. This is a first step toward an
electromagnetic follow-up program in the regime of routine detections with the
advanced GW instruments expected within this decade. In that regime
multi-wavelength observations will play a significant role in completing the
astrophysical identification of GW sources. We present the methods and results
from this first combined analysis and discuss its implications in terms of
sensitivity for the present and future instruments.
-
The Richardson-Lucy method is the most popular deconvolution method in
astronomy because it preserves the number of counts and the non-negativity of
the original object. Regularization is, in general, obtained by an early
stopping of Richardson-Lucy iterations. In the case of point-wise objects such
as binaries or open star clusters, iterations can be pushed to convergence.
However, it is well-known that Richardson-Lucy is an inefficient method. In
most cases, acceptable solutions are obtained at the cost of hundreds or
thousands of iterations. A general optimization method, referred to as the
scaled gradient projection method, has been proposed for the constrained
minimization of continuously differentiable convex functions. It is applicable
to the non-negative minimization of the Kullback-Leibler divergence. If the
scaling suggested by Richardson-Lucy is used in this method, then it provides a
considerable increase in the efficiency of Richardson-Lucy. Therefore the aim
of this paper is to apply the scaled gradient projection method to a number of
imaging problems in astronomy such as single image deconvolution, multiple
image deconvolution, and boundary effect correction. The corresponding
algorithms are derived and implemented in interactive data language. To attempt
to achieve a further increase in efficiency, we also consider an implementation
on graphic processing units. The proposed algorithms are tested on simulated
images. The acceleration of scaled gradient projection methods achieved with
respect to the corresponding Richardson-Lucy methods strongly depends on both
the problem and the specific object to be reconstructed, and in our simulations
the improvement achieved ranges from about a factor of 4 to more than 30.
Moreover, significant accelerations of up to two orders of magnitude have been
observed between the serial and parallel implementations of the algorithms.
-
We present the results of a search for gravitational waves associated with
154 gamma-ray bursts (GRBs) that were detected by satellite-based gamma-ray
experiments in 2009-2010, during the sixth LIGO science run and the second and
third Virgo science runs. We perform two distinct searches: a modeled search
for coalescences of either two neutron stars or a neutron star and black hole;
and a search for generic, unmodeled gravitational-wave bursts. We find no
evidence for gravitational-wave counterparts, either with any individual GRB in
this sample or with the population as a whole. For all GRBs we place lower
bounds on the distance to the progenitor, under the optimistic assumption of a
gravitational-wave emission energy of 10^-2 M c^2 at 150 Hz, with a median
limit of 17 Mpc. For short hard GRBs we place exclusion distances on binary
neutron star and neutron star-black hole progenitors, using astrophysically
motivated priors on the source parameters, with median values of 16 Mpc and 28
Mpc respectively. These distance limits, while significantly larger than for a
search that is not aided by GRB satellite observations, are not large enough to
expect a coincidence with a GRB. However, projecting these exclusions to the
sensitivities of Advanced LIGO and Virgo, which should begin operation in 2015,
we find that the detection of gravitational waves associated with GRBs will
become quite possible.
-
This paper presents results of an all-sky searches for periodic gravitational
waves in the frequency range [50, 1190] Hz and with frequency derivative ranges
of [-2 x 10^-9, 1.1 x 10^-10] Hz/s for the fifth LIGO science run (S5). The
novelty of the search lies in the use of a non-coherent technique based on the
Hough-transform to combine the information from coherent searches on timescales
of about one day. Because these searches are very computationally intensive,
they have been deployed on the Einstein@Home distributed computing project
infrastructure. The search presented here is about a factor 3 more sensitive
than the previous Einstein@Home search in early S5 LIGO data. The
post-processing has left us with eight surviving candidates. We show that
deeper follow-up studies rule each of them out. Hence, since no statistically
significant gravitational wave signals have been detected, we report upper
limits on the intrinsic gravitational wave amplitude h0. For example, in the
0.5 Hz-wide band at 152.5 Hz, we can exclude the presence of signals with h0
greater than 7.6 x 10^-25 with a 90% confidence level.
-
Between 2007 and 2010 Virgo collected data in coincidence with the LIGO and
GEO gravitational-wave (GW) detectors. These data have been searched for GWs
emitted by cataclysmic phenomena in the universe, by non-axisymmetric rotating
neutron stars or from a stochastic background in the frequency band of the
detectors. The sensitivity of GW searches is limited by noise produced by the
detector or its environment. It is therefore crucial to characterize the
various noise sources in a GW detector. This paper reviews the Virgo detector
noise sources, noise propagation, and conversion mechanisms which were
identified in the three first Virgo observing runs. In many cases, these
investigations allowed us to mitigate noise sources in the detector, or to
selectively flag noise events and discard them from the data. We present
examples from the joint LIGO-GEO-Virgo GW searches to show how well noise
transients and narrow spectral lines have been identified and excluded from the
Virgo data. We also discuss how detector characterization can improve the
astrophysical reach of gravitational-wave searches.
-
The advanced interferometer network will herald a new era in observational
astronomy. There is a very strong science case to go beyond the advanced
detector network and build detectors that operate in a frequency range from 1
Hz-10 kHz, with sensitivity a factor ten better in amplitude. Such detectors
will be able to probe a range of topics in nuclear physics, astronomy,
cosmology and fundamental physics, providing insights into many unsolved
problems in these areas.
-
We present the results of a weakly modeled burst search for gravitational
waves from mergers of non-spinning intermediate mass black holes (IMBH) in the
total mass range 100--450 solar masses and with the component mass ratios
between 1:1 and 4:1. The search was conducted on data collected by the LIGO and
Virgo detectors between November of 2005 and October of 2007. No plausible
signals were observed by the search which constrains the astrophysical rates of
the IMBH mergers as a function of the component masses. In the most efficiently
detected bin centered on 88+88 solar masses, for non-spinning sources, the rate
density upper limit is 0.13 per Mpc^3 per Myr at the 90% confidence level.
-
We present results from a search for gravitational-wave bursts in the data
collected by the LIGO and Virgo detectors between July 7, 2009 and October 20,
2010: data are analyzed when at least two of the three LIGO-Virgo detectors are
in coincident operation, with a total observation time of 207 days. The
analysis searches for transients of duration < 1 s over the frequency band
64-5000 Hz, without other assumptions on the signal waveform, polarization,
direction or occurrence time. All identified events are consistent with the
expected accidental background. We set frequentist upper limits on the rate of
gravitational-wave bursts by combining this search with the previous LIGO-Virgo
search on the data collected between November 2005 and October 2007. The upper
limit on the rate of strong gravitational-wave bursts at the Earth is 1.3
events per year at 90% confidence. We also present upper limits on source rate
density per year and Mpc^3 for sample populations of standard-candle sources.
As in the previous joint run, typical sensitivities of the search in terms of
the root-sum-squared strain amplitude for these waveforms lie in the range 5
10^-22 Hz^-1/2 to 1 10^-20 Hz^-1/2. The combination of the two joint runs
entails the most sensitive all-sky search for generic gravitational-wave bursts
and synthesizes the results achieved by the initial generation of
interferometric detectors.
-
A stochastic background of gravitational waves is expected to arise from a
superposition of many incoherent sources of gravitational waves, of either
cosmological or astrophysical origin. This background is a target for the
current generation of ground-based detectors. In this article we present the
first joint search for a stochastic background using data from the LIGO and
Virgo interferometers. In a frequency band of 600-1000 Hz, we obtained a 95%
upper limit on the amplitude of $\Omega_{\rm GW}(f) = \Omega_3 (f/900
\mathrm{Hz})^3$, of $\Omega_3 < 0.33$, assuming a value of the Hubble parameter
of $h_{100}=0.72$. These new limits are a factor of seven better than the
previous best in this frequency band.
-
We report on a search for gravitational waves from coalescing compact
binaries using LIGO and Virgo observations between July 7, 2009 and October 20,
2010. We searched for signals from binaries with total mass between 2 and 25
solar masses; this includes binary neutron stars, binary black holes, and
binaries consisting of a black hole and neutron star. The detectors were
sensitive to systems up to 40 Mpc distant for binary neutron stars, and further
for higher mass systems. No gravitational-wave signals were detected. We report
upper limits on the rate of compact binary coalescence as a function of total
mass, including the results from previous LIGO and Virgo observations. The
cumulative 90%-confidence rate upper limits of the binary coalescence of binary
neutron star, neutron star- black hole and binary black hole systems are 1.3 x
10^{-4}, 3.1 x 10^{-5} and 6.4 x 10^{-6} Mpc^{-3}yr^{-1}, respectively. These
upper limits are up to a factor 1.4 lower than previously derived limits. We
also report on results from a blind injection challenge.
-
Aims. A transient astrophysical event observed in both gravitational wave
(GW) and electromagnetic (EM) channels would yield rich scientific rewards. A
first program initiating EM follow-ups to possible transient GW events has been
developed and exercised by the LIGO and Virgo community in association with
several partners. In this paper, we describe and evaluate the methods used to
promptly identify and localize GW event candidates and to request images of
targeted sky locations.
Methods. During two observing periods (Dec 17 2009 to Jan 8 2010 and Sep 2 to
Oct 20 2010), a low-latency analysis pipeline was used to identify GW event
candidates and to reconstruct maps of possible sky locations. A catalog of
nearby galaxies and Milky Way globular clusters was used to select the most
promising sky positions to be imaged, and this directional information was
delivered to EM observatories with time lags of about thirty minutes. A Monte
Carlo simulation has been used to evaluate the low-latency GW pipeline's
ability to reconstruct source positions correctly.
Results. For signals near the detection threshold, our low-latency algorithms
often localized simulated GW burst signals to tens of square degrees, while
neutron star/neutron star inspirals and neutron star/black hole inspirals were
localized to a few hundred square degrees. Localization precision improves for
moderately stronger signals. The correct sky location of signals well above
threshold and originating from nearby galaxies may be observed with ~50% or
better probability with a few pointings of wide-field telescopes.
-
X-radiation from energetic electrons is the prime diagnostic of
flare-accelerated electrons. The observed X-ray flux (and polarization state)
is fundamentally a convolution of the cross-section for the hard X-ray emission
process(es) in question with the electron distribution function, which is in
turn a function of energy, direction, spatial location and time. To address the
problems of particle propagation and acceleration one needs to infer as much
information as possible on this electron distribution function, through a
deconvolution of this fundamental relationship. This review presents recent
progress toward this goal using spectroscopic, imaging and polarization
measurements, primarily from the Reuven Ramaty High Energy Solar Spectroscopic
Imager (RHESSI). Previous conclusions regarding the energy, angular (pitch
angle) and spatial distributions of energetic electrons in solar flares are
critically reviewed. We discuss the role and the observational evidence of
several radiation processes: free-free electron-ion, free-free
electron-electron, free-bound electron-ion bremsstrahlung, photoelectric
absorption and Compton back-scatter (albedo), using both spectroscopic and
imaging techniques. This unprecedented quality of data allows for the first
time inference of the angular distributions of the X-ray-emitting electrons
using albedo, improved model-independent inference of electron energy spectra
and emission measures of thermal plasma. Moreover, imaging spectroscopy has
revealed hitherto unknown details of solar flare morphology and detailed
spectroscopy of coronal, footpoint and extended sources in flaring regions.
Additional attempts to measure hard X-ray polarization were not sufficient to
put constraints on the degree of anisotropy of electrons, but point to the
importance of obtaining good quality polarization data.
-
We report on an all-sky search for periodic gravitational waves in the
frequency band 50-800 Hz and with the frequency time derivative in the range of
0 through -6e-9 Hz/s. Such a signal could be produced by a nearby spinning and
slightly non-axisymmetric isolated neutron star in our galaxy. After recent
improvements in the search program that yielded a 10x increase in computational
efficiency, we have searched in two years of data collected during LIGO's fifth
science run and have obtained the most sensitive all-sky upper limits on
gravitational wave strain to date. Near 150 Hz our upper limit on worst-case
linearly polarized strain amplitude $h_0$ is 1e-24, while at the high end of
our frequency range we achieve a worst-case upper limit of 3.8e-24 for all
polarizations and sky locations. These results constitute a factor of two
improvement upon previously published data. A new detection pipeline utilizing
a Loosely Coherent algorithm was able to follow up weaker outliers, increasing
the volume of space where signals can be detected by a factor of 10, but has
not revealed any gravitational wave signals. The pipeline has been tested for
robustness with respect to deviations from the model of an isolated neutron
star, such as caused by a low-mass or long-period binary companion.
-
The gravitational-wave (GW) sky may include nearby pointlike sources as well
as astrophysical and cosmological stochastic backgrounds. Since the relative
strength and angular distribution of the many possible sources of GWs are not
well constrained, searches for GW signals must be performed in a
model-independent way. To that end we perform two directional searches for
persistent GWs using data from the LIGO S5 science run: one optimized for
pointlike sources and one for arbitrary extended sources. The latter result is
the first of its kind. Finding no evidence to support the detection of GWs, we
present 90% confidence level (CL) upper-limit maps of GW strain power with
typical values between 2-20x10^-50 strain^2 Hz^-1 and 5-35x10^-49 strain^2
Hz^-1 sr^-1 for pointlike and extended sources respectively. The limits on
pointlike sources constitute a factor of 30 improvement over the previous best
limits. We also set 90% CL limits on the narrow-band root-mean-square GW strain
from interesting targets including Sco X-1, SN1987A and the Galactic Center as
low as ~7x10^-25 in the most sensitive frequency range near 160 Hz. These
limits are the most constraining to date and constitute a factor of 5
improvement over the previous best limits.
-
The Virgo gravitational wave detector is an interferometer (ITF) with 3km
arms located in Pisa, Italy. From July to October 2010, Virgo performed its
third science run (VSR3) in coincidence with the LIGO detectors. Despite
several techniques adopted to isolate the interferometer from the environment,
seismic noise remains an important issue for Virgo. Vibrations produced by the
detector infrastructure (such as air conditioning units, water
chillers/heaters, pumps) are found to affect Virgo's sensitivity, with the main
coupling mechanisms being through beam jitter and scattered light processes.
The Advanced Virgo (AdV) design seeks to reduce ITF couplings to environmental
noise by having most vibration-sensitive components suspended and in-vacuum, as
well as muffle and relocate loud machines. During the months of June and July
2010, a Guralp-3TD seismometer was stationed at various locations around the
Virgo site hosting major infrastructure machines. Seismic data were examined
using spectral and coherence analysis with seismic probes close to the
detector. The primary aim of this study was to identify noisy machines which
seismically affect the ITF environment and thus require mitigation attention.
Analyzed machines are located at various distances from the experimental halls,
ranging from 10m to 100m. An attempt is made to measure the attenuation of
emitted noise at the ITF and correlate it to the distance from the source and
to seismic attenuation models in soil.
-
We present direct upper limits on continuous gravitational wave emission from
the Vela pulsar using data from the Virgo detector's second science run. These
upper limits have been obtained using three independent methods that assume the
gravitational wave emission follows the radio timing. Two of the methods
produce frequentist upper limits for an assumed known orientation of the star's
spin axis and value of the wave polarization angle of, respectively,
$1.9\ee{-24}$ and $2.2\ee{-24}$, with 95% confidence. The third method, under
the same hypothesis, produces a Bayesian upper limit of $2.1\ee{-24}$, with 95%
degree of belief. These limits are below the indirect {\it spin-down limit} of
$3.3\ee{-24}$ for the Vela pulsar, defined by the energy loss rate inferred
from observed decrease in Vela's spin frequency, and correspond to a limit on
the star ellipticity of $\sim 10^{-3}$. Slightly less stringent results, but
still well below the spin-down limit, are obtained assuming the star's spin
axis inclination and the wave polarization angles are unknown.
-
Soft gamma repeaters (SGRs) and anomalous X-ray pulsars (AXPs) are thought to
be magnetars: neutron stars powered by extreme magnetic fields. These rare
objects are characterized by repeated and sometimes spectacular gamma-ray
bursts. The burst mechanism might involve crustal fractures and excitation of
non-radial modes which would emit gravitational waves (GWs). We present the
results of a search for GW bursts from six galactic magnetars that is sensitive
to neutron star f-modes, thought to be the most efficient GW emitting
oscillatory modes in compact stars. One of them, SGR 0501+4516, is likely ~1
kpc from Earth, an order of magnitude closer than magnetars targeted in
previous GW searches. A second, AXP 1E 1547.0-5408, gave a burst with an
estimated isotropic energy >10^{44} erg which is comparable to the giant
flares. We find no evidence of GWs associated with a sample of 1279
electromagnetic triggers from six magnetars occurring between November 2006 and
June 2009, in GW data from the LIGO, Virgo, and GEO600 detectors. Our lowest
model-dependent GW emission energy upper limits for band- and time-limited
white noise bursts in the detector sensitive band, and for f-mode ringdowns (at
1090 Hz), are 3.0x10^{44} d_1^2 erg and 1.4x10^{47} d_1^2 erg respectively,
where d_1 = d_{0501} / 1 kpc and d_{0501} is the distance to SGR 0501+4516.
These limits on GW emission from f-modes are an order of magnitude lower than
any previous, and approach the range of electromagnetic energies seen in SGR
giant flares for the first time.
-
We present the first modeled search for gravitational waves using the
complete binary black hole gravitational waveform from inspiral through the
merger and ringdown for binaries with negligible component spin. We searched
approximately 2 years of LIGO data taken between November 2005 and September
2007 for systems with component masses of 1-99 solar masses and total masses of
25-100 solar masses. We did not detect any plausible gravitational-wave signals
but we do place upper limits on the merger rate of binary black holes as a
function of the component masses in this range. We constrain the rate of
mergers for binary black hole systems with component masses between 19 and 28
solar masses and negligible spin to be no more than 2.0 per Mpc^3 per Myr at
90% confidence.
-
The Virgo detector is a kilometer-length interferometer for gravitational
wave detection located near Pisa (Italy). During its second science run (VSR2)
in 2009, six months of data were accumulated with a sensitivity close to its
design. In this paper, the methods used to determine the parameters for
sensitivity estimation and gravitational wave reconstruction are described. The
main quantities to be calibrated are the frequency response of the mirror
actuation and the sensing of the output power. Focus is also put on their
absolute timing. The monitoring of the calibration data as well as the
parameter estimation with independent techniques are discussed to provide an
estimation of the calibration uncertainties. Finally, the estimation of the
Virgo sensitivity in the frequency-domain is described and typical
sensitivities measured during VSR2 are shown.
In the case of ground-based telescopes equipped with adaptive optics systems, the point spread function (PSF) is only poorly known or completely unknown. Moreover, an accurate modeling of the PSF is in general not available. Therefore in several imaging situations the so-called blind deconvolution methods, aiming at estimating both the scientific target and the PSF from the detected image, can be useful. A blind deconvolution problem is severely ill-posed and, in order to reduce the extremely large number of possible solutions, it is necessary to introduce sensible constraints on both the scientific target and the PSF. In a previous paper we proposed a sound mathematical approach based on a suitable inexact alternating minimization strategy for minimizing the generalized Kullback-Leibler divergence, assuring global convergence. In the framework of this method we showed that an important constraint on the PSF is the upper bound which can be derived from the knowledge of its Strehl ratio. The efficacy of the approach was demonstrated by means of numerical simulations. In this paper, besides improving the previous approach by the use of a further constraint on the unknown scientific target, we extend it to the case of multiple images of the same target obtained with different PSFs. The main application we have in mind is to Fizeau interferometry. As it is known this is a special feature of the Large Binocular Telescope (LBT). The method is applied to realistic simulations of imaging both by single mirrors and Fizeau interferometers. Successes and failures of the method in the imaging of stellar fields are demonstrated in simple cases. These preliminary results look promising at least in specific situations. The IDL code of the proposed method is available on request and will be included in the forthcoming version of the Software Package AIRY (v.6.1).
In this paper we present new optimization strategies for the reconstruction of X-ray images of solar flares by means of the data collected by the Reuven Ramaty High Energy Solar Spectroscopic Imager (RHESSI). The imaging concept of the satellite is based of rotating modulation collimator instruments, which allow the use of both Fourier imaging approaches and reconstruction techniques based on the straightforward inversion of the modulated count profiles. Although in the last decade a greater attention has been devoted to the former strategies due to their very limited computational cost, here we consider the latter model and investigate the effectiveness of different accelerated gradient methods for the solution of the corresponding constrained minimization problem. Moreover, regularization is introduced through either an early stopping of the iterative procedure, or a Tikhonov term added to the discrepancy function, by means of a discrepancy principle accounting for the Poisson nature of the noise affecting the data.
Compact binary systems with neutron stars or black holes are one of the most promising sources for ground-based gravitational wave detectors. Gravitational radiation encodes rich information about source physics; thus parameter estimation and model selection are crucial analysis steps for any detection candidate events. Detailed models of the anticipated waveforms enable inference on several parameters, such as component masses, spins, sky location and distance that are essential for new astrophysical studies of these sources. However, accurate measurements of these parameters and discrimination of models describing the underlying physics are complicated by artifacts in the data, uncertainties in the waveform models and in the calibration of the detectors. Here we report such measurements on a selection of simulated signals added either in hardware or software to the data collected by the two LIGO instruments and the Virgo detector during their most recent joint science run, including a "blind injection" where the signal was not initially revealed to the collaboration. We exemplify the ability to extract information about the source physics on signals that cover the neutron star and black hole parameter space over the individual mass range 1 Msun - 25 Msun and the full range of spin parameters. The cases reported in this study provide a snap-shot of the status of parameter estimation in preparation for the operation of advanced detectors.
In this paper we propose a blind deconvolution method which applies to data perturbed by Poisson noise. The objective function is a generalized Kullback-Leibler divergence, depending on both the unknown object and unknown point spread function (PSF), without the addition of regularization terms; constrained minimization, with suitable convex constraints on both unknowns, is considered. The problem is nonconvex and we propose to solve it by means of an inexact alternating minimization method, whose global convergence to stationary points of the objective function has been recently proved in a general setting. The method is iterative and each iteration, also called outer iteration, consists of alternating an update of the object and the PSF by means of fixed numbers of iterations, also called inner iterations, of the scaled gradient projection (SGP) method. The use of SGP has two advantages: first, it allows to prove global convergence of the blind method; secondly, it allows the introduction of different constraints on the object and the PSF. The specific constraint on the PSF, besides non-negativity and normalization, is an upper bound derived from the so-called Strehl ratio, which is the ratio between the peak value of an aberrated versus a perfect wavefront. Therefore a typical application is the imaging of modern telescopes equipped with adaptive optics systems for partial correction of the aberrations due to atmospheric turbulence. In the paper we describe the algorithm and we recall the results leading to its convergence. Moreover we illustrate its effectiveness by means of numerical experiments whose results indicate that the method, pushed to convergence, is very promising in the reconstruction of non-dense stellar clusters. The case of more complex astronomical targets is also considered, but in this case regularization by early stopping of the outer iterations is required.
We report a search for gravitational waves from the inspiral, merger and ringdown of binary black holes (BBH) with total mass between 25 and 100 solar masses, in data taken at the LIGO and Virgo observatories between July 7, 2009 and October 20, 2010. The maximum sensitive distance of the detectors over this period for a (20,20) Msun coalescence was 300 Mpc. No gravitational wave signals were found. We thus report upper limits on the astrophysical coalescence rates of BBH as a function of the component masses for non-spinning components, and also evaluate the dependence of the search sensitivity on component spins aligned with the orbital angular momentum. We find an upper limit at 90% confidence on the coalescence rate of BBH with non-spinning components of mass between 19 and 28 Msun of 3.3 \times 10^-7 mergers /Mpc^3 /yr.
We present the first multi-wavelength follow-up observations of two candidate gravitational-wave (GW) transient events recorded by LIGO and Virgo in their 2009-2010 science run. The events were selected with low latency by the network of GW detectors and their candidate sky locations were observed by the Swift observatory. Image transient detection was used to analyze the collected electromagnetic data, which were found to be consistent with background. Off-line analysis of the GW data alone has also established that the selected GW events show no evidence of an astrophysical origin; one of them is consistent with background and the other one was a test, part of a "blind injection challenge". With this work we demonstrate the feasibility of rapid follow-ups of GW transients and establish the sensitivity improvement joint electromagnetic and GW observations could bring. This is a first step toward an electromagnetic follow-up program in the regime of routine detections with the advanced GW instruments expected within this decade. In that regime multi-wavelength observations will play a significant role in completing the astrophysical identification of GW sources. We present the methods and results from this first combined analysis and discuss its implications in terms of sensitivity for the present and future instruments.
The Richardson-Lucy method is the most popular deconvolution method in astronomy because it preserves the number of counts and the non-negativity of the original object. Regularization is, in general, obtained by an early stopping of Richardson-Lucy iterations. In the case of point-wise objects such as binaries or open star clusters, iterations can be pushed to convergence. However, it is well-known that Richardson-Lucy is an inefficient method. In most cases, acceptable solutions are obtained at the cost of hundreds or thousands of iterations. A general optimization method, referred to as the scaled gradient projection method, has been proposed for the constrained minimization of continuously differentiable convex functions. It is applicable to the non-negative minimization of the Kullback-Leibler divergence. If the scaling suggested by Richardson-Lucy is used in this method, then it provides a considerable increase in the efficiency of Richardson-Lucy. Therefore the aim of this paper is to apply the scaled gradient projection method to a number of imaging problems in astronomy such as single image deconvolution, multiple image deconvolution, and boundary effect correction. The corresponding algorithms are derived and implemented in interactive data language. To attempt to achieve a further increase in efficiency, we also consider an implementation on graphic processing units. The proposed algorithms are tested on simulated images. The acceleration of scaled gradient projection methods achieved with respect to the corresponding Richardson-Lucy methods strongly depends on both the problem and the specific object to be reconstructed, and in our simulations the improvement achieved ranges from about a factor of 4 to more than 30. Moreover, significant accelerations of up to two orders of magnitude have been observed between the serial and parallel implementations of the algorithms.
We present the results of a search for gravitational waves associated with 154 gamma-ray bursts (GRBs) that were detected by satellite-based gamma-ray experiments in 2009-2010, during the sixth LIGO science run and the second and third Virgo science runs. We perform two distinct searches: a modeled search for coalescences of either two neutron stars or a neutron star and black hole; and a search for generic, unmodeled gravitational-wave bursts. We find no evidence for gravitational-wave counterparts, either with any individual GRB in this sample or with the population as a whole. For all GRBs we place lower bounds on the distance to the progenitor, under the optimistic assumption of a gravitational-wave emission energy of 10^-2 M c^2 at 150 Hz, with a median limit of 17 Mpc. For short hard GRBs we place exclusion distances on binary neutron star and neutron star-black hole progenitors, using astrophysically motivated priors on the source parameters, with median values of 16 Mpc and 28 Mpc respectively. These distance limits, while significantly larger than for a search that is not aided by GRB satellite observations, are not large enough to expect a coincidence with a GRB. However, projecting these exclusions to the sensitivities of Advanced LIGO and Virgo, which should begin operation in 2015, we find that the detection of gravitational waves associated with GRBs will become quite possible.
This paper presents results of an all-sky searches for periodic gravitational waves in the frequency range [50, 1190] Hz and with frequency derivative ranges of [-2 x 10^-9, 1.1 x 10^-10] Hz/s for the fifth LIGO science run (S5). The novelty of the search lies in the use of a non-coherent technique based on the Hough-transform to combine the information from coherent searches on timescales of about one day. Because these searches are very computationally intensive, they have been deployed on the Einstein@Home distributed computing project infrastructure. The search presented here is about a factor 3 more sensitive than the previous Einstein@Home search in early S5 LIGO data. The post-processing has left us with eight surviving candidates. We show that deeper follow-up studies rule each of them out. Hence, since no statistically significant gravitational wave signals have been detected, we report upper limits on the intrinsic gravitational wave amplitude h0. For example, in the 0.5 Hz-wide band at 152.5 Hz, we can exclude the presence of signals with h0 greater than 7.6 x 10^-25 with a 90% confidence level.
Between 2007 and 2010 Virgo collected data in coincidence with the LIGO and GEO gravitational-wave (GW) detectors. These data have been searched for GWs emitted by cataclysmic phenomena in the universe, by non-axisymmetric rotating neutron stars or from a stochastic background in the frequency band of the detectors. The sensitivity of GW searches is limited by noise produced by the detector or its environment. It is therefore crucial to characterize the various noise sources in a GW detector. This paper reviews the Virgo detector noise sources, noise propagation, and conversion mechanisms which were identified in the three first Virgo observing runs. In many cases, these investigations allowed us to mitigate noise sources in the detector, or to selectively flag noise events and discard them from the data. We present examples from the joint LIGO-GEO-Virgo GW searches to show how well noise transients and narrow spectral lines have been identified and excluded from the Virgo data. We also discuss how detector characterization can improve the astrophysical reach of gravitational-wave searches.
The advanced interferometer network will herald a new era in observational astronomy. There is a very strong science case to go beyond the advanced detector network and build detectors that operate in a frequency range from 1 Hz-10 kHz, with sensitivity a factor ten better in amplitude. Such detectors will be able to probe a range of topics in nuclear physics, astronomy, cosmology and fundamental physics, providing insights into many unsolved problems in these areas.
A stochastic background of gravitational waves is expected to arise from a superposition of many incoherent sources of gravitational waves, of either cosmological or astrophysical origin. This background is a target for the current generation of ground-based detectors. In this article we present the first joint search for a stochastic background using data from the LIGO and Virgo interferometers. In a frequency band of 600-1000 Hz, we obtained a 95% upper limit on the amplitude of $\Omega_{\rm GW}(f) = \Omega_3 (f/900 \mathrm{Hz})^3$, of $\Omega_3 < 0.33$, assuming a value of the Hubble parameter of $h_{100}=0.72$. These new limits are a factor of seven better than the previous best in this frequency band.
We report on a search for gravitational waves from coalescing compact binaries using LIGO and Virgo observations between July 7, 2009 and October 20, 2010. We searched for signals from binaries with total mass between 2 and 25 solar masses; this includes binary neutron stars, binary black holes, and binaries consisting of a black hole and neutron star. The detectors were sensitive to systems up to 40 Mpc distant for binary neutron stars, and further for higher mass systems. No gravitational-wave signals were detected. We report upper limits on the rate of compact binary coalescence as a function of total mass, including the results from previous LIGO and Virgo observations. The cumulative 90%-confidence rate upper limits of the binary coalescence of binary neutron star, neutron star- black hole and binary black hole systems are 1.3 x 10^{-4}, 3.1 x 10^{-5} and 6.4 x 10^{-6} Mpc^{-3}yr^{-1}, respectively. These upper limits are up to a factor 1.4 lower than previously derived limits. We also report on results from a blind injection challenge.
Aims. A transient astrophysical event observed in both gravitational wave (GW) and electromagnetic (EM) channels would yield rich scientific rewards. A first program initiating EM follow-ups to possible transient GW events has been developed and exercised by the LIGO and Virgo community in association with several partners. In this paper, we describe and evaluate the methods used to promptly identify and localize GW event candidates and to request images of targeted sky locations. Methods. During two observing periods (Dec 17 2009 to Jan 8 2010 and Sep 2 to Oct 20 2010), a low-latency analysis pipeline was used to identify GW event candidates and to reconstruct maps of possible sky locations. A catalog of nearby galaxies and Milky Way globular clusters was used to select the most promising sky positions to be imaged, and this directional information was delivered to EM observatories with time lags of about thirty minutes. A Monte Carlo simulation has been used to evaluate the low-latency GW pipeline's ability to reconstruct source positions correctly. Results. For signals near the detection threshold, our low-latency algorithms often localized simulated GW burst signals to tens of square degrees, while neutron star/neutron star inspirals and neutron star/black hole inspirals were localized to a few hundred square degrees. Localization precision improves for moderately stronger signals. The correct sky location of signals well above threshold and originating from nearby galaxies may be observed with ~50% or better probability with a few pointings of wide-field telescopes.
X-radiation from energetic electrons is the prime diagnostic of flare-accelerated electrons. The observed X-ray flux (and polarization state) is fundamentally a convolution of the cross-section for the hard X-ray emission process(es) in question with the electron distribution function, which is in turn a function of energy, direction, spatial location and time. To address the problems of particle propagation and acceleration one needs to infer as much information as possible on this electron distribution function, through a deconvolution of this fundamental relationship. This review presents recent progress toward this goal using spectroscopic, imaging and polarization measurements, primarily from the Reuven Ramaty High Energy Solar Spectroscopic Imager (RHESSI). Previous conclusions regarding the energy, angular (pitch angle) and spatial distributions of energetic electrons in solar flares are critically reviewed. We discuss the role and the observational evidence of several radiation processes: free-free electron-ion, free-free electron-electron, free-bound electron-ion bremsstrahlung, photoelectric absorption and Compton back-scatter (albedo), using both spectroscopic and imaging techniques. This unprecedented quality of data allows for the first time inference of the angular distributions of the X-ray-emitting electrons using albedo, improved model-independent inference of electron energy spectra and emission measures of thermal plasma. Moreover, imaging spectroscopy has revealed hitherto unknown details of solar flare morphology and detailed spectroscopy of coronal, footpoint and extended sources in flaring regions. Additional attempts to measure hard X-ray polarization were not sufficient to put constraints on the degree of anisotropy of electrons, but point to the importance of obtaining good quality polarization data.
We report on an all-sky search for periodic gravitational waves in the frequency band 50-800 Hz and with the frequency time derivative in the range of 0 through -6e-9 Hz/s. Such a signal could be produced by a nearby spinning and slightly non-axisymmetric isolated neutron star in our galaxy. After recent improvements in the search program that yielded a 10x increase in computational efficiency, we have searched in two years of data collected during LIGO's fifth science run and have obtained the most sensitive all-sky upper limits on gravitational wave strain to date. Near 150 Hz our upper limit on worst-case linearly polarized strain amplitude $h_0$ is 1e-24, while at the high end of our frequency range we achieve a worst-case upper limit of 3.8e-24 for all polarizations and sky locations. These results constitute a factor of two improvement upon previously published data. A new detection pipeline utilizing a Loosely Coherent algorithm was able to follow up weaker outliers, increasing the volume of space where signals can be detected by a factor of 10, but has not revealed any gravitational wave signals. The pipeline has been tested for robustness with respect to deviations from the model of an isolated neutron star, such as caused by a low-mass or long-period binary companion.
The gravitational-wave (GW) sky may include nearby pointlike sources as well as astrophysical and cosmological stochastic backgrounds. Since the relative strength and angular distribution of the many possible sources of GWs are not well constrained, searches for GW signals must be performed in a model-independent way. To that end we perform two directional searches for persistent GWs using data from the LIGO S5 science run: one optimized for pointlike sources and one for arbitrary extended sources. The latter result is the first of its kind. Finding no evidence to support the detection of GWs, we present 90% confidence level (CL) upper-limit maps of GW strain power with typical values between 2-20x10^-50 strain^2 Hz^-1 and 5-35x10^-49 strain^2 Hz^-1 sr^-1 for pointlike and extended sources respectively. The limits on pointlike sources constitute a factor of 30 improvement over the previous best limits. We also set 90% CL limits on the narrow-band root-mean-square GW strain from interesting targets including Sco X-1, SN1987A and the Galactic Center as low as ~7x10^-25 in the most sensitive frequency range near 160 Hz. These limits are the most constraining to date and constitute a factor of 5 improvement over the previous best limits.
The Virgo gravitational wave detector is an interferometer (ITF) with 3km arms located in Pisa, Italy. From July to October 2010, Virgo performed its third science run (VSR3) in coincidence with the LIGO detectors. Despite several techniques adopted to isolate the interferometer from the environment, seismic noise remains an important issue for Virgo. Vibrations produced by the detector infrastructure (such as air conditioning units, water chillers/heaters, pumps) are found to affect Virgo's sensitivity, with the main coupling mechanisms being through beam jitter and scattered light processes. The Advanced Virgo (AdV) design seeks to reduce ITF couplings to environmental noise by having most vibration-sensitive components suspended and in-vacuum, as well as muffle and relocate loud machines. During the months of June and July 2010, a Guralp-3TD seismometer was stationed at various locations around the Virgo site hosting major infrastructure machines. Seismic data were examined using spectral and coherence analysis with seismic probes close to the detector. The primary aim of this study was to identify noisy machines which seismically affect the ITF environment and thus require mitigation attention. Analyzed machines are located at various distances from the experimental halls, ranging from 10m to 100m. An attempt is made to measure the attenuation of emitted noise at the ITF and correlate it to the distance from the source and to seismic attenuation models in soil.
We present direct upper limits on continuous gravitational wave emission from the Vela pulsar using data from the Virgo detector's second science run. These upper limits have been obtained using three independent methods that assume the gravitational wave emission follows the radio timing. Two of the methods produce frequentist upper limits for an assumed known orientation of the star's spin axis and value of the wave polarization angle of, respectively, $1.9\ee{-24}$ and $2.2\ee{-24}$, with 95% confidence. The third method, under the same hypothesis, produces a Bayesian upper limit of $2.1\ee{-24}$, with 95% degree of belief. These limits are below the indirect {\it spin-down limit} of $3.3\ee{-24}$ for the Vela pulsar, defined by the energy loss rate inferred from observed decrease in Vela's spin frequency, and correspond to a limit on the star ellipticity of $\sim 10^{-3}$. Slightly less stringent results, but still well below the spin-down limit, are obtained assuming the star's spin axis inclination and the wave polarization angles are unknown.
Soft gamma repeaters (SGRs) and anomalous X-ray pulsars (AXPs) are thought to be magnetars: neutron stars powered by extreme magnetic fields. These rare objects are characterized by repeated and sometimes spectacular gamma-ray bursts. The burst mechanism might involve crustal fractures and excitation of non-radial modes which would emit gravitational waves (GWs). We present the results of a search for GW bursts from six galactic magnetars that is sensitive to neutron star f-modes, thought to be the most efficient GW emitting oscillatory modes in compact stars. One of them, SGR 0501+4516, is likely ~1 kpc from Earth, an order of magnitude closer than magnetars targeted in previous GW searches. A second, AXP 1E 1547.0-5408, gave a burst with an estimated isotropic energy >10^{44} erg which is comparable to the giant flares. We find no evidence of GWs associated with a sample of 1279 electromagnetic triggers from six magnetars occurring between November 2006 and June 2009, in GW data from the LIGO, Virgo, and GEO600 detectors. Our lowest model-dependent GW emission energy upper limits for band- and time-limited white noise bursts in the detector sensitive band, and for f-mode ringdowns (at 1090 Hz), are 3.0x10^{44} d_1^2 erg and 1.4x10^{47} d_1^2 erg respectively, where d_1 = d_{0501} / 1 kpc and d_{0501} is the distance to SGR 0501+4516. These limits on GW emission from f-modes are an order of magnitude lower than any previous, and approach the range of electromagnetic energies seen in SGR giant flares for the first time.
We present the first modeled search for gravitational waves using the complete binary black hole gravitational waveform from inspiral through the merger and ringdown for binaries with negligible component spin. We searched approximately 2 years of LIGO data taken between November 2005 and September 2007 for systems with component masses of 1-99 solar masses and total masses of 25-100 solar masses. We did not detect any plausible gravitational-wave signals but we do place upper limits on the merger rate of binary black holes as a function of the component masses in this range. We constrain the rate of mergers for binary black hole systems with component masses between 19 and 28 solar masses and negligible spin to be no more than 2.0 per Mpc^3 per Myr at 90% confidence.
The Virgo detector is a kilometer-length interferometer for gravitational wave detection located near Pisa (Italy). During its second science run (VSR2) in 2009, six months of data were accumulated with a sensitivity close to its design. In this paper, the methods used to determine the parameters for sensitivity estimation and gravitational wave reconstruction are described. The main quantities to be calibrated are the frequency response of the mirror actuation and the sensing of the output power. Focus is also put on their absolute timing. The monitoring of the calibration data as well as the parameter estimation with independent techniques are discussed to provide an estimation of the calibration uncertainties. Finally, the estimation of the Virgo sensitivity in the frequency-domain is described and typical sensitivities measured during VSR2 are shown.