title: Memory-Efficient Backpropagation through Large Linear Layers

abstract: In modern neural networks like Transformers, linear layers require significant memory to store activations during backward pass. This study proposes a memory reduction approach to perform backpropagation through linear layers. Since the gradients of linear layers are computed by matrix multiplications, we consider methods for randomized matrix multiplications and demonstrate that they require less memory with a moderate decrease of the test accuracy. Also, we investigate the variance of the gradient estimate induced by the randomized matrix multiplication. We compare this variance with the variance coming from gradient estimation based on the batch of samples. We demonstrate the benefits of the proposed method on the fine-tuning of the pre-trained RoBERTa model on GLUE tasks.


title: A bang-bang principle for the conditional expectation vector measure

abstract: We prove a conditional expectation bang-bang principle. Based on properties of the conditional expectation vector measure, we establish that the conditional expectation of a set-valued mapping coincides with the conditional expectation of the set of selections of its extreme points part. As a by-product, we obtain straightforwardly a purification principle.


title: Exciton-Coupled Coherent Magnons in a 2D Semiconductor

abstract: Two-dimensional (2D) magnetic semiconductors feature both tightly-bound excitons with large oscillator strength and potentially long-lived coherent magnons due to the presence of bandgap and spatial confinement. While magnons and excitons are energetically mismatched by orders of magnitude, their coupling can lead to efficient optical access to spin information. Here we report strong magnon-exciton coupling in the 2D van der Waals (vdW) antiferromagnetic (AFM) semiconductor CrSBr. Coherent magnons launched by above-gap excitation modulate the interlayer hybridization, which leads to dynamic modulation of excitonic energies. Time-resolved exciton sensing reveals magnons that can coherently travel beyond 7 micrometer, with coherence time above 5 ns. We observe this exciton-coupled coherent magnons in both even and odd number of layers, with and without compensated magnetization, down to the bilayer limit. Given the versatility of vdW heterostructures, these coherent 2D magnons may be basis for optically accessible magnonics and quantum interconnects.


title: A subsampling approach for Bayesian model selection

abstract: It is common practice to use Laplace approximations to compute marginal likelihoods in Bayesian versions of generalised linear models (GLM). Marginal likelihoods combined with model priors are then used in different search algorithms to compute the posterior marginal probabilities of models and individual covariates. This allows performing Bayesian model selection and model averaging. For large sample sizes, even the Laplace approximation becomes computationally challenging because the optimisation routine involved needs to evaluate the likelihood on the full set of data in multiple iterations. As a consequence, the algorithm is not scalable for large datasets. To address this problem, we suggest using a version of a popular batch stochastic gradient descent (BSGD) algorithm for estimating the marginal likelihood of a GLM by subsampling from the data. We further combine the algorithm with Markov chain Monte Carlo (MCMC) based methods for Bayesian model selection and provide some theoretical results on the convergence of the estimates. Finally, we report results from experiments illustrating the performance of the proposed algorithm.


title: Passage of test particles through oscillating spherically-symmetric dark matter configurations

abstract: Applying the perturbative approach to geodesic equations, we study motion of the test particles in time-dependent spherically symmetric spacetimes created by oscillating dark matter. Assuming the weakness of the gravitational field, we derive general formulas that describe infinite trajectories of the test particles and determine the total deflection angle in the leading order approximation. The obtained formulas are valid for both time-dependent and static matter configurations. Using these results, we calculate the deflection angle of a test particle passing through a spherically symmetric oscillating distribution of a self-gravitating scalar field with a logarithmic potential. It turned out that, in a wide range of amplitudes, oscillations in the deflection angle are sinusoidal and become small for ultrarelativistic particles.


title: Complex Langevin studies of the emergent space-time in the type IIB matrix model

abstract: The type IIB matrix model has been proposed as a non-perturbative definition of superstring theory since 1996. We study a simplified model that describes the late time behavior of the type IIB matrix model non-perturbatively using Monte Carlo methods, and we use the complex Langevin method to overcome the sign problem. We investigate a scenario where the space-time signature changes dynamically from Euclidean at early times to Lorentzian at late times. We discuss the possibility of the emergence of the (3+1)D expanding universe.


title: The Donor of the Black-Hole X-Ray Binary MAXI J1820+070

abstract: We estimate the parameters of the donor of the accreting black-hole binary MAXI J1820+070. The measured values of the binary period, rotational and radial velocities and constraints on the orbital inclination imply the donor is a subgiant with the mass of $M_2\approx 0.49^{+0.10}_{-0.10}M_\odot$ and the radius of $R_2\approx 1.19^{+0.08}_{-0.08}R_\odot$. We re-analyze the previously obtained optical spectrum from the Gran Telescopio Canarias, and found it yields a strict lower limit on the effective temperature of $T>4200$ K. We compile optical and infrared fluxes observed during the quiescence of this system. From the minima $r$ and $i$-band fluxes found in Pan-STARSS1 Data Release 2 pre-discovery imaging and for a distance of $D\approx3$ kpc, reddening of $E(B$--$V)=0.23$ and $R_2\approx{1.11R_\odot}$, we find $T\lesssim4230$ K, very close to the above lower limit. For a larger distance, the temperature can be higher, up to about 4500 K (corresponding to a K5 spectral type, preferred by previous studies) at $D=3.5$ kpc, allowed by the Gaia parallax. We perform evolutionary calculations for the binary system and compare them to the observational constraints. Our model fitting the above temperature and radius constraints at $D\approx 3$ kpc has the mass of $0.4M_\odot$, $T\approx4200$ K and solar metallicity. Two alternative models require $D\gtrsim 3.3$--3.4 kpc at $0.4 M_\odot$, $T\approx4500$ K and half solar metallicity, and $0.5M_\odot$, $T\approx4300$ K and solar metallicity. These models yield mass transfer rates of $\sim\!\!10^{-10}M_\odot$/yr, compatible with those based on the estimated accreted mass of $\approx\!2\times 10^{25}$ g and the time between the 2018 discovery and the 1934 historical outburst.


title: ODSearch: Fast and Resource Efficient On-device Natural Language Search for Fitness Trackers' Data

abstract: Mobile and wearable technologies have promised significant changes to the healthcare industry. Although cutting-edge communication and cloud-based technologies have allowed for these upgrades, their implementation and popularization in low-income countries have been challenging. We propose "ODSearch", an On-device Search framework equipped with a natural language interface for mobile and wearable devices. To implement search, "ODSearch" employs compression and Bloom filter, it provides near real-time search query responses without network dependency. In particular, the Bloom filter reduces the temporal scope of the search and compression reduces the size of the data to be searched. Our experiments were conducted on a mobile phone and smartwatch. We compared "ODSearch" with current state-of-the-art search mechanisms, and it outperformed them on average by 53 times in execution time, 26 times in energy usage, and 2.3% in memory utilization.


title: Comment on the thermal Hall effect in cuprates

abstract: Insights from stripe incommensurabilities and antiferromagnetic stability indicate that the magnetic moments of both host Cu^2+ ions and Cu atoms from electron doping support the thermal Hall effect in cuprates, whereas those of O atoms from hole doping oppose it.


title: Additive polynomial superfunctors and cohomology

abstract: We prove a classification of additive polynomial superfunctors, which allows us to compute some extensions of a superfunctor of the form $F \circ A$ where $F$ is a classical polynomial functor and $A$ is additive. We get a formula which relates these extensions to the classical ones of $F$. A possible generalisation is conjectured at the end.


title: Vibronic mass computation for the $EF$-$GK$-$H\bar{H}$ $^1\Sigma_\text{g}^+$ manifold of molecular hydrogen

abstract: A variational procedure is described for the computation of the non-adiabatic mass-correction tensor applicable for multi-dimensional electronic manifolds. The 30-year-old computations of Wolniewicz, Dressler, and their co-workers are appended with the computed vibronic mass-correction functions corresponding to the $EF$-$GK$-$H\bar{H}$-$S5$-$S6$ $^1\Sigma_\text{g}^+$ manifold of the hydrogen molecule. Initial results are reported for the vibronic energies. Necessary further improvements and further developments are discussed.


title: About the Evolution of the Concept of Heat and Thermal Energy

abstract: The concept of thermal energy was only fully formulated in the late 19th and early 20th centuries. The search for the understanding of calorific phenomena was what provided the construction and development of the concept of thermal energy. In this work we present the evolution of the concept of heat and thermal energy over time, following a chronological order and the evolution of the ideas and principles in question. For this, a vast bibliographical research on the subject was carried out. The theme has been addressed from the ancient Greek philosophers to the advent of modern physics with the elucidation of the phenomenon of thermal radiation.


title: Non-Diophantine arithmetic as the mathematical foundation for quantum field theory

abstract: The problem of infinities in quantum field theory (QRT) is a long standing problem in physics.For solving this problem, different renormalization techniques have been suggested but the problem still persists. Here we suggest another approachto the elimination of infinities in QFT, which is based on non-Diophantine arithmetics - a novel mathematical area that already found useful applications in physics. To achieve this goal, new non-Diophantine arithmetics are constructed and their properties are studied. This allows using these arithmetics for computing integrals describing Feynman diagrams. Although in the conventional QFT these integrals diverge, their non-Diophantine counterparts are convergent and rigorously defined.


title: Elliptic Integrals for Calculation of the Propagation Time of a Signal, Emitted and Received by Satellites on One Orbit

abstract: The propagation time of a signal, emitted by a moving along an elliptical orbit satellite from the GPS (or GLONASS) satellite confi gurations is a very important ingredient of the theory, based on the formalism of the null cone and accounting for the effects of the General Relativity Theory. For the case of satellites, orbiting along a plane elliptic orbit, it has been proved that the propagation time for the signal between the satellites is given by a combination of elliptic integrals of the first, second and third kind. For the more general case of satellites on a space-distributed elliptic orbit, the propagation time is expressed by higher (fourth) order elliptic integrals, which according to the standard theory can be expressed recurrently by means of lower-order elliptic integrals. In the concrete case, the elliptic integrals of the second and the fourth order are expressed by means of a combination of irrational functions and the zero-order elliptic integral in the Legendre form. It has been proved that for the investigated case, second-order elliptic integrals can be expressed by elementary functions.


title: A Hypothesis about the Physical Nature of the Phenomenon of Life (A Contribution to the Discussion of the Paper by G.R. Ivanitskii 21st Century: What Is Life from the Perspective of Physics)

abstract: A macroscopic quantum state that is characterized by an energy gap in the electronic spectrum and occurs in organic structural components of the cell is hypothesized to provide a basis for the phenomenon of life. The width of the band gap has been estimated. The possibility of detecting this energy gap in an experimental study is discussed.


title: The background field method applied to cosmological phase transition

abstract: In this paper, the cosmic phase transition is investigated by background gauge field method. As a continuation of previous our work, some numerical results and graphic solutions at $T\neq 0$ are presented. Hence the mechanism of cosmic phase transition in the early Universe is considered. It is shown that the breaking of symmetry significantly depend on the nonzero temperature and chemical potential. Furthermore, it is the first order of phase transition. Non - restoration of symmetry in hot gauge theories for Cosmology.


title: Three-dimensional AdS black holes in massive-power-Maxwell theor

abstract: Recently, it was shown that the power-Maxwell (PM) theory could remove the singularity of the electric field \cite{PM2}. Motivated by a great interest in three-dimensional black holes and a surge of success in studying massive gravity from both the cosmological and astrophysical points of view, we investigate three-dimensional black hole solutions in de Rham, Gabadadze, and Tolley (dRGT) massive theory of gravity in the presence of PM electrodynamics. First, we extract exact three-dimensional solutions in this theory of gravity. Then we study the geometrical properties of these solutions. Calculating conserved and thermodynamic quantities, we check the validity of the first law of thermodynamics for these black holes. We also examine the stability of these black holes in the context of the canonical ensemble. We continue calculating this kind of black hole's optical features, such as the photon orbit radius, the energy emission rate, and the deflection angle. Considering these optical quantities, finally, we analyze the effective role of the parameters of models on them.


title: Application of Ionizing Radiation for the Enhancement of Hydrogen Yield in Electrolytic Water Decomposition

abstract: During the past decay many experimental configurations to improve the yield of hydrogen by electrolysis. These attempts include tests of different materials for the electrodes. This study proposes the use of {\gamma}-radiation of waste sources of nuclear reactors. The related experimental configuration can be reduced to boxes filled with water exposed to ionizing radiation and capacitors as electrodes as well-known in other technical disciplies. Test measurements based on radiation of a linear accelerator (6 MV beam) provided significant differences of the yield of hydrogen between with and without radiation.


title: Fermi surface segmentation in the helical state of a Rashba superconductor

abstract: We investigate the quasiparticle excitations in the FFLO- type helical state of a superconductor with inversion-symmetry breaking and strong Rashba spin-orbit coupling. We restrict to a state with single finite momentum of Cooper pairs in the helical phase that is determined by minimization of the condensation energy. We derive the dependence of quasiparticle dispersions on the Rashba coupling strength and external field. It leads to a peculiar momentum-space segmentation of the corresponding Rashba Fermi surface sheets which has not yet been observed experimentally. We show that it may be directly visualized by the method of quasiparticle interference that identifies the critical points of the segmented sheets and can map their evolution with field strength, bias voltage and Rashba coupling. We also indicate a strategy how to determine the finite Cooper-pair momentum from experimental quantities. This investigation has the potential for a more detailed microscopic understanding of the helical superconducting state under the influence of Rashba spin-orbit coupling.


title: Variational vibrational states of HCOOH

abstract: Vibrational states of the formic acid molecule are converged using the GENIUSH-Smolyak approach and the potential energy surface taken from [D. Tew and W. Mizukami, J. Phys. Chem. A 120, 9815 (2016)]. The quantum nuclear motion is described by using the $cis$-$trans$ torsional coordinate and eight curvilinear normal coordinates defined with respect to an instantaneous reference configuration changing as a function of the torsional degree of freedom. Harmonic oscillator basis functions are used for the curvilinear normal coordinates, a Fourier basis for the torsional coordinate, and a simple basis pruning condition is combined with a Smolyak integration grid. $Trans$, $cis$, and $delocalized$ vibrational states are reported up to and slightly beyond the isomerization barrier.


title: The Curvature Effect in Gaussian Random Fields

abstract: Random field models are mathematical structures used in the study of stochastic complex systems. In this paper, we compute the shape operator of Gaussian random field manifolds using the first and second fundamental forms (Fisher information matrices). Using Markov Chain Monte Carlo techniques, we simulate the dynamics of these random fields and compute the Gaussian curvature of the parametric space, analyzing how this quantity changes along phase transitions. During the simulation, we have observed an unexpected phenomenon that we called the \emph{curvature effect}, which indicates that a highly asymmetric geometric deformation happens in the underlying parametric space when there are significant increase/decrease in the system's entropy. This asymmetric pattern relates to the emergence of hysteresis, leading to an intrinsic arrow of time along the dynamics.

categories: cs.IT cond-mat.stat-mech cs.CE math.IT physics.comp-ph


title: Weyl and Majorana Spinors as Pure Goldstone Bosons

abstract: General spinors in polar form display a structure that is made up by real scalar degrees of freedom plus six components which can be recognized as Goldstone bosons: in the present paper we show that of all singular spinors, Weyl and Majorana spinors have no real degree of freedom and so that they can be interpreted as pure Goldstone states.


title: Fast Distributed k-Means with a Small Number of Rounds

abstract: We propose a new algorithm for k-means clustering in a distributed setting, where the data is distributed across many machines, and a coordinator communicates with these machines to calculate the output clustering. Our algorithm guarantees a cost approximation factor and a number of communication rounds that depend only on the computational capacity of the coordinator. Moreover, the algorithm includes a built-in stopping mechanism, which allows it to use fewer communication rounds whenever possible. We show both theoretically and empirically that in many natural cases, indeed 1-4 rounds suffice. In comparison with the popular k-means|| algorithm, our approach allows exploiting a larger coordinator capacity to obtain a smaller number of rounds. Our experiments show that the k-means cost obtained by the proposed algorithm is usually better than the cost obtained by k-means||, even when the latter is allowed a larger number of rounds. Moreover, the machine running time in our approach is considerably smaller than that of k-means||. Code for running the algorithm and experiments is available at https://github.com/selotape/distributed_k_means.


title: Powering Stellar Magnetism: Energy Transfers in Cyclic Dynamos of Sun-like Stars

abstract: We use the ASH code to model the convective dynamo of solar-type stars. Based on a series of 15 3-D MHD simulations spanning 4 bins in rotation and mass, we show what mechanisms are at work in these stellar dynamos with and without magnetic cycles and how global stellar parameters affect the outcome. We also derive scaling laws for the differential rotation and magnetic field based on these simulations. We find a weaker trend between differential rotation and stellar rotation rate, ($\Delta \Omega \sim (|\Omega|/\Omega_{\odot})^{0.46}$) in the MHD solutions than in their HD counterpart $(|\Omega|/\Omega_{\odot})^{0.66})$, yielding a better agreement with the observational trends based on power laws. We find that for a fluid Rossby number between $0.15 \lesssim Ro_f \lesssim 0.65$ the solutions possess long magnetic cycle, if $Ro_f \lesssim 0.42$ a short cycle and if $Ro_f \gtrsim 1$ (anti-solar-like differential rotation) a statistically steady state. We show that short-cycle dynamos follow the classical Parker-Yoshimura rule whereas the long-cycle period ones do not. We further demonstrate that the Rossby number dependency of the large-scale surface magnetic field in the simulation ($B_{L,surf} \sim Ro_{f}^{-1.26}$) agrees better with observations ($B_{V} \sim Ro_{s}^{-1.4 \pm 0.1}$) and differs from dynamo scaling based on the global magnetic energy ($B_{bulk} \sim Ro_{f}^{-0.5}$). We also show that up to few percents of the stellar luminosity can be channelled into the star's magnetism, hence providing a large energy reservoir for possible surface eruptive events.


title: The Extreme-ultraviolet Stellar Characterization for Atmospheric Physics and Evolution (ESCAPE) Mission: Motivation and Overview

abstract: The Extreme-ultraviolet Stellar Characterization for Atmospheric Physics and Evolution (ESCAPE) mission is an astrophysics Small Explorer employing ultraviolet spectroscopy (EUV: 80 - 825 \AA\ and FUV: 1280 - 1650 \AA) to explore the high-energy radiation environment in the habitable zones around nearby stars. ESCAPE provides the first comprehensive study of the stellar EUV and coronal mass ejection environments which directly impact the habitability of rocky exoplanets. In a 20 month science mission, ESCAPE will provide the essential stellar characterization to identify exoplanetary systems most conducive to habitability and provide a roadmap for NASA's future life-finder missions. ESCAPE accomplishes this goal with roughly two-order-of-magnitude gains in EUV efficiency over previous missions. ESCAPE employs a grazing incidence telescope that feeds an EUV and FUV spectrograph. The ESCAPE science instrument builds on previous ultraviolet and X-ray instrumentation, grazing incidence optical systems, and photon-counting ultraviolet detectors used on NASA astrophysics, heliophysics, and planetary science missions. The ESCAPE spacecraft bus is the versatile and high-heritage Ball Aerospace BCP Small spacecraft. Data archives will be housed at the Mikulski Archive for Space Telescopes (MAST).


title: The criterion for infinite positron feedback in dynamics of relativistic runaway electron avalanches

abstract: Relativistic runaway electron avalanches (RREA) accelerated by thunderstorm large-scale electric fields are one of the sources of atmospheric gamma radiation. In strong electric fields, RREAs can multiply by the relativistic feedback. Infinite relativistic feedback makes avalanches self-sustainable and hypothetically can cause a terrestrial gamma-ray flash (TGF). This paper introduces a kinetic approach to study the relativistic feedback caused by positrons since positron feedback dominates for the directly observed electric field strengths. With this approach, the criterion for infinite positron feedback within thunderstorms is derived. Discovered criterion allows obtaining the thunderstorm electric field parameters required for infinite positron feedback for any altitude. The possibility of derived thunderstorm conditions is discussed.


title: Risk-based Design of Regular Plane Frames Subject to Damage by Abnormal Events: a Conceptual Study

abstract: Constructed facilities should be robust with respect to the loss of load-bearing elements due to abnormal events. Yet, strengthening structures to withstand such damage has a significant impact on construction costs. Strengthening costs should be justified by the threat and should result in smaller expected costs of progressive collapse. In regular frame structures, beams and columns compete for the strengthening budget. In this paper, we present a risk-based formulation to address the optimal design of regular plane frames under element loss conditions. We address the threat probabilities for which strengthening has better cost-benefit than usual design, for different frame configurations, and study the impacts of strengthening extent and cost. The risk-based optimization reveals optimum points of compromise between competing failure modes: local bending of beams, local crushing of columns, and global pancake collapse, for frames of different aspect ratios. The conceptual study is based on a simple analytical model for progressive collapse, but it provides relevant insight for the design and strengthening of real structures.


title: Perspective on Code Submission and Automated Evaluation Platforms for University Teaching

abstract: We present a perspective on platforms for code submission and automated evaluation in the context of university teaching. Due to the COVID-19 pandemic, such platforms have become an essential asset for remote courses and a reasonable standard for structured code submission concerning increasing numbers of students in computer sciences. Utilizing automated code evaluation techniques exhibits notable positive impacts for both students and teachers in terms of quality and scalability. We identified relevant technical and non-technical requirements for such platforms in terms of practical applicability and secure code submission environments. Furthermore, a survey among students was conducted to obtain empirical data on general perception. We conclude that submission and automated evaluation involves continuous maintenance yet lowers the required workload for teachers and provides better evaluation transparency for students.


title: A Complete Helmholtz Decomposition on Multiply Connected Subdivision Surfaces and Its Application to Integral Equations

abstract: The analysis of electromagnetic scattering in the isogeometric analysis (IGA) framework based on Loop subdivision has long been restricted to simply-connected geometries. The inability to analyze multiply-connected objects is a glaring omission. In this paper, we address this challenge. IGA provides seamless integration between the geometry and analysis by using the same basis set to represent both. In particular, IGA methods using subdivision basis sets exploit the fact that the basis functions used for surface description are smooth (with continuous second derivatives) almost everywhere. On simply-connected surfaces, this permits the definition of basis sets that are divergence-free and curl-free. What is missing from this suite is a basis set that is both divergence-free and curl-free, a necessary ingredient for a complete Helmholtz decomposition of currents on multiply-connected structures. In this paper, we effect this missing ingredient numerically using random polynomial vector fields. We show that this basis set is analytically divergence-free and curl-free. Furthermore, we show that these basis recovers curl-free, divergence-free, and curl-free and divergence-free fields. Finally, we use this basis set to discretize a well-conditioned integral equation for analyzing perfectly conducting objects and demonstrate excellent agreement with other methods.


title: Evaluating a Methodology for Increasing AI Transparency: A Case Study

abstract: In reaction to growing concerns about the potential harms of artificial intelligence (AI), societies have begun to demand more transparency about how AI models and systems are created and used. To address these concerns, several efforts have proposed documentation templates containing questions to be answered by model developers. These templates provide a useful starting point, but no single template can cover the needs of diverse documentation consumers. It is possible in principle, however, to create a repeatable methodology to generate truly useful documentation. Richards et al. [25] proposed such a methodology for identifying specific documentation needs and creating templates to address those needs. Although this is a promising proposal, it has not been evaluated. This paper presents the first evaluation of this user-centered methodology in practice, reporting on the experiences of a team in the domain of AI for healthcare that adopted it to increase transparency for several AI models. The methodology was found to be usable by developers not trained in user-centered techniques, guiding them to creating a documentation template that addressed the specific needs of their consumers while still being reusable across different models and use cases. Analysis of the benefits and costs of this methodology are reviewed and suggestions for further improvement in both the methodology and supporting tools are summarized.


title: A note on the determinant 308 in Proskuryakov's linear algebra book

abstract: We put in evidence and correct a mistake in the formula for the determinant 308 in Proskuryakov's linear algebra book. We apply this formula to reprove the well-known fact that the Fubini-Study metric on the complex projective space is Einstein.


title: Online Assessment Misconduct Detection using Internet Protocol and Behavioural Classification

abstract: With the recent prevalence of remote education, academic assessments are often conducted online, leading to further concerns surrounding assessment misconducts. This paper investigates the potentials of online assessment misconduct (e-cheating) and proposes practical countermeasures against them. The mechanism for detecting the practices of online cheating is presented in the form of an e-cheating intelligent agent, comprising of an internet protocol (IP) detector and a behavioural monitor. The IP detector is an auxiliary detector which assigns randomised and unique assessment sets as an early procedure to reduce potential misconducts. The behavioural monitor scans for irregularities in assessment responses from the candidates, further reducing any misconduct attempts. This is highlighted through the proposal of the DenseLSTM using a deep learning approach. Additionally, a new PT Behavioural Database is presented and made publicly available. Experiments conducted on this dataset confirm the effectiveness of the DenseLSTM, resulting in classification accuracies of up to 90.7%.


title: A Proof Procedure For Separation Logic With Inductive Definitions and Theory Reasoning

abstract: A proof procedure, in the spirit of the sequent calculus, is proposed to check the validity of entailments between Separation Logic formulas combining inductively defined predicates denoted structures of bounded tree width and theory reasoning. The calculus is sound and complete, in the sense that a sequent is valid iff it admits a (possibly infinite) proof tree. We show that the procedure terminates in the two following cases: (i) When the inductive rules that define the predicates occurring on the left-hand side of the entailment terminate, in which case the proof tree is always finite. (ii) When the theory is empty, in which case every valid sequent admits a rational proof tree, where the total number of pairwise distinct sequents occurring in the proof tree is doubly exponential w.r.t.\ the size of the end-sequent.


title: Commuting Tuple of Multiplication Operators Homogeneous under the Unitary Group

abstract: Let $\mathcal U(d)$ be the group of $d\times d$ unitary matrices. We find conditions to ensure that a $\mathcal U(d)$-homogeneous $d$-tuple $\boldsymbol T$ is unitarily equivalent to multiplication by the coordinate functions on some reproducing kernel Hilbert space $\mathcal H_K(\mathbb B_d, \mathbb C^n) \subseteq \mbox{\rm Hol}(\mathbb B_d, \mathbb C^n)$, $n= \dim \cap_{j=1}^d \ker T^*_{j}.$ We describe this class of $\mathcal U(d)$-homogeneous operators, equivalently, non-negative kernels $K$ quasi-invariant under the action of $\mathcal U(d)$. We classify quasi-invariant kernels $K$ transforming under $\mathcal U(d)$ with two specific choice of multipliers. A crucial ingredient of the proof is that the group $SU(d)$ has exactly two inequivalent irreducible unitary representations of dimension $d$ and none in dimensions $2, \ldots , d-1$, $d\geq 3$. We obtain explicit criterion for boundedness, reducibility and mutual unitary equivalence among these operators.


title: Network-level Safety Metrics for Overall Traffic Safety Assessment: A Case Study

abstract: Driving safety analysis has recently experienced unprecedented improvements thanks to technological advances in precise positioning sensors, artificial intelligence (AI)-based safety features, autonomous driving systems, connected vehicles, high-throughput computing, and edge computing servers. Particularly, deep learning (DL) methods empowered volume video processing to extract safety-related features from massive videos captured by roadside units (RSU). Safety metrics are commonly used measures to investigate crashes and near-conflict events. However, these metrics provide limited insight into the overall network-level traffic management. On the other hand, some safety assessment efforts are devoted to processing crash reports and identifying spatial and temporal patterns of crashes that correlate with road geometry, traffic volume, and weather conditions. This approach relies merely on crash reports and ignores the rich information of traffic videos that can help identify the role of safety violations in crashes. To bridge these two perspectives, we define a new set of network-level safety metrics (NSM) to assess the overall safety profile of traffic flow by processing imagery taken by RSU cameras. Our analysis suggests that NSMs show significant statistical associations with crash rates. This approach is different than simply generalizing the results of individual crash analyses, since all vehicles contribute to calculating NSMs, not only the ones involved in crash incidents. This perspective considers the traffic flow as a complex dynamic system where actions of some nodes can propagate through the network and influence the crash risk for other nodes. We also provide a comprehensive review of surrogate safety metrics (SSM) in the Appendix A.


title: POTATO: exPlainable infOrmation exTrAcTion framewOrk

abstract: We present POTATO, a task- and languageindependent framework for human-in-the-loop (HITL) learning of rule-based text classifiers using graph-based features. POTATO handles any type of directed graph and supports parsing text into Abstract Meaning Representations (AMR), Universal Dependencies (UD), and 4lang semantic graphs. A streamlit-based user interface allows users to build rule systems from graph patterns, provides real-time evaluation based on ground truth data, and suggests rules by ranking graph features using interpretable machine learning models. Users can also provide patterns over graphs using regular expressions, and POTATO can recommend refinements of such rules. POTATO is applied in projects across domains and languages, including classification tasks on German legal text and English social media data. All components of our system are written in Python, can be installed via pip, and are released under an MIT License on GitHub.


title: Solving particle-antiparticle and cosmological constant problems

abstract: We argue that fundamental objects in particle theory are not elementary particles and antiparticles but objects described by irreducible representations (IRs) of the de Sitter (dS) algebra. One might ask why, then, experimental data give the impression that particles and antiparticles are fundamental and there are conserved additive quantum numbers (electric charge, baryon quantum number and others). The matter is that, at the present stage of the universe, the contraction parameter $R$ from the dS to the Poincare algebra is very large and, in the formal limit $R\to\infty$, one IR of the dS algebra splits into two IRs of the Poincare algebra corresponding to a particle and its antiparticle with the same masses. The problem why the quantities $(c,\hbar,R)$ are as are does not arise because they are contraction parameters for transitions from more general Lie algebras to less general ones. Then the baryon asymmetry of the universe problem does not arise. At the present stage of the universe, the phenomenon of cosmological acceleration (PCA) is described without uncertainties as an inevitable {\it kinematical} consequence of quantum theory in semiclassical approximation. In particular, it is not necessary to involve dark energy the physical meaning of which is a mystery. In our approach, background space and its geometry are not used and $R$ has nothing to do with the radius of dS space. In semiclassical approximation, the results for PCA are the same as in General Relativity if $\Lambda=3/R^2$, i.e., $\Lambda>0$ and there is no freedom in choosing the value of $\Lambda$.


title: EXSeQETIC: Expert System to Support the Implementation of eQETIC Model

abstract: The digital educational solutions are increasingly used demanding high quality functionalities. In this sense, standards and models are made available by governments, associations, and researchers being most used in quality control and assessment sessions. The eQETIC model was built according to the approach of continuous process improvement favoring the quality management for development and maintenance of digital educational solutions. This article presents two expert systems to support the implementation of eQETIC model and demonstrates that such systems are able to support users during the model implementation. Developed according to two types of shells (SINTA/UFC and e2gLite/eXpertise2go), the systems were used by a professional who develops these type of solutions and showed positive results regarding the support offered by them in implementing the rules proposed by eQETIC model.


title: "It's A Blessing and A Curse": Unpacking Creators' Practices with Non-Fungible Tokens (NFTs) and Their Communities

abstract: NFTs (Non-Fungible Tokens) are blockchain-based cryptographic tokens to represent ownership of unique content such as images, videos, or 3D objects. Despite NFTs' increasing popularity and skyrocketing trading prices, little is known about people's perceptions of and experiences with NFTs. In this work, we focus on NFT creators and present results of an exploratory qualitative study in which we interviewed 15 NFT creators from nine different countries. Our participants had nuanced feelings about NFTs and their communities. We found that most of our participants were enthusiastic about the underlying technologies and how they empower individuals to express their creativity and pursue new business models of content creation. Our participants also gave kudos to the NFT communities that have supported them to learn, collaborate, and grow in their NFT endeavors. However, these positivities were juxtaposed by their accounts of the many challenges that they encountered and thorny issues that the NFT ecosystem is grappling with around ownership of digital content, low-quality NFTs, scams, possible money laundering, and regulations. We discuss how the built-in properties (e.g., decentralization) of blockchains and NFTs might have contributed to some of these issues. We present design implications on how to improve the NFT ecosystem (e.g., making NFTs even more accessible to newcomers and the broader population).


title: A simple and fast algorithm for computing discrete Voronoi, Johnson-Mehl or Laguerre diagrams of points

abstract: This article presents an algorithm to compute digital images of Voronoi, Johnson-Mehl or Laguerre diagrams of a set of punctual sites, in a domain of a Euclidean space of any dimension. The principle of the algorithm is, in a first step, to investigate the voxels in balls centred around the sites, and, in a second step, to process the voxels remaining outside the balls. The optimal choice of ball radii can be determined analytically or numerically, which allows a performance of the algorithm in $O(N_v \ln N_s$), where $N_v$ is the total number of voxels of the domain and $N_s$ the number of sites of the tessellation. Periodic and non-periodic boundary conditions are considered. A major advantage of the algorithm is its simplicity which makes it very easy to implement. This makes the algorithm suitable for creating high resolution images of microstructures containing a large number of cells, in particular when calculations using FFT-based homogenisation methods are then to be applied to the simulated materials.


title: A hybrid deep learning approach for purchasing strategy of carbon emission rights -- Based on Shanghai pilot market

abstract: The price of carbon emission rights play a crucial role in carbon trading markets. Therefore, accurate prediction of the price is critical. Taking the Shanghai pilot market as an example, this paper attempted to design a carbon emission purchasing strategy for enterprises, and establish a carbon emission price prediction model to help them reduce the purchasing cost. To make predictions more precise, we built a hybrid deep learning model by embedding Generalized Autoregressive Conditional Heteroskedastic (GARCH) into the Gate Recurrent Unit (GRU) model, and compared the performance with those of other models. Then, based on the Iceberg Order Theory and the predicted price, we proposed the purchasing strategy of carbon emission rights. As a result, the prediction errors of the GARCH-GRU model with a 5-day sliding time window were the minimum values of all six models. And in the simulation, the purchasing strategy based on the GARCH-GRU model was executed with the least cost as well. The carbon emission purchasing strategy constructed by the hybrid deep learning method can accurately send out timing signals, and help enterprises reduce the purchasing cost of carbon emission permits.


title: Mass as a dynamical quantity

abstract: The Standard Model (SM) ascribes the observed mass of elementary particles to an effective interaction between basis states defined without mass terms and a scalar potential associated with the Higgs boson. In the relativistic field theory that underlies the SM, mass itself, understood as the Lorentz-invariant squared 4-momentum of a particle or field, is fixed a priori, imposing a constraint on possible momentum states. Stueckelberg introduced an alternative approach, positing antiparticles as particles evolving backward in time, thus relaxing the mass shell constraint for individual particles. Further work by Piron and Horwitz established a covariant Hamiltonian mechanics on an unconstrained 8D phase space, leading to a gauge field theory that mediates the exchange of mass between particles, while the total mass of particles and fields remains conserved. In a recently developed extension of general relativity, consistent with this approach, the spacetime metric evolves in a manner that permits the exchange of mass across spacetime through the gravitational field. Mechanisms that restrict mass exchange between particles have also been identified. Nevertheless, mass exchange remains possible under certain circumstances and may have phenomenological implications in particle physics and cosmology.


title: CDNNs: The coupled deep neural networks for coupling of the Stokes and Darcy-Forchheimer problems

abstract: In this article, we present an efficient deep learning method called coupled deep neural networks (CDNNs) for coupled physical problems. Our method compiles the interface conditions of the coupled PDEs into the networks properly and can be served as an efficient alternative to the complex coupled problems. To impose energy conservation constraints, the CDNNs utilize simple fully connected layers and a custom loss function to perform the model training process as well as the physical property of the exact solution. The approach can be beneficial for the following reasons: Firstly, we sampled randomly and only input spatial coordinates without being restricted by the nature of samples. Secondly, our method is meshfree which makes it more efficient than the traditional methods. Finally, our method is parallel and can solve multiple variables independently at the same time. We give the theory to guarantee the convergence of the loss function and the convergence of the neural networks to the exact solution. Some numerical experiments are performed and discussed to demonstrate the performance of the proposed method.


title: Fourth Generation fermions providing new candidates for Dark Matter and Dark Energy

abstract: Clifford Unification describes all known elementary fermions and their interactions, providing a useful tool in the analysis of experimental data. It also predicts the existence of a stable fourth generation (G4) of fermions, with electric charges different to their first, second and third generation counterparts. This makes neutral G4 baryon and neutral G4 lepton composites possible, respectively providing candidates for Dark Matter and Dark Energy, which are examined in the light of experimental evidence.


abstract: Quantum mechanics with quaternionic mass is considered. The momentum eigen-value equation with quaternionic mass yields the Klein-Gordon equation with a mass consisting of longitudinal and traverse masses. The scalar field total mass is found to be a sum of these masses. This field appears to be connected with two subfields conserving linear momentum. It is found that a particle with real mass satisfies the quantum Telegraph equation, whereas that one with quaternionic mass satisfies the Klein-Gordon equation. A quantum force acting on the particle is found to be proportional to its velocity. When the particle field is coupled to an electromagnetic field, an additional term in the particle's energy appears reflecting the interaction of the particle's angular momentum with the magnetic field.


title: Grid-Free Monte Carlo for PDEs with Spatially Varying Coefficients

abstract: Partial differential equations (PDEs) with spatially-varying coefficients arise throughout science and engineering, modeling rich heterogeneous material behavior. Yet conventional PDE solvers struggle with the immense complexity found in nature, since they must first discretize the problem -- leading to spatial aliasing, and global meshing/sampling that is costly and error-prone. We describe a method that approximates neither the domain geometry, the problem data, nor the solution space, providing the exact solution (in expectation) even for problems with extremely detailed geometry and intricate coefficients. Our main contribution is to extend the walk on spheres (WoS) algorithm from constant- to variable-coefficient problems, by drawing on techniques from volumetric rendering. In particular, an approach inspired by null-scattering yields unbiased Monte Carlo estimators for a large class of 2nd-order elliptic PDEs, which share many attractive features with Monte Carlo rendering: no meshing, trivial parallelism, and the ability to evaluate the solution at any point without solving a global system of equations.


title: Observation of magnetization surface textures of the van der Waals antiferromagnet FePS3 by spin Hall magnetoresistance

abstract: Van der Waals materials are a new platform to study two-dimensional systems, including magnetic order. Since the number of spins is relatively small, measuring the magnetization is challenging. Here we report spin Hall magnetoresistance (SMR) up to room temperature caused by the magnetic surface texture of exfoliated flakes of magnetic van der Waals materials. For the antiferromagnet FePS3 the SMR amounts to 0.1 % for an applied magnetic field of 7 T at 5 K which implies a substantial canting of the magnetic moments relative to the colinear antiferromagnetic order. The canting is substantial even for a magnetic field along the N\'eel vector, which illustrates the unique power of the SMR to detect magnetic surface textures in van der Waals magnets.


title: Correcting diacritics and typos with a ByT5 transformer model

abstract: Due to the fast pace of life and online communications and the prevalence of English and the QWERTY keyboard, people tend to forgo using diacritics, make typographical errors (typos) when typing in other languages. Restoring diacritics and correcting spelling is important for proper language use and the disambiguation of texts for both humans and downstream algorithms. However, both of these problems are typically addressed separately: the state-of-the-art diacritics restoration methods do not tolerate other typos, but classical spellcheckers also cannot deal adequately with all the diacritics missing. In this work, we tackle both problems at once by employing the newly-developed universal ByT5 byte-level seq2seq transformer model that requires no language-specific model structures. For a comparison, we perform diacritics restoration on benchmark datasets of 12 languages, with the addition of Lithuanian. The experimental investigation proves that our approach is able to achieve results (> 98%) comparable to the previous state-of-the-art, despite being trained less and on fewer data. Our approach is also able to restore diacritics in words not seen during training with > 76% accuracy. Our simultaneous diacritics restoration and typos correction approach reaches > 94% alpha-word accuracy on the 13 languages. It has no direct competitors and strongly outperforms classical spell-checking or dictionary-based approaches. We also demonstrate all the accuracies to further improve with more training. Taken together, this shows the great real-world application potential of our suggested methods to more data, languages, and error classes.


title: FastIoT -- A framework and holistic approach for rapid development of IIoT systems

abstract: While lots of research has been conducted on the architecture of Industrial Internet of Things (IIoT) systems, concepts of structuring their development processes are missing. Therefore, we propose a holistic approach supporting organizations in rapid development of IIoT systems. It includes the structuring of the development process into multiple projects sharing project conventions. Utilizing a single configurable build script for all projects, our goal is to make the integration of code from various projects into broader IIoT systems easy and the systems decomposable with minimal effort.


title: A generalization of a question asked by B. H. Neumann

abstract: Let $w \in F_2$ be a word and let $m$ and $n$ be two positive integers. We say that a finite group $G$ has the $w_{m,n}$-property if however a set $M$ of $m$ elements and a set $N$ of $n$ elements of the group is chosen, there exist at least one element of $x \in M$ and at least one element of $y \in M$ such that $w(x,y)=1.$ Assume that there exists a constant $\gamma < 1$ such that whenever $w$ is not an identity in a finite group $X$, then the probability that $w(x_1,x_2)=1$ in $X$ is at most $\gamma.$ If $m\leq n$ and $G$ satisfies the $w_{m,n}$-property, then either $w$ is an identity in $G$ or $|G|$ is bounded in terms of $\gamma, m$ and $n$. We apply this result to the 2-Engel word.


title: CP Violation for the Heavens and the Earth

abstract: Electroweak baryogenesis can be driven by the top quark in a general two Higgs doublet model with extra Yukawa couplings. Higgs quartics provide the first order phase transition, while extra top Yukawa coupling $\rho_{tt}$ can fuel the cosmic baryon asymmetry through the $\lambda_t\,{\rm Im}\,\rho_{tt}$ product, with flavor-changing $\rho_{tc}$ coupling as backup. The impressive ACME 2018 bound on the electron electric dipole moment calls for an extra electron coupling $\rho_{ee}$ for exquisite cancellation among dangerous diagrams, broadening the baryogenesis solution space. The mechanism suggests that extra Yukawa couplings echo the hierarchical structure of standard Yukawa couplings. Phenomenological consequences in Higgs search and flavor physics are discussed, with $\mu$ and $\tau$ EDM touched upon.


title: Impact of Naturalistic Field Acoustic Environments on Forensic Text-independent Speaker Verification System

abstract: Audio analysis for forensic speaker verification offers unique challenges in system performance due in part to data collected in naturalistic field acoustic environments where location/scenario uncertainty is common in the forensic data collection process. Forensic speech data as potential evidence can be obtained in random naturalistic environments resulting in variable data quality. Speech samples may include variability due to vocal efforts such as yelling over 911 emergency calls, whereas others might be whisper or situational stressed voice in a field location or interview room. Such speech variability consists of intrinsic and extrinsic characteristics and makes forensic speaker verification a complicated and daunting task. Extrinsic properties include recording equipment such as microphone type and placement, ambient noise, room configuration including reverberation, and other environmental scenario-based issues. Some factors, such as noise and non-target speech, will impact the verification system performance by their mere presence. To investigate the impact of field acoustic environments, we performed a speaker verification study based on the CRSS-Forensic corpus with audio collected from 8 field locations including police interviews. This investigation includes an analysis of the impact of seven unseen acoustic environments on speaker verification system performance using an x-Vector system.


title: The Sturdiness of the Shell Model: Informal Review

abstract: With shell model codes being able to encompassing larger and larger spaces we find that the percentage occupancy of the leading spaces becomes smaller and smaller. How can the shell model survive in such circumstances? We will not solve this puzzle here but rather will show examples where, with some explanations, the shell model holds fast. We will use nuclear moments as an example.


title: SafeAPT: Safe Simulation-to-Real Robot Learning using Diverse Policies Learned in Simulation

abstract: The framework of Simulation-to-real learning, i.e, learning policies in simulation and transferring those policies to the real world is one of the most promising approaches towards data-efficient learning in robotics. However, due to the inevitable reality gap between the simulation and the real world, a policy learned in the simulation may not always generate a safe behaviour on the real robot. As a result, during adaptation of the policy in the real world, the robot may damage itself or cause harm to its surroundings. In this work, we introduce a novel learning algorithm called SafeAPT that leverages a diverse repertoire of policies evolved in the simulation and transfers the most promising safe policy to the real robot through episodic interaction. To achieve this, SafeAPT iteratively learns a probabilistic reward model as well as a safety model using real-world observations combined with simulated experiences as priors. Then, it performs Bayesian optimization on the repertoire with the reward model while maintaining the specified safety constraint using the safety model. SafeAPT allows a robot to adapt to a wide range of goals safely with the same repertoire of policies evolved in the simulation. We compare SafeAPT with several baselines, both in simulated and real robotic experiments and show that SafeAPT finds high-performance policies within a few minutes in the real world while minimizing safety violations during the interactions.


title: Measured distribution of cloud chamber tracks from radioactive decay: a new empirical approach to investigating the quantum measurement problem

abstract: Using publicly available video of a diffusion cloud chamber with a very small radioactive source, I measure the spatial distribution of where tracks start, and consider possible implications. This is directly relevant to the quantum measurement problem and its possible resolution, and appears never to have been done before. The raw data are relatively uncontrolled, leading to caveats that should guide future, more tailored experiments. Aspects of the results may suggest a modification to Born's rule at very small wavefunction, with possibly profound implications for the detection of extremely rare events such as proton decay, but other explanations are not ruled out. Speculatively, I introduce two candidate small-wavefunction Born rule modifications, a hard cutoff, and an offset model with a stronger underlying physical rationale. Track distributions from decays in cloud chambers represent a previously unappreciated way to probe the foundations of quantum mechanics, and a novel case of wavefunctions with macroscopic signatures.


title: Differentiating and Integrating ZX Diagrams with Applications to Quantum Machine Learning

abstract: ZX-calculus has proved to be a useful tool for quantum technology with a wide range of successful applications. Most of these applications are of an algebraic nature. However, other tasks that involve differentiation and integration remain unreachable with current ZX techniques. Here we elevate ZX to an analytical perspective by realising differentiation and integration entirely within the framework of ZX-calculus. We explicitly illustrate the new analytic framework of ZX-calculus by applying it in context of quantum machine learning for the analysis of barren plateaus.


title: Stability conditions on fibred threefolds

abstract: We give a conjectural construction of Bridgeland stability conditions on the derived category of fibred threefolds. The construction depends on a conjectural Bogomolov-Gieseker type inequality for certain stable complexes. It can be considered as a relative version of the construction of Bayer, Macr\`i and Toda. We prove the conjectural Bogomolov-Gieseker type inequality in the case of relative projective planes over curves. This gives the the existence of Bridgeland stability conditions on such threefolds.


title: Homological properties of 0-Hecke modules for dual immaculate quasisymmetric functions

abstract: Let $n$ be a nonnegative integer. For each composition $\alpha$ of $n$, Berg $\textit{et al.}$ introduced a cyclic indecomposable $H_n(0)$-module $\mathcal{V}_\alpha$ with a dual immaculate quasisymmetric function as the image of the quasisymmetric characteristic. In this paper, we study $\mathcal{V}_\alpha$'s from the homological viewpoint. To be precise, we construct a minimal projective presentation of $\mathcal{V}_\alpha$ and a minimal injective presentation of $\mathcal{V}_\alpha$ as well. Using them, we compute ${\rm Ext}^1_{H_n(0)}(\mathcal{V}_\alpha, {\bf F}_\beta)$ and ${\rm Ext}^1_{H_n(0)}( {\bf F}_\beta, \mathcal{V}_\alpha)$, where ${\bf F}_\beta$ is the simple $H_n(0)$-module attached to a composition $\beta$ of $n$. We also compute ${\rm Ext}_{H_n(0)}^i(\mathcal{V}_\alpha,\mathcal{V}_{\beta})$ when $i=0,1$ and $\beta \le_l \alpha$, where $\le_l$ represents the lexicographic order on compositions.


title: On Two Families of Generalizations of Pascal's Triangle

abstract: We consider two families of Pascal-like triangles that have all ones on the left side and ones separated by $m-1$ zeros on the right side. The $m=1$ cases are Pascal's triangle and the two families also coincide when $m=2$. Members of the first family obey Pascal's recurrence everywhere inside the triangle. We show that the $m$-th triangle can also be obtained by reversing the elements up to and including the main diagonal in each row of the $(1/(1-x^m),x/(1-x))$ Riordan array. Properties of this family of triangles can be obtained quickly as a result. The $(n,k)$-th entry in the $m$-th member of the second family of triangles is the number of tilings of an $(n+k)\times1$ board that use $k$ $(1,m-1)$-fences and $n-k$ unit squares. A $(1,g)$-fence is composed of two unit square sub-tiles separated by a gap of width $g$. We show that the entries in the antidiagonals of these triangles are coefficients of products of powers of two consecutive Fibonacci polynomials and give a bijective proof that these coefficients give the number of $k$-subsets of $\{1,2,\ldots,n-m\}$ such that no two elements of a subset differ by $m$. Other properties of the second family of triangles are also obtained via a combinatorial approach. Finally, we give necessary and sufficient conditions for any Pascal-like triangle (or its row-reversed version) derived from tiling $(n\times1)$-boards to be a Riordan array.


title: Learning Hamiltonians of constrained mechanical systems

abstract: Recently, there has been an increasing interest in modelling and computation of physical systems with neural networks. Hamiltonian systems are an elegant and compact formalism in classical mechanics, where the dynamics is fully determined by one scalar function, the Hamiltonian. The solution trajectories are often constrained to evolve on a submanifold of a linear vector space. In this work, we propose new approaches for the accurate approximation of the Hamiltonian function of constrained mechanical systems given sample data information of their solutions. We focus on the importance of the preservation of the constraints in the learning strategy by using both explicit Lie group integrators and other classical schemes.


title: Multidimensional examples of the Metropolis algorithm

abstract: Consider the problem of approximating a given probability distribution on the cube $[0,1]^n$ via the use of a square lattice discretization with mesh-size $1/N$ and the Metropolis algorithm. Here the dimension $n$ is fixed and we focus for the most part on the case $n=2$. In order to understand the speed of convergence of such a procedure, one needs to control the spectral gap, $\lambda$, of the associated finite Markov chain, and how it depends on the parameter $N$. In this work, we study basic examples for which good upper-bounds and lower-bounds on $\lambda$ can be obtained via appropriate application of path techniques.


title: Proximal Denoiser for Convergent Plug-and-Play Optimization with Nonconvex Regularization

abstract: Plug-and-Play (PnP) methods solve ill-posed inverse problems through iterative proximal algorithms by replacing a proximal operator by a denoising operation. When applied with deep neural network denoisers, these methods have shown state-of-the-art visual performance for image restoration problems. However, their theoretical convergence analysis is still incomplete. Most of the existing convergence results consider nonexpansive denoisers, which is non-realistic, or limit their analysis to strongly convex data-fidelity terms in the inverse problem to solve. Recently, it was proposed to train the denoiser as a gradient descent step on a functional parameterized by a deep neural network. Using such a denoiser guarantees the convergence of the PnP version of the Half-Quadratic-Splitting (PnP-HQS) iterative algorithm. In this paper, we show that this gradient denoiser can actually correspond to the proximal operator of another scalar function. Given this new result, we exploit the convergence theory of proximal algorithms in the nonconvex setting to obtain convergence results for PnP-PGD (Proximal Gradient Descent) and PnP-ADMM (Alternating Direction Method of Multipliers). When built on top of a smooth gradient denoiser, we show that PnP-PGD and PnP-ADMM are convergent and target stationary points of an explicit functional. These convergence results are confirmed with numerical experiments on deblurring, super-resolution and inpainting.


title: A skew logistic distribution with application to modelling COVID-19 epidemic waves

abstract: A novel yet simple extension of the symmetric logistic distribution is proposed by introducing a skewness parameter. It is shown how the three parameters of the ensuing skew logistic distribution may be estimated using maximum likelihood. The skew logistic distribution is then extended to the skew bi-logistic distribution to allow the modelling of multiple waves in epidemic time series data. The proposed skew-logistic model is validated on COVID-19 data from the UK, and is evaluated for goodness-of-fit against the logistic and normal distributions using the recently formulated empirical survival Jensen-Shannon divergence (${\cal E}SJS$) and the Kolmogorov-Smirnov two-sample test statistic ($KS2$). We employ 95\% bootstrap confidence intervals to assess the improvement in goodness-of-fit of the skew logistic distribution over the other distributions. The obtained confidence intervals for the ${\cal E}SJS$ are narrower than those for the $KS2$ on using this data set, implying that the ${\cal E}SJS$ is more powerful than the $KS2$.


title: Liouville closed $H_T$-fields

abstract: Let $T$ be an o-minimal theory extending the theory of real closed ordered fields. An $H_T$-field is a model $K$ of $T$ equipped with a $T$-derivation such that the underlying ordered differential field of $K$ is an $H$-field. We study $H_T$-fields and their extensions. Our main result is that if $T$ is power bounded, then every $H_T$-field $K$ has either exactly one or exactly two minimal Liouville closed $H_T$-field extensions up to $K$-isomorphism. The assumption of power boundedness can be relaxed to allow for certain exponential cases, such as $T = \operatorname{Th}(\mathbb{R}_{\operatorname{an},\exp})$.


title: Trajectory balance: Improved credit assignment in GFlowNets

abstract: Generative flow networks (GFlowNets) are a method for learning a stochastic policy for generating compositional objects, such as graphs or strings, from a given unnormalized density by sequences of actions, where many possible action sequences may lead to the same object. We find previously proposed learning objectives for GFlowNets, flow matching and detailed balance, which are analogous to temporal difference learning, to be prone to inefficient credit propagation across long action sequences. We thus propose a new learning objective for GFlowNets, trajectory balance, as a more efficient alternative to previously used objectives. We prove that any global minimizer of the trajectory balance objective can define a policy that samples exactly from the target distribution. In experiments on four distinct domains, we empirically demonstrate the benefits of the trajectory balance objective for GFlowNet convergence, diversity of generated samples, and robustness to long action sequences and large action spaces.


title: Dynamical evolution of the Schottky barrier as a determinant contribution to electron-hole pair stabilization and photocatalysis of plasmon-induced hot carriers

abstract: The harnessing of plasmon-induced hot carriers promises to open new avenues for the development of clean energies and chemical catalysis. The extraction of carriers before thermalization and recombination is of primordial importance to obtain appealing conversion yields. Here, hot carrier injection in the paradigmatic Au-TiO$_{2}$ system is studied by means of electronic and electron-ion dynamics. Our results show that pure electronic features (without considering many-body interactions or dissipation to the environment) contribute to the electron-hole separation stability. These results reveal the existence of a dynamic contribution to the interfacial potential barrier (Schottky barrier) that arises at the charge injection pace, impeding electronic back transfer. Furthermore, we show that this charge separation stabilization provides the time needed for the charge to leak to capping molecules placed over the TiO$_{2}$ surface triggering a coherent bond oscillation that will lead to a photocatalytic dissociation. We expect that our results will add new perspectives to the interpretation of the already detected long-lived hot carrier lifetimes, their catalytical effect, and concomitantly to their technological applications.


title: A note on the use of equidistant contours for presenting scientific data

abstract: The passionate plea for the use of scientific colour maps misses some aspects in the visual presentation of scientific data. While a linear colour map based on scientific human colour perception is useful for the presentation of some images, like the three examples given of the topography of the earth, an apple and a passport photograph, scientific data are not presented. In this note it will be shown that there is more in scientific oceanographic data as they are presented in forms varying from historic equidistant contours, via a linear black-(grey)-white b&w map, a linear colour map and a nonlinear colour map. From an objective perspective, equidistant contouring is the best means for presenting scientific information in a relatively unbiased way. Nonlinear colour maps may add information to that by highlighting certain aspects also by varying the colour range if needed. Such information is not available from linear colour maps. Finally, the aesthetic aspect of visual data presentation is discussed.


title: A geometric model in 3+1D space-time for electrodynamic phenomena

abstract: With the idea to find geometric formulations of particle physics we investigate the predictions of a three dimensional generalisation of the Sine-Gordon model, very close to the Skyrme model and to the Wu-Yang description of Dirac monopoles. With three rotational degrees of freedom of spatial Dreibeins we formulate a Lagrangian and confront the predictions to electromagnetic phenomena. Stable solitonic excitations we compare with the lightest fundamental electric charges, electrons and positrons. Two Goldstone bosons we relate to the properties of photons. These particles are characterised by three topological quantum numbers, which we compare to charge, spin and photon number. Finally we conjecture some ideas for further comparisons with experiments.


title: Bootstrap percolation on the stochastic block model

abstract: We analyze the bootstrap percolation process on the stochastic block model (SBM), a natural extension of the Erd\H{o}s--R\'{e}nyi random graph that incorporates the community structure observed in many real systems. In the SBM, nodes are partitioned into two subsets, which represent different communities, and pairs of nodes are independently connected with a probability that depends on the communities they belong to. Under mild assumptions on the system parameters, we prove the existence of a sharp phase transition for the final number of active nodes and characterize the sub-critical and the super-critical regimes in terms of the number of initially active nodes, which are selected uniformly at random in each community.


title: X-ray diffraction with micrometer spatial resolution for highly absorbing samples

abstract: X-ray diffraction with high spatial resolution is commonly used to characterize (poly-)crystalline samples with, for example, respect to local strain, residual stress, grain boundaries and texture. However, the investigation of highly absorbing samples or the simultaneous assessment of high-Z materials by X-ray fluorescence have been limited due to the utilisation of low photon energies. Here, we report on a goniometer-based setup implemented at the P06 beamline of PETRA III that allows for micrometer spatial resolution with a photon energy of 35 keV and above. A highly focused beam was achieved by using compound refractive lenses and high precision sample manipulation was enabled by a goniometer that allows for up to 5D scans (3 rotations & 2 translations). As experimental examples, we demonstrate the determination of local strain variations in martensitic steel samples with micrometer spatial resolution as well as the simultaneous elemental distribution for high-Z materials in a thin film solar cell. Our proposed approach allows users from the materials science community to determine micro-structural properties even in highly absorbing samples.


title: Local existence of strong solutions to micro-macro models for reactive transport in evolving porous media

abstract: Two-scale models pose a promising approach in simulating reactive flow and transport in evolving porous media. Classically, homogenized flow and transport equations are solved on the macroscopic scale, while effective parameters are obtained from auxiliary cell problems on possibly evolving reference geometries (micro-scale). Despite their perspective success in rendering lab/field-scale simulations computationally feasible, analytic results regarding the arising two-scale bilaterally coupled system often restrict to simplified models. In this paper, we first derive smooth-dependence results concerning the partial coupling from the underlying geometry to macroscopic quantities. Therefore, alterations of the representative fluid domain are described by smooth paths of diffeomorphisms. Exploiting the gained regularity of the effective space- and time-dependent macroscopic coefficients, we present local-in-time existence results for strong solutions to the partially coupled micro-macro system using fixed-point arguments. What is more, we extend our results to the bilaterally coupled diffusive transport model including a level-set description of the evolving geometry.


title: Aggregation and Transformation of Vector-Valued Messages in the Shuffle Model of Differential Privacy

abstract: Advances in communications, storage and computational technology allow significant quantities of data to be collected and processed by distributed devices. Combining the information from these endpoints can realize significant societal benefit but presents challenges in protecting the privacy of individuals, especially important in an increasingly regulated world. Differential privacy (DP) is a technique that provides a rigorous and provable privacy guarantee for aggregation and release. The Shuffle Model for DP has been introduced to overcome challenges regarding the accuracy of local-DP algorithms and the privacy risks of central-DP. In this work we introduce a new protocol for vector aggregation in the context of the Shuffle Model. The aim of this paper is twofold; first, we provide a single message protocol for the summation of real vectors in the Shuffle Model, using advanced composition results. Secondly, we provide an improvement on the bound on the error achieved through using this protocol through the implementation of a Discrete Fourier Transform, thereby minimizing the initial error at the expense of the loss in accuracy through the transformation itself. This work will further the exploration of more sophisticated structures such as matrices and higher-dimensional tensors in this context, both of which are reliant on the functionality of the vector case.


title: Micro-level Reserving for General Insurance Claims using a Long Short-Term Memory Network

abstract: Detailed information about individual claims are completely ignored when insurance claims data are aggregated and structured in development triangles for loss reserving. In the hope of extracting predictive power from the individual claims characteristics, researchers have recently proposed to move away from these macro-level methods in favor of micro-level loss reserving approaches. We introduce a discrete-time individual reserving framework incorporating granular information in a deep learning approach named Long Short-Term Memory (LSTM) neural network. At each time period, the network has two tasks: first, classifying whether there is a payment or a recovery, and second, predicting the corresponding non-zero amount, if any. We illustrate the estimation procedure on a simulated and a real general insurance dataset. We compare our approach with the chain-ladder aggregate method using the predictive outstanding loss estimates and their actual values. Based on a generalized Pareto model for excess payments over a threshold, we adjust the LSTM reserve prediction to account for extreme payments.


title: PreDefense: Defending Underserved AI Students and Researchers from Predatory Conferences

abstract: Mentorship in the AI community is crucial to maintaining and increasing diversity, especially with respect to fostering the academic growth of underserved students. While the research process itself is important, there is not sufficient emphasis on the submission, presentation, and publication process, which is a cause for concern given the meteoric rise of predatory scientific conferences, which are based on profit only and have little to no peer review. These conferences are a direct threat to integrity in science by promoting work with little to no scientific merit. However, they also threaten diversity in the AI community by marginalizing underrepresented groups away from legitimate conferences due to convenience and targeting mechanisms like e-mail invitations. Due to the importance of conference presentation in AI research, this very specific problem must be addressed through direct mentorship. In this work, we propose PreDefense, a mentorship program that seeks to guide underrepresented students through the scientific conference and workshop process, with an emphasis on choosing legitimate venues that align with the specific work that the students are focused in and preparing students of all backgrounds for future successful, integrous AI research careers.


title: Towards machine learning aided real-time range imaging in proton therapy

abstract: In this work, we report on the advantageous aspects of the i-TED Compton imager for proton-range monitoring, based on the results of the first Monte Carlo study of its applicability to this field. i-TED is an array of Compton cameras, that have been designed for neutron-capture nuclear physics experiments, which are characterized by $\gamma$-ray energies spanning up to 5-6 MeV, rather low $\gamma$-ray emission yields and intense neutron induced $\gamma$-ray backgrounds. Our developments to cope with these three aspects are concomitant with those required in the field of hadron therapy, especially in terms of high efficiency for real-time monitoring, low sensitivity to neutron backgrounds and reliable performance at the high $\gamma$-ray energies. We find that signal-to-background ratios can be appreciably improved with i-TED thanks to its light-weight design and the low neutron-capture cross sections of its LaCl$_{3}$ crystals, when compared to other similar systems based on LYSO, CdZnTe or LaBr$_{3}$. Its high time-resolution (CRT$\sim$500 ps) represents an additional advantage for background suppression when operated in pulsed HT mode. Each i-TED module features two detection planes of very large LaCl$_{3}$ monolithic crystals, thereby achieving a high efficiency in coincidence of 0.2% for a point-like 1MeV $\gamma$-ray source at 5 cm distance. This leads to sufficient statistics for reliable image reconstruction with an array of four i-TED detectors assuming clinical intensities of 10$^{8}$ protons per treatment point. The use of a two-plane design instead of three-planes has been preferred owing to the higher attainable efficiency for double time-coincidences than for threefold events. The loss of full-energy events for high energy $\gamma$-rays is compensated by means of Machine-Learning algorithms, which allow one to enhance the signal-to-total ratio up to a factor of 2.


title: On Ternary Diophantine Equations of Signature $(p,p,3)$ over Number Fields

abstract: In this paper, we prove results about solutions of the Diophantine equation $x^p+y^p=z^3$ over various number fields using the modular method. Firstly, by assuming some standard modularity conjecture we prove an asymptotic result for general number fields of narrow class number one satisfying some technical conditions. Secondly, we show that there is an explicit bound such that the equation $x^p+y^p=z^3$ does not have a particular type of solution over $K=\Q(\sqrt{-d})$ where $d=1,7,19,43,67$ whenever $p$ is bigger than this bound. During the course of the proof we prove various results about the irreducibility of Galois representations, image of inertia groups and Bianchi newforms.


title: StRegA: Unsupervised Anomaly Detection in Brain MRIs using a Compact Context-encoding Variational Autoencoder

abstract: Expert interpretation of anatomical images of the human brain is the central part of neuro-radiology. Several machine learning-based techniques have been proposed to assist in the analysis process. However, the ML models typically need to be trained to perform a specific task, e.g., brain tumour segmentation or classification. Not only do the corresponding training data require laborious manual annotations, but a wide variety of abnormalities can be present in a human brain MRI - even more than one simultaneously, which renders representation of all possible anomalies very challenging. Hence, a possible solution is an unsupervised anomaly detection (UAD) system that can learn a data distribution from an unlabelled dataset of healthy subjects and then be applied to detect out of distribution samples. Such a technique can then be used to detect anomalies - lesions or abnormalities, for example, brain tumours, without explicitly training the model for that specific pathology. Several Variational Autoencoder (VAE) based techniques have been proposed in the past for this task. Even though they perform very well on controlled artificially simulated anomalies, many of them perform poorly while detecting anomalies in clinical data. This research proposes a compact version of the "context-encoding" VAE (ceVAE) model, combined with pre and post-processing steps, creating a UAD pipeline (StRegA), which is more robust on clinical data, and shows its applicability in detecting anomalies such as tumours in brain MRIs. The proposed pipeline achieved a Dice score of 0.642$\pm$0.101 while detecting tumours in T2w images of the BraTS dataset and 0.859$\pm$0.112 while detecting artificially induced anomalies, while the best performing baseline achieved 0.522$\pm$0.135 and 0.783$\pm$0.111, respectively.


title: Output-Feedback Control of Viscous Liquid-Tank System and its Numerical Approximation

abstract: We solve the output-feedback stabilization problem for a tank with a liquid modeled by the viscous Saint-Venant PDE system. The control input is the acceleration of the tank and a Control Lyapunov Functional methodology is used. The measurements are the tank position and the liquid level at the tank walls. The control scheme is a combination of a state feedback law with functional observers for the tank velocity and the liquid momentum. Four different types of output feedback stabilizers are proposed. A full-order observer and a reduced-order observer are used in order to estimate the tank velocity while the unmeasured liquid momentum is either estimated by using an appropriate scalar filter or is ignored. The reduced order observer differs from the full order observer because it omits the estimation of the measured tank position. Exponential convergence of the closed-loop system to the desired equilibrium point is achieved in each case. An algorithm is provided that guarantees that a robotic arm can move a glass of water to a pre-specified position no matter how full the glass is, without spilling water out of the glass, without residual end point sloshing and without measuring the water momentum and the glass velocity. Finally, the efficiency of the proposed output feedback laws is validated by numerical examples, obtained by using a simple finite-difference numerical scheme. The properties of the proposed, explicit, finite-difference scheme are determined.

categories: math.OC cs.NA cs.SY eess.SY math.AP math.NA physics.flu-dyn


title: Some asymptotic results for time series model selection

abstract: We consider the model selection problem for a large class of time series models, including, multivariate count processes, causal processes with exogenous covariates. A procedure based on a general penalized contrast is proposed. Some asymptotic results for weak and strong consistency are established. The non consistency issue is addressed, and a class of penalty term, that does not ensure consistency is provided. Examples of continuous valued and multivariate count autoregressive time series are considered.


title: TOI-1442 b and TOI-2445 b: two potentially rocky ultra-short period planets around M dwarfs

abstract: Context. Exoplanets with orbital periods of less than one day are known as ultra-short period (USP) planets. They are relatively rare products of planetary formation and evolution processes, but especially favourable for characterisation with current planet detection methods. At the time of writing, 125 USP planets have already been confirmed. Aims. Our aim is to validate the planetary nature of two new transiting planet candidates around M dwarfs announced by the NASA Transiting Exoplanet Survey Satellite (TESS), registered as TESS Objects of Interest (TOIs) TOI-1442.01 and TOI-2445.01. Methods. We used TESS data, ground-based photometric light curves, and Subaru/IRD spectrograph radial velocity (RV) measurements to validate both planetary candidates and to establish their physical properties. Results. TOI-1442 b is a validated exoplanet with an orbital period of $P$=0.4090682+/-0.0000004 d, a radius of $R_p$=1.15+/-0.06$R_{\oplus}$, and equilibrium temperature of $T_{p,eq}$=1357$_{-42}^{+49}$K. TOI-2445 b is also validated with an orbital period of $P$=0.3711286+/-0.0000004 d, a radius of $R_p$= 1.33+/-0.09$R_{\oplus}$, and equilibrium temperature of $T_{p,eq}$=1330$_{-56}^{+61}$K. Their physical properties align with current empirical trends and formation theories of USP planets. Based on the RV measurements, we set 3$\sigma$ upper mass limits of 8$M_{\oplus}$ and 20$M_{\oplus}$, thus confirming the non-stellar, sub-Jovian nature of both transiting objects. More RV measurements will be needed to constrain the planetary masses and mean densities, and the predicted presence of outer planetary companions. These targets extend the small sample of USP planets orbiting around M dwarfs up to 21 members. They are also among the 20 most suitable terrestrial planets for atmospheric characterisation via secondary eclipse with the JWST, according to a widespread emission spectroscopy metric.


title: Spectroscopy of the massive interacting binary UU Cassiopeiae

abstract: Context. The eclipsing close binary UU Cas is an interacting massive double-periodic system with a gainer star partly hidden in an accretion disk. Aims. In order to study the physics of the accretion process in greater detail, along with the structure and dynamics of the circumstellar matter in the system, we supplement our previous results obtained from photometry with an analysis of the spectra of UU Cas. Methods. We collected all available spectra used in previous publications on UU Cas and we acquired new ones. The method of disentangling was applied to this set of spectra spanning the years 2008-2021. The orbital parameters were disentangled and a fit of the separated component spectra by synthetic ones has been used to determine the physical parameters of the component stars. We compared the results to models of the evolution of interacting binaries. Results. We found that in addition to the dominant role of the donor star and a weak contribution of the gainer, the line profiles are strongly influenced by the circumstellar matter. The absorption lines reveal a presence of a disk wind emanating above the orbital plane. The variability of H{\alpha} emission yields evidence of changes in the structure of the circumstellar matter on a timescale of several orbital periods.


title: A Tunable Monolithic SQUID in Twisted Bilayer Graphene

abstract: Magic-angle twisted bilayer graphene (MATBG) hosts a number of correlated states of matter that can be tuned by electrostatic doping. Superconductivity has drawn considerable attention and the mechanism behind it is a topic of active discussion. MATBG has been experimentally characterized by numerous transport and scanning-probe experiments. The material has also emerged as a versatile platform for superconducting electronics, as proven by the realization of monolithic Josephson junctions. However, even though phase-coherent phenomena have been measured, no control of the superconducting phase has been demonstrated so far. Here, we present a Superconducting Quantum Interference Device (SQUID) in MATBG, where the superconducting phase difference is controlled through the magnetic field. We observe magneto-oscillations of the critical current, demonstrating long-range coherence agreeing with an effective charge of 2e for the superconducting charge carriers. We tune to both asymmetric and symmetric SQUID configurations by electrostatically controlling the critical currents through the junctions. With this tunability, we study the inductances in the device, finding values of up to 2{\mu}H. Furthermore, we directly observe the current-phase relation of one of the Josephson junctions of the device. Our results show that superconducting devices in MATBG can be scaled up and used to reveal properties of the material. We expect this to foster a more systematic realization of devices of this type, increasing the accuracy with which microscopic characteristics of the material are extracted. We also envision more complex devices to emerge, such as phase-slip junctions or high kinetic inductance detectors.


title: On the Hofer-Zehnder conjecture on weighted projective spaces

abstract: We prove an extension of the homology version of the Hofer-Zehnder conjecture proved by Shelukhin to the weighted projective spaces which are symplectic orbifolds. In particular, we prove that if the number of fixed points counted with their isotropy order as multiplicity of a non-degenerate Hamiltonian diffeomorphism of such a space is larger than the minimum number possible, then there are infinitely many periodic points.


title: Combining Local and Global Pose Estimation for Precise Tracking of Similar Objects

abstract: In this paper, we present a multi-object 6D detection and tracking pipeline for potentially similar and non-textured objects. The combination of a convolutional neural network for object classification and rough pose estimation with a local pose refinement and an automatic mismatch detection enables direct application in real-time AR scenarios. A new network architecture, trained solely with synthetic images, allows simultaneous pose estimation of multiple objects with reduced GPU memory consumption and enhanced performance. In addition, the pose estimates are further improved by a local edge-based refinement step that explicitly exploits known object geometry information. For continuous movements, the sole use of local refinement reduces pose mismatches due to geometric ambiguities or occlusions. We showcase the entire tracking pipeline and demonstrate the benefits of the combined approach. Experiments on a challenging set of non-textured similar objects demonstrate the enhanced quality compared to the baseline method. Finally, we illustrate how the system can be used in a real AR assistance application within the field of construction.


title: UQGAN: A Unified Model for Uncertainty Quantification of Deep Classifiers trained via Conditional GANs

abstract: We present an approach to quantifying both aleatoric and epistemic uncertainty for deep neural networks in image classification, based on generative adversarial networks (GANs). While most works in the literature that use GANs to generate out-of-distribution (OoD) examples only focus on the evaluation of OoD detection, we present a GAN based approach to learn a classifier that produces proper uncertainties for OoD examples as well as for false positives (FPs). Instead of shielding the entire in-distribution data with GAN generated OoD examples which is state-of-the-art, we shield each class separately with out-of-class examples generated by a conditional GAN and complement this with a one-vs-all image classifier. In our experiments, in particular on CIFAR10, CIFAR100 and Tiny ImageNet, we improve over the OoD detection and FP detection performance of state-of-the-art GAN-training based classifiers. Furthermore, we also find that the generated GAN examples do not significantly affect the calibration error of our classifier and result in a significant gain in model accuracy.


title: Hierarchical clustering of mixed-type data based on barycentric coding

abstract: Clustering of mixed-type datasets can be a particularly challenging task as it requires taking into account the associations between variables with different level of measurement, i.e., nominal, ordinal and/or interval. In some cases, hierarchical clustering is considered a suitable approach, as it makes few assumptions about the data and its solution can be easily visualized. Since most hierarchical clustering approaches assume variables are measured on the same scale, a simple strategy for clustering mixed-type data is to homogenize the variables before clustering. This would mean either recoding the continuous variables as categorical ones or vice versa. However, typical discretization of continuous variables implies loss of information. In this work, an agglomerative hierarchical clustering approach for mixed-type data is proposed, which relies on a barycentric coding of continuous variables. The proposed approach minimizes information loss and is compatible with the framework of correspondence analysis. The utility of the method is demonstrated on real and simulated data.


title: Convergence of level sets in fractional Laplacian regularization

abstract: The use of the fractional Laplacian in image denoising and regularization of inverse problems has enjoyed a recent surge in popularity, since for discontinuous functions it can behave less aggressively than methods based on $H^1$ norms, while being linear and computable with fast spectral numerical methods. In this work, we examine denoising and linear inverse problems regularized with fractional Laplacian in the vanishing noise and regularization parameter regime. The clean data is assumed piecewise constant in the first case, and continuous and satisfying a source condition in the second. In these settings, we prove results of convergence of level set boundaries with respect to Hausdorff distance, and additionally convergence rates in the case of denoising and indicatrix clean data. The main technical tool for this purpose is a family of barriers constructed by Savin and Valdinoci for studying the fractional Allen-Cahn equation. To help put these fractional methods in context, comparisons with the total variation and classical Laplacian are provided throughout.


title: On Tusi's Classification of Cubic Equations and its Connections to Cardano's Formula and Khayyam's Geometric Solution

abstract: Omar Khayyam's studies on cubic equations inspired the 12th century Persian mathematician Sharaf al-Din Tusi to investigate the number of positive roots. According to the noted mathematical historian Rashed, Tusi analyzed the problem for five different types of equations. In fact all cubic equations are reducible to a form {\it Tusi form} $x^2-x^3=c$. Tusi determined that the maximum of $x^2-x^3$ on $(0,1)$ occurs at $\frac{2}{3}$ and concluded when $c=\frac{4}{27} \delta$, $\delta \in (0,1)$, there are roots in $(0, \frac{2}{3})$ and $(\frac{2}{3},1)$, ignoring the root in $(-\frac{1}{3},0)$. Given a {\it reduced form} $x^3+px+q=0$, when $p <0$, we show it is reducible to a Tusi form with $\delta = \frac{1}{2} + {3\sqrt{3} q}/{4\sqrt{-p^3}}$. It follows there are three real roots if and only if $\Delta =-(\frac{q^2}{4}+\frac{p^3}{27})$ is positive. This gives an explicit connection between $\delta$ in Tusi form and $\Delta$ in Cardano's formula. Thus when $\delta \in (0,1)$, rather than using Cardano's formula in complex numbers one can approximate the roots iteratively. On the other hand, for a reduced form with $p >0$ we give a novel proof of Cardono's formula. While Rashed attributes Tusi's computation of the maximum to the use of derivatives, according to Hogendijk, Tusi was probably influenced by Euclid. Here we show the maximizer in Tusi form is computable via elementary algebraic manipulations. Indeed for a {\it quadratic Tusi form}, $x-x^2=\delta/4$, Tusi's approach results in a simple derivation of the quadratic formula, comparable with the pedagogical approach of Po-Shen Loh. Moreover, we derive analogous results for the {\it general Tusi form}. Finally, we present a novel derivation of Khayyam's geometric solution. The results complement previous findings on Tusi's work and reveal further facts on history, mathematics and pedagogy in solving cubic equations.


title: On reversible asynchronous non-uniform cellular automata

abstract: We study the class of asynchronous non-uniform cellular automata (ANUCA) over an arbitrary group universe with multiple local transition rules. We introduce the notion of stable injectivity, stable reversibility, stable post-surjectivity and investigate several dynamical properties of such automata. In particular, we establish the equivalence between reversibility, stable reversibility, and stable injectivity for ANUCA. We also prove the invertibility of several classes of injective and stably injective ANUCA. Counter-examples are given to highlight the differences between cellular automata and ANUCA.