title: Energy-Efficient Offloading in Delay-Constrained Massive MIMO Enabled Edge Network Using Data Partitioning

abstract: We study a wireless edge-computing system which allows multiple users to simultaneously offload computation-intensive tasks to multiple massive-MIMO access points, each with a collocated multi-access edge computing (MEC) server. Massive-MIMO enables simultaneous uplink transmissions from all users, significantly shortening the data offloading time compared to sequential protocols, and makes the three phases of data offloading, computing, and downloading have comparable durations. Based on this three-phase structure, we formulate a novel problem to minimize a weighted sum of the energy consumption at both the users and the MEC server under a round-trip latency constraint, using a combination of data partitioning, transmit power control and CPU frequency scaling at both the user and server ends. We design a novel nested primal-dual algorithm using two different methods to solve this problem efficiently. Optimized solutions show that for larger requests, more data is offloaded to the MECs to reduce local computation time in order to meet the latency constraint, despite higher energy cost of wireless transmissions. Massive-MIMO channel estimation errors under pilot contamination also causes more data to be offloaded to the MECs. Compared to binary offloading, partial offloading with data partitioning is superior and leads to significant reduction in the overall energy consumption.


title: Hidden conformal symmetry for Kerr-Newman-NUT-AdS black holes

abstract: We show the generic non-extremal Kerr-Newman-NUT-AdS black holes are holographically dual to the hidden conformal field theories in two different pictures. The two pictures can be merged together to the dual CFTs in general picture, that are generated by the $SL(2,\mathbb{Z})$ modular group. We also extend the calculation to the extremal limit and find the corresponding quantities. Using the extended central charges from the Kerr/CFT correspondence, we also find agreement between the macroscopic and microscopic entropies. We also find the absorption cross-section of the scalar probes for the generic and extremal Kerr-Newman-NUT-AdS black holes, to further support the different dual CFTs to the black holes.


title: Computation of effective front form Hamiltonians for massive Abelian gauge theory

abstract: Renormalization group procedure for effective particles (RGPEP) is applied in terms of a second-order perturbative computation to an Abelian gauge theory, as an example of application worth studying on the way toward derivation of a dynamical connection between the spectroscopy of bound states and their parton-model picture in the front form of Hamiltonian dynamics. In addition to the ultraviolet transverse divergences that are handled using the RGPEP in previously known ways, the small-x divergences are handled by introducing a mass parameter and a third polarization state for gauge bosons using a mechanism analogous to spontaneous breaking of global gauge symmetry, in a special limit that simplifies the theory to Soper's front form of massive QED. The resulting orders of magnitude of scales involved in the dynamics of effective constituents or partons in the simplified theory are identified for the fermion and boson mass counter terms, effective masses and self-interactions, as well as for the Coulomb-like effective interactions in bound states of fermions. Computations in orders higher than second are mentioned but not described in this article.


title: Uniform propagation of chaos for the thermostated Kac model

abstract: We consider Kac's 1D N-particle system coupled to an ideal thermostat at temperature T, introduced by Bonetto, Loss, and Vaidyanathan in 2014. We obtain a propagation of chaos result for this system, with explicit and uniform-in-time rates of order N^(-1/3) in the 2-Wasserstein metric. We also show well-posedness and equilibration for the limit kinetic equation in the space of probability measures. The proofs use a coupling argument previously introduced by Cortez and Fontbona in 2016.


title: X-CIGALE: fitting AGN/galaxy SEDs from X-ray to infrared

abstract: CIGALE is a powerful multiwavelength spectral energy distribution (SED) fitting code for extragalactic studies. However, the current version of CIGALE is not able to fit X-ray data, which often provide unique insights into AGN intrinsic power. We develop a new X-ray module for CIGALE, allowing it to fit SEDs from the X-ray to infrared (IR). We also improve the AGN fitting of CIGALE from UV-to-IR wavelengths. We implement a modern clumpy two-phase torus model, SKIRTOR. To account for moderately extincted type 1 AGNs, we implement polar-dust extinction. We publicly release the source code (named X-CIGALE). We test X-CIGALE with X-ray detected AGNs in SDSS, COSMOS, and AKARI-NEP. The fitting quality (as indicated by reduced $\chi^2$) is good in general, indicating that X-CIGALE is capable of modelling the observed SED from X-ray to IR. We discuss constrainability and degeneracy of model parameters in the fitting of AKARI-NEP, for which excellent mid-IR photometric coverage is available. We also test fitting a sample of AKARI-NEP galaxies for which only X-ray upper limits are available from Chandra observations, and find that the upper limit can effectively constrain the AGN SED contribution for some systems. Finally, using X-CIGALE, we assess the ability of Athena to constrain the AGN activity in future extragalactic studies.


title: Heart rate variability code: Does it exist and can we hack it?

abstract: Heart rate variability (HRV) has been studied for over 50 years, yet an integrative concept is missing on what HRV's mathematical properties represent physiologically. Here I introduce the notion of HRV code as an attempt to address this challenge systematically. I review the existing evidence from physiological studies in various species to support this concept and propose experiments to help validate and expand this notion further.


title: Lipschitz regularity of the invariant measure of Random Dynamical Systems

abstract: We consider a class of Random Dynamical Systems and prove that its invariant measure has a Lipschitz disintegration. As a consequence, we prove exponential decay of correlations over regular observables. Both results are obtained as an application of spectral gap for its transfers operator which is also proved in the paper.


title: Searching for central stars of planetary nebulae in Gaia DR2

abstract: Context: Accurate distance measurements are fundamental to the study of Planetary Nebulae (PNe) but have long been elusive. The most accurate and model-independent distance measurements for galactic PNe come from the trigonometric parallaxes of their central stars, which were only available for a few tens of objects prior to the Gaia mission. Aims: Accurate identification of PN central stars in the Gaia source catalogues is a critical prerequisite for leveraging the unprecedented scope and precision of the trigonometric parallaxes measured by Gaia. Our aim is to build a complete sample of PN central star detections with minimal contamination. Methods: We develop and apply an automated technique based on the likelihood ratio method to match candidate central stars in Gaia Data Release 2 (DR2) to known PNe in the Hong Kong/AAO/Strasbourg H$\alpha$ (HASH) PN catalogue, taking into account the BP--RP colours of the \emph{Gaia} sources as well as their positional offsets from the nebula centres. These parameter distributions for both true central stars and background sources are inferred directly from the data. Results: We present a catalogue of over 1000 Gaia sources that our method has automatically identified as likely PN central stars. We demonstrate how the best matches enable us to trace nebula and central star evolution and to validate existing statistical distance scales, and discuss the prospects for further refinement of the matching based on additional data. We also compare the accuracy of our catalogue to that of previous works.


title: Tailoring a pair of pants

abstract: We show how to deform the map $\operatorname{Log}\colon (\mathbb{C}^*)^n \to \mathbb{R}^n$ such that the image of the complex pair of pants $P^\circ \subset {(\mathbb{C}^*)^n}$ is the tropical hyperplane by showing an (ambient) isotopy between $P^\circ \subset {(\mathbb{C}^*)^n}$ and a natural polyhedral subcomplex of the product of the two skeleta $S\times \Sigma \subset \mathcal{A} \times \mathcal{C}$ of the amoeba $\mathcal{A}$ and the coamoeba $\mathcal{C}$ of $P^\circ$. This lays the groundwork for having the discriminant to be of codimension 2 in topological Strominger-Yau-Zaslow torus fibrations.


title: An asymptotic preserving semi-implicit multiderivative solver

abstract: In this work we construct a multiderivative implicit-explicit (IMEX) scheme for a class of stiff ordinary differential equations. Our solver is high-order accurate and has an asymptotic preserving (AP) property. The proposed method is based upon a two-derivative backward Taylor series base solver, which we show has an AP property. Higher order accuracies are found by iterating the result over a high-order multiderivative interpolant of the right hand side function, which we again prove has an AP property. Theoretical results showcasing the asymptotic consistency as well as the high-order accuracy of the solver are presented. In addition, an extension of the solver to an arbitrarily split right hand side function is also offered. Numerical results for a collection of standard test cases from the literature are presented that support the theoretical findings of the paper.


title: Representation Learning for Medical Data

abstract: We propose a representation learning framework for medical diagnosis domain. It is based on heterogeneous network-based model of diagnostic data as well as modified metapath2vec algorithm for learning latent node representation. We compare the proposed algorithm with other representation learning methods in two practical case studies: symptom/disease classification and disease prediction. We observe a significant performance boost in these task resulting from learning representations of domain data in a form of heterogeneous network.


title: Cartan subalgebras for non-principal twisted groupoid $C^*$-algebras

abstract: Renault proved in 2008 that if $G$ is a topologically principal groupoid, then $C_0(G^{(0)})$ is a Cartan subalgebra in $C^*_r(G, \Sigma)$ for any twist $\Sigma$ over $G$. However, there are many groupoids which are not topologically principal, yet their (twisted) $C^*$-algebras admit Cartan subalgebras. This paper gives a dynamical description of a class of such Cartan subalgebras, by identifying conditions on a 2-cocycle $c$ on $G$ and a subgroupoid $S \subseteq G$ under which $C^*_r(S, c)$ is Cartan in $C^*_r(G, c)$. When $G$ is a discrete group, we also describe the Weyl groupoid and twist associated to these Cartan pairs, under mild additional hypotheses.


title: To quantum or not to quantum: towards algorithm selection in near-term quantum optimization

abstract: The Quantum Approximate Optimization Algorithm (QAOA) constitutes one of the often mentioned candidates expected to yield a quantum boost in the era of near-term quantum computing. In practice, quantum optimization will have to compete with cheaper classical heuristic methods, which have the advantage of decades of empirical domain-specific enhancements. Consequently, to achieve optimal performance we will face the issue of algorithm selection, well-studied in practical computing. Here we introduce this problem to the quantum optimization domain. Specifically, we study the problem of detecting those problem instances of where QAOA is most likely to yield an advantage over a conventional algorithm. As our case study, we compare QAOA against the well-understood approximation algorithm of Goemans and Williamson (GW) on the Max-Cut problem. As exactly predicting the performance of algorithms can be intractable, we utilize machine learning to identify when to resort to the quantum algorithm. We achieve cross-validated accuracy well over 96\%, which would yield a substantial practical advantage. In the process, we highlight a number of features of instances rendering them better suited for QAOA. While we work with simulated idealised algorithms, the flexibility of ML methods we employed provides confidence that our methods will be equally applicable to broader classes of classical heuristics, and to QAOA running on real-world noisy devices.


title: Humps and bumps: The effects of shocks on the optical light curves of fundamental-mode RR Lyrae stars

abstract: We present the most extended and homogeneous study carried out so far of the main and early shocks in 1485 RR~Lyrae stars in the Galactic bulge observed by the Optical Gravitational Lensing Experiment (OGLE). We selected non-modulated fundamental-mode RR~Lyrae stars with good-quality photometry. Using a self-developed method, we determined the centers and strengths of main and early shock features in the phased light curves. We found that the position of both humps and bumps are highly correlated with the pulsation properties of the studied variables. Pulsators with a pronounced main shock are concentrated in the low-amplitude regime of the period-amplitude diagram, while stars with a strong early shock have average and above-average pulsation amplitudes. A connection between the main and early shocks and the Fourier coefficients is also observed. In the color-magnitude diagram (CMD), we see a separation between stars with strong and weak shocks. Variables with a pronounced main shock cluster close to the fundamental red edge of the instability strip (IS), while stars with a strong early shock tend to clump in the center and near the fundamental blue edge of the IS. The appearance of shocks and their properties seem independent of the direction of evolution estimated from the period change rate of the studied stars. In addition, the differences in the period change rate between the two main Oosterhoff groups found in the Galactic bulge suggest that stars of Oosterhoff type I are located close to the zero-age horizontal branch while Oosterhoff type II variables are on their way toward the fundamental red edge of the instability strip, thus having already left the zero-age horizontal branch.


title: Gradient Tomography of Jet Quenching in Heavy-Ion Collisions

abstract: Transverse momentum broadening and energy loss of a propagating parton are dictated by the space-time profile of the jet transport coefficient $\hat q$ in a dense QCD medium. The spatial gradient of $\hat q$ perpendicular to the propagation direction can lead to a drift and asymmetry in parton transverse momentum distribution. Such an asymmetry depends on both the spatial position along the transverse gradient and path length of a propagating parton as shown by numerical solutions of the Boltzmann transport in the simplified form of a drift-diffusion equation. In high-energy heavy-ion collisions, this asymmetry with respect to a plane defined by the beam and trigger particle (photon, hadron or jet) with a given orientation relative to the event plane is shown to be closely related to the transverse position of the initial jet production in full event-by-event simulations within the linear Boltzmann transport model. Such a gradient tomography can be used to localize the initial jet production position for more detailed study of jet quenching and properties of the quark-gluon plasma along a given propagation path in heavy-ion collisions.


title: Gauge boson mass as regulator

abstract: Small-x divergences of Abelian gauge theory in the front form of Hamiltonian dynamics are regulated using a mass parameter for gauge bosons, introduced through a mechanism analogous to the spontaneous breaking of global gauge symmetry. A corresponding family of ultraviolet and infrared finite scale-dependent renormalized Hamiltonians, is calculable order-by-order using the renormalization group procedure for effective particles. The second-order terms described here suggest the magnitude of mass corrections that may be involved in resolving the small-x parton and front-form vacuum and zero-mode problems, assuming that the gauge boson mass that counts does not exceed the current upper bound on the photon mass.


title: Learning Discontinuous Piecewise Affine Fitting Functions using Mixed Integer Programming for Segmentation and Denoising

abstract: Piecewise affine functions are widely used to approximate nonlinear and discontinuous functions. However, most, if not all existing models only deal with fitting continuous functions. In this paper, we investigate the problem of fitting a discontinuous piecewise affine function to given data that lie in an orthogonal grid, where no restriction on the partition is enforced (i.e., its geometric shape can be nonconvex). This is useful for segmentation and denoising when data corresponding to images. We propose a novel Mixed Integer Program (MIP) formulation for the piecewise affine fitting problem, where binary variables determine the location of break-points. To obtain consistent partitions (i.e. image segmentation), we include multi-cut constraints in the formulation. Since the resulting problem is $\mathcal{NP}$-hard, two techniques are introduced to improve the computation. One is to add facet-defining inequalities to the formulation and the other to provide initial integer solutions using a special heuristic algorithm. We conduct extensive experiments by some synthetic images as well as real depth images, and the results demonstrate the feasibility of our model.


title: On the convergence of multi-level Hermite-Pad\'e approximants for a class of meromorphic functions

abstract: The present paper deals with the convergence properties of multi-level Hermite-Pad\'e approximants for a class of meromorphic functions given by rational perturbations with real coefficients of a Nikishin system of functions, and study the zero location of the corresponding multiple orthogonal polynomials.


title: Intermittent Pulling with Local Compensation for Communication-Efficient Federated Learning

abstract: Federated Learning is a powerful machine learning paradigm to cooperatively train a global model with highly distributed data. A major bottleneck on the performance of distributed Stochastic Gradient Descent (SGD) algorithm for large-scale Federated Learning is the communication overhead on pushing local gradients and pulling global model. In this paper, to reduce the communication complexity of Federated Learning, a novel approach named Pulling Reduction with Local Compensation (PRLC) is proposed. Specifically, each training node intermittently pulls the global model from the server in SGD iterations, resulting in that it is sometimes unsynchronized with the server. In such a case, it will use its local update to compensate the gap between the local model and the global model. Our rigorous theoretical analysis of PRLC achieves two important findings. First, we prove that the convergence rate of PRLC preserves the same order as the classical synchronous SGD for both strongly-convex and non-convex cases with good scalability due to the linear speedup with respect to the number of training nodes. Second, we show that PRLC admits lower pulling frequency than the existing pulling reduction method without local compensation. We also conduct extensive experiments on various machine learning models to validate our theoretical results. Experimental results show that our approach achieves a significant pulling reduction over the state-of-the-art methods, e.g., PRLC requiring only half of the pulling operations of LAG.


title: LECX: a cubesat experiment to detect and localize cosmic explosions in hard X rays

abstract: With the advent of the nanosat/cubesat revolution, new opportunities have appeared to develop and launch small ($\sim$\ts 1000 cm$^3$), low-cost ($\sim$\ts US\$ 1M) experiments in space in very short timeframes ($\sim$ 2\ts years). In the field of high-energy astrophysics, in particular, it is a considerable challenge to design instruments with compelling science and competitive capabilities that can fit in very small satellite buses such as a cubesat platform, and operate them with very limited resources. Here we describe a hard X-ray (30--200\ts keV) experiment, LECX ("Localizador de Explos\~oes C\'osmicas de Raios X" -- Locator of X-Ray Cosmic Explosions), that is capable of detecting and localizing within a few degrees events like Gamma-Ray Bursts and other explosive phenomena in a 2U-cubesat platform, at a rate of $\sim${\bf 5 events year$^{-1}$.} In the current gravitational wave era of astronomy, a constellation or swarm of small spacecraft carrying instruments such as LECX can be a very cost-effective way to search for electromagnetic counterparts of gravitational wave events produced by the coalescence of compact objects.


title: Transition-Based Dependency Parsing using Perceptron Learner

abstract: Syntactic parsing using dependency structures has become a standard technique in natural language processing with many different parsing models, in particular data-driven models that can be trained on syntactically annotated corpora. In this paper, we tackle transition-based dependency parsing using a Perceptron Learner. Our proposed model, which adds more relevant features to the Perceptron Learner, outperforms a baseline arc-standard parser. We beat the UAS of the MALT and LSTM parsers. We also give possible ways to address parsing of non-projective trees.


title: A Recent Starbust in the Low Surface Brightness Galaxy UGC 628

abstract: We present the star-formation history of the low surface brightness (LSB) galaxy UGC 628 as part of the MUSCEL program (MUltiwavelength observations of the Structure, Chemistry, and Evolution of LSB galaxies). The star-formation histories of LSB galaxies represent a significant gap in our knowledge of galaxy assembly, with implications for dark matter / baryon feedback, IGM gas accretion, and the physics of star formation in low metallicity environments. Our program uses ground-based IFU spectra in tandem with space-based UV and IR imaging to determine the star-formation histories of LSB galaxies in a spatially resolved fashion. In this work we present the fitted history of our first target to demonstrate our techniques and methodology. Our technique splits the history of this galaxy into 15 semi-logarithmically spaced timesteps. Within each timestep the star-formation rate of each spaxel is assumed constant. We then determine the set of 15 star-formation rates that best recreate the spectra and photometry measured in each spaxel. Our main findings with respect to UGC 628 are: a) the visible properties of UGC 628 have varied over time, appearing as a high surface brightness spiral earlier than 8 Gyr ago and a starburst galaxy during a recent episode of star formation several tens of Myr ago, b) the central bar/core region was established early, around 8-10 Gyr ago, but has been largely inactive since, and c) star formation in the past 3 Gyr is best characterised as patchy and sporadic.


title: Convolutional Codes

abstract: The article provides a survey on convolutional codes stressing the connections to module theory and systems theory. Constructions of codes with maximal possible distance and distance profile are provided. The article will appear as book chapter in "A Concise Encyclopedia of Coding Theory" to be published by CRC Press.


title: Mutual Inclination Excitation by Stellar Oblateness

abstract: Ultra-short-period planets (USPs) provide important clues to planetary formation and migration. Recently, it is found that the mutual inclinations of the planetary systems are larger if the inner orbits are closer ($\lesssim 5R_*$) and if the planetary period ratios are larger ($P_2/P_1 \gtrsim 5$) (Dai et al. 2018). This suggests that the USPs experienced both inclination excitation and orbital shrinkage. Here we investigate the increase in the mutual inclination due to stellar oblateness. We find that the stellar oblateness (within $\sim 1$Gyr) is sufficient to enhance the mutual inclination to explain the observed signatures. This suggests that the USPs can migrate closer to the host star in a near coplanar configuration with their planetary companions (e.g., disk migration+tides or in-situ+tides), before mutual inclination gets excited due to stellar oblateness.


title: Engineering Weyl phases and nonlinear Hall effects in T$_d$-MoTe$_2$

abstract: MoTe$_2$ has recently attracted much attention due to the observation of pressure-induced superconductivity, exotic topological phase transitions, and nonlinear quantum effects. However, there has been debate on the intriguing structural phase transitions among various observed phases of MoTe$_2$, and their connection to the underlying topological electronic properties. In this work, by means of density-functional theory (DFT+U) calculations, we investigate the structural phase transition between the polar T$_d$ and nonpolar 1T$'$ phases of MoTe$_2$ in reference to a hypothetical high-symmetry T$_0$ phase that exhibits higher-order topological features. In the T$_d$ phase we obtain a total of 12 Weyl points, which can be created/annihilated, dynamically manipulated, and switched by tuning a polar phonon mode. We also report the existence of a tunable nonlinear Hall effect in T$_d$-MoTe$_2$, and propose the use of this effect as a probe for the detection of polarity orientation in polar (semi)metals. By studying the role of dimensionality, we identify a configuration in which a nonlinear surface response current emerges. The potential technological applications of the tunable Weyl phase and the nonlinear Hall effect are discussed.


title: Dynamics of extended Schelling models

abstract: We explore extensions of Schelling's model of social dynamics, in which two types of agents live on a checkerboard lattice and move in order to optimize their own satisfaction, which depends on how many agents among their neighbors are of their same type. For each number $n$ of same-type nearest neighbors we independently assign a binary satisfaction variable $s_{k}$ which is equal to one only if the agent is satisfied with that condition, and is equal to zero otherwise. This defines 32 different satisfaction rules, which we investigate in detail, focusing on pattern formation and measuring segregation with the help of an "energy" function which is related to the number of neighboring agents of different types and plays no role in the dynamics. We consider the checkerboard lattice to be fully occupied and the dynamics consists of switching the locations of randomly selected unsatisfied agents of opposite types. We show that, starting from a random distribution of agents, only a small number of rules lead to (nearly) fully segregated patterns in the long run, with many rules leading to chaotic steady-state behavior. Nevertheless, other interesting patterns may also be dynamically generated, such as "anti-segregate d" patterns as well as patterns resembling sponges.


title: The kinetic space of multistationarity in dual phosphorylation

abstract: Multistationarity in molecular systems underlies switch-like responses in cellular decision making. Determining whether and when a system displays multistationarity is in general a difficult problem. In this work we completely determine the set of kinetic parameters that enable multistationarity in a ubiquitous motif involved in cell signaling, namely a dual phosphorylation cycle. In addition we show that the regions of multistationarity and monostationarity are both path connected. We model the dynamics of the concentrations of the proteins over time by means of a parametrized polynomial ordinary differential equation (ODE) system arising from the mass-action assumption. Since this system has three linear first integrals defined by the total amounts of the substrate and the two enzymes, we study for what parameter values the ODE system has at least two positive steady states after suitably choosing the total amounts. We employ a suite of techniques from (real) algebraic geometry, which in particular concern the study of the signs of a multivariate polynomial over the positive orthant and sums of nonnegative circuit polynomials.


title: A Multi-Scale Tensor Network Architecture for Classification and Regression

abstract: We present an algorithm for supervised learning using tensor networks, employing a step of preprocessing the data by coarse-graining through a sequence of wavelet transformations. We represent these transformations as a set of tensor network layers identical to those in a multi-scale entanglement renormalization ansatz (MERA) tensor network, and perform supervised learning and regression tasks through a model based on a matrix product state (MPS) tensor network acting on the coarse-grained data. Because the entire model consists of tensor contractions (apart from the initial non-linear feature map), we can adaptively fine-grain the optimized MPS model backwards through the layers with essentially no loss in performance. The MPS itself is trained using an adaptive algorithm based on the density matrix renormalization group (DMRG) algorithm. We test our methods by performing a classification task on audio data and a regression task on temperature time-series data, studying the dependence of training accuracy on the number of coarse-graining layers and showing how fine-graining through the network may be used to initialize models with access to finer-scale features.


title: Wild Galois representations: a family of hyperelliptic curves with large inertia image

abstract: In this work we generalise the main result of arXiv:1812.05651 to the family of hyperelliptic curves with potentially good reduction over a $p$-adic field which have degree $p$ and the largest possible image of inertia under the $\ell$-adic Galois representation associated to its Jacobian. We will prove that this Galois representation factors as the tensor product of an unramified character and an irreducible representation of a finite group, which can be either equal to the inertia image (in which case the representation is easily determined) or a $C_2$-extension of it. In this second case, there are two suitable representations and we will describe the Galois action explicitly in order to determine the correct one.


title: Characterizing Smart Home IoT Traffic in the Wild

abstract: As the smart home IoT ecosystem flourishes, it is imperative to gain a better understanding of the unique challenges it poses in terms of management, security, and privacy. Prior studies are limited because they examine smart home IoT devices in testbed environments or at a small scale. To address this gap, we present a measurement study of smart home IoT devices in the wild by instrumenting home gateways and passively collecting real-world network traffic logs from more than 200 homes across a large metropolitan area in the United States. We characterize smart home IoT traffic in terms of its volume, temporal patterns, and external endpoints along with focusing on certain security and privacy concerns. We first show that traffic characteristics reflect the functionality of smart home IoT devices such as smart TVs generating high volume traffic to content streaming services following diurnal patterns associated with human activity. While the smart home IoT ecosystem seems fragmented, our analysis reveals that it is mostly centralized due to its reliance on a few popular cloud and DNS services. Our findings also highlight several interesting security and privacy concerns in smart home IoT ecosystem such as the need to improve policy-based access control for IoT traffic, lack of use of application layer encryption, and prevalence of third-party advertising and tracking services. Our findings have important implications for future research on improving management, security, and privacy of the smart home IoT ecosystem.


title: Substitution of subspace collections with nonorthogonal subspaces to accelerate Fast Fourier Transform methods applied to conducting composites

abstract: We show the power of the algebra of subspace collections developed in Chapter 7 of the book "Extending the Theory of Composites to Other Areas of Science (edited by Milton, 2016). Specifically we accelerate the Fast Fourier Transform schemes of Moulinec and Suquet and Eyre and Milton (1994, 1998) for computing the fields and effective tensor in a conducting periodic medium by substituting a subspace collection with nonorthogonal subspaces inside one with orthogonal subspaces. This can be done when the effective conductivity as a function of the conductivity $\sigma_1$ of the inclusion phase (with the matrix phase conductivity set to $1$) has its singularities confined to an interval $[-\beta,-\alpha]$ of the negative real $\sigma_1$ axis. Numerical results of Moulinec and Suquet show accelerated convergence for the model example of a square array of squares at $25\%$ volume fraction. For other problems we show how $Q^*_C$-convex functions can be used to restrict the region where singularities of the effective tensor as a function of the component tensors might be found.


title: Transformer-based Online CTC/attention End-to-End Speech Recognition Architecture

abstract: Recently, Transformer has gained success in automatic speech recognition (ASR) field. However, it is challenging to deploy a Transformer-based end-to-end (E2E) model for online speech recognition. In this paper, we propose the Transformer-based online CTC/attention E2E ASR architecture, which contains the chunk self-attention encoder (chunk-SAE) and the monotonic truncated attention (MTA) based self-attention decoder (SAD). Firstly, the chunk-SAE splits the speech into isolated chunks. To reduce the computational cost and improve the performance, we propose the state reuse chunk-SAE. Sencondly, the MTA based SAD truncates the speech features monotonically and performs attention on the truncated features. To support the online recognition, we integrate the state reuse chunk-SAE and the MTA based SAD into online CTC/attention architecture. We evaluate the proposed online models on the HKUST Mandarin ASR benchmark and achieve a 23.66% character error rate (CER) with a 320 ms latency. Our online model yields as little as 0.19% absolute CER degradation compared with the offline baseline, and achieves significant improvement over our prior work on Long Short-Term Memory (LSTM) based online E2E models.


title: Machine learning models for the secondary Bjerknes force between two insonated bubbles

abstract: The secondary Bjerknes force plays a significant role in the evolution of bubble clusters. However, due to the complex dependence of the force on multiple parameters, it is highly non-trivial to include the effects of this force in the simulations of bubble clusters. In this paper, machine learning is used to develop a data-driven model for the secondary Bjerknes force between two insonated bubbles as a function of the equilibrium radii of the bubbles, the distance between the bubbles, the amplitude and the frequency of the pressure. The force varies over several orders of magnitude, which poses a serious challenge for the usual machine learning models. To overcome this difficulty, the magnitudes and the signs of the force are separated and modelled separately. A nonlinear regression is obtained with a feed-forward network model for the logarithm of the magnitude, whereas the sign is modelled by a support-vector machine model. The principle, the practical aspects related to the training and validation of the machine models are introduced. The predictions from the models are checked against the values computed from the Keller-Miksis equations. The results show that the models are extremely efficient while providing accurate estimate of the force. The models make it computationally feasible for the future simulations of the bubble clusters to include the effects of the secondary Bjerknes force.


title: How many simplices are needed to triangulate a Grassmannian?

abstract: We compute a lower bound for the number of simplices that are needed to triangulate the Grassmann manifold $G_k(\mathbb{R}^n)$. In particular, we show that the number of top-dimensional simplices grows exponentially with $n$. More precise estimates are given for $k=2,3,4$. Our method can be used to estimate the minimal size of triangulations for other spaces, like Lie groups, flag manifolds, Stiefel manifolds etc.


title: "How over is it?" Understanding the Incel Community on YouTube

abstract: YouTube is by far the largest host of user-generated video content worldwide. Alas, the platform has also come under fire for hosting inappropriate, toxic, and hateful content. One community that has often been linked to sharing and publishing hateful and misogynistic content are the Involuntary Celibates (Incels), a loosely defined movement ostensibly focusing on men's issues. In this paper, we set out to analyze the Incel community on YouTube by focusing on this community's evolution over the last decade and understanding whether YouTube's recommendation algorithm steers users towards Incel-related videos. We collect videos shared on Incel communities within Reddit and perform a data-driven characterization of the content posted on YouTube. Among other things, we find that the Incel community on YouTube is getting traction and that, during the last decade, the number of Incel-related videos and comments rose substantially. We also find that users have a 6.3% chance of being suggested an Incel-related video by YouTube's recommendation algorithm within five hops when starting from a non Incel-related video. Overall, our findings paint an alarming picture of online radicalization: not only Incel activity is increasing over time, but platforms may also play an active role in steering users towards such extreme content.


title: On the structure of invertible elements in certain Fourier-Stieltjes algebras

abstract: For a locally compact abelian group $G$, J. L. Taylor (1971) gave a complete characterization of invertible elements in the measure algebra $M(G)$. Using the Fourier-Stieltjes transform, this characterization can be carried out in the context of Fourier-Stieltjes algebras $B(G)$. We obtain this latter characterization for the Fourier-Stieltjes algebra $B(G)$ of certain classes of locally compact groups, in particular, many totally minimal groups and the $ax+b$-group.


title: Models of Lubin-Tate spectra via Real bordism theory

abstract: We study certain formal group laws equipped with an action of the cyclic group of order a power of $2$. We construct $C_{2^n}$-equivariant Real oriented models of Lubin-Tate spectra $E_h$ at heights $h=2^{n-1}m$ and give explicit formulas of the $C_{2^n}$-action on their coefficient rings. Our construction utilizes equivariant formal group laws associated with the norms of the Real bordism theory $MU_{\mathbb{R}}$, and our work examines the height of the formal group laws of the Hill-Hopkins-Ravenel norms of $MU_{\mathbb{R}}$.


title: Flare-like Variability of the Mg~II $\lambda$2798 \AA\ Emission Line and UV Fe~II band in the Blazar CTA~102

abstract: We report on the detection of a statistically significant flare-like event in the Mg~II~$\lambda$ 2798~\AA\ emission line and the UV~Fe~II band of CTA~102 during the outburst of autumn 2017. The ratio between the maximum and minimum of $\lambda$3000~\AA\ continuum flux for the observation period ($2010-2017$) is 179$\pm$15. Respectively, the max/min ratios 8.1$\pm$10.5 and 34.0$\pm$45.5 confirmed the variability of the Mg~II emission line and of the Fe~II band. The highest levels of emission lines fluxes recorded coincide with a superluminal jet component traversing through a stationary component located at $\sim$0.1 mas from the 43 GHz core. Additionally, comparing the Mg~II line profile in the minimum of activity against the one in the maximum, we found that the latter is broader and blue-shifted. As a result of these findings, we can conclude that the non-thermal continuum emission produced by material in the jet moving at relativistic speeds is related to the broad emission line fluctuations. In consequence, these fluctuations are also linked to the presence of broad-line region (BLR) clouds located at $\sim$25 pc from the central engine, outside from the inner parsec, where the canonical BLR is located. Our results suggest that during strong activity in CTA~102, the source of non-thermal emission and broad-line clouds outside the inner parsec introduces uncertainties in the estimates of black hole (BH) mass. Therefore, it is important to estimate the BH mass, using single-epoch or reverberation mapping techniques, only with spectra where the continuum luminosity is dominated by the accretion disk.


title: The establishment of the Standard Cosmological Model through observations

abstract: Over the last decades, observations with increasing quality have revolutionized our understanding of the general properties of the Universe. Questions posed for millenia by mankind about the origin, evolution and structure of the cosmos have found an answer. This has been possible mainly thanks to observations of the Cosmic Microwave Background, of the large-scale distribution of matter structure in the local Universe, and of type Ia supernovae that have revealed the accelerated expansion of the Universe. All these observations have successfully converged into the so-called "concordance model". In spite of all these observational successes, there are still some important open problems, the most obvious of which are what generated the initial matter inhomogeneities that led to the structure observable in today's Universe, and what is the nature of dark matter, and of the dark energy that drives the accelerated expansion. In this chapter I will expand on the previous aspects. I will present a general description of the Standard Cosmological Model of the Universe, with special emphasis on the most recent observations that have us allowed to consolidate this model. I will also discuss the shortfalls of this model, its most pressing open questions, and will briefly describe the observational programmes that are being planned to tackle these issues.


title: Proxy Tasks and Subjective Measures Can Be Misleading in Evaluating Explainable AI Systems

abstract: Explainable artificially intelligent (XAI) systems form part of sociotechnical systems, e.g., human+AI teams tasked with making decisions. Yet, current XAI systems are rarely evaluated by measuring the performance of human+AI teams on actual decision-making tasks. We conducted two online experiments and one in-person think-aloud study to evaluate two currently common techniques for evaluating XAI systems: (1) using proxy, artificial tasks such as how well humans predict the AI's decision from the given explanations, and (2) using subjective measures of trust and preference as predictors of actual performance. The results of our experiments demonstrate that evaluations with proxy tasks did not predict the results of the evaluations with the actual decision-making tasks. Further, the subjective measures on evaluations with actual decision-making tasks did not predict the objective performance on those same tasks. Our results suggest that by employing misleading evaluation methods, our field may be inadvertently slowing its progress toward developing human+AI teams that can reliably perform better than humans or AIs alone.


title: GLIB: Efficient Exploration for Relational Model-Based Reinforcement Learning via Goal-Literal Babbling

abstract: We address the problem of efficient exploration for transition model learning in the relational model-based reinforcement learning setting without extrinsic goals or rewards. Inspired by human curiosity, we propose goal-literal babbling (GLIB), a simple and general method for exploration in such problems. GLIB samples relational conjunctive goals that can be understood as specific, targeted effects that the agent would like to achieve in the world, and plans to achieve these goals using the transition model being learned. We provide theoretical guarantees showing that exploration with GLIB will converge almost surely to the ground truth model. Experimentally, we find GLIB to strongly outperform existing methods in both prediction and planning on a range of tasks, encompassing standard PDDL and PPDDL planning benchmarks and a robotic manipulation task implemented in the PyBullet physics simulator. Video: https://youtu.be/F6lmrPT6TOY Code: https://git.io/JIsTB


title: Overcoming Noisy and Irrelevant Data in Federated Learning

abstract: Many image and vision applications require a large amount of data for model training. Collecting all such data at a central location can be challenging due to data privacy and communication bandwidth restrictions. Federated learning is an effective way of training a machine learning model in a distributed manner from local data collected by client devices, which does not require exchanging the raw data among clients. A challenge is that among the large variety of data collected at each client, it is likely that only a subset is relevant for a learning task while the rest of data has a negative impact on model training. Therefore, before starting the learning process, it is important to select the subset of data that is relevant to the given federated learning task. In this paper, we propose a method for distributedly selecting relevant data, where we use a benchmark model trained on a small benchmark dataset that is task-specific, to evaluate the relevance of individual data samples at each client and select the data with sufficiently high relevance. Then, each client only uses the selected subset of its data in the federated learning process. The effectiveness of our proposed approach is evaluated on multiple real-world image datasets in a simulated system with a large number of clients, showing up to $25\%$ improvement in model accuracy compared to training with all data.


title: How to Support Users in Understanding Intelligent Systems? Structuring the Discussion

abstract: The opaque nature of many intelligent systems violates established usability principles and thus presents a challenge for human-computer interaction. Research in the field therefore highlights the need for transparency, scrutability, intelligibility, interpretability and explainability, among others. While all of these terms carry a vision of supporting users in understanding intelligent systems, the underlying notions and assumptions about users and their interaction with the system often remain unclear. We review the literature in HCI through the lens of implied user questions to synthesise a conceptual framework integrating user mindsets, user involvement, and knowledge outcomes to reveal, differentiate and classify current notions in prior work. This framework aims to resolve conceptual ambiguity in the field and enables researchers to clarify their assumptions and become aware of those made in prior work. We thus hope to advance and structure the dialogue in the HCI research community on supporting users in understanding intelligent systems.


title: A B\'{e}koll\`{e}-Bonami Class of Weights for Certain Pseudoconvex Domains

abstract: We prove the weighted $L^p$ regularity of the ordinary Bergman projection on certain pseudoconvex domains where the weight belongs to an appropriate generalization of the B\'{e}koll\`{e}-Bonami class. The main tools used are estimates on the Bergman kernel obtained by McNeal and B\'{e}koll\`{e}'s original approach of proving a good-lambda inequality.


title: Regularity for a class of quasilinear degenerate parabolic equations in the Heisenberg group

abstract: We extend to the parabolic setting some of the ideas originated with Xiao Zhong's proof in \cite{Zhong} of the H\"older regularity of $p-$harmonic functions in the Heisenberg group $\Hn$. Given a number $p\ge 2$, in this paper we establish the $C^{\infty}$ smoothness of weak solutions of quasilinear pde's in $\Hn$ modelled on the equation $$\p_t u= \sum_{i=1}^{2n} X_i \bigg((1+|\nabla_0 u|^2)^{\frac{p-2}{2}} X_i u\bigg).$$


title: Cloud Atlas: High-precision HST/WFC3/IR Time-Resolved Observations of Directly-Imaged Exoplanet HD106906b

abstract: HD106906b is an ~11$M_{\mathrm{Jup}}$, ~15Myr old directly-imaged exoplanet orbiting at an extremely large distance from its host star. The wide separation (7.11 arcsec) between HD106906b and its host star greatly reduces the difficulty in direct-imaging observations, making it one of the most favorable directly-imaged exoplanets for detailed characterization. In this paper, we present HST/WFC3/IR time-resolved observations of HD106906b in the F127M, F139M, and F153M bands. We have achieved ~1% precision in the lightcurves in all three bands. The F127M lightcurve demonstrates marginally-detectable ($2.7\sigma$ significance) variability with a best-fitting period of 4 hr, while the lightcurves in the other two bands are consistent with flat lines. We construct primary-subtracted deep images and use these images to exclude additional companions to HD106906 that are more massive than 4$M_{\mathrm{Jup}}$ and locate at projected distances of more than ~500 au. We measure the astrometry of HD106906b in two HST/WFC3 epochs and achieve precisions better than 2.5 mas. The position angle and separation measurements do not deviate from those in the 2004 HST/ACS/HRC images for more than $1\sigma$ uncertainty. We provide the HST/WFC3 astrometric results for 25 background stars that can be used as reference sources in future precision astrometry studies. Our observations also provide the first 1.4-micron water band photometric measurement for HD106906b. HD106906b's spectral energy distribution and the best-fitting BT-Settl model have an inconsistency in the 1.4-micron water absorption band, which highlights the challenges in modeling atmospheres of young planetary-mass objects.


title: Inverse spectral results for non-abelian group actions

abstract: In this paper we will extend to non-abelian groups inverse spectral results, proved by us in an earlier paper, for compact abelian groups, i.e. tori. More precisely, Let $\mathsf G$ be a compact Lie group acting isometrically on a compact Riemannian manifold $X$. We will show that for the Schr\"odinger operator $-\hbar^2 \Delta+V$ with $V \in C^\infty(X)^{\mathsf G}$, the potential function $V$ is, in some interesting examples, determined by the $\mathsf G$-equivariant spectrum. The key ingredient in this proof is a generalized Legendrian relation between the Lagrangian manifolds $\mathrm{Graph}(dV)$ and $\mathrm{Graph}(dF)$, where $F$ is a spectral invariant defined on an open subset of the positive Weyl chamber.


title: NGC 6712: the variable star population of a tidally disrupted globular cluster

abstract: We present an analysis of $\emph{VI}$ CCD time series photometry of globular cluster NGC 6712. Our main goal is to study the variable star population as indicators of the cluster mean physical parameters. We employed the Fourier decomposition of RR Lyrae light curves to confirm that $[\rm Fe/H]_{UVES} = -1.0 \pm 0.05$ is a solid estimate. We estimated the reddening to the cluster as $E(B-V)$ = 0.35 $\pm$ 0.04 from the RRab stars colour curves. The distance to the cluster was estimated using three independent methods which yielded a weighted mean distance $\big < d \big>$ = 8.1 $\pm$ 0.2 kpc. The distribution of RRab and RRc stars on the HB shows a clear segregation around the first overtone red edge of the instability strip, which seems to be a common feature in OoI-type cluster with a very red horizontal branch. We carried out a membership analysis of 60,447 stars in our FoV using the data from $Gaia$-DR2 and found 1529 likely members; we possess the light curves of 1100 among the member stars. This allowed us to produce a clean colour-magnitude diagram, consistent with an age of 12 Gyrs, and enabled us to discover close unresolved contaminants for several variable stars. From the proper motion analysis we found evidence of non-member stars in the FoV of the cluster being tidally affected by the gravitational pull of the bulge of the Galaxy. We found that the RRab variable V6, shows a previously undetected Blazhko effect. Finally, we report sixteen new variables of the EW-type (9) and SR-type (7).


title: Model-Based Deep Learning for Reconstruction of Joint k-q Under-sampled High Resolution Diffusion MRI

abstract: We propose a model-based deep learning architecture for the reconstruction of highly accelerated diffusion magnetic resonance imaging (MRI) that enables high resolution imaging. The proposed reconstruction jointly recovers all the diffusion weighted images in a single step from a joint k-q under-sampled acquisition in a parallel MRI setting. We propose the novel use of a pre-trained denoiser as a regularizer in a model-based reconstruction for the recovery of highly under-sampled data. Specifically, we designed the denoiser based on a general diffusion MRI tissue microstructure model for multi-compartmental modeling. By using a wide range of biologically plausible parameter values for the multi-compartmental microstructure model, we simulated diffusion signal that spans the entire microstructure parameter space. A neural network was trained in an unsupervised manner using an autoencoder to learn the diffusion MRI signal subspace. We employed the autoencoder in a model-based reconstruction and show that the autoencoder provides a strong denoising prior to recover the q-space signal. We show reconstruction results on a simulated brain dataset that shows high acceleration capabilities of the proposed method.


title: Bayesian design for minimising prediction uncertainty in bivariate spatial responses with applications to air quality monitoring

abstract: Model-based geostatistical design involves the selection of locations to collect data to minimise an expected loss function over a set of all possible locations. The loss function is specified to reflect the aim of data collection, which, for geostatistical studies, could be to minimise the prediction uncertainty at unobserved locations. In this paper, we propose a new approach to design such studies via a loss function derived through considering the entropy about the model predictions and the parameters of the model. The approach also includes a multivariate extension to generalised linear spatial models, and thus can be used to design experiments with more than one response. Unfortunately, evaluating our proposed loss function is computationally expensive so we provide an approximation such that our approach can be adopted to design realistically sized geostatistical studies. This is demonstrated through a simulated study and through designing an air quality monitoring program in Queensland, Australia. The results show that our designs remain highly efficient in achieving each experimental objective individually, providing an ideal compromise between the two objectives. Accordingly, we advocate that our approach could be adopted more generally in model-based geostatistical design.


title: Positive factorizations of pseudoperiodic homeomorphisms

abstract: We generalize a classical result concerning smooth germs of surfaces, by proving that monodromies on links of isolated complex surface singularities associated with reduced holomorphic map germs admit a positive factorization. As a consequence of this and a topological characterization of these monodromies by Anne Pichon, we conclude that a pseudoperiodic homeomorphism on an oriented surface with boundary with positive fractional Dehn twist coefficients and screw numbers, admits a positive factorization. We use the main theorem to give a sufficiency criterion for certain pseudoperiodic homeomorphisms with negative screw numbers to admit a positive factorization.


title: Quantum phases of collective SU(3) spin systems with bipartite symmetry

abstract: We study a bipartite collective spin-$1$ model with exchange interaction between the spins. The bipartite nature of the model manifests itself by the spins being divided into two equal-sized subsystems; within each subsystem the spin-spin interactions are of the same strength, across the subsystems they are also equal, but the two coupling values within and across the subsystem are different. Such a set-up is inspired by recent experiments with ultracold atoms. Using the $\mathrm{SU}(3)$-symmetry of the exchange interaction and the permutation symmetry within the subsystems, we can employ representation theoretic methods to diagonalize the Hamiltonian of the system in the entire parameter space of the two coupling-strengths. These techniques then allow us to explicitly construct and explore the ground-state phase diagram. The phase diagram turns out to be rich containing both gapped and gapless phases. An interesting observation is that one of the five phases features a strong bipartite symmetry breaking, meaning that the two subsystems in the ground states are in different $\mathrm{SU}(3)$ representations.


title: Learning to adapt class-specific features across domains for semantic segmentation

abstract: Recent advances in unsupervised domain adaptation have shown the effectiveness of adversarial training to adapt features across domains, endowing neural networks with the capability of being tested on a target domain without requiring any training annotations in this domain. The great majority of existing domain adaptation models rely on image translation networks, which often contain a huge amount of domain-specific parameters. Additionally, the feature adaptation step often happens globally, at a coarse level, hindering its applicability to tasks such as semantic segmentation, where details are of crucial importance to provide sharp results. In this thesis, we present a novel architecture, which learns to adapt features across domains by taking into account per class information. To that aim, we design a conditional pixel-wise discriminator network, whose output is conditioned on the segmentation masks. Moreover, following recent advances in image translation, we adopt the recently introduced StarGAN architecture as image translation backbone, since it is able to perform translations across multiple domains by means of a single generator network. Preliminary results on a segmentation task designed to assess the effectiveness of the proposed approach highlight the potential of the model, improving upon strong baselines and alternative designs.


title: Arithmetic Combinatorics on Vinogradov systems

abstract: In this paper, we present a variant of the Balog-Szemer\'edi-Gowers theorem for the Vinogradov system. We then use our result to deduce a higher degree analogue of the sum-product phenomenon.


title: Analytic spread and integral closure of integrally decomposable modules

abstract: We relate the analytic spread of a module expressed as the direct sum of two submodules with the analytic spread of its components. We also study a class of submodules whose integral closure can be obtained by means of a simple computer algebra procedure once the integral closure of each row ideal is known. In particular, we analyze a class of modules, not necessarily of maximal rank, whose integral closure is determined by the family of Newton polyhedra of their row ideals.


title: Studies in Astronomical Time Series Analysis: VII. An Enquiry Concerning Non-Linearity, the RMS-Mean Flux Relation, and log-Normal Flux Distributions

abstract: A broad and widely used class of stationary, linear, additive time series models can have statistical properties which many authors have asserted imply that the underlying process must be non-linear, non-stationary, multiplicative, or inconsistent with shot noise. This result is demonstrated with exact and numerical evaluation of the model flux distribution function and dependence of flux standard deviation on mean flux (here and in the literature called the \emph{rms-flux relation}). These models can: (1) exhibit normal, log-normal or other flux distributions; (2) show linear or slightly non-linear rms-mean flux dependencies; as well as (3) match arbitrary second order statistics of the time series data. Accordingly the above assertions cannot be made on the basis of statistical time series analysis alone. Also discussed are ambiguities in the meaning of terms relevant to this study -- \emph{linear}, \emph{stationary} and \emph{multiplicative} -- and functions that can transform observed fluxes to a normal distribution as well or better than the logarithm.


title: Weight bounds for $(3,\gamma)$-hyperelliptic curves

abstract: {\it $(N,\gamma)$-hyperelliptic} semigroups were introduced by Fernando Torres to encapsulate the most salient properties of Weierstrass semigroups associated to totally-ramified points of $N$-fold covers of curves of genus $\gamma$. Torres characterized $(2,\gamma)$-hyperelliptic semigroups of maximal weight whenever their genus is large relative to $\gamma$. Here we do the same for $(3,\gamma)$-hyperelliptic semigroups, and we formulate a conjecture about the general case whenever $N \geq 3$ is prime.


title: Constraining the $^{22}$Ne($\alpha$,$\gamma$)$^{26}$Mg and $^{22}$Ne($\alpha$,n)$^{25}$Mg reaction rates using sub-Coulomb $\alpha$-transfer reactions

abstract: The $^{22}$Ne($\alpha$,$\gamma$)$^{26}$Mg and $^{22}$Ne($\alpha$,n)$^{25}$Mg reactions play an important role in astrophysics because they have significant influence on the neutron flux during the weak branch of the s-process. We constrain the astrophysical rates for these reactions by measuring partial $\alpha$-widths of resonances in $^{26}$Mg located in the Gamow window for the $^{22}$Ne+$\alpha$ capture. These resonances were populated using $^{22}$Ne($^6$Li,d)$^{26}$Mg and $^{22}$Ne($^7$Li,t)$^{26}$Mg reactions at energies near the Coulomb barrier. At these low energies $\alpha$-transfer reactions favor population of low spin states and the extracted partial $\alpha$-widths for the observed resonances exhibit only minor dependence on the model parameters. The astrophysical rates for both the $^{22}$Ne($\alpha$,$\gamma$)$^{26}$Mg and the $^{22}$Ne($\alpha$,n)$^{25}$Mg reactions are shown to be significantly different than the previously suggested values.


title: Deep Transformer Models for Time Series Forecasting: The Influenza Prevalence Case

abstract: In this paper, we present a new approach to time series forecasting. Time series data are prevalent in many scientific and engineering disciplines. Time series forecasting is a crucial task in modeling time series data, and is an important area of machine learning. In this work we developed a novel method that employs Transformer-based machine learning models to forecast time series data. This approach works by leveraging self-attention mechanisms to learn complex patterns and dynamics from time series data. Moreover, it is a generic framework and can be applied to univariate and multivariate time series data, as well as time series embeddings. Using influenza-like illness (ILI) forecasting as a case study, we show that the forecasting results produced by our approach are favorably comparable to the state-of-the-art.


title: Consistent Discretization of a Class of Predefined-Time Stable Systems

abstract: As the main contribution, this document provides a consistent discretization of a class of fixed-time stable systems, namely predefined-time stable systems. In the unperturbed case, the proposed approach allows obtaining not only a consistent but exact discretization of the considered class of predefined-time stable systems, whereas in the perturbed case, the consistent discretization preserves the predefined-time stability property. All the results are validated through simulations and compared with the conventional explicit Euler scheme, highlighting the advantages of this proposal.


title: Fano-Resonant, Asymmetric, Metamaterial-Assisted Tweezers for Single Nanoparticle Trapping

abstract: Plasmonic nanostructures can overcome Abbe's diffraction limit to generate strong gradient fields, enabling efficient optical trapping of nano-sized particles. However, it remains challenging to achieve stable trapping with low incident laser intensity. Here, we demonstrate a Fano resonance-assisted plasmonic optical tweezers (FAPOT), for single nanoparticle trapping in an array of asymmetrical split nano-apertures, milled on a 50 nm gold thin film. Stable trapping is achieved by tuning the trapping wavelength and varying the incident trapping laser intensity. A very large normalized trap stiffness of 8.65 fN/nm/mW for 20 nm polystyrene particles at a near-resonance trapping wavelength of 930 nm was achieved. We show that trap stiffness on resonance is enhanced by a factor of 63 compared to off-resonance conditions. This can be attributed to the ultra-small mode volume, which enables large near-field strengths and a cavity Purcell effect contribution. These results should facilitate strong trapping with low incident trapping laser intensity, thereby providing new options for studying transition paths of single molecules, such as proteins, DNA, or viruses.


title: The theory of direct laser excitation of nuclear transitions

abstract: A comprehensive theoretical study of direct laser excitation of a nuclear state based on the density matrix formalism is presented. The nuclear clock isomer $^{229\text{m}}$Th is discussed in detail, as it could allow for direct laser excitation using existing technology and provides the motivation for this work. The optical Bloch equations are derived for the simplest case of a pure nuclear two-level system and for the more complex cases taking into account the presence of magnetic sub-states, hyperfine-structure and Zeeman splitting in external fields. Nuclear level splitting for free atoms and ions as well as for nuclei in a solid-state environment is discussed individually. Based on the obtained equations, nuclear population transfer in the low-saturation limit is reviewed. Further, nuclear Rabi oscillations, power broadening and nuclear two-photon excitation are considered. Finally, the theory is applied to the special cases of $^{229\text{m}}$Th and $^{235\text{m}}$U, being the nuclear excited states of lowest known excitation energies. The paper aims to be a didactic review with many calculations given explicitly.


title: Bulk medium evolution has considerable effects on jet observables

abstract: We consider the case, in QCD, of a single jet propagating within a strongly interacting fluid, of finite extent. Interactions lead to the appearance of a source of energy-momentum within the fluid. The remnant jet that escapes the container is analyzed along with portions of the medium excited by the jet. We study the effect of a static versus a semi-realistic expanding medium, with jets traveling inward versus outward. We consider the medium response via recoils in partonic scatterings based on a weakly-coupled description and its combination with hydrodynamical medium response based on a strongly-coupled description, followed by incorporation into a jet. The effects of these limits on the reconstructed energy, momentum and mass of the jet, as a function of the angle away from the original parton direction are studied. It is demonstrated that different flow velocity configurations in the medium produce considerable differences in jet observables. This work highlights the importance of accurate dynamical modeling of the soft medium as a foundation on which to calculate jet modification, and casts skepticism on results obtained without such modeling.


title: FsNet: Feature Selection Network on High-dimensional Biological Data

abstract: Biological data including gene expression data are generally high-dimensional and require efficient, generalizable, and scalable machine-learning methods to discover their complex nonlinear patterns. The recent advances in machine learning can be attributed to deep neural networks (DNNs), which excel in various tasks in terms of computer vision and natural language processing. However, standard DNNs are not appropriate for high-dimensional datasets generated in biology because they have many parameters, which in turn require many samples. In this paper, we propose a DNN-based, nonlinear feature selection method, called the feature selection network (FsNet), for high-dimensional and small number of sample data. Specifically, FsNet comprises a selection layer that selects features and a reconstruction layer that stabilizes the training. Because a large number of parameters in the selection and reconstruction layers can easily result in overfitting under a limited number of samples, we use two tiny networks to predict the large, virtual weight matrices of the selection and reconstruction layers. Experimental results on several real-world, high-dimensional biological datasets demonstrate the efficacy of the proposed method.


title: Lattice results on dibaryons and baryon--baryon interactions

abstract: We present the recent study on dibaryons at the almost physical pion mass in lattice QCD by the HAL QCD potential method.


title: Optimal fermionic swap networks for Hubbard models

abstract: We propose an efficient variation of the fermionic swap network scheme used to efficiently simulate n-dimensional Fermi-Hubbard-model Hamiltonians encoded using the Jordan-Wigner transform. For the two-dimensional versions, we show that our choices minimize swap depth and number of Hamiltonian interaction layers. The proofs, along with the choice of swap network, rely on isoperimetric inequality results from the combinatorics literature, and are closely related to graph bandwidth problems. The machinery has the potential to be extended to maximize swap network efficiency for other types of lattices.


title: A population of dust-enshrouded objects orbiting the Galactic black hole

abstract: The central 0.1 parsecs of the Milky Way host a supermassive black hole identified with the position of the radio and infrared source Sagittarius A*, a cluster of young, massive stars (the S stars) and various gaseous features. Recently, two unusual objects have been found to be closely orbiting Sagittarius A*: the so-called G sources, G1 and G2. These objects are unresolved (having a size of the order of 100 astronomical units, except at periapse, where the tidal interaction with the black hole stretches them along the orbit) and they show both thermal dust emission and line emission from ionized gas. G1 and G2 have generated attention because they appear to be tidally interacting with the supermassive Galactic black hole, possibly enhancing its accretion activity. No broad consensus has yet been reached concerning their nature: the G objects show the characteristics of gas and dust clouds but display the dynamical properties of stellar-mass objects. Here we report observations of four additional G objects, all lying within 0.04 parsecs of the black hole and forming a class that is probably unique to this environment. The widely varying orbits derived for the six G objects demonstrate that they were commonly but separately formed.


title: On the Polarization of Ligands by Proteins

abstract: Although ligand-binding sites in many proteins contain a high number density of charged side chains that can polarize small organic molecules and influence binding, the magnitude of this effect has not been studied in many systems. Here, we use a quantum mechanics/molecular mechanics (QM/MM) approach in which the ligand is the QM region to compute the ligand polarization energy of 286 protein-ligand complexes from the PDBBind Core Set (release 2016). We observe that the ligand polarization energy is linearly correlated with the magnitude of the electric field acting on the ligand, the magnitude of the induced dipole moment, and the classical polarization energy. The influence of protein and cation charges on the ligand polarization diminishes with the distance and is below 2 kcal/mol at 9 $\unicode{x212B}$ and 1 kcal/mol at 12 $\unicode{x212B}$. Considering both polarization and solvation appears essential to computing negative binding energies in some crystallographic complexes. Solvation, but not polarization, is essential for achieving moderate correlation with experimental binding free energies.


title: The Reciprocal Bayesian LASSO

abstract: A reciprocal LASSO (rLASSO) regularization employs a decreasing penalty function as opposed to conventional penalization approaches that use increasing penalties on the coefficients, leading to stronger parsimony and superior model selection relative to traditional shrinkage methods. Here we consider a fully Bayesian formulation of the rLASSO problem, which is based on the observation that the rLASSO estimate for linear regression parameters can be interpreted as a Bayesian posterior mode estimate when the regression parameters are assigned independent inverse Laplace priors. Bayesian inference from this posterior is possible using an expanded hierarchy motivated by a scale mixture of double Pareto or truncated normal distributions. On simulated and real datasets, we show that the Bayesian formulation outperforms its classical cousin in estimation, prediction, and variable selection across a wide range of scenarios while offering the advantage of posterior inference. Finally, we discuss other variants of this new approach and provide a unified framework for variable selection using flexible reciprocal penalties. All methods described in this paper are publicly available as an R package at: https://github.com/himelmallick/BayesRecipe.


title: A Deep Learning Approach to Behavior-Based Learner Modeling

abstract: The increasing popularity of e-learning has created demand for improving online education through techniques such as predictive analytics and content recommendations. In this paper, we study learner outcome predictions, i.e., predictions of how they will perform at the end of a course. We propose a novel Two Branch Decision Network for performance prediction that incorporates two important factors: how learners progress through the course and how the content progresses through the course. We combine clickstream features which log every action the learner takes while learning, and textual features which are generated through pre-trained GloVe word embeddings. To assess the performance of our proposed network, we collect data from a short online course designed for corporate training and evaluate both neural network and non-neural network based algorithms on it. Our proposed algorithm achieves 95.7% accuracy and 0.958 AUC score, which outperforms all other models. The results also indicate the combination of behavior features and text features are more predictive than behavior features only and neural network models are powerful in capturing the joint relationship between user behavior and course content.


title: Leveraging Neighborhood Summaries for Efficient RDF Queries on RDBMS

abstract: Using structural informations to summarize graph-structured RDF data is helpful in tackling query performance issues. However, leveraging structural indexes needs to revise or even redesign the internal of RDF systems. Given an RDF dataset that have already been bulk loaded into a relational RDF system, we aim at improving the query performance on such systems. We do so by summarizing neighborhood structures and encoding them into triples which can be managed along side the exist instance data. At query time, we optimally select the effective structural patterns, and adding these patterns to the existing queries to gain an improved query performance. Empirical evaluations shown the effectiveness of our method.


title: Maximizing the $p$-th moment of exit time of planar Brownian motion from a given domain

abstract: In this paper we address the question of finding the point which maximizes the $p$-th moment of the exit time of planar Brownian motion from a given domain. We present a geometrical method of excluding parts of the domain from consideration which makes use of a coupling argument and the conformal invariance of Brownian motion. In many cases the maximizing point can be localized to a relatively small region. Several illustrative examples are presented.


title: On the probability of fast exits and long stays of planar Brownian motion in simply connected domains

abstract: Let $T^D$ denote the first exit time of a planar Brownian motion from a domain $D$. Given two simply connected planar domains $U,W \neq \SC$ containing $0$, we investigate the cases in which we are more likely to have fast exits (meaning for instance ${\bf P}(T^U<t) > {\bf P}(T^W<t)$ for $t$ small) from $U$ than from $W$, or long stays (meaning ${\bf P}(T^U>t) > {\bf P}(T^W>t)$ for $t$ large). We prove several results on these questions. In particular, we show that the primary factor in the probability of fast exits is the proximity of the boundary to the origin, while for long stays an important factor is the moments of the exit time. The complex analytic theory that motivated our inquiry is also discussed.


title: Nonlinear control in the nematode C. elegans

abstract: Recent whole-brain calcium imaging recordings of the nematode C. elegans have demonstrated that neural activity is dominated by dynamics on a low-dimensional manifold that can be clustered according to behavioral states. Despite progress in modeling the dynamics with linear or locally linear models, it remains unclear how a single network of neurons can produce the observed features. In particular, there are multiple clusters, or fixed points, observed in the data which cannot be characterized by a single linear model. We propose a nonlinear control model which is global and parameterized by only four free parameters that match the features displayed by the low-dimensional C. elegans neural activity. In addition to reproducing the average probability distribution of the data, long and short time-scale changes in transition statistics can be characterized via changes in a single parameter. Some of these macro-scale transitions have experimental correlates to single neuro-modulators that seem to act as biological controls, allowing this model to generate testable hypotheses about the effect of these neuro-modulators on the global dynamics. The theory provides an elegant characterization of the neuron population dynamics in C. elegans. Moreover, the mathematical structure of the nonlinear control framework provides a paradigm that can be generalized to more complex systems with an arbitrary number of behavioral states.


title: Applying Recent Innovations from NLP to MOOC Student Course Trajectory Modeling

abstract: This paper presents several strategies that can improve neural network-based predictive methods for MOOC student course trajectory modeling, applying multiple ideas previously applied to tackle NLP (Natural Language Processing) tasks. In particular, this paper investigates LSTM networks enhanced with two forms of regularization, along with the more recently introduced Transformer architecture.


title: Concentration phenomena to a higher order Liouville equation

abstract: We study blow-up and quantization phenomena for a sequence of solutions $(u_k)$ to the prescribed $Q$-curvature problem $$ (-\Delta)^nu_k= Q_ke^{2nu_k}\quad \text{in }\Omega\subset\mathbb{R}^{2n},\quad \int_{\Omega}e^{2nu_k}dx\leq C,$$ under natural assumptions on $Q_k$. It is well-known that, up to a subsequence, either $(u_k)$ is bounded in a suitable norm, or there exists $\beta_k\to\infty$ such that $ u_k=\beta_k(\varphi+o(1))$ in $\Omega\setminus (S_1\cup S_\varphi)$ for some non-trivial non-positive $n$-harmonic function $\varphi$ and for a finite set $S_1$, where $S_\varphi$ is the zero set of $\varphi$. We prove quantization of the total curvature $\int_{\tilde\Omega}Q_ke^{2nu_k}dx$ on the region $\tilde\Omega\Subset(\Omega\setminus S_\varphi)$. We also consider a non-local case in dimension three.


title: Numerical Abstract Persuasion Argumentation for Expressing Concurrent Multi-Agent Negotiations

abstract: A negotiation process by 2 agents e1 and e2 can be interleaved by another negotiation process between, say, e1 and e3. The interleaving may alter the resource allocation assumed at the inception of the first negotiation process. Existing proposals for argumentation-based negotiations have focused primarily on two-agent bilateral negotiations, but scarcely on the concurrency of multi-agent negotiations. To fill the gap, we present a novel argumentation theory, basing its development on abstract persuasion argumentation (which is an abstract argumentation formalism with a dynamic relation). Incorporating into it numerical information and a mechanism of handshakes among members of the dynamic relation, we show that the extended theory adapts well to concurrent multi-agent negotiations over scarce resources.


title: Geometric Conditions for the Discrepant Posterior Phenomenon and Connections to Simpson's Paradox

abstract: The discrepant posterior phenomenon (DPP) is a counter-intuitive phenomenon that can frequently occur in a Bayesian analysis of multivariate parameters. It refers to the phenomenon that a parameter estimate based on a posterior is more extreme than both of those inferred based on either the prior or the likelihood alone. Inferential claims that exhibit DPP defy the common intuition that the posterior is a prior-data compromise, and the phenomenon can be surprisingly ubiquitous in well-behaved Bayesian models. In this paper we revisit this phenomenon and, using point estimation as an example, derive conditions under which the DPP occurs in Bayesian models with exponential quadratic likelihoods and conjugate multivariate Gaussian priors. The family of exponential quadratic likelihood models includes Gaussian models and those models with local asymptotic normality property. We provide an intuitive geometric interpretation of the phenomenon and show that there exists a nontrivial space of marginal directions such that the DPP occurs. We further relate the phenomenon to the Simpson's paradox and discover their deep-rooted connection that is associated with marginalization. We also draw connections with Bayesian computational algorithms when difficult geometry exists. Our discovery demonstrates that DPP is more prevalent than previously understood and anticipated. Theoretical results are complemented by numerical illustrations. Scenarios covered in this study have implications for parameterization, sensitivity analysis, and prior choice for Bayesian modeling.


title: Relaxation of ferromagnetic domains in a disordered lattice in 2D

abstract: We investigate the relaxation process of ferromagnetic domains in 2D subjected to the influence of both, static disorder of variable strength and weak interactions. The domains are represented by a two species bosonic mixture of $^{87}$Rb ultracold atoms, such that initially each specie lies on left and right halves of a square lattice. The dynamics of the double domain is followed by describing the two-component superfluid, at mean field level, through the time dependent Gross-Pitaevskii coupled equations, considering values of the intra and inter-species interaction, reachable in current experimental setups, that guaranty miscibility of the components. A robust analysis for several values inter-species interaction leads us to conclude that the presence of structural disorder leads to slowdown the relaxation process of the initial ferromagnetic order. As shown by our numerical experiments, magnetization is maintained up to 60 percent of its initial value for the largest disorder amplitude.


title: Planar Heyting Algebras for Children 2: Local Operators, J-Operators, and Slashings

abstract: Choose a topos $E$. There are several different "notions of sheafness" on $E$. How do we visualize them? Let's refer to the classifier object of $E$ as $\Omega$, and to its Heyting Algebra of truth-values, $Sub(1_E)$, as $H$; we will sometimes call $H$ the "logic" of the topos. There is a well-known way of representing notions of sheafness as morphisms $j:\Omega\to \Omega$, but these `$j$'s yield big diagrams when we draw them explicitly; here we will see a way to represent these `$j$'s as maps $J:H\to H$ in a way that is much more manageable. In the previous paper of this series we showed how certain toy models of Heyting Algebras, called "ZHAs", can be used to develop visual intuition for how Heyting Algebras and Intuitionistic Propositional Logic work; here we will extend that to sheaves. The full idea is this: notions of sheafness correspond to local operators and vice-versa; local operators correspond to J-operators and vice-versa; if our Heyting Algebra $H$ is a ZHA then J-operators correspond to slashings on $H$, and vice-versa; slashings on $H$ correspond to "sets of question marks" and vice-versa, and each set of question marks induces a notion of erasing and reconstructing, which induces a sheaf. Also, every ZHA $H$ corresponds to an (acyclic) 2-column graph, and vice-versa, and for any two-column graph $(P,A)$ the logic of the topos $\mathbf{Set}^{(P,A)}$ is exactly the ZHA $H$ associated to $(P,A)$.


title: Cobordism invariance of topological edge-following states

abstract: We prove that a spectral gap-filling phenomenon occurs whenever a Hamiltonian operator encounters a coarse index obstruction upon compression to a domain with boundary. Furthermore, the gap-filling spectra contribute to quantised current channels, which follow and are localised at the possibly complicated boundary. This index obstruction is shown to be insensitive to deformations of the domain boundary, so the phenomenon is generic for magnetic Laplacians modelling quantum Hall systems and Chern topological insulators. A key construction is a quasi-equivariant version of Roe's algebra of locally compact finite propagation operators.

categories: math-ph cond-mat.mes-hall math.KT math.MP math.OA


title: The Radial Acceleration Relation in CLASH Galaxy Clusters

abstract: The radial acceleration relation (RAR) in galaxies describes a tight empirical scaling law between the total acceleration $g_\mathrm{tot}(r)=GM_\mathrm{tot}(<r)/r^2$ observed in galaxies and that expected from their baryonic mass $g_\mathrm{bar}(r)=GM_\mathrm{bar}(<r)/r^2$, with a characteristic acceleration scale of $g_\dagger\simeq 1.2\times 10^{-10}$ms$^{-2}$. Here, we examine if such a correlation exists in galaxy clusters using weak-lensing, strong-lensing, and X-ray data sets available for 20 high-mass clusters targeted by the CLASH survey. By combining our CLASH data with stellar mass estimates for the brightest cluster galaxies (BCGs) and accounting for the stellar baryonic component in clusters, we determine, for the first time, an RAR on BCG--cluster scales. The resulting RAR is well described by a tight power-law relation, $g_\mathrm{tot}\propto g_\mathrm{bar}^{0.51^{+0.04}_{-0.05}}$, with lognormal intrinsic scatter of $14.7^{+2.9}_{-2.8}\%$. The slope is consistent with the low acceleration limit of the RAR in galaxies, $g_\mathrm{tot}=\sqrt{g_\dagger\,g_\mathrm{bar}}$, whereas the intercept implies a much higher acceleration scale of $g_\ddagger = (2.02\pm0.11)\times 10^{-9}$ms$^{-2}$, indicating that there is no universal RAR that holds on all scales from galaxies to clusters. We find that the observed RAR in CLASH clusters is consistent with predictions from a semi-analytical model developed in the standard $\Lambda$CDM framework. Our results also predict the presence of a baryonic Faber--Jackson relation ($\sigma_v^4\propto M_\mathrm{bar}$) on cluster scales.


title: Reconfigurable multiport switch in coupled mode devices

abstract: We consider a coupled mode system where the effective propagation constants of localized modes are amenable to modulation. Starting from an unmodulated system where power transfer is heavily suppressed, we demonstrate that small, periodic modulation of the propagation constants enhance power transfer using a slowly varying envelope approximation for the field mode amplitudes. We calculate an approximate modulation frequency enabling complete transfer between otherwise negligibly coupled elements. The ability to control these modulations by electrical or thermal effects allows for reconfigurable multiport switching. We use an array of coupled silica waveguides and the thermo-optic effect to test our predictions in the telecom C-band. However, this requires a refractive index modulation with period of the order $~10^{-3}~\mathrm{m}$ and yields total power transfer with a propagation distance of the order of $10^{-1}~\mathrm{m}$, which might make it unattractive for integrated photonic applications. Nevertheless, our results are valid for devices described by an equivalent coupled mode matrix for space or time propagation; for example, arrays of microring or terahertz resonators, microwave cavities, radio frequency antennas, or RLC circuits.


title: KMTNet Nearby Galaxy Survey II. Searching for Dwarf Galaxies in Deep and Wide-field Images of the NGC 1291 system

abstract: We present newly discovered dwarf galaxy candidates in deep and wide-field images of NGC 1291 obtained with the Korea Microlensing Telescope Network. We identify 15 dwarf galaxy candidates by visual inspection. Using imaging simulations, we demonstrate that the completeness rate of our detection is greater than 70% for the central surface brightness value of ${\mu}_{0,R}\lesssim26$ mag arcsec$^{-2}$ and for magnitudes $M_R\lesssim-10$ mag. The structural and photometric properties of the dwarf galaxy candidates appear to be broadly consistent with those of ordinary dwarf galaxies in nearby groups and clusters, with ${\mu}_{0,R}\sim$ 22.5 to 26.5 mag arcsec$^{-2}$ and effective radii of 200 pc to 1 kpc. The dwarf galaxy candidates show a concentration towards NGC 1291 and tend to be redder the closer they are to the center, possibly indicating that they are associated with NGC 1291. The dwarf candidates presented in this paper appear to be bluer than those in denser environments, revealing that the quenching of star formation in dwarf galaxies is susceptible to the environment, while the morphology shaping is not.


title: Demonstrating a Continuous Set of Two-qubit Gates for Near-term Quantum Algorithms

abstract: Quantum algorithms offer a dramatic speedup for computational problems in machine learning, material science, and chemistry. However, any near-term realizations of these algorithms will need to be heavily optimized to fit within the finite resources offered by existing noisy quantum hardware. Here, taking advantage of the strong adjustable coupling of gmon qubits, we demonstrate a continuous two-qubit gate set that can provide a 3x reduction in circuit depth as compared to a standard decomposition. We implement two gate families: an iSWAP-like gate to attain an arbitrary swap angle, $\theta$, and a CPHASE gate that generates an arbitrary conditional phase, $\phi$. Using one of each of these gates, we can perform an arbitrary two-qubit gate within the excitation-preserving subspace allowing for a complete implementation of the so-called Fermionic Simulation, or fSim, gate set. We benchmark the fidelity of the iSWAP-like and CPHASE gate families as well as 525 other fSim gates spread evenly across the entire fSim($\theta$, $\phi$) parameter space achieving purity-limited average two-qubit Pauli error of $3.8 \times 10^{-3}$ per fSim gate.


title: Minimizing the Discounted Probability of Exponential Parisian Ruin via Reinsurance

abstract: We study the problem of minimizing the discounted probability of exponential Parisian ruin, that is, the discounted probability that an insurer's surplus exhibits an excursion below zero in excess of an exponentially distributed clock. The insurer controls its surplus via reinsurance priced according to the mean-variance premium principle, as in Liang, Liang, and Young (2019). We consider the classical risk model and apply stochastic Perron's method, as introduced by Bayraktar and Sirbu (2012,2013,2014), to show that the minimum discounted probability of exponential Parisian ruin is the unique viscosity solution of its Hamilton-Jacobi-Bellman equation with boundary conditions at $\pm \infty$. A major difficulty in proving the comparison principle arises from the discontinuity of the Hamiltonian.


title: Target-Embedding Autoencoders for Supervised Representation Learning

abstract: Autoencoder-based learning has emerged as a staple for disciplining representations in unsupervised and semi-supervised settings. This paper analyzes a framework for improving generalization in a purely supervised setting, where the target space is high-dimensional. We motivate and formalize the general framework of target-embedding autoencoders (TEA) for supervised prediction, learning intermediate latent representations jointly optimized to be both predictable from features as well as predictive of targets---encoding the prior that variations in targets are driven by a compact set of underlying factors. As our theoretical contribution, we provide a guarantee of generalization for linear TEAs by demonstrating uniform stability, interpreting the benefit of the auxiliary reconstruction task as a form of regularization. As our empirical contribution, we extend validation of this approach beyond existing static classification applications to multivariate sequence forecasting, verifying their advantage on both linear and nonlinear recurrent architectures---thereby underscoring the further generality of this framework beyond feedforward instantiations.


title: Unsupervised learning of control signals and their encodings in $\textit{C. elegans}$ whole-brain recordings

abstract: Recent whole brain imaging experiments on $\textit{C. elegans}$ has revealed that the neural population dynamics encode motor commands and stereotyped transitions between behaviors on low dimensional manifolds. Efforts to characterize the dynamics on this manifold have used piecewise linear models to describe the entire state space, but it is unknown how a single, global dynamical model can generate the observed dynamics. Here, we propose a control framework to achieve such a global model of the dynamics, whereby underlying linear dynamics is actuated by sparse control signals. This method learns the control signals in an unsupervised way from data, then uses $\textit{ Dynamic Mode Decomposition with control}$ (DMDc) to create the first global, linear dynamical system that can reconstruct whole-brain imaging data. These control signals are shown to be implicated in transitions between behaviors. In addition, we analyze the time-delay encoding of these control signals, showing that these transitions can be predicted from neurons previously implicated in behavioral transitions, but also additional neurons previously unidentified. Moreover, our decomposition method allows one to understand the observed nonlinear global dynamics instead as linear dynamics with control. The proposed mathematical framework is generic and can be generalized to other neurosensory systems, potentially revealing transitions and their encodings in a completely unsupervised way.


title: On a conjectural congruence of Guo

abstract: Let $p\equiv3\pmod{4}$ be a prime and $r$ a positive integer. We show that $$ \prod_{k=1}^{(p^{2r}-1)/2}\frac{4k-1}{4k+1}\equiv1\pmod{p^2}. $$ This confirms a recent conjecture of Guo.


title: Fermi liquid theory for nonlinear transport through a multilevel Anderson impurity

abstract: We present a microscopic Fermi-liquid view on the low-energy transport through an Anderson impurity with $N$ discrete levels, at arbitrary electron filling $N_d$. It is applied to nonequilibrium current fluctuations, for which the two-quasiparticle collision integral and the three-body correlations that determine the quasiparticle energy shift play important roles. Using the numerical renormalization group up to $N=6$, we find that for strong interactions the three-body fluctuations are determined by a single parameter other than the Kondo energy scale in a wide filling range $1 \lesssim N_d \lesssim N-1$. It significantly affects the current noise for $N>2$ and the behavior of noise in magnetic fields.


title: Investigating naturalistic hand movements by behavior mining in long-term video and neural recordings

abstract: Recent technological advances in brain recording and artificial intelligence are propelling a new paradigm in neuroscience beyond the traditional controlled experiment. Rather than focusing on cued, repeated trials, naturalistic neuroscience studies neural processes underlying spontaneous behaviors performed in unconstrained settings. However, analyzing such unstructured data lacking a priori experimental design remains a significant challenge, especially when the data is multi-modal and long-term. Here we describe an automated approach for analyzing simultaneously recorded long-term, naturalistic electrocorticography (ECoG) and naturalistic behavior video data. We take a behavior-first approach to analyzing the long-term recordings. Using a combination of computer vision, discrete latent-variable modeling, and string pattern-matching on the behavioral video data, we find and annotate spontaneous human upper-limb movement events. We show results from our approach applied to data collected for 12 human subjects over 7--9 days for each subject. Our pipeline discovers and annotates over 40,000 instances of naturalistic human upper-limb movement events in the behavioral videos. Analysis of the simultaneously recorded brain data reveals neural signatures of movement that corroborate prior findings from traditional controlled experiments. We also prototype a decoder for a movement initiation detection task to demonstrate the efficacy of our pipeline as a source of training data for brain-computer interfacing applications. Our work addresses the unique data analysis challenges in studying naturalistic human behaviors, and contributes methods that may generalize to other neural recording modalities beyond ECoG. We publicly release our curated dataset, providing a resource to study naturalistic neural and behavioral variability at a scale not previously available.