-
This Chapter, "High-dimensional ABC", is to appear in the forthcoming
Handbook of Approximate Bayesian Computation (2018). It details the main ideas
and concepts behind extending ABC methods to higher dimensions, with supporting
examples and illustrations.
-
This Chapter, "Overview of Approximate Bayesian Computation", is to appear as
the first chapter in the forthcoming Handbook of Approximate Bayesian
Computation (2018). It details the main ideas and concepts behind ABC methods
with many examples and illustrations.
-
This Chapter, "ABC Samplers", is to appear in the forthcoming Handbook of
Approximate Bayesian Computation (2018). It details the main ideas and
algorithms used to sample from the ABC approximation to the posterior
distribution, including methods based on rejection/importance sampling, MCMC
and sequential Monte Carlo.
-
In the context of engineered quantum systems, there is a demand for
superconducting tunable devices able to operate with high Q-factors at power
levels equivalent to only a few photons. In this work, we developed a 3D
microwave reentrant cavity with such characteristics ready to provide a very
fine-tuning of a high-Q resonant mode over a large dynamic range. This system
has an electronic tuning mechanism based on a mechanically amplified
piezoelectric actuator, which controls the resonator dominant mode frequency by
changing the cavity narrow gap by very small displacements. Experiments were
conducted at room and dilution refrigerator temperatures showing a large
dynamic range up to 4 GHz and 1 GHz, respectively, and were compared to a FEM
model simulated data. At elevated microwave power input, nonlinear thermal
effects were observed to destroy the superconductivity of the cavity due to the
large electric fields generated in the small gap of the reentrant cavity.
-
We consider the application of approximate Bayesian Computation (ABC) in the
context of medical imaging data. We consider the parameter estimation of
compartmental models in PET imaging analysis, and provide a simple ABC
algorithm for its estimation. We demonstrate the utility of the proposed
estimation methods on a neurotransmitter response model, and compare our
approach to existing methods.
-
This Supplement provides supporting material for arXiv:1602.08492 . We
briefly summarize past electromagnetic (EM) follow-up efforts as well as the
organization and policy of the current EM follow-up program. We compare the
four probability sky maps produced for the gravitational-wave transient
GW150914, and provide additional details of the EM follow-up observations that
were performed in the different bands.
-
A gravitational-wave (GW) transient was identified in data recorded by the
Advanced Laser Interferometer Gravitational-wave Observatory (LIGO) detectors
on 2015 September 14. The event, initially designated G184098 and later given
the name GW150914, is described in detail elsewhere. By prior arrangement,
preliminary estimates of the time, significance, and sky location of the event
were shared with 63 teams of observers covering radio, optical, near-infrared,
X-ray, and gamma-ray wavelengths with ground- and space-based facilities. In
this Letter we describe the low-latency analysis of the GW data and present the
sky localization of the first observed compact binary merger. We summarize the
follow-up observations reported by 25 teams via private Gamma-ray Coordinates
Network circulars, giving an overview of the participating facilities, the GW
sky localization coverage, the timeline and depth of the observations. As this
event turned out to be a binary black hole merger, there is little expectation
of a detectable electromagnetic (EM) signature. Nevertheless, this first
broadband campaign to search for a counterpart of an Advanced LIGO source
represents a milestone and highlights the broad capabilities of the transient
astronomy community and the observing strategies that have been developed to
pursue neutron star binary merger events. Detailed investigations of the EM
data and results of the EM follow-up campaign are being disseminated in papers
by the individual teams.
-
The Coronal Multichannel Polarimeter (CoMP) routinely performs coronal
polarimetric measurements using the Fe XIII 10747 $\AA$ and 10798 $\AA$ lines,
which are sensitive to the coronal magnetic field. However, inverting such
polarimetric measurements into magnetic field data is a difficult task because
the corona is optically thin at these wavelengths and the observed signal is
therefore the integrated emission of all the plasma along the line of sight. To
overcome this difficulty, we take on a new approach that combines a
parameterized 3D magnetic field model with forward modeling of the polarization
signal. For that purpose, we develop a new, fast and efficient, optimization
method for model-data fitting: the Radial-basis-functions Optimization
Approximation Method (ROAM). Model-data fitting is achieved by optimizing a
user-specified log-likelihood function that quantifies the differences between
the observed polarization signal and its synthetic/predicted analogue. Speed
and efficiency are obtained by combining sparse evaluation of the magnetic
model with radial-basis-function (RBF) decomposition of the log-likelihood
function. The RBF decomposition provides an analytical expression for the
log-likelihood function that is used to inexpensively estimate the set of
parameter values optimizing it. We test and validate ROAM on a synthetic test
bed of a coronal magnetic flux rope and show that it performs well with a
significantly sparse sample of the parameter space. We conclude that our
optimization method is well-suited for fast and efficient model-data fitting
and can be exploited for converting coronal polarimetric measurements, such as
the ones provided by CoMP, into coronal magnetic field data.
-
Quantile regression models provide a wide picture of the conditional
distributions of the response variable by capturing the effect of the
covariates at different quantile levels. In most applications, the parametric
form of those conditional distributions is unknown and varies across the
covariate space, so fitting the given quantile levels simultaneously without
relying on parametric assumptions is crucial. In this work we propose a
Bayesian model for simultaneous linear quantile regression. More specifically,
we propose to model the conditional distributions by using random probability
measures known as quantile pyramids. Unlike many existing approaches, our
framework allows us to specify meaningful priors on the conditional
distributions, whilst retaining the flexibility afforded by the nonparametric
error distribution formulation. Simulation studies demonstrate the flexibility
of the proposed approach in estimating diverse scenarios, generally
outperforming other competitive methods. The method is particularly promising
for modelling the extremal quantiles. Applications to linear splines and
extreme value analysis are also explored through real data examples.
-
Statistical extreme value theory is concerned with the use of asymptotically
motivated models to describe the extreme values of a process. A number of
commonly used models are valid for observed data that exceed some high
threshold. However, in practice a suitable threshold is unknown and must be
determined for each analysis. While there are many threshold selection methods
for univariate extremes, there are relatively few that can be applied in the
multivariate setting. In addition, there are only a few Bayesian-based methods,
which are naturally attractive in the modelling of extremes due to data
scarcity. The use of Bayesian measures of surprise to determine suitable
thresholds for extreme value models is proposed. Such measures quantify the
level of support for the proposed extremal model and threshold, without the
need to specify any model alternatives. This approach is easily implemented for
both univariate and multivariate extremes.
-
Gamma-ray bursts (GRBs) are most probably powered by collimated relativistic
outflows (jets) from accreting black holes at cosmological distances. Bright
afterglows are produced when the outflow collides with the ambient medium.
Afterglow polarization directly probes the magnetic properties of the jet, when
measured minutes after the burst, and the geometric properties of the jet and
the ambient medium when measured hours to days after the burst. High values of
optical polarization detected minutes after burst in GRB 120308A indicate the
presence of large-scale ordered magnetic fields originating from the central
engine (the power source of the GRB). Theoretical models predict low degrees of
linear polarization and negligable circular polarization at late times, when
the energy in the original ejecta is quickly transferred to the ambient medium
and propagates farther into the medium as a blastwave. Here we report the
detection of circularly polarized optical light in the afterglow of GRB
121024A, measured 0.15 days after the burst. We show that the circular
polarization is intrinsic to the afterglow and unlikely to be produced by dust
scattering or plasma propagation effects. A possible explanation is to invoke
anisotropic (rather than the commonly assumed isotropic) electron pitch angle
distributions, and we suggest that new models are required to produce the
complex microphysics of realistic shocks in relativistic jets.
-
We study the dark nature of GRB 130528A through multi-wavelength observations
and conclude that the main reason for the optical darkness is local extinction
inside of the host galaxy. Automatic observations were performed at
BOOTES-4/MET robotic telescope. We also triggered target of opportunity (ToO)
observation at the OSN, IRAM PdBI and the GTC+OSIRIS. The host galaxy
photometric observations in optical to near-infrared (nIR) wavelengths were
achieved through large ground-based aperture telescopes, such as the 10.4m GTC,
the 4.2m WHT, 6m BTA, and the 2m LT. Based on these observations, spectral
energy distributions (SED) for the host galaxy and afterglow were constructed.
Thanks to mm observations at PdBI, we confirm the presence of a mm source
within the XRT error circle that faded over the course of our observations and
identify the host galaxy. However, we do not find any credible optical source
within early observations with BOOTES-4/MET and 1.5m OSN telescopes.
Spectroscopic observation of this galaxy by GTC showed a single faint emission
line that likely corresponds to [OII] 3727\{AA} at a redshift of 1.250+/-0.001
implying a SFR(M_sun/yr) > 6.18 M_sun/yr without correcting for dust
extinction. The probable extinction was revealed through analysis of the
afterglow SED, resulting in a value of AV >= ~ 0.9 at the rest frame, this is
comparable to extinction levels found among other dark GRBs. The SED of the
host galaxy is explained well (chi2/d.o.f.=0.564) by a luminous (MB=-21.16),
low-extinction (AV =0, rest frame), and aged (2.6 Gyr) stellar population. We
can explain this apparent contradiction in global and line-of-sight extinction
if the GRB birth place happened to lie in a local dense environment. In light
of having relatively small specific SFR (SSFR) ~ 5.3 M_sun/yr (L/L_star)-1,
this also could explain the age of the old stellar population of host galaxy.
-
Analysis of the properties of resonant modes in a reentrant cavity structure
comprising of a post and a ring was undertaken and verified experimentally. In
particular we show the existence of higher order reentrant cavity modes in such
a structure. Results show that the new cavity has two re-entrant modes, one
which has a better displacement sensitivity than the single post resonator and
the other with a reduced sensitivity. The more sensitive mode is better than
the single post resonator by a factor of 2 to 1.5 when the gap spacing is below
100 $ \mu$m. This type of cavity has the potential to operate as a highly
sensitive transducer for a variety of precision measurement applications, in
particular applications which require coupling to more than one sensitive
transducer mode.
-
Cylindrical re-entrant cavities are unique three-dimensional structures that
resonate with their electric and magnetic fields in separate parts of the
cavity. To further understand these devices, we undertake rigorous analysis of
the properties of the resonance using in-house developed Finite Element Method
(FEM) software capable of dealing with small gap structures of extreme aspect
ratio. Comparisons between the FEM method and experiments are consistent and we
illustrate where predictions using established lumped element models work well
and where they are limited. With the aid of the modeling we design a highly
tunable cavity that can be tuned from 2 GHz to 22 GHz just by inserting a post
into a fixed dimensioned cylindrical cavity. We show this is possible as the
mode structure transforms from a re-entrant mode during the tuning process to a
standard cylindrical Transverse Magnetic (TM) mode.
-
A new method is proposed for the correction of confidence intervals when the
original interval does not have the correct nominal coverage probabilities in
the frequentist sense. The proposed method is general and does not require any
distributional assumptions. It can be applied to both frequentist and Bayesian
inference where interval estimates are desired. We provide theoretical results
for the consistency of the proposed estimator, and give two complex examples,
on confidence interval correction for composite likelihood estimators and in
approximate Bayesian computation (ABC), to demonstrate the wide applicability
of the new method. Comparison is made with the double-bootstrap and other
methods of improving confidence interval coverage.
-
The coronal magnetic field is the primary driver of solar dynamic events.
Linear and circular polarization signals of certain infrared coronal emission
lines contain information about the magnetic field, and to access this
information, either a forward or an inversion method must be used. We study
three coronal magnetic configurations that are applicable to polar-crown
filament cavities by doing forward calculations to produce synthetic
polarization data. We analyze these forward data to determine the
distinguishing characteristics of each model. We conclude that it is possible
to distinguish between cylindrical flux ropes, spheromak flux ropes, and
sheared arcades using coronal polarization measurements. If one of these models
is found to be consistent with observational measurements, it will mean
positive identification of the magnetic morphology that surrounds certain
quiescent filaments, which will lead to a greater understanding of how they
form and why they erupt.
-
The Coronal Multi-channel Polarimeter (CoMP) obtains daily full-Sun
above-the-limb coronal observations in linear polarization, allowing for the
first time a diagnostic of the coronal magnetic field direction in quiescent
prominence cavities. We find that these cavities consistently possess a
characteristic "lagomorphic" signature in linear polarization indicating twist
or shear extending up into the cavity above the neutral line. We demonstrate
that such a signature may be explained by a magnetic flux-rope model, a
topology with implications for solar eruptions. We find corroborating evidence
for a flux rope structure in the pattern of concentric rings within cavities
seen in CoMP line-of-sight velocity.
-
Bayes linear analysis and approximate Bayesian computation (ABC) are
techniques commonly used in the Bayesian analysis of complex models. In this
article we connect these ideas by demonstrating that regression-adjustment ABC
algorithms produce samples for which first and second order moment summaries
approximate adjusted expectation and variance for a Bayes linear analysis. This
gives regression-adjustment methods a useful interpretation and role in
exploratory analysis in high-dimensional problems. As a result, we propose a
new method for combining high-dimensional, regression-adjustment ABC with
lower-dimensional approaches (such as using MCMC for ABC). This method first
obtains a rough estimate of the joint posterior via regression-adjustment ABC,
and then estimates each univariate marginal posterior distribution separately
in a lower-dimensional analysis. The marginal distributions of the initial
estimate are then modified to equal the separately estimated marginals, thereby
providing an improved estimate of the joint posterior. We illustrate this
method with several examples. Supplementary materials for this article are
available online.
-
Approximate Bayesian computation (ABC) methods, which are applicable when the
likelihood is difficult or impossible to calculate, are an active topic of
current research. Most current ABC algorithms directly approximate the
posterior distribution, but an alternative, less common strategy is to
approximate the likelihood function. This has several advantages. First, in
some problems, it is easier to approximate the likelihood than to approximate
the posterior. Second, an approximation to the likelihood allows reference
analyses to be constructed based solely on the likelihood. Third, it is
straightforward to perform sensitivity analyses for several different choices
of prior once an approximation to the likelihood is constructed, which needs to
be done only once. The contribution of the present paper is to consider
regression density estimation techniques to approximate the likelihood in the
ABC setting. Our likelihood approximations build on recently developed marginal
adaptation density estimators by extending them for conditional density
estimation. Our approach facilitates reference Bayesian inference, as well as
frequentist inference. The method is demonstrated via a challenging problem of
inference for stereological extremes, where we perform both frequentist and
Bayesian inference.
-
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.
-
Using a 3D MHD simulation, we model the quasi-static evolution and the onset
of eruption of a coronal flux rope. The simulation begins with a twisted flux
rope emerging at the lower boundary and pushing into a pre-existing coronal
potential arcade field. At a chosen time the emergence is stopped with the
lower boundary taken to be rigid. Then the coronal flux rope settles into a
quasi-static rise phase during which an underlying, central sigmoid-shaped
current layer forms along the so called hyperbolic flux tube (HFT), a
generalization of the X-line configuration. Reconnections in the dissipating
current layer effectively add twisted flux to the flux rope and thus allow it
to rise quasi-statically, even though the magnetic energy is decreasing as the
system relaxes. We examine the thermal features produced by the current layer
formation and the associated "tether-cutting" reconnections as a result of
heating and field aligned thermal conduction. It is found that a central hot,
low-density channel containing reconnected, twisted flux threading under the
flux rope axis forms on top of the central current layer. When viewed in the
line of sight roughly aligned with the hot channel (which is roughly along the
neutral line), the central current layer appears as a high-density vertical
column with upward extensions as a "U" shaped dense shell enclosing a central
hot, low-density void. Such thermal features have been observed within coronal
prominence cavities. Our MHD simulation suggests that they are the signatures
of the development of the HFT topology and the associated tether-cutting
reconnections, and that the central void grows and rises with the
reconnections, until the flux rope reaches the critical height for the onset of
the torus instability and dynamic eruption ensues.
-
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.
-
We present the results of a LIGO search for gravitational waves (GWs)
associated with GRB 051103, a short-duration hard-spectrum gamma-ray burst
(GRB) whose electromagnetically determined sky position is coincident with the
spiral galaxy M81, which is 3.6 Mpc from Earth. Possible progenitors for
short-hard GRBs include compact object mergers and soft gamma repeater (SGR)
giant flares. A merger progenitor would produce a characteristic GW signal that
should be detectable at the distance of M81, while GW emission from an SGR is
not expected to be detectable at that distance. We found no evidence of a GW
signal associated with GRB 051103. Assuming weakly beamed gamma-ray emission
with a jet semi-angle of 30 deg we exclude a binary neutron star merger in M81
as the progenitor with a confidence of 98%. Neutron star-black hole mergers are
excluded with > 99% confidence. If the event occurred in M81 our findings
support the the hypothesis that GRB 051103 was due to an SGR giant flare,
making it the most distant extragalactic magnetar observed to date.
-
Quantile regression has received increased attention in the statistics
community in recent years. This article adapts an auxiliary variable method,
commonly used in Bayesian variable selection for mean regression models, to the
fitting of quantile regression curves. We focus on the fitting of regression
splines, with unknown number and location of knots. We provide an efficient
algorithm with Metropolis-Hastings updates whose tuning is fully automated. The
method is tested on simulated and real examples and its extension to additive
models is described. Finally we propose a simple postprocessing procedure to
deal with the problem of the crossing of multiple separately estimated quantile
curves.
In the context of engineered quantum systems, there is a demand for superconducting tunable devices able to operate with high Q-factors at power levels equivalent to only a few photons. In this work, we developed a 3D microwave reentrant cavity with such characteristics ready to provide a very fine-tuning of a high-Q resonant mode over a large dynamic range. This system has an electronic tuning mechanism based on a mechanically amplified piezoelectric actuator, which controls the resonator dominant mode frequency by changing the cavity narrow gap by very small displacements. Experiments were conducted at room and dilution refrigerator temperatures showing a large dynamic range up to 4 GHz and 1 GHz, respectively, and were compared to a FEM model simulated data. At elevated microwave power input, nonlinear thermal effects were observed to destroy the superconductivity of the cavity due to the large electric fields generated in the small gap of the reentrant cavity.
We consider the application of approximate Bayesian Computation (ABC) in the context of medical imaging data. We consider the parameter estimation of compartmental models in PET imaging analysis, and provide a simple ABC algorithm for its estimation. We demonstrate the utility of the proposed estimation methods on a neurotransmitter response model, and compare our approach to existing methods.
This Supplement provides supporting material for arXiv:1602.08492 . We briefly summarize past electromagnetic (EM) follow-up efforts as well as the organization and policy of the current EM follow-up program. We compare the four probability sky maps produced for the gravitational-wave transient GW150914, and provide additional details of the EM follow-up observations that were performed in the different bands.
A gravitational-wave (GW) transient was identified in data recorded by the Advanced Laser Interferometer Gravitational-wave Observatory (LIGO) detectors on 2015 September 14. The event, initially designated G184098 and later given the name GW150914, is described in detail elsewhere. By prior arrangement, preliminary estimates of the time, significance, and sky location of the event were shared with 63 teams of observers covering radio, optical, near-infrared, X-ray, and gamma-ray wavelengths with ground- and space-based facilities. In this Letter we describe the low-latency analysis of the GW data and present the sky localization of the first observed compact binary merger. We summarize the follow-up observations reported by 25 teams via private Gamma-ray Coordinates Network circulars, giving an overview of the participating facilities, the GW sky localization coverage, the timeline and depth of the observations. As this event turned out to be a binary black hole merger, there is little expectation of a detectable electromagnetic (EM) signature. Nevertheless, this first broadband campaign to search for a counterpart of an Advanced LIGO source represents a milestone and highlights the broad capabilities of the transient astronomy community and the observing strategies that have been developed to pursue neutron star binary merger events. Detailed investigations of the EM data and results of the EM follow-up campaign are being disseminated in papers by the individual teams.
Quantile regression models provide a wide picture of the conditional distributions of the response variable by capturing the effect of the covariates at different quantile levels. In most applications, the parametric form of those conditional distributions is unknown and varies across the covariate space, so fitting the given quantile levels simultaneously without relying on parametric assumptions is crucial. In this work we propose a Bayesian model for simultaneous linear quantile regression. More specifically, we propose to model the conditional distributions by using random probability measures known as quantile pyramids. Unlike many existing approaches, our framework allows us to specify meaningful priors on the conditional distributions, whilst retaining the flexibility afforded by the nonparametric error distribution formulation. Simulation studies demonstrate the flexibility of the proposed approach in estimating diverse scenarios, generally outperforming other competitive methods. The method is particularly promising for modelling the extremal quantiles. Applications to linear splines and extreme value analysis are also explored through real data examples.
Statistical extreme value theory is concerned with the use of asymptotically motivated models to describe the extreme values of a process. A number of commonly used models are valid for observed data that exceed some high threshold. However, in practice a suitable threshold is unknown and must be determined for each analysis. While there are many threshold selection methods for univariate extremes, there are relatively few that can be applied in the multivariate setting. In addition, there are only a few Bayesian-based methods, which are naturally attractive in the modelling of extremes due to data scarcity. The use of Bayesian measures of surprise to determine suitable thresholds for extreme value models is proposed. Such measures quantify the level of support for the proposed extremal model and threshold, without the need to specify any model alternatives. This approach is easily implemented for both univariate and multivariate extremes.
Gamma-ray bursts (GRBs) are most probably powered by collimated relativistic outflows (jets) from accreting black holes at cosmological distances. Bright afterglows are produced when the outflow collides with the ambient medium. Afterglow polarization directly probes the magnetic properties of the jet, when measured minutes after the burst, and the geometric properties of the jet and the ambient medium when measured hours to days after the burst. High values of optical polarization detected minutes after burst in GRB 120308A indicate the presence of large-scale ordered magnetic fields originating from the central engine (the power source of the GRB). Theoretical models predict low degrees of linear polarization and negligable circular polarization at late times, when the energy in the original ejecta is quickly transferred to the ambient medium and propagates farther into the medium as a blastwave. Here we report the detection of circularly polarized optical light in the afterglow of GRB 121024A, measured 0.15 days after the burst. We show that the circular polarization is intrinsic to the afterglow and unlikely to be produced by dust scattering or plasma propagation effects. A possible explanation is to invoke anisotropic (rather than the commonly assumed isotropic) electron pitch angle distributions, and we suggest that new models are required to produce the complex microphysics of realistic shocks in relativistic jets.
We study the dark nature of GRB 130528A through multi-wavelength observations and conclude that the main reason for the optical darkness is local extinction inside of the host galaxy. Automatic observations were performed at BOOTES-4/MET robotic telescope. We also triggered target of opportunity (ToO) observation at the OSN, IRAM PdBI and the GTC+OSIRIS. The host galaxy photometric observations in optical to near-infrared (nIR) wavelengths were achieved through large ground-based aperture telescopes, such as the 10.4m GTC, the 4.2m WHT, 6m BTA, and the 2m LT. Based on these observations, spectral energy distributions (SED) for the host galaxy and afterglow were constructed. Thanks to mm observations at PdBI, we confirm the presence of a mm source within the XRT error circle that faded over the course of our observations and identify the host galaxy. However, we do not find any credible optical source within early observations with BOOTES-4/MET and 1.5m OSN telescopes. Spectroscopic observation of this galaxy by GTC showed a single faint emission line that likely corresponds to [OII] 3727\{AA} at a redshift of 1.250+/-0.001 implying a SFR(M_sun/yr) > 6.18 M_sun/yr without correcting for dust extinction. The probable extinction was revealed through analysis of the afterglow SED, resulting in a value of AV >= ~ 0.9 at the rest frame, this is comparable to extinction levels found among other dark GRBs. The SED of the host galaxy is explained well (chi2/d.o.f.=0.564) by a luminous (MB=-21.16), low-extinction (AV =0, rest frame), and aged (2.6 Gyr) stellar population. We can explain this apparent contradiction in global and line-of-sight extinction if the GRB birth place happened to lie in a local dense environment. In light of having relatively small specific SFR (SSFR) ~ 5.3 M_sun/yr (L/L_star)-1, this also could explain the age of the old stellar population of host galaxy.
The coronal magnetic field is the primary driver of solar dynamic events. Linear and circular polarization signals of certain infrared coronal emission lines contain information about the magnetic field, and to access this information, either a forward or an inversion method must be used. We study three coronal magnetic configurations that are applicable to polar-crown filament cavities by doing forward calculations to produce synthetic polarization data. We analyze these forward data to determine the distinguishing characteristics of each model. We conclude that it is possible to distinguish between cylindrical flux ropes, spheromak flux ropes, and sheared arcades using coronal polarization measurements. If one of these models is found to be consistent with observational measurements, it will mean positive identification of the magnetic morphology that surrounds certain quiescent filaments, which will lead to a greater understanding of how they form and why they erupt.
The Coronal Multi-channel Polarimeter (CoMP) obtains daily full-Sun above-the-limb coronal observations in linear polarization, allowing for the first time a diagnostic of the coronal magnetic field direction in quiescent prominence cavities. We find that these cavities consistently possess a characteristic "lagomorphic" signature in linear polarization indicating twist or shear extending up into the cavity above the neutral line. We demonstrate that such a signature may be explained by a magnetic flux-rope model, a topology with implications for solar eruptions. We find corroborating evidence for a flux rope structure in the pattern of concentric rings within cavities seen in CoMP line-of-sight velocity.
Bayes linear analysis and approximate Bayesian computation (ABC) are techniques commonly used in the Bayesian analysis of complex models. In this article we connect these ideas by demonstrating that regression-adjustment ABC algorithms produce samples for which first and second order moment summaries approximate adjusted expectation and variance for a Bayes linear analysis. This gives regression-adjustment methods a useful interpretation and role in exploratory analysis in high-dimensional problems. As a result, we propose a new method for combining high-dimensional, regression-adjustment ABC with lower-dimensional approaches (such as using MCMC for ABC). This method first obtains a rough estimate of the joint posterior via regression-adjustment ABC, and then estimates each univariate marginal posterior distribution separately in a lower-dimensional analysis. The marginal distributions of the initial estimate are then modified to equal the separately estimated marginals, thereby providing an improved estimate of the joint posterior. We illustrate this method with several examples. Supplementary materials for this article are available online.
Approximate Bayesian computation (ABC) methods, which are applicable when the likelihood is difficult or impossible to calculate, are an active topic of current research. Most current ABC algorithms directly approximate the posterior distribution, but an alternative, less common strategy is to approximate the likelihood function. This has several advantages. First, in some problems, it is easier to approximate the likelihood than to approximate the posterior. Second, an approximation to the likelihood allows reference analyses to be constructed based solely on the likelihood. Third, it is straightforward to perform sensitivity analyses for several different choices of prior once an approximation to the likelihood is constructed, which needs to be done only once. The contribution of the present paper is to consider regression density estimation techniques to approximate the likelihood in the ABC setting. Our likelihood approximations build on recently developed marginal adaptation density estimators by extending them for conditional density estimation. Our approach facilitates reference Bayesian inference, as well as frequentist inference. The method is demonstrated via a challenging problem of inference for stereological extremes, where we perform both frequentist and Bayesian inference.
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.
Using a 3D MHD simulation, we model the quasi-static evolution and the onset of eruption of a coronal flux rope. The simulation begins with a twisted flux rope emerging at the lower boundary and pushing into a pre-existing coronal potential arcade field. At a chosen time the emergence is stopped with the lower boundary taken to be rigid. Then the coronal flux rope settles into a quasi-static rise phase during which an underlying, central sigmoid-shaped current layer forms along the so called hyperbolic flux tube (HFT), a generalization of the X-line configuration. Reconnections in the dissipating current layer effectively add twisted flux to the flux rope and thus allow it to rise quasi-statically, even though the magnetic energy is decreasing as the system relaxes. We examine the thermal features produced by the current layer formation and the associated "tether-cutting" reconnections as a result of heating and field aligned thermal conduction. It is found that a central hot, low-density channel containing reconnected, twisted flux threading under the flux rope axis forms on top of the central current layer. When viewed in the line of sight roughly aligned with the hot channel (which is roughly along the neutral line), the central current layer appears as a high-density vertical column with upward extensions as a "U" shaped dense shell enclosing a central hot, low-density void. Such thermal features have been observed within coronal prominence cavities. Our MHD simulation suggests that they are the signatures of the development of the HFT topology and the associated tether-cutting reconnections, and that the central void grows and rises with the reconnections, until the flux rope reaches the critical height for the onset of the torus instability and dynamic eruption ensues.
We present the results of a LIGO search for gravitational waves (GWs) associated with GRB 051103, a short-duration hard-spectrum gamma-ray burst (GRB) whose electromagnetically determined sky position is coincident with the spiral galaxy M81, which is 3.6 Mpc from Earth. Possible progenitors for short-hard GRBs include compact object mergers and soft gamma repeater (SGR) giant flares. A merger progenitor would produce a characteristic GW signal that should be detectable at the distance of M81, while GW emission from an SGR is not expected to be detectable at that distance. We found no evidence of a GW signal associated with GRB 051103. Assuming weakly beamed gamma-ray emission with a jet semi-angle of 30 deg we exclude a binary neutron star merger in M81 as the progenitor with a confidence of 98%. Neutron star-black hole mergers are excluded with > 99% confidence. If the event occurred in M81 our findings support the the hypothesis that GRB 051103 was due to an SGR giant flare, making it the most distant extragalactic magnetar observed to date.
Quantile regression has received increased attention in the statistics community in recent years. This article adapts an auxiliary variable method, commonly used in Bayesian variable selection for mean regression models, to the fitting of quantile regression curves. We focus on the fitting of regression splines, with unknown number and location of knots. We provide an efficient algorithm with Metropolis-Hastings updates whose tuning is fully automated. The method is tested on simulated and real examples and its extension to additive models is described. Finally we propose a simple postprocessing procedure to deal with the problem of the crossing of multiple separately estimated quantile curves.