title: Cellular Fourier analysis for geometrically disordered materials

abstract: Many media are divided into elementary units with irregular shape and size, as exemplified by domains in magnetic materials, bubbles in foams, or cells in biological tissues. Such media are essentially characterized by geometrical disorder of their elementary units, which we term cells. Cells set a reference scale at which are often assessed parameters and fields reflecting material properties and state. Here, we consider the spectral analysis of spatially varying fields. Such analysis is difficult in geometrically disordered media, because space discretization based on standard coordinate systems is not commensurate with the natural discretization into geometrically disordered cells. Indeed, we found that two classical spectral methods, the Fast Fourier Transform and the Graph Fourier transform, fail to reproduce all expected properties of spectra of plane waves and of white noise. We therefore built a method, which we call Cellular Fourier Transform (CFT), to analyze cell-scale fields, which comprise both discrete fields defined only at cell level and continuous fields smoothed out from their sub-cell variations. Our approach is based on the construction of a discrete operator suited to the disordered geometry and on the computation of its eigenvectors, which, respectively, play the same role as the Laplace operator and sine waves in Euclidean coordinate systems. We show that CFT has the expected behavior for sinusoidal fields and for random fields with long-range correlations. Our approach for spectral analysis is suited to any geometrically disordered material, such as a biological tissue with complex geometry, opening the path to systematic multiscale analyses of material behavior.


title: A Contour Stochastic Gradient Langevin Dynamics Algorithm for Simulations of Multi-modal Distributions

abstract: We propose an adaptively weighted stochastic gradient Langevin dynamics algorithm (SGLD), so-called contour stochastic gradient Langevin dynamics (CSGLD), for Bayesian learning in big data statistics. The proposed algorithm is essentially a \emph{scalable dynamic importance sampler}, which automatically \emph{flattens} the target distribution such that the simulation for a multi-modal distribution can be greatly facilitated. Theoretically, we prove a stability condition and establish the asymptotic convergence of the self-adapting parameter to a {\it unique fixed-point}, regardless of the non-convexity of the original energy function; we also present an error analysis for the weighted averaging estimators. Empirically, the CSGLD algorithm is tested on multiple benchmark datasets including CIFAR10 and CIFAR100. The numerical results indicate its superiority to avoid the local trap problem in training deep neural networks.


title: Spread of Tweets in Climate Discussions

abstract: Characterising the spreading of ideas within echo chambers is essential for understanding polarisation. In this paper, we explore the characteristics of popular and viral content in climate change discussions on Twitter around the 2019 announcement of the Nobel Peace Prize, where we find the retweet network of users to be polarised into two well-separated groups of activists and sceptics. Operationalising popularity as the number of retweets and virality as the spreading probability inferred using an independent cascade model, we find that the viral themes echo and differ from the popular themes in interesting ways. Most importantly, we find that the most viral themes in the two groups reflect different types of bonds that tie the community together, yet both function to enhance ingroup connections while repulsing outgroup engagement. With this, our study sheds light, from an information spreading perspective, on the formation and upkeep of echo chambers in climate discussions.


title: A Differentiable Newton Euler Algorithm for Multi-body Model Learning

abstract: In this work, we examine a spectrum of hybrid model for the domain of multi-body robot dynamics. We motivate a computation graph architecture that embodies the Newton Euler equations, emphasizing the utility of the Lie Algebra form in translating the dynamical geometry into an efficient computational structure for learning. We describe the used virtual parameters that enable unconstrained physical plausible dynamics and the used actuator models. In the experiments, we define a family of 26 grey-box models and evaluate them for system identification of the simulated and physical Furuta Pendulum and Cartpole. The comparison shows that the kinematic parameters, required by previous white-box system identification methods, can be accurately inferred from data. Furthermore, we highlight that models with guaranteed bounded energy of the uncontrolled system generate non-divergent trajectories, while more general models have no such guarantee, so their performance strongly depends on the data distribution. Therefore, the main contributions of this work is the introduction of a white-box model that jointly learns dynamic and kinematics parameters and can be combined with black-box components. We then provide extensive empirical evaluation on challenging systems and different datasets that elucidates the comparative performance of our grey-box architecture with comparable white- and black-box models.


title: Adversarial Training for Code Retrieval with Question-Description Relevance Regularization

abstract: Code retrieval is a key task aiming to match natural and programming languages. In this work, we propose adversarial learning for code retrieval, that is regularized by question-description relevance. First, we adapt a simple adversarial learning technique to generate difficult code snippets given the input question, which can help the learning of code retrieval that faces bi-modal and data-scarce challenges. Second, we propose to leverage question-description relevance to regularize adversarial learning, such that a generated code snippet should contribute more to the code retrieval training loss, only if its paired natural language description is predicted to be less relevant to the user given question. Experiments on large-scale code retrieval datasets of two programming languages show that our adversarial learning method is able to improve the performance of state-of-the-art models. Moreover, using an additional duplicate question prediction model to regularize adversarial learning further improves the performance, and this is more effective than using the duplicated questions in strong multi-task learning baselines


title: Computational Aspects of Speed-Dependent Voigt and Rautian Profiles

abstract: For accurate line-by-line modeling of molecular cross sections several physical processes "beyond Voigt" have to be considered. For the speed-dependent Voigt and Rautian profiles (SDV, SDR) and the Hartmann-Tran profile the difference $w(i z_-)-w(i z_+)$ of two complex error functions (essentially Voigt functions) has to be evaluated where the function arguments $z_\pm$ are given by the sum and difference of two square roots. These two terms describing $z_\pm$ can be huge and the default implementation of the difference can lead to large cancellation errors. First we demonstrate that these problems can be avoided by a simple reformulation of $z_-$. Furthermore we show that a single rational approximation of the complex error function valid in the whole complex plane (e.g. by Humlicek, 1979 or Weideman, 1994) enables computation of the SDV and SDR with four significant digits or better. Our benchmarks indicate that the SDV and SDR functions are about a factor 2.2 slower compared to the Voigt function, but for evaluation of molecular cross sections this time lag does not significantly prolong the overall program execution because speed-dependent parameters are available only for a fraction of strong lines.


title: Number of points of certain Artin-Schreier curves

abstract: We prove a conjecture of Johansen, Helleseth and Kholosha concerning equality of exponential sums related to the cross-correlation of $m$-sequences. In the proof we show that certain Artin-Schreier curves have the same number of points over finite fields. This has a consequence regarding the L-polynomials of these curves.


title: Non-existence of a universal zero entropy system for non-periodic amenable group actions

abstract: Let $G$ be a non-periodic amenable group. We prove that there does not exist a topological action of $G$ for which the set of ergodic invariant measures coincides with the set of all ergodic measure-theoretic $G$-systems of entropy zero. Previously J. Serafin, answering a question by B. Weiss, proved the same for $G = \mathbb{Z}$.


title: Improving Prediction of Real-Time Loneliness and Companionship Type Using Geosocial Features of Personal Smartphone Data

abstract: Loneliness is a widely affecting mental health symptom and can be mediated by and co-vary with patterns of social exposure. Using momentary survey and smartphone sensing data collected from 129 Android-using college student participants over three weeks, we (1) investigate and uncover the relations between momentary loneliness experience and companionship type and (2) propose and validate novel geosocial features of smartphone-based Bluetooth and GPS data for predicting loneliness and companionship type in real time. We base our features on intuitions characterizing the quantity and spatiotemporal predictability of an individual's Bluetooth encounters and GPS location clusters to capture personal significance of social exposure scenarios conditional on their temporal distribution and geographic patterns. We examine our features' statistical correlation with momentary loneliness through regression analyses and evaluate their predictive power using a sliding window prediction procedure. Our features achieved significant performance improvement compared to baseline for predicting both momentary loneliness and companionship type, with the effect stronger for the loneliness prediction task. As such we recommend incorporation and further evaluation of our geosocial features proposed in this study in future mental health sensing and context-aware computing applications.


title: Imitation with Neural Density Models

abstract: We propose a new framework for Imitation Learning (IL) via density estimation of the expert's occupancy measure followed by Maximum Occupancy Entropy Reinforcement Learning (RL) using the density as a reward. Our approach maximizes a non-adversarial model-free RL objective that provably lower bounds reverse Kullback-Leibler divergence between occupancy measures of the expert and imitator. We present a practical IL algorithm, Neural Density Imitation (NDI), which obtains state-of-the-art demonstration efficiency on benchmark control tasks.


title: Parametrized topological complexity of collision-free motion planning in the plane

abstract: Parametrized motion planning algorithms have high degrees of universality and flexibility, as they are designed to work under a variety of external conditions, which are viewed as parameters and form part of the input of the underlying motion planning problem. In this paper, we analyze the parameterized motion planning problem for the motion of many distinct points in the plane, moving without collision and avoiding multiple distinct obstacles with a priori unknown positions. This complements our prior work [arXiv:2009.06023], where parameterized motion planning algorithms were introduced, and the obstacle-avoiding collision-free motion planning problem in three-dimensional space was fully investigated. The planar case requires different algebraic and topological tools than its spatial analog.


title: Connections between Relational Event Model and Inverse Reinforcement Learning for Characterizing Group Interaction Sequences

abstract: In this paper we explore previously unidentified connections between relational event model (REM) from the field of network science and inverse reinforcement learning (IRL) from the field of machine learning with respect to their ability to characterize sequences of directed social interaction events in group settings. REM is a conventional approach to tackle such a problem whereas the application of IRL is a largely unbeaten path. We begin by examining the mathematical components of both REM and IRL and find straightforward analogies between the two methods as well as unique characteristics of the IRL approach. We demonstrate the special utility of IRL in characterizing group social interactions with an empirical experiment, in which we use IRL to infer individual behavioral preferences based on a sequence of directed communication events from a group of virtual-reality game players interacting and cooperating to accomplish a shared goal. Our comparison and experiment introduce fresh perspectives for social behavior analytics and help inspire new research opportunities at the nexus of social network analysis and machine learning.


title: Pair production due to an electric field in 1+1 dimensions and the validity of the semiclassical approximation

abstract: Solutions to the backreaction equation in 1+1-dimensional semiclassical electrodynamics are obtained and analyzed when considering a time-varying homogeneous electric field initially generated by a classical electric current, coupled to either a quantized scalar field or a quantized spin-$\frac{1}{2}$ field. Particle production by way of the Schwinger effect leads to backreaction effects that modulate the electric field strength. Details of the particle production process are investigated along with the transfer of energy between the electric field and the particles. The validity of the semiclassical approximation is also investigated using a criterion previously implemented for chaotic inflation and, in an earlier form, semiclassical gravity. The criterion states that the semiclassical approximation will break down if any linearized gauge-invariant quantity constructed from solutions to the linear response equation, with finite nonsingular data, grows rapidly for some period of time. Approximations to homogeneous solutions of the linear response equation are computed and it is found that the criterion is violated when the maximum value, $E_{\rm max}$, obtained by the electric field is of the order of the critical scale for the Schwinger effect, $E_{\rm max} \sim E_{\rm crit}\equiv m^2/q$, where $m$ is the mass of the quantized field and $q$ is its electric charge. For these approximate solutions the criterion appears to be satisfied in the extreme limits $\frac{qE_{\rm max}}{m^2} \ll 1$ and $\frac{qE_{\rm max}}{m^2} \gg 1$.


title: Weakly Mixing Systems with Dense Prime Orbits

abstract: We show existence of smooth, weakly mixing reparametrizations of some linear flows on $\mathbb{T}^2$ for which all orbits sampled at prime times are dense.


title: Finite Temperature Auxiliary Field Quantum Monte Carlo in the Canonical Ensemble

abstract: Finite temperature auxiliary field-based Quantum Monte Carlo methods, including Determinant Quantum Monte Carlo (DQMC) and Auxiliary Field Quantum Monte Carlo (AFQMC), have historically assumed pivotal roles in the investigation of the finite temperature phase diagrams of a wide variety of multidimensional lattice models and materials. Despite their utility, however, these techniques are typically formulated in the grand canonical ensemble, which makes them difficult to apply to condensates like superfluids and difficult to benchmark against alternative methods that are formulated in the canonical ensemble. Working in the grand canonical ensemble is furthermore accompanied by the increased overhead associated with having to determine the chemical potentials that produce desired fillings. Given this backdrop, in this work, we present a new recursive approach for performing AFQMC simulations in the canonical ensemble that does not require knowledge of chemical potentials. To derive this approach, we exploit the convenient fact that AFQMC solves the many-body problem by decoupling many-body propagators into integrals over one-body problems to which non-interacting theories can be applied. We benchmark the accuracy of our technique on illustrative Bose and Fermi Hubbard models and demonstrate that it can converge more quickly to the ground state than grand canonical AFQMC simulations. We believe that our novel use of HS-transformed operators to implement algorithms originally derived for non-interacting systems will motivate the development of a variety of other methods and anticipate that our technique will enable direct performance comparisons against other many-body approaches formulated in the canonical ensemble.


title: Social Hierarchy-based Distributed Economic Model Predictive Control of Floating Offshore Wind Farms

abstract: This paper implements a recently developed social hierarchy-based distributed economic model predictive control (DEMPC) algorithm in floating offshore wind farms for the purpose of power maximization. The controller achieves this objective using the concept of yaw and induction-based turbine repositioning (YITuR), which minimizes the overlap areas between adjacent floating wind turbine rotors in real-time to minimize the wake effect. Floating wind farm dynamics and performance are predicted numerically using FOWFSim-Dyn. To ensure fast decision-making by the DEMPC algorithm, feed-forward neural networks are used to estimate floating wind turbine dynamics during the process of dynamic optimization. For simulated wind farms with layouts ranging from 1-by-2 to 1-by-5, an increase of 20% in energy production is predicted when using YITuR instead of greedy operation. Increased variability in wind speed and direction is also studied and is shown to diminish controller performance due to rising errors in neural network predictions.


title: Free Energy of a Knotted Polymer Confined to Narrow Cylindrical and Conical Channels

abstract: Monte Carlo simulations are used to study the conformational behavior of a semiflexible polymer confined to cylindrical and conical channels. The channels are sufficiently narrow that the conditions for the Odijk regime are marginally satisfied. For cylindrical confinement, we examine polymers with a single knot of topology $3_1$, $4_1$, or $5_1$, as well as unknotted polymers that are capable of forming S-loops. We measure the variation of the free energy $F$ with the end-to-end polymer extension length $X$ and examine the effect of varying the polymer topology, persistence length $P$ and cylinder diameter $D$ on the free energy functions. Similarly, we characterize the behavior of the knot span along the channel. We find that increasing the knot complexity increases the typical size of the knot. In the regime of low $X$, where the knot/S-loop size is large, the conformational behavior is independent of polymer topology. In addition, the scaling properties of the free energy and knot span are in agreement with predictions from a theoretical model constructed using known properties of interacting polymers in the Odijk regime. We also examine the variation of $F$ with position of a knot in conical channels for various values of the cone angle $\alpha$. The free energy decreases as the knot moves in a direction where the cone widens, and it also decreases with increasing $\alpha$ and with increasing knot complexity. The behavior is in agreement with predictions from a theoretical model in which the dominant contribution to the change in $F$ is the change in the size of the hairpins as the knot moves to the wider region of the channel.


title: Essential self-adjointness of symmetric first-order differential systems and confinement of Dirac particles on bounded domains in $\mathbb{R}^d$

abstract: We prove essential self-adjointness of Dirac operators with Lorentz scalar potentials which grow sufficiently fast near the boundary $\partial\Omega$ of the spatial domain $\Omega\subset\mathbb R^d$. On the way, we first consider general symmetric first order differential systems, for which we identify a new, large class of potentials, called scalar potentials, ensuring essential self-adjointness. Furthermore, using the supersymmetric structure of the Dirac operator in the two dimensional case, we prove confinement of Dirac particles, i.e. essential self-adjointness of the operator, solely by magnetic fields $\mathcal{B}$ assumed to grow, near $\partial\Omega$, faster than $1/\big(2\text{dist} (x, \partial\Omega)^2\big)$.


title: Strain Relaxation via Phase Transformation in SrSnO3

abstract: SrSnO3 (SSO) is an emerging ultra-wide bandgap (UWBG) semiconductor with potential for highpower applications. In-plane compressive strain was recently shown to stabilize the high temperature tetragonal phase of SSO at room temperature (RT) which exists at T > 1062 K in bulk. Here, we report on the study of strain relaxation in epitaxial, tetragonal phase of Nd-doped SSO films grown on GdScO3 (110) (GSO) substrates using radical-based hybrid molecular beam epitaxy. The thinnest SSO film (thickness, t = 12 nm) yielded a fully coherent tetragonal phase at RT. At 12 nm < t < 110 nm, the tetragonal phase first transformed into orthorhombic phase and then at t > 110 nm, the orthorhombic phase began to relax by forming misfit dislocations. Remarkably, the tetragonal phase remained fully coherent until it completely transformed into the orthorhombic phase. Using thickness- and temperature-dependent electronic transport measurements, we discuss the important roles of the surface, phase coexistence, and misfit dislocations on carrier density and mobility in Nd-doped SSO. This study provides unprecedented insights into the strain relaxation behavior and its consequences for electronic transport in doped SSO with implications in the development of high-power electronic devices.


title: Online Active Model Selection for Pre-trained Classifiers

abstract: Given $k$ pre-trained classifiers and a stream of unlabeled data examples, how can we actively decide when to query a label so that we can distinguish the best model from the rest while making a small number of queries? Answering this question has a profound impact on a range of practical scenarios. In this work, we design an online selective sampling approach that actively selects informative examples to label and outputs the best model with high probability at any round. Our algorithm can be used for online prediction tasks for both adversarial and stochastic streams. We establish several theoretical guarantees for our algorithm and extensively demonstrate its effectiveness in our experimental studies.


title: Comparative Analysis of Control Barrier Functions and Artificial Potential Fields for Obstacle Avoidance

abstract: Artificial potential fields (APFs) and their variants have been a staple for collision avoidance of mobile robots and manipulators for almost 40 years. Its model-independent nature, ease of implementation, and real-time performance have played a large role in its continued success over the years. Control barrier functions (CBFs), on the other hand, are a more recent development, commonly used to guarantee safety for nonlinear systems in real-time in the form of a filter on a nominal controller. In this paper, we address the connections between APFs and CBFs. At a theoretic level, we prove that APFs are a special case of CBFs: given a APF one obtains a CBFs, while the converse is not true. Additionally, we prove that CBFs obtained from APFs have additional beneficial properties and can be applied to nonlinear systems. Practically, we compare the performance of APFs and CBFs in the context of obstacle avoidance on simple illustrative examples and for a quadrotor, both in simulation and on hardware using onboard sensing. These comparisons demonstrate that CBFs outperform APFs.


title: Solving relativistic three-body integral equations in the presence of bound states

abstract: We present a systematically improvable method for numerically solving relativistic three-body integral equations for the partial-wave projected amplitudes. The method consists of a discretization procedure in momentum space, which approximates the continuum problem with a matrix equation. It is solved for different matrix sizes, and in the end, an extrapolation is employed to restore the continuum limit. Our technique is tested by solving a three-body problem of scalar particles with an $S$ wave two-body bound state. We discuss two methods of incorporating the pole contribution in the integral equations, both of them leading to agreement with previous results obtained using finite-volume spectra of the same theory. We provide an analytic and numerical estimate of the systematic errors. Although we focus on kinematics below the three-particle threshold, we provide numerical evidence that the methods presented allow for determination of amplitude above this threshold as well.


title: Isotropy and Combination Problems

abstract: In a previous paper, the author and his collaborators studied the phenomenon of isotropy in the context of single-sorted equational theories, and showed that the isotropy group of the category of models of any such theory encodes a notion of inner automorphism for the theory. Using results from the treatment of combination problems in term rewriting theory, we show in this article that if $\mathbb{T}_1$ and $\mathbb{T}_2$ are (disjoint) equational theories satisfying minimal assumptions, then any free, finitely generated model of the disjoint union theory $\mathbb{T}_1 + \mathbb{T}_2$ has trivial isotropy group, and hence the only inner automorphisms of such models, i.e. the only automorphisms of such models that are coherently extendible, are the identity automorphisms. As a corollary, we show that the global isotropy group of the category of models $(\mathbb{T}_1 + \mathbb{T}_2)\mathsf{mod}$, i.e. the group of invertible elements of the centre of this category, is the trivial group.


title: Is the new model better? One metric says yes, but the other says no. Which metric do I use?

abstract: Incremental value (IncV) evaluates the performance change from an existing risk model to a new model. It is one of the key considerations in deciding whether a new risk model performs better than the existing one. Problems arise when different IncV metrics contradict each other. For example, compared with a prescribed-dose model, an ovarian-dose model for predicting acute ovarian failure has a slightly lower area under the receiver operating characteristic curve (AUC) but increases the area under the precision-recall curve (AP) by 48%. This phenomenon of conflicting conclusions is not uncommon, and it creates a dilemma in medical decision making. In this article, we examine the analytical connections and differences between two IncV metrics: IncV in AUC (IncV-AUC) and IncV in AP (IncV-AP). Additionally, since they are both semi-proper scoring rules, we compare them with a strictly proper scoring rule: the IncV of the scaled Brier score (IncV-sBrS), via a numerical study. We demonstrate that both IncV-AUC and IncV-AP are weighted averages of the changes (from the existing model to the new one) in separating the risk score distributions between events and non-events. However, IncV-AP assigns heavier weights to the changes in the high-risk group, whereas IncV-AUC weights the changes equally. In the numerical study, we find that IncV-AP has a wide range, from negative to positive, but the size of IncV-AUC is much smaller. In addition, IncV-AP and IncV-sBr Sare highly consistent, but IncV-AUC is negatively correlated with IncV-sBrS and IncV-AP at a low event rate. IncV-AUC and IncV-AP are the least consistent among the three pairs, and their differences are more pronounced as the event rate decreases.


title: Calibration of the MaGIXS experiment I: Calibration of the X-ray source at the X-ray and Cryogenic Facility (XRCF)

abstract: The Marshall Grazing Incidence Spectrometer {\it MaGIXS} is a sounding rocket experiment that will observe the soft X-ray spectrum of the Sun from 24 - 6.0 \AA\ (0.5 - 2.0 keV) and is scheduled for launch in 2021. Component and instrument level calibrations for the {\it MaGIXS} instrument are carried out using the X-ray and Cryogenic Facility (XRCF) at NASA Marshall Space Flight Center. In this paper, we present the calibration of the incident X-ray flux from the electron impact source with different targets at the XRCF using a CCD camera; the photon flux at the CCD was low enough to enable its use as a "photon counter" i.e. the ability to identify individual photon hits and calculate their energy. The goal of this paper is two-fold: 1) to confirm that the flux measured by the XRCF beam normalization detectors is consistent with the values reported in the literature and therefore reliable for {\it MaGIXS} calibration and 2) to develop a method of counting photons in CCD images that best captures their number and energy


title: Multi-Objective Bayesian Optimization for Accelerator Tuning

abstract: Particle accelerators require constant tuning during operation to meet beam quality, total charge and particle energy requirements for use in a wide variety of physics, chemistry and biology experiments. Maximizing the performance of an accelerator facility often necessitates multi-objective optimization, where operators must balance trade-offs between multiple objectives simultaneously, often using limited, temporally expensive beam observations. Usually, accelerator optimization problems are solved offline, prior to actual operation, with advanced beamline simulations and parallelized optimization methods (NSGA-II, Swarm Optimization). Unfortunately, it is not feasible to use these methods for online multi-objective optimization, since beam measurements can only be done in a serial fashion, and these optimization methods require a large number of measurements to converge to a useful solution.Here, we introduce a multi-objective Bayesian optimization scheme, which finds the full Pareto front of an accelerator optimization problem efficiently in a serialized manner and is thus a critical step towards practical online multi-objective optimization in accelerators.This method uses a set of Gaussian process surrogate models, along with a multi-objective acquisition function, which reduces the number of observations needed to converge by at least an order of magnitude over current methods.We demonstrate how this method can be modified to specifically solve optimization challenges posed by the tuning of accelerators.This includes the addition of optimization constraints, objective preferences and costs related to changing accelerator parameters.


title: An evidence-based assumption that helps to reduce the discrepancy between the observed and predicted $^7$Be abundances in novae

abstract: Recent spectroscopic measurements of the equivalent widths of the resonant Be II doublet and Ca II K lines and their ratios in expanding nova ejecta indicate surprisingly high abundances of $^7$Be with a typical mass fraction $X_\mathrm{obs}(^7\mathrm{Be}) = 10^{-4}$. This is an order of magnitude larger than theoretically predicted values of $X_\mathrm{theor}(^7\mathrm{Be})\sim 10^{-5}$ for novae. We use an analytical solution of the $^7$Be production equations to demonstrate that $X_\mathrm{theor}(^7\mathrm{Be})$ is proportional to the $^4$He mass fraction $Y$ in the nova accreted envelope and then we perform computations of 1D hydrostatic evolution of the $1.15\,M_\odot$ CO nova model that confirm our conclusion based on the analytical solution. Our assumption of enhanced $^4$He abundances helps to reduce, although not completely eliminate, the discrepancy between $X_\mathrm{obs}(^7\mathrm{Be})$ and $X_\mathrm{theor}(^7\mathrm{Be})$. It is supported by UV, optical and IR spectroscopy data that reveal unusually high values of $Y$ in nova ejecta. We also show that a significantly increased abundance of $^3$He in nova accreted envelopes does not lead to higher values of $X_\mathrm{theor}(^7\mathrm{Be})$ because this assumption affects the evolution of nova models resulting in a decrease of both their peak temperatures and accreted masses and, as a consequence, in a reduced production of $^7$Be.


title: Lattice contribution to the unconventional charge density wave transition in $2H$-NbSe$_2$: a non-equilibrium optical approach

abstract: The complex Fermi surfaces of transition-metal dichalcogenides (TMDCs) challenge the standard Peierls-instability-driven charge-density-wave (CDW) formation. Recently, evidence has been accumulating of a prominent role of ionic thermal fluctuations, which frozen out below $T_{CDW}$ inducing a periodic lattice distortion (PLD). We focus on $2H$-NbSe$_2$, displaying a quasi-commensurate CDW below $T_{CDW}$ $\approx$33 K, and use time-resolved optical spectroscopy (TR-OS) to detect and disentangle the electronic and lattice degrees of freedom. We reveal a fingerprint of the ordered phase at h$\nu$ $\sim$1 eV: at $T_{CDW}$, a divergent relaxation timescale and a sign-change of the differential reflectivity indicate that CDW gap opening and PLD formation occur at the same temperature. However, we show that these effects can be decoupled under moderate photoexcitation, forming a long-lived state in which the electronic order is destroyed, but the lattice distortion is not. Our results and computations suggest an unconventional CDW mechanims in $2H$-NbSe$_2$, highlighting the dominant role of the lattice in driving the ordered-phase formation.


title: On the Shape Fields Finiteness Principle

abstract: In this paper, we improve the finiteness constant for the finiteness principles for $C^m(\mathbb{R}^n,\mathbb{R}^d)$ and $C^{m-1,1}(\mathbb{R}^n,\mathbb{R}^D)$ selection proven by Fefferman, Israel, and the second author and extend the more general shape fields finiteness principle to the vector-valued case.


title: Cross-Lingual Transfer in Zero-Shot Cross-Language Entity Linking

abstract: Cross-language entity linking grounds mentions in multiple languages to a single-language knowledge base. We propose a neural ranking architecture for this task that uses multilingual BERT representations of the mention and the context in a neural network. We find that the multilingual ability of BERT leads to robust performance in monolingual and multilingual settings. Furthermore, we explore zero-shot language transfer and find surprisingly robust performance. We investigate the zero-shot degradation and find that it can be partially mitigated by a proposed auxiliary training objective, but that the remaining error can best be attributed to domain shift rather than language transfer.


title: Synthetic high-order PT symmetry in a single coil resonator

abstract: The exploration of non-Hermitian systems with parity-time (PT) symmetry has witnessed immense research interest both fundamentally and technologically in a wide range of subject areas in physics and engineering. One significant example of the principal emerging fields in this context is the PT symmetric wireless applications using multiple coils that are spatially separated but mutually coupled with position-dependent coupling strength. Such a spatial PT configuration limits the flexibility and miniaturization of the PT symmetric designs. As far as this is concerned, inspired by scattering induced two opposite whispering-gallery (WG) modes in an optical resonator, analogously here we experimentally demonstrate a specially constructed second-order (2-nd order) PT symmetry in a single coil resonator, whose currents with two different directions are excited by internal bypass capacitor. Our proposed structure has the following peculiar feature: First, the bypass capacitor induces coupling in spectral resonances allow us to observe a 2-nd order phase transition between symmetry regimes, without the need of a second coil in the spatial PT case. Under this circumstance, this specially constructed PT symmetry can be regarded as synthetic PT symmetry, which is enabled by coupling modes with different directions. Second, by introducing two or more internal bypass capacitors, the synthetic high-order PT symmetric system bearing such as third-order exceptional point in a single coil resonator can be realized. These results will provide a new paradigm to realize higher-order PT symmetry towards the investigation of non-Hermitian physics in a synthetic perspective, which can be extended to other physical platforms such as optics and acoustics.


title: Remarks on multivariate Gaussian Process

abstract: Gaussian processes occupy one of the leading places in modern statistics and probability theory due to their importance and a wealth of strong results. The common use of Gaussian processes is in connection with problems related to estimation, detection, and many statistical or machine learning models. With the fast development of Gaussian process applications, it is necessary to consolidate the fundamentals of vector-valued stochastic processes, in particular multivariate Gaussian processes, which is the essential theory for many applied problems with multiple correlated responses. In this paper, we propose a precise definition of multivariate Gaussian processes based on Gaussian measures on vector-valued function spaces, and provide an existence proof. In addition, several fundamental properties of multivariate Gaussian processes, such as strict stationarity and independence, are introduced. We further derive multivariate Brownian motion including It\^o lemma as a special case of a multivariate Gaussian process, and present a brief introduction to multivariate Gaussian process regression as a useful statistical learning method for multi-output prediction problems.


title: Miniaturized backward coupler realized by the circuit-based magnetic hyperbolic waveguide

abstract: Planar waveguides can limit the transmission of electromagnetic waves in a specific direction and have a wide range of applications in filters, sensors, and energy-transfer devices. However, given the increasing demand for planar integrated photonics, new waveguides are required with excellent characteristics such as more functionality, greater efficiency, and smaller size. In this work, we report the experimental results for a subwavelength planar microwave-regime waveguide fabricated from circuit-based magnetic hyperbolic metamaterial (HMM) and develop the associated theory. HMM is a special type of anisotropic metamaterial whose isofrequency contour (IFC) takes the form of an open hyperboloid. Because of the open-dispersion IFCs, HMMs support propagating high-k modes with large effective refractive indices, which allows planar hyperbolic waveguides to be miniaturized. In particular, as opposed to the traditional dielectric slab waveguide, the group velocity of the hyperbolic guided modes is negative-a characteristic that can be exploited to realize the backward propagation of electromagnetic waves. Furthermore, we study the abnormal photonic spin Hall effect and experimentally demonstrate an ultracompact backward coupler based on hyperbolic waveguides. The experimental results are consistent with numerical simulations. This work not only reveals the significant potential of circuit-based platforms for the experimental study of the propagation and coupling of guided modes but also promotes the use of HMMs in the microwave regime for numerous integrated functional devices.


title: Dream and Search to Control: Latent Space Planning for Continuous Control

abstract: Learning and planning with latent space dynamics has been shown to be useful for sample efficiency in model-based reinforcement learning (MBRL) for discrete and continuous control tasks. In particular, recent work, for discrete action spaces, demonstrated the effectiveness of latent-space planning via Monte-Carlo Tree Search (MCTS) for bootstrapping MBRL during learning and at test time. However, the potential gains from latent-space tree search have not yet been demonstrated for environments with continuous action spaces. In this work, we propose and explore an MBRL approach for continuous action spaces based on tree-based planning over learned latent dynamics. We show that it is possible to demonstrate the types of bootstrapping benefits as previously shown for discrete spaces. In particular, the approach achieves improved sample efficiency and performance on a majority of challenging continuous-control benchmarks compared to the state-of-the-art.


title: Note on local mixing techniques for stochastic differential equations

abstract: This paper discusses several techniques which may be used for applying the coupling method to solutions of stochastic differential equations (SDEs). They all work in dimension $d\ge 1$, although, in $d=1$ the most natural way is to use intersections of trajectories, which requires nothing but strong Markov property and non-degeneracy of the diffusion coefficient. In dimensions $d>1$ it is possible to use embedded Markov chains either by considering discrete times $n=0,1,\ldots$, or by arranging special stopping time sequences and to use local Markov -- Dobrushin's (MD) condition. Further applications may be based on one or another version of the MD condition. For studies of convergence and mixing rates the (Markov) process must be strong Markov and recurrent; however, recurrence is a separate issue which is not discussed in this paper.


title: Topology-Aware Graph Pooling Networks

abstract: Pooling operations have shown to be effective on computer vision and natural language processing tasks. One challenge of performing pooling operations on graph data is the lack of locality that is not well-defined on graphs. Previous studies used global ranking methods to sample some of the important nodes, but most of them are not able to incorporate graph topology. In this work, we propose the topology-aware pooling (TAP) layer that explicitly considers graph topology. Our TAP layer is a two-stage voting process that selects more important nodes in a graph. It first performs local voting to generate scores for each node by attending each node to its neighboring nodes. The scores are generated locally such that topology information is explicitly considered. In addition, graph topology is incorporated in global voting to compute the importance score of each node globally in the entire graph. Altogether, the final ranking score for each node is computed by combining its local and global voting scores. To encourage better graph connectivity in the sampled graph, we propose to add a graph connectivity term to the computation of ranking scores. Results on graph classification tasks demonstrate that our methods achieve consistently better performance than previous methods.


title: Conditions on the regularity of balanced $c$-partite tournaments for the existence of strong subtournaments with high minimum degree

abstract: We consider the following problem posed by Volkmann in 2007: How close to regular must a c-partite tournament be, to secure a strongly connected subtournament of order $c$? We give sufficient conditions on the regularity of balanced $c$-partite tournaments to assure the existence of strong maximal subtournament with minimum degree at least $\left\lfloor \frac{c-2}{4}\right\rfloor+1$. We obtain this result as an application of counting the number of subtournaments of order $c$ for which a vertex has minimum out-degree (resp. in-degree) at most $q\geq 0$.


title: Tuning resonance energy transfer with magneto-optical properties of graphene

abstract: We investigate the resonance energy transfer (RET) rate between two quantum emitters near a suspended graphene sheet in vacuum under the influence of an external magnetic field. We perform the analysis for low and room temperatures and show that, due to the extraordinary magneto-optical response of graphene, it allows for an active control and tunability of the RET even in the case of room temperature. We also demonstrate that the RET rate is extremely sensitive to small variations of the applied magnetic field, and can be tuned up to a striking six orders of magnitude for quite realistic values of magnetic field. Moreover, we evidence the fundamental role played by the magnetoplasmon polaritons supported by the graphene monolayer as the dominant channel for the RET within a certain distance range. Our results suggest that magneto-optical media may take the manipulation of energy transfer between quantum emitters to a whole new level, and broaden even more its great spectrum of applications.


title: Isotropy Groups of Free Racks and Quandles

abstract: In this article, we characterize the (covariant) isotropy groups of free, finitely generated racks and quandles. As a consequence, we show that the usual inner automorphisms of such racks and quandles are precisely those automorphisms that are "coherently extendible". We then use this result to compute the global isotropy groups of the categories of racks and quandles, i.e. the automorphism groups of the identity functors of these categories.


title: Open quantum systems beyond Fermi's golden rule: Diagrammatic expansion of the steady-state time-convolutionless master equation

abstract: Steady-state observables, such as occupation numbers and currents, are crucial experimental signatures in open quantum systems. The time-convolutionless (TCL) master equation, which is both exact and time-local, is an ideal candidate for the perturbative computation of such observables. We develop a diagrammatic approach to evaluate the steady-state TCL generator based on operators rather than superoperators. We obtain the steady-state occupation numbers, extend our formulation to the calculation of currents, and provide a simple physical interpretation of the diagrams. We further benchmark our method on a single non-interacting level coupled to Fermi reservoirs, where we recover the exact expansion to next-to-leading order. The low number of diagrams appearing in our formulation makes the extension to higher orders accessible. Combined, these properties make the steady-state time-convolutionless master equation an effective tool for the calculation of steady-state properties in open quantum systems.


title: New Properties of the Data Distillation Method When Working With Tabular Data

abstract: Data distillation is the problem of reducing the volume oftraining data while keeping only the necessary information. With thispaper, we deeper explore the new data distillation algorithm, previouslydesigned for image data. Our experiments with tabular data show thatthe model trained on distilled samples can outperform the model trainedon the original dataset. One of the problems of the considered algorithmis that produced data has poor generalization on models with differenthyperparameters. We show that using multiple architectures during distillation can help overcome this problem.


title: Domain Wall Constraints on Two Higgs Doublet Models with $Z_2$ Symmetry

abstract: The Two Higgs Doublet Model (2HDM) with spontaneously broken $Z_2$ symmetry predicts a production of domain walls at the electroweak scale. We derive cosmological constraints on model parameters for both Type-I and Type-II 2HDMs from the requirement that domain walls do not dominate the Universe by the present day. For Type-I 2HDMs, we deduce the lower bound on the key parameter $\tan\beta > 10^5$ for a wide range of Higgs-boson masses $\sim$ 100 GeV or greater close to the Standard Model alignment limit. In addition, we perform numerical simulations of the 2HDM with an approximate as well as an exact $Z_2$ symmetry but biased initial conditions. In both cases, we find that domain wall networks are unstable and, hence, do not survive at late times. The domain walls experience an exponential suppression of scaling in these models which can help ameliorate the stringent constraints found in the case of an exact discrete symmetry. For a 2HDM with softly-broken $Z_2$ symmetry, we relate the size of this exponential suppression to the soft-breaking bilinear parameter $m_{12}$ allowing limits to be placed on this parameter of order $\mu$eV, such that domain wall domination can be avoided. In particular, for Type-II 2HDMs, we obtain a corresponding lower limit on the CP-odd phase $\theta$ generated by QCD instantons, $\theta \ \stackrel{>}{{}_\sim}\ 10^{-11}/(\sin\beta \cos\beta)$, which is in some tension with the upper limit of $\theta \ \stackrel{<}{{}_\sim}\ 10^{-11}$--$10^{-10}$, as derived from the non-observation of a non-zero neutron electric dipole moment. For a $Z_2$-symmetric 2HDM with biased initial conditions, we are able to relate the size of the exponential suppression to a biasing parameter $\varepsilon$ so as to avoid domain wall domination.


title: Counting Dimensions of Tangent Spaces to Hilbert Schemes of Points

abstract: In this paper we prove two results which further classify smoothness properties of Hilbert schemes of points. This is done by counting classes of arrows on Young diagrams corresponding to monomial ideals, building on the approach taken by Jan Cheah to show smoothness in the 2 dimensional case. We prove sufficient conditions for points to be smooth on Hilbert schemes of points in 3 dimensions and on nested schemes in 2 dimensions in terms of the geometry of the Young diagram. In particular, we proved that when the region between the two diagrams at a point of the nested scheme is rectangular, the corresponding point is smooth. We also proved that if the three dimensional Young diagram at a point can be oriented such that its horizontal layers are rectangular, then the point is smooth.


title: Robot Design With Neural Networks, MILP Solvers and Active Learning

abstract: Central to the design of many robot systems and their controllers is solving a constrained blackbox optimization problem. This paper presents CNMA, a new method of solving this problem that is conservative in the number of potentially expensive blackbox function evaluations; allows specifying complex, even recursive constraints directly rather than as hard-to-design penalty or barrier functions; and is resilient to the non-termination of function evaluations. CNMA leverages the ability of neural networks to approximate any continuous function, their transformation into equivalent mixed integer linear programs (MILPs) and their optimization subject to constraints with industrial strength MILP solvers. A new learning-from-failure step guides the learning to be relevant to solving the constrained optimization problem. Thus, the amount of learning is orders of magnitude smaller than that needed to learn functions over their entire domains. CNMA is illustrated with the design of several robotic systems: wave-energy propelled boat, lunar lander, hexapod, cartpole, acrobot and parallel parking. These range from 6 real-valued dimensions to 36. We show that CNMA surpasses the Nelder-Mead, Gaussian and Random Search optimization methods against the metric of number of function evaluations.


title: What breach? Measuring online awareness of security incidents by studying real-world browsing behavior

abstract: Awareness about security and privacy risks is important for developing good security habits. Learning about real-world security incidents and data breaches can alert people to the ways in which their information is vulnerable online, thus playing a significant role in encouraging safe security behavior. This paper examines 1) how often people read about security incidents online, 2) of those people, whether and to what extent they follow up with an action, e.g., by trying to read more about the incident, and 3) what influences the likelihood that they will read about an incident and take some action. We study this by quantitatively examining real-world internet-browsing data from 303 participants. Our findings present a bleak view of awareness of security incidents. Only 16% of participants visited any web pages related to six widely publicized large-scale security incidents; few read about one even when an incident was likely to have affected them (e.g., the Equifax breach almost universally affected people with Equifax credit reports). We further found that more severe incidents as well as articles that constructively spoke about the incident inspired more action. We conclude with recommendations for specific future research and for enabling useful security incident information to reach more people.


title: Degenerate solutions to the Dirac equation for massive particles and their applications in quantum tunneling

abstract: In a recent work we have proven the existence of degenerate solutions to the Dirac equation, corresponding to an infinite number of different electromagnetic fields, providing also some examples regarding massless particles. In the present article our results are extended significantly, providing degenerate solutions to the Dirac equation for particles with arbitrary mass, which, under certain conditions, could be interpreted as pairs of particles (or antiparticles) moving in a potential barrier with energy equal to the height of the barrier and spin opposite to each other. We calculate the electromagnetic fields corresponding to these solutions, providing also some examples regarding both spatially constant electromagnetic fields and electromagnetic waves. Further, we discuss some potential applications of our work, mainly regarding the control of the particles outside the potential barrier, without affecting their state inside the barrier. Finally, we study the effect of small perturbations to the degenerate solutions, showing that our results are still valid, in an approximate sense, provided that the amplitude of the electromagnetic fields corresponding to the exact degenerate solutions is sufficiently small.


title: Criniferous entire maps with absorbing Cantor bouquets

abstract: It is known that, for many transcendental entire functions in the Eremenko-Lyubich class $\mathcal{B}$, every escaping point can eventually be connected to infinity by a curve of escaping points. When this is the case, we say that the functions are criniferous. In this paper, we extend this result to a new class of maps in $\mathcal{B}$. Furthermore, we show that if a map belongs to this class, then its Julia set contains a Cantor bouquet; in other words, it is a subset of $\mathbb{C}$ ambiently homeomorphic to a straight brush.


title: Degenerate solutions to the massless Dirac and Weyl equations and a proposed method for controlling the quantum state of Weyl particles

abstract: In a recent work, we have shown that all solutions to the Weyl equation and a special class of solutions to the Dirac equation are degenerate in the sense that they remain unaltered under the influence of a wide variety of different electromagnetic fields. In this study, our previous work is significantly extended, providing a wide class of degenerate solutions to the Dirac equation for massless particles. The electromagnetic fields corresponding to these solutions are calculated and examples regarding both spatially constant electromagnetic fields and electromagnetic waves are also provided. Furthermore, some general solutions to the Weyl equation are presented, and the corresponding electromagnetic fields are calculated. Based on these results, a method for fully controlling the quantum state of Weyl particles using appropriate electromagnetic fields is proposed. Finally, the transition from degenerate to non-degenerate solutions as the particles acquire mass is discussed.


title: Idle Vehicle Relocation Strategy through Deep Learning for Shared Autonomous Electric Vehicle System Optimization

abstract: In optimization of a shared autonomous electric vehicle (SAEV) system, idle vehicle relocation strategies are important to reduce operation costs and customers' wait time. However, for an on-demand service, continuous optimization for idle vehicle relocation is computationally expensive, and thus, not effective. This study proposes a deep learning-based algorithm that can instantly predict the optimal solution to idle vehicle relocation problems under various traffic conditions. The proposed relocation process comprises three steps. First, a deep learning-based passenger demand prediction model using taxi big data is built. Second, idle vehicle relocation problems are solved based on predicted demands, and optimal solution data are collected. Finally, a deep learning model using the optimal solution data is built to estimate the optimal strategy without solving relocation. In addition, the proposed idle vehicle relocation model is validated by applying it to optimize the SAEV system. We present an optimal service system including the design of SAEV vehicles and charging stations. Further, we demonstrate that the proposed strategy can drastically reduce operation costs and wait times for on-demand services.


title: Critical Scaling of Compression-Driven Jamming of Athermal Frictionless Spheres in Suspension

abstract: We study numerically a system of athermal, overdamped, frictionless spheres, as in a non-Brownian suspension, in two and three dimensions. Compressing the system isotropically at a fixed rate $\dot\epsilon$, we investigate the critical behavior at the jamming transition. The finite compression rate introduces a control timescale, which allows one to probe the critical timescale associated with jamming. As was found previously for steady-state shear-driven jamming, we find for compression-driven jamming that pressure obeys a critical scaling relation as a function of packing fraction $\phi$ and compression rate $\dot\epsilon$, and that the bulk viscosity $p/\dot\epsilon$ diverges upon jamming. A scaling analysis determines the critical exponents associated with the compression-driven jamming transition. Our results suggest that stress-isotropic, compression-driven, jamming may be in the same universality class as stress-anisotropic, shear-driven, jamming.


title: Facial Emotion Recognition with Noisy Multi-task Annotations

abstract: Human emotions can be inferred from facial expressions. However, the annotations of facial expressions are often highly noisy in common emotion coding models, including categorical and dimensional ones. To reduce human labelling effort on multi-task labels, we introduce a new problem of facial emotion recognition with noisy multi-task annotations. For this new problem, we suggest a formulation from the point of joint distribution match view, which aims at learning more reliable correlations among raw facial images and multi-task labels, resulting in the reduction of noise influence. In our formulation, we exploit a new method to enable the emotion prediction and the joint distribution learning in a unified adversarial learning game. Evaluation throughout extensive experiments studies the real setups of the suggested new problem, as well as the clear superiority of the proposed method over the state-of-the-art competing methods on either the synthetic noisy labeled CIFAR-10 or practical noisy multi-task labeled RAF and AffectNet. The code is available at https://github.com/sanweiliti/noisyFER.


title: The Anatomical Edutainer

abstract: Physical visualizations (i.e., data representations by means of physical objects) have been used for many centuries in medical and anatomical education. Recently, 3D printing techniques started also to emerge. Still, other medical physicalizations that rely on affordable and easy-to-find materials are limited, while smart strategies that take advantage of the optical properties of our physical world have not been thoroughly investigated. We propose the Anatomical Edutainer, a workflow to guide the easy, accessible, and affordable generation of physicalizations for tangible, interactive anatomical edutainment. The Anatomical Edutainer supports 2D printable and 3D foldable physicalizations that change their visual properties (i.e., hues of the visible spectrum) under colored lenses or colored lights, to reveal distinct anatomical structures through user interaction.


title: Can I Trust My Fairness Metric? Assessing Fairness with Unlabeled Data and Bayesian Inference

abstract: We investigate the problem of reliably assessing group fairness when labeled examples are few but unlabeled examples are plentiful. We propose a general Bayesian framework that can augment labeled data with unlabeled data to produce more accurate and lower-variance estimates compared to methods based on labeled data alone. Our approach estimates calibrated scores for unlabeled examples in each group using a hierarchical latent variable model conditioned on labeled examples. This in turn allows for inference of posterior distributions with associated notions of uncertainty for a variety of group fairness metrics. We demonstrate that our approach leads to significant and consistent reductions in estimation error across multiple well-known fairness datasets, sensitive attributes, and predictive models. The results show the benefits of using both unlabeled data and Bayesian inference in terms of assessing whether a prediction model is fair or not.


title: Evaluating the Cost of Atomic Operations on Modern Architectures

abstract: Atomic operations (atomics) such as Compare-and-Swap (CAS) or Fetch-and-Add (FAA) are ubiquitous in parallel programming. Yet, performance tradeoffs between these operations and various characteristics of such systems, such as the structure of caches, are unclear and have not been thoroughly analyzed. In this paper we establish an evaluation methodology, develop a performance model, and present a set of detailed benchmarks for latency and bandwidth of different atomics. We consider various state-of-the-art x86 architectures: Intel Haswell, Xeon Phi, Ivy Bridge, and AMD Bulldozer. The results unveil surprising performance relationships between the considered atomics and architectural properties such as the coherence state of the accessed cache lines. One key finding is that all the tested atomics have comparable latency and bandwidth even if they are characterized by different consensus numbers. Another insight is that the hardware implementation of atomics prevents any instruction-level parallelism even if there are no dependencies between the issued operations. Finally, we discuss solutions to the discovered performance issues in the analyzed architectures. Our analysis enables simpler and more effective parallel programming and accelerates data processing on various architectures deployed in both off-the-shelf machines and large compute systems.


title: (How) Do people change their passwords after a breach?

abstract: To protect against misuse of passwords compromised in a breach, consumers should promptly change affected passwords and any similar passwords on other accounts. Ideally, affected companies should strongly encourage this behavior and have mechanisms in place to mitigate harm. In order to make recommendations to companies about how to help their users perform these and other security-enhancing actions after breaches, we must first have some understanding of the current effectiveness of companies' post-breach practices. To study the effectiveness of password-related breach notifications and practices enforced after a breach, we examine---based on real-world password data from 249 participants---whether and how constructively participants changed their passwords after a breach announcement. Of the 249 participants, 63 had accounts on breached domains; only 33% of the 63 changed their passwords and only 13% (of 63) did so within three months of the announcement. New passwords were on average 1.3x stronger than old passwords (when comparing log10-transformed strength), though most were weaker or of equal strength. Concerningly, new passwords were overall more similar to participants' other passwords, and participants rarely changed passwords on other sites even when these were the same or similar to their password on the breached domain. Our results highlight the need for more rigorous password-changing requirements following a breach and more effective breach notifications that deliver comprehensive advice.


title: High-Performance Distributed RMA Locks

abstract: We propose a topology-aware distributed Reader-Writer lock that accelerates irregular workloads for supercomputers and data centers. The core idea behind the lock is a modular design that is an interplay of three distributed data structures: a counter of readers/writers in the critical section, a set of queues for ordering writers waiting for the lock, and a tree that binds all the queues and synchronizes writers with readers. Each structure is associated with a parameter for favoring either readers or writers, enabling adjustable performance that can be viewed as a point in a three dimensional parameter space. We also develop a distributed topology-aware MCS lock that is a building block of the above design and improves state-of-the-art MPI implementations. Both schemes use non-blocking Remote Memory Access (RMA) techniques for highest performance and scalability. We evaluate our schemes on a Cray XC30 and illustrate that they outperform state-of-the-art MPI-3 RMA locking protocols by 81% and 73%, respectively. Finally, we use them to accelerate a distributed hashtable that represents irregular workloads such as key-value stores or graph processing.


title: Combinatorics of criniferous entire maps with escaping critical values

abstract: A transcendental entire function is called criniferous if every point in its escaping set can eventually be connected to infinity by a curve of escaping points. Many transcendental entire functions with bounded singular set have this property, and this class has recently attracted much attention in complex dynamics. In the presence of escaping critical values, these curves break or split at (preimages of) critical points. In this paper, we develop combinatorial tools that allow us to provide a complete description of the escaping set of any criniferous function without asymptotic values on its Julia set. In particular, our description precisely reflects the splitting phenomenon. This combinatorial structure provides the foundation for further study of this class of functions. For example, we use these results in [arXiv:1905.03778] to give the first full description of the topological dynamics of a class of transcendental entire maps with unbounded postsingular set.


title: Anomaly Detection on X-Rays Using Self-Supervised Aggregation Learning

abstract: Deep anomaly detection models using a supervised mode of learning usually work under a closed set assumption and suffer from overfitting to previously seen rare anomalies at training, which hinders their applicability in a real scenario. In addition, obtaining annotations for X-rays is very time consuming and requires extensive training of radiologists. Hence, training anomaly detection in a fully unsupervised or self-supervised fashion would be advantageous, allowing a significant reduction of time spent on the report by radiologists. In this paper, we present SALAD, an end-to-end deep self-supervised methodology for anomaly detection on X-Ray images. The proposed method is based on an optimization strategy in which a deep neural network is encouraged to represent prototypical local patterns of the normal data in the embedding space. During training, we record the prototypical patterns of normal training samples via a memory bank. Our anomaly score is then derived by measuring similarity to a weighted combination of normal prototypical patterns within a memory bank without using any anomalous patterns. We present extensive experiments on the challenging NIH Chest X-rays and MURA dataset, which indicate that our algorithm improves state-of-the-art methods by a wide margin.


title: The STDyn-SLAM: A stereo vision and semantic segmentation approach for SLAM in dynamic outdoor environments

abstract: Commonly, SLAM algorithms are focused on a static environment, however, there are several scenes where dynamic objects are present. This work presents the STDyn-SLAM an image feature-based SLAM system working on dynamic environments using a series of sub-systems, like optic flow, orb features extraction, visual odometry, and convolutional neural networks to discern moving objects in the scene. The neural network is used to support object detection and segmentation to avoid erroneous maps and wrong system localization. The STDyn-SLAM employs a stereo pair and is developed for outdoor environments. Moreover, the processing time of the proposed system is fast enough to run in real-time as it was demonstrated through the experiments given in real dynamic outdoor environments. Further, we compare our SLAM with state-of-the-art methods achieving promising results.


title: Vehicular Visible Light Positioning for Collision Avoidance and Platooning: A Survey

abstract: Relative vehicle positioning methods can contribute to safer and more efficient autonomous driving by enabling collision avoidance and platooning applications. For full automation, these applications require cm-level positioning accuracy and greater than 50 Hz update rate. Since sensor-based methods (e.g., LIDAR, cameras) have not been able to reliably satisfy these requirements under all conditions so far, complementary methods are sought. Recently, positioning based on visible light communication signals from vehicle head/tail LED lights (VLP) has shown significant promise as a complementary method attaining cm-level accuracy and near-kHz rate in realistic driving scenarios. Vehicular VLP methods measure relative bearing (angle) or range (distance) of transmitters (i.e., head/tail lights) based on received signals from on-board photodiodes and estimate transmitter relative positions based on those measurements. In this survey, we first review existing vehicular VLP methods and propose a new method that advances the state-of-the-art in positioning performance. Next, we analyze the theoretical and simulated performance of all methods in realistic driving scenarios under challenging noise and weather conditions, real asymmetric light beam patterns and different vehicle dimensions and light placements. Our simulation results show that the newly proposed VLP method is the overall best performer, and can indeed satisfy the accuracy and rate requirements for localization in collision avoidance and platooning applications within practical constraints. Finally, we discuss remaining open challenges that are faced for the deployment of VLP solutions in the automotive sector and further research questions.


title: A Unified Approach for Autonomous Volumetric Exploration of Large Scale Environments under Severe Odometry Drift

abstract: Exploration is a fundamental problem in robot autonomy. A major limitation, however, is that during exploration robots oftentimes have to rely on on-board systems alone for state estimation, accumulating significant drift over time in large environments. Drift can be detrimental to robot safety and exploration performance. In this work, a submap-based, multi-layer approach for both mapping and planning is proposed to enable safe and efficient volumetric exploration of large scale environments despite odometry drift. The central idea of our approach combines local (temporally and spatially) and global mapping to guarantee safety and efficiency. Similarly, our planning approach leverages the presented map to compute global volumetric frontiers in a changing global map and utilizes the nature of exploration dealing with partial information for efficient local and global planning. The presented system is thoroughly evaluated and shown to outperform state of the art methods even under drift-free conditions. Our system, termed GLoca}, will be made available open source.


title: The Polylogarithm Function in Julia

abstract: The polylogarithm function is one of the constellation of important mathematical functions. It has a long history, and many connections to other special functions and series, and many applications, for instance in statistical physics. However, the practical aspects of its numerical evaluation have not received the type of comprehensive treatments lavished on its siblings. Only a handful of formal publications consider the evaluation of the function, and most focus on a specific domain and/or presume arbitrary precision arithmetic will be used. And very little of the literature contains any formal validation of numerical performance. In this paper we present an algorithm for calculating polylogarithms for both complex parameter and argument and evaluate it thoroughly in comparison to the arbitrary precision implementation in Mathematica. The implementation was created in a new scientific computing language Julia, which is ideal for the purpose, but also allows us to write the code in a simple, natural manner so as to make it easy to port the implementation to other such languages.


title: On the compressive nature of turbulence driven by ionising feedback in the pillars of the Carina Nebula

abstract: The ionizing radiation of massive stars sculpts the surrounding neutral gas into pillar-like structures. Direct signatures of star formation through outflows and jets are observed in these structures, typically at their tips. Recent numerical simulations have suggested that this star formation could potentially be triggered by photoionising radiation, driving compressive modes of turbulence in the pillars. In this study we use recent high-resolution ALMA observations of $^{12}\mathrm{CO}$, $^{13}\mathrm{CO}$, and $\mathrm{C}^{18} \mathrm{O}, \; J = 2-1$ emission to test this hypothesis for pillars in the Carina Nebula. We analyse column density and intensity-weighted velocity maps, and subtract any large-scale bulk motions in the plane of the sky to isolate the turbulent motions. We then reconstruct the dominant turbulence driving mode in the pillars, by computing the turbulence driving parameter $b$, characterised by the relation $\sigma_{\rho/\rho_0} = b \mathcal{M}$ between the standard deviation of the density contrast $\sigma_{\rho/\rho_0}$ (with gas density $\rho$ and its average $\rho_0$) and the turbulent Mach number $\mathcal{M}$. We find values of \mbox{$b\sim0.7$--$1.0$} for most of the pillars, suggesting that predominantly compressive modes of turbulence are driven in the pillars by the ionising radiation from nearby massive stars. We find that this range of $b$ values can produce star formation rates in the pillars that are a factor $\sim 3$ greater than with $b \sim 0.5$, a typical average value of $b$ for spiral-arm molecular clouds. Our results provide further evidence for the potential triggering of star formation in pillars through compressive turbulent motions.


title: Higher Derivatives of the Tangent and Inverse Tangent Functions and Chebyshev Polynomials

abstract: The higher derivatives of the tangent and hyperbolic tangent functions are determined. Formulas for the higher derivatives of the inverse tangent and inverse hyperbolic tangent functions as polynomials are stated and proved. Using another formula for the higher derivatives of the inverse tangent function from literature, two known formulas for the Chebyshev polynomials of the first and second kind are proved. From these formulas the higher derivatives of the inverse tangent and inverse hyperbolic tangent functions in terms of the Chebyshev polynomial of the second kind are provided.


title: $k$-Minkowski-deformation of $U(1)$ gauge theory

abstract: We construct a non-commutative kappa-Minkowski deformation of U(1) gauge theory, following a general approach, recently proposed in JHEP 2008 (2020) 041. We obtain an exact (all orders in the non-commutativity parameter) expression for both the deformed gauge transformations and the deformed field strength, which is covariant under these transformations. The corresponding Yang-Mills Lagrangian is gauge covariant and reproduces the Maxwell Lagrangian in the commutative limit. Gauge invariance of the action functional requires a non-trivial integration measure which, in the commutative limit, does not reduce to the trivial one. We discuss the physical meaning of such a nontrivial commutative limit, relating it to a nontrivial space-time curvature of the undeformed theory. Moreover, we propose a rescaled kappa-Minkowski non-commutative structure, which exhibits a standard flat commutative limit.


title: On an equichordal property of a pair of convex bodies

abstract: Let $d\ge 2$ and let $K$ and $L$ be two convex bodies in ${\mathbb R^d}$ such that $L\subset \textrm{int}\,K$ and the boundary of $L$ does not contain a segment. If $K$ and $L$ satisfy the $(d+1)$-equichordal property, i.e., for any line $l$ supporting the boundary of $L$ and the points $\{\zeta_{\pm}\}$ of the intersection of the boundary of $K$ with $l$, $$ \textrm{dist}^{d+1}(L\cap l, \zeta_+)+\textrm{dist}^{d+1}(L\cap l, \zeta_-)=2\sigma^{d+1} $$ holds, where the constant $\sigma$ is independent of $l$, does it follow that $K$ and $L$ are concentric Euclidean balls? We prove that if $K$ and $L$ have $C^2$-smooth boundaries and $L$ is a body of revolution, then $K$ and $L$ are concentric Euclidean balls.


title: Failure Prediction by Confidence Estimation of Uncertainty-Aware Dirichlet Networks

abstract: Reliably assessing model confidence in deep learning and predicting errors likely to be made are key elements in providing safety for model deployment, in particular for applications with dire consequences. In this paper, it is first shown that uncertainty-aware deep Dirichlet neural networks provide an improved separation between the confidence of correct and incorrect predictions in the true class probability (TCP) metric. Second, as the true class is unknown at test time, a new criterion is proposed for learning the true class probability by matching prediction confidence scores while taking imbalance and TCP constraints into account for correct predictions and failures. Experimental results show our method improves upon the maximum class probability (MCP) baseline and predicted TCP for standard networks on several image classification tasks with various network architectures.


title: Compressing Colour Images with Joint Inpainting and Prediction

abstract: Inpainting-based codecs store sparse, quantised pixel data directly and decode by interpolating the discarded image parts. This interpolation can be used simultaneously for efficient coding by predicting pixel data to be stored. Such joint inpainting and prediction approaches yield good results with simple components such as regular grids and Shepard interpolation on grey value images, but they lack a dedicated mode for colour images. Therefore, we evaluate different approaches for inpainting-based colour compression. Inpainting operators are able to reconstruct a large range of colours from a small colour palette of the known pixels. We exploit this with a luma preference mode which uses higher sparsity in YCbCr colour channels than in the brightness channel. Furthermore, we propose the first full vector quantisation mode for an inpainting-based codec that stores only a small codebook of colours. Our experiments reveal that both colour extensions yield significant improvements.


title: Diffusion-induced blowup solutions for the shadow limit model of a singular Gierer-Meinhardt system

abstract: In the current paper, we provide a thorough investigation of the blowing up behaviour induced via diffusion of the solution of the following non local problem \begin{equation*} \left\{\begin{array}{rcl} \partial_t u &=& \Delta u - u + \displaystyle{\frac{u^p}{ \left(\mathop{\,\rlap{-}\!\!\int}\nolimits_\Omega u^r dr \right)^\gamma }}\quad\text{in}\quad \Omega \times (0,T), \\[0.2cm] \frac{ \partial u}{ \partial \nu} & = & 0 \text{ on } \Gamma = \partial \Omega \times (0,T),\\ u(0) & = & u_0, \end{array} \right. \end{equation*} where $\Omega$ is a bounded domain in $\mathbb{R}^N$ with smooth boundary $\partial \Omega;$ such problem is derived as the shadow limit of a singular Gierer-Meinhardt system, cf. \cite{KSN17, NKMI2018}. Under the Turing type condition $$ \frac{r}{p-1} < \frac{N}{2}, \gamma r \ne p-1, $$ we construct a solution which blows up in finite time and only at an interior point $x_0$ of $\Omega,$ i.e. $$ u(x_0, t) \sim (\theta^*)^{-\frac{1}{p-1}} \left[\kappa (T-t)^{-\frac{1}{p-1}} \right], $$ where $$ \theta^* := \lim_{t \to T} \left(\mathop{\,\rlap{-}\!\!\int}\nolimits_\Omega u^r dr \right)^{- \gamma} \text{ and } \kappa = (p-1)^{-\frac{1}{p-1}}. $$ More precisely, we also give a description on the final asymptotic profile at the blowup point $$ u(x,T) \sim ( \theta^* )^{-\frac{1}{p-1}} \left[ \frac{(p-1)^2}{8p} \frac{|x-x_0|^2}{ |\ln|x-x_0||} \right]^{ -\frac{1}{p-1}} \text{ as } x \to 0, $$ and thus we unveil the form of the Turing patterns occurring in that case due to driven-diffusion instability. The applied technique for the construction of the preceding blowing up solution mainly relies on the approach developed in \cite{MZnon97} and \cite{DZM3AS19}.


title: Nitric Oxide and other molecules: Molecular Modelling and Low FrequencyExploration using the Murchison Widefield Array

abstract: We present new molecular modelling for 14NO and 15NO and a deep, blind molecular line survey at low radio frequencies (99-129 MHz). This survey is the third in a series completed with the Murchison Widefield Array (MWA), but in comparison with the previous surveys, uses four times more data (17 hours vs. 4 hours) and is three times better in angular resolution (1' vs. 3'). The new molecular modelling for nitric oxide and its main isotopologue has seven transitions within the MWA frequency band (although we also present the higher frequency transitions). Although we did not detect any new molecular lines at a limit of 0.21 Jy beam^-1, this work is an important step in understanding the data processing challenges for the future Square Kilometre Array (SKA) and places solid limits on what is expected in the future of low-frequency surveys. The modelling can be utilised for future searches of nitric oxide.


title: Inertia groups in the metastable range

abstract: We prove that the inertia groups of all sufficiently-connected, high-dimensional $(2n)$-manifolds are trivial. This is a key step toward a general classification of manifolds in the metastable range. Specifically, for $m \gg 0$ and $k>5/12$, suppose $M$ is a $\lfloor km \rfloor$-connected, smooth, closed, oriented $m$-manifold and $\Sigma$ is an exotic $m$-sphere. We prove that, if $M \sharp \Sigma$ is diffeomorphic to $M$, then $\Sigma$ bounds a parallelizable manifold. Our proof is built on an understanding of the second extended power functor in Pstr\k{a}gowski's category of synthetic spectra.


title: DeepApple: Deep Learning-based Apple Detection using a Suppression Mask R-CNN

abstract: Robotic apple harvesting has received much research attention in the past few years due to growing shortage and rising cost in labor. One key enabling technology towards automated harvesting is accurate and robust apple detection, which poses great challenges as a result of the complex orchard environment that involves varying lighting conditions and foliage/branch occlusions. This letter reports on the development of a novel deep learning-based apple detection framework named DeepApple. Specifically, we first collect a comprehensive apple orchard dataset for 'Gala' and 'Blondee' apples, using a color camera, under different lighting conditions (sunny vs. overcast and front lighting vs. back lighting). We then develop a novel suppression Mask R-CNN for apple detection, in which a suppression branch is added to the standard Mask R-CNN to suppress non-apple features generated by the original network. Comprehensive evaluations are performed, which show that the developed suppression Mask R-CNN network outperforms state-of-the-art models with a higher F1-score of 0.905 and a detection time of 0.25 second per frame on a standard desktop computer.


title: Evidence for Two Distinct Broad-Line Regions from Reverberation Mapping of PG 0026+129

abstract: We report on the results of a new spectroscopic monitoring campaign of the quasar PG 0026+129 at the Calar Alto Observatory 2.2m telescope from July 2017 to February 2020. Significant variations in the fluxes of the continuum and broad-emission lines, including H$\beta$ and He II, were observed in the first and third years, and clear time lags between them are measured. The broad H$\beta$ line profile consists of two Gaussian components: an intermediate-width H$\beta_{\rm IC}$ with a full width at half-maximum (FWHM) of 1964$\pm$18 $\rm km~s^{-1}$ and another very broad H$\beta_{\rm VBC}$ with a FWHM of 7570$\pm$83 $\rm km~s^{-1}$. H$\beta_{\rm IC}$ has long time lags of $\sim$40--60 days in the rest frame, while H$\beta_{\rm VBC}$ shows nearly zero time delay with respect to the optical continuum at 5100 \AA. The velocity-resolved delays show consistent results: lags of $\sim$30--50 days at the core of the broad H$\beta$ line and roughly zero lags at the wings. H$\beta_{\rm IC}$ has a redshift of $\sim$400 $\rm km~s^{-1}$ which seems to be stable for nearly 30 years by comparing with archived spectra, and may originate from an infall. The root mean square (rms) spectrum of H$\beta_{\rm VBC}$ shows a double-peaked profile with brighter blue peak and extended red wing in the first year, which matches the signature of a thin disk. Both the double-peaked profile and the near-zero lag suggest that H$\beta_{\rm VBC}$ comes from a region associated with the part of the accretion disc that emits the optical continuum. Adopting the FWHM (in the rms spectrum) and the time lag measured for the total H$\beta$ line, and a virial factor of 1.5, we obtain a virial mass of $2.89_{-0.69}^{+0.60} \times10^7 M_{\odot}$ for the central black hole in this quasar.


title: Kohn-Sham Density Functional Theory of Abelian Anyons

abstract: We develop a density functional treatment of non-interacting abelian anyons, which is capable, in principle, of dealing with a system of a large number of anyons in an external potential. Comparison with exact results for few particles shows that the model captures the behavior qualitatively and semi-quantitatively, especially in the vicinity of the fermionic statistics. We then study anyons with statistics parameter $1+1/n$, which are thought to condense into a superconducting state. An indication of the superconducting behavior is the mean-field result that, for uniform density systems, the ground state energy increases under the application of an external magnetic field independent of its direction. Our density-functional-theory based analysis does not find that to be the case for finite systems of anyons, which can accommodate a weak external magnetic field through density transfer between the bulk and the boundary rather than through transitions across effective Landau levels, but the "Meissner repulsion" of the external magnetic field is recovered in the thermodynamic limit as the effect of the boundary becomes negligible. We also consider the quantum Hall effect of anyons, and show that its topological properties, such as the charge and statistics of the excitations and the quantized Hall conductance, arise in a self-consistent fashion.


title: Evolution and ion kinetics of a XUV-induced nanoplasma in ammonia clusters

abstract: High-intensity extreme ultraviolet (XUV) pulses from a free-electron laser can be used to create a nanoplasma in clusters. In Ref. [Michiels et al. PCCP, 2020; 22: 7828-7834] we investigated the formation of excited states in an XUV-induced nanoplasma in ammonia clusters. In the present article we expand our previous study with a detailed analysis of the nanoplasma evolution and ion kinetics. We use a time-delayed UV laser as probe to ionize excited states of H and H$_2^+$ in the XUV-induced plasma. Employing covariance mapping techniques, we show that the correlated emission of protons plays an important role in the plasma dynamics. The time-dependent kinetic energy of the ions created by the probe laser is measured, revealing the charge neutralization of the cluster happens on a sub-picosecond timescale. Furthermore, we observe ro-vibrationally excited molecular hydrogen ions H$_2^{+*}$ being ejected from the clusters. We rationalize our data through a qualitative model of a finite-size non-thermal plasma.


title: Combined synchrotron X-ray diffraction and NV diamond magnetic microscopy measurements at high pressure

abstract: We report the possibility to simultaneously perform wide-field nitrogen-vacancy (NV) diamond magnetic microscopy and synchrotron X-ray diffraction (XRD) measurements at high pressure. NV color centers are created on the culet of a diamond anvil which is integrated in a diamond anvil cell for static compression of the sample. The optically detected spin resonance of the NV centers is used to map the stray magnetic field produced by the sample magnetization. Using this combined scheme, the magnetic and structural behaviors can be simultaneously measured. As a proof-of-principle, we record the correlated {\alpha}-Fe to {\epsilon}-Fe structural and magnetic transitions of iron that occur here between 15 and 20 GPa at 300 K.


title: Combining Ensembles and Data Augmentation can Harm your Calibration

abstract: Ensemble methods which average over multiple neural network predictions are a simple approach to improve a model's calibration and robustness. Similarly, data augmentation techniques, which encode prior information in the form of invariant feature transformations, are effective for improving calibration and robustness. In this paper, we show a surprising pathology: combining ensembles and data augmentation can harm model calibration. This leads to a trade-off in practice, whereby improved accuracy by combining the two techniques comes at the expense of calibration. On the other hand, selecting only one of the techniques ensures good uncertainty estimates at the expense of accuracy. We investigate this pathology and identify a compounding under-confidence among methods which marginalize over sets of weights and data augmentation techniques which soften labels. Finally, we propose a simple correction, achieving the best of both worlds with significant accuracy and calibration gains over using only ensembles or data augmentation individually. Applying the correction produces new state-of-the art in uncertainty calibration across CIFAR-10, CIFAR-100, and ImageNet.


title: Cusped spaces and quasi-isometries of relatively hyperbolic groups

abstract: A group $\Gamma$ with a family of subgroups $\mathbb{P}$ is relatively hyperbolic if $\Gamma$ admits a cusp-uniform action on a proper $\delta$--hyperbolic space. We show that any two such spaces for a given group pair are quasi-isometric, provided the spaces have "constant horospherical distortion," a condition satisfied by Groves--Manning's cusped Cayley graph and by all negatively curved symmetric spaces. Consequently the Bowditch boundary admits a canonical quasisymmetric structure, which coincides with the "naturally occurring" quasisymmetric structure of the symmetric space when considering lattices in rank one symmetric spaces. We show that a group $\Gamma$ is a lattice in a negatively curved symmetric space $X$ if and only if a cusped space for $\Gamma$ is quasi-isometric to the symmetric space. We also prove an ideal triangle characterization of the $\delta$--hyperbolic spaces with uniformly perfect boundary due to Meyer and Bourdon--Kleiner. An appendix concerns the equivalence of several definitions of conical limit point found in the literature.


title: A Concentration of Measure Framework to study convex problems and other implicit formulation problems in machine learning

abstract: This paper provides a framework to show the concentration of solutions $Y^*$ to convex minimizing problem where the objective function $\phi(X)(Y)$ depends on some random vector $X$ satisfying concentration of measure hypotheses. More precisely, the convex problem translates into a contractive fixed point equation that ensure the transmission of the concentration from $X$ to $Y^*$. This result is of central interest to characterize many machine learning algorithms which are defined through implicit equations (e.g., logistic regression, lasso, boosting, etc.). Based on our framework, we provide precise estimations for the first moments of the solution $Y^*$, when $X= (x_1,\ldots, x_n)$ is a data matrix of independent columns and $\phi(X)(y)$ writes as a sum $\frac{1}{n}\sum_{i=1}^n h_i(x_i^TY)$. That allows to describe the behavior and performance (e.g., generalization error) of a wide variety of machine learning classifiers.


title: Exoplanets and the Sun

abstract: We review the recent progress in understanding the jet structures on exoplanets as well as on and inside the Sun. The emphasis is on the more robust aspects of observation and numerical modeling that relate directly to jets. For the exoplanets, the primary focus is on hot-Jupiters since many more observations are available for them presently than other types of exoplanets. Because not much is known about the morphology and strength of the jets on exoplanets, there is currently not much agreement. In contrast, the picture is very different for the Sun. In fact, the jet structure of the Sun is arguably one of the best known jet structures of all the planets and stars, due to the fact that Sun's disk is resolved and its interior can be probed with helioseismology. A discussion of several critical issues pertaining to the modeling of jets on exoplanets and the Sun is presented, along with a brief outlook on the subject.


title: Modelling Complex Survey Data Using R, SAS, SPSS and Stata: A Comparison Using CLSA Datasets

abstract: The R software has become popular among researchers due to its flexibility and open-source nature. However, researchers in the fields of public health and epidemiological studies are more customary to commercial statistical softwares such as SAS, SPSS and Stata. This paper provides a comprehensive comparison on analysis of health survey data using the R survey package, SAS, SPSS and Stata. We describe detailed R codes and procedures for other software packages on commonly encountered statistical analyses, such as estimation of population means and regression analysis, using datasets from the Canadian Longitudinal Study on Aging (CLSA). It is hoped that the paper stimulates interest among health science researchers to carry data analysis using R and also serves as a cookbook for statistical analysis using different software packages.


title: New Solutions for Rotating Boson Stars

abstract: It has been shown that scalar fields can form gravitationally bound compact objects called boson stars. In this study, we analyze boson star configurations where the scalar fields contain a small amount of angular momentum and find two new classes of solutions. In the first case all particles are in the same slowly rotating state and in the second case the majority of particles are in the non-rotating ground state and a small number of particles are in an excited rotating state. In both cases, we solve the underlying Gross-Pitaevskii-Poisson equations that describe the profile of these compact objects both numerically as well as analytically through series expansions.


title: Parity of the coefficients of certain eta-quotients

abstract: We investigate the parity of the coefficients of certain eta-quotients, extensively examining the case of $m$-regular partitions. Our theorems concern the density of their odd values, in particular establishing lacunarity modulo 2 for specified coefficients; self-similarities modulo 2; and infinite families of congruences in arithmetic progressions. For all $m \leq 28$, we either establish new results of these types where none were known, extend previous ones, or conjecture that such results are impossible. All of our work is consistent with a new, overarching conjecture that we present for arbitrary eta-quotients, greatly extending Parkin-Shanks' classical conjecture for the partition function. We pose several other open questions throughout the paper, and conclude by suggesting a list of specific research directions for future investigations in this area.


title: Circular Convolution and Product Theorem for Affine Discrete Fractional Fourier Transform

abstract: The Fractional Fourier Transform is a ubiquitous signal processing tool in basic and applied sciences. The Fractional Fourier Transform generalizes every property and application of the Fourier Transform. Despite the practical importance of the discrete fractional Fourier transform, its applications in digital communications have been elusive. The convolution property of the discrete Fourier transform plays a vital role in designing multi-carrier modulation systems. Here we report a closed-form affine discrete fractional Fourier transform and we show the circular convolution property for it. The proposed approach is versatile and generalizes the discrete Fourier transform and can find applications in Fourier based signal processing tools.


title: Impact of invasive metal probes on Hall measurements in semiconductor nanostructures

abstract: Recent advances in bottom-up growth are giving rise to a range of new two-dimensional nanostructures. Hall effect measurements play an important role in their electrical characterization. However, size constraints can lead to device geometries that deviate significantly from the ideal of elongated Hall bars with currentless contacts. Many devices using these new materials have a low aspect ratio and feature metal Hall probes that overlap with the semiconductor channel. This can lead to a significant distortion of the current flow. We present experimental data from InAs 2D nanofin devices with different Hall probe geometries to study the influence of Hall probe length and width. We use finite-element simulations to further understand the implications of these aspects and expand the scope to contact resistance and sample aspect ratios. Our key finding is that invasive probes lead to a significant underestimation in the measured Hall voltage, typically of the order of 40-80%. This in turn leads to a subsequent proportional overestimation of carrier concentration and an underestimation of mobility