title: $\mathcal W^{\alpha, p}$ and $C^{0,\gamma}$ regularity of solutions to $(\mu - \Delta + b \cdot \nabla)u=f$ with form-bounded vector fields

abstract: We consider the operator $-\Delta +b \cdot \nabla$ with $b:\mathbb R^d \rightarrow \mathbb R^d$ ($d \geq 3$) in the class of form-bounded vector fields (containing vector fields having critical-order singularities), and characterize quantitative dependence of the $\mathcal W^{1+\frac{2}{q},p}$ ($2 \leq p < q$) and the $C^{0,\gamma}$ regularity of solutions to the corresponding elliptic equation in $L^p$ on the value of the form-bound of $b$.


title: Spectroscopic fingerprints of many-body renormalizations in$1T$-TiSe$_2$

abstract: We have employed high resolution angle resolved photoemission spectroscopy (ARPES) measurements to investigate many-body renormalizations of the single-particle excitations in $1T$-TiSe$_2$. The energy distribution curves of the ARPES data reveal intrinsic peak-dip-hump feature, while the electronic dispersion derived from the momentum distribution curves of the data highlights, for the first time, multiple kink structures. These are canonical signatures of a coupling between the electronic degrees of freedom and some Bosonic mode in the system. We demonstrate this using a model calculation of the single-particle spectral function at the presence of an electron-Boson coupling. From the self-energy analysis of our ARPES data, we discern some of the critical energy scales of the involved Bosonic mode, which are $\sim$15 and 26 meV. Based on a comparison between these energies and the characteristic energy scales of our Raman scattering data, we identify these Bosonic modes as Raman active breathing (${\text{A}}_\text{1g}$) and shear (${\text{E}}_\text{g}$) modes, respectively. Direct observation of the band-renormalization due to electron-phonon coupling increases the possibility that electron-phonon interactions are central to the collective quantum states such as charge density wave (CDW) and superconductivity in the compounds based on $1T$-TiSe$_2$.


title: Multiwavelength Light Curves of Two Remarkable Sagittarius A* Flares

abstract: Sgr A*, the supermassive black hole (SMBH) at the center of our Milky Way Galaxy, is known to be a variable source of X-ray, near-infrared (NIR), and submillimeter (submm) radiation and therefore a prime candidate to study the electromagnetic radiation generated by mass accretion flow onto a black hole and/or a related jet. Disentangling the power source and emission mechanisms of this variability is a central challenge to our understanding of accretion flows around SMBHs. Simultaneous multiwavelength observations of the flux variations and their time correlations can play an important role in obtaining a better understanding of possible emission mechanisms and their origin. This paper presents observations of two flares that both apparently violate the previously established patterns in the relative timing of submm/NIR/X-ray flares from Sgr A*. One of these events provides the first evidence of coeval structure between NIR and submm flux increases, while the second event is the first example of the sequence of submm/X-ray/NIR flux increases all occurring within ~1 hr. Each of these two events appears to upend assumptions that have been the basis of some analytic models of flaring in Sgr A*. However, it cannot be ruled out that these events, even though unusual, were just coincidental. These observations demonstrate that we do not fully understand the origin of the multiwavelength variability of Sgr A*, and show that there is a continued and important need for long-term, coordinated, and precise multiwavelength observations of Sgr A* to characterize the full range of variability behavior.


title: Relative Amenability and Relative Soficity

abstract: We define a notion of relative soficity for countable groups with respect to a family of groups. A group is sofic if and only if it is relative sofic with respect to the family consisting only of the trivial group. If a group is relatively sofic with respect to a family of sofic groups, then the group is sofic. Using this notion we generalize a theorem of Elek and Szabo on the soficity of an extension of sofic groups by amenable groups. In particular we prove that groups that are relatively sofic with respect to a family of sofic groups are sofic and more importantly extensions of sofic groups by residually amenable groups are sofic. We also construct examples of relatively amenable groups with respect to an infinite family of subgroups but not with respect to any finite subfamily of the subgroups. As an application of these constructions we show that Deligne's central extension groups are sofic.


title: The Extended Solar Cycle: Muddying the Waters of Solar/Stellar Dynamo Modeling Or Providing Crucial Observational Constraints?

abstract: In 1844 Schwabe discovered that the number of sunspots increased and decreased over a period of about 11 years, that variation became known as the sunspot cycle. Almost eighty years later, Hale described the nature of the Sun's magnetic field, identifying that it takes about 22 years for the Sun's magnetic polarity to cycle. It was also identified that the latitudinal distribution of sunspots resembles the wings of a butterfly showing migration of sunspots in each hemisphere that abruptly start at mid-latitudes towards the Sun's equator over the next 11 years. These sunspot patterns were shown to be asymmetric across the equator. In intervening years, it was deduced that the Sun (and sun-like stars) possess magnetic activity cycles that are assumed to be the physical manifestation of a dynamo process that results from complex circulatory transport processes in the star's interior. Understanding the Sun's magnetism, its origin and its variation, has become a fundamental scientific objective \-- the distribution of magnetism, and its interaction with convective processes, drives various plasma processes in the outer atmosphere. In the past few decades, a range of diagnostic techniques have been employed to systematically study finer scale magnetized objects, and associated phenomena. The patterns discerned became known as the ``Extended Solar Cycle'' (ESC). The patterns of the ESC appeared to extend the wings of the activity butterfly back in time, nearly a decade before the formation of the sunspot pattern, and to much higher solar latitudes. In this short review, we describe their observational patterns of the ESC and discuss possible connections to the solar dynamo as we depart on a multi-national collaboration to investigate the origins of solar magnetism through a blend of archived and contemporary data analysis with the goal of improving solar dynamo understanding and modeling.


title: Quantum evolution of black hole initial data sets: Foundations

abstract: We construct a formalism for evolving spherically symmetric black hole initial data sets within a canonical approach to quantum gravity. This problem can be formulated precisely in quantum reduced loop gravity, a framework which has been successfully applied to give a full theory derivation of loop quantum cosmology. We extend this setting by implementing a particular choice of partial gauge which is then used to select a kinematical Hilbert space where the symmetry reduction is imposed through semiclassical states. The main result of this investigation is an effective Hamiltonian that can be used to solve for quantum black hole geometries by evolving classical black hole initial data sets.


title: Doubly Stochastic Adversarial Autoencoder

abstract: Any autoencoder network can be turned into a generative model by imposing an arbitrary prior distribution on its hidden code vector. Variational Autoencoder (VAE) [2] uses a KL divergence penalty to impose the prior, whereas Adversarial Autoencoder (AAE) [1] uses {\it generative adversarial networks} GAN [3]. GAN trades the complexities of {\it sampling} algorithms with the complexities of {\it searching} Nash equilibrium in minimax games. Such minimax architectures get trained with the help of data examples and gradients flowing through a generator and an adversary. A straightforward modification of AAE is to replace the adversary with the maximum mean discrepancy (MMD) test [4-5]. This replacement leads to a new type of probabilistic autoencoder, which is also discussed in our paper. We propose a novel probabilistic autoencoder in which the adversary of AAE is replaced with a space of {\it stochastic} functions. This replacement introduces a new source of randomness, which can be considered as a continuous control for encouraging {\it explorations}. This prevents the adversary from fitting too closely to the generator and therefore leads to a more diverse set of generated samples.


title: On the Mordell-Weil ranks of supersingular abelian varieties in cyclotomic extensions

abstract: Let $F$ be a number field unramified at an odd prime $p$ and $F_\infty$ be the $\mathbf{Z}_p$-cyclotomic extension of $F$. Let $A$ be an abelian variety defined over $F$ with good supersingular reduction at all primes of $F$ above $p$. B\"uy\"ukboduk and the first named author have defined modified Selmer groups associated to $A$ over $F_\infty$. Assuming that the Pontryagin dual of these Selmer groups are torsion $\mathbf{Z}_p[[\mathrm{Gal}(F_\infty/F)]]$-modules, we give an explicit sufficient condition for the rank of the Mordell-Weil group $A(F_n)$ to be bounded as $n$ varies.


title: Avoiding 3-Term Geometric Progressions in Hurwitz Quaternions

abstract: Several recent papers have considered the problem of how large a subset of integers can be without containing any 3-term geometric progressions. This problem has also recently been generalized to rings of integers in quadratic number fields and polynomial rings over finite fields. We study the analogous problem in the Hurwitz quaternion order to see how non-commutativity affects the problem. We compute an exact formula for the density of a 3-term geometric-progression-free set of Hurwitz quaternions arising from a greedy algorithm and derive upper and lower bounds for the supremum of upper densities of 3-term geometric-progression-free sets of Hurwitz quaternions.


title: Simple correction to bandgap in IV and III-V semiconductors: an improved first-principles local density functional theory

abstract: We report results from a fast, efficient, and first-principles full-potential N$^{th}$-order muffin-tin orbital (FP-NMTO) method combined with van Leeuwen-Baerends correction to local density exchange-correlation potential. We show that more complete and compact basis set is critical in improving the electronic and structural properties. We exemplify the self-consistent FP-NMTO calculations on group IV and III-V semiconductors. Notably, predicted bandgaps, lattice constants, and bulk moduli are in good agreement with experiments (e.g., we find for Ge $0.86~e$V, $5.57$~\AA, $75$~GPa vs. measured $0.74~e$V, $5.66$~\AA, $77.2$~GPa). We also showcase its application to the electronic properties of 2-dimensional $h-$BN and $h-$SiC, again finding good agreement with experiments.


title: On the structure of signed Selmer groups

abstract: Let $F$ be a number field unramified at an odd prime $p$ and $F_\infty$ be the $\mathbf{Z}_p$-cyclotomic extension of $F$. Generalizing Kobayashi plus/minus Selmer groups for elliptic curves, B\"uy\"ukboduk and Lei have defined modified Selmer groups, called signed Selmer groups, for certain non-ordinary $\mathrm{Gal}(\overline{F}/F)$-representations. In particular, their construction applies to abelian varieties defined over $F$ with good supersingular reduction at primes of $F$ dividing $p$. Assuming that these Selmer groups are cotorsion $\mathbf{Z}_p[[\mathrm{Gal}(F_\infty/F)]]$-modules, we show that they have no proper sub-$\mathbf{Z}_p[[\mathrm{Gal}(F_\infty/F)]]$-module of finite index. We deduce from this a number of arithmetic applications. On studying the Euler-Poincar\'e characteristic of these Selmer groups, we obtain an explicit formula on the size of the Bloch-Kato Selmer group attached to these representations. Furthermore, for two such representations that are isomorphic modulo $p$, we compare the Iwasawa-invariants of their signed Selmer groups.


title: Hydrodynamic modes in a magnetized chiral plasma with vorticity

abstract: By making use of a covariant formulation of the chiral kinetic theory in the relaxation-time approximation, we derive the first-order dissipative hydrodynamics equations for a charged chiral plasma with background electromagnetic fields. We identify the global equilibrium state for a rotating chiral plasma confined to a cylindrical region with realistic boundary conditions. Then, by using linearized hydrodynamic equations, supplemented by the Maxwell equations, we study hydrodynamic modes of magnetized rotating chiral plasma in the regimes of high temperature and high density. We find that, in both regimes, dynamical electromagnetism has profound effects on the spectrum of propagating modes. In particular, there are only the sound and Alfv\'en waves in the regime of high temperature, and the plasmons and helicons at high density. We also show that the chiral magnetic wave is universally overdamped because of high electrical conductivity in plasma that causes an efficient screening of charge fluctuations. The physics implications of the main results are briefly discussed.


title: Atmospheric Escape and the Evolution of Close-in Exoplanets

abstract: Exoplanets with substantial Hydrogen/Helium atmospheres have been discovered in abundance, many residing extremely close to their parent stars. The extreme irradiation levels these atmospheres experience causes them to undergo hydrodynamic atmospheric escape. Ongoing atmospheric escape has been observed to be occurring in a few nearby exoplanet systems through transit spectroscopy both for hot Jupiters and lower-mass super-Earths/mini-Neptunes. Detailed hydrodynamic calculations that incorporate radiative transfer and ionization chemistry are now common in one-dimensional models, and multi-dimensional calculations that incorporate magnetic-fields and interactions with the interstellar environment are cutting edge. However, there remains very limited comparison between simulations and observations. While hot Jupiters experience atmospheric escape, the mass-loss rates are not high enough to affect their evolution. However, for lower mass planets atmospheric escape drives and controls their evolution, sculpting the exoplanet population we observe today.


title: Unsupervised Metric Learning in Presence of Missing Data

abstract: For many machine learning tasks, the input data lie on a low-dimensional manifold embedded in a high dimensional space and, because of this high-dimensional structure, most algorithms are inefficient. The typical solution is to reduce the dimension of the input data using standard dimension reduction algorithms such as ISOMAP, LAPLACIAN EIGENMAPS or LLES. This approach, however, does not always work in practice as these algorithms require that we have somewhat ideal data. Unfortunately, most data sets either have missing entries or unacceptably noisy values. That is, real data are far from ideal and we cannot use these algorithms directly. In this paper, we focus on the case when we have missing data. Some techniques, such as matrix completion, can be used to fill in missing data but these methods do not capture the non-linear structure of the manifold. Here, we present a new algorithm MR-MISSING that extends these previous algorithms and can be used to compute low dimensional representation on data sets with missing entries. We demonstrate the effectiveness of our algorithm by running three different experiments. We visually verify the effectiveness of our algorithm on synthetic manifolds, we numerically compare our projections against those computed by first filling in data using nlPCA and mDRUR on the MNIST data set, and we also show that we can do classification on MNIST with missing data. We also provide a theoretical guarantee for MR-MISSING under some simplifying assumptions.


title: Multimodal Functions as Flow Signatures in Complex Porous Media

abstract: This study refutes the premise that the distribution of flow speeds in complex porous media can be described by a simple function such as a normal or exponential variation. In many complex porous media, including those relevant for subsurface storage and recovery applications, a separation of scales exists between larger inter-granular pores and micro-porosity inside the grains, leading to different flow signatures that need to be described by multimodal functions with distinct flow field characteristics. We demonstrate this finding by devising a novel methodology to simulate fluid flow in carbonate rock based on a representation of the pore space obtained by differential X-ray imaging. The permeability assignment in the micro-porous space is estimated from the pore size inferred from mercury injection porosimetry. Model predictions considering micro-porosity agree well with experimentally measured porosities and permeabilities. The micro-porosity can contribute significantly to the overall permeability, particularly in the more heterogeneous media.


title: Adaptive Variational Particle Filtering in Non-stationary Environments

abstract: Online convex optimization is a sequential prediction framework with the goal to track and adapt to the environment through evaluating proper convex loss functions. We study efficient particle filtering methods from the perspective of such a framework. We formulate an efficient particle filtering methods for the non-stationary environment by making connections with the online mirror descent algorithm which is known to be a universal online convex optimization algorithm. As a result of this connection, our proposed particle filtering algorithm proves to achieve optimal particle efficiency.


title: Derivation degree sequences of non-free arrangements

abstract: In this note we study the logarithmic derivation module of a non-free arrangement. We prove a generalized addition theorem for all arrangements. This addition theorem allows us to find various relationships between non-free arrangements, free arrangements and restriction counts. For graphic arrangements we can use these results to find a lower bound for the maximal degree generator in terms of triangles in the associated graph. We also apply these results to the case of hypersolvable arrangements where we define hyperexponents and use them to find a lower bound for their maximal degree generator.


title: A topological approach to Soergel theory

abstract: We develop a "Soergel theory" for Bruhat-constructible perverse sheaves on the flag variety $G/B$ of a complex reductive group $G$, with coefficients in an arbitrary field $\Bbbk$. Namely, we describe the endomorphisms of the projective cover of the skyscraper sheaf in terms of a "multiplicative" coinvariant algebra, and then establish an equivalence of categories between projective (or tilting) objects in this category and a certain category of "Soergel modules" over this algebra. We also obtain a description of the derived category of $T$-monodromic $\Bbbk$-sheaves on $G/U$ (where $U$, $T\subset B$ are the unipotent radical and the maximal torus), as a monoidal category, in terms of coherent sheaves on the formal neighborhood of the base point in $T^\vee_\Bbbk \times_{(T^\vee_\Bbbk)^W} T^\vee_\Bbbk$, where $T^\vee_\Bbbk$ is the $\Bbbk$-torus dual to $T$.


title: Sparse space-time models: Concentration Inequalities and Lasso

abstract: Inspired by Kalikow-type decompositions, we introduce a new stochastic model of infinite neuronal networks, for which we establish sharp oracle inequalities for Lasso methods and restricted eigenvalue properties for the associated Gram matrix with high probability. These results hold even if the network is only partially observed. The main argument rely on the fact that concentration inequalities can easily be derived whenever the transition probabilities of the underlying process admit a sparse space-time representation.


title: Dynamical generation of field mixing via flavor vacuum condensate

abstract: In this paper we study dynamical chiral symmetry breaking of a generic quantum field theoretical model with global chiral flavor symmetry. By purely algebraic means we analyze the vacuum structure for different symmetry breaking schemes and show explicitly how the ensuing non-trivial flavor vacuum condensate characterizes the phenomenon of field mixing. In addition, with the help of Ward-Takahashi identities we demonstrate the emergence of a correct number of Nambu-Goldstone modes in the physical spectrum.


title: SoniControl - A Mobile Ultrasonic Firewall

abstract: The exchange of data between mobile devices in the near-ultrasonic frequency band is a new promising technology for near field communication (NFC) but also raises a number of privacy concerns. We present the first ultrasonic firewall that reliably detects ultrasonic communication and provides the user with effective means to prevent hidden data exchange. This demonstration showcases a new media-based communication technology ("data over audio") together with its related privacy concerns. It enables users to (i) interactively test out and experience ultrasonic information exchange and (ii) shows how to protect oneself against unwanted tracking.


title: IC4499 revised: spectro-photometric evidence of small light-element variations

abstract: It has been suggested that IC4499 is one of the very few old globulars to not host multiple populations with light-element variations. To follow-up on this very interesting result, here we make use of accurate HST photometry and FLAMES@VLT high-resolution spectroscopy to investigate in more detail the stellar population properties of this system. We find that the red giant branch of the cluster is clearly bimodal in near-UV -- optical colour-magnitude diagrams, thus suggesting that IC4499 is actually composed by two sub-populations of stars with different nitrogen abundances. This represents the first detection of multiple populations in IC4499. Consistently, we also find that one star out of six is Na-rich to some extent, while we do not detect any evidence of intrinsic spread in both Mg and O. The number ratio between stars with normal and enriched nitrogen is in good agreement with the number ratio - mass trend observed in Galactic globular clusters. Also, as typically found in other systems, nitrogen rich stars are more centrally concentrated than normal stars, although this result cannot be considered conclusive because of the limited field of view covered by our observations (~$1r_h$). On the contrary, we observe that both the RGB UV color spread, which is a proxy of N variations, and Na abundance variations, are significantly smaller than those observed in Milky Way globular clusters with mass and metallicity comparable to IC4499. The modest N and Na spreads observed in this system can be tentatively connected to the fact that IC4499 likely formed in a disrupted dwarf galaxy orbiting the Milky Way, as previously proposed based on its orbit.


title: Generalized Metric Repair on Graphs

abstract: Many modern data analysis algorithms either assume that or are considerably more efficient if the distances between the data points satisfy a metric. These algorithms include metric learning, clustering, and dimensionality reduction. Because real data sets are noisy, the similarity measures often fail to satisfy a metric. For this reason, Gilbert and Jain [11] and Fan, et al. [8] introduce the closely related problems of $\textit{sparse metric repair}$ and $\textit{metric violation distance}$. The goal of each problem is to repair as few distances as possible to ensure that the distances between the data points satisfy a metric. We generalize these problems so as to no longer require all the distances between the data points. That is, we consider a weighted graph $G$ with corrupted weights w and our goal is to find the smallest number of modifications to the weights so that the resulting weighted graph distances satisfy a metric. This problem is a natural generalization of the sparse metric repair problem and is more flexible as it takes into account different relationships amongst the input data points. As in previous work, we distinguish amongst the types of repairs permitted (decrease, increase, and general repairs). We focus on the increase and general versions and establish hardness results and show the inherent combinatorial structure of the problem. We then show that if we restrict to the case when $G$ is a chordal graph, then the problem is fixed parameter tractable. We also present several classes of approximation algorithms. These include and improve upon previous metric repair algorithms for the special case when $G = K_n$


title: Border and skewness functions from a leading order fit to DVCS data

abstract: We propose new parameterizations for the border and skewness functions appearing in the description of 3D nucleon structure in the language of Generalized Parton Distributions (GPDs). These parameterizations are constructed in a way to fulfill the basic properties of GPDs, like their reduction to Parton Density Functions and Elastic Form Factors. They also rely on the power behavior of GPDs in the $x \to 1$ limit and the propounded analyticity property of Mellin moments of GPDs. We evaluate Compton Form Factors (CFFs), the sub-amplitudes of the Deeply Virtual Compton Scattering (DVCS) process, at the leading order and leading twist accuracy. We constrain the restricted number of free parameters of these new parameterizations in a global CFF analysis of almost all existing proton DVCS measurements. The fit is performed within the PARTONS framework, being the modern tool for generic GPD studies. A distinctive feature of this CFF fit is the careful propagation of uncertainties based on the replica method. The fit results genuinely permit nucleon tomography and may give some insight into the distribution of forces acting on partons.


title: Approximate Collapsed Gibbs Clustering with Expectation Propagation

abstract: We develop a framework for approximating collapsed Gibbs sampling in generative latent variable cluster models. Collapsed Gibbs is a popular MCMC method, which integrates out variables in the posterior to improve mixing. Unfortunately for many complex models, integrating out these variables is either analytically or computationally intractable. We efficiently approximate the necessary collapsed Gibbs integrals by borrowing ideas from expectation propagation. We present two case studies where exact collapsed Gibbs sampling is intractable: mixtures of Student-t's and time series clustering. Our experiments on real and synthetic data show that our approximate sampler enables a runtime-accuracy tradeoff in sampling these types of models, providing results with competitive accuracy much more rapidly than the naive Gibbs samplers one would otherwise rely on in these scenarios.


title: Distributed Power Control Schemes for In-Band Full-Duplex Energy Harvesting Wireless Networks

abstract: This paper studies two power control problems in energy harvesting wireless networks where one hybrid base station (HBS) and all user equipments (UEs) are operating in in-band full-duplex mode. We consider minimizing the aggregate power subject to the quality of service requirement constraint, and maximizing the aggregate throughput. We address these two problems by proposing two distributed power control schemes for controlling the uplink transmit power by the UEs and the downlink energy harvesting signal power by the HBS. In our proposed schemes, the HBS updates the downlink transmit power level of the energy-harvesting signal so that each UE is enabled to harvest its required energy for powering the operating circuit and transmitting its uplink information signal with the power level determined by the proposed schemes. We show that our proposed power control schemes converge to their corresponding unique fixed points starting from any arbitrary initial transmit power. We will show that our proposed schemes well address the stated problems, which is also demonstrated by our extensive simulation results.


title: Tsallis-INF: An Optimal Algorithm for Stochastic and Adversarial Bandits

abstract: We derive an algorithm that achieves the optimal (within constants) pseudo-regret in both adversarial and stochastic multi-armed bandits without prior knowledge of the regime and time horizon. The algorithm is based on online mirror descent (OMD) with Tsallis entropy regularization with power $\alpha=1/2$ and reduced-variance loss estimators. More generally, we define an adversarial regime with a self-bounding constraint, which includes stochastic regime, stochastically constrained adversarial regime (Wei and Luo), and stochastic regime with adversarial corruptions (Lykouris et al.) as special cases, and show that the algorithm achieves logarithmic regret guarantee in this regime and all of its special cases simultaneously with the adversarial regret guarantee.} The algorithm also achieves adversarial and stochastic optimality in the utility-based dueling bandit setting. We provide empirical evaluation of the algorithm demonstrating that it significantly outperforms UCB1 and EXP3 in stochastic environments. We also provide examples of adversarial environments, where UCB1 and Thompson Sampling exhibit almost linear regret, whereas our algorithm suffers only logarithmic regret. To the best of our knowledge, this is the first example demonstrating vulnerability of Thompson Sampling in adversarial environments. Last, but not least, we present a general stochastic analysis and a general adversarial analysis of OMD algorithms with Tsallis entropy regularization for $\alpha\in[0,1]$ and explain the reason why $\alpha=1/2$ works best.


title: The Limiting Eigenvalue Distribution of Iterated k-Regular Graph Cylinders

abstract: We explore the limiting empirical eigenvalue distributions arising from matrices of the form \[A_{n+1} = \begin{bmatrix} A_n & I\\ I & A_n \end{bmatrix} , \]where $A_0$ is the adjacency matrix of a $k$-regular graph. We find that for bipartite graphs, the distributions are centered symmetric binomial distributions, and for non-bipartite graphs, the distributions are asymmetric. This research grew out of our work on neural networks in $k$-regular graphs. Our original question was whether or not the graph cylinder construction would produce an expander graph that is a suitable candidate for the neural networks being developed at Nousot. This question is answered in the negative for our computational purposes. However, the limiting distribution is still of theoretical interest to us; thus, we present our results here.


title: A comprehensive comparison between APOGEE and LAMOST: Radial Velocities and Atmospheric Stellar Parameters

abstract: We undertake a critical and comprehensive comparison of the radial velocities and the main stellar atmospheric parameters for stars in common between the latest data releases from the APOGEE and the LAMOST surveys. There is a total of 42,420 dwarfs/giants stars in common between the APOGEE DR14 - LAMOST DR3 stellar catalogs. A comparison between the RVs shows an offset of 4.54 $\pm$ 0.03 km/s, with a dispersion of 5.8 km/s, in the sense that APOGEE RVs are larger. We observe a small offset in the Teff of about 13 K, with a scatter of 155 K. Small offset in [Fe/H] of about 0.06 dex together with a scatter of 0.13 dex is also observed. We notice that the largest offset between the surveys occurs in the surface gravities. Using only surface gravities in calibrated red giants from APOGEE DR14, where there are 24,074 stars in common, a deviation of 0.14 dex is found with substantial scatter. There are 17,482 red giant stars in common between APOGEE DR14 and those in LAMOST tied to APOGEE DR12 via the Cannon. There is generally good agreement between the two data-sets. However, we find dependencies of the differences of the stellar parameters on effective temperature. For metal-rich stars, a different trend for the [Fe/H] discrepancies is found. Surprisingly, we see no correlation between the internal APOGEE DR14 - DR12 differences in Teff and those in DR14 - LAMOST tied to DR12, where a correlation should be expected since LAMOST has been calibrated to APOGEE DR12. We also find no correlation between the [Fe/H] discrepancies, suggesting that LAMOST/Cannon is not well coupled to the APOGEE DR12 stellar parameters scale. A [Fe/H] dependence between the stellar parameters in APOGEE DR12 and those in DR14 is reported. We find a weak correlation in the differences between APOGEE DR14 - DR12 and LAMOST on DR12 surface gravity for stars hotter than 4800 K and in the log g range between 2.0 and 2.8.


title: Quasi-polynomial time approximation schemes for packing and covering problems in planar graphs

abstract: We consider two optimization problems in planar graphs. In Maximum Weight Independent Set of Objects we are given a graph $G$ and a family $\mathcal{D}$ of objects, each being a connected subgraph of $G$ with a prescribed weight, and the task is to find a maximum-weight subfamily of $\mathcal{D}$ consisting of pairwise disjoint objects. In Minimum Weight Distance Set Cover we are given an edge-weighted graph $G$, two sets $\mathcal{D},\mathcal{C}$ of vertices of $G$, where vertices of $\mathcal{D}$ have prescribed weights, and a nonnegative radius $r$. The task is to find a minimum-weight subset of $\mathcal{D}$ such that every vertex of $\mathcal{C}$ is at distance at most $r$ from some selected vertex. Via simple reductions, these two problems generalize a number of geometric optimization tasks, notably Maximum Weight Independent Set for polygons in the plane and Weighted Geometric Set Cover for unit disks and unit squares. We present quasi-polynomial time approximation schemes (QPTASs) for both of the above problems in planar graphs: given an accuracy parameter $\epsilon>0$ we can compute a solution whose weight is within multiplicative factor of $(1+\epsilon)$ from the optimum in time $2^{\mathrm{poly}(1/\epsilon,\log |\mathcal{D}|)}\cdot n^{\mathcal{O}(1)}$, where $n$ is the number of vertices of the input graph. Our main technical contribution is to transfer the techniques used for recursive approximation schemes for geometric problems due to Adamaszek, Har-Peled, and Wiese to the setting of planar graphs. In particular, this yields a purely combinatorial viewpoint on these methods.


title: Rapid Time Series Prediction with a Hardware-Based Reservoir Computer

abstract: Reservoir computing is a neural network approach for processing time-dependent signals that has seen rapid development in recent years. Physical implementations of the technique using optical reservoirs have demonstrated remarkable accuracy and processing speed at benchmark tasks. However, these approaches require an electronic output layer to maintain high performance, which limits their use in tasks such as time-series prediction, where the output is fed back into the reservoir. We present here a reservoir computing scheme that has rapid processing speed both by the reservoir and the output layer. The reservoir is realized by an autonomous, time-delay, Boolean network configured on a field-programmable gate array. We investigate the dynamical properties of the network and observe the fading memory property that is critical for successful reservoir computing. We demonstrate the utility of the technique by training a reservoir to learn the short- and long-term behavior of a chaotic system. We find accuracy comparable to state-of-the-art software approaches of similar network size, but with a superior real-time prediction rate up to 160 MHz.


title: Particle Detectors, Cavities, and the Weak Equivalence Principle

abstract: We analyze a quantum version of the weak equivalence principle, in which we compare the response of a static particle detector crossed by an accelerated cavity with the response of an accelerated detector crossing a static cavity in (1+1)-dimensional flat spacetime. We show, for both massive and massless scalar fields, that the non-locality of the field is enough for the detector to distinguish the two scenarios. We find this result holds for vacuum and excited field states of different kinds and we clarify the role of field mass in this setup.


title: Energy-chirp compensation of laser-driven ion beams enabled by structured targets

abstract: We show using 3D simulations that the challenge of generating dense mono-energetic laser-driven ion beams with low angular divergence can be overcome by utilizing structured targets with a relativistically transparent channel and an overdense wall. In contrast to a uniform target that produces a chirped ion beam, the target structure facilitates formation of a dense electron bunch whose longitudinal electric field reverses the energy chirp. This approach works in conjunction with existing acceleration mechanisms, augmenting the ion spectra. For example, our 3D simulations predict a significant improvement for a 2 PW laser pulse with a peak intensity of $5 \times 10^{22}$ W/cm$^2$. The simulations show a mono-energetic proton peak in a highly desirable energy range of 200 MeV with an unprecedented charge of several nC and relatively low divergence that is below 10$^{\circ}$.


title: Renormalisation and locality: branched zeta values

abstract: Multivariate renormalisation techniques are implemented in order to build, study and then renormalise at the poles, branched zeta functions associated with trees. For this purpose, we first prove algebraic results and develop analytic tools, which we then combine to study branched zeta functions. The algebraic aspects concern universal properties for locality algebraic structures, some of which had been discussed in previous work; we "branch/ lift" to trees operators acting on the decoration set of trees, and factorise branched maps through words by means of universal properties for words which we prove in the locality setup. The analytic tools are multivariate meromorphic germs of pseudodifferential symbols with linear poles which generalise the meromorphic germs of functions with linear poles studied in previous work. Multivariate meromorphic germs of pseudodifferential symbols form a locality algebra on which we build various locality maps in the framework of locality structures. We first show that the finite part at infinity defines a locality character from the latter symbol valued meromorphic germs to the scalar valued ones. We further equip the locality algebra of germs of pseudodifferential symbols with locality Rota-Baxter operators given by regularised sums and integrals. By means of the universal properties in the framework of locality structures we can lift Rota-Baxter operators to trees, and use the lifted discrete sums in order to build and study renormalised branched zeta values associated with trees. By construction these renormalised branched zeta values factorise on mutually independent (for the locality relation) trees.


title: Spin Nernst effect in the Paramagnetic Regime of An Antiferromagnetic Insulator

abstract: We theoretically investigate a pure spin Hall current driven by a longitudinal temperature gradient, i.e., the spin Nernst effect (SNE), in a paramagnetic state of a collinear antiferromagnetic insulator with the Dzyaloshinskii-Moriya interaction. The SNE in a magnetic ordered state in such an insulator was proposed by Cheng {\it et al.} [R. Cheng, S. Okamoto, and D. Xiao, Phys. Rev. Lett. 117, 217202 (2016)]. Here we show that the Dzyaloshinskii-Moriya interaction can generate a pure spin Hall current even without magnetic ordering. By using a Schwinger boson mean-field theory, we calculate the temperature dependence of SNE in a disordered phase. We also discuss the implication of our results to experimental realizations.


title: Critique of the Fox-Lu model for Hodgkin-Huxley fluctuations in neuron ion channels

abstract: Using a well known result that every FP equation has an antecedent Langevin equation LE Fox and Lu proposed such a description for ion channels in 1994. Their contraction followed the works of van Kampen and of T. Kurtz. The contraction produces a diffusion term with a state dependent diffusion matrix, D, that arises from the coupling matrix, S, in the LE. This S connected the noise terms to the channel subunit variables in the LE. Fox and Lu and many others later on observed that SS = D. Since D was determined by the contraction of the MC equations into the FP equation, this left the problem of determining the square root matrix, S, for every time step of the simulation. Since this is time consuming, Fox and Lu introduced simplified models not requiring the square root of a matrix. Subsequently, numerous studies were published that showed the several shortcomings of these simplified models. In 2011, Goldwyn et al. [6] rediscovered the overlooked original matrix dependent approach in the Fox-Lu 1994 paper. They showed that it produced results in very good agreement with the MC results. In 1991, Fox and Keizer [7] wrote a paper on an unrelated topic that utilized the work of van Kampen and of Kurtz. In that work the connection between D and S is SST = D. ST is the adjoint (transpose) of S. D remains a positive definite symmetric matrix but S need not be. Fox has reproduced the 2012 results of Orio and Soudry for potassium channels and has also found in closed form the solution for the more complicated sodium channels. The square root problem generally must be done numerically, but the Cholesky is always doable in closed form. Thereby the S matrix for sodium is given explicitly for the first time in this paper


title: Polarisation oscillations in birefringent emitter-cavity systems

abstract: We present the effects of resonator birefringence on the cavity-enhanced interfacing of quantum states of light and matter, including the first observation of single photons with a time-dependent polarisation state that evolves within their coherence time. A theoretical model is introduced and experimentally verified by the modified polarisation of temporally-long single photons emitted from a $^{87}$Rb atom coupled to a high-finesse optical cavity by a vacuum-stimulated Raman adiabatic passage (V-STIRAP) process. Further theoretical investigation shows how a change in cavity birefringence can both impact the atom-cavity coupling and engender starkly different polarisation behaviour in the emitted photons. With polarisation a key resource for encoding quantum states of light and modern micron-scale cavities particularly prone to birefringence, the consideration of these effects is vital to the faithful realisation of efficient and coherent emitter-photon interfaces for distributed quantum networking and communications.


title: Rectification of heat currents across nonlinear quantum chains: A versatile approach beyond weak thermal contact

abstract: Within the emerging field of quantum thermodynamics the issues of heat transfer and heat rectification are basic ingredients for the understanding and design of heat engines or refrigerators at nanoscales. Here, a consistent and versatile approach for mesoscopic devices operating with continuous degrees of freedom is developed valid from low up to strong system-reservoir couplings and over the whole temperature range. It allows to cover weak to moderate nonlinearities and is applicable to various scenarios including the presence of disorder and external time-dependent fields. As a particular application coherent one-dimensional chains of anharmonic oscillators terminated by thermal reservoirs are analyzed with particular focus on rectification. The efficiency of the method opens a door to treat also rather long chains and extensions to higher dimensions and geometries.


title: Multi-Resolution Hashing for Fast Pairwise Summations

abstract: A basic computational primitive in the analysis of massive datasets is summing simple functions over a large number of objects. Modern applications pose an additional challenge in that such functions often depend on a parameter vector $y$ (query) that is unknown a priori. Given a set of points $X\subset \mathbb{R}^{d}$ and a pairwise function $w:\mathbb{R}^{d}\times \mathbb{R}^{d}\to [0,1]$, we study the problem of designing a data-structure that enables sublinear-time approximation of the summation $Z_{w}(y)=\frac{1}{|X|}\sum_{x\in X}w(x,y)$ for any query $y\in \mathbb{R}^{d}$. By combining ideas from Harmonic Analysis (partitions of unity and approximation theory) with Hashing-Based-Estimators [Charikar, Siminelakis FOCS'17], we provide a general framework for designing such data structures through hashing that reaches far beyond what previous techniques allowed. A key design principle is a collection of $T\geq 1$ hashing schemes with collision probabilities $p_{1},\ldots, p_{T}$ such that $\sup_{t\in [T]}\{p_{t}(x,y)\} = \Theta(\sqrt{w(x,y)})$. This leads to a data-structure that approximates $Z_{w}(y)$ using a sub-linear number of samples from each hash family. Using this new framework along with Distance Sensitive Hashing [Aumuller, Christiani, Pagh, Silvestri PODS'18], we show that such a collection can be constructed and evaluated efficiently for any log-convex function $w(x,y)=e^{\phi(\langle x,y\rangle)}$ of the inner product on the unit sphere $x,y\in \mathcal{S}^{d-1}$. Our method leads to data structures with sub-linear query time that significantly improve upon random sampling and can be used for Kernel Density or Partition Function Estimation. We provide extensions of our result from the sphere to $\mathbb{R}^{d}$ and from scalar functions to vector functions.


title: Sparse Averages of Partial Sums of Fourier Series

abstract: We study convergence properties of sparse averages of partial sums of Fourier series of continuous functions. By sparse averages, we are considering an increasing sequences of integers $n_0 < n_1 < n_2 < ...$ and looking at \begin{equation*} \tilde{\sigma}_N(f)(t) = \frac{1}{N+1}\sum_{k=0}^{N}s_{n_k}(f)(t) \end{equation*} to determine the necessary conditions on the sequence $\{n_k\}$ for uniform convergence. Among our results, we find that convergence is dependent on the sequence: we give a proof of convergence for the linear case, $n_k = pk$, for $p$ a positive integer, and present strong experimental evidence for convergence of the quadratic $n_k = k^2$ and cubic $n_k = k^3$ cases, but divergence for the exponential case, $n_k = 2^k$. We also present experimental evidence that if we replace the deterministic rules above by random processes with the same asymptotic behavior then almost surely the answer is the same.


title: Inventor collaboration and its persistence across European regions

abstract: Inventor collaborations that span across regions facilitate the combination and diffusion of innovation-related knowledge. While remote partnerships are gaining importance; most R&D cooperations remain embedded in local environments and little is known about how spatial patterns of creation and persistence of ties in large-scale co-inventor networks differ. We use the publicly available OECD REGPAT database to construct a co-inventor network of the EU27 and continental EFTA countries from patents filed between 2006 and 2010, and identify those ties that had been persisted from before 2006. A community detection exercise reveals that persisted collaborations cluster at a smaller geographical scale than what is observed in the full network. We find that in general the estimated probability of persisted collaboration does not differ from the complete collaboration network when geographical distance, technological similarity or the number of common third partners across regions are assumed to influence tie probability separately. However, persistent collaboration gains exceptional likelihood if regions are proximate in all three dimensions. Our results therefore offer evidence that repeated inventor collaboration drives regional innovation systems towards isolation, which is a threat for European innovation policy.


title: Collision system and beam energy dependence of anisotropic flow fluctuations

abstract: New measurements of two- and four-particle elliptic flow are used to investigate flow fluctuations in collisions of U+U at $\sqrt{s_{NN}}$ = 193~GeV, Cu+Au at $\sqrt{s_{NN}}$ = 200~GeV and Au+Au at several beam energies. These measurements highlight the dependence of these fluctuations on the event-shape, system-size and beam energy and indicate a dominant role for initial-state-driven fluctuations. These measurements could provide further constraints for initial-state models, as well as for precision extraction of the temperature-dependent specific shear viscosity $\frac{\eta}{s}(T)$.


title: Simulations of splashing high and low viscosity droplets

abstract: In this work simulations are presented of low viscosity ethanol and high viscosity silicone oil droplets impacting on a dry solid surface at atmospheric and reduced ambient pressure. The simulations are able to capture both the effect of the ambient gas pressure and liquid viscosity on droplet impact and breakup. The results suggests that the early time droplet impact and gas film behavior for both low and high viscosity liquids share the same physics. However, for later time liquid sheet formation and breakup high and low viscosity liquids behave differently. These results explain why for both kinds of liquids the pressure effect can be observed, while at the same time different high and low viscosity splashing regimes have been identified experimentally.


title: Coloring in Graph Streams

abstract: In this paper, we initiate the study of the vertex coloring problem of a graph in the semi streaming model. In this model, the input graph is defined by a stream of edges, arriving in adversarial order and any algorithm must process the edges in the order of arrival using space linear (up to polylogarithmic factors) in the number of vertices of the graph. In the offline settings, there is a simple greedy algorithm for $(\Delta+1)$-vertex coloring of a graph with maximum degree $\Delta$. We design a one pass randomized streaming algorithm for $(1+\varepsilon)\Delta$-vertex coloring problem for any constant $\varepsilon >0$ using $O(\varepsilon^{-1} n ~\mathrm{ poly} \log n)$ space where $n$ is the number of vertices in the graph. Much more color efficient algorithms are known for graphs with bounded arboricity in the offline settings. Specifically, there is a simple $2\alpha$-vertex coloring algorithm for a graph with arboricity $\alpha$. We present a $O(\varepsilon^{-1}\log n)$ pass randomized vertex coloring algorithm that requires at most $(2+\varepsilon)\alpha$ many colors for any constant $\varepsilon>0$ for a graph with arboricity $\alpha$ in the semi streaming model.


title: An X-ray view of central engines of low luminosity quasars (LLQSO) in the local Universe

abstract: We have carried out a systematic X-ray spectral analysis of a sample of low luminosity quasars (LLQSO) to investigate the nature of the central engines of these sources. The optically-selected LLQSO sample consists of close, known bright active galactic nuclei (AGN) which serves as an important link between the powerful quasars at higher redshift and local Seyfert galaxies. We find that although the bolometric luminosities of the LLQSOs are lower than those of the higher redshift quasars by almost an order of magnitude, their distribution of the Eddington rate $\lambda_{\rm Edd}$ is similar. We detect a strong anti-correlation between $\alpha_{\rm OX}$ and $L_{2500 \rm \AA}$, as has also been detected in several other quasar studies with large sample sizes, indicating that as the UV luminosity of the source increases, the X-ray luminosity decreases. We do not detect any significant neutral obscuration ($N_{\rm H} \ge10^{22}\, \rm cm^{-2}$) in the X-ray spectra of the LLQSOs, and hence rule out obscuration as a possible cause for their lower luminosity. We conclude that the central engines of the LLQSOs function similarly to those of the higher redshift quasars, and the difference is possibly because of the fact that the LLQSOs have lower black hole masses. We do not find any correlation between the molecular gas in the host galaxies and accretion states of the AGN. This indicates that the presence of molecular gas in the host galaxies of the LLQSOs does not significantly influence the instantaneous accretion rates of their SMBHs.


title: Explicit inverse of nonsingular Jacobi matrices

abstract: We present here the necessary and sufficient conditions for the invertibility of tridiagonal matrices, commonly named Jacobi matrices, and explicitly compute their inverse. The techniques we use are related with the solution of Sturm-Liouville boundary value problems associated to second order linear difference equations. These boundary value problems can be expressed throughout a discrete Schr\"odinger operator and their solutions can be computed using recent advances in the study of linear difference equations. The conditions that ensure the uniqueness solution of the boundary value problem lead us to the invertibility conditions for the matrix, whereas the solutions of the boundary value problems provides the entries of the inverse matrix.


title: Physical-type correctness in scientific Python

abstract: The representation of units and dimensions in informatics systems is barely codified and often ignored. For instance, the major languages used in scientific computing (Fortran, C and Python), have no type for dimension or unit, and so physical quantities are represented in a program by variables of type real, resulting in the possibility of unit or dimensional errors. In view of this danger, many authors have proposed language schemes for unit-checking and conversion. However, since many physical quantities have the same units, it is possible for a block of code to be unit-compatible, but still physically meaningless. We demonstrate the limitations of three Python unit-libraries and present a justification and method for checking kind-of-quantity.


title: High-Order Implicit Hybridizable Discontinuous Galerkin Method for the Boltzmann Equation

abstract: The high-order hybridizable discontinuous Galerkin (HDG) method combining with an implicit iterative scheme is used to find the steady-state solution of the Boltzmann equation with full collision integral on two-dimensional triangular meshes. The velocity distribution function and its trace are approximated in the piecewise polynomial space of degree up to 4. The fast spectral method (FSM) is incorporated into the DG discretization to evaluate the collision operator. Specific polynomial approximation is proposed for the collision term to reduce the computational cost. The proposed scheme is proved to be accurate and efficient.


title: Distributed local approximation algorithms for maximum matching in graphs and hypergraphs

abstract: We describe approximation algorithms in Linial's classic LOCAL model of distributed computing to find maximum-weight matchings in a hypergraph of rank $r$. Our main result is a deterministic algorithm to generate a matching which is an $O(r)$-approximation to the maximum weight matching, running in $\tilde O(r \log \Delta + \log^2 \Delta + \log^* n)$ rounds. (Here, the $\tilde O()$ notations hides $\text{polyloglog } \Delta$ and $\text{polylog } r$ factors). This is based on a number of new derandomization techniques extending methods of Ghaffari, Harris & Kuhn (2017). As a main application, we obtain nearly-optimal algorithms for the long-studied problem of maximum-weight graph matching. Specifically, we get a $(1+\epsilon)$ approximation algorithm using $\tilde O(\log \Delta / \epsilon^3 + \text{polylog}(1/\epsilon, \log \log n))$ randomized time and $\tilde O(\log^2 \Delta / \epsilon^4 + \log^*n / \epsilon)$ deterministic time. The second application is a faster algorithm for hypergraph maximal matching, a versatile subroutine introduced in Ghaffari et al. (2017) for a variety of local graph algorithms. This gives an algorithm for $(2 \Delta - 1)$-edge-list coloring in $\tilde O(\log^2 \Delta \log n)$ rounds deterministically or $\tilde O( (\log \log n)^3 )$ rounds randomly. Another consequence (with additional optimizations) is an algorithm which generates an edge-orientation with out-degree at most $\lceil (1+\epsilon) \lambda \rceil$ for a graph of arboricity $\lambda$; for fixed $\epsilon$ this runs in $\tilde O(\log^6 n)$ rounds deterministically or $\tilde O(\log^3 n )$ rounds randomly.


title: Socio-Material Network Analysis: A Mixed Method Study of Five European Artistic Collectives

abstract: In this paper, I argue that we can better understand the relationship between social structure and materiality by combining qualitative analysis of practices in shared physical space with statistical analysis. Drawing on the two-mode approach, I treat social and material structures together with the relationship between them as a two-level socio-material network. In a mixed method study, formalized ethnographic data on such networks in five European artistic collectives are subjected to multilevel exponential random graph modelling. It sheds light on how different types of interpersonal ties condition the engagement of individuals with similar materiality over time.


title: A general fluid-sediment mixture model and constitutive theory

abstract: We present a thermodynamically consistent constitutive model for fluid-saturated sediments, spanning dense to dilute regimes, developed from the basic balance laws for two phase-mixtures. The model can represent various limiting cases, such as pure fluid and dry grains. It is formulated to capture a number of key behaviors such as: (i) viscous inertial rheology of submerged wet grains under steady shearing flows, (ii) the critical state behavior of grains, which causes granular Reynolds dilation/contraction due to shear, (iii) the viscous thickening of the fluid response due to the presence of suspended grains, and (iv) the Darcy-like drag interaction observed in both dense and dilute mixtures, which gives rise to complex fluid-grain interactions under dilation and flow. The full constitutive model is combined with the basic equations of motion for each mixture phase and implemented in the material point method (MPM) to accurately model the coupled dynamics of the mixed system. Qualitative results show the breadth of problems which this model can address. Quantitative results demonstrate the accuracy of this model as compared with analytical limits and experimental observations of fluid and grain behaviors in inhomogeneous geometries.


title: An improved uncertainty principle for functions with symmetry

abstract: Chebotar\"ev proved that every minor of a discrete Fourier matrix of prime order is nonzero. We prove a generalization of this result that includes analogues for discrete cosine and discrete sine matrices as special cases. We establish these results via a generalization of the Bir\'o-Meshulam-Tao uncertainty principle to functions with symmetries that arise from certain group actions, with some of the simplest examples being even and odd functions. We show that our result is best possible and in some cases is stronger than that of Bir\'o-Meshulam-Tao. Some of these results hold in certain circumstances for non-prime fields; Gauss sums play a central role in such investigations.


title: Ergoregion instability of a rotating quantum system

abstract: Using the analogy between acoustic perturbations in an ideal fluid and the description of a Klein-Gordon scalar field in a curved spacetime, we study the quasinormal modes of a quantum system: the rotating Bose-Einstein condensate. To compute quasinormal frequencies, we use two different numerical techniques, namely the direct integration and the continued-fraction methods. We study in detail the ergoregion instability of this linearly perturbed system, comparing the results with different setup configurations.


title: Near-Optimal Distributed Estimation for a Network of Sensing Units Operating Under Communication Constraints

abstract: We study the problem of distributed state estimation in a network of sensing units that can exchange their measurements but the rate of communication between the units is constrained. The units collect noisy, possibly only partial observations of the unknown state; they are assisted by a relay center which can communicate at a higher rate and schedules the exchange of measurements between the units. We consider the task of minimizing the total mean-square estimation error of the network while promoting balance between the individual units' performances. This problem is formulated as the maximization of a monotone objective function subject to a cardinality constraint. By leveraging the notion of weak submodularity, we develop an efficient greedy algorithm for the proposed formulation and show that the greedy algorithm achieves a constant factor approximation of the optimal objective. Our extensive simulation studies illustrate the efficacy of the proposed formulation and the greedy algorithm.


title: First-Passage Duality

abstract: We show that the distribution of times for a diffusing particle to first hit an absorber is \emph{independent} of the direction of an external flow field, when we condition on the event that the particle reaches the target for flow away from the target. Thus, in one dimension, the average time for a particle to travel to an absorber a distance $\ell$ away is $\ell/|v|$, independent of the sign of $v$. This duality extends to all moments of the hitting time. In two dimensions, the distribution of first-passage times to an absorbing circle in the radial velocity field $v(r)=Q/(2\pi r)$ again exhibits duality. Our approach also gives a new perspective on how varying the radial velocity is equivalent to changing the spatial dimension, as well as the transition between transience and strong transience in diffusion.


title: Twisted quantum affinizations and their vertex representations

abstract: In this paper we generalize Drinfeld's twisted quantum affine algebras to construct twisted quantum algebras for all simply-laced generalized Cartan matrices and present their vertex representation realizations.


title: Three-boson spectrum in the presence of 1D spin-orbit coupling: Efimov's generalized radial scaling law

abstract: Spin-orbit coupled cold atom systems, governed by Hamiltonians that contain quadratic kinetic energy terms typical for a particle's motion in the usual Schr\"odinger equation and linear kinetic energy terms typical for a particle's motion in the usual Dirac equation, have attracted a great deal of attention recently since they provide an alternative route for realizing fractional quantum Hall physics, topological insulators, and spintronics physics. The present work focuses on the three-boson system in the presence of 1D spin-orbit coupling, which is most relevant to ongoing cold atom experiments. In the absence of spin-orbit coupling terms, the three-boson system exibits the Efimov effect: the entire energy spectrum is uniquely determined by the $s$-wave scattering length and a single three-body parameter, i.e., using one of the energy levels as input, the other energy levels can be obtained via Efimov's radial scaling law, which is intimately tied to a discrete scaling symmetry. It is demonstrated that the discrete scaling symmetry persists in the presence of 1D spin-orbit coupling, implying the validity of a generalized radial scaling law in five-dimensional space. The dependence of the energy levels on the scattering length, spin-orbit coupling parameters, and center-of-mass momentum is discussed. It is conjectured that three-body systems with other types of spin-orbit coupling terms are also governed by generalized radial scaling laws, provided the system exhibits the Efimov effect in the absence of spin-orbit coupling.


title: Inflation in Multi-field Modified DBM Potentials

abstract: We study multi-field inflation in random potentials generated via a non-equilibrium random matrix theory process. We make a novel modification of the process to include correlations between the elements of the Hessian and the height of the potential, similar to a Random Gaussian Field (RGF). We present the results of over 50,000 inflationary simulations involving 5-100 fields. For the first time, we present results of $\mathcal{O}(100)$ fields using the full `transport method', without slow-roll approximation. We conclude that Planck compatibility is a common prediction of such models, however significant isocurvature power at the end of inflation is possible.


title: Characteristics of unconventional Rb magneto-optical traps

abstract: We study several new magneto-optical trapping configurations in $^{87}$Rb. These unconventional MOTs all use type-II transitions, where the angular momentum of the ground state is greater than or equal to that of the excited state, and they may use either red-detuned or blue-detuned light. We describe the conditions under which each new MOT forms. The various MOTs exhibit an enormous range of lifetimes, temperatures and density distributions. At the detunings where they are maximized, the lifetimes of the various MOTs vary from 0.1 to 15 s. One MOT forms large ring-like structures with no density at the centre. The temperature in the red-detuned MOTs can be three orders of magnitude higher than in the blue-detuned MOTs. We present measurements of the capture velocity of a blue-detuned MOT, and we study how the loss rate due to ultracold collisions depends on laser intensity and detuning.


title: Striving for Realism, not for Determinism: Historical Misconceptions on Einstein and Bohm

abstract: In this paper I show that, while Einstein and Bohm both pursued a deterministic description of quantum mechanics, their philosophical concern was in fact primarily realism and not determinism. Their alleged firm adherence to determinism is based on a long-lasting misleading narrative and it should therefore be greatly scaled down.


title: Performance Enhancement of Soil-Structure Systems Using a Controlled Rocking

abstract: Application of performance-based design (PBD) in earthquake geotechnical design codes has been gaining attention for the past decade. Strong vibrations (e.g., earthquake loading) may generate an uplift or a partial separation of shallow foundations from the underneath soil. To minimize the damage and keep the superstructure safe from the vibrations, a controlled rocking can be exploited. In this study, a finite element method is used to investigate the effect of nonlinear behavior of foundation geomaterial on the rocking behavior of soil-structure systems. To get a better insight about the rocking behavior of the soil-structure system, nonlinear elastic-perfectly plastic behavior is considered for the simulated geomaterial of the foundation. Next, a parametric study, using a wide range of geomaterials with different stiffness values, is conducted in accordance with high-rise structures with different configurations. The results indicate that a significant proportion of the input energy is dissipated by deploying rotational behavior of shallow foundations especially for the structures on softer geomaterials.


title: Deriving star cluster parameters with convolutional neural networks. I. Age, mass, and size

abstract: Context. Convolutional neural networks (CNNs) have been proven to perform fast classification and detection on natural images and have potential to infer astrophysical parameters on the exponentially increasing amount of sky survey imaging data. The inference pipeline can be trained either from real human-annotated data or simulated mock observations. Until now star cluster analysis was based on integral or individual resolved stellar photometry. This limits the amount of information that can be extracted from cluster images. Aims. Develop a CNN-based algorithm aimed to simultaneously derive ages, masses, and sizes of star clusters directly from multi-band images. Demonstrate CNN capabilities on low mass semi-resolved star clusters in a low signal-to-noise ratio regime. Methods. A CNN was constructed based on the deep residual network (ResNet) architecture and trained on simulated images of star clusters with various ages, masses, and sizes. To provide realistic backgrounds, M31 star fields taken from the PHAT survey were added to the mock cluster images. Results. The proposed CNN was verified on mock images of artificial clusters and has demonstrated high precision and no significant bias for clusters of ages $\lesssim$3Gyr and masses between 250 and 4,000 ${\rm M_\odot}$. The pipeline is end-to-end, starting from input images all the way to the inferred parameters; no hand-coded steps have to be performed: estimates of parameters are provided by the neural network in one inferential step from raw images.


title: A consistent estimate for Gravitational Wave and Electromagnetic transient rates

abstract: Gravitational wave transients, resulting from the merger of two stellar remnants, are now detectable. The properties and rates of these directly relates to the stellar population which gave rise to their progenitors, and thus to other, electromagnetic transients which result from stellar death. We aim to estimate simultaneously the event rates and delay time distribution of gravitational wave-driven compact object mergers together with the rates of core collapse and thermonuclear supernovae within a single consistent stellar population synthesis paradigm. We combine event delay-time distributions at different metallicities from the Binary Population and Spectral Synthesis (BPASS) models with an analytic model of the volume-averaged cosmic star formation rate density and chemical evolution to determine the volume-averaged rates of each event rate at the current time. We estimate rates in excellent agreement with extant observational constraints on core-collapse supernovae, thermonuclear supernovae and long GRBs. We predict rates for gravitational wave mergers based on the same stellar populations, and find rates consistent with current LIGO estimates. We note that tighter constraints on the rates of these events will be required before it is possible to determine their redshift evolution, progenitor metallicity dependence or constrain uncertain aspects of stellar evolution.


title: An independent search of gravitational waves in the first observation run of advanced LIGO using cross-correlation

abstract: This work describes a template-free method to search gravitational waves (GW) using data from the LIGO observatories simultaneously. The basic idea of this method is that a GW signal is present in a short-duration data segment if the maximum correlation-coefficient between the strain signals is higher than a significant threshold and its time difference is lower than the 10 ms of inter-observatory light propagation time. Hence, this method can be used to carry out blind searches of any types of GW irrespective of the waveform and of the source type and sky location. An independent search of injected and real GW signals from compact binary coalescences (CBC) contained in the first observation run (O1) of advanced LIGO was carried out to asses its performance. On the basis of the results, the proposed method was able to detect GW produced by binary systems without making any assumption about them.


title: Field theory actions for ambitwistor string and superstring

abstract: We analyze the free ambitwistor string field theory action for the bosonic string, heterotic string and both GSO sectors of the Type II string. The spectrum contains non-unitary states and provides an interesting consistency test for one-loop ambitwistor string computations.


title: Loud and Interactive Paper Prototyping in Requirements Elicitation: What is it Good for?

abstract: Requirements Engineering is a multidisciplinary and a human-centered process, therefore, the artifacts produced from RE are always error-prone. The most significant of these errors are missing or misunderstanding requirements. Information loss in RE could result in omitted logic in the software, which will be onerous to correct at the later stages of development. In this paper, we demonstrate and investigate how interactive and Loud Paper Prototyping (LPP) can be integrated to collect stakeholders' needs and expectations than interactive prototyping or face-to-face meetings alone. To this end, we conducted a case study of (1) 31 mobile application (App) development teams who applied either of interactive or loud prototyping and (2) 19 mobile App development teams who applied only the face-to-face meetings. From this study, we found that while using Silent Paper Prototyping (SPP) rather than No Paper Prototyping (NPP) is a more efficient technique to capture Non-Functional Requirements (NFRs), User Interface (UI) requirements, and existing requirements, LPP is more applicable to manage NFRs, UI requirements, as well as adding new requirements and removing/modifying the existing requirements. We also found that among LPP and SPP, LPP is more efficient to capture and influence Functional Requirements (FRs).


title: Automatically Designing CNN Architectures for Medical Image Segmentation

abstract: Deep neural network architectures have traditionally been designed and explored with human expertise in a long-lasting trial-and-error process. This process requires huge amount of time, expertise, and resources. To address this tedious problem, we propose a novel algorithm to optimally find hyperparameters of a deep network architecture automatically. We specifically focus on designing neural architectures for medical image segmentation task. Our proposed method is based on a policy gradient reinforcement learning for which the reward function is assigned a segmentation evaluation utility (i.e., dice index). We show the efficacy of the proposed method with its low computational cost in comparison with the state-of-the-art medical image segmentation networks. We also present a new architecture design, a densely connected encoder-decoder CNN, as a strong baseline architecture to apply the proposed hyperparameter search algorithm. We apply the proposed algorithm to each layer of the baseline architectures. As an application, we train the proposed system on cine cardiac MR images from Automated Cardiac Diagnosis Challenge (ACDC) MICCAI 2017. Starting from a baseline segmentation architecture, the resulting network architecture obtains the state-of-the-art results in accuracy without performing any trial-and-error based architecture design approaches or close supervision of the hyperparameters changes.


title: Enhanced dynamic nuclear polarization via swept microwave frequency combs

abstract: Dynamic Nuclear Polarization (DNP) has enabled enormous gains in magnetic resonance signals and led to vastly accelerated NMR/MRI imaging and spectroscopy. Unlike conventional cw-techniques, DNP methods that exploit the full electron spectrum are appealing since they allow direct participation of all electrons in the hyperpolarization process. Such methods typically entail sweeps of microwave radiation over the broad electron linewidth to excite DNP, but are often inefficient because the sweeps, constrained by adiabaticity requirements, are slow. In this paper we develop a technique to overcome the DNP bottlenecks set by the slow sweeps, employing a swept microwave frequency comb that increases the effective number of polarization transfer events while respecting adiabaticity constraints. This allows a multiplicative gain in DNP enhancement, scaling with the number of comb frequencies and limited only by the hyperfine-mediated electron linewidth. We demonstrate the technique for the optical hyperpolarization of 13C nuclei in powdered microdiamonds at low fields, increasing the DNP enhancement from 30 to 100 measured with respect to the thermal signal at 7T. For low concentrations of broad linewidth electron radicals, e.g. TEMPO, these multiplicative gains could exceed an order of magnitude.

categories: physics.app-ph physics.atom-ph physics.chem-ph quant-ph


title: Hierarchical Reinforcement Learning for Zero-shot Generalization with Subtask Dependencies

abstract: We introduce a new RL problem where the agent is required to generalize to a previously-unseen environment characterized by a subtask graph which describes a set of subtasks and their dependencies. Unlike existing hierarchical multitask RL approaches that explicitly describe what the agent should do at a high level, our problem only describes properties of subtasks and relationships among them, which requires the agent to perform complex reasoning to find the optimal subtask to execute. To solve this problem, we propose a neural subtask graph solver (NSGS) which encodes the subtask graph using a recursive neural network embedding. To overcome the difficulty of training, we propose a novel non-parametric gradient-based policy, graph reward propagation, to pre-train our NSGS agent and further finetune it through actor-critic method. The experimental results on two 2D visual domains show that our agent can perform complex reasoning to find a near-optimal way of executing the subtask graph and generalize well to the unseen subtask graphs. In addition, we compare our agent with a Monte-Carlo tree search (MCTS) method showing that our method is much more efficient than MCTS, and the performance of NSGS can be further improved by combining it with MCTS.


title: Anomaly constraint on massless QCD and the role of Skyrmions in chiral symmetry breaking

abstract: We discuss consequences of the 't Hooft anomaly matching condition for Quantum Chromodynamics (QCD) with massless fundamental quarks. We derive the new discrete 't Hooft anomaly of massless QCD for generic numbers of color $N_\mathrm{c}$ and flavor $N_\mathrm{f}$, and an exotic chiral-symmetry broken phase without quark-bilinear condensate is ruled out from possible QCD vacua. We show that the $U(1)_\mathrm{B}$ baryon number symmetry is anomalously broken when the $(\mathbb{Z}_{2N_\mathrm{f}})_\mathrm{A}$ discrete axial symmetry and the flavor symmetry are gauged. In the ordinary chiral symmetry breaking, the Skyrmion current turns out to reproduce this 't Hooft anomaly of massless QCD. In the exotic chiral symmetry breaking, however, the anomalous breaking of $U(1)_\mathrm{B}$ does not take the correct form, and it is inconsistent with anomaly matching. This no-go theorem is based only on symmetries and anomalies, and thus has a wider range of applicability to the QCD phase diagram than the previous one obtained by QCD inequalities. Lastly, as another application, we check that duality of $\mathcal{N}=1$ supersymmetric QCD with $N_\mathrm{f}\ge N_\mathrm{c}+1$ satisfies the new anomaly matching.


title: Simple security proof of twin-field type quantum key distribution protocol

abstract: Twin-field (TF) quantum key distribution (QKD) was conjectured to beat the private capacity of a point-to-point QKD link by using single-photon interference in a central measuring station. This remarkable conjecture has recently triggered an intense research activity to prove its security. Here, we introduce a TF-type QKD protocol which is conceptually simpler than the original proposal. It relies on local phase randomization, instead of global phase randomization, which significantly simplifies its security analysis and is arguably less demanding experimentally. We demonstrate that the secure key rate of our protocol has a square-root improvement over the point-to-point private capacity, as conjectured by the original TF-QKD scheme.


title: Survival of non-coplanar, closely-packed planetary systems after a close encounter

abstract: Planetary systems with more than two bodies will experience orbital crossings at a time related to the initial orbital separations of the planets. After a crossing, the system enters a period of chaotic evolution ending in the reshaping of the system's architecture via planetary collisions or ejections. We carry out N-body integrations on a large number of systems with equally-spaced planets (in units of the Hill radius) to determine the distribution of instability times for a given planet separation. We investigate both the time to the initiation of instability through a close encounter and the time to a planet-planet collision. We find that a significant portion of systems with non-zero mutual inclinations survive after a close encounter and do not promptly experience a planet-planet collision. Systems with significant inclinations can continue to evolve for over 1,000 times longer than the encounter time. The fraction of long lived systems is dependent on the absolute system scale and the initial inclination of the planets. These results have implications to the assumed stability of observed planetary systems.


title: Modelling of consciousness and interpretation of quantum mechanics

abstract: I start from the fundamental principles of non-relativistic quantum mechanics, without probability, and interpret them using the notion of coexistence: a quantum state can be read, not uniquely, as a coexistence of other quantum states, which are pairwise orthogonal. In this formalism, I prove that a conscious observer is necessarily a physical object that can memorize local events by setting one of its parts in an exactly specified constant quantum state (hypotheses H1, H2 and H3). Then I define the probability of a future event as the proportion of initial observers, all identical, who will actually experience that event. It then becomes possible to establish the usual results of quantum mechanics. Furthermore, I detail the link between probabilities and relative frequencies. Additionally, I study the biological feasibility of this modelling of observer's mind. The second part of this paper completes the neuronal description of the mind functions, based on current neuroscientific knowledge. It provides a model that is compatible with the assumptions of the first part and consistent with our daily conscious experience. In particular, it develops a model of self-consciousness based on an explicit use of the random component of neuron behaviour; according to the first part, that random is in fact the coexistence of a multiplicity of possibilities. So, when the mind measures the random part of certain neurons in the brain, he goes himself within each of these possibilities. The mind has a decision-making component that is active in this situation, appearing then as the cause of the choice of this possibility among all the others. This models the self-consciousness which then ensures the unity of our conscious experience by equating this experience with ``what the ego is conscious about''. The conclusion details the points that remain to be developed.


title: Statistical generalized derivative applied to the profile likelihood estimation in a mixture of semiparametric models

abstract: There is a difficulty in finding an estimate of variance of the profile likelihood estimator in the joint model of longitudinal and survival data. We solve the difficulty by introducing the ``statistical generalized derivative''. The derivative is used to show the asymptotic normality of the estimator without assuming the second derivative of the density function in the model exists.


title: $\mathcal{N}=2$ Supersymmetry with Central Charge: A Twofold Implementation

abstract: In this work, we analyze an extended $\mathcal{N}=2$ supersymmetry with central charge and develop its superspace formulation under two distinct viewpoints. Initially, in the context of classical mechanics, we discuss the introduction of deformed supersymmetric derivatives and their consequence on the deformation of one-dimensional non-linear sigma model. After that, considering a field-theoretical framework, we present an implementation of this superalgebra in two dimensions, such that one of the coordinates is related to the central charge. As an application, in this two-dimensional scenario, we consider topological (bosonic) configurations of a special self-coupled matter model and present a non-trivial fermionic solution.


title: Quantum Circuit Model of Black Hole Evaporation

abstract: We consider a quantum circuit model describing the evaporation process of black holes. We specifically examine the behavior of the multipartite entanglement represented by this model, and find that the entanglement structure depends on the black hole mass $M$ and the frequency of the Hawking radiation $\omega$. For sufficiently small values of $M\omega$, the black hole and the radiation system becomes a separable state after the Page time and a firewall-like structure appears. On the contrary, for larger values of $M\omega$, the entanglement between the black hole and the radiation is not lost. These behaviors imply that owing to the monogamy property of the multipartite entanglement, low frequency modes of the Hawking radiation destroys the quantum correlation between the black hole and the emitted Hawking radiation.


title: Generalized cosmological term in D=4 supergravity from a new AdS-Lorentz superalgebra

abstract: A new supersymmetrization of the so-called AdS-Lorentz algebra is presented. It involves two fermionic generators and is obtained by performing an abelian semigroup expansion of the superalgebra osp(4|1). The peculiar properties of the aforesaid expansion method are then exploited to construct a D=4 supergravity action involving a generalized supersymmetric cosmological term in a geometric way, only from the curvatures of the novel superalgebra. The action obtained with this procedure is a MacDowell-Mansouri like action. Gauge invariance and supersymmetry of the action are also analyzed.


title: Bounding Box Embedding for Single Shot Person Instance Segmentation

abstract: We present a bottom-up approach for the task of object instance segmentation using a single-shot model. The proposed model employs a fully convolutional network which is trained to predict class-wise segmentation masks as well as the bounding boxes of the object instances to which each pixel belongs. This allows us to group object pixels into individual instances. Our network architecture is based on the DeepLabv3+ model, and requires only minimal extra computation to achieve pixel-wise instance assignments. We apply our method to the task of person instance segmentation, a common task relevant to many applications. We train our model with COCO data and report competitive results for the person class in the COCO instance segmentation task.


title: Hard bremsstrahlung from a high-voltage atmospheric discharge and its anisotropy

abstract: The results of the experiments on recording hard gamma radiation and measurements of its angular distribution at the initial stage of a laboratory high-voltage atmospheric discharge are presented. The experiments were performed on an ERG installation at a voltage of $\sim 1$ MV, an atmospheric discharge current of up to 12 kA, and a gap of 0.55 m. The duration of the voltage pulse was about 1~$\mu$s with a pulse rise time of 150-200 ns. The radiation was recorded by an assembly of 10 identical scintillation detectors installed each 10$^\circ$ around the circumference of a quarter of a circle with a curvature of 1 m. In order to separate the radiation with energies from 20 keV to 1.5 MeV, Al and Pb filters of different thicknesses were used. The obtained results show that, as a rule, a multi-beam radiation pattern and several bursts of radiation (each with a directional pattern) are recorded in each shot. In a considerable number of "shots", hard radiation with photon energies comparable to or exceeding the maximum electron energy corresponding to the applied voltage is recorded. In these cases, a needle-like radiation pattern is observed, including at large angles to the axis of the discharge. This may indicate the acceleration of electrons in different plasma channels.


title: Poster: Improving Bug Localization with Report Quality Dynamics and Query Reformulation

abstract: Recent findings from a user study suggest that IR-based bug localization techniques do not perform well if the bug report lacks rich structured information such as relevant program entity names. On the contrary, excessive structured information such as stack traces in the bug report might always not be helpful for the automated bug localization. In this paper, we conduct a large empirical study using 5,500 bug reports from eight subject systems and replicating three existing studies from the literature. Our findings (1) empirically demonstrate how quality dynamics of bug reports affect the performances of IR-based bug localization, and (2) suggest potential ways (e.g., query reformulations) to overcome such limitations.


title: Circuit Complexity for Coherent States

abstract: We examine the circuit complexity of coherent states in a free scalar field theory, applying Nielsen's geometric approach as in [1]. The complexity of the coherent states have the same UV divergences as the vacuum state complexity and so we consider the finite increase of the complexity of these states over the vacuum state. One observation is that generally, the optimal circuits introduce entanglement between the normal modes at intermediate stages even though our reference state and target states are not entangled in this basis. We also compare our results from Nielsen's approach with those found using the Fubini-Study method of [2]. For general coherent states, we find that the complexities, as well as the optimal circuits, derived from these two approaches, are different.


title: Arithmetic aspects of symmetric edge polytopes

abstract: We investigate arithmetic, geometric and combinatorial properties of symmetric edge polytopes. We give a complete combinatorial description of their facets. By combining Gr\"obner basis techniques, half-open decompositions and methods for interlacing polynomials we provide an explicit formula for the $h^\ast$-polynomial in case of complete bipartite graphs. In particular, we show that the $h^\ast$-polynomial is $\gamma$-positive and real-rooted. This proves Gal's conjecture for arbitrary flag unimodular triangulations in this case, and, beyond that, we prove a strengthing due to Nevo and Petersen (2011).


title: Blue straggler stars beyond the Milky Way. III. Detection of evolved blue straggler candidates in Large Magellanic Cloud clusters

abstract: We analyze {\sl Hubble Space Telescope} observations of nine Large Magellanic Cloud star clusters with ages of 1--2 Gyr to search for evolved counterparts of blue straggler stars. Near the red clump regions in the clusters' color--magnitude diagrams, we find branches of evolved stars that are much brighter than normal evolved stars. We examine the effects of photometric artifacts, differential reddening, and field contamination. We conclude that these bright evolved stars cannot be explained by any of these effects. Our statistical tests show that the contributions of photometric uncertainties and crowding effects, as well as that owing to differential reddening, to these bright evolved stars are insufficient to fully explain the presence of these bright evolved stars. Based on isochrone fitting, we also ruled out the possibility that these bright evolved stars could be reproduced by an internal chemical abundance spread. The spatial distributions of the bright evolved stars exhibit clear concentrations that cannot be explained by homogeneously distributed field stars. This is further confirmed based on Monte Carlo-based tests. By comparing our observations with stellar evolution models, we find that the masses of most of bright evolved stars do not exceed twice the average mass of normal evolved stars. We suggest that these bright evolved stars are, in fact, evolved blue straggler stars.


title: Generalized Stochastic Frank-Wolfe Algorithm with Stochastic "Substitute" Gradient for Structured Convex Optimization

abstract: The stochastic Frank-Wolfe method has recently attracted much general interest in the context of optimization for statistical and machine learning due to its ability to work with a more general feasible region. However, there has been a complexity gap in the guaranteed convergence rate for stochastic Frank-Wolfe compared to its deterministic counterpart. In this work, we present a new generalized stochastic Frank-Wolfe method which closes this gap for the class of structured optimization problems encountered in statistical and machine learning characterized by empirical loss minimization with a certain type of ``linear prediction'' property (formally defined in the paper), which is typically present loss minimization problems in practice. Our method also introduces the notion of a ``substitute gradient'' that is a not-necessarily-unbiased sample of the gradient. We show that our new method is equivalent to a particular randomized coordinate mirror descent algorithm applied to the dual problem, which in turn provides a new interpretation of randomized dual coordinate descent in the primal space. Also, in the special case of a strongly convex regularizer our generalized stochastic Frank-Wolfe method (as well as the randomized dual coordinate descent method) exhibits linear convergence. Furthermore, we present computational experiments that indicate that our method outperforms other stochastic Frank-Wolfe methods consistent with the theory developed herein.


title: Optimal Co-design of Industrial Networked Control Systems with State-dependent Correlated Fading Channels

abstract: This paper examines a co-design problem for industrial networked control systems (NCS) whereby physical systems are controlled over wireless fading channels. In particular, the considered wireless channels are also stochastically dependent on the physical states of moving machineries in an industrial working space. In this paper, the moving machineries are modeled as Markov decision processes whereas the characteristics of the correlated fading channels are modeled as a binary random process whose probability measure is a function of both the moving machineries' physical states and the channels' transmission power. Under such state-dependent fading channel models, sufficient conditions which ensure the stochastic safety of the NCS are first derived. Using the derived safety conditions, the co-design problem is then formulated as a constrained joint optimization problem that seeks for optimal control and transmission power policies which simultaneously minimize both the communication and control costs in an infinite time horizon. This paper shows that such optimal co-design policies can be obtained in an efficient manner from the solution of convex programs. Simulation results from industrial NCS models consisting forklift truck and networked DC motor systems are also presented to illustrate and verify both the advantages and efficacy of the proposed co-design solution framework.


title: Cross-Modality Distillation: A case for Conditional Generative Adversarial Networks

abstract: In this paper, we propose to use a Conditional Generative Adversarial Network (CGAN) for distilling (i.e. transferring) knowledge from sensor data and enhancing low-resolution target detection. In unconstrained surveillance settings, sensor measurements are often noisy, degraded, corrupted, and even missing/absent, thereby presenting a significant problem for multi-modal fusion. We therefore specifically tackle the problem of a missing modality in our attempt to propose an algorithm based on CGANs to generate representative information from the missing modalities when given some other available modalities. Despite modality gaps, we show that one can distill knowledge from one set of modalities to another. Moreover, we demonstrate that it achieves better performance than traditional approaches and recent teacher-student models.


title: The Rickman-Picard Theorem

abstract: We give a new and conceptually simple proof of the Rickman-Picard theorem for quasiregular maps based on potential-theoretic methods.


title: The Surface Magnetic Activity of the Weak-Line T Tauri Stars TWA 9A and V1095 Sco

abstract: We present a detailed analysis of high-resolution spectropolarimetric observations of the weak-line T Tauri stars (wTTSs) TWA 9A and V1095 Sco as part of a wider sur- vey of magnetic properties and activity in weak-line T Tauri stars, called MaTYSSE (Magnetic Topologies of Young Stars and the Survival of close-in giant Exoplanets). Our targets have similar masses but differing ages which span the stage of radiative core formation in solar-mass stars. We use the intensity line profiles to reconstruct surface brightness maps for both stars. The reconstructed brightness maps for TWA 9A and V1095 Sco were used to model and subtract the activity-induced jitter, reducing the RMS in the radial velocity measurements of TWA 9A by a factor of $\sim$7, and for V1095 Sco by a factor of $\sim$3. We detect significant circular polarisation for both stars, but only acquired a high quality circular polarisation time-series for TWA 9A. Our reconstructed large-scale magnetic field distribution of TWA 9A indicates a strong, non-axisymmetric field. We also analyse the chromospheric activity of both stars by investigating their H$\alpha$ emission, finding excess blue-ward emission for most observations of V1095 Sco, and symmetric, double-peaked emission for TWA 9A, with enhanced emission at one epoch likely indicating a flaring event.