title: HCNAF: Hyper-Conditioned Neural Autoregressive Flow and its Application for Probabilistic Occupancy Map Forecasting

abstract: We introduce Hyper-Conditioned Neural Autoregressive Flow (HCNAF); a powerful universal distribution approximator designed to model arbitrarily complex conditional probability density functions. HCNAF consists of a neural-net based conditional autoregressive flow (AF) and a hyper-network that can take large conditions in non-autoregressive fashion and outputs the network parameters of the AF. Like other flow models, HCNAF performs exact likelihood inference. We conduct a number of density estimation tasks on toy experiments and MNIST to demonstrate the effectiveness and attributes of HCNAF, including its generalization capability over unseen conditions and expressivity. Finally, we show that HCNAF scales up to complex high-dimensional prediction problems of the magnitude of self-driving and that HCNAF yields a state-of-the-art performance in a public self-driving dataset.


title: A learning-based algorithm to quickly compute good primal solutions for Stochastic Integer Programs

abstract: We propose a novel approach using supervised learning to obtain near-optimal primal solutions for two-stage stochastic integer programming (2SIP) problems with constraints in the first and second stages. The goal of the algorithm is to predict a "representative scenario" (RS) for the problem such that, deterministically solving the 2SIP with the random realization equal to the RS, gives a near-optimal solution to the original 2SIP. Predicting an RS, instead of directly predicting a solution ensures first-stage feasibility of the solution. If the problem is known to have complete recourse, second-stage feasibility is also guaranteed. For computational testing, we learn to find an RS for a two-stage stochastic facility location problem with integer variables and linear constraints in both stages and consistently provide near-optimal solutions. Our computing times are very competitive with those of general-purpose integer programming solvers to achieve a similar solution quality.


title: Improved Surrogates in Inertial Confinement Fusion with Manifold and Cycle Consistencies

abstract: Neural networks have become very popular in surrogate modeling because of their ability to characterize arbitrary, high dimensional functions in a data driven fashion. This paper advocates for the training of surrogates that are consistent with the physical manifold -- i.e., predictions are always physically meaningful, and are cyclically consistent -- i.e., when the predictions of the surrogate, when passed through an independently trained inverse model give back the original input parameters. We find that these two consistencies lead to surrogates that are superior in terms of predictive performance, more resilient to sampling artifacts, and tend to be more data efficient. Using Inertial Confinement Fusion (ICF) as a test bed problem, we model a 1D semi-analytic numerical simulator and demonstrate the effectiveness of our approach. Code and data are available at https://github.com/rushilanirudh/macc/


title: Monotone Catenary Degree in Numerical Monoids

abstract: Recent investigations on the catenary degrees of numerical monoids have demonstrated that this invariant is a powerful tool in understanding the factorization theory of this class of monoids. Although useful, the catenary degree is largely not sensitive to the lengths of factorizations of an element. In this paper, we study the monotone catenary degree of numerical monoids, which is a variant of catenary degree that requires chains run through factorization lengths monotonically. In general, the monotone catenary is greater than or equal to the catenary degree. We begin by providing an important class of monoids (arithmetical numerical monoids) for which monotone catenary degree is equal to the catenary degree. Conversely, we provide several classes of embedding dimension 3 numerical monoids where monotone catenary degree is strictly greater. We conclude by showing that this difference can grow arbitrarily large.


title: On the configurations of nine points on a cubic curve

abstract: We study the reciprocal position of nine points in the plane, according to their collinearities. In particular, we consider the case in which the nine points are contained in an irreducible cubic curve and we give their classification. If we consider two configurations different when the associated incidence structures are not isomorphic, we see that there are 131 configurations that can be realized in $\mathbb{P}^2_{\mathbb{Q}}$, and there are two more in $\mathbb{P}^2_K$, where $K =\mathbb{Q}[\sqrt{-3}]$ (one of the two is the Hesse configuration given by the nine inflection points of a cubic curve). Finally, we compute the possible Hilbert functions of the ideals of the nine points.


title: When Your Robot Breaks: Active Learning During Plant Failure

abstract: Detecting and adapting to catastrophic failures in robotic systems requires a robot to learn its new dynamics quickly and safely to best accomplish its goals. To address this challenging problem, we propose probabilistically-safe, online learning techniques to infer the altered dynamics of a robot at the moment a failure (e.g., physical damage) occurs. We combine model predictive control and active learning within a chance-constrained optimization framework to safely and efficiently learn the new plant model of the robot. We leverage a neural network for function approximation in learning the latent dynamics of the robot under failure conditions. Our framework generalizes to various damage conditions while being computationally light-weight to advance real-time deployment. We empirically validate within a virtual environment that we can regain control of a severely damaged aircraft in seconds and require only 0.1 seconds to find safe, information-rich trajectories, outperforming state-of-the-art approaches.


title: New K\"ahler invariant Fayet-Iliopoulos terms in supergravity and cosmological applications

abstract: Recently, a new type of constant Fayet-Iliopoulos (FI) terms was introduced in $\mathcal{N}=1$ supergravity, which do not require the gauging of the $R$-symmetry. We revisit and generalise these constructions, building a new class of K\"ahler invariant FI terms parametrised by a function of the gravitino mass as functional of the chiral superfields, which is then used to describe new models of inflation. They are based on a no-scale supergravity model of the inflaton chiral multiplet, supplemented by an abelian vector multiplet with the new FI-term. We show that the inflaton potential is compatible with the CMB observational data, with a vacuum energy at the minimum that can be tuned to a tiny positive value. Finally, the axionic shift symmetry can be gauged by the $U(1)$ which becomes massive. These models offer a mechanism for fixing the gravitino mass in no-scale supergravities, that corresponds to a flat direction of the scalar potential in the absence of the new FI-term; its origin in string theory is an interesting open problem.


title: Full-spectrum high resolution modeling of the dielectric function of water

abstract: In view of the vital role of water in chemical and physical processes, an exact knowledge of its dielectric function over a large frequency range is important. In this article we report on currently available measurements of the dielectric function of water at room temperature (25$^{\circ}$C) across the full electromagnetic spectrum: microwave, IR, UV and X-ray (up to 100 eV). We provide parameterisations of the complex dielectric function of water with two Debye (microwave) oscillators and high resolution of IR and UV/X-ray oscillators. We also report dielectric parameters for ice-cold water with a microwave/IR spectrum measured at $0.4^\circ$C, while taking the UV spectrum from 25$^{\circ}$C (assuming negligible temperature dependence in UV). We illustrate the consequences of the model via calculations of van der Waals interactions of gas molecules near water surfaces, and an assessment of the thickness of water films on ice and ice films on water. In contrast to earlier models of ice-cold water, we predict that a micron-scale layer of ice is stabilised on a bulk water surface. Similarly, the van der Waals interaction promotes complete freezing rather than supporting a thin premelting layer of water on a bulk ice surface. Density-based extrapolation from warm to cold water of the dielectric function at imaginary frequencies is found to be satisfactory in the microwave but poor (40% error) at IR frequencies.


title: NOMA versus OMA in Finite Blocklength Regime: Link-Layer Rate Performance

abstract: In this paper, we investigate the latency performance of non-orthogonal multiple access (NOMA) and orthogonal multiple access (OMA) technologies in finite blocklength regime. In the comparative study, we derive the achievable effective capacity of two-user NOMA and its OMA counterpart under delay quality-of-service constraints. We then obtain closed-form expressions for the achievable effective capacity of the weak and strong users in both scenarios considering transmissions over Rayleigh fading channels. Numerical results are provided. In particular, it is shown that at low signal-to-noise ratios (SNRs), the OMA user with better channel condition outperforms both NOMA users. \blue{We also evaluate the impact of fixed power allocation scheme on the achievable effective capacity of two-user NOMA.} The comparative analysis of the total link-layer rate shows that at high SNRs, the total link-layer rate of NOMA with finite blocklength outperforms the one of OMA when the delay exponent is loose.


title: Thermodynamic and ordering kinetics in asymmetric PS-b-PMMA block copolymer thin films

abstract: The ordering kinetics of standing cylinder-forming polystyrene-block-poly (methyl methacrylate) block copolymers (molecular weight: 39 kg/mol) close to the order-disorder transition is experimentally investigated following the temporal evolution of the correlation length at different annealing temperatures. The grow exponent of the grain coarsening process is determined to be 1/2, signature of a curvature driven ordering mechanism. The measured activation enthalpy and the resulting Meyer-Neldel temperature for this specific copolymer along with the data already known for PS-b-PMMA block copolymers in strong segregation limit, allows investigating the interplay between the ordering kinetics and the thermodynamic driving force during the grain coarsening. These findings unveil various phenomena concomitantly occurring during the thermally activated ordering kinetics at segmental, single chain, and collective level.


title: Modeling Escape from a One-Dimensional Potential Well at Zero or Very Low Temperatures

abstract: The process of activation out a one-dimensional potential is investigated systematically in zero and nonzero temperature conditions. The features of the potential are traced through statistical escape out of its wells whose depths are tuned in time by a forcing term. The process is carried out on the damped pendulum system imposing specific initial conditions on the potential variable. While for relatively high values of the dissipation the statistical properties follow a behavior that can be derived from the standard Kramers model, decreasing the dissipation we observe responses/deviations which have regular dependencies on initial conditions, temperature, and loss parameter itself. It is shown that failures of the thermal activation model are originated at low temperatures, and very low dissipation, by the initial conditions and intrinsic, namely T=0, characteristic oscillations of the potential-generated dynamical equation.


title: Estimating the maximum gravitational mass of nonrotating neutron stars from the GW170817/GRB 170817A/AT2017gfo observation

abstract: Assuming that the differential rotation of the massive neutron star (NS) formed in the double NS (DNS) mergers has been effectively terminated by the magnetic braking and a uniform rotation has been subsequently established (i.e., a supramassive NS is formed), we analytically derive in this work an approximated expression for the critical total gravitational mass ($M_{\rm tot,c}$) to form supramassive NS (SMNS) in the DNS mergers, benefited from some equation of state (EoS) insensitive relationships. The maximum gravitational mass of the nonrotating NSs ($M_{\rm TOV}$) as well as the dimensionless angular momentum of the remnant ($j$) play the dominant roles in modifying $M_{\rm tot,c}$, while the radius and mass differences of the premerger NSs do not. The GW170817/GRB 170817A/AT2017gfo observations have provided so far the best opportunity to quantitatively evaluate $M_{\rm TOV}$. Supposing the central engine for GRB 170817A is a black hole quickly formed in the collapse of an SMNS, we find $M_{\rm TOV}=2.13^{+0.09}_{-0.08}M_\odot$ (68.3% credibility interval, including also the uncertainties of the EoS insensitive relationships), which is consistent with the constraints set by current NS mass measurements.


title: Microscopic theory of the friction force exerted on a quantum impurity in one-dimensional quantum liquids

abstract: We study the motion of a slow quantum impurity in one-dimensional environments focusing on systems of strongly interacting bosons and weakly interacting fermions. While at zero temperature the impurity motion is frictionless, at low temperatures finite friction appears. The dominant process is the scattering of the impurity off two fermionic quasiparticles. We evaluate the friction force and show that, at low temperatures, it scales either as the fourth or the sixth power of temperature, depending on the system parameters. This is a result of the scattering of the impurity off two fermionic quasiparticles that are situated around different Fermi points. It is the dominant process at low temperatures. We also evaluate the contribution to the friction force originating from the scattering of the impurity off two fermionic quasiparticles that are situated around different Fermi points. It behaves as the tenth power of temperature.


title: SpaRCe: Improved Learning of Reservoir Computing Systems through Sparse Representations

abstract: "Sparse" neural networks, in which relatively few neurons or connections are active, are common in both machine learning and neuroscience. Whereas in machine learning, "sparsity" is related to a penalty term that leads to some connecting weights becoming small or zero, in biological brains, sparsity is often created when high spiking thresholds prevent neuronal activity. Here we introduce sparsity into a reservoir computing network via neuron-specific learnable thresholds of activity, allowing neurons with low thresholds to contribute to decision-making but suppressing information from neurons with high thresholds. This approach, which we term "SpaRCe", optimises the sparsity level of the reservoir without affecting the reservoir dynamics. The read-out weights and the thresholds are learned by an on-line gradient rule that minimises an error function on the outputs of the network. Threshold learning occurs by the balance of two opposing forces: reducing inter-neuronal correlations in the reservoir by deactivating redundant neurons, while increasing the activity of neurons participating in correct decisions. We test SpaRCe on classification problems and find that threshold learning improves performance compared to standard reservoir computing. SpaRCe alleviates the problem of catastrophic forgetting, a problem most evident in standard echo state networks and recurrent neural networks in general, due to increasing the number of task-specialised neurons that are included in the network decisions.


title: Quartic monoid surfaces with maximum number of lines

abstract: In 1884 the German mathematician Karl Rohn published a substantial paper on \cite{ROH} on the properties of quartic surfaces with triple points, proving (among many other things) that the maximum number of lines contained in a quartic monoid surface is $31$. In this paper we study in details this class of surfaces. We prove that there exists an open subset $A \subseteq \mathbb{P}^1_K$ ($K$ is a characteristic zero field) that parametrizes (up to a projectivity) all the quartic monoid surfaces with $31$ lines; then we study the action of $\mathrm{PGL}(4,K)$ on these surfaces, we show that the stabiliser of each of them is a group isomorphic to $S_3$ except for one surface of the family, whose stabiliser is a group isomorphic to $S_3 \times C_3$. Finally we show that the $j$-invariant allows one to decide, also in this situation, when two elements of $A$ give the same surface up to a projectivity. To get our results, several computational tools, available in computer algebra systems, are used.


title: Theta lifts for Lorentzian lattices and coefficients of mock theta functions

abstract: We evaluate regularized theta lifts for Lorentzian lattices in three different ways. In particular, we obtain formulas for their values at special points involving coefficients of mock theta functions. By comparing the different evaluations, we derive recurrences for the coefficients of mock theta functions, such as Hurwitz class numbers, Andrews' spt-function, and Ramanujan's mock theta functions.


title: A weighted central limit theorem for $\log|\zeta(1/2+it)|$

abstract: Under the Riemann Hypothesis, we show that as $t$ varies in $T\leq t \leq 2T$, the distribution of $\log|\zeta(1/2+it)|$ with respect to the measure $|\zeta(1/2+it)|^2dt$ is approximately normal with mean $\log\log T$ and variance $\frac{1}{2}\log\log T$.


title: Form class groups and class fields of CM-fields

abstract: Let $F$ be a totally real number field of class number one, and let $K$ be a CM-field with $F$ as its maximal real subfield. For each positive integer $N$, we construct a class group of certain binary quadratic forms over $F$ which is isomorphic to the ray class group of $K$ modulo $N$. Assuming further that the narrow class number of $F$ is one, we construct a class field of the reflex field of $K$ in terms of the singular values of Hilbert modular functions.


title: Deep XMM-Newton Observations Reveal the Origin of Recombining Plasma in the Supernova Remnant W44

abstract: Recent X-ray studies revealed over-ionized recombining plasmas (RPs) in a dozen mixed-morphology (MM) supernova remnants (SNRs). However, the physical process of the over-ionization has not been fully understood yet. Here we report on spatially resolved spectroscopy of X-ray emission from W44, one of the over-ionized MM-SNRs, using XMM-Newton data from deep observations, aiming to clarify the physical origin of the over-ionization. We find that combination of low electron temperature and low recombination timescale is achieved in the region interacting with dense molecular clouds. Moreover, a clear anti-correlation between the electron temperature and the recombining timescale is obtained from each of the regions with and without the molecular clouds. The results are well explained if the plasma was over-ionized by rapid cooling through thermal conduction with the dense clouds hit by the blast wave of W44. Given that a few other over-ionized SNRs show evidence for adiabatic expansion as the major driver of the rapid cooling, our new result indicates that both processes can contribute to over-ionization in SNRs, with the dominant channel depending on the evolutionary stage.


title: General multilevel adaptations for stochastic approximation algorithms

abstract: In this article we establish central limit theorems for multilevel Polyak-Ruppert averaged stochastic approximation schemes. We work under very mild technical assumptions and consider the slow regime in wich typical errors decay like $N^{-\delta}$ with $\delta\in(0,\frac 12)$ and the critical regime in which errors decay of order $N^{-1/2}\sqrt{\log N}$ in the runtime $N$ of the algorithm.


title: Universal algebraic growth of entanglement entropy in many-body localized systems with power-law interactions

abstract: Power-law interactions play a key role in a large variety of physical systems. In the presence of disorder, these systems may undergo many-body localization for a sufficiently large disorder. Within the many-body localized phase the system presents in time an algebraic growth of entanglement entropy, $S_{vN}(t)\propto t^{\gamma}$. Whereas the critical disorder for many-body localization depends on the system parameters, we find by extensive numerical calculations that the exponent $\gamma$ acquires a universal value $\gamma_c\simeq 0.33$ at the many-body localization transition, for different lattice models and decay powers. Moreover, our results suggest an intriguing relation between $\gamma_c$ and the critical minimal decay power of interactions necessary for the observation of many-body localization.


title: Bernoulli property for certain skew products over hyperbolic systems

abstract: We study the Bernoulli property for a class of partially hyperbolic systems arising from skew products. More precisely, we consider a hyperbolic map $(T,M,\mu)$, where $\mu$ is a Gibbs measure, an aperiodic H\"older continuous cocycle $\phi:M\to \mathbb R$ with zero mean and a zero-entropy flow $(K_t,N,\nu)$. We then study the skew product $$ T_\phi(x,y)=(Tx,K_{\phi(x)}y), $$ acting on $(M\times N,\mu \times \nu)$. We show that if $(K_t)$ is of slow growth and has good equidistribution properties, then $T_\phi$ remains Bernoulli. In particular, our main result applies to $(K_t)$ being a typical translation flow on a surface of genus $\geq 1$ or a smooth reparametrization of isometric flows on $\mathbb T^2$. This provides examples of non-algebraic, partially hyperbolic systems which are Bernoulli and for which the center is non-isometric (in fact might be weakly mixing).


title: Modeling tomographic measurements of photoelectron vortices in counter-rotating circularly polarized laser pulses

abstract: Recent experiments [D. Pengel, S. Kerbstadt, L. Englert, T. Bayer, and M. Wollenhaupt, \href{https://journals.aps.org/pra/abstract/10.1103/PhysRevA.96.043426}{{\PRA} {\bf 96} 043426 (2017)}] have measured the photoelectron momentum distribution for three-photon ionization of potassium by counter-rotating circularly polarized 790-nm laser pulses. The distribution displays spiral vortices, arising from the interference of ionizing wavepackets with different magnetic quantum numbers. The high level of multidimensional detail observed in the distribution makes this an ideal case in which to demonstrate the accuracy of emerging theoretical techniques applicable to such problems. We use the \(R\)-matrix with time dependence approach to investigate this process. We calculate the full-dimensional photoelectron momentum distribution, and compare against a set of planar projections of this distribution previously measured in experiment.


title: Particle Ratios in a Multi Component Non-Ideal Hadron Resonance Gas

abstract: We have considered formation of a multi component non-ideal hot and dense gas of hadronic resonances in the ultra-relativistic heavy-ion collisions. In the statistical thermal model approach the equation of state (EoS) of the non interacting ideal hadron resonance gas (IHRG) does not incorporate either the attractive part or the short range repulsive part of the baryonic interaction. On the other hand in the non-ideal hadron resonance gas (NIHRG) model we can incorporate these interactions using the Van der Waals (VDW) type approach. Studies have been made to see its effect on the critical parameters of the quark-hadron phase transition. However, it can also lead to modifications in the calculated relative particle yields. In this paper we have attempted to understand the effect of such Van der Waals type interactions on the relative particle yields and also studied their dependences on the system thermal parameters, such as the temperature and baryon chemical potential ({\mu}_B). We have also taken into account the decay contributions of the heavier resonances. These results on particle ratios are compared with the corresponding results obtained from the point-like i.e. non interacting IHRG model. It is found that the particle ratios get modified by incorporating the Van der Waals type interactions, especially in a baryon rich system which is expected to be formed at lower RHIC energies, SPS energies and in the forthcoming CBM experiments due to high degree of nuclear stopping in these experiments.


title: On Exact Proper Time Wilsonian RG Flows

abstract: We discuss the possibility to define exact RG equations for a UV regulated Wilsonian action based on a proper time (PT) regulator function. We start from a functional mapping which shows how each particular flow equation (and RG scheme) is associated to infinitely many scale dependent field redefinitions, which are related to specific coarse-graining procedures. On specializing to a sub-family of one parameter PT regulators we briefly analyze few results for the Ising Universality class in three dimensions, obtained within a second order truncation in the derivative expansion of the Wilsonian action.


title: Direction Concentration Learning: Enhancing Congruency in Machine Learning

abstract: One of the well-known challenges in computer vision tasks is the visual diversity of images, which could result in an agreement or disagreement between the learned knowledge and the visual content exhibited by the current observation. In this work, we first define such an agreement in a concepts learning process as congruency. Formally, given a particular task and sufficiently large dataset, the congruency issue occurs in the learning process whereby the task-specific semantics in the training data are highly varying. We propose a Direction Concentration Learning (DCL) method to improve congruency in the learning process, where enhancing congruency influences the convergence path to be less circuitous. The experimental results show that the proposed DCL method generalizes to state-of-the-art models and optimizers, as well as improves the performances of saliency prediction task, continual learning task, and classification task. Moreover, it helps mitigate the catastrophic forgetting problem in the continual learning task. The code is publicly available at https://github.com/luoyan407/congruency.


title: The relative trace formula in analytic number theory

abstract: We discuss a variety of applications of the relative trace formula of Kuznetsov type in analytic number theory and the theory of automorphic forms.


title: Detection of a Source Code Plagiarism in a Student Programming Competition

abstract: The article presents a system for testing the independence of solutions to algorithmic problems sent by students as part of the student programming competition. First, the context was discussed, as well as the need to organize programming competitions resulting from this context. Then, an algorithm was proposed to study the mutual similarity of source codes of programs sent as part of a programming competition. Since, after implementation, the algorithm was used in practice, examples of its application for detecting the plagiarism of source codes of solutions in two programming competitions conducted as part of classes on Algorithms and Numerical Methods were also presented. Finally, the effectiveness of the solutions used in the work was discussed.


title: Test of lepton universality with $\Lambda^{0}_{b} \to p K^- \ell^+ \ell^-$ decays

abstract: The ratio of branching fractions of the decays $\Lambda^{0}_{b}\to pK^{-}e^{+}e^{-}$ and $\Lambda^{0}_{b}\to pK^{-}\mu^{+}\mu^{-}$, $R^{-1}_{pK}$, is measured for the first time using proton-proton collision data corresponding to an integrated luminosity of 4.7 $fb^{-1}$ recorded with the LHCb experiment at center-of-mass energies of 7, 8 and 13 TeV. In the dilepton mass-squared range $0.1 < q^{2} < 6.0$ $GeV^{2}/c^{4}$ and the $pK^{-}$ mass range $m(pK^{-}) < 2600$ $MeV/c^{2}$, the ratio of branching fractions is measured to be $R^{-1}_{pK} = 1.17 ^{+0.18}_{-0.16} \pm 0.07$, where the first uncertainty is statistical and the second systematic. This is the first test of lepton universality with b baryons and the first observation of the decay $\Lambda^{0}_{b}\to pK^{-}e^{+}e^{-}$.


title: On-the-fly Global Embeddings Using Random Projections for Extreme Multi-label Classification

abstract: The goal of eXtreme Multi-label Learning (XML) is to automatically annotate a given data point with the most relevant subset of labels from an extremely large vocabulary of labels (e.g., a million labels). Lately, many attempts have been made to address this problem that achieve reasonable performance on benchmark datasets. In this paper, rather than coming-up with an altogether new method, our objective is to present and validate a simple baseline for this task. Precisely, we investigate an on-the-fly global and structure preserving feature embedding technique using random projections whose learning phase is independent of training samples and label vocabulary. Further, we show how an ensemble of multiple such learners can be used to achieve further boost in prediction accuracy with only linear increase in training and prediction time. Experiments on three public XML benchmarks show that the proposed approach obtains competitive accuracy compared with many existing methods. Additionally, it also provides around 6572x speed-up ratio in terms of training time and around 14.7x reduction in model-size compared to the closest competitors on the largest publicly available dataset.


title: Performance of regression models as a function of experiment noise

abstract: A challenge in developing machine learning regression models is that it is difficult to know whether maximal performance has been reached on a particular dataset, or whether further model improvement is possible. In biology this problem is particularly pronounced as sample labels (response variables) are typically obtained through experiments and therefore have experiment noise associated with them. Such label noise puts a fundamental limit to the performance attainable by regression models. We address this challenge by deriving a theoretical upper bound for the coefficient of determination (R2) for regression models. This theoretical upper bound depends only on the noise associated with the response variable in a dataset as well as its variance. The upper bound estimate was validated via Monte Carlo simulations and then used as a tool to bootstrap performance of regression models trained on biological datasets, including protein sequence data, transcriptomic data, and genomic data. Although we study biological datasets in this work, the new upper bound estimates will hold true for regression models from any research field or application area where response variables have associated noise.


title: Causality matters in medical imaging

abstract: This article discusses how the language of causality can shed new light on the major challenges in machine learning for medical imaging: 1) data scarcity, which is the limited availability of high-quality annotations, and 2) data mismatch, whereby a trained algorithm may fail to generalize in clinical practice. Looking at these challenges through the lens of causality allows decisions about data collection, annotation procedures, and learning strategies to be made (and scrutinized) more transparently. We discuss how causal relationships between images and annotations can not only have profound effects on the performance of predictive models, but may even dictate which learning strategies should be considered in the first place. For example, we conclude that semi-supervision may be unsuitable for image segmentation---one of the possibly surprising insights from our causal analysis, which is illustrated with representative real-world examples of computer-aided diagnosis (skin lesion classification in dermatology) and radiotherapy (automated contouring of tumours). We highlight that being aware of and accounting for the causal relationships in medical imaging data is important for the safe development of machine learning and essential for regulation and responsible reporting. To facilitate this we provide step-by-step recommendations for future studies.


title: An Extreme Metastable Model for Quantum Measurement

abstract: Quantum experiments are observed as probability density functions. We often encounter multi-peaked densities which we model in this paper by a metastable dynamical system. The dynamics can be regarded in a thought experiment where a mouse is in either one of two disjoint traps which possess little outlets. The mouse spends most of its time in one or the other trap, and once in awhile makes its way to the other trap. Hence, a bi-peaked density function. When the experiment is observed - such as a light shone on the traps - the mouse stays in one trap. This measurement process results in a single peaked density that, we prove, can be modelled by the metastable process selecting an extreme density function, which is single peaked.


title: Constructing a control-ready model of EEG signal during general anesthesia in humans

abstract: Significant effort toward the automation of general anesthesia has been made in the past decade. One open challenge is in the development of control-ready patient models for closed-loop anesthesia delivery. Standard depth-of-anesthesia tracking does not readily capture inter-individual differences in response to anesthetics, especially those due to age, and does not aim to predict a relationship between a control input (infused anesthetic dose) and system state (commonly, a function of electroencephalography (EEG) signal). In this work, we developed a control-ready patient model for closed-loop propofol-induced anesthesia using data recorded during a clinical study of EEG during general anesthesia in ten healthy volunteers. We used principal component analysis to identify the low-dimensional state-space in which EEG signal evolves during anesthesia delivery. We parameterized the response of the EEG signal to changes in propofol target-site concentration using logistic models. We note that inter-individual differences in anesthetic sensitivity may be captured by varying a constant cofactor of the predicted effect-site concentration. We linked the EEG dose-response with the control input using a pharmacokinetic model. Finally, we present a simple nonlinear model predictive control in silico demonstration of how such a closed-loop system would work.


title: A noncommutative calculus on the cyclic dual of Ext

abstract: We show that if the cochain complex computing Ext groups (in the category of modules over Hopf algebroids) admits a cocyclic structure, then the noncommutative Cartan calculus structure on Tor over Ext dualises in a cyclic sense to a calculus on Coext over Cotor. More precisely, the cyclic duals of the chain resp. cochain spaces computing the two classical derived functors lead to complexes that compute the more exotic ones, giving a cyclic opposite module over an operad with multiplication that induce operations such as a Lie derivative, a cap product (or contraction), and a (cyclic) differential, along with higher homotopy operators defining a noncommutative Cartan calculus up to homotopy. In particular, this allows to recover the classical Cartan calculus from differential geometry or the Chevalley-Eilenberg calculus for Lie(-Rinehart) algebras without any finiteness conditions or the use of topological tensor products.


title: On nilpotent automorphism groups of function fields

abstract: We study the automorphisms of a function field of genus $g\geq 2$ over an algebraically closed field of characteristic $p>0$. More precisely, we show that the order of a nilpotent subgroup $G$ of its automorphism group is bounded by $16 (g-1)$ when G is not a $p$-group. We show that if $|G|=16(g-1) $, then $g-1$ is a power of $2$. Furthermore, we provide an infinite family of function fields attaining the bound.


title: New Bounds on Antipowers in Words

abstract: Fici et al. defined a word to be a k-power if it is the concatenation of k consecutive identical blocks, and an r-antipower if it is the concatenation of r pairwise distinct blocks of the same size. They defined N (k, r) as the smallest l such that every binary word of length l contains either a k-power or an r-antipower. In this note we obtain some new upper and lower bounds on N (k, r). We also consider avoiding 3-antipowers and 4-antipowers over larger alphabets, and obtain a lower bound for N (k, 5) in the binary case.


title: Enhanced LMMSE Estimation Capable of Selecting Parameters

abstract: In the linear minimum mean square error (LMMSE) estimation for orthogonal frequency division multiplexing (OFDM) systems, the problem about the determination of the algorithm's parameters, especially those related with channel frequency response (CFR) correlation, has not been readily solved yet. Although many approaches have been proposed to determine the statistic parameters, it is hard to choose the best one within those approaches in the design phase, since every approach has its own most suitable application conditions and the real channel condition is unpredictable. In this paper, we propose an enhance LMMSE estimation capable of selecting parameters by itself. To this end, sampled noise MSE is first proposed to evaluate the practical performance of interpolation. Based on this evaluation index, a novel parameter comparison scheme is proposed to determine the parameters which can endow LMMSE estimation best performance within a parameter set. After that, the structure of the enhanced LMMSE is illustrated, and it is applied in OFDM systems. Besides, the issues about theoretical analysis on accuracy of the parameter comparison scheme, the parameter set design and algorithm complexity are explained in detail. At last, our analyses and performance of the proposed estimation method are demonstrated by simulation experiments.


title: Estimation of Residential Radon Concentration in Pennsylvania Counties by Data Fusion

abstract: A data fusion method for the estimation of residential radon level distribution in any Pennsylvania county is proposed. The method is based on a multi-sample density ratio model with variable tilts and is applied to combined radon data from a reference county of interest and its neighboring counties. Beaver county and its four immediate neighbors are taken as a case in point. The distribution of radon concentration is estimated in each of six periods, and then the analysis is repeated combining the data from all the periods to obtain estimates of Beaver threshold probabilities and the corresponding confidence intervals.


title: Nonlinear scattering of Fast Radio Bursts

abstract: Nonlinear radio waves modulate the plasma, scatter on the modulations, and develop an intermittent power spectrum -- perhaps. Rudiments of theory, numerical simulations, and qualitative modeling of nonlinear scattering are presented. Nonlinear radio waves in electron-positron plasma are shown to be unstable, the instability growth rate is calculated. The long-term outcome of the instability is studied numerically; an intermittent power spectrum is found. The frequency coherence scale of the intermittent power spectrum is interpreted as a measure of the nonlinear scattering screen.


title: Stellar populations and physical properties of starbursts in the Antennae galaxy from self-consistent modelling of MUSE spectra

abstract: We have modelled the stellar and nebular continua and emission-line intensity ratios of massive stellar populations in the Antennae galaxy using high resolution and self-consistent libraries of model HII regions around central clusters of aging stars. The model libraries are constructed using the stellar population synthesis code, Starburst99, and photoionisation model, Cloudy. The Geneva and PARSEC stellar evolutionary models are plugged into Starburst99 to allow comparison between the two models. Using a spectrum-fitting methodology that allows the spectral features in the stellar and nebular continua (e.g. Wolf-Rayet features, Paschen jump), and emission-line diagnostics to constrain the models, we apply the libraries to the high-resolution MUSE spectra of the starbursting regions in the Antennae galaxy. Through this approach, we were able to model the continuum emission from Wolf-Rayet stars and extract stellar and gas metallicities, ages, electron temperatures and densities of starbursts by exploiting the full spectrum. From the application to the Antennae galaxy, we find that (1) the starbursts in the Antennae galaxy are characterised by stellar and gas metallicities of around solar, (2) the star-forming gas in starbursts in the Western loop of NGC 4038 appear to be more enriched, albeit slightly, than the rest of galaxy, (3) the youngest starbursts are found across the overlap region and over parts of the western-loop, though in comparison, the regions in the western-loop appear to be at a slightly later stage in star-formation than the overlap region, and (4) the results obtained from fitting the Geneva and Parsec models are largely consistent.


title: On the conditions for the breaking of oscillations in a cold plasma

abstract: The Cauchy problem for a quasilinear system of hyperbolic equations describing plane one-dimensional relativistic oscillations of electrons in a cold plasma is considered. For some simplified formulation of the problem, a criterion for the existence of a global in time solutions is obtained. For the original problem, a sufficient condition for the loss of smoothness is found, as well as a sufficient condition for the solution to remain smooth at least for time $ 2 \pi $. In addition, it is shown that in the general case, arbitrarily small perturbations of the trivial state lead to the formation of singularities in a finite time. It is further proved that there are special initial data such that the respective solution remains smooth for all time, even in the relativistic case. Periodic in space traveling wave gives an example of such a solution. In order for such a wave to be smooth, the velocity of the wave must be greater than a certain constant that depends on the initial data. Nevertheless, arbitrary small perturbation of general form destroys these global in time smooth solutions. The nature of the singularities of the solutions is illustrated by numerical examples.

categories: math-ph math.AP math.MP physics.plasm-ph


title: Peculiar solar sources and geospace disturbances on 20-26 August 2018

abstract: On the approach to minimum of Solar Cycle 24, on 26 August 2018, an unexpectedly strong geomagnetic storm (GMS) suddenly occurred. Its Dst index reached -174 nT, that is the third of the most intense storms during the cycle. The analysis showed that it was initiated by a two-step long filament eruption, which occurred on 20 August in the central sector of the solar disk. The eruptions were accompanied by two large-scale divergent ribbons and dimmings of a considerable size and were followed by relatively weak but evident Earth-directed coronal mass ejections. In the inner corona, their estimated speed was very low of about 200-360 km/s. The respective interplanetary transients apparently propagated between two high-speed solar wind streams originated from a two-component coronal hole and therefore their expansion was limited. The resulting ejecta arrived at the Earth only on 25 August and brought an unexpectedly strong field of Bt ~ 18.2 nT with a predominantly negative Bz component of almost the same strength. The geospace storm also manifested itself in the form of a peculiar Forbush decrease (FD). Its magnitude was about 1.5%, which is rather small for the observed G3-class GMS. The main unusual feature of the event is that large positive bursts with an enhancement up to 3% above the pre-event level were recorded on the FD background. We argue that these bursts were mainly caused by an unusually large and changeable cosmic ray anisotropy combined with lowering of the geomagnetic cutoff rigidity in the perturbed Earth's magnetosphere under conditions of the cycle minimum.


title: $\mathrm{SL}_2(\mathbb{R})$-dynamics on the moduli space of one-holed tori

abstract: We study the $\mathrm{SL}_2(\mathbb{R})$-action on the moduli space of (triangulable) dilation tori with one boundary component. We prove that every orbit is either closed or dense, and that every orbit of the Teichmuller flow escapes to infinity.


title: Twisted Dirac operator on quantum SU(2) in disc coordinates

abstract: The quantum disc is used to define a noncommutative analogue of a dense coordinate chart and of left-invariant vector fields on quantum SU(2). This yields two twisted Dirac operators for different twists that are related by a gauge transformation and have bounded twisted commutators with a suitable algebra of differentiable functions on quantum SU(2).


title: New constraints for nuclear parton distribution functions from hadron-nucleus collision processes

abstract: This work studies collinearly factorizable nuclear parton distribution functions (nPDFs) in perturbative Quantum Chromodynamics (QCD) at next-to-leading order in the light of hadron-nucleus collision data which have not been included in nPDF analyses previously. The aim is at setting new constraints on the nuclear modifications of the gluon distribution and on the flavour separation of quark nuclear modifications. The introductory part provides an outline of the theoretical framework of QCD collinear factorization and the used statistical methods and relates the work presented here to other similar contemporary analyses. As a result, a new set of nPDFs, EPPS16, is presented, including for the first time electroweak-boson and dijet production data from CERN-LHC proton-lead collisions and allowing a full flavour separation in the fit. The flavour separation is constrained with Drell-Yan dilepton-production data from fixed target pion-nucleus experiments and neutrino-nucleus deep-inelastic scattering data, which are shown to give evidence for the similarity of the u and d valence-quark nuclear modifications. For studying the gluon degrees of freedom, collider data are essential and in the EPPS16 analysis new constraints are derived from the dijet production at the LHC. Possible further constraints for the gluons are investigated in terms of the LHC data on nuclear modification ratios of dijet and D-meson production. Using a non-quadratically improved Hessian reweighting method, these measurements are found to put stringent constraints on the gluon modifications in the lead nucleus, reaching smaller values of the nucleon momentum fraction than previously accessible. A study on the future prospects of constraining nPDFs within a multi-observable approach with the BNL-RHIC is also given.


title: Generalized Perron Roots and Solvability of the Absolute Value Equation

abstract: Let $A$ be a $n\times n$ real matrix. The piecewise linear equation system $z-A\vert z\vert =b$ is called an absolute value equation (AVE). It is well-known to be equivalent to the linear complementarity problem. Unique solvability of the AVE is known to be characterized in terms of a generalized Perron root called the sign-real spectral radius of $A$. For mere, possibly non-unique, solvability no such characterization exists. We narrow this gap in the theory. That is, we define the concept of the aligned spectrum of $A$ and prove, under some mild genericity assumptions on $A$, that the mapping degree of the piecewise linear function $F_A:\mathbb{R}^n\to\mathbb{R}^n\,, z\mapsto z-A\lvert z\rvert$ is congruent to $(k+1)\mod 2$, where $k$ is the number of aligned values of $A$ which are larger than $1$. We also derive an exact--but more technical--formula for the degree of $F_A$ in terms of the aligned spectrum. Finally, we derive the analogous quantities and results for the LCP.


title: Exploring Bistability in the Cycles of the Solar Dynamo through Global Simulations

abstract: The calling card of solar magnetism is the sunspot cycle, during which sunspots regularly reverse their polarity sense every 11 years. However, a number of more complicated time-dependent behaviors have also been identified. In particular, there are temporal modulations associated with active longitudes and hemispheric asymmetry, when sunspots appear at certain solar longitudes or else in one hemisphere preferentially. So far, a direct link between between this asymmetric temporal behavior and the underlying solar dynamo has remained elusive. In this work, we present results from global, 3D magnetohydrodynamic (MHD) simulations, which for the first time display both behavior reminiscent of the sunspot cycle (regular polarity reversals and equatorward migration of internal magnetic field) and asymmetric, irregular behavior that in the simulations we interpret as active longitudes and hemispheric asymmetry. The simulations are thus bistable, in that the turbulent convection can stably support two distinct flavors of magnetism at different times, in superposition, or with smooth transitions from one state to the other. We discuss this new family of dynamo models in the context of the extensive observations of the Sun's surface magnetic field with the Solar and Heliospheric Observatory (SOHO) and the Solar Dynamics Observatory (SDO), as well as earlier observations of sunspot number and synoptic maps. We suggest that the solar dynamo itself may be bistable in nature, exhibiting two types of temporal behavior in the magnetic field.


title: A perturbative algorithm for rotational decoherence

abstract: Recent advances in levitated optomechanics provide new perspectives for the use of rotational degrees of freedom for the development of quantum technologies as well as for testing fundamental physics. As for the translational case, their use, especially in the quantum regime, is limited by environmental noises, whose characterization is fundamental in order to assess, control and minimize their effect, in particular decoherence. Here, we present a general perturbative approach to compute decoherence for a quantum system in a superposition of its rotational degrees of freedom. The specific cases of the dipole-dipole and quadrupole-quadrupole interactions are solved explicitly, and we show that the rotational degrees of freedom decohere on a time scale that can be longer than the translational one.


title: Quantitative Theory of Triplet Pairing in the Unconventional Superconductor LaNiGa$_2$

abstract: The exceptionally low-symmetry crystal structures of the time-reversal symmetry breaking superconductors LaNiC$_2$ and LaNiGa$_2$ lead to an internally-antisymmetric non-unitary triplet (INT) state as the only possibility compatible with experiments. We argue that this state has a distinct signature: a double-peak structure in the Density of States (DOS) which resolves in the spin channel in a particular way. We construct a detailed model of LaNiGa$_2$ capturing its electronic band structure and magnetic properties ab initio. The pairing mechanism is described via a single adjustable parameter. The latter is fixed by the critical temperature $T_c$ allowing parameter-free predictions. We compute the electronic specific heat and find excellent agreement with experiment. The size of the ordered moment in the superconducting state is compatible with zero-field muon spin relaxation experiments and the predicted spin-resolved DOS suggests the spin-splitting is within the reach of present experimental technology.


title: Second powers of cover ideals of paths

abstract: We show that the second power of the cover ideal of a path graph has linear quotients. To prove our result we construct a recursively defined order on the generators of the ideal which yields linear quotients. Our construction has a natural generalization to the larger class of chordal graphs. This generalization allows us to raise some questions that are related to some open problems about powers of cover ideals of chordal graphs.


title: Efficiency of Observed Information Adaptive Designs

abstract: In this work the primary objective is to maximize the precision of the maximum likelihood estimate in a linear regression model through the efficient design of the experiment. One common measure of precision is the unconditional mean square error. Unconditional mean square error has been a primary motivator for optimal designs; commonly, defined as the design that maximizes a concave function of the expected Fisher information. The inverse of expected Fisher information is asymptotically equal to the mean square error of the maximum likelihood estimate. There is a substantial amount of existing literature that argues the mean square error conditioned on an appropriate ancillary statistic better represents the precision of the maximum likelihood estimate. Despite evidence in favor of conditioning, limited effort has been made to find designs that are optimal with respect to conditional mean square error. The inverse of observed Fisher information is a higher order approximation of the conditional mean square error than the inverse of expected Fisher information [Efron and Hinkley (1978)]. In light of this, a more relevant objective is to find designs that optimize observed Fisher information. Unlike expected Fisher information, observed Fisher information depends on the observed data and cannot be used to design an experiment completely in advance of data collection. In a sequential experiment the observed Fisher information from past observations is available to inform the design of the next observation. In this work an adaptive design that incorporates observed Fisher information is proposed for linear regression models. It is shown that the proposed design is more efficient, at the limit, than any fixed design, including the optimal design, with respect to conditional mean square error.


title: Simple mathematical model of aging

abstract: A simple mathematical model of the aging process for long-lived organisms is considered. The key point in this model is the assumption that the body does not have internal clocks that count out the chronological time at scales of decades. At these scales, we may limit ourselves by empirical consideration only the background (smoothed, averaged) processes. The body is dealing with internal biological factors, which can be considered as the biological clocks in suitable parameterization of corresponding variables. The dynamics of these variables is described using a system of autonomous ODEs. A particular representation of the right-hand side of equations in the form of quadratic polynomials is considered. In the simplest case of one variable we deal with a logistic equation, which has an analytical solution. Such quadratic model is justified if it is used to predict the dynamic of aging process for relatively small time intervals. However, since a well-defined biological interpretation can be given for the quadratic right-hand side of the equations, we can expect that the area of applicability of this simplified empirical model can extend to relatively large time intervals. The considered model is, in our opinion, a good and simple mathematical framework for organizing experimental data. The model may be useful to quantify and arrange the chains of causal relationships that may appear significant for the development of the aging process.


title: Weakly compact sets and weakly compact pointwise multipliers in Banach function lattices

abstract: We prove that the class of Banach function lattices in which all relatively weakly compact sets are equi-integrable sets (i.e. spaces satisfying the Dunford-Pettis criterion) coincides with the class of 1-disjointly homogeneous Banach lattices. A new examples of such spaces are provided. Furthermore, it is shown that Dunford-Pettis criterion is equivalent to de la Vallee Poussin criterion in all rearrangement invariant spaces on the interval. Finally, the results are applied to characterize weakly compact pointwise multipliers between Banach function lattices.


title: Cyanure: An Open-Source Toolbox for Empirical Risk Minimization for Python, C++, and soon more

abstract: Cyanure is an open-source C++ software package with a Python interface. The goal of Cyanure is to provide state-of-the-art solvers for learning linear models, based on stochastic variance-reduced stochastic optimization with acceleration mechanisms. Cyanure can handle a large variety of loss functions (logistic, square, squared hinge, multinomial logistic) and regularization functions (l_2, l_1, elastic-net, fused Lasso, multi-task group Lasso). It provides a simple Python API, which is very close to that of scikit-learn, which should be extended to other languages such as R or Matlab in a near future.


title: APRICOT: A Dataset of Physical Adversarial Attacks on Object Detection

abstract: Physical adversarial attacks threaten to fool object detection systems, but reproducible research on the real-world effectiveness of physical patches and how to defend against them requires a publicly available benchmark dataset. We present APRICOT, a collection of over 1,000 annotated photographs of printed adversarial patches in public locations. The patches target several object categories for three COCO-trained detection models, and the photos represent natural variation in position, distance, lighting conditions, and viewing angle. Our analysis suggests that maintaining adversarial robustness in uncontrolled settings is highly challenging, but it is still possible to produce targeted detections under white-box and sometimes black-box settings. We establish baselines for defending against adversarial patches through several methods, including a detector supervised with synthetic data and unsupervised methods such as kernel density estimation, Bayesian uncertainty, and reconstruction error. Our results suggest that adversarial patches can be effectively flagged, both in a high-knowledge, attack-specific scenario, and in an unsupervised setting where patches are detected as anomalies in natural images. This dataset and the described experiments provide a benchmark for future research on the effectiveness of and defenses against physical adversarial objects in the wild.


title: Towards personalized diagnosis of Glioblastoma in Fluid-attenuated inversion recovery (FLAIR) by topological interpretable machine learning

abstract: Glioblastoma multiforme (GBM) is a fast-growing and highly invasive brain tumour, it tends to occur in adults between the ages of 45 and 70 and it accounts for 52 percent of all primary brain tumours. Usually, GBMs are detected by magnetic resonance images (MRI). Among MRI, Fluid-attenuated inversion recovery (FLAIR) sequence produces high quality digital tumour representation. Fast detection and segmentation techniques are needed for overcoming subjective medical doctors (MDs) judgment. In the present investigation, we intend to demonstrate by means of numerical experiments that topological features combined with textural features can be enrolled for GBM analysis and morphological characterization on FLAIR. To this extent, we have performed three numerical experiments. In the first experiment, Topological Data Analysis (TDA) of a simplified 2D tumour growth mathematical model had allowed to understand the bio-chemical conditions that facilitate tumour growth: the higher the concentration of chemical nutrients the more virulent the process. In the second experiment topological data analysis was used for evaluating GBM temporal progression on FLAIR recorded within 90 days following treatment (e.g., chemo-radiation therapy - CRT) completion and at progression. The experiment had confirmed that persistent entropy is a viable statistics for monitoring GBM evolution during the follow-up period. In the third experiment we had developed a novel methodology based on topological and textural features and automatic interpretable machine learning for automatic GBM classification on FLAIR. The algorithm reached a classification accuracy up to the 97%.


title: Differentiable programming and its applications to dynamical systems

abstract: Differentiable programming is the combination of classical neural networks modules with algorithmic ones in an end-to-end differentiable model. These new models, that use automatic differentiation to calculate gradients, have new learning capabilities (reasoning, attention and memory). In this tutorial, aimed at researchers in nonlinear systems with prior knowledge of deep learning, we present this new programming paradigm, describe some of its new features such as attention mechanisms, and highlight the benefits they bring. Then, we analyse the uses and limitations of traditional deep learning models in the modeling and prediction of dynamical systems. Here, a dynamical system is meant to be a set of state variables that evolve in time under general internal and external interactions. Finally, we review the advantages and applications of differentiable programming to dynamical systems.


title: $O(N)$ Models with Boundary Interactions and their Long Range Generalizations

abstract: We study the critical properties of scalar field theories in $d+1$ dimensions with $O(N)$ invariant interactions localized on a $d$-dimensional boundary. By a combination of large $N$ and epsilon expansions, we provide evidence for the existence of non-trivial $O(N)$ BCFTs in $1<d<4$. Due to having free fields in the bulk, these models possess bulk higher-spin currents which are conserved up to terms localized on the boundary. We suggest that this should lead to a set of protected spinning operators on the boundary, and give evidence that their anomalous dimensions vanish. We also discuss the closely related long-range $O(N)$ models in $d$ dimensions, and in particular study a weakly coupled description of the $d=1$ long range $O(N)$ model near the upper critical value of the long range parameter, which is given in terms of a non-local non-linear sigma model. By combining the known perturbative descriptions, we provide some estimates of critical exponents in $d=1$.


title: Consensus seeking gradient descent flows on boundaries of convex sets

abstract: Consensus on nonlinear spaces is of use in many control applications. This paper proposes a gradient descent flow algorithm for consensus on hypersurfaces. We show that if an inequality holds, then the system converges for almost all initial conditions and all connected graphs. The inequality involves the hypersurface Gauss map and the gradient and Hessian of the implicit equation. Moreover, for the inequality to hold, it is necessary that the manifold is the boundary of a convex set. The literature already contains an algorithm for consensus on hypersurfaces. That algorithm on any ellipsoid is equivalent to our algorithm on the unit sphere. In particular, that algorithm achieves almost global synchronization on ellipsoids. These findings suggest that strong convergence results for consensus seeking gradient descent flows may be established on manifolds that are the boundaries of convex sets.


title: Two-sided optimal stopping for L\'evy processes

abstract: Infinite horizon optimal stopping problems for a L\'evy processes with a two-sided reward function are considered. A two-sided verification theorem is presented in terms of the overall supremum and the overall infimum of the process. A result to compute the angle of the value function at the optimal thresholds of the stopping region is given. To illustrate the results, the optimal stopping problem of a compound Poisson process with two-sided exponential jumps and a two-sided payoff function is solved. In this example, the smooth-pasting condition does not hold.


title: Shapes of fluid membranes with chiral edges

abstract: We carry out Monte Carlo simulations of a colloidal fluid membrane composed of chiral rod-like viruses. The membrane is modeled by a triangular mesh of beads connected by bonds in which the bonds and beads are free to move at each Monte Carlo step. Since the constituent viruses are experimentally observed to twist only near the membrane edge, we use an effective energy that favors a particular sign of the geodesic torsion of the edge. The effective energy also includes membrane bending stiffness, edge bending stiffness, and edge tension. We find three classes of membrane shapes resulting from the competition of the various terms in the free energy: branched shapes, chiral disks, and vesicles. Increasing the edge bending stiffness smooths the membrane edge, leading to correlations among the membrane normal at different points along the edge. We also consider membrane shapes under an external force by fixing the distance between two ends of the membrane, and find the shape for increasing values of the distance between the two ends. As the distance increases, the membrane twists into a ribbon, with the force eventually reaching a plateau.


title: Function Approximation via The Subsampled Poincar\' e Inequality

abstract: Function approximation and recovery via some sampled data have long been studied in a wide array of applied mathematics and statistics fields. Analytic tools, such as the Poincar\'e inequality, have been handy for estimating the approximation errors in different scales. The purpose of this paper is to study a generalized Poincar\' e inequality, where the measurement function is of subsampled type, with a small but non-zero lengthscale that will be made precise. Our analysis identifies this inequality as a basic tool for function recovery problems. We discuss and demonstrate the optimality of the inequality concerning the subsampled lengthscale, connecting it to existing results in the literature. In application to function approximation problems, the approximation accuracy using different basis functions and under different regularity assumptions is established by using the subsampled Poincar\'e inequality. We observe that the error bound blows up as the subsampled lengthscale approaches zero, due to the fact that the underlying function is not regular enough to have well-defined pointwise values. A weighted version of the Poincar\' e inequality is proposed to address this problem; its optimality is also discussed.


title: Effect of Electromagnetic Interaction on Galactic Center Flare Components

abstract: Recently, near-infrared GRAVITY@ESO observations at $2.2\,\mu{\rm m}$ have announced the detection of three bright "flares" in the vicinity of the Galactic center supermassive black hole (SMBH) that exhibited orbital motion at a distance of about $6 - 11$ gravitational radii from an $\sim 4\times 10^6\, M_{\odot}$ black hole. There are indications of the presence of a large-scale, organized component of the magnetic field at the Galactic center. Electromagnetic effects on the flare dynamics were previously not taken into account despite the relativistic motion of a plasma in magnetic field leading to the charge separation and nonnegligible net charge density in the plasma. Applying various approaches, we find the net charge number density of the flare components of the order of $10^{-3} - 10^{-4}$ cm$^{-3}$, while the particles' total number density is of the order of $10^{6} - 10^{8}$ cm$^{-3}$. However, even such a tiny excess of charged particles in the quasi-neutral plasma can significantly affect the dynamics of flare components, which can then lead to the degeneracy in the measurements of spin of the SMBH. Analyzing the dynamics of recent flares in the case of the rapidly rotating black hole, we also constrain the inclination angle between the magnetic field and spin axis to $\alpha < 50^{\circ}$, as for larger angles, the motion of the hot spot is strongly chaotic.


title: Towards a quantum Oppenheimer-Snyder model

abstract: We present a consistent canonical formulation of the flat Oppenheimer-Snyder model, including the Schwarzschild exterior. The switching between comoving and stationary observer is realized by promoting the coordinate transformation between dust proper time and Schwarzschild-Killing time to a canonical one. This leads to two different forms of the Hamiltonian constraint, both (almost) deparameterizable with regard to one of these times. A preliminary quantization of these constraints reveals a consistent picture for both observers: the singularity is avoided by a bounce.


title: Density scaling of generalized Lennard-Jones fluids in different dimensions

abstract: Liquids displaying strong virial-potential energy correlations conform to an approximate density scaling of their structural and dynamical observables. This scaling property does not extend to the entire phase diagram, in general. The validity of the scaling can be quantified by a correlation coefficient. In this work a simple scheme to predict the correlation coefficient and the density-scaling exponent is presented. Although this scheme is exact only in the dilute gas regime or in high dimension d, a comparison with results from molecular dynamics simulations in d = 1 to 4 shows that it reproduces well the behavior of generalized Lennard-Jones systems in a large portion of the fluid phase.


title: Lift & Learn: Physics-informed machine learning for large-scale nonlinear dynamical systems

abstract: We present Lift & Learn, a physics-informed method for learning low-dimensional models for large-scale dynamical systems. The method exploits knowledge of a system's governing equations to identify a coordinate transformation in which the system dynamics have quadratic structure. This transformation is called a lifting map because it often adds auxiliary variables to the system state. The lifting map is applied to data obtained by evaluating a model for the original nonlinear system. This lifted data is projected onto its leading principal components, and low-dimensional linear and quadratic matrix operators are fit to the lifted reduced data using a least-squares operator inference procedure. Analysis of our method shows that the Lift & Learn models are able to capture the system physics in the lifted coordinates at least as accurately as traditional intrusive model reduction approaches. This preservation of system physics makes the Lift & Learn models robust to changes in inputs. Numerical experiments on the FitzHugh-Nagumo neuron activation model and the compressible Euler equations demonstrate the generalizability of our model.


title: AeroRIT: A New Scene for Hyperspectral Image Analysis

abstract: We investigate applying convolutional neural network (CNN) architecture to facilitate aerial hyperspectral scene understanding and present a new hyperspectral dataset-AeroRIT-that is large enough for CNN training. To date the majority of hyperspectral airborne have been confined to various sub-categories of vegetation and roads and this scene introduces two new categories: buildings and cars. To the best of our knowledge, this is the first comprehensive large-scale hyperspectral scene with nearly seven million pixel annotations for identifying cars, roads, and buildings. We compare the performance of three popular architectures - SegNet, U-Net, and Res-U-Net, for scene understanding and object identification via the task of dense semantic segmentation to establish a benchmark for the scene. To further strengthen the network, we add squeeze and excitation blocks for better channel interactions and use self-supervised learning for better encoder initialization. Aerial hyperspectral image analysis has been restricted to small datasets with limited train/test splits capabilities and we believe that AeroRIT will help advance the research in the field with a more complex object distribution to perform well on. The full dataset, with flight lines in radiance and reflectance domain, is available for download at https://github.com/aneesh3108/AeroRIT. This dataset is the first step towards developing robust algorithms for hyperspectral airborne sensing that can robustly perform advanced tasks like vehicle tracking and occlusion handling.


title: Path integral quantization of a spinning particle

abstract: Following the idea of Alekseev and Shatashvili we derive the path integral quantization of a modified relativistic particle action that results in the Feynman propagator of a free field with arbitrary spin. This propagator can be associated with the Duffin, Kemmer, and Petiau (DKP) form of a free field theory. We show explicitly that the obtained DKP propagator is equivalent to the standard one, for spins 0 and 1. We argue that this equivalence holds also for higher spins.


title: Deep Radar Waveform Design for Efficient Automotive Radar Sensing

abstract: In radar systems, unimodular (or constant-modulus) waveform design plays an important role in achieving better clutter/interference rejection, as well as a more accurate estimation of the target parameters. The design of such sequences has been studied widely in the last few decades, with most design algorithms requiring sophisticated a priori knowledge of environmental parameters which may be difficult to obtain in real-time scenarios. In this paper, we propose a novel hybrid model-driven and data-driven architecture that adapts to the ever changing environment and allows for adaptive unimodular waveform design. In particular, the approach lays the groundwork for developing extremely low-cost waveform design and processing frameworks for radar systems deployed in autonomous vehicles. The proposed model-based deep architecture imitates a well-known unimodular signal design algorithm in its structure, and can quickly infer statistical information from the environment using the observed data. Our numerical experiments portray the advantages of using the proposed method for efficient radar waveform design in time-varying environments.


title: A remark on sets with few distances in $\mathbb{R}^{d}$

abstract: A celebrated theorem due to Bannai-Bannai-Stanton says that if $A$ is a set of points in $\mathbb{R}^{d}$, which determines $s$ distinct distances, then $$|A| \leq {d+s \choose s}.$$ In this note, we give a new simple proof of this result by combining Sylvester's Law of Inertia for quadratic forms with the proof of the so-called Croot-Lev-Pach Lemma from additive combinatorics.


title: Probing the Sivers function with an unpolarized target: GTMD distributions and the Odderons

abstract: It is commonly believed that the Sivers function has uniquely to do with processes involving a transversely polarized nucleon. In this paper we show that it is not necessarily the case. We demonstrate that exclusive pion production in $un$polarized electron-proton scattering in the forward region is a direct probe of the gluon Sivers function due to its connection to the QCD Odderon.


title: Mean field flow equations and asymptotically free scalar fields

abstract: The flow equations of the renormalisation group permit to analyse the perturbative $n$-point functions of renormalisable quantum field theories. Rigorous bounds implying renormalisablility allow to control large momentum behaviour, infrared singularities and large order behaviour in the number of loops and the number of arguments $n\,$. Gauge symmetry which is broken by the flow in momentum or position space, can be shown to be restored in the renormalised theory. In this paper we want to do a first but important step towards a rigorous nonperturbative analysis of the flow equations (FEs). We restrict to massive scalar fields and analyse the {\it mean field limit} where the Schwinger or 1PI functions are considered to be momentum independent or, otherwise stated, are replaced by their zero momentum values. We regard smooth solutions of the system of FEs for the $n$-point functions for different sets of boundary conditions. We will realise that allowing for non\-vanishing irrelevant terms permits to construct {\it asymptotically free} and thus nontrivial {\it scalar field theories} in the mean field approximation. We will also analyse the so-called trivial solution so far generally believed to exhaust four-dimensional scalar field theory. The method paves the way to a study of the system of FEs beyond the mean field limit.


title: On Fano threefolds with $\mathbb{C}^*$-action

abstract: We study three-dimensional Fano varieties with $\mathbb{C}^*$-action. Complementing recent results [13], we give classification results in the canonical case, where the maximal orbit quotient is $\mathbb{P}_2$ having a line arrangement of five lines in special position as its critical values.


title: A classification of the finite non-solvable minimal non-CA-groups

abstract: A group is called a CA-group if the centralizer of every non-central element is abelian. Furthermore, a group is called a minimal non-CA-group if it is not a CA-group itself, but all of its proper subgroups are. In this paper, we give a classification of the finite non-solvable minimal non-CA-groups.


title: Short proof that Kneser graphs are Hamiltonian for $n\geq 4k$

abstract: For integers $n\geq k\geq 1$, the Kneser graph $K(n,k)$ is the graph with vertex set $V=[n]^{(k)}$ and edge set $E=\{\{x,y\} \in V^{(2)}: x\cap y=\emptyset\}$. Chen proved that for $n\geq 3k$, Kneser graphs are Hamiltonian and later improved this to $n\geq 2.62k+1$. Furthermore, Chen and F\"uredi gave a short proof that if $k | n$, Kneser graphs are Hamiltonian for $n\geq 3k$. In this note, we present a short proof that does not need the divisibility condition, i.e., we give a short proof that $K(n,k)$ is Hamiltonian for $n\geq 4k$.


title: Underlying event measurements in $p$+$p$ collisions at $\sqrt{s}= 200 $ GeV at RHIC

abstract: Particle production sensitive to non-factorizable and non-perturbative processes that contribute to the underlying event associated with a high transverse momentum ($p_{T}$) jet in proton+proton collisions at $\sqrt{s}$=200 GeV is studied with the STAR detector. Each event is divided into three regions based on the azimuthal angle with respect to the highest-$p_{T}$ jet direction: in the leading jet direction ("Toward"), opposite to the leading jet ("Away"), and perpendicular to the leading jet ("Transverse"). In the Transverse region, the average charged particle density is found to be between 0.4 and 0.6 and the mean transverse momentum, $\langle p_{T}\rangle$, between 0.5-0.7 GeV/$c$ for particles with $p_{T}$$>$0.2 GeV/$c$ at mid-pseudorapidity ($|\eta|$$<$1) and jet $p_{T}$$>$15 GeV/$c$. Both average particle density and $\langle p_{T}\rangle$ depend weakly on the leading jet $p_{T}$. Closer inspection of the Transverse region hints that contributions to the underlying event from initial- and final-state radiation are significantly smaller in these collisions than at the higher energies, up to 13 TeV, recorded at the LHC. Underlying event measurements associated with a high-$p_{T}$ jet will contribute to our understanding of QCD processes at hard and soft scales at RHIC energies, as well as provide constraints to modeling of underlying event dynamics.


title: Opening the 1$\,$Hz axion window

abstract: An axion-like particle (ALP) with mass $m_\phi \sim 10^{-15}$eV oscillates with frequency $\sim$1 Hz. This mass scale lies in an open window of astrophysical constraints, and appears naturally as a consequence of grand unification (GUT) in string/M-theory. However, with a GUT-scale decay constant such an ALP overcloses the Universe, and cannot solve the strong CP problem. In this paper, we present a two axion model in which the 1 Hz ALP constitutes the entirety of the dark matter (DM) while the QCD axion solves the strong CP problem but contributes negligibly to the DM relic density. The mechanism to achieve the correct relic densities relies on low-scale inflation ($m_\phi \lesssim H_{\rm inf}\lesssim 1$ MeV), and we present explicit realisations of such a model. The scale in the axion potential leading to the 1 Hz axion generates a value for the strong CP phase which oscillates around $\bar{\theta}_{\rm QCD}\sim 10^{-12}$, within reach of the proton storage ring electric dipole moment experiment. The 1 Hz axion is also in reach of near future laboratory and astrophysical searches.


title: Learning from Discriminatory Training Data

abstract: Supervised learning systems are trained using historical data and, if the data was tainted by discrimination, they may unintentionally learn to discriminate against protected groups. We propose that fair learning methods, despite training on potentially discriminatory datasets, shall perform well on fair test datasets. Such dataset shifts crystallize application scenarios for specific fair learning methods. For instance, the removal of direct discrimination can be represented as a particular dataset shift problem. For this scenario, we propose a learning method that provably minimizes model error on fair datasets, while blindly training on datasets poisoned with direct additive discrimination. The method is compatible with existing legal systems and provides a solution to the widely discussed issue of protected groups' intersectionality by striking a balance between the protected groups. Technically, the method applies probabilistic interventions, has causal and counterfactual formulations, and is computationally lightweight - it can be used with any supervised learning model to prevent discrimination via proxies while maximizing model accuracy for business necessity.


title: Violent Buckling Benefits Galactic Bars

abstract: Galactic bars are unstable to a vertical buckling instability which heats the disk and in some cases forms a boxy/peanut shaped bulge. We analyze the buckling instability as an application of classical Euler buckling followed by nonlinear gravitational Landau damping in the collisionless system. We find that the buckling instability is dictated by the kinematic properties and geometry of the bar. The analytical result is compared to simulations of isolated galaxies containing the disk and dark matter components. Our results demonstrate that violent buckling does not destroy bars while a less energetic buckling can dissolve the bar. The disks that undergo gentle buckling remain stable to bar formation which may explain the observed bar fraction in the local universe. Our results align with the results from recent surveys.


title: Securing Big Data systems, A cybersecurity management discussion

abstract: This paper explores the essential areas of cybersecurity management for big data systems. Big data platform stems its complexity from being a collection of interrelated non-standardized systems that interact with each other to process large data-sets. This complexity increases the chances of overlooking critical cybersecurity management practices. The paper discusses various security management for each part of big data systems. This includes security measures, standards, and methodologies. The primary objective is to highlight the essential segments of security management that are expected to be implemented flawlessly in big data systems.


title: First-Passage Percolation under extreme disorder: from bond-percolation to Kardar-Parisi-Zhang universality

abstract: We consider the statistical properties of arrival times and balls on first-passage percolation (FPP) $2D$ square lattices with strong disorder in the link-times. A previous work showed a crossover in the weak disorder regime, between Gaussian and Kardar-Parisi-Zhang (KPZ) universality, with the crossover length decreasing as the noise amplitude grows. On the other hand, this work presents a very different behavior in the strong-disorder regime. A new crossover length appears below which the model is described by bond-percolation universality class. This characteristic length scale grows with the noise amplitude and diverges at the infinite-disorder limit. We provide a thorough characterization of the bond-percolation phase, reproducing its associated critical exponents through a careful scaling analysis of the balls. This is carried out through to a continuous mapping of the FPP passage time into the occupation probability of the bond-percolation problem. Moreover, the crossover length can be explained merely in terms of properties of the link-time distribution. The interplay between the new characteristic length and the correlation length intrinsic to bond-percolation determines the crossover between the initial percolation-like growth and theasymptotic KPZ scaling.


title: PointRend: Image Segmentation as Rendering

abstract: We present a new method for efficient high-quality image segmentation of objects and scenes. By analogizing classical computer graphics methods for efficient rendering with over- and undersampling challenges faced in pixel labeling tasks, we develop a unique perspective of image segmentation as a rendering problem. From this vantage, we present the PointRend (Point-based Rendering) neural network module: a module that performs point-based segmentation predictions at adaptively selected locations based on an iterative subdivision algorithm. PointRend can be flexibly applied to both instance and semantic segmentation tasks by building on top of existing state-of-the-art models. While many concrete implementations of the general idea are possible, we show that a simple design already achieves excellent results. Qualitatively, PointRend outputs crisp object boundaries in regions that are over-smoothed by previous methods. Quantitatively, PointRend yields significant gains on COCO and Cityscapes, for both instance and semantic segmentation. PointRend's efficiency enables output resolutions that are otherwise impractical in terms of memory or computation compared to existing approaches. Code has been made available at https://github.com/facebookresearch/detectron2/tree/master/projects/PointRend.


title: Comment on "Nonlocal Coherent Perfect Absorption"

abstract: We show that, contrary to what is claimed in Ref. [J. Jeffers, Phys. Rev. Lett. 123, 143602 (2019)], there is no nonlocal coherent absorption by the two lossy beam splitters considered in detail in the treatment of the setup presented in this reference. What happens is a local coherent absorption by the other two lossy beam splitters of the setup.


title: Learning Generalizable Visual Representations via Interactive Gameplay

abstract: A growing body of research suggests that embodied gameplay, prevalent not just in human cultures but across a variety of animal species including turtles and ravens, is critical in developing the neural flexibility for creative problem solving, decision making, and socialization. Comparatively little is known regarding the impact of embodied gameplay upon artificial agents. While recent work has produced agents proficient in abstract games, these environments are far removed from the real world and thus these agents can provide little insight into the advantages of embodied play. Hiding games, such as hide-and-seek, played universally, provide a rich ground for studying the impact of embodied gameplay on representation learning in the context of perspective taking, secret keeping, and false belief understanding. Here we are the first to show that embodied adversarial reinforcement learning agents playing Cache, a variant of hide-and-seek, in a high fidelity, interactive, environment, learn generalizable representations of their observations encoding information such as object permanence, free space, and containment. Moving closer to biologically motivated learning strategies, our agents' representations, enhanced by intentionality and memory, are developed through interaction and play. These results serve as a model for studying how facets of vision develop through interaction, provide an experimental framework for assessing what is learned by artificial agents, and demonstrates the value of moving from large, static, datasets towards experiential, interactive, representation learning.


title: Relativistic non-Fermi liquid from interacting birefringent fermions: A robust superuniversality

abstract: We address the emergent quantum critical phenomena for (pseudo)spin-3/2 birefringent fermions, featuring two effective Fermi velocities, when they reside close to itinerant Mott transitions realized through spontaneous symmetry breaking and triggered by strong local or Hubbardlike repulsive interactions. Irrespective of the nature of the mass orderings that produce fully gapped quasiparticle spectra in the ordered phase, which otherwise can be grouped into three classes, the system always possesses a \emph{unique} terminal velocity near the corresponding quantum critical point. The associated critical regime accommodates a relativistic non-Fermi liquid of strongly coupled collective bosonic and spin-1/2 Dirac excitations with vanishing weight of the quasiparticle pole. These conclusions are also operative near superconducting critical points. Therefore, relativistic non-Fermi liquid possibly constitutes a robust superuniversal description for the entire family of strongly correlated arbitrary half-integer spin Dirac materials.


title: Lightweight and Robust Representation of Economic Scales from Satellite Imagery

abstract: Satellite imagery has long been an attractive data source that provides a wealth of information on human-inhabited areas. While super resolution satellite images are rapidly becoming available, little study has focused on how to extract meaningful information about human habitation patterns and economic scales from such data. We present READ, a new approach for obtaining essential spatial representation for any given district from high-resolution satellite imagery based on deep neural networks. Our method combines transfer learning and embedded statistics to efficiently learn critical spatial characteristics of arbitrary size areas and represent them into a fixed-length vector with minimal information loss. Even with a small set of labels, READ can distinguish subtle differences between rural and urban areas and infer the degree of urbanization. An extensive evaluation demonstrates the model outperforms the state-of-the-art in predicting economic scales, such as population density for South Korea (R^2=0.9617), and shows a high potential use for developing countries where district-level economic scales are not known.


title: srlearn: A Python Library for Gradient-Boosted Statistical Relational Models

abstract: We present srlearn, a Python library for boosted statistical relational models. We adapt the scikit-learn interface to this setting and provide examples for how this can be used to express learning and inference problems.


title: Lieb, Entropy and Logarithmic uncertainty principles for the multivariate continuous quaternion Shearlet Transform

abstract: In this paper, we generalize the continuous quaternion shearlet transform on $\mathbb{R}^{2}$ to $\mathbb{R}^{2d}$, called the multivariate two sided continuous quaternion shearlet transform. Using the two sided quaternion Fourier transform, we derive several important properties such as (reconstruction formula, reproducing kernel, plancherel's formula, etc.). We present several example of the multivariate two sided continuous quaternion shearlet transform. We apply the multivariate two sided continuous quaternion shearlet transform properties and the two sided quaternion Fourier transform to establish Lieb uncertainty principle and the Logarithmic uncertainty principle. Last we study the Beckner's uncertainty principle in term of entropy for the multivariate two sided continuous quaternion shearlet transform.


title: Visualisation of Brain Statistics with R-packages ggseg and ggseg3d

abstract: There is an increased emphasis on visualizing neuroimaging results in more intuitive ways. Common statistical tools for dissemination, such as bar charts, lack the spatial dimension that is inherent in neuroimaging data. Here we present two packages for the statistical software R, ggseg and ggseg3d, that integrate this spatial component. The ggseg and ggseg3d packages visualize pre-defined brain segmentations as both 2D polygons and 3D meshes, respectively. Both packages are integrated with other well-established R-packages, allowing great flexibility. In this tutorial, we present the main data and functions in the ggseg and ggseg3d packages for brain atlas visualization. The main highlighted functions are able to display brain segmentation plots in R. Further, the accompanying ggsegExtra-package includes a wider collection of atlases, and is intended for community-based efforts to develop more compatible atlases to ggseg and ggseg3d. Overall, the ggseg-packages facilitate parcellation-based visualizations in R, improve and ease the dissemination of the results, and increase the efficiency of the workflows.


title: Microwave power absorption to high energy electrons in the ECR ion thruster

abstract: The microwave power absorption efficiency of the mu10 ECR ion thruster, utilized in the Japanese asteroid explorers Hayabusa and Hayabusa2, is investigated in order to allow performance measurement and provide information for its improvement. A model detailing the local electron behavior in a real ECR plasma discharge, based the magnetic field characteristics, is presented. Three methods to evaluate the microwave power absorption efficiency are proposed: an estimation based on the chamber geometry and magnetic field characteristics, a measuremen based on performance parameters and a measurement performed with Langmuir probes. The equations used for each method are analytically derived. The local electron behavior model is confirmed with a Langmuir probe experiment. Measurement of the microwave power absorption efficiency is performed with the two independent methods proposed. Results from the two experiments show good agreement with each other and with the theory. Finally, a diffusion model explaining the different electron temperature distributions observed in the chamber is proposed. The model and experiments clarify the physics behind previously observed performance variations and give valuable hints for future chamber improvement.


title: Extrinsic Kernel Ridge Regression Classifier for Planar Kendall Shape Space

abstract: Kernel methods have had great success in Statistics and Machine Learning. Despite their growing popularity, however, less effort has been drawn towards developing kernel based classification methods on Riemannian manifolds due to difficulty in dealing with non-Euclidean geometry. In this paper, motivated by the extrinsic framework of manifold-valued data analysis, we propose a new positive definite kernel on planar Kendall shape space $\Sigma_2^k$, called extrinsic Veronese Whitney Gaussian kernel. We show that our approach can be extended to develop Gaussian kernels on any embedded manifold. Furthermore, kernel ridge regression classifier (KRRC) is implemented to address the shape classification problem on $\Sigma_2^k$, and their promising performances are illustrated through the real data analysis.


title: Three dimensional waveguide-interconnects for scalable integration of photonic neural networks

abstract: Photonic waveguides are prime candidates for integrated and parallel photonic interconnects. Such interconnects correspond to large-scale vector matrix products, which are at the heart of neural network computation. However, parallel interconnect circuits realized in two dimensions, for example by lithography, are strongly limited in size due to disadvantageous scaling. We use three dimensional (3D) printed photonic waveguides to overcome this limitation. 3D optical-couplers with fractal topology efficiently connect large numbers of input and output channels, and we show that the substrate's footprint area scales linearly. Going beyond simple couplers, we introduce functional circuits for discrete spatial filters identical to those used in deep convolutional neural networks.


title: Atomic responses to general dark matter-electron interactions

abstract: In the leading paradigm of modern cosmology, about 80% of our Universe's matter content is in the form of hypothetical, as yet undetected particles. These do not emit or absorb radiation at any observable wavelengths, and therefore constitute the so-called Dark Matter (DM) component of the Universe. Detecting the particles forming the Milky Way DM component is one of the main challenges for astroparticle physics and basic science in general. One promising way to achieve this goal is to search for rare DM-electron interactions in low-background deep underground detectors. Key to the interpretation of this search is the response of detectors' materials to elementary DM-electron interactions defined in terms of electron wave functions' overlap integrals. In this work, we compute the response of atomic argon and xenon targets used in operating DM search experiments to general, so far unexplored DM-electron interactions. We find that the rate at which atoms can be ionized via DM-electron scattering can in general be expressed in terms of four independent atomic responses, three of which we identify here for the first time. We find our new atomic responses to be numerically important in a variety of cases, which we identify and investigate thoroughly using effective theory methods. We then use our atomic responses to set 90% confidence level (C.L.) exclusion limits on the strength of a wide range of DM-electron interactions from the null result of DM search experiments using argon and xenon targets.


title: Automated label flows for excited states of correlation functions in lattice gauge theory

abstract: Extracting excited states from lattice gauge theory correlation functions can be achieved through chi-squared minimization fits or algebraic approaches such as the variational method and Prony's method. Performing any kind of error analysis, such as bootstrap resampling, often leads to overlapping confidence regions of model parameters, even when the spectrum is not particularly dense. In order to correctly estimate errors, one must beware of mislabeling the states. In this work, we provide an algorithm that we call automated label flows which consistently and systematically identifies a deterministic labeling of states. In the context of Prony's method, we analyze lattice correlation functions by using automated label flows, and compare the results to fits obtained from chi-square minimization fits to exponentials.