title: The central spheroids of Milky Way mass-sized galaxies

abstract: We study the properties of the central spheroids located within 10 kpc of the centre of mass of Milky Way mass-sized galaxies simulated in a cosmological context. The simulated central regions are dominated by stars older than 10 Gyr, mostly formed in situ, with a contribution of ~30 per cent from accreted stars. These stars formed in well-defined starbursts, although accreted stars exhibit sharper and earlier ones. The fraction of accreted stars increases with galactocentric distance, so that at a radius of ~8-10 kpc a fraction of ~40 per cent, on average, are detected. Accreted stars are slightly younger, lower metallicity, and more $\alpha$-enhanced than in situ stars. A significant fraction of old stars in the central regions come from a few ($2-3$) massive satellites ($\sim 10^{10}{\rm M}_\odot$). The bulge components receive larger contributions of accreted stars formed in dwarfs smaller than $\sim 10^{9.5}{\rm M}_\odot$. The difference between the distributions of ages and metallicities of old stars is thus linked to the accretion histories -- those central regions with a larger fraction of accreted stars are those with contributions from more massive satellites. The kinematical properties of in situ and accreted stars are consistent with the latter being supported by their velocity dispersions, while the former exhibit clear signatures of rotational support. Our simulations demonstrate a range of characteristics, with some systems exhibiting a co-existing bar and spheroid in their central regions, resembling in some respect the central region of the Milky Way.


title: Optimal projection of observations in a Bayesian setting

abstract: Optimal dimensionality reduction methods are proposed for the Bayesian inference of a Gaussian linear model with additive noise in presence of overabundant data. Three different optimal projections of the observations are proposed based on information theory: the projection that minimizes the Kullback-Leibler divergence between the posterior distributions of the original and the projected models, the one that minimizes the expected Kullback-Leibler divergence between the same distributions, and the one that maximizes the mutual information between the parameter of interest and the projected observations. The first two optimization problems are formulated as the determination of an optimal subspace and therefore the solution is computed using Riemannian optimization algorithms on the Grassmann manifold. Regarding the maximization of the mutual information, it is shown that there exists an optimal subspace that minimizes the entropy of the posterior distribution of the reduced model; a basis of the subspace can be computed as the solution to a generalized eigenvalue problem; an a priori error estimate on the mutual information is available for this particular solution; and that the dimensionality of the subspace to exactly conserve the mutual information between the input and the output of the models is less than the number of parameters to be inferred. Numerical applications to linear and nonlinear models are used to assess the efficiency of the proposed approaches, and to highlight their advantages compared to standard approaches based on the principal component analysis of the observations.


title: High-dimensional posterior consistency for hierarchical non-local priors in regression

abstract: The choice of tuning parameters in Bayesian variable selection is a critical problem in modern statistics. In particular, for Bayesian linear regression with non-local priors, the scale parameter in the non-local prior density is an important tuning parameter which reflects the dispersion of the non-local prior density around zero, and implicitly determines the size of the regression coefficients that will be shrunk to zero. Current approaches treat the scale parameter as given, and suggest choices based on prior coverage/asymptotic considerations. In this paper, we consider the fully Bayesian approach introduced in (Wu, 2016) with the pMOM non-local prior and an appropriate Inverse-Gamma prior on the tuning parameter to analyze the underlying theoretical property. Under standard regularity assumptions, we establish strong model selection consistency in a high-dimensional setting, where $p$ is allowed to increase at a polynomial rate with n$or even at a sub-exponential rate with n. Through simulation studies, we demonstrate that our model selection procedure can outperform other Bayesian methods which treat the scale parameter as given, and commonly used penalized likelihood methods, in a range of simulation settings.


title: Clifford algebras and their applications to Lie groups and spinors

abstract: In these lectures, we discuss some well-known facts about Clifford algebras: matrix representations, Cartan's periodicity of 8, double coverings of orthogonal groups by spin groups, Dirac equation in different formalisms, spinors in $n$ dimensions, etc. We also present our point of view on some problems. Namely, we discuss the generalization of the Pauli theorem, the basic ideas of the method of averaging in Clifford algebras, the notion of quaternion type of Clifford algebra elements, the classification of Lie subalgebras of specific type in Clifford algebra, etc.


title: The Toric SO(10) F-Theory Landscape

abstract: Supergravity theories in more than four dimensions with grand unified gauge symmetries are an important intermediate step towards the ultraviolet completion of the Standard Model in string theory. Using toric geometry, we classify and analyze six-dimensional F-theory vacua with gauge group SO(10) taking into account Mordell-Weil U(1) and discrete gauge factors. We determine the full matter spectrum of these models, including charged and neutral SO(10) singlets. Based solely on the geometry, we compute all matter multiplicities and confirm the cancellation of gauge and gravitational anomalies independent of the base space. Particular emphasis is put on symmetry enhancements at the loci of matter fields and to the frequent appearance of superconformal points. They are linked to non-toric K\"ahler deformations which contribute to the counting of degrees of freedom. We compute the anomaly coefficients for these theories as well by using a base-independent blow-up procedure and superconformal matter transitions. Finally, we identify six-dimensional supergravity models which can yield the Standard Model with high-scale supersymmetry by further compactification to four dimensions in an Abelian flux background.


title: Yaglom limits for R-transient chains with non-trivial Martin boundary

abstract: We give conditions for the existence of a Yaglom limit for R-transient Markov chains with non-trivial rho-Martin entrance boundary (rho=1/R) and we characterize the rho-invariant limiting quasistationary distribution.


title: Demagnetizing fields in all-optical switching

abstract: Time-resolved pump-probe measurements show ultrafast and heat accumulation demagnetization in Co/Pd superlattices on glass substrates. A model of demagnetizing fields and micromagnetic simulations are applied to examine the evolution of a demagnetized cylinder into a switched state.


title: Finite searches, Chowla's cosine problem, and large Newman polynomials

abstract: A length $n$ cosine sum is an expression of the form $\cos a_1\theta + \cdots + \cos a_n\theta$ where $a_1 < \cdots < a_n$ are positive integers, and a length $n$ Newman polynomial is an expression of the form $z^{a_1} + \cdots + z^{a_n}$ where $a_1 < \cdots < a_n$ are nonnegative integers. We define $-\lambda(n)$ to be the largest minimum of a length $n$ cosine sum as $\{a_1,\ldots,a_n\}$ ranges over all sets of $n$ positive integers, and we define $\mu(n)$ to be the largest minimum modulus on the unit circle of a length $n$ Newman polynomial as $\{a_1,\ldots,a_n\}$ ranges over all sets of $n$ nonnegative integers. Since there are infinitely many possibilities for the $a_j$, it is not obvious how to compute $\lambda(n)$ or $\mu(n)$ for a given $n$ in finitely many steps. Campbell et al. found the value of $\mu(3)$ in 1983, and Goddard found the value of $\mu(4)$ in 1992. In this paper, we find the values of $\lambda(2)$ and $\lambda(3)$ and nontrivial bounds on $\mu(5)$. We also include further remarks on the seemingly difficult general task of reducing the computation of $\lambda(n)$ or $\mu(n)$ to a finite problem.


title: A Server-based Approach for Predictable GPU Access with Improved Analysis

abstract: We propose a server-based approach to manage a general-purpose graphics processing unit (GPU) in a predictable and efficient manner. Our proposed approach introduces a GPU server that is a dedicated task to handle GPU requests from other tasks on their behalf. The GPU server ensures bounded time to access the GPU, and allows other tasks to suspend during their GPU computation to save CPU cycles. By doing so, we address the two major limitations of the existing real-time synchronization-based GPU management approach: busy waiting and long priority inversion. We have implemented a prototype of the server-based approach on a real embedded platform. This case study demonstrates the practicality and effectiveness of the server-based approach. Experimental results indicate that the server-based approach yields significant improvements in task schedulability over the existing synchronization-based approach in most practical settings. Although we focus on a GPU in this paper, the server-based approach can also be used for other types of computational accelerators.


title: An Analog Neural Network Computing Engine using CMOS-Compatible Charge-Trap-Transistor (CTT)

abstract: An analog neural network computing engine based on CMOS-compatible charge-trap transistor (CTT) is proposed in this paper. CTT devices are used as analog multipliers. Compared to digital multipliers, CTT-based analog multiplier shows significant area and power reduction. The proposed computing engine is composed of a scalable CTT multiplier array and energy efficient analog-digital interfaces. Through implementing the sequential analog fabric (SAF), the engine mixed-signal interfaces are simplified and hardware overhead remains constant regardless of the size of the array. A proof-of-concept 784 by 784 CTT computing engine is implemented using TSMC 28nm CMOS technology and occupied 0.68mm2. The simulated performance achieves 76.8 TOPS (8-bit) with 500 MHz clock frequency and consumes 14.8 mW. As an example, we utilize this computing engine to address a classic pattern recognition problem -- classifying handwritten digits on MNIST database and obtained a performance comparable to state-of-the-art fully connected neural networks using 8-bit fixed-point resolution.


title: Permutational symmetries for coincidence rates in multi-mode multi-photonic interferometry

abstract: We obtain coincidence rates for passive optical interferometry by exploiting the permutational symmetries of partially distinguishable input photons, and our approach elucidates qualitative features of multi-photon coincidence landscapes. We treat the interferometer input as a product state of any number of photons in each input mode with photons distinguished by their arrival time. Detectors at the output of the interferometer count photons from each output mode over a long integration time. We generalize and prove the claim of Tillmann et al. [Phys. Rev. X 5 041015 (2015)] that coincidence rates can be elegantly expressed in terms of immanants. Immanants are functions of matrices that exhibit permutational symmetries and the immanants appearing in our coincidence-rate expressions share permutational symmetries with the input state. Our results are obtained by employing representation theory of the symmetric group to analyze systems of arbitrary number of photons in arbitrarily sized interferometers.


title: On Collaborative Compressive Sensing Systems: The Framework, Design and Algorithm

abstract: Based on the maximum likelihood estimation principle, we derive a collaborative estimation framework that fuses several different estimators and yields a better estimate. Applying it to compressive sensing (CS), we propose a collaborative CS (CCS) scheme consisting of a bank of $K$ CS systems that share the same sensing matrix but have different sparsifying dictionaries. This CCS system is expected to yield better performance than each individual CS system, while requiring the same time as that needed for each individual CS system when a parallel computing strategy is used. We then provide an approach to designing optimal CCS systems by utilizing a measure that involves both the sensing matrix and dictionaries and hence allows us to simultaneously optimize the sensing matrix and all the $K$ dictionaries. An alternating minimization-based algorithm is derived for solving the corresponding optimal design problem. With a rigorous convergence analysis, we show that the proposed algorithm is convergent. Experiments are carried out to confirm the theoretical results and show that the proposed CCS system yields significant improvements over the existing CS systems in terms of the signal recovery accuracy.


title: A PAC-Bayesian Analysis of Randomized Learning with Application to Stochastic Gradient Descent

abstract: We study the generalization error of randomized learning algorithms -- focusing on stochastic gradient descent (SGD) -- using a novel combination of PAC-Bayes and algorithmic stability. Importantly, our generalization bounds hold for all posterior distributions on an algorithm's random hyperparameters, including distributions that depend on the training data. This inspires an adaptive sampling algorithm for SGD that optimizes the posterior at runtime. We analyze this algorithm in the context of our generalization bounds and evaluate it on a benchmark dataset. Our experiments demonstrate that adaptive sampling can reduce empirical risk faster than uniform sampling while also improving out-of-sample accuracy.


title: Boundary regularity and non-transversal intersection for the fully nonlinear obstacle problem

abstract: In this paper non-transversal intersection of the free and fixed boundary is shown to hold in any dimension for obstacle problems generated by fully nonlinear uniformly elliptic operators. Moreover, $C^1$ regularity results of the free boundary are obtained and a classification of blow-up solutions is given.


title: On Sinc Quadrature Approximations of Fractional Powers of Regularly Accretive Operators

abstract: We consider the finite element approximation of fractional powers of regularly accretive operators via the Dunford-Taylor integral approach. We use a sinc quadrature scheme to approximate the Balakrishnan representation of the negative powers of the operator as well as its finite element approximation. We improve the exponentially convergent error estimates from [A. Bonito, J. E. Pasciak, IMA J. Numer. Anal. (2016) 00, 1-29] by reducing the regularity required on the data. Numerical experiments illustrating the new theory are provided.


title: Learning of Coordination Policies for Robotic Swarms

abstract: Inspired by biological swarms, robotic swarms are envisioned to solve real-world problems that are difficult for individual agents. Biological swarms can achieve collective intelligence based on local interactions and simple rules; however, designing effective distributed policies for large-scale robotic swarms to achieve a global objective can be challenging. Although it is often possible to design an optimal centralized strategy for smaller numbers of agents, those methods can fail as the number of agents increases. Motivated by the growing success of machine learning, we develop a deep learning approach that learns distributed coordination policies from centralized policies. In contrast to traditional distributed control approaches, which are usually based on human-designed policies for relatively simple tasks, this learning-based approach can be adapted to more difficult tasks. We demonstrate the efficacy of our proposed approach on two different tasks, the well-known rendezvous problem and a more difficult particle assignment problem. For the latter, no known distributed policy exists. From extensive simulations, it is shown that the performance of the learned coordination policies is comparable to the centralized policies, surpassing state-of-the-art distributed policies. Thereby, our proposed approach provides a promising alternative for real-world coordination problems that would be otherwise computationally expensive to solve or intangible to explore.


title: Localization in the Disordered Holstein model

abstract: The Holstein model describes the motion of a tight-binding tracer particle interacting with a field of quantum harmonic oscillators. We consider this model with an on-site random potential. Provided the hopping amplitude for the particle is small, we prove localization for matrix elements of the resolvent, in particle position and in the field Fock space. These bounds imply a form of dynamical localization for the particle position that leaves open the possibility of resonant tunneling in Fock space between equivalent field configurations.


title: Distributed Training Large-Scale Deep Architectures

abstract: Scale of data and scale of computation infrastructures together enable the current deep learning renaissance. However, training large-scale deep architectures demands both algorithmic improvement and careful system configuration. In this paper, we focus on employing the system approach to speed up large-scale training. Via lessons learned from our routine benchmarking effort, we first identify bottlenecks and overheads that hinter data parallelism. We then devise guidelines that help practitioners to configure an effective system and fine-tune parameters to achieve desired speedup. Specifically, we develop a procedure for setting minibatch size and choosing computation algorithms. We also derive lemmas for determining the quantity of key components such as the number of GPUs and parameter servers. Experiments and examples show that these guidelines help effectively speed up large-scale deep learning training.


title: Secure Beamforming in Full-Duplex SWIPT Systems

abstract: Physical layer security is a key issue in the full duplex (FD) communication systems due to the broadcast nature of wireless channels. In this paper, the joint design of information and artificial noise beamforming vectors is proposed for the FD simultaneous wireless information and power transferring (FD-SWIPT) systems. To guarantee high security and energy harvesting performance of the FD-SWIPT system, the proposed design is formulated as a sum information transmission rate (SITR) maximization problem under information-leakage and energy constraints. In addition, we consider the fairness issue between the uplink and downlink information transmission rates by formulating a \mbox{fairness-aware} SITR-maximization problem. Although the formulated \mbox{SITR-maximization} and \mbox{fairness-aware} \mbox{SITR-maximization} problems are non-convex, we solve them via semidefinite relaxation and one-dimensional search. The optimality of our proposed algorithms is theoretically proved, and the computation complexities are established. Moreover, we propose two suboptimal solutions to the formulated optimization problems. In terms of the SITR-maximization problem, numerical results show that the performance achieved by one of the two suboptimal algorithms is close to the performance of the optimal algorithm with increasing maximum transmission power of the FD-BST.


title: On the multiplicity of isolated roots of sparse polynomial systems

abstract: We give formulas for the multiplicity of any affine isolated zero of a generic polynomial system of n equations in n unknowns with prescribed sets of monomials. First, we consider sets of supports such that the origin is an isolated root of the corresponding generic system and prove formulas for its multiplicity. Then, we apply these formulas to solve the problem in the general case, by showing that the multiplicity of an arbitrary affine isolated zero of a generic system with given supports equals the multiplicity of the origin as a common zero of a generic system with an associated family of supports. The formulas obtained are in the spirit of the classical Bernstein's theorem, in the sense that they depend on the combinatorial structure of the system, namely, geometric numerical invariants associated to the supports, such as mixed volumes of convex sets and, alternatively, mixed integrals of convex functions.


title: Dynamic Cross-Layer Beamforming in Hybrid Powered Communication Systems With Harvest-Use-Trade Strategy

abstract: The application of renewable energy is a promising solution to realize the Green Communications. However, if the cellular systems are solely powered by the renewable energy, the weather dependence of the renewable energy arrival makes the systems unstable. On the other hand, the proliferation of the smart grid facilitates the loads with two-way energy trading capability. Hence, a hybrid powered cellular system, which combines the smart grid with the base stations, can reduce the grid energy expenditure and improve the utilization efficiency of the renewable energy. In this paper, the long-term grid energy expenditure minimization problem is formulated as a stochastic optimization model. By leveraging the stochastic optimization theory, we reformulate the stochastic optimization problem as a \mbox{per-frame} grid energy plus weighted penalized packet rate minimization problem, which is NP-hard. As a result, two suboptimal algorithms, which jointly consider the effects of the channel quality and the packet reception failure, are proposed based on the successive approximation beamforming (SABF) technique and the \mbox{zero-forcing} beamforming (ZFBF) technique. The convergence properties of the proposed suboptimal algorithms are established, and the corresponding computational complexities are analyzed. Simulation results show that the proposed SABF algorithm outperforms the ZFBF algorithm in both grid energy expenditure and packet delay. By tuning a control parameter, the grid energy expenditure can be traded for the packet delay under the proposed stochastic optimization model.


title: Ultracompact amplitude modulator by coupling hyperbolic polaritons over a graphene-covered gap

abstract: The hyperbolic phonon-polaritons within the Reststrahlen band of hBN are of great interest for applications in nanophotonics as they are capable of propagating light signals with low losses over large distances. However, due to the phononic nature of the polaritons in hBN, amplitude modulation of its signal proves to be difficult and has been underexplored. In this paper, we propose theoretically a broadband efficient amplitude modulator for hyperbolic rays in hBN operating in the frequency range between 1450 cm$^{-1}$ and 1550 cm$^{-1}$. The modulating region comprises a few tens of nanometers wide gap carved within the hBN slab and covered by a graphene layer, where electrostatically gated graphene serves as a mediator that facilitates the coupling between phonon-polaritons on each side of the gap through plasmonic modes within graphene. We demonstrate that such an ultra compact modulator has insertion losses as low as 3 dB and provides modulation depth varying between 14 and 20 dB within the type-II hyperbolicity region of hBN.


title: On the inertial effects of density variation in stratified shear flows

abstract: In this paper, we first revisit the celebrated Boussinesq approximation in stratified flows. Using scaling arguments we show that when the background shear is weak, the Boussinesq approximation yields either (i) $A_t\ll \mathcal{O}(1)$ or (ii) $Fr_c^2 \ll \mathcal{O}(1)$, where $A_t$ is the ratio of density variation to the mean density and $Fr_c$ is the ratio of the phase speed to the long wave speed. The second clause implies, contrary to the commonly accepted notion, that a flow with large density variations can also be Boussinesq. Indeed, we show that deepwater surface gravity waves are Boussinesq while shallow water surface gravity waves are not. However, in the presence of moderate/strong shear, Boussinesq approximation implies the conventionally accepted $A_t\ll \mathcal{O}(1)$. To understand the inertial effects of density variation, our second objective is to explore various non-Boussinesq shear flows and study different kinds of stably propagating waves that can be present at an interface between two fluids of different background densities and vorticities. Furthermore, three kinds of density interfaces - neutral, stable and unstable - embedded in a background shear layer, are investigated. Instabilities ensuing from these configurations, which includes Kelvin-Helmholtz, Holmboe, Rayleigh-Taylor and triangular-jet, are studied in terms of resonant wave interactions. The effects of density stratification and the shear on the stability of each of these flow configurations are explored. Some of the results, e.g. the destabilizing role of density stratification, stabilizing role of shear, etc. are apparently counter-intuitive, but physical explanations are possible if the instabilities are interpreted from wave interactions perspective.


title: Homogeneous rank one perturbations and inverse square potentials

abstract: Following [D,BDG,DR], I describe several exactly solvable families of closed operators. Some of these families are defined by the theory of singular rank one perturbations. The remaining families are Schrodinger operators with inverse square potentials and various boundary conditions. I describe a close relationship between these families. In all of them one can observe interesting renormalization group flows (action of the group of dilations).


title: Composite quasianalytic functions

abstract: We prove two main results on Denjoy-Carleman classes: (1) a composite function theorem which asserts that a function f(x) in a quasianalytic Denjoy-Carleman class Q, which is formally composite with a generically submersive mapping y=h(x) of class Q, at a single given point in the source (or in the target) of h, can be written locally as f(x) = g(h(x)), where g(y) belongs to a shifted Denjoy-Carleman class Q' ; (2) a statement on a similar loss of regularity for functions definable in the o-minimal structure given by expansion of the real field by restricted functions of quasianalytic class Q. Both results depend on an estimate for the regularity of an infinitely differentiable solution g of the equation f(x) = g(h(x)), with f and h as above. The composite function result depends also on a quasianalytic continuation theorem, which shows that the formal assumption at a given point in (1) propagates to a formal composition condition at every point in a neighbourhood.


title: On Lagrange polynomials and the rate of approximation of planar sets by polynomial Julia sets

abstract: We revisit the approximation of nonempty compact planar sets by filled-in Julia sets of polynomials developed by Lindsey and Younsi and analyze the rate of approximation. We use slightly modified fundamental Lagrange interpolation polynomials and show that taking certain classes of nodes with subexponential growth of Lebesgue constants improves the approximation rate. To this end we investigate properties of some arrays of points in $\mathbb{C}$. In particular we prove subexponential growth of Lebesgue constants for pseudo Leja sequences with bounded Edrei growth on finite unions of quasiconformal arcs. Finally, for some classes of sets we estimate more precisely the rate of approximation by filled-in Julia sets in Hausdorff and Klimek metrics.


title: Comparing the fractal basins of attraction in the Hill problem with oblateness and radiation

abstract: The basins of convergence, associated with the roots (attractors) of a complex equation, are revealed in the Hill problem with oblateness and radiation, using a large variety of numerical methods. Three cases are investigated, regarding the values of the oblateness and radiation. In all cases, a systematic and thorough scan of the complex plane is performed in order to determine the basins of attraction of the several iterative schemes. The correlations between the attracting domains and the corresponding required number of iterations are also illustrated and discussed. Our numerical analysis strongly suggests that the basins of convergence, with the highly fractal basin boundaries, produce extraordinary and beautiful formations on the complex plane.


title: Biomechanical analysis of a cranial Patient Specific Implant on the interface with the bone using the Finite Element Method

abstract: - New advance technologies based on reverse engineering , design and additive manufacturing, have expanded design capabilities for biomedical applications to include Patient Specific Implants (PSI). This change in design paradigms needs advanced tools to assess the mechanical performance of the product, and simulate the impact on the patient. In this work, we perform a structural analysis on the interface of a cranial PSI under static loading conditions. Based on those simulations, we have identified the regions with high stress and strain and checked the failure criteria both in the implant and the skull. We evaluate the quality of the design of the implant and determine their response given different materials, in order to ensure optimality of the final product to be manufactured .


title: Multilevel mixed effects parametric survival analysis: Estimation, simulation and application

abstract: In this article, I present the user written stmixed command for the fitting of multilevel survival models, which serves as both an alternative to Stata's official mestreg, and a complimentary program with substantial extensions. stmixed can fit multilevel survival models with any number of levels and random effects at each level, including flexible spline-based approaches (such as Royston-Parmar and the log hazard equivalent) or user-defined hazard models. Simple or complex time-dependent effects can be included, as well as the addition of expected mortality for a relative survival model. Left-truncation/delayed entry can be used and t-distributed random effects are provided as an alternative to Gaussian random effects. The methods are illustrated with a commonly used dataset of patients with kidney disease suffering recurrent infections, and a simulated example, illustrating a simple approach to simulating clustered survival data using survsim (Crowther and Lambert 2012, 2013). stmixed is part of the merlin family (Crowther 2017, 2018).


title: Neutron star mergers as a probe of modifications of general relativity with finite-range scalar forces

abstract: Observations of gravitational radiation from compact binary systems provide an unprecedented opportunity to test General Relativity in the strong field dynamical regime. In this paper, we investigate how future observations of gravitational radiation from binary neutron star mergers might provide constraints on finite-range forces from a universally coupled massive scalar field. Such scalar degrees of freedom are a characteristic feature of many extensions of General Relativity. For concreteness, we work in the context of metric $f(R)$ gravity, which is equivalent to General Relativity and a universally coupled scalar field with a non-linear potential whose form is fixed by the choice of $f(R)$. In theories where neutron stars (or other compact objects) obtain a significant scalar charge, the resulting attractive finite-range scalar force has implications for both the inspiral and merger phases of binary systems. We first present an analysis of the inspiral dynamics in Newtonian limit, and forecast the constraints on the mass of the scalar and charge of the compact objects for the Advanced LIGO gravitational wave observatory. We then perform a comparative study of binary neutron star mergers in General Relativity with those of a one-parameter model of $f(R)$ gravity using fully relativistic hydrodynamical simulations. These simulations elucidate the effects of the scalar on the merger and post-merger dynamics. We comment on the utility of the full waveform (inspiral, merger, post-merger) to probe different regions of parameter space for both the particular model of $f(R)$ gravity studied here and for finite-range scalar forces more generally.


title: Estimating model evidence using ensemble-based data assimilation with localization - The model selection problem

abstract: IIn recent years, there has been a growing interest in applying data assimilation (DA) methods, originally designed for state estimation, to the model selection problem. In this setting, Carrassi et al. (2017) introduced the contextual formulation of model evidence (CME) and showed that CME can be efficiently computed using a hierarchy of ensemble-based DA procedures. Although Carrassi et al. (2017) analyzed the DA methods most commonly used for operational atmospheric and oceanic prediction worldwide, they did not study these methods in conjunction with localization to a specific domain. Yet any application of ensemble DA methods to realistic geophysical models requires the implementation of some form of localization. The present study extends the theory for estimating CME to ensemble DA methods with domain localization. The domain-localized CME (DL-CME) developed herein is tested for model selection with two models: (i) the Lorenz 40-variable mid-latitude atmospheric dynamics model (L95); and (ii) the simplified global atmospheric SPEEDY model. The CME is compared to the root-mean-square-error (RMSE) as a metric for model selection. The experiments show that CME improves systematically over the RMSE, and that this skill improvement is further enhanced by applying localization in the estimate of the CME, using the DL-CME. The potential use and range of applications of the CME and DL-CME as a model selection metric are also discussed.


title: An Attention-based Collaboration Framework for Multi-View Network Representation Learning

abstract: Learning distributed node representations in networks has been attracting increasing attention recently due to its effectiveness in a variety of applications. Existing approaches usually study networks with a single type of proximity between nodes, which defines a single view of a network. However, in reality there usually exists multiple types of proximities between nodes, yielding networks with multiple views. This paper studies learning node representations for networks with multiple views, which aims to infer robust node representations across different views. We propose a multi-view representation learning approach, which promotes the collaboration of different views and lets them vote for the robust representations. During the voting process, an attention mechanism is introduced, which enables each node to focus on the most informative views. Experimental results on real-world networks show that the proposed approach outperforms existing state-of-the-art approaches for network representation learning with a single view and other competitive approaches with multiple views.


title: Structure of free semigroupoid algebras

abstract: A free semigroupoid algebra is the closure of the algebra generated by a TCK family of a graph in the weak operator topology. We obtain a structure theory for these algebras analogous to that of free semigroup algebra. We clarify the role of absolute continuity and wandering vectors. These results are applied to obtain a Lebesgue-von Neumann-Wold decomposition of TCK families, along with reflexivity, a Kaplansky density theorem and classification for free semigroupoid algebras. Several classes of examples are discussed and developed, including self-adjoint examples and a classification of atomic free semigroupoid algebras up to unitary equivalence.


title: Time-reversal symmetry breaking phase and gapped surface states in d-wave nanoscale superconductors

abstract: We solve the Bogoliubov-de Gennes equations self-consistently for d-wave superconductors with [110] surfaces. We find spontaneous breaking of time-reversal symmetry (TRS) caused by the spontaneous occurrence of a complex order parameter with an extended s-wave symmetry along the [110] surfaces. In the TRS-breaking phase, the d-wave order parameter itself becomes complex. We also show that vortex-antivortex pairs can appear along the [110] surfaces of nanoribbons. These pairs are detectable with surface sensitive probes.


title: Different interpretations of quantum mechanics make different predictions in non-linear quantum mechanics, and some do not violate the no-signaling condition

abstract: Nonlinear modifications of quantum mechanics have a troubled history. They were initially studied for many promising reasons: resolving the measurement problem, testing the limits of standard quantum mechanics, and reconciling it with gravity. Two results substantially undermined the credibility of non-linear theories. Some have been experimentally refuted, and more importantly, all deterministic non-linear theories can be used for superluminal communication. However, these results are unconvincing because they overlook the fact that the distribution of measurement results predicted by non-linear quantum mechanics depends on the interpretation of quantum mechanics that one uses. For instance, although the Everett and Copenhagen interpretations agree on the expression of Born's rule for the outcomes of multiple measurements in linear quantum mechanics, they disagree in non-linear quantum mechanics. We present the range of expressions of Born's rule that can be obtained by applying different formulations of quantum mechanics to a class of non-linear quantum theories. We then determine that many do not allow for superluminal communication but only two seem to have a reasonable justification. The first is the Everett interpretation, and the second, which we name causal-conditional, states that a measurement broadcasts its outcome to degrees of freedom in its future light-cone, who update the wavefunction that their non-linear Hamiltonian depends on according to this new information.


title: Construction C*: an inter-level coded version of Construction C

abstract: Besides all the attention given to lattice constructions, it is common to find some very interesting nonlattice constellations, as Construction C, for example, which also has relevant applications in communication problems (multi-level coding, multi-stage decoding, good quantization efficieny). In this work we present a constellation which is a subset of Construction C, based on inter-level coding, which we call Construction C*. This construction may have better immunity to noise and it also provides a simple way of describing the Leech lattice $\Lambda_{24}.$ A condition under which Construction C* is a lattice constellation is given.


title: Dead Alphas as Risk Factors

abstract: We give an explicit algorithm and source code for extracting equity risk factors from dead (a.k.a. "flatlined" or "hockey-stick") alphas and using them to improve performance characteristics of good (tradable) alphas. In a nutshell, we use dead alphas to extract directions in the space of stock returns along which there is no money to be made (and/or those bets are too volatile). In practice the number of dead alphas can be large compared with the number of underlying stocks and care is required in identifying the aforesaid directions.


title: Entropy of random symbolic high-order bilinear Markov chains

abstract: The main goal of this paper is to develop an estimate for the entropy of random stationary ergodic symbolic sequences with elements belonging to a finite alphabet. We present here the detailed analytical study of the entropy for the high-order Markov chain in the bilinear approximation. The appendix contains a short comprehensive introduction into the subject of study.


title: Radial oscillations of strange quark stars admixed with condensed dark matter

abstract: We compute the 20 lowest frequency radial oscillation modes of strange stars admixed with condensed dark matter. We assume a self-interacting bosonic dark matter, and we model dark matter inside the star as a Bose-Einstein condensate. In this case the equation of state is a polytropic one with index $1+1/n=2$ and a constant $K$ that is computed in terms of the mass of the dark matter particle and the scattering length. Assuming a mass and a scattering length compatible with current observational bounds for self-interacting dark matter, we have integrated numerically first the Tolman-Oppenheimer-Volkoff equations for the hydrostatic equilibrium, and then the equations for the perturbations $\xi=\Delta r/r$ and $\eta=\Delta P/P$. For a compact object with certain mass and radius we have considered here three cases, namely no dark matter at all and two different dark matter scenarios. Our results show that i) the separation between consecutive modes increases with the amount of dark matter, and ii) the effect is more pronounced for higher order modes. These effects are relevant even for a strange star made of 5\% dark matter.


title: The anisotropic non-equilibrium hydrodynamic attractor

abstract: We determine the dynamical attractors associated with anisotropic hydrodynamics (aHydro) and the DNMR equations for a 0+1d conformal system using kinetic theory in the relaxation time approximation. We compare our results to the non-equilibrium attractor obtained from exact solution of the 0+1d conformal Boltzmann equation, Navier-Stokes theory, and second-order Mueller-Israel-Stewart theory. We demonstrate that the aHydro attractor equation resums an infinite number of terms in the inverse Reynolds number. The resulting resummed aHydro attractor possesses a positive longitudinal to transverse pressure ratio and is virtually indistinguishable from the exact attractor. This suggests that kinetic theory involves not only a resummation in gradients (Knudsen number) but also a novel resummation in inverse Reynolds number. We also demonstrate that the DNMR result provides a better approximation to the exact kinetic theory attractor than Mueller-Israel-Stewart theory. Finally, we introduce a new method for obtaining approximate aHydro equations which relies solely on an expansion in inverse Reynolds number, carry out this expansion to third order, and compare these third-order results to the exact kinetic theory solution.


title: Closed-Loop Statistical Verification of Stochastic Nonlinear Systems Subject to Parametric Uncertainties

abstract: This paper proposes a statistical verification framework using Gaussian processes (GPs) for simulation-based verification of stochastic nonlinear systems with parametric uncertainties. Given a small number of stochastic simulations, the proposed framework constructs a GP regression model and predicts the system's performance over the entire set of possible uncertainties. Included in the framework is a new metric to estimate the confidence in those predictions based on the variance of the GP's cumulative distribution function. This variance-based metric forms the basis of active sampling algorithms that aim to minimize prediction error through careful selection of simulations. In three case studies, the new active sampling algorithms demonstrate up to a 35% improvement in prediction error over other approaches and are able to correctly identify regions with low prediction confidence through the variance metric.


title: The long-range spin-singlet proximity effect for the Josephson system with single-crystal ferromagnet due to its band structure features

abstract: The possible explanation of the long-range proximity effect observed in the single-crystalline cobalt nanowire sandwiched between two tungsten superconducting electrodes [Wang, M. \textit{et al}. \textit{ Nat. Phys}. \textbf{6}, 389 (2010)] is proposed. The theoretical approach is based on the features of band structure of a ferromagnet. To connect the exchange field with the momentum of quasiparticles the distinction between their effective masses in majority and minority spin bands and the Fermi surface anisotropy are taken into account. The derived Eilenberger-like equations allow to obtain the renormalized effective exchange interaction that can be completely compensated for some crystallographic direction under certain conditions. The proposed theoretical model is also compared with previous approaches.


title: Dust in the first galaxies

abstract: Using cosmological volume simulations and a custom built sub-grid model for Pop~III star formation, we examine the baseline dust extinction in the first galaxies due to Pop~III metal enrichment in the first billion years of cosmic history. We find that while the most enriched, high-density lines of sight in primordial galaxies can experience a measurable amount of extinction from Pop~III dust ($E(B-V)_{\rm max}=0.07,\ A_{\rm V,max}\approx0.28$), the average extinction is very low with $\left< E(B-V) \right> \lesssim 10^{-3}$. We derive a power-law relationship between dark matter halo mass and extinction of $E(B-V)\propto M_{\rm halo}^{0.80}$. Performing a Monte Carlo parameter study, we establish the baseline reddening of the UV spectra of dwarf galaxies at high redshift due to Pop~III enrichment only. With this method, we find $\left<\beta_{\rm UV}\right>-2.51\pm0.07$, which is both nearly halo mass and redshift independent.


title: Halving the cost of quantum addition

abstract: We improve the number of T gates needed to perform an n-bit adder from 8n + O(1) to 4n + O(1). We do so via a "temporary logical-AND" construction which uses four T gates to store the logical-AND of two qubits into an ancilla and zero T gates to later erase the ancilla. This construction is equivalent to one by Jones, except that our framing makes it clear that the technique is far more widely applicable than previously realized. Temporary logical-ANDs can be applied to integer arithmetic, modular arithmetic, rotation synthesis, the quantum Fourier transform, Shor's algorithm, Grover oracles, and many other circuits. Because T gates dominate the cost of quantum computation based on the surface code, and temporary logical-ANDs are widely applicable, this represents a significant reduction in projected costs of quantum computation. In addition to our n-bit adder, we present an n-bit controlled adder circuit with T-count of 8n + O(1), a temporary adder that can be computed for the same cost as the normal adder but whose result can be kept until it is later uncomputed without using T gates, and discuss some other constructions whose T-count is improved by the temporary logical-AND.


title: Orbital decomposition of CALIFA spiral galaxies

abstract: Schwarzschild orbit-based dynamical models are widely used to uncover the internal dynamics of early-type galaxies and globular clusters. Here we present for the first time the Schwarzschild models of late-type galaxies: an SBb galaxy NGC 4210 and an S0 galaxy NGC 6278 from the CALIFA survey. The mass profiles within $2\,R_e$ are constrained well with $1\sigma$ statistical error of $\sim 10\%$. The luminous and dark mass can be disentangled with uncertainties of $\sim 20\%$ and $\sim 50\%$ respectively. From $R_e$ to $2\,R_e$, the dark matter fraction increases from $14\pm10\%$ to $18\pm10\%$ for NGC 4210 and from $15\pm10\%$ to $30\pm20\%$ for NGC 6278. The velocity anisotropy profiles of both $\sigma_r/\sigma_t$ and $\sigma_z/\sigma_R$ are well constrained. The inferred internal orbital distributions reveal clear substructures. The orbits are naturally separated into three components: a cold component with near circular orbits, a hot component with near radial orbits, and a warm component in between. The photometrically-identified exponential disks are predominantly made up of cold orbits only beyond $\sim 1\,R_e$, while they are constructed mainly with the warm orbits inside. Our dynamical hot components are concentrated in the inner regions, similar to the photometrically-identified bulges. The reliability of the results, especially the orbit distribution, are verified by applying the model to mock data.


title: On Graphs and the Gotsman-Linial Conjecture for d = 2

abstract: We give an infinite class of counterexamples to the Gotsman-Linial conjecture when d = 2. On the other hand, we establish an asymptotic form of the conjecture for quadratic threshold functions whose non-zero quadratic terms define a graph with either low fractional chromatic number or few edges. Our techniques are elementary and our exposition is self-contained, if you're into that.


title: On the Kottwitz conjecture for local shtuka spaces

abstract: Kottwitz's conjecture describes the contribution of a supercuspidal represention to the cohomology of a local Shimura variety in terms of the local Langlands correspondence. A natural extension of this conjecture concerns Scholze's more general spaces of local shtukas. Using a new Lefschetz-Verdier trace formula for v-stacks, we prove the extended conjecture, disregarding the action of the Weil group, and modulo a virtual representation whose character vanishes on the locus of elliptic elements. As an application, we show that for an irreducible smooth representation of an inner form of $\mathrm{GL}_n$, the $L$-parameter constructed by Fargues-Scholze agrees with the usual semisimplified parameter arising from local Langlands.


title: Distributed event-triggered control for multi-agent formation stabilization and tracking

abstract: This paper addresses the problem of formation control and tracking a of desired trajectory by an Euler-Lagrange multi-agent systems. It is inspired by recent results by Qingkai et al. and adopts an event-triggered control strategy to reduce the number of communications between agents. For that purpose, to evaluate its control input, each agent maintains estimators of the states of the other agents. Communication is triggered when the discrepancy between the actual state of an agent and the corresponding estimate reaches some threshold. The impact of additive state perturbations on the formation control is studied. A condition for the convergence of the multi-agent system to a stable formation is studied. Simulations show the effectiveness of the proposed approach.


title: Unique Information via Dependency Constraints

abstract: The partial information decomposition (PID) is perhaps the leading proposal for resolving information shared between a set of sources and a target into redundant, synergistic, and unique constituents. Unfortunately, the PID framework has been hindered by a lack of a generally agreed-upon, multivariate method of quantifying the constituents. Here, we take a step toward rectifying this by developing a decomposition based on a new method that quantifies unique information. We first develop a broadly applicable method---the dependency decomposition---that delineates how statistical dependencies influence the structure of a joint distribution. The dependency decomposition then allows us to define a measure of the information about a target that can be uniquely attributed to a particular source as the least amount which the source-target statistical dependency can influence the information shared between the sources and the target. The result is the first measure that satisfies the core axioms of the PID framework while not satisfying the Blackwell relation, which depends on a particular interpretation of how the variables are related. This makes a key step forward to a practical PID.

categories: cond-mat.stat-mech cs.IT cs.LG math.IT math.ST stat.TH


title: Keeping Context In Mind: Automating Mobile App Access Control with User Interface Inspection

abstract: Recent studies observe that app foreground is the most striking component that influences the access control decisions in mobile platform, as users tend to deny permission requests lacking visible evidence. However, none of the existing permission models provides a systematic approach that can automatically answer the question: Is the resource access indicated by app foreground? In this work, we present the design, implementation, and evaluation of COSMOS, a context-aware mediation system that bridges the semantic gap between foreground interaction and background access, in order to protect system integrity and user privacy. Specifically, COSMOS learns from a large set of apps with similar functionalities and user interfaces to construct generic models that detect the outliers at runtime. It can be further customized to satisfy specific user privacy preference by continuously evolving with user decisions. Experiments show that COSMOS achieves both high precision and high recall in detecting malicious requests. We also demonstrate the effectiveness of COSMOS in capturing specific user preferences using the decisions collected from 24 users and illustrate that COSMOS can be easily deployed on smartphones as a real-time guard with a very low performance overhead.


title: Low loss QKD optical scheme for fast polarization encoding

abstract: We present a new optical scheme for BB84 protocol quantum key distribution (QKD). The proposed setup consists of a compact all-fiber polarization encoding optical scheme based on LiNbO3 phase modulators, single laser source and two single-photon detectors. Optical scheme consists of standard telecommunication components and is suitable for both fiber and free-space quantum communication channels. Low losses (~2dB) in Bob's device increase both the key generation rate and distance limit. A new technique for solving polarization mode dispersion (PMD) issue in LiNbO3 is implemented, allowing two crystals to neutralize the effect of each other. Several proof-of-concept experiments have been conducted at 10 MHz repetition frequency over 50 km of standard optical fiber under laboratory conditions and over 30 km of urban fiber with high losses (13dB), which is a link within a QKD network. To achieve this, calibration algorithms have been developed, allowing the system to work autonomously and making it promising for practical applications.


title: Deep Reinforcement Learning for Event-Driven Multi-Agent Decision Processes

abstract: The incorporation of macro-actions (temporally extended actions) into multi-agent decision problems has the potential to address the curse of dimensionality associated with such decision problems. Since macro-actions last for stochastic durations, multiple agents executing decentralized policies in cooperative environments must act asynchronously. We present an algorithm that modifies generalized advantage estimation for temporally extended actions, allowing a state-of-the-art policy optimization algorithm to optimize policies in Dec-POMDPs in which agents act asynchronously. We show that our algorithm is capable of learning optimal policies in two cooperative domains, one involving real-time bus holding control and one involving wildfire fighting with unmanned aircraft. Our algorithm works by framing problems as "event-driven decision processes," which are scenarios in which the sequence and timing of actions and events are random and governed by an underlying stochastic process. In addition to optimizing policies with continuous state and action spaces, our algorithm also facilitates the use of event-driven simulators, which do not require time to be discretized into time-steps. We demonstrate the benefit of using event-driven simulation in the context of multiple agents taking asynchronous actions. We show that fixed time-step simulation risks obfuscating the sequence in which closely separated events occur, adversely affecting the policies learned. In addition, we show that arbitrarily shrinking the time-step scales poorly with the number of agents.


title: Particle number dependence in the non-linear evolution of N-body self-gravitating systems

abstract: Simulations of purely self-gravitating N-body systems are often used in astrophysics and cosmology to study the collisionless limit of such systems. Their results for macroscopic quantities should then converge well for sufficiently large N. Using a study of the evolution from a simple space of spherical initial conditions - including a region characterised by so-called "radial orbit instability" - we illustrate that the values of N at which such convergence is obtained can vary enormously. In the family of initial conditions we study, good convergence can be obtained up to a few dynamical times with N $ \sim 10^3$ - just large enough to suppress two body relaxation - for certain initial conditions, while in other cases such convergence is not attained at this time even in our largest simulations with N $\sim 10^5$. The qualitative difference is due to the stability properties of fluctuations introduced by the N-body discretisation, of which the initial amplitude depends on N. We discuss briefly why the crucial role which such fluctuations can potentially play in the evolution of the N-body system could, in particular, constitute a serious problem in cosmological simulations of dark matter.


title: An Introduction to Rule-based Modeling of Immune Receptor Signaling

abstract: Cells process external and internal signals through chemical interactions. Cells that constitute the immune system (e.g., antigen presenting cell, T-cell, B-cell, mast cell) can have different functions (e.g., adaptive memory, inflammatory response) depending on the type and number of receptor molecules on the cell surface and the specific intracellular signaling pathways activated by those receptors. Explicitly modeling and simulating kinetic interactions between molecules allows us to pose questions about the dynamics of a signaling network under various conditions. However, the application of chemical kinetics to biochemical signaling systems has been limited by the complexity of the systems under consideration. Rule-based modeling (BioNetGen, Kappa, Simmune, PySB) is an approach to address this complexity. In this chapter, by application to the Fc$\varepsilon$RI receptor system, we will explore the origins of complexity in macromolecular interactions, show how rule-based modeling can be used to address complexity, and demonstrate how to build a model in the BioNetGen framework. Open source BioNetGen software and documentation are available at http://bionetgen.org.


title: Benchmarking Numerical Methods for Lattice Equations with the Toda Lattice

abstract: We compare performances of well-known numerical time-stepping methods that are widely used to compute solutions of the doubly-infinite Fermi-Pasta-Ulam-Tsingou (FPUT) lattice equations. The methods are benchmarked according to (1) their accuracy in capturing the soliton peaks and (2) in capturing highly-oscillatory parts of the solutions of the Toda lattice resulting from a variety of initial data. The numerical inverse scattering transform method is used to compute a reference solution with high accuracy. We find that benchmarking a numerical method on pure-soliton initial data can lead one to overestimate the accuracy of the method.


title: Space Technology for Directly Imaging and Characterizing Exo-Earths

abstract: The detection of Earth-like exoplanets in the habitable zone of their stars, and their spectroscopic characterization in a search for biosignatures, requires starlight suppression that exceeds the current best ground-based performance by orders of magnitude. The required planet/star brightness ratio of order 1e-10 at visible wavelengths can be obtained by blocking stellar photons with an occulter, either externally (a starshade) or internally (a coronagraph) to the telescope system, and managing diffracted starlight, so as to directly image the exoplanet in reected starlight. Coronagraph instruments require advancement in telescope aperture (either monolithic or segmented), aperture obscurations (obscured by secondary mirror and its support struts), and wavefront error sensitivity (e.g. line-of-sight jitter, telescope vibration, polarization). The starshade, which has never been used in a science application, benefits a mission by being decoupled from the telescope, allowing a loosening of telescope stability requirements. In doing so, it transfers the difficult technology from the telescope system to a large deployable structure (tens of meters to greater than 100 m in diameter) that must be positioned precisely at a distance of tens of thousands of kilometers from the telescope. We describe in this paper a roadmap to achieving the technological capability to search for biosignatures on an Earth-like exoplanet from a future space telescope. Two of these studies, HabEx and LUVOIR, include the direct imaging of Earth-sized habitable exoplanets as a central science theme.


title: Simulation-based reachability analysis for nonlinear systems using componentwise contraction properties

abstract: A shortcoming of existing reachability approaches for nonlinear systems is the poor scalability with the number of continuous state variables. To mitigate this problem we present a simulation-based approach where we first sample a number of trajectories of the system and next establish bounds on the convergence or divergence between the samples and neighboring trajectories. We compute these bounds using contraction theory and reduce the conservatism by partitioning the state vector into several components and analyzing contraction properties separately in each direction. Among other benefits this allows us to analyze the effect of constant but uncertain parameters by treating them as state variables and partitioning them into a separate direction. We next present a numerical procedure to search for weighted norms that yield a prescribed contraction rate, which can be incorporated in the reachability algorithm to adjust the weights to minimize the growth of the reachable set.


title: Verifying Properties of Binarized Deep Neural Networks

abstract: Understanding properties of deep neural networks is an important challenge in deep learning. In this paper, we take a step in this direction by proposing a rigorous way of verifying properties of a popular class of neural networks, Binarized Neural Networks, using the well-developed means of Boolean satisfiability. Our main contribution is a construction that creates a representation of a binarized neural network as a Boolean formula. Our encoding is the first exact Boolean representation of a deep neural network. Using this encoding, we leverage the power of modern SAT solvers along with a proposed counterexample-guided search procedure to verify various properties of these networks. A particular focus will be on the critical property of robustness to adversarial perturbations. For this property, our experimental results demonstrate that our approach scales to medium-size deep neural networks used in image classification tasks. To the best of our knowledge, this is the first work on verifying properties of deep neural networks using an exact Boolean encoding of the network.


title: Linear Computer-Music through Sequences over Galois Fields

abstract: It is shown how binary sequences can be associated with automatic composition of monophonic pieces. We are concerned with the composition of e-music from finite field structures. The information at the input may be either random or information from a black-and-white, grayscale or color picture. New e-compositions and music score are made available, including a new piece from the famous Lenna picture: the score of the e-music <<Between Lenna's eyes in C major.>> The corresponding stretch of music score are presented. Some particular structures, including clock arithmetic (mod 12), GF(7), GF(8), GF(13) and GF(17) are addressed. Further, multilevel block-codes are also used in a new approach of e-music composition, engendering a particular style as an e-composer. As an example, Pascal multilevel block codes recently introduced are handled to generate a new style of electronic music over GF(13).


title: Curriculum Learning of Visual Attribute Clusters for Multi-Task Classification

abstract: Visual attributes, from simple objects (e.g., backpacks, hats) to soft-biometrics (e.g., gender, height, clothing) have proven to be a powerful representational approach for many applications such as image description and human identification. In this paper, we introduce a novel method to combine the advantages of both multi-task and curriculum learning in a visual attribute classification framework. Individual tasks are grouped after performing hierarchical clustering based on their correlation. The clusters of tasks are learned in a curriculum learning setup by transferring knowledge between clusters. The learning process within each cluster is performed in a multi-task classification setup. By leveraging the acquired knowledge, we speed-up the process and improve performance. We demonstrate the effectiveness of our method via ablation studies and a detailed analysis of the covariates, on a variety of publicly available datasets of humans standing with their full-body visible. Extensive experimentation has proven that the proposed approach boosts the performance by 4% to 10%.


title: Inverse mean curvature evolution of entire graphs

abstract: We study the evolution of strictly mean-convex entire graphs over $R^n$ by Inverse Mean Curvature flow. First we establish the global existence of starshaped entire graphs with superlinear growth at infinity. The main result in this work concerns the critical case of asymptotically conical entire convex graphs. In this case we show that there exists a time $ T < +\infty$, which depends on the growth at infinity of the initial data, such that the unique solution of the flow exists for all $t < T$. Moreover, as $t \to T$ the solution converges to a flat plane. Our techniques exploit the ultra-fast diffusion character of the fully-nonlinear flow, a property that implies that the asymptotic behavior at spatial infinity of our solution plays a crucial influence on the maximal time of existence, as such behavior propagates infinitely fast towards the interior.


title: Colored Khovanov-Rozansky homology for infinite braids

abstract: We show that the limiting unicolored $\mathfrak{sl}(N)$ Khovanov-Rozansky chain complex of any infinite positive braid categorifies a highest-weight projector. This result extends an earlier result of Cautis categorifying highest-weight projectors using the limiting complex of infinite torus braids. Additionally, we show that the results hold in the case of colored HOMFLY-PT Khovanov-Rozansky homology as well. An application of this result is given in finding a partial isomorphism between the HOMFLY-PT homology of any braid positive link and the stable HOMFLY-PT homology of the infinite torus knot as computed by Hogancamp.


title: Phosphorous dimerization in GaP at high pressure

abstract: Using combined experimental and computational approaches we show that at 43 GPa and 1300 K gallium phosphide adopts the super-Cmcm structure, here indicated with its Pearson notation oS24. First-principles enthalpy calculations demonstrate that this oS24 structure is more thermodynamically stable above about 20 GPa than other proposed polymorphs, such as cubic ZB and SC16 or the ubiquitous orthorhombic Cmcm phase (oS8). It is found that oS24-structured GaP exhibits short P-P bonds, according to our high-resolution structural analysis. Such phosphorous dimerization in GaP, observed for the first time and confirmed by ab initio calculations, sheds light on the nature of the super-Cmcm structure and provides critical new insights into the high-pressure polymorphism of octet semiconductors. Bond directionality and anisotropy, in addition to 5-fold phosphorous coordination, explain the relatively low symmetry of this high-pressure structure.


title: Fast and Reliable Autonomous Surgical Debridement with Cable-Driven Robots Using a Two-Phase Calibration Procedure

abstract: Automating precision subtasks such as debridement (removing dead or diseased tissue fragments) with Robotic Surgical Assistants (RSAs) such as the da Vinci Research Kit (dVRK) is challenging due to inherent non-linearities in cable-driven systems. We propose and evaluate a novel two-phase coarse-to-fine calibration method. In Phase I (coarse), we place a red calibration marker on the end effector and let it randomly move through a set of open-loop trajectories to obtain a large sample set of camera pixels and internal robot end-effector configurations. This coarse data is then used to train a Deep Neural Network (DNN) to learn the coarse transformation bias. In Phase II (fine), the bias from Phase I is applied to move the end-effector toward a small set of specific target points on a printed sheet. For each target, a human operator manually adjusts the end-effector position by direct contact (not through teleoperation) and the residual compensation bias is recorded. This fine data is then used to train a Random Forest (RF) to learn the fine transformation bias. Subsequent experiments suggest that without calibration, position errors average 4.55mm. Phase I can reduce average error to 2.14mm and the combination of Phase I and Phase II can reduces average error to 1.08mm. We apply these results to debridement of raisins and pumpkin seeds as fragment phantoms. Using an endoscopic stereo camera with standard edge detection, experiments with 120 trials achieved average success rates of 94.5%, exceeding prior results with much larger fragments (89.4%) and achieving a speedup of 2.1x, decreasing time per fragment from 15.8 seconds to 7.3 seconds. Source code, data, and videos are available at https://sites.google.com/view/calib-icra/.


title: A textual transform of multivariate time-series for prognostics

abstract: Prognostics or early detection of incipient faults is an important industrial challenge for condition-based and preventive maintenance. Physics-based approaches to modeling fault progression are infeasible due to multiple interacting components, uncontrolled environmental factors and observability constraints. Moreover, such approaches to prognostics do not generalize to new domains. Consequently, domain-agnostic data-driven machine learning approaches to prognostics are desirable. Damage progression is a path-dependent process and explicitly modeling the temporal patterns is critical for accurate estimation of both the current damage state and its progression leading to total failure. In this paper, we present a novel data-driven approach to prognostics that employs a novel textual representation of multivariate temporal sensor observations for predicting the future health state of the monitored equipment early in its life. This representation enables us to utilize well-understood concepts from text-mining for modeling, prediction and understanding distress patterns in a domain agnostic way. The approach has been deployed and successfully tested on large scale multivariate time-series data from commercial aircraft engines. We report experiments on well-known publicly available benchmark datasets and simulation datasets. The proposed approach is shown to be superior in terms of prediction accuracy, lead time to prediction and interpretability.


title: Dex-Net 3.0: Computing Robust Robot Vacuum Suction Grasp Targets in Point Clouds using a New Analytic Model and Deep Learning

abstract: Vacuum-based end effectors are widely used in industry and are often preferred over parallel-jaw and multifinger grippers due to their ability to lift objects with a single point of contact. Suction grasp planners often target planar surfaces on point clouds near the estimated centroid of an object. In this paper, we propose a compliant suction contact model that computes the quality of the seal between the suction cup and local target surface and a measure of the ability of the suction grasp to resist an external gravity wrench. To characterize grasps, we estimate robustness to perturbations in end-effector and object pose, material properties, and external wrenches. We analyze grasps across 1,500 3D object models to generate Dex-Net 3.0, a dataset of 2.8 million point clouds, suction grasps, and grasp robustness labels. We use Dex-Net 3.0 to train a Grasp Quality Convolutional Neural Network (GQ-CNN) to classify robust suction targets in point clouds containing a single object. We evaluate the resulting system in 350 physical trials on an ABB YuMi fitted with a pneumatic suction gripper. When evaluated on novel objects that we categorize as Basic (prismatic or cylindrical), Typical (more complex geometry), and Adversarial (with few available suction-grasp points) Dex-Net 3.0 achieves success rates of 98$\%$, 82$\%$, and 58$\%$ respectively, improving to 81$\%$ in the latter case when the training set includes only adversarial objects. Code, datasets, and supplemental material can be found at http://berkeleyautomation.github.io/dex-net .


title: Think Globally, Embed Locally --- Locally Linear Meta-embedding of Words

abstract: Distributed word embeddings have shown superior performances in numerous Natural Language Processing (NLP) tasks. However, their performances vary significantly across different tasks, implying that the word embeddings learnt by those methods capture complementary aspects of lexical semantics. Therefore, we believe that it is important to combine the existing word embeddings to produce more accurate and complete \emph{meta-embeddings} of words. For this purpose, we propose an unsupervised locally linear meta-embedding learning method that takes pre-trained word embeddings as the input, and produces more accurate meta embeddings. Unlike previously proposed meta-embedding learning methods that learn a global projection over all words in a vocabulary, our proposed method is sensitive to the differences in local neighbourhoods of the individual source word embeddings. Moreover, we show that vector concatenation, a previously proposed highly competitive baseline approach for integrating word embeddings, can be derived as a special case of the proposed method. Experimental results on semantic similarity, word analogy, relation classification, and short-text classification tasks show that our meta-embeddings to significantly outperform prior methods in several benchmark datasets, establishing a new state of the art for meta-embeddings.


title: Reasoning about proof and knowledge

abstract: In previous work [Lewitzka, Log. J. IGPL 2017], we presented a hierarchy of classical modal systems, along with algebraic semantics, for the reasoning about intuitionistic truth, belief and knowledge. Deviating from G\"odel's interpretation of IPC in S4, our modal systems contain IPC in the way established in [Lewitzka, J. Log. Comp. 2015]. The modal operator can be viewed as a predicate for intuitionistic truth, i.e. proof. Epistemic principles are partially adopted from Intuitionistic Epistemic Logic IEL [Artemov and Protopopescu, Rev. Symb. Log. 2016]. In the present paper, we show that the S5-style systems of our hierarchy correspond to an extended Brouwer-Heyting-Kolmogorov interpretation and are complete w.r.t. a relational semantics based on intuitionistic general frames. In this sense, our S5-style logics are adequate and complete systems for the reasoning about proof combined with belief or knowledge. The proposed relational semantics is a uniform framework in which also IEL can be modeled. Verification-based intuitionistic knowledge formalized in IEL turns out to be a special case of the kind of knowledge described by our S5-style systems.


title: Why PairDiff works? -- A Mathematical Analysis of Bilinear Relational Compositional Operators for Analogy Detection

abstract: Representing the semantic relations that exist between two given words (or entities) is an important first step in a wide-range of NLP applications such as analogical reasoning, knowledge base completion and relational information retrieval. A simple, yet surprisingly accurate method for representing a relation between two words is to compute the vector offset (\PairDiff) between their corresponding word embeddings. Despite the empirical success, it remains unclear as to whether \PairDiff is the best operator for obtaining a relational representation from word embeddings. We conduct a theoretical analysis of generalised bilinear operators that can be used to measure the $\ell_{2}$ relational distance between two word-pairs. We show that, if the word embeddings are standardised and uncorrelated, such an operator will be independent of bilinear terms, and can be simplified to a linear form, where \PairDiff is a special case. For numerous word embedding types, we empirically verify the uncorrelation assumption, demonstrating the general applicability of our theoretical result. Moreover, we experimentally discover \PairDiff from the bilinear relation composition operator on several benchmark analogy datasets.


title: Microscopic Dynamics of Li$^{+}$ in Rutile TiO$_{2}$ Revealed by $^{8}$Li $\beta$-detected NMR

abstract: We report measurements of the dynamics of isolated $^{8}$Li$^{+}$ in single crystal rutile TiO$_{2}$ using $\beta$-detected NMR. From spin-lattice relaxation and motional narrowing, we find two sets of thermally activated dynamics: one below 100 K; and one at higher temperatures. At low temperature, the activation barrier is $26.8(6)$ meV with prefactor $1.23(5) \times 10^{10}$ s$^{-1}$. We suggest this is unrelated to Li$^{+}$ motion, and rather is a consequence of electron polarons in the vicinity of the implanted $^{8}$Li$^{+}$ that are known to become mobile in this temperature range. Above 100 K, Li$^{+}$ undergoes long-range diffusion as an isolated uncomplexed cation, characterized by an activation energy and prefactor of $0.32(2)$ eV and $1.0(5) \times 10^{16}$ s$^{-1}$, in agreement with macroscopic diffusion measurements. These results in the dilute limit from a microscopic probe indicate that Li$^{+}$ concentration does not limit the diffusivity even up to high concentrations, but that some key ingredient is missing in the calculations of the migration barrier. The anomalous prefactors provide further insight into both Li$^{+}$ and polaron motion.


title: Talk Resource-Efficiently to Me: Optimal Communication Planning for Distributed Loop Closure Detection

abstract: Due to the distributed nature of cooperative simultaneous localization and mapping (CSLAM), detecting inter-robot loop closures necessitates sharing sensory data with other robots. A na\"{\i}ve approach to data sharing can easily lead to a waste of mission-critical resources. This paper investigates the logistical aspects of CSLAM. Particularly, we present a general resource-efficient communication planning framework that takes into account both the total amount of exchanged data and the induced division of labor between the participating robots. Compared to other state-of-the-art approaches, our framework is able to verify the same set of potential inter-robot loop closures while exchanging considerably less data and influencing the induced workloads. We develop a fast algorithm for finding globally optimal communication policies, and present theoretical analysis to characterize the necessary and sufficient conditions under which simpler strategies are optimal. The proposed framework is extensively evaluated with data from the KITTI odometry benchmark datasets.


title: Evolution of Interfaces for the Nonlinear Double Degenerate Parabolic Equation of Turbulent Filtration with Absorption

abstract: We prove the short-time asymptotic formula for the interfaces and local solutions near the interfaces for the nonlinear double degenerate reaction-diffusion equation of turbulent filtration with strong absorption \[ u_t=\Big(|(u^{m})_x|^{p-1}(u^{m})_x\Big)_x-bu^{\beta}, \, mp>1, \, \beta >0. \] Full classification is pursued in terms of the nonlinearity parameters $m, p,\beta$ and asymptotics of the initial function near its support. Numerical analysis using a weighted essentially nonoscillatory (WENO) scheme with interface capturing is implemented, and comparison of numerical and analytical results is presented.


title: Balanced truncation for model order reduction of linear dynamical systems with quadratic outputs

abstract: We investigate model order reduction (MOR) for linear dynamical systems, where a quadratic output is defined as a quantity of interest. The system can be transformed into a linear dynamical system with many linear outputs. MOR is feasible by the method of balanced truncation, but suffers from the large number of outputs in approximate methods. To ameliorate this shortcoming we derive an equivalent quadratic-bilinear system with a single linear output and analyze the properties of this system. We examine MOR for this system via the technique of balanced truncation, which requires a stabilization of the system. Therein, the solution of two quadratic Lyapunov equations is traced back to the solution of just two linear Lyapunov equations. We present numerical results for several test examples comparing the two MOR approaches.


title: A blueprint for demonstrating quantum supremacy with superconducting qubits

abstract: Fundamental questions in chemistry and physics may never be answered due to the exponential complexity of the underlying quantum phenomena. A desire to overcome this challenge has sparked a new industry of quantum technologies with the promise that engineered quantum systems can address these hard problems. A key step towards demonstrating such a system will be performing a computation beyond the capabilities of any classical computer, achieving so-called quantum supremacy. Here, using 9 superconducting qubits, we demonstrate an immediate path towards quantum supremacy. By individually tuning the qubit parameters, we are able to generate thousands of unique Hamiltonian evolutions and probe the output probabilities. The measured probabilities obey a universal distribution, consistent with uniformly sampling the full Hilbert-space. As the number of qubits in the algorithm is varied, the system continues to explore the exponentially growing number of states. Combining these large datasets with techniques from machine learning allows us to construct a model which accurately predicts the measured probabilities. We demonstrate an application of these algorithms by systematically increasing the disorder and observing a transition from delocalized states to localized states. By extending these results to a system of 50 qubits, we hope to address scientific questions that are beyond the capabilities of any classical computer.


title: Controllability and data-driven identification of bipartite consensus on nonlinear signed networks

abstract: Nonlinear networked systems are of interest in several areas of research, such as multi-agent systems and social networks. In this paper, we examine the controllability of several classes of nonlinear networked dynamics on which the underlying graph admits negative weights. Such signed networks exhibit bipartite clustering when the underlying graph is structurally balanced. We show that structural balance is the key ingredient inducing uncontrollability when combined with a leader-node symmetry and a certain type of dynamical symmetry. We then examine the problem of extracting the bipartite structure of such graphs from data using Extended Dynamic Mode Decomposition to approximate the corresponding Koopman operator.


title: Deep Lattice Networks and Partial Monotonic Functions

abstract: We propose learning deep models that are monotonic with respect to a user-specified set of inputs by alternating layers of linear embeddings, ensembles of lattices, and calibrators (piecewise linear functions), with appropriate constraints for monotonicity, and jointly training the resulting network. We implement the layers and projections with new computational graph nodes in TensorFlow and use the ADAM optimizer and batched stochastic gradients. Experiments on benchmark and real-world datasets show that six-layer monotonic deep lattice networks achieve state-of-the art performance for classification and regression with monotonicity guarantees.


title: Environmental effects on galaxy evolution. II: quantifying the tidal features in NIR-images of the cluster Abell 85

abstract: This work is part of a series of papers devoted to investigate the evolution of cluster galaxies during their infall. In the present article we imaged in NIR a selected sample of galaxies through- out the massive cluster Abell 85 (z = 0.055). We obtained (JHK) photometry for 68 objects, reaching 1 mag/arcsec^2 deeper than 2MASS. We use these images to unveil asymmetries in the outskirts of a sample of bright galaxies and develop a new asymmetry index, alpha_An, which allows to quantify the degree of disruption by the relative area occupied by the tidal features on the plane of the sky. We measure the asymmetries for a subsample of 41 large area objects finding clear asymmetries in ten galaxies, most of them being in groups and pairs projected at different clustercentric distances, some of them located beyond R500 . Combining information on the Hi-gas content of blue galaxies and the distribution of sub-structures across Abell 85, with the present NIR asymmetry analysis, we obtain a very powerful tool to confirm that tidal mechanisms are indeed present and are currently affecting a fraction of galaxies in Abell 85. However, when comparing our deep NIR images with UV-blue images of two very disrupted (jellyfish) galaxies in this cluster, we discard the presence of tidal 1 interactions down to our detection limit. Our results suggest that ram-pressure stripping is at the origin of such spectacular disruptions. We conclude that across a complex cluster like Abell 85, environment mechanisms, both gravitational and hydrodynamical, are playing an active role in driving galaxy evolution.


title: Random matrices: overcrowding estimates for the spectrum

abstract: We address overcrowding estimates for the singular values of random iid matrices, as well as for the eigenvalues of random Wigner matrices. We show evidence of long range separation under arbitrary perturbation even in matrices of discrete entry distributions. In many cases our method yields nearly optimal bounds


title: OptionGAN: Learning Joint Reward-Policy Options using Generative Adversarial Inverse Reinforcement Learning

abstract: Reinforcement learning has shown promise in learning policies that can solve complex problems. However, manually specifying a good reward function can be difficult, especially for intricate tasks. Inverse reinforcement learning offers a useful paradigm to learn the underlying reward function directly from expert demonstrations. Yet in reality, the corpus of demonstrations may contain trajectories arising from a diverse set of underlying reward functions rather than a single one. Thus, in inverse reinforcement learning, it is useful to consider such a decomposition. The options framework in reinforcement learning is specifically designed to decompose policies in a similar light. We therefore extend the options framework and propose a method to simultaneously recover reward options in addition to policy options. We leverage adversarial methods to learn joint reward-policy options using only observed expert states. We show that this approach works well in both simple and complex continuous control tasks and shows significant performance increases in one-shot transfer learning.


title: The Classification Problem for Simple Unital Finite Rank Dimension Groups

abstract: The Borel complexity of the isomorphism problem for finite-rank unital simple dimension groups increases with rank. This implies that the isomorphism problems for the corresponding classes of Bratteli diagrams and LDA-groups also increase with rank.


title: Concentration of distances in Wigner matrices

abstract: It is well-known that distances in random iid matrices are highly concentrated around their mean. In this note we extend this concentration phenomenon to Wigner matrices. Exponential bounds for the lower tail are also included.


title: Sieve: Actionable Insights from Monitored Metrics in Microservices

abstract: Major cloud computing operators provide powerful monitoring tools to understand the current (and prior) state of the distributed systems deployed in their infrastructure. While such tools provide a detailed monitoring mechanism at scale, they also pose a significant challenge for the application developers/operators to transform the huge space of monitored metrics into useful insights. These insights are essential to build effective management tools for improving the efficiency, resiliency, and dependability of distributed systems. This paper reports on our experience with building and deploying Sieve - a platform to derive actionable insights from monitored metrics in distributed systems. Sieve builds on two core components: a metrics reduction framework, and a metrics dependency extractor. More specifically, Sieve first reduces the dimensionality of metrics by automatically filtering out unimportant metrics by observing their signal over time. Afterwards, Sieve infers metrics dependencies between distributed components of the system using a predictive-causality model by testing for Granger Causality. We implemented Sieve as a generic platform and deployed it for two microservices-based distributed systems: OpenStack and ShareLatex. Our experience shows that (1) Sieve can reduce the number of metrics by at least an order of magnitude (10 - 100$\times$), while preserving the statistical equivalence to the total number of monitored metrics; (2) Sieve can dramatically improve existing monitoring infrastructures by reducing the associated overheads over the entire system stack (CPU - 80%, storage - 90%, and network - 50%); (3) Lastly, Sieve can be effective to support a wide-range of workflows in distributed systems - we showcase two such workflows: orchestration of autoscaling, and Root Cause Analysis (RCA).


title: Phase-field modeling and electronic structural analysis of flexoelectric effect at 180{\deg} domain walls in ferroelectric PbTiO3

abstract: Flexoelectric effect is the coupling between strain, polarization and their gradients, which are prominent at the nanoscale. Although this effect is important to understand nanostructures, such as domain walls in ferroelectrics, its electronic mechanism is not clear. In this work, we combined phase-field simulations and first-principles calculations to study the 180{\deg} domain walls in tetragonal ferroelectric PbTiO3, and found that the ultimate source of N\'eel components is the gradient of the square of spontaneous polarizations. Electronic structural analysis reveals that there is a redistribution of electronic charge density and potential around domain walls, which produces the electric field and N\'eel components. This work thus sheds light on the electronic mechanism of the flexoelectric effect around 180{\deg} domain walls in tetragonal ferroelectrics.


title: Property Testing in High Dimensional Ising models

abstract: This paper explores the information-theoretic limitations of graph property testing in zero-field Ising models. Instead of learning the entire graph structure, sometimes testing a basic graph property such as connectivity, cycle presence or maximum clique size is a more relevant and attainable objective. Since property testing is more fundamental than graph recovery, any necessary conditions for property testing imply corresponding conditions for graph recovery, while custom property tests can be statistically and/or computationally more efficient than graph recovery based algorithms. Understanding the statistical complexity of property testing requires the distinction of ferromagnetic (i.e., positive interactions only) and general Ising models. Using combinatorial constructs such as graph packing and strong monotonicity, we characterize how target properties affect the corresponding minimax upper and lower bounds within the realm of ferromagnets. On the other hand, by studying the detection of an antiferromagnetic (i.e., negative interactions only) Curie-Weiss model buried in Rademacher noise, we show that property testing is strictly more challenging over general Ising models. In terms of methodological development, we propose two types of correlation based tests: computationally efficient screening for ferromagnets, and score type tests for general models, including a fast cycle presence test. Our correlation screening tests match the information-theoretic bounds for property testing in ferromagnets.


title: $n$-exangulated categories

abstract: For each positive integer $n$ we introduce the notion of $n$-exangulated categories as higher dimensional analogues of extriangulated categories defined by Nakaoka-Palu. We characterize which $n$-exangulated categories are $n$-exact in the sense of Jasso and which are $(n+2)$-angulated in the sense of Geiss-Keller-Oppermann. For extriangulated categories with enough projectives and injectives we introduce the notion of $n$-cluster tilting subcategories and show that under certain conditions such $n$-cluster tilting subcategories are $n$-exangulated.


title: Social Network Differences of Chronotypes Identified from Mobile Phone Data

abstract: Human activity follows an approximately 24-hour day-night cycle, but there is significant individual variation in awake and sleep times. Individuals with circadian rhythms at the extremes can be categorized into two chronotypes: "larks", those who wake up and go to sleep early, and "owls", those who stay up and wake up late. It is well established that a person's chronotype can affect their activities and health. However, less is known on the effects of chronotypes on the social behavior, even though it is evident that social interactions require coordinated timings. To study how chronotypes relate to social behavior, we use data collected using a smartphone app on a population of more than seven hundred volunteer students to simultaneously determine their chronotypes and social network structure. We find that owls maintain larger personal networks, albeit with less time spent per contact. On average, owls are more central in the social network of students than larks, frequently occupying the dense core of the network. Owls also display strong homophily, as seen in an unexpectedly large number of social ties connecting owls to owls.


title: Presentations for cusped arithmetic hyperbolic lattices

abstract: We present a general method to compute a presentation for any cusped arithmetic hyperbolic lattice $\Gamma$, applying a classical result of Macbeath to a suitable $\Gamma$-invariant horoball cover of the corresponding symmetric space. As applications we compute presentations for the Picard modular groups ${\rm PU}(2,1,\mathcal{O}_d)$ for $d=1,3,7$ and the quaternion hyperbolic lattice ${\rm PU}(2,1,\mathcal{H})$ with entries in the Hurwitz integer ring $\mathcal{H}$. The implementation of the method for these groups is computer-assisted.


title: A Voting-Based System for Ethical Decision Making

abstract: We present a general approach to automating ethical decisions, drawing on machine learning and computational social choice. In a nutshell, we propose to learn a model of societal preferences, and, when faced with a specific ethical dilemma at runtime, efficiently aggregate those preferences to identify a desirable choice. We provide a concrete algorithm that instantiates our approach; some of its crucial steps are informed by a new theory of swap-dominance efficient voting rules. Finally, we implement and evaluate a system for ethical decision making in the autonomous vehicle domain, using preference data collected from 1.3 million people through the Moral Machine website.


title: Flavor structure from misalignment of inner products in noncommutative geometry

abstract: In this letter, we consider an idea that induces flavor structure from inner products in noncommutative geometry. Assuming proper components of vectors $v_{(L,R) i}$ in enlarged representation space for fermions, we can induce the waterfall texture for Yukawa matrices retaining gauge interactions are universal. The hierarchy of the Yukawa interactions is a consequence of "misalignment" between the vectors $v_{Li}$ and $v_{Rj}$.


title: Spectral discretization errors in filtered subspace iteration

abstract: We consider filtered subspace iteration for approximating a cluster of eigenvalues (and its associated eigenspace) of a (possibly unbounded) selfadjoint operator in a Hilbert space. The algorithm is motivated by a quadrature approximation of an operator-valued contour integral of the resolvent. Resolvents on infinite dimensional spaces are discretized in computable finite-dimensional spaces before the algorithm is applied. This study focuses on how such discretizations result in errors in the eigenspace approximations computed by the algorithm. The computed eigenspace is then used to obtain approximations of the eigenvalue cluster. Bounds for the Hausdorff distance between the computed and exact eigenvalue clusters are obtained in terms of the discretization parameters within an abstract framework. A realization of the proposed approach for a model second-order elliptic operator using a standard finite element discretization of the resolvent is described. Some numerical experiments are conducted to gauge the sharpness of the theoretical estimates.


title: Effective dimension of some weighted pre-Sobolev spaces with dominating mixed partial derivatives

abstract: This paper considers two notions of effective dimension for quadrature in weighted pre-Sobolev spaces with dominating mixed partial derivatives. We begin by finding a ball in those spaces just barely large enough to contain a function with unit variance. If no function in that ball has more than $\varepsilon$ of its variance from ANOVA components involving interactions of order $s$ or more, then the space has effective dimension at most $s$ in the superposition sense. A similar truncation sense notion replaces the cardinality of the ANOVA component by the largest index it contains. Some Poincar\'e type inequalities are used to bound variance components by multiples of these space's squared norm and those in turn provide bounds on effective dimension. Very low effective dimension in the superposition sense holds for some spaces defined by product weights in which quadrature is strongly tractable. The superposition dimension is $O( \log(1/\varepsilon)/\log(\log(1/\varepsilon)))$ just like the superposition dimension used in the multidimensional decomposition method. Surprisingly, even spaces where all subset weights are equal, regardless of their cardinality or included indices, have low superposition dimension in this sense. This paper does not require periodicity of the integrands.