title: A GN-model closed-form formula supporting ultra-low fiber loss and short fiber spans

abstract: We report on a fully closed-form approximation of the GN model for Nyquist WDM systems that extends the range of applicability of previously available formulas to values of fiber loss and span loss that were not previously covered. In particular, so far, a closed-form formula was available for zero loss, and another for span loss of at least about 7 dB. The new formula is accurate over any fiber loss and any span loss value. The interest for this new formula is varied, ranging from hollow-core fibers to improving modeling of standard systems that present lumped loss along installed spans (such as splices or connectors).


title: Fast global spectral methods for three-dimensional partial differential equations

abstract: Global spectral methods offer the potential to compute solutions of partial differential equations numerically to very high accuracy. In this work, we develop a novel global spectral method for linear partial differential equations on cubes by extending ideas of Chebop2 [Townsend and Olver, J. Comput. Phys., 299 (2015)] to the three-dimensional setting utilizing expansions in tensorized polynomial bases. Solving the discretized PDE involves a linear system that can be recast as a linear tensor equation. Under suitable additional assumptions, the structure of these equations admits for an efficient solution via the blocked recursive solver [Chen and Kressner, Numer. Algorithms, 84 (2020)]. In the general case, when these assumptions are not satisfied, this solver is used as a preconditioner to speed up computations.


title: Recommendations on using VPN over SATCOM

abstract: VPN are a secured tunnel that help service providers to exchange data over non-secured networks. There is a large variety of VPN solutions that have variable deployment impacts on the target architecture as well as performance limitations or opportunities. This technical report compares Wireguard and OpenVPN for various SATCOM deployment scenarios and topologies.


title: Configurational Approaches to Urban Form: Empirical Test on the City of Nice (France)

abstract: Configurational analysis can be seen as the adaptation to urban space of complex network approaches already developed in social network analysis (Freeman 1979). Its central idea is that elements of urban form should be analysed with respect to the relationships that they establish with all other form elements within a given scale of analysis. In twenty years of research on urban form using configurational analysis, a certain variety of approaches and of more specific techniques has been developed. The main goal of this paper is to give a systematic overview of the different techniques and approaches and to test selected techniques in their ability to recognise different urban forms within the city of Nice, France.


title: Long Short-Term Memory Implementation Exploiting Passive RRAM Crossbar Array

abstract: The ever-increasing demand to extract temporal correlations across sequential data and perform context-based learning in this era of big data has led to the development of long short-term memory (LSTM) networks. Furthermore, there is an urgent need to perform these time-series data-dependent applications including speech/video processing and recognition, language modelling and translation, etc. on compact internet-of-things (IoT) edge devices with limited energy. To this end, in this work, for the first time, we propose an extremely area- and energy-efficient LSTM network implementation exploiting the passive resistive random access memory (RRAM) crossbar array. We developed a hardware-aware LSTM network simulation framework and performed an extensive analysis of the proposed LSTM implementation considering the non-ideal hardware artifacts such as spatial (device-to-device) and temporal variations, non-linearity, noise, etc. utilizing an experimentally calibrated comprehensive phenomenological model for passive RRAM crossbar array. Our results indicate that the proposed passive RRAM crossbar-based LSTM network implementation not only outperforms the prior digital and active 1T-1R crossbar-based LSTM implementations by more than three orders of magnitude in terms of area and two orders of magnitude in terms of training energy for identical network accuracy, but also exhibits robustness against spatial and temporal variations and noise, and a faster convergence rate. Our work may provide the incentive for experimental realization of LSTM networks on passive RRAM crossbar arrays.


title: An Improved Local Search Algorithm for k-Median

abstract: We present a new local-search algorithm for the $k$-median clustering problem. We show that local optima for this algorithm give a $(2.836+\epsilon)$-approximation; our result improves upon the $(3+\epsilon)$-approximate local-search algorithm of Arya et al. [STOC 01]. Moreover, a computer-aided analysis of a natural extension suggests that this approach may lead to an improvement over the best-known approximation guarantee for the problem. The new ingredient in our algorithm is the use of a potential function based on both the closest and second-closest facilities to each client. Specifically, the potential is the sum over all clients, of the distance of the client to its closest facility, plus (a small constant times) the truncated distance to its second-closest facility. We move from one solution to another only if the latter can be obtained by swapping a constant number of facilities, and has a smaller potential than the former. This refined potential allows us to avoid the bad local optima given by Arya et al. for the local-search algorithm based only on the cost of the solution.


title: Rashba and Dresselhaus effects in 2D Pb-I-based perovskites

abstract: Bulk hybride halide perovskites are governed by significant Rashba and Dresselhaus splitting. This indicates that such effects will not only affect their optoelectronic properties but also those of their two dimensional layered relatives. This work aims at understanding how different ways of symmetry breaking influence these effects in those materials. For this purpose, model structures are adopted where the organic compounds are replaced by Cs atoms. Disregarding possible distortions in the inorganic layers, results in structures with composition Cs$_{n+1}$Pb$_n$I$_{3n+1}$. Using the all-electron full-potential density-functional-theory code \texttt{exciting}, the impact of atomic displacement on the band structure is systematically studied for $n=1$, 2, 3 and $\infty$. The displacement patterns that yield Rashba or Dresselhaus splitting are identified, and the amount of the splitting is determined as a function of displacement. Furthermore, the spin textures in the electronic states around the band gap are analyzed to differentiate between Rashba and Dresselhaus effects. This study reveals in-plane Pb displacements as the origin of the strongest effects.


title: A note on the representations of $\text{SO}(1,d+1)$

abstract: $\text{SO}(1, d+1)$ is the isometry group of $(d+1)$-dimensional de Sitter spacetime $\text{dS}_{d+1}$ and the conformal group of $\mathbb{R}^{d}$. This note gives a pedagogical introduction to the representation theory of $\text{SO}(1, d+1)$, from the perspective of de Sitter quantum field theory and using tools from conformal field theory. Topics include (1) the construction and classification of all unitary irreducible representations (UIRs) of $\text{SO}(1,2)$ and $\text{SL}(2,\mathbb R)$, (2) the construction and classification of all UIRs of $\text{SO}(1,d+1)$ that describe integer-spin fields in $\text{dS}_{d+1}$, (3) a physical framework for understanding these UIRs, (4) the definition and derivation of Harish-Chandra group characters of $\text{SO}(1,d+1)$, and (5) a comparison between UIRs of $\text{SO}(1, d+1)$ and $\text{SO}(2,d)$.


title: Can Satellite Mega-Constellations Justify their Impact on Astronomy?

abstract: Private space companies are increasingly commercializing near-earth space. Some estimates show that more than 100,000 satellites may orbit the Earth by 2030. This increase in the number of satellites has many consequences for science, particularly astronomy. Satellite mega-constellations for satellite Internet connectivity are one of the main drivers behind the explosion in the number of satellites. Here, we briefly note whether such satellite mega-constellations can justify their impact on astronomy.


title: Frequency-Dependent F-Number Increases the Contrast and the Spatial Resolution in Fast Pulse-Echo Ultrasound Imaging

abstract: Fixed $F$-numbers reduce grating lobe artifacts in fast pulse-echo ultrasound imaging. Such $F$-numbers result in dynamic receive subapertures whose widths vary with the focal position. These subapertures, however, ignore useful low-frequency components in the excluded radio frequency (RF) signals and, thus, reduce the lateral resolution. Here, we propose a frequency-dependent $F$-number to simultaneously suppress grating lobe artifacts and maintain the lateral resolution. This $F$-number, at high frequencies, reduces the receive subaperture to remove spatially undersampled components of the RF signals and suppress grating lobes. The $F$-number, at low frequencies, enlarges the receive subaperture to use the components of all RF signals and maintain the lateral resolution. Experiments validated the proposed $F$-number and demonstrated improvements in the contrast and the widths of wire targets of up to 3.2 % and 12.8 %, respectively.


title: Migration of Jupiter mass planets in discs with laminar accretion flows

abstract: Migration of giant planets in discs with low viscosity has been studied recently. The proportionality between migration speed and the disc's viscosity is broken by the presence of vortices that appear at the edges of the planet-induced gap. Our goal is to investigate vortex-driven migration in low-viscosity discs in the presence of radial advection of gas, as expected from angular momentum removal by magnetised disc winds. We performed three dimensional simulations using the grid-based code FARGOCA. We mimicked the effects of a disc wind by applying a synthetic torque on a surface layer of the disc characterised by a prescribed column density Sigma_A so that it results in a disc accretion rate of 10^-8 Solar masses per year. Discs with this structure are called 'layered' and the layer where the torque is applied is denoted as 'active'. We also consider the case of accretion focussed near the disc midplane to mimic transport properties induced by a large Hall effect or by weak Ohmic diffusion. We observe two migration phases: in the first phase, the migration of the planet is driven by the vortex and is directed inwards. This phase ends when the vortex disappear. Migration depends on the ability of the torque from the planet to block the accretion flow. When the flow is fast and unimpeded (small Sigma_A) migration is very slow. When the accretion flow is completely blocked, migration is faster and the speed is controlled by the rate at which the accretion flow refills the gap behind the migrating planet. The migration speed of a giant planet in a layered protoplanetary disc depends on the thickness of the accreting layer. The lack of large-scale migration apparently experienced by the majority of giant exoplanets can be explained if the accreting layer is sufficiently thin to allow unimpeded accretion through the disc.


title: Graphs can be succinctly indexed for pattern matching in $ O(|E|^2 + |V|^{5 / 2}) $ time

abstract: For the first time we provide a succinct pattern matching index for arbitrary graphs that can be built in polynomial time, which requires less space and answers queries more efficiently than the one in [SODA 2021]. We show that, given an edge-labeled graph $ G = (V, E) $, there exists a data structures of $|E /_{\le_G}|(\lceil \log|\Sigma|\rceil + \lceil\log q\rceil + 2)\cdot (1+o(1)) + |V /_{\le_G}|\cdot (1+o(1))$ bits which can be built in $ O(|E|^2 + |V /_{\le_G}|^{5 / 2}) $ time and supports pattern matching on $ G $ in $O(|P| \cdot q^2 \cdot \log(q\cdot |\Sigma|))$ time, where $ G /_{\le_G} = (V /_{\le_G}, E /_{\le_G}) $ is a quotient graph obtained by collapsing some nodes in $ G $ (so $ |V /_{\le_G}| \le |V| $ and $ |E /_{\le_G}| \le |E| $) and $ q $ is the width of the maximum co-lex relation on $ G $. Our results have relevant applications in automata theory. First, we can build a succinct data structure to decide whether a string is accepted by a given automaton. Second, starting from an automaton $ \mathcal{A} $, one can define a relation $ \preceq_\mathcal{A} $ and a quotient automaton that capture the nondeterminism of $ \mathcal{A} $, improving the results in [SODA 2021].


title: Inertial Newton Algorithms Avoiding Strict Saddle Points

abstract: We study the asymptotic behavior of second-order algorithms mixing Newton's method and inertial gradient descent in non-convex landscapes. We show that, despite the Newtonian behavior of these methods, they almost always escape strict saddle points. We also evidence the role played by the hyper-parameters of these methods in their qualitative behavior near critical points. The theoretical results are supported by numerical illustrations.


title: Neyman-Pearson Multi-class Classification via Cost-sensitive Learning

abstract: Most existing classification methods aim to minimize the overall misclassification error rate. However, in applications such as loan default prediction, different types of errors can have varying consequences. To address this asymmetry issue, two popular paradigms have been developed: the Neyman-Pearson (NP) paradigm and the cost-sensitive (CS) paradigm. Previous studies on the NP paradigm have primarily focused on the binary case, while the multi-class NP problem poses a greater challenge due to its unknown feasibility. In this work, we tackle the multi-class NP problem by establishing a connection with the CS problem via strong duality and propose two algorithms. We extend the concept of NP oracle inequalities, crucial in binary classifications, to NP oracle properties in the multi-class context. Our algorithms satisfy these NP oracle properties under certain conditions. Furthermore, we develop practical algorithms to assess the feasibility and strong duality in multi-class NP problems, which can offer practitioners the landscape of a multi-class NP problem with various target error levels. Simulations and real data studies validate the effectiveness of our algorithms. To our knowledge, this is the first study to address the multi-class NP problem with theoretical guarantees. The proposed algorithms have been implemented in the R package \texttt{npcs}, which is available on CRAN.


title: 3D Hypersound Microscopy of Van der Waals Heterostructures

abstract: We employ here a picosecond ultrasonic technique to study Van der Waals heterostructures. Temporal variation of the reflection coefficient of the Al film that covers Van der Waals hBN/WSe$_2$/hBN heterostructures on a sapphire substrate after the femtosecond laser pulse excitation is carefully measured using an interferometric technique with spatial resolution. The laser pulse generates a broadband sound wave packet in aluminum film propagating perpendicular to the plane direction and partially reflecting from the heterostructural interfaces. The demonstrated technique has enough sensitivity to resolve a WSe$_2$ monolayer embedded in hBN. We apply a multilayered model of the optical and acoustical response that allows to evaluate the mechanical parameters, in particular, rigidity of interfaces, inaccessible from the other measurements. Mapping of the Fourier spectra of the response clearly visualizes different composition regions and can therefore serve as an acoustic tomography tool. Our findings demonstrate almost zero acoustic phonon dissipation below 150 GHz at the interfaces and in the layers that makes Van der Waals heterostructures perspective for nano-acoustical applications.


title: Thermodynamics of Precision in Markovian Open Quantum Dynamics

abstract: The thermodynamic and kinetic uncertainty relations indicate trade-offs between the relative fluctuation of observables and thermodynamic quantities such as dissipation and dynamical activity. Although these relations have been well studied for classical systems, they remain largely unexplored in the quantum regime. In this paper, we investigate such trade-off relations for Markovian open quantum systems whose underlying dynamics are quantum jumps, such as thermal processes and quantum measurement processes. Specifically, we derive finite-time lower bounds on the relative fluctuation of both dynamical observables and their first passage times for arbitrary initial states. The bounds imply that the precision of observables is constrained not only by thermodynamic quantities but also by quantum coherence. We find that the product of the relative fluctuation and entropy production or dynamical activity is enhanced by quantum coherence in a generic class of dissipative processes of systems with nondegenerate energy levels. Our findings provide insights into the survival of the classical uncertainty relations in quantum cases.


title: Rational Framing Motions and Spatial Rational Pythagorean Hodograph Curves

abstract: We propose a new method for constructing rational spatial Pythagorean Hodograph (PH) curves based on determining a suitable rational framing motion. While the spherical component of the framing motion is arbitrary, the translation part is determined be a modestly sized and nicely structured system of linear equations. Rather surprisingly, generic input data will only result in polynomial PH curves. We provide a complete characterization of all cases that admit truly rational (non-polynomial) solutions. Examples illustrate our ideas and relate them to existing literature.


title: On the Stochastic Stability of Deep Markov Models

abstract: Deep Markov models (DMM) are generative models that are scalable and expressive generalization of Markov models for representation, learning, and inference problems. However, the fundamental stochastic stability guarantees of such models have not been thoroughly investigated. In this paper, we provide sufficient conditions of DMM's stochastic stability as defined in the context of dynamical systems and propose a stability analysis method based on the contraction of probabilistic maps modeled by deep neural networks. We make connections between the spectral properties of neural network's weights and different types of used activation functions on the stability and overall dynamic behavior of DMMs with Gaussian distributions. Based on the theory, we propose a few practical methods for designing constrained DMMs with guaranteed stability. We empirically substantiate our theoretical results via intuitive numerical experiments using the proposed stability constraints.


title: Self-consistent description of high-spin states in doubly magic $^{208}$Pb

abstract: We analyze recent data on a long series of high-spin states in $^{208}$Pb with a self-consistent phonon-coupling model for nuclear excitations based on the Skyrme functionals. The model is the renormalized time-blocking approximation (RenTBA) which takes the coherent one-particle-one-hole (1p1h) states of the random-phase approximation (RPA) as starting point and develops from that more complex configurations beyond RPA. To the best of our knowledge, this is the first investigation of high spin states in $^{208}$Pb using self-consistent nuclear models. The interesting point here is that complex configurations are compulsory to describe the upper end of the long spin series at all. The data thus provide an ideal testing ground for phonon-coupling models as they give direct access to complex configurations. We find that standard Skyrme functionals which perform well in ground state properties and giant resonance excitations deliver at once an agreeable description of these high spin states.


title: Nuclear spin dynamics and noise in anisotropic large box model

abstract: We consider the central spin model in the box approximation taking into account external magnetic field and anisotropy of the hyperfine interaction. From the exact Hamiltonian diagonalization we obtain analytical expressions for the nuclear spin dynamics in the limit of many nuclear spins. We predict the nuclear spin precession in zero magnetic field for the case of anisotropic interaction between electron and nuclear spins. We calculate and describe the nuclear spin noise spectra in the thermodynamic equilibrium. The obtained results can be used for the analysis of the nuclear spin induced current fluctuations in organic semiconductors.


title: On the conjectured gravity-related collapse rate $E_\Delta/\hbar$ of massive quantum superpositions

abstract: Roger Penrose and the author share the proposal that the spatial superposition $|x_1\rangle+|x_2\rangle$ of a massive object collapses into its localized components $|x_1\rangle$ or $|x_2\rangle$ with the characteristic time $\hbar/E_\Delta$ where $E_\Delta$ is the gravitational self-energy excess of the superposition versus the localized states. Underlying arguments of such radical departure from standard quantum mechanics and different derivations of the rate equation are briefly recapitulated and discussed.


title: Study of dijet events with large rapidity separation in proton-proton collisions at $\sqrt{s} =$ 2.76 TeV

abstract: The cross sections for inclusive and Mueller-Navelet dijet production are measured as a function of the rapidity separation between the jets in proton-proton collisions at $\sqrt{s} =$ 2.76 TeV for jets with transverse momentum $p_\mathrm{T}$ $\gt$ 35 GeV and rapidity $\vert y\vert$ $\lt$ 4.7. Various dijet production cross section ratios are also measured. A veto on additional jets with $p_\mathrm{T}$ $\gt$ 20 GeV is introduced to improve the sensitivity to the effects of the Balitsky-Fadin-Kuraev-Lipatov (BFKL) evolution. The measurement is compared with the predictions of various Monte Carlo models based on leading-order and next-to-leading-order calculations including the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi leading-logarithm (LL) parton shower as well as the LL BFKL resummation.


title: Microstructure, grain boundary evolution and anisotropic Fe segregation in (0001) textured Ti thin films

abstract: The structure and chemistry of grain boundaries (GBs) are crucial in determining polycrystalline materials' properties. Faceting and solute segregation to minimize the GB energy is a commonly observed phenomenon. In this paper, a deposition process to obtain pure tilt GBs in titanium (Ti) thin films is presented. By increasing the power density, a transition from polycrystalline film growth to a maze bicrystalline Ti film on SrTiO$_3$ (001) substrate is triggered. All the GBs in the bicrystalline thin film are characterized to be $\Sigma$13 [0001] coincident site lattice (CSL) boundaries. The GB planes are seen to distinctly facet into symmetric {$\bar{7}520$} and asymmetric {$10\bar{1}0$} // {$11\bar{2}0$} segments of 20-50~nm length. Additionally, EDS reveals preferential segregation of iron (Fe) in every alternate symmetric {$\bar{7}520$} segment. Both the faceting and the segregation are explained by a difference in the CSL density between the facet planes. Furthermore, in the GB plane containing Fe segregation, atom probe tomography is used to experimentally determine the GB excess solute to be 1.25~atoms/nm$^{2}$. In summary, the study reveals for the first time a methodology to obtain bicrystalline Ti thin films with strong faceting and anisotropy in iron (Fe) segregation behaviour within the same family of planes.


title: A simple upper bound for Lebesgue constants associated with Leja points on the real line

abstract: Let $K\subset \mathbb R$ be a regular compact set and let $g(z)=g_{\overline{\mathbb C}\setminus K}(z,\infty)$ be the Green function for $\overline{\mathbb C}\setminus K$ with pole at infinity. For $\delta>0$, define $$ G(\delta):=\max\{ g(z): z\in \mathbb C, \,\operatorname{dist}(z,K)\le 2\delta\}. $$ Let $\{ x_n\}_{n=0}^\infty$ be a Leja sequence of points of $K$. Then the uniform norm $\|T_n\|=\Lambda_n, n=1,2,\ldots$ of the associated interpolation operator $T_n$, i.e., the $n$-th Lebesgue constant, is bounded from above by $$ \min_{\delta>0}2n\left[\frac{\operatorname{diam}( K)}{\delta}e^{nG(\delta)}\right]^{9/8}. $$ In particular, when $K$ is a uniformly perfect subset of $\mathbb R$, the Lebesgue constants grow at most polynomially in $n$. To the best of our knowledge, the result is new even when $K$ is a finite union of intervals.


title: Variance asymptotics and central limit theory for geometric functionals of Poisson cylinder processes

abstract: This paper deals with the union set of a stationary Poisson process of cylinders in $\mathbb{R}^n$ having an $(n-m)$-dimensional base and an $m$-dimensional direction space, where $m\in\{0,1,\ldots,n-1\}$ and $n\geq 2$. The concept simultaneously generalises those of a Boolean model and a Poisson hyperplane or $m$-flat process. Under very general conditions on the typical cylinder base a Berry-Esseen bound for the volume of the union set within a sequence of growing test sets is derived. Assuming convexity of the cylinder bases and of the window a similar result is shown for a broad class of geometric functionals, including the intrinsic volumes. In this context the asymptotic variance constant is analysed in detail, which in contrast to the Boolean model leads to a new degeneracy phenomenon. A quantitative central limit theory is developed in a multivariate set-up as well.


title: A Private and Computationally-Efficient Estimator for Unbounded Gaussians

abstract: We give the first polynomial-time, polynomial-sample, differentially private estimator for the mean and covariance of an arbitrary Gaussian distribution $\mathcal{N}(\mu,\Sigma)$ in $\mathbb{R}^d$. All previous estimators are either nonconstructive, with unbounded running time, or require the user to specify a priori bounds on the parameters $\mu$ and $\Sigma$. The primary new technical tool in our algorithm is a new differentially private preconditioner that takes samples from an arbitrary Gaussian $\mathcal{N}(0,\Sigma)$ and returns a matrix $A$ such that $A \Sigma A^T$ has constant condition number.


title: Sum-of-squares hierarchies for polynomial optimization and the Christoffel-Darboux kernel

abstract: Consider the problem of minimizing a polynomial $f$ over a compact semialgebraic set ${\mathbf{X} \subseteq \mathbb{R}^n}$. Lasserre introduces hierarchies of semidefinite programs to approximate this hard optimization problem, based on classical sum-of-squares certificates of positivity of polynomials due to Putinar and Schm\"udgen. When $\mathbf{X}$ is the unit ball or the standard simplex, we show that the hierarchies based on the Schm\"udgen-type certificates converge to the global minimum of $f$ at a rate in $O(1/r^2)$, matching recently obtained convergence rates for the hypersphere and hypercube $[-1,1]^n$. For our proof, we establish a connection between Lasserre's hierarchies and the Christoffel-Darboux kernel, and make use of closed form expressions for this kernel derived by Xu.


title: Safety Validation of Autonomous Vehicles using Assertion Checking

abstract: Safety and mission performance validation of autonomous vehicles (AVs) is a major challenge. In this paper we describe a methodology for constructing and applying assertion checks to validate the behaviour of an AV operating either in simulation or in the real world. We have identified a taxonomy of assertion types and the general format of their specification, and we have developed procedures for translating driving codes of practice to yield formal logical expressions that can be monitored automatically by computer, either by direct translation or by physical modelling. We have developed examples of assertions derived from the UK Highway Code (UKHC), as an example of a code of practice. We illustrate the approach with an example of assertion checking for vehicle overtaking, using a geospatial information system in an SQL database for validation and performance assessment. We present initial simulation and runtime monitoring experiments that apply assertions relevant in this overtaking scenario together with an analysis of the safety and mission performance characteristics measured.


title: Machine Learning Guided 3D Image Recognition for Carbonate Pore and Mineral Volumes Determination

abstract: Automated image processing algorithms can improve the quality, efficiency, and consistency of classifying the morphology of heterogeneous carbonate rock and can deal with a massive amount of data and images seamlessly. Geoscientists face difficulties in setting the direction of the optimum method for determining petrophysical properties from rock images, Micro-Computed Tomography (uCT), or Magnetic Resonance Imaging (MRI). Most of the successful work is from the homogeneous rocks focusing on 2D images with less focus on 3D and requiring numerical simulation. Currently, image analysis methods converge to three approaches: image processing, artificial intelligence, and combined image processing with artificial intelligence. In this work, we propose two methods to determine the porosity from 3D uCT and MRI images: an image processing method with Image Resolution Optimized Gaussian Algorithm (IROGA); advanced image recognition method enabled by Machine Learning Difference of Gaussian Random Forest (MLDGRF). We have built reference 3D micro models and collected images for calibration of IROGA and MLDGRF methods. To evaluate the predictive capability of these calibrated approaches, we ran them on 3D uCT and MRI images of natural heterogeneous carbonate rock. We measured the porosity and lithology of the carbonate rock using three and two industry-standard ways, respectively, as reference values. Notably, IROGA and MLDGRF have produced porosity results with an accuracy of 96.2% and 97.1% on the training set and 91.7% and 94.4% on blind test validation, respectively, in comparison with the three experimental measurements. We measured limestone and pyrite reference values using two methods, X-ray powder diffraction, and grain density measurements. MLDGRF has produced lithology (limestone and Pyrite) volumes with 97.7% accuracy.


title: Variational Automatic Curriculum Learning for Sparse-Reward Cooperative Multi-Agent Problems

abstract: We introduce a curriculum learning algorithm, Variational Automatic Curriculum Learning (VACL), for solving challenging goal-conditioned cooperative multi-agent reinforcement learning problems. We motivate our paradigm through a variational perspective, where the learning objective can be decomposed into two terms: task learning on the current task distribution, and curriculum update to a new task distribution. Local optimization over the second term suggests that the curriculum should gradually expand the training tasks from easy to hard. Our VACL algorithm implements this variational paradigm with two practical components, task expansion and entity progression, which produces training curricula over both the task configurations as well as the number of entities in the task. Experiment results show that VACL solves a collection of sparse-reward problems with a large number of agents. Particularly, using a single desktop machine, VACL achieves 98% coverage rate with 100 agents in the simple-spread benchmark and reproduces the ramp-use behavior originally shown in OpenAI's hide-and-seek project. Our project website is at https://sites.google.com/view/vacl-neurips-2021.


title: Learning Filterbanks for End-to-End Acoustic Beamforming

abstract: Recent work on monaural source separation has shown that performance can be increased by using fully learned filterbanks with short windows. On the other hand it is widely known that, for conventional beamforming techniques, performance increases with long analysis windows. This applies also to most hybrid neural beamforming methods which rely on a deep neural network (DNN) to estimate the spatial covariance matrices. In this work we try to bridge the gap between these two worlds and explore fully end-to-end hybrid neural beamforming in which, instead of using the Short-Time-Fourier Transform, also the analysis and synthesis filterbanks are learnt jointly with the DNN. In detail, we explore two different types of learned filterbanks: fully learned and analytic. We perform a detailed analysis using the recent Clarity Challenge data and show that by using learnt filterbanks it is possible to surpass oracle-mask based beamforming for short windows.


title: Safe Control of Arbitrary Nonlinear Systems using Dynamic Extension

abstract: Safe control for control-affine systems has been extensively studied. However, due to the complexity of system dynamics, it is challenging and time-consuming to apply these methods directly to non-control-affine systems, which cover a large group of dynamic systems, such as UAVs and systems with data-driven Neural Network Dynamic Models (NNDMs). Although all dynamic systems can be written in control-affine forms through dynamic extension, it remains unclear how to optimally design a computationally efficient algorithm to safely control the extended system. This paper addresses this challenge by proposing an optimal approach to synthesize safe control for the extended system under the framework of energy-function-based safe control. The proposed method first extends the energy function and then performs hyperparameter optimization to maximize performance while guaranteeing safety. It has been theoretically proved that our method guarantees safety (forward invariance of the safe set) and performance (bounded tracking error and smoother trajectories). It has been numerically validated that the proposed method is computationally efficient for non-control-affine systems.


title: Character Vectors of Strongly Regular Vertex Operator Algebras

abstract: We summarize interactions between vertex operator algebras and number theory through the lens of Zhu theory. The paper begins by recalling basic facts on vertex operator algebras (VOAs) and modular forms, and then explains Zhu's theorem on characters of VOAs in a slightly new form. We then axiomatize the desirable properties of modular forms that have played a role in Zhu's theorem and related classification results of VOAs. After this we summarize known classification results in rank two, emphasizing the geometric theory of vector-valued modular forms as a means for simplifying the discussion. We conclude by summarizing some known examples, and by providing some new examples, in higher ranks. In particular, the paper contains a number of potential character vectors that could plausibly correspond to a VOA, but such that the existence of a corresponding hypothetical VOA is presently unknown.


title: Isometry invariant shape recognition of projectively perturbed point clouds by the mergegram extending 0D persistence

abstract: Rigid shapes should be naturally compared up to rigid motion or isometry, which preserves all inter-point distances. The same rigid shape can be often represented by noisy point clouds of different sizes. Hence the isometry shape recognition problem requires methods that are independent of a cloud size. This paper studies stable-under-noise isometry invariants for the recognition problem stated in the harder form when given clouds can be related by affine or projective transformations. The first contribution is the stability proof for the invariant mergegram, which completely determines a single-linkage dendrogram in general position. The second contribution is the experimental demonstration that the mergegram outperforms other invariants in recognizing isometry classes of point clouds extracted from perturbed shapes in images.


title: A global shadow lemma and logarithm law for geometrically finite Hilbert geometries

abstract: For geometrically finite group actions on hyperbolic metric spaces and under certain assumptions on the growth of parabolic subgroups, we prove a global shadow lemma for Patterson-Sullivan measures, as well as a Dirichlet-type theorem and a logarithm law for excursion of geodesics into cusps. We then apply these results to geometrically finite quotients of strictly convex Hilbert geometries with C^1 boundary.


title: Efficient multi-partition topology optimization

abstract: In topology optimization, the state of structures is typically obtained by numerically evaluating a discretized PDE-based model. The degrees of freedom of such a model can be partitioned in free and prescribed sets to define the boundary conditions. A multi-partition problem involves multiple partitions of the same discretization, typically corresponding to different loading scenarios. As a result, solving multi-partition problems involves multiple factorization/preconditionings of the system matrix, requiring a high computational effort. In this paper, a novel method is proposed to efficiently calculate the responses and accompanying design sensitivities in such multi-partition problems using static condensation for use in gradient-based topology optimization. A main problem class that benefits from the proposed method is the topology optimization of small-displacement multi-input-multi-output compliant mechanisms. However, the method is applicable to any linear problem. We present its formulation and an algorithmic complexity analysis to estimate computational advantages for both direct and iterative solution methods to solve the system of equations, verified by numerical experiments. It is demonstrated that substantial gains are achievable for large-scale multi-partition problems. This is especially true for problems with both a small set of number of degrees of freedom that fully describes the performance of the structure and with large similarities between the different partitions. A major contribution to the gain is the lack of large adjoint analyses required to obtain the sensitivities of the performance measure.


title: A simple and versatile topology optimization formulation for flexure synthesis

abstract: High-tech equipment critically relies on flexures for precise manipulation and measurement. Through elastic deformation, flexures offer extreme position repeatability within a limited range of motion in their degrees of freedom, while constraining motion in the degrees of constraint. Topology optimization proves a prospective tool for the design of short-stroke flexures, providing maximum design freedom and allowing for application-specific requirements. State-of-the-art topology optimization formulations for flexure synthesis are subject to challenges like ease of use, versatility, implementation complexity, and computational cost, leaving a generally accepted formulation absent. This study proposes a novel topology optimization formulation for the synthesis of short-stroke flexures uniquely based on strain energy measures under prescribed displacement scenarios. The resulting self-adjoint optimization problem resembles great similarity to classic compliance minimization and inherits similar implementation simplicity, computational efficiency, and convergence properties. Numerical examples demonstrate the versatility in flexure types and the extendability of additional design requirements. The provided source code encourages the formulation to be explored and applied in academia and industry.


title: General-relativistic neutrino-radiation magnetohydrodynamic simulation of seconds-long black hole-neutron star mergers

abstract: Seconds-long numerical-relativity simulations for black hole-neutron star mergers are performed for the first time to obtain a self-consistent picture of the merger and post-merger evolution processes. To investigate the case that tidal disruption takes place, we choose the initial mass of the black hole to be $5.4M_\odot$ or $8.1M_\odot$ with the dimensionless spin of 0.75. The neutron-star mass is fixed to be $1.35M_\odot$. We find that after the tidal disruption, dynamical mass ejection takes place spending $\lesssim 10$ ms together with the formation of a massive accretion disk. Subsequently, the magnetic field in the disk is amplified by the magnetic winding and magnetorotational instability, establishing a turbulent state and inducing the angular momentum transport. The post-merger mass ejection by the magnetically-induced viscous effect sets in at $\sim 300$-500\,ms after the tidal disruption, at which the neutrino luminosity drops below $\sim 10^{51.5}\,{\rm erg/s}$, and continues for several hundreds ms. A magnetosphere near the rotational axis of the black hole is developed after the matter and magnetic flux fall into the black hole from the accretion disk, and high-intensity Poynting flux generation sets in at a few hundreds ms after the tidal disruption. The intensity of the Poynting flux becomes low after the significant post-merger mass ejection, because the opening angle of the magnetosphere increases. The lifetime for the stage with the strong Poynting flux is $1$-2 s, which agrees with the typical duration of short-hard gamma-ray bursts.


title: On $V$-filtration, Hodge filtration and Fourier transform

abstract: For $i : Z \to X$ a closed immersion of smooth varieties, we study how the $V$-filtration along $Z$ and the Hodge filtration on a mixed Hodge module $M$ on $X$ interact with each other. We also give a formula for the functors $i^*$, $i^!$ in terms of this $V$-filtration. As applications, we obtain results on the Hodge filtration of monodromic mixed Hodge modules and we give a Hodge theoretic proof of Skoda's theorem on multiplier ideals. Finally, we use the results to study the Fourier-Laplace transform of a monodromic mixed Hodge module.


title: Stripes in the extended $t-t^\prime$ Hubbard model: A Variational Monte Carlo analysis

abstract: By using variational quantum Monte Carlo techniques, we investigate the instauration of stripes (i.e., charge and spin inhomogeneities) in the Hubbard model on the square lattice at hole doping $\delta=1/8$, with both nearest- ($t$) and next-nearest-neighbor hopping ($t^\prime$). Stripes with different wavelengths $\lambda$ (denoting the periodicity of the charge inhomogeneity) and character (bond- or site-centered) are stabilized for sufficiently large values of the electron-electron interaction $U/t$. The general trend is that $\lambda$ increases going from negative to positive values of $t^\prime/t$ and decreases by increasing $U/t$. In particular, the $\lambda=8$ stripe obtained for $t^\prime=0$ and $U/t=8$ [L.F. Tocchio, A. Montorsi, and F. Becca, SciPost Phys. {\bf 7}, 21 (2019)] shrinks to $\lambda=6$ for $U/t\gtrsim 10$. For $t^\prime/t<0$, the stripe with $\lambda=5$ is found to be remarkably stable, while for $t^\prime/t>0$, stripes with wavelength $\lambda=12$ and $\lambda=16$ are also obtained. In all these cases, pair-pair correlations are highly suppressed with respect to the uniform state (obtained for large values of $|t^\prime/t|$), suggesting that striped states are not superconducting at $\delta=1/8$.


title: Optimal purification of a spin ensemble by quantum-algorithmic feedback

abstract: Purifying a high-temperature ensemble of quantum particles towards a known state is a key requirement to exploit quantum many-body effects. An alternative to passive cooling, which brings a system to its ground state, is based on feedback to stabilise the system actively around a target state. This alternative, if realised, offers additional control capabilities for the design of quantum states. Here we present a quantum feedback algorithm capable of stabilising the collective state of an ensemble from an infinite-temperature state to the limit of single quanta. We implement this on ~50,000 nuclei in a semiconductor quantum dot, and show that the nuclear-spin fluctuations are reduced 83-fold down to 10 spin macrostates. While our algorithm can purify a single macrostate, system-specific inhomogeneities prevent reaching this limit. Our feedback algorithm further engineers classically correlated ensemble states via macrostate tuning, weighted bimodality, and latticed multistability, constituting a pre-cursor towards quantum-correlated macrostates.


title: DeepSteal: Advanced Model Extractions Leveraging Efficient Weight Stealing in Memories

abstract: Recent advancements of Deep Neural Networks (DNNs) have seen widespread deployment in multiple security-sensitive domains. The need of resource-intensive training and use of valuable domain-specific training data have made these models a top intellectual property (IP) for model owners. One of the major threats to the DNN privacy is model extraction attacks where adversaries attempt to steal sensitive information in DNN models. Recent studies show hardware-based side channel attacks can reveal internal knowledge about DNN models (e.g., model architectures) However, to date, existing attacks cannot extract detailed model parameters (e.g., weights/biases). In this work, for the first time, we propose an advanced model extraction attack framework DeepSteal that effectively steals DNN weights with the aid of memory side-channel attack. Our proposed DeepSteal comprises two key stages. Firstly, we develop a new weight bit information extraction method, called HammerLeak, through adopting the rowhammer based hardware fault technique as the information leakage vector. HammerLeak leverages several novel system-level techniques tailed for DNN applications to enable fast and efficient weight stealing. Secondly, we propose a novel substitute model training algorithm with Mean Clustering weight penalty, which leverages the partial leaked bit information effectively and generates a substitute prototype of the target victim model. We evaluate this substitute model extraction method on three popular image datasets (e.g., CIFAR-10/100/GTSRB) and four DNN architectures (e.g., ResNet-18/34/Wide-ResNet/VGG-11). The extracted substitute model has successfully achieved more than 90 % test accuracy on deep residual networks for the CIFAR-10 dataset. Moreover, our extracted substitute model could also generate effective adversarial input samples to fool the victim model.


title: Procurements with Bidder Asymmetry in Cost and Risk-Aversion

abstract: We propose an empirical method to analyze data from first-price procurements where bidders are asymmetric in their risk-aversion (CRRA) coefficients and distributions of private costs. Our Bayesian approach evaluates the likelihood by solving type-symmetric equilibria using the boundary-value method and integrates out unobserved heterogeneity through data augmentation. We study a new dataset from Russian government procurements focusing on the category of printing papers. We find that there is no unobserved heterogeneity (presumably because the job is routine), but bidders are highly asymmetric in their cost and risk-aversion. Our counterfactual study shows that choosing a type-specific cost-minimizing reserve price marginally reduces the procurement cost; however, inviting one more bidder substantially reduces the cost, by at least 5.5%. Furthermore, incorrectly imposing risk-neutrality would severely mislead inference and policy recommendations, but the bias from imposing homogeneity in risk-aversion is small.


title: Test-beam demonstration of a TORCH prototype module

abstract: The TORCH time-of-flight detector is designed to provide a 15 ps timing resolution for charged particles, resulting in $\pi$/$K$ particle identification up to 10 GeV/c momentum over a 10 m flight path. Cherenkov photons, produced in a quartz plate of 10 mm thickness, are focused onto an array of micro-channel plate photomultipliers (MCP-PMTs) which measure the photon arrival times and spatial positions. A half-scale ($660\times1250\times10$ mm$^3$) TORCH demonstrator module has been tested in an 8 GeV/c mixed proton-pion beam at CERN. Customised square MCP-PMTs of active area $53\times53$ mm$^2$ and granularity $64\times64$ pixels have been employed, which have been developed in collaboration with an industrial partner. The single-photon timing performance and photon yields have been measured as a function of beam position in the radiator, giving measurements which are consistent with expectations. The expected performance of TORCH for high luminosity running of the LHCb Upgrade II has been simulated.


title: Accelerating GAN training using highly parallel hardware on public cloud

abstract: With the increasing number of Machine and Deep Learning applications in High Energy Physics, easy access to dedicated infrastructure represents a requirement for fast and efficient R&D. This work explores different types of cloud services to train a Generative Adversarial Network (GAN) in a parallel environment, using Tensorflow data parallel strategy. More specifically, we parallelize the training process on multiple GPUs and Google Tensor Processing Units (TPU) and we compare two algorithms: the TensorFlow built-in logic and a custom loop, optimised to have higher control of the elements assigned to each GPU worker or TPU core. The quality of the generated data is compared to Monte Carlo simulation. Linear speed-up of the training process is obtained, while retaining most of the performance in terms of physics results. Additionally, we benchmark the aforementioned approaches, at scale, over multiple GPU nodes, deploying the training process on different public cloud providers, seeking for overall efficiency and cost-effectiveness. The combination of data science, cloud deployment options and associated economics allows to burst out heterogeneously, exploring the full potential of cloud-based services.


title: A Pixelated Approach to Galaxy Catalogue Incompleteness: Improving the Dark Siren Measurement of the Hubble Constant

abstract: The use of gravitational wave standard sirens for cosmological analyses is becoming well known, with particular interest in measuring the Hubble constant, $H_0$, and in shedding light on the current tension between early- and late-time measurements. The current tension is over $4\sigma$ and standard sirens will be able to provide a completely independent measurement. Dark sirens (binary black hole or neutron star mergers with no electromagnetic counterparts) can be informative if the missing redshift information is provided through the use of galaxy catalogues to identify potential host galaxies of the merger. However, galaxy catalogue incompleteness affects this analysis, and accurate modelling of it is essential for obtaining an unbiased measurement of $H_0$. Previously this has been done by assuming uniform completeness within the sky area of a gravitational wave event. This paper presents an updated methodology in which the completeness of the galaxy catalogue is estimated in a directionally-dependent matter, by pixelating the sky and computing the completeness of the galaxy catalogue along each line of sight. The $H_0$ inference for a single event is carried out on a pixel-by-pixel basis, and the pixels are combined for the final result. A reanalysis of the events in the first gravitational wave transient catalogue (GWTC-1) leads to an improvement on the measured value of $H_0$ of approximately 5% compared to the 68.3% highest-density interval of the equivalent LIGO and Virgo result, with $H_0= 68.8^{+15.9}_{-7.8}$ km s$^{-1}$ Mpc$^{-1}$.


title: Strong convergence rate of Euler-Maruyama approximations in temporal-spatial H\"older-norms

abstract: Classical approximation results for stochastic differential equations analyze the $L^p$-distance between the exact solution and its Euler-Maruyama approximations. In this article we measure the error with temporal-spatial H\"older-norms. Our motivation for this are multigrid approximations of the exact solution viewed as a function of the starting point. We establish the classical strong convergence rate $0.5$ with respect to temporal-spatial H\"older-norms if the coefficient functions have bounded derivatives of first and second order.


title: SPT-SLIM: A Line Intensity Mapping Pathfinder for the South Pole Telescope

abstract: The South Pole Telescope Summertime Line Intensity Mapper (SPT-SLIM) is a pathfinder experiment that will demonstrate the use of on-chip filter-bank spectrometers for mm-wave line intensity mapping (LIM). The SPT-SLIM focal plane consists of 18 dual-polarization R=300 filter-bank spectrometers covering 120-180 GHz, coupled to aluminum kinetic inductance detectors. A compact cryostat holds the detectors at 100 mK and performs observations without removing the SPT-3G receiver. SPT-SLIM will be deployed to the 10-m South Pole Telescope for observations during the 2023-24 austral summer. We discuss the overall instrument design, expected detector performance and sensitivity to the LIM signal from CO at 0.5 < z < 2. The technology and observational techniques demonstrated by SPT-SLIM will enable next-generation LIM experiments that constrain cosmology beyond the redshift reach of galaxy surveys.


title: The Simulation and Design of an On-Chip Superconducting Millimetre Filter-Bank Spectrometer

abstract: Superconducting on-chip filter-banks provide a scalable, space saving solution to create imaging spectrometers at millimetre and sub-millimetre wavelengths. We present an easy to realise, lithographed superconducting filter design with a high tolerance to fabrication error. Using a capacitively coupled $\lambda/2$ microstrip resonator to define a narrow ($\lambda/\Delta\lambda = 300$) spectral pass band, the filtered output of a given spectrometer channel directly connects to a Lumped Element Kinetic Inductance Detector (LEKID). We show the tolerance analysis of our design, demonstrating $<11\%$ change in filter quality factor to any one realistic fabrication errors and a full filter-bank efficiency forecast to be 60\% after accounting for fabrication errors and dielectric loss tangent.


title: Design of SPT-SLIM focal plane; a spectroscopic imaging array for the South Pole Telescope

abstract: The Summertime Line Intensity Mapper (SLIM) is a mm-wave line-intensity mapping (mm-LIM) experiment for the South Pole Telescope (SPT). The goal of SPT-SLIM is to serve as a technical and scientific pathfinder for the demonstration of the suitability and in-field performance of multi-pixel superconducting filterbank spectrometers for future mm-LIM experiments. Scheduled to deploy in the 2023-24 austral summer, the SPT-SLIM focal plane will include 18 dual-polarization pixels, each coupled to an $R = \lambda/\Delta\lambda$ = 300 thin- film microstrip filterbank spectrometer that spans the 2 mm atmospheric window (120-180 GHz). Each individual spectral channel feeds a microstrip-coupled lumped-element kinetic inductance detector, which provides the highly multiplexed readout for the 10k detectors needed for SPT-SLIM. Here we present an overview of the preliminary design of key aspects of the SPT-SLIM the focal plane array, a description of the detector architecture and predicted performance, and initial test results that will be used to inform the final design of the SPT- SLIM spectrometer array.


title: Existence of Rotating Stars with Variable Entropy

abstract: We model a rotating star as a compressible fluid subject to gravitational forces. In almost all the mathematical literature the entropy is considered to be constant. Here we allow it to be variable. We consider a star that steadily rotates differentially around a fixed axis, say the $z$-axis. We prove the existence of a family of such stars with small angular velocity $\omega$ and small entropy variation $s$ and with an equation of state $p=Ke^s\rho^\gamma$. Our analysis reduces to a hyperbolic equation for the modified entropy coupled to an elliptic equation for the modified density, together with a mass constraint. Due to the variable entropy and the consequent loss of both regularity and variational structure, all the methods in the previous literature fail. We develop a new ad hoc perturbative strategy that allows us to construct rotating stars that bifurcate from the non-rotating ones.


title: CORE: a Complex Event Recognition Engine

abstract: Complex Event Recognition (CER) systems are a prominent technology for finding user-defined query patterns over large data streams in real time. CER query evaluation is known to be computationally challenging, since it requires maintaining a set of partial matches, and this set quickly grows super-linearly in the number of processed events. We present CORE, a novel COmplex event Recognition Engine that focuses on the efficient evaluation of a large class of complex event queries, including time windows as well as the partition-by event correlation operator. This engine uses a novel automaton-based evaluation algorithm that circumvents the super-linear partial match problem: under data complexity, it takes constant time per input event to maintain a data structure that compactly represents the set of partial matches and, once a match is found, the query results may be enumerated from the data structure with output-linear delay. We experimentally compare CORE against state-of-the-art CER systems on real-world data. We show that (1) CORE's performance is stable with respect to both query and time window size, and (2) CORE outperforms the other systems by up to five orders of magnitude on different workloads.


title: Improving the utility of locally differentially private protocols for longitudinal and multidimensional frequency estimates

abstract: This paper investigates the problem of collecting multidimensional data throughout time (i.e., longitudinal studies) for the fundamental task of frequency estimation under Local Differential Privacy (LDP) guarantees. Contrary to frequency estimation of a single attribute, the multidimensional aspect demands particular attention to the privacy budget. Besides, when collecting user statistics longitudinally, privacy progressively degrades. Indeed, the "multiple" settings in combination (i.e., many attributes and several collections throughout time) impose several challenges, for which this paper proposes the first solution for frequency estimates under LDP. To tackle these issues, we extend the analysis of three state-of-the-art LDP protocols (Generalized Randomized Response -- GRR, Optimized Unary Encoding -- OUE, and Symmetric Unary Encoding -- SUE) for both longitudinal and multidimensional data collections. While the known literature uses OUE and SUE for two rounds of sanitization (a.k.a. memoization), i.e., L-OUE and L-SUE, respectively, we analytically and experimentally show that starting with OUE and then with SUE provides higher data utility (i.e., L-OSUE). Also, for attributes with small domain sizes, we propose Longitudinal GRR (L-GRR), which provides higher utility than the other protocols based on unary encoding. Last, we also propose a new solution named Adaptive LDP for LOngitudinal and Multidimensional FREquency Estimates (ALLOMFREE), which randomly samples a single attribute to be sent with the whole privacy budget and adaptively selects the optimal protocol, i.e., either L-GRR or L-OSUE. As shown in the results, ALLOMFREE consistently and considerably outperforms the state-of-the-art L-SUE and L-OUE protocols in the quality of the frequency estimates.


title: A hemispheric two-channel code accounts for binaural unmasking in humans

abstract: Sound in noise is better detected or understood if target and masking sources originate from different locations. Mammalian physiology suggests that the neurocomputational process that underlies this binaural unmasking is based on two hemispheric channels that encode interaural differences in their relative neuronal activity. Here, we introduce a mathematical formulation of the two-channel model - the complex-valued correlation coefficient. We show that this formulation quantifies the amount of temporal fluctuations in interaural differences, which we suggest underlie binaural unmasking. We applied this model to an extensive library of psychoacoustic experiments, accounting for 98% of the variance across eight studies. Combining physiological plausibility with its success in explaining behavioral data, the proposed mechanism is a significant step towards a unified understanding of binaural unmasking and the encoding of interaural differences in general.


title: Evolution of the parent body of enstatite (EL) chondrites

abstract: Chondrites stem from undifferentiated asteroidal parent bodies that nevertheless experienced a certain degree of metamorphism after their formation in the early solar system. Maximum temperatures of metamorphism depend mainly on formation time and the abundance of the main heating source, which is short-lived 26Al. Enstatite chondrites formed under reducing conditions and include many strongly metamorphosed members of petrologic type 6. We model the thermal evolution of the parent body of the low metal enstatite chondrite class (EL). The model takes into account accretion, heating, sintering and compaction by hot pressing of the initially porous material, temperature dependent heat conductivity, and insulation effects by the remaining regolith layer. Key parameters of the parent body (formation time, radius, and porosity) are estimated by fitting thermal histories of EL6 chondrites (LON 94100, Neuschwanstein, Khairpur, Blithfield, Daniel's Kuil), constrained mainly by I-Xe and Ar-Ar ages and their respective closure temperatures. Viable fits are obtained for parent bodies with 120 -- 210 km radius, formed c. 1.8 -- 2.1 Ma after Ca,Al-rich inclusions (CAIs), and an initial porosity of 30%, relatively independent on initial disk temperatures. Thermal histories of the different EL6 chondrites are indistinguishable and so are their burial depths. While the exact layering depth is somewhat model dependent (c. 12--20 km), the thickness of the layer from which all five EL6 chondrites stem is < 1 km. Hence, an origin from a quite small asteroidal fragment is possible, particular as most excavation ages inferred from cosmic ray exposure data are compatible with a separation as meter sized meteoroids from a small Apollo asteroid 33 Ma ago.


title: S3RP: Self-Supervised Super-Resolution and Prediction for Advection-Diffusion Process

abstract: We present a super-resolution model for an advection-diffusion process with limited information. While most of the super-resolution models assume high-resolution (HR) ground-truth data in the training, in many cases such HR dataset is not readily accessible. Here, we show that a Recurrent Convolutional Network trained with physics-based regularizations is able to reconstruct the HR information without having the HR ground-truth data. Moreover, considering the ill-posed nature of a super-resolution problem, we employ the Recurrent Wasserstein Autoencoder to model the uncertainty.


title: Experiments conducted in the burning plasma regime with inertial fusion implosions

abstract: An experimental program is currently underway at the National Ignition Facility (NIF) to compress deuterium and tritium (DT) fuel to densities and temperatures sufficient to achieve fusion and energy gain. The primary approach being investigated is indirect drive inertial confinement fusion (ICF), where a high-Z radiation cavity (a hohlraum) is heated by lasers, converting the incident energy into x-ray radiation which in turn drives the DT fuel filled capsule causing it to implode. Previous experiments reported DT fuel gain exceeding unity [O.A. Hurricane et al., Nature 506, 343 (2014)] and then exceeding the kinetic energy of the imploding fuel [S. Le Pape et al., Phys. Rev. Lett. 120, 245003 (2018)]. We report on recent experiments that have achieved record fusion neutron yields on NIF, greater than 100 kJ with momentary fusion powers exceeding 1PW, and have for the first time entered the burning plasma regime where fusion alpha-heating of the fuel exceeds the energy delivered to the fuel via compression. This was accomplished by increasing the size of the high-density carbon (HDC) capsule, increasing energy coupling, while controlling symmetry and implosion design parameters. Two tactics were successful in controlling the radiation flux symmetry and therefore the implosion symmetry: transferring energy between laser cones via plasma waves, and changing the shape of the hohlraum. In conducting these experiments, we controlled for known sources of degradation. Herein we show how these experiments were performed to produce record performance, and demonstrate the data fidelity leading us to conclude that these shots have entered the burning plasma regime.


title: Statistical properties of large data sets with linear latent features

abstract: Analytical understanding of how low-dimensional latent features reveal themselves in large-dimensional data is still lacking. We study this by defining a linear latent feature model with additive noise constructed from probabilistic matrices, and analytically and numerically computing the statistical distributions of pairwise correlations and eigenvalues of the correlation matrix. This allows us to resolve the latent feature structure across a wide range of data regimes set by the number of recorded variables, observations, latent features and the signal-to-noise ratio. We find a characteristic imprint of latent features in the distribution of correlations and eigenvalues and provide an analytic estimate for the boundary between signal and noise even in the absence of a clear spectral gap.

categories: cond-mat.dis-nn math.ST physics.data-an stat.TH


title: Privacy-Preserving Distributed Average Consensus in Finite Time using Random Gossip

abstract: In this paper, we develop and analyze a gossip-based average consensus algorithm that enables all of the components of a distributed system, each with some initial value, to reach (approximate) average consensus on their initial values after executing a finite number of iterations, and without having to reveal the specific value they contribute to the average calculation. We consider a fully-connected (undirected) network in which each pair of components (nodes) can be randomly selected to perform pairwise standard gossip averaging of their values, and propose an enhancement that can be followed by each node that does not want to reveal its initial value to other (curious) nodes. We assume that curious nodes try to identify the initial values of other nodes but do not interfere in the computation in any other way; however, as a worst-case assumption, curious nodes are allowed to collaborate arbitrarily and are assumed to know the privacy-preserving strategy (but not the actual parameters chosen by the nodes that want to preserve their privacy). We characterize precisely conditions on the information exchange that guarantee privacy-preservation for a specific node. The protocol also provides a criterion that allows the nodes to determine, in a distributed manner (while running the enhanced gossip protocol), when to terminate their operation because approximate average consensus has been reached, i.e., all nodes have obtained values that are within a small distance from the exact average of their initial values.


title: (super)Schwarzian mechanics

abstract: In this paper we revisit the construction of supersymmetric Schwarzians using nonlinear realizations. We show that ${\cal N}=0,1,2,3,4$ supersymmetric Schwarzians can be systematically obtained as certain projections of Maurer-Cartan forms of superconformal groups after imposing simple conditions on them. We also present the supersymmetric Schwarzian actions defined as the integrals of products of Cartan forms. In contrast with the previous attempts to obtain the super-Schwarzians within nonlinear realizations technique, our set of constraints does not impose any restriction on the super-Schwarzians.


title: Surface Quasi-Geostrophic Equation driven by Space-Time White Noise

abstract: We consider the Surface Quasi-Geostrophic equation (SQG) driven by space-time white noise and show the existence of a local in time solution by applying the theory of regularity structures. A main difficulty is the presence of Riesz-transforms in the non-linearity. We show how to lift singular integral operators with a particular structure to the level of regularity structures and using this result we deduce the existence of a solution to SQG by a renormalisation procedure. The fact that the Riesz-transforms act only in the spatial variable makes it necessary to use inhomogeneous models instead of the standard ones.


title: Bayesian modelling of statistical region- and family-level clustered ordinal outcome data from Turkey

abstract: This study is concerned with the analysis of three-level ordinal outcome data with polytomous logistic regression in the presence of random-effects. It is assumed that the random-effects follow a Bridge distribution for the logit link, which allows one to obtain marginal interpretations of the regression coefficients. The data are obtained from the Turkish Income and Living Conditions Study, where the outcome variable is self-rated health (SRH), which is ordinal in nature. The analysis of these data is to compare covariate sub-groups and draw region- and family-level inferences in terms of SRH. Parameters and random-effects are sampled from the joint posterior densities following a Bayesian paradigm. Three criteria are used for model selection: Watenable information criterion, log pseudo marginal likelihood, and deviance information criterion. All three suggest that we need to account for both region- and family-level variabilities in order to model SRH. The extent to which the models replicate the observed data is examined by posterior predictive checks. Differences in SRH are found between levels of economic and demographic variables, regions of Turkey, and families who participated in the survey. Some of the interesting findings are that unemployed people are 19% more likely to report poorer health than employed people, and rural Aegean is the region that has the least probability of reporting poorer health.


title: Intelligent Reflecting Surfaces for Enhanced NOMA-based Visible Light Communications

abstract: The emerging intelligent reflecting surface (IRS) technology introduces the potential of controlled light propagation in visible light communication (VLC) systems. This concept opens the door for new applications in which the channel itself can be altered to achieve specific key performance indicators. In this paper, for the first time in the open literature, we investigate the role that IRSs can play in enhancing the link reliability in VLC systems employing non-orthogonal multiple access (NOMA). We propose a framework for the joint optimisation of the NOMA and IRS parameters and show that it provides significant enhancements in link reliability. The enhancement is even more pronounced when the VLC channel is subject to blockage and random device orientation.


title: Composition and Style Attributes Guided Image Aesthetic Assessment

abstract: The aesthetic quality of an image is defined as the measure or appreciation of the beauty of an image. Aesthetics is inherently a subjective property but there are certain factors that influence it such as, the semantic content of the image, the attributes describing the artistic aspect, the photographic setup used for the shot, etc. In this paper we propose a method for the automatic prediction of the aesthetics of an image that is based on the analysis of the semantic content, the artistic style and the composition of the image. The proposed network includes: a pre-trained network for semantic features extraction (the Backbone); a Multi Layer Perceptron (MLP) network that relies on the Backbone features for the prediction of image attributes (the AttributeNet); a self-adaptive Hypernetwork that exploits the attributes prior encoded into the embedding generated by the AttributeNet to predict the parameters of the target network dedicated to aesthetic estimation (the AestheticNet). Given an image, the proposed multi-network is able to predict: style and composition attributes, and aesthetic score distribution. Results on three benchmark datasets demonstrate the effectiveness of the proposed method, while the ablation study gives a better understanding of the proposed network.


title: Nanotubes as Sinks for Quantum Particles

abstract: Nanotubes with proper thickness, size and texture make ultra-efficient sinks for the quantum particles traveling into specific background media. Several optimal semiconducting cylindrical layers are reported to achieve enhancement in the trapping of matter waves by 2-3 orders of magnitude. The identified shells can be used as pieces in quantum devices that involve the focusing of incident beams from charge pumps and superconducting capacitors to radiation pattern controllers and matter-wave lenses.


title: Tidal evolution and diffusive growth during high-eccentricity planet migration: revisiting the eccentricity distribution of hot Jupiters

abstract: High-eccentricity tidal migration is a potential formation channel for hot Jupiters. During this process, the planetary f-mode may experience a phase of diffusive growth, allowing its energy to quickly build up to large values. In Yu et al. (2021, ApJ, 917, 31), we demonstrated that nonlinear mode interactions between a parent f-mode and daughter f- and p-modes expand the parameter space over which the diffusive growth of the parent is triggered. We extend that study by incorporating (1) the angular momentum transfer between the orbit and the mode, and consequently the evolution of the pericenter distance, (2) a phenomenological correction to the nonlinear frequency shift at high parent mode energies, and (3) dissipation of the parent's energy due to both turbulent convective damping of the daughter modes and strongly nonlinear wave-breaking events. The new ingredients allow us to follow the coupled evolution of the mode and orbit over $\gtrsim 10^4$ years, covering the diffusive evolution from its onset to its termination. We find that the semi-major axis shrinks by a factor of nearly ten over $10^4$ years, corresponding to a tidal quality factor $\mathcal{Q}\sim10$. The f-mode's diffusive growth terminates while the eccentricity is still high, at around e=0.8-0.95. Using these results, we revisit the eccentricity distribution of proto-hot Jupiters. We estimate that less than 1 proto-HJ with eccentricity >0.9 should be expected in Kepler's data once the diffusive regime is accounted for, explaining the observed paucity of this population.


title: Minimum-lap-time Control Strategies for All-wheel Drive Electric Race Cars via Convex Optimization

abstract: This paper presents a convex optimization framework to compute the minimum-lap-time control strategies of all-wheel drive (AWD) battery electric race cars, accounting for the grip limitations of the individual tyres. Specifically, we first derive the equations of motion (EOM) of the race car and simplify them to a convex form. Second, we leverage convex models of the electric motors (EMs) and battery, and frame the time-optimal final-drives design and EMs control problem in space domain. The resulting optimization problem is fully convex and can be efficiently solved with global optimality guarantees using second-order conic programming algorithms. Finally, we validate our modeling assumptions via the original non-convex EOM, and simulate our framework on the Formula Student Netherlands endurance race track. Thereby, we compare a torque vectoring with a fixed power split configuration, showing that via torque vectoring we can make a better use of the individual tyre grip, and significantly improve the achievable lap time by more than 4%. Finally, we present a design study investigating the respective impact of the front and rear EM size on lap time, revealing that the rear motor sizing is predominant due to the higher vertical rear tyre load caused by the center of pressure position and rearwards load transfer under acceleration.


title: Saturation of the asymmetric current filamentation instability under conditions relevant to relativistic shock precursors

abstract: The current filamentation instability, which generically arises in the counterstreaming of supersonic plasma flows, is known for its ability to convert the free energy associated with anisotropic momentum distributions into kinetic-scale magnetic fields. The saturation of this instability has been extensively studied in symmetric configurations where the interpenetrating plasmas share the same properties (velocity, density, temperature). In many physical settings, however, the most common configuration is that of asymmetric plasma flows. For instance, the precursor of relativistic collisionless shock waves involves a hot, dilute beam of accelerated particles reflected at the shock front and a cold, dense inflowing background plasma. To determine the appropriate criterion for saturation in this case, we have performed large-scale 2D particle-in-cell simulations of counterstreaming electron-positron pair and electron-ion plasmas. We show that, in interpenetrating pair plasmas, the relevant criterion is that of magnetic trapping as applied to the component (beam or plasma) that carries the larger inertia of the two; namely, the instability growth suddenly slows down once the quiver frequency of those particles equals or exceeds the instability growth rate. We present theoretical approximations for the saturation level. These findings remain valid for electron-ion plasmas provided that electrons and ions are close to equipartition in the plasma flow of larger inertia. Our results can be directly applied to the physics of relativistic, weakly magnetized shock waves, but they can also be generalized to other cases of study.


title: Optimal convex lifted sparse phase retrieval and PCA with an atomic matrix norm regularizer

abstract: We present novel analysis and algorithms for solving sparse phase retrieval and sparse principal component analysis (PCA) with convex lifted matrix formulations. The key innovation is a new mixed atomic matrix norm that, when used as regularization, promotes low-rank matrices with sparse factors. We show that convex programs with this atomic norm as a regularizer provide near-optimal sample complexity and error rate guarantees for sparse phase retrieval and sparse PCA. While we do not know how to solve the convex programs exactly with an efficient algorithm, for the phase retrieval case we carefully analyze the program and its dual and thereby derive a practical heuristic algorithm. We show empirically that this practical algorithm performs similarly to existing state-of-the-art algorithms.


title: Adaptive Steganography Based on bargain Game

abstract: The capacity and security of the confidential message on the channel are two important challenges in steganography. In this paper, a new block steganography model is presented using the bargain method so that a competitive model is introduced. In this game, the blocks are the same players. The bargain is provided with the aim of embedding information without reducing capacity as well as increasing security. The proposed model shows that it can be used both of the special domain and the transform domain, which are two important methods of steganography. For this purpose, an example of a special domain model is introduced in which, In the first step, the image is divided into $n \times n$ blocks, and in the second step using the graph coloring algorithm, pixels are considered to embed confidential information in each block. In the third step, regarding the bargaining method in game theory, each block plays the role of a player, that the competition between players is based on the defined goal function, and in the best blocks in terms of two criteria of capacity and security, which here means each block has a higher security-to-capacity ratio, so it has a higher priority, which is determined based on the bargaining model. Also, information embedded in LSB two bits. An example of a conversion domain method is also shows that security increases without decreasing in capacity. The conclusion is evaluated by three criteria: PSNR, histogram, and $\epsilon-secure$ also, 2000 standard images were evaluated and observed that the proposed method improves the block methods of embedding information.


title: Search for Quantum Gravity Using Astrophysical Neutrino Flavour with IceCube

abstract: Along their long propagation from production to detection, neutrino states undergo quantum interference which converts their types, or flavours. High-energy astrophysical neutrinos, first observed by the IceCube Neutrino Observatory, are known to propagate unperturbed over a billion light years in vacuum. These neutrinos act as the largest quantum interferometer and are sensitive to the smallest effects in vacuum due to new physics. Quantum gravity (QG) aims to describe gravity in a quantum mechanical framework, unifying matter, forces and space-time. QG effects are expected to appear at the ultra-high-energy scale known as the Planck energy, $E_{P}\equiv 1.22\times 10^{19}$~giga-electronvolts (GeV). Such a high-energy universe would have existed only right after the Big Bang and it is inaccessible by human technologies. On the other hand, it is speculated that the effects of QG may exist in our low-energy vacuum, but are suppressed by the Planck energy as $E_{P}^{-1}$ ($\sim 10^{-19}$~GeV$^{-1}$), $E_{P}^{-2}$ ($\sim 10^{-38}$~GeV$^{-2}$), or its higher powers. The coupling of particles to these effects is too small to measure in kinematic observables, but the phase shift of neutrino waves could cause observable flavour conversions. Here, we report the first result of neutrino interferometry~\cite{Aartsen:2017ibm} using astrophysical neutrino flavours to search for new space-time structure. We did not find any evidence of anomalous flavour conversion in IceCube astrophysical neutrino flavour data. We place the most stringent limits of any known technologies, down to $10^{-42}$~GeV$^{-2}$, on the dimension-six operators that parameterize the space-time defects for preferred astrophysical production scenarios. For the first time, we unambiguously reach the signal region of quantum-gravity-motivated physics.


title: Phase separation dynamics in deformable droplets

abstract: Phase separation can drive spatial organization of multicomponent mixtures. For instance in developing animal embryos, effective phase separation descriptions have been used to account for the spatial organization of different tissue types. Similarly, separation of different tissue types and the emergence of a polar organization is also observed in cell aggregates mimicking early embryonic axis formation. Here, we describe such aggregates as deformable two-phase fluid droplets, which are suspended in a fluid environment (third phase). Using hybrid finite-volume Lattice-Boltzmann simulations, we numerically explore the out-of-equilibrium routes that can lead to the polar equilibrium state of such a droplet (Janus droplet). We focus on the interplay between spinodal decomposition and advection with hydrodynamic flows driven by interface tensions, which we characterize by a Peclet number $Pe$. Consistent with previous work, for large $Pe$ the coarsening process is generally accelerated. However, for intermediate $Pe$ we observe long-lived, strongly elongated droplets, where both phases form an alternating stripe pattern. We show that these ``croissant'' states are close to mechanical equilibrium and coarsen only slowly through diffusive fluxes in an Ostwald-ripening-like process. Finally, we show that a surface tension asymmetry between both droplet phases leads to transient, rotationally symmetric states whose resolution leads to flows reminiscent of Marangoni flows. Our work highlights the importance of advection for the phase separation process in finite, deformable systems.


title: Newly born extragalactic millisecond pulsars as efficient emitters of PeV neutrinos

abstract: The origins of the diffuse flux of cosmogenic PeV neutrinos detected in the IceCube experiment during 2010 - 2017 remain unidentified. A population of extragalactic newly born fast spinning pulsars are investigated as possible candidates for generating the PeV energy scale neutrinos. A two-step mechanism of particle acceleration is adopted for transferring energy from star rotation to high energy electrons. Electrons might be boosted up to $\approx 0.01$ EeV energies and above, and produce ultra-high-energy (UHE) neutrinos and gamma rays when these electrons interact with low energy positrons and soft radiations in the acceleration zone. The theoretically derived extragalactic diffuse muon neutrino flux in the energy range [1-10] PeV is found consistent with the IceCube level if only a fraction ($\eta_{k} \approx 0.31\%$) of the total bolometric luminosity of pulsars are transferred to power the PeV neutrinos. Using the above value of the cosmic ray electron loading parameter $\eta_{k}$, the diffuse gamma ray flux from inverse Compton scattering off UHE electrons in the soft radiation field might be predicted from the model.


title: Distribution efficiency and structure of complex networks

abstract: Flow networks efficiently transport nutrients and other solutes in many physical systems, such as plant and animal vasculature. In the case of the animal circulatory system, an adequate oxygen and nutrient supply is not guaranteed everywhere: as nutrients travel through the microcirculation and get absorbed, they become less available at the venous side of the vascular network. Ensuring that the nutrient distribution is homogeneous provides a fitness advantage, as all tissue gets enough supply to survive while waste is minimized. How do animals build such a uniform perfusing flow system? We propose a local adaptation rule for the vessel radii that is able to equalize perfusion, while minimizing energy dissipation to circulate the flow and the material cost. The rule is a combination of different objective cost functions that compete to produce complex network morphologies ranging from hierarchical architectures to uniform mesh grids, depending on how each cost is weighted. We find that our local adaptation rules are consistent with experimental data of the rat mesentery vasculature.


title: Consistent Sufficient Explanations and Minimal Local Rules for explaining regression and classification models

abstract: To explain the decision of any model, we extend the notion of probabilistic Sufficient Explanations (P-SE). For each instance, this approach selects the minimal subset of features that is sufficient to yield the same prediction with high probability, while removing other features. The crux of P-SE is to compute the conditional probability of maintaining the same prediction. Therefore, we introduce an accurate and fast estimator of this probability via random Forests for any data $(\boldsymbol{X}, Y)$ and show its efficiency through a theoretical analysis of its consistency. As a consequence, we extend the P-SE to regression problems. In addition, we deal with non-discrete features, without learning the distribution of $\boldsymbol{X}$ nor having the model for making predictions. Finally, we introduce local rule-based explanations for regression/classification based on the P-SE and compare our approaches w.r.t other explainable AI methods. These methods are available as a Python package at \url{www.github.com/salimamoukou/acv00}.


title: XXL-HSC: An updated catalogue of high-redshift ($z\geqslant3.5$) X-ray AGN in the XMM-XXL northern field. Constraints on the bright end of the soft logN-logS

abstract: X-rays offer a reliable method to identify Active Galactic Nuclei (AGN). However, in the high-redshift Universe X-ray AGN are poorly sampled due to their relatively low space density and the small areas covered by X-ray surveys. In addition to wide-area X-ray surveys, it is important to have deep optical data in order to locate the optical counterparts and determine their redshifts. In this work, we build a high-redshift ($z\geqslant3.5$) X-ray selected AGN sample in the XMM-XXL northern field using the most updated [0.5-2 keV] catalogue along with a plethora of new spectroscopic and multi-wavelength catalogues, including the deep optical Subaru Hyper Suprime-Cam (HSC) data reaching magnitude limits i$\sim$26 mag. We select all the spectroscopically confirmed AGN and complement this sample with high-redshift candidates that are HSC g- and r-band dropouts. To confirm the dropouts, we derive their photometric redshifts using spectral energy distribution techniques. We end up with a sample of 54 high-z sources (28 with spec-z), the largest in this field so far (almost three times larger than in previous studies), and we estimate the possible contamination and completeness. We calculate the number counts (logN-logS) in different redshift bins and compare our results with previous studies and models. We provide the strongest high-redshift AGN constraints yet at bright fluxes ($\rm f_{0.5-2~keV} > 10^{-15} erg~s^{-1}~cm^{-2}$). The samples of $z\geqslant3.5$, $z\geqslant4$ and $z\geqslant5$ are in agreement with an exponential decline model similar to that witnessed at optical wavelengths. Our work emphasizes the importance of using wide area X-ray surveys with deep optical data to uncover high redshift AGN.


title: First-degree prime ideals of composite extensions

abstract: Let $\mathbb{Q}(\alpha)$ and $\mathbb{Q}(\beta)$ be linearly disjoint number fields and let $\mathbb{Q}(\theta)$ be their compositum. We prove that the first-degree prime ideals of $\mathbb{Z}[\theta]$ may almost always be constructed in terms of the first-degree prime ideals of $\mathbb{Z}[\alpha]$ and $\mathbb{Z}[\beta]$, and vice-versa. We also classify the cases in which this correspondence does not hold, by providing explicit counterexamples. We show that for every pair of coprime integers $d,e \in \mathbb{Z}$, such a correspondence almost always respects the divisibility of principal ideals of the form $(e+d\theta)\mathbb{Z}[\theta]$, with a few exceptions that we characterize. Finally, we discuss the computational improvement of such an approach, and we verify the reduction in time needed for computing such primes for certain concrete cases.


title: Higher Order $c$-Differentials

abstract: EFRST20, the notion of $c$-differentials was introduced as a potential expansion of differential cryptanalysis against block ciphers utilizing substitution boxes. Drawing inspiration from the technique of higher order differential cryptanalysis, in this paper we propose the notion of higher order $c$-derivatives and differentials and investigate their properties. Additionally, we consider how several classes of functions, namely the multiplicative inverse function and the Gold function, perform under higher order $c$-differential uniformity.


title: Genus-zero Permutation-twisted Conformal Blocks for Tensor Product Vertex Operator Algebras: The Tensor-factorizable Case

abstract: For a vertex operator algebra $V$, we construct an explicit isomorphism between the space of genus-0 conformal blocks associated to permutation-twisted $V^{\otimes n}$-modules and the space of conformal blocks associated to untwisted $V$-modules and a branched covering C of the Riemann sphere. As a consequence, when V is CFT-type, rational, and C2 cofinite, the fusion rules for permutation-twisted modules are determined. We also relate the sewing and factorization of permutation-twisted $V^{\otimes n}$-conformal blocks and untwisted $V$-conformal blocks. Various applications are discussed.


title: Wasserstein Distributionally Robust Optimization with Expected Value Constraints

abstract: We investigate a stochastic program with expected value constraints, addressing the problem in a general context through Distributionally Robust Optimization (DRO) approach using Wasserstein distances, where the ambiguity set depends on the decision variable. We demonstrate that this approach can be reformulated into a finite-dimensional optimization problem, which, in certain instances, can be convex. Moreover, we establish criteria for determining the feasibility of the problem concerning the Wasserstein radius and the parameter governing the constraint. Finally, we present numerical results within the context of portfolio optimization. In particular, we highlight the distinctions between our approach and several existing non-robust methods, using both simulated data and real financial market data.


title: Drag anisotropy of cylindrical solids in fluid-saturated granular beds

abstract: We study the direction-dependent drag acting on a cylindrical solid intruder with length $L$ and diameter $D$ as it moves in water-saturated granular beds at constant depth. Polysterene and hydrogel grains with diameter $d$ are used to investigate materials which have high contact friction and those which are nearly frictionless, respectively. The drag on the intruder is measured while oriented perpendicular $F_\perp$ and parallel $F_\parallel$ to its axis as a function of speed $U$ from the quasi-static to the rate-dependent regime. We find that the drag anisotropy $\xi = F_\perp/F_\parallel$ is not constant, and increases significantly with driving rate and $L/D$ in both mediums. In particular for $L/D = 40$, $\xi$ increases from 2.6 to 4.5, and from 7.0 to 8.2 in the high and low friction beds, respectively, as the nondimensional Froude number $Fr =U/\sqrt{g(D+d)}$ is varied between $10^{-4}$ and $2 \times 10^{-2}$. On average, $\xi$ is observed to increase logarithmically with $L/D$ for $L/D \gg 1$. Exploiting the index-matched nature of hydrogel grains in water, we further show that the sediment flow around the cylinder in the two orientations is consistent with skin friction dominated drag. The relative contributions of the cylindrical side and the circular flat-ends on $\xi$ are estimated with thin disks to understand the observed variation of drag with aspect ratio and surface friction.


title: Evaluating Predictive Uncertainty and Robustness to Distributional Shift Using Real World Data

abstract: Most machine learning models operate under the assumption that the training, testing and deployment data is independent and identically distributed (i.i.d.). This assumption doesn't generally hold true in a natural setting. Usually, the deployment data is subject to various types of distributional shifts. The magnitude of a model's performance is proportional to this shift in the distribution of the dataset. Thus it becomes necessary to evaluate a model's uncertainty and robustness to distributional shifts to get a realistic estimate of its expected performance on real-world data. Present methods to evaluate uncertainty and model's robustness are lacking and often fail to paint the full picture. Moreover, most analysis so far has primarily focused on classification tasks. In this paper, we propose more insightful metrics for general regression tasks using the Shifts Weather Prediction Dataset. We also present an evaluation of the baseline methods using these metrics.


title: Machine Learning-based Test Selection for Simulation-based Testing of Self-driving Cars Software

abstract: Abstract Simulation platforms facilitate the development of emerging cyber-physical systems (CPS) like self-driving cars (SDC) because they are more efficient and less dangerous than field operational tests. Despite this, thoroughly testing SDCs in simulated environments remains challenging because SDCs must be tested in a sheer amount of long-running test scenarios. Past results on software testing optimization have shown that not all the tests contribute equally to establishing confidence in test subjects' quality and reliability, with some \uninformative" tests that can be skipped (or removed) to reduce testing effort. However, this problem was partially addressed in the context of SDC simulation platforms. In this paper, we investigate test selection strategies to increase the cost-effectiveness of simulation-based testing in the context of SDCs. We propose an approach called SDC-Scissor (SDC coSt-effeCtIve teSt SelectOR), which leverages machine learning (ML) strategies to identify and skip tests that are unlikely to detect faults in SDCs before executing them. Specifically, SDC-Scissor extract features concerning the characteristics of the test scenarios being executed in the simulation environment and via ML strategies predict tests that lead to faults before executing them. Our evaluation shows that SDC-Scissor achieved high classification accuracy (up to 93.4%) in classifying tests leading to a fault which allows improving testing cost-effectiveness: SDC-Scissor was able to reduce (ca. 170%) the time spent in running irrelevant tests as well as identified 33% more failure triggering tests compared to a randomized baseline. Interestingly, SDC-Scissor does not introduce significant computational overhead in the SDCs testing process, which is critical to SDC development in industrial settings.


title: Comment on "Using an atom interferometer to infer gravitational entanglement generation''

abstract: Our paper arXiv:2101.11629 contains a technical error which changes some of the conclusions. We thank Streltsov, Pedernales, and Plenio for bringing the essence of this error to our attention. Here we explain the error, examine its consequences, and suggest methods to overcome the resulting weakness in the proposed experiment.


title: $\pi-\pi$ scattering, QED and finite-volume quantization

abstract: Using the Coulomb gauge formulation of QED we present a lattice QCD procedure to calculate the $\pi^+\pi^+$ scattering phase shift including the effects of the Coulomb potential which appears in this formulation. The approach described here incorporates the effects of relativity and avoids finite-volume corrections that vanish as a power of the volume in which the lattice calculation is performed. This is the first step in developing a complete lattice QCD calculation of the electromagnetic and isospin-breaking light-quark mass contributions to $\varepsilon'$, the parameter describing direct CP violating effects in $K_L\to\pi\pi$ decay.