title: Physics Program for the STAR/CBM eTOF Upgrade

abstract: The STAR Collaboration and the CBM Collaboration institutions: Heidelberg, Darmstadt, CCNU, Tsinghua, and USTC propose to install an end-cap time-of-flight upgrade (eTOF) to the STAR detector for the RHIC beam energy scan phase II (BES-II) program in 2019 and 2020. BES-II will cover the collision energy range 3.0 to 19.6 GeV. This is the region of interest in the search for a critical point and first-order phase transition, identified by the results from BES-I and by model calculations. For the collider-mode portion of the energy scan, 7.7 to 19.6 GeV, eTOF will extend particle identification (PID) for pions, kaons, and protons to a rapidity of 1.2, complementing the inner Time Projection Chamber (iTPC) upgrade to the forward tracking. The rapidity coverage for PID would extend to only 0.8 without the eTOF upgrade. The eTOF upgrade will enable precision studies of the key bulk property observables, essential to the BES-II search. An internal fixed-target program will allow the energy scan to cover 3.0 to 7.7 GeV. The eTOF upgrade will provide essential mid-rapidity PID for the 4.5 to 7.7 GeV portion of the scan in fixed-target mode. Otherwise there would be a large energy gap in the middle of the BES-II program. A full description of the physics provided by the eTOF upgrade to STAR is presented in this note.


title: Learning Tuple Probabilities

abstract: Learning the parameters of complex probabilistic-relational models from labeled training data is a standard technique in machine learning, which has been intensively studied in the subfield of Statistical Relational Learning (SRL), but---so far---this is still an under-investigated topic in the context of Probabilistic Databases (PDBs). In this paper, we focus on learning the probability values of base tuples in a PDB from labeled lineage formulas. The resulting learning problem can be viewed as the inverse problem to confidence computations in PDBs: given a set of labeled query answers, learn the probability values of the base tuples, such that the marginal probabilities of the query answers again yield in the assigned probability labels. We analyze the learning problem from a theoretical perspective, cast it into an optimization problem, and provide an algorithm based on stochastic gradient descent. Finally, we conclude by an experimental evaluation on three real-world and one synthetic dataset, thus comparing our approach to various techniques from SRL, reasoning in information extraction, and optimization.


title: Intrinsic normalization and extrinsic denormalization of formant data of vowels

abstract: Using a known speaker-intrinsic normalization procedure, formant data are scaled by the reciprocal of the geometric mean of the first three formant frequencies. This reduces the influence of the talker but results in a distorted vowel space. The proposed speaker-extrinsic procedure re-scales the normalized values by the mean formant values of vowels. When tested on the formant data of vowels published by Peterson and Barney, the combined approach leads to well separated clusters by reducing the spread due to talkers. The proposed procedure performs better than two top-ranked normalization procedures based on the accuracy of vowel classification as the objective measure.


title: STAX. An Axion-like Particle Search with Microwave Photons

abstract: We discuss an improved detection scheme for a light-shining-through-wall (LSW) experiment for axion-like particle searches. We propose to use: gyrotrons or klystrons, which can provide extremely intense photon fluxes at frequencies around 30 GHz; transition-edge-sensors (TES) single photon detectors in this frequency domain, with efficiency $\approx1$; high quality factor Fabry-Perot cavities in the microwave domain, both on the photon-axion conversion and photon regeneration sides. We compute that present laboratory exclusion limits on axion-like particles might be improved by at least four orders of magnitude for axion masses $\lesssim 0.02$ meV.


title: Near-Surface Boundary Layer Turbulence Along a Horizontally-Moving, Surface-Piercing Vertical Wall

abstract: The complex interaction between turbulence and the free surface in boundary layer shear flow created by a vertical surface-piercing wall is considered. A laboratory-scale device was built that utilizes a surface-piercing stainless steel belt that travels in a loop around two vertical rollers, with one length of the belt between the rollers acting as a horizontally-moving flat wall. The belt is accelerated suddenly from rest until reaching constant speed in order to create a temporally-evolving boundary layer analogous to the spatially-evolving boundary layer that would exist along a surface-piercing towed flat plate. Surface profiles are measured with a cinematic laser-induced fluorescence system and sub-surface velocity fields are recorded using a high-speed planar particle image velocimetry system. It is found that the belt initially travels through the water without creating any significant waves, before the free surface bursts with activity close to the belt surface. These free surface ripples travel away from the belt before appearing to become freely-propagating waves. From sub-surface velocity measurements, it is found that close to the surface, transition to turbulence happens sooner than far from the surface, leading to an overall thicker boundary layer in the vicinity of the free surface. A secondary peak in streamwise velocity fluctuations accompanies this transition to turbulence and this peak reaches a maximum value a short time later before smoothing outward. Using momentum thickness as a length scale and the streamwise velocity fluctuations at the location of this outer peak as a velocity scale, free surface bursting and air entrainment onset are found to depend in some way on Weber number and agreement is found with scaling arguments for air entrainment presented by Brocchini and Pregrine (J. Fluid Mech., 499, 225-254, 2001).


title: Data assimilation for the heat equation using stabilized finite element methods

abstract: We consider data assimilation for the heat equation using a finite element space semi-discretization. The approach is optimization based, but the design of regularization operators and parameters rely on techniques from the theory of stabilized finite elements. The space semi-discretized system is shown to admit a unique solution. Combining sharp estimates of the numerical stability of the discrete scheme and conditional stability estimates of the ill-posed continuous pde-model we then derive error estimates that reflect the approximation order of the finite element space and the stability of the continuous model. Two different data assimilation situations with different stability properties are considered to illustrate the framework. Full detail on how to adapt known stability estimates for the continuous model to work with the numerical analysis framework is given in appendix.


title: A multi-sensor multi-Bernoulli filter

abstract: In this paper we derive a multi-sensor multi-Bernoulli (MS-MeMBer) filter for multi-target tracking. Measurements from multiple sensors are employed by the proposed filter to update a set of tracks modeled as a multi-Bernoulli random finite set. An exact implementation of the MS-MeMBer update procedure is computationally intractable. We propose an efficient approximate implementation by using a greedy measurement partitioning mechanism. The proposed filter allows for Gaussian mixture or particle filter implementations. Numerical simulations conducted for both linear-Gaussian and non-linear models highlight the improved accuracy of the MS-MeMBer filter and its reduced computational load with respect to the multi-sensor cardinalized probability hypothesis density filter and the iterated-corrector cardinality-balanced multi-Bernoulli filter especially for low probabilities of detection.


title: Preferential sampling of helicity by isotropic helicoids

abstract: We present a theoretical and numerical study on the motion of isotropic helicoids in complex flows. These are particles whose motion is invariant under rotations but not under mirror reflections of the particle. This is the simplest, yet unexplored, extension of the much studied case of small spherical particles. We show that heavy isotropic helicoids, due to the coupling between translational and rotational degrees of freedom, preferentially sample different helical regions in laminar or chaotic advecting flows. This opens the way to control and engineer particles able to track complex flow structures with potential applications to microfluidics and turbulence.


title: Parameterized and Approximation Complexity of Partial VC Dimension

abstract: We introduce the problem Partial VC Dimension that asks, given a hypergraph $H=(X,E)$ and integers $k$ and $\ell$, whether one can select a set $C\subseteq X$ of $k$ vertices of $H$ such that the set $\{e\cap C, e\in E\}$ of distinct hyperedge-intersections with $C$ has size at least $\ell$. The sets $e\cap C$ define equivalence classes over $E$. Partial VC Dimension is a generalization of VC Dimension, which corresponds to the case $\ell=2^k$, and of Distinguishing Transversal, which corresponds to the case $\ell=|E|$ (the latter is also known as Test Cover in the dual hypergraph). We also introduce the associated fixed-cardinality maximization problem Max Partial VC Dimension that aims at maximizing the number of equivalence classes induced by a solution set of $k$ vertices. We study the algorithmic complexity of Partial VC Dimension and Max Partial VC Dimension both on general hypergraphs and on more restricted instances, in particular, neighborhood hypergraphs of graphs.


title: Detection with Multimodal Dependent Data Using Low Dimensional Random Projections

abstract: Performing likelihood ratio based detection with high dimensional multimodal data is a challenging problem since the computation of the joint probability density functions (pdfs) in the presence of inter-modal dependence is difficult. While some computationally expensive approaches have been proposed for dependent multimodal data fusion (e.g., based on copula theory), a commonly used tractable approach is to compute the joint pdf as the product of marginal pdfs ignoring dependence. However, this method leads to poor performance when the data is strongly dependent. In this paper, we consider the problem of detection when dependence among multimodal data is modeled in a compressed domain where compression is obtained using low dimensional random projections. We employ a Gaussian approximation while modeling inter-modal dependence in the compressed domain which is computationally more efficient. We show that, under certain conditions, detection with multimodal dependent data in the compressed domain with a small number of compressed measurements yields enhanced performance compared to detection with high dimensional data via either the product approach or other suboptimal fusion approaches proposed in the literature.


title: Barcodes for Medical Image Retrieval Using Autoencoded Radon Transform

abstract: Using content-based binary codes to tag digital images has emerged as a promising retrieval technology. Recently, Radon barcodes (RBCs) have been introduced as a new binary descriptor for image search. RBCs are generated by binarization of Radon projections and by assembling them into a vector, namely the barcode. A simple local thresholding has been suggested for binarization. In this paper, we put forward the idea of "autoencoded Radon barcodes". Using images in a training dataset, we autoencode Radon projections to perform binarization on outputs of hidden layers. We employed the mini-batch stochastic gradient descent approach for the training. Each hidden layer of the autoencoder can produce a barcode using a threshold determined based on the range of the logistic function used. The compressing capability of autoencoders apparently reduces the redundancies inherent in Radon projections leading to more accurate retrieval results. The IRMA dataset with 14,410 x-ray images is used to validate the performance of the proposed method. The experimental results, containing comparison with RBCs, SURF and BRISK, show that autoencoded Radon barcode (ARBC) has the capacity to capture important information and to learn richer representations resulting in lower retrieval errors for image retrieval measured with the accuracy of the first hit only.


title: Bleach: A Distributed Stream Data Cleaning System

abstract: In this paper we address the problem of rule-based stream data cleaning, which sets stringent requirements on latency, rule dynamics and ability to cope with the unbounded nature of data streams. We design a system, called Bleach, which achieves real-time violation detection and data repair on a dirty data stream. Bleach relies on efficient, compact and distributed data structures to maintain the necessary state to repair data, using an incremental version of the equivalence class algorithm. Additionally, it supports rule dynamics and uses a "cumulative" sliding window operation to improve cleaning accuracy. We evaluate a prototype of Bleach using a TPC-DS derived dirty data stream and observe its high throughput, low latency and high cleaning accuracy, even with rule dynamics. Experimental results indicate superior performance of Bleach compared to a baseline system built on the micro-batch streaming paradigm.


title: Staggered domain wall fermions

abstract: We construct domain wall fermions with a staggered kernel and investigate their spectral and chiral properties numerically in the Schwinger model. In some relevant cases we see an improvement of chirality by more than an order of magnitude as compared to usual domain wall fermions. Moreover, we present first results for four-dimensional quantum chromodynamics, where we also observe significant reductions of chiral symmetry violations for staggered domain wall fermions.


title: Dense Wide-Baseline Scene Flow From Two Handheld Video Cameras

abstract: We propose a new technique for computing dense scene flow from two handheld videos with wide camera baselines and different photometric properties due to different sensors or camera settings like exposure and white balance. Our technique innovates in two ways over existing methods: (1) it supports independently moving cameras, and (2) it computes dense scene flow for wide-baseline scenarios.We achieve this by combining state-of-the-art wide-baseline correspondence finding with a variational scene flow formulation. First, we compute dense, wide-baseline correspondences using DAISY descriptors for matching between cameras and over time. We then detect and replace occluded pixels in the correspondence fields using a novel edge-preserving Laplacian correspondence completion technique. We finally refine the computed correspondence fields in a variational scene flow formulation. We show dense scene flow results computed from challenging datasets with independently moving, handheld cameras of varying camera settings.


title: Adaptive Windowing for ICI Mitigation in Doubly Selective Channels with Unknown Statistics

abstract: In doubly selective channels, receiver windowing constitutes an effective technique for enhancing the banded structure of the frequency-domain channel matrix, and thus improving the effectiveness of a banded equalizer for intercarrier interference (ICI) mitigation. A common window design technique, which performs close to optimal, is based on the criterion of maximum average signal-to-interference-plus-noise ratio (SINR). The optimality of this technique has been verified for stationary channels with perfectly known statistics. However, in cases where this assumption does not hold, a near optimal performance can be achieved at the expense of high complexity cost. To overcome these limitations, an adaptive windowing technique is proposed that is able to track the optimal receiver window offering low-complexity requirements. Through simulation experiments it has been verified that the proposed technique is able to adapt to the varying channel statistics with increased robustness to channel modeling errors.


title: Troisi\`eme groupe de cohomologie non ramifi\'ee des torseurs universels sur les surfaces rationnelles

abstract: Let $k$ a field of characteristic zero. Let $X$ be a smooth, projective, geometrically rational $k$-surface. Let $\mathcal{T}$ be a universal torsor over $X$ with a $k$-point et $\mathcal{T}^c$ a smooth compactification of $\mathcal{T}$. There is an open question: is $\mathcal{T}^c$ $k$-birationally equivalent to a projective space? We know that the unramified cohomology groups of degree 1 and 2 of $\mathcal{T}$ and $\mathcal{T}^c$ are reduced to their constant part. For the analogue of the third cohomology groups, we give a sufficient condition using the Galois structure of the geometrical Picard group of $X$. This enables us to show that $H^{3}_{nr}(\mathcal{T}^{c},\mathbb{Q}/\mathbb{Z}(2))/H^3(k,\mathbb{Q}/\mathbb{Z}(2))$ vanishes if $X$ is a generalised Ch\^atelet surface and that this group is reduced to its $2$-primary part if $X$ is a del Pezzo surface of degree at least 2.


title: Radon-Gabor Barcodes for Medical Image Retrieval

abstract: In recent years, with the explosion of digital images on the Web, content-based retrieval has emerged as a significant research area. Shapes, textures, edges and segments may play a key role in describing the content of an image. Radon and Gabor transforms are both powerful techniques that have been widely studied to extract shape-texture-based information. The combined Radon-Gabor features may be more robust against scale/rotation variations, presence of noise, and illumination changes. The objective of this paper is to harness the potentials of both Gabor and Radon transforms in order to introduce expressive binary features, called barcodes, for image annotation/tagging tasks. We propose two different techniques: Gabor-of-Radon-Image Barcodes (GRIBCs), and Guided-Radon-of-Gabor Barcodes (GRGBCs). For validation, we employ the IRMA x-ray dataset with 193 classes, containing 12,677 training images and 1,733 test images. A total error score as low as 322 and 330 were achieved for GRGBCs and GRIBCs, respectively. This corresponds to $\approx 81\%$ retrieval accuracy for the first hit.


title: Deep Impression: Audiovisual Deep Residual Networks for Multimodal Apparent Personality Trait Recognition

abstract: Here, we develop an audiovisual deep residual network for multimodal apparent personality trait recognition. The network is trained end-to-end for predicting the Big Five personality traits of people from their videos. That is, the network does not require any feature engineering or visual analysis such as face detection, face landmark alignment or facial expression recognition. Recently, the network won the third place in the ChaLearn First Impressions Challenge with a test accuracy of 0.9109.


title: Triangular reductions of 2D Toda hierarchy

abstract: New reductions of the 2D Toda equations associated with low-triangular difference operators are proposed. Their explicit Hamiltonian description is obtained.


title: Cooper pairs and bipolarons

abstract: It is shown that Cooper pairs are a solution of the bipolaron problem for model Fr\"{o}hlich Hamiltonian. The total energy of a pair for the initial Fr\"{o}hlich Hamiltonian is found. Differences between the solutions for the model and initial two-particle problems are discussed.


title: Search for anomalous electroweak production of $WW/WZ$ in association with a high-mass dijet system in $pp$ collisions at $\sqrt{s}=8$ TeV with the ATLAS detector

abstract: A search is presented for anomalous quartic gauge boson couplings in vector-boson scattering. The data for the analysis correspond to $20.2$ fb$^{-1}$ of $\sqrt{s}=8$ TeV $pp$ collisions, and were collected in 2012 by the ATLAS experiment at the Large Hadron Collider. The search looks for the production of $WW$ or $WZ$ boson pairs accompanied by a high-mass dijet system, with one $W$ decaying leptonically, and a $W$ or $Z$ decaying hadronically. The hadronically decaying $W/Z$ is reconstructed as either two small-radius jets or one large-radius jet using jet substructure techniques. Constraints on the anomalous quartic gauge boson coupling parameters $\alpha_4$ and $\alpha_5$ are set by fitting the transverse mass of the diboson system, and the resulting 95% confidence intervals are $-0.024<\alpha_4<0.030$ and $-0.028<\alpha_5<0.033$.


title: Learning Opposites Using Neural Networks

abstract: Many research works have successfully extended algorithms such as evolutionary algorithms, reinforcement agents and neural networks using "opposition-based learning" (OBL). Two types of the "opposites" have been defined in the literature, namely \textit{type-I} and \textit{type-II}. The former are linear in nature and applicable to the variable space, hence easy to calculate. On the other hand, type-II opposites capture the "oppositeness" in the output space. In fact, type-I opposites are considered a special case of type-II opposites where inputs and outputs have a linear relationship. However, in many real-world problems, inputs and outputs do in fact exhibit a nonlinear relationship. Therefore, type-II opposites are expected to be better in capturing the sense of "opposition" in terms of the input-output relation. In the absence of any knowledge about the problem at hand, there seems to be no intuitive way to calculate the type-II opposites. In this paper, we introduce an approach to learn type-II opposites from the given inputs and their outputs using the artificial neural networks (ANNs). We first perform \emph{opposition mining} on the sample data, and then use the mined data to learn the relationship between input $x$ and its opposite $\breve{x}$. We have validated our algorithm using various benchmark functions to compare it against an evolving fuzzy inference approach that has been recently introduced. The results show the better performance of a neural approach to learn the opposites. This will create new possibilities for integrating oppositional schemes within existing algorithms promising a potential increase in convergence speed and/or accuracy.


title: Critical Rayleigh number of for error function temperature profile with a quasi-static assumption

abstract: When a semi-infinite body is heated from below by a sudden increase in temperature (or cooled from above) an error function temperature profile grows as the heat diffuses into the fluid. The stability of such a profile is investigated using a large-wavelength asymptotic expansion under the quasi-static, or frozen-time, approximation. The critical Rayleigh number for this layer is found to be $Ra=\pi^{1/2}$ based on the length-scale $(\kappa t)^{1/2}$ where $\kappa$ is the thermal diffusivity and $t$ the time since the onset of heating.


title: Counting unstable eigenvalues in Hamiltonian spectral problems via commuting operators

abstract: We present a general counting result for the unstable eigenvalues of linear operators of the form $JL$ in which $J$ and $L$ are skew- and self-adjoint operators, respectively. Assuming that there exists a self-adjoint operator $K$ such that the operators $JL$ and $JK$ commute, we prove that the number of unstable eigenvalues of $JL$ is bounded by the number of nonpositive eigenvalues of~$K$. As an application, we discuss the transverse stability of one-dimensional periodic traveling waves in the classical KP-II (Kadomtsev--Petviashvili) equation. We show that these one-dimensional periodic waves are transversely spectrally stable with respect to general two-dimensional bounded perturbations, including periodic and localized perturbations in either the longitudinal or the transverse direction, and that they are transversely linearly stable with respect to doubly periodic perturbations.

categories: nlin.SI math-ph math.AP math.DS math.MP nlin.PS


title: On regular Stein neighborhoods of a union of two maximal totally real subspaces in $\mathbb{C}^n$

abstract: We present a construction of regular Stein neighborhoods of a union of maximally totally real subspaces $M=(A+iI)\mathbb{R}^n$ and $N=\mathbb{R}^n$ in $\mathbb{C}^n$, provided that the entries of a real $n \times n$ matrix $A$ are sufficiently small. Our proof is based on a local construction of a suitable plurisubharmonic function $\rho$ near the origin, such that the sublevel sets of $\rho$ are strongly pseudoconvex and admit strong deformation retraction to $M\cup N$. We also give the application of this result to totally real immersions of real $n$-manifolds in $\mathbb{C}^n$ with only finitely many double points, and such that the union of the tangent spaces at each intersection in some local coordinates coincides with $M\cup N$, described above.


title: Skew plane partitions according to the $m$th largest and $m$th smallest parts

abstract: We extend recent results by G. E. Andrews and G. Simay on the $m$th largest and $m$th smallest parts of a partition to the more general context of skew plane partitions. In order to do this, we introduce new objects called skew plane overpartitions.


title: Epidemic Processes over Time-Varying Networks

abstract: The spread of viruses in biological networks, computer networks, and human contact networks can have devastating effects; developing and analyzing mathematical models of these systems can be insightful and lead to societal benefits. Prior research has focused mainly on network models with static graph structures, however the systems being modeled typically have dynamic graph structures. Therefore to better understand and analyze virus spread, further study is required. In this paper, we consider virus spread models over networks with dynamic graph structures, and investigate the behavior of diseases in these systems. A stability analysis of epidemic processes over time-varying networks is performed, examining conditions for the disease free equilibrium, in both the deterministic and stochastic cases. We present simulation results, propose a number of corollaries based on these simulations, and discuss quarantine control via simulation.


title: A geometric invariant characterising initial data for the Kerr-Newman spacetime

abstract: We describe the construction of a geometric invariant characterising initial data for the Kerr-Newman spacetime. This geometric invariant vanishes if and only if the initial data set corresponds to exact Kerr-Newman initial data, and so characterises this type of data. We first illustrate the characterisation of the Kerr-Newman spacetime in terms of Killing spinors. The space spinor formalism is then used to obtain a set of four independent conditions on an initial Cauchy hypersurface that guarantee the existence of a Killing spinor on the development of the initial data. Following a similar analysis in the vacuum case, we study the properties of solutions to the approximate Killing spinor equation and use them to construct the geometric invariant.


title: SemanticFusion: Dense 3D Semantic Mapping with Convolutional Neural Networks

abstract: Ever more robust, accurate and detailed mapping using visual sensing has proven to be an enabling factor for mobile robots across a wide variety of applications. For the next level of robot intelligence and intuitive user interaction, maps need extend beyond geometry and appearence - they need to contain semantics. We address this challenge by combining Convolutional Neural Networks (CNNs) and a state of the art dense Simultaneous Localisation and Mapping (SLAM) system, ElasticFusion, which provides long-term dense correspondence between frames of indoor RGB-D video even during loopy scanning trajectories. These correspondences allow the CNN's semantic predictions from multiple view points to be probabilistically fused into a map. This not only produces a useful semantic 3D map, but we also show on the NYUv2 dataset that fusing multiple predictions leads to an improvement even in the 2D semantic labelling over baseline single frame predictions. We also show that for a smaller reconstruction dataset with larger variation in prediction viewpoint, the improvement over single frame segmentation increases. Our system is efficient enough to allow real-time interactive use at frame-rates of approximately 25Hz.


title: Global Sensitivity Analysis in a Mathematical Model of the Renal Interstitium

abstract: The pressure in the renal interstitium is an important factor for normal kidney function. Here we develop a computational model of the rat kidney and use it to investigate the relationship between arterial blood pressure and interstitial fluid pressure. In addition, we investigate how tissue flexibility influences this relationship. Due to the complexity of the model, the large number of parameters, and the inherent uncertainty of the experimental data, we utilize Monte Carlo sampling to study the model's behavior under a wide range of parameter values and to compute first- and total-order sensitivity indices. Characteristically, at elevated arterial blood pressure, the model predicts cases with increased or reduced interstitial pressure. The transition between the two cases is controlled mostly by the compliance of the blood vessels located before the afferent arterioles.


title: Low Complexity Multiply Accumulate Unit for Weight-Sharing Convolutional Neural Networks

abstract: Convolutional Neural Networks (CNNs) are one of the most successful deep machine learning technologies for processing image, voice and video data. CNNs require large amounts of processing capacity and memory, which can exceed the resources of low power mobile and embedded systems. Several designs for hardware accelerators have been proposed for CNNs which typically contain large numbers of Multiply Accumulate (MAC) units. One approach to reducing data sizes and memory traffic in CNN accelerators is "weight sharing", where the full range of values in a trained CNN are put in bins and the bin index is stored instead of the original weight value. In this paper we propose a novel MAC circuit that exploits binning in weight-sharing CNNs. Rather than computing the MAC directly we instead count the frequency of each weight and place it in a bin. We then compute the accumulated value in a subsequent multiply phase. This allows hardware multipliers in the MAC circuit to be replaced with adders and selection logic. Experiments show that for the same clock speed our approach results in fewer gates, smaller logic, and reduced power.


title: On finite determinacy for matrices of power series

abstract: Let $R=K[[x_1,...,x_s]]$ be the ring of formal power series with maximal ideal $\mathfrak{m}$ over a field $K$ of arbitrary characteristic. On the ring $M_{m,n}$ of $m\times n$ matrices $A$ with entries in $R$ we consider several equivalence relations given by the action on $M_{m,n}$ of a group $G$. $G$ can be the group of automorphisms of $R$, combined with the multiplication of invertible matrices from the left, from the right, or from both sides, respectively. We call $A$ finitely $G$-determined if $A$ is $G$-equivalent to any matrix $B$ with ${A-B} \in {\mathfrak{m}^k M_{m,n}}$ for some finite integer $k$, which implies in particular that $A$ is $G$--equivalent to a matrix with polynomial entries. The classical criterion for analytic or differential map germs $f:(K^s,0) \to (K^m,0)$, $K = \mathbb{R}, \mathbb{C}$, says that $f \in M_{m,1}$ is finitely determined (with respect to various group actions) iff the tangent space to the orbit of $f$ has finite codimension in $M_{m,1}$. We extend this criterion to arbitrary matrices in $M_{m,n}$ if the characteristic of K is 0 or, more general, if the orbit map is separable. In positive characteristic however, the problem is more subtle since the orbit map is in general not separable, as we show by an example. This fact had been overlooked in previous papers. Our main result is a general sufficient criterion for finite $G$-determinacy in $M_{m,n}$ in arbitrary characteristic in terms of the tangent image of the orbit map, which we introduce in this paper. This criterion provides a computable bound for the $G$-determinacy of a matrix $A$ in $M_{m,n}$, which is new even in characteristic 0.


title: Intrinsic coherence in assisted sub-state discrimination

abstract: We study intrinsic coherence in the tripartite process to unambiguously discriminate two nonorthogonal states of a qubit, entangled with another one, and assisted by an auxiliary system. The optimal success probability is found to be benefited by initial intrinsic coherence, but no extra one is required. The transformations among different contributions of intrinsic coherence are necessary in this procedure, which increase with the overlap between the states to recognize. Such state discrimination is a key step of the probabilistic teleportation protocol. Entanglement of the quantum channel decreases the coherence characterizing the reliance on an ancilla.


title: Virtual Micromagnetics: A Framework for Accessible and Reproducible Micromagnetic Simulation

abstract: Computational micromagnetics requires numerical solution of partial differential equations to resolve complex interactions in magnetic nanomaterials. The Virtual Micromagnetics project described here provides virtual machine simulation environments to run open-source micromagnetic simulation packages. These environments allow easy access to simulation packages that are often difficult to compile and install, and enable simulations and their data to be shared and stored in a single virtual hard disk file, which encourages reproducible research. Virtual Micromagnetics can be extended to automate the installation of micromagnetic simulation packages on non-virtual machines, and to support closed-source and new open-source simulation packages, including packages from disciplines other than micromagnetics, encouraging reuse. Virtual Micromagnetics is stored in a public GitHub repository under a three-clause Berkeley Software Distribution (BSD) license.


title: Hardness Results for Consensus-Halving

abstract: We study the consensus-halving problem of dividing an object into two portions, such that each of $n$ agents has equal valuation for the two portions. The $\epsilon$-approximate consensus-halving problem allows each agent to have an $\epsilon$ discrepancy on the values of the portions. We prove that computing $\epsilon$-approximate consensus-halving solution using $n$ cuts is in PPA, and is PPAD-hard, where $\epsilon$ is some positive constant; the problem remains PPAD-hard when we allow a constant number of additional cuts. It is NP-hard to decide whether a solution with $n-1$ cuts exists for the problem. As a corollary of our results, we obtain that the approximate computational version of the Continuous Necklace Splitting Problem is PPAD-hard when the number of portions $t$ is two.


title: A unifying framework for fast randomization of ecological networks with fixed (node) degrees

abstract: The switching model is a Markov chain approach to sample graphs with fixed degree sequence uniformly at random. The recently invented Curveball algorithm for bipartite graphs applies several switches simultaneously (`trades'). Here, we introduce Curveball algorithms for simple (un)directed graphs which use single or simultaneous trades. We show experimentally that these algorithms converge magnitudes faster than the corresponding switching models.


title: Heavy flavor Production and Interactions in Relativistic Heavy-Ion Collisions in CMS Experiment

abstract: This paper presents the CMS measurements of quarkonia and open heavy flavor production in \pp, \pPb, and \PbPb collisions at \sqrtsnn = 2.76 and 5.02 TeV. A brief outlook of the near-future CMS heavy flavor physics analyses is provided at the end.


title: Existence of weak solutions for a general porous medium equation with nonlocal pressure

abstract: We study the general nonlinear diffusion equation $u_t=\nabla\cdot (u^{m-1}\nabla (-\Delta)^{-s}u)$ that describes a flow through a porous medium which is driven by a nonlocal pressure. We consider constant parameters $m>1$ and $0<s<1$, we assume that the solutions are non-negative and the problem is posed in the whole space. In this paper we prove existence of weak solutions for all integrable initial data $u_0 \ge 0$ and for all exponents $m>1$ by developing a new approximation method that allows to treat the range $m\ge 3$ that could not be covered by previous works. We also extend the class of initial data to include any non-negative measure $\mu$ with finite mass. In passing from bounded initial data to measure data we make strong use of an $L^1$-$L^\infty$ smoothing effect and other functional estimates. Finite speed of propagation is established for all $m\ge 2$, and this property implies the existence of free boundaries. The authors had already proved that finite propagation does not hold for $m<2$.


title: The Option-Critic Architecture

abstract: Temporal abstraction is key to scaling up learning and planning in reinforcement learning. While planning with temporally extended actions is well understood, creating such abstractions autonomously from data has remained challenging. We tackle this problem in the framework of options [Sutton, Precup & Singh, 1999; Precup, 2000]. We derive policy gradient theorems for options and propose a new option-critic architecture capable of learning both the internal policies and the termination conditions of options, in tandem with the policy over options, and without the need to provide any additional rewards or subgoals. Experimental results in both discrete and continuous environments showcase the flexibility and efficiency of the framework.


title: STEM: A Scheme for Two-phase Evaluation of Majority Logic

abstract: The switching time of a magnet in a spin current based majority gate depends on the input vector combination, and this often restricts the speed of majority-based circuits. To address this issue, this work proposes a novel two-phase scheme to implement majority logic and evaluates it on an all-spin logic (ASL) majority-based logic structures. In Phase 1, the output is initialized to a preset value. Next in Phase 2, the inputs are evaluated to switch the output magnet to its correct value. The time window for the output to switch in Phase 2 is fixed. Using such a scheme, an n-input AND gate which requires a total of (2n-1) inputs in the conventional implementation can now be implemented with only (n+1) inputs. When applied to standard logic functions, it is demonstrated that the proposed method of designing ASL gates are 1.6-3.4X faster and 1.9-6.9X more energy-efficient than the conventional method, and for a five-magnet full adder, it is shown that the proposed ASL implementation is 1.5X faster, 2.2X more energy-efficient, and provides a 16% improvement in area.


title: Spectral geometry of the Steklov problem on orbifolds

abstract: We consider how the geometry and topology of a compact $n$-dimensional Riemannian orbifold with boundary relates to its Steklov spectrum. In two dimensions, motivated by work of A. Girouard, L. Parnovski, I. Polterovich and D. Sher in the manifold setting, we compute the precise asymptotics of the Steklov spectrum in terms of only boundary data. As a consequence, we prove that the Steklov spectrum detects the presence and number of orbifold singularities on the boundary of an orbisurface and it detects the number each of smooth and singular boundary components. Moreover, we find that the Steklov spectrum also determines the lengths of the boundary components modulo an equivalence relation, and we show by examples that this result is the best possible. We construct various examples of Steklov isospectral Riemannian orbifolds which demonstrate that these two-dimensional results do not extend to higher dimensions. In addition, we give two-dimensional examples which show that the Steklov spectrum does \emph{not} detect the presence of interior singularities nor does it determine the orbifold Euler characteristic. In fact, a flat disk is Steklov isospectral to a cone. In another direction, we obtain upper bounds on the Steklov eigenvalues of a Riemannian orbifold in terms of the isoperimetric ratio and a conformal invariant. We generalize results of B. Colbois, A. El Soufi and A. Girouard, and the fourth author to the orbifold setting; in the process, we gain a sharpness result on these bounds that was not evident in the manifold setting. In dimension two, our eigenvalue bounds are solely in terms of the orbifold Euler characteristic and the number each of smooth and singular boundary components.


title: Target-driven Visual Navigation in Indoor Scenes using Deep Reinforcement Learning

abstract: Two less addressed issues of deep reinforcement learning are (1) lack of generalization capability to new target goals, and (2) data inefficiency i.e., the model requires several (and often costly) episodes of trial and error to converge, which makes it impractical to be applied to real-world scenarios. In this paper, we address these two issues and apply our model to the task of target-driven visual navigation. To address the first issue, we propose an actor-critic model whose policy is a function of the goal as well as the current state, which allows to better generalize. To address the second issue, we propose AI2-THOR framework, which provides an environment with high-quality 3D scenes and physics engine. Our framework enables agents to take actions and interact with objects. Hence, we can collect a huge number of training samples efficiently. We show that our proposed method (1) converges faster than the state-of-the-art deep reinforcement learning methods, (2) generalizes across targets and across scenes, (3) generalizes to a real robot scenario with a small amount of fine-tuning (although the model is trained in simulation), (4) is end-to-end trainable and does not need feature engineering, feature matching between frames or 3D reconstruction of the environment. The supplementary video can be accessed at the following link: https://youtu.be/SmBxMDiOrvs.


title: Photoproduction of a large invariant mass $\gamma \rho$ pair at small momentum transfer

abstract: The collinear factorization framework allows to describe the exclusive photoproduction of a $\gamma\,\rho$ pair in the generalized Bjorken regime in terms of a perturbatively calculable coefficient function and universal generalized parton distributions. The kinematics are defined by a large invariant mass of the $\gamma \rho$ pair and a small transverse momentum of the final nucleon. We calculate the scattering amplitude at leading order in $\alpha_s$ and the differential cross sections for the process where the $\rho-$meson is either longitudinally or transversely polarized, in the kinematics of the near future Jlab experiments. Our estimate of the cross section demonstrates that this process is measurable at JLab 12-GeV.


title: The $A^2$ asymmetry and propagators in lattice $SU(2)$ gluodynamics at $T>T_c$

abstract: We study numerically the chromoelectric-chromomagnetic asymmetry of the dimension two $A^2$ gluon condensate as well as the transverse and longitudinal gluon propagators at $T>T_c$ in the Landau-gauge $SU(2)$ lattice gauge theory with a particular emphasis on finite-volume effects. We show that previously found so called symmetric point at which asymmetry changes sign is an artifact of the finite volume effects. We find that with increasing temperature the asymmetry decreases approaching zero value from above in agreement with perturbative result. Instead of the asymmetry we suggest the ratio of the transverse to longitudinal propagator taken at zero momentum as an indicator of the boundary of the postconfinement domain and find it at $T \simeq 1.7 T_c$.


title: On asymptotic dynamics for $L^2$ critical generalized KdV equations with a saturated perturbation

abstract: In this paper, we consider the $L^2$ critical gKdV equation with a saturated perturbation: $\partial_t u+(u_{xx}+u^5-\gamma u|u|^{q-1})_x=0$, where $q>5$ and $0<\gamma\ll1$. For any initial data $u_0\in H^1$, the corresponding solution is always global and bounded in $H^1$. This equation has a family of solitons, and our goal is to classify the dynamics near soliton. Together with a suitable decay assumption, there are only three possibilities: (i) the solution converges asymptotically to a solitary wave, whose $H^1$ norm is of size $\gamma^{-2/(q-1)}$, as $\gamma\rightarrow0$; (ii) the solution is always in a small neighborhood of the modulated family of solitary waves, but blows down at $+\infty$; (iii) the solution leaves any small neighborhood of the modulated family of the solitary waves. This extends the classification of the rigidity dynamics near the ground state for the unperturbed $L^2$ critical gKdV (corresponding to $\gamma=0$) by Martel, Merle and Rapha\"el. However, the blow-down behavior (ii) is completely new, and the dynamics of the saturated equation cannot be viewed as a perturbation of the $L^2$ critical dynamics of the unperturbed equation. This is the first example of classification of the dynamics near ground state for a saturated equation in this context. The cases of $L^2$ critical NLS and $L^2$ supercritical gKdV, where similar classification results are expected, are completely open.


title: High Contrast Imaging in the Visible: First Experimental Results at the Large Binocular Telescope

abstract: In February 2014, the SHARK-VIS (System for High contrast And coronography from R to K at VISual bands) Forerunner, a high contrast experimental imager operating at visible wavelengths, was installed at LBT (Large Binocular Telescope). Here we report on the first results obtained by recent on-sky tests. These results show the extremely good performance of the LBT ExAO (Extreme Adaptive Optics) system at visible wavelengths, both in terms of spatial resolution and contrast achieved. Similarly to what was done by (Amara et al. 2012), we used the SHARK-VIS Forerunner data to quantitatively assess the contrast enhancement. This is done by injecting several different synthetic faint objects in the acquired data and applying the ADI (angular differential imaging) technique. A contrast of the order of $5 \times 10^{-5}$ is obtained at 630 nm for angular separations from the star larger than 100 mas. These results are discussed in light of the future development of SHARK-VIS and compared to those obtained by other high contrast imagers operating at similar wavelengths.


title: Discovering and Deciphering Relationships Across Disparate Data Modalities

abstract: Understanding the relationships between different properties of data, such as whether a connectome or genome has information about disease status, is becoming increasingly important in modern biological datasets. While existing approaches can test whether two properties are related, they often require unfeasibly large sample sizes in real data scenarios, and do not provide any insight into how or why the procedure reached its decision. Our approach, "Multiscale Graph Correlation" (MGC), is a dependence test that juxtaposes previously disparate data science techniques, including k-nearest neighbors, kernel methods (such as support vector machines), and multiscale analysis (such as wavelets). Other methods typically require double or triple the number samples to achieve the same statistical power as MGC in a benchmark suite including high-dimensional and nonlinear relationships - spanning polynomial (linear, quadratic, cubic), trigonometric (sinusoidal, circular, ellipsoidal, spiral), geometric (square, diamond, W-shape), and other functions, with dimensionality ranging from 1 to 1000. Moreover, MGC uniquely provides a simple and elegant characterization of the potentially complex latent geometry underlying the relationship, providing insight while maintaining computational efficiency. In several real data applications, including brain imaging and cancer genetics, MGC is the only method that can both detect the presence of a dependency and provide specific guidance for the next experiment and/or analysis to conduct.


title: Hierarchical Quantum Master Equation Approach to Electronic-Vibrational Coupling in Nonequilibrium Transport through Nanosystems

abstract: Within the hierarchical quantum master equation (HQME) framework, an approach is presented, which allows a numerically exact description of nonequilibrium charge transport in nanosystems with strong electronic-vibrational coupling. The method is applied to a generic model of vibrationally coupled transport considering a broad spectrum of parameters ranging from the nonadiabatic to the adiabatic regime and including both resonant and off-resonant transport. We show that nonequilibrium effects are important in all these regimes. In particular in the off-resonant transport regime, the inelastic co-tunneling signal is analyzed for a vibrational mode in full nonequilibrium, revealing a complex interplay of different transport processes and deviations from the commonly used $G_0/2$-thumb-rule. In addition, the HQME-approach is used to benchmark approximate master equation and nonequilibrium Green's function methods.


title: Lamda*-Closed sets and new separation axioms in Alexandroff spaces

abstract: Here we have studied the ideas of g*-closed sets, g^tou -sets and Lamda*-closed sets and investigate some of their properties in the spaces of A. D. Alexandroff [1]. We have also studied few separation axioms like T-omega/4,T-3omega/8, T-omega in Alexandroff spaces and have introduced a new separation axiom namely T-5omega/8 axiom in this space.


title: Random Search with Memory in Patchy Media: Exploration-Exploitation Tradeoff

abstract: How to best exploit patchy resources? This long-standing question belongs to the extensively studied class of explore/exploit problems that arise in a wide range of situations, from animal foraging, to robotic exploration, and to human decision processes. Despite its broad relevance, the issue of optimal exploitation has previously only been tackled through two paradigmatic limiting models---patch-use and random search---that do not account for the interplay between searcher motion within patches and resource depletion. Here, we bridge this gap by introducing a minimal patch exploitation model that incorporates this coupling: the searcher depletes the resources along its random-walk trajectory within a patch and travels to a new patch after it takes $\mathcal{S}$ consecutive steps without finding resources. We compute the distribution of the amount of resources $F_t$ consumed by time $t$ for this non-Markovian random walker and show that exploring multiple patches is beneficial. In one dimension, we analytically derive the optimal strategy to maximize $F_t$. We show that this strategy is robust with respect to the distribution of resources within patches and the criterion for leaving a given patch. We also show that $F_t$ can be optimized in the ecologically-relevant case of two-dimensional patchy environments.


title: Style Imitation and Chord Invention in Polyphonic Music with Exponential Families

abstract: Modeling polyphonic music is a particularly challenging task because of the intricate interplay between melody and harmony. A good model should satisfy three requirements: statistical accuracy (capturing faithfully the statistics of correlations at various ranges, horizontally and vertically), flexibility (coping with arbitrary user constraints), and generalization capacity (inventing new material, while staying in the style of the training corpus). Models proposed so far fail on at least one of these requirements. We propose a statistical model of polyphonic music, based on the maximum entropy principle. This model is able to learn and reproduce pairwise statistics between neighboring note events in a given corpus. The model is also able to invent new chords and to harmonize unknown melodies. We evaluate the invention capacity of the model by assessing the amount of cited, re-discovered, and invented chords on a corpus of Bach chorales. We discuss how the model enables the user to specify and enforce user-defined constraints, which makes it useful for style-based, interactive music generation.


title: The Emergence of Solar Supergranulation as a Natural Consequence of Rotationally-Constrained Interior Convection

abstract: We investigate how rotationally-constrained, deep convection might give rise to supergranulation, the largest distinct spatial scale of convection observed in the solar photosphere. While supergranulation is only weakly influenced by rotation, larger spatial scales of convection sample the deep convection zone and are presumably rotationally influenced. We present numerical results from a series of nonlinear, 3-D simulations of rotating convection and examine the velocity power distribution realized under a range of Rossby numbers. When rotation is present, the convective power distribution possesses a pronounced peak, at characteristic wavenumber $\ell_\mathrm{peak}$, whose value increases as the Rossby number is decreased. This distribution of power contrasts with that realized in non-rotating convection, where power increases monotonically from high to low wavenumbers. We find that spatial scales smaller than $\ell_\mathrm{peak}$ behave in analogy to non-rotating convection. Spatial scales larger than $\ell_\mathrm{peak}$ are rotationally-constrained and possess substantially reduced power relative to the non-rotating system. We argue that the supergranular scale emerges due to a suppression of power on spatial scales larger than $\ell\approx100$ owing to the presence of deep, rotationally-constrained convection. Supergranulation thus represents the largest non-rotationally-constrained mode of solar convection. We conclude that the characteristic spatial scale of supergranulation bounds that of the deep convective motions from above, making supergranulation an indirect measure of the deep-seated dynamics at work in the solar dynamo. Using the spatial scale of supergranulation in conjunction with our numerical results, we estimate an upper bound of 10 m s$^{-1}$ for the Sun's bulk rms convective velocity.


title: Relatively hyperbolic groups with fixed peripherals

abstract: We build quasi--isometry invariants of relatively hyperbolic groups which detect the hyperbolic parts of the group; these are variations of the stable dimension constructions previously introduced by the authors. We prove that, given any finite collection of finitely generated groups $\mathcal{H}$ each of which either has finite stable dimension or is non-relatively hyperbolic, there exist infinitely many quasi--isometry types of one--ended groups which are hyperbolic relative to $\mathcal{H}$. The groups are constructed using small cancellation theory over free products.


title: Ultrafast Magnetic Switching of GdFeCo with Electronic Heat Currents

abstract: We report the magnetic response of Au/GdFeCo bilayers to optical irradiation of the Au surface. For bilayers with Au thickness greater than 50 nm, the great majority of energy is absorbed by the Au electrons, creating an initial temperature differential of thousands of Kelvin between the Au and GdFeCo layers. The resulting electronic heat currents between the Au and GdFeCo layers last for several picoseconds with energy flux in excess of 2 TW m-2, and provide sufficient heating to the GdFeCo electrons to induce deterministic reversal of the magnetic moment.


title: Finite dimensional thermo-mechanical systems and second order constraints

abstract: In this paper we study a class of physical systems that combine a finite number of mechanical and thermodynamic observables. We call them finite dimensional thermo-mechanical systems. We introduce these systems by means of simple examples. The evolution equations of the involved observables are obtained in each example by using, essentially, the Newton's law and the First Law of Thermodynamics only. We show that such equations are similar to those defining certain mechanical systems with higher order constraints. Moreover, we show that all of the given examples can be described in a variational formalism in terms of second order constrained systems.

categories: math-ph math.DG math.MP physics.class-ph


title: On the spatial distribution of neutral hydrogen in the Universe: bias and shot-noise of the HI Power Spectrum

abstract: The spatial distribution of neutral hydrogen (HI) in the Universe contains a wealth of cosmological information. The 21 cm emission line can be used to map the HI up to very high redshift and therefore reveal us something about the evolution of the large scale structures in the Universe. However little is known about the abundance and clustering properties of the HI over cosmic time. Motivated by this, we build an analytic framework where the relevant parameters that govern how the HI is distributed among dark matter halos can be fixed using observations. At the same time we provide tools to study the column density distribution function of the HI absorbers together with their clustering properties. Our formalism is the first one able to account for all observations at a single redshift, $z = 2.3$. The linear bias of the HI and the mean number density of HI sources, two main ingredients in the calculation of the signal-to-noise ratio of a cosmological survey, are then discussed in detail, also extrapolating the results to low and high redshift. We find that HI bias is relatively higher than the value reported in similar studies, but the shot noise level is always sub dominant, making the HI Power Spectrum always a high signal-to-noise measurements up to $z\simeq5$ in the limit of no instrumental noise and foreground contamination.


title: Real-Time Single Image and Video Super-Resolution Using an Efficient Sub-Pixel Convolutional Neural Network

abstract: Recently, several models based on deep neural networks have achieved great success in terms of both reconstruction accuracy and computational performance for single image super-resolution. In these methods, the low resolution (LR) input image is upscaled to the high resolution (HR) space using a single filter, commonly bicubic interpolation, before reconstruction. This means that the super-resolution (SR) operation is performed in HR space. We demonstrate that this is sub-optimal and adds computational complexity. In this paper, we present the first convolutional neural network (CNN) capable of real-time SR of 1080p videos on a single K2 GPU. To achieve this, we propose a novel CNN architecture where the feature maps are extracted in the LR space. In addition, we introduce an efficient sub-pixel convolution layer which learns an array of upscaling filters to upscale the final LR feature maps into the HR output. By doing so, we effectively replace the handcrafted bicubic filter in the SR pipeline with more complex upscaling filters specifically trained for each feature map, whilst also reducing the computational complexity of the overall SR operation. We evaluate the proposed approach using images and videos from publicly available datasets and show that it performs significantly better (+0.15dB on Images and +0.39dB on Videos) and is an order of magnitude faster than previous CNN-based methods.


title: Density distribution of a dust cloud in three-dimensional complex plasmas

abstract: We propose a novel method of determination of the dust particle spatial distribution in dust clouds that form in three-dimensional (3D) complex plasmas under microgravity conditions. The method utilizes the data obtained during the 3D scanning of a cloud and provides a reasonably good accuracy. Based on this method, we investigate the particle density in a dust cloud realized in gas discharge plasma in the PK-3 Plus setup onboard the International Space Station. We find that the treated dust clouds are both anisotropic and inhomogeneous. One can isolate two regimes, in which a stationary dust cloud can be observed. At low pressures, the particle density decreases monotonically with the increase of the distance from the discharge center; at higher pressures, the density distribution has a shallow minimum. Regardless of the regime, we detect a cusp of the distribution at the void boundary and a slowly varying density at larger distances (in the foot region). A theoretical interpretation of obtained results is developed that leads to reasonable estimates of the densities for both the cusp and foot. The modified ionization equation of state, which allows for violation of the local quasineutrality in the cusp region, predicts the spatial distributions of ion and electron densities to be measured in future experiments.


title: Energy-Efficient Resource Allocation for SWIPT in Multiple Access Channels

abstract: In this paper, we study optimal resource allocation strategies for simultaneous information and power transfer (SWIPT) focusing on the system energy efficiency. We consider two-user multiple access channels in which energy harvesting (EH) and information decoding (ID) nodes are spatially separated. We formulate optimization problems that maximize system energy efficiency while taking harvested energy constraints into account. These are concave-linear fractional problems, and hence Karush-Kuhn-Tucker (KKT) conditions are necessary and sufficient to obtain globally optimal solution. Solving these optimization problems, we provide analytical expressions for optimal transmit power allocation among the source nodes, and identify the corresponding energy efficiency. We confirm the theoretical analysis via numerical results. Furthermore, we also characterize the effect of circuit power consumption on the system's efficiency as the harvested energy demand varies.


title: Rational Whitney tower filtration of links

abstract: We present complete classifications of links in the 3-sphere modulo framed and twisted Whitney towers in a rational homology 4-ball. This provides a geometric characterization of the vanishing of the Milnor invariants of links in terms of Whitney towers. Our result also says that the higher order Arf invariants, which are conjectured to be nontrivial, measure the potential difference between the Whitney tower theory in rational homology 4-balls and that in the 4-ball extensively developed by Conant, Schneiderman and Teichner.


title: No-Regret Replanning under Uncertainty

abstract: This paper explores the problem of path planning under uncertainty. Specifically, we consider online receding horizon based planners that need to operate in a latent environment where the latent information can be modeled via Gaussian Processes. Online path planning in latent environments is challenging since the robot needs to explore the environment to get a more accurate model of latent information for better planning later and also achieves the task as quick as possible. We propose UCB style algorithms that are popular in the bandit settings and show how those analyses can be adapted to the online robotic path planning problems. The proposed algorithm trades-off exploration and exploitation in near-optimal manner and has appealing no-regret properties. We demonstrate the efficacy of the framework on the application of aircraft flight path planning when the winds are partially observed.


title: Quantum Thermal Rectification to design thermal diodes and transistors

abstract: We study in this article how heat can be exchanged between two level systems (TLS) each of them being coupled to a thermal reservoir. Calculation are performed solving a master equation for the density matrix using the Born markov-approximation. We analyse the conditions for which a thermal diode and a thermal transistor can be obtained as well as their optimization.


title: Dispersive estimates for Dirac Operators in dimension three with obstructions at threshold energies

abstract: We investigate $L^1\to L^\infty$ dispersive estimates for the three dimensional Dirac equation with a potential. We also classify the structure of obstructions at the thresholds of the essential spectrum as being composed of a two dimensional space of resonances and finitely many eigenfunctions. We show that, as in the case of the Schr\"odinger evolution, the presence of a threshold obstruction generically leads to a loss of the natural $t^{-\frac32}$ decay rate. In this case we show that the solution operator is composed of a finite rank operator that decays at the rate $t^{-\frac12}$ plus a term that decays at the rate $t^{-\frac32}$.


title: ADC Bit Allocation under a Power Constraint for MmWave Massive MIMO Communication Receivers

abstract: Millimeter wave (mmWave) systems operating over a wide bandwidth and using a large number of antennas impose a heavy burden on power consumption. In a massive multiple-input multiple-output (MIMO) uplink, analog-to-digital con- verters (ADCs) would be the primary consumer of power in the base station receiver. This paper proposes a bit allocation (BA) method for mmWave multi-user (MU) massive MIMO systems under a power constraint. We apply ADCs to the outputs of an analog phased array for beamspace projection to exploit mmWave channel sparsity. We relax a mean square quantization error (MSQE) minimization problem and map the closed-form solution to non-negative integer bits at each ADC. In link-level simulations, the proposed method gives better communication performance than conventional low-resolution ADCs for the same or less power. Our contribution is a near optimal low-complexity BA method that minimizes total MSQE under a power constraint.


title: On fixed gain recursive estimators with discontinuity in the parameters

abstract: In this paper we estimate the tracking error of a fixed gain stochastic approximation scheme. The underlying process is not assumed Markovian, a mixing condition is required instead. Furthermore, the updating function may be discontinuous in the parameter.


title: Improving the semidefinite programming bound for the kissing number by exploiting polynomial symmetry

abstract: The kissing number of $\mathbb{R}^n$ is the maximum number of pairwise-nonoverlapping unit spheres that can simultaneously touch a central unit sphere. Mittelmann and Vallentin (2010), based on the semidefinite programming bound of Bachoc and Vallentin (2008), computed the best known upper bounds for the kissing number for several values of $n \leq 23$. In this paper, we exploit the symmetry present in the semidefinite programming bound to provide improved upper bounds for $n = 9, \ldots, 23$.


title: Lasing Action with Gold Nanorod Hyperbolic Metamaterials

abstract: Coherent nanoscale photon sources are of paramount importance to achieving all-optical communication. Several nanolasers smaller than the diffraction limit have been theoretically proposed and experimentally demonstrated using plasmonic cavities to confine optical fields. Such compact cavities exhibit large Purcell factors, thereby enhancing spontaneous emission, which feeds into the lasing mode. However, most plasmonic nanolasers reported so far have employed resonant nanostructures and therefore had the lasing restricted to the proximity of the resonance wavelength. Here, we report on an approach based on gold nanorod hyperbolic metamaterials for lasing. Hyperbolic metamaterials provide broadband Purcell enhancement due to large photonic density of optical states, while also supporting surface plasmon modes to deliver optical feedback for lasing due to nonlocal effects in nanorod media. We experimentally demonstrate the advantage of hyperbolic metamaterials in achieving lasing action by its comparison with that obtained in a metamaterial with elliptic dispersion. The conclusions from the experimental results are supported with numerical simulations comparing the Purcell factors and surface plasmon modes for the metamaterials with different dispersions. We show that although the metamaterials of both types support lasing, emission with hyperbolic samples is about twice as strong with 35% lower threshold vs. the elliptic ones. Hence, hyperbolic metamaterials can serve as a convenient platform of choice for nanoscale coherent photon sources in a broad wavelength range.


title: An integral Relationship for a new Fractional One-phase Stefan Problem

abstract: A one-dimensional fractional one-phase Stefan problem with a temperature boundary condition at the fixed face is considered. An integral relationship between the temperature and the free boundary is obtained which is equivalent to the fractional Stefan condition. Moreover, an exact solution of similarity type expressed in terms of Wright functions is given.


title: Should Terminology Principles be re-examined?

abstract: Operationalization of terminology for IT applications has revived the Wusterian approach. The conceptual dimension once more prevails after taking back seat to specialised lexicography. This is demonstrated by the emergence of ontology in terminology. While the Terminology Principles as defined in Felber manual and the ISO standards remain at the core of traditional terminology, their computational implementation raises some issues. In this article, while reiterating their importance, we will be re-examining these Principles from a dual perspective: that of logic in the mathematical sense of the term and that of epistemology as in the theory of knowledge. We will thus be clarifying and describing some of them so as to take into account advances in knowledge engineering (ontology) and formal systems (logic). The notion of ontoterminology, terminology whose conceptual system is a formal ontology, results from this approach.


title: Cumulants and nonlinear response of high $p_T$ harmonic flow at $\sqrt{s_{NN}}=5.02$ TeV

abstract: Event-by-event fluctuations caused by quantum mechanical fluctuations in the wave function of colliding nuclei in ultrarelativistic heavy ion collisions were recently shown to be necessary for the simultaneous description of $R_{AA}$ as well as the elliptic and triangular flow harmonics at high $p_T$ in PbPb collisions at the Large Hadron Collider. In fact, the presence of a finite triangular flow as well as cumulants of the flow harmonic distribution that differ from the mean are only possible when these event-by-event fluctuations are considered. In this paper we combine event-by-event viscous hydrodynamics and jet quenching to make predictions for high $p_T$ $R_{AA}$, $v_2\{2\}$, $v_3\{2\}$, and $v_2\{4\}$ in PbPb collisions at $\sqrt{s_{NN}}=5.02$ TeV. With an order of magnitude larger statistics we find that high $p_T$ elliptic flow does not scale linearly with the soft elliptical flow, as originally thought, but has deviations from perfectly linear scaling. A new experimental observable, which involves the difference between the ratio of harmonic flow cumulants at high and low $p_T$, is proposed to investigate the fluctuations of high $p_T$ flow harmonics and measure this nonlinear response. By varying the path length dependence of the energy loss and the viscosity of the evolving medium we find that $R_{AA}(p_T)$ and $v_2\{2\}(p_T)$ strongly depend on the choice for the path length dependence of the energy loss, which can be constrained using the new LHC run 2 data.


title: Long term stability and reproducibility of magnetic colloids are key issues for steady values of Specific Power Absorption through time

abstract: Virtually all clinical applications of magnetic nanoparticles (MNPs) require the formulation of biocompatible, water-based magnetic colloids. For magnetic hyperthermia, the requirements also include a high colloidal stability against precipitation and agglomeration of the constituent MNPs, in order to keep the heating efficiency of the ferrofluid in the long term. The specific power absorption (SPA) of single-domain MNPs depends critically on the average particle size and size distribution width, therefore first-rate reproducibility among different batches regarding these parameters are also needed. We have studied the evolution of the SPA of highly reproducible and stable water-based colloid composed of polymer coated $Fe_{3}O_{4}$ magnetic nanoparticles. By measuring the specific power absorption (SPA) values along one year as a function of field amplitude and frequencies ($H \leq 24 kA/m$; $260 \leq f \leq 830 kHz$), we demonstrated that SPA in these samples can be made reproducible between successive synthesis, and stable along several months, due to the in situ polymer coating that provides colloidal stability and keeps dipolar interactions negligible.


title: Simulating device-to-device communications in OMNeT++ with SimuLTE: scenarios and configurations

abstract: SimuLTE is a tool that enables system-level simulations of LTE/LTE-Advanced networks within OMNeT++. It is designed such that it can be plugged within network elements as an additional Network Interface Card (NIC) to those already provided by the INET framework (e.g. Wi-Fi). Recently, device-to-device (D2D) technology has been widely studied by the research community, as a mechanism to allow direct communications between devices of a LTE cellular network. In this work, we present how SimuLTE can be employed to simulate both one-to-one and one-to-many D2D communications, so that the latter can be exploited as a new communication opportunity in several research fields, like vehicular networks, IoT and machine-to-machine (M2M) applications.


title: On groups, slow heat kernel decay yields Liouville property and sharp entropy bounds

abstract: Let $\mu$ be a symmetric probability measure of finite entropy on a group $G$. We show that if $-\log \mu^{(2n)}(id)=o(n^{1/2})$, then the pair $(G,\mu)$ has the Liouville property (all bounded $\mu$-harmonic functions on $G$ are constant). Furthermore, if $-\log \mu^{(2n)}(id)=O(n^{\beta})$ where $\beta\in(0,1/2)$, then the entropy of the $n$-fold convolution power $\mu^{(n)}$ satisfies $H(\mu^{(n)})=O\left(n^{\frac{\beta}{1-\beta}}\right)$. This improves earlier results of Gournay and of Saloff-Coste and the second author. We extend the bounds to transitive graphs and illustrate their sharpness on a family of groups.


title: A test of Gaia Data Release 1 parallaxes: implications for the local distance scale

abstract: We present a comparison of Gaia Data Release 1 (DR1) parallaxes with photometric parallaxes for a sample of 212 Galactic Cepheids at a median distance of 2~kpc, and explore their implications on the distance scale and the local value of the Hubble constant H_0. The Cepheid distances are estimated from a recent calibration of the near-infrared Period-Luminosity P-L relation. The comparison is carried out in parallax space, where the DR1 parallax errors, with a median value of half the median parallax, are expected to be well-behaved. With the exception of one outlier, the DR1 parallaxes are in remarkably good global agreement with the predictions, and the published errors may be conservatively overestimated by about 20%. The parallaxes of 9 Cepheids brighter than G = 6 may be systematically underestimated, trigonometric parallaxes measured with the HST FGS for three of these objects confirm this trend. If interpreted as an independent calibration of the Cepheid luminosities and assumed to be otherwise free of systematic uncertainties, DR1 parallaxes would imply a decrease of 0.3% in the current estimate of the local Hubble constant, well within their statistical uncertainty, and corresponding to a value 2.5 sigma (3.5 sigma if the errors are scaled) higher than the value inferred from Planck CMB data used in conjunction with Lambda-CDM. We also test for a zeropoint error in Gaia parallaxes and find none to a precision of ~20 muas. We caution however that with this early release, the complete systematic properties of the measurements may not be fully understood at the statistical level of the Cepheid sample mean, a level an order of magnitude below the individual uncertainties. The early results from DR1 demonstrate again the enormous impact that the full mission will likely have on fundamental questions in astrophysics and cosmology.


title: Components of affine Springer fibers

abstract: Let $\mathbf{G}$ be a connected split reductive group over a field of characteristic zero or sufficiently large characteristic, $\gamma_0\in(\operatorname{Lie}\mathbf{G})((t))$ be any topologically nilpotent regular semisimple element, and $\gamma=t\gamma_0$. Using methods from $p$-adic orbital integrals, we show that the number of components of the Iwahori affine Springer fiber over $\gamma$ modulo $Z_{\mathbf{G}((t))}(\gamma)$ is equal to the order of the Weyl group.


title: The microstructural foundations of leverage effect and rough volatility

abstract: We show that typical behaviors of market participants at the high frequency scale generate leverage effect and rough volatility. To do so, we build a simple microscopic model for the price of an asset based on Hawkes processes. We encode in this model some of the main features of market microstructure in the context of high frequency trading: high degree of endogeneity of market, no-arbitrage property, buying/selling asymmetry and presence of metaorders. We prove that when the first three of these stylized facts are considered within the framework of our microscopic model, it behaves in the long run as a Heston stochastic volatility model, where leverage effect is generated. Adding the last property enables us to obtain a rough Heston model in the limit, exhibiting both leverage effect and rough volatility. Hence we show that at least part of the foundations of leverage effect and rough volatility can be found in the microstructure of the asset.


title: Privacy Preserving Distance Computation using Somewhat-trusted Third Parties

abstract: A critically important component of most signal processing procedures is that of computing the distance between signals. In multi-party processing applications where these signals belong to different parties, this introduces privacy challenges. The signals may themselves be private, and the parties to the computation may not be willing to expose them. Solutions proposed to the problem in the literature generally invoke homomorphic encryption schemes, secure multi-party computation, or other cryptographic methods which introduce significant computational complexity into the proceedings, often to the point of making more complex computations requiring repeated computations unfeasible. Other solutions invoke third parties, making unrealistic assumptions about their trustworthiness. In this paper we propose an alternate approach, also based on third party computation, but without assuming as much trust in the third party. Individual participants to the computation "secure" their data through a proposed secure hashing scheme with shared keys, prior to sharing it with the third party. The hashing ensures that the third party cannot recover any information about the individual signals or their statistics, either from analysis of individual computations or their long-term aggregate patterns. We provide theoretical proof of these properties and empirical demonstration of the feasibility of the computation.


title: Extending OMNeT++ Towards a Platform for the Design of Future In-Vehicle Network Architectures

abstract: In-vehicle communication technologies are evolving. While today's cars are equipped with fieldbusses to interconnect the various electronic control units, next generation vehicles have timing and bandwidth requirements that exceed the capacities. In particular Advanced Driver Assistance Systems (ADAS) and automated driving using high bandwidth sensors such as cameras, LIDAR or radar will challenge the in-car network. Automotive Ethernet is the most promising candidate to solve the upcoming challenges. But to design and evaluate new protocols, concepts, and architectures suitable analysis tools are required. Especially in the interim period with architectures using automotive Ethernet and legacy fieldbusses together, careful planning and design is of vital importance. Simulation can provide a good understanding of the expectable network metrics in an early development phase. This paper contributes a workflow as well as the required toolchain to evaluate new real-time Ethernet communication architectures using event based simulation in OMNeT++. We introduce a domain specific language (DSL) - the Abstract Network Description Language (ANDL) - to describe and configure the simulation and present the required simulation models for real-time Ethernet and fieldbus technologies such as CAN and FlexRay. We further introduce new analysis tools for special in-vehicle network use-cases and the interaction of the simulation with third-party applications established in the automotive domain.


title: Grammatical Templates: Improving Text Difficulty Evaluation for Language Learners

abstract: Language students are most engaged while reading texts at an appropriate difficulty level. However, existing methods of evaluating text difficulty focus mainly on vocabulary and do not prioritize grammatical features, hence they do not work well for language learners with limited knowledge of grammar. In this paper, we introduce grammatical templates, the expert-identified units of grammar that students learn from class, as an important feature of text difficulty evaluation. Experimental classification results show that grammatical template features significantly improve text difficulty prediction accuracy over baseline readability features by 7.4%. Moreover, we build a simple and human-understandable text difficulty evaluation approach with 87.7% accuracy, using only 5 grammatical template features.


title: Information Theoretic Limits of Data Shuffling for Distributed Learning

abstract: Data shuffling is one of the fundamental building blocks for distributed learning algorithms, that increases the statistical gain for each step of the learning process. In each iteration, different shuffled data points are assigned by a central node to a distributed set of workers to perform local computations, which leads to communication bottlenecks. The focus of this paper is on formalizing and understanding the fundamental information-theoretic trade-off between storage (per worker) and the worst-case communication overhead for the data shuffling problem. We completely characterize the information theoretic trade-off for $K=2$, and $K=3$ workers, for any value of storage capacity, and show that increasing the storage across workers can reduce the communication overhead by leveraging coding. We propose a novel and systematic data delivery and storage update strategy for each data shuffle iteration, which preserves the structural properties of the storage across the workers, and aids in minimizing the communication overhead in subsequent data shuffling iterations.


title: Emission of Photons and Relativistic Axions from Axion Stars

abstract: The number of nonrelativistic axions can be changed by inelastic reactions that produce photons or relativistic axions. Any odd number of axions can annihilate into two photons. Any even number of nonrelativistic axions can scatter into two relativistic axions. We calculate the rate at which axions are lost from axion stars from these inelastic reactions. In dilute systems of axions, the dominant inelastic reaction is axion decay into two photons. In sufficiently dense systems of axions, the dominant inelastic reaction is the scattering of four nonrelativistic axions into two relativistic axions. The scattering of odd numbers of axions into two photons produces monochromatic radio-frequency signals at odd-integer harmonics of the fundamental frequency set by the axion mass. This provides a unique signature for dense systems of axions, such as a dense axion star or a collapsing dilute axion star.


title: Relativistic Aberration and the Cosmological Constant in Gravitational Lensing I: Introduction

abstract: An analysis of null geodesics in Schwarzschild de Sitter space is presented with special attention to their global `bending angles', local measurable angles, and the involvement of the cosmological constant. We make use of a general technique which allows for finding observable intersection angles of null trajectories analytically. A general relativistic aberration relationship is established as one of its applications. The question of whether or not the cosmological constant, $\Lambda$, contributes to orbits of light and to related observable quantities is addressed in detail. We also discuss the ongoing debate on this issue and respond to some recent papers on the topic. The dependence of measurable quantities on the motion of observers is stressed throughout. Exact formulas for measurable intersection angles, as well as gravitational lens equations for observers in the Schwarzschild de Sitter background are provided.


title: Estimation of Sobolev embedding constant on a domain dividable into bounded convex domains

abstract: This paper is concerned with an explicit value of the embedding constant from $W^{1,q}(\Omega)$ to $L^{p}(\Omega)$ for a bounded domain $\Omega\subset\mathbb{R}^N~(N\in\mathbb{N})$, where $1\leq q\leq p\leq \infty$. To obtain this value, we previously proposed a formula for estimating the embedding constant on bounded and unbounded Lipschitz domains by estimating the norm of Stein's extension operator, in the article (K. Tanaka, K. Sekine, M. Mizuguchi, and S. Oishi, Estimation of Sobolev-type embedding constant on domains with minimally smooth boundary using extension operator, Journal of Inequalities and Applications, Vol. 389, pp. 1-23, 2015). This formula is also applicable to a domain that can be divided into Lipschitz domains. However, the values computed by the previous formula are very large. In this paper, we propose several sharper estimations of the embedding constant on a bounded domain that can be divided into convex domains.


title: The wild monodromy of the Fifth Painlev\'e equation and its action on wild character variety: an approach of confluence

abstract: The article studies the Fifth Painlev\'e equation and of the nonlinear Stokes phenomenon at its irregular singularity at infinity from the point of view of confluence from the Sixth Painlev\'e equation. This approach is developped separately on both sides of the Riemann-Hilbert correspondance. On the side of the nonlinear Painlev\'e-Okamoto foliation the relation between the nonlinear monodromy group of Painlev\'e VI and the "nonlinear wild monodromy pseudogroup" of Painlev\'e V (that is the pseudogroup generated by nonlinear monodromy, nonlinear Stokes operators and nonlinear exponential torus) is explained in detail. On the side of the corresponding linear isomonodromic problem, the "wild" character variety (the space of the linear monodromy and Stokes data) associated to Painlev\'e V is constructed through a birational transformation from the character variety (the space of the linear monodromy data) associated to Painlev\'e VI. This allows to transport the known description of the action of the nonlinear monodromy of Painlev\'e VI on its character variety to that of Painlev\'e V, and to provide explicit formulas for the action of the "nonlinear wild monodromy" of Painlev\'e V on its character variety.


title: Spatial Multiresolution Analysis of the Effect of PM2.5 on Birth Weights

abstract: Fine particulate matter (PM2.5) measured at a given location is a mix of pollution generated locally and pollution traveling long distances in the atmosphere. Therefore, the identification of spatial scales associated with health effects can inform on pollution sources responsible for these effects, resulting in more targeted regulatory policy. Recently, prediction methods that yield high-resolution spatial estimates of PM2.5 exposures allow one to evaluate such scale-specific associations. We propose a two-dimensional wavelet decomposition that alleviates restrictive assumptions required for standard wavelet decompositions. Using this method we decompose daily surfaces of PM2.5 to identify which scales of pollution are most associated with adverse health outcomes. A key feature of the approach is that it can remove the purely temporal component of variability in PM2.5 levels and calculate effect estimates derived solely from spatial contrasts. This eliminates the potential for unmeasured confounding of the exposure - outcome associations by temporal factors, such as season. We apply our method to a study of birth weights in Massachusetts, U.S.A from 2003-2008 and find that both local and urban sources of pollution are strongly negatively associated with birth weight. Results also suggest that failure to eliminate temporal confounding in previous analyses attenuated the overall effect estimate towards zero, with the effect estimate growing in magnitude once this source of variability is removed.


title: Optimizing isotope and vacancy engineering in graphene ribbons to enhance the thermoelectric performance without degrading the electronic properties

abstract: The enhancement of thermoelectric figure of merit ZT requires to either increase the power factor or reduce the phonon conductance, or even both. In graphene, the high phonon thermal conductivity is the main factor limiting the thermoelectric conversion. The common strategy to enhance ZT is therefore to introduce phonon scatterers to suppress the phonon conductance while retaining high electrical conductance and Seebeck coefficient. Although thermoelectric performance is eventually enhanced, all studies based on this strategy show a significant reduction of the electrical conductance, most often leading to a lower electronic performance. In this study we show that appropriate sources of disorder, including isotopes and vacancies at lowest electron density positions, can be used as phonon scatterers to reduce the phonon conductance in graphene ribbons without degrading the electrical conductance, particularly in the low-energy region which is the most important range for device operation. By means of atomistic calculations using semi-empirical Tight-Binding and Force Constant models in combination with Non-Equilibrium Green function formalism, we show that the natural electronic properties of graphene ribbons can be fully preserved while their thermoelectric efficiency is strongly enhanced. For ribbons of width M = 5 dimer lines, room-temperature ZT is enhanced from less than 0.26 for defect-free ribbons to more than 2.5. This study is likely to set the milestones of a new generation of nano-devices with dual electronic, thermoelectric functionalities.


title: Mode testing, critical bandwidth and excess mass

abstract: The identification of peaks or maxima in probability densities, by mode testing or bump hunting, has become an important problem in applied fields. This task has been approached in the statistical literature from different perspectives, with the proposal of testing procedures which are based on kernel density estimators or on the quantification of excess mass. However, none of the existing proposals provides a satisfactory performance in practice. In this work, a new procedure which combines the previous approaches (smoothing and excess mass) is presented and compared with the existing methods, showing a superior behaviour. A real data example on philatelic data is also included for illustration purposes.


title: Higher order selfdual toric varieties

abstract: The notion of higher order dual varieties of a projective variety, introduced in \cite{P83}, is a natural generalization of the classical notion of projective duality. In this paper we present geometric and combinatorial characterizations of those equivariant projective toric embeddings that satisfy higher order selfduality. We also give several examples and general constructions. In particular, we highlight the relation with Cayley-Bacharach questions and with Cayley configurations.


title: Theory News Higgs

abstract: I review new developments in Higgs physics, with a focus on Yukawa couplings in and beyond the standard model. In particular, I discuss different methods of measuring the light Yukawas, new sources of CP violation in the Higgs sector, and lepton flavor violation.


title: Gradient Descent Learns Linear Dynamical Systems

abstract: We prove that stochastic gradient descent efficiently converges to the global optimizer of the maximum likelihood objective of an unknown linear time-invariant dynamical system from a sequence of noisy observations generated by the system. Even though the objective function is non-convex, we provide polynomial running time and sample complexity bounds under strong but natural assumptions. Linear systems identification has been studied for many decades, yet, to the best of our knowledge, these are the first polynomial guarantees for the problem we consider.


title: Extended X-ray Monitoring of Gravitational Lenses with Chandra and Joint Constraints on X-ray Emission Regions

abstract: We present an X-ray photometric analysis of six gravitationally lensed quasars, with observation campaigns spanning from 5 to 14 years, measuring the total (0.83 - 21.8 keV restframe), soft (0.83 - 3.6 keV), and hard (3.6 - 21.8 keV) band image flux ratios or each epoch. Using the ratios of the model-predicted macro-magnifications as baselines, we build differential microlensing light curves and obtain joint likelihood functions for the average X-ray emission region sizes. Our analysis yields a Probability Distribution Function for the average half-light radius of the X-Ray emission region in the sample that peaks slightly above 1 gravitational radius and with nearly indistinguishable 68% confidence (one-sided) upper limits of 17.8 and 18.9 gravitational radii for the soft and hard X-ray emitting regions, assuming a mean stellar mass of 0.3 solar masses. We see hints of energy dependent microlensing between the soft and hard bands in two of the objects. In a separate analysis on the root-mean-square of the microlensing variability, we find significant differences between the soft and hard bands but the sign of the difference is not consistent across the sample. This suggests the existence of some kind of spatial structure to the X-ray emission in an otherwise extremely compact source. We also discover a correlation between the RMS microlensing variability and the average microlensing amplitude.


title: Casimir Force Phase Transitions in the Graphene Family

abstract: The Casimir force is a universal interaction induced by electromagnetic quantum fluctuations between any types of objects. The expansion of the graphene family by adding silicene, germanene, and stanene, 2D allotropes of Si, Ge, and Sn, lands itself as a platform to probe Dirac-like physics in honeycomb staggered systems in such a ubiquitous interaction. We discover Casimir force phase transitions between these staggered 2D materials induced by the complex interplay between Dirac physics, spin-orbit coupling, and externally applied fields. In particular, we find that the interaction energy experiences different power law distance decays, magnitudes, and dependences on characteristic physical constants. Furthermore, due to the topological properties of these materials, repulsive and quantized Casimir interactions become possible.


title: The Bradley--Terry condition is $L_1$--testable

abstract: We provide an algorithm with constant running time that given a weighted tournament $T$, distinguishes with high probability of success between the cases that $T$ can be represented by a Bradley--Terry model, or cannot even be approximated by one. The same algorithm tests whether the corresponding Markov chain is reversible.


title: Mathematical Properties of Numerical Inversion for Jet Calibrations

abstract: Numerical inversion is a general detector calibration technique that is independent of the underlying spectrum. This procedure is formalized and important statistical properties are presented, using high energy jets at the Large Hadron Collider as an example setting. In particular, numerical inversion is inherently biased and common approximations to the calibrated jet energy tend to over-estimate the resolution. Analytic approximations to the closure and calibrated resolutions are demonstrated to effectively predict the full forms under realistic conditions. Finally, extensions of numerical inversion are presented which can reduce the inherent biases. These methods will be increasingly important to consider with degraded resolution at low jet energies due to a much higher instantaneous luminosity in the near future.


abstract: We study Jordan-Lie inner ideals of finite dimensional associative algebras and the corresponding Lie algebras and prove that they admit Levi decompositions. Moreover, we classify Jordan-Lie inner ideals satisfying a certain minimality condition and show that they are generated by pairs of idempotents.


title: Parametric Resonance in the Early Universe - A Fitting Analysis

abstract: Particle production via parametric resonance in the early Universe, is a nonperturbative, non-linear and out-of-equilibrium phenomenon. Although it is a well studied topic, whenever a new scenario exhibits parametric resonance, a full re-analysis is normally required. To avoid this tedious task, many works present often only a simplified linear treatment of the problem. In order to surpass this circumstance in the future, we provide a fitting analysis of parametric resonance through all its relevant stages: initial linear growth, non-linear evolution, and relaxation towards equilibrium. Using lattice simulations in an expanding grid in 3 + 1 dimensions, we parametrize the dynamics outcome scanning over the relevant ingredients: role of the oscillatory field, particle coupling strength, initial conditions, and background expansion rate. We emphasize the inaccuracy of the linear calculation of the decay time of the oscillatory field, and propose a more appropriate definition of this scale based on the subsequent non-linear dynamics. We provide simple fits to the relevant time scales and particle energy fractions at each stage. Our fits can be applied to post-inflationary preheating scenarios, where the oscillatory field is the inflaton, or to spectator-field scenarios, where the oscillatory field can be e.g. a curvaton, or the Standard Model Higgs.