The paper deals with point-wise estimates for the heat kernel of a nonlocal
convolution type operator with a kernel that decays at least exponentially at
infinity. It is shown that the large time behaviour of the heat kernel depends
essentially on whether $|x|\ll t$, or $|x|\sim t$, or $|x|\gg t$. We obtain
sharp point-wise upper bounds for the heat kernel in all these regions. In the
first region the nonlocal heat kernel behaves like the classical one, while in
the other regions we observe an essential difference.
We present an analysis on the potential prevalence of mid infrared excesses
in stars with and without planetary companions. Based on an extended database
of stars detected with the ${\it WISE}$ satellite, we studied two stellar
samples: one with 236 planet hosts and another with 986 objects for which
planets have been searched but not found. We determined the presence of an
excess over the photosphere by comparing the observed flux ratio at 22 $\mu$m
and 12 $\mu$m ($f_{22}/f_{12}$) with the corresponding synthetic value, derived
from results of classical model photospheres. We found a detection rate of
0.85$\%$ at 22 $\mu$m (2 excesses) in the sample of stars with planets and
0.1$\%$ (1 detection) for the stars without planets. The difference of the
detection rate between the two samples is not statistically significant, a
result that is independent of the different approaches found in the literature
to define an excess in the wavelength range covered by ${\it WISE}$
observations. As an additional result, we found that the ${\it WISE}$ fluxes
required a normalisation procedure to make them compatible with synthetic data,
probably pointing out a revision of the ${\it WISE}$ data calibration.
Authors' comments: 10 pages, 6 figures, 3 tables. Accepted for publication in MNRAS
Wide-angle photometric surveys of previously uncharted sky areas or
wavelength regimes will always bring in unexpected sources whose existence and
properties cannot be easily predicted from earlier observations: novelties or
even anomalies. Such objects can be efficiently sought for with novelty
detection algorithms. Here we present an application of such a method, called
one-class support vector machines (OCSVM), to search for anomalous patterns
among sources preselected from the mid-infrared AllWISE catalogue covering the
whole sky. To create a model of expected data we train the algorithm on a set
of objects with spectroscopic identifications from the SDSS DR13 database,
present also in AllWISE. OCSVM detects as anomalous those sources whose
patterns - WISE photometric measurements in this case - are inconsistent with
the model. Among the detected anomalies we find artefacts, such as objects with
spurious photometry due to blending, but most importantly also real sources of
genuine astrophysical interest. Among the latter, OCSVM has identified a sample
of heavily reddened AGN/quasar candidates distributed uniformly over the sky
and in a large part absent from other WISE-based AGN catalogues. It also
allowed us to find a specific group of sources of mixed types, mostly stars and
compact galaxies. By combining the semi-supervised OCSVM algorithm with
standard classification methods it will be possible to improve the latter by
accounting for sources which are not present in the training sample but are
otherwise well-represented in the target set. Anomaly detection adds
flexibility to automated source separation procedures and helps verify the
reliability and representativeness of the training samples. It should be thus
considered as an essential step in supervised classification schemes to ensure
completeness and purity of produced catalogues.
Authors' comments: 14 pages, 15 figures
We present the first results of an ongoing variability monitoring program of
active galactic nuclei (AGNs) using the 46\,cm telescope of the Wise
observatory in Israel. The telescope has a field of view of $1.25^{\circ}
\times 0.84^{\circ}$ and is specially equipped with five narrow band filters at
4300, 5200, 5700, 6200 and 7000\,\AA\, in order to perform photometric
reverberation mapping studies of the central engine of AGNs. The program aims
to observe a sample of 27 AGNs (V $<17$ mag) selected according to tentative
continuum and line time delay measurements obtained in previous works. We
describe the autonomous operation of the telescope together with the fully
automatic pipeline used to achieve high-performance unassisted observations,
data reduction, and light curves extraction using different photometric
methods. The science verification data presented here, demonstrates the
performance of the monitoring program in particular for efficiently photometric
reverberation mapping of AGNs with additional capabilities to carry out
complementary studies of other transient and variable phenomena such as
variable stars studies.
Authors' comments: 14 pages, 13 figures, accepted for publication in PASP
The interplay of rotation and manifested magnetic activity on ultracool
dwarfs (UCDs) is of key importance for gathering clues as to the operation of
the dynamos within these objects. A number of magnetized UCDs host kG-strength
magnetic fields. It was recently reported that the L8 dwarf WISEP
J060738.65+242953.4 is a radio-emitting UCD that is likely observed pole-on,
due to its lack of photometric variability and narrow spectral lines. Follow-up
radio observations at Arecibo Observatory, together with careful analysis of
previously published details, however, suggest that the scientific and
statistical significance of the radio and spectroscopic data have been
overstated. If the UCD is observed along its aligned spin/magnetic axis, the
absence of observed H$\alpha$ activity may present challenges to the auroral
model of UCD magnetism, although short or long-term cyclic magnetic activity
may explain this behavior. Monte Carlo simulations presented here suggest that
the source probably rotates with v sin i=6-12 km s$^{-1}$, indicating that its
inclination angle and rotational velocity are unexceptional and that its
angular momentum has evolved as expected for brown dwarfs observed in ~1
Myr-old clusters. The discovery and verification of the most rapidly and
slowest rotating brown dwarfs place valuable constraints on the angular
momentum evolution and magnetic activity history of these objects.
Authors' comments: Published in ApJ; 7 pages, 1 figure
The two-component broad-band spectral energy distributions of blazars were
suggested to form a sequence in which (1) the peak frequency of the low-energy
(synchrotron) component nu_syn is anticorrelated with the synchrotron
luminosity L_syn, and (2) the luminosity ratio of the high-energy (inverse
Compton) to synchrotron components q = L_IC / L_syn (Compton dominance)
increases with L_syn from the BL Lac objects (BL Lacs) to the flat spectrum
radio quasars (FSRQs). The Compton dominance parameter is an important probe of
plasma magnetisation in the blazar zones within relativistic jets. We
investigate a sample of blazars detected by WISE in the mid-infrared (MIR) band
and by Fermi/LAT in the GeV gamma-ray band, with the focus on the distribution
of luminosities and photon indices. Our findings are the following: (1) the MIR
photon index Gamma_W12 is a useful probe of the blazar sequence, with the
exception of low-luminosity BL Lacs that are most likely contaminated by their
host galaxies (L_W1 ~ 10^44 erg/s and Gamma_W12 < 1); (2) Gamma_W12 is
correlated with the gamma-ray photon index Gamma_1-100GeV, with the MIR
luminosity L_W1, and with the Fermi/WISE Compton dominance q_FW = L_1GeV/L_W1;
(3) a clean separation between FSRQs and BL Lacs can be seen in the parameter
space of Gamma_W12 and q_FW; (4) the observed distribution of MIR luminosity
L_W1 vs. Compton dominance q_FW for the entire sample of blazars can be modeled
as a sequence of lepto-magnetic jet powers in the range (log_10 P_eB) in
[42:45] with the preference for sub-equipartition magnetic fields P_B/P_e in
[0.05:1], assuming fixed bulk Lorentz factor Gamma_j = 15, fixed jet opening
angle Gamma_j Theta_j = 0.3, fixed radiative efficiency of jet electrons
epsilon_em = 50%, and that external radiation luminosity scales like L_ext
propto (P_eB)^(1.6) (parameter degeneracies are discussed).
Authors' comments: 9 pages, 5 figures, submitted to A&A
Context. Mars-crossing asteroids (MCs) are a dynamically unstable group
between the main belt and the near-Earth populations. Characterising the
physical properties of a large sample of MCs can help to understand the
original sources of many near-Earth asteroids, some of which may produce
meteorites on Earth. Aims. Our aim is to provide diameters and albedos of MCs
with available WISE/NEOWISE data. Methods. We used the near-Earth asteroid
thermal model to find the best-fitting values of equivalent diameter and,
whenever possible, the infrared beaming parameter. With the diameter and
tabulated asteroid absolute magnitudes we also computed the visible geometric
albedos. Results. We determined the diameters and beaming parameters of 404
objects observed during the fully cryogenic phase of the WISE mission, most of
which have not been published elsewhere. We also obtained 1572 diameters from
data from the 3-Band and posterior non-cryogenic phases using a default value
of beaming parameter. The average beaming parameter is 1.2 +/- 0.2 for objects
smaller than 10 km, which constitute most of our sample. This is higher than
the typical value of 1.0 found for the whole main belt and is possibly related
to the fact that WISE is able to observe many more small objects at shorter
heliocentric distances, i.e. at higher phase angles. We argue that this is a
better default value for modelling Mars-crossing asteroids from the
WISE/NEOWISE catalogue and discuss the effects of this choice on the diameter
and albedo distributions. We find a double-peaked distribution for the visible
geometric albedos, which is expected since this population is compositionally
diverse and includes objects in the major spectral complexes. However, the
distribution of beaming parameters is homogeneous for both low- and high-albedo
objects.
Authors' comments: 8 pages, 6 figures, accepted for publication in Astronomy &
Astrophysics
Many real datasets contain values missing not at random (MNAR). In this
scenario, investigators often perform list-wise deletion, or delete samples
with any missing values, before applying causal discovery algorithms. List-wise
deletion is a sound and general strategy when paired with algorithms such as
FCI and RFCI, but the deletion procedure also eliminates otherwise good samples
that contain only a few missing values. In this report, we show that we can
more efficiently utilize the observed values with test-wise deletion while
still maintaining algorithmic soundness. Here, test-wise deletion refers to the
process of list-wise deleting samples only among the variables required for
each conditional independence (CI) test used in constraint-based searches.
Test-wise deletion therefore often saves more samples than list-wise deletion
for each CI test, especially when we have a sparse underlying graph. Our
theoretical results show that test-wise deletion is sound under the justifiable
assumption that none of the missingness mechanisms causally affect each other
in the underlying causal graph. We also find that FCI and RFCI with test-wise
deletion outperform their list-wise deletion and imputation counterparts on
average when MNAR holds in both synthetic and real data.
How to develop slim and accurate deep neural networks has become crucial for
real- world applications, especially for those employed in embedded systems.
Though previous work along this research line has shown some promising results,
most existing methods either fail to significantly compress a well-trained deep
network or require a heavy retraining process for the pruned deep network to
re-boost its prediction performance. In this paper, we propose a new layer-wise
pruning method for deep neural networks. In our proposed method, parameters of
each individual layer are pruned independently based on second order
derivatives of a layer-wise error function with respect to the corresponding
parameters. We prove that the final prediction performance drop after pruning
is bounded by a linear combination of the reconstructed errors caused at each
layer. Therefore, there is a guarantee that one only needs to perform a light
retraining process on the pruned network to resume its original prediction
performance. We conduct extensive experiments on benchmark datasets to
demonstrate the effectiveness of our pruning method compared with several
state-of-the-art baseline methods.
We have reprocessed over 100 terabytes of single-exposure WISE/NEOWISE images
to create the deepest ever full-sky maps at 3-5 microns. We incorporate all
publicly available W1 and W2 imaging - a total of ~8 million exposures in each
band - from ~37 months of observations spanning 2010 January to 2015 December.
Our coadds preserve the native WISE resolution and feature depth of coverage ~3
times greater than that of the AllWISE Atlas stacks. Our coadds are designed to
enable deep forced photometry, in particular for the Dark Energy Camera Legacy
Survey (DECaLS) and Mayall z-Band Legacy Survey (MzLS), both of which are being
used to select targets for the Dark Energy Spectroscopic Instrument (DESI). We
describe newly introduced processing steps aimed at leveraging added redundancy
to remove artifacts, with the intent of facilitating uniform target selection
and searches for rare/exotic objects (e.g. high-redshift quasars and distant
galaxy clusters). Forced photometry depths achieved with these coadds extend
0.56 (0.46) magnitudes deeper in W1 (W2) than is possible with only
pre-hibernation WISE imaging.
Authors' comments: data release available at http://unwise.me
Recent work has found evidence for a difference in the bias and dark matter
halo masses of WISE-selected obscured and unobscured quasars, implying a
distinction between these populations beyond random line-of-sight effects.
However, the significance of this difference in the most up-to-date
measurements is relatively weak, at $\sim$2$\sigma$ for individual measurements
but bolstered by agreement from different techniques, including angular
clustering and cross-correlations with cosmic microwave background (CMB)
lensing maps. Here, we expand the footprint of previous work, aiming to improve
the precision of both methods. In this larger area we correct for position
dependent selection effects, in particular fluctuations of the WISE-selected
quasar density as a function of Galactic latitude. We also measure the
cross-correlation of the obscured and unobscured samples and confirm that they
are well-matched in redshift, both centred at $z=1$. Combined with very similar
detection fractions and magnitude distributions in the long-wavelength WISE
bands, this redshift match strongly supports the fact that IR selection
identifies obscured and unobscured quasars of similar bolometric luminosity.
Finally, we perform cross-correlations with confirmed spectroscopic quasars,
again confirming the results from other methods --- obscured quasars reside in
haloes a factor of 3 times more massive than unobscured quasars. This
difference is significant at the $\sim$5$\sigma$ level when the measurements
are combined, strong support for the idea that obscuration in at least some
quasars is tied to the larger environment, and may have an evolutionary
component.
Authors' comments: 15 pages, 12 figures, accepted for publication in MNRAS
We investigate extremely luminous dusty galaxies in the environments around
WISE-selected hot dust obscured galaxies (Hot DOGs) and WISE/radio-selected
active galactic nuclei (AGNs) at average redshifts of z = 2.7 and z = 1.7,
respectively. Previous observations have detected overdensities of companion
submillimetre-selected sources around 10 Hot DOGs and 30 WISE/radio AGNs, with
overdensities of ~ 2 - 3 and ~ 5 - 6 , respectively. We find that the space
densities in both samples to be overdense compared to normal star-forming
galaxies and submillimetre galaxies (SMGs) in the SCUBA-2 Cosmology Legacy
Survey (S2CLS). Both samples of companion sources have consistent mid-IR
colours and mid-IR to submm ratios as SMGs. The brighter population around
WISE/radio AGNs could be responsible for the higher overdensity reported. We
also find the star formation rate density (SFRDs) are higher than the field,
but consistent with clusters of dusty galaxies. WISE-selected AGNs appear to be
good signposts for protoclusters at high redshift on arcmin scales. The results
reported here provide an upper limit to the strength of angular clustering
using the two-point correlation function. Monte Carlo simulations show no
angular correlation, which could indicate protoclusters on scales larger than
the SCUBA-2 1.5arcmin scale maps.
Authors' comments: 10 pages, 7 figures, 6 tables, MNRAS accepted
CTR prediction in real-world business is a difficult machine learning problem
with large scale nonlinear sparse data. In this paper, we introduce an
industrial strength solution with model named Large Scale Piece-wise Linear
Model (LS-PLM). We formulate the learning problem with $L_1$ and $L_{2,1}$
regularizers, leading to a non-convex and non-smooth optimization problem.
Then, we propose a novel algorithm to solve it efficiently, based on
directional derivatives and quasi-Newton method. In addition, we design a
distributed system which can run on hundreds of machines parallel and provides
us with the industrial scalability. LS-PLM model can capture nonlinear patterns
from massive sparse data, saving us from heavy feature engineering jobs. Since
2012, LS-PLM has become the main CTR prediction model in Alibaba's online
display advertising system, serving hundreds of millions users every day.
The eShel, an off-the-shelf, fiber-fed echelle spectrograph ($R \approx
10,000$), was installed on the 1m telescope at the Wise observatory in Israel.
We report the installation of the multi-order spectrograph, and describe our
pipeline to extract stellar radial velocity from the obtained spectra. We also
introduce a new algorithm---UNICOR, to remove radial-velocity systematics that
can appear in some of the observed orders. We show that the system performance
is close to the photon-noise limit for exposures with more than $10^7$ counts,
with a precision that can get better than 200 m/s for F--K stars, for which the
eShel spectral response is optimal. This makes the eShel at Wise a useful tool
for studying spectroscopic binaries brighter than $m_V=11$. We demonstrate this
capability with orbital solutions of two binaries from projects being performed
at Wise.
Authors' comments: 30 pages, 12 figures. Accepted for publication in PASP
Hot, Dust-Obscured Galaxies, or "Hot DOGs", are a rare, dusty, hyperluminous
galaxy population discovered by the WISE mission. Predominantly at redshifts
2-3, they include the most luminous known galaxies in the universe. Their high
luminosities likely come from accretion onto highly obscured super massive
black holes (SMBHs). We have conducted a pilot survey to measure the SMBH
masses of five z~2 Hot DOGs via broad H_alpha emission lines, using
Keck/MOSFIRE and Gemini/FLAMINGOS-2. We detect broad H_alpha emission in all
five Hot DOGs. We find substantial corresponding SMBH masses for these Hot DOGs
(~ 10^{9} M_sun), and their derived Eddington ratios are close to unity. These
z~2 Hot DOGs are the most luminous AGNs at given BH masses, suggesting they are
accreting at the maximum rates for their BHs. A similar property is found for
known z~6 quasars. Our results are consistent with scenarios in which Hot DOGs
represent a transitional, high-accretion phase between obscured and unobscured
quasars. Hot DOGs may mark a special evolutionary stage before the red quasar
and optical quasar phases, and they may be present at other cosmic epochs.
Authors' comments: 15 pages, 9 figures. Accepted by ApJ
The stochastic block model is widely used for detecting community structures
in network data. How to test the goodness-of-fit of the model is one of the
fundamental problems and has gained growing interests in recent years. In this
article, we propose a novel goodness-of-fit test based on the maximum entry of
the centered and re-scaled adjacency matrix for the stochastic block model. One
noticeable advantage of the proposed test is that the number of communities can
be allowed to grow linearly with the number of nodes ignoring a logarithmic
factor. We prove that the null distribution of the test statistic converges in
distribution to a Gumbel distribution, and we show that both the number of
communities and the membership vector can be tested via the proposed method.
Further, we show that the proposed test has asymptotic power guarantee against
a class of alternatives. We also demonstrate that the proposed method can be
extended to the degree-corrected stochastic block model. Both simulation
studies and real-world data examples indicate that the proposed method works
well.
Authors' comments: 31 pages
Characterization of classes of switching signals that ensure stability of
switched systems occupies a significant portion of the switched systems
literature. This article collects a multitude of stabilizing switching signals
under an umbrella framework. We achieve this in two steps: Firstly, given a
family of systems, possibly containing unstable dynamics, we propose a new and
general class of stabilizing switching signals. Secondly, we demonstrate that
prior results based on both point-wise and asymptotic characterizations follow
our result. This is the first attempt in the switched systems literature where
these switching signals are unified under one banner.
Authors' comments: 7 pages
WISE J080822.18$-$644357.3 was recently identified as a new M dwarf debris
disc system and a candidate member of the 45 Myr-old Carina association. Given
that the strength of its infrared excess ($L_{\rm IR}/L_{\rm \star}\simeq 0.1$)
appears to be more consistent with a young protoplanetary disc, we present the
first optical spectra of the star and reassess its evolutionary and membership
status. We find WISE J0808-6443 to be a Li-rich M5 star with strong H$\alpha$
emission ($-125 < \textrm{EW} < -65$ \AA\ over 4 epochs) whose strength and
broad width are consistent with accretion at a low level ($\sim$10$^{-10}$
$M_{\odot}$ yr$^{-1}$) from its disc. The spectral energy distribution of the
star is consistent with a primordial disc and is well-fit using a
two-temperature blackbody model with $T_{\rm inner}\simeq$ 1100 K and $T_{\rm
outer}\simeq$ 240 K. AllWISE multi-epoch photometry shows the system exhibits
significant variability in the 3.4 $\mu$m and 4.6 $\mu$m bands. We calculate an
improved proper motion based on archival astrometry, and combined with a new
radial velocity, the kinematics of the star are consistent with membership in
Carina at a kinematic distance of $90\pm9$ pc. The spectroscopic and
photometric data are consistent with WISE J0808-6443 being a $\sim$0.1
$M_{\odot}$ Classical T-Tauri star and one of the oldest known accreting M-type
stars. These results provide further evidence that the upper limit on the
lifetimes of gas-rich discs - and hence the timescales to form and evolve
protoplanetary systems - around the lowest mass stars may be longer than
previously recognised, or some mechanism may be responsible for regenerating
short-lived discs at later stages of pre-main sequence evolution.
Authors' comments: MNRAS Accepted; 12 pages, 10 figures and 5 tables
Permutation testing is a non-parametric method for obtaining the max null
distribution used to compute corrected $p$-values that provide strong control
of false positives. In neuroimaging, however, the computational burden of
running such an algorithm can be significant. We find that by viewing the
permutation testing procedure as the construction of a very large permutation
testing matrix, $T$, one can exploit structural properties derived from the
data and the test statistics to reduce the runtime under certain conditions. In
particular, we see that $T$ is low-rank plus a low-variance residual. This
makes $T$ a good candidate for low-rank matrix completion, where only a very
small number of entries of $T$ ($\sim0.35\%$ of all entries in our experiments)
have to be computed to obtain a good estimate. Based on this observation, we
present RapidPT, an algorithm that efficiently recovers the max null
distribution commonly obtained through regular permutation testing in
voxel-wise analysis. We present an extensive validation on a synthetic dataset
and four varying sized datasets against two baselines: Statistical
NonParametric Mapping (SnPM13) and a standard permutation testing
implementation (referred as NaivePT). We find that RapidPT achieves its best
runtime performance on medium sized datasets ($50 \leq n \leq 200$), with
speedups of 1.5x - 38x (vs. SnPM13) and 20x-1000x (vs. NaivePT). For larger
datasets ($n \geq 200$) RapidPT outperforms NaivePT (6x - 200x) on all
datasets, and provides large speedups over SnPM13 when more than 10000
permutations (2x - 15x) are needed. The implementation is a standalone toolbox
and also integrated within SnPM13, able to leverage multi-core architectures
when available.
Authors' comments: 36 pages, 16 figures
Answering a question asked by Agol and Wise, we show that a desired stronger
form of Wise's malnormal special quotient theorem does not hold. The
counterexamples are generalizations of triangle groups, built using the
Ramanujan graphs constructed by Lubotzky--Phillips--Sarnak.
Authors' comments: 12 pages, 2 figures. V2 incorporates changes suggested by the
referee. This is the final version accepted for publication
