title: Riemann surfaces of genus $1+q^2$ with $3q^2$ automorphisms

abstract: In this article we classify compact Riemann surfaces of genus $1+q^2$ with a group of automorphisms of order $3q^2,$ where $q$ is a prime number. We also study decompositions of the corresponding Jacobian varieties.


title: First passage leapovers of L\'evy flights and the proper formulation of absorbing boundary conditions

abstract: An important open problem in the theory of L\'evy flights concerns the analytically tractable formulation of absorbing boundary conditions. Although numerical studies using the correctly defined nonlocal approach have yielded substantial insights regarding the statistics of first passage, the resultant modifications to the dynamical equations hinder the detailed analysis possible in the absence of these conditions. In this study it is demonstrated that using the first-hit distribution, related to the first passage leapover, as the absorbing sink preserves the tractability of the dynamical equations for a particle undergoing L\'evy flight. In particular, knowledge of the first-hit distribution is sufficient to fully determine the first passage time and position density of the particle, without requiring integral truncation or numerical simulations. In addition, we report on the first-hit and leapover properties of first passages and arrivals for L\'evy flights of arbitrary skew parameter, and extend these results to L\'evy flights in a certain ubiquitous class of potentials satisfying an integral condition.


title: Clustering of Local Group distances: publication bias or correlated measurements? VI. Extending to Virgo cluster distances

abstract: We have established an internally consistent Local Group distance framework, using the Galactic Center, the Large Magellanic Cloud, and Messier 31 (M31) as important stepping stones. At greater distances, few distance benchmarks are available. As a consequence, M87 and/or Virgo cluster distances are often invoked as the next rung on the ladder to more distant objects such as the Fornax and Coma clusters. Therefore, we extensively mined the published literature for independently derived distance estimates to either M87 or the center of the Virgo cluster. Based on our newly compiled, comprehensive database of 213 such distances, published between 1929 and 2017 July, we recommend an outward extension to our distance framework, $(m-M)_0^{\rm M87} = 31.03 \pm 0.14$ mag ($D = 16.07 \pm 1.03$ Mpc; where the uncertainty represents the Gaussian $\sigma$ of the distribution), based on a subset of recent (post-1990) M87/Virgo cluster distance measurements. The most stable distance tracers employed here were derived from analysis of both primary and secondary distance indicators. Among the former, we preferentially rely on Cepheid period--luminosity relations and red-giant-branch terminal magnitudes; our preferred secondary distance tracers are surface brightness fluctuations. Our updated distance modulus to M87 implies a slightly reduced black hole mass of $(5.9 \pm 0.6) \times 10^9 M_\odot$ with respect to that determined by the Event Horizon Telescope collaboration.


title: Mean-field calculations of charge radii in ground and isomeric states of Cd isotopes

abstract: Quadrupole moments and charge radii in cadmium isotopes are studied from a microscopic perspective. The results obtained from self-consistent deformed Hartree-Fock+BCS calculations with Skyrme forces are compared with isomer and isotope shifts measured from high resolution laser spectroscopy experiments. The microscopic calculations reproduce fairly well the main features observed in those isotopes that include the linear increase with the neutron number of the quadrupole moments of the 11/2- isomers, as well as the parabolic behavior of the mean-square charge radii difference between isomers and ground states.


title: Low power, fast and broadband ESR quantum control using a stripline resonator

abstract: Using a home-built Ku band ESR spectrometer equipped with an arbitrary waveform generator and a stripline resonator, we implement two types of pulses that would benefit quantum computers: BB1 composite pulse and a microwave frequency comb. Broadband type 1 (BB1) composite pulse is commonly used to combat systematic errors but previous experiments were carried out only on extremely narrow linewidth samples. Using a sample with a linewidth of 9.35 MHz, we demonstrate that BB1 composite pulse is still effective against pulse length errors at a Rabi frequency of 38.46 MHz. The fast control is realized with low microwave power which is required for initialization of electron spin qubits at 0.6 T. We also digitally design and implement a microwave frequency comb to excite multiple spin packets of a different sample. Using this pulse, we demonstrate coherent and well resolved excitations spanning over the entire spectrum of the sample (ranging from -20 to 20 MHz). In anticipation of scaling up to a system with large number of qubits, this approach provides an efficient technique to selectively and simultaneously control multiple qubits defined in the frequency-domain.


title: Incompressible limit of the Ericksen-Leslie hyperbolic liquid crystal model in compressible flow

abstract: Ericksen and Leslie proposed a hydrodynamic model for liquid crystals in the format of conservation laws in the 1960s. Their original model includes inertial and compressibility effects, which makes the model a coupled parabolic-hyperbolic system. In this paper we build up the connection between the compressible and incompressible parabolic-hyperbolic liquid crystal model in the framework of classical solutions. We first derive the scaled Ericksen-Leslie system with dimensionless numbers, including Mach, Reynolds, and Ericksen numbers. In particular, we introduce the so-called inertial constant $\chi$ which characterizes the inertial effect of the liquid crystal molecular. Next, we establish the energy estimates uniform in the Mach number $\epsilon$ for both the compressible system and its time-derivative system with small data. Then, we pass to the limit $\epsilon \rightarrow 0$ in the compressible system, so that we establish the global classical solution of the incompressible system by the compactness arguments. Moreover, we also obtain the convergence rates associated with $L^2$-norm in the case of well-prepared initial data. This is the first result on the incompressible limit of the compressible parabolic-hyperbolic liquid crystal model and confirms the relations of different parabolic-hyperbolic liquid crystal model rigorously.


title: Exoskeleton Knee Compliance Improves Gait Velocity and Stability in a Spinal Cord Injured User: A Case Report

abstract: Spinal cord injuries frequently impair the ability to walk. Powered lower limb exoskeletons offer a promising solution to restore walking ability. However, they are currently restricted to even ground. We hypothesized that compliant exoskeleton knees could decrease required effort to maneuver on uneven terrain, and increase gait velocity and stability. We describe a case study of a motor-complete spinal cord injury user (AIS A, Th12) walking with a powered exoskeleton on even and uneven ground over multiple sessions after extensive training. Measurements with compliant or rigid exoskeleton knee joints were performed on three different days for each configuration. Body motion and crutch ground interaction forces were recorded to assess gait performance. We observed higher walking speeds with a compliant exoskeleton knee configuration (mean: 0.116 m/s on uneven and 0.145 m/s on even ground) compared to a rigid configuration (mean: 0.083 m/s and 0.100 m/s). Crutch force impulse was significantly reduced in the compliant configuration. Lastly, gait was more symmetric when the knee joints were compliant. In conclusion, compliant exoskeleton knee joints can help maneuver uneven ground faster and with less user effort than rigid joints. Based on our findings, exoskeleton designers should consider introducing compliance into their design to increase gait robustness and performance, and render exoskeletons more suitable for daily life use.


title: Machine Learning for high speed channel optimization

abstract: Design of printed circuit board (PCB) stack-up requires the consideration of characteristic impedance, insertion loss and crosstalk. As there are many parameters in a PCB stack-up design, the optimization of these parameters needs to be efficient and accurate. A less optimal stack-up would lead to expensive PCB material choices in high speed designs. In this paper, an efficient global optimization method using parallel and intelligent Bayesian optimization is proposed for the stripline design.


title: The drifting subpulses of PSR B0031-07 and its synchronously modulated radio polarization

abstract: We establish that for PSR B0031-07 the orthogonal polarization modes switch at a single pulse level synchronously with the periodic drifting subpulses seen in total intensity. There are only four other pulsars known for which this phenomenon is observed. PSR B0031-07 is unique as it is the only source in this group which has multiple stable drift modes. For both drift modes visible at our observing frequency centered at 1369 MHz, the modulation of polarization modes is synchronous with the drifting subpulses. In one of the drift modes, a discontinuity in the modulation pattern of polarization properties occurs halfway through the pulse, coinciding with a slight change in the slope of the intensity drift band. In contrast to what has been suggested for this pulsar in the past, this, plus other differences in the polarization of the modulated emission observed for the two drift modes, suggests that a drift mode change is more than a change in the underlying carousel radius and that magnetospheric propagation effects play an important role. The ellipticity evolves asymmetrically in time during the modulation cycle, which in the framework of a carousel model implies that the polarized sub-beams are asymmetric with respect to the sense of circulation, something which is not observed for other pulsars. Birefringence in the magnetosphere, resulting in the orthogonal polarization modes to spatially separate, is not enough to explain these results. It is argued that more complex magnetospheric processes, which possibly allow conversion between orthogonal polarization modes, play a role.


title: Coupled Magnetic and Superconducting Transitions in Sr$_{2}$VO$_{3}$FeAs Under Pressure

abstract: We report muon spin rotation ($\mu$SR) and magnetization measurements on superconducting Sr$_{2}$VO$_{3}$FeAs under pressure. At ambient pressure, Sr$_{2}$VO$_{3}$FeAs undergoes an antiferromagnetic transition of the V moments at $T_{N}$ and becomes superconducting at $T_{c}<T_{N}$. As a function of pressure, $T_{N}$ initially decreases while $T_{c}$ increases. Surprisingly, once $T_{N}\approx T_{c}$ at 0.6~GPa, $T_{N}$ reverses its trend and increases together with $T_{c}$ suggesting that the static V magnetism is a prerequisite for superconductivity. We explain this cooperative coupling by a possible localization of the V 3$d$ states below the magnetic transition which enables the nesting of the Fermi surface necessary for superconductivity.


title: Identifying Strong Lenses with Unsupervised Machine Learning using Convolutional Autoencoder

abstract: In this paper we develop a new unsupervised machine learning technique comprised of a feature extractor, a convolutional autoencoder (CAE), and a clustering algorithm consisting of a Bayesian Gaussian mixture model (BGM). We apply this technique to visual band space-based simulated imaging data from the Euclid Space Telescope using data from the Strong Gravitational Lenses Finding Challenge. Our technique promisingly captures a variety of lensing features such as Einstein rings with different radii, distorted arc structures, etc, without using predefined labels. After the clustering process, we obtain several classification clusters separated by different visual features which are seen in the images. Our method successfully picks up $\sim$63\ percent of lensing images from all lenses in the training set. With the assumed probability proposed in this study, this technique reaches an accuracy of $77.25\pm 0.48$\% in binary classification using the training set. Additionally, our unsupervised clustering process can be used as the preliminary classification for future surveys of lenses to efficiently select targets and to speed up the labelling process. As the starting point of the astronomical application using this technique, we not only explore the application to gravitationally lensed systems, but also discuss the limitations and potential future uses of this technique.


title: Sobolev spaces in extended metric-measure spaces

abstract: These lecture notes contain an extended version of the material presented in the C.I.M.E. summer course in 2017, aiming to give a detailed introduction to the metric Sobolev theory. The notes are divided in four main parts. The first one is devoted to a preliminary study of the underlying topological, metric, and measure-theoretic aspects of a general extended metric-topological measure space $\mathbb X=(X,\tau,\mathsf d,\mathfrak m)$. The second part is devoted to the construction of the Cheeger energy, initially defined on a distinguished unital algebra of Lipschitz functions. The third part deals with the basic tools needed for the dual characterization of the Sobolev spaces: the notion of $p$-Modulus of a collection of (nonparametric) rectifiable arcs and its duality with the class of nonparametric dynamic plans, i.e.~Radon measures on the space of rectifiable arcs with finite $q$-barycentric entropy with respect to $\mathfrak m$. The final part of the notes is devoted to the dual/weak formulation of the Sobolev spaces $W^{1,p}(\mathbb X)$ in terms of nonparametric dynamic plans and to their relations with the Newtonian spaces $N^{1,p}(\mathbb X)$ and with the spaces $H^{1,p}(\mathbb X)$ obtained by the Cheeger construction. In particular, when $(X,\mathsf d)$ is complete, a new proof of the equivalence between these approaches is given by a direct duality argument. A substantial part of these Lecture notes relies on well established theories. New contributions concern the extended metric setting, the role of general compatible algebras of Lipschitz functions and their density w.r.t.~the Sobolev energy, a compactification trick, the study of reflexivity and infinitesimal Hilbertianity inherited from the underlying space, and the use of nonparametric dynamic plans for the definition of weak upper gradients.


title: Kernel Dependence Regularizers and Gaussian Processes with Applications to Algorithmic Fairness

abstract: Current adoption of machine learning in industrial, societal and economical activities has raised concerns about the fairness, equity and ethics of automated decisions. Predictive models are often developed using biased datasets and thus retain or even exacerbate biases in their decisions and recommendations. Removing the sensitive covariates, such as gender or race, is insufficient to remedy this issue since the biases may be retained due to other related covariates. We present a regularization approach to this problem that trades off predictive accuracy of the learned models (with respect to biased labels) for the fairness in terms of statistical parity, i.e. independence of the decisions from the sensitive covariates. In particular, we consider a general framework of regularized empirical risk minimization over reproducing kernel Hilbert spaces and impose an additional regularizer of dependence between predictors and sensitive covariates using kernel-based measures of dependence, namely the Hilbert-Schmidt Independence Criterion (HSIC) and its normalized version. This approach leads to a closed-form solution in the case of squared loss, i.e. ridge regression. Moreover, we show that the dependence regularizer has an interpretation as modifying the corresponding Gaussian process (GP) prior. As a consequence, a GP model with a prior that encourages fairness to sensitive variables can be derived, allowing principled hyperparameter selection and studying of the relative relevance of covariates under fairness constraints. Experimental results in synthetic examples and in real problems of income and crime prediction illustrate the potential of the approach to improve fairness of automated decisions.


title: Control of photodissociation with dynamic Stark effect induced by Thz pulses

abstract: We demonstrate how dynamic Stark control (DSC) can be achieved on molecular photodissociation in the dipole limit, using single-cycle (FWHM) laser pulses in the terahertz (THz) regime. As the laser-molecule interaction follows the instantaneous electric field through the permanent dipoles, the molecular potentials dynamically oscillate and so does the crossings between them. In this paper, we consider rotating-vibrating diatomic molecules (2D description) and reveal the interplay between the dissociating wave packet and the dynamically fluctuating crossing seam located in the configuration space of the molecules spanned by the R vibrational and $\theta$ rotational coordinates. Our showcase example is the widely studied lithium-fluoride (LiF) molecule for which the two lowest $\Sigma$ states are nonadiabatically coupled at an avoided crossing (AC), furthermore a low-lying pure repulsive $\Pi$ state is energetically close. Optical pumping of the system in the ground state thus results in two dissociation channels: one indirect route via the AC in the ground $\Sigma$ state and one direct path in the $\Pi$ state. We show that applying THz control pulses with specific time delays relative to the pumping, can significantly alter the population dynamics, as well as, the kinetic energy and angular distribution of the photofragments.


title: Scale-up of room-temperature constructive quantum interference from single molecules to self-assembled molecular-electronic films

abstract: The realization of self-assembled molecular-electronic films, whose room-temperature transport properties are controlled by quantum interference (QI), is an essential step in the scale-up QI effects from single molecules to parallel arrays of molecules. Recently, the effect of destructive QI (DQI) on the electrical conductance of self-assembled monolayers (SAMs) has been investigated. Here, through a combined experimental and theoretical investigation, we demonstrate chemical control of different forms of constructive QI (CQI) in cross-plane transport through SAMs and assess its influence on cross-plane thermoelectricity in SAMs. It is known that the electrical conductance of single molecules can be controlled in a deterministic manner, by chemically varying their connectivity to external electrodes. Here, by employing synthetic methodologies to vary the connectivity of terminal anchor groups around aromatic anthracene cores, and by forming SAMs of the resulting molecules, we clearly demonstrate that this signature of CQI can be translated into SAM-on-gold molecular films. We show that the conductance of vertical molecular junctions formed from anthracene-based molecules with two different connectivities differ by a factor of approximately 16, in agreement with theoretical predictions for their conductance ratio based on constructive QI effects within the core. We also demonstrate that for molecules with thiol anchor groups, the Seebeck coefficient of such films is connectivity dependent and with an appropriate choice of connectivity can be boosted by ~50%. This demonstration of QI and its influence on thermoelectricity in SAMs represents a critical step towards functional ultra-thin-film devices for future thermoelectric and molecular-scale electronics applications.


title: Muon spin rotation measurements on RbEuFe$_{4}$As$_{4}$ under pressure

abstract: We report muon spin rotation and magnetization measurements on the magnetic superconductor RbEuFe$_{4}$As$_{4}$ under hydrostatic pressures up to 3.8 GPa. At ambient pressure, RbEuFe$_{4}$As$_{4}$ exhibits a superconducting transition at $T_{c} \approx$ 36.5 K and a magnetic transition at $T_{m} \approx$ 15 K below which the magnetic and the superconducting order coexist. With increasing pressure, $T_{c}$ decreases while $T_{m}$ and the ordered Eu magnetic moment increase. In contrast to iron-based superconductors with ordering Fe moments, the size of the ordered Eu moment is not proportional to $T_{m}$. The muon spin rotation signal is dominated by the magnetic response impeding the determination of the superconducting properties.


title: ASP-Core-2 Input Language Format

abstract: Standardization of solver input languages has been a main driver for the growth of several areas within knowledge representation and reasoning, fostering the exploitation in actual applications. In this document we present the ASP-Core-2 standard input language for Answer Set Programming, which has been adopted in ASP Competition events since 2013.


title: Far-Field Asymptotics for Multiple-Pole Solitons in the Large-Order Limit

abstract: The integrable focusing nonlinear Schrodinger equation admits soliton solutions whose associated spectral data consist of a single pair of conjugate poles of arbitrary order. We study families of such multiple-pole solitons generated by Darboux transformations as the pole order tends to infinity. We show that in an appropriate scaling, there are four regions in the space-time plane where solutions display qualitatively distinct behaviors: an exponential-decay region, an algebraic-decay region, a non-oscillatory region, and an oscillatory region. Using the nonlinear steepest-descent method for analyzing Riemann-Hilbert problems, we compute the leading-order asymptotic behavior in the algebraic-decay, non-oscillatory, and oscillatory regions.


title: Eliashberg equations for an electron-phonon version of the Sachdev-Ye-Kitaev model: Pair Breaking in non-Fermi liquid superconductors

abstract: We present a theory that is a non-Fermi-liquid counterpart of the Abrikosov-Gor'kov pair-breaking theory due to paramagnetic impurities in superconductors. To this end we analyze a model of interacting electrons and phonons that is a natural generalization of the Sachdev-Ye-Kitaev-model. In the limit of large numbers of degrees of freedom, the Eliashberg equations of superconductivity become exact and emerge as saddle-point equations of a field theory with fluctuating pairing fields. In its normal state the model is governed by two non-Fermi liquid fixed points, characterized by distinct universal exponents. At low temperatures a superconducting state emerges from the critical normal state. We study the role of pair-breaking on $T_{c}$, where we allow for disorder that breaks time-reversal symmetry. For small Bogoliubov quasi-particle weight, relevant for systems with strongly incoherent normal state, $T_{c}$ drops rapidly as function of the pair breaking strength and reaches a small but finite value before it vanishes at a critical pair-breaking strength via an essential singularity. The latter signals a breakdown of the emergent conformal symmetry of the non-Fermi liquid normal state.


title: Tumor stabilization induced by T-cell recruitment fluctuations

abstract: The influence of random fluctuations on the recruitment of effector cells towards a tumor is studied by means of a stochastic mathematical model. Aggressively growing tumors are confronted against varying intensities of the cell-mediated immune response for which chaotic and periodic oscillations coexist together with stable tumor dynamics. A thorough parametric analysis of the noise-induced transition from this oscillatory regime to complete tumor dominance is carried out. A hysteresis phenomenon is uncovered, which stabilizes the tumor at its carrying capacity and drives the healthy and the immune cell populations to their extinction. Furthermore, it is shown that near a crisis bifurcation such transitions occur under weak noise intensities. Finally, the corresponding noise-induced chaos-order transformation is analyzed and discussed in detail.


title: Security Issues in WBANs

abstract: Wireless Body Area Network (WBAN) refers to short-range, wireless communications near or inside a human body. WBAN is emerging solution to cater the needs of local and remote health care related facility. Medical and non-medical applications have been revolutionarily under consideration for providing a healthy and gratify service to the humanity. Being very critical in communication from body it faces a lot of challenges which are to be tackled for the safety of life and benefit of the user. There is variety of challenges faced by WBANs. WBAN is favorite playground for attackers due to its usability in various applications. This article provides systematic overview of main challenges in WBANs in security perspectives.


title: Poisson Cohomology of holomorphic toric Poisson manifolds. II

abstract: In this paper, we give a description of holomorphic multi-vector fields on smooth compact toric varieties, which generalizes Demazure's result of holomorphic vector fields on toric varieties. Based on the result, we compute the Poisson cohomology groups of holomorphic toric Poisson manifolds, i.e., toric varieties endowed with $T$-invariant holomorphic Poisson structures.


title: Compound-tunable embedding potential method and its application to ytterbium fluoride crystals YbF$_2$ and YbF$_3$

abstract: Compound-tunable embedding potential (CTEP) method developed in previous works to describe electronic structure of fragments in materials is applied to crystals containing periodically arranged lanthanide atoms, which can have open $4f$ shell. We consider YbF$_2$ and YbF$_3$ as examples such that $4f$ shell is excluded from both the crystal and cluster stages of generating the CTEP. Instead, 10 and 11 valence-electron pseudopotentials for Yb, correspondingly, are applied and the latter treats the $4f$-hole implicitly. At the next stage of the two-component embedded cluster studies of the YbF$_{2,3}$ crystals we apply the 42 valence-electron relativistic pseudopotential for Yb and, thus, $4f$ shell is treated explicitly. A remarkable agreement of the electronic density and interatomic distances within the fragment with those of the original periodic crystal calculation is attained.


title: Photoexcited organic molecules $en~route$ to highly efficient autoionization

abstract: The conversion of optical and electrical energy in novel materials is key to modern optoelectronic and light-harvesting applications. Here, we investigate the equilibration dynamics of photoexcited 2,7-bis(biphenyl-4-yl)-2,7-ditertbutyl-9,9-spirobiuorene (SP6) molecules adsorbed on ZnO(10-10) using femtosecond time-resolved two-photon photoelectron (2PPE) and optical spectroscopy. We find that, after initial ultrafast relaxation on fs and ps timescales, an optically dark state is populated, likely the SP6 triplet (T) state, that undergoes Dexter-type energy transfer ($r_{\mathrm{Dex}} = 1.3~\mathrm{nm}$) and exhibits a long decay time of 0.1 s. Because of this long lifetime a photostationary state with average T-T distances below 2 nm is established at excitation densities in the $10^{20}~\mathrm{cm}^{-2}~\mathrm{s}^{-1}$ range. This large density enables decay by T-T annihilation (TTA) mediating autoionization despite an extremely low TTA rate of $k_{\mathrm{TTA}} = 4.5~10^{-26}~\mathrm{m}^3~\mathrm{s}^{-1}$. The large external quantum efficiency of the autoionization process (up to 15 %) and photocurrent densities in the \mathrm{mA}~\mathrm{cm}^{-2}$ range offer great potential for light-harvesting applications.


title: Stochastic Difference-of-Convex Algorithms for Solving nonconvex optimization problems

abstract: The paper deals with stochastic difference-of-convex functions (DC) programs, that is, optimization problems whose the cost function is a sum of a lower semicontinuous DC function and the expectation of a stochastic DC function with respect to a probability distribution. This class of nonsmooth and nonconvex stochastic optimization problems plays a central role in many practical applications. Although there are many contributions in the context of convex and/or smooth stochastic optimization, algorithms dealing with nonconvex and nonsmooth programs remain rare. In deterministic optimization literature, the DC Algorithm (DCA) is recognized to be one of the few algorithms to solve effectively nonconvex and nonsmooth optimization problems. The main purpose of this paper is to present some new stochastic DCAs for solving stochastic DC programs. The convergence analysis of the proposed algorithms is carefully studied, and numerical experiments are conducted to justify the algorithms' behaviors.


title: Systematic Comparison of the Influence of Different Data Preprocessing Methods on the Performance of Gait Classifications Using Machine Learning

abstract: Human movements are characterized by highly non-linear and multi-dimensional interactions within the motor system. Recently, an increasing emphasis on machine-learning applications has led to a significant contribution to the field of gait analysis, e.g., in increasing the classification performance. In order to ensure the generalizability of the machine-learning models, different data preprocessing steps are usually carried out to process the measured raw data before the classifications. In the past, various methods have been used for each of these preprocessing steps. However, there are hardly any standard procedures or rather systematic comparisons of these different methods and their impact on the classification performance. Therefore, the aim of this analysis is to compare different combinations of commonly applied data preprocessing steps and test their effects on the classification performance of gait patterns. A publicly available dataset on intra-individual changes of gait patterns was used for this analysis. Forty-two healthy participants performed 6 sessions of 15 gait trials for 1 day. For each trial, two force plates recorded the 3D ground reaction forces (GRFs). The data was preprocessed with the following steps: GRF filtering, time derivative, time normalization, data reduction, weight normalization and data scaling. Subsequently, combinations of all methods from each preprocessing step were analyzed by comparing their prediction performance in a six-session classification using Support Vector Machines, Random Forest Classifiers, Multi-Layer Perceptrons, and Convolutional Neural Networks. In conclusion, the present results provide first domain-specific recommendations for commonly applied data preprocessing methods and might help to build more comparable and more robust classification models based on machine learning that are suitable for a practical application.


title: Fair Meta-Learning: Learning How to Learn Fairly

abstract: Data sets for fairness relevant tasks can lack examples or be biased according to a specific label in a sensitive attribute. We demonstrate the usefulness of weight based meta-learning approaches in such situations. For models that can be trained through gradient descent, we demonstrate that there are some parameter configurations that allow models to be optimized from a few number of gradient steps and with minimal data which are both fair and accurate. To learn such weight sets, we adapt the popular MAML algorithm to Fair-MAML by the inclusion of a fairness regularization term. In practice, Fair-MAML allows practitioners to train fair machine learning models from only a few examples when data from related tasks is available. We empirically exhibit the value of this technique by comparing to relevant baselines.


title: Bayesian Non-Parametric Factor Analysis for Longitudinal Spatial Surfaces

abstract: We introduce a Bayesian non-parametric spatial factor analysis model with spatial dependency induced through a prior on factor loadings. For each column of the loadings matrix, spatial dependency is encoded using a probit stick-breaking process (PSBP) and a multiplicative gamma process shrinkage prior is used across columns to adaptively determine the number of latent factors. By encoding spatial information into the loadings matrix, meaningful factors are learned that respect the observed neighborhood dependencies, making them useful for assessing rates over space. Furthermore, the spatial PSBP prior can be used for clustering temporal trends, allowing users to identify regions within the spatial domain with similar temporal trajectories, an important task in many applied settings. In the manuscript, we illustrate the model's performance in simulated data, but also in two real-world examples: longitudinal monitoring of glaucoma and malaria surveillance across the Peruvian Amazon. The R package spBFA, available on CRAN, implements the method.


title: Active Learning for Black-Box Adversarial Attacks in EEG-Based Brain-Computer Interfaces

abstract: Deep learning has made significant breakthroughs in many fields, including electroencephalogram (EEG) based brain-computer interfaces (BCIs). However, deep learning models are vulnerable to adversarial attacks, in which deliberately designed small perturbations are added to the benign input samples to fool the deep learning model and degrade its performance. This paper considers transferability-based black-box attacks, where the attacker trains a substitute model to approximate the target model, and then generates adversarial examples from the substitute model to attack the target model. Learning a good substitute model is critical to the success of these attacks, but it requires a large number of queries to the target model. We propose a novel framework which uses query synthesis based active learning to improve the query efficiency in training the substitute model. Experiments on three convolutional neural network (CNN) classifiers and three EEG datasets demonstrated that our method can improve the attack success rate with the same number of queries, or, in other words, our method requires fewer queries to achieve a desired attack performance. To our knowledge, this is the first work that integrates active learning and adversarial attacks for EEG-based BCIs.


title: Slingshot prominences: coronal structure, mass loss and spin down

abstract: The structure of a star's coronal magnetic field is a fundamental property that governs the high-energy emission from the hot coronal gas and the loss of mass and angular momentum in the stellar wind. It is, however, extremely difficult to measure. We report a new method to trace this structure in rapidly-rotating young convective stars, using the cool gas trapped on coronal field lines as markers. This gas forms "slingshot prominences" which appear as transient absorption features in H-$\alpha$. By using different methods of extrapolating this field from the surface measurements, we determine locations for prominence support and produce synthetic H-$\alpha$ stacked spectra. The absorption features produced with a potential field extrapolation match well this those observed, while those from a non-potential field do not. In systems where the rotation and magnetic axes are well aligned, up to $50\%$ of the prominence mass may transit the star and so produces a observable feature. This fraction may fall as low as $~2\%$ in very highly inclined systems. Ejected prominences carry away mass and angular momentum at rates that vary by two orders of magnitude, but which may approach those carried by the stellar wind.


title: The systemic growth constants of climate change: From its origin in 1780 to its major post-WWII acceleration

abstract: The recent discovery of long term growth constants in the accumulation of atmospheric CO2, confirmed by two methods, enables analog methods for dating the beginning of climate change at ~1780 and projecting its near term future. Here we show that the preceding wavy variation in CO2 PPM abruptly shifts to exponential, moving symmetrically around a growth constant of 1.48 %/yr until WWII, and after a pause rises to hover about a higher constant growth constant from 1960 on of 2.0 %/yr. Such long term steady states of global environmental change suggest transitions between stable states of global self-organization. A method of analog curve fitting to project the current steady state of acceleration is tested, suggesting a 2040 earth temperature rise of 1.89 deg C above the IPCC baseline. A very brief following discussion of what a systemic growth constant for atmospheric CO2 implies and model strategies for responding to it.


title: A Minimal Contrast Estimator for the Linear Fractional Stable Motion

abstract: In this paper we present an estimator for the three-dimensional parameter $(\sigma, \alpha, H)$ of the linear fractional stable motion, where $H$ represents the self-similarity parameter, and $(\sigma, \alpha)$ are the scaling and stability parameters of the driving symmetric L\'evy process $L$. Our approach is based upon a minimal contrast method associated with the empirical characteristic function combined with a ratio type estimator for the selfsimilarity parameter $H$. The main result investigates the strong consistency and weak limit theorems for the resulting estimator. Furthermore, we propose several ideas to obtain feasible confidence regions in various parameter settings. Our work is mainly related to [16, 18], in which parameter estimation for the linear fractional stable motion and related L\'evy moving average processes has been studied.


title: Rydberg-positronium velocity and self-ionization studies in 1T magnetic field and cryogenic environment

abstract: We characterized the pulsed Rydberg-positronium production inside the AEgIS (Antimatter Experiment: Gravity, Interferometry, Spectroscopy) apparatus in view of antihydrogen formation by means of a charge exchange reaction between cold antiprotons and slow Rydberg-positronium atoms. Velocity measurements on positronium along two axes in a cryogenic environment (10K) and in 1T magnetic field were performed. The velocimetry was done by MCP-imaging of photoionized positronium previously excited to the $n=3$ state. One direction of velocity was measured via Doppler-scan of this $n=3$-line, another direction perpendicular to the former by delaying the exciting laser pulses in a time-of-flight measurement. Self-ionization in the magnetic field due to motional Stark effect was also quantified by using the same MCP-imaging technique for Rydberg positronium with an effective principal quantum number $n_{eff}$ ranging between 14 and 22. We conclude with a discussion about the optimization of our experimental parameters for creating Rydberg-positronium in preparation for an efficient pulsed production of antihydrogen.


title: On the Existence of Heterotic-String and Type-II-Superstring Field Theory Vertices

abstract: We consider two problems associated to the space of Riemann surfaces that are closely related to string theory. We first consider the problem of existence of heterotic-string and type-II-superstring field theory vertices in the product of spaces of bordered surfaces parameterizing left- and right-moving sectors of these theories. It turns out that this problem can be solved by proving the existence of a solution to the BV quantum master equation in moduli spaces of bordered spin-Riemann surfaces. We first prove that for arbitrary genus $g$, $n_{\text{NS}}$ Neveu-Schwarz boundary components, and $n_{\text{R}}$ Ramond boundary components such solutions exist. We also prove that these solutions are unique up to homotopies in the category of BV algebras and are mapped to fundamental classes of Deligne-Mumford stacks of associated punctured spin curves. These results generalize the work of Costello on the existence of a solution to the BV quantum master equations in moduli spaces of bordered ordinary Riemann surfaces. Using this, we prove that heterotic-string and type-II-superstring field theory vertices exist. Furthermore, we prove the existence of a solution to the BV quantum master equation in spaces of bordered $\mathscr{N}=1$ super-Riemann surfaces for arbitrary genus $g$, $n_{\text{NS}}$ Neveu-Schwarz boundary components, and $n_{\text{R}}$ Ramond boundary components. Secondly, we turn to the moduli problem of stable $\mathscr{N}=1$ SUSY curves. Using Vaintrob's deformation theory of complex superspaces and Deligne's model for the compactification of the space of $\mathscr{N}=1$ SUSY curves, we prove the representability of this moduli problem by a proper smooth Deligne-Mumford superstack. This is a generalization of the work of Deligne and Mumford on the representability of the moduli problem of stable ordinary curves by a Deligne-Mumford stack.


title: Bunch theory: working notes on applications, axioms and models

abstract: In his book A Practical Theory of Programming, Eric Hehner proposes and applies a remarkably radical reformulation of set theory, in which the collection and packaging of elements are seen as separate activities. This provides for unpackaged collections, referred to as bunches. Bunches allow us to reason about non-determinism at the level of terms, and, very remarkably, allow us to reason about the conceptual entity nothing, which is just an empty bunch (and very different from an empty set). This eliminates mathematical gaps caused by undefined terms. We compare the use of bunches with other approaches to this problem, and we illustrate the use of Bunch Theory in formulating program semantics combining non-deterministic, preferential, and probabilistic choice to provide a guarded command language whose exceptional expressivity we illustrate with a short case study. We show how an existing axiomatisation of set theory can be extended to incorporate bunches, and we provide and validate a model.


title: Thom spectra, higher $THH$ and tensors in $\infty$-categories

abstract: Let $f:G\to \mathrm{Pic}(R)$ be a map of $E_\infty$-groups, where $\mathrm{Pic}(R)$ denotes the Picard space of an $E_\infty$-ring spectrum $R$. We determine the tensor $X\otimes_R Mf$ of the Thom $E_\infty$-$R$-algebra $Mf$ with a space $X$; when $X$ is the circle, the tensor with $X$ is topological Hochschild homology over $R$. We use the theory of localizations of $\infty$-categories as a technical tool: we contribute to this theory an $\infty$-categorical analogue of Day's reflection theorem about closed symmetric monoidal structures on localizations, and we prove that for a smashing localization $L$ of the $\infty$-category of presentable $\infty$-categories, the free $L$-local presentable $\infty$-category on a small simplicial set $K$ is given by presheaves on $K$ valued on the $L$-localization of the $\infty$-category of spaces. If $X$ is a pointed space, a map $g: A\to B$ of $E_\infty$-ring spectra satisfies $X$-base change if $X\otimes B$ is the pushout of $A\to X\otimes A$ along $g$. Building on a result of Mathew, we prove that if $g$ is \'etale then it satisfies $X$-base change provided $X$ is connected. We also prove that $g$ satisfies $X$-base change provided the multiplication map of $B$ is an equivalence. Finally, we prove that, under some hypotheses, the Thom isomorphism of Mahowald cannot be an instance of $S^0$-base change.


title: A note on Franel numbers and $SU(3)$

abstract: We study weight multiplicities in tensor powers of the adjoint representation of $SU(3)$ and relate them to Franel numbers.


title: A structure-preserving FEM for the uniaxially constrained $\mathbf{Q}$-tensor model of nematic liquid crystals

abstract: We consider the one-constant Landau - de Gennes model for nematic liquid crystals. The order parameter is a traceless tensor field $\mathbf{Q}$, which is constrained to be uniaxial: $\mathbf{Q} = s (\mathbf{n}\otimes\mathbf{n} - d^{-1} \mathbf{I})$ where $\mathbf{n}$ is a director field, $s\in\mathbb{R}$ is the degree of orientation, and $d\ge2$ is the dimension. Building on similarities with the one-constant Ericksen energy, we propose a structure-preserving finite element method for the computation of equilibrium configurations. We prove stability and consistency of the method without regularization, and $\Gamma$-convergence of the discrete energies towards the continuous one as the mesh size goes to zero. We design an alternating direction gradient flow algorithm for the solution of the discrete problems, and we show that such a scheme decreases the energy monotonically. Finally, we illustrate the method's capabilities by presenting some numerical simulations in two and three dimensions including non-orientable line fields.


title: Semi-discrete optimal transport

abstract: In the current book I suggest an off-road path to the subject of optimal transport. I tried to avoid prior knowledge of analysis, PDE theory and functional analysis, as much as possible. Thus I concentrate on discrete and semi-discrete cases, and always assume compactness for the underlying spaces. However, some fundamental knowledge of measure theory and convexity is unavoidable. In order to make it as self-contained as possible I included an appendix with some basic definitions and results. I believe that any graduate student in mathematics, as well as advanced undergraduate students, can read and understand this book. Some chapters (in particular in Parts II\&III ) can also be interesting for experts. Starting with the the most fundamental, fully discrete problem I attempted to place optimal transport as a particular case of the celebrated stable marriage problem. From there we proceed to the partition problem, which can be formulated as a transport from a continuous space to a discrete one. Applications to information theory and game theory (cooperative and non-cooperative) are introduced as well. Finally, the general case of transport between two compact measure spaces is introduced as a coupling between two semi-discrete transports.


title: Unconditional uniqueness of solutions for nonlinear dispersive equations

abstract: When a solution to the Cauchy problem for nonlinear dispersive equations is obtained by a fixed point argument using auxiliary function spaces, it is non-trivial to ensure uniqueness of solutions in a natural space such as the class of continuous curves in the data space. This property is called unconditional uniqueness, and proving it often requires some additional work. In the last decade, unconditional uniqueness has been shown for some canonical nonlinear dispersive equations by an integration-by-parts technique, which can be regarded as a variant of the (Poincar\'e-Dulac) normal form reduction. In this article, we aim to provide an abstract framework for establishing unconditional uniqueness as well as existence of certain weak solutions via infinite iteration of the normal form reduction. In particular, in an abstract setting we find two sets of fundamental estimates, each of which can be used repeatedly to generate all multilinear estimates of arbitrarily high degrees required in this scheme. Then, we confirm versatility of the framework by applying it to various equations, including the cubic nonlinear Schr\"odinger equation (NLS) in higher dimension, the cubic NLS with fractional Laplacians, the cubic derivative NLS, and the Zakharov system, for which new results on unconditional uniqueness are obtained under the periodic boundary condition.


title: Recurrence versus Transience for Weight-Dependent Random Connection Models

abstract: We investigate random graphs on the points of a Poisson process in $d$-dimensional space, which combine scale-free degree distributions and long-range effects. Every Poisson point carries an independent random mark and given marks and positions of the points we form an edge between two points independently with a probability depending via a kernel on the two marks and the distance of the points. Different kernels allow the mark to play different roles, like weight, radius or birth time of a vertex. The kernels depend on a parameter~$\gamma$, which determines the power-law exponent of the degree distributions. A further independent parameter $\delta$ characterises the decay of the connection probabilities of vertices as their distance increases. We prove transience of the infinite cluster in the entire supercritical phase in regimes given by the parameters $\gamma$ and~$\delta$, and complement these results by recurrence results if $d=2$. Our results are particularly interesting for the soft Boolean graph model discussed in the preprint [arXiv:2108:11252] and the age-dependent random connection model recently introduced by Gracar et al.\ [Queueing Syst. 93.3-4 (2019)]}


title: Stronger Convergence Results for Deep Residual Networks: Network Width Scales Linearly with Training Data Size

abstract: Deep neural networks are highly expressive machine learning models with the ability to interpolate arbitrary datasets. Deep nets are typically optimized via first-order methods and the optimization process crucially depends on the characteristics of the network as well as the dataset. This work sheds light on the relation between the network size and the properties of the dataset with an emphasis on deep residual networks (ResNets). Our contribution is that if the network Jacobian is full rank, gradient descent for the quadratic loss and smooth activation converges to the global minima even if the network width $m$ of the ResNet scales linearly with the sample size $n$, and independently from the network depth. To the best of our knowledge, this is the first work which provides a theoretical guarantee for the convergence of neural networks in the $m=\Omega(n)$ regime.


title: A novel class of stabilized greedy kernel approximation algorithms: Convergence, stability & uniform point distribution

abstract: Kernel based methods provide a way to reconstruct potentially high-dimensional functions from meshfree samples, i.e., sampling points and corresponding target values. A crucial ingredient for this to be successful is the distribution of the sampling points. Since the computation of an optimal selection of sampling points may be an infeasible task, one promising option is to use greedy methods. Although these methods may be very effective, depending on the specific greedy criterion the chosen points might quickly lead to instabilities in the computation. To circumvent this problem, we introduce and investigate a new class of \textit{stabilized} greedy kernel algorithms, which can be used to create a scale of new selection strategies. We analyze these algorithms, and in particular we prove convergence results and quantify in a precise way the distribution of the selected points. These results allow to prove, in the case of certain Sobolev kernels, that the algorithms have optimal stability and optimal convergence rates, including for functions outside the native space of the kernel. The results also apply to the case of the usual $P$-greedy algorithm, significantly improving state-of-the-art results available in the literature. Illustrative experiments are presented that support the theoretical findings.


title: On Quantifying Large Lattice Relaxations in Photovoltaic Devices

abstract: Temporal variations of Cu(In,Ga)Se$_2$ photovoltaic device properties during light exposure at various temperatures and voltage biases for times up to 100 h were analyzed using the kinetic theory of large lattice relaxations. Open-circuit voltage and p-type doping increased with charge injection and decreased with temperature at low injection conditions. Lattice relaxation can account for both trends and activation energies extracted from the data were approximately 0.9 and 1.2 eV for devices with lower and higher sodium content, respectively. In these devices, increased sodium content resulted in higher initial p-type doping with greater stability. First principles calculations providing revised activation energies for the ($V_{Se}-V_{Cu}$) complex suggest that this defect does not account for the metastability observed here.


title: Almost-Killing equation: Stability, hyperbolicity, and black hole Gauss law

abstract: We examine the Hamiltonian formulation and hyperbolicity of the almost-Killing equation (AKE). We find that for all but one parameter choice, the Hamiltonian is unbounded, and in some cases, the AKE has ghost degrees of freedom. We also show the AKE is only strongly hyperbolic for one parameter choice, which corresponds to a case in which the AKE has ghosts. Fortunately, one finds that the AKE reduces to the homogeneous Maxwell equation in a vacuum, so that with the addition of the divergence-free constraint (a "Lorenz gauge"), one can still obtain a well-posed problem that is stable in the sense that the corresponding Hamiltonian is positive definite. An analysis of the resulting Komar currents reveals an exact Gauss law for a system of black holes in vacuum spacetimes and suggests a possible measure of matter content in asymptotically flat spacetimes.


title: The Crisanti-Sommers Formula for Spherical Spin Glasses with Vector Spins

abstract: We obtain the analogue of the Crisanti-Sommers variational formula for spherical spin glasses with vector spins. This formula is derived from the discrete Parisi variational formula for the limit of the free energy of constrained copies of spherical spin glasses. In vector spin models, the variations of the functional order parameters must preserve the monotonicity of matrix paths which introduces a new challenge in contrast to the derivation of the classical Crisanti-Sommers formula.


title: Constraining the Local Burst Rate Density of Primordial Black Holes with HAWC

abstract: Primordial Black Holes (PBHs) may have been created by density fluctuations in the early Universe and could be as massive as $> 10^9$ solar masses or as small as the Planck mass. It has been postulated that a black hole has a temperature inversely-proportional to its mass and will thermally emit all species of fundamental particles via Hawking Radiation. PBHs with initial masses of $\sim 5 \times 10^{14}$ g (approximately one gigaton) should be expiring today with bursts of high-energy gamma radiation in the GeV--TeV energy range. The High Altitude Water Cherenkov (HAWC) Observatory is sensitive to gamma rays with energies of $\sim$300 GeV to past 100 TeV, which corresponds to the high end of the PBH burst spectrum. With its large instantaneous field-of-view of $\sim 2$ sr and a duty cycle over 95%, the HAWC Observatory is well suited to perform an all-sky search for PBH bursts. We conducted a search using 959 days of HAWC data and exclude the local PBH burst rate density above $3400~\mathrm{pc^{-3}~yr^{-1}}$ at 99% confidence, the strongest limit on the local PBH burst rate density from any existing electromagnetic measurement.


title: Limited View and Sparse Photoacoustic Tomography for Neuroimaging with Deep Learning

abstract: Photoacoustic tomography (PAT) is a nonionizing imaging modality capable of acquiring high contrast and resolution images of optical absorption at depths greater than traditional optical imaging techniques. Practical considerations with instrumentation and geometry limit the number of available acoustic sensors and their view of the imaging target, which result in significant image reconstruction artifacts degrading image quality. Iterative reconstruction methods can be used to reduce artifacts but are computationally expensive. In this work, we propose a novel deep learning approach termed pixelwise deep learning (PixelDL) that first employs pixelwise interpolation governed by the physics of photoacoustic wave propagation and then uses a convolution neural network to directly reconstruct an image. Simulated photoacoustic data from synthetic vasculature phantom and mouse-brain vasculature were used for training and testing, respectively. Results demonstrated that PixelDL achieved comparable performance to iterative methods and outperformed other CNN-based approaches for correcting artifacts. PixelDL is a computationally efficient approach that enables for realtime PAT rendering and for improved image quality, quantification, and interpretation.


title: Minimal colorings for properly colored subgraphs in complete graphs

abstract: Let $pr(K_{n}, G)$ be the maximum number of colors in an edge-coloring of $K_{n}$ with no properly colored copy of $G$. In this paper, we show that $pr(K_{n}, G)-ex(n, \mathcal{G'})=o(n^{2}), $ where $\mathcal{G'}=\{G-M: M \text{ is a matching of }G\}$. Furthermore, we determine the value of $pr(K_{n}, P_{l})$ for $l\ge 27$ and $n\ge 2l^{3}$ and the exact value of $pr(K_{n}, G)$, where $G$ is $C_{5}, C_{6}$ and $K_{4}^{-}$, respectively. Also, we give an upper bound and a lower bound of $pr(K_{n}, K_{2,3})$.


title: Relativistic and non-Gaussianity contributions to the one-loop power spectrum

abstract: We compute the one-loop density power spectrum including Newtonian and relativistic contributions, as well as the primordial non-Gaussianity contributions from $f_{\rm NL}$ and $g_{\rm NL}$ in the local configuration. To this end we take solutions to the Einstein equations in the long-wavelength approximation and provide expressions for the matter density perturbation at second and third order. These solutions have shown to be complementary to the usual Newtonian cosmological perturbations. We confirm a sub-dominant effect from pure relativistic terms, manifested at scales dominated by cosmic variance, but find that a sizable effect of order one comes from $g_{\rm NL}$ values allowed by Planck-2018 constraints, manifested at scales probed by forthcoming galaxy surveys like DESI and Euclid. As a complement, we present the matter bispectrum at the tree-level including the mentioned contributions.


title: Quantum random access codes and incompatibility of measurements

abstract: We prove that quantum random access code (QRAC) performs better than its classical counterpart only when incompatible quantum measurements are used in the decoding task. As a consequence, evaluating the average success probability for QRAC provides a semi-device independent test for the detection of quantum incompatibility. We further demonstrate that any incompatible pair of projective measurements gives an advantage over all classical strategies. Finally, we establish a connection between the maximal average success probability for QRAC and earlier quantities introduced to assess incompatibility.


title: Attending to Entities for Better Text Understanding

abstract: Recent progress in NLP witnessed the development of large-scale pre-trained language models (GPT, BERT, XLNet, etc.) based on Transformer (Vaswani et al. 2017), and in a range of end tasks, such models have achieved state-of-the-art results, approaching human performance. This demonstrates the power of the stacked self-attention architecture when paired with a sufficient number of layers and a large amount of pre-training data. However, on tasks that require complex and long-distance reasoning where surface-level cues are not enough, there is still a large gap between the pre-trained models and human performance. Strubell et al. (2018) recently showed that it is possible to inject knowledge of syntactic structure into a model through supervised self-attention. We conjecture that a similar injection of semantic knowledge, in particular, coreference information, into an existing model would improve performance on such complex problems. On the LAMBADA (Paperno et al. 2016) task, we show that a model trained from scratch with coreference as auxiliary supervision for self-attention outperforms the largest GPT-2 model, setting the new state-of-the-art, while only containing a tiny fraction of parameters compared to GPT-2. We also conduct a thorough analysis of different variants of model architectures and supervision configurations, suggesting future directions on applying similar techniques to other problems.


title: To Populate is To Regulate

abstract: We examine the effects of instantiating Lewis signaling games within a population of speaker and listener agents with the aim of producing a set of general and robust representations of unstructured pixel data. Preliminary experiments suggest that the set of representations associated with languages generated within a population outperform those generated between a single speaker-listener pair on this objective, making a case for the adoption of population-based approaches in emergent communication studies. Furthermore, post-hoc analysis reveals that population-based learning induces a number of novel factors to the conventional emergent communication setup, inviting a wide range of future research questions regarding communication dynamics and the flow of information within them.


title: Global, local and dense non-mixing of the 3D Euler equation

abstract: We prove a non-mixing property of the flow of the 3D Euler equation which has a local nature: in any neighbourhood of a "typical" steady solution there is a generic set of initial conditions, such that the corresponding Euler flows will never enter a vicinity of the original steady one. More precisely, we establish that there exist stationary solutions $u_0$ of the Euler equation on $\mathbb S^3$ and divergence-free vector fields $v_0$ arbitrarily close to $u_0$, whose (non-steady) evolution by the Euler flow cannot converge in the $C^k$ H\"older norm ($k>10$ non-integer) to any stationary state in a small (but fixed a priori) $C^k$-neighbourhood of $u_0$. The set of such initial conditions $v_0$ is open and dense in the vicinity of $u_0$. A similar (but weaker) statement also holds for the Euler flow on $\mathbb T^3$. Two essential ingredients in the proof of this result are a geometric description of all steady states near certain nondegenerate stationary solutions, and a KAM-type argument to generate knotted invariant tori from elliptic orbits.


title: The Lagrangian Mechanics and Pseudo-Periodicity of The Many-Body Planar Pendulum System

abstract: We present the Euler--Langrage equations for a many-body system of coupled planar pendulums. Hence, imposing initial condition data, the equations of motion are linearized and later developed in an idealized model for the pseudo-periodicity of the system as a function of the number of pendulums $N$. The result is empirically corroborated by comparing the model with data obtained via a numerical simulation, and by employing Kane's Method integrator in Python.


title: Conditionally Learn to Pay Attention for Sequential Visual Task

abstract: Sequential visual task usually requires to pay attention to its current interested object conditional on its previous observations. Different from popular soft attention mechanism, we propose a new attention framework by introducing a novel conditional global feature which represents the weak feature descriptor of the current focused object. Specifically, for a standard CNN (Convolutional Neural Network) pipeline, the convolutional layers with different receptive fields are used to produce the attention maps by measuring how the convolutional features align to the conditional global feature. The conditional global feature can be generated by different recurrent structure according to different visual tasks, such as a simple recurrent neural network for multiple objects recognition, or a moderate complex language model for image caption. Experiments show that our proposed conditional attention model achieves the best performance on the SVHN (Street View House Numbers) dataset with / without extra bounding box; and for image caption, our attention model generates better scores than the popular soft attention model.


title: MuSCAT2 multicolour validation of TESS candidates: an ultra-short-period substellar object around an M dwarf

abstract: We report the discovery of TOI 263.01 (TIC 120916706), a transiting substellar object (R = 0.87 RJup) orbiting a faint M3.5~V dwarf (V=18.97) on a 0.56~d orbit. We set out to determine the nature of the TESS planet candidate TOI 263.01 using ground-based multicolour transit photometry. The host star is faint, which makes RV confirmation challenging, but the large transit depth makes the candidate suitable for validation through multicolour photometry. Our analysis combines three transits observed simultaneously in r', i', and z_s bands using the MuSCAT2 multicolour imager, three LCOGT-observed transit light curves in g, r', and i' bands, a TESS light curve from Sector 3, and a low-resolution spectrum for stellar characterisation observed with the ALFOSC spectrograph. We model the light curves with PyTransit using a transit model that includes a physics-based light contamination component that allows us to estimate the contamination from unresolved sources from the multicolour photometry. This allows us to derive the true planet-star radius ratio marginalised over the contamination allowed by the photometry, and, combined with the stellar radius, gives us a reliable estimate of the object's absolute radius. The ground-based photometry excludes contamination from unresolved sources with a significant colour difference to TOI 263. Further, contamination from sources of same stellar type as the host is constrained to levels where the true radius ratio posterior has a median of 0.217. The median radius ratio corresponds to an absolute planet radius of 0.87 RJup, which confirms the substellar nature of the planet candidate. The object is either a giant planet or a brown dwarf (BD) located deep inside the so-called "brown dwarf desert". Both possibilities offer a challenge to current planet/BD formation models and makes 263.01 an object deserving of in-depth follow-up studies.


title: Nonlocal Thermoelectricity in a S-TI-S Junction in Contact with a N-Metal Probe: Evidence for Helical Edge States

abstract: We consider a Josephson junction hosting a Kramers pair of helical edge states of a quantum spin Hall bar in contact with a normal-metal probe. In this hybrid system, the orbital phase induced by a small magnetic field threading the junction known as Doppler shift (DS), combines with the conventional Josephson phase difference and originates an effect akin to a Zeeman field in the spectrum. As a consequence, when a temperature bias is applied to the superconducting terminals, a thermoelectric current is established in the normal probe. We argue that this purely non-local thermoelectric effect is a unique signature of the helical nature of the edge states coupled to superconducting leads and it can constitute a useful tool for probing the helical nature of the edge states in systems where the Hall bar configuration is difficult to achieve. We fully characterize thermoelectric response and performance of this hybrid junction in a wide range of parameters, demonstrating that the external magnetic flux inducing the DS can be used as a knob to control the thermoelectric response and the heat flow in a novel device based on topological junctions.


title: Deep-glassy ice VI confirmed with a combination of neutron spectroscopy and diffraction

abstract: The recent discovery of a low-temperature endotherm upon heating hydrochloric-acid doped ice VI has sparked a vivid controversy. The two competing explanations aiming to explain its origin range from a new distinct crystalline phase of ice to deep-glassy states of the well-known ice VI. Problems with the slow kinetics of deuterated phases have been raised, which we circumvent here entirely by simultaneously measuring the inelastic neutron spectra and neutron diffraction data of H2O samples. These measurements clearly confirm the deep-glassy ice VI scenario and rule out alternative explanations. Additionally, we show that the crystallographic model of D2O ice XV, the ordered counterpart of ice VI, also applies to the corresponding H2O phase. The discovery of deep-glassy ice VI now provides a fascinating new example of ultra-stable glasses which are encountered across a wide range of other materials.


title: Probing quantum walks through coherent control of high-dimensionally entangled photons

abstract: Quantum walks in atomic systems, owing to their continuous nature, are especially well-suited for the simulation of many-body physics and can potentially offer an exponential speedup in solving certain black box problems. Photonics offers an alternate route to simulating such nonclassical behavior in a more robust platform. However, in photonic implementations to date, an increase to the depth of a continuous quantum walk requires modifying the footprint of the system. Here we report continuous walks of a two-photon quantum frequency comb with entanglement across multiple dimensions. The coupling between frequency modes is mediated by electro-optic phase modulation, which makes the evolution of the state completely tunable over a continuous range. With arbitrary control of the phase across different modes, we demonstrate a rich variety of behavior: from walks exhibiting ballistic transport or strong energy confinement, to subspaces featuring bosonic or fermionic character. We also explore the role of entanglement dimensionality and demonstrate biphoton energy bound states, which are only possible with multilevel entanglement. This suggests the potential for such walks to quantify entanglement in high-dimensional systems.


title: Mapping Spatial Variations of HI Turbulent Properties in the Small and Large Magellanic Cloud

abstract: We developed methods for mapping spatial variations of the spatial power spectrum (SPS) and structure function (SF) slopes, with a goal of connecting neutral hydrogen (HI) statistical properties with the turbulent drivers. The new methods were applied on the HI observations of the Small and Large Magellanic Clouds (SMC and LMC). In the case of the SMC, we find highly uniform turbulent properties of HI, with no evidence for local enhancements of turbulence due to stellar feedback. Such properties could be caused by a significant turbulent driving on large-scales. Alternatively, a significant line-of-sight depth of the SMC could be masking out localized regions with a steeper SPS slope caused by stellar feedback. In contrast to the SMC, the LMC HI shows a large diversity in terms of its turbulent properties. Across most of the LMC, the small-scale SPS slope is steeper than the large-scale slope due to the presence of the HI disk. On small spatial scales, we find several areas of localized steepening of the SPS slope around major HII regions, with the 30 Doradus region being the most prominent. This is in agreement with predictions from numerical simulations which suggest steepening of the SPS slope due to stellar feedback eroding and destroying interstellar clouds. We also find localized steepening of the large-scale SPS slope in the outskirts of the LMC. This is likely caused by the flaring of the HI disk, or alternatively ram-pressure stripping of the LMC disk due to the interactions with the surrounding halo gas.


title: Bottom of spectra and coverings

abstract: We discuss the behaviour of the bottom of the spectrum of scalar Schr\"odinger operators under Riemannian coverings.


title: Information carefull worstcase DecreaseKey heaps with simple nonMeld variant

abstract: We analyze priority queues including DecreaseKey method in its interface. The paper is inspired by Strict Fibonacci Heaps [2], where G. S. Brodal, G. Lagogiannis, and R. E. Tarjan implemented the heap with DecreaseKey and Meld interface in assymptotically optimal worst case times (based on key comparisons). At the end of the paper there are mentioned possible variants of other structural properties an violations than they have used in the analysis. In the main variant a lot of information is wasted during violation reduction steps. Our goal is to concentrate on other variants and to invent natural strategy not losing that much in the information value. In other words we try to choose among them one which corresponds to superexpensive comparision principle as much as possible. The principle was described in [5] of myself, but after publication I have found these ideas in [4] of H. Kaplan, R. E. Tarjan, and U. Zwick.


title: A Combinatorial Formula for Kazhdan-Lusztig Polynomials of $\rho$-Removed Uniform Matroids

abstract: Let $\rho$ be a non-negative integer. A $\rho$-removed uniform matroid is a matroid obtained from a uniform matroid by removing a collection of $\rho$ disjoint bases. We present a combinatorial formula for Kazhdan-Lusztig polynomials of $\rho$-removed uniform matroids, using skew Young Tableaux. Even for uniform matroids, our formula is new, gives manifestly positive integer coefficients, and is more manageable than known formulas.


title: Combining human cell line transcriptome analysis and Bayesian inference to build trustworthy machine learning models for prediction of animal toxicity in drug development

abstract: Biomedical data, particularly in the field of genomics, has characteristics which make it challenging for machine learning applications - it can be sparse, high dimensional and noisy. Biomedical applications also present challenges to model selection - whilst powerful, accurate predictions are necessary, they alone are not sufficient for a model to be deemed useful. Due to the nature of the predictions, a model must also be trustworthy and transparent, empowering a practitioner with confidence that its use is appropriate and reliable. In this paper, we propose that this can be achieved through the use of judiciously built feature sets coupled with Bayesian models, specifically Gaussian processes. We apply Gaussian processes to drug discovery, using inexpensive transcriptomic profiles from human cell lines to predict animal kidney and liver toxicity after treatment with specific chemical compounds. This approach has the potential to reduce invasive and expensive animal testing during clinical trials if in vitro human cell line analysis can accurately predict model animal phenotypes. We compare results across a range of feature sets and models, to highlight model importance for medical applications.


title: Quasisymmetric orbit-flexibility of multicritical circle maps

abstract: Two given orbits of a minimal circle homeomorphism $f$ are said to be geometrically equivalent if there exists a quasisymmetric circle homeomorphism identifying both orbits and commuting with $f$. By a well-known theorem due to Herman and Yoccoz, if $f$ is a smooth diffeomorphism with Diophantine rotation number, then any two orbits are geometrically equivalent. As it follows from the a-priori bounds of Herman and Swiatek, the same holds if $f$ is a critical circle map with rotation number of bounded type. By contrast, we prove in the present paper that if $f$ is a critical circle map whose rotation number belongs to a certain full Lebesgue measure set in $(0,1)$, then the number of equivalence classes is uncountable (Theorem A). The proof of this result relies on the ergodicity of a two-dimensional skew product over the Gauss map. As a by-product of our techniques, we construct topological conjugacies between multicritical circle maps which are not quasisymmetric, and we show that this phenomenon is abundant, both from the topological and measure-theoretical viewpoints (Theorems B and C).


title: Determining the phase diagram of atomically thin layered antiferromagnet CrCl$_3$

abstract: Changes in the spin configuration of atomically-thin, magnetic van-der-Waals multilayers can cause drastic modifications in their opto-electronic properties. Conversely, the opto-electronic response of these systems provides information about the magnetic state, very difficult to obtain otherwise. Here we show that in CrCl$_3$ multilayers, the dependence of the tunnelling conductance on applied magnetic field ($H$), temperature ($T$), and number of layers ($N$) tracks the evolution of the magnetic state, enabling the magnetic phase diagram of these systems to be determined experimentally. Besides a high-field spin-flip transition occurring for all thicknesses, the in-plane magnetoconductance exhibits an even-odd effect due to a low-field spin-flop transition. If the layer number $N$ is even, the transition occurs at $\mu_0 H \sim 0$ T due to the very small in-plane magnetic anisotropy, whereas for odd $N$ the net magnetization of the uncompensated layer causes the transition to occur at finite $H$. Through a quantitative analysis of the phenomena, we determine the interlayer exchange coupling as well as the staggered magnetization, and show that in CrCl$_3$ shape anisotropy dominates. Our results reveal the rich behaviour of atomically-thin layered antiferromagnets with weak magnetic anisotropy.


title: Markov chains in random environment with applications in queueing theory and machine learning

abstract: We prove the existence of limiting distributions for a large class of Markov chains on a general state space in a random environment. We assume suitable versions of the standard drift and minorization conditions. In particular, the system dynamics should be contractive on the average with respect to the Lyapunov function and large enough small sets should exist with large enough minorization constants. We also establish that a law of large numbers holds for bounded functionals of the process. Applications to queuing systems, to machine learning algorithms and to autoregressive processes are presented.

categories: math.PR math.ST physics.data-an stat.ML stat.TH


title: DRAB-LOCUS: An Area-Efficient AES Architecture for Hardware Accelerator Co-Location on FPGAs

abstract: Advanced Encryption Standard (AES) implementations on Field Programmable Gate Arrays (FPGA) commonly focus on maximizing throughput at the cost of utilizing high volumes of FPGA slice logic. High resource usage limits systems' abilities to implement other functions (such as video processing or machine learning) that may want to share the same FPGA resources. In this paper, we address the shared resource challenge by proposing and evaluating a low-area, but high-throughput, AES architecture. In contrast to existing work, our DSP/RAM-Based Low-CLB Usage (DRAB-LOCUS) architecture leverages block RAM tiles and Digital Signal Processing (DSP) slices to implement the AES Sub Bytes, Mix Columns, and Add Round Key sub-round transformations, reducing resource usage by a factor of 3 over traditional approaches. To achieve area-efficiency, we built an inner-pipelined architecture using the internal registers of block RAM tiles and DSP slices. Our DRAB-LOCUS architecture features a 12-stage pipeline capable of producing 7.055 Gbps of interleaved encrypted or decrypted data, and only uses 909 Look Up tables, 593 Flip Flops, 16 block RAMs, and 18 DSP slices in the target device.


title: Modeling EEG data distribution with a Wasserstein Generative Adversarial Network to predict RSVP Events

abstract: Electroencephalography (EEG) data are difficult to obtain due to complex experimental setups and reduced comfort with prolonged wearing. This poses challenges to train powerful deep learning model with the limited EEG data. Being able to generate EEG data computationally could address this limitation. We propose a novel Wasserstein Generative Adversarial Network with gradient penalty (WGAN-GP) to synthesize EEG data. This network addresses several modeling challenges of simulating time-series EEG data including frequency artifacts and training instability. We further extended this network to a class-conditioned variant that also includes a classification branch to perform event-related classification. We trained the proposed networks to generate one and 64-channel data resembling EEG signals routinely seen in a rapid serial visual presentation (RSVP) experiment and demonstrated the validity of the generated samples. We also tested intra-subject cross-session classification performance for classifying the RSVP target events and showed that class-conditioned WGAN-GP can achieve improved event-classification performance over EEGNet.


title: On the Cosmological Models with Matter Creation

abstract: The matter creation model of Prigogine--G\'eh\'eniau--Gunzig--Nardone is revisited in terms of a redefined creation pressure which does not lead to irreversible adiabatic evolution at constant specific entropy. With the resulting freedom to choose a particular gas process, a flat FRWL cosmological model is proposed based on three input characteristics: (i) a perfect fluid comprising of an ideal gas, (ii) a quasi-adiabatic polytropic process, and (iii) a particular rate of particle creation. Such model leads to the description of the late-time acceleration of the expanding Universe with a natural transition from decelerating to accelerating regime. Only the Friedmann equations and the laws of thermodynamics are used and no assumptions of dark energy component is made. The model also allows the explicit determination as functions of time of all variables, including the entropy, the non-conserved specific entropy and the time the accelerating phase begins. A form of correspondence with the dark energy models (quintessence, in particular) is established via the $Om$ diagnostics. Parallels with the concordance cosmological $\Lambda$CDM model for the matter-dominated epoch and the present epoch of accelerated expansion are also established via slight modifications of both models.


title: Beyond Social Fragmentation: Coexistence of Cultural Diversity and Structural Connectivity Is Possible with Social Constituent Diversity

abstract: Social fragmentation caused by widening differences among constituents has recently become a highly relevant issue to our modern society. Theoretical models of social fragmentation using the adaptive network framework have been proposed and studied in earlier literature, which are known to either converge to a homogeneous, well-connected network or fragment into many disconnected sub-networks with distinct states. Here we introduced the diversities of behavioral attributes among social constituents and studied their effects on social network evolution. We investigated, using a networked agent-based simulation model, how the resulting network states and topologies would be affected when individual constituents' cultural tolerance, cultural state change rate, and edge weight change rate were systematically diversified. The results showed that the diversity of cultural tolerance had the most direct effect to keep the cultural diversity within the society high and simultaneously reduce the average shortest path length of the social network, which was not previously reported in the earlier literature. Diversities of other behavioral attributes also had effects on final states of the social network, with some nonlinear interactions. Our results suggest that having a broad distribution of cultural tolerance levels within society can help promote the coexistence of cultural diversity and structural connectivity.


title: GRASS: Graph Spectral Sparsification Leveraging Scalable Spectral Perturbation Analysis

abstract: Spectral graph sparsification aims to find ultra-sparse subgraphs whose Laplacian matrix can well approximate the original Laplacian eigenvalues and eigenvectors. In recent years, spectral sparsification techniques have been extensively studied for accelerating various numerical and graph-related applications. Prior nearly-linear-time spectral sparsification methods first extract low-stretch spanning tree from the original graph to form the backbone of the sparsifier, and then recover small portions of spectrally-critical off-tree edges to the spanning tree to significantly improve the approximation quality. However, it is not clear how many off-tree edges should be recovered for achieving a desired spectral similarity level within the sparsifier. Motivated by recent graph signal processing techniques, this paper proposes a similarity-aware spectral graph sparsification framework that leverages efficient spectral off-tree edge embedding and filtering schemes to construct spectral sparsifiers with guaranteed spectral similarity (relative condition number) level. An iterative graph densification scheme is also introduced to facilitate efficient and effective filtering of off-tree edges for highly ill-conditioned problems. The proposed method has been validated using various kinds of graphs obtained from public domain sparse matrix collections relevant to VLSI CAD, finite element analysis, as well as social and data networks frequently studied in many machine learning and data mining applications. For instance, a sparse SDD matrix with 40 million unknowns and 180 million nonzeros can be solved (1E-3 accuracy level) within two minutes using a single CPU core and about 6GB memory.


title: RAD: On-line Anomaly Detection for Highly Unreliable Data

abstract: Classification algorithms have been widely adopted to detect anomalies for various systems, e.g., IoT, cloud and face recognition, under the common assumption that the data source is clean, i.e., features and labels are correctly set. However, data collected from the wild can be unreliable due to careless annotations or malicious data transformation for incorrect anomaly detection. In this paper, we present a two-layer on-line learning framework for robust anomaly detection (RAD) in the presence of unreliable anomaly labels, where the first layer is to filter out the suspicious data, and the second layer detects the anomaly patterns from the remaining data. To adapt to the on-line nature of anomaly detection, we extend RAD with additional features of repetitively cleaning, conflicting opinions of classifiers, and oracle knowledge. We on-line learn from the incoming data streams and continuously cleanse the data, so as to adapt to the increasing learning capacity from the larger accumulated data set. Moreover, we explore the concept of oracle learning that provides additional information of true labels for difficult data points. We specifically focus on three use cases, (i) detecting 10 classes of IoT attacks, (ii) predicting 4 classes of task failures of big data jobs, (iii) recognising 20 celebrities faces. Our evaluation results show that RAD can robustly improve the accuracy of anomaly detection, to reach up to 98% for IoT device attacks (i.e., +11%), up to 84% for cloud task failures (i.e., +20%) under 40% noise, and up to 74% for face recognition (i.e., +28%) under 30% noisy labels. The proposed RAD is general and can be applied to different anomaly detection algorithms.


title: Provably Convergent Two-Timescale Off-Policy Actor-Critic with Function Approximation

abstract: We present the first provably convergent two-timescale off-policy actor-critic algorithm (COF-PAC) with function approximation. Key to COF-PAC is the introduction of a new critic, the emphasis critic, which is trained via Gradient Emphasis Learning (GEM), a novel combination of the key ideas of Gradient Temporal Difference Learning and Emphatic Temporal Difference Learning. With the help of the emphasis critic and the canonical value function critic, we show convergence for COF-PAC, where the critics are linear and the actor can be nonlinear.


title: Visualizing and Understanding Self-attention based Music Tagging

abstract: Recently, we proposed a self-attention based music tagging model. Different from most of the conventional deep architectures in music information retrieval, which use stacked 3x3 filters by treating music spectrograms as images, the proposed self-attention based model attempted to regard music as a temporal sequence of individual audio events. Not only the performance, but it could also facilitate better interpretability. In this paper, we mainly focus on visualizing and understanding the proposed self-attention based music tagging model.


title: Fault Detection and Identification using Bayesian Recurrent Neural Networks

abstract: In processing and manufacturing industries, there has been a large push to produce higher quality products and ensure maximum efficiency of processes. This requires approaches to effectively detect and resolve disturbances to ensure optimal operations. While the control system can compensate for many types of disturbances, there are changes to the process which it still cannot handle adequately. It is therefore important to further develop monitoring systems to effectively detect and identify those faults such that they can be quickly resolved by operators. In this paper, a novel probabilistic fault detection and identification method is proposed which adopts a newly developed deep learning approach using Bayesian recurrent neural networks~(BRNNs) with variational dropout. The BRNN model is general and can model complex nonlinear dynamics. Moreover, compared to traditional statistic-based data-driven fault detection and identification methods, the proposed BRNN-based method yields uncertainty estimates which allow for simultaneous fault detection of chemical processes, direct fault identification, and fault propagation analysis. The outstanding performance of this method is demonstrated and contrasted to (dynamic) principal component analysis, which are widely applied in the industry, in the benchmark Tennessee Eastman process~(TEP) and a real chemical manufacturing dataset.


title: Analyzing second order stochasticity of neural spiking under stimuli-bundle exposure

abstract: Conventional analysis of neuroscience data involves computing average neural activity over a group of trials and/or a period of time. This approach may be particularly problematic when assessing the response patterns of neurons to more than one simultaneously presented stimulus. In such cases, the brain must represent each individual component of the stimuli bundle, but trial-and-time-pooled averaging methods are fundamentally unequipped to address the means by which multi-item representation occurs. We introduce and investigate a novel statistical analysis framework that relates the firing pattern of a single cell, exposed to a stimuli bundle, to the ensemble of its firing patterns under each constituent stimulus. Existing statistical tools focus on what may be called "first order stochasticity" in trial-to-trial variation in the form of unstructured noise around a fixed firing rate curve associated with a given stimulus. Our analysis is based upon the theoretical premise that exposure to a stimuli bundle induces additional stochasticity in the cell's response pattern, in the form of a stochastically varying recombination of its single stimulus firing rate curves. We discuss challenges to statistical estimation of such "second order stochasticity" and address them with a novel dynamic admixture Poisson process (DAPP) model. DAPP is a hierarchical point process model that decomposes second order stochasticity into a Gaussian stochastic process and a random vector of interpretable features, and, facilitates borrowing of information on the latter across repeated trials through latent clustering. We present empirical evidence of the utility of the DAPP analysis with synthetic and real neural recordings.


title: Finite bulk Josephson currents and chirality blockade removal from inter-orbital pairing in magnetic Weyl semimetals

abstract: Magnetic Weyl semimetals (WSMs) have been presumed to be immune to proximity-induced spin-singlet superconducting pairing due to a chirality blockade of the regular Andreev reflection. In this work we study all possible superconducting pairing induced in a WSM sandwiched between two conventional superconductors in a Josephson junction setup. We confirm that, while conventional intra-orbital $s$-wave pairing is present on the surface of the WSM, it cannot propagate into the bulk due to the chirality blockade. However, inter-orbital $s$-wave pairing, with both even-frequency spin-singlet and odd-frequency mixed-spin triplet symmetry, propagates into the bulk of the WSM, as do several $p$-wave symmetries. To demonstrate the importance of these finite inter-orbital and $p$-wave pair amplitudes in an experimental setup, we calculate the Josephson current and find a finite and even increasing current when the chirality blockade effect for the conventional intra-orbital pairing is enhanced.


title: A Biologically Plausible Benchmark for Contextual Bandit Algorithms in Precision Oncology Using in vitro Data

abstract: Precision oncology, the genetic sequencing of tumors to identify druggable targets, has emerged as the standard of care in the treatment of many cancers. Nonetheless, due to the pace of therapy development and variability in patient information, designing effective protocols for individual treatment assignment in a sample-efficient way remains a major challenge. One promising approach to this problem is to frame precision oncology treatment as a contextual bandit problem and to apply sequential decision-making algorithms designed to minimize regret in this setting. However, a clear prerequisite for considering this methodology in high-stakes clinical decisions is careful benchmarking to understand realistic costs and benefits. Here, we propose a benchmark dataset to evaluate contextual bandit algorithms based on real in vitro drug response of approximately 900 cancer cell lines. Specifically, we curated a dataset of complete treatment responses for a subset of 7 treatments from prior in vitro studies. This allows us to compute the regret of proposed decision policies using biologically plausible counterfactuals. We ran a suite of Bayesian bandit algorithms on our benchmark, and found that the methods accumulate less regret over a sequence of treatment assignment tasks than a rule-based baseline derived from current clinical practice. This effect was more pronounced when genomic information was included as context. We expect this work to be a starting point for evaluation of both the unique structural requirements and ethical implications for real-world testing of bandit based clinical decision support.


title: Least-squares solutions of boundary-value problems in hybrid systems

abstract: This paper looks to apply the mathematical framework of the Theory of Functional Connections to the solution of boundary-value problems arising from hybrid systems. The Theory of Functional Connections is a technique to derive constrained expressions which are analytical expressions with embedded constraints. These expressions are particularly suitable to transform a large class of constrained optimization problems into unconstrained problems. The initial and most useful application of this technique is in the solution of differential equations where the problem can be posed as an unconstrained optimization problem and solved with simple numerical techniques (i.e. least-squares). A hybrid system is simply a sequence of different differential equations. The approach developed in this work derives an analytical constrained expression for the entire range of a hybrid system, enforcing both the boundary conditions as well as the continuity conditions across the sequence of differential equations. This reduces the searched solution space of the hybrid system to only admissible solutions. The transformation allows for a least-squares solution of the sequence for linear differential equations and a iterative least-squares solution for nonlinear differential equations. Lastly, the general formulation for "n" segments is developed and validation is provided through numerical tests for three differential equation sequences: linear/linear, linear/nonlinear, and nonlinear/nonlinear. The accuracy level obtained are all at machine-error, which is consistent with the accuracy experienced in past studies on the application of the Theory of Functional Connections to solve single ordinary differential equations.


title: Hybrid $P_{\ell}(k)$: general, unified, non-linear matter power spectrum in redshift space

abstract: Constraints on gravity and cosmology will greatly benefit from performing joint clustering and weak lensing analyses on large-scale structure data sets. Utilising non-linear information coming from small physical scales can greatly enhance these constraints. At the heart of these analyses is the matter power spectrum. Here we employ a simple method, dubbed "Hybrid $P_\ell(k)$", based on the Gaussian Streaming Model (GSM), to calculate the quasi non-linear redshift space matter power spectrum multipoles. This employs a fully non-linear and theoretically general prescription for the matter power spectrum. We test this approach against comoving Lagrangian acceleration simulation measurements performed in GR, DGP and $f(R)$ gravity and find that our method performs comparably or better to the dark matter TNS redshift space power spectrum model for dark matter. When comparing the redshift space multipoles for halos, we find that the Gaussian approximation of the GSM with a linear bias and a free stochastic term, $N$, is competitive to the TNS model. Our approach offers many avenues for improvement in accuracy as well as further unification under the halo model.


title: Discrete variants of Brunn-Minkowski type inequalities

abstract: We present an alternative, short proof of a recent discrete version of the Brunn-Minkowski inequality due to Lehec and the second named author. Our proof also yields the four functions theorem of Ahlswede and Daykin as well as some new variants.


title: Simplifying Random Forests: On the Trade-off between Interpretability and Accuracy

abstract: We analyze the trade-off between model complexity and accuracy for random forests by breaking the trees up into individual classification rules and selecting a subset of them. We show experimentally that already a few rules are sufficient to achieve an acceptable accuracy close to that of the original model. Moreover, our results indicate that in many cases, this can lead to simpler models that clearly outperform the original ones.


title: On the log minimal model program for $3$-folds over imperfect fields of characteristic $p>5$

abstract: We prove that many of the results of the LMMP hold for $3$-folds over fields of characteristic $p>5$ which are not necessarily perfect. In particular, the existence of flips, the cone theorem, the contraction theorem for birational extremal rays, and the existence of log minimal models. As well as pertaining to the geometry of fibrations of relative dimension $3$ over algebraically closed fields, they have applications to tight closure in dimension $4$.