title: On Estimating Multi-Attribute Choice Preferences using Private Signals and Matrix Factorization

abstract: Revealed preference theory studies the possibility of modeling an agent's revealed preferences and the construction of a consistent utility function. However, modeling agent's choices over preference orderings is not always practical and demands strong assumptions on human rationality and data-acquisition abilities. Therefore, we propose a simple generative choice model where agents are assumed to generate the choice probabilities based on latent factor matrices that capture their choice evaluation across multiple attributes. Since the multi-attribute evaluation is typically hidden within the agent's psyche, we consider a signaling mechanism where agents are provided with choice information through private signals, so that the agent's choices provide more insight about his/her latent evaluation across multiple attributes. We estimate the choice model via a novel multi-stage matrix factorization algorithm that minimizes the average deviation of the factor estimates from choice data. Simulation results are presented to validate the estimation performance of our proposed algorithm.


title: Actions Speak Louder Than Goals: Valuing Player Actions in Soccer

abstract: Assessing the impact of the individual actions performed by soccer players during games is a crucial aspect of the player recruitment process. Unfortunately, most traditional metrics fall short in addressing this task as they either focus on rare actions like shots and goals alone or fail to account for the context in which the actions occurred. This paper introduces (1) a new language for describing individual player actions on the pitch and (2) a framework for valuing any type of player action based on its impact on the game outcome while accounting for the context in which the action happened. By aggregating soccer players' action values, their total offensive and defensive contributions to their team can be quantified. We show how our approach considers relevant contextual information that traditional player evaluation metrics ignore and present a number of use cases related to scouting and playing style characterization in the 2016/2017 and 2017/2018 seasons in Europe's top competitions.


title: Local Differential Privacy for Evolving Data

abstract: There are now several large scale deployments of differential privacy used to collect statistical information about users. However, these deployments periodically recollect the data and recompute the statistics using algorithms designed for a single use. As a result, these systems do not provide meaningful privacy guarantees over long time scales. Moreover, existing techniques to mitigate this effect do not apply in the "local model" of differential privacy that these systems use. In this paper, we introduce a new technique for local differential privacy that makes it possible to maintain up-to-date statistics over time, with privacy guarantees that degrade only in the number of changes in the underlying distribution rather than the number of collection periods. We use our technique for tracking a changing statistic in the setting where users are partitioned into an unknown collection of groups, and at every time period each user draws a single bit from a common (but changing) group-specific distribution. We also provide an application to frequency and heavy-hitter estimation.


title: Deep BCD-Net Using Identical Encoding-Decoding CNN Structures for Iterative Image Recovery

abstract: In "extreme" computational imaging that collects extremely undersampled or noisy measurements, obtaining an accurate image within a reasonable computing time is challenging. Incorporating image mapping convolutional neural networks (CNN) into iterative image recovery has great potential to resolve this issue. This paper 1) incorporates image mapping CNN using identical convolutional kernels in both encoders and decoders into a block coordinate descent (BCD) signal recovery method and 2) applies alternating direction method of multipliers to train the aforementioned image mapping CNN. We refer to the proposed recurrent network as BCD-Net using identical encoding-decoding CNN structures. Numerical experiments show that, for a) denoising low signal-to-noise-ratio images and b) extremely undersampled magnetic resonance imaging, the proposed BCD-Net achieves significantly more accurate image recovery, compared to BCD-Net using distinct encoding-decoding structures and/or the conventional image recovery model using both wavelets and total variation.


title: Universal qudit Hamiltonians

abstract: A family of quantum Hamiltonians is said to be universal if any other finite-dimensional Hamiltonian can be approximately encoded within the low-energy space of a Hamiltonian from that family. If the encoding is efficient, universal families of Hamiltonians can be used as universal analogue quantum simulators and universal quantum computers, and the problem of approximately determining the ground-state energy of a Hamiltonian from a universal family is QMA-complete. One natural way to categorise Hamiltonians into families is in terms of the interactions they are built from. Here we prove universality of some important classes of interactions on qudits ($d$-level systems): (1) We completely characterise the $k$-qudit interactions which are universal, if augmented with arbitrary 1-local terms. We find that, for all $k \geqslant 2$ and all local dimensions $d \geqslant 2$, almost all such interactions are universal aside from a simple stoquastic class. (2) We prove universality of generalisations of the Heisenberg model that are ubiquitous in condensed-matter physics, even if free 1-local terms are not provided. We show that the $SU(d)$ and $SU(2)$ Heisenberg interactions are universal for all local dimensions $d \geqslant 2$ (spin $\geqslant 1/2$), implying that a quantum variant of the Max-$d$-Cut problem is QMA-complete. We also show that for $d=3$ all bilinear-biquadratic Heisenberg interactions are universal. One example is the general AKLT model. (3) We prove universality of any interaction proportional to the projector onto a pure entangled state.


title: Semi-direct products involving $Sp_{2n}$ or $Spin_n$ with free algebras of symmetric invariants

abstract: This is a part of an ongoing project, the goal of which is to classify all semi-direct products $\mathfrak s=\mathfrak g{\ltimes} V$ such that $\mathfrak g$ is a simple Lie algebra, $V$ is a $\mathfrak g$-module, and $\mathfrak s$ has a free algebra of symmetric invariants. In this paper, we obtain such a classification for the representations of the orthogonal and symplectic algebras.


title: Capstone: Mobility Modeling on Smartphones to Achieve Privacy by Design

abstract: Sharing location traces with context-aware service providers has privacy implications. Location-privacy preserving mechanisms, such as obfuscation, anonymization and cryptographic primitives, have been shown to have impractical utility/privacy tradeoff. Another solution for enhancing user privacy is to minimize data sharing by executing the tasks conventionally carried out at the service providers' end on the users' smartphones. Although the data volume shared with the untrusted entities is significantly reduced, executing computationally demanding server-side tasks on resource-constrained smartphones is often impracticable. To this end, we propose a novel perspective on lowering the computational complexity by treating spatiotemporal trajectories as space-time signals. Lowering the data dimensionality facilitates offloading the computational tasks onto the digital-signal processors and the usage of the non-blocking signal-processing pipelines. While focusing on the task of user mobility modeling, we achieve the following results in comparison to the state of the art techniques: (i) mobility models with precision and recall greater than 80%, (ii) reduction in computational complexity by a factor of 2.5, and (iii) reduction in power consumption by a factor of 0.5. Furthermore, our technique does not rely on users' behavioral parameters that usually result in privacy-leakage and conclusive bias in the existing techniques. Using three real-world mobility datasets, we demonstrate that our technique addresses these weaknesses while formulating accurate user mobility models.


title: Towards Deep Representation Learning with Genetic Programming

abstract: Genetic Programming (GP) is an evolutionary algorithm commonly used for machine learning tasks. In this paper we present a method that allows GP to transform the representation of a large-scale machine learning dataset into a more compact representation, by means of processing features from the original representation at individual level. We develop as a proof of concept of this method an autoencoder. We tested a preliminary version of our approach in a variety of well-known machine learning image datasets. We speculate that this method, used in an iterative manner, can produce results competitive with state-of-art deep neural networks.


title: Characterization of generalized Petersen graphs that are Kronecker covers

abstract: The family of generalized Petersen graphs $G(n, k)$, introduced by Coxeter et al. [4] and named by Mark Watkins (1969), is a family of cubic graphs formed by connecting the vertices of a regular polygon to the corresponding vertices of a star polygon. The Kronecker cover $\mathrm{KC}(G)$ of a simple undirected graph $G$ is a a special type of bipartite covering graph of $G$, isomorphic to the direct (tensor) product of $G$ and $K_2$. We characterize all the members of generalized Petersen graphs that are Kronecker covers, and describe the structure of their respective quotients. We observe that some of such quotients are again generalized Petersen graphs, and describe all such pairs.


title: Analytic Analysis of Irregular Discrete Universes

abstract: In this work we investigate the dynamics of cosmological models with spherical topology containing up to 600 Schwarzschild black holes arranged in an irregular manner. We solve the field equations by tessellating the 3-sphere into eight identical cells, each having a single edge which is shared by all cells. The shared edge is enforced to be locally rotationally symmetric (LRS), thereby allowing for solving the dynamics to high accuracy along this edge. Each cell will then carry an identical (up to parity) configuration which can however have an arbitrarily random distribution. The dynamics of such models is compared to that of previous works on regularly distributed black holes as well as with the standard isotropic dust models of the FLRW type. The irregular models are shown to have richer dynamics than that of the regular models. The randomization of the distribution of the black holes is done both without bias and also with a certain clustering bias. The geometry of the initial configuration of our models is shown to be qualitatively different from the regular case in the way it approaches the isotropic model


title: Height of rational points on congruent number elliptic curves

abstract: We prove that a positive proportion of squarefree integers are congruent numbers such that the canonical height of the lowest non-torsion rational point on the corresponding elliptic curve satisfies a strong lower bound.


title: On the cosmological constant in the deformed Einstein-Cartan gauge gravity in De Donder-Weyl Hamiltonian formulation

abstract: A modification of the Einstein-Hilbert theory, the Covariant Canonical Gauge Gravity (CCGG), leads to a cosmological constant that represents the energy of the space-time continuum when deformed from its (A)dS ground state to a flat geometry. CCGG is based on the canonical transformation theory in the De Donder-Weyl (DW) Hamiltonian formulation. That framework modifies the Einstein-Hilbert Lagrangian of the free gravitational field by a quadratic Riemann-Cartan concomitant. The theory predicts a total energy-momentum of the system of space-time and matter to vanish, in line with the conjecture of a "Zero-Energy-Universe" going back to Lorentz (1916) and Levi-Civita (1917). Consequently a flat geometry can only exist in presence of matter where the bulk vacuum energy of matter, regardless of its value, is eliminated by the vacuum energy of space-time.% $\lambda_0$. The observed cosmological constant $\Lambda_{\mathrm{obs}}$ is found to be merely a small correction %of the order $10^{-120} \,\lambda_0$ attributable to deviations from a flat geometry and effects of complex dynamical geometry of space-time, namely torsion and possibly also vacuum fluctuations. That quadratic extension of General Relativity, anticipated already in 1918 by Einstein \cite{einstein18}, thus provides a significant and natural contribution to resolving the "cosmological constant problem".


title: BRST-BV quantization of gauge theories with global symmetries

abstract: We consider the general gauge theory with a closed irreducible gauge algebra possessing the non-anomalous global (super)symmetry in the case when the gauge fixing procedure violates the global invariance of classical action. The theory is quantized in the framework of BRST-BV approach in the form of functional integral over all fields of the configuration space. It is shown that the global symmetry transformations are deformed in the process of quantization and the full quantum action is invariant under such deformed global transformations in the configuration space. The deformed global transformations are calculated in an explicit form in the one-loop approximation.


title: Superconductivity of Bi-III phase of elemental Bismuth: insights from Muon-Spin Rotation and Density Functional Theory

abstract: Using muon-spin rotation the pressure-induced superconductivity in the Bi-III phase of elemental Bismuth (transition temperature $T_{\rm c}\simeq7.05$ K) was investigated. The Ginzburg-Landau parameter $\kappa=\lambda/\xi=30(6)$ ($\lambda$ is the magnetic penetration depth, $\xi$ is the coherence length) was estimated which is the highest among single element superconductors. The temperature dependence of the superconducting energy gap [$\Delta(T)$] reconstructed from $\lambda^{-2}(T)$ deviates from the weak-coupled BCS prediction. The coupling strength $2\Delta/k_{\rm B}T_{\rm c}\simeq 4.34$ was estimated thus implying that Bi-III stays within the strong coupling regime. The Density Functional Theory calculations suggest that superconductivity in Bi-III could be described within the Eliashberg approach with the characteristic phonon frequency $\omega_{\rm ln}\simeq 5.5$ meV. An alternative pairing mechanism to the electron-phonon coupling involves the possibility of Cooper pairing induced by the Fermi surface nesting.


title: A Testability Analysis Framework for Non-Functional Properties

abstract: This paper presents background, the basic steps and an example for a testability analysis framework for non-functional properties.


title: Arrival Time Distributions of Spin-1/2 Particles

abstract: The arrival time statistics of spin-1/2 particles governed by Pauli's equation, and defined by their Bohmian trajectories, show unexpected and very well articulated features. Comparison with other proposed statistics of arrival times that arise from either the usual (convective) quantum flux or from semiclassical considerations suggest testing the notable deviations in an arrival time experiment, thereby probing the predictive power of Bohmian trajectories. The suggested experiment, including the preparation of the wave functions, could be done with present-day experimental technology.


title: Mallows permutations as stable matchings

abstract: We show that the Mallows measure on permutations of $1,\ldots,n$ arises as the law of the unique Gale-Shapley stable matching of the random bipartite graph conditioned to be perfect, where preferences arise from a total ordering of the vertices but are restricted to the (random) edges of the graph. We extend this correspondence to infinite intervals, for which the situation is more intricate. We prove that almost surely every stable matching of the random bipartite graph obtained by performing Bernoulli percolation on the complete bipartite graph $K_{\mathbb{Z},\mathbb{Z}}$ falls into one of two classes: a countable family $(\sigma_n)_{n\in\mathbb{Z}}$ of tame stable matchings, in which the length of the longest edge crossing $k$ is $O(\log |k|)$ as $k\to\pm \infty$, and an uncountable family of wild stable matchings, in which this length is $\exp \Omega(k)$ as $k\to +\infty$. The tame stable matching $\sigma_n$ has the law of the Mallows permutation of $\mathbb{Z}$ (as constructed by Gnedin and Olshanski) composed with the shift $k\mapsto k+n$. The permutation $\sigma_{n+1}$ dominates $\sigma_{n}$ pointwise, and the two permutations are related by a shift along a random strictly increasing sequence.


title: Breaking the Loop: Recursive Proofs for Coinductive Predicates in Fibrations

abstract: The purpose of this paper is to develop and study recursive proofs of coinductive predicates. Such recursive proofs allow one to discover proof goals in the construction of a proof of a coinductive predicate, while still allowing the use of up-to techniques. This approach lifts the burden to guess invariants, like bisimulation relations, beforehand. Rather, they allow one to start with the sought-after proof goal and develop the proof from there until a point is reached, at which the proof can be closed through a recursion step. Proofs given in this way are both easier to construct and to understand, similarly to proofs given in cyclic proof systems or by appealing parameterised coinduction. In this paper, we develop a framework for recursive proofs of coinductive predicates that are given through fibrational predicate liftings. This framework is built on the so-called later modality, which has made its appearance in type theoretic settings before. In particular, we show the soundness and completeness of recursive proofs, we prove that compatible up-to techniques can be used as inference rules in recursive proofs, and provide some illustrating examples.


title: ILP-based Local Search for Graph Partitioning

abstract: Computing high-quality graph partitions is a challenging problem with numerous applications. In this paper, we present a novel meta-heuristic for the balanced graph partitioning problem. Our approach is based on integer linear programs that solve the partitioning problem to optimality. However, since those programs typically do not scale to large inputs, we adapt them to heuristically improve a given partition. We do so by defining a much smaller model that allows us to use symmetry breaking and other techniques that make the approach scalable. For example, in Walshaw's well-known benchmark tables we are able to improve roughly half of all entries when the number of blocks is high.


title: Coarse-Grained Simulation of DNA using LAMMPS

abstract: During the last decade coarse-grained nucleotide models have emerged that allow us to DNA and RNA on unprecedented time and length scales. Among them is oxDNA, a coarse-grained, sequence-specific model that captures the hybridisation transition of DNA and many structural properties of single- and double-stranded DNA. oxDNA was previously only available as standalone software, but has now been implemented into the popular LAMMPS molecular dynamics code. This article describes the new implementation and analyses its parallel performance. Practical applications are presented that focus on single-stranded DNA, an area of research which has been so far under-investigated. The LAMMPS implementation of oxDNA lowers the entry barrier for using the oxDNA model significantly, facilitates future code development and interfacing with existing LAMMPS functionality as well as other coarse-grained and atomistic DNA models.

categories: cond-mat.soft physics.bio-ph physics.comp-ph q-bio.BM


title: Stability and convergence of second order backward differentiation schemes for parabolic Hamilton-Jacobi-Bellman equations

abstract: We study a second order BDF (Backward Differentiation Formula) scheme for the numerical approximation of parabolic HJB (Hamilton-Jacobi-Bellman) equations. The scheme under consideration is implicit, non-monotone, and second order accurate in time and space. The lack of monotonicity prevents the use of well-known convergence results for solutions in the viscosity sense. In this work, we establish rigorous stability results in a general nonlinear setting as well as convergence results for some particular cases with additional regularity assumptions. While most results are presented for one-dimensional, linear parabolic and non-linear HJB equations, some results are also extended to multiple dimensions and to Isaacs equations. Numerical tests are included to validate the method.


title: Practical Pulse Engineering: Gradient Ascent Without Matrix Exponentiation

abstract: Since 2005 there has been a huge growth in the use of engineered control pulses to perform desired quantum operations in systems such as NMR quantum information processors. These approaches, which build on the original gradient ascent pulse engineering (GRAPE) algorithm, remain computationally intensive because of the need to calculate matrix exponentials for each time step in the control pulse. Here we discuss how the propagators for each time step can be approximated using the Trotter--Suzuki formula, and a further speed up achieved by avoiding unnecessary operations. The resulting procedure can give a substantial speed gain with negligible cost in propagator error, providing a more practical approach to pulse engineering.


title: Correlated pseudo-marginal schemes for time-discretised stochastic kinetic models

abstract: The challenging problem of conducting fully Bayesian inference for the reaction rate constants governing stochastic kinetic models (SKMs) is considered. Given the challenges underlying this problem, the Markov jump process representation is routinely replaced by an approximation based on a suitable time discretisation of the system of interest. Improving the accuracy of these schemes amounts to using an ever finer discretisation level, which in the context of the inference problem, requires integrating over the uncertainty in the process at a predetermined number of intermediate times between observations. Pseudo-marginal Metropolis-Hastings schemes are increasingly used, since for a given discretisation level, the observed data likelihood can be unbiasedly estimated using a particle filter. When observations are particularly informative an auxiliary particle filter can be implemented, by employing an appropriate construct to push the state particles towards the observations in a sensible way. Recent work in state-space settings has shown how the pseudo-marginal approach can be made much more efficient by correlating the underlying pseudo-random numbers used to form the likelihood estimate at the current and proposed values of the unknown parameters. We extend this approach to the time-discretised SKM framework by correlating the innovations that drive the auxiliary particle filter. We find that the resulting approach offers substantial gains in efficiency over a standard implementation.


title: On the search of the elusive Intermediate Mass Black-Holes

abstract: Ultra-Luminous X-ray sources (ULXs) are accreting black holes for which their X-ray properties have been seen to be different to the case of stellar-mass black hole binaries. For most of the cases their intrinsic energy spectra are well described by a cold accretion disc (thermal) plus a curved high-energy emission components. The mass of the black hole (BH) derived from the thermal disc component is usually in the range of 100-10^5 solar masses, which have led to the idea that this can represent strong evidence of the Intermediate Mass Black Holes (IMBH), proposed to exist by theoretical studies but with no firm detection (as a class) so far. Recent theoretical and observational developments are leading towards the idea that these sources are instead compact objects accreting at an unusual super-Eddington regime instead. On the other hand, gravitational waves have been seen to be a useful tool for finding (some of these) IMBHs. We give a brief overview about the recent advent of the discovery of gravitational waves and their relationship with these so far elusive IMBHs.


title: The Algebraic Approach to Duality: An Introduction

abstract: This survey article gives an elementary introduction to the algebraic approach to Markov process duality, as opposed to the pathwise approach. In the algebraic approach, a Markov generator is written as the sum of products of simpler operators, which each have a dual with respect to some duality function. We discuss at length the recent suggestion by Giardin\`a, Redig, and others, that it may be a good idea to choose these simpler operators in such a way that they form an irreducible representation of some known Lie algebra. In particular, we collect the necessary background on representations of Lie algebras that is crucial for this approach. We also discuss older work by Lloyd and Sudbury on duality functions of product form and the relation between intertwining and duality.


title: Exact results for the $O(N)$ model with quenched disorder

abstract: We use scale invariant scattering theory to exactly determine the lines of renormalization group fixed points for $O(N)$-symmetric models with quenched disorder in two dimensions. Random fixed points are characterized by two disorder parameters: a modulus that vanishes when approaching the pure case, and a phase angle. The critical lines fall into three classes depending on the values of the disorder modulus. Besides the class corresponding to the pure case, a second class has maximal value of the disorder modulus and includes Nishimori-like multicritical points as well as zero temperature fixed points. The third class contains critical lines that interpolate, as $N$ varies, between the first two classes. For positive $N$, it contains a single line of infrared fixed points spanning the values of $N$ from $\sqrt{2}-1$ to $1$. The symmetry sector of the energy density operator is superuniversal (i.e. $N$-independent) along this line. For $N=2$ a line of fixed points exists only in the pure case, but accounts also for the Berezinskii-Kosterlitz-Thouless phase observed in presence of disorder.


title: On an inverse source problem for the Biot's equations in electro-seismic imaging

abstract: Electro-seismic imaging is a novel hybrid imaging modality in geophysical exploration. This paper concerns an inverse source problem for Biot's equations that arise in electro-seismic imaging. Using the time reversal method, we derive an explicit reconstruction formula, which immediately gives the uniqueness and stability of the reconstructed solution.


title: Chow rings and gonality of general abelian varieties

abstract: We study the (covering) gonality of abelian varieties and their orbits of zero-cycles for rational equivalence. We show that any orbit for rational equivalence of zero-cycles of degree $k$ has dimension at most $k-1$. Building on the work of Pirola, we show that very general abelian varieties of dimension $g$ have covering gonality $k\geq f(g)$ where $f(g)$ grows like ${\rm log}\,g$. This answers a question asked by Bastianelli, De Poi, Ein, Lazarsfeld and B. Ullery. We also obtain results on the Chow ring of very general abelian varieties, eg. if $g\geq 2k-1$, for any divisor $D\in {\rm Pic}^0(A)$, $D^k$ is not a torsion cycle.


title: On the support of matrix coefficients of supercuspidal representations of the general linear group over a local non-archimedean field

abstract: We derive an upper bound on the support of matrix coefficients of suprecuspidal representations of the general linear group over a non-archimedean local field. The results are in par with those which can be obtained from the Bushnell--Kutzko classification of supercuspidal representations, but they are proved independently.


title: The lithium depletion boundary and the age of the Hyades cluster

abstract: Determination of the lithium depletion boundary (LDB), i.e., the observational limit below which the cores of very low-mass objects do not reach high enough temperature for Li destruction, has been used to obtain ages for several open clusters and stellar associations younger than 200 Myr, which until now has been considered as the practical upper limit on the range of applicability of this method. In this work we show that the LDB method can be extended to significant older ages than previously thought. Intermediate resolution optical spectra of six L-type candidate members in the Hyades cluster obtained using OSIRIS at the 10.4-m Gran Telescopio Canarias are presented. The Li~I 670.8~nm resonance doublet is clearly detected only in the two faintest and coolest of these objects, which are classified as L3.5 to L4 brown dwarf cluster members with luminosities around 10$^{-4}$ solar. Lithium depletion factors are estimated for our targets with the aid of synthetic spectra and they are compared with predictions from evolutionary models. A LDB age of 650$\pm$70 Myr for the Hyades provides a consistent description of our data using a set of state-of-the-art evolutionary models for brown dwarfs calculated by \citet{2015A&A..577..42}.


title: $B$ decay anomalies and dark matter from strong dynamics

abstract: Indications of lepton flavor universality violation in semileptonic $B$ decays to $K$ or $K^*$ and muons or electrons can be explained by leptoquark exchange. I present a model in which the leptoquark is a bound state of constituents charged under a new confining SU($N_{\rm HC}$) hypercolor interaction. The lightest neutral bound state in the theory is an asymmetric dark matter candidate, that might be directly detectable through its magnetic dipole moment interaction.


title: Comment on "A structural test for the conformal invariance of the critical 3d Ising model" by S. Meneses, S. Rychkov, J. M. Viana Parente Lopes and P. Yvernay. arXiv:1802.02319

abstract: In a recent preprint [ArXiv 1802.02319], Meneses et al. challenge our proof that scale invariance implies conformal invariance for the three-dimensional Ising model [B. Delamotte, M. Tissier and N. Wschebor, Phys. Rev. E 93 (2016), 012144.]. We refute their arguments. We also point out a mistake in their one-loop calculation of the dimension of the vector operator $V_\mu$ of lowest dimension which is not a total derivative.


title: Friction and Radiative Heat Exchange in a System of Two Parallel Plates Moving Sideways : Levin-Polevoy-Rytov Theory Revisited

abstract: It is shown that the fundamental results obtained in the works by Levine, Polevoi, Rytov (1980) and Polevoi (1990), based on the fluctuation-electromagnetic theory by Levine and Rytov, adequately describe the rate of radiative heat exchange and the frictional force in a system of two parallel thick plates in relative lateral motion. A numerically calculated friction force for good metals and thin gaps turns out to be by a factor 107 higher than earlier obtained by Polevoi and increases with increasing conductivity of the metals.


title: A Closed-loop controller to improve the Stability of Cascaded DC/DC Converters

abstract: Study of the buck converter and cascaded system considering the voltage mode controller has been done. First the small signal analysis of a buck dc/dc converter is presented and its mathematical representation has been showed. Then, the cascaded converter model regarding close loop impedances and voltage gain has been studied. The controller for this converter is proposed to stabilize the performance of the plant. The effectiveness of the proposed controller has been tested on a typical buck converter.


title: Performance analysis of a novel hybrid FSO / RF communication system

abstract: In this paper, a novel dual-hop relay-assisted hybrid Free Space Optical / Radio Frequency (FSO / RF) communication system is presented. In this structure an access point connects users within the building to the Base Station via a hybrid parallel FSO / RF link, this link is proposed firstly. Parallel combination of FSO and RF links and use of an access point, will increase capacity, reliability and data rate of the system. It is the first time that the effect of number of users on the performance of a dual-hop relay-assisted hybrid parallel FSO / RF system is investigated. FSO link is considered in Gamma-Gamma atmospheric turbulence with the effect of pointing error and RF link is considered in Rayleigh fading. For the first time, closed-form expressions are derived for Bit Error Rate (BER) and Outage Probability (P_out) of the proposed system. Derived expressions are verified through MATLAB simulations. It is shown that the performance of the proposed system is almost independent of atmospheric turbulence intensity, thereby when atmospheric turbulence strengthens, low power consumption is required for maintenance of the system performance. Hence the proposed structure is particularly suitable for mobile communication systems in which a small mobile battery supplies transmitter power. Also the proposed system performance of the system is preferable even at low signal to noise ratio (SNR). Therefore, proposed structure significantly reduces power consumption while maintaining performance of the system.


title: Fourier optics: basic concepts

abstract: Based on diffraction theory and the propagation of the light, Fourier optics is a powerful tool allowing the estimation of a visible-range imaging system to transfer the spatial frequency components of an object. The analyses of the imaging systems can thus be performed and the the performance retrieved. For a better understanding of the optical study, we present a short introduction of the Fourier optics and we review the mathematical treatment depending on the illumination conditions of the imaging system. Furthermore, resolution criteria based on Fourier optics are detailed. Also, the impact of aberrations on the imaging quality are discussed.


title: Marginal Stability of Sweet-Parker Type Current Sheets at Low Lundquist Numbers

abstract: Magnetohydrodynamic simulations have shown that a non-unique critical Lundquist number $S_c$ exists, hovering around $S_c \sim 10^4$, above which threshold Sweet-Parker type stationary reconnecting configurations become unstable to a fast tearing mode dominated by plasmoid generation. It is known that the flow along the sheet plays a stabilizing role, though a satisfactory explanation of the non-universality and variable critical Lundquist numbers observed is still lacking. Here we discuss this question using 2D linear MHD simulations and linear stability analyses of Sweet-Parker type current sheets in the presence of background stationary inflows and outflows at low Lundquist numbers ($S\le 10^4$). Simulations show that the inhomogeneous outflow stabilizes the current sheet by stretching the growing magnetic islands and at the same time evacuating the magnetic islands out of the current sheet. This limits the time during which fluctuations which begin at any given wave-length can remain unstable, rendering the instability non-exponential. We find that the linear theory based on the expanding-wavelength assumption works well for $S$ larger than $\sim 1000$. However we also find that the inflow and location of the initial perturbation also affect the stability threshold.


title: Transport-Based Pattern Theory: A Signal Transformation Approach

abstract: In many scientific fields imaging is used to relate a certain physical quantity to other dependent variables. Therefore, images can be considered as a map from a real-world coordinate system to the non-negative measurements being acquired. In this work we describe an approach for simultaneous modeling and inference of such data, using the mathematics of optimal transport. To achieve this, we describe a numerical implementation of the linear optimal transport transform, based on the solution of the Monge-Ampere equation, which uses Brenier's theorem to characterize the solution of the Monge functional as the derivative of a convex potential function. We use our implementation of the transform to compute a curl-free mapping between two images, and show that it is able to match images with lower error that existing methods. Moreover, we provide theoretical justification for properties of the linear optimal transport framework observed in the literature, including a theorem for the linear separation of data classes. Finally, we use our optimal transport method to empirically demonstrate that the linear separability theorem holds, by rendering non-linearly separable data as linearly separable following transform to transport space.


title: Cubic graphs, their Ehrhart quasi-polynomials, and a scissors congruence phenomenon

abstract: The scissors congruence conjecture for the unimodular group is an analogue of Hilbert's third problem, for the equidecomposability of polytopes. Liu and Osserman studied the Ehrhart quasi-polynomials of polytopes naturally associated to graphs whose vertices have degree one or three. In this paper, we prove the scissors congruence conjecture, posed by Haase and McAllister, for this class of polytopes. The key ingredient in the proofs is the nearest neighbor interchange on graphs and a naturally arising piecewise unimodular transformation.


title: On a certain identity involving the Gamma function

abstract: The goal of this paper is to prove the identity \begin{align}\sum \limits_{j=0}^{\lfloor s\rfloor}\frac{(-1)^j}{s^j}\eta_s(j)+\frac{1}{e^{s-1}s^s}\sum \limits_{j=0}^{\lfloor s\rfloor}(-1)^{j+1}\alpha_s(j)+\bigg(\frac{1-((-1)^{s-\lfloor s\rfloor +2})^{1/(s-\lfloor s\rfloor +2)}}{2}\bigg)\nonumber \\ \bigg(\sum \limits_{j=\lfloor s\rfloor +1}^{\infty}\frac{(-1)^j}{s^j}\eta_s(j)+\frac{1}{e^{s-1}s^s}\sum \limits_{j=\lfloor s\rfloor +1}^{\infty}(-1)^{j+1}\alpha_s(j)\bigg)=\frac{1}{\Gamma(s+1)},\nonumber \end{align}where \begin{align}\eta_s(j):=\bigg(e^{\gamma (s-j)}\prod \limits_{m=1}^{\infty}\bigg(1+\frac{s-j}{m}\bigg)\nonumber \\e^{-(s-j)/m}\bigg)\bigg(2+\log s-\frac{j}{s}+\sum \limits_{m=1}^{\infty}\frac{s}{m(s+m)}-\sum \limits_{m=1}^{\infty}\frac{s-j}{m(s-j+m)}\bigg), \nonumber \end{align}and \begin{align}\alpha_s(j):=\bigg(e^{\gamma (s-j)}\prod \limits_{m=1}^{\infty}\bigg(1+\frac{s-j}{m}\bigg)e^{-(s-j)/m}\bigg)\bigg(\sum \limits_{m=1}^{\infty}\frac{s}{m(s+m)}-\sum \limits_{m=1}^{\infty}\frac{s-j}{m(s-j+m)}\bigg),\nonumber \end{align}where $\Gamma(s+1)$ is the Gamma function defined by $\Gamma(s):=\int \limits_{0}^{\infty}e^{-t}t^{s-1}dt$ and $\gamma =\lim \limits_{n\longrightarrow \infty}\bigg(\sum \limits_{k=1}^{n}\frac{1}{k}-\log n\bigg)=0.577215664\cdots $ is the Euler-Mascheroni constant.


title: Closed, Two Dimensional Surface Dynamics

abstract: We present dynamic equations for two dimensional closed surfaces and analytically solve it for some simplified cases. We derive final equations for surface normal motions by two different ways. The solution of the equations of motions in normal direction indicates that any closed, two dimensional, homogeneous surface with time invariable surface energy density adopts constant mean curvature shape when it comes in equilibrium with environment. As an example, we apply the formalism to analyze equilibrium shapes of micelles and explain why they adopt spherical, lamellar and cylindrical shapes. We show that theoretical calculation for micellar optimal radius is in good agreement with all atom simulations and experiments.


title: High-Quality Prediction Intervals for Deep Learning: A Distribution-Free, Ensembled Approach

abstract: This paper considers the generation of prediction intervals (PIs) by neural networks for quantifying uncertainty in regression tasks. It is axiomatic that high-quality PIs should be as narrow as possible, whilst capturing a specified portion of data. We derive a loss function directly from this axiom that requires no distributional assumption. We show how its form derives from a likelihood principle, that it can be used with gradient descent, and that model uncertainty is accounted for in ensembled form. Benchmark experiments show the method outperforms current state-of-the-art uncertainty quantification methods, reducing average PI width by over 10%.


title: Some results relating to (p,q)-th relative Gol'dberg order and (p,q)-relative Gol'dberg type of entire functions of several variables

abstract: In this paper we introduce the notions of (p,q)-th relative Gol'dberg order and (p,q)-th relative Gol'dberg type of entire functions of several complex variables where p,q are any positive integers. Then we study some growth properties of entire functions of several complex variables on the basis of their (p,q)-th relative Gol'dberg order and (p,q)-th relative Gol'dberg type.


title: Self-organization principles of intracellular pattern formation

abstract: Dynamic patterning of specific proteins is essential for the spatiotemporal regulation of many important intracellular processes in procaryotes, eucaryotes, and multicellular organisms. The emergence of patterns generated by interactions of diffusing proteins is a paradigmatic example for self-organization. In this article we review quantitative models for intracellular Min protein patterns in E. coli, Cdc42 polarization in S. cerevisiae, and the bipolar PAR protein patterns found in C. elegans. By analyzing the molecular processes driving these systems we derive a theoretical perspective on general principles underlying self-organized pattern formation. We argue that intracellular pattern formation is not captured by concepts such as "activators"', "inhibitors", or "substrate-depletion". Instead, intracellular pattern formation is based on the redistribution of proteins by cytosolic diffusion, and the cycling of proteins between distinct conformational states. Therefore, mass-conserving reaction-diffusion equations provide the most appropriate framework to study intracellular pattern formation. We conclude that directed transport, e.g. cytosolic diffusion along an actively maintained cytosolic gradient, is the key process underlying pattern formation. Thus the basic principle of self-organization is the establishment and maintenance of directed transport by intracellular protein dynamics.


title: CytonMT: an Efficient Neural Machine Translation Open-source Toolkit Implemented in C++

abstract: This paper presents an open-source neural machine translation toolkit named CytonMT (https://github.com/arthurxlw/cytonMt). The toolkit is built from scratch only using C++ and NVIDIA's GPU-accelerated libraries. The toolkit features training efficiency, code simplicity and translation quality. Benchmarks show that CytonMT accelerates the training speed by 64.5% to 110.8% on neural networks of various sizes, and achieves competitive translation quality.


title: Constrained minimum Riesz and Green energy problems for vector measures associated with a generalized condenser

abstract: For a finite collection $\mathbf A=(A_i)_{i\in I}$ of locally closed sets in $\mathbb R^n$, $n\geqslant3$, with the sign $\pm1$ prescribed such that the oppositely charged plates are mutually disjoint, we consider the minimum energy problem relative to the $\alpha$-Riesz kernel $|x-y|^{\alpha-n}$, $\alpha\in(0,2]$, over positive vector Radon measures $\boldsymbol\mu=(\mu^i)_{i\in I}$ such that each $\mu^i$, $i\in I$, is carried by $A_i$ and normalized by $\mu^i(A_i)=a_i\in(0,\infty)$. We show that, though the closures of oppositely charged plates may intersect each other even in a set of nonzero capacity, this problem has a solution $\boldsymbol\lambda^{\boldsymbol\xi}_{\mathbf A}=(\lambda^i_{\mathbf A})_{i\in I}$ (also in the presence of an external field) if we restrict ourselves to $\boldsymbol\mu$ with $\mu^i\leqslant\xi^i$, $i\in I$, where the constraint $\boldsymbol\xi=(\xi^i)_{i\in I}$ is properly chosen. We establish the sharpness of the sufficient conditions on the solvability thus obtained, provide descriptions of the weighted vector $\alpha$-Riesz potentials of the solutions, single out their characteristic properties, and analyze the supports of the $\lambda^i_{\mathbf A}$, $i\in I$. Our approach is based on the simultaneous use of the vague topology and an appropriate semimetric structure defined in terms of the $\alpha$-Riesz energy on a set of vector measures associated with $\mathbf A$, as well as on the establishment of an intimate relationship between the constrained minimum $\alpha$-Riesz energy problem and a constrained minimum $\alpha$-Green energy problem, suitably formulated. The results are illustrated by examples.


title: Spin-1/2 anisotropic Heisenberg antiferromagnet with Dzyaloshinskii-Moriya interaction via mean-field approximation

abstract: The spin-1/2 anisotropic Heisenberg model with antiferromagnetic exchange interactions in the presence of a external magnetic field and a Dzyaloshinskii-Moriya interaction is studied by employing the usual mean-field approximation. The magnetic properties are obtained and it is shown that only second-order phase transitions take place for any values of the theoretical Hamiltonian parameters. Contrary to previous results from effective field theory, no anomalies have been observed at low temperatures. However, some re-entrancies still persist in some region of the phase diagram.


title: Existence of the gauge for fractional Laplacian Schr\"{o}dinger operators

abstract: Let $\Omega \subseteq \mathbb{R}^n$ be an open set, where $n \geq 2$. Suppose $\omega $ is a locally finite Borel measure on $\Omega$. For $\alpha \in (0,2)$, define the fractional Laplacian $(-\triangle )^{\alpha/2}$ via the Fourier transform on $\mathbb{R}^n$, and let $G $ be the corresponding Green's operator of order $\alpha$ on $\Omega$. Define $T(u) = G(u \omega).$ If $\Vert T \Vert_{L^2(\omega) \rightarrow L^2 (\omega)} <1$, we obtain a representation for the unique weak solution $u$ in the homogeneous Sobolev space $L^{\alpha/2, 2}_0 (\Omega)$ of \[ (-\triangle)^{\alpha/2} u = u \omega + \nu \,\,\, \mbox{on} \,\,\, \Omega, \,\,\, u=0 \,\,\, \mbox{on} \,\,\, \Omega^c, \] for $\nu$ in the dual Sobolev space $L^{-\alpha/2, 2} (\Omega)$. If $\Omega$ is a bounded $C^{1,1}$ domain, this representation yields matching exponential upper and lower pointwise estimates for the solution when $\nu = \chi_{\Omega}$. These estimates are used to study the existence of a solution $u_1$ (called the "gauge") of the integral equation $u_1=1+G(u_1 \omega)$ corresponding to the problem \[ (-\triangle)^{\alpha/2} u = u \omega \,\,\, \mbox{on} \,\,\, \Omega, \,\,\, u \geq 0 \,\,\, \mbox{on} \,\,\, \Omega, \,\,\, u=1 \,\,\, \mbox{on} \,\,\, \Omega^c . \] We show that if $\Vert T \Vert <1$, then $u_1$ always exists if $0<\alpha <1$. For $1 \leq \alpha <2$, a solution exists if the norm of $T$ is sufficiently small. We also show that the condition $\Vert T \Vert <1$ does not imply the existence of a solution if $1 < \alpha <2$.


title: Soft-photon corrections to the Bethe-Heitler process in the $\gamma p\rightarrow l^+l^-p$ reaction

abstract: We report on the calculation of first-order QED corrections for the $\gamma p\rightarrow l^+l^-p^{\prime}$ process. An upcoming experiment at MAMI (Mainz) aims to compare the cross sections of muon- and electron-pair production in this reaction to test lepton universality. Precise knowledge of the electromagnetic radiative corrections is needed for these measurements. As a first step, we present the leading QED radiative corrections in the soft-photon approximation when accounting for the finite lepton mass. For the kinematics at MAMI, we find corrections of the percent level for muons, and of order $10\%$ for electrons.


title: The parameterized complexity of finding a 2-sphere in a simplicial complex

abstract: We consider the problem of finding a subcomplex K' of a simplicial complex K such that K' is homeomorphic to the 2-dimensional sphere, S^2. We study two variants of this problem. The first asks if there exists such a K' with at most k triangles, and we show that this variant is W[1]-hard and, assuming ETH, admits no O(n^{o(sqrt(k))}) time algorithm. We also give an algorithm that is tight with regards to this lower bound. The second problem is the dual of the first, and asks if K' can be found by removing at most k triangles from K. This variant has an immediate O(3^k poly(|K|)) time algorithm, and we show that it admits a polynomial kernelization to O(k^2) triangles, as well as a polynomial compression to a weighted version with bit-size O(k log k).


title: Adaptive Sampling for Coarse Ranking

abstract: We consider the problem of active coarse ranking, where the goal is to sort items according to their means into clusters of pre-specified sizes, by adaptively sampling from their reward distributions. This setting is useful in many social science applications involving human raters and the approximate rank of every item is desired. Approximate or coarse ranking can significantly reduce the number of ratings required in comparison to the number needed to find an exact ranking. We propose a computationally efficient PAC algorithm LUCBRank for coarse ranking, and derive an upper bound on its sample complexity. We also derive a nearly matching distribution-dependent lower bound. Experiments on synthetic as well as real-world data show that LUCBRank performs better than state-of-the-art baseline methods, even when these methods have the advantage of knowing the underlying parametric model.


title: Wireless Expanders

abstract: This paper introduces an extended notion of expansion suitable for radio networks. A graph $G=(V,E)$ is called an $(\alpha_w, \beta_w)$-{wireless expander} if for every subset $S \subseteq V$ s.t. $|S|\leq \alpha_w \cdot |V|$, there exists a subset $S'\subseteq S$ s.t. there are at least $\beta_w \cdot |S|$ vertices in $V\backslash S$ adjacent in $G$ to exactly one vertex in $S'$. The main question we ask is the following: to what extent are ordinary expanders also good {wireless} expanders? We answer this question in a nearly tight manner. On the positive side, we show that any $(\alpha, \beta)$-expander with maximum degree $\Delta$ and $\beta\geq 1/\Delta$ is also a $(\alpha_w, \beta_w)$ wireless expander for $\beta_w = \Omega(\beta / \log (2 \cdot \min\{\Delta / \beta, \Delta \cdot \beta\}))$. Thus the wireless expansion is smaller than the ordinary expansion by at most a factor logarithmic in $\min\{\Delta / \beta, \Delta \cdot \beta\}$, which depends on the graph \emph{average degree} rather than maximum degree; e.g., for low arboricity graphs, the wireless expansion matches the ordinary expansion up to a constant. We complement this positive result by presenting an explicit construction of a "bad" $(\alpha, \beta)$-expander for which the wireless expansion is $\beta_w = O(\beta / \log (2 \cdot \min\{\Delta / \beta, \Delta \cdot \beta\})$. We also analyze the theoretical properties of wireless expanders and their connection to unique neighbor expanders, and demonstrate their applicability: Our results yield improved bounds for the {spokesmen election problem} that was introduced in the seminal paper of Chlamtac and Weinstein (1991) to devise efficient broadcasting for multihop radio networks. Our negative result yields a significantly simpler proof than that from the seminal paper of Kushilevitz and Mansour (1998) for a lower bound on the broadcast time in radio networks.


title: Primitive abundant and weird numbers with many prime factors

abstract: We give an algorithm to enumerate all primitive abundant numbers (briefly, PANs) with a fixed $\Omega$ (the number of prime factors counted with their multiplicity), and explicitly find all PANs up to $\Omega=6$, count all PANs and square-free PANs up to $\Omega=7$ and count all odd PANs and odd square-free PANs up to $\Omega=8$. We find primitive weird numbers (briefly, PWNs) with up to 16 prime factors, improving the previous results of [Amato-Hasler-Melfi-Parton] where PWNs with up to 6 prime factors have been given. The largest PWN we find has 14712 digits: as far as we know, this is the largest example existing, the previous one being 5328 digits long [Melfi]. We find hundreds of PWNs with exactly one square odd prime factor: as far as we know, only five were known before. We find all PWNs with at least one odd prime factor with multiplicity greater than one and $\Omega = 7$ and prove that there are none with $\Omega < 7$. Regarding PWNs with a cubic (or higher) odd prime factor, we prove that there are none with $\Omega\le 7$, and we did not find any with larger $\Omega$. Finally, we find several PWNs with 2 square odd prime factors, and one with 3 square odd prime factors. These are the first such examples.


title: Finite temperature electronic structure of Diamond and Silicon

abstract: The electron-phonon interaction contribution to the electronic energies is included in density functional total energy calculations with ab initio pseudopotentials via the Allen formalism [Phys. Rev. B 18, 5217 (1978)] to obtain temperature dependent electronic structure of diamond and silicon. This method allows us to obtain the thermally-averaged ab initio electronic structure in a straightforward and computationally inexpensive way. Our investigations on the finite temperature electronic structure of diamond and silicon lead to a new criterion, temperature transferability, which is required in the ab initio pseudopotentials for temperature dependent studies. The temperature transferability of the Troullier-Martins pseudopotentials used in this work is strongly dependent on the cut-off radius and the inclusion of the unbound 3d$^0$ state. The finite temperature indirect band gaps are highly sensitive to the choice of cut-off radius used in the pseudopotentials. The finite temperature band structures and density of states show that thermal vibrations affect the electron energies throughout the valence and conduction band. We compare our results on the band gap shifts with that due to the Debye-Waller term in the Allen-Heine theory and discuss the observed differences in the zero point and high temperature band gap shifts. Although, the electron energy shifts in the highest occupied valence band and lowest unoccupied conduction band enable to obtain the changes in the indirect and direct band gaps at finite temperatures, the shifts in other electronic levels with temperature enable investigations into the finite temperature valence charge distribution in the bonding region. Thus, we demonstrate that the Allen theory provides a simple and theoretically justified formalism to obtain finite temperature valence electron charge densities that go beyond the rigid pseudo-atom approximation.


title: Comparison of threshold-based algorithms for sparse signal recovery

abstract: Intensively growing approach in signal processing and acquisition, the Compressive Sensing approach, allows sparse signals to be recovered from small number of randomly acquired signal coefficients. This paper analyses some of the commonly used threshold-based algorithms for sparse signal reconstruction. Signals satisfy the conditions required by the Compressive Sensing theory. The Orthogonal Matching Pursuit, Iterative Hard Thresholding and Single Iteration Reconstruction algorithms are observed. Comparison in terms of reconstruction error and execution time is performed within the experimental part of the paper.


title: Algorithmic Information Dynamics of Persistent Patterns and Colliding Particles in the Game of Life

abstract: Without loss of generalisation to other systems, including possibly non-deterministic ones, we demonstrate the application of methods drawn from algorithmic information dynamics to the characterisation and classification of emergent and persistent patterns, motifs and colliding particles in Conway's Game of Life (GoL), a cellular automaton serving as a case study illustrating the way in which such ideas can be applied to a typical discrete dynamical system. We explore the issue of local observations of closed systems whose orbits may appear open because of inaccessibility to the global rules governing the overall system. We also investigate aspects of symmetry related to complexity in the distribution of patterns that occur with high frequency in GoL (which we thus call motifs) and analyse the distribution of these motifs with a view to tracking the changes in their algorithmic probability over time. We demonstrate how the tools introduced are an alternative to other computable measures that are unable to capture changes in emergent structures in evolving complex systems that are often too small or too subtle to be properly characterised by methods such as lossless compression and Shannon entropy.


title: The Gaussian Process Autoregressive Regression Model (GPAR)

abstract: Multi-output regression models must exploit dependencies between outputs to maximise predictive performance. The application of Gaussian processes (GPs) to this setting typically yields models that are computationally demanding and have limited representational power. We present the Gaussian Process Autoregressive Regression (GPAR) model, a scalable multi-output GP model that is able to capture nonlinear, possibly input-varying, dependencies between outputs in a simple and tractable way: the product rule is used to decompose the joint distribution over the outputs into a set of conditionals, each of which is modelled by a standard GP. GPAR's efficacy is demonstrated on a variety of synthetic and real-world problems, outperforming existing GP models and achieving state-of-the-art performance on established benchmarks.


title: Unitary Scheme Model Calculations of the Ground- and Excited-State Characteristics of 3H and 4He

abstract: The ground and excited states of the 3H and 4He nuclei are studied in the framework of group-theoretical methods. Basis functions of the unitary scheme model corresponding to even numbers of quanta of excitation in the range from zero to twenty are constructed for the even-parity states of these two nuclei, and bases from one to nineteen are constructed for the odd-parity states of the 4He nucleus. The ground-state and first-excited-state energies and wave functions, the ground S-, P- and D-state probabilities, the root-mean-square radius and the magnetic dipole moment of Triton are calculated. Furthermore, for the 4He nucleus, the spectrum and the wave functions, the ground S-, P- and D-state probabilities, the root-mean-square radius and the total integral cross section of the dipole electric transition accompanying the photoabsorption of gamma quanta by this nucleus are calculated. The GPT and AV8' two-body interactions and the Urbana nine three-nucleon interaction are used in these investigations. Moreover, the convergence of the calculations is examined by incrementally extrapolating the nuclear characteristics calculated for N less than or equal to twenty to N equals thirty.


title: Comment on the Bekenstein bound

abstract: We propose a rigorous derivation of the Bekenstein upper limit for the entropy/information that can be contained by a physical system in a given finite region of space with given finite energy. The starting point is the observation that the derivation of such a bound provided by Casini [6] is similar to the description of the black hole incremental free energy that had been given by the first named author [23]. The approach here is different but close in the spirit to [6]. Our bound is obtained by operator algebraic methods, in particular Connes' bimodules, Tomita-Takesaki modular theory and Jones' index are essential ingredients inasmuch as the von Neumann algebras in question are typically of type III. We rely on the general mathematical framework, recently set up in [26], concerning quantum information of infinite systems.


title: The Planetary Nebula Luminosity Function (PNLF): Contamination from Supernova Remnants

abstract: The planetary nebula luminosity function (PNLF) has been used as an extragalactic distance indicator since the 1980's, but there are still unsolved problems associated with its use. One of the most serious involves PNLF distances beyond ~ 10 Mpc, which tend to be slightly smaller than those of other methods. We consider the implications of previous spectroscopic investigations that found that several of the brightest planetary nebula (PN) candidates in M74 are actually compact supernova remnants (SNRs). Using narrow-band imaging data from the KPNO 4-m telescope, we measure the [O III] $\lambda$5007 and H$\alpha$ fluxes of all the known SNRs in M31 and M33, and test whether those objects could be misidentified as bright PNe at distances beyond ~ 10 Mpc. Our results suggest that compact SNRs are not an important source of contamination in photometric surveys for extragalactic PNe.


title: Stochastic compressible Euler equations and inviscid limits

abstract: We prove the existence of a unique local strong solution to the stochastic compressible Euler system with nonlinear multiplicative noise. This solution exists up to a positive stopping time and is strong in both the PDE and probabilistic sense. Based on this existence result, we study the inviscid limit of the stochastic compressible Navier--Stokes system. As the viscosity tends to zero, any sequence of finite energy weak martingale solutions converges to the compressible Euler system.


title: Leader-follower based Coalition Formation in Large-scale UAV Networks, A Quantum Evolutionary Approach

abstract: The problem of decentralized multiple Point of Interests (PoIs) detection and associated task completion in an unknown environment with multiple resource-constrained and self-interested Unmanned Aerial Vehicles (UAVs) is studied. The UAVs form several coalitions to efficiently complete the compound tasks which are impossible to be performed individually. The objectives of such coalition formation are to firstly minimize resource consumption in completing the encountered tasks on time, secondly to enhance the reliability of the coalitions, and lastly in segregating the most trusted UAVs amid the self interested of them. As many previous publications have merely focused on minimizing costs, this study considers a multi-objective optimization coalition formation problem that considers the three aforementioned objectives. In doing so, a leader-follower- inspired coalition formation algorithm amalgamating the three objectives to address the problem of the computational complexity of coalition formation in large-scale UAV networks is proposed. This algorithm attempts to form the coalitions with minimally exceeding the required resources for the encountered tasks while maximizing the number of completed tasks. The proposed algorithm is based on Quantum Evolutionary Algorithms(QEA) which are a combination of quantum computing and evolutionary algorithms. Results from simulations show that the proposed algorithm significantly outperforms the existing coalition formation algorithms such as merge-and-split and a famous multi-objective genetic algorithm called NSGA-II.


title: Adjoint Sensitivity Analysis of Hybrid Multibody Dynamical Systems

abstract: Sensitivity analysis of multibody systems computes the derivatives of general cost functions that depend on the system solution with respect to parameters or initial conditions. This work develops adjoint sensitivity analysis for hybrid multibody dynamic systems. Hybrid systems are characterized by trajectories that are piecewise continuous in time, with finitely-many discontinuities being caused by events such as elastic/inelastic impacts or sudden changes in constraints. The corresponding direct and adjoint sensitivity variables are also discontinuous at the time of events. The framework discussed herein uses a jump sensitivity matrix to relate the jump conditions for the direct and adjoint sensitivities before and after the time event, and provides analytical jump equations for the adjoint variables. The theoretical framework for sensitivities for hybrid systems is validated on a five-bar mechanism with non-smooth contacts.


title: On configuration space, Born's rule and ontological states

abstract: It is shown how configuration space, possibly encompassing ordinary spatial structures, Born's rule, and ontological states aiming to address an underlying reality beyond Quantum Mechanics relate to each other in models of Hamiltonian cellular automata.


title: Evolutionary games: natural selection of strategies

abstract: In this paper, I model and study the process of natural selection between all possible mixed strategies in classical two-player two-strategy games. I derive and solve an equation that is a natural generalization of the Taylor-Jonker replicator equation that describes dynamics of pure strategy frequencies. I then investigate the evolution of not only frequencies of pure strategies but also of total distribution of mixed strategies. I show that the process of natural selection of strategies for all games obeys the dynamical Principle of minimum of information gain.


title: Neural Architecture Search with Bayesian Optimisation and Optimal Transport

abstract: Bayesian Optimisation (BO) refers to a class of methods for global optimisation of a function $f$ which is only accessible via point evaluations. It is typically used in settings where $f$ is expensive to evaluate. A common use case for BO in machine learning is model selection, where it is not possible to analytically model the generalisation performance of a statistical model, and we resort to noisy and expensive training and validation procedures to choose the best model. Conventional BO methods have focused on Euclidean and categorical domains, which, in the context of model selection, only permits tuning scalar hyper-parameters of machine learning algorithms. However, with the surge of interest in deep learning, there is an increasing demand to tune neural network \emph{architectures}. In this work, we develop NASBOT, a Gaussian process based BO framework for neural architecture search. To accomplish this, we develop a distance metric in the space of neural network architectures which can be computed efficiently via an optimal transport program. This distance might be of independent interest to the deep learning community as it may find applications outside of BO. We demonstrate that NASBOT outperforms other alternatives for architecture search in several cross validation based model selection tasks on multi-layer perceptrons and convolutional neural networks.


title: Robust Optimal Eco-driving Control with Uncertain Traffic Signal Timing

abstract: This paper proposes a robust optimal eco-driving control strategy considering multiple signalized intersections with uncertain traffic signal timing. A spatial vehicle velocity profile optimization formulation is developed to minimize the global fuel consumption, with driving time as one state variable. We introduce the concept of effective red-light duration (ERD), formulated as a random variable, to describe the feasible passing time through signalized intersections. A chance constraint is appended to the optimal control problem to incorporate robustness with respect to uncertain signal timing. The optimal eco-driving control problem is solved via dynamic programming (DP). Simulation results demonstrate that the optimal eco-driving can save fuel consumption by 50-57% while maintaining arrival time at the same level, compared with a modified intelligent driver model as the benchmark. The robust formulation significantly reduces traffic intersection violations, in the face of uncertain signal timing, with small sacrifice on fuel economy compared to a non-robust approach.


title: Construction of Milnorian representations

abstract: We prove a partial converse to the main theorem of the author's previous paper "Proper affine actions: a sufficient criterion" (submitted; available at arXiv:1612.08942). More precisely, let $G$ be a semisimple real Lie group with a representation $\rho$ on a finite-dimensional real vector space $V$, that does not satisfy the criterion from the previous paper. Assuming that $\rho$ is irreducible and under some additional assumptions on $G$ and $\rho$, we then prove that there does not exist a group of affine transformations acting properly discontinuously on $V$ whose linear part is Zariski-dense in $\rho(G)$.


title: A Note on "Quantum Algorithm for Linear Systems of Equations"

abstract: Recently, an efficient quantum algorithm for linear systems of equations introduced by Harrow, Hassidim, and Lloyd, has received great concern from the academic community. However, the error and complexity analysis for this algorithm seems so complicated that it may not be applicable to other filter functions for other tasks. In this note, a concise proof is proposed. We hope that it may inspire some novel HHL-based algorithms that can compute $F(A)|b\rangle$ for any computable $F$.


title: Ultrafast magnetization dynamics in pure and doped Heusler and inverse Heusler alloys

abstract: By using a multiscale approach based on first-principles density functional theory combined with atomistic spin dynamics, we investigate the electronic structure and magnetization dynamics of an inverse Heusler and a Heusler compound and their alloys, i. e. Mn$_{2-x}Z_x$CoAl and Mn$_{2-x}Z_x$VAl, where $Z$ = Mo, W, Os and Ru, respectively. A signature of the ferrimagnetic ordering of Mn$_{2}$CoAl and Mn$_{2}$VAl Heusler alloys is reflected in the calculated Heisenberg exchange constants. They decay very rapidly with the interatomic distance and have short range, which is a consequence of the existence of the finite gap in the minority spin band. The calculated Gilbert damping parameter of both Mn$_2$CoAl and Mn$_2$VAl is high compared to other half-metals, but interestingly in the particular case of the inverse Mn$_{2}$CoAl alloys and due to the spin-gapless semiconducting property, the damping parameters decrease with the doping concentration in clear contradiction to the general trend. Atomistic spin dynamics simulations predict ultrafast magnetisation switching in Mn$_{2}$CoAl and Mn$_{2}$VAl under the influence of an external magnetic field, starting from a threshold field of $2\text{T}$. Our overall finding extends with Heusler and inverse Heusler alloys, the class of materials that exhibits laser induced magnetic switching.


title: On Davenport constant of finite abelian groups

abstract: $G$ be an additive finite abelian group. The Davenport constant $\mathsf D(G)$ is the smallest integer $t$ such that every sequence (multiset) $S$ over $G$ of length $|S|\ge t$ has a non-empty zero-sum subsequence. Recently, B. Girard proved that for every fixed integer $r > 1$ the Davenport constant $\mathsf D(C_n^r)$ is asymptotic to $rn$ when $n$ tends to infinity. In this paper, for every fixed positive integer $r$, we prove that $$\mathsf D(C_n^r)=rn+O(\frac{n}{\ln n}).$$ This is an explicit version of the above result of B. Girard. Furthermore, we can get better estimates of the error term for some $n$ of special types. Finally, we get an asymptotic result for some finite abelian groups of special types. Our proof combines a classical argument in the zero-sum theory together with some basic tools and results from analytic number theory.


title: Topological nature of spinons and holons: Elementary excitations from matrix product states with conserved symmetries

abstract: We develop variational matrix product state (MPS) methods with symmetries to determine dispersion relations of one dimensional quantum lattices as a function of momentum and preset quantum number. We test our methods on the XXZ spin chain, the Hubbard model and a non-integrable extended Hubbard model, and determine the excitation spectra with a precision similar to the one of the ground state. The formulation in terms of quantum numbers makes the topological nature of spinons and holons very explicit. In addition, the method also enables an easy and efficient direct calculation of the necessary magnetic field or chemical potential required for a certain ground state magnetization or particle density.


title: DarkSide-50 532-day Dark Matter Search with Low-Radioactivity Argon

abstract: The DarkSide-50 direct-detection dark matter experiment is a dual-phase argon time projection chamber operating at Laboratori Nazionali del Gran Sasso. This paper reports on the blind analysis of a (16,660+-270) kg d exposure using a target of low-radioactivity argon extracted from underground sources. We find no events in the dark matter selection box and set a 90% C.L. upper limit on the dark matter-nucleon spin-independent cross section of 1.14E-44 cm^2 (3.78E-44 cm^2, 3.43E-43 cm^2) for a WIMP mass of 100 GeV/c^2 (1 TeV/c^2, 10 TeV/c^2).


title: A Parametric MPC Approach to Balancing the Cost of Abstraction for Differential-Drive Mobile Robots

abstract: When designing control strategies for differential-drive mobile robots, one standard tool is the consideration of a point at a fixed distance along a line orthogonal to the wheel axis instead of the full pose of the vehicle. This abstraction supports replacing the non-holonomic, three-state unicycle model with a much simpler two-state single-integrator model (i.e., a velocity-controlled point). Yet this transformation comes at a performance cost, through the robot's precision and maneuverability. This work contains derivations for expressions of these precision and maneuverability costs in terms of the transformation's parameters. Furthermore, these costs show that only selecting the parameter once over the course of an application may cause an undue loss of precision. Model Predictive Control (MPC) represents one such method to ameliorate this condition. However, MPC typically realizes a control signal, rather than a parameter, so this work also proposes a Parametric Model Predictive Control (PMPC) method for parameter and sampling horizon optimization. Experimental results are presented that demonstrate the effects of the parameterization on the deployment of algorithms developed for the single-integrator model on actual differential-drive mobile robots.


title: Asymptotics of Hitchin's metric on the Hitchin section

abstract: We consider Hitchin's hyperk\"ahler metric $g$ on the moduli space $\mathcal{M}$ of degree zero $\mathrm{SL}(2)$-Higgs bundles over a compact Riemann surface. It has been conjectured that, when one goes to infinity along a generic ray in $\mathcal{M}$, $g$ converges to an explicit "semiflat" metric $g^{\mathrm{sf}}$, with an exponential rate of convergence. We show that this is indeed the case for the restriction of $g$ to the tangent bundle of the Hitchin section $\mathcal{B} \subset \mathcal{M}$.


title: Inference of gene loss rates after whole genome duplications at early vertebrates through ancient genome reconstructions

abstract: The famous 2R hypothesis was first proposed by Susumu Ohno in 1970. It states that the two whole genome duplications had shaped the genome of early vertebrates. The most convincing evidence for 2R hypothesis comes from the 4:1 ratio chromosomal regions that have preserved both gene content and order in vertebrates compared with closely related. However, due to the shortage of such strict evidence, the 2R hypothesis is still under debates. Here, we present a combined perspective of phylogenetic and genomic homology to revisit the hypothesis of 2R whole genome duplications. Ancestral vertebrate genomes as well as ancient duplication events were created from 17 extant vertebrate species. Extant descendants from the duplication events at early vertebrates were extracted and reorganized to partial genomes. We then examined the gene order based synteny, and projected back to phylogenetic gene trees for examination of synteny evidence of the reconstructed early vertebrate genes. We identified 7877 ancestral genes that were created from 3026 duplication events at early vertebrates, and more than 50% of the duplication events show synteny evidence. Thus, our reconstructions provide very strong evidence for the 2R hypothesis. We also reconstructed the genome of early vertebrates, and built a model of the gene gains and losses in early vertebrates. We estimated that there were about 12,000 genes in early vertebrates before 2R, and the probability of a random gene get lost after the first round of whole genome duplication is around 0.45, and the probability of a random gene get lost after the second round of whole genome duplication is around 0.55. This research provides convincing evidence for the 2R hypothesis, and may provide further insights in vertebral evolution. Data availability: https://github.com/haimingt/Ohnologs-and-2R-WGD


title: Regularity of bicyclic Graphs and their powers

abstract: Let $I(G)$ be the edge ideal of a bicyclic graph. In this paper, we characterize the Castelnuovo-Mumford regularity of $I(G)$ in terms of the induced matching number of $G$. For the base case of this family of graphs, i.e. dumbbell graph, we explicitly compute the induced matching number. Moreover, we prove that $ \operatorname{reg}(I(G)^q)=2q+\operatorname{reg}(I(G))-2$, for all $ q\geq 1 $, when $ G $ is a dumbbell graph with a connecting path having no more than two vertices.


title: Endomorphisms of Koszul complexes: formality and application to deformation theory

abstract: We study the differential graded Lie algebra of endomorphisms of the Koszul resolution of a regular sequence on a unitary commutative $K$-algebra $R$ and we prove that it is homotopy abelian over $K$, while it is generally not formal over $R$. We apply this result to prove an annihilation theorem for obstructions of (derived) deformations of locally complete intersection ideal sheaves on projective schemes.


title: Automation of gene function prediction through modeling human curators' decisions in GO phylogenetic annotation project

abstract: The Gene Ontology Consortium launched the GO-PAINT project (Phylogenetic Annotation and INference Tool) 9 years ago and is currently being used in the GO Reference Genome Annotation Project to support inference of GO function terms (molecular function, cellular component and biological process) by homology. PAINT uses a phylogenetic model to infer gene function by homology, a process that requires manual curation of experienced biocurators. Tremendous amount of time and efforts have been spent on the GO-PAINT project yielding more than 4000 fully annotated phylogenetic families with more than 170,000 annotations. These preliminary data have thus enabled potential algorithmic representation and automatic solvation of the additional 9000 unannoated phylogenetic families. Here we present an automated pipeline for phylogenetic annotation and inference, which simulates the standard annotation procedures of curators and models the curators' decisions during the manual curation process. The pipeline has been built into the newest version of PAINT software available at http://www.pantherdb.org/downloads/index.jsp. The standalone automation pipeline and datasets are available at https://github.com/haimingt/GO-PAINT-automation


title: Information gain and loss for a quantum Maxwell's demon

abstract: We use continuous weak measurements of a driven superconducting qubit to experimentally study the information dynamics of a quantum Maxwell's demon. We show how information gained by a demon who can track single quantum trajectories of the qubit can be converted into work using quantum coherent feedback. We verify the validity of a quantum fluctuation theorem with feedback by utilizing information obtained along single trajectories. We demonstrate, in particular, that quantum backaction can lead to a loss of information in imperfect measurements. We furthermore probe the transition between information gain and loss by varying the initial purity of the qubit.


title: The Fifth Force in the Local Cosmic Web

abstract: Extensions of the standard models of particle physics and cosmology often lead to long-range fifth forces with properties dependent on gravitational environment. Fifth forces on astrophysical scales are best studied in the cosmic web where perturbation theory breaks down. We present constraints on chameleon- and symmetron-screened fifth forces with Yukawa coupling and megaparsec range -- as well as unscreened fifth forces with differential coupling to galactic mass components -- by searching for the displacements they predict between galaxies' stars and gas. Taking data from the Alfalfa HI survey, identifying galaxies' gravitational environments with the maps of Desmond et al. (2018a) and forward-modelling with a Bayesian likelihood framework, we set upper bounds on fifth-force strength relative to Newtonian gravity from $\Delta G/G_N < \text{few} \: \times 10^{-4}$ for range $\lambda_C = 50$ Mpc, to $\Delta G/G_N \lesssim 0.1$ for $\lambda_C = 500$ kpc. In $f(R)$ gravity this requires $f_{R0} < \text{few} \: \times \: 10^{-8}$. The analogous bounds without screening are $\Delta G/G_N < \text{few} \: \times 10^{-4}$ and $\Delta G/G_N < \text{few} \times 10^{-3}$. These are the tightest and among the only fifth-force constraints on galaxy scales. We show how our results may be strengthened with future survey data and identify the key features of an observational programme for furthering fifth-force tests beyond the Solar System.


title: AutoPrognosis: Automated Clinical Prognostic Modeling via Bayesian Optimization with Structured Kernel Learning

abstract: Clinical prognostic models derived from largescale healthcare data can inform critical diagnostic and therapeutic decisions. To enable off-theshelf usage of machine learning (ML) in prognostic research, we developed AUTOPROGNOSIS: a system for automating the design of predictive modeling pipelines tailored for clinical prognosis. AUTOPROGNOSIS optimizes ensembles of pipeline configurations efficiently using a novel batched Bayesian optimization (BO) algorithm that learns a low-dimensional decomposition of the pipelines high-dimensional hyperparameter space in concurrence with the BO procedure. This is achieved by modeling the pipelines performances as a black-box function with a Gaussian process prior, and modeling the similarities between the pipelines baseline algorithms via a sparse additive kernel with a Dirichlet prior. Meta-learning is used to warmstart BO with external data from similar patient cohorts by calibrating the priors using an algorithm that mimics the empirical Bayes method. The system automatically explains its predictions by presenting the clinicians with logical association rules that link patients features to predicted risk strata. We demonstrate the utility of AUTOPROGNOSIS using 10 major patient cohorts representing various aspects of cardiovascular patient care.


title: Experimental single-photon pulse characterization by electro-optic shearing interferometry

abstract: The ability to characterize the complete quantum state of light is essential for both fundamental and applied science. For single photons the quantum state is provided by the mode that it occupies. The spectral temporal mode structure of light has recently emerged as an essential means for quantum information science. Here we experimentally demonstrate a self-referencing technique to completely determine the pulse-mode structure of single photons by means of spectral shearing interferometry. We detail the calibration and resolution of the measurement and discuss challenges and critical requirements for future advances of this method.


title: Distributed Symmetry-Breaking Algorithms for Congested Cliques

abstract: The {Congested Clique} is a distributed-computing model for single-hop networks with restricted bandwidth that has been very intensively studied recently. It models a network by an $n$-vertex graph in which any pair of vertices can communicate one with another by transmitting $O(\log n )$ bits in each round. Various problems have been studied in this setting, but for some of them the best-known results are those for general networks. In this paper we devise significantly improved algorithms for various symmetry-breaking problems, such as forests-decompositions, vertex-colorings, and maximal independent set. We analyze the running time of our algorithms as a function of the arboricity $a$ of a clique subgraph that is given as input. Our algorithms are especially efficient in Trees, planar graphs, graphs with constant genus, and many other graphs that have bounded arboricity, but unbounded size. We obtain $O(a)$-forest-decomposition algorithm with $O(\log a)$ time that improves the previously-known $O(\log n)$ time, $O(a^{2 + \epsilon})$-coloring in $O(\log^* n)$ time that improves upon an $O(\log n)$-time algorithm, $O(a)$-coloring in $O(a^{\epsilon})$-time that improves upon several previous algorithms, and a maximal independent set algorithm with $O(\sqrt a)$ time that improves at least quadratically upon the state-of-the-art for small and moderate values of $a$. Those results are achieved using several techniques. First, we produce a forest decomposition with a helpful structure called {$H$-partition} within $O(\log a)$ rounds. In general graphs this structure requires $\Theta(\log n)$ time, but in Congested Cliques we are able to compute it faster. We employ this structure in conjunction with partitioning techniques that allow us to solve various symmetry-breaking problems efficiently.


title: Real-Time Dense Stereo Matching With ELAS on FPGA Accelerated Embedded Devices

abstract: For many applications in low-power real-time robotics, stereo cameras are the sensors of choice for depth perception as they are typically cheaper and more versatile than their active counterparts. Their biggest drawback, however, is that they do not directly sense depth maps; instead, these must be estimated through data-intensive processes. Therefore, appropriate algorithm selection plays an important role in achieving the desired performance characteristics. Motivated by applications in space and mobile robotics, we implement and evaluate a FPGA-accelerated adaptation of the ELAS algorithm. Despite offering one of the best trade-offs between efficiency and accuracy, ELAS has only been shown to run at 1.5-3 fps on a high-end CPU. Our system preserves all intriguing properties of the original algorithm, such as the slanted plane priors, but can achieve a frame rate of 47fps whilst consuming under 4W of power. Unlike previous FPGA based designs, we take advantage of both components on the CPU/FPGA System-on-Chip to showcase the strategy necessary to accelerate more complex and computationally diverse algorithms for such low power, real-time systems.


title: An algorithm for determining the rotation count of pulsars

abstract: We present here a simple, systematic method for determining the correct global rotation count of a radio pulsar; an essential step for the derivation of an accurate phase-coherent ephemeris. We then build on this method by developing a new algorithm for determining the global rotational count for pulsars with sparse timing data sets. This makes it possible to obtain phase-coherent ephemerides for pulsars for which this has been impossible until now. As an example, we do this for PSR J0024-7205aa, an extremely faint MSP recently discovered in the globular cluster 47 Tucanae. This algorithm has the potential to significantly reduce the number of observations and the amount of telescope time needed to follow up on new pulsar discoveries.


title: Acoustic microstreaming near a plane wall due to a pulsating free/coated bubble: Nyborg result revisited and vortices

abstract: Acoustic microstreaming due to an oscillating microbubble is analytically investigated to obtain the circular streaming motion adjacent to a nearby wall. Classical theory due to Nyborg is carefully derived in the radial coordinates. The theory is used to obtain the flow field and the vortical motion caused by the microbubble motion. The length of the vertices are decreasing when the microbubble is excited at distances close to the rigid wall, while the maximum shear stress is increasing.


title: Heegaard splittings of graph manifolds

abstract: In this paper we give a method to construct Heegaard splittings of oriented graph manifolds with orientable bases. A graph manifold is a closed $3$-manifold admitting only Seifert-fibered pieces in its Jaco-Shalen decomposition; for technical reasons, we restrict our attention to the fully oriented case, i.e. both the pieces and the bases are oriented.


title: Inverse Design of Discrete Mechanical Metamaterials

abstract: Mechanical and phononic metamaterials exhibiting negative elastic moduli, gapped vibrational spectra, or topologically protected modes enable precise control of structural and acoustic functionalities. While much progress has been made in their experimental and theoretical characterization, the inverse design of mechanical metamaterials with arbitrarily programmable spectral properties and mode localization remains an unsolved problem. Here, we present a flexible computational inverse-design framework that allows the efficient tuning of one or more gaps at nearly arbitrary positions in the spectrum of discrete phononic metamaterial structures. The underlying algorithm optimizes the linear response of elastic networks directly, is applicable to ordered and disordered structures, scales efficiently in 2D and 3D, and can be combined with a wide range of numerical optimization schemes. We illustrate the broad practical potential of this approach by designing mechanical bandgap switches that open and close pre-programmed spectral gaps in response to an externally applied stimulus such as shear or compression. We further show that the designed structures can host topologically protected edge modes, and validate the numerical predictions through explicit 3D finite element simulations of continuum elastica with experimentally relevant material parameters. Generally, this network-based inverse design paradigm offers a direct pathway towards manufacturing phononic metamaterials, DNA origami structures and topolectric circuits that can realize a wide range of static and dynamic target functionalities.


title: Preparation contextuality as an essential feature underlying quantum communication advantage

abstract: The study of ontology (hidden variables) provides for a vital ground on which significant non-classical features of quantum theory are revealed. One such non-classical ontic-feature is preparation contextuality (PC) and advantage in oblivious communication tasks is one of its operational signatures. This article primarily pursues the ontic-feature underlying quantum advantage in communication complexity (CC). We construct oblivious communication tasks tailored to given CC problems. We upper-bound the classical success probability of these oblivious communication tasks, obtaining preparation non-contextual inequalities. We use the very states and measurements responsible for advantage in CC problems along with the orthogonal mixtures of these states to orchestrate an advantageous protocol for the associated oblivious communication tasks and the violation of the associated inequalities, thereby unveiling PC. To showcase the vitality of our results, we find a criterion for unbounded violation of these inequalities and demonstrate the same for two widely studied CC problems.


title: Cosmic Visions Dark Energy: Small Projects Portfolio

abstract: Understanding cosmic acceleration is one of the key science drivers for astrophysics and high-energy physics in the coming decade (2014 P5 Report). With the Large Synoptic Survey Telescope (LSST) and the Dark Energy Spectroscopic Instrument (DESI) and other new facilities beginning operations soon, we are entering an exciting phase during which we expect an order of magnitude improvement in constraints on dark energy and the physics of the accelerating Universe. This is a key moment for a matching Small Projects portfolio that can (1) greatly enhance the science reach of these flagship projects, (2) have immediate scientific impact, and (3) lay the groundwork for the next stages of the Cosmic Frontier Dark Energy program. In this White Paper, we outline a balanced portfolio that can accomplish these goals through a combination of observational, experimental, and theory and simulation efforts.


title: Measuring frequency fluctuations in nonlinear nanomechanical resonators

abstract: Advances in nanomechanics within recent years have demonstrated an always expanding range of devices, from top-down structures to appealing bottom-up MoS$_2$ and graphene membranes, used for both sensing and component-oriented applications. One of the main concerns in all of these devices is frequency noise, which ultimately limits their applicability. This issue has attracted a lot of attention recently, and the origin of this noise remains elusive up to date. In this Letter we present a very simple technique to measure frequency noise in nonlinear mechanical devices, based on the presence of bistability. It is illustrated on silicon-nitride high-stress doubly-clamped beams, in a cryogenic environment. We report on the same $T/f$ dependence of the frequency noise power spectra as reported in the literature. But we also find unexpected {\it damping fluctuations}, amplified in the vicinity of the bifurcation points; this effect is clearly distinct from already reported nonlinear dephasing, and poses a fundamental limit on the measurement of bifurcation frequencies. The technique is further applied to the measurement of frequency noise as a function of mode number, within the same device. The relative frequency noise for the fundamental flexure $\delta f/f_0$ lies in the range $0.5 - 0.01~$ppm (consistent with literature for cryogenic MHz devices), and decreases with mode number in the range studied. The technique can be applied to {\it any types} of nano-mechanical structures, enabling progresses towards the understanding of intrinsic sources of noise in these devices.


title: A 3D Game Theoretical Framework for the Evaluation of Unmanned Aircraft Systems Airspace Integration Concepts

abstract: Predicting the outcomes of integrating Unmanned Aerial Systems (UAS) into the National Airspace System (NAS) is a complex problem which is required to be addressed by simulation studies before allowing the routine access of UAS into the NAS. This paper focuses on providing a 3-dimensional (3D) simulation framework using a game theoretical methodology to evaluate integration concepts using scenarios where manned and unmanned air vehicles co-exist. In the proposed method, human pilot interactive decision making process is incorporated into airspace models which can fill the gap in the literature where the pilot behavior is generally assumed to be known a priori. The proposed human pilot behavior is modeled using dynamic level-k reasoning concept and approximate reinforcement learning. The level-k reasoning concept is a notion in game theory and is based on the assumption that humans have various levels of decision making. In the conventional "static" approach, each agent makes assumptions about his or her opponents and chooses his or her actions accordingly. On the other hand, in the dynamic level-k reasoning, agents can update their beliefs about their opponents and revise their level-k rule. In this study, Neural Fitted Q Iteration, which is an approximate reinforcement learning method, is used to model time-extended decisions of pilots with 3D maneuvers. An analysis of UAS integration is conducted using an example 3D scenario in the presence of manned aircraft and fully autonomous UAS equipped with sense and avoid algorithms.


title: Leibniz algebras as non-associative algebras

abstract: In this paper we define the basic concepts for left or right Leibniz algebras and prove some of the main results. Our proofs are often variations of the known proofs and several results seem to be new.


title: Thermoelectricity in correlated narrow-gap semiconductors

abstract: We review many-body effects, their microscopic origin, as well as their impact onto thermoelectricity in correlated narrow-gap semiconductors. Members of this class---such as FeSi and FeSb$_2$---display an unusual temperature dependence in various observables: insulating with large thermopowers at low temperatures, they turn bad metals at temperatures much smaller than the size of their gaps. This insulator-to-metal crossover is accompanied by spectral weight-transfers over large energies in the optical conductivity and by a gradual transition from activated to Curie-Weiss-like behaviour in the magnetic susceptibility. We show a retrospective of the understanding of these phenomena, discuss the relation to heavy-fermion Kondo insulators---such as Ce$_3$Bi$_4$Pt$_3$ for which we present new results---and propose a general classification of paramagnetic insulators. From the latter FeSi emerges as an orbital-selective Kondo insulator. Focussing on intermetallics such as silicides, antimonides, skutterudites, and Heusler compounds we showcase successes and challenges for the realistic simulation of transport properties in the presence of electronic correlations. Further, we advert to new avenues in which electronic correlations may contribute to the improvement of thermoelectric performance.


title: Geometric structures and representations of discrete groups

abstract: We describe recent links between two topics: geometric structures on manifolds in the sense of Ehresmann and Thurston, and dynamics "at infinity" for representations of discrete groups into Lie groups.


title: 007: Democratically Finding The Cause of Packet Drops

abstract: Network failures continue to plague datacenter operators as their symptoms may not have direct correlation with where or why they occur. We introduce 007, a lightweight, always-on diagnosis application that can find problematic links and also pinpoint problems for each TCP connection. 007 is completely contained within the end host. During its two month deployment in a tier-1 datacenter, it detected every problem found by previously deployed monitoring tools while also finding the sources of other problems previously undetected.