title: A high-order wideband direct solver for electromagnetic scattering from bodies of revolution

abstract: The generalized Debye source representation of time-harmonic electromagnetic fields yields well-conditioned second-kind integral equations for a variety of boundary value problems, including the problems of scattering from perfect electric conductors and dielectric bodies. Furthermore, these representations, and resulting integral equations, are fully stable in the static limit as $\omega \to 0$ in multiply connected geometries. In this paper, we present the first high-order accurate solver based on this representation for bodies of revolution. The resulting solver uses a Nystr\"om discretization of a one-dimensional generating curve and high-order integral equation methods for applying and inverting surface differentials. The accuracy and speed of the solvers are demonstrated in several numerical examples.


title: Dualism between Optical and Difference Parametric Amplification

abstract: Breaking the symmetry in a coupled wave system can result in unusual amplification behavior. In the case of difference parametric amplification the resonant pump frequency is equal to the difference, instead of the sum, frequency of the normal modes. We show that sign reversal in the symmetry relation of parametric coupling give rise to difference parametric amplification as a dual of optical parametric amplification. For optical systems, our result can potentially be used for efficient XUV amplification.


title: Lectures on the Spin and Loop $O(n)$ Models

abstract: The classical spin $O(n)$ model is a model on a $d$-dimensional lattice in which a vector on the $(n-1)$-dimensional sphere is assigned to every lattice site and the vectors at adjacent sites interact ferromagnetically via their inner product. Special cases include the Ising model ($n=1$), the XY model ($n=2$) and the Heisenberg model ($n=3$). We discuss questions of long-range order and decay of correlations in the spin $O(n)$ model for different combinations of the lattice dimension $d$ and the number of spin components $n$. The loop $O(n)$ model is a model for a random configuration of disjoint loops. We discuss its properties on the hexagonal lattice. The model is parameterized by a loop weight $n\ge0$ and an edge weight $x\ge 0$. Special cases include self-avoiding walk ($n=0$), the Ising model ($n=1$), critical percolation ($n=x=1$), dimer model ($n=1,x=\infty$), proper $4$-coloring ($n=2, x=\infty)$, integer-valued ($n=2$) and tree-valued (integer $n>=3$) Lipschitz functions and the hard hexagon model ($n=\infty$). The object of study in the model is the typical structure of loops. We review the connection of the model with the spin $O(n)$ model and discuss its conjectured phase diagram, emphasizing the many open problems remaining.


title: Asymptotically optimal private estimation under mean square loss

abstract: We consider the minimax estimation problem of a discrete distribution with support size $k$ under locally differential privacy constraints. A privatization scheme is applied to each raw sample independently, and we need to estimate the distribution of the raw samples from the privatized samples. A positive number $\epsilon$ measures the privacy level of a privatization scheme. In our previous work (arXiv:1702.00610), we proposed a family of new privatization schemes and the corresponding estimator. We also proved that our scheme and estimator are order optimal in the regime $e^{\epsilon} \ll k$ under both $\ell_2^2$ and $\ell_1$ loss. In other words, for a large number of samples the worst-case estimation loss of our scheme was shown to differ from the optimal value by at most a constant factor. In this paper, we eliminate this gap by showing asymptotic optimality of the proposed scheme and estimator under the $\ell_2^2$ (mean square) loss. More precisely, we show that for any $k$ and $\epsilon,$ the ratio between the worst-case estimation loss of our scheme and the optimal value approaches $1$ as the number of samples tends to infinity.


title: Application of Bayesian Networks for Estimation of Individual Psychological Characteristics

abstract: An accurate qualitative and comprehensive assessment of human potential is one of the most important challenges in any company or collective. We apply Bayesian networks for developing more accurate overall estimations of psychological characteristics of an individual, based on psychological test results, which identify how much an individual possesses a certain trait. Examples of traits could be a stress resistance, the readiness to take a risk, the ability to concentrate on certain complicated work. The most common way of studying psychological characteristics of each individual is testing. Additionally, the overall estimation is usually based on personal experiences and the subjective perception of a psychologist or a group of psychologists about the investigated psychological personality traits.


title: A Seifert algorithm for lamination links

abstract: We generalize H. Seifert's algorithm for finding a Seifert surface for a knot or link. The generalization applies to "framed oriented measured lamination links." For knots, a Seifert surface determines a unique framing. In our setting, we analyze the set of framed lamination links which bound Seifert laminations and are carried by an $S^1$-fibered tube neighborhood of an oriented train track embedded in a 3-manifold.


title: Sparse Structural Approach for Rating Transitions

abstract: In banking practice, rating transition matrices have become the standard approach of deriving multi-year probabilities of default (PDs) from one-year PDs, the latter normally being available from Basel ratings. Rating transition matrices have gained in importance with the newly adopted IFRS 9 accounting standard. Here, the multi-year PDs can be used to calculate the so-called expected credit losses (ECL) over the entire lifetime of relevant credit assets. A typical approach for estimating the rating transition matrices relies on calculating empirical rating migration counts and frequencies from rating history data. For small portfolios, however, this approach often leads to zero counts and high count volatility, which makes the estimations unreliable and unstable, and can also produce counter-intuitive prediction patterns such as non-parallel/crossing forward PD patterns. This paper proposes a structural model which overcomes these problems. We make a plausible assumption of an underlying autoregressive mean-reverting ability-to-pay process. With only three parameters, this sparse process can well describe an entire typical rating transition matrix, provided the one-year PDs of the rating classes are specified. The transition probabilities produced by the structural approach are well-behaved by design. The approach significantly reduces the statistical degrees of freedom of the estimated transition probabilities, which makes the rating transition matrix more reliable for small portfolios. The approach can be applied to data with as few as 50 observed rating transitions. Moreover, the approach can be efficiently applied to data consisting of continuous PDs (prior to rating discretization). In the IFRS 9 context, the approach offers an additional merit: it can easily account for the macroeconomic adjustments, which are required by the IFRS 9 accounting standard.


title: Compressibility anomalies in stretched water and their interplay with density anomalies

abstract: Water keeps puzzling scientists because of its numerous properties which behave oppositely to usual liquids: for instance, water expands upon cooling, and liquid water is denser than ice. To explain this anomalous behaviour, several theories have been proposed, with different predictions for the properties of supercooled water (liquid at conditions where ice is more stable). However, discriminating between those theories with experiments has remained elusive because of spontaneous ice nucleation. Here we measure the sound velocity in liquid water stretched to negative pressure, and derive an experimental equation of state, which reveals compressibility anomalies. We show by rigorous thermodynamic relations how these anomalies are intricately linked with the density anomaly. Some features we observe are necessary conditions for the validity of two theories of water.


title: The inverse eigenvalue problem of a graph: Multiplicities and minors

abstract: The inverse eigenvalue problem of a given graph $G$ is to determine all possible spectra of real symmetric matrices whose off-diagonal entries are governed by the adjacencies in $G$. Barrett et al. introduced the Strong Spectral Property (SSP) and the Strong Multiplicity Property (SMP) in [8]. In that paper it was shown that if a graph has a matrix with the SSP (or the SMP) then a supergraph has a matrix with the same spectrum (or ordered multiplicity list) augmented with simple eigenvalues if necessary, that is, subgraph monotonicity. In this paper we extend this to a form of minor monotonicity, with restrictions on where the new eigenvalues appear. These ideas are applied to solve the inverse eigenvalue problem for all graphs of order five, and to characterize forbidden minors of graphs having at most one multiple eigenvalue.


title: Time-Dependent Representation for Neural Event Sequence Prediction

abstract: Existing sequence prediction methods are mostly concerned with time-independent sequences, in which the actual time span between events is irrelevant and the distance between events is simply the difference between their order positions in the sequence. While this time-independent view of sequences is applicable for data such as natural languages, e.g., dealing with words in a sentence, it is inappropriate and inefficient for many real world events that are observed and collected at unequally spaced points of time as they naturally arise, e.g., when a person goes to a grocery store or makes a phone call. The time span between events can carry important information about the sequence dependence of human behaviors. In this work, we propose a set of methods for using time in sequence prediction. Because neural sequence models such as RNN are more amenable for handling token-like input, we propose two methods for time-dependent event representation, based on the intuition on how time is tokenized in everyday life and previous work on embedding contextualization. We also introduce two methods for using next event duration as regularization for training a sequence prediction model. We discuss these methods based on recurrent neural nets. We evaluate these methods as well as baseline models on five datasets that resemble a variety of sequence prediction tasks. The experiments revealed that the proposed methods offer accuracy gain over baseline models in a range of settings.


title: HESS J0632+057: hydrodynamics and nonthermal emission

abstract: HESS J0632+057 is an eccentric gamma-ray Be binary that produces non-thermal radio, X-rays, GeV, and very high-energy gamma rays. The non-thermal emission of HESS J0632+057 is modulated with the orbital period, with a dominant maximum before apastron passage. The nature of the compact object in HESS J0632+057 is not known, although it has been proposed to be a young pulsar as in PSR B1259-63, the only gamma-ray emitting high-mass binary known to host a non-accreting pulsar. In this Letter, we present hydrodynamical simulations of HESS J0632+057 in the context of a pulsar and a stellar wind interacting in an eccentric binary, and propose a scenario for the non-thermal phenomenology of the source. In this scenario, the non-thermal activity before and around apastron is linked to the accumulation of non-thermal particles in the vicinity of the binary, and the sudden drop of the emission before apastron is produced by the disruption of the two-wind interaction structure, allowing these particles to efficiently escape. In addition to providing a framework to explain the non-thermal phenomenology of the source, this scenario predicts extended, moving X-ray emitting structures similar to those observed in PSR B1259-63.


title: On $A_p$ weights and the Landau equation

abstract: In this manuscript we investigate the regularization of solutions for the spatially homogeneous Landau equation. For moderately soft potentials, it is shown that weak solutions become smooth instantaneously and stay so over all times, and the estimates depend only on the initial mass, energy, and entropy. For very soft potentials we obtain a conditional regularity result, hinging on what may be described as a nonlinear Morrey space bound, assumed to hold uniformly over time. This bound always holds in the case of moderately soft potentials, and nearly holds for general potentials, including Coulomb. This latter phenomenon captures the intuition that for moderately soft potentials, the dissipative term in the equation is of the same order as the quadratic term driving the growth (and potentially, singularities). In particular, for the Coulomb case, the conditional regularity result shows a rate of regularization much stronger than what is usually expected for regular parabolic equations. The main feature of our proofs is the analysis of the linearized Landau operator around an arbitrary and possibly irregular distribution. This linear operator is shown to be a degenerate elliptic Schr\"odinger operator whose coefficients are controlled by $A_p$-weights.


title: Boundary streaming by internal waves

abstract: Damped internal wave beams in stratified fluids have long been known to generate strong mean flows through a mechanism analogous to acoustic streaming. While the role of viscous boundary layers in acoustic streaming has thoroughly been addressed, it remains largely unexplored in the case of internal waves. Here we compute the mean flow generated close to an undulating wall that emits internal waves in a viscous, linearly stratified two-dimensional Boussinesq fluid. Using a quasi-linear approach, we demonstrate that the form of the boundary conditions dramatically impacts the generated boundary streaming. In the no-slip scenario, the early time Reynolds stress divergence within the viscous boundary layer is much stronger than within the bulk while also driving flow in the opposite direction. Whatever the boundary condition, boundary streaming is however dominated by bulk streaming at a larger time. Using a WKB approach, we investigate the consequences of adding boundary streaming effects to an idealised model of wave-mean flow interactions known to reproduce the salient features of the quasi-biennial oscillation. The presence of wave boundary layers has a quantitative impact on the flow reversals.


title: Learning Robust Representations for Computer Vision

abstract: Unsupervised learning techniques in computer vision often require learning latent representations, such as low-dimensional linear and non-linear subspaces. Noise and outliers in the data can frustrate these approaches by obscuring the latent spaces. Our main goal is deeper understanding and new development of robust approaches for representation learning. We provide a new interpretation for existing robust approaches and present two specific contributions: a new robust PCA approach, which can separate foreground features from dynamic background, and a novel robust spectral clustering method, that can cluster facial images with high accuracy. Both contributions show superior performance to standard methods on real-world test sets.


title: Noise driven neuromorphic tuned amplifier

abstract: Living systems implement and execute an extraordinary plethora of computational tasks. The inherent degree of large scale coordination emerges as a global property, from the intricate sea of microscopic interactions. The brain, with its structural and functional architecture, represents an emblematic example of hierarchic self-organization: elementary units, the neurons, act much like instruments of an orchestra, which combine diverse timbres to create harmonious symphonies. Neurons come indeed in different types, varying in shapes, connections and electrical properties. They all team up to process external stimuli from a number of sources and integrate the information to yield, from neurons to mind, different cognitive faculties. Identifying the coarse grained modules that exert, from bottom to up, pivotal neuronal functions constitutes a goal of paramount importance. On the other hand, the brain and its unraveled secrets could inspire novel biomimetics technologies to adaptively handle complex problems. Here, we investigate the intertwined stochastic dynamics of two populations of excitatory and inhibitory units, arranged in a directed lattice. The endogenous noise seeds a coherent amplification across the chain generating giant oscillations with tunable frequencies, a process that the brain could exploit to enhance, and eventually encode, different signals. The system works as an out-equilibrium thermal device under stationary operating conditions, the associated entropy production rate being analytically quantified. The same scheme could be invoked to design a novel family of manmade detectors capable of reacting to spatially distributed low intensity alerts.

categories: cond-mat.dis-nn nlin.AO nlin.PS physics.bio-ph


title: On plane quartics with a Galois invariant Steiner hexad

abstract: We describe a construction of plane quartics with prescribed Galois operation on the 28 bitangents, in the particular case of a Galois invariant Steiner hexad. As an application, we solve the inverse Galois problem for degree two del Pezzo surfaces in the corresponding particular case.


title: A Component-oriented Framework for Autonomous Agents

abstract: The design of a complex system warrants a compositional methodology, i.e., composing simple components to obtain a larger system that exhibits their collective behavior in a meaningful way. We propose an automaton-based paradigm for compositional design of such systems where an action is accompanied by one or more preferences. At run-time, these preferences provide a natural fallback mechanism for the component, while at design-time they can be used to reason about the behavior of the component in an uncertain physical world. Using structures that tell us how to compose preferences and actions, we can compose formal representations of individual components or agents to obtain a representation of the composed system. We extend Linear Temporal Logic with two unary connectives that reflect the compositional structure of the actions, and show how it can be used to diagnose undesired behavior by tracing the falsification of a specification back to one or more culpable components.


title: Progress with Xenon Liquid Hole Multipliers

abstract: The bubble-assisted Liquid Hole Multiplier (LHM) is a recently-proposed concept for the combined detection of ionization electrons and primary scintillation photons in noble-liquid time projection chambers. The LHM comprises a perforated micro-pattern electrode (e.g. Thick Gas Electron Multiplier - THGEM, or Gas Electron Multiplier - GEM) immersed in the liquid, with a bubble of the noble gas supported underneath. Ionization electrons and scintillation-induced photoelectrons extracted from a cesium iodide photocathode drift through the electrode's holes and induce electroluminescence (EL) signals in the bubble; these are recorded by photon detectors located closely below the electrode. We present recent results in the development of LHMs, comparing the response of different electrodes to ionization and photon-induced electrons.


title: Point Transformations and the Relationships Among Anomalous Diffusion, Normal Diffusion and the Central Limit Theorem

abstract: We present new connections among anomalous diffusion (AD), normal diffusion (ND) and the Central Limit Theorem. This is done by defining a point transformation to a new position variable, which we postulate to be Cartesian, motivated by considerations from super-symmetric quantum mechanics. Canonically quantizing in the new position and momentum variables according to Dirac gives rise to generalized negative semi-definite and self-adjoint Laplacian operators. These lead to new generalized Fourier transformations and associated probability distributions, which are form invariant under the corresponding transform. The new Laplacians also lead us to generalized diffusion equations, which imply a connection to the CLT. We show that the derived diffusion equations capture all of the Fractal and Non-Fractal Diffusion equations of O'Shaughnessy and Procaccia. However, we also obtain new equations that cannot (so far as we are able to tell) be expressed as examples of the O'Shaughnessy and Procaccia equations. These equations also possess asymptotics that are related to a CLT but with bi-modal distributions as limits. The results show, in part, that experimentally measuring the diffusion scaling law can determine the point transformation (for monomial point transformations). We also show that AD in the original, physical position is actually ND when viewed in terms of displacements in an appropriately transformed position variable. Finally, we show that there is a new, anomalous diffusion possible for bi-modal probability distributions that also display attractor behavior which is the consequence of an underlying CLT.


title: Efficient Regret Minimization in Non-Convex Games

abstract: We consider regret minimization in repeated games with non-convex loss functions. Minimizing the standard notion of regret is computationally intractable. Thus, we define a natural notion of regret which permits efficient optimization and generalizes offline guarantees for convergence to an approximate local optimum. We give gradient-based methods that achieve optimal regret, which in turn guarantee convergence to equilibrium in this framework.


title: Capturing the Connections: Unboxing Internet of Things Devices

abstract: Based upon a study of how to capture data from Internet of Things (IoT) devices, this paper explores the challenges for data centric design ethnography. Often purchased to perform specific tasks, IoT devices exist in a complex ecosystem. This paper describes a study that used a variety of methods to capture the interactions an IoT device engaged in when it was first setup. The complexity of the study that is explored through the annotated documentation across video and router activity, presents the ethnographic challenges that designers face in an age of connected things.


title: Bayesian Sparsification of Recurrent Neural Networks

abstract: Recurrent neural networks show state-of-the-art results in many text analysis tasks but often require a lot of memory to store their weights. Recently proposed Sparse Variational Dropout eliminates the majority of the weights in a feed-forward neural network without significant loss of quality. We apply this technique to sparsify recurrent neural networks. To account for recurrent specifics we also rely on Binary Variational Dropout for RNN. We report 99.5% sparsity level on sentiment analysis task without a quality drop and up to 87% sparsity level on language modeling task with slight loss of accuracy.


title: Bayesian Dyadic Trees and Histograms for Regression

abstract: Many machine learning tools for regression are based on recursive partitioning of the covariate space into smaller regions, where the regression function can be estimated locally. Among these, regression trees and their ensembles have demonstrated impressive empirical performance. In this work, we shed light on the machinery behind Bayesian variants of these methods. In particular, we study Bayesian regression histograms, such as Bayesian dyadic trees, in the simple regression case with just one predictor. We focus on the reconstruction of regression surfaces that are piecewise constant, where the number of jumps is unknown. We show that with suitably designed priors, posterior distributions concentrate around the true step regression function at a near-minimax rate. These results do not require the knowledge of the true number of steps, nor the width of the true partitioning cells. Thus, Bayesian dyadic regression trees are fully adaptive and can recover the true piecewise regression function nearly as well as if we knew the exact number and location of jumps. Our results constitute the first step towards understanding why Bayesian trees and their ensembles have worked so well in practice. As an aside, we discuss prior distributions on balanced interval partitions and how they relate to an old problem in geometric probability. Namely, we relate the probability of covering the circumference of a circle with random arcs whose endpoints are confined to a grid, a new variant of the original problem.


title: Towards the Success Rate of One: Real-time Unconstrained Salient Object Detection

abstract: In this work, we propose an efficient and effective approach for unconstrained salient object detection in images using deep convolutional neural networks. Instead of generating thousands of candidate bounding boxes and refining them, our network directly learns to generate the saliency map containing the exact number of salient objects. During training, we convert the ground-truth rectangular boxes to Gaussian distributions that better capture the ROI regarding individual salient objects. During inference, the network predicts Gaussian distributions centered at salient objects with an appropriate covariance, from which bounding boxes are easily inferred. Notably, our network performs saliency map prediction without pixel-level annotations, salient object detection without object proposals, and salient object subitizing simultaneously, all in a single pass within a unified framework. Extensive experiments show that our approach outperforms existing methods on various datasets by a large margin, and achieves more than 100 fps with VGG16 network on a single GPU during inference.


title: A topological realization of the congruence subgroup Kernel A

abstract: A number of years ago, Kumar Murty pointed out to me that the computation of the fundamental group of a Hilbert modular surface ([7],IV,${\S}$6), and the computation of the congruence subgroup kernel of SL(2) ([6]) were surprisingly similar. We puzzled over this, in particular over the role of elementary matrices in both computations. We formulated a very general result on the fundamental group of a Satake compactification of a locally symmetric space. This lead to our joint paper [1] with Lizhen Ji and Les Saper on these fundamental groups. Although the results in it were intriguingly similar to the corresponding calculations of the congruence subgroup kernel of the underlying algebraic group in [5], we were not able to demonstrate a direct connection (cf. [1], ${\S}$7). The purpose of this note is to explain such a connection. A covering space is constructed from inverse limits of reductive Borel-Serre compactifications. The congruence subgroup kernel then appears as the group of deck transformations of this covering. The key to this is the computation of the fundamental group in [1].


title: Relating centrality to impact parameter in nucleus-nucleus collisions

abstract: In ultrarelativistic heavy-ion experiments, one estimates the centrality of a collision by using a single observable, say $n$, typically given by the transverse energy or the number of tracks observed in a dedicated detector. The correlation between $n$ and the impact parameter, $b$, of the collision is then inferred by fitting a specific model of the collision dynamics, such as the Glauber model, to experimental data. The goal of this paper is to assess precisely which information about $b$ can be extracted from data without any specific model of the collision. Under the sole assumption that the probability distribution of $n$ for a fixed $b$ is Gaussian, we show that the probability distribution of the impact parameter in a narrow centrality bin can be accurately reconstructed up to $5\%$ centrality. We apply our methodology to data from the Relativistic Heavy Ion Collider and the Large Hadron Collider. We propose a simple measure of the precision of the centrality determination, which can be used to compare different experiments.


title: On the zeros of the spectrogram of white noise

abstract: In a recent paper, Flandrin [2015] has proposed filtering based on the zeros of a spectrogram, using the short-time Fourier transform and a Gaussian window. His results are based on empirical observations on the distribution of the zeros of the spectrogram of white Gaussian noise. These zeros tend to be uniformly spread over the time-frequency plane, and not to clutter. Our contributions are threefold: we rigorously define the zeros of the spectrogram of continuous white Gaussian noise, we explicitly characterize their statistical distribution, and we investigate the computational and statistical underpinnings of the practical implementation of signal detection based on the statistics of spectrogram zeros. In particular, we stress that the zeros of spectrograms of white Gaussian noise correspond to zeros of Gaussian analytic functions, a topic of recent independent mathematical interest [Hough et al., 2009].


title: Hybrid Beamforming with Selection for Multi-user Massive MIMO Systems

abstract: This work studies a variant of hybrid beamforming, namely, hybrid beamforming with selection (HBwS), as an attractive solution to reduce the hardware cost of multi-user Massive Multiple-Input-Multiple-Output systems, while retaining good performance. Unlike conventional hybrid beamforming, in a transceiver with HBwS, the antenna array is fed by an analog beamforming matrix with $\bar{L}$ input ports, where $\bar{L}$ is larger than the number of up/down-conversion chains $\bar{K}$. A bank of switches connects the instantaneously best $\bar{K}$ out of the $\bar{L}$ input ports to the up/down-conversion chains. The analog beamformer is designed based on average channel statistics and therefore needs to be updated only infrequently, while the switches operate based on instantaneous channel knowledge. HBwS allows use of simpler hardware in the beamformer that only need to adjust to the statistics, while also enabling the effective analog beams to adapt to the instantaneous channel variations via switching. This provides better user separability, beamforming gain, and/or simpler hardware than some conventional hybrid schemes. In this work, a novel design for the analog beamformer is derived and approaches to reduce the hardware and computational cost of a multi-user HBwS system are explored. In addition, we study how $\bar{L}$, the switch bank architecture, the number of users and the channel estimation overhead impact system performance.


title: Compressed channeled spectropolarimetry

abstract: Channeled spectropolarimetry measures the spectrally resolved Stokes parameters. A key aspect of this technique is to accurately reconstruct the Stokes parameters from a modulated measurement of the channeled spectropolarimeter. The state-of-the-art reconstruction algorithm uses the Fourier transform to extract the Stokes parameters from channels in the Fourier domain. While this approach is straightforward, it can be sensitive to noise and channel cross-talk, and it imposes bandwidth limitations that cut off high frequency details. To overcome these drawbacks, we present a reconstruction method called compressed channeled spectropolarimetry. In our proposed framework, reconstruction in channeled spectropolarimetry is an underdetermined problem, where we take N measurements and solve for 3N unknown Stokes parameters. We formulate an optimization problem by creating a mathematical model of the channeled spectropolarimeter with inspiration from compressed sensing. We show that our approach offers greater noise robustness and reconstruction accuracy compared with the Fourier transform technique in simulations and experimental measurements. By demonstrating more accurate reconstructions, we push performance to the native resolution of the sensor, allowing more information to be recovered from a single measurement of a channeled spectropolarimeter.


title: Dynamic Variable Selection with Spike-and-Slab Process Priors

abstract: We address the problem of dynamic variable selection in time series regression with unknown residual variances, where the set of active predictors is allowed to evolve over time. To capture time-varying variable selection uncertainty, we introduce new dynamic shrinkage priors for the time series of regression coefficients. These priors are characterized by two main ingredients: smooth parameter evolutions and intermittent zeroes for modeling predictive breaks. More formally, our proposed Dynamic Spike-and-Slab (DSS) priors are constructed as mixtures of two processes: a spike process for the irrelevant coefficients and a slab autoregressive process for the active coefficients. The mixing weights are themselves time-varying and depend on lagged values of the series. Our DSS priors are probabilistically coherent in the sense that their stationary distribution is fully known and characterized by spike-and-slab marginals. For posterior sampling over dynamic regression coefficients, model selection indicators as well as unknown dynamic residual variances, we propose a Dynamic SSVS algorithm based on forward-filtering and backward-sampling. To scale our method to large data sets, we develop a Dynamic EMVS algorithm for MAP smoothing. We demonstrate, through simulation and a topical macroeconomic dataset, that DSS priors are very effective at separating active and noisy coefficients. Our fast implementation significantly extends the reach of spike-and-slab methods to large time series data.


title: A priori estimates for the free-boundary Euler equations with surface tension in three dimensions

abstract: We derive a priori estimates for the incompressible free-boundary Euler equations with surface tension in three spatial dimensions. Working in Lagrangian coordinates, we provide a priori estimates for the local existence when the initial velocity, which is rotational, belongs to $H^3$ and the trace of initial velocity on the free boundary to $H^{3.5}$, thus lowering the requirement on the regularity of initial data in the Lagrangian setting. Our methods are direct and involve three key elements: estimates for the pressure, the boundary regularity provided by the mean curvature, and the Cauchy invariance.


title: Investigating quantum wireless multihop teleportation under decoherence

abstract: This research work scrutinizes quantum routing protocol with multihop teleportation for wireless mesh backbone networks, in amplitude and phase damping channels. After analyzing the quantum multihop protocol, we select a four-qubit cluster state as the quantum channel for the protocol. The quantum channel linking the intermediate nodes has been established via entanglement swapping based on four-qubit cluster state. Also, we established the classical and the quantum route in a distributed manner. We show that from the source node to the destination node, quantum information can be teleported hop-by-hop through an amplitude damping channel. We show that the quantum teleportation could be successful if the sender node performs Bell state measurements (BSM), and the receiver introduces auxiliary particles, applies positive operative value measure and then utilizes corresponding unitary transformation to recover the transmitted state. We scrutinize the success probability of transferring the quantum state through a noisy channel. We found that optimum probability would be attained if decoherence rate of amplitude damping channel ($\xi_a$) is zero or the number of hops ($N$) is above $75$. Our numerical results evince susceptibility of success probability to $\xi_a$ and $N$. It has been shown that as the decoherence increases, the fidelity exponentially decays until it vanishes. This decay is as a consequence of information loss from the system to the surrounding. However, the fidelity can be enhanced by considering fewer hops.


title: Learning Algorithms for Active Learning

abstract: We introduce a model that learns active learning algorithms via metalearning. For a distribution of related tasks, our model jointly learns: a data representation, an item selection heuristic, and a method for constructing prediction functions from labeled training sets. Our model uses the item selection heuristic to gather labeled training sets from which to construct prediction functions. Using the Omniglot and MovieLens datasets, we test our model in synthetic and practical settings.


title: Formal Embeddings between $\mathcal{BSD}$-Models

abstract: It is studied the Classification Problem for Formal (Holomorphic) Embeddings between Shilov Boundaries of Bounded Symmetric Domains of First Type.


title: A Next Generation Low Band Observatory: A Community Study Exploring Low Frequency Options for ngVLA

abstract: We present a community study exploring the low frequency (5 - 800 MHz) options and opportunities for the ngVLA project and its infrastructure. We describe a Next Generation LOw Band Observatory (ngLOBO) that will provide access to the low frequency sky in a commensal fashion, operating independently from the ngVLA, but leveraging common infrastructure. This approach provides continuous coverage through an aperture array (called ngLOBO-Low) below 150 MHz and by accessing the primary focus of the ngVLA antennas (called ngLOBO-High) above 150 MHz. ngLOBO preconditions include a) non-interference and b) low relative cost (<5%) with respect to ngVLA. ngLOBO has three primary scientific missions: (1) Radio Large Synoptic Survey Telescope (Radio-LSST): one naturally wide beam, commensal with ngVLA, will conduct a continuous synoptic survey of large swaths of the sky for both slow and fast transients; (2) This same commensal beam will provide complementary low frequency images of all ngVLA targets and their environment {\it when such data enhances their value}. (3) Independent beams from the ngLOBO-Low aperture array will conduct research in astrophysics, Earth science and space weather applications, engaging new communities and attracting independent resources. If ngVLA operates down to 2 GHz or lower, ngLOBO data will enhance ngVLA calibration and dynamic scheduling. Finally, non-variable field sources outside the ngVLA field of view can be harvested for serendipitous science, e.g. population studies for thermal and non-thermal continuum sources. The ngVLA will be a superb, high frequency instrument; ngLOBO will allow it to participate in the worldwide renaissance in low frequency science as well.


title: Discrete probabilistic and algebraic dynamics: a stochastic commutative Gelfand-Naimark Theorem

abstract: We introduce a category of stochastic maps (certain Markov kernels) on compact Hausdorff spaces, construct a stochastic analogue of the Gelfand spectrum functor, and prove a stochastic version of the commutative Gelfand-Naimark Theorem. This relates concepts from algebra and operator theory to concepts from topology and probability theory. For completeness, we review stochastic matrices, their relationship to positive maps on commutative $C^*$-algebras, and the Gelfand-Naimark Theorem. No knowledge of probability theory nor $C^*$-algebras is assumed and several examples are drawn from physics.


title: Conditional Expectation Bounds with Applications in Cryptography

abstract: We derive two conditional expectation bounds, which we use to simplify cryptographic security proofs. The first bound relates the expectation of a bounded random variable and the average of its conditional expectations with respect to a set of i.i.d. random objects. It shows, under certain conditions, that the conditional expectation average has a small tail probability when the expectation of the random variable is sufficiently large. It is used to simplify the proof that the existence of weakly one-way functions implies the existence of strongly one-way functions. The second bound relaxes the independence requirement on the random objects to give a result that has applications to expander graph constructions in cryptography. It is used to simplify the proof that there is a security preserving reduction from weakly one-way functions to strongly one-way functions. To satisfy the hypothesis for this bound, we prove a hitting property for directed graphs that are expander-permutation hybrids.


title: Superadiabatic driving of a three-level quantum system

abstract: We study superadiabatic quantum control of a three-level quantum system whose energy spectrum exhibits multiple avoided crossings. In particular, we investigate the possibility of treating the full control task in terms of independent two-level Landau-Zener problems. We first show that the time profiles of the elements of the full control Hamiltonian are characterized by peaks centered around the crossing times. These peaks decay algebraically for large times. In principle, such a power-law scaling invalidates the hypothesis of perfect separability. Nonetheless, we address the problem from a pragmatic point of view by studying the fidelity obtained through separate control as a function of the intercrossing separation. This procedure may be a good approach to achieve approximate adiabatic driving of a specific instantaneous eigenstate in realistic implementations.


title: On facial unique-maximum (edge-)coloring

abstract: A facial unique-maximum coloring of a plane graph is a vertex coloring where on each face $\alpha$ the maximal color appears exactly once on the vertices of $\alpha$. If the coloring is required to be proper, then the upper bound for the minimal number of colors required for such a coloring is set to $5$. Fabrici and G\"oring [Fabrici and Goring 2016] even conjectured that $4$ colors always suffice. Confirming the conjecture would hence give a considerable strengthening of the Four Color Theorem. In this paper, we prove that the conjecture holds for subcubic plane graphs, outerplane graphs and plane quadrangulations. Additionally, we consider the facial edge-coloring analogue of the aforementioned coloring and prove that every $2$-connected plane graph admits such a coloring with at most $4$ colors.


title: The Projective Planarity Question for Matroids of $3$-Nets and Biased Graphs

abstract: A biased graph is a graph with a class of selected circles ("cycles", "circuits"), called "balanced", such that no theta subgraph contains exactly two balanced circles. A biased graph has two natural matroids, the frame matroid and the lift matroid. A classical question in matroid theory is whether a matroid can be embedded in a projective geometry. There is no known general answer, but for matroids of biased graphs it is possible to give algebraic criteria. Zaslavsky has previously given such criteria for embeddability of biased-graphic matroids in Desarguesian projective spaces; in this paper we establish criteria for the remaining case, that is, embeddability in an arbitrary projective plane that is not necessarily Desarguesian. The criteria depend on the embeddability of a quasigroup associated to the graph into the additive or multiplicative loop of a ternary coordinate ring for the plane. A 3-node biased graph is equivalent to an abstract partial 3-net; thus, we have a new algebraic criterion for an abstract 3-net to be realized in a non-Desarguesian projective plane. We work in terms of a special kind of 3-node biased graph called a biased expansion of a triangle. Our results apply to all finite 3-node biased graphs because, as we prove, every such biased graph is a subgraph of a finite biased expansion of a triangle. A biased expansion of a triangle, in turn, is equivalent to an isostrophe class of quasigroups, which is equivalent to a $3$-net. Much is not known about embedding a quasigroup into a ternary ring, so we do not say our criteria are definitive. For instance, it is not even known whether there is a finite quasigroup that cannot be embedded in any finite ternary ring. If there is, then there is a finite rank-3 matroid (of the corresponding biased expansion) that cannot be embedded in any finite projective plane---a presently unsolved problem.


title: Structure of transition metal clusters: A force-biased Monte Carlo approach

abstract: We present a force-biased Monte Carlo (FMC) method for structural modeling of transition metal clusters of Fe, Ni, and Cu with 5 to 60 atoms. By employing the Finnis-Sinclair potential for Fe and the Sutton-Chen potential for Ni and Cu, the total energy of the clusters is minimized using a method that utilizes atomic forces in Monte Carlo simulations. The structural configurations of the clusters obtained from this biased Monte Carlo approach are analyzed and compared with the same from the Cambridge Cluster Database (CCD). The results show that the total-energy of the FMC clusters is very close to the corresponding value of the CCD clusters as listed in the Cambridge Cluster Database. A comparison of the FMC and CCD clusters is presented by computing the pair-correlation function, the bond-angle distribution, and the distribution of atomic-coordination numbers in the first-coordination shell, which provide information about the two-body and three-body correlation functions, the local atomic structure, and the bonding environment of the atoms in the clusters.


title: Geometry of logarithmic forms and deformations of complex structures

abstract: We present a new method to solve certain $\bar{\partial}$-equations for logarithmic differential forms by using harmonic integral theory for currents on Kahler manifolds. The result can be considered as a $\bar{\partial}$-lemma for logarithmic forms. As applications, we generalize the result of Deligne about closedness of logarithmic forms, give geometric and simpler proofs of Deligne's degeneracy theorem for the logarithmic Hodge to de Rham spectral sequences at $E_1$-level, as well as certain injectivity theorem on compact Kahler manifolds. Furthermore, for a family of logarithmic deformations of complex structures on Kahler manifolds, we construct the extension for any logarithmic $(n,q)$-form on the central fiber and thus deduce the local stability of log Calabi-Yau structure by extending an iteration method to the logarithmic forms. Finally we prove the unobstructedness of the deformations of a log Calabi-Yau pair and a pair on a Calabi-Yau manifold by differential geometric method.


title: The Code2Text Challenge: Text Generation in Source Code Libraries

abstract: We propose a new shared task for tactical data-to-text generation in the domain of source code libraries. Specifically, we focus on text generation of function descriptions from example software projects. Data is drawn from existing resources used for studying the related problem of semantic parser induction (Richardson and Kuhn, 2017b; Richardson and Kuhn, 2017a), and spans a wide variety of both natural languages and programming languages. In this paper, we describe these existing resources, which will serve as training and development data for the task, and discuss plans for building new independent test sets.


title: Mixture Data-Dependent Priors

abstract: We propose a two-component mixture of a noninformative (diffuse) and an informative prior distribution, weighted through the data in such a way to prefer the first component if a prior-data conflict arises. The data-driven approach for computing the mixture weights makes this class data-dependent. Although rarely used with any theoretical motivation, data-dependent priors are often used for different reasons, and their use has been a lot debated over the last decades. However, our approach is justified in terms of Bayesian inference as an approximation of a hierarchical model and as a conditioning on a data statistic. This class of priors turns out to provide less information than an informative prior, perhaps it represents a suitable option for not dominating the inference in presence of small samples. First evidences from simulation studies show that this class could also be a good proposal for reducing mean squared errors.


title: Analysis of the Polya-Gamma block Gibbs sampler for Bayesian logistic linear mixed models

abstract: In this article, we construct a two-block Gibbs sampler using Polson et al. (2013) data augmentation technique with Polya-Gamma latent variables for Bayesian logistic linear mixed models under proper priors. Furthermore, we prove the uniform ergodicity of this Gibbs sampler, which guarantees the existence of the central limit theorems for MCMC based estimators.


title: MOSCAB: A geyser-concept bubble chamber to be used in a dark matter search

abstract: The MOSCAB experiment (Materia OSCura A Bolle) uses the "geyser technique", a variant of the superheated liquid technique of extreme simplicity. Operating principles of the new dark matter detector and technical solutions of the device are reported in detail. First results obtained in a series of test runs taken in laboratory demonstrate that we have successfully built and tested a geyser-concept bubble chamber that can be used in particle physics, especially in dark matter searches, and that we are ready to move underground for extensive data taking.


title: Advantages and Limitations of using Successor Features for Transfer in Reinforcement Learning

abstract: One question central to Reinforcement Learning is how to learn a feature representation that supports algorithm scaling and re-use of learned information from different tasks. Successor Features approach this problem by learning a feature representation that satisfies a temporal constraint. We present an implementation of an approach that decouples the feature representation from the reward function, making it suitable for transferring knowledge between domains. We then assess the advantages and limitations of using Successor Features for transfer.


title: Quantum triple point and quantum critical endpoints in metallic magnets

abstract: In low-temperature metallic magnets, ferromagnetic (FM) and antiferromagnetic (AFM) orders can exist in a single system in different parts of the phase diagram as a function of some control parameter. These phases can be adjacent, or exist concurrently, resulting in a phase transition between a FM phase and an AFM one, or between a phase of concurrent FM and AFM order and either of the pure phases. We show that universal quantum fluctuations qualitatively alter the known phase diagrams for classical magnets: They shrink the region of concurrent FM and AFM order, change various transitions from second to first order, and, in the presence of a magnetic field, lead to either a quantum triple point where the FM, AFM and paramagnetic phases coexist, or to a quantum critical end point.


title: High-power, high repetition rate, tunable, ultrafast vortex beam in the near-infrared

abstract: We report on experimental demonstration of high power, ultrafast, high repetition rate vortex beam source tunable in the near-IR wavelength range. Based on single-pass optical parametric generation of Yb-fiber laser of vortex order lp=1 in a 50 mm long MgO doped periodically poled LiNbO3 crystal, the source produces signal beam in vortex profile of order ls=1 across 1433-1553 nm. Additionally, the source produces broadband idler radiation tunable across 3379-4132 nm in the Gaussian beam profile. We observed that the vortex profile of the pump beam is always transferred to the signal beam due to the highest overlapping integral among the interacting beams and the idler maintains a Gaussian spatial profile owing to conservation of orbital angular momentum in optical parametric processes. For pump power of 4.72 W, the signal and idler beams have maximum power of 1.7 W at 1509 nm and 0.48 W at 3625 nm respectively. The signal vortex beam has output pulses of width 637 fs at a repetition rate of 78 MHz. The signal (idler) has spectral width of 4.3 nm (129.5 nm) and passive peak-to-peak power fluctuation better than 3% (1.1%) over 30 minutes, respectively.


title: Representations on Partially Holomorphic Cohomology Spaces, Revisited

abstract: This is a semi--expository update and rewrite of my 1974 AMS AMS Memoir describing Plancherel formulae and partial Dolbeault cohomology realizations for standard tempered representations for general real reductive Lie groups. Even after so many years, much of that Memoir is up to date, but of course there have been a number of refinements, advances and new developments, most of which have applied to smaller classes of real reductive Lie groups. Here we rewrite that AMS Memoir in in view of these advances and indicate the ties with some of the more recent (or at least less classical) approaches to geometric realization of unitary representations.


title: Material Editing Using a Physically Based Rendering Network

abstract: The ability to edit materials of objects in images is desirable by many content creators. However, this is an extremely challenging task as it requires to disentangle intrinsic physical properties of an image. We propose an end-to-end network architecture that replicates the forward image formation process to accomplish this task. Specifically, given a single image, the network first predicts intrinsic properties, i.e. shape, illumination, and material, which are then provided to a rendering layer. This layer performs in-network image synthesis, thereby enabling the network to understand the physics behind the image formation process. The proposed rendering layer is fully differentiable, supports both diffuse and specular materials, and thus can be applicable in a variety of problem settings. We demonstrate a rich set of visually plausible material editing examples and provide an extensive comparative study.


title: Learned in Translation: Contextualized Word Vectors

abstract: Computer vision has benefited from initializing multiple deep layers with weights pretrained on large supervised training sets like ImageNet. Natural language processing (NLP) typically sees initialization of only the lowest layer of deep models with pretrained word vectors. In this paper, we use a deep LSTM encoder from an attentional sequence-to-sequence model trained for machine translation (MT) to contextualize word vectors. We show that adding these context vectors (CoVe) improves performance over using only unsupervised word and character vectors on a wide variety of common NLP tasks: sentiment analysis (SST, IMDb), question classification (TREC), entailment (SNLI), and question answering (SQuAD). For fine-grained sentiment analysis and entailment, CoVe improves performance of our baseline models to the state of the art.


title: Experimental characterization of electron beam driven wakefield modes in a dielectric woodpile Cartesian symmetric structure

abstract: The use of photonic structures in the terahertz (THz) spectral region may enable the essential characteristics of confinement, modal control, and electric field shielding for very high gradient accelerators based on wakefields in dielectrics. We report here an experimental investigation of THz wakefield modes in a 3D photonic woodpile structure. Selective control in exciting or suppressing of wakefield modes with non-zero transverse wave vector is demonstrated by using drive beams of varying transverse ellipticity. Additionally, we show that the wakefield spectrum is insensitive to the offset position of the strongly elliptical beam. These results are consistent with analytic theory and 3D simulations, and illustrate a key advantage of wakefield systems with Cartesian symmetry, the suppression of transverse wakes by elliptical beams.


title: A Labelling Framework for Probabilistic Argumentation

abstract: The combination of argumentation and probability paves the way to new accounts of qualitative and quantitative uncertainty, thereby offering new theoretical and applicative opportunities. Due to a variety of interests, probabilistic argumentation is approached in the literature with different frameworks, pertaining to structured and abstract argumentation, and with respect to diverse types of uncertainty, in particular the uncertainty on the credibility of the premises, the uncertainty about which arguments to consider, and the uncertainty on the acceptance status of arguments or statements. Towards a general framework for probabilistic argumentation, we investigate a labelling-oriented framework encompassing a basic setting for rule-based argumentation and its (semi-) abstract account, along with diverse types of uncertainty. Our framework provides a systematic treatment of various kinds of uncertainty and of their relationships and allows us to back or question assertions from the literature.


title: Naturally occurring $^{32}$Si and low-background silicon dark matter detectors

abstract: The naturally occurring radioisotope $^{32}$Si represents a potentially limiting background in future dark matter direct-detection experiments. We investigate sources of $^{32}$Si and the vectors by which it comes to reside in silicon crystals used for fabrication of radiation detectors. We infer that the $^{32}$Si concentration in commercial single-crystal silicon is likely variable, dependent upon the specific geologic and hydrologic history of the source (or sources) of silicon "ore" and the details of the silicon-refinement process. The silicon production industry is large, highly segmented by refining step, and multifaceted in terms of final product type, from which we conclude that production of $^{32}$Si-mitigated crystals requires both targeted silicon material selection and a dedicated refinement-through-crystal-production process. We review options for source material selection, including quartz from an underground source and silicon isotopically reduced in $^{32}$Si. To quantitatively evaluate the $^{32}$Si content in silicon metal and precursor materials, we propose analytic methods employing chemical processing and radiometric measurements. Ultimately, it appears feasible to produce silicon detectors with low levels of $^{32}$Si, though significant assay method development is required to validate this claim and thereby enable a quality assurance program during an actual controlled silicon-detector production cycle.


title: A Continuous Relaxation of Beam Search for End-to-end Training of Neural Sequence Models

abstract: Beam search is a desirable choice of test-time decoding algorithm for neural sequence models because it potentially avoids search errors made by simpler greedy methods. However, typical cross entropy training procedures for these models do not directly consider the behaviour of the final decoding method. As a result, for cross-entropy trained models, beam decoding can sometimes yield reduced test performance when compared with greedy decoding. In order to train models that can more effectively make use of beam search, we propose a new training procedure that focuses on the final loss metric (e.g. Hamming loss) evaluated on the output of beam search. While well-defined, this "direct loss" objective is itself discontinuous and thus difficult to optimize. Hence, in our approach, we form a sub-differentiable surrogate objective by introducing a novel continuous approximation of the beam search decoding procedure. In experiments, we show that optimizing this new training objective yields substantially better results on two sequence tasks (Named Entity Recognition and CCG Supertagging) when compared with both cross entropy trained greedy decoding and cross entropy trained beam decoding baselines.


title: Retrofitting Distributional Embeddings to Knowledge Graphs with Functional Relations

abstract: Knowledge graphs are a versatile framework to encode richly structured data relationships, but it can be challenging to combine these graphs with unstructured data. Methods for retrofitting pre-trained entity representations to the structure of a knowledge graph typically assume that entities are embedded in a connected space and that relations imply similarity. However, useful knowledge graphs often contain diverse entities and relations (with potentially disjoint underlying corpora) which do not accord with these assumptions. To overcome these limitations, we present Functional Retrofitting, a framework that generalizes current retrofitting methods by explicitly modeling pairwise relations. Our framework can directly incorporate a variety of pairwise penalty functions previously developed for knowledge graph completion. Further, it allows users to encode, learn, and extract information about relation semantics. We present both linear and neural instantiations of the framework. Functional Retrofitting significantly outperforms existing retrofitting methods on complex knowledge graphs and loses no accuracy on simpler graphs (in which relations do imply similarity). Finally, we demonstrate the utility of the framework by predicting new drug--disease treatment pairs in a large, complex health knowledge graph.


title: Prediction and Generation of Binary Markov Processes: Can a Finite-State Fox Catch a Markov Mouse?

abstract: Understanding the generative mechanism of a natural system is a vital component of the scientific method. Here, we investigate one of the fundamental steps toward this goal by presenting the minimal generator of an arbitrary binary Markov process. This is a class of processes whose predictive model is well known. Surprisingly, the generative model requires three distinct topologies for different regions of parameter space. We show that a previously proposed generator for a particular set of binary Markov processes is, in fact, not minimal. Our results shed the first quantitative light on the relative (minimal) costs of prediction and generation. We find, for instance, that the difference between prediction and generation is maximized when the process is approximately independently, identically distributed.

categories: cond-mat.stat-mech cs.CC cs.IT math.IT nlin.CD


title: Thermodynamic Framework for Compact q-Gaussian Distributions

abstract: Recent works have associated systems of particles, characterized by short-range repulsive interactions and evolving under overdamped motion, to a nonlinear Fokker-Planck equation within the class of nonextensive statistical mechanics, with a nonlinear diffusion contribution whose exponent is given by $\nu=2-q$. The particular case $\nu=2$ applies to interacting vortices in type-II superconductors, whereas $\nu>2$ covers systems of particles characterized by short-range power-law interactions, where correlations among particles are taken into account. In the former case, several studies presented a consistent thermodynamic framework based on the definition of an effective temperature $\theta$ (presenting experimental values much higher than typical room temperatures $T$, so that thermal noise could be neglected), conjugated to a generalized entropy $s_{\nu}$ (with $\nu=2$). Herein, the whole thermodynamic scheme is revisited and extended to systems of particles interacting repulsively, through short-ranged potentials, described by an entropy $s_{\nu}$, with $\nu>1$, covering the $\nu=2$ (vortices in type-II superconductors) and $\nu>2$ (short-range power-law interactions) physical examples. The main results achieved are: (a) The definition of an effective temperature $\theta$ conjugated to the entropy $s_{\nu}$; (b) The construction of a Carnot cycle, whose efficiency is shown to be $\eta=1-(\theta_2/\theta_1)$, where $\theta_1$ and $\theta_2$ are the effective temperatures associated with two isothermal transformations, with $\theta_1>\theta_2$; (c) Thermodynamic potentials, Maxwell relations, and response functions. The present thermodynamic framework, for a system of interacting particles under the above-mentioned conditions, and associated to an entropy $s_{\nu}$, with $\nu>1$, certainly enlarges the possibility of experimental verifications.


title: On Arithmetic Series involving the fractional part function

abstract: We present new relationships between the work of H. Davenport and A. I. Popov. A new general formula involving the Von Mangoldt function is presented, as well as a criteria for the Riemann Hypothesis.


title: Solving an abstract nonlinear eigenvalue problem by the inverse iteration method

abstract: Let $\left( X,\left\Vert \cdot\right\Vert_{X}\right) $ and $\left( Y,\left\Vert \cdot\right\Vert_{Y}\right) $ be Banach spaces over $\mathbb{R},$ with $X$ uniformly convex and compactly embedded into $Y.$ The inverse iteration method is applied to solve the abstract eigenvalue problem $A(w)=\lambda\left\Vert w\right\Vert_{Y}^{p-q}B(w),$ where the maps $A:X\rightarrow X^{\star}$ and $B:Y\rightarrow Y^{\star}$ are homogeneous of degrees $p-1$ and $q-1,$ respectively.


title: Compiling Deep Learning Models for Custom Hardware Accelerators

abstract: Convolutional neural networks (CNNs) are the core of most state-of-the-art deep learning algorithms specialized for object detection and classification. CNNs are both computationally complex and embarrassingly parallel. Two properties that leave room for potential software and hardware optimizations for embedded systems. Given a programmable hardware accelerator with a CNN oriented custom instructions set, the compiler's task is to exploit the hardware's full potential, while abiding with the hardware constraints and maintaining generality to run different CNN models with varying workload properties. Snowflake is an efficient and scalable hardware accelerator implemented on programmable logic devices. It implements a control pipeline for a custom instruction set. The goal of this paper is to present Snowflake's compiler that generates machine level instructions from Torch7 model description files. The main software design points explored in this work are: model structure parsing, CNN workload breakdown, loop rearrangement for memory bandwidth optimizations and memory access balancing. The performance achieved by compiler generated instructions matches against hand optimized code for convolution layers. Generated instructions also efficiently execute AlexNet and ResNet18 inference on Snowflake. Snowflake with $256$ processing units was synthesized on Xilinx's Zynq XC7Z045 FPGA. At $250$ MHz, AlexNet achieved in $93.6$ frames/s and $1.2$ GB/s of off-chip memory bandwidth, and $21.4$ frames/s and $2.2$ GB/s for ResNet18. Total on-chip power is $5$ W.


title: Anomaly Detection Using Optimally-Placed Micro-PMU Sensors in Distribution Grids

abstract: As the distribution grid moves toward a tightly-monitored network, it is important to automate the analysis of the enormous amount of data produced by the sensors to increase the operators situational awareness about the system. In this paper, focusing on Micro-Phasor Measurement Unit ($\mu$PMU) data, we propose a hierarchical architecture for monitoring the grid and establish a set of analytics and sensor fusion primitives for the detection of abnormal behavior in the control perimeter. Due to the key role of the $\mu$PMU devices in our architecture, a source-constrained optimal $\mu$PMU placement is also described that finds the best location of the devices with respect to our rules. The effectiveness of the proposed methods are tested through the synthetic and real $\mu$PMU data.


title: Classifying Families of Character Degree Graphs of Solvable Groups

abstract: We investigate prime character degree graphs of solvable groups. In particular, we consider a family of graphs $\Gamma_{k,t}$ constructed by adjoining edges between two complete graphs in a one-to-one fashion. In this paper we determine completely which graphs $\Gamma_{k,t}$ occur as the prime character degree graph of a solvable group.


title: Multiple Stakeholders in Music Recommender Systems

abstract: Music recommendation services collectively spin billions of songs for millions of listeners on a daily basis. Users can typically listen to a variety of songs tailored to their personal tastes and preferences. Music is not the only type of content encountered in these services, however. Advertisements are generally interspersed throughout the music stream to generate revenue for the business. Additional content may include artist messaging, ticketing, sports, news and weather. In this paper, we discuss issues that arise when multiple content providers are stakeholders in the recommendation process. These stakeholders each have their own objectives and must work in concert to sustain a healthy music recommendation service.


title: Computing the Margin of Victory in Preferential Parliamentary Elections

abstract: We show how to use automated computation of election margins to assess the number of votes that would need to change in order to alter a parliamentary outcome for single-member preferential electorates. In the context of increasing automation of Australian electoral processes, and accusations of deliberate interference in elections in Europe and the USA, this work forms the basis of a rigorous statistical audit of the parliamentary election outcome. Our example is the New South Wales Legislative Council election of 2015, but the same process could be used for any similar parliament for which data was available, such as the Australian House of Representatives given the proposed automatic scanning of ballots.


title: Application of Support Vector Machine Modeling and Graph Theory Metrics for Disease Classification

abstract: Disease classification is a crucial element of biomedical research. Recent studies have demonstrated that machine learning techniques, such as Support Vector Machine (SVM) modeling, produce similar or improved predictive capabilities in comparison to the traditional method of Logistic Regression. In addition, it has been found that social network metrics can provide useful predictive information for disease modeling. In this study, we combine simulated social network metrics with SVM to predict diabetes in a sample of data from the Behavioral Risk Factor Surveillance System. In this dataset, Logistic Regression outperformed SVM with ROC index of 81.8 and 81.7 for models with and without graph metrics, respectively. SVM with a polynomial kernel had ROC index of 72.9 and 75.6 for models with and without graph metrics, respectively. Although this did not perform as well as Logistic Regression, the results are consistent with previous studies utilizing SVM to classify diabetes.


title: Quantum Projective Simulation with Hamiltonian Evolution: A study in reinforcement learning

abstract: Projective Simulation was introduced as a novel approach to Artificial Intelligence. It involves a deliberation procedure that consists of a random walk on a graph of clips and allows for the learning agent to project itself into the future before committing to an action. Here we study and analyze a quantum mechanical version in which the random walk is performed by two kinds of Hamiltonians. The first kind is implemented by naively embedding the classical model in a quantum model by turning the clips into qubits. The other allows for storing clips in superpositions of qubits allowing for a potentially purely quantum mechanical learning procedure in which the perception of the environment is purely quantum mechanical but the action is classical. We lastly introduce the concept of interacting projective agents for both the classical and quantum mechanical case.


title: Kosterlitz-Thouless transition and vortex-antivortex lattice melting in two-dimensional Fermi gases with $p$- or $d$-wave pairing

abstract: We present a theoretical study of the finite-temperature Kosterlitz-Thouless (KT) and vortex-antivortex lattice (VAL) melting transitions in two-dimensional Fermi gases with $p$- or $d$-wave pairing. For both pairings, when the interaction is tuned from weak to strong attractions, we observe a quantum phase transition from the Bardeen-Cooper-Schrieffer (BCS) superfluidity to the Bose-Einstein condensation (BEC) of difermions. The KT and VAL transition temperatures increase during this BCS-BEC transition and approach constant values in the deep BEC region. The BCS-BEC transition is characterized by the non-analyticities of the chemical potential, the superfluid order parameter, and the sound velocities as functions of the interaction strength at both zero and finite temperatures; however, the temperature effect tends to weaken the non-analyticities comparing to the zero temperature case. The effect of mismatched Fermi surfaces on the $d$-wave pairing is also studied.


title: A Note on Upper Bounds for Some Generalized Folkman Numbers

abstract: We present some new constructive upper bounds based on product graphs for generalized vertex Folkman numbers. They lead to new upper bounds for some special cases of generalized edge Folkman numbers, including $F_e(K_3,K_4-e; K_5) \leq 27$ and $F_e(K_4-e,K_4-e; K_5) \leq 51$. The latter bound follows from a construction of a $K_5$-free graph on 51 vertices, for which every coloring of its edges with two colors contains a monochromatic $K_4-e$.


title: Excitonic Resonance Effects and Davydov Splitting in Circularly Polarized Raman Spectra of Few-Layer WSe2

abstract: Few-layer tungsten diselenide (WSe2) is investigated using circularly polarized Raman spectroscopy with up to eight excitation energies. The main E2g1 and A1g modes near 250 cm-1 appear as a single peak in the Raman spectrum taken without consideration of polarization but are resolved by using circularly polarized Raman scattering. The resonance behaviors of the E2g1 and A1g modes are examined. Firstly, both the E2g1 and A1g modes are enhanced near resonances with the exciton states. The A1g mode exhibits Davydov splitting for trilayers or thicker near some of the exciton resonances. The low-frequency Raman spectra show shear and breathing modes involving rigid vibrations of the layers and also exhibit strong dependence on the excitation energy. An unidentified peak at ~19 cm-1 that does not depend on the number of layers appears near resonance with the B exciton state at 1.96 eV (632.8 nm). The strengths of the intra- and inter-layer interactions are estimated by comparing the mode frequencies and Davydov splitting with the linear chain model, and the contribution of the next-nearest-neighbor interaction to the inter-layer interaction turns out to be about 34% of the nearest-neighbor interaction. Fano resonance is observed for 1.58-eV excitation, and its origin is found to be the interplay between two-phonon scattering and indirect band transition.


title: Periodic fourth-order cubic NLS: Local well-posedness and Non-squeezing property

abstract: In this paper, we consider the cubic fourth-order nonlinear Schr\"odinger equation (4NLS) under the periodic boundary condition. We prove two results. One is the local well-posedness in $H^s$ with $-1/3 \le s < 0$ for the Cauchy problem of the Wick ordered 4NLS. The other one is the non-squeezing property for the flow map of 4NLS in the symplectic phase space $L^2(\mathbb{T})$. To prove the former we used the ideas introduced in [Takaoka and Tsutsumi 2004] and [Nakanish et al 2010], and to prove the latter we used the ideas in [Colliander et al 2005].


title: Imaging from the Inside Out: Inverse Scattering with Photoactivated Internal Sources

abstract: We propose a method to reconstruct the optical properties of a scattering medium with subwavelength resolution. The method is based on the solution to the inverse scattering problem with photoactivated internal sources. Numerical simulations of three-dimensional structures demonstrate that a resolution of approximately $\lambda/25$ is achievable.


title: Deep Generative Adversarial Neural Networks for Realistic Prostate Lesion MRI Synthesis

abstract: Generative Adversarial Neural Networks (GANs) are applied to the synthetic generation of prostate lesion MRI images. GANs have been applied to a variety of natural images, is shown show that the same techniques can be used in the medical domain to create realistic looking synthetic lesion images. 16mm x 16mm patches are extracted from 330 MRI scans from the SPIE ProstateX Challenge 2016 and used to train a Deep Convolutional Generative Adversarial Neural Network (DCGAN) utilizing cutting edge techniques. Synthetic outputs are compared to real images and the implicit latent representations induced by the GAN are explored. Training techniques and successful neural network architectures are explained in detail.


title: Predicting Session Length in Media Streaming

abstract: Session length is a very important aspect in determining a user's satisfaction with a media streaming service. Being able to predict how long a session will last can be of great use for various downstream tasks, such as recommendations and ad scheduling. Most of the related literature on user interaction duration has focused on dwell time for websites, usually in the context of approximating post-click satisfaction either in search results, or display ads. In this work we present the first analysis of session length in a mobile-focused online service, using a real world data-set from a major music streaming service. We use survival analysis techniques to show that the characteristics of the length distributions can differ significantly between users, and use gradient boosted trees with appropriate objectives to predict the length of a session using only information available at its beginning. Our evaluation on real world data illustrates that our proposed technique outperforms the considered baseline.


title: Reconfiguration of quantum states in $\mathcal PT$-symmetric quasi-one dimensional lattices

abstract: We demonstrate mesoscopic transport through quantum states in quasi-1D lattices maintaining the combination of parity and time-reversal symmetries by controlling energy gain and loss. We investigate the phase diagram of the non-Hermitian system where transitions take place between unbroken and broken $\mathcal{PT}$-symmetric phases via exceptional points. Quantum transport in the lattice is measured only in the unbroken phases in the energy band-but not in the broken phases. The broken phase allows for spontaneous symmetry-broken states where the cross-stitch lattice is separated into two identical single lattices corresponding to conditionally degenerate eigenstates. These degeneracies show a lift-up in the complex energy plane, caused by the non-Hermiticity with $\mathcal{PT}$-symmetry.


title: On Tensor Train Rank Minimization: Statistical Efficiency and Scalable Algorithm

abstract: Tensor train (TT) decomposition provides a space-efficient representation for higher-order tensors. Despite its advantage, we face two crucial limitations when we apply the TT decomposition to machine learning problems: the lack of statistical theory and of scalable algorithms. In this paper, we address the limitations. First, we introduce a convex relaxation of the TT decomposition problem and derive its error bound for the tensor completion task. Next, we develop an alternating optimization method with a randomization technique, in which the time complexity is as efficient as the space complexity is. In experiments, we numerically confirm the derived bounds and empirically demonstrate the performance of our method with a real higher-order tensor.


title: Grounding Language for Transfer in Deep Reinforcement Learning

abstract: In this paper, we explore the utilization of natural language to drive transfer for reinforcement learning (RL). Despite the wide-spread application of deep RL techniques, learning generalized policy representations that work across domains remains a challenging problem. We demonstrate that textual descriptions of environments provide a compact intermediate channel to facilitate effective policy transfer. Specifically, by learning to ground the meaning of text to the dynamics of the environment such as transitions and rewards, an autonomous agent can effectively bootstrap policy learning on a new domain given its description. We employ a model-based RL approach consisting of a differentiable planning module, a model-free component and a factorized state representation to effectively use entity descriptions. Our model outperforms prior work on both transfer and multi-task scenarios in a variety of different environments. For instance, we achieve up to 14% and 11.5% absolute improvement over previously existing models in terms of average and initial rewards, respectively.


title: Repulsive Fermi polarons with negative effective mass

abstract: Recent LENS experiment on a 3D Fermi gas has reported a negative effective mass ($m^*<0$) of Fermi polarons in the strongly repulsive regime. There naturally arise a question whether the negative $m^*$ is a precursor of the instability towards phase separation (or itinerant ferromagnetism). In this work, we make use of the exact solutions to study the ground state and excitation properties of repulsive Fermi polarons in 1D, which can also exhibit a negative $m^*$ in the super Tonks-Girardeau regime. By analyzing the total spin, quasi-momentum distribution and pair correlations, we conclude that the negative $m^*$ is irrelevant to the instability towards ferromagnetism or phase separation, but rather an intrinsic feature of collective excitations for fermions in the strongly repulsive regime. Surprisingly, for large and negative $m^*$, such excitation is accompanied with a spin density modulation when the majority fermions move closer to the impurity rather than being repelled far away, contrary to the picture of phase separation. These results suggest an alternative interpretation of negative $m^*$ as observed in recent LENS experiment.


title: $s$-wave Contacts of Quantum Gases in Quasi-one and Quasi-two Dimensions

abstract: In quasi-one- or quasi-two-dimensional traps with strong transverse confinements, quantum gases behave like strictly one- or two-dimensional systems at large length scales. However, at short distance, the two-body scattering intrinsically has three-dimensional characteristics such that an exact description of any universal thermodynamic relation requires three-dimensional contacts, no matter how strong the confinement is. A fundamental question arises as to whether one- or two-dimensional contacts, which were originally defined for strictly one or two dimensions, are capable of describing quantum gases in quasi-one- or quasi-two-dimensional traps. Here, we point out an exact relation between the three- and low-dimensional contacts in these highly anisotropic traps. Such relation allows us to directly connect physical quantities at different length scales, and to characterise the quasi-one- or quasi-two-dimensional traps using universal thermodynamic relations that were derived for strict one or two dimensions.


title: TNSPackage: A Fortran2003 library designed for tensor network state methods

abstract: Recently, the tensor network states (TNS) methods have proven to be very powerful tools to investigate the strongly correlated many-particle physics in one and two dimensions. The implementation of TNS methods depends heavily on the operations of tensors, including contraction, permutation, reshaping tensors, SVD and son on. Unfortunately, the most popular computer languages for scientific computation, such as Fortran and C/C++ do not have a standard library for such operations, and therefore make the coding of TNS very tedious. We develop a Fortran2003 package that includes all kinds of basic tensor operations designed for TNS. It is user-friendly and flexible for different forms of TNS, and therefore greatly simplifies the coding work for the TNS methods.


title: Causation does not explain contextuality

abstract: Realist interpretations of quantum mechanics presuppose the existence of elements of reality that are independent of the actions used to reveal them. Such a view is challenged by several no-go theorems that show quantum correlations cannot be explained by non-contextual ontological models, where physical properties are assumed to exist prior to and independently of the act of measurement. However, all such contextuality proofs assume a traditional notion of causal structure, where causal influence flows from past to future according to ordinary dynamical laws. This leaves open the question of whether the apparent contextuality of quantum mechanics is simply the signature of some exotic causal structure, where the future might affect the past or distant systems might get correlated due to non-local constraints. Here we show that quantum predictions require a deeper form of contextuality: even allowing for arbitrary causal structure, no model can explain quantum correlations from non-contextual ontological properties of the world, be they initial states, dynamical laws, or global constraints.


title: The differential geometry of perceptual similarity

abstract: Human similarity judgments are inconsistent with Euclidean, Hamming, Mahalanobis, and the majority of measures used in the extensive literatures on similarity and dissimilarity. From intrinsic properties of brain circuitry, we derive principles of perceptual metrics, showing their conformance to Riemannian geometry. As a demonstration of their utility, the perceptual metrics are shown to outperform JPEG compression. Unlike machine-learning approaches, the outperformance uses no statistics, and no learning. Beyond the incidental application to compression, the metrics offer broad explanatory accounts of empirical perceptual findings such as Tverskys triangle inequality violations, contradictory human judgments of identical stimuli such as speech sounds, and a broad range of other phenomena on percepts and concepts that may initially appear unrelated. The findings constitute a set of fundamental principles underlying perceptual similarity.


title: Citizen Astronomy in China: Present and Future

abstract: Citizen science refers to scientific research conducted or participated by non-professional scientists (such as hobbyists or members from the general public). Citizen astronomy is a classic example of citizen science. Citizen astronomers has benefited from technological advancements in the recent decades as they fill the scientific gaps left by professional astronomers, in the areas such as time domain observations, visual classification and data mining. Chinese citizen astronomers have made a visible contribution in the discoveries of new objects; however, comparing to their counterparts in the western world, they appear to be less interested in researches that do not involve making new discovery, such as visual classification, long-term monitoring of objects, and data mining. From a questionnaire survey that aimed to investigate the motivation of Chinese citizen astronomers, we find that this population is predominantly male ($92\%$) who mostly reside in economically developed provinces. A large fraction ($69\%$) of the respondents are students and young professionals younger than the age of 25, which differs significantly from the occupation and age distribution of typical Chinese Internet users as well as the user distribution of large international citizen science projects such as the Galaxy Zoo. This suggests that youth generation in China is more willing to participate citizen astronomy research than average generation. Additionally, we find that interests in astronomy, desire to learn new knowledges, have a fun experience and meet new friends in the community are all important driving factors for Chinese citizen astronomers to participate research. This also differs from their counterparts in western countries. (full abstract in PDF)


title: Digit Serial Methods with Applications to Division and Square Root (with mechanically checked correctness proofs)

abstract: We present a generic digit serial method (DSM) to compute the digits of a real number $V$ . Bounds on these digits, and on the errors in the associated estimates of $V$ formed from these digits, are derived. To illustrate our results, we derive such bounds for a parameterized family of high-radix algorithms for division and square root. These bounds enable a DSM designer to determine, for example, whether a given choice of parameters allows rapid formation and rounding of its approximation to $V$. All our claims are mechanically verified using the HOL-Light theorem prover, and are included in the appendix with commentary.


title: Chern-Ricci flows on noncompact complex manifolds

abstract: In this work, we obtain existence criteria for Chern-Ricci flows on noncompact manifolds. We generalize a result by Tossati-Wienkove on Chern-Ricci flows to noncompact manifolds and at the same time generalize a result for Kahler-Ricci flows by Lott-Zhang to Chern-Ricci flows. Using the existence results, we prove that any complete noncollapsed Kahler metric with nonnegative bisectional curvature on a noncompact complex manifold can be deformed to a complete Kahler metric with nonnegative and bounded bisectional curvature which will have maximal volume growth if the initial metric has maximal volume. Combining this result with the result of Chau-Tam, we give another proof that a complete noncompact Kahler manifold with nonnegative bisectional curvature (not necessarily bounded) and maximal volume growth is biholomorphic to the complex Euclidean space. This last result has already been proved by Gang Liu recently using other methods. This last result is partial confirmation of a uniformization conjecture of Yau.


title: Triply heavy $QQ\bar Q\bar q$ tetraquark states

abstract: Within the framework of QCD sum rules, we have investigated the tetraquark states with three heavy quarks. We systematically construct the interpolating currents for the possible $cc\bar{c}\bar{q}$, $cc\bar{b}\bar{q}$, $bc\bar{b}\bar{q}$, $bb\bar{b}\bar{q}$ tetraquark states with quantum numbers $J^{P}=0^{+}$ and $J^{P}=1^{+}$. Using these interpolating currents, we have calculated the two-point correlation functions and extracted the mass spectra for the above tetraquark states. We also discuss the decay patterns of these tetraquarks, and notice that the $cc\bar{c}\bar{q}$, $cc\bar{b}\bar{q}$, $bc\bar{b}\bar{q}$ may decay quickly with a narrow width due to their mass spectra. The $bb\bar{b}\bar{q}$ tetraquarks are expected to be very narrow resonances since their OZI-allowed decay modes are kinematically forbidden. These states may be searched for in the final states with a $B$ meson plus a light meson or photon.


title: Improved Algorithms for Scheduling Unsplittable Flows on Paths

abstract: In this paper, we investigate offline and online algorithms for rufpp, the problem of minimizing the number of rounds required to schedule a set of unsplittable flows of non-uniform sizes on a given path with non-uniform edge capacities. rufpp is NP-hard and constant-factor approximation algorithms are known under the no bottleneck assumption (NBA), which stipulates that maximum size of a flow is at most the minimum edge capacity. We study rufpp without the NBA, and present improved online and offline algorithms. We first study offline rufpp for a restricted class of instances called $\alpha$-small, where the size of each flow is at most $\alpha$ times the capacity of its bottleneck edge, and present an $O(\log(1/(1-\alpha)))$-approximation algorithm. Our main result is an online $O(\log\log c_{\max})$-competitive algorithm for rufpp for general instances, where $c_{\max}$ is the largest edge capacities, improving upon the previous best bound of $O(\log c_{\max})$ due to Epstein et al. Our result leads to an offline $O(\min(\log n, \log m, \log\log c_{\max}))$-approximation algorithm and an online $O(\min(\log m, \log\log c_{\max}))$-competitive algorithm for rufpp, where $n$ is the number of flows and $m$ is the number of edges.


title: On A Conjecture Regarding Permutations Which Destroy Arithmetic Progressions

abstract: Hegarty conjectured for $n\neq 2, 3, 5, 7$ that $\mathbb{Z}/n\mathbb{Z}$ has a permutation which destroys all arithmetic progressions mod $n$. For $n\ge n_0$, Hegarty and Martinsson demonstrated that $\mathbb{Z}/n\mathbb{Z}$ has an arithmetic-progression destroying permutation. However $n_0\approx 1.4\times 10^{14}$ and thus resolving the conjecture in full remained out of reach of any computational techniques. However, this paper using constructions modeled after those used by Elkies and Swaminathan for the case of $\mathbb{Z}/p\mathbb{Z}$ with $p$ being prime, establish the conjecture in full. Furthermore our results do not rely on the fact that it suffices to study when $n<n_0$ and thus our results completely independent of the proof given by Hegarty and Martinsson.


title: Semiparametric Efficiency in Convexity Constrained Single Index Model

abstract: We consider estimation and inference in a single index regression model with an unknown convex link function. We introduce a convex and Lipschitz constrained least squares estimator (CLSE) for both the parametric and the nonparametric components given independent and identically distributed observations. We prove the consistency and find the rates of convergence of the CLSE when the errors are assumed to have only $q \ge 2$ moments and are allowed to depend on the covariates. When $q\ge 5$, we establish $n^{-1/2}$-rate of convergence and asymptotic normality of the estimator of the parametric component. Moreover, the CLSE is proved to be semiparametrically efficient if the errors happen to be homoscedastic. {We develop and implement a numerically stable and computationally fast algorithm to compute our proposed estimator in the R package~\texttt{simest}}. We illustrate our methodology through extensive simulations and data analysis. Finally, our proof of efficiency is geometric and provides a general framework that can be used to prove efficiency of estimators in a wide variety of semiparametric models even when they do not satisfy the efficient score equation directly.


title: Large-Scale Low-Rank Matrix Learning with Nonconvex Regularizers

abstract: Low-rank modeling has many important applications in computer vision and machine learning. While the matrix rank is often approximated by the convex nuclear norm, the use of nonconvex low-rank regularizers has demonstrated better empirical performance. However, the resulting optimization problem is much more challenging. Recent state-of-the-art requires an expensive full SVD in each iteration. In this paper, we show that for many commonly-used nonconvex low-rank regularizers, a cutoff can be derived to automatically threshold the singular values obtained from the proximal operator. This allows such operator being efficiently approximated by power method. Based on it, we develop a proximal gradient algorithm (and its accelerated variant) with inexact proximal splitting and prove that a convergence rate of O(1/T) where T is the number of iterations is guaranteed. Furthermore, we show the proposed algorithm can be well parallelized, which achieves nearly linear speedup w.r.t the number of threads. Extensive experiments are performed on matrix completion and robust principal component analysis, which shows a significant speedup over the state-of-the-art. Moreover, the matrix solution obtained is more accurate and has a lower rank than that of the nuclear norm regularizer.


title: Adiabatic control of surface plasmon-polaritons in a 3-layers graphene curved configuration

abstract: In this paper, we utilize coupled mode theory (CMT) to model the coupling between surface plasmon-polaritons (SPPs) between multiple graphene sheets. By using the Stimulated Raman Adiabatic Passage (STIRAP) Quantum Control Technique, we propose a novel directional coupler based on SPPs evolution in three layers of graphene sheets in some curved configuration. Our calculated results show that the SPP can be transferred efficiently from the input graphene sheet to the output graphene sheet, and the coupling is also robust that it is not sensitive to the length of the device.


title: Tameness in least fixed-point logic and McColm's conjecture

abstract: We investigate four model-theoretic tameness properties in the context of least fixed-point logic over a family of finite structures. We find that each of these properties depends only on the elementary (i.e., first-order) limit theory, and we completely determine the valid entailments among them. In contrast to the context of first-order logic on arbitrary structures, the order property and independence property are equivalent in this setting. McColm conjectured that least fixed-point definability collapses to first-order definability exactly when proficiency fails. McColm's conjecture is known to be false in general. However, we show that McColm's conjecture is true for any family of finite structures whose limit theory is model-theoretically tame.


title: Adaptive Hierarchical Clustering Using Ordinal Queries

abstract: In many applications of clustering (for example, ontologies or clusterings of animal or plant species), hierarchical clusterings are more descriptive than a flat clustering. A hierarchical clustering over $n$ elements is represented by a rooted binary tree with $n$ leaves, each corresponding to one element. The subtrees rooted at interior nodes capture the clusters. In this paper, we study active learning of a hierarchical clustering using only ordinal queries. An ordinal query consists of a set of three elements, and the response to a query reveals the two elements (among the three elements in the query) which are "closer" to each other than to the third one. We say that elements $x$ and $x'$ are closer to each other than $x"$ if there exists a cluster containing $x$ and $x'$, but not $x"$. When all the query responses are correct, there is a deterministic algorithm that learns the underlying hierarchical clustering using at most $n \log_2 n$ adaptive ordinal queries. We generalize this algorithm to be robust in a model in which each query response is correct independently with probability $p > \frac{1}{2}$, and adversarially incorrect with probability $1 - p$. We show that in the presence of noise, our algorithm outputs the correct hierarchical clustering with probability at least $1 - \delta$, using $O(n \log n + n \log(1/\delta))$ adaptive ordinal queries. For our results, adaptivity is crucial: we prove that even in the absence of noise, every non-adaptive algorithm requires $\Omega(n^3)$ ordinal queries in the worst case.