title: Minimal generating set of Sylow 2-subgroups commutator subgroup of alternating group. Commutator width in Sylow $p$-subgroups of alternating, symmetric groups and in the wreath product of groups

abstract: The size of minimal generating set for commutator of Sylow 2-subgroup of alternating group was found. Given a permutational wreath product of finite cyclic groups sequence we prove that the commutator width of such groups is 1 and we research some properties of its commutator subgroup. It was shown that $(Syl_2 A_{2^k})^2 = Syl'_2 (A_{2^k}), \, k>2$. A new approach to presentation of Sylow 2-subgroups of alternating group ${A_{{2^{k}}}}$ was applied. As a result the short proof that the commutator width of Sylow 2-subgroups of alternating group ${A_{{2^{k}}}}$, permutation group ${S_{{2^{k}}}}$ and Sylow $p$-subgroups of $Syl_2 A_{p^k}$ ($Syl_2 S_{p^k}$) are equal to 1 was obtained. Commutator width of permutational wreath product $B \wr C_n$ were investigated. It was proven that the commutator length of an arbitrary element of commutator of the wreath product of cyclic groups $C_{p_i}, \, p_i\in \mathbb{N} $ equals to 1. The commutator width of direct limit of wreath product of cyclic groups are found. As a corollary, it was shown that the commutator width of Sylows $p$-subgroups $Syl_2(S_{{p^{k}}})$ of symmetric $S_{{p^{k}}}$ and alternating groups $A_{{p^{k}}}$ $p \geq 2$ are also equal to 1. A recursive presentation of Sylows $2$-subgroups $Syl_2(A_{{2^{k}}})$ of $A_{{2^{k}}}$ was introduced. The structure of Sylows $2$-subgroups commutator of symmetric and alternating groups were investigated. For an arbitrary group $B$ an upper bound of commutator width of $C_p \wr B$ was founded.


title: A Method for Measuring Network Effects of One-to-One Communication Features in Online A/B Tests

abstract: A/B testing is an important decision making tool in product development because can provide an accurate estimate of the average treatment effect of a new features, which allows developers to understand how the business impact of new changes to products or algorithms. However, an important assumption of A/B testing, Stable Unit Treatment Value Assumption (SUTVA), is not always a valid assumption to make, especially for products that facilitate interactions between individuals. In contexts like one-to-one messaging we should expect network interference; if an experimental manipulation is effective, behavior of the treatment group is likely to influence members in the control group by sending them messages, violating this assumption. In this paper, we propose a novel method that can be used to account for network effects when A/B testing changes to one-to-one interactions. Our method is an edge-based analysis that can be applied to standard Bernoulli randomized experiments to retrieve an average treatment effect that is not influenced by network interference. We develop a theoretical model, and methods for computing point estimates and variances of effects of interest via network-consistent permutation testing. We then apply our technique to real data from experiments conducted on the messaging product at LinkedIn. We find empirical support for our model, and evidence that the standard method of analysis for A/B tests underestimates the impact of new features in one-to-one messaging contexts.


title: Search for production of an invisible dark photon in $\pi^0$ decays

abstract: The results of a search for $\pi^0$ decays to a photon and an invisible massive dark photon at the NA62 experiment at the CERN SPS are reported. From a total of $4.12\times10^8$ tagged $\pi^0$ mesons, no signal is observed. Assuming a kinetic-mixing interaction, limits are set on the dark photon coupling to the ordinary photon as a function of the dark photon mass, improving on previous searches in the mass range 60--110 MeV/$c^2$. The present results are interpreted in terms of an upper limit of the branching ratio of the electro-weak decay $\pi^0 \to \gamma \nu \overline{\nu}$, improving the current limit by more than three orders of magnitude.


title: Anomaly Flow and T-Duality

abstract: In this paper, we study the dual Anomaly flow, which is a dual version of the Anomaly flow under T-duality. A family of monotone functionals is introduced and used to estimate the dilaton function along the flow. Many examples and reductions of the dual Anomaly flow are worked out in detail.


title: An analytical model of radial dust trapping in protoplanetary disks

abstract: We study dust concentration in axisymmetric gas rings in protoplanetary disks. Given the gas surface density, we derived an analytical total dust surface density by taking into account the differential concentration of all the grain sizes. This model allows us to predict the local dust-to-gas mass ratio and the slope of the particle size distribution, as a function of radius. We test this analytical model comparing it with a 3D magneto-hydrodynamical simulation of dust evolution in an accretion disk. The model is also applied to the disk around HD 169142. By fitting the disk continuum observations simultaneously at $\lambda = 0.87$, 1.3, 3.0 mm, we obtain a global dust-to-gas mass ratio $\epsilon_{\rm global} = 1.05 \times 10^{-2}$ and a viscosity coefficient $\alpha = 1.35 \times 10^{-2}$. This model can be easily implemented in numerical simulations of accretion disks.


title: Syzygies in Hilbert schemes of complete intersections

abstract: Let $ e_1, ..., e_c $ be positive integers and let $ Y \subseteq \mathbb{P}^n$ be the monomial complete intersection defined by the vanishing of $x_1^{e_1}, ..., x_c^{e_c}$. In this paper we study sharp upper bounds on the number of equations and syzygies of subschemes parametrized by the Hilbert scheme of points $Hilb^d(Y)$, and discuss applications to the Hilbert scheme of points $Hilb^d(X)$ of arbitrary complete intersections $X \subseteq \mathbb{P}^n$.


title: Monotonicity-based inversion of the fractional Schr\"odinger equation II. General potentials and stability

abstract: In this work, we use monotonicity-based methods for the fractional Schr\"odinger equation with general potentials $q\in L^\infty(\Omega)$ in a Lipschitz bounded open set $\Omega\subset \mathbb R^n$ in any dimension $n\in \mathbb N$. We demonstrate that if-and-only-if monotonicity relations between potentials and the Dirichlet-to-Neumann map hold up to a finite dimensional subspace. Based on these if-and-only-if monotonicity relations, we derive a constructive global uniqueness results for the fractional Calder\'on problem and its linearized version. We also derive a reconstruction method for unknown obstacles in a given domain that only requires the background solution of the fractional Schr\"odinger equation, and we prove uniqueness and Lipschitz stability from finitely many measurements for potentials lying in an a-priori known bounded set in a finite dimensional subset of $L^\infty(\Omega)$.


title: ToyArchitecture: Unsupervised Learning of Interpretable Models of the World

abstract: Research in Artificial Intelligence (AI) has focused mostly on two extremes: either on small improvements in narrow AI domains, or on universal theoretical frameworks which are usually uncomputable, incompatible with theories of biological intelligence, or lack practical implementations. The goal of this work is to combine the main advantages of the two: to follow a big picture view, while providing a particular theory and its implementation. In contrast with purely theoretical approaches, the resulting architecture should be usable in realistic settings, but also form the core of a framework containing all the basic mechanisms, into which it should be easier to integrate additional required functionality. In this paper, we present a novel, purposely simple, and interpretable hierarchical architecture which combines multiple different mechanisms into one system: unsupervised learning of a model of the world, learning the influence of one's own actions on the world, model-based reinforcement learning, hierarchical planning and plan execution, and symbolic/sub-symbolic integration in general. The learned model is stored in the form of hierarchical representations with the following properties: 1) they are increasingly more abstract, but can retain details when needed, and 2) they are easy to manipulate in their local and symbolic-like form, thus also allowing one to observe the learning process at each level of abstraction. On all levels of the system, the representation of the data can be interpreted in both a symbolic and a sub-symbolic manner. This enables the architecture to learn efficiently using sub-symbolic methods and to employ symbolic inference.


title: Robust Image Segmentation Quality Assessment

abstract: Deep learning based image segmentation methods have achieved great success, even having human-level accuracy in some applications. However, due to the black box nature of deep learning, the best method may fail in some situations. Thus predicting segmentation quality without ground truth would be very crucial especially in clinical practice. Recently, people proposed to train neural networks to estimate the quality score by regression. Although it can achieve promising prediction accuracy, the network suffers robustness problem, e.g. it is vulnerable to adversarial attacks. In this paper, we propose to alleviate this problem by utilizing the difference between the input image and the reconstructed image, which is conditioned on the segmentation to be assessed, to lower the chance to overfit to the undesired image features from the original input image, and thus to increase the robustness. Results on ACDC17 dataset demonstrated our method is promising.


title: Permutation patterns in genome rearrangement problems: the reversal model

abstract: In the context of the genome rearrangement problem, we analyze two well known models, namely the reversal and the prefix reversal models, by exploiting the connection with the notion of permutation pattern. More specifically, for any $k$, we provide a characterization of the set of permutations having distance $\leq k$ from the identity (which is known to be a permutation class) in terms of what we call generating peg permutations and we describe some properties of its basis, which allow to compute such a basis for small values of $k$.


title: Hidden Asymptotic Symmetry in Long Elastic Beams on Softening Foundations

abstract: Transverse wrinkles are known to appear in thin rectangular elastic sheets when stretched in the long direction. Numerically computed bifurcation diagrams for extremely thin, highly stretched films indicate entire orbits of wrinkling solutions, cf. Healey, et. al. [J. Nonlinear Sci., 23 (2013), pp.~777--805]. These correspond to arbitrary phase shifts of the wrinkled pattern in the transverse direction. While such behavior is normally associated with problems in the presence of a continuous symmetry group, an unloaded rectangular sheet possesses only a finite symmetry group. In order to understand this phenomenon, we consider a simpler problem more amenable to analysis -- a finite-length beam on a nonlinear softening foundation under axial compression. We first obtain asymptotic results via amplitude equations, that are valid as a certain non-dimensional beam length becomes sufficiently large. We deduce that any two phase-shifts of a solution differ from one another by exponentially small terms in that length. We validate this observation with numerical computations, indicating the presence of solution orbits for sufficiently long beams. We refer to this as "hidden asymptotic symmetry".


title: Linear quadratic mean field games: Asymptotic solvability and relation to the fixed point approach

abstract: Mean field game theory has been developed largely following two routes. One of them, called the direct approach, starts by solving a large-scale game and next derives a set of limiting equations as the population size tends to infinity. The second route is to apply mean field approximations and formalize a fixed point problem by analyzing the best response of a representative player. This paper addresses the connection and difference of the two approaches in a linear quadratic (LQ) setting. We first introduce an asymptotic solvability notion for the direct approach, which means for all sufficiently large population sizes, the corresponding game has a set of feedback Nash strategies in addition to a mild regularity requirement. We provide a necessary and sufficient condition for asymptotic solvability and show that in this case the solution converges to a mean field limit. This is accomplished by developing a re-scaling method to derive a low dimensional ordinary differential equation (ODE) system, where a non-symmetric Riccati ODE has a central role. We next compare with the fixed point approach which determines a two point boundary value (TPBV) problem, and show that asymptotic solvability implies feasibility of the fixed point approach, but the converse is not true. We further address non-uniqueness in the fixed point approach and examine the long time behavior of the non-symmetric Riccati ODE in the asymptotic solvability problem.


title: The Disk Gas Mass and the Far-IR Revolution

abstract: The gaseous mass of protoplanetary disks is a fundamental quantity in planet formation. The presence of gas is necessary to assemble planetesimals, it determines timescales of giant planet birth, and it is an unknown factor for a wide range of properties of planet formation, from chemical abundances (X/H) to the mass efficiency of planet formation. The gas mass obtained from traditional tracers, such as dust thermal continuum and CO isotopologues, are now known to have significant (1 - 2 orders of magnitude) discrepancies. Emission from the isotopologue of H2, hydrogen deuteride (HD), offers an alternative measurement of the disk gas mass. Of all of the regions of the spectrum, the far-infrared stands out in that orders of magnitude gains in sensitivity can be gleaned by cooling a large aperture telescope to 8 K. Such a facility can open up a vast new area of the spectrum to exploration. One of the primary benefits of this far-infrared revolution would be the ability to survey hundreds of planet-forming disks in HD emission to derive their gaseous masses. For the first time, we will have statistics on the gas mass as a function of evolution, tracing birth to dispersal as a function of stellar spectral type. These measurements have broad implications for our understanding of the time scale during which gas is available to form giant planets, the dynamical evolution of the seeds of terrestrial worlds, and the resulting chemical composition of pre-planetary embryos carrying the elements needed for life. Measurements of the ground-state line of HD requires a space-based observatory operating in the far-infrared at 112 microns.


title: Provable Certificates for Adversarial Examples: Fitting a Ball in the Union of Polytopes

abstract: We propose a novel method for computing exact pointwise robustness of deep neural networks for all convex $\ell_p$ norms. Our algorithm, GeoCert, finds the largest $\ell_p$ ball centered at an input point $x_0$, within which the output class of a given neural network with ReLU nonlinearities remains unchanged. We relate the problem of computing pointwise robustness of these networks to that of computing the maximum norm ball with a fixed center that can be contained in a non-convex polytope. This is a challenging problem in general, however we show that there exists an efficient algorithm to compute this for polyhedral complices. Further we show that piecewise linear neural networks partition the input space into a polyhedral complex. Our algorithm has the ability to almost immediately output a nontrivial lower bound to the pointwise robustness which is iteratively improved until it ultimately becomes tight. We empirically show that our approach generates distance lower bounds that are tighter compared to prior work, under moderate time constraints.


title: Effective upper bound of analytic torsion under Arakelov metric

abstract: Given a choice of metric on the Riemann surface, the regularized determinant of Laplacian (analytic torsion) is defined via the complex power of elliptic operators: $$ \det(\Delta)=\exp(-\zeta'(0)) $$ In this paper we gave an asymptotic effective estimate of analytic torsion under Arakelov metric. In particular, after taking the logarithm it is asymptotically upper bounded by $g$ for $g>1$. The construction of a cohomology theory for arithmetic surfaces in Arakelov theory has long been an open problem. In particular, it is not known if $h^{1}(X,L)\ge 0$. We view this as an indirect piece of evidence that if such a cohomology theory exists, the $h^{1}$ term may be effectively estimated.


title: Linear quadratic mean field games with a major player: The multi-scale approach

abstract: This paper considers linear quadratic (LQ) mean field games with a major player and analyzes an asymptotic solvability problem. It starts with a large-scale system of coupled dynamic programming equations and applies a re-scaling technique introduced in Huang and Zhou (2018a, 2018b) to derive a set of Riccati equations in lower dimensions, the solvability of which determines the necessary and sufficient condition for asymptotic solvability. We next derive the mean field limit of the strategies and the value functions. Finally, we show that the two decentralized strategies can be interpreted as the best responses of a major player and a representative minor player embedded in an infinite population, which have the property of consistent mean field approximations.


title: Fault-Tolerant Nanosatellite Computing on a Budget

abstract: Micro- and nanosatellites have become popular platforms for a variety of commercial and scientific applications, but today are considered suitable mainly for short and low-priority space missions due to their low reliability. In part, this can be attributed to their reliance upon cheap, low-feature size, COTS components originally designed for embedded and mobile-market applications, for which traditional hardware-voting concepts are ineffective. Software-fault-tolerance concepts have been shown effective for such systems, but have largely been ignored by the space industry due to low maturity, as most have only been researched in theory. In practice, designers of payload instruments and miniaturized satellites are usually forced to sacrifice reliability in favor deliver the level of performance necessary for cutting-edge science and innovative commercial applications. Thus, we developed a software-fault-tolerance-approach based upon thread-level coarse-grain lockstep, which was validated using fault-injection. To offer strong long-term fault coverage, our architecture is implemented as tiled MPSoC on an FPGA, utilizing partial reconfiguration, as well as mixed criticality. This architecture can satisfy the high performance requirements of current and future scientific and commercial space missions at very low cost, while offering the strong fault-coverage guarantees necessary for platform control even for missions with a long duration. This architecture was developed for a 4-year ESA project. Together with two industrial partners, we are developing a prototype to then undergo radiation testing.


title: Epstein-Zin Utility Maximization on a Random Horizon

abstract: This paper solves the consumption-investment problem under Epstein-Zin preferences on a random horizon. In an incomplete market, we take the random horizon to be a stopping time adapted to the market filtration, generated by all observable, but not necessarily tradable, state processes. Contrary to prior studies, we do not impose any fixed upper bound for the random horizon, allowing for truly unbounded ones. Focusing on the empirically relevant case where the risk aversion and the elasticity of intertemporal substitution are both larger than one, we characterize the optimal consumption and investment strategies using backward stochastic differential equations with superlinear growth on unbounded random horizons. This characterization, compared with the classical fixed-horizon result, involves an additional stochastic process that serves to capture the randomness of the horizon. As demonstrated in two concrete examples, changing from a fixed horizon to a random one drastically alters the optimal strategies.


title: A Reciprocal Formulation of Nonexponential Radiative Transfer. 3: Binary Mixtures

abstract: We derive the form of reciprocal generalized radiative transfer (RGRT) that includes the Levermore-Pomraning attenuation law for paths leaving a deterministic origin. The resulting model describes linear transport within multi-dimensional stochastic binary mixtures with Markovian mixing statistics and nonstochastic albedo and phase function. The derivation includes a new attenuation law and related free-path distribution between collisions, which was previously estimated using a Monte Carlo approach. We also show that previous stationary descriptions of binary mixtures with binomial mixing statistics reduce to exponential attenuations and, thus, have an exact homogenization within classical radiative transfer.


title: Real-time Ecological Velocity Planning for Plug-in Hybrid Vehicles with Partial Communication to Traffic Lights

abstract: This paper presents the design of an ecological adaptive cruise controller (ECO-ACC) for a plug-in hybrid vehicle (PHEV) which exploits automated driving and connectivity. Most existing papers for ECO-ACC focus on a short-sighted control scheme. A two-level control framework for long-sighted ECO-ACC was only recently introduced. However, that work is based on a deterministic traffic signal phase and timing (SPaT) over the entire route. In practice, connectivity with traffic lights may be limited by communication range, e.g. just one upcoming traffic light. We propose a two-level receding-horizon control framework for long-sighted ECO-ACC that exploits deterministic SPaT for the upcoming traffic light, and utilizes historical SPaT for other traffic lights within a receding control horizon. We also incorporate a powertrain control mechanism to enhance PHEV energy prediction accuracy. Hardware-in-the-loop simulation results validate the energy savings of the receding-horizon control framework in various traffic scenarios.


title: Efficient demultiplexed single-photon source with a quantum dot coupled to a nanophotonic waveguide

abstract: Planar nanostructures allow near-ideal extraction of emission from a quantum emitter embedded within, thereby realizing deterministic single-photon sources. Such a source can be transformed into M single-photon sources by implementing active temporal-to-spatial mode demultiplexing. We report on the realization of such a demultiplexed source based on a quantum dot embedded in a nanophotonic waveguide. Efficient outcoupling (>60%) from the waveguide into a single mode optical fiber is obtained with high-efficiency grating couplers. As a proof-of-concept, active demultiplexing into M=4 spatial channels is demonstrated by the use of electro-optic modulators with an end-to-end efficiency of >81% into single-mode fibers. Overall we demonstrate four-photon coincidence rates of >1 Hz even under non-resonant excitation of the quantum dot. The main limitation of the current source is the residual population of other exciton transitions that corresponds to a finite preparation efficiency of the desired transition. We quantitatively extract a preparation efficiency of 15% using the second-order correlation function measurements. The experiment highlights the applicability of planar nanostructures as efficient multiphoton sources through temporal-to-spatial demultiplexing and lays out a clear path way of how to scale up towards demonstrating quantum advantages with the quantum dot sources.


title: Positive solutions for non-variational fractional elliptic systems with negative exponents

abstract: In this paper, we study strongly coupled elliptic systems in non-variational form with negative exponents involving fractional Laplace operators. We investigate the existence, nonexistence, and uniqueness of the positive classical solution. The results obtained here are a natural extension of the results obtained by Ghergu (2010), for the fractional case.


title: Quot schemes of curves and surfaces: virtual classes, integrals, Euler characteristics

abstract: We compute tautological integrals over Quot schemes on curves and surfaces. After obtaining several explicit formulas over Quot schemes of dimension 0 quotients on curves (and finding a new symmetry), we apply the results to tautological integrals against the virtual fundamental classes of Quot schemes of dimension 0 and 1 quotients on surfaces (using also universality, torus localization, and cosection localization). The virtual Euler characteristics of Quot schemes of surfaces, a new theory parallel to the Vafa-Witten Euler characteristics of the moduli of bundles, is defined and studied. Complete formulas for the virtual Euler characteristics are found in the case of dimension 0 quotients on surfaces. Dimension 1 quotients are studied on K3 surfaces and surfaces of general type with connections to the Kawai-Yoshioka formula and the Seiberg-Witten invariants respectively. The dimension 1 theory is completely solved for minimal surfaces of general type admitting a nonsingular canonical curve. Along the way, we find a new connection between weighted tree counting and multivariate Fuss-Catalan numbers which is of independent interest.


title: Selective Attention for Context-aware Neural Machine Translation

abstract: Despite the progress made in sentence-level NMT, current systems still fall short at achieving fluent, good quality translation for a full document. Recent works in context-aware NMT consider only a few previous sentences as context and may not scale to entire documents. To this end, we propose a novel and scalable top-down approach to hierarchical attention for context-aware NMT which uses sparse attention to selectively focus on relevant sentences in the document context and then attends to key words in those sentences. We also propose single-level attention approaches based on sentence or word-level information in the context. The document-level context representation, produced from these attention modules, is integrated into the encoder or decoder of the Transformer model depending on whether we use monolingual or bilingual context. Our experiments and evaluation on English-German datasets in different document MT settings show that our selective attention approach not only significantly outperforms context-agnostic baselines but also surpasses context-aware baselines in most cases.


title: Interpreting Neural Networks Using Flip Points

abstract: Neural networks have been criticized for their lack of easy interpretation, which undermines confidence in their use for important applications. Here, we introduce a novel technique, interpreting a trained neural network by investigating its flip points. A flip point is any point that lies on the boundary between two output classes: e.g. for a neural network with a binary yes/no output, a flip point is any input that generates equal scores for "yes" and "no". The flip point closest to a given input is of particular importance, and this point is the solution to a well-posed optimization problem. This paper gives an overview of the uses of flip points and how they are computed. Through results on standard datasets, we demonstrate how flip points can be used to provide detailed interpretation of the output produced by a neural network. Moreover, for a given input, flip points enable us to measure confidence in the correctness of outputs much more effectively than softmax score. They also identify influential features of the inputs, identify bias, and find changes in the input that change the output of the model. We show that distance between an input and the closest flip point identifies the most influential points in the training data. Using principal component analysis (PCA) and rank-revealing QR factorization (RR-QR), the set of directions from each training input to its closest flip point provides explanations of how a trained neural network processes an entire dataset: what features are most important for classification into a given class, which features are most responsible for particular misclassifications, how an adversary might fool the network, etc. Although we investigate flip points for neural networks, their usefulness is actually model-agnostic.


title: Constraining the Rastall parameters in static spacetimes with galaxy-scale strong gravitational lensing

abstract: Recently, Rastall gravity is undergoing a significant surge in popularity. We obtain a power-law total mass-density profile for the inner region (within several effective radius) of early-type galaxies (ETGs) from the space-time structures which are described by the static spherically-symmetric solutions of Rastall gravity under the assumption of perfect fluid matter. We find that in the inner region of ETGs, the Rastall dimensionless parameter $\beta=\kappa\lambda$ determines the mass distribution. We then use 118 galaxy-galaxy strong gravitational lensing systems to constrain the Rastall dimensionless parameter $\beta$. We find that the mean value of $\beta$ for total 118 ETGs is $\beta=0.163\pm0.001$(68\% CL) with a minor intrinsic scatter of $\delta=0.020\pm 0.001$. Our work observationally illustrates the physical meaning of the Rastall dimensionless parameter in galaxy scale. From the Newtonian approximation of Rastall gravity, we also find that an absolute isothermal mass distribution for ETGs is not allowed in the framework of Rastall gravity.


title: On generalized near-group fusion categories

abstract: In this paper, we study the structure of a generalized near-group fusion category and classified it when it is slightly degenerate.


title: End-to-End Safe Reinforcement Learning through Barrier Functions for Safety-Critical Continuous Control Tasks

abstract: Reinforcement Learning (RL) algorithms have found limited success beyond simulated applications, and one main reason is the absence of safety guarantees during the learning process. Real world systems would realistically fail or break before an optimal controller can be learned. To address this issue, we propose a controller architecture that combines (1) a model-free RL-based controller with (2) model-based controllers utilizing control barrier functions (CBFs) and (3) on-line learning of the unknown system dynamics, in order to ensure safety during learning. Our general framework leverages the success of RL algorithms to learn high-performance controllers, while the CBF-based controllers both guarantee safety and guide the learning process by constraining the set of explorable polices. We utilize Gaussian Processes (GPs) to model the system dynamics and its uncertainties. Our novel controller synthesis algorithm, RL-CBF, guarantees safety with high probability during the learning process, regardless of the RL algorithm used, and demonstrates greater policy exploration efficiency. We test our algorithm on (1) control of an inverted pendulum and (2) autonomous car-following with wireless vehicle-to-vehicle communication, and show that our algorithm attains much greater sample efficiency in learning than other state-of-the-art algorithms and maintains safety during the entire learning process.


title: Slightly trivial extensions of a fusion category

abstract: We introduce and study the notion of slightly trivial extensions of a fusion category which can be viewed as the first level of complexity of extensions. We also provide two examples of slightly trivial extensions which arise from rank $3$ fusion categories.


title: Parallel Batch-Dynamic Graph Connectivity

abstract: In this paper, we study batch parallel algorithms for the dynamic connectivity problem, a fundamental problem that has received considerable attention in the sequential setting. The most well known sequential algorithm for dynamic connectivity is the elegant level-set algorithm of Holm, de Lichtenberg and Thorup (HDT), which achieves $O(\log^2 n)$ amortized time per edge insertion or deletion, and $O(\log n / \log\log n)$ time per query. We design a parallel batch-dynamic connectivity algorithm that is work-efficient with respect to the HDT algorithm for small batch sizes, and is asymptotically faster when the average batch size is sufficiently large. Given a sequence of batched updates, where $\Delta$ is the average batch size of all deletions, our algorithm achieves $O(\log n \log(1 + n / \Delta))$ expected amortized work per edge insertion and deletion and $O(\log^3 n)$ depth w.h.p. Our algorithm answers a batch of $k$ connectivity queries in $O(k \log(1 + n/k))$ expected work and $O(\log n)$ depth w.h.p. To the best of our knowledge, our algorithm is the first parallel batch-dynamic algorithm for connectivity.


title: Largest 2-regular subgraphs in 3-regular graphs

abstract: For a graph $G$, let $f_2(G)$ denote the largest number of vertices in a $2$-regular subgraph of $G$. We determine the minimum of $f_2(G)$ over $3$-regular $n$-vertex simple graphs $G$. To do this, we prove that every $3$-regular multigraph with exactly $c$ cut-edges has a $2$-regular subgraph that omits at most $\max\{0,\lfloor (c-1)/2\rfloor\}$ vertices. More generally, every $n$-vertex multigraph with maximum degree $3$ and $m$ edges has a $2$-regular subgraph that omits at most $\max\{0,\lfloor (3n-2m+c-1)/2\rfloor\}$ vertices. These bounds are sharp; we describe the extremal multigraphs.


title: Evidence for a Young Stellar Population in Nearby Type 1 Active Galaxies

abstract: To understand the physical origin of the close connection between supermassive black holes and their host galaxies, it is vital to investigate star formation properties in active galaxies. Using a large dataset of nearby type 1 active galactic nuclei (AGNs) with detailed structural decomposition based on high-resolution optical images obtained with the Hubble Space Telescope, we study the correlation between black hole mass and bulge luminosity and the (Kormendy) relation between bulge effective radius and surface brightness. In both relations, the bulges of type 1 AGNs tend to be more luminous than those of inactive galaxies with the same black hole mass or the same bulge size. This suggests that the central regions of AGN host galaxies have characteristically lower mass-to-light ratios than inactive galaxies, most likely due to the presence of a younger stellar population in active systems. In addition, the degree of luminosity excess appears to be proportional to the accretion rate of the AGN, revealing a physical connection between stellar growth and black hole growth. Adopting a simple toy model for the increase of stellar mass and black hole mass, we show that the fraction of young stellar population flattens out toward high accretion rates, possibly reflecting the influence of AGN-driven feedback.


title: Magnetic Skyrmions in Atomic Thin CrI$_3$ Monolayer

abstract: In this letter, we report the visualization of topologically protected spin textures, in the form of magnetic skyrmions, in recently discovered monoatomic-thin two-dimensional CrI$_3$. By combining density functional theory and atomistic spin dynamic simulation, we demonstrate that an application of out-of-plane electric field to CrI$_3$ lattice favors the formation of sub-10 nm skyrmions at 0 K temperature. The spin texture arises due to a strong correlation between magneto-crystalline anisotropy, Dzyaloshinskii-Moriya interaction and the vertical electric field, whose shape and size could be tuned with the magnetic field. Such finding will open new avenues for atomic-scale quantum engineering and precision sensing.


title: Emulators for the non-linear matter power spectrum beyond $\Lambda$CDM

abstract: Accurate predictions for the non-linear matter power spectrum are needed to confront theory with observations in current and near future weak lensing and galaxy clustering surveys. We propose a computationally cheap method to create an emulator for modified gravity models by utilizing existing emulators for $\Lambda$CDM. Using a suite of $N$-body simulations we construct a fitting function for the enhancement of both the linear and non-linear matter power spectrum in the commonly studied Hu-Sawicki $f(R)$ gravity model valid for wave-numbers $k \lesssim 5-10\, h\text{Mpc}^{-1}$ and redshifts $z \lesssim 3$. We show that the cosmology dependence of this enhancement is relatively weak so that our fit, using simulations coming from only one cosmology, can be used to get accurate predictions for other cosmological parameters. We also show that the cosmology dependence can, if needed, be included by using linear theory, approximate $N$-body simulations (such as COLA) and semi-analytical tools like the halo model. Our final fit can easily be combined with any emulator or semi-analytical models for the non-linear $\Lambda$CDM power spectrum to accurately, and quickly, produce a non-linear power spectrum for this particular modified gravity model. The method we use can be applied to fairly cheaply construct an emulator for other modified gravity models. As an application of our fitting formula we use it to compute Fisher-forecasts for how well galaxy clustering and weak lensing in a Euclid-like survey will be at constraining modifications of gravity.


title: The pure cohomology of multiplicative quiver varieties

abstract: To a quiver $Q$ and choices of nonzero scalars $q_i$, non-negative integers $\alpha_i$, and integers $\theta_i$ labeling each vertex $i$, Crawley-Boevey--Shaw associate a "multiplicative quiver variety" $\mathcal{M}_\theta^q(\alpha)$, a trigonometric analogue of the Nakajima quiver variety associated to $Q$, $\alpha$, and $\theta$. We prove that the pure cohomology, in the Hodge-theoretic sense, of the stable locus $\mathcal{M}_\theta^q(\alpha)^s$ is generated as a $\mathbb{Q}$-algebra by the tautological characteristic classes. In particular, the pure cohomology of genus $g$ twisted character varieties of $GL_n$ is generated by tautological classes.


title: Impact of Temperature-Dependent Rattling Phonons on Lattice Dynamics and Thermal Transport in Ag$_{6}$Ge$_{10}$P$_{12}$

abstract: Crystalline compounds exhibiting low-frequency rattling phonons constitute an important class of high-performance thermoelectrics owning to their intrinsically very low lattice thermal conductivity ($\kappa_{l}$). Theoretical approach that is capable of revealing the physical origin and accurately predicting $\kappa_{l}$ is of particular interest, which, however, still remains an outstanding challenge. In this study, we perform a case study of lattice dynamics and thermal transport properties of Ag$_{6}$Ge$_{10}$P$_{12}$, which has recently been identified as a high-performance thermoelectric phosphide due to low $\kappa_{l}$, arising from rattling vibrations associated with Ag$_{6}$ clusters. Analysis within a first-principles-based lattice-dynamics framework based on self-consistent phonon theory reveals a strong temperature dependence of rattling phonons due to high-order anharmonic interactions. Anharmonic hardening of the rattling optical modes has a strong effect on the lifetimes of heat-carrying acoustic phonons by decreasing the rate of three-phonon combination processes. This mechanism results in a significant increase in $\kappa_l$ and changes its temperature dependence to $1/T^{0.64}$.


title: A Unified Analytical Framework for Trustable Machine Learning and Automation Running with Blockchain

abstract: Traditional machine learning algorithms use data from databases that are mutable, and therefore the data cannot be fully trusted. Also, the machine learning process is difficult to automate. This paper proposes building a trustable machine learning system by using blockchain technology, which can store data in a permanent and immutable way. In addition, smart contracts are used to automate the machine learning process. This paper makes three contributions. First, it establishes a link between machine learning technology and blockchain technology. Previously, machine learning and blockchain have been considered two independent technologies without an obvious link. Second, it proposes a unified analytical framework for trustable machine learning by using blockchain technology. This unified framework solves both the trustability and automation issues in machine learning. Third, it enables a computer to translate core machine learning implementation from a single thread on a single machine to multiple threads on multiple machines running with blockchain by using a unified approach. The paper uses association rule mining as an example to demonstrate how trustable machine learning can be implemented with blockchain, and it shows how this approach can be used to analyze opioid prescriptions to help combat the opioid crisis.


title: Calculations of Bethe logarithms for Hydrogen and Helium atoms using B-splines in different gauges

abstract: The efficient and simple B-splines variational method of calculating the hydrogen atom Bethe logarithms in the acceleration gauge [Y-B Tang et.al., Phys. Rev. A $\mathbf{87}$, 022510 (2013)] is successfully applied to other gauges. The ground state Bethe logarithm of H with fourteen accurate figures is obtained in the velocity gauge, and in the length gauge the ground state value has eleven accurate figures. Present velocity- and length-gauge results for the $ns$, $np$, $nd$, and $nf$ states up to $n = 200$ of H are at the 10$^{-10}$ level of accuracy. This represents the successful variational attempt to calculate Bethe logarithms of the hydrogen atom in the velocity and length gauges. Furthermore, the B-splines variational method is successfully extended to the helium atom by the aid of configuration interaction. Bethe logarithms for $n\,^{1,3}S\,\,(n=1-8)$ states are calculated in three different gauges. The acceleration-gauge value of 4.364 036 7(2) a.u. and the velocity-acceleration gauge value of 4.364 036 4(2) a.u. for the $2\,^3S$ state are achieved, which are both at the 10$^{-7}$ level of accuracy. For other triplet states, present results in three different gauges all have five to seven accurate figures. While for the singlet states, the best convergent values are obtained in the pa-gauge, of which the numerical precision is also at the $10^{-5}$ to $10^{-7}$ level of accuracy.


title: Prevention and Mitigation of Catastrophic Failures in Demand-Supply Interdependent Networks

abstract: We propose a generic system model for a special category of interdependent networks, demand-supply networks, in which the demand and the supply nodes are associated with heterogeneous loads and resources, respectively. Our model sheds a light on a unique cascading failure mechanism induced by resource/load fluctuations, which in turn opens the door to conducting stress analysis on interdependent networks. Compared to the existing literature mainly concerned with the node connectivity, we focus on developing effective resource allocation methods to prevent these cascading failures from happening and to mitigate/confine them upon occurrence in the network. To prevent cascading failures, we identify some dangerous stress mechanisms, based on which we quantify the robustness of the network in terms of the resource configuration scheme. Afterward, we identify the optimal resource configuration under two resource/load fluctuations scenarios: uniform and proportional fluctuations. We further investigate the optimal resource configuration problem considering heterogeneous resource sharing costs among the nodes. To mitigate/confine ongoing cascading failures, we propose two network adaptations mechanisms: intentional failure and resource re-adjustment, based on which we propose an algorithm to mitigate an ongoing cascading failure while reinforcing the surviving network with a high robustness to avoid further failures.


title: RAIRE: Risk-Limiting Audits for IRV Elections

abstract: Risk-limiting post election audits guarantee a high probability of correcting incorrect election results, independent of why the result was incorrect. Ballot-polling audits select ballots at random and interpret those ballots as evidence for and against the reported result, continuing this process until either they support the recorded result, or they fall back to a full manual recount. For elections with digitised scanning and counting of ballots, a comparison audit compares randomly selected digital ballots with their paper versions. Discrepancies are referred to as errors, and are used to build evidence against or in support of the recorded result. Risk-limiting audits for first-past-the-post elections are well understood, and used in some US elections. We define a number of approaches to ballot-polling and comparison risk-limiting audits for Instant Runoff Voting (IRV) elections. We show that for almost all real elections we found, we can perform a risk-limiting audit by looking at only a small fraction of the total ballots (assuming no errors were made in the tallying and distribution of votes).


title: Research topics in finite groups and vertex algebras

abstract: We suggest a few projects for studying vertex algebras with emphasis on finite group viewpoints.


title: Robust Contraction Analysis of Nonlinear Systems via Differential IQC

abstract: We present a new approach to verifying contraction and $L_2$-gain of uncertain nonlinear systems, extending the well-known method of integral quadratic constraints. The uncertain system consists of a feedback interconnection of a nonlinear nominal system and uncertainties satisfying differential integral quadratic constraints. A pointwise linear matrix inequality condition is formulated to verify the closed-loop differential $ L_2 $ gain, which can lead to global reference-independent $ L_2 $ gain performance of the nonlinear uncertain system. For a polynomial nominal system, the convex verification conditions can be solved via sum-of-squares programming. A simple computational example based on jet-engine surge with input delays illustrates the approach.


title: General Methods for Digital Quantum Simulation of Gauge Theories

abstract: A general scheme is presented for simulating gauge theories, with matter fields, on a digital quantum computer. A Trotterized time-evolution operator that respects gauge symmetry is constructed, and a procedure for obtaining time-separated, gauge-invariant operators is detailed. We demonstrate the procedure on small lattices, including the simulation of a 2+1D non-Abelian gauge theory.

categories: hep-lat cond-mat.str-el nucl-th quant-ph


title: Bidirectional Recurrent Models for Offensive Tweet Classification

abstract: In this paper we propose four deep recurrent architectures to tackle the task of offensive tweet detection as well as further classification into targeting and subject of said targeting. Our architectures are based on LSTMs and GRUs, we present a simple bidirectional LSTM as a baseline system and then further increase the complexity of the models by adding convolutional layers and implementing a split-process-merge architecture with LSTM and GRU as processors. Multiple pre-processing techniques were also investigated. The validation F1-score results from each model are presented for the three subtasks as well as the final F1-score performance on the private competition test set. It was found that model complexity did not necessarily yield better results. Our best-performing model was also the simplest, a bidirectional LSTM; closely followed by a two-branch bidirectional LSTM and GRU architecture.


title: Rotating black string in dRGT massive gravity

abstract: One of the solutions of Einstein field equation with cylindrical symmetry is known as black string solution. In this work, the rotating black string solution in dRGT massive gravity is obtained and then called rotating-dRGT black string solution. This solution is a kind of generalized version of rotating AdS/dS black sting solution containing an additional two more terms characterizing the structure of graviton mass. The horizon structures of the black string are explored. The thermodynamical properties of the black string are investigated. We found that it is possible to obtain the Hawking-Page phase transition depending on the additional structure of the graviton mass, while it is not possible for rotating-AdS/dS black string. By analyzing the free energy, we also found that the stable rotating black string is bigger than the non-rotating one.


title: Link prediction with hyperbolic geometry

abstract: Link prediction is a paradigmatic problem in network science with a variety of applications. In latent space network models this problem boils down to ranking pairs of nodes in the order of increasing latent distances between them. The network model with hyperbolic latent spaces has a number of attractive properties suggesting it must be a powerful tool to predict links, but the past work in this direction reported mixed results. Here we perform systematic investigation of the utility of latent hyperbolic geometry for link prediction in networks. We first show that some measures of link prediction accuracy are extremely sensitive with respect to inaccuracies in the inference of latent hyperbolic coordinates of nodes, so that we develop a new coordinate inference method that maximizes the accuracy of such inference. Applying this method to synthetic and real networks, we then find that while there exists a multitude of competitive methods to predict obvious easy-to-predict links, among which hyperbolic link prediction is rarely the best but often competitive, it is the best, often by far, when the task is to predict less obvious missing links that are really hard to predict. These links include missing links in incomplete networks with large fractions of missing links, missing links between nodes that do not have any common neighbors, and missing links between dissimilar nodes at large latent distances. Overall these results suggest that the harder a specific link prediction task is, the more seriously one should consider using hyperbolic geometry.


title: Networks for Joint Affine and Non-parametric Image Registration

abstract: We introduce an end-to-end deep-learning framework for 3D medical image registration. In contrast to existing approaches, our framework combines two registration methods: an affine registration and a vector momentum-parameterized stationary velocity field (vSVF) model. Specifically, it consists of three stages. In the first stage, a multi-step affine network predicts affine transform parameters. In the second stage, we use a Unet-like network to generate a momentum, from which a velocity field can be computed via smoothing. Finally, in the third stage, we employ a self-iterable map-based vSVF component to provide a non-parametric refinement based on the current estimate of the transformation map. Once the model is trained, a registration is completed in one forward pass. To evaluate the performance, we conducted longitudinal and cross-subject experiments on 3D magnetic resonance images (MRI) of the knee of the Osteoarthritis Initiative (OAI) dataset. Results show that our framework achieves comparable performance to state-of-the-art medical image registration approaches, but it is much faster, with a better control of transformation regularity including the ability to produce approximately symmetric transformations, and combining affine and non-parametric registration.


title: Magnetic and electrical anisotropy with correlation and orbital effects in dimerized honeycomb ruthenate Li$_2$RuO$_3$

abstract: Li2RuO3 undergoes a structural transition at a relatively high temperature of 550 K with a distinct dimerization of Ru-Ru bonds on the otherwise isotropic honeycomb lattice. It exhibits a unique herringbone dimerization pattern with a largest ever reported value of the bond shrinkage of about ~ 0.5 \r{A}. Despite extensive studies, both theoretical and experimental, however, its origin and its effect on physical properties still remain to be understood. In this work, using high quality single crystals we investigated the anisotropy of resistivity ($\rho$) and magnetic susceptibility ($\chi$) to find a very clear anisotropy: $\rho$$_c*$ > $\rho$$_b$ > $\rho$$_a$ and $\chi$$_b$ > $\chi$$_a$ > $\chi$$_c*$. For possible theoretical interpretations, we carried out density functional calculations to conclude that these anisotropic behavior is due to the correlation effects combined with the unique orbital structure and the dimerization of Ru 4d bands.


title: Stellar mass as a galaxy cluster mass proxy: application to the Dark Energy Survey redMaPPer clusters

abstract: We introduce a galaxy cluster mass observable, $\mu_\star$, based on the stellar masses of cluster members, and we present results for the Dark Energy Survey (DES) Year 1 observations. Stellar masses are computed using a Bayesian Model Averaging method, and are validated for DES data using simulations and COSMOS data. We show that $\mu_\star$ works as a promising mass proxy by comparing our predictions to X-ray measurements. We measure the X-ray temperature-$\mu_\star$ relation for a total of 150 clusters matched between the wide-field DES Year 1 redMaPPer catalogue, and Chandra and XMM archival observations, spanning the redshift range $0.1<z<0.7$. For a scaling relation which is linear in logarithmic space, we find a slope of $\alpha = 0.488\pm0.043$ and a scatter in the X-ray temperature at fixed $\mu_\star$ of $\sigma_{{\rm ln} T_X|\mu_\star}=0.266^{+0.019}_{-0.020}$ for the joint sample. By using the halo mass scaling relations of the X-ray temperature from the Weighing the Giants program, we further derive the $\mu_\star$-conditioned scatter in mass, finding $\sigma_{{\rm ln} M|\mu_\star}=0.26^{+ 0.15}_{- 0.10}$. These results are competitive with well-established cluster mass proxies used for cosmological analyses, showing that $\mu_\star$ can be used as a reliable and physically motivated mass proxy to derive cosmological constraints.


title: Prostate Segmentation from Ultrasound Images using Residual Fully Convolutional Network

abstract: Medical imaging based prostate cancer diagnosis procedure uses intra-operative transrectal ultrasound (TRUS) imaging to visualize the prostate shape and location to collect tissue samples. Correct tissue sampling from prostate requires accurate prostate segmentation in TRUS images. To achieve this, this study uses a novel residual connection based fully convolutional network. The advantage of this segmentation technique is that it requires no pre-processing of TRUS images to perform the segmentation. Thus, it offers a faster and straightforward prostate segmentation from TRUS images. Results show that the proposed technique can achieve around 86% Dice Similarity accuracy using only few TRUS datasets.


title: Seiberg-Witten equation on a manifold with rank $2$-foliation

abstract: Let $M$ be a closed oriented $4$-manifold admitting a rank-$2$ oriented foliation with a metric of leafwise positive scalar curvature. If $b^+>1$, then we will show that the Seiberg-Witten invariant vanishes for all \spinc structures.


title: Empirical Evaluations of Seed Set Selection Strategies for Predictive Coding

abstract: Training documents have a significant impact on the performance of predictive models in the legal domain. Yet, there is limited research that explores the effectiveness of the training document selection strategy - in particular, the strategy used to select the seed set, or the set of documents an attorney reviews first to establish an initial model. Since there is limited research on this important component of predictive coding, the authors of this paper set out to identify strategies that consistently perform well. Our research demonstrated that the seed set selection strategy can have a significant impact on the precision of a predictive model. Enabling attorneys with the results of this study will allow them to initiate the most effective predictive modeling process to comb through the terabytes of data typically present in modern litigation. This study used documents from four actual legal cases to evaluate eight different seed set selection strategies. Attorneys can use the results contained within this paper to enhance their approach to predictive coding.


title: Dual Residual Networks Leveraging the Potential of Paired Operations for Image Restoration

abstract: In this paper, we study design of deep neural networks for tasks of image restoration. We propose a novel style of residual connections dubbed "dual residual connection", which exploits the potential of paired operations, e.g., up- and down-sampling or convolution with large- and small-size kernels. We design a modular block implementing this connection style; it is equipped with two containers to which arbitrary paired operations are inserted. Adopting the "unraveled" view of the residual networks proposed by Veit et al., we point out that a stack of the proposed modular blocks allows the first operation in a block interact with the second operation in any subsequent blocks. Specifying the two operations in each of the stacked blocks, we build a complete network for each individual task of image restoration. We experimentally evaluate the proposed approach on five image restoration tasks using nine datasets. The results show that the proposed networks with properly chosen paired operations outperform previous methods on almost all of the tasks and datasets.


title: Contingency Model Predictive Control for Automated Vehicles

abstract: We present Contingency Model Predictive Control (CMPC), a novel and implementable control framework which tracks a desired path while simultaneously maintaining a contingency plan -- an alternate trajectory to avert an identified potential emergency. In this way, CMPC anticipates events that might take place, instead of reacting when emergencies occur. We accomplish this by adding an additional prediction horizon in parallel to the classical receding MPC horizon. The contingency horizon is constrained to maintain a feasible avoidance solution; as such, CMPC is selectively robust to this emergency while tracking the desired path as closely as possible. After defining the framework mathematically, we demonstrate its effectiveness experimentally by comparing its performance to a state-of-the-art deterministic MPC. The controllers drive an automated research platform through a left-hand turn which may be covered by ice. Contingency MPC prepares for the potential loss of friction by purposefully and intuitively deviating from the prescribed path to approach the turn more conservatively; this deviation significantly mitigates the consequence of encountering ice.


title: On Measures and Measurements: a Fibre Bundle approach to Contextuality

abstract: Contextuality is the failure of "local" probabilistic models to become global ones. In this paper we introduce the notions of \emph{measurable fibre bundles}, \emph{probability fibre bundles}, and \emph{sample fibre bundle} which capture and make precise the former statement. The central notions of contextuality are discussed under this formalism, examples worked out, and some new aspects pointed out.


title: Stabilization and heteroepitaxial growth of metastable tetragonal FeS thin films by pulsed laser deposition

abstract: Pulsed laser deposition, a non-equilibrium thin-film growth technique, was used to stabilize metastable tetragonal iron sulfide (FeS), the bulk state of which is known as a superconductor with a critical temperature of 4 K. Comprehensive experiments revealed four important factors to stabilize tetragonal FeS epitaxial thin films: (i) an optimum growth temperature of 300 {\deg}C followed by thermal quenching, (ii) an optimum growth rate of ~7 nm/min, (iii) use of a high-purity bulk target, and (iv) use of a single-crystal substrate with small in-plane lattice mismatch (CaF2). Electrical resistivity measurements indicated that none of all the films exhibited superconductivity. Although an electric double-layer transistor structure was fabricated using the tetragonal FeS epitaxial film as a channel layer to achieve high-density carrier doping, no phase transition was observed. Possible reasons for the lack of superconductivity include lattice strain, off-stoichiometry of the film, electrochemical etching by the ionic liquid under gate bias, and surface degradation during device fabrication.


title: Analytical low-rank compression via proxy point selection

abstract: It has been known in potential theory that, for some kernels matrices corresponding to well-separated point sets, fast analytical low-rank approximation can be achieved via the use of proxy points. This proxy point method gives a surprisingly convenient way of explicitly writing out approximate basis matrices for a kernel matrix. However, this elegant strategy is rarely known or used in the numerical linear algebra community. It still needs clear algebraic understanding of the theoretical background. Moreover, rigorous quantifications of the approximation errors and reliable criteria for the selection of the proxy points are still missing. In this work, we use contour integration to clearly justify the idea in terms of a class of important kernels. We further provide comprehensive accuracy analysis for the analytical compression and show how to choose nearly optimal proxy points. The analytical compression is then combined with fast rank-revealing factorizations to get compact low-rank approximations and also to select certain representative points. We provide the error bounds for the resulting overall low-rank approximation. This work thus gives a fast and reliable strategy for compressing those kernel matrices. Furthermore, it provides an intuitive way of understanding the proxy point method and bridges the gap between this useful analytical strategy and practical low-rank approximations. Some numerical examples help to further illustrate the ideas.


title: Semantic Security on Wiretap Channels using Universal Hashing with Fading Applications

abstract: We furnish a procedure based on universal hash families (UHFs) that can convert an error correcting coding scheme (ECC) of rate $R$ into a semantically secure wiretap coding scheme of rate $R - \xi$ where $\xi$ is a parameter derived from the eavesdropper's point-to-point channel. This conversion is shown to be polynomial time efficient with block length and is applicable to any channel. When an ECC is chosen, our procedure induces a wiretap coding scheme that is concrete and efficient as long as the ECC is also such. To prove this induced wiretap coding scheme is semantically secure, we have constructed bounds on the information leaked to the eavesdropper. Our construction is an upgrade of bounds from recent literature: the novelty here being that our leakage bounds hold for any message distribution. Indeed, our wiretap procedure using UHFs and our characterization of its semantic leakage is the first main contribution of this work. The other main contribution of this work is as follows. We apply the aforementioned procedure to a variety of wiretap channels in order to show the procedure's efficacy, and as a result of such applications, our results establish new achievable semantically secure rates.


title: Finite size corrections at the hard edge for the Laguerre $\beta$ ensemble

abstract: A fundamental question in random matrix theory is to quantify the optimal rate of convergence to universal laws. We take up this problem for the Laguerre $\beta$ ensemble, characterised by the Dyson parameter $\beta$, and the Laguerre weight $x^a e^{-\beta x/2}$, $x > 0$ in the hard edge limit. The latter relates to the eigenvalues in the vicinity of the origin in the scaled variable $x \mapsto x/4N$. Previous work has established the corresponding functional form of various statistical quantities --- for example the distribution of the smallest eigenvalue, provided that $a \in \mathbb Z_{\ge 0}$. We show, using the theory of multidimensional hypergeometric functions based on Jack polynomials, that with the modified hard edge scaling $x \mapsto x/4(N+a/\beta)$, the rate of convergence to the limiting distribution is $O(1/N^2)$, which is optimal. In the case $\beta = 2$, general $a> -1$ the explicit functional form of the distribution of the smallest eigenvalue at this order can be computed, as it can for $a=1$ and general $\beta > 0$. An iterative scheme is presented to numerically approximate the functional form for general $a \in \mathbb Z_{\ge 2}$.


title: Perfect codes from PGL(2,5) in Star graphs

abstract: The Star graph $S_n$ is the Cayley graph of the symmetric group $Sym_n$ with the generating set $\{(1\mbox{ }i): 2\leq i\leq n \}$. Arumugam and Kala proved that $\{\pi\in Sym_n: \pi(1)=1\}$ is a perfect code in $S_n$ for any $n, n\geq 3$. In this note we show that for any $n, n\geq 6$ the Star graph $S_n$ contains a perfect code which is a union of cosets of the embedding of $PGL(2,5)$ into $Sym_6$.


title: Onset of giant planet migration before 4480 million years ago

abstract: Immediately after their formation, the terrestrial planets experienced intense impact bombardment by comets, leftover planetesimals from primary accretion, and asteroids. This temporal interval in solar system evolution, termed late accretion, thermally and chemically modified solid planetary surfaces and may have impeded the emergence of life on the Hadean Earth. The sources and tempo of late accretion are, however, vague. Here, we present a timeline that relates variably retentive radiometric ages from asteroidal meteorites, to new dynamical models of late accretion that invokes giant planet migration. Reconciliation of the geochronological data with dynamical models shows that giant planet migration immediately leads to an intense 30 Myr influx of comets to the entire solar system. The absence of whole-sale crustal reset ages after 4450 Ma for the most resilient chronometers from Earth, Moon, Mars, Vesta and various meteorite parent bodies confines the onset of giant planet migration to no later than ca. 4480 Ma. Waning impacts from planetesimals, asteroids (and a minor cometary component) continue to strike the inner planets through a protracted monotonic decline in impactor flux; this is in agreement with predictions from crater chronology. Amended global 3-D thermal analytical bombardment models derived from our new impact mass-production functions show that persistent niches for prebiotic chemistry on the early Hadean Earth could endure late accretion for at least the last 4400 Myr.


title: On Preempting Advanced Persistent Threats Using Probabilistic Graphical Models

abstract: This paper presents PULSAR, a framework for pre-empting Advanced Persistent Threats (APTs). PULSAR employs a probabilistic graphical model (specifically a Factor Graph) to infer the time evolution of an attack based on observed security events at runtime. PULSAR (i) learns the statistical significance of patterns of events from past attacks; (ii) composes these patterns into FGs to capture the progression of the attack; and (iii) decides on preemptive actions. PULSAR's accuracy and its performance are evaluated in three experiments at SystemX: (i) a study with a dataset containing 120 successful APTs over the past 10 years (PULSAR accurately identifies 91.7%); (ii) replaying of a set of ten unseen APTs (PULSAR stops 8 out of 10 replayed attacks before system integrity violation, and all ten before data exfiltration); and (iii) a production deployment of PULSAR (during a month-long deployment, PULSAR took an average of one second to make a decision).


title: What's in a frequency: new tools for graph Fourier Transform visualization

abstract: Recent progress in graph signal processing (GSP) has addressed a number of problems, including sampling and filtering. Proposed methods have focused on generic graphs and defined signals with certain characteristics, e.g., bandlimited signals, based on t he graph Fourier transform (GFT). However, the effect of GFT properties (e.g., vertex localization) on the behavior of such methods is not as well understood. In this paper, we propose novel GFT visualization tools and provide some examples to illustrate certain GFT properties and their impact on sampling or wavelet transforms.


title: Convolutional Neural Network on Semi-Regular Triangulated Meshes and its Application to Brain Image Data

abstract: We developed a convolution neural network (CNN) on semi-regular triangulated meshes whose vertices have 6 neighbours. The key blocks of the proposed CNN, including convolution and down-sampling, are directly defined in a vertex domain. By exploiting the ordering property of semi-regular meshes, the convolution is defined on a vertex domain with strong motivation from the spatial definition of classic convolution. Moreover, the down-sampling of a semi-regular mesh embedded in a 3D Euclidean space can achieve a down-sampling rate of 4, 16, 64, etc. We demonstrated the use of this vertex-based graph CNN for the classification of mild cognitive impairment (MCI) and Alzheimer's disease (AD) based on 3169 MRI scans of the Alzheimer's Disease Neuroimaging Initiative (ADNI). We compared the performance of the vertex-based graph CNN with that of the spectral graph CNN.


title: Exact slice sampler for Hierarchical Dirichlet Processes

abstract: We propose an exact slice sampler for Hierarchical Dirichlet process (HDP) and its associated mixture models (Teh et al., 2006). Although there are existing MCMC algorithms for sampling from the HDP, a slice sampler has been missing from the literature. Slice sampling is well-known for its desirable properties including its fast mixing and its natural potential for parallelization. On the other hand, the hierarchical nature of HDPs poses challenges to adopting a full-fledged slice sampler that automatically truncates all the infinite measures involved without ad-hoc modifications. In this work, we adopt the powerful idea of Bayesian variable augmentation to address this challenge. By introducing new latent variables, we obtain a full factorization of the joint distribution that is suitable for slice sampling. Our algorithm has several appealing features such as (1) fast mixing; (2) remaining exact while allowing natural truncation of the underlying infinite-dimensional measures, as in (Kalli et al., 2011), resulting in updates of only a finite number of necessary atoms and weights in each iteration; and (3) being naturally suited to parallel implementations. The underlying principle for joint factorization of the full likelihood is simple and can be applied to many other settings, such as designing sampling algorithms for general dependent Dirichlet process (DDP) models.


title: The Blueshift Of Civ Broad Emission Line In Qsos

abstract: For the sample from Ge et al. of 87 low-$z$ Palomar--Green (PG) quasi-stellar objects (QSOs) and 130 high-$z$ QSOs ($0<z<5$) with $\hb$-based single-epoch supermassive black hole (SMBH) masses, we performed a uniform decomposition of the \civ\ $\lambda$1549 broad-line profile. Based on the rest frame defined by the \oiii $\lambda$5007 narrow emission line, a medium-strong positive correlation is found between the \civ\ blueshift and the luminosity at 5100\AA\ or the Eddington ratio \leddR. A medium-strong negative relationship is found between the \civ\ blueshift and \civ\ equivalent width. These results support the postulation where the radiation pressure may be the driver of \civ\ blueshift. There is a medium strong correlation between the mass ratio of \civ-based to $\hb$-based \mbh and the \civ\ blueshift, which indicates that the bias for \civ-based \mbh is affected by the \civ\ profile.


title: Non-target Structural Displacement Measurement Using Reference Frame Based Deepflow

abstract: Structural displacement is crucial for structural health monitoring, although it is very challenging to measure in field conditions. Most existing displacement measurement methods are costly, labor intensive, and insufficiently accurate for measuring small dynamic displacements. Computer vision (CV) based methods incorporate optical devices with advanced image processing algorithms to accurately, cost-effectively, and remotely measure structural displacement with easy installation. However, non-target based CV methods are still limited by insufficient feature points, incorrect feature point detection, occlusion, and drift induced by tracking error accumulation. This paper presents a reference frame based Deepflow algorithm integrated with masking and signal filtering for non-target based displacement measurements. The proposed method allows the user to select points of interest for images with a low gradient for displacement tracking and directly calculate displacement without drift accumulated by measurement error. The proposed method is experimentally validated on a cantilevered beam under ambient and occluded test conditions. The accuracy of the proposed method is compared with that of a reference laser displacement sensor for validation. The significant advantage of the proposed method is its flexibility in extracting structural displacement in any region on structures that do not have distinct natural features.


title: EHT-HOPS pipeline for millimeter VLBI data reduction

abstract: We present the design and implementation of an automated data calibration and reduction pipeline for very-long-baseline interferometric (VLBI) observations taken at millimeter wavelengths. These short radio-wavelengths provide the best imaging resolution available from ground-based VLBI networks such as the Event Horizon Telescope (EHT) and the Global Millimeter VLBI Array (GMVA), but require specialized processing due to the strong effects from atmospheric opacity and turbulence as well as the heterogeneous nature of existing global arrays. The pipeline builds upon a calibration suite (HOPS) originally designed for precision geodetic VLBI. To support the reduction of data for astronomical observations, we have developed an additional framework for global phase and amplitude calibration which provides output in a standard data format for astronomical imaging and analysis. We test the pipeline on observations taken at 3.5 mm (86 GHz) by the GMVA joined by the phased Atacama Large Millimeter/submillimeter Array in April 2017, and demonstrate the benefits from the specialized processing of high frequency VLBI data with respect to classical analysis techniques.


title: Enhanced photoresponse in MoTe2 photodetectors with asymmetric graphene contacts

abstract: Atomically thin two dimensional (2D) materials are promising candidates for miniaturized high-performance optoelectronic devices. Here, we report on multilayer MoTe2 photodetectors contacted with asymmetric electrodes based on n- and p-type graphene layers. The asymmetry in the graphene contacts creates a large (Ebi ~100 kV cm-1) built-in electric field across the short (l = 15 nm) MoTe2 channel, causing a high and broad (400 to 1400 nm) photoresponse even without any externally applied voltage. Spatially resolved photovoltage maps reveal an enhanced photoresponse and larger built-in electric field in regions of the MoTe2 layer between the two graphene contacts. Furthermore, a fast (~0.01 ms) photoresponse is achieved in both the photovoltaic and photoconductive operation modes of the junction. Our findings could be extended to other 2D materials and offer prospects for the implementation of asymmetric graphene contacts in future low-power optoelectronic applications.


title: Exterior powers in Iwasawa theory

abstract: The Iwasawa theory of CM fields has traditionally concerned Iwasawa modules that are abelian pro-p Galois groups with ramification allowed at a maximal set of primes over p such that the module is torsion. A main conjecture for such an Iwasawa module describes its codimension one support in terms of a p-adic L-function attached to the primes of ramification. In this paper, we study more general and potentially much smaller Iwasawa modules that are quotients of exterior powers of Iwasawa modules with ramification at a set of primes over p by sums of exterior powers of inertia subgroups. We show that the higher codimension support of such quotients can be measured by finite collections of p-adic L-functions under the relevant CM main conjectures.


title: Ultra Low-Power System for Remote ECG Monitoring

abstract: A complete system solution extracting signals from the patient chest with three leads including motion artifact removal in both analog and digital implementations are described. The resulting ECG signal is transferred via Bluetooth low energy to a mobile phone. Using deep sleep modes, the overall power consumption is less than 300 $\mu$A and the device can operate for more than 20 days using a 150mAh battery. The screening software looks for suspicious traces such as those with missing pulses, tachycardia, bradycardia, etc. The mobile phone software also eliminates any remaining motion artifact. The traces are subsequently processed in detail in a cloud server and to a physicians dashboard for long-term monitoring.


title: Towards Robust Curve Text Detection with Conditional Spatial Expansion

abstract: It is challenging to detect curve texts due to their irregular shapes and varying sizes. In this paper, we first investigate the deficiency of the existing curve detection methods and then propose a novel Conditional Spatial Expansion (CSE) mechanism to improve the performance of curve text detection. Instead of regarding the curve text detection as a polygon regression or a segmentation problem, we treat it as a region expansion process. Our CSE starts with a seed arbitrarily initialized within a text region and progressively merges neighborhood regions based on the extracted local features by a CNN and contextual information of merged regions. The CSE is highly parameterized and can be seamlessly integrated into existing object detection frameworks. Enhanced by the data-dependent CSE mechanism, our curve text detection system provides robust instance-level text region extraction with minimal post-processing. The analysis experiment shows that our CSE can handle texts with various shapes, sizes, and orientations, and can effectively suppress the false-positives coming from text-like textures or unexpected texts included in the same RoI. Compared with the existing curve text detection algorithms, our method is more robust and enjoys a simpler processing flow. It also creates a new state-of-art performance on curve text benchmarks with F-score of up to 78.4$\%$.


title: Coalgebraic Geometric Logic: Basic Theory

abstract: Using the theory of coalgebra, we introduce a uniform framework for adding modalities to the language of propositional geometric logic. Models for this logic are based on coalgebras for an endofunctor on some full subcategory of the category of topological spaces and continuous functions. We investigate derivation systems, soundness and completeness for such geometric modal logics, and we specify a method of lifting an endofunctor on Set, accompanied by a collection of predicate liftings, to an endofunctor on the category of topological spaces, again accompanied by a collection of (open) predicate liftings. Furthermore, we compare the notions of modal equivalence, behavioural equivalence and bisimulation on the resulting class of models, and we provide a final object for the corresponding category.


title: 3D non-LTE line formation of neutral carbon in the Sun

abstract: Carbon abundances in late-type stars are important in a variety of astrophysical contexts. However C i lines, one of the main abundance diagnostics, are sensitive to departures from local thermodynamic equilibrium (LTE). We present a model atom for non-LTE analyses of C i lines, that uses a new, physically-motivated recipe for the rates of neutral hydrogen impact excitation. We analyse C i lines in the solar spectrum, employing a three-dimensional (3D) hydrodynamic model solar atmosphere and 3D non-LTE radiative transfer. We find negative non-LTE abundance corrections for C i lines in the solar photosphere, in accordance with previous studies, reaching up to around 0.1 dex in the disk-integrated flux. We also present the first fully consistent 3D non-LTE solar carbon abundance determination: we infer log $\epsilon_{\text{C}}$ = $8.44\pm0.02$, in good agreement with the current standard value. Our models reproduce the observed solar centre-to-limb variations of various C i lines, without any adjustments to the rates of neutral hydrogen impact excitation, suggesting that the proposed recipe may be a solution to the long-standing problem of how to reliably model inelastic collisions with neutral hydrogen in late-type stellar atmospheres.


title: Weakly-Supervised Discovery of Geometry-Aware Representation for 3D Human Pose Estimation

abstract: Recent studies have shown remarkable advances in 3D human pose estimation from monocular images, with the help of large-scale in-door 3D datasets and sophisticated network architectures. However, the generalizability to different environments remains an elusive goal. In this work, we propose a geometry-aware 3D representation for the human pose to address this limitation by using multiple views in a simple auto-encoder model at the training stage and only 2D keypoint information as supervision. A view synthesis framework is proposed to learn the shared 3D representation between viewpoints with synthesizing the human pose from one viewpoint to the other one. Instead of performing a direct transfer in the raw image-level, we propose a skeleton-based encoder-decoder mechanism to distil only pose-related representation in the latent space. A learning-based representation consistency constraint is further introduced to facilitate the robustness of latent 3D representation. Since the learnt representation encodes 3D geometry information, mapping it to 3D pose will be much easier than conventional frameworks that use an image or 2D coordinates as the input of 3D pose estimator. We demonstrate our approach on the task of 3D human pose estimation. Comprehensive experiments on three popular benchmarks show that our model can significantly improve the performance of state-of-the-art methods with simply injecting the representation as a robust 3D prior.


title: A Cold Stellar Stream in Pegasus

abstract: We report the serendipitous discovery of a stellar stream in the constellation Pegasus in the south Galactic hemisphere. The stellar stream was detected using the SDSS Data Release 14 by means of a matched filter in the color--magnitude diagram that is optimised for a stellar population that is 8 Gyr old with [Fe/H] = $-$0.46 dex, and located at heliocentric distance of 18 kpc. The candidate stream is faint (turnoff point at $r_0 \sim$ 19.6), sparse and barely visible in SDSS photometry. It is also detected in the (shallower) Pan-STARRs data. The residual stellar density in the $(u-g)_0$, $(g-r)_0$ color--color diagram gives the same estimate for the age and [Fe/H] of this stellar population. The stream is located at a Galactic coordinates $(l,b) = (79.4,-24.6)$ and extends over 9$^\circ$ (2.5 kpc), with a width of 112 pc. The narrow width suggests a globular cluster progenitor.


title: Some multiplicative equations in finite fields

abstract: In this paper we consider estimating the number of solutions to multiplicative equations in finite fields when the variables run through certain sets with high additive structure. In particular, we consider estimating the multiplicative energy of generalized arithmetic progressions in prime fields and of boxes in arbitrary finite fields and obtain sharp bounds in more general scenarios than previously known. Our arguments extend some ideas of Konyagin and Bourgain and Chang into new settings.


title: Gravitational Production of Superheavy Dark Matter and Associated Cosmological Signatures

abstract: We study the gravitational production of super-Hubble-mass dark matter in the very early universe. We first review the simplest scenario where dark matter is produced mainly during slow roll inflation. Then we move on to consider the cases where dark matter is produced during the transition period between inflation and the subsequent cosmological evolution. The limits of smooth and sudden transitions are studied, respectively. The relic abundances and the cosmological collider signals are calculated.


title: Ramond-Ramond fields and twisted differential K-theory

abstract: We provide a systematic approach to describing the Ramond-Ramond (RR) fields as elements in twisted differential K-theory. This builds on a series of constructions by the authors on geometric and computational aspects of twisted differential K-theory, which to a large extent were originally motivated by this problem. In addition to providing a new conceptual framework and a mathematically solid setting, this allows us to uncover interesting and novel effects. Explicitly, we use our recently constructed Atiyah-Hirzebruch spectral sequence (AHSS) for twisted differential K-theory to characterize the RR fields and their quantization, which involves interesting interplay between geometric and topological data. We illustrate this with the examples of spheres, tori, and Calabi-Yau threefolds.


title: Directional limits on persistent gravitational waves using data from Advanced LIGO's first two observing runs

abstract: We perform an unmodeled search for persistent, directional gravitational wave (GW) sources using data from the first and second observing runs of Advanced LIGO. We do not find evidence for any GW signals. We place limits on the broadband GW flux emitted at 25~Hz from point sources with a power law spectrum at $F_{\alpha,\Theta} <(0.05-25)\times 10^{-8} ~{\rm erg\,cm^{-2}\,s^{-1}\,Hz^{-1}}$ and the (normalized) energy density spectrum in GWs at 25 Hz from extended sources at $\Omega_{\alpha}(\Theta) <(0.19-2.89)\times 10^{-8} ~{\rm sr^{-1}}$ where $\alpha$ is the spectral index of the energy density spectrum. These represent improvements of $2.5-3\times$ over previous limits. We also consider point sources emitting GWs at a single frequency, targeting the directions of Sco X-1, SN 1987A, and the Galactic Center. The best upper limits on the strain amplitude of a potential source in these three directions range from $h_0 < (3.6-4.7)\times 10^{-25}$, 1.5$\times$ better than previous limits set with the same analysis method. We also report on a marginally significant outlier at 36.06~Hz. This outlier is not consistent with a persistent gravitational-wave source as its significance diminishes when combining all of the available data.


title: Transverse lines to surfaces over finite fields

abstract: We prove that if $S$ is a smooth reflexive surface in $\mathbb{P}^3$ defined over a finite field $\mathbb{F}_q$, then there exists an $\mathbb{F}_q$-line meeting $S$ transversely provided that $q\geq c\operatorname{deg}(S)$, where $c=\frac{3+\sqrt{17}}{4}\approx 1.7808$. Without the reflexivity hypothesis, we prove the existence of a transverse $\mathbb{F}_q$-line for $q\geq \operatorname{deg}(S)^2$.


title: The Kepler problem: polynomial algebra of non-polynomial first integrals

abstract: The sum of elliptic integrals simultaneously determines orbits in thr Kepler problem and the addition of divisors on elliptic curves. Periodic motion of a body in physical space is defined by symmetries, whereas periodic motion of divisors is defined by a fixed point on the curve. Algebra of the first integrals associated with symmetries is a well-known mathematical object, whereas algebra of the first integrals associated with coordinates of fixed points is unknown. In this paper, we discuss polynomial algebras of non-polynomial first integrals of superintegrable systems associated with elliptic curves.


title: Parametic Classification of Handvein Patterns Based on Texture Features

abstract: In this paper, we have developed Biometric recognition system adopting hand based modality Handvein, which has the unique pattern for each individual and it is impossible to counterfeit and fabricate as it is an internal feature. We have opted in choosing feature extraction algorithms such as LBP-visual descriptor ,LPQ-blur insensitive texture operator, Log-Gabor-Texture descriptor. We have chosen well known classifiers such as KNN and SVM for classification. We have experimented and tabulated results of single algorithm recognition rate for Handvein under different distance measures and kernel options. The feature level fusion is carried out which increased the performance level.


title: Exploring neutrino mass and mass hierarchy in interacting dark energy models

abstract: We investigate how the dark energy properties impact the constraints on the total neutrino mass in interacting dark energy (IDE) models. In this study, we focus on two typical interacting dynamical dark energy models, i.e., the interacting $w$ cold dark matter (I$w$CDM) model and the interacting holographic dark energy (IHDE) model. To avoid the large-scale instability problem in IDE models, we apply the parameterized post-Friedmann approach to calculate the perturbation of dark energy. We employ the Planck 2015 cosmic microwave background temperature and polarization data, combined with low-redshift measurements on baryon acoustic oscillation distance scales, type Ia supernovae, and the Hubble constant, to constrain the cosmological parameters. We find that the dark energy properties could influence the constraint limits on the total neutrino mass. Once dynamical dark energy is considered in the IDE models, the upper bounds of $\sum m_\nu$ will be changed. By considering the values of $\chi^2_{\rm min}$, we find that in these IDE models the normal hierarchy case is slightly preferred over the inverted hierarchy case; for example, $\Delta\chi^2=2.720$ is given in the IHDE+$\sum m_\nu$ model. In addition, we also find that in the I$w$CDM+$\sum m_\nu$ model $\beta=0$ is consistent with current observational data inside the 1$\sigma$ range, and in the IHDE+$\sum m_\nu$ model $\beta>0$ is favored at more than 2$\sigma$ level.


title: Neural Networks in Hybrid Precoding for Millimeter Wave Massive MIMO Systems

abstract: Neural networks have been applied to the physical layer of wireless communication systems to solve complex problems. In millimeter wave (mmWave) massive multiple-input multiple-output (MIMO) systems, hybrid precoding has been considered as an energy-efficient technology to replace fully-digital precoding. The way of designing hybrid precoding in mmWave massive MIMO systems by multi-layer neural networks has not been investigated. Based on further decomposing the baseband precoding matrix, an idea is proposed in this paper to map hybrid precoding structure to a multi-layer neural network. Considering the deterioration in the throughput and energy efficiency of mmWave massive MIMO systems, the feasibility of the proposed idea is analyzed. Moreover, a singular value decomposition (SVD) based decomposing (SVDDE) algorithm is proposed to evaluate the feasibility of the proposed idea. Simulation results indicate that there is an optimal number of users which can minimize the performance deterioration. Moreover, the simulation results also show that slight deterioration in the throughput and energy efficiency of mmWave massive MIMO systems is caused by further decomposing the baseband precoding matrix. In other words, further decomposing the baseband precoding matrix is a feasible way to map the hybrid precoding structure to a multi-layer neural network.


title: Stochastic Optimization of Sorting Networks via Continuous Relaxations

abstract: Sorting input objects is an important step in many machine learning pipelines. However, the sorting operator is non-differentiable with respect to its inputs, which prohibits end-to-end gradient-based optimization. In this work, we propose NeuralSort, a general-purpose continuous relaxation of the output of the sorting operator from permutation matrices to the set of unimodal row-stochastic matrices, where every row sums to one and has a distinct arg max. This relaxation permits straight-through optimization of any computational graph involve a sorting operation. Further, we use this relaxation to enable gradient-based stochastic optimization over the combinatorially large space of permutations by deriving a reparameterized gradient estimator for the Plackett-Luce family of distributions over permutations. We demonstrate the usefulness of our framework on three tasks that require learning semantic orderings of high-dimensional objects, including a fully differentiable, parameterized extension of the k-nearest neighbors algorithm.


title: On the structure of Hamiltonian impact systems

abstract: Near-integrability is usually associated with smooth small perturbations of smooth integrable systems. Studying integrable mechanical Hamiltonian flows with impacts that respect the symmetries of the integrable structure provides an additional rich class of non-smooth systems that can be analyzed. Such systems exhibit rich dynamics, as, in addition to the underlying integrable structure, some of the trajectories may undergo only transverse impact, others may undergo also tangential impacts, and some trajectories do not impact at all. Under perturbations, each of these classes of orbits behaves differently. Tools for classifying these different types of dynamics in 2 degrees-of-freedom Hamiltonian impact systems with underlying separable integrable dynamics are presented. Moreover, some of these tools may be extended to far from integrable cases. In particular, a generalization of the energy momentum bifurcation diagram, Fomenko graphs and the hierarchy of bifurcations framework to impact systems is constructed. It is shown that such representations classify dynamically different regions in phase space. For the integrable and near integrable (small perturbations) cases these provide global information on the dynamics whereas for the far from integrable (non-perturbative) regimes, these provide rough classification of the first impact map. The interpretation of these results in terms of projections of solutions to the configuration space as well as the relations of these to the Hill region are presented.


title: A novel energy factorization approach for the diffuse-interface model with Peng-Robinson equation of state

abstract: The Peng-Robinson equation of state (PR-EoS) has become one of the most extensively applied equations of state in chemical engineering and petroleum industry due to its excellent accuracy in predicting the thermodynamic properties of a wide variety of materials, especially hydrocarbons. Although great efforts have been made to construct efficient numerical methods for the diffuse interface models with PR-EoS, there is still not a linear numerical scheme that can be proved to preserve the original energy dissipation law. In order to pursue such a numerical scheme, we propose a novel energy factorization (EF) approach, which first factorizes an energy function into a product of several factors and then treats the factors using their properties to obtain the semi-implicit linear schemes. We apply the EF approach to deal with the Helmholtz free energy density determined by PR-EoS, and then propose a linear semi-implicit numerical scheme that inherits the original energy dissipation law. Moreover, the proposed scheme is proved to satisfy the maximum principle in both the time semi-discrete form and the cell-centered finite difference fully discrete form under certain conditions. Numerical results are presented to demonstrate the stability and efficiency of the proposed scheme.


title: A convex programming approach for discrete-time Markov decision processes under the expected total reward criterion

abstract: In this work, we study discrete-time Markov decision processes (MDPs) under constraints with Borel state and action spaces and where all the performance functions have the same form of the expected total reward (ETR) criterion over the infinite time horizon. One of our objective is to propose a convex programming formulation for this type of MDPs. It will be shown that the values of the constrained control problem and the associated convex program coincide and that if there exists an optimal solution to the convex program then there exists a stationary randomized policy which is optimal for the MDP. It will be also shown that in the framework of constrained control problems, the supremum of the expected total rewards over the set of randomized policies is equal to the supremum of the expected total rewards over the set of stationary randomized policies. We consider standard hypotheses such as the so-called continuity-compactness conditions and a Slater-type condition. Our assumptions are quite weak to deal with cases that have not yet been addressed in the literature. An example is presented to illustrate our results with respect to those of the literature.


title: Manifold constrained non-uniformly elliptic problems

abstract: We consider the problem of minimizing variational integrals defined on \cc{nonlinear} Sobolev spaces of competitors taking values into the sphere. The main novelty is that the underlying energy features a non-uniformly elliptic integrand exhibiting different polynomial growth conditions and no homogeneity. We develop a few intrinsic methods aimed at proving partial regularity of minima and providing techniques for treating larger classes of similar constrained non-uniformly elliptic variational problems. In order to give estimates for the singular sets we use a general family of Hausdorff type measures following the local geometry of the integrand. A suitable comparison is provided with respect to the naturally associated capacities.