We introduce a pathwise integration for Volterra processes driven by L\'evy
noise or martingale noise. These processes are widely used in applications to
turbulence, signal processes, biology, and in environmental finance. Indeed
they constitute a very flexible class of models, which include fractional
Brownian and L\'evy motions and it is part of the so-called ambit fields. A
pathwise integration with respect of such Volterra processes aims at producing
a framework where modelling is easily understandable from an information
perspective. The techniques used are based on fractional calculus and in this
there is a bridging of the stochastic and deterministic techniques. The present
paper aims at setting the basis for a framework in which further computational
rules can be devised. Our results are general in the choice of driving noise.
Additionally we propose some further details in the relevant context
subordinated Wiener processes.
Route & Wolszczan (2016) recently detected five radio bursts from the T6
dwarf WISEP J112254.73+255021.5 and used the timing of these events to propose
that this object rotates with an ultra-short period of ~17.3 minutes. We
conducted follow-up observations with the Very Large Array and Gemini-North but
found no evidence for this periodicity. We do, however, observe variable,
highly circularly polarized radio emission possibly with a period of 116
minutes, although our observation lasted only 162 minutes and so more data are
needed to confirm it. Our proposed periodicity is typical of other radio-active
ultracool dwarfs. The handedness of the circular polarization alternates with
time and there is no evidence for any unpolarized emission component, the first
time such a phenomenology has been observed in radio studies of very low-mass
stars and brown dwarfs. We suggest that the object's magnetic dipole axis may
be highly misaligned relative to its rotation axis.
Authors' comments: 9 pages, 5 figures, AAS submitted
We address Gillis' recent criticism [arXiv:1506.05795] of a series of papers
(by different combinations of the present authors) on formulations of Bell's
theorem. Those papers intended to address an unfortunate gap of communication
between two broad camps in the quantum foundations community that we identify
as "operationalists" and "realists". Here, we once again urge the readers to
approach the question from an unbiased standpoint, and explain that Gillis'
criticism draws too heavily on the philosophical inclinations of one side of
that debate -- the realist camp. As part of that explanation we discuss
intuition versus proof, look again at Bell's formalizations of locality, and
correct misstatements by Gillis of our views, and those of Bell and Einstein.
Authors' comments: 10 pages. To be published in International Journal of Quantum
Foundations
Derandomization of Chernoff bound with union bound is already proven in many
papers. We here give another explicit version of it that obtains a construction
of size that is arbitrary close to the probabilistic nonconstructive size.
We apply this to give a new simple polynomial time constructions of almost
$k$-wise independent sets. We also give almost tight lower bounds for the size
of $k$-wise independent sets.
We determine the 22$\mu$m luminosity evolution and luminosity function for
quasars from a data set of over 20,000 objects obtained by combining
flux-limited Sloan Digital Sky Survey optical and Wide field Infrared Survey
Explorer mid-infrared data. We apply methods developed in previous works to
access the intrinsic population distributions non-parametrically, taking into
account the truncations and correlations inherent in the data. We find that the
population of quasars exhibits positive luminosity evolution with redshift in
the mid-infrared, but with considerably less mid-infrared evolution than in the
optical or radio bands. With the luminosity evolutions accounted for, we
determine the density evolution and local mid-infrared luminosity function. The
latter displays a sharp flattening at local luminosities below $\sim 10^{31}$
erg sec$^{-1}$ Hz$^{-1}$, which has been reported previously at 15 $\mu$m for
AGN classified as both type-1 and type-2. We calculate the integrated total
emission from quasars at 22 $\mu$m and find it to be a small fraction of both
the cosmic infrared background light and the integrated emission from all
sources at this wavelength.
Authors' comments: 13 pages, 11 figures, published in ApJ. arXiv admin note: text
overlap with arXiv:1207.3396
We have monitored photometrically the Y0 brown dwarf WISEP
J173835.52+273258.9 (W1738) at both near- and mid-infrared wavelengths. This ~1
Gyr-old 400K dwarf is at a distance of 8pc and has a mass around 5 M_Jupiter.
We observed W1738 using two near-infrared filters at lambda~1um, Y and J, on
Gemini observatory, and two mid-infrared filters at lambda~4um, [3.6] and
[4.5], on the Spitzer observatory. Twenty-four hours were spent on the source
by Spitzer on each of June 30 and October 30 2013 UT. Between these
observations, around 5 hours were spent on the source by Gemini on each of July
17 and August 23 2013 UT. The mid-infrared light curves show significant
evolution between the two observations separated by four months. We find that a
double sinusoid can be fit to the [4.5] data, where one sinusoid has a period
of 6.0 +/- 0.1 hours and the other a period of 3.0 +/- 0.1 hours. The
near-infrared observations suggest variability with a ~3.0 hour period,
although only at a <~2 sigma confidence level. We interpret our results as
showing that the Y dwarf has a 6.0 +/- 0.1 hour rotation period, with one or
more large-scale surface features being the source of variability. The
peak-to-peak amplitude of the light curve at [4.5] is 3%. The amplitude of the
near-infrared variability, if real, may be as high as 5 to 30%. Intriguingly,
this size of variability and the wavelength dependence can be reproduced by
atmospheric models that include patchy KCl and Na_2S clouds and associated
small changes in surface temperature. The small number of large features, and
the timescale for evolution of the features, is very similar to what is seen in
the atmospheres of the solar system gas giants.
Authors' comments: Accepted by ApJ July 26 2016. Twenty-six pages include 8 Figures and
5 Tables
We present an analysis of the mid-infrared WISE sources seen within the
equatorial GAMA G12 field, located in the North Galactic Cap. Our motivation is
to study and characterize the behavior of WISE source populations in
anticipation of the deep multi-wavelength surveys that will define the next
decade, with the principal science goal of mapping the 3D large scale
structures and determining the global physical attributes of the host galaxies.
In combination with cosmological redshifts, we identify galaxies from their
WISE W1 3.4um resolved emission, and by performing a star-galaxy separation
using apparent magnitude, colors and statistical modeling of star-counts. The
resultant galaxy catalog has ~590,000 sources in 60 deg^2, reaching a W1
5-sigma depth of 31 uJy. At the faint end, where redshifts are not available,
we employ a luminosity function analysis to show that approximately 27% of all
WISE extragalactic sources to a limit of 17.5 mag (31 uJy) are at high
redshift, z > 1. The spatial distribution is investigated using two-point
correlation functions and a 3D source density characterization at 5 Mpc and 20
Mpc scales. For angular distributions, we find brighter and more massive
sources are strongly clustered relative to fainter and lower mass source;
likewise, based on WISE colors, spheroidal galaxies have the strongest
clustering, while late-type disk galaxies have the lowest clustering
amplitudes. Along the radial direction, the strongest clustering is in the
largest redshift shell, while the weakest is in the nearest redshift shell,
consistent with the stellar mass and morphological type dependency results. In
three dimensions, we find a number of distinct groupings, often bridged by
filaments and super-structures. Using special visualization tools, we map these
structures, exploring how clustering may play a role with stellar mass and
galaxy type.
Authors' comments: accepted for publication in the ApJ (Dec 30, 2016); animations and
additional graphics: http://www.ast.uct.ac.za/~jarrett/G12/
We cross-match the two currently largest all-sky photometric catalogs,
mid-infrared WISE and SuperCOSMOS scans of UKST/POSS-II photographic plates, to
obtain a new galaxy sample that covers 3pi steradians. In order to characterize
and purify the extragalactic dataset, we use external GAMA and SDSS
spectroscopic information to define quasar and star loci in multicolor space,
aiding the removal of contamination from our extended-source catalog. After
appropriate data cleaning we obtain a deep wide-angle galaxy sample that is
approximately 95% pure and 90% complete at high Galactic latitudes. The catalog
contains close to 20 million galaxies over almost 70% of the sky, outside the
Zone of Avoidance and other confused regions, with a mean surface density of
over 650 sources per square degree. Using multiwavelength information from two
optical and two mid-IR photometric bands, we derive photometric redshifts for
all the galaxies in the catalog, using the ANNz framework trained on the final
GAMA-II spectroscopic data. Our sample has a median redshift of z_{med} = 0.2
but with a broad dN/dz reaching up to z>0.4. The photometric redshifts have a
mean bias of |delta_z|~10^{-3}, normalized scatter of sigma_z = 0.033 and less
than 3% outliers beyond 3sigma_z. Comparison with external datasets shows no
significant variation of photo-z quality with sky position. Together with the
overall statistics, we also provide a more detailed analysis of photometric
redshift accuracy as a function of magnitudes and colors. The final catalog is
appropriate for `all-sky' 3D cosmology to unprecedented depths, in particular
through cross-correlations with other large-area surveys. It should also be
useful for source pre-selection and identification in forthcoming surveys such
as TAIPAN or WALLABY.
Authors' comments: Data available from http://ssa.roe.ac.uk/WISExSCOS
The two currently largest all-sky photometric datasets, WISE and SuperCOSMOS,
were cross-matched by Bilicki et al. (2016) (B16) to construct a novel
photometric redshift catalogue on 70% of the sky. Galaxies were therein
separated from stars and quasars through colour cuts, which may leave
imperfections because of mixing different source types which overlap in colour
space. The aim of the present work is to identify galaxies in the
WISExSuperCOSMOS catalogue through an alternative approach of machine learning.
This allows us to define more complex separations in the multi-colour space
than possible with simple colour cuts, and should provide more reliable source
classification. For the automatised classification we use the support vector
machines learning algorithm, employing SDSS spectroscopic sources cross-matched
with WISExSuperCOSMOS as the training and verification set. We perform a number
of tests to examine the behaviour of the classifier (completeness, purity and
accuracy) as a function of source apparent magnitude and Galactic latitude. We
then apply the classifier to the full-sky data and analyse the resulting
catalogue of candidate galaxies. We also compare thus produced dataset with the
one presented in B16. The tests indicate very high accuracy, completeness and
purity (>95%) of the classifier at the bright end, deteriorating for the
faintest sources, but still retaining acceptable levels of 85%. No significant
variation of classification quality with Galactic latitude is observed.
Application of the classifier to all-sky WISExSuperCOSMOS data gives 15 million
galaxies after masking problematic areas. The resulting sample is purer than
the one in B16, at a price of lower completeness over the sky. The automatic
classification gives a successful alternative approach to defining a reliable
galaxy sample as compared to colour cuts.
Authors' comments: 12 pages, 15 figures, accepted for publication in A&A. Obtained
catalogue will be included in the public release of the WISExSuperCOSMOS
galaxy catalogue available from http://ssa.roe.ac.uk/WISExSCOS
We present a simultaneous, multi-wavelength campaign targeting the nearby
(7.2 pc) L8/L9 (optical/near-infrared) dwarf WISEP J060738.65+242953.4 in the
mid-infrared, radio, and optical. Spitzer Space Telescope observations show no
variability at the 0.2% level over 10 hours each in the 3.6 and 4.5 micron
bands. Kepler K2 monitoring over 36 days in Campaign 0 rules out stable
periodic signals in the optical with amplitudes great than 1.5% and periods
between 1.5 hours and 2 days. Non-simultaneous Gemini optical spectroscopy
detects lithium, constraining this L dwarf to be less than ~2 Gyr old, but no
Balmer emission is observed. The low measured projected rotation velocity (v
sin i < 6 km/s) and lack of variability are very unusual compared to other
brown dwarfs, and we argue that this substellar object is likely viewed
pole-on. We detect quiescent (non-bursting) radio emission with the VLA.
Amongst radio detected L and T dwarfs, it has the lowest observed L_nu and the
lowest v sin i. We discuss the implications of a pole-on detection for various
proposed radio emission scenarios.
Authors' comments: Accepted to A.J
H{\sc i} absorption studies of active galaxies enable us to probe their
circumnuclear regions and the general interstellar medium, and study the supply
of gas which may trigger the nuclear activity. In this paper, we investigate
the detection rate of H{\sc i} absorption on the nature of radio galaxies based
on their emission-line spectra, nature of the host galaxies based on the
\textit{WISE} colours and their radio structure, which may help understand the
different accretion modes. We find significant difference in distributions of
W2$-$W3 colour for sources with H{\sc i} absorption detections and
non-detections. We report a high detection rate of H{\sc i} absorption in the
galaxies with \textit{WISE} infrared colours W2$-$W3 $>$ 2, which is typical of
gas-rich systems, along with a compact radio structure. The H{\sc i} detection
rate for low-excitation radio galaxies (LERGs) with W2$-$W3 $>$ 2 and compact
radio structure is high (70.6$\pm$20.4 \%). In HERGs, compact radio structure
in the nuclear or circumnuclear region could give rise to absorption by gas in
the dusty torus in addition to gas in the interstellar medium. However, higher
specific star formation rate (sSFR) for the LERGs with W2$-$W3 $>$ 2 suggests
that H{\sc i} absorption may be largely due to star-forming gas in their hosts.
LERGs with extended radio structure tend to have significantly lower values of
W2$-$W3 compared to those with compact structure. Extended radio sources and
those with W2$-$W3 $<$ 2 have low H{\sc i} detection rates.
Authors' comments: 10 pages, 7 figures, this updated version is accepted for publication
in MNRAS
This technical report revisits the analysis of family-wise error rates in
statistical parametric mapping - using random field theory - reported in
(Eklund et al., 2015). Contrary to the understandable spin that these sorts of
analyses attract, a review of their results suggests that they endorse the use
of parametric assumptions - and random field theory - in the analysis of
functional neuroimaging data. We briefly rehearse the advantages parametric
analyses offer over nonparametric alternatives and then unpack the implications
of (Eklund et al., 2015) for parametric procedures.
In this paper, we tackle the problem of temporally consistent boundary
detection and hierarchical segmentation in videos. While finding the best
high-level reasoning of region assignments in videos is the focus of much
recent research, temporal consistency in boundary detection has so far only
rarely been tackled. We argue that temporally consistent boundaries are a key
component to temporally consistent region assignment. The proposed method is
based on the point-wise mutual information (PMI) of spatio-temporal voxels.
Temporal consistency is established by an evaluation of PMI-based point
affinities in the spectral domain over space and time. Thus, the proposed
method is independent of any optical flow computation or previously learned
motion models. The proposed low-level video segmentation method outperforms the
learning-based state of the art in terms of standard region metrics.
There is mounting evidence that supermassive black holes form and grow in
bulgeless galaxies. However, a robust determination of the fraction of AGNs in
bulgeless galaxies, an important constraint to models of supermassive black
hole seed formation and merger-free models of AGN fueling, is unknown, since
optical studies have been shown to be incomplete for low mass AGNs. In a recent
study using the Wide-field Infrared Survey Explorer, we discovered hundreds of
bulgeless galaxies that display mid-infrared signatures of extremely hot dust
suggestive of powerful accreting massive black holes, despite having no
signatures of black hole activity at optical wavelengths. Here we report X-ray
follow-up observations of J122434.66+555522.3, a nearby (z=0.052) isolated
bulgeless galaxy that contains an unresolved X-ray source detected at the 3
sigma level by XMM-Newton with an observed luminosity uncorrected for intrinsic
absorption of L{2-10~keV}=1.1+/-0.4 10^40 ergs/s. Ground-based near-infrared
spectroscopy with the Large Binocular Telescope together with multiwavelength
observations from ultraviolet to millimeter wavelengths together suggest that
J1224+5555 harbors a highly absorbed AGN with an intrinsic absorption of ~NH
>10^24 cm^-2. The hard X-ray luminosity of the putative AGN corrected for
absorption is L{2-10~keV}~3x10^42 ergs/s, which, depending on the bolometric
correction factor, corresponds to a bolometric luminosity of the AGN of
6x10^43ergs/s - 3x10^44 erg/s, and a lower mass limit for the black hole of
MBH~2x10^6 Msun, based on the Eddington limit.
Authors' comments: Accepted to ApJ, 14 pages, 7 figures
The recent past years have seen a noticeable increase of interest in the
correlation analysis of electrohysterographic (EHG) signals in the perspective
of improving the pregnancy monitoring. Here we propose a new approach based on
the functional connectivity between multichannel (4x4 matrix) EHG signals
recorded from the women abdomen. The proposed pipeline includes i) the
computation of the statistical couplings between the multichannel EHG signals,
ii) the characterization of the connectivity matrices, computed by using the
imaginary part of the coherence, based on the graph-theory analysis and iii)
the use of these measures for pregnancy monitoring. The method was evaluated on
a dataset of EHGs, in order to track the correlation between EHGs collected by
each electrode of the matrix (called node-wise analysis) and follow their
evolution along weeks before labor. Results showed that the strength of each
node significantly increases from pregnancy to labor. Electrodes located on the
median vertical axis of the uterus seemed to be the more discriminant. We
speculate that the network-based analysis can be a very promising tool to
improve pregnancy monitoring.
Authors' comments: 4 pages, 3 figures, accepted in the IEEE EMBC conferance
The best arm identification problem (BEST-1-ARM) is the most basic pure
exploration problem in stochastic multi-armed bandits. The problem has a long
history and attracted significant attention for the last decade. However, we do
not yet have a complete understanding of the optimal sample complexity of the
problem: The state-of-the-art algorithms achieve a sample complexity of
$O(\sum_{i=2}^{n} \Delta_{i}^{-2}(\ln\delta^{-1} + \ln\ln\Delta_i^{-1}))$
($\Delta_{i}$ is the difference between the largest mean and the $i^{th}$
mean), while the best known lower bound is $\Omega(\sum_{i=2}^{n}
\Delta_{i}^{-2}\ln\delta^{-1})$ for general instances and $\Omega(\Delta^{-2}
\ln\ln \Delta^{-1})$ for the two-arm instances. We propose to study the
instance-wise optimality for the BEST-1-ARM problem. Previous work has proved
that it is impossible to have an instance optimal algorithm for the 2-arm
problem. However, we conjecture that modulo the additive term
$\Omega(\Delta_2^{-2} \ln\ln \Delta_2^{-1})$ (which is an upper bound and worst
case lower bound for the 2-arm problem), there is an instance optimal algorithm
for BEST-1-ARM. Moreover, we introduce a new quantity, called the gap entropy
for a best-arm problem instance, and conjecture that it is the instance-wise
lower bound. Hence, resolving this conjecture would provide a final answer to
the old and basic problem.
Authors' comments: To appear in COLT 2016 Open Problems
Most of machine learning approaches have stemmed from the application of
minimizing the mean squared distance principle, based on the computationally
efficient quadratic optimization methods. However, when faced with
high-dimensional and noisy data, the quadratic error functionals demonstrated
many weaknesses including high sensitivity to contaminating factors and
dimensionality curse. Therefore, a lot of recent applications in machine
learning exploited properties of non-quadratic error functionals based on $L_1$
norm or even sub-linear potentials corresponding to quasinorms $L_p$ ($0<p<1$).
The back side of these approaches is increase in computational cost for
optimization. Till so far, no approaches have been suggested to deal with {\it
arbitrary} error functionals, in a flexible and computationally efficient
framework. In this paper, we develop a theory and basic universal data
approximation algorithms ($k$-means, principal components, principal manifolds
and graphs, regularized and sparse regression), based on piece-wise quadratic
error potentials of subquadratic growth (PQSQ potentials). We develop a new and
universal framework to minimize {\it arbitrary sub-quadratic error potentials}
using an algorithm with guaranteed fast convergence to the local or global
error minimum. The theory of PQSQ potentials is based on the notion of the cone
of minorant functions, and represents a natural approximation formalism based
on the application of min-plus algebra. The approach can be applied in most of
existing machine learning methods, including methods of data approximation and
regularized and sparse regression, leading to the improvement in the
computational cost/accuracy trade-off. We demonstrate that on synthetic and
real-life datasets PQSQ-based machine learning methods achieve orders of
magnitude faster computational performance than the corresponding
state-of-the-art methods.
Authors' comments: Edited and extended version with algortihms of regularized regression
We report the detection of a significant infrared variability of the nearest
tidal disruption event (TDE) ASASSN-14li using Wide-field Infrared Survey
Explorer} and newly released Near-Earth Object WISE Reactivation} data. In
comparison with the quiescent state, the infrared flux is brightened by 0.12
and 0.16 magnitude in the W1 ($3.4\mu$m) and W2 ($4.6\mu$m) bands at 36 days
after the optical discovery (or $\sim110$ days after the peak disruption date).
The flux excess is still detectable $\sim170$ more days later. Assuming that
the flare-like infrared emission is from the dust around the black hole, its
blackbody temperature is estimated to be $\sim2.1\times10^3$~K, slightly higher
than the dust sublimation temperature, indicating that the dust is likely
located close to the dust sublimation radius. The equilibrium between the
heating and radiation of the dust claims a bolometric luminosity of
$\sim10^{43}-10^{45}$~\lum, comparable with the observed peak luminosity. This
result has for the first time confirmed the detection of infrared emission from
the dust echoes of TDEs.
Authors' comments: published by ApJL
The Stokes system with constant viscosity can be cast into different
formulations by exploiting the incompressibility constraint. For instance the
strain in the weak formulation can be replaced by the gradient to decouple the
velocity components in the different coordinate directions. Thus the
discretization of the simplified problem leads to fewer nonzero entries in the
stiffness matrix. This is of particular interest in large scale simulations
where a reduced memory bandwidth requirement can help to significantly
accelerate the computations. In the case of a piecewise constant viscosity, as
it typically arises in multi-phase flows, or when the boundary conditions
involve traction, the situation is more complex, and one has to treat the cross
derivatives in the original Stokes system with care. A naive application of the
standard vectorial Laplacian results in a physically incorrect solution, while
formulations based on the strain increase the computational effort everywhere,
even when the inconsistencies arise only from an incorrect treatment in a small
fraction of the computational domain. Here we propose a new approach that is
consistent with the strain-based formulation and preserves the decoupling
advantages of the gradient-based formulation in isoviscous subdomains. The
modification is equivalent to locally changing the discretization stencils,
hence the more expensive discretization is restricted to a lower dimensional
interface, making the additional computational cost asymptotically negligible.
We demonstrate the consistency and convergence properties of the method and
show that in a massively parallel setup, the multigrid solution of the
resulting discrete systems is faster than for the classical strain-based
formulation. Moreover, we give an application example which is inspired by
geophysical research.
The astronomy community has at its disposal a large back catalog of public
spectroscopic galaxy redshift surveys that can be used for the measurement of
luminosity functions. Utilizing the back catalog with new photometric surveys
to maximum efficiency requires modeling the color selection bias imposed on
selection of target galaxies by flux limits at multiple wavelengths. The
likelihood derived herein can address, in principle, all possible color
selection biases through the use of a generalization of the luminosity
function, $\Phi(L)$, over the space of all spectra: the spectro-luminosity
functional, $\Psi[L_\nu]$. It is, therefore, the first estimator capable of
simultaneously analyzing multiple redshift surveys in a consistent way. We also
propose a new way of parametrizing the evolution of the classic Shechter
function parameters, $L_\star$ and $\phi_\star$, that improves both the
physical realism and statistical performance of the model. The techniques
derived in this work will be used in an upcoming paper to measure the
luminosity function of galaxies at the rest frame wavelength of
$2.4\operatorname{\mu m}$ using the Widefield Infrared Survey Explorer (WISE).
Authors' comments: Accepted to AJ 14 pages, 6 figures
