title: The FluidFlower International Benchmark Study: Process, Modeling Results, and Comparison to Experimental Data

abstract: Successful deployment of geological carbon storage (GCS) requires an extensive use of reservoir simulators for screening, ranking and optimization of storage sites. However, the time scales of GCS are such that no sufficient long-term data is available yet to validate the simulators against. As a consequence, there is currently no solid basis for assessing the quality with which the dynamics of large-scale GCS operations can be forecasted. To meet this knowledge gap, we have conducted a major GCS validation benchmark study. To achieve reasonable time scales, a laboratory-size geological storage formation was constructed (the "FluidFlower"), forming the basis for both the experimental and computational work. A validation experiment consisting of repeated GCS operations was conducted in the FluidFlower, providing what we define as the true physical dynamics for this system. Nine different research groups from around the world provided forecasts, both individually and collaboratively, based on a detailed physical and petrophysical characterization of the FluidFlower sands. The major contribution of this paper is a report and discussion of the results of the validation benchmark study, complemented by a description of the benchmarking process and the participating computational models. The forecasts from the participating groups are compared to each other and to the experimental data by means of various indicative qualitative and quantitative measures. By this, we provide a detailed assessment of the capabilities of reservoir simulators and their users to capture both the injection and post-injection dynamics of the GCS operations.


title: Towards a responsible machine learning approach to identify forced labor in fisheries

abstract: Many fishing vessels use forced labor, but identifying vessels that engage in this practice is challenging because few are regularly inspected. We developed a positive-unlabeled learning algorithm using vessel characteristics and movement patterns to estimate an upper bound of the number of positive cases of forced labor, with the goal of helping make accurate, responsible, and fair decisions. 89% of the reported cases of forced labor were correctly classified as positive (recall) while 98% of the vessels certified as having decent working conditions were correctly classified as negative. The recall was high for vessels from different regions using different gears, except for trawlers. We found that as much as ~28% of vessels may operate using forced labor, with the fraction much higher in squid jiggers and longlines. This model could inform risk-based port inspections as part of a broader monitoring, control, and surveillance regime to reduce forced labor. * Translated versions of the English title and abstract are available in five languages in S1 Text: Spanish, French, Simplified Chinese, Traditional Chinese, and Indonesian.


title: Travel times and ray paths for acoustic and elastic waves in generally anisotropic media

abstract: Wavefield travel time tomography is used for a variety of purposes in acoustics, geophysics and non-destructive testing. Since the problem is non-linear, assessing uncertainty in the results requires many forward evaluations. It is therefore important that the forward evaluation of travel times and ray paths is efficient, which is challenging in anisotropic media. Given a computed travel time field, ray tracing can be performed to obtain the fastest ray path from any point in the medium to the source of the travel time field. These rays can then be used to speed up gradient based inversion methods. We present a forward modeller for calculating travel time fields by localised estimation of wavefronts, and a novel approach to ray tracing through travel time fields. These methods have been tested in a complex anisotropic weld and give travel times comparable to those obtained using finite element modelling while being computationally cheaper.


title: On the weak-gravity bound for a shift-symmetric scalar field

abstract: The weak-gravity bound has been discovered in several asymptotically safe gravity-matter systems. It limits the strength of gravitational fluctuations that are compatible with an ultraviolet-complete matter sector, and results from the collision of two partial fixed points of the matter system as a function of the strength of the gravitational interactions. In this paper, we will investigate this mechanism in detail for a shift-symmetric scalar field. First, we will study the fixed point structure of the scalar system without gravity. We find indications that the Gaussian fixed point is the only viable fixed point, suggesting that a weak-gravity bound resulting from the collision of two partial fixed points is a truncation artefact. We will then couple the scalar system to gravity and perform different expansions to track the Gaussian fixed point as gravitational fluctuations become stronger. We also introduce a new notion of the weak-gravity bound that is based on the number of relevant operators.


title: A characterization of Rieffel's deformed algebra as Heisenberg smooth operators

abstract: Let $\mathcal{C}$ be a unital C$^*$-algebra and $E_n$ be the Hilbert $\mathcal{C}$-module defined as the completion of the $\mathcal{C}$-valued Schwartz function space $\mathcal{S}^\mathcal{C}(\mathbb{R}^n)$ with respect to the norm $\|f\|_2 := \left\| \int_{\mathbb{R}^n} f(x)^*f(x) \, dx \right\|_\mathcal{C}^{1 / 2}$. Also, let $\text{Ad }\mathcal{U}$ be the canonical action of the $(2n + 1)$-dimensional Heisenberg group by conjugation on the algebra of adjointable operators on $E_n$ and let $J$ be a skew-symmetric linear transformation on $\mathbb{R}^n$. We characterize the smooth vectors under $\text{Ad }\mathcal{U}$ which commute with a certain algebra of right multiplication operators $R_h$, with $h \in \mathcal{S}^\mathcal{C}(\mathbb{R}^n)$, where the product is "twisted" with respect to $J$ according to a deformation quantization procedure introduced by M.A. Rieffel. More precisely, we establish that they coincide with the corresponding algebra of left multiplication operators, as conjectured by Rieffel.


title: Enabling Versatile Privacy Interfaces Using Machine-Readable Transparency Information

abstract: Transparency regarding the processing of personal data in online services is a necessary precondition for informed decisions on whether or not to share personal data. In this paper, we argue that privacy interfaces shall incorporate the context of display, personal preferences, and individual competences of data subjects following the principles of universal design and usable privacy. Doing so requires -- among others -- to consciously decouple the provision of transparency information from their ultimate presentation. To this end, we provide a general model of how transparency information can be provided from a data controller to data subjects, effectively leveraging machine-readable transparency information and facilitating versatile presentation interfaces. We contribute two actual implementations of said model: 1) a GDPR-aligned privacy dashboard and 2) a chatbot and virtual voice assistant enabled by conversational AI. We evaluate our model and implementations with a user study and find that these approaches provide effective and time-efficient transparency. Consequently, we illustrate how transparency can be enhanced using machine-readable transparency information and how data controllers can meet respective regulatory obligations.


title: A model of anisotropic winds from rotating stars for evolutionary calculations

abstract: Context: The surface properties of rotating stars can vary from pole to equator, resulting in anisotropic stellar winds which are not included in the currently available evolutionary models. Aims: We develop a formalism to describe the mass and angular momentum loss of rotating stars which takes into account both the varying surface properties and distortion due to rotation. Methods: Adopting the mass-loss recipe for non-rotating stars, we assigned to each point on the surface of a rotating star an equivalent non-rotating star, for which the surface mass flux is given by the recipe. The global mass-loss and angular momentum loss rates are then given by integrating over the deformed stellar surface as appropriate. Evolutionary models were computed and our prescription is compared to the currently used simple mass-loss enhancement recipes for rotating stars. Results: We find that mass-loss rates are largely insensitive to rotation for models not affected by the bi-stability jump. For those affected by the bi-stability jump, the increase in mass-loss rates with respect to time is smoothed. As our prescription considers the variation of physical conditions over the stellar surface, the region affected by the bi-stability jump is able to grow gradually instead of the whole star suddenly being affected. Conclusion: We have provided an easy to implement and flexible, yet physically meaningful prescription for calculating mass and angular momentum loss rates of rotating stars in a one-dimensional stellar evolution code which compares favourably to more physically comprehensive models. The implementation of our scheme in the stellar evolution code MESA is available online: https://zenodo.org/record/7437006


title: A convergent finite-volume scheme for nonlocal cross-diffusion systems for multi-species populations

abstract: An implicit Euler finite-volume scheme for a nonlocal cross-diffusion system on the one-dimensional torus, arising in population dynamics, is proposed and analyzed. The kernels are assumed to be in detailed balance and satisfy a weak cross-diffusion condition. The latter condition allows for negative off-diagonal coefficients and for kernels defined by an indicator function. The scheme preserves the nonnegativity of the densities, conservation of mass, and production of the Boltzmann and Rao entropies. The key idea is to ``translate'' the entropy calculations for the continuous equations to the finite-volume scheme, in particular to design discretizations of the mobilities, which guarantee a discrete chain rule even in the presence of nonlocal terms. Based on this idea, the existence of finite-volume solutions and the convergence of the scheme are proven. As a by-product, we deduce the existence of weak solutions to the continuous cross-diffusion system. Finally, we present some numerical experiments illustrating the behavior of the solutions to the nonlocal and associated local models.


title: Impact of artificial topological changes on flow and transport through fractured media due to mesh resolution

abstract: We performed a set of numerical simulations to characterize the interplay of fracture network topology, upscaling, and mesh refinement on flow and transport properties in fractured porous media. We generated a set of generic three-dimensional discrete fracture networks at various densities, where the radii of the fractures were sampled from a truncated power-law distribution, and whose parameters were loosely based on field site characterizations. We also considered five network densities, which were defined using a dimensionless version of density based on percolation theory. Once the networks were generated, we upscaled them into a single continuum model using the upscaled discrete fracture matrix model presented by Sweeney et al. We considered steady, isothermal pressure-driven flow through each domain and then simulated conservative, decaying, and adsorbing tracers using a pulse injection into the domain. For each simulation, we calculated the effective permeability and solute breakthrough curves as quantities of interest to compare between network realizations. We found that selecting a mesh resolution such that the global topology of the upscaled mesh matches the fracture network is essential. If the upscaled mesh has a connected pathway of fracture (higher permeability) cells but the fracture network does not, then the estimates for effective permeability and solute breakthrough will be incorrect. False connections cannot be eliminated entirely, but they can be managed by choosing appropriate mesh resolution and refinement for a given network. Adopting octree meshing to obtain sufficient levels of refinement leads to fewer computational cells (up to a 90% reduction in overall cell count) when compared to using a uniform resolution grid and can result in a more accurate continuum representation of the true fracture network.

categories: cs.CE cs.NA math.NA physics.comp-ph physics.flu-dyn


title: Vandermonde Trajectory Bounds for Linear Companion Systems

abstract: Fast and accurate safety assessment and collision checking are essential for motion planning and control of highly dynamic autonomous robotic systems. Informative, intuitive, and explicit motion trajectory bounds enable explainable and time-critical safety verification of autonomous robot motion. In this paper, we consider feedback linearization of nonlinear systems in the form of proportional-and-higher-order-derivative (PhD) control corresponding to companion dynamics. We introduce a novel analytic convex trajectory bound, called $\textit{Vandermonde simplex}$, for high-order companion systems, that is given by the convex hull of a finite weighted combination of system position, velocity, and other relevant higher-order state variables. Our construction of Vandermonde simplexes is built based on expressing the solution trajectory of companion dynamics in a newly introduced family of $\textit{Vandermonde basis functions}$ that offer new insights for understanding companion system motion compared to the classical exponential basis functions. In numerical simulations, we demonstrate that Vandermonde simplexes offer significantly more accurate motion prediction (e.g., at least an order of magnitude improvement in estimated motion volume) for describing the motion trajectory of companion systems compared to the standard invariant Lyapunov ellipsoids as well as exponential simplexes built based on exponential basis functions.


title: A two-stage stochastic programming model for electric substation flood mitigation prior to an imminent hurricane

abstract: We present a stochastic programming model for informing the deployment of ad hoc flood mitigation measures to protect electrical substations prior to an imminent and uncertain hurricane. The first stage captures the deployment of a fixed number of mitigation resources, and the second stage captures grid operation in response to a contingency. The primary objective is to minimize expected load shed. We develop methods for simulating flooding induced by extreme rainfall and construct two geographically realistic case studies, one based on Tropical Storm Imelda and the other on Hurricane Harvey. Applying our model to those case studies, we investigate the effect of the mitigation budget on the optimal objective value and solutions. Our results highlight the sensitivity of the optimal mitigation to the budget, a consequence of those decisions being discrete. We additionally assess the value of having better mitigation options and the spatial features of the optimal mitigation.


title: Robust Auto-landing Control of an agile Regional Jet Using Fuzzy Q-learning

abstract: A robust auto-landing problem of a Truss-braced Wing (TBW) regional jet aircraft with poor stability characteristics is presented in this study employing a Fuzzy Reinforcement Learning scheme. Reinforcement Learning (RL) has seen a recent surge in practical uses in control systems. In contrast to many studies implementing Deep Learning in RL algorithms to generate continuous actions, the methodology of this study is straightforward and avoids complex neural network architectures by applying Fuzzy rules. An innovative, agile civil aircraft is selected not only to meet future aviation community expectations but also to demonstrate the robustness of the suggested method. In order to create a multi-objective RL environment, a Six-degree-of-freedom (6-DoF) simulation is first developed. By transforming the auto-landing problem of the aircraft into a Markov Decision Process (MDP) formulation, the problem is solved by designing a low-level Fuzzy Q-learning (FQL) controller. More specifically, the well-known Q-learning method, which is a discrete RL algorithm, is supplemented by Fuzzy rules to provide continuous actions with no need to complex learning structures. The performance of the proposed system is then evaluated by extensive flight simulations in different flight conditions considering severe wind gusts, measurement noises, actuator faults, and model uncertainties. Besides, the controller effectiveness would be compared with existing competing techniques such as Dynamic Inversion (DI) and Q-learning. The simulation results indicate the superior performance of the proposed control system as a reliable and robust control method to be employed in real applications.


title: On the LS-category of homomorphisms of groups with torsion

abstract: We prove the equality $\cat(\phi)=\cd(\phi)$ for homomorphisms $\phi:\Gamma\to \Lambda$ between finitely generated abelian groups $\Gamma$ and $\Lambda$, where $\phi(T(\Gamma))=0$ for the torsion subgroups $T(\Gamma)$ of $\Gamma$.


title: On validity of the quasi-static approximation in scalar-tensor theories

abstract: The discovery of cosmic acceleration motivated extensive studies of dynamical dark energy and modified gravity models. Of particular interest are the scalar-tensor theories, with a scalar field dark energy non-minimally coupled to matter. Cosmological constraints on these models often employ the quasi-static approximation (QSA), in which the dynamics of the scalar field perturbations is proportional to the perturbation in the matter density. Using the QSA simplifies the physical interpretation of the phenomenology of scalar-tensor theories, and results in substantial savings of computing time when deriving parameter constraints. Focusing on the symmetron model, which is a well-motivated scalar-tensor theory with a screening mechanism, we compare the exact solution of the linearly perturbed field equations to those obtained under the QSA and identify the range of the model parameters for which the QSA is valid. We find that the evolution of background scalar field is most important, namely, whether it is dominated by the Hubble friction or the scalar field potential. This helps us derive a criterion for the symmetron model, but same argument can be applied to other scalar-tensor theories of generalized Brans-Dicke type. We consider two scenarios, one where the scalar field is only coupled to dark matter and where it couples to all of the matter.


title: CHA2: CHemistry Aware Convex Hull Autoencoder Towards Inverse Molecular Design

abstract: Optimizing molecular design and discovering novel chemical structures to meet certain objectives, such as quantitative estimates of the drug-likeness score (QEDs), is NP-hard due to the vast combinatorial design space of discrete molecular structures, which makes it near impossible to explore the entire search space comprehensively to exploit de novo structures with properties of interest. To address this challenge, reducing the intractable search space into a lower-dimensional latent volume helps examine molecular candidates more feasibly via inverse design. Autoencoders are suitable deep learning techniques, equipped with an encoder that reduces the discrete molecular structure into a latent space and a decoder that inverts the search space back to the molecular design. The continuous property of the latent space, which characterizes the discrete chemical structures, provides a flexible representation for inverse design in order to discover novel molecules. However, exploring this latent space requires certain insights to generate new structures. We propose using a convex hall surrounding the top molecules in terms of high QEDs to ensnare a tight subspace in the latent representation as an efficient way to reveal novel molecules with high QEDs. We demonstrate the effectiveness of our suggested method by using the QM9 as a training dataset along with the Self- Referencing Embedded Strings (SELFIES) representation to calibrate the autoencoder in order to carry out the Inverse molecular design that leads to unfold novel chemical structure.


title: Tensor enriched categorical generalization of the Eilenberg-Watts theorem

abstract: Given monoids $\mathfrak{b}$, $\mathfrak{b}'$ in a B\'{e}nabou cosmos, we prove that the category of cocontinuous enriched functors between enriched categories of right modules is equivalent to the category of $(\mathfrak{b},\mathfrak{b}')$-bimodules. When the given monoids $\mathfrak{b}$, $\mathfrak{b}'$ are commutative, we further show that the category of cocontinuous lax tensor enriched functors between tensor enriched categories of right modules is equivalent to the category of commutative monoids over $\mathfrak{b}\otimes\mathfrak{b}'$.


title: Learning Physical Models that Can Respect Conservation Laws

abstract: Recent work in scientific machine learning (SciML) has focused on incorporating partial differential equation (PDE) information into the learning process. Much of this work has focused on relatively "easy" PDE operators (e.g., elliptic and parabolic), with less emphasis on relatively "hard" PDE operators (e.g., hyperbolic). Within numerical PDEs, the latter problem class requires control of a type of volume element or conservation constraint, which is known to be challenging. Delivering on the promise of SciML requires seamlessly incorporating both types of problems into the learning process. To address this issue, we propose ProbConserv, a framework for incorporating conservation constraints into a generic SciML architecture. To do so, ProbConserv combines the integral form of a conservation law with a Bayesian update. We provide a detailed analysis of ProbConserv on learning with the Generalized Porous Medium Equation (GPME), a widely-applicable parameterized family of PDEs that illustrates the qualitative properties of both easier and harder PDEs. ProbConserv is effective for easy GPME variants, performing well with state-of-the-art competitors; and for harder GPME variants it outperforms other approaches that do not guarantee volume conservation. ProbConserv seamlessly enforces physical conservation constraints, maintains probabilistic uncertainty quantification (UQ), and deals well with shocks and heteroscedasticities. In each case, it achieves superior predictive performance on downstream tasks.


title: Chow-Witt rings and topology of flag varieties

abstract: The paper computes the Witt-sheaf cohomology rings of partial flag varieties in type A in terms of the Pontryagin classes of the subquotient bundles. The proof is based on a Leray-Hirsch-type theorem for Witt-sheaf cohomology for the maximal rank cases, and a detailed study of cohomology ring presentations and annihilators of characteristic classes for the general case. The computations have consequences for the topology of real flag manifolds: we show that all torsion in the integral cohomology is 2-torsion, which was not known in full generality previously. This allows for example to compute the Poincar\'e polynomials of complete flag varieties for cohomology with twisted integer coefficients. The computations also allow to describe the Chow-Witt rings of flag varieties, and we sketch an enumerative application to counting flags satisfying multiple incidence conditions to given hypersurfaces.


title: On the Convergence of an IEQ-based first-order Numerical Scheme for the Beris-Edwards System

abstract: We present a convergence analysis of an unconditionally energy-stable first-order semi-discrete numerical scheme designed for a hydrodynamic Q-tensor model, the so-called Beris-Edwards system, based on the Invariant Energy Quadratization Method (IEQ). The model consists of the Navier-Stokes equations for the fluid flow, coupled to the Q-tensor gradient flow describing the liquid crystal molecule alignment. By using the Invariant Energy Quadratization Method, we obtain a linearly implicit scheme, accelerating the computational speed. However, this introduces an auxiliary variable to replace the bulk potential energy and it is a priori unclear whether the reformulated system is equivalent to the Beris-Edward system. In this work, we prove stability properties of the scheme and show its convergence to a weak solution of the coupled liquid crystal system. We also demonstrate the equivalence of the reformulated and original systems in the weak sense.


title: A topological space associated to corank 1 tropical phased matroids

abstract: A consequence of the Folkman-Lawrence topological representation theorem is that the geometric realization of the order complex of the poset of non-zero covectors of a loopless rank $n-1$ oriented matroid on $[n]$ is homeomorphic to an $(n-2)$-sphere. In this paper, we begin the study of an analogous theorem for tropical phased matroids by proving that the topological order complex for a loopless rank $n-1$ tropical phased matroid on $[n]$ is homeomorphic to a $(2n-3)$-sphere.


title: Data-driven reduced-order modelling for blood flow simulations with geometry-informed snapshots

abstract: Parametric reduced-order modelling often serves as a surrogate method for hemodynamics simulations to improve the computational efficiency in many-query scenarios or to perform real-time simulations. However, the snapshots of the method require to be collected from the same discretisation, which is a straightforward process for physical parameters, but becomes challenging for geometrical problems, especially for those domains featuring unparameterised and unique shapes, e.g. patient-specific geometries. In this work, a data-driven surrogate model is proposed for the efficient prediction of blood flow simulations on similar but distinct domains. The proposed surrogate model leverages group surface registration to parameterise those shapes and formulates corresponding hemodynamics information into geometry-informed snapshots by the diffeomorphisms constructed between a reference domain and original domains. A non-intrusive reduced-order model for geometrical parameters is subsequently constructed using proper orthogonal decomposition, and a radial basis function interpolator is trained for predicting the reduced coefficients of the reduced-order model based on compressed geometrical parameters of the shape. Two examples of blood flowing through a stenosis and a bifurcation are presented and analysed. The proposed surrogate model demonstrates its accuracy and efficiency in hemodynamics prediction and shows its potential application toward real-time simulation or uncertainty quantification for complex patient-specific scenarios.

categories: cs.CE cs.LG physics.bio-ph physics.flu-dyn


title: Unification of popular artificial neural network activation functions

abstract: We present a unified representation of the most popular neural network activation functions. Adopting Mittag-Leffler functions of fractional calculus, we propose a flexible and compact functional form that is able to interpolate between various activation functions and mitigate common problems in training neural networks such as vanishing and exploding gradients. The presented gated representation extends the scope of fixed-shape activation functions to their adaptive counterparts whose shape can be learnt from the training data. The derivatives of the proposed functional form can also be expressed in terms of Mittag-Leffler functions making it a suitable candidate for gradient-based backpropagation algorithms. By training multiple neural networks of different complexities on various datasets with different sizes, we demonstrate that adopting a unified gated representation of activation functions offers a promising and affordable alternative to individual built-in implementations of activation functions in conventional machine learning frameworks.


title: Model adaptation for hyperbolic balance laws

abstract: In this work, we devise a model adaptation strategy for a class of model hierarchies consisting of two levels of model complexity. In particular, the fine model consists of a system of hyperbolic balance laws with stiff reaction terms and the coarse model consists of a system of hyperbolic conservation laws. We employ the relative entropy stability framework to obtain an a posteriori modeling error estimator. The efficiency of the model adaptation strategy is demonstrated by conducting simulations for chemically reacting fluid mixtures in one space dimension.


title: Challenges in the Design and Implementation of IoT Testbeds in Smart-Cities: A Systematic Review

abstract: Advancements in wireless communication and the increased accessibility to low-cost sensing and data processing IoT technologies have increased the research and development of urban monitoring systems. Most smart city research projects rely on deploying proprietary IoT testbeds for indoor and outdoor data collection. Such testbeds typically rely on a three-tier architecture composed of the Endpoint, the Edge, and the Cloud. Managing the system's operation whilst considering the security and privacy challenges that emerge, such as data privacy controls, network security, and security updates on the devices, is challenging. This work presents a systematic study of the challenges of developing, deploying and managing urban monitoring testbeds, as experienced in a series of urban monitoring research projects, followed by an analysis of the relevant literature. By identifying the challenges in the various projects and organising them under the V-model development lifecycle levels, we provide a reference guide for future projects. Understanding the challenges early on will facilitate current and future smart-cities IoT research projects to reduce implementation time and deliver secure and resilient testbeds.


title: Formality in the Deligne-Langlands correspondence

abstract: The Deligne-Langlands correspondence parametrizes irreducible representations of the affine Hecke algebra $\mathcal{H}^{\text{aff}}$ by certain perverse sheaves. We show that this can be lifted to an equivalence of triangulated categories. More precisely, we construct for each central character $\chi$ of $\mathcal{H}^{\text{aff}}$ an equivalence of triangulated categories between a perfect derived category of dg-modules $D_{\text{perf}}(\mathcal{H}^{\text{aff}}/(\text{ker}(\chi)) - \text{dgMod})$ and the triangulated category generated by the corresponding perverse sheaves. The main step in this construction is a formality result that we prove for a wide range of `Springer sheaves'.


title: On an Instance of the Small Cohen-Macaulay Conjecture

abstract: We provide a simplified proof of a theorem proved by Tavanfar and Shimomoto which states that a quasi-Gorenstein deformation of a $3$-dimensional quasi-Gorenstein local ring $(A,m,k)$ with $H^2_m(A)=k$ admits a small Cohen-Macaulay module.


title: Likelihood Annealing: Fast Calibrated Uncertainty for Regression

abstract: Recent advances in deep learning have shown that uncertainty estimation is becoming increasingly important in applications such as medical imaging, natural language processing, and autonomous systems. However, accurately quantifying uncertainty remains a challenging problem, especially in regression tasks where the output space is continuous. Deep learning approaches that allow uncertainty estimation for regression problems often converge slowly and yield poorly calibrated uncertainty estimates that can not be effectively used for quantification. Recently proposed post hoc calibration techniques are seldom applicable to regression problems and often add overhead to an already slow model training phase. This work presents a fast calibrated uncertainty estimation method for regression tasks called Likelihood Annealing, that consistently improves the convergence of deep regression models and yields calibrated uncertainty without any post hoc calibration phase. Unlike previous methods for calibrated uncertainty in regression that focus only on low-dimensional regression problems, our method works well on a broad spectrum of regression problems, including high-dimensional regression.Our empirical analysis shows that our approach is generalizable to various network architectures, including multilayer perceptrons, 1D/2D convolutional networks, and graph neural networks, on five vastly diverse tasks, i.e., chaotic particle trajectory denoising, physical property prediction of molecules using 3D atomistic representation, natural image super-resolution, and medical image translation using MRI.


title: The Complexity of Learning (Pseudo)random Dynamics of Black Holes and Other Chaotic Systems

abstract: It has been recently proposed that the naive semiclassical prediction of non-unitary black hole evaporation can be understood in the fundamental description of the black hole as a consequence of ignorance of high-complexity information. Validity of this conjecture implies that any algorithm which is polynomially bounded in computational complexity cannot accurately reconstruct the black hole dynamics. In this work, we prove that such bounded quantum algorithms cannot accurately predict (pseudo)random unitary dynamics, even if they are given access to an arbitrary set of polynomially complex observables under this time evolution; this shows that "learning" a (pseudo)random unitary is computationally hard. We use the common simplification of modeling black holes and more generally chaotic systems via (pseudo)random dynamics. The quantum algorithms that we consider are completely general, and their attempted guess for the time evolution of black holes is likewise unconstrained: it need not be a linear operator, and may be as general as an arbitrary (e.g. decohering) quantum channel.


title: Assessment of Reinforcement Learning for Macro Placement

abstract: We provide open, transparent implementation and assessment of Google Brain's deep reinforcement learning approach to macro placement and its Circuit Training (CT) implementation in GitHub. We implement in open source key "blackbox" elements of CT, and clarify discrepancies between CT and Nature paper. New testcases on open enablements are developed and released. We assess CT alongside multiple alternative macro placers, with all evaluation flows and related scripts public in GitHub. Our experiments also encompass academic mixed-size placement benchmarks, as well as ablation and stability studies. We comment on the impact of Nature and CT, as well as directions for future research.


title: Water Dynamics in Highly Concentrated Protein Systems -- Insight from Nuclear Magnetic Resonance Relaxometry

abstract: 1H spin-lattice relaxation experiments have been performed for water - Bovine Serum Al-bumin (BSA) mixtures including 20% wt and 40% wt of BSA. The experiments have been carried out in a frequency range encompassing three orders of magnitude, from 10 kHz to 10 MHz, versus temperature. The relaxation data have been thoroughly analyzed in terms of several relaxation models with the purpose to reveal the mechanisms of water motion. For this purpose four relaxation models have been used: the data have been decomposed into relaxation contributions expressed in terms of Lorentzian spectral densities, then three dimensional translation diffusion has been assumed, next two dimensional surface diffusion has been considered and eventually, a model of surface diffusion mediated by acts of adsorption to the surface has been employed. In this way it has been demonstrated that the last concept is most plausible. Parameters describing the dynamics in a quantita-tive manner have been determined and discussed.


title: Radiation Reaction for Non-Spinning Bodies at 4.5PN in the Effective Field Theory Approach

abstract: We calculate the 2 post-Newtonian correction to the radiation reaction acceleration for non-spinning binary systems, which amounts to the 4.5 post-Newtonian correction to Newtonian acceleration. The calculation is carried out completely using the effective field theory approach. The center-of-mass corrections to the results are complicated and are discussed in detail. Non-trivial consistency checks are performed and we compare with corresponding results in the literature. Analytic results are supplied in the supplementary materials.


title: Enhancing Energy System Models Using Better Load Forecasts

abstract: Energy system models require a large amount of technical and economic data, the quality of which significantly influences the reliability of the results. Some of the variables on the important data source ENTSO-E transparency platform, such as transmission system operators' day-ahead load forecasts, are known to be biased. These biases and high errors affect the quality of energy system models. We propose a simple time series model that does not require any input variables other than the load forecast history to significantly improve the transmission system operators' load forecast data on the ENTSO-E transparency platform in real-time, i.e., we successively improve each incoming data point. We further present an energy system model developed specifically for the short-term day-ahead market. We show that the improved load data as inputs reduce pricing errors of the model, with strong reductions particularly in times when prices are high and the market is tight.


title: Aerodynamic loads on groups of offshore wind turbine towers stored on quaysides during the pre-assembly phase

abstract: Offshore wind turbine towers are pre-assembled and temporarily held in close proximity to each other in group arrangements on port quaysides. Accurate estimates of the aerodynamic loads on the individual towers and on the overall group are essential for the safe and economic design of the quayside's supporting structures and foundations. Wind tunnel tests represent the main way to address the problem, but results may lead to overconservative designs due to inevitable mismatches in the Reynolds number. This crucial issue is dealt with here using an original engineering solution based on concentrated but discontinuous surface roughness, which allows, for the first time in the case of finite-height towers arranged in groups and subjected to an atmospheric boundary layer flow, the successful simulation of the target high Reynolds number regime. This case study assumes slender wind turbine towers with a height of 115 m, and for the sake of generality, the investigations focus principally on a cylindrical shape rather than the more complex real-world geometry. The constant diameter of the towers is determined based on the theoretical equivalence of the mean overturning moment. The rationality of this procedure is verified a posteriori using a set of measurements on the real-shape towers. The experiments show a regular behavior of the maximum mean base shear force and moment for towers arranged in double-row groups, while the results are more complicated for single-row groups, since biased flow sometimes occurs in symmetric or nearly-symmetric configurations. Dynamic loads are also inspected, and gust factors in good agreement with Eurocode 1's prescriptions are found. Several parametric studies are carried out, the most extensive of which is devoted to assessing the role of tower height. A complicated non-monotonic pattern of the load coefficients with the tower height is encountered.


title: Using Semantic Information for Defining and Detecting OOD Inputs

abstract: As machine learning models continue to achieve impressive performance across different tasks, the importance of effective anomaly detection for such models has increased as well. It is common knowledge that even well-trained models lose their ability to function effectively on out-of-distribution inputs. Thus, out-of-distribution (OOD) detection has received some attention recently. In the vast majority of cases, it uses the distribution estimated by the training dataset for OOD detection. We demonstrate that the current detectors inherit the biases in the training dataset, unfortunately. This is a serious impediment, and can potentially restrict the utility of the trained model. This can render the current OOD detectors impermeable to inputs lying outside the training distribution but with the same semantic information (e.g. training class labels). To remedy this situation, we begin by defining what should ideally be treated as an OOD, by connecting inputs with their semantic information content. We perform OOD detection on semantic information extracted from the training data of MNIST and COCO datasets and show that it not only reduces false alarms but also significantly improves the detection of OOD inputs with spurious features from the training data.


title: Comparisons of the core and mantle compositions of earth analogs from different terrestrial planet formation scenarios

abstract: The chemical compositions of Earth's core and mantle provide insight into the processes that led to their formation. N-body simulations, on the other hand, generally do not contain chemical information, and seek to only reproduce the masses and orbits of the terrestrial planets. These simulations can be grouped into four potentially viable scenarios of Solar System formation (Classical, Annulus, Grand Tack, and Early Instability) for which we compile a total of 433 N-body simulations. We relate the outputs of these simulations to the chemistry of Earth's core and mantle using a melt-scaling law combined with a multi-stage model of core formation. We find the compositions of Earth analogs to be largely governed by the fraction of equilibrating embryo cores and the initial embryo masses in N-body simulations. Simulation type may be important when considering magma ocean lifetimes, where Grand Tack simulations have the largest amounts of material accreted after the last giant impact. However, we cannot rule out any accretion scenarios or initial embryo masses due to the sensitivity of Earth's mantle composition to different parameters and the stochastic nature of N-body simulations. Comparing the last embryo impacts experienced by Earth analogs to specific Moon-forming scenarios, we find the characteristics of the Moon-forming impact are dependent on the initial conditions in N-body simulations where larger initial embryo masses promote larger and slower Moon-forming impactors. Mars-sized initial embryos are most consistent with the canonical hit-and-run scenario onto a solid mantle. Our results suggest that constraining the fraction of equilibrating impactor core and the initial embryo masses in N-body simulations could be significant for understanding both Earth's accretion history and characteristics of the Moon-forming impact.


title: MVMTnet: A Multi-variate Multi-modal Transformer for Multi-class Classification of Cardiac Irregularities Using ECG Waveforms and Clinical Notes

abstract: Deep learning provides an excellent avenue for optimizing diagnosis and patient monitoring for clinical-based applications, which can critically enhance the response time to the onset of various conditions. For cardiovascular disease, one such condition where the rising number of patients increasingly outweighs the availability of medical resources in different parts of the world, a core challenge is the automated classification of various cardiac abnormalities. Existing deep learning approaches have largely been limited to detecting the existence of an irregularity, as in binary classification, which has been achieved using networks such as CNNs and RNN/LSTMs. The next step is to accurately perform multi-class classification and determine the specific condition(s) from the inherently noisy multi-variate waveform, which is a difficult task that could benefit from (1) a more powerful sequential network, and (2) the integration of clinical notes, which provide valuable semantic and clinical context from human doctors. Recently, Transformers have emerged as the state-of-the-art architecture for forecasting and prediction using time-series data, with their multi-headed attention mechanism, and ability to process whole sequences and learn both long and short-range dependencies. The proposed novel multi-modal Transformer architecture would be able to accurately perform this task while demonstrating the cross-domain effectiveness of Transformers, establishing a method for incorporating multiple data modalities within a Transformer for classification tasks, and laying the groundwork for automating real-time patient condition monitoring in clinical and ER settings.


title: An\'alisis del argumento cosmol\'ogico kalam desde la f\'isica contempor\'anea

abstract: The Kalam cosmological argument, which is logical-philosophical based on its two premises, affirms that the universe had a cause. We carry out an analysis from the perspective of current physics to contrast or object to the premises of said argument. For this purpose, we will use the results found in contemporary physics literature related to the premises of the Kalam argument. Analyzing quantum theory and the latest experimental results that prove that quantum phenomena are indeterministic and therefore not causal. And even in deterministic theories, such as a space-time with negative curvature, causality is in doubt. On the other hand, the first moments of the universe require models that describe states before it to contrast experimental observations, such as cosmic inflation theories. Thus, we will conclude that both premises of the Kalam argument are not supported by contemporary physics.


title: Learning signatures of decision making from many individuals playing the same game

abstract: Human behavior is incredibly complex and the factors that drive decision making--from instinct, to strategy, to biases between individuals--often vary over multiple timescales. In this paper, we design a predictive framework that learns representations to encode an individual's 'behavioral style', i.e. long-term behavioral trends, while simultaneously predicting future actions and choices. The model explicitly separates representations into three latent spaces: the recent past space, the short-term space, and the long-term space where we hope to capture individual differences. To simultaneously extract both global and local variables from complex human behavior, our method combines a multi-scale temporal convolutional network with latent prediction tasks, where we encourage embeddings across the entire sequence, as well as subsets of the sequence, to be mapped to similar points in the latent space. We develop and apply our method to a large-scale behavioral dataset from 1,000 humans playing a 3-armed bandit task, and analyze what our model's resulting embeddings reveal about the human decision making process. In addition to predicting future choices, we show that our model can learn rich representations of human behavior over multiple timescales and provide signatures of differences in individuals.


title: Gradient Flows for Sampling: Mean-Field Models, Gaussian Approximations and Affine Invariance

abstract: Sampling a probability distribution with an unknown normalization constant is a fundamental problem in computational science and engineering. This task may be cast as an optimization problem over all probability measures, and an initial distribution can be evolved to the desired minimizer dynamically via gradient flows. Mean-field models, whose law is governed by the gradient flow in the space of probability measures, may also be identified; particle approximations of these mean-field models form the basis of algorithms. The gradient flow approach is also the basis of algorithms for variational inference, in which the optimization is performed over a parameterized family of probability distributions such as Gaussians, and the underlying gradient flow is restricted to the parameterized family. By choosing different energy functionals and metrics for the gradient flow, different algorithms with different convergence properties arise. In this paper, we concentrate on the Kullback-Leibler divergence after showing that, up to scaling, it has the unique property that the gradient flows resulting from this choice of energy do not depend on the normalization constant. For the metrics, we focus on variants of the Fisher-Rao, Wasserstein, and Stein metrics; we introduce the affine invariance property for gradient flows, and their corresponding mean-field models, determine whether a given metric leads to affine invariance, and modify it to make it affine invariant if it does not. We study the resulting gradient flows in both probability density space and Gaussian space. The flow in the Gaussian space may be understood as a Gaussian approximation of the flow. We demonstrate that the Gaussian approximation based on the metric and through moment closure coincide, establish connections between them, and study their long-time convergence properties showing the advantages of affine invariance.


title: Asteroseismology of $\delta$ Scuti stars: emulating model grids using a neural network

abstract: Young $\delta$ Scuti stars have proven to be valuable asteroseismic targets but obtaining robust uncertainties on their inferred properties is challenging. We aim to quantify the random uncertainties in grid-based modelling of $\delta$ Sct stars. We apply Bayesian inference using nested sampling and a neural network emulator of stellar models, testing our method on both simulated and real stars. Based on results from simulated stars we demonstrate that our method can recover plausible posterior probability density estimates while accounting for both the random uncertainty from the observations and neural network emulation. We find that the posterior distributions of the fundamental parameters can be significantly non-Gaussian, multi-modal, and have strong covariance. We conclude that our method reliably estimates the random uncertainty in the modelling of $\delta$ Sct stars and paves the way for the investigation and quantification of the systematic uncertainty.


title: Unsourced Multiple Access with Common Alarm Messages: Network Slicing for Massive and Critical IoT

abstract: We investigate the coexistence of massive and critical Internet of Things (IoT) services in the context of the unsourced multiple access (UMA) framework introduced by Polyanskiy (2017), where all users employ a common codebook and the receiver returns an unordered list of decoded codewords. This setup is suitably modified to introduce heterogeneous traffic. Specifically, to model the massive IoT service, a standard message originates independently from each IoT device as in the standard UMA setup. To model the critical IoT service, we assume the generation of alarm messages that are common for all devices. This setup requires a significant redefinition of the error events, i.e., misdetections and false positives. We further assume that the number of active users in each transmission attempt is random and unknown. We derive a random-coding achievability bound on the misdetection and false positive probabilities of both standard and alarm messages on the Gaussian multiple access channel. Using our bound, we demonstrate that orthogonal network slicing enables massive and critical IoT to coexist under the requirement of high energy efficiency. On the contrary, we show that nonorthogonal network slicing is energy inefficient due to the residual interference from the alarm signal when decoding the standard messages.


title: Analysis of Real-Time Hostile Activitiy Detection from Spatiotemporal Features Using Time Distributed Deep CNNs, RNNs and Attention-Based Mechanisms

abstract: Real-time video surveillance, through CCTV camera systems has become essential for ensuring public safety which is a priority today. Although CCTV cameras help a lot in increasing security, these systems require constant human interaction and monitoring. To eradicate this issue, intelligent surveillance systems can be built using deep learning video classification techniques that can help us automate surveillance systems to detect violence as it happens. In this research, we explore deep learning video classification techniques to detect violence as they are happening. Traditional image classification techniques fall short when it comes to classifying videos as they attempt to classify each frame separately for which the predictions start to flicker. Therefore, many researchers are coming up with video classification techniques that consider spatiotemporal features while classifying. However, deploying these deep learning models with methods such as skeleton points obtained through pose estimation and optical flow obtained through depth sensors, are not always practical in an IoT environment. Although these techniques ensure a higher accuracy score, they are computationally heavier. Keeping these constraints in mind, we experimented with various video classification and action recognition techniques such as ConvLSTM, LRCN (with both custom CNN layers and VGG-16 as feature extractor) CNNTransformer and C3D. We achieved a test accuracy of 80% on ConvLSTM, 83.33% on CNN-BiLSTM, 70% on VGG16-BiLstm ,76.76% on CNN-Transformer and 80% on C3D.


title: Kinematics signature of a giant planet in the disk of AS 209

abstract: [abridged] ALMA observations of dust in protoplanetary disks are revealing the existence of sub-structures such as rings, gaps and cavities. Such morphology are expected to be the outcome of dynamical interaction between the disk and planets. However, other mechanisms are able to produce similar dust sub-structures. A solution is to look at the perturbation induced by the planet to the gas surface density and/or to the kinematics. In the case of the disk around AS 209, a prominent gap has been reported in the surface density of CO at $r \sim 100\,$au. Recently, Bae et al. (2022) detected a localized velocity perturbation in the $^{12}$CO $J=2-1$ emission along with a clump in $^{13}$CO $J=2-1$ at nearly 200 au, interpreted as a gaseous circumplanetary disk. We report a new analysis of ALMA archival observations of $^{12}$CO and $^{13}$CO J=2-1. A clear kinematics perturbation (kink) is detected in multiple channels and over a wide azimuth range in both dataset. We compared the observed perturbation with a semi-analytic model of velocity perturbations due to planet-disk interaction. The observed kink is not consistent with a planet at 200\,au as this would require a low gas disk scale height ($< 0.05$) in contradiction with previous estimate ($h/r \sim 0.118$ at $r = 100$ au). When we fix the disk scale height to 0.118 (at $r = 100$ au) we find instead that a planet of 3-5 M$_{\rm Jup}$ at 100 au induces a kinematics perturbation similar to the observed one. Thus, we conclude that a giant protoplanet orbiting at $r \sim 100\,$au is responsible of the large scale kink as well as of the perturbed dust and gas surface density previously detected. The position angle of the planet is constrained to be between 60$^{\circ}$-100$^{\circ}$. Future observations with high contrast imaging technique in the near- and mid- infrared are needed to confirm the presence and position of such a planet.


title: Ponderomotive Recoil for Electromagnetic Waves

abstract: When waves damp or amplify on resonant particles in a plasma, the nonresonant particles experience a recoil force that conserves the total momentum between the particles and electromagnetic fields. This force is important to understand, as it can completely negate current drive and rotation drive mechanisms that are predicted on the basis of only the resonant particles. Here, the existing electrostatic theory of this recoil force is extended to electromagnetic waves. While the result bears close similarity to historical fluid theories of laser-plasma interactions, it now incorporates both resonant and nonresonant particles, allowing momentum conservation to be self-consistently proven. Furthermore, the result is shown to be generally valid for kinetic plasmas, which is verified through single-particle hot-plasma simulations. The new form of the force provides physical insight into the nature of the generalized Minkowski (plasmon) momentum of geometrical optics, which is shown to correspond to the momentum gained by the field and nonresonant particles as the wave is self-consistently ramped up from vanishing amplitude.


title: Log $p$-divisible groups and semi-stable representations

abstract: Let $\mathscr{O}_K$ be a henselian DVR with field of fractions $K$ and residue field of characteristic $p>0$. Let $S$ denote $\mathop{\mathrm{Spec}} \mathscr{O}_K$ endowed with the canonical log structure. We show that the generic fiber functor $\mathbf{BT}_{S, {\mathrm{d}}}^{\log}\to \mathbf{BT}^{\mathrm{st}}_K$ between the category of dual representable log $p$-divisible groups over $S$ and the category of $p$-divisible groups with semistable reduction over $K$ is an equivalence. If $\mathscr{O}_K$ is further complete and of mixed characteristic, we show that $\mathbf{BT}_{S, {\mathrm{d}}}^{\log}$ is also equivalent to the category of semistable Galois $\mathbb{Z}_p$-representations with Hodge-Tate weights in $\{0,1\}$. Finally, we show that the above equivalences respect monodromies.


title: Two-parabolic-generator subgroups of hyperbolic 3-manifold groups

abstract: We give a detailed account of Agol's theorem and his proof concerning two-meridional-generator subgroups of hyperbolic 2-bridge link groups, which is included in the slide of his talk at the Bolyai conference 2001. We also give a generalization of the theorem to two-parabolic-generator subgroups of hyperbolic 3-manifold groups, which gives a refinement of a result due to Boileau-Weidmann.


title: Boosting Nystr\"{o}m Method

abstract: The Nystr\"{o}m method is an effective tool to generate low-rank approximations of large matrices, and it is particularly useful for kernel-based learning. To improve the standard Nystr\"{o}m approximation, ensemble Nystr\"{o}m algorithms compute a mixture of Nystr\"{o}m approximations which are generated independently based on column resampling. We propose a new family of algorithms, boosting Nystr\"{o}m, which iteratively generate multiple ``weak'' Nystr\"{o}m approximations (each using a small number of columns) in a sequence adaptively - each approximation aims to compensate for the weaknesses of its predecessor - and then combine them to form one strong approximation. We demonstrate that our boosting Nystr\"{o}m algorithms can yield more efficient and accurate low-rank approximations to kernel matrices. Improvements over the standard and ensemble Nystr\"{o}m methods are illustrated by simulation studies and real-world data analysis.


title: MultiVehicle Simulator (MVSim): lightweight dynamics simulator for multiagents and mobile robotics research

abstract: Development of applications related to closed-loop control requires either testing on the field or on a realistic simulator, with the latter being more convenient, inexpensive, safe, and leading to shorter development cycles. To address that need, the present work introduces MVSim, a simulator for multiple vehicles or robots capable of running dozens of agents in simple scenarios, or a handful of them in complex scenarios. MVSim employs realistic physics-grounded friction models for tire-ground interaction, and aims at accurate and GPU-accelerated simulation of most common modern sensors employed in mobile robotics and autonomous vehicle research, such as depth and RGB cameras, or 2D and 3D LiDAR scanners. All depth-related sensors are able to accurately measure distances to 3D models provided by the user to define custom world elements. Efficient simulation is achieved by means of focusing on ground vehicles, which allows the use of a simplified 2D physics engine for body collisions while solving wheel-ground interaction forces separately. The core parts of the system are written in C++ for maximum efficiency, while Python, ROS 1, and ROS 2 wrappers are also offered for easy integration into user systems. A custom publish/subscribe protocol based on ZeroMQ (ZMQ) is defined to allow for multiprocess applications to access or modify a running simulation. This simulator enables and makes easier to do research and development on vehicular dynamics, autonomous navigation algorithms, and simultaneous localization and mapping (SLAM) methods.


title: Counterfeit Chip Detection using Scattering Parameter Analysis

abstract: The increase in the number of counterfeit and recycled microelectronic chips in recent years has created significant security and safety concerns in various applications. Hence, detecting such counterfeit chips in electronic systems is critical before deployment in the field. Unfortunately, the conventional verification tools using physical inspection and side-channel methods are costly, unscalable, error-prone, and often incompatible with legacy systems. This paper introduces a generic non-invasive and low-cost counterfeit chip detection based on characterizing the impedance of the system's power delivery network (PDN). Our method relies on the fact that the impedance of the counterfeit and recycled chips differs from the genuine ones. To sense such impedance variations confidently, we deploy scattering parameters, frequently used for impedance characterization of RF/microwave circuits. Our proposed approach can directly be applied to soldered chips on the system's PCB and does not require any modifications on the legacy systems. To validate our claims, we perform extensive measurements on genuine and aged samples from two families of STMicroelectronics chips to assess the effectiveness of the proposed approach.


title: Color-avoiding connected spanning subgraphs with minimum number of edges

abstract: We call a (not necessarily properly) edge-colored graph edge-color-avoiding connected if after the removal of edges of any single color, the graph remains connected. For vertex-colored graphs, similar definitions of color-avoiding connectivity can be given. In this article, we investigate the problem of determining the maximum number of edges that can be removed from a color-avoiding connected graph so that it remains color-avoiding connected. First, we prove that this problem is NP-hard, then we give a polynomial-time approximation algorithm for it. To analyze the approximation factor of this algorithm, we determine the minimum number of edges of color-avoiding connected graphs on a given number of vertices and with a given number of colors. Furthermore, we also consider a generalization of edge-color-avoiding connectivity to matroids.


title: CrowdLogo: crowd simulation in NetLogo

abstract: Planning the evacuation of people from crowded places, such as squares, stadiums, or indoor arenas during emergency scenarios is a fundamental task that authorities must deal with. This article summarizes the work of the authors to simulate an emergency scenario in a square using NetLogo, a multi-agent programmable modeling environment. The emergency scenario is based on a real event, which took place in Piazza San Carlo, Turin, on the 3rd of June 2017. The authors have developed a model and conducted various experiments, the results of which are presented, discussed and analyzed. The article concludes by offering suggestions for further research and summarizing the key takeaways.


title: Sharp estimates for imaginary powers of Bessel operators

abstract: Let $L f(x):=-\frac{d^2}{dx^2}f(x)-\frac{ r}{x}\frac{d}{dx}f(x),\quad x>0, r>0$ be the Bessel operator on $((0,\infty), |\cdot|, x^rdx)$. In this paper, we prove the sharp weak type $(1,1)$ estimate for the imaginary power $L^{i\alpha}, \alpha\in \mathbb R$, of the Bessel operator.


title: Preliminary Examination of Guardian Cap Head Impact Data Using Instrumented Mouthguards

abstract: Purpose The objective of this study is to present preliminary on-field head kinematics data for NCAA Division I American football players through closely matched pre-season workouts both with and without Guardian Caps (GCs). Methods 42 NCAA Division I American football players wore instrumented mouthguards (iMMs) for 6 closely matched workouts, 3 in traditional helmets (PRE) and 3 with GCs (POST) affixed to the exterior of their helmets. This includes 7 players who had consistent data through all workouts. Results There was no significant difference between the collapsed mean values for the entire sample between PRE and POST for peak linear acceleration (PLA) (PRE=16.3, POST=17.2Gs; p=0.20), Peak Angular Acceleration (PAA) (PRE=992.1, POST=1029.4rad/s2; p=0.51 and the total amount of impacts (PRE=9.3, POST=9.7; p=0.72). Similarly, no difference was observed between PRE and POST for PLA (PRE=16.1, POST=17.2Gs; p=0.32), PAA (PRE=951.2, POST=1038.0rad/s2; p=0.29 and total impacts (PRE=9.6, POST=9.7; p=0.32) between sessions for the 7 repeated players. Conclusion These data suggest no difference in head kinematics data (PLA, PAA and total impacts) when GCs are worn. This study suggests GCs are not effective in reducing the magnitude of head impacts experienced by NCAA Division I American football players.


title: The Lefschetz property for an algebra defined by matchings

abstract: In this article, we consider the weighted generating function of matchings in the complete graph. We define an Artinian Gorenstein algebra as the quotient ring of a polynomial ring by the annihilator of the generating function. We show the strong Lefschetz property of the algebra.


title: An Accelerated Asynchronous Distributed Method for Convex Constrained Optimization Problems

abstract: We consider a class of multi-agent cooperative consensus optimization problems with local nonlinear convex constraints where only those agents connected by an edge can directly communicate, hence, the optimal consensus decision lies in the intersection of these private sets. We develop an asynchronous distributed accelerated primal-dual algorithm to solve the considered problem. The proposed scheme is the first asynchronous method with an optimal convergence guarantee for this class of problems, to the best of our knowledge. In particular, we provide an optimal convergence rate of $\mathcal{O(1/K)}$ for suboptimality, infeasibility, and consensus violation.


title: Counterexamples in rotundity of norms in Banach spaces

abstract: We study several classical concepts in the topic of strict convexity of norms in infinite dimensional Banach spaces. Specifically, and in descending order of strength, we deal with Uniform Rotundity (UR), Weak Uniform Rotundity (WUR) and Uniform Rotundity in Every Direction (URED). Our first three results show that we may distinguish between all of these three properties in every Banach space where such renormings are possible. Specifically, we show that in every infinite dimensional Banach space which admits a WUR (resp. URED) renorming, we can find a norm with the same condition and which moreover fails to be UR (resp. WUR). We prove that these norms can be constructed to be Locally Uniformly Rotund (LUR) in Banach spaces admitting such renormings. Additionally, we obtain that in every Banach space with a LUR norm we can find a LUR renorming which is not URED. These results solve three open problems posed by A.J. Guirao, V. Montesinos and V. Zizler. The norms we construct in this first part are dense. In the last part of this note, we solve a fourth question posed by the same three authors by constructing a $C^\infty$-smooth norm in $c_0$ whose dual norm is not strictly convex.


title: In-context Example Selection with Influences

abstract: In-context learning (ICL) is a powerful paradigm emerged from large language models (LLMs). Despite its promises, ICL performance is known to be highly sensitive to input examples. In this work, we use $\textit{in-context influences}$ to analyze few-shot ICL performance directly from the in-context examples. Our proposed influence-based example selection method can identify both positive and negative examples, outperforming several baselines when evaluated on 9 SuperGLUE tasks. Our analysis uncovers up to a $16.3\%$ performance gap between using the most negative in-context examples compared to the most positive. In a case study, we apply our influence-based framework to quantify the phenomena of recency bias in example ordering for few-shot ICL.


title: Constructing Deterministic Parity Automata from Positive and Negative Examples

abstract: We present a polynomial time algorithm that constructs a deterministic parity automaton (DPA) from a given set of positive and negative ultimately periodic example words. We show that this algorithm is complete for the class of $\omega$-regular languages, that is, it can learn a DPA for each regular $\omega$-language. For use in the algorithm, we give a definition of a DPA, that we call the precise DPA of a language, and show that it can be constructed from the syntactic family of right congruences for that language (introduced by Maler and Staiger in 1997). Depending on the structure of the language, the precise DPA can be of exponential size compared to a minimal DPA, but it can also be a minimal DPA. The upper bound that we obtain on the number of examples required for our algorithm to find a DPA for $L$ is therefore exponential in the size of a minimal DPA, in general. However we identify two parameters of regular $\omega$-languages such that fixing these parameters makes the bound polynomial.


title: The Target-Charging Technique for Privacy Accounting across Interactive Computations

abstract: We propose the \emph{Target Charging Technique} (TCT), a unified privacy analysis framework for interactive settings where a sensitive dataset is accessed multiple times using differentially private algorithms. Unlike traditional composition, where privacy guarantees deteriorate quickly with the number of accesses, TCT allows computations that don't hit a specified \emph{target}, often the vast majority, to be essentially free (while incurring instead a small overhead on those that do hit their targets). TCT generalizes tools such as the sparse vector technique and top-$k$ selection from private candidates and extends their remarkable privacy enhancement benefits from noisy Lipschitz functions to general private algorithms.


title: A Montgomery-Hooley theorem for the k-fold divisor function

abstract: Let $d_k(n)$ denote the $k$-fold divisor function. For a wide range of large $q$ the expected bound $$\sum_{n\leq x\atop {n\equiv a(q)}}d_k(n)-\text { main term }\approx \sqrt {x/q}$$ is shown to be true in an average sense -- for all $k$. This generalises the work of Pongsriiam and Vaughan [15] who studied $k=2$, and answers the work of Rodgers and Soundararajan [17], who used the asymptotic large sieve to study a smoothed version of the problem. We use a circle method approach as developed by Goldston and Vaughan [7] to study the unsmoothed problem.


title: Teachable Reality: Prototyping Tangible Augmented Reality with Everyday Objects by Leveraging Interactive Machine Teaching

abstract: This paper introduces Teachable Reality, an augmented reality (AR) prototyping tool for creating interactive tangible AR applications with arbitrary everyday objects. Teachable Reality leverages vision-based interactive machine teaching (e.g., Teachable Machine), which captures real-world interactions for AR prototyping. It identifies the user-defined tangible and gestural interactions using an on-demand computer vision model. Based on this, the user can easily create functional AR prototypes without programming, enabled by a trigger-action authoring interface. Therefore, our approach allows the flexibility, customizability, and generalizability of tangible AR applications that can address the limitation of current marker-based approaches. We explore the design space and demonstrate various AR prototypes, which include tangible and deformable interfaces, context-aware assistants, and body-driven AR applications. The results of our user study and expert interviews confirm that our approach can lower the barrier to creating functional AR prototypes while also allowing flexible and general-purpose prototyping experiences.


title: Eight-node solid brick element high-order stiffness matrix template

abstract: In this paper, the template will be developed from an assumed Stress Method, which its formulation is based on the Hellinger-Reissner principle developed according to Kang's study in 1986. The element stiffness is decomposed into a basic part that takes care of consistency and mix-ability, and a HO element stiffness part that takes care of stability (also known as rank sufficient) and accuracy. In the FE method, the HO stiffness is based on a displacement formulation, whereas the basis stiffness is method independent. To start, one should be familiar with the definition of a solid brick element. Solid brick element is three-dimensional finite elements that can model solid bodies and structures without any a priori geometric simplification.


abstract: We propose a novel model-based offline Reinforcement Learning (RL) framework, called Adversarial Model for Offline Reinforcement Learning (ARMOR), which can robustly learn policies to improve upon an arbitrary reference policy regardless of data coverage. ARMOR is designed to optimize policies for the worst-case performance relative to the reference policy through adversarially training a Markov decision process model. In theory, we prove that ARMOR, with a well-tuned hyperparameter, can compete with the best policy within data coverage when the reference policy is supported by the data. At the same time, ARMOR is robust to hyperparameter choices: the policy learned by ARMOR, with "any" admissible hyperparameter, would never degrade the performance of the reference policy, even when the reference policy is not covered by the dataset. To validate these properties in practice, we design a scalable implementation of ARMOR, which by adversarial training, can optimize policies without using model ensembles in contrast to typical model-based methods. We show that ARMOR achieves competent performance with both state-of-the-art offline model-free and model-based RL algorithms and can robustly improve the reference policy over various hyperparameter choices.


title: Framework for Certification of AI-Based Systems

abstract: The current certification process for aerospace software is not adapted to "AI-based" algorithms such as deep neural networks. Unlike traditional aerospace software, the precise parameters optimized during neural network training are as important as (or more than) the code processing the network and they are not directly mathematically understandable. Despite their lack of explainability such algorithms are appealing because for some applications they can exhibit high performance unattainable with any traditional explicit line-by-line software methods. This paper proposes a framework and principles that could be used to establish certification methods for neural network models for which the current certification processes such as DO-178 cannot be applied. While it is not a magic recipe, it is a set of common sense steps that will allow the applicant and the regulator increase their confidence in the developed software, by demonstrating the capabilities to bring together, trace, and track the requirements, data, software, training process, and test results.


title: Edgeformers: Graph-Empowered Transformers for Representation Learning on Textual-Edge Networks

abstract: Edges in many real-world social/information networks are associated with rich text information (e.g., user-user communications or user-product reviews). However, mainstream network representation learning models focus on propagating and aggregating node attributes, lacking specific designs to utilize text semantics on edges. While there exist edge-aware graph neural networks, they directly initialize edge attributes as a feature vector, which cannot fully capture the contextualized text semantics of edges. In this paper, we propose Edgeformers, a framework built upon graph-enhanced Transformers, to perform edge and node representation learning by modeling texts on edges in a contextualized way. Specifically, in edge representation learning, we inject network information into each Transformer layer when encoding edge texts; in node representation learning, we aggregate edge representations through an attention mechanism within each node's ego-graph. On five public datasets from three different domains, Edgeformers consistently outperform state-of-the-art baselines in edge classification and link prediction, demonstrating the efficacy in learning edge and node representations, respectively.


title: Fusion of Global and Local Knowledge for Personalized Federated Learning

abstract: Personalized federated learning, as a variant of federated learning, trains customized models for clients using their heterogeneously distributed data. However, it is still inconclusive about how to design personalized models with better representation of shared global knowledge and personalized pattern. To bridge the gap, we in this paper explore personalized models with low-rank and sparse decomposition. Specifically, we employ proper regularization to extract a low-rank global knowledge representation (GKR), so as to distill global knowledge into a compact representation. Subsequently, we employ a sparse component over the obtained GKR to fuse the personalized pattern into the global knowledge. As a solution, we propose a two-stage proximal-based algorithm named \textbf{Fed}erated learning with mixed \textbf{S}parse and \textbf{L}ow-\textbf{R}ank representation (FedSLR) to efficiently search for the mixed models. Theoretically, under proper assumptions, we show that the GKR trained by FedSLR can at least sub-linearly converge to a stationary point of the regularized problem, and that the sparse component being fused can converge to its stationary point under proper settings. Extensive experiments also demonstrate the superior empirical performance of FedSLR. Moreover, FedSLR reduces the number of parameters, and lowers the down-link communication complexity, which are all desirable for federated learning algorithms. Source code is available in \url{https://github.com/huangtiansheng/fedslr}.


title: Que2Engage: Embedding-based Retrieval for Relevant and Engaging Products at Facebook Marketplace

abstract: Embedding-based Retrieval (EBR) in e-commerce search is a powerful search retrieval technique to address semantic matches between search queries and products. However, commercial search engines like Facebook Marketplace Search are complex multi-stage systems optimized for multiple business objectives. At Facebook Marketplace, search retrieval focuses on matching search queries with relevant products, while search ranking puts more emphasis on contextual signals to up-rank the more engaging products. As a result, the end-to-end searcher experience is a function of both relevance and engagement, and the interaction between different stages of the system. This presents challenges to EBR systems in order to optimize for better searcher experiences. In this paper we presents Que2Engage, a search EBR system built towards bridging the gap between retrieval and ranking for end-to-end optimizations. Que2Engage takes a multimodal & multitask approach to infuse contextual information into the retrieval stage and to balance different business objectives. We show the effectiveness of our approach via a multitask evaluation framework and thorough baseline comparisons and ablation studies. Que2Engage is deployed on Facebook Marketplace Search and shows significant improvements in searcher engagement in two weeks of A/B testing.


title: ChameleonControl: Teleoperating Real Human Surrogates through Mixed Reality Gestural Guidance for Remote Hands-on Classrooms

abstract: We present ChameleonControl, a real-human teleoperation system for scalable remote instruction in hands-on classrooms. In contrast to existing video or AR/VR-based remote hands-on education, ChameleonControl uses a real human as a surrogate of a remote instructor. Building on existing human-based telepresence approaches, we contribute a novel method to teleoperate a human surrogate through synchronized mixed reality hand gestural navigation and verbal communication. By overlaying the remote instructor's virtual hands in the local user's MR view, the remote instructor can guide and control the local user as if they were physically present. This allows the local user/surrogate to synchronize their hand movements and gestures with the remote instructor, effectively teleoperating a real human. We deploy and evaluate our system in classrooms of physiotherapy training, as well as other application domains such as mechanical assembly, sign language and cooking lessons. The study results confirm that our approach can increase engagement and the sense of co-presence, showing potential for the future of remote hands-on classrooms.


title: Conversational Text-to-SQL: An Odyssey into State-of-the-Art and Challenges Ahead

abstract: Conversational, multi-turn, text-to-SQL (CoSQL) tasks map natural language utterances in a dialogue to SQL queries. State-of-the-art (SOTA) systems use large, pre-trained and finetuned language models, such as the T5-family, in conjunction with constrained decoding. With multi-tasking (MT) over coherent tasks with discrete prompts during training, we improve over specialized text-to-SQL T5-family models. Based on Oracle analyses over n-best hypotheses, we apply a query plan model and a schema linking algorithm as rerankers. Combining MT and reranking, our results using T5-3B show absolute accuracy improvements of 1.0% in exact match and 3.4% in execution match over a SOTA baseline on CoSQL. While these gains consistently manifest at turn level, context dependent turns are considerably harder. We conduct studies to tease apart errors attributable to domain and compositional generalization, with the latter remaining a challenge for multi-turn conversations, especially in generating SQL with unseen parse trees.


title: SGD learning on neural networks: leap complexity and saddle-to-saddle dynamics

abstract: We investigate the time complexity of SGD learning on fully-connected neural networks with isotropic data. We put forward a complexity measure -- the leap -- which measures how "hierarchical" target functions are. For $d$-dimensional uniform Boolean or isotropic Gaussian data, our main conjecture states that the time complexity to learn a function $f$ with low-dimensional support is $\tilde\Theta (d^{\max(\mathrm{Leap}(f),2)})$. We prove a version of this conjecture for a class of functions on Gaussian isotropic data and 2-layer neural networks, under additional technical assumptions on how SGD is run. We show that the training sequentially learns the function support with a saddle-to-saddle dynamic. Our result departs from [Abbe et al. 2022] by going beyond leap 1 (merged-staircase functions), and by going beyond the mean-field and gradient flow approximations that prohibit the full complexity control obtained here. Finally, we note that this gives an SGD complexity for the full training trajectory that matches that of Correlational Statistical Query (CSQ) lower-bounds.


title: Exchange-Driven Intermixing of Bulk and Topological Surface State by Chiral Excitons in Bi2Se3

abstract: Topological surface states (TSS) in the prototypical topological insulator (TI) Bi2Se3 are frequently characterized using optical probes, but electron-hole interactions and their effect on surface localization and optical response of the TSS remain unexplored. Here, we use ab initio calculations to understand excitonic effects in the bulk and surface of Bi2Se3. We identify multiple series of chiral excitons that exhibit both bulk and TSS character, due to exchange-driven mixing. Our results address fundamental questions about the degree to which electron-hole interactions can relax the topological protection of surface states and dipole selection rules for circularly polarized light in TIs by elucidating the complex intermixture of bulk and surface states excited in optical measurements and their coupling to light.


title: Indium-Bond-And-Stop-Etch (IBASE) Technique for Dual-side Processing of Thin High-mobility GaAs/AlGaAs Epitaxial Layers

abstract: We present a reliable flip-chip technique for dual-side processing of thin (<1 micron) high-mobility GaAs/AlGaAs epitaxial layers. The technique allows the fabrication of small (micron-scale with standard UV photolithography) patterned back gates and dual-gate structures on the thin GaAs/AlGaAs films with good alignment accuracy using only frontside alignment. The technique preserves the high-mobility (>10^6 cm^2 /V-s at 2 K) and most (>95%) of the charge density of the 2-dimensional electron gas (2DEG) systems, and allows linear control of the charge density with small (< 1 V) electrostatic gate bias. Our technique is motivated by a novel THz quantum-well detector based on intersubband transitions in a single, wide GaAs/AlGaAs quantum well, in which a symmetric, well-aligned dual-gate structure (with a typical gate dimension of ~5 micron by 5 micron) is required for accurate and precise tuning of the THz detection frequency. Using our Indium-Bond-And-Stop-Etch (IBASE) technique, we realize such dual-gate structure on 660-nm thick GaAs/AlGaAs epitaxial layers that contain a modulation-doped, 40-nm wide, single square quantum well. By independently controlling the charge density and the DC electric field set between the gates, we demonstrate robust tuning of the intersubband absorption behavior of the 40-nm quantum well near 3.44 THz at 30 K.


title: Bayesian Age Category Reconciliation for Age- and Cause-specific Under-five Mortality Estimates

abstract: Age-disaggregated health data is crucial for effective public health planning and monitoring. Monitoring under-five mortality, for example, requires highly detailed age data since the distribution of potential causes of death varies substantially within the first few years of life. Comparative researchers often have to rely on multiple data sources yet, these sources often have ages aggregated at different levels, making it difficult to combine the data into a single, coherent picture. To address this challenge in the context of under-five cause-specific mortality, we propose a Bayesian approach, that calibrates data with different age structures to produce unified and accurate estimates of the standardized age group distributions. We consider age-disaggregated death counts as fully-classified multinomial data and show that by incorporating partially-classified aggregated data, we can construct an improved Bayes estimator of the multinomial parameters under the Kullback-Leibler (KL) loss. We illustrate the method using both synthetic and real data, demonstrating that the proposed method achieves adequate performance in imputing incomplete classification. Finally, we present the results of numerical studies examining the conditions necessary for obtaining improved estimators. These studies provide insights and interpretations that can be used to aid future research and inform guidance for practitioners on appropriate levels of age disaggregation, with the aim of improving the accuracy and reliability of under-five cause-specific mortality estimates.


title: Loop corrections in a solvable UV-finite model and its effective field theory

abstract: We examine some features of the non-renormalizability induced through the use of low-energy effective Lagrangians in loop diagrams, in the context of a scalar model which is ultraviolet finite and partially soluble. In this framework, one can directly demonstrate the mechanism leading to the non-renormalizability of the effective theory. This behavior is generated by approximations that are applicable at low energies but are generally inappropriate for evaluating loop diagrams that contain virtual high-energy particles. However, it is explicitly shown that one can match the results obtained in the renormalized effective theory with those found in the full theory at low energy. We argue that the infrared sectors of these theories are inherently similar, independently of the matching procedure. A closed-form expression is obtained, to leading order in the energy expansion, for the complete effective Lagrangian at all orders in perturbation theory. The model may be useful to clarify certain aspects of realistic, but more complex, effective field theories.


title: Defect-induced states, defect-induced phase transition and excitonic states in bent transition metal dichalcogenide (TMD) nanoribbons: density functional vs. many body theory

abstract: Two-dimensional (2D) transition metal dichalcogenide (TMD) materials have versatile electronic and optical properties. TMD nanoribbons show interesting properties due to reduced dimensionality, quantum confinement, and edge states. Tang et al. showed that the edge bands evolved with bending can tune the optical properties for various widths of TMD nanoribbons. Defects are commonly present in 2D TMD materials, and can dramatically change the material properties. In this following work, we investigate the interaction between the edge and the defect states in WS2 nanoribbons with line defects under different bending conditions, using density functional theory (DFT). We reveal interesting semiconducting-to-metallic phase transitions, suggesting potential applications in nano-electronics or molecular electronics. We also calculate the optical absorption of the nanoribbons with different defect positions with the many-body GW-BSE (Bethe-Salpeter equation) approach, revealing a tunable optical spectrum and diverse exciton states in the defected TMD nanoribbons.


title: Physics-informed Spectral Learning: the Discrete Helmholtz--Hodge Decomposition

abstract: In this work, we further develop the Physics-informed Spectral Learning (PiSL) by Espath et al. \cite{Esp21} based on a discrete $L^2$ projection to solve the discrete Hodge--Helmholtz decomposition from sparse data. Within this physics-informed statistical learning framework, we adaptively build a sparse set of Fourier basis functions with corresponding coefficients by solving a sequence of minimization problems where the set of basis functions is augmented greedily at each optimization problem. Moreover, our PiSL computational framework enjoys spectral (exponential) convergence. We regularize the minimization problems with the seminorm of the fractional Sobolev space in a Tikhonov fashion. In the Fourier setting, the divergence- and curl-free constraints become a finite set of linear algebraic equations. The proposed computational framework combines supervised and unsupervised learning techniques in that we use data concomitantly with the projection onto divergence- and curl-free spaces. We assess the capabilities of our method in various numerical examples including the `Storm of the Century' with satellite data from 1993.


title: A Log-linear Gradient Descent Algorithm for Unbalanced Binary Classification using the All Pairs Squared Hinge Loss

abstract: Receiver Operating Characteristic (ROC) curves are plots of true positive rate versus false positive rate which are used to evaluate binary classification algorithms. Because the Area Under the Curve (AUC) is a constant function of the predicted values, learning algorithms instead optimize convex relaxations which involve a sum over all pairs of labeled positive and negative examples. Naive learning algorithms compute the gradient in quadratic time, which is too slow for learning using large batch sizes. We propose a new functional representation of the square loss and squared hinge loss, which results in algorithms that compute the gradient in either linear or log-linear time, and makes it possible to use gradient descent learning with large batch sizes. In our empirical study of supervised binary classification problems, we show that our new algorithm can achieve higher test AUC values on imbalanced data sets than previous algorithms, and make use of larger batch sizes than were previously feasible.


title: About the quantum Talbot effect on the sphere

abstract: The Schr\"odinger equation on a circle with an initially localized profile of the wave function is known to give rise to revivals or replications, where the probability density of the particle is partially reproduced at rational times. As a consequence of the convolutional form of the general solution it is deduced that a piecewise constant initial wave function remains piecewise constant at rational times as well. For a sphere instead, it is known that this piecewise revival does not necessarily occur, indeed the wave function becomes singular at some specific locations at rational times. It may be desirable to study the same problem, but with an initial condition being a localized Dirac delta instead of a piecewise constant function, and this is the purpose of the present work. By use of certain summation formulas for the Legendre polynomials together with properties of Gaussian sums, it is found that revivals on the sphere occur at rational times for some specific locations, and the structure of singularities of the resulting wave function is characterized in detail. In addition, a partial study of the regions where the density vanishes, named before valley of shadows in the context of the circle, is initiated here. It is suggested that, differently from the circle case, these regions are not lines but instead some specific set of points along the sphere. A conjecture about the precise form of this set is stated and the intuition behind it is clarified.


title: Task-Oriented Prediction and Communication Co-Design for Haptic Communications

abstract: Prediction has recently been considered as a promising approach to meet low-latency and high-reliability requirements in long-distance haptic communications. However, most of the existing methods did not take features of tasks and the relationship between prediction and communication into account. In this paper, we propose a task-oriented prediction and communication co-design framework, where the reliability of the system depends on prediction errors and packet losses in communications. The goal is to minimize the required radio resources subject to the low-latency and high-reliability requirements of various tasks. Specifically, we consider the just noticeable difference (JND) as a performance metric for the haptic communication system. We collect experiment data from a real-world teleoperation testbed and use time-series generative adversarial networks (TimeGAN) to generate a large amount of synthetic data. This allows us to obtain the relationship between the JND threshold, prediction horizon, and the overall reliability including communication reliability and prediction reliability. We take 5G New Radio as an example to demonstrate the proposed framework and optimize bandwidth allocation and data rates of devices. Our numerical and experimental results show that the proposed framework can reduce wireless resource consumption up to 77.80% compared with a task-agnostic benchmark.


title: Probing the incompressibility of dense hadronic matter near QCD phase transition in relativistic heavy-ion collisions

abstract: Based on the extended hadronic transport model of relativistic heavy-ion collisions, the incompressibility of dense hadronic matter created in relativistic Au+Au heavy-ion collisions at $\sqrt{s_{NN}} = 3$ GeV is studied. By comparing experimental proton directed flow, productions of strange hadrons $\phi$, $K^{-}$ as well as their ratio $\phi/K^{-}$, proton high-order cumulants to the model calculations, a large incompressibility of dense hadronic matter is obtained from nucleon observabels while a rather small incompressibility is needed to fit the data of strange hadrons. This may indicate hadronic matter possesses different incompressibilities in different density regions, i.e., the incompressibility may become stiffer from saturation density to a certain baryon density and then turn to soft before reaching hadron-quark phase transition. The study also shows that the incompressibility significantly affects the critical baryon density of hadron-quark phase transition.


title: Reinforcement Learning for Block Decomposition of CAD Models

abstract: We present a novel AI-assisted method for decomposing (segmenting) planar CAD (computer-aided design) models into well shaped rectangular blocks as a proof-of-principle of a general decomposition method applicable to complex 2D and 3D CAD models. The decomposed blocks are required for generating good quality meshes (tilings of quadrilaterals or hexahedra) suitable for numerical simulations of physical systems governed by conservation laws. The problem of hexahedral mesh generation of general CAD models has vexed researchers for over 3 decades and analysts often spend more than 50% of the design-analysis cycle time decomposing complex models into simpler parts meshable by existing techniques. Our method uses reinforcement learning to train an agent to perform a series of optimal cuts on the CAD model that result in a good quality block decomposition. We show that the agent quickly learns an effective strategy for picking the location and direction of the cuts and maximizing its rewards as opposed to making random cuts. This paper is the first successful demonstration of an agent autonomously learning how to perform this block decomposition task effectively thereby holding the promise of a viable method to automate this challenging process.


title: Ending States of a Special Variant of the Chip-Firing Algorithm

abstract: We investigate a special variant of chip-firing, in which we consider an infinite set of rooms on a number line, some of which are occupied by violinists. In a move, we take two violinists in adjacent rooms, and send one of them to the closest unoccupied room to the left and the other to the closest unoccupied room to the right. We classify the different possible final states from repeatedly performing this operation. We introduce numbers $R(N,\ell,x)$ that count labeled recursive rooted trees with $N$ vertices, $\ell$ leaves, and the smallest rooted path ending in $x$. We describe the properties of these numbers and connect them to permutations. We conjecture that these numbers describe the probabilities ending with different final states when the moves are chosen uniformly.