title: On the mass and luminosity functions of tidal disruption flares: rate suppression due to black hole event horizons

abstract: The tidal disruption of a star by a massive black hole is expected to yield a luminous flare of thermal emission. About two dozen of these stellar tidal disruption flares (TDFs) may have been detected in optical transient surveys. However, explaining the observed properties of these events within the tidal disruption paradigm is not yet possible. This theoretical ambiguity has led some authors to suggest that optical TDFs are due to a different process, such as a nuclear supernova or accretion disk instabilities. Here we present a test of a fundamental prediction of the tidal disruption event scenario: a suppression of the flare rate due to the direct capture of stars by the black hole. Using a recently compiled sample of candidate TDFs with black hole mass measurements, plus a careful treatment of selection effects in this flux-limited sample, we confirm that the dearth of observed TDFs from high-mass black holes is statistically significant. All the TDF impostor models we consider fail to explain the observed mass function; the only scenario that fits the data is a suppression of the rate due to direct captures. We find that this suppression can explain the low volumetric rate of the luminous TDF candidate ASASSN-15lh, thus supporting the hypothesis that this flare belongs to the TDF family. Our work is the first to present the optical TDF luminosity function. A steep power law is required to explain the observed rest-frame g-band luminosity, $dN/dL_{g} \propto L_{g}^{-2.5}$. The mean event rate of the flares in our sample is about $1 \times10^{-4}$ per galaxy per year, consistent with the theoretically expected tidal disruption rate.


title: The Tutte q-Polynomial

abstract: $q$-Matroids are defined on complemented modular support lattices. Minors of length 2 are of four types as in a "classical" matroid. Tutte polynomials $\tau(x,y)$ of matroids are calculated either by recursion over deletion/contraction of single elements, by an enumeration of bases with respect to internal/external activities, or by substitution $x \to (x-1),\; y \to (y-1)$ in their rank generating functions $\rho(x,y)$. The $q$-analogue of the passage from a Tutte polynomial to its corresponding RGF is straight-forward, but the analogue of the reverse process $x \to (x-1),\; y \to (y-1)$ is more delicate. For matroids $M(S)$ on a set $S$, and relative to any linear order on the points, the concept of internal/external activity of a point relative to a basis gives rise to a partition of the underlying Boolean algebra $B(S)$ into a set of "prime-free" (or "structureless") minors, such minors being direct sums of loops and isthmi (coloops), with one such prime-free minor for each basis. What usually goes unnoticed is that each prime-free minor has a unique clopen flat. The latter property carries over to $q$-matroids, but each prime-free minor will contain many bases. So internal and external activity in $q$-matroids must be defined not for points relative to bases, but rather for coverings in the underlying complemented modular lattice. Following lattice paths from arbitrary subspaces $A$ along active coverings (downward for internally active, upward for externally active) will lead to the unique clopen subspace in the prime-free minor containing the subspace $A$. There are a number of interesting questions concerning $q$-matroids that remain unsolved.


title: Tuning Perpendicular Magnetic Anisotropy by Oxygen Octahedral Rotations in (La$_{1-x}$Sr$_{x}$MnO$_{3}$)/(SrIrO$_{3}$) Superlattices

abstract: Perpendicular magnetic anisotropy (PMA) plays a critical role in the development of spintronics, thereby demanding new strategies to control PMA. Here we demonstrate a conceptually new type of interface induced PMA that is controlled by oxygen octahedral rotation. In superlattices comprised of La$_{1-x}$Sr$_{x}$MnO$_{3}$ and SrIrO$_{3}$, we find that all superlattices (0$\leq$x$\leq$1) exhibit ferromagnetism despite the fact that La$_{1-x}$Sr$_{x}$MnO$_{3}$ is antiferromagnetic for x$>$0.5. PMA as high as 4$\times$10$^6$ erg/cm$^3$ is observed by increasing x and attributed to a decrease of oxygen octahedral rotation at interfaces. We also demonstrate that oxygen octahedral deformation cannot explain the trend in PMA. These results reveal a new degree of freedom to control PMA, enabling discovery of emergent magnetic textures and topological phenomena.


title: Galaxy properties from J-PAS narrow-band photometry

abstract: We study the consistency of the physical properties of galaxies retrieved from SED-fitting as a function of spectral resolution and signal-to-noise ratio (SNR). Using a selection of physically motivated star formation histories, we set up a control sample of mock galaxy spectra representing observations of the local universe in high-resolution spectroscopy, and in 56 narrow-band and 5 broad-band photometry. We fit the SEDs at these spectral resolutions and compute their corresponding the stellar mass, the mass- and luminosity-weighted age and metallicity, and the dust extinction. We study the biases, correlations, and degeneracies affecting the retrieved parameters and explore the r\^ole of the spectral resolution and the SNR in regulating these degeneracies. We find that narrow-band photometry and spectroscopy yield similar trends in the physical properties derived, the former being considerably more precise. Using a galaxy sample from the SDSS, we compare more realistically the results obtained from high-resolution and narrow-band SEDs (synthesized from the same SDSS spectra) following the same spectral fitting procedures. We use results from the literature as a benchmark to our spectroscopic estimates and show that the prior PDFs, commonly adopted in parametric methods, may introduce biases not accounted for in a Bayesian framework. We conclude that narrow-band photometry yields the same trend in the age-metallicity relation in the literature, provided it is affected by the same biases as spectroscopy; albeit the precision achieved with the latter is generally twice as large as with the narrow-band, at SNR values typical of the different kinds of data.


title: Optimal design for high-throughput screening via false discovery rate control

abstract: High-throughput screening (HTS) is a large-scale hierarchical process in which a large number of chemicals are tested in multiple stages. Conventional statistical analyses of HTS studies often suffer from high testing error rates and soaring costs in large-scale settings. This article develops new methodologies for false discovery rate control and optimal design in HTS studies. We propose a two-stage procedure that determines the optimal numbers of replicates at different screening stages while simultaneously controlling the false discovery rate in the confirmatory stage subject to a constraint on the total budget. The merits of the proposed methods are illustrated using both simulated and real data. We show that the proposed screening procedure effectively controls the error rate and the design leads to improved detection power. This is achieved at the expense of a limited budget.


title: Division subspaces and integrable kernels

abstract: In this note we prove that the reproducing kernel of a Hilbert space satisfying the division property has integrable form, is locally of trace class, and the Hilbert space itself is a Hilbert space of holomorphic functions.


title: Uncertainty relations and joint numerical ranges

abstract: Noncommutativity lies in the heart of quantum theory and provides rich set of interesting questions in physics and mathematics. In this work I present some of them through the concept of Joint Numerical Range (JNR) - the set of simultaneously attainable averages of quantum observables, which in general need not commute. Full description of JNR geometry quickly becomes complex as the dimensionality of quantum system grows. In this thesis I present already existing results in the theory of JNR, generalizations of the formalism and applications of JNR in the theory of phase transitions and uncertainty relations.


title: Algebraic operads up to homotopy

abstract: This paper deals with the homotopy theory of differential graded operads. We endow the Koszul dual category of curved conilpotent cooperads, where the notion of quasi-isomorphism barely makes sense, with a model category structure Quillen equivalent to that of operads. This allows us to describe the homotopy properties of differential graded operads in a simpler and richer way, using obstruction methods.


title: Control strategy to limit duty cycle impact of earthquakes on the LIGO gravitational-wave detectors

abstract: Advanced gravitational-wave detectors such as the Laser Interferometer Gravitational-Wave Observatories (LIGO) require an unprecedented level of isolation from the ground. When in operation, they are expected to observe changes in the space-time continuum of less than one thousandth of the diameter of a proton. Strong teleseismic events like earthquakes disrupt the proper functioning of the detectors, and result in a loss of data until the detectors can be returned to their operating states. An earthquake early-warning system, as well as a prediction model have been developed to help understand the impact of earthquakes on LIGO. This paper describes a control strategy to use this early-warning system to reduce the LIGO downtime by 30%. It also presents a plan to implement this new earthquake configuration in the LIGO automation system.


title: Individual Recognition in Schizophrenia using Deep Learning Methods with Random Forest and Voting Classifiers: Insights from Resting State EEG Streams

abstract: Recently, there has been a growing interest in monitoring brain activity for individual recognition system. So far these works are mainly focussing on single channel data or fragment data collected by some advanced brain monitoring modalities. In this study we propose new individual recognition schemes based on spatio-temporal resting state Electroencephalography (EEG) data. Besides, instead of using features derived from artificially-designed procedures, modified deep learning architectures which aim to automatically extract an individual's unique features are developed to conduct classification. Our designed deep learning frameworks are proved of a small but consistent advantage of replacing the $softmax$ layer with Random Forest. Additionally, a voting layer is added at the top of designed neural networks in order to tackle the classification problem arisen from EEG streams. Lastly, various experiments are implemented to evaluate the performance of the designed deep learning architectures; Results indicate that the proposed EEG-based individual recognition scheme yields a high degree of classification accuracy: $81.6\%$ for characteristics in high risk (CHR) individuals, $96.7\%$ for clinically stable first episode patients with schizophrenia (FES) and $99.2\%$ for healthy controls (HC).


title: Recovering Dense Tissue Multispectral Signal from in vivo RGB Images

abstract: Hyperspectral/multispectral imaging (HSI/MSI) contains rich information clinical applications, such as 1) narrow band imaging for vascular visualisation; 2) oxygen saturation for intraoperative perfusion monitoring and clinical decision making [1]; 3) tissue classification and identification of pathology [2]. The current systems which provide pixel-level HSI/MSI signal can be generally divided into two types: spatial scanning and spectral scanning. However, the trade-off between spatial/spectral resolution, the acquisition time, and the hardware complexity hampers implementation in real-world applications, especially intra-operatively. Acquiring high resolution images in real-time is important for HSI/MSI in intra-operative imaging, to alleviate the side effect caused by breathing, heartbeat, and other sources of motion. Therefore, we developed an algorithm to recover a pixel-level MSI stack using only the captured snapshot RGB images from a normal camera. We refer to this technique as "super-spectral-resolution". The proposed method enables recovery of pixel-level-dense MSI signals with 24 spectral bands at ~11 frames per second (FPS) on a GPU. Multispectral data captured from porcine bowel and sheep/rabbit uteri in vivo has been used for training, and the algorithm has been validated using unseen in vivo animal experiments.


title: Machine Learning in Appearance-based Robot Self-localization

abstract: An appearance-based robot self-localization problem is considered in the machine learning framework. The appearance space is composed of all possible images, which can be captured by a robot's visual system under all robot localizations. Using recent manifold learning and deep learning techniques, we propose a new geometrically motivated solution based on training data consisting of a finite set of images captured in known locations of the robot. The solution includes estimation of the robot localization mapping from the appearance space to the robot localization space, as well as estimation of the inverse mapping for modeling visual image features. The latter allows solving the robot localization problem as the Kalman filtering problem.


title: Place recognition: An Overview of Vision Perspective

abstract: Place recognition is one of the most fundamental topics in computer vision and robotics communities, where the task is to accurately and efficiently recognize the location of a given query image. Despite years of wisdom accumulated in this field, place recognition still remains an open problem due to the various ways in which the appearance of real-world places may differ. This paper presents an overview of the place recognition literature. Since condition invariant and viewpoint invariant features are essential factors to long-term robust visual place recognition system, We start with traditional image description methodology developed in the past, which exploit techniques from image retrieval field. Recently, the rapid advances of related fields such as object detection and image classification have inspired a new technique to improve visual place recognition system, i.e., convolutional neural networks (CNNs). Thus we then introduce recent progress of visual place recognition system based on CNNs to automatically learn better image representations for places. Eventually, we close with discussions and future work of place recognition.


title: Proceedings of the 2017 AdKDD & TargetAd Workshop

abstract: Proceedings of the 2017 AdKDD and TargetAd Workshop held in conjunction with the 23rd ACM SIGKDD Conference on Knowledge Discovery and Data Mining Halifax, Nova Scotia, Canada.


title: Life in the coffee-ring: how evaporation-driven density gradients dictate the outcome of inter-bacterial competition

abstract: When a drop dries, it often leaves a ring-shaped stain through a ubiquitous phenomenon known as the coffee-ring effect. This also occurs when the liquid contains suspended microbes; evaporation leaves cells at higher concentrations in the ring than the drop interior. Using biofilm experiments and cellular automata simulations, we show that the physical structure created by the coffee-ring effect can dramatically alter the outcome of inter-bacterial competition. We experimentally study this effect using two strains of $\textit{Vibrio cholerae}$ that compete using a contact-dependent killing mechanism termed the type 6 secretion system. By creating a heterogeneous density profile, the coffee-ring effect changes the outcome of competition: the bacterial strain that wins in the low-density interior loses in the far denser ring. Through simulations parameterized with experimentally-determined density profiles, we recapitulate our experimental findings. We examine the role of a $\textit{V. cholerae}$ strain's frequency, its relative efficacy at killing competitors, and the initial concentration of cells in the droplet in determining the outcome of competition. By scaling from individual cellular interactions to overall changes in strain frequency, our work demonstrates how the coffee-ring effect plays a powerful role in structuring microbial communities, indirectly driving ecological changes in community composition.


title: Stacco: Differentially Analyzing Side-Channel Traces for Detecting SSL/TLS Vulnerabilities in Secure Enclaves

abstract: Intel Software Guard Extension (SGX) offers software applications enclave to protect their confidentiality and integrity from malicious operating systems. The SSL/TLS protocol, which is the de facto standard for protecting transport-layer network communications, has been broadly deployed for a secure communication channel. However, in this paper, we show that the marriage between SGX and SSL may not be smooth sailing. Particularly, we consider a category of side-channel attacks against SSL/TLS implementations in secure enclaves, which we call the control-flow inference attacks. In these attacks, the malicious operating system kernel may perform a powerful man-in-the-kernel attack to collect execution traces of the enclave programs at page, cacheline, or branch level, while positioning itself in the middle of the two communicating parties. At the center of our work is a differential analysis framework, dubbed Stacco, to dynamically analyze the SSL/TLS implementations and detect vulnerabilities that can be exploited as decryption oracles. Surprisingly, we found exploitable vulnerabilities in the latest versions of all the SSL/TLS libraries we have examined. To validate the detected vulnerabilities, we developed a man-in-the-kernel adversary to demonstrate Bleichenbacher attacks against the latest OpenSSL library running in the SGX enclave (with the help of Graphene) and completely broke the PreMasterSecret encrypted by a 4096-bit RSA public key with only 57286 queries. We also conducted CBC padding oracle attacks against the latest GnuTLS running in Graphene-SGX and an open-source SGX-implementation of mbedTLS (i.e., mbedTLS-SGX) that runs directly inside the enclave, and showed that it only needs 48388 and 25717 queries, respectively, to break one block of AES ciphertext. Empirical evaluation suggests these man-in-the-kernel attacks can be completed within 1 or 2 hours.


title: Quantum spaces, central extensions of Lie groups and related quantum field theories

abstract: Quantum spaces with $\frak{su}(2)$ noncommutativity can be modelled by using a family of $SO(3)$-equivariant differential $^*$-representations. The quantization maps are determined from the combination of the Wigner theorem for $SU(2)$ with the polar decomposition of the quantized plane waves. A tracial star-product, equivalent to the Kontsevich product for the Poisson manifold dual to $\mathfrak{su}(2)$ is obtained from a subfamily of differential $^*$-representations. Noncommutative (scalar) field theories free from UV/IR mixing and whose commutative limit coincides with the usual $\phi^4$ theory on $\mathbb{R}^3$ are presented. A generalization of the construction to semi-simple possibly non simply connected Lie groups based on their central extensions by suitable abelian Lie groups is discussed.


title: Stability of partially locked states in the Kuramoto model through Landau damping with Sobolev regularity

abstract: The Kuramoto model is a mean-field model for the synchronisation behaviour of oscillators, which exhibits Landau damping. In a recent work, the nonlinear stability of a class of spatially inhomogeneous stationary states was shown under the assumption of analytic regularity. This paper proves the nonlinear Landau damping under the assumption of Sobolev regularity. The weaker regularity required the construction of a different more robust bootstrap argument, which focuses on the nonlinear Volterra equation of the order parameter.


title: Vibrationally resolved electronic spectra including vibrational pre-excitation: Theory and application to VIPER spectroscopy

abstract: Vibrationally resolved electronic absorption spectra including the effect of vibrational pre-excitation are computed in order to interpret and predict vibronic transitions that are probed in the Vibrationally Promoted Electronic Resonance (VIPER) experiment [L. J. G. W. van Wilderen et al., Angew. Chem. Int. Ed. 53, 2667 (2014)]. To this end, we employ time-independent and time-dependent methods based on the evaluation of Franck-Condon overlap integrals and Fourier transformation of time-domain wavepacket autocorrelation functions, respectively. The time-independent approach uses a generalized version of the FCclasses method [F. Santoro et al., J. Chem. Phys. 126, 084509 (2007)]. In the time-dependent approach, autocorrelation functions are obtained by wavepacket propagation and by evaluation of analytic expressions, within the harmonic approximation including Duschinsky rotation effects. For several medium-sized polyatomic systems, it is shown that selective pre-excitation of particular vibrational modes leads to a red-shift of the low-frequency edge of the electronic absorption spectrum, which is a prerequisite for the VIPER experiment. This effect is typically most pronounced upon excitation of ring distortion modes within an aromatic pi-system. Theoretical predictions as to which modes show the strongest VIPER effect are found to be in excellent agreement with experiment.


title: A Grazing Gaussian Beam

abstract: We consider Friedlander's wave equation in two space dimensions in the half-space x > 0 with the boundary condition u(x,y,t)=0 when x=0. For a Gaussian beam w(x,y,t;k) concentrated on a ray path that is tangent to x=0 at (x,y,t)=(0,0,0) we calculate the "reflected" wave z(x,y,t;k) in t > 0 such that w(x,y,t;k)+z(x,y,t;k) satisfies Friedlander's wave equation and vanishes on x=0. These computations are done to leading order in k on the ray path. The interaction of beams with boundaries has been studied for non-tangential beams and for beams gliding along the boundary. We find that the amplitude of the solution on the central ray for large k after leaving the boundary is very nearly one half of that of the incoming beam.


title: Round Compression for Parallel Matching Algorithms

abstract: For over a decade now we have been witnessing the success of {\em massive parallel computation} (MPC) frameworks, such as MapReduce, Hadoop, Dryad, or Spark. One of the reasons for their success is the fact that these frameworks are able to accurately capture the nature of large-scale computation. In particular, compared to the classic distributed algorithms or PRAM models, these frameworks allow for much more local computation. The fundamental question that arises in this context is though: can we leverage this additional power to obtain even faster parallel algorithms? A prominent example here is the {\em maximum matching} problem---one of the most classic graph problems. It is well known that in the PRAM model one can compute a 2-approximate maximum matching in $O(\log{n})$ rounds. However, the exact complexity of this problem in the MPC framework is still far from understood. Lattanzi et al. showed that if each machine has $n^{1+\Omega(1)}$ memory, this problem can also be solved $2$-approximately in a constant number of rounds. These techniques, as well as the approaches developed in the follow up work, seem though to get stuck in a fundamental way at roughly $O(\log{n})$ rounds once we enter the near-linear memory regime. It is thus entirely possible that in this regime, which captures in particular the case of sparse graph computations, the best MPC round complexity matches what one can already get in the PRAM model, without the need to take advantage of the extra local computation power. In this paper, we finally refute that perplexing possibility. That is, we break the above $O(\log n)$ round complexity bound even in the case of {\em slightly sublinear} memory per machine. In fact, our improvement here is {\em almost exponential}: we are able to deliver a $(2+\epsilon)$-approximation to maximum matching, for any fixed constant $\epsilon>0$, in $O((\log \log n)^2)$ rounds.


title: A MacDonald formula for zeta functions of varieties over finite fields

abstract: We provide a formula for the generating series of the zeta function $Z(X,t)$ of symmetric powers $Sym^n X$ of varieties over finite fields. This realizes $Z(X,t)$ as an exponentiable motivic measure whose associated Kapranov motivic zeta function takes values in $W(R)$ the big Witt ring of $R=W(\mathbb{Z})$. We apply our formula to compute $Z(Sym^n X,t)$ in a number of explicit cases. Moreover, we show that all $\lambda$-ring motivic measures have zeta functions which are exponentiable. In this setting, the formula for $Z(X,t)$ takes the form of a MacDonald formula for the zeta function.


title: Topological order in the Haldane model with spin-spin on-site interactions

abstract: Ultracold atom experiments allow the study of topological insulators, such as the noninteracting Haldane model. In this work we study a generalization of the Haldane model with spin-spin on-site interactions that can be implemented on such experiments. We focus on measuring the winding number, a topological invariant, of the ground state, which we compute using a mean-field calculation that effectively captures long range correlations and a matrix product state computation in a lattice with 64 sites. Our main result is that we show how the topological phases present in the noninteracting model survive until the interactions are comparable to the kinetic energy. We also demonstrate the accuracy of our mean-field approach in efficiently capturing long-range correlations. Based on state-of-the-art ultracold atom experiments, we propose an implementation of our model that can give information about the topological phases.


title: Magnetic pseudo-fields in a rotating electron-nuclear spin system

abstract: A precessing spin observed in a rotating frame of reference appears frequency-shifted, an effect analogous to the precession of a Foucault pendulum observed on the rotating Earth. This frequency shift can be understood as arising from a magnetic pseudo-field in the rotating frame that nevertheless has physically significant consequences, such as the Barnett effect. Detecting these pseudo-fields is experimentally challenging, as a rotating-frame sensor is required. Previous work has realised classical rotating-frame detectors. Here we use quantum sensors, nitrogen-vacancy (NV) centres, in a rapidly rotating diamond to detect pseudo-fields in the rotating frame. While conventional magnetic fields induce precession at a rate proportional to the gyromagnetic ratio, rotation shifts the precession of all spins equally, and thus primarily affect nearby $^{13}$C nuclear spins. We are thus able to explore these effects via quantum sensing in a rapidly rotating frame, and define a new approach to quantum control using rotationally-induced nuclear spin-selective magnetic fields. This work provides an integral step towards realising precision rotation sensing and quantum spin gyroscopes.


title: Discrete Bethe-Sommerfeld Conjecture

abstract: In this paper, we prove a discrete version of the Bethe-Sommerfeld conjecture. Namely, we show that the spectra of multi-dimensional discrete periodic Schr\"odinger operators on $\mathbb{Z}^d$ lattice with sufficiently small potentials contain at most two intervals. Moreover, the spectrum is a single interval, provided one of the periods is odd, and can have a gap whenever all periods are even.


title: A class of rotating hairy black holes in arbitrary dimensions

abstract: A class of exact rotating black hole solutions of gravity nonminimally coupled to a self-interacting scalar field in arbitrary dimensions is presented. These spacetimes are asymptotically locally anti-de Sitter manifolds and have a Ricci-flat event horizon hiding a curvature singularity at the origin. The scalar field is real and regular everywhere and its effective mass, coming from the nonminimal coupling with the scalar curvature, saturates the Breitenlohner-Freedman bound for the corresponding spacetime dimension. The rotating black hole is obtained by applying an improper coordinate transformation to the static one. Although both spacetimes are locally equivalent, they are globally different, as it is confirmed by the nonvanishing angular momentum of the rotating black hole. It is found that the mass is bounded from below by the angular momentum in agreement with the existence of an event horizon. The thermodynamical analysis is carried out in the grand canonical ensemble. The first law is satisfied and a Smarr formula is exhibited. The thermodynamical local stability of the rotating hairy black holes is established from their Gibbs free energy. However, the global stability analysis establishes that the vacuum spacetime is always preferred over the hairy black hole. Thus, the hairy black hole is likely to decay into the vacuum one for any temperature.


title: Anthropics of Aluminum-26 Decay and Biological Homochirality

abstract: Results of recent experiment reinstate feasibility to the hypothesis that biomolecular homochirality originates from beta decay. Coupled with hints that this process occurred extraterrestrially suggests aluminum-26 as the most likely source. If true, then its appropriateness is highly dependent on the half-life and energy of this decay. Demanding that this mechanism hold places new constraints on the anthropically allowed range for multiple parameters, including the electron mass, difference between up and down quark masses, the fine structure constant, and the electroweak scale. These new constraints on particle masses are tighter than those previously found. However, one edge of the allowed region is nearly degenerate with an existing bound, which, using what is termed here as `the principle of noncoincident peril', is argued to be a strong indicator that the fine structure constant must be an environmental parameter in the multiverse.


title: Coefficient groups inducing nonbranched optimal transport

abstract: In this work we consider an optimal transport problem with coefficients in a normed Abelian group $G$, and extract a purely intrinsic condition on $G$ that guarantees that the optimal transport (or the corresponding minimum filling) is not branching. The condition turns out to be equivalent to the nonbranching of minimum fillings in geodesic metric spaces. We completely characterize finitely generated normed groups and finite-dimensional normed vector spaces of coefficients that induce nonbranching optimal transport plans. We also provide a complete classification of normed groups for which the optimal transport plans, besides being nonbranching, have acyclic support. This seems to initiate a new geometric classifications of certain normed groups. In the nonbranching case we also provide a global version of calibration, i.e. a generalization of Monge-Kantorovich duality.


title: Bounded pregeometries and pairs of fields

abstract: A definable set in a pair (K, k) of algebraically closed fields is co-analyzable relative to the subfield k of the pair if and only if it is almost internal to k. To prove this and some related results for tame pairs of real closed fields we introduce a certain kind of bounded pregeometry for such pairs.


title: Robust Estimation from Multiple Graphs under Gross Error Contamination

abstract: Estimation of graph parameters based on a collection of graphs is essential for a wide range of graph inference tasks. In practice, weighted graphs are generally observed with edge contamination. We consider a weighted latent position graph model contaminated via an edge weight gross error model and propose an estimation methodology based on robust Lq estimation followed by low-rank adjacency spectral decomposition. We demonstrate that, under appropriate conditions, our estimator both maintains Lq robustness and wins the bias-variance tradeoff by exploiting low-rank graph structure. We illustrate the improvement offered by our estimator via both simulations and a human connectome data experiment.


title: The Spectral Position of Neumann domains on the torus

abstract: Neumann domains of Laplacian eigenfunctions form a natural counterpart of nodal domains. The restriction of an eigenfunction to one of its nodal domains is the first Dirichlet eigenfunction of that domain. This simple observation is fundamental in many works on nodal domains. We consider a similar property for Neumann domains. Namely, given a Laplacian eigenfunction $f$ and its Neumann domain $\Omega$, what is the position of $\ .f\ |_{\Omega}$ in the Neumann spectrum of $\Omega$? The current paper treats this spectral position problem on the two-dimensional torus. We fully solve it for separable eigenfunctions on the torus and complement our analytic solution with numerics for random waves on the torus. These results answer questions from [37, 8] and arouse new ones.


title: Global stability properties of renewal epidemic models

abstract: We investigate the global dynamics of a general Kermack-McKendrick-type epidemic model formulated in terms of a system of renewal equations. Specifically, we consider a renewal model for which both the force of infection and the infected removal rates are arbitrary functions of the infection age, $\tau$, and use the direct Lyapunov method to establish the global asymptotic stability of the equilibrium solutions. In particular, we show that the basic reproduction number, $R_0$, represents a sharp threshold parameter such that for $R_0\leq 1$, the infection-free equilibrium is globally asymptotically stable; whereas the endemic equilibrium becomes globally asymptotically stable when $R_0 > 1$, i.e. when it exists.


title: Detecting Policy Preferences and Dynamics in the UN General Debate with Neural Word Embeddings

abstract: Foreign policy analysis has been struggling to find ways to measure policy preferences and paradigm shifts in international political systems. This paper presents a novel, potential solution to this challenge, through the application of a neural word embedding (Word2vec) model on a dataset featuring speeches by heads of state or government in the United Nations General Debate. The paper provides three key contributions based on the output of the Word2vec model. First, it presents a set of policy attention indices, synthesizing the semantic proximity of political speeches to specific policy themes. Second, it introduces country-specific semantic centrality indices, based on topological analyses of countries' semantic positions with respect to each other. Third, it tests the hypothesis that there exists a statistical relation between the semantic content of political speeches and UN voting behavior, falsifying it and suggesting that political speeches contain information of different nature then the one behind voting outcomes. The paper concludes with a discussion of the practical use of its results and consequences for foreign policy analysis, public accountability, and transparency.


title: Creatism: A deep-learning photographer capable of creating professional work

abstract: Machine-learning excels in many areas with well-defined goals. However, a clear goal is usually not available in art forms, such as photography. The success of a photograph is measured by its aesthetic value, a very subjective concept. This adds to the challenge for a machine learning approach. We introduce Creatism, a deep-learning system for artistic content creation. In our system, we break down aesthetics into multiple aspects, each can be learned individually from a shared dataset of professional examples. Each aspect corresponds to an image operation that can be optimized efficiently. A novel editing tool, dramatic mask, is introduced as one operation that improves dramatic lighting for a photo. Our training does not require a dataset with before/after image pairs, or any additional labels to indicate different aspects in aesthetics. Using our system, we mimic the workflow of a landscape photographer, from framing for the best composition to carrying out various post-processing operations. The environment for our virtual photographer is simulated by a collection of panorama images from Google Street View. We design a "Turing-test"-like experiment to objectively measure quality of its creations, where professional photographers rate a mixture of photographs from different sources blindly. Experiments show that a portion of our robot's creation can be confused with professional work.


title: Fault-Induced Dynamics of Oblivious Robots on a Line

abstract: The study of computing in presence of faulty robots in the Look-Compute-Move model has been the object of extensive investigation, typically with the goal of designing algorithms tolerant to as many faults as possible. In this paper, we initiate a new line of investigation on the presence of faults, focusing on a rather different issue. We are interested in understanding the dynamics of a group of robots when they execute an algorithm designed for a fault-free environment, in presence of some undetectable crashed robots. We start this investigation focusing on the classic point-convergence algorithm by Ando et al. for robots with limited visibility, in a simple setting (which already presents serious challenges): the robots operate fully synchronously on a line, and at most two of them are faulty. Interestingly, and perhaps surprisingly, the presence of faults induces the robots to perform some form of scattering, rather than point-convergence. In fact, we discover that they arrange themselves inside the segment delimited by the two faults in interleaved sequences of equidistant robots. The structure that they form has a hierarchical nature: robots organize themselves in groups where a group of some level converges to an equidistant distribution only after all groups of lower levels have converged. This is the first study on the unintended dynamics of oblivious robots induced by the presence of faults.


title: Numerical integration of blow-up problems on the basis of non-local transformations and differential constraints

abstract: Several new methods of numerical integration of Cauchy problems with blow-up solutions for nonlinear ordinary differential equations of the first- and second-order are described. Solutions of such problems have singularities whose positions are unknown a priori (the standard numerical methods for solving problems with blow-up solutions can lead to significant errors). The first proposed method is based on the transition to an equivalent system of equations by introducing a new independent variable chosen as the first derivative. The second method is based on introducing a new auxiliary non-local variable with the subsequent transformation to the Cauchy problem for the corresponding system of ODEs. The third method is based on adding to the original equation of a differential constraint, which is an auxiliary ODE connecting the given variables and a new variable. The proposed methods lead to problems whose solutions are represented in parametric form and do not have blowing-up singular points; therefore the transformed problems admit the application of standard fixed-step numerical methods. The efficiency of these methods is illustrated by solving a number of test problems that admit an exact analytical solution. It is shown that: (i) the methods based on non-local transformations of a special kind are more efficient than several other methods, (ii) among the proposed methods, the most general method is the method based on the differential constraints. Some examples of nonclassical blow-up problems are considered. Simple theoretical estimates are derived for the critical value of an independent variable. It is shown that the method based on a non-local transformation of the general form as well as the method based on the differential constraints admit generalizations to the $n$th-order ODEs and systems of coupled ODEs.


title: Unsupervised robust nonparametric learning of hidden community properties

abstract: We consider learning of fundamental properties of communities in large noisy networks, in the prototypical situation where the nodes or users are split into two classes according to a binary property, e.g., according to their opinions or preferences on a topic. For learning these properties, we propose a nonparametric, unsupervised, and scalable graph scan procedure that is, in addition, robust against a class of powerful adversaries. In our setup, one of the communities can fall under the influence of a knowledgeable adversarial leader, who knows the full network structure, has unlimited computational resources and can completely foresee our planned actions on the network. We prove strong consistency of our results in this setup with minimal assumptions. In particular, the learning procedure estimates the baseline activity of normal users asymptotically correctly with probability 1; the only assumption being the existence of a single implicit community of asymptotically negligible logarithmic size. We provide experiments on real and synthetic data to illustrate the performance of our method, including examples with adversaries.


title: Lengthening and Extending Binary Private Information Retrieval Codes

abstract: It was recently shown by Fazeli et al. that the storage overhead of a traditional $t$-server private information retrieval (PIR) protocol can be significantly reduced using the concept of a $t$-server PIR code. In this work, we show that a family of $t$-server PIR codes (with increasing dimensions and blocklengths) can be constructed from an existing $t$-server PIR code through lengthening by a single information symbol and code extension by at most $\bigl\lceil t/2\bigr\rceil$ code symbols. Furthermore, by extending a code construction notion from Steiner systems by Fazeli et al., we obtain a specific family of $t$-server PIR codes. Based on a code construction technique that lengthens and extends a $t$-server PIR code simultaneously, a basic algorithm to find good (i.e., small blocklength) $t$-server PIR codes is proposed. For the special case of $t=5$, we find provably optimal PIR codes for code dimensions $k\leq 6$, while for all $7\leq k\leq 32$ we find codes of smaller blocklength than the best known codes from the literature. Furthermore, in the case of $t = 8$, we also find better codes for $k = 5, 6, 11, 12$. Numerical results show that most of the best found $5$-server PIR codes can be constructed from the proposed family of codes connected to Steiner systems.


title: Conservation laws by symmetries and adjoint symmetries

abstract: Conservation laws are formulated for systems of differential equations by using symmetries and adjoint symmetries, and an application to systems of evolution equations is made, together with illustrative examples. The formulation does not require the existence of a Lagrangian for a given system, and the presented examples include computations of conserved densities for the heat equation, Burgers' equation and the Korteweg-de Vries equation.


title: Value Prediction Network

abstract: This paper proposes a novel deep reinforcement learning (RL) architecture, called Value Prediction Network (VPN), which integrates model-free and model-based RL methods into a single neural network. In contrast to typical model-based RL methods, VPN learns a dynamics model whose abstract states are trained to make option-conditional predictions of future values (discounted sum of rewards) rather than of future observations. Our experimental results show that VPN has several advantages over both model-free and model-based baselines in a stochastic environment where careful planning is required but building an accurate observation-prediction model is difficult. Furthermore, VPN outperforms Deep Q-Network (DQN) on several Atari games even with short-lookahead planning, demonstrating its potential as a new way of learning a good state representation.


title: Mean Reversion Trading with Sequential Deadlines and Transaction Costs

abstract: We study the optimal timing strategies for trading a mean-reverting price process with afinite deadline to enter and a separate finite deadline to exit the market. The price process is modeled by a diffusion with an affine drift that encapsulates a number of well-known models,including the Ornstein-Uhlenbeck (OU) model, Cox-Ingersoll-Ross (CIR) model, Jacobi model,and inhomogeneous geometric Brownian motion (IGBM) model.We analyze three types of trading strategies: (i) the long-short (long to open, short to close) strategy; (ii) the short-long(short to open, long to close) strategy, and (iii) the chooser strategy whereby the trader has the added flexibility to enter the market by taking either a long or short position, and subsequently close the position. For each strategy, we solve an optimal double stopping problem with sequential deadlines, and determine the optimal timing of trades. Our solution methodology utilizes the local time-space calculus of Peskir (2005) to derive nonlinear integral equations of Volterra-type that uniquely characterize the trading boundaries. Numerical implementation ofthe integral equations provides examples of the optimal trading boundaries.


title: The Engineering Development Array: A low frequency radio telescope utilising SKA precursor technology

abstract: We describe the design and performance of the Engineering Development Array (EDA), which is a low frequency radio telescope comprising 256 dual-polarisation dipole antennas working as a phased-array. The EDA was conceived of, developed, and deployed in just 18 months via re-use of Square Kilometre Array (SKA) precursor technology and expertise, specifically from the Murchison Widefield Array (MWA) radio telescope. Using drift scans and a model for the sky brightness temperature at low frequencies, we have derived the EDA's receiver temperature as a function of frequency. The EDA is shown to be sky-noise limited over most of the frequency range measured between 60 and 240 MHz. By using the EDA in interferometric mode with the MWA, we used calibrated visibilities to measure the absolute sensitivity of the array. The measured array sensitivity matches very well with a model based on the array layout and measured receiver temperature. The results demonstrate the practicality and feasibility of using MWA-style precursor technology for SKA-scale stations. The modular architecture of the EDA allows upgrades to the array to be rolled out in a staged approach. Future improvements to the EDA include replacing the second stage beamformer with a fully digital system, and to transition to using RF-over-fibre for the signal output from first stage beamformers.


title: Banking risk as an epidemiological model: an optimal control approach

abstract: The process of contagiousness spread modelling is well-known in epidemiology. However, the application of spread modelling to banking market is quite recent. In this work, we present a system of ordinary differential equations, simulating data from the largest European banks. Then, an optimal control problem is formulated in order to study the impact of a possible measure of the Central Bank in the economy. The proposed approach enables qualitative specifications of contagion in banking obtainment and an adequate analysis and prognosis within the financial sector development and macroeconomic as a whole. We show that our model describes well the reality of the largest European banks. Simulations were done using MATLAB and BOCOP optimal control solver, and the main results are taken for three distinct scenarios.


title: NO Need to Worry about Adversarial Examples in Object Detection in Autonomous Vehicles

abstract: It has been shown that most machine learning algorithms are susceptible to adversarial perturbations. Slightly perturbing an image in a carefully chosen direction in the image space may cause a trained neural network model to misclassify it. Recently, it was shown that physical adversarial examples exist: printing perturbed images then taking pictures of them would still result in misclassification. This raises security and safety concerns. However, these experiments ignore a crucial property of physical objects: the camera can view objects from different distances and at different angles. In this paper, we show experiments that suggest that current constructions of physical adversarial examples do not disrupt object detection from a moving platform. Instead, a trained neural network classifies most of the pictures taken from different distances and angles of a perturbed image correctly. We believe this is because the adversarial property of the perturbation is sensitive to the scale at which the perturbed picture is viewed, so (for example) an autonomous car will misclassify a stop sign only from a small range of distances. Our work raises an important question: can one construct examples that are adversarial for many or most viewing conditions? If so, the construction should offer very significant insights into the internal representation of patterns by deep networks. If not, there is a good prospect that adversarial examples can be reduced to a curiosity with little practical impact.


title: Deep Learning for Sensor-based Activity Recognition: A Survey

abstract: Sensor-based activity recognition seeks the profound high-level knowledge about human activities from multitudes of low-level sensor readings. Conventional pattern recognition approaches have made tremendous progress in the past years. However, those methods often heavily rely on heuristic hand-crafted feature extraction, which could hinder their generalization performance. Additionally, existing methods are undermined for unsupervised and incremental learning tasks. Recently, the recent advancement of deep learning makes it possible to perform automatic high-level feature extraction thus achieves promising performance in many areas. Since then, deep learning based methods have been widely adopted for the sensor-based activity recognition tasks. This paper surveys the recent advance of deep learning based sensor-based activity recognition. We summarize existing literature from three aspects: sensor modality, deep model, and application. We also present detailed insights on existing work and propose grand challenges for future research.


title: THz-frequency cavity magnon-phonon-polaritons in the strong coupling regime

abstract: We demonstrate the strong coupling of both magnons and phonons to terahertz (THz) frequency electromagnetic (EM) waves confined to a photonic crystal (PhC) cavity. Our cavity consists of a two-dimensional array of air-holes cut into a hybrid slab of ferroelectric lithium niobate (LiNbO$_3$) and erbium orthoferrite (ErFeO$_3$), a canted antiferromagnetic crystal. The phonons in LiNbO$_3$ and the magnons in ErFeO$_3$ are strongly coupled to the electric and magnetic field components of the confined EM wave, respectively. This leads to the formation of new cavity magnon-phonon-polariton modes, which we experimentally observe as a normal-mode splitting in the frequency spectrum and an avoided crossing in the temperature-frequency plot. The cavity also has a mode volume of $V=3.4\times10^{-3}\lambda^3\simeq0.5(\lambda/n)^3$ $\mu$m$^3$ and can achieve a Q-factor as high as 1000. These factors facilitate the pursuit of the fields of THz cavity spintronics and quantum electrodynamics.


title: Strongly screening electron capture rates of chromium isotopes in presupernova

abstract: Taking into account the effect of electron screening on the electron energy and electron capture threshold energy, by using the method of Shell-Model Monte Carlo and Random Phase Approximation theory, we investigate the strong electron screening capture rates of chromium isotopes according to the linear response theory screening model. The strong screening rates can decrease by about 40.43\% (e.g., for $^{60}$Cr at $T_9=3.44,Y_e=0.43$). Our conclusions may be helpful to the researches of supernova explosion and numerical simulation.


title: Proximally Guided Stochastic Subgradient Method for Nonsmooth, Nonconvex Problems

abstract: In this paper, we introduce a stochastic projected subgradient method for weakly convex (i.e., uniformly prox-regular) nonsmooth, nonconvex functions---a wide class of functions which includes the additive and convex composite classes. At a high-level, the method is an inexact proximal point iteration in which the strongly convex proximal subproblems are quickly solved with a specialized stochastic projected subgradient method. The primary contribution of this paper is a simple proof that the proposed algorithm converges at the same rate as the stochastic gradient method for smooth nonconvex problems. This result appears to be the first convergence rate analysis of a stochastic (or even deterministic) subgradient method for the class of weakly convex functions.


title: Phase correction for ALMA - Investigating water vapour radiometer scaling:The long-baseline science verification data case study

abstract: The Atacama Large millimetre/submillimetre Array (ALMA) makes use of water vapour radiometers (WVR), which monitor the atmospheric water vapour line at 183 GHz along the line of sight above each antenna to correct for phase delays introduced by the wet component of the troposphere. The application of WVR derived phase corrections improve the image quality and facilitate successful observations in weather conditions that were classically marginal or poor. We present work to indicate that a scaling factor applied to the WVR solutions can act to further improve the phase stability and image quality of ALMA data. We find reduced phase noise statistics for 62 out of 75 datasets from the long-baseline science verification campaign after a WVR scaling factor is applied. The improvement of phase noise translates to an expected coherence improvement in 39 datasets. When imaging the bandpass source, we find 33 of the 39 datasets show an improvement in the signal-to-noise ratio (S/N) between a few to ~30 percent. There are 23 datasets where the S/N of the science image is improved: 6 by <1%, 11 between 1 and 5%, and 6 above 5%. The higher frequencies studied (band 6 and band 7) are those most improved, specifically datasets with low precipitable water vapour (PWV), <1mm, where the dominance of the wet component is reduced. Although these improvements are not profound, phase stability improvements via the WVR scaling factor come into play for the higher frequency (>450 GHz) and long-baseline (>5km) observations. These inherently have poorer phase stability and are taken in low PWV (<1mm) conditions for which we find the scaling to be most effective. A promising explanation for the scaling factor is the mixing of dry and wet air components, although other origins are discussed. We have produced a python code to allow ALMA users to undertake WVR scaling tests and make improvements to their data.


title: Antimagic orientations of even regular graphs

abstract: A labeling of a digraph $D$ with $m$ arcs is a bijection from the set of arcs of $D$ to $\{1, \ldots, m\}$. A labeling of $D$ is antimagic if no two vertices in $D$ have the same vertex-sum, where the vertex-sum of a vertex $u\in V(D)$ for a labeling is the sum of labels of all arcs entering $u$ minus the sum of labels of all arcs leaving $u$. Motivated by the conjecture of Hartsfield and Ringel from 1990 on antimagic labelings of graphs, Hefetz, M\"utze, and Schwartz [On antimagic directed graphs, J Graph Theory 64 (2010) 219--232] initiated the study of antimagic labelings of digraphs, and conjectured that every connected graph admits an antimagic orientation, where an orientation $D$ of a graph $G$ is antimagic if $D$ has an antimagic labeling. It remained unknown whether every disjoint union of cycles admits an antimagic orientation. In this paper, we first answer this question in the positive by proving that every $2$-regular graph has an antimagic orientation. We then show that for any integer $d\ge2$, every connected, $2d$-regular graph has an antimagic orientation. Our technique is new.


title: Doping and Defect-Induced Germanene: A Superior Media for Sensing H2S, SO2, and CO2 gas molecules

abstract: First-principles calculations based on density functional theory (DFT) have been employed to investigate the structural, electronic, and gas-sensing properties of pure, defected, and doped germanene nanosheets. Our calculations have revealed that while a pristine germanene nanosheet adsorbs CO2 weakly, H2S moderately, and SO2 strongly, the introduction of vacancy defects increases the sensitivity significantly which is promising for future gas-sensing applications.Mulliken population analysis imparts that an appreciable amount of charge transfer occurs between gas molecules and a germanene nanosheet which supports our results for adsorption energies of the systems. The enhancement of the interactions between gas molecules and the germanene nanosheet has been further investigated by density of states.Projected density of states provides detailed insight of the gas molecules contribution in the gas-sensing system.Additionally, the influences of substituted dopant atoms such as B, N, and Al in the germanene nanosheet have also been considered to study the impact on its gas sensing ability. There was no significant improvement found in doped gas sensing capability of germanene over the vacancy defects, except for CO2 upon adsorption on N-doped germanene.


title: Polarity-tunable magnetic tunnel junctions based on ferromagnetism at oxide heterointerfaces

abstract: Complex oxide systems have attracted considerable attention because of their fascinating properties, including the magnetic ordering at the conducting interface between two band insulators, such as LaAlO3 (LAO) and SrTiO3 (STO). However, the manipulation of the spin degree of freedom at the LAO/STO heterointerface has remained elusive. Here, we have fabricated hybrid magnetic tunnel junctions consisting of Co and LAO/STO ferromagnets with the insertion of a Ti layer in between, which clearly exhibit magnetic switching and the tunnelling magnetoresistance (TMR) effect below 10 K. The magnitude and the of the TMR are strongly dependent on the direction of the rotational magnetic field parallel to the LAO/STO plane, which is attributed to a strong Rashba-type spin orbit coupling in the LAO/STO heterostructure. Our study provides a further support for the existence of the macroscopic ferromagnetism at LAO/STO heterointerfaces and opens a novel route to realize interfacial spintronics devices.


title: Association of Networked Flying Platforms with Small Cells for Network Centric 5G+ C-RAN

abstract: 5G+ systems expect enhancement in data rate and coverage area under limited power constraint. Such requirements can be fulfilled by the densification of small cells (SCs). However, a major challenge is the management of fronthaul links connected with an ultra dense network of SCs. A cost effective and scalable idea of using network flying platforms (NFPs) is employed here, where the NFPs are used as fronthaul hubs that connect the SCs to the core network. The association problem of NFPs and SCs is formulated considering a number of practical constraints such as backhaul data rate limit, maximum supported links and bandwidth by NFPs and quality of service requirement of the system. The network centric case of the system is considered that aims to maximize the number of associated SCs without any biasing, i.e., no preference for high priority SCs. Then, two new efficient greedy algorithms are designed to solve the presented association problem. Numerical results show a favorable performance of our proposed methods in comparison to exhaustive search.


title: Multisoliton solutions of the two-component Camassa-Holm system and their reductions

abstract: We develop a systematic procedure for constructing soliton solutions of an integrable two-component Camassa-Holm (CH2) system. The parametric representation of the multisoliton solutions is obtained by using a direct method combined with a reciprocal transformation. The properties of the solutions are then investigated in detail focusing mainly on the smooth one- and two-soliton solutions. The general $N$-soliton case is described shortly. Subsequently, we show that the CH2 system reduces to the CH equation and the two-component Hunter-Saxton (HS2) system by means of appropriate limiting procedures. The corresponding expressions of the multisoliton solutions are presented in parametric forms, reproducing the existing results for the reduced equations. Last, we discuss the reduction from the HS2 system to the HS equation.


title: Stretching p-wave molecules by transverse confinements

abstract: We revisit the confinement-induced p-wave resonance in quasi-one-dimensional (quasi-1D) atomic gases and study the induced molecules near resonance. We derive the reduced 1D interaction parameters and show that they can well predict the binding energy of shallow molecules in quasi-1D system. Importantly, these shallow molecules are found to be much more spatially extended compared to those in three dimensions (3D) without transverse confinement. Our results strongly indicate that a p-wave interacting atomic gas can be much more stable in quasi-1D near the induced p-wave resonance, where most weight of the molecule lies outside the short-range regime and thus the atom loss could be suppressed.


title: Detecting primordial gravitational waves with circular polarization of the redshifted 21 cm line: I. Formalism

abstract: We propose a new method to measure the tensor-to-scalar ratio $r$ using the circular polarization of the 21 cm radiation from the pre-reionization epoch. Our method relies on the splitting of the $F = 1$ hyperfine level of neutral hydrogen due to the quadrupole moment of the CMB. We show that unlike the Zeeman effect, where $M_{F}=\pm 1$ have opposite energy shifts, the CMB quadrupole shifts $M_{F}=\pm 1$ together relative to $M_{F}= 0$. This splitting leads to a small circular polarization of the emitted 21 cm radiation. In this paper (Paper I in a series on this effect), we present calculations on the microphysics behind this effect, accounting for all processes that affect the hyperfine transition. We conclude with an analytic formula for the circular polarization from the Dark Ages as a function of pre-reionization parameters and the value of the remote quadrupole of the CMB. We also calculate the splitting of the $F = 1$ hyperfine level due to other anisotropic radiation sources and show that they are not dominant. In a companion paper (Paper II) we make forecasts for measuring the tensor-to-scalar ratio $r$ using future radio arrays.


title: Detecting primordial gravitational waves with circular polarization of the redshifted 21 cm line: II. Forecasts

abstract: In the first paper of this series, we showed that the CMB quadrupole at high redshifts results in a small circular polarization of the emitted 21 cm radiation. In this paper we forecast the sensitivity of future radio experiments to measure the CMB quadrupole during the era of first cosmic light ($z\sim 20$). The tomographic measurement of 21 cm circular polarization allows us to construct a 3D remote quadrupole field. Measuring the $B$-mode component of this remote quadrupole field can be used to put bounds on the tensor-to-scalar ratio $r$. We make Fisher forecasts for a future Fast Fourier Transform Telescope (FFTT), consisting of an array of dipole antennas in a compact grid configuration, as a function of array size and observation time. We find that a FFTT with a side length of 100 km can achieve $\sigma(r)\sim 4\times 10^{-3}$ after ten years of observation and with a sky coverage $f_{\mathrm{sky}}\sim 0.7$. The forecasts are dependent on the evolution of the Lyman-$\alpha$ flux in the pre-reionization era, that remains observationally unconstrained. Finally, we calculate the typical order of magnitudes for circular polarization foregrounds and comment on their mitigation strategies. We conclude that detection of primordial gravitational waves with 21 cm observations is in principle possible, so long as the primordial magnetic field amplitude is small, but would require a very futuristic experiment with corresponding advances in calibration and foreground suppression techniques.


title: Benchmarking Data Analysis and Machine Learning Applications on the Intel KNL Many-Core Processor

abstract: Knights Landing (KNL) is the code name for the second-generation Intel Xeon Phi product family. KNL has generated significant interest in the data analysis and machine learning communities because its new many-core architecture targets both of these workloads. The KNL many-core vector processor design enables it to exploit much higher levels of parallelism. At the Lincoln Laboratory Supercomputing Center (LLSC), the majority of users are running data analysis applications such as MATLAB and Octave. More recently, machine learning applications, such as the UC Berkeley Caffe deep learning framework, have become increasingly important to LLSC users. Thus, the performance of these applications on KNL systems is of high interest to LLSC users and the broader data analysis and machine learning communities. Our data analysis benchmarks of these application on the Intel KNL processor indicate that single-core double-precision generalized matrix multiply (DGEMM) performance on KNL systems has improved by ~3.5x compared to prior Intel Xeon technologies. Our data analysis applications also achieved ~60% of the theoretical peak performance. Also a performance comparison of a machine learning application, Caffe, between the two different Intel CPUs, Xeon E5 v3 and Xeon Phi 7210, demonstrated a 2.7x improvement on a KNL node.


title: Portfolio Risk Assessment using Copula Models

abstract: In the paper, we use and investigate copulas models to represent multivariate dependence in financial time series. We propose the algorithm of risk measure computation using copula models. Using the optimal mean-$CVaR$ portfolio we compute portfolio's Profit and Loss series and corresponded risk measures curves. Value-at-risk and Conditional-Value-at-risk curves were simulated by three copula models: full Gaussian, Student's $t$ and regular vine copula. These risk curves are lower than historical values of the risk measures curve. All three models have superior prediction ability than a usual empirical method. Further directions of research are described.


title: Diffuse galaxy cluster emission at 168 MHz within the Murchison Widefield Array Epoch of Reionization 0-hour field

abstract: We detect and characterise extended, diffuse radio emission from galaxy clusters at 168 MHz within the Epoch of Reionization 0-hour field: a $45^\circ \times 45^\circ$ region of the southern sky centred on R.~A.${}= 0^\circ$, decl.${}=-27^\circ$. We detect 29 sources of interest; a newly detected halo in Abell 0141; a newly detected relic in Abell 2751; 4 new halo candidates and a further 4 new relic candidates; and a new phoenix candidate in Abell 2556. Additionally, we find 9 clusters with unclassifiable, diffuse steep-spectrum emission as well as a candidate double relic system associated with RXC J2351.0-1934. We present measured source properties such as their integrated flux densities, spectral indices ($\alpha$, where $S_\nu \propto \nu^\alpha$), and sizes where possible. We find several of the diffuse sources to have ultra-steep spectra including the halo in Abell 0141, if confirmed, showing $\alpha \leq -2.1 \pm 0.1$ with the present data making it one of the steepest-spectrum haloes known. Finally, we compare our sample of haloes with previously detected haloes and revisit established scaling relations of the radio halo power ($P_{1.4}$) with the cluster X-ray luminosity ($L_{\mathrm{X}}$) and mass ($M_{500}$). We find that the newly detected haloes and candidate haloes are consistent with the $P_{1.4}$-$L_{\mathrm{X}}$ and $P_{1.4}$-$M_{500}$ relations, and see an increase in scatter in the previously found relations with increasing sample size likely caused by inhomogeneous determination of $P_{1.4}$ across the full halo sample. We show that the MWA is capable of detecting haloes and relics within most of the galaxy clusters within the Planck catalogue of Sunyaev-Zel'dovich sources depending on exact halo or relic properties.


title: Two-channel Kondo Effect Emerging from Nd Ions

abstract: We discuss Kondo phenomena in a seven-orbital impurity Anderson model hybridized with $\Gamma_8$ conduction electrons by employing a numerical renormalization group method. In particular, we focus on the case with three local $f$ electrons, corresponding to a Nd$^{3+}$ ion. For realistic values of Coulomb interactions, spin-orbit coupling, cubic crystalline electric field potentials, and hybridization, we find a residual entropy of $0.5\log 2$, a characteristic of two-channel Kondo phenomena, for the wide range of parameters of the local $\Gamma_6$ ground state. This is considered to be the magnetic two-channel Kondo effect, consistent with the result from an extended $s$-$d$ model constructed on the basis of the $j$-$j$ coupling scheme. Finally, we briefly discuss candidates of Nd compounds to observe the two-channel Kondo effect.


title: Affine Hecke Algebras via DAHA

abstract: A method is suggested for obtaining the Plancherel measure for Affine Hecke Algebras as a limit of integral-type formulas for inner products in the polynomial and related modules of Double Affine Hecke Algebras. The analytic continuation necessary here is a generalization of "picking up residues" due to Arthur, Heckman, Opdam and others, which can be traced back to Hermann Weyl. Generally, it is a finite sum of integrals over double affine residual subtori; a complete formula is presented for $A_1$ in the spherical case.


title: Inflation in multi-field random Gaussian landscapes

abstract: We investigate slow-roll inflation in a multi-field random Gaussian landscape. The landscape is assumed to be small-field, with a correlation length much smaller than the Planck scale. Inflation then typically occurs in small patches of the landscape, localized near inflection or saddle points. We find that the inflationary track is typically close to a straight line in the field space, and the statistical properties of inflation are similar to those in a one-dimensional landscape. This picture of multi-field inflation is rather different from that suggested by the Dyson Brownian motion model; we discuss the reasons for this difference. We also discuss tunneling from inflating false vacua to the neighborhood of inflection and saddle points and show that the tunneling endpoints tend to concentrate along the flat direction in the landscape.


title: Crossover in the log-gamma polymer from the replica coordinate Bethe Ansatz

abstract: The coordinate Bethe Ansatz solution of the log-gamma polymer is extended to boundary conditions with one fixed end and the other attached to one half of a one-dimensional lattice. The large-time limit is studied using a saddle-point approximation,and the cumulative distribution function of the rescaled free energy of a long polymer is expressed as a Fredholm determinant. Scaling limits of the kernel are identified, leading to a crossover from the GUE to the GOE Tracy--Widom distributions. The continuum limit reproduces the crossover from droplet to flat initial conditions of the Kardar--Parisi--Zhang equation.


title: An iterative aggregation and disaggregation approach to the calculation of steady-state distributions of continuous processes

abstract: A mapping of the process on a continuous configuration space to the symbolic representation of the motion on a discrete state space will be combined with an iterative aggregation and disaggregation (IAD) procedure to obtain steady state distributions of the process. The IAD speeds up the convergence to the unit eigenvector, which is the steady state distribution, by forming smaller aggregated matrices whose unit eigenvector solutions are used to refine approximations of the steady state vector until convergence is reached. This method works very efficiently and can be used together with distributed or parallel computing methods to obtain high resolution images of the steady state distribution of complex atomistic or energy landscape type problems. The method is illustrated in two numerical examples. In the first example the transition matrix is assumed to be known. The second example represents an overdamped Brownian motion process subject to a dichotomously changing external potential.


title: Broad line-width of antiferromagnetic spinwave due to electrons correlation

abstract: We study magnetic excitations in a bilayer of an antiferromagnetic (AF) insulator and a correlated metal, in which double occupancy is forbidden. The effective action of the AF spin wave in the AF insulator is derived by using the path integral formula within the second order of interplane coupling. The electron correlation in the correlated metal is treated by the Gutzwiller approximation, which renormalizes the hopping integrals by g_t as proportional to the hole density. The linewidth of the AF spin wave excitations originates from particle-hole excitations in the correlated metal. By increasing the correlation effect, i.e., by decreasing g_t, it is found that the linewidth at low energies increases inversely proportional to g_t. The present results will also be useful for bilayers of a metal and ferrimagnet.


title: A mathematical account of the NEGF formalism

abstract: The main goal of this paper is to put on solid mathematical grounds the so-called Non-Equilibrium Green's Function (NEGF) transport formalism for open systems. In particular, we derive the Jauho-Meir-Wingreen formula for the time-dependent current through an interacting sample coupled to non-interacting leads. Our proof is non-perturbative and uses neither complex-time Keldysh contours, nor Langreth rules of 'analytic continuation'. We also discuss other technical identities (Langreth, Keldysh) involving various many body Green's functions. Finally, we study the Dyson equation for the advanced/retarded interacting Green's function and we rigorously construct its (irreducible) self-energy, using the theory of Volterra operators.

categories: math-ph cond-mat.mes-hall cond-mat.stat-mech math.MP quant-ph


title: Entanglement entropy from SU(2) Chern-Simons theory and symmetric webs

abstract: A path integral on a link complement of a three-sphere fixes a vector (the "link state") in Chern-Simons theory. The link state can be written in a certain basis with the colored link invariants as its coefficients. We use symmetric webs to systematically compute the colored link invariants, by which we can write down the multi-partite entangled state of any given link. It is still unknown if a product state necessarily implies that the corresponding components are unlinked, and we leave it as a conjecture.


title: Generalized Ensemble Theory with Non-extensive Statistics

abstract: The non-extensive canonical ensemble theory is reconsidered with the method of Lagrange multipliers by maximizing Tsallis entropy, with the constraint that the normalized term of Tsallis' $q-$average of physical quantities, the sum $\sum p_j^q$, is independent of the probability $p_i$ for Tsallis parameter $q$. The self-referential problem in the deduced probability and thermal quantities in non-extensive statistics is thus avoided, and thermodynamical relationships are obtained in a consistent and natural way. We also extend the study to the non-extensive grand canonical ensemble theory and obtain the $q$-deformed Bose-Einstein distribution as well as the $q$-deformed Fermi-Dirac distribution. The theory is further applied to the generalized Planck law to demonstrate the distinct behaviors of the various generalized $q$-distribution functions discussed in literature.

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


title: Oseba: Optimization for Selective Bulk Analysis in Big Data Processing

abstract: Selective bulk analyses, such as statistical learning on temporal/spatial data, are fundamental to a wide range of contemporary data analysis. However, with the increasingly larger data-sets, such as weather data and marketing transactions, the data organization/access becomes more challenging in selective bulk data processing with the use of current big data processing frameworks such as Spark or keyvalue stores. In this paper, we propose a method to optimize selective bulk analysis in big data processing and referred to as Oseba. Oseba maintains a super index for the data organization in memory to support fast lookup through targeting the data involved with each selective analysis program. Oseba is able to save memory as well as computation in comparison to the default data processing frameworks.


title: Low frame coherence via zero-mean tensor embeddings

abstract: This paper is concerned with achieving optimal coherence for highly redundant real unit-norm frames. As the redundancy grows, the number of vectors in the frame becomes too large to admit equiangular arrangements. In this case, other geometric optimality criteria need to be identified. To this end, we use an iteration of the embedding technique by Conway, Hardin and Sloane. As a consequence of their work, a quadratic mapping embeds equiangular lines into a simplex in a real Euclidean space. Here, higher degree polynomial maps embed highly redundant unit-norm frames to simplices in high-dimensional Euclidean spaces. We focus on the lowest degree case in which the embedding is quartic.


title: Tunable Reactive Synthesis for Lipschitz-Bounded Systems with Temporal Logic Specifications

abstract: We address the problem of synthesizing reactive controllers for cyber-physical systems subject to Signal Temporal Logic (STL) specifications in the presence of adversarial inputs. Given a finite horizon, we define a reactive hierarchy of control problems that differ in the degree of information available to the system about the adversary's actions over the horizon. We show how to construct reactive controllers at various levels of the hierarchy, leveraging the existence of Lipschitz bounds on system dynamics and the quantitative semantics of STL. Our approach, a counterexample-guided inductive synthesis (CEGIS) scheme based on optimization and satisfiability modulo theories (SMT) solving, builds a strategy tree representing the interaction between the system and its environment. In every iteration of the CEGIS loop, we use a mix of optimization and SMT to maximally discard controllers falsified by a given counterexample. Our approach can be applied to any system with local Lipschitz-bounded dynamics, including linear, piecewise-linear and differentially-flat systems. Finally we show an application in the autonomous car domain.


title: A Cluster Elastic Net for Multivariate Regression

abstract: We propose a method for estimating coefficients in multivariate regression when there is a clustering structure to the response variables. The proposed method includes a fusion penalty, to shrink the difference in fitted values from responses in the same cluster, and an L1 penalty for simultaneous variable selection and estimation. The method can be used when the grouping structure of the response variables is known or unknown. When the clustering structure is unknown the method will simultaneously estimate the clusters of the response and the regression coefficients. Theoretical results are presented for the penalized least squares case, including asymptotic results allowing for p >> n. We extend our method to the setting where the responses are binomial variables. We propose a coordinate descent algorithm for both the normal and binomial likelihood, which can easily be extended to other generalized linear model (GLM) settings. Simulations and data examples from business operations and genomics are presented to show the merits of both the least squares and binomial methods.


title: Mutual Interpretability of Robinson Arithmetic and Adjunctive Set Theory with Extensionality

abstract: An elementary rheory of concatenation is introduced and used to establish mutual interpretability of Robinson arithmetic, Minimal Predicative Set Theory, the quantifier-free part of Kirby's finitary set theory, and Adjunctive Set Theory, with or without extensionality.


title: Cell-to-cell variation sets a tissue-rheology-dependent bound on collective gradient sensing

abstract: When a single cell senses a chemical gradient and chemotaxes, stochastic receptor-ligand binding can be a fundamental limit to the cell's accuracy. For clusters of cells responding to gradients, however, there is a critical difference: even genetically identical cells have differing responses to chemical signals. With theory and simulation, we show collective chemotaxis is limited by cell-to-cell variation in signaling. We find that when different cells cooperate the resulting bias can be much larger than the effects of ligand-receptor binding. Specifically, when a strongly-responding cell is at one end of a cell cluster, cluster motion is biased toward that cell. These errors are mitigated if clusters average measurements over times long enough for cells to rearrange. In consequence, fluid clusters are better able to sense gradients: we derive a link between cluster accuracy, cell-to-cell variation, and the cluster rheology. Because of this connection, increasing the noisiness of individual cell motion can actually increase the collective accuracy of a cluster by improving fluidity.

categories: q-bio.CB cond-mat.soft physics.bio-ph q-bio.TO


title: Qualitative types of cosmological evolution in hydrodynamic models with barotropic equation of state

abstract: We study solutions of the Friedmann equations in case of the homogeneous isotropic Universe filled with a perfect fluid. The main points concern the monotony properties of the solutions, the possibility to extend the solutions on all times and occurrence of singularities. We present a qualitative classification of all possible solutions in case of the general smooth barotropic equation of state of the fluid, provided the speed of sound is finite. The list of possible scenarios includes analogs of the "Big Rip" in the future and/or in the past as well as singularity free solutions and oscillating Universes. Extensions of the results to the multicomponent fluids are discussed.


title: $ZZ\gamma$ production in the NLO QCD+EW accuracy at the LHC

abstract: In this paper we present the first study of the impact of the $\mathcal{O}(\alpha)$ EW correction to the $pp \to ZZ \gamma+X$ process at the CERN Large Hadron Collider (LHC). The subsequent $Z$-boson leptonic decays are considered at the leading order using the MadSpin method, which takes into account the spin-correlation and off-shell effects from the $Z$-boson decays. We provide numerical results of the integrated cross section and the kinematic distributions for this process. In coping with final-state photon-jet separation in the QCD real emission and photon-induced processes, we adopt both the Frixione isolated-photon plus jets algorithm and the phenomenological quark-to-photon fragmentation function method for comparison. We find that the next-to-leading order (NLO) EW correction to the $ZZ\gamma$ production can be sizeable and amounts to about $-7\%$ of the integrated cross section, and provides a non$-$negligible contribution to the kinematic distributions, particularly in the high energy region. We conclude that the NLO EW correction should be included in precision theoretical predictions in order to match future experimental accuracy.


title: Influence of growth flux solvent on anneal-tuning of ground states in CaFe2As2

abstract: The effects of anneal-tuning of single-crystalline samples of CaFe$_2$As$_2$ synthesized via a molten Sn flux method are investigated using x-ray diffraction, chemical composition, electrical transport and magnetic susceptibility measurements in order to understand the role of growth conditions on the resultant phase diagram. Previous studies of CaFe$_2$As$_2$ crystals synthesized using a self-flux (FeAs) method revealed an ability to tune the structural and magnetic properties of this system by control of post-synthesis annealing conditions, resulting in an ambient pressure phase diagram that spans from tetragonal/orthorhombic antiferromagnetism to the collapsed tetragonal phase of this system. In this work, we compare previous results to those obtained on crystals synthesized via Sn flux, finding similar tunability in both self- and Sn-flux cases, but less sensitivity to annealing temperatures in the latter case, resulting in a temperature-shifted phase diagram.


title: Application of Dictionary Learning in Alleviating Computational Burden of EEG Source Localization

abstract: Two techniques are proposed to alleviate the computational burden of MUltiple SIgnal Classification (MUSIC) algorithm applied to Electroencephalogram (EEG) source localization. A significant reduction was achieved by parsing the cortex surface into smaller regions and nominating only a few regions for the exhaustive search inherent in the MUSIC algorithm. The nomination procedure involves a dictionary learning phase in which each region is assigned an atom matrix. Moreover, a dimensionality reduction step provided by excluding some of the electrodes is designed such that the Cramer-Rao bound of localization is maintained. It is shown by simulation that computational complexity of the MUSIC-based localization can be reduced by up to $80\%$.


title: Black Hole Growth in Disk Galaxies Mediated by the Secular Evolution of Short Bars

abstract: The growth of black holes (BHs) in disk galaxies lacking classical bulges, which implies an absence of significant mergers, appears to be driven by secular processes. Short bars of sub-kiloparsec radius have been hypothesized to be an important mechanism for driving gas inflows to small scale, feeding central BHs. In order to quantify the maximum BH mass allowed by this mechanism, we examine the robustness of short bars to the dynamical influence of BHs. Large-scale bars are expected to be robust, long-lived structures; extremely massive BHs, which are rare, are needed to completely destroy such bars. However, we find that short bars, which are generally embedded in large-scale outer bars, can be destroyed quickly when BHs of mass $M_{\rm bh}\sim0.05-0.2\%$ of the total stellar mass ($M_\star$) are present. In agreement with this prediction, all galaxies observed to host short bars have BHs with a mass fraction less than $0.2\%M_\star$. Thus the dissolution of short inner bars is possible, perhaps even frequent, in the universe. An important implication of this result is that inner bar-driven gas inflows may be terminated when BHs grow to $\sim0.1\%M_\star$. We predict that $0.2\%M_\star$ is the maximum mass of BHs allowed if they are fed predominately via inner bars. This value matches well the maximum ratio of BH-to-host galaxy stellar mass observed in galaxies with pseudo-bulges and most nearby AGN host galaxies. This hypothesis provides a novel explanation for the lower $M_{\rm bh}/M_\star$ in galaxies which have avoided significant mergers compared with galaxies with classical bulges.


title: An Introduction to the Practical and Theoretical Aspects of Mixture-of-Experts Modeling

abstract: Mixture-of-experts (MoE) models are a powerful paradigm for modeling of data arising from complex data generating processes (DGPs). In this article, we demonstrate how different MoE models can be constructed to approximate the underlying DGPs of arbitrary types of data. Due to the probabilistic nature of MoE models, we propose the maximum quasi-likelihood (MQL) estimator as a method for estimating MoE model parameters from data, and we provide conditions under which MQL estimators are consistent and asymptotically normal. The blockwise minorization-maximizatoin (blockwise-MM) algorithm framework is proposed as an all-purpose method for constructing algorithms for obtaining MQL estimators. An example derivation of a blockwise-MM algorithm is provided. We then present a method for constructing information criteria for estimating the number of components in MoE models and provide justification for the classic Bayesian information criterion (BIC). We explain how MoE models can be used to conduct classification, clustering, and regression and we illustrate these applications via a pair of worked examples.


title: Impact of Angular Spread in Moderately Large MIMO Systems under Pilot Contamination

abstract: Pilot contamination is known to be one of the main bottlenecks for massive multi-input multi-output (MIMO) networks. For moderately large antenna arrays (of importance to recent/emerging deployments) and correlated MIMO, pilot contamination may not be the dominant limiting factor in certain scenarios. To the best of our knowledge, a rigorous characterization of the achievable rates and their explicit dependence on the angular spread (AS) is not available in the existing literature for moderately large antenna array regime. In this paper, considering eigen-beamforming (EBF) precoding, we derive an exact analytical expression for achievable rates in multi-cell MIMO systems under pilot contamination, and characterize the relation between the AS, array size, and respective user rates. Our analytical and simulation results reveal that the achievable rates for both the EBF and the regularized zero-forcing (RZF) precoders follow a non-monotonic behavior for increasing AS when the antenna array size is moderate. We argue that knowledge of this non-monotonic behavior can be exploited to develop effective user-cell pairing techniques.


title: VMEXT: A Visualization Tool for Mathematical Expression Trees

abstract: Mathematical expressions can be represented as a tree consisting of terminal symbols, such as identifiers or numbers (leaf nodes), and functions or operators (non-leaf nodes). Expression trees are an important mechanism for storing and processing mathematical expressions as well as the most frequently used visualization of the structure of mathematical expressions. Typically, researchers and practitioners manually visualize expression trees using general-purpose tools. This approach is laborious, redundant, and error-prone. Manual visualizations represent a user's notion of what the markup of an expression should be, but not necessarily what the actual markup is. This paper presents VMEXT - a free and open source tool to directly visualize expression trees from parallel MathML. VMEXT simultaneously visualizes the presentation elements and the semantic structure of mathematical expressions to enable users to quickly spot deficiencies in the Content MathML markup that does not affect the presentation of the expression. Identifying such discrepancies previously required reading the verbose and complex MathML markup. VMEXT also allows one to visualize similar and identical elements of two expressions. Visualizing expression similarity can support support developers in designing retrieval approaches and enable improved interaction concepts for users of mathematical information retrieval systems. We demonstrate VMEXT's visualizations in two web-based applications. The first application presents the visualizations alone. The second application shows a possible integration of the visualizations in systems for mathematical knowledge management and mathematical information retrieval. The application converts LaTeX input to parallel MathML, computes basic similarity measures for mathematical expressions, and visualizes the results using VMEXT.


title: A Scalable Semidefinite Relaxation Approach to Grid Scheduling

abstract: Determination of the most economic strategies for supply and transmission of electricity is a daunting computational challenge. Due to theoretical barriers, so-called NP-hardness, the amount of effort to optimize the schedule of generating units and route of power, can grow exponentially with the number of decision variables. Practical approaches to this problem involve legacy approximations and ad-hoc heuristics that may undermine the efficiency and reliability of power system operations, that are ever growing in scale and complexity. Therefore, the development of powerful optimization methods for detailed power system scheduling is critical to the realization of smart grids and has received significant attention recently. In this paper, we propose for the first time a computational method, which is capable of solving large-scale power system scheduling problems with thousands of generating units, while accurately incorporating the nonlinear equations that govern the flow of electricity on the grid. The utilization of this accurate nonlinear model, as opposed to its linear approximations, results in a more efficient and transparent market design, as well as improvements in the reliability of power system operations. We design a polynomial-time solvable third-order semidefinite programming (TSDP) relaxation, with the aim of finding a near globally optimal solution for the original NP-hard problem. The proposed method is demonstrated on the largest available benchmark instances from real-world European grid data, for which provably optimal or near-optimal solutions are obtained.


title: Modeling Financial System with Interbank Flows, Borrowing, and Investing

abstract: In our model, private actors with interbank cash flows similar to, but nore general than (Carmona, Fouque, Sun, 2013) borrow from the outside economy at a certain interest rate, controlled by the central bank, and invest in risky assets. Each private actor aims to maximize its expected terminal logarithmic utility. The central bank, in turn, aims to control the overall economy by means of an exponential utility function. We solve all stochastic optimal control problems explicitly. We are able to recreate occasions such as liquidity trap. We study distribution of the number of defaults (net worth of a private actor going below a certain threshold).


title: Computing Entropies With Nested Sampling

abstract: The Shannon entropy, and related quantities such as mutual information, can be used to quantify uncertainty and relevance. However, in practice, it can be difficult to compute these quantities for arbitrary probability distributions, particularly if the probability mass functions or densities cannot be evaluated. This paper introduces a computational approach, based on Nested Sampling, to evaluate entropies of probability distributions that can only be sampled. I demonstrate the method on three examples: a simple gaussian example where the key quantities are available analytically; (ii) an experimental design example about scheduling observations in order to measure the period of an oscillating signal; and (iii) predicting the future from the past in a heavy-tailed scenario.

categories: stat.CO astro-ph.IM cs.IT math.IT physics.data-an


title: Score-informed syllable segmentation for a cappella singing voice with convolutional neural networks

abstract: This paper introduces a new score-informed method for the segmentation of jingju a cappella singing phrase into syllables. The proposed method estimates the most likely sequence of syllable boundaries given the estimated syllable onset detection function (ODF) and its score. Throughout the paper, we first examine the jingju syllables structure and propose a definition of the term "syllable onset". Then, we identify which are the challenges that jingju a cappella singing poses. Further, we investigate how to improve the syllable ODF estimation with convolutional neural networks (CNNs). We propose a novel CNN architecture that allows to efficiently capture different time-frequency scales for estimating syllable onsets. In addition, we propose using a score-informed Viterbi algorithm -instead of thresholding the onset function-, because the available musical knowledge we have (the score) can be used to inform the Viterbi algorithm in order to overcome the identified challenges. The proposed method outperforms the state-of-the-art in syllable segmentation for jingju a cappella singing. We further provide an analysis of the segmentation errors which points possible research directions.


title: Quantum coherence of planar spin models with Dzyaloshinsky-Moriya interaction

abstract: The quantum coherence of one dimensional planar spin models with the Dzyaloshinsky-Moriya interaction is investigated. The anisotropic XY model, the isotropic XX model and the transverse field model are studied in the large N-limit using the two qubit reduced density matrices and the two point correlation functions. From our investigations we find that the coherence as measured using the Jensen-Shannon divergence can be used to detect the quantum phase transitions and the quantum critical points. The derivative of coherence shows non-analytic behavior at the critical points leading to the conclusions that these transitions are of second order. Further we show that the presence of the Dzyaloshinsky-Moriya coupling suppresses the phase transition due to the residual ferromagnetism which is caused by spin canting.


title: Stein's method using approximate zero bias couplings with applications to combinatorial central limit theorems under the Ewens distribution

abstract: We generalize the well-known zero bias distribution and the $\lambda$-Stein pair to an approximate zero bias distribution and an approximate $\lambda,R$-Stein pair, respectively. Berry Esseen type bounds to the normal, based on approximate zero bias couplings and approximate $\lambda,R$-Stein pairs, are obtained using Stein's method. The bounds are then applied to combinatorial central limit theorems where the random permutation has the Ewens $\mathcal{E}_\theta$ distribution with $\theta>0$ which can be specialized to the uniform distribution by letting $\theta=1$. The family of the Ewens distributions appears in the context of population genetics in biology.