title: Predicting the properties of black holes merger remnants with Deep Neural Networks

abstract: We present the first estimation of the mass and spin magnitude of Kerr black holes resulting from the coalescence of binary black holes using a deep neural network. The network is trained on a dataset containing 80\% of the full publicly available catalog of numerical simulations of gravitational waves emission by binary black hole systems, including full precession effects for spinning binaries. The network predicts the remnant black holes mass and spin with an error less than 0.04\% and 0.3\% respectively for 90\% of the values in the non-precessing test dataset, it is 0.1\% and 0.3\% respectively in the precessing test dataset. When compared to existing fits in the LIGO algorithm software library, the network enables to reduce the remnant mass root mean square error to one half in the non-precessing case. In the precessing case, both remnant mass and spin mean square errors are decreased to one half, and the network corrects the bias observed in available fits.


title: On Compositionality in Neural Machine Translation

abstract: We investigate two specific manifestations of compositionality in Neural Machine Translation (NMT) : (1) Productivity - the ability of the model to extend its predictions beyond the observed length in training data and (2) Systematicity - the ability of the model to systematically recombine known parts and rules. We evaluate a standard Sequence to Sequence model on tests designed to assess these two properties in NMT. We quantitatively demonstrate that inadequate temporal processing, in the form of poor encoder representations is a bottleneck for both Productivity and Systematicity. We propose a simple pre-training mechanism which alleviates model performance on the two properties and leads to a significant improvement in BLEU scores.


title: Studies for New Experiments at the CERN M2 Beamline within "Physics Beyond Colliders": AMBER/COMPASS++, NA64mu, MuonE

abstract: The "Physics Beyond Colliders (PBC)" study explores fundamental physics opportunities at the CERN accelerator complex complementary to collider experiments. Three new collaborations aim to exploit the M2 beamline in the North Area with existing high-intensity muon and hadron beams, but also aspire to go beyond the current M2 capabilities with a RF-separated, high intensity hadron beam, under study. The AMBER/COMPASS++ collaboration proposes an ambitious program with a measurement of the proton radius with muon beams, as well as QCD-related studies from pion PDFs / Drell-Yan to cross section measurements for dark sector searches. Assuming feasibility of the RF-separated beam, the spectrum of strange mesons would enter a high precision era while kaon PDFs as well as nucleon TMDs would be accessible via Drell-Yan reactions. The NA64mu collaboration proposes to search for dark sector mediators such as a dark scalar A' or a hypothetical Z_mu using the M2 muon beam and complementing their on-going A' searches with electron beams. The MuonE collaboration intends to assess the hadronic component of the vacuum polarization via elastic mu-e scattering, the dominant uncertainty in the determination of (g-2)_mu. An overview of the three new experimental programs will be presented together with implications for the M2 beamline and the experimental area EHN2, based on the studies of the PBC "Conventional Beams" Working Group.


title: Scalable and Data Privacy Conserving Controller Tuning for Large-Scale Power Networks

abstract: The increasing share of renewable generation leads to new challenges in reliable power system operation, such as the rising volatility of power generation, which leads to time-varying dynamics and behavior of the system. To counteract the changing dynamics, we propose to adapt the parameters of existing controllers to the changing conditions. Doing so, however, is challenging, as large power systems often involve multiple subsystem operators, which, for safety and privacy reasons, do not want to exchange detailed information about their subsystems. Furthermore, centralized tuning of structured controllers for large-scale systems, such as power networks, is often computationally very challenging. For this reason, we present a hierarchical decentralized approach for controller tuning, which increases data security and scalability. The proposed method is based on the exchange of structured reduced models of subsystems, which conserves data privacy and reduces computational complexity. For this purpose, suitable methods for model reduction and model matching are introduced. Furthermore, we demonstrate how increased renewable penetration leads to time-varying dynamics on the IEEE 68 bus power system, which underlines the importance of the problem. Then, we apply the proposed approach on simulation studies to show its effectiveness. As shown, a similar system performance as with a centralized method can be obtained. Finally, we show the scalability of the approach on a large power system with more than 2500 states and about 1500 controller parameters.


title: Transmission Lines Positive Sequence Parameters Estimation and Instrument Transformers Calibration Based on PMU Measurement Error Model

abstract: Phasor Measurement Unit measurement data have been widely used in nowadays power system applications both in steady state and dynamic analysis. The performance of these applications running in utilities' energy management system depends heavily on an accurate positive sequence power system model. However, it is impractical to find this accurate model with transmission line parameters calculated directly with the PMU measurements due to ratio errors brought by instrument transformers and communication errors brought by PMUs. Therefore, a methodology is proposed in this paper to estimate the actual transmission lines parameters throughout the whole system and, at the same time, calibrate the corresponding instrument transformers. A PMU positive sequence measurement error model is proposed targeting at the aforementioned errors, which is applicable to both transposed and un-transposed transmission lines. A single line parameters estimation method is designed based on Least Squares Estimation and this error model. This method requires only one set of reference measurements and the accuracy can be propagated throughout the whole network along with the topology acquired by the introduced Edge-based Breadth-first Search algorithm. The IEEE 118-bus system and the Texas 2000-bus system are used to demonstrate the effectiveness and efficiency of the proposed method. The potential for deployment in reality is also discussed.


title: The effect of particle gas composition and boundary conditions on triboplasma generation: a computational study using the particle-in-cell method

abstract: Two dimensional particle in cell simulations of free charge creation by collisional ionization of C12 and C60 molecules immersed in plasma for the parameters of relevance to plasma gasification are presented. Our main findings are that (i) in uniform plasmas with smooth walls two optimal values which emerge for free electron production by collisional ionization (i.e. a most efficient discharge condition creation) are $C60:C12$ fractions of $10:90$ and $80:20$, (ii) in plasmas with rough walls, modelled by comb-like electric field at the boundary, the case of tangential electric field creates significant charge localization in C12+ and C60+ species, again creating most favorable discharge condition for tribo-electrically generated plasma. The numerical simulation results are discussed with reference to recent triboelectric plasma experiments and are corroborated by suitable analytical models.


title: Combinatorial list-decoding of Reed-Solomon codes beyond the Johnson radius

abstract: List-decoding of Reed-Solomon (RS) codes beyond the so called Johnson radius has been one of the main open questions since the work of Guruswami and Sudan. It is now known by the work of Rudra and Wootters, using techniques from high dimensional probability, that over large enough alphabets most RS codes are indeed list-decodable beyond this radius. In this paper we take a more combinatorial approach which allows us to determine the precise relation (up to the exact constant) between the decoding radius and the list size. We prove a generalized Singleton bound for a given list size, and conjecture that the bound is tight for most RS codes over large enough finite fields. We also show that the conjecture holds true for list sizes $2 \text{ and }3$, and as a by product show that most RS codes with a rate of at least $1/9$ are list-decodable beyond the Johnson radius. Lastly, we give the first explicit construction of such RS codes. The main tools used in the proof are a new type of linear dependency between codewords of a code that are contained in a small Hamming ball, and the notion of cycle space from Graph Theory. Both of them have not been used before in the context of list-decoding.


title: Metropolized Forest Recombination for Monte Carlo Sampling of Graph Partitions

abstract: We develop a new Markov chain on graph partitions that makes relatively global moves yet is computationally feasible to be used as the proposal in the Metropolis-Hastings method. Our resulting algorithm can be made reversible and able to sample from a specified measure on partitions. Both of these properties are critical to some important applications and computational Bayesian statistics in general. Our proposal chain modifies the recently developed method called Recombination (ReCom), which draws spanning trees on joined partitions and then randomly cuts them to repartition. We improve the computational efficiency by augmenting the state space from partitions to spanning forests. The extra information accelerates the computation of the forward and backward proposal probabilities. We demonstrate this method by sampling redistricting plans and find promising convergence results on several key observables of interest.


title: Statistical physics approaches to Unique Games

abstract: We show how two techniques from statistical physics can be adapted to solve a variant of the notorious Unique Games problem, potentially opening new avenues towards the Unique Games Conjecture. The variant, which we call Count Unique Games, is a promise problem in which the "yes" case guarantees a certain number of highly satisfiable assignments to the Unique Games instance. In the standard Unique Games problem, the "yes" case only guarantees at least one such assignment. We exhibit efficient algorithms for Count Unique Games based on approximating a suitable partition function for the Unique Games instance via (i) a zero-free region and polynomial interpolation, and (ii) the cluster expansion. We also show that a modest improvement to the parameters for which we give results would refute the Unique Games Conjecture.


title: Conformal field theory and the non-abelian $SU(2)_k$ chiral spin liquid

abstract: We construct a family of 1D and 2D long-range SU(2) spin models as parent Hamiltonians associated with infinite dimensional matrix product states that arise from simple current correlation functions in $SU(2)_k$ WZW models. Our results provide a conformal field theory foundation for recent proposals by Greiter and coauthors regarding the realization of non-abelian chiral spin liquids. We explain, in particular, how the symmetrization procedure for the amplitudes of the ground state wave function suggested in Ref. arXiv:0903.4547 originates from the conformal field theory description.


title: Aerial Access and Backhaul in mmWave B5G Systems: Performance Dynamics and Optimization

abstract: The use of unmanned aerial vehicle (UAV)-based communication in millimeter-wave (mmWave) frequencies to provide on-demand radio access is a promising approach to improve capacity and coverage in beyond-5G (B5G) systems. There are several design aspects to be addressed when optimizing for the deployment of such UAV base stations. As traffic demand of mobile users varies across time and space, dynamic algorithms that correspondingly adjust the UAV locations are essential to maximize performance. In addition to careful tracking of spatio-temporal user/traffic activity, such optimization needs to account for realistic backhaul constraints. In this work, we first review the latest 3GPP activities behind integrated access and backhaul system design, support for UAV base stations, and mmWave radio relaying functionality. We then compare static and mobile UAV-based communication options under practical assumptions on the mmWave system layout, mobility and clusterization of users, antenna array geometry, and dynamic backhauling. We demonstrate that leveraging the UAV mobility to serve moving users may improve the overall system performance even in the presence of backhaul capacity limitations.


title: Minimal Solvers for Rectifying from Radially-Distorted Conjugate Translations

abstract: This paper introduces minimal solvers that jointly solve for radial lens undistortion and affine-rectification using local features extracted from the image of coplanar translated and reflected scene texture, which is common in man-made environments. The proposed solvers accommodate different types of local features and sampling strategies, and three of the proposed variants require just one feature correspondence. State-of-the-art techniques from algebraic geometry are used to simplify the formulation of the solvers. The generated solvers are stable, small and fast. Synthetic and real-image experiments show that the proposed solvers have superior robustness to noise compared to the state of the art. The solvers are integrated with an automated system for rectifying imaged scene planes from coplanar repeated texture. Accurate rectifications on challenging imagery taken with narrow to wide field-of-view lenses demonstrate the applicability of the proposed solvers.


title: Verifying Visibility-Based Weak Consistency

abstract: Multithreaded programs generally leverage efficient and thread-safe concurrent objects like sets, key-value maps, and queues. While some concurrent-object operations are designed to behave atomically, each witnessing the atomic effects of predecessors in a linearization order, others forego such strong consistency to avoid complex control and synchronization bottlenecks. For example, contains (value) methods of key-value maps may iterate through key-value entries without blocking concurrent updates, to avoid unwanted performance bottlenecks, and consequently overlook the effects of some linearization-order predecessors. While such weakly-consistent operations may not be atomic, they still offer guarantees, e.g., only observing values that have been present. In this work we develop a methodology for proving that concurrent object implementations adhere to weak-consistency specifications. In particular, we consider (forward) simulation-based proofs of implementations against relaxed-visibility specifications, which allow designated operations to overlook some of their linearization-order predecessors, i.e., behaving as if they never occurred. Besides annotating implementation code to identify linearization points, i.e., points at which operations' logical effects occur, we also annotate code to identify visible operations, i.e., operations whose effects are observed; in practice this annotation can be done automatically by tracking the writers to each accessed memory location. We formalize our methodology over a general notion of transition systems, agnostic to any particular programming language or memory model, and demonstrate its application, using automated theorem provers, by verifying models of Java concurrent object implementations.


title: Understanding racial bias in health using the Medical Expenditure Panel Survey data

abstract: Over the years, several studies have demonstrated that there exist significant disparities in health indicators in the United States population across various groups. Healthcare expense is used as a proxy for health in algorithms that drive healthcare systems and this exacerbates the existing bias. In this work, we focus on the presence of racial bias in health indicators in the publicly available, and nationally representative Medical Expenditure Panel Survey (MEPS) data. We show that predictive models for care management trained using this data inherit this bias. Finally, we demonstrate that this inherited bias can be reduced significantly using simple mitigation techniques.


title: Deployment Architectures for Cyber-Physical Control Systems

abstract: We consider the problem of how to deploy a controller to a (networked) cyber-physical system (CPS). Controlling a CPS is an involved task, and synthesizing a controller to respect sensing, actuation, and communication constraints is only part of the challenge. In addition to controller synthesis, one should also consider how the controller will work in the CPS. Put another way, the cyber layer and its interaction with the physical layer need to be taken into account. In this work, we aim to bridge the gap between theoretical controller synthesis and practical CPS deployment. We adopt the system level synthesis (SLS) framework to synthesize a state-feedback controller and provide a deployment architecture for the standard SLS controller. Furthermore, we derive a new controller realization for open-loop stable systems and introduce four different architectures for deployment, ranging from fully centralized to fully distributed. Finally, we compare the trade-offs among them in terms of robustness, memory, computation, and communication overhead.


title: Model uncertainty estimation using the expectation maximization algorithm and a particle flow filter

abstract: Model error covariances play a central role in the performance of data assimilation methods applied to nonlinear state-space models. However, these covariances are largely unknown in most of the applications. A misspecification of the model error covariance has a strong impact on the computation of the posterior probability density function, leading to unreliable estimations and even to a total failure of the assimilation procedure. In this work, we propose the combination of the Expectation-Maximization algorithm (EM) with an efficient particle filter to estimate the model error covariance, using a batch of observations. Based on the EM algorithm principles, the proposed method encompasses two stages: the expectation stage, in which a particle filter is used with the present estimate of the model error covariance as given to find the probability density function that maximizes the likelihood, followed by a maximization stage in which the expectation under the probability density function found in the expectation step is maximized as a function of the elements of the model error covariance. This novel algorithm here presented combines the EM with a fixed point algorithm and does not require a particle smoother to approximate the posterior densities. We demonstrate that the new method accurately and efficiently solves the linear model problem. Furthermore, for the chaotic nonlinear Lorenz-96 model the method is stable even for observation error covariance 10 times larger than the estimated model error covariance matrix, and also that it is successful in high-dimensional situations where the dimension of the estimated matrix is 1600.


title: Universal finite size scaling around tricriticality between topologically ordered, SPT, and trivial phases

abstract: A quantum tricritical point is shown to exists in coupled time-reversal symmetry (TRS) broken Majorana chains. The tricriticality separates topologically ordered, symmetry protected topological (SPT), and trivial phases of the system. Here we demonstrate that the breaking of the TRS manifests itself in an emergence of a new dimensionless scale, $g = \alpha(\xi) B \sqrt{N}$, where $N$ is the system size, $B$ is a generic TRS breaking field, and $\alpha(\xi)$, with $\alpha(0)\equiv 1$, is a model-dependent function of the localization length, $\xi$, of boundary Majorana zero modes at the tricriticality. This scale determines the scaling of the finite size corrections around the tricriticality, which are shown to be {\it universal}, and independent of the nature of the breaking of the TRS. We show that the single variable scaling function, $f(w)$, $w\propto m N$, where $m$ is the excitation gap, that defines finite-size corrections to the ground state energy of the system around topological phase transition at $B=0$, becomes double-scaling, $f=f(w,g)$, at finite $B$. We realize TRS breaking through three different methods with completely different lattice details and find the same universal behavior of $f(w,g)$. In the critical regime, $m=0$, the function $f(0,g)$ is nonmonotonic, and reproduces the Ising conformal field theory scaling only in limits $g=0$ and $g\rightarrow \infty$. The obtained result sets a scale of $N \gg 1/(\alpha B)^2$ for the system to reach the thermodynamic limit in the presence of the TRS breaking. We derive the effective low-energy theory describing the tricriticality and analytically find the asymptotic behavior of the finite-size scaling function. Our results show that the boundary entropy around the tricriticality is also a universal function of $g$ at $m=0$.


title: Local on-surface radiation condition for multiple scattering of waves from convex obstacles

abstract: We propose a novel on-surface radiation condition to approximate the outgoing solution to the Helmholtz equation in the exterior of several impenetrable convex obstacles. Based on a local approximation of the Dirichlet-to-Neumann operator and a local formula for wave propagation, this new method simultaneously accounts for the outgoing behavior of the solution as well as the reflections arising from the multiple obstacles. The method involves tangential derivatives only, avoiding the use of integration over the surfaces of the obstacles. As a consequence, the method leads to sparse matrices and O(N) complexity. Numerical results are presented to illustrate the performance of the proposed method. Possible improvements and extensions are also discussed.


title: On the Form of Solutions of Fuchsian differential Equations with n regular singular Points

abstract: The form of the coefficients of power series expressions corresponding to solutions of Fuchsian differential equations (or their associated degenerated confluent forms) with n regular singular points is determined by solving the corresponding n-term recurrence relations in full generality. Some important special cases are discussed in which the solutions coincide with special functions of mathematical physics.


title: Can the Elliptic Billiard Still Surprise Us?

abstract: Can any secrets still be shed by that much studied, uniquely integrable, Elliptic Billiard? Starting by examining the family of 3-periodic trajectories and the loci of their Triangular Centers, one obtains a beautiful and variegated gallery of curves: ellipses, quartics, sextics, circles, and even a stationary point. Secondly, one notices this family conserves an intriguing ratio: Inradius-to-Circumradius. In turn this implies three conservation corollaries: (i) the sum of bounce angle cosines, (ii) the product of excentral cosines, and (iii) the ratio of excentral-to-orbit areas. Monge's Orthoptic Circle's close relation to 4-periodic Billiard trajectories is well-known. Its geometry provided clues with which to generalize 3-periodic invariants to trajectories of an arbitrary number of edges. This was quite unexpected. Indeed, the Elliptic Billiard did surprise us!


title: Latent Theme Dictionary Model for Finding Co-occurrent Patterns in Process Data

abstract: Process data, temporally ordered categorical observations, are of recent interest due to its increasing abundance and the desire to extract useful information. A process is a collection of time-stamped events of different types, recording how an individual behaves in a given time period. The process data are too complex in terms of size and irregularity for the classical psychometric models to be applicable, at least directly, and, consequently, it is desirable to develop new ways for modeling and analysis. We introduce herein a latent theme dictionary model (LTDM) for processes that identifies co-occurrent event patterns and individuals with similar behavioral patterns. Theoretical properties are established under certain regularity conditions for the likelihood based estimation and inference. A non-parametric Bayes LTDM algorithm using the Markov Chain Monte Carlo method is proposed for computation. Simulation studies show that the proposed approach performs well in a range of situations. The proposed method is applied to an item in the 2012 Programme for International Student Assessment with interpretable findings.


title: The Next Generation Very Large Array

abstract: The next generation Very Large Array (ngVLA) is a transformational radio observatory being designed by the U.S. National Radio Astronomy Observatory (NRAO). It will provide order of magnitude improvements in sensitivity, resolution, and uv coverage over the current Jansky Very Large Array (VLA) at ~1.2-50 GHz and extend the frequency range up to 70-115 GHz. This document is a white paper written by members of the Canadian community for the 2020 Long Range Plan panel, which will be making recommendations on Canada's future directions in astronomy. Since Canadians have been historically major users of the VLA and have been valued partners with NRAO for ALMA, Canada's participation in ngVLA is welcome. Canadians have been actually involved in ngVLA discussions for the past five years, and have played leadership roles in the ngVLA Science and Technical Advisory Councils. Canadian technologies are also very attractive for the ngVLA, in particular our designs for radio antennas, receivers, correlates, and data archives, and our industrial capacities to realize them. Indeed, the Canadian designs for the ngVLA antennas and correlator/beamformer are presently the baseline models for the project. Given the size of Canada's radio community and earlier use of the VLA (and ALMA), we recommend Canadian participation in the ngVLA at the 7% level. Such participation would be significant enough to allow Canadian leadership in gVLA's construction and usage. Canada's participation in ngVLA should not preclude its participation in SKA; access to both facilities is necessary to meet Canada's radio astronomy needs. Indeed, ngVLA will fill the gap between those radio frequencies observable with the SKA and ALMA at high sensitivities and resolutions. Canada's partnership in ngVLA will give it access to cutting-edge facilities together covering approximately three orders of magnitude in frequency.


title: Logjams are not jammed: measurements of log motions in Big Creek, Idaho

abstract: Colloquially, a "logjam" indicates a kinematic arrest of movement. Taken literally, it refers to a type of dense accumulation of wood in rivers widely recognized as bestowing numerous biological and physical benefits to the system but also present serious hazards to infrastructure. Despite this, no in-situ field measurements have assessed the degree of arrest in a naturally-formed logjam. Using time-lapse photography, repeat total station surveys and water level loggers, we provide an unprecedented perspective on the evolution of a logjam in central Idaho. Despite the namesake, we find that the logjam is not jammed. The ensemble of logs progressively deforms in response to shear and buoyant lift of flowing water, modulated by the rising limb, peak and falling limb of the snowmelt hydrograph. As water rises and log drag against the bed and banks decreases, they collectively translate downstream, generating a heterogeneous pattern of deformation. As streamflow recedes and the logs reconnect with the bed and banks, the coherent deformation pattern degrades as logs settle opportunistically amongst their neighbors. Field observations of continuous movement at a low rate are qualitatively similar to creep and clogging, behaviors that are common to a wide class of disordered materials. These similarities open the possibility to inform future studies of environmental clogging, wood-laden flows, logjams, hazard mitigation and the design of engineered logjams by bridging these practices with frontier research efforts in soft matter physics and granular rheology.


title: SAMM Long Videos: A Spontaneous Facial Micro- and Macro-Expressions Dataset

abstract: With the growth of popularity of facial micro-expressions in recent years, the demand for long videos with micro- and macro-expressions remains high. Extended from SAMM, a micro-expressions dataset released in 2016, this paper presents SAMM Long Videos dataset for spontaneous micro- and macro-expressions recognition and spotting. SAMM Long Videos dataset consists of 147 long videos with 343 macro-expressions and 159 micro-expressions. The dataset is FACS-coded with detailed Action Units (AUs). We compare our dataset with Chinese Academy of Sciences Macro-Expressions and Micro-Expressions (CAS(ME)2) dataset, which is the only available fully annotated dataset with micro- and macro-expressions. Furthermore, we preprocess the long videos using OpenFace, which includes face alignment and detection of facial AUs. We conduct facial expression spotting using this dataset and compare it with the baseline of MEGC III. Our spotting method outperformed the baseline result with F1-score of 0.3299.


title: Role of Coherent Eddies in Potential Vorticity Transport in Two-layer Quasigeostrophic Turbulence

abstract: The transport by materially coherent eddies is studied in a two-layer quasigeostrophic model of geophysical turbulence. The coherent eddies are identified by closed contours of the Lagrangian-averaged vorticity deviation obtained from Lagrangian particles advected by the flow. A series of flow regimes with different bottom friction strengths are considered---it is found that coherent eddies become more prevalent and longer-lasting as the strength of the bottom drag increases. These coherent eddies, with average core radius close to the deformation radius, propagate zonally with speeds close to the long baroclinic Rossby wave speed and meridionally with a preference for cyclones to propagate poleward and anticyclones to propagate equatorward. The meridional propagation preference of the coherent eddies gives rise to a systematic upgradient potential vorticity (PV) transport, which is in the opposite direction as the background PV transport and not captured by standard Lagrangian diffusivity estimates. The upgradient PV transport by coherent eddy cores is less than 10% of the total PV transport, but the PV transport by the periphery flow induced by the PV inside coherent eddies is significant and downgradient. This clarifies the distinct roles of the trapping and stirring effect of coherent eddies in PV transport in geophysical turbulence.


title: Multilateration of Random Networks with Community Structure

abstract: The minimal number of nodes required to multilaterate a network endowed with geodesic distance (i.e., to uniquely identify all nodes based on shortest path distances to the selected nodes) is called its metric dimension. This quantity is related to a useful technique for embedding graphs in low-dimensional Euclidean spaces and representing the nodes of a graph numerically for downstream analyses such as vertex classification via machine learning. While metric dimension has been studied for many kinds of graphs, its behavior on the Stochastic Block Model (SBM) ensemble has not. The simple community structure of graphs in this ensemble make them interesting in a variety of contexts. Here we derive probabilistic bounds for the metric dimension of random graphs generated according to the SBM, and describe algorithms of varying complexity to find---with high probability---subsets of nodes for multilateration. Our methods are tested on SBM ensembles with parameters extracted from real-world networks. We show that our methods scale well with increasing network size as compared to the state-of-the-art Information Content Heuristic algorithm for metric dimension approximation.


title: Generalized Self-concordant Hessian-barrier algorithms

abstract: Many problems in statistical learning, imaging, and computer vision involve the optimization of a non-convex objective function with singularities at the boundary of the feasible set. For such challenging instances, we develop a new interior-point technique building on the Hessian-barrier algorithm recently introduced in Bomze, Mertikopoulos, Schachinger and Staudigl, [SIAM J. Opt. 2019 29(3), pp. 2100-2127], where the Riemannian metric is induced by a generalized self-concordant function. This class of functions is sufficiently general to include most of the commonly used barrier functions in the literature of interior point methods. We prove global convergence to an approximate stationary point of the method, and in cases where the feasible set admits an easily computable self-concordant barrier, we verify worst-case optimal iteration complexity of the method. Applications in non-convex statistical estimation and $L^{p}$-minimization are discussed to given the efficiency of the method.


title: Counterexample-Guided Synthesis of Perception Models and Control

abstract: Recent advances in learning-based perception systems have led to drastic improvements in the performance of robotic systems like autonomous vehicles and surgical robots. These perception systems, however, are hard to analyze and errors in them can propagate to cause catastrophic failures. In this paper, we consider the problem of synthesizing safe and robust controllers for robotic systems which rely on complex perception modules for feedback. We propose a counterexample-guided synthesis framework that iteratively builds simple surrogate models of the complex perception module and enables us to find safe control policies. The framework uses a falsifier to find counterexamples, or traces of the systems that violate a safety property, to extract information that enables efficient modeling of the perception modules and errors in it. These models are then used to synthesize controllers that are robust to errors in perception. If the resulting policy is not safe, we gather new counterexamples. By repeating the process, we eventually find a controller which can keep the system safe even when there is a perception failure. We demonstrate our framework on two scenarios in simulation, namely lane keeping and automatic braking, and show that it generates controllers that are safe, as well as a simpler model of a deep neural network-based perception system that can provide meaningful insight into operations of the perception system.


title: Efficiency Improvement of Commercially Available MPPT Controllers Using Boost Converter

abstract: The use of charge controllers in photovoltaic generation systems increases the energy harvested significantly. The most efficient type of charge controller is the Maximum Power Point Tracking (MPPT) controller. The commercially available MPPT controller is limited to the input voltage of the solar panel and to the charging voltage of the battery making it inefficient to very low irradiance weather condition. In this paper, it is proposed to increase the efficiency of commercially available MPPT controllers by adding Power Management System that has boost converter and switching circuit that will operate at low irradiance weather conditions like cloudy, sunrise, and sunset. The boost converter is controlled by the PWM function of Arduino Uno to regulate the voltage required of the MPPT controller. The switching circuit will automatically change the direction of current depending on the output of the voltage sensors attached to the PV panel and input of MPPT controller. The results showed that there is 8.77% increase in average power to the system with Power Management System compared to the MPPT controller alone.


title: Statistical Inference in Mean-Field Variational Bayes

abstract: We conduct non-asymptotic analysis on the mean-field variational inference for approximating posterior distributions in complex Bayesian models that may involve latent variables. We show that the mean-field approximation to the posterior can be well-approximated relative to the Kullback-Leibler divergence discrepancy measure by a normal distribution whose center is the maximum likelihood estimator (MLE). In particular, our results imply that the center of the mean-field approximation matches the MLE up to higher-order terms and there is essentially no loss of efficiency in using it as a point estimator for the parameter in any regular parametric model with latent variables. We also propose a new class of variational weighted likelihood bootstrap (VWLB) methods for quantifying the uncertainty in the mean-field variational inference. The proposed VWLB can be viewed as a new sampling scheme that produces independent samples for approximating the posterior. Comparing with traditional sampling algorithms such Markov Chain Monte Carlo, VWLB can be implemented in parallel and is free of tuning.

categories: math.ST stat.AP stat.CO stat.ME stat.ML stat.TH


title: Beryllium Polyhydride Be4H8(H2)2 Synthesized at HP/HT

abstract: We report the XRD and Raman scattering measurements in combination with DFT calculations that reveal the formation of beryllium polyhydride Be4H8(H2)2 by laser heating Be/H2 mixture to above 1700 K at pressures between 5 GPa and 8 GPa. The Be4H8(H2)2 crystallizes in a P6_3/mmc structure and consists of corner-sharing BeH4 tetrahedrons and H2 molecules that are in an interstitial site. The Be4H8(H2)2 is stable at least to 14 GPa on compression and stable down to 4 GPa at room temperature. Our ab-initio calculations suggest that the Be4H8(H2)2 is a meta-stable phase of Be-H system.


title: Same Stats, Different Graphs: Exploring the Space of Graphs in Terms of Graph Properties

abstract: Data analysts commonly utilize statistics to summarize large datasets. While it is often sufficient to explore only the summary statistics of a dataset (e.g., min/mean/max), Anscombe's Quartet demonstrates how such statistics can be misleading. We consider a similar problem in the context of graph mining. To study the relationships between different graph properties, we examine low-order non-isomorphic graphs and provide a simple visual analytics system to explore correlations across multiple graph properties. However, for larger graphs, studying the entire space quickly becomes intractable. We use different random graph generation methods to further look into the distribution of graph properties for higher order graphs and investigate the impact of various sampling methodologies. We also describe a method for generating many graphs that are identical over a number of graph properties and statistics yet are clearly different and identifiably distinct.


title: BAS: An Answer Selection Method Using BERT Language Model

abstract: In recent years, Question Answering systems have become more popular and widely used by users. Despite the increasing popularity of these systems, the their performance is not even sufficient for textual data and requires further research. These systems consist of several parts that one of them is the Answer Selection component. This component detects the most relevant answer from a list of candidate answers. The methods presented in previous researches have attempted to provide an independent model to undertake the answer-selection task. An independent model cannot comprehend the syntactic and semantic features of questions and answers with a small training dataset. To fill this gap, language models can be employed in implementing the answer selection part. This action enables the model to have a better understanding of the language in order to understand questions and answers better than previous works. In this research, we will present the "BAS" (BERT Answer Selection) that uses the BERT language model to comprehend language. The empirical results of applying the model on the TrecQA Raw, TrecQA Clean, and WikiQA datasets demonstrate that using a robust language model such as BERT can enhance the performance. Using a more robust classifier also enhances the effect of the language model on the answer selection component. The results demonstrate that language comprehension is an essential requirement in natural language processing tasks such as answer-selection.


title: Closing the Reality Gap with Unsupervised Sim-to-Real Image Translation

abstract: Deep learning approaches have become the standard solution to many problems in computer vision and robotics, but obtaining sufficient training data in high enough quality is challenging, as human labor is error prone, time consuming, and expensive. Solutions based on simulation have become more popular in recent years, but the gap between simulation and reality is still a major issue. In this paper, we introduce a novel method for augmenting synthetic image data through unsupervised image-to-image translation by applying the style of real world images to simulated images with open source frameworks. The generated dataset is combined with conventional augmentation methods and is then applied to a neural network model running in real-time on autonomous soccer robots. Our evaluation shows a significant improvement compared to models trained on images generated entirely in simulation.


title: The effects of disk self-gravity and radiative cooling on the formation of gaps and spirals by young planets

abstract: We have carried out two-dimensional hydrodynamical simulations to study the effects of disk self-gravity and radiative cooling on the formation of gaps and spirals. (1) With disk self-gravity included, we find stronger, more tightly-wound spirals and deeper gaps in more massive disks. The deeper gaps are due to the larger Angular Momentum Flux (AMF) of the waves excited in more massive disks, as expected from the linear theory. The position of the secondary gap does not change, provided that the disk is not extremely massive ($Q \gtrsim 2$). (2) With radiative cooling included, the excited spirals become monotonically more open (less tightly-wound) as the disk's cooling timescale increases. On the other hand, the amplitude and strength of the spirals decrease when the cooling time increases from a small value to $\sim 1/\Omega$, but then the amplitude starts to increase again when the cooling time continues to increase. This indicates that radiative dissipation becomes important for waves with $T_{cool}\sim$ 1. Consequently, the induced primary gap is narrower and the secondary gap becomes significantly shallower when the cooling time becomes $\sim 1/\Omega$. When the secondary gap is present, the position of it moves to the inner disk from the fast cooling cases to the slow cooling cases. The dependence of gap properties on the cooling timescale (e.g. in AS 209) provides a new way to constrain the disk optical depth and thus disk surface density.


title: Real-Time Sensor Anomaly Detection and Recovery in Connected Automated Vehicle Sensors

abstract: In this paper we propose a novel observer-based method to improve the safety and security of connected and automated vehicle (CAV) transportation. The proposed method combines model-based signal filtering and anomaly detection methods. Specifically, we use adaptive extended Kalman filter (AEKF) to smooth sensor readings of a CAV based on a nonlinear car-following motion model. Under the assumption of a car-following model, the subject vehicle utilizes its leading vehicle's information to detect sensor anomalies by employing previously-trained One Class Support Vector Machine (OCSVM) models. This approach allows the AEKF to estimate the state of a vehicle not only based on the vehicle's location and speed, but also by taking into account the state of the surrounding traffic. A communication time delay factor is considered in the car-following model to make it more suitable for real-world applications. Our experiments show that compared with the AEKF with a traditional $\chi^2$-detector, our proposed method achieves a better anomaly detection performance. We also demonstrate that a larger time delay factor has a negative impact on the overall detection performance.


title: Real-time Funnel Generation for Restricted Motion Planning

abstract: In autonomous systems, a motion planner generates reference trajectories which are tracked by a low-level controller. For safe operation, the motion planner should account for inevitable controller tracking error when generating avoidance trajectories. In this article we present a method for generating provably safe tracking error bounds, while reducing over-conservatism that exists in existing methods. We achieve this goal by restricting possible behaviors for the motion planner. We provide an algebraic method based on sum-of-squares programming to define restrictions on the motion planner and find small bounds on the tracking error. We demonstrate our method on two case studies and show how we can integrate the method into already developed motion planning techniques. Results suggest that our method can provide acceptable tracking error wherein previous work were not applicable.


title: Speaker-invariant Affective Representation Learning via Adversarial Training

abstract: Representation learning for speech emotion recognition is challenging due to labeled data sparsity issue and lack of gold standard references. In addition, there is much variability from input speech signals, human subjective perception of the signals and emotion label ambiguity. In this paper, we propose a machine learning framework to obtain speech emotion representations by limiting the effect of speaker variability in the speech signals. Specifically, we propose to disentangle the speaker characteristics from emotion through an adversarial training network in order to better represent emotion. Our method combines the gradient reversal technique with an entropy loss function to remove such speaker information. Our approach is evaluated on both IEMOCAP and CMU-MOSEI datasets. We show that our method improves speech emotion classification and increases generalization to unseen speakers.


title: Another method to solve the grasshopper problem (the International Mathematical Olympiad)

abstract: The 6th problem of the 50th International Mathematical Olympiad (IMO), held in Germany, 2009, is called 'the grasshopper problem'. To this problem Kos developed theory from unique viewpoints by reference of Noga Alon's combinatorial Nullstellensatz. We have tried to solve this problem by an original method inspired by a polynomial function that Kos defined, then examined for n=3, 4 and 5. For almost cases the claim of this problem follows, but there remains imperfection due to 'singularity'.


title: Scalable Deep Generative Relational Models with High-Order Node Dependence

abstract: We propose a probabilistic framework for modelling and exploring the latent structure of relational data. Given feature information for the nodes in a network, the scalable deep generative relational model (SDREM) builds a deep network architecture that can approximate potential nonlinear mappings between nodes' feature information and the nodes' latent representations. Our contribution is two-fold: (1) We incorporate high-order neighbourhood structure information to generate the latent representations at each node, which vary smoothly over the network. (2) Due to the Dirichlet random variable structure of the latent representations, we introduce a novel data augmentation trick which permits efficient Gibbs sampling. The SDREM can be used for large sparse networks as its computational cost scales with the number of positive links. We demonstrate its competitive performance through improved link prediction performance on a range of real-world datasets.


title: Spectral Representations of Graphons in Very Large Network Systems Control

abstract: Graphon-based control has recently been proposed and developed to solve control problems for dynamical systems on networks which are very large or growing without bound (see Gao and Caines, CDC 2017, CDC 2018). In this paper, spectral representations, eigenfunctions and approximations of graphons, and their applications to graphon-based control are studied. First, spectral properties of graphons are presented and then approximations based on Fourier approximated eigenfunctions are analyzed. Within this framework, two classes of graphons with simple spectral representations are given. Applications to graphon-based control analysis are next presented; in particular, the controllability of systems distributed over very large networks is expressed in terms of the properties of the corresponding graphon dynamical systems. Moreover, spectral analysis based upon real-world network data is presented, which demonstrates that low-dimensional spectral approximations of networks are possible. Finally, an initial, exploratory investigation of the utility of the spectral analysis methodology in graphon systems control to study the control of epidemic spread is presented.


title: Verification and Parameter Synthesis for Stochastic Systems using Optimistic Optimization

abstract: We present an algorithm for formal verification and parameter synthesis of continuous state-space Markov chains. This class of problems captures the design and analysis of a wide variety of autonomous and cyber-physical systems defined by nonlinear and black-box modules. In order to solve these problems, one has to maximize certain probabilistic objective functions overall choices of initial states and parameters. In this paper, we identify the assumptions that make it possible to view this problem as a multi-armed bandit problem. Based on this fresh perspective, we propose an algorithm (HOO-MB) for solving the problem that carefully instantiates an existing bandit algorithm -- Hierarchical Optimistic Optimization -- with appropriate parameters. As a consequence, we obtain theoretical regret bounds on sample efficiency of our solution that depends on key problem parameters like smoothness, near-optimality dimension, and batch size. The batch size parameter enables us to strike a balance between the sample efficiency and the memory usage of the algorithm. Our experiments, using the tool HooVer, suggest that the approach scales to realistic-sized problems and is often more sample-efficient compared to PlasmaLab -- a leading tool for verification of stochastic systems. Specifically, HooVer has distinct advantages in analyzing models in which the objective function has sharp slopes. In addition, HooVer shows promising behavior in parameter synthesis for a linear quadratic regulator (LQR) example.


title: Online tuning and light source control using a physics-informed Gaussian process Adi

abstract: Operating large-scale scientific facilities often requires fast tuning and robust control in a high dimensional space. In this paper we introduce a new physics-informed optimization algorithm based on Gaussian process regression. Our method takes advantage of the existing domain knowledge in the form of realizations of a physics model of the observed system. We have applied a physics-informed Gaussian Process method experimentally at the SPEAR3 storage ring to demonstrate online accelerator optimization. This method outperforms Gaussian Process trained on data as well as the standard approach routinely used for operation, in terms of convergence speed and optimal point. The proposed method could be applicable to automatic tuning and control of other complex systems, without a prerequisite for any observed data.


title: A Girsanov type representation of quadratic-exponential cost functionals for linear quantum stochastic systems

abstract: This paper is concerned with multimode open quantum harmonic oscillators and quadratic-exponential functionals (QEFs) as quantum risk-sensitive performance criteria. Such systems are described by linear quantum stochastic differential equations driven by multichannel bosonic fields. We develop a finite-horizon expansion for the system variables using the eigenbasis of their two-point commutator kernel with noncommuting position-momentum pairs as coefficients. This quantum Karhunen-Loeve expansion is used in order to obtain a Girsanov type representation for the quadratic-exponential functions of the system variables. This representation is valid regardless of a particular system-field state and employs the averaging over an auxiliary classical Gaussian random process whose covariance operator is defined in terms of the quantum commutator kernel. We use this representation in order to relate the QEF to the moment-generating functional of the system variables. This result is also specified for the invariant multipoint Gaussian quantum state when the oscillator is driven by vacuum fields.

categories: quant-ph cs.SY eess.SY math-ph math.MP math.OC math.PR


title: Motivic Galois coaction and one-loop Feynman graphs

abstract: Following the work of Brown, we can canonically associate a family of motivic periods -- called the motivic Feynman amplitude -- to any convergent Feynman integral, viewed as a function of the kinematic variables. The motivic Galois theory of motivic Feynman amplitudes provides an organizing principle, as well as strong constraints, on the space of amplitudes in general, via Brown's "small graphs principle". This serves as motivation for explicitly computing the motivic Galois action, or, dually, the coaction of the Hopf algebra of functions on the motivic Galois group. In this paper, we study the motivic Galois coaction on the motivic Feynman amplitudes associated to one-loop Feynman graphs. We study the associated variations of mixed Hodge structures, and provide an explicit formula for the coaction on the four-edge cycle graph -- the box graph -- with non-vanishing generic kinematics, which leads to a formula for all one-loop graphs with non-vanishing generic kinematics in four-dimensional space-time. We also show how one computes the coaction in some degenerate configurations -- when defining the motive of the graph requires blowing up the underlying family of varieties -- on the example of the three-edge cycle graph.


title: Limitations of the Hyperplane Separation Technique for Bounding the Extension Complexity of Polytopes

abstract: We illustrate the limitations of the hyperplane separation bound, a non-combinatorial lower bound on the extension complexity of a polytope. Most notably, this bounding technique is used by Rothvo{\ss} (J ACM 64.6:41, 2017) to establish an exponential lower bound for the perfect matching polytope. We point out that the technique is sensitive to the particular choice of slack matrix. For the canonical slack matrices of the spanning tree polytope and the completion time polytope, we show that the lower bounds produced by the hyperplane separation method are trivial. These bounds may, however, be strengthened by normalizing rows and columns of the slack matrices.


title: embComp: Visual Interactive Comparison of Vector Embeddings

abstract: This paper introduces embComp, a novel approach for comparing two embeddings that capture the similarity between objects, such as word and document embeddings. We survey scenarios where comparing these embedding spaces is useful. From those scenarios, we derive common tasks, introduce visual analysis methods that support these tasks, and combine them into a comprehensive system. One of embComp's central features are overview visualizations that are based on metrics for measuring differences in the local structure around objects. Summarizing these local metrics over the embeddings provides global overviews of similarities and differences. Detail views allow comparison of the local structure around selected objects and relating this local information to the global views. Integrating and connecting all of these components, embComp supports a range of analysis workflows that help understand similarities and differences between embedding spaces. We assess our approach by applying it in several use cases, including understanding corpora differences via word vector embeddings, and understanding algorithmic differences in generating embeddings.


title: Physics driven reduced order model for real time blood flow simulations

abstract: Predictive modeling of blood flow and pressure have numerous applications ranging from non-invasive assessment of functional significance of disease to planning invasive procedures. While several such predictive modeling techniques have been proposed, their use in the clinic has been limited due in part to the significant time required to perform virtual interventions and compute the resultant changes in hemodynamic conditions. We propose a fast hemodynamic assessment method based on first constructing an exploration space of geometries, tailored to each patient, and subsequently building a physics driven reduced order model in this space. We demonstrate that this method can predict fractional flow reserve derived from coronary computed tomography angiography in response to changes to a patient-specific lumen geometry in real time while achieving high accuracy when compared to computational fluid dynamics simulations. We validated this method on over 1300 patients that received a coronary CT scan and demonstrated a correlation coefficient of 0.98 with an error of 0.005 +- 0.015 (95% confidence interval: (-0.020, 0.031)) as compared to three-dimensional blood flow calculations.

categories: math.NA cs.NA physics.comp-ph physics.med-ph


title: The generalization error of max-margin linear classifiers: Benign overfitting and high dimensional asymptotics in the overparametrized regime

abstract: Modern machine learning classifiers often exhibit vanishing classification error on the training set. They achieve this by learning nonlinear representations of the inputs that maps the data into linearly separable classes. Motivated by these phenomena, we revisit high-dimensional maximum margin classification for linearly separable data. We consider a stylized setting in which data $(y_i,{\boldsymbol x}_i)$, $i\le n$ are i.i.d. with ${\boldsymbol x}_i\sim\mathsf{N}({\boldsymbol 0},{\boldsymbol \Sigma})$ a $p$-dimensional Gaussian feature vector, and $y_i \in\{+1,-1\}$ a label whose distribution depends on a linear combination of the covariates $\langle {\boldsymbol \theta}_*,{\boldsymbol x}_i \rangle$. While the Gaussian model might appear extremely simplistic, universality arguments can be used to show that the results derived in this setting also apply to the output of certain nonlinear featurization maps. We consider the proportional asymptotics $n,p\to\infty$ with $p/n\to \psi$, and derive exact expressions for the limiting generalization error. We use this theory to derive two results of independent interest: $(i)$ Sufficient conditions on $({\boldsymbol \Sigma},{\boldsymbol \theta}_*)$ for `benign overfitting' that parallel previously derived conditions in the case of linear regression; $(ii)$ An asymptotically exact expression for the generalization error when max-margin classification is used in conjunction with feature vectors produced by random one-layer neural networks.


title: Compositional Generalization with Tree Stack Memory Units

abstract: We study compositional generalization, viz., the problem of zero-shot generalization to novel compositions of concepts in a domain. Standard neural networks fail to a large extent on compositional learning. We propose Tree Stack Memory Units (Tree-SMU) to enable strong compositional generalization. Tree-SMU is a recursive neural network with Stack Memory Units (\SMU s), a novel memory augmented neural network whose memory has a differentiable stack structure. Each SMU in the tree architecture learns to read from its stack and to write to it by combining the stacks and states of its children through gating. The stack helps capture long-range dependencies in the problem domain, thereby enabling compositional generalization. Additionally, the stack also preserves the ordering of each node's descendants, thereby retaining locality on the tree. We demonstrate strong empirical results on two mathematical reasoning benchmarks. We use four compositionality tests to assess the generalization performance of Tree-SMU and show that it enables accurate compositional generalization compared to strong baselines such as Transformers and Tree-LSTMs.


title: Being Optimistic to Be Conservative: Quickly Learning a CVaR Policy

abstract: While maximizing expected return is the goal in most reinforcement learning approaches, risk-sensitive objectives such as conditional value at risk (CVaR) are more suitable for many high-stakes applications. However, relatively little is known about how to explore to quickly learn policies with good CVaR. In this paper, we present the first algorithm for sample-efficient learning of CVaR-optimal policies in Markov decision processes based on the optimism in the face of uncertainty principle. This method relies on a novel optimistic version of the distributional Bellman operator that moves probability mass from the lower to the upper tail of the return distribution. We prove asymptotic convergence and optimism of this operator for the tabular policy evaluation case. We further demonstrate that our algorithm finds CVaR-optimal policies substantially faster than existing baselines in several simulated environments with discrete and continuous state spaces.


title: On the Measure of Intelligence

abstract: To make deliberate progress towards more intelligent and more human-like artificial systems, we need to be following an appropriate feedback signal: we need to be able to define and evaluate intelligence in a way that enables comparisons between two systems, as well as comparisons with humans. Over the past hundred years, there has been an abundance of attempts to define and measure intelligence, across both the fields of psychology and AI. We summarize and critically assess these definitions and evaluation approaches, while making apparent the two historical conceptions of intelligence that have implicitly guided them. We note that in practice, the contemporary AI community still gravitates towards benchmarking intelligence by comparing the skill exhibited by AIs and humans at specific tasks such as board games and video games. We argue that solely measuring skill at any given task falls short of measuring intelligence, because skill is heavily modulated by prior knowledge and experience: unlimited priors or unlimited training data allow experimenters to "buy" arbitrary levels of skills for a system, in a way that masks the system's own generalization power. We then articulate a new formal definition of intelligence based on Algorithmic Information Theory, describing intelligence as skill-acquisition efficiency and highlighting the concepts of scope, generalization difficulty, priors, and experience. Using this definition, we propose a set of guidelines for what a general AI benchmark should look like. Finally, we present a benchmark closely following these guidelines, the Abstraction and Reasoning Corpus (ARC), built upon an explicit set of priors designed to be as close as possible to innate human priors. We argue that ARC can be used to measure a human-like form of general fluid intelligence and that it enables fair general intelligence comparisons between AI systems and humans.


title: The GALAH Survey: Chemically tagging the Fimbulthul stream to the globular cluster $\omega$ Centauri

abstract: Using kinematics from Gaia and the large elemental abundance space of the second data release of the GALAH survey, we identify two new members of the Fimbulthul stellar stream, and chemically tag them to massive, multi-metallic globular cluster $\omega$ Centauri. Recent analysis of the second data release of Gaia had revealed the Fimbulthul stellar stream in the halo of the Milky Way. It had been proposed that the stream is associated with the $\omega$ Centauri, but this proposition relied exclusively upon the kinematics and metallicities of the stars to make the association. In this work, we find our two new members of the stream to be metal-poor stars that are enhanced in sodium and aluminium, typical of second population globular cluster stars, but not otherwise seen in field stars. Furthermore, the stars share the s-process abundance pattern seen in $\omega$ Centauri, which is rare in field stars. Apart from one star within 1.5 deg of $\omega$ Centauri, we find no other stars observed by GALAH spatially near $\omega$ Centauri or the Fimbulthul stream that could be kinematically and chemically linked to the cluster. Chemically tagging stars in the Fimbulthul stream to $\omega$ Centauri confirms the earlier work, and further links this tidal feature in the Milky Way halo to $\omega$ Centauri.


title: Expansive systems on lattices

abstract: We study expansive dynamical systems in the setting of distributive lattices and their automorphisms, the usual notion of expansiveness for a homeomorphism of a compact metric space being the particular case when the lattice is the topology of the phase space ordered by inclusion and the automorphism the one induced by the homeomorphism, mapping open sets to open sets. We prove in this context generalizations of Ma\~n\'e's Theorem and Utz's Theorem about the finite dimensionality of the phase space of an expansive system, and the finiteness of a space supporting a positively expansive homeomorphism, respectively. We also discuss the notion of entropy in this setting calculating it for the non-Hausdorff shifts.


title: Global regularity of solutions for the 3D non-resistive and non-diffusive MHD-Boussinesq system with axisymmetric data

abstract: In this paper, we will show that solutions of the three-dimensional non-resistive and non-diffusive MHD-Boussinesq system are globally regular if the initial data is axisymmetric and the swirl components of the velocity and the magnetic vorticity are zero. Our main result extends previous ones on the three-dimensional non-resistive MHD system and non-diffusive Boussinesq system, and the method used here can also be applied to the magnetic Rayleigh-B\'{e}nard convection system.


title: Dynamic Graph Embedding via LSTM History Tracking

abstract: Many real world networks are very large and constantly change over time. These dynamic networks exist in various domains such as social networks, traffic networks and biological interactions. To handle large dynamic networks in downstream applications such as link prediction and anomaly detection, it is essential for such networks to be transferred into a low dimensional space. Recently, network embedding, a technique that converts a large graph into a low-dimensional representation, has become increasingly popular due to its strength in preserving the structure of a network. Efficient dynamic network embedding, however, has not yet been fully explored. In this paper, we present a dynamic network embedding method that integrates the history of nodes over time into the current state of nodes. The key contribution of our work is 1) generating dynamic network embedding by combining both dynamic and static node information 2) tracking history of neighbors of nodes using LSTM 3) significantly decreasing the time and memory by training an autoencoder LSTM model using temporal walks rather than adjacency matrices of graphs which are the common practice. We evaluate our method in multiple applications such as anomaly detection, link prediction and node classification in datasets from various domains.


title: A Conceptual Blend Analysis of Physics Quantitative Literacy Reasoning Inventory Items

abstract: Mathematical reasoning flexibility across physics contexts is a desirable learning outcome of introductory physics, where the math world and physical world meet. Physics Quantitative Literacy (PQL) is a set of interconnected skills and habits of mind that support quantitative reasoning about the physical world. The Physics Inventory of Quantitative Literacy (PIQL), which we are currently refining and validating, assesses students proportional reasoning, covariational reasoning, and reasoning with signed quantities in physics contexts. In this paper, we apply a Conceptual Blending Theory analysis of two exemplar PIQL items to demonstrate how we are using this theory to help develop an instrument that represents the kind of blended reasoning that characterizes expertise in physics. A Conceptual Blending Theory analysis allows for assessment of hierarchical partially correct reasoning patterns, and thereby holds potential to map the emergence of mathematical reasoning flexibility throughout the introductory physics sequence.


title: Global symbolic calculus of pseudo-differential operators on homogeneous vector bundles

abstract: A symbolic calculus for a pseudo-differential operators acting on sections of a homogeneous vector bundle over a compact homogeneous space $G/H$ with compact $G$ and $H$ is developed. We realize the symbol of a pseudo-differential operator as a linear operator acting on corresponding irreducible unitary representations of $H$ valued in the algebra $C^\infty(G)$ of smooth functions. We write down how left invariant vector fields of $SU(2)$ act on the sections of homogeneous vector bundles associated to the fibration $\mathbb{T}\hookrightarrow SU(2)\rightarrow\mathbb{C}P^1$, which is known as the Hopf fibration. Lastly, we outline how functional calculus of a pseudo-differential operator can be computed using our calculus.


title: Global Entanglement Distribution with Multi-mode Non-Gaussian Operations

abstract: Non-Gaussian operations have been studied intensively in recent years due to their ability to enhance the entanglement of quantum states. However, most previous studies on such operations are carried out in a single-mode setting, even though in reality any quantum state contains multi-mode components in frequency space. Whilst there have been general frameworks developed for multi-mode photon subtraction (PS) and photon addition (PA), an important gap exists in that no framework has thus far been developed for multi-mode photon catalysis (PC). In this work we close that gap. We then apply our newly developed PC framework to the problem of continuous variable (CV) entanglement distribution via quantum-enabled satellites. Due to the high pulse rate envisioned for such systems, multi-mode effects will be to the fore in space-based CV deployments. After determining the entanglement distribution possible via multi-mode PC, we then compare our results with the entanglement distribution possible using multi-mode PS and PA. Our results show that multi-mode PC carried out at the transmitter is the superior non-Gaussian operation when the initial squeezing is below some threshold. When carried out at the receiver, multi-mode PC is found to be the superior non-Gaussian operation when the mean channel attenuation is above some threshold. Our new results should prove valuable for next-generation deployments of CV quantum-enabled satellites.


title: Properly colored short cycles in edge-colored graphs

abstract: Properly colored cycles in edge-colored graphs are closely related to directed cycles in oriented graphs. As an analogy of the well-known Caccetta-H\"{a}ggkvist Conjecture, we study the existence of properly colored cycles of bounded length in an edge-colored graph. We first prove that for all integers $s$ and $t$ with $t\geq s\geq2$, every edge-colored graph $G$ with no properly colored $K_{s,t}$ contains a spanning subgraph $H$ which admits an orientation $D$ such that every directed cycle in $D$ is a properly colored cycle in $G$. Using this result, we show that for $r\geq4$, if the Caccetta-H\"{a}ggkvist Conjecture holds , then every edge-colored graph of order $n$ with minimum color degree at least $n/r+2\sqrt{n}+1$ contains a properly colored cycle of length at most $r$. In addition, we also obtain an asymptotically tight total color degree condition which ensures a properly colored (or rainbow) $K_{s,t}$.


title: Review-based Question Generation with Adaptive Instance Transfer and Augmentation

abstract: Online reviews provide rich information about products and service, while it remains inefficient for potential consumers to exploit the reviews for fulfilling their specific information need. We propose to explore question generation as a new way of exploiting review information. One major challenge of this task is the lack of review-question pairs for training a neural generation model. We propose an iterative learning framework for handling this challenge via adaptive transfer and augmentation of the training instances with the help of the available user-posed question-answer data. To capture the aspect characteristics in reviews, the augmentation and generation procedures incorporate related features extracted via unsupervised learning. Experiments on data from 10 categories of a popular E-commerce site demonstrate the effectiveness of the framework, as well as the usefulness of the new task.


title: Benchmarking Simulated Robotic Manipulation through a Real World Dataset

abstract: We present a benchmark to facilitate simulated manipulation; an attempt to overcome the obstacles of physical benchmarks through the distribution of a real world, ground truth dataset. Users are given various simulated manipulation tasks with assigned protocols having the objective of replicating the real world results of a recorded dataset. The benchmark comprises of a range of metrics used to characterise the successes of submitted environments whilst providing insight into their deficiencies. We apply our benchmark to two simulation environments, PyBullet and V-Rep, and publish the results. All materials required to benchmark an environment, including protocols and the dataset, can be found at the benchmarks' website https://research.csiro.au/robotics/manipulation-benchmark/.


title: Modelling synchrotron and synchrotron self-Compton emission of gamma-ray burst afterglows from radio to very-high energies

abstract: Synchrotron radiation from a decelerating blastwave is a widely accepted model of radio to X-ray afterglow emission from gamma-ray bursts (GRBs). GeV gamma-ray emission detected by the Fermi Large Area Telescope (LAT) and the duration of which extends beyond the prompt gamma-ray emission phase, is also compatible with broad features of afterglow emission. We revisit the synchrotron self-Compton (SSC) emission model from a decelerating blastwave to fit multiwavelength data from three bright GRBs, namely GRB~190114C, GRB~130427A and GRB~090510. We constrain the afterglow model parameters using the simultaneous fit of the spectral energy distributions at different times and light curves at different frequencies for these bursts. We find that a constant density interstellar medium is favored for the short GRB~090510, while a wind-type environment is favored for the long GRB~130427A and GRB~190114C. The sub-TeV component in GRB~190114C detected by MAGIC is the SSC emission in our modelling. Furthermore we find that the SSC emission in the Thomson regime is adequate to fit the spectra and light curves of GRB~190114C. For the other two GRBs, lacking sub-TeV detection, the SSC emissions are also modeled in the Thomson regime. For the model parameters we have used, the $\gamma\gamma$ attenuation in the blastwave is negligible in the sub-TeV range compared to the redshift-dependent $\gamma\gamma$ attenuation in the extragalactic background light.


title: A Tale of Evil Twins: Adversarial Inputs versus Poisoned Models

abstract: Despite their tremendous success in a range of domains, deep learning systems are inherently susceptible to two types of manipulations: adversarial inputs -- maliciously crafted samples that deceive target deep neural network (DNN) models, and poisoned models -- adversely forged DNNs that misbehave on pre-defined inputs. While prior work has intensively studied the two attack vectors in parallel, there is still a lack of understanding about their fundamental connections: what are the dynamic interactions between the two attack vectors? what are the implications of such interactions for optimizing existing attacks? what are the potential countermeasures against the enhanced attacks? Answering these key questions is crucial for assessing and mitigating the holistic vulnerabilities of DNNs deployed in realistic settings. Here we take a solid step towards this goal by conducting the first systematic study of the two attack vectors within a unified framework. Specifically, (i) we develop a new attack model that jointly optimizes adversarial inputs and poisoned models; (ii) with both analytical and empirical evidence, we reveal that there exist intriguing "mutual reinforcement" effects between the two attack vectors -- leveraging one vector significantly amplifies the effectiveness of the other; (iii) we demonstrate that such effects enable a large design spectrum for the adversary to enhance the existing attacks that exploit both vectors (e.g., backdoor attacks), such as maximizing the attack evasiveness with respect to various detection methods; (iv) finally, we discuss potential countermeasures against such optimized attacks and their technical challenges, pointing to several promising research directions.


title: Memory in nonmonotonic stress relaxation of a granular system

abstract: We demonstrate experimentally that a granular packing of glass spheres is capable of storing memory of multiple strain states in the dynamic process of stress relaxation. Modeling the system as a non-interacting population of relaxing elements, we find that the functional form of the predicted relaxation requires a quantitative correction which grows in severity with each additional memory and is suggestive of interactions between elements. Our findings have implications for the broad class of soft matter systems that display memory and anomalous relaxation.


title: Almost-sure enhanced dissipation and uniform-in-diffusivity exponential mixing for advection-diffusion by stochastic Navier-Stokes

abstract: We study the mixing and dissipation properties of the advection-diffusion equation with diffusivity $0 < \kappa \ll 1$ and advection by a class of random velocity fields on $\mathbb T^d$, $d=\{2,3\}$, including solutions of the 2D Navier-Stokes equations forced by sufficiently regular-in-space, non-degenerate white-in-time noise. We prove that the solution almost surely mixes exponentially fast uniformly in the diffusivity $\kappa$. Namely, that there is a deterministic, exponential rate (independent of $\kappa$) such that all mean-zero $H^1$ initial data decays exponentially fast in $H^{-1}$ at this rate with probability one. This implies almost-sure enhanced dissipation in $L^2$. Specifically that there is a deterministic, uniform-in-$\kappa$, exponential decay in $L^2$ after time $t \gtrsim |\log \kappa|$. Both the $O(|\log \kappa|)$ time-scale and the uniform-in-$\kappa$ exponential mixing are optimal for Lipschitz velocity fields and, to our knowledge, are the first rigorous examples of velocity fields satisfying these properties (deterministic or stochastic). This work is also a major step in our program on scalar mixing and Lagrangian chaos necessary for a rigorous proof of the Batchelor power spectrum of passive scalar turbulence.


title: DeepRacer: Educational Autonomous Racing Platform for Experimentation with Sim2Real Reinforcement Learning

abstract: DeepRacer is a platform for end-to-end experimentation with RL and can be used to systematically investigate the key challenges in developing intelligent control systems. Using the platform, we demonstrate how a 1/18th scale car can learn to drive autonomously using RL with a monocular camera. It is trained in simulation with no additional tuning in physical world and demonstrates: 1) formulation and solution of a robust reinforcement learning algorithm, 2) narrowing the reality gap through joint perception and dynamics, 3) distributed on-demand compute architecture for training optimal policies, and 4) a robust evaluation method to identify when to stop training. It is the first successful large-scale deployment of deep reinforcement learning on a robotic control agent that uses only raw camera images as observations and a model-free learning method to perform robust path planning. We open source our code and video demo on GitHub: https://git.io/fjxoJ.


title: Sub-cycle coherent control of ionic dynamics via transient ionization injection

abstract: We investigate the interwoven dynamic evolutions of neutral nitrogen molecules together with nitrogen ions created through transient tunnel ionization in an intense laser field. By treating the molecules as open quantum systems, it is found that considering real-time injection of ions and strong couplings among their electronic states, nitrogen molecular ions are primarily populated in the electronically excited states, rather than staying in the ground state as predicted by the well-known tunneling theory. The unexpected result is attributed to sub-cycle switch-on of time-dependent polarization by transient ionization and dynamic Stark shift mediated near-resonant multiphoton transitions. Their combined contribution also causes that the vibrational distribution of N$_2^+$ does not comply with Franck-Condon principle. These findings corroborate the mechanism of nitrogen molecular ion lasing and are likely to be universal. The present work opens a new route to explore the important role of transient ionization injection in strong-field induced non-equilibrium dynamics.


title: Existence and Distributional Chaos of Points that are Recurrent but Not Banach Recurrent

abstract: According to the recurrent frequency, many levels of recurrent points are found, such as periodic points, almost periodic points, weakly almost periodic points, quasi-weakly almost periodic points and Banach recurrent points. In this paper, we consider symbolic dynamics and show the existence of six refined levels between Banach recurrence and general recurrence. Despite the fact that these refined levels are all null-measure under any invariant measure, we further show they carry strong topological complexity. Each refined level of recurrent points is dense in the whole space and contains an uncountable distributionally chaotic subset. For a wide range of dynamical systems, such as expansive systems with the shadowing property, we also show the distributional chaos of the points that are recurrent but not Banach recurrent.


title: Deep Collaborative Discrete Hashing with Semantic-Invariant Structure

abstract: Existing deep hashing approaches fail to fully explore semantic correlations and neglect the effect of linguistic context on visual attention learning, leading to inferior performance. This paper proposes a dual-stream learning framework, dubbed Deep Collaborative Discrete Hashing (DCDH), which constructs a discriminative common discrete space by collaboratively incorporating the shared and individual semantics deduced from visual features and semantic labels. Specifically, the context-aware representations are generated by employing the outer product of visual embeddings and semantic encodings. Moreover, we reconstruct the labels and introduce the focal loss to take advantage of frequent and rare concepts. The common binary code space is built on the joint learning of the visual representations attended by language, the semantic-invariant structure construction and the label distribution correction. Extensive experiments demonstrate the superiority of our method.


title: Action minimizing orbits in the 2-center problems with simple choreography constraint

abstract: The aim of this paper is to study the motion of $2+n$-body problem where two equal masses are assumed to be fixed. We assume that the value of each fixed mass is equal to $M>0$ and the remaining $n$ moving particles have equal masses $m>0$. According to Newton's second law and the universal gravitation law, the $n$ particles move under the interaction of each other and the affection of the two fixed particles. Also, this motion has a natural variational structure. Under the simple choreography constraint, we show that the Lagrangian action functional attains its absolute minimum on a uniform circular motion.


title: Spin chains and classical strings in two parameters $q$-deformed AdS$_3\times$S$^3$

abstract: In this paper, we study the spin chain and string excitation in the two-parameters $q$-deformed AdS$_3\times$S$^3$ proposed by Hoare. We obtain the deformed spin chain model at the fast spin limit for choices of deformed parameters. General ansatz for giant magnons are studied in great detail and complicated dispersion relation is treated perturbatively. We also study several types of hanging string solutions and their charges and spins are analyzed numerically. At last, we explore its pp-wave limit and find its solution only depends on the difference of deformed parameters.


title: Engel's law in the commodity composition of exports

abstract: Different shares of distinct commodity sectors in production, trade, and consumption illustrate how resources and capital are allocated and invested. Economic progress has been claimed to change the share distribution in a universal manner as exemplified by the Engel's law for the household expenditure and the shift from primary to manufacturing and service sector in the three sector model. Searching for large-scale quantitative evidence of such correlation, we analyze the gross-domestic product (GDP) and international trade data based on the standard international trade classification (SITC) in the period 1962 to 2000. Three categories, among ten in the SITC, are found to have their export shares significantly correlated with the GDP over countries and time; The machinery category has positive and food and crude materials have negative correlations. The export shares of commodity categories of a country are related to its GDP by a power-law with the exponents characterizing the GDP-elasticity of their export shares. The distance between two countries in terms of their export portfolios is measured to identify several clusters of countries sharing similar portfolios in 1962 and 2000. We show that the countries whose GDP is increased significantly in the period are likely to transit to the clusters displaying large share of the machinery category.

categories: q-fin.GN econ.GN physics.soc-ph q-fin.EC


title: PAPR Reduction Using Iterative Clipping/Filtering and ADMM Approaches for OFDM-Based Mixed-Numerology Systems

abstract: Mixed-numerology transmission is proposed to support a variety of communication scenarios with diverse requirements. However, as the orthogonal frequency division multiplexing (OFDM) remains as the basic waveform, the peak-to average power ratio (PAPR) problem is still cumbersome. In this paper, based on the iterative clipping and filtering (ICF) and optimization methods, we investigate the PAPR reduction in the mixed-numerology systems. We first illustrate that the direct extension of classical ICF brings about the accumulation of inter-numerology interference (INI) due to the repeated execution. By exploiting the clipping noise rather than the clipped signal, the noise-shaped ICF (NS-ICF) method is then proposed without increasing the INI. Next, we address the in-band distortion minimization problem subject to the PAPR constraint. By reformulation, the resulting model is separable in both the objective function and the constraints, and well suited for the alternating direction method of multipliers (ADMM) approach. The ADMM-based algorithms are then developed to split the original problem into several subproblems which can be easily solved with closed-form solutions. Furthermore, the applications of the proposed PAPR reduction methods combined with filtering and windowing techniques are also shown to be effective.


title: Capturing Multireference Excited States by Constrained DFT

abstract: The computation of excited electronic states with commonly employed (approximate) methods is challenging, typically yielding states of lower quality than the corresponding ground state for a higher computational cost. In this work, we present a mean field method that extends the previously proposed eXcited Constrained DFT (XCDFT) from single Slater determinants to ensemble 1-RDMs for computing low-lying excited states. The method still retains an associated computational complexity comparable to a semilocal DFT calculation while at the same time is capable of approaching states with multireference character. We benchmark the quality of this method on well-established test sets, finding good descriptions of the electronic structure of multireference states and maintaining an overall accuracy for the predicted excitation energies comparable to semilocal TDDFT.


title: A Liouville Theorem for Axi-symmetric Navier-Stokes Equations on $\mathbb{R}^2 \times \mathbb{T}^1$

abstract: We establish a Liouville theorem for bounded mild ancient solutions to the axi-symmetric incompressible Navier-Stokes equations on $(-\infty, 0] \times (\mathbb{R}^2 \times \mathbb{T}^1)$. This is a step forward to completely solve the conjecture on $(-\infty, 0] \times \mathbb{R}^3$ which was made in \cite{KNSS} to describe the potential singularity structures of the Cauchy problem.


title: Search for the field-induced magnetic instability around the upper critical field of superconductivity in H || c in CeCoIn5

abstract: We present nuclear spin-lattice relaxation rate (1/T1) at the Co site and ac-susceptibility results in the normal and superconducting (SC) states of CeCoIn5 for H || c near the SC upper critical field Hc2 above 0.1 K. At 4.2 T, 1/T1 rapidly decreases below the SC transition temperature, consistent with the previous reports. Although the field dependence of 1/T1T at 0.1 K shows a peak at 5.2 T above Hc2, the temperature dependence of 1/T1T at 5.2 T is independent of temperature below 0.2 K, showing a Fermi-liquid behavior. In addition, we found no NMR-spectrum broadening by the appearance of internal fields around Hc2 at 0.1 K. We could not detect any field-induced magnetic instability around Hc2 down to 0.1 K although the remarkable non-Fermi-liquid behavior towards Hc2 was observed in various physical quantities.


title: Improving Operational Feasibility of Low-voltage Distribution Network by Phase-switching Devices

abstract: High penetration of residential photovoltaic (PV) in low-voltage distribution networks (LVDN) makes the unbalance issue more sever due to the asymmetry of generation/load characteristic in different phases, and may lead to infeasible operation of the whole network. Phase switching device (PSD), which can switch the connected phase of a residential load as required, is viable and efficient equipment to help address this issue. This paper, based on three-phase power flow (TUPF) formulation, aims to investigate the benefit of PSD on improving the operation feasibility of LVDN. The linear model of TUPF with PSD is presented to effectively take the flexible device into account, which can be conveniently used to seek the PSD positions that lead to a viable operation strategy of the original problem when infeasibility is reported by the traditional iteration-based algorithm. Case study based on a practical LVDN in Australia demonstrates the efficacy of the proposed method.


title: A Novel, Finite-Element Model for Spin-Exchange Optical Pumping Using an Open-Source Code

abstract: A new model is presented for spin-exchange optical pumping using an open-source code, ElmerFEM-CSC. The model builds on previous models by adding the effects of alkali-vapor heterogeneity in optical pumping cells and by modeling the effects of hyperpolarized-gas wall-relaxation using a diffusion model. The code supports full, three-dimensional solutions to optical-pumping models, and solves for (1) laser absorption, (2) alkali vapor concentration, (3) fluid flow parameters, (4) thermal affects due to the pumping laser, and (5) noble gas polarization. The source code for the model is available for researchers to utilize and modify.


title: Importance Sampling via Local Sensitivity

abstract: Given a loss function $F:\mathcal{X} \rightarrow \R^+$ that can be written as the sum of losses over a large set of inputs $a_1,\ldots, a_n$, it is often desirable to approximate $F$ by subsampling the input points. Strong theoretical guarantees require taking into account the importance of each point, measured by how much its individual loss contributes to $F(x)$. Maximizing this importance over all $x \in \mathcal{X}$ yields the \emph{sensitivity score} of $a_i$. Sampling with probabilities proportional to these scores gives strong guarantees, allowing one to approximately minimize of $F$ using just the subsampled points. Unfortunately, sensitivity sampling is difficult to apply since (1) it is unclear how to efficiently compute the sensitivity scores and (2) the sample size required is often impractically large. To overcome both obstacles we introduce \emph{local sensitivity}, which measures data point importance in a ball around some center $x_0$. We show that the local sensitivity can be efficiently estimated using the \emph{leverage scores} of a quadratic approximation to $F$ and that the sample size required to approximate $F$ around $x_0$ can be bounded. We propose employing local sensitivity sampling in an iterative optimization method and analyze its convergence when $F$ is smooth and convex.


title: Correcting for attenuation due to measurement error

abstract: I present a frequentist method for quantifying uncertainty when correcting correlations for attenuation due to measurement error. The method is conservative but has far better coverage properties than the methods currently used when sample sizes are small. I recommend the use of confidence curves in favor of confidence intervals when this method is used. I introduce the R package "attenuation" which can be used to calculate and visualize the methods described in this paper.


title: Improving Long Handwritten Text Line Recognition with Convolutional Multi-way Associative Memory

abstract: Convolutional Recurrent Neural Networks (CRNNs) excel at scene text recognition. Unfortunately, they are likely to suffer from vanishing/exploding gradient problems when processing long text images, which are commonly found in scanned documents. This poses a major challenge to goal of completely solving Optical Character Recognition (OCR) problem. Inspired by recently proposed memory-augmented neural networks (MANNs) for long-term sequential modeling, we present a new architecture dubbed Convolutional Multi-way Associative Memory (CMAM) to tackle the limitation of current CRNNs. By leveraging recent memory accessing mechanisms in MANNs, our architecture demonstrates superior performance against other CRNN counterparts in three real-world long text OCR datasets.


title: Singular principal bundles on reducible nodal curves

abstract: Studying degenerations of moduli spaces of semistable principal bundles on smooth curves leads to the problem of constructing and studying moduli spaces on singular curves. In this note, we will see that the moduli spaces of $\delta$-semistable pseudo bundles on a nodal curve constructed by the first author become, for large values of $\delta$, the moduli spaces for semistable singular principal bundles. The latter are reasonable candidates for degenerations and a potential basis of further developments as on irreducible nodal curves. In particular, we find a notion of semistability for principal bundles on reducible nodal curves. The understanding of the asymptotic behavior of $\delta$-semistability rests on a lemma from geometric invariant theory. The results will allow the construction of a universal moduli space of semistable singular principal bundles relative to the moduli space $\overline{\mathcal M}_g$ of stable curves of genus $g$.


title: A note on the lower bound of representation functions

abstract: For a set $A$ of nonnegative integers, let $R_2(A,n)$ denote the number of solutions to $n=a+a'$ with $a,a'\in A$, $a<a'$. Let $A_0$ be the Thue-Morse sequence and $B_0=\mathbb{N}\setminus A_0$. Let $A\subset \mathbb{N}$ and $N$ be a positive integer such that $R_2(A,n)=R_2(\mathbb{N}\setminus A,n)$ for all $n\geq 2N-1$. Previously, the first author proved that if $|A\cap A_0|=+\infty$ and $|A\cap B_0|=+\infty$, then $R_2(A,n)\geq \frac{n+3}{56N-52}-1$ for all $n\geq 1$. In this paper, we prove that the above lower bound is nearly best possible. We also get some other results.


title: Polarized primordial gravitational waves in the ghost-free parity-violating gravity

abstract: The tests of parity symmetry in the gravitational interaction is an attractive issue in the gravitational-wave astronomy. In the general theories of gravity with parity violation, one of the fundamental results is that the primordial gravitational waves (PGWs) produced during the slow-roll inflation is circularly polarized. In this article, we investigate the polarization of PGWs in the recently proposed ghost-free parity-violating gravity, which generalizes the Chern-Simons gravity by including the higher derivatives of the coupling scalar field. For this purpose, we first construct the approximate analytical solution to the mode function of the PGWs during the slow-roll inflation by using the uniform asymptotic approximation. With the approximate solution, we calculate explicitly the power spectrum and the corresponding circular polarization of the PGWs analytically, and find that the contributions of the higher derivatives of the coupling scalar field to the circular polarization is of the same order of magnitude as that of Chern-Simons gravity. The degree of circular polarization of PGWs is suppressed by the energy scale of parity violation in gravity, which is hardly detected on the basis solely of two-point statistics from future cosmic microwave background data.