title: The best of both worlds: Combining LOFAR and Apertif to derive resolved radio spectral index images

abstract: Supermassive black holes at the centres of galaxies can cycle through periods of activity and quiescence. Characterising the duty cycle of active galactic nuclei is crucial for understanding the impact of the energy they release on the host galaxy. For radio AGN, this can be done by identifying dying (remnant) and restarted radio galaxies from their radio spectral properties. Using the combination of images at 1400 MHz produced by Apertif, the new phased-array feed receiver installed on the Westerbork Synthesis Radio Telescope, and images at 150 MHz provided by LOFAR, we have derived resolved spectral index images (at a resolution of ~15 arcsec) for all the sources within ~6 deg^2 area of the Lockman Hole region. We were able to select 15 extended radio sources with emission (partly or entirely) characterised by extremely steep spectral indices (steeper than 1.2). These objects represent radio sources in the remnant or the restarted phases of their life cycle. Our findings suggest this cycle to be relatively fast. They also show a variety of properties relevant for modelling the evolution of radio galaxies. For example, the restarted activity can occur while the remnant structure from a previous phase of activity is still visible. This provides constraints on the duration of the 'off' (dying) phase. In extended remnants with ultra-steep spectra at low frequencies, the activity likely stopped a few hundred megayears ago, and they correspond to the older tail of the age distribution of radio galaxies, in agreement with simulations of radio source evolution. We find remnant radio sources with a variety of structures (from double-lobed to amorphous), suggesting different types of progenitors. The present work sets the stage for exploiting low-frequency spectral index studies of extended sources by taking advantage of the large areas common to the LOFAR and the Apertif surveys.


title: Chiral gauge leptoquark mass limits and branching ratios of $ K_L^0, B^0, B_s \to l^+_i l^-_j $ decays with account of the general fermion mixing in leptoquark currents

abstract: The contributions of the chiral gauge leptoquarks $V^{L,R}$ induced by the chiral four color quark-lepton symmetry to the branching ratios of $ K_L^0, B^0, B_s \to l_1 \, l_2 $ decays are calculated and analysed using the general parametrizations of the fermion mixing matrices in the leptoguark currents. From the current experimental data on these decays under assumption $ m_{V^L} \ll m_{V^R} $ the lower mass limit $ m_{V^L} \cos{\gamma_L} > 5.68 \,\, \mbox{TeV}$ is found, which in particular case of equal gauge coupling constants gives $ m_{V^L} > 8.03 \,\, \mbox{TeV} $. The branching ratios of the decays under consideration predicted by the chiral gauge leptoquarks are calculated and analysed in dependence on the leptoquark masses and the mixing parameters. It is shown that in consistency with the current experimental data these branching ratios for $ B_s, B^0 \to \mu e $ decays can be close to their experimental limits and those for $ B_s, B^0 \to \tau e, \tau \mu $ decays can be of order of~$10^{-7}$. The calculated branching ratios will be useful in the further experimental searches for these decays.


title: Diffusion with Local Resetting and Exclusion

abstract: Stochastic resetting models diverse phenomena across numerous scientific disciplines. Current understanding stems from the renewal framework, which relates systems subject to global resetting to their non-resetting counterparts. Yet, in interacting many-body systems, even the simplest scenarios involving resetting give rise to the notion of local resetting, whose analysis falls outside the scope of the renewal approach. A prime example is that of diffusing particles with excluded volume interactions that independently attempt to reset their position to the origin of a 1D lattice. With renewal rendered ineffective, we instead employ a mean-field approach whose validity is corroborated via extensive numerical simulations. The emerging picture sheds first light on the non-trivial interplay between interactions and resetting in many-body systems.


title: Opportunities and Challenges for Circuit Board Level Hardware Description Languages

abstract: Board-level hardware description languages (HDLs) are one approach to increasing automation and raising the level of abstraction for designing electronics. These systems borrow programming languages concepts like generators and type systems, but also must be designed with human factors in mind to serve existing hardware engineers. In this work, we look at one recent prototype system, and discuss open questions spanning from fundamental models through usable interfaces.


title: Dialog Simulation with Realistic Variations for Training Goal-Oriented Conversational Systems

abstract: Goal-oriented dialog systems enable users to complete specific goals like requesting information about a movie or booking a ticket. Typically the dialog system pipeline contains multiple ML models, including natural language understanding, state tracking and action prediction (policy learning). These models are trained through a combination of supervised or reinforcement learning methods and therefore require collection of labeled domain specific datasets. However, collecting annotated datasets with language and dialog-flow variations is expensive, time-consuming and scales poorly due to human involvement. In this paper, we propose an approach for automatically creating a large corpus of annotated dialogs from a few thoroughly annotated sample dialogs and the dialog schema. Our approach includes a novel goal-sampling technique for sampling plausible user goals and a dialog simulation technique that uses heuristic interplay between the user and the system (Alexa), where the user tries to achieve the sampled goal. We validate our approach by generating data and training three different downstream conversational ML models. We achieve 18 ? 50% relative accuracy improvements on a held-out test set compared to a baseline dialog generation approach that only samples natural language and entity value variations from existing catalogs but does not generate any novel dialog flow variations. We also qualitatively establish that the proposed approach is better than the baseline. Moreover, several different conversational experiences have been built using this method, which enables customers to have a wide variety of conversations with Alexa.


title: On the precision of full-spectrum fitting of simple stellar populations. II. The dependence on star cluster mass in the wavelength range 0.3-5.0 $\mu$m

abstract: In this second paper of a series on the accuracy and precision of the determination of age and metallicity of simple stellar populations (SSPs) by means of the full spectrum fitting technique, we study the influence of star cluster mass through stochastic fluctuations of the number of stars near the top of the stellar mass function, which dominate the flux in certain wavelength regimes depending on the age. We consider SSP models based on the Padova isochrones, spanning the age range 7.0 $\leq$ log (age/yr) $\leq$ 10.1. Simulated spectra of star clusters in the mass range $10^4 \leq M/M_{\odot} < 10^6$ are compared with SSP model spectra to determine best-fit ages and metallicities using a full-spectrum fitting routine in four wavelength regimes: the blue optical (0.35-0.70 $\mu$m), the red optical (0.6-1.0 $\mu$m), the near-IR (1.0-2.5 $\mu$m), and the mid-IR (2.5-5.0 $\mu$m). We compare the power of each wavelength regime in terms of both the overall precision of age and metallicity determination, and of its dependence on cluster mass. We also study the relevance of spectral resolution in this context by utilizing two different spectral libraries (BaSeL and BT-Settl). We highlight the power of the mid-IR regime in terms of identifying young massive clusters in dusty star forming regions in external galaxies. The spectra of the simulated star clusters and SSPs are made available online to enable follow-up studies by the community.


title: Melvin's "magnetic universe", the role of the magnetic tension and the implications for gravitational collapse

abstract: Electromagnetism appears to have the potential to alter the fate of relativistic gravitational implosion. Coulomb forces, for example, can in principle prevent the ultimate collapse of charged matter to a singularity. Also, half a century ago, Mael Melvin noted that magnetic forcelines do not collapse under their own gravity, no matter how strong the latter may be. Instead, magnetic fields seemed capable of stabilising themselves against gravitational self-contraction. Intrigued, Melvin wondered whether magnetism could also support against the realistic collapse of charged matter. Here, we look back into these issues by means of two complementary scenarios. The first reconsiders Melvin's highly idealised "magnetic universe" and reaches the same conclusion, namely that magnetic fields do not self-gravitate. The second is a realistic collapse scenario, like that of magnetised massive compact stars. Looking for answers to Melvin's question, we find that magnetic fields do show a generic tendency to support against relativistic gravitational implosion. The magnetic effect in question involves the curvature of the host space, it grows stronger as the collapse progresses and the triggering agent is the field's tension. In other words, the reason is the elasticity of the field lines and their inherent tendency to react against anything that distorts them from equilibrium. In this case, the distorting agent is the ever increasing curvature of the collapsing host space.


title: Interplay of electron correlations, spin-orbit couplings, and structural effects for Cu centers in the quasi-two-dimensional magnet InCu$_{2/3}$V$_{1/3}$O$_3$

abstract: Less common ligand coordination of transition-metal centers is often associated with peculiar valence-shell electron configurations and outstanding physical properties. One example is the Fe$^+$ ion with linear coordination, actively investigated in the research area of single-molecule magnetism. Here we address the nature of 3$d^9$ states for Cu$^{2+}$ ions sitting in the center of trigonal bipyramidal ligand cages in the quasi-two-dimensional honeycomb compound InCu$_{2/3}$V$_{1/3}$O$_3$, whose unusual magnetic properties were intensively studied in the recent past. In particular, we discuss the interplay of structural effects, electron correlations, and spin-orbit couplings in this material. A relevant computational finding is a different sequence of the Cu ($xz$, $yz$) and ($xy$, $x^2\!-\!y^2$) levels as compared to existing electronic-structure models, which has implications for the interpretation of various excitation spectra. Spin-orbit interactions, both first- and second-order, turn out to be stronger than previously assumed, suggesting that rather rich single-ion magnetic properties can be in principle achieved also for the 3$d^9$ configuration by properly adjusting the sequence of crystal-field states for such less usual ligand coordination.


title: Stability of rotating scalar boson stars with nonlinear interactions

abstract: We study the stability of rotating scalar boson stars, comparing those made from a simple massive complex scalar (referred to as mini boson stars), to those with several different types of nonlinear interactions. To that end, we numerically evolve the nonlinear Einstein-Klein-Gordon equations in 3D, beginning with stationary boson star solutions. We show that the linear, non-axisymmetric instability found in mini boson stars with azimuthal number $m=1$ persists across the entire parameter space for these stars, though the timescale diverges in the Newtonian limit. Therefore, any boson star with $m=1$ that is sufficiently far into the non-relativistic regime, where the leading order mass term dominates, will be unstable, independent of the nonlinear scalar self-interactions. However, we do find regions of $m=1$ boson star parameter space where adding nonlinear interactions to the scalar potential quenches the non-axisymmetric instability, both on the non-relativistic, and the relativistic branches of solutions. We also consider select boson stars with $m=2$, finding instability in all cases. For the cases exhibiting instability, we follow the nonlinear development, finding a range of dynamics including fragmentation into multiple unbound non-rotating stars, and formation of binary black holes. Finally, we comment on the relationship between stability and criteria based on the rotating boson star's frequency in relation to that of a spherical boson star or the existence of a co-rotation point. The boson stars that we find not to exhibit instability when evolved for many dynamical times include rapidly rotating cases where the compactness is comparable to that of a black hole or neutron star.


title: Sufficient Conditions for Feasibility of Optimal Control Problems Using Control Barrier Functions

abstract: It has been shown that satisfying state and control constraints while optimizing quadratic costs subject to desired (sets of) state convergence for affine control systems can be reduced to a sequence of quadratic programs (QPs) by using Control Barrier Functions (CBFs) and Control Lyapunov Functions (CLFs). One of the main challenges in this approach is ensuring the feasibility of these QPs, especially under tight control bounds and safety constraints of high relative degree. In this paper, we provide sufficient conditions for guranteed feasibility. The sufficient conditions are captured by a single constraint that is enforced by a CBF, which is added to the QPs such that their feasibility is always guaranteed. The additional constraint is designed to be always compatible with the existing constraints, therefore, it cannot make a feasible set of constraints infeasible - it can only increase the overall feasibility. We illustrate the effectiveness of the proposed approach on an adaptive cruise control problem.


title: Diffuse Radio Emission from Galaxy Clusters in the LOFAR Two-metre Sky Survey Deep Fields

abstract: Low-frequency radio observations are revealing an increasing number of diffuse synchrotron sources from galaxy clusters, dominantly in the form of radio halos or radio relics. The existence of this diffuse synchrotron emission indicates the presence of relativistic particles and magnetic fields. It is still an open question what mechanisms exactly are responsible for the population of relativistic electrons driving this synchrotron emission. The LOFAR Two-metre Sky Survey Deep Fields offer a unique view of this problem. Reaching noise levels below 30 $\mu$Jy/beam, these are the deepest images made at the low frequency of 144 MHz. This paper presents a search for diffuse emission in galaxy clusters in the first data release of the LOFAR Deep Fields. We detect a new high-redshift radio halo with a flux density of $8.9 \pm 1.0$ mJy and corresponding luminosity of $P_{144\mathrm{MHz}}=(3.6 \pm 0.6)\times10^{25}$ W Hz$^{-1}$ in an X-ray detected cluster at $z=0.77$ with a mass estimate of $M_{500} = 3.3_{-1.7}^{+1.1} \times 10^{14} M_\odot.$ Deep upper limits are placed on clusters with non-detections. We compare the results to the correlation between halo luminosity and cluster mass derived for radio halos found in the literature. This study is one of few to find diffuse emission in low mass ($M_{500} < 5\times10^{14} M_\odot$) systems and shows that deep low-frequency observations of galaxy clusters are fundamental for opening up a new part of parameter space in the study of non-thermal phenomena in galaxy clusters.


title: Improved Load Balancing in Large Scale Systems using Attained Service Time Reporting

abstract: Our interest lies in load balancing jobs in large scale systems consisting of multiple dispatchers and FCFS servers. In the absence of any information on job sizes, dispatchers typically use queue length information reported by the servers to assign incoming jobs. When job sizes are highly variable, using only queue length information is clearly suboptimal and performance can be improved if some indication can be provided to the dispatcher about the size of an ongoing job. In a FCFS server measuring the attained service time of the ongoing job is easy and servers can therefore report this attained service time together with the queue length when queried by a dispatcher. In this paper we propose and analyse a variety of load balancing policies that exploit both the queue length and attained service time to assign jobs, as well as policies for which only the attained service time of the job in service is used. We present a unified analysis for all these policies in a large scale system under the usual asymptotic independence assumptions. The accuracy of the proposed analysis is illustrated using simulation. We present extensive numerical experiments which clearly indicate that a significant improvement in waiting (and thus also in response) time may be achieved by using the attained service time information on top of the queue length of a server. Moreover, the policies which do not make use of the queue length still provide an improved waiting time for moderately loaded systems.


title: Cosmic microwave background lensing with optimal convergence and shear estimators

abstract: We present the optimal convergence and shear estimators for lensing reconstruction from the cosmic microwave background temperature field. This generalizes the deflection estimator, is sensitive to non-lensing modes, provides internal consistency checks, and is always at least as optimal. Previously, these estimators were only known in the squeezed limit. This paper decomposes convergence and shear fields into cosine and sine waves and the lensed correlation function is then Taylor expanded in the wave amplitudes. Maximizing the likelihood function gives the optimal estimators for the convergence and shear fields. This method has the potential to improve the lensing reconstruction of the cosmic microwave background polarization field: the shear and convergence can be optimally combined to form a deflection estimator, or used separately to separate non-lensing modes, or utilize lensing of non-Gaussian secondary foregrounds.


title: Modelling heterogeneous dust particles. An application to cometary polarization

abstract: In this work, we introduce a comet dust model that incorporates multiple dust morphologies along with inhomogeneous mixture of silicate minerals and carbonaceous materials under power-law size distribution, to replicate the standard polarization-phase curve observed in several comets in the narrow-band continuum. Following the results from Rosetta/MIDAS and COSIMA, we create high porosity Hierarchical Aggregates (HA) and low porosity (< 10$\%$) Solids in the form of agglomerated debris. We also introduce a moderate porosity structure with solids in the core, surrounded by fluffy aggregates called Fluffy Solids (FS). We study the mixing combinations, (HA and Solids), (HA and FS) and (HA, FS and Solids) for a range of power-law index n=2.0 to 3.0 for different sets of mixing percentage of silicate minerals and carbonaceous materials. Polarimetry of the short period comets 1P/Halley and 67P/Churyumov-Gerasimenko match best with the polarisation resulting from the combination of HA and Solids while the combinations (HA and FS) and (HA, FS and Solids) provide the best fit results for the long period comets C/1995 O1 (Hale-Bopp) and C/1996 B2 (Hyakutake). The best fit model results also recreate the observed wavelength dependence of polarization. Our dust model agree with the idea that the long period comets may have high percentage of loose particles (HA and FS) compared to those in the case of short period comets as the short period comets experience more frequent and/or higher magnitude of weathering.


title: Heterogeneous Paxos: Technical Report

abstract: In distributed systems, a group of $\textit{learners}$ achieve $\textit{consensus}$ when, by observing the output of some $\textit{acceptors}$, they all arrive at the same value. Consensus is crucial for ordering transactions in failure-tolerant systems. Traditional consensus algorithms are homogeneous in three ways: - all learners are treated equally, - all acceptors are treated equally, and - all failures are treated equally. These assumptions, however, are unsuitable for cross-domain applications, including blockchains, where not all acceptors are equally trustworthy, and not all learners have the same assumptions and priorities. We present the first consensus algorithm to be heterogeneous in all three respects. Learners set their own mixed failure tolerances over differently trusted sets of acceptors. We express these assumptions in a novel $\textit{Learner Graph}$, and demonstrate sufficient conditions for consensus. We present $\textit{Heterogeneous Paxos}$: an extension of Byzantine Paxos. Heterogeneous Paxos achieves consensus for any viable Learner Graph in best-case three message sends, which is optimal. We present a proof-of-concept implementation, and demonstrate how tailoring for heterogeneous scenarios can save resources and latency.


title: Personalized Cardiovascular Disease Risk Mitigation via Longitudinal Inverse Classification

abstract: Cardiovascular disease (CVD) is a serious illness affecting millions world-wide and is the leading cause of death in the US. Recent years, however, have seen tremendous growth in the area of personalized medicine, a field of medicine that places the patient at the center of the medical decision-making and treatment process. Many CVD-focused personalized medicine innovations focus on genetic biomarkers, which provide person-specific CVD insights at the genetic level, but do not focus on the practical steps a patient could take to mitigate their risk of CVD development. In this work we propose longitudinal inverse classification, a recommendation framework that provides personalized lifestyle recommendations that minimize the predicted probability of CVD risk. Our framework takes into account historical CVD risk, as well as other patient characteristics, to provide recommendations. Our experiments show that earlier adoption of the recommendations elicited from our framework produce significant CVD risk reduction.


title: Learning differential equation models from stochastic agent-based model simulations

abstract: Agent-based models provide a flexible framework that is frequently used for modelling many biological systems, including cell migration, molecular dynamics, ecology, and epidemiology. Analysis of the model dynamics can be challenging due to their inherent stochasticity and heavy computational requirements. Common approaches to the analysis of agent-based models include extensive Monte Carlo simulation of the model or the derivation of coarse-grained differential equation models to predict the expected or averaged output from the agent-based model. Both of these approaches have limitations, however, as extensive computation of complex agent-based models may be infeasible, and coarse-grained differential equation models can fail to accurately describe model dynamics in certain parameter regimes. We propose that methods from the equation learning field provide a promising, novel, and unifying approach for agent-based model analysis. Equation learning is a recent field of research from data science that aims to infer differential equation models directly from data. We use this tutorial to review how methods from equation learning can be used to learn differential equation models from agent-based model simulations. We demonstrate that this framework is easy to use, requires few model simulations, and accurately predicts model dynamics in parameter regions where coarse-grained differential equation models fail to do so. We highlight these advantages through several case studies involving two agent-based models that are broadly applicable to biological phenomena: a birth-death-migration model commonly used to explore cell biology experiments and a susceptible-infected-recovered model of infectious disease spread.


title: Excitation of Zonal Flow by Intermediate-Scale Toroidal Electron Temperature Gradient Turbulence

abstract: We show that zonal flow can be preferentially excited by intermediate-scale toroidal electron temperature gradient (ETG) turbulence in tokamak plasmas. Previous theoretical studies that yielded an opposite conclusion assumed a fluid approximation for ETG modes. Here, we carry out a gyrokinetic analysis which ultimately yields a nonlinear Schr\"{o}dinger equation for the ETG dynamics with a Navier-Stokes type nonlinearity. For typical tokamak parameters, it is found that zonal flow generation plays an important role in the intermediate-scale ETG turbulence. This finding offers an explanation for recent multi-scale gyrokinetic simulations.


title: Oort cloud Ecology II: Extra-solar Oort clouds and the origin of asteroidal interlopers

abstract: We simulate the formation and evolution of Oort clouds around the 200 nearest stars (within 16pc according to the Gaia DR2) database. This study is performed by numerically integrating the planets and minor bodies in orbit around the parent star and in the Galactic potential. The calculations start 1\,Gyr ago and continue for 100Myr into the future. In this time frame, we simulate how asteroids (and planets) are ejected from the star's vicinity and settle in an Oort cloud and how they escape the local stellar gravity to form tidal steams. A fraction of 0.0098 to 0.026 of the asteroids remain bound to their parent star. The orbits of these asteroids isotropizes and circularizes due to the influence of the Galactic tidal field to eventually form an Oort cloud between 10^4 and 2 10^5au. We estimate that 6% of the nearby stars may have a planet in its Oort cloud. The majority of asteroids (and some of the planets) become unbound from the parent star to become free floating in the Galactic potential. These soli lapides remain in a similar orbit around the Galactic center as their host star, forming dense streams of rogue interstellar asteroids and planets. The Solar system occasionally passes through such tidal streams, potentially giving rise to occasional close encounters with object in this stream. The two recently discovered objects, 1I/(2017 Q3) 'Oumuamua and 2I/(2019 Q4) Borisov, may be such objects. Although the direction from which an individual solus lapis originated cannot easily be traced back to the original host, multiple such objects coming from the same source might help to identify their origin. At the moment the Solar system is in the bow or wake of the tidal stream of 10 of the nearby stars which might contribute considerably to the interaction rate. (abridged)


title: Extending Charged Holographic R\'enyi Entropy

abstract: Motivated by extended black hole thermodynamics, we generalize the R\'enyi entropy of charged holographic conformal field theories (CFTs) in $d$-dimensions. Specifically, following (1807.09215), we extend the quench description of the R\'enyi entropy of globally charged holographic CFTs by including pressure variations of charged hyperbolically sliced anti de Sitter black holes. We provide an exhaustive analysis of the new type of charged R\'enyi entropy, where we find an interesting interplay between a parameter controlling the pressure of the black hole and its charge. A field theoretic interpretation of this extended charged R\'enyi entropy is given. In particular, in $d=2$, where the bulk geometry becomes the charged Ba\~nados, Teitelboim, Zanelli black hole, we write down the extended charged R\'enyi entropy in terms of the twist operators of the charged field theory. An area law prescription for the extended R\'enyi entropy is formulated. We comment on several avenues for future work, including how global charge conservation relates to black hole super-entropicity.


title: On the Bezrukavnikov-Kaledin quantization of symplectic varieties in characteristic $p$

abstract: We prove that after inverting the Planck constant $h$ the Bezrukavnikov-Kaledin quantization $(X, \mathcal{O}_h)$ of symplectic variety $X$ in characteristic $p$ is Morita equivalent to a certain central reduction of the algebra of differential operators on $X$.


title: Inferring the properties of the sources of reionization using the morphological spectra of the ionized regions

abstract: High-redshift 21-cm observations will provide crucial insights into the physical processes of the Epoch of Reionization. Next-generation interferometers such as the Square Kilometer Array will have enough sensitivity to directly image the 21-cm fluctuations and trace the evolution of the ionizing fronts. In this work, we develop an inferential approach to recover the sources and IGM properties of the process of reionization using the number and, in particular, the morphological pattern spectra of the ionized regions extracted from realistic mock observations. To do so, we extend the Markov Chain Monte Carlo analysis tool 21CMMC by including these 21-cm tomographic statistics and compare this method to only using the power spectrum. We demonstrate that the evolution of the number-count and morphology of the ionized regions as a function of redshift provides independent information to disentangle multiple reionization scenarios because it probes the average ionizing budget per baryon. Although less precise, we find that constraints inferred using 21-cm tomographic statistics are more robust to the presence of contaminants such as foreground residuals. This work highlights that combining power spectrum and tomographic analyses more accurately recovers the astrophysics of reionization.


title: Large-scale kernelized GRANGER causality to infer topology of directed graphs with applications to brain networks

abstract: Graph topology inference of network processes with co-evolving and interacting time-series is crucial for network studies. Vector autoregressive models (VAR) are popular approaches for topology inference of directed graphs; however, in large networks with short time-series, topology estimation becomes ill-posed. The present paper proposes a novel nonlinearity-preserving topology inference method for directed networks with co-evolving nodal processes that solves the ill-posedness problem. The proposed method, large-scale kernelized Granger causality (lsKGC), uses kernel functions to transform data into a low-dimensional feature space and solves the autoregressive problem in the feature space, then finds the pre-images in the input space to infer the topology. Extensive simulations on synthetic datasets with nonlinear and linear dependencies and known ground-truth demonstrate significant improvement in the Area Under the receiver operating characteristic Curve ( AUC ) of the receiver operating characteristic for network recovery compared to existing methods. Furthermore, tests on real datasets from a functional magnetic resonance imaging (fMRI) study demonstrate 96.3 percent accuracy in diagnosis tasks of schizophrenia patients, which is the highest in the literature with only brain time-series information.


title: A Probabilistic Approach in Historical Linguistics Word Order Change in Infinitival Clauses: from Latin to Old French

abstract: This research offers a new interdisciplinary approach to the field of Linguistics by using Computational Linguistics, NLP, Bayesian Statistics and Sociolinguistics methods. This thesis investigates word order change in infinitival clauses from Object-Verb (OV) to Verb-Object (VO) in the history of Latin and Old French. By applying a variationist approach, I examine a synchronic word order variation in each stage of language change, from which I infer the character, periodization and constraints of diachronic variation. I also show that in discourse-configurational languages, such as Latin and Early Old French, it is possible to identify pragmatically neutral contexts by using information structure annotation. I further argue that by mapping pragmatic categories into a syntactic structure, we can detect how word order change unfolds. For this investigation, the data are extracted from annotated corpora spanning several centuries of Latin and Old French and from additional resources created by using computational linguistic methods. The data are then further codified for various pragmatic, semantic, syntactic and sociolinguistic factors. This study also evaluates previous factors proposed to account for word order alternation and change. I show how information structure and syntactic constraints change over time and propose a method that allows researchers to differentiate a stable word order alternation from alternation indicating a change. Finally, I present a three-stage probabilistic model of word order change, which also conforms to traditional language change patterns.


title: Minimal Zeeman field requirement for a topological transition in superconductors

abstract: Platforms for creating Majorana quasiparticles rely on superconductivity and breaking of time-reversal symmetry. By studying continuous deformations to known trivial states, we find that the relationship between superconducting pairing and time reversal breaking imposes rigorous bounds on the topology of the system. Applying these bounds to $s$-wave systems with a Zeeman field, we conclude that a topological phase transition requires that the Zeeman energy at least locally exceed the superconducting pairing by the energy gap of the full Hamiltonian. Our results are independent of the geometry and dimensionality of the system.


title: Spectra of locally matrix algebras

abstract: We describe spectra of associative (not necessarily unital and not necessarily countable-dimensional) locally matrix algebras. We determine all possible spectra of locally matrix algebras and give a new proof of Dixmier-Baranov Theorem. As an application of our description of spectra we determine embeddings of locally matrix algebras.


title: Observation-infused simulations of high-speed boundary layer transition

abstract: High-speed boundary-layer transition is extremely sensitive to the free-stream disturbances which are often uncertain. This uncertainty compromises predictions of models and simulations. To enhance the fidelity of simulations, we directly infuse them with available observations. Our methodology is general and can be adopted with any simulation tool, and is herein demonstrated using direct numerical simulations. An ensemble variational (EnVar) optimization is performed, whereby we determine the upstream flow that optimally reproduces the observations. The cost functional accounts for our relative confidence in the model and the observations, and judicious choice of the ensemble members improves convergence and reduces the prediction uncertainty. We demonstrate our observation-infused predictions for boundary-layer transition at Mach 4.5. Without prior knowledge of the free-stream condition, and using only observations of wall pressure at isolated locations from an independent computation (true flow), all of the relevant inflow disturbances are identified. We then evaluate the entire flow field, beyond the original limited wall observations, and interpret simulations consistently with data and vice versa. Our predicted flow compares favorably to the true `unknown' state, and discrepancies are analyzed in detail. We also examine the impact of weighting of observations. Improved convergence of the inverse problem and accuracy of the inflow amplitudes and phases are demonstrated, and are explained by aid of a simple example from two-dimensional unstable, chaotic convection.


title: Individual Dynamical Masses of DENIS J063001.4$-$184014AB Reveal A Likely Young Brown Dwarf Triple

abstract: The binary nature of the M8.5 dwarf DENIS J063001.4$-$184014AB (DE0630$-$18) was discovered with astrometric monitoring from the ground, which determined the unresolved photocentric orbit and the trigonometric parallax of the system. Here we present radial-velocity monitoring and resolved observations in the near-infrared with Keck aperture masking that allow us to measure the system's relative separation and brightness. By combining all available information, we determine the individual dynamical masses of the binary components to be $M_1 = 0.052^{+0.009}_{-0.008} M_\mathrm{Sun}$ and $M_2 = 0.052^{+0.005}_{-0.004} M_\mathrm{Sun}$, both firmly in the substellar regime. These masses are surprising given the object's M8.5 optical spectral type and equivalent absolute magnitude, and the significant difference in brightness between the components ($\Delta{K}$ = 1.74$\pm$0.06 mag). Our results suggest that DE0630$-$18 is a relatively young system ($\sim$200 Myr) with a secondary component that is itself a potentially unresolved binary.


title: Hausdorff dimension of escaping sets of meromorphic functions II

abstract: A function which is transcendental and meromorphic in the plane has at least two singular values. On one hand, if a meromorphic function has exactly two singular values, it is known that the Hausdorff dimension of the escaping set can only be either $2$ or $1/2$. On the other hand, the Hausdorff dimension of escaping sets of Speiser functions can attain every number in $[0,2]$ (cf. \cite{ac1}). In this paper, we show that number of singular values which is needed to attain every Hausdorff dimension of escaping sets is not more than $4$.


title: ACCORDANT: A Domain Specific Model and DevOpsApproach for Big Data Analytics Architectures

abstract: Big data analytics (BDA) applications use machine learning algorithms to extract valuable insights from large, fast, and heterogeneous data sources. New software engineering challenges for BDA applications include ensuring performance levels of data-driven algorithms even in the presence of large data volume, velocity, and variety (3Vs). BDA software complexity frequently leads to delayed deployments, longer development cycles and challenging performance assessment. This paper proposes a Domain-Specific Model (DSM), and DevOps practices to design, deploy, and monitor performance metrics in BDA applications. Our proposal includes a design process, and a framework to define architectural inputs, software components, and deployment strategies through integrated high-level abstractions to enable QS monitoring. We evaluate our approach with four use cases from different domains to demonstrate a high level of generalization. Our results show a shorter deployment and monitoring times, and a higher gain factor per iteration compared to similar approaches.


title: Inter-regional correlation estimators for functional magnetic resonance imaging

abstract: Functional magnetic resonance imaging (fMRI) functional connectivity between brain regions is often computed using parcellations defined by functional or structural atlases. Typically, some kind of voxel averaging is performed to obtain a single temporal correlation estimate per region pair. However, several estimators can be defined for this task, with various assumptions and degrees of robustness to local noise, global noise, and region size. In this paper, we systematically present and study the properties of 9 different functional connectivity estimators taking into account the spatial structure of fMRI data, based on a simple fMRI data spatial model. These include 3 existing estimators and 6 novel estimators. We demonstrate the empirical properties of the estimators using synthetic, animal, and human data, in terms of graph structure, repeatability and reproducibility, discriminability, dependence on region size, as well as local and global noise robustness. We prove analytically the link between regional intra-correlation and inter-region correlation, and show that the choice of estimator has a strong influence on inter-correlation values.


title: Gravitational decoherence of photons

abstract: Models of gravitational decoherence are not commonly applied to ultra-relativistic systems, including photons. As a result, few quantum optical tests of gravitational decoherence have been developed. In this paper, we generalize the gravitational decoherence model of Anastopoulos and Hu (arXiv:1305.5231) to photons. In this model, decoherence originates from a bath of stochastic gravitational perturbations, possibly of fundamental origin. We derive a master equation for general states of the electromagnetic (EM) field; the only free parameter is a noise temperature $\Theta$ of the gravitational fluctuations. We find that interference experiments with long baselines, accessible in near-future experiments, can, in principle, lead to strong constraints in $\Theta$.


title: Probabilistic Mapping of Dark Matter by Neural Score Matching

abstract: The Dark Matter present in the Large-Scale Structure of the Universe is invisible, but its presence can be inferred through the small gravitational lensing effect it has on the images of far away galaxies. By measuring this lensing effect on a large number of galaxies it is possible to reconstruct maps of the Dark Matter distribution on the sky. This, however, represents an extremely challenging inverse problem due to missing data and noise dominated measurements. In this work, we present a novel methodology for addressing such inverse problems by combining elements of Bayesian statistics, analytic physical theory, and a recent class of Deep Generative Models based on Neural Score Matching. This approach allows to do the following: (1) make full use of analytic cosmological theory to constrain the 2pt statistics of the solution, (2) learn from cosmological simulations any differences between this analytic prior and full simulations, and (3) obtain samples from the full Bayesian posterior of the problem for robust Uncertainty Quantification. We present an application of this methodology on the first deep-learning-assisted Dark Matter map reconstruction of the Hubble Space Telescope COSMOS field.


title: NLPGym -- A toolkit for evaluating RL agents on Natural Language Processing Tasks

abstract: Reinforcement learning (RL) has recently shown impressive performance in complex game AI and robotics tasks. To a large extent, this is thanks to the availability of simulated environments such as OpenAI Gym, Atari Learning Environment, or Malmo which allow agents to learn complex tasks through interaction with virtual environments. While RL is also increasingly applied to natural language processing (NLP), there are no simulated textual environments available for researchers to apply and consistently benchmark RL on NLP tasks. With the work reported here, we therefore release NLPGym, an open-source Python toolkit that provides interactive textual environments for standard NLP tasks such as sequence tagging, multi-label classification, and question answering. We also present experimental results for 6 tasks using different RL algorithms which serve as baselines for further research. The toolkit is published at https://github.com/rajcscw/nlp-gym


title: Machine Learning and Soil Humidity Sensing: Signal Strength Approach

abstract: The IoT vision of ubiquitous and pervasive computing gives rise to future smart irrigation systems comprising physical and digital world. Smart irrigation ecosystem combined with Machine Learning can provide solutions that successfully solve the soil humidity sensing task in order to ensure optimal water usage. Existing solutions are based on data received from the power hungry/expensive sensors that are transmitting the sensed data over the wireless channel. Over time, the systems become difficult to maintain, especially in remote areas due to the battery replacement issues with large number of devices. Therefore, a novel solution must provide an alternative, cost and energy effective device that has unique advantage over the existing solutions. This work explores a concept of a novel, low-power, LoRa-based, cost-effective system which achieves humidity sensing using Deep learning techniques that can be employed to sense soil humidity with the high accuracy simply by measuring signal strength of the given underground beacon device.


title: A simple way to compute structure constants of semi-simple Lie algebras

abstract: The standard way to compute the structure constants of semi-simple Lie algebras involves the additive structure of the roots. In earlier work, I described how ideas of Jacques Tits could be applied to do this by using the structure of the Weyl group. In this work, I explain how to combine some of the earlier ideas with some suggestions of Kottwitz to present a very simple algorithm which comes close to constructing a canonical basis of the Lie algebra. Unfortunately, unlike previous methods, it does not work for arbitrary Kac-Moody algebra.


title: Fat tails arise endogenously in asset prices from supply/demand, with or without jump processes

abstract: We show that the quotient of Levy processes of jump-diffusion type has a fat-tailed distribution. An application is to price theory in economics. We show that fat tails arise endogenously from modeling of price change based on an excess demand analysis resulting in a quotient of arbitrarily correlated demand and supply whether or not jump discontinuities are present. The assumption is that supply and demand are described by drift terms, Brownian (i.e., Gaussian) and compound Poisson jump processes. If $P^{-1}dP/dt$ (the relative price change in an interval $dt$) is given by a suitable function of relative excess demand, $\left( \mathcal{D}% -\mathcal{S}\right) /\mathcal{S}$ (where $\mathcal{D}$ and $\mathcal{S}$ are demand and supply), then the distribution has tail behavior $F\left( x\right) \sim x^{-\zeta}$ for a power $\zeta$ that depends on the function $G$ in $P^{-1}dP/dt=G\left( \mathcal{D}/\mathcal{S}\right) $. For $G\left( x\right) \sim\left\vert x\right\vert ^{1/q}$ one has $\zeta=q.$ The empirical data for assets typically yields a value, $\zeta\tilde{=}3,$ or $\ \zeta \in\left[ 3,5\right] $ for some markets. The discrepancy between the empirical result and theory never arises if one models price dynamics using basic economics methodology, i.e., generalized Walrasian adjustment, rather than the usual starting point for classical finance which assumes a normal distribution of price changes. The function $G$ is deterministic, and can be calibrated with a smaller data set. The results establish a simple link between the decay exponent of the density function and the price adjustment function, a feature that can improve methodology for risk assessment. The mathematical results can be applied to other problems involving the relative difference or quotient of Levy processes of jump-diffusion type.


title: Void Galaxies Follow a Distinct Evolutionary Path in the Environmental COntext Catalog

abstract: We measure the environmental dependence, where environment is defined by the distance to the third nearest neighbor, of multiple galaxy properties inside the Environmental COntext (ECO) catalog. We focus primarily on void galaxies, which we define as the $10 \%$ of galaxies having the lowest local density. We compare the properties of void and non-void galaxies: baryonic mass, color, fractional stellar mass growth rate (FSMGR), morphology, and gas-to-stellar-mass ratio (estimated from a combination of HI data and photometric gas fractions calibrated with the RESOLVE survey). Our void galaxies typically have lower baryonic masses than galaxies in denser environments, and they display the properties expected of a lower mass population: they have more late-types, are bluer, have higher FSMGR, and are more gas rich. We control for baryonic mass and investigate the extent to which void galaxies are different at fixed mass. Void galaxies are bluer, more gas-rich, and more star forming at fixed mass than non-void galaxies, which is a possible signature of galaxy assembly bias. Furthermore, we show that these trends persist even at fixed mass and morphology, and we find that voids host a distinct population of early-types that are bluer and more star-forming than the typical red and quenched early-types. In addition to these empirical observational results, we also present theoretical results from mock catalogs with built-in galaxy assembly bias. We show that a simple matching of galaxy properties to (sub)halo properties, such as mass and age, can recover the observed environmental trends in ECO galaxies.


title: Where Are You? Localization from Embodied Dialog

abstract: We present Where Are You? (WAY), a dataset of ~6k dialogs in which two humans -- an Observer and a Locator -- complete a cooperative localization task. The Observer is spawned at random in a 3D environment and can navigate from first-person views while answering questions from the Locator. The Locator must localize the Observer in a detailed top-down map by asking questions and giving instructions. Based on this dataset, we define three challenging tasks: Localization from Embodied Dialog or LED (localizing the Observer from dialog history), Embodied Visual Dialog (modeling the Observer), and Cooperative Localization (modeling both agents). In this paper, we focus on the LED task -- providing a strong baseline model with detailed ablations characterizing both dataset biases and the importance of various modeling choices. Our best model achieves 32.7% success at identifying the Observer's location within 3m in unseen buildings, vs. 70.4% for human Locators.


title: A bibliometric methodology to unveil territorial inequities in the scientific wealth to combat COVID-19

abstract: In this paper we develop a methodology to assess the scientific wealth of territories at field level. Our methodology uses a bibliometric approach based on the observation of academic research performance and overall scientific production in each territory. We apply it to assess disparities in the Italian territories in the medical specialties at the front line of the COVID-19 emergency. Italy has been the first among western countries to be severely affected by the onset of the COVID-19 pandemic. The study reveals remarkable inequities across territories, with scientific weaknesses concentrated in the south. Policies for rebalancing the north-south divide should also consider, in addition to tangible assets, the gap in production and availability of quality medical knowledge.


title: Imaging the Dusty Substructures due to Terrestrial Planets in Planet-forming Disks with ALMA and the Next Generation Very Large Array

abstract: We present simulations of the capabilities of the Atacama Large Millimeter/submillimeter Array (ALMA) and of a Next Generation Very Large Array (ngVLA) to detect and resolve substructures due to terrestrial planets and Super-Earths in nearby planet-forming disks. We adopt the results of global 2-D hydrodynamical planet-disk simulations that account for the dynamics of gas and dust in a disk with an embedded planet. Our simulations follow the combined evolution of gas and dust for several thousand planetary orbits. We show that long integrations (several tens of hours) with the ngVLA can detect and spatially resolve dust structures due to low-mass rocky planets in the terrestrial planet formation regions of nearby disks (stellocentric radii $r = 1 - 3$ au), under the assumption that the disk viscosity in those regions is low ($\alpha \le 10^{-5}$). ALMA is instead unable to resolve these structures in these disk regions. We also show that high-resolution ngVLA observations separated by several days to few weeks would allow to detect the proper motion of the azimuthally asymmetric structures expected in the disk regions of terrestrial planet formation.


title: Measurements of Higgs Bosons Decaying to Bottom Quarks from Vector Boson Fusion Production with the ATLAS Experiment at $\sqrt{s}=$ 13 TeV

abstract: The paper presents a measurement of the Standard Model Higgs Boson decaying to $b$-quark pairs in the vector boson fusion (VBF) production mode. A sample corresponding to 126 fb$^{-1}$ of $\sqrt{s} =$ 13 TeV proton-proton collision data, collected with the ATLAS experiment at the Large Hadron Collider, is analyzed utilizing an adversarial neural network for event classification. The signal strength, defined as the ratio of the measured signal yield to that predicted by the Standard Model for VBF Higgs production, is measured to be $0.95^{+0.38}_{-0.36}$, corresponding to an observed (expected) significance of 2.6 (2.8) standard deviations from the background only hypothesis. The results are additionally combined with an analysis of Higgs bosons decaying to $b$-quarks, produced via VBF in association with a photon.


title: Avoiding Communication in Logistic Regression

abstract: Stochastic gradient descent (SGD) is one of the most widely used optimization methods for solving various machine learning problems. SGD solves an optimization problem by iteratively sampling a few data points from the input data, computing gradients for the selected data points, and updating the solution. However, in a parallel setting, SGD requires interprocess communication at every iteration. We introduce a new communication-avoiding technique for solving the logistic regression problem using SGD. This technique re-organizes the SGD computations into a form that communicates every $s$ iterations instead of every iteration, where $s$ is a tuning parameter. We prove theoretical flops, bandwidth, and latency upper bounds for SGD and its new communication-avoiding variant. Furthermore, we show experimental results that illustrate that the new Communication-Avoiding SGD (CA-SGD) method can achieve speedups of up to $4.97\times$ on a high-performance Infiniband cluster without altering the convergence behavior or accuracy.


title: Covariance matrix testing in high dimension using random projections

abstract: Estimation and hypothesis tests for the covariance matrix in high dimensions is a challenging problem as the traditional multivariate asymptotic theory is no longer valid. When the dimension is larger than or increasing with the sample size, standard likelihood based tests for the covariance matrix have poor performance. Existing high dimensional tests are either computationally expensive or have very weak control of type I error. In this paper, we propose a test procedure, CRAMP, for testing hypotheses involving one or more covariance matrices using random projections. Projecting the high dimensional data randomly into lower dimensional subspaces alleviates of the curse of dimensionality, allowing for the use of traditional multivariate tests. An extensive simulation study is performed to compare CRAMP against asymptotics-based high dimensional test procedures. An application of the proposed method to two gene expression data sets is presented.


title: Center of Poisson and skein algebras associated to loops on surfaces

abstract: We discuss and develop a systematic method to compute the Poisson center (Casimir) of various Poisson algebras associated to loops on orientable surfaces (possibly with boundary and punctures) introduced by Goldman and Wolpert in 80's while studying Thurston's earthquakes deformations. Our computation extends a result of Etingof to all finite type hyperbolic surfaces. We use these methods to compute the center of various skein algebras introduced by Turaev for the quantization of these Poisson algebras. As another application of our results we compute the center of homotopy skein algebra introduced by Hoste and Przytycki.


title: Tsirelson's Bound and the Quantum Monogamy Bound from Global Determinism

abstract: We demonstrate that in a globally deterministic universe, all spatiotemporally symmetric processes must obey counterfactual parameter independence. We show that the Tsirelson bound can be derived from counterfactual parameter independence. We show that the quantum monogamy bound can also be obtained from global determinism, and then then use these results to propose a novel solution to the horizon problem. We also explain how global determinism relates to contextuality in quantum mechanics.


title: Dual simulation of a Polyakov loop model at finite baryon density: phase diagram and local observables

abstract: Many Polyakov loop models can be written in a dual formulation which is free of sign problem even when a non-vanishing baryon chemical potential is introduced in the action. Here, results of numerical simulations of a dual representation of one such effective Polyakov loop model at finite baryon density are presented. We compute various local observables such as energy density, baryon density, quark condensate and describe in details the phase diagram of the model. The regions of the first order phase transition and the crossover, as well as the line of the second order phase transition, are established. We also compute several correlation functions of the Polyakov loops.


title: Single-particle entanglement gives rise to truly nonlocal effects like single-particle steering

abstract: In 1927, at the Solvay conference, Einstein posed a thought experiment with the primary intention of showing the incompleteness of quantum mechanics; to prove it, he uses the instantaneous nonlocal effects caused by the collapse of the wave function of a single particle -the spooky action at a distance-, when a measurement is done. This historical event precede the well-know Einstein-Podolsk-Rosen criticism over the incompleteness of quantum mechanics. Here, by using the Stern-Gerlach experiment (SGE), we demonstrate how the instantaneous nonlocal feature of the collapse of the wave function together with the single-particle entanglement can be used to produce the nonlocal effect of steering. In the steering process Bob gets a quantum state depending on which observable Alice decides to measure. To accomplish this, we fully exploit the spreading (over large distances) of the entangled wave function of the single-particle. In particular, we demonstrate that the nonlocality of the single-particle entanglement allows the particle to know which detector Alice is using to steer Bob's state. Therefore, notwithstanding strong counterarguments, we prove that the single-particle entanglement gives rise to truly nonlocal effects at two far a away places. This open the possibility of using the single-particle entanglement for implementing truly nonlocal task.


title: On inner automorphisms preserving subspaces of Clifford algebras

abstract: In this paper, we consider inner automorphisms that leave invariant fixed subspaces of real and complex Clifford algebras -- subspaces of fixed grades and subspaces determined by the reversion and the grade involution. We present groups of elements that define such inner automorphisms and study their properties. Some of these Lie groups can be interpreted as generalizations of Clifford, Lipschitz, and spin groups. We study the corresponding Lie algebras. Some of the results can be reformulated for the case of more general algebras -- graded central simple algebras or graded central simple algebras with involution.


title: Spherical objects, transitivity and auto-equivalences of Kodaira cycles via gentle algebras

abstract: This paper studies the class of spherical objects over any Kodaira $n$-cycle of projective lines and provides a parametrization of their isomorphism classes in terms of closed curves on the $n$-punctured torus without self-intersections. Employing recent results on gentle algebras, we derive a topological model for the bounded derived category of any Kodaira cycle. The groups of triangle auto-equivalences of these categories are computed and are shown to act transitively on isomorphism classes of spherical objects. This answers a question by Polishchuk and extends earlier results by Burban-Kreussler and Lekili-Polishchuk. The description of auto-equivalences is further used to establish faithfulness of a mapping class group action defined by Sibilla. The final part describes the closed curves which correspond to vector bundles and simple vector bundles. This leads to an alternative proof of a result by Bodnarchuk-Drozd-Greuel which states that simple vector bundles on cycles of projective lines are uniquely determined by their multi-degree, rank and determinant. As a by-product we obtain a closed formula for the cyclic sequence of any simple vector bundle on $C_n$ as introduced by Burban-Drozd-Greuel.


title: Clifford Analysis with Indefinite Signature

abstract: We extend constructions of classical Clifford analysis to the case of indefinite non-degenerate quadratic forms. We define (p,q)-left- and right-monogenic functions by means of Dirac operators that factor a certain wave operator. We prove two different versions of Cauchy's integral formulas for these functions. The two formulas arise from dealing with singularities in distinct ways, and are inspired by the methods of [L, FL]. These results indicate the merit of these methods for dealing with singularities.


title: Estimating the correlation in network disturbance models

abstract: The Network Disturbance Model of Doreian (1989) expresses the dependency between observations taken at the vertices of a network by modelling the correlation between neighbouring vertices, using a single correlation parameter $\rho$. It has been observed that estimation of $\rho$ in dense graphs, using the method of Maximum Likelihood, leads to results that can be both biased and very unstable. In this paper, we sketch why this is the case, showing that the variability cannot be avoided, no matter how large the network. We also propose a more intuitive estimator of $\rho$, which shows little bias. The related Network Effects Model is briefly discussed.


title: A Two-Phase Approach for Abstractive Podcast Summarization

abstract: Podcast summarization is different from summarization of other data formats, such as news, patents, and scientific papers in that podcasts are often longer, conversational, colloquial, and full of sponsorship and advertising information, which imposes great challenges for existing models. In this paper, we focus on abstractive podcast summarization and propose a two-phase approach: sentence selection and seq2seq learning. Specifically, we first select important sentences from the noisy long podcast transcripts. The selection is based on sentence similarity to the reference to reduce the redundancy and the associated latent topics to preserve semantics. Then the selected sentences are fed into a pre-trained encoder-decoder framework for the summary generation. Our approach achieves promising results regarding both ROUGE-based measures and human evaluations.


title: LOFAR Deep Fields: Probing a broader population of polarized radio galaxies in ELAIS-N1

abstract: We present deep polarimetric observations of the European Large Area ISO Survey-North 1 (ELAIS-N1) field using the Low Frequency Array (LOFAR) at 114.9-177.4 MHz. The ELAIS-N1 field is part of the LOFAR Two-metre Sky Survey deep fields data release I. For six eight-hour observing epochs, we align the polarization angles and stack the 20"-resolution Stokes $Q$, $U$-parameter data cubes. This produces a 16 deg$^2$ image with 1$\sigma_{\rm QU}$ sensitivity of 26 $\mu$Jy/beam in the central area. In this paper, we demonstrate the feasibility of the stacking technique, and we generate a catalog of polarized sources in ELAIS-N1 and their associated Faraday rotation measures (RMs). While in a single-epoch observation we detect three polarized sources, this number increases by a factor of about three when we consider the stacked data, with a total of ten sources. This yields a surface density of polarized sources of one per 1.6 deg$^2$. The Stokes $I$ images of three of the ten detected polarized sources have morphologies resembling those of FR I radio galaxies. This represents a greater fraction of this type of source than previously found, which suggests that more sensitive observations may help with their detection.


title: Kinematic Analysis of a Protostellar Multiple System: Measuring the Protostar Masses and Assessing Gravitational Instability in the Disks of L1448 IRS3B and L1448 IRS3A

abstract: We present new Atacama Large Millimeter/submillimeter Array (ALMA) observations towards a compact (230~au separation) triple protostar system, L1448 IRS3B, at 879~\micron with \contbeam~resolution. Spiral arm structure within the circum-multiple disk is well resolved in dust continuum toward IRS3B, and we detect the known wide (2300~au) companion, IRS3A, also resolving possible spiral substructure. Using dense gas tracers, C17O, H13CO$+$, and H13CN, we resolve the Keplerian rotation for both the circum-triple disk in IRS3B and the disk around IRS3A. Furthermore, we use the molecular line kinematic data and radiative transfer modeling of the molecular line emission to confirm that the disks are in Keplerian rotation with fitted masses of $1.19^{+0.13}_{-0.07}$ for IRS3B-ab, $1.51^{+0.06}_{-0.07}$~Msun for IRS3A, and place an upper limit on the central protostar mass for the tertiary IRS3B-c of 0.2~Msun. We measure the mass of the fragmenting disk of IRS3B to be 0.29~Msun from the dust continuum emission of the circum-multiple disk and estimate the mass of the clump surrounding IRS3B-c to be 0.07~Msun. We also find that the disk around IRS3A has a mass of 0.04~Msun. By analyzing the Toomre~Q parameter, we find the IRS3A circumstellar disk is gravitationally stable (Q$>$5), while the IRS3B disk is consistent with a gravitationally unstable disk (Q$<$1) between the radii 200-500~au. This coincides with the location of the spiral arms and the tertiary companion IRS3B-c, supporting the hypothesis that IRS3B-c was formed in situ via fragmentation of a gravitationally unstable disk.


title: The contribution of discrete sources to the sky temperature at 144 MHz

abstract: In recent years, the level of the extragalactic radio background has become a point of considerable interest, with some lines of argument pointing to an entirely new cosmological synchrotron background. The contribution of the known discrete source population to the sky temperature is key to this discussion. Because of the steep spectral index of the excess over the cosmic microwave background, it is best studied at low frequencies where the signal is strongest. The Low-Frequency Array (LOFAR) wide and deep sky surveys give us the best constraints yet on the contribution of discrete extragalactic sources at 144 MHz, and in particular allow us to include contributions from diffuse, low-surface-brightness emission that could not be fully accounted for in previous work. We show that, even with these new data, known sources can still only account for around a quarter of the estimated extragalactic sky temperature at LOFAR frequencies.


title: Real-Time Radio Technology and Modulation Classification via an LSTM Auto-Encoder

abstract: Identification of the type of communication technology and/or modulation scheme based on detected radio signal are challenging problems encountered in a variety of applications including spectrum allocation and radio interference mitigation. They are rendered difficult due to a growing number of emitter types and varied effects of real-world channels upon the radio signal. Existing spectrum monitoring techniques are capable of acquiring massive amounts of radio and real-time spectrum data using compact sensors deployed in a variety of settings. However, state-of-the-art methods that use such data to classify emitter types and detect communication schemes struggle to achieve required levels of accuracy at a computational efficiency that would allow their implementation on low-cost computational platforms. In this paper, we present a learning framework based on an LSTM denoising auto-encoder designed to automatically extract stable and robust features from noisy radio signals, and infer modulation or technology type using the learned features. The algorithm utilizes a compact neural network architecture readily implemented on a low-cost computational platform while exceeding state-of-the-art accuracy. Results on realistic synthetic as well as over-the-air radio data demonstrate that the proposed framework reliably and efficiently classifies received radio signals, often demonstrating superior performance compared to state-of-the-art methods.


title: Discrete logarithm problem in some families of sandpile groups

abstract: Biggs proposed the sandpile group of certain modified wheel graphs for cryptosystems relying on the difficulty of the discrete logarithm problem. Blackburn and independently Shokrieh showed that the discrete logarithm problem is efficiently solvable. We study Shokrieh's method in cases of graphs such that the sandpile group is not cyclic, namely the square cycle graphs and the wheel graphs. Knowing generators of the group or the form of the pseudoinverse of the Laplacian matrix makes the problem more vulnerable. We also consider the discrete logarithm problem in case of the so-called subdivided banana graphs. In certain cases the sandpile group is cyclic and a generator is known and one can solve the discrete logarithm problem without computing the pseudoinverse of the Laplacian matrix.


title: Clinical Micro-CT Empowered by Interior Tomography, Robotic Scanning, and Deep Learning

abstract: While micro-CT systems are instrumental in preclinical research, clinical micro-CT imaging has long been desired with cochlear implantation as a primary example. The structural details of the cochlear implant and the temporal bone require a significantly higher image resolution than that (about 0.2 mm) provided by current medical CT scanners. In this paper, we propose a clinical micro-CT (CMCT) system design integrating conventional spiral cone-beam CT, contemporary interior tomography, deep learning techniques, and technologies of micro-focus X-ray source, photon-counting detector (PCD), and robotic arms for ultrahigh resolution localized tomography of a freely-selected volume of interest (VOI) at a minimized radiation dose level. The whole system consists of a standard CT scanner for a clinical CT exam and VOI specification, and a robotic-arm based micro-CT scanner for a local scan at much higher spatial and spectral resolution as well as much reduced radiation dose. The prior information from global scan is also fully utilized for background compensation to improve interior tomography from local data for accurate and stable VOI reconstruction. Our results and analysis show that the proposed hybrid reconstruction algorithm delivers superior local reconstruction, being insensitive to the misalignment of the isocenter position and initial view angle in the data/image registration while the attenuation error caused by scale mismatch can be effectively addressed with bias correction. These findings demonstrate the feasibility of our system design. We envision that deep learning techniques can be leveraged for optimized imaging performance. With high resolution imaging, high dose efficiency and low system cost synergistically, our proposed CMCT system has great potentials in temporal bone imaging as well as various other clinical applications.


title: Facebook AI's WMT20 News Translation Task Submission

abstract: This paper describes Facebook AI's submission to WMT20 shared news translation task. We focus on the low resource setting and participate in two language pairs, Tamil <-> English and Inuktitut <-> English, where there are limited out-of-domain bitext and monolingual data. We approach the low resource problem using two main strategies, leveraging all available data and adapting the system to the target news domain. We explore techniques that leverage bitext and monolingual data from all languages, such as self-supervised model pretraining, multilingual models, data augmentation, and reranking. To better adapt the translation system to the test domain, we explore dataset tagging and fine-tuning on in-domain data. We observe that different techniques provide varied improvements based on the available data of the language pair. Based on the finding, we integrate these techniques into one training pipeline. For En->Ta, we explore an unconstrained setup with additional Tamil bitext and monolingual data and show that further improvement can be obtained. On the test set, our best submitted systems achieve 21.5 and 13.7 BLEU for Ta->En and En->Ta respectively, and 27.9 and 13.0 for Iu->En and En->Iu respectively.


title: Foundations of Bayesian Learning from Synthetic Data

abstract: There is significant growth and interest in the use of synthetic data as an enabler for machine learning in environments where the release of real data is restricted due to privacy or availability constraints. Despite a large number of methods for synthetic data generation, there are comparatively few results on the statistical properties of models learnt on synthetic data, and fewer still for situations where a researcher wishes to augment real data with another party's synthesised data. We use a Bayesian paradigm to characterise the updating of model parameters when learning in these settings, demonstrating that caution should be taken when applying conventional learning algorithms without appropriate consideration of the synthetic data generating process and learning task. Recent results from general Bayesian updating support a novel and robust approach to Bayesian synthetic-learning founded on decision theory that outperforms standard approaches across repeated experiments on supervised learning and inference problems.


title: Strict hierarchy between parallel, sequential, and indefinite-causal-order strategies for channel discrimination

abstract: We present an instance of a task of minimum-error discrimination of two qubit-qubit quantum channels for which a sequential strategy outperforms any parallel strategy. We then establish two new classes of strategies for channel discrimination that involve indefinite causal order and show that there exists a strict hierarchy among the performance of all four strategies. Our proof technique employs a general method of computer-assisted proofs. We also provide a systematic method for finding pairs of channels that showcase this phenomenon, demonstrating that the hierarchy between the strategies is not exclusive to our main example.


title: ABCD of Kondo effect

abstract: We explore the Kondo effect incorporating the localized impurity transforming under generic symmetry group $G$, that we call the $G$-Kondo effect. We derive the one-dimensional effective model coupled with the impurity, and studied the thermodynamic properties based on the boundary conformal field theory approach. We in particular study the impurity entropy at the conformal fixed point, and the low temperature behavior of the specific heat and the susceptibility. We also consider the Wilson ratio based on these results, and mention the case with exceptional group symmetry.


title: Detecting Receptivity for mHealth Interventions in the Natural Environment

abstract: JITAI is an emerging technique with great potential to support health behavior by providing the right type and amount of support at the right time. A crucial aspect of JITAIs is properly timing the delivery of interventions, to ensure that a user is receptive and ready to process and use the support provided. Some prior works have explored the association of context and some user-specific traits on receptivity, and have built post-study machine-learning models to detect receptivity. For effective intervention delivery, however, a JITAI system needs to make in-the-moment decisions about a user's receptivity. To this end, we conducted a study in which we deployed machine-learning models to detect receptivity in the natural environment, i.e., in free-living conditions. We leveraged prior work regarding receptivity to JITAIs and deployed a chatbot-based digital coach~-- Ally~-- that provided physical-activity interventions and motivated participants to achieve their step goals. We extended the original Ally~app to include two types of machine-learning model that used contextual information about a person to predict when a person is receptive: a \textit{static model\/} that was built before the study started and remained constant for all participants and an \textit{adaptive model\/} that continuously learned the receptivity of individual participants and updated itself as the study progressed. For comparison, we included a \textit{control model\/} that sent intervention messages at random times. The app randomly selected a delivery model for each intervention message. We observed that the machine-learning models led up to a 40\% improvement in receptivity as compared to the control model. Further, we evaluated the temporal dynamics of the different models and observed that receptivity to messages from the adaptive model increased over the course of the study.


title: Asymptotic Rate Analysis for Impairments-Aware Multi-Carrier FD Massive MIMO Relay Networks utilizing MRT/MRC Strategy

abstract: In this paper, we analyze the asymptotic rate for a multi-carrier (MC) full-duplex (FD) massive multiple input multiple output (mMIMO) decode and forward (DF) relay system which serves multiple MC single-antenna half-duplex (HD) nodes. We take into account the impact of hardware distortions resulting in residual self-interference (SI) and inter-carrier leakage (ICL) as well as the impact of imperfect channel state information (CSI). We derive the asymptotic rate expression of our system employed with maximum ratio transmitting (MRT)/ maximum ratio combining (MRC) strategy when the number of the antenna becomes large (goes to $\infty$). It is noticed that the impact of hardware distortion becomes remarkable in a large-scale antenna array regime. On contrary to the effect of multi-user interference and receiver noise, which vanishes as the number of relay antenna goes to infinity, the residual SI and ICL caused by the hardware impairments remains in the MC system.


title: An Interpretable Machine Learning Model for Deformation of Multi-Walled Carbon Nanotubes

abstract: We present a novel interpretable machine learning model to accurately predict complex rippling deformations of Multi-Walled Carbon Nanotubes(MWCNTs) made of millions of atoms. Atomistic-physics-based models are accurate but computationally prohibitive for such large systems. To overcome this bottleneck, we have developed a machine learning model that consists of a novel dimensionality reduction technique and a deep neural network-based learning in the reduced dimension. The proposed nonlinear dimensionality reduction technique extends the functional principal component analysis to satisfy the constraint of deformation. Its novelty lies in designing a function space that satisfies the constraint exactly, which is crucial for efficient dimensionality reduction. Owing to the dimensionality reduction and several other strategies adopted in the present work, learning through deep neural networks is remarkably accurate. The proposed model accurately matches an atomistic-physics-based model while being orders of magnitude faster. It extracts universally dominant patterns of deformation in an unsupervised manner. These patterns are comprehensible and elucidate how the model predicts, yielding interpretability. The proposed model can form a basis for the exploration of machine learning toward the mechanics of one and two-dimensional materials.


title: Dynamically Tagged Groups of Very Metal-poor Halo Stars from the HK and Hamburg/ESO Surveys

abstract: We analyze the dynamical properties of $\sim$1500 very metal-poor (VMP; [Fe/H] $\lesssim -2.0$) halo stars, based primarily on medium-resolution spectroscopic data from the HK and Hamburg/ESO surveys. These data, collected over the past thirty years, are supplemented by a number of calibration stars and other small samples, along with astrometric information from $Gaia$ DR2. We apply a clustering algorithm to the 4-D energy-action space of the sample, and identify a set of 38 Dynamically Tagged Groups (DTGs), containing between 5 and 30 member stars. Many of these DTGs can be associated with previously known prominent substructures such as $Gaia$-Sausage/Enceladus (GSE), Sequoia, the Helmi Stream (HStr), and Thamnos. Others are associated with previously identified smaller dynamical groups of stars and streams. We identify 10 new DTGs as well, many of which have strongly retrograde orbits. We also investigate possible connections between our DTGs and $\sim$300 individual $r$-process-enhanced (RPE) stars from a recent literature compilation. We find that several of these objects have similar dynamical properties to GSE (5), the HStr (4), Sequoia (1), and Rg5 (1), indicating that their progenitors might have been important sources of RPE stars in the Galaxy. Additionally, a number of our newly identified DTGs are shown to be associated with at least two RPE stars each (DTG-2: 3, DTG-7: 2; DTG-27: 2). Taken as a whole, these results are consistent with ultra-faint and/or dwarf spheroidal galaxies as birth environments in which $r$-process nucleosynthesis took place, and then were disrupted by the Milky Way.


abstract: Deep Subspace Clustering Networks (DSC) provide an efficient solution to the problem of unsupervised subspace clustering by using an undercomplete deep auto-encoder with a fully-connected layer to exploit the self expressiveness property. This method uses undercomplete representations of the input data which makes it not so robust and more dependent on pre-training. To overcome this, we propose a simple yet efficient alternative method - Overcomplete Deep Subspace Clustering Networks (ODSC) where we use overcomplete representations for subspace clustering. In our proposed method, we fuse the features from both undercomplete and overcomplete auto-encoder networks before passing them through the self-expressive layer thus enabling us to extract a more meaningful and robust representation of the input data for clustering. Experimental results on four benchmark datasets show the effectiveness of the proposed method over DSC and other clustering methods in terms of clustering error. Our method is also not as dependent as DSC is on where pre-training should be stopped to get the best performance and is also more robust to noise. Code - \href{https://github.com/jeya-maria-jose/Overcomplete-Deep-Subspace-Clustering}{https://github.com/jeya-maria-jose/Overcomplete-Deep-Subspace-Clustering


title: Rotational bands beyond the Elliott model

abstract: Rotational bands are commonplace in the spectra of atomic nuclei. Inspired by early descriptions of these bands by quadrupole deformations of a liquid drop, Elliott constructed a discrete nucleon representations of $\mathrm{SU}(3)$ from fermionic creation and annihilation operators. Ever since, Elliott's model has been foundational to descriptions of rotation in nuclei. Later work, however, suggested the symplectic extension $\mathrm{Sp}(3,R)$ provides a more unified picture. We decompose no-core shell-model nuclear wave functions into symmetry-defined subspaces for several beryllium isotopes, as well as $^{20}$Ne, using the quadratic Casimirs of both Elliott's $\mathrm{SU}(3)$ and $\mathrm{Sp}(3,R)$. The band structure, delineated by strong $B(E2)$ values, has a more consistent description in $\mathrm{Sp}(3,R)$ rather than $\mathrm{SU}(3)$. {In particular, we confirm previous work finding in some nuclides strongly connected upper and lower bands with the same underlying symplectic structure.


title: Impacts of density-dependent predation, cannibalism and fishing in a stage-structured population model of the blue crab in Chesapeake Bay

abstract: The blue crab (Callinectes sapidus) is a dominant ecological species of high commercial value. Spawning stock and recruitment of the Chesapeake Bay population declined by 80% in the 1990s. After severe management actions were implemented in 2008, female abundance rebounded to pre-1994 levels and stabilized. The stepwise decline in the early 1990s, followed by a consistently low level of abundance for 15 y and a jump to high abundance after 2008, suggested the existence of alternative stable states. Alternatively, high fishing pressure combined with low recruitment in 1992 could have triggered a proportional decline in the population, followed by a population increase in 2008 due to rigorous management actions that reduced fishing. We evaluated these alternatives with a stage-structured dynamic population model using ordinary differential equations. In addition, stock assessment models assume that fishing and mortality are independent of density. Hence, we also investigated the role of density-dependent predation, cannibalism and fishing in blue crab population dynamics. We conclude that for the blue crab population in Chesapeake Bay: (1) bistable positive states are not likely with biologically realistic parameter values; (2) hyperbolic (depensatory) fishing will not produce extinction at the range of population densities observed in the bay; and (3) crabs can survive a higher fishing rate under the more realistic assumption of sigmoidal (density-dependent) predation and cannibalism than under constant (density-independent) predation and cannibalism. These collectively indicate that the blue crab population in Chesapeake Bay is resilient to a range of biotic and abiotic disturbances.


title: Change point detection for COVID-19 excess deaths in Belgium

abstract: Emerging at the end of 2019, COVID-19 has become a public health threat to people worldwide. Apart from the deaths who tested positive for COVID-19, many others have died from causes indirectly related to COVID-19. Therefore, the COVID-19 confirmed deaths underestimate the influence of the pandemic on the society; instead, the measure of `excess deaths' is a more objective and comparable way to assess the scale of the epidemic and formulate lessons. One common practical issue in analyzing the impact of COVID-19 is to determine the `pre-COVID-19' period and the `post-COVID-19' period. We apply a change point detection method to identify any change points using the excess deaths in Belgium.


title: New radial velocity observations of AH Her: evidence for material outside the tidal radius

abstract: Spectroscopic observations of AH Herculis during a deep quiescent state are put forward. We found the object in a rare long minima, allowing us to derive accurately the semi-amplitudes: $K_1 =121 \pm \, 4$ \kms\ and $K_2 =152 \pm 2$ \kms\ and its mass functions $M_W{ \sin }^{ 3 }i=0.30 \pm 0.01$ M$_{\odot}$ and $M_R{ \sin }^{ 3 }i=0.24 \pm 0.02$ M$_{\odot}$, while its binary separation is given by $a \sin i =1.39 \pm 0.02$~R$_{\odot}$. The orbital period $P_{orb}$~=~ 0.25812~$\pm~0.00032$~days was found from a power spectrum analysis of the radial velocities of the secondary star. These values are consistent with those determined by \citet{Horne:1986}. Our observations indicate that K5 is the most likely spectral type of the secondary. We discuss why we favour the assumption that the donor in AH Her is a slightly evolved star, in which case we find that the best solution for the inclination yields $i = 48^\circ \pm 2^\circ$. Nonetheless, should the donor be a ZAMS star, we obtain that the inclination is between $ i = 43^\circ$ and $i = 44^\circ$. We also present Doppler tomography of H$\alpha$ and H$\beta$, and found that the emission in both lines is concentrated in a large asymmetric region at low velocities, but at an opposite position to the secondary star, outside the tidal radius and therefore at an unstable position. We also analyse the H$\alpha$ and H$\beta$ line profiles, which show a single broad peak and compare it with the previous quiescent state study which shows a double-peaked profile, providing evidence for its transient nature.


title: Educational robotics for children and their teachers

abstract: This paper describes a Google Educator funded project devoted to the training of teachers (primary and secondary school) through an e-learning platform that will introduce them to educational robotics using Wolly, a social, educational and affective robot.


title: Sparse handlebody decompositions and non-finiteness of $g_3=0$

abstract: We prove that a PL manifold admits a handle decomposition into handles of index $\le k$ if and only if $M$ is $k$-stacked, i.e., it admits a PL triangulation in which all $(d-k-1)$-faces are on $\partial M$. We use this to solve a problem posed in 2008 by Kalai: In any dimension higher than four, there are infinitely many homology-spheres with $g_3 =0$.


title: IR-improved DGLAP parton shower effects for associated production of a W boson and jets in pp collisions at $\sqrt{s}=$8 and 13 TeV

abstract: In a previous paper, hereafter referred to as I, we have analyzed the 7 TeV LHC data on W + jets events from the standpoint of IR-improved DGLAP parton shower effects, using the IR-improved Herwiri1.031 parton shower MC in comparison with the Herwig6.5 parton shower MC in the context of the exact $O(\alpha_s)$ matrix element matched parton shower framework provided by MG5\_aMC@NLO. In the current paper, we extend this analysis to the LHC 8 and 13 TeV data to investigate the energy dependence of the results obtained in I. Specifically, W~+ jet events are generated in the MADGRAPH5\_aMC@NLO framework and showered by HERWIG6.521 and HERWIRI1.031 with $\mathtt{PTRMS}=2.2$ and 0 GeV, respectively. The differential cross sections are reported as functions of jet multiplicity, transverse linear momenta ($P_{T}$), the jet pesudo-rapidity ($\eta$) and the scalar sum of jet transverse momenta ($H_{T}$) for different jet multiplicities 1--3. The dijet cross sections as functions of transverse linear momenta, invariant mass of the dijet and the jet separation are shown as well. Distributions of angular correlations between the jets and the muon are examined as well and the corresponding cross sections are presented. The respective measured cross sections are compared with the exact next-to-leading-order (NLO) matrix element matched parton shower theoretical predictions provided by MADGRAPH5\_aMC@NLO/HERWIRI1.031~($\mathtt{PTRMS}=0)$ and MADGRAPH5\_aMC@NLO/HERWIG6.521~($\mathtt{PTRMS}=2.2~\mathrm{GeV})$ and the phenomenological consequences are discussed with an eye toward their energy dependence.


title: Quantifying information via Shannon entropy in spatially structured optical beams

abstract: While information is ubiquitously generated, shared, and analyzed in a modern-day life, there is still some controversy around the ways to asses the amount and quality of information inside a noisy optical channel. A number of theoretical approaches based on, e.g., conditional Shannon entropy and Fisher information have been developed, along with some experimental validations. Some of these approaches are limited to a certain alphabet, while others tend to fall short when considering optical beams with non-trivial structure, such as Hermite-Gauss, Laguerre-Gauss and other modes with non-trivial structure. Here, we propose a new definition of classical Shannon information via the Wigner distribution function, while respecting the Heisenberg inequality. Following this definition, we calculate the amount of information in a Gaussian, Hermite-Gaussian, and Laguerre-Gaussian laser modes in juxtaposition and experimentally validate it by reconstruction of the Wigner distribution function from the intensity distribution of structured laser beams. We experimentally demonstrate the technique that allows to infer field structure of the laser beams in singular optics to assess the amount of contained information. Given the generality, this approach of defining information via analyzing the beam complexity is applicable to laser modes of any topology that can be described by 'well-behaved' functions. Classical Shannon information defined in this way is detached from a particular alphabet, i.e. communication scheme, and scales with the structural complexity of the system. Such a synergy between the Wigner distribution function encompassing the information in both real and reciprocal space, and information being a measure of disorder, can contribute into future coherent detection algorithms and remote sensing.


title: Anonymizing Sensor Data on the Edge: A Representation Learning and Transformation Approach

abstract: The abundance of data collected by sensors in Internet of Things (IoT) devices, and the success of deep neural networks in uncovering hidden patterns in time series data have led to mounting privacy concerns. This is because private and sensitive information can be potentially learned from sensor data by applications that have access to this data. In this paper, we aim to examine the tradeoff between utility and privacy loss by learning low-dimensional representations that are useful for data obfuscation. We propose deterministic and probabilistic transformations in the latent space of a variational autoencoder to synthesize time series data such that intrusive inferences are prevented while desired inferences can still be made with sufficient accuracy. In the deterministic case, we use a linear transformation to move the representation of input data in the latent space such that the reconstructed data is likely to have the same public attribute but a different private attribute than the original input data. In the probabilistic case, we apply the linear transformation to the latent representation of input data with some probability. We compare our technique with autoencoder-based anonymization techniques and additionally show that it can anonymize data in real time on resource-constrained edge devices.


title: Perturbation theory of the quadratic Lotka-Volterra double center

abstract: We revisit the bifurcation theory of the Lotka-Volterra quadratic system \begin{eqnarray} X_0 :\left\{\begin{aligned} \dot{x}=& - y -x^2+y^2 ,\\ \dot{y}= &\;\;\;\;x - 2xy \end{aligned} \right. \end{eqnarray} with respect to arbitrary quadratic deformations. The system $X_0$ has a double center, which is moreover isochronous. We show that the deformed system $X_0$ can have at most two limit cycles on the finite plane, with possible distribution $(i,j)$, where $i+j\leq2$. Our approach is based on the study of pairs of bifurcation functions associated to the centers, expressed in terms of iterated path integrals of length two.


title: Feature Sharing and Integration for Cooperative Cognition and Perception with Volumetric Sensors

abstract: The recent advancement in computational and communication systems has led to the introduction of high-performing neural networks and high-speed wireless vehicular communication networks. As a result, new technologies such as cooperative perception and cognition have emerged, addressing the inherent limitations of sensory devices by providing solutions for the detection of partially occluded targets and expanding the sensing range. However, designing a reliable cooperative cognition or perception system requires addressing the challenges caused by limited network resources and discrepancies between the data shared by different sources. In this paper, we examine the requirements, limitations, and performance of different cooperative perception techniques, and present an in-depth analysis of the notion of Deep Feature Sharing (DFS). We explore different cooperative object detection designs and evaluate their performance in terms of average precision. We use the Volony dataset for our experimental study. The results confirm that the DFS methods are significantly less sensitive to the localization error caused by GPS noise. Furthermore, the results attest that detection gain of DFS methods caused by adding more cooperative participants in the scenes is comparable to raw information sharing technique while DFS enables flexibility in design toward satisfying communication requirements.


title: Mode Localization and Suppressed Heat Transport in Amorphous Alloys

abstract: Glasses usually represent the lower limit for the thermal conductivity of solids, but a fundamental understanding of lattice heat transport in amorphous materials can provide design rules to beat such a limit. Here we investigate the role of mass disorder in glasses by studying amorphous silicon-germanium alloy (a-Si$_{1-x}$Ge$_x$) over the full range of atomic concentration from $x=0$ to $x=1$, using molecular dynamics and the quasi-harmonic Green-Kubo lattice dynamics formalism. We find that the thermal conductivity of a-Si$_{1-x}$Ge$_x$ as a function of $x$ exhibits a smoother U-shape than in crystalline mass-disordered alloys. The main contribution to the initial drop of thermal conductivity at low Ge concentration stems from the localization of otherwise extended modes that make up the lowest 8\% of the population by frequency. Contributions from intermediate frequency modes are decreased more gradually with increasing Ge to reach a broad minimum thermal conductivity between concentrations of Ge from $x=0.25$ to $0.75$.


title: Multi-Step Recurrent Q-Learning for Robotic Velcro Peeling

abstract: Learning object manipulation is a critical skill for robots to interact with their environment. Even though there has been significant progress in robotic manipulation of rigid objects, interacting with non-rigid objects remains challenging for robots. In this work, we introduce velcro peeling as a representative application for robotic manipulation of non-rigid objects in complex environments. We present a method of learning force-based manipulation from noisy and incomplete sensor inputs in partially observable environments by modeling long term dependencies between measurements with a multi-step deep recurrent network. We present experiments on a real robot to show the necessity of modeling these long term dependencies and validate our approach in simulation and robot experiments. Our results show that using tactile input enables the robot to overcome geometric uncertainties present in the environment with high fidelity in ~90% of all cases, outperforming the baselines by a large margin.