title: The AGN Fraction in Dwarf Galaxies from eROSITA: First Results and Future Prospects

abstract: Determining the fraction of nearby dwarf galaxies hosting massive black holes (BHs) can inform our understanding of the origin of "seed" black holes at high redshift. Here we search for signatures of accreting massive BHs in a sample of dwarf galaxies ($M_\star \le 3 \times 10^9~M_\odot$, $z \leq 0.15$) selected from the NASA-Sloan Atlas (NSA) using X-ray observations from the eROSITA Final Equatorial Depth Survey (eFEDS). On average, our search is sensitive to active galactic nuclei (AGNs) in dwarf galaxies that are accreting at $\gtrsim 1%$ of their Eddington luminosity. Of the ${\sim}28,000$ X-ray sources in eFEDS and the 495 dwarf galaxies in the NSA within the eFEDS footprint, we find six galaxies hosting possible active massive BHs. If the X-ray sources are indeed associated with the dwarf galaxies, the X-ray emission is above that expected from star formation, with X-ray source luminosities of $L_{0.5-8~\textrm{keV}} \sim 10^{39\textrm{-}40}$ erg s$^{-1}$. Additionally, after accounting for chance alignments of background AGNs with dwarf galaxies, we estimate there are between 0-9 real associations between dwarf galaxies and X-ray sources in the eFEDS field at the 95% confidence level. From this we find an upper limit on the eFEDS-detected dwarf galaxy AGN fraction of $\le 1.8%$, which is broadly consistent with similar studies at other wavelengths. We extrapolate these findings from the eFEDS sky coverage to the planned eROSITA All-Sky Survey and estimate that upon completion, the all-sky survey could yield as many as ${\sim}1350$ AGN candidates in dwarf galaxies at low redshift.


title: Sandpile groups of supersingular isogeny graphs

abstract: Let $p$ and $q$ be distinct primes, and let $X_{p,q}$ be the $(q+1)$-regular graph whose nodes are supersingular elliptic curves over $\overline{\mathbb{F}}_p$ and whose edges are $q$-isogenies. For fixed $p$, we compute the distribution of the $\ell$-Sylow subgroup of the sandpile group (i.e.\ Jacobian) of $X_{p,q}$ as $q \to \infty$. We find that the distribution disagrees with the Cohen-Lenstra heuristic in this context. Our proof is via Galois representations attached to modular curves. As a corollary of our result, we give an upper bound on the probability that the Jacobian is cyclic, which we conjecture to be sharp.


title: The ComMA Dataset V0.2: Annotating Aggression and Bias in Multilingual Social Media Discourse

abstract: In this paper, we discuss the development of a multilingual dataset annotated with a hierarchical, fine-grained tagset marking different types of aggression and the "context" in which they occur. The context, here, is defined by the conversational thread in which a specific comment occurs and also the "type" of discursive role that the comment is performing with respect to the previous comment. The initial dataset, being discussed here (and made available as part of the ComMA@ICON shared task), consists of a total 15,000 annotated comments in four languages - Meitei, Bangla, Hindi, and Indian English - collected from various social media platforms such as YouTube, Facebook, Twitter and Telegram. As is usual on social media websites, a large number of these comments are multilingual, mostly code-mixed with English. The paper gives a detailed description of the tagset being used for annotation and also the process of developing a multi-label, fine-grained tagset that can be used for marking comments with aggression and bias of various kinds including gender bias, religious intolerance (called communal bias in the tagset), class/caste bias and ethnic/racial bias. We also define and discuss the tags that have been used for marking different the discursive role being performed through the comments, such as attack, defend, etc. We also present a statistical analysis of the dataset as well as results of our baseline experiments with developing an automatic aggression identification system using the dataset developed.


title: Data Excellence for AI: Why Should You Care

abstract: The efficacy of machine learning (ML) models depends on both algorithms and data. Training data defines what we want our models to learn, and testing data provides the means by which their empirical progress is measured. Benchmark datasets define the entire world within which models exist and operate, yet research continues to focus on critiquing and improving the algorithmic aspect of the models rather than critiquing and improving the data with which our models operate. If "data is the new oil," we are still missing work on the refineries by which the data itself could be optimized for more effective use.


title: Perturbatively Exact $w_{1+\infty}$ Asymptotic Symmetry of Quantum Self-Dual Gravity

abstract: The infinite tower of positive-helicity soft gravitons in any minimally coupled, tree-level, asymptotically flat four-dimensional (4D) gravity was recently shown to generate a $w_{1+\infty}$ asymptotic symmetry algebra. It is natural to ask whether this classical algebra acquires quantum corrections at loop level. We explore this in quantum self-dual gravity, whose amplitudes acquire known one-loop exact all-plus helicity quantum corrections. We show using collinear splitting formulae that, remarkably, the $w_{1+\infty}$ algebra persists in quantum self-dual gravity without corrections.


title: The $r$-coloring and maximum stable set problem in hypergraphs with bounded matching number and edge size

abstract: Motivated by the analogous questions in graphs, we study the complexity of coloring and stable set problems in hypergraphs with forbidden substructures and bounded edge size. Letting $\nu(G)$ denote the maximum size of a matching in $H$, we obtain complete dichotomies for the complexity of the following problems parametrized by fixed $r, k, s \in \mathbb{N}$: $r$-Coloring in hypergraphs $G$ with edge size at most $k$ and $\nu(G) \leq s$; $r$-Precoloring Extension in $k$-uniform hypergraphs $G$ with $\nu(G) \leq s$; $r$-Precoloring Extension in hypergraphs $G$ with edge size at most $k$ and $\nu(G) \leq s$; Maximum Stable Set in $k$-uniform hypergraphs $G$ with $\nu(G) \leq s$; Maximum Weight Stable Set in $k$-uniform hypergraphs with $\nu(G) \leq s$; as well as partial results for $r$-Coloring in $k$-uniform hypergraphs $\nu(G) \leq s$. We then turn our attention to $2$-Coloring in 3-uniform hypergraphs with forbidden induced subhypergraphs, and give a polynomial-time algorithm when restricting the input to hypergraphs excluding a fixed one-edge hypergraph. Finally, we consider linear 3-uniform hypergraphs (in which every two edges share at most one vertex), and show that excluding an induced matching in $G$ implies that $\nu(G)$ is bounded by a constant; and that $3$-coloring linear $3$-uniform hypergraphs $G$ with $\nu(G) \leq 532$ is NP-hard.


title: A template method to measure the $t \bar t$ polarisation

abstract: We develop a template method for the measurement of the polarisation of $t \bar t$ pairs produced in hadron collisions. The method would allow to extract the individual fractions of $t_L \bar t_L$, $t_R \bar t_R$, $t_L \bar t_R$ and $t_R \bar t_L$ pairs with a fit to data, where $L,R$ refer to the polarisation along any axis. These polarisation fractions have not been independently measured at present. Secondarily, the method also provides the net polarisation of $t$ and $\bar t$, as well as their spin correlation for arbitrary axes.


title: Very high-energy collective states of partons in fractional quantum Hall liquids

abstract: The low energy physics of fractional quantum Hall (FQH) states -- a paradigm of strongly correlated topological phases of matter -- to a large extent is captured by weakly interacting quasiparticles known as composite fermions (CFs). In this paper, based on numerical simulations and effective field theory, we argue that some \emph{high energy} states in the FQH spectra necessitate a different description based on \emph{parton} quasiparticles. We show that Jain states at filling factor $\nu{=}n/(2pn\pm1)$ with integers $n,p{\geq}2$, support two kinds of collective modes: in addition to the well-known Girvin-MacDonald-Platzman (GMP) mode, they host a high energy collective mode, which is interpreted as the GMP mode of partons. We elucidate observable signatures of the parton mode in the dynamics following a geometric quench. We construct a microscopic wave function for the parton mode, and demonstrate agreement between its variational energy and exact diagonalization. Using the parton construction, we derive a field theory of the Jain states and show that the previously proposed effective theories follow from our approach. Our results point to partons being "real" quasiparticles which, in a way reminiscent of quarks, only become observable at sufficiently high energies.


title: Evidence for Nonlinear Gluon Effects in QCD and their $A$ Dependence at STAR

abstract: The STAR Collaboration reports measurements of back-to-back azimuthal correlations of di-$\pi^0$s produced at forward pseudorapidities ($2.6<\eta<4.0$) in $p$+$p$, $p+$Al, and $p+$Au collisions at a center-of-mass energy of 200 GeV. We observe a clear suppression of the correlated yields of back-to-back $\pi^0$ pairs in $p+$Al and $p+$Au collisions compared to the $p$+$p$ data. The observed suppression of back-to-back pairs as a function of transverse momentum suggests nonlinear gluon dynamics arising at high parton densities. The larger suppression found in $p+$Au relative to $p+$Al collisions exhibits a dependence of the saturation scale, $Q_s^2$, on the mass number, $A$. A linear scaling of the suppression with $A^{1/3}$ is observed with a slope of $-0.09$ $\pm$ $0.01$.


title: A Radon Transform on the Cylinder

abstract: We define a parametric Radon transform $R$ that assigns to a Sobolev function on the cylinder $\mathbb{S}\times \mathbb{R}$ in $\mathbb{R}^3$ its mean values along sets $E_\zeta$ formed by the intersections of planes through the origin and the cylinder. We show that $R$ is a continuous operator, prove an inversion formula, provide a support theorem, as well as a characterization of its null space. We conclude by presenting a formula for the dual transform $R^*$. We show that $R$ and its dual $R^*$ are related to the right-sided and left-sided Chebyshev fractional integrals. Using this relationship, we characterize the null space of $R$ and $R^*$ and provide an inversion formula for $R^*$.


title: Fast Discovery of Nested Dependencies on JSON Data

abstract: Functional and inclusion dependencies are the most widely used classes of data dependencies in data profiling due to their ability to identify relationships in data such as primary and foreign keys. These relationships are equally important when dealing with nested data formats such as JSON. However, the definition of functional and inclusion dependencies makes use of a flat, unnested relational model which misses many useful types of dependencies on data which involve nested data models. In this work, we identify types of dependencies which are not captured by traditional functional and inclusion dependencies but which nevertheless capture meaningful relationships among nested data. We also demonstrate how algorithms for mining these traditional dependencies can be adapted to also mine nested dependencies. The first strategy simply flattens the input data and feeds into unmodified existing algorithms. We present a second strategy which instead adapts the algorithm to efficiently process JSON data as input which in some cases leads to a reduction in runtime by multiple orders of magnitude on real-world datasets. We further show how these algorithms can be adapted to produce useful results in the presence of a percentage of incomplete or invalid data.


title: What Stops Learning-based 3D Registration from Working in the Real World?

abstract: Much progress has been made on the task of learning-based 3D point cloud registration, with existing methods yielding outstanding results on standard benchmarks, such as ModelNet40, even in the partial-to-partial matching scenario. Unfortunately, these methods still struggle in the presence of real data. In this work, we identify the sources of these failures, analyze the reasons behind them, and propose solutions to tackle them. We summarise our findings into a set of guidelines and demonstrate their effectiveness by applying them to different baseline methods, DCP and IDAM. In short, our guidelines improve both their training convergence and testing accuracy. Ultimately, this translates to a best-practice 3D registration network (BPNet), constituting the first learning-based method able to handle previously-unseen objects in real-world data. Despite being trained only on synthetic data, our model generalizes to real data without any fine-tuning, reaching an accuracy of up to 67% on point clouds of unseen objects obtained with a commercial sensor.


title: RacketStore: Measurements of ASO Deception in Google Play via Mobile and App Usage

abstract: Online app search optimization (ASO) platforms that provide bulk installs and fake reviews for paying app developers in order to fraudulently boost their search rank in app stores, were shown to employ diverse and complex strategies that successfully evade state-of-the-art detection methods. In this paper we introduce RacketStore, a platform to collect data from Android devices of participating ASO providers and regular users, on their interactions with apps which they install from the Google Play Store. We present measurements from a study of 943 installs of RacketStore on 803 unique devices controlled by ASO providers and regular users, that consists of 58,362,249 data snapshots collected from these devices, the 12,341 apps installed on them and their 110,511,637 Google Play reviews. We reveal significant differences between ASO providers and regular users in terms of the number and types of user accounts registered on their devices, the number of apps they review, and the intervals between the installation times of apps and their review times. We leverage these insights to introduce features that model the usage of apps and devices, and show that they can train supervised learning algorithms to detect paid app installs and fake reviews with an F1-measure of 99.72% (AUC above 0.99), and detect devices controlled by ASO providers with an F1-measure of 95.29% (AUC = 0.95). We discuss the costs associated with evading detection by our classifiers and also the potential for app stores to use our approach to detect ASO work with privacy.


title: Community-Detection via Hashtag-Graphs for Semi-Supervised NMF Topic Models

abstract: Extracting topics from large collections of unstructured text-documents has become a central task in current NLP applications and algorithms like NMF, LDA as well as their generalizations are the well-established current state of the art. However, especially when it comes to short text documents like Tweets, these approaches often lead to unsatisfying results due to the sparsity of the document-feature matrices. Even though, several approaches have been proposed to overcome this sparsity by taking additional information into account, these are merely focused on the aggregation of similar documents and the estimation of word-co-occurrences. This ultimately completely neglects the fact that a lot of topical-information can be actually retrieved from so-called hashtag-graphs by applying common community detection algorithms. Therefore, this paper outlines a novel approach on how to integrate topic structures of hashtag graphs into the estimation of topic models by connecting graph-based community detection and semi-supervised NMF. By applying this approach on recently streamed Twitter data it will be seen that this procedure actually leads to more intuitive and humanly interpretable topics.


title: Influence Dynamics Among Narratives: A Case Study of the Venezuelan Presidential Crisis

abstract: It is widely understood that diffusion of and simultaneous interactions between narratives -- defined here as persistent point-of-view messaging -- significantly contributes to the shaping of political discourse and public opinion. In this work, we propose a methodology based on Multi-Variate Hawkes Processes and our newly-introduced Process Influence Measures for quantifying and assessing how such narratives influence (Granger-cause) each other. Such an approach may aid social scientists enhance their understanding of socio-geopolitical phenomena as they manifest themselves and evolve in the realm of social media. In order to show its merits, we apply our methodology on Twitter narratives during the 2019 Venezuelan presidential crisis. Our analysis indicates a nuanced, evolving influence structure between 8 distinct narratives, part of which could be explained by landmark historical events.


title: Towards Integrative Multi-Modal Personal Health Navigation Systems: Framework and Application

abstract: It is well understood that an individual's health trajectory is influenced by choices made in each moment, such as from lifestyle or medical decisions. With the advent of modern sensing technologies, individuals have more data and information about themselves than any other time in history. How can we use this data to make the best decisions to keep the health state optimal? We propose a generalized Personal Health Navigation (PHN) framework. PHN takes individuals towards their personal health goals through a system which perpetually digests data streams, estimates current health status, computes the best route through intermediate states utilizing personal models, and guides the best inputs that carry a user towards their goal. In addition to describing the general framework, we test the PHN system in two experiments within the field of cardiology. First, we prospectively test a knowledge-infused cardiovascular PHN system with a pilot clinical trial of 41 users. Second, we build a data-driven personalized model on cardiovascular exercise response variability on a smartwatch data-set of 33,269 real-world users. We conclude with critical challenges in health computing for PHN systems that require deep future investigation.


title: Tweet Analysis for Enhancement of COVID-19 Epidemic Simulation: A Case Study in Japan

abstract: The COVID-19 pandemic, which began in December 2019, progressed in a complicated manner and thus caused problems worldwide. Seeking clues to the reasons for the complicated progression is necessary but challenging in the fight against the pandemic. We sought clues by investigating the relationship between reactions on social media and the COVID-19 epidemic in Japan. Twitter was selected as the social media platform for study because it has a large user base in Japan and because it quickly propagates short topic-focused messages ("tweets"). Analysis using Japanese Twitter data suggests that reactions on social media and the progression of the COVID-19 pandemic may have a close relationship. Experiments to evaluate the potential of using tweets to support the prediction of how an epidemic will progress demonstrated the value of using epidemic-related social media data. Our findings provide insights into the relationship between user reactions on social media, particularly Twitter, and epidemic progression, which can be used to fight pandemics.


title: Optical Vortex Braiding with Bessel Beams

abstract: We propose the braiding of optical vortices in a laser beam with more than 2{\pi} rotation by superposing Bessel modes with a plane wave. We experimentally demonstrate this by using a Bessel-Gaussian beam and a coaxial Gaussian, and we present measurement of three complete braids. The amount of braiding is fundamentally limited only by the numerical aperture of the system and we discuss how braiding can be controlled experimentally for any number of vortices.


title: Exact Convergence Analysis for Metropolis-Hastings Independence Samplers in Wasserstein Distances

abstract: Under mild assumptions, we show the exact convergence rate in total variation is also exact in weaker Wasserstein distances for the Metropolis-Hastings independence sampler. We develop a new upper and lower bound on the worst-case Wasserstein distance when initialized from points. For an arbitrary point initialization, we show the convergence rate is the same and matches the convergence rate in total variation. We derive exact convergence expressions for more general Wasserstein distances when initialization is at a specific point. Using optimization, we construct a novel centered independent proposal to develop exact convergence rates in Bayesian quantile regression and many generalized linear model settings. We show the exact convergence rate can be upper bounded in Bayesian binary response regression (e.g. logistic and probit) when the sample size and dimension grow together.


title: Taming Triangulation Dependence of T6/Z2xZ2 Resolutions

abstract: Resolutions of certain toroidal orbifolds, like T6/Z2xZ2, are far from unique, due to triangulation dependence of their resolved local singularities. This leads to an explosion of the number of topologically distinct smooth geometries associated to a single orbifold. By introducing a parameterisation to keep track of the triangulations used at all resolved singularities simultaneously, (self-)intersection numbers and integrated Chern classes can be determined for any triangulation configuration. Using this method the consistency conditions of line bundle models and the resulting chiral spectra can be worked out for any choice of triangulation. Moreover, by superimposing the Bianchi identities for all triangulation options much simpler though stronger conditions are uncovered. When these are satisfied, flop--transitions between all different triangulations are admissible. Various methods are exemplified by a number of concrete models on resolutions of the T6/Z2xZ2 orbifold.


title: Scalar correlators and normal modes in holographic neutron stars

abstract: We study a scalar field in the background of a holographic neutron star at finite temperature and analyze its asymptotic behavior to compute the two-point correlator at the boundary. From the normal modes of the field we determine the resonances and decay constants of the boundary field theory. We show that the correlator is dominated by the normal modes in the stable regions of phase space of the neutron star, becoming a power law as we move into the unstable zones.


title: The Acrobatics of BQP

abstract: One can fix the randomness used by a randomized algorithm, but there is no analogous notion of fixing the quantumness used by a quantum algorithm. Underscoring this fundamental difference, we show that, in the black-box setting, the behavior of quantum polynomial-time ($\mathsf{BQP}$) can be remarkably decoupled from that of classical complexity classes like $\mathsf{NP}$. Specifically: -There exists an oracle relative to which $\mathsf{NP^{BQP}}\not\subset\mathsf{BQP^{PH}}$, resolving a 2005 problem of Fortnow. As a corollary, there exists an oracle relative to which $\mathsf{P}=\mathsf{NP}$ but $\mathsf{BQP}\neq\mathsf{QCMA}$. -Conversely, there exists an oracle relative to which $\mathsf{BQP^{NP}}\not\subset\mathsf{PH^{BQP}}$. -Relative to a random oracle, $\mathsf{PP}=\mathsf{PostBQP}$ is not contained in the "$\mathsf{QMA}$ hierarchy" $\mathsf{QMA}^{\mathsf{QMA}^{\mathsf{QMA}^{\cdots}}}$. -Relative to a random oracle, $\mathsf{\Sigma}_{k+1}^\mathsf{P}\not\subset\mathsf{BQP}^{\mathsf{\Sigma}_{k}^\mathsf{P}}$ for every $k$. -There exists an oracle relative to which $\mathsf{BQP}=\mathsf{P^{\# P}}$ and yet $\mathsf{PH}$ is infinite. -There exists an oracle relative to which $\mathsf{P}=\mathsf{NP}\neq\mathsf{BQP}=\mathsf{P^{\# P}}$. To achieve these results, we build on the 2018 achievement by Raz and Tal of an oracle relative to which $\mathsf{BQP}\not \subset \mathsf{PH}$, and associated results about the Forrelation problem. We also introduce new tools that might be of independent interest. These include a "quantum-aware" version of the random restriction method, a concentration theorem for the block sensitivity of $\mathsf{AC^0}$ circuits, and a (provable) analogue of the Aaronson-Ambainis Conjecture for sparse oracles.


title: Unconventional Flat Chern Bands and 2$e$ Charges in Skyrmionic Moir\'e Superlattices

abstract: The interplay of topological characteristics in real space and reciprocal space can lead to the emergence of unconventional topological phases. In this Letter, we implement a novel mechanism for generating higher-Chern flat bands on the basis of twisted bilayer graphene (TBG) coupled to topological magnetic structures in the form of the skyrmion lattice. In particular, we discover a scenario for generating $|C|=2$ dispersionless electronic bands when the skyrmion periodicity and the moir\'e periodicity are matched. Following the Wilczek argument, the statistics of the charge-carrying excitations in this case is \textit{bosonic}, characterized by electronic charge $Q =2e$, that is \textit{even} in units of electron charge $e$. The required skyrmion coupling strength triggering the topological phase transition is realistic, with its threshold estimated as low as 4~meV. The Hofstadter butterfly spectrum of this phase is different resulting in an unexpected quantum Hall conductance sequence $\pm \frac{2 e^2}{h}, \ \pm \frac{4 e^2}{h}, ...$ for TBG with skyrmion order.


title: A Transient Semantics for Typed Racket

abstract: Mixed-typed languages enable programmers to link typed and untyped components in various ways. Some offer rich type systems to facilitate the smooth migration of untyped code to the typed world; others merely provide a convenient form of type Dynamic together with a conventional structural type system. Orthogonal to this dimension, Natural systems ensure the integrity of types with a sophisticated contract system, while Transient systems insert simple first-order checks at strategic places within typed code. Furthermore, each method of ensuring type integrity comes with its own blame-assignment strategy. Typed Racket has a rich migratory type system and enforces the types with a Natural semantics. Reticulated Python has a simple structural type system extended with Dynamic and enforces types with a Transient semantics. While Typed Racket satisfies the most stringent gradual-type soundness properties at a significant performance cost, Reticulated Python seems to limit the performance penalty to a tolerable degree and is nevertheless type sound. This comparison raises the question of whether Transient checking is applicable to and beneficial for a rich migratory type system. This paper reports on the surprising difficulties of adapting the Transient semantics of Reticulated Python to the rich migratory type system of Typed Racket. The resulting implementation, Shallow Typed Racket, is faster than the standard Deep Typed Racket but only when the Transient blame assignment strategy is disabled. For language designers, this report provides valuable hints on how to equip an existing compiler to support a Transient semantics. For theoreticians, the negative experience with Transient blame calls for a thorough investigation of this strategy.


title: Types for Tables: A Language Design Benchmark

abstract: Context: Tables are ubiquitous formats for data. Therefore, techniques for writing correct programs over tables, and debugging incorrect ones, are vital. Our specific focus in this paper is on rich types that articulate the properties of tabular operations. We wish to study both their expressive power and _diagnostic quality_. Inquiry: There is no "standard library" of table operations. As a result, every paper (and project) is free to use its own (sub)set of operations. This makes artifacts very difficult to compare, and it can be hard to tell whether omitted operations were left out by oversight or because they cannot actually be expressed. Furthermore, virtually no papers discuss the quality of type error feedback. Approach: We combed through several existing languages and libraries to create a "standard library" of table operations. Each entry is accompanied by a detailed specification of its "type," expressed independent of (and hence not constrained by) any type language. We also studied and categorized a corpus of (student) program edits that resulted in table-related errors. We used this to generate a suite of erroneous programs. Finally, we adapted the concept of a datasheet to facilitate comparisons of different implementations. Knowledge: Our benchmark creates a common ground to frame work in this area. Language designers who claim to support typed programming over tables have a clear suite against which to demonstrate their system's expressive power. Our family of errors also gives them a chance to demonstrate the quality of feedback. Researchers who improve one aspect -- especially error reporting -- without changing the other can demonstrate their improvement, as can those who engage in trade-offs between the two. The net result should be much better science in both expressiveness and diagnostics. We also introduce a datasheet format for presenting this knowledge in a methodical way. Grounding: We have generated our benchmark from real languages, libraries, and programs, as well as personal experience conducting and teaching data science. We have drawn on experience in engineering and, more recently, in data science to generate the datasheet. Importance: Claims about type support for tabular programming are hard to evaluate. However, tabular programming is ubiquitous, and the expressive power of type systems keeps growing. Our benchmark and datasheet can help lead to more orderly science. It also benefits programmers trying to choose a language.


title: Continuation-Passing Style, Defunctionalization, Accumulations, and Associativity

abstract: Context: Reynolds showed us how to use continuation-passing style and defunctionalization to transform a recursive interpreter for a language into an abstract machine for programs in that language. The same techniques explain other programming tricks, including zippers and accumulating parameters. Inquiry: Buried within all those applications there is usually a hidden appeal to the algebraic property of associativity. Approach: Our purpose in this paper is to entice associativity out of the shadows and into the limelight. Knowledge: We revisit some well-known applications (factorial, fast reverse, tree flattening, and a compiler for a simple expression language) to spotlight their dependence on associativity. Grounding: We replay developments of these programs through a series of program transformations and data refinements, justified by equational reasoning. Importance: Understanding the crucial role played by associativity clarifies when continuation-passing style and defunctionalization can help and when they cannot, and may prompt other applications of these techniques.


title: Automated, Targeted Testing of Property-Based Testing Predicates

abstract: Context: This work is based on property-based testing (PBT). PBT is an increasingly important form of software testing. Furthermore, it serves as a concrete gateway into the abstract area of formal methods. Specifically, we focus on students learning PBT methods. Inquiry: How well do students do at PBT? Our goal is to assess the quality of the predicates they write as part of PBT. Prior work introduced the idea of decomposing the predicate's property into a conjunction of independent subproperties. Testing the predicate against each subproperty gives a "semantic" understanding of their performance. Approach: The notion of independence of subproperties both seems intuitive and was an important condition in prior work. First, we show that this condition is overly restrictive and might hide valuable information: it both undercounts errors and makes it hard to capture misconceptions. Second, we introduce two forms of automation, one based on PBT tools and the other on SAT-solving, to enable testing of student predicates. Third, we compare the output of these automated tools against manually-constructed tests. Fourth, we also measure the performance of those tools. Finally, we re-assess student performance reported in prior work. Knowledge: We show the difficulty caused by the independent subproperty requirement. We provide insight into how to use automation effectively to assess PBT predicates. In particular, we discuss the steps we had to take to beat human performance. We also provide insight into how to make the automation work efficiently. Finally, we present a much richer account than prior work of how students did. Grounding: Our methods are grounded in mathematical logic. We also make use of well-understood principles of test generation from more formal specifications. This combination ensures the soundness of our work. We use standard methods to measure performance. Importance: As both educators and programmers, we believe PBT is a valuable tool for students to learn, and its importance will only grow as more developers appreciate its value. Effective teaching requires a clear understanding of student knowledge and progress. Our methods enable a rich and automated analysis of student performance on PBT that yields insight into their understanding and can capture misconceptions. We therefore expect these results to be valuable to educators.


title: Rotating Bose gas dynamically entering the lowest Landau level

abstract: Motivated by recent experiments, we model the dynamics of a condensed Bose gas in a rotating anisotropic trap, where the equations of motion are analogous to those of charged particles in a magnetic field. As the rotation rate is ramped from zero to the trapping frequency, the condensate stretches along one direction and is squeezed along another, becoming long and thin. When the trap anisotropy is slowly switched off on a particular timescale, the condensate is left in the lowest Landau level. We use a time dependent variational approach to quantify these dynamics and give intuitive arguments about the structure of the condensate wavefunction. This preparation of a lowest Landau level condensate can be an important first step in realizing bosonic analogs of quantum Hall states.

categories: cond-mat.quant-gas cond-mat.mes-hall cond-mat.str-el physics.atom-ph quant-ph


title: The return of the spin period in DW Cnc and evidence of new high state outbursts

abstract: DW Cnc is an intermediate polar which has previously been observed in both high and low states. Observations of the high state of DW Cnc have previously revealed a spin period at ~ 38.6 min, however observations from the 2018/19 low state showed no evidence of the spin period. We present results from our analysis of 12 s cadence photometric data collected by NGTS of DW Cnc during the high state which began in 2019. Following the previously reported suppression of the spin period signal we identify the return of this signal during the high state, consistent with previous observations of it. We identify this as the restarting of accretion during the high state. We further identified three short outbursts lasting ~ 1 d in DW Cnc with a mean recurrence time of ~ 60 d and an amplitude of ~ 1 mag. These are the first outbursts identified in DW Cnc since 2008. Due to the short nature of these events we identify them not as a result of accretion instabilities but instead either from instabilities originating from the interaction of the magnetorotational instability in the accretion disc and the magnetic field generated by the white dwarf or the result of magnetic gating.


title: The Close AGN Reference Survey (CARS): IFU survey data and the BH mass dependence of long-term AGN variability

abstract: [Abridged] AGN are thought to be intimately connected with their host galaxies through feeding and feedback processes. A spatially resolved multiwavelength survey is required to map the interaction of AGN with their host galaxies on different spatial scales and different phases of the ISM. The goal of CARS is to obtain the necessary spatially resolved multiwavelength observations for an unbiased sample of local unobscured luminous AGN. We present the overall CARS survey design and the associated wide-field optical IFU spectroscopy for all 41 CARS targets at z<0.06 randomly selected from the Hamburg/ESO survey of luminous unobscured AGN. This data set provides the backbone of CARS and allows us to characterize host galaxy morphologies, AGN parameters, precise systemic redshifts, and ionized gas distributions including excitation conditions, kinematics, and metallicities in unprecedented detail. We focus our study on the size of the ENLR which has been traditionally connected to AGN luminosity. Given the large scatter in the ENLR size-luminosity relation, we performed a large parameter search to identify potentially more fundamental relations. Remarkably, we identified the strongest correlation between the maximum projected ENLR size and the black hole mass, consistent with an $R_\mathrm{ENLR,max}\sim M_\mathrm{BH}^{0.5}$ relationship. We interpret the maximum ENLR size as a timescale indicator of a single BH radiative-efficient accretion episode for which we inferred log(t_AGN) = (0.45+- 0.08)log(M_BH)+1.78 using forward modeling. The extrapolation of our inferred relation toward higher BH masses is consistent with an independent lifetime estimate from the HeII proximity zones around luminous AGN at z~3. While our proposed link between the BH mass and AGN lifetime might be a secondary correlation itself or impacted by unknown biases, it has a few relevant implications if confirmed.


title: The Close AGN Reference Survey (CARS): Locating the [O III] wing component in luminous local Type 1 AGN

abstract: [Abridged]The strong asymmetry in the optical [O III]$\lambda$5007 emission line is one of the best signatures of AGN-driven warm (~10$^4$ K) ionized gas outflows on host galaxy scales. While large spectroscopic surveys like SDSS have characterized the kinematics of [O III] for large samples of AGN, estimating the associated energetics requires spatially resolving these outflows with, for example, IFU studies. As part of CARS we obtained spatially-resolved IFU spectroscopy for a representative sample of 39 luminous type 1 AGN at 0.01<z<0.06 with MUSE and VIMOS IFUs at the VLT to infer the spatial location of the ionized gas outflows. We compare the light distributions of the [O III] wing to that of the H$\beta$ broad emission line region, a classical point source (PSF). We then use the PSF to distinguish between the unresolved and resolved [O III] wing emission. We further determine its location using spectro-astrometry for the point-like sources. The [O III] wing is spatially unresolved in 23 out of the 36 AGN with >80 % of the flux associated with a point-like source. We measure <100 pc offsets in the spatial location of the outflow from the AGN nucleus using the spectro-astrometry technique for these sources. For the other 13 AGN, the [O III] wing emission is resolved and possibly extended on kpc scale. We conclude that [O III] wing emission can be compact or extended in an unbiased luminous AGN sample, where both cases are likely to appear. Electron density in the compact [O III] wing regions (median $n_e$~1900 cm$^{-3}$) is nearly a magnitude higher than in the extended ones (median $n_e$~500 cm$^{-3}$). The presence of spatially extended and compact [O III] wing emission is unrelated to the AGN bolometric luminosity and to inclination effects, which means other features such as time delays, or mechanical feedback/radio jets may shape the ionized gas outflow properties.


title: The Close AGN Reference Survey (CARS): No obvious signature of AGN feedback on star formation, but subtle trends

abstract: [Abridged] Active Galactic Nuclei (AGN) are thought to be responsible for the suppression of star formation in massive ~10$^{10}$ M$_\odot$ galaxies. While this process is a key feature in numerical simulations, it is not yet unambiguously confirmed in observational studies. Characterization of the star formation rate (SFR) in AGN host galaxies is challenging as AGN light contaminates most SFR tracers. We aim to obtain and compare SFR estimates from different tracers for AGN host galaxies in the Close AGN Reference Survey (CARS) to provide new observational insights. We construct integrated panchromatic spectral energy distributions (SED) to measure the FIR luminosity as a tracer for the recent (< 100 Myr) SFR. In addition, we use integral-field unit observation of the CARS targets to employ the H$\alpha$ luminosity decontaminated by AGN excitation as a proxy for the current (< 5 Myr) SFR. We find that significant differences in specific SFR of the AGN host galaxies as compared with the larger galaxy population disappear once cold gas mass, in addition to stellar mass, is used to predict the SFR. We identify individual galaxies with a significant difference in their SFR which can be related to a recent enhancement or decline in their SFR history that might be related to various processes including interactions, gas consumption, outflows and AGN feedback. AGN can occur in various stages of galaxy evolution which makes it difficult to relate the SFR solely to the impact of the AGN. We do not find any strong evidence for global positive or negative AGN feedback in the CARS sample. However, there is tentative evidence that 1) the relative orientation of the AGN engine with respect to the host galaxies might alter the efficiency of AGN feedback and 2) the recent SFH is an additional tool to identify rapid changes in galaxy growth driven by the AGN or other processes.


title: Fermionic Wave Functions from Neural-Network Constrained Hidden States

abstract: We introduce a systematically improvable family of variational wave functions for the simulation of strongly correlated fermionic systems. This family consists of Slater determinants in an augmented Hilbert space involving "hidden" additional fermionic degrees of freedom. These determinants are projected onto the physical Hilbert space through a constraint which is optimized, together with the single-particle orbitals, using a neural network parametrization. This construction draws inspiration from the success of hidden particle representations but overcomes the limitations associated with the mean-field treatment of the constraint often used in this context. Our construction provides an extremely expressive family of wave functions, which is proven to be universal. We apply this construction to the ground state properties of the Hubbard model on the square lattice, achieving levels of accuracy which are competitive with state-of-the-art variational methods.


title: Dispersion relation of a polaron in the Yang-Gaudin Bose gas

abstract: We study a one-dimensional Bose gas with two internal states described by the Yang-Gaudin model and calculate analytically the dispersion relation of a polaron quasiparticle, which is the lowest excitation branch. We study the dispersion in the thermodynamic limit in the regimes of weak and strong interaction without limitations on the momentum. At weak interaction the polaron dispersion is in the vicinity of the dark soliton one; we calculate the leading deviation in the parametric form. At strong interaction we find an ansatz for the explicit form of the polaron dispersion. It has the form of a power series of the sine function of the momentum with interaction-dependent coefficients. By increasing the power of the series, the corresponding coefficients show faster decay and thus one practically needs only a few of them; we give the results for the first three. The coefficients of the series are connected to the Maclaurin series of the polaron dispersion and thus it is sufficient to calculate latter quantity to automatically find the power series result for the dispersion at all momenta. The derived results for the dispersion can be used to obtain explicit expressions for the exponents of the power-law singularities in the response functions at the spectral edge.


title: Modeling Presymptomatic Spread in Epidemics via Mean-Field Games

abstract: This paper is concerned with developing mean-field game models for the evolution of epidemics. Specifically, an agent's decision -- to be socially active in the midst of an epidemic -- is modeled as a mean-field game with health-related costs and activity-related rewards. By considering the fully and partially observed versions of this problem, the role of information in guiding an agent's rational decision is highlighted. The main contributions of the paper are to derive the equations for the mean-field game in both fully and partially observed settings of the problem, to present a complete analysis of the fully observed case, and to present some analytical results for the partially observed case.


title: Notes on Confinement on $\mathbf{R^3 \times S^1}$: From Yang-Mills, super-Yang-Mills, and QCD(adj) to QCD(F)

abstract: This is a pedagogical introduction to the physics of confinement on $R^3 \times S^1$, using $SU(2)$ Yang-Mills with massive or massless adjoint fermions as the prime example; at the end, we also add fundamental flavours. The small-$S^1$ limit is remarkable, allowing for controlled semiclassical determination of the nonperturbative physics in these, mostly non-supersymmetric, theories. We begin by reviewing the Polyakov confinement mechanism on $R^3$. Moving on to $R^3 \times S^1$, we show how introducing adjoint fermions stabilizes center symmetry, leading to abelianization and semiclassical calculability. We explain how monopole-instantons and twisted monopole-instantons arise. We describe the role of various novel topological excitations in extending Polyakov's confinement to the locally four-dimensional case, discuss the nature of the confining string, and the $\theta$-angle dependence.~We study the global symmetry realization and, when available, present evidence for the absence of phase transitions as a function of the $S^1$ size. As our aim is not to cover all work on the subject, but to prepare the interested reader for its study, we also include brief descriptions of topics not covered in detail: the necessity for analytic continuation of path integrals, the study of more general theories, and the 't Hooft anomalies involving higher-form symmetries.


title: Shadowing, recurrence, and rigidity in dynamical systems

abstract: In this paper we examine the interplay between recurrence properties and the shadowing property in dynamical systems on compact metric spaces. In particular, we demonstrate that if the dynamical system $(X,f)$ has shadowing, then it is recurrent if and only if it is minimal. Furthermore, we show that a uniformly rigid system $(X,f)$ has shadowing if and only if $X$ is totally disconnected and use this to demonstrate the existence of a space $X$ for which no surjective system $(X,f)$ has shadowing. We further refine these results to discuss the dynamics that can occur in spaces with compact space of self-maps.


title: Flexible Inference of Optimal Individualized Treatment Strategy in Covariate Adjusted Randomization with Multiple Covariates

abstract: To maximize clinical benefit, clinicians routinely tailor treatment to the individual characteristics of each patient, where individualized treatment rules are needed and are of significant research interest to statisticians. In the covariate-adjusted randomization clinical trial with many covariates, we model the treatment effect with an unspecified function of a single index of the covariates and leave the baseline response completely arbitrary. We devise a class of estimators to consistently estimate the treatment effect function and its associated index while bypassing the estimation of the baseline response, which is subject to the curse of dimensionality. We further develop inference tools to identify predictive covariates and isolate effective treatment region. The usefulness of the methods is demonstrated in both simulations and a clinical data example.


title: Modeling and Analysis of the Landing Gear System with the Generalized Contracts

abstract: Nowadays, there are several complex systems in different sectors such as aviation, air traffic control ...etc. These systems do not have a precise perimeter, they are open and made of various specific components built with different languages and environments. The modeling, assembly and analysis of such open and complex heterogeneous systems are challenges in software engineering. This paper describes how the Minarets method decreases the difficulty of modeling, composition and analysis of the well known case study of the landing gear system. The method consists in: equipping individual components with generalized contracts that integrate various facets related to different concerns, composing these components according to their facets and verifying the resulting system with respect to the involved facets as well. The proposed method may be used or extended to cover more facets, and by strengthening assistance tool through proactive aspects in modeling, composing multi-facets contracts and finally the verification of the heterogeneous systems.


title: DIVeR: Real-time and Accurate Neural Radiance Fields with Deterministic Integration for Volume Rendering

abstract: DIVeR builds on the key ideas of NeRF and its variants -- density models and volume rendering -- to learn 3D object models that can be rendered realistically from small numbers of images. In contrast to all previous NeRF methods, DIVeR uses deterministic rather than stochastic estimates of the volume rendering integral. DIVeR's representation is a voxel based field of features. To compute the volume rendering integral, a ray is broken into intervals, one per voxel; components of the volume rendering integral are estimated from the features for each interval using an MLP, and the components are aggregated. As a result, DIVeR can render thin translucent structures that are missed by other integrators. Furthermore, DIVeR's representation has semantics that is relatively exposed compared to other such methods -- moving feature vectors around in the voxel space results in natural edits. Extensive qualitative and quantitative comparisons to current state-of-the-art methods show that DIVeR produces models that (1) render at or above state-of-the-art quality, (2) are very small without being baked, (3) render very fast without being baked, and (4) can be edited in natural ways.


title: Examining Two-Dimensional Luminosity-time Correlations for Gamma Ray Burst Radio Afterglows with VLA and ALMA

abstract: Gamma-ray burst (GRB) afterglow emission can be observed from sub-TeV to radio wavelengths, though only 6.6\% of observed GRBs present radio afterglows. We examine GRB radio light curves (LCs) to look for the presence of radio plateaus, resembling the plateaus observed in X-ray and optical. We analyze 404 GRBs from the literature with observed radio afterglow and fit 82 GRBs with at least 5 data points with a broken power law (BPL) model, requiring 4 parameters. From these, we find 18 GRBs that present a break feature resembling a plateau. We conduct the first multi-wavelength study of the Dainotti correlation between the luminosity $L_a$ and the rest-frame time of break $T_a^*$ for those 18 GRBs, concluding that the correlation exists and resembles the corresponding correlation in X-ray and optical wavelengths after correction for evolutionary effects. We compare the $T_a^*$ for the radio sample with $T_a^*$ values in X-ray and optical data (Dainotti et al. 2013, 2020b), finding significantly later break times in radio. We propose that this late break time and compatibility in slope suggests either a long-lasting plateau or the passage of a spectral break in the radio band. We also correct the distribution of the isotropic energy $E_{\rm iso}$ vs. the rest-frame burst duration $T^*_{90}$ for evolutionary effects and conclude that there is no significant difference between the $T^*_{90}$ distribution for the radio LCs with a break and those without.


title: Electro-osmosis at surfactant-laden liquid-gas interfaces: beyond standard models

abstract: Electro-osmosis (EO) is a powerful tool to manipulate liquids in micro and nanofluidic systems. While EO has been studied extensively at liquid-solid interfaces, the case of liquid-vapor interfaces, found e.g. in foam films and bubbles, remains to be explored. Here we perform molecular dynamics (MD) simulations of EO in a film of aqueous electrolyte covered with fluid layers of ionic surfactants and surrounded by gas. Following the experimental procedure, we compute the zeta potential from the EO velocity, defined as the velocity difference between the middle of the liquid film and the surrounding gas. We show that the zeta potential can be smaller or larger than existing predictions depending on the surfactant coverage. We explain the failure of previous descriptions by the fact that surfactants and bound ions move as rigid bodies and do not transmit the electric driving force to the liquid locally. Considering the reciprocal streaming current effect, we then develop an extended model, which can be used to predict the experimental zeta potential of surfactant-laden liquid-gas interfaces.


title: Some Error Analysis for the Quantum Phase Estimation Algorithms

abstract: This paper is concerned with the phase estimation algorithm in quantum computing algorithms, especially the scenarios where (1) the input vector is not an eigenvector; (2) the unitary operator is not exactly implemented; (3) random approximations are used for the unitary operator, e.g., the QDRIFT method. We characterize the probability of computing the phase values in terms of the consistency error, including the residual error, Trotter splitting error, or statistical mean-square error. In the first two cases, we show that in order to obtain the phase value with {error less or equal to $2^{-n}$ } and probability at least $1-\epsilon$, the required number of qubits is $ t \geq n + \log \big(2 + \frac{\delta^2 }{2 \epsilon \Delta\!E^2 } \big).$ The parameter $\delta$ quantifies the error associated with the inexact eigenvector and/or the unitary operator, and $\Delta\! E$ characterizes the spectral gap, i.e., the separation from the rest of the phase values. For the third case, we found a similar estimate, but the number of random steps has to be sufficiently large.


title: Measurement and QCD analysis of double-differential inclusive jet cross sections in proton-proton collisions at $\sqrt{s}$ = 13 TeV

abstract: A measurement of the inclusive jet production in proton-proton collisions at the LHC at $\sqrt{s}$ = 13 TeV is presented. The double-differential cross sections are measured as a function of the jet transverse momentum $p_\mathrm{T}$ and the absolute jet rapidity $\lvert y \rvert$. The anti-$k_\mathrm{T}$ clustering algorithm is used with distance parameter of 0.4 (0.7) in a phase space region with jet $p_\mathrm{T}$ from 97 GeV up to 3.1 TeV and $\lvert y \rvert$ $\lt$ 2.0. Data collected with the CMS detector are used, corresponding to an integrated luminosity of 36.3 fb$^{-1}$ (33.5 fb$^{-1}$). The measurement is used in a comprehensive QCD analysis at next-to-next-to-leading order, which results in significant improvement in the accuracy of the parton distributions in the proton. Simultaneously, the value of the strong coupling constant at the Z boson mass is extracted as $\alpha_\mathrm{S}$(Z) = 0.1170 $\pm$ 0.0019. For the first time, these data are used in a standard model effective field theory analysis at next-to-leading order, where parton distributions and the QCD parameters are extracted simultaneously with imposed constraints on the Wilson coefficient $c_1$ of 4-quark contact interactions. Note added: in the Addendum to this paper, available as Appendix B in this document, an improved value of $\alpha_\mathrm{S}$(Z) = 0.1166 $\pm$ 0.0017 has been extracted. This result supersedes the number in the above abstract of the original publication.


title: Probing Modified Gravity with Integrated Sachs-Wolfe CMB and Galaxy Cross-correlations

abstract: We use the cross-correlation power spectrum of the integrated Sachs-Wolfe (ISW) effect in the cosmic microwave background (CMB) temperature anisotropy and galaxy fluctuations to probe the physics of late-time cosmic acceleration. For this purpose, we focus on three models of dark energy that belong to a sub-class of Horndeski theories with the speed of gravity equivalent to that of light: Galileon Ghost Condensate (GGC), Generalized Cubic Covariant Galileon (GCCG), and K-mouflage. In the GGC and GCCG models, the existence of cubic-order scalar self-interactions allows a possibility for realizing negative ISW-galaxy cross-correlations, while the K-mouflage model predicts a positive correlation similar to the $\Lambda$-cold-dark-matter (LCDM) model. In our analysis, we fix the parameters of each model to their best-fit values derived from a baseline likelihood analysis with observational data from CMB, baryon acoustic oscillations, and supernovae type Ia. Then we fit those best-fit models to the ISW-galaxy cross-correlation power spectrum extracted from a collection of photometric redshift surveys. We find that both GGC and GCCG best-fit models degrade the fit to the ISW-galaxy cross-correlation data compared to LCDM best-fit model. This is attributed to the fact that, for their best-fit values constrained from the baseline likelihood, the cubic-order scalar self-interaction gives rise to suppressed ISW tails relative to $\Lambda$CDM. The K-mouflage best-fit model is largely degenerate with the LCDM best-fit model and has a positively correlated ISW-galaxy power close to that of LCDM.


title: Qubit Decoherence and Symmetry Restoration through Real-Time Instantons

abstract: A parametrically driven quantum oscillator, stabilized by a nonlinear dissipation, exhibits a spontaneous breaking of the parity symmetry. It results in the quantum bi-stability, corresponding to a Bloch sphere of dark states. This makes such a driven-dissipative system an attractive candidate for a qubit. The parity symmetry breaking is exact both on the classical level and within the quantum mechanical perturbation theory. Here we show that non-perturbative quantum effects lead to the symmetry restoration and result in exponentially small but finite qubit decoherence rate. Technically the symmetry restoration is due to real time instanton trajectories of the Keldysh path integral, which represents the Lindbladian evolution of the driven-dissipative oscillator.


title: Machine Learning for Mechanical Ventilation Control (Extended Abstract)

abstract: Mechanical ventilation is one of the most widely used therapies in the ICU. However, despite broad application from anaesthesia to COVID-related life support, many injurious challenges remain. We frame these as a control problem: ventilators must let air in and out of the patient's lungs according to a prescribed trajectory of airway pressure. Industry-standard controllers, based on the PID method, are neither optimal nor robust. Our data-driven approach learns to control an invasive ventilator by training on a simulator itself trained on data collected from the ventilator. This method outperforms popular reinforcement learning algorithms and even controls the physical ventilator more accurately and robustly than PID. These results underscore how effective data-driven methodologies can be for invasive ventilation and suggest that more general forms of ventilation (e.g., non-invasive, adaptive) may also be amenable.


title: Accounting for Selection Bias and Redshift Evolution in GRB Radio Afterglow Data

abstract: Gamma-ray Bursts (GRBs) are highly energetic events that can be observed at extremely high redshift. However, inherent bias in GRB data due to selection effects and redshift evolution can significantly skew any subsequent analysis. We correct for important variables related to the GRB emission, such as the burst duration, $T_{90}^*$, the prompt isotropic energy, $E_{\rm iso}$, the rest-frame end time of the plateau emission, $T_{\rm a, radio}^*$, and its correspondent luminosity $L_{\rm a, radio}$, for radio afterglow. In particular, we use the Efron-Petrosian method presented in 1992 for the correction of our variables of interest. Specifically, we correct $E_{\rm iso}$ and $T_{90}^*$ for 80 GRBs, and $L_{\rm a, radio}$ and $T_{\rm a, radio}^*$ for a subsample of 18 GRBs that present a plateau-like flattening in their light curve. Upon application of this method, we find strong evolution with redshift in most variables, particularly in $L_{\rm a, radio}$, with values similar to those found in past and current literature in radio, X-ray and optical wavelengths, indicating that these variables are susceptible to observational bias. This analysis emphasizes the necessity of correcting observational data for evolutionary effects to obtain the intrinsic behavior of correlations to use them as discriminators among the most plausible theoretical models and as reliable cosmological tools.


title: A counter-example to the probabilistic universal graph conjecture via randomized communication complexity

abstract: We refute the Probabilistic Universal Graph Conjecture of Harms, Wild, and Zamaraev, which states that a hereditary graph property admits a constant-size probabilistic universal graph if and only if it is stable and has at most factorial speed. Our counter-example follows from the existence of a sequence of $n \times n$ Boolean matrices $M_n$, such that their public-coin randomized communication complexity tends to infinity, while the randomized communication complexity of every $\sqrt{n}\times \sqrt{n}$ submatrix of $M_n$ is bounded by a universal constant.


title: Modeling the diffusion-erosion crossover dynamics in drug release

abstract: A computational model is proposed to investigate drug delivery systems in which erosion and diffusion mechanisms are participating in the drug release process. Our approach allowed us to analytically estimate the crossover point between those mechanisms through the value of the parameter $b$ ($b_c = 1$) and the scaling behavior of parameter $\tau$ on the Weibull function, $\exp[-(t/\tau)^b]$, used to adjust drug release data in pharmaceutical literature. Numerical investigations on the size dependence of the characteristic release time $\tau$ found it to satisfy either linear or quadratic scaling relations on either erosive or diffusive regimes. Along the crossover the characteristic time scales with the average coefficient observed on the extreme regimes ({\it i.e.}, $\tau \sim L^{3/2}$), and we show that this result can be derived analytically by assuming an Arrhenius relation for the diffusion coefficient inside the capsule. Based on these relations a phenomenological expression for the characteristic release in terms of size $L$ and erosion rate $\kappa$ is proposed, which can be useful for predicting the crossover erosion rate $\kappa_c$. We applied this relation to the experimental literature data for the release of acetaminophen immersed in a wax matrix and found them to be consistent with our numerical results.


title: Multi-Sensory HMI for Human-Centric Industrial Digital Twins: A 6G Vision of Future Industry

abstract: The next revolution of industry will turn the industries as well as the entire society into a human-centric shape. The human presence in industrial environment and the human participation in industrial processes will be magnified more than ever before. To cope with the emerging challenges raised by this revolution, 6G ambitions to bridge the three domains of digital information, physical assets and humans into one merged cyber-physical-human world. This proposes not only an unprecedented demand for digital twin solutions, but also new technical requirements. Especially, aiming at a human-centric industrial DT system, novel multi-sensory human-machine interfaces will play a key role in this paradigm shift.


title: Modular-topology optimization of structures and mechanisms with free material design and clustering

abstract: Topology optimization of modular structures and mechanisms enables balancing the performance of automatically-generated individualized designs, as required by Industry 4.0, with enhanced sustainability by means of component reuse. For optimal modular design, two key questions must be answered: (i) what should the topology of individual modules be like and (ii) how should modules be arranged at the product scale? We address these challenges by proposing a bi-level sequential strategy that combines free material design, clustering techniques, and topology optimization. First, using free material optimization enhanced with post-processing for checkerboard suppression, we determine the distribution of elasticity tensors at the product scale. To extract the sought-after modular arrangement, we partition the obtained elasticity tensors with a novel deterministic clustering algorithm and interpret its outputs within Wang tiling formalism. Finally, we design interiors of individual modules by solving a single-scale topology optimization problem with the design space reduced by modular mapping, conveniently starting from an initial guess provided by free material optimization. We illustrate these developments with three benchmarks first, covering compliance minimization of modular structures, and, for the first time, the design of non-periodic compliant modular mechanisms. Furthermore, we design a set of modules reusable in an inverter and in gripper mechanisms, which ultimately pave the way towards the rational design of modular architectured (meta)materials.


title: Thermopower and thermophase in a $d$-wave superconductor

abstract: In an unconventional superconductor, the interplay of scattering off impurities and Andreev processes may lead to different scattering times for electronlike and holelike quasiparticles. Such electron-hole asymmetry appears when the impurity scattering phase shift is intermediate between the Born and unitary limits and leads to an expectation for large thermoelectric effects. Here, we examine the thermoelectric response of a $d$-wave superconductor connected to normal-metal reservoirs under a temperature bias using a fully self-consistent quasiclassical theory. The thermoelectrically induced quasiparticle current is cancelled by superflow in an open circuit setup, but at the cost of a charge imbalance induced at the contacts and extending across the structure. We investigate the resulting thermopower and thermophase and their dependencies on scattering phase shift, mean free path, and interface transparency. For crystal-axis orientations such that surface-bound zero-energy Andreev states are formed, the thermoelectric effect is reduced as a result of locally reduced electron-hole asymmetry. For a semiballistic superconductor with good contacts, we find thermopowers of order several ${\mu}V/K$, suggesting a thermovoltage measurement as a promising path to investigate thermoelectricity in unconventional superconductors.


title: Lightcurve Evolution of the nearest Tidal Disruption Event: A late-time, radio-only flare

abstract: Tidal disruption events (TDEs) occur when a star passes close enough to a galaxy's supermassive black hole to be disrupted by tidal forces. We discuss new observations of IGRJ12580+0134, a TDE observed in NGC 4845 (d=17 Mpc) in November 2010, with the Karl G. Jansky Very Large Array (JVLA). We also discuss a reanalysis of 2010-2011 Swift and XMM-Newton observations, as well as new, late-time Swift observations. Our JVLA observations show a decay of the nuclear radio flux until 2015, when a plateau was seen, and then a significant (~factor 3) radio flare during 2016. The 2016 radio flare was also accompanied by radio spectral changes, but was not seen in the X-rays. We model the flare as resulting from the interaction of the nuclear jet with a cloud in the interstellar medium. This is distinct from late-time X-ray flares in a few other TDEs where changes in the accretion state and/or a fallback event were suggested, neither of which appears possible in this case. Our reanalysis of the Swift and XMM-Newton data from 2011 shows significant evidence for thermal emission from a disk, as well as a very soft power-law. This, in addition to the extreme X-ray flux increase seen in 2010 (a factor of >$100) bolsters the identification of IGRJ12580+0134 as a TDE, not an unusual AGN variability event.


title: Thermodynamics of Equilibrium Alkali Plasma. Simple and Accurate Analytical Model for Non-Trivial Case

abstract: Analytical equation of state for pure alkali metals (lithium, sodium, potassium, rubidium and cesium) in equilibrium gas phase is proposed. This equation has a simple form, generalizes the equation of state for a perfect gas and is universal for all alkali elements. It correctly reproduces the experimental data for the equilibrium gas phase over a wide range of pressures (up to $\sim 10^5$~Pa) and temperatures (up to $\sim 3 \cdot 10^3$~K). On the basis of this equation of state, expressions for thermal and caloric coefficients as well as for \textit{other physical characteristics} are obtained. The results of this study confirm feasibility of the principle of corresponding states in relation to the group of alkali elements.

categories: cond-mat.mtrl-sci cond-mat.dis-nn cond-mat.stat-mech physics.comp-ph physics.plasm-ph


title: Evaluation of automated airway morphological quantification for assessing fibrosing lung disease

abstract: Abnormal airway dilatation, termed traction bronchiectasis, is a typical feature of idiopathic pulmonary fibrosis (IPF). Volumetric computed tomography (CT) imaging captures the loss of normal airway tapering in IPF. We postulated that automated quantification of airway abnormalities could provide estimates of IPF disease extent and severity. We propose AirQuant, an automated computational pipeline that systematically parcellates the airway tree into its lobes and generational branches from a deep learning based airway segmentation, deriving airway structural measures from chest CT. Importantly, AirQuant prevents the occurrence of spurious airway branches by thick wave propagation and removes loops in the airway-tree by graph search, overcoming limitations of existing airway skeletonisation algorithms. Tapering between airway segments (intertapering) and airway tortuosity computed by AirQuant were compared between 14 healthy participants and 14 IPF patients. Airway intertapering was significantly reduced in IPF patients, and airway tortuosity was significantly increased when compared to healthy controls. Differences were most marked in the lower lobes, conforming to the typical distribution of IPF-related damage. AirQuant is an open-source pipeline that avoids limitations of existing airway quantification algorithms and has clinical interpretability. Automated airway measurements may have potential as novel imaging biomarkers of IPF severity and disease extent.


title: Infrared Absolute Calibration I: Comparison of Sirius with Fainter Calibration Stars

abstract: A challenge in absolute calibration is to relate very bright stars with physical flux measurements to faint ones within range of modern instruments, e.g. those on large groundbased telescopes or on the James Webb Space Telescope (JWST). We propose Sirius as the fiducial color standard: it is an A0V star that is slowly rotating and does not have infrared excesses either due to hot dust or a planetary debris disk; it also has a number of accurate (~ 1 - 2%) absolute flux measurements. We transfer the near infrared flux from Sirius accurately to BD +60 1753, an unobscured early A-type star (A1V, V ~ 9.6, E(B-V) ~ 0.009) that is faint enough to serve as a primary absolute flux calibrator for JWST. Its near-infrared spectral energy distribution and that of Sirius should be virtually identical. We have determined its output relative to that of Sirius in a number of different ways, all of which give consistent results within ~ 1%. We also transfer the calibration to GSPC P330-E, a well calibrated close solar analog (G2V). We have emphasized the 2MASS K_S band since it represents a large number and long history of measurements, but theoretical spectra (i.e., from CALSPEC) of these stars can be used to extend this result throughout the near- and mid-infrared.


title: Multifractal Magnetoconductance Fluctuations in Mesoscopic Systems

abstract: We perform a multifractal detrended fluctuation analysis of the magnetoconductance data of two standard types of mesoscopic systems: a disordered nanowire and a ballistic chaotic billiard, with two different lattice structures. We observe in all cases that multifractality is generally present and that it becomes stronger in the quantum regime of conduction, i.e., when the number of open scattering channels is small. We argue that this behavior originates from correlations induced by the magnetic field, which can be characterized through the distribution of conductance increments in the corresponding "stochastic time series," with the magnetic field playing the role of a fictitious time. More specifically, we show that the distributions of conductance increments are well fitted by $q$ Gaussians and that the value of the parameter $q$ is a useful quantitative measure of multifractality in magnetoconductance fluctuations.


title: Multiplicity structure of the arc space of a fat point

abstract: The equation $x^m = 0$ defines a fat point on a line. The algebra of regular functions on the arc space of this scheme is the quotient of $k[x, x', x^{(2)}, \ldots]$ by all differential consequences of $x^m = 0$. This infinite-dimensional algebra admits a natural filtration by finite dimensional algebras corresponding to the truncations of arcs. We show that the generating series for their dimensions equals $\frac{m}{1 - mt}$. We also determine the lexicographic initial ideal of the defining ideal of the arc space. These results are motivated by nonreduced version of the geometric motivic Poincar\'e series, multiplicities in differential algebra, and connections between arc spaces and the Rogers-Ramanujan identities. We also prove a recent conjecture put forth by Afsharijoo in the latter context.


title: DyFormer: A Scalable Dynamic Graph Transformer with Provable Benefits on Generalization Ability

abstract: Transformers have achieved great success in several domains, including Natural Language Processing and Computer Vision. However, its application to real-world graphs is less explored, mainly due to its high computation cost and its poor generalizability caused by the lack of enough training data in the graph domain. To fill in this gap, we propose a scalable Transformer-like dynamic graph learning method named Dynamic Graph Transformer (DyFormer) with spatial-temporal encoding to effectively learn graph topology and capture implicit links. To achieve efficient and scalable training, we propose temporal-union graph structure and its associated subgraph-based node sampling strategy. To improve the generalization ability, we introduce two complementary self-supervised pre-training tasks and show that jointly optimizing the two pre-training tasks results in a smaller Bayesian error rate via an information-theoretic analysis. Extensive experiments on the real-world datasets illustrate that DyFormer achieves a consistent 1%-3% AUC gain (averaged over all time steps) compared with baselines on all benchmarks.


title: Parallel algorithms for computing the tensor-train decomposition

abstract: The tensor-train (TT) decomposition expresses a tensor in a data-sparse format used in molecular simulations, high-order correlation functions, and optimization. In this paper, we propose four parallelizable algorithms that compute the TT format from various tensor inputs: (1) Parallel-TTSVD for traditional format, (2) PSTT and its variants for streaming data, (3) Tucker2TT for Tucker format, and (4) TT-fADI for solutions of Sylvester tensor equations. We provide theoretical guarantees of accuracy, parallelization methods, scaling analysis, and numerical results. For example, for a $d$-dimension tensor in $\mathbb{R}^{n\times\dots\times n}$, a two-sided sketching algorithm PSTT2 is shown to have a memory complexity of $\mathcal{O}(n^{\lfloor d/2 \rfloor})$, improving upon $\mathcal{O}(n^{d-1})$ from previous algorithms.


title: Neutrino spin operator and dispersion in moving matter

abstract: We found a spin integral of motion for neutrino, which propagates in moving and polarized matter. The operator obtained opens up the possibility of consistent classification of neutrino states in such a medium and, as a consequence, a systematic description of the related physical phenomena. Using the operator, we obtain a dispersion relation for neutrinos and consider its particular cases.


title: Birth-death chains on a spider: spectral analysis and reflecting-absorbing factorization

abstract: We consider discrete-time birth-death chains on a spider, i.e. a graph consisting of $N$ discrete half lines on the plane that are joined at the origin. This process can be identified with a discrete-time quasi-birth-death process on the state space $\mathbb{N}_0 \times \{1, 2, \dots, N\},$ represented by a block tridiagonal transition probability matrix. We prove that we can analyze this process by using spectral methods and obtain the $n$-step transition probabilities in terms of a weight matrix and the corresponding matrix-valued orthogonal polynomials (the so-called Karlin-McGregor formula). We also study under what conditions we can get a reflecting-absorbing factorization of the birth-death chain on a spider which can be seen as a stochastic UL block factorization of the transition probability matrix of the quasi-birth-death process. With this factorization we can perform a discrete Darboux transformation and get new families of "almost" birth-death chains on a spider. The spectral matrix associated with the Darboux transformation will be a Geronimus transformation of the original spectral matrix. Finally, we apply our results to the random walk on a spider, i.e. with constant transition probabilities.


title: A Swift study of long-term changes in the X-ray flaring properties of Sagittarius A*

abstract: The radiative counterpart of the supermassive black hole at the Galactic Centre, Sagittarius A*, displays flaring emission in the X-ray band atop a steady, quiescent level. Flares are also observed in the near-infrared band. The physical process producing the flares is not fully understood and it is unclear if the flaring rate varies, although some recent works suggest it has reached unprecedented variability in recent years. Using over a decade of regular X-ray monitoring of Neil Gehrels Swift Observatory, we studied the variations in count rate of Sgr A* on time scales of years. We decomposed the X-ray emission into quiescent and flaring emission, modelled as a constant and power law process, respectively. We found that the complete, multi-year dataset cannot be described by a stationary distribution of flare fluxes, while individual years follow this model better. In three of the ten studied years, the data is consistent with a purely Poissonian quiescent distribution, while for five years only an upper limit of the flare flux distribution parameter could be determined. We find that these possible changes cannot be explained fully by the different number of observations per year. Combined, these results are instead consistent with a changing flaring rate of Sgr A*, appearing more active between 2006-2007 and 2017-2019, than between 2008-2012. Finally, we discuss this result in the context of flare models and the passing of gaseous objects, and discuss the extra statistical steps taken, for instance to deal with the background in the Swift observations.


title: MURAL: An Unsupervised Random Forest-Based Embedding for Electronic Health Record Data

abstract: A major challenge in embedding or visualizing clinical patient data is the heterogeneity of variable types including continuous lab values, categorical diagnostic codes, as well as missing or incomplete data. In particular, in EHR data, some variables are {\em missing not at random (MNAR)} but deliberately not collected and thus are a source of information. For example, lab tests may be deemed necessary for some patients on the basis of suspected diagnosis, but not for others. Here we present the MURAL forest -- an unsupervised random forest for representing data with disparate variable types (e.g., categorical, continuous, MNAR). MURAL forests consist of a set of decision trees where node-splitting variables are chosen at random, such that the marginal entropy of all other variables is minimized by the split. This allows us to also split on MNAR variables and discrete variables in a way that is consistent with the continuous variables. The end goal is to learn the MURAL embedding of patients using average tree distances between those patients. These distances can be fed to nonlinear dimensionality reduction method like PHATE to derive visualizable embeddings. While such methods are ubiquitous in continuous-valued datasets (like single cell RNA-sequencing) they have not been used extensively in mixed variable data. We showcase the use of our method on one artificial and two clinical datasets. We show that using our approach, we can visualize and classify data more accurately than competing approaches. Finally, we show that MURAL can also be used to compare cohorts of patients via the recently proposed tree-sliced Wasserstein distances.


title: Enhancing the sensitivity and selectivity of pyrene-based sensors for detection of small gaseous molecules via destructive quantum interference

abstract: Graphene-based sensors are exceptionally sensitive with high carrier mobility and low intrinsic noise, and have been intensively investigated in the past decade. The detection of individual gas molecules has been reported, albeit the underlying sensing mechanism is not yet well understood. We focus on the adsorption of NO$_2$, H$_2$O, and NH$_3$ on a molecular junction with a pyrene core, which can be considered as a minimal graphene-like unit. We systematically investigate the chemiresistive response within the framework of density functional theory and non-equilibrium Green's functions. We highlight the fundamental role of quantum interference (QI) in the sensing process, and we propose it as a paradigmatic mechanism for sensing. Owing to the open-shell character of NO$_2$, its interaction with pyrene gives rise to a Fano resonance thereby triggering the strongest chemiresistive response, while the weaker interactions with H$_2$O and NH$_3$ result in lower sensitivity. We demonstrate that by exploiting destructive QI arising in the meta-substituted pyrene, it is possible to calibrate the sensor to enhance both its sensitivity and chemical selectivity by almost two orders of magnitude so that individual molecules can be detected and distinguished. These results provide a fundamental strategy to design high-performance chemical sensors with graphene functional blocks.


title: Spherical harmonic shape descriptors of nodal force demands for quantifying spatial truss connection complexity

abstract: The connections of a spatial truss structure play a critical role in the safe and efficient transfer of axial forces between members. For discrete connections, they can also improve construction efficiency by acting as registration devices that lock members in precise orientations. As more geometrically complex spatial trusses are enabled by computational workflows and the demand for material-efficient spanning systems, there is a need to understand the effects of global form on the demands at the connections. For large-scale structures with irregular geometry, customizing individual nodes to meet exact member orientations and force demands may be infeasible; conversely, standardizing all connections results in oversized nodes and a compromise in registration potential. We propose a method for quantifying the complexity of spatial truss designs by the variation in nodal force demands. By representing nodal forces as a geometric object, we leverage the spherical harmonic shape descriptor, developed for applications in computational geometry, to characterize each node by a rotation and translation-invariant fixed-length vector. We define a complexity score for spatial truss design by the variance in the positions of the feature vectors in higher-dimensional space, providing an additional performance metric during early stage design exploration. We then develop a pathway towards reducing complexity by clustering nodes with respect to their feature vectors to reduce the number of unique connectors for design while minimizing the effects of mass standardization.


title: Evidence for subdominant multipole moments and precession in merging black-hole-binaries from GWTC-2.1

abstract: The LIGO--Virgo--KAGRA collaborations (LVK) produced a catalogue containing gravitational-wave (GW) observations from the first half of the third GW observing run (O3a). This catalogue, GWTC-2.1, includes for the first time a number of \emph{exceptional} GW candidates produced from merging black-hole-binaries with unequivocally unequal component masses. Since subdominant multipole moments and spin-induced orbital precession are more likely to leave measurable imprints on the emitted GW from unequal component mass binaries, these general relativistic phenomena may now be measurable. Indeed, both GW190412 and GW190814 have already shown conclusive evidence for subdominant multipole moments. This provides valuable insights into the dynamics of the binary. We calculate the evidence for subdominant multipole moments and spin-induced orbital precession for all merging black-hole-binaries in GWTC-2.1 that were observed during O3a and show that (a) no gravitational-wave candidate has measurable higher order multipole content beyond $\ell = 3$, (b) in addition to the confident subdominant multipole measurements in GW190412 and GW190814, GW190519\_153544 and GW190929\_012149 show marginal evidence for the $(\ell, |m|) = (3, 3)$ subdominant multipole, (c) GW190521 may have measurable subdominant multipole content and (d) GW190412 may show evidence for spin-induced orbital precession.


title: A lattice structure for ancestral configurations arising from the relationship between gene trees and species trees

abstract: To a given gene tree topology $G$ and species tree topology $S$ with leaves labeled bijectively from a fixed set $X$, one can associate a set of ancestral configurations, each of which encodes a set of gene lineages that can be found at a given node of the species tree. We introduce a lattice structure on ancestral configurations, studying the directed graphs that provide graphical representations of lattices of ancestral configurations. For a matching gene tree topology and species tree topology $G=S$, we present a method for defining the digraph of ancestral configurations from the tree topology by using iterated cartesian products of graphs. We show that a specific set of paths on the digraph of ancestral configurations is in bijection with the set of labeled histories -- a well-known phylogenetic object that enumerates possible temporal orderings of the coalescences of a tree. For each of a series of tree families, we obtain closed-form expressions for the number of labeled histories by using this bijection to count paths on associated digraphs. Finally, we prove that our lattice construction extends to nonmatching tree pairs, and we use it to characterize pairs $(G,S)$ having the maximal number of ancestral configurations for a fixed $G$. We discuss how the construction provides new methods for performing enumerations of combinatorial aspects of gene and species trees.


title: Model selection applied to reconstructions of the Dark Energy

abstract: The main aim of this paper is to perform a model comparison for some reconstructions of the key properties that describe the dark energy of the Universe i.e. energy density and the equation of state (EoS). We carry out this process by using a binning and a linear interpolation methodologies, and on top of that, we incorporate a correlation function mechanism. An extension of the two of them was also introduced, where internal amplitudes are allowed to vary in height as well as in position. The reconstructions were made with data from the Hubble parameter, Supernovae Type Ia and Baryon Acoustic Oscillations (H+SN+BAO), all of which span a range from $z=0.01$ to $z=2.34$. First we perform the parameter estimation for each of the reconstructions to then provide a model selection through the Bayesian Evidence. Throughout our process we found a better fit to the data, up to $4\sigma$ compared to $\Lambda$CDM, and the presence of some interesting features, i.e. an oscillatory behaviour at late times, a decrease in the dark energy density component at early times and a transition to the phantom divide-line in the EoS. To discern these features from noisy contributions, we include a principal component analysis and found that some of these characteristics should be taken into account to satisfy current observations.


title: INCHE: High-Performance Encoding for Relational Databases through Incrementally Homomorphic Encryption

abstract: Homomorphic encryption (HE) offers data confidentiality by executing queries directly on encrypted fields in the database-as-a-service (DaaS) paradigm. While fully HE exhibits great expressiveness but prohibitive performance overhead, a better balance between flexibility and efficiency can be achieved by partially HE schemes. Performance-wise, however, the encryption rate of state-of-the-art HE schemes is still orders of magnitude lower than the I/O throughput, rendering the HE scheme the performance bottleneck. This paper proposes INCHE, an incrementally homomorphic encryption scheme, which aims to boost the performance of HE schemes by incrementally encrypting fields in relational databases. The key idea of INCHE is to explore the intrinsic correlation between plaintexts and cache them for future reuse such that expensive HE primitives from plaintexts to ciphertexts are avoided. We prove the semantic security of INCHE under the chosen-plaintext attack (CPA) model and show that its time complexity is linear in the plaintext length. We implement an INCHE prototype by extending the Symmetria cryptosystem and verify its effectiveness on both randomly-generated data and the TPC-H benchmark.


title: Identifying Population Movements with Non-Negative Matrix Factorization from Wi-Fi User Counts in Smart and Connected Cities

abstract: Non-Negative Matrix Factorization (NMF) is a valuable matrix factorization technique which produces a "parts-based" decomposition of data sets. Wi-Fi user counts are a privacy-preserving indicator of population movements in smart and connected urban environments. In this paper, we apply NMF with a novel matrix embedding to Wi-Fi user count data from the University of Colorado at Boulder Campus for the purpose of automatically identifying patterns of human movement in a Smart and Connected infrastructure environment.


title: Compactness of Fixed Point Maps and the Ball-Marsden-Slemrod Conjecture

abstract: Given a parameter dependent fixed point equation $x = F(x,u)$, we derive an abstract compactness principle for the fixed point map $u \mapsto x^*(u)$ under the assumptions that (i) the fixed point equation can be solved by the contraction principle and (ii) the map $u \mapsto F(x,u)$ is compact for fixed $x$. This result is applied to infinite-dimensional, semi-linear control systems and their reachable sets. More precisely, we extend a non-controllability result of Ball, Marsden, and Slemrod [1] to semi-linear systems. First we consider $L^p$-controls, $p>1$. Subsequently we analyze the case $p=1$.


title: Gaussian Process Inference Using Mini-batch Stochastic Gradient Descent: Convergence Guarantees and Empirical Benefits

abstract: Stochastic gradient descent (SGD) and its variants have established themselves as the go-to algorithms for large-scale machine learning problems with independent samples due to their generalization performance and intrinsic computational advantage. However, the fact that the stochastic gradient is a biased estimator of the full gradient with correlated samples has led to the lack of theoretical understanding of how SGD behaves under correlated settings and hindered its use in such cases. In this paper, we focus on hyperparameter estimation for the Gaussian process (GP) and take a step forward towards breaking the barrier by proving minibatch SGD converges to a critical point of the full log-likelihood loss function, and recovers model hyperparameters with rate $O(\frac{1}{K})$ for $K$ iterations, up to a statistical error term depending on the minibatch size. Our theoretical guarantees hold provided that the kernel functions exhibit exponential or polynomial eigendecay which is satisfied by a wide range of kernels commonly used in GPs. Numerical studies on both simulated and real datasets demonstrate that minibatch SGD has better generalization over state-of-the-art GP methods while reducing the computational burden and opening a new, previously unexplored, data size regime for GPs.


title: Online Coverage Planning for an Autonomous Weed Mowing Robot with Curvature Constraints

abstract: The land used for grazing cattle takes up about one-third of the land in the United States. These areas can be highly rugged. Yet, they need to be maintained to prevent weeds from taking over the nutritious grassland. This can be a daunting task especially in the case of organic farming since herbicides cannot be used. In this paper, we present the design of Cowbot, an autonomous weed mowing robot for pastures. Cowbot is an electric mower designed to operate in the rugged environments on cow pastures and provide a cost-effective method for weed control in organic farms. Path planning for the Cowbot is challenging since weed distribution on pastures is unknown. Given a limited field of view, online path planning is necessary to detect weeds and plan paths to mow them. We study the general online path planning problem for an autonomous mower with curvature and field of view constraints. We develop two online path planning algorithms that are able to utilize new information about weeds to optimize path length and ensure coverage. We deploy our algorithms on the Cowbot and perform field experiments to validate the suitability of our methods for real-time path planning. We also perform extensive simulation experiments which show that our algorithms result in up to 60 % reduction in path length as compared to baseline boustrophedon and random-search based coverage paths.


title: Correlations between Panoramic Imagery and Gamma-Ray Background in an Urban Area

abstract: When searching for radiological sources in an urban area, a vehicle-borne detector system will often measure complex, varying backgrounds primarily from natural gamma-ray sources. Much work has been focused on developing spectral algorithms that retain sensitivity and minimize the false positive rate even in the presence of such spectral and temporal variability. However, information about the environment surrounding the detector system might also provide useful clues about the expected background, which if incorporated into an algorithm, could improve performance. Recent work has focused on extensive measuring and modeling of urban areas with the goal of understanding how these complex backgrounds arise. This work presents an analysis of panoramic video images and gamma-ray background data collected in Oakland, California by the Radiological Multi-sensor Analysis Platform (RadMAP) vehicle. Features were extracted from the panoramic images by semantically labeling the images and then convolving the labeled regions with the detector response. A linear model was used to relate the image-derived features to gamma-ray spectral features obtained using Non-negative Matrix Factorization (NMF) under different regularizations. We find some gamma-ray background features correlate strongly with image-derived features that measure the response-adjusted solid angle subtended by sky and buildings, and we discuss the implications for the development of future, contextually-aware detection algorithms.


title: Likelihood analysis of the flavour anomalies and $g-2$ in the general two Higgs doublet model

abstract: We present a likelihood analysis of the general two Higgs doublet model, using the most important currently measured flavour observables, in view of the anomalies in charged current tree-level and neutral current one-loop rare decays of $B$ mesons in $b\to c l \overline{\nu}$ and $b\to s\mu^{+}\mu^{-}$ transitions, respectively. We corroborate that the model explains the latter and it is able to simultaneously fit the experimental values of the $R(D)$ ratio at $1\sigma$, but it can not accommodate the $R(D^{*})$ and $F_{L}(D^{*})$ observables. We find that the fitted values for the angular observables in $b\to s\mu^{+}\mu^{-}$ transitions exhibit better agreement with the general two Higgs double model in comparison to the SM. We also make predictions for future collider observables $\mathrm{BR}(t\to ch)$, $\mathrm{BR}(h\to bs)$, $\mathrm{BR}(h\to \tau\mu)$, $\mathrm{BR}(B_{s}\rightarrow\tau^{+}\tau^{-})$, $\mathrm{BR}(B^{+}\rightarrow K^{+}\tau^{+}\tau^{-})$ and the flavour violating decays of the $\tau$ lepton, $\mathrm{BR}(\tau\rightarrow3\mu)$ and $\mathrm{BR}(\tau\to\mu\gamma)$. The model predicts values of $\mathrm{BR}(t\to ch)$, $\mathrm{BR}(B_{s}\rightarrow\tau^{+}\tau^{-})$ and $\mathrm{BR}(B^{+}\rightarrow K^{+}\tau^{+}\tau^{-})$ that are out of reach of future experiments, but its predictions for $\mathrm{BR}(h\to bs)$ and $\mathrm{BR}(h\to \tau\mu)$ are within the future sensitivity of the HL-LHC or the ILC. We also find that the predictions for the $\tau\rightarrow3\mu$ and $\tau\to\mu\gamma$ decays are well within the projected limits of the Belle II experiment. Finally, using the latest measurement of the Fermilab Muon $g-2$ Collaboration, we performed a simultaneous fit to $\Delta a_{\mu}$ constrained by the charged anomalies, finding solutions at the $1\sigma$ level. Once the neutral anomalies are included, however, a simultaneous explanation is unfeasible.