title: The application of sub-seasonal to seasonal (S2S) predictions for hydropower forecasting

abstract: Inflow forecasts play an essential role in the management of hydropower reservoirs. Forecasts help operators schedule power generation in advance to maximise economic value, mitigate downstream flood risk, and meet environmental requirements. The horizon of operational inflow forecasts is often limited in range to ~2 weeks ahead, marking the predictability barrier of deterministic weather forecasts. Reliable inflow forecasts in the sub-seasonal to seasonal (S2S) range would allow operators to take proactive action to mitigate risks of adverse weather conditions, thereby improving water management and increasing revenue. This study outlines a method of deriving skilful S2S inflow forecasts using a case study reservoir in the Scottish Highlands. We generate ensemble inflow forecasts by training a linear regression model for the observed inflow onto S2S ensemble precipitation predictions from the European Centre for Medium-range Weather Forecasting (ECMWF). Subsequently, post-processing techniques from Ensemble Model Output Statistics are applied to derive calibrated S2S probabilistic inflow forecasts, without the application of a separate hydrological model. We find the S2S probabilistic inflow forecasts hold skill relative to climatological forecasts up to 6 weeks ahead. The inflow forecasts hold greater skill during winter compared with summer. The forecasts, however, struggle to predict high summer inflows, even at short lead-times. The potential for the S2S probabilistic inflow forecasts to improve water management and deliver increased economic value is confirmed using a stylised cost model. While applied to hydropower forecasting, the results and methods presented here are relevant to broader fields of water management and S2S forecasting applications.


title: Enhancing audio quality for expressive Neural Text-to-Speech

abstract: Artificial speech synthesis has made a great leap in terms of naturalness as recent Text-to-Speech (TTS) systems are capable of producing speech with similar quality to human recordings. However, not all speaking styles are easy to model: highly expressive voices are still challenging even to recent TTS architectures since there seems to be a trade-off between expressiveness in a generated audio and its signal quality. In this paper, we present a set of techniques that can be leveraged to enhance the signal quality of a highly-expressive voice without the use of additional data. The proposed techniques include: tuning the autoregressive loop's granularity during training; using Generative Adversarial Networks in acoustic modelling; and the use of Variational Auto-Encoders in both the acoustic model and the neural vocoder. We show that, when combined, these techniques greatly closed the gap in perceived naturalness between the baseline system and recordings by 39% in terms of MUSHRA scores for an expressive celebrity voice.


title: Playing with a Quantum Computer

abstract: The high public attention given to quantum computing shows that it is perceived as an interesting topic. We want to utilize this motivating effect for the teaching and learning of quantum physics. Specifically, we want to take advantage of the access to real quantum computers, which various providers make available free of charge. We show a direct and straightforward way to use quantum computers in an introductory course on quantum physics. We use an algorithm that solves a simple and easily understandable problem while providing a quantum advantage. The algorithm we propose is a simple game in which the use of quantum physics offers a winning advantage. The game is called Quantum Penny Flip and was proposed by David A. Meyer back in 1999. It can be easily reformulated to be described by quantum gates. We can therefore use it to teach the programming of a quantum computer. We demonstrate its implementation in IBM's Quantum Composer.


title: Assessing time series irreversibility through micro-scale trends

abstract: Time irreversibility, defined as the lack of invariance of the statistical properties of a system or time series under the operation of time reversal, has received an increasing attention during the last decades, thanks to the information it provides about the mechanisms underlying the observed dynamics. Following the need of analysing real-world time series, many irreversibility metrics and tests have been proposed, each one associated with different requirements in terms of e.g. minimum time series length or computational cost. We here build upon previously proposed tests based on the concept of permutation patterns, but deviating from them through the inclusion of information about the amplitude of the signal and how this evolves over time. We show, by means of synthetic time series, that the results yielded by this method are complementary to the ones obtained by using permutation patterns alone, thus suggesting that "one irreversibility metric does not fit all" We further apply the proposed metric to the analysis of two real-world data sets.


title: A simple lower bound for ARRIVAL

abstract: The ARRIVIAL problem introduced by Dohrau, G\"artner, Kohler, Matou\v{s}ek and Welzl concerns a train moving on a directed graph proceeding along outward edges according to the position of 'switches' at each vertex, which in turn are toggled whenever the train passes through them. The problem asks whether the train every reaches a designated destination vertex. It is known that ARRIVAL is contained in UP $\cap$ coUP, while the previously best published lower bound is that it is NL-hard. In this note we provide a simple reduction to the $\mathsf{DIGICOMP}_\mathsf{EXP}$ problem considered by Aaronson. It follows in particular that ARRIVAL is both CC-hard and PL-hard.


title: Towards artificially intelligent recycling Improving image processing for waste classification

abstract: The ever-increasing amount of global refuse is overwhelming the waste and recycling management industries. The need for smart systems for environmental monitoring and the enhancement of recycling processes is thus greater than ever. Amongst these efforts lies IBM's Wastenet project which aims to improve recycling by using artificial intelligence for waste classification. The work reported in this paper builds on this project through the use of transfer learning and data augmentation techniques to ameliorate classification accuracy. Starting with a convolutional neural network (CNN), a systematic approach is followed for selecting appropriate splitting ratios and for tuning multiple training parameters including learning rate schedulers, layers freezing, batch sizes and loss functions, in the context of the given scenario which requires classification of waste into different recycling types. Results are compared and contrasted using 10-fold cross validation and demonstrate that the model developed achieves a 91.21% test accuracy. Subsequently, a range of data augmentation techniques are then incorporated into this work including flipping, rotation, shearing, zooming, and brightness control. Results show that these augmentation techniques further improve the test accuracy of the final model to 95.40%. Unlike other work reported in the field, this paper provides full details regarding the training of the model. Furthermore, the code for this work has been made open-source and we have demonstrated that the model can perform successful real-time classification of recycling waste items using a standard computer webcam.


title: Leading Edge Vortex Dynamics of Airfoils, Pitching Continuously at High Amplitudes

abstract: Airfoils pitching in the stalled regime have been of keen interest in recent years due to their desirable aerodynamic force characteristics. In this study, we are numerically investigating the unsteady flow past a NACA 0012 airfoil under sinusoidally pitching motion, using a finite volume based sharp interface immersed boundary solver. The flow is investigated in the low Reynolds number regime (Re=3000) for reduced frequencies of 0.1 and 0.5, at three different pivot locations (c/3, c/2 and 2c/3 from the leading edge). The airfoil is subjected to sinusoidal oscillations with its incidence angle varying from 15o to 45o. Leading edge vortices (LEVs) that are formed during the pitching motion dictate the transient aerodynamic characteristics. The flow field data is used to identify individual LEVs in the flow field. The spatio-temporal evolution of LEVs in terms of their strengths are traced throughout the pitching cycle to obtain a quantitative estimate of its evolution. The effect of pivot location and pitching frequency on the vortex evolution and aerodynamic forces is investigated. Change in pitching frequency is found to have a drastic effect on the vortex dynamics. LEV formation is delayed to higher fractional times and LEVs are found to interact and merge into each other at higher frequencies. The maximum lift and drag are also found to increase with pitching frequency. Pivot location is found to have a higher influence on the magnitude and phase of aerodynamic coefficients at higher frequencies. However, at both frequencies, moving the pivot axis aftward is found to delay LEV evolution and cause a phase lag on the aerodynamic forces.


title: NMPC-Based Cooperative Strategy For A Target Pair To Lure Two Attackers Into Collision

abstract: This paper presents a cooperative target defense strategy using nonlinear model-predictive control (NMPC) framework for a two--targets two--attackers (2T2A) game. The 2T2A game consists of two attackers and two targets. Each attacker needs to capture a designated target individually. However, the two targets cooperate to lure the attackers into a collision. We assume that the cooperative target pair do not have perfect knowledge of the attacker states, and hence they estimate the attacker states using an extended Kalman filter (EKF). The NMPC scheme computes closed- loop optimal control commands for the targets while respecting imposed state and control constraints. Theoretical analysis is carried out to determine regions that will lead to the targets' survival, given the initial positions of the attacker and target agents. Numerical simulations are carried out to evaluate the performance of the proposed NMPC- based strategy for different scenarios.


title: Towards Structured Dynamic Sparse Pre-Training of BERT

abstract: Identifying algorithms for computational efficient unsupervised training of large language models is an important and active area of research. In this work, we develop and study a straightforward, dynamic always-sparse pre-training approach for BERT language modeling task, which leverages periodic compression steps based on magnitude pruning followed by random parameter re-allocation. This approach enables us to achieve Pareto improvements in terms of the number of floating-point operations (FLOPs) over statically sparse and dense models across a broad spectrum of network sizes. Furthermore, we demonstrate that training remains FLOP-efficient when using coarse-grained block sparsity, making it particularly promising for efficient execution on modern hardware accelerators.


title: Self-resonant Dark Matter

abstract: We present a novel mechanism for Sommerfeld enhancement for dark matter interactions without the need for light mediators. Considering a model for two-component dark matter with a triple coupling, we find that one of dark matter particles leads to an $u$-channel resonance in dark matter elastic scattering. From the sum of the $u$-channel ladder diagrams, we obtain a Bethe-Salpeter equation with a delay term and identify the Sommerfeld factor from the elastic scattering of two dark matter components for the first time. We discuss the implications of our results for enhancing dark matter self-scattering and annihilation.


title: On Single and Multiple Representations in Dense Passage Retrieval

abstract: The advent of contextualised language models has brought gains in search effectiveness, not just when applied for re-ranking the output of classical weighting models such as BM25, but also when used directly for passage indexing and retrieval, a technique which is called dense retrieval. In the existing literature in neural ranking, two dense retrieval families have become apparent: single representation, where entire passages are represented by a single embedding (usually BERT's [CLS] token, as exemplified by the recent ANCE approach), or multiple representations, where each token in a passage is represented by its own embedding (as exemplified by the recent ColBERT approach). These two families have not been directly compared. However, because of the likely importance of dense retrieval moving forward, a clear understanding of their advantages and disadvantages is paramount. To this end, this paper contributes a direct study on their comparative effectiveness, noting situations where each method under/over performs w.r.t. each other, and w.r.t. a BM25 baseline. We observe that, while ANCE is more efficient than ColBERT in terms of response time and memory usage, multiple representations are statistically more effective than the single representations for MAP and MRR@10. We also show that multiple representations obtain better improvements than single representations for queries that are the hardest for BM25, as well as for definitional queries, and those with complex information needs.


title: Understanding Structural Vulnerability in Graph Convolutional Networks

abstract: Recent studies have shown that Graph Convolutional Networks (GCNs) are vulnerable to adversarial attacks on the graph structure. Although multiple works have been proposed to improve their robustness against such structural adversarial attacks, the reasons for the success of the attacks remain unclear. In this work, we theoretically and empirically demonstrate that structural adversarial examples can be attributed to the non-robust aggregation scheme (i.e., the weighted mean) of GCNs. Specifically, our analysis takes advantage of the breakdown point which can quantitatively measure the robustness of aggregation schemes. The key insight is that weighted mean, as the basic design of GCNs, has a low breakdown point and its output can be dramatically changed by injecting a single edge. We show that adopting the aggregation scheme with a high breakdown point (e.g., median or trimmed mean) could significantly enhance the robustness of GCNs against structural attacks. Extensive experiments on four real-world datasets demonstrate that such a simple but effective method achieves the best robustness performance compared to state-of-the-art models.


title: Modal-Adaptive Gated Recoding Network for RGB-D Salient Object Detection

abstract: The multi-modal salient object detection model based on RGB-D information has better robustness in the real world. However, it remains nontrivial to better adaptively balance effective multi-modal information in the feature fusion phase. In this letter, we propose a novel gated recoding network (GRNet) to evaluate the information validity of the two modes, and balance their influence. Our framework is divided into three phases: perception phase, recoding mixing phase and feature integration phase. First, A perception encoder is adopted to extract multi-level single-modal features, which lays the foundation for multi-modal semantic comparative analysis. Then, a modal-adaptive gate unit (MGU) is proposed to suppress the invalid information and transfer the effective modal features to the recoding mixer and the hybrid branch decoder. The recoding mixer is responsible for recoding and mixing the balanced multi-modal information. Finally, the hybrid branch decoder completes the multi-level feature integration under the guidance of an optional edge guidance stream (OEGS). Experiments and analysis on eight popular benchmarks verify that our framework performs favorably against 9 state-of-art methods.


title: Identification of Incomplete Preferences

abstract: We provide a sharp identification region for discrete choice models where consumers' preferences are not necessarily complete and only aggregate choice data is available. Behavior is modeled using an upper and a lower utility for each alternative so that non-comparability can arise. The identification region places intuitive bounds on the probability distribution of upper and lower utilities. We show that the existence of an instrumental variable can be used to reject the hypothesis that the preferences of all consumers are complete. We apply our methods to data from the 2018 mid-term elections in Ohio.


title: Random Subspace Mixture Models for Interpretable Anomaly Detection

abstract: We present a new subspace-based method to construct probabilistic models for high-dimensional data and highlight its use in anomaly detection. The approach is based on a statistical estimation of probability density using densities of random subspaces combined with geometric averaging. In selecting random subspaces, equal representation of each attribute is used to ensure correct statistical limits. Gaussian mixture models (GMMs) are used to create the probability densities for each subspace with techniques included to mitigate singularities allowing for the ability to handle both numerical and categorial attributes. The number of components for each GMM is determined automatically through Bayesian information criterion to prevent overfitting. The proposed algorithm attains competitive AUC scores compared with prominent algorithms against benchmark anomaly detection datasets with the added benefits of being simple, scalable, and interpretable.


title: On the Cauchy problem of defocusing mKdV equation with finite density initial data: long time asymptotics in soliton-less regions

abstract: We investigate the long-time asymptotics for the solutions to the Cauchy problem of defocusing modified Kortweg-de Vries (mKdV) equation with finite density initial data. The present paper is the subsequent work of our previous paper [arXiv:2108.03650], which gives the soliton resolution for the defocusing mKdV equation in the central asymptotic sector $\{(x,t): \vert \xi \vert<6\}$ with $\xi:=x/t$. In the present paper, via the Riemann-Hilbert (RH) problem associated to the Cauchy problem, the long-time asymptotics in the soliton-less regions $\{(x,t): \vert \xi \vert>6, |\xi|=\mathcal{O}(1)\}$ for the defocusing mKdV equation are further obtained. It is shown that the leading term of the asymptotics are in compatible with the ``background solution'' and the error terms are derived via rigorous analysis.


title: The geometry of spectral interlacing

abstract: We provide a detailed description of the maps associated with spectral interlacing, for rank one perturbations and bordering of symmetric and Hermitian matrices. The arguments rely on standard techniques of nonlinear analysis.


title: Dark energy interactions near the galactic centre

abstract: We investigate scalar-tensor theories, motivated by dark energy models, in the strong gravity regime around the black hole at the centre of our galaxy. In such theories general relativity is modified since the scalar field couples to matter. We consider the most general conformal and disformal couplings of the scalar field to matter to study the orbital behavior of the nearby stars around the galactic star center $Sgr A^{*}$. Markov Chain Monte Carlo (MCMC) simulation yields a bound on the parameters of the couplings of the scalar field to matter. Using Bayesian Analysis yields the first constraints on such theories in the strong gravity regime.


title: Two-level Systems Coupled to Graphene plasmons: A Lindblad equation approach

abstract: In this paper we review the theory of open quantum systems and macroscopic quantum electrodynamics, providing a self-contained account of many aspects of these two theories. The former is presented in the context of a qubit coupled to a electromagnetic thermal bath, the latter is presented in the context of a quantization scheme for surface-plasmon polaritons (SPPs) in graphene based on Langevin noise currents. This includes a calculation of the dyadic Green's function (in the electrostatic limit) for a Graphene sheet between two semi-infinite linear dieletric media, and its subsequent application to the construction of SPP creation and annihilation operators. We then bring the two fields together and discuss the entanglement of two qubits in the vicinity of a graphene sheet which supports SPPs. The two qubits communicate with each other via the emission and absorption of SPPs. We find that a Sch\"odinger cat state involving the two qubits can be partially protected from decoherence by taking advantage of the dissipative dynamics in graphene. A comparison is also drawn between the dynamics at zero temperature, obtained via Schrodinger's equation, and at finite temperature, obtained using the Lindblad equation.


title: Mathematical modelling by help of category theory: models and relations between them

abstract: The growing complexity of modern practical problems puts high demands on the mathematical modelling. Given that various models can be used for modelling one physical phenomenon, the role of model comparison and model choice becomes particularly important. Methods for model comparison and model choice typically used in practical applications nowadays are computation-based, and thus, time consuming and computationally costly. Therefore, it is necessary to develop other approaches for working abstractly, i.e. without computations, with mathematical models. The abstract description of mathematical models can be achieved by help of abstract mathematics, implying formalisation of models and relations between them. In this paper, a category theory-based approach to mathematical modelling is proposed. On this way, mathematical models are formalised in the language of categories, relations between the models are formally defined, as well as several practically relevant properties are introduced on the level of categories. Finally, an illustrative example is presented underlying how the category-theory based approach can be used in practice. Further, all constructions presented in this paper are also discussed from the modelling point of view by making explicit the link to concrete modelling scenarios.


title: Code Perfumes: Reporting Good Code to Encourage Learners

abstract: Block-based programming languages like Scratch enable children to be creative while learning to program. Even though the block-based approach simplifies the creation of programs, learning to program can nevertheless be challenging. Automated tools such as linters therefore support learners by providing feedback about potential bugs or code smells in their programs. Even when this feedback is elaborate and constructive, it still represents purely negative criticism and by construction ignores what learners have done correctly in their programs. In this paper we introduce an orthogonal approach to linting: We complement the criticism produced by a linter with positive feedback. We introduce the concept of code perfumes as the counterpart to code smells, indicating the correct application of programming practices considered to be good. By analysing not only what learners did wrong but also what they did right we hope to encourage learners, to provide teachers and students a better understanding of learners' progress, and to support the adoption of automated feedback tools. Using a catalogue of 25 code perfumes for Scratch, we empirically demonstrate that these represent frequent practices in Scratch, and we find that better programs indeed contain more code perfumes.


title: Sklyanin algebras and a cubic root of 1

abstract: We consider Sklyanin algebras $S$ with 3 generators, which are quadratic algebras over a field $\K$ with $3$ generators $x,y,z$ given by $3$ relations $pxy+qyx+rzz=0$, $pyz+qzy+rxx=0$ and $pzx+qxz+ryy=0$, where $p,q,r\in\K$. This class of algebras has enjoyed much attention. In particular, using tools from algebraic geometry Artin, Tate and Van Den Berg \cite{ATV2} showed that if at least two of the parameters $p$, $q$ and $r$ are non-zero and at least two of three numbers $p^3$, $q^3$ and $r^3$ are distinct, then $S$ is Artin--Schelter regular. More specifically, $S$ is Koszul and has the same Hilbert series as the algebra of commutative polynomials in 3 indeterminates. It has became commonly accepted that it is impossible to achieve the same objective by purely algebraic and combinatorial means like the Gr\"obner basis technique. The authors have previously dispelled this belief. However our previous proof was no less complicated than the one based on algebraic geometry. It used a construcion of a Gr\"obner basis in a suitable one-sided module over $S$ and had quite a number of cases to consider. In this paper we exhibit a linear substitution after which it becomes possible to determine the leading monomials of a reduced Gr\"obner basis for the ideal of relations of $S$ itself (without passing to a module). We also find out explicitly (in terms of parameters) which Sklyanin algebras are isomorphic. The only drawback of the new technique is that it fails if the characteristic of the ground field equals 3.

categories: math.RA math-ph math.AG math.GR math.MP math.RT


title: On the cohomology of the Ree groups and kernels of exceptional isogenies

abstract: Let $G$ be a simple, simply connected algebraic group over an algebraically closed field $k$ of characteristic $p>0$. Let $\sigma : G \rightarrow G$ be a surjective endomorphism of $G$ such that the fixed point set $G(\sigma)$ is a Suzuki or Ree group. Then, let $G_{\sigma}$ denote the scheme-theoretic kernel of $\sigma.$ Using methods of Jantzen and Bendel-Nakano-Pillen, we compute the $1$-cohomology for the Frobenius kernels with coefficients in the induced modules, $H^{1}(G_{\sigma}, H^{0}(\lambda))$, and the $1$-cohomology for the Frobenius kernels with coefficients in the simple modules, $H^{1}(G_{\sigma}, L(\lambda))$ for the Suzuki and Ree groups. Moreover, we improve the known bounds for identifying extensions for the Ree groups of type $F_4$ with the ones for the algebraic group.


title: The Lang-Trotter Conjecture for the elliptic curve $y^2=x^3+Dx$

abstract: Let $E$ be an elliptic curve over $\mathbb{Q}.$ Let $a_p$ denote the trace of the Frobenius endomorphism at a rational prime $p$. For a fixed integer $r,$ define the prime-counting function as $\pi_{E,r}(x):=\sum_{p\leq x,p\nmid \Delta_E,a_p=r}1$. The Lang-Trotter Conjecture predicts that $$\pi_{E,r}(x)=C_{E,r}\cdot \frac{\sqrt{x}}{{\rm log}x}+o(\frac{\sqrt{x}}{{\rm log}x})$$ as $x\longrightarrow \infty,$ where $C_{E,r}$ is a specific non-negative constant. The Hardy-Littlewood Conjecture gives a similar asymptotic formula as above for the number of primes of the form $ax^2+bx+c$. We establish a relationship between the Hardy-Littlewood Conjecture and the Lang-Trotter Conjecture for the elliptic curve $y^2=x^3+Dx.$ We show that the Hardy-Littlewood Conjecture implies the Lang-Trotter Conjecture for $y^2=x^3+Dx.$ Conversely, if the Lang-Trotter Conjecture holds for some $D$ and $2r$ (for $y^2=x^3+Dx, p\nmid D, a_p$ is always even) with positive constant $C_{E,2r},$ then the polynomial $x^2+r^2$ represents infinitely many primes. For a prime $p$, if $a_p=2r$, then $p$ is necessarily of the form $x^2+r^2$. Fixing $r$ and $D$, and assuming that the Hardy-Littlewood Conjecture holds, we obtain the density of the primes with $a_p=2r$ inside the set of primes of the form $x^2+r^2$. In some cases, the density is $1/4$, which is a natural expectation, but it fails to be true for all $D$. In particular, we give a full list of $D$ and $r$ when there is no prime $p$ for $a_p=2r$.


title: Combined sensitivity of JUNO and KM3NeT/ORCA to the neutrino mass ordering

abstract: This article presents the potential of a combined analysis of the JUNO and KM3NeT/ORCA experiments to determine the neutrino mass ordering. This combination is particularly interesting as it significantly boosts the potential of either detector, beyond simply adding their neutrino mass ordering sensitivities, by removing a degeneracy in the determination of $\Delta m_{31}^2$ between the two experiments when assuming the wrong ordering. The study is based on the latest projected performances for JUNO, and on simulation tools using a full Monte Carlo approach to the KM3NeT/ORCA response with a careful assessment of its energy systematics. From this analysis, a $5\sigma$ determination of the neutrino mass ordering is expected after 6 years of joint data taking for any value of the oscillation parameters. This sensitivity would be achieved after only 2 years of joint data taking assuming the current global best-fit values for those parameters for normal ordering.


title: Ghosts without runaway

abstract: We present a simple class of mechanical models where a canonical degree of freedom interacts with another one with a negative kinetic term, i.e. with a ghost. We prove analytically that the classical motion of the system is completely stable for all initial conditions, notwithstanding that the conserved Hamiltonian is unbounded from below and above. This is fully supported by numerical computations. Systems with negative kinetic terms often appear in modern cosmology, quantum gravity and high energy physics and are usually deemed as unstable. Our result demonstrates that for mechanical systems this common lore can be too naive and that living with ghosts can be stable.


title: Diachronic Analysis of German Parliamentary Proceedings: Ideological Shifts through the Lens of Political Biases

abstract: We analyze bias in historical corpora as encoded in diachronic distributional semantic models by focusing on two specific forms of bias, namely a political (i.e., anti-communism) and racist (i.e., antisemitism) one. For this, we use a new corpus of German parliamentary proceedings, DeuPARL, spanning the period 1867--2020. We complement this analysis of historical biases in diachronic word embeddings with a novel measure of bias on the basis of term co-occurrences and graph-based label propagation. The results of our bias measurements align with commonly perceived historical trends of antisemitic and anti-communist biases in German politics in different time periods, thus indicating the viability of analyzing historical bias trends using semantic spaces induced from historical corpora.


title: An ML-style Record Calculus with Extensible Records

abstract: In this work, we develop a polymorphic record calculus with extensible records. Extensible records are records that can have new fields added to them, or preexisting fields removed from them. We also develop a static type system for this calculus and a sound and complete type inference algorithm. Most ML-style polymorphic record calculi that support extensible records are based on row variables. We present an alternative construction based on the polymorphic record calculus developed by Ohori. Ohori based his polymorphic record calculus on the idea of kind restrictions. This allowed him to express polymorphic operations on records such as field selection and modification. With the addition of extensible types, we were able to extend Ohori's original calculus with other powerful operations on records such as field addition and removal.


title: Sequential Subspace Optimization for Quasar-Convex Optimization Problems with Inexact Gradient

abstract: It is well-known that accelerated gradient first-order methods possess optimal complexity estimates for the class of convex smooth minimization problems. In many practical situations it makes sense to work with inexact gradient information. However, this can lead to an accumulation of corresponding inexactness in the theoretical estimates of the rate of convergence. We propose some modification of the Sequential Subspace Optimization Method (SESOP) for minimization problems with quasar-convex functions with inexact gradient. A theoretical result is obtained indicating the absence of accumulation of gradient inexactness. A numerical implementation of the proposed version of the SESOP method and its comparison with the known Similar Triangle Method with an inexact gradient is carried out.


title: Analysis of birefringence and dispersion effects from spacetime-symmetry breaking in gravitational waves

abstract: In this work, we review the effective field theory framework to search for Lorentz and CPT symmetry breaking during the propagation of gravitational waves. The article is written so as to bridge the gap between the theory of spacetime-symmetry breaking and the analysis of gravitational-waves signals detected by ground-based interferometers. The primary physical effects beyond General Relativity that we explore here are dispersion and birefringence of gravitational waves. We discuss their implementation in the open-source LIGO-Virgo algorithm library suite, as well as the statistical method used to perform a Bayesian inference of the posterior probability of the coefficients for symmetry-breaking. We present preliminary results of this work in the form of simulations of modified gravitational waveforms, together with sensitivity studies of the measurements of the coefficients for Lorentz and CPT violation. The findings show the high potential of gravitational wave sources across the sky to probe sensitively for these signals of new physics.


title: Criterion for the functional dissipativity of the Lam\'e operator

abstract: After introducing the concept of functional dissipativity of the Dirichlet problem in a domain $\Omega\subset {\mathbb R}^N$ for systems of partial differential operators of the form $\partial_{h}({\mathscr A}^{hk}(x)\partial_{k})$ (${\mathscr A}^{hk}(x)$ being $m \times m$ matrices with complex valued $L^\infty$ entries), we find necessary and sufficient conditions for the functional dissipativity of the two-dimensional Lam\'e system. As an application of our theory we provide two regularity results for the displacement vector in the $N$-dimensional equilibrium problem, when the body is fixed along its boundary.


title: Multiplicity one theorem over characteristic 2

abstract: It is shown for all local fields $\mathbb{F}$ which are of characteristic different from $2$ that any distribution on $GL_{n+1}(\mathbb{F})$ which is invariant under conjugation by $GL_n(\mathbb{F})$ is also invariant under transposition. In this paper we give an adaptation of the proof of this theorem to fields of characteristic 2.


title: Existence and Non-existence of Ground states of bi-harmonic equations involving constant and degenerate Rabinowitz potentials

abstract: Recently, the authors of the current paper established in [9] the existence of a ground-state solution to the following bi-harmonic equation with the constant potential or Rabinowitz potential: \begin{equation} (-\Delta)^{2}u+V(x)u=f(u)\ \text{in}\ \mathbb{R}^{4}, \end{equation} when the nonlinearity has the special form $f(t)=t(\exp(t^2)-1)$ and $V(x)\geq c>0$ is a constant or the Rabinowitz potential. One of the crucial elements used in [9] is the Fourier rearrangement argument. However, this argument is not applicable if $f(t)$ is not an odd function. Thus, it still remains open whether the above equation with the general critical exponential nonlinearity $f(u)$ admits a ground-state solution even when $V(x)$ is a positive constant. The first purpose of this paper is to develop a Fourier rearrangement-free approach to solve the above problem. More precisely, we will prove that there is a threshold $\gamma^{*}$ such that for any $\gamma\in (0,\gamma^*)$, the above equation with the constant potential $V(x)=\gamma>0$ admits a ground-state solution, while does not admit any ground-state solution for any $\gamma\in (\gamma^{*},+\infty)$. The second purpose of this paper is to establish the existence of a ground-state solution to the above equation with any degenerate Rabinowitz potential $V$ vanishing on some bounded open set. Among other techniques, the proof also relies on a critical Adams inequality involving the degenerate potential which is of its own interest.


title: Context Aware Object Geotagging

abstract: Localization of street objects from images has gained a lot of attention in recent years. We propose an approach to improve asset geolocation from street view imagery by enhancing the quality of the metadata associated with the images using Structure from Motion. The predicted object geolocation is further refined by imposing contextual geographic information extracted from OpenStreetMap. Our pipeline is validated experimentally against the state of the art approaches for geotagging traffic lights.


title: On the uniqueness of Gibbs distributions with a non-negative and subcritical pair potential

abstract: We prove that the distribution of a Gibbs process with non-negative pair potential is uniquely determined as soon as an associated Poisson-driven random connection model (RCM) does not percolate. Our proof combines disagreement coupling in continuum with a coupling of a Gibbs process and a RCM. The improvement over previous uniqueness results is illustrated both in theory and simulations.


title: Finite groups whose real classes have prime-power size

abstract: We present a characterization of the finite groups in which all real classes have prime powers size.


title: "That Star Is Not on the Map": The German Side of the discovery

abstract: Neptune was telescopically discovered by Johan Gottfried Galle and Heinrich Louis d'Arrest in Berlin on 23 September 1846 based on the prediction by Urbain Jean Joseph Le Verrier. The role German astronomers played in the discovery has often been overshadowed by the controversies that erupted in England and France after the discovery. However, their role was crucial, not only in bringing about the discovery in the first place, but also in resolving some of the post-discovery controversies that erupted around priority for the prediction and naming of the planet. German astronomers in Central and Eastern Europe possessed some of the best telescopes of the day and had established themselves over several decades as being at the forefront of observational astronomy. They had produced the star charts that in the end proved indispensable for allowing identification of the planet, while a German publication, Astronomische Nachrichten, published by H. C. Schumacher in Altona, in the vicinity of Hamburg, then part of the Kingdom of Denmark, served as the journal of record. The general neglect of the German part of the story is most strongly attested by the inaccuracies concerning what actually happened on discovery night, which were long propagated in the English-speaking literature. Notably, for 30 years after the discovery, it was not appreciated that there were two observers, while the role of the sky map used was often exaggerated. This chapter sets forth a more complete picture, and in particular emphasises more than previous accounts the critical role of d'Arrest, arguing that he should be celebrated as a co-discoverer. In addition, evidence presented here, much of it not previously available in the English-speaking literature, shows that the name "Neptune" was eventually accepted throughout the scientific community based on German precedents.


title: 3D point cloud segmentation using GIS

abstract: In this paper we propose an approach to perform semantic segmentation of 3D point cloud data by importing the geographic information from a 2D GIS layer (OpenStreetMap). The proposed automatic procedure identifies meaningful units such as buildings and adjusts their locations to achieve best fit between the GIS polygonal perimeters and the point cloud. Our processing pipeline is presented and illustrated by segmenting point cloud data of Trinity College Dublin (Ireland) campus constructed from optical imagery collected by a drone.


title: Skateboard Tricks and Topological Flips

abstract: We study the motion of skateboard flip tricks by modeling them as continuous curves in the group $SO(3)$ of special orthogonal matrices. We show that up to continuous deformation there are only four flip tricks. The proof relies on an analysis of the lift of such curves to the unit 3-sphere. We also derive explicit formulas for a number of tricks and continuous deformations between them.


title: Directional mean dimension and continuum-wise expansive $\mathbb{Z}^k$-actions

abstract: We study directional mean dimension of $\mathbb{Z}^k$-actions (where $k$ is a positive integer). On the one hand, we show that there is a $\mathbb{Z}^2$-action whose directional mean dimension (considered as a $[0,+\infty]$-valued function on the torus) is not continuous. On the other hand, we prove that if a $\mathbb{Z}^k$-action is continuum-wise expansive, then the values of its $(k-1)$-dimensional directional mean dimension are bounded. This is a generalization (with a view towards Meyerovitch and Tsukamoto's theorem on mean dimension and expansive multiparameter actions) of a classical result due to Ma\~n\'e: Any compact metrizable space admitting an expansive homeomorphism (with respect to a compatible metric) is finite-dimensional.


title: Spatio-Temporal Split Learning

abstract: This paper proposes a novel split learning framework with multiple end-systems in order to realize privacypreserving deep neural network computation. In conventional split learning frameworks, deep neural network computation is separated into multiple computing systems for hiding entire network architectures. In our proposed framework, multiple computing end-systems are sharing one centralized server in split learning computation, where the multiple end-systems are with input and first hidden layers and the centralized server is with the other hidden layers and output layer. This framework, which is called as spatio-temporal split learning, is spatially separated for gathering data from multiple end-systems and also temporally separated due to the nature of split learning. Our performance evaluation verifies that our proposed framework shows nearoptimal accuracy while preserving data privacy.


title: MeetSum: Transforming Meeting Transcript Summarization using Transformers!

abstract: Creating abstractive summaries from meeting transcripts has proven to be challenging due to the limited amount of labeled data available for training neural network models. Moreover, Transformer-based architectures have proven to beat state-of-the-art models in summarizing news data. In this paper, we utilize a Transformer-based Pointer Generator Network to generate abstract summaries for meeting transcripts. This model uses 2 LSTMs as an encoder and a decoder, a Pointer network which copies words from the inputted text, and a Generator network to produce out-of-vocabulary words (hence making the summary abstractive). Moreover, a coverage mechanism is used to avoid repetition of words in the generated summary. First, we show that training the model on a news summary dataset and using zero-shot learning to test it on the meeting dataset proves to produce better results than training it on the AMI meeting dataset. Second, we show that training this model first on out-of-domain data, such as the CNN-Dailymail dataset, followed by a fine-tuning stage on the AMI meeting dataset is able to improve the performance of the model significantly. We test our model on a testing set from the AMI dataset and report the ROUGE-2 score of the generated summary to compare with previous literature. We also report the Factual score of our summaries since it is a better benchmark for abstractive summaries since the ROUGE-2 score is limited to measuring word-overlaps. We show that our improved model is able to improve on previous models by at least 5 ROUGE-2 scores, which is a substantial improvement. Also, a qualitative analysis of the summaries generated by our model shows that these summaries and human-readable and indeed capture most of the important information from the transcripts.


title: The Integrated Density of States of the 1D Discrete Anderson-Bernoulli Model at Rational Energies

abstract: We show there is a countable dense set of energies at which the integrated density of states of the 1D discrete Anderson-Bernoulli model can be given explicitly and does not depend on the disorder parameter, provided the latter is above an energy-dependent threshold.


title: Matrix Concentration Inequalities and Free Probability

abstract: A central tool in the study of nonhomogeneous random matrices, the noncommutative Khintchine inequality, yields a nonasymptotic bound on the spectral norm of general Gaussian random matrices $X=\sum_i g_i A_i$ where $g_i$ are independent standard Gaussian variables and $A_i$ are matrix coefficients. This bound exhibits a logarithmic dependence on dimension that is sharp when the matrices $A_i$ commute, but often proves to be suboptimal in the presence of noncommutativity. In this paper, we develop nonasymptotic bounds on the spectrum of arbitrary Gaussian random matrices that can capture noncommutativity. These bounds quantify the degree to which the spectrum of $X$ is captured by that of a noncommutative model $X_{\rm free}$ that arises from free probability theory. This "intrinsic freeness" phenomenon provides a powerful tool for the study of various questions that are outside the reach of classical methods of random matrix theory. Our nonasymptotic bounds are easily applicable in concrete situations, and yield sharp results in examples where the noncommutative Khintchine inequality is suboptimal. When combined with a linearization argument, our bounds imply strong asymptotic freeness for a remarkably general class of Gaussian random matrix models that may be very sparse, have dependent entries, and lack any special symmetries. When combined with a universality principle, our bounds extend beyond the Gaussian setting to general sums of independent random matrices.


title: Accelerating Approximate Aggregation Queries with Expensive Predicates

abstract: Researchers and industry analysts are increasingly interested in computing aggregation queries over large, unstructured datasets with selective predicates that are computed using expensive deep neural networks (DNNs). As these DNNs are expensive and because many applications can tolerate approximate answers, analysts are interested in accelerating these queries via approximations. Unfortunately, standard approximate query processing techniques to accelerate such queries are not applicable because they assume the result of the predicates are available ahead of time. Furthermore, recent work using cheap approximations (i.e., proxies) do not support aggregation queries with predicates. To accelerate aggregation queries with expensive predicates, we develop and analyze a query processing algorithm that leverages proxies (ABae). ABae must account for the key challenge that it may sample records that do not satisfy the predicate. To address this challenge, we first use the proxy to group records into strata so that records satisfying the predicate are ideally grouped into few strata. Given these strata, ABae uses pilot sampling and plugin estimates to sample according to the optimal allocation. We show that ABae converges at an optimal rate in a novel analysis of stratified sampling with draws that may not satisfy the predicate. We further show that ABae outperforms on baselines on six real-world datasets, reducing labeling costs by up to 2.3x.


title: A Dataset for Answering Time-Sensitive Questions

abstract: Time is an important dimension in our physical world. Lots of facts can evolve with respect to time. For example, the U.S. President might change every four years. Therefore, it is important to consider the time dimension and empower the existing QA models to reason over time. However, the existing QA datasets contain rather few time-sensitive questions, hence not suitable for diagnosing or benchmarking the model's temporal reasoning capability. In order to promote research in this direction, we propose to construct a time-sensitive QA dataset. The dataset is constructed by 1) mining time-evolving facts from WikiData and aligning them to their corresponding Wikipedia page, 2) employing crowd workers to verify and calibrate these noisy facts, 3) generating question-answer pairs based on the annotated time-sensitive facts. Our dataset poses challenges in the aspect of both temporal understanding and temporal reasoning. We evaluate different SoTA long-document QA systems like BigBird and FiD on our dataset. The best-performing model FiD can only achieve 46\% accuracy, still far behind the human performance of 87\%. We demonstrate that these models are still lacking the ability to perform consistent temporal reasoning. Therefore, we believe that our dataset could serve as a benchmark to develop NLP models more sensitive to temporal shifts. The dataset and code are released in~\url{https://github.com/wenhuchen/Time-Sensitive-QA}.


title: Mean dimension of Bernstein spaces and universal real flows

abstract: We study the action of translation on the spaces of uniformly bounded continuous functions on the real line which are uniformly band-limited in a compact interval. We prove that two intervals themselves will decide if two spaces are topologically conjugate, while the length of an interval tells the mean dimension of a space. We also investigate universal real flows. We construct a sequence of compact invariant subsets of a space consisting of uniformly bounded smooth one-Lipschitz functions on the real line, which have mean dimension equal to one, such that all real flows can be equivariantly embedded in the translation on their product space. Moreover, we show that the countable self-product of any among them does not satisfy such a universal property. This, on the one hand, presents a more reasonable choice of a universal real flow with a view towards mean dimension, and on the other hand, clarifies a seemingly plausible impression; meanwhile, it refines the previously known results in this direction. Our proof goes through an approach of harmonic analysis. Furthermore, both the universal space that we provide and an embedding mapping which we build for any real flow are explicit.


title: Downlink Resource Allocation in Multiuser Cell-free MIMO Networks with User-centric Clustering

abstract: In this paper, we optimize user scheduling, power allocation and beamforming in distributed multiple-input multiple-output (MIMO) networks implementing user-centric clustering. We study both the coherent and non-coherent transmission modes, formulating a weighted sum rate maximization problem for each; finding the optimal solution to these problems is known to be NP-hard. We use tools from fractional programming, block coordinate descent, and compressive sensing to construct an algorithm that optimizes the beamforming weights and user scheduling and converges in a smooth non-decreasing pattern. Channel state information (CSI) being crucial for optimization, we highlight the importance of employing a low-overhead pilot assignment policy for scheduling problems. In this regard, we use a variant of hierarchical agglomerative clustering, which provides a suboptimal, but feasible, pilot assignment scheme; for our cell-free case, we formulate an area-based pilot reuse factor. Our results show that our scheme provides large gains in the long-term network sum spectral efficiency compared to benchmark schemes such as zero-forcing and conjugate beamforming (with round-robin scheduling) respectively. Furthermore, the results show the superiority of coherent transmission compared to the non-coherent mode under ideal and imperfect CSI for the area-based pilot-reuse factors we consider.

categories: cs.IT cs.NI cs.PF cs.SY eess.SP eess.SY math.IT


title: Towards Efficient Point Cloud Graph Neural Networks Through Architectural Simplification

abstract: In recent years graph neural network (GNN)-based approaches have become a popular strategy for processing point cloud data, regularly achieving state-of-the-art performance on a variety of tasks. To date, the research community has primarily focused on improving model expressiveness, with secondary thought given to how to design models that can run efficiently on resource constrained mobile devices including smartphones or mixed reality headsets. In this work we make a step towards improving the efficiency of these models by making the observation that these GNN models are heavily limited by the representational power of their first, feature extracting, layer. We find that it is possible to radically simplify these models so long as the feature extraction layer is retained with minimal degradation to model performance; further, we discover that it is possible to improve performance overall on ModelNet40 and S3DIS by improving the design of the feature extractor. Our approach reduces memory consumption by 20$\times$ and latency by up to 9.9$\times$ for graph layers in models such as DGCNN; overall, we achieve speed-ups of up to 4.5$\times$ and peak memory reductions of 72.5%.


title: Neuromorphic Processing: A Unifying Tutorial

abstract: All systolic or distributed neuromorphic architectures require power-efficient processing nodes. In this paper, a unifying tutorial is presented which implements multiple neuromorphic processing elements using a systematic analog approach including synapse, neuron and astrocyte models. It is shown that the proposed approach can successfully synthesize multidimensional dynamical systems into analog circuitry with minimum effort.


title: Modular Proofs of Gosper's Identities

abstract: We give unified modular proofs to all of Gosper's identities on the $q$-constant $\Pi_q$. We also confirm Gosper's observation that for any distinct positive integers $n_1,\cdots,n_m$ with $m\geq 3$, $\Pi_{q^{n_1}}$, $\cdots$, $\Pi_{q^{n_m}}$ satisfy a nonzero homogeneous polynomial. Our proofs provide a method to rediscover Gosper's identities. Meanwhile, several results on $\Pi_q$ found by El Bachraoui have been corrected. Furthermore, we illustrate a strategy to construct some of Gosper's identities using hauptmoduls for genus zero congruence subgroups.


title: Newton, entanglement, and the graviton

abstract: Many experiments have recently been proposed to test whether non-relativistic gravitational interactions can generate entanglement. In this note, I consider the extent to which these experiments can test if the graviton exists. Assuming unitarity and Lorentz invariance of the $S$-matrix, I demonstrate that this "Newtonian entanglement" requires the existence of massless bosons, universally coupled to mass, in the Hilbert space of low-energy scattering states. These bosons could be the usual spin-2 gravitons, but in principle there are other possibilities like spin-0 scalar gravitons. I suggest a concept for a more refined experiment to rule these out. The special role of $d=3+1$ spacetime dimensions and the possibility that unitarity is violated by gravity are highlighted.


title: On the relative logarithmic connections and relative residue formula

abstract: We investigate the relative logarithmic connections on a holomorphic vector bundle over a complex analytic family. We give a sufficient condition for the existence of a relative logarithmic connection on a holomorphic vector bundle singular over a relative simple normal crossing divisor. We define the relative residue of relative logarithmic connection and express relative Chern classes of a holomorphic vector bundle in terms of relative residues.


title: Quantifying and Improving Performance of Distributed Deep Learning with Cloud Storage

abstract: Cloud computing provides a powerful yet low-cost environment for distributed deep learning workloads. However, training complex deep learning models often requires accessing large amounts of data, which can easily exceed the capacity of local disks. Prior research often overlooks this training data problem by implicitly assuming that data is available locally or via low latency network-based data storage. Such implicit assumptions often do not hold in a cloud-based training environment, where deep learning practitioners create and tear down dedicated GPU clusters on demand, or do not have the luxury of local storage, such as in serverless workloads. In this work, we investigate the performance of distributed training that leverages training data residing entirely inside cloud storage buckets. These buckets promise low storage costs, but come with inherent bandwidth limitations that make them seem unsuitable for an efficient training solution. To account for these bandwidth limitations, we propose the use of two classical techniques, namely caching and pre-fetching, to mitigate the training performance degradation. We implement a prototype, DELI, based on the popular deep learning framework PyTorch by building on its data loading abstractions. We then evaluate the training performance of two deep learning workloads using Google Cloud's NVIDIA K80 GPU servers and show that we can reduce the time that the training loop is waiting for data by 85.6%-93.5% compared to loading directly from a storage bucket - thus achieving comparable performance to loading data directly from disk - while only storing a fraction of the data locally at a time. In addition, DELI has the potential of lowering the cost of running a training workload, especially on models with long per-epoch training times.


title: Long-Lived Spin-Polarized Intermolecular Exciplex States in Thermally Activated Delayed Fluorescence-Based Organic Light-Emitting Diodes

abstract: Spin-spin interactions in organic light-emitting diodes (OLEDs) based on thermally activated delayed fluorescence (TADF) are pivotal because radiative recombination is largely determined by triplet-to-singlet conversion, also called reverse intersystem crossing (RISC). To explore the underlying process, we apply a spin-resonance spectral hole-burning technique to probe electroluminescence. We find that the triplet exciplex states in OLEDs are highly spin-polarized and show that these states can be decoupled from the heterogeneous nuclear environment as a source of spin dephasing and can even be coherently manipulated on a spin-spin relaxation time scale T2* of 30 ns. Crucially, we obtain the characteristic triplet exciplex spin-lattice relaxation time T1 in the range of 50 us, which far exceeds the RISC time. We conclude that slow spin relaxation rather than RISC is an efficiency-limiting step for intermolecular donor:acceptor systems. Finding TADF emitters with faster spin relaxation will benefit this type of TADF OLEDs.


title: Non-parametric estimation of cumulative (residual) extropy

abstract: Extropy and its properties are explored to quantify the uncertainty. In this paper, we obtain alternative expressions for cumulative residual extropy and negative cumulative extropy. We obtain simple estimators of cumulative (residual) extropy. Asymptotic properties of the proposed estimators are studied. We also present new estimators of cumulative (residual) extropy when the data is right censored. The finite sample performance of the estimators is evaluated through Monte Carlo simulation studies. We use the proposed estimators to analyse different real data sets. Finally, we obtain the relationship between different dynamic and weighted extropy measures and reliability concepts, which leads to several open problems associated with these measures.


title: Continual Backprop: Stochastic Gradient Descent with Persistent Randomness

abstract: The Backprop algorithm for learning in neural networks utilizes two mechanisms: first, stochastic gradient descent and second, initialization with small random weights, where the latter is essential to the effectiveness of the former. We show that in continual learning setups, Backprop performs well initially, but over time its performance degrades. Stochastic gradient descent alone is insufficient to learn continually; the initial randomness enables only initial learning but not continual learning. To the best of our knowledge, ours is the first result showing this degradation in Backprop's ability to learn. To address this degradation in Backprop's plasticity, we propose an algorithm that continually injects random features alongside gradient descent using a new generate-and-test process. We call this the \textit{Continual Backprop} algorithm. We show that, unlike Backprop, Continual Backprop is able to continually adapt in both supervised and reinforcement learning (RL) problems. Continual Backprop has the same computational complexity as Backprop and can be seen as a natural extension of Backprop for continual learning.


title: The exact SL(K+3,C) symmetry of string theory

abstract: By using on-shell recursion relation of string scattering amplitudes (SSA), we show that all n-point SSA of the open bosonic string theory can be expressed in terms of the Lauricella functions. This result extends the previous exact SL(K+3,C) symmetry of the 4-point Lauricella SSA (LSSA) of three tachyons and one arbitrary string states to the whole tree-level open bosonic string theory. Moreover, we present three applications of the SL(K+3,C) symmetry on the SSA. They are the solvability of all SSA in terms of one amplitude, the existence of iteration relations among residues of a given SSA so as to soften its hard scattering behavior and finally the re-derivation of infinite linear relations among hard SSA [12].


title: On the history of the Nekrasov equation

abstract: Appearing in 1921 as an equation for small-amplitude waves on the surface of an infinitely deep liquid, the Nekrasov equation quickly became a source of new results. This manifested itself both in the field of mathematics (theory of nonlinear integral equations of A.I. Nekrasov; 1922, later - of N.N. Nazarov; 1941), and in the field of mechanics (transition to a fluid of finite depth - A.I. Nekrasov; 1927 and refusal on the smallness of the wave amplitude - Yu.P. Krasovskii; 1960).The main task of the author is to find out the prehistory of the Nekrasov equation and to trace the change in approaches to its solution in the context of the nonlinear functional analysis development in the 1940s - 1960s. Close attention will be paid to the contribution of European and Russian mathematicians and mechanics: A.M. Lyapunov, E. Schmidt, T. Levi-Civita, A. Villat, L. Lichtenstein, M.A. Krasnoselskii, N.N. Moiseev, V.V. Pokornyi, etc. In the context of the development of qualitative methods for the Nekrasov equation investigating, the question of the interaction between Voronezh school of nonlinear functional analysis under the guidance of Professor M.A. Krasnoselskii and Rostov school of nonlinear mechanics under the guidance of Professor I.I. Vorovich.


title: Heat conduction in an irregular magnetic field: Part II. Heat transport as a measure of the effective non-integrable volume

abstract: Given the large anisotropy of transport processes in magnetized plasmas, the magnetic field structure can strongly impact heat diffusion: magnetic surfaces and cantori form barriers to transport while chaotic layers and island structures can degrade confinement. When a small but non-zero amount of perpendicular diffusion is included, the structure of the magnetic field becomes less important, allowing pressure gradients to be supported across chaotic regions and island chains. We introduce a metric for the effective volume over which the local parallel diffusion dominates based on the solution to the anisotropic heat diffusion equation. To validate this metric, we consider model fields with a single island chain and a strongly chaotic layer for which analytic predictions of the relative parallel and perpendicular transport can be made. We also analyze critically chaotic fields produced from different sets of perturbations, highlighting the impact of the mode number spectrum on the heat transport. Our results indicate that this metric coincides with the effective volume of non-integrability in the limit $\kappa_{\perp} \rightarrow 0$. We propose that this metric be used to assess the impact of non-integrability on the heat transport in stellarator equilibria.


title: Low-Resource Adaptation of Open-Domain Generative Chatbots

abstract: Recent work building open-domain chatbots has demonstrated that increasing model size improves performance. On the other hand, latency and connectivity considerations dictate the move of digital assistants on the device. Giving a digital assistant like Siri, Alexa, or Google Assistant the ability to discuss just about anything leads to the need for reducing the chatbot model size such that it fits on the user's device. We demonstrate that low parameter models can simultaneously retain their general knowledge conversational abilities while improving in a specific domain. Additionally, we propose a generic framework that accounts for variety in question types, tracks reference throughout multi-turn conversations, and removes inconsistent and potentially toxic responses. Our framework seamlessly transitions between chatting and performing transactional tasks, which will ultimately make interactions with digital assistants more human-like. We evaluate our framework on 1 internal and 4 public benchmark datasets using both automatic (Perplexity) and human (SSA - Sensibleness and Specificity Average) evaluation metrics and establish comparable performance while reducing model parameters by 90%.


title: Discovery of a New Population of Galactic HII Regions with Ionized Gas Velocity Gradients

abstract: We investigate the kinematic properties of Galactic HII regions using radio recombination line (RRL) emission detected by the Australia Telescope Compact Array (ATCA) at 4-10 GHz and the Jansky Very Large Array (VLA) at 8-10 GHz. Our HII region sample consists of 425 independent observations of 374 nebulae that are relatively well isolated from other, potentially confusing sources and have a single RRL component with a high signal-to-noise ratio. We perform Gaussian fits to the RRL emission in position-position-velocity data cubes and discover velocity gradients in 178 (42%) of the nebulae with magnitudes between 5 and 200 m/s/arcsec. About 15% of the sources also have a RRL width spatial distribution that peaks toward the center of the nebula. The velocity gradient position angles appear to be random on the sky with no favored orientation with respect to the Galactic Plane. We craft HII region simulations that include bipolar outflows or solid body rotational motions to explain the observed velocity gradients. The simulations favor solid body rotation since, unlike the bipolar outflow kinematic models, they are able to produce both the large, > 40 m/s/arcsec, velocity gradients and also the RRL width structure that we observe in some sources. The bipolar outflow model, however, cannot be ruled out as a possible explanation for the observed velocity gradients for many sources in our sample. We nevertheless suggest that most HII region complexes are rotating and may have inherited angular momentum from their parent molecular clouds.


title: Discovery of Methanimine Megamasers Toward Compact Obscured Galaxy Nuclei

abstract: We present the first search for the 5.29 GHz methanimine($\rm{CH}_2\rm{NH}$) $1_{10}-1_{11}$ transition toward a sample of galaxy nuclei. We target seven galaxies that host Compact Obscured Nuclei (CONs) with the Karl G. Jansky Very Large Array. These galaxies are characterized by Compton-thick cores. $\rm{CH}_2\rm{NH}$ emission is detected toward six CONs. The brightness temperatures measured toward Arp220 indicate maser emission. Isotropic luminosities of the $\rm{CH}_2\rm{NH}$ transition, from all sources where it is detected, exceed 1 L$_{\odot}$ and thus may be considered megamasers. We also detect formaldehyde ($\rm{H}_2\rm{CO}$) emission toward three CONs. The isotropic $\rm{CH}_2\rm{NH}$ luminosities are weakly correlated with the infrared luminosity of the host galaxy and strongly correlated with OH megamaser luminosities from the same galaxies. Non-LTE radiative transfer models suggest that the maser is pumped by the intense mm/submm radiation field of the CONs. Our study suggests that $\rm{CH}_2\rm{NH}$ megamasers are linked to the nuclear processes within 100 pc of the Compton Thick nucleus within CONs.


title: FlipDA: Effective and Robust Data Augmentation for Few-Shot Learning

abstract: Most previous methods for text data augmentation are limited to simple tasks and weak baselines. We explore data augmentation on hard tasks (i.e., few-shot natural language understanding) and strong baselines (i.e., pretrained models with over one billion parameters). Under this setting, we reproduced a large number of previous augmentation methods and found that these methods bring marginal gains at best and sometimes degrade the performance much. To address this challenge, we propose a novel data augmentation method FlipDA that jointly uses a generative model and a classifier to generate label-flipped data. Central to the idea of FlipDA is the discovery that generating label-flipped data is more crucial to the performance than generating label-preserved data. Experiments show that FlipDA achieves a good tradeoff between effectiveness and robustness -- it substantially improves many tasks while not negatively affecting the others.


title: Self-Calibrating the Look-Elsewhere Effect: Fast Evaluation of the Statistical Significance Using Peak Heights

abstract: In experiments where one searches a large parameter space for an anomaly, one often finds many spurious noise-induced peaks in the likelihood. This is known as the look-elsewhere effect, and must be corrected for when performing statistical analysis. This paper introduces a method to calibrate the false alarm probability (FAP), or $p$-value, for a given dataset by considering the heights of the highest peaks in the likelihood. In the simplest form of self-calibration, the look-elsewhere-corrected $\chi^2$ of a physical peak is approximated by the $\chi^2$ of the peak minus the $\chi^2$ of the highest noise-induced peak. Generalizing this concept to consider lower peaks provides a fast method to quantify the statistical significance with improved accuracy. In contrast to alternative methods, this approach has negligible computational cost as peaks in the likelihood are a byproduct of every peak-search analysis. We apply to examples from astronomy, including planet detection, periodograms, and cosmology.

categories: astro-ph.IM astro-ph.CO astro-ph.EP physics.data-an stat.AP


title: Natural Selection Rules: New Positivity Bounds for Massive Spinning Particles

abstract: We derive new effective field theory (EFT) positivity bounds on the elastic $2\to2$ scattering amplitudes of massive spinning particles from the standard UV properties of unitarity, causality, locality and Lorentz invariance. By bounding the $t$ derivatives of the amplitude (which can be represented as angular momentum matrix elements) in terms of the total ingoing helicity, we derive stronger unitarity bounds on the $s$- and $u$-channel branch cuts which determine the dispersion relation. In contrast to previous positivity bounds, which relate the $t$-derivative to the forward-limit EFT amplitude with no $t$ derivatives, our bounds establish that the $t$-derivative alone must be strictly positive for sufficiently large helicities. Consequently, they provide stronger constraints beyond the forward limit and can be used to constrain dimension-6 interactions with a milder assumption about the high-energy growth of the UV amplitude.


title: Chaos in the vicinity of a singularity in the Three-Body Problem: The equilateral triangle experiment in the zero angular momentum limit

abstract: We present numerical simulations of the gravitational three-body problem, in which three particles lie at rest close to the vertices of an equilateral triangle. In the unperturbed problem, the three particles fall towards the center of mass of the system to form a three-body collision, or singularity, where the particles overlap in space and time. By perturbing the initial positions of the particles, we are able to study chaos in the vicinity of the singularity. Here we cover both the singular region close to the unperturbed configuration and the binary-single scattering regime where one side of the triangle is very short compared to the other two. We make phase space plots to study the regular and ergodic subsets of our simulations and compare them with the outcomes expected from the statistical escape theory of the three-body problem. We further provide fits to the ergodic subset to characterize the properties of the left-over binaries. We identify the discrepancy between the statistical theory and the simulations in the regular subset of interactions, which only exhibits weak chaos. As we decrease the scale of the perturbations in the initial positions, the phase space becomes entirely dominated by regular interactions, according to our metric for chaos. Finally, we show the effect of general relativity corrections by simulating the same scenario with the inclusion of post-Newtonian corrections to the equations of motion.

categories: nlin.CD astro-ph.EP astro-ph.HE astro-ph.SR


title: Parametrization of the Relative Amplitude of Geomagnetic and Askaryan Radio Emission from Cosmic-Ray Air Showers using CORSIKA/CoREAS Simulations (ICRC2021)

abstract: Cosmic rays are messengers from highly energetic events in the Universe. These rare ultra-high-energy particles can be detected efficiently and in an affordable way using large arrays of radio antennas. Linearly polarized geomagnetic emission is the dominant emission mechanism produced when charged particles in air showers get deflected in the Earth's magnetic field. The sub-dominant Askaryan emission is radially polarized and produced due to the time-varying negative-charge excess in the shower front. The relative amplitude of these two emission components depends on various air shower parameters, such as the arrival direction and the depth of the shower maximum. We studied these dependencies using CoREAS simulations of the radio emission from air showers at the South Pole using a star-shaped antenna layout. On the one hand, the parametrization of the Askaryan-to-geomagnetic ratio can be used as input for a more accurate reconstruction of the shower energy. On the other hand, if measured precisely enough, this ratio may provide a new method to reconstruct the atmospheric depth of the shower maximum.


title: Dual Path Learning for Domain Adaptation of Semantic Segmentation

abstract: Domain adaptation for semantic segmentation enables to alleviate the need for large-scale pixel-wise annotations. Recently, self-supervised learning (SSL) with a combination of image-to-image translation shows great effectiveness in adaptive segmentation. The most common practice is to perform SSL along with image translation to well align a single domain (the source or target). However, in this single-domain paradigm, unavoidable visual inconsistency raised by image translation may affect subsequent learning. In this paper, based on the observation that domain adaptation frameworks performed in the source and target domain are almost complementary in terms of image translation and SSL, we propose a novel dual path learning (DPL) framework to alleviate visual inconsistency. Concretely, DPL contains two complementary and interactive single-domain adaptation pipelines aligned in source and target domain respectively. The inference of DPL is extremely simple, only one segmentation model in the target domain is employed. Novel technologies such as dual path image translation and dual path adaptive segmentation are proposed to make two paths promote each other in an interactive manner. Experiments on GTA5$\rightarrow$Cityscapes and SYNTHIA$\rightarrow$Cityscapes scenarios demonstrate the superiority of our DPL model over the state-of-the-art methods. The code and models are available at: \url{https://github.com/royee182/DPL}


title: Graph2MDA: a multi-modal variational graph embedding model for predicting microbe-drug associations

abstract: Accumulated clinical studies show that microbes living in humans interact closely with human hosts, and get involved in modulating drug efficacy and drug toxicity. Microbes have become novel targets for the development of antibacterial agents. Therefore, screening of microbe-drug associations can benefit greatly drug research and development. With the increase of microbial genomic and pharmacological datasets, we are greatly motivated to develop an effective computational method to identify new microbe-drug associations. In this paper, we proposed a novel method, Graph2MDA, to predict microbe-drug associations by using variational graph autoencoder (VGAE). We constructed multi-modal attributed graphs based on multiple features of microbes and drugs, such as molecular structures, microbe genetic sequences, and function annotations. Taking as input the multi-modal attribute graphs, VGAE was trained to learn the informative and interpretable latent representations of each node and the whole graph, and then a deep neural network classifier was used to predict microbe-drug associations. The hyperparameter analysis and model ablation studies showed the sensitivity and robustness of our model. We evaluated our method on three independent datasets and the experimental results showed that our proposed method outperformed six existing state-of-the-art methods. We also explored the meaningness of the learned latent representations of drugs and found that the drugs show obvious clustering patterns that are significantly consistent with drug ATC classification. Moreover, we conducted case studies on two microbes and two drugs and found 75\%-95\% predicted associations have been reported in PubMed literature. Our extensive performance evaluations validated the effectiveness of our proposed method.\


title: Optimality and complexity of classification by random projection

abstract: The generalization error of a classifier is related to the complexity of the set of functions among which the classifier is chosen. We study a family of low-complexity classifiers consisting of thresholding a random one-dimensional feature. The feature is obtained by projecting the data on a random line after embedding it into a higher-dimensional space parametrized by monomials of order up to k. More specifically, the extended data is projected n-times and the best classifier among those n, based on its performance on training data, is chosen. We show that this type of classifier is extremely flexible, as it is likely to approximate, to an arbitrary precision, any continuous function on a compact set as well as any boolean function on a compact set that splits the support into measurable subsets. In particular, given full knowledge of the class conditional densities, the error of these low-complexity classifiers would converge to the optimal (Bayes) error as k and n go to infinity. On the other hand, if only a training dataset is given, we show that the classifiers will perfectly classify all the training points as k and n go to infinity. We also bound the generalization error of our random classifiers. In general, our bounds are better than those for any classifier with VC dimension greater than O (ln n) . In particular, our bounds imply that, unless the number of projections n is extremely large, there is a significant advantageous gap between the generalization error of the random projection approach and that of a linear classifier in the extended space. Asymptotically, as the number of samples approaches infinity, the gap persists for any such n. Thus, there is a potentially large gain in generalization properties by selecting parameters at random, rather than optimization.


title: yupi: Generation, Tracking and Analysis of Trajectory data in Python

abstract: The study of trajectories is often a core task in several research fields. In environmental modelling, trajectories are crucial to study fluid pollution, animal migrations, oil slick patterns or land movements. In this contribution, we address the lack of standardization and integration existing in current approaches to handle trajectory data. Within this scenario, challenges extend from the extraction of a trajectory from raw sensor data to the application of mathematical tools for modeling or making inferences about populations and their environments. This work introduces a generic framework that addresses the problem as a whole, i.e., a software library to handle trajectory data. It contains a robust tracking module aiming at making data acquisition handy, artificial generation of trajectories powered by different stochastic models to aid comparisons among experimental and theoretical data, a statistical kit for analyzing patterns in groups of trajectories and other resources to speed up pre-processing of trajectory data. It is worth emphasizing that this library does not make assumptions about the nature of trajectories (e.g., those from GPS), which facilitates its usage across different disciplines. We validate the software by reproducing key results when modelling dynamical systems related to environmental modelling applications. An example script to facilitate reproduction is presented for each case.


title: Interpolation between noncommutative martingale Hardy and BMO spaces: the case $0<p<1$

abstract: Let $\mathcal{M}$ be a semifinite von Nemann algebra equipped with an increasing filtration $(\mathcal{M}_n)_{n\geq 1}$ of (semifinite) von Neumann subalgebras of $\mathcal{M}$. For $0<p <\infty$, let $\mathsf{h}_p^c(\mathcal{M})$ denote the noncommutative column conditioned martingale Hardy space and $\bmo^c(\M)$ denote the column \lq\lq little\rq\rq \ martingale BMO space associated with the filtration $(\mathcal{M}_n)_{n\geq 1}$. We prove the following real interpolation identity: if $0<p <\infty$ and $0<\theta<1$, then for $1/r=(1-\theta)/p$, \[ \big(\mathsf{h}_p^c(\mathcal{M}), \bmo^c(\mathcal{M})\big)_{\theta, r}=\mathsf{h}_{r}^c(\mathcal{M}), \] with equivalent quasi norms. For the case of complex interpolation, we obtain that if $0<p<q<\infty$ and $0<\theta<1$, then for $1/r =(1-\theta)/p +\theta/q$, \[ \big[\mathsf{h}_p^c(\mathcal{M}), \mathsf{h}_q^c(\mathcal{M})\big]_{\theta}=\mathsf{h}_{r}^c(\mathcal{M}) \] with equivalent quasi norms. These extend previously known results from $p\geq 1$ to the full range $0<p<\infty$. Other related spaces such as spaces of adapted sequences and Junge's noncommutative conditioned $L_p$-spaces are also shown to form interpolation scale for the full range $0<p<\infty$ when either the real method or the complex method is used. Our method of proof is based on a new algebraic atomic decomposition for Orlicz space version of Junge's noncommutative conditioned $L_p$-spaces. We apply these results to derive various inequalities for martingales in noncommutative symmetric quasi-Banach spaces.


title: Giant magnetoanisotropy in the Josephson effect and switching of staggered order in antiferromagnets

abstract: We predict that the amplitude of the Josephson current through an antiferromagnetic weak link changes by several orders of magnitude upon rotation of the N\'eel order parameter characterizing the staggered magnetic order. This occurs due the presence of spin-orbit coupling arising from structural inversion asymmetry which makes the band gap in the antiferromagnet highly sensitive to the staggered order parameter direction. We also demonstrate that when phase-biasing the junction, magnetization dynamics is induced which switches the direction of the N\'eel vector in the antiferromagnet. These results reveal an interesting versatility of antiferromagnetic Josephson junctions as they offer both a large tunability of the supercurrent magnitude via the staggerered magnetization and phase-coherent control over the N\'eel order parameter.


title: Morphology of radio relics II: Properties of polarised emission

abstract: Radio relics are diffuse radio sources in galaxy clusters that are associated with merger shock waves. Detailed observations of radio relics in total intensity and in polarisation show complex structures on kiloparsec scales. The relation between the observed features and the underlying morphology of the magnetic field is not clear. Using three-dimensional magneto-hydrodynamical simulations, we study the polarised emission produced by a shock wave that propagates through a turbulent medium that resembles the intracluster medium. We model the polarised synchrotron emission on the basis of diffusive shock acceleration of cosmic-ray electrons. We find that the synchrotron emission produced in a shocked turbulent medium can reproduce some of the observed features in radio relics. Shock compression can give rise to a high polarisation fraction at the shock front and a partial alignment of the polarisation $E$-vectors with the shock normal. Our work confirms that radio relics can also be formed in an environment with a tangled magnetic field. We also discuss the effect of Faraday Rotation intrinsic to the source, and how our results depend on the angular resolution of observations.


title: Gravitational perturbations of rotating black holes in Lorenz gauge

abstract: Perturbations of Kerr spacetime are typically studied with the Teukolsky formalism, in which a pair of invariant components of the perturbed Weyl tensor are expressed in terms of separable modes that satisfy ordinary differential equations. However, for certain applications it is desirable to construct the full metric perturbation in the Lorenz gauge, in which the linearized Einstein field equations take a manifestly hyperbolic form. Here we obtain a set of Lorenz-gauge solutions to the vacuum field equations in terms of homogeneous solutions to the spin-2, spin-1 and spin-0 Teukolsky equations; and completion pieces that represent perturbations to the mass and angular momentum of the spacetime. The solutions are valid in vacuum Petrov type-D spacetimes that admit a conformal Killing-Yano tensor.


title: Cooling-Induced Vortex Decay in Keplerian Disks

abstract: Vortices are readily produced by hydrodynamical instabilities, such as the Rossby wave instability, in protoplanetary disks. However, large-scale asymmetries indicative of dust-trapping vortices are uncommon in sub-millimeter continuum observations. One possible explanation is that vortices have short lifetimes. In this paper, we explore how radiative cooling can lead to vortex decay. Elliptical vortices in Keplerian disks go through adiabatic heating and cooling cycles. Radiative cooling modifies these cycles and generates baroclinicity that changes the potential vorticity of the vortex. We show that the net effect is typically a spin down, or decay, of the vortex for a sub-adiabatic radial stratification. We perform a series of two-dimensional shearing box simulations, varying the gas cooling (or relaxation) time, $t_{\rm cool}$, and initial vortex strength. We measure the vortex decay half-life, $t_{\rm half}$, and find that it can be roughly predicted by the timescale ratio $t_{\rm cool}/t_{\rm turn}$, where $t_{\rm turn}$ is the vortex turnaround time. Decay is slow in both the isothermal ($t_{\rm cool}\ll t_{\rm turn}$) and adiabatic ($t_{\rm cool}\gg t_{\rm turn}$) limits; it is fastest when $t_{\rm cool}\sim0.1\,t_{\rm turn}$, where $t_{\rm half}$ is as short as $\sim300$ orbits. At tens of au where disk rings are typically found, $t_{\rm turn}$ is likely much longer than $t_{\rm cool}$, potentially placing vortices in the fast decay regime.


title: Metallic and Deconfined Quantum Criticality in Dirac Systems

abstract: Motivated by the physics of spin-orbital liquids, we study a model of interacting Dirac fermions on a bilayer honeycomb lattice at half filling, featuring an explicit global SO(3)$\times$U(1) symmetry. Using large-scale auxiliary-field quantum Monte Carlo (QMC) simulations, we locate two zero-temperature phase transitions as function of increasing interaction strength. First, we observe a continuous transition from the weakly-interacting semimetal to a different semimetallic phase in which the SO(3) symmetry is spontaneously broken and where two out of three Dirac cones acquire a mass gap. The associated quantum critical point can be understood in terms of a Gross-Neveu-SO(3) theory. Second, we subsequently observe a transition towards an insulating phase in which the SO(3) symmetry is restored and the U(1) symmetry is spontaneously broken. While strongly first order at the mean-field level, the QMC data is consistent with a direct and continuous transition. It is thus a candidate for a new type of deconfined quantum critical point that features gapless fermionic degrees of freedom.


title: Dijet impact factor in DIS at next-to-leading order in the Color Glass Condensate

abstract: We compute the next-to-leading order impact factor for inclusive dijet production in deeply inelastic electron-nucleus scattering at small $x_{\rm Bj}$. Our computation, performed in the framework of the Color Glass Condensate effective field theory, includes all real and virtual contributions in the gluon shock wave background of all-twist lightlike Wilson line correlators. We demonstrate explicitly that the rapidity evolution of these correlators, to leading logarithmic accuracy, is described by the JIMWLK Hamiltonian. When combined with the next-to-leading order JIMWLK Hamiltonian, our results for the impact factor improve the accuracy of the inclusive dijet cross section to $\mathcal{O}(\alpha_s^2 \ln(x_f/x_{\rm Bj}))$, where $x_f$ is a rapidity factorization scale. These results are an essential ingredient in assessing the discovery potential of inclusive dijets to uncover the physics of gluon saturation at the Electron-Ion Collider.


title: The cumulative hierarchy in Homotopy Type Theory

abstract: We explore the cumulative hierarchy $V$ defined in Chapter 10 of the HoTT book. We begin by showing how to translate formulas of set theory in HoTT, and proceed to examine which axioms are satisfied in $V$. In particular, we show that $V$ models ZF$^-$ in HoTT+PR, while LEM is required to obtain full ZF. Finally, we attempt to model constructive set theories in V, and although this is achieved for ECST, we only obtain IZF and CZF with LEM.


title: Criticality-Enhanced Quantum Sensing via Continuous Measurement

abstract: Present protocols of criticality-enhanced sensing with open quantum sensors assume direct measurement of the sensor and omit the radiation quanta emitted to the environment, thereby potentially missing valuable information. Here we propose a protocol for criticality-enhanced sensing via continuous observation of the emitted radiation quanta. Under general assumptions, we establish a scaling theory for the global quantum Fisher information of the joint system and environment state at dissipative critical points. We derive universal scaling laws featuring transient and long-time behavior governed by the underlying critical exponents. Importantly, such scaling laws exceed the standard quantum limit and can in principle saturate the Heisenberg limit. To harness such advantageous scaling, we propose a practical sensing scheme based on continuous detection of the emitted quanta as realized experimentally in various quantum-optical setups. In such a scheme a single interrogation corresponds to a (stochastic) quantum trajectory of the open system evolving under the nonunitary dynamics dependent on the parameter to be sensed and the backaction of the continuous measurement. Remarkably, we demonstrate that the associated precision scaling significantly exceeds that based on direct measurement of the critical steady state, thereby establishing the metrological value of the continuous detection of the emitted quanta at dissipative criticality. We illustrate our protocol via counting the photons emitted by the open Rabi model, a paradigmatic model for the study of dissipative phase transition with finite components. Our protocol is applicable to generic quantum-optical open sensors permitting continuous readout, and may find applications at the frontier of quantum sensing, such as the human-machine interface, magnetic diagnosis of heart disease, and zero-field nuclear magnetic resonance.


title: Computing with jets

abstract: We introduce a Macaulay2 package for working with jet schemes. The main method constructs jets of ideals, polynomial rings and their quotients, ring homomorphisms, affine varieties, and (hyper)graphs. The package also includes additional methods to compute principal components and radicals of jets of monomial ideals.


title: q-Fibonacci bicomplex quaternions

abstract: In the paper, we define the $q$-Fibonacci bicomplex quaternions and the $q$-Lucas bicomplex quaternions, respectively. Then, we give some algebraic properties of $q$-Fibonacci bicomplex quaternions and the $q$-Lucas bicomplex quaternions.


title: Comparison Between International Space Station and Airplane Flying

abstract: Inertial forces are not intuitive; therefore, interesting examples are great help for learners. In this paper, we examine simple models for International Space Station (ISS) and airplane flying and identify the inertial forces involved and their effect. In both cases there are centrifugal forces but in the case of the ISS, the centrifugal force is the only major force opposing the gravitational attraction that is pulling it towards the Earth at 0.885 the gravity at the surface of the Earth. While in the case of an airplane, it is a minor component amounting to less than 1% in the best case of flying along the equator west to east. Other inertial forces including the Coriolis force affect the airplane. Another thing we see is that flying with the Earth direction west to east gives a markedly higher forces than flying in the opposite direction. The maintenance of the ISS orientation with respect to Earth requires rotation of ISS around its center of gravity (cg) at the same rate as it revolves around the Earth (4 degrees per minute) while for an airplane, maintaining level flight at a constant speed and height does not require direct intervention because the amount of adjustment is so small and gradual that it is handled by the pilot system as a part of the disturbances that affects the airplane.


title: Chiral Symmetry Breaking and Confinement from an Interacting Ensemble of Instanton-dyons in Two-flavor Massless QCD

abstract: In this work we present the results from numerical simulations of an interacting ensemble of instanton-dyons in the $SU(3)$ gauge group with $N_f=2$ flavors of massless quarks. Dynamical quarks are included via the effective interactions induced by the fermionic determinant evaluated in the subspace of topological zero modes. The eigenvalue spectrum of the Dirac operator is studied at different volumes to extract the chiral condensate and eigenvalue gap, with both observables providing consistent values of the chiral transition temperature $T_c$. We find that a sufficient density of dyons is responsible for generating the confining potential and breaking the chiral symmetry, both of which are compatible with second-order transitions.


title: A Generalized Definition of Fractional Derivative with Applications

abstract: A generalized fractional derivative (GFD) definition is proposed in this work. For a differentiable function that can be expanded by Taylor series, we show that D^Elafa*D^Beta f(t)=D^(Elafa+Beta)f(t). GFD is applied for some functions in which we investigate that GFD coincides with Caputo and Riemann-Liouville fractional derivatives' results. The solutions of Riccati fractional differential equation are simply obtained via GFD. A comparison with other definitions is also discussed. The results show that the proposed definition in this work gives better accuracy than the commonly known conformable derivative definition. Therefore, GFD has some advantages in comparison with other definitions in which a new path is provided for simple analytical solutions of many problems in the context of fractional calculus.


title: Boundary dynamics in competing critical black hole formation

abstract: Expanding upon our previous study of competing critical phenomena in black hole formation, we numerically investigate the behavior of dominant exponents across the boundary separating asymptotically dispersing and collapsing regions in a two-dimensional configuration space of initial data. We find that across the Type II boundary section the dominant exponent remains constant, equal to the reciprocal of Choptuik's well-known quasi-universal value, whereas across the Type I section the exponent noticeably varies. We postulate that this change reflects the existence of a third critical solution in addition to the two primary competing solutions, possibly another member of the family of metastable soliton stars constituting the Type I attractor.


title: Neutrino transport with Monte Carlo method: II. Quantum Kinetic Equations

abstract: Neutrinos have an unique quantum feature as flavor conversions. Recent studies suggested that collective neutrino oscillations play important roles in high-energy astrophysical phenomena. Quantum kinetic equation (QKE) is capable of describing the neutrino flavor conversion, transport and matter collision self-consistently. However, we have experienced many technical difficulties in their numerical implementation. In this paper, we present a new QKE solver based on Monte Carlo (MC) approach. This is an upgraded version of our classical MC neutrino transport solver; in essence, a flavor degree of freedom including mixing state is added into each MC particle. This extension requires updating numerical treatments of collision terms, in particular for scattering processes. We deal with the technical problem by generating a new MC particle at each scattering event. To reduce statistical noise inherent in MC methods, we develop the effective mean free path method. This suppresses a sudden change of flavor state due to collisions without increasing the number of MC particles. We present a suite of code tests to validate these new modules with comparing to the results reported in previous studies. Our QKE-MC solver is developed with fundamentally different philosophy and design from other deterministic- and mesh methods, suggesting that it will be complementary to others, and potentially provide new insights into physical processes of neutrino dynamics.