title: Evaluating the feasibility of interpretable machine learning for globular cluster detection

abstract: Extragalactic globular clusters (GCs) are important tracers of galaxy formation and evolution. Obtaining GC catalogues from photometric data involves several steps which will likely become too time-consuming to perform on the large data volumes that are expected from upcoming wide-field imaging projects such as Euclid. In this work, we explore the feasibility of various machine learning (ML) methods to aid the search for GCs. We use archival Hubble Space Telescope data in the F475W and F850LP bands of 141 early-type galaxies in the Fornax and Virgo galaxy clusters. Using existing GC catalogues to label the data, we obtain an extensive data set of 84929 sources containing 18556 GCs and we train several ML methods both on image and tabular data containing physically relevant features extracted from the images. We find that our evaluated ML models are capable of producing catalogues of similar quality as the existing ones. The best performing methods, ensemble-based models like random forests and convolutional neural networks, recover ~ 90-94 % of GCs while producing an acceptable amount of false detections (~ 6-8 %) - with some falsely detected sources being identifiable as GCs that have not been labelled as such in the used catalogues. In the magnitude range 22 < m4_g < 24.5 mag, 98 - 99 % of GCs are recovered. We even find such high performance levels when training on Virgo and evaluating on Fornax data (and vice versa), illustrating that the models are transferable to environments with different conditions such as different distances than in the used training data. Additionally, we demonstrate how interpretable methods can be used to better understand model predictions, recovering that magnitudes, colours, and sizes are important for identifying GCs. These are encouraging results, indicating that similar methods can be applied for creating GC catalogues for a large number of galaxies.


title: Chaos enhancement in large-spin chains

abstract: We study the chaotic properties of a large-spin XXZ chain with onsite disorder and a small number of excitations above the fully polarized state. We show that while the classical limit, which is reached for large spins, is chaotic, enlarging the spin suppresses quantum chaos features. We attribute this suppression to the occurrence of large and slow clusters of onsite excitations, and propose a way to facilitate their fragmentation by introducing additional decay channels. We numerically verify that the introduction of such relaxation channels restores chaoticity for large spins, so that only three excitations are required to achieve strong level repulsion and ergodic eigenstates.


title: Kepler-1656b's Extreme Eccentricity: Signature of a Gentle Giant

abstract: Highly eccentric orbits are one of the major surprises of exoplanets relative to the Solar System and indicate rich and tumultuous dynamical histories. One system of particular interest is Kepler-1656, which hosts a sub-Jovian planet with an eccentricity of 0.8. Sufficiently eccentric orbits will shrink in semi-major axis due to tidal dissipation of orbital energy during periastron passage. Here our goal was to assess whether Kepler-1656b is currently undergoing such high-eccentricity migration, and to further understand the system's origins and architecture. We confirm a second planet in the system with $M_{\rm c}= 0.40 \pm 0.09M_{\rm jup}$ and P$_{\rm c}= 1919\pm 27\,$days. We simulated the dynamical evolution of planet b in the presence of planet c and find a variety of possible outcomes for the system, such as tidal migration and engulfment. The system is consistent with an in situ dynamical origin of planet b followed by subsequent Eccentric Kozai Lidov (EKL) perturbations that excite Kepler-1656b's eccentricity gently, i.e. without initiating tidal migration. Thus, despite its high eccentricity, we find no evidence that planet b is or has migrated through the high-eccentricity channel. Finally, we predict the outer orbit to be mutually inclined in a nearly perpendicular configuration with respect to the inner planet orbit based on the outcomes of our simulations, and make observable predictions for the inner planet's spin-orbit angle. Our methodology can be applied to other eccentric or tidally locked planets to constrain their origins, orbital configurations and properties of a potential companion.


title: Spectacular 240 kpc double-sided relativistic jets in a spiral-hosted narrow-line Seyfert 1 galaxy

abstract: Narrow-line Seyfert 1 (NLS1) galaxies are a peculiar subclass of active galactic nuclei (AGN). They have demonstrated that the presence of relativistic jets in an AGN is not strictly related to its radio-loudness, the black hole mass, or the host galaxy type. Here we present a remarkable example of a radio-quiet NLS1, 6dFGS gJ035432.8-134008 (J0354-1340). In our Karl G. Jansky Very Large Array observations at 5.5 GHz the source shows a bright core with a flat spectral index, and extended emission corresponding to very elongated jets. These are the largest double-sided radio jets found to date in an NLS1, with a deprojected linear size of almost 250 kpc. We also analysed near-infrared and optical images obtained by the Magellan Baade and the European Southern Observatory New Technology Telescope. By means of photometric decomposition and colour maps, we determined that J0354-1340 is hosted by a spiral/disk galaxy. Fully evolved relativistic jets have traditionally been associated with high-mass elliptical galaxies hosting the most massive black holes, but our results confirm that also less massive black holes in spiral galaxies can launch and sustain powerful jets, implying that the launching of the jets is governed by factors other than previously believed.


title: Looking for structure in the cobordism conjecture

abstract: The cobordism conjecture of the swampland program states that the bordism group of quantum gravity must be trivial. We investigate this statement in several directions, on both the mathematical and physical side. We consider the Whitehead tower construction as a possible organising principle for the topological structures entering the formulation of the conjecture. We discuss why and how to include geometric structures in bordism groups, such as higher U(1)-bundles with connection. The inclusion of magnetic defects is also addressed in some detail. We further elaborate on how the conjecture could predict Kaluza--Klein monopoles, and we study the gravity decoupling limit in the cobordism conjecture, with a few observations on NSNS string backgrounds. We end with comments in relation to T-duality, as well as the finiteness conjecture.


title: Hubble constant and nuclear equation of state from kilonova spectro-photometric light curves

abstract: The merger of two compact objects of which at least one is a neutron star is signalled by transient electromagnetic emission in a kilonova (KN). This event is accompanied by gravitational waves and possibly other radiation messengers such as neutrinos or cosmic rays. The electromagnetic emission arises from the radioactive decay of heavy $r-$process elements synthesized in the material ejected during and after the merger. In this paper we show that the analysis of KNe light curves can provide cosmological distance measurements and constrain the properties of the ejecta. In this respect, MAAT, the new Integral Field Unit in the OSIRIS spectrograph on the $10.4$ m Gran Telescopio CANARIAS (GTC), is well suited for the study of KNe by performing absolute spectro-photometry over the entire 3600-10000 Angstron spectral range. Here, we study the most representative cases regarding the scientific interest of KNe from binary neutron stars, and we evaluate the observational prospects and performance of MAAT on the GTC to do the following: a) study the impact of the equation of state on the KN light curve, and determine to what extent bounds on neutron star (NS) radii or compactness deriving from KN peak magnitudes can be identified and b) measure the Hubble constant, $H_0$, with precision improved by up to 40$\%$, when both gravitational wave data and photometric-light curves are used. In this context we discuss how the equation of state, the viewing angle, and the distance affect the precision and estimated value of $H_0$.

categories: astro-ph.CO astro-ph.HE astro-ph.IM gr-qc


title: Field-Based Physical Inference From Peculiar Velocity Tracers

abstract: We present a Bayesian hierarchical modelling approach to reconstruct the initial cosmic matter density field constrained by peculiar velocity observations. As our approach features a model for the gravitational evolution of dark matter to connect the initial conditions to late-time observations, it reconstructs the final density and velocity fields as natural byproducts. We implement this field-based physical inference approach by adapting the Bayesian Origin Reconstruction from Galaxies (BORG) algorithm, which explores the high-dimensional posterior through the use of Hamiltonian Monte Carlo sampling. We test the self-consistency of the method using random sets of mock tracers, and assess its accuracy in a more complex scenario where peculiar velocity tracers are non-linearly evolved mock haloes. We find that our framework self-consistently infers the initial conditions, density and velocity fields, and shows some robustness to model mis-specification. As compared to the state-of-the-art approach of constrained Gaussian random fields/Wiener filtering, our method produces more accurate final density and velocity field reconstructions. It also allows us to constrain the initial conditions by peculiar velocity observations, complementing in this aspect previous field-based approaches based on other cosmological observables.


title: Darker and brighter branes: Suppression and enhancement of photon production in a strongly coupled magnetized plasma

abstract: We extend our holographic analysis of the emission of photons by a strongly coupled plasma subject to a very intense external magnetic field. We previously showed that in a particular model, any photon produced by the plasma had to be in its only polarization state parallel to the reaction plane. In this paper we consider a construction that relaxes a formerly imposed constraint, permitting the emission of photons with either out-plane or in-plane polarization. This constitutes a completion of our former study because the fully back-reacted equations decouple for these two polarization states in such a manner that those involving the in-plane are identical to the ones we explored previously. In view of the above, part of the details concerning the calculations and of the numerical results for the differential rate of emitted photons that we present here correspond exactly to those omitted in our preceding letter. One of our main results is that the production of photons is increased by the introduction of a non-vanishing magnetic field with an intensity up to a value $B_{\vartheta}$, above which the effect is reversed and said production becomes lower than the $B=0$ case. The characteristic intensity $B_{\vartheta}$ depends on the propagation direction and tends to zero as the photon momentum becomes aligned with the magnetic field. Additionally, we also show that the magnetic field has the effect of increasing the value of the elliptic flow, providing a possible explanation for the excess measured in collision experiments. The holographic model is constructed using an effective five-dimensional action that includes a scalar field in addition to the constant magnetic field.


title: Stripped-envelope stars in different metallicity environments. II. Type I supernovae and compact remnants

abstract: Stripped-envelope stars can be observed as Wolf-Rayet (WR) stars, or as less luminous hydrogen-poor stars with low mass loss rates and transparent winds. Both types are potential progenitors of Type I core-collapse supernovae (SNe). We use grids of core-collapse models obtained from helium stars at different metallicities to study the effects of metallicity on the transients and remnants these stars produce. We characterise the surface and core properties of our core collapse models, and investigate their explodability employing three criteria. In cases where explosions are predicted, we estimate the ejecta mass, explosion energy, nickel mass and neutron star (NS) mass. Otherwise, we predict the mass of the resulting black hole (BH). We construct a simplified population model, and find that the properties SNe and compact objects depend strongly on metallicity. Ejecta masses and explosion energies for Type Ic SNe are best reproduced by models with Z=0.04 which exhibit strong winds during core helium burning. This implies that either their mass loss rates are underestimated, or that Type Ic SN progenitors experience mass loss through other mechanisms before exploding. The distributions of ejecta masses, explosion energies and nickel mass for Type Ib SNe are not well reproduced by progenitor models with WR mass loss, but are better reproduced if we assume no mass loss in progenitors with luminosities below the minimum WR star luminosity. We find that Type Ic SNe become more common as metallicity increases, and that the vast majority of progenitors of Type Ib SNe must be transparent-wind stripped-envelope stars. We find several models with pre-collapse CO-masses of up to $\sim 30 M_{\odot}$ may form $\sim 3 M_{\odot}$ BHs in fallback SNe. This may carry important consequences for our understanding of SNe, binary BH and NS systems, X-ray binary systems and gravitational wave transients.


title: Constraining black-hole binary spin precession and nutation with sequential prior conditioning

abstract: We investigate the detectability of sub-dominant spin effects in merging black-hole binaries using current gravitational-wave data. Using a phenomenological model that separates the spin dynamics into precession (azimuthal motion) and nutation (polar motion), we present constraints on the resulting amplitudes and frequencies. We also explore current constraints on the spin morphologies, indicating if binaries are trapped near spin-orbit resonances. We dissect such weak effects from the signals using a sequential prior conditioning approach, where parameters are progressively re-sampled from their posterior distribution. This allows us to investigate whether the data contain additional information beyond what is already provided by quantities that are better measured, namely the masses and the effective spin. For the current catalog of events, we find no significant measurements of weak spin effects such as nutation and spin-orbit locking. We synthesize a source with a high nutational amplitude and show that near-future detections will allow us to place powerful constraints, hinting that we may be at the cusp of detecting spin nutations in gravitational-wave data.


title: Soft gamma-ray polarimetry with COSI using maximum likelihood analysis

abstract: Measurements of the linear polarization of high-energy emission from pulsars, accreting black holes, and gamma-ray bursts (GRBs) provide an opportunity for constraining the emission mechanisms and geometries (e.g., of the accretion disk, jet, magnetic field, etc.) in the sources. For photons in the soft (MeV) gamma-ray band, Compton scattering is the most likely interaction to occur in detectors. Compton telescopes detect multiple interactions from individual incoming photons, allowing for scattering angles to be measured. After many photons are detected from a source, the distribution of azimuthal angles provides polarization information. While the standard method relies on binning the photons to produce and fit an azimuthal scattering angle distribution, improved polarization sensitivity is obtained by using additional information to more accurately weight each event's contribution to the likelihood statistic. In this chapter, we describe the Compton Spectrometer and Imager (COSI) and its capabilities for polarization measurements. We also describe the maximum likelihood technique, its application to COSI data analysis, and plans for its future use.


title: Gravitational Wave Signals from Finite Size Effects in Spinning Binary Inspirals Including Parity Violating Constituents

abstract: We generalize the world line EFT formalism to account for parity violating finite size effects. Results are presented for potentials and radiating moments of a binary inspiral for the parity conserving sector, and agreement is found with, previous calculations. Furthermore, we generate new results in this sector, calculating the current quadrupole moment induced by finite size gravitomagnetic effects. We also present novel results for parity violating sources, which might be due to beyond standard model physics, and show that they generate GW signals with the unique signature that the current-moment appears at 0.5PN order earlier relative to the mass-moment in the PN expansion. Parity violation also induces a new type of potential, which is proportional to the $\textbf{S}\cdot \textbf{r}$. Finally, we present new results for the dissipative force for parity violating constituents, which leads to the curious signature of a force normal to the orbit.


title: Spin chain on a metallic surface: Dissipation-induced order vs. Kondo entanglement

abstract: We explore the physics of a spin-1/2 Heisenberg chain with Kondo interaction, $J_k$, to a two-dimensional electron gas. At weak $J_k$ the problem maps onto a Heisenberg chain locally coupled to a dissipative Ohmic bath. At the decoupled fixed point, the dissipation is a marginally relevant perturbation and drives long-range antiferromagnetic order along the chain. In the dynamical spin structure factor we observe a quadratic low-energy dispersion akin to Landau-damped Goldstone modes. At large $J_k$ Kondo screening dominates, and the spin correlations of the chain inherit the power law of the host metal, akin to a paramagnetic heavy Fermi liquid. In both phases we observe heavy bands near the Fermi energy in the composite-fermion spectral function. Our results, obtained from auxiliary-field quantum Monte Carlo simulations, provide a unique negative-sign-free realization of a quantum transition between an antiferromagnetic metal and a heavy-fermion metal. We discuss the relevance of our results in the context of scanning tunneling spectroscopy experiments of magnetic adatom chains on metallic surfaces.


title: A misfired outburst in the neutron star X-ray binary Centaurus X-4

abstract: We report on a long-term optical monitoring of the neutron star X-ray binary Centaurus X-4 performed during the last 13.5 years. This source has been in quiescence since its outburst in 1979. Our monitoring reveals the overall evolution of the accretion disc; we detect short-duration flares, likely originating also in the disc, superimposed with a small-amplitude (< 0.1 mag) ellipsoidal modulation from the companion star due to geometrical effects. A long-term (~2300 days) downward trend, followed by a shorter (~1000 days) upward one, is observed in the disc light curve. Such a rise in the optical has been observed for other X-ray binaries preceding outbursts, as predicted by the disc instability model. For Cen X-4, the rise of the optical flux proceeded for ~3 years, and culminated in a flux increase at all wavelengths (optical-UV-X-rays) at the end of 2020. This increase faded after ~2 weeks, without giving rise to a full outburst. We suggest that the propagation of an inside-out heating front was ignited due to a partial ionization of hydrogen in the inner disc. The propagation might have stalled soon after the ignition due to the increasing surface density in the disc that the front encountered while propagating outwards. The stall was likely eased by the low level irradiation of the outer regions of the large accretion disc, as shown by the slope of the optical/X-ray correlation, suggesting that irradiation does not play a strong role in the optical, compared to other sources of emission.


title: Conformal Dirac-Einstein equations on manifolds with boundary

abstract: In this paper we study Dirac-Einstein equations on manifolds with boundary, restricted to a conformal class with constant boundary volume, under chiral bag boundary conditions for the Dirac operator. We characterize the bubbling phenomenon, also classifying ground state bubbles. Finally, we prove an Aubin-type inequality and a related existence result.


title: Truth Serum: Poisoning Machine Learning Models to Reveal Their Secrets

abstract: We introduce a new class of attacks on machine learning models. We show that an adversary who can poison a training dataset can cause models trained on this dataset to leak significant private details of training points belonging to other parties. Our active inference attacks connect two independent lines of work targeting the integrity and privacy of machine learning training data. Our attacks are effective across membership inference, attribute inference, and data extraction. For example, our targeted attacks can poison <0.1% of the training dataset to boost the performance of inference attacks by 1 to 2 orders of magnitude. Further, an adversary who controls a significant fraction of the training data (e.g., 50%) can launch untargeted attacks that enable 8x more precise inference on all other users' otherwise-private data points. Our results cast doubts on the relevance of cryptographic privacy guarantees in multiparty computation protocols for machine learning, if parties can arbitrarily select their share of training data.


title: Self-interacting dark baryons

abstract: Using results from lattice QCD, it is possible to quantitatively design models of dark baryons leading to velocity-dependent self-interaction cross sections that match the values needed for solving small-scale structure problems of standard cold dark matter. However it is not obvious that the main dark matter component in such models will be nucleons rather than large nuclei, or dark pions or atoms, whose scattering properties would be different. We first identify the parameters of a dark SU(3) sector analogous to QCD -- the confinement scale $\Lambda$ and pion mass $m_\pi$ -- needed to reproduce desired self-interaction cross sections. Then we show that these values can generically be compatible with the absence of a sufficiently stable deuteron bound state, and hence leading to no heavier nuclei, thus establishing the consistency of the scenario for self-interacting dark nucleons. The range of dark photon masses needed to avoid dominant pion or atomic dark matter is determined, as well as allowed values for the kinetic mixing parameter. The dark proton might be detected directly in future searches, by dark photon exchange.


title: Parallel Proof-of-Work with Concrete Bounds

abstract: Authorization is challenging in distributed systems that cannot rely on the identification of nodes. Proof-of-work offers an alternative gate-keeping mechanism, but its probabilistic nature is incompatible with conventional security definitions. Recent related work establishes concrete bounds for the failure probability of Bitcoin's sequential proof-of-work mechanism. We propose a family of state replication protocols using parallel proof-of-work. Our bottom-up design from an agreement sub-protocol allows us to give concrete bounds for the failure probability in adversarial synchronous networks. After the typical interval of 10 minutes, parallel proof-of-work offers two orders of magnitude more security than sequential proof-of-work. This means that state updates can be sufficiently secure to support commits after one block (i.e., after 10 minutes), removing the risk of double-spending in many applications. We offer guidance on the optimal choice of parameters for a wide range of network and attacker assumptions. Simulations show that the proposed construction is robust against violations of design assumptions.


title: Curiosity Driven Self-supervised Tactile Exploration of Unknown Objects

abstract: Intricate behaviors an organism can exhibit is predicated on its ability to sense and effectively interpret complexities of its surroundings. Relevant information is often distributed between multiple modalities, and requires the organism to exhibit information assimilation capabilities in addition to information seeking behaviors. While biological beings leverage multiple sensing modalities for decision making, current robots are overly reliant on visual inputs. In this work, we want to augment our robots with the ability to leverage the (relatively under-explored) modality of touch. To focus our investigation, we study the problem of scene reconstruction where touch is the only available sensing modality. We present Tactile Slam (tSLAM) -- which prepares an agent to acquire information seeking behavior and use implicit understanding of common household items to reconstruct the geometric details of the object under exploration. Using the anthropomorphic `ADROIT' hand, we demonstrate that tSLAM is highly effective in reconstructing objects of varying complexities within 6 seconds of interactions. We also established the generality of tSLAM by training only on 3D Warehouse objects and testing on ContactDB objects.


title: A Statistical Decision-Theoretical Perspective on the Two-Stage Approach to Parameter Estimation

abstract: One of the most important problems in system identification and statistics is how to estimate the unknown parameters of a given model. Optimization methods and specialized procedures, such as Empirical Minimization (EM) can be used in case the likelihood function can be computed. For situations where one can only simulate from a parametric model, but the likelihood is difficult or impossible to evaluate, a technique known as the Two-Stage (TS) Approach can be applied to obtain reliable parametric estimates. Unfortunately, there is currently a lack of theoretical justification for TS. In this paper, we propose a statistical decision-theoretical derivation of TS, which leads to Bayesian and Minimax estimators. We also show how to apply the TS approach on models for independent and identically distributed samples, by computing quantiles of the data as a first step, and using a linear function as the second stage. The proposed method is illustrated via numerical simulations.


title: Higher-order interactions promote chimera states

abstract: Since the discovery of chimera states, the presence of a nonzero phase lag parameter turns out to be an essential attribute for the emergence of chimeras in a nonlocally coupled identical Kuramoto phase oscillators' network with pairwise interactions. In this letter, we report the emergence of chimeras without phase lag in nonlocally coupled identical Kuramoto network owing to the introduction of non-pairwise interactions. The influence of added nonlinearity in the coupled system dynamics in the form of simplicial complexes mitigates the requisite of a nonzero phase lag parameter for the emergence of chimera states. Chimera states stimulated by the reciprocity of the pairwise and non-pairwise interaction strengths and their multistable nature are characterized with appropriate measures and are demonstrated in the parameter spaces.


title: Optimal Resource Scheduling and Allocation under Allowable Over-Scheduling

abstract: This paper studies optimal scheduling and resource allocation under allowable over-scheduling. Formulating an optimisation problem where over-scheduling is embedded, we derive an optimal solution that can be implemented by means of a new additive increase multiplicative decrease (AIMD) algorithm. After describing the AIMD-like scheduling mechanism as a switching system, we show convergence of the scheme, based on the joint spectral radius of symmetric matrices, and propose two methods for fitting an optimal AIMD tuning to the optimal solution derived. Finally, we demonstrate the overall optimal design strategy via an illustrative example.


title: Lattice QCD and the Computational Frontier

abstract: The search for new physics requires a joint experimental and theoretical effort. Lattice QCD is already an essential tool for obtaining precise model-free theoretical predictions of the hadronic processes underlying many key experimental searches, such as those involving heavy flavor physics, the anomalous magnetic moment of the muon, nucleon-neutrino scattering, and rare, second-order electroweak processes. As experimental measurements become more precise over the next decade, lattice QCD will play an increasing role in providing the needed matching theoretical precision. Achieving the needed precision requires simulations with lattices with substantially increased resolution. As we push to finer lattice spacing we encounter an array of new challenges. They include algorithmic and software-engineering challenges, challenges in computer technology and design, and challenges in maintaining the necessary human resources. In this white paper we describe those challenges and discuss ways they are being dealt with. Overcoming them is key to supporting the community effort required to deliver the needed theoretical support for experiments in the coming decade.


title: Integrating Biological Knowledge in Kernel-Based Analyses of Environmental Mixtures and Health

abstract: A key goal of environmental health research is to assess the risk posed by mixtures of pollutants. As epidemiologic studies of mixtures can be expensive to conduct, it behooves researchers to incorporate prior knowledge about mixtures into their analyses. This work extends the Bayesian multiple index model (BMIM), which assumes the exposure-response function is a non-parametric function of a set of linear combinations of pollutants formed with a set of exposure-specific weights. The framework is attractive because it combines the flexibility of response-surface methods with the interpretability of linear index models. We propose three strategies to incorporate prior toxicological knowledge into construction of indices in a BMIM: (a) constraining index weights, (b) structuring index weights by exposure transformations, and (c) placing informative priors on the index weights. We propose a novel prior specification that combines spike-and-slab variable selection with informative Dirichlet distribution based on relative potency factors often derived from previous toxicological studies. In simulations we show that the proposed priors improve inferences when prior information is correct and can protect against misspecification suffered by naive toxicological models when prior information is incorrect. Moreover, different strategies may be mixed-and-matched for different indices to suit available information (or lack thereof). We demonstrate the proposed methods on an analysis of data from the National Health and Nutrition Examination Survey and incorporate prior information on relative chemical potencies obtained from toxic equivalency factors available in the literature.


title: Diffusion of photo-excited holes in viscous electron fluid

abstract: The diffusion of photo-generated holes is studied in a high-mobility mesoscopic GaAs\ channel where electrons exhibit hydrodynamic properties. It is shown that the injection of holes into such an electron system leads to the formation of a hydrodynamic three-component mixture consisted of electrons and photo-generated heavy and light holes. The obtained results are analyzed within the framework of ambipolar diffusion, which reveals characteristics of a viscous flow. Both hole types exhibit similar hydrodynamic characteristics. In such a way the diffusion lengths, ambipolar diffusion coefficient and the effective viscosity of the electron-hole system are determined.


title: Efficient algorithms for computing ground states of the 2D random-field Ising model

abstract: We investigate the application of graph-cut methods for the study of the critical behaviour of the two-dimensional random-field Ising model. We focus on exact ground-state calculations, crossing the phase boundary of the model at zero temperature and varying the disorder strength. For this purpose we employ two different minimum-cut--maximum-flow algorithms, one of augmenting-path and another of push-relabel style. We implement these approaches for the square and triangular lattice problems and compare their computational efficiency.


title: Efficient Active Learning with Abstention

abstract: The goal of active learning is to achieve the same accuracy achievable by passive learning, while using much fewer labels. Exponential savings in terms of label complexity have been proved in very special cases, but fundamental lower bounds show that such improvements are impossible in general. This suggests a need to explore alternative goals for active learning. Learning with abstention is one such alternative. In this setting, the active learning algorithm may abstain from prediction and incur an error that is marginally smaller than random guessing. We develop the first computationally efficient active learning algorithm with abstention. Our algorithm provably achieves $\mathsf{polylog}(\frac{1}{\varepsilon})$ label complexity, without any low noise conditions. Such performance guarantee reduces the label complexity by an exponential factor, relative to passive learning and active learning that is not allowed to abstain. Furthermore, our algorithm is guaranteed to only abstain on hard examples (where the true label distribution is close to a fair coin), a novel property we term \emph{proper abstention} that also leads to a host of other desirable characteristics (e.g., recovering minimax guarantees in the standard setting, and avoiding the undesirable ``noise-seeking'' behavior often seen in active learning). We also provide novel extensions of our algorithm that achieve \emph{constant} label complexity and deal with model misspecification.


title: Corrections to scaling in geometrical clusters of the 2D Ising model

abstract: We study the scaling of the average cluster size and percolation strength of geometrical clusters for the two-dimensional Ising model. By means of Monte Carlo simulations and a finite-size scaling analysis we discuss the appearance of corrections to scaling for different definitions of cluster sets. We find that including all percolating clusters, or excluding only clusters that percolate in one but not the other direction, leads to smaller corrections to scaling for the average cluster size as compared to the other definitions considered. The percolation strength is less sensitive to the definitions used.


title: dS$\boldsymbol{_2}$ as excitation of AdS$\boldsymbol{_2}$

abstract: We introduce a family of 2D dilaton gravity models with state-dependent constant curvature so that dS$_2$ emerges as an excitation of AdS$_2$. Curiously, the strong coupling region corresponds to the asymptotic region geometrically. Apart from these key differences, many features resemble the Almheiri--Polchinski model. We discuss perturbative and non-perturbative thermodynamical stability, bubble nucleation through matter shockwaves, and semiclassical backreaction effects. In some of these models, we find that low temperatures are dominated by AdS$_2$ but high temperatures are dominated by dS$_2$, concurrent with a recent proposal by Susskind.


title: Geometric numerical methods for Lie systems and their application in optimal control

abstract: A Lie system is a non-autonomous system of first-order ordinary differential equations whose general solution can be written via an autonomous function, a so-called (nonlinear) superposition rule of a finite number of particular solutions and some parameters to be related to initial conditions. Even if the superposition rules for some Lie systems are known, the explicit analytic expression of their solutions frequently is not. This is why this article focuses on a novel geometric attempt to integrate Lie systems analytically and numerically. We focus on two families of methods: those based on Magnus expansions and the Runge-Kutta-Munthe-Kaas method, which are here adapted to the geometric properties of Lie systems. To illustrate the accuracy of our techniques we propose examples based on the SL$(n,\mathbb{R})$ Lie group, which plays a very relevant role in mechanics. In particular, we depict an optimal control problem for a vehicle with quadratic cost function. Particular numerical solutions of the studied examples are given.


title: Multicanonical simulations of the 2D spin-$1$ Baxter-Wu model in a crystal field

abstract: We investigate aspects of universality in the two-dimensional (2D) spin-$1$ Baxter-Wu model in a crystal field $\Delta$ using a parallel version of the multicanonical algorithm employed at constant temperature $T$. A detailed finite-size scaling analysis in the continuous regime of the $\Delta-T$ phase diagram of the model indicates that the transition belongs to the universality class of the $4$-state Potts model. The presence of first-order-like finite-size effects that become more pronounced as one approaches the pentacritical point of the model is highlighted and discussed.


title: Improving Biomarker Based HIV Incidence Estimation in the Treatment Era

abstract: Estimating HIV-1 incidence using biomarker assays in cross-sectional surveys is important for understanding the HIV pandemic. However, the utility of these estimates has been limited by uncertainty about what input parameters to use for False Recency Rate (FRR) and Mean Duration of Recent Infection (MDRI) after applying recent infection testing algorithm (RITA). This article shows how testing and diagnosis in a population reduce both FRR and MDRI compared to a treatment-naive population. Using self-reported testing history, a new method is proposed for calculating appropriate context-specific estimates of FRR and MDRI. The result of this is a new formula for incidence that depends only on reference FRR and MDRI parameters derived in an undiagnosed, treatment-naive, non-elite controller, non-AIDS-progressed population.


title: AKF-SR: Adaptive Kalman Filtering-based Successor Representation

abstract: Recent studies in neuroscience suggest that Successor Representation (SR)-based models provide adaptation to changes in the goal locations or reward function faster than model-free algorithms, together with lower computational cost compared to that of model-based algorithms. However, it is not known how such representation might help animals to manage uncertainty in their decision-making. Existing methods for SR learning do not capture uncertainty about the estimated SR. In order to address this issue, the paper presents a Kalman filter-based SR framework, referred to as Adaptive Kalman Filtering-based Successor Representation (AKF-SR). First, Kalman temporal difference approach, which is a combination of the Kalman filter and the temporal difference method, is used within the AKF-SR framework to cast the SR learning procedure into a filtering problem to benefit from the uncertainty estimation of the SR, and also decreases in memory requirement and sensitivity to model's parameters in comparison to deep neural network-based algorithms. An adaptive Kalman filtering approach is then applied within the proposed AKF-SR framework in order to tune the measurement noise covariance and measurement mapping function of Kalman filter as the most important parameters affecting the filter's performance. Moreover, an active learning method that exploits the estimated uncertainty of the SR to form the behaviour policy leading to more visits to less certain values is proposed to improve the overall performance of an agent in terms of received rewards while interacting with its environment.


title: Leakage Localization in Water Distribution Networks: A Model-Based Approach

abstract: The paper studies the problem of leakage localization in water distribution networks. For the case of a single pipe that suffers from a single leak, by taking recourse to pressure and flow measurements, and assuming those are noiseless, we provide a closed-form expression for leak localization, leak exponent and leak constant. For the aforementioned setting, but with noisy pressure and flow measurements, an expression for estimating the location of the leak is provided. Finally, assuming the existence of a single leak, for a network comprising of more than one pipe and assuming that the network has a tree structure, we provide a systematic procedure for determining the leak location, the leak exponent, and the leak constant


title: Phase Field Crystal model for particles with n-fold rotational symmetry in two dimensions

abstract: We introduce a Phase Field Crystal (PFC) model for particles with n-fold rotational symmetry in two dimensions. Our approach is based on a free energy functional that depends on the reduced one-particle density, the strength of the orientation, and the direction of the orientation, where all these order parameters depend on the position. The functional is constructed such that for particles with axial symmetry (i. e. n = 2) the PFC model for liquid crystals as introduced by H. L\"owen [J. Phys.: Condens. Matter 22, 364105 (2010)] is recovered. We discuss the stability of the functional and explore phases that occur for $1 \leq n \leq 6$. In addition to isotropic, nematic, stripe, and triangular order, we also observe cluster crystals with square, rhombic, honeycomb, and even quasicrystalline symmetry. The n-fold symmetry of the particles corresponds to the one that can be realized for colloids with symmetrically arranged patches. We explain how both, repulsive as well as attractive patches, are described in our model.


title: Digitizing Historical Balance Sheet Data: A Practitioner's Guide

abstract: This paper discusses how to successfully digitize large-scale historical micro-data by augmenting optical character recognition (OCR) engines with pre- and post-processing methods. Although OCR software has improved dramatically in recent years due to improvements in machine learning, off-the-shelf OCR applications still present high error rates which limit their applications for accurate extraction of structured information. Complementing OCR with additional methods can however dramatically increase its success rate, making it a powerful and cost-efficient tool for economic historians. This paper showcases these methods and explains why they are useful. We apply them against two large balance sheet datasets and introduce quipucamayoc, a Python package containing these methods in a unified framework.


title: Ultra-high terahertz index in deep subwavelength coupled bi-layer free-standing flexible metamaterials

abstract: We report extensive enhancement in the refractive index of artificially designed metamaterials by exploiting the deep subwavelength coupling in a free-standing, thin-film metal-dielectric-metal checkboard structure. A record high refractive index of 77.02+43.22i is obtained at terahertz frequencies. The detailed investigations reveal that the enhancement of the effective refractive index of the structure via deep subwavelength coupling in a bilayer design is governed by a power law, which is an effective and simpler approach to design high index metamaterial. The approach relies on deep subwavelength coupling to obtain extremely high refractive index values that can lead to many practical applications in the field of imaging, lithography, design of delay lines and interferometers.


title: Distributed Robust Geocast Multicast Routing for Inter-Vehicle Communication

abstract: Numerous protocols for geocast have been proposed in literature. It has been shown that explicit route setup approaches perform poorly with VANETs due to limited route lifetime and frequent network fragmentation. The broadcast based approaches have considerable redundancy and add significantly to the overhead of the protocol. A completely distributed and robust geocast approach is presented in this paper, that is resilient to frequent topology changes and network fragmentation. A distance-based backoff algorithm is used to reduce the number of hops and a novel mechanism to reduce redundant broadcasts is introduced. The performance of the proposed protocol is evaluated for various scenarios and compared with simple flooding and a protocol based on explicit route setup.


title: Inequalities for higher order differences of the logarithm of the overpartition function and a problem of Wang-Xie-Zhang

abstract: Let $\overline{p}(n)$ denote the overpartition function. In this paper, our primary goal is to study the asymptotic behavior of the finite differences of the logarithm of the overpartition function, i.e., $(-1)^{r-1}\Delta^r \log \p(n)$, by studying the inequality of the following form $$\log \Bigl(1+\dfrac{C(r)}{n^{r-1/2}}-\dfrac{C_1(r)}{n^{r}}\Bigr)<(-1)^{r-1}\Delta^r \log \p(n) <\log \Bigl(1+\dfrac{C(r)}{n^{r-1/2}}\Bigr)\ \text{for}\ n \geq N(r),$$ where $C(r), C_1(r), \text{and}\ N(r)$ are computable constants depending on the positive integer $r$, determined explicitly. This solves a problem posed by Wang, Xie and Zhang in the context of searching for a better lower bound of $(-1)^{r-1}\Delta^r \log \p(n)$ than $0$. By settling the problem, we are able to show that \begin{equation*} \lim_{n\rightarrow \infty}(-1)^{r-1}\Delta^r \log \p(n) =\dfrac{\pi}{2}\Bigl(\dfrac{1}{2}\Bigr)_{r-1}n^{\frac{1}{2}-r}. \end{equation*}


title: SELFIES and the future of molecular string representations

abstract: Artificial intelligence (AI) and machine learning (ML) are expanding in popularity for broad applications to challenging tasks in chemistry and materials science. Examples include the prediction of properties, the discovery of new reaction pathways, or the design of new molecules. The machine needs to read and write fluently in a chemical language for each of these tasks. Strings are a common tool to represent molecular graphs, and the most popular molecular string representation, SMILES, has powered cheminformatics since the late 1980s. However, in the context of AI and ML in chemistry, SMILES has several shortcomings -- most pertinently, most combinations of symbols lead to invalid results with no valid chemical interpretation. To overcome this issue, a new language for molecules was introduced in 2020 that guarantees 100\% robustness: SELFIES (SELF-referencIng Embedded Strings). SELFIES has since simplified and enabled numerous new applications in chemistry. In this manuscript, we look to the future and discuss molecular string representations, along with their respective opportunities and challenges. We propose 16 concrete Future Projects for robust molecular representations. These involve the extension toward new chemical domains, exciting questions at the interface of AI and robust languages and interpretability for both humans and machines. We hope that these proposals will inspire several follow-up works exploiting the full potential of molecular string representations for the future of AI in chemistry and materials science.


title: ElectAnon: A Blockchain-Based, Anonymous, Robust and Scalable Ranked-Choice Voting Protocol

abstract: Remote voting has become more critical in recent years, especially after the Covid-19 outbreak. Blockchain technology and its benefits like decentralization, security, and transparency have encouraged remote voting systems to use blockchains. Analysis of existing solutions reveals that anonymity, robustness, and scalability are common problems in blockchain-based election systems. In this work, we propose ElectAnon, a blockchain-based, ranked-choice election protocol focusing on anonymity, robustness, and scalability. ElectAnon achieves anonymity by enabling voters to cast their votes via zero-knowledge proofs anonymously. Robustness is realized by removing the direct control of the authorities in the voting process by using timed-state machines. Results show that ElectAnon is scalable amongst existing works as it reduces the gas consumption up to 89% compared to previous works. The proposed protocol includes a candidate proposal system and swappable tallying libraries. An extension is also proposed to minimize the trust assumption on election authorities. Our code is available on https://github.com/ceyonur/electanon.


title: The Multilinear Spherical Maximal Function in one dimension

abstract: In dimension $n=1$ we obtain $L^{p_1}(\mathbb R) \times\dots\times L^{p_m}(\mathbb R)$ to $L^p(\mathbb R)$ boundedness for the multilinear spherical maximal function in the largest possible open set of indices and we provide counterexamples that indicate the optimality of our results.


title: Lyapunov based Stochastic Stability of Human-Machine Interaction: A Quantum Decision System Approach

abstract: In mathematical psychology, decision makers are modeled using the Lindbladian equations from quantum mechanics to capture important human-centric features such as order effects and violation of the sure thing principle. We consider human-machine interaction involving a quantum decision maker (human) and a controller (machine). Given a sequence of human decisions over time, how can the controller dynamically provide input messages to adapt these decisions so as to converge to a specific decision? We show via novel stochastic Lyapunov arguments how the Lindbladian dynamics of the quantum decision maker can be controlled to converge to a specific decision asymptotically. Our methodology yields a useful mathematical framework for human-sensor decision making. The stochastic Lyapunov results are also of independent interest as they generalize recent results in the literature.


title: Rogue waves in discrete-time quantum walks

abstract: Rogue waves are rapid and unpredictable events of exceptional amplitude reported in various fields, such as oceanography and optics, with much of the interest being targeted towards their physical origins and likelihood of occurrence. Here, we use the all-round framework of discrete-time quantum walks to study the onset of those events due to a random phase modulation, unveiling its long-tailed statistics, distribution profile, and dependence upon the degree of randomness. We find that those rogue waves belong the Gumbel family of extreme value distributions.


title: Data-augmented cross-lingual synthesis in a teacher-student framework

abstract: Cross-lingual synthesis can be defined as the task of letting a speaker generate fluent synthetic speech in another language. This is a challenging task, and resulting speech can suffer from reduced naturalness, accented speech, and/or loss of essential voice characteristics. Previous research shows that many models appear to have insufficient generalization capabilities to perform well on every of these cross-lingual aspects. To overcome these generalization problems, we propose to apply the teacher-student paradigm to cross-lingual synthesis. While a teacher model is commonly used to produce teacher forced data, we propose to also use it to produce augmented data of unseen speaker-language pairs, where the aim is to retain essential speaker characteristics. Both sets of data are then used for student model training, which is trained to retain the naturalness and prosodic variation present in the teacher forced data, while learning the speaker identity from the augmented data. Some modifications to the student model are proposed to make the separation of teacher forced and augmented data more straightforward. Results show that the proposed approach improves the retention of speaker characteristics in the speech, while managing to retain high levels of naturalness and prosodic variation.


title: SimPO: Simultaneous Prediction and Optimization

abstract: Many machine learning (ML) models are integrated within the context of a larger system as part of a key component for decision making processes. Concretely, predictive models are often employed in estimating the parameters for the input values that are utilized for optimization models as isolated processes. Traditionally, the predictive models are built first, then the model outputs are used to generate decision values separately. However, it is often the case that the prediction values that are trained independently of the optimization process produce sub-optimal solutions. In this paper, we propose a formulation for the Simultaneous Prediction and Optimization (SimPO) framework. This framework introduces the use of a joint weighted loss of a decision-driven predictive ML model and an optimization objective function, which is optimized end-to-end directly through gradient-based methods.


title: On the Generalised Ricci Solitons and Sasakian Manifolds

abstract: In this note, we find a necessary condition on odd-dimensional Riemannian manifolds under which both of Sasakian structure and the generalised Ricci soliton equation are satisfied, and we give some examples.


title: Simplistic approach for Determining Center Shift of M\"ossbauer Spectrum: exemplified on two cases

abstract: A very simple method for determining the center (isomer) shift, CS, of a M\"ossbauer spectrum is outlined. Its applicability is demonstrated on two examples viz. pyrite and a ternary sigma-phase Fe-Cr-Ni compound. Sets of the spectra recorded in the temperature interval of 78-295 K for the former, and 5-293 K for the latter were analyzed with the simple and a commonly used methods. CS(T)-values obtained with both ways of the fitting procedures were analyzed in terms of the Debye model. The determined therefrom values of the Debye temperatures agree within the error limit with each other proving thereby that this very simple method gives correct values of CS.


title: Importance of Different Temporal Modulations of Speech: A Tale of Two Perspectives

abstract: How important are different temporal speech modulations for speech recognition? We answer this question from two complementary perspectives. Firstly, we quantify the amount of phonetic \textit{information} in the modulation spectrum of speech by computing the mutual information between temporal modulations with frame-wise phoneme labels. Looking from another perspective, we ask - which speech modulations an Automatic Speech Recognition (ASR) system prefers for its operation. Data-driven weights are learned over the modulation spectrum and optimized for an end-to-end ASR task. Both methods unanimously agree that speech information is mostly contained in slow modulation. Maximum mutual information occurs around 3-6 Hz which also happens to be the range of modulations most preferred by the ASR. In addition, we show that the incorporation of this knowledge into ASRs significantly reduces their dependency on the amount of training data.


title: Jordan types of triangular matrices over a finite field

abstract: Let $\lambda$ be a partition of an integer $n$ and ${\mathbb F}_q$ be a finite field of order $q$. Let $P_\lambda(q)$ be the number of strictly upper triangular $n\times n$ matrices of the Jordan type $\lambda$. It is known that the polynomial $P_\lambda$ has a tendency to be divisible by high powers of $q$ and $Q=q-1$, and we put $P_\lambda(q)=q^{d(\lambda)}Q^{e(\lambda)}R_\lambda(q)$, where $R_\lambda(0)\neq0$ and $R_\lambda(1)\neq0$. In this article, we study the polynomials $P_\lambda(q)$ and $R_\lambda(q)$. Our main results: an explicit formula for $d(\lambda)$ (an explicit formula for $e(\lambda)$ is known, see Proposition 3.3 below), a recursive formula for $R_\lambda(q)$ (a similar formula for $P_\lambda(q)$ is known, see Proposition 3.1 below), the stabilization of $R_\lambda$ with respect to extending $\lambda$ by adding strings of 1's, and an explicit formula for the limit series $R_{\lambda1^\infty}$. Our studies are motivated by projected applications to the orbit method in the representation theory of nilpotent algebraic groups over finite fields.


title: rfPhen2Gen: A machine learning based association study of brain imaging phenotypes to genotypes

abstract: Imaging genetic studies aim to find associations between genetic variants and imaging quantitative traits. Traditional genome-wide association studies (GWAS) are based on univariate statistical tests, but when multiple traits are analyzed together they suffer from a multiple-testing problem and from not taking into account correlations among the traits. An alternative approach to multi-trait GWAS is to reverse the functional relation between genotypes and traits, by fitting a multivariate regression model to predict genotypes from multiple traits simultaneously. However, current reverse genotype prediction approaches are mostly based on linear models. Here, we evaluated random forest regression (RFR) as a method to predict SNPs from imaging QTs and identify biologically relevant associations. We learned machine learning models to predict 518,484 SNPs using 56 brain imaging QTs. We observed that genotype regression error is a better indicator of permutation p-value significance than genotype classification accuracy. SNPs within the known Alzheimer disease (AD) risk gene APOE had lowest RMSE for lasso and random forest, but not ridge regression. Moreover, random forests identified additional SNPs that were not prioritized by the linear models but are known to be associated with brain-related disorders. Feature selection identified well-known brain regions associated with AD,like the hippocampus and amygdala, as important predictors of the most significant SNPs. In summary, our results indicate that non-linear methods like random forests may offer additional insights into phenotype-genotype associations compared to traditional linear multi-variate GWAS methods.


title: Automatic Classification of Alzheimer's Disease using brain MRI data and deep Convolutional Neural Networks

abstract: Alzheimer's disease (AD) is one of the most common public health issues the world is facing today. This disease has a high prevalence primarily in the elderly accompanying memory loss and cognitive decline. AD detection is a challenging task which many authors have developed numerous computerized automatic diagnosis systems utilizing neuroimaging and other clinical data. MRI scans provide high-intensity visible features, making these scans the most widely used brain imaging technique. In recent years deep learning has achieved leading success in medical image analysis. But a relatively little investigation has been done to apply deep learning techniques for the brain MRI classification. This paper explores the construction of several deep learning architectures evaluated on brain MRI images and segmented images. The idea behind segmented images investigates the influence of image segmentation step on deep learning classification. The image processing presented a pipeline consisting of pre-processing to enhance the MRI scans and post-processing consisting of a segmentation method for segmenting the brain tissues. The results show that the processed images achieved a better accuracy in the binary classification of AD vs. CN (Cognitively Normal) across four different architectures. ResNet architecture resulted in the highest prediction accuracy amongst the other architectures (90.83% for the original brain images and 93.50% for the processed images).


title: A One-Dimensional Model for Investigating Scale-separated Approaches to the Interaction of Oceanic Internal Waves

abstract: High-frequency wave propagation in near-inertial wave shear has been considered fundamental in setting the spectral character of the oceanic internal wave continuum and for transporting energy to wave-breaking. We compare idealized ray tracing numerical results with metrics derived using a wave turbulence derivation for the kinetic equation and a path integral to study such scale-separated interactions. At small inertial wave amplitudes, all three provide consistent descriptions for the mean drift of wavepackets in the spectral domain and dispersion about that mean drift. Extrapolating our results to the background internal wavefield over-predicts downscale energy transports by an order of magnitude. At oceanic amplitudes, however, the numerics support diminished transport and dispersion that coincide with the mean drift time scale becoming similar to the lagged correlation time scale. Despite this decrease, numerical estimates of downscale energy transfer are still significantly larger than oceanic derived metrics. Our results support replacing the long standing interpretive paradigm of extreme scale-separated interactions with a more nuanced slate of 'local' interactions in the kinetic equation.


title: KPZ physics and phase transition in a classical single random walker under continuous measurement

abstract: We introduce and study a new model consisting of a single classical random walker undergoing continuous monitoring at rate $\gamma$ on a discrete lattice. Although such a continuous measurement cannot affect physical observables, it has a non-trivial effect on the probability distribution of the random walker. At small $\gamma$, we show analytically that the time-evolution of the latter can be mapped to the Stochastic Heat Equation (SHE). In this limit, the width of the log probability thus follows a Family-Vicsek scaling law, $N^{\alpha}f(t/N^{\alpha/\beta})$, with roughness and growth exponents corresponding to the Kardar-Parisi-Zhang (KPZ) universality class, i.e $\alpha^{\rm{1D}}_{\rm{KPZ}}=1/2$ and $\beta^{\rm{1D}}_{\rm{KPZ}}=1/3$ respectively. When $\gamma$ is increased outside this regime, we find numerically in 1D a crossover from the KPZ class to a new universality class characterized by exponents $\alpha^{1\rm{D}}_{\text{M}}\approx 1$ and $\beta^{1\rm{D}}_{\text{M}}\approx 1.4$. In 3D, varying $\gamma$ beyond a critical value $\gamma^c_{\rm{M}}$ leads to a phase transition from a smooth phase that we identify as the Edwards-Wilkinson (EW) class to a new universality class with $\alpha^{3\rm{D}}_{\text{M}}\approx1$.


title: Numerical Solution of the Steady-State Network Flow Equations for a Non-Ideal Gas

abstract: We formulate a steady-state network flow problem for non-ideal gas that relates injection rates and nodal pressures in the network to flows in pipes. For this problem, we present and prove a theorem on uniqueness of generalized solution for a broad class of non-ideal pressure-density relations that satisfy a monotonicity property. Further, we develop a Newton-Raphson algorithm for numerical solution of the steady-state problem, which is made possible by a systematic non-dimensionalization of the equations. The developed algorithm has been extensively tested on benchmark instances and shown to converge robustly to a generalized solution. Previous results indicate that the steady-state network flow equations for an ideal gas are difficult to solve by the Newton-Raphson method because of its extreme sensitivity to the initial guess. In contrast, we find that non-dimensionalization of the steady-state problem is key to robust convergence of the Newton-Raphson method. We identify criteria based on the uniqueness of solutions under which the existence of a non-physical generalized solution found by a non-linear solver implies non-existence of a physical solution, i.e., infeasibility of the problem. Finally, we compare pressure and flow solutions based on ideal and non-ideal equations of state to demonstrate the need to apply the latter in practice. The solver developed in this article is open-source and is made available for both the academic and research communities as well as the industry.


title: Cyclic independence: Boolean and monotone

abstract: The present paper introduces a modified version of cyclic-monotone independence which originally arose in the context of random matrices, and also introduces its natural analogy called cyclic-Boolean independence. We investigate formulas for convolutions, limit theorems for sums of independent random variables, and also classify infinitely divisible distributions with respect to cyclic-Boolean convolution. Finally, we provide applications to the eigenvalues of the adjacency matrices of iterated star products of graphs and also iterated comb products of graphs.


title: On-line Estimation of Stability and Passivity Metrics

abstract: We consider the problem of on-line evaluation of critical characteristic parameters such as the L_2-gain (L2G), input feedforward passivity index (IFP) and output feedback passivity index (OFP) of non-linear systems using their input-output data. Typically, having an accurate measure of such "system indices" enables the application of systematic control design techniques. Moreover, if such system indices can efficiently be evaluated on-line, they can be exploited to device intelligent controller reconfiguration and fault-tolerant control techniques. However, the existing estimation methods of such system indices (i.e., L2G, IFP and OFP) are predominantly off-line, computationally inefficient, and require a large amount of actual or synthetically generated input-output trajectory data under some specific initial/terminal conditions. On the other hand, the existing on-line estimation methods take an averaging-based approach, which may be sub-optimal, computationally inefficient and susceptible to estimate saturation. In this paper, to overcome these challenges (in the on-line estimation of system indices), we establish and exploit several interesting theoretical results on a particular class of fractional function optimization problems. For comparison purposes, the details of an existing averaging-based approach are provided for the same on-line estimation problem. Finally, several numerical examples are discussed to demonstrate the proposed on-line estimation approach and to highlight our contributions.


title: When Artificial Parameter Evolution Gets Real: Particle Filtering for Time-Varying Parameter Estimation in Deterministic Dynamical Systems

abstract: Estimating and quantifying uncertainty in unknown system parameters from limited data remains a challenging inverse problem in a variety of real-world applications. While many approaches focus on estimating constant parameters, a subset of these problems includes time-varying parameters with unknown evolution models that often cannot be directly observed. This work develops a systematic particle filtering approach that reframes the idea behind artificial parameter evolution to estimate time-varying parameters in nonstationary inverse problems arising from deterministic dynamical systems. Focusing on systems modeled by ordinary differential equations, we present two particle filter algorithms for time-varying parameter estimation: one that relies on a fixed value for the noise variance of a parameter random walk; another that employs online estimation of the parameter evolution noise variance along with the time-varying parameter of interest. Several computed examples demonstrate the capability of the proposed algorithms in estimating time-varying parameters with different underlying functional forms and different relationships with the system states (i.e., additive vs. multiplicative).


title: The holographic entropy cone from marginal independence

abstract: The holographic entropy cone characterizes the relations between entanglement entropies for a spatial partitioning of the boundary spacetime of a holographic CFT in any state describing a classical bulk geometry. We argue that the holographic entropy cone, for an arbitrary number of parties, can be reconstructed from more fundamental data determined solely by subadditivity of quantum entropy. We formulate certain conjectures about graph models of holographic entanglement, for which we provide strong evidence, and rigorously prove that they all imply that such a reconstruction is possible. Our conjectures (except only for the weakest) further imply that the necessary data is remarkably simple. In essence, all one needs to know to reconstruct the holographic entropy cone, is a certain subset of the extreme rays of this simpler "subadditivity cone", namely those which can be realized in holography. This recasting of the bewildering entanglement structure of geometric states into primal building blocks paves the way to distilling the essence of holography for the emergence of a classical bulk spacetime.


title: Reachability Logic for Low-Level Programs

abstract: Automatic exploit generation is a relatively new area of research. Work in this area aims to automate the manual and labor intensive task of finding exploits in software. In this paper we present a novel program logic to support automatic exploit generation. We develop a program logic called Reachability Logic, which formally defines the relation between reachability of an assertion and the preconditions which allow them to occur. This relation is then used to calculate the search space of preconditions. We show that Reachability Logic is a powerful tool in automatically finding evidence that an assertion is reachable. We verify that the system works for small litmus tests, as well as real-world algorithms. An implementation has been developed, and the entire system is proven to be sound and complete in a theorem prover. This work represents an important step towards formally verified automatic exploit generation.


title: Efficient Maximal Coding Rate Reduction by Variational Forms

abstract: The principle of Maximal Coding Rate Reduction (MCR$^2$) has recently been proposed as a training objective for learning discriminative low-dimensional structures intrinsic to high-dimensional data to allow for more robust training than standard approaches, such as cross-entropy minimization. However, despite the advantages that have been shown for MCR$^2$ training, MCR$^2$ suffers from a significant computational cost due to the need to evaluate and differentiate a significant number of log-determinant terms that grows linearly with the number of classes. By taking advantage of variational forms of spectral functions of a matrix, we reformulate the MCR$^2$ objective to a form that can scale significantly without compromising training accuracy. Experiments in image classification demonstrate that our proposed formulation results in a significant speed up over optimizing the original MCR$^2$ objective directly and often results in higher quality learned representations. Further, our approach may be of independent interest in other models that require computation of log-determinant forms, such as in system identification or normalizing flow models.


title: The DESC Stellarator Code Suite Part III: Quasi-symmetry optimization

abstract: The DESC stellarator optimization code takes advantage of advanced numerical methods to search the full parameter space much faster than conventional tools. Only a single equilibrium solution is needed at each optimization step thanks to automatic differentiation, which efficiently provides exact derivative information. A Gauss-Newton trust-region optimization method uses second-order derivative information to take large steps in parameter space and converges rapidly. With just-in-time compilation and GPU portability, high-dimensional stellarator optimization runs take orders of magnitude less computation time with DESC compared to other approaches. This paper presents the theory of the DESC fixed-boundary local optimization algorithm along with demonstrations of how to easily implement it in the code. Example quasi-symmetry optimizations are shown and compared to results from conventional tools. Three different forms of quasi-symmetry objectives are available in DESC, and their relative advantages are discussed in detail. In the examples presented, the triple product formulation yields the best optimization results in terms of minimized computation time and particle transport. This paper concludes with an explanation of how the modular code suite can be extended to accommodate other types of optimization problems.


title: Sensitive label-free and compact ultrasonic sensor based on double silicon-on-insulator slot micro-ring resonators

abstract: We propose a new label-free ultrasonic sensor, which comprises a slot wave-guide and double silicon-on-insulator (SOI) slot micro-ring resonators. The all-optical sensors do not suffer from electromagnetic interference. We choose to integrate a silicon slot double micro-ring (SDMR) resonators in an acoustically resonant membrane. Optimization of the several key structural parameters is investigated to achieve the mode-field distributions of transmission spectrum based on Comsol Multiphysics software. Our numerical studies show that the proposed ultrasonic sensor offers higher sensitivity and a larger detection frequency range than conventional piezoelectric based ultrasound transducer. For a SDMR system with an area of 15um*30um, sensitivity as high as 2453.7mv/kpa, and over a bandwidth range of 1-150MHz. The sensitivity value is 36 times higher than that of single slot micro-ring ultrasonic sensor. The theoretical-factor of the SDMR can be approximately 1.24*10^6 with bending radius of 5um. The investigation on the SDMR system is a valuable exploration of the photo-acoustic microscopy for the ultra-high Q factor and large frequency range.


title: 4Weed Dataset: Annotated Imagery Weeds Dataset

abstract: Weeds are a major threat to crops and are responsible for reducing crop yield worldwide. To mitigate their negative effect, it is advantageous to accurately identify them early in the season to prevent their spread throughout the field. Traditionally, farmers rely on manually scouting fields and applying herbicides for different weeds. However, it is easy to confuse between crops with weeds during the early growth stages. Recently, deep learning-based weed identification has become popular as deep learning relies on convolutional neural networks that are capable of learning important distinguishable features between weeds and crops. However, training robust deep learning models requires access to large imagery datasets. Therefore, an early-season weeds dataset was acquired under field conditions. The dataset consists of 159 Cocklebur images, 139 Foxtail images, 170 Redroot Pigweed images and 150 Giant Ragweed images corresponding to four common weed species found in corn and soybean production systems.. Bounding box annotations were created for each image to prepare the dataset for training both image classification and object detection deep learning networks capable of accurately locating and identifying weeds within corn and soybean fields. (https://osf.io/w9v3j/)


title: Support-vector-machine with Bayesian optimization for lithofacies classification using elastic properties

abstract: We investigate an applicability of Bayesian-optimization (BO) to optimize hyperparameters associated with support-vector-machine (SVM) in order to classify facies using elastic properties derived from well data in the East Central Graben, UKCS. The cross-plot products of the field dataset appear to be successfully classified with non-linear boundaries. Although there are a few factors to be predetermined in the BO scheme such as an iteration number to deal with a trade-off between the prediction accuracy and the computational cost, this approach effectively reduces possible human subjectivity connected to the architecture of the SVM. Our proposed workflow might be beneficial in resource-exploration and development in terms of subsurface objective technical evaluations.


title: Quantitative analysis of diaphragm motion during fluoroscopic sniff test to assist in diagnosis of hemidiaphragm paralysis

abstract: The current imaging gold standard for detecting paradoxical diaphragm motion and diagnosing hemidiaphragm paralysis is to perform the fluoroscopic sniff test. The images are visually examined by an experienced radiologist, and if one hemidiaphragm ascends while the other descends, then it is described as paradoxical motion, which is highly suggestive of hemidiaphragm paralysis. However, diagnosis can be challenging because diaphragm motion during sniffing is fast, paradoxical motion can be subtle, and the analysis is based on a 2-dimensional projection of a 3-dimensional surface. This paper presents a case of chronic left hemidiaphragm elevation that was initially reported as mild paradoxical motion on fluoroscopy. After measuring the elevations of the diaphragms and modeling their temporal correlation using Gaussian process regression, the systematic trend of the hemidiaphragmatic motion along with its stochastic properties was determined. When analyzing the trajectories of the hemidiaphragms, no statistically significant paradoxical motion was detected. This could potentially change the prognosis if the patient was to consider diaphragm plication as treatment. The presented method provides a more objective analysis of hemidiaphragm motions and can potentially improve diagnostic accuracy.


title: Irreducibility of Chebyshev-Lissajous polynomials

abstract: We study certain kind of polynomials associated with Lissajous curves, called Chebyshev-Lissajous polynomials. We investigate their irreducibilities over the real numbers and complex numbers, thus comfirming two conjectures proposed by Merino.


title: Evolution of barchan dune interactions investigated by a downscaled water tunnel experiment: the temporal characteristics and a soliton-like behavior

abstract: This paper reports a downscaled water tunnel experiment to study the temporal characteristics of a double dune interaction system and the new pattern of dune interaction when the initial mass ratio of the two dunes is large. These topics are useful for a comprehensive understanding of the dune interaction system but were rarely covered before. The turnover time scale under dune interaction is defined, and its time averaged value is found to have a nonmonotonic relationship with the initial mass ratio. A nonmonotonic relationship is also found between the convexity of the downstream dune tip and the initial mass ratio. The stationary points of the two nonmonotonic curves above correspond to the same dune interaction pattern named 'exchange-chasing', which is considered indispensable in the classification map of dune interactions. The upstream dune acts as an energy transmitter between fluid flow and the downstream dune. A soliton-like behavior occurs when the downstream dune enlarges, where a small dune is detached from the downstream dune tip and gets passed by the upstream dune approximately without mass exchange. The activity of such temporary soliton is found to be negatively related with the initial dune spacing and positively related with the initial mass ratio.


title: Bifurcation of limit cycles in piecewise quadratic differential systems with an invariant straight line

abstract: We solve the center-focus problem in a class of piecewise quadratic polynomial differential systems with an invariant straight line. The separation curve is also a straight line which is not invariant. We provide families having at the origin a weak-foci of maximal order. In the continuous class, the cyclicity problem is also solved, being $3$ such maximal number. Moreover, for the discontinuous class but without sliding segment, we prove the existence of $7$ limit cycles of small amplitude.


title: Discoveries and Timing of Pulsars in NGC 6440

abstract: Using the MeerKAT radio telescope, a series of observations have been conducted to time the known pulsars and search for new pulsars in the globular cluster NGC 6440. As a result, two pulsars have been discovered, NGC 6440G and NGC 6440H, one of which is isolated and the other a non-eclipsing (at frequencies above 962 MHz) "Black Widow", with a very low mass companion (M$_{\rm c}$ > 0.006 M$_{\odot}$). It joins the other binary pulsars discovered so far in this cluster which all have low companion masses (M$_{\rm c}$ < 0.30 M$_{\odot}$). We present the results of long-term timing solutions obtained using data from both Green Bank and MeerKAT telescopes for these two new pulsars and an analysis of the pulsars NGC 6440C and NGC 6440D. For the isolated pulsar NGC 6440C, we searched for planets using a Markov Chain Monte Carlo technique. We find evidence for significant unmodelled variations but they cannot be well modelled as planets nor as part of a power-law red-noise process. Studies of the eclipses of the "Redback" pulsar NGC 6440D at two different frequency bands reveal a frequency dependence with longer and asymmetric eclipses at lower frequencies (962-1283 MHz).


title: A quantum learning approach based on Hidden Markov Models for failure scenarios generation

abstract: Finding the failure scenarios of a system is a very complex problem in the field of Probabilistic Safety Assessment (PSA). In order to solve this problem we will use the Hidden Quantum Markov Models (HQMMs) to create a generative model. Therefore, in this paper, we will study and compare the results of HQMMs and classical Hidden Markov Models HMM on a real datasets generated from real small systems in the field of PSA. As a quality metric we will use Description accuracy DA and we will show that the quantum approach gives better results compared with the classical approach, and we will give a strategy to identify the probable and no-probable failure scenarios of a system.


title: Speech and the n-Back task as a lens into depression. How combining both may allow us to isolate different core symptoms of depression

abstract: Embedded in any speech signal is a rich combination of cognitive, neuromuscular and physiological information. This richness makes speech a powerful signal in relation to a range of different health conditions, including major depressive disorders (MDD). One pivotal issue in speech-depression research is the assumption that depressive severity is the dominant measurable effect. However, given the heterogeneous clinical profile of MDD, it may actually be the case that speech alterations are more strongly associated with subsets of key depression symptoms. This paper presents strong evidence in support of this argument. First, we present a novel large, cross-sectional, multi-modal dataset collected at Thymia. We then present a set of machine learning experiments that demonstrate that combining speech with features from an n-Back working memory assessment improves classifier performance when predicting the popular eight-item Patient Health Questionnaire depression scale (PHQ-8). Finally, we present a set of experiments that highlight the association between different speech and n-Back markers at the PHQ-8 item level. Specifically, we observe that somatic and psychomotor symptoms are more strongly associated with n-Back performance scores, whilst the other items: anhedonia, depressed mood, change in appetite, feelings of worthlessness and trouble concentrating are more strongly associated with speech changes.


title: Investigating Top-$k$ White-Box and Transferable Black-box Attack

abstract: Existing works have identified the limitation of top-$1$ attack success rate (ASR) as a metric to evaluate the attack strength but exclusively investigated it in the white-box setting, while our work extends it to a more practical black-box setting: transferable attack. It is widely reported that stronger I-FGSM transfers worse than simple FGSM, leading to a popular belief that transferability is at odds with the white-box attack strength. Our work challenges this belief with empirical finding that stronger attack actually transfers better for the general top-$k$ ASR indicated by the interest class rank (ICR) after attack. For increasing the attack strength, with an intuitive interpretation of the logit gradient from the geometric perspective, we identify that the weakness of the commonly used losses lie in prioritizing the speed to fool the network instead of maximizing its strength. To this end, we propose a new normalized CE loss that guides the logit to be updated in the direction of implicitly maximizing its rank distance from the ground-truth class. Extensive results in various settings have verified that our proposed new loss is simple yet effective for top-$k$ attack. Code is available at: \url{https://bit.ly/3uCiomP}


title: Gravitational Waves in Einstein-Cartan Theory: On the Effects of Dark Matter Spin Tensor

abstract: This article studies the effects of an arbitrary dark matter spin tensor on the propagation of gravitational wave amplitude in the context of Einstein-Cartan theory. We choose to work with an arbitrary spin tensor because, given our ignorance of the nature of dark matter, it is sensible not to make further hypotheses on its spin and not to assume any particular model for its spin tensor (or its vanishing). The analysis focuses on a \textquotedblleft weak-torsion regime,\textquotedblright\ such that gravitational wave emission, at leading and subleading orders, does not deviate from standard General Relativity. We show that, in principle, the background torsion induced by an eventual dark matter spin component could lead to an anomalous dampening or amplification of the gravitational wave amplitude, after going across a long cosmological distance. We assess the importance of this torsion-induced anomalous amplitude propagation for binary black hole mergers in a way as model-free as possible in Einstein-Cartan gravity. It is possible to prove that at its best, for realistic late-universe cosmological scenarios, the effect is tiny and falls below detection thresholds, even for near-future interferometers such as LISA. Therefore, detecting this effect may not be impossible, but it is still beyond our technological capabilities. As a model-independent result in the Einstein-Cartan context, it also implies that mergers are robust standard sirens without considering any potential dark-matter-induced torsional effects.

categories: gr-qc astro-ph.CO hep-th math-ph math.MP


title: Atomic-scale origin of the low grain-boundary resistance in perovskite solid electrolytes

abstract: Oxide solid electrolytes (OSEs) have the potential to achieve improved safety and energy density for lithium-ion batteries, but their high grain-boundary (GB) resistance is a general bottleneck. In the most well studied perovskite OSE, Li3xLa2/3-xTiO3 (LLTO), the ionic conductivity of GBs is about three orders of magnitude lower than that of the bulk. In contrast, the related Li0.375Sr0.4375Ta0.75Zr0.25O3 (LSTZ0.75) perovskite exhibits low GB resistance for reasons yet unknown. Here, we used aberration-corrected scanning transmission electron microscopy and spectroscopy, along with an active learning moment tensor potential, to reveal the atomic scale structure and composition of LSTZ0.75 GBs. Vibrational electron energy loss spectroscopy is applied for the first time to characterize the otherwise unmeasurable Li distribution in GBs of LSTZ0.75. We found that Li depletion, which is a major reason for the low GB ionic conductivity of LLTO, is absent for the GBs of LSTZ0.75. Instead, the low GB resistivity of LSTZ0.75 is attributed to the formation of a unique defective cubic perovskite interfacial structure that contained abundant vacancies. Our study provides insights into the atomic scale mechanisms of low GB resistivity and sheds light on possible paths for designing OSEs with high total ionic conductivity.


title: Machine learning for a finite size correction in periodic coupled cluster theory calculations

abstract: We introduce a straightforward Gaussian process regression (GPR) model for the transition structure factor of metal periodic coupled cluster singles and doubles (CCSD) calculations. This is inspired by the method introduced by Liao and Gr\"uneis for interpolating over the transition structure factor to obtain a finite size correction for CCSD [J. Chem. Phys. 145, 141102 (2016)], and by our own prior work using the transition structure factor to efficiently converge CCSD for metals to the thermodynamic limit [Nat. Comput. Sci. 1, 801 (2021)]. In our CCSD-FS-GPR method to correct for finite size errors, we fit the structure factor to a 1D function in the momentum transfer, $G$. We then integrate over this function by projecting it onto a k-point mesh to obtain comparisons with extrapolated results. Results are shown for lithium, sodium, and the uniform electron gas.


title: The Spectrum of Triangle-free Graphs

abstract: Denote by $q_n(G)$ the smallest eigenvalue of the signless Laplacian matrix of an $n$-vertex graph $G$. Brandt conjectured in 1997 that for regular triangle-free graphs $q_n(G) \leq \frac{4n}{25}$. We prove a stronger result: If $G$ is a triangle-free graph then $q_n(G) \leq \frac{15n}{94}< \frac{4n}{25}$. Brandt's conjecture is a subproblem of two famous conjectures of Erd\H{o}s: (1) Sparse-Half-Conjecture: Every $n$-vertex triangle-free graph has a subset of vertices of size $\lceil\frac{n}{2}\rceil$ spanning at most $n^2/50$ edges. (2) Every $n$-vertex triangle-free graph can be made bipartite by removing at most $n^2/25$ edges. In our proof we use linear algebraic methods to upper bound $q_n(G)$ by the ratio between the number of induced paths with 3 and 4 vertices. We give an upper bound on this ratio via the method of flag algebras.


title: Perceptive, non-linear Speech Processing and Spiking Neural Networks

abstract: Source separation and speech recognition are very difficult in the context of noisy and corrupted speech. Most conventional techniques need huge databases to estimate speech (or noise) density probabilities to perform separation or recognition. We discuss the potential of perceptive speech analysis and processing in combination with biologically plausible neural network processors. We illustrate the potential of such non-linear processing of speech on a source separation system inspired by an Auditory Scene Analysis paradigm. We also discuss a potential application in speech recognition.


title: Tooth Instance Segmentation on Panoramic Dental Radiographs Using U-Nets and Morphological Processing

abstract: Automatic teeth segmentation in panoramic x-ray images is an important research subject of the image analysis in dentistry. In this study, we propose a post-processing stage to obtain a segmentation map in which the objects in the image are separated, and apply this technique to tooth instance segmentation with U-Net network. The post-processing consists of grayscale morphological and filtering operations, which are applied to the sigmoid output of the network before binarization. A dice overlap score of 95.4 - 0.3% is obtained in overall teeth segmentation. The proposed post-processing stages reduce the mean error of tooth count to 6.15%, whereas the error without post-processing is 26.81%. The performances of both segmentation and tooth counting are the highest in the literature, to our knowledge. Moreover, this is achieved by using a relatively small training dataset, which consists of 105 images. Although the aim in this study is to segment tooth instances, the presented method is applicable to similar problems in other domains, such as separating the cell instances


title: Iterative Reconstruction of the Electron Density and Effective Atomic Number using a Non-Linear Forward Model

abstract: For material identification, characterization, and quantification, it is useful to estimate system-independent material properties that do not depend on the detailed specifications of the X-ray computed tomography (CT) system such as spectral response. System independent rho-e and Z-e (SIRZ) refers to a suite of methods for estimating the system independent material properties of electron density, rho-e, and effective atomic number, Z-e, of an object scanned using dual-energy X-ray CT (DECT). The current state-of-the-art approach, SIRZ-2, makes certain approximations that lead to inaccurate estimates for large atomic numbered (Z-e) materials. In this paper, we present an extension, SIRZ-3, which iteratively reconstructs the unknown rho-e and Z-e while avoiding the limiting approximations made by SIRZ-2. Unlike SIRZ-2, this allows SIRZ-3 to accurately reconstruct rho-e and Z-e even at large values of Z-e. SIRZ-3 relies on the use of a new non-linear differentiable forward measurement model that expresses the DECT measurement data as a direct analytical function of rho-e and Z-e. Leveraging this new forward model, we use an iterative optimization algorithm to reconstruct (or solve for) rho-e and Z-e directly from the DECT data. Compared to SIRZ-2, we show that the magnitude of performance improvement using SIRZ-3 increases with increasing values for Z-e.


title: TransGeo: Transformer Is All You Need for Cross-view Image Geo-localization

abstract: The dominant CNN-based methods for cross-view image geo-localization rely on polar transform and fail to model global correlation. We propose a pure transformer-based approach (TransGeo) to address these limitations from a different perspective. TransGeo takes full advantage of the strengths of transformer related to global information modeling and explicit position information encoding. We further leverage the flexibility of transformer input and propose an attention-guided non-uniform cropping method, so that uninformative image patches are removed with negligible drop on performance to reduce computation cost. The saved computation can be reallocated to increase resolution only for informative patches, resulting in performance improvement with no additional computation cost. This "attend and zoom-in" strategy is highly similar to human behavior when observing images. Remarkably, TransGeo achieves state-of-the-art results on both urban and rural datasets, with significantly less computation cost than CNN-based methods. It does not rely on polar transform and infers faster than CNN-based methods. Code is available at https://github.com/Jeff-Zilence/TransGeo2022.


title: Readout for Calorimetry at Future Colliders: A Snowmass 2021 White Paper

abstract: Calorimeters will provide critical measurements at future collider detectors. As the traditional challenge of high dynamic range, high precision, and high readout rates for signal amplitudes is compounded by increasing granularity and precision timing the readout systems will become increasingly complex. This white paper reviews the challenges and opportunities in calorimeter readout at future collider detectors.


title: Decidability bounds for Presburger arithmetic extended by sine

abstract: We consider Presburger arithmetic extended by the sine function, call this extension sine-Presburger arithmetic ($\sin$-PA), and systematically study decision problems for sets of sentences in $\sin$-PA. In particular, we detail a decision algorithm for existential $\sin$-PA sentences under assumption of Schanuel's conjecture. This procedure reduces decisions to the theory of the ordered additive group of real numbers extended by sine, which is decidable under Schanuel's conjecture. On the other hand, we prove that four alternating quantifier blocks suffice for undecidability of $\sin$-PA sentences. To do so, we explicitly interpret the weak monadic second-order theory of the grid, which is undecidable, in $\sin$-PA.


title: Distributed Stochastic Nash Equilibrium Learning in Locally Coupled Network Games with Unknown Parameters

abstract: In stochastic Nash equilibrium problems (SNEPs), it is natural for players to be uncertain about their complex environments and have multi-dimensional unknown parameters in their models. Among various SNEPs, this paper focuses on locally coupled network games where the objective of each rational player is subject to the aggregate influence of its neighbors. We propose a distributed learning algorithm based on the proximal-point iteration and ordinary least-square estimator, where each player repeatedly updates the local estimates of neighboring decisions, makes its augmented best-response decisions given the current estimated parameters, receives the realized objective values, and learns the unknown parameters. Leveraging the Robbins-Siegmund theorem and the law of large deviations for M-estimators, we establish the almost sure convergence of the proposed algorithm to solutions of SNEPs when the updating step sizes decay at a proper rate.


title: Higher-order magnetohydrodynamic numerics

abstract: In this chapter, we aim at presenting the basic techniques necessary to go beyond the widely accepted paradigm of second-order numerics. We specifically focus on finite-volume schemes for hyperbolic conservation laws occuring in fluid approximations such as the equations of ideal magnetohydrodynamics or the Euler equations of gas dynamics. For the sake of clarity, a simple fourth-order ideal magnetohydrodynamic (MHD) solver which allows to simulate strongly shocked systems serves as an instructive example. Issues that only or mainly arise in the world of higher-order numerics are given specific focus. Alternative algorithms as well as refinements and improvements are dicussed and are referenced to in the literature. As an example of application, some results on decaying compressible turbulence are presented.

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


title: Dynamic Multimodal Fusion

abstract: Deep multimodal learning has achieved great progress in recent years. However, current fusion approaches are static in nature, i.e., they process and fuse multimodal inputs with identical computation, without accounting for diverse computational demands of different multimodal data. In this work, we propose dynamic multimodal fusion (DynMM), a new approach that adaptively fuses multimodal data and generates data-dependent forward paths during inference. To this end, we propose a gating function to provide modality-level or fusion-level decisions on-the-fly based on multimodal features and a resource-aware loss function that encourages computational efficiency. Results on various multimodal tasks demonstrate the efficiency and wide applicability of our approach. For instance, DynMM can reduce the computation costs by 46.5% with only a negligible accuracy loss (CMU-MOSEI sentiment analysis) and improve segmentation performance with over 21% savings in computation (NYU Depth V2 semantic segmentation) when compared with static fusion approaches. We believe our approach opens a new direction towards dynamic multimodal network design, with applications to a wide range of multimodal tasks.