
The properties of precessing, coalescing binary black holes are presently
inferred through comparison with two approximate models of compact binary
coalescence. In this work we show these two models often disagree substantially
when binaries have modestly large spins ($a\gtrsim 0.4$) and modest mass ratios
($q\gtrsim 2$). We demonstrate these disagreements using standard figures of
merit and the parameters inferred for recent detections of binary black holes.
By comparing to numerical relativity, we confirm these disagreements reflect
systematic errors. We provide concrete examples to demonstrate that these
systematic errors can significantly impact inferences about astrophysically
significant binary parameters. For the immediate future, parameter inference
for binary black holes should be performed with multiple models (including
numerical relativity), and carefully validated by performing inference under
controlled circumstances with similar synthetic events.

We present constraints on the annihilation cross section of WIMP dark matter
based on the joint statistical analysis of four dwarf galaxies with VERITAS.
These results are derived from an optimized photon weighting statistical
technique that improves on standard imaging atmospheric Cherenkov telescope
(IACT) analyses by utilizing the spectral and spatial properties of individual
photon events. We report on the results of $\sim$230 hours of observations of
five dwarf galaxies and the joint statistical analysis of four of the dwarf
galaxies. We find no evidence of gammaray emission from any individual dwarf
nor in the joint analysis. The derived upper limit on the dark matter
annihilation cross section from the joint analysis is $1.35\times 10^{23}
{\mathrm{ cm^3s^{1}}}$ at 1 TeV for the bottom quark ($b\bar{b}$) final state,
$2.85\times 10^{24}{\mathrm{ cm^3s^{1}}}$ at 1 TeV for the tau lepton
($\tau^{+}\tau^{}$) final state and $1.32\times 10^{25}{\mathrm{
cm^3s^{1}}}$ at 1 TeV for the gauge boson ($\gamma\gamma$) final state.

We present an extensive study of the BL Lac object Mrk 501 based on a data
set collected during the multiinstrument campaign spanning from 2009 March 15
to 2009 August 1 which includes, among other instruments, MAGIC, VERITAS,
Whipple 10m, FermiLAT, RXTE, Swift, GASPWEBT and VLBA. We find an increase
in the fractional variability with energy, while no significant interband
correlations of flux changes are found in the acquired data set. The higher
variability in the very high energy (>100 GeV, VHE) gammaray emission and the
lack of correlation with the Xray emission indicate that the highestenergy
electrons that are responsible for the VHE gammarays do not make a dominant
contribution to the ~1 keV emission. Alternatively, there could be a very
variable component contributing to the VHE gammaray emission in addition to
that coming from the synchrotron selfCompton (SSC) scenarios. The space of SSC
model parameters is probed following a dedicated gridscan strategy, allowing
for a wide range of models to be tested and offering a study of the degeneracy
of modeltodata agreement in the individual model parameters. We find that
there is some degeneracy in both the onezone and the twozone SSC scenarios
that were probed, with several combinations of model parameters yielding a
similar modeltodata agreement, and some parameters better constrained than
others. The SSC model gridscan shows that the flaring activity around 2009 May
22 cannot be modeled adequately with a onezone SSC scenario, while it can be
suitably described within a twoindependentzone SSC scenario. The observation
of an electric vector polarization angle rotation coincident with the gammaray
flare from 2009 May 1 resembles those reported previously for low frequency
peaked blazars, hence suggesting that there are many similarities in the
flaring mechanisms of blazars with different jet properties.

A flare from the TeV blazar Mrk 421, occurring in March 2010, was observed
for 13 consecutive days from radio to very high energy (VHE, E > 100 GeV)
gammarays with MAGIC, VERITAS, Whipple, FermiLAT, MAXI, RXTE, Swift,
GASPWEBT, and several optical and radio telescopes. We model the dayscale
SEDs with onezone and twozone synchrotron selfCompton (SSC) models,
investigate the physical parameters, and evaluate whether the observed
broadband SED variability can be associated to variations in the relativistic
particle population.
Flux variability was remarkable in the Xray and VHE bands while it was minor
or not significant in the other bands. The onezone SSC model can describe
reasonably well the SED of each day for the 13 consecutive days. This flaring
activity is also very well described by a twozone SSC model, where one zone is
responsible for the quiescent emission while the other smaller zone, which is
spatially separated from the first one, contributes to the dailyvariable
emission occurring in Xrays and VHE gammarays.
Both the onezone SSC and the twozone SSC models can describe the daily SEDs
via the variation of only four or five model parameters, under the hypothesis
that the variability is associated mostly to the underlying particle
population. This shows that the particle acceleration and cooling mechanism
producing the radiating particles could be the main one responsible for the
broadband SED variations during the flaring episodes in blazars. The twozone
SSC model provides a better agreement to the observed SED at the narrow peaks
of the low and highenergy bumps during the highest activity, although the
reported onezone SSC model could be further improved by the variation of the
parameters related to the emitting region itself ($\delta$, $B$ and $R$), in
addition to the parameters related to the particle population.

We performed a 4.5month multiinstrument campaign (from radio to VHE gamma
rays) on Mrk421 between January 2009 and June 2009, which included VLBA,
FGAMMA, GASPWEBT, Swift, RXTE, FermiLAT, MAGIC, and Whipple, among other
instruments and collaborations. Mrk421 was found in its typical (nonflaring)
activity state, with a VHE flux of about half that of the Crab Nebula, yet the
light curves show significant variability at all wavelengths, the highest
variability being in the Xrays. We determined the power spectral densities
(PSD) at most wavelengths and found that all PSDs can be described by
powerlaws without a break, and with indices consistent with pink/rednoise
behavior. We observed a harderwhenbrighter behavior in the Xray spectra and
measured a positive correlation between VHE and Xray fluxes with zero time
lag. Such characteristics have been reported many times during flaring
activity, but here they are reported for the first time in the nonflaring
state. We also observed an overall anticorrelation between optical/UV and
Xrays extending over the duration of the campaign.
The harderwhenbrighter behavior in the Xray spectra and the measured
positive Xray/VHE correlation during the 2009 multiwavelength campaign
suggests that the physical processes dominating the emission during nonflaring
states have similarities with those occurring during flaring activity. In
particular, this observation supports leptonic scenarios as being responsible
for the emission of Mrk421 during nonflaring activity. Such a temporally
extended Xray/VHE correlation is not driven by any single flaring event, and
hence is difficult to explain within the standard hadronic scenarios. The
highest variability is observed in the Xray band, which, within the onezone
synchrotron selfCompton scenario, indicates that the electron energy
distribution is most variable at the highest energies.

We present the results of 71.6 hours of observations of the Geminga pulsar
(PSR J0633+1746) with the VERITAS veryhighenergy gammaray telescope array.
Data taken with VERITAS between November 2007 and February 2013 were
phasefolded using a Geminga pulsar timing solution derived from data recorded
by the XMM\emph{Newton} and \emph{Fermi}LAT space telescopes. No significant
pulsed emission above 100 GeV is observed, and we report upper limits at the
95% confidence level on the integral flux above 135 GeV (spectral analysis
threshold) of 4.0$\times10^{13}$ s$^{1}$ cm$^{2}$ and 1.7$\times10^{13}$
s$^{1}$ cm$^{2}$ for the two principal peaks in the emission profile. These
upper limits, placed in context with phaseresolved spectral energy
distributions determined from five years of data from the \emph{Fermi}LAT,
constrain possible hardening of the Geminga pulsar emission spectra above
$\sim$50 GeV.

Multiple epidemiological models have been proposed to predict the spread of
Ebola in West Africa. These models include consideration of countermeasures
meant to slow and, eventually, stop the spread of the disease. Here, we examine
one component of Ebola dynamics that is of growing concern  the transmission
of Ebola from the dead to the living. We do so by applying the toolkit of
mathematical epidemiology to analyze the consequences of postdeath
transmission. We show that underlying disease parameters cannot be inferred
with confidence from earlystage incidence data (that is, they are not
"identifiable") because different parameter combinations can produce virtually
the same epidemic trajectory. Despite this identifiability problem, we find
robustly that inferences that don't account for postdeath transmission tend to
underestimate the basic reproductive number  thus, given the observed rate of
epidemic growth, larger amounts of postdeath transmission imply larger
reproductive numbers. From a control perspective, we explain how improvements
in reducing postdeath transmission of Ebola may reduce the overall epidemic
spread and scope substantially. Increased attention to the proportion of
postdeath transmission has the potential to aid both in projecting the course
of the epidemic and in evaluating a portfolio of control strategies.

Mrk 501 is one of the brightest blazars at TeV energies and has been
extensively studied since its first VHE detection in 1996. Our goal is to
characterize in detail the source gammaray emission, together with the
radiotoXray emission, during the nonflaring (low) activity, which is less
often studied than the occasional flaring (high) activity. We organized a
multiwavelength (MW) campaign on Mrk 501 between March and May 2008. This
multiinstrument effort included the most sensitive VHE gammaray instruments
in the northern hemisphere, namely the imaging atmospheric Cherenkov telescopes
MAGIC and VERITAS, as well as Swift, RXTE, the FGAMMA, GASPWEBT, and other
collaborations and instruments. Mrk 501 was found to be in a low state of
activity during the campaign, with a VHE flux in the range of 10%20% of the
Crab nebula flux. Nevertheless, significant flux variations were detected with
various instruments, with a trend of increasing variability with energy. The
broadband spectral energy distribution during the two different emission states
of the campaign can be adequately described within the homogeneous onezone
synchrotron selfCompton model, with the (slightly) higher state described by
an increase in the electron number density. This agrees with previous studies
of the broadband emission of this source during flaring and nonflaring states.
We report for the first time a tentative XraytoVHE correlation during a low
VHE activity. Although marginally significant, this positive correlation
between Xray and VHE, which has been reported many times during flaring
activity, suggests that the mechanisms that dominate the Xray/VHE emission
during nonflaringactivity are not substantially different from those that are
responsible for the emission during flaring activity.

The NumericalRelativityAnalyticalRelativity (NRAR) collaboration is a
joint effort between members of the numerical relativity, analytical relativity
and gravitationalwave data analysis communities. The goal of the NRAR
collaboration is to produce numericalrelativity simulations of compact
binaries and use them to develop accurate analytical templates for the
LIGO/Virgo Collaboration to use in detecting gravitationalwave signals and
extracting astrophysical information from them. We describe the results of the
first stage of the NRAR project, which focused on producing an initial set of
numerical waveforms from binary black holes with moderate mass ratios and
spins, as well as one nonspinning binary configuration which has a mass ratio
of 10. All of the numerical waveforms are analysed in a uniform and consistent
manner, with numerical errors evaluated using an analysis code created by
members of the NRAR collaboration. We compare previouslycalibrated,
nonprecessing analytical waveforms, notably the effectiveonebody (EOB) and
phenomenological template families, to the newlyproduced numerical waveforms.
We find that when the binary's total mass is ~100200 solar masses, current EOB
and phenomenological models of spinning, nonprecessing binary waveforms have
overlaps above 99% (for advanced LIGO) with all of the nonprecessingbinary
numerical waveforms with mass ratios <= 4, when maximizing over binary
parameters. This implies that the loss of event rate due to modelling error is
below 3%. Moreover, the nonspinning EOB waveforms previously calibrated to
five nonspinning waveforms with mass ratio smaller than 6 have overlaps above
99.7% with the numerical waveform with a mass ratio of 10, without even
maximizing on the binary parameters.

We address the statistical mechanics of randomly and permanently crosslinked
networks. We develop a theoretical framework (vulcanization theory) which can
be used to systematically analyze the correlation between the statistical
properties of random networks and their histories of formation. Generalizing
the original idea of Deam and Edwards, we consider an instantaneous
crosslinking process, where all crosslinkers (modeled as Gaussian springs) are
introduced randomly at once in an equilibrium liquid state, referred to as the
preparation state. The probability that two functional sites are crosslinked by
a spring exponentially decreases with their distance squared. After formally
averaging over network connectivity, we obtained an effective theory with all
degrees of freedom replicated 1 + n times. Two thermodynamic ensembles, the
preparation ensemble and the measurement ensemble, naturally appear in this
theory. The former describes the thermodynamic fluctuations in the state of
preparation, while the latter describes the thermodynamic fluctuations in the
state of measurement. We classify various correlation functions and discuss
their physical significances. In particular, the memory correlation functions
characterize how the properties of networks depend on their history of
formation, and are the hallmark properties of all randomly crosslinked
materials. We clarify the essential difference between our approach and that of
DeamEdwards, discuss the saddlepoint order parameters and its physical
significance. Finally we also discuss the connection between saddlepoint
approximation of vulcanization theory, and the classical theory of rubber
elasticity as well as the neoclassical theory of nematic elastomers.

We discuss the prospects of doing tests of Lorentz invariance with gammarays
observed with present and future ground based gammaray observatories.

Toeplitzstructured linear systems arise often in practical engineering
problems. Correspondingly, a number of algorithms have been developed that
exploit Toeplitz structure to gain computational efficiency when solving these
systems. The earliest "fast" algorithms for Toeplitz systems required O(n^2)
operations, while more recent "superfast" algorithms reduce the cost to O(n
(log n)^2) or below. In this work, we present a superfast algorithm for
Tikhonov regularization of Toeplitz systems. Using an
"extensionandtransformation" technique, our algorithm translates a
Tikhonovregularized Toeplitz system into a type of specialized polynomial
problem known as tangential interpolation. Under this formulation, we can
compute the solution in only O(n (log n)^2) operations. We use numerical
simulations to demonstrate our algorithm's complexity and verify that it
returns stable solutions.

We discuss the interdependence of resource state, measurement setting and
temporal order in measurementbased quantum computation. The possible temporal
orders of measurement events are constrained by the principle that the
randomness inherent in quantum measurement should not affect the outcome of the
computation. We provide a classification for all temporal relations among
measurement events compatible with a given initial stabilizer state and
measurement setting, in terms of a matroid. Conversely, we show that classical
processing relations necessary for turning the local measurement outcomes into
computational output determine the resource state and measurement setting up to
local equivalence. Further, we find a symmetry transformation related to local
complementation that leaves the temporal relations invariant.

In this paper, we use semidefinite programming and generalized principal
component analysis (GPCA) to distinguish between two or more different facial
expressions. In the first step, semidefinite programming is used to reduce the
dimension of the image data and "unfold" the manifold which the data points
(corresponding to facial expressions) reside on. Next, GPCA is used to fit a
series of subspaces to the data points and associate each data point with a
subspace. Data points that belong to the same subspace are claimed to belong to
the same facial expression category. An example is provided.

We report the observation of a Fano resonance between continuum Mie
scattering and a narrow Bragg band in synthetic opal photonic crystals. The
resonance leads to a transmission spectrum exhibiting a Bragg dip with an
asymmetric profile, which can be tunably reversed to a Bragg rise. The Fano
asymmetry parameter is linked with the dielectric contrast between the
permittivity of the filler and the specific value determined by the opal
matrix. The existence of the Fano resonance is directly related to disorder due
to nonuniformity of aSiO2 opal spheres. Proposed theoretical "quasi3D" model
produces results in excellent agreement with the experimental data.

Vibrational energy flows unevenly in molecules, repeatedly going back and
forth between trapping and roaming. We identify bottlenecks between diffusive
and chaotic behavior, and describe generic mechanisms of these transitions,
taking the carbonyl sulphide molecule OCS as a case study. The bottlenecks are
found to be lowerdimensional tori; their bifurcations and unstable manifolds
govern the transition mechanisms.

Segal proposed ultraquantum commutation relations with two ultraquantum
constants hbar' and hbar'' besides Planck's quantum constant hbar with a
variable i. The Heisenberg quantum algebra is a contraction  in a more general
sense than that of Inonu and Wigner  of the Segal ultraquantum algebra. The
usual constant i arises as a vacuum orderparameter in the quantum limit where
hbar' and hbar'' approach zero. One physical consequence is a discrete spectrum
for canonical variable and spacetime coordinates. Another is an
interconvention of time and energy accompnying spacetime meltdown (disorder),
with a fundamental conversion factor of some kilograms of energy per second.

Formation and evolution of fragmentation instabilities in fractal islands,
obtained by deposition of silver clusters on graphite, are studied. The
fragmentation dynamics and subsequent relaxation to the equilibrium shapes are
controlled by the deposition conditions and cluster composition. Sharing common
features with other materials' breakup phenomena, the fragmentation instability
is governed by the lengthtowidth ratio of the fractal arms.

We study analytically the effect of a constant magnetic field on the dynamics
of a two dimensional Josephson array. The magnetic field induces spatially
dependent states and coupling between rows, even in the absence of an external
load. Numerical simulations support these conclusions.