title: Toward fully compressible numerical simulations of stellar magneto-convection with the RAMSES code

abstract: Numerical simulations of magneto-convection have greatly expanded our understanding of stellar interiors and stellar magnetism. Recently, fully compressible hydrodynamical simulations of full-star models have demonstrated the feasibility of studying the excitation and propagation of pressure and internal gravity waves in stellar interiors, which would allow for a direct comparison with asteroseismological measurements. However, the impact of magnetic fields on such waves has not been taken into account yet in three-dimensional simulations. We conduct a proof of concept for the realization of three-dimensional, fully compressible, magneto-hydrodynamical numerical simulations of stellar interiors with the RAMSES code. We adapted the RAMSES code to deal with highly subsonic turbulence, typical of stellar convection, by implementing a well-balanced scheme in the numerical solver. We then ran and analyzed three-dimensional hydrodynamical and magneto-hydrodynamical simulations with different resolutions of a plane-parallel convective envelope on a Cartesian grid. Both hydrodynamical and magneto-hydrodynamical simulations develop a quasi-steady, turbulent convection layer from random density perturbations introduced over the initial profiles. The convective flows are characterized by small-amplitude fluctuations around the hydrodynamical equilibrium of the stellar interior, which is preserved over the whole simulation time. Using our compressible well-balanced scheme, we were able to model flows with Mach numbers as low as $\mathcal{M} \sim 10^{-3}$, but even lower Mach number flows are possible in principle. In the magneto-hydrodynamical runs, we observe an exponential growth of magnetic energy consistent with the action of a small-scale dynamo. (Abridged)


title: reXcor: A Model of the X-ray Spectrum of Active Galactic Nuclei that Combines Ionized Reflection and a Warm Corona

abstract: The X-ray spectra of active galactic nuclei (AGNs) often exhibit an excess of emission above the primary power-law at energies <~ $2$ keV. Two models for the origin of this `soft excess' are ionized relativistic reflection from the inner accretion disc and Comptonization of thermal emission in a warm corona. Here, we introduce reXcor, a new AGN X-ray ($0.3$-$100$ keV) spectral fitting model that self-consistently combines the effects of both ionized relativistic reflection and the emission from a warm corona. In this model, the accretion energy liberated in the inner disc is distributed between a warm corona, a lamppost X-ray source, and the accretion disc. The emission and ionized reflection spectrum from the inner $400$ $r_g$ of the disc is computed, incorporating the effects of relativistic light-bending and blurring. The resulting spectra predict a variety of soft excess shapes and sizes that depend on the fraction of energy dissipated in the warm corona and lamppost. We illustrate the use of reXcor by fitting to the joint XMM-Newton and NuSTAR observations of the Seyfert 1 galaxies HE 1143-1820 and NGC 4593, and find that both objects require a warm corona contribution to the soft excess. Eight reXcor table models, covering different values of accretion rate, lamppost height and black hole spin, are publicly available through the XSPEC website. Systematic use of reXcor will provide insight into the distribution of energy in AGN accretion flows.


title: Tailored max-out networks for learning convex PWQ functions

abstract: Convex piecewise quadratic (PWQ) functions frequently appear in control and elsewhere. For instance, it is well-known that the optimal value function (OVF) as well as Q-functions for linear MPC are convex PWQ functions. Now, in learning-based control, these functions are often represented with the help of artificial neural networks (NN). In this context, a recurring question is how to choose the topology of the NN in terms of depth, width, and activations in order to enable efficient learning. An elegant answer to that question could be a topology that, in principle, allows to exactly describe the function to be learned. Such solutions are already available for related problems. In fact, suitable topologies are known for piecewise affine (PWA) functions that can, for example, reflect the optimal control law in linear MPC. Following this direction, we show in this paper that convex PWQ functions can be exactly described by max-out-NN with only one hidden layer and two neurons.


title: Variance Reduction for Policy-Gradient Methods via Empirical Variance Minimization

abstract: Policy-gradient methods in Reinforcement Learning(RL) are very universal and widely applied in practice but their performance suffers from the high variance of the gradient estimate. Several procedures were proposed to reduce it including actor-critic(AC) and advantage actor-critic(A2C) methods. Recently the approaches have got new perspective due to the introduction of Deep RL: both new control variates(CV) and new sub-sampling procedures became available in the setting of complex models like neural networks. The vital part of CV-based methods is the goal functional for the training of the CV, the most popular one is the least-squares criterion of A2C. Despite its practical success, the criterion is not the only one possible. In this paper we for the first time investigate the performance of the one called Empirical Variance(EV). We observe in the experiments that not only EV-criterion performs not worse than A2C but sometimes can be considerably better. Apart from that, we also prove some theoretical guarantees of the actual variance reduction under very general assumptions and show that A2C least-squares goal functional is an upper bound for EV goal. Our experiments indicate that in terms of variance reduction EV-based methods are much better than A2C and allow stronger variance reduction.


title: Multivariate Normality test for colored data

abstract: Performances of the Multivariate Kurtosis are investigated when applied to colored data, with or without Auto-Regressive pre-whitening, and with or without projection onto a lower-dimensional random subspace. Computer experiments demonstrate the importance of taking into account the possible color of the process in calculating the power of the normality test, in all the scenarios.


title: Efficient Decoder-free Object Detection with Transformers

abstract: Vision transformers (ViTs) are changing the landscape of object detection approaches. A natural usage of ViTs in detection is to replace the CNN-based backbone with a transformer-based backbone, which is straightforward and effective, with the price of bringing considerable computation burden for inference. More subtle usage is the DETR family, which eliminates the need for many hand-designed components in object detection but introduces a decoder demanding an extra-long time to converge. As a result, transformer-based object detection can not prevail in large-scale applications. To overcome these issues, we propose a novel decoder-free fully transformer-based (DFFT) object detector, achieving high efficiency in both training and inference stages, for the first time. We simplify objection detection into an encoder-only single-level anchor-based dense prediction problem by centering around two entry points: 1) Eliminate the training-inefficient decoder and leverage two strong encoders to preserve the accuracy of single-level feature map prediction; 2) Explore low-level semantic features for the detection task with limited computational resources. In particular, we design a novel lightweight detection-oriented transformer backbone that efficiently captures low-level features with rich semantics based on a well-conceived ablation study. Extensive experiments on the MS COCO benchmark demonstrate that DFFT_SMALL outperforms DETR by 2.5% AP with 28% computation cost reduction and more than $10$x fewer training epochs. Compared with the cutting-edge anchor-based detector RetinaNet, DFFT_SMALL obtains over 5.5% AP gain while cutting down 70% computation cost.


title: Mixed-Mode Chimera states in Pendula Networks

abstract: We report the emergence of peculiar chimera states in networks of identical pendula with global phase-lagged coupling. The states reported include both rotating and quiescent modes, i.e. with non-zero and zero average frequencies. This kind mixed-mode chimeras may be interpreted as images of bump states known in neuroscience in the context of modelling the working memory. We illustrate this striking phenomenon for a network of N=100 coupled pendula, followed by a detailed description of the minimal non-trivial case of N=3. Parameter regions for five characteristic types of the system behavior are identified consisting: two mixed-mode chimeras with one and two rotating pendula, classical weak chimera with all three pendula rotating, synchronous rotation and quiescent state. The network dynamics is multistable: up to four of the states can coexist in the system phase state as demonstrated through the basins of attraction. The analysis suggests that the robust mixed-mode chimera states can generically describe the complex dynamics of diverse pendula-like systems widespread in nature.


title: Never bet against Einstein

abstract: Among the geniuses of mankind, Einstein was probably one of those who made more erroneous claims, and often changed his opinion during the years on important scientific subjects. However, it is important to bear in mind that his mistakes were always due to a misinterpretation of the results obtained from his equations, since he was often biased by his own view of the universe. Einstein's equations have always given the correct answer and, till now, there is no evidence that they might not be correct.


title: PaleoJump: A database for abrupt transitions in past climates

abstract: Tipping points (TPs) in the Earth system have been studied with growing interest and concern in recent years due to the potential risk of anthropogenic forcing causing abrupt, and possibly irreversible, climate transitions. Paleoclimate records are essential for identifying TPs in Earth's past and for properly understanding the climate system's underlying nonlinearities and bifurcation mechanisms. Due to the variations in quality, resolution, and dating methods, it is crucial to select the records that give the best representation of past climates. Furthermore, as paleoclimate time series vary in their origin, time spans, and periodicities, an objective, automated methodology is crucial for identifying and comparing TPs. To reach this goal, we present here the PaleoJump database of carefully selected, high-resolution records originating in ice cores, marine sediments, speleothems, terrestrial records, and lake sediments. These records describe climate variability on centennial, millennial, or longer time scales and cover all the continents and ocean basins. We provide an overview of their spatial distribution and discuss the gaps in coverage. Our statistical methodology includes an augmented Kolmogorov-Smirnov test and Recurrence Quantification Analysis; it is applied here to selected records to automatically detect abrupt transitions therein and to investigate the presence of potential tipping elements. These transitions are shown in the PaleoJump database together with other essential information, including location, temporal scale and resolution, along with temporal plots. This database represents, therefore, a valuable resource for researchers investigating TPs in past climates.


title: Density Regression with Conditional Support Points

abstract: Density regression characterizes the conditional density of the response variable given the covariates, and provides much more information than the commonly used conditional mean or quantile regression. However, it is often computationally prohibitive in applications with massive data sets, especially when there are multiple covariates. In this paper, we develop a new data reduction approach for the density regression problem using conditional support points. After obtaining the representative data, we exploit the penalized likelihood method as the downstream estimation strategy. Based on the connections among the continuous ranked probability score, the energy distance, the $L_2$ discrepancy and the symmetrized Kullback-Leibler distance, we investigate the distributional convergence of the representative points and establish the rate of convergence of the density regression estimator. The usefulness of the methodology is illustrated by modeling the conditional distribution of power output given multivariate environmental factors using a large scale wind turbine data set. Supplementary materials for this article are available online.


title: Distributed Coordination of Charging Stations Considering Aggregate EV Power Flexibility

abstract: In recent years, electric vehicle (EV) charging stations have witnessed a rapid growth. However, effective management of charging stations is challenging due to individual EV owners' privacy concerns, competing interests of different stations, and the coupling distribution network constraints. To cope with this challenge, this paper proposes a two-stage scheme. In the first stage, the aggregate EV power flexibility region is derived by solving an optimization problem. We prove that any trajectory within the obtained region corresponds to at least one feasible EV dispatch strategy. By submitting this flexibility region instead of the detailed EV data to the charging station operator, EV owners' privacy can be preserved and the computational burden can be reduced. In the second stage, a distributed coordination mechanism with a clear physical interpretation is developed with consideration of AC power flow constraints. We prove that the proposed mechanism is guaranteed to converge to the centralized optimum. Case studies validate the theoretical results. Comprehensive performance comparisons are carried out to demonstrate the advantages of the proposed scheme.


title: A result on the $c_2$ invariant for powers of primes

abstract: The $c_2$ invariant is an arithmetic graph invariant related to quantum field theory. We give a relation modulo $p$ between the $c_2$ invariant at $p$ and the $c_2$ invariant at $p^s$ by proving a relation modulo $p$ between certain coefficients of powers of products of particularly nice polynomials. The relation at the level of the $c_2$ invariant provides evidence for a conjecture of Schnetz.


title: "hasSignification()": une nouvelle fonction de distance pour soutenir la d\'etection de donn\'ees personnelles

abstract: Today with Big Data and data lakes, we are faced of a mass of data that is very difficult to manage it manually. The protection of personal data in this context requires an automatic analysis for data discovery. Storing the names of attributes already analyzed in a knowledge base could optimize this automatic discovery. To have a better knowledge base, we should not store any attributes whose name does not make sense. In this article, to check if the name of an attribute has a meaning, we propose a solution that calculate the distances between this name and the words in a dictionary. Our studies on the distance functions like N-Gram, Jaro-Winkler and Levenshtein show limits to set an acceptance threshold for an attribute in the knowledge base. In order to overcome these limitations, our solution aims to strengthen the score calculation by using an exponential function based on the longest sequence. In addition, a double scan in dictionary is also proposed in order to process the attributes which have a compound name.


title: A dynamical theory for one-dimensional fermions with strong two-body losses: universal non-Hermitian Zeno physics and spin-charge separation

abstract: We study an interacting one-dimensional gas of spin-1/2 fermions with two-body losses. The dynamical phase diagram that characterises the approach to the stationary state displays a wide quantum-Zeno region, identified by a peculiar behaviour of the lowest eigenvalues of the associated non-Hermitian Hamiltonian. We characterise the universal dynamics of this Zeno regime using an approximation scheme based on a effective decoupling of charge and spin degrees of freedom, where the latter effectively evolve in a non-Markovian way according to a non-Hermitian Heisenberg Hamiltonian. We present detailed results for the time evolution from initial states with one particle per site with either incoherent and antiferromagnetic spin order, showing how the two non-equilibrium evolutions build up drastically different charge correlations.


title: Architectural patterns for handling runtime uncertainty of data-driven models in safety-critical perception

abstract: Data-driven models (DDM) based on machine learning and other AI techniques play an important role in the perception of increasingly autonomous systems. Due to the merely implicit definition of their behavior mainly based on the data used for training, DDM outputs are subject to uncertainty. This poses a challenge with respect to the realization of safety-critical perception tasks by means of DDMs. A promising approach to tackling this challenge is to estimate the uncertainty in the current situation during operation and adapt the system behavior accordingly. In previous work, we focused on runtime estimation of uncertainty and discussed approaches for handling uncertainty estimations. In this paper, we present additional architectural patterns for handling uncertainty. Furthermore, we evaluate the four patterns qualitatively and quantitatively with respect to safety and performance gains. For the quantitative evaluation, we consider a distance controller for vehicle platooning where performance gains are measured by considering how much the distance can be reduced in different operational situations. We conclude that the consideration of context information of the driving situation makes it possible to accept more or less uncertainty depending on the inherent risk of the situation, which results in performance gains.


title: Filtrations and torsion pairs in Abramovich Polishchuk's heart

abstract: We study some abelian subcategories and torsion pairs in Abramovich Polishchuk's heart. And we construct stability conditions on a full triangulated subcategory $\mathcal{D}^{\leq 1}_S$ in $D(X\times S)$, for an arbitrary smooth projective variety S. We also define a notion of $l$-th level stability, which is a generalization of the slope stability and the Gieseker stability. We show that for any object E in Abramovich Polishchuk's heart, there is a unique filtration whose factors are $l$-th level semistable, and the phase vectors are decreasing in a lexicographic order.


title: On the number of non-isomorphic choices on four elements

abstract: We use a combinatorial approach to compute the number of non-isomorphic choices on four elements that can be explained by models of bounded rationality.


title: Robust Reinforcement Learning with Distributional Risk-averse formulation

abstract: Robust Reinforcement Learning tries to make predictions more robust to changes in the dynamics or rewards of the system. This problem is particularly important when the dynamics and rewards of the environment are estimated from the data. In this paper, we approximate the Robust Reinforcement Learning constrained with a $\Phi$-divergence using an approximate Risk-Averse formulation. We show that the classical Reinforcement Learning formulation can be robustified using standard deviation penalization of the objective. Two algorithms based on Distributional Reinforcement Learning, one for discrete and one for continuous action spaces are proposed and tested in a classical Gym environment to demonstrate the robustness of the algorithms.


title: On neighborhoods of embedded complex tori

abstract: The goal of the article is to show that an n-dimensional complex torus embedded in a complex manifold of dimensional n+d, with a split tangent bundle, has neighborhood biholomorphic a neighborhood of the zero section in its normal bundle, provided the latter has (locally constant) Hermitian transition functions and satisfies a non-resonant Diophantine condition.


title: Role of particle local curvature in cellular wrapping

abstract: Cellular uptake through the lipid membranes plays an important role in adsorbing nutrients and fighting infection and can be used for drug delivery and nanomedicine developments. Endocytosis is one of the known pathways of the cellular uptake which associate with elastic deformation of the membrane wrapping around the foreign particle. The deformability of the membrane itself is strongly regulated by the presence of a cortical cytoskeleton placed underneath the membrane. It has been shown that size, shape, and orientation of the particles influence on their entry into the cell. Here, we study the role of particle local curvature in the cellular uptake by investigating the wrapping of an elastic membrane around a long cylindrical object with an elliptical cross section. The membrane itself is adhered to a substrate mimicking the cytoskeleton. Membrane wrapping proceeds differently whether the initial contact occurs at the particle's highly curved tip (prolate) or along its long side (oblate). We obtain a wrapping phase diagram as a function of the membrane-cytoskeleton and the membrane-target adhesion energy, which includes three distinct regimes of engulfment(unwrapped, partially wrapped, and fully wrapped), separated by two phase transitions. We also provide analytical expressions for the boundary between the different regimes which confirm that the transitions strongly depend on the orientation and aspect ratio of the particle.


title: Fully Automated Assessment of Cardiac Coverage in Cine Cardiovascular Magnetic Resonance Images using an Explainable Deep Visual Salient Region Detection Model

abstract: Cardiovascular magnetic resonance (CMR) imaging has become a modality with superior power for the diagnosis and prognosis of cardiovascular diseases. One of the essential basic quality controls of CMR images is to investigate the complete cardiac coverage, which is necessary for the volumetric and functional assessment. This study examines the full cardiac coverage using a 3D convolutional model and then reduces the number of false predictions using an innovative salient region detection model. Salient regions are extracted from the short-axis cine CMR stacks using a three-step proposed algorithm. Combining the 3D CNN baseline model with the proposed salient region detection model provides a cascade detector that can reduce the number of false negatives of the baseline model. The results obtained on the images of over 6,200 participants of the UK Biobank population cohort study show the superiority of the proposed model over the previous state-of-the-art studies. The dataset is the largest regarding the number of participants to control the cardiac coverage. The accuracy of the baseline model in identifying the presence/absence of basal/apical slices is 96.25\% and 94.51\%, respectively, which increases to 96.88\% and 95.72\% after improving using the proposed salient region detection model. Using the salient region detection model by forcing the baseline model to focus on the most informative areas of the images can help the model correct misclassified samples' predictions. The proposed fully automated model's performance indicates that this model can be used in image quality control in population cohort datasets and also real-time post-imaging quality assessments.


title: Effect of architecture disorder on the elastic response of two-dimensional lattice materials

abstract: We examine how disordering joint position influences the linear elastic behavior of lattice materials via numerical simulations in two-dimensional beam networks. Three distinct initial crystalline geometries are selected as representative of mechanically isotropic materials low connectivity, mechanically isotropic materials with high connectivity, and mechanically anisotropic materials with intermediate connectivity. Introducing disorder generates spatial fluctuations in the elasticity tensor at the local (joint) scale. Proper coarse-graining reveals a well-defined continuum-level scale elasticity tensor. Increasing disorder aids in making initially anisotropic materials more isotropic. The disorder impact on the material stiffness depends on the lattice connectivity: Increasing the disorder softens lattices with high connectivity and stiffens those with low connectivity, without modifying the scaling between elastic modulus and density (linear scaling for high connectivity and cubic scaling for low connectivity). Introducing disorder in lattices with intermediate fixed connectivity reveals both scaling: the linear scaling occurs for low density, the cubic one at high density, and the crossover density increases with disorder. Contrary to classical formulations, this work demonstrates that connectivity is not the sole parameter governing elastic modulus scaling. It offers a promising route to access novel mechanical properties in lattice materials via disordering the architectures.


title: Combining Image Space and q-Space PDEs for Lossless Compression of Diffusion MR Images

abstract: Diffusion MRI is a modern neuroimaging modality with a unique ability to acquire microstructural information by measuring water self-diffusion at the voxel level. However, it generates huge amounts of data, resulting from a large number of repeated 3D scans. Each volume samples a location in q-space, indicating the direction and strength of a diffusion sensitizing gradient during the measurement. This captures detailed information about the self-diffusion, and the tissue microstructure that restricts it. Lossless compression with GZIP is widely used to reduce the memory requirements. We introduce a novel lossless codec for diffusion MRI data. It reduces file sizes by more than 30% compared to GZIP, and also beats lossless codecs from the JPEG family. Our codec builds on recent work on lossless PDE-based compression of 3D medical images, but additionally exploits smoothness in q-space. We demonstrate that, compared to using only image space PDEs, q-space PDEs further improve compression rates. Moreover, implementing them with Finite Element Methods and a custom acceleration significantly reduces computational expense. Finally, we show that our codec clearly benefits from integrating subject motion correction, and slightly from optimizing the order in which the 3D volumes are coded.


title: On the Finite-Time Performance of the Knowledge Gradient Algorithm

abstract: The knowledge gradient (KG) algorithm is a popular and effective algorithm for the best arm identification (BAI) problem. Due to the complex calculation of KG, theoretical analysis of this algorithm is difficult, and existing results are mostly about the asymptotic performance of it, e.g., consistency, asymptotic sample allocation, etc. In this research, we present new theoretical results about the finite-time performance of the KG algorithm. Under independent and normally distributed rewards, we derive bounds for the sample allocation of the algorithm. With these bounds, existing asymptotic results become simple corollaries. Furthermore, we derive upper and lower bounds for the probability of error and simple regret of the algorithm, and show the performance of the algorithm for the multi-armed bandit (MAB) problem. These developments not only extend the existing analysis of the KG algorithm, but can also be used to analyze other improvement-based algorithms. Last, we use numerical experiments to compare the bounds we derive and the performance of the KG algorithm.


title: The lifetime of charged dust in the atmosphere

abstract: Windblown dust plays critical roles in numerous geophysical and biological systems, yet current models fail to explain the transport of coarse-mode particles (>5 $\mu$m) to great distances from their sources. For particles larger than a few microns, electrostatic effects have been invoked to account for longer-than-predicted atmospheric residence times. Although much effort has focused on elucidating the charging processes, comparatively little effort has been expended understanding the stability of charge on particles once electrified. Overall, electrostatic-driven transport requires that charge remain present on particles for days to weeks. Here, we present a set of experiments designed to explore the longevity of electrostatic charge on levitated airborne particles after a single charging event. Using an acoustic levitator, we measured the charge on particles of different material compositions suspended in atmospheric conditions for long periods of time. In dry environments, the total charge on particles decayed in over 1 week. The decay timescale decreased to days in humid environments. These results were independent of particle material and charge polarity. However, exposure to UV radiation could both increase or decrease the decay time depending on polarity. Our work suggests that the rate of charge decay on airborne particles is solely determined by ion capture from the air. Furthermore, using a one-dimensional sedimentation model, we predict that atmospheric dust of order 10 $\mu$m will experience the largest change in residence time due to electrostatic forces.


title: On the Gauss-Manin Connection and Real Singularities

abstract: We prove that to each real singularity $f: (\mathbb{R}^{n+1}, 0) \to (\mathbb{R}, 0)$ one can associate two systems of differential equations $\mathfrak{g}^{k\pm}_f$ which are pushforwards in the category of $\mathcal{D}$-modules over $\mathbb{R}^{\pm}$, of the sheaf of real analytic functions on the total space of the positive, respectively negative, Milnor fibration. We prove that for $k=0$ if $f$ is an isolated singularity then $\mathfrak{g}^{\pm}$ determines the the $n$-th homology groups of the positive, respectively negative, Milnor fibre. We then calculate $\mathfrak{g}^{+}$ for ordinary quadratic singularities and prove that under certain conditions on the choice of morsification, one recovers the top homology groups of the Milnor fibers of any isolated singularity $f$. As an application we construct a public-key encryption scheme based on morsification of singularities.


title: Four-field Hamiltonian fluid closures of the one-dimensional Vlasov-Poisson equation

abstract: We consider a reduced dynamics for the first four fluid moments of the onedimensional Vlasov-Poisson equation, namely, the fluid density, fluid velocity, pressure and heat flux. This dynamics depends on an equation of state to close the system. This equation of state (closure) connects the fifth order moment-related to the kurtosis in velocity of the Vlasov distribution-with the first four moments. By solving the Jacobi identity, we derive an equation of state which ensures that the resulting reduced fluid model is Hamiltonian. We show that this Hamiltonian closure allows symmetric homogeneous equilibria of the reduced fluid model to be stable.


title: Hanbury Brown-Twiss effect with electromagnetic scattered field generated by a collection of particles of $\mathcal{L}$ types

abstract: A theoretical framework in the spherical polar coordinate system is developed to systematically treat the correlation between intensity fluctuations (CIF) of electromagnetic light waves on scattering from a collection of particles of $\mathcal{L}$ types. Two $\mathcal{L}\times \mathcal{L}$ matrices called pair-potential matrix (PPM) and pair-structure matrix (PSM) are introduced to jointly formulate the normalized CIF of the scattered field for the first time. We build a closed-form relation that associates the normalized CIF with the PPM and the PSM as well as the spectral degree of polarization $\mathcal{P}$ of the incident field, showing that the normalized CIF is closely related to the trace of the product of the PSM and the transpose of the PPM, and its dependence on $\mathcal{P}$ is completely determined by the scattering polar angle and azimuth angle. For a special case where the spatial distributions of scattering potentials of particles of different types are similar and the same is true of their density distributions, the PPM and the PSM will reduce to two new matrices whose elements separately quantify the degree of angular correlation of the scattering potentials of particles and their density distributions, and the number of species of particles in this special case appears as a scaled factor to ensure the normalization of the CIF. The effects of the off-diagonal elements of the PPM and the PSM on the normalized CIF and its dependence on $\mathcal{P}$ are illustrated by two numerical examples.


title: Probing non-standard neutrino interactions with a light boson from next galactic and diffuse supernova neutrinos

abstract: Non-standard neutrino interactions with a massive boson can produce the bosons in the core of core-collapse supernovae (SNe). After the emission of the bosons from the SN core, their subsequent decays into neutrinos can modify the SN neutrino flux. We show future observations of neutrinos from a next galactic SN in Super-Kamiokande (SK) and Hyper-Kamiokande (HK) can probe flavor-universal non-standard neutrino couplings to a light boson, improving the previous limit from the SN 1987A neutrino burst by several orders of magnitude. We also discuss sensitivity of the flavor-universal non-standard neutrino interactions in future observations of diffuse neutrinos from all the past SNe, known as the diffuse supernova neutrino background (DSNB). According to our analysis, observations of the DSNB in HK, JUNO and DUNE experiments can probe such couplings by a factor of $\sim 2$ beyond the SN 1987A constraint. However, our result is also subject to a large uncertainty concerning the precise estimation of the DSNB.


title: Fast convergence of inertial dynamics with Hessian-driven damping under geometry assumptions

abstract: First-order optimization algorithms can be considered as a discretization of ordinary differential equations (ODEs) \cite{su2014differential}. In this perspective, studying the properties of the corresponding trajectories may lead to convergence results which can be transfered to the numerical scheme. In this paper we analyse the following ODE introduced by Attouch et al. in \cite{attouch2016fast}: \begin{equation*} \forall t\geqslant t_0,~\ddot{x}(t)+\frac{\alpha}{t}\dot{x}(t)+\beta H_F(x(t))\dot{x}(t)+\nabla F(x(t))=0,\end{equation*} where $\alpha>0$, $\beta>0$ and $H_F$ denotes the Hessian of $F$. This ODE can be derived to build numerical schemes which do not require $F$ to be twice differentiable as shown in \cite{attouch2020first,attouch2021convergence}. We provide strong convergence results on the error $F(x(t))-F^*$ and integrability properties on $\|\nabla F(x(t))\|$ under some geometry assumptions on $F$ such as quadratic growth around the set of minimizers. In particular, we show that the decay rate of the error for a strongly convex function is $O(t^{-\alpha-\varepsilon})$ for any $\varepsilon>0$. These results are briefly illustrated at the end of the paper.


title: On the explainable properties of 1-Lipschitz Neural Networks: An Optimal Transport Perspective

abstract: Input gradients have a pivotal role in a variety of applications, including adversarial attack algorithms for evaluating model robustness, explainable AI techniques for generating Saliency Maps, and counterfactual explanations.However, Saliency Maps generated by traditional neural networks are often noisy and provide limited insights. In this paper, we demonstrate that, on the contrary, the Saliency Maps of 1-Lipschitz neural networks, learned with the dual loss of an optimal transportation problem, exhibit desirable XAI properties:They are highly concentrated on the essential parts of the image with low noise, significantly outperforming state-of-the-art explanation approaches across various models and metrics. We also prove that these maps align unprecedentedly well with human explanations on ImageNet.To explain the particularly beneficial properties of the Saliency Map for such models, we prove this gradient encodes both the direction of the transportation plan and the direction towards the nearest adversarial attack. Following the gradient down to the decision boundary is no longer considered an adversarial attack, but rather a counterfactual explanation that explicitly transports the input from one class to another. Thus, Learning with such a loss jointly optimizes the classification objective and the alignment of the gradient, i.e. the Saliency Map, to the transportation plan direction.These networks were previously known to be certifiably robust by design, and we demonstrate that they scale well for large problems and models, and are tailored for explainability using a fast and straightforward method.


title: A new convergence proof for approximations of the Stefan problem

abstract: We consider the Stefan problem, firstly with regular data and secondly with irregular data. In both cases is given a proof for the convergence of an approximation obtained by regularising the problem. These proofs are based on weak formulations and on compactness results in some Sobolev spaces with negative exponents.


title: Finite time extinction for a diffusion equation with spatially inhomogeneous strong absorption

abstract: The phenomenon of finite time extinction of bounded and non-negative solutions to the diffusion equation with strong absorption $$\partial_t u-\Delta u^m+|x|^{\sigma}u^q=0, \qquad (t,x)\in(0,\infty)\times\mathbb{R}^N,$$ with $m\geq1$, $q\in(0,1)$ and $\sigma>0$, is addressed. Introducing the critical exponent $\sigma^* := 2(1-q)/(m-1)$ for $m>1$ and $\sigma_*=\infty$ for $m=1$, extinction in finite time is known to take place for $\sigma\in [0,\sigma^*)$ and an alternative proof is provided therein. When $m>1$ and $\sigma\ge \sigma^*$, the occurrence of finite time extinction is proved for a specific class of initial conditions, thereby supplementing results on non-extinction that are available in that range of $\sigma$ and showing their sharpness.


title: Model Predictive Control of a Tandem-Rotor Helicopter With a Nonuniformly Spaced Prediction Horizon

abstract: This letter considers model predictive control of a tandem-rotor helicopter. The error is formulated using the matrix Lie group $SE_2(3)$. A reference trajectory to a target is calculated using a quartic guidance law, leveraging the differentially flat properties of the system, and refined using a finite-horizon linear quadratic regulator. The nonlinear system is linearized about the reference trajectory enabling the formulation of a quadratic program with control input, attitude keep-in zone, and attitude error constraints. A non-uniformly spaced prediction horizon is leveraged to capture the multi-timescale dynamics while keeping the problem size tractable. Monte-Carlo simulations demonstrate robustness of the proposed control structure to initial conditions, model uncertainty, and environmental disturbances.


title: Monoidal Kleisli Bicategories and the Arithmetic Product of Coloured Symmetric Sequences

abstract: We extend the arithmetic product of species of structures and symmetric sequences studied by Maia and Mendez and by Dwyer and Hess to coloured symmetric sequences and show that it determines a normal oplax monoidal structure on the bicategory of coloured symmetric sequences. In order to do this, we establish general results on extending monoidal structures to Kleisli bicategories. Our approach uses monoidal double categories, which help us to attack the difficult problem of verifying the coherence conditions for a monoidal bicategory in an efficient way.


title: A recursive coloring function without $\Pi_3^0$ solution for Hindman theorem

abstract: We show that there exists a recursive coloring function $c$ such that any $\Pi^0_3$ set is not a solution to $c$ for Hindman's theorem. We also show that there exists a recursive coloring function $c$ such that any $\Delta^0_3$ set is not a solution to $c$ for Hindman's theorem restricted on sum of at most three numbers.


title: The local Universe in the era of large surveys -- II. Multi-wavelength characterisation of activity in nearby S0 galaxies

abstract: This is the second paper in a series using data from tens of thousands S0 galaxies of the local Universe ($z\lesssim 0.1$) retrieved from the NASA-Sloan Atlas. It builds on the outcomes of the previous work, which introduced a new classification scheme for these objects based on the principal component analysis (PCA) of their optical spectrum and its projections on to the first two eigenvectors or principal components (the PC1$\unicode{x2013}$PC2 diagram). We provide a comprehensive characterization of the activity of present-day S0s throughout both the broad-band PC1$\unicode{x2013}$PC2 spectral classifier and the conventional narrow-line BPT/WHAN ones, contrasting the different types of activity classes they define, and present an alternative diagram that exploits the concordance between WHAN and PCA demarcations. The analysis is extended to the mid-infrared, radio and X-ray wavelengths by crossmatching our core sample with data from the WISE, FIRST, XMM$\unicode{x2013}$Newton, and Chandra surveys. This has allowed us to carry out a thorough comparison of the most important activity diagnostics in the literature over different wavebands, discuss their similarities and differences, and explore the connections between them and with parameters related to star formation and black hole accretion. In particular, we find evidence that the bulk of nebular emission from radio and X-ray detected S0$\unicode{x2013}$Seyfert and LINER systems is not driven by star birth, while the dominant ionising radiation for a number of LINERs might come from post-AGB stars. These and other outcomes from the present work should be transferable to other morphologies.


title: The Stieltjes--Fekete problem and degenerate orthogonal polynomials

abstract: A famous result of Stieltjes relates the zeroes of the classical orthogonal polynomials with the configurations of points on the line that minimize a suitable energy. The energy has logarithmic interactions and an external field whose exponential is related to the weight of the classical orthogonal polynomials. The optimal configuration satisfies an algebraic set of equations: we call this set of algebraic equations the Stieltjes--Fekete problem or equivalently the Stieltjes--Bethe equations. In this work we consider the Stieltjes-Fekete problem when the derivative of the external field is an arbitrary rational complex function. We show that its solutions are in one-to-one correspondence with the zeroes of certain non-hermitean orthogonal polynomials that satisfy an excess of orthogonality conditions and are thus termed "degenerate". This generalizes the original result of Stieltjes.


title: Evaluating histopathology transfer learning with ChampKit

abstract: Histopathology remains the gold standard for diagnosis of various cancers. Recent advances in computer vision, specifically deep learning, have facilitated the analysis of histopathology images for various tasks, including immune cell detection and microsatellite instability classification. The state-of-the-art for each task often employs base architectures that have been pretrained for image classification on ImageNet. The standard approach to develop classifiers in histopathology tends to focus narrowly on optimizing models for a single task, not considering the aspects of modeling innovations that improve generalization across tasks. Here we present ChampKit (Comprehensive Histopathology Assessment of Model Predictions toolKit): an extensible, fully reproducible benchmarking toolkit that consists of a broad collection of patch-level image classification tasks across different cancers. ChampKit enables a way to systematically document the performance impact of proposed improvements in models and methodology. ChampKit source code and data are freely accessible at https://github.com/kaczmarj/champkit .


title: How are policy gradient methods affected by the limits of control?

abstract: We study stochastic policy gradient methods from the perspective of control-theoretic limitations. Our main result is that ill-conditioned linear systems in the sense of Doyle inevitably lead to noisy gradient estimates. We also give an example of a class of stable systems in which policy gradient methods suffer from the curse of dimensionality. Our results apply to both state feedback and partially observed systems.


title: Searching for Afterglow: Light Dark Matter Boosted by Supernova Neutrinos

abstract: A novel analysis is performed, incorporating time-of-flight (TOF) information to study the interactions of dark matter (DM) with standard model particles. After supernova (SN) explosions, DM with mass $m_\chi\lesssim\mathcal{O}({\rm MeV})$ in the halo can be boosted by SN neutrinos (SN$\nu$) to relativistic speed. The SN$\nu$ boosted DM (BDM) arrives on Earth with TOF which depends only on $m_\chi$ and is independent of the cross section. These BDMs can interact with detector targets in low-background experiments and manifest as afterglow events after the arrival of SN$\nu$. The characteristic TOF spectra of the BDM events can lead to large background suppression and unique determination of $m_\chi$. New cross section constraints on $\sqrt{\sigma_{\chi e} \sigma_{\chi\nu}}$ are derived from SN1987a in the Large Magellanic Cloud with data from the Kamiokande and Super-Kamiokande experiments. Potential sensitivities for the next galactic SN with Hyper-Kamiokande are projected. This analysis extends the existing bounds on $\sqrt{\sigma_{\chi e}\sigma_{\chi \nu}}$ over a broad range of $r_\chi=\sigma_{\chi \nu}/\sigma_{\chi e}$. In particular, the improvement is by 1-3 orders of magnitude for $m_\chi<\mathcal{O}(100\,{\rm keV})$ for $\sigma_{\chi e}\sim\sigma_{\chi \nu}$. Prospects of exploiting TOF information in other astrophysical systems to probe exotic physics with other DM candidates are discussed.


title: Rough paths and symmetric-Stratonovich integrals driven by singular covariance gaussian processes

abstract: We examine the relation between a stochastic version of the rough path integral with the symmetric-Stratonovich integral in the sense of regularization. Under mild regularity conditions in the sense of Malliavin calculus, we establish equality between stochastic rough path and symmetric-Stratonovich integrals driven by a class of Gaussian processes. As a by-product, we show that solutions of multi-dimensional rough differential equations driven by a large class of Gaussian rough paths they are actually solutions to Stratonovich stochastic differential equations. We obtain almost sure convergence rates of the first-order Stratonovich scheme to rough paths integrals in the sense of Gubinelli. In case the time-increment of the Malliavin derivative of the integrands is regular enough, the rates are essentially sharp. The framework applies to a large class of Gaussian processes whose the second-order derivative of the covariance function is a sigma-finite non-positive measure on ${\mathbb R}^2$ + off diagonal.


title: EXPTIME-hardness of higher-dimensional Minkowski spacetime

abstract: We prove the EXPTIME-hardness of the validity problem for the basic temporal logic on Minkowski spacetime with more than one space dimension. We prove this result for both the lightspeed-or-slower and the slower-than-light accessibility relations (and for both the irreflexive and the reflexive versions of these relations). As an auxiliary result, we prove the EXPTIME-hardness of validity on any frame for which there exists an embedding of the infinite complete binary tree satisfying certain conditions. The proof is by a reduction from the two-player corridor-tiling game.


title: Optomechanically induced transparency and directional amplification in a non-Hermitian optomechanical lattice

abstract: Cavity optomechanics is important in both quantum information processing and basic physics research. In this paper, we propose an optomechanical lattice which manifests non-Hermitian physics . We first use the non-Bloch band theory to investigate the energy spectrum and transmission property of an optomechanical lattice. The generalized Brillouin zone of the system is calculated with the help of the resultant. And the periodical boundary condition (PBC) and open boundary condition energy spectrum are given, subsequently. By introducing probe laser on different sites we observed the directional amplification of the system. The direction of the amplification is analyzed combined with the non-Hermitian skin effect. The frequency that supports the amplification is analyzed by considering the PBC energy spectrum. By introducing probe laser on one site we investigate the onsite transmission property. Optomechanically induced transparency (OMIT) can be achieved in our system. By varying the parameters and size of the system, the OMIT peak can be effectively modulated or even turned into optomechanically induced amplification . Our system shows its potential as the function of a single-way signal filter. And our model can be extended to other non-Hermitian Bosonic model which may possess topological features and bipolar non-Hermitian skin effect.


title: Natural Language Sentence Generation from API Specifications

abstract: APIs are everywhere; they provide access to automation solutions that could help businesses automate some of their tasks. Unfortunately, they may not be accessible to the business users who need them but are not equipped with the necessary technical skills to leverage them. Wrapping these APIs with chatbot capabilities is one solution to make these automation solutions interactive. In this work, we propose a system to generate sentences to train intent recognition models, a crucial component within chatbots to understand natural language utterances from users. Evaluation of our approach based on deep learning models showed promising and inspiring results, and the human-in-the-loop interaction will provide further improvement on the system.


title: k-Deformed Fourier Transform

abstract: We present a new formulation of Fourier transform in the picture of the $\kappa$-algebra derived in the framework of the $\kappa$-generalized statistical mechanics. The $\kappa$-Fourier transform is obtained from a $\kappa$-Fourier series recently introduced by us [2013 Entropy {\bf15} 624]. The kernel of this transform, that reduces to the usual exponential phase in the $\kappa\to0$ limit, is composed by a $\kappa$-deformed phase and a damping factor that gives a wavelet-like behavior. We show that the $\kappa$-Fourier transform is isomorph to the standard Fourier transform through a changing of time and frequency variables. Nevertheless, the new formalism is useful to study, according to Fourier analysis, those functions defined in the realm of the $\kappa$-algebra. As a relevant application, we discuss the central limit theorem for the $\kappa$-sum of $n$-iterate statistically independent random variables.


title: A Novel RIS-Aided EMF Exposure Aware Approach using an Angularly Equalized Virtual Propagation Channel

abstract: Massive Multiple-Input Multiple-Output systems with beamforming are key components of the 5th and the future 6th generation of networks. However, in some cases, where the BS serves the same user for a long period, and in some propagation conditions, such systems reduce their transmit power to avoid creating unwanted regions of electromagnetic field exposure exceeding the regulatory threshold, beyond the circle around the BS that limits the distance between people and the BS antenna. Such power reduction strongly degrades the received power at the target user. Recently, exposition aware beamforming schemes aided by self-tuning reconfigurable intelligent surfaces derived from maximum ratio transmission beamforming, have been proposed: truncated beamforming. However, such scheme is highly complex. In this paper, we propose a novel and low complexity reconfigurable intelligent surface aided beamforming scheme called Equalized beamforming, which applies maximum ratio transmission to an angularly equalized virtual propagation channel. Our simulations show that our proposed scheme outperforms the reduced beamforming scheme, whilst complying with the exposition regulation.


title: An ion-independent phenomenological relative biological effectiveness (RBE) model for proton therapy

abstract: Background: A relative biological effectiveness (RBE) of 1.1 is used for proton therapy though clinical evidence of varying RBE was raised. Clinical studies on RBE variability have been conducted for decades for carbon radiation, which could advance the understanding of the clinical proton RBE given an ion-independent RBE model. In this work, such a model, linear and simple, using the beam quantity Q = Z^2/E (Z = ion charge, E = kinetic energy per nucleon) was tested and compared to the commonly used, proton-specific and linear energy transfer (LET) based Wedenberg RBE model. Material and methods: The Wedenberg and Q models, both predicting RBEmax and RBEmin (i.e., RBE at vanishing and very high dose, respectively), are compared in terms of ion-dependence and prediction power. An experimental in-vitro data ensemble covering 115 publications for various ions was used as dataset. Results: The model parameter of the Q model was observed to be similar for different ions (in contrast to LET). The Q model was trained without any prior knowledge of proton data. For proton RBE, the differences between experimental data and corresponding predictions of the Wedenberg or the Q model were highly comparable. Conclusions: A simple linear RBE model using Q instead of LET was proposed and tested to be able to predict proton RBE using model parameter trained based on only RBE data of other particles in a clinical proton energy range for a large in-vitro dataset. Adding (pre)clinical knowledge from carbon ion therapy may, therefore, reduce the dominating biological uncertainty in proton RBE modelling. This would translate in reduced RBE related uncertainty in proton therapy treatment planning.


title: On Provably Robust Meta-Bayesian Optimization

abstract: Bayesian optimization (BO) has become popular for sequential optimization of black-box functions. When BO is used to optimize a target function, we often have access to previous evaluations of potentially related functions. This begs the question as to whether we can leverage these previous experiences to accelerate the current BO task through meta-learning (meta-BO), while ensuring robustness against potentially harmful dissimilar tasks that could sabotage the convergence of BO. This paper introduces two scalable and provably robust meta-BO algorithms: robust meta-Gaussian process-upper confidence bound (RM-GP-UCB) and RM-GP-Thompson sampling (RM-GP-TS). We prove that both algorithms are asymptotically no-regret even when some or all previous tasks are dissimilar to the current task, and show that RM-GP-UCB enjoys a better theoretical robustness than RM-GP-TS. We also exploit the theoretical guarantees to optimize the weights assigned to individual previous tasks through regret minimization via online learning, which diminishes the impact of dissimilar tasks and hence further enhances the robustness. Empirical evaluations show that (a) RM-GP-UCB performs effectively and consistently across various applications, and (b) RM-GP-TS, despite being less robust than RM-GP-UCB both in theory and in practice, performs competitively in some scenarios with less dissimilar tasks and is more computationally efficient.


title: A Data-Driven Simulation of the New York State Foster Care System

abstract: We introduce an analytic pipeline to model and simulate youth trajectories through the New York state foster care system. Our goal in doing so is to forecast how proposed interventions may impact the foster care system's ability to achieve it's stated goals \emph{before these interventions are actually implemented and impact the lives of thousands of youth}. Here, we focus on two specific stated goals of the system: racial equity, and, as codified most recently by the 2018 Family First Prevention Services Act (FFPSA), a focus on keeping all youth out of foster care. We also focus on one specific potential intervention -- a predictive model, proposed in prior work and implemented elsewhere in the U.S., which aims to determine whether or not a youth is in need of care. We use our method to explore how the implementation of this predictive model in New York would impact racial equity and the number of youth in care. While our findings, as in any simulation model, ultimately rely on modeling assumptions, we find evidence that the model would not necessarily achieve either goal. Primarily, then, we aim to further promote the use of data-driven simulation to help understand the ramifications of algorithmic interventions in public systems.


title: Does Open Access Really Increase Impact? A Large-Scale Randomized Analysis

abstract: The Open Access Citation Advantage (OACA) has been a major topic of discussion in the literature over the past twenty years. In this paper, we propose a method to constitute a control group to isolate the OACA effect. Thus, we compared citation impact (MNCS) of 2,458,378 publications in fully OA journals to that (weighted MNCS) of a control group of non-OA publications (\#10,310,842). Similarly, we did the same exercise for OA publications in hybrid journals (\#1,024,430) and their control group (\#11,533,001), over the period 2010-2020. The results showed that there is no OACA for publications in fully OA journals, and that there is rather a disadvantage. Conversely, the OACA seems to be a reality in hybrid journals, suggesting that a better accessibility in this context tends to improve the visibility of publications. The lack of OACA for publications in fully OA journals is to be expected, as a great proportion of OA journals are newly created and less attractive to high-impact senior researchers.Another striking result of this paper is the fall of the OACA from 2016. The citation advantage fell from 70% to 9% between 2016 and 2020 (for publications in hybrid journals). We wonder if this fall is linked to the increase in the notoriety of pirate sites (eg Sci-Hub) from 2016. In other words, the democratization of pirate sites instantly cancels the positive effect of OA publication insofar as the question of access to scientific content no longer arises.


title: A Lagrangian free-stream boundary condition for weakly compressible smoothed particle hydrodynamics

abstract: This paper proposes the first free-stream boundary condition in a purely Lagrangian framework for weakly-compressible smoothed particle hydrodynamics (WCSPH). The boundary condition is implemented based on several numerical techniques, including a spatio-temporal method and an ensuring technique for surface particle identification, far-field corrections of density and velocity near the free-stream boundary, and improved in- and outlet conditions which guarantee smooth inflow profile and free outflow particles without the issue of the truncated kernel. The accuracy and effectiveness of the method have been validated by the reference solutions of flow past a flat plate and circular cylinder, respectively. The method has also been tested with more complex flows to demonstrate its applicability, such as fluid-structure interaction (FSI) and 3D flow past a sphere.


title: Comparing the hardness of MAX 2-SAT problem instances for quantum and classical algorithms

abstract: An algorithm for a particular problem may find some instances of the problem easier and others harder to solve, even for a fixed input size. We numerically analyse the relative hardness of MAX 2-SAT problem instances for various continuous-time quantum algorithms and a comparable classical algorithm. This has two motivations: to investigate whether small-sized problem instances, which are commonly used in numerical simulations of quantum algorithms for benchmarking purposes, are a good representation of larger instances in terms of their hardness to solve, and to determine the applicability of continuous-time quantum algorithms in a portfolio approach, where we take advantage of the variation in the hardness of instances between different algorithms by running them in parallel. We find that, while there are correlations in instance hardness between all of the algorithms considered, they appear weak enough that a portfolio approach would likely be desirable in practice. Our results also show a widening range of hardness of randomly generated instances as the problem size is increased, which demonstrates both the difference in the distribution of hardness at small sizes and the value of a portfolio approach that can reduce the number of extremely hard instances. We identify specific weaknesses of these quantum algorithms that can be overcome with a portfolio approach, such their inability to efficiently solve satisfiable instances (which is easy classically).


title: Intrinsically projectively linked graphs

abstract: A graph is intrinsically projectively linked (IPL) if its every embedding in projective space contains a nonsplit link. Some minor-minimal IPL graphs have been found previously. We determine that no minor-minimal IPL graphs on 16 edges exists and identify new minor-minimal IPL graphs by applying $\Delta-Y$ exchanges to $K_{7}-2e$. We prove that for a nonouter-projective-planar graph $G$, $G+\bar{K}_{2}$ is IPL and describe the necessary and sufficient conditions on a projective planar graph $G$ such that $G+\bar{K}_{2}$ is IPL. Lastly, we deduce conditions for $f(G + \bar{K_{2}})$ to have no nonsplit link, where $G$ is projective planar, $\bar{K_{2}} = \{w_{0},w_{1}\}$, and $f(G + \bar{K_{2}})$ is the embedding onto $\mathbb{R}P^{3}$ with $f(G)$ in $z=0$, $w_{0}$ above $z=0$, and $w_{1}$ below $z=0$ such that every edge connecting ${w_{0},w_{1}}$ to $G$ avoids the boundary of the 3-ball, whose antipodal points are identified to obtain projective space.


title: Temporal Multimodal Multivariate Learning

abstract: We introduce temporal multimodal multivariate learning, a new family of decision making models that can indirectly learn and transfer online information from simultaneous observations of a probability distribution with more than one peak or more than one outcome variable from one time stage to another. We approximate the posterior by sequentially removing additional uncertainties across different variables and time, based on data-physics driven correlation, to address a broader class of challenging time-dependent decision-making problems under uncertainty. Extensive experiments on real-world datasets ( i.e., urban traffic data and hurricane ensemble forecasting data) demonstrate the superior performance of the proposed targeted decision-making over the state-of-the-art baseline prediction methods across various settings.


title: Creating a Secure Underlay for the Internet

abstract: Adversaries can exploit inter-domain routing vulnerabilities to intercept communication and compromise the security of critical Internet applications. Meanwhile the deployment of secure routing solutions such as Border Gateway Protocol Security (BGPsec) and Scalability, Control and Isolation On Next-generation networks (SCION) are still limited. How can we leverage emerging secure routing backbones and extend their security properties to the broader Internet? We design and deploy an architecture to bootstrap secure routing. Our key insight is to abstract the secure routing backbone as a virtual Autonomous System (AS), called Secure Backbone AS (SBAS). While SBAS appears as one AS to the Internet, it is a federated network where routes are exchanged between participants using a secure backbone. SBAS makes BGP announcements for its customers' IP prefixes at multiple locations (referred to as Points of Presence or PoPs) allowing traffic from non-participating hosts to be routed to a nearby SBAS PoP (where it is then routed over the secure backbone to the true prefix owner). In this manner, we are the first to integrate a federated secure non-BGP routing backbone with the BGP-speaking Internet. We present a real-world deployment of our architecture that uses SCIONLab to emulate the secure backbone and the PEERING framework to make BGP announcements to the Internet. A combination of real-world attacks and Internet-scale simulations shows that SBAS substantially reduces the threat of routing attacks. Finally, we survey network operators to better understand optimal governance and incentive models.


title: Creating and Operating Areas With Reduced Electromagnetic Field Exposure Thanks to Reconfigurable Intelligent Surfaces

abstract: Mobile network operators must provide a target quality of service within a target coverage area. Each generation of networks from the 2nd to the 5th has reached higher quality-of-service (QoS) targets and coverage area sizes. However, the deployment of new networks is sometimes challenged by Electromagnetic Field (EMF) exposure constraints. In this paper, to take into account these constraints, we assess the novel and recent concept of Reduced EMF Exposure Area. Such an area would be created and operated by a mobile network operator upon the request of its customers. In such an area, customers keep enjoying high data rate internet access while getting a reduced EMF exposure. To ensure this EMF exposure reduction, we propose to deploy Reconfigurable Intelligent Surfaces (RIS) connected to the mobile network and exploit a joint RIS-M-MIMO uplink beamforming scheme. We use our ray-tracing-based simulation tool to visualize and characterize the Reduced EMF Exposure Area in a challenging environment in terms of propagation. Our simulations show that an operator can create and operate such an area under the condition that it carefully places the RIS in the environment.


title: Combinatorial Derived Matroids

abstract: Let $M$ be an arbitrary matroid with circuits $\mathcal{C}(M)$. We propose a definition of a derived matroid $\delta M$ that has as its ground set $\mathcal{C}(M)$. Unlike previous attempts of such a definition, our definition applies to arbitrary matroids, and is completely combinatorial. We prove that the rank of $\delta M$ is bounded from above by $|M|-r(M)$, that it is connected if and only if $M$ is connected. We compute examples including the derived matroids of uniform matroids, the V\'amos matroid and the graphical matroid $M(K_4)$. We formulate conjectures relating our construction to previous definitions of derived matroids.


title: An Accurate HDDL Domain Learning Algorithm from Partial and Noisy Observations

abstract: The Hierarchical Task Network ({\sf HTN}) formalism is very expressive and used to express a wide variety of planning problems. In contrast to the classical {\sf STRIPS} formalism in which only the action model needs to be specified, the {\sf HTN} formalism requires to specify, in addition, the tasks of the problem and their decomposition into subtasks, called {\sf HTN} methods. For this reason, hand-encoding {\sf HTN} problems is considered more difficult and more error-prone by experts than classical planning problem. To tackle this problem, we propose a new approach (HierAMLSI) based on grammar induction to acquire {\sf HTN} planning domain knowledge, by learning action models and {\sf HTN} methods with their preconditions. Unlike other approaches, HierAMLSI is able to learn both actions and methods with noisy and partial inputs observation with a high level or accuracy.


title: Charge instabilities of the extended attractive Hubbard Model on the cubic lattice

abstract: The paramagnetic phase of the extended attractive Hubbard model on the cubic lattice is studied within the spin rotation invariant Kotliar-Ruckenstein slave-boson representation at zero temperature. It is obtained that the quasiparticle residue of the Fermi liquid phase vanishes for all densities at an interaction strength slightly smaller than $U_c$ that signals the Brinkman-Rice transition, and that it weakly depends on density. While for vanishing non-local interaction parameters homogeneous static charge instabilities are found in a rather narrow window centered around quarter filling and $U \simeq 0.8\ U_c$, increasing them to $V = -0.2\ U$ results into a severe narrowing of this window. On the contrary, when all interaction parameters are attractive, for example for $V = 0.2\ U$, a large parameter range in which homogeneous static charge instabilities is found. Yet, this systematically happens inside the Fermi liquid phase.


title: Influence of impurities on electronic structure in cuprate superconductors

abstract: The impurity is inherently manifest in cuprate superconductors, as cation substitution or intercalation is necessary for the introduction of charge carriers, and its influence on the electronic state is at the heart of a great debate in physics. Here based on the microscopic octet scattering model, the influence of the impurity scattering on the electronic structure of cuprate superconductors is investigated in terms of the self-consistent T-matrix approach. The impurity scattering self-energy is evaluated firstly in the Fermi-arc-tip approximation of the quasiparticle excitations and scattering processes, and the obtained results show that the decisive role played by the impurity scattering self-energy in the particle-hole channel is the further renormalization of the quasiparticle band structure with a reduced quasiparticle lifetime, while the impurity scattering self-energy in the particle-particle channel induces a strong deviation from the d-wave behaviour of the superconducting gap, leading to the existence of a finite gap over the entire electron Fermi surface. Moreover, these impurity scattering self-energies are employed to study the exotic features of the line-shape in the quasiparticle excitation spectrum and the autocorrelation of the quasiparticle excitation spectra, and the obtained results are then compared with the corresponding experimental data. The theory therefore also indicates that the unconventional features of the electronic structure in cuprate superconductors is generated by both the strong electron correlation and impurity scattering.


title: Neural interval-censored Cox regression with feature selection

abstract: The classical Cox model emerged in 1972 promoting breakthroughs in how patient prognosis is quantified using time-to-event analysis in biomedicine. One of the most useful characteristics of the model for practitioners is the interpretability of the variables in the analysis. However, this comes at the price of introducing strong assumptions concerning the functional form of the regression model. To break this gap, this paper aims to exploit the explainability advantages of the classical Cox model in the setting of interval-censoring using a new Lasso neural network that simultaneously selects the most relevant variables while quantifying non-linear relations between predictors and survival times. The gain of the new method is illustrated empirically in an extensive simulation study with examples that involve linear and non-linear ground dependencies. We also demonstrate the performance of our strategy in the analysis of physiological, clinical and accelerometer data from the NHANES 2003-2006 waves to predict the effect of physical activity on the survival of patients. Our method outperforms the prior results in the literature that use the traditional Cox model.


title: Space-efficient Quantization Method for Reversible Markov Chains

abstract: In a seminal paper, Szegedy showed how to construct a quantum walk $W(P)$ for any reversible Markov chain $P$ such that its eigenvector with eigenphase $0$ is a quantum sample of the limiting distribution of the random walk and its eigenphase gap is quadratically larger than the spectral gap of $P$. The standard construction of Szegedy's quantum walk requires an ancilla register of Hilbert-space dimension equal to the size of the state space of the Markov chain. We show that it is possible to avoid this doubling of state space for certain Markov chains that employ a symmetric proposal probability and a subsequent accept/reject probability to sample from the Gibbs distribution. For such Markov chains, we give a quantization method which requires an ancilla register of dimension equal to only the number of different energy values, which is often significantly smaller than the size of the state space. To accomplish this, we develop a technique for block encoding Hadamard products of matrices which may be of wider interest.


title: Objectivity and Rigor in Classical Italian Algebraic Geometry

abstract: The classification of algebraic surfaces by the Italian School of algebraic geometry is universally recognized as a breakthrough in 20th-century mathematics. The methods by which it was achieved do not, however, meet the modern standard of rigor and therefore appear dubious from a contemporary viewpoint. In this article, we offer a glimpse into the mathematical practice of the three leading exponents of the Italian School of algebraic geometry: Castelnuovo, Enriques, and Severi. We then bring into focus their distinctive conception of rigor and intuition. Unlike what is often assumed today, from their perspective, rigor is neither opposed to intuition nor understood as a unitary phenomenon - Enriques distinguishes between small-scale rigor and large-scale rigor and Severi between formal rigor and substantial rigor. Finally, we turn to the notion of mathematical objectivity. We draw from our case study in order to advance a multi-dimensional analysis of objectivity. Specifically, we suggest that various types of rigor may be associated with different conceptions of objectivity: namely objectivity as faithfulness to facts and objectivity as intersubjectivity.


title: CERT: Continual Pre-Training on Sketches for Library-Oriented Code Generation

abstract: Code generation is a longstanding challenge, aiming to generate a code snippet based on a natural language description. Usually, expensive text-code paired data is essential for training a code generation model. Recently, thanks to the success of pre-training techniques, large language models are trained on large-scale unlabelled code corpora and perform well in code generation. In this paper, we investigate how to leverage an unlabelled code corpus to train a model for library-oriented code generation. Since it is a common practice for programmers to reuse third-party libraries, in which case the text-code paired data are harder to obtain due to the huge number of libraries. We observe that library-oriented code snippets are more likely to share similar code sketches. Hence, we present CERT with two steps: a sketcher generates the sketch, then a generator fills the details in the sketch. Both the sketcher and the generator are continually pre-trained upon a base model using unlabelled data. Furthermore, we craft two benchmarks named PandasEval and NumpyEval to evaluate library-oriented code generation. Experimental results demonstrate the impressive performance of CERT. For example, it surpasses the base model by an absolute 15.67% improvement in terms of pass@1 on PandasEval. Our work is available at https://github.com/microsoft/PyCodeGPT.


title: Symplectic leaves of Calogero-Moser spaces of type $G(\ell,1,n)$

abstract: We study symplectic leaves of Calogero-Moser spaces of type $G(\ell,1,n)$. We prove that the normalization of the closure of each symplectic leaf is isomorphic to some Calogero-Moser space. We also give a nice combinatorial parameterization of the symplectic leaves.


title: RDU: A Region-based Approach to Form-style Document Understanding

abstract: Key Information Extraction (KIE) is aimed at extracting structured information (e.g. key-value pairs) from form-style documents (e.g. invoices), which makes an important step towards intelligent document understanding. Previous approaches generally tackle KIE by sequence tagging, which faces difficulty to process non-flatten sequences, especially for table-text mixed documents. These approaches also suffer from the trouble of pre-defining a fixed set of labels for each type of documents, as well as the label imbalance issue. In this work, we assume Optical Character Recognition (OCR) has been applied to input documents, and reformulate the KIE task as a region prediction problem in the two-dimensional (2D) space given a target field. Following this new setup, we develop a new KIE model named Region-based Document Understanding (RDU) that takes as input the text content and corresponding coordinates of a document, and tries to predict the result by localizing a bounding-box-like region. Our RDU first applies a layout-aware BERT equipped with a soft layout attention masking and bias mechanism to incorporate layout information into the representations. Then, a list of candidate regions is generated from the representations via a Region Proposal Module inspired by computer vision models widely applied for object detection. Finally, a Region Categorization Module and a Region Selection Module are adopted to judge whether a proposed region is valid and select the one with the largest probability from all proposed regions respectively. Experiments on four types of form-style documents show that our proposed method can achieve impressive results. In addition, our RDU model can be trained with different document types seamlessly, which is especially helpful over low-resource documents.


title: Homomorphisms on algebras of analytic functions on non-symmetrically regular spaces

abstract: We study homomorphisms on the algebra of analytic functions of bounded type on a Banach space. When the domain space lacks symmetric regularity, we show that in every fiber of the spectrum there are evaluations (in higher duals) which do not coincide with evaluations in the second dual. We also consider the commutativity of convolutions between evaluations. We show that in some Banach spaces $X$ (for example, $X=\ell_1$) the only evaluations that commute with every other evaluation in $X''$ are those in $X$. Finally, we establish conditions ensuring the symmetry of the canonical extension of a symmetric multilinear operator (on a non-symmetrically regular space) and present some applications.


title: A new algorithm for structural restrictions in Bayesian vector autoregressions

abstract: A comprehensive methodology for inference in vector autoregressions (VARs) using sign and other structural restrictions is developed. The reduced-form VAR disturbances are driven by a few common factors and structural identification restrictions can be incorporated in their loadings in the form of parametric restrictions. A Gibbs sampler is derived that allows for reduced-form parameters and structural restrictions to be sampled efficiently in one step. A key benefit of the proposed approach is that it allows for treating parameter estimation and structural inference as a joint problem. An additional benefit is that the methodology can scale to large VARs with multiple shocks, and it can be extended to accommodate non-linearities, asymmetries, and numerous other interesting empirical features. The excellent properties of the new algorithm for inference are explored using synthetic data experiments, and by revisiting the role of financial factors in economic fluctuations using identification based on sign restrictions.


title: Possible Coexistence of Antiferromagnetic and Ferromagnetic Spin Fluctuations in the Spin-triplet Superconductor UTe2 Revealed by 125Te NMR under Pressure

abstract: A spin-triplet superconducting state mediated by ferromagnetic (FM) spin fluctuations has been suggested to occur in the newly discovered heavy-fermion superconductor UTe$_2$. However, the recent neutron scattering measurements revealed the presence of antiferromagnetic (AFM) spin fluctuations in UTe$_2$. Here, we report the $^{125}$Te nuclear magnetic resonance (NMR) studies of a single-crystal UTe$_2$, suggesting the coexistence of FM and AFM spin fluctuations in UTe$_2$. Owing to the two different Te sites in the compound, we conclude that the FM spin fluctuations are dominant within ladders and the AFM spin fluctuations originate from the inter-ladder magnetic coupling. Although AFM spin fluctuations exist in the system, the FM spin fluctuations in the ladders may play an important role in the appearance of the spin-triplet superconducting state of UTe$_2$.


title: Comparative Study of Magnetic Moment of Leptons in Hot and Dense Media

abstract: We study the magnetic moment of leptons in extremely hot universe and superdense media of stars at high temperatures. Anomalous magnetic moment of charged leptons is inversely proportional to its mass, whereas, the induced dipole moment of neutral leptons is directly proportional to their mass. Neutral massive point particles exert nonzero magnetic moment as a higher order effect, which is smaller than the anomalous magnetic moment of a charged particle of their same flavor partner. All leptons acquire some extra mass due to their interaction with the medium and affect the magnetic moment accordingly. We compare the contribution to the magnetic moment of various leptons due to their temperature and chemical potential dependent masses. It is shown that the magnetic moment contributions are non-ignorable for lighter leptons and heavy neutrinos. These calculations are very important to study particle propagation in the early universe and in superdense stellar media.


title: The cortical V1 transform as a heterogeneous Poisson problem

abstract: Receptive profiles of V1 cortical cells are very heterogeneous and act by differentiating the stimulus image as operators changing from point to point. A lightness and color constancy image can be reconstructed as the solution of the associated inverse problem, that is a Poisson equation with heterogeneous differential operators. At the neural level the weights of short range connectivity constitute the fundamental solution of the Poisson problem adapted point by point. A first demonstration of convergence of the result towards homogeneous reconstructions is proposed by means of homogenisation techniques.


title: Model reduction for second-order systems with inhomogeneous initial conditions

abstract: In this paper, we consider the problem of finding surrogate models for large-scale second-order linear time-invariant systems with inhomogeneous initial conditions. For this class of systems, the superposition principle allows us to decompose the system behavior into three independent components. The first behavior corresponds to the transfer between the input and output having zero initial conditions. In contrast, the other two correspond to the transfer between the initial position and the initial velocity conditions having zero input, respectively. Based on this superposition of systems, our goal is to propose model reduction schemes allowing to preserve the second-order structure in the surrogate models. To this aim, we introduce tailored second-order Gramians for each system component and compute them numerically, solving Lyapunov equations. As a consequence, two methodologies are proposed. The first one consists in reducing each of the components independently using a suitable balanced truncation procedure. The sum of these reduced systems provides an approximation of the original system. This methodology allows flexibility on the order of the reduced-order model. The second proposed methodology consists in extracting the dominant subspaces from the sum of Gramians to build the projection matrices leading to a surrogate model. Additionally, we discuss error bounds for the overall output approximation. Finally, the proposed methods are illustrated by means of benchmark problems.


title: On the Message Passing Efficiency of Polar and Low-Density Parity-Check Decoders

abstract: This study focuses on the efficiency of message-passing-based decoding algorithms for polar and low-density parity-check (LDPC) codes. Both successive cancellation (SC) and belief propagation (BP) decoding algorithms are studied {in} the message-passing framework. Counter-intuitively, SC decoding demonstrates the highest decoding efficiency, although it was considered a weak decoder {in terms of} error-correction performance. We analyze the complexity-performance tradeoff to dynamically track the decoding efficiency, where the complexity is measured by the number of messages passed (NMP), and the performance is measured by the statistical distance to the maximum a posteriori (MAP) estimate. This study offers a new insight into the contribution of each message passed in decoding, and compares various decoding algorithms on a message-by-message level. The analysis corroborates recent results on terabits-per-second polar SC decoders, and might shed light on better scheduling strategies.


title: Anti-Ramsey theory problems, lattice point counts on polytopes, and Hodge structures on the cohomology of toric varieties

abstract: We find families of graphs $G$ and subgraphs $H$ of $G$ such that the number of edge colorings of $G$ avoiding a monochromatic coloring of $H$ is determined by lattice point counts or a Hodge structure on the cohomology of a certain toric variety. In general, this gives a class of ``anti-Ramsey theory problems'' with a geometric structure. For example, we find one for Ramsey numbers of classes of such graphs. The key observation is that our previous result expressing simplicial chromatic polynomials in terms of $h$-vectors of auxiliary simplicial complexes can be reinterpreted as one on edge colorings of graphs avoiding monochromatic colorings of specified forbidden subgraphs. Specializing to simplicial complexes arising from triangulations of polytopes (e.g. unimodular triangulations), we obtain families of graphs and forbidden subgraphs where edge colorings avoiding monochromatic colorings of the forbidden subgraphs depend on lattice point counts or Hodge structures on the cohomology of toric varieties.


title: Signatures of regular black holes from the quasar continuum spectrum

abstract: Regular black holes arising in Einstein gravity coupled to non-linear electrodynamics are worth studying as they can circumvent the r = 0 curvature singularity arising in general relativity. In this work we explore the signatures of regular black holes with a Minkowski core from the quasar continuum spectrum. We use thin-disk approximation to derive the theoretical luminosity from the accretion disk and compare it with the optical data of eighty Palomar Green quasars. Our analysis based on error estimators like the chi-square, the Nash-Sutcliffe efficiency, the index of agreement etc. reveal that optical observations of quasars favor the Kerr scenario compared to black holes in non-linear electrodynamics. The implications are discussed.


title: Grad-GradaGrad? A Non-Monotone Adaptive Stochastic Gradient Method

abstract: The classical AdaGrad method adapts the learning rate by dividing by the square root of a sum of squared gradients. Because this sum on the denominator is increasing, the method can only decrease step sizes over time, and requires a learning rate scaling hyper-parameter to be carefully tuned. To overcome this restriction, we introduce GradaGrad, a method in the same family that naturally grows or shrinks the learning rate based on a different accumulation in the denominator, one that can both increase and decrease. We show that it obeys a similar convergence rate as AdaGrad and demonstrate its non-monotone adaptation capability with experiments.


title: Conservative Hamiltonian Monte Carlo

abstract: We introduce a new class of Hamiltonian Monte Carlo (HMC) algorithm called Conservative Hamiltonian Monte Carlo (CHMC), where energy-preserving integrators, derived from the Discrete Multiplier Method, are used instead of symplectic integrators. Due to the volume being no longer preserved under such a proposal map, a correction involving the determinant of the Jacobian of the proposal map is introduced within the acceptance probability of HMC. For a $p$-th order accurate energy-preserving integrator using a time step size $\tau$, we show that CHMC satisfies stationarity without detailed balance. Moreover, we show that CHMC satisfies approximate stationarity with an error of $\mathcal{O}(\tau^{(m+1)p})$ if the determinant of the Jacobian is truncated to its first $m+1$ terms of its Taylor polynomial in $\tau^p$. We also establish a lower bound on the acceptance probability of CHMC which depends only on the desired tolerance $\delta$ for the energy error and approximate determinant. In particular, a cost-effective and gradient-free version of CHMC is obtained by approximating the determinant of the Jacobian as unity, leading to an $\mathcal{O}(\tau^p)$ error to the stationary distribution and a lower bound on the acceptance probability depending only on $\delta$. Furthermore, numerical experiments show increased performance in acceptance probability and convergence to the stationary distribution for the Gradient-free CHMC over HMC in high dimensional problems.


title: Three-point correlation functions in the $\mathfrak{sl}_3$ Toda theory I: Reflection coefficients

abstract: Toda Conformal Field Theories (CFTs) form a family of 2d CFTs indexed by semisimple and complex Lie algebras. They are natural generalizations of the Liouville CFT in that they enjoy an enhanced level of symmetry encoded by W-algebras. These theories can be rigorously defined using a probabilistic framework that involves the consideration of correlated Gaussian Multiplicative Chaos measures. This document provides a first step towards the computation of a class of three-point correlation functions, that generalize the celebrated DOZZ formula and whose expressions were predicted in the physics literature by Fateev-Litvinov, within the probabilistic framework associated to the $\mathfrak{sl}_3$ Toda CFT. Namely this first article of a two-parts series is dedicated to the probabilistic derivation of the reflection coefficients of general Toda CFTs, which are essential building blocks in the understanding of Toda correlation functions. Along the computations of these reflection coefficients a new path decomposition for diffusion processes in Euclidean spaces, based on a suitable notion of minimum and that generalizes the celebrated one-dimensional result of Williams, will be unveiled. As a byproduct we describe the joint tail expansion of correlated Gaussian Multiplicative Chaos measures together with an asymptotic expansion of class one Whitakker functions.


title: A Local Optima Network Analysis of the Feedforward Neural Architecture Space

abstract: This study investigates the use of local optima network (LON) analysis, a derivative of the fitness landscape of candidate solutions, to characterise and visualise the neural architecture space. The search space of feedforward neural network architectures with up to three layers, each with up to 10 neurons, is fully enumerated by evaluating trained model performance on a selection of data sets. Extracted LONs, while heterogeneous across data sets, all exhibit simple global structures, with single global funnels in all cases but one. These results yield early indication that LONs may provide a viable paradigm by which to analyse and optimise neural architectures.


title: On cohomological and formal properties of Strong K\"ahler with torsion and astheno-K\"ahler metrics

abstract: We provide families of compact astheno-K\"ahler nilmanifolds and we study the behaviour of the complex blowup of such manifolds. We prove that the existence of an astheno-K\"ahler metric satisfying an extra differential condition is not preserved by blowup. We also study the interplay between Strong K\"ahler with torsion metrics and geometrically Bott-Chern metrics. We show that Fino-Parton-Salamon nilmanifolds are geometrically-Bott-Chern-formal, whereas we obtain negative results on the product of two copies of primary Kodaira surface, Inoue surface of type $\mathcal{S}_M$ and on the product of a Kodaira surface with an Inoue surface.


title: A note on the representation theory of central extensions of reductive $p$-adic groups

abstract: In this note, we verify that several fundamental results from the theory of representations of reductive $p$-adic groups, extend to finite central extensions of these groups.


title: Giant mid-IR resonant coupling to molecular vibrations in sub-nm gaps of plasmonic multilayer metafilms

abstract: Nanomaterials capable of confining light are desirable for enhancing spectroscopies such as Raman scattering, infrared absorption, and nonlinear optical processes. Plasmonic superlattices have shown the ability to host collective resonances in the mid-infrared, but require stringent fabrication processes to create well-ordered structures. Here, we demonstrate how short-range-ordered Au nanoparticle multilayers on a mirror, self-assembled by a sub-nm molecular spacer, support collective plasmon-polariton resonances in the visible and infrared, continuously tunable beyond 11 $\mu$m by simply varying the nanoparticle size and number of layers. The resulting molecule-plasmon system approaches vibrational strong coupling, and displays giant Fano dip strengths, SEIRA enhancement factors ~10$^6$, light-matter coupling strengths g~100 cm$^{-1}$, Purcell factors ~10$^6$, and mode volume compression factors ~10$^8$. The collective plasmon-polariton mode is highly robust to nanoparticle vacancy disorder and is sustained by the consistent gap size defined by the molecular spacer. Structural disorder efficiently couples light into the gaps between the multilayers and mirror, enabling Raman and infrared sensing of sub-picolitre sample volumes.


title: Bounds on higher graph gonality

abstract: We prove new lower and upper bounds on the higher gonalities of finite graphs. These bounds are generalizations of known upper and lower bounds for first gonality to higher gonalities, including upper bounds on gonality involving independence number, and lower bounds on gonality by scramble number. We apply our bounds to study the computational complexity of computing higher gonalities, proving that it is NP-hard to compute the second gonality of a graph when restricting to multiplicity-free divisors.


title: LPCSE: Neural Speech Enhancement through Linear Predictive Coding

abstract: The increasingly stringent requirement on quality-of-experience in 5G/B5G communication systems has led to the emerging neural speech enhancement techniques, which however have been developed in isolation from the existing expert-rule based models of speech pronunciation and distortion, such as the classic Linear Predictive Coding (LPC) speech model because it is difficult to integrate the models with auto-differentiable machine learning frameworks. In this paper, to improve the efficiency of neural speech enhancement, we introduce an LPC-based speech enhancement (LPCSE) architecture, which leverages the strong inductive biases in the LPC speech model in conjunction with the expressive power of neural networks. Differentiable end-to-end learning is achieved in LPCSE via two novel blocks: a block that utilizes the expert rules to reduce the computational overhead when integrating the LPC speech model into neural networks, and a block that ensures the stability of the model and avoids exploding gradients in end-to-end training by mapping the Linear prediction coefficients to the filter poles. The experimental results show that LPCSE successfully restores the formants of the speeches distorted by transmission loss, and outperforms two existing neural speech enhancement methods of comparable neural network sizes in terms of the Perceptual evaluation of speech quality (PESQ) and Short-Time Objective Intelligibility (STOI) on the LJ Speech corpus.


title: Krylov subspace residual and restarting for certain second order differential equations

abstract: We propose algorithms for efficient time integration of large systems of oscillatory second order ordinary differential equations (ODEs) whose solution can be expressed in terms of trigonometric matrix functions. Our algorithms are based on a residual notion for second order ODEs, which allows to extend the ``residual-time restarting'' Krylov subspace framework -- which was recently introduced for exponential and $\varphi$-functions occurring in time integration of first order ODEs -- to our setting. We then show that the computational cost can be further reduced in many cases by using our restarting in the Gautschi cosine scheme. We analyze residual convergence in terms of Faber and Chebyshev series and supplement these theoretical results by numerical experiments illustrating the efficiency of the proposed methods.


title: Stability of solitary waves in nonlinear Klein-Gordon equations

abstract: The stability of topological solitary waves and pulses in one-dimensional nonlinear Klein-Gordon systems is revisited. The linearized equation describing small deviations around the static solution leads to a Sturm-Liouville problem, which is solved in a systematic way for the $-l\,(l+1)\,\sech^2(x)$-potential, showing the orthogonality and completeness relations fulfilled by the set of its solutions for all values $l\in\mathbb{N}$. This approach allows to determine the linear stability of kinks and pulses of certain nonlinear Klein-Gordon equations. Two families of novel nonlinear Klein-Gordon potentials are introduced. The exact solutions (kinks and pulses) for these potentials are exactly calculated, even when the nonlinear potential is not explicitly known. The kinks of the novel models are found to be stable, whereas the pulses are unstable. The stability of the pulses is achieved by introducing certain spatial inhomogeneities.