title: Simulation of Systematics in Future Single-Dish HI Intensity Mapping Experiments

abstract: HI intensity mapping (IM) is an exciting new probe that could revolutionize the future of cosmology. However, the relative faintness of the HI signal when compared to foregrounds of astrophysical or terrestrial origin will make HI IM extremely challenging. The imprint of these foregrounds may result in systematic errors in the recovered cosmological signal. We discuss an IM simulation pipeline developed at Manchester that can introduce systematic errors at the TOD level in order to help assess their impact. We will present results for two potential sources of systematics for HI IM surveys: 1/f noise and the integrated emission from global navigation satellites.


title: Deleting edges to restrict the size of an epidemic in temporal networks

abstract: Spreading processes on graphs are a natural model for a wide variety of real-world phenomena, including information spread over social networks and biological diseases spreading over contact networks. Often, the networks over which these processes spread are dynamic in nature, and can be modeled with temporal graphs. Here, we study the problem of deleting edges from a given temporal graph in order to reduce the number of vertices (temporally) reachable from a given starting point. This could be used to control the spread of a disease, rumour, etc. in a temporal graph. In particular, our aim is to find a temporal subgraph in which a process starting at any single vertex can be transferred to only a limited number of other vertices using a temporally-feasible path. We introduce a natural edge-deletion problem for temporal graphs and provide positive and negative results on its computational complexity and approximability.


title: A fast algorithm with minimax optimal guarantees for topic models with an unknown number of topics

abstract: We propose a new method of estimation in topic models, that is not a variation on the existing simplex finding algorithms, and that estimates the number of topics K from the observed data. We derive new finite sample minimax lower bounds for the estimation of A, as well as new upper bounds for our proposed estimator. We describe the scenarios where our estimator is minimax adaptive. Our finite sample analysis is valid for any number of documents (n), individual document length (N_i), dictionary size (p) and number of topics (K), and both p and K are allowed to increase with n, a situation not handled well by previous analyses. We complement our theoretical results with a detailed simulation study. We illustrate that the new algorithm is faster and more accurate than the current ones, although we start out with a computational and theoretical disadvantage of not knowing the correct number of topics K, while we provide the competing methods with the correct value in our simulations.


title: Sums of Weighted Differentiation Composition Operators

abstract: We solve an interpolation problem in $A^p_\alpha$ involving specifying a set of (possibly not distinct) $n$ points, where the $k^{\textrm{th}}$ derivative at the $k^{\textrm{th}}$ point is up to a constant as large as possible for functions of unit norm. The solution obtained has norm bounded by a constant independent of the points chosen. As a direct application, we obtain a characterization of the order-boundedness of a sum of products of weighted composition and differentiation operators acting between weighted Bergman spaces. We also characterize the compactness of such operators that map a weighted Bergman space into the space of bounded analytic functions.


title: Methods for the inclusion of real world evidence in network meta-analysis

abstract: Background: Network Meta-Analysis (NMA) is a key component of submissions to reimbursement agencies world-wide, especially when there is limited direct head-to-head evidence for multiple technologies from randomised controlled trials (RCTs). Many NMAs include only data from RCTs. However, real-world evidence (RWE) is also becoming widely recognised as a valuable source of clinical data. We investigate methods for the inclusion of RWE in NMA and its impact on the uncertainty around the effectiveness estimates. Methods: A range of methods for inclusion of RWE in evidence synthesis, including Bayesian hierarchical and power prior models, were investigated by applying them to an example in relapsing remitting multiple sclerosis. The effect of the inclusion of RWE was investigated by varying the degree of down weighting of this part of evidence by the use of a power prior. Results: Whilst the inclusion of the RWE led to an increase in the level of uncertainty surrounding effect estimates in this example, this depended on the method of inclusion adopted for the RWE. Power prior NMA model resulted in stable effect estimates for fingolimod yet increasing the width of the credible intervals with increasing weight given to RWE data. The hierarchical NMA models were effective in allowing for heterogeneity between study designs; however, this also increased the level of uncertainty. Conclusion: The power prior approach for the inclusion of RWE in NMAs indicates that the degree to which RWE is taken into account can have a significant impact on the overall level of uncertainty. The hierarchical modelling approach further allowed for accommodating differences between study types. Consequently, further work investigating both empirical evidence for biases associated with individual RWE studies and methods of elicitation from experts on the extent of such biases is warranted.


title: Extending functions from a neighborhood of the sphere to the ball

abstract: In this article, we are interested in the problem of extending the germ of a smooth function $\tilde{f}$ defined along the standard sphere of dimension $n$ to a function defined on the ball which has no critical points. The article gives a necessary condition using the Morse chain complex associated to the function $f$, restriction of $\tilde{f}$ to the sphere $\mathbb{S}^n$, which is assumed to be a Morse function.


title: Is the IMF in ellipticals bottom-heavy? Clues from their chemical abundances

abstract: We tested the implementation of different IMFs in our model for the chemical evolution of ellipticals, with the aim of reproducing the observed relations of [Fe/H] and [Mg/Fe] abundances with galaxy mass in a sample of early-type galaxies selected from the SPIDER-SDSS catalog. Abundances in the catalog were derived from averaged spectra, obtained by stacking individual spectra according to central velocity dispersion, as a proxy of galaxy mass. We tested initial mass functions already used in a previous work, as well as two new models, based on low-mass tapered ("bimodal") IMFs, where the IMF becomes either (1) bottom-heavy in more massive galaxies, or (2) is time-dependent, switching from top-heavy to bottom-heavy in the course of galactic evolution. We found that observations could only be reproduced by models assuming either a constant, Salpeter IMF, or a time-dependent distribution, as other IMFs failed. We further tested the models by calculating their M/L ratios. We conclude that a constant, time-independent bottom-heavy IMF does not reproduce the data, especially the increase of the $[\alpha/Fe]$ ratio with galactic stellar mass, whereas a variable IMF, switching from top to bottom-heavy, can match observations. For the latter models, the IMF switch always occurs at the earliest possible considered time, i.e. $t_{\text{switch}}= 0.1$ Gyr.


title: Weakly singular shock profiles for a non-dispersive regularization of shallow-water equations

abstract: We study a regularization of the classical Saint-Venant (shallow-water) equations, recently introduced by D. Clamond and D. Dutykh (Commun. Nonl. Sci. Numer. Simulat. 55 (2018) 237-247). This regularization is non-dispersive and formally conserves mass, momentum and energy. We show that for every classical shock wave, the system admits a corresponding non-oscillatory traveling wave solution which is continuous and piecewise smooth, having a weak singularity at a single point where energy is dissipated as it is for the classical shock. The system also admits cusped solitary waves of both elevation and depression.


title: Decay behavior of the strange and beauty partners of $P_c$ hadronic molecules

abstract: We extend our previous study on the decay behavior of $P_c$ hadronic molecules to their strange and beauty partners. While $P_c(4380)$ and $P_c(4450)$ locate just below the $\bar D\Sigma_c^*$ and $\bar D^*\Sigma_c$ thresholds, respectively, their proposed strange partners $N(1875)$ and $N(2080)$ sit just below the $K\Sigma^*$ and $K^*\Sigma$ thresholds, respectively. Using the effective Lagrangian approach as the same as for the study of $P_c$ hadronic molecular states and with the couplings determined by the $\mathrm{SU}(3)$ flavor symmetry, the decay patterns of $N(1875)$ and $N(2080)$ are obtained assuming them to be the $S$-wave $K\Sigma^*$ and $K^*\Sigma$ molecules, respectively. It is found that the measured decay properties of $N(1875)$ and $N(2080)$ are reproduced well. Our results suggest that both $N(1875)$ and $N(2080)$ can be ascribed as the hadronic molecular pentaquark states in the hidden strangeness sector with quantum number $I(J^P)=1/2({3/2}^-)$. Further checks on their hadronic molecule nature by various experiments are proposed. With the same approach, we also give our predictions for the decay behavior of the possible beauty partners of the $P_c$ hadronic molecular states. The $B^*\Lambda_b$ is found to be the largest decay mode for both $B\Sigma_b^*$ and $B^*\Sigma_b$ hadronic molecules, hence can be used to look for the hidden beauty pentaquark states in forthcoming experiments.


title: Time Fractional Schro\"odinger Equation

abstract: We propose a time fractional extension of the Schr{\"o}dinger equation that keeps the main mechanical and quantum properties of the classical Schr{\"o}dinger equation. This extension is shown to be equivalent to another well identified time first order PDE with fractional hamiltonian.


title: Two-dimensional nonlinear modes and frequency combs in bottle microresonators

abstract: We investigate theoretically frequency comb generation in a bottle microresonator accounting for the azimuthal and axial degrees of freedom. We first identify a discrete set of the axial nonlinear modes of a bottle microresonator that appear as tilted resonances bifurcating from the spectrum of linear axial modes. We then study azimuthal modulational instability of these modes and show that families of 2D soliton states localized both azimuthally and axially bifurcate from them at critical pump frequencies. Depending on detuning, 2D solitons can be either stable, or form persistent breathers, chaotic spatio-temporal patterns, or exhibit collapse-like evolution.


title: RotDCF: Decomposition of Convolutional Filters for Rotation-Equivariant Deep Networks

abstract: Explicit encoding of group actions in deep features makes it possible for convolutional neural networks (CNNs) to handle global deformations of images, which is critical to success in many vision tasks. This paper proposes to decompose the convolutional filters over joint steerable bases across the space and the group geometry simultaneously, namely a rotation-equivariant CNN with decomposed convolutional filters (RotDCF). This decomposition facilitates computing the joint convolution, which is proved to be necessary for the group equivariance. It significantly reduces the model size and computational complexity while preserving performance, and truncation of the bases expansion serves implicitly to regularize the filters. On datasets involving in-plane and out-of-plane object rotations, RotDCF deep features demonstrate greater robustness and interpretability than regular CNNs. The stability of the equivariant representation to input variations is also proved theoretically under generic assumptions on the filters in the decomposed form. The RotDCF framework can be extended to groups other than rotations, providing a general approach which achieves both group equivariance and representation stability at a reduced model size.


title: Quantitative structure of stable sets in finite abelian groups

abstract: We prove an arithmetic regularity lemma for stable subsets of finite abelian groups, generalising our previous result for high-dimensional vector spaces over finite fields of prime order. A qualitative version of this generalisation was recently obtained by the first author in joint work with Conant and Pillay, using model-theoretic techniques. In contrast, the approach in the present paper is highly quantitative and relies on several key ingredients from arithmetic combinatorics.


title: Edge-statistics on large graphs

abstract: The inducibility of a graph $H$ measures the maximum number of induced copies of $H$ a large graph $G$ can have. Generalizing this notion, we study how many induced subgraphs of fixed order $k$ and size $\ell$ a large graph $G$ on $n$ vertices can have. Clearly, this number is $\binom{n}{k}$ for every $n$, $k$ and $\ell \in \left \{0, \binom{k}{2} \right\}$. We conjecture that for every $n$, $k$ and $0 < \ell < \binom{k}{2}$ this number is at most $\left(1/e + o_k(1) \right) \binom{n}{k}$. If true, this would be tight for $\ell \in \{1, k-1\}$. In support of our `Edge-statistics conjecture' we prove that the corresponding density is bounded away from $1$ by an absolute constant. Furthermore, for various ranges of the values of $\ell$ we establish stronger bounds. In particular, we prove that for `almost all' pairs $(k, \ell)$ only a polynomially small fraction of the $k$-subsets of $V(G)$ has exactly $\ell$ edges, and prove an upper bound of $(1/2 + o_k(1))\binom{n}{k}$ for $\ell = 1$. Our proof methods involve probabilistic tools, such as anti-concentration results relying on fourth moment estimates and Brun's sieve, as well as graph-theoretic and combinatorial arguments such as Zykov's symmetrization, Sperner's theorem and various counting techniques.


title: A More Accurate and Competitive Estimative of $H_0$ in Intermediate Redshifts

abstract: In order to clarify the tension between estimates of the Hubble Constant ($H_0$) from local ($z \ll 1$) and global ($z \gg 1$) measurements, Lima and Cunha (LC) proposed a new method to measure $H_0$ in intermediate redshifts ($z \approx 1$), which were obtained $H_0 = 74.1 \pm 2.2$ km s$^{-1} $Mpc$^{-1}$ ($1 \sigma$), in full agreement to local measurements via Supernovae/Cepheid dataset. However, Holanda \textit{et al.} (2014) affirm that a better understanding of the morphology of galaxy clusters in LC framework is needed to a more robust and accurate determination of $H_0$. Moreover, that kind of sample has been strongly questioned in the literature. In this context, (i) we investigated if the sample of galaxy clusters used by LC has a relevant role in their results, then (ii) we perform a more accurate and competitive determination of $H_0$ in intermediate redshifts, free of unknown systematic uncertainties. First, we found that the exclusion of the sample of galaxy clusters from the determination initially proposed by LC leads to significantly different results. Finally, we performed a new determination in $H_0$, where we obtained $H_0 = 68.00 \pm 2.20$ km s$^{-1}$ Mpc$^{-1}$ ($1 \sigma$) with statistical and systematic errors and $H_0 = 68.71^{+1.37}_{-1.45}$ km s$^{-1}$ Mpc$^{-1}$ ($1 \sigma$) with statistical errors only. Contrary to those obtained by LC, these values are in full harmony with the global measurements via Cosmic Microwave Background (CMB) radiation and to the other recent estimates of $H_0$ in intermediate redshifts.


title: A NICER look at the Aql X-1 hard state

abstract: We report on a spectral-timing analysis of the neutron star low-mass X-ray binary Aql~X-1 with the Neutron Star Interior Composition Explorer (NICER) on the International Space Station. Aql~X-1 was observed with NICER during a dim outburst in 2017 July, collecting approximately $50$ ks of good exposure. The spectral and timing properties of the source correspond to that of an (hard) extreme island state in the atoll classification. We find that the fractional amplitude of the low frequency ($<0.3$ Hz) band-limited noise shows a dramatic turnover as a function of energy: it peaks at 0.5 keV with nearly 5\% rms, drops to $12\%$ rms at 2 keV, and rises to $15\%$ rms at 10 keV. Through the analysis of covariance spectra, we demonstrate that band-limited noise exists in both the soft thermal emission and the power-law emission. Additionally, we measure hard time lags, indicating the thermal emission at $0.5$ keV leads the power-law emission at 10 keV on a timescale of $\sim100$ ms at $0.3$ Hz to $\sim10$ ms at $3$ Hz. Our results demonstrate that the thermal emission in the hard state is intrinsically variable, and driving the modulation of the higher energy power-law. Interpreting the thermal spectrum as disk emission, we find our results are consistent with the disk propagation model proposed for accretion onto black holes.


title: Do Bi-Stable Poisson-Nernst-Planck Models Describe Single Channel Gating?

abstract: Experiments measuring currents through single protein channels show unstable currents, a phenomena called the gating of a single channel. Channels switch between an 'open' state with a well defined single amplitude of current and 'closed' states with nearly zero current. The existing mean-field theory of ion channels focuses almost solely on the open state. The physical modeling of the dynamical features of ion channels is still in its infancy, and does not describe the transitions between open and closed states, nor the distribution of the duration times of open states. One hypothesis is that gating corresponds to noise-induced fast transitions between multiple steady (equilibrium) states of the underlying system. In this work, we aim to test this hypothesis. Particularly, our study focuses on the (high order) steric Poisson-Nernst-Planck-Cahn-Hilliard model since it has been successful in predicting permeability and selectivity of ionic channels in their open state, and since it gives rise to multiple steady states. We show that this system gives rise to a gating-like behavior, but that important features of this switching behavior are different from the defining features of gating in biological systems. Furthermore, we show that noise prohibits switching in the system of study. The above phenomena are expected to occur in other PNP-type models, strongly suggesting that one has to go beyond over-damped (gradient flow) Nernst-Planck type dynamics to explain the spontaneous gating of single channels.


title: Saddle Point Least Squares Preconditioning of Mixed Methods

abstract: We present a simple way to discretize and precondition mixed variational formulations. Our theory connects with, and takes advantage of, the classical theory of symmetric saddle point problems and the theory of preconditioning symmetric positive definite operators. Efficient iterative processes for solving the discrete mixed formulations are proposed and choices for discrete spaces that are always compatible are provided. For the proposed discrete spaces and solvers, a basis is needed only for the test spaces and assembly of a global saddle point system is avoided. We prove sharp approximation properties for the discretization and iteration errors and also provide a sharp estimate for the convergence rate of the proposed algorithm in terms of the condition number of the elliptic preconditioner and the discrete $\inf-\sup$ and $\sup-\sup$ constants of the pair of discrete spaces.


title: Unusual interplay between superconductivity and field-induced charge order in YBa2Cu3Oy

abstract: We present a detailed study of the temperature (T) and magnetic field (H) dependence of the electronic density of states (DOS) at the Fermi level, as deduced from specific heat and Knight shift measurements in underdoped YBa2Cu3Oy. We find that the DOS becomes field-independent above a characteristic field H_{DOS} and that the H_{DOS}(T) line displays an unusual inflection near the onset of the long range 3D charge-density wave order. The unusual S-shape of H_{DOS}(T) is suggestive of two mutually-exclusive orders that eventually establish a form of cooperation in order to coexist at low T. On theoretical grounds, such a collaboration could result from the stabilisation of a pair-density wave state, which calls for further investigations in this region of the phase diagram


title: Free $(\mathbb{Z}/p)^n$-complexes and $p$-DG modules

abstract: We reformulate the problem of bounding the total rank of the homology of perfect chain complexes over the group ring $\mathbb{F}_p[G]$ of an elementary abelian $p$-group $G$ in terms of commutative algebra. This extends results of Carlsson for $p=2$ to all primes. As an intermediate step, we construct an embedding of the derived category of perfect chain complexes over $\mathbb{F}_p[G]$ into the derived category of $p$-DG modules over a polynomial ring.


title: Learning non-smooth models: instrumental variable quantile regressions and related problems

abstract: This paper proposes computationally efficient methods that can be used for instrumental variable quantile regressions (IVQR) and related methods with statistical guarantees. This is much needed when we investigate heterogenous treatment effects since interactions between the endogenous treatment and control variables lead to an increased number of endogenous covariates. We prove that the GMM formulation of IVQR is NP-hard and finding an approximate solution is also NP-hard. Hence, solving the problem from a purely computational perspective seems unlikely. Instead, we aim to obtain an estimate that has good statistical properties and is not necessarily the global solution of any optimization problem. The proposal consists of employing $k$-step correction on an initial estimate. The initial estimate exploits the latest advances in mixed integer linear programming and can be computed within seconds. One theoretical contribution is that such initial estimators and Jacobian of the moment condition used in the k-step correction need not be even consistent and merely $k=4\log n$ fast iterations are needed to obtain an efficient estimator. The overall proposal scales well to handle extremely large sample sizes because lack of consistency requirement allows one to use a very small subsample to obtain the initial estimate and the k-step iterations on the full sample can be implemented efficiently. Another contribution that is of independent interest is to propose a tuning-free estimation for the Jacobian matrix, whose definition nvolves conditional densities. This Jacobian estimator generalizes bootstrap quantile standard errors and can be efficiently computed via closed-end solutions. We evaluate the performance of the proposal in simulations and an empirical example on the heterogeneous treatment effect of Job Training Partnership Act.


title: Projections with fixed difference: a Hopf-Rinow theorem

abstract: The set $D_{A_0}$, of pairs of orthogonal projections $(P,Q)$ in generic position with fixed difference $P-Q=A_0$, is shown to be a homogeneus smooth manifold: it is the quotient of the unitary group of the commutant $\{A_0\}'$ divided by the unitary subgroup of the commutant $\{P_0, Q_0\}'$, where $(P_0,Q_0)$ is any fixed pair in $D_{A_0}$. Endowed with a natural reductive structure (a linear connection) and the quotient Finsler metric of the operator norm, it behaves as a classic Riemannian space: any two pairs in $D_{A_0}$ are joined by a geodesic of minimal length. Given a base pair $(P_0,Q_0)$, pairs in an open dense subset of $D_{A_0}$ can be joined to $(P_0,Q_0)$ by a {\it unique} minimal geodesic.


title: Topological Dirac Semimetal Na3Bi Films in the Ultrathin Limit via Alternating Layer Molecular Beam Epitaxy

abstract: Ultrathin films of Na3Bi on insulating substrates are desired for opening a bulk band gap and generating the quantum spin Hall effect from a topological Dirac semimetal, though continuous films in the few nanometer regime have been difficult to realize. Here, we utilize alternating layer molecular beam epitaxy (MBE) to achieve uniform and continuous single crystal films of Na3Bi(0001) on insulating Al2O3(0001) substrates and demonstrate electrical transport on films with 3.8 nm thickness (4 unit cells). The high material quality is confirmed through in situ reflection high-energy electron diffraction (RHEED), scanning tunneling microscopy (STM), x-ray diffraction (XRD), and x-ray photoelectron spectroscopy (XPS). In addition, these films are employed as seed layers for subsequent growth by codeposition, leading to atomic layer-by-layer growth as indicated by RHEED intensity oscillations. These material advances facilitate the pursuit of quantum phenomena in thin films of Dirac semimetals.


title: Phonon Quantum Nondemolition Measurements in Nonlinearly Coupled Optomechanical Cavities

abstract: In the field of cavity optomechanics, proposals for quantum nondemolition (QND) measurements of phonon number provide a promising avenue by which one can study the quantum nature of nanoscale mechanical resonators. Here, we investigate these QND measurements for an optomechanical system whereby quadratic coupling arises due to shared symmetries between a single optical resonance and a mechanical mode. We establish a relaxed limit on the amount of linear coupling that can exist in this type of system while still allowing for a QND measurement of Fock states. This new condition enables optomechanical QND measurements, which can be used to probe the decoherence of mesoscopic mechanical Fock states, providing an experimental testbed for quantum collapse theories.


title: On the total volume of the double hyperbolic space

abstract: Let the \emph{double hyperbolic space} $\mathbb{DH}^n$, proposed in this paper as an extension of the hyperbolic space $\mathbb{H}^n$, contain a two-sheeted hyperboloid with the two sheets connected to each other along the boundary at infinity. We propose to extend the volume of convex polytopes in $\mathbb{H}^n$ to polytopes in $\mathbb{DH}^n$, where the volume is invariant under isometry but can possibly be complex valued. We show that the total volume of $\mathbb{DH}^n$ is equal to $i^n V_n(\mathbb{S}^n)$ for both even and odd dimensions, and prove a Schl\"{a}fli differential formula (\SDF{}) for $\mathbb{DH}^n$. For $n$ odd, the volume of a polytope in $\mathbb{DH}^n$ is shown to be completely determined by its intersection with $\partial\mathbb{H}^n$ and induces a new intrinsic \emph{volume} on $\partial\mathbb{H}^n$ that is invariant under M\"{o}bius transformations.


title: Quantum Transport in a Low-Density Periodic Potential: Homogenisation via Homogeneous Flows

abstract: We show that the time evolution of a quantum wavepacket in a periodic potential converges in a combined high-frequency/Boltzmann-Grad limit, up to second order in the coupling constant, to terms that are compatible with the linear Boltzmann equation. This complements results of Eng and Erd\"os for low-density random potentials, where convergence to the linear Boltzmann equation is proved in all orders. We conjecture, however, that the linear Boltzmann equation fails in the periodic setting for terms of order four and higher. Our proof uses Floquet-Bloch theory, multi-variable theta series and equidistribution theorems for homogeneous flows. Compared with other scaling limits traditionally considered in homogenisation theory, the Boltzmann-Grad limit requires control of the quantum dynamics for longer times, which are inversely proportional to the total scattering cross section of the single-site potential.


title: Answer Set Programming Modulo `Space-Time'

abstract: We present ASP Modulo `Space-Time', a declarative representational and computational framework to perform commonsense reasoning about regions with both spatial and temporal components. Supported are capabilities for mixed qualitative-quantitative reasoning, consistency checking, and inferring compositions of space-time relations; these capabilities combine and synergise for applications in a range of AI application areas where the processing and interpretation of spatio-temporal data is crucial. The framework and resulting system is the only general KR-based method for declaratively reasoning about the dynamics of `space-time' regions as first-class objects. We present an empirical evaluation (with scalability and robustness results), and include diverse application examples involving interpretation and control tasks.


title: Design Identification of Curve Patterns on Cultural Heritage Objects: Combining Template Matching and CNN-based Re-Ranking

abstract: The surfaces of many cultural heritage objects were embellished with various patterns, especially curve patterns. In practice, most of the unearthed cultural heritage objects are highly fragmented, e.g., sherds of potteries or vessels, and each of them only shows a very small portion of the underlying full design, with noise and deformations. The goal of this paper is to address the challenging problem of automatically identifying the underlying full design of curve patterns from such a sherd. Specifically, we formulate this problem as template matching: curve structure segmented from the sherd is matched to each location with each possible orientation of each known full design. In this paper, we propose a new two-stage matching algorithm, with a different matching cost in each stage. In Stage 1, we use a traditional template matching, which is highly computationally efficient, over the whole search space and identify a small set of candidate matchings. In Stage 2, we derive a new matching cost by training a dual-source Convolutional Neural Network (CNN) and apply it to re-rank the candidate matchings identified in Stage 1. We collect 600 pottery sherds with 98 full designs from the Woodland Period in Southeastern North America for experiments and the performance of the proposed algorithm is very competitive.


title: A remark on constant scalar curvature K\"ahler metrics on minimal models

abstract: In this short note, we prove the existence of constant scalar curvature K\"ahler metrics on compact K\"ahler manifolds with semi-ample canonical bundles.


title: Resource allocation under uncertainty: an algebraic and qualitative treatment

abstract: We use an algebraic viewpoint, namely a matrix framework to deal with the problem of resource allocation under uncertainty in the context of a qualitative approach. Our basic qualitative data are a plausibility relation over the resources, a hierarchical relation over the agents and of course the preference that the agents have over the resources. With this data we propose a qualitative binary relation $\unrhd$ between allocations such that $\mathcal{F}\unrhd \mathcal{G}$ has the following intended meaning: the allocation $\mathcal{F}$ produces more or equal social welfare than the allocation $\mathcal{G}$. We prove that there is a family of allocations which are maximal with respect to $\unrhd$. We prove also that there is a notion of simple deal such that optimal allocations can be reached by sequences of simple deals. Finally, we introduce some mechanism for discriminating {optimal} allocations.


title: Optimal Scheduling and Exact Response Time Analysis for Multistage Jobs

abstract: Scheduling to minimize mean response time in an M/G/1 queue is a classic problem. The problem is usually addressed in one of two scenarios. In the perfect-information scenario, the scheduler knows each job's exact size, or service requirement. In the zero-information scenario, the scheduler knows only each job's size distribution. The well-known shortest remaining processing time (SRPT) policy is optimal in the perfect-information scenario, and the more complex Gittins policy is optimal in the zero-information scenario. In real systems the scheduler often has partial but incomplete information about each job's size. We introduce a new job model, that of multistage jobs, to capture this partial-information scenario. A multistage job consists of a sequence of stages, where both the sequence of stages and stage sizes are unknown, but the scheduler always knows which stage of a job is in progress. We give an optimal algorithm for scheduling multistage jobs in an M/G/1 queue and an exact response time analysis of our algorithm.


title: On the projections of mutual multifractal spectra

abstract: The aim of this article is to study the behaviour of the relative multifractal spectrum under projections. First of all, we depict a relationship between the mutual multifractal spectra of a couple of measures $(\mu, \nu)$ and its orthogonal projections in Euclidean space. As an application, we improve Svetova's result (Tr. Petrozavodsk. Gos. Univ. Ser. Mat., 11 (2004), 41-46) and study the mutual multifractal analysis of the projections of measures.


title: Superconductivity from Valley Fluctuations and Approximate SO(4) Symmetry in a Weak Coupling Theory of Twisted Bilayer Graphene

abstract: We develop a weak coupling approach to superconductivity in twisted bilayer graphene, starting from the Fermi liquid regime. A key observation is that near half filling, the fermiology consists of well nested Fermi pockets derived from opposite valleys, leading to enhanced valley fluctuation, which in turn can mediate superconductivity. This scenario is studied within the random phase approximation. We find that inter-valley electron pairing with either chiral ($d+i d$ mixed with $p-i p$) or helical form factor is the dominant instability. An approximate SO(4) spin-valley symmetry implies a near degeneracy of spin-singlet and triplet pairing. On increasing interactions, commensurate inter-valley coherence wave (IVCW) order can arise, with simultaneous condensation at the three M points in the Brillouin Zone, and a $2\times2$ pattern in real space. In simple treatments though, this leads to a full gap at fillings $\pm (1/2+1/8)$, slightly away from half-filling. An SO(4) symmetry breaking "anti-Hunds" coupling favors the spin-singlet order both for the IVCW and the superconductor, consistent with observations. Mott insulators derived from phase fluctuating superconductors are also discussed, which exhibit both symmetry protected and intrinsic topological orders.


title: Gaussian functions are optimal for waveguided nonlinear-quantum-optical processes

abstract: Many nonlinear optical technologies require the two-mode spectral amplitude function that describes them---the \emph{joint spectral amplitude} (JSA)---to be separable. We prove that the JSA factorizes \emph{only} when the incident pump field and phase-matching function are Gaussian functions. We show this by mapping our problem to a known result, in continuous variable quantum information, that only squeezed states remain unentangled when combined on a beam splitter. We then conjecture that only a squeezed state minimizes entanglement when sent through a beam splitter with another pre-specified ket. This implies that to maximize JSA separability when one of the (pump or nonlinear medium) functions is non-Gaussian, the other function \emph{must} be Gaussian. This answers an outstanding question about optimal design of certain nonlinear processes, and is of practical interest to researchers using waveguide nonlinear optics to generate and manipulate quantum light.


title: Revisiting the tree edit distance and its backtracing: A tutorial

abstract: Almost 30 years ago, Zhang and Shasha (1989) published a seminal paper describing an efficient dynamic programming algorithm computing the tree edit distance, that is, the minimum number of node deletions, insertions, and replacements that are necessary to transform one tree into another. Since then, the tree edit distance has been widely applied, for example in biology and intelligent tutoring systems. However, the original paper of Zhang and Shasha can be challenging to read for newcomers and it does not describe how to efficiently infer the optimal edit script. In this contribution, we provide a comprehensive tutorial to the tree edit distance algorithm of Zhang and Shasha. We further prove metric properties of the tree edit distance, and describe efficient algorithms to infer the cheapest edit script, as well as a summary of all cheapest edit scripts between two trees.


title: Quantum-enabled electrometer measures field transients and correlation function

abstract: We realize a non-invasive electrometer based on state engineering in a Rydberg hydrogenic manifold. A quantum interference process involving states with very different dipoles measures directly the time correlation of a stochastic electric field, with a 200~ns time resolution and with a $83.2\pm1.4$~mV/m single-atom sensitivity, beyond the Standard Quantum Limit in this context. This quantum-enabled correlation electrometer opens the way to applications in mesoscopic physics, e.g. for the study of individual charge transits in nano-structures.


title: Visiting Newton's Atelier before the Principia, 1679-1684

abstract: The manuscripts that presumably contained Newton's early development of the fundamental concepts that led to his Principia have been lost. A plausible reconstruction of this development is presened based on Newton's exchange of letters with Robert Hooke in 1679, with Edmund Halley in 1686, and on some clues in the diagram associated with Proposition1 in Book1 of the Principia that have been ignored in the past. The graphical method associated with this proposition leads to a rapidly convergent method to obtain orbital curves for central forces, and elucidates how Newton may have have been led to formulate some of his other propositions in the Principia.


title: Coding for Interactive Communication with Small Memory and Applications to Robust Circuits

abstract: Classically, coding theory has been concerned with the problem of transmitting a single message in a format which is robust to noise. Recently, researchers have turned their attention to designing coding schemes to make two-way conversations robust to noise. That is, given an interactive communication protocol $\Pi$, an \emph{interactive coding scheme} converts $\Pi$ into another communication protocol $\Pi'$ such that, even if errors are introduced during the execution of $\Pi'$, the parties are able to determine what the outcome of running $\Pi$ would be in a noise-free setting. We consider the problem of designing interactive coding schemes which allow the parties to simulate the original protocol using little memory. Specifically, given any communication protocol $\Pi$ we construct robust simulating protocols which tolerate a constant noise rate and require the parties to use only $O(\log d \log s)$ memory, where $d$ is the depth of $\Pi$ and $s$ is a measure of the size of $\Pi$. Prior to this work, all known coding schemes required the parties to use at least $\Omega(d)$ memory, as the parties were required to remember the transcript of the conversation thus far. Moreover, our coding scheme achieves a communication rate of $1-O(\sqrt{\varepsilon})$ over oblivious channels and $1-O(\sqrt{\varepsilon\log\log\tfrac{1}{\varepsilon}})$ over adaptive adversarial channels, matching the conjecturally optimal rates. Lastly, we point to connections between fault-tolerant circuits and coding for interactive communication with small memory.


title: Modular forms invariant under non-split Cartan subgroups

abstract: In this paper we describe a method for computing a basis for the space of weight $2$ cusp forms invariant under a non-split Cartan subgroup of prime level $p$. As an application we compute, for certain small values of $p$, explicit equations over $\bf Q$ for the canonical embeddings of the associated modular curves.


title: Structural properties of the Jagla fluid

abstract: The structural properties of the Jagla fluid are studied by Monte Carlo (MC) simulations, numerical solutions of integral equation theories, and the (semi-analytical) rational-function approximation (RFA) method. In the latter case, the results are obtained from the assumption (supported by our MC simulations) that the Jagla potential and a potential with a hard core plus an appropriate piecewise constant function lead to practically the same cavity function. The predictions obtained for the radial distribution function, $g(r)$, from this approach are compared against MC simulations and integral equations for the Jagla model, and also for the limiting cases of the triangle-well potential and the ramp potential, with a general good agreement. The analytical form of the RFA in Laplace space allows us to describe the asymptotic behavior of $g(r)$ in a clean way and compare it with MC simulations for representative states with oscillatory or monotonic decay. The RFA predictions for the Fisher--Widom and Widom lines of the Jagla fluid are obtained.


title: NeuralNetwork-Viterbi: A Framework for Weakly Supervised Video Learning

abstract: Video learning is an important task in computer vision and has experienced increasing interest over the recent years. Since even a small amount of videos easily comprises several million frames, methods that do not rely on a frame-level annotation are of special importance. In this work, we propose a novel learning algorithm with a Viterbi-based loss that allows for online and incremental learning of weakly annotated video data. We moreover show that explicit context and length modeling leads to huge improvements in video segmentation and labeling tasks andinclude these models into our framework. On several action segmentation benchmarks, we obtain an improvement of up to 10% compared to current state-of-the-art methods.


title: Composite Asymmetric Dark Matter with a Dark Photon Portal

abstract: Asymmetric dark matter (ADM) is an attractive framework relating the observed baryon asymmetry of the Universe to the dark matter density. A composite particle in a new strong dynamics is a promising candidate for ADM as the strong dynamics naturally explains the ADM mass in the GeV range. Its large annihilation cross section due to the strong dynamics leaves the asymmetric component to be dominant over the symmetric component. In such composite ADM scenarios, the dark sector has a relatively large entropy density in the early Universe. The large dark sector entropy results in the overclosure of the Universe or at best contradicts with the observations of the cosmic microwave background and the successful Big-Bang Nucleosynthesis. Thus, composite ADM models generically require some portal to transfer the entropy of the dark sector into the Standard Model sector. In this paper, we consider a dark photon portal with a mass in the sub-GeV range and kinetic mixing with the Standard Model photon. We investigate the viable parameter space of the dark photon in detail, which can find broad applications to dark photon portal models. We also provide a simple working example of composite ADM with a dark photon portal. Our model is compatible with thermal leptogenesis and $B - L$ symmetry. By taking into account the derived constraints, we show that the parameter space is largely tested by direct detection experiments.


title: Quantum Zeno effect: A possible resolution to the leakage problem in superconducting quantum computing architectures

abstract: We propose to use the continuous version of the quantum Zeno effect to eliminate leakage to higher energy states in superconducting quantum computing architectures based on Josephson phase and flux qubits. We are particularly interested in the application of this approach to the single-step Greenberger-Horne-Zeilinger (GHZ) state protocol described in [A. Galiautdinov and J. M. Martinis, Phys. Rev. A 78, 010305(R) (2008)]. While being conceptually appealing, the protocol was found to be plagued with a number of spectral crowding and leakage problems. Here we argue that by coupling the qubits to a measuring device which continuously monitors leakage to higher energy states (say, to a very lossy resonator of frequency omega = E3(qubit)-E1(qubit), with 1 labeling the ground state of the qubit), we could potentially restrict the multi-qubit system's evolution to its computational subspace, thus circumventing the above mentioned problems.


title: Anisotropic light-shift and magic-polarization of the intercombination line of Dysprosium atoms in a far-detuned dipole trap

abstract: We characterize the anisotropic differential ac-Stark shift for the Dy $626$ nm intercombination transition, induced in a far-detuned $1070$ nm optical dipole trap, and observe the existence of a "magic polarization" for which the polarizabilities of the ground and excited states are equal. From our measurements we extract both the scalar and tensorial components of the dynamic dipole polarizability for the excited state, $\alpha_E^\text{s} = 188 (12)\,\alpha_\text{0}$ and $\alpha_E^\text{t} = 34 (12)\,\alpha_\text{0}$, respectively, where $\alpha_\text{0}$ is the atomic unit for the electric polarizability. We also provide a theoretical model allowing us to predict the excited state polarizability and find qualitative agreement with our observations. Furthermore, we utilize our findings to optimize the efficiency of Doppler cooling of a trapped gas, by controlling the sign and magnitude of the inhomogeneous broadening of the optical transition. The resulting initial gain of the collisional rate allows us, after forced evaporation cooling, to produce a quasi-pure Bose-Einstein condensate of $^{162}$Dy with $3\times 10^4$ atoms.


title: Neural language representations predict outcomes of scientific research

abstract: Many research fields codify their findings in standard formats, often by reporting correlations between quantities of interest. But the space of all testable correlates is far larger than scientific resources can currently address, so the ability to accurately predict correlations would be useful to plan research and allocate resources. Using a dataset of approximately 170,000 correlational findings extracted from leading social science journals, we show that a trained neural network can accurately predict the reported correlations using only the text descriptions of the correlates. Accurate predictive models such as these can guide scientists towards promising untested correlates, better quantify the information gained from new findings, and has implications for moving artificial intelligence systems from predicting structures to predicting relationships in the real world.


title: It's all Relative: Monocular 3D Human Pose Estimation from Weakly Supervised Data

abstract: We address the problem of 3D human pose estimation from 2D input images using only weakly supervised training data. Despite showing considerable success for 2D pose estimation, the application of supervised machine learning to 3D pose estimation in real world images is currently hampered by the lack of varied training images with corresponding 3D poses. Most existing 3D pose estimation algorithms train on data that has either been collected in carefully controlled studio settings or has been generated synthetically. Instead, we take a different approach, and propose a 3D human pose estimation algorithm that only requires relative estimates of depth at training time. Such training signal, although noisy, can be easily collected from crowd annotators, and is of sufficient quality for enabling successful training and evaluation of 3D pose algorithms. Our results are competitive with fully supervised regression based approaches on the Human3.6M dataset, despite using significantly weaker training data. Our proposed algorithm opens the door to using existing widespread 2D datasets for 3D pose estimation by allowing fine-tuning with noisy relative constraints, resulting in more accurate 3D poses.


title: Changing Observations in Epistemic Temporal Logic

abstract: We study dynamic changes of agents' observational power in logics of knowledge and time. We consider CTL*K, the extension of CTL* with knowledge operators, and enrich it with a new operator that models a change in an agent's way of observing the system. We extend the classic semantics of knowledge for perfect-recall agents to account for changes of observation, and we show that this new operator strictly increases the expressivity of CTL*K. We reduce the model-checking problem for our logic to that for CTL*K, which is known to be decidable. This provides a solution to the model-checking problem for our logic, but its complexity is not optimal. Indeed we provide a direct decision procedure with better complexity.


title: Time-reversal symmetry-breaking nematic superconductivity in FeSe

abstract: FeSe is a unique member of the family of iron-based superconductors, not only because of the high values of $T_{c}$ in FeSe monolayer, but also because in bulk FeSe superconductivity emerges inside a nematic phase without competing with long-range magnetic order. Near $T_{c}$, superconducting order necessarily has $s+d$ symmetry, because nematic order couples linearly the $s$-wave and $d$-wave harmonics of the superconducting order parameter. Here we argue that the near-degeneracy between $s$-wave and $d$-wave pairing instabilities in FeSe, combined with the sign-change of the nematic order parameter between hole and electron pockets, allows the superconducting order to break time-reversal symmetry at a temperature $T^{*}<T_{c}$. The transition from an $s+d$ state to an $s+\mathrm{e}^{i\alpha}d$ state should give rise to a peak in the specific heat and to the emergence of a soft collective mode that can be potentially detected by Raman spectroscopy.


title: Particle acceleration in relativistic magnetic flux-merging events

abstract: Using analytical and numerical methods (fluid and particle-in-cell simulations) we study a number of model problems involving merger of magnetic flux tubes in relativistic magnetically-dominated plasma. Mergers of current-carrying flux tubes (exemplified by the two dimensional `ABC' structures) and zero total current magnetic flux tubes are considered. In all cases regimes of spontaneous and driven evolution are investigated. We identify two stages of particle acceleration during flux mergers: (i) fast explosive prompt X-point collapse and (ii) ensuing island merger. The fastest acceleration occurs during the initial catastrophic X-point collapse, with the reconnection electric field of the order of the magnetic field. During the X-point collapse particles are accelerated by charge-starved electric fields, which can reach (and even exceed) values of the local magnetic field. The explosive stage of reconnection produces non-thermal power-law tails with slopes that depend on the average magnetization $\sigma$. For plasma magnetization $\sigma \leq 10^2$ the spectrum power law index is $p> 2$; in this case the maximal energy depends linearly on the size of the reconnecting islands. For higher magnetization, $\sigma \geq 10^2$, the spectra are hard, $p< 2$, yet the maximal energy $\gamma_{max}$ can still exceed the average magnetic energy per particle, $ \sim \sigma$, by orders of magnitude (if $p$ is not too close to unity). The X-point collapse stage is followed by magnetic island merger that dissipates a large fraction of the initial magnetic energy in a regime of forced magnetic reconnection, further accelerating the particles, but proceeds at a slower reconnection rate.


title: Individual Control and Readout of Qubits in a Sub-Diffraction Volume

abstract: Medium-scale ensembles of coupled qubits offer a platform for near-term quantum technologies including computing, sensing, and the study of mesoscopic quantum systems. Atom-like emitters in solids have emerged as promising quantum memories, with demonstrations of spin-spin entanglement by optical and magnetic interactions. Magnetic coupling in particular is attractive for efficient and deterministic entanglement gates, but raises the problem of individual spin addressing at the necessary nanometer-scale separation. Current super-resolution techniques can reach this resolution, but are destructive to the states of nearby qubits. Here, we demonstrate the measurement of individual qubit states in a sub-diffraction cluster by selectively exciting spectrally distinguishable nitrogen vacancy (NV) centers. We demonstrate super-resolution localization of single centers with nanometer spatial resolution, as well as individual control and readout of spin populations. These measurements indicate a readout-induced crosstalk on non-addressed qubits below $4\times10^{-2}$. This approach opens the door to high-speed control and measurement of qubit registers in mesoscopic spin clusters, with applications ranging from entanglement-enhanced sensors to error-corrected qubit registers to multiplexed quantum repeater nodes.


title: Anomaly matching and symmetry-protected critical phases in $SU(N)$ spin systems in 1+1 dimensions

abstract: We study $(1+1)$-dimensional $SU(N)$ spin systems in the presence of the global $SU(N)$ rotation and lattice translation symmetries. By matching the mixed anomaly of the $PSU(N)\times\mathbb{Z}$ symmetry in the continuum limit, we identify a topological index for spin model evaluated as the total number of Young-tableaux boxes of spins per unit cell modulo $N$, which characterizes the "ingappability" of the system. A nontrivial index implies either a ground-state degeneracy in a gapped phase, which can be regarded as a field-theory version of the Lieb-Schultz-Mattis theorem, or a restriction of the possible universality classes in a critical phase -- the symmetry-protected critical phase, e.g. only a class of $SU(N)$ Wess-Zumino-Witten theories can be realized in the low-energy limit of the given lattice model in the presence of the symmetries. Similar constraints also apply when a higher global symmetry emerges in the model with a lower symmetry. Our prediction agrees with several examples known in previous studies of $SU(N)$ models.

categories: cond-mat.str-el cond-mat.quant-gas cond-mat.stat-mech hep-th


title: Chiral symmetry restoration by parity doubling and the structure of neutron stars

abstract: We investigate the equation of state for a recently developed hybrid quark-meson-nucleon model under neutron star conditions of $\beta-$equilibrium and charge neutrality. The model has the characteristic feature that at increasing baryon density chiral symmetry is restored in a first order transition within the hadronic phase by lifting the mass splitting between chiral partner states, before quark deconfinement takes place. Most important for this study are the nucleon (neutron, proton) and $N(1535)$ states. We present three sets for the two free parameters which result in compact star mass-radius relations in accordance with modern constraints on the mass from PSR~J0437-4715 and on the compactness from GW170817. We also consider the threshold for the direct URCA process for which a new relationship is given and suggest as an additional constraint on the parameter choice of the model that this process shall become operative at best for stars with masses above the range for binary radio pulsars, $M>1.4~M_\odot$.


title: Comparison of \nu\mu-Ar multiplicity distributions observed by MicroBooNE to GENIE model predictions

abstract: We measure a large set of observables in inclusive charged current muon neutrino scattering on argon with the MicroBooNE liquid argon time projection chamber operating at Fermilab. We evaluate three neutrino interaction models based on the widely used GENIE event generator using these observables. The measurement uses a data set consisting of neutrino interactions with a final state muon candidate fully contained within the MicroBooNE detector. These data were collected in 2016 with the Fermilab Booster Neutrino Beam, which has an average neutrino energy of 800 MeV, using an exposure corresponding to 5E19 protons-on-target. The analysis employs fully automatic event selection and charged particle track reconstruction and uses a data-driven technique to separate neutrino interactions from cosmic ray background events. We find that GENIE models consistently describe the shapes of a large number of kinematic distributions for fixed observed multiplicity.


title: Star Cluster Disruption by a Supermassive Black Hole Binary

abstract: Supermassive black hole binaries (BHBs) are expected to be one of the most powerful sources of low-frequency gravitational waves (GWs) for future space-borne detectors. Prior to the GW emission stage, BHBs evolving in gas-poor nuclei shrink primarily through the slingshot ejection of stars approaching the BHB from sufficiently close distances. Here we address the possibility that the BHB shrinking rate is enhanced through the infall of a star cluster (SC) onto the BHB. We present the results of direct summation N-body simulations exploring different orbits for the SC infall, and we show that SCs reaching the BHB on non-zero angular momentum orbits (with eccentricity 0.75) fail to enhance the BHB hardening, while SCs approaching the BHB on radial orbits reduce the BHB separation by ~10% in less than 10 Myr, effectively shortening the BHB path towards GWs.


title: Probing star formation and ISM properties using galaxy disk inclination II: Testing typical FUV attenuation corrections out to z$\sim$0.7

abstract: We evaluate dust-corrected far ultraviolet (FUV) star formation rates (SFRs) for samples of star-forming galaxies at $z\sim0$ and $z\sim0.7$ and find significant differences between values obtained through corrections based on UV colour, from a hybrid mid-infrared (MIR) plus FUV relation, and from a radiative transfer based attenuation correction method. The performances of the attenuation correction methods are assessed by their ability to remove the dependency of the corrected SFR on inclination, as well as returning, on average, the expected population mean SFR. We find that combining MIR (rest-frame $\sim$13$\mu$m) and FUV luminosities gives the most inclination independent SFRs and reduces the intrinsic SFR scatter out of the methods tested. However, applying the radiative transfer based method of Tuffs et al. gives corrections to the FUV SFR that are inclination independent and in agreement with the expected SFRs at both $z\sim0$ and $z\sim0.7$. SFR corrections based on the UV-slope perform worse than the other two methods tested. For our local sample, the UV-slope method works on average but does not remove inclination biases. At $z\sim$0.7 we find that the UV-slope correction used locally flattens the inclination dependence compared to the raw FUV measurements but was not sufficient to correct for the large attenuation observed at $z\sim$0.7.


title: Walter Talbot's thesis

abstract: Walter Richard Talbot was the fourth African American to earn a PhD in Mathematics. His doctoral degree is from the University of Pittsburgh in 1934 in geometric group theory. A contemporary research program was the determination of fundamental domains of finite group actions on complex vector spaces. His thesis is not widely available, and this note gives a brief synopsis of the main results of his thesis, expressed using modern mathematical methods and language, and placed in general context.


title: Analyzing the love affair of Romeo and Juliet with modern mathematical tools

abstract: We facetiously suggest that the romance between Romeo and Juliet can be interpreted using modern terminology and include current temptations. Using this model, we consider various factors such as the time that they might spend consulting social networks, the time that they could spend alone together and along with friends as well as their tolerance of being able to waste the couple's money. The model consists of a set of differential equations which describes the relationship between them. Finally, we analyze the eigenvalues the mathematical equations in order to determine if the critical point in this model is stable or not in four different hypothetical scenarios.


title: Empirically determined dilution factors of stripped-envelope, core-collapse SNe: Paper II - Using GRB-SNe to determine the Hubble Constant

abstract: The aim of this work is to use gamma-ray burst supernovae (GRB-SNe) as cosmological probes to measure the Hubble constant, $H_0$, in the local Universe. In the context of the Expanding Photosphere Method (EPM), I use empirically derived dilution factors of a sample of nearby SNe Ic, which were derived in Paper I of a two-paper series, as a proxy for the dilution factors of GRB-SNe. It is seen that the dilution factors as a function of temperature in $VI$ display the least amount of scatter, relative to $BVI$ and $BV$. A power-law function is fit to the former, and is used to derive model dilution factors which are then used to derive EPM distances to GRB-SNe 1998bw and 2003lw: $36.7\pm9.6$ and $372.2\pm137.1$ Mpc, respectively. In linear Hubble diagrams in filters $BVR$, I determine the offset of the Hubble ridge line, and armed with the peak absolute magnitudes in these filters for the two aforementioned GRB-SNe, I find a (weighted average) Hubble constant of $\bar{H_{0,\rm w}} = 61.9\pm12.3$ km s$^{-1}$ Mpc$^{-1}$ for GRB-SNe located at redshifts $z\le0.1$. The 20\% error is consistent with the value of $H_0$ calculated by Planck and SNe Ia within 1$\sigma$. I tested the fitting method on five nearby SNe Ic, and found that their EPM distances varied by 18-50\%, with smaller errors found for those SNe which had more numerous usable observations. For SN 2002ap, its EPM distance was overestimated by 18\%, and if the distance to SN 1998bw was similarly over-estimated by the same amount, the resultant value of the Hubble constant is $H_0 = 72$ km s$^{-1}$ Mpc$^{-1}$, which perfectly matches that obtained using SNe Ia. [abridged]


title: Java Subtyping as an Infinite Self-Similar Partial Graph Product

abstract: Due to supporting variance annotations, such as wildcard types, the subtyping relation in Java and other generic nominally-typed OO programming languages is both interesting and intricate. In these languages, the subtyping relation between ground object types, i.e., ones with no type variables, is the basis for defining the full OO subtyping relation, i.e., that includes type variables. As an ordering relation over the set of types, the subtyping relation in object-oriented programming languages can always be represented as a directed graph. In order to better understand some of the subtleties of the subtyping relation in Java, in this paper we present how the subtyping relation between ground Java types can be precisely constructed using two new operations (a binary operation and a unary one) on directed graphs. The binary operation we use, called a partial Cartesian graph product, is similar in its essence to standard graph products and group products. Its definition is based in particular on that of the standard Cartesian graph product. We believe the use of graph operations in constructing the ground generic Java subtyping relation reveals some of the not-immediately-obvious structure of the subtyping relation not only in Java but, more generally, also in mainstream generic nominally-typed OO programming languages such as C#, Scala and Kotlin. Accordingly, we believe that describing precisely how graph operations can be used to explicitly construct the subtyping relation in these languages, as we do in this paper, may significantly improve our understanding of features of the type systems of these languages such as wildcard types and variance annotations, and of the dependency of these features on nominal subtyping in nominally-typed OOP.


title: Kinematics of outer halo globular clusters: M 75 and NGC 6426

abstract: Globular clusters (GCs) and their dynamic interactions with the Galactic components provide an important insight into the structure and formation of the early Milky Way. Here, we present a kinematic study of two outer halo GCs based on a combination of VLT/FORS2, VLT/FLAMES, and Magellan/MIKE low- and high-resolution spectroscopy of 32 and 27 member stars, respectively. Although both clusters are located at Galactocentric distances of 15 kpc, they have otherwise very different properties. M 75 is a luminous and metal-rich system at [Fe/H] = $-1.2$ dex, a value that we confirm from the calcium triplet region. This GC shows mild evidence for rotation with an amplitude of A$_{\rm rot}\sim$5 km s$^{-1}$. One of the most metal-poor GCs in the Milky Way (at [FeII/H] = $-2.3$ dex), NGC 6426 exhibits marginal evidence of internal rotation at the 2 km s$^{-1}$ level. Both objects have velocity dispersions that are consistent with their luminosity. Although limited by small-number statistics, the resulting limits on their $A_{\rm rot}/\sigma_0$ ratios suggest that M 75 is a slow rotator driven by internal dynamics rather than being effected by the weak Galactic tides at its large distances. Here, M 75 ($A_{\rm rot}/\sigma=0.31$) is fully consistent with the properties of other, younger halo clusters. At $A_{\rm rot}/\sigma_0=0.8\pm0.4$, NGC 6426 appears to have a remarkably ordered internal motion for its low metallicity, but the large uncertainty does not allow for an unambiguous categorization as a fast rotator. An accretion origin of M 75 cannot be excluded, based on the eccentric orbit, which we derived from the recent data release 2 of Gaia, and considering its younger age.


title: Emergent locality in systems with power-law interactions

abstract: Locality imposes stringent constraints on the spreading of information in nonrelativistic quantum systems, which is reminiscent of a "light-cone," a casual structure arising in their relativistic counterparts. Long-range interactions can potentially soften such constraints, allowing almost instantaneous long jumps of particles, thus defying causality. Since interactions decaying as a power-law with distance, $r^{-\alpha}$, are ubiquitous in nature, it is pertinent to understand what is the fate of causality and information spreading in such systems. Using a numerically exact technique we address these questions by studying the out-of-time-order correlation function of a representative generic system in one-dimension. We show that while the interactions are long-range, their effect on information spreading is asymptotically negligible as long as $\alpha>1$. In this range we find a complex compound behavior, where after a short transient a fully local behavior emerges, yielding asymptotic "light-cones" virtually indistinguishable from "light-cones" in corresponding local models. The long-range nature of the interaction is only expressed in the power-law leaking of information from the "light-cone," with the same exponent as the exponent of the interaction, $\alpha$. Our results directly imply that all previously obtained rigorous bounds on information spreading in long-range interacting systems are not tight, and thus could be improved.


title: The Formation of Ultra Diffuse Galaxies in Cored Dark Matter Halos Through Tidal Stripping and Heating

abstract: We propose that the Ultra-Diffuse Galaxy (UDG) population represents a set of satellite galaxies born in $\sim10^{10}-10^{11}$ M$_\odot$ halos, similar to field dwarfs, which suffer a dramatic reduction in surface brightness due to tidal stripping and heating. This scenario is observationally motivated by the radial alignment of UDGs in Coma as well as the significant dependence of UDG abundance on cluster mass. As a test of this formation scenario, we apply a semi-analytic model describing the change in stellar mass and half-light radius of dwarf satellites, occupying either cored or cuspy halos, to cluster subhalos in the Illustris-dark simulation. Key to this model are results from simulations which indicate that galaxies in cored dark-matter halos expand significantly in response to tidal stripping and heating, whereas galaxies in cuspy halos experience limited size evolution. Our analysis indicates that a population of tidally-stripped dwarf galaxies, residing in cored halos (like those hosting low-surface brightness field dwarfs), is able to reproduce the observed sizes and stellar masses of UDGs in clusters remarkably well.


title: SNe Kicks in Hierarchical Triple Systems

abstract: Most massive stars, if not all, are in binary configuration or higher multiples. These massive stars undergo supernova explosions and end their lives as either black holes or neutron stars. Recent observations have suggested that neutron stars and perhaps even black holes receive large velocity kicks at birth. Such natal kicks and the sudden mass loss can significantly alter the orbital configuration of the system. Here we derive general analytical expressions that describe the effects of natal kicks in binaries on hierarchical triple systems. We explore several proof-of-concept applications such as black hole and neutron stars binaries and X-ray binaries with either stellar or Supermassive Black Hole (SMBH) companions on a wide orbit. Kicks can disrupt the hierarchical configuration, although it is harder to escape the potential well of an SMBH. Some binary systems do escape the SMBH system resulting in hyper-velocity binary system. Furthermore, kicks can result in increasing or decreasing the orbital separations. Decreasing the orbital separation may have significant consequences in these astrophysical systems. For example, shrinking the separation post-supernova kick can lead to the shrinking of an inner compact binary that then may merge via Gravitational Wave (GW) emission. This process yields a supernova that is shortly followed by a possible GW-LIGO event. Interestingly, we find that the natal kick can result in shrinking the outer orbit, and the binary may cross the tertiary Roche limit, breaking up the inner binary. Thus, in the case of SMBH companion, this process can lead to either a tidal disruption event or a GW-LISA detection event with a supernova precursor.


title: Gas-Assisted Growth of Protoplanets in a Turbulent Medium

abstract: Pebble accretion is a promising process for decreasing growth timescales of planetary cores, allowing gas giants to form at wide orbital separations. However, nebular turbulence can reduce the efficiency of this gas-assisted growth. We present an order of magnitude model of pebble accretion, which calculates the impact of turbulence on the average velocity of small bodies, the radius for binary capture, and the sizes of the small bodies that can be accreted. We also include the effect of turbulence on the particle scale height, which has been studied in previous works. We find that turbulence does not prevent rapid growth in the high-mass regime: the last doubling time to the critical mass to trigger runaway gas accretion ($M \sim 10 \, M_\oplus$) is well within the disk lifetime even for strong ($\alpha \gtrsim 10^{-2}$) turbulence. We find that while the growth timescale is quite sensitive to the local properties of the protoplanetary disk, over large regimes of parameter space large cores grow in less than the disk lifetime if appropriately-sized small bodies are present. Instead, the effects of turbulence are most pronounced for low planetary masses. For strong turbulence the growth timescale is longer than the gas disk lifetime until the core reaches masses $\gtrsim 10^{-2}-10^{-1} M_\oplus$. A "Flow Isolation Mass", at which binary capture ceases, emerges naturally from our model framework. We comment that the dependence of this mass on orbital separation is similar to the semi-major axis distribution of solar-system cores.


title: Twisted Fracton Models in Three Dimensions

abstract: We study novel three-dimensional gapped quantum phases of matter which support quasiparticles with restricted mobility, including immobile "fracton" excitations. So far, most existing fracton models may be instructively viewed as generalized Abelian lattice gauge theories. Here, by analogy with Dijkgraaf-Witten topological gauge theories, we discover a natural generalization of fracton models, obtained by twisting the gauge symmetries. Introducing generalized gauge transformation operators carrying an extra phase factor depending on local configurations, we construct a plethora of exactly solvable three-dimensional models, which we dub "twisted fracton models." A key result of our approach is to demonstrate the existence of rich non-Abelian fracton phases of distinct varieties in a three-dimensional system with finite-range interactions. For an accurate characterization of these novel phases, the notion of being inextricably non-Abelian is introduced for fractons and quasiparticles with one-dimensional mobility, referring to their new behavior of displaying braiding statistics that is, and remains, non-Abelian regardless of which quasiparticles with higher mobility are added to or removed from them. We also analyze these models by embedding them on a three-torus and computing their ground state degeneracies, which exhibit a surprising and novel dependence on the system size in the non-Abelian fracton phases. Moreover, as an important advance in the study of fracton order, we develop a general mathematical framework which systematically captures the fusion and braiding properties of fractons and other quasiparticles with restricted mobility.

categories: cond-mat.str-el cond-mat.stat-mech hep-th quant-ph


title: Quantum gravity: a quantum-first approach

abstract: A "quantum-first" approach to gravity is described, where rather than quantizing general relativity, one seeks to formulate the physics of gravity within a quantum-mechanical framework with suitably general postulates. Important guides are the need for appropriate mathematical structure on Hilbert space, and correspondence with general relativity and quantum field theory in weak-gravity situations. A basic physical question is that of "Einstein separability:" how to define mutually independent subsystems, e.g. through localization. Standard answers via tensor products or operator algebras conflict with properties of gravity, as is seen in the correspondence limit; this connects with discussions of "soft hair." Instead, gravitational behavior suggests a networked Hilbert space structure. This structure plus unitarity provide important clues towards a quantum formulation of gravity.


title: Gravitational Wave Signals from the First Massive Black Hole Seeds

abstract: Recent numerical simulations reveal that the isothermal collapse of pristine gas in atomic cooling haloes may result in stellar binaries of supermassive stars with $M_* \gtrsim 10^4\ \mathrm{M}_{\odot}$. For the first time, we compute the in-situ merger rate for such massive black hole remnants by combining their abundance and multiplicity estimates. For black holes with initial masses in the range $10^{4-6} \ \mathrm{M}_{\odot}$ merging at redshifts $z \gtrsim 15$ our optimistic model predicts that LISA should be able to detect 0.6 mergers per year. This rate of detection can be attributed, without confusion, to the in-situ mergers of seeds from the collapse of very massive stars. Equally, in the case where LISA observes no mergers from heavy seeds at $z \gtrsim 15$ we can constrain the combined number density, multiplicity, and coalesence times of these high-redshift systems. This letter proposes gravitational wave signatures as a means to constrain theoretical models and processes that govern the abundance of massive black hole seeds in the early Universe.


title: Cosmological simulations of black hole growth II: how (in)significant are merger events for fuelling nuclear activity?

abstract: Which mechanism(s) are mainly driving nuclear activity in the centres of galaxies is a major unsettled question. In this study, we investigate the statistical relevance of galaxy mergers for fuelling gas onto the central few kpc of a galaxy, potentially resulting in an active galactic nucleus (AGN). To robustly address that, we employ large-scale cosmological hydrodynamic simulations from the Magneticum Pathfinder set, including models for BH accretion and AGN feedback. Our simulations predict that for luminous AGN ($L_{\rm AGN} > 10^{45} {\rm erg/s}$) at $z = 2$, more than 50 per cent of their host galaxies have experienced a merger in the last 0.5~Gyr. These high merger fractions, however, merely reflect the intrinsically high merger fractions of massive galaxies at $z=2$, in which luminous AGN preferentially occur. Apart from that, our simulations suggest that merger events are not the statistically dominant fuelling mechanism for nuclear activity over a redshift range $z=0-2$: irrespective of AGN luminosity, less than 20 per cent of AGN hosts have on average undergone a recent merger, in agreement with a number of observational studies. The central ISM conditions required for inducing AGN activity can be, but are not necessarily caused by a merger. Despite the statistically minor relevance of mergers, at a given AGN luminosity and stellar mass, the merger fractions of AGN hosts can be by up to three times higher than that of inactive galaxies. Such elevated merger fractions still point towards an intrinsic connection between AGN and mergers, consistent with our traditional expectation.


title: A UV Complete Framework of Freeze-in Massive Particle Dark Matter

abstract: We propose a way to generate tiny couplings of freeze-in massive particle dark matter with the Standard Model particles dynamically by considering an extension of the electroweak gauge symmetry. The dark matter is considered to be a singlet under this extended gauge symmetry which we have assumed to be the one in a very widely studied scenario called left-right symmetric model. Several heavy particles, that can be thermally inaccessible in the early Universe due to their masses being greater than the reheat temperature after inflation, can play the role of portals between dark matter and Standard Model particles through one loop couplings. Due to the loop suppression, one can generate the required non-thermal dark matter couplings without any need of highly fine tuned Yukawa couplings beyond that of electron Yukawa with the Standard Model like Higgs boson. We show that generic values of Yukawa couplings as large as $\mathcal{O}(0.01)$ to $\mathcal{O}(1)$ can keep the dark matter out of thermal equilibrium in the early Universe and produce the correct relic abundance later through the freeze-in mechanism. Though the radiative couplings of dark matter are tiny as required by the freeze-in scenario, the associated rich particle sector of the model can be probed at ongoing and near future experiments. The allowed values of dark matter mass can remain in a wide range from keV to TeV order keeping the possibilities of warm and cold dark matter equally possible.


title: Asymptotic Correlations in Gapped and Critical Topological Phases of 1D Quantum Systems

abstract: Topological phases protected by symmetry can occur in gapped and---surprisingly---in critical systems. We consider non-interacting fermions in one dimension with spinless time-reversal symmetry. It is known that the phases are classified by a topological invariant $\omega$ and a central charge $c$. We investigate the correlations of string operators, giving insight into the interplay between topology and criticality. In the gapped phases, these non-local string order parameters allow us to extract $\omega$. Remarkably, ratios of correlation lengths are universal. In the critical phases, the scaling dimensions of these operators serve as an order parameter, encoding $\omega$ and $c$. We derive exact asymptotics of these correlation functions using Toeplitz determinant theory. We include physical discussion, e.g., relating lattice operators to the conformal field theory. Moreover, we discuss the dual spin chains. Using the aforementioned universality, the topological invariant of the spin chain can be obtained from correlations of local observables.

categories: cond-mat.str-el cond-mat.stat-mech math-ph math.MP quant-ph


title: A Theory for the Variation of Dust Attenuation Laws in Galaxies

abstract: In this paper, we provide a physical model for the origin of variations in the shapes and bump strengths of dust attenuation laws in galaxies by combining a large suite of cosmological "zoom-in" galaxy formation simulations with 3D Monte Carlo dust radiative transfer calculations. We model galaxies over 3 orders of magnitude in stellar mass, ranging from Milky Way like systems through massive galaxies at high-redshift. Critically, for these calculations we employ a constant underlying dust extinction law in all cases, and examine how the role of geometry and radiative transfer effects impact the resultant attenuation curves. Our main results follow. Despite our usage of a constant dust extinction curve, we find dramatic variations in the derived attenuation laws. The slopes of normalized attenuation laws depend primarily on the complexities of star-dust geometry. Increasing fractions of unobscured young stars flatten normalized curves, while increasing fractions of unobscured old stars steepen curves. Similar to the slopes of our model attenuation laws, we find dramatic variation in the 2175 Angstrom ultraviolet (UV) bump strength, including a subset of curves with little to no bump. These bump strengths are primarily influenced by the fraction of unobscured O and B stars in our model, with the impact of scattered light having only a secondary effect. Taken together, these results lead to a natural relationship between the attenuation curve slope and 2175 Angstrom bump strength. Finally, we apply these results to a 25 Mpc/h box cosmological hydrodynamic simulation in order to model the expected dispersion in attenuation laws at integer redshifts from z=0-6. A significant dispersion is expected at low redshifts, and decreases toward z=6. We provide tabulated results for the best fit median attenuation curve at all redshifts.


title: Emergent Superconductivity in the weak Mott insulator phase of bilayer Graphene Moir\'e Superlattice

abstract: We propose a phenomenological understanding of the recently discovered weak Mott insulator in the moir\'e superlattice of twisted bilayer graphene, especially the emergent superconductivity at low temperature within the weak Mott insulator phase, namely while lowering temperature, the longitudinal resistivity first grows below temperature $T_m$, but then rapidly drops to zero at even lower temperature $T_c$. An emergent superconductor in an insulator phase is very unusual. Here we propose that this phenomenon is due to the pure two-dimensional nature of the bilayer graphene moir\'e superlattice. We also compare our results with other theories proposed so far.


title: Type Ia supernovae have two physical width-luminosity relations and they favor sub-Chandrasekhar and direct collision models - II. Color evolution

abstract: While the width-luminosity relation (WLR) among type Ia supernovae (slower is brighter) is one of the best studied properties of this type of events, its physical basis has not been identified convincingly. The 'luminosity' is known to be related to a clear physical quantity -- the amount of $^{56}$Ni synthesized, but the 'width' has not been quantitatively linked yet to a physical time scale. We show that the recombination time of $^{56}$Fe and $^{56}$Co from doubly to singly ionized states causes the typical observed break in the color curve B-V due to a cliff in the mean opacities, and is a robust width measure of the light curve, which is insensitive to radiation transfer uncertainties. A simple photospheric model is shown to predict the recombination time to an accuracy of $\sim5$ days, allowing a quantitative understanding of the color WLR. Two physical times scales of the width luminosity relation are shown to be set by two column densities -- the total column density which sets the gamma-ray escape time $t_0$ (previous Paper I) and the $^{56}$Ni column density which sets the recombination time (this Paper II). Central detonations of sub-$\rm M_{ch}$ WDs and direct WD collision models have gamma-ray escape times and recombination times which are consistent with observations across the luminosity range of type Ia's. Delayed detonation Chandrasekhar mass models have recombination times that are broadly consistent with observations, with tension at the bright end of the luminosity range and inconsistent gamma-ray escape times at the faint end.


title: Towards Large-scale Functional Verification of Universal Quantum Circuits

abstract: We introduce a framework for the formal specification and verification of quantum circuits based on the Feynman path integral. Our formalism, built around exponential sums of polynomial functions, provides a structured and natural way of specifying quantum operations, particularly for quantum implementations of classical functions. Verification of circuits over all levels of the Clifford hierarchy with respect to either a specification or reference circuit is enabled by a novel rewrite system for exponential sums with free variables. Our algorithm is further shown to give a polynomial-time decision procedure for checking the equivalence of Clifford group circuits. We evaluate our methods by performing automated verification of optimized Clifford+T circuits with up to 100 qubits and thousands of T gates, as well as the functional verification of quantum algorithms using hundreds of qubits. Our experiments culminate in the automated verification of the Hidden Shift algorithm for a class of Boolean functions in a fraction of the time it has taken recent algorithms to simulate.


title: Fully Convolutional Model for Variable Bit Length and Lossy High Density Compression of Mammograms

abstract: Early works on medical image compression date to the 1980's with the impetus on deployment of teleradiology systems for high-resolution digital X-ray detectors. Commercially deployed systems during the period could compress 4,096 x 4,096 sized images at 12 bpp to 2 bpp using lossless arithmetic coding, and over the years JPEG and JPEG2000 were imbibed reaching upto 0.1 bpp. Inspired by the reprise of deep learning based compression for natural images over the last two years, we propose a fully convolutional autoencoder for diagnostically relevant feature preserving lossy compression. This is followed by leveraging arithmetic coding for encapsulating high redundancy of features for further high-density code packing leading to variable bit length. We demonstrate performance on two different publicly available digital mammography datasets using peak signal-to-noise ratio (pSNR), structural similarity (SSIM) index and domain adaptability tests between datasets. At high density compression factors of >300x (~0.04 bpp), our approach rivals JPEG and JPEG2000 as evaluated through a Radiologist's visual Turing test.


title: Investigating the scalar sector of left-right symmetric models with leptonic probes

abstract: We investigate the potential collider signatures of singly-charged and doubly-charged Higgs bosons such as those arising in minimal left-right symmetric models. Focusing on multileptonic probes in the context of the high-luminosity run of the Large Hadron Collider, we separately assess the advantages of the four-leptonic and trileptonic final states for a representative benchmark setup designed by considering a large set of experimental constraints. Our study establishes possibilities of identifying singly-charged and doubly-charged scalars at the Large Hadron Collider with a large significance, for luminosity goals expected to be reached during the high-luminosity phase of the Large Hadron Collider. We generalise our results and demonstrate that existing limits can in principle be pushed much further in the heavy mass regime.


title: On the Capacity of MIMO Broadband Power Line Communications Channels

abstract: Communications over power lines in the frequency range above 2 MHz, commonly referred to as broadband (BB) power line communications (PLC), has been the focus of increasing research attention and standardization efforts in recent years. BB-PLC channels are characterized by a dominant colored non-Gaussian additive noise, as well as by periodic variations of the channel impulse response and the noise statistics. In this work we study the fundamental rate limits for BB-PLC channels by bounding their capacity while accounting for the unique properties of these channels. We obtain explicit expressions for the derived bounds for several BB-PLC noise models, and illustrate the resulting fundamental limits in a numerical analysis.


title: Fast reinforcement learning for decentralized MAC optimization

abstract: In this paper, we propose a novel decentralized framework for optimizing the transmission strategy of Irregular Repetition Slotted ALOHA (IRSA) protocol in sensor networks. We consider a hierarchical communication framework that ensures adaptivity to changing network conditions and does not require centralized control. The proposed solution is inspired by the reinforcement learning literature, and, in particular, Q-learning. To deal with sensor nodes' limited lifetime and communication range, we allow them to decide how many packet replicas to transmit considering only their own buffer state. We show that this information is sufficient and can help avoiding packets' collisions and improving the throughput significantly. We solve the problem using the decentralized partially observable Markov Decision Process (Dec-POMDP) framework, where we allow each node to decide independently of the others how many packet replicas to transmit. We enhance the proposed Q-learning based method with the concept of virtual experience, and we theoretically and experimentally prove that convergence time is, thus, significantly reduced. The experiments prove that our method leads to large throughput gains, in particular when network traffic is heavy, and scales well with the size of the network. To comprehend the effect of the problem's nature on the learning dynamics and vice versa, we investigate the waterfall effect, a severe degradation in performance above a particular traffic load, typical for codes-on-graphs and prove that our algorithm learns to alleviate it.


title: Spectroscopic decomposition of the galaxy and halo of the cD galaxy NGC 3311

abstract: Information on the star-formation histories of cD galaxies and their extended stellar haloes lie in their spectra. Therefore, to determine whether these structures evolved together or through a two-phase formation, we need to spectroscopically separate the light from each component. We present a pilot study to use BUDDI to fit and extract the spectra of the cD galaxy NGC 3311 and its halo in an Integral Field Spectroscopy datacube, and carry out a simple stellar populations analysis to study their star-formation histories. Using MUSE data, we were able to isolate the light of the galaxy and its halo throughout the datacube, giving spectra representing purely the light from each of these structures. The stellar populations analysis of the two components indicates that, in this case, the bulk of the stars in both the halo and the central galaxy are very old, but the halo is more metal poor and less $\alpha$-enriched than the galaxy. This result is consistent with the halo forming through the accretion of much smaller satellite galaxies with more extended star formation. It is noteworthy that the apparent gradients in age and metallicity indicators across the galaxy are entirely consistent with the radially-varying contributions of galaxy and halo components, which individually display no gradients. The success of this study is promising for its application to a larger sample of cD galaxies that are currently being observed by IFU surveys.


title: Newton Polygons Arising for Special Families of Cyclic Covers of the Projective Line

abstract: By a result of Moonen, there are exactly 20 positive-dimensional families of cyclic covers of the projective line for which the Torelli image is open and dense in the associated Shimura variety. For each of these, we compute the Newton polygons, and the $\mu$-ordinary Ekedahl--Oort type, occurring in the characteristic $p$ reduction of the Shimura variety. We prove that all but a few of the Newton polygons appear on the open Torelli locus. As an application, we produce multiple new examples of Newton polygons and Ekedahl--Oort types of Jacobians of smooth curves in characteristic $p$. Under certain congruence conditions on $p$, these include: the supersingular Newton polygon for genus $5,6,7$; fourteen new non-supersingular Newton polygons for genus $5-7$; eleven new Ekedahl--Oort types for genus $4-7$ and, for all $g \geq 6$, the Newton polygon with $p$-rank $g-6$ with slopes $1/6$ and $5/6$.


title: A Note on Coding and Standardization of Categorical Variables in (Sparse) Group Lasso Regression

abstract: Categorical regressor variables are usually handled by introducing a set of indicator variables, and imposing a linear constraint to ensure identifiability in the presence of an intercept, or equivalently, using one of various coding schemes. As proposed in Yuan and Lin [J. R. Statist. Soc. B, 68 (2006), 49-67], the group lasso is a natural and computationally convenient approach to perform variable selection in settings with categorical covariates. As pointed out by Simon and Tibshirani [Stat. Sin., 22 (2011), 983-1001], "standardization" by means of block-wise orthonormalization of column submatrices each corresponding to one group of variables can substantially boost performance. In this note, we study the aspect of standardization for the special case of categorical predictors in detail. The main result is that orthonormalization is not required; column-wise scaling of the design matrix followed by re-scaling and centering of the coefficients is shown to have exactly the same effect. Similar reductions can be achieved in the case of interactions. The extension to the so-called sparse group lasso, which additionally promotes within-group sparsity, is considered as well. The importance of proper standardization is illustrated via extensive simulations.


title: Morphology of AGN Emission Line Regions in SDSS-IV MaNGA Survey

abstract: Extended narrow-line regions (NLRs) around active galactic nuclei (AGN) are shaped by the distribution of gas in the host galaxy and by the geometry of the circumnuclear obscuration, and thus they can be used to test the AGN unification model. In this work, we quantify the morphologies of the narrow-line regions in 308 nearby AGNs ($z=0-0.14$, \lbol $\sim 10^{42.4-44.1}$ \erg{}) from the MaNGA survey. Based on the narrow-line region maps, we find that a large fraction (81\%) of these AGN have bi-conical NLR morphology. The distribution of their measured opening angles suggests that the intrinsic opening angles of the ionization cones has a mean value of 85--98$^\circ$ with a finite spread of 39-44$^\circ$ (1-$\sigma$). Our inferred opening angle distribution implies a number ratio of type I to type II AGN of 1:1.6--2.3, consistent with other measurements of the type I / type II ratio at low AGN luminosities. Combining these measurements with the WISE photometry data, we find that redder mid-IR color (lower effective temperature of dust) corresponds to stronger and narrower photo-ionized bicones. This relation is in agreement with the unification model that suggests that the bi-conical narrow-line regions are shaped by a toroidal dusty structure within a few pc from the AGN. Furthermore, we find a significant alignment between the minor axis of host galaxy disks and AGN ionization cones. Together, these findings suggest that obscuration on both circumnuclear ($\sim $pc) and galactic ($\sim$ kpc) scales are important in shaping and orienting the AGN narrow-line regions.


title: Sensitivity Analysis for Rare Events based on R\'enyi Divergence

abstract: Rare events play a key role in many applications and numerous algorithms have been proposed for estimating the probability of a rare event. However, relatively little is known on how to quantify the sensitivity of the probability with respect to model parameters. In this paper, instead of the direct statistical estimation of rare event sensitivities, we develop novel and general uncertainty quantification and sensitivity bounds which are not tied to specific rare event simulation methods and which apply to families of rare events. Our method is based on a recently derived variational representation for the family of R\'enyi divergences in terms of risk sensitive functionals associated with the rare events under consideration. Based on the derived bounds, we propose new sensitivity indices for rare events and relate them to the moment generating function of the score function. The bounds scale in such a way that we additionally develop sensitivity indices for large deviation rate functions.


title: C/2016 R2 (PANSTARRS): A comet rich in CO and depleted in HCN

abstract: We observed comet C/2016 R2 (PANSTARRS) with the ARO 10-m SMT, and report the first detection of CO emission from this comet with amounts high enough to be the primary driver of activity. We obtained spectra and maps of the CO J=2-1 rotational line at 230 GHz between 2017 December and 2018 January. We calculated an average production rate of Q(CO)=(4.6+/-0.4)x10$^{28}$ mol s$^{-1}$ at r ~2.9 au and delta ~2.1 au. The CO line is thin FWHM ~ 0.8 km s$^{-1}$ with a slight blue-shift ~ -0.1 km s$^{-1}$ from the ephemeris velocity, and we derive a gas expansion velocity of V$_{exp}$ = 0.50+/-0.15 km s$^{-1}$. This comet produced approximately half the CO that comet C/1995 O1 (Hale-Bopp) did at 3 au. If CO production scales with nucleus surface area, then the radius need not exceed ~15 km. The spectra and mapping data are consistent with CO arising from a combination of a sunward-side active area and an isotropic source. For HCN, we calculated a 3-sigma upper limit production rate of Q(HCN) < 8x10$^{24}$ molecules s$^{-1}$, which corresponds to an extraordinarily high abundance ratio limit of Q(CO)/Q(HCN) > 5000. We inferred a production rate of molecular nitrogen of Q(N$_2$) ~2.8x10$^{27}$ molecules s$^{-1}$ using our CO data and the reported N$_2$/CO column density ratio (Cochran & MacKay 2018a,b). The comet does not show the typical nitrogen depletion seen in comets, and the CO-rich, N$_2$-rich and HCN-depleted values are consistent with formation in a cold environment of T < 50 K that may have provided significant N$_2$ shielding.


title: Luzin-type properties and the difference quotient of a real function

abstract: We introduce the notion of the difference quotient set of a real valued function $f$ on a set $E\subset[0,1]$, and compare this set to the range of $f$ on $E$. We discuss the measure theoretic properties of both the range and the difference quotient set of $f$ over $E$ under different assumptions on $f$ and $E$.