title: Towards a Theory of Parameterized Streaming Algorithms

abstract: Parameterized complexity attempts to give a more fine-grained analysis of the complexity of problems: instead of measuring the running time as a function of only the input size, we analyze the running time with respect to additional parameters. This approach has proven to be highly successful in delineating our understanding of \NP-hard problems. Given this success with the TIME resource, it seems but natural to use this approach for dealing with the SPACE resource. First attempts in this direction have considered a few individual problems, with some success: Fafianie and Kratsch [MFCS'14] and Chitnis et al. [SODA'15] introduced the notions of streaming kernels and parameterized streaming algorithms respectively. For example, the latter shows how to refine the $\Omega(n^2)$ bit lower bound for finding a minimum Vertex Cover (VC) in the streaming setting by designing an algorithm for the parameterized $k$-VC problem which uses $O(k^{2}\log n)$ bits. In this paper, we initiate a systematic study of graph problems from the paradigm of parameterized streaming algorithms. We first define a natural hierarchy of space complexity classes of FPS, SubPS, SemiPS, SupPS and BrutePS, and then obtain tight classifications for several well-studied graph problems such as Longest Path, Feedback Vertex Set, Dominating Set, Girth, Treewidth, etc. into this hierarchy. (see paper for full abstract)


title: A class of non-weight modules over the super-BMS$_3$ algebra

abstract: In the present paper, a class of non-weight modules over the super-BMS$_3$ algebras $\S^{\epsilon}$ ($\epsilon=0$ or $\frac{1}{2}$) are constructed. Assume that $\mathfrak{t}=\C L_0\oplus\C W_0\oplus\C G_0$ and $\mathfrak{T}=\C L_0\oplus\C W_0$ are the Cartan subalgebra (modulo center) of $\S^{0}$ and $\S^{\frac{1}{2}}$, respectively. These modules over $\S^{0}$ when restricted to the $\mathfrak{t}$ are free of rank $1$, while these modules over $\S^{\frac{1}{2}}$ when restricted to the $\mathfrak{T}$ are free of rank $2$. Then we determine the necessary and sufficient conditions for these modules being simple, as well as determining the necessary and sufficient conditions for two $\S^{\epsilon}$-modules being isomorphic. %Moreover, we see that the category of free $U(\mathfrak{t})$-modules of rank $1$ over $\S^0$ is %equivalent to the category of free $U(\mathfrak{T})$-modules of rank $2$ over %$\S^{\frac{1}{2}}$.


title: Unsupervised Domain Adaptation by Optical Flow Augmentation in Semantic Segmentation

abstract: It is expensive to generate real-life image labels and there is a domain gap between real-life and simulated images, hence a model trained on the latter cannot adapt to the former. Solving this can totally eliminate the need for labeling real-life datasets completely. Class balanced self-training is one of the existing techniques that attempt to reduce the domain gap. Moreover, augmenting RGB with flow maps has improved performance in simple semantic segmentation and geometry is preserved across domains. Hence, by augmenting images with dense optical flow map, domain adaptation in semantic segmentation can be improved.


title: Thermal wave crystals based on dual-phase-lag model

abstract: Thermal wave crystals based on the dual-phase-lag model are investigated in this paper by both theoretical analysis and numerical simulation to control the non-Fourier heat conduction process. The transfer matrix method is used to calculate the complex dispersion curves. The temperature field is calculated by the finite difference time domain method. The results show that thermal band-gaps exist due to the Bragg-scattering. The key parameters characterizing the band-gaps are analyzed. The thermal wave impedance and mid-gap frequencies are introduced to predict band-gaps theoretically. Our results show that the larger the difference in the thermal wave impedances is, the wider of the thermal band-gaps will be. This study demonstrates a type of the thermal metamaterials which have potential innovative applications such as thermal imagining, thermal diodes and thermal waveguides for energy transmission.


title: Possible instabilities in quadratic and cubic nodal line fermion systems with correlated interactions

abstract: Influence of short-range four-fermion interactions on quadratic and cubic nodal line fermion systems is studied by renormalization group theory. It is found that arbitrarily weak four-fermion interaction could drive quadratic or cubic nodal line fermion system to a new phase. According to the initial conditions and value of fermion flavor, the system may appear three kinds of instabilities. First, quadratic or cubic nodal line is split into conventional nodal lines, thus the system becomes nodal line semimetal. Second, finite excitonic gap is generated, and the system becomes an excitonic insulator. Third, the system is driven into superconducting phase. Thus, quadratic and cubic nodal line fermion systems are rare strong correlated fermion systems in three dimension under the influence of four-fermion interactions. These theoretical results may be verified in the candidates for quadratic and cubic nodal line fermion systems.


title: Temporal Reasoning via Audio Question Answering

abstract: Multimodal question answering tasks can be used as proxy tasks to study systems that can perceive and reason about the world. Answering questions about different types of input modalities stresses different aspects of reasoning such as visual reasoning, reading comprehension, story understanding, or navigation. In this paper, we use the task of Audio Question Answering (AQA) to study the temporal reasoning abilities of machine learning models. To this end, we introduce the Diagnostic Audio Question Answering (DAQA) dataset comprising audio sequences of natural sound events and programmatically generated questions and answers that probe various aspects of temporal reasoning. We adapt several recent state-of-the-art methods for visual question answering to the AQA task, and use DAQA to demonstrate that they perform poorly on questions that require in-depth temporal reasoning. Finally, we propose a new model, Multiple Auxiliary Controllers for Linear Modulation (MALiMo) that extends the recent Feature-wise Linear Modulation (FiLM) model and significantly improves its temporal reasoning capabilities. We envisage DAQA to foster research on AQA and temporal reasoning and MALiMo a step towards models for AQA.


title: Bayesian forecasting of multivariate time series: Scalability, structure uncertainty and decisions

abstract: I overview recent research advances in Bayesian state-space modeling of multivariate time series. A main focus is on the decouple/recouple concept that enables application of state-space models to increasingly large-scale data, applying to continuous or discrete time series outcomes. The scope includes large-scale dynamic graphical models for forecasting and multivariate volatility analysis in areas such as economics and finance, multi-scale approaches for forecasting discrete/count time series in areas such as commercial sales and demand forecasting, and dynamic network flow models for areas including internet traffic monitoring. In applications, explicit forecasting, monitoring and decision goals are paramount and should factor into model assessment and comparison, a perspective that is highlighted.


title: Stochastic Modelling in Fluid Dynamics: It\^o vs Stratonovich

abstract: Suppose the observations of Lagrangian trajectories for fluid flow in some physical situation can be modelled sufficiently accurately by a spatially correlated It\^o stochastic process (with zero mean) obtained from data which is taken in fixed Eulerian space. Suppose we also want to apply Hamilton's principle to derive the stochastic fluid equations for this situation. Now, the variational calculus for applying Hamilton's principle requires the Stratonovich process, so we must transform from It\^o noise in the \emph{data frame} to the equivalent Stratonovich noise. However, the transformation from the It\^o process in the data frame to the corresponding Stratonovich process shifts the drift velocity of the transformed Lagrangian fluid trajectory out of the data frame into a non-inertial frame obtained from the It\^o correction. The issue is, "Will non-inertial forces arising from this transformation of reference frames make a difference in the interpretation of the solution behaviour of the resulting stochastic equations?" This issue will be resolved by elementary considerations.


title: 3D MHD Modeling of the Impact of Subsurface Stratification on the Solar Dynamo

abstract: Various models of solar subsurface stratification are tested in the global EULAG-MHD solver to simulate diverse regimes of near-surface convective transport. Sub- and superadiabacity are altered at the surface of the model ($ r > 0.95~R_{\odot}$) to either suppress or enhance convective flow speeds in an effort to investigate the impact of the near-surface layer on global dynamics. A major consequence of increasing surface convection rates appears to be a significant alteration of the distribution of angular momentum, especially below the tachocline where the rotational frequency predominantly increases at higher latitudes. These hydrodynamic changes correspond to large shifts in the development of the current helicity in this stable layer ($r<0.72R_{\odot}$), significantly altering its impact on the generation of poloidal and toroidal fields at the tachocline and below, acting as a major contributor towards transitions in the dynamo cycle. The enhanced near-surface flow speed manifests in a global shift of the toroidal field ($B_{\phi}$) in the butterfly diagram - from a North-South symmetric pattern to a staggered anti-symmetric emergence.


title: AdaFilter: Adaptive Filter Fine-tuning for Deep Transfer Learning

abstract: There is an increasing number of pre-trained deep neural network models. However, it is still unclear how to effectively use these models for a new task. Transfer learning, which aims to transfer knowledge from source tasks to a target task, is an effective solution to this problem. Fine-tuning is a popular transfer learning technique for deep neural networks where a few rounds of training are applied to the parameters of a pre-trained model to adapt them to a new task. Despite its popularity, in this paper, we show that fine-tuning suffers from several drawbacks. We propose an adaptive fine-tuning approach, called AdaFilter, which selects only a part of the convolutional filters in the pre-trained model to optimize on a per-example basis. We use a recurrent gated network to selectively fine-tune convolutional filters based on the activations of the previous layer. We experiment with 7 public image classification datasets and the results show that AdaFilter can reduce the average classification error of the standard fine-tuning by 2.54%.


title: Estimating uncertainty of earthquake rupture using Bayesian neural network

abstract: Bayesian neural networks (BNN) are the probabilistic model that combines the strengths of both neural network (NN) and stochastic processes. As a result, BNN can combat overfitting and perform well in applications where data is limited. Earthquake rupture study is such a problem where data is insufficient, and scientists have to rely on many trial and error numerical or physical models. Lack of resources and computational expenses, often, it becomes hard to determine the reasons behind the earthquake rupture. In this work, a BNN has been used (1) to combat the small data problem and (2) to find out the parameter combinations responsible for earthquake rupture and (3) to estimate the uncertainty associated with earthquake rupture. Two thousand rupture simulations are used to train and test the model. A simple 2D rupture geometry is considered where the fault has a Gaussian geometric heterogeneity at the center, and eight parameters vary in each simulation. The test F1-score of BNN (0.8334), which is 2.34% higher than plain NN score. Results show that the parameters of rupture propagation have higher uncertainty than the rupture arrest. Normal stresses play a vital role in determining rupture propagation and are also the highest source of uncertainty, followed by the dynamic friction coefficient. Shear stress has a moderate role, whereas the geometric features such as the width and height of the fault are least significant and uncertain.


title: Paraphrasing with Large Language Models

abstract: Recently, large language models such as GPT-2 have shown themselves to be extremely adept at text generation and have also been able to achieve high-quality results in many downstream NLP tasks such as text classification, sentiment analysis and question answering with the aid of fine-tuning. We present a useful technique for using a large language model to perform the task of paraphrasing on a variety of texts and subjects. Our approach is demonstrated to be capable of generating paraphrases not only at a sentence level but also for longer spans of text such as paragraphs without needing to break the text into smaller chunks.


title: Ideal Solar Cell Efficiencies

abstract: In a recent paper, Guillemoles et al [J.F. Guillemoles, T. Kirchartz, D. Cahen and U Rau, Guide for the perplexed to the Shockley-Queisser model for solar cells, Nat. Photonics, 13, 501 (2019)] attempt to clarify and explain the often cited paper by Shockley and Queisser (SQ) which defines the limits to photovoltaic conversion by a single-junction solar cell. The SQ paper is not easy to read and is therefore easily misunderstood. As modern solar cells approach theoretical efficiency limits, the fundamentals become particularly important and the effort by Guillemoles et al is therefore to be welcome. However, in doing so, the authors have fallen into several pitfalls and the aim of the present note is to clarify a number of misconceptions and correct some errors in that paper.


title: Entangleability of cones

abstract: We solve a long-standing conjecture by Barker, proving that the minimal and maximal tensor products of two finite-dimensional proper cones coincide if and only if one of the two cones is generated by a linearly independent set. Here, given two proper cones $C_1$, $C_2$, their minimal tensor product is the cone generated by products of the form $x_1 \otimes x_2$, where $x_1 \in C_1$ and $x_2 \in C_2$, while their maximal tensor product is the set of tensors that are positive under all product functionals $f_1 \otimes f_2$, where $f_1$ is positive on $C_1$ and $f_2$ is positive on $C_2$. Our proof techniques involve a mix of convex geometry, elementary algebraic topology, and computations inspired by quantum information theory. Our motivation comes from the foundations of physics: as an application, we show that any two non-classical systems modelled by general probabilistic theories can be entangled.


title: Time dependence of few-body F\"{o}rster interactions among ultracold Rydberg atoms

abstract: Rubidium Rydberg atoms in either $|m_j|$-sublevel of the $36p_{3/2}$ state can exchange energy via Stark-tuned F\"{o}rster resonances, including two-, three-, and four-body dipole-dipole interactions. Three-body interactions of this type were first reported and categorized by Faoro, \textit{et al.}~[Nat.\ Commun.\ \textbf{6}, 8173 (2015)] and their Borromean nature was confirmed by Tretyakov, \textit{et al.}~[Phys.\ Rev.\ Lett. \textbf{119}, 173402 (2017)]. We report the time dependence of the $N$-body F\"{o}rster resonance $N\times 36p_{3/2,|m_j|=1/2}\rightarrow 36s_{1/2}+37s_{1/2}+(N-2)\times 36p_{3/2,|m_j|=3/2}$, for $N=2,3$, and 4, by measuring the fraction of initially excited atoms that end up in the $37s_{1/2}$ state as a function of time. The essential features of these interactions are captured in an analytical model that includes only the many-body matrix elements and neighboring atom distribution. A more sophisticated simulation reveals the importance of beyond-nearest-neighbor interactions and of always-resonant interactions.


title: Adversarial Examples Improve Image Recognition

abstract: Adversarial examples are commonly viewed as a threat to ConvNets. Here we present an opposite perspective: adversarial examples can be used to improve image recognition models if harnessed in the right manner. We propose AdvProp, an enhanced adversarial training scheme which treats adversarial examples as additional examples, to prevent overfitting. Key to our method is the usage of a separate auxiliary batch norm for adversarial examples, as they have different underlying distributions to normal examples. We show that AdvProp improves a wide range of models on various image recognition tasks and performs better when the models are bigger. For instance, by applying AdvProp to the latest EfficientNet-B7 [28] on ImageNet, we achieve significant improvements on ImageNet (+0.7%), ImageNet-C (+6.5%), ImageNet-A (+7.0%), Stylized-ImageNet (+4.8%). With an enhanced EfficientNet-B8, our method achieves the state-of-the-art 85.5% ImageNet top-1 accuracy without extra data. This result even surpasses the best model in [20] which is trained with 3.5B Instagram images (~3000X more than ImageNet) and ~9.4X more parameters. Models are available at https://github.com/tensorflow/tpu/tree/master/models/official/efficientnet.


title: Entanglement islands in higher dimensions

abstract: It has been suggested in recent work that the Page curve of Hawking radiation can be recovered using computations in semi-classical gravity provided one allows for "islands" in the gravity region of quantum systems coupled to gravity. The explicit computations so far have been restricted to black holes in two-dimensional Jackiw-Teitelboim gravity. In this note, we numerically construct a five-dimensional asymptotically AdS geometry whose boundary realizes a four-dimensional Hartle-Hawking state on an eternal AdS black hole in equilibrium with a bath. We also numerically find two types of extremal surfaces: ones that correspond to having or not having an island. The version of the information paradox involving the eternal black hole exists in this setup, and it is avoided by the presence of islands. Thus, recent computations exhibiting islands in two-dimensional gravity generalize to higher dimensions as well.


title: First Resolved Scattered-Light Images of Four Debris Disks in Scorpius-Centaurus with the Gemini Planet Imager

abstract: We present the first spatially resolved scattered-light images of four debris disks around members of the Scorpius-Centaurus (Sco-Cen) OB Association with high-contrast imaging and polarimetry using the Gemini Planet Imager (GPI). All four disks are resolved for the first time in polarized light and one disk is also detected in total intensity. The three disks imaged around HD 111161, HD 143675, and HD 145560 are symmetric in both morphology and brightness distribution. The three systems span a range of inclinations and radial extents. The disk imaged around HD 98363 shows indications of asymmetries in morphology and brightness distribution, with some structural similarities to the HD 106906 planet-disk system. Uniquely, HD 98363 has a wide co-moving stellar companion Wray 15-788 with a recently resolved disk with very different morphological properties. HD 98363 A/B is the first binary debris disk system with two spatially resolved disks. All four targets have been observed with ALMA, and their continuum fluxes range from one non-detection to one of the brightest disks in the region. With the new results, a total of 15 A/F-stars in Sco-Cen have resolved scattered light debris disks, and approximately half of these systems exhibit some form of asymmetry. Combining the GPI disk structure results with information from the literature on millimeter fluxes and imaged planets reveals a diversity of disk properties in this young population. Overall, the four newly resolved disks contribute to the census of disk structures measured around A/F-stars at this important stage in the development of planetary systems.


title: Visualization by Example

abstract: While visualizations play a crucial role in gaining insights from data, generating useful visualizations from a complex dataset is far from an easy task. Besides understanding the functionality provided by existing visualization libraries, generating the desired visualization also requires reshaping and aggregating the underlying data as well as composing different visual elements to achieve the intended visual narrative. This paper aims to simplify visualization tasks by automatically synthesizing the required program from simple visual sketches provided by the user. Specifically, given an input data set and a visual sketch that demonstrates how to visualize a very small subset of this data, our technique automatically generates a program that can be used to visualize the entire data set. Automating visualization poses several challenges. First, because many visualization tasks require data wrangling in addition to generating plots, we need to decompose the end-to-end synthesis task into two separate sub-problems. Second, because the intermediate specification that results from the decomposition is necessarily imprecise, this makes the data wrangling task particularly challenging in our context. In this paper, we address these problems by developing a new compositional visualization-by-example technique that (a) decomposes the end-to-end task into two different synthesis problems over different DSLs and (b) leverages bi-directional program analysis to deal with the complexity that arises from having an imprecise intermediate specification. We implemented our visualization-by-example algorithm and evaluate it on 83 visualization tasks collected from on-line forums and tutorials. Viser can solve 84% of these benchmarks within a 600 second time limit, and, for those tasks that can be solved, the desired visualization is among the top-5 generated by Viser in 70% of the cases.


title: Regularizing Neural Networks by Stochastically Training Layer Ensembles

abstract: Dropout and similar stochastic neural network regularization methods are often interpreted as implicitly averaging over a large ensemble of models. We propose STE (stochastically trained ensemble) layers, which enhance the averaging properties of such methods by training an ensemble of weight matrices with stochastic regularization while explicitly averaging outputs. This provides stronger regularization with no additional computational cost at test time. We show consistent improvement on various image classification tasks using standard network topologies.


title: M51 ULX-7: super-orbital periodicity and constraints on the neutron star magnetic field

abstract: In the current work we explore the applicability of standard theoretical models of accretion to the observed properties of M51 ULX-7. The spin-up rate and observed X-ray luminosity are evidence of a neutron star with a surface magnetic field of $2-7\times10^{13}$ G, rotating near equilibrium. Analysis of the X-ray light-curve of the system (Swift/XRT data) reveals the presence of a $\sim$39 d super-orbital period. We argue that the super-orbital periodicity is due to disc precession, and that material is accreted onto the neutron star at a constant rate throughout it. Moreover, by attributing this modulation to the free precession of the neutron star we estimate a surface magnetic field strength of $3-4\times10^{13}$ G. The agreement of these two independent estimates provide strong constraints on the surface polar magnetic field strength of the neutron star.


title: LL/SC and Atomic Copy: Constant Time, Space Efficient Implementations using only pointer-width CAS

abstract: When designing concurrent algorithms, Load-Link/Store-Conditional (LL/SC) is often the ideal primitive to have because unlike Compare and Swap (CAS), LL/SC is immune to the ABA problem. However, the full semantics of LL/SC are not supported by any modern machine, so there has been a significant amount of work on simulations of LL/SC using Compare and Swap (CAS), a synchronization primitive that enjoys widespread hardware support. All of the algorithms so far that are constant time either use unbounded sequence numbers (and thus base objects of unbounded size), or require $\Omega(MP)$ space for $M$ LL/SC object (where $P$ is the number of processes). We present a constant time implementation of $M$ LL/SC objects using $\Theta(M+kP^2)$ space, where $k$ is the maximum number of overlapping LL/SC operations per process (usually a constant), and requiring only pointer-sized CAS objects. Our implementation can also be used to implement $L$-word $LL/SC$ objects in $\Theta(L)$ time (for both $LL$ and $SC$) and $\Theta((M+kP^2)L)$ space. To achieve these bounds, we begin by implementing a new primitive called Single-Writer Copy which takes a pointer to a word sized memory location and atomically copies its contents into another object. The restriction is that only one process is allowed to write/copy into the destination object at a time. We believe this primitive will be very useful in designing other concurrent algorithms as well.


title: Euclidean operator growth and quantum chaos

abstract: We consider growth of local operators under Euclidean time evolution in lattice systems with local interactions. We derive rigorous bounds on the operator norm growth and then proceed to establish an analog of the Lieb-Robinson bound for the spatial growth. In contrast to the Minkowski case when ballistic spreading of operators is universal, in the Euclidean case spatial growth is system-dependent and indicates if the system is integrable or chaotic. In the integrable case, the Euclidean spatial growth is at most polynomial. In the chaotic case, it is the fastest possible: exponential in 1D, while in higher dimensions and on Bethe lattices local operators can reach spatial infinity in finite Euclidean time. We use bounds on the Euclidean growth to establish constraints on individual matrix elements and operator power spectrum. We show that one-dimensional systems are special with the power spectrum always being superexponentially suppressed at large frequencies. Finally, we relate the bound on the Euclidean growth to the bound on the growth of Lanczos coefficients. To that end, we develop a path integral formalism for the weighted Dyck paths and evaluate it using saddle point approximation. Using a conjectural connection between the growth of the Lanczos coefficients and the Lyapunov exponent controlling the growth of OTOCs, we propose an improved bound on chaos valid at all temperatures.


title: Exploring Whether Super-Puffs Can Be Explained as Ringed Exoplanets

abstract: An intriguing, growing class of planets are the "super-puffs," objects with exceptionally large radii for their masses and thus correspondingly low densities ($\lesssim0.3\rm\,g\,cm^{-3}$). Here we consider whether they could have large inferred radii because they are in fact ringed. This would naturally explain why super-puffs have thus far only shown featureless transit spectra. We find that this hypothesis can work in some cases but not all. The close proximity of the super-puffs to their parent stars necessitates rings with a rocky rather than icy composition. This limits the radius of the rings, and makes it challenging to explain the large size of Kepler 51b, 51c, 51d, and 79d unless the rings are composed of porous material. Furthermore, the short tidal locking timescales for Kepler 18d, 223d, and 223e mean that these planets may be spinning too slowly, resulting in a small oblateness and rings that are warped by their parent star. Kepler 87c and 177c have the best chance of being explained by rings. Using transit simulations, we show that testing this hypothesis requires photometry with a precision of somewhere between ~10 ppm and ~50 ppm, which roughly scales with the ratio of the planet and star's radii. We conclude with a note about the recently discovered super-puff HIP 41378f.


title: The dynamics of stratified horizontal shear flows at low P\'eclet number

abstract: We consider the dynamics of a vertically stratified, horizontally-forced Kolmogorov flow. Motivated by astrophysical systems where the Prandtl number is often asymptotically small, our focus is the little-studied limit of high Reynolds number but low P\'eclet number (which is defined to be the product of the Reynolds number and the Prandtl number). Through a linear stability analysis, we demonstrate that the stability of two-dimensional modes to infinitesimal perturbations is independent of the stratification, whilst three-dimensional modes are always unstable in the limit of strong stratification and strong thermal diffusion. The subsequent nonlinear evolution and transition to turbulence is studied numerically using direct numerical simulations. For sufficiently large Reynolds numbers, four distinct dynamical regimes naturally emerge, depending upon the strength of the background stratification. By considering dominant balances in the governing equations, we derive scaling laws for each regime which explain the numerical data.


title: HST/COS Spectra of the Wind Lines of VFTS 102 and 285

abstract: Rapid rotation in massive stars imposes a latitudinal variation in the mass loss from radiatively driven winds that can lead to enhanced mass loss at the poles (with little angular momentum loss) and/or equator (with maximal angular momentum loss). Here we present an examination of the stellar wind lines of the two O-type stars with the fastest known equatorial velocities, VFTS 102 ($V\sin i = 610 \pm 30$ km/s; O9:Vnnne+) and VFTS 285 ($V\sin i = 609 \pm 29$ km/s; O7.5 Vnnn) in the Large Magellanic Cloud. Ultraviolet spectra of both stars were obtained with the Hubble Space Telescope Cosmic Origins Spectrograph. The spectrum of VFTS 285 displays a fast outflow in N V and a much slower wind in Si IV, and we argue that there is a two-wind regime in which mass loss is strong at the poles (fast and tenuous wind) but dominant at the equator (slow and dense winds). These ions and wind lines are not present in the spectrum of the cooler star VFTS 102, but the double-peaked H$\alpha$ emission in its spectrum implies equatorial mass loss into a circumstellar disk. The results suggest that in the fastest rotating O-stars, most mass is lost as an equatorial outflow, promoting angular momentum loss that contributes to a spin down over time.


title: Third-Person Visual Imitation Learning via Decoupled Hierarchical Controller

abstract: We study a generalized setup for learning from demonstration to build an agent that can manipulate novel objects in unseen scenarios by looking at only a single video of human demonstration from a third-person perspective. To accomplish this goal, our agent should not only learn to understand the intent of the demonstrated third-person video in its context but also perform the intended task in its environment configuration. Our central insight is to enforce this structure explicitly during learning by decoupling what to achieve (intended task) from how to perform it (controller). We propose a hierarchical setup where a high-level module learns to generate a series of first-person sub-goals conditioned on the third-person video demonstration, and a low-level controller predicts the actions to achieve those sub-goals. Our agent acts from raw image observations without any access to the full state information. We show results on a real robotic platform using Baxter for the manipulation tasks of pouring and placing objects in a box. Project video and code are at https://pathak22.github.io/hierarchical-imitation/


title: The EDGE-CALIFA Survey: Using Optical Extinction to Probe the Spatially-Resolved Distribution of Gas in Nearby Galaxies

abstract: We present an empirical relation between the cold gas surface density ($\Sigma_{\rm gas}$) and the optical extinction (${\rm A_V}$) in a sample of 103 galaxies from the Extragalactic Database for Galaxy Evolution (EDGE) survey. This survey provides CARMA interferometric CO observations for 126 galaxies included in the Calar Alto Legacy Integral Field Area (CALIFA) survey. The matched, spatially resolved nature of these data sets allows us to derive the $\Sigma_{\rm gas}$-${\rm A_V}$ relation on global, radial, and kpc (spaxel) scales. We determine ${\rm A_V}$ from the Balmer decrement (H$\alpha$/H$\beta$). We find that the best fit for this relation is $\Sigma_{\rm gas} ({\rm M_\odot pc^{-2}})\sim~26~\times~ {\rm A_V}({\rm mag})$, and that it does not depend on the spatial scale used for the fit. However, the scatter in the fits increases as we probe smaller spatial scales, reflecting the complex relative spatial distributions of stars, gas, and dust. We investigate the $\Sigma_{\rm gas}$/ ${\rm A_V}$ ratio on radial and spaxel scales as a function of ${\rm EW(H\alpha)}$. We find that at larger values of ${\rm EW(H\alpha)}$ (i.e., actively star-forming regions) this ratio tend to converge to the value expected for dust-star mixed geometries ($\sim$ 30 $\mathrm{M_{\odot} \,pc^{-2}\,mag^{-1}}$). On radial scales, we do not find a significant relation between the $\Sigma_{\rm gas}$/${\rm A_V}$ ratio and the ionized gas metallicity. We contrast our estimates of $\Sigma_{\rm gas}$ using ${\rm A_V}$ with compilations in the literature of the gas fraction on global and radial scales as well as with well known scaling relations such as the radial star-formation law and the $\Sigma_{\rm gas}$-$\Sigma_*$ relation. These tests show that optical extinction is a reliable proxy for estimating $\Sigma_{\rm gas}$ in the absence of direct sub/millimeter observations of the cold gas.


title: Intermediate-Mass Black Holes

abstract: We describe ongoing searches for intermediate-mass black holes with M_BH ~ 100-10^5 M_sun. We review a range of search mechanisms, both dynamical and those that rely on accretion signatures. We find that dynamical and accretion signatures alike point to a high fraction of 10^9-10^10 M_sun galaxies hosting black holes with M_BH<10^5 M_sun. In contrast, there are no solid detections of black holes in globular clusters. There are few observational constraints on black holes in any environment with M_BH ~ 100-10^4 M_sun. Considering low-mass galaxies with dynamical black hole masses and constraining limits, we find that the M_BH-sigma_* relation continues unbroken to M_BH~10^5 M_sun, albeit with large scatter. We believe the scatter is at least partially driven by a broad range in black hole mass, since the occupation fraction appears to be relatively high in these galaxies. We fold the observed scaling relations with our empirical limits on occupation fraction and the galaxy mass function to put observational bounds on the black hole mass function in galaxy nuclei. We are pessimistic that local demographic observations of galaxy nuclei alone could constrain seeding mechanisms, although either high-redshift luminosity functions or robust measurements of off-nuclear black holes could begin to discriminate models.


title: Two-dimensional partially ionized magnetohydrodynamic turbulence

abstract: Ionization occurs in the upper atmospheres of hot Jupiters and in the interiors of Gas Giant Planets, leading to magnetohydrodynamic (MHD) effects which couple the momentum and the magnetic field, thereby significantly altering the dynamics. In regions of moderate temperatures the gas is only partially ionized, which also leads to interactions with neutral molecules. To explore the turbulent dynamics of these regions we utilize Partially-Ionized MHD (PIMHD), a two-fluid model -- one neutral and one ionized -- coupled by a collision term proportional to the difference in velocities. Motivated by planetary settings where rotation constrains the large-scale motions to be mostly two-dimensional, we perform a suite of simulations to examine the parameter space of 2D PIMHD turbulence and pay particular attention to collisions and their role in the dynamics, dissipation, and energy exchange between the two species. We arrive at, and numerically confirm, an expression for the energy loss due to collisions in both the weakly and strongly collisional limits, and show that, in the latter limit, the neutral fluid couples to the ions and behaves as an MHD fluid. Finally, we discuss some implications of our findings to current understanding of gas giant planet atmospheres.


title: On comparison of estimators for proportional error nonlinear regression models in the limit of small measurement error

abstract: In this paper, we compare maximum likelihood (ML), quasi likelihood (QL) and weighted least squares (WLS) estimators for proportional error nonlinear regression models. Literature on thermoluminescence sedimentary dating revealed another estimator similar to weighted least squares but observed responses used as weights. This estimator that we refer to as data weighted least squares (DWLS) is also included in the comparison. We show that on the order $\sigma, $ all four estimators behave similar to ordinary least squares estimators for standard linear regression models. On the order of $\sigma^2, $ the estimators have biases. Formulae that are valid in the limit of small measurement error are derived for the biases and the variances of the four estimators. The maximum likelihood estimator has less bias compared to the quasi likelihood estimator. Conditions are derived under which weighted least squares and maximum likelihood estimators have similar biases. On the order of $\sigma^ {2} $, all estimators have similar standard errors. On higher order of $\sigma$, maximum likelihood estimator has smaller variance compared to quasi likelihood estimator, provided that the random errors have the same first four moments as the normal distribution. The maximum likelihood and quasi-likelihood estimating equations are unbiased. In large samples, these two estimators are distributed as multivariate normal. The estimating equations for weighted least squares and data weighted least squares are biased. However, in the limit of $\sigma \to 0$ and $n \to \infty, $ if $ n^ {1/2} \sigma$ remains bounded, these two estimators are also distributed as multivariate normal. A simulation study justified the applicability of the derived formulae in the presence of measurement errors typical in sedimentary data. Results are illustrated with a data set from thermoluminescence sedimentary dating.


title: The artefact of the Natural Resources Curse

abstract: This paper reexamines the validity of the natural resource curse hypothesis, using the database of mineral exporting countries. Our findings are as follows: (i) Resource-rich countries (RRCs) do not necessarily exhibit poor political, economic and social performance; (ii) RRCs that perform poorly have a low diversified exports portfolio; (iii) In contrast, RRCs with a low diversified exports portfolio do not necessarily perform poorly. Then, we develop a model of strategic interaction from a Bayesian game setup to study the role of leadership and governance in the management of natural resources. We show that an improvement in the leadership-governance binomial helps to discipline the behavior of lobby groups (theorem 1) and generate a Pareto improvement in the management of natural resources (theorem 2). Evidence from the World Bank Group's CPIA data confirms the later finding. Our results remain valid after some robustness checks.


title: Quantum Observables for continuous control of the Quantum Approximate Optimization Algorithm via Reinforcement Learning

abstract: We present a classical control mechanism for Quantum devices using Reinforcement Learning. Our strategy is applied to the Quantum Approximate Optimization Algorithm (QAOA) in order to optimize an objective function that encodes a solution to a hard combinatorial problem. This method provides optimal control of the Quantum device following a reformulation of QAOA as an environment where an autonomous classical agent interacts and performs actions to achieve higher rewards. This formulation allows a hybrid classical-Quantum device to train itself from previous executions using a continuous formulation of deep Q-learning to control the continuous degrees of freedom of QAOA. Our approach makes a selective use of Quantum measurements to complete the observations of the Quantum state available to the agent. We run tests of this approach on MAXCUT instances of size up to N = 21 obtaining optimal results. We show how this formulation can be used to transfer the knowledge from shorter training episodes to reach a regime of longer executions where QAOA delivers higher results.


title: In situ measurements of the variable slow solar wind near sector boundaries

abstract: The release of density structures at the tip of the coronal helmet streamers, likely as a consequence of magnetic reconnection, contributes to the mass flux of the slow solar wind. In situ measurements in the vicinity of the heliospheric plasma sheet of the magnetic field, protons and suprathermal electrons reveal details of the processes at play during the formation of density structures near the Sun. In a previous article, we exploited remote-sensing observations to derive a 3-D picture of the dynamic evolution of a streamer. We found evidence of the recurrent and continual release of dense blobs from the tip of the streamers. In the present paper, we interpret in situ measurements of the slow solar wind during solar maximum. Through both case and statistical analysis, we show that in situ signatures (magnetic field magnitude, smoothness and rotation, proton density and suprathermal electrons, in the first place) are consistent with the helmet streamers producing, in alternation, high-density regions (mostly disconnected) separated by magnetic flux ropes (mostly connected to the Sun). This sequence of emission of dense blobs and flux ropes also seems repeated at smaller scales inside each of the high-density regions. These properties are further confirmed with in situ measurements much closer to the Sun using Helios observations. We conclude on a model for the formation of dense blobs and flux ropes that explains both the in situ measurements and the remote-sensing observations presented in our previous studies.


title: How do atomic code uncertainties affect abundance measurements in the intracluster medium?

abstract: Accurate chemical abundance measurements of X-ray emitting atmospheres pervading massive galaxies, galaxy groups, and clusters provide essential information on the star formation and chemical enrichment histories of these large scale structures. Although the collisionally ionised nature of the intracluster medium (ICM) makes these abundance measurements relatively easy to derive, underlying spectral models can rely on different atomic codes, which brings additional uncertainties on the inferred abundances. Here, we provide a simple, yet comprehensive comparison between the codes SPEXACT v3.0.5 (cie model) and AtomDB v3.0.9 (vapec model) in the case of moderate, CCD-like resolution spectroscopy. We show that, in cool plasmas ($kT \lesssim 2$ keV), systematic differences up to $\sim$20% for the Fe abundance and $\sim$45% for the O/Fe, Mg/Fe, Si/Fe, and S/Fe ratios may still occur. Importantly, these discrepancies are also found to be instrument-dependent, at least for the absolute Fe abundance. Future improvements in these two codes will be necessary to better address questions on the ICM enrichment.


title: The ASAS-SN Catalog of Variable Stars VII: Contact Binaries are Different Above and Below the Kraft Break

abstract: We characterize ${\sim} 71,200$ W UMa type (EW) contact binaries, including ${\sim} 12,600$ new discoveries, using ASAS-SN $V$-band all-sky light curves along with archival data from Gaia, 2MASS, AllWISE, LAMOST, GALAH, RAVE, and APOGEE. There is a clean break in the EW period-luminosity relation at $\rm \log (\rm P/d){\simeq}-0.30$, separating the longer period early-type EW binaries from the shorter period, late-type systems. The two populations are even more cleanly separated in the space of period and effective temperature, by $\rm T_{eff}=6710\,K-1760\,K\,\log(P/0.5\,d)$. Early-type and late-type EW binaries follow opposite trends in $\rm T_{eff}$ with orbital period. For longer periods, early-type EW binaries are cooler, while late-type systems are hotter. We derive period-luminosity relationships (PLRs) in the $W_{JK}$, $V$, Gaia DR2 $G$, $J$, $H$, $K_s$ and $W_1$ bands for the late-type and early-type EW binaries separated both by period and effective temperature, and by period alone. The dichotomy of contact binaries is almost certainly related to the Kraft break and the related changes in envelope structure, winds and angular momentum loss.


title: Luminosity Models and Density Profiles for Nuclear Star Clusters for a Nearby Volume-Limited Sample of 29 Galaxies

abstract: Nuclear Star Clusters (NSCs) are dense clusters of stars that reside in the centers of a majority of the galaxies. In this paper, we study the density profiles for 29 galaxies in a volume-limited survey within 10 Mpc to characterize their NSCs. These galaxies span a 3$\times$10$^8$ - 8$\times$10$^{10}$ M$_\odot$ and a wide range of Hubble types. We use high-resolution Hubble space telescope archival data to create luminosity models for the galaxies using S\'ersic profiles to parameterize the NSCs. We also provide estimates for photometric masses of NSCs and their host galaxies using color-M/L relationships and examine their correlation. We use the Multi-Gaussian Expansion (MGE) to derive the NSC densities and their 3-D mass density profiles. The 3-D density profiles characterize the NSC densities on the scales as small as $\sim$1 pc, approaching the likely spheres of influence for black holes in these objects. We find that these densities correlate with galaxy mass, with galaxies most NSCs being denser than typical globular clusters. We parametrize the 3-D NSC density profiles, and their scatter and slope as a function of galaxy stellar mass to enable the construction of realistic nuclear mass profiles. These are useful in predicting the rate of tidal disruption events in galaxies. We will verify the results of this paper in a follow-up paper that presents the dynamical modeling of the same sample of NSCs.


title: Searching for further evidence for cloud-cloud collisions in L1188

abstract: In order to search for further observational evidence of cloud-cloud collisions in one of the promising candidates, L1188, we carried out observations of multiple molecular lines toward the intersection region of the two nearly orthogonal filamentary molecular clouds in L1188. Based on these observations, we find two parallel filamentary structures, both of which have at least two velocity components being connected with broad bridging features. We also found a spatially complementary distribution between the two molecular clouds, as well as enhanced $^{13}$CO emission and $^{12}$CO self-absorption toward their abutting regions. At the most blueshifted velocities, we unveil a 1~pc-long arc ubiquitously showing $^{12}$CO line wings. We discover two 22 GHz water masers, which are the first maser detections in L1188. An analysis of line ratios at a linear resolution of 0.2 pc suggests that L1188 is characterised by kinetic temperatures of 13--23~K and H$_{2}$ number densities of 10$^{3}$--10$^{3.6}$ cm$^{-3}$. On the basis of previous theoretical predictions and simulations, we suggest that these observational features can be naturally explained by the scenario of a cloud-cloud collision in L1188, although an additional contribution of stellar feedback from low-mass young stellar objects cannot be ruled out.


title: Coupled spin-charge dynamics in magnetic van der Waals heterostructures

abstract: We present a phenomenological theory for coupled spin-charge dynamics in magnetic van der Waals heterostructures. The system studied consists of a layered antiferromagnet inserted into a capacitive vdW heterostructure. It has been recently demonstrated that charge doping in such layered antiferromagnets can modulate the strength, and even the sign, of exchange coupling between the layer magnetizations. This provides a mechanism for electrically generating magnetization dynamics. The central result we predict here is that the magnetization dynamics reciprocally results in inducing charge dynamics. Such dynamics makes magnetic van der Waals heterostructures interesting candidates for spintronics applications. To this end, we also show that these systems can be used to convert sub-THz radiation induced magnetization dynamics into electrical signals.


title: On the Gauge Invariance of Cosmological Gravitational Waves

abstract: The issue of the gauge invariance of gravitational waves arises if they are produced in the early universe at second-order in perturbation theory. We address it by dividing the discussion about the gauge invariance in three parts: the production of gravitational waves, their propagation in the real universe, and their measurement.


title: Canonical analysis of Holst action without second-class constraints

abstract: We perform the canonical analysis of the Holst action for general relativity with a cosmological constant without introducing second-class constraints. Our approach consists in identifying the dynamical and nondynamical parts of the involved variables from the very outset. After integrating out the nondynamical variables associated with the connection, we obtain the description of phase space in terms of manifestly $SO(3,1)$ [or $SO(4)$, depending on the signature] covariant canonical variables and first-class constraints only. We impose the time gauge on them and show that the Ashtekar-Barbero formulation of general relativity emerges. Later, we discuss a family of canonical transformations that allows us to construct new $SO(3,1)$ [or $SO(4)$] covariant canonical variables for the phase space of the theory and compare them with the ones already reported in the literature, pointing out the presence of a set of canonical variables not considered before. Finally, we resort to the time gauge again and find that the theory, when written in terms of the new canonical variables, either collapses to the $SO(3)$ ADM formalism or to the Ashtekar-Barbero formalism with a rescaled Immirzi parameter.


title: Modeling the magnetized Local Bubble from dust data

abstract: The Sun is embedded in the so-called Local Bubble (LB) -- a cavity of hot plasma created by supernova explosions and surrounded by a shell of cold, dusty gas. Knowing the local distortion of the Galactic magnetic field associated with the LB is critical for the modeling of interstellar polarization data at high Galactic latitudes. In this his paper, we relate the structure of the Galactic magnetic field on the LB scale to three-dimensional (3D) maps of the local interstellar medium (ISM). First, we extracted the geometry of the LB shell, its inner surface, in particular from 3D dust extinction maps of the local ISM. We expanded the shell inner surface in spherical harmonics, up to a variable maximum multipole degree, which enabled us to control the level of complexity for the modeled surface. Next, we applied an analytical model for the ordered magnetic field in the shell to the modeled shell surface. This magnetic field model was successfully fitted to the \textit{Planck} 353~GHz dust polarized emission maps over the Galactic polar caps. For each polar cap, the direction of the mean magnetic field derived from dust polarization (together with the prior that the field points toward longitude $90^\circ \pm 90^\circ$) is found to be consistent with the Faraday spectra of the nearby diffuse synchrotron emission. Our work presents a new approach to modeling the local structure of the Galactic magnetic field. We expect our methodology and our results to be useful both in modeling the local ISM as traced by its different components and in modeling the dust polarized emission, which is a long-awaited input for studies of the polarized foregrounds for cosmic microwave background.


title: Spiral arms and instability within the AFGL 4176 mm1 disc

abstract: We present high-resolution (30 mas or 130 au at 4.2 kpc) Atacama Large Millimeter/submillimeter Array observations at 1.2 mm of the disc around the forming O-type star AFGL 4176 mm1. The disc (AFGL 4176 mm1-main) has a radius of ~1000 au and contains significant structure, most notably a spiral arm on its redshifted side. We fitted the observed spiral with logarithmic and Archimedean spiral models. We find that both models can describe its structure, but the Archimedean spiral with a varying pitch angle fits its morphology marginally better. As well as signatures of rotation across the disc, we observe gas arcs in CH$_3$CN that connect to other millimetre continuum sources in the field, supporting the picture of interactions within a small cluster around AFGL 4176 mm1-main. Using local thermodynamic equilibrium modelling of the CH$_3$CN K-ladder, we determine the temperature and velocity field across the disc, and thus produce a map of the Toomre stability parameter. Our results indicate that the outer disc is gravitationally unstable and has already fragmented or is likely to fragment in the future, possibly producing further companions. These observations provide evidence that disc fragmentation is one possible pathway towards explaining the high fraction of multiple systems around high-mass stars.


title: Two-dimensional quantum-link lattice Quantum Electrodynamics at finite density

abstract: We present an unconstrained tree tensor network approach to the study of lattice gauge theories in two spatial dimensions showing how to perform numerical simulations of theories in presence of fermionic matter and four-body magnetic terms, at zero and finite density, with periodic and open boundary conditions. We exploit the quantum link representation of the gauge fields and demonstrate that a fermionic rishon representation of the quantum links allows us to efficiently handle the fermionic matter while finite densities are naturally enclosed in the tensor network description. We explicit perform calculations for quantum electrodynamics in the spin-one quantum link representation on lattice sizes of up to 16x16 sites, detecting and characterizing different quantum regimes. In particular, at finite density, we detect signatures of a phase separation as a function of the bare mass values at different filling densities. The presented approach can be extended straightforwardly to three spatial dimensions.


title: Neutral Hydrogen Observations of Low Surface Brightness Galaxies around M101 and NGC 5485

abstract: We present atomic hydrogen (HI) observations using the Robert C. Byrd Green Bank Telescope along the lines-of-sight to 27 low surface brightness (LSB) dwarf galaxy candidates discovered in optical searches around M101. We detect HI reservoirs in 5 targets and place stringent upper limits on the remaining 22, implying that they are gas poor. The distances to our HI detections range from 7 Mpc --150 Mpc, demonstrating the utility of wide-bandpass HI observations as a follow-up tool. The systemic velocities of 3 detections are consistent with that of the NGC~5485 group behind M101, and we suggest that our 15 non-detections with lower distance limits from the optical are associated with and have been stripped by that group. We find that the gas richnesses of confirmed M101 satellites are broadly consistent with those of the Milky Way satellites, as well as with those of satellites around other hosts of comparable mass, when survey completeness is taken into account. This suggests that satellite quenching and gas stripping proceeds similarly around halos of similar mass, in line with theoretical expectations.


title: The Massive and Distant Clusters of WISE Survey VII: The Environments and Properties of Radio Galaxies in Clusters at z~1

abstract: We present the results from a study with NSF's Karl G. Jansky Very Large Array (VLA) to determine the radio morphologies of extended radio sources and the properties of their host galaxies in 50 massive galaxy clusters at z~1. We find a majority of the radio morphologies to be Fanaroff-Riley (FR) type IIs. By analyzing the infrared counterparts of the radio sources, we find that ~40% of the host galaxies are the candidate brightest cluster galaxy (BCG) and ~83% are consistent with being one of the top six most massive galaxies in the cluster. We investigate the role of environmental factors on the radio-loud AGN population by examining correlations between environmental and radio-galaxy properties. We find that the highest stellar mass hosts ($M_{*} \gtrsim$ 4$\times 10^{11} M_{\odot}$) are confined to the cluster center and host compact jets. There is evidence for an increase in the size of the jets with cluster-centric radius, which may be attributed to the decreased ICM pressure confinement with increasing radius. Besides this correlation, there are no other significant correlations between the properties of the radio-AGN (luminosity, morphology, or size) and environmental properties (cluster richness and location within the cluster). The fact that there are more AGN in the cluster environment than the field at this epoch, combined with the lack of strong correlation between galaxy and environmental properties, argues that the cluster environment fosters radio activity but does not solely drive the evolution of these sources at this redshift.


title: Generalized probability rules from a timeless formulation of Wigner's friend scenarios

abstract: The quantum measurement problem can be regarded as the tension between the two alternative dynamics prescribed by quantum mechanics: the unitary evolution of the wave function and the state-update rule (or "collapse") at the instant a measurement takes place. The notorious Wigner's friend gedankenexperiment constitutes the paradoxical scenario in which different observers (one of whom is observed by the other) describe one and the same interaction differently, one --the Friend-- via state-update and the other --Wigner-- unitarily. This can lead to Wigner and his friend assigning different probabilities to the outcome of the same subsequent measurement. In this paper, we apply the Page-Wootters mechanism (PWM) as a timeless description of Wigner's friend-like scenarios. We show that the standard rules to assign two-time conditional probabilities within the PWM need to be modified to deal with the Wigner's friend gedankenexperiment. We identify three main definitions of such modified rules to assign two-time conditional probabilities, all of which reduce to standard quantum theory for non-Wigner's friend scenarios. However, when applied to the Wigner's friend setup each rule assigns different conditional probabilities, potentially resolving the probability-assignment paradox in a different manner. Moreover, one rule imposes strict limits on when a joint probability distribution for the measurement outcomes of Wigner and his Friend is well-defined, which single out those cases where Wigner's measurement does not disturb the Friend's memory and such a probability has an operational meaning in terms of collectible statistics. Interestingly, the same limits guarantee that said measurement outcomes fulfill the consistency condition of the consistent histories framework.


title: Heterotic strings on $(K3\times T^2)/\mathbb{Z}_3$ and their dual Calabi-Yau threefolds

abstract: In this paper we study compactifications of the ${\cal N}=2$ heterotic $E_8\times E_8$ string on $(K3\times T^2)/\mathbb{Z}_3$ with various gauge backgrounds and calculate the topological couplings in the effective supergravity action that arise from one-loop amplitudes. We then identify candidates for dual type IIA compactifications on Calabi-Yau threefolds and compare the heterotic results with the corresponding topological string amplitudes. We find that the dual Calabi-Yau geometries are $K3$ fibrations that are also genus one fibered with three-sections. Moreover, we show that the intersection form on the polarization lattice of the $K3$ fibration has to be three times the intersection form on the Narain lattice $\Gamma^{1,1}$.


title: Parallelising MCMC via Random Forests

abstract: For Bayesian computation in big data contexts, the divide-and-conquer MCMC concept splits the whole data set into batches, runs MCMC algorithms separately over each batch to produce samples of parameters, and combines them to produce an approximation of the target distribution. In this article, we embed random forests into this framework and use each subposterior/partial-posterior as a proposal distribution to implement importance sampling. Unlike the existing divide-and-conquer MCMC, our methods are based on scaled subposteriors, whose scale factors are not necessarily restricted to being equal to one or to the number of subsets. Through several experiments, we show that our methods work well with models ranging from Gaussian cases to strongly non-Gaussian cases, and include model misspecification.


title: Topic Model for four-top at the LHC

abstract: We study the implementation of a Topic Model algorithm in four-top searches at the LHC as a test-probe of a not ideal system for applying this technique. We study this Topic Model behavior as its different hypotheses such as mutual reducibility and equal distribution in all samples shift from true. The four-top final state at the LHC is not only relevant because it does not fulfill these conditions, but also because it is a difficult and inefficient system to reconstruct and current Monte Carlo modeling of signal and backgrounds suffers from non-negligible uncertainties. We implement this Topic Model algorithm in the Same-Sign lepton channel where S/B is of order one and all backgrounds cannot have more than two b-jets at parton level. We define different mixtures according to the number of b-jets and we use the total number of jets to demix. Since only the background has an anchor bin, we find that we can reconstruct the background in the signal region independently of Monte Carlo. We propose to use this information to tune the Monte Carlo in the signal region and then compare signal prediction with data. We also explore Machine Learning techniques applied to this Topic Model algorithm and find slight improvements as well as potential roads to investigate. Although our findings indicate that still with the full LHC run 3 data the implementation would be challenging, we pursue through this work to find ways to reduce the impact of Monte Carlo simulations in four-top searches at the LHC.


title: Algebraic solution to constrained bi-criteria decision problem of rating alternatives through pairwise comparisons

abstract: We consider a decision-making problem to evaluate absolute ratings of alternatives from the results of their pairwise comparisons according to two criteria, subject to constraints on the ratings. We formulate the problem as a bi-objective optimization problem of constrained matrix approximation in the Chebyshev sense in logarithmic scale. The problem is to approximate the pairwise comparison matrices for each criterion simultaneously by a common consistent matrix of unit rank, which determines the vector of ratings. We represent and solve the optimization problem in the framework of tropical (idempotent) algebra, which deals with the theory and applications of idempotent semirings and semifields. The solution involves the introduction of two parameters that represent the minimum values of approximation error for each matrix and thereby describe the Pareto frontier for the bi-objective problem. The optimization problem then reduces to a parametrized vector inequality. The necessary and sufficient conditions for solutions of the inequality serve to derive the Pareto frontier for the problem. All solutions of the inequality, which correspond to the Pareto frontier, are taken as a complete Pareto-optimal solution to the problem. We apply these results to the decision problem of interest and present illustrative examples.


title: Correlating photons using the collective nonlinear response of atoms weakly coupled to an optical mode

abstract: Photons in a nonlinear medium can repel or attract each other, resulting in a strongly correlated quantum many-body system. Typically, such strongly correlated states of light arise from the extreme nonlinearity granted by quantum emitters that are strongly coupled to a photonic mode. However, in these approaches, unavoidable dissipation, like photon loss, blurs nonlinear quantum effects. Here, we generate strongly correlated photon states using only weak coupling and taking advantage of dissipation. We launch light through an ensemble of non-interacting waveguide-coupled atoms, which induce correlations between simultaneously arriving photons through collectively enhanced nonlinear interactions. These correlated photons then experience less dissipation than the uncorrelated ones. Depending on the number of atoms, we experimentally observe strong photon bunching or anti-bunching of the transmitted light. This realization of a collectively enhanced nonlinearity may turn out transformational for quantum information science and opens new avenues for generating nonclassical light, covering frequencies from the microwave to the X-ray regime.


title: Mechanisms of Spatiotemporal Mode-Locking

abstract: Mode-locking is a process in which different modes of an optical resonator establish, through nonlinear interactions, stable synchronization. This self-organization underlies light sources that enable many modern scientific applications, such as ultrafast and high-field optics and frequency combs. Despite this, mode-locking has almost exclusively referred to self-organization of light in a single dimension - time. Here we present a theoretical approach, attractor dissection, for understanding three-dimensional (3D) spatiotemporal mode-locking (STML). The key idea is to find, for each distinct type of 3D pulse, a specific, minimal reduced model, and thus to identify the important intracavity effects responsible for its formation and stability. An intuition for the results follows from the 'minimum loss principle,' the idea that a laser strives to find the configuration of intracavity light that minimizes loss (maximizes gain extraction). Through this approach, we identify and explain several distinct forms of STML. These novel phases of coherent laser light have no analogues in 1D and are supported by experimental measurements of the three-dimensional field, revealing STML states comprising more than $10^7$ cavity modes. Our results should facilitate the discovery and understanding of new higher-dimensional forms of coherent light which, in turn, may enable new applications.


title: Maurice Janet's algorithms on systems of linear partial differential equations

abstract: This article presents the emergence of formal methods in theory of partial differential equations (PDE) in the french school of mathematics through Janet's work in the period 1913-1930. In his thesis and in a series of articles published during this period, M. Janet introduced an original formal approach to deal with the solvability of the problem of initial conditions for finite linear PDE systems. His constructions implicitly used an interpretation of a monomial PDE system as a generating family of a multiplicative set of monomials. He introduced an algorithmic method on multiplicative sets to compute compatibility conditions, and to study the problem of the existence and the unicity of a solution to a linear PDE system with given initial conditions. The compatibility conditions are formulated using a refinement of the division operation on monomials defined with respect to a partition of the set of variables into multiplicative and non-multiplicative variables. M. Janet was a pioneer in the development of these algorithmic methods, and the completion procedure that he introduced on polynomials was the first one in a long and rich series of works on completion methods which appeared independently throughout the 20th century in various algebraic contexts.


title: A Conceptual Framework for Lifelong Learning

abstract: Humans can learn a variety of concepts and skills incrementally over the course of their lives while exhibiting many desirable properties, such as continual learning without forgetting, forward transfer and backward transfer of knowledge, and learning a new concept or task with only a few examples. Several lines of machine learning research, such as lifelong learning, few-shot learning, and transfer learning, attempt to capture these properties. However, most previous approaches can only demonstrate subsets of these properties, often by different complex mechanisms. In this work, we propose a simple yet powerful unified framework that supports almost all of these properties and approaches through one central mechanism. We also draw connections between many peculiarities of human learning (such as memory loss and "rain man") and our framework. While we do not present any state-of-the-art results, we hope that this conceptual framework provides a novel perspective on existing work and proposes many new research directions.


title: NGTS clusters survey. I. Rotation in the young benchmark open cluster Blanco 1

abstract: We determine rotation periods for 127 stars in the ~115 Myr old Blanco 1 open cluster using ~200 days of photometric monitoring with the Next Generation Transit Survey (NGTS). These stars span F5-M3 spectral types (1.2 $\gtrsim M \gtrsim$ 0.3 M$_{\odot}$) and increase the number of known rotation periods in Blanco 1 by a factor of four. We determine rotation periods using three methods: Gaussian process (GP) regression, generalised autocorrelation (G-ACF) and Lomb-Scargle (LS) periodograms, and find that GPs and G-ACF are more applicable to evolving spot modulation patterns. Between mid-F and mid-K spectral types, single stars follow a well-defined rotation sequence from ~2 to 10 days, whereas stars in photometric multiple systems typically rotate faster. This may suggest that the presence of a moderate-to-high mass ratio companion inhibits angular momentum loss mechanisms during the early pre-main sequence, and this signature has not been erased at ~100 Myr. The majority of mid-F to mid-K stars display evolving modulation patterns, whereas most M stars show stable modulation signals. This morphological change coincides with the shift from a well-defined rotation sequence (mid-F to mid-K stars) to a broad rotation period distribution (late-K and M stars). Finally, we compare our rotation results for Blanco 1 to the similarly-aged Pleiades: the single star populations in both clusters possess consistent rotation period distributions, which suggests that the angular momentum evolution of stars follows a well-defined pathway that is, at least for mid-F to mid-K stars, strongly imprinted by ~100 Myr.


title: Density and temperature of cosmic-web filaments on scales of tens of megaparsecs

abstract: We studied physical properties of matter in 24,544 filaments ranging from 30 to 100 Mpc in length, identified in the Sloan Digital Sky Survey (SDSS). We stacked the Comptonization y map produced by the Planck Collaboration around the filaments, excluding the resolved galaxy groups and clusters above a mass of ~3*10^13 Msun. We detected the thermal Sunyaev-Zel'dovich signal for the first time at a significance of 4.4 sigma in filamentary structures on such a large scale. We also stacked the Planck cosmic microwave background (CMB) lensing convergence map in the same manner and detected the lensing signal at a significance of 8.1 sigma. To estimate physical properties of the matter, we considered an isothermal cylindrical filament model with a density distribution following a beta-model (beta=2/3). Assuming that the gas distribution follows the dark matter distribution, we estimate that the central gas and matter overdensity and gas temperature are overdensity = (19.0 +27.3 -12.1) and temperature = (1.2 +- 0.4)*10^6 K, which results in a measured baryon fraction of (0.080 +0.116 -0.051) * Omega_b.


title: High Resolution Photoexcitation Measurements Exacerbate the Long-Standing Fe XVII Oscillator Strength Problem

abstract: For more than 40 years, most astrophysical observations and laboratory studies of two key soft x-ray diagnostic $2p-3d$ transitions, $3C$ and $3D$, in Fe XVII ions found oscillator strength ratios $f(3C)/f(3D)$ disagreeing with theory, but uncertainties had precluded definitive statements on this much studied conundrum. Here, we resonantly excite these lines using synchrotron radiation at PETRA III, and reach, at a millionfold lower photon intensities, a 10 times higher spectral resolution, and 3 times smaller uncertainty than earlier work. Our final result of $f(3C)/f(3D) = 3.09(8)(6)$ supports many of the earlier clean astrophysical and laboratory observations, while departing by five sigmas from our own newest large-scale ab initio calculations, and excluding all proposed explanations, including those invoking nonlinear effects and population transfers.

categories: physics.atom-ph astro-ph.HE astro-ph.IM astro-ph.SR physics.plasm-ph


title: On the number of vertices of Newton--Okounkov polygons

abstract: The Newton--Okounkov body of a big divisor D on a smooth surface is a numerical invariant in the form of a convex polygon. We study the geometric significance of the shape of Newton--Okounkov polygons of ample divisors, showing that they share several important properties of Newton polygons on toric surfaces. In concrete terms, sides of the polygon are associated to some particular irreducible curves, and their lengths are determined by the intersection numbers of these curves with D. As a consequence of our description we determine the numbers k such that D admits some k-gon as a Newton--Okounkov body, elucidating the relationship of these numbers with the Picard number of the surface, which was first hinted at by work of K\"uronya, Lozovanu and Maclean.


title: Automatically Neutralizing Subjective Bias in Text

abstract: Texts like news, encyclopedias, and some social media strive for objectivity. Yet bias in the form of inappropriate subjectivity - introducing attitudes via framing, presupposing truth, and casting doubt - remains ubiquitous. This kind of bias erodes our collective trust and fuels social conflict. To address this issue, we introduce a novel testbed for natural language generation: automatically bringing inappropriately subjective text into a neutral point of view ("neutralizing" biased text). We also offer the first parallel corpus of biased language. The corpus contains 180,000 sentence pairs and originates from Wikipedia edits that removed various framings, presuppositions, and attitudes from biased sentences. Last, we propose two strong encoder-decoder baselines for the task. A straightforward yet opaque CONCURRENT system uses a BERT encoder to identify subjective words as part of the generation process. An interpretable and controllable MODULAR algorithm separates these steps, using (1) a BERT-based classifier to identify problematic words and (2) a novel join embedding through which the classifier can edit the hidden states of the encoder. Large-scale human evaluation across four domains (encyclopedias, news headlines, books, and political speeches) suggests that these algorithms are a first step towards the automatic identification and reduction of bias.


title: Multi-scale observations of thermal non-equilibrium cycles in coronal loops

abstract: Thermal non-equilibrium (TNE) is a phenomenon that can occur in solar coronal loops when the heating is quasi-constant and highly-stratified. Under such heating conditions, coronal loops undergo cycles of evaporation and condensation. The recent observations of ubiquitous long-period intensity pulsations in coronal loops and their relationship with coronal rain have demonstrated that understanding the characteristics of TNE cycles is an essential step in constraining the circulation of mass and energy in the corona. We report unique observations with the Solar Dynamics Observatory (SDO) and the Swedish 1-m Solar Telescope (SST) that link the captured thermal properties across the extreme spatiotemporal scales covered by TNE processes. Within the same coronal loop bundle, we captured 6 hr period coronal intensity pulsations in SDO/AIA and coronal rain observed off-limb in the chromospheric Halpha and Ca II K spectral lines with SST/CRISP and SST/CHROMIS. We combined a multi-thermal analysis of the cycles with AIA and an extensive spectral characterisation of the rain clumps with the SST. We find clear evidence of evaporation-condensation cycles in the corona which are linked with periodic coronal rain showers. The high-resolution spectroscopic instruments at the SST reveal the fine-structured rain strands and allow us to probe the cooling phase of one of the cycles down to chromospheric temperatures. These observations reinforce the link between long-period intensity pulsations and coronal rain. They also demonstrate the capability of TNE to shape the dynamics of active regions on the large scales as well as on the smallest scales currently resolvable.


title: Superconducting Contact and Quantum Interference Between Two-Dimensional van der Waals and Three-Dimensional Conventional Superconductors

abstract: Two-dimensional (2D) transition-metal dichalcogenide superconductors have unique and desirable properties for integration with conventional superconducting circuits. However, fully superconducting contact must be made between the 2D material and three-dimensional (3D) superconductors in order to employ the standard microwave drive and readout of qubits in such circuits. Here, we present a method for creating zero-resistance contacts between 2D NbSe$_2$ and 3D aluminum that behave as Josephson junctions (JJs) with large effective areas compared to 3D-3D JJs. We present a model for the supercurrent flow in a 2D-3D superconducting structure by numerical solution of the Ginzburg-Landau equations and find good agreement with experiment. These results demonstrate a crucial step towards a new generation of hybrid superconducting quantum circuits.


title: RefinedMPL: Refined Monocular PseudoLiDAR for 3D Object Detection in Autonomous Driving

abstract: In this paper, we strive for solving the ambiguities arisen by the astoundingly high density of raw PseudoLiDAR for monocular 3D object detection for autonomous driving. Without much computational overhead, we propose a supervised and an unsupervised sparsification scheme of PseudoLiDAR prior to 3D detection. Both the strategies assist the standard 3D detector gain better performance over the raw PseudoLiDAR baseline using only ~5% of its points on the KITTI object detection benchmark, thus making our monocular framework and LiDAR-based counterparts computationally equivalent (Figure 1). Moreover, our architecture agnostic refinements provide state-of-the-art results on KITTI3D test set for "Car" and "Pedestrian" categories with 54% relative improvement for "Pedestrian". Finally, exploratory analysis is performed on the discrepancy between monocular and LiDAR-based 3D detection frameworks to guide future endeavours.


title: Using Socially Expressive Mixed Reality Arms for Enhancing Low-Expressivity Robots

abstract: Expressivity--the use of multiple modalities to convey internal state and intent of a robot--is critical for interaction. Yet, due to cost, safety, and other constraints, many robots lack high degrees of physical expressivity. This paper explores using mixed reality to enhance a robot with limited expressivity by adding virtual arms that extend the robot's expressiveness. The arms, capable of a range of non-physically-constrained gestures, were evaluated in a between-subject study ($n=34$) where participants engaged in a mixed reality mathematics task with a socially assistive robot. The study results indicate that the virtual arms added a higher degree of perceived emotion, helpfulness, and physical presence to the robot. Users who reported a higher perceived physical presence also found the robot to have a higher degree of social presence, ease of use, usefulness, and had a positive attitude toward using the robot with mixed reality. The results also demonstrate the users' ability to distinguish the virtual gestures' valence and intent.


title: Statistical Guarantees for Local Spectral Clustering on Random Neighborhood Graphs

abstract: We study the Personalized PageRank (PPR) algorithm, a local spectral method for clustering, which extracts clusters using locally-biased random walks around a given seed node. In contrast to previous work, we adopt a classical statistical learning setup, where we obtain samples from an unknown nonparametric distribution, and aim to identify sufficiently salient clusters. We introduce a trio of population-level functionals -- the normalized cut, conductance, and local spread, analogous to graph-based functionals of the same name -- and prove that PPR, run on a neighborhood graph, recovers clusters with small population normalized cut and large conductance and local spread. We apply our general theory to establish that PPR identifies connected regions of high density (density clusters) that satisfy a set of natural geometric conditions. We also show a converse result, that PPR can fail to recover geometrically poorly-conditioned density clusters, even asymptotically. Finally, we provide empirical support for our theory.


title: Efficient Drone Mobility Support Using Reinforcement Learning

abstract: Flying drones can be used in a wide range of applications and services from surveillance to package delivery. To ensure robust control and safety of drone operations, cellular networks need to provide reliable wireless connectivity to drone user equipments (UEs). To date, existing mobile networks have been primarily designed and optimized for serving ground UEs, thus making the mobility support in the sky challenging. In this paper, a novel handover (HO) mechanism is developed for a cellular-connected drone system to ensure robust wireless connectivity and mobility support for drone-UEs. By leveraging tools from reinforcement learning, HO decisions are dynamically optimized using a Q-learning algorithm to provide an efficient mobility support in the sky. The results show that the proposed approach can significantly reduce (e.g., by 80%) the number of HOs, while maintaining connectivity, compared to the baseline HO scheme in which the drone always connects to the strongest cell.


title: Too Quiet in the Library: An Empirical Study of Security Updates in Android Apps' Native Code

abstract: Android apps include third-party native libraries to increase performance and to reuse functionality. Native code is directly executed from apps through the Java Native Interface or the Android Native Development Kit. Android developers add precompiled native libraries to their projects, enabling their use. Unfortunately, developers often struggle or simply neglect to update these libraries in a timely manner. This results in the continuous use of outdated native libraries with unpatched security vulnerabilities years after patches became available. To further understand such phenomena, we study the security updates in native libraries in the most popular 200 free apps on Google Play from Sept. 2013 to May 2020. A core difficulty we face in this study is the identification of libraries and their versions. Developers often rename or modify libraries, making their identification challenging. We create an approach called LibRARIAN (LibRAry veRsion IdentificAtioN) that accurately identifies native libraries and their versions as found in Android apps based on our novel similarity metric bin2sim. LibRARIAN leverages different features extracted from libraries based on their metadata and identifying strings in read-only sections. We discovered 53/200 popular apps (26.5%) with vulnerable versions with known CVEs between Sept. 2013 and May 2020, with 14 of those apps remaining vulnerable. We find that app developers took, on average, 528.71 days to apply security patches, while library developers release a security patch after 54.59 days - a 10 times slower rate of update.


title: Escalating core formation with dark matter self-heating

abstract: Exothermic scatterings of dark matter (DM) produce DM particles with significant kick velocities inside DM halos. In collaboration with DM self-interaction, the excess kinetic energy of the produced DM particles is distributed to the others, which self-heats the DM particles as a whole. The DM self-heating is efficient towards the halo center, and the heat injection is used to enhance the formation of a uniform density core inside halos. The effect of DM self-heating is expected to be more significant in smaller halos for two reasons: 1) the exothermic cross section times the relative velocity, $\left\langle\sigma_{\rm exo}v_{\rm rel}\right\rangle$, is constant; 2) and the ratio of the injected heat to the velocity dispersion squared gets larger towards smaller halos. For the first time, we quantitatively investigate the core formation from DM self-heating for halos in a wide mass range ($10^{9}$-$10^{15}\,{\rm M}_\odot$) using the gravothermal fluid formalism. Notably, we demonstrate that the core formation is sharply escalating towards smaller halos by taking the self-heating DM (i.e., DM that semi-annihilates and self-interacts) as an example. We show that the sharp escalation of core formation may cause a tension in simultaneously explaining the observed central mass deficit of Milky Way satellites, and field dwarf/low surface brightness spiral galaxies. While the details of the self-heating effect may differ among models, we expect that the sharp halo-mass dependence of the core formation is a general feature of exothermic DM.


title: Harmonic analysis on the rank-$2$ value group of a two-dimensional local field

abstract: In this work we construct harmonic analysis on free Abelian groups of rank $2$, namely: we construct and investigate spaces of functions and distributions, Fourier transforms, actions of discrete and extended discrete Heisenberg groups. In case of the rank-$2$ value group of a two-dimensional local field with finite last residue field we connect this harmonic analysis with harmonic analysis on the two-dimensional local field, where the latter harmonic analysis was constructed in earlier works by the authors.


title: Low frequency view of GRB 190114C reveals time varying shock micro-physics

abstract: We present radio and optical afterglow observations of the TeV-bright long Gamma Ray Burst (GRB) 190114C at a redshift of $z=0.425$, which was detected by the MAGIC telescope. Our observations with ALMA, ATCA, and uGMRT were obtained by our low frequency observing campaign and range from $\sim1$ to $\sim140$ days after the burst and the optical observations were done with three optical telescopes spanning up to $\sim25$ days after the burst. Long term radio/mm observations reveal the complex nature of the afterglow, which does not follow the spectral and temporal closure relations expected from the standard afterglow model. We find that the microphysical parameters of the external forward shock, representing the share of shock-created energy in the non-thermal electron population and magnetic field, are evolving with time. The inferred kinetic energy in the blast-wave depends strongly on the assumed ambient medium density profile, with a constant density medium demanding almost an order of magnitude higher energy than in the prompt emission, while a stellar wind-driven medium requires approximately the same amount energy as in prompt emission.


title: Universal structure of dark matter haloes over a mass range of 20 orders of magnitude

abstract: Cosmological models in which dark matter consists of cold elementary particles predict that the dark halo population should extend to masses many orders of magnitude below those at which galaxies can form. Here we report a cosmological simulation of the formation of present-day haloes over the full range of observed halo masses (20 orders of magnitude) when dark matter is assumed to be in the form of weakly interacting massive particles of mass approximately 100 gigaelectronvolts. The simulation has a full dynamic range of 30 orders of magnitude in mass and resolves the internal structure of hundreds of Earth-mass haloes in as much detail as it does for hundreds of rich galaxy clusters. We find that halo density profiles are universal over the entire mass range and are well described by simple two-parameter fitting formulae. Halo mass and concentration are tightly related in a way that depends on cosmology and on the nature of the dark matter. For a fixed mass, the concentration is independent of the local environment for haloes less massive than those of typical galaxies. Haloes over the mass range of 10^3 to 10^11 solar masses contribute about equally (per logarithmic interval) to the luminosity produced by dark matter annihilation, which we find to be smaller than all previous estimates by factors ranging up to one thousand.


title: Communication-Efficient and Byzantine-Robust Distributed Learning with Error Feedback

abstract: We develop a communication-efficient distributed learning algorithm that is robust against Byzantine worker machines. We propose and analyze a distributed gradient-descent algorithm that performs a simple thresholding based on gradient norms to mitigate Byzantine failures. We show the (statistical) error-rate of our algorithm matches that of Yin et al.~\cite{dong}, which uses more complicated schemes (coordinate-wise median, trimmed mean). Furthermore, for communication efficiency, we consider a generic class of $\delta$-approximate compressors from Karimireddi et al.~\cite{errorfeed} that encompasses sign-based compressors and top-$k$ sparsification. Our algorithm uses compressed gradients and gradient norms for aggregation and Byzantine removal respectively. We establish the statistical error rate for non-convex smooth loss functions. We show that, in certain range of the compression factor $\delta$, the (order-wise) rate of convergence is not affected by the compression operation. Moreover, we analyze the compressed gradient descent algorithm with error feedback (proposed in \cite{errorfeed}) in a distributed setting and in the presence of Byzantine worker machines. We show that exploiting error feedback improves the statistical error rate. Finally, we experimentally validate our results and show good performance in convergence for convex (least-square regression) and non-convex (neural network training) problems.


title: EvAn: Neuromorphic Event-based Anomaly Detection

abstract: Event-based cameras are bio-inspired novel sensors that asynchronously record changes in illumination in the form of events, thus resulting in significant advantages over conventional cameras in terms of low power utilization, high dynamic range, and no motion blur. Moreover, such cameras, by design, encode only the relative motion between the scene and the sensor (and not the static background) to yield a very sparse data structure, which can be utilized for various motion analytics tasks. In this paper, for the first time in event data analytics community, we leverage these advantages of an event camera towards a critical vision application - video anomaly detection. We propose to model the motion dynamics in the event domain with dual discriminator conditional Generative adversarial Network (cGAN) built on state-of-the-art architectures. To adapt event data for using as input to cGAN, we also put forward a deep learning solution to learn a novel representation of event data, which retains the sparsity of the data as well as encode the temporal information readily available from these sensors. Since there is no existing dataset for anomaly detection in event domain, we also provide an anomaly detection event dataset with an exhaustive set of anomalies. Careful analysis reveals that the proposed method results in huge reduction in computational complexity as compared to previous state-of-the-art conventional anomaly detection networks.


title: Fast Sparse ConvNets

abstract: Historically, the pursuit of efficient inference has been one of the driving forces behind research into new deep learning architectures and building blocks. Some recent examples include: the squeeze-and-excitation module, depthwise separable convolutions in Xception, and the inverted bottleneck in MobileNet v2. Notably, in all of these cases, the resulting building blocks enabled not only higher efficiency, but also higher accuracy, and found wide adoption in the field. In this work, we further expand the arsenal of efficient building blocks for neural network architectures; but instead of combining standard primitives (such as convolution), we advocate for the replacement of these dense primitives with their sparse counterparts. While the idea of using sparsity to decrease the parameter count is not new, the conventional wisdom is that this reduction in theoretical FLOPs does not translate into real-world efficiency gains. We aim to correct this misconception by introducing a family of efficient sparse kernels for ARM and WebAssembly, which we open-source for the benefit of the community as part of the XNNPACK library. Equipped with our efficient implementation of sparse primitives, we show that sparse versions of MobileNet v1, MobileNet v2 and EfficientNet architectures substantially outperform strong dense baselines on the efficiency-accuracy curve. On Snapdragon 835 our sparse networks outperform their dense equivalents by $1.3-2.4\times$ -- equivalent to approximately one entire generation of MobileNet-family improvement. We hope that our findings will facilitate wider adoption of sparsity as a tool for creating efficient and accurate deep learning architectures.


title: Information-Theoretic Confidence Bounds for Reinforcement Learning

abstract: We integrate information-theoretic concepts into the design and analysis of optimistic algorithms and Thompson sampling. By making a connection between information-theoretic quantities and confidence bounds, we obtain results that relate the per-period performance of the agent with its information gain about the environment, thus explicitly characterizing the exploration-exploitation tradeoff. The resulting cumulative regret bound depends on the agent's uncertainty over the environment and quantifies the value of prior information. We show applicability of this approach to several environments, including linear bandits, tabular MDPs, and factored MDPs. These examples demonstrate the potential of a general information-theoretic approach for the design and analysis of reinforcement learning algorithms.


title: Influence of sub- and super-solar metallicities on the compositions of solid planetary building blocks

abstract: The composition of the protoplanetary disc is linked to the composition of the host star, where a higher overall metallicity of the host star provides more building blocks for planets. However, most planet formation simulations only link the stellar iron abundance [Fe/H] to planet formation and [Fe/H] in itself is used as a proxy to scale all elements. But large surveys of stellar abundances show that this is not true. We use here stellar abundances from the GALAH surveys to determine the average detailed abundances of Fe, Si, Mg, O, and C for a broad range of [Fe/H] spanning from -0.4 to +0.4. Using an equilibrium chemical model that features the most important rock forming molecules as well as volatile contributions of H$_2$O, CO$_2$, CH$_4$ and CO, we calculate the chemical composition of solid planetary building blocks. Solid building blocks that are formed entirely interior to the water ice line (T>150K) only show an increase in Mg$_2$SiO$_4$ and a decrease in MgSiO$_3$ for increasing host star metallicity, related to the increase of Mg/Si for higher [Fe/H]. Solid planetary building blocks forming exterior to the water ice line (T<150K) show dramatic changes in their composition. The water ice content decreases from around $\sim$50\% at [Fe/H]=-0.4 to $\sim$6\% at [Fe/H]=0.4 in our chemical model. This is mainly caused by the increasing C/O ratio with increasing [Fe/H], which binds most of the oxygen in gaseous CO and CO$_2$, resulting in a small water ice fraction. Planet formation simulations coupled with the chemical model confirm these results by showing that the water ice content of super-Earths decreases with increasing host star metallicity due to the increased C/O ratio. This decrease of the water ice fraction has important consequences for planet formation, planetary composition and the eventual habitability of planetary systems formed around these high metallicity stars.


title: Large Effects from Small QCD Instantons: Making Soft Bombs at Hadron Colliders

abstract: It is a common belief that the last missing piece of the Standard Model of particles physics was found with the discovery of the Higgs boson at the Large Hadron Collider. However, there remains a major prediction of quantum tunnelling processes mediated by instanton solutions in the Yang-Mills theory, that is still untested in the Standard Model. The direct experimental observation of instanton-induced processes, which are a consequence of the non-trivial vacuum structure of the Standard Model and of quantum tunnelling in QFT, would be a major breakthrough in modern particle physics. In this paper, we present for the first time a full calculation of QCD instanton-induced processes in proton-proton collisions accounting for quantum corrections due to both initial and final state gluon interactions, a first implementation in an MC event generator as well as a basic strategy how to observe these effects experimentally.


title: On the closure requirement for VOF simulations with RANS modeling

abstract: The volume of fluid (VOF) method is increasingly used in computational fluid dynamics (CFD) simulations of turbulent two-phase flows. The Reynolds-Averaged Navier-Stokes (RANS) approach is an economic and practical way for turbulent VOF simulations. Even though RANS-VOF simulations are widely conducted, the underlying physics is barely discussed. This study reveals the very basic closure requirement for RANS-VOF simulations.