title: The eclipse of the V773 Tau B circumbinary disk

abstract: A deep (~70%) and extended (~150 days) eclipse was seen towards the young multiple stellar system V773 Tau in 2010. We interpret it as due to the passage of a circumbinary disk around the B components moving in front of the A components. Our aim is to characterise the orientation and structure of the disk, to refine the orbits of the subcomponents, and to predict when the next eclipse will occur. We combine the photometry from several ground based surveys, construct a model for the light curve of the eclipse, and use high angular resolution imaging to refine the orbits of the three components of the system, A, B and C. Frequency analysis of the light curves, including from the TESS satellite, enables characterisation of the rotational periods of the Aa and Ab stars. A toy model of the circumbinary disk shows that it extends out to approximately 5 au around the B binary and has an inclination of 73 degrees with respect to the orbital plane of AB, where the lower bound of the radius of the disk is constrained by the geometry of the AB orbit and the upper bound is set by the stability of the disk. We identify several frequencies in the photometric data that we attribute to rotational modulation of the Aa and Ab stellar companions. We produce the first determination of the orbit of the more distant C component around the AB system and limit its inclination to 93 degrees. The high inclination and large diameter of the disk, together with the match from theory suggest that B is an almost equal mass, moderately eccentric binary. We identify the rotational periods of the Aa and Ab stars, identify a third frequency in the light curve that we attribute to the orbital period of the stars in the B binary. We predict that the next eclipse will be around 2037, during which both detailed photometric and spectroscopic monitoring will characterise the disk in greater detail.


title: Hypergraph Tur\'{a}n densities can have arbitrarily large algebraic degree

abstract: Grosu [On the algebraic and topological structure of the set of Tur\'{a}n densities. \emph{J. Combin. Theory Ser. B} \textbf{118} (2016) 137--185] asked if there exist an integer $r\ge 3$ and a finite family of $r$-graphs whose Tur\'{a}n density, as a real number, has (algebraic) degree greater than~$r-1$. In this note we show that, for all integers $r\ge 3$ and $d$, there exists a finite family of $r$-graphs whose Tur\'{a}n density has degree at least~$d$, thus answering Grosu's question in a strong form.


title: Learning from Label Relationships in Human Affect

abstract: Human affect and mental state estimation in an automated manner, face a number of difficulties, including learning from labels with poor or no temporal resolution, learning from few datasets with little data (often due to confidentiality constraints) and, (very) long, in-the-wild videos. For these reasons, deep learning methodologies tend to overfit, that is, arrive at latent representations with poor generalisation performance on the final regression task. To overcome this, in this work, we introduce two complementary contributions. First, we introduce a novel relational loss for multilabel regression and ordinal problems that regularises learning and leads to better generalisation. The proposed loss uses label vector inter-relational information to learn better latent representations by aligning batch label distances to the distances in the latent feature space. Second, we utilise a two-stage attention architecture that estimates a target for each clip by using features from the neighbouring clips as temporal context. We evaluate the proposed methodology on both continuous affect and schizophrenia severity estimation problems, as there are methodological and contextual parallels between the two. Experimental results demonstrate that the proposed methodology outperforms all baselines. In the domain of schizophrenia, the proposed methodology outperforms previous state-of-the-art by a large margin, achieving a PCC of up to 78% performance close to that of human experts (85%) and much higher than previous works (uplift of up to 40%). In the case of affect recognition, we outperform previous vision-based methods in terms of CCC on both the OMG and the AMIGOS datasets. Specifically for AMIGOS, we outperform previous SoTA CCC for both arousal and valence by 9% and 13% respectively, and in the OMG dataset we outperform previous vision works by up to 5% for both arousal and valence.


title: The NNLO quark beam function for jet-veto resummation

abstract: We consider the quark beam function that describes collinear initial-state radiation that is constrained by a veto on reconstructed jets. As the veto is imposed on the transverse momenta of the jets, the beam function is subject to rapidity divergences, and we use the collinear-anomaly framework to extract the perturbative matching kernels to next-to-next-to-leading order (NNLO) in the strong-coupling expansion. Our calculation is based on a novel framework that automates the computation of beam functions in Mellin space and it provides the ingredients to extend jet-veto resummations for quark-initiated processes to NNLL$'$ accuracy.


title: Are We Building on the Rock? On the Importance of Data Preprocessing for Code Summarization

abstract: Code summarization, the task of generating useful comments given the code, has long been of interest. Most of the existing code summarization models are trained and validated on widely-used code comment benchmark datasets. However, little is known about the quality of the benchmark datasets built from real-world projects. Are the benchmark datasets as good as expected? To bridge the gap, we conduct a systematic research to assess and improve the quality of four benchmark datasets widely used for code summarization tasks. First, we propose an automated code-comment cleaning tool that can accurately detect noisy data caused by inappropriate data preprocessing operations from existing benchmark datasets. Then, we apply the tool to further assess the data quality of the four benchmark datasets, based on the detected noises. Finally, we conduct comparative experiments to investigate the impact of noisy data on the performance of code summarization models. The results show that these data preprocessing noises widely exist in all four benchmark datasets, and removing these noisy data leads to a significant improvement on the performance of code summarization. We believe that the findings and insights will enable a better understanding of data quality in code summarization tasks, and pave the way for relevant research and practice.


title: Online Video Instance Segmentation via Robust Context Fusion

abstract: Video instance segmentation (VIS) aims at classifying, segmenting and tracking object instances in video sequences. Recent transformer-based neural networks have demonstrated their powerful capability of modeling spatio-temporal correlations for the VIS task. Relying on video- or clip-level input, they suffer from high latency and computational cost. We propose a robust context fusion network to tackle VIS in an online fashion, which predicts instance segmentation frame-by-frame with a few preceding frames. To acquire the precise and temporal-consistent prediction for each frame efficiently, the key idea is to fuse effective and compact context from reference frames into the target frame. Considering the different effects of reference and target frames on the target prediction, we first summarize contextual features through importance-aware compression. A transformer encoder is adopted to fuse the compressed context. Then, we leverage an order-preserving instance embedding to convey the identity-aware information and correspond the identities to predicted instance masks. We demonstrate that our robust fusion network achieves the best performance among existing online VIS methods and is even better than previously published clip-level methods on the Youtube-VIS 2019 and 2021 benchmarks. In addition, visual objects often have acoustic signatures that are naturally synchronized with them in audio-bearing video recordings. By leveraging the flexibility of our context fusion network on multi-modal data, we further investigate the influence of audios on the video-dense prediction task, which has never been discussed in existing works. We build up an Audio-Visual Instance Segmentation dataset, and demonstrate that acoustic signals in the wild scenarios could benefit the VIS task.


title: Breaking of internal waves parametrically excited by ageostrophic anticyclonic instability

abstract: A gradient-wind balanced flow with an elliptic streamline parametrically excites internal inertia-gravity waves through ageostrophic anticyclonic instability (AAI). This study numerically investigates the breaking of internal waves and the following turbulence generation resulting from the AAI. In our simulation, we periodically distort the calculation domain following the streamlines of an elliptic vortex and integrate the equations of motion using a Fourier spectral method. This technique enables us to exclude the overall structure of the large-scale vortex from the computation and concentrate on resolving the small-scale waves and turbulence. From a series of experiments, we identify two different scenarios of wave breaking conditioned on the magnitude of the instability growth rate scaled by the buoyancy frequency, $\lambda/N$. First, when $\lambda/N\gtrsim0.008$, the primary wave amplitude excited by AAI quickly goes far beyond the overturning threshold and directly breaks. The resulting state is thus strongly nonlinear turbulence. Second, if $\lambda/N\lesssim0.008$, weak wave-wave interactions begin to redistribute energy across frequency space before the primary wave reaches a breaking limit. Then, after a sufficiently long time, the system approaches a Garrett-Munk-like stationary spectrum, in which wave breaking occurs at finer vertical scales. Throughout the experimental conditions, the growth and decay time scales of the primary wave energy are well correlated. However, since the primary wave amplitude reaches a prescribed limit in one scenario but not in the other, the energy dissipation rates exhibit two types of scaling properties. This scaling classification has similarities and differences with D'Asaro and Lien's (2000) wave-turbulence transition model.


title: Investigating the Impact of Independent Rule Fitnesses in a Learning Classifier System

abstract: Achieving at least some level of explainability requires complex analyses for many machine learning systems, such as common black-box models. We recently proposed a new rule-based learning system, SupRB, to construct compact, interpretable and transparent models by utilizing separate optimizers for the model selection tasks concerning rule discovery and rule set composition.This allows users to specifically tailor their model structure to fulfil use-case specific explainability requirements. From an optimization perspective, this allows us to define clearer goals and we find that -- in contrast to many state of the art systems -- this allows us to keep rule fitnesses independent. In this paper we investigate this system's performance thoroughly on a set of regression problems and compare it against XCSF, a prominent rule-based learning system. We find the overall results of SupRB's evaluation comparable to XCSF's while allowing easier control of model structure and showing a substantially smaller sensitivity to random seeds and data splits. This increased control can aid in subsequently providing explanations for both training and final structure of the model.


title: A galaxy-driven model of type Ia supernova luminosity variations

abstract: Type Ia supernovae (SNe Ia) are used as standardisable candles to measure cosmological distances, but differences remain in their corrected luminosities which display a magnitude step as a function of host galaxy properties such as stellar mass and rest-frame $U-R$ colour. Identifying the cause of these steps is key to cosmological analyses and provides insight into SN physics. Here we investigate the effects of SN progenitor ages on their light curve properties using a galaxy-based forward model that we compare to the Dark Energy Survey 5-year SN Ia sample. We trace SN Ia progenitors through time and draw their light-curve width parameters from a bimodal distribution according to their age. We find that an intrinsic luminosity difference between SNe of different ages cannot explain the observed trend between step size and SN colour. The data split by stellar mass are better reproduced by following recent work implementing a step in total-to-selective dust extinction ratio $(R_V)$ between low- and high-mass hosts, although an additional intrinsic luminosity step is still required to explain the data split by host galaxy $U-R$. Modelling the $R_V$ step as a function of galaxy age provides a better match overall. Additional age vs. luminosity steps marginally improve the match to the data, although most of the step is absorbed by the width vs. luminosity coefficient $\alpha$. Furthermore, we find no evidence that $\alpha$ varies with SN age.


title: Multi-Behavior Hypergraph-Enhanced Transformer for Sequential Recommendation

abstract: Learning dynamic user preference has become an increasingly important component for many online platforms (e.g., video-sharing sites, e-commerce systems) to make sequential recommendations. Previous works have made many efforts to model item-item transitions over user interaction sequences, based on various architectures, e.g., recurrent neural networks and self-attention mechanism. Recently emerged graph neural networks also serve as useful backbone models to capture item dependencies in sequential recommendation scenarios. Despite their effectiveness, existing methods have far focused on item sequence representation with singular type of interactions, and thus are limited to capture dynamic heterogeneous relational structures between users and items (e.g., page view, add-to-favorite, purchase). To tackle this challenge, we design a Multi-Behavior Hypergraph-enhanced Transformer framework (MBHT) to capture both short-term and long-term cross-type behavior dependencies. Specifically, a multi-scale Transformer is equipped with low-rank self-attention to jointly encode behavior-aware sequential patterns from fine-grained and coarse-grained levels. Additionally, we incorporate the global multi-behavior dependency into the hypergraph neural architecture to capture the hierarchical long-range item correlations in a customized manner. Experimental results demonstrate the superiority of our MBHT over various state-of-the-art recommendation solutions across different settings. Further ablation studies validate the effectiveness of our model design and benefits of the new MBHT framework. Our implementation code is released at: https://github.com/yuh-yang/MBHT-KDD22.


title: Two results on $x^r + y^r = dz^p$

abstract: This note proves two theorems regarding Fermat-type equation $x^r + y^r = dz^p$ where $r \geq 5$ is a prime. Our main result shows that, for infinitely many integers~$d$, the previous equation has no non-trivial primitive solutions such that $2 \mid x+y$ or $r \mid x+y$, for a set of exponents $p$ of positive density. We use the modular method with a symplectic argument to prove this result.


title: Photonic heat transport from weak to strong coupling

abstract: Superconducting circuits provide a favorable platform for quantum thermodynamic experiments. An important component for such experiments is a heat valve, i.e. a device which allows one to control the heat power flowing through the system. Here we theoretically study the heat valve based on a superconducting quantum interference device (SQUID) coupled to two heat baths via two resonators. The heat current in such system can be tuned by magnetic flux. We investigate how does the heat current modulation depend on the coupling strength g between the SQUID and the resonators. In the weak coupling regime the heat current modulation grows as g2, but, surprisingly, at the intermediate coupling it can be strongly suppressed. This effect is linked to the resonant nature of the heat transport at weak coupling, where the heat current dependence on the magnetic flux is a periodic set of narrow peaks. At the intermediate coupling, the peaks become broader and overlap, thus reducing the heat modulation. At very strong coupling the heat modulation grows again and finally saturates at a constant value.


title: Rigidity results on Liouville equation

abstract: We give a complete classification of solutions bounded from above of the Liouville equation $$-\Delta u=e^{2u}\quad\mbox{in}\quad {\mathbf{R}}^2.$$ More generally, solutions in the class $$N:=\{ u:\limsup_{z\to\infty} u(z)/\log|z|:=k(u)<\infty\}$$ are described. As a consequence, we obtain five rigidity results. First, $k(u)$ can take only a discrete set of values: either $k=-2$, or $2k$ is a non-negative integer. Second, $u\to-\infty$ as $z\to\infty$, if and only if $u$ is radial about some point. Third, if $u$ is symmetric with respect to $x$ and $y$ axes and $u_x<0,\; u_y<0$ in the first quadrant then $u$ is radially symmetric. Fourth, if $u$ is concave and bounded from above, then $u$ is one-dimensional. Fifth, if $u$ is bounded from above, and the diameter of ${\mathbf{R}}^2$ with the metric $e^{2u}\delta$ is $\pi$, where $\delta$ is the Euclidean metric, then $u$ is either radial about a point or one-dimensional. In addition, we extend the concavity rigidity result on Liouville equation in higher dimensions.


title: Ego-motion Estimation Based on Fusion of Images and Events

abstract: Event camera is a novel bio-inspired vision sensor that outputs event stream. In this paper, we propose a novel data fusion algorithm called EAS to fuse conventional intensity images with the event stream. The fusion result is applied to some ego-motion estimation frameworks, and is evaluated on a public dataset acquired in dim scenes. In our 3-DoF rotation estimation framework, EAS achieves the highest estimation accuracy among intensity images and representations of events including event slice, TS and SITS. Compared with original images, EAS reduces the average APE by 69%, benefiting from the inclusion of more features for tracking. The result shows that our algorithm effectively leverages the high dynamic range of event cameras to improve the performance of the ego-motion estimation framework based on optical flow tracking in difficult illumination conditions.


title: MultiShape: A Spectral Element Method, with Applications to Dynamic Density Functional Theory and PDE-Constrained Optimization

abstract: A numerical framework is developed to solve various types of PDEs on complicated domains, including steady and time-dependent, non-linear and non-local PDEs, with different boundary conditions that can also include non-linear and non-local terms. This numerical framework, called MultiShape, is a class in Matlab, and the software is open source. We demonstrate that MultiShape is compatible with other numerical methods, such as differential--algebraic equation solvers and optimization algorithms. The numerical implementation is designed to be user-friendly, with most of the set-up and computations done automatically by MultiShape and with intuitive operator definition, notation, and user-interface. Validation tests are presented, before we introduce three examples motivated by applications in Dynamic Density Functional Theory and PDE-constrained optimization, illustrating the versatility of the method.


title: Integer programming models for the semi-obnoxious p-median problem

abstract: The p-median problem concerns the location of facilities so that the sum of distances between the demand points and their nearest facility is minimized. We study a variant of this classic location problem where minimum distance constraints exist both between the facilities and between the facilities and the demand points. This specific type of problem can be used to model situations where the facilities to be located are semi-obnoxious. But despite its relevance to real life scenarios, it has received little attention within the vast literature on location problems. We present twelve ILP models for this problem, coupling three formulations of the p-median problem with four formulations of the distance constraints. We utilize Gurobi Optimizer v9.0.3 in order to compare these ILP models on a large dataset of problems. Experimental results demonstrate that the classic p-median model proposed by ReVelle \& Swain and the model proposed by Rosing et al. are the best performers.


title: Self-calibration and improving image fidelity for ALMA and other radio interferometers

abstract: This manual is intended to help ALMA and other interferometer users improve images by recognising limitations and how to overcome them and deciding when and how to use self-calibration. The images provided by the ALMA Science Archive are calibrated using standard observing and data processing routines, including a quality assurance process to make sure that the observations meet the proposer's science requirements. This may not represent the full potential of the data, since any interferometry observation can be imaged with a range of resolutions and surface brightness sensitivity. The separation between phase calibration source and target usually limits the target dynamic range to a few hundred (or 50--100 for challenging conditions) but if the noise in the target field has not reached the thermal limit, improvements may be possible using self-calibration. This often requires judgements based on the target properties and is not yet automated for all situations. This manual provides background on the instrumental and atmospheric causes of visibility phase and amplitude errors, their effects on imaging and how to improve the signal to noise ratio and image fidelity by self-calibration. We introduce the conditions for self-calibration to be useful and how to estimate calibration parameter values for a range of observing modes (continuum, spectral line etc.). We also summarise more general error recognition and other techniques to tackle imaging problems. The examples are drawn from ALMA interferometric data processed using CASA, but the principles are generally applicable to most similar cm to sub-mm imaging.


title: Squarefree values of polynomial discriminants II

abstract: We determine the density of integral binary forms of given degree that have squarefree discriminant, proving for the first time that the lower density is positive. Furthermore, we determine the density of integral binary forms that cut out maximal orders in number fields. The latter proves, in particular, an ``arithmetic Bertini theorem'' conjectured by Poonen for $\mathbb{P}^1_\mathbb{Z}$. Our methods also allow us to prove that there are $\gg X^{1/2+1/(n-1)}$ number fields of degree~$n$ having associated Galois group~$S_n$ and absolute discriminant less than $X$, improving the best previously known lower bound of $\gg X^{1/2+1/n}$. Finally, our methods correct an error in and thus resurrect earlier (retracted) results of Nakagawa on lower bounds for the number of totally unramified $A_n$-extensions of quadratic number fields of bounded discriminant.


title: Metal Borohydrides as high-$T_{c}$ ambient pressure superconductors

abstract: The extreme pressures required to stabilize the recently discovered superhydrides represent a major obstacle to their practical application. In this paper, we propose a novel route to attain high-temperature superconductivity in hydrides at ambient pressure, by doping commercial metal borohydrides. Using first-principles calculations based on Density Functional Theory and Migdal-Eliashberg theory, we demonstrate that in Ca(BH$_4$)$_2$ a moderate hole doping of 0.03 holes per formula unit, obtained through a partial replacement of Ca with monovalent K, is sufficient to achieve $T_c$'s as high as 110 K. The high-$T_c$ arises because of the strong electron-phonon coupling between the B-H $\sigma$ molecular orbitals and bond-stretching phonons. Using a random sampling of large supercells to estimate the local effects of doping, we show that the required doping can be achieved without significant disruption of the electronic structure and at moderate energetic cost. Given the wide commercial availability of metal borohydrides, the ideas presented here can find prompt experimental confirmation. If successful, the synthesis of high-$T_c$ doped borohydrides will represent a formidable advancement towards technological exploitation of conventional superconductors.


title: Bound orbit domains in the phase space of the Kerr geometry

abstract: We derive the conditions for a non-equatorial eccentric bound orbit to exist around a Kerr black hole in two-parameter spaces: the energy, angular momentum of the test particle, spin of the black hole, and Carter's constant space ($E$, $L$, $a$, $Q$), and eccentricity, inverse-latus rectum space ($e$, $\mu$, $a$, $Q$). These conditions distribute various kinds of bound orbits in different regions of the ($E$, $L$) and ($e$, $\mu$) planes, depending on which pair of roots of the effective potential forms a bound orbit. We provide a prescription to select these parameters for bound orbits, which are useful inputs to study bound trajectory evolution in various astrophysical applications like simulations of gravitational wave emission from extreme-mass ratio inspirals, relativistic precession around black holes, and the study of gyroscope precession as a test of general relativity.


title: Energy Dependence of the Breit-Wheeler process in Heavy-Ion Collisions and its Application to Nuclear Charge Radius Measurements

abstract: The collision energy dependence of the cross section and the transverse momentum distribution of dielectrons from the Breit-Wheeler process in heavy-ion collisions are computed in the lowest-order QED and found to be sensitive to the nuclear charge distribution and the infrared-divergence of the ultra-Lorentz boosted Coulomb field. Within a given experimental kinematic acceptance, the cross section is found to increase while the pair transverse momentum ($\sqrt{\langle p_{T}^{2} \rangle}$) decreases with increasing beam energy. We demonstrate that the transverse-momentum component of Weizsacker-Williams photons is due to the finite extent of the charge source and electric field component in the longitudinal direction. We further clarify the connection between the nuclear charge distribution and the kinematics of produced $e^+e^-$ from the Breit-Wheeler process, and propose a criterion for the validity of the Breit-Wheeler process in relativistic heavy-ion collisions. Following this approach we demonstrate that the experimental measurements of the Breit-Wheeler process in ultra-relativistic heavy-ion collisions can be used to quantitatively constrain the nuclear charge radius. The extracted parameters show sensitivity to the impact parameter dependence, and can be used to study the initial-state and final-state effects in hadronic interactions.


title: Quantum modulation of a coherent state wavepacket with a single electron spin

abstract: The interaction of quantum objects lies at the heart of fundamental quantum physics and is key to a wide range of quantum information technologies. Photon-quantum-emitter interactions are among the most widely studied. Two-qubit interactions are generally simplified into two quantum objects in static well-defined states . In this work we explore a fundamentally new dynamic type of spin-photon interaction. We demonstrate modulation of a coherent narrowband wavepacket with another truly quantum object, a quantum dot with ground state spin degree of freedom. What results is a quantum modulation of the wavepacket phase (either 0 or {\pi} but no values in between), a new quantum state of light that cannot be described classically.


title: On globally solving nonconvex trust region subproblem via projected gradient method

abstract: The trust region subproblem (TRS) is to minimize a possibly nonconvex quadratic function over a Euclidean ball. There are typically two cases for (TRS), the so-called ``easy case'' and ``hard case''. Even in the ``easy case'', the sequence generated by the classical projected gradient method (PG) may converge to a saddle point at a sublinear local rate, when the initial point is arbitrarily selected from a nonzero measure feasible set. To our surprise, when applying (PG) to solve a cheap and possibly nonconvex reformulation of (TRS), the generated sequence initialized with {\it any} feasible point almost always converges to its global minimizer. The local convergence rate is at least linear for the ``easy case'', without assuming that we have possessed the information that the ``easy case'' holds. We also consider how to use (PG) to globally solve equality-constrained (TRS).


title: Analysis of $\Omega(2012)$ as a molecule in the chiral quark model

abstract: Inspired by the updated information on $\Omega(2012)$ by the Belle Collaboration, we conduct a study of all possible $S$-wave pentaquark systems with quark contents $sssq\bar{q},q=u,d$ in a chiral quark model with the help of Gaussian expansion method. Channel coupling is also considered. The real-scaling method (stabilization method) is employed to identify and check the bound states and the genuine resonances. In addition, the decay widths of all resonances are given. The results show that $\Omega(2012)$ can be interpreted as a $\Xi^*K$ molecular state with quantum number of $IJ^P=0(\frac{3}{2})^-$. Other resonances are obtained: $\Xi^* K^*$ with $IJ^P=0(\frac{1}{2})^-$ and$0(\frac{3}{2})^-$, $\Omega\pi$ with $IJ^P=1(\frac{3}{2})^-$. These pentaquark states is expected to be further verified in future experiments.


title: Sprague-Grundy values and complexity for LCTR

abstract: Given an integer partition of $n$, we consider the impartial combinatorial game LCTR in which moves consist of removing either the left column or top row of its Young diagram. We show that for both normal and mis\`ere play, the optimal strategy can consist mostly of mirroring the opponent's moves. We also establish that both LCTR and Downright are domestic as well as returnable, and on the other hand neither tame nor forced. For both games, those structural observations allow for computing the Sprague-Grundy value any position in $O(\log(n))$ time, assuming that the time unit allows for reading an integer, or performing a basic arithmetic operation. This improves on the previously known bound of $O(n)$ due to Ili\'c (2019). We also cover some other complexity measures of both games, such as state-space complexity, and number of leaves and nodes in the corresponding game tree.


title: A stochastic model to reproduce the star-formation history of individual galaxies in hydrodynamic simulations

abstract: The star formation history (SFH) of galaxies is critical for understanding galaxy evolution. Hydrodynamical simulations enable us to precisely reconstruct the SFH of galaxies and establish a link to the underlying physical processes. In this work, we present a model to describe individual galaxies' SFHs from three simulations: TheThreeHundred, Illustris-1 and TNG100-1. This model divides the galaxy SFH into two distinct components: the "main sequence" and the "variation". The "main sequence" part is generated by tracing the history of the $SFR-M_*$ main sequence of galaxies across time. The "variation" part consists of the scatter around the main sequence, which is reproduced by fractional Brownian motions. We find that: 1) The evolution of the main sequence varies between simulations; 2) fractional Brownian motions can reproduce many features of SFHs, however, discrepancies still exist; 3) The variations and mass-loss rate are crucial for reconstructing the SFHs of the simulations. This model provides a fair description of the SFHs in simulations. On the other hand, by correlating the fractional Brownian motion model to simulation data, we provide a 'standard' against which to compare simulations.


title: Square Root Compression and Noise Effects in Digitally Transformed Images

abstract: We report on a particular example of noise and data representation interacting to introduce systematic error. Many instruments collect integer digitized values and appy nonlinear coding, in particular square-root coding, to compress the data for transfer or downlink; this can introduce surprising systematic errors when they are decoded for analysis. Square root coding and subsequent decoding typically introduces a variable, $\pm 1$ count value-dependent systematic bias in the data after reconstitution. This is significant when large numbers of measurements (e.g., image pixels) are averaged together. Using direct modeling of the probabiliity distribution of particular coded values in the presence of instrument noise, one may apply Bayes' Theorem to construct a decoding table that reduces this error source to a very small fraction of a digitizer step; in our example, systematic error from square root coding is reduced by a factor of 20 from 0.23 count RMS to 0.013 count RMS. The method is suitable both for new experiments such as the upcoming PUNCH mission, and also for post facto application to existing data sets -- even if the instrument noise properties are only loosely known. Further, the method does not depend on the specifics of the coding formula, and may be applied to other forms of nonlinear coding or representation of data values.


title: Nanosecond machine learning regression with deep boosted decision trees in FPGA for high energy physics

abstract: We present a novel application of the machine learning / artificial intelligence method called boosted decision trees to estimate physical quantities on field programmable gate arrays (FPGA). The software package fwXmachina features a new architecture called parallel decision paths that allows for deep decision trees with arbitrary number of input variables. It also features a new optimization scheme to use different numbers of bits for each input variable, which produces optimal physics results and ultraefficient FPGA resource utilization. Problems in high energy physics of proton collisions at the Large Hadron Collider (LHC) are considered. Estimation of missing transverse momentum (ETmiss) at the first level trigger system at the High Luminosity LHC (HL-LHC) experiments, with a simplified detector modeled by Delphes, is used to benchmark and characterize the firmware performance. The firmware implementation with a maximum depth of up to 10 using eight input variables of 16-bit precision gives a latency value of O(10) ns, independent of the clock speed, and O(0.1)% of the available FPGA resources without using digital signal processors.


title: Identification of Substation Configurations in Modern Power Systems using Artificial Intelligence

abstract: Power system transmission network topology is utilized in energy management system applications. Substation configurations are fundamental to transmission network topology processing. Modern power systems consisting of renewable energy sources require reliable and fast network topology processing due to the variable nature of wind and solar power plants. Currently used transmission network topology processing, which is based on the relay signals communicated through SCADA is not highly reliable or highly accurate. Substation configuration identification (SCI) for different substation arrangements including main and transfer bus arrangement (MTBA), ring bus arrangement (RBA), and single bus arrangement (SBA) is investigated. Synchrophasor measurement based SCI for functional arrangements (FA) using artificial intelligence (AI) approaches is proposed in this paper. This method improves monitoring FA. Typical results for MTBA, RBA and SBA substation configuration identification is presented. A modified two-area four-machine power system model with two grid connected solar PV plants consisting of MTBA, RBA and SBA is simulated on real-time digital simulator. AI based SCI is shown to accurately identify all possible FAs for the three substation arrangements under any operating condition.


title: Time-Optimal Trajectory Planning with Interaction with the Environment

abstract: Optimal motion planning along prescribed paths can be solved with several techniques, but most of them do not take into account the wrenches exerted by the end-effector when in contact with the environment. When a dynamic model of the environment is not available, no consolidated methodology exists to consider the effect of the interaction. Regardless of the specific performance index to optimize, this article proposes a strategy to include external wrenches in the optimal planning algorithm, considering the task specifications. This procedure is instantiated for minimum-time trajectories and validated on a real robot performing an interaction task under admittance control. The results prove that the inclusion of end-effector wrenches affect the planned trajectory, in fact modifying the manipulator's dynamic capability.


title: Towards Real-time High-Definition Image Snow Removal: Efficient Pyramid Network with Asymmetrical Encoder-decoder Architecture

abstract: In winter scenes, the degradation of images taken under snow can be pretty complex, where the spatial distribution of snowy degradation is varied from image to image. Recent methods adopt deep neural networks to directly recover clean scenes from snowy images. However, due to the paradox caused by the variation of complex snowy degradation, achieving reliable High-Definition image desnowing performance in real time is a considerable challenge. We develop a novel Efficient Pyramid Network with asymmetrical encoder-decoder architecture for real-time HD image desnowing. The general idea of our proposed network is to utilize the multi-scale feature flow fully and implicitly mine clean cues from features. Compared with previous state-of-the-art desnowing methods, our approach achieves a better complexity-performance trade-off and effectively handles the processing difficulties of HD and Ultra-HD images. The extensive experiments on three large-scale image desnowing datasets demonstrate that our method surpasses all state-of-the-art approaches by a large margin both quantitatively and qualitatively, boosting the PSNR metric from 31.76 dB to 34.10 dB on the CSD test dataset and from 28.29 dB to 30.87 dB on the SRRS test dataset.


title: Incommensurate and multiple-$\boldsymbol{q}$ magnetic misfit order in the frustrated quantum-spin-ladder material antlerite, Cu$_3$SO$_4$(OH)$_4$

abstract: In frustrated magnetic systems, the competition amongst interactions can introduce extremely high degeneracy and prevent the system from readily selecting a unique ground state. In such cases, the magnetic order is often exquisitely sensitive to the balance among the interactions, allowing tuning among novel magnetically ordered phases. In antlerite, Cu$_3$SO$_4$(OH)$_4$, Cu$^{2+}$ ($S=1/2$) quantum spins populate three-leg zigzag ladders in a highly frustrated quasi-one-dimensional structural motif. We demonstrate that at zero applied field, in addition to its recently reported low-temperature phase of coupled ferromagnetic and antiferromagnetic spin chains, this mineral hosts an incommensurate helical+cycloidal state, an idle-spin state, and a multiple-$q$ phase which is the magnetic analog of misfit crystal structures. The antiferromagnetic order on the central leg is reentrant. The high tunability of the magnetism in antlerite makes it a particularly promising platform for pursuing exotic magnetic order.


title: Lower bounds for eigenfunction restrictions in lacunary regions

abstract: Let $(M,g)$ be a compact, smooth Riemannian manifold and $\{u_h\}$ be a sequence of $L^2$-normalized Laplace eigenfunctions that has a localized defect measure $\mu$ in the sense that $ M \setminus \text{supp}(\pi_* \mu) \neq \emptyset$ where $\pi:T^*M \to M$ is the canonical projection. Using Carleman estimates we prove that for any real-smooth closed hypersurface $H \subset (M\setminus \text{supp} (\pi_* \mu))$ sufficiently close to $ \text{supp}(\pi_* \mu),$ and for all $\delta >0,$ $$ \int_{H} |u_h|^2 d\sigma \geq C_{\delta}\, e^{- [\, d(H, \text{supp}(\pi_* \mu)) + \,\delta] /h} $$ as $h \to 0^+$. We also show that the result holds for eigenfunctions of Schr\"odinger operators and give applications to eigenfunctions on warped products and joint eigenfunctions of quantum completely integrable (QCI) systems.


title: Inner Monologue: Embodied Reasoning through Planning with Language Models

abstract: Recent works have shown how the reasoning capabilities of Large Language Models (LLMs) can be applied to domains beyond natural language processing, such as planning and interaction for robots. These embodied problems require an agent to understand many semantic aspects of the world: the repertoire of skills available, how these skills influence the world, and how changes to the world map back to the language. LLMs planning in embodied environments need to consider not just what skills to do, but also how and when to do them - answers that change over time in response to the agent's own choices. In this work, we investigate to what extent LLMs used in such embodied contexts can reason over sources of feedback provided through natural language, without any additional training. We propose that by leveraging environment feedback, LLMs are able to form an inner monologue that allows them to more richly process and plan in robotic control scenarios. We investigate a variety of sources of feedback, such as success detection, scene description, and human interaction. We find that closed-loop language feedback significantly improves high-level instruction completion on three domains, including simulated and real table top rearrangement tasks and long-horizon mobile manipulation tasks in a kitchen environment in the real world.


title: Eccentric equatorial trajectories around a Kerr black hole as a QPO model for M82X-1

abstract: We study the bound orbit conditions for equatorial and eccentric orbits around a Kerr black hole both in the parameter space ($E$, $L$, $a$) representing the energy, angular momentum of the test particle, and spin of the black hole, and also ($e$, $\mu$, $a$) space representing the eccentricity, inverse-latus rectum of the orbit, and spin. We apply these conditions and implement the relativistic precession (RP) model to M82X-1, which is an Intermediate-mass black hole (IMBH) system, where two high-frequency Quasi-Periodic Oscillations (HFQPOs) and a low-frequency QPO were simultaneously observed. Assuming that the QPO frequencies can also be generated by equatorial and eccentric trajectories, we calculate the probability distributions to infer $e$, $a$, and periastron distance, $r_p$, of the orbit giving rise to simultaneous QPOs. We find that an eccentric orbit solution is possible in the region between innermost stable circular orbit (ISCO) and the marginally bound circular orbit (MBCO) for $e=0.2768^{+0.0657}_{-0.0451}$, $a=0.2897 \pm 0.0087$, and $r_p=4.6164^{+0.0694}_{-0.1259}$.


title: Abstraction Logic: A New Foundation for (Computer) Mathematics

abstract: Abstraction logic is a new logic, serving as a foundation of mathematics. It combines features of both predicate logic and higher-order logic: abstraction logic can be viewed both as higher-order logic minus static types as well as predicate logic plus operators and variable binding. We argue that abstraction logic is the best foundational logic possible because it maximises both simplicity and practical expressivity. This argument is supported by the observation that abstraction logic has simpler terms and a simpler notion of proof than all other general logics. At the same time, abstraction logic can formalise both intuitionistic and classical abstraction logic, and is sound and complete for these logics and all other logics extending deduction logic with equality.


title: Intelligent Reflecting Surface Enabled Sensing: Cram\'er-Rao Bound Optimization

abstract: This paper investigates intelligent reflecting surface (IRS) enabled non-line-of-sight (NLoS) wireless sensing, in which an IRS is dedicatedly deployed to assist an access point (AP) to sense a target at its NLoS region. It is assumed that the AP is equipped with multiple antennas and the IRS is equipped with a uniform linear array. We consider two types of target models, namely the point and extended targets, for which the AP aims to estimate the target's direction-of-arrival (DoA) and the target response matrix with respect to the IRS, respectively, based on the echo signals from the AP-IRS-target-IRS-AP link. Under this setup, we jointly design the transmit beamforming at the AP and the reflective beamforming at the IRS to minimize the Cram\'er-Rao bound (CRB) on the estimation error. Towards this end, we first obtain the CRB expressions for the two target models in closed form. It is shown that in the point target case, the CRB for estimating the DoA depends on both the transmit and reflective beamformers; while in the extended target case, the CRB for estimating the target response matrix only depends on the transmit beamformers. Next, for the point target case, we optimize the joint beamforming design to minimize the CRB, via alternating optimization, semi-definite relaxation, and successive convex approximation. For the extended target case, we obtain the optimal transmit beamforming solution to minimize the CRB in closed form. Finally, numerical results show that for both cases, the proposed designs based on CRB minimization achieve improved sensing performance in terms of mean squared error, as compared to other traditional schemes.


title: A density-matrix renormalization group algorithm for simulating quantum circuits with a finite fidelity

abstract: We develop a density-matrix renormalization group (DMRG) algorithm for the simulation of quantum circuits. This algorithm can be seen as the extension of time-dependent DMRG from the usual situation of hermitian Hamiltonian matrices to quantum circuits defined by unitary matrices. For small circuit depths, the technique is exact and equivalent to other matrix product state (MPS) based techniques. For larger depths, it becomes approximate in exchange for an exponential speed up in computational time. Like an actual quantum computer, the quality of the DMRG results is characterized by a finite fidelity. However, unlike a quantum computer, the fidelity depends strongly on the quantum circuit considered. For the most difficult possible circuit for this technique, the so-called "quantum supremacy" benchmark of Google Inc. , we find that the DMRG algorithm can generate bit strings of the same quality as the seminal Google experiment on a single computing core. For a more structured circuit used for combinatorial optimization (Quantum Approximate Optimization Algorithm or QAOA), we find a drastic improvement of the DMRG results with error rates dropping by a factor of 100 compared with random quantum circuits. Our results suggest that the current bottleneck of quantum computers is their fidelities rather than the number of qubits.


title: Making Python Code Idiomatic by Automatic Refactoring Non-Idiomatic Python Code with Pythonic Idioms

abstract: Compared to other programming languages (e.g., Java), Python has more idioms to make Python code concise and efficient. Although pythonic idioms are well accepted in the Python community, Python programmers are often faced with many challenges in using them, for example, being unaware of certain pythonic idioms or do not know how to use them properly. Based on an analysis of 7,638 Python repositories on GitHub, we find that non-idiomatic Python code that can be implemented with pythonic idioms occurs frequently and widely. Unfortunately, there is no tool for automatically refactoring such non-idiomatic code into idiomatic code. In this paper, we design and implement an automatic refactoring tool to make Python code idiomatic. We identify nine pythonic idioms by systematically contrasting the abstract syntax grammar of Python and Java. Then we define the syntactic patterns for detecting non-idiomatic code for each pythonic idiom. Finally, we devise atomic AST-rewriting operations and refactoring steps to refactor non-idiomatic code into idiomatic code. We test and review over 4,115 refactorings applied to 1,065 Python projects from GitHub, and submit 90 pull requests for the 90 randomly sampled refactorings to 84 projects. These evaluations confirm the high-accuracy, practicality and usefulness of our refactoring tool on real-world Python code. Our refactoring tool can be accessed at 47.242.131.128:5000.


title: Max-Min Fairness of Rate-Splitting Multiple Access with Finite Blocklength Communications

abstract: Rate-Splitting Multiple Access (RSMA) has emerged as a flexible and powerful framework for wireless networks. In this paper, we investigate the user fairness of downlink multi-antenna RSMA in short-packet communications with/without cooperative (user-relaying) transmission. We design optimal time allocation and linear precoders that maximize the Max-Min Fairness (MMF) rate with Finite Blocklength (FBL) constraints. The relation between the MMF rate and blocklength of RSMA, as well as the impact of cooperative transmission are investigated for a wide range of network loads. Numerical results demonstrate that RSMA can achieve the same MMF rate as Non-Orthogonal Multiple Access (NOMA) and Space Division Multiple Access (SDMA) with smaller blocklengths (and therefore lower latency), especially in cooperative transmission deployment. Hence, we conclude that RSMA is a promising multiple access for guaranteeing user fairness in low-latency communications.


title: Contrastive Learning for Online Semi-Supervised General Continual Learning

abstract: We study Online Continual Learning with missing labels and propose SemiCon, a new contrastive loss designed for partly labeled data. We demonstrate its efficiency by devising a memory-based method trained on an unlabeled data stream, where every data added to memory is labeled using an oracle. Our approach outperforms existing semi-supervised methods when few labels are available, and obtain similar results to state-of-the-art supervised methods while using only 2.6% of labels on Split-CIFAR10 and 10% of labels on Split-CIFAR100.


title: Set input-to-state stability under input delays for nonlinear systems with disturbances

abstract: The study proposes new results on the set input-to-state stability (ISS) subject to a small input time delay for compact, invariant sets that contains the origin. First, using the nonlinear small-gain theory, we prove a Razumikhin-type theorem that ensures ISS for sets in the setup of functional differential equations (FDEs) with disturbances. Next we show how this theorem can be used to ensure set ISS for nonlinear systems with input delays and disturbances. In comparison to the existing research on set ISS robustness with respect to small time delays at the input, our results are rather broad, retaining the ISS gain and without any constraints on time delayed states. The effectiveness of the method is illustrated through two case-studies addressing set stability for classes of nonlinear oscillators of practical interest.


title: Mechanizing Refinement Types (extended)

abstract: Practical checkers based on refinement types use the combination of implicit semantic sub-typing and parametric polymorphism to simplify the specification and automate the verification of sophisticated properties of programs. However, a formal meta-theoretic accounting of the soundness of refinement type systems using this combination has proved elusive. We present \lambda_RF a core refinement calculus that combines semantic sub-typing and parametric polymorphism. We develop a meta-theory for this calculus and prove soundness of the type system. Finally, we give a full mechanization of our meta-theory using the refinement-type based LiquidHaskell as a proof checker, showing how refinements can be used for mechanization.


title: Htex: Per-Halfedge Texturing for Arbitrary Mesh Topologies

abstract: We introduce per-halfedge texturing (Htex) a GPU-friendly method for texturing arbitrary polygon-meshes without an explicit parameterization. Htex builds upon the insight that halfedges encode an intrinsic triangulation for polygon meshes, where each halfedge spans a unique triangle with direct adjacency information. Rather than storing a separate texture per face of the input mesh as is done by previous parameterization-free texturing methods, Htex stores a square texture for each halfedge and its twin. We show that this simple change from face to halfedge induces two important properties for high performance parameterization-free texturing. First, Htex natively supports arbitrary polygons without requiring dedicated code for, e.g, non-quad faces. Second, Htex leads to a straightforward and efficient GPU implementation that uses only three texture-fetches per halfedge to produce continuous texturing across the entire mesh. We demonstrate the effectiveness of Htex by rendering production assets in real time.


title: GRB 190829A -- A Showcase of Binary Late Evolution

abstract: GRB 190829A is the fourth closest gamma-ray burst (GRB) to date ($z=0.0785$). Owing to its wide range of radio, optical, X-ray, and the very-high-energy (VHE) observations by H.E.S.S., it has become an essential new source examined by various models with complementary approaches. We here show in GRB 190829A the double-prompt pulses and the three-multiwavelength afterglows are consistent with the type II binary-driven hypernova (BdHN II) model. The progenitor is a binary composed of a carbon-oxygen (CO) star and a neutron star (NS) companion. The gravitational collapse of the iron core of the CO star produces a supernova (SN) explosion and leaves behind a new neutron star ($\nu$NS) at its center. The accretion of the SN ejecta onto the NS companion and onto the $\nu$NS via matter fallback spins up the NSs and produces the double-peak prompt emission. The synchrotron emission from the expanding SN ejecta with the energy injection from the rapidly spinning $\nu$NS and its subsequently spin-down leads to the afterglow in the radio, optical, and X-ray. We model the sequence of physical and related radiation processes in BdHNe and focus on individuating the binary properties that play the relevant roles.


title: LudVision -- Remote Detection of Exotic Invasive Aquatic Floral Species using Drone-Mounted Multispectral Data

abstract: Remote sensing is the process of detecting and monitoring the physical characteristics of an area by measuring its reflected and emitted radiation at a distance. It is being broadly used to monitor ecosystems, mainly for their preservation. Ever-growing reports of invasive species have affected the natural balance of ecosystems. Exotic invasive species have a critical impact when introduced into new ecosystems and may lead to the extinction of native species. In this study, we focus on Ludwigia peploides, considered by the European Union as an aquatic invasive species. Its presence can negatively impact the surrounding ecosystem and human activities such as agriculture, fishing, and navigation. Our goal was to develop a method to identify the presence of the species. We used images collected by a drone-mounted multispectral sensor to achieve this, creating our LudVision data set. To identify the targeted species on the collected images, we propose a new method for detecting Ludwigia p. in multispectral images. The method is based on existing state-of-the-art semantic segmentation methods modified to handle multispectral data. The proposed method achieved a producer's accuracy of 79.9% and a user's accuracy of 95.5%.


title: MSP-Former: Multi-Scale Projection Transformer for Single Image Desnowing

abstract: Snow removal causes challenges due to its characteristic of complex degradations. To this end, targeted treatment of multi-scale snow degradations is critical for the network to learn effective snow removal. In order to handle the diverse scenes, we propose a multi-scale projection transformer (MSP-Former), which understands and covers a variety of snow degradation features in a multi-path manner, and integrates comprehensive scene context information for clean reconstruction via self-attention operation. For the local details of various snow degradations, the local capture module is introduced in parallel to assist in the rebuilding of a clean image. Such design achieves the SOTA performance on three desnowing benchmark datasets while costing the low parameters and computational complexity, providing a guarantee of practicality.


title: A Dynamic Programming Framework for Optimal Planning of Redundant Robots Along Prescribed Paths With Kineto-Dynamic Constraints

abstract: Offline optimal planning of trajectories for redundant robots along prescribed task space paths is usually broken down into two consecutive processes: first, the task space path is inverted to obtain a joint space path, then, the latter is parametrized with a time law. If the two processes are separated, they cannot optimize the same objective function, ultimately providing sub-optimal results. In this paper, a unified approach is presented where dynamic programming is the underlying optimization technique. Its flexibility allows accommodating arbitrary constraints and objective functions, thus providing a generic framework for optimal planning of real systems. To demonstrate its applicability to a real world scenario, the framework is instantiated for time-optimality on Franka Emika's Panda robot. The well-known issues associated with the execution of non-smooth trajectories on a real controller are partially addressed at planning level, through the enforcement of constraints, and partially through post-processing of the optimal solution. The experiments show that the proposed framework is able to effectively exploit kinematic redundancy to optimize the performance index defined at planning level and generate feasible trajectories that can be executed on real hardware with satisfactory results.


title: Attacking (and defending) the Maritime Radar System

abstract: Operation of radar equipment is one of the key facilities used by navigators to gather situational awareness about their surroundings. With an ever increasing need for always-running logistics and tighter shipping schedules, operators are relying more and more on computerized instruments and their indications. As a result, modern ships have become a complex cyber-physical system in which sensors and computers constantly communicate and coordinate. In this work, we discuss novel threats related to the radar system, which is one of the most security-sensitive component on a ship. In detail, we first discuss some new attacks capable of compromising the integrity of data displayed on a radar system, with potentially catastrophic impacts on the crew' situational awareness or even safety itself. Then, we present a detection system aimed at highlighting anomalies in the radar video feed, requiring no modifications to the target ship configuration. Finally, we stimulate our detection system by performing the attacks inside of a simulated environment. The experimental results clearly indicate that the attacks are feasible, rather easy to carry out, and hard-to-detect. Moreover, they prove that the proposed detection technique is effective.


title: DWTCP: Ultra Low Latency Congestion Control Protocol for Data Centers

abstract: Congestion control algorithms rely on a variety of congestion signals (packet loss, Explicit Congestion Notification, delay, etc.) to achieve fast convergence, high utilization, and fairness among flows. A key limitation of these congestion signals is that they are either late in feedback or they incur significant overheads. An ideal congestion control must discover any available bandwidth in the network, detect congestion as soon as link utilization approaches full capacity, and react timely to avoid queuing and packet drops, without significant overheads. To this end, this work proposes Scout service that leverages priority queues to infer bandwidth availability and link busyness at the host. The key observation here is that as the high priority queue (HPQ) gets busier, the low priority queue (LPQ) is served less. Therefore, the state of the link can be observed from the LPQ and any congestion can be detected several RTTs earlier than observing the HPQ. We propose a new transport protocol, Double-Window Transmission Control Protocol (DWTCP) that builds upon the Scout service to dynamically adjust its congestion window. Our testbed and simulation-based evaluation demonstrates that Scout enables a data center transport to achieve high throughput, near-zero queues, lower latency, and high fairness.


title: Nonlinear wave interactions in geochemical modeling

abstract: This paper is concerned with the study of the main wave interactions in a system of conservation laws in geochemical modeling. We study the modeling of the chemical complexes on the rock surface. The presence of stable surface complexes affects the relative permeability. We add terms representing surface complexes to the accumulation function in the model presented in \cite{lambert2019nonlinear1}. This addition allows to take into account the interaction of ions with the rock surface in the modeling of the oil recovery by the injection of carbonated water. Compatibility hypotheses with the modeling are made on the coefficients of the system to obtain meaningful solutions. We developed a Riemann solver taking into account the complexity of the interactions and bifurcations of nonlinear waves. Such bifurcations occur at the inflection and resonance surfaces. We present the solution of a generalized eigenvalue problem in a (n+1)-dimensional case, which allows the construction of rarefaction curves. A method to find the discontinuous solutions is also presented. We find the solution path for some examples.


title: TreeExplorer: a coding algorithm for rooted trees with application to wireless and ad hoc routing

abstract: Routing tables in ad hoc and wireless routing protocols can be represented using rooted trees. The constant need for communication and storage of these trees in routing protocols demands an efficient rooted tree coding algorithm. This efficiency is defined in terms of the average code length, and the optimality of the algorithm is measured by comparing the average code length with the entropy of the source. In this work, TreeExplorer is introduced as an easy-to-implement and nearly optimal algorithm for coding rooted tree structures. This method utilizes the number of leaves of the tree as an indicator for choosing the best method of coding. We show how TreeExplorer can improve existing routing protocols for ad hoc and wireless systems, which normally entails a significant communication overhead.


title: On the role of initial coherence in the spin phase-space entropy production rate

abstract: Recent studies have pointed out the intrinsic dependence of figures of merit of thermodynamic relevance -- such as work, heat and entropy production -- on the amount of quantum coherences that is made available to a system. However, whether coherences hinder or enhance the value taken by such quantifiers of thermodynamic performance is yet to be ascertained. We show that, when considering entropy production generated in a process taking a finite-size bipartite quantum system out of equilibrium through local non-unitary channels, no general monotonicity relationship exists between the entropy production and degree of quantum coherence in the state of the system. A direct correspondence between such quantities can be retrieved when considering specific forms of open-system dynamics applied to suitably chosen initial states. Our results call for a systematic study of the role of genuine quantum features in the non-equilibrium thermodynamics of quantum processes.


title: Non-uniqueness theory in sampled STFT phase retrieval

abstract: The reconstruction of a function from its spectrogram (i.e., the absolute value of its short-time Fourier transform (STFT)) arises as a key problem in several important applications, including coherent diffraction imaging and audio processing. It is a classical result that for suitable windows any function can, in principle, be uniquely recovered up to a global phase factor from its spectrogram. However, for most practical applications only discrete samples - typically from a lattice - of the spectrogram are available. This raises the question of whether lattice samples of the spectrogram contain sufficient information for determining a function $f\in L^2(\mathbb{R}^d)$ up to a global phase factor. In the present paper, we answer this question in the negative by providing general non-identifiability results which lead to a non-uniqueness theory for the sampled STFT phase retrieval problem. Precisely, given any dimension $d$, any window function $g$ and any (symplectic or separable) lattice $\mathcal{L} \subseteq \mathbb{R}^d$, we construct pairs of functions $f,h\in L^2(\mathbb{R}^d)$ that do not agree up to a global phase factor, but whose spectrograms agree on $\mathcal{L}$. Our techniques are sufficiently flexible to produce counterexamples to unique recoverability under even more stringent assumptions; for example, if the window function is real-valued, the functions $f,h$ can even be chosen to satisfy $|f|=|h|$. Our results thus reveal the non-existence of a critical sampling density in the absence of phase information, a property which is in stark contrast to uniqueness results in time-frequency analysis.


title: On the Bernstein--Zelevinsky classification in families

abstract: We study the variation of admissible representations of $p$-adic $GL_n$ in families from the point of view of the Bernstein--Zelevinsky classification and show that the ramified parts of these families are rigid. We explain how to apply our results in the context of $GL_n$-eigenvarieties.


title: Totally disconnected semigroup compactifications of topological groups

abstract: We introduce the notion of an introverted Boolean algebra $\cal B$ of closed-and-open subsets of a topological group $G$, show that the associated Stone space $(\nu_{\cal B} G, \nu_{\cal B})$ is a totally disconnected semigroup compactification of $G$, and show that every totally disconnected semigroup compactification of $G$ takes this form. We identify and study the universal totally disconnected semigroup compactification, the universal totally disconnected semitopological semigroup compactification and the universal totally disconnected group compactification of $G$. Our main results are obtained independently of Gelfand theory and well-known properties of the (typically non-totally disconnected) universal compactifications $G^{LUC}$, $G^{WAP}$ and $G^{AP}$, though we do employ Gelfand theory to clarify the relationship between these familiar universal compactifications and their totally disconnected counterparts.


title: DGPO: Discovering Multiple Strategies with Diversity-Guided Policy Optimization

abstract: Most reinforcement learning algorithms seek a single optimal strategy that solves a given task. However, it can often be valuable to learn a diverse set of solutions, for instance, to make an agent's interaction with users more engaging, or improve the robustness of a policy to an unexpected perturbance. We propose Diversity-Guided Policy Optimization (DGPO), an on-policy algorithm that discovers multiple strategies for solving a given task. Unlike prior work, it achieves this with a shared policy network trained over a single run. Specifically, we design an intrinsic reward based on an information-theoretic diversity objective. Our final objective alternately constraints on the diversity of the strategies and on the extrinsic reward. We solve the constrained optimization problem by casting it as a probabilistic inference task and use policy iteration to maximize the derived lower bound. Experimental results show that our method efficiently discovers diverse strategies in a wide variety of reinforcement learning tasks. Compared to baseline methods, DGPO achieves comparable rewards, while discovering more diverse strategies, and often with better sample efficiency.


title: The Science Performance of JWST as Characterized in Commissioning

abstract: This paper characterizes the actual science performance of the James Webb Space Telescope (JWST), as determined from the six month commissioning period. We summarize the performance of the spacecraft, telescope, science instruments, and ground system, with an emphasis on differences from pre-launch expectations. Commissioning has made clear that JWST is fully capable of achieving the discoveries for which it was built. Moreover, almost across the board, the science performance of JWST is better than expected; in most cases, JWST will go deeper faster than expected. The telescope and instrument suite have demonstrated the sensitivity, stability, image quality, and spectral range that are necessary to transform our understanding of the cosmos through observations spanning from near-earth asteroids to the most distant galaxies.


title: Bochner formulas, functional inequalities and generalized Ricci flow

abstract: As a consequence of the Bochner formula for the Bismut connection acting on gradients, we show sharp universal Poincar\'e and log-Sobolev inequalities along solutions to generalized Ricci flow. Using the two-form potential we define a twisted connection on spacetime which determines an adapted Brownian motion on the frame bundle, yielding an adapted Malliavin gradient on path space. We show a Bochner formula for this operator, leading to characterizations of generalized Ricci flow in terms of universal Poincar\'e and log-Sobolev type inequalities for the associated Malliavin gradient and Ornstein-Uhlenbeck operator.


title: GANzzle: Reframing jigsaw puzzle solving as a retrieval task using a generative mental image

abstract: Puzzle solving is a combinatorial challenge due to the difficulty of matching adjacent pieces. Instead, we infer a mental image from all pieces, which a given piece can then be matched against avoiding the combinatorial explosion. Exploiting advancements in Generative Adversarial methods, we learn how to reconstruct the image given a set of unordered pieces, allowing the model to learn a joint embedding space to match an encoding of each piece to the cropped layer of the generator. Therefore we frame the problem as a R@1 retrieval task, and then solve the linear assignment using differentiable Hungarian attention, making the process end-to-end. In doing so our model is puzzle size agnostic, in contrast to prior deep learning methods which are single size. We evaluate on two new large-scale datasets, where our model is on par with deep learning methods, while generalizing to multiple puzzle sizes.


title: Martini 3 Coarse-Grained Force Field for Carbohydrates

abstract: The Martini 3 force field is a full re-parametrization of the Martini coarse-grained model for biomolecular simulations. Due to the improved interaction balance it allows for more accurate description of condensed phase systems. In the present work we develop a consistent strategy to parametrize carbohydrate molecules accurately within the framework of Martini 3. In particular, we develop a canonical mapping scheme that decomposes arbitrarily large carbohydrates into a limited number of fragments. Bead types for these fragments have been assigned by matching physicochemical properties of mono- and disaccharides. In addition, guidelines for assigning bonds, angles, and dihedrals are developed. These guidelines enable a more accurate description of carbohydrate conformations than in the Martini 2 force field. We show that models obtained with this approach are able to accurately reproduce osmotic pressures of carbohydrate water solutions. Furthermore, we provide evidence that the model differentiates correctly the solubility of the poly-glucoses dextran (water soluble) and cellulose (water insoluble, but soluble in ionic-liquids). Finally, we demonstrate that the new building blocks can be applied to glycolipids, being able to reproduce membrane properties and to induce binding of peripheral membrane proteins. These test cases demonstrate the validity and transferability of our approach.


title: Neural Posterior Estimation with Differentiable Simulators

abstract: Simulation-Based Inference (SBI) is a promising Bayesian inference framework that alleviates the need for analytic likelihoods to estimate posterior distributions. Recent advances using neural density estimators in SBI algorithms have demonstrated the ability to achieve high-fidelity posteriors, at the expense of a large number of simulations ; which makes their application potentially very time-consuming when using complex physical simulations. In this work we focus on boosting the sample-efficiency of posterior density estimation using the gradients of the simulator. We present a new method to perform Neural Posterior Estimation (NPE) with a differentiable simulator. We demonstrate how gradient information helps constrain the shape of the posterior and improves sample-efficiency.


title: Anomalous Cyclic in the Neutrino Oscillations

abstract: Neutrino physics is one of the most important topics in particle physics and cosmology. Despite the many physical properties of neutrinos that are understood theoretically and experimentally, it is known that there are many unsolved problems in this area. In this study, we suppose that the deformed space-time caused by the gravitational perturbation can play an important role in neutrino oscillations. We analytically analyzed these effects on the neutrino oscillation and showed that this effect leads to an anomalous cyclic in the neutrino oscillation. The results clearly indicate that these anomalous cyclics depend on the degree of deformation of space-time. The role of this anomaly in neutrino oscillation may be important. There may be a relation between cyclic anomaly and other anomalies of neutrinos such as mass limits, energy limits, oscillation lengths, mixing angles, and speed of neutrinos. The cyclic anomalies might be detected experimentally. If these anomalies are confirmed, it will appear that we need to think more about neutrino physics.


title: Representations by sexternary quadratic forms with coefficients 1,2,5 and 10

abstract: In this article, we find bases for the spaces of modular forms $M_{3}(\Gamma _{0}(40),\left( \frac{d}{\cdot }\right) )$ for $d=-4,-8,-20\text{ and }-40.$ We then derive formulas for the number of representations of a positive integer by all the diagonal sextenary quadratic forms with coefficients $% 1,2,5$ and $10$.


title: Positive co-degree density of hypergraphs

abstract: The \emph{minimum positive co-degree} of a non-empty $r$-graph ${H}$, denoted $\delta_{r-1}^+( {H})$, is the maximum $k$ such that if $S$ is an $(r-1)$-set contained in a hyperedge of $ {H}$, then $S$ is contained in at least $k$ distinct hyperedges of $ {H}$. Given an $r$-graph ${F}$, we introduce the \emph{positive co-degree Tur\'an number} $\mathrm{co^+ex}(n, {F})$ as the maximum positive co-degree $\delta_{r-1}^+(H)$ over all $n$-vertex $r$-graphs $H$ that do not contain $F$ as a subhypergraph. In this paper we concentrate on the behavior of $\mathrm{co^+ex}(n, {F})$ for $3$-graphs $F$. In particular, we determine asymptotics and bounds for several well-known concrete $3$-graphs $F$ (e.g.\ $K_4^-$ and the Fano plane). We also show that, for $r$-graphs, the limit \[ \gamma^+(F) := \lim_{n \rightarrow \infty} \frac{\mathrm{co^+ex}(n, {F})}{n} \] exists, and ``jumps'' from $0$ to $1/r$, i.e., it never takes on values in the interval $(0,1/r)$. Moreover, we characterize which $r$-graphs $F$ have $\gamma^+(F)=0$. Our motivation comes primarily from the study of (ordinary) co-degree Tur\'an numbers where a number of results have been proved that inspire our results.


title: Measurement of Charmonium Production in $p + p$ and $p + d$ Interactions in the Fermilab SeaQuest Experiment

abstract: The Fermilab SeaQuest experiment has measured dimuon events from the interactions of 120GeV proton beam on liquid hydrogen and deuterium targets with dimuon mass between 2 and 9GeV. These dimuon events contain both the Drell-Yan process and the charmonium ($J/\psi$ and $\psi^\prime$) production. Unlike the Drell-Yan process which probes the antiquark distributions in the nucleons, the charmonium production is sensitive to both quark and gluon distributions. SeaQuest has extracted the $\sigma^{pd}/2\sigma^{pp}$ ratio as well as the differential cross sections for charmonium production in the kinematic region of $0.4 < x_F < 0.9$. The $\sigma^{pd}/2\sigma^{pp}$ ratio for charmonium production are found to be significantly different from that of the Drell-Yan process. The measured differential cross sections for charmonium production are compared with theoretical calculations using Color Evaporation Model and Non-Relativistic QCD model.


title: Backdoor Attacks on Crowd Counting

abstract: Crowd counting is a regression task that estimates the number of people in a scene image, which plays a vital role in a range of safety-critical applications, such as video surveillance, traffic monitoring and flow control. In this paper, we investigate the vulnerability of deep learning based crowd counting models to backdoor attacks, a major security threat to deep learning. A backdoor attack implants a backdoor trigger into a target model via data poisoning so as to control the model's predictions at test time. Different from image classification models on which most of existing backdoor attacks have been developed and tested, crowd counting models are regression models that output multi-dimensional density maps, thus requiring different techniques to manipulate. In this paper, we propose two novel Density Manipulation Backdoor Attacks (DMBA$^{-}$ and DMBA$^{+}$) to attack the model to produce arbitrarily large or small density estimations. Experimental results demonstrate the effectiveness of our DMBA attacks on five classic crowd counting models and four types of datasets. We also provide an in-depth analysis of the unique challenges of backdooring crowd counting models and reveal two key elements of effective attacks: 1) full and dense triggers and 2) manipulation of the ground truth counts or density maps. Our work could help evaluate the vulnerability of crowd counting models to potential backdoor attacks.


title: Scalable Bayesian Inference for Detection and Deblending in Astronomical Images

abstract: We present a new probabilistic method for detecting, deblending, and cataloging astronomical sources called the Bayesian Light Source Separator (BLISS). BLISS is based on deep generative models, which embed neural networks within a Bayesian model. For posterior inference, BLISS uses a new form of variational inference known as Forward Amortized Variational Inference. The BLISS inference routine is fast, requiring a single forward pass of the encoder networks on a GPU once the encoder networks are trained. BLISS can perform fully Bayesian inference on megapixel images in seconds, and produces highly accurate catalogs. BLISS is highly extensible, and has the potential to directly answer downstream scientific questions in addition to producing probabilistic catalogs.


title: SafeDrones: Real-Time Reliability Evaluation of UAVs using Executable Digital Dependable Identities

abstract: The use of Unmanned Arial Vehicles (UAVs) offers many advantages across a variety of applications. However, safety assurance is a key barrier to widespread usage, especially given the unpredictable operational and environmental factors experienced by UAVs, which are hard to capture solely at design-time. This paper proposes a new reliability modeling approach called SafeDrones to help address this issue by enabling runtime reliability and risk assessment of UAVs. It is a prototype instantiation of the Executable Digital Dependable Identity (EDDI) concept, which aims to create a model-based solution for real-time, data-driven dependability assurance for multi-robot systems. By providing real-time reliability estimates, SafeDrones allows UAVs to update their missions accordingly in an adaptive manner.


title: Stability of a point charge for the repulsive Vlasov-Poisson system

abstract: We consider solutions of the repulsive Vlasov-Poisson system which are a combination of a point charge and a small gas, i.e.\ measures of the form $\delta_{(\mathcal{X}(t),\mathcal{V}(t))}+\mu^2d{\bf x}d{\bf v}$ for some $(\mathcal{X}, \mathcal{V}):\mathbb{R}\to\mathbb{R}^6$ and a small gas distribution $\mu:\mathbb{R}\to L^2_{{\bf x},{\bf v}}$, and study asymptotic dynamics in the associated initial value problem. If initially suitable moments on $\mu_0=\mu(t=0)$ are small, we obtain a global solution of the above form, and the electric field generated by the gas distribution $\mu$ decays at an almost optimal rate. Assuming in addition boundedness of suitable derivatives of $\mu_0$, the electric field decays at an optimal rate and we derive a modified scattering dynamics for the motion of the point charge and the gas distribution. Our proof makes crucial use of the Hamiltonian structure. The linearized system is transport by the Kepler ODE, which we integrate exactly through an asymptotic action-angle transformation. Thanks to a precise understanding of the associated kinematics, moment and derivative control is achieved via a bootstrap analysis that relies on the decay of the electric field associated to $\mu$. The asymptotic behavior can then be deduced from the properties of Poisson brackets in asymptotic action coordinates.


title: Speed limits on correlations in bipartite quantum systems

abstract: Quantum speed limit is bound on the minimum time a quantum system requires to evolve from an initial state to final state under a given dynamical process. It sheds light on how fast a desired state transformation can take place which is pertinent for design and control of quantum technologies. In this paper, we derive speed limits on correlations such as entanglement, Bell-CHSH correlation, and quantum mutual information of quantum systems evolving under dynamical processes. Our main result is speed limit on an entanglement monotone called negativity which holds for arbitrary dimensional bipartite quantum systems and processes. Another entanglement monotone which we consider is the concurrence. To illustrate efficacy of our speed limits, we analytically and numerically compute the speed limits on the negativity, concurrence, and Bell-CHSH correlation for various quantum processes of practical interest. We are able to show that for practical examples we have considered, some of the speed limits we derived are actually attainable and hence these bounds can be considered to be tight.


title: Resonant Multilevel Amplitude Damping Channels

abstract: We introduce a new set of quantum channels: resonant multilevel amplitude damping (ReMAD) channels. Among other instances, they can describe energy dissipation effects in multilevel atomic systems induced by the interaction with a zero-temperature bosonic environment. At variance with the already known class of multilevel amplitude damping (MAD) channels, this new class of maps allows the presence of an environment unable to discriminate transitions with identical energy gaps. After characterizing the algebra of their composition rules, by analyzing the qutrit case, we show that this new set of channels can exhibit degradability and antidegradability in vast regions of the allowed parameter space. There we compute their quantum capacity and private classical capacity. We show that these capacities can be computed exactly also in regions of the parameter space where the channels aren't degradable nor antidegradable.


title: How Much Entanglement Does a Quantum Code Need?

abstract: In the setting of entanglement-assisted quantum error-correcting codes (EAQECCs), the sender and the receiver have access to pre-shared entanglement. Such codes promise better information rates or improved error handling properties. Entanglement incurs costs and must be judiciously calibrated in designing quantum codes with good performance, relative to their deployment parameters. Revisiting known constructions, we devise tools from classical coding theory to better understand how the amount of entanglement can be varied. We present three new propagation rules and discuss how each of them affects the error handling. Tables listing the parameters of the best performing qubit and qutrit EAQECCs that we can explicitly construct are supplied for reference and comparison.


title: Docent: A content-based recommendation system to discover contemporary art

abstract: Recommendation systems have been widely used in various domains such as music, films, e-shopping etc. After mostly avoiding digitization, the art world has recently reached a technological turning point due to the pandemic, making online sales grow significantly as well as providing quantitative online data about artists and artworks. In this work, we present a content-based recommendation system on contemporary art relying on images of artworks and contextual metadata of artists. We gathered and annotated artworks with advanced and art-specific information to create a completely unique database that was used to train our models. With this information, we built a proximity graph between artworks. Similarly, we used NLP techniques to characterize the practices of the artists and we extracted information from exhibitions and other event history to create a proximity graph between artists. The power of graph analysis enables us to provide an artwork recommendation system based on a combination of visual and contextual information from artworks and artists. After an assessment by a team of art specialists, we get an average final rating of 75% of meaningful artworks when compared to their professional evaluations.


title: Elastic turbulence homogenizes fluid transport in stratified porous media

abstract: Many key environmental, industrial, and energy processes rely on controlling fluid transport within subsurface porous media. These media are typically structurally heterogeneous, often with vertically-layered strata of distinct permeabilities -- leading to uneven partitioning of flow across strata, which can be undesirable. Here, using direct in situ visualization, we demonstrate that polymer additives can homogenize this flow by inducing a purely-elastic flow instability that generates random spatiotemporal fluctuations and excess flow resistance in individual strata. In particular, we find that this instability arises at smaller imposed flow rates in higher-permeability strata, diverting flow towards lower-permeability strata and helping to homogenize the flow. Guided by the experiments, we develop a parallel-resistor model that quantitatively predicts the flow rate at which this homogenization is optimized for a given stratified medium. Thus, our work provides a new approach to homogenizing fluid and passive scalar transport in heterogeneous porous media.

categories: physics.flu-dyn cond-mat.dis-nn cond-mat.soft nlin.CD physics.app-ph


title: A Single-Loop Gradient Descent and Perturbed Ascent Algorithm for Nonconvex Functional Constrained Optimization

abstract: Nonconvex constrained optimization problems can be used to model a number of machine learning problems, such as multi-class Neyman-Pearson classification and constrained Markov decision processes. However, such kinds of problems are challenging because both the objective and constraints are possibly nonconvex, so it is difficult to balance the reduction of the loss value and reduction of constraint violation. Although there are a few methods that solve this class of problems, all of them are double-loop or triple-loop algorithms, and they require oracles to solve some subproblems up to certain accuracy by tuning multiple hyperparameters at each iteration. In this paper, we propose a novel gradient descent and perturbed ascent (GDPA) algorithm to solve a class of smooth nonconvex inequality constrained problems. The GDPA is a primal-dual algorithm, which only exploits the first-order information of both the objective and constraint functions to update the primal and dual variables in an alternating way. The key feature of the proposed algorithm is that it is a single-loop algorithm, where only two step-sizes need to be tuned. We show that under a mild regularity condition GDPA is able to find Karush-Kuhn-Tucker (KKT) points of nonconvex functional constrained problems with convergence rate guarantees. To the best of our knowledge, it is the first single-loop algorithm that can solve the general nonconvex smooth problems with nonconvex inequality constraints. Numerical results also showcase the superiority of GDPA compared with the best-known algorithms (in terms of both stationarity measure and feasibility of the obtained solutions).


title: Crystal and magnetic structures of $R_2$Ni$_2$In compounds ($R$ = Tb and Ho)

abstract: Crystal and magnetic structures of $R_2$Ni$_2$In ($R$ = Tb and Ho) have been studied by powder neutron diffraction at low temperatures. The compounds crystallize in an orthorhombic crystal structure of the Mn$_2$AlB$_2$-type. At low temperatures, the magnetic moments localized solely on the rare earth atoms form antiferromagnetic structures. The Tb magnetic moments, equal to 8.65(6) $\mu_B$ and parallel to the $c$-axis, form a collinear magnetic structure described by the propagation vector $\boldsymbol{k} = [\frac{1}{2}, \frac{1}{2}, \frac{1}{2}]$. This magnetic structure is stable up to the N\'eel temperature equal to 40 K. For Ho$_2$Ni$_2$In a complex, temperature-dependent magnetic structure is detected. In the temperature range 3.5-8.6 K, an incommensurate magnetic structure, described by the propagation vector $\boldsymbol{k}_1 = [0.76, 0, 0.52]$ is observed, while in the temperature interval 2.2-3.1 K the magnetic order is described by two propagation vectors, namely, $\boldsymbol{k}_2 = [\frac{5}{6}, 0.16, \frac{1}{2}]$ and its third harmonics $3\boldsymbol{k}_2 = [\frac{5}{2}, 0.48, \frac{3}{2}]$. Below 2 K, a coexistence of all magnetic structures detected at higher temperatures is observed. For all magnetic phases, the Ho magnetic moments are parallel to the $c$-axis. The low temperature heat capacity data confirm a first order transition near 3 K.


title: Accelerating astronomical and cosmological inference with Preconditioned Monte Carlo

abstract: We introduce Preconditioned Monte Carlo (PMC), a novel Monte Carlo method for Bayesian inference that facilitates efficient sampling of probability distributions with non-trivial geometry. PMC utilises a Normalising Flow (NF) in order to decorrelate the parameters of the distribution and then proceeds by sampling from the preconditioned target distribution using an adaptive Sequential Monte Carlo (SMC) scheme. The results produced by PMC include samples from the posterior distribution and an estimate of the model evidence that can be used for parameter inference and model comparison respectively. The aforementioned framework has been thoroughly tested in a variety of challenging target distributions achieving state-of-the-art sampling performance. In the cases of primordial feature analysis and gravitational wave inference, PMC is approximately 50 and 25 times faster respectively than Nested Sampling (NS). We found that in higher dimensional applications the acceleration is even greater. Finally, PMC is directly parallelisable, manifesting linear scaling up to thousands of CPUs.


title: A spectral surrogate model for stochastic simulators computed from trajectory samples

abstract: Stochastic simulators are non-deterministic computer models which provide a different response each time they are run, even when the input parameters are held at fixed values. They arise when additional sources of uncertainty are affecting the computer model, which are not explicitly modeled as input parameters. The uncertainty analysis of stochastic simulators requires their repeated evaluation for different values of the input variables, as well as for different realizations of the underlying latent stochasticity. The computational cost of such analyses can be considerable, which motivates the construction of surrogate models that can approximate the original model and its stochastic response, but can be evaluated at much lower cost. We propose a surrogate model for stochastic simulators based on spectral expansions. Considering a certain class of stochastic simulators that can be repeatedly evaluated for the same underlying random event, we view the simulator as a random field indexed by the input parameter space. For a fixed realization of the latent stochasticity, the response of the simulator is a deterministic function, called trajectory. Based on samples from several such trajectories, we approximate the latter by sparse polynomial chaos expansion and compute analytically an extended Karhunen-Lo\`eve expansion (KLE) to reduce its dimensionality. The uncorrelated but dependent random variables of the KLE are modeled by advanced statistical techniques such as parametric inference, vine copula modeling, and kernel density estimation. The resulting surrogate model approximates the marginals and the covariance function, and allows to obtain new realizations at low computational cost. We observe that in our numerical examples, the first mode of the KLE is by far the most important, and investigate this phenomenon and its implications.


title: Spinless Mirror Chern Insulator from Projective Symmetry Algebra

abstract: It was commonly believed that a mirror Chern insulator (MCI) must require spin-orbital coupling, since time-reversal symmetry for spinless systems contradicts with the mirror Chern number. So MCI cannot be realized in spinless systems which include the large field of topological artificial crystals. Here, we disprove this common belief. The first point to clarify is that the fundamental constraint is not from spin-orbital coupling but the symmetry algebra of time reversal and mirror operations. Then, our theory is based on the conceptual transformation that the symmetry algebras will be projectively modified under gauge fields. Particularly, we show that the symmetry algebra of mirror reflection and time-reversal required for MCI can be achieved projectively in spinless systems with lattice $\mathbb{Z}_2$ gauge fields, i.e., by allowing real hopping amplitudes to take $\pm$ signs. Moreover, we propose the basic structure, the twisted $\pi$-flux blocks, to fulfill the projective symmetry algebra, and develop a general approach to construct spinless MCIs based on these building blocks. Two concrete spinless MCI models are presented, which can be readily realized in artificial systems such as acoustic crystals.


title: Exact solutions for the insulated and perfect conductivity problems with concentric balls

abstract: The insulated and perfect conductivity problems arising from high-contrast composite materials are considered in all dimensions. The solution and its gradient, respectively, represent the electric potential and field. The novelty of this paper lies in finding exact solutions for the insulated and perfect conductivity problems with concentric balls. Our results show that there appears no electric field concentration for the insulated conductivity problem, while the electric field for the perfect conductivity problem exhibits sharp singularity with respect to the small distance between interfacial boundaries of the interior and exterior balls. This discrepancy reveals that concentric balls is the optimal structure of insulated composites, but not for superconducting composites.


title: Phase estimation with limited coherence

abstract: We investigate the ultimate precision limits for quantum phase estimation in terms of the coherence, $C$, of the probe. For pure states, we give the minimum estimation variance attainable, $V(C)$, and the optimal state, in the asymptotic limit when the probe system size, $n$, is large. We prove that pure states are optimal only if $C$ scales as $n$ with a sufficiently large proportionality factor, and that the rank of the optimal state increases with decreasing $C$, eventually becoming full-rank. We show that the variance exhibits a Heisenberg-like scaling, $V(C) \sim a_n/C^2$, where $a_n$ decreases to $\pi^2/3$ as $n$ increases, leading to a dimension-independent relation.


title: Topological models to infer multiphase interstellar medium properties

abstract: (abridged) Spectroscopic observations of high-redshift galaxies slowly reveal the same complexity of the interstellar medium (ISM) as expected from resolved observations in nearby galaxies. While providing a wealth of diagnostics, such high-z spectra are often unresolved, and inferring reliable diagnostics represents a major obstacle. Bright, nearby, unresolved galaxies provide many constraints to design methods to infer ISM properties, but they have so far been limited to deterministic methods and/or with simple topological assumptions. Multicomponent ISM models associate a complex enough distribution of matter and phases with a simple enough topological description to be constrained with probabilistic frameworks. Our goal is to provide a flexible statistical framework that is agnostic to the model grid and that considers either a few discrete components defined by their parameter values and/or statistical distributions of parameters. In this study we present a the MULTIGRIS framework and use it to model a suite of spectral lines in the star-forming regions of the metal-poor dwarf galaxy IZw18. We are able to calculate accurate values for the metallicity, number of ionizing photons, masses of ionized and neutral hydrogen, as well as the dust mass and the dust-to-gas mass ratio. We find a relatively modest amount of H$_2$ which is predominantly CO-dark and traced by C$^+$ rather than C$^0$. Nevertheless, more than $90$% of the [CII] emission is associated with the neutral atomic gas. Our models confirm the necessity to consider an X-ray source. Finally, we investigate the escape fraction of ionizing photons for different energy ranges. While the escape fraction for the main HII region lies around $50-65$%, we show that most of the soft X-ray photons are able to escape and may play a role in the ionization and heating of the circumgalactic or intergalactic medium.