• Systematic challenges for future gravitational wave measurements of precessing binary black holes(1709.03095)

Sept. 29, 2017 gr-qc, astro-ph.HE
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 gamma-ray 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 multi-instrument campaign spanning from 2009 March 15 to 2009 August 1 which includes, among other instruments, MAGIC, VERITAS, Whipple 10-m, Fermi-LAT, RXTE, Swift, GASP-WEBT 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) gamma-ray emission and the lack of correlation with the X-ray emission indicate that the highest-energy electrons that are responsible for the VHE gamma-rays do not make a dominant contribution to the ~1 keV emission. Alternatively, there could be a very variable component contributing to the VHE gamma-ray emission in addition to that coming from the synchrotron self-Compton (SSC) scenarios. The space of SSC model parameters is probed following a dedicated grid-scan strategy, allowing for a wide range of models to be tested and offering a study of the degeneracy of model-to-data agreement in the individual model parameters. We find that there is some degeneracy in both the one-zone and the two-zone SSC scenarios that were probed, with several combinations of model parameters yielding a similar model-to-data agreement, and some parameters better constrained than others. The SSC model grid-scan shows that the flaring activity around 2009 May 22 cannot be modeled adequately with a one-zone SSC scenario, while it can be suitably described within a two-independent-zone SSC scenario. The observation of an electric vector polarization angle rotation coincident with the gamma-ray 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) gamma-rays with MAGIC, VERITAS, Whipple, FermiLAT, MAXI, RXTE, Swift, GASP-WEBT, and several optical and radio telescopes. We model the day-scale SEDs with one-zone and two-zone synchrotron self-Compton (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 X-ray and VHE bands while it was minor or not significant in the other bands. The one-zone 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 two-zone 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 daily-variable emission occurring in X-rays and VHE gamma-rays. Both the one-zone SSC and the two-zone 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 two-zone SSC model provides a better agreement to the observed SED at the narrow peaks of the low- and high-energy bumps during the highest activity, although the reported one-zone 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.5-month multi-instrument campaign (from radio to VHE gamma rays) on Mrk421 between January 2009 and June 2009, which included VLBA, F-GAMMA, GASP-WEBT, Swift, RXTE, Fermi-LAT, MAGIC, and Whipple, among other instruments and collaborations. Mrk421 was found in its typical (non-flaring) 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 X-rays. We determined the power spectral densities (PSD) at most wavelengths and found that all PSDs can be described by power-laws without a break, and with indices consistent with pink/red-noise behavior. We observed a harder-when-brighter behavior in the X-ray spectra and measured a positive correlation between VHE and X-ray 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 non-flaring state. We also observed an overall anti-correlation between optical/UV and X-rays extending over the duration of the campaign. The harder-when-brighter behavior in the X-ray spectra and the measured positive X-ray/VHE correlation during the 2009 multi-wavelength campaign suggests that the physical processes dominating the emission during non-flaring 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 non-flaring activity. Such a temporally extended X-ray/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 X-ray band, which, within the one-zone synchrotron self-Compton 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 very-high-energy gamma-ray telescope array. Data taken with VERITAS between November 2007 and February 2013 were phase-folded 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 phase-resolved 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.
• Post-death Transmission of Ebola: Challenges for Inference and Opportunities for Control(1411.3435)

Nov. 13, 2014 q-bio.PE
Multiple epidemiological models have been proposed to predict the spread of Ebola in West Africa. These models include consideration of counter-measures 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 post-death transmission. We show that underlying disease parameters cannot be inferred with confidence from early-stage 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 post-death transmission tend to underestimate the basic reproductive number -- thus, given the observed rate of epidemic growth, larger amounts of post-death transmission imply larger reproductive numbers. From a control perspective, we explain how improvements in reducing post-death transmission of Ebola may reduce the overall epidemic spread and scope substantially. Increased attention to the proportion of post-death 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 gamma-ray emission, together with the radio-to-X-ray emission, during the non-flaring (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 multi-instrument effort included the most sensitive VHE gamma-ray instruments in the northern hemisphere, namely the imaging atmospheric Cherenkov telescopes MAGIC and VERITAS, as well as Swift, RXTE, the F-GAMMA, GASP-WEBT, 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 one-zone synchrotron self-Compton 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 non-flaring states. We report for the first time a tentative X-ray-to-VHE correlation during a low VHE activity. Although marginally significant, this positive correlation between X-ray and VHE, which has been reported many times during flaring activity, suggests that the mechanisms that dominate the X-ray/VHE emission during non-flaring-activity are not substantially different from those that are responsible for the emission during flaring activity.
• Error-analysis and comparison to analytical models of numerical waveforms produced by the NRAR Collaboration(1307.5307)

Dec. 11, 2013 gr-qc
The Numerical-Relativity-Analytical-Relativity (NRAR) collaboration is a joint effort between members of the numerical relativity, analytical relativity and gravitational-wave data analysis communities. The goal of the NRAR collaboration is to produce numerical-relativity simulations of compact binaries and use them to develop accurate analytical templates for the LIGO/Virgo Collaboration to use in detecting gravitational-wave 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 non-spinning 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 previously-calibrated, non-precessing analytical waveforms, notably the effective-one-body (EOB) and phenomenological template families, to the newly-produced numerical waveforms. We find that when the binary's total mass is ~100-200 solar masses, current EOB and phenomenological models of spinning, non-precessing binary waveforms have overlaps above 99% (for advanced LIGO) with all of the non-precessing-binary 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 non-spinning EOB waveforms previously calibrated to five non-spinning 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.
• Generalized Deam-Edwards Approach to the Statistical Mechanics of Randomly Crosslinked Systems(1305.0718)

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 Deam-Edwards, discuss the saddle-point order parameters and its physical significance. Finally we also discuss the connection between saddle-point approximation of vulcanization theory, and the classical theory of rubber elasticity as well as the neo-classical theory of nematic elastomers.
• Tests of Lorentz Invariance Violation with Gamma Rays to probe Quantum Gravity(1305.0264)

May 1, 2013 astro-ph.HE
We discuss the prospects of doing tests of Lorentz invariance with gamma-rays observed with present and future ground based gamma-ray observatories.
• Superfast Tikhonov Regularization of Toeplitz Systems(1302.6945)

Feb. 28, 2013 math.NA
Toeplitz-structured 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 "extension-and-transformation" technique, our algorithm translates a Tikhonov-regularized 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.
• Symmetry constraints on temporal order in measurement-based quantum computation(1210.0620)

Oct. 2, 2012 quant-ph
We discuss the interdependence of resource state, measurement setting and temporal order in measurement-based 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.
• Segmentation of Facial Expressions Using Semi-Definite Programming and Generalized Principal Component Analysis(0906.1763)

June 10, 2009 cs.CV
In this paper, we use semi-definite programming and generalized principal component analysis (GPCA) to distinguish between two or more different facial expressions. In the first step, semi-definite 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.
• Fano Resonance Between Mie and Bragg Scattering in Photonic Crystals(0902.4053)

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 non-uniformity of a-SiO2 opal spheres. Proposed theoretical "quasi-3D" model produces results in excellent agreement with the experimental data.
• Sticky obstacles to intramolecular energy flow(0708.1893)

Aug. 14, 2007 nlin.CD, physics.chem-ph
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 lower-dimensional tori; their bifurcations and unstable manifolds govern the transition mechanisms.
• Ultraquantum Dynamics(hep-th/0304031)

April 3, 2003 hep-th
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 order-parameter in the quantum limit where hbar' and hbar'' approach zero. One physical consequence is a discrete spectrum for canonical variable and space-time coordinates. Another is an interconvention of time and energy accompnying space-time meltdown (disorder), with a fundamental conversion factor of some kilograms of energy per second.
• Instability driven fragmentation of nanoscale fractal islands(cond-mat/0204327)

April 15, 2002 cond-mat.mtrl-sci
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 length-to-width ratio of the fractal arms.
• Magnetic Field Effect in a Two-dimensional Array of Short Josephson Junctions(chao-dyn/9504017)

May 3, 1995 nlin.CD, chao-dyn, cond-mat
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.