title: Offset-Based Beamforming: A New Approach to Robust Downlink Transmission

abstract: The design of a set of beamformers for the multiuser multiple-input single-output (MISO) downlink that provides the receivers with prespecified levels of quality-of-service (QoS) can be quite challenging when the channel state information is not perfectly known at the base station. The constraint of having the SINR meet or exceed a given threshold with high probability is intractable in general, which results in problems that are fundamentally hard to solve. In this paper, we will develop a high quality approximation of the SINR outage constraint that, along with a semidefinite relaxation, enables us to formulate the beamformer design problem as a convex optimization problem that can be efficiently solved. For systems in which the uncertainty size is small, a further approximation yields algorithms based on iterative evaluations of closed-form expressions that have substantially lower computational cost. Since finding the beamforming directions incurs most of the computational load of these algorithms, analogous power loading algorithms for predefined beamforming directions are developed and their performance is shown to be close to optimal. When the system contains a large number of antennas, the proposed power loading can be obtained at a computational cost that grows only linearly in the number of antennas. The proposed power loading algorithm provides an explicit relationship between the outage probability required and the power consumed, which allows us to precisely control the power consumption, and automatically identifies users who are consuming most of the power resources. The flexibility of the proposed approach is illustrated by developing a power loading technique that minimizes an average notion of outage.


title: Optimal Shrinkage of Singular Values Under Random Data Contamination

abstract: A low rank matrix X has been contaminated by uniformly distributed noise, missing values, outliers and corrupt entries. Reconstruction of X from the singular values and singular vectors of the contaminated matrix Y is a key problem in machine learning, computer vision and data science. In this paper we show that common contamination models (including arbitrary combinations of uniform noise,missing values, outliers and corrupt entries) can be described efficiently using a single framework. We develop an asymptotically optimal algorithm that estimates X by manipulation of the singular values of Y , which applies to any of the contamination models considered. Finally, we find an explicit signal-to-noise cutoff, below which estimation of X from the singular value decomposition of Y must fail, in a well-defined sense.


title: FashionBrain Project: A Vision for Understanding Europe's Fashion Data Universe

abstract: A core business in the fashion industry is the understanding and prediction of customer needs and trends. Search engines and social networks are at the same time a fundamental bridge and a costly middleman between the customer's purchase intention and the retailer. To better exploit Europe's distinctive characteristics e.g., multiple languages, fashion and cultural differences, it is pivotal to reduce retailers' dependence to search engines. This goal can be achieved by harnessing various data channels (manufacturers and distribution networks, online shops, large retailers, social media, market observers, call centers, press/magazines etc.) that retailers can leverage in order to gain more insight about potential buyers, and on the industry trends as a whole. This can enable the creation of novel on-line shopping experiences, the detection of influencers, and the prediction of upcoming fashion trends. In this paper, we provide an overview of the main research challenges and an analysis of the most promising technological solutions that we are investigating in the FashionBrain project.


title: Evaluating Password Advice

abstract: Password advice is constantly circulated by standards agencies, companies, websites and specialists. But there appears to be great diversity in terms of the advice that is given. Users have noticed that different websites are enforcing different restrictions. For example, requiring different combinations of uppercase and lowercase letters, numbers and special characters. We collected password advice and found that the advice distributed by one organization can directly contradict advice given by another. Our paper aims to illuminate interesting characteristics for a sample of the password advice distributed. We also create a framework for identifying the costs associated with implementing password advice. In doing so we identify a reason for why password advice is often both derided and ignored.


title: Process interpretation of current entropic bounds

abstract: We show for Markov diffusion processes that the quadratic entropic bound, recently derived for the rate functions of nonequilibrium currents, can be seen as being produced by an effective process that creates current fluctuations in a sub-optimal way by modifying only the non-reversible part of the drift or force of the process considered while keeping its reversible part constant. This provides a clear interpretation of the bound in terms of a physical process, which explains, among other things, its relation to the fluctuation relation, linear response, and reversible limits. The existence of more general quadratic bounds, and related uncertainty relations, for physical quantities other than currents is also discussed.


title: Structural Variability from Noisy Tomographic Projections

abstract: In cryo-electron microscopy, the 3D electric potentials of an ensemble of molecules are projected along arbitrary viewing directions to yield noisy 2D images. The volume maps representing these potentials typically exhibit a great deal of structural variability, which is described by their 3D covariance matrix. Typically, this covariance matrix is approximately low-rank and can be used to cluster the volumes or estimate the intrinsic geometry of the conformation space. We formulate the estimation of this covariance matrix as a linear inverse problem, yielding a consistent least-squares estimator. For $n$ images of size $N$-by-$N$ pixels, we propose an algorithm for calculating this covariance estimator with computational complexity $\mathcal{O}(nN^4+\sqrt{\kappa}N^6 \log N)$, where the condition number $\kappa$ is empirically in the range $10$--$200$. Its efficiency relies on the observation that the normal equations are equivalent to a deconvolution problem in 6D. This is then solved by the conjugate gradient method with an appropriate circulant preconditioner. The result is the first computationally efficient algorithm for consistent estimation of 3D covariance from noisy projections. It also compares favorably in runtime with respect to previously proposed non-consistent estimators. Motivated by the recent success of eigenvalue shrinkage procedures for high-dimensional covariance matrices, we introduce a shrinkage procedure that improves accuracy at lower signal-to-noise ratios. We evaluate our methods on simulated datasets and achieve classification results comparable to state-of-the-art methods in shorter running time. We also present results on clustering volumes in an experimental dataset, illustrating the power of the proposed algorithm for practical determination of structural variability.


title: Supersymmetric Gauge Theory with Space-time Dependent Couplings

abstract: We study deformations of ${\cal N}=4$ supersymmetric Yang-Mills theory with couplings and masses depending on space-time. The conditions to preserve part of the supersymmetry are derived and a lot of solutions of these conditions are found. The main example is the case with $ISO(1,1)\times SO(3)\times SO(3)$ symmetry, in which couplings, as well as masses and the theta parameter, can depend on two spatial coordinates. In the case $ISO(1,1)$ is enhanced to $ISO(1,2)$, it reproduces the supersymmetric Janus configuration found by Gaiotto and Witten. When $SO(3)\times SO(3)$ is enhanced to $SO(6)$, it agrees with the world-volume theory of D3-branes embedded in F-theory (a background with 7-branes in type IIB string theory). We also find the general solution of the supersymmetry conditions for the cases with $ISO(1,1)\times SO(2)\times SO(4)$ symmetry. The cases with time dependent couplings and/or masses are also considered.


title: Statistical Inference on Tree Swallow Migrations with Random Forests

abstract: Bird species' migratory patterns have typically been studied through individual observations and historical records. In recent years however, the eBird citizen science project, which solicits observations from thousands of bird watchers around the world, has opened the door for a data-driven approach to understanding the large-scale geographical movements. Here, we focus on the North American Tree Swallow (\textit{Tachycineta bicolor}) occurrence patterns throughout the eastern United States. Migratory departure dates for this species are widely believed by both ornithologists and casual observers to vary substantially across years, but the reasons for this are largely unknown. In this work, we present evidence that maximum daily temperature is a major factor influencing Tree Swallow occurrence. Because it is generally understood that species occurrence is a function of many complex, high-order interactions between ecological covariates, we utilize the flexible modeling approach offered by random forests. Making use of recent asymptotic results, we provide formal hypothesis tests for predictive significance various covariates and also develop and implement a permutation-based approach for formally assessing interannual variations by treating the prediction surfaces generated by random forests as functional data. Each of these tests suggest that maximum daily temperature has a significant effect on migration patterns.


title: A strong collapse increasing the geometric simplicial Lusternik-Schnirelmann category

abstract: In [3], after defining notions of LS category in the simplicial context, the authors show that the geometric simplicial LS category is non-decreasing under strong collapses. However, they do not give examples where it increases strictly, but they conjecture that such an example should exist, and thus that the geometric simplicial LS category is not strong homotopy invariant. The purpose of this note is to provide with such an example. We construct a simplicial complex whose simplicial and geometric simplicial LS categories are different, and using this, we provide an example of a strong collapse that increases the geometric simplicial LS category, thus settling the geometric simplicial LS category not being strong homotopy invariant.


title: Synchronization Strings: Explicit Constructions, Local Decoding, and Applications

abstract: This paper gives new results for synchronization strings, a powerful combinatorial object that allows to efficiently deal with insertions and deletions in various communication settings: $\bullet$ We give a deterministic, linear time synchronization string construction, improving over an $O(n^5)$ time randomized construction. Independently of this work, a deterministic $O(n\log^2\log n)$ time construction was just put on arXiv by Cheng, Li, and Wu. We also give a deterministic linear time construction of an infinite synchronization string, which was not known to be computable before. Both constructions are highly explicit, i.e., the $i^{th}$ symbol can be computed in $O(\log i)$ time. $\bullet$ This paper also introduces a generalized notion we call long-distance synchronization strings that allow for local and very fast decoding. In particular, only $O(\log^3 n)$ time and access to logarithmically many symbols is required to decode any index. We give several applications for these results: $\bullet$ For any $\delta<1$ and $\epsilon>0$ we provide an insdel correcting code with rate $1-\delta-\epsilon$ which can correct any $O(\delta)$ fraction of insdel errors in $O(n\log^3n)$ time. This near linear computational efficiency is surprising given that we do not even know how to compute the (edit) distance between the decoding input and output in sub-quadratic time. We show that such codes can not only efficiently recover from $\delta$ fraction of insdel errors but, similar to [Schulman, Zuckerman; TransInf'99], also from any $O(\delta/\log n)$ fraction of block transpositions and replications. $\bullet$ We show that highly explicitness and local decoding allow for infinite channel simulations with exponentially smaller memory and decoding time requirements. These simulations can be used to give the first near linear time interactive coding scheme for insdel errors.


title: Topology in the $SU(N_f)$ chiral symmetry restored phase of unquenched QCD and axion cosmology

abstract: The axion is one of the more interesting candidates to make the dark matter of the universe, and the axion potential plays a fundamental role in the determination of the dynamics of the axion field. Moreover, the way in which the $U(1)_A$ anomaly manifests itself in the chiral symmetry restored phase of $QCD$ at high temperature could be tested when probing the $QCD$ phase transition in relativistic heavy ion collisions. With these motivations, we investigate the physical consequences of the survival of the effects of the $U(1)_A$ anomaly in the chiral symmetric phase of $QCD$, and show that the free energy density is a singular function of the quark mass $m$, in the chiral limit, and that the $\sigma$ and $\bar\pi$ susceptibilities diverge in this limit at any $T\ge T_c$. We also show that the difference between the $\bar\pi$ and $\bar\delta$ susceptibilities diverges in the chiral limit at any $T\ge T_c$, a result that can be contrasted with the existing lattice calculations; and discuss on the generalization of these results to the $N_f\ge 3$ model.


title: Interference Queueing Networks on Grids

abstract: Consider a countably infinite collection of interacting queues, with a queue located at each point of the $d$-dimensional integer grid, having independent Poisson arrivals, but dependent service rates. The service discipline is of the processor sharing type,with the service rate in each queue slowed down, when the neighboring queues have a larger workload. The interactions are translation invariant in space and is neither of the Jackson Networks type, nor of the mean-field type. Coupling and percolation techniques are first used to show that this dynamics has well defined trajectories. Coupling from the past techniques are then proposed to build its minimal stationary regime. The rate conservation principle of Palm calculus is then used to identify the stability condition of this system, where the notion of stability is appropriately defined for an infinite dimensional process. We show that the identified condition is also necessary in certain special cases and conjecture it to be true in all cases. Remarkably, the rate conservation principle also provides a closed form expression for the mean queue size. When the stability condition holds, this minimal solution is the unique translation invariant stationary regime. In addition, there exists a range of small initial conditions for which the dynamics is attracted to the minimal regime. Nevertheless, there exists another range of larger though finite initial conditions for which the dynamics diverges, even though stability criterion holds.


title: Lip2AudSpec: Speech reconstruction from silent lip movements video

abstract: In this study, we propose a deep neural network for reconstructing intelligible speech from silent lip movement videos. We use auditory spectrogram as spectral representation of speech and its corresponding sound generation method resulting in a more natural sounding reconstructed speech. Our proposed network consists of an autoencoder to extract bottleneck features from the auditory spectrogram which is then used as target to our main lip reading network comprising of CNN, LSTM and fully connected layers. Our experiments show that the autoencoder is able to reconstruct the original auditory spectrogram with a 98% correlation and also improves the quality of reconstructed speech from the main lip reading network. Our model, trained jointly on different speakers is able to extract individual speaker characteristics and gives promising results of reconstructing intelligible speech with superior word recognition accuracy.


title: Numerical Method for the Maxwell-Liouville-von Neumann Equations using Efficient Matrix Exponential Computations

abstract: We present a novel method to solve the Maxwell-Liouville-von Neumann (MLN) equations in an accurate and efficient way without invoking the rotating wave approximation (RWA). The method is a combination of two established concepts, namely the operator splitting method as well as the adjoint representation of the Lie algebra SU(N) (or pseudospin representation). The former concept ensures the accuracy of the approach, but is computationally expensive. The latter concept provides an efficient representation of the problem and two optimization possibilities. We have implemented and verified both optimization approaches and demonstrate that substantial speedup can be achieved.


title: Correction of the orbital mass of double galaxies estimation

abstract: We obtain a more accurate statistical estimation of the mass of double galaxies moving in circular orbits, including confidence intervals for different confidence levels.


abstract: Morse foliations of codimension one on the sphere S^3 are studied and the existence of special components for these foliations is derived. As a corollary the instability of Morse foliations can be proven in almost all cases.


title: Invariant linear functionals on $L^{\infty}(\mathbb{R}_+)$

abstract: We consider a continuous version of the classical notion of Banach limits, namely, positive linear functionals on $L^{\infty}(\mathbb{R}_+)$ invariant under translations $f(x) \mapsto f(x+s)$ of $L^{\infty}(\mathbb{R}_+)$ for every $s \ge 0$. We give its characterization in terms of the invariance under the operation of a certain linear mapping on $L^{\infty}(\mathbb{R}_+)$. Applications to summability methods are provided in the last section.


title: Subnanometer localization accuracy in widefield optical microscopy

abstract: The common assumption that precision is the limit of accuracy in localization microscopy and the typical absence of comprehensive calibration of optical microscopes lead to a widespread issue - overconfidence in measurement results with nanoscale statistical uncertainties that can be invalid due to microscale systematic errors. In this article, we report a comprehensive solution to this underappreciated problem. We develop arrays of subresolution apertures into the first reference materials that enable localization errors approaching the atomic scale across a submillimeter field. We present novel methods for calibrating our microscope system using aperture arrays and develop aberration corrections that reach the precision limit of our reference materials. We correct and register localization data from multiple colors and test different sources of light emission with equal accuracy, indicating the general applicability of our reference materials and calibration methods. In a first application of our new measurement capability, we introduce the concept of critical dimension localization microscopy, facilitating tests of nanofabrication processes and quality control of aperture arrays. In a second application, we apply these stable reference materials to answer open questions about the apparent instability of fluorescent nanoparticles that commonly serve as fiducial markers. Our study establishes a foundation for subnanometer localization accuracy in widefield optical microscopy.


title: Loop-TNR analysis of CP(1) model with theta term

abstract: The phase structure of the two dimensional lattice CP(1) model in the presence of the $\theta$ term is analyzed by tensor network methods. The tensor renormalization group, which is a standard renormalization method of tensor networks, is used for the regions $\theta =0$ and $\theta \neq 0$. Loop-TNR, which is more suitable for the analysis of near criticality, is also implemented for the region $\theta =0$. The application of Loop-TNR for the region $\theta \neq 0$ is left for future work.


title: Improving Negative Sampling for Word Representation using Self-embedded Features

abstract: Although the word-popularity based negative sampler has shown superb performance in the skip-gram model, the theoretical motivation behind oversampling popular (non-observed) words as negative samples is still not well understood. In this paper, we start from an investigation of the gradient vanishing issue in the skipgram model without a proper negative sampler. By performing an insightful analysis from the stochastic gradient descent (SGD) learning perspective, we demonstrate that, both theoretically and intuitively, negative samples with larger inner product scores are more informative than those with lower scores for the SGD learner in terms of both convergence rate and accuracy. Understanding this, we propose an alternative sampling algorithm that dynamically selects informative negative samples during each SGD update. More importantly, the proposed sampler accounts for multi-dimensional self-embedded features during the sampling process, which essentially makes it more effective than the original popularity-based (one-dimensional) sampler. Empirical experiments further verify our observations, and show that our fine-grained samplers gain significant improvement over the existing ones without increasing computational complexity.


title: Minimum Circuit Size, Graph Isomorphism, and Related Problems

abstract: We study the computational power of deciding whether a given truth-table can be described by a circuit of a given size (the Minimum Circuit Size Problem, or MCSP for short), and of the variant denoted as MKTP where circuit size is replaced by a polynomially-related Kolmogorov measure. All prior reductions from supposedly-intractable problems to MCSP / MKTP hinged on the power of MCSP / MKTP to distinguish random distributions from distributions produced by hardness-based pseudorandom generator constructions. We develop a fundamentally different approach inspired by the well-known interactive proof system for the complement of Graph Isomorphism (GI). It yields a randomized reduction with zero-sided error from GI to MKTP. We generalize the result and show that GI can be replaced by any isomorphism problem for which the underlying group satisfies some elementary properties. Instantiations include Linear Code Equivalence, Permutation Group Conjugacy, and Matrix Subspace Conjugacy. Along the way we develop encodings of isomorphism classes that are efficiently decodable and achieve compression that is at or near the information-theoretic optimum; those encodings may be of independent interest.


title: Mixed analytical-stochastic simulation method for the recovery of a Brownian gradient source from probability fluxes to small windows

abstract: Is it possible to recover the position of a source from the steady-state fluxes of Brownian particles to small absorbing windows located on the boundary of a domain? To address this question, we develop a numerical procedure to avoid tracking Brownian trajectories in the entire infinite space. Instead, we generate particles near the absorbing windows, computed from the analytical expression of the exit probability. When the Brownian particles are generated by a steady-state gradient at a single point, we compute asymptotically the fluxes to small absorbing holes distributed on the boundary of half-space and on a disk in two dimensions, which agree with stochastic simulations. We also derive an expression for the splitting probability between small windows using the matched asymptotic method. Finally, when there are more than two small absorbing windows, we show how to reconstruct the position of the source from the diffusion fluxes. The present approach provides a computational first principle for the mechanism of sensing a gradient of diffusing particles, a ubiquitous problem in cell biology.


title: CMB constraints on $\beta$-exponential inflationary models

abstract: We analyze a class of generalized inflationary models proposed in Ref. [1], known as $\beta$-exponential inflation. We show that this kind of potential can arise in the context of brane cosmology, where the field describing the size of the extra-dimension is interpreted as the inflaton. We discuss the observational viability of this class of model in light of the latest Cosmic Microwave Background (CMB) data from the Planck Collaboration through a Bayesian analysis, and impose tight constraints on the model parameters. We find that the CMB data alone prefer \textit{weakly} the minimal standard model ($\Lambda$CDM) over the $\beta$-exponential inflation. However, when current local measurements of the Hubble parameter, $H_0$, are considered, the $\beta$-inflation model is {\it moderately} preferred over the $\Lambda$CDM cosmology, making the study of this class of inflationary models interesting in the context of the current $H_0$ tension.


title: Joint Screening Tests for LASSO

abstract: This paper focusses on "safe" screening techniques for the LASSO problem. Motivated by the need for low-complexity algorithms, we propose a new approach, dubbed "joint" screening test, allowing to screen a set of atoms by carrying out one single test. The approach is particularized to two different sets of atoms, respectively expressed as sphere and dome regions. After presenting the mathematical derivations of the tests, we elaborate on their relative effectiveness and discuss the practical use of such procedures.


title: The Experimental Program for High Pressure Gas Filled Radio Frequency Cavities for Muon Cooling Channels

abstract: An intense beam of muons is needed to provide a luminosity on the order of 10$^{34}$ cm$^{-2}$s$^{-1}$ for a multi-TeV collider. Because muons produced by colliding a multi-MW proton beam with a target made of carbon or mercury have a large phase space, significant six dimensional cooling is required. Through ionization cooling - the only cooling method that works within the lifetime of the muon - and emittance exchange, the desired emittances for a Higgs Factory or higher energy collider are attainable. A cooling channel utilizing gas filled radio frequency cavities has been designed to deliver the requisite cool muon beam. Technology development of these RF cavities has progressed from breakdown studies, through beam tests, to dielectric loaded and reentrant cavity designs. The results of these experiments are summarized.


title: Spatial behavior in a Mott insulator near the voltage-driven resistive transition

abstract: We develop a real space theory of the voltage bias driven transition from a Mott insulator to a correlated metal. Within our Keldysh mean field approach the problem reduces to a self-consistency scheme for the charge and spin profiles in this open system. We solve this problem for a two dimensional antiferromagnetic Mott insulator at zero temperature. The charge and spin magnitude is uniform over the system at zero bias, but a bias $V$ leads to spatial modulation over a lengthscale $\xi(V)$ near the edges. $\xi(V)$ grows rapidly and becomes comparable to system size as $V$ increases towards a threshold scale $V_c$. The linear response conductance of the insulator is zero with the current being exponentially small for $V \ll V_c$. The current increases rapidly as $V \rightarrow V_c$. Beyond $V_c$, we observe an inhomogeneous low moment antiferromagnetic metal, and at even larger bias a current saturated paramagnetic metal. We suggest an approximate scheme for the spectral features of this nonequilibrium system.


title: Coulomb effects on the photo-excited quantum dynamics of electrons in a plasmonic nanosphere

abstract: With recent experiments investigating the optical properties of progressively smaller plasmonic particles, quantum effects become increasingly more relevant, requiring a microscopic description. Using the density matrix formalism we analyze the photo-excited few-electron dynamics of a small nanosphere. Following the standard derivation of the bulk plasmon we particularly aim on elucidating the role of the Coulomb interaction. Calculating the dielectric susceptibility spectrum in the linear optical response we find discrete resonances resulting from a collective response mediated by the Coulomb interaction between the electrons. In the nonlinear regime, the occupations of the system exhibit oscillations between the interacting eigenstates. Our approach provides an ideal platform to study and explain nonlinear and quantum plasmonics, revealing that the photo-excited dynamics of plasmonic nano spheres has similarities with and combines characteristics of both, the well-known two-level Rabi dynamics, and the collective many-electron behavior typical of plasmons.


title: Sparse Diffusion-Convolutional Neural Networks

abstract: The predictive power and overall computational efficiency of Diffusion-convolutional neural networks make them an attractive choice for node classification tasks. However, a naive dense-tensor-based implementation of DCNNs leads to $\mathcal{O}(N^2)$ memory complexity which is prohibitive for large graphs. In this paper, we introduce a simple method for thresholding input graphs that provably reduces memory requirements of DCNNs to O(N) (i.e. linear in the number of nodes in the input) without significantly affecting predictive performance.


title: Linear convergence of the generalized Douglas-Rachford algorithm for feasibility problems

abstract: In this paper, we study the generalized Douglas-Rachford algorithm and its cyclic variants which include many projection-type methods such as the classical Douglas-Rachford algorithm and the alternating projection algorithm. Specifically, we establish several local linear convergence results for the algorithm in solving feasibility problems with finitely many closed possibly nonconvex sets under different assumptions. Our findings not only relax some regularity conditions but also improve linear convergence rates in the literature. In the presence of convexity, the linear convergence is global.


title: Coherent transport structures in magnetized plasmas, I : Theory

abstract: In a pair of linked articles (called Article I and II respectively) we apply the concept of Lagrangian Coherent Structures (LCSs) borrowed from the study of Dynamical Systems to magnetic field configurations in order to separate regions where field lines have different kind of behaviour. In the present article, article I, after recalling the definition and the properties of the LCSs, we show how this conceptual framework can be applied to the study of particle transport in a magnetized plasma. Futhermore we introduce a simplified model that allows us to consider explicitly the case where the magnetic configuration evolves in time on timescales comparable to the particle transit time through the configuration. In contrast with previous works on this topic, this analysis requires that a system that is aperiodic in time be investigated. In this case the Poincar\'e map technique cannot be applied and LCSs remain the only viable tool.


title: Acoustic response in a a one-dimensional layered pseudo-Hermitian metamaterial containing defects

abstract: Using transfer-matrix methods, we investigate the response of a multilayered metamaterial system containing defects to an incident acoustic plane wave at normal or oblique incidence. The transmission response is composed of pass-bands with oscillatory behaviour, separated by band gaps and covers a wide frequency range. The presence of gain and loss in the layers leads to the emergence of symmetry breaking and re-entrant phases. In the general case, a system containing defects will display a more general property, pseudo-Hermiticity (PH), of which $\mathcal{PT}$ systems are a subset. In the PH-symmetric phase, unidirectional responses of the reflection, accomplished by reversing the parity $\mathcal{P}$, can be found but the response sometimes deviates from the predictions of simple scattering theory which call for a pseudo-unitarity relation relating the transmission and the two directions of reflections to hold. The converse of reversing the parity, reversing the time operator $\mathcal{T}$ in a spatially-asymmetric system within the PH-symmetric regime can lead to different transmissions: a pass-band versus a stop-band. As regions of stable PH-symmetric pass-band transmission oscillations occur over a wide spectral range, there is a large flexibility in system parameters such as layer thicknesses, for leading to the desired unidirectional traits. In addition, we find that while defects in general lead to a near or complete loss of PH symmetry at all frequencies, they can be exploited to produce highly-sensitive responses, making such systems good candidates for sensor applications.


title: $\overline{B}\rightarrow D^\ast\ell\overline{\nu}$ at non-zero recoil

abstract: We present preliminary results from our analysis of the form factors for the $\overline{B}\rightarrow D^\ast\ell\overline{\nu}$ decay at non-zero recoil. Our analysis includes 15 MILC asqtad ensembles with $N_f=2+1$ flavors of sea quarks and lattice spacings ranging from $a\approx 0.15$ fm down to $0.045$ fm. The valence light quarks employ the asqtad action, whereas the heavy quarks are treated using the Fermilab action. We conclude with a discussion of future plans and phenomenological implications. When combined with experimental measurements of the decay rate, our calculation will enable a determination of the CKM matrix element $|V_{cb}|$.


title: Signatures of graviton masses on the CMB

abstract: The impact of the existence of gravitons with non-vanishing masses on the B modes of the Cosmic Microwave Background (CMB) is investigated. We also focus on putative modifications to the speed of the gravitational waves. We find that a change of the graviton speed shifts the acoustic peaks of the CMB and then could be easily constrained. For the case of massive gravity, we show analytically how the B modes are sourced in a manner differing from the massless case leading to a plateau at low $l$ in the CMB spectrum. We also study the case when there are more than one graviton, and when pressure instabilities are present. The latter would occur in doubly coupled bigravity in the radiation era. We focus on the case where a massless graviton becomes tachyonic in the radiation era whilst a massive one remains stable. As the unstable mode decouples from matter in the radiation era, we find that the effects of the instability is largely reduced on the spectrum of B-modes as long as the unstable graviton does not grow into the non-linear regime. In all cases when both massless and massive gravitons are present, we find that the B-mode CMB spectrum is characterised by a low $l$ plateau together with a shifted position for the first few peaks compared to a purely massive graviton spectrum, a shift which depends on the mixing between the gravitons in their coupling to matter and could serve as a hint in favour of the existence of multiple gravitons.


title: Algorithms for Deforming and Contracting Simply Connected Discrete Closed Manifolds (III)

abstract: In a recent paper, {\it Algorithms for Deforming and Contracting Simply Connected Discrete Closed Manifolds (II)}, we discussed two algorithms for deforming and contracting a simply connected discrete closed manifold into a discrete sphere. The first algorithm was a continuation of work that began in {\it Algorithms for Deforming and Contracting Simply Connected Discrete Closed Manifolds (I)}, the second algorithm contained a more direct treatment of contraction for discrete manifolds. In this paper, we clarify that we can use this same method on standard piecewise linear (PL) complexes on the triangulation of general smooth manifolds. Our discussion is based on triangulation techniques invented by Cairns, Whitehead, and Whitney more than half of a century ago. In this paper, we use PL or simplicial complexes to replace certain concepts of discrete manifolds in previous papers. Note that some details in the original papers related to discrete manifolds may also need to be slightly modified for this purpose. In this paper, we use the algorithmic procedure (a contraction process) to prove the following theorem: For a finite triangulation of a simply-connected closed and orientable 3-manifold $M$ in Euclidean space, if a (simply-connected closed) 2-cycle which was made by 2-cells of this triangulation separates $M$ into two connected components, then each of the components will also be simply-connected. In addition, we can algorithmically make $M$ to be homeomorphic to a 3-sphere. The relationship of the theorem to a generalized Jordan separation problem, the general Jordan-Schoenflies theorem, and other important problems are also discussed. This revision adds more figures and explanations as well as two more special cases. We will post the detailed algorithm/procedure for practical triangulations in the following papers.


title: Spiking Optical Flow for Event-based Sensors Using IBM's TrueNorth Neurosynaptic System

abstract: This paper describes a fully spike-based neural network for optical flow estimation from Dynamic Vision Sensor data. A low power embedded implementation of the method which combines the Asynchronous Time-based Image Sensor with IBM's TrueNorth Neurosynaptic System is presented. The sensor generates spikes with sub-millisecond resolution in response to scene illumination changes. These spike are processed by a spiking neural network running on TrueNorth with a 1 millisecond resolution to accurately determine the order and time difference of spikes from neighboring pixels, and therefore infer the velocity. The spiking neural network is a variant of the Barlow Levick method for optical flow estimation. The system is evaluated on two recordings for which ground truth motion is available, and achieves an Average Endpoint Error of 11% at an estimated power budget of under 80mW for the sensor and computation.


title: Cell Line Classification Using Electric Cell-substrate Impedance Sensing (ECIS)

abstract: We consider cell line classification using multivariate time series data obtained from electric cell-substrate impedance sensing (ECIS) technology. The ECIS device, which monitors the attachment and spreading of mammalian cells in real time through the collection of electrical impedance data, has historically been used to study one cell line at a time. However, we show that if applied to data from multiple cell lines, ECIS can be used to classify unknown or potentially mislabeled cells, which may help to mitigate the current crisis of reproducibility in the biological literature. We assess a range of approaches to this new problem, testing different classification methods and deriving a dictionary of 29 features to characterize ECIS data. Our analysis also makes use of simultaneous multi-frequency ECIS data, where previous studies have focused on only one frequency. In classification tests on fifteen mammalian cell lines, we obtain very high out-of-sample accuracy. These preliminary findings provide a baseline for future large-scale studies in this field.


title: The Brown-Peterson spectrum is not $\mathbb{E}_{2(p^2+2)}$ at odd primes

abstract: We show that the odd-primary Brown-Peterson spectrum $\mathrm{BP}$ does not admit the structure of an $\mathbb{E}_{2(p^2+2)}$ ring spectrum and that there can be no map $\mathrm{MU} \to \mathrm{BP}$ of $\mathbb{E}_{2p+3}$ ring spectra. We also prove the same results for truncated Brown-Peterson spectra $\mathrm{BP} \langle n \rangle$ of height $n \geq 4$. This extends results of Lawson at the prime $2$.


title: Angle Dependent Interlayer Magnetoresistance (ILMR) in Multilayer Graphene Stacks

abstract: Interlayer magnetoresistance (ILMR) effect is explored in a vertical stack of weakly coupled multilayer graphene as grown by chemical vapor deposition (CVD). This effect has been characterized as a function of temperature and tilt angle of the magnetic field with respect to the interlayer current. To our knowledge, this is the first experimental report on angle dependent ILMR effect in graphitic systems. Our data agrees qualitatively with the existing theories of ILMR in multilayer massless Dirac Fermion systems. However, a sharper change in ILMR has been observed as the tilt angle of the magnetic field is varied. A physical explanation of this effect is proposed, which is consistent with our experimental scenario.


title: Klout Topics for Modeling Interests and Expertise of Users Across Social Networks

abstract: This paper presents Klout Topics, a lightweight ontology to describe social media users' topics of interest and expertise. Klout Topics is designed to: be human-readable and consumer-friendly; cover multiple domains of knowledge in depth; and promote data extensibility via knowledge base entities. We discuss why this ontology is well-suited for text labeling and interest modeling applications, and how it compares to available alternatives. We show its coverage against common social media interest sets, and examples of how it is used to model the interests of over 780M social media users on Klout.com. Finally, we open the ontology for external use.


title: On the role of synaptic stochasticity in training low-precision neural networks

abstract: Stochasticity and limited precision of synaptic weights in neural network models are key aspects of both biological and hardware modeling of learning processes. Here we show that a neural network model with stochastic binary weights naturally gives prominence to exponentially rare dense regions of solutions with a number of desirable properties such as robustness and good generalization performance, while typical solutions are isolated and hard to find. Binary solutions of the standard perceptron problem are obtained from a simple gradient descent procedure on a set of real values parametrizing a probability distribution over the binary synapses. Both analytical and numerical results are presented. An algorithmic extension aimed at training discrete deep neural networks is also investigated.


title: Supersymmetric $\mathrm{AdS}_3$ supergravity backgrounds and holography

abstract: We analyse the conditions for $\mathrm{AdS}_3 \times \mathcal{M}_7$ backgrounds with pure NS-NS flux to be supersymmetric. We find that a necessary condition is that $\mathcal{M}_7$ is a $\mathrm{U}(1)$-fibration over a balanced manifold. We subsequently classify all $\mathcal{N}=(2,2)$ solutions where $\mathcal{M}_7$ satisfies the stronger condition of being a $\mathrm{U}(1)$-fibration over a Kahler manifold. We compute the BPS spectrum of all the backgrounds in this classification. We assign a natural dual CFT to the backgrounds and confirm that the BPS spectra agree, thus providing evidence in favour of the proposal.


title: Diffusive hydrodynamics of out-of-time-ordered correlators with charge conservation

abstract: The scrambling of quantum information in closed many-body systems, as measured by out-of-time-ordered correlation functions (OTOCs), has lately received considerable attention. Recently, a hydrodynamical description of OTOCs has emerged from considering random local circuits, aspects of which are conjectured to be universal to ergodic many-body systems, even without randomness. Here we extend this approach to systems with locally conserved quantities (e.g., energy). We do this by considering local random unitary circuits with a conserved U$(1)$ charge and argue, with numerical and analytical evidence, that the presence of a conservation law slows relaxation in both time ordered {\textit{and}} out-of-time-ordered correlation functions, both can have a diffusively relaxing component or "hydrodynamic tail" at late times. We verify the presence of such tails also in a deterministic, peridocially driven system. We show that for OTOCs, the combination of diffusive and ballistic components leads to a wave front with a specific, asymmetric shape, decaying as a power law behind the front. These results also explain existing numerical investigations in non-noisy ergodic systems with energy conservation. Moreover, we consider OTOCs in Gibbs states, parametrized by a chemical potential $\mu$, and apply perturbative arguments to show that for $\mu\gg 1$ the ballistic front of information-spreading can only develop at times exponentially large in $\mu$ -- with the information traveling diffusively at earlier times. We also develop a new formalism for describing OTOCs and operator spreading, which allows us to interpret the saturation of OTOCs as a form of thermalization on the Hilbert space of operators.

categories: cond-mat.stat-mech cond-mat.str-el hep-th quant-ph


title: Gaussian Process Regression for Generalized Frequency Response Function Estimation

abstract: Kernel-based modeling of dynamic systems has garnered a significant amount of attention in the system identification literature since its introduction to the field. While the method was originally applied to linear impulse response estimation in the time domain, the concepts have since been extended to the frequency domain for estimation of frequency response functions (FRFs), as well as to the estimation of the Volterra series in time domain. In the latter case, smoothness and exponential decay was imposed along the hypersurfaces of the multidimensional impulse responses, allowing lower variance estimates than could be obtained in a simple least squares framework. The Volterra series can also be expressed in a frequency domain context, however there are several competing representations which all possess some unique advantages. Perhaps the most natural representation is the generalized frequency response function (GFRF), which is defined as the multidimensional Fourier transform of the corresponding Volterra kernel in the time-domain series. The representation leads to a series of frequency domain functions with increasing dimension.


title: Dynamic Routing Between Capsules

abstract: A capsule is a group of neurons whose activity vector represents the instantiation parameters of a specific type of entity such as an object or an object part. We use the length of the activity vector to represent the probability that the entity exists and its orientation to represent the instantiation parameters. Active capsules at one level make predictions, via transformation matrices, for the instantiation parameters of higher-level capsules. When multiple predictions agree, a higher level capsule becomes active. We show that a discrimininatively trained, multi-layer capsule system achieves state-of-the-art performance on MNIST and is considerably better than a convolutional net at recognizing highly overlapping digits. To achieve these results we use an iterative routing-by-agreement mechanism: A lower-level capsule prefers to send its output to higher level capsules whose activity vectors have a big scalar product with the prediction coming from the lower-level capsule.


title: Computations over Local Rings in Macaulay2

abstract: Local rings are ubiquitous in algebraic geometry. Not only are they naturally meaningful in a geometric sense, but also they are extremely useful as many problems can be attacked by first reducing to the local case and taking advantage of their nice properties. Any localization of a ring $R$, for instance, is flat over $R$. Similarly, when studying finitely generated modules over local rings, projectivity, flatness, and freeness are all equivalent. We introduce the packages PruneComplex, Localization and LocalRings for Macaulay2. The first package consists of methods for pruning chain complexes over polynomial rings and their localization at prime ideals. The second package contains the implementation of such local rings. Lastly, the third package implements various computations for local rings, including syzygies, minimal free resolutions, length, minimal generators and presentation, and the Hilbert--Samuel function. The main tools and procedures in this paper involve homological methods. In particular, many results depend on computing the minimal free resolution of modules over local rings.


title: A positivity preserving and conservative variational scheme for phase-field modeling of two-phase flows

abstract: We present a positivity preserving variational scheme for the phase-field modeling of incompressible two-phase flows with high density ratio and using meshes of arbitrary topology. The variational finite element technique relies on the Allen-Cahn phase-field equation for capturing the phase interface on a fixed mesh with a mass conservative and energy-stable discretization. Mass is conserved by enforcing a Lagrange multiplier which has both temporal and spatial dependence on the solution of the phase-field equation. The spatial part of the Lagrange multiplier is written as a mid-point approximation to make the scheme energy-stable. This enables us to form a conservative, energy-stable and positivity preserving scheme. The proposed variational technique reduces spurious and unphysical oscillations in the solution while maintaining second-order spatial accuracy. To model a generic two-phase free-surface flow, we couple the Allen-Cahn phase-field equation with the Navier-Stokes equations. Comparison of results between standard linear stabilized finite element method and the present variational formulation shows a remarkable reduction of oscillations in the solution while retaining the boundedness of the phase-indicator field. We perform a standalone test to verify the accuracy and stability of the Allen-Cahn two-phase solver. Standard two-phase flow benchmarks such as Laplace-Young law and sloshing tank problem are carried out to assess the convergence and accuracy of the coupled Navier-Stokes and Allen-Cahn solver. Two- and three-dimensional dam break problem are then solved to assess the scheme for the problem with topological changes of the air-water interface on unstructured meshes. Finally, we demonstrate the phase-field solver for a practical problem of wave-structure interaction in offshore engineering using general three-dimensional unstructured meshes.


title: Effective chemical potential for non-equilibrium systems and its application to molecular beam epitaxy of Bi2Se3

abstract: First-principles studies often rely on the assumption of equilibrium, which can be a poor approximation, e.g., for growth. Here, an effective chemical potential method for non-equilibrium systems is developed. A salient feature of the theory is that it maintains the equilibrium limits as the correct limit. In application to molecular beam epitaxy, rate equations are solved for the concentrations of small clusters, which serve as feedstock for growth. We find that the effective chemical potential is determined by the most probable, rather than by the lowest-energy, cluster. In the case of Bi2Se3, the chemical potential is found to be highly supersaturated, leading to a high nucleus concentration in agreement with experiment.


title: ADHM and the 4d Quantum Hall Effect

abstract: Yang-Mills instantons are solitonic particles in d=4+1 dimensional gauge theories. We construct and analyse the quantum Hall states that arise when these particles are restricted to the lowest Landau level. We describe the ground state wavefunctions for both Abelian and non-Abelian quantum Hall states. Although our model is purely bosonic, we show that the excitations of this 4d quantum Hall state are governed by the Nekrasov partition function of a certain five dimensional supersymmetric gauge theory with Chern-Simons term. The partition function can also be interpreted as a variant of the Hilbert series of the instanton moduli space, counting holomorphic sections rather than holomorphic functions. It is known that the Hilbert series of the instanton moduli space can be rewritten using mirror symmetry of 3d gauge theories in terms of Coulomb branch variables. We generalise this approach to include the effect of a five dimensional Chern-Simons term. We demonstrate that the resulting Coulomb branch formula coincides with the corresponding Higgs branch Molien integral which, in turn, reproduces the standard formula for the Nekrasov partition function.


title: Deep Illumination: Approximating Dynamic Global Illumination with Generative Adversarial Network

abstract: We present Deep Illumination, a novel machine learning technique for approximating global illumination (GI) in real-time applications using a Conditional Generative Adversarial Network. Our primary focus is on generating indirect illumination and soft shadows with offline rendering quality at interactive rates. Inspired from recent advancement in image-to-image translation problems using deep generative convolutional networks, we introduce a variant of this network that learns a mapping from Gbuffers (depth map, normal map, and diffuse map) and direct illumination to any global illumination solution. Our primary contribution is showing that a generative model can be used to learn a density estimation from screen space buffers to an advanced illumination model for a 3D environment. Once trained, our network can approximate global illumination for scene configurations it has never encountered before within the environment it was trained on. We evaluate Deep Illumination through a comparison with both a state of the art real-time GI technique (VXGI) and an offline rendering GI technique (path tracing). We show that our method produces effective GI approximations and is also computationally cheaper than existing GI techniques. Our technique has the potential to replace existing precomputed and screen-space techniques for producing global illumination effects in dynamic scenes with physically-based rendering quality.


title: Operator spreading and the emergence of dissipative hydrodynamics under unitary evolution with conservation laws

abstract: We study the scrambling of local quantum information in chaotic many-body systems in the presence of a locally conserved quantity like charge or energy that moves diffusively. The interplay between conservation laws and scrambling sheds light on the mechanism by which unitary quantum dynamics, which is reversible, gives rise to diffusive hydrodynamics, which is a dissipative process. We obtain our results in a random quantum circuit model that is constrained to have a conservation law. We find that a generic spreading operator consists of two parts: (i) a conserved part which comprises the weight of the spreading operator on the local conserved densities, whose dynamics is described by diffusive charge spreading. This conserved part also acts as a source that steadily emits a flux of (ii) non-conserved operators. This emission leads to dissipation in the operator hydrodynamics, with the dissipative process being the conversion of operator weight from local conserved operators to nonconserved, at a rate set by the local diffusion current. The emitted nonconserved parts then spread ballistically at a butterfly speed, thus becoming highly nonlocal and hence essentially non-observable, thereby acting as the "reservoir" that facilitates the dissipation. In addition, we find that the nonconserved component develops a power law tail behind its leading ballistic front due to the slow dynamics of the conserved components. This implies that the out-of-time-order commutator (OTOC) between two initially separated operators grows sharply upon the arrival of the ballistic front but, in contrast to systems with no conservation laws, it develops a diffusive tail and approaches its asymptotic late-time value only as a power of time instead of exponentially. We also derive these results within an effective hydrodynamic description which contains multiple coupled diffusion equations.


title: Estimation of Rain Attenuation at EHF bands for Earth-to-Satellite Links in Bangladesh

abstract: Due to heavy congestion in lower frequency bands, engineers are looking for new frequency bands to support new services that require higher data rates, which in turn needs broader bandwidths. To meet this requirement, extremely high frequency (EHF), particularly Q (36 to 46 GHz) and V (46 to 56 GHz) bands, is the best viable solution because of its complete availability. The most serious challenge the EHF band poses is the attenuation caused by rain. This paper investigates the effect of the rain on Q and V bands' performances in Bangladeshi climatic conditions. The rain attenuations of the two bands are predicted for the four main regions of Bangladesh using ITU rain attenuation model. The measured rain statistics is used for this prediction. It is observed that the attenuation due to rain in the Q/V band reaches up to 150 dB which is much higher than that of the currently used Ka band. The variability of the rain attenuation is also investigated over different sessions of Bangladesh. The attenuation varies from 40 dB to 170 dB depending on the months. Finally, the amount of rain fade required to compensate the high rain attenuation is also predicted for different elevation angles.


title: Saturation of the Quantum Null Energy Condition in Far-From-Equilibrium Systems

abstract: The Quantum Null Energy Condition (QNEC) is a new local energy condition that a general Quantum Field Theory (QFT) is believed to satisfy, relating the classical null energy condition (NEC) to the second functional derivative of the entanglement entropy in the corresponding null direction. We present the first series of explicit computations of QNEC in a strongly coupled QFT, using holography. We consider the vacuum, thermal equilibrium, a homogeneous far-from-equilibrium quench as well as a colliding system that violates NEC. For vacuum and the thermal phase QNEC is always weaker than NEC. While for the homogeneous quench QNEC is satisfied with a finite gap, we find the interesting result that the colliding system can saturate QNEC, depending on the null direction.


title: An elevation of 0.1 light-seconds for the optical jet base in an accreting Galactic black hole system

abstract: Relativistic plasma jets are observed in many accreting black holes. According to theory, coiled magnetic fields close to the black hole accelerate and collimate the plasma, leading to a jet being launched. Isolating emission from this acceleration and collimation zone is key to measuring its size and understanding jet formation physics. But this is challenging because emission from the jet base cannot be easily disentangled from other accreting components. Here, we show that rapid optical flux variations from a Galactic black-hole binary are delayed with respect to X-rays radiated from close to the black hole by ~0.1 seconds, and that this delayed signal appears together with a brightening radio jet. The origin of these sub-second optical variations has hitherto been controversial. Not only does our work strongly support a jet origin for the optical variations, it also sets a characteristic elevation of <~10$^3$ Schwarzschild radii for the main inner optical emission zone above the black hole, constraining both internal shock and magnetohydrodynamic models. Similarities with blazars suggest that jet structure and launching physics could potentially be unified under mass-invariant models. Two of the best-studied jetted black hole binaries show very similar optical lags, so this size scale may be a defining feature of such systems.


title: Detecting Black Hole Binaries by Gaia

abstract: We study the prospect of the Gaia satellite to identify black hole binary systems by detecting the orbital motion of the companion stars. Taking into account the initial mass function, mass transfer, common envelope phase, interstellar absorption and identifiability of black holes, we estimate the number of black hole binaries detected by Gaia and their distributions with respect to the black hole mass for several models with different parameters. We find that $\sim 200-1000$ black hole binaries will be detected by Gaia during its 5 years operation. The shape of distribution function of the black hole mass is affected most severely by the relation between the zero-age main sequence (ZAMS) stellar mass and the black hole mass in parameters adopted in this paper, which implies that black holes detected with Gaia enables us to constrain the mass relation.


title: Fast preparation of critical ground states using superluminal fronts

abstract: We propose a spatio-temporal quench protocol that allows for the fast preparation of ground states of gapless models with Lorentz invariance. Assuming the system initially resides in the ground state of a corresponding massive model, we show that a superluminally-moving `front' that $\textit{locally}$ quenches the mass, leaves behind it (in space) a state $\textit{arbitrarily close}$ to the ground state of the gapless model. Importantly, our protocol takes time $\mathcal{O} \left( L \right)$ to produce the ground state of a system of size $\sim L^d$ ($d$ spatial dimensions), while a fully adiabatic protocol requires time $\sim \mathcal{O} \left( L^2 \right)$ to produce a state with exponential accuracy in $L$. The physics of the dynamical problem can be understood in terms of relativistic rarefaction of excitations generated by the mass front. We provide proof-of-concept by solving the proposed quench exactly for a system of free bosons in arbitrary dimensions, and for free fermions in $d = 1$. We discuss the role of interactions and UV effects on the free-theory idealization, before numerically illustrating the usefulness of the approach via simulations on the quantum Heisenberg spin-chain.


title: The impact of Gaia on our understanding of the Vast Polar Structure of the Milky Way

abstract: The Milky Way (MW) is surrounded by a Vast Polar Structure (VPOS) of satellite galaxies, star clusters, and streams. Proper motion measurements for the brightest MW satellites indicate that they are predominantly co-orbiting along the VPOS. This is consistent with a dynamically stable structure. Assuming that all satellites that are aligned with the VPOS also co-orbit along this structure allows to empirically predict their systemic proper motions. Testing predictions for individual satellite galaxies at large distances requires high-accuracy proper motion measurements such as with the Hubble Space Telescope. However, for nearby MW satellites, Gaia will allow to test these predictions, in particular for a statistical sample of satellites since proper motion predictions exist for almost all of them. This will clarify how rotationally supported the VPOS is. In addition, Gaia will discover Galactic substructure, in particular stellar streams. The degree of their alignment with the VPOS might further constrain its stability and nature.


title: Machine learning vortices at the Kosterlitz-Thouless transition

abstract: Efficient and automated classification of phases from minimally processed data is one goal of machine learning in condensed matter and statistical physics. Supervised algorithms trained on raw samples of microstates can successfully detect conventional phase transitions via learning a bulk feature such as an order parameter. In this paper, we investigate whether neural networks can learn to classify phases based on topological defects. We address this question on the two-dimensional classical XY model which exhibits a Kosterlitz-Thouless transition. We find significant feature engineering of the raw spin states is required to convincingly claim that features of the vortex configurations are responsible for learning the transition temperature. We further show a single-layer network does not correctly classify the phases of the XY model, while a convolutional network easily performs classification by learning the global magnetization. Finally, we design a deep network capable of learning vortices without feature engineering. We demonstrate the detection of vortices does not necessarily result in the best classification accuracy, especially for lattices of less than approximately 1000 spins. For larger systems, it remains a difficult task to learn vortices.


title: Absence of Thermalization in Finite Isolated Interacting Floquet Systems

abstract: Conventional wisdom suggests that the long time behavior of isolated interacting periodically driven (Floquet) systems is a featureless maximal entropy state characterized by an infinite temperature. Efforts to thwart this uninteresting fixed point include adding sufficient disorder to realize a Floquet many-body localized phase or working in a narrow region of drive frequencies to achieve glassy non-thermal behavior at long time. Here we show that in clean systems the Floquet eigenstates can exhibit non-thermal behavior due to finite system size. We consider a one-dimensional system of spinless fermions with nearest-neighbor interactions where the interaction term is driven. Interestingly, even with no static component of the interaction, the quasienergy spectrum contains gaps and a significant fraction of the Floquet eigenstates, at all quasienergies, have non-thermal average doublon densities. We show that this non-thermal behavior arises due to emergent integrability at large interaction strength and discuss how the integrability breaks down with power-law dependence on system size.


title: Alone Together: Compositional Reasoning and Inference for Weak Isolation

abstract: Serializability is a well-understood correctness criterion that simplifies reasoning about the behavior of concurrent transactions by ensuring they are isolated from each other while they execute. However, enforcing serializable isolation comes at a steep cost in performance and hence database systems in practice support, and often encourage, developers to implement transactions using weaker alternatives. Unfortunately, the semantics of weak isolation is poorly understood, and usually explained only informally in terms of low-level implementation artifacts. Consequently, verifying high-level correctness properties in such environments remains a challenging problem. To address this issue, we present a novel program logic that enables compositional reasoning about the behavior of concurrently executing weakly-isolated transactions. Recognizing that the proof burden necessary to use this logic may dissuade application developers, we also describe an inference procedure based on this foundation that ascertains the weakest isolation level that still guarantees the safety of high-level consistency invariants associated with such transactions. The key to effective inference is the observation that weakly-isolated transactions can be viewed as functional (monadic) computations over an abstract database state, allowing us to treat their operations as state transformers over the database. This interpretation enables automated verification using off-the-shelf SMT solvers. Case studies and experiments of real-world applications (written in an embedded DSL in OCaml) demonstrate the utility of our approach.


title: Pre-inflationary dynamics in loop quantum cosmology: Power-law potentials

abstract: In this paper, we study the pre-inflationary dynamics for the power-law potential $(V(\phi) \propto \phi^n)$ with $n<2$ in the framework of loop quantum cosmology. In the case where the kinetic energy of the inflaton dominates at the initial, the evolution of the universe can always be divided into three different phases prior to preheating: {\em bouncing, transition and slow-roll inflation}. During the bouncing phase, the evolution of the expansion factor is independent not only on the initial conditions but also the inflationary potentials, and is given explicitly by an analytical solution. In contrast, for the potential energy dominated initial conditions, this universality is lost. We also obtain total number of e-folds during the slow-roll inflation, whereby physically viable models are identified. In addition, we present phase space analysis for the inflationary potentials under consideration and compare our results with the ones obtained previously for different potentials.


title: Cosmological singlet diagnostics of neutrinophilic dark matter

abstract: The standard model of particle physics is extended by adding a purely neutrinophilic dark sector. It is shown that theories which accommodate standard model neutrinos as dark radiation are resurrected. Sterile neutrinos bridge the visible and the dark sector and can keep their mutual interactions effective even after the epoch of big bang nucleosythesis. This prolonged contact between dark matter and standard model neutrinos solves the well-known small-scale structure problems of the cosmological standard model. The dark sector presented in this work satisfies all experimental and observational constraints from particle physics and cosmology.


title: Confirming chemical clocks: asteroseismic age dissection of the Milky Way disk(s)

abstract: Investigations of the origin and evolution of the Milky Way disk have long relied on chemical and kinematic identification of its components to reconstruct our Galactic past. Difficulties in determining precise stellar ages have restricted most studies to small samples, normally confined to the solar neighbourhood. Here we break this impasse with the help of asteroseismic inference and perform a chronology of the evolution of the disk throughout the age of the Galaxy. We chemically dissect the Milky Way disk population using a sample of red giant stars spanning out to 2~kpc in the solar annulus observed by the {\it Kepler} satellite, with the added dimension of asteroseismic ages. Our results reveal a clear difference in age between the low- and high-$\alpha$ populations, which also show distinct velocity dispersions in the $V$ and $W$ components. We find no tight correlation between age and metallicity nor [$\alpha$/Fe] for the high-$\alpha$ disk stars. Our results indicate that this component formed over a period of more than 2~Gyr with a wide range of [M/H] and [$\alpha$/Fe] independent of time. Our findings show that the kinematic properties of young $\alpha$-rich stars are consistent with the rest of the high-$\alpha$ population and different from the low-$\alpha$ stars of similar age, rendering support to their origin being old stars that went through a mass transfer or stellar merger event, making them appear younger, instead of migration of truly young stars formed close to the Galactic bar.


title: Implementation and characterization of a controllable dephasing channel based on coupling polarization and spatial degrees of freedom of light

abstract: We present the experimental implementation and theoretical model of a controllable dephasing quantum channel using photonic systems. The channel is implemented by coupling the polarization and the spatial distribution of light, that play, in the perspective of Open Quantum Systems, the role of quantum system and environment, respectively. The capability of controlling our channel allows us to visualize its effects in a quantum system. Differently from standard dephasing channels, our channel presents an exotic behavior in the sense that the evolution of a state, from a pure to a mixed state, shows an oscillatory behavior if tracked in the Bloch sphere. Additionally, we report the evolution of the purity and perform a quantum process tomography to obtain the \chi matrix associated to our channel.


title: Resolving ambiguities in model-independent partial-wave analysis of three-body decay

abstract: Partial wave analysis is an important tool for analyzing large data sets in hadronic decays of light and heavy mesons. It commonly relies on the isobar model, which assumes multihadron final states originate from successive two-body decays of well-known undisturbed intermediate states. Recently, analyses of heavy-meson decays and diffractively produced states have attempted to overcome the strong model dependences of the isobar model. These analyses have overlooked that model-independent, or freed-isobar, partial-wave analysis can break the orthogonality of partial waves and introduce mathematical ambiguities into results. We show how these ambiguities arise and present general techniques for identifying their presence and for correcting for them. We demonstrate these techniques with specific examples in both heavy-meson decay and pion--proton scattering.


title: A Note on Lower Bounds for Induced Ramsey Numbers

abstract: We say that a graph $F$ strongly arrows a pair of graphs $(G,H)$ if any 2-colouring of its edges with red and blue leads to either a red $G$ or a blue $H$ appearing as induced subgraphs of $F$. The induced Ramsey number, $IR(G,H)$ is defined as the minimum number of vertices of a graph $F$ which strongly arrows a pair $(G,H)$. We will consider two aspects of induced Ramsey numbers. Firstly there will be shown that the lower bound of the induced Ramsey number for a connected graph $G$ with independence number $\alpha$ and a graph $H$ with clique number $\omega$ roughly $\frac{\omega^2\alpha}{2}$. This bounds is sharp. Moreover we discuss also the case when $G$ is not connected providing also a sharp lower bound which is linear in both parameters


title: Inefficient jet-induced star formation in Centaurus A: High resolution ALMA observations of the northern filaments

abstract: NGC 5128 is one of the best targets to study AGN-feedback in the local Universe. At 13.5 kpc from the galaxy, optical filaments with recent star formation lie along the radio-jet direction. It is a testbed region for positive feedback (jet-induced star formation). APEX revealed strong CO emission in star-forming regions but also in regions with no detected tracers of star formation. When observed, star formation appears to be inefficient compared to the Kennicutt-Schmidt relation. We used ALMA to map the 12CO(1-0) emission all along the filaments at 1.3"~ 23.8 pc resolution. The CO emission is clumpy and distributed in two main structures: (i) the Horseshoe complex, outside the HI cloud, where gas is mostly excited by shocks and no star formation is observed; (ii) the Vertical filament, at the edge of the HI shell, which is a region of moderate star formation. We identified 140 molecular clouds. A statistical study reveals that they have very similar physical properties that in the inner Milky Way. However, the range of radius available with the present observations does not enable to investigate whether the clouds follow the Larson relation or not. The large virial parameter of the clouds suggests that gravity is not dominant. Finally, the total energy injection in the filaments is of the same order as in the inner part of the Milky Way. The strong CO emission detected in the filaments is an indication that the energy injected by the jet acts positively in the formation of dense molecular gas. The relatively high virial parameter of the molecular clouds suggests that the injected kinetic energy is too strong for star formation to be efficient. This is particularly the case in the Horseshoe complex where the virial parameter is the largest and where strong CO is detected with no associated star formation. This is the first evidence of inefficient AGN positive feedback.


title: BayesVP: a Bayesian Voigt profile fitting package

abstract: We introduce a Bayesian approach for modeling Voigt profiles in absorption spectroscopy and its implementation in the python package, BayesVP, publicly available at https://github.com/cameronliang/BayesVP. The code fits the absorption line profiles within specified wavelength ranges and generates posterior distributions for the column density, Doppler parameter, and redshifts of the corresponding absorbers. The code uses publicly available efficient parallel sampling packages to sample posterior and thus can be run on parallel platforms. BayesVP supports simultaneous fitting for multiple absorption components in high-dimensional parameter space. We provide other useful utilities in the package, such as explicit specification of priors of model parameters, continuum model, Bayesian model comparison criteria, and posterior sampling convergence check.


title: Characterization of Invariant subspaces in the polydisc

abstract: We give a complete characterization of invariant subspaces for $(M_{z_1}, \ldots, M_{z_n})$ on the Hardy space $H^2(\mathbb{D}^n)$ over the unit polydisc $\mathbb{D}^n$ in $\mathbb{C}^n$, $n >1$. In particular, this yields a complete set of unitary invariants for invariant subspaces for $(M_{z_1}, \ldots, M_{z_n})$ on $H^2(\mathbb{D}^n)$, $n > 1$. As a consequence, we classify a large class of $n$-tuples, $n > 1$, of commuting isometries. All of our results hold for vector-valued Hardy spaces over $\mathbb{D}^n$, $n > 1$. Our invariant subspace theorem solves the well-known open problem on characterizations of invariant subspaces of the Hardy space over the unit polydisc.


title: Gradient Sparsification for Communication-Efficient Distributed Optimization

abstract: Modern large scale machine learning applications require stochastic optimization algorithms to be implemented on distributed computational architectures. A key bottleneck is the communication overhead for exchanging information such as stochastic gradients among different workers. In this paper, to reduce the communication cost we propose a convex optimization formulation to minimize the coding length of stochastic gradients. To solve the optimal sparsification efficiently, several simple and fast algorithms are proposed for approximate solution, with theoretical guaranteed for sparseness. Experiments on $\ell_2$ regularized logistic regression, support vector machines, and convolutional neural networks validate our sparsification approaches.


title: Enhanced laser-driven ion acceleration by superponderomotive electrons generated from near-critical-density plasma

abstract: We report on the experimental studies of laser driven ion acceleration from double-layer target where a near-critical density target with a few-micron thickness is coated in front of a nanometer thin diamond-like carbon foil. A significant enhancement of proton maximum energies from 12 to ~30 MeV is observed when relativistic laser pulse impinge on the double-layer target under linear polarization. We attributed the enhanced acceleration to superponderomotive electrons that were simultaneously measured in the experiments with energies far beyond the free-electron ponderomotive limit. Our interpretation is supported by two-dimensional simulation results.


title: Software Engineering Modeling Applied to English Verb Classification (and Poetry)

abstract: In requirements specification, software engineers create a textual description of the envisioned system as well as develop conceptual models using such tools as Universal Modeling Language (UML) and System Modeling Language (SysML). One such tool, called FM, has recently been developed as an extension of the INPUT-PROCESS-OUTPUT (IPO) model. IPO has been used extensively in many interdisciplinary applications and is described as one of the most fundamental and important of all descriptive tools. This paper is an attempt to understanding the PROCESS in IPO. The fundamental way to describe PROCESS is in verbs. This use of language has an important implication for systems modeling since verbs express the vast range of actions and movements of all things. It is clear that modeling needs to examine verbs. Accordingly, this paper involves a study of English verbs as a bridge to learn about processes, not as linguistic analysis but rather to reveal the semantics of processes, particularly the five verbs that form the basis of FM states: create, process, receive, release, and transfer. The paper focuses on verb classification, and specifically on how to model the action of verbs diagrammatically. From the linguistics point of view, according to some researchers, further exploration of the notion of verb classes is needed for real-world tasks such as machine translation, language generation, and document classification. Accordingly, this non-linguistics study may benefit linguistics.


title: Spectral invariants and playing hide-and-seek on surfaces

abstract: We prove the expected duration of a game of hide-and-seek played on a Riemannian manifold under the laws of Brownian Motion is a spectral invariant: it is a zeta-regularized version of the `trace' of the Laplacian. An analogous hide-and-seek game may be played on Markov chains, where the spectral invariant that emerges is a classical quantity known as Kemeny's constant. We develop the analogies between the two settings in order to highlight the connections between the regularized trace and Kemeny's constant. Our proof relies on the connections between Green's functions and expected hitting times, and the fact that the regularized trace may be approached via the Green's function.


title: First Monte Carlo global analysis of nucleon transversity with lattice QCD constraints

abstract: We report on the first global QCD analysis of the quark transversity distributions in the nucleon from semi-inclusive deep-inelastic scattering (SIDIS), using a new Monte Carlo method based on nested sampling and constraints on the isovector tensor charge $g_T$ from lattice QCD. A simultaneous fit to the available SIDIS Collins asymmetry data is compatible with $g_T$ values extracted from a comprehensive reanalysis of existing lattice simulations, in contrast to previous analyses which found significantly smaller $g_T$ values. The contributions to the nucleon tensor charge from $u$ and $d$ quarks are found to be $\delta u = 0.3(2)$ and $\delta d = -0.7(2)$ at a scale $Q^2 = 2$ GeV$^2$.


title: Kernel k-Groups via Hartigan's Method

abstract: Energy statistics was proposed by Sz\' ekely in the 80's inspired by Newton's gravitational potential in classical mechanics and it provides a model-free hypothesis test for equality of distributions. In its original form, energy statistics was formulated in Euclidean spaces. More recently, it was generalized to metric spaces of negative type. In this paper, we consider a formulation for the clustering problem using a weighted version of energy statistics in spaces of negative type. We show that this approach leads to a quadratically constrained quadratic program in the associated kernel space, establishing connections with graph partitioning problems and kernel methods in machine learning. To find local solutions of such an optimization problem, we propose kernel k-groups, which is an extension of Hartigan's method to kernel spaces. Kernel k-groups is cheaper than spectral clustering and has the same computational cost as kernel k-means (which is based on Lloyd's heuristic) but our numerical results show an improved performance, especially in higher dimensions. Moreover, we verify the efficiency of kernel k-groups in community detection in sparse stochastic block models which has fascinating applications in several areas of science.

categories: stat.ML cs.CV cs.DS cs.LG math.ST stat.TH


title: DoShiCo Challenge: Domain Shift in Control Prediction

abstract: Training deep neural network policies end-to-end for real-world applications so far requires big demonstration datasets in the real world or big sets consisting of a large variety of realistic and closely related 3D CAD models. These real or virtual data should, moreover, have very similar characteristics to the conditions expected at test time. These stringent requirements and the time consuming data collection processes that they entail, are currently the most important impediment that keeps deep reinforcement learning from being deployed in real-world applications. Therefore, in this work we advocate an alternative approach, where instead of avoiding any domain shift by carefully selecting the training data, the goal is to learn a policy that can cope with it. To this end, we propose the DoShiCo challenge: to train a model in very basic synthetic environments, far from realistic, in a way that it can be applied in more realistic environments as well as take the control decisions on real-world data. In particular, we focus on the task of collision avoidance for drones. We created a set of simulated environments that can be used as benchmark and implemented a baseline method, exploiting depth prediction as an auxiliary task to help overcome the domain shift. Even though the policy is trained in very basic environments, it can learn to fly without collisions in a very different realistic simulated environment. Of course several benchmarks for reinforcement learning already exist - but they never include a large domain shift. On the other hand, several benchmarks in computer vision focus on the domain shift, but they take the form of a static datasets instead of simulated environments. In this work we claim that it is crucial to take the two challenges together in one benchmark.


title: Machine learning for crystal identification and discovery

abstract: As computers get faster, researchers -- not hardware or algorithms -- become the bottleneck in scientific discovery. Computational study of colloidal self-assembly is one area that is keenly affected: even after computers generate massive amounts of raw data, performing an exhaustive search to determine what (if any) ordered structures occur in a large parameter space of many simulations can be excruciating. We demonstrate how machine learning can be applied to discover interesting areas of parameter space in colloidal self assembly. We create numerical fingerprints -- inspired by bond orientational order diagrams -- of structures found in self-assembly studies and use these descriptors to both find interesting regions in a phase diagram and identify characteristic local environments in simulations in an automated manner for simple and complex crystal structures. Utilizing these methods allows analysis methods to keep up with the data generation ability of modern high-throughput computing environments.


title: Description of pp forward elastic scattering at 7 and 8 TeV

abstract: We analyse the recent LHC data at 7 and 8 TeV for pp elastic scattering with special attention for the structure of the real part, which is shown to be crucial to describe the differential cross section in the forward region. We determine accurately the position of the zero of the real amplitude, which corresponds to the zero of a theorem by A. Martin.


title: Growth and structural characterization of large superconducting crystals of La$_{2-x}$Ca$_{1+x}$Cu$_2$O$_{6}$

abstract: Large crystals of La$_{2-x}$Ca$_{1+x}$Cu$_2$O$_{6}$ (La-Ca-2126) with $x=0.10$ and 0.15 have been grown and converted to bulk superconductors by high-pressure oxygen annealing. The superconducting transition temperature, $T_c$, is as high as 55~K; this can be raised to 60~K by post-annealing in air. Here we present structural and magnetic characterizations of these crystals using neutron scattering and muon spin rotation techniques. While the as-grown, non-superconducting crystals are single phase, we find that the superconducting crystals contain 3 phases forming coherent domains stacked along the $c$ axis: the dominant La-Ca-2126 phase, very thin (1.5 unit-cell) intergrowths of La$_2$CuO$_4$, and an antiferromagnetic La$_8$Cu$_8$O$_{20}$ phase. We propose that the formation and segregation of the latter phases increases the Ca concentration of the La-Ca-2126, thus providing the hole-doping that supports superconductivity.


title: Recursive functions and existentially closed structures

abstract: The purpose of this paper is to clarify the relationship between various conditions implying essential undecidability: our main result is that there exists a theory $T$ in which all partially recursive functions are representable, yet $T$ does not interpret Robinson's theory $R$. To this end, we borrow tools from model theory--specifically, we investigate model-theoretic properties of the model completion of the empty theory in a language with function symbols. We obtain a certain characterization of $\exists\forall$ theories interpretable in existential theories in the process.


title: Blowing bubbles on the torus

abstract: We consider the regularized trace of the inverse of the Laplacian on a skinny torus. With its flat metric, a skinny torus has large trace, but we show that there are conformally equivalent metrics making the trace close to that of a sphere of the same area. This behavior is in sharp contrast to that of the log-determinant, a well-known spectral invariant which is extremized at the flat metric on any torus. Our examples are bubbled tori, where you take a sphere, discard polar regions, and glue top to bottom. In a addendum, we belatedly notice that our bubbled tori have trace less than the sphere, and outline how to exploit this to get Okikiolu's result that by means of a conformal factor depending only on longitude, any torus can be made to have trace less than the sphere.


title: Aerosols and Methane in the Ice Giant Atmospheres Inferred from Spatially Resolved, Near-Infrared Spectra: I. Uranus, 2001-2007

abstract: We present a radiative transfer analysis of latitudinally resolved H (1.487-1.783 micron) and K (2.028-2.364 micron) band spectra of Uranus, from which we infer the distributions of aerosols and methane in the planet's atmosphere. Data were acquired in 2001, 2002, 2004, 2005, and 2007 using the 200-inch (5.1 m) Hale Telescope and the Palomar High Angular Resolution Observer (PHARO) near-infrared adaptive optics (AO) camera system (Hayward, 2001). Observations sample a range of latitudes between 80-deg S and 60-deg N on the Uranian disk. At each latitude, a vertical distributions of aerosols was retrieved using a custom non-linear constrained retrieval algorithm. Two layers of aerosols are needed to match the observations: a thin upper layer peaking just below the 100-mb tropopause and a lower clouds at ~1.9 bars. Latitudinal variations in aerosols are interpreted in context of notional circulation models, while temporal changes suggest potential seasonal effects. We infer significant reduction in aerosol scattering optical thickness in southern latitudes between 2001 and 2007, in agreement with trends reported in studies covering part of the same period using different data and retrieval algorithms (e.g., Irwin et al., 2009, 2010, 2012; Sromovsky et al., 2009). Best fits to the data are consistent with proposed models of polar depletion of methane (e.g., Karkoschka and Tomasko, 2011). Finally, a discrete cloud from 2007 is analyzed in context of simple parcel theory, with the goal of identifying likely formation mechanisms. The low scattering optical thicknesses of the discrete high cloud are consistent with formation associated with vortices and shallow lift rather than deep convection.


title: Understanding Early Word Learning in Situated Artificial Agents

abstract: Neural network-based systems can now learn to locate the referents of words and phrases in images, answer questions about visual scenes, and execute symbolic instructions as first-person actors in partially-observable worlds. To achieve this so-called grounded language learning, models must overcome challenges that infants face when learning their first words. While it is notable that models with no meaningful prior knowledge overcome these obstacles, researchers currently lack a clear understanding of how they do so, a problem that we attempt to address in this paper. For maximum control and generality, we focus on a simple neural network-based language learning agent, trained via policy-gradient methods, which can interpret single-word instructions in a simulated 3D world. Whilst the goal is not to explicitly model infant word learning, we take inspiration from experimental paradigms in developmental psychology and apply some of these to the artificial agent, exploring the conditions under which established human biases and learning effects emerge. We further propose a novel method for visualising semantic representations in the agent.


title: How far did we get in face spoofing detection?

abstract: The growing use of control access systems based on face recognition shed light over the need for even more accurate systems to detect face spoofing attacks. In this paper, an extensive analysis on face spoofing detection works published in the last decade is presented. The analyzed works are categorized by their fundamental parts, i.e., descriptors and classifiers. This structured survey also brings the temporal evolution of the face spoofing detection field, as well as a comparative analysis of the works considering the most important public data sets in the field. The methodology followed in this work is particularly relevant to observe trends in the existing approaches, to discuss still opened issues, and to propose new perspectives for the future of face spoofing detection.


title: Bounds for traces of Hecke operators and applications to modular and elliptic curves over a finite field

abstract: We give an upper bound for the trace of a Hecke operator acting on the space of holomorphic cusp forms with respect to certain congruence subgroups. Such an estimate has applications to the analytic theory of elliptic curves over a finite field, going beyond the Riemann hypothesis over finite fields. As the main tool to prove our bound on traces of Hecke operators, we develop a Petersson formula for newforms for general nebentype characters.


title: Exoplanet Atmosphere Retrieval from Multifractal Analysis of Secondary Eclipse Spectra

abstract: We extend a data-based model-free multifractal method of exoplanet detection to probe exoplanetary atmospheres. Whereas the transmission spectrum is studied during the primary eclipse, we analyze the emission spectrum during the secondary eclipse, thereby probing the atmospheric limb. In addition to the spectral structure of exoplanet atmospheres, the approach provides information to study phenomena such as atmospheric flows, tidal-locking behavior, and the dayside-nightside redistribution of energy. The approach is demonstrated using Spitzer data for exoplanet HD189733b. The central advantage of the method is the lack of model assumptions in the detection and observational schemes.

categories: astro-ph.EP astro-ph.IM physics.data-an stat.AP


title: Trajectory Deformations from Physical Human-Robot Interaction

abstract: Robots are finding new applications where physical interaction with a human is necessary: manufacturing, healthcare, and social tasks. Accordingly, the field of physical human-robot interaction (pHRI) has leveraged impedance control approaches, which support compliant interactions between human and robot. However, a limitation of traditional impedance control is that---despite provisions for the human to modify the robot's current trajectory---the human cannot affect the robot's future desired trajectory through pHRI. In this paper, we present an algorithm for physically interactive trajectory deformations which, when combined with impedance control, allows the human to modulate both the actual and desired trajectories of the robot. Unlike related works, our method explicitly deforms the future desired trajectory based on forces applied during pHRI, but does not require constant human guidance. We present our approach and verify that this method is compatible with traditional impedance control. Next, we use constrained optimization to derive the deformation shape. Finally, we describe an algorithm for real time implementation, and perform simulations to test the arbitration parameters. Experimental results demonstrate reduction in the human's effort and improvement in the movement quality when compared to pHRI with impedance control alone.


title: Addendum: New approach to the resummation of logarithms in Higgs-boson decays to a vector quarkonium plus a photon [Phys. Rev. D 95, 054018 (2017)]

abstract: In this addendum to Phys.\ Rev.\ D {\bf 95}, 054018 (2017) we recompute the rates for the decays of the Higgs boson to a vector quarkonium plus a photon, where the vector quarkonium is $J/\psi$, $\Upsilon(1S)$, $\Upsilon(2S)$, or $\Upsilon(3S)$. We correct an error in the Abel-Pad\'e summation formula that was used to carry out the evolution of the quarkonium light-cone distribution amplitude in Phys.\ Rev.\ D {\bf 95}, 054018 (2017). We also correct an error in the scale of quarkonium wave function at the origin in Phys.\ Rev.\ D {\bf 95}, 054018 (2017) and introduce several additional refinements in the calculation.


title: Joint distribution in residue classes of the base-$q$ and Ostrowski digital sums

abstract: Let $q$ be an integer $\geq 2$ and let $S_q(n)$ denote the sum of digits of $n$ in base $q$. For \[ \alpha=[0;\overline{1,m}],\ m\geq 2, \] let $S_{\alpha}(n)$ denote the sum of digits in the Ostrowski $\alpha$-representation of $n$. Let $m_1,m_2\geq 2$ be integers with $$\gcd(q-1,m_1)=\gcd(m,m_2)=1.$$ We prove that there exists $\delta>0$ such that for all integers $a_1,a_2$, \begin{eqnarray*} &&|\{0\leq n<N: S_{q}(n)\equiv a_1\pmod{m_1},\ S_{\alpha}(n)\equiv a_2\pmod{m_2}\}| &=&\frac{N}{m_1m_2}+O(N^{1-\delta}). \end{eqnarray*} The asymptotic relation implied by this equality was proved by Coquet, Rhin & Toffin and the equality was proved for the case $\alpha=[\ \overline{1}\ ]$ by Spiegelhofer.


title: Loss Compensation in Time-Dependent Elastic Metamaterials

abstract: Materials with properties that are modulated in time are known to display wave phenomena showing energy increasing with time, with the rate mediated by the modulation. Until now there has been no accounting for material dissipation, which clearly counteracts energy growth. This paper provides an exact expression for the amplitude of elastic or acoustic waves propagating in lossy materials with properties that are periodically modulated in time. It is found that these materials can support a special propagation regime in which waves travel at constant amplitude, with temporal modulation compensating for the normal energy dissipation. We derive a general condition under which amplification due to time-dependent properties offsets the material dissipation. This identity relates band-gap properties associated with the temporal modulation and the average of the viscosity coefficient, thereby providing a simple recipe for the design of loss-compensated mechanical metamaterials.