title: Sup-norms of eigenfunctions in the level aspect for compact arithmetic surfaces

abstract: Let $D$ be an indefinite quaternion division algebra over $\mathbb{Q}$. We approach the problem of bounding the sup-norms of automorphic forms $\phi$ on $D^\times(\mathbb{A})$ that belong to irreducible automorphic representations and transform via characters of unit groups of orders of $D$. We obtain a non-trivial upper bound for $\|\phi\|_\infty$ in the level aspect that is valid for arbitrary orders. This generalizes and strengthens previously known upper bounds for $\|\phi\|_\infty$ in the setting of newforms for Eichler orders. In the special case when the index of the order in a maximal order is a squarefull integer $N$, our result specializes to $\|\phi\|_\infty \ll_{\pi_\infty, \epsilon} N^{1/3 + \epsilon} \|\phi\|_2$. A key application of our result is to automorphic forms $\phi$ which correspond at the ramified primes to either minimal vectors (in the sense of Hu-Nelson-Saha), or $p$-adic microlocal lifts (in the sense of Nelson). For such forms, our bound specializes to $\| \phi\|_\infty \ll_{\epsilon} C^{\frac16 + \epsilon}\|\phi\|_2$ where $C$ is the conductor of the representation $\pi$ generated by $\phi$. This improves upon the previously known local bound $\|\phi\|_\infty \ll_{\lambda, \epsilon} C^{\frac14 + \epsilon}\|\phi\|_2$ in these cases.


title: Schwarz reflections and the Tricorn

abstract: We continue our exploration of the family $\mathcal{S}$ of Schwarz reflection maps with respect to the cardioid and a circle which was initiated in our earlier work. We prove that there is a natural combinatorial bijection between the geometrically finite maps of this family and those of the basilica limb of the Tricorn, which is the connectedness locus of quadratic anti-holomorphic polynomials. We also show that every geometrically finite map in $\mathcal{S}$ arises as a conformal mating of a unique geometrically finite quadratic anti-holomorphic polynomial and a reflection map arising from the ideal triangle group. We then follow up with a combinatorial mating description for the periodically repelling maps in $\mathcal{S}$. Finally, we show that the locally connected topological model of the connectedness locus of $\mathcal{S}$ is naturally homeomorphic to such a model of the basilica limb of the Tricorn.


title: Optimal Sensor and Actuator Selection using Balanced Model Reduction

abstract: Optimal sensor and actuator selection is a central challenge in high-dimensional estimation and control. Nearly all subsequent control decisions are affected by these sensor/actuator locations, and optimal placement amounts to an intractable brute-force search among the combinatorial possibilities. In this work, we exploit balanced model reduction and greedy optimization to efficiently determine sensor and actuator selections that optimize observability and controllability. In particular, we determine locations that optimize scalar measures of observability and controllability via greedy matrix QR pivoting on the dominant modes of the direct and adjoint balancing transformations. Pivoting runtime scales linearly with the state dimension, making this method tractable for high-dimensional systems. The results are demonstrated on the linearized Ginzburg-Landau system, for which our algorithm approximates known optimal placements computed using costly gradient descent methods.


title: New Galactic multi-mode Cepheids from the ASAS-SN Survey

abstract: A systematic search for multi-mode Cepheids using the database of the ASAS-SN survey has led to the detection of thirteen new double-mode and two triple-mode Cepheids in the Galactic disk. These discoveries have increased the number of Galactic disk multi-mode Cepheids by 33%. One of the new triple-mode variables pulsates simultaneously in the fundamental and in the first and the second radial overtone modes and the other in the first three radial overtone modes. Overtone triple-mode Cepheids were identified only in the Galactic bulge and in the Large and Small Magellanic Clouds previously.


title: Direct measurement of the Higgs self-coupling in $e^+e^- \to ZH$

abstract: A new method to measure the trilinear Higgs self-coupling $\lambda_3^{}$ in a single Higgs production process is proposed. Time-reversal-odd (T-odd) asymmetries in the process $e^+_{}e^-_{} \to ZH$, $Z \to f\bar{f}$ are computed from the absorptive part of the electroweak one-loop amplitude. Since the T-odd asymmetries measure the tree-level $t$-channel $ZH \to ZH$ scattering, they can be direct probes of $\lambda_3^{}$. The proposed method is quite challenging; a relatively large statistics and polarized $e^+_{}e^-_{}$ beams are demanded. However, this is probably the only approach to directly measure $\lambda_3^{}$ in $e^+_{}e^-_{}$ collisions, when a beam energy above the $ZHH$ production threshold is not available.


title: A Duality-Based Unified Approach to Bayesian Mechanism Design

abstract: We provide a unified view of many recent developments in Bayesian mechanism design, including the black-box reductions of Cai et al. [CDW13b], simple auctions for additive buyers [HN12], and posted-price mechanisms for unit-demand bidders [CHK07]. Additionally, we show that viewing these three previously disjoint lines of work through the same lens leads to new developments as well. First, we provide a duality framework for Bayesian mechanism design, which naturally accommodates multiple agents and arbitrary objectives/feasibility constraints. Using this, we prove that either a posted-price mechanism or the Vickrey-Clarke-Groves auction with per-bidder entry fees achieves a constant-factor of the optimal revenue achievable by a Bayesian Incentive Compatible mechanism whenever buyers are unit-demand or additive, unifying previous breakthroughs of Chawla et al. [CHMS10] and Yao [Yao15], and improving both approximation ratios (from 30 to 24 and 69 to 8, respectively). Finally, we show that this view also leads to improved structural characterizations in the Cai et al. framework.


title: Joint analysis of two-way laser range and one-way frequency measurements for the gravitational redshift experiment with the RadioAstron spacecraft

abstract: We consider the problem of joint analysis of two-way laser range and one-way frequency measurements in high-precision tests of general relativity with spacecrafts. Of main interest to such tests is the accuracy of the computed values of the one-way frequency observables. We identify the principal sources of error in these observables to be the errors in the modeled corrections due to various `small' effects, such as that of the troposphere, the error in the reflection time of the laser pulse from the spacecraft, and the error of fitting the spacecraft trajectory to the laser data. We suggest ways to evaluating these errors.


title: A Search for Optical Emission from Binary-Black-Hole Merger GW170814 with the Dark Energy Camera

abstract: Binary black hole (BBH) mergers found by the LIGO and Virgo detectors are of immense scientific interest to the astrophysics community, but are considered unlikely to be sources of electromagnetic emission. To test whether they have rapidly fading optical counterparts, we used the Dark Energy Camera to perform an $i$-band search for the BBH merger GW170814, the first gravitational wave detected by three interferometers. The 87-deg$^2$ localization region (at 90\% confidence) centered in the Dark Energy Survey (DES) footprint enabled us to image 86\% of the probable sky area to a depth of $i\sim 23$ mag and provide the most comprehensive dataset to search for EM emission from BBH mergers. To identify candidates, we perform difference imaging with our search images and with templates from pre-existing DES images. The analysis strategy and selection requirements were designed to remove supernovae and to identify transients that decline in the first two epochs. We find two candidates, each of which is spatially coincident with a star or a high-redshift galaxy in the DES catalogs, and they are thus unlikely to be associated with GW170814. Our search finds no candidates associated with GW170814, disfavoring rapidly declining optical emission from BBH mergers brighter than $i\sim 23$ mag ($L_{\rm optical} \sim 5\times10^{41}$ erg/s) 1-2 days after coalescence. In terms of GW sky map coverage, this is the most complete search for optical counterparts to BBH mergers to date


title: Renormalon-motivated evaluation of QCD observables

abstract: A method of evaluation of spacelike QCD observables ${\cal D}(Q^2)$ is developed, motivated by the renormalon structure of these quantities. A related auxiliary quantity ${\widetilde {\cal D}}(Q^2)$ is introduced, which is renomalization scale independent only at the one-loop level, and a large-$\beta_0$-type renormalon motivated ansatz is made for the Borel transform of this quantity. This leads to a `dressed' Borel transform of the considered observable ${\cal D}(Q^2)$. From there, a Neubert-type distribution is obtained for the observable. The described method is then applied to the massless Adler function and the related decay ratio of the $\tau$ lepton semihadronic decays. Comparisons are then made with an evaluation method at higher truncated orders, developed in our earlier works, which is a renormalization scale invariant extension of the diagonal Pad\'e approximants.


title: On Tur\'{a}n problems for Cartesian products of graphs

abstract: Let $A,B$ be disjoint sets of sizes $n$ and $m$. Let ${\mathcal Q}$ be a family of quadruples, having $2$ elements from $A$ and $2$ from $B$, such that any subset $S \subseteq A \cup B$ with $|S|=7$, $|S \cap A| \geq 2$ and $|S \cap B| \geq 2$ contains one of the quadruples. We prove that the smallest size of ${\mathcal Q}$ is $(1/16 + O(1/n) + O(1/m)) n^2 m^2$ as $n,m\to\infty$. We also solve asymptotically a more general two-partite Tur\'{a}n problem for quadruples.


title: Astronomy summer camp "Beli Brezi", Bulgaria - building the astronomical community of the future

abstract: Why study astronomy, why teach astronomy? We give answers to these fundamental questions based on our experience with the Astronomical Camp "Beli Brezi" (White Aspens; Kardzhali, Bulgaria). It has been a place for teaching astronomy to high schools kids for nearly half a century. We describe shortly the history of the camp and draw some conclusions based on nearly five decades of experience. Major among them is that the camp has gone further than just distributing astronomical knowledge - while this is an important and worthy task, the main achievement has been the cultivation of critical thinking among the pupils and we think that that is the main motivation to give positively reassuring answers the questions we asked at the beginning.


title: On the lifespan of three-dimensional gravity water waves with vorticity

abstract: We prove a long-term regularity result for three-dimensional gravity water waves with small initial data but nonzero initial vorticity. We consider solutions whose vorticity vanishes on the free boundary and use this to derive a system for the evolution of the free boundary which reduces to the Zakharov/Craig-Sulem formulation in the irrotational case. We are able to continue the solution until a time determined by the size of the initial vorticity in such a way that if the vorticity is zero, one recovers a lifespan $T\sim \epsilon^{-N}$ where $N$ can be taken arbitrarily large if the initial data is taken to be arbitrarily smooth.


title: Detect-to-Retrieve: Efficient Regional Aggregation for Image Search

abstract: Retrieving object instances among cluttered scenes efficiently requires compact yet comprehensive regional image representations. Intuitively, object semantics can help build the index that focuses on the most relevant regions. However, due to the lack of bounding-box datasets for objects of interest among retrieval benchmarks, most recent work on regional representations has focused on either uniform or class-agnostic region selection. In this paper, we first fill the void by providing a new dataset of landmark bounding boxes, based on the Google Landmarks dataset, that includes $86k$ images with manually curated boxes from $15k$ unique landmarks. Then, we demonstrate how a trained landmark detector, using our new dataset, can be leveraged to index image regions and improve retrieval accuracy while being much more efficient than existing regional methods. In addition, we introduce a novel regional aggregated selective match kernel (R-ASMK) to effectively combine information from detected regions into an improved holistic image representation. R-ASMK boosts image retrieval accuracy substantially with no dimensionality increase, while even outperforming systems that index image regions independently. Our complete image retrieval system improves upon the previous state-of-the-art by significant margins on the Revisited Oxford and Paris datasets. Code and data available at the project webpage: https://github.com/tensorflow/models/tree/master/research/delf.


title: Constraints on New Scalar and Vector Mediators from LHC Dijet Searches

abstract: We present a reanalysis of the latest results from CMS dijet searches for an integrated luminosity of $36\,$fb$^{-1}$ together with preliminary results for $78\,$fb$^{-1}$ in the framework of simplified models for dark matter interacting with quarks through the exchange of a scalar, pseudoscalar, vector or pseudovector mediator particle. Within the same framework, we also project the sensitivity of dijet searches in future LHC runs and study how well the parameters of a simplified model could be reconstructed in case of a future discovery at the high luminosity (HL) LHC. Finally, we explore the possibility of discriminating different mediator scenarios by extending the sensitivity of dijet searches for simplified models through the use of angular information. It is the first time that these studies are performed systematically for the case of spin 0 mediators. Among other results we find: 1) no evidence for a dijet signal in the simplified model framework; 2) improvements due to an increased luminosity at the HL-LHC are significant, but mostly for heavy mediators, where dijet searches are limited by statistical, rather than systematical uncertainties; 3) Information on the angular separation of dijets could be used at the HL-LHC to discriminate different mediator scenarios.


title: Comment on "Hubble flow variations as a test for inhomogeneous cosmology"

abstract: Saulder et al [2019, A&A, 622, A83; arXiv:1811.11976] have performed a novel observational test of the local expansion of the Universe for the standard cosmology as compared to an alternative model with differential cosmic expansion. Their analysis employs mock galaxy samples from the Millennium Simulation, a Newtonian $N$-body simulation on a $\Lambda$CDM background. For the differential expansion case the simulation has been deformed in an attempt to incorporate features of a particular inhomogeneous cosmology: the timescape model. It is shown that key geometrical features of the timescape cosmology have been omitted in this rescaling. Consequently, the differential expansion model tested by Saulder et al (2019) cannot be considered to approximate the timescape cosmology.


title: A Family of Projective Representations of the Thompson Group and Lifting Problems

abstract: The Thompson group F has a natural unitary representation on $H=L^2[0,1]$. With some projections, we construct a family of projective unitary representations on a Fermionic Fock space associated with $H$. It comes from the representation of the associated CAR algebra. After $H^2(F;S^1)$ is obtained, we mainly study whether any of these projective unitary representations can be lifted to an ordinary one. We will discuss the lifting problem of these projective representations.


title: Two-Step Many-Objective Optimal Power Flow Based on Knee Point-Driven Evolutionary Algorithm

abstract: To coordinate the economy, security and environment protection in the power system operation, a two-step many-objective optimal power flow (MaOPF) solution method is proposed. In step 1, it is the first time that knee point-driven evolutionary algorithm (KnEA) is introduced to address the MaOPF problem, and thereby the Pareto-optimal solutions can be obtained. In step 2, an integrated decision analysis technique is utilized to provide decision makers with decision supports by combining fuzzy c-means (FCM) clustering and grey relational projection (GRP) method together. In this way, the best compromise solutions (BCSs) that represent decision makers' different, even conflicting, preferences can be automatically determined from the set of Pareto-optimal solutions. The primary contribution of the proposal is the innovative application of many-objective optimization together with decision analysis for addressing MaOPF problems. Through examining the two-step method via the IEEE 118-bus system and the real-world Hebei provincial power system, it is verified that our approach is suitable for addressing the MaOPF problem of power systems.


title: $2$-stratifolds with fundamental group $\mathbb{Z}$

abstract: $2$-stratifolds are a generalization of $2$-manifolds that occur as objects in applications such as in TDA. These spaces can be described by an associated bicoloured labelled graph. In previous papers we obtained a classification of 1-connected trivalent $2$-stratifolds. In this paper we classify trivalent $2$-stratifolds that have fundamental group $\mathbb{Z}$. This classification implicitly gives an efficient algorithm that applies to a bicoloured labelled graph to decide wether the associated $2$-stratifold has fundamental group $\mathbb{Z}$.


title: Gravitational solitons on Kasner background revisited: The simplest solitons with physical context

abstract: We revise the one-pair complex poles soliton solutions on a Kasner background. In the literature, these were rejected as solutions with no cosmological interest due to singularities that supposedly show up at space-like infinity. The only accepted solutions of this kind were those with background metric parameter $d=\pm 1$. By computing the scalars $I, J$ we find that there are no scalar singularities at all, for a wide range of the background parameter $d$. This means that there are actually an enormous number of acceptable simple complex-poles solutions, besides the $d=\pm1$ cases. These solutions are interesting, because they are much simpler than the two-pairs complex poles solutions and, consequently, it is easier to draw conclusions and relate physical phenomena to them.


title: A Parallel Double Greedy Algorithm for Submodular Maximization

abstract: We study parallel algorithms for the problem of maximizing a non-negative submodular function. Our main result is an algorithm that achieves a nearly-optimal $1/2 -\epsilon$ approximation using $O(\log(1/\epsilon) / \epsilon)$ parallel rounds of function evaluations. Our algorithm is based on a continuous variant of the double greedy algorithm of Buchbinder et al. that achieves the optimal $1/2$ approximation in the sequential setting. Our algorithm applies more generally to the problem of maximizing a continuous diminishing-returns (DR) function.


title: Debris disks among Kepler solar rotational analog stars

abstract: Observations of circumstellar disks provide a powerful tool for our understanding of planetary systems dynamics. Analogs to the Solar System asteroid belts, debris disks result from the collision of the remaining solid material of the planet formation process. Even if the presence of disk is now reported for hundreds of stars, its detection around stars similar to the Sun is still very sparse. We report the results of a search for debris disks around Kepler stars with surface physical parameters close to solar values, including rotation period, using observations by the Wide-field infrared Survey Explorer (WISE). From the entire sample of Kepler stars, 881 targets were identified with these parameters and only six of them (KIC 1868785, 7267949, 7435796, 10533222, 11352643, and KIC 11666436) show unambiguous infrared excess, for which we determined debris disk physical parameters. Interestingly, the present study reveals traces of debris disks much more massive and brighter than the Solar System zodiacal dust, probably resulting from recent violent collisional events, orbiting stars with ages around the solar values.


title: Improving Semantic Segmentation via Video Propagation and Label Relaxation

abstract: Semantic segmentation requires large amounts of pixel-wise annotations to learn accurate models. In this paper, we present a video prediction-based methodology to scale up training sets by synthesizing new training samples in order to improve the accuracy of semantic segmentation networks. We exploit video prediction models' ability to predict future frames in order to also predict future labels. A joint propagation strategy is also proposed to alleviate mis-alignments in synthesized samples. We demonstrate that training segmentation models on datasets augmented by the synthesized samples leads to significant improvements in accuracy. Furthermore, we introduce a novel boundary label relaxation technique that makes training robust to annotation noise and propagation artifacts along object boundaries. Our proposed methods achieve state-of-the-art mIoUs of 83.5% on Cityscapes and 82.9% on CamVid. Our single model, without model ensembles, achieves 72.8% mIoU on the KITTI semantic segmentation test set, which surpasses the winning entry of the ROB challenge 2018. Our code and videos can be found at https://nv-adlr.github.io/publication/2018-Segmentation.


title: Impact of a XENONnT Signal on LHC Dijet Searches

abstract: It is well-known that dark matter (DM) direct detection experiments and the LHC are complementary, since they probe physical processes occurring at different energy scales. And yet, there are aspects of this complementarity which are still not fully understood, or exploited. For example, what is the impact that the discovery of DM at XENONnT would have on present and future searches for DM in LHC final states involving a pair of hadronic jets? In this work we investigate the impact of a XENONnT signal on the interpretation of current dijet searches at the LHC, and on the prospects for dijet signal discovery at the High-Luminosity (HL) LHC in the framework of simplified models. Specifically, we focus on a general class of simplified models where DM can have spin 0, 1/2 or 1, and interacts with quarks through the exchange of a scalar, pseudo-scalar, vector, or pseudo-vector mediator. We find that exclusion limits on the mediator's mass and its coupling to quarks from dijet searches at the LHC are significantly affected by a signal at XENONnT, and that $\mathcal{O}(100)$ signal events at XENONnT would drastically narrow the region in the parameter space of simplified models where a dijet signal can be discovered at $5\sigma$ C.L. at the HL-LHC.


title: From Crossing-Free Graphs on Wheel Sets to Embracing Simplices and Polytopes with Few Vertices

abstract: A set $P = H \cup \{w\}$ of $n+1$ points in general position in the plane is called a wheel set if all points but $w$ are extreme. We show that for the purpose of counting crossing-free geometric graphs on such a set $P$, it suffices to know the frequency vector of $P$. While there are roughly $2^n$ distinct order types that correspond to wheel sets, the number of frequency vectors is only about $2^{n/2}$. We give simple formulas in terms of the frequency vector for the number of crossing-free spanning cycles, matchings, triangulations, and many more. Based on that, the corresponding numbers of graphs can be computed efficiently. In particular, we rediscover an already known formula for $w$-embracing triangles spanned by $H$. Also in higher dimensions, wheel sets turn out to be a suitable model to approach the problem of computing the simplicial depth of a point $w$ in a set $H$, i.e., the number of $w$-embracing simplices. While our previous arguments in the plane do not generalize easily, we show how to use similar ideas in $\mathbb{R}^d$ for any fixed $d$. The result is an $O(n^{d-1})$ time algorithm for computing the simplicial depth of a point $w$ in a set $H$ of $n$ points, improving on the previously best bound of $O(n^d\log n)$. Based on our result about simplicial depth, we can compute the number of facets of the convex hull of $n=d+k$ points in general position in $\mathbb{R}^d$ in time $O(n^{\max\{\omega,k-2\}})$ where $\omega \approx 2.373$, even though the asymptotic number of facets may be as large as $n^k$.


title: Gravitational Thermodynamics of Causal Diamonds in (A)dS

abstract: The static patch of de Sitter spacetime and the Rindler wedge of Minkowski spacetime are causal diamonds admitting a true Killing field, and they behave as thermodynamic equilibrium states under gravitational perturbations. We explore the extension of this gravitational thermodynamics to all causal diamonds in maximally symmetric spacetimes. Although such diamonds generally admit only a conformal Killing vector, that seems in all respects to be sufficient. We establish a Smarr formula for such diamonds and a "first law" for variations to nearby solutions. The latter relates the variations of the bounding area, spatial volume of the maximal slice, cosmological constant, and matter Hamiltonian. The total Hamiltonian is the generator of evolution along the conformal Killing vector that preserves the diamond. To interpret the first law as a thermodynamic relation, it appears necessary to attribute a negative temperature to the diamond, as has been previously suggested for the special case of the static patch of de Sitter spacetime. With quantum corrections included, for small diamonds we recover the "entanglement equilibrium" result that the generalized entropy is stationary at the maximally symmetric vacuum at fixed volume, and we reformulate this as the stationarity of free conformal energy with the volume not fixed.


title: A Study of Data Store-based Home Automation

abstract: Home automation platforms provide a new level of convenience by enabling consumers to automate various aspects of physical objects in their homes. While the convenience is beneficial, security flaws in the platforms or integrated third-party products can have serious consequences for the integrity of a user's physical environment. In this paper we perform a systematic security evaluation of two popular smart home platforms, Google's Nest platform and Philips Hue, that implement home automation "routines" (i.e., trigger-action programs involving apps and devices) via manipulation of state variables in a centralized data store. Our semi-automated analysis examines, among other things, platform access control enforcement, the rigor of non-system enforcement procedures, and the potential for misuse of routines. This analysis results in ten key findings with serious security implications. For instance, we demonstrate the potential for the misuse of smart home routines in the Nest platform to perform a lateral privilege escalation, illustrate how Nest's product review system is ineffective at preventing multiple stages of this attack that it examines, and demonstrate how emerging platforms may fail to provide even bare-minimum security by allowing apps to arbitrarily add/remove other apps from the user's smart home. Our findings draw attention to the unique security challenges of platforms that execute routines via centralized data stores and highlight the importance of enforcing security by design in emerging home automation platforms.


title: Monocular Total Capture: Posing Face, Body, and Hands in the Wild

abstract: We present the first method to capture the 3D total motion of a target person from a monocular view input. Given an image or a monocular video, our method reconstructs the motion from body, face, and fingers represented by a 3D deformable mesh model. We use an efficient representation called 3D Part Orientation Fields (POFs), to encode the 3D orientations of all body parts in the common 2D image space. POFs are predicted by a Fully Convolutional Network (FCN), along with the joint confidence maps. To train our network, we collect a new 3D human motion dataset capturing diverse total body motion of 40 subjects in a multiview system. We leverage a 3D deformable human model to reconstruct total body pose from the CNN outputs by exploiting the pose and shape prior in the model. We also present a texture-based tracking method to obtain temporally coherent motion capture output. We perform thorough quantitative evaluations including comparison with the existing body-specific and hand-specific methods, and performance analysis on camera viewpoint and human pose changes. Finally, we demonstrate the results of our total body motion capture on various challenging in-the-wild videos. Our code and newly collected human motion dataset will be publicly shared.


title: Dark matter stability and Dirac neutrinos using only Standard Model symmetries

abstract: We provide a generic framework to obtain stable dark matter along with naturally small Dirac neutrino masses generated at the loop level. This is achieved through the spontaneous breaking of the global $U(1)_{B-L}$ symmetry already present in Standard Model. The $U(1)_{B-L}$ symmetry is broken down to a residual even $\mathcal{Z}_n$; $n \geq 4$ subgroup. The residual $\mathcal{Z}_n$ symmetry simultaneously guarantees dark matter stability and protects the Dirac nature of neutrinos. The $U(1)_{B-L}$ symmetry in our setup is anomaly free and can also be gauged in a straightforward way. Finally, we present an explicit example using our framework to show the idea in action.


title: AutoFocus: Efficient Multi-Scale Inference

abstract: This paper describes AutoFocus, an efficient multi-scale inference algorithm for deep-learning based object detectors. Instead of processing an entire image pyramid, AutoFocus adopts a coarse to fine approach and only processes regions which are likely to contain small objects at finer scales. This is achieved by predicting category agnostic segmentation maps for small objects at coarser scales, called FocusPixels. FocusPixels can be predicted with high recall, and in many cases, they only cover a small fraction of the entire image. To make efficient use of FocusPixels, an algorithm is proposed which generates compact rectangular FocusChips which enclose FocusPixels. The detector is only applied inside FocusChips, which reduces computation while processing finer scales. Different types of error can arise when detections from FocusChips of multiple scales are combined, hence techniques to correct them are proposed. AutoFocus obtains an mAP of 47.9% (68.3% at 50% overlap) on the COCO test-dev set while processing 6.4 images per second on a Titan X (Pascal) GPU. This is 2.5X faster than our multi-scale baseline detector and matches its mAP. The number of pixels processed in the pyramid can be reduced by 5X with a 1% drop in mAP. AutoFocus obtains more than 10% mAP gain compared to RetinaNet but runs at the same speed with the same ResNet-101 backbone.


title: Learning 3D Human Dynamics from Video

abstract: From an image of a person in action, we can easily guess the 3D motion of the person in the immediate past and future. This is because we have a mental model of 3D human dynamics that we have acquired from observing visual sequences of humans in motion. We present a framework that can similarly learn a representation of 3D dynamics of humans from video via a simple but effective temporal encoding of image features. At test time, from video, the learned temporal representation give rise to smooth 3D mesh predictions. From a single image, our model can recover the current 3D mesh as well as its 3D past and future motion. Our approach is designed so it can learn from videos with 2D pose annotations in a semi-supervised manner. Though annotated data is always limited, there are millions of videos uploaded daily on the Internet. In this work, we harvest this Internet-scale source of unlabeled data by training our model on unlabeled video with pseudo-ground truth 2D pose obtained from an off-the-shelf 2D pose detector. Our experiments show that adding more videos with pseudo-ground truth 2D pose monotonically improves 3D prediction performance. We evaluate our model, Human Mesh and Motion Recovery (HMMR), on the recent challenging dataset of 3D Poses in the Wild and obtain state-of-the-art performance on the 3D prediction task without any fine-tuning. The project website with video, code, and data can be found at https://akanazawa.github.io/human_dynamics/.


title: Diameter Spanners, Eccentricity Spanners, and Approximating Extremal Distances

abstract: The diameter of a graph is one if its most important parameters, being used in many real-word applications. In particular, the diameter dictates how fast information can spread throughout data and communication networks. Thus, it is a natural question to ask how much can we sparsify a graph and still guarantee that its diameter remains preserved within an approximation $t$. This property is captured by the notion of extremal-distance spanners. Given a graph $G=(V,E)$, a subgraph $H=(V,E_H)$ is defined to be a $t$-diameter spanner if the diameter of $H$ is at most $t$ times the diameter of $G$. We show that for any $n$-vertex and $m$-edges directed graph $G$, we can compute a sparse subgraph $H$ that is a $(1.5)$-diameter spanner of $G$, such that $H$ contains at most $\tilde O(n^{1.5})$ edges. We also show that the stretch factor cannot be improved to $(1.5-\epsilon)$. For a graph whose diameter is bounded by some constant, we show the existence of $\frac{5}{3}$-diameter spanner that contains at most $\tilde O(n^\frac{4}{3})$ edges. We also show that this bound is tight. Additionally, we present other types of extremal-distance spanners, such as $2$-eccentricity spanners and $2$-radius spanners, both contain only $\tilde O(n)$ edges and are computable in $\tilde O(m)$ time. Finally, we study extremal-distance spanners in the dynamic and fault-tolerant settings. An interesting implication of our work is the first $\tilde O(m)$-time algorithm for computing $2$-approximation of vertex eccentricities in general directed weighted graphs. Backurs et al. [STOC 2018] gave an $\tilde O(m\sqrt{n})$ time algorithm for this problem, and also showed that no $O(n^{2-o(1)})$ time algorithm can achieve an approximation factor better than $2$ for graph eccentricities, unless SETH fails; this shows that our approximation factor is essentially tight.


title: Left-Right SU(4) Vector Leptoquark Model for Flavor Anomalies

abstract: Building on our recent proposal to explain the experimental hints of new physics in $B$ meson decays within the framework of Pati-Salam quark-lepton unification, through the interactions of the $(3,1)_{2/3}$ vector leptoquark, we construct a realistic model of this type based on the gauge group ${\rm SU}(4)_L \times {\rm SU}(4)_R \times {\rm SU}(2)_L \times {\rm U}(1)'$ and consistent with all experimental constraints. The key feature of the model is that ${\rm SU}(4)_R$ is broken at a high scale, which suppresses right-handed lepton flavor changing currents at the low scale and evades the stringent bounds from searches for lepton flavor violation. The mass of the leptoquark can be as low as $10 \ {\rm TeV}$ without the need to introduce mixing of quarks or leptons with new vector-like fermions. We provide a comprehensive list of model-independent bounds from low energy processes on the couplings in the effective Hamiltonian that arises from generic leptoquark interactions, and then apply these to the model presented here. We discuss various meson decay channels that can be used to probe the model and we investigate the prospects for discovering the new gauge boson at future colliders.


title: Chandra centres for COSMOS X-ray galaxy groups: Differences in stellar properties between central dominant and offset brightest group galaxies

abstract: We present the results of a search for galaxy clusters and groups in the $\sim2$ square degree of the COSMOS field using all available X-ray observations from the XMM-Newton and Chandra observatories. We reach an X-ray flux limit of $3\times10^{-16}\;ergs\;cm^{-2}\;s^{-1}$ in 0.5--2 keV range, and identify 247 X-ray groups with $M_{200c}=8\times10^{12}-3\times10^{14}\;M_{\odot}$ at a redshift range of $0.08\leq z<1.53$, using the multiband photometric redshift and the master spectroscopic redshift catalogues of the COSMOS. The X-ray centres of groups are determined using high-resolution Chandra imaging. We investigate the relations between the offset of the brightest group galaxies (BGGs) from halo X-ray centre and group properties and compare with predictions from semi-analytic models and hydrodynamical simulations. We find that BGG offset decreases with both increasing halo mass and decreasing redshift with no strong dependence on the X-ray flux and SNR. We show that the BGG offset decreases as a function of increasing magnitude gap with no considerable redshift dependent trend. The stellar mass of BGGs in observations extends over a wider dynamic range compared to model predictions. At $z<0.5$, the central dominant BGGs become more massive than those with large offsets by up to 0.3dex, in agreement with model prediction. The observed and predicted lognormal scatter in the stellar mass of both low- and large-offset BGGs at fixed halo mass is $\sim0.3$dex.


title: proEQUIB: IDL Library for Plasma Diagnostics and Abundance Analysis

abstract: The emission lines emitted from gaseous nebulae carry valuable information about the physical conditions and chemical abundances of ionized gases in these objects, as well as the interstellar extinction. "proEQUIB" is a library containing several application programming interface (API) functions developed in the Interactive Data Language (IDL), which can be used for plasma diagnostics and abundance analysis of nebular spectra. This IDL library can also be used by the GNU Data Language (GDL), which is a free and open-source IDL compiler. This package includes several API functions to determine physical conditions and chemical abundances from collisionally excited lines (CEL) and recombination lines (RL), derive interstellar extinctions from Balmer lines, and deredden the observed fluxes. This IDL library heavily relies on the IDL Astronomy User's library and the IDL "AtomNeb" library. The API functions of this IDL library can easily be utilized for spatially-resolved studies of ionized gaseous nebulae observed using integral field spectroscopy.


title: qpower2 - a fast and accurate algorithm for the computation of exoplanet transit light curves with the power-2 limb-darkening law

abstract: The power-2 law, $I_{\lambda}(\mu)=1- c (1 - \mu^{\alpha})$, accurately represents the limb-darkening profile for cool stars. It has been implemented in a few transit models to-date using numerical integration but there is as-yet no implementation of the power-2 law in analytic form that is generally available. An algorithm to implement the power-2 law is derived using a combination of an approximation to the required integral and a Taylor expansion of the power-2 law. The accuracy of stellar and planetary radii derived by fitting transit light curves with this approximation is tested using light curves computed by numerical integration of limb-darkening profiles from 3D stellar model atmospheres. Our algorithm (qpower2) is accurate to about 100 ppm for broad-band optical light curves of systems with a star-planet radius ratio p = 0.1. The implementation requires less than 40 lines of python code so can run extremely fast on graphical processing units (GPUs; $\sim$ 1 million models per second for the analysis of 1000 data points). Least-squares fits to simulated light curves show that the star and planet radius are recovered to better than 1% for p < 0.2. The qpower2 algorithm can be used to efficiently and accurately analyse large numbers of high-precision transit light curves using Monte Carlo methods.


title: On the structure of the set of positive maps

abstract: The full description of the set of positive maps $T: \qA \to \cB(\cH)$ ($\qA$ a $C^*$-algebra) is given. The approach is based on the simple prescription for selecting various types of positive maps. This prescription stems from the Grothendieck theory of projective tensor products complemented by the theory of tensor connes. In particular, the origin of non-decomposable maps is clarified.


title: Generating High Fidelity Images with Subscale Pixel Networks and Multidimensional Upscaling

abstract: The unconditional generation of high fidelity images is a longstanding benchmark for testing the performance of image decoders. Autoregressive image models have been able to generate small images unconditionally, but the extension of these methods to large images where fidelity can be more readily assessed has remained an open problem. Among the major challenges are the capacity to encode the vast previous context and the sheer difficulty of learning a distribution that preserves both global semantic coherence and exactness of detail. To address the former challenge, we propose the Subscale Pixel Network (SPN), a conditional decoder architecture that generates an image as a sequence of sub-images of equal size. The SPN compactly captures image-wide spatial dependencies and requires a fraction of the memory and the computation required by other fully autoregressive models. To address the latter challenge, we propose to use Multidimensional Upscaling to grow an image in both size and depth via intermediate stages utilising distinct SPNs. We evaluate SPNs on the unconditional generation of CelebAHQ of size 256 and of ImageNet from size 32 to 256. We achieve state-of-the-art likelihood results in multiple settings, set up new benchmark results in previously unexplored settings and are able to generate very high fidelity large scale samples on the basis of both datasets.


title: Approximating the solution to wave propagation using deep neural networks

abstract: Humans gain an implicit understanding of physical laws through observing and interacting with the world. Endowing an autonomous agent with an understanding of physical laws through experience and observation is seldom practical: we should seek alternatives. Fortunately, many of the laws of behaviour of the physical world can be derived from prior knowledge of dynamical systems, expressed through the use of partial differential equations. In this work, we suggest a neural network capable of understanding a specific physical phenomenon: wave propagation in a two-dimensional medium. We define `understanding' in this context as the ability to predict the future evolution of the spatial patterns of rendered wave amplitude from a relatively small set of initial observations. The inherent complexity of the wave equations -- together with the existence of reflections and interference -- makes the prediction problem non-trivial. A network capable of making approximate predictions also unlocks the opportunity to speed-up numerical simulations for wave propagation. To this aim, we created a novel dataset of simulated wave motion and built a predictive deep neural network comprising of three main blocks: an encoder, a propagator made by 3 LSTMs, and a decoder. Results show reasonable predictions for as long as 80 time steps into the future on a dataset not seen during training. Furthermore, the network is able to generalize to an initial condition that is qualitatively different from those seen during training.


title: Imaging at the mesoscale (LEEM, PEEM)

abstract: *To be published in Springer Handbook of Surface Science (Springer Verlag) [Preprint]* The capability to display images containing chemical, magnetic and structural information and to perform spectroscopy and diffraction from a {\mu}m-sized area makes cathode lens electron microscopy one of the most used and reliable techniques to analyze surfaces at the mesoscale. Thanks to its versatility, LEEM/PEEM systems are currently employed to study model systems in the fields of nanotechnology, nanomagnetism, material science, catalysis, energy storage, thin films and 2D materials. In the following chapter, we will present a brief but complete review of this class of instruments. After an historical digression in the introducing section, we will show first the basic operating principles of a simple setup and then the elements that can be added to improve the performances. Later, two sections will be dedicated to LEEM and PEEM respectively. In both cases, a theoretical discussion on the contrast mechanisms will prelude to a showcase of the operating modes of the instrument, with clear examples that will show the best performances available nowadays. Finally, a brief discussion about the future developments of cathode lens electron microscopy will close the chapter.


title: Radiative AGN feedback on a moving mesh: the impact of the galactic disc and dust physics on outflow properties

abstract: Feedback from accreting supermassive black holes, active galactic nuclei (AGN), is now a cornerstone of galaxy formation models. In this work, we present radiation-hydrodynamic simulations of radiative AGN feedback using the novel Arepo-RT. A central black hole emits radiation at a constant luminosity and drives an outflow via radiation pressure on dust grains. Utilising an isolated NFW halo we validate our setup in the single and multi-scattering regimes, with the simulated shock front propagation in excellent agreement with the expected analytic result. For a spherically symmetric NFW halo, an examination of the simulated outflow properties generated by radiative feedback demonstrates that they are lower than typically observed at a fixed AGN luminosity, regardless of the collimation of the radiation. We then explore the impact of a central disc galaxy and the assumed dust model on the outflow properties. The contraction of the halo during the galaxy's formation and modelling the production of dust grains results in a factor $100$ increase in the halo's optical depth. Radiation is then able to couple momentum more efficiently to the gas, driving a stronger shock and producing a mass-loaded $\sim10^{3}\,\mathrm{M}_{\odot}\,\mathrm{yr}^{-1}$ outflow with a velocity of $\sim2000\,\mathrm{km}\,\mathrm{s}^{-1}$, in agreement with observations. However, the inclusion of dust destruction mechanisms, like thermal sputtering, leads to the rapid destruction of dust grains within the outflow, reducing its properties below typically observed values. We conclude that radiative AGN feedback can drive outflows, but a thorough numerical and physical treatment is required to assess its true impact.


title: The active flux scheme on Cartesian grids and its low Mach number limit

abstract: Finite volume schemes for hyperbolic conservation laws require a numerical intercell flux. In one spatial dimension the numerical flux can be successfully obtained by solving (exactly or approximately) Riemann problems that are introduced at cell interfaces. This is more challenging in multiple spatial dimensions. The active flux scheme is a finite volume scheme that considers continuous reconstructions instead. The intercell flux is obtained using additional degrees of freedom distributed along the cell boundary. For their time evolution an exact evolution operator is employed, which naturally ensures the correct direction of information propagation and provides stability. This paper presents an implementation of active flux for the acoustic equations on two-dimensional Cartesian grids and demonstrates its ability to simulate discontinuous solutions with an explicit time stepping in a stable manner. Additionally, it is shown that the active flux scheme for linear acoustics is low Mach number compliant without the need for any fix.


title: Foliated Field Theory and String-Membrane-Net Condensation Picture of Fracton Order

abstract: Foliated fracton order is a qualitatively new kind of phase of matter. It is similar to topological order, but with the fundamental difference that a layered structure, referred to as a foliation, plays an essential role and determines the mobility restrictions of the topological excitations. In this work, we introduce a new kind of field theory to describe these phases: a foliated field theory. We also introduce a new lattice model and string-membrane-net condensation picture of these phases, which is analogous to the string-net condensation picture of topological order.


title: Efficacy of early stellar feedback in low gas surface density environments

abstract: We present a suite of high resolution radiation hydrodynamic simulations of a small patch ($1 \ {\rm kpc}^2$) of the inter-stellar medium (ISM) performed with Arepo-RT, with the aim to quantify the efficacy of various feedback processes like supernovae explosions (SNe), photoheating and radiation pressure in low gas surface density galaxies ($\Sigma_{\rm gas} \simeq 10 \ {\rm M}_\odot \ {\rm pc}^{-2}$). We show that radiation fields decrease the star formation rate and therefore the total stellar mass formed by a factor of $\sim 2$. This increases the gas depletion timescale and brings the simulated Kennicutt-Schmidt relation closer to the observational estimates. Radiation feedback coupled with SNe is more efficient at driving outflows with the mass and energy loading increasing by a factor of $\sim 10$. This increase is mainly driven by the additional entrainment of medium density ($10^{-2} \leq n< 1 \ {\rm cm}^{-3}$), warm ($300 \ {\rm K}\leq T<8000 \ {\rm K}$) material. Therefore including radiation fields tends to launch colder, denser and higher mass and energy loaded outflows. This is because photoheating of the high density gas around a newly formed star over-pressurises the region, causing it to expand. This reduces the ambient density in which the SNe explode by a factor of $10-100$ which in turn increases their momentum output by a factor of $\sim 1.5-2.5$. Finally, we note that in these low gas surface density environments, radiation fields primarily impact the ISM via photoheating and radiation pressure has only a minimal role in regulating star formation.


title: Star Clusters Across Cosmic Time

abstract: Star clusters stand at the intersection of much of modern astrophysics: the interstellar medium, gravitational dynamics, stellar evolution, and cosmology. Here we review observations and theoretical models for the formation, evolution, and eventual disruption of star clusters. Current literature suggests a picture of this life cycle with several phases: (1) Clusters form in hierarchically-structured, accreting molecular clouds that convert gas into stars at a low rate per dynamical time until feedback disperses the gas. (2) The densest parts of the hierarchy resist gas removal long enough to reach high star formation efficiency, becoming dynamically-relaxed and well-mixed. These remain bound after gas removal. (3) In the first $\sim 100$ Myr after gas removal, clusters disperse moderately fast, through a combination of mass loss and tidal shocks by dense molecular structures in the star-forming environment. (4) After $\sim 100$ Myr, clusters lose mass via two-body relaxation and shocks by giant molecular clouds, processes that preferentially affect low-mass clusters and cause a turnover in the cluster mass function to appear on $\sim 1-10$ Gyr timescales. (5) Even after dispersal, some clusters remain coherent and thus detectable in chemical or action space for multiple galactic orbits. In the next decade a new generation of space- and AO-assisted ground-based telescopes will enable us to test and refine this picture.


title: Large Stark Tuning of InAs/InP Quantum Dots

abstract: InAs/InP quantum dots are excellent sources of telecom single-photon emission and are among the most promising candidates for scalable quantum photonic circuits. However, geometric differences in each quantum dot leads to slightly different emission wavelengths and hinders the possibility of generating multiple identical quantum emitters on the same chip. Stark tuning is an efficient technique to overcome this issue as it can control the emission energy of individual quantum dots through the quantum-confined Stark effect. Realizing this technique in InAs/InP quantum dots has previously been limited to shifts of less than 0.8 meV due to jumps in the emission energy because of additional charges at high electric field intensities. We demonstrate up to 5.1 meV of Stark tuning in the emission wavelength of InAs/InP quantum dots. To eliminate undesirable jumps to charged state, we use a thin oxide insulator to prevent carrier injection from the contacts, thereby significantly improves the tuning range of the Stark effect. Moreover, the single-photon nature and narrow linewidth of the quantum dot emission is preserved under a wide range of applied electric fields. Using photoluminescence intensity measurements and time-resolved lifetime spectroscopy we confirmed that this Stark tuning range is limited by carrier tunneling at high electric fields. This result is an important step toward integrating multiple identical quantum emitters at telecom wavelengths on-a-chip, which is crucial for realizing complex quantum photonic circuits for quantum information processing.

categories: cond-mat.mes-hall physics.app-ph physics.optics quant-ph


title: Oracles for Gauss's law on digital quantum computers

abstract: Formulating a lattice gauge theory using only physical degrees of freedom generically leads to non-local interactions. A local Hamiltonian is desirable for quantum simulation, and this is possible by treating the Hilbert space as a subspace of a much larger Hilbert space in which Gauss's law is not automatic. Digital quantum simulations of this local formulation will wander into unphysical sectors due to errors from Trotterization or from quantum noise. In this work, oracles are constructed that use local Gauss law constraints to projectively distinguish physical and unphysical wave functions in Abelian lattice gauge theories. Such oracles can be used to detect errors that break Gauss's law.


title: Shadow Imaging of Transiting Objects

abstract: We consider the problem of inferring the shape of a transiting object's silhouette from its light curve alone, without assuming a physical model for the object. We model the object as a grid of pixels which transits the star; each pixel has an opacity, ranging from transparent to opaque, which we infer from the light curve. We explore three interesting degeneracies inherent to this problem, in which markedly different transiting shapes can produce identical light curves: (i) the "flip" degeneracy, by which two pixels transiting at the same impact parameter on opposite sides of the star's horizontal midplane generate the same light curve; (ii) the "arc" degeneracy, by which opacity can be redistributed along the semicircular arc of pixels which undergoes ingress or egress at the same time without consequence to the light curve, and (iii) the "stretch" degeneracy, by which a wide shape moving fast can produce the same light curve as a narrow shape moving more slowly. By understanding these degeneracies and adopting some additional assumptions, we are able to numerically recover informative shadow images of transiting objects, and we explore a number of different algorithmic approaches to this problem. We apply our methods to real data, including the TRAPPIST-1c,e,f triple transit and two dips of Boyajian's Star. We provide Python code to calculate the transit light curve of any grid and, conversely, infer the image grid which generates any light curve in a software package accompanying this paper, EightBitTransit.


title: Excited mesons, baryons, glueballs and tetraquarks: Predictions of the Holography Inspired Stringy Hadron model

abstract: In this note we collect and summarize the predictions of the Holography Inspired Stringy Hadron (HISH) model. We list the masses and widths of predicted excited states across the spectrum, based on placing the different hadrons on the non-linear Regge trajectories of a string with massive endpoints. Our predicted states include: (i) Light, heavy-light and heavy-heavy mesons. (ii) Baryons, including charmed, doubly charmed and bottom baryons. (iii) Glueballs, together with a method to disentangle them from flavorless mesons. (iv) Genuine tetraquarks, which are not "molecules" of hadrons, and are characterized by their decay into a baryon and an anti-baryon.


title: Follow-up signals from superradiant instabilities of black hole merger remnants

abstract: Superradiant instabilities can trigger the formation of bosonic clouds around rotating black holes. If the bosonic field growth is sufficiently fast, these clouds could form shortly after a binary black hole merger. Such clouds are continuous sources of gravitational waves whose detection (or lack thereof) can probe the existence of ultralight bosons (such as axion-like particles) and their properties. Motivated by the binary black hole mergers seen by Advanced LIGO so far, we investigate in detail the parameter space that can be probed with continuous gravitational wave signals from ultralight scalar field clouds around black hole merger remnants with particular focus on future ground-based detectors (A+, Voyager and Cosmic Explorer). We also study the impact that the confusion noise from a putative stochastic gravitational-wave background from unresolved sources would have on such searches and we estimate, under different astrophysical priors, the number of binary black-hole merger events that could lead to an observable post-merger signal. Under our most optimistic assumptions, Cosmic Explorer could detect dozens of post-merger signals.


title: Ultra-slow dynamics in a translationally invariant spin model for multiplication and factorization

abstract: We construct a model of short-range interacting Ising spins on a translationally invariant two-dimensional lattice that mimics a reversible circuit that multiplies or factorizes integers, depending on the choice of boundary conditions. We prove that, for open boundary conditions, the model exhibits no finite-temperature phase transition. Yet we find that it displays glassy dynamics with astronomically slow relaxation times, numerically consistent with a double exponential dependence on the inverse temperature. The slowness of the dynamics arises due to errors that occur during thermal annealing that cost little energy but flip an extensive number of spins. We argue that the energy barrier that needs to be overcome in order to heal such defects scales linearly with the correlation length, which diverges exponentially with inverse temperature, thus yielding the double exponential behavior of the relaxation time.


title: A special Calabi-Yau degeneration with trivial monodromy

abstract: A well-known theorem of Kulikov, Persson and Pinkham states that a degeneration of a family of K3-surfaces with trivial monodromy can be completed to a smooth family. We give a simple example that an analogous statement does not hold for Calabi-Yau threefolds.


title: RadioAstron orbit determination and evaluation of its results using correlation of space-VLBI observations

abstract: A crucial part of a space mission for very-long baseline interferometery (VLBI), which is the technique capable of providing the highest resolution images in astronomy, is orbit determination of the mission's space radio telescope(s). In order to successfully detect interference fringes that result from correlation of the signals recorded by a ground-based and a space-borne radio telescope, the propagation delays experienced in the near-Earth space by radio waves emitted by the source and the relativity effects on each telescope's clock need to be evaluated, which requires accurate knowledge of position and velocity of the space radio telescope. In this paper we describe our approach to orbit determination (OD) of the RadioAstron spacecraft of the RadioAstron space-VLBI mission. Determining RadioAstron's orbit is complicated due to several factors: strong solar radiation pressure, a highly eccentric orbit, and frequent orbit perturbations caused by the attitude control system. We show that in order to maintain the OD accuracy required for processing space-VLBI observations at cm-wavelengths it is required to take into account the additional data on thruster firings, reaction wheel rotation rates, and attitude of the spacecraft. We also investigate into using the unique orbit data available only for a space-VLBI spacecraft, i.e. the residual delays and delay rates that result from VLBI data processing, as a means to evaluate the achieved OD accuracy. We present the results of the first experience of OD accuracy evaluation of this kind, using more than 5,000 residual values obtained as a result of space-VLBI observations performed over 7 years of the RadioAstron mission operations.


title: Obliquity Tides May Drive WASP-12b's Rapid Orbital Decay

abstract: Recent analyses have revealed a mystery. The orbital period of the highly inflated hot Jupiter, WASP-12b, is decreasing rapidly. The rate of inspiral, however, is too fast to be explained by either eccentricity tides or equilibrium stellar tides. While dynamical stellar tides are possible, they require a subgiant structure for the star, whereas stellar models point toward a main sequence host. Here, we show that these hitherto irreconcilable observations might be explained by planetary obliquity tides if planet b's spin vector is trapped in a high-obliquity state maintained by a secular spin-orbit resonance with an unseen exterior perturbing planet. We derive constraints on the obliquity ($\epsilon\gtrsim50^{\circ}$), reduced tidal quality factor ($Q^{\prime}\sim10^{6}-10^{7}$), and perturbing planet parameters ($M_{2}\sim10-20M_{\oplus}$, $a_2\lesssim0.04\,{\rm AU}$) required to generate the observed orbital decay. Direct N-body simulations that include tidal and spin dynamics reinforce the plausibility of the scenario. Furthermore, we show that the resonance could have been captured when planet b's obliquity was small, making the proposed sequence of events easy to explain. The hypothetical perturbing planet is within the limits of current radial velocity constraints on the system yet is also detectable. If it exists, it could provide evidence in favor of the in situ formation hypothesis for hot Jupiters.


title: Nontrivial Quantum Cellular Automata in Higher Dimensions

abstract: We construct a three-dimensional quantum cellular automaton (QCA), an automorphism of the local operator algebra on a lattice of qubits, which disentangles the ground state of the Walker-Wang three fermion model. We show that if this QCA can be realized by a quantum circuit of constant depth, then there exists a two-dimensional commuting projector Hamiltonian which realizes the three fermion topological order which is widely believed not to be possible. We conjecture in accordance with this belief that this QCA is not a quantum circuit of constant depth, and we provide two further pieces of evidence to support the conjecture. We show that this QCA maps every local Pauli operator to a local Pauli operator, but is not a Clifford circuit of constant depth. Further, we show that if the three-dimensional QCA can be realized by a quantum circuit of constant depth, then there exists a two-dimensional QCA acting on fermionic degrees of freedom which cannot be realized by a quantum circuit of constant depth; i.e., we prove the existence of a nontrivial QCA in either three or two dimensions. The square of our three-dimensional QCA can be realized by a quantum circuit of constant depth, and this suggests the existence of a $\mathbb{Z}_2$ invariant of a QCA in higher dimensions, totally distinct from the classification by positive rationals (i.e., by one integer index for each prime) in one dimension. In an appendix, unrelated to the main body of this paper, we give a fermionic generalization of a result of Bravyi and Vyalyi on ground states of 2-local commuting Hamiltonians.

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


title: Structure of a protobinary system: an asymmetric circumbinary disk and spiral arms

abstract: We investigate the gas structures around young binary stars by using three-dimensional numerical simulations. Each model exhibits circumstellar disks, spiral arms, and a circumbinary disk with an inner gap or cavity. The circumbinary disk has an asymmetric pattern rotating at an angular velocity of approximately one-fourth of the binary orbit of the moderate-temperature models. Because of this asymmetry, the circumbinary disk has a density bump and a vortex, both of which continue to exist until the end of our calculation. The density bump and vortex are attributed to enhanced angular momentum, which is promoted by the gravitational torque of the stars. In a hot model ($c \ge 2.0$), the asymmetry rotates considerably more slowly than in the moderate-temperature models. The cold models ($c \le 0.02$) exhibit eccentric circumbinary disks, the precession of which is approximated by a secular motion of the ballistic particles. The asymmetry in the circumbinary disk does not depend on the mass ratio, but it becomes less clear as the specific angular momentum of the infalling envelope increases. The relative accretion rate onto the stars is sensitive to the angular momentum of the infalling envelope. For envelopes with constant angular momentum, the secondary tends to have a higher accretion rate than the primary, except in very low angular momentum cases. For envelopes with a constant angular velocity, the primary has a higher accretion rate than the secondary because gas with low specific angular momentum falls along the polar directions.


title: A response function approach for rapid far-field tsunami forecasting

abstract: Predicting tsunami impacts at remote coasts largely relies on tsunami en-route measurements in an open ocean. In this work, these measurements are used to generate instant tsunami predictions in deep water and near the coast. The predictions are generated as a response or a combination of responses to one or more tsunameters, with each response obtained as a convolution of real-time tsunameter measurements and a pre-computed Pulse Response Function (PRF). Practical implementation of this method requires tables of PRFs in a 3D parameter space: earthquake location - tsunameter - forecasted site. Examples of hindcasting the 2010 Chilean and the 2011 Tohoku-Oki tsunamis along the US West Coast and beyond demonstrated high accuracy of the suggested technology in application to trans-Pacific seismically-generated tsunamis.


title: Playing Text-Adventure Games with Graph-Based Deep Reinforcement Learning

abstract: Text-based adventure games provide a platform on which to explore reinforcement learning in the context of a combinatorial action space, such as natural language. We present a deep reinforcement learning architecture that represents the game state as a knowledge graph which is learned during exploration. This graph is used to prune the action space, enabling more efficient exploration. The question of which action to take can be reduced to a question-answering task, a form of transfer learning that pre-trains certain parts of our architecture. In experiments using the TextWorld framework, we show that our proposed technique can learn a control policy faster than baseline alternatives. We have also open-sourced our code at https://github.com/rajammanabrolu/KG-DQN.


title: Star Clusters in the Elliptical Galaxy NGC 4589 Hosting a Calcium-rich SN Ib (SN 2005CZ)

abstract: _NGC 4589, a bright E2 merger-remnant galaxy, hosts the peculiar fast and faint calcium-rich Type Ib supernova (SN) SN 2005cz. The progenitor of Ca-rich SNe Ib has been controversial: it could be a) a young massive star with 6-12 M$\odot$ in a binary system, or b) an old low-mass star in a binary system that was kicked out from the galaxy center. Moreover, previous distance estimates for this galaxy have shown a large spread, ranging from 20 Mpc to 60 Mpc. Thus, using archival $Hubble$ $Space$ $Telescope$/ACS $F435W$, $F555W$, and $F814W$ images, we search for star clusters in NGC 4589 in order to help resolve these issues. We find a small population of young star clusters with $25<V\leq27$ ($-7.1<M_V\leq-5.1$) mag and age $< 1$ Gyr in the central region at $R<0.5'$ ($<3.8$ kpc), thus supporting the massive-star progenitor scenario for SN 2005cz. In addition to young star clusters, we also find a large population of old globular clusters. In contrast to previous results in the literature, we find that the color distribution of the globular clusters is clearly bimodal. The turnover (Vega) magnitude in the $V$-band luminosity functions of the blue (metal-poor) globular clusters is determined to be $V_0{(\rm max)}=24.40\pm0.10$ mag. We derive the total number of globular clusters, $N_{\rm GC} =640\pm50$, and the specific frequency, $S_N =1.7\pm0.2$. Adopting a calibration for the metal-poor globular clusters, $M_V({\rm max})=-7.66\pm0.14$ mag, we derive a distance to this galaxy: $(m-M)_0=32.06\pm0.10({\rm ran})\pm0.15({\rm sys})$ ($d=25.8\pm2.2$ Mpc).


title: Probing Neutrino Dirac Mass in Left-Right Symmetric Models at the LHC and Next Generation Colliders

abstract: We assess the sensitivity of the LHC, its high energy upgrade, and a prospective 100 TeV hadronic collider to the Dirac Yukawa coupling of the heavy neutrinos in left-right symmetric models (LRSMs). We focus specifically on the trilepton final state in regions of parameter space yielding prompt decays of the right-handed gauge bosons ($W_R$) and neutrinos ($N_R$). In the minimal LRSM, the Dirac Yukawa couplings are completely fixed in terms of the mass matrices for the heavy and light neutrinos. In this case, the trilepton signal provides a direct probe of the Dirac mass term for a fixed $W_R$ and $N_R$ mass. We find that while it is possible to discover the $W_R$ at the LHC, probing the Dirac Yukawa couplings will require a 100 TeV $pp$ collider. We also show that the observation of the trilepton signal at the LHC would indicate the presence of a non-minimal LRSM scenario.


title: 4D ensembles of percolating center vortices and chains

abstract: In this work, we review a recently proposed measure to compute center-element averages in a mixed ensemble of center vortices and chains with non-Abelian d.o.f. and monopole fusion. When center vortices percolate and monopoles condense, the average is captured by a saddle point and collective modes in a YMH model. In this manner, the L\"uscher term, confining flux tubes with N-ality and confined gluons were accommodated in an ensemble picture.


title: Discovery of a Damped Ly-alpha System in a Low-z Galaxy Group: Possible Evidence for Gas Inflow and Nuclear Star Formation

abstract: We present a low-redshift (z=0.029) Damped Lyman-alpha (DLA) system in the spectrum of a background Quasi-Stellar Object (QSO). The DLA is associated with an interacting galaxy pair within a galaxy group. We detected weak Lyman-alpha emission centered at the absorption trough of the DLA. The emission was likely tracing the neutral HI reservoir around the galaxies in the interacting pair, which scattered the Lyman-alpha generated by star formation within those galaxies. We also found that the interacting pair is enveloped by a large HI cloud with $M(HI)=2\times 10^{10}M_{\odot}$. We discovered blueshifted 21cm HI emission, corresponding to M(HI)=$ 2\times10^{9}~M_{\odot}$, associated with J151225.15+012950.4 - one of the galaxies in the interacting pair. The blueshifted HI was tracing gas flowing into the galaxy from behind and towards us. Gas at similar blueshifted velocities was seen in the QSO sightline thus suggesting the presence of a filamentary structure of the order of 100kpc feeding the galaxy. We estimated a mass inflow rate of $2 M_{\odot}~yr^{-1}$ into the galaxy, which matches the star formation rate estimated from H-alpha emission. It is likely that the inflow of enormous amounts of gas has triggered star formation in this galaxy. The sudden acquisition of cold gas may lead to a starburst in this galaxy like those commonly seen in simulations.


title: Exotic Higgs Decays in Type-II 2HDMs at the LHC and Future 100 TeV Hadron Colliders

abstract: The exotic decay modes of non-Standard Model (SM) Higgses in models with extended Higgs sectors have the potential to serve as powerful search channels to explore the space of Two-Higgs Doublet Models (2HDMs) that cannot be studied effectively using conventional decay channels. Once kinematically allowed, heavy Higgses could decay into pairs of light non-SM Higgses, or a non-SM Higgs and a SM gauge boson, with branching fractions that dominate those of the conventional decay modes to SM particles. In this study, we focus on the prospects of probing exotic decay channels at the LHC and a future 100 TeV \emph{pp} collider in the context of Type-II 2HDMs. We study the three prominent exotic decay channels, A -> HZ, A -> H^+ W and H^+ -> HW, and find that a 100-TeV pp collider can probe the entire region of the Type-II 2HDM parameter space that survives current theoretical and experimental constraints with exotic decay branching fraction > 20%.


title: A Semismooth Predictor Corrector Method for Real-Time Constrained Parametric Optimization with Applications in Model Predictive Control

abstract: Real-time optimization problems are ubiquitous in control and estimation, and are typically parameterized by incoming measurement data and/or operator commands. This paper proposes solving parameterized constrained nonlinear programs using a semismooth predictor-corrector (SSPC) method. Nonlinear complementarity functions are used to reformulate the first order necessary conditions of the optimization problem into a parameterized non-smooth root-finding problem. Starting from an approximate solution, a semismooth Euler-Newton algorithm is proposed for tracking the trajectory of the primal-dual solution as the parameter varies over time. Active set changes are naturally handled by the SSPC method, which only requires the solution of linear systems of equations. The paper establishes conditions under which the solution trajectories of the root-finding problem are well behaved and provides sufficient conditions for ensuring boundedness of the tracking error. Numerical case studies featuring the application of the SSPC method to nonlinear model predictive control are reported and demonstrate the advantages of the proposed method.


title: Helium Suppression in Impulsive Solar Energetic-Particle Events

abstract: We have studied the element abundances and energy spectra of the small "He-poor" impulsive solar energetic-particle (SEP) events, comparing them with other impulsive SEP events with more-normal abundances of He. He-poor events can have abundances as low as He/O ~ 2, while both impulsive and gradual SEP events usually have source abundances of 30 < He/O < 100 with mean values of 50 - 60. He/C ratios are not only low, but often decrease with energy in He-poor events. Abundance enhancement patterns of other elements with atomic numbers 6 < Z < 56, and likely values of their mass-to-charge ratios A/Q, are generally unaltered in He-poor events, as are the probable source-plasma temperatures of 2.5 - 3.2 MK for all impulsive SEP events. One He-poor event is also an example of a rarer C-poor event with C/O = 0.08 +- 0.04, sup-pressed by a factor over 5 from the mean. We discuss suggestions of a possible A/Q threshold during acceleration and of the sluggish ionization of He entering the corona, because of its uniquely high first ionization potential (FIP), but the suppression of He and the decline of He/C with energy is difficult to explain if both He and C are fully ionized with A/Q = 2 as expected at 2.5 - 3.2 MK. Although less dramatic, a possible excess enhancement of Ne in some impulsive SEP events is also considered. Possible causes of the large ~30% spectral and abundance variations in impulsive events are also discussed. However, the physics of the He-poor events remains a mystery.


title: Measuring ISW with next-generation radio surveys

abstract: The late-time integrated Sachs-Wolfe (ISW) signal in the CMB temperature anisotropies is an important probe of dark energy when it can be detected by cross-correlation with large-scale structure surveys. Because of their huge sky area, surveys in the radio are well-suited to ISW detection. We show that 21cm intensity mapping and radio continuum surveys with the SKA in Phase~1 promise a $\sim5\sigma$ detection if we use tomography, with a similar forecast for the precursor EMU survey. In SKA Phase~2, the 21cm galaxy redshift survey and the continuum survey could deliver a $\sim6\sigma$ detection. Our analysis of the radio surveys aims for theoretical accuracy on large scales. Firstly, we include all the effects on the radio surveys from observing on the past lightcone: redshift-space distortions and lensing magnification can have a significant impact on the ISW signal to noise ratio, while Doppler and other relativistic distortions are not significant. Secondly, we use the full information in the observable galaxy angular power spectra $C_\ell(z,z')$, by avoiding the Limber approximation and by including all cross-correlations between redshift bins in the covariance. Without these cross-bin correlations, the ISW signal to noise ratio is biased.


title: Vertex-Connectivity Measures for Node Failure Identification in Boolean Network Tomography

abstract: In this paper we study the node failure identification problem in undirected graphs by means of Boolean Network Tomography. We argue that vertex connectivity plays a central role. We show tight bounds on the maximal identifiability in a particular class of graphs, the Line of Sight networks. We prove slightly weaker bounds on arbitrary networks. Finally we initiate the study of maximal identifiability in random networks. We focus on two models: the classical Erd\H{o}s-R\'enyi model, and that of Random Regular graphs. The framework proposed in the paper allows a probabilistic analysis of the identifiability in random networks giving a tradeoff between the number of monitors to place and the maximal identifiability.


title: A Navigational Approach to Health: Actionable Guidance for Improved Quality of Life

abstract: Lifestyle and environment interacting with our biological machine are primarily responsible for shaping our health and wellbeing. Continuous, multi-modal, and quantitative approaches to understanding and controlling these factors will allow each person to better reach their desired quality of life. A navigational paradigm can help users towards a specified health goal by using constantly captured measurements to feed estimations of how a user's health is continuously changing in order to provide actionable guidance. As various actions are taken by the user, measurements of the resulting effects loop back into the estimation and the next step of guidance. This perpetual cycle of measuring, estimating, guiding, and acting articulates a Personal Health Navigation information and actuation framework. Personal Health Navigation focuses on fulfilling a user's health goals by ensuring minimal deviation from healthy states, rather than treating disease or symptoms after derailment from proper biological function.


title: Adsorption and dissociation of iron phthalocyanine on H/Si(111): Impact of van-der-Waals interactions and perspectives for subsurface doping

abstract: The adsorption of iron phthalocyanine (FePc) on the passivated H/Si(111) surface is explored from first principles. We find that the organic molecule is predominantly physisorbed with a distance to the surface of $2.6 \pm 0.1$ Angstrom, but also exhibits sizable resonance with the underlying substrate. This establishes the present system as interesting mixed covalent-van-der-Waals-bound test case, which we use to compare the impact of different approaches to van-der-Waals interactions. (Spin-polarized) scanning tunneling microscopy (SP STM) images are simulated, selectively accessing different molecular orbitals via the applied bias voltage in the spirit of scanning tunneling spectroscopy. Comparison with experimental STM images reveals very good agreement. We find a significant magnetic contrast exceeding $\pm 1$ Angstrom in the SP STM images for $-2$ and $+1.5$ V. Binding energies of different (transition metal) atoms in the center of the Pc ring are presented, which particularly show that Fe is strongly bound in the molecule (about $9.6$ eV). Finally, we discuss different reactions for subsurface doping by room-temperature FePc deposition and point out two feasible reactions.


title: Overcoming Catastrophic Forgetting by Soft Parameter Pruning

abstract: Catastrophic forgetting is a challenge issue in continual learning when a deep neural network forgets the knowledge acquired from the former task after learning on subsequent tasks. However, existing methods try to find the joint distribution of parameters shared with all tasks. This idea can be questionable because this joint distribution may not present when the number of tasks increase. On the other hand, It also leads to "long-term" memory issue when the network capacity is limited since adding tasks will "eat" the network capacity. In this paper, we proposed a Soft Parameters Pruning (SPP) strategy to reach the trade-off between short-term and long-term profit of a learning model by freeing those parameters less contributing to remember former task domain knowledge to learn future tasks, and preserving memories about previous tasks via those parameters effectively encoding knowledge about tasks at the same time. The SPP also measures the importance of parameters by information entropy in a label free manner. The experiments on several tasks shows SPP model achieved the best performance compared with others state-of-the-art methods. Experiment results also indicate that our method is less sensitive to hyper-parameter and better generalization. Our research suggests that a softer strategy, i.e. approximate optimize or sub-optimal solution, will benefit alleviating the dilemma of memory. The source codes are available at https://github.com/lehaifeng/Learning_by_memory.


title: New dual gravitational charges

abstract: We show that there are a further infinite number of, previously unknown, supertranslation charges. These can be viewed as duals of the known BMS charges corresponding to supertranslations. In Newman-Penrose language, these new supertranslation charges roughly correspond to the imaginary part of the leading term in $\psi_2$. We find these charges by dualising the Barnich-Brandt asymptotic charges and argue that this prescription gives rise to new bona fide charges at null infinity.


title: Paired parton quantum Hall states: a coupled wire construction

abstract: The Pfaffian fractional quantum Hall (FQH) states are incompressible non-Abelian topological fluids present in a half-filled electron Landau level, where there is a balanced population of electrons and holes. They give rise to half-integral quantum Hall plateaus that divide critical transitions between integer quantum Hall (IQH) states. On the other hand, there are Abelian FQH states, such as the Laughlin state, that can be understood using partons, which are fermionic divisions of the electron. In this paper, we propose a new family of incompressible paired parton FQH states at filling $\nu=1/6$ (modulo 1) that emerge from critical transitions between IQH states and Abelian FQH states at filling $\nu=1/3$ (modulo 1). These paired parton states are originated from a half-filled parton Landau level, where there is an equal amount of partons and holes. They generically support Ising-like anyonic quasiparticle excitations and carry non-Abelian Pfaffian topological orders (TO) for partons. We prove the principle existence of these paired parton states using exactly solvable interacting arrays of electronic wires under a magnetic field. Moreover, we establish a new notion of particle-hole (PH) symmetry for partons and relate the PH symmetric parton Pfaffian TO with the gapped symmetric surface TO of a fractional topological insulator in three dimension.


title: Five-dimensional electromagnetism

abstract: The usual pre-metric Maxwell equations of electromagnetism in Minkowski space are extended to five dimensions, and the equations of linear five-dimensional electromagnetic waves are obtained from them by using the simplest electromagnetic constitutive law. The role of the dispersion law for five-dimensional plane waves in dictating the most appropriate extension of the four-dimensional Minkowski space metric to the five-dimensional one is discussed. The fact that the five-dimensional dispersion law includes both massless waves and massive ones is emphasized, as well as the fact that one can obtain the Klein-Gordon equation or the Proca equation by separating the fifth coordinate.


title: Updated CLIC luminosity staging baseline and Higgs coupling prospects

abstract: An updated luminosity staging baseline for CLIC is presented. Assuming accelerator ramp-up and up-time scenarios that are harmonized with those of other potential future colliders, CLIC will deliver 1 ab^-1 at sqrt(s)=380 GeV, 2.5 ab^-1 at sqrt(s)=1.5 TeV, and 5 ab^-1 at sqrt(s)=3 TeV. The complete programme will take 25-30 years. The baseline scenario for luminosity sharing between the two electron beam polarisation states is also discussed. Updated Higgs coupling sensitivities are given for this new luminosity staging baseline.


title: Effect of a laser field in the confinement potential of two electrons in a double quantum dot

abstract: We have studied a system consisted of two coupled quantum dots containing two electrons subjected by a laser field. The effect of the laser is described by the dressed-band approach involving the concept of the conduction/valence effective mass, valid far from resonance. The interaction between the electrons and the quantum dots is described by a phenomenological tridimensional potential, which simulates quantum dots in GaAs heterostructure. In this study we have employed the approach already presented in a previous work [Olavo et al., J. Phys. B: At. Mol. Opt. Phys. {\bf 49}, 145004 (2016)]. We have used a code based on the full interaction configuration method. We have employed as basis set the {\it Cartesian anisotropic Gaussian-type} orbitals which allows one to explore the confining characteristics of a potential due to their flexibility of using different exponents for each direction space. We present an analysis based on the energy levels of the singlet and triplet as function of the confinement parameters.


title: From spatio-temporal data to chronological networks: An application to wildfire analysis

abstract: Network theory has established itself as an appropriate tool for complex systems analysis and pattern recognition. In the context of spatiotemporal data analysis, correlation networks are used in the vast majority of works. However, the Pearson correlation coefficient captures only linear relationships and does not correctly capture recurrent events. This missed information is essential for temporal pattern recognition. In this work, we propose a chronological network construction process that is capable of capturing various events. Similar to the previous methods, we divide the area of study into grid cells and represent them by nodes. In our approach, links are established if two consecutive events occur in two different nodes. Our method is computationally efficient, adaptable to different time windows and can be applied to any spatiotemporal data set. As a proof-of-concept, we evaluated the proposed approach by constructing chronological networks from the MODIS dataset for fire events in the Amazon basin. We explore two data analytic approaches: one static and another temporal. The results show some activity patterns on the fire events and a displacement phenomenon over the year. The validity of the analyses in this application indicates that our data modeling approach is very promising for spatio-temporal data mining.

categories: physics.soc-ph cs.CY physics.data-an physics.geo-ph


title: Rigorous Agent Evaluation: An Adversarial Approach to Uncover Catastrophic Failures

abstract: This paper addresses the problem of evaluating learning systems in safety critical domains such as autonomous driving, where failures can have catastrophic consequences. We focus on two problems: searching for scenarios when learned agents fail and assessing their probability of failure. The standard method for agent evaluation in reinforcement learning, Vanilla Monte Carlo, can miss failures entirely, leading to the deployment of unsafe agents. We demonstrate this is an issue for current agents, where even matching the compute used for training is sometimes insufficient for evaluation. To address this shortcoming, we draw upon the rare event probability estimation literature and propose an adversarial evaluation approach. Our approach focuses evaluation on adversarially chosen situations, while still providing unbiased estimates of failure probabilities. The key difficulty is in identifying these adversarial situations -- since failures are rare there is little signal to drive optimization. To solve this we propose a continuation approach that learns failure modes in related but less robust agents. Our approach also allows reuse of data already collected for training the agent. We demonstrate the efficacy of adversarial evaluation on two standard domains: humanoid control and simulated driving. Experimental results show that our methods can find catastrophic failures and estimate failures rates of agents multiple orders of magnitude faster than standard evaluation schemes, in minutes to hours rather than days.


title: A Characterization of Controllability for Discrete-Time Linear Systems with Convex Constraints

abstract: We give a characterization of the controllability for discrete-time linear systems with convex output constraints. It extends all previously known characterizations in the literature, as well as our previous results on controllability of discrete-time linear systems with conic output constraints.


title: Instrument Simulator and Data Reduction Pipeline for the iLocater Spectrograph

abstract: iLocater is a near-infrared (NIR) radial velocity (RV) spectrograph that is being developed for the Large Binocular Telescope in Arizona. Unlike seeing limited designs, iLocater uses adaptive optics to inject starlight directly into a single mode fiber. This feature offers high spectral resolution while simultaneously maintaining a compact optical design. Although this approach shows promise to generate extremely precise RV measurements, it differs from conventional Doppler spectrographs, and therefore carries additional risk. To aid with the design of the instrument, we have developed a comprehensive simulator and data reduction pipeline. In this paper, we describe the simulation code and quantify its performance in the context of understanding terms in a RV error budget. We find that the program has an intrinsic precision of $\sigma < 5$ cm/s, thereby justifying its use in a number of instrument trade studies. The code is written in Matlab and available for download on GitHub.


title: An efficient approach to extract parameters from star cluster CMDs: fitCMD

abstract: This work presents an approach (fitCMD) designed to obtain a comprehensive set of astrophysical parameters from colour-magnitude diagrams (CMDs) of star clusters. Based on initial mass function (IMF) properties taken from isochrones, fitCMD searches for the values of total (or cluster) stellar mass, age, global metallicity, foreground reddening, distance modulus, and magnitude-dependent photometric completeness that produce the artificial CMD that best reproduces the observed one; photometric scatter is also taken into account in the artificial CMDs. Inclusion of photometric completeness proves to be an important feature of fitCMD, something that becomes apparent especially when luminosity functions are considered. These parameters are used to build a synthetic CMD that also includes photometric scatter. Residual minimization between the observed and synthetic CMDs leads to the best-fit parameters. When tested against artificial star clusters, fitCMD shows to be efficient both in terms of computational time and ability to recover the input values.


title: $D_5^{(1)}$- Geometric Crystal corresponding to the Dynkin spin node $i=5$ and its ultra-discretization

abstract: Let $g$ be an affine Lie algebra with index set $I = \{0, 1, 2, \cdots , n\}$ and $g^L$ be its Langlands dual. It is conjectured that for each Dynkin node $i \in I \setminus \{0\}$ the affine Lie algebra $g$ has a positive geometric crystal whose ultra-discretization is isomorphic to the limit of certain coherent family of perfect crystals for $g^L$. In this paper we construct a positive geometric crystal $V(D_5^{(1)})$ in the level zero fundamental spin $D_5^{(1)}$- module $W(\varpi_5)$. Then we define explicit $0$-action on the level $l$ known $D_5^{(1)}$- perfect crystal $B^{5, l}$ and show that $\{B^{5, l}\}_{l \geq 1}$ is a coherent family of perfect crystals with limit $B^{5, \infty}$. Finally we show that the ultra-discretization of $V(D_5^{(1)})$ is isomorphic to $B^{5, \infty}$ as crystals which proves the conjecture in this case.


title: Sub-micrometer Nanostructure-based RGB Filters for CMOS Image Sensors

abstract: Digital color imaging relies on spectral filters on top of a pixelated sensor, such as a CMOS image sensor. An important parameter of imaging devices is their resolution, which depends on the size of the pixels. For many applications, a high resolution is desirable, consequently requiring small spectral filters. Dielectric nanostructures, due to their resonant behavior and its tunability, offer the possibility to be assembled into flexible and miniature spectral filters, which could potentially replace conventional pigmented and dye-based color filters. In this paper, we demonstrate the generation of transmissive structural colors based on uniform-height amorphous silicon nanostructures. We optimize the structures for the primary RGB colors and report the construction of sub-micrometer RGB filter arrays for a pixel size down to 0.5 {\mu}m.


title: Model theory and metric convergence II: Averages of unitary polynomial actions

abstract: We use model theory of metric structures to prove the pointwise convergence, with a uniform metastability rate, of averages of a polynomial sequence $\{T_n\}$ (in Leibman's sense) of unitary transformations of a Hilbert space. As a special case, this applies to unitary sequences $\{U^{p(n)}\}$ where $p$ is a polynomial $\mathbb{Z}\to\mathbb{Z}$ and $U$ a fixed unitary operator; however, our convergence results hold for arbitrary Leibman sequences. As a case study, we show that the non-nilpotent "lamplighter group" $\mathbb{Z}\wr\mathbb{Z}$ is realized as the range of a suitable quadratic Leibman sequence. We also indicate how these convergence results generalize to arbitrary Folner averages of unitary polynomial actions of any abelian group $\mathbb{G}$ in place of $\mathbb{Z}$.


title: Tate cohomology of connected k-theory for elementary abelian groups revisited

abstract: Tate cohomology (as well as Borel homology and cohomology) of connective K-theory for $G=(\mathbb{Z}/2)^n$ was completely calculated by Bruner and Greenlees. In this note, we essentially redo the calculation by a different, more elementary method, and we extend it to $p>2$ prime. We also identify the resulting spectra, which are products of Eilenberg-Mac Lane spectra, and finitely many finite Postnikov towers. For $p=2$, we also reconcile our answer completely with the result of Bruner and Greenlees, which is in a different form, and hence the comparison involves some non-trivial combinatorics.


title: A probabilistic incremental proximal gradient method

abstract: In this paper, we propose a probabilistic optimization method, named probabilistic incremental proximal gradient (PIPG) method, by developing a probabilistic interpretation of the incremental proximal gradient algorithm. We explicitly model the update rules of the incremental proximal gradient method and develop a systematic approach to propagate the uncertainty of the solution estimate over iterations. The PIPG algorithm takes the form of Bayesian filtering updates for a state-space model constructed by using the cost function. Our framework makes it possible to utilize well-known exact or approximate Bayesian filters, such as Kalman or extended Kalman filters, to solve large-scale regularized optimization problems.


title: Chemical and kinematic analysis of CN-strong Metal-poor Field Stars in LAMOST DR3

abstract: The large amount of chemical and kinematic information available in large spectroscopic surveys have inspired the search for chemically peculiar stars in the field. Though these metal-poor field stars ([Fe/H$]<-1$) are commonly enriched in nitrogen, their detailed spatial, kinematic, and chemical distributions suggest that various groups may exist, and thus their origin is still a mystery. To study these stars statistically, we increase the sample size by identifying new CN-strong stars with LAMOST DR3 for the first time. We use CN-CH bands around 4000 \AA~to find CN-strong stars, and further separate them into CH-normal stars (44) and CH-strong (or CH) stars (35). The chemical abundances from our data-driven software and APOGEE DR 14 suggest that most CH-normal stars are N-rich, and it cannot be explained by only internal mixing process. The kinematics of our CH-normal stars indicate a substantial fraction of these stars are retrograding, pointing to an extragalactic origin. The chemistry and kinematics of CH-normal stars imply that they may be GC-dissolved stars, or accreted halo stars, or both.


title: Extension of reilly formula for a class of elliptic differential operator in divergence form

abstract: We prove the Reilly formula for a class of elliptic divergence differential operator $L_Au=div(A\nabla u)$, where $A$ is a (1,1)-Codazzi tensor field. Then we get some estimates for the first positive eigenvalue of the operator.