title: VNect: Real-time 3D Human Pose Estimation with a Single RGB Camera

abstract: We present the first real-time method to capture the full global 3D skeletal pose of a human in a stable, temporally consistent manner using a single RGB camera. Our method combines a new convolutional neural network (CNN) based pose regressor with kinematic skeleton fitting. Our novel fully-convolutional pose formulation regresses 2D and 3D joint positions jointly in real time and does not require tightly cropped input frames. A real-time kinematic skeleton fitting method uses the CNN output to yield temporally stable 3D global pose reconstructions on the basis of a coherent kinematic skeleton. This makes our approach the first monocular RGB method usable in real-time applications such as 3D character control---thus far, the only monocular methods for such applications employed specialized RGB-D cameras. Our method's accuracy is quantitatively on par with the best offline 3D monocular RGB pose estimation methods. Our results are qualitatively comparable to, and sometimes better than, results from monocular RGB-D approaches, such as the Kinect. However, we show that our approach is more broadly applicable than RGB-D solutions, i.e. it works for outdoor scenes, community videos, and low quality commodity RGB cameras.


title: $q$-Generalized representation of the $d$-dimensional Dirac delta and $q$-Fourier transform

abstract: We discuss a generalized representation of the Dirac delta function in $d$ dimensions in terms of $q$-exponential functions. We apply this new representation to the study of the so-called $q$-Fourier transform, proving its invertibility for any value of $d$. We finally illustrate the effect of the $q$-deformation on the Gibbs phenomenon.


title: Reality from maximizing overlap in the future-included real action theory

abstract: In the future-included real action theory whose path runs over not only past but also future, we demonstrate a theorem, which states that the normalized matrix element of a Hermitian operator $\hat{\cal O}$ defined in terms of the future state at the final time $T_B$ and the fixed past state at the initial time $T_A$ becomes real for the future state selected such that the absolute value of the transition amplitude from the past state to the future state is maximized. This is a special version of our previously proposed theorem for the future-included complex action theory. We find that though the maximization principle leads to the reality of the normalized matrix element in the future-included real action theory, it does not specify the future and past states so much as in the case of the future-included complex action theory. In addition, we argue that the normalized matrix element seems to be more natural than the usual expectation value. Thus we speculate that the functional integral formalism of quantum theory could be most elegant in the future-included complex action theory.


title: The power graph of a torsion-free group

abstract: The \emph{power graph} $P(G)$ of a group $G$ is the graph whose vertex set is $G$, with $x$ and $y$ joined if one is a power of the other; the \emph{directed power graph} $\vec{P}(G)$ has the same vertex set, with an arc from $x$ to $y$ if $y$ is a power of $x$. It is known that, for finite groups, the power graph determines the directed power graph up to isomorphism. However, it is not true that any isomorphism between power graphs induces an isomorphism between directed power graphs. Moreover, for infinite groups the power graph may fail to determine the directed power graph. In this paper, we consider power graphs of torsion-free groups. Our main results are that, for torsion-free nilpotent groups of class at most $2$, and for groups in which every non-identity element lies in a unique maximal cyclic subgroup, the power graph determines the directed power graph up to isomorphism. For specific groups such as $\mathbb{Z}$ and $\mathbb{Q}$, we obtain more precise results. Any isomorphism $P(\mathbb{Z})\to P(G)$ preserves orientation, so induces an isomorphism between directed power graphs; in the case of $\mathbb{Q}$, the orientations are either all preserved or all reversed. We also obtain results about groups in which every element is contained in a unique maximal cyclic subgroup (this class includes the free and free abelian groups), and about subgroups of the additive group of $\mathbb{Q}$ and about $\mathbb{Q}^n$.


title: Convergence of Positive Operator Semigroups

abstract: We present new conditions for semigroups of positive operators to converge strongly as time tends to infinity. Our proofs are based on a novel approach combining the well-known splitting theorem by Jacobs, de Leeuw and Glicksberg with a purely algebraic result about positive group representations. Thus we obtain convergence theorems not only for one-parameter semigroups but for a much larger class of semigroup representations. Our results allow for a unified treatment of various theorems from the literature that, under technical assumptions, a bounded positive $C_0$-semigroup containing or dominating a kernel operator converges strongly as $t \to \infty$. We gain new insights into the structure theoretical background of those theorems and generalise them in several respects; especially we drop any kind of continuity or regularity assumption with respect to the time parameter.


title: Legacy ExtraGalactic UV Survey with The Hubble Space Telescope. Stellar cluster catalogues and first insights into cluster formation and evolution in NGC 628

abstract: We report the large effort which is producing comprehensive high-level young star cluster (YSC) catalogues for a significant fraction of galaxies observed with the Legacy ExtraGalactic UV Survey (LEGUS) Hubble treasury program. We present the methodology developed to extract cluster positions, verify their genuine nature, produce multiband photometry (from NUV to NIR), and derive their physical properties via spectral energy distribution fitting analyses. We use the nearby spiral galaxy NGC628 as a test case for demonstrating the impact that LEGUS will have on our understanding of the formation and evolution of YSCs and compact stellar associations within their host galaxy. Our analysis of the cluster luminosity function from the UV to the NIR finds a steepening at the bright end and at all wavelengths suggesting a dearth of luminous clusters. The cluster mass function of NGC628 is consistent with a power-law distribution of slopes $\sim -2$ and a truncation of a few times $10^5$ M$_\odot$. After their formation YSCs and compact associations follow different evolutionary paths. YSCs survive for a longer timeframe, confirming their being potentially bound systems. Associations disappear on time scales comparable to hierarchically organized star-forming regions, suggesting that they are expanding systems. We find mass-independent cluster disruption in the inner region of NGC628, while in the outer part of the galaxy there is little or no disruption. We observe faster disruption rates for low mass ($\leq$ $10^4$ M$_\odot$) clusters suggesting that a mass-dependent component is necessary to fully describe the YSC disruption process in NGC628.


title: Stable Secretaries

abstract: We define and study a new variant of the secretary problem. Whereas in the classic setting multiple secretaries compete for a single position, we study the case where the secretaries arrive one at a time and are assigned, in an on-line fashion, to one of multiple positions. Secretaries are ranked according to talent, as in the original formulation, and in addition positions are ranked according to attractiveness. To evaluate an online matching mechanism, we use the notion of blocking pairs from stable matching theory: our goal is to maximize the number of positions (or secretaries) that do not take part in a blocking pair. This is compared with a stable matching in which no blocking pair exists. We consider the case where secretaries arrive randomly, as well as that of an adversarial arrival order, and provide corresponding upper and lower bounds.


title: First spectropolarimetric measurement of a brown dwarf magnetic field in molecular bands

abstract: We present the first measurement of the surface magnetic field of a late-M dwarf, LSR J1835+3259, with the help of the full-Stokes spectropolarimetry in the bands of diatomic molecules. Our magnetic field measurements at different rotational phases of the dwarf yielded one 5$\sigma$ and two 3$\sigma$ detections. The observational data have been obtained with the LRISp polarimeter at the Keck observatory on August 22nd and 23rd, 2012. These data have been compared against synthetic full-Stokes spectra in the bands of the molecules CrH, FeH and TiO, which have been calculated for a range of the stellar parameters and magnetic field strengths. Making use of $\chi^2$-minimization and maximum likelihood estimation, we determine the net magnetic field strength $B$ (and not flux $Bf$) of LSR J1835+3259 to $\sim5$ kG, mainly with the help of the Paschen-Back effect in the CrH lines. Our measurements at different rotational phases suggest that the dwarf's surface might be covered with strong small-scale magnetic fields. In addition, recent findings of the hydrogen emission and the Stokes V signal in the Na I near-infrared doublet indicate that the dwarf's surface magnetic field might be changing rapidly giving rise to flare activity, similar to young dMe dwarfs. We substantiate the substellar origin of LSR J1835+3259 by making use of our own data as well as the photometric data from the all-sky surveys 2MASS and WISE.


title: Demonstrating research subcommunities in mathematical networks

abstract: We propose a method for demonstrating sub community structure in scientific networks of relatively small size from analyzing databases of publications. Research relationships between the network members can be visualized as a graph with vertices corresponding to authors and with edges indicating joint authorship. Using a fast clustering algorithm combined with a graph layout algorithm, we demonstrate how to display these clustering results in an attractive and informative way. The small size of the graph allows us to develop tools that keep track of how these research sub communities evolve in time, as well as to present the research articles that create the links between the network members. These tools are included in a web app, where the visitor can easily identify the various sub communities, providing also valuable information for administrational purposes. Our method was developed for the GEAR mathematical network and it can be applied to other networks.


title: Short- and long-baseline sterile neutrino phenomenology

abstract: Several anomalies observed in short-baseline neutrino experiments indicate that the standard 3-flavor framework may be incomplete and point towards the existence of light sterile neutrinos. Here, we present a concise review of the status of the neutrino oscillations within the 3+1 scheme, which is a minimal extension of the standard 3-flavor framework with one sterile neutrino species. We emphasize the potential role of LBL experiments in the searches of CP violation connected to sterile neutrinos and their complementarity with the SBL experiments.


title: A Bound on the Spectral Radius of Hypergraphs with $e$ Edges

abstract: For $r\geq 3$, let $f_r\colon [0,\infty)\to [1,\infty)$ be the unique analytic function such that $f_r({k\choose r})={k-1\choose r-1}$ for any $k\geq r-1$. We prove that the spectral radius of an $r$-uniform hypergraph $H$ with $e$ edges is at most $f_r(e)$. The equality holds if and only if $e={k\choose r}$ for some positive integer $k$ and $H$ is the union of a complete $r$-uniform hypergraph $K_k^r$ and some possible isolated vertices. This result generalizes the classical Stanley's theorem on graphs.


title: Some homological properties of ideals with cohomological dimension one

abstract: Let R denote a commutative Noetherian ring and let I be an ideal of R such that H_i^I(R) = 0, for all integers i greater than or equal to 2. In this paper we shall prove some results concerning the homological properties of I.


title: Homomorphisms Are a Good Basis for Counting Small Subgraphs

abstract: We introduce graph motif parameters, a class of graph parameters that depend only on the frequencies of constant-size induced subgraphs. Classical works by Lov\'asz show that many interesting quantities have this form, including, for fixed graphs $H$, the number of $H$-copies (induced or not) in an input graph $G$, and the number of homomorphisms from $H$ to $G$. Using the framework of graph motif parameters, we obtain faster algorithms for counting subgraph copies of fixed graphs $H$ in host graphs $G$: For graphs $H$ on $k$ edges, we show how to count subgraph copies of $H$ in time $k^{O(k)}\cdot n^{0.174k + o(k)}$ by a surprisingly simple algorithm. This improves upon previously known running times, such as $O(n^{0.91k + c})$ time for $k$-edge matchings or $O(n^{0.46k + c})$ time for $k$-cycles. Furthermore, we prove a general complexity dichotomy for evaluating graph motif parameters: Given a class $\mathcal C$ of such parameters, we consider the problem of evaluating $f\in \mathcal C$ on input graphs $G$, parameterized by the number of induced subgraphs that $f$ depends upon. For every recursively enumerable class $\mathcal C$, we prove the above problem to be either FPT or #W[1]-hard, with an explicit dichotomy criterion. This allows us to recover known dichotomies for counting subgraphs, induced subgraphs, and homomorphisms in a uniform and simplified way, together with improved lower bounds. Finally, we extend graph motif parameters to colored subgraphs and prove a complexity trichotomy: For vertex-colored graphs $H$ and $G$, where $H$ is from a fixed class $\mathcal H$, we want to count color-preserving $H$-copies in $G$. We show that this problem is either polynomial-time solvable or FPT or #W[1]-hard, and that the FPT cases indeed need FPT time under reasonable assumptions.


title: Capacity of Burst Noise-Erasure Channels With and Without Feedback and Input Cost

abstract: A class of burst noise-erasure channels which incorporate both errors and erasures during transmission is studied. The channel, whose output is explicitly expressed in terms of its input and a stationary ergodic noise-erasure process, is shown to have a so-called "quasi-symmetry" property under certain invertibility conditions. As a result, it is proved that a uniformly distributed input process maximizes the channel's block mutual information, resulting in a closed-form formula for its non-feedback capacity in terms of the noise-erasure entropy rate and the entropy rate of an auxiliary erasure process. The feedback channel capacity is also characterized, showing that feedback does not increase capacity and generalizing prior related results. The capacity-cost function of the channel with and without feedback is also investigated. A sequence of finite-letter upper bounds for the capacity-cost function without feedback is derived. Finite-letter lower bonds for the capacity-cost function with feedback are obtained using a specific encoding rule. Based on these bounds, it is demonstrated both numerically and analytically that feedback can increase the capacity-cost function for a class of channels with Markov noise-erasure processes.


title: Testing different approaches to quantum gravity with cosmology: An overview

abstract: Among the available quantum gravity proposals, string theory, loop quantum gravity, non-commutative geometry, group field theory, causal sets, asymptotic safety, causal dynamical triangulation, emergent gravity are among the best motivated models. As an introductory summary to this special issue of Comptes Rendus Physique, I explain how those different theories can be tested or constrained by cosmological observations.


title: cuTT: A High-Performance Tensor Transpose Library for CUDA Compatible GPUs

abstract: We introduce the CUDA Tensor Transpose (cuTT) library that implements high-performance tensor transposes for NVIDIA GPUs with Kepler and above architectures. cuTT achieves high performance by (a) utilizing two GPU-optimized transpose algorithms that both use a shared memory buffer in order to reduce global memory access scatter, and by (b) computing memory positions of tensor elements using a thread-parallel algorithm. We evaluate the performance of cuTT on a variety of benchmarks with tensor ranks ranging from 2 to 12 and show that cuTT performance is independent of the tensor rank and that it performs no worse than an approach based on code generation. We develop a heuristic scheme for choosing the optimal parameters for tensor transpose algorithms by implementing an analytical GPU performance model that can be used at runtime without need for performance measurements or profiling. Finally, by integrating cuTT into the tensor algebra library TAL-SH, we significantly reduce the tensor transpose overhead in tensor contractions, achieving as low as just one percent overhead for arithmetically intensive tensor contractions.


title: Characterization and Photometric Performance of the Hyper Suprime-Cam Software Pipeline

abstract: The Subaru Strategic Program (SSP) is an ambitious multi-band survey using the Hyper Suprime-Cam (HSC) on the Subaru telescope. The Wide layer of the SSP is both wide and deep, reaching a detection limit of i~26.0 mag. At these depths, it is challenging to achieve accurate, unbiased, and consistent photometry across all five bands. The HSC data are reduced using a pipeline that builds on the prototype pipeline for the Large Synoptic Survey Telescope. We have developed a Python-based, flexible framework to inject synthetic galaxies into real HSC images called SynPipe. Here we explain the design and implementation of SynPipe and generate a sample of synthetic galaxies to examine the photometric performance of the HSC pipeline. For stars, we achieve 1% photometric precision at i~19.0 mag and 6% precision at i~25.0 in the i-band. For synthetic galaxies with single-Sersic profiles, forced CModel photometry achieves 13% photometric precision at i~20.0 mag and 18% precision at i~25.0 in the i-band. We show that both forced PSF and CModel photometry yield unbiased color estimates that are robust to seeing conditions. We identify several caveats that apply to the version of HSC pipeline used for the first public HSC data release (DR1) that need to be taking into consideration. First, the degree to which an object is blended with other objects impacts the overall photometric performance. This is especially true for point sources. Highly blended objects tend to have larger photometric uncertainties, systematically underestimated fluxes and slightly biased colors. Second, >20% of stars at 22.5< i < 25.0 mag can be misclassified as extended objects. Third, the current CModel algorithm tends to strongly underestimate the half-light radius and ellipticity of galaxy with i>21.5 mag.


title: Coupling Polynomial Stratonovich Integrals: the two-dimensional Brownian case

abstract: We show how to build an immersion coupling of a two-dimensional Brownian motion $(W_1, W_2)$ along with $\binom{n}{2} + n= \tfrac12n(n+1)$ integrals of the form $\int W_1^iW_2^j \circ dW_2$, where $j=1,\ldots,n$ and $i=0, \ldots, n-j$ for some fixed $n$. The resulting construction is applied to the study of couplings of certain hypoelliptic diffusions (driven by two-dimensional Brownian motion using polynomial vector fields). This work follows up previous studies concerning coupling of Brownian stochastic areas and time integrals (Ben Arous, Cranston and Kendall (1995), Kendall and Price (2004), Kendall (2007), Kendall (2009), Kendall (2013), Banerjee and Kendall (2015), Banerjee, Gordina and Mariano (2016)), and is part of an ongoing research programme aimed at gaining a better understanding of when it is possible to couple not only diffusions but also multiple selected integral functionals of the diffusions.


title: Semi-supervised cross-entropy clustering with information bottleneck constraint

abstract: In this paper, we propose a semi-supervised clustering method, CEC-IB, that models data with a set of Gaussian distributions and that retrieves clusters based on a partial labeling provided by the user (partition-level side information). By combining the ideas from cross-entropy clustering (CEC) with those from the information bottleneck method (IB), our method trades between three conflicting goals: the accuracy with which the data set is modeled, the simplicity of the model, and the consistency of the clustering with side information. Experiments demonstrate that CEC-IB has a performance comparable to Gaussian mixture models (GMM) in a classical semi-supervised scenario, but is faster, more robust to noisy labels, automatically determines the optimal number of clusters, and performs well when not all classes are present in the side information. Moreover, in contrast to other semi-supervised models, it can be successfully applied in discovering natural subgroups if the partition-level side information is derived from the top levels of a hierarchical clustering.


title: Inflation in an external four-form gauge field

abstract: We consider the possibility that the universe, viewed as a three-brane, originated in a region of strong external field strength due to a four-form gauge field in the bulk. It is shown that in a scenario of this kind inflation is generic for a wide range of initial conditions. This is true even for small field inflation with a simple quadratic inflaton potential, as well as for Higgs potentials with the initial field well away from the local maximum, not necessarily starting from rest. The power spectrum, spectral index, and $r$ parameter help to constrain parameters in this scenario, and $r<0.1$ favors Higgs potentials.


title: Vortex sheets and diffeomorphism groupoids

abstract: In 1966 V.Arnold suggested a group-theoretic approach to ideal hydrodynamics in which the motion of an inviscid incompressible fluid is described as the geodesic flow of the right-invariant $L^2$-metric on the group of volume-preserving diffeomorphisms of the flow domain. Here we propose geodesic, group-theoretic, and Hamiltonian frameworks to include fluid flows with vortex sheets. It turns out that the corresponding dynamics is related to a certain groupoid of pairs of volume-preserving diffeomorphisms with common interface. We also develop a general framework for Euler-Arnold equations for geodesics on groupoids equipped with one-sided invariant metrics.


title: Emerging dynamics arising from coarse-grained quantum systems

abstract: The purpose of physics is to describe nature from elementary particles all the way up to cosmological objects like cluster of galaxies and black holes. Although a unified description for all this spectrum of events is desirable, this would be highly impractical. To not get lost in unnecessary details, effective descriptions are mandatory. Here we analyze the dynamics that may emerge from a full quantum description when one does not have access to all the degrees of freedom of a system. More concretely, we describe the properties of the dynamics that arise from quantum mechanics if one has access only to a coarse-grained description of the system. We obtain that the effective maps are not necessarily of Kraus form, due to correlations between accessible and nonaccessible degrees of freedom, and that the distance between two effective states may increase under the action of the effective map. We expect our framework to be useful for addressing questions such as the thermalization of closed quantum systems, as well as the description of measurements in quantum mechanics.


title: Implications of the complex Singlet field for Noncommutative Geometry model

abstract: We consider a complex singlet scalar in the spectral action approach to the standard model. It is shown that there is a range of initial values at the unification scale which is able to produce Higgs and top quark masses at low energies. The stability of the vacuum and the deviation of gauge couplings from experimental values are discussed and compared at the two-loop level with a real scalar singlet and the pure standard model.


title: Gemini Observations of Galaxies in Rich Early Environments (GOGREEN) I: Survey Description

abstract: We describe a new Large Program in progress on the Gemini North and South telescopes: Gemini Observations of Galaxies in Rich Early Environments (GOGREEN). This is an imaging and deep spectroscopic survey of 21 galaxy systems at $1<z<1.5$, selected to span a factor $>10$ in halo mass. The scientific objectives include measuring the role of environment in the evolution of low-mass galaxies, and measuring the dynamics and stellar contents of their host haloes. The targets are selected from the SpARCS, SPT, COSMOS and SXDS surveys, to be the evolutionary counterparts of today's clusters and groups. The new red-sensitive Hamamatsu detectors on GMOS, coupled with the nod-and-shuffle sky subtraction, allow simultaneous wavelength coverage over $\lambda\sim 0.6$--$1.05\mu$m, and this enables a homogeneous and statistically complete redshift survey of galaxies of all types. The spectroscopic sample targets galaxies with AB magnitudes $z^{\prime}<24.25$ and [3.6]$\mu$m$<22.5$, and is therefore statistically complete for stellar masses $M_\ast\gtrsim10^{10.3}M_\odot$, for all galaxy types and over the entire redshift range. Deep, multiwavelength imaging has been acquired over larger fields for most systems, spanning $u$ through $K$, in addition to deep IRAC imaging at 3.6$\mu$m. The spectroscopy is $\sim 50$ per cent complete as of semester 17A, and we anticipate a final sample of $\sim 500$ new cluster members. Combined with existing spectroscopy on the brighter galaxies from GCLASS, SPT and other sources, GOGREEN will be a large legacy cluster and field galaxy sample at this redshift that spectroscopically covers a wide range in stellar mass, halo mass, and clustercentric radius.


title: Edge Shear Flows and Particle Transport near the Density Limit in the HL-2A Tokamak

abstract: Edge shear flow and its effect on regulating turbulent transport have long been suspected to play an important role in plasmas operating near the Greenwald density limit $ n_G $. In this study, equilibrium profiles as well as the turbulent particle flux and Reynolds stress across the separatrix in the HL-2A tokamak are examined as $ n_G $ is approached in ohmic L-mode discharges. As the normalized line-averaged density $ \bar{n}_e/n_G $ is raised, the shearing rate of the mean poloidal flow $ \omega_{\rm sh} $ drops, and the turbulent drive for the low-frequency zonal flow (the Reynolds power $ \mathcal{P}_{Re} $) collapses. Correspondingly, the turbulent particle transport increases drastically with increasing collision rates. The geodesic acoustic modes (GAMs) gain more energy from the ambient turbulence at higher densities, but have smaller shearing rate than low-frequency zonal flows. The increased density also introduces decreased adiabaticity which not only enhances the particle transport but is also related to a reduction in the eddy-tilting and the Reynolds power. Both effects may lead to the cooling of edge plasmas and therefore the onset of MHD instabilities that limit the plasma density.


title: Towards the solution of the many-electron problem in real materials: equation of state of the hydrogen chain with state-of-the-art many-body methods

abstract: We present numerical results for the equation of state of an infinite chain of hydrogen atoms. A variety of modern many-body methods are employed, with exhaustive cross-checks and validation. Approaches for reaching the continuous space limit and the thermodynamic limit are investigated, proposed, and tested. The detailed comparisons provide a benchmark for assessing the current state of the art in many-body computation, and for the development of new methods. The ground-state energy per atom in the linear chain is accurately determined versus bondlength, with a confidence bound given on all uncertainties.


title: Essential Dimension of Generic Symbols in Characteristic p

abstract: In this article the $p$-essential dimension of generic symbols over fields of characteristic $p$ is studied. In particular, the $p$-essential dimension of the length $\ell$ generic $p$-symbol of degree $n+1$ is bounded below by $n+\ell$ when the base field is algebraically closed of characteristic $p$. The proof uses new techniques for working with residues in Milne-Kato $p$-cohomology and builds on work of Babic and Chernousov in the Witt group in characteristic 2. Two corollaries on $p$-symbol algebras (i.e, degree 2 symbols) result from this work. The generic $p$-symbol algebra of length $\ell$ is shown to have $p$-essential dimension equal to $\ell+1$ as a $p$-torsion Brauer class. The second is a lower bound of $\ell+1$ on the $p$-essential dimension of the functor $\mathrm{Alg}_{p^\ell,p}$. Roughly speaking this says that you will need at least $\ell+1$ independent parameters to be able to specify any given algebra of degree $p^{\ell}$ and exponent $p$ over a field of characteristic $p$ and improves on the previously established lower bound of 3.


title: Memory Effects in Active Particles with Exponentially Correlated Propulsion

abstract: The Ornstein--Uhlenbeck Particle (OUP) model imagines a microscopic swimmer propelled by an active force which is correlated with itself on a finite time-scale. Here we investigate the influence of external potentials on an ideal suspension of OUPs, in both one and two spatial dimensions, with particular attention paid to the pressure exerted on "confining walls". We employ a mathematical connection between the local density of OUPs and the statistics of their propulsion force to demonstrate the existence of an equation of state in one dimension. In higher dimensions we show that active particles generate a non-conservative force field in the surrounding medium. A simplified far-from-equilibrium model is proposed to account for OUP behaviour in the vicinity of potentials. Building on this, we interpret simulations of OUPs in more complicated situations involving asymmetrical and spatially curved potentials, characterising the inhomogeneous local stresses which result in terms of competing active length-scales.


title: Diamonds, Compactness, and Measure Sequences

abstract: We establish the consistency of the failure of the diamond principle on a cardinal $\kappa$ which satisfies a strong simultaneous reflection property. The result is based on an analysis of Radin forcing, and further leads to a characterization of weak compactness of $\kappa$ in a Radin generic extension.


title: Evidence for an intermediate-mass black hole in the globular cluster NGC 6624

abstract: PSR B1820$-$30A is located in the globular cluster NGC 6624 and is the closest known pulsar to the centre of any globular cluster. We present more than 25 years of high-precision timing observations of this millisecond pulsar and obtain four rotational frequency time derivative measurements. Modelling these higher-order derivatives as being due to orbital motion, we find solutions which indicate that the pulsar is in either a low-eccentricity ($0.33\lesssim e\lesssim0.4$) smaller orbit with a low mass companion (such as a main sequence star, white dwarf, neutron star, or stellar mass black hole) or a high-eccentricity ($e\gtrsim0.9$) larger orbit with a massive companion. The cluster mass properties and the observed properties of 4U 1820$-$30 and the other pulsars in the cluster argue against the low-eccentricity possibility. The high-eccentricity solution reveals that the pulsar is most likely orbiting around an intermediate-mass black hole (IMBH) of mass $> 7,500$~M$_\odot$ located at the cluster centre. A gravitational model for the globular cluster, which includes such a central black hole (BH), predicts an acceleration that is commensurate with that measured for the pulsar. It further predicts that the model-dependent minimum mass of the IMBH is $\sim60,000$~M$_\odot$. Accounting for the associated contribution to the observed period derivative indicates that the $\gamma$-ray efficiency of the pulsar should be between 0.08 and 0.2. Our results suggest that other globular clusters may also contain central black holes and they may be revealed by the study of new pulsars found sufficiently close to their centres. Note that we found an erratum in Section 5 and thus, the $\sim$60,000~M$_\odot$ mass mentioned above has to be replaced by the correct model-dependent mass limit of $\sim$20,000~M$_\odot$. See the erratum appended.


title: Automatically Identifying Fake News in Popular Twitter Threads

abstract: Information quality in social media is an increasingly important issue, but web-scale data hinders experts' ability to assess and correct much of the inaccurate content, or `fake news,' present in these platforms. This paper develops a method for automating fake news detection on Twitter by learning to predict accuracy assessments in two credibility-focused Twitter datasets: CREDBANK, a crowdsourced dataset of accuracy assessments for events in Twitter, and PHEME, a dataset of potential rumors in Twitter and journalistic assessments of their accuracies. We apply this method to Twitter content sourced from BuzzFeed's fake news dataset and show models trained against crowdsourced workers outperform models based on journalists' assessment and models trained on a pooled dataset of both crowdsourced workers and journalists. All three datasets, aligned into a uniform format, are also publicly available. A feature analysis then identifies features that are most predictive for crowdsourced and journalistic accuracy assessments, results of which are consistent with prior work. We close with a discussion contrasting accuracy and credibility and why models of non-experts outperform models of journalists for fake news detection in Twitter.


title: A Generalized Labeled Multi-Bernoulli Filter with Object Spawning

abstract: Previous labeled random finite set filter developments use a motion model that only accounts for survival and birth. While such a model provides the means for a multi-object tracking filter such as the Generalized Labeled Multi-Bernoulli (GLMB) filter to capture object births and deaths in a wide variety of applications, it lacks the capability to capture spawned tracks and their lineages. In this paper, we propose a new GLMB based filter that formally incorporates spawning, in addition to birth. This formulation enables the joint estimation of a spawned object's state and information regarding its lineage. Simulations results demonstrate the efficacy of the proposed formulation.


title: Rhythmic production of consonant-vowel syllables synchronizes traveling waves in speech-processing brain regions

abstract: Nature is abundant in oscillatory activity, with oscillators that have the remarkable ability of synchronizing to external events. Using electrocorticographic (ECoG) recordings from a subject rhythmically producing consonant-vowel syllables (CVSs) we show that neural oscillators recorded at individual ECoG electrodes become precisely synchronized to initiations of the production of CVSs (i.e., that these initiations occur at precise phases of bandpassed-filtered voltages recorded at most ECoG electrodes). This synchronization is not a trivial consequence of the rhythmic production of CVSs, since it takes several minutes to be fully established and is observed at the frequency of CVS production and at its second harmonic. The phase of filtered voltages at which CVSs are produced varies systematically across the grid of electrodes, consistently with the propagation of traveling waves (TWs). Using these synchronized phases we isolate a first TW in voltages (filtered at the median CVS-production frequency) moving from primary auditory to premotor cortex, and a second TW in high-gamma amplitude (coupled to phase at the CVS-production frequency) moving along the same path but in opposite direction. To our knowledge, this is the first report of rhythmic motor acts synchronizing spatio-temporally organized cortical activity in the human brain.


title: Strong solutions of SDE's with generalized drift and multidimensional fractional Brownian initial noise

abstract: In this paper we prove the existence of strong solutions to a SDE with a generalized drift driven by a multidimensional fractional Brownian motion for small Hurst parameters H<1/2. Here the generalized drift is given as the local time of the unknown solution process, which can be considered an extension of the concept of a skew Brownian motion to the case of fractional Brownian motion. Our approach for the construction of strong solutions is new and relies on techniques from Malliavin calculus combined with a "local time variational calculus" argument.


title: Wallpaper Fermions and the Nonsymmorphic Dirac Insulator

abstract: Recent developments in the relationship between bulk topology and surface crystal symmetry have led to the discovery of materials whose gapless surface states are protected by crystal symmetries. In fact, there exists only a very limited set of possible surface crystal symmetries, captured by the 17 "wallpaper groups." We show that a consideration of symmetry-allowed band degeneracies in the wallpaper groups can be used to understand previous topological crystalline insulators, as well as to predict new examples. In particular, the two wallpaper groups with multiple glide lines, $pgg$ and $p4g$, allow for a new topological insulating phase, whose surface spectrum consists of only a single, fourfold-degenerate, true Dirac fermion. Like the surface state of a conventional topological insulator, the surface Dirac fermion in this "nonsymmorphic Dirac insulator" provides a theoretical exception to a fermion doubling theorem. Unlike the surface state of a conventional topological insulator, it can be gapped into topologically distinct surface regions while keeping time-reversal symmetry, allowing for networks of topological surface quantum spin Hall domain walls. We report the theoretical discovery of new topological crystalline phases in the A$_2$B$_3$ family of materials in SG 127, finding that Sr$_2$Pb$_3$ hosts this new topological surface Dirac fermion. Furthermore, (100)-strained Au$_2$Y$_3$ and Hg$_2$Sr$_3$ host related topological surface hourglass fermions. We also report the presence of this new topological hourglass phase in Ba$_5$In$_2$Sb$_6$ in SG 55. For orthorhombic space groups with two glides, we catalog all possible bulk topological phases by a consideration of the allowed non-abelian Wilson loop connectivities, and we develop topological invariants for these systems. Finally, we show how in a particular limit, these crystalline phases reduce to copies of the SSH model.


title: Magnetic properties of Dy nano-islands on graphene

abstract: We have determined the magnetic properties of epitaxially grown Dy islands on graphene/SiC(0001) that are passivated by a gold film (deposited in the ultra-high vacuum growth chamber) for {\it ex-situ} X-ray magnetic circular dichroism (XMCD). Our sum-rule analysis of the Dy $M_{4,5}$ XMCD spectra at low temperatures ($T=15$ K) as a function of magnetic field assuming Dy$^{3+}$ (spin configuration $^6H_{15/2}$) indicate that the projection of the magnetic moment along an applied magnetic field of 5 T is 3.5(3) $\mu_B$. Temperature dependence of the magnetic moment (extracted from the $M_5$ XMCD spectra) shows an onset of a change in magnetic moment at about 175 K in proximity of the transition from paramagnetic to helical magnetic structure at $T_{\rm H} =179$ K in bulk Dy. No feature at the vicinity of the ferromagnetic transition of hcp bulk Dy at $T_{\rm c}$ = 88 K is observed. However, below $\sim$130 K, the inverse magnetic moment (extracted from the XMCD) is linear in temperature as commonly expected from a paramagnetic system suggesting different behavior of Dy nano-island than bulk Dy.


title: Complex dimensions and their observability

abstract: We show that the dimension of spacetime becomes complex-valued when its short-scale geometry is invariant under a discrete scaling symmetry. This characteristic can generically arise in quantum gravities, for instance, in those based on combinatorial or multifractal structures or as the partial breaking of continuous dilation symmetry in any conformal-invariant theory. With its infinite scale hierarchy, discrete scale invariance overlaps with the traditional separation between ultraviolet and infrared physics and it can leave an all-range observable imprint, such as a pattern of log oscillations and sharp features in the cosmic microwave background primordial power spectrum.


title: Improving Nanopore Reads Raw Signal Alignment

abstract: We investigate usage of dynamic time warping (DTW) algorithm for aligning raw signal data from MinION sequencer. DTW is mostly using for fast alignment for selective sequencing to quickly determine whether a read comes from sequence of interest. We show that standard usage of DTW has low discriminative power mainly due to problem with accurate estimation of scaling parameters. We propose a simple variation of DTW algorithm, which does not suffer from scaling problems and has much higher discriminative power.


title: Existence of Infinite Product Measures

abstract: A construction of product measures is given for an arbitrary sequence of measure spaces via outer measure techniques without imposing any condition on the underlying measure spaces. This approach concludes finally the problem of the existence of product measures in an elementary manner. Moreover, the $L_{p}$ spaces of this measures are simplified in terms of finite product measures following the approach of [21]. This decomposition simplifies infinite dimensional integration and gives to this theory a computational framework.


title: Stabilization of the wave equation with moving boundary

abstract: We deal with the wave equation with assigned moving boundary ($0<x<a(t)$) upon which Dirichlet-Neuman boundary conditions are satisfied, here $a(t)$ is assumed to move slower than the light and periodically. We give a feedback which guarantees the exponential decay of the energy. The proof relies on a reduction theorem of Yoccoz. At the end we give a remark on the moving-pointwise stabilization problem.


title: Dual numbers, weighted quivers, and extended Somos and Gale-Robinson sequences

abstract: We investigate a general method that allows one to construct new integer sequences extending existing ones. We apply this method to the classic Somos-4 and Somos-5, and the Gale-Robinson sequences, as well as to more general class of sequences introduced by Fordy and Marsh, and produce a great number of new sequences. The method is based on the notion of "weighted quiver", a quiver with a $\mathbb Z$-valued function on the set of vertices that obeys very special rules of mutation.


title: Distributed generalized Nash equilibria computation of monotone games via a preconditioned proximal point algorithm

abstract: In this paper, we investigate distributed generalized Nash equilibrium (GNE) computation of monotone games with affine coupling constraints. Each player can only utilize its local objective function, local feasible set and a local block of the coupling constraint, and can only communicate with its neighbours. We assume the game has monotone pseudo-subdifferential without Lipschitz continuity restrictions. We design novel center-free distributed GNE seeking algorithms for equality and inequality affine coupling constraints, respectively. A proximal alternating direction method of multipliers(ADMM) is proposed for the equality case, while for the inequality case, a parallel splitting type algorithm is proposed. In both algorithms, the GNE seeking task is decomposed into a sequential NE computation of regularized subgames and distributed update of multipliers and auxiliary variables, based on local data and local communication. Our two double-layer GNE algorithms need not specify the inner-loop NE seeking algorithm and moreover, only require that the strongly monotone subgames are inexactly solved. We prove their convergence by showing that the two algorithms can be seen as specific instances of preconditioned proximal point algorithms} (PPPA) for finding zeros of monotone operators. Applications and numerical simulations are given for illustration.


title: SIMP J013656.5+093347 is Likely a Planetary-Mass Object in the Carina-Near Moving Group

abstract: We report the discovery that the nearby ($\sim$ 6 pc) photometrically variable T2.5 dwarf SIMP~J013656.5+093347 is a likely member of the $\sim$ 200 Myr-old Carina-Near moving group with a probability of $>$ 99.9% based on its full kinematics. Our $v\sin i$ measurement of $50.9 \pm 0.8$ km s$^{-1}$ combined with the known rotation period inferred from variability measurements provide a lower limit of $1.01 \pm 0.02$ $R_{\mathrm{Jup}}$ on the radius of SIMP0136+0933, an independent verification that it must be younger than $\sim$ 950 Myr according to evolution models. We estimate a field interloper probability of 0.2% based on the density of field T0-T5 dwarfs. At the age of Carina-Near, SIMP0136+0933 has an estimated mass of $12.7 \pm 1.0$ $M_{\mathrm{Jup}}$ and is predicted to have burned roughly half of its original deuterium. SIMP0136+0933 is the closest known young moving group member to the Sun, and is one of only a few known young T dwarfs, making it an important benchmark for understanding the atmospheres of young planetary-mass objects.


title: Compressing DMA Engine: Leveraging Activation Sparsity for Training Deep Neural Networks

abstract: Popular deep learning frameworks require users to fine-tune their memory usage so that the training data of a deep neural network (DNN) fits within the GPU physical memory. Prior work tries to address this restriction by virtualizing the memory usage of DNNs, enabling both CPU and GPU memory to be utilized for memory allocations. Despite its merits, virtualizing memory can incur significant performance overheads when the time needed to copy data back and forth from CPU memory is higher than the latency to perform the computations required for DNN forward and backward propagation. We introduce a high-performance virtualization strategy based on a "compressing DMA engine" (cDMA) that drastically reduces the size of the data structures that are targeted for CPU-side allocations. The cDMA engine offers an average 2.6x (maximum 13.8x) compression ratio by exploiting the sparsity inherent in offloaded data, improving the performance of virtualized DNNs by an average 32% (maximum 61%).


title: Two dimensional solitary waves in shear flows

abstract: In this paper we study existence and asymptotic behavior of solitary-wave solutions for the generalized Shrira equation, a two-dimensional model appearing in shear flows. The method used to show the existence of such special solutions is based on the mountain pass theorem. One of the main difficulties consists in showing the compact embedding of the energy space in the Lesbesgue spaces; this is dealt with interpolation theory. Regularity and decay properties of the solitary waves are also established.


title: Quasi-Automorphism Groups of Type F-infinity

abstract: The groups QF, QT, and QV are groups of quasi-automorphisms of the infinite binary tree. Their names indicate a similarity with Thompson's well-known groups F, T, and V. We will use the theory of diagram groups over semigroup presentations to prove that all of the above groups (and several generalizations) have type F-infinity. Our proof uses certain types of hybrid diagrams, which have properties in common with both planar diagrams and braided diagrams. The diagram groups defined by hybrid diagrams also act properly and isometrically on CAT(0) cubical complexes.


title: A Formal Semantics for Data Analytics Pipelines

abstract: In this report, we present a new programming model based on Pipelines and Operators, which are the building blocks of programs written in PiCo, a DSL for Data Analytics Pipelines. In the model we propose, we use the term Pipeline to denote a workflow that processes data collections -- rather than a computational process -- as is common in the data processing community. The novelty with respect to other frameworks is that all PiCo operators are polymorphic with respect to data types. This makes it possible to 1) re-use the same algorithms and pipelines on different data models (e.g., streams, lists, sets, etc); 2) reuse the same operators in different contexts, and 3) update operators without affecting the calling context, i.e., the previous and following stages in the pipeline. Notice that in other mainstream frameworks, such as Spark, the update of a pipeline by changing a transformation with another is not necessarily trivial, since it may require the development of an input and output proxy to adapt the new transformation for the calling context. In the same line, we provide a formal framework (i.e., typing and semantics) that characterizes programs from the perspective of how they transform the data structures they process -- rather than the computational processes they represent. This approach allows to reason about programs at an abstract level, without taking into account any aspect from the underlying execution model or implementation.


title: Cusp-shaped Elastic Creases and Furrows

abstract: The surfaces of growing biological tissues, swelling gels, and compressed rubbers do not remain smooth, but frequently exhibit highly localized inward folds. We reveal the morphology of this surface folding in a novel experimental setup, which permits to deform the surface of a soft gel in a controlled fashion. The interface first forms a sharp furrow, whose tip size decreases rapidly with deformation. Above a critical deformation, the furrow bifurcates to an inward folded crease of vanishing tip size. We show experimentally and numerically that both creases and furrows exhibit a universal cusp-shape, whose width scales like $y^{3/2}$ at a distance $y$ from the tip. We provide a similarity theory that captures the singular profiles before and after the self-folding bifurcation, and derive the length of the fold from large deformation elasticity.


title: Pressure and Flow of Exponentially Self-Correlated Active Particles

abstract: Microscopic swimming particles, which dissipate energy to execute persistent directed motion, are a classic example of a non-equilibrium system. We investigate the non-interacting Ornstein--Uhlenbeck Particle (OUP), which is propelled through a viscous medium by a force which is correlated over a finite time. We obtain an exact expression for the steady state phase-space density of a single OUP confined by a quadratic potential, and use the result to explore more complex geometries, both through analytical approximations and numerical simulations. In a "Casimir"-style setup involving two narrowly-spaced walls, we describe a particle-trapping phenomenon, which leads to a repulsive effective interaction between the walls; while in a two-dimensional annulus geometry, we observe net stresses which resemble the Laplace pressure.


title: Electron doped layered nickelates: spanning the phase diagram of the cuprates

abstract: Pr$_4$Ni$_3$O$_8$ is an overdoped analog of hole-doped layered cuprates. Here we show via ab initio calculations that Ce-doped Pr$_4$Ni$_3$O$_8$ (Pr$_3$CeNi$_3$O$_8$) has the same electronic structure as the antiferromagnetic insulating phase of parent cuprates. We find that substantial Ce-doping should be thermodynamically stable and that other 4+ cations would yield a similar antiferromagnetic insulating state, arguing this configuration is robust for layered nickelates of low enough valence. The analogies with cuprates at different $d$ fillings suggest that intermediate Ce-doping concentrations near 1/8 should be an appropriate place to search for superconductivity in these low-valence Ni oxides.


title: Validity of the Weizs\"acker-Williams Approximation and the Analysis of Beam Dump Experiments: Production of an axion, a dark photon, or a new axial-vector boson

abstract: Beam dump experiments have been used to search for new particles, $\phi$, with null results interpreted in terms of limits on masses $m_\phi$ and coupling constants $\epsilon$. However these limits have been obtained by using approximations [including the Weizs\"{a}cker-Williams (WW) approximation] or Monte-Carlo simulations. We display methods to obtain the cross section and the resulting particle production rates without using approximations on the phase space integral or Monte-Carlo simulations. In our previous work we examined the case of the new scalar boson production; in this paper we explore all possible new spin-0 and spin-1 particles. We show that the approximations cannot be used to obtain accurate values of cross sections. The corresponding exclusion plots differ by substantial amounts when seen on a linear scale. Furthermore, a new region ($m_\phi<2m_e$) of parameter space can be explored without using one of the common approximations, $m_\phi\gg m_e$. We derive new expressions for the three photon decays of dark photon and four photon decays of new axial-vector bosons. As a result, the production cross section and exclusion region of different low mass ($m_\phi<2m_e$) bosons are very different. Moreover, our method can be used as a consistency check for Monte-Carlo simulations.


title: Galaxy Protoclusters as Drivers of Cosmic Star-Formation History in the First 2 Gyr

abstract: Present-day clusters are massive halos containing mostly quiescent galaxies, while distant protoclusters are extended structures containing numerous star-forming galaxies. We investigate the implications of this fundamental change in a cosmological context using a set of N-body simulations and semi-analytic models. We find that the fraction of the cosmic volume occupied by all (proto)clusters increases by nearly three orders of magnitude from z=0 to z=7. We show that (proto)cluster galaxies are an important, and even dominant population at high redshift, as their expected contribution to the cosmic star-formation rate density rises (from 1% at z=0) to 20% at z=2 and 50% at z=10. Protoclusters thus provide a significant fraction of the cosmic ionizing photons, and may have been crucial in driving the timing and topology of cosmic reionization. Internally, the average history of cluster formation can be described by three distinct phases: at z~10-5, galaxy growth in protoclusters proceeded in an inside-out manner, with centrally dominant halos that are among the most active regions in the Universe; at z~5-1.5, rapid star formation occurred within the entire 10-20 Mpc structures, forming most of their present-day stellar mass; at z<~1.5, violent gravitational collapse drove these stellar contents into single cluster halos, largely erasing the details of cluster galaxy formation due to relaxation and virialization. Our results motivate observations of distant protoclusters in order to understand the rapid, extended stellar growth during Cosmic Noon, and their connection to reionization during Cosmic Dawn.


title: Reconstructing rational stable motivic homotopy theory

abstract: Using a recent computation of the rational minus part of $SH(k)$ by Ananyevskiy-Levine-Panin, a theorem of Cisinski-Deglise and a version of the Roendigs-Ostvaer theorem, rational stable motivic homotopy theory over an infinite perfect field of characteristic different from 2 is recovered in this paper from finite Milnor-Witt correspondences in the sense of Calmes-Fasel.


title: Study of secondary neutron interactions with $^{232}$Th, $^{129}$I, and $^{127}$I nuclei with the uranium assembly "QUINTA" at 2, 4, and 8 GeV deuteron beams of the JINR Nuclotron accelerator

abstract: The natural uranium assembly, "QUINTA", was irradiated with 2, 4, and 8 GeV deuterons. The $^{232}$Th, $^{127}$I, and $^{129}$I samples have been exposed to secondary neutrons produced in the assembly at a 20-cm radial distance from the deuteron beam axis. The spectra of gamma rays emitted by the activated $^{232}$Th, $^{127}$I, and $^{129}$I samples have been analyzed and several tens of product nuclei have been identified. For each of those products, neutron-induced reaction rates have been determined. The transmutation power for the $^{129}$I samples is estimated. Experimental results were compared to those calculated with well-known stochastic and deterministic codes.


title: Asymptotically optimal bound on the adjacent vertex distinguishing edge choice number

abstract: An adjacent vertex distinguishing edge colouring of a graph $G$ without isolated edges is its proper edge colouring such that no pair of adjacent vertices meets the same set of colours in $G$. We show that such colouring can be chosen from any set of lists associated to the edges of $G$ as long as the size of every list is at least $\Delta+C\Delta^{\frac{1}{2}}(\log\Delta)^4$, where $\Delta$ is the maximum degree of $G$ and $C$ is a constant. The proof is probabilistic. The same is true in the environment of total colourings.


title: Longitudinal negative magnetoresistance and magneto-transport phenomena in conventional and topological conductors

abstract: Recently a large negative longitudinal (parallel to the magnetic field) magnetoresistance was observed in Weyl and Dirac semimetals. It is believed to be related to the chiral anomaly associated with topological electron band structure of these materials. We show that in a certain range of parameters such a phenomenon can also exist in conventional centrosymmetric and time reversal conductors, lacking topological protection of the electron spectrum and the chiral anomaly. We also discuss the magnetic field enhancement of the longitudinal components of the thermal conductivity and thermoelectric tensors.


title: Gaiotto's Lagrangian subvarieties via loop groups

abstract: The purpose of this note is to give a simple proof of the fact that a certain substack, defined in [2], of the moduli stack $T^{\ast}Bun_G(\Sigma)$ of Higgs bundles over a curve $\Sigma$, for a connected, simply connected semisimple group $G$, possesses a Lagrangian structure. The substack, roughly speaking, consists of images under the moment map of global sections of principal $G$-bundles over $\Sigma$ twisted by a smooth symplectic variety with a Hamiltonian $G$-action.


title: State-Dependent Gaussian Multiple Access Channels: New Outer Bounds and Capacity Results

abstract: This paper studies a two-user state-dependent Gaussian multiple-access channel (MAC) with state noncausally known at one encoder. Two scenarios are considered: i) each user wishes to communicate an independent message to the common receiver, and ii) the two encoders send a common message to the receiver and the non-cognitive encoder (i.e., the encoder that does not know the state) sends an independent individual message (this model is also known as the MAC with degraded message sets). For both scenarios, new outer bounds on the capacity region are derived, which improve uniformly over the best known outer bounds. In the first scenario, the two corner points of the capacity region as well as the sum rate capacity are established, and it is shown that a single-letter solution is adequate to achieve both the corner points and the sum rate capacity. Furthermore, the full capacity region is characterized in situations in which the sum rate capacity is equal to the capacity of the helper problem. The proof exploits the optimal-transportation idea of Polyanskiy and Wu (which was used previously to establish an outer bound on the capacity region of the interference channel) and the worst-case Gaussian noise result for the case in which the input and the noise are dependent.


title: A combinatorial approach to Rauzy-type dynamics I: permutations and the Kontsevich--Zorich--Boissy classification theorem

abstract: Rauzy-type dynamics are group actions on a collection of combinatorial objects. The first and best known example concerns an action on permutations, associated to interval exchange transformations (IET) for the Poincar\'e map on compact orientable translation surfaces. The equivalence classes on the objects induced by the group action are related to components of the moduli spaces of Abelian differentials with prescribed singularities, and, in two variants of the problem, have been classified by Kontsevich and Zorich, and by Boissy, through methods involving both combinatorics and algebraic geometry. We provide here a purely combinatorial proof of both classification theorems, and in passing establish a few previously unnoticed features. As will be shown elsewhere, our methods extend also to other Rauzy-type dynamics, both on labeled and unlabeled structures. Some of these dynamics have a geometrical interpretation (e.g., matchings, related to IET on non-orientable surfaces), while some others do not have one so far.


title: When is the Chernoff Exponent for Quantum Operations finite?

abstract: We consider the problem of testing two hypotheses of quantum operations in a setting of many uses where an arbitrary prior probability distribution is given. The Chernoff exponent for quantum operations is investigated to track the minimal average error probability of discriminating two quantum operations asymptotically. We answer the question, "When is the Chernoff exponent for quantum operations finite?" We show that either two quantum operations can be perfectly distinguished with finite uses, or the minimal discrimination error decays exponentially with respect to the number of uses asymptotically. That is, the Chernoff exponent is finite if and only if the quantum operations can not be perfectly distinguished with finite uses. This rules out the possibility of super-exponential decay of error probability. Upper bounds of the Chernoff exponent for quantum operations are provided.


title: Prefer Nested Segmentation to Compound Segmentation

abstract: Introduction: Intra-organ radiation dose sensitivity is becoming increasingly relevant in clinical radiotherapy. One method for assessment involves partitioning delineated regions of interest and comparing the relative contributions or importance to clinical outcomes. We show that an intuitive method for dividing organ contours, compound (sub-)segmentation, can unintentionally lead to sub-segments with inconsistent volumes, which will bias relative importance assessment. An improved technique, nested segmentation, is introduced and compared. Methods: Clinical radiotherapy planning parotid contours from 510 patients were segmented. Counts of radiotherapy dose matrix voxels interior to sub-segments were used to determine the equivalency of sub-segment volumes. The distribution of voxel counts within sub-segments were compared using Kolmogorov-Smirnov tests and characterized by their dispersion. Analytical solutions for 2D/3D analogues were derived and sub-segment area/volume were compared directly. Results: Both parotid and 2D/3D region of interest analogue segmentation confirmed compound segmentation intrinsically produces sub-segments with volumes that depend on the region of interest shape and selection location. Significant volume differences were observed when sub-segmenting parotid contours into 18ths, and vanishingly small sub-segments were observed when sub-segmenting into 96ths. Central sub-segments were considerably smaller than sub-segments on the periphery. Nested segmentation did not exhibit these shortcomings and produced sub-segments with equivalent volumes when dose grid and contour collinearity was addressed, even when dividing the parotid into 96ths. Nested segmentation was always faster or equivalent in runtime to compound segmentation. Conclusions: Nested segmentation is more suited than compound segmentation for analyses requiring equal weighting of sub-segments.


title: Combinatorial Auctions Do Need Modest Interaction

abstract: We study the necessity of interaction for obtaining efficient allocations in subadditive combinatorial auctions. This problem was originally introduced by Dobzinski, Nisan, and Oren (STOC'14) as the following simple market scenario: $m$ items are to be allocated among $n$ bidders in a distributed setting where bidders valuations are private and hence communication is needed to obtain an efficient allocation. The communication happens in rounds: in each round, each bidder, simultaneously with others, broadcasts a message to all parties involved and the central planner computes an allocation solely based on the communicated messages. Dobzinski et.al. showed that no non-interactive ($1$-round) protocol with polynomial communication (in the number of items and bidders) can achieve approximation ratio better than $\Omega(m^{{1}/{4}})$, while for any $r \geq 1$, there exists $r$-round protocols that achieve $\widetilde{O}(r \cdot m^{{1}/{r+1}})$ approximation with polynomial communication; in particular, $O(\log{m})$ rounds of interaction suffice to obtain an (almost) efficient allocation. A natural question at this point is to identify the "right" level of interaction (i.e., number of rounds) necessary to obtain an efficient allocation. In this paper, we resolve this question by providing an almost tight round-approximation tradeoff for this problem: we show that for any $r \geq 1$, any $r$-round protocol that uses polynomial communication can only approximate the social welfare up to a factor of $\Omega(\frac{1}{r} \cdot m^{{1}/{2r+1}})$. This in particular implies that $\Omega(\frac{\log{m}}{\log\log{m}})$ rounds of interaction are necessary for obtaining any efficient allocation in these markets. Our work builds on the recent multi-party round-elimination technique of Alon, Nisan, Raz, and Weinstein (FOCS'15) and settles an open question posed by Dobzinski et.al. and Alon et. al.


title: Vertex algebras and 4-manifold invariants

abstract: We propose a way of computing 4-manifold invariants, old and new, as chiral correlation functions in half-twisted 2d $\mathcal{N}=(0,2)$ theories that arise from compactification of fivebranes. Such formulation gives a new interpretation of some known statements about Seiberg-Witten invariants, such as the basic class condition, and gives a prediction for structural properties of the multi-monopole invariants and their non-abelian generalizations.


title: Recursive Integral Method with Cayley Transformation

abstract: Recently, a non-classical eigenvalue solver, called RIM, was proposed to compute (all) eigenvalues in a region on the complex plane. Without solving any eigenvalue problem, it tests if a region contains eigenvalues using an approximate spectral projection. Regions that contain eigenvalues are subdivided and tested recursively until eigenvalues are isolated with a specified precision. This makes RIM an eigensolver distinct from all existing methods. Furthermore, it requires no a priori spectral information. In this paper, we propose an improved version of {\bf RIM} for non-Hermitian eigenvalue problems. Using Cayley transformation and Arnoldi's method, the computation cost is reduced significantly. Effectiveness and efficiency of the new method are demonstrated by numerical examples and compared with 'eigs' in Matlab.


title: Integral Gassman equivalence of algebraic and hyperbolic manifolds

abstract: In this paper we construct arbitrarily large families of smooth projective varieties and closed Riemannian manifolds that share many algebraic and analytic invariants. For instance, every non-arithmetic, closed hyperbolic $3$--manifold admits arbitrarily large collections of non-isometric finite covers which are strongly isospectral, length isospectral, and have isomorphic integral cohomology where the isomorphisms commute with restriction and co-restriction. We can also construct arbitrarily large collections of pairwise non-isomorphic smooth projective surfaces where these isomorphisms in cohomology are natural with respect to Hodge structure or as Galois modules. In particular, the projective varieties have isomorphic Picard and Albanese varieties, and they also have isomorphic effective Chow motives. Our construction employs an integral refinement of the Gassman--Sunada construction that has recently been utilized by D. Prasad. One application of our work shows the non-injectivity of the map from the Grothendieck group of varieties over $\overline{\mathbf{Q}}$ to the Grothendieck group of the category of effective Chow motives. We also answer a question of D. Prasad.


title: Decomposing graphs into forests

abstract: We give a simple graph-theoretic proof of a classical result due to C. St. J. A. Nash-Williams on covering graphs by forests. Moreover we derive a slight generalisation of this statement where some edges are preassigned to distinct forests.


title: Solvent Fluctuations and Nuclear Quantum Effects Modulate the Molecular Hyperpolarizability of Water

abstract: Second-Harmonic Scatteringh (SHS) experiments provide a unique approach to probe non-centrosymmetric environments in aqueous media, from bulk solutions to interfaces, living cells and tissue. A central assumption made in analyzing SHS experiments is that the each molecule scatters light according to a constant molecular hyperpolarizability tensor $\boldsymbol{\beta}^{(2)}$. Here, we investigate the dependence of the molecular hyperpolarizability of water on its environment and internal geometric distortions, in order to test the hypothesis of constant $\boldsymbol{\beta}^{(2)}$. We use quantum chemistry calculations of the hyperpolarizability of a molecule embedded in point-charge environments obtained from simulations of bulk water. We demonstrate that both the heterogeneity of the solvent configurations and the quantum mechanical fluctuations of the molecular geometry introduce large variations in the non-linear optical response of water. This finding has the potential to change the way SHS experiments are interpreted: in particular, isotopic differences between H$_2$O and D$_2$O could explain recent second-harmonic scattering observations. Finally, we show that a simple machine-learning framework can predict accurately the fluctuations of the molecular hyperpolarizability. This model accounts for the microscopic inhomogeneity of the solvent and represents a first step towards quantitative modelling of SHS experiments.


title: Chiral vortical effect in pionic superfluid vs spin alignment of baryons

abstract: We consider chiral fluids, with (nearly) massless fermionic constituents, in the confining phase. Chiral vortical effect (CVE) is the flow of axial current along the axis of rotation of the fluid while the spin alignment is a non-vanishing correlation of polarizations of baryons with the axis of rotation. As the theoretical framework we use the model of pionic superfluidity induced by a non-vanishing isotopic chemical potential. We note that the average value of spin of virtual baryons reproduces the CVE. The role of defects, or vortices is crucial. The model does not apply directly to the quark-gluon plasma but might indicate existence of a mechanism to produce baryons with relatively large polarization in heavy-ion collisions.


title: Modeling temporal constraints for a system of interactive scores

abstract: In this chapter we explain briefly the fundamentals of the interactive scores formalism. Then we develop a solution for implementing the ECO machine by mixing petri nets and constraints propagation. We also present another solution for implementing the ECO machine using concurrent constraint programming. Finally, we present an extension of interactive score with conditional branching.


title: Polluted Bootstrap Percolation with Threshold Two in All Dimensions

abstract: In the polluted bootstrap percolation model, the vertices of a graph are independently declared initially occupied with probability p or closed with probability q. At subsequent steps, a vertex becomes occupied if it is not closed and it has at least r occupied neighbors. On the cubic lattice Z^d of dimension d>=3 with threshold r=2, we prove that the final density of occupied sites converges to 1 as p and q both approach 0, regardless of their relative scaling. Our result partially resolves a conjecture of Morris, and contrasts with the d=2 case, where Gravner and McDonald proved that the critical parameter is q/{p^2}.


title: Eigenspaces of symmetric graphs are not typically irreducible

abstract: We construct rich families of Schr\"odinger operators on symmetric graphs, both quantum and combinatorial, whose spectral degeneracies are persistently larger than the maximal dimension of an irreducible representations of the symmetry group.


title: Are Unobservables Separable?

abstract: It is common to assume in empirical research that observables and unobservables are additively separable, especially, when the former are endogenous. This is done because it is widely recognized that identification and estimation challenges arise when interactions between the two are allowed for. Starting from a nonseparable IV model, where the instrumental variable is independent of unobservables, we develop a novel nonparametric test of separability of unobservables. The large-sample distribution of the test statistics is nonstandard and relies on a novel Donsker-type central limit theorem for the empirical distribution of nonparametric IV residuals, which may be of independent interest. Using a dataset drawn from the 2015 US Consumer Expenditure Survey, we find that the test rejects the separability in Engel curves for most of the commodities.


title: Zitterbewegung in time-reversal Weyl Semimetals

abstract: We perform a systematic study of the Zitterbewegung effect of fermions, which are described by a Gaussian wave with broken spatial-inversion symmetry in a three-dimensional low-energy Weyl semimetal. Our results show that the motion of fermions near the Weyl points is characterized by rectilinear motion and Zitterbewegung oscillation. The ZB oscillation is affected by the width of the Gaussian wave packet, the position of the Weyl node, and the chirality and anisotropy of the fermions. By introducing a one-dimensional cosine potential, the new generated massless fermions have lower Fermi Velocities, which results in a robust relativistic oscillation. Modulating the height and periodicity of periodic potential demonstrates that the ZB effect of fermions in the different Brillouin zones exhibits quasi-periodic behavior. These results may provide an appropriate system for probing the Zitterbewegung effect experimentally.


title: The portrait of eikonal instability in Lovelock theories

abstract: Perturbations and eikonal instabilities of black holes and branes in the Einstein-Gauss-Bonnet theory and its Lovelock generalization were considered in the literature for several particular cases, where the asymptotic conditions (flat, dS, AdS), the number of spacetime dimensions $D$, non-vanishing coupling constants ($\alpha_1$, $\alpha_2$, $\alpha_3$ etc.) and other parameters have been chosen in a specific way. Here we give a comprehensive analysis of the eikonal instabilities of black holes and branes for the \emph{most general} Lovelock theory, not limited by any of the above cases. Although the part of the stability analysis is performed here purely analytically and formulated in terms of the inequalities for the black hole parameters, the most general case is treated numerically and the accurate regions of instabilities are presented. The shared Mathematica(R) code allows the reader to construct the regions of eikonal instability for any desired values of the parameters.


title: The Brauer group of $\mathscr{M}_{1,1}$ over algebraically closed fields of characteristic $2$

abstract: We prove that the Brauer group of the moduli stack of elliptic curves $\mathscr{M}_{1,1,k}$ over an algebraically closed field $k$ of characteristic $2$ is isomorphic to $\mathbb{Z}/(2)$. We also compute the Brauer group of $\mathscr{M}_{1,1,k}$ where $k$ is a finite field of characteristic $2$.


title: The Stability of Full Dimensional KAM tori for Nonlinear Schr\"odinger equation

abstract: In this paper, it is proved that the full dimensional invariant tori obtained by Bourgain [J. Funct. Anal., \textbf{229} (2005), no. 1, 62-94.] is stable in a very long time for 1D nonlinear Schr\"{o}dinger equation with periodic boundary conditions.


title: Graded hypoellipticity of BGG sequences

abstract: This article studies hypoellipticity on general filtered manifolds. We extend the Rockland criterion to a pseudodifferential calculus on filtered manifolds, construct a parametrix and describe its precise analytic structure. We use this result to study Rockland sequences, a notion generalizing elliptic sequences to filtered manifolds. The main application that we present is to the analysis of the Bernstein--Gelfand--Gelfand (BGG) sequences over regular parabolic geometries. We do this by generalizing the BGG machinery to more general filtered manifolds (in a non-canonical way) and show that the generalized BGG sequences are Rockland in a graded sense.


title: MapReduce Particle Filtering with Exact Resampling and Deterministic Runtime

abstract: Particle filtering is a numerical Bayesian technique that has great potential for solving sequential estimation problems involving non-linear and non-Gaussian models. Since the estimation accuracy achieved by particle filters improves as the number of particles increases, it is natural to consider as many particles as possible. MapReduce is a generic programming model that makes it possible to scale a wide variety of algorithms to Big data. However, despite the application of particle filters across many domains, little attention has been devoted to implementing particle filters using MapReduce. In this paper, we describe an implementation of a particle filter using MapReduce. We focus on a component that what would otherwise be a bottleneck to parallel execution, the resampling component. We devise a new implementation of this component, which requires no approximations, has $O\left(N\right)$ spatial complexity and deterministic $O\left(\left(\log N\right)^2\right)$ time complexity. Results demonstrate the utility of this new component and culminate in consideration of a particle filter with $2^{24}$ particles being distributed across $512$ processor cores.


title: Learning Hierarchical Shape Segmentation and Labeling from Online Repositories

abstract: We propose a method for converting geometric shapes into hierarchically segmented parts with part labels. Our key idea is to train category-specific models from the scene graphs and part names that accompany 3D shapes in public repositories. These freely-available annotations represent an enormous, untapped source of information on geometry. However, because the models and corresponding scene graphs are created by a wide range of modelers with different levels of expertise, modeling tools, and objectives, these models have very inconsistent segmentations and hierarchies with sparse and noisy textual tags. Our method involves two analysis steps. First, we perform a joint optimization to simultaneously cluster and label parts in the database while also inferring a canonical tag dictionary and part hierarchy. We then use this labeled data to train a method for hierarchical segmentation and labeling of new 3D shapes. We demonstrate that our method can mine complex information, detecting hierarchies in man-made objects and their constituent parts, obtaining finer scale details than existing alternatives. We also show that, by performing domain transfer using a few supervised examples, our technique outperforms fully-supervised techniques that require hundreds of manually-labeled models.


title: Execution Templates: Caching Control Plane Decisions for Strong Scaling of Data Analytics

abstract: Control planes of cloud frameworks trade off between scheduling granularity and performance. Centralized systems schedule at task granularity, but only schedule a few thousand tasks per second. Distributed systems schedule hundreds of thousands of tasks per second but changing the schedule is costly. We present execution templates, a control plane abstraction that can schedule hundreds of thousands of tasks per second while supporting fine-grained, per-task scheduling decisions. Execution templates leverage a program's repetitive control flow to cache blocks of frequently-executed tasks. Executing a task in a template requires sending a single message. Large-scale scheduling changes install new templates, while small changes apply edits to existing templates. Evaluations of execution templates in Nimbus, a data analytics framework, find that they provide the fine-grained scheduling flexibility of centralized control planes while matching the strong scaling of distributed ones. Execution templates support complex, real-world applications, such as a fluid simulation with a triply nested loop and data dependent branches.


title: Supercongruences for rigid hypergeometric Calabi--Yau threefolds

abstract: We establish the supercongruences for the fourteen rigid hypergeometric Calabi--Yau threefolds over $\mathbb Q$ conjectured by Rodriguez-Villegas in 2003. Our first method is based on Dwork's theory of $p$-adic unit roots and it allows us to establish the supercongruences between the truncated hypergeometric series and the corresponding unit roots for ordinary primes. The other method makes use of the theory of hypergeometric motives, in particular, adapts the techniques from the recent work of Beukers, Cohen and Mellit on finite hypergeometric sums over $\mathbb Q$. Essential ingredients in executing the both approaches are the modularity of the underlying Calabi--Yau threefolds and a $p$-adic perturbation method applied to hypergeometric functions.


title: Picometer-scale atom position analysis in annular bright-field STEM imaging

abstract: We study the effects of specimen mistilt on the picometer-scale measurement of local structure by combing experiment and simulation in annular bright-field scanning transmission electron microscopy (ABF-STEM). A relative distance measurement method is proposed to separate the tilt effects from the scan noise and scan distortion. We find that under a typical experimental condition a small specimen tilt (~6 mrad) in 25 nm thick SrTiO3 along [001] causes 11.9 pm artificial displacement between O and Sr/TiO columns in ABF image, which is more than 3 times of scan noise and sample drift induced image distortion ~3.2 pm, suggesting the tilt effect could be dominant for the quantitative analysis of ABF images. The artifact depends the crystal mistilt angle, specimen thickness, defocus, convergence angle and uncorrected aberration. Our study provides useful insights into detecting and correcting tilt effects during both experiment operation and data analysis to extract the real structure information and avoid mis-interpretations of atomic structure as well as the properties such as oxygen octahedral distortion/shift.


title: GalRotpy: an educational tool to understand and parametrize the rotation curve and gravitational potential of disk-like galaxies

abstract: \textbf{GalRotpy} is an educational \verb+Python3+-based visual tool, which is useful to undestand how is the contribution of each mass component to the gravitational potential of disc-like galaxies by means of their rotation curve. Besides, \textbf{GalRotpy} allows the user to perform a parametric fit of a given rotation curve, which relies on a MCMC procedure implemented by using \verb+emcee+ package. Here the gravitational potential of disc-like galaxies is built from the contribution of a Miyamoto-Nagai potential model for the bulge/core and the thin/thick disc, an exponential disc, together with the NFW (Navarro-Frenk- White) potential or the Burkert (cored density profile) potential for the Dark Matter halo, where each contribution is implemented by using \verb+galpy+ package. We summarize the properties of each contribution to the rotation curve involved, and then describe how \textbf{GalRotpy} is implemented along with its capabilities. Finally we present the characterization of two galaxies, NGC6361 and M33, and show that the results for M33 provided by \textbf{GalRotpy} are consistent with those found in the literature.


title: EUV and Magnetic Activities Associated with Type-I Solar Radio Bursts

abstract: Type-I bursts (i.e. noise storms) are the earliest-known type of solar radio emission at the metre wavelength. They are believed to be excited by non-thermal energetic electrons accelerated in the corona. The underlying dynamic process and exact emission mechanism still remain unresolved. Here, with a combined analysis of extreme ultraviolet (EUV), radio and photospheric magnetic field data of unprecedented quality recorded during a type-I storm on 30 July 2011, we identify a good correlation between the radio bursts and the co-spatial EUV and magnetic activities. The EUV activities manifest themselves as three major brightening stripes above a region adjacent to a compact sunspot, while the magnetic field there presents multiple moving magnetic features (MMFs) with persistent coalescence or cancelation and a morphologically similar three-part distribution. We find that the type-I intensities are correlated with those of the EUV emissions at various wavelengths with a correlation coefficient of 0.7-0.8. In addition, in the region between the brightening EUV stripes and the radio sources there appear consistent dynamic motions with a series of bi-directional flows, suggesting ongoing small-scale reconnection there. Mainly based on the induced connection between the magnetic motion at the photosphere and the EUV and radio activities in the corona, we suggest that the observed type-I noise storms and the EUV brightening activities are the consequence of small-scale magnetic reconnection driven by MMFs. This is in support of the original proposal made by Bentely et al. (Solar Phys. 193, 227, 2000).


title: Link Mining for Kernel-based Compound-Protein Interaction Predictions Using a Chemogenomics Approach

abstract: Virtual screening (VS) is widely used during computational drug discovery to reduce costs. Chemogenomics-based virtual screening (CGBVS) can be used to predict new compound-protein interactions (CPIs) from known CPI network data using several methods, including machine learning and data mining. Although CGBVS facilitates highly efficient and accurate CPI prediction, it has poor performance for prediction of new compounds for which CPIs are unknown. The pairwise kernel method (PKM) is a state-of-the-art CGBVS method and shows high accuracy for prediction of new compounds. In this study, on the basis of link mining, we improved the PKM by combining link indicator kernel (LIK) and chemical similarity and evaluated the accuracy of these methods. The proposed method obtained an average area under the precision-recall curve (AUPR) value of 0.562, which was higher than that achieved by the conventional Gaussian interaction profile (GIP) method (0.425), and the calculation time was only increased by a few percent.


title: On the Necessity of Superparametric Geometry Representation for Discontinuous Galerkin Methods on Domains with Curved Boundaries

abstract: We provide numerical evidence demonstrating the necessity of employing a superparametric geometry representation in order to obtain optimal convergence orders on two-dimensional domains with curved boundaries when solving the Euler equations using Discontinuous Galerkin methods. However, concerning the obtention of optimal convergence orders for the Navier-Stokes equations, we demonstrate numerically that the use of isoparametric geometry representation is sufficient for the case considered here.


title: Compact Resolved Ejecta in the Nearest Tidal Disruption Event

abstract: Tidal disruption events (TDEs) occur when a star or sub-stellar object passes close enough to a galaxy's supermassive black hole to be disrupted by tidal forces. NGC 4845 (d=17 Mpc) was host to a TDE, IGR J12580+0134, detected in November 2010. Its proximity offers us a unique close-up of the TDE and its aftermath. We discuss new Very Long Baseline Array (VLBA) and Karl G. Jansky Very Large Array (JVLA) observations, which show that the radio flux from the active nucleus created by the TDE has decayed in a manner consistent with predictions from a jet-circumnuclear medium interaction model. This model explains the source's broadband spectral evolution, which shows a spectral peak that has moved from the submm (at the end of 2010) to GHz radio frequencies (in 2011-2013) to <1 GHz in 2015. The milliarcsecond-scale core is circularly polarized at 1.5 GHz but not at 5 GHz, consistent with the model. The VLBA images show a complex structure at 1.5 GHz that includes an east west extension ~40 milliarcsec (3 pc) long as well as a resolved component 52 milliarcsec (4.1 pc) northwest of the flat-spectrum core, which is all that can be seen at 5 GHz. If ejected in 2010, the NW component must have had v=0.96 c over five years. However, this is unlikely, as our model suggests strong deceleration to speeds < 0.5c within months and a much smaller, sub-parsec size. In this interpretation, the northwest component could have either a non-nuclear origin or be from an earlier event.


title: Qubit-loss-free fusion of W states employing weak cross-Kerr nonlinearities

abstract: With the assistance of weak cross-Kerr nonlinearities, we introduce an optical scheme to fuse two small-size polarization entangled W states into a large-scale W state without qubit loss, i.e.,$\mathrm{W}_{n+m}$ state can be generated from an $n$-qubit W state and a $m$-qubit W state. To complete the fusion task, two polarization entanglement processes and one spatial entanglement process are applied. The fulfillments of the above processes are contributed by a cross-Kerr nonlinear interaction between the signal photons and a coherent state via Kerr media. We analyze the resource cost and the success probability of the scheme. There is no complete failure output in our fusion mechanism, and all the garbage states are recyclable. In addition, there is no need for any controlled quantum gate and any ancillary photon, so it is simple and feasible under the current experiment technology.


title: Towards Simulation and Risk Assessment of Weather-Related Cascading Outages

abstract: Weather and environmental factors are verified to have played significant roles in historical major cascading outages and blackouts. Therefore, in the simulation and risk assessment of cascading outages in power systems, it is necessary to consider the weather and environmental effects. This paper proposes a method for the risk assessment of weather-related cascading outages. Based on the analysis of historical outage records and temperature-dependent physical outage mechanisms of transmission lines, an outage rate model considering weather condition and conductor temperature is proposed, and the analytical form of outage probability of lines are derived. With the weather-dependent outage model, a two-stage risk assessment method based on Markovian tree (MT) search is proposed, which consists of offline full assessment, and online efficient update of risk assessment results and continued MT search using updated NWP data. The test cases on NPCC 140-bus test system model in winter and summer scenarios verify the advantages of the proposed risk assessment method in both accuracy and efficiency.


title: Infinite time blow-up for the 3-dimensional energy critical heat equation

abstract: We construct globally defined in time, unbounded positive solutions to the energy-critical heat equation in dimension three $$ u_t = \Delta u + u^5 , \quad {\mbox {in}} \quad \R^3 \times (0,\infty), \ \ u(x, 0)= u_0 (x)\inn \R^3. $$ For each $\gamma>1$ we find initial data (not necessarily radially symmetric) with $\lim\limits_{r \to \infty} |x|^\gamma u_0 (x) >0$ such that as $t \to \infty$ $$ \| u(\cdot ,t ) \|_\infty \sim t^{\gamma-1 \over 2} , \quad {\mbox {if}} \quad 1<\gamma <2, \quad \| u(\cdot ,t ) \|_\infty \sim \sqrt{t}, \quad {\mbox {if}} \quad \gamma >2, \quad $$ and $$ \| u(\cdot , t)\|_\infty \sim \sqrt{t}\, (\ln t )^{-1} , \quad {\mbox {if}} \quad \gamma = 2. $$ Furthermore we show that this infinite time blow-up is co-dimensional one stable. The existence of such solutions was conjectured by Fila and King.


title: Ab initio study of sodium cointercalation with diglyme molecule into graphite

abstract: The cointercalation of sodium with the solvent organic molecule into graphite can resolve difficulty of forming the stage-I Na-graphite intercalation compound, which is a predominant anode of Na-ion battery. To clarify the mechanism of such cointercalation, we investigate the atomistic structure, energetics, electrochemical properties, ion and electron conductance, and charge transferring upon de/intercalation of the solvated Na-diglyme ion into graphite with {\it ab initio} calculations. It is found that the Na(digl)$_2$C$_n$ compound has the negatively lowest intercalation energy at $n\approx$21, the solvated Na(digl)$_2$ ion diffuses fast in the interlayer space, and their electronic conductance can be enhanced compared to graphite. The calculations reveal that the diglyme molecules as well as Na atom donates electrons to the graphene layer, resulting in the formation of ionic bonding between the graphene layer and the moiety of diglyme molecule. This work will contribute to the development of innovative anode materials for alkali-ion battery applications.


title: Semi-Global Weighted Least Squares in Image Filtering

abstract: Solving the global method of Weighted Least Squares (WLS) model in image filtering is both time- and memory-consuming. In this paper, we present an alternative approximation in a time- and memory- efficient manner which is denoted as Semi-Global Weighed Least Squares (SG-WLS). Instead of solving a large linear system, we propose to iteratively solve a sequence of subsystems which are one-dimensional WLS models. Although each subsystem is one-dimensional, it can take two-dimensional neighborhood information into account due to the proposed special neighborhood construction. We show such a desirable property makes our SG-WLS achieve close performance to the original two-dimensional WLS model but with much less time and memory cost. While previous related methods mainly focus on the 4-connected/8-connected neighborhood system, our SG-WLS can handle a more general and larger neighborhood system thanks to the proposed fast solution. We show such a generalization can achieve better performance than the 4-connected/8-connected neighborhood system in some applications. Our SG-WLS is $\sim20$ times faster than the WLS model. For an image of $M\times N$, the memory cost of SG-WLS is at most at the magnitude of $max\{\frac{1}{M}, \frac{1}{N}\}$ of that of the WLS model. We show the effectiveness and efficiency of our SG-WLS in a range of applications. The code is publicly available at: https://github.com/wliusjtu/Semi-Global-Weighted-Least-Squares-in-Image-Filtering.


title: Pricing in non-convex markets with quadratic deliverability costs

abstract: The problem of obtaining market-clearing prices for markets with non-convexities has been widely studied in the literature. This is particularly the case in electricity markets, where worldwide deregulation leads to markets in which non-convexities arise from the decisions of market operators regarding which generators are committed to provide electricity power. Here, we extend seminal results in this area to address the problem of obtaining market-clearing prices for markets in which beyond non-convexities, it is relevant to account for convex quadratic market costs. In a general market, such costs arise from quadratic commodity costs or transactions costs. In an electricity market, such quadratic costs arise when ramping costs need to be considered due to the presence of renewable energy sources, which continue to increase their participation in electricity markets. To illustrate our results, we compute and analyze the clearing prices of a classical market problem with the addition of ramping costs.


title: Destroying a near-extremal Kerr-Newman-AdS black hole

abstract: If radiative and self-force effects are neglected, we find that feeding a test particle into a near-extremal Kerr-Newman-AdS black hole could lead to destroy their event horizon, giving rise to naked singularities. Hence radiative and self-force effects must be taken into account to further test cosmic censorship. Moreover, the allowed parameter range for this test particle is very narrow, this leaves the possibility of radiative and self-force effects considerations to cure the problem of WCC violation in Kerr-Newman-AdS spacetime.


title: Optimized Regression Discontinuity Designs

abstract: The increasing popularity of regression discontinuity methods for causal inference in observational studies has led to a proliferation of different estimating strategies, most of which involve first fitting non-parametric regression models on both sides of a treatment assignment boundary and then reporting plug-in estimates for the effect of interest. In applications, however, it is often difficult to tune the non-parametric regressions in a way that is well calibrated for the specific target of inference; for example, the model with the best global in-sample fit may provide poor estimates of the discontinuity parameter. In this paper, we propose an alternative method for estimation and statistical inference in regression discontinuity designs that uses numerical convex optimization to directly obtain the finite-sample-minimax linear estimator for the regression discontinuity parameter, subject to bounds on the second derivative of the conditional response function. Given a bound on the second derivative, our proposed method is fully data-driven, and provides uniform confidence intervals for the regression discontinuity parameter with both discrete and continuous running variables. The method also naturally extends to the case of multiple running variables.