title: A Free Boundary Problem Related to Thermal Insulation: Flat Implies Smooth

abstract: We study the regularity of the interface for a new free boundary problem introduced by Caffarelli and Kriventsov. We show that for minimizers of the functional \[ F_1(A,u) = \int_A |\nabla u|^2 d\mathcal{L}^n + \int_{\partial A} u^2 + \bar{C} \mathcal{L}^n(A) \] over all pairs $(A,u)$ of open sets $A$ containing a fixed set $\Omega$ and functions $u\in H^1(A)$ which equal $1$ on $\Omega$, the boundary $\partial A$ locally coincides with the union of the graphs of two $C^{1,\alpha}$ functions near most points. Specifically, this happens at all points where the interface is trapped between two planes which are sufficiently close together. The proof combines ideas introduced by Ambrosio, Fusco, and Pallara for the Mumford-Shah functional with new arguments specific to the problem considered.


title: Staleness-aware Async-SGD for Distributed Deep Learning

abstract: Deep neural networks have been shown to achieve state-of-the-art performance in several machine learning tasks. Stochastic Gradient Descent (SGD) is the preferred optimization algorithm for training these networks and asynchronous SGD (ASGD) has been widely adopted for accelerating the training of large-scale deep networks in a distributed computing environment. However, in practice it is quite challenging to tune the training hyperparameters (such as learning rate) when using ASGD so as achieve convergence and linear speedup, since the stability of the optimization algorithm is strongly influenced by the asynchronous nature of parameter updates. In this paper, we propose a variant of the ASGD algorithm in which the learning rate is modulated according to the gradient staleness and provide theoretical guarantees for convergence of this algorithm. Experimental verification is performed on commonly-used image classification benchmarks: CIFAR10 and Imagenet to demonstrate the superior effectiveness of the proposed approach, compared to SSGD (Synchronous SGD) and the conventional ASGD algorithm.


title: Small symplectic Calabi-Yau surfaces and exotic 4-manifolds via genus-3 pencils

abstract: We explicitly produce symplectic genus-3 Lefschetz pencils (with base points), whose total spaces are homeomorphic but not diffeomorphic to rational surfaces CP^2 # p (-CP^2) for p= 7, 8, 9. We then give a new construction of an infinite family of symplectic Calabi-Yau surfaces with first Betti number b_1=2,3, along with a surface with b_1=4 homeomorphic to the 4-torus. These are presented as the total spaces of symplectic genus-3 Lefschetz pencils we construct via new positive factorizations in the mapping class group of a genus-3 surface. Our techniques in addition allow us to answer in the negative a question of Korkmaz regarding the upper bound on b_1 of a genus-g fibration.


title: Prioritized Experience Replay

abstract: Experience replay lets online reinforcement learning agents remember and reuse experiences from the past. In prior work, experience transitions were uniformly sampled from a replay memory. However, this approach simply replays transitions at the same frequency that they were originally experienced, regardless of their significance. In this paper we develop a framework for prioritizing experience, so as to replay important transitions more frequently, and therefore learn more efficiently. We use prioritized experience replay in Deep Q-Networks (DQN), a reinforcement learning algorithm that achieved human-level performance across many Atari games. DQN with prioritized experience replay achieves a new state-of-the-art, outperforming DQN with uniform replay on 41 out of 49 games.


title: The Bogoliubov free energy functional II. The dilute limit

abstract: We analyse the canonical Bogoliubov free energy functional at low temperatures in the dilute limit. We prove existence of a first order phase transition and, in the limit $a_0\to a$, we determine the critical temperature to be $T_{\rm{c}}=T_{\rm{fc}}(1+1.49(\rho^{1/3}a))$ to leading order. Here, $T_{\rm{fc}}$ is the critical temperature of the free Bose gas, $\rho$ is the density of the gas, $a$ is the scattering length of the pair-interaction potential $V$, and $a_0=(8\pi)^{-1}\widehat{V}(0)$ its first order approximation. We also prove asymptotic expansions for the free energy. In particular, we recover the Lee-Huang-Yang formula in the limit $a_0\to a$.


title: WASP-135b: a highly irradiated, inflated hot Jupiter orbiting a G5V star

abstract: We report the discovery of a new transiting planet from the WASP survey. WASP-135b is a hot Jupiter with a radius of 1.30 pm 0.09 Rjup, a mass of 1.90 pm 0.08 Mjup and an orbital period of 1.401 days. Its host is a Sun-like star, with a G5 spectral type and a mass and radius of 0.98 pm 0.06 Msun and 0.96 pm 0.05 Rsun respectively. The proximity of the planet to its host means that WASP-135b receives high levels of insolation, which may be the cause of its inflated radius. Additionally, we find weak evidence of a transfer of angular momentum from the planet to its star.


title: An information-theoretic approach to the gravitational-wave burst detection problem

abstract: The observational era of gravitational-wave astronomy began in the Fall of 2015 with the detection of GW150914. One potential type of detectable gravitational wave is short-duration gravitational-wave bursts, whose waveforms can be difficult to predict. We present the framework for a new detection algorithm for such burst events -- \textit{oLIB} -- that can be used in low-latency to identify gravitational-wave transients independently of other search algorithms. This algorithm consists of 1) an excess-power event generator based on the Q-transform -- \textit{Omicron} --, 2) coincidence of these events across a detector network, and 3) an analysis of the coincident events using a Markov chain Monte Carlo Bayesian evidence calculator -- \textit{LALInferenceBurst}. These steps compress the full data streams into a set of Bayes factors for each event; through this process, we use elements from information theory to minimize the amount of information regarding the signal-versus-noise hypothesis that is lost. We optimally extract this information using a likelihood-ratio test to estimate a detection significance for each event. Using representative archival LIGO data, we show that the algorithm can detect gravitational-wave burst events of astrophysical strength in realistic instrumental noise across different burst waveform morphologies. We also demonstrate that the combination of Bayes factors by means of a likelihood-ratio test can improve the detection efficiency of a gravitational-wave burst search. Finally, we show that oLIB's performance is robust against the choice of gravitational-wave populations used to model the likelihood-ratio test likelihoods.


title: Fundamental Limits to Moore's Law

abstract: The theoretical and practical aspects of the fundamental, ultimate, physical limits to scaling, or Moore-s law, is presented.


title: Supplementary Materials for "How to Avoid Reidentification with Proper Anonymization"- Comment on "Unique in the shopping mall: on the reidentifiability of credit card metadata"

abstract: The study by De Montjoye et al. ("Science", 30 January 2015, p. 536) claimed that most individuals can be reidentified from a deidentified credit card transaction database and that anonymization mechanisms are not effective against reidentification. Such claims deserve detailed quantitative scrutiny, as they might seriously undermine the willingness of data owners and subjects to share data for research. In a recent Technical Comment published in "Science" (18 March 2016, p. 1274), we demonstrate that the reidentification risk reported by De Montjoye et al. was significantly overestimated (due to a misunderstanding of the reidentification attack) and that the alleged ineffectiveness of anonymization is due to the choice of poor and undocumented methods and to a general disregard of 40 years of anonymization literature. The technical comment also shows how to properly anonymize data, in order to reduce unequivocal reidentifications to zero while retaining even more analytical utility than with the poor anonymization mechanisms employed by De Montjoye et al. In conclusion, data owners, subjects and users can be reassured that sound privacy models and anonymization methods exist to produce safe and useful anonymized data. Supplementary materials detailing the data sets, algorithms and extended results of our study are available here. Moreover, unlike the De Montjoye et al.'s data set, which was never made available, our data, anonymized results, and anonymization algorithms can be freely downloaded from http://crises-deim.urv.cat/opendata/SPD_Science.zip


title: Comparative Design, Scaling, and Control of Appendages for Inertial Reorientation

abstract: This paper develops a comparative framework for the design of actuated inertial appendages for planar, aerial reorientation. We define the Inertial Reorientation template, the simplest model of this behavior, and leverage its linear dynamics to reveal the design constraints linking a task with the body designs capable of completing it. As practicable inertial appendage designs lead to morphology that is generally more complex, we advance a notion of "anchoring" whereby a judicious choice of physical design in concert with an appropriate control policy yields a system whose closed loop dynamics are sufficiently captured by the template as to permit all further design to take place in its far simpler parameter space. This approach is effective and accurate over the diverse design spaces afforded by existing platforms, enabling performance comparison through the shared task space. We analyze examples from the literature and find advantages to each body type, but conclude that tails provide the highest potential performance for reasonable designs. Thus motivated, we build a physical example by retrofitting a tail to a RHex robot and present empirical evidence of its efficacy.


title: Convergent Perturbation Theory for the lattice $\phi^4$-model

abstract: The standard lattice perturbation theory leads to the asymptotic series because of the incorrect interchange of the summation and integration. However, changing the initial approximation of the perturbation theory, one can generate the convergent series. We study the lattice $\phi^4$-model and compare the operator $\langle\phi_n^2\rangle$ calculated using the convergent series and obtained by Monte Carlo simulations.


title: ABC-CNN: An Attention Based Convolutional Neural Network for Visual Question Answering

abstract: We propose a novel attention based deep learning architecture for visual question answering task (VQA). Given an image and an image related natural language question, VQA generates the natural language answer for the question. Generating the correct answers requires the model's attention to focus on the regions corresponding to the question, because different questions inquire about the attributes of different image regions. We introduce an attention based configurable convolutional neural network (ABC-CNN) to learn such question-guided attention. ABC-CNN determines an attention map for an image-question pair by convolving the image feature map with configurable convolutional kernels derived from the question's semantics. We evaluate the ABC-CNN architecture on three benchmark VQA datasets: Toronto COCO-QA, DAQUAR, and VQA dataset. ABC-CNN model achieves significant improvements over state-of-the-art methods on these datasets. The question-guided attention generated by ABC-CNN is also shown to reflect the regions that are highly relevant to the questions.


title: Mapping out spin and particle conductances in a quantum point contact

abstract: We study particle and spin transport in a single mode quantum point contact using a charge neutral, quantum degenerate Fermi gas with tunable, attractive interactions. This yields the spin and particle conductance of the point contact as a function of chemical potential or confinement. The measurements cover a regime from weak attraction, where quantized conductance is observed, to the resonantly interacting superfluid. Spin conductance exhibits a broad maximum when varying the chemical potential at moderate interactions, which signals the emergence of Cooper pairing. In contrast, the particle conductance is unexpectedly enhanced even before the gas is expected to turn into a superfluid, continuously rising from the plateau at 1/h for weak interactions to plateaux-like features at non-universal values as high as 4/h for intermediate interactions. For strong interactions, the particle conductance plateaux disappear and the spin conductance gets suppressed, confirming the spin-insulating character of a superfluid. Our observations document the breakdown of universal conductance quantization as many-body correlations appear. The observed anomalous quantization challenges a Fermi liquid description of the normal phase, shedding new light on the nature of the strongly attractive Fermi gases.


title: Searching for scalar gravitational interactions in current and future cosmological data

abstract: Modified gravity theories often contain a scalar field of gravitational strength which interacts with matter. We examine constraints on the range and the coupling strength of a scalar gravitational degree of freedom using a subset of current data that can be safely analyzed within the linear perturbation theory. Using a model-independent implementation of scalar-tensor theories in MGCAMB in terms of two functions of the scale factor describing the mass and the coupling of the scalar degree of freedom, we derive constraints on the $f(R)$, generalized chameleon, Symmetron and Dilaton models. Since most of the large scale structure data available today is from relatively low redshifts, only a limited range of observed scales is in the linear regime, leading to relatively weak constraints. We then perform a forecast for a future large scale structure survey, such as Large Synoptic Survey Telescope (LSST), which will map a significant volume at higher redshifts, and show that it will produce much stronger constraints on scalar interactions in specific models. We also perform a principal component analysis and find that future surveys should be able to provide tight constraints on several eigenmodes of the scalar mass evolution.


title: Global equilibrium and local thermodynamics in stationary spacetimes

abstract: In stationary spacetimes global equilibrium states can be defined, applying the maximum entropy principle, by the introduction of local thermodynamic fields determined solely by geometry. As an example, we study a class of equilibrium states for a scalar field in Einstein's static universe, characterized by inhomogeneous thermodynamic properties and nonvanishing angular momentum.


title: The Bino Variations: Effective Field Theory Methods for Dark Matter Direct Detection

abstract: We apply effective field theory methods to compute bino-nucleon scattering, in the case where tree-level interactions are suppressed and the leading contribution is at loop order via heavy flavor squarks or sleptons. We find that leading log corrections to fixed-order calculations can increase the bino mass reach of direct detection experiments by a factor of two in some models. These effects are particularly large for the bino-sbottom coannihilation region, where bino dark matter as heavy as 5-10 TeV may be detected by near future experiments. For the case of stop- and selectron-loop mediated scattering, an experiment reaching the neutrino background will probe thermal binos as heavy as 500 and 300 GeV, respectively. We present three key examples that illustrate in detail the framework for determining weak scale coefficients, and for mapping onto a low energy theory at hadronic scales, through a sequence of effective theories and renormalization group evolution. For the case of a squark degenerate with the bino, we extend the framework to include a squark degree of freedom at low energies using heavy particle effective theory, thus accounting for large logarithms through a "heavy-light current." Benchmark predictions for scattering cross sections are evaluated, including complete leading order matching onto quark and gluon operators, and a systematic treatment of perturbative and hadronic uncertainties.


title: The frequency of accretion disks around single stars: Chamaeleon I

abstract: It is well known that stellar companions can influence the evolution of a protoplanetary disk. Nevertheless, previous disk surveys did not - and could not - consistently exclude binaries from their samples. We present a study dedicated to investigating the frequency of ongoing disk accretion around single stars in a star-forming region. We obtained near-infrared spectroscopy of 54 low-mass stars selected from a high-angular resolution survey in the 2-3 Myr-old Chamaeleon I region to determine the presence of Brackett-$\gamma$ emission, taking the residual chance of undetected multiplicity into account, which we estimate to be on the order of 30%. The result is compared with previous surveys of the same feature in binary stars of the same region to provide a robust estimate of the difference between the accretor fractions of single stars and individual components of binary systems. We find Br$\gamma$ emission among $39.5^{+14.0}_{-9.9}$% of single stars, which is a significantly higher fraction than for binary stars in Chamaeleon I. In particular, close binary systems with separations <100 AU show emission in only $6.5^{+16.5}_{-3.0}$% of the cases according to the same analysis. The emitter frequency of wider binaries appears consistent with the single star value. Interpreting Br$\gamma$ emission as a sign of ongoing accretion and correcting for sensitivity bias, we infer an accretor fraction of single stars of F_acc=$47.8^{+14.0}_{-9.9}$%. This is slightly higher but consistent with previous estimates that do not clearly exclude binaries from their samples. Through our robust and consistent analysis, we confirm that the fraction of young single stars harboring accretion disks is much larger than that of close binaries at the same age. Our findings have important implications for the timescales of disk evolution and planet formation.


title: Penetrative internally heated convection in two and three dimensions

abstract: Convection of an internally heated fluid, confined between top and bottom plates of equal temperature, is studied by direct numerical simulation in two and three dimensions. The unstably stratified upper region drives convection that penetrates into the stably stratified lower region. The fraction of produced heat escaping across the bottom plate, which is one half without convection, initially decreases as convection strengthens. Entering the turbulent regime, this decrease reverses in two dimensions but continues monotonically in three dimensions. The mean fluid temperature, which grows proportionally to the heating rate ($H$) without convection, grows proportionally to $H^{4/5}$ when convection is strong in both two and three dimensions. The ratio of the heating rate to the fluid temperature is likened to the Nusselt number of Rayleigh-B\'enard convection. Simulations are reported for Prandtl numbers between 0.1 and 10 and for Rayleigh numbers (defined in terms of the heating rate) up to $5\times10^{10}$.


title: The CTA aims at the Inert Doublet Model

abstract: We show that the Cherenkov Telescope Array (CTA) can realistically challenge the Inert Doublet Model, excluding its heavy regime up to dark matter masses of 800 GeV and probing a large fraction of the remaining viable parameter space at even higher masses. Two features of the Inert Doublet Model make it particularly suitable for CTA searches. First, the dark matter mass (in the heavy regime) must be larger than 500 GeV. Second, the dark matter annihilation cross section, $\sigma v$, is always larger than the thermal one, reaching values as high as $10^{-25} \mathrm{cm^3s^{-1}}$. This higher value of $\sigma v$ is the result of the unavoidable coannhilation effects that determine the relic density via thermal freeze-out in the early Universe. We find that with 100 hours of Galactic Center exposure, CTA's expected limit widely surpasses, even after the inclusion of systematic errors, current and projected bounds from Fermi-LAT and HESS on this model.


title: Star Formation Suppression Due to Jet Feedback in Radio Galaxies with Shocked Warm Molecular Gas

abstract: We present Herschel observations of 22 radio galaxies, selected for the presence of shocked, warm molecular hydrogen emission. We measured and modeled spectral energy distributions (SEDs) in 33 bands from the ultraviolet to the far-infrared to investigate the impact of jet feedback on star formation activity. These galaxies are massive, early-type galaxies with normal gas-to-dust ratios, covering a range of optical and infrared colors. We find that the star formation rate (SFR) is suppressed by a factor of ~3-6, depending on how molecular gas mass is estimated. We suggest this suppression is due to the shocks driven by the radio jets injecting turbulence into the interstellar medium (ISM), which also powers the luminous warm H2 line emission. Approximately 25% of the sample shows suppression by more than a factor of 10. However, the degree of SFR suppression does not correlate with indicators of jet feedback including jet power, diffuse X-ray emission, or intensity of warm molecular H2 emission, suggesting that while injected turbulence likely impacts star formation, the process is not purely parameterized by the amount of mechanical energy dissipated into the ISM. Radio galaxies with shocked warm molecular gas cover a wide range in SFR-stellar mass space, indicating that these galaxies are in a variety of evolutionary states, from actively star-forming and gas-rich to quiescent and gas-poor. SFR suppression appears to have the largest impact on the evolution of galaxies that are moderately gas-rich.


title: Parameter inference with estimated covariance matrices

abstract: When inferring parameters from a Gaussian-distributed data set by computing a likelihood, a covariance matrix is needed that describes the data errors and their correlations. If the covariance matrix is not known a priori, it may be estimated and thereby becomes a random object with some intrinsic uncertainty itself. We show how to infer parameters in the presence of such an estimated covariance matrix, by marginalising over the true covariance matrix, conditioned on its estimated value. This leads to a likelihood function that is no longer Gaussian, but rather an adapted version of a multivariate t-distribution, which has the same numerical complexity as the multivariate Gaussian. As expected, marginalisation over the true covariance matrix improves inference when compared with Hartlap et al.'s method, which uses an unbiased estimate of the inverse covariance matrix but still assumes that the likelihood is Gaussian.


title: Incoherent thermal transport from dirty black holes

abstract: We study thermal transport in strongly disordered, strongly interacting quantum field theories without quasiparticles using gauge-gravity duality. We analyze linear perturbations of black holes with broken translational symmetry in Einstein-Maxwell-dilaton theories of gravity. Using general geometric arguments in the bulk, we derive bounds on thermal conductivity for the dual disordered field theories in one and two spatial dimensions. In the latter case, the thermal conductivity is always non-zero at finite temperature, so long as the dilaton potential is bounded from below. Hence, generic holographic models make non-trivial predictions about the thermal conductivity in a strongly disordered, strongly coupled metal in two spatial dimensions.


title: Lagrangian Methods Of Cosmic Web Classification

abstract: The cosmic web defines the large scale distribution of matter we see in the Universe today. Classifying the cosmic web into voids, sheets, filaments and nodes allows one to explore structure formation and the role environmental factors have on halo and galaxy properties. While existing studies of cosmic web classification concentrate on grid based methods, this work explores a Lagrangian approach where the V-web algorithm proposed by Hoffman et al. (2012) is implemented with techniques borrowed from smoothed particle hydrodynamics. The Lagrangian approach allows one to classify individual objects (e.g. particles or halos) based on properties of their nearest neighbours in an adaptive manner. It can be applied directly to a halo sample which dramatically reduces computational cost and potentially allows an application of this classification scheme to observed galaxy samples. Finally, the Lagrangian nature admits a straight forward inclusion of the Hubble flow negating the necessity of a visually defined threshold value which is commonly employed by grid based classification methods.


title: Design Considerations of Biaxially Tensile-Strained Germanium-on-Silicon Lasers

abstract: Physical models of Ge energy band structure and material loss were implemented in LASTIP(TM), a 2D simulation tool for edge emitting laser diodes. The model calculation is able to match experimental data available. Important design parameters of a Fabry-Perot Ge laser, such as the cavity length, thickness, width, polycrystalline Si cladding layer thickness were studied and optimized. The laser structure optimizations alone were shown to reduce Ith by 22-fold and increase the differential efficiency by 11 times. The simulations also showed that improving the defect limited carrier lifetime is critical for achieving an efficient and low-threshold Ge laser. With the optimized structure design (300 micron for the cavity length, 0.4 micron for the cavity width, 0.3 micron for the cavity thickness, and 0.6 micron for the polycrystalline Si cladding layer thickness) and a defect limited carrier lifetime of 100 ns, a wall-plug efficiency of 14.6% at 1mW output is predicted, where Jth of 2.8 kA/cm2, Ith of 3.3 mA, I_1mA of 9 mA, and differential efficiency of 23.6% can also be achieved. These are tremendous improvements from the available experimental values at 280 kA/cm2, 756 mA, 837 mA and 1.9%, respectively.


title: Evidence for Gas from a Disintegrating Extrasolar Asteroid

abstract: We report high-resolution spectroscopic observations of WD 1145+017 -- a white dwarf that recently has been found to be transitted by multiple asteroid-sized objects within its tidal radius. We have discovered numerous circumstellar absorption lines with linewidths of $\sim$ 300 km s$^{-1}$ from Mg, Ca, Ti, Cr, Mn, Fe and Ni, possibly from several gas streams produced by collisions among the actively disintegrating objects. The atmosphere of WD 1145+017 is polluted with 11 heavy elements, including O, Mg, Al, Si, Ca, Ti, V:, Cr, Mn, Fe and Ni. Evidently, we are witnessing the active disintegration and subsequent accretion of an extrasolar asteroid.


title: Affleck-Dine baryogenesis just after inflation

abstract: We propose a new scenario of Affleck-Dine baryogenesis where a flat direction in the MSSM generates B-L asymmetry just after the end of inflation. The resulting amount of baryon asymmetry is independent of low-energy supersymmetric models but is dependent on inflation models. We consider the hybrid and chaotic inflation models and find that reheating temperature is required to be higher than that in the conventional scenario of Affleck-Dine baryogenesis. In particular, non-thermal gravitino-overproduction problem is naturally avoided in the hybrid inflation model. Our results imply that Affleck-Dine baryogenesis can be realized in a broader range of supersymmetry and inflation models than expected in the literature.


title: Rainbow Vacua of Colored Higher Spin Gravity in Three Dimensions

abstract: We study the color-decoration of higher-spin (anti)-de Sitter gravity in three dimensions. We show that the rainbow vacua, which we found recently for the colored gravity theory, also pertain in the colored higher-spin theory. The color singlet spin-two plays the role of first fundamental form (metric). The difference is that when spontaneous breaking of color symmetry takes place, the Goldstone modes of massless spin-two combine with all other spins and become the maximal-depth partially massless fields of the highest spin in the theory, forming a Regge trajectory.


title: Stratifying systems over the hereditary path algebra with quiver $\mathbb{A}_{p,q}$

abstract: The authors have proved in [J. Algebra Appl. 14 (2015), no. 6] that the size of a stratifying system over a finite-dimensional hereditary path algebra $A$ is at most $n$, where $n$ is the number of isomorphism classes of simple $A$-modules. Moreover, if $A$ is of Euclidean type a stratifying system over $A$ has at most $n-2$ regular modules. In this work, we construct a family of stratifying systems of size $n$ with a maximal number of regular elements, over the hereditary path algebra with quiver $\widetilde{\mathbb {A}}_{p,q} $, canonically oriented.


title: What is QFT? Resurgent trans-series, Lefschetz thimbles, and new exact saddles

abstract: This is an introductory level review of recent applications of resurgent trans-series and Picard-Lefschetz theory to quantum mechanics and quantum field theory. Resurgence connects local perturbative data with global topological structure. In quantum mechanical systems, this program provides a constructive relation between different saddles. For example, in certain cases it has been shown that all information around the instanton saddle is encoded in perturbation theory around the perturbative saddle. In quantum field theory, such as sigma models compactified on a circle, neutral bions provide a semi-classical interpretation of the elusive IR-renormalon, and fractional kink instantons lead to the non-perturbatively induced gap, of order of the strong scale. In the path integral formulation of quantum mechanics, saddles must be found by solving the holomorphic Newton's equation in the inverted (holomorphized) potential. Some saddles are complex, multi-valued, and even singular, but of finite action, and their inclusion is strictly necessary to prevent inconsistencies. The multi-valued saddles enter either via resurgent cancellations, or their phase is tied with a hidden topological angle. We emphasize the importance of the destructive/constructive interference effects between equally dominant saddles in the Lefschetz thimble decomposition. This is especially important in the context of the sign problem.


title: Fossil groups in the Millennium simulation. From the brightest to the faintest galaxies during the past 8 Gyr

abstract: We investigate the evolution of bright and faint galaxies in fossil and non-fossil groups. We used mock galaxies constructed based on the Millennium run simulation II. We identified fossil groups at redshift zero according to two different selection criteria, and then built reliable control samples of non-fossil groups that reproduce the fossil virial mass and assembly time distributions. The faint galaxies were defined as having r-band absolute magnitudes in the range [-16,-11]. We analysed the properties of the bright and faint galaxies in fossil and non-fossil groups during the past 8 Gyr. We observed that the brightest galaxy in fossil groups is typically brighter and more massive than their counterparts in control groups. Fossil groups developed their large magnitude gap between the brightest galaxies around 3.5 Gyr ago. The brightest galaxy stellar masses of all groups show a notorious increment at that time. By analysing the behaviour of the magnitude gap between the first and the second, third, and fourth ranked galaxies, we found that at earlier times, fossil groups comprised two large brightest galaxies with similar magnitudes surrounded by much fainter galaxies, while in control groups these magnitude gaps were never as large as in fossils. At early times, fossil groups in the faint population were denser than non-fossil groups, then this trend reversed, and finally they became similar at the present day. The mean number of faint galaxies in non-fossil systems increases in an almost constant rate towards later times, while this number in fossil groups reaches a plateau at $z\sim0.6$ that lasts $\sim 2$ Gyr, and then starts growing again more rapidly. The formation of fossil groups is defined at the very beginning of the groups according to their galaxy luminosity sampling, which could be determined by their merging rate at early times.


title: Infinite irredundant equational axiomatisability for a finite monoid

abstract: It is shown that a finite monoid can have an infinite irredundant basis of equations.


title: Evidence for single top-quark production in the $s$-channel in proton-proton collisions at $\sqrt{s}=$8 TeV with the ATLAS detector using the Matrix Element Method

abstract: This Letter presents evidence for single top-quark production in the $s$-channel using proton-proton collisions at a centre-of-mass energy of 8 TeV with the ATLAS detector at the CERN Large Hadron Collider. The analysis is performed on events containing one isolated electron or muon, large missing transverse momentum and exactly two $b$-tagged jets in the final state. The analysed data set corresponds to an integrated luminosity of 20.3 fb$^{-1}$. The signal is extracted using a maximum-likelihood fit of a discriminant which is based on the matrix element method and optimized in order to separate single-top-quark $s$-channel events from the main background contributions, which are top-quark pair production and $W$ boson production in association with heavy-flavour jets. The measurement leads to an observed signal significance of 3.2 standard deviations and a measured cross-section of $\sigma_s=4.8\pm0.8$(stat.)$^{+1.6}_{-1.3}$(syst.) pb, which is consistent with the Standard Model expectation. The expected significance for the analysis is 3.9 standard deviations.


title: Computation of the Madelung constant for hypercubic crystal structures in any dimension

abstract: A new method of computing the Madelung constants for hypercubic crystal structures in any dimension $n\geq 2$ is given. It is shown for $n\geq 3$ that the Madelung constant may be obtained in a simple, efficient and unambiguous way as the Hadamard finite part of the integral representation of the potential within the crystal which is divergent at any point charge location. Such a regularization method fails in the bidimensional case due to the logarithmic nature of singularities for the potential. In that case, a specific approach is proposed taking in account the scale invariance of the Poisson equation and the existence of a finite horizon for each point charge in the plane. Since a closed-form exact solution for the 2D electrostatic potential may be derived, one shows that the Madelung constant may be defined via an appropriate limit calculation as the mean value of potential energies of charges composing the unit cell.


title: Mobility Measurements Probe Conformational Changes in Membrane Proteins due to Tension

abstract: The function of membrane-embedded proteins such as ion channels depends crucially on their conformation. We demonstrate how conformational changes in asymmetric membrane proteins may be inferred from measurements of their diffusion. Such proteins cause local deformations in the membrane, which induce an extra hydrodynamic drag on the protein. Using membrane tension to control the magnitude of the deformations and hence the drag, measurements of diffusivity can be used to infer--- via an elastic model of the protein--- how conformation is changed by tension. Motivated by recent experimental results [Quemeneur et al., Proc. Natl. Acad. Sci. USA, 111 5083 (2014)] we focus on KvAP, a voltage-gated potassium channel. The conformation of KvAP is found to change considerably due to tension, with its `walls', where the protein meets the membrane, undergoing significant angular strains. The torsional stiffness is determined to be 26.8 kT at room temperature. This has implications for both the structure and function of such proteins in the environment of a tension-bearing membrane.


title: Neutrino mass limits: robust information from the power spectrum of galaxy surveys

abstract: We present cosmological upper limits on the sum of active neutrino masses using large-scale power spectrum data from the WiggleZ Dark Energy Survey and from the Sloan Digital Sky Survey - Data Release 7 (SDSS-DR7) sample of Luminous Red Galaxies (LRG). Combining measurements on the Cosmic Microwave Background temperature and polarisation anisotropies by the Planck satellite together with WiggleZ power spectrum results in a neutrino mass bound of 0.37 eV at 95% C.L., while replacing WiggleZ by the SDSS-DR7 LRG power spectrum, the 95% C.L. bound on the sum of neutrino masses is 0.38 eV. Adding Baryon Acoustic Oscillation (BAO) distance scale measurements, the neutrino mass upper limits greatly improve, since BAO data break degeneracies in parameter space. Within a $\Lambda$CDM model, we find an upper limit of 0.13 eV (0.14 eV) at 95% C.L., when using SDSS-DR7 LRG (WiggleZ) together with BAO and Planck. The addition of BAO data makes the neutrino mass upper limit robust, showing only a weak dependence on the power spectrum used. We also quantify the dependence of neutrino mass limit reported here on the CMB lensing information. The tighter upper limit (0.13 eV) obtained with SDSS-DR7 LRG is very close to that recently obtained using Lyman-alpha clustering data, yet uses a completely different probe and redshift range, further supporting the robustness of the constraint. This constraint puts under some pressure the inverted mass hierarchy and favours the normal hierarchy.


title: Universal collisionless transport of graphene

abstract: The impact of the electron-electron Coulomb interaction on the optical conductivity of graphene has led to a controversy that calls into question the universality of collisionless transport in this and other Dirac materials. Using a lattice calculation that avoids divergences present in previous nodal Dirac approaches, our work settles this controversy and obtains results in quantitative agreement with experiment over a wide frequency range. We also demonstrate that dimensional regularization methods agree, as long as the scaling properties of the conductivity and the regularization of the theory in modified dimension are correctly implemented. Tight-binding lattice and nodal Dirac theory calculations are shown to coincide at low energies even when the non-zero size of the atomic orbital wave function is included, conclusively demonstrating the universality of the optical conductivity of graphene.


title: Iron-Group Abundances in the Metal-Poor Main Sequence Turnoff Star HD~84937

abstract: We have derived new very accurate abundances of the Fe-group elements Sc through Zn (Z = 21-30) in the bright main-sequence turnoff star HD 84937, based on high-resolution spectra covering the visible and ultraviolet spectral regions. New or recent laboratory transition data for 14 species of seven elements have been used. Abundances from more than 600 lines of non-Fe species have been combined with about 550 Fe lines in HD 84937 to yield abundance ratios of high precision. The abundances have been determined from both neutral and ionized transitions, which generally are in agreement with each other. We find no substantial departures from standard LTE Saha ionization balance in this [Fe/H] = -2.32 star. Noteworthy among the abundances are: [Co/Fe] = 0.14 and [Cu/Fe] = -0.83, in agreement with past studies abundance trends in this and other low metallicity stars; and <[Sc,Ti,V/Fe]> = 0.31, which has not been noted previously. A detailed examination of scandium, titanium, and vanadium abundances in large-sample spectroscopic surveys reveals that they are positively correlated in stars with [Fe/H] < -2; HD 84937 lies at the high end of this correlation. These trends constrain the synthesis mechanisms of Fe-group elements. We also examine the GCE abundance trends of the Fe-group elements, including a new nucleosynthesis model with jet-like explosion effects.


title: Computing with Harmonic Functions

abstract: This document is the manual for a free Mathematica package for computing with harmonic functions. This package allows the user to make calculations that would take a prohibitive amount of time if done without a computer. For example, the Poisson integral of any polynomial can be computed exactly. This software can find exact solutions to Dirichlet, Neumann, and biDirichlet problems in R^n with polynomial data on balls, ellipsoids, and annular regions. It can also find bases for spaces of spherical harmonics, compute projections onto the harmonic Bergman space, and perform other manipulations with harmonic functions.


title: Algorithms for Communication Problems for Mobile Agents Exchanging Energy

abstract: We consider communication problems in the setting of mobile agents deployed in an edge-weighted network. The assumption of the paper is that each agent has some energy that it can transfer to any other agent when they meet (together with the information it holds). The paper deals with three communication problems: data delivery,convergecast and broadcast. These problems are posed for a centralized scheduler which has full knowledge of the instance. It is already known that, without energy exchange, all three problems are NP-complete even if the network is a line. Surprisingly, if we allow the agents to exchange energy, we show that all three problems are polynomially solvable on trees and have linear time algorithms on the line. On the other hand for general undirected and directed graphs we show that these problems, even if energy exchange is allowed, are still NP-complete.


title: Spectral Model of Non-Stationary, Inhomogeneous Turbulence

abstract: We compare results from a spectral model for non-stationary, inhomogeneous turbulence (Besnard et al., Theor. Comp. Fluid. Dyn., vol. 8, pp 1-35, 1996) with Direct Numerical Simulation (DNS) data of a shear-free mixing layer (SFML) (Tordella et al., Phys. Rev. E, vol. 77, 016309, 2008). The SFML is used as a test case in which the efficacy of the model closure for the physical-space transport of the fluid velocity field can be tested in a flow with inhomogeneity, without the additional complexity of mean-flow coupling. The model is able to capture certain features of the SFML quite well for intermediate to long-times, including the evolution of the mixing-layer width and turbulent kinetic energy. At short-times, and for more sensitive statistics such as the generation of the velocity field anisotropy, the model is less accurate. We present arguments, supported by the DNS data, that a significant cause of the discrepancies is the local approximation to the intrinsically non-local pressure-transport in physical-space that was made in the model, the effects of which would be particularly strong at short-times when the inhomogeneity of the SFML is strongest.


title: A Model of Quark and Lepton Mixing and Mass Hierarchy

abstract: It is shown that an idea proposed in 1996 that relates in a qualitatively correct way the inter-family mass hierarchies of the up quarks, down quarks, charged leptons, and neutrinos, can be combined with a predictive scheme recently proposed for relating quark mixing and neutrino mixing. In the resulting model, the entire flavor structure of the quarks and leptons is expressible in terms of two "master matrices": a diagonal matrix that gives the inter-family mass ratios, and an off-diagonal matrix that controls all flavor mixing.


title: Mu_x boosted-bottom-jet tagging and Z-prime boson searches

abstract: We present a new technique for tagging heavy-flavor jets with p_T > 500 GeV called "mu_x tagging." Current track-based methods of b-jet tagging lose efficiency and experience a large rise in fake rate in the boosted regime. Using muons from B hadron decay, we combine angular information and jet substructure to tag b jets, c jets, light jets, and "light-heavy" jets (those containing B hadrons from gluon splitting). We find tagging efficiencies of epsilon_b = 14%, epsilon_c = 6.5%, epsilon_{light-light} = 0.14%, and epsilon_{light-heavy} = 0.5%, respectively, that are nearly independent of transverse momentum at high energy. We demonstrate the usefulness of this new scheme by examining the discovery potential for multi-TeV leptophobic Z-prime bosons in the boosted-b-tagged dijet channel at the Large Hadron Collider.


title: Supergravity Solutions with $AdS_4$ from Non-Abelian T-Dualities

abstract: We present a large class of new backgrounds that are solutions of type II supergravity with a warped AdS${}_4$ factor, non-trivial axion-dilaton, B-field, and three- and five-form Ramond-Ramond fluxes. We obtain these solutions by applying non-Abelian T-dualities with respect to SU(2) or SU(2)/U(1) isometries to reductions to 10d IIA of 11d sugra solutions of the form AdS${}_4 \times Y^7$, with $Y^7 = S^7/\mathbb{Z}_k, S^7, M^{1,1,1}, Q^{1,1,1}$ and $N(1,1)$. The main class of reductions to IIA is along the Hopf fiber and leads to solutions of the form $AdS_4 \times K_6$, where $K_6 $ is K\"ahler Einstein with $K_6=\mathbb{CP}^3, S^2\times \mathbb{CP}^2, S^2\times S^2 \times S^2$; the first member of this class is dual to the ABJM field theory in the 't Hooft limit. We also consider other less symmetric but susy preserving reductions along circles that are not the Hopf fiber. In the case of $N(1,1)$ we find an additional breaking of isometries in the NAT-dual background. To initiate the study of some properties of the field theory dual, we explicitly compute the central charge holographically.


title: Models of the Thermal Evolution of the Intergalactic Medium After Reionization

abstract: Recent years have brought more precise temperature measurements of the low-density intergalactic medium (IGM). These new measurements constrain the processes that heated the IGM, such as the reionization of H I and of He II. We present a semi-analytical model for the thermal history of the IGM that follows the photoheating history of primordial gas. Our model adopts a multizone approach that, compared to previous models, more accurately captures the inhomogeneous heating and cooling of the IGM during patchy reionization. We compare our model with recent temperature measurements spanning z = 1.6-4.8, finding that these measurements are consistent with scenarios in which the He II was reionized at z = 3-4 by quasars. Significantly longer duration or higher redshift He II reionization scenarios are ruled out by the measurements. For hydrogen reionization, we find that only low-redshift and high-temperature scenarios are excluded. For example, a model in which the IGM was heated to 30 000 K when an ionization front passed, and with hydrogen reionization occurring over 6 < z < 9, is ruled out. Finally, we place constraints on how much heating could owe to TeV blazars, cosmic rays, and other non-standard mechanisms. We find that by z = 2, a maximum of 1 eV of additional heat could be injected per baryon over standard photoheating-only models, with this limit becoming <~0.5 eV at z > 3.


title: An algebraic model for rational G-spectra over an exceptional subgroup

abstract: We give a simple algebraic model for rational G-spectra over an exceptional subgroup, for any compact Lie group G. Moreover, all our Quillen equivalences are symmetric monoidal, so as a corollary we obtain a monoidal algebraic model for rational G-spectra when G is finite. We also present a study of the relationship between induction - restriction - coinduction adjunctions and left Bousfield localizations at idempotents of the rational Burnside ring.


title: Gravitational-wave cosmology across 29 decades in frequency

abstract: Quantum fluctuations of the gravitational field in the early Universe, amplified by inflation, produce a primordial gravitational-wave background across a broad frequency band. We derive constraints on the spectrum of this gravitational radiation, and hence on theories of the early Universe, by combining experiments that cover 29 orders of magnitude in frequency. These include Planck observations of cosmic microwave background temperature and polarization power spectra and lensing, together with baryon acoustic oscillations and big bang nucleosynthesis measurements, as well as new pulsar timing array and ground-based interferometer limits. While individual experiments constrain the gravitational-wave energy density in specific frequency bands, the combination of experiments allows us to constrain cosmological parameters, including the inflationary spectral index, $n_t$, and the tensor-to-scalar ratio, $r$. Results from individual experiments include the most stringent nanohertz limit of the primordial background to date from the Parkes Pulsar Timing Array, $\Omega_{\rm gw}(f)<2.3\times10^{-10}$. Observations of the cosmic microwave background alone limit the gravitational-wave spectral index at 95\% confidence to $n_t\lesssim5$ for a tensor-to-scalar ratio of $r = 0.11$. However, the combination of all the above experiments limits $n_t<0.36$. Future Advanced LIGO observations are expected to further constrain $n_t<0.34$ by 2020. When cosmic microwave background experiments detect a non-zero $r$, our results will imply even more stringent constraints on $n_t$ and hence theories of the early Universe.


title: On melting of boron phosphide under pressure

abstract: Melting of cubic boron phosphide, BP has been studied at pressures to 9 GPa using synchrotron X-ray diffraction and electrical resistivity measurements. It has been found that above 2.6 GPa BP melts congruently, and the melting curve exhibits negative slope -60(7) K/GPa, which is indicative of a higher density of the melt as compared to the solid phase.


title: Uncertainty-based Arbitration of Human-Machine Shared Control

abstract: Manufacturing requires consistent production rate and task success for sustainable operation. Some manufacturing tasks require a semi-autonomous approach, exploiting the combination of human adaptability and machine precision and speed, to be cost effective. The main contribution of this paper is a new approach to determine the level of autonomy for human-machine shared control based on the automation uncertainty. Moreover, the haptic feedback is scaled by the level of autonomy to indicate machine confidence to the operator. Experimentation results, with a human-robot peg-in-a-hole testbed, show more than 5 times improvement in the error tolerance for task completion with the shared control approach when compared to a purely autonomous method.


title: Anomalous Microwave Emission from Spinning Dust and its Polarization Spectrum

abstract: Nearly twenty years after the discovery of anomalous microwave emission (AME) that contaminates to the cosmic microwave background (CMB) radiation, its origin remains inconclusive. Observational results from numerous experiments have revealed that AME is most consistent with spinning dust emission from rapidly spinning ultrasmall interstellar grains. In this paper, I will first review our improved model of spinning dust, which treats realistic dynamics of wobbling non-spherical grains, impulsive interactions of grains with ions in the ambient plasma, and some other important effects. I will then discuss recent progress in quantifying the polarization of spinning dust emission from polycyclic aromatic hydrocarbons. I will finish with a brief discussion on remaining issues about the origins of AME.


title: Theory of the Structural Glass Transition: A Pedagogical Review

abstract: The random first-order transition (RFOT) theory of the structural glass transition is reviewed in a pedagogical fashion. The rigidity that emerges in crystals and glassy liquids is of the same fundamental origin. In both cases, it corresponds with a breaking of the translational symmetry; analogies with freezing transitions in spin systems can also be made. The common aspect of these seemingly distinct phenomena is a spontaneous emergence of the molecular field, a venerable and well-understood concept. In crucial distinction from periodic crystallisation, the free energy landscape of a glassy liquid is vastly degenerate, which gives rise to new length and time scales while rendering the emergence of rigidity gradual. We obviate the standard notion that to be mechanically stable a structure must be essentially unique; instead, we show that bulk degeneracy is perfectly allowed but should not exceed a certain value. The present microscopic description thus explains both crystallisation and the emergence of the landscape regime followed by vitrification in a unified, thermodynamics-rooted fashion. The article contains a self-contained exposition of the basics of the classical density functional theory and liquid theory, which are subsequently used to quantitatively estimate, without using adjustable parameters, the key attributes of glassy liquids, viz., the relaxation barriers, glass transition temperature, and cooperativity size. These results are then used to quantitatively discuss many diverse glassy phenomena, including: the intrinsic connection between the excess liquid entropy and relaxation rates, the non-Arrhenius temperature dependence of $\alpha$-relaxation, the dynamic heterogeneity, ... (see Comments for the remainder of Abstract.)


title: Method for detection and reconstruction of gravitational wave transients with networks of advanced detectors

abstract: We present a method for detection and reconstruction of the gravitational-wave (GW) transients with the networks of advanced detectors. Originally designed to search for transients with the initial GW detectors, it uses significantly improved algorithms, which enhances both the low-latency searches with rapid localization of GW events for the electro-magnetic followup and high confidence detection of a broad range of the transient GW sources. In the paper we present the analytic framework of the method. Following a short description of the core analysis algorithms, we introduce a novel approach to the reconstruction of the GW polarization from a pattern of detector responses to a GW signal. This polarization pattern is a unique signature of an arbitrary GW signal that can be measured independent from the other source parameters. The polarization measurements enable rapid reconstruction of the GW waveforms, sky localization and helps identification of the source origin.


title: A Duality Based 2-Approximation Algorithm for Maximum Agreement Forest

abstract: We give a 2-approximation algorithm for the Maximum Agreement Forest problem on two rooted binary trees. This NP-hard problem has been studied extensively in the past two decades, since it can be used to compute the Subtree Prune-and-Regraft (SPR) distance between two phylogenetic trees. Our result improves on the very recent 2.5-approximation algorithm due to Shi, Feng, You and Wang (2015). Our algorithm is the first approximation algorithm for this problem that uses LP duality in its analysis.


title: A pilot study on the daily control capability of s-EMG prosthetic hands by amputees

abstract: Surface electromyography is a valid tool to gather muscular contraction signals from intact and amputated subjects. Electromyographic signals can be used to control prosthetic devices in a noninvasive way distinguishing the movements performed by the particular EMG electrodes activity. According to the literature, several algorithms have been used to control prosthetic hands through s-EMG signals. The main issue is to correctly classify the signals acquired as the movement actually performed. This work presents a study on the Support Vector Machine's performance in a short-time period, gained using two different feature representation (Mean Absolute Value and Waveform Length) of the sEMG signals. In particular, we paid close attention to the repeatability problem, that is the capability to achieve a stable and satisfactory level of accuracy in repeated experiments. Results on a limited setting are encouraging, as they show an average accuracy above 73% even in the worst case scenario.


title: Neutral Higgs Boson Production at e+e- Colliders in the Complex MSSM: A Full One-Loop Analysis

abstract: For the search for additional Higgs bosons in the Minimal Supersymmetric Standard Model (MSSM) as well as for future precision analyses in the Higgs sector a precise knowledge of their production properties is mandatory. We evaluate the cross sections for the neutral Higgs boson production at e^+e^- colliders in the MSSM with complex parameters (cMSSM). The evaluation is based on a full one-loop calculation of the production channels e^+e^- -> h_i Z, h_i gamma, h_i h_j (i,j = 1,2,3), including soft and hard QED radiation. The dependence of the Higgs boson production cross sections on the relevant cMSSM parameters is analyzed numerically. We find sizable contributions to many cross sections. They are, depending on the production channel, roughly of 10-20% of the tree-level results, but can go up to 50% or higher. The full one-loop contributions are important for a future linear e^+e^- collider such as the ILC or CLIC. There are plans to implement the evaluation of the Higgs boson production cross sections into the code FeynHiggs.


title: Measurement of the transverse single-spin asymmetry in $p^\uparrow+p \to W^{\pm}/Z^0$ at RHIC

abstract: We present the measurement of the transverse single-spin asymmetry of weak boson production in transversely polarized proton-proton collisions at $\sqrt{s} = 500~\text{GeV}$ by the STAR experiment at RHIC. The measured observable is sensitive to the Sivers function, one of the transverse momentum dependent parton distribution functions, which is predicted to have the opposite sign in proton-proton collisions from that observed in deep inelastic lepton-proton scattering. These data provide the first experimental investigation of the non-universality of the Sivers function, fundamental to our understanding of QCD.


title: Studying the control of non invasive prosthetic hands over large time spans

abstract: The electromyography (EMG) signal is the electrical manifestation of a neuromuscular activation that provides access to physiological processes which cause the muscle to generate force and produce movement. Non invasive prostheses use such signals detected by the electrodes placed on the user's stump, as input to generate hand posture movements according to the intentions of the prosthesis wearer. The aim of this pilot study is to explore the repeatability issue, i.e. the ability to classify 17 different hand postures, represented by EMG signal, across a time span of days by a control algorithm. Data collection experiments lasted four days and signals were collected from the forearm of a single subject. We find that Support Vector Machine (SVM) classification results are high enough to guarantee a correct classification of more than 10 postures in each moment of the considered time span.


title: Distinguished graduates in mathematics of Jagiellonian University in the interwar period. Part II: 1928-1939

abstract: In this study, we continue to present profiles of some distinguished graduates in mathematics of the Jagiellonian University from the years 1918-1939. We discuss their professional paths and scholarly achievements, instances of scientific collaboration, connections with other academic centers in Poland and worldwide, involvement in mathematical education and teacher training, as well as their later roles in Polish scientific and academic life. We also try to understand in what way they were shaped by their studies and how much of Krak\'ow scientific traditions they continued. We find strong support for the claim that there was a distinct, diverse and deep mathematical stream in Krak\'ow between the wars, rooted in classical disciplines like differential equations and geometry, but also open to new trends in mathematics. Part II concerns people who graduated after the university reform, getting the newly introduced master's degree in mathematics, in 1928-39.


title: HI and CO Velocity Dispersions in Nearby Galaxies

abstract: We analyze the velocity dispersions of individual HI and CO profiles in a number of nearby galaxies from the high-resolution HERACLES CO and THINGS HI surveys. Focusing on regions with bright CO emission, we find a CO dispersion value: 7.3 $\pm$ 1.7 km/s. The corresponding HI dispersion is 11.7 $\pm$ 2.3 km/s, yielding a mean HI/CO dispersion ratio of 1.4 $\pm$ 0.2, independent of radius. We find that the CO velocity dispersion increases towards lower peak fluxes. This is consistent with previous work where we showed that when using spectra averaged ("stacked") over large areas, larger values for the CO dispersion are found, and a lower dispersion ratio: 1.0 $\pm$ 0.2. The stacking method is more sensitive to low-level diffuse emission, whereas individual profiles trace narrow-line, GMC-dominated, bright emission. These results provide further evidence that disk galaxies contain not only a thin, low velocity dispersion, high density CO disk that is dominated by GMCs, but also a fainter, higher dispersion, diffuse disk component.


title: Mutual distance dependence drives the observed jet power - radio luminosity scaling relations in radio galaxies

abstract: The kinetic power of radio jets is a quantity of fundamental importance to studies of the AGN feedback process and radio galaxy physics. A widely used proxy for jet power is the extended radio luminosity. A number of empirical methods have been used to calibrate a scaling relationship between jet power (Q) and radio luminosity (L) of the form log(Q) = beta_L * log(L) + C. The regression slope has typically been found to be beta_L ~ 0.7 - 0.8. Here we show that the previously reported scaling relations are strongly affected by the confounding variable, distance. We find that in a sample of FRI X-ray cavity systems, after accounting for the mutual distance dependence, the jet power and radio luminosity are only weakly correlated, with slope beta_L ~ 0.3: significantly flatter than previously reported. We also find that in previously used samples of high-power sources, no evidence for an intrinsic correlation is present when the effect of distance is accounted for. Using a simple model we show that beta_L is expected to be significantly lower in samples of FRI radio galaxies than it is for FRIIs, due to the differing dynamics for these two classes of radio source. For FRI X-ray cavity systems the model predicts beta_L (FRI) ~ 0.5 in contrast to FRII radio galaxies, for which beta_L(FRII) ~ 0.8. We discuss the implications of our finding for studies of radio mode feedback, and radio galaxy physics.


title: Glacial Cycles and Milankovitch Forcing

abstract: Using a recent conceptual model of the glacial-interglacial cycles we present more evidence of Milankovitch cycles being the trigger for retreat and forming of ice sheets in the cycles. This model is based on a finite approximation of an infinite dimensional model which has three components: Budyko's energy balance model describing the annual mean temperatures at latitudes, Widiasih's ODE which describes the behavior of the edge of the ice sheet, and Walsh et al. who introduced a snow line to account for glacial accumulation and ablation zones. Certain variables in the model are made to depend on the Milankovitch cycles, in particular, the obliquity of the Earth's axis and the eccentricity of the Earth's orbit. We see as a result that deglaciation and glaciation do occur mostly due to obliquity and to some extent eccentricity.


title: Robust self-trapping of vortex beams in a saturable optical medium

abstract: We report the first observation of robust self-trapping of vortex beams propagating in a uniform condensed medium featuring local saturable self-focusing nonlinearity. Optical vortices with topological charge m=1, that remain self-trapped over ~ 5 Rayleigh lengths, are excited in carbon disulfide using a helical light beam at 532 nm and intensities from 8 to 10 GW/cm^2. At larger intensities, the vortex beams lose their stability, spontaneously breaking into two fragments. Numerical simulations based on the nonlinear Schr\"odinger equation including the three-photon absorption and nonpolynomial saturation of the refractive nonlinearity demonstrate close agreement with the experimental findings.


title: A Roth type theorem for dense subsets of $\mathbb{R}^d$

abstract: Let $1 < p < \infty$, $p\neq 2$. We prove that if $d\geq d_p$ is sufficiently large, and $A\subs\R^d$ is a measurable set of positive upper density then there exists $\la_0=\la_0(A)$ such for all $\la\geq\la_0$ there are $x,y\in\R^d$ such that $\{x,x+y,x+2y\}\subs A$ and $|y|_p=\la$, where $||y||_p=(\sum_i |y_i|^p)^{1/p}$ is the $l^p(\mathbb R^d)$-norm of a point $y=(y_1,\ldots,y_d)\in\R^d$. This means that dense subsets of $\R^d$ contain 3-term progressions of all sufficiently large gaps when the gap size is measured in the $l^p$-metric. This statement is known to be false in the Euclidean $l^2$-metric as well as in the $l^1$ and $\ell^{\infty}$-metrics. One of the goals of this note is to understand this phenomenon. A distinctive feature of the proof is the use of multilinear singular integral operators, widely studied in classical time-frequency analysis, in the estimation of forms counting configurations.


title: A passive THz video camera based on lumped element kinetic inductance detectors

abstract: We have developed a passive 350 GHz (850 {\mu}m) video-camera to demonstrate lumped element kinetic inductance detectors (LEKIDs) -- designed originally for far-infrared astronomy -- as an option for general purpose terrestrial terahertz imaging applications. The camera currently operates at a quasi-video frame rate of 2 Hz with a noise equivalent temperature difference per frame of $\sim$0.1 K, which is close to the background limit. The 152 element superconducting LEKID array is fabricated from a simple 40 nm aluminum film on a silicon dielectric substrate and is read out through a single microwave feedline with a cryogenic low noise amplifier and room temperature frequency domain multiplexing electronics.


title: Discovery of polarized line emission in SN1006

abstract: Laming (1990) predicted that the narrow Balmer line core of the ~3000 km/s shock in the SN 1006 remnant would be significantly polarized due to electron and proton impact polarization. Here, based on deep spectrally resolved polarimetry obtained with the European Southern Observatory (ESO)'s Very Large Telescope (VLT), we report the discovery of polarized line emission of polarization degree approx 1.3 percent with position angle orthogonal to the SNR filament. Correcting for an unpolarized broad line component, the implied narrow line polarization is approx 2.0 percent, close to the predictions of Laming (1990). The predicted polarization is primarily sensitive to shock velocity and post-shock temperature equilibration. By measuring polarization for the SN1006 remnant, we validate and enable a new diagnostic that has important applications in a wide variety of astrophysical situations, such as shocks, intense radiation fields, high energy particle streams and conductive interfaces.


title: Statistical Engineering: An Idea Whose Time Has Come?

abstract: Several authors, including the American Statistician (ASA), have noted the challenges facing statisticians when attacking large, complex, unstructured problems, as opposed to well-defined textbook problems. Clearly, the standard paradigm of selecting the one "correct" statistical method for such problems is not sufficient; a new paradigm is needed. Statistical engineering has been proposed as a discipline that can provide a viable paradigm to attack such problems, used in conjunction with sound statistical science. Of course, in order to develop as a true discipline, statistical engineering needs a well-developed theory, not just a formal definition and successful case studies. This article documents and disseminates the current state of the underlying theory of statistical engineering. Our purpose is to provide a vehicle for applied statisticians to further enhance the practice of statistics, and for academics so interested to continue development of the underlying theory of statistical engineering.


title: Regret Analysis of the Finite-Horizon Gittins Index Strategy for Multi-Armed Bandits

abstract: I analyse the frequentist regret of the famous Gittins index strategy for multi-armed bandits with Gaussian noise and a finite horizon. Remarkably it turns out that this approach leads to finite-time regret guarantees comparable to those available for the popular UCB algorithm. Along the way I derive finite-time bounds on the Gittins index that are asymptotically exact and may be of independent interest. I also discuss some computational issues and present experimental results suggesting that a particular version of the Gittins index strategy is a modest improvement on existing algorithms with finite-time regret guarantees such as UCB and Thompson sampling.


title: Active Object Localization with Deep Reinforcement Learning

abstract: We present an active detection model for localizing objects in scenes. The model is class-specific and allows an agent to focus attention on candidate regions for identifying the correct location of a target object. This agent learns to deform a bounding box using simple transformation actions, with the goal of determining the most specific location of target objects following top-down reasoning. The proposed localization agent is trained using deep reinforcement learning, and evaluated on the Pascal VOC 2007 dataset. We show that agents guided by the proposed model are able to localize a single instance of an object after analyzing only between 11 and 25 regions in an image, and obtain the best detection results among systems that do not use object proposals for object localization.


title: Necessary and Sufficient Condition on the Lindblad Equation to Prevent Entropy Increase

abstract: It is shown that in order for the solutions of the Lindblad equation never to give a decreasing von Neumann entropy, it is necessary and sufficient that the operators appearing in this equation should be unitary linear combinations of their adjoints. In this case, these operators may be replaced with Hermitian operators, without changing the evolution of density matrices.


title: Rate of Price Discovery in Iterative Combinatorial Auctions

abstract: We study a class of iterative combinatorial auctions which can be viewed as subgradient descent methods for the problem of pricing bundles to balance supply and demand. We provide concrete convergence rates for auctions in this class, bounding the number of auction rounds needed to reach clearing prices. Our analysis allows for a variety of pricing schemes, including item, bundle, and polynomial pricing, and the respective convergence rates confirm that more expressive pricing schemes come at the cost of slower convergence. We consider two models of bidder behavior. In the first model, bidders behave stochastically according to a random utility model, which includes standard best-response bidding as a special case. In the second model, bidders behave arbitrarily (even adversarially), and meaningful convergence relies on properly designed activity rules.


title: Segmental Recurrent Neural Networks

abstract: We introduce segmental recurrent neural networks (SRNNs) which define, given an input sequence, a joint probability distribution over segmentations of the input and labelings of the segments. Representations of the input segments (i.e., contiguous subsequences of the input) are computed by encoding their constituent tokens using bidirectional recurrent neural nets, and these "segment embeddings" are used to define compatibility scores with output labels. These local compatibility scores are integrated using a global semi-Markov conditional random field. Both fully supervised training -- in which segment boundaries and labels are observed -- as well as partially supervised training -- in which segment boundaries are latent -- are straightforward. Experiments on handwriting recognition and joint Chinese word segmentation/POS tagging show that, compared to models that do not explicitly represent segments such as BIO tagging schemes and connectionist temporal classification (CTC), SRNNs obtain substantially higher accuracies.


title: Decoupling of the Antiferromagnetic and Insulating States in Tb doped Sr2IrO4

abstract: Sr2IrO4 is a spin-orbit coupled insulator with an antiferromagnetic (AFM) transition at TN=240 K. We report results of a comprehensive study of single-crystal Sr2Ir1-xTbxO4. This study found that mere 3% (x=0.03) tetravalent Tb4+(4f7) substituting for Ir4+ (rather than Sr2+) completely suppresses the long-range collinear AFM transition but retains the insulating state, leading to a phase diagram featuring a decoupling of magnetic interactions and charge gap. The insulating state at x=0.03 is characterized by an unusually large specific heat at low temperatures and an incommensurate magnetic state having magnetic peaks at (0.95, 0, 0) and (0, 0.95, 0) in the neutron diffraction, suggesting a spiral or spin density wave order. It is apparent that Tb doping effectively changes the relative strength of the SOI and the tetragonal CEF and enhances the Hund's rule coupling that competes with the SOI, and destabilizes the AFM state. However, the disappearance of the AFM accompanies no metallic state chiefly because an energy level mismatch for the Ir and Tb sites weakens charge carrier hopping and renders a persistent insulating state. This work highlights an unconventional correlation between the AFM and insulating states in which the magnetic transition plays no critical role in the formation of the charge gap in the iridate.


title: A First Szeg\H{o}'s Limit Theorem for a class of non-Toeplitz matrices

abstract: We compute the limiting statistical distribution of the eigenvalues of sequences of matrices whose entries satisfy what we call a vanishing mean variation condition and are $\mu$-distributed for some probability measure. As an application of our results, we extend the well-known class of Kac-Murdock-Szeg\H{o} generalized Toeplitz matrices to sequences of matrices whose diagonal entries are modeled by Riemann integrable functions.


title: Hubble Tarantula Treasury Project. III. Photometric Catalog and Resulting Constraints on the Progression of Star Formation in the 30 Doradus Region

abstract: We present and describe the astro-photometric catalog of more than 800,000 sources found in the Hubble Tarantula Treasury Project (HTTP). HTTP is a Hubble Space Telescope (HST) Treasury program designed to image the entire 30 Doradus region down to the sub-solar (~0.5 solar masses) mass regime using the Wide Field Camera 3 (WFC3) and the Advanced Camera for Surveys (ACS). We observed 30 Doradus in the near ultraviolet (F275W, F336W), optical (F555W, F658N, F775W), and near infrared (F110W, F160W) wavelengths. The stellar photometry was measured using point-spread function (PSF) fitting across all the bands simultaneously. The relative astrometric accuracy of the catalog is 0.4 mas. The astro-photometric catalog, results from artificial star experiments and the mosaics for all the filters are available for download. Color-magnitude diagrams are presented showing the spatial distributions and ages of stars within 30 Dor as well as in the surrounding fields. HTTP provides the first rich and statistically significant sample of intermediate and low mass pre-main sequence candidates and allows us to trace how star formation has been developing through the region. The depth and high spatial resolution of our analysis highlight the dual role of stellar feedback in quenching and triggering star formation on the giant HII region scale. Our results are consistent with stellar sub-clustering in a partially filled gaseous nebula that is offset towards our side of the Large Magellanic Cloud.


title: Simulating Branching Programs with Edit Distance and Friends or: A Polylog Shaved is a Lower Bound Made

abstract: A recent and active line of work achieves tight lower bounds for fundamental problems under the Strong Exponential Time Hypothesis (SETH). A celebrated result of Backurs and Indyk (STOC'15) proves that the Edit Distance of two sequences of length n cannot be computed in strongly subquadratic time under SETH. The result was extended by follow-up works to simpler looking problems like finding the Longest Common Subsequence (LCS). SETH is a very strong assumption, asserting that even linear size CNF formulas cannot be analyzed for satisfiability with an exponential speedup over exhaustive search. We consider much safer assumptions, e.g. that such a speedup is impossible for SAT on much more expressive representations, like NC circuits. Intuitively, this seems much more plausible: NC circuits can implement complex cryptographic primitives, while CNFs cannot even approximately compute an XOR of bits. Our main result is a surprising reduction from SAT on Branching Programs to fundamental problems in P like Edit Distance, LCS, and many others. Truly subquadratic algorithms for these problems therefore have consequences that we consider to be far more remarkable than merely faster CNF SAT algorithms. For example, SAT on arbitrary o(n)-depth bounded fan-in circuits (and therefore also NC-Circuit-SAT) can be solved in (2-eps)^n time. A very interesting feature of our work is that we can prove major consequences even from mildly subquadratic algorithms for Edit Distance or LCS. For example, we show that if we can shave an arbitrarily large polylog factor from n^2 for Edit Distance then NEXP does not have non-uniform NC^1 circuits. A more fine-grained examination shows that even shaving a $\log^c{n}$ factor, for a specific constant $c \approx 10^3$, already implies new circuit lower bounds.


title: Bounds on the Burning Number

abstract: Motivated by a graph theoretic process intended to measure the speed of the spread of contagion in a graph, Bonato, Janssen, and Roshanbin [Burning a Graph as a Model of Social Contagion, Lecture Notes in Computer Science 8882 (2014) 13-22] define the burning number $b(G)$ of a graph $G$ as the smallest integer $k$ for which there are vertices $x_1,\ldots,x_k$ such that for every vertex $u$ of $G$, there is some $i\in \{ 1,\ldots,k\}$ with ${\rm dist}_G(u,x_i)\leq k-i$, and ${\rm dist}_G(x_i,x_j)\geq j-i$ for every $i,j\in \{ 1,\ldots,k\}$. For a connected graph $G$ of order $n$, they prove that $b(G)\leq 2\left\lceil\sqrt{n}\right\rceil-1$, and conjecture $b(G)\leq \left\lceil\sqrt{n}\right\rceil$. We show that $b(G)\leq \sqrt{\frac{32}{19}\cdot \frac{n}{1-\epsilon}}+\sqrt{\frac{27}{19\epsilon}}$ and $b(G)\leq \sqrt{\frac{12n}{7}}+3\approx 1.309 \sqrt{n}+3$ for every connected graph $G$ of order $n$ and every $0<\epsilon<1$. For a tree $T$ of order $n$ with $n_2$ vertices of degree $2$, and $n_{\geq 3}$ vertices of degree at least $3$, we show $b(T)\leq \left\lceil\sqrt{(n+n_2)+\frac{1}{4}}+\frac{1}{2}\right\rceil$ and $b(T)\leq \left\lceil\sqrt{n}\right\rceil+n_{\geq 3}$. Furthermore, we characterize the binary trees of depth $r$ that have burning number $r+1$.


title: Vector leptoquark resolution of $R_K$ and $R_{D^{(*)}}$ puzzles

abstract: We propose that three recent anomalies in $B$ meson decays, $R_{D^{(*)}}$, $R_K$, and $P_5'$, might be explained by a single vector leptoquark weak triplet state. The constraints on the parameter space are obtained by considering $t \to b \tau^+ \nu$ data, lepton flavor universality tests in the kaon sector, bounds on $B \to K^{(*)} \bar \nu \nu$, bound on the lepton flavor violating decay $B \to K \mu \tau$, and measurements of $b \to c \mu^- \bar\nu$ decays. The presence of such vector leptoquark could be exposed in precise measurements of $t \to b \tau \nu$ and $B \to K^{(*)} \bar \nu \nu$ decays. The model also predicts approximate equality of lepton flavor universality ratios $R_{K^*}$, $R_{K}$, and suppressed branching fraction of $B_s \to \mu^+ \mu^-$.


title: Seeing relativity -- I. Ray tracing in a Schwarzschild metric to explore the maximal analytic extension of the metric and making a proper rendering of the stars

abstract: We present an implementation of a ray tracing code in the Schwarzschild metric. We aim at building a numerical code with a correct implementation of both special (aberration, amplification, Doppler) and general (deflection of light, lensing, gravitational redshift) relativistic effects so as to simulate what an observer with arbitrary velocity would see near, or possibly within, the black hole. We also pay some specific attention to perform a satisfactory rendering of stars. Using this code, we then show several unexplored features of the maximal analytical extension of the metric. In particular, we study the aspect of the second asymptotic region of the metric as seen by an observer crossing the horizon. We also address several aspects related to the white hole region (i.e., past singularity) seen both from outside the black hole, inside the future horizon and inside the past horizon, which gives rise to the most counter-intuitive effects.


title: Packings of 3D stars: Stability and structure

abstract: We describe a series of experiments involving the creation of cylindrical packings of star-shaped particles, and an exploration of the stability of these packings. The stars cover a broad range of arm sizes and frictional properties. We carried out three different kinds of experiments, all of which involve columns that are prepared by raining star particles one-by-one into hollow cylinders. As an additional part of the protocol, we sometimes vibrated the column before removing the confining cylinder. We rate stability in terms of r, the ratio of the mass of particles that fall off a pile when it collapsed, to the total particle mass. The first experiment involved the intrinsic stability of the pile when the confining cylinder was removed. The second kind of experiment involved adding a uniform load to the top of the column, and then determining the collapse properties. A third experiment involved testing stability to tipping of the piles. We find a stability diagram relating the pile height, h, vs. pile diameter, delta, where the stable and unstable regimes are separated by a boundary that is roughly a power-law in h vs. delta with an exponent that is less than one. Increasing friction and vibration both tend to stabilize piles, while increasing particle size can destabilize the system under certain conditions.


title: On the Refracted-Reflected Spectrally Negative L\'evy Processes

abstract: We study a combination of the refracted and reflected L\'evy processes. Given a spectrally negative L\'evy process and two boundaries, it is reflected at the lower boundary while, whenever it is above the upper boundary, a linear drift at a constant rate is subtracted from the increments of the process. Using the scale functions, we compute the resolvent measure, the Laplace transform of the occupation times as well as other fluctuation identities that will be useful in applied probability including insurance, queues, and inventory management.


title: Optimal inference in a class of regression models

abstract: We consider the problem of constructing confidence intervals (CIs) for a linear functional of a regression function, such as its value at a point, the regression discontinuity parameter, or a regression coefficient in a linear or partly linear regression. Our main assumption is that the regression function is known to lie in a convex function class, which covers most smoothness and/or shape assumptions used in econometrics. We derive finite-sample optimal CIs and sharp efficiency bounds under normal errors with known variance. We show that these results translate to uniform (over the function class) asymptotic results when the error distribution is not known. When the function class is centrosymmetric, these efficiency bounds imply that minimax CIs are close to efficient at smooth regression functions. This implies, in particular, that it is impossible to form CIs that are tighter using data-dependent tuning parameters, and maintain coverage over the whole function class. We specialize our results to inference on the regression discontinuity parameter, and illustrate them in simulations and an empirical application.


title: A Tensor-Train accelerated solver for integral equations in complex geometries

abstract: We present a framework using the Quantized Tensor Train (QTT) decomposition to accurately and efficiently solve volume and boundary integral equations in three dimensions. We describe how the QTT decomposition can be used as a hierarchical compression and inversion scheme for matrices arising from the discretization of integral equations. For a broad range of problems, computational and storage costs of the inversion scheme are extremely modest $O(\log N)$ and once the inverse is computed, it can be applied in $O(N \log N)$. We analyze the QTT ranks for hierarchically low rank matrices and discuss its relationship to commonly used hierarchical compression techniques such as FMM and HSS. We prove that the QTT ranks are bounded for translation-invariant systems and argue that this behavior extends to non-translation invariant volume and boundary integrals. For volume integrals, the QTT decomposition provides an efficient direct solver requiring significantly less memory compared to other fast direct solvers. We present results demonstrating the remarkable performance of the QTT-based solver when applied to both translation and non-translation invariant volume integrals in 3D. For boundary integral equations, we demonstrate that using a QTT decomposition to construct preconditioners for a Krylov subspace method leads to an efficient and robust solver with a small memory footprint. We test the QTT preconditioners in the iterative solution of an exterior elliptic boundary value problem (Laplace) formulated as a boundary integral equation in complex, multiply connected geometries.


title: BIRDNEST: Bayesian Inference for Ratings-Fraud Detection

abstract: Review fraud is a pervasive problem in online commerce, in which fraudulent sellers write or purchase fake reviews to manipulate perception of their products and services. Fake reviews are often detected based on several signs, including 1) they occur in short bursts of time; 2) fraudulent user accounts have skewed rating distributions. However, these may both be true in any given dataset. Hence, in this paper, we propose an approach for detecting fraudulent reviews which combines these 2 approaches in a principled manner, allowing successful detection even when one of these signs is not present. To combine these 2 approaches, we formulate our Bayesian Inference for Rating Data (BIRD) model, a flexible Bayesian model of user rating behavior. Based on our model we formulate a likelihood-based suspiciousness metric, Normalized Expected Surprise Total (NEST). We propose a linear-time algorithm for performing Bayesian inference using our model and computing the metric. Experiments on real data show that BIRDNEST successfully spots review fraud in large, real-world graphs: the 50 most suspicious users of the Flipkart platform flagged by our algorithm were investigated and all identified as fraudulent by domain experts at Flipkart.


abstract: We study Fourier--Mukai partners of elliptic ruled surfaces. We also describe the autoequivalence group of the derived categories of ruled surfaces with an elliptic fibration, by using \cite{Ue15}.


title: The Entropic Measure Transform

abstract: We introduce the entropic measure transform (EMT) problem for a general process and prove the existence of a unique optimal measure characterizing the solution. The density process of the optimal measure is characterized using a semimartingale BSDE under general conditions. The EMT is used to reinterpret the conditional entropic risk-measure and to obtain a convenient formula for the conditional expectation of a process which admits an affine representation under a related measure. The entropic measure transform is then used provide a new characterization of defaultable bond prices, forward prices, and futures prices when the asset is driven by a jump diffusion. The characterization of these pricing problems in terms of the EMT provides economic interpretations as a maximization of returns subject to a penalty for removing financial risk as expressed through the aggregate relative entropy. The EMT is shown to extend the optimal stochastic control characterization of default-free bond prices of Gombani and Runggaldier (Math. Financ. 23(4):659-686, 2013). These methods are illustrated numerically with an example in the defaultable bond setting.


title: EigenRec: Generalizing PureSVD for Effective and Efficient Top-N Recommendations

abstract: We introduce EigenRec; a versatile and efficient Latent-Factor framework for Top-N Recommendations that includes the well-known PureSVD algorithm as a special case. EigenRec builds a low dimensional model of an inter-item proximity matrix that combines a similarity component, with a scaling operator, designed to control the influence of the prior item popularity on the final model. Seeing PureSVD within our framework provides intuition about its inner workings, exposes its inherent limitations, and also, paves the path towards painlessly improving its recommendation performance. A comprehensive set of experiments on the MovieLens and the Yahoo datasets based on widely applied performance metrics, indicate that EigenRec outperforms several state-of-the-art algorithms, in terms of Standard and Long-Tail recommendation accuracy, exhibiting low susceptibility to sparsity, even in its most extreme manifestations -- the Cold-Start problems. At the same time EigenRec has an attractive computational profile and it can apply readily in large-scale recommendation settings.


title: Binary Locally Repairable Codes ---Sequential Repair for Multiple Erasures

abstract: Locally repairable codes (LRC) for distribute storage allow two approaches to locally repair multiple failed nodes: 1) parallel approach, by which each newcomer access a set of $r$ live nodes $(r$ is the repair locality$)$ to download data and recover the lost packet; and 2) sequential approach, by which the newcomers are properly ordered and each newcomer access a set of $r$ other nodes, which can be either a live node or a newcomer ordered before it. An $[n,k]$ linear code with locality $r$ and allows local repair for up to $t$ failed nodes by sequential approach is called an $(n,k,r,t)$-exact locally repairable code (ELRC). In this paper, we present a family of binary codes which is equivalent to the direct product of $m$ copies of the $[r+1,r]$ single-parity-check code. We prove that such codes are $(n,k,r,t)$-ELRC with $n=(r+1)^m,k=r^m$ and $t=2^m-1$, which implies that they permit local repair for up to $2^m-1$ erasures by sequential approach. Our result shows that the sequential approach has much bigger advantage than parallel approach.


title: Malthusian Locks

abstract: Applications running in modern multithreaded environments are sometimes \emph{over-threaded}. The excess threads do not improve performance, and in fact may act to degrade performance via \emph{scalability collapse}. Often, such software also has highly contended locks. We opportunistically leverage the existence of such locks by modifying the lock admission policy so as to intentionally limit the number of threads circulating over the lock in a given period. Specifically, if there are more threads circulating than are necessary to keep the lock saturated, our approach will selectively cull and passivate some of those threads. We borrow the concept of \emph{swapping} from the field of memory management and intentionally impose \emph{concurrency restriction} (CR) if a lock is oversubscribed. In the worst case CR does no harm, but it often yields performance benefits. The resultant admission order is unfair over the short term but we explicitly provide long-term fairness by periodically shifting threads between the set of passivated threads and those actively circulating. Our approach is palliative, but often effective.


title: Stochastic gradient method with accelerated stochastic dynamics

abstract: In this paper, we propose a novel technique to implement stochastic gradient methods, which are beneficial for learning from large datasets, through accelerated stochastic dynamics. A stochastic gradient method is based on mini-batch learning for reducing the computational cost when the amount of data is large. The stochasticity of the gradient can be mitigated by the injection of Gaussian noise, which yields the stochastic Langevin gradient method; this method can be used for Bayesian posterior sampling. However, the performance of the stochastic Langevin gradient method depends on the mixing rate of the stochastic dynamics. In this study, we propose violating the detailed balance condition to enhance the mixing rate. Recent studies have revealed that violating the detailed balance condition accelerates the convergence to a stationary state and reduces the correlation time between the samplings. We implement this violation of the detailed balance condition in the stochastic gradient Langevin method and test our method for a simple model to demonstrate its performance.

categories: stat.ML cond-mat.dis-nn cond-mat.stat-mech cs.CV


title: Enumeration and Random Generation of Unlabeled Classes of Graphs: A Practical Study of Cycle Pointing and the Dissymmetry Theorem

abstract: Our work studies the enumeration and random generation of unlabeled combinatorial classes of unrooted graphs. While the technique of vertex pointing provides a straightforward procedure for analyzing a labeled class of unrooted graphs by first studying its rooted counterpart, the existence of nontrivial symmetries in the unlabeled case causes this technique to break down. Instead, techniques such as the dissymmetry theorem (of Otter) and cycle pointing (of Bodirsky et al.) have emerged in the unlabeled case, with the former providing an enumeration of the class and the latter providing both an enumeration and an unbiased sampler. In this work, we extend the power of the dissymmetry theorem by showing that it in fact provides a Boltzmann sampler for the class in question. We then present an exposition of the cycle pointing technique, with a focus on the enumeration and random generation of the underlying unpointed class. Finally, we apply cycle pointing to enumerate and implement samplers for the classes of distance-hereditary graphs and three-leaf power graphs.


title: Neural Variational Inference for Text Processing

abstract: Recent advances in neural variational inference have spawned a renaissance in deep latent variable models. In this paper we introduce a generic variational inference framework for generative and conditional models of text. While traditional variational methods derive an analytic approximation for the intractable distributions over latent variables, here we construct an inference network conditioned on the discrete text input to provide the variational distribution. We validate this framework on two very different text modelling applications, generative document modelling and supervised question answering. Our neural variational document model combines a continuous stochastic document representation with a bag-of-words generative model and achieves the lowest reported perplexities on two standard test corpora. The neural answer selection model employs a stochastic representation layer within an attention mechanism to extract the semantics between a question and answer pair. On two question answering benchmarks this model exceeds all previous published benchmarks.


title: Combining $Q_T$ and small-$x$ resummations

abstract: We analyze transverse momentum ($Q_T$) resummation of a colorless final state, e.g. Higgs production in gluon fusion or the production of a lepton pair via the Drell-Yan mechanism, in the limit where the invariant mass of the final state is much less then the center-of-mass energy, i.e. $Q^2\ll s$. We show how the traditional resummation of logarithms of $Q_T/Q$ can be supplemented with the resummation of the leading logarithmic contributions at small $x=Q^2/s$ and we compute the necessary ingredients to perform such joint resummation.


title: A Hierarchical Deep Temporal Model for Group Activity Recognition

abstract: In group activity recognition, the temporal dynamics of the whole activity can be inferred based on the dynamics of the individual people representing the activity. We build a deep model to capture these dynamics based on LSTM (long-short term memory) models. To make use of these ob- servations, we present a 2-stage deep temporal model for the group activity recognition problem. In our model, a LSTM model is designed to represent action dynamics of in- dividual people in a sequence and another LSTM model is designed to aggregate human-level information for whole activity understanding. We evaluate our model over two datasets: the collective activity dataset and a new volley- ball dataset. Experimental results demonstrate that our proposed model improves group activity recognition perfor- mance with compared to baseline methods.