title: Vectors of type II Hermite-Pad\'e approximations and a new linear independence criterion

abstract: We propose a linear independence criterion, and outline an application of it. Down to its simplest case, it aims at solving this problem: given three real numbers, typically as special values of analytic functions, how to prove that the $\mathbb{Q}$-vector space spanned by $1$ and those three numbers has dimension at least 3, whenever we are unable to achieve full linear independence, by using simultaneous approximations, i.e. those usually arising from Hermite-Pad\'e approximations of type II and their suitable generalizations. It should be recalled that approximations of type I and II are related, at least in principle: when the numerical application consists in specializing actual functional constructions of the two types, they can be obtained, one from the other, as explained in a well-known paper by K.Mahler [34]. That relation is reflected in a relation between the asymptotic behavior of the approximations at the infinite place of $\mathbb{Q}$. Rather interestingly, the two view-points split away regarding the asymptotic behaviors at finite places (i.e. primes) of $\mathbb{Q}$, and this makes the use of type II more convenient for particular purposes. In addition, sometimes we know type II approximations to a given set of functions, for which type I approximations are not known explicitly. Our approach can be regarded as a dual version of the standard linear independence criterion, which goes back to Siegel.


title: Hasse--Witt matrices and mirror toric pencils

abstract: Mirror symmetry suggests unexpected relationships between arithmetic properties of distinct families of algebraic varieties. For example, Wan and others have shown that for some mirror pairs, the number of rational points over a finite field matches modulo the order of the field. In this paper, we obtain a similar result for certain mirror pairs of toric hypersurfaces. We use recent results by Huang, Lian, Yau and Yu describing the relationship between the Picard-Fuchs equations of these varieties and their Hasse--Witt matrices, which encapsulate information about the number of points. The result allows us to compute the number of points modulo the order of the field explicitly. We illustrate this by computing K3 surface examples related to hypergeometric functions.


title: Universal Geometric Graphs

abstract: We introduce and study the problem of constructing geometric graphs that have few vertices and edges and that are universal for planar graphs or for some sub-class of planar graphs; a geometric graph is \emph{universal} for a class $\mathcal H$ of planar graphs if it contains an embedding, i.e., a crossing-free drawing, of every graph in $\mathcal H$. Our main result is that there exists a geometric graph with $n$ vertices and $O(n \log n)$ edges that is universal for $n$-vertex forests; this extends to the geometric setting a well-known graph-theoretic result by Chung and Graham, which states that there exists an $n$-vertex graph with $O(n \log n)$ edges that contains every $n$-vertex forest as a subgraph. Our $O(n \log n)$ bound on the number of edges cannot be improved, even if more than $n$ vertices are allowed. We also prove that, for every positive integer $h$, every $n$-vertex convex geometric graph that is universal for $n$-vertex outerplanar graphs has a near-quadratic number of edges, namely $\Omega_h(n^{2-1/h})$; this almost matches the trivial $O(n^2)$ upper bound given by the $n$-vertex complete convex geometric graph. Finally, we prove that there exists an $n$-vertex convex geometric graph with $n$ vertices and $O(n \log n)$ edges that is universal for $n$-vertex caterpillars.


title: Two-timescale evolution of extreme-mass-ratio inspirals: waveform generation scheme for quasicircular orbits in Schwarzschild spacetime

abstract: Extreme-mass-ratio inspirals, in which a stellar-mass compact object spirals into a supermassive black hole in a galactic core, are expected to be key sources for LISA. Modelling these systems with sufficient accuracy for LISA science requires going to second (or {\em post-adiabatic}) order in gravitational self-force theory. Here we present a practical two-timescale framework for achieving this and generating post-adiabatic waveforms. The framework comprises a set of frequency-domain field equations that apply on the fast, orbital timescale, together with a set of ordinary differential equations that determine the evolution on the slow, inspiral timescale. Our analysis is restricted to the special case of quasicircular orbits around a Schwarzschild black hole, but its general structure carries over to the realistic case of generic (inclined and eccentric) orbits in Kerr spacetime. In our restricted context, we also develop a tool that will be useful in all cases: a formulation of the frequency-domain field equations using hyperboloidal slicing, which significantly improves the behavior of the sources near the boundaries. We give special attention to the slow evolution of the central black hole, examining its impact on both the two-timescale evolution and the earlier self-consistent evolution scheme.


title: Boosted Dark Matter Interpretation of the XENON1T Excess

abstract: We propose boosted dark matter (BDM) as a possible explanation for the excess of keV electron recoil events observed by XENON1T. BDM particles have velocities much larger than those typical of virialized dark matter, and, as such, BDM-electron scattering can naturally produce keV electron recoils. We show that the required BDM-electron scattering cross sections can be easily realized in a simple model with a heavy vector mediator. Though these cross sections are too large for BDM to escape from the Sun, the BDM flux can originate from the Galactic Center or from halo dark matter annihilations. Furthermore, a daily modulation of the BDM signal will be present, which could not only be used to differentiate it from various backgrounds, but would also provide important directional information for the BDM flux.


title: Proper scoring rules for evaluating asymmetry in density forecasting

abstract: This paper proposes a novel asymmetric continuous probabilistic score (ACPS) for evaluating and comparing density forecasts. It extends the proposed score and defines a weighted version, which emphasizes regions of interest, such as the tails or the center of a variable's range. A test is also introduced to statistically compare the predictive ability of different forecasts. The ACPS is of general use in any situation where the decision maker has asymmetric preferences in the evaluation of the forecasts. In an artificial experiment, the implications of varying the level of asymmetry in the ACPS are illustrated. Then, the proposed score and test are applied to assess and compare density forecasts of macroeconomic relevant datasets (US employment growth) and of commodity prices (oil and electricity prices) with particular focus on the recent COVID-19 crisis period.


title: An operator view of policy gradient methods

abstract: We cast policy gradient methods as the repeated application of two operators: a policy improvement operator $\mathcal{I}$, which maps any policy $\pi$ to a better one $\mathcal{I}\pi$, and a projection operator $\mathcal{P}$, which finds the best approximation of $\mathcal{I}\pi$ in the set of realizable policies. We use this framework to introduce operator-based versions of traditional policy gradient methods such as REINFORCE and PPO, which leads to a better understanding of their original counterparts. We also use the understanding we develop of the role of $\mathcal{I}$ and $\mathcal{P}$ to propose a new global lower bound of the expected return. This new perspective allows us to further bridge the gap between policy-based and value-based methods, showing how REINFORCE and the Bellman optimality operator, for example, can be seen as two sides of the same coin.


title: Fast Matrix Square Roots with Applications to Gaussian Processes and Bayesian Optimization

abstract: Matrix square roots and their inverses arise frequently in machine learning, e.g., when sampling from high-dimensional Gaussians $\mathcal{N}(\mathbf 0, \mathbf K)$ or whitening a vector $\mathbf b$ against covariance matrix $\mathbf K$. While existing methods typically require $O(N^3)$ computation, we introduce a highly-efficient quadratic-time algorithm for computing $\mathbf K^{1/2} \mathbf b$, $\mathbf K^{-1/2} \mathbf b$, and their derivatives through matrix-vector multiplication (MVMs). Our method combines Krylov subspace methods with a rational approximation and typically achieves $4$ decimal places of accuracy with fewer than $100$ MVMs. Moreover, the backward pass requires little additional computation. We demonstrate our method's applicability on matrices as large as $50,\!000 \times 50,\!000$ - well beyond traditional methods - with little approximation error. Applying this increased scalability to variational Gaussian processes, Bayesian optimization, and Gibbs sampling results in more powerful models with higher accuracy.


title: High-$T$$_\textrm{C}$ Superconductivity in Hydrogen Clathrates Mediated by Coulomb Interactions between Hydrogen and Central-Atom Electrons

abstract: The uniquely characteristic macrostructures of binary hydrogen-clathrate compounds $M$H$_\textrm{n}$ formed at high pressure, a cage of hydrogens surrounding a central-atom host, is theoretically predicted in various studies to include structurally stable phonon-mediated superconductors. High superconductive transition temperatures $T$$_\textrm{C}$ have thus far been measured for syntheses with $M$ = La, Y, and Th. In compressed LaH$_\textrm{10}$, independent studies report $T$$_\textrm{C}$ of 250 K and over 260 K, a maximum in $T$$_\textrm{C}$ with pressure $P$, and normal-state resistance scaling with temperature (suggesting unconventional pairing). According to reported band structure calculations of $Fm$$\bar{3}$$m$-phase LaH$_\textrm{10}$, the La is anionic, with the chemical valence electrons appearing evenly split between La and H$_\textrm{10}$. Thus, compressed LaH$_\textrm{10}$ contains the combination of structure, charge separation, and optimal balanced allocation of valence electrons for supporting unconventional high-$T$$_\textrm{C}$ superconductivity mediated by Coulomb interactions between electronic charges associated with La and H$_\textrm{10}$. A general expression for the optimal superconducting transition temperature for $M$H$_\textrm{n}$ clathrates is derived as $T$$_\textrm{C0}$ = $k$$_\textrm{B}$$^{-1}$$\Lambda$[(n + $v$)/2$A$]$^{1/2}$$e$$^{2}$/$\zeta$, where $\Lambda$ is a universal constant, (n + $v$) is the chemical valence sum per formula unit, taking unity for H and $v$ for atom $M$, $A$ is the surface area of the H-polyhedron cage, and $\zeta$ is the mean distance between the $M$ site and the centroids of the polyhedron faces. Applied to $Fm$$\bar{3}$$m$ LaH$_\textrm{10}$, $T$$_\textrm{C0}$ values of 249.8(1.3) K and 260.7(2.0) K are found for the two experiments. Associated attributes of charge allocation, structure, effective Coulomb potential, . . .


title: The Lang-Trotter Conjecture for products of non-CM elliptic curves

abstract: Inspired by the work of Lang-Trotter on the densities of primes with fixed Frobenius traces for elliptic curves defined over $\mathbb{Q}$ and by the subsequent generalization of Cojocaru-Davis-Silverberg-Stange to generic abelian varieties, we study the analogous question for abelian surfaces isogenous to products of non-CM elliptic curves over $\mathbb{Q}$. We formulate the corresponding conjectural asymptotic, provide upper bounds, and explicitly compute (when the elliptic curves lie outside a thin set) the arithmetically significant constants appearing in the asymptotic. This allows us to provide computational evidence for the conjecture.


title: Computational model on COVID-19 Pandemic using Probabilistic Cellular Automata

abstract: Coronavirus disease (COVID-19) which is caused by SARS-COV2 has become a pandemic. This disease is highly infectious and potentially fatal, causing a global public health concern. To contain the spread of COVID-19, governments are adopting nationwide interventions, like lockdown, containment and quarantine, restrictions on travel, cancelling social events and extensive testing. To understand the effects of these measures on the control of the epidemic in a data-driven manner, we propose a probabilistic cellular automata (PCA) based modified SEIQR model. The transitions associated with the model is driven by data available on chronology, symptoms, pathogenesis and transmissivity of the virus. By arguing that the lattice-based model captures the features of the dynamics along with the existing fluctuations, we perform rigorous computational analyses of the model to take into account of the spatial dynamics of social distancing measures imposed on the people. Considering the probabilistic behavioural aspects associated with mitigation strategies, we study the model considering factors like population density and testing efficiency. Using the model, we focus on the variability of epidemic dynamics data for different countries and point out the reasons behind these contrasting observations. To the best of our knowledge, this is the first attempt to model COVID-19 spread using PCA that gives us both spatial and temporal variations of the infection spread with the insight about the contributions of different infection parameters.

categories: physics.soc-ph cs.SI nlin.CG q-bio.PE q-bio.QM


title: Quantification of competing magnetic states and switching pathways in curved nanowires by direct dynamic imaging

abstract: For viable applications, spintronic devices based e.g. on domain wall motion need to be highly reliable with stable magnetization states and highly reproducible switching pathways transforming one state to another. The existence of multiple stable states and switching pathways in a system is a definitive barrier for device operation, yet rare and stochastic events are difficult to detect and understand. We demonstrate an approach to quantify competing magnetic states and stochastic switching pathways based on time-resolved scanning electron microscopy with polarization analysis, applied to the technologically relevant control of vortex domain wall chirality via field and curvature in curved wires. While being a pump-probe technique, our analysis scheme nonetheless allows for the disentanglement of different occurring dynamic pathways and we can even identify the rare events leading to changes from one magnetization switching pathway to another pathway via temperature- and geometry-dependent measurements. The experimental imaging is supported by micromagnetic simulations to reveal the mechanisms responsible for the change of the pathway. Together the results allow us to explain the origin and details of the domain wall chirality control and to quantify the frequency and the associated energy barriers of thermally activated changes of the states and switching pathways.


title: Uncertainty Principles for Quaternion Windowed Offset Linear Canonical Transform of Two Dimensional Signals

abstract: The offset linear canonical transform encompassing the numerous integral transforms, is a promising tool for analyzing non-stationary signals with more degrees of freedom. In this paper, we generalize the windowed offset linear canonical transform for quaternion-valued signals, by introducing a novel time-frequency transform namely the quaternion windowed offset linear canonical transform of 2D quaternion-valued signals. We initiate our investigation by studying some fundamental properties of the proposed transform including inner product relation, energy conservation, and reproducing formula by employing the machinery of quaternion offset linear canonical transforms. Some uncertainty principles such as Heisenberg-Weyl, logarithmic and local uncertainty principle are also derived for quaternion windowed offset linear canonical transform. Finally, we gave an example of quaternion windowed offset linear canonical transform.


title: Priming reasoning increases intentions to wear a face covering to slow down COVID-19 transmission

abstract: Finding mechanisms to promote the use of face masks is fundamental during the second phase of the COVID-19 pandemic response, when shelter-in-place rules are relaxed and some segments of the population are allowed to circulate more freely. Here we report three pre-registered studies (total N = 1,920), using an heterogenous sample of people living in the USA, showing that priming people to "rely on their reasoning" rather than to "rely on their emotions" significantly increases their intentions to wear a face covering. Compared to the baseline, priming reasoning promotes intentions to wear a face covering, whereas priming emotion has no significant effect. These findings have theoretical and practical implications. Practically, they offer a simple and scalable intervention to promote intentions to wear a face mask. Theoretically, they shed light on the cognitive basis of intentions to wear a face covering.


title: On Reward-Free Reinforcement Learning with Linear Function Approximation

abstract: Reward-free reinforcement learning (RL) is a framework which is suitable for both the batch RL setting and the setting where there are many reward functions of interest. During the exploration phase, an agent collects samples without using a pre-specified reward function. After the exploration phase, a reward function is given, and the agent uses samples collected during the exploration phase to compute a near-optimal policy. Jin et al. [2020] showed that in the tabular setting, the agent only needs to collect polynomial number of samples (in terms of the number states, the number of actions, and the planning horizon) for reward-free RL. However, in practice, the number of states and actions can be large, and thus function approximation schemes are required for generalization. In this work, we give both positive and negative results for reward-free RL with linear function approximation. We give an algorithm for reward-free RL in the linear Markov decision process setting where both the transition and the reward admit linear representations. The sample complexity of our algorithm is polynomial in the feature dimension and the planning horizon, and is completely independent of the number of states and actions. We further give an exponential lower bound for reward-free RL in the setting where only the optimal $Q$-function admits a linear representation. Our results imply several interesting exponential separations on the sample complexity of reward-free RL.


title: Center-based 3D Object Detection and Tracking

abstract: Three-dimensional objects are commonly represented as 3D boxes in a point-cloud. This representation mimics the well-studied image-based 2D bounding-box detection but comes with additional challenges. Objects in a 3D world do not follow any particular orientation, and box-based detectors have difficulties enumerating all orientations or fitting an axis-aligned bounding box to rotated objects. In this paper, we instead propose to represent, detect, and track 3D objects as points. Our framework, CenterPoint, first detects centers of objects using a keypoint detector and regresses to other attributes, including 3D size, 3D orientation, and velocity. In a second stage, it refines these estimates using additional point features on the object. In CenterPoint, 3D object tracking simplifies to greedy closest-point matching. The resulting detection and tracking algorithm is simple, efficient, and effective. CenterPoint achieved state-of-the-art performance on the nuScenes benchmark for both 3D detection and tracking, with 65.5 NDS and 63.8 AMOTA for a single model. On the Waymo Open Dataset, CenterPoint outperforms all previous single model method by a large margin and ranks first among all Lidar-only submissions. The code and pretrained models are available at https://github.com/tianweiy/CenterPoint.


title: An Improved Spin-Down Rate for the Proposed White-Dwarf Pulsar AR Scorpii

abstract: We analyze rapid-cadence, multiwavelength photometry of AR Scorpii from three observatories, covering five observing seasons. We measure the arrival times of the system's beat pulses and use them to compute an updated ephemeris. The white dwarf spin-down rate is estimated with an uncertainty of only 4%. These results confirm, beyond any doubt, that the white dwarf's spin period is increasing at the rate consistent with by that of Stiller et al. (2018). We study the evolution of the beat pulse's color index across the orbit. The color of the primary pulse maxima varies significantly across the orbit, with the peaks being bluer after superior conjunction than in the first half of the orbit. Specifically, at orbital phase 0.5, the color index of the primary pulse shows a very sharp discontinuity towards bluer indices. This supports the Potter & Buckley (2018b) synchrotron emission model where the two emitting poles differ significantly in color. However, no corresponding jump in the color of the secondary pulses is seen. Furthermore, our analysis reveals that the arrival times of the pulses can differ by as much as 6s in simultaneous $u$ and $r$ photometry, depending on the binary orbital phase. If left uncorrected, this wavelength-dependent timing offset could lead to erroneous measurements of the spin-period derivative, particularly with heterogeneous datasets.


title: Comments on: "A new additive decomposition of velocity gradient", by B. Sun

abstract: Comments on "A new additive decomposition of velocity gradient" [Phys. Fluids 31, 061702 (2019), arXiv:1908.01638] is presented


title: The MCC-F1 curve: a performance evaluation technique for binary classification

abstract: Many fields use the ROC curve and the PR curve as standard evaluations of binary classification methods. Analysis of ROC and PR, however, often gives misleading and inflated performance evaluations, especially with an imbalanced ground truth. Here, we demonstrate the problems with ROC and PR analysis through simulations, and propose the MCC-F1 curve to address these drawbacks. The MCC-F1 curve combines two informative single-threshold metrics, MCC and the F1 score. The MCC-F1 curve more clearly differentiates good and bad classifiers, even with imbalanced ground truths. We also introduce the MCC-F1 metric, which provides a single value that integrates many aspects of classifier performance across the whole range of classification thresholds. Finally, we provide an R package that plots MCC-F1 curves and calculates related metrics.


title: Distributionally Robust Profit Opportunities

abstract: This paper expands the notion of robust profit opportunities in financial markets to incorporate distributional uncertainty using Wasserstein distance as the ambiguity measure. Financial markets with risky and risk-free assets are considered. The infinite dimensional primal problems are formulated, leading to their simpler finite dimensional dual problems. A principal motivating question is how does distributional uncertainty help or hurt the robustness of the profit opportunity. Towards answering this question, some theory is developed and computational experiments are conducted. Finally some open questions and suggestions for future research are discussed.


title: Self-PU: Self Boosted and Calibrated Positive-Unlabeled Training

abstract: Many real-world applications have to tackle the Positive-Unlabeled (PU) learning problem, i.e., learning binary classifiers from a large amount of unlabeled data and a few labeled positive examples. While current state-of-the-art methods employ importance reweighting to design various risk estimators, they ignored the learning capability of the model itself, which could have provided reliable supervision. This motivates us to propose a novel Self-PU learning framework, which seamlessly integrates PU learning and self-training. Self-PU highlights three "self"-oriented building blocks: a self-paced training algorithm that adaptively discovers and augments confident positive/negative examples as the training proceeds; a self-calibrated instance-aware loss; and a self-distillation scheme that introduces teacher-students learning as an effective regularization for PU learning. We demonstrate the state-of-the-art performance of Self-PU on common PU learning benchmarks (MNIST and CIFAR-10), which compare favorably against the latest competitors. Moreover, we study a real-world application of PU learning, i.e., classifying brain images of Alzheimer's Disease. Self-PU obtains significantly improved results on the renowned Alzheimer's Disease Neuroimaging Initiative (ADNI) database over existing methods. The code is publicly available at: https://github.com/TAMU-VITA/Self-PU.


title: Self-interacting multistate boson stars

abstract: In this paper, we consider rotating multistate boson stars with quartic self-interactions. In contrast to the nodeless quartic-boson stars in \cite{Herdeiro:2015tia}, the self-interacting multistate boson stars (SIMBSs) have two types of nodes, including the $^1S^2S$ and $^1S^2P$ states. We show the mass $M$ of SIMBSs as a function of the synchronized frequency $\omega$, and the nonsynchronized frequency $\omega_2$ for three different cases. Moreover, for the case of two coexisting states with self-interacting potential, we study the mass $M$ of SIMBSs versus the angular momentum $J$ for the synchronized frequency $\omega$ and the nonsynchronized frequency $\omega_2$. Furthermore, for three different cases, we analyze the coexisting phase with both the ground and first excited states for SIMBSs. We also calculate the maximum value of coupling parameter $\Lambda$, and find the coupling parameter $\Lambda$ exists the finite range.


title: Efficient implementations of echo state network cross-validation

abstract: Background/introduction: Cross-Validation (CV) is still uncommon in time series modeling. Echo State Networks (ESNs), as a prime example of Reservoir Computing (RC) models, are known for their fast and precise one-shot learning, that often benefit from good hyper-parameter tuning. This makes them ideal to change the status quo. Methods: We discuss CV of time series for predicting a concrete time interval of interest, suggest several schemes for cross-validating ESNs and introduce an efficient algorithm for implementing them. This algorithm is presented as two levels of optimizations of doing $k$-fold CV. Training an RC model typically consists of two stages: (i) running the reservoir with the data and (ii) computing the optimal readouts. The first level of our optimization addresses the most computationally expensive part (i) and makes it remain constant irrespective of $k$. It dramatically reduces reservoir computations in any type of RC system and is enough if $k$ is small. The second level of optimization also makes the (ii) part remain constant irrespective of large $k$, as long as the dimension of the output is low. We discuss when the proposed validation schemes for ESNs could be beneficial, three options for producing the final model and empirically investigate them on six different real-world datasets, as well as do empirical computation time experiments. We provide the code in an online repository. Results: Proposed CV schemes give better and more stable test performance in all the six different real-world datasets, three task types. Empirical run times confirm our complexity analysis. Conclusions: In most situations $k$-fold CV of ESNs and many other RC models can be done for virtually the same time and space complexity as a simple single-split validation. This enables CV to become a standard practice in RC.


title: Spatial Concentration of Caching in Wireless Heterogeneous Networks

abstract: We propose a decentralized caching policy for wireless heterogeneous networks that makes content placement decisions based on pairwise interactions between cache nodes. We call our proposed scheme {\gamma}-exclusion cache placement (GEC), where a parameter {\gamma} controls an exclusion radius that discourages nearby caches from storing redundant content. GEC takes into account item popularity and the nodes' caching priorities and leverages negative dependence to relax the classic 0-1 knapsack problem to yield spatially balanced sampling across caches. We show that GEC guarantees a better concentration (reduced variance) of the required cache storage size than the state of the art, and that the cache size constraints can be satisfied with high probability. Given a cache hit probability target, we compare the 95\% confidence intervals of the required cache sizes for three caching schemes: (i) independent placement, (ii) hard exclusion caching (HEC), and (iii) the proposed GEC approach. For uniform spatial traffic, we demonstrate that GEC provides approximately a 3x and 2x reduction in required cache size over (i) and (ii), respectively. For non-uniform spatial traffic based on realistic peak-hour variations in urban scenarios, the gains are even greater.


title: Improving Locality Sensitive Hashing by Efficiently Finding Projected Nearest Neighbors

abstract: Similarity search in high-dimensional spaces is an important task for many multimedia applications. Due to the notorious curse of dimensionality, approximate nearest neighbor techniques are preferred over exact searching techniques since they can return good enough results at a much better speed. Locality Sensitive Hashing (LSH) is a very popular random hashing technique for finding approximate nearest neighbors. Existing state-of-the-art Locality Sensitive Hashing techniques that focus on improving performance of the overall process, mainly focus on minimizing the total number of IOs while sacrificing the overall processing time. The main time-consuming process in LSH techniques is the process of finding neighboring points in projected spaces. We present a novel index structure called radius-optimized Locality Sensitive Hashing (roLSH). With the help of sampling techniques and Neural Networks, we present two techniques to find neighboring points in projected spaces efficiently, without sacrificing the accuracy of the results. Our extensive experimental analysis on real datasets shows the performance benefit of roLSH over existing state-of-the-art LSH techniques.


title: Experimental Analysis of Locality Sensitive Hashing Techniques for High-Dimensional Approximate Nearest Neighbor Searches

abstract: Finding nearest neighbors in high-dimensional spaces is a fundamental operation in many multimedia retrieval applications. Exact tree-based indexing approaches are known to suffer from the notorious curse of dimensionality for high-dimensional data. Approximate searching techniques sacrifice some accuracy while returning good enough results for faster performance. Locality Sensitive Hashing (LSH) is a very popular technique for finding approximate nearest neighbors in high-dimensional spaces. Apart from providing theoretical guarantees on the query results, one of the main benefits of LSH techniques is their good scalability to large datasets because they are external memory based. The most dominant costs for existing LSH techniques are the algorithm time and the index I/Os required to find candidate points. Existing works do not compare both of these dominant costs in their evaluation. In this experimental survey paper, we show the impact of both these costs on the overall performance of the LSH technique. We compare three state-of-the-art techniques on four real-world datasets, and show that, in contrast to recent works, C2LSH is still the state-of-the-art algorithm in terms of performance while achieving similar accuracy as its recent competitors.


title: It has to be cool: on supergiant progenitors of binary black hole mergers from common-envelope evolution

abstract: Common-envelope (CE) evolution in massive binary systems is thought to be one of the most promising channels for the formation of compact binary mergers. In the case of merging binary black holes (BBHs), the essential CE phase takes place at a stage when the first BH is already formed and the companion star expands as a supergiant. We study which BH binaries with supergiant companions will evolve through and potentially survive a CE phase. To this end, we compute envelope binding energies from detailed massive stellar models at different evolutionary stages and metallicities. We make multiple physically extreme choices of assumptions that favor easier CE ejection as well as account for recent advancements in mass transfer stability criteria. We find that even with the most optimistic assumptions, a successful CE ejection in BH (and also NS) binaries is only possible if the donor is a massive convective-envelope giant, a red supergiant (RSG). In other words, pre-CE progenitors of BBH mergers are BH binaries with RSG companions. We find that due to its influence on the radial expansion of massive giants, metallicity has an indirect but a very strong effect on the envelope structure and binding energies of RSGs. Our results suggest that merger rates from population synthesis models could be severely overestimated, especially at low metallicity. Additionally, the lack of observed RSGs with luminosities above log($L/L_{\odot}$) = 5.6-5.8, corresponding to stars with $M > 40 M_{\odot}$, puts into question the viability of the CE channel for the formation of the most massive BBH mergers. Either such RSGs elude detection due to very short lifetimes, or they do not exist and the CE channel can only produce BBH systems with total mass $< 50 M_{\odot}$. We discuss an alternative CE scenario, in which a partial envelope ejection is followed by a phase of possibly long and stable mass transfer.


title: Discovering Symbolic Models from Deep Learning with Inductive Biases

abstract: We develop a general approach to distill symbolic representations of a learned deep model by introducing strong inductive biases. We focus on Graph Neural Networks (GNNs). The technique works as follows: we first encourage sparse latent representations when we train a GNN in a supervised setting, then we apply symbolic regression to components of the learned model to extract explicit physical relations. We find the correct known equations, including force laws and Hamiltonians, can be extracted from the neural network. We then apply our method to a non-trivial cosmology example-a detailed dark matter simulation-and discover a new analytic formula which can predict the concentration of dark matter from the mass distribution of nearby cosmic structures. The symbolic expressions extracted from the GNN using our technique also generalized to out-of-distribution data better than the GNN itself. Our approach offers alternative directions for interpreting neural networks and discovering novel physical principles from the representations they learn.

categories: cs.LG astro-ph.CO astro-ph.IM physics.comp-ph stat.ML


title: Recurrent Neutrino Emission from Supermassive Black Hole Mergers

abstract: The recent detection of possible neutrino emission from the blazar TXS 0506+056 was the first high-energy neutrino associated with an astrophysical source, making this special type of active galaxies promising neutrino emitters. The fact that two distinct episodes of neutrino emission were detected with a separation of around 3 years suggests that emission could be periodic. Periodic emission is expected from supermassive binary black hole systems due to jet precession close to the binary's merger. Here we show that if TXS 0506+056 is a binary source then the next neutrino flare could occur before the end of 2021. We derive the binary properties that would lead to the detection of gravitational waves from this system by LISA. Our results for the first time quantify the time scale of these correlations for the example of TXS 0506+056, providing clear predictions for both the neutrino and gravitational-wave signatures of such sources.


title: Shadow cast by a rotating charged black hole in quintessential dark energy

abstract: The existence of quintessential dark energy around a black hole has considerable consequences on its spacetime geometry. Hence, in this article, we explore its effect on horizons and the silhouette generated by a Kerr-Newman black hole in quintessential dark energy. Moreover, to analyze the deflection angle of light, we utilize the Gauss-Bonnet theorem. The obtained result demonstrates that, due to the dragging effect, the black hole spin elongates its shadow in the direction of the rotational axis, while increases the deflection angle. On the other hand, the black hole charge diminishing its shadow, as well as the angle of lights deflection. Besides, both spin and charge significantly increase the distortion effect in the black hole's shadow. The quintessence parameter \gamma, increases the shadow radius, while decreases the distortion effect at higher values of charge and spin parameters.


title: One residue to rule them all: Electroweak symmetry breaking, inflation and field-space geometry

abstract: We point out that the successful generation of the electroweak scale via gravitational instanton configurations in certain scalar-tensor theories can be viewed as the aftermath of a simple requirement: the existence of a quadratic pole with a sufficiently small residue in the Einstein-frame kinetic term for the Higgs field. In some cases, the inflationary dynamics may also be controlled by this residue and therefore related to the Fermi-to-Planck mass ratio, up to possible uncertainties associated with the instanton regularization. We present here a unified framework for this hierarchy generation mechanism, showing that the aforementioned residue can be associated with the curvature of the Einstein-frame target manifold in models displaying spontaneous breaking of dilatations. Our findings are illustrated through examples previously considered in the literature.


title: Entanglement harvesting with coherently delocalized matter

abstract: We study entanglement harvesting for matter systems such as atoms, ions or molecules, whose center of mass degrees of freedom are quantum delocalized and which couple to a relativistic quantum field. We employ a generalized Unruh-deWitt detector model for the light-matter interaction, and we investigate how the coherent spreading of the quantum center of mass wave function of two delocalized detector systems impacts their ability to become entangled with one another, via their respective interaction with a quantum field. For very massive detectors with initially highly localized centers of mass, we recover the results of entanglement harvesting for pointlike Unruh-deWitt detectors with classical center of mass degrees of freedom. We find that entanglement harvesting is Gaussian suppressed in the initial center of mass delocalization of the detectors. We further find that spatial smearing profiles, which are commonly employed to model the finite size of atoms due to their atomic orbitals, are not suited to model center of mass delocalization. Finally, for coherently delocalized detectors, we compare entanglement harvesting in the vacuum to entanglement harvesting in media. We find that entanglement harvesting is significantly suppressed in media in which the wave propagation speed is much smaller than the vacuum speed of light.


title: The Dual Conformal Box Integral in Minkowski Space

abstract: The dual conformal box integral in Minkowski space is not fully determined by the conformal invariants $z$ and $\bar{z}$. Depending on the kinematic region its value is on a 'branch' of the Bloch-Wigner function which occurs in the Euclidean case. Dual special conformal transformations in Minkowski space can change the kinematic region in such a way that the value of the integral jumps to another branch of this function, encoding a subtle breaking of dual conformal invariance for the integral. We classify conformally equivalent configurations of four points in compactified Minkowski space. We show that starting with any configuration, one can reach up to four branches of the integral using dual special conformal transformations. We also show that most configurations with real $z$ and $\bar{z}$ can be conformally mapped to a configuration in the same kinematic region with two points at infinity, where the box integral can be calculated directly in Minkowski space using only the residue theorem.


title: You Can't Always Get What You Want: The Impact of Prior Assumptions on Interpreting GW190412

abstract: GW190412 is the first observation of a black hole binary with definitively unequal masses. GW190412's mass asymmetry, along with the measured positive effective inspiral spin, allowed for inference of a component black hole spin: the primary black hole in the system was found to have a dimensionless spin magnitude between 0.17 and 0.59 (90% credible range). We investigate how the choice of priors for the spin magnitudes and tilts of the component black holes affect the robustness of parameter estimates for GW190412, and report Bayes factors across a suite of prior assumptions. Depending on the waveform family used to describe the signal, we find either marginal to moderate (2:1-6:1) or strong ($\gtrsim$ 20:1) support for the primary black hole being spinning compared to cases where only the secondary is allowed to have spin. We show how these choices influence parameter estimates, and find the asymmetric masses and positive effective inspiral spin of GW190412 to be qualitatively, but not quantitatively, robust to prior assumptions. Our results highlight the importance of both considering astrophysically motivated or population-based priors in interpreting observations and considering their relative support from the data.


title: Curvature homogeneous manifolds in dimension 4

abstract: We classify complete curvature homogeneous metrics on simply connected four dimensional manifolds which are invariant under a cohomogeneity one action. We show that they are either isometric to a symmetric space with one of its cohomogeneity one actions, or to a complete example by Tsukada on the normal bundle of the Veronese surface in CP^2. Along the way we show (in any dimension) that via an equivariant diffeomorphism the functions describing the metric can be partially diagonalized, a fact that may be useful for other problems as well


title: First principles modeling of exciton-polaritons in polydiacetylene chains

abstract: Exciton-polaritons in organic materials are hybrid states that result from the strong interaction of photons and the bound excitons that these materials host. Organic polaritons hold great interest for optoelectronic applications, however progress towards this end has been impeded by the lack of a first principles approach that quantifies light-matter interactions in these systems, and which would allow the formulation of molecular design rules. Here we develop such a first principles approach, quantifying light-matter interactions. We exemplify our approach by studying variants of the conjugated polymer polydiacetylene, and we show that a large polymer conjugation length is critical towards strong exciton-photon coupling, hence underlying the importance of pure structures without static disorder. By comparing to our experimental reflectivity measurements, we show that the coupling of excitons to vibrations, manifested by phonon side bands in the absorption, has a strong impact on the magnitude of light-matter coupling over a range of frequencies. Our approach opens the way towards a deeper understanding of polaritons in organic materials, and we highlight that a quantitatively accurate calculation of the exciton-photon interaction would require accounting for all sources of disorder self-consistently.


title: Three-component modelling of C-rich AGB star winds V. Effects of frequency-dependent radiative transfer including drift

abstract: Stellar winds of cool carbon stars enrich the interstellar medium with significant amounts of carbon and dust. We present a study of the influence of two-fluid flow on winds where we add descriptions of frequency-dependent radiative transfer. Our radiation hydrodynamic models in addition include stellar pulsations, grain growth and ablation, gas-to-dust drift using one mean grain size, dust extinction based on both the small particle limit and Mie scattering, and an accurate numerical scheme. We calculate models at high spatial resolution using 1024 gridpoints and solar metallicities at 319 frequencies, and we discern effects of drift by comparing drift models to non-drift models. Our results show differences of up to 1000 per cent in comparison to extant results. Mass-loss rates and wind velocities of drift models are typically, but not always, lower than in non-drift models. Differences are larger when Mie scattering is used instead of the small particle limit. Amongst other properties, the mass-loss rates of the gas and dust, dust-to-gas density ratio, and wind velocity show an exponential dependence on the dust-to-gas speed ratio. Yields of dust in the least massive winds increase by a factor four when drift is used. We find drift velocities in the range 10-67 km/s, which is drastically higher than in our earlier works that use grey radiative transfer. It is necessary to include an estimate of drift velocities to reproduce high yields of dust and low wind velocities.


title: Compressing combined probes: redshift weights for joint lensing and clustering analyses

abstract: Combining different observational probes, such as galaxy clustering and weak lensing, is a promising technique for unveiling the physics of the Universe with upcoming dark energy experiments. Whilst this strategy significantly improves parameter constraints, decreasing the degeneracies of individual analyses and controlling the systematics, processing data from tens of millions of galaxies is not a trivial task. In this work we derive and test a new estimator for joint clustering and lensing data analysis, maximising the scientific return and decreasing the computational cost. Our estimator compresses the data by up-weighting the components most sensitive to the parameters of interest, with no loss of information, taking into account information from the cross-correlation between the two probes. We derive optimal redshift weights which may be applied to individual galaxies when testing a given statistic and cosmological model.


title: Upper boundaries of AGN regions in optical diagnostic diagrams

abstract: The distribution of galaxies in optical diagnostic diagrams can provide information about their physical parameters when compared with ionization models under proper assumptions. By using a sample of central emitting regions from the MaNGA survey, we find evidence of the existence of upper boundaries for narrow-line regions (NLRs) of active galactic nuclei (AGN) in optical BPT diagrams, especially in the diagrams involving [S II]$\lambda \lambda$6716, 6731/H$\alpha$. Photoionization models can well reproduce the boundaries as a consequence of the decrease of [S II]$\lambda \lambda$6716, 6731/H$\alpha$ and [O III]$\lambda$5007/H$\beta$ ratios at very high metallicity. Whilst the exact location of the upper boundary in the [S II] BPT diagram only weakly depends on the electron density of the ionized cloud and the secondary nitrogen prescription, its dependence on the shapes of the input SEDs is much stronger. This allows us to constrain the power-law index of the AGN SED between 1 Ryd and $\sim100$ Ryd to be less than or equal to $-1.40\pm 0.05$. The coverage of the photoionization models in the [N II] BPT diagram has a stronger dependence on the electron density and the secondary nitrogen prescription. With the density constrained by the [S II] doublet ratio and the input SED constrained by the [S II] BPT diagram, we find that the extent of the data in the [N II] BPT diagram favors those prescriptions with high N/O ratios. Although shock-ionized cloud can produce similar line ratios as those by photoionization, the resulting shapes of the upper boundaries, if exist, would likely be different from those of a photoionizing origin.


title: Dynamics of one-dimensional quantum many-body systems in time-periodic linear potentials

abstract: We study a system of one-dimensional interacting quantum particles subjected to a time-periodic potential linear in space. After discussing the cases of driven one- and two-particles systems, we derive the analogous results for the many-particles case in presence of a general interaction two-body potential and the corresponding Floquet Hamiltonian. When the undriven model is integrable, the Floquet Hamitlonian is shown to be integrable too. We determine the micro-motion operator and the expression for a generic time evolved state of the system. We discuss various aspects of the dynamics of the system both at stroboscopic and intermediate times, in particular the motion of the center of mass of a generic wavepacket and its spreading over time. We also discuss the case of accelerated motion of the center of mass, obtained when the integral of the coeffcient strenght of the linear potential on a time period is non-vanishing, and we show that the Floquet Hamiltonian gets in this case an additional static linear potential. We also discuss the application of the obtained results to the Lieb-Liniger model.


title: From Demonstrations to Task-Space Specifications: Using Causal Analysis to Extract Rule Parameterization from Demonstrations

abstract: Learning models of user behaviour is an important problem that is broadly applicable across many application domains requiring human-robot interaction. In this work, we show that it is possible to learn generative models for distinct user behavioural types, extracted from human demonstrations, by enforcing clustering of preferred task solutions within the latent space. We use these models to differentiate between user types and to find cases with overlapping solutions. Moreover, we can alter an initially guessed solution to satisfy the preferences that constitute a particular user type by backpropagating through the learned differentiable models. An advantage of structuring generative models in this way is that we can extract causal relationships between symbols that might form part of the user's specification of the task, as manifested in the demonstrations. We further parameterize these specifications through constraint optimization in order to find a safety envelope under which motion planning can be performed. We show that the proposed method is capable of correctly distinguishing between three user types, who differ in degrees of cautiousness in their motion, while performing the task of moving objects with a kinesthetically driven robot in a tabletop environment. Our method successfully identifies the correct type, within the specified time, in 99% [97.8 - 99.8] of the cases, which outperforms an IRL baseline. We also show that our proposed method correctly changes a default trajectory to one satisfying a particular user specification even with unseen objects. The resulting trajectory is shown to be directly implementable on a PR2 humanoid robot completing the same task.


title: Gravitational waves affect vacuum entanglement

abstract: The entanglement harvesting protocol is an operational way to probe vacuum entanglement. This protocol relies on two atoms, modelled by Unruh-DeWitt detectors, that are initially unentangled. These atoms then interact locally with the field and become entangled. If the atoms remain spacelike separated, any entanglement between them is a result of entanglement that is `harvested' from the field. Thus, quantifying this entanglement serves as a proxy for how entangled the field is across the regions in which the atoms interacted. Using this protocol, it is demonstrated that while the transition probability of an individual inertial atom is unaffected by the presence of a gravitational wave, the entanglement harvested by two atoms depends sensitively on the frequency of the gravitational wave, exhibiting novel resonance effects when the energy gap of the detectors is tuned to the frequency of the gravitational wave. This suggests that the entanglement signature left by a gravitational wave may be useful in characterizing its properties, and potentially useful in exploring the gravitational-wave memory effect and gravitational-wave induced decoherence.


title: Eigenstate Entanglement Entropy in Random Quadratic Hamiltonians

abstract: The eigenstate entanglement entropy has been recently shown to be a powerful tool to distinguish integrable from generic quantum-chaotic models. In integrable models, a unique feature of the average eigenstate entanglement entropy (over all Hamiltonian eigenstates) is that the volume-law coefficient depends on the subsystem fraction. Hence, it deviates from the maximal (subsystem fraction independent) value encountered in quantum-chaotic models. Using random matrix theory for quadratic Hamiltonians, we obtain a closed-form expression for the average eigenstate entanglement entropy as a function of the subsystem fraction. We test its correctness against numerical results for the quadratic Sachdev-Ye-Kitaev model. We also show that it describes the average entanglement entropy of eigenstates of the power-law random banded matrix model (in the delocalized regime), and that it is close but not the same as the result for quadratic models that exhibit localization in quasimomentum space.

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


title: Constraining the energy spectrum of neutral pions in ultra-high-energy proton-air interactions

abstract: Fluctuations in the muon content of extensive air showers are anticorrelated to the fluctuations of the energy taken by the neutral pions which emerge from the first interaction of the cosmic ray in the atmosphere. We show that the high-energy tail of the neutral pion spectrum produced in the first proton-air interaction can be constrained, within the uncertainties of present cosmic ray experiments, through the analysis of the shower-to-shower distribution of the muon content of the air showers, $P(N_{\mu})$.


title: One of Everything: The Breakthrough Listen Exotica Catalog

abstract: We present Breakthrough Listen's "Exotica" Catalog as the centerpiece of our efforts to expand the diversity of targets surveyed in the Search for Extraterrestrial Intelligence (SETI). As motivation, we introduce the concept of survey breadth, the diversity of objects observed during a program. Several reasons for pursuing a broad program are given, including increasing the chance of a positive result in SETI, commensal astrophysics, and characterizing systematics. The Exotica Catalog is a 963 entry collection of 816 distinct targets intended to include "one of everything" in astronomy. It contains four samples: the Prototype sample, with an archetype of every known major type of non-transient celestial object; the Superlative sample of objects with the most extreme properties; the Anomaly sample of enigmatic targets that are in some way unexplained; and the Control sample with sources not expected to produce positive results. As far as we are aware, this is the first object list in recent times with the purpose of spanning the breadth of astrophysics. We share it with the community in hopes that it can guide treasury surveys and as a general reference work. Accompanying the catalog is extensive discussion of classification of objects and a new classification system for anomalies. Extensive notes on the objects in the catalog are available online. We discuss how we intend to proceed with observations in the catalog, contrast it with our extant Exotica efforts, and suggest similar tactics may be applied to other programs.


title: Generalization of Agent Behavior through Explicit Representation of Context

abstract: In order to deploy autonomous agents in digital interactive environments, they must be able to act robustly in unseen situations. The standard machine learning approach is to include as much variation as possible into training these agents. The agents can then interpolate within their training, but they cannot extrapolate much beyond it. This paper proposes a principled approach where a context module is coevolved with a skill module in the game. The context module recognizes the temporal variation in the game and modulates the outputs of the skill module so that the action decisions can be made robustly even in previously unseen situations. The approach is evaluated in the Flappy Bird and LunarLander video games, as well as in the CARLA autonomous driving simulation. The Context+Skill approach leads to significantly more robust behavior in environments that require extrapolation beyond training. Such a principled generalization ability is essential in deploying autonomous agents in real-world tasks, and can serve as a foundation for continual adaptation as well.


title: Kilonova Luminosity Function Constraints based on Zwicky Transient Facility Searches for 13 Neutron Star Mergers

abstract: We present a systematic search for optical counterparts to 13 gravitational wave (GW) triggers involving at least one neutron star during LIGO/Virgo's third observing run. We searched binary neutron star (BNS) and neutron star black hole (NSBH) merger localizations with the Zwicky Transient Facility (ZTF) and undertook follow-up with the Global Relay of Observatories Watching Transients Happen (GROWTH) collaboration. The GW triggers had a median localization of 4480 deg^2, median distance of 267 Mpc and false alarm rates ranging from 1.5 to 1e-25 per yr. The ZTF coverage had a median enclosed probability of 39%, median depth of 20.8mag, and median response time of 1.5 hr. The O3 follow-up by the GROWTH team comprised 340 UVOIR photometric points, 64 OIR spectra, and 3 radio. We find no promising kilonova (radioactivity-powered counterpart) and we convert the upper limits to constrain the underlying kilonova luminosity function. Assuming that all kilonovae are at least as luminous as GW170817 at discovery (-16.1mag), we calculate our joint probability of detecting zero kilonovae is only 4.2%. If we assume that all kilonovae are brighter than -16.6mag (extrapolated peak magnitude of GW170817) and fade at 1 mag/day (similar to GW170817), the joint probability of zero detections is 7%. If we separate the NSBH and BNS populations, the joint probability of zero detections, assuming all kilonovae are brighter than -16.6mag, is 9.7% for NSBH and 7.9% for BNS mergers. Moreover, <57% (<89%) of putative kilonovae could be brighter than -16.6mag assuming flat (fading) evolution, at 90% confidence. If we further account for the online terrestrial probability for each GW trigger, we find that <68% of putative kilonovae could be brighter than -16.6mag. Comparing to model grids, we find that some kilonovae must have Mej < 0.03 Msun or Xlan>1e-4 or phi>30deg to be consistent with our limits. (Abridged)


title: SkyHaul: An Autonomous Gigabit Network Fabric in the Sky

abstract: We design and build SKYHAUL, the first large-scale, autonomous, self-organizing network of Unmanned Aerial Vehicles (UAVs) that are connected using a mmWave wireless mesh backhaul. While the use of a mmWave backhaul paves the way for a new class of bandwidth-intensive, latency-sensitive cooperative applications (e.g., LTE coverage during disasters, surveillance during rescue in challenging terrains), the network of UAVs allows these applications to be executed at operating ranges that are far beyond the line-of-sight distances that limit individual UAVs today. To realize the challenging vision of deploying and maintaining an airborne mmWave mesh backhaul to cater to dynamic applications, SKYHAUL's design incorporates various elements: (1) Role-specific UAV operations that simultaneously address application tracking and backhaul connectivity (2) Novel algorithms to jointly address the problem of deployment (position, yaw of UAVs) and traffic routing across the UAV network; and (3) A provably optimal solution for fast and safe reconfiguration of UAV backhaul during application dynamics. We implement SKYHAUL on four DJI Matrice 600 Pros to demonstrate its practicality and performance through autonomous flight operations, complemented by large scale simulations.


title: Hot Exoplanet Atmospheres Resolved with Transit Spectroscopy (HEARTS) IV. A spectral inventory of atoms and molecules in the high-resolution transmission spectrum of WASP-121 b

abstract: Aims: We survey the transmission spectrum of WASP-121 b for line-absorption by metals and molecules at high spectral resolution, and elaborate on existing interpretations of the optical transmission spectrum observed with HST/STIS and WFC3. Methods: We use the cross-correlation technique and direct differential spectroscopy to search for sodium and other neutral and ionised atoms, TiO, VO and SH in high-resolution transit spectra obtained with the HARPS spectrograph. We inject models assuming chemical and hydrostatic equilibrium with varying temperature and composition to enable model comparison, and employ two bootstrap methods to test the robustness of our detections. Results: We detect neutral Mg, Na, Ca, Cr, Fe, Ni and V, which we predict exists in equilibrium with a significant quantity of VO, supporting earlier observations by HST/WFC3. Non-detections of Ti and TiO support the hypothesis that Ti is depleted via a cold-trap mechanism as has been proposed in the literature. Atomic line depths are under-predicted by hydrostatic models by a factor of 1.5 to 8, confirming recent findings that the atmosphere is extended. We predict the existence of significant concentrations of gas-phase TiO$_2$, VO$_2$ and TiS, which could be important absorbers at optical and NIR wavelengths in hot Jupiter atmospheres, but for which accurate line-list data is currently not available. We find no evidence for absorption by SH, and find that inflated atomic lines can plausibly explain the slope of the transmission spectrum observed in the NUV with HST/STIS. The Na D lines are significantly broadened and show a difference in their respective depths of 15 scale heights, which is not expected from isothermal hydrostatic theory.


title: Modelling Agent Policies with Interpretable Imitation Learning

abstract: As we deploy autonomous agents in safety-critical domains, it becomes important to develop an understanding of their internal mechanisms and representations. We outline an approach to imitation learning for reverse-engineering black box agent policies in MDP environments, yielding simplified, interpretable models in the form of decision trees. As part of this process, we explicitly model and learn agents' latent state representations by selecting from a large space of candidate features constructed from the Markov state. We present initial promising results from an implementation in a multi-agent traffic environment.


abstract: Deschamps and Mottez (hereafter DM) argue that the Gauss-Matuyama terrestrial magnetic field reversal may have left a vanishing main dipole moment to the field for a time of order 10,000 years. They say this may have allowed an enhanced cosmic ray flux, boosting the effect we proposed in Melott and Thomas (2019). We point out that the bulk of the cosmic ray flux from a nearby supernova should be too energetic, up to a million times more energetic than the limits of deflection by the terrestrial magnetic field. In fact, only those highly energetic ones will directly reach the troposphere, relevant for cloud-to-ground lightning. From Cosmic Explosions to Terrestrial Fires?: A Discussion. F. Deschamps and F. Mottez. J. Geology 128, online ahead of print. (2020) From Cosmic Explosions to Terrestrial Fires?: A Reply A.L. Melott and B.C. Thomas. J. Geology 128, online ahead of print. (2020) From cosmic explosions to terrestrial fires? (A.L. Melott and B.C. Thomas) Journal of Geology, 127, 475-481 10.1086/703418 (2019) [arXiv:1903.01501]


title: Remarks on blow-up phenomena in p-Laplacian heat equation with inhomogeneous nonlinearity

abstract: We investigate the $p-$Laplace heat equation $u_t-\Delta_p u=\zeta(t)f(u)$ on a bounded smooth domain $\Omega\subset\mathbb{R}^N$. Using differential inequalities arguments, we prove blow-up results under suitable conditions on $\zeta, f$, and the initial data $u_0$. We also give an upper bound for the blow-up time in each case.


title: Envy-freeness up to one item: Shall we add or remove resources?

abstract: We consider a fair division model in which agents have general valuations for bundles of indivisible items. We propose two new axiomatic properties for allocations in this model: EF1+- and EFX+-. We compare these with the existing EF1 and EFX. Although EF1 and EF1+- allocations often exist, our results assert eloquently that EFX+- and PO allocations exist in each case where EFX and PO allocations do not exist. Additionally, we prove several new impossibility and incompatibility results.


title: Information theoretic limits of learning a sparse rule

abstract: We consider generalized linear models in regimes where the number of nonzero components of the signal and accessible data points are sublinear with respect to the size of the signal. We prove a variational formula for the asymptotic mutual information per sample when the system size grows to infinity. This result allows us to derive an expression for the minimum mean-square error (MMSE) of the Bayesian estimator when the signal entries have a discrete distribution with finite support. We find that, for such signals and suitable vanishing scalings of the sparsity and sampling rate, the MMSE is nonincreasing piecewise constant. In specific instances the MMSE even displays an all-or-nothing phase transition, that is, the MMSE sharply jumps from its maximum value to zero at a critical sampling rate. The all-or-nothing phenomenon has previously been shown to occur in high-dimensional linear regression. Our analysis goes beyond the linear case and applies to learning the weights of a perceptron with general activation function in a teacher-student scenario. In particular, we discuss an all-or-nothing phenomenon for the generalization error with a sublinear set of training examples.


title: On the Chow ring of Fano varieties of type $S2$

abstract: We show that certain Fano eightfolds (obtained as hyperplane sections of an orthogonal Grassmannian, and studied by Ito-Miura-Okawa-Ueda and by Fatighenti-Mongardi) have a multiplicative Chow-K\"unneth decomposition. As a corollary, the Chow ring of these eightfolds behaves like that of K3 surfaces.


title: Classifying groups with a small number of subgroups

abstract: We provide lower bounds on the number of subgroups of a group $G$ as a function of the primes and exponents appearing in the prime factorization of $|G|$. Using these bounds, we classify all abelian groups with 22 or fewer subgroups, and all non-abelian groups with 19 or fewer subgroups. This allows us to extend the integer sequence A274847 \cite{OEIS} introduced by Slattery in \cite{Slattery}.


title: SqueezeBERT: What can computer vision teach NLP about efficient neural networks?

abstract: Humans read and write hundreds of billions of messages every day. Further, due to the availability of large datasets, large computing systems, and better neural network models, natural language processing (NLP) technology has made significant strides in understanding, proofreading, and organizing these messages. Thus, there is a significant opportunity to deploy NLP in myriad applications to help web users, social networks, and businesses. In particular, we consider smartphones and other mobile devices as crucial platforms for deploying NLP models at scale. However, today's highly-accurate NLP neural network models such as BERT and RoBERTa are extremely computationally expensive, with BERT-base taking 1.7 seconds to classify a text snippet on a Pixel 3 smartphone. In this work, we observe that methods such as grouped convolutions have yielded significant speedups for computer vision networks, but many of these techniques have not been adopted by NLP neural network designers. We demonstrate how to replace several operations in self-attention layers with grouped convolutions, and we use this technique in a novel network architecture called SqueezeBERT, which runs 4.3x faster than BERT-base on the Pixel 3 while achieving competitive accuracy on the GLUE test set. The SqueezeBERT code will be released.


title: The path integral of 3D gravity near extremality; or, JT gravity with defects as a matrix integral

abstract: We propose that a class of new topologies, for which there is no classical solution, should be included in the path integral of three-dimensional pure gravity, and that their inclusion solves pathological negativities in the spectrum, replacing them with a nonperturbative shift of the BTZ extremality bound. We argue that a two-dimensional calculation using a dimensionally reduced theory captures the leading effects in the near extremal limit. To make this argument, we study a closely related two-dimensional theory of Jackiw-Teitelboim gravity with dynamical defects. We show that this theory is equivalent to a matrix integral.


title: Coupled skyrmion breathing modes in synthetic ferri- and antiferromagnets

abstract: We present micromagnetic simulations of the dynamic GHz-range resonance modes of skyrmions excited by either out-of-plane ac magnetic fields or spin torques in prototypical synthetic ferri- and antiferromagnetic trilayer structures. The observed features in the calculated power spectra exhibit a systematic dependence on the coupling strength between the individual magnetic layers and are related to pure in-phase and anti-phase breathing modes as well as to hybridizations of breathing and spin-wave modes that are characteristic for the considered circular-shaped geometry. The experimental detection of these resonant oscillation modes may provide a means for skyrmion sensing applications and for the general characterization of skyrmion states in multilayer stacks with antiferromagnetic interlayer exchange coupling.


title: The Challenge of Stochastic St{\o}rmer-Verlet Thermostats Generating Correct Statistics

abstract: In light of the recently developed complete GJ set of single random variable stochastic, discrete-time St{\o}rmer-Verlet algorithms for statistically accurate simulations of Langevin equations, we investigate two outstanding questions: 1) Are there any algorithmic or statistical benefits from including multiple random variables per time-step, and 2) are there objective reasons for using one or more methods from the available set of statistically correct algorithms? To address the first question, we assume a general form for the discrete-time equations with two random variables and then follow the systematic, brute-force GJ methodology by enforcing correct thermodynamics in linear systems. It is concluded that correct configurational Boltzmann sampling of a particle in a harmonic potential implies correct configurational free-particle diffusion, and that these requirements only can be accomplished if the two random variables per time step are identical. We consequently submit that the GJ set represents all possible stochastic St{\o}rmer-Verlet methods that can reproduce time-step-independent statistics of linear systems. The second question is thus addressed within the GJ set. Based in part on numerical simulations of complex molecular systems, and in part on analytic scaling of time, we analyze the apparent difference in stability between different methods. We attribute this difference to the inherent time scaling in each method, and suggest that this scaling may lead to inconsistencies in the interpretation of dynamical and statistical simulation results. We therefore suggest that the method with the least inherent time-scaling, the GJ-I/GJF-2GJ method, be preferred for statistical applications where spurious rescaling of time is undesirable.


title: Magnetism and its coexistence with superconductivity in CaK(Fe$_{0.949}$Ni$_{0.051}$)$_4$As$_4$: muon spin rotation/relaxation studies

abstract: The magnetic response of CaK(Fe$_{0.949}$Ni$_{0.051}$)$_4$As$_4$ was investigated by means of the muon-spin rotation/relaxation. The long-range commensurate magnetic order sets in below the N\'{e}el temperature $T_{\rm N}= 50.0(5)$~K. The density-functional theory calculations have identified three possible muon stopping sites. The experimental data were found to be consistent with only one type of magnetic structure, namely, the long-range magnetic spin-vortex-crystal order with the hedgehog motif within the $ab-$plane and the antiferromagnetic stacking along the $c-$direction. The value of the ordered magnetic moment at $T\approx3$ K was estimated to be $m_{\rm Fe}=0.38(11)$ $\mu_{\rm B}$ ($\mu_{\rm B}$ is the Bohr magneton). A microscopic coexistence of magnetic and superconducting phases accompanied by a reduction of the magnetic order parameter below the superconducting transition temperature $T_{\rm c}\simeq 9$ K is observed. Comparison with 11, 122, and 1144 families of Fe-based pnictides points to existence of correlation between the reduction of the magnetic order parameter at $T\rightarrow 0$ and the ratio of the transition temperatures $T_{\rm c}/T_{\rm N}$. Such correlations were found to be described by Machida's model for coexistence of itinerant spin-density wave magnetism and superconductivity [Machida, J. Phys. Soc. Jpn. 50, 2195 (1981) and Bud'ko et al., Phys. Rev. B 98, 144520 (2018)].


title: AutoOD: Automated Outlier Detection via Curiosity-guided Search and Self-imitation Learning

abstract: Outlier detection is an important data mining task with numerous practical applications such as intrusion detection, credit card fraud detection, and video surveillance. However, given a specific complicated task with big data, the process of building a powerful deep learning based system for outlier detection still highly relies on human expertise and laboring trials. Although Neural Architecture Search (NAS) has shown its promise in discovering effective deep architectures in various domains, such as image classification, object detection, and semantic segmentation, contemporary NAS methods are not suitable for outlier detection due to the lack of intrinsic search space, unstable search process, and low sample efficiency. To bridge the gap, in this paper, we propose AutoOD, an automated outlier detection framework, which aims to search for an optimal neural network model within a predefined search space. Specifically, we firstly design a curiosity-guided search strategy to overcome the curse of local optimality. A controller, which acts as a search agent, is encouraged to take actions to maximize the information gain about the controller's internal belief. We further introduce an experience replay mechanism based on self-imitation learning to improve the sample efficiency. Experimental results on various real-world benchmark datasets demonstrate that the deep model identified by AutoOD achieves the best performance, comparing with existing handcrafted models and traditional search methods.


title: Two-photon Doppler-free ultraviolet laser spectroscopy on sulphur atoms

abstract: The $3p^{4}$ $^{3}$P$_{J}$ - $3p^{3}4p$ $^{3}$P$_{J}$ transition in the sulphur atom is investigated in a precision two-photon excitation scheme under Doppler-free and collision-free circumstances yielding an absolute accuracy of 0.0009 cm$^{-1}$, using a narrowband pulsed laser. This verifies and improves the level separations between amply studied odd parity levels with even parity levels in S I. An improved value for the $^{3}$P$_{2}$ - $^{3}$P$_{1}$ ground state fine structure splitting is determined at $396.0564$ (7) cm$^{-1}$. A $^{34}$S - $^{32}$S atomic isotope shift was measured from combining time-of-flight mass spectrometry with laser spectroscopy.


title: Semiclassical theory of the circular photogalvanic effect in gyrotropic systems

abstract: We develop a theory of circular photogalvanic effect (CPGE) for classically high photon energies which exceed the electron scattering rate but are small compared to the average electron kinetic energy. In this frequency range one can calculate the CPGE by using two different approaches. In the fully quantum-mechanical approach we find the photocurrent density by applying Fermi's golden rule for indirect intraband optical transitions with virtual intermediate states both in the conduction and valence bands. In the framework of the semiclassical approach, we apply a generalized Boltzmann equation with accounts for the Berry-curvature induced anomalous velocity, side jumps and skew scattering. The calculation is carried out for a wurtzite symmetry crystal. Both methods yield the same results for the CPGE current demonstrating consistency between the two approaches and applicability of the semiclassical theory for the description of nonlinear high-frequency transport.


title: The effect of metric behavior at spatial infinity on pointwise wave decay in the asymptotically flat stationary setting

abstract: The current work considers solutions to the wave equation on asymptotically flat, stationary, Lorentzian spacetimes in (1+3) dimensions. We investigate the relationship between the rate at which the geometry tends to flat and the pointwise decay rate of solutions. The case where the spacetime tends toward flat at a rate of $|x|^{-1}$ was studied in \cite{tat2013}, where a $t^{-3}$ pointwise decay rate was established. Here we extend the result to geometries tending toward flat at a rate of $|x|^{-\kappa}$ and establish a pointwise decay rate of $t^{-\kappa-2}$ for $\kappa \in \mathbb{N}$ with $\kappa \ge 2$. We assume a weak local energy decay estimate holds, which restricts the geodesic trapping allowed on the underlying geometry. We use the resolvent to connect the time Fourier Transform of a solution to the Cauchy data. Ultimately the rate of pointwise wave decay depends on the low frequency behavior of the resolvent, which is sensitive to the rate at which the background geometry tends to flat.


title: Self-Supervised Prototypical Transfer Learning for Few-Shot Classification

abstract: Most approaches in few-shot learning rely on costly annotated data related to the goal task domain during (pre-)training. Recently, unsupervised meta-learning methods have exchanged the annotation requirement for a reduction in few-shot classification performance. Simultaneously, in settings with realistic domain shift, common transfer learning has been shown to outperform supervised meta-learning. Building on these insights and on advances in self-supervised learning, we propose a transfer learning approach which constructs a metric embedding that clusters unlabeled prototypical samples and their augmentations closely together. This pre-trained embedding is a starting point for few-shot classification by summarizing class clusters and fine-tuning. We demonstrate that our self-supervised prototypical transfer learning approach ProtoTransfer outperforms state-of-the-art unsupervised meta-learning methods on few-shot tasks from the mini-ImageNet dataset. In few-shot experiments with domain shift, our approach even has comparable performance to supervised methods, but requires orders of magnitude fewer labels.


title: An effective field theory approach to monopolium

abstract: In this work we investigate the interaction between spin-zero and spin-one monopoles by making use of an effective field theory based on two-body and four-body interaction parts. In particular, we analyze the formation of bound state of monopole-antimonopole (i.e. monopolium). The magnetic-charge conjugation symmetry is studied in analogy to the usual charge conjugation to define a particle basis, for which we find bound-state solutions with relatively small binding energies and which allows us to identify the bounds on the parameters in the effective Lagrangians. Estimations of their masses, binding energies and scattering lengths are performed as functions of monopole masses and interaction strength in a specific renormalization scheme. We also examine the general validity of the approach and the feasibility of detecting the monopolium.


title: Teaching Programming to Novices: A Large-scale Analysis of App Inventor Projects

abstract: Teaching programming to K-12 students has become essential. In this context, App Inventor is a popular block-based programming environment used by a wide audience, from K-12 to higher education, including end-users to create mobile applications to support their primary job or hobbies. Although learning programming with App Inventor has been investigated, a question that remains is which programming concepts are typically used and how this compares to other block-based programming environments. Therefore, we explore the characteristics of App Inventor projects through a large-scale analysis of 88,606 apps from the App Inventor Gallery. We discovered that the size of App Inventor projects varies from projects with very few blocks to some surprisingly large projects with more than 60,000 blocks. In general, much fewer design components are used than programming blocks, as typically, to work properly, several programming blocks are necessary for each design component in an App Inventor project. In addition, we also compare our results with the analysis of 233,491 Scratch projects reported by Aivaloglou and Hermans [4]. Several differences can be observed, as in App Inventor projects events are more predominant, with lesser use of conditionals and loops. These findings may guide the decision on the adoption of App Inventor for teaching computing depending on the specific learning objectives or indicate the need for tailoring the curricula.


title: Class Normalization for (Continual)? Generalized Zero-Shot Learning

abstract: Normalization techniques have proved to be a crucial ingredient of successful training in a traditional supervised learning regime. However, in the zero-shot learning (ZSL) world, these ideas have received only marginal attention. This work studies normalization in ZSL scenario from both theoretical and practical perspectives. First, we give a theoretical explanation to two popular tricks used in zero-shot learning: normalize+scale and attributes normalization and show that they help training by preserving variance during a forward pass. Next, we demonstrate that they are insufficient to normalize a deep ZSL model and propose Class Normalization (CN): a normalization scheme, which alleviates this issue both provably and in practice. Third, we show that ZSL models typically have more irregular loss surface compared to traditional classifiers and that the proposed method partially remedies this problem. Then, we test our approach on 4 standard ZSL datasets and outperform sophisticated modern SotA with a simple MLP optimized without any bells and whistles and having ~50 times faster training speed. Finally, we generalize ZSL to a broader problem -- continual ZSL, and introduce some principled metrics and rigorous baselines for this new setup. The project page is located at https://universome.github.io/class-norm.


title: Quantum and classical approaches in statistical physics: some basic inequalities

abstract: We present some basic inequalities between the classical and quantum values of free energy, entropy and mean energy. We investigate the transition from the deterministic case (classical mechanics) to the probabilistic case (quantum mechanics). In the first part of the paper, we assume that the reduced Planck constant $\hbar$, the absolute temperature $T$, the frequency of an oscillator $\omega$, and the degree of freedom of a system $N$ are fixed. This approach to the problem of comparing quantum and classical mechanics is new (see [35]--[37]). In the second part of the paper, we simultaneously derive the semiclassical limits for four cases, that is, for $\hbar{\to}0$, $T{\to}\infty$, $\omega{\to}0$, and $N{\to}\infty$. We note that only the case $\hbar{\to}0$ is usually considered in quantum mechanics (see [21]). The cases $T{\to}\infty$ and $\omega{\to}0$ in quantum mechanics were initially studied by M. Planck and by A. Einstein, respectively.


title: Pycro-manager: open-source software for integrated microscopy hardware control and image processing

abstract: {\mu}Manager, an open-source microscopy acquisition software, has been an essential tool for many microscopy experiments over the past 15 years, but is not easy to use for experiments in which image acquisition and analysis are closely coupled. This is because {\mu}Manager libraries are written in C++ and Java, whereas image processing is increasingly carried out with data science and machine learning tools most easily accessible through the Python programming language. We present Pycro-Manager, a tool that enables rapid development of such experiments, while also providing access to the wealth of existing tools within {\mu}Manager through Python.


title: Polymer escape through a three dimensional Double-Nanopore System

abstract: We study escape dynamics of a double-stranded DNA (dsDNA) through an idealized double nanopore (DNP) geometry subject to two equal and opposite forces (tug-of-war) using Brownian dynamics (BD) simulation. In addition to the geometrical restrictions imposed on the cocaptured dsDNA segment in between the pores, the presence of tug-of-war forces at each pore results in a variation of the local chain stiffness for the segment of the chain in between the pores which increases the overall stiffness of the chain. We use BD simulation results to understand how the intrinsic chain stiffness and the TOW forces affect the escape dynamics by monitoring the local chain persistence length $\ell_p$, the residence time of the individual monomers $W(m)$ in the nanopores, and the chain length dependence of the escape time $\langle \tau \rangle$ and its distribution. Finally, we generalize the scaling theory for the unbiased single nanopore translocation for a fully flexible chain for the escape of a semi-flexible chain through a DNP in presence of TOW forces. We establish that the stiffness dependent part of the escape time is approximately independent of the translocation mechanism so that $\langle \tau \rangle \sim \ell_p^{2/D+2}$, and therefore the generalized escape time for a semi-flexible chain can be written as $\langle \tau \rangle = AN^\alpha\ell_p^{2/D+2}$. We use BD simulation results to compare the predictions of the scaling theory. Our numerical studies supplemented by scaling analysis provide fundamental insights to design new experiments where a dsDNA moves slowly through a series of graphene nanopores.


title: Quantum Geometric Machine Learning for Quantum Circuits and Control

abstract: The application of machine learning techniques to solve problems in quantum control together with established geometric methods for solving optimisation problems leads naturally to an exploration of how machine learning approaches can be used to enhance geometric approaches to solving problems in quantum information processing. In this work, we review and extend the application of deep learning to quantum geometric control problems. Specifically, we demonstrate enhancements in time-optimal control in the context of quantum circuit synthesis problems by applying novel deep learning algorithms in order to approximate geodesics (and thus minimal circuits) along Lie group manifolds relevant to low-dimensional multi-qubit systems, such as SU(2), SU(4) and SU(8). We demonstrate the superior performance of greybox models, which combine traditional blackbox algorithms with prior domain knowledge of quantum mechanics, as means of learning underlying quantum circuit distributions of interest. Our results demonstrate how geometric control techniques can be used to both (a) verify the extent to which geometrically synthesised quantum circuits lie along geodesic, and thus time-optimal, routes and (b) synthesise those circuits. Our results are of interest to researchers in quantum control and quantum information theory seeking to combine machine learning and geometric techniques for time-optimal control problems.


title: Energy Reflection and Transmission at 2D Holographic Interfaces

abstract: Scattering from conformal interfaces in two dimensions is universal in that the flux of reflected and transmitted energy does not depend on the details of the initial state. In this letter, we present the first gravitational calculation of energy reflection and transmission coefficients for interfaces with thin-brane holographic duals. Our result for the reflection coefficient depends monotonically on the tension of the dual string anchored at the interface, and obeys the lower bound recently derived from the ANEC in conformal field theory. The B(oundary)CFT limit is recovered for infinite ratio of the central charges.


title: The computational strength of matchings in countable graphs

abstract: In a 1977 paper, Steffens identified an elegant criterion for determining when a countable graph has a perfect matching. In this paper, we will investigate the proof-theoretic strength of this result and related theorems. We show that a number of natural variants of these theorems are equivalent, or closely related, to the ``big five'' subsystems of reverse mathematics. The results of this paper explore the relationship between graph theory and logic by showing the way in which specific changes to a single graph-theoretic principle impact the corresponding proof-theoretical strength. Taken together, the results and questions of this paper suggest that the existence of matchings in countable graphs provides a rich context for understanding reverse mathematics more broadly.


title: Spin-orbit-proximitized ferromagnetic metal by monolayer transition metal dichalcogenide: Atlas of spectral functions, spin textures and spin-orbit torques in Co/MoSe$_2$, Co/WSe$_2$ and Co/TaSe$_2$ heterostructures

abstract: The bilayer heterostructures composed of an ultrathin ferromagnetic metal (FM) and a material hosting strong spin-orbit (SO) coupling are principal resource for SO torque and spin-to-charge conversion nonequilibrium effects in spintronics. We demonstrate how hybridization of wavefunctions of Co layer and a monolayer of transition metal dichalcogenides (TMDs)---such as semiconducting MoSe$_2$ and WSe$_2$ or metallic TaSe$_2$---can lead to dramatic transmutation of electronic and spin structure of Co within some distance away from its interface with TMD, when compared to the bulk of Co or its surface in contact with vacuum. This is due to proximity induced SO splitting of Co bands encoded in the spectral functions and spin textures on its monolayers, which we obtain using noncollinear density functional theory (ncDFT) combined with equilibrium Green function (GF) calculations. In fact, SO splitting is present due to structural inversion asymmetry of the bilayer even if SO coupling within TMD monolayer is artificially switched off in ncDFT calculations, but switching it on makes the effects associated with proximity SO coupling within Co layer about five times larger. Injecting spin-unpolarized charge current through SO-proximitized monolayers of Co generates nonequilibrium spin density over them, so that its cross product with the magnetization of Co determines SO torque. The SO torque computed via first-principles quantum transport methodology, which combines ncDFT with nonequilibrium GF calculations, can be used as the screening parameter to identify optimal combination of materials and their interfaces for applications in spintronics. In particular, we identify heterostructure two-monolayer-Co/monolayer-WSe$_2$ as the most optimal.


title: Behavior of Limited Memory BFGS when Applied to Nonsmooth Functions and their Nesterov Smoothings

abstract: The motivation to study the behavior of limited-memory BFGS (L-BFGS) on nonsmooth optimization problems is based on two empirical observations: the widespread success of L-BFGS in solving large-scale smooth optimization problems, and the effectiveness of the full BFGS method in solving small to medium-sized nonsmooth optimization problems, based on using a gradient, not a subgradient, oracle paradigm. We first summarize our theoretical results on the behavior of the scaled L-BFGS method with one update applied to a simple convex nonsmooth function that is unbounded below, stating conditions under which the method converges to a non-optimal point regardless of the starting point. We then turn to empirically investigating whether the same phenomenon holds more generally,focusing on a difficult problem of Nesterov, as well as eigenvalue optimization problems arising in semidefinite programming applications. We find that when applied to a nonsmooth function directly, L-BFGS, especially its scaled variant, often breaks down with a poor approximation to an optimal solution, in sharp contrast to full BFGS. Unscaled L-BFGS is less prone to breakdown but conducts far more function evaluations per iteration than scaled L-BFGS does, and thus it is slow. Nonetheless, it is often the case that both variants obtain better results than the provably convergent, but slow, subgradient method. On the other hand, when applied to Nesterov's smooth approximation of a nonsmooth function, scaled L-BFGS is generally much more efficient than unscaled L-BFGS, often obtaining good results even when the problem is quite ill-conditioned. Summarizing, for large-scale nonsmooth optimization problems for which full BFGS and other methods for nonsmooth optimization are not practical, it is often better to apply L-BFGS to a smooth approximation of a nonsmooth problem than to apply it directly to the nonsmooth problem.


title: Image Sentiment Transfer

abstract: In this work, we introduce an important but still unexplored research task -- image sentiment transfer. Compared with other related tasks that have been well-studied, such as image-to-image translation and image style transfer, transferring the sentiment of an image is more challenging. Given an input image, the rule to transfer the sentiment of each contained object can be completely different, making existing approaches that perform global image transfer by a single reference image inadequate to achieve satisfactory performance. In this paper, we propose an effective and flexible framework that performs image sentiment transfer at the object level. It first detects the objects and extracts their pixel-level masks, and then performs object-level sentiment transfer guided by multiple reference images for the corresponding objects. For the core object-level sentiment transfer, we propose a novel Sentiment-aware GAN (SentiGAN). Both global image-level and local object-level supervisions are imposed to train SentiGAN. More importantly, an effective content disentanglement loss cooperating with a content alignment step is applied to better disentangle the residual sentiment-related information of the input image. Extensive quantitative and qualitative experiments are performed on the object-oriented VSO dataset we create, demonstrating the effectiveness of the proposed framework.


title: Distinguishing the Rigidity Dependences of Acceleration and Transport in Solar Energetic Particles

abstract: In solar energetic particle (SEP) events, the power-law dependence of element abundance enhancements on their mass-to-charge ratios A/Q provides a new tool that measures the combined rigidity dependences from both acceleration and transport. Distinguishing these two processes can be more challenging. However, the effects of acceleration dominate when SEP events are small or when the ions even propagate scatter-free, and transport can dominate the time evolution of large events with streaming-limited intensities. Magnetic reconnection in solar jets produces positive powers of A/Q from +2 to +7 and shock acceleration produces mostly negative powers from -2 to +1 in small and moderate SEP events where transport effects are minimal. This variation in the rigidity dependence of shock acceleration may reflect the non-planer structure, complexity, and time variation of coronal shocks themselves. Wave amplification by streaming protons in the largest SEP events suppresses the escape of ions with low A/Q, creating observed powers of A/Q from +1 to +3 upstream of the accelerating shock, decreasing to small negative powers downstream. Of course, the powers of A/Q are correlated with the spectral indices of He, O, and Fe, yet unexplained departures exist.


title: Video Panoptic Segmentation

abstract: Panoptic segmentation has become a new standard of visual recognition task by unifying previous semantic segmentation and instance segmentation tasks in concert. In this paper, we propose and explore a new video extension of this task, called video panoptic segmentation. The task requires generating consistent panoptic segmentation as well as an association of instance ids across video frames. To invigorate research on this new task, we present two types of video panoptic datasets. The first is a re-organization of the synthetic VIPER dataset into the video panoptic format to exploit its large-scale pixel annotations. The second is a temporal extension on the Cityscapes val. set, by providing new video panoptic annotations (Cityscapes-VPS). Moreover, we propose a novel video panoptic segmentation network (VPSNet) which jointly predicts object classes, bounding boxes, masks, instance id tracking, and semantic segmentation in video frames. To provide appropriate metrics for this task, we propose a video panoptic quality (VPQ) metric and evaluate our method and several other baselines. Experimental results demonstrate the effectiveness of the presented two datasets. We achieve state-of-the-art results in image PQ on Cityscapes and also in VPQ on Cityscapes-VPS and VIPER datasets. The datasets and code are made publicly available.