title: Neutral Carbon Emission in luminous infrared galaxies The \CI\ Lines as Total Molecular Gas Tracers

abstract: We present a statistical study on the [C I]($^{3} \rm P_{1} \rightarrow {\rm ^3 P}_{0}$), [C I] ($^{3} \rm P_{2} \rightarrow {\rm ^3 P}_{1}$) lines (hereafter [C I] (1$-$0) and [C I] (2$-$1), respectively) and the CO (1$-$0) line for a sample of (ultra)luminous infrared galaxies [(U)LIRGs]. We explore the correlations between the luminosities of CO (1$-$0) and [C I] lines, and find that $L'_\mathrm{CO(1-0)}$ correlates almost linearly with both $L'_ \mathrm{[CI](1-0)}$ and $L'_\mathrm{[CI](2-1)}$, suggesting that [C I] lines can trace total molecular gas mass at least for (U)LIRGs. We also investigate the dependence of $L'_\mathrm{[CI](1-0)}$/$L'_\mathrm{CO(1-0)}$, $L'_\mathrm{[CI](2-1)}$/$L'_\mathrm{CO(1-0)}$ and $L'_\mathrm{[CI](2-1)}$/$L'_\mathrm{[CI](1-0)}$ on the far-infrared color of 60-to-100 $\mu$m, and find non-correlation, a weak correlation and a modest correlation, respectively. Under the assumption that these two carbon transitions are optically thin, we further calculate the [C I] line excitation temperatures, atomic carbon masses, and the mean [C I] line flux-to-H$_2$ mass conversion factors for our sample. The resulting $\mathrm{H_2}$ masses using these [C I]-based conversion factors roughly agree with those derived from $L'_\mathrm{CO(1-0)}$ and CO-to-H$_2$ conversion factor.


title: Search for supersymmetry in multijet events with missing transverse momentum in proton-proton collisions at 13 TeV

abstract: A search for supersymmetry is presented based on multijet events with large missing transverse momentum produced in proton-proton collisions at a center-of-mass energy of sqrt(s) = 13 TeV. The data, corresponding to an integrated luminosity of 35.9 inverse femtobarns, were collected with the CMS detector at the CERN LHC in 2016. The analysis utilizes four-dimensional exclusive search regions defined in terms of the number of jets, the number of tagged bottom quark jets, the scalar sum of jet transverse momenta, and the magnitude of the vector sum of jet transverse momenta. No evidence for a significant excess of events is observed relative to the expectation from the standard model. Limits on the cross sections for the pair production of gluinos and squarks are derived in the context of simplified models. Assuming the lightest supersymmetric particle to be a weakly interacting neutralino, 95% confidence level lower limits on the gluino mass as large as 1800 to 1960 GeV are derived, and on the squark mass as large as 960 to 1390 GeV, depending on the production and decay scenario.


title: Mass and p-factor of the type II Cepheid OGLE-LMC-T2CEP-098 in a binary system

abstract: We present the results of a study of the type II Cepheid ($P_{puls} = 4.974 d$) in the eclipsing binary system OGLE-LMC-T2CEP-098 ($P_{orb} = 397.2 d$). The Cepheid belongs to the peculiar W Vir group, for which the evolutionary status is virtually unknown. It is the first single-lined system with a pulsating component analyzed using the method developed by Pilecki et al. (2013). We show that the presence of a pulsator makes it possible to derive accurate physical parameters of the stars even if radial velocities can be measured for only one of the components. We have used four different methods to limit and estimate the physical parameters, eventually obtaining precise results by combining pulsation theory with the spectroscopic and photometric solutions. The Cepheid radius, mass and temperature are $25.3 \pm 0.2 R_\odot$, $1.51 \pm 0.09 M_\odot$ and $5300 \pm 100 K$, respectively, while its companion has similar size ($26.3 R_\odot$), but is more massive ($6.8 M_\odot$) and hotter ($9500 K$). Our best estimate for the p-factor of the Cepheid is $1.30 \pm 0.03$. The mass, position on the period-luminosity diagram, and pulsation amplitude indicate that the pulsating component is very similar to the Anomalous Cepheids, although it has a much longer period and is redder in color. The very unusual combination of the components suggest that the system has passed through a mass transfer phase in its evolution. More complicated internal structure would then explain its peculiarity.


title: Extreme submillimetre starburst galaxies

abstract: We use two catalogues, a Herschel catalogue selected at 500 mu (HerMES) and an IRAS catalogue selected at 60 mu (RIFSCz), to contrast the sky at these two wavelengths. Both surveys demonstrate the existence of extreme starbursts, with star-formation rates (SFRs) > 5000 Msun/yr. The maximum intrinsic star-formation rate appears to be ~30,000 Msun/yr. The sources with apparent SFR estimates higher than this are in all cases either lensed systems, blazars, or erroneous photometric redshifts. At redshifts of 3 to 5, the time-scale for the Herschel galaxies to make their current mass of stars at their present rate of formation ~ 10^8 yrs, so these galaxies are making a significant fraction of their stars in the current star-formation episode. Using dust mass as a proxy for gas mass, the Herschel galaxies at redshift 3 to 5 have gas masses comparable to their mass in stars. Of the 38 extreme starbursts in our Herschel survey for which we have more complete SED information, over 50% show evidence for QSO-like optical emission, or exhibit AGN dust tori in the mid-infrared SEDs. In all cases however the infrared luminosity is dominated by a starburst component. We derive a mean covering factor for AGN dust as a function of redshift and derive black hole masses and black hole accretion rates. There is a universal ratio of black-hole mass to stellar mass, ~ 10^{-3}, driven by the strong period of star-formation and black-hole growth at z = 1-5.


title: Tight bounds on the coefficients of partition functions via stability

abstract: Partition functions arise in statistical physics and probability theory as the normalizing constant of Gibbs measures and in combinatorics and graph theory as graph polynomials. For instance the partition functions of the hard-core model and monomer-dimer model are the independence and matching polynomials respectively. We show how stability results follow naturally from the recently developed occupancy method for maximizing and minimizing physical observables over classes of regular graphs, and then show these stability results can be used to obtain tight extremal bounds on the individual coefficients of the corresponding partition functions. As applications, we prove new bounds on the number of independent sets and matchings of a given size in regular graphs. For large enough graphs and almost all sizes, the bounds are tight and confirm the Upper Matching Conjecture of Friedland, Krop, and Markstr\"om and a conjecture of Kahn on independent sets for a wide range of parameters. Additionally we prove tight bounds on the number of $q$-colorings of cubic graphs with a given number of monochromatic edges, and tight bounds on the number of independent sets of a given size in cubic graphs of girth at least $5$.


title: Thinking Fast and Slow: Optimization Decomposition Across Timescales

abstract: Many real-world control systems, such as the smart grid and human sensorimotor control systems, have decentralized components that react quickly using local information and centralized components that react slowly using a more global view. This paper seeks to provide a theoretical framework for how to design controllers that are decomposed across timescales in this way. The framework is analogous to how the network utility maximization framework uses optimization decomposition to distribute a global control problem across independent controllers, each of which solves a local problem; except our goal is to decompose a global problem temporally, extracting a timescale separation. Our results highlight that decomposition of a multi-timescale controller into a fast timescale, reactive controller and a slow timescale, predictive controller can be near-optimal in a strong sense. In particular, we exhibit such a design, named Multi-timescale Reflexive Predictive Control (MRPC), which maintains a per-timestep cost within a constant factor of the offline optimal in an adversarial setting.


title: A new approach to exact solutions construction in scalar cosmology with a Gauss-Bonnet term

abstract: We study the cosmological model based on Einstein-Gauss-Bonnet gravity with non-minimal coupling of a scalar field to a Gauss-Bonnet term in 4D Friedmann universe. We show how constructing the exact solutions by the method based on a confrontation of the Hubble parameter in the model under consideration with that in a standard scalar field inflationary cosmology.


title: Finding Exogenous Variation in Data

abstract: We reconsider the classic problem of recovering exogenous variation from an endogenous regressor. Two-stage least squares recovers exogenous variation through presuming the existence of an instrumental variable. We rely instead on the assumption that the regressor is a mixture of exogenous and endogenous observations--say as the result of temporary natural experiments. With this assumption, we propose an alternative two-stage method based on nonparametrically estimating a mixture model to recover a subset of the exogenous observations. We demonstrate that our method recovers exogenous observations in simulation and can be used to find pricing experiments hidden in grocery store scanner data.


title: Properly embedded minimal annuli in $\mathbb{H}^2 \times \mathbb{R}$

abstract: In this paper we study the moduli space of properly Alexandrov-embedded, minimal annuli in $\mathbb{H}^2 \times \mathbb{R}$ with horizontal ends. We say that the ends are horizontal when they are graphs of $\mathcal{C}^{2, \alpha}$ functions over $\partial_\infty \mathbb{H}^2$. Contrary to expectation, we show that one can not fully prescribe the two boundary curves at infinity, but rather, one can prescribe the bottom curve, but the top curve only up to a translation and a tilt, along with the position of the neck and the vertical flux of the annulus. We also prove general existence theorems for minimal annuli with discrete groups of symmetries.


title: On variance estimation for generalizing from a trial to a target population

abstract: Randomized controlled trials (RCTs) provide strong internal validity compared with observational studies. However, selection bias threatens the external validity of randomized trials. Thus, RCT results may not apply to either broad public policy populations or narrow populations, such as specific insurance pools. Some researchers use propensity scores (PSs) to generalize results from an RCT to a target population. In this scenario, a PS is defined as the probability of participating in the trial conditioning on observed covariates. We study a model-free inverse probability weighted estimator (IPWE) of the average treatment effect in a target population with data from a randomized trial. We present variance estimators and compare the performance of our method with that of model-based approaches. We examine the robustness of the model-free estimators to heterogeneous treatment effects.


title: FWDA: a Fast Wishart Discriminant Analysis with its Application to Electronic Health Records Data Classification

abstract: Linear Discriminant Analysis (LDA) on Electronic Health Records (EHR) data is widely-used for early detection of diseases. Classical LDA for EHR data classification, however, suffers from two handicaps: the ill-posed estimation of LDA parameters (e.g., covariance matrix), and the "linear inseparability" of EHR data. To handle these two issues, in this paper, we propose a novel classifier FWDA -- Fast Wishart Discriminant Analysis, that makes predictions in an ensemble way. Specifically, FWDA first surrogates the distribution of inverse covariance matrices using a Wishart distribution estimated from the training data, then "weighted-averages" the classification results of multiple LDA classifiers parameterized by the sampled inverse covariance matrices via a Bayesian Voting scheme. The weights for voting are optimally updated to adapt each new input data, so as to enable the nonlinear classification. Theoretical analysis indicates that FWDA possesses a fast convergence rate and a robust performance on high dimensional data. Extensive experiments on large-scale EHR dataset show that our approach outperforms state-of-the-art algorithms by a large margin.


title: Concave Flow on Small Depth Directed Networks

abstract: Small depth networks arise in a variety of network related applications, often in the form of maximum flow and maximum weighted matching. Recent works have generalized such methods to include costs arising from concave functions. In this paper we give an algorithm that takes a depth $D$ network and strictly increasing concave weight functions of flows on the edges and computes a $(1 - \epsilon)$-approximation to the maximum weight flow in time $mD \epsilon^{-1}$ times an overhead that is logarithmic in the various numerical parameters related to the magnitudes of gradients and capacities. Our approach is based on extending the scaling algorithm for approximate maximum weighted matchings by [Duan-Pettie JACM`14] to the setting of small depth networks, and then generalizing it to concave functions. In this more restricted setting of linear weights in the range $[w_{\min}, w_{\max}]$, it produces a $(1 - \epsilon)$-approximation in time $O(mD \epsilon^{-1} \log( w_{\max} /w_{\min}))$. The algorithm combines a variety of tools and provides a unified approach towards several problems involving small depth networks.


title: The Gordian distance of handlebody-knots and Alexander biquandle colorings

abstract: We give lower bounds for the Gordian distance and the unknotting number of handlebody-knots by using Alexander biquandle colorings. We construct handlebody-knots with Gordian distance $n$ and unknotting number $n$ for any positive integer $n$.


title: Arabidopsis roots segmentation based on morphological operations and CRFs

abstract: Arabidopsis thaliana is a plant species widely utilized by scientists to estimate the impact of genetic differences in root morphological features. For this purpose, images of this plant after genetic modifications are taken to study differences in the root architecture. This task requires manual segmentations of radicular structures, although this is a particularly tedious and time-consuming labor. In this work, we present an unsupervised method for Arabidopsis thaliana root segmentation based on morphological operations and fully-connected Conditional Random Fields. Although other approaches have been proposed to this purpose, all of them are based on more complex and expensive imaging modalities. Our results prove that our method can be easily applied over images taken using conventional scanners, with a minor user intervention. A first data set, our results and a fully open source implementation are available online.


title: Refinement of Bolide Characteristics from Infrasound measurements

abstract: We have detected and performed signal measurements on 78 individual bolide events as recorded at 179 infrasound stations between 2006 and 2015. We compared period-yield relations with AFTAC nuclear period-yield data, finding these to be similar with a slight offset. Scatter in period measurements for individual bolide is found to be caused in part by station noise levels and by attenuation effects with range. No correlation was found between the infrasound signal period and any of bolide height at peak brightness, entry speed or impact angle. We examined in detail three well constrained bolides having energy deposition curves, known trajectories and infrasound detections finding some evidence at shorter ranges that a component of station period scatter is due to varying source heights sampled by different stations. However, for longer-range stations in these three case studies, we were not able to assign unique source heights using raytracing due to large uncertainties in atmospheric conditions. Our results suggest that while source height contributes to the observed variance in infrasound signal periods from a given bolide, range and station noise play a larger role.


title: In-plane spectrum in superlattices

abstract: We show that the existing theory does not give correct in-plane spectrum of superlattices at small in-plane momentum. Magneto-absorption experiments demonstrate that the energy range of the parabolic region of the spectrum near the electron subband bottom is by the order of magnitude lower than the value predicted by the traditional approach. We developed a modified theory according to which the energy range of the parabolic region and carrier in-plane effective masses are determined by the effective bandgap of the superlattice rather than by the bulk bandgaps of the superlattice layers. The results of the new theory are consistent with the experiment.


title: Ribbon graphs and the fundamental group of surfaces

abstract: In the present work we are going to give a formal exposition of the ribbon graphs topic based on notes of Labourie \cite{Lab}, since is difficult to find as such in the literature.


title: Gentle heating by mixing in cooling flow clusters

abstract: We analyze three-dimensional hydrodynamical simulations of the interaction of jets and the bubbles they inflate with the intra-cluster medium (ICM), and show that the heating of the ICM by mixing hot bubble gas with the ICM operates over tens of millions of years, and hence can smooth the sporadic activity of the jets. The inflation process of hot bubbles by propagating jets forms many vortices, and these vortices mix the hot bubble gas with the ICM. The mixing, hence the heating of the ICM, starts immediately after the jets are launched, but continues for tens of millions of years. We suggest that the smoothing of the active galactic nucleus (AGN) sporadic activity by the long-lived vortices accounts for the recent finding of a gentle energy coupling between AGN heating and the ICM.


title: Limitations on Transversal Computation through Quantum Homomorphic Encryption

abstract: Transversality is a simple and effective method for implementing quantum computation fault-tolerantly. However, no quantum error-correcting code (QECC) can transversally implement a quantum universal gate set (Eastin and Knill, Phys. Rev. Lett., 102, 110502). Since reversible classical computation is often a dominating part of useful quantum computation, whether or not it can be implemented transversally is an important open problem. We show that, other than a small set of non-additive codes that we cannot rule out, no binary QECC can transversally implement a classical reversible universal gate set. In particular, no such QECC can implement the Toffoli gate transversally. We prove our result by constructing an information theoretically secure (but inefficient) quantum homomorphic encryption (ITS-QHE) scheme inspired by Ouyang et al. (arXiv:1508.00938). Homomorphic encryption allows the implementation of certain functions directly on encrypted data, i.e. homomorphically. Our scheme builds on almost any QECC, and implements that code's transversal gate set homomorphically. We observe a restriction imposed by Nayak's bound (FOCS 1999) on ITS-QHE, implying that any ITS quantum fully homomorphic scheme (ITS-QFHE) implementing the full set of classical reversible functions must be highly inefficient. While our scheme incurs exponential overhead, any such QECC implementing Toffoli transversally would still violate this lower bound through our scheme.


title: Invariant Measures for TASEP with a Slow Bond

abstract: Totally Asymmetric Simple Exclusion Process (TASEP) on $\mathbb{Z}$ is one of the classical exactly solvable models in the KPZ universality class. We study the "slow bond" model, where TASEP on $\mathbb{Z}$ is imputed with a slow bond at the origin. The slow bond increases the particle density immediately to its left and decreases the particle density immediately to its right. Whether or not this effect is detectable in the macroscopic current started from the step initial condition has attracted much interest over the years and this question was settled recently (Basu, Sidoravicius, Sly (2014)) by showing that the current is reduced even for arbitrarily small strength of the defect. Following non-rigorous physics arguments (Janowsky, Lebowitz (1992, 1994)) and some unpublished works by Bramson, a conjectural description of properties of invariant measures of TASEP with a slow bond at the origin was provided in Liggett's 1999 book. We establish Liggett's conjectures and in particular show that TASEP with a slow bond at the origin, started from step initial condition, converges in law to an invariant measure that is asymptotically close to product measures with different densities far away from the origin towards left and right. Our proof exploits the correspondence between TASEP and the last passage percolation on $\mathbb{Z}^2$ with exponential weights and uses the understanding of geometry of maximal paths in those models.

categories: math.PR cond-mat.stat-mech math-ph math.MP


title: A data-driven analysis of the heavy quark transport coefficient

abstract: Using a Bayesian model-to-data analysis, we estimate the temperature dependence of the heavy quark diffusion coefficients by calibrating to the experimental data of $D$-meson $R_{\mathrm{AA}}$ and $v_2$ in AuAu collisions ($\sqrt{s_{NN}}=200$ GeV) and PbPb collisions ($\sqrt{s_{NN}}=2.76$ TeV)~\cite{Xie:2016iwq}. The spatial diffusion coefficient $D_s2\pi T$ is found to be mostly constraint around $(1.3-1.5) T_c$ and is compatible with lattice QCD calculations. We demonstrate the capability of our improved Langevin model to simultaneously describe the $R_{\mathrm{AA}}$ and $v_2$ at both RHIC and the LHC energies, as well as the feasibility to apply a Bayesian analysis to quantitatively study the heavy flavor transport in heavy-ion collisions.


title: Existence results for minimizers of parametric elliptic functionals

abstract: We prove a compactness principle for the anisotropic formulation of the Plateau problem in any codimension, in the same spirit of the previous works of the authors \cite{DelGhiMag,DePDeRGhi,DeLDeRGhi16}. In particular, we perform a new strategy for the proof of the rectifiability of the minimal set, based on the new anisotropic counterpart of the Allard rectifiability theorem proved by the authors in \cite{DePDeRGhi2}. As a consequence we provide a new proof of Reifenberg existence theorem.


title: Variation of ionic conductivity in a plastic-crystalline mixture

abstract: Ionically-conducting plastic crystals are possible candidates for solid-state electrolytes in energy-storage devices. Interestingly, the admixture of larger molecules to the most prominent molecular PC electrolyte, succinonitrile, was shown to drastically enhance its ionic conductivity. Therefore, binary mixtures seem to be a promising way to tune the conductivity of such solid-state electrolytes. However, to elucidate the general mechanisms of ionic charge transport in plastic crystals and the influence of mixing, a much broader data base is needed. In the present work, we investigate mixtures of two well-known plastic-crystalline systems, cyclohexanol and cyclooctanol, to which 1 mol% of Li ions were added. Applying differential scanning calorimetry and dielectric spectroscopy, we present a thorough investigation of the phase behavior and the ionic and dipolar dynamics of this system. All mixtures reveal plastic-crystalline phases with corresponding orientational glass-transitions. Moreover, their conductivity seems to be dominated by the "revolving-door" mechanism, implying a close coupling between the ionic translational and the molecular reorientational dynamics of the surrounding plastic-crystalline matrix. In contrast to succinonitrile-based mixtures, there is no strong variation of this coupling with the mixing ratio.


title: Big-bang nucleosynthesis constraints on dark-matter sectors revisited

abstract: We extend previous studies of big bang nucleosynthesis, with the assumption that ordinary matter and dark matter sectors are entangled through the number of degrees of freedom entering the Friedmann equations. This conjecture allows us to find a relation between the temperatures of the weakly interacting matter and dark-matter sectors. The constraints imposed by observations are studied by comparison with calculations of big bang nucleosynthesis for the abundance of light elements. It is found that these constraints are compatible with cold dark matter and a wide range number of dark sectors.


title: SfM-Net: Learning of Structure and Motion from Video

abstract: We propose SfM-Net, a geometry-aware neural network for motion estimation in videos that decomposes frame-to-frame pixel motion in terms of scene and object depth, camera motion and 3D object rotations and translations. Given a sequence of frames, SfM-Net predicts depth, segmentation, camera and rigid object motions, converts those into a dense frame-to-frame motion field (optical flow), differentiably warps frames in time to match pixels and back-propagates. The model can be trained with various degrees of supervision: 1) self-supervised by the re-projection photometric error (completely unsupervised), 2) supervised by ego-motion (camera motion), or 3) supervised by depth (e.g., as provided by RGBD sensors). SfM-Net extracts meaningful depth estimates and successfully estimates frame-to-frame camera rotations and translations. It often successfully segments the moving objects in the scene, even though such supervision is never provided.


title: Two famous results of Einstein derived from the Jarzynski equality

abstract: The Jarzynski equality (JE) is a remarkable statement relating transient irreversible processes to infinite-time free energy differences. Although twenty years old, the JE remains unfamiliar to many; nevertheless it is a robust and powerful law. We examine two of Einstein's most simple and well-known discoveries, one classical and one quantum, and show how each of these follows from the JE. Our first example is Einstein's relation between the drag and diffusion coefficients of a particle in Brownian motion. In this context we encounter a paradox in the macroscopic limit of the JE which is fascinating, but also warns us against using the JE too freely outside of the microscopic domain. Our second example is the equality of Einstein's $B$ coefficients for absorption and stimulated emission of quanta. Here resonant light does irreversible work on a sample, and the argument differs from Einstein's equilibrium reasoning using the Planck black-body spectrum. We round out our examples with a brief derivation and discussion of Jarzynski's remarkable equality.


title: A Note on Experiments and Software For Multidimensional Order Statistics

abstract: In this note we describe experiments on an implementation of two methods proposed in the literature for computing regions that correspond to a notion of order statistics for multidimensional data. Our implementation, which works for any dimension greater than one, is the only that we know of to be publicly available. Experiments run using the software confirm that half-space peeling generally gives better results than directly peeling convex hulls, but at a computational cost.


title: A decentralized proximal-gradient method with network independent step-sizes and separated convergence rates

abstract: This paper proposes a novel proximal-gradient algorithm for a decentralized optimization problem with a composite objective containing smooth and non-smooth terms. Specifically, the smooth and nonsmooth terms are dealt with by gradient and proximal updates, respectively. The proposed algorithm is closely related to a previous algorithm, PG-EXTRA \cite{shi2015proximal}, but has a few advantages. First of all, agents use uncoordinated step-sizes, and the stable upper bounds on step-sizes are independent of network topologies. The step-sizes depend on local objective functions, and they can be as large as those of the gradient descent. Secondly, for the special case without non-smooth terms, linear convergence can be achieved under the strong convexity assumption. The dependence of the convergence rate on the objective functions and the network are separated, and the convergence rate of the new algorithm is as good as one of the two convergence rates that match the typical rates for the general gradient descent and the consensus averaging. We provide numerical experiments to demonstrate the efficacy of the introduced algorithm and validate our theoretical discoveries.

categories: math.OC cs.DC cs.LG cs.NA math.NA stat.ML


title: Pulsations of intermediate-mass stars on the asymptotic giant branch

abstract: Evolutionary tracks from the zero age main sequence to the asymptotic giant branch were computed for stars with initial masses 2M_\odot <= Mzams <= 5M_\odot and metallicity Z=0.02. Some models of evolutionary sequences were used as initial conditions for equations of radiation hydrodynamics and turbulent convection describing radial stellar pulsations. The early asymptotic giant branch stars are shown to pulsate in the fundamental mode with periods from 30 to 400 day. The rate of period change gradually increases as the star evolves but is too small to be detected (d ln P/dt < 1.e-5 yr^-1). Pulsation properties of thermally pulsing AGB stars are investigated on time intervals comprising 17 thermal pulses for evolutionary sequences with intial masses Mzams=2M_\odot and 3M_\odot and 6 thermal pulses for Mzams=4M_\odot and 5M_\odot. TP-AGB stars with initial masses Mzams<=3M_\odot pulsate either in the fundamental mode or in the first overtone, whereas more massive red giants (Mzams>=4M_\odot) pulsate in the fundamental mode with periods as long as 1.e3 day. Most rapid pulsation period change with rate -0.02 < d ln P/dt < -0.01 yr^-1 occurs during decrease of the surface luminosity after the maximum helium luminosity. The rate of subsequent increase of the period is d ln P/dt <= 5.e-3 yr^-1.


title: Hand Keypoint Detection in Single Images using Multiview Bootstrapping

abstract: We present an approach that uses a multi-camera system to train fine-grained detectors for keypoints that are prone to occlusion, such as the joints of a hand. We call this procedure multiview bootstrapping: first, an initial keypoint detector is used to produce noisy labels in multiple views of the hand. The noisy detections are then triangulated in 3D using multiview geometry or marked as outliers. Finally, the reprojected triangulations are used as new labeled training data to improve the detector. We repeat this process, generating more labeled data in each iteration. We derive a result analytically relating the minimum number of views to achieve target true and false positive rates for a given detector. The method is used to train a hand keypoint detector for single images. The resulting keypoint detector runs in realtime on RGB images and has accuracy comparable to methods that use depth sensors. The single view detector, triangulated over multiple views, enables 3D markerless hand motion capture with complex object interactions.


title: Sparse domination of Hilbert transforms along curves

abstract: We obtain sharp sparse bounds for Hilbert transforms along curves in $\mathbb{R}^n$, and derive as corollaries weighted norm inequalities for such operators. The curves that we consider include monomial curves and arbitrary $C^n$ curves with nonvanishing torsion.


title: $R$-triviality of some exceptional groups

abstract: The main aim of this paper is to prove $R$-triviality for simple, simply connected algebraic groups with Tits index $E_{8,2}^{78}$ or $E_{7,1}^{78}$, defined over a field $k$ of arbitrary characteristic. Let $G$ be such a group. We prove that there exists a quadratic extension $K$ of $k$ such that $G$ is $R$-trivial over $K$, i.e., for any extension $F$ of $K$, $G(F)/R=\{1\}$, where $G(F)/R$ denotes the group of $R$-equivalence classes in $G(F)$, in the sense of Manin (see \cite{M}). As a consequence, it follows that the variety $G$ is retract $K$-rational and that the Kneser-Tits conjecture holds for these groups over $K$. Moreover, $G(L)$ is projectively simple as an abstract group for any field extension $L$ of $K$. In their monograph (\cite{TW}) J. Tits and Richard Weiss conjectured that for an Albert division algebra $A$ over a field $k$, its structure group $Str(A)$ is generated by scalar homotheties and its $U$-operators. This is known to be equivalent to the Kneser-Tits conjecture for groups with Tits index $E_{8,2}^{78}$. We settle this conjecture for Albert division algebras which are first constructions, in affirmative. These results are obtained as corollaries to the main result, which shows that if $A$ is an Albert division algebra which is a first construction and $\Gamma$ its structure group, i.e., the algebraic group of the norm similarities of $A$, then $\Gamma(F)/R=\{1\}$ for any field extension $F$ of $k$, i.e., $\Gamma$ is $R$-trivial.


title: Discrete Symmetries of Calabi-Yau Hypersurfaces in Toric Four-Folds

abstract: We analyze freely-acting discrete symmetries of Calabi-Yau three-folds defined as hypersurfaces in ambient toric four-folds. An algorithm which allows the systematic classification of such symmetries which are linearly realised on the toric ambient space is devised. This algorithm is applied to all Calabi-Yau manifolds with $h^{1,1}(X)\leq 3$ obtained by triangulation from the Kreuzer-Skarke list, a list of some $350$ manifolds. All previously known freely-acting symmetries on these manifolds are correctly reproduced and we find five manifolds with freely-acting symmetries. These include a single new example, a manifold with a $\mathbb{Z}_2\times\mathbb{Z}_2$ symmetry where only one of the $\mathbb{Z}_2$ factors was previously known. In addition, a new freely-acting $\mathbb{Z}_2$ symmetry is constructed for a manifold with $h^{1,1}(X)=6$. While our results show that there are more freely-acting symmetries within the Kreuzer-Skarke set than previously known, it appears that such symmetries are relatively rare.


title: Unsupervised Learning of Depth and Ego-Motion from Video

abstract: We present an unsupervised learning framework for the task of monocular depth and camera motion estimation from unstructured video sequences. We achieve this by simultaneously training depth and camera pose estimation networks using the task of view synthesis as the supervisory signal. The networks are thus coupled via the view synthesis objective during training, but can be applied independently at test time. Empirical evaluation on the KITTI dataset demonstrates the effectiveness of our approach: 1) monocular depth performing comparably with supervised methods that use either ground-truth pose or depth for training, and 2) pose estimation performing favorably with established SLAM systems under comparable input settings.


title: Disaggregating Input--Output Tables by the Multidimensional RAS Method: A Case Study of the Czech Republic

abstract: The RAS method is an iterative procedure that bi-proportionally scales an input-output table to be consistent with given row and column sums. It can be used to disaggregate an annual national table to more detailed tables, such as regional, quarterly, and domestic/imported tables. However, the regular two-dimensional RAS method does not ensure the consistency of the disaggregated tables with the original table. For this problem, we use the multidimensional RAS method, which besides input and output totals, also ensures regional, quarterly, and domestic/imported totals. Our analysis of Czech industries shows that the multidimensional RAS method increases the accuracy of table estimation as well as the accuracy of the Leontief inverse, the quarterly value added, and (to some degree) the regional Isard's model. We also rigorously demonstrate its relation to the cross-entropy model.


title: Supercoiling DNA locates mismatches

abstract: We present a method of detecting sequence defects by supercoiling DNA with magnetic tweezers. The method is sensitive to a single mismatched base pair in a DNA sequence of several thousand base pairs. We systematically compare DNA molecules with 0 to 16 adjacent mismatches at 1 M monovalent salt and 3.5 pN force and show that, under these conditions, a single plectoneme forms and is stably pinned at the defect. We use these measurements to estimate the energy and degree of end-loop kinking at defects. From this, we calculate the relative probability of plectoneme pinning at the mismatch under physiologically relevant conditions. Based on this estimate, we propose that DNA supercoiling could contribute to mismatch and damage sensing in vivo.


title: Introspective Classification with Convolutional Nets

abstract: We propose introspective convolutional networks (ICN) that emphasize the importance of having convolutional neural networks empowered with generative capabilities. We employ a reclassification-by-synthesis algorithm to perform training using a formulation stemmed from the Bayes theory. Our ICN tries to iteratively: (1) synthesize pseudo-negative samples; and (2) enhance itself by improving the classification. The single CNN classifier learned is at the same time generative --- being able to directly synthesize new samples within its own discriminative model. We conduct experiments on benchmark datasets including MNIST, CIFAR-10, and SVHN using state-of-the-art CNN architectures, and observe improved classification results.


title: $A$-dependence of quasielastic charged-current neutrino-nucleus cross sections

abstract: Background: 12C has been and is still widely used in neutrino-nucleus scattering and oscillation experiments. More recently, 40Ar has emerged as an important nuclear target for current and future experiments. Liquid argon time projection chambers (LArTPCs) possess various advantages in measuring electroweak neutrino-nucleus cross sections. Concurrent theoretical research is an evident necessity. Purpose: 40Ar is larger than 12C, and one expects nuclear effects to play a bigger role in reactions. We present inclusive differential and total cross section results for charged-current neutrino scattering on 40Ar and perform a comparison with 12C, 16O and 56Fe targets, to find out about the A-dependent behavior of model predictions. Method: Our model starts off with a Hartree-Fock description of the nucleus, with the nucleons interacting through a mean field generated by an effective Skyrme force. Long-range correlations are introduced by means of a continuum random phase approximation (CRPA) approach. Further methods to improve the accuracy of model predictions are also incorporated in the calculations. Results: We present calculations for 12C, 16O, 40Ar and 56Fe, showcasing differential cross sections over a broad range of kinematic values in the quasielastic regime. We furthermore show flux-folded results for 40Ar and we discuss the differences between nuclear responses. Conclusions: At low incoming energies and forward scattering we identify an enhancement in the 40Ar cross section compared to 12C, as well as in the high $\omega$ (low $T_\mu$ ) region across the entire studied $E_\nu$ range. The contribution to the folded cross section of the reaction strength at values of $\omega$ lower than 50 MeV for forward scattering is sizeable.


title: Skyrmion defects and competing singlet orders in a half-filled antiferromagnetic Kondo-Heisenberg model on the honeycomb lattice

abstract: Due to the interaction between topological defects of an order parameter and underlying fermions, the defects can possess induced fermion numbers, leading to several exotic phenomena of fundamental importance to both condensed matter and high energy physics. One of the intriguing outcome of induced fermion number is the presence of fluctuating competing orders inside the core of topological defect. In this regard, the interaction between fermions and skyrmion excitations of antiferromagnetic phase can have important consequence for understanding the global phase diagrams of many condensed matter systems where antiferromagnetism and several singlet orders compete. We critically investigate the relation between fluctuating competing orders and skyrmion excitations of the antiferromagnetic insulating phase of a half-filled Kondo-Heisenberg model on honeycomb lattice. By combining analytical and numerical methods we obtain exact eigenstates of underlying Dirac fermions in the presence of a single skyrmion configuration, which are used for computing induced chiral charge. Additionally, by employing this nonperturbative eigenbasis we calculate the susceptibilities of different translational symmetry breaking charge, bond and current density wave orders and translational symmetry preserving Kondo singlet formation. Based on the computed susceptibilities we establish spin Peierls and Kondo singlets as dominant competing orders of antiferromagnetism. We show favorable agreement between our findings and field theoretic predictions based on perturbative gradient expansion scheme which crucially relies on adiabatic principle and plane wave eigenstates for Dirac fermions. The methodology developed here can be applied to many other correlated systems supporting competition between spin-triplet and spin-singlet orders in both lower and higher spatial dimensions.


title: Notes on $G_2$: The Lie algebra and the Lie group

abstract: These notes have been prepared for the Workshop on "(Non)-existence of complex structures on $\mathbb{S}^6$", to be celebrated in Marburg in March, 2017. The material is not intended to be original. It contains a survey about the smallest of the exceptional Lie groups: $G_2$, its definition and different characterizations joint with its relationship with $\mathbb{S}^6$ and with $\mathbb{S}^7$. With the exception of the summary of the Killing-Cartan classification, this survey is self-contained, and all the proofs are given, mainly following linear algebra arguments. Although these proofs are well-known, they are spread and some of them are difficult to find. The approach is algebraical, working at the Lie algebra level most of times. We analyze the complex Lie algebra (and group) of type $G_2$ as well as the two real Lie algebras of type $G_2$, the split and the compact one. Octonions will appear, but it is not the starting point. Also, 3-forms approach and spinorial approach are viewed and related.


title: Introspective Generative Modeling: Decide Discriminatively

abstract: We study unsupervised learning by developing introspective generative modeling (IGM) that attains a generator using progressively learned deep convolutional neural networks. The generator is itself a discriminator, capable of introspection: being able to self-evaluate the difference between its generated samples and the given training data. When followed by repeated discriminative learning, desirable properties of modern discriminative classifiers are directly inherited by the generator. IGM learns a cascade of CNN classifiers using a synthesis-by-classification algorithm. In the experiments, we observe encouraging results on a number of applications including texture modeling, artistic style transferring, face modeling, and semi-supervised learning.


title: Chemodynamical Modelling of the Galactic Bulge and Bar

abstract: We present the first self-consistent chemodynamical model fitted to reproduce data for the galactic bulge, bar and inner disk. We extend the Made-to-Measure method to an augmented phase-space including the metallicity of stars, and show its first application to the bar region of the Milky Way. Using data from the ARGOS and APOGEE (DR12) surveys, we adapt the recent dynamical model from Portail et al. to reproduce the observed spatial and kinematic variations as a function of metallicity, thus allowing the detailed study of the 3D density distributions, kinematics and orbital structure of stars in different metallicity bins. We find that metal-rich stars with [Fe/H] > -0.5 are strongly barred and have dynamical properties that are consistent with a common disk origin. Metal-poor stars with [Fe/H] < -0.5 show strong kinematic variations with metallicity, indicating varying contributions from the underlying stellar populations. Outside the central kpc, metal-poor stars are found to have the density and kinematics of a thick disk while in the inner kpc, evidence for an extra concentration of metal-poor stars is found. Finally, the combined orbit distributions of all metallicities in the model naturally reproduce the observed vertex deviations in the bulge. This paper demonstrates the power of Made-to-Measure chemodynamical models, that when extended to other chemical dimensions will be very powerful tools to maximize the information obtained from large spectroscopic surveys such as APOGEE, GALAH and MOONS.


title: Evolution of long-lived globular cluster stars IV. Initial helium content and white-dwarf properties

abstract: Context. Globular clusters host stars with chemical peculiarities. The associated helium enrichment is expected to affect the evolution of stars, in general, and of low-mass stars, and in particular the progenitors of white dwarfs (WDs). Aims. We investigate the effects of different initial helium contents on the properties of white dwarfs such as their masses, compositions, and the time since their formation. Results. In a globular cluster at a given age, the He enrichment implies lower initial masses for stars at a given phase. Thus it leads to a decrease of the masses of white dwarfs reaching the cooling sequence. In addition the He enrichment increases the total mass and number of white dwarfs and eventually allows the presence of He white dwarf from single progenitors. Conclusions. The low He enrichment determined in most globular clusters with different methods results in negligible effects on the white dwarf properties. However, in the few globular clusters that display a high He enrichment, this may significantly affect the characteristics of the white dwarfs. In NGC 2808 and {\omega} Centauri the high He enrichment even leads to the formation of He white dwarfs from single He-rich progenitors. Therefore investigating the white dwarf mass domain in globular clusters with a high He enrichment would provide an additional indirect way to measure and constrain the He enrichment degree.


title: Towards Understanding the Impact of Human Mobility on Police Allocation

abstract: Motivated by recent findings that human mobility is proxy for crime behavior in big cities and that there is a superlinear relationship between the people's movement and crime, this article aims to evaluate the impact of how these findings influence police allocation. More precisely, we shed light on the differences between an allocation strategy, in which the resources are distributed by clusters of floating population, and conventional allocation strategies, in which the police resources are distributed by an Administrative Area (typically based on resident population). We observed a substantial difference in the distributions of police resources allocated following these strategies, what evidences the imprecision of conventional police allocation methods.


title: Metrically Ramsey ultrafilters

abstract: Given a metric space $(X,d)$, we say that a mapping $\chi: [X]^{2}\longrightarrow\{0.1\}$ is an isometric coloring if $d(x,y)=d(z,t)$ implies $\chi(\{x,y\})=\chi(\{z,t\})$. A free ultrafilter $\mathcal{U}$ on an infinite metric space $(X,d)$ is called metrically Ramsey if, for every isometric coloring $\chi$ of $[X]^{2}$, there is a member $U\in\mathcal{U}$ such that the set $[U]^{2}$ is $\chi$-monochrome. We prove that each infinite ultrametric space $(X,d)$ has a countable subset $Y$ such that each free ultrafilter $\mathcal{U}$ on $X$ satisfying $Y\in\mathcal{U}$ is metrically Ramsey. On the other hand, it is an open question whether every metrically Ramsey ultrafilter on the natural numbers $\mathbb{N}$ with the metric $|x-y|$ is a Ramsey ultrafilter. We prove that every metrically Ramsey ultrafilter $\mathcal{U}$ on $\mathbb{N}$ has a member with no arithmetic progression of length 2, and if $\mathcal{U}$ has a thin member then there is a mapping $f:\mathbb{N}\longrightarrow\omega $ such that $f(\mathcal{U})$ is a Ramsey ultrafilter.


title: Dynamical system analysis of interacting models

abstract: We perform a dynamical system analysis of a cosmological model with linear dependence between the vacuum density and the Hubble parameter, with constant-rate creation of dark matter. We show that the de Sitter spacetime is an asymptotically stable critical point, future limit of any expanding solution. Our analysis also shows that the Minkowski spacetime is an unstable critical point, which eventually collapses to a singularity. In this way, such a prescription for the vacuum decay not only predicts the correct future de Sitter limit, but also forbids the existence of a stable Minkowski universe. We also study the effect of matter creation on the growth of structures and their peculiar velocities, showing that it is inside the current errors of redshift space distortions observations.


title: Predicting Financial Crime: Augmenting the Predictive Policing Arsenal

abstract: Financial crime is a rampant but hidden threat. In spite of this, predictive policing systems disproportionately target "street crime" rather than white collar crime. This paper presents the White Collar Crime Early Warning System (WCCEWS), a white collar crime predictive model that uses random forest classifiers to identify high risk zones for incidents of financial crime.


title: Accurate Location and Manipulation of Nano-Scaled Objects Buried under Spin-Coated Films

abstract: Detection and precise localization of nano-scale structures buried beneath spin coated films are highly valuable additions to nano-fabrication technology. In principle, the topography of the final film contains information about the location of the buried features. However, it is generally believed that the relation is masked by flow effects, which lead to an upstream shift of the dry film's topography and render precise localization impossible. Here we demonstrate, theoretically and experimentally, that the flow-shift paradigm does not apply at the sub-micron scale. Specifically, we show that the resist topography is accurately obtained from a convolution operation with a symmetric Gaussian Kernel whose parameters solely depend on the resist characteristics. We exploit this finding for a 3 nm precise overlay fabrication of metal contacts to an InAs nanowire with a diameter of 27 nm using thermal scanning probe lithography.


title: Friendships, Rivalries, and Trysts: Characterizing Relations between Ideas in Texts

abstract: Understanding how ideas relate to each other is a fundamental question in many domains, ranging from intellectual history to public communication. Because ideas are naturally embedded in texts, we propose the first framework to systematically characterize the relations between ideas based on their occurrence in a corpus of documents, independent of how these ideas are represented. Combining two statistics --- cooccurrence within documents and prevalence correlation over time --- our approach reveals a number of different ways in which ideas can cooperate and compete. For instance, two ideas can closely track each other's prevalence over time, and yet rarely cooccur, almost like a "cold war" scenario. We observe that pairwise cooccurrence and prevalence correlation exhibit different distributions. We further demonstrate that our approach is able to uncover intriguing relations between ideas through in-depth case studies on news articles and research papers.


title: Wide Binaries in Tycho-{\it Gaia}: Search Method and the Distribution of Orbital Separations

abstract: We mine the Tycho-{\it Gaia} astrometric solution (TGAS) catalog for wide stellar binaries by matching positions, proper motions, and astrometric parallaxes. We separate genuine binaries from unassociated stellar pairs through a Bayesian formulation that includes correlated uncertainties in the proper motions and parallaxes. Rather than relying on assumptions about the structure of the Galaxy, we calculate Bayesian priors and likelihoods based on the nature of Keplerian orbits and the TGAS catalog itself. We calibrate our method using radial velocity measurements and obtain 6196 high-confidence candidate wide binaries with projected separations $s\lesssim1$ pc. The normalization of this distribution suggests that at least 0.6\% of TGAS stars have an associated, distant TGAS companion in a wide binary. We demonstrate that {\it Gaia}'s astrometry is precise enough that it can detect projected orbital velocities in wide binaries with orbital periods as large as 10$^6$ yr. For pairs with $s\ \lesssim\ 4\times10^4$~AU, characterization of random alignments indicate our contamination to be $\approx$5\%. For $s \lesssim 5\times10^3$~AU, our distribution is consistent with \"{O}pik's Law. At larger separations, the distribution is steeper and consistent with a power-law $P(s)\propto s^{-1.6}$; there is no evidence in our data of any bimodality in this distribution for $s \lesssim$ 1 pc. Using radial velocities, we demonstrate that at large separations, i.e., of order $s \sim$ 1 pc and beyond, any potential sample of genuine wide binaries in TGAS cannot be easily distinguished from ionized former wide binaries, moving groups, or contamination from randomly aligned stars.


title: zBEAMS: A unified solution for supernova cosmology with redshift uncertainties

abstract: Supernova cosmology without spectra will be an important component of future surveys such as LSST. This lack of supernova spectra results in uncertainty in the redshifts which, if ignored, leads to significantly biased estimates of cosmological parameters. Here we present a hierarchical Bayesian formalism -- zBEAMS -- that addresses this problem by marginalising over the unknown or uncertain supernova redshifts to produce unbiased cosmological estimates that are competitive with supernova data with spectroscopically confirmed redshifts. zBEAMS provides a unified treatment of both photometric redshifts and host galaxy misidentification (occurring due to chance galaxy alignments or faint hosts), effectively correcting the inevitable contamination in the Hubble diagram. Like its predecessor BEAMS, our formalism also takes care of non-Ia supernova contamination by marginalising over the unknown supernova type. We illustrate this technique with simulations of supernovae with photometric redshifts and host galaxy misidentification. A novel feature of the photometric redshift case is the important role played by the redshift distribution of the supernovae.


title: Disassembling the Clockwork Mechanism

abstract: The clockwork mechanism is a means of naturally generating exponential hierarchies in theories without significant hierarchies among fundamental parameters. We emphasize the role of interactions in the clockwork mechanism, demonstrating that clockwork is an intrinsically abelian phenomenon precluded in non-abelian theories such as Yang-Mills, non-linear sigma models, and gravity. We also show that clockwork is not realized in extra-dimensional theories through purely geometric effects, but may be generated by appropriate localization of zero modes.


title: Mapping Exoplanets

abstract: The varied surfaces and atmospheres of planets make them interesting places to live, explore, and study from afar. Unfortunately, the great distance to exoplanets makes it impossible to resolve their disk with current or near-term technology. It is still possible, however, to deduce spatial inhomogeneities in exoplanets provided that different regions are visible at different times -- this can be due to rotation, orbital motion, and occultations by a star, planet, or moon. Astronomers have so far constructed maps of thermal emission and albedo for short period giant exoplanets. These maps constrain atmospheric dynamics and cloud patterns in exotic atmospheres. In the future, exo-cartography could yield surface maps of terrestrial planets, hinting at the geophysical and geochemical processes that shape them.


title: The Gaia-ESO Survey: Dynamical models of flattened, rotating globular clusters

abstract: We present a family of self-consistent axisymmetric rotating globular cluster models which are fitted to spectroscopic data for NGC 362, NGC 1851, NGC 2808, NGC 4372, NGC 5927 and NGC 6752 to provide constraints on their physical and kinematic properties, including their rotation signals. They are constructed by flattening Modified Plummer profiles, which have the same asymptotic behaviour as classical Plummer models, but can provide better fits to young clusters due to a slower turnover in the density profile. The models are in dynamical equilibrium as they depend solely on the action variables. We employ a fully Bayesian scheme to investigate the uncertainty in our model parameters (including mass-to-light ratios and inclination angles) and evaluate the Bayesian evidence ratio for rotating to non-rotating models. We find convincing levels of rotation only in NGC 2808. In the other clusters, there is just a hint of rotation (in particular, NGC 4372 and NGC 5927), as the data quality does not allow us to draw strong conclusions. Where rotation is present, we find that it is confined to the central regions, within radii of $R \leq 2 r_h$. As part of this work, we have developed a novel q-Gaussian basis expansion of the line-of-sight velocity distributions, from which general models can be constructed via interpolation on the basis coefficients.


title: Double and cyclic $\lambda$-deformations and their canonical equivalents

abstract: We prove that the doubly lambda-deformed sigma-models, which include integrable cases, are canonically equivalent to the sum of two single lambda-deformed models. This explains the equality of the exact beta-functions and current anomalous dimensions of the doubly lambda-deformed sigma-models to those of two single lambda-deformed models. Our proof is based upon agreement of their Hamiltonian densities and of their canonical structure. Subsequently, we show that it is possible to take a well defined non-Abelian type limit of the doubly-deformed action. Last, but not least, by extending the above, we construct multi-matrix integrable deformations of an arbitrary number of WZW models.


title: Nuclear effects are relevant to the calorimetric reconstruction of neutrino energy

abstract: As the calorimetric method of neutrino-energy reconstruction is generally considered to be largely insensitive to nuclear effects, its application seems to be an effective way for reducing systematic uncertainties in oscillation experiments. To verify the validity of this opinion, we quantitatively study the sensitivity of the calorimetric energy reconstruction to the effect of final-state interactions in an ideal detector and in a realistic scenario. We find that when particles escaping detection carry away a non-negligible fraction of neutrino energy, the calorimetric reconstruction method becomes sensitive to nuclear effects which, in turn, affects the outcome of the oscillation analysis. These findings suggest that the best strategy for reduction of systematic uncertainties in future neutrino-oscillation studies---such as the Deep Underground Neutrino Experiment---is to increase their sensitivity to particles of low energy. The ambitious precision goals appear also to require an extensive development of theoretical models capable of providing an accurate predictions for exclusive cross sections of well-controlled uncertainties.


title: Avoiding resonance capture in multi-planet extrasolar systems

abstract: A commonly noted feature of the population of multi-planet extrasolar systems is the rarity of planet pairs in low-order mean-motion resonances. We revisit the physics of resonance capture via convergent disk-driven migration. We point out that for planet spacings typical of stable configurations for Kepler systems, the planets can routinely maintain a small but nonzero eccentricity due to gravitational perturbations from their neighbors. Together with the upper limit on the migration rate needed for capture, the finite eccentricity can make resonance capture difficult or impossible in Sun-like systems for planets smaller than ~Neptune-sized. This mass limit on efficient capture is broadly consistent with observed exoplanet pairs that have mass determinations: of pairs with the heavier planet exterior to the lighter planet -- which would have been undergoing convergent migration in their disks -- those in or nearly in resonance are much more likely to have total mass greater than two Neptune masses than to have smaller masses. The agreement suggests that the observed paucity of resonant pairs around sun-like stars may simply arise from a small resonance capture probability for lower-mass planets. Planet pairs that thereby avoid resonance capture are much less likely to collide in an eventual close approach than to simply migrate past one another to become a divergently migrating pair with the lighter planet exterior. For systems around M stars we expect resonant pairs to be much more common, since there the minimum mass threshhold for efficient capture is about an Earth mass.


title: Biases from neutrino bias: to worry or not to worry?

abstract: The relation between the halo field and the matter fluctuations (halo bias), in the presence of massive neutrinos depends on the total neutrino mass, massive neutrinos introduce an additional scale-dependence of the bias which is usually neglected in cosmological analyses. We investigate the magnitude of the systematic effect on interesting cosmological parameters induced by neglecting this scale dependence, finding that while it is not a problem for current surveys, it is non-negligible for future, denser or deeper ones depending on the neutrino mass, the maximum scale used for the analyses and the details of the nuisance parameters considered. However there is a simple recipe to account for the bulk of the effect as to make it fully negligible, which we illustrate and advocate should be included in analysis of forthcoming large-scale structure surveys.


title: "Non-cold" dark matter at small scales: a general approach

abstract: Structure formation at small cosmological scales provides an important frontier for dark matter (DM) research. Scenarios with small DM particle masses, large momenta or hidden interactions tend to suppress the gravitational clustering at small scales. The details of this suppression depend on the DM particle nature, allowing for a direct link between DM models and astrophysical observations. However, most of the astrophysical constraints obtained so far refer to a very specific shape of the power suppression, corresponding to thermal warm dark matter (WDM), i.e., candidates with a Fermi-Dirac or Bose-Einstein momentum distribution. In this work we introduce a new analytical fitting formula for the power spectrum, which is simple yet flexible enough to reproduce the clustering signal of large classes of non-thermal DM models, which are not at all adequately described by the oversimplified notion of WDM. We show that the formula is able to fully cover the parameter space of sterile neutrinos (whether resonantly produced or from particle decay), mixed cold and warm models, fuzzy dark matter, as well as other models suggested by effective theory of structure formation (ETHOS). Based on this fitting formula, we perform a large suite of N-body simulations and we extract important nonlinear statistics, such as the matter power spectrum and the halo mass function. Finally, we present first preliminary astrophysical constraints, based on linear theory, from both the number of Milky Way satellites and the Lyman-$\alpha$ forest. This paper is a first step towards a general and comprehensive modeling of small-scale departures from the standard cold DM model.


title: Lyman-alpha Absorbers and the Coma Cluster

abstract: The spatial and kinematic distribution of warm gas in and around the Coma Cluster is presented through observations of Lyman-alpha absorbers using background QSOs. Updates to the Lyman-alpha absorber distribution found in Yoon et al. (2012) for the Virgo Cluster are also presented. At 0.2-2.0 R_vir of Coma we identify 14 Lyman-alpha absorbers (N_HI = 10^{12.8-15.9} cm^-2) towards 5 sightlines and no Lyman-alpha absorbers along 3 sightlines within 3\sigmav_coma. For both Coma and Virgo, most Lyman-alpha absorbers are found outside the virial radius or beyond 1\sigmav consistent with them largely representing the infalling intergalactic medium. The few exceptions in the central regions can be associated with galaxies. The Lyman-alpha absorbers avoid the hot ICM, consistent with the infalling gas being shock-heated within the cluster. The massive dark matter halos of clusters do not show the increasing column density with decreasing impact parameter relationship found for the smaller mass galaxy halos. In addition, while the covering fraction within R_vir is lower for clusters than galaxies, beyond R_vir the covering fraction is somewhat higher for clusters. The velocity dispersion of the absorbers compared to the galaxies is higher for Coma, consistent with the absorbers tracing additional turbulent gas motions in the cluster outskirts. The results are overall consistent with cosmological simulations, with the covering fraction being high in the observations standing out as the primary discrepancy.


title: Plethora of relativistic charged spheres: The full spectrum of Guilfoyle's static, electrically charged spherical solutions

abstract: We show that Guilfoyle's exact solutions of the Einstein-Maxwell equations for spherical symmetric static electrically charged matter with a Reissner-Nordstr\"om exterior possess a bewildering plethora of different types of solutions. For the parameter space of the solutions we use two normalized variables, $q^2/R^2$ and $r_0/R$, where $q$ is the total electric charge, $r_0$ is the radius of the object, and $R$ is a length representing the square root of the inverse energy density of the matter. The two other parameters, the mass $m$ and the Guilfoyle parameter $a$, both dependent on $q$, $r_0$ and $R$, are analyzed in detail. The full parameter space of solutions $q^2/R^2\times r_0/R$ is explored with the corresponding types of solutions being identified and analyzed. The different types of solutions are regular charged stars, including charged dust stars and stars saturating the Buchdahl-Andr\'easson bound, quasiblack holes, regular charged black holes with a de Sitter core, regular black holes with a core of phantom charged matter, other exotic regular black holes, Schwarzschild stars, Schwarzschild black holes, Kasner spacetimes, pointlike and planar naked singularities, and the Minkowski spacetime. Allowing for $q^2<0$, in which case it is not possible to interpret $q$ as electric charge, also yields new solutions, some of which are interesting and regular, others are singular. Some of these types of solutions as well as the matter properties have been previously found and studied, here the full spectrum being presented in a unified manner.


title: Transverse fields to tune an Ising-nematic quantum critical transition

abstract: The paradigmatic example of a continuous quantum phase transition is the transverse field Ising ferromagnet. In contrast to classical critical systems, whose properties depend only on symmetry and the dimension of space, the nature of a quantum phase transition also depends on the dynamics. In the transverse field Ising model, the order parameter is not conserved and increasing the transverse field enhances quantum fluctuations until they become strong enough to restore the symmetry of the ground state. Ising pseudo-spins can represent the order parameter of any system with a two-fold degenerate broken-symmetry phase, including electronic nematic order associated with spontaneous point-group symmetry breaking. Here, we show for the representative example of orbital-nematic ordering of a non-Kramers doublet that an orthogonal strain or a perpendicular magnetic field plays the role of the transverse field, thereby providing a practical route for tuning appropriate materials to a quantum critical point. While the transverse fields are conjugate to seemingly unrelated order parameters, their non-trivial commutation relations with the nematic order parameter, which can be represented by a Berry-phase term in an effective field theory, intrinsically intertwines the different order parameters.


title: Probing the spacetime evolution of heavy ion collisions with photons and dileptons

abstract: Photons and dileptons are emitted throughout the evolution of the deconfined nuclear medium produced in heavy ion collisions. As such they can provide valuable information about the different phases of the medium, and complement hadronic measurements and other observables. In this work, recent developments related to electromagnetic emissions at early time, in the cross-over region, and at late times are reviewed. The spacetime description of the nuclear medium from a hydrodynamic model of heavy ion collisions is used to provide context and guide the discussion.


title: The Velocity Dispersion Function for Quiescent Galaxies in the Local Universe

abstract: We investigate the distribution of central velocity dispersions for quiescent galaxies in the SDSS at $0.03 \leq z \leq 0.10$. To construct the field velocity dispersion function (VDF), we construct a velocity dispersion complete sample of quiescent galaxies with Dn4000$ > 1.5$. The sample consists of galaxies with central velocity dispersion larger than the velocity dispersion completeness limit of the SDSS survey. Our VDF measurement is consistent with previous field VDFs for $\sigma > 200$ km s$^{-1}$. In contrast with previous results, the VDF does not decline significantly for $\sigma < 200$ km s$^{-1}$. The field and the similarly constructed cluster VDFs are remarkably flat at low velocity dispersion ($\sigma < 250$ km s$^{-1}$). The cluster VDF exceeds the field for $\sigma > 250$ km s$^{-1}$ providing a measure of the relatively larger number of massive subhalos in clusters. The VDF is a probe of the dark matter halo distribution because the measured central velocity dispersion may be directly proportional to the dark matter velocity dispersion. Thus the VDF provides a potentially powerful test of simulations for models of structure formation.


title: Influence of the bound states in the Neumann Laplacian in a thin waveguide

abstract: We study the Neumann Laplacian operator $-\Delta_\Omega^N$ restricted to a twisted waveguide $\Omega$. The goal is to find the effective operator when the diameter of $\Omega$ tends to zero. However, when $\Omega$ is "squeezed" there are divergent eigenvalues due to the transverse oscillations. We show that each one of these eigenvalues influences the action of the effective operator in a different way. In the case where $\Omega$ is periodic and sufficiently thin, we find information about the absolutely continuous spectrum of $-\Delta_\Omega^N$ and the existence and location of band gaps in its structure.


title: Minimal Fundamental Partial Compositeness

abstract: Building upon the fundamental partial compositeness framework we provide consistent and complete composite extensions of the standard model. These are used to determine the effective operators emerging at the electroweak scale in terms of the standard model fields upon consistently integrating out the heavy composite dynamics. We exhibit the first effective field theories matching these complete composite theories of flavour and analyse their physical consequences for the third generation quarks. Relations with other approaches, ranging from effective analyses for partial compositeness to extra dimensions as well as purely fermionic extensions, are briefly discussed. Our methodology is applicable to any composite theory of dynamical electroweak symmetry breaking featuring a complete theory of flavour.


title: High energy behaviour of form factors

abstract: We solve renormalization group equations that govern infrared divergences of massless and massive form factors. By comparing to recent results for planar massive three-loop and massless four-loop form factors in QCD, we give predictions for the high-energy limit of massive form factors at the four- and for the massless form factor at five-loop order. Furthermore, we discuss the relation which connects infrared divergences regularized dimensionally and via a small quark mass and extend results present in the literature to higher order.


title: A first constraint on the average mass of ultra diffuse galaxies from weak gravitational lensing

abstract: The recent discovery of thousands of ultra diffuse galaxies (UDGs) in nearby galaxy clusters has opened a new window into the process of galaxy formation and evolution. Several scenarios have been proposed to explain the formation history of UDGs, and their ability to survive in the harsh cluster environments. A key requirement to distinguish between these scenarios is a measurement of their halo masses which, due to their low surface brightnesses, has proven difficult if one relies on stellar tracers of the potential. We exploit weak gravitational lensing, a technique that does not depend on these baryonic tracers, to measure the average subhalo mass of 784 UDGs selected in 18 clusters at $z\leq0.09$. Our sample of UDGs has a median stellar mass $\langle m_\star\rangle=2\times10^8\,\mathrm{M}_\odot$ and a median effective radius $\langle r_\mathrm{eff}\rangle=2.8$ kpc. We constrain the average mass of subhaloes within 30 kpc to $\log m_\mathrm{UDG}(r<30\,\mathrm{kpc})/\mathrm{M}_\odot\leq10.99$ at 95 per cent credibility, implying an effective virial mass $\log m_{200}/\mathrm{M}_\odot\leq11.80$, and a lower limit on the stellar mass fraction within 10 kpc of 1.0 per cent. Such mass is consistent with a simple extrapolation of the subhalo-to-stellar mass relation of typical satellite galaxies in massive clusters. However, our analysis is not sensitive to scatter about this mean mass; the possibility remains that extreme UDGs reside in haloes as massive as the Milky Way.


title: Spatial disease mapping using Directed Acyclic Graph Auto-Regressive (DAGAR) models

abstract: Hierarchical models for regionally aggregated disease incidence data commonly involve region specific latent random effects that are modeled jointly as having a multivariate Gaussian distribution. The covariance or precision matrix incorporates the spatial dependence between the regions. Common choices for the precision matrix include the widely used ICAR model, which is singular, and its nonsingular extension which lacks interpretability. We propose a new parametric model for the precision matrix based on a directed acyclic graph (DAG) representation of the spatial dependence. Our model guarantees positive definiteness and, hence, in addition to being a valid prior for regional spatially correlated random effects, can also directly model the outcome from dependent data like images and networks. Theoretical results establish a link between the parameters in our model and the variance and covariances of the random effects. Substantive simulation studies demonstrate that the improved interpretability of our model reaps benefits in terms of accurately recovering the latent spatial random effects as well as for inference on the spatial covariance parameters. Under modest spatial correlation, our model far outperforms the CAR models, while the performances are similar when the spatial correlation is strong. We also assess sensitivity to the choice of the ordering in the DAG construction using theoretical and empirical results which testify to the robustness of our model. We also present a large-scale public health application demonstrating the competitive performance of the model.


title: A structure theorem for product sets in extra special groups

abstract: Hegyv\'ari and Hennecart showed that if $B$ is a sufficiently large brick of a Heisenberg group, then the product set $B\cdot B$ contains many cosets of the center of the group. We give a new, robust proof of this theorem that extends to all extra special groups as well as to a large family of quasigroups.


title: Discovery prospects of a light Higgs boson at the LHC in type-I 2HDM

abstract: We present a comprehensive analysis of observing a light Higgs boson in the mass range $70$ -- $110$ GeV at the 13/14 TeV LHC, in the context of the type-I two-Higgs-doublet model. The decay of the light Higgs to a pair of bottom quarks is dominant in most parts of the parameter space, except in the fermiophobic limit. Here its decay to bosons, (mainly a pair of photons), becomes important. We perform an extensive collider analysis for the $b\bar{b}$ and $\gamma \gamma$ final states. The light scalar is tagged in the highly boosted regimes for the $b \bar{b}$ mode to reduce the enormous QCD background. This decay can be observed with a few thousand fb$^{-1}$ of integrated luminosity at the LHC. Near the fermiophobic limit, the decay of the light Higgs to a pair of photons can even be probed with a few hundred fb$^{-1}$ of integrated luminosity at the LHC.


title: Standard Model EFT and Extended Scalar Sectors

abstract: One of the simplest extensions of the Standard Model is the inclusion of an additional scalar multiplet, and we consider scalars in the $SU(2)_L$ singlet, triplet, and quartet representations. We examine models with heavy neutral scalars, $m_H\sim 1-2$ TeV, and the matching of the UV complete theories to the low energy effective field theory. We demonstrate the agreement of the kinematic distributions obtained in the singlet models for the gluon fusion of a Higgs pair with the predictions of the effective field theory. The restrictions on the extended scalar sectors due to unitarity and precision electroweak measurements are summarized and lead to highly restricted regions of viable parameter space for the triplet and quartet models.


title: Sub-string/Pattern Matching in Sub-linear Time Using a Sparse Fourier Transform Approach

abstract: We consider the problem of querying a string (or, a database) of length $N$ bits to determine all the locations where a substring (query) of length $M$ appears either exactly or is within a Hamming distance of $K$ from the query. We assume that sketches of the original signal can be computed off line and stored. Using the sparse Fourier transform computation based approach introduced by Pawar and Ramchandran, we show that all such matches can be determined with high probability in sub-linear time. Specifically, if the query length $M = O(N^\mu)$ and the number of matches $L=O(N^\lambda)$, we show that for $\lambda < 1-\mu$ all the matching positions can be determined with a probability that approaches 1 as $N \rightarrow \infty$ for $K \leq \frac{1}{6}M$. More importantly our scheme has a worst-case computational complexity that is only $O\left(\max\{N^{1-\mu}\log^2 N, N^{\mu+\lambda}\log N \}\right)$, which means we can recover all the matching positions in {\it sub-linear} time for $\lambda<1-\mu$. This is a substantial improvement over the best known computational complexity of $O\left(N^{1-0.359 \mu} \right)$ for recovering one matching position by Andoni {\em et al.} \cite{andoni2013shift}. Further, the number of Fourier transform coefficients that need to be computed, stored and accessed, i.e., the sketching complexity of this algorithm is only $O\left(N^{1-\mu}\log N\right)$. Several extensions of the main theme are also discussed.


title: Undecidability of the first order theories of free non-commutative Lie algebras

abstract: Let $R$ be a commutative integral unital domain and $L$ a free non-commutative Lie algebra over $R$. In this paper we show that the ring $R$ and its action on $L$ are 0-interpretable in $L$, viewed as a ring with the standard ring language $+, \cdot,0$. Furthermore, if $R$ has characteristic zero then we prove that the elementary theory $Th(L)$ of $L$ in the standard ring language is undecidable. To do so we show that the arithmetic ${\bf N} = \langle{\bf N}, +,\cdot,0 \rangle$ is 0-interpretable in $L$. This implies that the theory of $Th(L)$ has the independence property. These results answer some old questions on model theory of free Lie algebras.


title: Generating Liquid Simulations with Deformation-aware Neural Networks

abstract: We propose a novel approach for deformation-aware neural networks that learn the weighting and synthesis of dense volumetric deformation fields. Our method specifically targets the space-time representation of physical surfaces from liquid simulations. Liquids exhibit highly complex, non-linear behavior under changing simulation conditions such as different initial conditions. Our algorithm captures these complex phenomena in two stages: a first neural network computes a weighting function for a set of pre-computed deformations, while a second network directly generates a deformation field for refining the surface. Key for successful training runs in this setting is a suitable loss function that encodes the effect of the deformations, and a robust calculation of the corresponding gradients. To demonstrate the effectiveness of our approach, we showcase our method with several complex examples of flowing liquids with topology changes. Our representation makes it possible to rapidly generate the desired implicit surfaces. We have implemented a mobile application to demonstrate that real-time interactions with complex liquid effects are possible with our approach.


title: On faithfulness of the lifting for Hopf algebras and fusion categories

abstract: We use a version of Haboush's theorem over complete local Noetherian rings to prove faithfulness of the lifting for semisimple cosemisimple Hopf algebras and separable (braided, symmetric) fusion categories from characteristic $p$ to characteristic zero (arXiv/math:0203060, Section 9), showing that, moreover, any isomorphism between such structures can be reduced modulo $p$. This fills a gap in arXiv/math:0203060, Subsection 9.3. We also show that lifting of semisimple cosemisimple Hopf algebras is a fully faithful functor, and prove that lifting induces an isomorphism on Picard and Brauer-Picard groups. Finally, we show that a subcategory or quotient category of a separable multifusion category is separable (resolving an open question from arXiv/math:0203060, Subsection 9.4), and use this to show that certain classes of tensor functors between lifts of separable categories to characteristic zero can be reduced modulo $p$.


title: Separability by Piecewise Testable Languages is PTime-Complete

abstract: Piecewise testable languages form the first level of the Straubing-Th\'erien hierarchy. The membership problem for this level is decidable and testing if the language of a DFA is piecewise testable is NL-complete. The question has not yet been addressed for NFAs. We fill in this gap by showing that it is PSpace-complete. The main result is then the lower-bound complexity of separability of regular languages by piecewise testable languages. Two regular languages are separable by a piecewise testable language if the piecewise testable language includes one of them and is disjoint from the other. For languages represented by NFAs, separability by piecewise testable languages is known to be decidable in PTime. We show that it is PTime-hard and that it remains PTime-hard even for minimal DFAs.


title: Grouping-Based Random Access Collision Control for Massive Machine-Type Communication

abstract: Massive Machine-Type Communication (mMTC) is expected to be strongly supported by future 5G wireless networks. Its deployment, however, is seriously challenged by the high risk of random access (RA) collision. A popular concept is to reduce RA collisions by clustering mMTC devices into groups, and to connect group members with device-to-device (D2D) links. However, analytical models of this method and discussions about the reliability of D2D links are still absent. In this paper, existing grouping-based solutions are reviewed, an analytical model of grouped RA collision is proposed. Based on the analytical model, the impact of D2D reliability on the collision rate is also investigated. Afterwards, an efficient grouped RA procedure is designed to extend the state-of-the-art with an efficient local group update mechanism against D2D link exceptions.


title: Equation of state of a quark-meson mixture in the improved-PNJL model at finite chemical potential

abstract: We study the equation of state of QCD using an improved version of the three-flavor Polyakov-Nambu-Jona-Lasinio model beyond the mean-field approximation. It incorporates the effects of unquenched quarks into the Polyakov-loop effective potential, as well as mesonic contributions to the grand-canonical potential. We study in full detail the calculation of the thermodynamical potential in this approach and compare the resulting pressure and entropy density with the most-recent lattice-QCD calculations at zero baryochemical potential. Finally, we present some exploratory results at finite chemical potential which include the phase diagram of the model, the quark and meson masses, and finally, the thermodynamical pressure.


title: Higher-Spin Gauge Theories and Bulk Locality

abstract: We present a no-go result on consistent Noether interactions among higher-spin gauge fields on anti-de Sitter space-times. We show that there is a non-local obstruction at the classical level to consistent interacting field theory descriptions of massless higher-spin particles that are described in the free limit by the free Fronsdal action, under the assumption that such theories arise from the gauging of a global higher-spin symmetry. Our result suggests that the Fronsdal programme for introducing interactions among higher-spin gauge fields cannot be completed without introducing new guiding principles, which could lie beyond the framework of classical field theory.


title: Universality of Confluent, Self-Loop Deterministic Partially Ordered NFAs is Hard

abstract: An automaton is partially ordered if the only cycles in its transition diagram are self-loops. The expressivity of partially ordered NFAs (poNFAs) can be characterized by the Straubing-Th\'erien hierarchy. Level 3/2 is recognized by poNFAs, level 1 by confluent, self-loop deterministic poNFAs as well as by confluent poDFAs, and level 1/2 by saturated poNFAs. We study the universality problem for confluent, self-loop deterministic poNFAs. It asks whether an automaton accepts all words over its alphabet. Universality for both NFAs and poNFAs is a PSpace-complete problem. For confluent, self-loop deterministic poNFAs, the complexity drops to coNP-complete if the alphabet is fixed but is open if the alphabet may grow. We solve this problem by showing that it is PSpace-complete if the alphabet may grow polynomially. Consequently, our result provides a lower-bound complexity for some other problems, including inclusion, equivalence, and $k$-piecewise testability. Since universality for saturated poNFAs is easy, confluent, self-loop deterministic poNFAs are the simplest and natural kind of NFAs characterizing a well-known class of languages, for which deciding universality is as difficult as for general NFAs.


title: A Muon Source Proton Driver At JPARC-Based Parameters

abstract: An "ultimate" high intensity proton source for neutrino factories and/or muon colliders was projected to be a ~4 MW multi-GeV proton source providing short, intense proton pulses at ~15 Hz. The JPARC ~1 MW accelerators provide beam at parameters that in many respects overlap these goals. Proton pulses from the JPARC Main Ring can readily meet the pulsed intensity goals. We explore these parameters, describing the overlap and consider extensions that may take a JPARC-like facility toward this "ultimate" source. JPARC itself could serve as a stage 1 source for such a facility.


title: Role of multiorbital effects in the magnetic phase diagram of iron-pnictides

abstract: We elucidate the pivotal role of the bandstructure's orbital content in deciding the type of commensurate magnetic order stabilized within the itinerant scenario of iron-pnictides. Recent experimental findings in the tetragonal magnetic phase attest to the existence of the so-called charge and spin ordered density wave over the spin-vortex crystal phase, the latter of which tends to be favored in simplified band models of itinerant magnetism. Here we show that employing a multiorbital itinerant Landau approach based on realistic bandstructures can account for the experimentally observed magnetic phase, and thus shed light on the importance of the orbital content in deciding the magnetic order. In addition, we remark that the presence of a hole pocket centered at the Brillouin zone's ${\rm M}$-point favors a magnetic stripe rather than a tetragonal magnetic phase. For inferring the symmetry properties of the different magnetic phases, we formulate our theory in terms of magnetic order parameters transforming according to irreducible representations of the ensuing D$_{\rm 4h}$ point group. The latter method not only provides transparent understanding of the symmetry breaking schemes but also reveals that the leading instabilities always belong to the $\{A_{1g},B_{1g}\}$ subset of irreducible representations, independent of their C$_2$ or C$_4$ nature.


title: Multi-View Dynamic Facial Action Unit Detection

abstract: We propose a novel convolutional neural network approach to address the fine-grained recognition problem of multi-view dynamic facial action unit detection. We leverage recent gains in large-scale object recognition by formulating the task of predicting the presence or absence of a specific action unit in a still image of a human face as holistic classification. We then explore the design space of our approach by considering both shared and independent representations for separate action units, and also different CNN architectures for combining color and motion information. We then move to the novel setup of the FERA 2017 Challenge, in which we propose a multi-view extension of our approach that operates by first predicting the viewpoint from which the video was taken, and then evaluating an ensemble of action unit detectors that were trained for that specific viewpoint. Our approach is holistic, efficient, and modular, since new action units can be easily included in the overall system. Our approach significantly outperforms the baseline of the FERA 2017 Challenge, with an absolute improvement of 14% on the F1-metric. Additionally, it compares favorably against the winner of the FERA 2017 challenge. Code source is available at https://github.com/BCV-Uniandes/AUNets.


title: Radiation-driven turbulent accretion onto massive black holes

abstract: Accretion of gas and interaction of matter and radiation are at the heart of many questions pertaining to black hole (BH) growth and coevolution of massive BHs and their host galaxies. To answer them it is critical to quantify how the ionizing radiation that emanates from the innermost regions of the BH accretion flow couples to the surrounding medium and how it regulates the BH fueling. In this work we use high resolution 3-dimensional (3D) radiation-hydrodynamic simulations with the code Enzo, equipped with adaptive ray tracing module Moray, to investigate radiation-regulated BH accretion of cold gas. Our simulations reproduce findings from an earlier generation of 1D/2D simulations: the accretion powered UV and X-ray radiation forms a highly ionized bubble, which leads to suppression of BH accretion rate characterized by quasi-periodic outbursts. A new feature revealed by the 3D simulations is the highly turbulent nature of the gas flow in vicinity of the ionization front. During quiescent periods between accretion outbursts, the ionized bubble shrinks in size and the gas density that precedes the ionization front increases. Consequently, the 3D simulations show oscillations in the accretion rate of only ~2-3 orders of magnitude, significantly smaller than 1D/2D models. We calculate the energy budget of the gas flow and find that turbulence is the main contributor to the kinetic energy of the gas but corresponds to less than 10% of its thermal energy and thus does not contribute significantly to the pressure support of the gas.


title: Robust Estimators and Test-Statistics for One-Shot Device Testing Under the Exponential Distribution

abstract: This paper develops a new family of estimators, the minimum density power divergence estimators (MDPDEs), for the parameters of the one-shot device model as well as a new family of test statistics, Z-type test statistics based on MDPDEs, for testing the corresponding model parameters. The family of MDPDEs contains as a particular case the maximum likelihood estimator (MLE) considered in Balakrishnan and Ling (2012). Through a simulation study, it is shown that some MDPDEs have a better behavior than the MLE in relation to robustness. At the same time, it can be seen that some Z-type tests based on MDPDEs have a better behavior than the classical Z-test statistic also in terms of robustness.


title: Identifying spatiotemporal dynamics of Ebola in Sierra Leone using virus genomes

abstract: Containing the recent West African outbreak of Ebola virus (EBOV) required the deployment of substantial global resources. Operationally, health workers and surveillance teams treated cases, collected genetic samples, and tracked case contacts. Despite the substantial progress in analyzing and modeling EBOV epidemiological data, a complete characterization of the spatiotemporal spread of Ebola cases remains a challenge. In this work, we offer a novel perspective on the EBOV epidemic that utilizes virus genome sequences to inform population-level, spatial models. Calibrated to phylogenetic linkages, these dynamic spatial models provide unique insight into the disease mobility of EBOV in Sierra Leone. Further, we developed a model selection framework that identifies important epidemiological variables influencing the spatiotemporal propagation of EBOV. Consistent with other investigations, our results show that the spread of EBOV during the beginning and middle portions of the epidemic strongly depended on the size of and distance between populations. Our analysis also revealed a substantial decline in the dependence on population size at the end of the epidemic, coinciding with the large-scale intervention campaign: Operation Western Area Surge. More generally, we believe this framework, pairing molecular diagnostics with dynamic models, has the potential to be a powerful forecasting tool along with offering operationally-relevant guidance for surveillance and sampling strategies during an epidemic.


title: Nanoantenna-Microcavity Hybrids with Highly Cooperative Plasmonic-Photonic Coupling

abstract: Nanoantennas offer the ultimate spatial control over light by concentrating optical energy well below the diffraction limit, whereas their quality factor (Q) is constrained by large radiative and dissipative losses. Dielectric microcavities, on the other hand, are capable of generating a high Q-factor through an extended photon storage time but have a diffraction-limited optical mode volume. Here we bridge the two worlds, by studying an exemplary hybrid system integrating plasmonic gold nanorods acting as nanoantennas with an on-resonance dielectric photonic crystal (PC) slab acting as a low-loss microcavity and, more importantly, by synergistically combining their advantages to produce a much stronger local field enhancement than that of the separate entities. To achieve this synergy between the two polar opposite types of nanophotonic resonant elements, we show that it is crucial to coordinate both the dissipative loss of the nanoantenna and the Q-factor of the low-loss cavity. In comparison to the antenna-cavity coupling approach using a Fabry-Perot resonator, which has proved successful for resonant amplification of the antenna's local field intensity, we theoretically and experimentally show that coupling to a modest-Q PC guided resonance can produce a greater amplification by at least an order of magnitude. The synergistic nanoantenna-microcavity hybrid strategy opens new opportunities for further enhancing nanoscale light-matter interactions to benefit numerous areas such as nonlinear optics, nanolasers, plasmonic hot carrier technology, and surface-enhanced Raman and infrared absorption spectroscopies.


title: New robust statistical procedures for polytomous logistic regression models

abstract: This paper derives a new family of estimators, namely the minimum density power divergence estimators, as a robust generalization of the maximum likelihood estimator for the polytomous logistic regression model. Based on these estimators, a family of Wald-type test statistics for linear hypotheses is introduced. Robustness properties of both the proposed estimators and the test statistics are theoretically studied through the classical influence function analysis. Appropriate real life examples are presented to justify the requirement of suitable robust statistical procedures in place of the likelihood based inference for the polytomous logistic regression model. The validity of the theoretical results established in the paper are further confirmed empirically through suitable simulation studies. Finally, an approach for the data-driven selection of the robustness tuning parameter is proposed with empirical justifications.


title: On a free boundary problem and minimal surfaces

abstract: From minimal surfaces such as Simons' cone and catenoids, using refined Lyapunov-Schmidt reduction method, we construct new solutions for a free boundary problem whose free boundary has two components. In dimension $8$, using variational arguments, we also obtain solutions which are global minimizers of the corresponding energy functional. This shows that Savin's theorem is optimal.


title: Configurations of Linear Spaces of Codimension Two and the Containment Problem

abstract: We show that the union of certain codimension two flats in Fermat hyperplane arrangements does not have the symbolic cube of the defining ideal contained in the square of the ideal.