title: Homogeneous Sets in Graphs and a Chromatic Multisymmetric Function

abstract: In this paper, we extend the chromatic symmetric function $X$ to a chromatic $k$-multisymmetric function $X_k$, defined for graphs equipped with a partition of their vertex set into $k$ parts. We demonstrate that this new function retains the basic properties and basis expansions of $X$, and we give a method for systematically deriving new linear relationships for $X$ from previous ones by passing them through $X_k$. In particular, we show how to take advantage of homogeneous sets of $G$ (those $S \subseteq V(G)$ such that each vertex of $V(G) \backslash S$ is either adjacent to all of $S$ or is nonadjacent to all of $S$) to relate the chromatic symmetric function of $G$ to those of simpler graphs. Furthermore, we show how extending this idea to homogeneous pairs $S_1 \sqcup S_2 \subseteq V(G)$ generalizes the process used by Guay-Paquet to reduce the Stanley-Stembridge conjecture to unit interval graphs.


title: On the $(6,4)$-problem of Brown, Erd\H{o}s and S\'os

abstract: Let $f^{(r)}(n;s,k)$ be the maximum number of edges of an $r$-uniform hypergraph on $n$ vertices not containing a subgraph with $k$ edges and at most $s$ vertices. In 1973, Brown, Erd\H{o}s and S\'os conjectured that the limit $$\lim_{n\to \infty} n^{-2} f^{(3)}(n;k+2,k)$$ exists for all $k$ and confirmed it for $k=2$. Recently, Glock showed this for $k=3$. We settle the next open case, $k=4$, by showing that $f^{(3)}(n;6,4)=\left(\frac{7}{36}+o(1)\right)n^2$ as $n\to\infty$. More generally, for all $k\in \{3,4\}$, $r\ge 3$ and $t\in [2,r-1]$, we compute the value of the limit $\lim_{n\to \infty} n^{-t}f^{(r)}(n;k(r-t)+t,k)$, which settles a problem of Shangguan and Tamo.


title: Solving multiorbital dynamical mean-field theory using natural orbitals renormalization group

abstract: The natural orbitals renormalization group (NORG) has previously been proposed as an efficient numerical method for solving zero-temperature properties of multisite and multiorbital quantum impurity systems. Here, we implement the NORG as an impurity solver for dynamical mean-field theory (DMFT). In comparison with the exact diagonalization method, the NORG method can treat much more bath sites in an impurity model to which the DMFT maps a lattice model and can find accurate zero-temperature Matsubara and low-frequency retarded Green's functions. We demonstrate the effectiveness of this method on a two-orbital Hubbard model on the Bethe lattice and find successfully the orbital selective Mott transition with a Kondo resonance peak in the wide band and two holon-doublon bound state excitation peaks in the narrow band.


title: Unidirectional magnetic coupling

abstract: We show that interlayer Dzyaloshinskii-Moriya interaction in combination with non-local Gilbert damping gives rise to unidirectional magnetic coupling. That is, the coupling between two magnetic layers -- say the left and right layer -- is such that dynamics of the left layer leads to dynamics of the right layer, but not vice versa. We discuss the implications of this result for the magnetic susceptibility of a magnetic bilayer, electrically-actuated spin-current transmission, and unidirectional spin-wave packet generation and propagation. Our results may enable a route towards spin-current and spin-wave diodes and further pave the way to design spintronic devices via reservoir engineering.


title: Spectral and Timing properties of the recently discovered Be/X-ray pulsar eRASSUJ 052914.9-662446

abstract: We have presented NuSTAR and Swift observations of the newly discovered Be/X-ray pulsar eRASSU J052914.9-662446. This is the first detailed study of the temporal and spectral properties of the pulsar using 2020 observations. A coherent pulsation of 1411.5$\pm$0.5 s was detected from the source. The pulse profile was found to resemble a simple single peaked feature which may be due to emission from the surface of the neutron star only. Pulse profiles are highly energy dependent. The variation of the pulse fraction of the pulse profiles are found to be non-monotonic with energy. The 0.5-20 keV Swift and NuSTAR simultaneous can be fitted well with power-law modified by high energy cutoff of $\sim$ 5.7 keV. The NuSTAR luminosity in the 0.5-79 keV energy range was $\sim$ 7.9x10 35 erg/s. The spectral flux in 3-79 keV shows modulation with the pulse phase.


title: Linear estimation of average global effects

abstract: We study the problem of estimating the average causal effect of treating every member of a population, as opposed to none, using an experiment that treats only some. This is the policy-relevant estimand when deciding whether to scale up an intervention based on the results of an RCT, for example, but differs from the usual average treatment effect in the presence of spillovers. We consider both estimation and experimental design given a bound (parametrized by $\eta > 0$) on the rate at which spillovers decay with the ``distance'' between units, defined in a generalized way to encompass spatial and quasi-spatial settings, e.g. where the economically relevant concept of distance is a gravity equation. Over all estimators linear in the outcomes and all cluster-randomized designs the optimal geometric rate of convergence is $n^{-\frac{1}{2+\frac{1}{\eta}}}$, and this rate can be achieved using a generalized ``Scaling Clusters'' design that we provide. We then introduce the additional assumption, implicit in the OLS estimators used in recent applied studies, that potential outcomes are linear in population treatment assignments. These estimators are inconsistent for our estimand, but a refined OLS estimator is consistent and rate optimal, and performs better than IPW estimators when clusters must be small. Its finite-sample performance can be improved by incorporating prior information about the structure of spillovers. As a robust alternative to the linear approach we also provide a method to select estimator-design pairs that minimize a notion of worst-case risk when the data generating process is unknown. Finally, we provide asymptotically valid inference methods.


title: A Hochschild-Kostant-Rosenberg theorem and residue sequences for logarithmic Hochschild homology

abstract: This paper incorporates the theory of Hochschild homology into our program on log motives. We discuss a geometric definition of logarithmic Hochschild homology of derived pre-log rings and construct an Andr\'e-Quillen type spectral sequence. The latter degenerates for derived log smooth maps between discrete pre-log rings. We employ this to show a logarithmic version of the Hochschild-Kostant-Rosenberg theorem and that logarithmic Hochschild homology is representable in the category of log motives. Among the applications, we deduce a generalized residue sequence involving blow-ups of log schemes.


title: Primordial Black Hole Formation in Starobinsky's Linear Potential Model

abstract: We study the power spectrum of the comoving curvature perturbation $\cal R$ in the model that glues two linear potentials of different slopes, originally proposed by Starobinsky. We find that the enhanced power spectrum reaches its maximum at the wavenumber which is $\pi$ times the junction scale. The peak is $\sim2.61$ times larger than the ultraviolet plateau. We also show that its near-peak behavior can be well approximated by a constant-roll model, once we define the effective ultra-slow-roll $e$-folding number appropriately by considering the contribution from non-single-clock phase only. Such an abrupt transition to non-attractor phase can leave some interesting characteristic features in the energy spectrum of the scalar-induced gravitational waves, which are detectable in the space-borne interferometers if the primordial black holes generated at such a high peak are all the dark matter.


title: Possible points of blow-up in chemotaxis systems with spatially heterogeneous logistic source

abstract: We discuss the influence of possible spatial inhomogeneities in the coefficients of logistic source terms in parabolic-elliptic chemotaxis-growth systems of the form \begin{align*} u_t &= \Delta u - \nabla\cdot(u\nabla v) + \kappa(x)u-\mu(x)u^2, 0 &= \Delta v - v + u \end{align*} in smoothly bounded domains $\Omega\subset\mathbb{R}^2$. Assuming that the coefficient functions satisfy $\kappa,\mu\in C^0(\overline{\Omega})$ with $\mu\geq0$ we prove that finite-time blow-up of the classical solution can only occur in points where $\mu$ is zero, i.e.\ that the blow-up set $\mathcal{B}$ is contained in \begin{align*} \big\{x\in\overline{\Omega}\mid\mu(x)=0\big\}. \end{align*} Moreover, we show that whenever $\mu(x_0)>0$ for some $x_0\in\overline{\Omega}$, then one can find an open neighbourhood $U$ of $x_0$ in $\overline{\Omega}$ such that $u$ remains bounded in $U$ throughout evolution.


title: Classical and Quantum Iterative Optimization Algorithms Based on Matrix Legendre-Bregman Projections

abstract: We consider Legendre-Bregman projections defined on the Hermitian matrix space and design iterative optimization algorithms based on them. A general duality theorem is established for Bregman divergences on Hermitian matrices, and it plays a crucial role in proving the convergence of the iterative algorithms. We study both exact and approximate Bregman projection algorithms. In the particular case of Kullback-Leibler divergence, our approximate iterative algorithm gives rise to the non-commutative versions of both the generalized iterative scaling (GIS) algorithm for maximum entropy inference and the AdaBoost algorithm in machine learning as special cases. As the Legendre-Bregman projections are simple matrix functions on Hermitian matrices, quantum algorithmic techniques are applicable to achieve potential speedups in each iteration of the algorithm. We discuss several quantum algorithmic design techniques applicable in our setting, including the smooth function evaluation technique, two-phase quantum minimum finding, and NISQ Gibbs state preparation.


title: How unitizing affects annotation of cohesion

abstract: This paper investigates how unitizing affects external observers' annotation of group cohesion. We compared unitizing techniques belonging to these categories: interval coding, continuous coding, and a technique inspired by a cognitive theory on event perception. We applied such techniques for sampling coding units from a set of recordings of social interactions rich in behaviors related to cohesion. Then, we compared the cohesion scores the observers assigned to each coding unit. Results show that the three techniques can lead to suitable ratings and that the technique inspired to cognitive theories leads to scores reflecting variability in cohesion better than the other ones.


title: Object oriented data analysis of surface motion time series in peatland landscapes

abstract: Peatlands account for 10% of UK land area, 80% of which are degraded to some degree, emitting carbon at a similar magnitude to oil refineries or landfill sites. A lack of tools for rapid and reliable assessment of peatland condition has limited monitoring of vast areas of peatland and prevented targeting areas urgently needing action to halt further degradation. Measured using interferometric synthetic aperture radar (InSAR), peatland surface motion is highly indicative of peatland condition, largely driven by the eco-hydrological change in the peatland causing swelling and shrinking of the peat substrate. The computational intensity of recent methods using InSAR time series to capture the annual functional structure of peatland surface motion becomes increasingly challenging as the sample size increases. Instead, we utilize the behavior of the entire peatland surface motion time series using object oriented data analysis to assess peatland condition. In a Gibbs sampling scheme, our cluster analysis based on the functional behavior of the surface motion time series finds features representative of soft/wet peatlands, drier/shrubby peatlands and thin/modified peatlands align with the clusters. The posterior distribution of the assigned peatland types enables the scale of peatland degradation to be assessed, which will guide future cost-effective decisions for peatland restoration.


title: On fits to correlated and auto-correlated data

abstract: Observables in particle physics and specifically in lattice QCD calculations are often extracted from fits. Standard $\chi^2$ tests require a reliable determination of the covariance matrix and its inverse from correlated and auto-correlated data, a challenging task often leading to close-to-singular estimates. These motivate modifications of the definition of $\chi^2$ such as uncorrelated fits. We show how the goodness-of-fit measured by their p-value can still be estimated robustly for a broad class of such fits.

categories: hep-lat cond-mat.stat-mech hep-ex hep-ph physics.comp-ph


title: Class-Imbalanced Complementary-Label Learning via Weighted Loss

abstract: Complementary-label learning (CLL) is widely used in weakly supervised classification, but it faces a significant challenge in real-world datasets when confronted with class-imbalanced training samples. In such scenarios, the number of samples in one class is considerably lower than in other classes, which consequently leads to a decline in the accuracy of predictions. Unfortunately, existing CLL approaches have not investigate this problem. To alleviate this challenge, we propose a novel problem setting that enables learning from class-imbalanced complementary labels for multi-class classification. To tackle this problem, we propose a novel CLL approach called Weighted Complementary-Label Learning (WCLL). The proposed method models a weighted empirical risk minimization loss by utilizing the class-imbalanced complementary labels, which is also applicable to multi-class imbalanced training samples. Furthermore, we derive an estimation error bound to provide theoretical assurance. To evaluate our approach, we conduct extensive experiments on several widely-used benchmark datasets and a real-world dataset, and compare our method with existing state-of-the-art methods. The proposed approach shows significant improvement in these datasets, even in the case of multiple class-imbalanced scenarios. Notably, the proposed method not only utilizes complementary labels to train a classifier but also solves the problem of class imbalance.


title: Predicting the FCI energy of large systems to chemical accuracy from restricted active space density matrix renormalization group calculations

abstract: We theoretically derive and validate with large scale simulations a remarkably accurate power law scaling of errors for the restricted active space density matrix renormalization group (DMRG-RAS) method [arXiv:2111.06665] in electronic structure calculations. This yields a new extrapolation method, DMRG-RAS-X, which reaches chemical accuracy for strongly correlated systems such as the Chromium dimer, dicarbon up to a large cc-pVQZ basis, and even a large chemical complex like the FeMoco with significantly lower computational demands than previous methods. The method is free of empirical parameters, performed robustly and reliably in all examples we tested, and has the potential to become a vital alternative method for electronic structure calculations in quantum chemistry, and more generally for the computation of strong correlations in nuclear and condensed matter physics.


title: Parameter identification from single trajectory data: from linear to nonlinear

abstract: Our recent work lays out a general framework for inferring information about the parameters and associated dynamics of a differential equation model from a discrete set of data points collected from the system being modeled. Rigorous mathematical results have justified this approach and have identified some common features that arise for certain classes of integrable models. In this work we present a thorough numerical investigation that shows that several of these core features extend to a paradigmatic linear-in-parameters model, the Lotka-Volterra (LV) system, which we consider in the conservative case as well as under the addition of terms that perturb the system away from this regime. A central construct for this analysis is a concise representation of parameter features in the data space that we call the $P_n$-diagram, which is particularly useful for visualization of results for low-dimensional (small $n$) systems. Our work also exposes some new properties related to non-uniqueness that arise for these LV systems, with non-uniqueness manifesting as a multi-layered structure in the associated $P_2$-diagrams.


title: Atomistic modelling of the channeling process with radiation reaction force included

abstract: Methodology is developed that incorporates the radiation reaction force into the relativistic molecular dynamics framework implemented in the MBN Explorer software package. The force leads to a gradual decrease in the projectile's energy E due to the radiation emission. This effect is especially strong for ultra-relativistic projectiles passing through oriented crystals where they experience the action of strong electrostatic fields as has been shown in recent experiments. A case study has been carried out for the initial approbation of the methodology developed. Simulations of the processes of planar channeling and photon emission have been performed for 150 GeV positrons in a 200 microns thick single oriented Si(110) crystal. Several regimes for the decrease in E have been established and characterized. Further steps in developing the code to include the necessary quantum corrections are identified and possible algorithmic modifications are proposed.


title: Optimal Sobolev embeddings for the Ornstein-Uhlenbeck operator

abstract: A comprehensive analysis of Sobolev-type inequalities for the Ornstein-Uhlenbeck operator in the Gauss space is offered. A unified approach is proposed, providing one with criteria for their validity in the class of rearrangement-invariant function norms. Optimal target and domain norms in the relevant inequalities are characterized via a reduction principle to one-dimensional inequalities for a Calder\'on type integral operator patterned on the Gaussian isoperimetric function. Consequently, the best possible norms in a variety of specific families of spaces, including Lebesgue, Lorentz, Lorentz-Zygmund, Orlicz and Marcinkiewicz spaces, are detected. The reduction principle hinges on a preliminary discussion of the existence and uniqueness of generalized solutions to equations, in the Gauss space, for the Ornstein-Uhlenbeck operator, with a just integrable right-hand side. A decisive role is also played by a pointwise estimate, in rearrangement form, for these solutions.


title: Galaxy interactions in filaments and sheets: effects of the large-scale structures versus the local density

abstract: The major interactions are known to trigger star formation in galaxies and alter their colour. We study the major interactions in filaments and sheets using the SDSS data to understand the influence of large-scale environments on the galaxy interactions. We identify the galaxies in filaments and sheets using the local dimension and also find the major pairs residing in these environments. The star formation rate and colour of the interacting galaxies as a function of pair separation are separately analyzed in filaments and sheets. The analysis is repeated for three volume limited samples covering different magnitude ranges. The major pairs residing in the filaments show a significantly higher star formation rate (SFR) and bluer colour than those residing in the sheets up to the projected pair separation of $\sim 50$ kpc. We observe a complete reversal of this behaviour for both the SFR and colour of the galaxy pairs having a projected separation larger than 50 kpc. Some earlier studies report that the galaxy pairs align with the filament axis. Such alignment inside filaments indicates anisotropic accretion that may cause these differences. We do not observe these trends in the brighter galaxy samples. The pairs in filaments and sheets from the brighter galaxy samples trace relatively denser regions in these environments. The absence of these trends in the brighter samples may be explained by the dominant effect of the local density over the effects of the large-scale environment.


title: Secure Indoor Location for Airport Environments

abstract: This work presents a secure novel solution based on inertial measurement units to provide indoor location and positioning in airports. The use of different technologies allows to locate people with precision in this kind of indoor places where the use of the GPS is not possible. The system has been developed thinking in the low cost and in a possible future expansion of this kind of systems to improve the Quality of Service of the users in airports. The use of QR codes and low cost IMU devices through the use of people smartphones ensure this premise. An Android application has been developed to show the applicability and performance of the system. The security in this kind of systems is essential. This kind of systems needs to avoid the traceability of the IMU devices when users are using it. To solve this problem, the FourQ elliptic curve has been used to generate a shared key using the elliptic curve Diffie-Hellman protocol. The key generated with the FourQ is used then to cipher all communications through the use of the SNOW 3G stream cipher. The developed system offers promising results.


title: Quantitative Khintchine in Simultaneous Approximation

abstract: In a ground-breaking work \cite{BY}, Beresnevich and Yang recently proved Khintchine's theorem in simultaneous Diophantine approximation for nondegenerate manifolds resolving a long-standing problem in the theory of Diophantine approximation. In this paper, we prove an effective version of their result.


title: On Computing Exact Means of Time Series Using the Move-Split-Merge Metric

abstract: Computing an accurate mean of a set of time series is a critical task in applications like nearest-neighbor classification and clustering of time series. While there are many distance functions for time series, the most popular distance function used for the computation of time series means is the non-metric dynamic time warping (DTW) distance. A recent algorithm for the exact computation of a DTW-Mean has a running time of $\mathcal{O}(n^{2k+1}2^kk)$, where $k$ denotes the number of time series and $n$ their maximum length. In this paper, we study the mean problem for the move-split-merge (MSM) metric that not only offers high practical accuracy for time series classification but also carries of the advantages of the metric properties that enable further diverse applications. The main contribution of this paper is an exact and efficient algorithm for the MSM-Mean problem of time series. The running time of our algorithm is $\mathcal{O}(n^{k+3}2^k k^3 )$, and thus better than the previous DTW-based algorithm. The results of an experimental comparison confirm the running time superiority of our algorithm in comparison to the DTW-Mean competitor. Moreover, we introduce a heuristic to improve the running time significantly without sacrificing much accuracy.


title: Graph Universal Cycles: Compression and Connections to Universal Cycles

abstract: Universal cycles, such as De Bruijn cycles, are cyclic sequences of symbols that represent every combinatorial object from some family exactly once as a consecutive subsequence. Graph universal cycles are a graph analogue of universal cycles introduced in 2010. We introduce graph universal partial cycles, a more compact representation of graph classes, which use "do not know" edges. We show how to construct graph universal partial cycles for labeled graphs, threshold graphs, and permutation graphs. For threshold graphs and permutation graphs, we demonstrate that the graph universal cycles and graph universal partial cycles are closely related to universal cycles and compressed universal cycles, respectively. Using the same connection, for permutation graphs, we define and prove the existence of an $s$-overlap form of graph universal cycles. We also prove the existence of a generalized form of graph universal cycles for unlabeled graphs.


title: Active Transfer Prototypical Network: An Efficient Labeling Algorithm for Time-Series Data

abstract: The paucity of labeled data is a typical challenge in the automotive industry. Annotating time-series measurements requires solid domain knowledge and in-depth exploratory data analysis, which implies a high labeling effort. Conventional Active Learning (AL) addresses this issue by actively querying the most informative instances based on the estimated classification probability and retraining the model iteratively. However, the learning efficiency strongly relies on the initial model, resulting in the trade-off between the size of the initial dataset and the query number. This paper proposes a novel Few-Shot Learning (FSL)-based AL framework, which addresses the trade-off problem by incorporating a Prototypical Network (ProtoNet) in the AL iterations. The results show an improvement, on the one hand, in the robustness to the initial model and, on the other hand, in the learning efficiency of the ProtoNet through the active selection of the support set in each iteration. This framework was validated on UCI HAR/HAPT dataset and a real-world braking maneuver dataset. The learning performance significantly surpasses traditional AL algorithms on both datasets, achieving 90% classification accuracy with 10% and 5% labeling effort, respectively.


title: Building an Ethereum-Based Decentralized Vehicle Rental System

abstract: Blockchain technology, beyond cryptocurrencies, is called to be the new information exchange ecosystem due to its unique properties, such as immutability and transparency. The main objective of this work is to introduce the design of a decentralized rental system, which leverages smart contracts and the Ethereum public blockchain. The work started from an exhaustive investigation on the Ethereum platform, emphasizing the aspect of cryptography and all the technology behind this platform. In order to test the proposed scheme in a realistic use, the implementation of a web application for the rental of vehicles has been carried out. The application covers the entire vehicle rental process offered in traditional web applications, adding more autonomy and ease of use to users. Following Ethereum application development guidelines, all business logic is located in the smart contracts implemented in the Ethereum network, where these contracts control the entire vehicle rental system of customers. While this is a work in progress, the results obtained in the first proof of concept have been very promising.


title: Spatial Pruned Sparse Convolution for Efficient 3D Object Detection

abstract: 3D scenes are dominated by a large number of background points, which is redundant for the detection task that mainly needs to focus on foreground objects. In this paper, we analyze major components of existing sparse 3D CNNs and find that 3D CNNs ignore the redundancy of data and further amplify it in the down-sampling process, which brings a huge amount of extra and unnecessary computational overhead. Inspired by this, we propose a new convolution operator named spatial pruned sparse convolution (SPS-Conv), which includes two variants, spatial pruned submanifold sparse convolution (SPSS-Conv) and spatial pruned regular sparse convolution (SPRS-Conv), both of which are based on the idea of dynamically determining crucial areas for redundancy reduction. We validate that the magnitude can serve as important cues to determine crucial areas which get rid of the extra computations of learning-based methods. The proposed modules can easily be incorporated into existing sparse 3D CNNs without extra architectural modifications. Extensive experiments on the KITTI, Waymo and nuScenes datasets demonstrate that our method can achieve more than 50% reduction in GFLOPs without compromising the performance.


title: Making statistics work: a quantum engine in the BEC-BCS crossover

abstract: Heat engines convert thermal energy into mechanical work both in the classical and quantum regimes. However, quantum theory offers genuine nonclassical forms of energy, different from heat, which so far have not been exploited in cyclic engines to produce useful work. We here experimentally realize a novel quantum many-body engine fuelled by the energy difference between fermionic and bosonic ensembles of ultracold particles that follows from the Pauli exclusion principle. We employ a harmonically trapped superfluid gas of $^6$Li atoms close to a magnetic Feshbach resonance which allows us to effectively change the quantum statistics from Bose-Einstein to Fermi-Dirac. We replace the traditional heating and cooling strokes of a quantum Otto cycle by tuning the gas between a Bose- Einstein condensate of bosonic molecules and a unitary Fermi gas (and back) through a magnetic field. The quantum nature of such a Pauli engine is revealed by contrasting it to a classical thermal engine and to a purely interaction-driven device. We obtain a work output of several $10^6$ vibrational quanta per cycle with an efficiency of up to $25\%$. Our findings establish quantum statistics as a useful thermodynamic resource for work production, shifting the paradigm of energy-conversion devices to a new class of emergent quantum engines.


abstract: The unitary correspondence between Quantum Hall states in higher Landau levels and states in the lowest Landau level is discussed together with the resulting transformation formulas for particle densities and interaction potentials. This correspondence leads in particular to a representation of states in arbitrary Landau levels in terms of holomorphic functions. Some important special Quantum Hall states in higher Landau levels are also briefly discussed.


title: Determination of $\alpha_s(m_Z)$ by the non-perturbative decoupling method

abstract: We present the details and first results of a new strategy for the determination of $\alpha_s(m_Z)$. By simultaneously decoupling 3 fictitious heavy quarks we establish a relation between the $\Lambda$-parameters of three-flavor QCD and pure gauge theory. Very precise recent results in the pure gauge theory can thus be leveraged to obtain the three-flavour $\Lambda$-parameter in units of a common decoupling scale. Connecting this scale to hadronic physics in 3-flavour QCD leads to our result in physical units, $\Lambda^{(3)}_{\bar{\rm MS}} = 336(12)\, {\rm MeV}$, which translates to $\alpha_s(m_Z) = 0.11823(84)$. This is compatible with both the FLAG average and the previous ALPHA result, with a comparable, yet still statistics dominated, error. This constitutes a highly non-trivial check, as the decoupling strategy is conceptually very different from the 3-flavour QCD step-scaling method, and so are their systematic errors. These include the uncertainties of the combined decoupling and continuum limits, which we discuss in some detail. We also quantify the correlation between both results, due to some common elements, such as the scale determination in physical units and the definition of the energy scale where we apply decoupling.


title: Prompt-driven efficient Open-set Semi-supervised Learning

abstract: Open-set semi-supervised learning (OSSL) has attracted growing interest, which investigates a more practical scenario where out-of-distribution (OOD) samples are only contained in unlabeled data. Existing OSSL methods like OpenMatch learn an OOD detector to identify outliers, which often update all modal parameters (i.e., full fine-tuning) to propagate class information from labeled data to unlabeled ones. Currently, prompt learning has been developed to bridge gaps between pre-training and fine-tuning, which shows higher computational efficiency in several downstream tasks. In this paper, we propose a prompt-driven efficient OSSL framework, called OpenPrompt, which can propagate class information from labeled to unlabeled data with only a small number of trainable parameters. We propose a prompt-driven joint space learning mechanism to detect OOD data by maximizing the distribution gap between ID and OOD samples in unlabeled data, thereby our method enables the outliers to be detected in a new way. The experimental results on three public datasets show that OpenPrompt outperforms state-of-the-art methods with less than 1% of trainable parameters. More importantly, OpenPrompt achieves a 4% improvement in terms of AUROC on outlier detection over a fully supervised model on CIFAR10.


title: A multiscale and multiphysics framework to simulate radiation damage in nano-crystalline materials

abstract: This work presents a multiscale and multiphysics framework to investigate the radiation-induced damage in nano-crystalline materials. The framework combines two methodologies, including molecular dynamics simulations with electronic effects and long-term atomistic diffusion simulations in nano-crystalline materials. Using this framework, we investigated nano-crystalline materials' self-healing behavior under radiation events. We found that the number of defects generated in nano-crystals during the cascade simulations was less than in single crystals. This behavior was due to the fast absorption of interstitial atoms in the grain boundary network during the cascade simulations, while vacancies migrated to the boundaries in a much longer time scale than interstitial atoms. Thus, nano-crystalline materials showed a self-healing behavior where the number and size of the defects are drastically reduced with time. We found that the self-healing behavior of nano-crystalline materials is limited, and about 50% of vacancies survived. This effect resulted from clusters of vacancies' collective behavior, which are much more stable than individual vacancies.


title: A Fully Homomorphic Encryption Scheme for Real-Time Safe Control

abstract: Fully Homomorphic Encryption (FHE) has made it possible to perform addition and multiplication operations on encrypted data. Using FHE in control thus has the advantage that control effort for a plant can be calculated remotely without ever decrypting the exchanged information. FHE in its current form is however not practically applicable for real-time control as its computational load is very high compared to traditional encryption methods. In this paper a reformulation of the Gentry FHE scheme is proposed and applied on an FPGA to solve this problem. It is shown that the resulting FHE scheme can be implemented for real-time stabilization of an inverted double pendulum using discrete time control.


title: Optimality problems in Orlicz spaces

abstract: In mathematical modelling, the data and solutions are represented as measurable functions and their quality is oftentimes captured by the membership to a certain function space. One of the core questions for an analysis of a model is the mutual relationship between the data and solution quality. The optimality of the obtained results deserves a special focus. It requires a careful choice of families of function spaces balancing between their expressivity, i.e. the ability to capture fine properties of the model, and their accessibility, i.e. its technical difficulty for practical use. This paper presents a unified and general approach to optimality problems in Orlicz spaces. Orlicz spaces are parametrized by a single convex function and neatly balance the expressivity and accessibility. We prove a general principle that yields an easily verifiable necessary and sufficient condition for the existence or the non-existence of an optimal Orlicz space in various tasks. We demonstrate its use in specific problems, including the continuity of Sobolev embeddings and boundedness of integral operators such as the Hardy--Littlewood maximal operator and the Laplace transform.


title: Nature and origin of the operators entering the master equation of an open quantum system

abstract: By exploiting the peculiarities of a recently introduced formalism for describing open quantum systems (the Parametric Representation with Environmental Coherent States) we derive an equation of motion for the reduced density operator of an open quantum system that has the same structure of the celebrated Gorini-Kossakowski-Sudarshan-Lindblad equation, but holds regardless of markovianity being assumed. The operators in our result have explicit expressions in terms of the Hamiltonian describing the interactions with the environment, and can be computed once a specific model is considered. We find that, instead of a single set of Lindblad operators, in the general (non-markovian) case there one set of Lindblad-like operators for each and every point of a symplectic manifold associated to the environment. This intricacy disappears under some assumptions (which are related to markovianity and the classical limit of the environment), under which it is possible to recover the usual master equation formalism. Finally, we find such Lindblad-like operators for two different models of a qubit in a bosonic environment, and show that in the classical limit of the environment their renown master equations are recovered.


title: Understanding the impacts and the influences of Metro rail on urban environment Case studies and the Bengaluru scenario

abstract: Urban sprawl is brought on by migration brought on by the new and exciting prospects in the cities. Cities expand rapidly, which leads to increased population growth and exceeding vehicle expansion. To better understand the features and advantages of Metro Transit surveys on car ownership transit use frequency and trip purpose are crucial. These surveys are nothing more than case studies of specific transit operators from all throughout the country. Case studies explore the various city conditions the role of MRTS problems which are created by them problems solved by them and to envisage the future problems over a period of time. They also provide the clarity of the various issues and challenges faced by the city at different time of the day. Bengaluru is an outstanding example how a pensioners paradise has resulted now as a Meta city. It has undergone a huge metamorphosis trying to solve the various demands put forth by the economic activities. Transit oriented development TOD seems the logical way to address the traffic and transportation problems which results in a mixed-use development high rise high density. After Delhi and Hyderabad Bengaluru has entered the queue to embrace Namma Metro and implement Phases 2A and 2B after Phases 1 and 2. An inventory into selected case studies will relate and differentiate the conditions and systems in Bengaluru thus enriching our capability to tackle the Transit Oriented Development TOD issues at bay. Key words Public Transportation Metro rail Case study Namma Metro Bengaluru TOD


title: On Gelfand graded commutative rings

abstract: This paper deals with the graded commutative rings in which every homogeneous prime ideal is contained in a unique homogeneous maximal ideal called Gelfand graded ring. The purpose is to establish some topological and algebraic characterizations of these rings, one of which is the algebraic analogue of the Urysohn's lemma. Finally we look at a special class of those graded rings called pm$^+$ graded rings which can be viewed as graded ring with a Gelfand strong property.


title: Automated Quality Controlled Analysis of 2D Phase Contrast Cardiovascular Magnetic Resonance Imaging

abstract: Flow analysis carried out using phase contrast cardiac magnetic resonance imaging (PC-CMR) enables the quantification of important parameters that are used in the assessment of cardiovascular function. An essential part of this analysis is the identification of the correct CMR views and quality control (QC) to detect artefacts that could affect the flow quantification. We propose a novel deep learning based framework for the fully-automated analysis of flow from full CMR scans that first carries out these view selection and QC steps using two sequential convolutional neural networks, followed by automatic aorta and pulmonary artery segmentation to enable the quantification of key flow parameters. Accuracy values of 0.958 and 0.914 were obtained for view classification and QC, respectively. For segmentation, Dice scores were $>$0.969 and the Bland-Altman plots indicated excellent agreement between manual and automatic peak flow values. In addition, we tested our pipeline on an external validation data set, with results indicating good robustness of the pipeline. This work was carried out using multivendor clinical data consisting of 986 cases, indicating the potential for the use of this pipeline in a clinical setting.


title: On abelian and cyclic group codes

abstract: We determine a condition on the minimum Hamming weight of some special abelian group codes and, as a consequence of this result, we establish that any such code is, up to permutational equivalence, a subspace of the direct sum of $s$ copies of the repetition code of length $t$, for some suitable positive integers $s$ and $t$. Moreover, we provide a complete characterisation of permutation automorphisms of the linear code $C=\bigoplus_{i=1}^{s}Rep_{t}(\mathbb{F}_{q})$ and we establish that such a code is an abelian group code, for every pair of integers $s,t\geq1$. Finally, in a similar fashion as for abelian group codes, we give an equivalent characterisation of cyclic group codes.


title: Optical circular polarization induced by axionlike particles in blazars

abstract: We propose that the interaction between the axionlike particles (ALPs) and photons can be a possible origin of optical circular polarization (CP) in blazars. Given that there is no definite detection of optical CP at $\sim0.1\%$ level, a rough limit on ALP-photon coupling can be obtained, specifically $g_{a\gamma}\cdot B_\mathrm{T0}\lesssim7.9\times10^{-12}~\mathrm{G\cdot GeV}^{-1}$ for $m_{a}\lesssim 10^{-13}~\mathrm{eV}$, depending on the magnetic field configuration of the blazar jet. Obviously, for the blazar models with a larger magnetic field strength, such as hadronic radiation models, this constraint could be more stringent. We also perform a dedicated analysis of the tentative observations of optical CP in two blazars, namely 3C 66A and OJ 287, and we find that these observations could be explained by the ALP-photon mixing with $g_{a\gamma} \sim 10^{-11}~\mathrm{GeV}^{-1}$. As an outlook, our analysis can be improved by further research on the radiation models of blazars and high-precision joint measurements of optical CP and linear polarization.


title: Quantum Hall Phases of Cold Bose Gases

abstract: Cold atomic gases of interacting bosons subject to rapid rotation and confined in anharmonic traps can theoretically exhibit analogues of the fractional quantum Hall effect for electrons in strong magnetic fields. In this setting the Coriolis force due to the rotation mimics the Lorentz force on charged particles but artificial gauge fields can also be obtained by coupling the internal structure of the atoms to light fields. The chapter discusses mathematical aspects of transitions to different strongly correlated phases that appear when the parameters of a model Hamiltonian are varied.


title: Shining a Light on Forensic Black-Box Studies

abstract: Forensic science plays a critical role in the United States criminal justice system. For decades, many feature-based fields of forensic science, such as firearm and toolmark identification, developed outside the scientific community's purview. The results of these studies are widely relied on by judges nationwide. However, this reliance is misplaced. Black-box studies to date suffer from inappropriate sampling methods and high rates of missingness. Current black-box studies ignore both problems in arriving at the error rate estimates presented to courts. We explore the impact of each type of limitation using available data from black-box studies and court materials. We show that black-box studies rely on non-representative samples of examiners. Using a case study of a popular ballistics study, we find evidence that these unrepresentative samples may commit fewer errors than the wider population from which they came. We also find evidence that the missingness in black-box studies is non-ignorable. Using data from a recent latent print study, we show that ignoring this missingness likely results in systematic underestimates of error rates. Finally, we offer concrete steps to overcome these limitations.


title: Longitudinal Variability Analysis on Low-dose Abdominal CT with Deep Learning-based Segmentation

abstract: Metabolic health is increasingly implicated as a risk factor across conditions from cardiology to neurology, and efficiency assessment of body composition is critical to quantitatively characterizing these relationships. 2D low dose single slice computed tomography (CT) provides a high resolution, quantitative tissue map, albeit with a limited field of view. Although numerous potential analyses have been proposed in quantifying image context, there has been no comprehensive study for low-dose single slice CT longitudinal variability with automated segmentation. We studied a total of 1816 slices from 1469 subjects of Baltimore Longitudinal Study on Aging (BLSA) abdominal dataset using supervised deep learning-based segmentation and unsupervised clustering method. 300 out of 1469 subjects that have two year gap in their first two scans were pick out to evaluate longitudinal variability with measurements including intraclass correlation coefficient (ICC) and coefficient of variation (CV) in terms of tissues/organs size and mean intensity. We showed that our segmentation methods are stable in longitudinal settings with Dice ranged from 0.821 to 0.962 for thirteen target abdominal tissues structures. We observed high variability in most organ with ICC<0.5, low variability in the area of muscle, abdominal wall, fat and body mask with average ICC>0.8. We found that the variability in organ is highly related to the cross-sectional position of the 2D slice. Our efforts pave quantitative exploration and quality control to reduce uncertainties in longitudinal analysis.


title: DMAP: a Distributed Morphological Attention Policy for Learning to Locomote with a Changing Body

abstract: Biological and artificial agents need to deal with constant changes in the real world. We study this problem in four classical continuous control environments, augmented with morphological perturbations. Learning to locomote when the length and the thickness of different body parts vary is challenging, as the control policy is required to adapt to the morphology to successfully balance and advance the agent. We show that a control policy based on the proprioceptive state performs poorly with highly variable body configurations, while an (oracle) agent with access to a learned encoding of the perturbation performs significantly better. We introduce DMAP, a biologically-inspired, attention-based policy network architecture. DMAP combines independent proprioceptive processing, a distributed policy with individual controllers for each joint, and an attention mechanism, to dynamically gate sensory information from different body parts to different controllers. Despite not having access to the (hidden) morphology information, DMAP can be trained end-to-end in all the considered environments, overall matching or surpassing the performance of an oracle agent. Thus DMAP, implementing principles from biological motor control, provides a strong inductive bias for learning challenging sensorimotor tasks. Overall, our work corroborates the power of these principles in challenging locomotion tasks.


title: Strong-$CP$ with and without gravity

abstract: Conventionally, the strong-$CP$ problem is assumed to be a naturalness puzzle, with the axion solution sometimes viewed as an ad hoc fix. Gravity is either ignored or taken as a threat for the global Peccei-Quinn symmetry. We explain that the situation is fundamentally different. In gravity, axion is a matter of consistency imposed by the $S$-matrix: Each gauge sector must include axion with exact relaxation of the corresponding $\bar{\theta}$. We show that this favors an alternative and remarkably simple formulation of the axion, fully fixed by the gauge redundancy of QCD, without involvement of a global symmetry. The axion mechanism is a Higgs effect for the QCD $3$-form, ensuring that physics is independent of $\bar{\theta}$ to all orders in operator expansion. A near-future experimental detection of the neutron EDM will be an unambiguous signal of $CP$-violating physics beyond the Standard Model. The axion coupling is tied to the scale of gravity.


title: New Constraints on Warm Dark Matter from the Lyman-$\alpha$ Forest Power Spectrum

abstract: The forest of Lyman-$\alpha$ absorption lines detected in the spectra of distant quasars encodes information on the nature and properties of dark matter and the thermodynamics of diffuse baryonic material. Its main observable -- the 1D flux power spectrum (FPS) -- should exhibit a suppression on small scales and an enhancement on large scales in warm dark matter (WDM) cosmologies compared to standard $\Lambda$CDM. Here, we present an unprecedented suite of 1080 high-resolution cosmological hydrodynamical simulations run with the Graphics Processing Unit-accelerated code Cholla to study the evolution of the Lyman-$\alpha$ forest under a wide range of physically-motivated gas thermal histories along with different free-streaming lengths of WDM thermal relics in the early Universe. A statistical comparison of synthetic data with the forest FPS measured down to the smallest velocity scales ever probed at redshifts $4.0\lesssim z\lesssim 5.2$ (Boera et al. 2019) yields a lower limit $m_{\rm WDM}>3.1$ keV (95 percent CL) for the WDM particle mass and constrains the amplitude and spectrum of the photoheating and photoionizing background produced by star-forming galaxies and active galactic nuclei at these redshifts. Interestingly, our Bayesian inference analysis appears to weakly favor WDM models with a peak likelihood value at the thermal relic mass of $m_{\rm WDM}=4.5$ keV. We find that the suppression of the FPS from free-streaming saturates at $k\gtrsim 0.1\,\mathrm{s}\,\mathrm{km}^{-1}$ because of peculiar velocity smearing, and this saturated suppression combined with a slightly lower gas temperature provides a moderately better fit to the observed small-scale FPS for WDM cosmologies.


title: Exploiting oriented field projectors to open topological gaps in plasmonic nanoparticle arrays

abstract: In the last years there have been multiple proposals in nanophotonics to mimic topological condensed matter systems. However, nanoparticles have degrees of freedom that atoms lack of, like dimensions or shape, which can be exploited to explore topology beyond electronics. Elongated nanoparticles can act like projectors of the electric field in the direction of the major axis. Then, by orienting them in an array the coupling between them can be tuned, allowing to open a gap in an otherwise gapless system. As a proof of the potential of the use of orientation of nanoparticles for topology, we study 1D chains of prolate spheroidal silver nanoparticles. We show that in these arrays spatial modulation of the polarization allows to open gaps, engineer hidden crystalline symmetries and to switch on/off or left/right edge states depending on the polarization of the incident electric field. This opens a path towards exploiting features of nanoparticles for topology to go beyond analogues of condensed matter systems.


title: Online Subset Selection using $\alpha$-Core with no Augmented Regret

abstract: We revisit the classic problem of optimal subset selection in the online learning set-up. Assume that the set $[N]$ consists of $N$ distinct elements. On the $t$th round, an adversary chooses a monotone reward function $f_t: 2^{[N]} \to \mathbb{R}_+$ that assigns a non-negative reward to each subset of $[N].$ An online policy selects (perhaps randomly) a subset $S_t \subseteq [N]$ consisting of $k$ elements before the reward function $f_t$ for the $t$th round is revealed to the learner. As a consequence of its choice, the policy receives a reward of $f_t(S_t)$ on the $t$th round. Our goal is to design an online sequential subset selection policy to maximize the expected cumulative reward accumulated over a time horizon. In this connection, we propose an online learning policy called SCore (Subset Selection with Core) that solves the problem for a large class of reward functions. The proposed SCore policy is based on a new polyhedral characterization of the reward functions called $\alpha$-Core - a generalization of Core from the cooperative game theory literature. We establish a learning guarantee for the SCore policy in terms of a new performance metric called $\alpha$-augmented regret. In this new metric, the performance of the online policy is compared with an unrestricted offline benchmark that can select all $N$ elements at every round. We show that a large class of reward functions, including submodular, can be efficiently optimized with the SCore policy. We also extend the proposed policy to the optimistic learning set-up where the learner has access to additional untrusted hints regarding the reward functions. Finally, we conclude the paper with a list of open problems.


title: Subsonic and supersonic mechanisms in compressible turbulent boundary layers: a perspective from resolvent analysis

abstract: In this study we isolate the mechanisms that are separately responsible for amplifying the two distinct sets of modes that appear in a linearized Navier-Stokes based model for compressible turbulent boundary layer flows. These two sets of modes are: (i) the subsonic modes that have equivalents in the incompressible flow regime, and (ii) the supersonic modes that do not have such equivalents in the incompressible flow regime. The resolvent analysis framework is used to analyze the linear model, where the non-linear terms of the equations are taken to be a forcing to the linear terms. The most sensitive forcing and correspondingly the most amplified response are then analyzed. We find that, when considering the subsonic response modes, only the solenoidal component of the forcing to the momentum equations can amplify these modes. This is consistent with observations in the literature regarding incompressible flows. Next, when considering the supersonic response modes, we find that these are pressure fluctuations that radiate into the freestream. Within the freestream, these modes closely follow the trends of the inviscid Mach waves. There are three different forcing mechanisms that are found to amplify these modes. First, the forcing to the continuity and energy equations contribute. Second, the dilatational component of the forcing to the momentum equations has a significant role. Finally, although less significant for the resolvent amplification, there is also a contribution from the solenoidal component of the forcing. Using these identified forcing mechanisms, a decomposition of the forcing to the compressible resolvent operator, i.e. of the non-linear terms of the linearized Navier-Stokes equations, is proposed, that can approximately isolate the subsonic and the supersonic modes in turbulent boundary layers over adiabatic as well as cooled walls, and for a range of Mach numbers.


title: Gradient flow step-scaling function for SU(3) with $N_f$ = 6 or 4 fundamental flavors

abstract: Nonperturbative determinations of the renormalization group (RG) $\beta$ function are crucial to understand properties of gauge-fermion systems at strong coupling and connect lattice simulations and the perturbative ultraviolet regime. Choosing well-understood, QCD-like systems with SU(3) gauge group and either six or four fundamental flavors, we investigate their step-scaling $\beta$ function. In both cases we push the simulations to the boundary of chiral symmetry breaking and study the regime $g^2_{GF} \lesssim 8.2$ with six, and $g^2_{GF} \lesssim 6.6$ with four flavors. We carefully consider the lattice discretization errors by comparing three different gradient flows (GF), and for each flow three operators to estimate the renormalized finite volume coupling. We also consider the tree level improvement of the coupling. Noteworthy outcome is that nonperturbatively determined $\beta$ functions run much slower than perturbatively predicted.


title: Road Rutting Detection using Deep Learning on Images

abstract: Road rutting is a severe road distress that can cause premature failure of road incurring early and costly maintenance costs. Research on road damage detection using image processing techniques and deep learning are being actively conducted in the past few years. However, these researches are mostly focused on detection of cracks, potholes, and their variants. Very few research has been done on the detection of road rutting. This paper proposes a novel road rutting dataset comprising of 949 images and provides both object level and pixel level annotations. Object detection models and semantic segmentation models were deployed to detect road rutting on the proposed dataset, and quantitative and qualitative analysis of model predictions were done to evaluate model performance and identify challenges faced in the detection of road rutting using the proposed method. Object detection model YOLOX-s achieves mAP@IoU=0.5 of 61.6% and semantic segmentation model PSPNet (Resnet-50) achieves IoU of 54.69 and accuracy of 72.67, thus providing a benchmark accuracy for similar work in future. The proposed road rutting dataset and the results of our research study will help accelerate the research on detection of road rutting using deep learning.


title: Radio source-component association for the LOFAR Two-metre Sky Survey with region-based convolutional neural networks

abstract: Radio loud active galactic nuclei (RLAGNs) are often morphologically complex objects that can consist of multiple, spatially separated, components. Astronomers often rely on visual inspection to resolve radio component association. However, applying visual inspection to all the hundreds of thousands of well-resolved RLAGNs that appear in the images from the Low Frequency Array (LOFAR) Two-metre Sky Survey (LoTSS) at $144$ MHz, is a daunting, time-consuming process, even with extensive manpower. Using a machine learning approach, we aim to automate the radio component association of large ($> 15$ arcsec) radio components. We turned the association problem into a classification problem and trained an adapted Fast region-based convolutional neural network to mimic the expert annotations from the first LoTSS data release. We implemented a rotation data augmentation to reduce overfitting and simplify the component association by removing unresolved radio sources that are likely unrelated to the large and bright radio components that we consider using predictions from an existing gradient boosting classifier. For large ($> 15$ arcsec) and bright ($> 10$ mJy) radio components in the LoTSS first data release, our model provides the same associations for $85.3\%\pm0.6$ of the cases as those derived when astronomers perform the association manually. When the association is done through public crowd-sourced efforts, a result similar to that of our model is attained. Our method is able to efficiently carry out manual radio-component association for huge radio surveys and can serve as a basis for either automated radio morphology classification or automated optical host identification. This opens up an avenue to study the completeness and reliability of samples of radio sources with extended, complex morphologies.


title: UTC Time, Formally Verified

abstract: FV Time is a small-scale verification project developed in the Coq proof assistant using the Mathematical Components libraries. It is a library for managing conversions between time formats (UTC and timestamps), as well as commonly used functions for time arithmetic. As a library for time conversions, its novelty is the implementation of leap seconds, which are part of the UTC standard but usually not implemented in existing libraries. Since the verified functions of FV Time are reasonably simple yet non-trivial, it nicely illustrates our methodology for verifying software with Coq. In this paper we present a description of the project, emphasizing the main problems faced while developing the library, as well as some general-purpose solutions that were produced as by-products and may be used in other verification projects. These include a refinement package between proof-oriented MathComp numbers and computation-oriented primitive numbers from the Coq standard library, as well as a set of tactics to automatically prove certain decidable statements over finite ranges through brute-force computation.


title: Knowledge-Aware Bayesian Deep Topic Model

abstract: We propose a Bayesian generative model for incorporating prior domain knowledge into hierarchical topic modeling. Although embedded topic models (ETMs) and its variants have gained promising performance in text analysis, they mainly focus on mining word co-occurrence patterns, ignoring potentially easy-to-obtain prior topic hierarchies that could help enhance topic coherence. While several knowledge-based topic models have recently been proposed, they are either only applicable to shallow hierarchies or sensitive to the quality of the provided prior knowledge. To this end, we develop a novel deep ETM that jointly models the documents and the given prior knowledge by embedding the words and topics into the same space. Guided by the provided knowledge, the proposed model tends to discover topic hierarchies that are organized into interpretable taxonomies. Besides, with a technique for adapting a given graph, our extended version allows the provided prior topic structure to be finetuned to match the target corpus. Extensive experiments show that our proposed model efficiently integrates the prior knowledge and improves both hierarchical topic discovery and document representation.


title: Process-guidance improves predictive performance of neural networks for carbon turnover in ecosystems

abstract: Despite deep-learning being state-of-the-art for data-driven model predictions, it has not yet found frequent application in ecology. Given the low sample size typical in many environmental research fields, the default choice for the modelling of ecosystems and its functions remain process-based models. The process understanding coded in these models complements the sparse data and neural networks can detect hidden dynamics even in noisy data. Embedding the process model in the neural network adds information to learn from, improving interpretability and predictive performance of the combined model towards the data-only neural networks and the mechanism-only process model. At the example of carbon fluxes in forest ecosystems, we compare different approaches of guiding a neural network towards process model theory. Evaluation of the results under four classical prediction scenarios supports decision-making on the appropriate choice of a process-guided neural network.


title: Predictable trajectories of the reduced Collatz iteration and a possible pathway to the proof of the Collatz conjecture (Version 2)

abstract: I show here that there are three different kinds of iterations for the reduced Collatz algorithm; depending on whether the root of the number is odd or even. There is only one kind of iteration if the root is odd and two kinds if the root is even. I also show that iterations on numbers with odd roots will cause an increase in value and eventually lead to an even rooted number. The iterations on even rooted numbers will subsequently cause a decrease in values. Because increase in values during the odd root iterations are bounded, I conclude that the Collatz iteration cannot veer to infinity. Since the sequence generated by the Collatz iteration is infinite and the values of the numbers do not veer to infinity it must either cycle and/or converge. I postulate that any cycling must occur with alternating types of iterations: e.g. odd rooted iterations which cause the values of the numbers to increase followed by even rooted iterations which causes the values to decrease. I show here that for simpler types of cycles, valid values of odd rooted or even rooted numbers are only found in a narrow gap which closes as the number of iterations increase. I further generalize to all types of odd-even and even-odd iterations. Given that previous work has shown that only very large non-trivial cycles are feasible during the Collatz iteration and this study shows the low probability of large simple cycles, leads us to the conclusion most likely cycles other than the trivial cycle are not possible during the Collatz iteration.


title: Constant intensity acoustic propagation in the presence of non-uniform properties and impedance discontinuities: Hermitian and non-Hermitian solutions

abstract: Propagation of sound through a non-uniform medium without scattering is possible, in principle, if the density and acoustic compressibility assume complex values, requiring passive and active mechanisms, also known as Hermitian and non-Hermitian solutions, respectively. Two types of constant intensity wave conditions are identified: in the first the propagating acoustic pressure has constant amplitude, while in the second the energy flux remains constant. The fundamental problem of transmission across an impedance discontinuity without reflection or energy loss is solved using a combination of monopole and dipole resonators in parallel. The solution depends on an arbitrary phase angle which can be chosen to give a unique acoustic metamaterial with both resonators undamped and passive, requiring purely Hermitian acoustic elements. For other phase angles one of the two elements must be active and the other passive, resulting in a gain/loss non-Hermitian system. These results prove that uni-directional and reciprocal transmission through a slab separating two half spaces is possible using passive Hermitian acoustic elements without the need to resort to active gain/loss energetic mechanisms.


title: GDCP: A custom package for flexible detector setups in Geant4

abstract: Geant4 Detector Construction Pattern (GDCP) is a C++ package containing a set of utility classes for developers to utilize during the detector construction phase of a Geant4 application. The main objective of this study is to provide an interface that helps users develop a highly flexible, understandable, and maintainable detector construction code. The package also contains a couple of helper classes facilitating the implementation of optical components in scintillation detectors. Finally, an advanced Geant application has been provided to demonstrate the package's utility. Program summary Program title: GDCP CPC Library link to program files: Licensing provisions: GNU General Public License 3 Programming language: C++ External routines/libraries: Geant4 Nature of problem: Geant4 provides an abstract class, G4VUserDetectorConstruction, that is used for the whole detector construction. This is sufficient for a simple detector setup but implementing a single class for building a complex detector leads to complexity in the code. To solve this issue it is recommended to split the implementation into more methods or classes; however, there is no standard solution published yet and it is entirely left to the users' ability. Solution method: To address this issue, I propose to apply the "Builder design pattern" in the detector construction phase. For this, the elements of the Builder design pattern have been introduced and made available for developers. Keywords: G4VUserDetectorConstruction, Builder design pattern, Detector construction, Component builder, Optical photon, Scintillator


title: Obstacle Identification and Ellipsoidal Decomposition for Fast Motion Planning in Unknown Dynamic Environments

abstract: Collision avoidance in the presence of dynamic obstacles in unknown environments is one of the most critical challenges for unmanned systems. In this paper, we present a method that identifies obstacles in terms of ellipsoids to estimate linear and angular obstacle velocities. Our proposed method is based on the idea of any object can be approximately expressed by ellipsoids. To achieve this, we propose a method based on variational Bayesian estimation of Gaussian mixture model, the Kyachiyan algorithm, and a refinement algorithm. Our proposed method does not require knowledge of the number of clusters and can operate in real-time, unlike existing optimization-based methods. In addition, we define an ellipsoid-based feature vector to match obstacles given two timely close point frames. Our method can be applied to any environment with static and dynamic obstacles, including the ones with rotating obstacles. We compare our algorithm with other clustering methods and show that when coupled with a trajectory planner, the overall system can efficiently traverse unknown environments in the presence of dynamic obstacles.


title: A Jordan-like decomposition for linear relations in finite-dimensional spaces

abstract: A square matrix $A$ has the usual Jordan canonical form that describes the structure of $A$ via eigenvalues and the corresponding Jordan blocks. If $A$ is a linear relation in a finite-dimensional linear space ${\mathfrak H}$ (i.e., $A$ is a linear subspace of ${\mathfrak H} \times {\mathfrak H}$ and can be considered as a multivalued linear operator), then there is a richer structure. In addition to the classical Jordan chains (interpreted in the Cartesian product ${\mathfrak H} \times {\mathfrak H}$), there occur three more classes of chains: chains starting at zero (the chains for the eigenvalue infinity), chains starting at zero and also ending at zero (the singular chains), and chains with linearly independent entries (the shift chains). These four types of chains give rise to a direct sum decomposition (a Jordan-like decomposition) of the linear relation $A$. In this decomposition there is a completely singular part that has the extended complex plane as eigenvalues; a usual Jordan part that corresponds to the finite proper eigenvalues; a Jordan part that corresponds to the eigenvalue $\infty$; and a multishift, i.e., a part that has no eigenvalues at all. Furthermore, the Jordan-like decomposition exhibits a certain uniqueness, closing a gap in earlier results. The presentation is purely algebraic, only the structure of linear spaces is used. Moreover, the presentation has a uniform character: each of the above types is constructed via an appropriately chosen sequence of quotient spaces. The dimensions of the spaces are the Weyr characteristics, which uniquely determine the Jordan-like decomposition of the linear relation.


title: Exact and efficient multivariate two-sample tests through generalized linear rank statistics

abstract: So-called linear rank statistics provide a means for distribution-free (even in finite samples), yet highly flexible, two-sample testing in the setting of univariate random variables. Their flexibility derives from a choice of weights that can be adapted to any given (simple) alternative hypothesis to achieve efficiency in case of correct specification of said alternative, while their non-parametric nature guarantees well-calibrated $p$-values even under misspecification. By drawing connections to (generalized) maximum likelihood estimation, and expanding on recent work on ranks in multiple dimensions, we extend linear rank statistics both to multivariate random variables and composite alternatives. Doing so yields non-parametric, multivariate two-sample tests that mirror efficiency properties of likelihood ratio tests, while remaining robust against model misspecification and computationally tractable. We prove non-parametric versions of the classical Wald and score tests facilitating hypothesis testing in the asymptotic regime, and relate these generalized linear rank statistics to linear spacing statistics enabling exact $p$-value computations in the small to moderate sample setting. Moreover, viewing rank statistics through the lens of likelihood ratios affords applications beyond fully efficient two-sample testing, of which we demonstrate three: testing in the presence of nuisance alternatives, simultaneous detection of location and scale shifts, and $K$-sample testing.


title: On Primality Tests Grounded on Binomial Coefficients

abstract: In this paper, we introduce two primality tests based on new divisibility properties of binomial coefficients. These new properties were enunciated and proved in previous work. We also study two similar tests that can be obtained from well-known results in Number Theory. At the end we compare our results with the existing ones.


title: Accuracy, Fairness, and Interpretability of Machine Learning Criminal Recidivism Models

abstract: Criminal recidivism models are tools that have gained widespread adoption by parole boards across the United States to assist with parole decisions. These models take in large amounts of data about an individual and then predict whether an individual would commit a crime if released on parole. Although such models are not the only or primary factor in making the final parole decision, questions have been raised about their accuracy, fairness, and interpretability. In this paper, various machine learning-based criminal recidivism models are created based on a real-world parole decision dataset from the state of Georgia in the United States. The recidivism models are comparatively evaluated for their accuracy, fairness, and interpretability. It is found that there are noted differences and trade-offs between accuracy, fairness, and being inherently interpretable. Therefore, choosing the best model depends on the desired balance between accuracy, fairness, and interpretability, as no model is perfect or consistently the best across different criteria.


title: Set-Valued Shadow Matching Using Zonotopes for 3-D Map-Aided GNSS Localization

abstract: Unlike many urban localization methods that return point-valued estimates, a set-valued representation enables robustness by ensuring that a continuum of possible positions obeys safety constraints. One strategy with the potential for set-valued estimation is GNSS-based shadow matching~(SM), where one uses a three-dimensional (3-D) map to compute GNSS shadows (where line-of-sight is blocked). However, SM requires a point-valued grid for computational tractability, with accuracy limited by grid resolution. We propose zonotope shadow matching (ZSM) for set-valued 3-D map-aided GNSS localization. ZSM represents buildings and GNSS shadows using constrained zonotopes, a convex polytope representation that enables propagating set-valued estimates using fast vector concatenation operations. Starting from a coarse set-valued position, ZSM refines the estimate depending on the receiver being inside or outside each shadow as judged by received carrier-to-noise density. We demonstrated our algorithm's performance using simulated experiments on a simple 3-D example map and on a dense 3-D map of San Francisco.


title: How to solve a classification problem using a cooperative tiling Multi-Agent System?

abstract: Adaptive Multi-Agent Systems (AMAS) transform dynamic problems into problems of local cooperation between agents. We present smapy, an ensemble based AMAS implementation for mobility prediction, whose agents are provided with machine learning models in addition to their cooperation rules. With a detailed methodology, we propose a framework to transform a classification problem into a cooperative tiling of the input variable space. We show that it is possible to use linear classifiers for online non-linear classification on three benchmark toy problems chosen for their different levels of linear separability, if they are integrated in a cooperative Multi-Agent structure. The results obtained show a significant improvement of the performance of linear classifiers in non-linear contexts in terms of classification accuracy and decision boundaries, thanks to the cooperative approach.


title: Two examples of ungrading in higher education from the United States and from Germany

abstract: In this paper the authors discuss their experiences with ungrading at a small public university in the U.S. as well as a large public university in Germany. The courses described are Calculus 1, Mathematics for Liberal Arts, and courses for pre-service secondary teachers of mathematics. We outline and compare our approaches, discuss student performance and feedback and present some interesting patterns relating to gender. A shortened and revised version of this paper appeared in PRIMUS 33 (2023), no. 9, 1035-1054, DOI: 10.1080/10511970.2023.2229819.


title: Cross Ratio Geometry Advances for Four Co-Linear Points in the Desargues Affine Plane-Skew Field

abstract: This paper introduces advances in the geometry of the cross ratio of four co-linear points in in the Desargues affine plane. The cross-ratio of co-linear points of a skew field in the Desargues affine plane. The results given here have a clean rendition, based on Desargues affine plane axiomatics, skew field properties and the addition and multiplication of planar co-linear points.


title: Single Micro-hole per Pixel for Thin Ge-on-Si Image Sensor with Enhanced Sensitivity upto 1700 nm

abstract: We present a Ge-on_Si CMOS image sensor with backside illumination for the near-infrared electromagnetic waves, wavelengths range 300-1700nm, detection essential for optical sensor technology. The micro-holes help to enhance the optical efficiency and extend the range to the 1.7 microns wavelength. We demonstrate an optimization for the width and depth of the nano-holes for maximal absorption in the near infrared. We show a reduction in cross-talk by employing thin Si oxide deep trench isolation in between the pixels. Finally, we show a 26-50 percent reduction in the device capacitance with the introduction of a hole. Such CMOS-compatible Ge-on_Si sensor will enable high-density, ultra-fast and efficient near-infrared imaging.


title: A Closer Look at Evaluating the Bit-Flip Attack Against Deep Neural Networks

abstract: Deep neural network models are massively deployed on a wide variety of hardware platforms. This results in the appearance of new attack vectors that significantly extend the standard attack surface, extensively studied by the adversarial machine learning community. One of the first attack that aims at drastically dropping the performance of a model, by targeting its parameters (weights) stored in memory, is the Bit-Flip Attack (BFA). In this work, we point out several evaluation challenges related to the BFA. First of all, the lack of an adversary's budget in the standard threat model is problematic, especially when dealing with physical attacks. Moreover, since the BFA presents critical variability, we discuss the influence of some training parameters and the importance of the model architecture. This work is the first to present the impact of the BFA against fully-connected architectures that present different behaviors compared to convolutional neural networks. These results highlight the importance of defining robust and sound evaluation methodologies to properly evaluate the dangers of parameter-based attacks as well as measure the real level of robustness offered by a defense.


title: Lotka-Volterra Models for Extraterrestrial Self-Replicating Probes

abstract: A sufficiently advanced extraterrestrial civilization can send out a swarm of self-replicating probes for space exploration. Given the fast-growing number of such a probe, even if there is only one extraterrestrial civilization sending out such probes in the Milky Way galaxy, we should still expect to see them. The fact that we do not consists part of the Fermi paradox. The suggestion that self-replicating probes will eventually mutate to consume their progenitors and therefore significantly reduce the number of total probes has been investigated and dismissed in the literature. In this work, we re-visit this question with a more realistic Lotka-Volterra model, and show that mutated probes would drive the progenitor probes into "extinction", thereby replacing them to spread throughout the galaxy. Thus, the efficiency of mutated probes in reducing the total number of self-replicating probes is even less than previously thought. As part of the analysis, we also suggest that, somewhat counter-intuitively, in designing self-replicating probes, one should not program them to stop replicating when sufficient mutation causes the probes to fail to recognize the progenitor probes as "self".


title: Framework for Highway Traffic Profiling using Connected Vehicle Data

abstract: The connected vehicle (CV) data could potentially revolutionize the traffic monitoring landscape as a new source of CV data that are collected exclusively from original equipment manufactures (OEMs) have emerged in the commercial market in recent years. Compared to existing CV data that are used by agencies, the new-generation of CV data have certain advantages including nearly ubiquitous coverage, high temporal resolution, high spatial accuracy, and enriched vehicle telematics data (e.g., hard braking events). This paper proposed a traffic profiling framework that target vehicle-level performance indexes across mobility, safety, riding comfort, traffic flow stability, and fuel consumption. The proof-of-concept study of a major interstate highway (i.e., I-280 NJ), using the CV data, illustrates the feasibility of going beyond traditional aggregated traffic metrics. Lastly, potential applications for either historical analysis and even near real-time monitoring are discussed. The proposed framework can be easily scaled and is particularly valuable for agencies that wish to systemically monitoring regional or statewide roadways without substantial investment on infrastructure-based sensing (and the associated on-going maintenance costs)


title: The Type II Dirac Seesaw Portal to the mirror sector: Connecting neutrino masses and a solution to the strong CP problem

abstract: We present a version of the Type II Seesaw mechanism for parametrically small Dirac neutrino masses. Our model starts from an $\text{SU}(2)_\text{L} \otimes \text{SU}(2)'\otimes \text{U}(1)_\text{X}$ gauge extension of the Standard Model involving a sector of mirror fermions. A bidoublet scalar with a very small vacuum expectation value connects the SM leptons with their mirror counterparts and we can identify the mirror neutrino with the right-handed neutrino. Similar to the conventional Type II Seesaw, our particle spectrum features singly- and doubly-charged scalars. The strong CP problem is solved by a discrete exchange symmetry between the two sectors that forces the contributions of quarks and mirror quarks to the strong CP phase to cancel each other. We discuss the low-energy phenomenology, comment on the cosmological implications of this scenario and indicate how to realize successful Dirac leptogenesis.


title: Eigenvalue avoidance of structured matrices depending smoothly on a real parameter

abstract: We explore the concept of eigenvalue avoidance, which is well understood for real symmetric and Hermitian matrices, for other classes of structured matrices. We adopt a differential geometric perspective and study the generic behaviour of the eigenvalues of regular and injective curves $t \in ]a,b[ \mapsto A(t) \in \mathcal{N} $ where $\mathcal{N}$ is a smooth real Riemannian submanifold of either $\mathbb{R}^{n \times n}$ or $\mathbb{C}^{n \times n}$. We focus on the case where $\mathcal{N}$ corresponds to some class of (real or complex) structured matrices including skew-symmetric, skew-Hermitian, orthogonal, unitary, banded symmetric, banded Hermitian, banded skew-symmetric, and banded skew-Hermitian. We argue that for some structures eigenvalue avoidance always happens, whereas for other structures this may depend on the parity of the size, on the numerical value of the multiple eigenvalue, and possibly on the value of determinant. As a further application of our tools we also study singular value avoidance for unstructured matrices in $\mathbb{R}^{m \times n}$ or $\mathbb{C}^{m \times n}$.


title: Electronic Structure Calculation with the Exact Pseudopotential and Interpolating Wavelet Basis

abstract: Electronic structure calculations are mostly carried out with Coulomb potential singularity adapted basis sets like STO or contracted GTO. With other basis or for heavy elements the pseudopotentials may appear as a practical alternative. Here, we introduce the exact pseudopotential (EPP) to remove the Coulomb singularity and test it for orbitals of small atoms with the interpolating wave basis set. We apply EPP to the Galerkin method with a basis set consisting of Deslauriers--Dubuc scaling functions on the half-infinite real interval. We demonstrate the EPP--Galerkin method by computing the hydrogen atom 1s, 2s, and 2p orbitals and helium atom configurations $\mathrm{He\;1s^2}$, $\mathrm{He\;1s2s\;{}^1 S}$, and $\mathrm{He\;1s2s\;{}^3 S}$. We compare the method to the ordinary interpolating wavelet Galerkin method (OIW--Galerkin) handling the singularity at the nucleus by excluding the scaling function located at the origin from the basis. We also compare the performance of our approach to that of finite--difference approach, which is another practical method for spherical atoms. We find the accuracy of the EPP--Galerkin method better than both of the above mentioned methods.


title: Compositional Score Modeling for Simulation-based Inference

abstract: Neural Posterior Estimation methods for simulation-based inference can be ill-suited for dealing with posterior distributions obtained by conditioning on multiple observations, as they tend to require a large number of simulator calls to learn accurate approximations. In contrast, Neural Likelihood Estimation methods can handle multiple observations at inference time after learning from individual observations, but they rely on standard inference methods, such as MCMC or variational inference, which come with certain performance drawbacks. We introduce a new method based on conditional score modeling that enjoys the benefits of both approaches. We model the scores of the (diffused) posterior distributions induced by individual observations, and introduce a way of combining the learned scores to approximately sample from the target posterior distribution. Our approach is sample-efficient, can naturally aggregate multiple observations at inference time, and avoids the drawbacks of standard inference methods.


title: B2B Advertising: Joint Dynamic Scoring of Account and Users

abstract: When a business sells to another business (B2B), the buying business is represented by a group of individuals, termed account, who collectively decide whether to buy. The seller advertises to each individual and interacts with them, mostly by digital means. The sales cycle is long, most often over a few months. There is heterogeneity among individuals belonging to an account in seeking information and hence the seller needs to score the interest of each individual over a long horizon to decide which individuals must be reached and when. Moreover, the buy decision rests with the account and must be scored to project the likelihood of purchase, a decision that is subject to change all the way up to the actual decision, emblematic of group decision making. We score decision of the account and its individuals in a dynamic manner. Dynamic scoring allows opportunity to influence different individual members at different time points over the long horizon. The dataset contains behavior logs of each individual's communication activities with the seller; but, there are no data on consultations among individuals which result in the decision. Using neural network architecture, we propose several ways to aggregate information from individual members' activities, to predict the group's collective decision. Multiple evaluations find strong model performance.


title: The Y-Product

abstract: We present a topological construction that provides many examples of non-commutative Frobenius algebras that generalizes the well-known pair-of-pants. When applied to the solid torus, in conjunction with Crane-Yetter theory, we provide a topological proof of the Verlinde formula. We also apply the construction to a solid handlebody of higher genus, leading to a generalization of the Verlinde formula (not the higher genus Verlinde formula); in particular, we define a generalized $S$-matrix. Finally, we discuss the relation between our construction and Yetter's construction of a handle as a Hopf algebra, and give a generalization.


title: Physics-aware Differentiable Discrete Codesign for Diffractive Optical Neural Networks

abstract: Diffractive optical neural networks (DONNs) have attracted lots of attention as they bring significant advantages in terms of power efficiency, parallelism, and computational speed compared with conventional deep neural networks (DNNs), which have intrinsic limitations when implemented on digital platforms. However, inversely mapping algorithm-trained physical model parameters onto real-world optical devices with discrete values is a non-trivial task as existing optical devices have non-unified discrete levels and non-monotonic properties. This work proposes a novel device-to-system hardware-software codesign framework, which enables efficient physics-aware training of DONNs w.r.t arbitrary experimental measured optical devices across layers. Specifically, Gumbel-Softmax is employed to enable differentiable discrete mapping from real-world device parameters into the forward function of DONNs, where the physical parameters in DONNs can be trained by simply minimizing the loss function of the ML task. The results have demonstrated that our proposed framework offers significant advantages over conventional quantization-based methods, especially with low-precision optical devices. Finally, the proposed algorithm is fully verified with physical experimental optical systems in low-precision settings.


title: Universality in the tripartite information after global quenches

abstract: We consider macroscopically large 3-partitions $(A,B,C)$ of connected subsystems $A\cup B \cup C$ in infinite quantum spin chains and study the R\'enyi-$\alpha$ tripartite information $I_3^{(\alpha)}(A,B,C)$. At equilibrium in clean 1D systems with local Hamiltonians it generally vanishes. A notable exception is the ground state of conformal critical systems, in which $I_3^{(\alpha)}(A,B,C)$ is known to be a universal function of the cross ratio $x=|A||C|/[(|A|+|B|)(|C|+|B|)]$, where $|A|$ denotes $A$'s length. We identify different classes of states that, under time evolution with translationally invariant Hamiltonians, locally relax to states with a nonzero (R\'enyi) tripartite information, which furthermore exhibits a universal dependency on $x$. We report a numerical study of $I_3^{(\alpha)}$ in systems that are dual to free fermions, propose a field-theory description, and work out their asymptotic behaviour for $\alpha=2$ in general and for generic $\alpha$ in a subclass of systems. This allows us to infer the value of $I_3^{(\alpha)}$ in the scaling limit $x\rightarrow 1^-$, which we call ``residual tripartite information''. If nonzero, our analysis points to a universal residual value $-\log 2$ independently of the R\'enyi index $\alpha$, and hence applies also to the genuine (von Neumann) tripartite information.


title: Preempting to Minimize Age of Incorrect Information under Transmission Delay

abstract: We study the problem of optimizing the decisions of a preemptively capable transmitter to minimize the Age of Incorrect Information (AoII) when the communication channel has a random delay. We consider a slotted-time system where a transmitter observes a Markovian source and makes decisions based on the system status. In each time slot, the transmitter decides whether to preempt or skip when the channel is busy. When the channel is idle, the transmitter decides whether to send a new update. A remote receiver estimates the state of the Markovian source based on the update it receives. We consider a generic transmission delay and assume that the transmission delay is independent and identically distributed for each update. This paper aims to optimize the transmitter's decision in each time slot to minimize the AoII with generic time penalty functions. To this end, we first use the Markov decision process to formulate the optimization problem and derive the analytical expressions of the expected AoIIs achieved by two canonical preemptive policies. Then, we prove the existence of the optimal policy and provide a feasible value iteration algorithm to approximate the optimal policy. However, the value iteration algorithm will be computationally expensive if we want considerable confidence in the approximation. Therefore, we analyze the system characteristics under two canonical delay distributions and theoretically obtain the corresponding optimal policies using the policy improvement theorem. Finally, numerical results are presented to illustrate the performance improvements brought about by the preemption capability.


title: Discrete Mechanics and Optimal Control for Passive Walking with Foot Slippage

abstract: Forced variational integrators are given by the discretization of the Lagrange-d'Alembert principle for systems subject to external forces, and have proved useful for numerical simulation studies of complex dynamical systems. In this paper we model a passive walker with foot slip by using techniques of geometric mechanics, and we construct forced variational integrators for the system. Moreover, we present a methodology for generating (locally) optimal control policies for simple hybrid holonomically constrained forced Lagrangian systems, based on discrete mechanics, applied to a controlled walker with foot slip in a trajectory tracking problem.


title: Spiral instabilities: Mode saturation and decay

abstract: This paper continues a series reporting different aspects of the behaviour of disc galaxy simulations that support spiral instabilities. The focus in this paper is to demonstrate how linear spiral instabilities saturate and decay, and how the properties of the disc affect the limiting amplitude of the spirals. Once again, we employ idealized models that each possess a single instability that we follow until it has run its course. Remarkably, we find a tight correlation between the growth rate of the mode and its limiting amplitude, albeit from only six simulations. We show that non-linear orbit deflections near corotation cause the mode to saturate, and that the more time available in a slowly-growing mode creates the critical deflections at lower amplitude. We also find that scattering at the inner Lindblad resonance is insignificant until after the mode has saturated. Our objective in this series of papers, which we believe we have now achieved, has been to develop a convincing and well-documented account of the physical behaviour of the spiral patterns that have been observed in simulations by others, and by ourselves, for many decades. Understanding the simulations is an important step towards the greater objective, which is to find observational evidence from galaxies that could confront the identified mechanism.


title: Are Am stars and hot-Jupiter planets related?

abstract: Am stars are often components of short-period binary systems, where tidal interactions would result in low rotational velocities and help to develop the chemical peculiarities observed. However, the origin of single Am stars and Am stars that belong to wide binary systems is unclear. There is very recent evidence of an Am star hosting a hot-brown dwarf likely synchronized and other possible Am stars hosting hot-Jupiter planets. We wonder if these hot-low mass companions could play a role in the development of an Am star, that is to say, if they could help to mitigate the "single Am" problem. We studied a sample of 19 early-type stars, 7 of them hosting hot-brown dwarfs and 12 of them hosting hot-Jupiter planets. We detected 4 Am stars in our sample (KELT-19A, KELT-17, HATS-70 and TOI-503) and 2 possible Am stars (TOI-681 and HAT-P-69). In particular, we detected the new Am star HATS-70 which hosts a hot-brown dwarf, and rule out this class for the hot-Jupiter host WASP-189, both showing different composition than previously reported. We estimated the incidence of Am stars within stars hosting hot-brown dwarfs (50-75%) and within stars hosting hot-Jupiters (20-42%). The incidence of Am stars hosting hot-brown dwarfs resulted higher than the frequency of Am stars in general. This would imply that the presence of hot-brown dwarfs could play a role in the development of Am stars and possibly help to mitigate the "single Am" problem, different to the case of hot-Jupiter planets. Notably, these results would also indicate that the search for hot-brown dwarfs may be benefited by targeting single Am stars or Am stars in wide binary systems. We encourage the analysis off additional early-type stars hosting hot-companions in order to improve the significance of the initial trends found here. [abridged]


title: Helly-type theorems for the ordering of the vertices of a hypergraph

abstract: Let $H$ be a complete $r$-uniform hypergraph such that two vertices are marked in each edge as its `boundary' vertices. A linear ordering of the vertex set of $H$ is called an {\em agreeing linear order}, provided all vertices of each edge of $H$ lie between its two boundary vertices. We prove the following Helly-type theorem: if there is an {agreeing linear order} on the vertex set of every subhypergraph of $H$ with at most $2r-2$ vertices, then there is an agreeing linear order on the vertex set of $H$. We also show that the constant $2r-2$ cannot be reduced in the theorem. The case $r=3$ of the theorem has particular interest in the axiomatic theory of betweenness. Similar results are obtained for further $r$-uniform hypergraphs ($r\geq 3$), where one or two vertices are marked in each edge, and the linear orders need to satisfy various rules of agreement. In one of the cases we prove that no such Helly-type statement holds.