title: Effectiveness of Diffusing Information through a Social Network in Multiple Phases

abstract: We study the effectiveness of using multiple phases for maximizing the extent of information diffusion through a social network, and present insights while considering various aspects. In particular, we focus on the independent cascade model with the possibility of adaptively selecting seed nodes in multiple phases based on the observed diffusion in preceding phases, and conduct a detailed simulation study on real-world network datasets and various values of seeding budgets. We first present a negative result that more phases do not guarantee a better spread, however the adaptability advantage of more phases generally leads to a better spread in practice, as observed on real-world datasets. We study how diffusing in multiple phases affects the mean and standard deviation of the distribution representing the extent of diffusion. We then study how the number of phases impacts the effectiveness of multiphase diffusion, how the diffusion progresses phase-by-phase, and what is an optimal way to split the total seeding budget across phases. Our experiments suggest a significant gain when we move from single phase to two phases, and an appreciable gain when we further move to three phases, but the marginal gain thereafter is usually not very significant. Our main conclusion is that, given the number of phases, an optimal way to split the budget across phases is such that the number of nodes influenced in each phase is almost the same.


title: Topological nanophononic states by band inversion

abstract: Nanophononics is essential for the engineering of thermal transport in nanostructured electronic devices, it greatly facilitates the manipulation of mechanical resonators in the quantum regime, and could unveil a new route in quantum communications using phonons as carriers of information. Acoustic phonons also constitute a versatile platform for the study of fundamental wave dynamics, including Bloch oscillations, Wannier Stark ladders and other localization phenomena. Many of the phenomena studied in nanophononics were indeed inspired by their counterparts in optics and electronics. In these fields, the consideration of topological invariants to control wave dynamics has already had a great impact for the generation of robust confined states. Interestingly, the use of topological phases to engineer nanophononic devices remains an unexplored and promising field. Conversely, the use of acoustic phonons could constitute a rich platform to study topological states. Here, we introduce the concept of topological invariants to nanophononics and experimentally implement a nanophononic system supporting a robust topological interface state at 350 GHz. The state is constructed through band inversion, i.e. by concatenating two semiconductor superlattices with inverted spatial mode symmetries. The existence of this state is purely determined by the Zak phases of the constituent superlattices, i.e. that one-dimensional Berry phase. We experimentally evidenced the mode through Raman spectroscopy. The reported robust topological interface states could become part of nanophononic devices requiring resonant structures such as sensors or phonon lasers.


title: Integrating a Computational Perspective in Physics Courses

abstract: In this contribution we discuss how to develop a physics curriculum for undergraduate students that includes computing as a central element. Our contribution starts with a definition of computing and pertinent learning outcomes and assessment studies and programs. We end with a discussion on how to implement computing in various physics courses by presenting our experiences from Michigan State University in the USA and the University of Oslo in Norway.


title: Deep learning for conifer/deciduous classification of airborne LiDAR 3D point clouds representing individual trees

abstract: The purpose of this study was to investigate the use of deep learning for coniferous/deciduous classification of individual trees from airborne LiDAR data. To enable efficient processing by a deep convolutional neural network (CNN), we designed two discrete representations using leaf-off and leaf-on LiDAR data: a digital surface model with four channels (DSMx4) and a set of four 2D views (4x2D). A training dataset of labeled tree crowns was generated via segmentation of tree crowns, followed by co-registration with field data. Potential mislabels due to GPS error or tree leaning were corrected using a statistical ensemble filtering procedure. Because the training data was heavily unbalanced (~8% conifers), we trained an ensemble of CNNs on random balanced sub-samples of augmented data (180 rotational variations per instance). The 4x2D representation yielded similar classification accuracies to the DSMx4 representation (~82% coniferous and ~90% deciduous) while converging faster. The data augmentation improved the classification accuracies, but more real training instances (especially coniferous) likely results in much stronger improvements. Leaf-off LiDAR data were the primary source of useful information, which is likely due to the perennial nature of coniferous foliage. LiDAR intensity values also proved to be useful, but normalization yielded no significant improvements. Lastly, the classification accuracies of overstory trees (~90%) were more balanced than those of understory trees (~90% deciduous and ~65% coniferous), which is likely due to the incomplete capture of understory tree crowns via airborne LiDAR. Automatic derivation of optimal features via deep learning provide the opportunity for remarkable improvements in prediction tasks where captured data are not friendly to human visual system - likely yielding sub-optimal human-designed features.


title: On the Convergence Rates of GMsFEMs for Heterogeneous Elliptic Problems without Oversampling Techniques

abstract: This work is concerned with the rigorous analysis on the Generalized Multiscale Finite Element Methods (GMsFEMs) for elliptic problems with high-contrast heterogeneous coefficients. GMsFEMs are popular numerical methods for solving flow problems with heterogeneous high-contrast coefficients, and it has demonstrated extremely promising numerical results for a wide range of applications. However, the mathematical justification of the efficiency of the method is still largely missing. In this work, we analyze two types of multiscale basis functions, i.e., local spectral basis functions and basis functions of local harmonic extension type, within the GMsFEM framework. These constructions have found many applications in the past few years. We establish their optimal convergence in the energy norm under a very mild assumption that the source term belongs to some weighted $L^2$ space, and without the help of any oversampling technique. Furthermore, we analyze the model order reduction of the local harmonic extension basis and prove its convergence in the energy norm. These theoretical findings shed insights into the mechanism behind the efficiency of the GMsFEMs.


title: Phase-locked bi-frequency Raman lasing in a double-$\Lambda$ system

abstract: We show that it is possible to realize simultaneous Raman lasing at two different frequencies using a double-$\Lambda$ system pumped by a bi-frequency field. The Raman lasers are phase-locked to one another, and the beat-frequency matches the energy difference between the two meta-stable ground states. Akin to a conventional Raman laser, the phase-locked Raman laser pair is expected to be subluminal. As such, it is expected to be highly stable against perturbations in cavity length, and have a quantum noise limited linewidth that is far below that of a conventional laser. Because of these properties, the phase-locked Raman laser pair may find important applications in precision metrology, including atomic interferometry and magnetometry. To elucidate the behavior of this laser pair, we develop an analytical model that describes the stimulated Raman interaction in a double-$\Lambda$ system using an effective 2-level transition. The approximation is valid as long as the excited states adiabatically follow the ground states, as verified by numerical simulations. The effective model is used to identify the optimal operating conditions for the bi-frequency Raman lasing process. This model may also prove useful in other potential applications of the double-$\Lambda$ system, including generation of squeezed light and spatial solitons.


title: Dynamics of large-scale solar-wind streams obtained by the double superposed epoch analysis. 3. Deflection of speed vector

abstract: This work is a continuation of our previous articles (Yermolaev et al. in J. Geophys.Res. 120, 7094, 2015, Yermolaev et al. in Solar Phys. 292, 193, 2017), which describe the average temporal profiles of interplanetary plasma and field parameters in large-scale solar-wind (SW) streams: corotating interaction regions (CIRs), interplanetary coronal mass ejections (ICMEs including both magnetic clouds (MCs) and ejecta), and sheaths as well as interplanetary shocks (ISs). Changes of longitude angle {\phi} in CIRs from -2 to +2{\deg} agree with earlier observations. Besides we have for the first time analyzed the average temporal profiles of bulk velocity angles in Sheaths and ICMEs and found that the angle {\phi} in ICME changes from 2 to -2{\deg} while in Sheath it changes from -2 to 2{\deg} (similar to change in CIR), i.e., the streams in CIR/Sheath and ICME deflects in the opposite side. When averaging the latitude angle {\theta} on all intervals of the chosen SW type, the angle {\theta} is almost constant ~1{\deg}. We made selection of SW events with increasing and decreasing angle {\theta} for the first time and found that average temporal profiles for angle {\theta} in selected events have the same <<integral-like>> shape as for angle {\phi}. The difference in average profiles for angles {\phi} and {\theta} is explained by the fact that most events have increasing profiles for angle in ecliptic plane due to solar rotation while for angle in meridional plane numbers of increasing and decreasing profiles are equal.


title: Lov\'asz Meets Weisfeiler and Leman

abstract: In this paper, we relate a beautiful theory by Lov\'asz with a popular heuristic algorithm for the graph isomorphism problem, namely the color refinement algorithm and its k-dimensional generalization known as the Weisfeiler-Leman algorithm. We prove that two graphs G and H are indistinguishable by the color refinement algorithm if and only if, for all trees T, the number Hom(T,G) of homomorphisms from T to G equals the corresponding number Hom(T,H) for H. There is a natural system of linear equations whose nonnegative integer solutions correspond to the isomorphisms between two graphs. The nonnegative real solutions to this system are called fractional isomorphisms, and two graphs are fractionally isomorphic if and only if the color refinement algorithm cannot distinguish them (Tinhofer 1986, 1991). We show that, if we drop the nonnegativity constraints, that is, if we look for arbitrary real solutions, then a solution to the linear system exists if and only if, for all t, the two graphs have the same number of length-t walks. We lift the results for trees to an equivalence between numbers of homomorphisms from graphs of tree width k, the k-dimensional Weisfeiler-Leman algorithm, and the level-k Sherali-Adams relaxation of our linear program. We also obtain a partial result for graphs of bounded path width and solutions to our system where we drop the nonnegativity constraints. A consequence of our results is a quasi-linear time algorithm to decide whether, for two given graphs G and H, there is a tree T with Hom(T,G) = Hom(T,H).


title: Exponents of Bogomolov multipliers

abstract: We prove that if $G$ is a finite group, then the exponent of its Bogomolov multiplier divides the exponent of $G$ in the following four cases: (i) $G$ is metabelian, (ii) $\exp G=4$, (iii) $G$ is nilpotent of class $\le 5$, or (iv) $G$ is a $4$-Engel group.


title: Tales from the prehistory of Quantum Gravity. L\'eon Rosenfeld's earliest contribution

abstract: The main purpose of this paper is to analyse the earliest work of L\'eon Rosenfeld, one of the pioneers in the search of Quantum Gravity, the supposed theory unifying quantum theory and general relativity. We describe how and why Rosenfeld tried to face this problem in 1927, analysing the role of his mentors: Oskar Klein, Louis de Broglie and Th\'eophile De Donder. Rosenfeld asked himself how quantum mechanics should \textit{concretely} modify general relativity. In the context of a five-dimensional theory, Rosenfeld tried to construct a unifying framework for the gravitational and electromagnetic interaction and wave mechanics. Using a sort of "general relativistic quantum mechanics" Rosenfeld introduced a wave equation on a curved background. He investigated the metric created by what he called `quantum phenomena', represented by wave functions. Rosenfeld integrated Einstein equations in the weak field limit, with wave functions as source of the gravitational field. The author performed a sort of semi-classical approximation obtaining at the first order the Reissner-Nordstr\"om metric. We analyse how Rosenfeld's work is part of the history of Quantum Mechanics, because in his investigation Rosenfeld was guided by Bohr's correspondence principle. Finally we briefly discuss how his contribution is connected with the task of finding out which metric can be generated by a quantum field, a problem that quantum field theory on curved backgrounds will start to address 35 years later.


title: Multi-megawatt, self-seeded Mamyshev oscillator

abstract: We demonstrate a fiber oscillator that achieves 3 MW peak power, is easily started and is environmentally stable. The Mamyshev oscillator delivers 190-nJ pulses that can be compressed externally to 35 fs duration. Accurate numerical modeling of the gain medium provides insight into the behavior and performance of the device.


title: Asynchronous Stochastic Proximal Methods for Nonconvex Nonsmooth Optimization

abstract: We study stochastic algorithms for solving nonconvex optimization problems with a convex yet possibly nonsmooth regularizer, which find wide applications in many practical machine learning applications. However, compared to asynchronous parallel stochastic gradient descent (AsynSGD), an algorithm targeting smooth optimization, the understanding of the behavior of stochastic algorithms for nonsmooth regularized optimization problems is limited, especially when the objective function is nonconvex. To fill this theoretical gap, in this paper, we propose and analyze asynchronous parallel stochastic proximal gradient (Asyn-ProxSGD) methods for nonconvex problems. We establish an ergodic convergence rate of $O(1/\sqrt{K})$ for the proposed Asyn-ProxSGD, where $K$ is the number of updates made on the model, matching the convergence rate currently known for AsynSGD (for smooth problems). To our knowledge, this is the first work that provides convergence rates of asynchronous parallel ProxSGD algorithms for nonconvex problems. Furthermore, our results are also the first to show the convergence of any stochastic proximal methods without assuming an increasing batch size or the use of additional variance reduction techniques. We implement the proposed algorithms on Parameter Server and demonstrate its convergence behavior and near-linear speedup, as the number of workers increases, on two real-world datasets.


title: The effect of transmission-line dynamics on grid-forming dispatchable virtual oscillator control

abstract: In this work, we analyze the effect of transmission line dynamics on grid-forming control for inverter-based AC power systems. In particular, we investigate a dispatchable virtual oscillator control (dVOC) strategy that was recently proposed in the literature. When the dynamics of the transmission lines are neglected, i.e., if an algebraic model of the transmission network is used, dVOC ensures almost global asymptotic stability of a network of AC power inverters with respect to a pre-specified solution of the AC power-flow equations. While this approximation is typically justified for conventional power systems, the electromagnetic transients of the transmission lines can compromise the stability of an inverter-based power system. In this work, we establish explicit bounds on the controller setpoints, branch powers, and control gains that guarantee almost global asymptotic stability of dVOC in combination with a dynamic model of the transmission network.


title: A Block-wise, Asynchronous and Distributed ADMM Algorithm for General Form Consensus Optimization

abstract: Many machine learning models, including those with non-smooth regularizers, can be formulated as consensus optimization problems, which can be solved by the alternating direction method of multipliers (ADMM). Many recent efforts have been made to develop asynchronous distributed ADMM to handle large amounts of training data. However, all existing asynchronous distributed ADMM methods are based on full model updates and require locking all global model parameters to handle concurrency, which essentially serializes the updates from different workers. In this paper, we present a novel block-wise, asynchronous and distributed ADMM algorithm, which allows different blocks of model parameters to be updated in parallel. The lock-free block-wise algorithm may greatly speedup sparse optimization problems, a common scenario in reality, in which most model updates only modify a subset of all decision variables. We theoretically prove the convergence of our proposed algorithm to stationary points for non-convex general form consensus problems with possibly non-smooth regularizers. We implement the proposed ADMM algorithm on the Parameter Server framework and demonstrate its convergence and near-linear speedup performance as the number of workers increases.


title: Approximation of Kolmogorov-Smirnov Test Statistics

abstract: Motivated by the weak limit of the Kolmogorov-Smirnov test statistics, in this contribution, we concern the asymptotics of \begin{align*} \mathbb{P}\left\{\sup_{\boldsymbol{x}\in [0,1]^n}\left(W(\boldsymbol{x})\Big| W(\boldsymbol{1})=w\right)>u\right\}, \ w\in\mathbb{R}, \end{align*} for large $u$ where $W(\boldsymbol{x})$ is the multivariate Brownian sheet based on a distribution function $F$. The results related to general $F$ are investigated and some important examples are also showed.


title: The Duration of Energy Deposition on Unresolved Flaring Loops in the Solar Corona

abstract: Solar flares form and release energy across a large number of magnetic loops. The global parameters of flares, such as the total energy released, duration, physical size, etc., are routinely measured, and the hydrodynamics of a coronal loop subjected to intense heating have been extensively studied. It is not clear, however, how many loops comprise a flare, nor how the total energy is partitioned between them. In this work, we employ a hydrodynamic model to better understand the energy partition by synthesizing Si IV and Fe XXI line emission and comparing to observations of these lines with IRIS. We find that the observed temporal evolution of the Doppler shifts holds important information on the heating duration. To demonstrate this we first examine a single loop model, and find that the properties of chromospheric evaporation seen in Fe XXI can be reproduced by loops heated for long durations, while persistent red-shifts seen in Si IV cannot be reproduced by any single loop model. We then examine a multi-threaded model, assuming both a fixed heating duration on all loops, and a distribution of heating durations. For a fixed heating duration, we find that durations of 100 -- 200 s do a fair job of reproducing both the red- and blue-shifts, while a distribution of durations, with a median of about 50 -- 100 s, does a better job. Finally, we compare our simulations directly to observations of an M-class flare seen by IRIS, and find good agreement between the modeled and observed values given these constraints.


title: Importance of initial conditions in the polarization of complex networks

abstract: Currently used models of opinion formation use random initial conditions. In reality, most people in a social network, except for a small fraction of the population, are initially either unaware of, or indifferent to, the disputed issue. To explore the consequences of such specific initial conditions, we study the polarization of social networks when conflicting ideas arise on two different seed nodes and then spread according to a majority rule. Using the configuration model and the stochastic block model as examples, we show that this framework leads to substantially different outcomes than those which employ random initial conditions. Moreover, the empirically observed splits in the karate and the dolphins' networks naturally come out of this paradigm. Our work thus suggests that the existing opinion dynamics models should be reevaluated to incorporate the initial condition dependence.


title: Holomorphic torsion and geometric zeta functions for certain Hermitian locally symmetric manifolds

abstract: We give a dynamical description, in terms of a Weil-type zeta function, to the holomorphic torsion with coefficients for certain compact Hermitian locally symmetric manifolds, whose connected group G of isometries of the universal cover has only one conjugacy class of cuspidal maximal parabolic subgroup and satisfies a technical Ansatz relative to the given coefficients. A distinguishing feature of our zeta function is that its construction involves in an essential way the geometry of a standard compactification of the universal cover. The two senior authors are indebted to their junior colleague, Jan Frahm, for his laborious work shedding light on the scope of the validity of the Ansatz, and for writing up the attached Appendix. The results therein show that for real rank one groups G the Ansatz is satisfied with respect to any coefficients, for some rank two groups G it is satisfied with respect to certain coefficients, and also that there are groups G which do not obey the Ansatz.


title: Water from Two Rocks: Maximizing the Mutual Information

abstract: We build a natural connection between the learning problem, co-training, and forecast elicitation without verification (related to peer-prediction) and address them simultaneously using the same information theoretic approach. In co-training/multiview learning, the goal is to aggregate two views of data into a prediction for a latent label. We show how to optimally combine two views of data by reducing the problem to an optimization problem. Our work gives a unified and rigorous approach to the general setting. In forecast elicitation without verification we seek to design a mechanism that elicits high quality forecasts from agents in the setting where the mechanism does not have access to the ground truth. By assuming the agents' information is independent conditioning on the outcome, we propose mechanisms where truth-telling is a strict equilibrium for both the single-task and multi-task settings. Our multi-task mechanism additionally has the property that the truth-telling equilibrium pays better than any other strategy profile and strictly better than any other "non-permutation" strategy profile when the prior satisfies some mild conditions.


title: N-GCN: Multi-scale Graph Convolution for Semi-supervised Node Classification

abstract: Graph Convolutional Networks (GCNs) have shown significant improvements in semi-supervised learning on graph-structured data. Concurrently, unsupervised learning of graph embeddings has benefited from the information contained in random walks. In this paper, we propose a model: Network of GCNs (N-GCN), which marries these two lines of work. At its core, N-GCN trains multiple instances of GCNs over node pairs discovered at different distances in random walks, and learns a combination of the instance outputs which optimizes the classification objective. Our experiments show that our proposed N-GCN model improves state-of-the-art baselines on all of the challenging node classification tasks we consider: Cora, Citeseer, Pubmed, and PPI. In addition, our proposed method has other desirable properties, including generalization to recently proposed semi-supervised learning methods such as GraphSAGE, allowing us to propose N-SAGE, and resilience to adversarial input perturbations.


title: Remarks on the Ostrovsky's Theorem

abstract: In this paper we prove that the condition one-to-one of continuous open-resolvable mapping is necessary in the Ostrovsky's Theorem. Also we get that the Ostrovsky's Problem (Is every open-LCn function between Polish spaces piecewise open for n = 2, 3, ... ?) has a negative solution for any n > 1.


title: Scalar Resonant Relaxation of Stars Around a Massive Black Hole

abstract: In nuclear star clusters, the potential is governed by the central massive black hole, so that stars move on nearly Keplerian orbits and the total potential is almost stationary in time. Yet, the deviations of the potential from the Keplerian one, due to the enclosed stellar mass and general relativity, will cause the stellar orbits to precess. Moreover, as a result of the finite number of stars, small deviations of the potential from spherical symmetry induce residual torques that can change the stars' angular momentum faster than the standard two-body relaxation. The combination of these two effects drives a stochastic evolution of orbital angular momentum, a process named "resonant relaxation". Owing to recent developments in the description of the relaxation of self-gravitating systems, we can now fully describe scalar resonant relaxation (relaxation of the magnitude of the angular momentum) as a diffusion process. In this framework, the potential fluctuations due to the complex orbital motion of the stars are described by a random correlated noise, whose statistical properties are fully characterized by the stars' mean field motion. On long timescales, the cluster can be regarded as a diffusive system, whose diffusion coefficients depend explicitly on the mean field stellar distribution through the properties of the noise. We show here, for the first time, how the diffusion coefficients of scalar resonant relaxation, for a spherically symmetric system, can be fully calculated from first principles, without any free parameters. We also provide an open source code that evaluates these diffusion coefficients numerically.


title: Affine elliptic surfaces with type-A singularities and orbi-conifolds

abstract: Following the work of Castano-Bernard and Matessi on conifold transition in the Gross-Siebert program, we construct orbi-conifold transitions of the Shoen's Calabi-Yau threefold and their mirrors. The construction glues together the local models for orbi-conifold transitions in the previous work with Kanazawa.


title: Decompositions of linear spaces induced by $n$-linear maps

abstract: Let $\mathbb V$ be an arbitrary linear space and $f:\mathbb V \times \ldots \times \mathbb V \to \mathbb V$ an $n$-linear map. It is proved that, for each choice of a basis ${\mathcal B}$ of $\mathbb V$, the $n$-linear map $f$ induces a (nontrivial) decomposition $\mathbb V= \oplus V_j$ as a direct sum of linear subspaces of $\mathbb V$, with respect to ${\mathcal B}$. It is shown that this decomposition is $f$-orthogonal in the sense that $f(\mathbb V, \ldots, V_j, \ldots, V_k, \ldots, \mathbb V) =0$ when $j \neq k$, and in such a way that any $V_j$ is strongly $f$-invariant, meaning that $f(\mathbb V, \ldots, V_j, \ldots, \mathbb V) \subset V_j.$ A sufficient condition for two different decompositions of $\mathbb V$ induced by an $n$-linear map $f$, with respect to two different bases of $\mathbb V$, being isomorphic is deduced. The $f$-simplicity -- an analogue of the usual simplicity in the framework of $n$-liner maps -- of any linear subspace $V_j$ of a certain decomposition induced by $f$ is characterized. Finally, an application to the structure theory of arbitrary $n$-ary algebras is provided. This work is a close generalization the results obtained by A. J. Calder\'on (2018).


title: Accurate estimations of electromagnetic transitions of Sn IV for stellar and interstellar media

abstract: Here we report on accurate ab initio calculations to study astrophysically important electromagnetic transition parameters among different low-lying states of Sn IV. Our ab initio calculations are based on the sophisticated relativistic coupled-cluster theory, which almost exhausts many important electron correlations. To establish the accuracy of the calculations, we compare our results with the available experiments and estimates the transition amplitudes in length and velocity gauged forms. Most of these allowed and forbidden transition wavelengths lie in the infrared region, and they can be observed in the different cool stellar and interstellar media. For the improvement of uncertainty, we use experimental energies to the estimations of the above transition parameters. The presented data will be helpful to find the abundances of the ion in different astrophysical and laboratory plasma.


title: Improving Recall of In Situ Sequencing by Self-Learned Features and a Graphical Model

abstract: Image-based sequencing of mRNA makes it possible to see where in a tissue sample a given gene is active, and thus discern large numbers of different cell types in parallel. This is crucial for gaining a better understanding of tissue development and disease such as cancer. Signals are collected over multiple staining and imaging cycles, and signal density together with noise makes signal decoding challenging. Previous approaches have led to low signal recall in efforts to maintain high sensitivity. We propose an approach where signal candidates are generously included, and true-signal probability at the cycle level is self-learned using a convolutional neural network. Signal candidates and probability predictions are thereafter fed into a graphical model searching for signal candidates across sequencing cycles. The graphical model combines intensity, probability and spatial distance to find optimal paths representing decoded signal sequences. We evaluate our approach in relation to state-of-the-art, and show that we increase recall by $27\%$ at maintained sensitivity. Furthermore, visual examination shows that most of the now correctly resolved signals were previously lost due to high signal density. Thus, the proposed approach has the potential to significantly improve further analysis of spatial statistics in in situ sequencing experiments.


title: A Semi-Smooth Newton Algorithm for High-Dimensional Nonconvex Sparse Learning

abstract: The smoothly clipped absolute deviation (SCAD) and the minimax concave penalty (MCP) penalized regression models are two important and widely used nonconvex sparse learning tools that can handle variable selection and parameter estimation simultaneously, and thus have potential applications in various fields such as mining biological data in high-throughput biomedical studies. Theoretically, these two models enjoy the oracle property even in the high-dimensional settings, where the number of predictors $p$ may be much larger than the number of observations $n$. However, numerically, it is quite challenging to develop fast and stable algorithms due to their non-convexity and non-smoothness. In this paper we develop a fast algorithm for SCAD and MCP penalized learning problems. First, we show that the global minimizers of both models are roots of the nonsmooth equations. Then, a semi-smooth Newton (SSN) algorithm is employed to solve the equations. We prove that the SSN algorithm converges locally and superlinearly to the Karush-Kuhn-Tucker (KKT) points. Computational complexity analysis shows that the cost of the SSN algorithm per iteration is $O(np)$. Combined with the warm-start technique, the SSN algorithm can be very efficient and accurate. Simulation studies and a real data example suggest that our SSN algorithm, with comparable solution accuracy with the coordinate descent (CD) and the difference of convex (DC) proximal Newton algorithms, is more computationally efficient.


title: Vector-soliton storage and three-pulse-area theorem

abstract: In the present manuscript, we present a high-speed method to control, manipulate and retrieve an intense vector soliton stored in the ground state coherences of a four-level atomic system. Additionally, we show the importance of the pulse area in determining the evolution of the system and present a constant in the evolution defined as the three-pulse area, a surprising extension to previously defined pulse areas.


title: Applications of Picard and Magnus expansions to the Rabi model

abstract: We apply the Picard and Magnus expansions to both the semiclassical and the quantum Rabi model, with a switchable matter-field coupling. The case of the quantum Rabi model ia a paradigmatic example of finite-time quantum electrodynamics (QED), and in this case we build an intuitive diagrammatic representation of the Picard series. In particular, we show that regular oscillations in the mean number of photons, ascribed to the dynamical Casimir effect (DCE) for the the generation of photons and to the anti-DCE for their destruction, take place at twice the resonator frequency $\omega$. Such oscillations, which are a clear dynamical "smoking gun" of the ultrastrong coupling regime, can be predicted by first-order Picard expansion. We also show that the Magnus expansion can be used, through concatenation, as an efficient numerical integrator for both the semiclassical and the quantum Rabi model. In the first case, we find distinctive features in the Fourier spectrum of motion, with a single peak at the Rabi frequency $\Omega$ and doublets at frequencies $2n\omega\pm\Omega$, with $n$ positive integer. We explain these doublets, which are a feature beyond the rotating wave approximation (RWA), on the basis of the Picard series.


title: Dimensionally Tight Bounds for Second-Order Hamiltonian Monte Carlo

abstract: Hamiltonian Monte Carlo (HMC) is a widely deployed method to sample from high-dimensional distributions in Statistics and Machine learning. HMC is known to run very efficiently in practice and its popular second-order "leapfrog" implementation has long been conjectured to run in $d^{1/4}$ gradient evaluations. Here we show that this conjecture is true when sampling from strongly log-concave target distributions that satisfy a weak third-order regularity property associated with the input data. Our regularity condition is weaker than the Lipschitz Hessian property and allows us to show faster convergence bounds for a much larger class of distributions than would be possible with the usual Lipschitz Hessian constant alone. Important distributions that satisfy our regularity condition include posterior distributions used in Bayesian logistic regression for which the data satisfies an "incoherence" property. Our result compares favorably with the best available bounds for the class of strongly log-concave distributions, which grow like $d^{{1}/{2}}$ gradient evaluations with the dimension. Moreover, our simulations on synthetic data suggest that, when our regularity condition is satisfied, leapfrog HMC performs better than its competitors -- both in terms of accuracy and in terms of the number of gradient evaluations it requires.


title: Longitudinal Mapping Knot Invariant for SU(2)

abstract: The knot coloring polynomial defined by Eisermann for a finite pointed group is generalized to an infinite pointed group as the longitudinal mapping invariant of a knot. In turn this can be thought of as a generalization of the quandle 2-cocycle invariant for finite quandles. If the group is a topological group then this invariant can be thought of a topological generalization of the 2-cocycle invariant. The longitudinal mapping invariant is based on a meridian-longitude pair in the knot group. We also give an interpretation of the invariant in terms of quandle colorings of a 1-tangle for generalized Alexander quandles without use of a meridian-longitude pair in the knot group. The invariant values are concretely evaluated for the torus knots $T(2,n)$, their mirror images, and the figure eight knot for the group $SU(2)$.


title: Powers of tight Hamilton cycles in randomly perturbed hypergraphs

abstract: For $k\ge 2$ and $r\ge 1$ such that $k+r\ge 4$, we prove that, for any $\alpha>0$, there exists $\epsilon>0$ such that the union of an $n$-vertex $k$-graph with minimum codegree $\left(1-\binom{k+r-2}{k-1}^{-1}+\alpha\right)n$ and a binomial random $k$-graph $\mathbb{G}^{(k)}(n,p)$ with $p\ge n^{-\binom{k+r-2}{k-1}^{-1}-\epsilon}$ on the same vertex set contains the $r^{\text{th}}$ power of a tight Hamilton cycle with high probability. This result for $r=1$ was first proved by McDowell and Mycroft.


title: A quasi-physical dynamic reduced order model for thermospheric mass density via Hermitian Space Dynamic Mode Decomposition

abstract: Thermospheric mass density is a major driver of satellite drag, the largest source of uncertainty in accurately predicting the orbit of satellites in low Earth orbit (LEO) pertinent to space situational awareness. Most existing models for thermosphere are either physics-based or empirical. Physics-based models offer the potential for good predictive/forecast capabilities but require dedicated parallel resources for real-time evaluation and data assimilative capabilities that have yet to be developed. Empirical models are fast to evaluate, but offer very limited forecasting abilities. This paper presents a methodology of developing a reduced-order dynamic model from high-dimensional physics-based models by capturing the underlying dynamical behavior. This work develops a quasi-physical reduced order model (ROM) for thermospheric mass density using simulated output from NCAR's Thermosphere-Ionosphere-Electrodynamics General Circular Model (TIE-GCM). The ROM is derived using a dynamic system formulation from a large dataset of TIE-GCM simulations spanning 12 years and covering a complete solar cycle. Towards this end, a new reduced order modeling approach, based on Dynamic Mode Decomposition with control (DMDc), that uses the Hermitian space of the problem to derive the dynamics and input matrices in a tractable manner is developed. Results show that the ROM performs well in serving as a reduced order surrogate for TIE-GCM while almost always maintaining the forecast error to within 5\% of the simulated densities after 24 hours.


title: Galois descent for higher Brauer groups

abstract: For $X$ a smooth projective variety over a field $k$, we consider the problem of Galois descent for higher Brauer groups. More precisely, we extend a finiteness result of Colliot-Th\'el\`ene and Skorobogatov to higher Brauer groups.


title: Product Kernel Interpolation for Scalable Gaussian Processes

abstract: Recent work shows that inference for Gaussian processes can be performed efficiently using iterative methods that rely only on matrix-vector multiplications (MVMs). Structured Kernel Interpolation (SKI) exploits these techniques by deriving approximate kernels with very fast MVMs. Unfortunately, such strategies suffer badly from the curse of dimensionality. We develop a new technique for MVM based learning that exploits product kernel structure. We demonstrate that this technique is broadly applicable, resulting in linear rather than exponential runtime with dimension for SKI, as well as state-of-the-art asymptotic complexity for multi-task GPs.


title: Finding Long Lost Lexell's Comet: The Fate of the First Discovered Near-Earth Object

abstract: Jupiter-family Comet D/1770 L1 (Lexell) was the first discovered Near-Earth Object (NEO), and passed the Earth on 1770 Jul 1 at a recorded distance of 0.015 au. The comet was subsequently lost due to unfavorable observing circumstances during its next apparition followed by a close encounter with Jupiter in 1779. Since then, the fate of D/Lexell has attracted interest from the scientific community, and now we revisit this long-standing question. We investigate the dynamical evolution of D/Lexell based on a set of orbits recalculated using the observations made by Charles Messier, the comet's discoverer, and find that there is a $98\%$ chance that D/Lexell remains in the Solar System by the year of 2000. This finding remains valid even if a moderate non-gravitational effect is imposed. Messier's observations also suggest that the comet is one of the largest known near-Earth comets, with a nucleus of $\gtrsim 10$ km in diameter. This implies that the comet should have been detected by contemporary NEO surveys regardless of its activity level if it has remained in the inner Solar System. We identify asteroid 2010 JL$_{33}$ as a possible descendant of D/Lexell, with a $0.8\%$ probability of chance alignment, but a direct orbital linkage of the two bodies has not been successfully accomplished. We also use the recalculated orbit to investigate the meteors potentially originating from D/Lexell. While no associated meteors have been unambiguously detected, we show that meteor observations can be used to better constrain the orbit of D/Lexell despite the comet being long lost.


title: Bubble tree convergence for harmonic maps into compact locally CAT(1) spaces

abstract: We determine bubble tree convergence for a sequence of harmonic maps, with uniform energy bounds, from a compact Riemann surface into a compact locally CAT(1) space. In particular, we demonstrate energy quantization and the no-neck property for such a sequence. In the smooth setting, Jost and Parker respectively established these results by exploiting now classical arguments for harmonic maps. Our work demonstrates that these results can be reinterpreted geometrically. In the absence of a PDE, we take advantage of the local convexity properties of the target space. Included in this paper are an $\epsilon$-regularity theorem, an energy gap theorem, and a removable singularity theorem for harmonic maps for harmonic maps into metric spaces with upper curvature bounds. We also prove an isoperimetric inequality for conformal harmonic maps with small image.


title: Multi-frequency Raman amplifiers

abstract: In its usual implementation, the Raman amplifier features only one pump carrier frequency. However, pulses with well-separated frequencies can also be Raman amplified while compressed in time. Amplification with frequency-separated pumps is shown to hold even in the highly nonlinear, pump-depletion regime, as derived through a fluid model, and demonstrated via particle-in-cell (PIC) simulations. The resulting efficiency is similar to single-frequency amplifiers, but, due to the beat-wave waveform of both the pump lasers and the amplified seed pulses, these amplifiers feature higher seed intensities with a shorter spike duration. Advantageously, these amplifiers also suffer less noise backscattering, because the total fluence is split between the different spectral components.


title: Muon Hunter: a Zooniverse project

abstract: The large datasets and often low signal-to-noise inherent to the raw data of modern astroparticle experiments calls out for increasingly sophisticated event classification techniques. Machine learning algorithms, such as neural networks, have the potential to outperform traditional analysis methods, but come with the major challenge of identifying reliably classified training samples from real data. Citizen science represents an effective approach to sort through the large datasets efficiently and meet this challenge. Muon Hunter is a project hosted on the Zooniverse platform, wherein volunteers sort through pictures of data from the VERITAS cameras to identify muon ring images. Each image is classified multiple times to produce a "clean" dataset used to train and validate a convolutional neural network model both able to reject background events and identify suitable calibration data to monitor the telescope performance as a function of time.


title: Scalable Private Learning with PATE

abstract: The rapid adoption of machine learning has increased concerns about the privacy implications of machine learning models trained on sensitive data, such as medical records or other personal information. To address those concerns, one promising approach is Private Aggregation of Teacher Ensembles, or PATE, which transfers to a "student" model the knowledge of an ensemble of "teacher" models, with intuitive privacy provided by training teachers on disjoint data and strong privacy guaranteed by noisy aggregation of teachers' answers. However, PATE has so far been evaluated only on simple classification tasks like MNIST, leaving unclear its utility when applied to larger-scale learning tasks and real-world datasets. In this work, we show how PATE can scale to learning tasks with large numbers of output classes and uncurated, imbalanced training data with errors. For this, we introduce new noisy aggregation mechanisms for teacher ensembles that are more selective and add less noise, and prove their tighter differential-privacy guarantees. Our new mechanisms build on two insights: the chance of teacher consensus is increased by using more concentrated noise and, lacking consensus, no answer need be given to a student. The consensus answers used are more likely to be correct, offer better intuitive privacy, and incur lower-differential privacy cost. Our evaluation shows our mechanisms improve on the original PATE on all measures, and scale to larger tasks with both high utility and very strong privacy ($\varepsilon$ < 1.0).


title: Free-breathing cardiac MRI using bandlimited manifold modelling

abstract: We introduce a novel bandlimited manifold framework and an algorithm to recover freebreathing and ungated cardiac MR images from highly undersampled measurements. The image frames in the free breathing and ungated dataset are assumed to be points on a bandlimited manifold. We introduce a novel kernel low-rank algorithm to estimate the manifold structure (Laplacian) from a navigator-based acquisition scheme. The structure of the manifold is then used to recover the images from highly undersampled measurements. A computationally efficient algorithm, which relies on the bandlimited approximation of the Laplacian matrix, is used to recover the images. The proposed scheme is demonstrated on several patients with different breathing patterns and cardiac rates, without requiring the need for manually tuning the reconstruction parameters in each case. The proposed scheme enabled the recovery of free-breathing and ungated data, providing reconstructions that are qualitatively similar to breath-held scans performed on the same patients. This shows the potential of the technique as a clinical protocol for free-breathing cardiac scans.


title: Correlating Cellular Features with Gene Expression using CCA

abstract: To understand the biology of cancer, joint analysis of multiple data modalities, including imaging and genomics, is crucial. The involved nature of gene-microenvironment interactions necessitates the use of algorithms which treat both data types equally. We propose the use of canonical correlation analysis (CCA) and a sparse variant as a preliminary discovery tool for identifying connections across modalities, specifically between gene expression and features describing cell and nucleus shape, texture, and stain intensity in histopathological images. Applied to 615 breast cancer samples from The Cancer Genome Atlas, CCA revealed significant correlation of several image features with expression of PAM50 genes, known to be linked to outcome, while Sparse CCA revealed associations with enrichment of pathways implicated in cancer without leveraging prior biological understanding. These findings affirm the utility of CCA for joint phenotype-genotype analysis of cancer.

categories: eess.SP eess.IV q-bio.CB q-bio.QM stat.AP


title: Conformal bootstrap for percolation and polymers

abstract: The conformal bootstrap is applied to percolation and dilute self-avoiding polymers, two theories with Virasoro central charge $c=0$ in two dimensions. In both cases we propose a spectrum of operators motivated by Virasoro symmetry which is devoid of a stress energy tensor as an approximate means of enforcing $c=0$. Percolation is treated in $2\leq D \leq 6$ dimensions, and the self-avoiding walk in $2 \leq D \leq 4$.

categories: hep-th cond-mat.stat-mech math-ph math.MP


title: Large Scale Environment of a $z=6.61$ Luminous Quasar Probed by Ly$\alpha$ Emitters and Lyman Break Galaxies

abstract: Quasars (QSOs) hosting supermassive black holes are believed to reside in massive halos harboring galaxy overdensities. However, many observations revealed average or low galaxy densities around $z\gtrsim6$ QSOs. This could be partly because they measured galaxy densities in only tens of arcmin$^2$ around QSOs and might have overlooked potential larger scale galaxy overdensities. Some previous studies also observed only Lyman break galaxies (LBGs, massive older galaxies) and missed low mass young galaxies like Ly$\alpha$ emitters (LAEs) around QSOs. Here we present observations of LAE and LBG candidates in $\sim700$ arcmin$^2$ around a $z=6.61$ luminous QSO using Subaru Telescope Suprime-Cam with narrow/broadbands. We compare their sky distributions, number densities and angular correlation functions with those of LAEs/LBGs detected in the same manner and comparable data quality in our control blank field. In the QSO field, LAEs and LBGs are clustering in 4-20 comoving Mpc angular scales, but LAEs show mostly underdensity over the field while LBGs are forming $30\times60$ comoving Mpc$^2$ large scale structure containing 3-$7\sigma$ high density clumps. The highest density clump includes a bright (23.78 mag in the narrowband) extended ($\gtrsim 16$ kpc) Ly$\alpha$ blob candidate, indicative of a dense environment. The QSO could be part of the structure but is not located exactly at any of the high density peaks. Near the QSO, LAEs show underdensity while LBGs average to $4\sigma$ excess densities compared to the control field. If these environments reflect halo mass, the QSO may not be in the most massive halo, but still in a moderately massive one.


title: Constraining the $\bar{p}/p$ Ratio in TeV Cosmic Rays with Observations of the Moon Shadow by HAWC

abstract: An indirect measurement of the antiproton flux in cosmic rays is possible as the particles undergo deflection by the geomagnetic field. This effect can be measured by studying the deficit in the flux, or shadow, created by the Moon as it absorbs cosmic rays that are headed towards the Earth. The shadow is displaced from the actual position of the Moon due to geomagnetic deflection, which is a function of the energy and charge of the cosmic rays. The displacement provides a natural tool for momentum/charge discrimination that can be used to study the composition of cosmic rays. Using 33 months of data comprising more than 80 billion cosmic rays measured by the High Altitude Water Cherenkov (HAWC) observatory, we have analyzed the Moon shadow to search for TeV antiprotons in cosmic rays. We present our first upper limits on the $\bar{p}/p$ fraction, which in the absence of any direct measurements, provide the tightest available constraints of $\sim1\%$ on the antiproton fraction for energies between 1 and 10 TeV.


title: Breakdown of zero-energy quantum Hall state in graphene in the light of current fluctuations and shot noise

abstract: We have investigated the cross-over from Zener tunneling of single charge carriers to avalanche type of bunched electron transport in a suspended graphene Corbino disk in the zeroth Landau level. At low bias, we find a tunneling current that follows the gyrotropic Zener tunneling behavior. At larger bias, we find avalanche type of transport that sets in at a smaller current the larger the magnetic field is. The low-frequency noise indicates strong bunching of the electrons in the avalanches. On the basis of the measured low-frequency switching noise power, we deduce the characteristic switching rates of the avalanche sequence. The simultaneous microwave shot noise measurement also reveals intrinsic correlations within the avalanche pulses and indicate decrease of correlations with increasing bias.


title: Security: Doing Whatever is Needed... and Not a Thing More!

abstract: As malware, exploits, and cyber-attacks advance over time, so do the mitigation techniques available to the user. However, while attackers often abandon one form of exploitation in favor of a more lucrative one, mitigation techniques are rarely abandoned. Mitigations are rarely retired or disabled since proving they have outlived their usefulness is often impossible. As a result, performance overheads, maintenance costs, and false positive rates induced by the different mitigations accumulate, culminating in an outdated, inefficient, and costly security solution. We advocate for a new kind of tunable framework on which to base security mechanisms. This new framework enables a more reactive approach to security allowing us to optimize the deployment of security mechanisms based on the current state of attacks. Based on actual evidence of exploitation collected from the field, our framework can choose which mechanisms to enable/disable so that we can minimize the overall costs and false positive rates while maintaining a satisfactory level of security in the system. We use real-world Snort signatures to simulate the benefits of reactively disabling signatures when no evidence of exploitation is observed and compare them to the costs of the current state of deployment. Additionally, we evaluate the responsiveness of our framework and show that in case disabling a security mechanism triggers a reappearance of an attack we can respond in time to prevent mass exploitation. Through large-scale simulations that use integer linear and Bayesian solvers, we discover that our responsive strategy is both computationally affordable and results in significant reductions in false positives (~20% over traces that are about 9 years long), at the cost of introducing a moderate number of false negatives. Finding the optimal sampling strategy takes less than 2.5 minutes in the vast majority of cases.


title: SAT-based Reverse Engineering of Gate-Level Schematics using Fault Injection and Probing

abstract: Gate camouflaging is a known security enhancement technique that tries to thwart reverse engineering by hiding the functions of gates or the connections between them. A number of works on SAT-based attacks have shown that it is often possible to reverse engineer a circuit function by combining a camouflaged circuit model and the ability to have oracle access to the obfuscated combinational circuit. Especially in small circuits it is easy to reverse engineer the circuit function in this way, but SAT-based reverse engineering techniques provide no guarantees of recovering a circuit that is gate-by-gate equivalent to the original design. In this work we show that an attacker who does not know gate functions or connections of an aggressively camouflaged circuit cannot learn the correct gate-level schematic even if able to control inputs and probe all combinational nodes of the circuit. We then present a stronger attack that extends SAT-based reverse engineering with fault analysis to allow an attacker to recover the correct gate-level schematic. We analyze our reverse engineering approach on an S-Box circuit.


title: Permissive Barrier Certificates for Safe Stabilization Using Sum-of-squares

abstract: Motivated by the need to simultaneously guarantee safety and stability of safety-critical dynamical systems, we construct permissive barrier certificates in this paper that explicitly maximize the region where the system can be stabilized without violating safety constraints. An optimization strategy is developed to search for the maximum volume barrier certified region of safe stabilization. The barrier certified region, which is allowed to take any arbitrary shape, is proved to be strictly larger than safe regions generated with Lyapunov sublevel set based methods. The proposed approach effectively unites a Lyapunov function with multiple barrier functions that might not be compatible with each other. Iterative search algorithms are developed using sum-of-squares to compute the most permissive, that is, the maximum volume, barrier certificates. Simulation results of the iterative search algorithm demonstrate the effectiveness of the proposed method.


title: Anti-Ramsey number of edge-disjoint rainbow spanning trees

abstract: An edge-colored graph $G$ is called rainbow if every edge of $G$ receives a different color. The anti-Ramsey number of $t$ edge-disjoint rainbow spanning trees, denoted by $r(n,t)$, is defined as the maximum number of colors in an edge-coloring of $K_n$ containing no $t$ edge-disjoint rainbow spanning trees. Jahanbekam and West [J. Graph Theory, 2014] conjectured that for any fixed $t$, $r(n,t)=\binom{n-2}{2}+t$ whenever $n\geq 2t+2 \geq 6$. In this paper, we prove this conjecture. We also determine $r(n,t)$ when $n = 2t+1$. Together with previous results, this gives the anti-Ramsey number of $t$ edge-disjoint rainbow spanning trees for all values of $n$ and $t$.


title: Privacy Leakages in Approximate Adders

abstract: Approximate computing has recently emerged as a promising method to meet the low power requirements of digital designs. The erroneous outputs produced in approximate computing can be partially a function of each chip's process variation. We show that, in such schemes, the erroneous outputs produced on each chip instance can reveal the identity of the chip that performed the computation, possibly jeopardizing user privacy. In this work, we perform simulation experiments on 32-bit Ripple Carry Adders, Carry Lookahead Adders, and Han-Carlson Adders running at over-scaled operating points. Our results show that identification is possible, we contrast the identifiability of each type of adder, and we quantify how success of identification varies with the extent of over-scaling and noise. Our results are the first to show that approximate digital computations may compromise privacy. Designers of future approximate computing systems should be aware of the possible privacy leakages and decide whether mitigation is warranted in their application.


title: Geometric Surface-Based Tracking Control of a Quadrotor UAV

abstract: New quadrotor UAV control algorithms are developed, based on nonlinear surfaces composed of tracking errors that evolve directly on the nonlinear configuration manifold, thus inherently including in the control design the nonlinear characteristics of the SE(3) configuration space. In particular, geometric surface-based controllers are developed and are shown, through rigorous stability proofs, to have desirable almost global closed loop properties. For the first time in regards to the geometric literature, a region of attraction independent of the position error is identified and its effects are analyzed. The effectiveness of the proposed "surface based" controllers are illustrated by simulations of aggressive maneuvers in the presence of disturbances and motor saturation.


title: Rigidity of $\ell^p$ Roe-type algebras

abstract: We investigate the rigidity of the $\ell^p$ analog of Roe-type algebras. In particular, we show that if $p\in[1,\infty)\setminus\{2\}$, then an isometric isomorphism between the $\ell^p$ uniform Roe algebras of two metric spaces with bounded geometry yields a bijective coarse equivalence between the underlying metric spaces, while a stable isometric isomorphism yields a coarse equivalence. We also obtain similar results for other $\ell^p$ Roe-type algebras. In this paper, we do not assume that the metric spaces have Yu's property A or finite decomposition complexity.


title: Positively curved Killing foliations via deformations

abstract: We show that a compact manifold admitting a Killing foliation with positive transverse curvature fibers over finite quotients of spheres or weighted complex projective spaces, provided that the singular foliation defined by the closures of the leaves has maximal dimension. This result is obtained by deforming the foliation into a closed one while maintaining transverse geometric properties, which allows us to apply results from the Riemannian geometry of orbifolds to the space of leaves. We also show that the basic Euler characteristic is preserved by such deformations. Using this fact we prove that a Riemannian foliation of a compact manifold with finite fundamental group and nonvanishing Euler characteristic is closed. As another application we obtain that, for a positively curved Killing foliation of a compact manifold, if the structural algebra has sufficiently large dimension then the basic Euler characteristic is positive.


title: The Strong Trotter Property for Locally $\mu$-convex Lie Groups

abstract: We show that an infinite dimensional Lie group in Milnor's sense has the strong Trotter property if it is locally $\mu$-convex. This is a continuity condition imposed on the Lie group multiplication that generalizes the triangle inequality for locally convex vector spaces, and is equivalent to $C^0$-continuity of the evolution map on its domain. In particular, the result proven in this paper significantly extends the respective result obtained by Gl\"ockner in the context of measurable regularity.


title: Time Series Learning using Monotonic Logical Properties

abstract: Cyber-physical systems of today are generating large volumes of time-series data. As manual inspection of such data is not tractable, the need for learning methods to help discover logical structure in the data has increased. We propose a logic-based framework that allows domain-specific knowledge to be embedded into formulas in a parametric logical specification over time-series data. The key idea is to then map a time series to a surface in the parameter space of the formula. Given this mapping, we identify the Hausdorff distance between boundaries as a natural distance metric between two time-series data under the lens of the parametric specification. This enables embedding non-trivial domain-specific knowledge into the distance metric and then using off-the-shelf machine learning tools to label the data. After labeling the data, we demonstrate how to extract a logical specification for each label. Finally, we showcase our technique on real world traffic data to learn classifiers/monitors for slow-downs and traffic jams.


title: Generating retinal flow maps from structural optical coherence tomography with artificial intelligence

abstract: Despite significant advances in artificial intelligence (AI) for computer vision, its application in medical imaging has been limited by the burden and limits of expert-generated labels. We used images from optical coherence tomography angiography (OCTA), a relatively new imaging modality that measures perfusion of the retinal vasculature, to train an AI algorithm to generate vasculature maps from standard structural optical coherence tomography (OCT) images of the same retinae, both exceeding the ability and bypassing the need for expert labeling. Deep learning was able to infer perfusion of microvasculature from structural OCT images with similar fidelity to OCTA and significantly better than expert clinicians (P < 0.00001). OCTA suffers from need of specialized hardware, laborious acquisition protocols, and motion artifacts; whereas our model works directly from standard OCT which are ubiquitous and quick to obtain, and allows unlocking of large volumes of previously collected standard OCT data both in existing clinical trials and clinical practice. This finding demonstrates a novel application of AI to medical imaging, whereby subtle regularities between different modalities are used to image the same body part and AI is used to generate detailed and accurate inferences of tissue function from structure imaging.


title: Global existence and stability of nearly aligned flocks

abstract: We study regularity of a hydrodynamic singular model of collective behavior introduced in \cite{ST1}. In this note we address the question of global well-posedness in multi-dimensional settings. It is shown that any initial data $(u,\rho)$ with small velocity variations $|u(x) - u(y)| < \e$ relative to its higher order norms, gives rise to a unique global regular solution which aligns and flocks exponentially fast. Moreover, we prove that the limiting flocks are stable.


title: Evaluating Design Tradeoffs in Numeric Static Analysis for Java

abstract: Numeric static analysis for Java has a broad range of potentially useful applications, including array bounds checking and resource usage estimation. However, designing a scalable numeric static analysis for real-world Java programs presents a multitude of design choices, each of which may interact with others. For example, an analysis could handle method calls via either a top-down or bottom-up interprocedural analysis. Moreover, this choice could interact with how we choose to represent aliasing in the heap and/or whether we use a relational numeric domain, e.g., convex polyhedra. In this paper, we present a family of abstract interpretation-based numeric static analyses for Java and systematically evaluate the impact of 162 analysis configurations on the DaCapo benchmark suite. Our experiment considered the precision and performance of the analyses for discharging array bounds checks. We found that top-down analysis is generally a better choice than bottom-up analysis, and that using access paths to describe heap objects is better than using summary objects corresponding to points-to analysis locations. Moreover, these two choices are the most significant, while choices about the numeric domain, representation of abstract objects, and context-sensitivity make much less difference to the precision/performance tradeoff.


title: Curvature estimates and sheeting theorems for weakly stable CMC hypersurfaces

abstract: Weakly stable constant mean curvature (CMC) hypersurfaces are stable critical points of the area functional with respect to volume preserving deformations. We establish a pointwise curvature estimate (in the non-singular dimensions) and a sheeting theorem (in all dimensions) for weakly stable CMC hypersurfaces, giving an effective version of the compactness theorem for weakly stable CMC hypersurfaces established in the recent work of the first and third-named authors. Our results generalize the curvature estimate and the sheeting theorem proven respectively by Schoen--Simon--Yau and Schoen--Simon for strongly stable hypersurfaces.


title: Market based embedded Real Time Operation for Distributed Resources and Flexibility

abstract: We build upon previous work out of UC Berkeley's energy, controls, and applications laboratory (eCal) that developed a model for price prediction of the energy day-ahead market (DAM) and a stochastic load scheduling for distributed energy resources (DER) with DAM based objective\cite{Travacca}. In similar fashion to the work of Travacca et al., in this project we take the standpoint of a DER aggregator pooling a large number of electricity consumers - each of which have an electric vehicle and solar PV panels - to bid their pooled energy resources into the electricity markets. The primary contribution of this project is the optimization of an aggregated load schedule for participation in the California Independent System Operator (CAISO) real time (15-minute) electricity market. The goal of the aggregator is to optimally manage its pool of resources, particularly the flexible resources, in order to minimize its cost in the real time market. We achieves this through the use of a model predictive control scheme. A critical difference between the prior work in \cite{Travacca} is that the structure of the optimization problem is drastically different. Based upon our review of the current and public literature, no similar approaches exist. The main objective of this project were building methods. Nevertheless, to illustrate a simulation with 100 prosumers was realized. The results should therefore be taken with a grain of salt. We find that the Real Time operation does not substantially decrease or increase the total cost the aggregator faces in the RT market, but this is probably due to parameters that need further tuning and data, that need better processing.


title: Twisted Calabi-Yau ring spectra, string topology, and gauge symmetry

abstract: In this paper, we import the theory of "Calabi-Yau" algebras and categories from symplectic topology and topological field theories to the setting of spectra in stable homotopy theory. Twistings in this theory will be particularly important. There will be two types of Calabi-Yau structures in the setting of ring spectra: one that applies to compact algebras and one that applies to smooth algebras. The main application of twisted compact Calabi-Yau ring spectra that we will study is to describe, prove, and explain a certain duality phenomenon in string topology. This is a duality between the manifold string topology of Chas-Sullivan and the Lie group string topology of Chataur-Menichi. This will extend and generalize work of Gruher. Then, generalizing work of the first author and Jones, we show how the gauge group of the principal bundle acts on this compact Calabi-Yau structure, and compute some explicit examples. We then extend the notion of the Calabi-Yau structure to smooth ring spectra, and prove that Thom ring spectra of (virtual) bundles over the loop space, $\Omega M$, have this structure. In the case when $M$ is a sphere we will use these twisted smooth Calabi-Yau ring spectra to study Lagrangian immersions of the sphere into its cotangent bundle. We recast the work of Abouzaid-Kragh to show that the topological Hochschild homology of the Thom ring spectrum induced by the $h$-principle classifying map of the Lagrangian immersion, detects whether that immersion can be Lagrangian isotopic to an embedding. We then compute some examples. Finally, we interpret these Calabi-Yau structures directly in terms of topological Hochschild homology and cohomology.


title: Fast, accurate solutions for curvilinear earthquake faults and anelastic strain

abstract: Imaging the anelastic deformation within the crust and lithosphere using surface geophysical data remains a significant challenge in part due to the wide range of physical processes operating at different depths and to various levels of localization that they embody. Models of Earth's elastic properties from seismological imaging combined with geodetic modeling may form the basis of comprehensive rheological models of Earth's interior. However, representing the structural complexity of faults and shear zones in numerical models of deformation still constitutes a major difficulty. Here, we present numerical techniques for high-precision models of deformation and stress around both curvilinear faults and volumes undergoing anelastic (irreversible) strain in a heterogenous elastic half-space. To that end, we enhance the software Gamra to model triangular and rectangular fault patches and tetrahedral and cuboidal strain volumes. This affords a means of rapid and accurate calculations of elasto-static Green's functions for localized (e.g., faulting) and distributed (e.g., viscoelastic) deformation in Earth's crust and lithosphere. We demonstrate the correctness of the method with analytic tests, and we illustrate its practical performance by solving for coseismic and postseismic deformation following the 2015 Mw 7.8 Gorkha, Nepal earthquake to extremely high precision.


title: Auxiliary field diffusion Monte Carlo calculations of light and medium-mass nuclei with local chiral interactions

abstract: Quantum Monte Carlo methods have recently been employed to study properties of nuclei and infinite matter using local chiral effective field theory interactions. In this work, we present a detailed description of the auxiliary field diffusion Monte Carlo algorithm for nuclei in combination with local chiral two- and three-nucleon interactions up to next-to-next-to-leading order. We show results for the binding energy, charge radius, charge form factor, and Coulomb sum rule in nuclei with $3\le A\le16$. Particular attention is devoted to the effect of different operator structures in the three-body force for different cutoffs. The outcomes suggest that local chiral interactions fit to few-body observables give a very good description of the ground-state properties of nuclei up to $^{16}$O, with the exception of one fit for the softer cutoff which predicts overbinding in larger nuclei.


title: Circular support in random sorting networks

abstract: A sorting network is a shortest path from $12 \cdots n$ to $n \cdots 2 1$ in the Cayley graph of the symmetric group generated by adjacent transpositions. For a uniform random sorting network, we prove that in the global limit, particle trajectories are supported on $\pi$-Lipschitz paths. We show that the weak limit of the permutation matrix of a random sorting network at any fixed time is supported within a particular ellipse. This is conjectured to be an optimal bound on the support. We also show that in the global limit, trajectories of particles that start within distance $\epsilon$ of the edge are within $\sqrt{2\epsilon}$ of a sine curve in uniform norm.


title: The Archimedean limit of random sorting networks

abstract: A sorting network (also known as a reduced decomposition of the reverse permutation), is a shortest path from $12 \cdots n$ to $n \cdots 21$ in the Cayley graph of the symmetric group $S_n$ generated by adjacent transpositions. We prove that in a uniform random $n$-element sorting network $\sigma^n$, all particle trajectories are close to sine curves with high probability. We also find the weak limit of the time-$t$ permutation matrix measures of $\sigma^n$. As a corollary of these results, we show that if $S_n$ is embedded into $\mathbb{R}^n$ via the map $\tau \mapsto (\tau(1), \tau(2), \dots \tau(n))$, then with high probability, the path $\sigma^n$ is close to a great circle on a particular $(n-2)$-dimensional sphere in $\mathbb{R}^n$. These results prove conjectures of Angel, Holroyd, Romik, and Virag.


title: Identifying the occurrence or non occurrence of cognitive bias in situations resembling the Monty Hall problem

abstract: People reason heuristically in situations resembling inferential puzzles such as Bertrand's box paradox and the Monty Hall problem. The practical significance of that fact for economic decision making is uncertain because a departure from sound reasoning may, but does not necessarily, result in a "cognitively biased" outcome different from what sound reasoning would have produced. Criteria are derived here, applicable to both experimental and non-experimental situations, for heuristic reasoning in an inferential-puzzle situations to result, or not to result, in cognitively bias. In some situations, neither of these criteria is satisfied, and whether or not agents' posterior probability assessments or choices are cognitively biased cannot be determined.


title: A Dataset To Evaluate The Representations Learned By Video Prediction Models

abstract: We present a parameterized synthetic dataset called Moving Symbols to support the objective study of video prediction networks. Using several instantiations of the dataset in which variation is explicitly controlled, we highlight issues in an existing state-of-the-art approach and propose the use of a performance metric with greater semantic meaning to improve experimental interpretability. Our dataset provides canonical test cases that will help the community better understand, and eventually improve, the representations learned by such networks in the future. Code is available at https://github.com/rszeto/moving-symbols .


title: Detecting Comma-shaped Clouds for Severe Weather Forecasting using Shape and Motion

abstract: Meteorologists use shapes and movements of clouds in satellite images as indicators of several major types of severe storms. Satellite imaginary data are in increasingly higher resolution, both spatially and temporally, making it impossible for humans to fully leverage the data in their forecast. Automatic satellite imagery analysis methods that can find storm-related cloud patterns as soon as they are detectable are in demand. We propose a machine learning and pattern recognition based approach to detect "comma-shaped" clouds in satellite images, which are specific cloud distribution patterns strongly associated with the cyclone formulation. In order to detect regions with the targeted movement patterns, our method is trained on manually annotated cloud examples represented by both shape and motion-sensitive features. Sliding windows in different scales are used to ensure that dense clouds will be captured, and we implement effective selection rules to shrink the region of interest among these sliding windows. Finally, we evaluate the method on a hold-out annotated comma-shaped cloud dataset and cross-match the results with recorded storm events in the severe weather database. The validated utility and accuracy of our method suggest a high potential for assisting meteorologists in weather forecasting.


title: DID: Distributed Incremental Block Coordinate Descent for Nonnegative Matrix Factorization

abstract: Nonnegative matrix factorization (NMF) has attracted much attention in the last decade as a dimension reduction method in many applications. Due to the explosion in the size of data, naturally the samples are collected and stored distributively in local computational nodes. Thus, there is a growing need to develop algorithms in a distributed memory architecture. We propose a novel distributed algorithm, called \textit{distributed incremental block coordinate descent} (DID), to solve the problem. By adapting the block coordinate descent framework, closed-form update rules are obtained in DID. Moreover, DID performs updates incrementally based on the most recently updated residual matrix. As a result, only one communication step per iteration is required. The correctness, efficiency, and scalability of the proposed algorithm are verified in a series of numerical experiments.


title: Constraining the Collective Radio Emission of Large Scale Accretion Shocks

abstract: Accretion of gas onto already virialized structures like galaxy clusters should give rise to accretion shocks which can potentially accelerate cosmic rays. Here, we use the radio emission detected from Coma cluster and models of evolution of cosmic accretion shocks, to constrain the possible contribution of unresolved galaxy clusters to the cosmic radio background. We assume that Coma is a typical galaxy cluster and that its entire radio emission is produced by cosmic rays accelerated in accretion shocks, making our prediction an upper limit. Our models predict that at lower frequencies accretion shocks can have a potentially large contribution to the cosmic radio background, while on larger frequencies, e.g. 5 GHz, their contribution must be lower than < 2-35%, depending on the models of evolution of accretion shocks that we use.


title: M\"ossbauer spectroscopy study of magnetic fluctuations in superconducting RbGd$_2$Fe$_4$As$_4$O$_2$

abstract: $^{57}$Fe M\"ossbauer spectra were measured at different temperatures between 5.9 K and 300 K on the recently discovered self-doped superconducting RbGd$_2$Fe$_4$As$_4$O$_2$ with T$_c$ as high as 35 K. Singlet pattern was observed down to the lowest temperature measured in this work, indicating the absence of static magnetic order on the Fe site. The intermediate isomer shift in comparison with that of the samples RbFe$_2$As$_2$ and GdFeAsO confirms the self doping induced local electronic structure change. Surprisingly, we observe two magnetic fluctuation induced spectral broadenings below $\sim$15 K and $\sim$100 K which are believed to be originated from the transferred magnetic fluctuations of the Gd$^{3+}$ moments and that of the magnetic fluctuations of the Fe atoms, respectively.


title: Gradient Primal-Dual Algorithm Converges to Second-Order Stationary Solutions for Nonconvex Distributed Optimization

abstract: In this work, we study two first-order primal-dual based algorithms, the Gradient Primal-Dual Algorithm (GPDA) and the Gradient Alternating Direction Method of Multipliers (GADMM), for solving a class of linearly constrained non-convex optimization problems. We show that with random initialization of the primal and dual variables, both algorithms are able to compute second-order stationary solutions (ss2) with probability one. This is the first result showing that primal-dual algorithm is capable of finding ss2 when only using first-order information, it also extends the existing results for first-order, but primal-only algorithms. An important implication of our result is that it also gives rise to the first global convergence result to the ss2, for two classes of unconstrained distributed non-convex learning problems over multi-agent networks.


title: Local quantum uncertainty guarantees the measurement precision for two coupled two-level systems in non-Markovian environment

abstract: Quantum Fisher information (QFI) is an important feature for the precision of quantum parameter estimation based on the quantum Cram\'er-Rao inequality. When the quantum state satisfies the von Neumann-Landau equation, the local quantum uncertainty (LQU), as a kind of quantum correlation, present in a bipartite mixed state guarantees a lower bound on QFI in the optimal phase estimation protocol [Phys. Rev. Lett. 110 (2013) 240402]. However, in the open quantum systems, there is not an explicit relation between LQU and QFI generally. In this paper, we study the relation between LQU and QFI in open systems which is composed of two interacting two-level systems coupled to independent non-Markovian environments with the entangled initial state embedded by a phase parameter $\theta$. The analytical calculations show that the QFI does't depend on the phase parameter $\theta$, and its decay can be restrained through enhancing the coupling strength or non-Markovianity. Meanwhile, the LQU is related to the phase parameter $\theta$ and shows plentiful phenomena. In particular, we find that the LQU can well bound the QFI when the coupling between the two systems is switched off or the initial state is Bell state.


title: A Model for Innovation Diffusion with Intergroup Suppression

abstract: We present a new model for the diffusion of innovation. Here, the population is segmented into distinct groups. Adoption by a particular group of some cultural product may be inhibited both by large numbers of its own members already having adopted but also, in particular, by members of another group having adopted. Intergroup migration is also permitted. We determine the equilibrium points and carry out stability analysis for the model for a two-group population. We also simulate a discrete time version of the model. Lastly, we present data on tablet use in eight countries from 2012-2016 and show that the relationship between use in the "under 25" age group and "55+" age group conforms to the model.


title: Real Time Quantum Dynamics of Spontaneous Translational Symmetry Breakage in the Early Stage of Photo-induced Structural Phase Transitions

abstract: Real time quantum dynamics of the spontaneous translational symmetry breakage in the early stage of photoinduced structural phase transitions is reviewed and supplementally explained, under the guide of the Toyozawa theory, which is exactly in compliance with the conservation laws of the total momentum and energy. At the Franck Condon state, an electronic excitation just created by a visible light, is in a plane wave state, extended all over the crystal. While, after the lattice relaxation having been completed, it is localized around a certain lattice site of the crystal, as a new excitation. Is there a sudden shrinkage of the excitation wave function, in between. The wave function never shrinks, but only the spatial, or inter lattice site quantum coherence, interference of the excitation disappears, as the lattice relaxation proceeds. This is nothing but the spontaneous breakage of translational symmetry.


title: Validation and Initial Characterization of the Long Period Planet Kepler-1654 b

abstract: Fewer than 20 transiting Kepler planets have periods longer than one year. Our early search of the Kepler light curves revealed one such system, Kepler-1654 b (originally KIC~8410697b), which shows exactly two transit events and whose second transit occurred only 5 days before the failure of the second of two reaction wheels brought the primary Kepler mission to an end. A number of authors have also examined light curves from the Kepler mission searching for long period planets and identified this candidate. Starting in Sept. 2014 we began an observational program of imaging, reconnaissance spectroscopy and precision radial velocity measurements which confirm with a high degree of confidence that Kepler-1654 b is a {\it bona fide} transiting planet orbiting a mature G2V star (T$_{eff}= 5580$K, [Fe/H]=-0.08) with a semi-major axis of 2.03 AU, a period of 1047.84 days and a radius of 0.82$\pm$0.02 R$_{Jup}$. Radial Velocity (RV) measurements using Keck's HIRES spectrometer obtained over 2.5 years set a limit to the planet's mass of $<0.5\ (3\sigma$) M$_{Jup}$. The bulk density of the planet is similar to that of Saturn or possibly lower. We assess the suitability of temperate gas giants like Kepler-1654b for transit spectroscopy with the James Webb Space Telescope since their relatively cold equilibrium temperatures (T$_{pl}\sim 200$K) make them interesting from the standpoint of exo-planet atmospheric physics. Unfortunately, these low temperatures also make the atmospheric scale heights small and thus transmission spectroscopy challenging. Finally, the long time between transits can make scheduling JWST observations difficult---as is the case with Kepler-1654b.


title: Teacher Improves Learning by Selecting a Training Subset

abstract: We call a learner super-teachable if a teacher can trim down an iid training set while making the learner learn even better. We provide sharp super-teaching guarantees on two learners: the maximum likelihood estimator for the mean of a Gaussian, and the large margin classifier in 1D. For general learners, we provide a mixed-integer nonlinear programming-based algorithm to find a super teaching set. Empirical experiments show that our algorithm is able to find good super-teaching sets for both regression and classification problems.


title: $2$-groups behaving as automorphism groups of regular $3$-polytopes

abstract: In this paper, we classify regular polytopes with automorphism groups of order $2^n$ and Schl\"afli types $\{4, 2^{n-3}\}, \{4, 2^{n-4}\}$ and $\{4, 2^{n-5}\}$ for $n \geq 10$, therefore giving a partial answer to a problem proposed by Schulte and Weiss in [Problems on polytopes, their groups, and realizations, Periodica Math. Hungarica 53(2006) 231-255].


title: Multi-Oriented Scene Text Detection via Corner Localization and Region Segmentation

abstract: Previous deep learning based state-of-the-art scene text detection methods can be roughly classified into two categories. The first category treats scene text as a type of general objects and follows general object detection paradigm to localize scene text by regressing the text box locations, but troubled by the arbitrary-orientation and large aspect ratios of scene text. The second one segments text regions directly, but mostly needs complex post processing. In this paper, we present a method that combines the ideas of the two types of methods while avoiding their shortcomings. We propose to detect scene text by localizing corner points of text bounding boxes and segmenting text regions in relative positions. In inference stage, candidate boxes are generated by sampling and grouping corner points, which are further scored by segmentation maps and suppressed by NMS. Compared with previous methods, our method can handle long oriented text naturally and doesn't need complex post processing. The experiments on ICDAR2013, ICDAR2015, MSRA-TD500, MLT and COCO-Text demonstrate that the proposed algorithm achieves better or comparable results in both accuracy and efficiency. Based on VGG16, it achieves an F-measure of 84.3% on ICDAR2015 and 81.5% on MSRA-TD500.


title: OhioState at SemEval-2018 Task 7: Exploiting Data Augmentation for Relation Classification in Scientific Papers using Piecewise Convolutional Neural Networks

abstract: We describe our system for SemEval-2018 Shared Task on Semantic Relation Extraction and Classification in Scientific Papers where we focus on the Classification task. Our simple piecewise convolution neural encoder performs decently in an end to end manner. A simple inter-task data augmentation signifi- cantly boosts the performance of the model. Our best-performing systems stood 8th out of 20 teams on the classification task on noisy data and 12th out of 28 teams on the classification task on clean data.


title: Multi-Segment Reconstruction Using Invariant Features

abstract: Multi-segment reconstruction (MSR) problem consists of recovering a signal from noisy segments with unknown positions of the observation windows. One example arises in DNA sequence assembly, which is typically solved by matching short reads to form longer sequences. Instead of trying to locate the segment within the sequence through pair-wise matching, we propose a new approach that uses shift-invariant features to estimate both the underlying signal and the distribution of the positions of the segments. Using the invariant features, we formulate the problem as a constrained nonlinear least-squares. The non-convexity of the problem leads to its sensitivity to the initialization. However, with clean data, we show empirically that for longer segment lengths, random initialization achieves exact recovery. Furthermore, we compare the performance of our approach to the results of expectation maximization and demonstrate that the new approach is robust to noise and computationally more efficient.


title: On the discrete analog of gamma-Lomax distribution: properties and applications

abstract: This article represents how certain types of blockades in any industrial (heavy industries) production, in particular, industrial strikes can be modeled with the proposed discrete probabilistic distribution as a baseline distribution. We considered the number of outbreaks of strikes in the coal mining industry, the vehicle manufacturing industry, and the transpose industry in the UK obtained from Consul (1989). We fitted those data sets with the proposed discrete gamma-Lomax distribution and compared the fit with the discrete generalized Pareto distribution (Consul, 1989). For this purpose, we explore the basic properties of the discrete gamma-Lomax distribution including but not limited to: cumulative distribution, survival, probability mass, quantile and hazard functions, genesis and rth-order moments; consider maximum likelihood estimation under the normal set up as well as under the censored data set scenario. It is observed that the newly proposed model can be useful to describe strikes arising from various types of industries.


title: Efficient nonparametric causal inference with missing exposure information

abstract: Missing exposure information is a very common feature of many observational studies. Here we study identifiability and efficient estimation of causal effects on vector outcomes, in such cases where treatment is unconfounded but partially missing. We consider a missing at random setting where missingness in treatment can depend not only on complex covariates, but also on post-treatment outcomes. We give a new identifying expression for average treatment effects in this setting, along with the efficient influence function for this parameter in a nonparametric model, which yields a nonparametric efficiency bound. We use this latter result to construct nonparametric estimators that are less sensitive to the curse of dimensionality than usual, e.g., by having faster rates of convergence than the complex nuisance estimators they rely on. Further we show that these estimators can be root-n consistent and asymptotically normal under weak nonparametric conditions, even when constructed using flexible machine learning. Finally we apply these results to the problem of causal inference with a partially missing instrumental variable.


title: Robust Target-relative Localization with Ultra-Wideband Ranging and Communication

abstract: In this paper we propose a method to achieve relative positioning and tracking of a target by a quadcopter using Ultra-wideband (UWB) ranging sensors, which are strategically installed to help retrieve both relative position and bearing between the quadcopter and target. To achieve robust localization for autonomous flight even with uncertainty in the speed of the target, two main features are developed. First, an estimator based on Extended Kalman Filter (EKF) is developed to fuse UWB ranging measurements with data from onboard sensors including inertial measurement unit (IMU), altimeters and optical flow. Second, to properly handle the coupling of the target's orientation with the range measurements, UWB based communication capability is utilized to transfer the target's orientation to the quadcopter. Experiment results demonstrate the ability of the quadcopter to control its position relative to the target autonomously in both cases when the target is static and moving.


title: Polarization as a probe of dusty environments around Type Ia supernovae: radiative transfer models for SN 2012dn

abstract: Geometry of circumstellar (CS) medium around supernovae (SNe) provides important diagnostics to understand the nature of their progenitors. In this paper, properties of CS dust around SN 2012dn, a super-Chandrasekhar candidate Type Ia supernova (SC-SN), have been studied through detailed three dimensional radiation transfer simulations. With the detected near-infrared excess from SN 2012dn, we show that it has a disk-like dusty CS environment whose mass is roughly consistent with a branch of an accreting white dwarf system (the single degenerate scenario). We show that a similar system should produce polarization signals up to $\sim 8$ \% in the $B$ band, depending on the viewing direction if polarimetric observation is performed. We predict that the maximum polarization is reached around $\sim 60$ days after the $B$-band maximum. We show that the temporal and wavelength dependence of the polarization signals, together with other unique features, can be easily distinguished from the interstellar polarization and intrinsic SN polarization. Indeed, the small polarization degree observed for normal Type Ia SNe (SNe Ia) can constrain a parameter space in the CS dust mass and distribution. We thus encourage multi-band polarimetric observations for SNe Ia, especially for outliers including SC-SNe for which some arguments for the single degenerate scenario exist but the polarization data are very rare so far.


title: On the Broadcast Routing Problem in Computer Networks

abstract: Given an undirected graph $G = (V, E)$, and a vertex $r\in V$, an $r$-acyclic orientation of $G$ is an orientation $OE$ of the edges of $G$ such that the digraph $OG = (V, OE)$ is acyclic and $r$ is the unique vertex with indegree equal to 0. For $w\in \mathbb{R}^E_+$, $k(G, w)$ is the value of the $w$-maximum packing of $r$-arborescences for all $r\in V$ and all $r$-acyclic orientations $OE$ of $G$. In this case, the Broadcast Routing (in Computers Networks) Problem (BRP) is to compute $k(G, w)$, by finding an optimal $r$ and an optimal $r$-acyclic orientation. BRP is a mathematical formulation of multipath broadcast routing in computer networks. In this paper, we provide a polynomial time algorithm to solve BRP in outerplanar graphs. Outerplanar graphs are encountered in many applications such as computational geometry, robotics, etc.


title: Muliti-scale regularity of axisymmetric Navier-Stokes equations

abstract: By applying the delicate \textit{a priori} estimates for the equations of $(\Phi,\Gamma)$, which is introduced in the previous work, we obtain some multi-scale regularity criteria of the swirl component $u^{\theta}$ for the 3D axisymmetric Navier-Stokes equations. In particularly, the solution $\mathbf{u}$ can be continued beyond the time $T$, provided that $u^{\theta}$ satiesfies $$ u^{\theta} \in L^{p}_{T}L^{q_{v}}_{v}L^{q_{h},w}_{h},~~\frac{2}{p}+\frac{1}{q_{v}}+\frac{2}{q_{h}}\leq 1, ~2<q_{h}\leq\infty,~\frac{1}{q_{v}}+\frac{2}{q_{h}}<1. $$


title: Measuring quantum discord using the most distinguishable steered states

abstract: Any two-qubit state can be represented, geometrically, as an ellipsoid with a certain size and a center located within the Bloch sphere of one of the qubits. Points of this ellipsoid represent the post-measurement states when the other qubit is measured. Based on the most demolition concept in the definition of quantum discord, we study the amount of demolition when the two post-measurement states, represented as two points on the steering ellipsoid, have the most distinguishability. We use trace distance as a measure of distinguishability and obtain the maximum distinguishability for some classes of states, analytically. Using the optimum measurement that gives the most distinguishable steered states, we extract quantum correlation of the state and compare the result with the quantum discord. It is shown that there are some important classes of states for which the most demolition happens exactly at the most distinguished steered points. Correlations gathered from the most distinguished post-measurement states provide a faithful and tight upper bound touching the quantum discord in most of the cases.


title: Group theory approach to unification of gravity with internal symmetry gauge interactions: I. Canonical electrogravity

abstract: The infinite group of deformed diffeomorphisms of the spacetime continuum is put into the basis of the gauge theory of gravity. This gives rise to some new ways for unification of gravity with other gauge interactions.