title: Bolder is Better: Raising User Awareness through Salient and Concise Privacy Notices

abstract: This paper addresses the question whether the recently proposed approach of concise privacy notices in apps and on websites is effective in raising user awareness. To assess the effectiveness in a realistic setting, we included concise notices in a fictitious but realistic fitness tracking app and asked participants recruited from an online panel to provide their feedback on the usability of the app as a cover story. Importantly, after giving feedback, users were also asked to recall the data practices described in the notices. The experimental setup included the variation of different levels of saliency and riskiness of the privacy notices. Based on a total sample of 2,274 participants, our findings indicate that concise privacy notices are indeed a promising approach to raise user awareness for privacy information when displayed in a salient way, especially in case the notices describe risky data practices. Our results may be helpful for regulators, user advocates and transparency-oriented companies in creating or enforcing better privacy transparency towards average users that do not read traditional privacy policies.


title: Diffuse-Interface Blended Method for Imposing Physical Boundaries in Two-Fluid Flows

abstract: Multiphase flows are commonly found in chemical engineering processes such as distillation columns, bubble columns, fluidized beds and heat exchangers. The physical boundaries of domains in numerical simulations of multiphase flows are generally defined by a conformal unstructured mesh which, depending on the complexity of the physical system, results in time-consuming mesh generation which frequently requires user-intervention. Furthermore, the resulting conformal unstructured mesh could potentially contain a large number of skewed elements, which is undesirable for numerical stability and accuracy. The diffuse-interface approach allows for the use of a simple structured meshes to be used while still capturing the desired physical (e.g., solid-fluid) boundaries. In this work, a novel diffuse-interface method for the imposition of physical boundaries is developed for the incompressible two-fluid multiphase flow model. This model is appropriate for dispersed multiphase flows which are pervasive in chemical engineering processes, in that this flow regime results in high levels of mass and energy transfer between phases. A diffuse interface is used to define the physical boundaries and boundary conditions are imposed by blending the conservation equations from the two-fluid model with that of the nondeformable solid. The results from the diffuse-interface method are compared with results from a conformal unstructured mesh for different interface functions and widths. For small interface widths, the accuracy of the flow profile is unaffected by the choice of interface function and the phase fraction distribution and flow behavior are within 3% compared to those from a conformal mesh. As the interface width increases, the diffuse-interface solution deviates from the conformal mesh solution in both the localized gas fraction and the overall gas hold-up, resulting in a difference up to 30%.


title: Analysis of Bragg Curve Parameters and Lateral Straggle for Proton and Carbon Beams

abstract: Heavy ions have varying effects on the target. The most important factor in comparing this effect is Linear Energy Transfer (LET). Protons and carbons are heavy ions with high LET. Since these ions lose energy through collisions as they move through the tissue, their range is not long. This loss of energy increases along the way, and the maximum energy loss is reached at the end of the range. This whole process is represented by the Bragg curve. The input dose of the Bragg curve, full width at half maximum (FWHM) value, Bragg peak amplitude and position, and Penumbra thickness are important factors in determining which particle is advantageous in tumor treatment. While heavy ions move through the tissue, small deviations occur in their direction of travel due to Coulomb collisions. These small deviations cause lateral straggle in the dose profile. Lateral straggle is important in determining the type and energy of the particle used in tumor treatments close to critical organs. In our study, when the water phantom of protons and carbon beams with different energies is taken into consideration, the input dose, FWHM value, peak amplitude and position, penumbra thickness and lateral straggle are calculated using the TRIM code and the results are compared with Monte Carlo (MC) simulation. It was found that the proton has an average of 63% more FWHM and 53% more Penumbra than the carbon ion. The carbon ion has an average of 28-45 times greater Bragg peak amplitude at the same Bragg peak location than the proton. It was observed that the proton scattered approximately 70% more in lateral straggle. The difference was found to be around 1.32 mm. In line with all these results, the most commonly used proton and carbon heavy ions in hadron therapy applications were compared.


title: Scalable Deep Compressive Sensing

abstract: Deep learning has been used to image compressive sensing (CS) for enhanced reconstruction performance. However, most existing deep learning methods train different models for different subsampling ratios, which brings additional hardware burden. In this paper, we develop a general framework named scalable deep compressive sensing (SDCS) for the scalable sampling and reconstruction (SSR) of all existing end-to-end-trained models. In the proposed way, images are measured and initialized linearly. Two sampling masks are introduced to flexibly control the subsampling ratios used in sampling and reconstruction, respectively. To make the reconstruction model adapt to any subsampling ratio, a training strategy dubbed scalable training is developed. In scalable training, the model is trained with the sampling matrix and the initialization matrix at various subsampling ratios by integrating different sampling matrix masks. Experimental results show that models with SDCS can achieve SSR without changing their structure while maintaining good performance, and SDCS outperforms other SSR methods.


title: Observational constraints on the formation and evolution of Neptune-class exoplanets

abstract: Among exoplanets, the small-size population constitutes the dominant one, with a diversity of properties and compositions ranging from rocky to gas dominated envelope. While a large fraction of them have masses and radii similar to or smaller than Neptune, yet none share common properties in term of orbital period and insulation with our ice giants. These exoplanets belong to multi-planet systems where planets are closely packed within the first tenth of AU and often exposed to strong irradiation from their host star. Their formation process, subsequent evolution, and fate are still debated and trigger new developments of planet formation models. This paper reviews the characteristics and properties of this extended sample of planets with radii between $\sim$ 1.6 and 4.0$_\oplus$. Even though we still lack real Neptune/Uranus analogues, these exoplanets provide us with key observational constraints that allow the formation of our ice giants to be placed in a more general framework than the sole example of our solar system.


title: Size and shape fluctuations of ultrasoft colloids

abstract: Ultrasoft colloidal particle fluctuates due to its flexibility. Such fluctuation is essential for colloidal structure and dynamics, but is challenging to quantify experimentally. We use dendrimers as a model system to study the fluctuation of ultrasoft colloids. By considering the dynamic polydispersity in the small-angle neutron scattering (SANS) model, and introducing the fluctuation of invasive water into the contrast in SANS, we reveal the fluctuating amplitudes of the size and shape of the dendrimer of generation 6 at finite concentrations. The size fluctuation is suppressed while the shape fluctuation increases as the weight fraction of dendrimers passes 11%. With neutron spin echo data, we suggest that such crossover originates from the competition between the inter- and intra-particle dynamics. Further investigation on lower-generation samples shows a contrary result, which suggests a structural basis for these dynamic phenomena.


title: Data-Driven Distributionally Robust Optimization for Real-Time Economic Dispatch Considering Secondary Frequency Regulation Cost

abstract: With the large-scale integration of renewable power generation, frequency regulation resources (FRRs) are required to have larger capacities and faster ramp rates, which increases the cost of the frequency regulation ancillary service. Therefore, it is necessary to consider the frequency regulation cost and constraint along with real-time economic dispatch (RTED). In this paper, a data-driven distributionally robust optimization (DRO) method for RTED considering automatic generation control (AGC) is proposed. First, a Copula-based AGC signal model is developed to reflect the correlations among the AGC signal, load power and renewable generation variations. Secondly, samples of the AGC signal are taken from its conditional probability distribution under the forecasted load power and renewable generation variations. Thirdly, a distributionally robust RTED model considering the frequency regulation cost and constraint is built and transformed into a linear programming problem by leveraging the Wasserstein metric-based DRO technique. Simulation results show that the proposed method can reduce the total cost of power generation and frequency regulation.


title: An Alternative to Collective Coordinates

abstract: Collective coordinates provide a powerful tool for separating collective and elementary excitations, allowing both to be treated in the full quantum theory. The price is a canonical transformation which leads to a complicated starting point for subsequent calculations. Sometimes the collective behavior of a soliton is simple but nontrivial, and one is interested in the elementary excitations. We show that in this case an alternative prescription suffices, in which the canonical transformation is not necessary. The use of a nonperturbative operator which creates a soliton state allows the theory to be constructed perturbatively in terms of the soliton normal modes. We show how translation invariance may be perturbatively imposed. We apply this to construct the two-loop ground state of an arbitrary scalar kink.


title: Cavity and replica methods for the spectral density of sparse symmetric random matrices

abstract: We review the problem of how to compute the spectral density of sparse symmetric random matrices, i.e. weighted adjacency matrices of undirected graphs. Starting from the Edwards-Jones formula, we illustrate the milestones of this line of research, including the pioneering work of Bray and Rodgers using replicas. We focus first on the cavity method, showing that it quickly provides the correct recursion equations both for single instances and at the ensemble level. We also describe an alternative replica solution that proves to be equivalent to the cavity method. Both the cavity and the replica derivations allow us to obtain the spectral density via the solution of an integral equation for an auxiliary probability density function. We show that this equation can be solved using a stochastic population dynamics algorithm, and we provide its implementation. In this formalism, the spectral density is naturally written in terms of a superposition of local contributions from nodes of given degree, whose role is thoroughly elucidated. This paper does not contain original material, but rather gives a pedagogical overview of the topic. It is indeed addressed to students and researchers who consider entering the field. Both the theoretical tools and the numerical algorithms are discussed in detail, highlighting conceptual subtleties and practical aspects.


title: Adversarial Attacks for Tabular Data: Application to Fraud Detection and Imbalanced Data

abstract: Guaranteeing the security of transactional systems is a crucial priority of all institutions that process transactions, in order to protect their businesses against cyberattacks and fraudulent attempts. Adversarial attacks are novel techniques that, other than being proven to be effective to fool image classification models, can also be applied to tabular data. Adversarial attacks aim at producing adversarial examples, in other words, slightly modified inputs that induce the Artificial Intelligence (AI) system to return incorrect outputs that are advantageous for the attacker. In this paper we illustrate a novel approach to modify and adapt state-of-the-art algorithms to imbalanced tabular data, in the context of fraud detection. Experimental results show that the proposed modifications lead to a perfect attack success rate, obtaining adversarial examples that are also less perceptible when analyzed by humans. Moreover, when applied to a real-world production system, the proposed techniques shows the possibility of posing a serious threat to the robustness of advanced AI-based fraud detection procedures.


title: Observation of thermalization and information scrambling in a superconducting quantum processor

abstract: Understanding various phenomena in non-equilibrium dynamics of closed quantum many-body systems, such as quantum thermalization, information scrambling, and nonergodic dynamics, is a crucial for modern physics. Using a ladder-type superconducting quantum processor, we perform analog quantum simulations of both the $XX$ ladder and one-dimensional (1D) $XX$ model. By measuring the dynamics of local observables, entanglement entropy and tripartite mutual information, we signal quantum thermalization and information scrambling in the $XX$ ladder. In contrast, we show that the $XX$ chain, as free fermions on a 1D lattice, fails to thermalize, and local information does not scramble in the integrable channel. Our experiments reveal ergodicity and scrambling in the controllable qubit ladder, and opens the door to further investigations on the thermodynamics and chaos in quantum many-body systems.


title: Riemannian Manifold Optimization for Discriminant Subspace Learning

abstract: Linear discriminant analysis (LDA) is a widely used algorithm in machine learning to extract a low-dimensional representation of high-dimensional data, it features to find the orthogonal discriminant projection subspace by using the Fisher discriminant criterion. However, the traditional Euclidean-based methods for solving LDA are easily convergent to spurious local minima and hardly obtain an optimal solution. To address such a problem, in this paper, we propose a novel algorithm namely Riemannian-based discriminant analysis (RDA) for subspace learning. In order to obtain an explicit solution, we transform the traditional Euclidean-based methods to the Riemannian manifold space and use the trust-region method to learn the discriminant projection subspace. We compare the proposed algorithm to existing variants of LDA, as well as the unsupervised tensor decomposition methods on image classification tasks. The numerical results suggest that RDA achieves state-of-the-art performance in classification accuracy.


title: Lucky planets: how circum-binary planets survive the supernova in one of the inner-binary components

abstract: A planet hardly ever survives the supernova of the host star in a bound orbit, because mass loss in the supernova and the natal kick imparted to the newly formed compact object cause the planet to be ejected. A planet in orbit around a binary has a considerably higher probability to survive the supernova explosion of one of the inner binary stars. In those cases, the planet most likely remains bound to the companion of the exploding star, whereas the compact object is ejected. We estimate this to happen to $\sim 1/33$ the circum-binary planetary systems. These planetary orbits tend to be highly eccentric ($e \apgt 0.9$), and $\sim 20$\,\% of these planets have retrograde orbits compared to their former binary. The probability that the planet as well as the binary (now with a compact object) remains bound is about ten times smaller ($\sim 3\cdot 10^{-3}$). We then expect the Milky way Galaxy to host $\aplt 10$ x-ray binaries that are still orbited by a planet, and $\aplt 150$ planets that survived in orbit around the compact object's companion. These numbers should be convolved with the fraction of massive binaries that is orbited by a planet.


title: A Statistical Analysis of the Effect of Different Embolization Materials on Gamma Knife Arteriovenous Malformation Dose Distributions

abstract: Gamma Knife Radiosurgery (GKRS) is a treatment choice in newly diagnosed AVMs and residual AVMs following endovascular interventions. Aim of this study is to determine if commercially available liquid embolic agents reduce the radiation dose to the target due to the attenuation of the 60Co beam. Doses accumulated by three different embolization materials, namely: Onyx, cyanoacrylate and polyvinyl alcohol (PVA) were analysed. A collimator helmet size of 8mm was employed in a Monte-Carlo simulation implemented in Geant4 simulation toolkit. Obtained dose accumulations were evaluated and results demonstrated statistically significant differences in the dosimetries calculated.


title: Bias in ontologies -- a preliminary assessment

abstract: Logical theories in the form of ontologies and similar artefacts in computing and IT are used for structuring, annotating, and querying data, among others, and therewith influence data analytics regarding what is fed into the algorithms. Algorithmic bias is a well-known notion, but what does bias mean in the context of ontologies that provide a structuring mechanism for an algorithm's input? What are the sources of bias there and how would they manifest themselves in ontologies? We examine and enumerate types of bias relevant for ontologies, and whether they are explicit or implicit. These eight types are illustrated with examples from extant production-level ontologies and samples from the literature. We then assessed three concurrently developed COVID-19 ontologies on bias and detected different subsets of types of bias in each one, to a greater or lesser extent. This first characterisation aims contribute to a sensitisation of ethical aspects of ontologies primarily regarding representation of information and knowledge.


title: Weighted value distributions of the Riemann zeta function on the critical line

abstract: We prove a central limit theorem for $\log|\zeta(1/2+it)|$ with respect to the measure $|\zeta^{(m)}(1/2+it)|^{2k}dt$ ($k,m\in\mathbb N$), assuming RH and the asymptotic formula for twisted and shifted integral moments of zeta. Under the same hypotheses, we also study a shifted case, looking at the measure $|\zeta(1/2+it+i\alpha)|^{2k}dt$, with $\alpha\in(-1,1)$. Finally we prove unconditionally the analogue result in the random matrix theory context.


title: Effects of internal dynamics on chemotactic aggregation of bacteria

abstract: The effects of internal adaptation dynamics on the self-organized aggregation of chemotactic bacteria are investigated by Monte Carlo (MC) simulations based on a two-stream kinetic transport equation coupled with a reaction-diffusion equation of the chemoattractant that bacteria produce. A remarkable finding is a nonmonotonic behavior of the peak aggregation density with respect to the adaptation time; more specifically, aggregation is the most enhanced when the adaptation time is comparable to or moderately larger than the mean run time of bacteria. Another curious observation is the formation of a trapezoidal aggregation profile occurring at a very large adaptation time, where the biased motion of individual cells is rather hindered at the plateau regimes due to the boundedness of the tumbling frequency modulation. Asymptotic analysis of the kinetic transport system is also carried out, and a novel asymptotic equation is obtained at the large adaptation-time regime while the Keller-Segel type equations are obtained when the adaptation time is moderate. Numerical comparison of the asymptotic equations with MC results clarifies that trapezoidal aggregation is well described by the novel asymptotic equation, and the nonmonotonic behavior of the peak aggregation density is interpreted as the transient of the asymptotic solutions between different adaptation time regimes.


title: Infinitely many twin prime polynomials of odd degree

abstract: While the twin prime conjecture is still famously open, it holds true in the setting of finite fields: There are infinitely many pairs of monic irreducible polynomials over $\mathbb{F}_q$ that differ by a fixed constant, for each $q \geq 3$. Elementary, constructive proofs were given for different cases by Hall and Pollack. In the same spirit, we discuss the construction of a further infinite family of twin prime tuples of odd degree, and its relations to the existence of certain Wieferich primes and to arithmetic properties of the combinatorial Bell numbers.


title: Bridge the Vision Gap from Field to Command: A Deep Learning Network Enhancing Illumination and Details

abstract: With the goal of tuning up the brightness, low-light image enhancement enjoys numerous applications, such as surveillance, remote sensing and computational photography. Images captured under low-light conditions often suffer from poor visibility and blur. Solely brightening the dark regions will inevitably amplify the blur, thus may lead to detail loss. In this paper, we propose a simple yet effective two-stream framework named NEID to tune up the brightness and enhance the details simultaneously without introducing many computational costs. Precisely, the proposed method consists of three parts: Light Enhancement (LE), Detail Refinement (DR) and Feature Fusing (FF) module, which can aggregate composite features oriented to multiple tasks based on channel attention mechanism. Extensive experiments conducted on several benchmark datasets demonstrate the efficacy of our method and its superiority over state-of-the-art methods.


title: 1st Place Solution to ECCV-TAO-2020: Detect and Represent Any Object for Tracking

abstract: We extend the classical tracking-by-detection paradigm to this tracking-any-object task. Solid detection results are first extracted from TAO dataset. Some state-of-the-art techniques like \textbf{BA}lanced-\textbf{G}roup \textbf{S}oftmax (\textbf{BAGS}\cite{li2020overcoming}) and DetectoRS\cite{qiao2020detectors} are integrated during detection. Then we learned appearance features to represent any object by training feature learning networks. We ensemble several models for improving detection and feature representation. Simple linking strategies with most similar appearance features and tracklet-level post association module are finally applied to generate final tracking results. Our method is submitted as \textbf{AOA} on the challenge website. Code is available at https://github.com/feiaxyt/Winner_ECCV20_TAO.


title: One-dimensional game-theoretic differential equations

abstract: We provide a very brief introduction to typical paths and the corresponding It\^o type integration. Relying on this robust It\^o integration, we prove an existence and uniqueness result for one-dimensional differential equations driven by typical paths with non-Lipschitz continuous coefficients in the spirit of Yamada--Watanabe as well as an approximation result in the spirit of Doss--Sussmann.


title: The geodesic-transversal problem

abstract: A maximal geodesic in a graph is a geodesic (alias shortest path) which is not a subpath of a longer geodesic. The geodesic-transversal problem in a graph $G$ is introduced as the task to find a smallest set $S$ of vertices of $G$ such that each maximal geodesic has at least one vertex in $S$. The minimum cardinality of such a set is the geodesic-transversal number ${\rm gt}(G)$ of $G$. It is proved that ${\rm gt}(G) = 1$ if and only if $G$ is a subdivided star and that the geodesic-transversal problem is NP-complete. Fast algorithms to determine the geodesic-transversal number of trees and of spread cactus graphs are designed, respectively.


title: Long term X-Ray Observations of Seyfert 1 Galaxy Ark 120: On the origin of soft-excess

abstract: We present the long-term X-ray spectral and temporal analysis of a 'bare-type AGN' Ark 120. We consider the observations from XMM-Newton, Suzaku, Swift, and NuSTAR from 2003 to 2018. The spectral properties of this source are studied using various phenomenological and physical models present in the literature. We report (a) the variations of several physical parameters, such as the temperature and optical depth of the electron cloud, the size of the Compton cloud, and accretion rate for the last fifteen years. The spectral variations are explained from the change in the accretion dynamics; (b) the X-ray time delay between 0.2-2 keV and 3-10 keV light-curves exhibited zero-delay in 2003, positive delay of 4.71 \pm 2.1 ks in 2013, and negative delay of 4.15 \pm 1.5 ks in 2014. The delays are explained considering Comptonization, reflection, and light-crossing time; (c) the long term intrinsic luminosities obtained using nthcomp, of the soft-excess and the primary continuum show a correlation with a Pearson Correlation Coefficient of 0.922. This indicates that the soft-excess and the primary continuum are originated from the same physical process. From a physical model fitting, we infer that the soft excess for Ark 120 could be due to a small number of scatterings in the Compton cloud. Using Monte-Carlo simulations, we show that indeed the spectra corresponding to fewer scatterings could provide a steeper soft-excess power-law in the 0.2-3 keV range. Simulated luminosities are found to be in agreement with the observed values.


title: Bayesian Optimization Assisted Meal Bolus Decision Based on Gaussian Processes Learning and Risk-Sensitive Control

abstract: Effective postprandial glucose control is important to glucose management for subjects with diabetes mellitus. In this work, a data-driven meal bolus decision method is proposed without the need of subject-specific glucose management parameters. The postprandial glucose dynamics is learnt using Gaussian process regression. Considering the asymmetric risks of hyper- and hypoglycemia and the uncertainties in the predicted glucose trajectories, an asymmetric risk-sensitive cost function is designed. Bayesian optimization is utilized to solve the optimization problem, since the gradient of the cost function is unavailable. The proposed approach is evaluated using the 10-adult cohort of the FDA-accepted UVA/Padova T1DM simulator and compared with the standard insulin bolus calculator. For the case of announced meals, the proposed method achieves satisfactory and similar performance in terms of mean glucose and percentage time in [70, 180] mg/dL without increasing the risk of hypoglycemia. Similar results are observed for the case without the meal information (assuming that the patient follows a consistent diet) and the case of basal rate mismatches. In addition, advisory-mode analysis is performed based on clinical data, which indicates that the method can determine safe and reasonable meal boluses in real clinical settings. The results verify the effectiveness and robustness of the proposed method and indicate the feasibility of achieving improved postprandial glucose regulation through a data-driven optimal control method.


title: A class of Newton maps with Julia sets of Lebesgue measure zero

abstract: Let $g(z)=\int_0^zp(t)\exp(q(t))\,dt+c$ where $p,q$ are polynomials and $c\in\mathbb{C}$, and let $f$ be the function from Newton's method for $g$. We show that under suitable assumptions the Julia set of $f$ has Lebesgue measure zero. Together with a theorem by Bergweiler, our result implies that $f^n(z)$ converges to zeros of $g$ almost everywhere in $\mathbb{C}$ if this is the case for each zero of $g''$. In order to prove our result, we establish general conditions ensuring that Julia sets have Lebesgue measure zero.


title: Extended Reality (XR) Remote Research: a Survey of Drawbacks and Opportunities

abstract: Extended Reality (XR) technology - such as virtual and augmented reality - is now widely used in Human Computer Interaction (HCI), social science and psychology experimentation. However, these experiments are predominantly deployed in-lab with a co-present researcher. Remote experiments, without co-present researchers, have not flourished, despite the success of remote approaches for non-XR investigations. This paper summarises findings from a 30-item survey of 46 XR researchers to understand perceived limitations and benefits of remote XR experimentation. Our thematic analysis identifies concerns common with non-XR remote research, such as participant recruitment, as well as XR-specific issues, including safety and hardware variability. We identify potential positive affordances of XR technology, including leveraging data collection functionalities builtin to HMDs (e.g. hand, gaze tracking) and the portability and reproducibility of an experimental setting. We suggest that XR technology could be conceptualised as an interactive technology and a capable data-collection device suited for remote experimentation.


title: Two generalized strong convergence algorithms for the variational inequality problems in Banach spaces

abstract: In this paper, two generalized algorithms for solving the variational inequality problem in Banach spaces are proposed. Then the strong convergence of the sequences generated by these algorithms will be proved under the suitable conditions. Finally, using MATLAB software, we provide some numerical examples to illustrate our results.


title: Exploring Design and Governance Challenges in the Development of Privacy-Preserving Computation

abstract: Homomorphic encryption, secure multi-party computation, and differential privacy are part of an emerging class of Privacy Enhancing Technologies which share a common promise: to preserve privacy whilst also obtaining the benefits of computational analysis. Due to their relative novelty, complexity, and opacity, these technologies provoke a variety of novel questions for design and governance. We interviewed researchers, developers, industry leaders, policymakers, and designers involved in their deployment to explore motivations, expectations, perceived opportunities and barriers to adoption. This provided insight into several pertinent challenges facing the adoption of these technologies, including: how they might make a nebulous concept like privacy computationally tractable; how to make them more usable by developers; and how they could be explained and made accountable to stakeholders and wider society. We conclude with implications for the development, deployment, and responsible governance of these privacy-preserving computation techniques.


title: Evaluating uncertainties in electrochemical impedance spectra of solid oxide fuel cells

abstract: Electrochemical impedance spectroscopy (EIS) is a widely used tool for characterization of fuel cells and other electrochemical conversion systems. When applied to the on-line monitoring in the context of in-field applications, the disturbances, drifts and sensor noise may cause severe distortions in the evaluated spectra, especially in the low-frequency part. Failure to ignore the random effects can result in misinterpreted spectra and, consequently, in misleading diagnostic reasoning. This fact has not been often addressed in the research so far. In this paper, we propose an approach to the quantification of the spectral uncertainty, which relies on evaluating the uncertainty of the equivalent circuit model (ECM). We apply the computationally efficient variational Bayes (VB) method and compare the quality of the results with those obtained with the Markov chain Monte Carlo (MCMC) algorithm. Namely, MCMC algorithm returns accurate distributions of the estimated model parameters, while VB approach provides the approximate distributions. By using simulated and real data we show that approximate results provided by VB approach, although slightly over-optimistic, are still close to the more realistic MCMC estimates. A great advantage of the VB method for online monitoring is low computational load, which is several orders of magnitude lower compared to MCMC. The performance of VB algorithm is demonstrated on a case of ECM parameters estimation in a 6 cell solid oxide fuel cell (SOFC) stack. The complete numerical implementation for recreating the results can be found at https://repo.ijs.si/lznidaric/variational-bayes-supplementary-material.


title: Quantification and interpretation of the climate variability record

abstract: The spectral view of variability is a compelling and adaptable tool for understanding variability of the climate. In Mitchell (1976) seminal paper, it was used to express, on one graph with log scales, a very wide range of climate variations from millions of years to days. The spectral approach is particularly useful for suggesting causal links between forcing variability and climate response variability. However, a substantial degree of variability is intrinsic and the Earth system may respond to external forcing in a complex manner. There has been an enormous amount of work on understanding climate variability over the last decades. Hence in this paper, we address the question: Can we (after 40 years) update the Mitchell (1976) diagram and provide it with a better interpretation? By reviewing both the extended observations available for such a diagram and new methodological developments in the study of the interaction between internal and forced variability over a wide range of timescales, we give a positive answer to this question. In addition, we review alternative approaches to the spectral decomposition and pose some challenges for a more detailed quantification of climate variability.


title: Broadband Non-Geostationary Satellite Communication Systems: Research Challenges and Key Opportunities

abstract: Besides conventional geostationary (GSO) satellite broadband communication services, non-geostationary (NGSO) satellites are envisioned to support various new communication use cases from countless industries. These new scenarios bring many unprecedented challenges that will be discussed in this paper alongside with several potential future research opportunities. NGSO systems are known for various advantages, including their important features of low cost, lower propagation delay, smaller size, and lower losses in comparison to GSO satellites. However, there are still many deployment challenges to be tackled to ensure seamless integration not only with GSO systems but also with terrestrial networks. In this paper, we discuss several key challenges including satellite constellation and architecture designs, coexistence with GSO systems in terms of spectrum access and regulatory issues, resource management algorithms, and NGSO networking requirements. Additionally, the latest progress in provisioning secure communication via NGSO systems is discussed. Finally, this paper identifies multiple important open issues and research directions to inspire further studies towards the next generation of satellite networks.


title: Variational Autoencoders with a Structural Similarity Loss in Time of Flight MRAs

abstract: Time-of-Flight Magnetic Resonance Angiographs (TOF-MRAs) enable visualization and analysis of cerebral arteries. This analysis may indicate normal variation of the configuration of the cerebrovascular system or vessel abnormalities, such as aneurysms. A model would be useful to represent normal cerebrovascular structure and variabilities in a healthy population and to differentiate from abnormalities. Current anomaly detection using autoencoding convolutional neural networks usually use a voxelwise mean-error for optimization. We propose optimizing a variational-autoencoder (VAE) with structural similarity loss (SSIM) for TOF-MRA reconstruction. A patch-trained 2D fully-convolutional VAE was optimized for TOF-MRA reconstruction by comparing vessel segmentations of original and reconstructed MRAs. The method was trained and tested on two datasets: the IXI dataset, and a subset from the ADAM challenge. Both trained networks were tested on a dataset including subjects with aneurysms. We compared VAE optimization with L2-loss and SSIM-loss. Performance was evaluated between original and reconstructed MRAs using mean square error, mean-SSIM, peak-signal-to-noise-ratio and dice similarity index (DSI) of segmented vessels. The L2-optimized VAE outperforms SSIM, with improved reconstruction metrics and DSIs for both datasets. Optimization using SSIM performed best for visual image quality, but with discrepancy in quantitative reconstruction and vascular segmentation. The larger, more diverse IXI dataset had overall better performance. Reconstruction metrics, including SSIM, were lower for MRAs including aneurysms. A SSIM-optimized VAE improved the visual perceptive image quality of TOF-MRA reconstructions. A L2-optimized VAE performed best for TOF-MRA reconstruction, where the vascular segmentation is important. SSIM is a potential metric for anomaly detection of MRAs.


title: Fast formation and assembly of isogeometric Galerkin matrices for trimmed patches

abstract: This work explores the application of the fast assembly and formation strategy from [8, 17] to trimmed bi-variate parameter spaces. Two concepts for the treatment of basis functions cut by the trimming curve are investigated: one employs a hybrid Gauss-point-based approach, and the other computes discontinuous weighted quadrature rules. The concepts' accuracy and efficiency are examined for the formation of mass matrices and their application to L2-projection. Significant speed-ups compared to standard element by element finite element formation are observed. There is no clear preference between the concepts proposed. While the discontinuous weighted scheme scales favorably with the degree of the basis, it also requires additional effort for computing the quadrature weights. The hybrid Gauss approach does not have this overhead, which is determined by the complexity of the trimming curve. Hence, it is well-suited for moderate degrees, whereas discontinuous-weightedquadrature has potential for high degrees, in particular, if the related weights are computed in parallel.


title: Voltage Inference for and Coordination of Distributed Voltage Controls in Extremely-High DER-Penetration Distribution Networks

abstract: The unique problems and phenomena in the distributed voltage control of large-scale power distribution systems with extremely-high DER-penetration are targeted in this paper. First, a DER-explicit distribution network model and voltage sensitivity are derived. Based on that, a voltage inference method is implemented to fill the gap of measurement insufficiency in the grid-edge areas. Then, autonomous Q control being implemented in each local area, a $\overline{Q}$-coordinated P control is designed to coordinate the reactive and real power controls. All the algorithms have been tested in standard and synthetic systems, and have expected results. Moreover, an open-source software platform, which integrates the modeling of communication networks, DER controls, and power networks, is developed to enable the distributed control and optimization algorithms in the grid simulation of the large-scale distribution systems.


title: Effects of dynamical screening on the BCS-BEC crossover in double bilayer graphene: Density functional theory for exciton bilayers

abstract: We derive a gap equation for bilayer excitonic systems based on density functional theory and benchmark our results against quantum Monte-Carlo simulations and recent experiments on double bilayer graphene. The gap equation has a mean-field form but includes a consistent treatment of dynamical screening. We show that the gap survives at much higher densities than previously thought from mean-field estimates which gives strong indications that the double-bilayer graphene systems at zero magnetic field can be used as model systems to investigate the BCS-BEC crossover. Furthermore, we show that Josephson-like transfer of pairs can be substantial for small band gaps and densities.


title: On the Parametrization and Statistics of Propagation Graphs

abstract: Propagation graphs (PGs) serve as a frequency-selective, spatially consistent channel model suitable for fast channel simulations in a scattering environment. So far, however, the parametrization of the model, and its consequences, have received little attention. In this contribution, we propose a new parametrization for PGs that adheres to the doubly exponentially decaying cluster structure of the Saleh-Valenzuela (SV) model. We show how to compute the newly proposed internal model parameters based on an approximation of the $K$-factor and the two decay rates from the SV model. Furthermore, via the singular values of multiple-input multiple-output (MIMO) channels, we compare the degrees of freedom (DoF) between our new and another frequently used parametrization. Specifically, we compare the DoF loss when the distance between antennas within the transmitter and receiver arrays or the average distance between scatterers decreases. Based on this comparison, it is shown that, in contrast to the typical parametrization, our newly proposed parametrization loses DoF in both scenarios, as one would expect from a spatially consistent channel model.


title: Weak Convergence for Variational Inequalities with Inertial-Type Method

abstract: Weak convergence of inertial iterative method for solving variational inequalities is the focus of this paper. The cost function is assumed to be non-Lipschitz and monotone. We propose a projection-type method with inertial terms and give weak convergence analysis under appropriate conditions. Some test results are performed and compared with relevant methods in the literature to show the efficiency and advantages given by our proposed methods.


title: On the cubic Weyl sum

abstract: We obtain an estimate for the cubic Weyl sum which improves the bound obtained from Weyl differencing for short ranges of summation. In particular, we show that for any $\varepsilon>0$ there exists some $\delta>0$ such that for any coprime integers $a,q$ and real number $\gamma$ we have \begin{align*} \sum_{1\le n \le N}e\left(\frac{an^3}{q}+\gamma n\right)\ll (qN)^{1/4} q^{-\delta}, \end{align*} provided $q^{1/3+\varepsilon}\le N \le q^{1/2-\varepsilon}$. Our argument builds on some ideas of Enflo.


title: Description of longitudinal modes in moderately coupled Yukawa systems with the static local field correction

abstract: In moderately coupled Yukawa fluids, longitudinal mode dispersion is determined by the competition between kinetic and potential effects. In a recent paper [Khrapak and Cou\"edel, Phys. Rev. E 102, 033207 (2020)], a semi-phenomenological dispersion relation was constructed by the ad-hoc addition of the Bohm-Gross kinetic term to the generalized instantaneous excess bulk modulus, which showed very good agreement with simulations. In this paper, a nearly identical dispersion relation is derived in a rigorous manner based on a dielectric formulation with static local field corrections. At moderate coupling, this formalism is revealed to be more accurate than other successful theoretical approaches.


title: The Collatz Problem generalized to 3x+k

abstract: The Collatz problem with $3x+k$ is revisited. Positive and negative limit cycles are given up to k=9997 starting with $x_0=-2\cdot10^7...+2\cdot10^7$. A simple relation between the probability distribution for the Syracuse iterates for various k (not divisible by 2 and 3) is obtained. From this it follows that the oscillation considered by Tao 2019 ( arXiv:1909.03562v2 ) does not depend on k. Thus this piece of the proof of his theorem 1.3 "Almost all Collatz orbits attain almost bounded values" holds for all k not divisible by 2 and 3.


title: A Similarity Measure of Gaussian Process Predictive Distributions

abstract: Some scenarios require the computation of a predictive distribution of a new value evaluated on an objective function conditioned on previous observations. We are interested on using a model that makes valid assumptions on the objective function whose values we are trying to predict. Some of these assumptions may be smoothness or stationarity. Gaussian process (GPs) are probabilistic models that can be interpreted as flexible distributions over functions. They encode the assumptions through covariance functions, making hypotheses about new data through a predictive distribution by being fitted to old observations. We can face the case where several GPs are used to model different objective functions. GPs are non-parametric models whose complexity is cubic on the number of observations. A measure that represents how similar is one GP predictive distribution with respect to another would be useful to stop using one GP when they are modelling functions of the same input space. We are really inferring that two objective functions are correlated, so one GP is enough to model both of them by performing a transformation of the prediction of the other function in case of inverse correlation. We show empirical evidence in a set of synthetic and benchmark experiments that GPs predictive distributions can be compared and that one of them is enough to predict two correlated functions in the same input space. This similarity metric could be extremely useful used to discard objectives in Bayesian many-objective optimization.


title: Thread Evolution Kit for Optimizing Thread Operations on CE/IoT Devices

abstract: Most modern operating systems have adopted the one-to-one thread model to support fast execution of threads in both multi-core and single-core systems. This thread model, which maps the kernel-space and user-space threads in a one-to-one manner, supports quick thread creation and termination in high-performance server environments. However, the performance of time-critical threads is degraded when multiple threads are being run in low-end CE devices with limited system resources. When a CE device runs many threads to support diverse application functionalities, low-level hardware specifications often lead to significant resource contention among the threads trying to obtain system resources. As a result, the operating system encounters challenges, such as excessive thread context switching overhead, execution delay of time-critical threads, and a lack of virtual memory for thread stacks. This paper proposes a state-of-the-art Thread Evolution Kit (TEK) that consists of three primary components: a CPU Mediator, Stack Tuner, and Enhanced Thread Identifier. From the experiment, we can see that the proposed scheme significantly improves user responsiveness (7x faster) under high CPU contention compared to the traditional thread model. Also, TEK solves the segmentation fault problem that frequently occurs when a CE application increases the number of threads during its execution.


title: Component Tree Loss Function: Definition and Optimization

abstract: In this article, we propose a method to design loss functions based on component trees which can be optimized by gradient descent algorithms and which are therefore usable in conjunction with recent machine learning approaches such as neural networks. We show how the altitudes associated to the nodes of such hierarchical image representations can be differentiated with respect to the image pixel values. This feature is used to design a generic loss function that can select or discard image maxima based on various attributes such as extinction values. The possibilities of the proposed method are demonstrated on simulated and real image filtering.


title: Interpolation by multivariate polynomials in convex domains

abstract: Let $\Omega$ be a convex open set in $\mathbb R^n$ and let $\Lambda_k$ be a finite subset of $\Omega$. We find necessary geometric conditions for $\Lambda_k$ to be interpolating for the space of multivariate polynomials of degree at most $k$. Our results are asymptotic in $k$. The density conditions obtained match precisely the necessary geometric conditions that sampling sets are known to satisfy, and they are expressed in terms of the equilibrium potential of the convex set. Moreover, we prove that in the particular case of the unit ball, for $k$ large enough, there is no family of orthogonal reproducing kernels in the space of polynomials of degree at most $k$.


title: The functional form of Mahler conjecture for even log-concave functions in dimension $2$

abstract: Let $\Phi$ : R n $\rightarrow$ R $\cup$ {+$\infty$} be an even convex function and L$\Phi$ be its Legendre transform. We prove the functional form of Mahler conjecture concerning the functional volume product P ($\Phi$) = e --$\Phi$ e --L$\Phi$ in dimension 2: we give the sharp lower bound of this quantity and characterize the equality case. The proof uses the computation of the derivative in t of P (t$\Phi$) and ideas due to Meyer [M] for unconditional convex bodies, adapted to the functional case by Fradelizi-Meyer [FM2] and extended for symmetric convex bodies in dimension 3 by Iriyeh-Shibata [IS] (see also [FHMRZ]).


title: A note on Reidemeister Torsion of G-Anosov Representations

abstract: This article considers $G$-Anosov representations of a fixed closed oriented Riemann surface $\Sigma$ of genus at least $2$. Here, $G$ is the Lie group $\text{PSp}(2n,\mathbb{R}$), $\text{PSO}(n,n)$ or $\text{PSO}(n,n+1)$. It proves that Reidemeister torsion (R-torsion) associated to $\Sigma$ with coefficients in the adjoint bundle representations of such representations is well-defined. Moreover, by using symplectic chain complex method, it establishes a novel formula for R-torsion of such representations in terms of the Atiyah-Bott-Goldman symplectic form corresponding to the Lie group $G$. Furtermore, it applies the results to Hitchin components, in particular, Teichm\"{u}ller space.


title: Non-divisible point on a two-parameter family of elliptic curves

abstract: Let n be a positive integer and t a non-zero integer. We consider the elliptic curve over Q given by E : y 2 = x 3 + tx 2 -- n 2 (t + 3n 2)x + n 6. It is a special case of an elliptic surface studied recently by Bettin, David and Delaunay [2] and it generalizes Washington's family. The point (0, n 3) belongs to E(Q) and we obtain some results about its nondivisibility in E(Q). Our work extends to this two-parameter family of elliptic curves a previous study of Duquesne (mainly stated for n = 1 and t > 0).


title: Neural networks-based algorithms for stochastic control and PDEs in finance

abstract: This paper presents machine learning techniques and deep reinforcement learningbased algorithms for the efficient resolution of nonlinear partial differential equations and dynamic optimization problems arising in investment decisions and derivative pricing in financial engineering. We survey recent results in the literature, present new developments, notably in the fully nonlinear case, and compare the different schemes illustrated by numerical tests on various financial applications. We conclude by highlighting some future research directions.


title: The Tully-Fisher relation in dense groups at $z \sim 0.7$ in the MAGIC survey

abstract: Galaxies in dense environments are subject to interactions and mechanisms which directly affect their evolution by lowering their gas fractions and reducing their star-forming capacity earlier than their isolated counterparts. The aim of our project is to get new insights about the role of environment on the stellar and baryonic content of galaxies using a kinematic approach, through the study of the Tully-Fisher relation (TFR). We study a sample of galaxies in 8 groups spanning a redshift range of $0.5<z<0.8$ and located in 10 pointings of the MAGIC MUSE Guaranteed Time Observations program. We perform a morpho-kinematics analysis of this sample and set up a selection based on galaxy size, [OII] emission line doublet signal-to-noise ratio, bulge-to-disk ratio and nuclear activity to construct a robust kinematic sample of 67 star-forming galaxies. This selection considerably reduces the number of outliers in the TFR, which are predominantly dispersion-dominated galaxies. Our results suggest a significant offset of the TFR zero-point between galaxies in low- and high-density environments, whatever kinematics estimator is used. This can be interpreted as a decrease of either stellar mass by $\sim 0.05 - 0.3$ dex or an increase of rotation velocity by $\sim 0.02 - 0.06$ dex for galaxies in groups, depending on the samples used for comparison. We also studied the stellar and baryon mass fractions within stellar disks and found they both increase with stellar mass, the trend being more pronounced for the stellar component alone. These fractions do not exceed 50%. We show that this evolution of the TFR is consistent either with a decrease of star formation or with a contraction of the mass distribution due to the environment. These two effects probably act together with their relative contribution depending on the mass regime.


title: Towards the Right Direction in BiDirectional User Interfaces

abstract: Hundreds of millions of speakers of bidirectional (BiDi) languages rely on writing systems that mix the native right-to-left script with left-to-right strings. The global reach of interactive digital technologies requires special attention to these people, whose perception of interfaces is affected by this script mixture. However, empirical research on this topic is scarce. Although leading software vendors provide guidelines for BiDi design, bidirectional interfaces demonstrate inconsistent and incorrect directionality of UI elements, which may cause user confusion and errors. Through a websites' review, we identified problematic UI items and considered reasons for their existence. In an online survey with 234 BiDi speakers, we observed that in many cases, users' direction preferences were inconsistent with the guidelines. The findings provide potential insights for design rules and empirical evidence for the problem's complexity, suggesting the need for further empirical research and greater attention by the HCI community to the BiDi design problem.


title: The electric dipole moment of the tau lepton revisited

abstract: We reconsider the issue of the search for a nonzero electric dipole form factor (EDM) $d_\tau(s)$ using optimal observables in $\tau^+\tau^-$ production by $e^+ e^-$ collisions in the center-of-mass energy range from the $\tau$-pair threshold to about $\sqrt{s} \sim 15$ GeV. We discuss the general formalism of optimal observables and apply it to two $CP$-odd observables that are sensitive to the real and imaginary part of $d_\tau(s)$, respectively. We compute the expectation values and covariances of these optimal $CP$ observables for $\tau$-pair production at $\sqrt{s}=10.58$ GeV with subsequent decays of $\tau^\pm$ into major leptonic or semihadronic modes. For the $\tau$ decays to two pions and three charged pions we take the full kinematic information of the hadronic system into account. Assuming that the Belle II experiment at the KEKB accelerator will eventually analyze data corresponding to an integrated luminosity of 50 ab$^{-1}$ and applying acceptance cuts on the final-state pions we find that 1~s.d. sensitivities $\delta {\rm Re} d_\tau = 6.8 \times 10^{-20}$ e cm and $\delta {\rm} Im d_\tau = 4.0 \times 10^{-20}$ e cm can be obtained with events where both $\tau$'s decay semihadronically. In the ideal case where no cuts on the final-state particles are applied we find with 50 ab$^{-1}$ at $\sqrt{s}=10.58$ GeV corresponding to $4.5 \times 10^{10}$ $\tau^+ \tau^-$ events the 1 s.d. sensitivities $\delta {\rm Re} d_\tau = 5.8 \times 10^{-20}$ e cm and $\delta {\rm} Im d_\tau = 3.2 \times 10^{-20}$ e cm, again for events where both $\tau$ leptons decay semihadronically. Furthermore, we analyze the potential magnitude of the $\tau$ EDM form factor in the type-II two-Higgs doublet extension and in two scalar leptoquark extensions of the Standard Model.


title: A practical preconditioner for wide-field continuum imaging of radio interferometric data

abstract: The celebrated CLEAN algorithm has been the cornerstone of deconvolution algorithms in radio interferometry almost since its conception in the 1970s. For all its faults, CLEAN is remarkably fast, robust to calibration artefacts and in its ability to model point sources. We demonstrate how the same assumptions that afford CLEAN its speed can be used to accelerate more sophisticated deconvolution algorithms.


title: The stochastic Auxiliary Problem Principle in Banach spaces: measurability and convergence

abstract: The stochastic Auxiliary Problem Principle (APP) algorithm is a general Stochastic Approximation (SA) scheme that turns the resolution of an original optimization problem into the iterative resolution of a sequence of auxiliary problems. This framework has been introduced to design decomposition-coordination schemes but also encompasses many well-known SA algorithms such as stochastic gradient descent or stochastic mirror descent. We study the stochastic APP in the case where the iterates lie in a Banach space and we consider an additive error on the computation of the subgradient of the objective. In order to derive convergence results or efficiency estimates for a SA scheme, the iterates must be random variables. This is why we prove the measurability of the iterates of the stochastic APP algorithm. Then, we extend convergence results from the Hilbert space case to the Banach space case. Finally, we derive efficiency estimates for the function values taken at the averaged sequence of iterates or at the last iterate, the latter being obtained by adapting the concept of modified Fej{\'e}r monotonicity to our framework.


title: Flocking and Collision Avoidance for a Dynamic Squad of Fixed-Wing UAVs Using Deep Reinforcement Learning

abstract: Developing the flocking behavior for a dynamic squad of fixed-wing UAVs is still a challenge due to kinematic complexity and environmental uncertainty. In this paper, we deal with the decentralized flocking and collision avoidance problem through deep reinforcement learning (DRL). Specifically, we formulate a decentralized DRL-based decision making framework from the perspective of every follower, where a collision avoidance mechanism is integrated into the flocking controller. Then, we propose a novel reinforcement learning algorithm PS-CACER for training a shared control policy for all the followers. Besides, we design a plug-n-play embedding module based on convolutional neural networks and the attention mechanism. As a result, the variable-length system state can be encoded into a fixed-length embedding vector, which makes the learned DRL policy independent with the number and the order of followers. Finally, numerical simulation results demonstrate the effectiveness of the proposed method, and the learned policies can be directly transferred to semi-physical simulation without any parameter finetuning.


title: Asymptotic first boundary value problem for elliptic operators

abstract: In 1955, Lehto showed that, for every measurable function $\psi$ on the unit circle $\mathbb T,$ there is a function $f$ holomorphic in the unit disc, having $\psi$ as radial limit a.e. on $\mathbb T.$ We consider an analogous problem for solutions $f$ of homogenous elliptic equations $Pf=0$ and, in particular, for holomorphic functions on Riemann surfaces and harmonic functions on Riemannian manifolds.


title: Fluctuation theory for one-sided L\'evy processes with a matrix-exponential time horizon

abstract: There is an abundance of useful fluctuation identities for one-sided L\'evy processes observed up to an independent exponentially distributed time horizon. We show that all the fundamental formulas generalize to time horizons having matrix exponential distributions, and the structure is preserved. Essentially, the positive killing rate is replaced by a matrix with eigenvalues in the right half of the complex plane which, in particular, applies to the positive root of the Laplace exponent and the scale function. Various fundamental properties of thus obtained matrices and functions are established, resulting in an easy to use toolkit. An important application concerns deterministic time horizons which can be well approximated by concentrated matrix exponential distributions. Numerical illustrations are also provided.


title: Upstream mobility Finite Volumes for the Richards equation in heterogenous domains

abstract: This paper is concerned with the Richards equation in a heterogeneous domain, each subdomain of which is homogeneous and represents a rocktype. Our first contribution is to rigorously prove convergence toward a weak solution of cell-centered finite-volume schemes with upstream mobility and without Kirchhoff's transform. Our second contribution is to numerically demonstrate the relevance of locally refining the grid at the interface between subregions, where discontinuities occur, in order to preserve an acceptable accuracy for the results computed with the schemes under consideration.


title: Populations facing a nonlinear environmental gradient: steady states and pulsating fronts

abstract: We consider a population structured by a spacevariable and a phenotypical trait, submitted to dispersion,mutations, growth and nonlocal competition. This population is facing an {\it environmental gradient}: to survive at location $x$, an individual must have a trait close to some optimal trait $y_{opt}(x)$. Our main focus is to understand the effect of a {\it nonlinear} environmental gradient. We thus consider a nonlocal parabolic equation for the distribution of the population, with $y_{opt}(x) = \varepsilon\theta(x)$, $0<\vert \varepsilon \vert \ll 1$. We construct steady states solutions and, when $\theta$ is periodic, pulsating fronts. This requires the combination of rigorous perturbation techniques based on a careful application of the implicit function theorem in rather intricate function spaces. To deal with the phenotypic trait variable $y$ we take advantage of a Hilbert basis of $L^{2}(\mathbb{R})$ made of eigenfunctions of an underlying Schr\"odinger operator, whereas to deal with the space variable $x$ we use the Fourier series expansions. Our mathematical analysis reveals, in particular, how both the steady states solutions and the fronts (speed and profile) are distorted by the nonlinear environmental gradient, which are important biological insights.


title: NB-IoT Random Access for Non-Terrestrial Networks

abstract: The satellite component is recognized as a promising solution to complement and extend the coverage of future Internet of things (IoT) terrestrial networks (TNs). In this context, a study item to integrate satellites into narrowband IoT (NB-IoT) systems has been approved within the 3rd Generation Partnership Project (3GPP) standardization body. However, as NB-IoT systems were initially conceived for TNs, their basic design principles and operation might require some key modifications when incorporating the satellite component. These changes in NB-IoT systems, therefore, need to be carefully implemented in order to guarantee a seamless integration of both TN and non-terrestrial network (NTN) for a global coverage. This paper addresses this adaptation for the random access (RA) step in NB-IoT systems, which is in fact the most challenging aspect in the NTN context, for it deals with multi-user time-frequency synchronization and timing advance for data scheduling. In particular, we propose an RA technique which is robust to typical satellite channel impairments, including long delays, significant Doppler effects, and wide beams, without requiring any modification to the current NBIoT RA waveform. Performance evaluations demonstrate the proposal's capability of addressing different NTN configurations recently defined by 3GPP for the 5G new radio system.


title: A Damped Newton Algorithm for Generated Jacobian Equations

abstract: Generated Jacobian Equations have been introduced by Trudinger [Disc. cont. dyn. sys (2014), pp. 1663-1681] as a generalization of Monge-Amp{\`e}re equations arising in optimal transport. In this paper, we introduce and study a damped Newton algorithm for solving these equations in the semi-discrete setting, meaning that one of the two measures involved in the problem is finitely supported and the other one is absolutely continuous. We also present a numerical application of this algorithm to the near-field parallel refractor problem arising in non-imaging problems.


title: Recursive Trellis Processing of Large Polarization Kernels

abstract: A reduced complexity algorithm is presented for computing the log-likelihood ratios arising in the successive cancellation decoder for polar codes with large kernels of arbitrary dimension. The proposed algorithm exploits recursive trellis representation of the codes generated by submatrices of the polarization kernel, and enables codes based on large kernels to provide better performance compared to the codes based on Arikan kernel with the same decoding complexity.


title: Matrix Kesten Recursion, Inverse-Wishart Ensemble and Fermions in a Morse Potential

abstract: The random variable $1+z_1+z_1z_2+\dots$ appears in many contexts and was shown by Kesten to exhibit a heavy tail distribution. We consider natural extensions of this variable and its associated recursion to $N \times N$ matrices either real symmetric $\beta=1$ or complex Hermitian $\beta=2$. In the continuum limit of this recursion, we show that the matrix distribution converges to the inverse-Wishart ensemble of random matrices. The full dynamics is solved using a mapping to $N$ fermions in a Morse potential, which are non-interacting for $\beta=2$. At finite $N$ the distribution of eigenvalues exhibits heavy tails, generalizing Kesten's results in the scalar case. The density of fermions in this potential is studied for large $N$, and the power-law tail of the eigenvalue distribution is related to the properties of the so-called determinantal Bessel process which describes the hard edge universality of random matrices. For the discrete matrix recursion, using free probability in the large $N$ limit, we obtain a self-consistent equation for the stationary distribution. The relation of our results to recent works of Rider and Valk\'o, Grabsch and Texier, as well as Ossipov, is discussed.

categories: cond-mat.stat-mech cond-mat.dis-nn math-ph math.MP math.PR


title: Developments and applications of Shapley effects to reliability-oriented sensitivity analysis with correlated inputs

abstract: Reliability-oriented sensitivity analysis methods have been developed for understanding the influence of model inputs relative to events which characterize the failure of a system (e.g., a threshold exceedance of the model output). In this field, the target sensitivity analysis focuses primarily on capturing the influence of the inputs on the occurrence of such a critical event. This paper proposes new target sensitivity indices, based on the Shapley values and called "target Shapley effects", allowing for interpretable sensitivity measures under dependent inputs. Two algorithms (one based on Monte Carlo sampling, and a given-data algorithm based on a nearest-neighbors procedure) are proposed for the estimation of these target Shapley effects based on the $\ell^2$ norm. Additionally, the behavior of these target Shapley effects are theoretically and empirically studied through various toy-cases. Finally, the application of these new indices in two real-world use-cases (a river flood model and a COVID-19 epidemiological model) is discussed.


title: Absorptive laser threshold magnetometry: combining visible diamond Raman lasers and nitrogen-vacancy centres

abstract: We propose a high-sensitivity magnetometry scheme based on a diamond Raman laser with visible pump absorption by an ensemble of coherently microwave driven negatively charged nitrogen-vacancy centres (NV) in the same diamond crystal. The NV centres' absorption and emission are spin-dependent. We show how the varying absorption of the NV centres changes the Raman laser output. A shift in the diamond Raman laser threshold and output occurs with the external magnetic-field and microwave driving. We develop a theoretical framework including steady-state solutions to describe the effects of coherently driven NV centres in a diamond Raman laser. We discuss that such a laser working at the threshold can be employed for magnetic-field sensing. In contrast to previous studies on NV magnetometry with visible laser absorption, the laser threshold magnetometry method is expected to have low technical noise, due to low background light in the measurement signal. For magnetic-field sensing, we project a shot-noise limited DC sensitivity of a few $\mathrm{pT}/\sqrt{\mathrm{Hz}}$ in a well-calibrated cavity with realistic parameters. This sensor employs the broad visible absorption of NV centres and unlike previous laser threshold magnetometry proposals it does not rely on active NV centre lasing or an infrared laser medium at the specific wavelength of the NV centre's infrared absorption line.


title: Few-shot Action Recognition with Prototype-centered Attentive Learning

abstract: Few-shot action recognition aims to recognize action classes with few training samples. Most existing methods adopt a meta-learning approach with episodic training. In each episode, the few samples in a meta-training task are split into support and query sets. The former is used to build a classifier, which is then evaluated on the latter using a query-centered loss for model updating. There are however two major limitations: lack of data efficiency due to the query-centered only loss design and inability to deal with the support set outlying samples and inter-class distribution overlapping problems. In this paper, we overcome both limitations by proposing a new Prototype-centered Attentive Learning (PAL) model composed of two novel components. First, a prototype-centered contrastive learning loss is introduced to complement the conventional query-centered learning objective, in order to make full use of the limited training samples in each episode. Second, PAL further integrates a hybrid attentive learning mechanism that can minimize the negative impacts of outliers and promote class separation. Extensive experiments on four standard few-shot action benchmarks show that our method clearly outperforms previous state-of-the-art methods, with the improvement particularly significant (10+\%) on the most challenging fine-grained action recognition benchmark.


title: Qudits for Witnessing Quantum Gravity Induced Entanglement of Masses Under Decoherence

abstract: Recently a theoretical and an experimental protocol known as quantum gravity induced entanglement of masses (QGEM) has been proposed to test the quantum nature of gravity using two mesoscopic masses each placed in a superposition of two locations. If, after eliminating all non-gravitational interactions between them, the particles become entangled, one can conclude that the gravitational potential is induced via a quantum mediator, i.e. a virtual graviton. In this paper, we examine a range of different experimental set-ups, considering different geometries and the number of spatially superposed states taken, in order to determine which would generate entanglement faster. We conclude that without decoherence, and given a maximum distance $\Delta x$ between any two spatial states of a superposition, a set of two qubits placed in spatial superposition parallel to one another will outperform all other models given realistic experimental parameters. Furthermore, when a sufficiently high decoherence rate is introduced, multi-component superpositions can outperform the two-qubit set-up. This is further verified with an experimental simulation, showing that $O(10^3)$ measurements are required to reject the no entanglement hypothesis with a parallel qubits set-up without decoherence at a 99.9$\%$ confidence level. The number of measurements increases when decoherence is introduced. When the decoherence rate reaches $0.125$~Hz, 6-dimensional qudits are required as the two-qubit system entanglement cannot be witnessed anymore. However, in this case, $O(10^6)$ measurements will be required. One can group the witness operators to measure in order to reduce the number of measurements (up to ten-fold). However, this may be challenging to implement experimentally.


title: The Challenges of Persian User-generated Textual Content: A Machine Learning-Based Approach

abstract: Over recent years a lot of research papers and studies have been published on the development of effective approaches that benefit from a large amount of user-generated content and build intelligent predictive models on top of them. This research applies machine learning-based approaches to tackle the hurdles that come with Persian user-generated textual content. Unfortunately, there is still inadequate research in exploiting machine learning approaches to classify/cluster Persian text. Further, analyzing Persian text suffers from a lack of resources; specifically from datasets and text manipulation tools. Since the syntax and semantics of the Persian language is different from English and other languages, the available resources from these languages are not instantly usable for Persian. In addition, recognition of nouns and pronouns, parts of speech tagging, finding words' boundary, stemming or character manipulations for Persian language are still unsolved issues that require further studying. Therefore, efforts have been made in this research to address some of the challenges. This presented approach uses a machine-translated datasets to conduct sentiment analysis for the Persian language. Finally, the dataset has been rehearsed with different classifiers and feature engineering approaches. The results of the experiments have shown promising state-of-the-art performance in contrast to the previous efforts; the best classifier was Support Vector Machines which achieved a precision of 91.22%, recall of 91.71%, and F1 score of 91.46%.


title: Frequency-comb-induced radiation pressure force in dense atomic clouds

abstract: We investigate the frequency comb induced radiation pressure force acting on a cloud of cold $^{87}$Rb atoms. Reduction and spectral broadening of the frequency comb force are observed as the cloud's optical thickness is increased. Since the radiation pressure force is uniquely determined by light scattered by an atomic cloud, we discuss different scattering mechanisms, and point to the shadow effect as the dominant mechanism affecting FC-induced force in resonantly excited dense atomic clouds. Our results improve the understanding of the interaction of frequency comb light with many-atom ensembles, which is essential for novel frequency comb applications in simultaneous multi-species cooling, multi-mode quantum memories, and multi-mode atom-light interfaces.


title: Generalization properties of restricted Boltzmann machine for short-range order

abstract: The restricted Boltzmann machine (RBM) is used to investigate short-range order in binary alloys. The network is trained on the data collected by Monte Carlo simulations for a simple Ising-like binary alloy model and used to calculate the Warren--Cowley short-range order parameter and other thermodynamic properties. We demonstrate that RBM not only reproduces the order parameters for the alloy concentration at which it was trained, but can also predict them for any other concentrations.


title: Discriminant groups of wild cyclic quotient singularities

abstract: Let p be prime. We describe explicitly the resolution of singularities of several families of wild Z/pZ-quotient singularities in dimension two, including families that generalize the quotient singularities of type E_6, E_7, and E_8 from p=2 to arbitrary characteristics. We prove that for odd primes, any power of p can appear as the determinant of the intersection matrix of a wild Z/pZ-quotient singularity. We also provide evidence towards the conjecture that in this situation one may choose the wild action to be ramified precisely at the origin.


title: A survey of joint intent detection and slot-filling models in natural language understanding

abstract: Intent classification and slot filling are two critical tasks for natural language understanding. Traditionally the two tasks have been deemed to proceed independently. However, more recently, joint models for intent classification and slot filling have achieved state-of-the-art performance, and have proved that there exists a strong relationship between the two tasks. This article is a compilation of past work in natural language understanding, especially joint intent classification and slot filling. We observe three milestones in this research so far: Intent detection to identify the speaker's intention, slot filling to label each word token in the speech/text, and finally, joint intent classification and slot filling tasks. In this article, we describe trends, approaches, issues, data sets, evaluation metrics in intent classification and slot filling. We also discuss representative performance values, describe shared tasks, and provide pointers to future work, as given in prior works. To interpret the state-of-the-art trends, we provide multiple tables that describe and summarise past research along different dimensions, including the types of features, base approaches, and dataset domain used.


title: An effective field theory of holographic dark energy

abstract: A general covariant local field theory of the holographic dark energy model is presented. It turns out the low energy effective theory of the holographic dark energy is the massive gravity theory whose graviton has 3 polarisations, including one scalar mode and two tensor modes. The Compton wavelength is the size of the future event horizon of the universe. The UV-IR correspondence in the holographic dark energy model stems from the scalar graviton's strong coupling at the energy scale that marks the breaking down of the effective field theory.


title: A NEAT Quantum Error Decoder

abstract: We investigate the use of the evolutionary NEAT algorithm for the optimization of a policy network that performs quantum error decoding on the toric code, with bitflip and depolarizing noise, one qubit at a time. We find that these NEAT-optimized network decoders have similar performance to previously reported machine-learning based decoders, but use roughly three to four orders of magnitude fewer parameters to do so.


title: Generalized Tur\'an problems for complete bipartite graphs

abstract: For graph $G$, $F$ and integer $n$, the generalized Tu\'an number $ex(n,G,F)$ denotes the maximum number of copies of $G$ that an $F$-free $n$-vertex graph can have. We study this parameter when both $G$ and $F$ are complete bipartite graphs.


title: Automatic Differentiation via Effects and Handlers: An Implementation in Frank

abstract: Automatic differentiation (AD) is an important family of algorithms which enables derivative based optimization. We show that AD can be simply implemented with effects and handlers by doing so in the Frank language. By considering how our implementation behaves in Frank's operational semantics, we show how our code performs the dynamic creation of programs during evaluation.


title: The role of reflections in the generation of a time delay in strong field ionization

abstract: The problem of time delay in tunneling ionization is revisited. The origin of time delay at the tunnel exit is analysed, underlining the two faces of the concept of the tunnelling time delay: the time delay around the tunnel exit and the asymptotic time delay at a detector. We show that the former time delay, in the sense of a delay in the peak of the wavefunction, exists as a matter of principle and arises due to the sub-barrier interference of the reflected and transmitted components of the tunneling electronic wavepacket. We exemplify this by describing the tunnelling ionization of an electron bound by a short-range potential within the strong field approximation in a "deep tunnelling" regime. If sub-barrier reflections are extracted from this wavefunction, then the time delay of the peak is shown to vanish. Thus, we assert that the disturbance of the tunnelling wavepacket by the reflection from the surface of the barrier causes a time delay in the neighbourhood of the tunnel exit.


title: Difference-frequency generation in an AlGaAs Bragg-reflection waveguide using an on-chip electrically-pumped quantum dot laser

abstract: Nonlinear frequency conversion is ubiquitous in laser engineering and quantum information technology. A long-standing goal in photonics is to integrate on-chip semiconductor laser sources with nonlinear optical components. Engineering waveguide lasers with spectra that phase-match to nonlinear processes on the same device is a formidable challenge. Here, we demonstrate difference-frequency generation in an AlGaAs Bragg reflection waveguide which incorporates the gain medium for the pump laser in its core. We include quantum dot layers in the AlGaAs waveguide that generate electrically driven laser light at ~790 nm, and engineer the structure to facilitate nonlinear processes at this wavelength. We perform difference-frequency generation between 1540 nm and 1630 nm using the on-chip laser, which is enabled by the broad modal phase-matching of the AlGaAs waveguide, and measure normalized conversion efficiencies up to $(0.64\pm0.21)$ %/W/cm$^2$. Our work demonstrates a pathway towards devices that utilize on-chip active elements and strong optical nonlinearities to enable highly integrated photonic systems-on-chip.


title: Noncommutative Henselizations

abstract: In this paper, the familiar notion of a Henselian pair is extended to the noncommutative case. Furthermore, the problem of Henselizations is studied in the noncommutative context, and it is shown that every (not necessarily commutative) pair which is Hausdorff with respect to a certain topology has a left (and right) Henselization.


title: Comparison results for solutions to p-Laplace equations with Robin boundary conditions

abstract: In the last decades comparison results of Talenti type for Elliptic Problems with Dirichlet boundary conditions have been widely investigated. In this paper, we generalize the results obtained in arXiv:1909.11950 to the case of p-Laplace operator with Robin boundary conditions. The point-wise comparison, obtained in arXiv:1909.11950 only in the planar case, holds true in any dimension if p is sufficiently small.


title: Active Model Learning using Informative Trajectories for Improved Closed-Loop Control on Real Robots

abstract: Model-based controllers on real robots require accurate knowledge of the system dynamics to perform optimally. For complex dynamics, first-principles modeling is not sufficiently precise, and data-driven approaches can be leveraged to learn a statistical model from real experiments. However, the efficient and effective data collection for such a data-driven system on real robots is still an open challenge. This paper introduces an optimization problem formulation to find an informative trajectory that allows for efficient data collection and model learning. We present a sampling-based method that computes an approximation of the trajectory that minimizes the prediction uncertainty of the dynamics model. This trajectory is then executed, collecting the data to update the learned model. In experiments we demonstrate the capabilities of our proposed framework when applied to a complex omnidirectional flying vehicle with tiltable rotors. Using our informative trajectories results in models which outperform models obtained from non-informative trajectory by 13.3\% with the same amount of training data. Furthermore, we show that the model learned from informative trajectories generalizes better than the one learned from non-informative trajectories, achieving better tracking performance on different tasks.


title: Unicritical Laminations

abstract: Thurston introduced \emph{invariant (quadratic) laminations} in his 1984 preprint as a vehicle for understanding the connected Julia sets and the parameter space of quadratic polynomials. Important ingredients of his analysis of the angle doubling map $\sigma_2$ on the unit circle $\mathbb{S}^1$ were the Central Strip Lemma, non-existence of wandering polygons, the transitivity of the first return map on vertices of periodic polygons, and the non-crossing of minors of quadratic invariant laminations. We use Thurston's methods to prove similar results for \emph{unicritical} laminations of arbitrary degree $d$ and to show that the set of so-called \emph{minors} of unicritical laminations themselves form a \emph{Unicritical Minor Lamination} $\mathrm{UML}_d$. In the end we verify the \emph{Fatou conjecture} for the unicritical laminations and extend the \emph{Lavaurs algorithm} onto $\mathrm{UML}_d$.


title: Designing Accountable Systems

abstract: Accountability is an often called for property of technical systems. It is a requirement for algorithmic decision systems, autonomous cyber-physical systems, and for software systems in general. As a concept, accountability goes back to the early history of Liberalism and is suggested as a tool to limit the use of power. This long history has also given us many, often slightly differing, definitions of accountability. The problem that software developers now face is to understand what accountability means for their systems and how to reflect it in a system's design. To enable the rigorous study of accountability in a system, we need models that are suitable for capturing such a varied concept. In this paper, we present a method to express and compare different definitions of accountability using Structural Causal Models. We show how these models can be used to evaluate a system's design and present a small use case based on an autonomous car.


title: Bottomonium precision tests from full lattice QCD: hyperfine splitting, $\Upsilon$ leptonic width and $b$ quark contribution to $e^+e^- \rightarrow$ hadrons

abstract: We calculate the mass difference between the $\Upsilon$ and $\eta_b$ and the $\Upsilon$ leptonic width from lattice QCD using the Highly Improved Staggered Quark formalism for the $b$ quark and including $u$, $d$, $s$ and $c$ quarks in the sea. We have results for lattices with lattice spacing as low as 0.03 fm and multiple heavy quark masses, enabling us to map out the heavy quark mass dependence and determine values at the $b$ quark mass. Our results are: $M_{\Upsilon} -M_{\eta_b} = 57.5(2.3)(1.0) \,\mathrm{MeV}$ (where the second uncertainty comes from neglect of quark-line disconnected correlation functions) and decay constants, $f_{\eta_b}=724(12)$ MeV and $f_{\Upsilon} =677.2(9.7)$ MeV, giving $\Gamma(\Upsilon \rightarrow e^+e^-) = 1.292(37)(3) \,\mathrm{keV}$. The hyperfine splitting and leptonic width are both in good agreement with experiment, and provide the most accurate lattice QCD results to date for these quantities by some margin. At the same time results for the time moments of the vector-vector correlation function can be compared to values for the $b$ quark contribution to $\sigma(e^+e^- \rightarrow \mathrm{hadrons})$ determined from experiment. Moments 4--10 provide a 2\% test of QCD and yield a $b$ quark contribution to the anomalous magnetic moment of the muon of 0.300(15)$\times 10^{-10}$. Our results, covering a range of heavy quark masses, may also be useful to constrain QCD-like composite theories for beyond the Standard Model physics.


title: A Generalized Weisfeiler-Lehman Graph Kernel

abstract: The Weisfeiler-Lehman graph kernels are among the most prevalent graph kernels due to their remarkable time complexity and predictive performance. Their key concept is based on an implicit comparison of neighborhood representing trees with respect to equality (i.e., isomorphism). This binary valued comparison is, however, arguably too rigid for defining suitable similarity measures over graphs. To overcome this limitation, we propose a generalization of Weisfeiler-Lehman graph kernels which takes into account the similarity between trees rather than equality. We achieve this using a specifically fitted variation of the well-known tree edit distance which can efficiently be calculated. We empirically show that our approach significantly outperforms state-of-the-art methods in terms of predictive performance on datasets containing structurally more complex graphs beyond the typically considered molecular graphs.


title: Holographic Schwinger Effect in Anisotropic Media

abstract: According to gauge/gravity correspondence, we study the holographic Schwinger effect within an anisotropic background. Firstly, the separate length of the particle-antiparticle pairs is computed within the context of an anisotropic background which is parameterized by dynamical exponent $\nu$. It is found that the maximum separate length $x$ increases with the increase of dynamical exponent $\nu$. By analyzing the potential energy, we find that the potential barrier increases with the dynamical exponent $\nu$ at a small separate distance. This observation implies that the Schwinger effect within an anisotropic background is comparatively weaker when contrasted with its manifestation in an isotropic background.Finally, we also find that the Schwinger effect in the transverse direction is weakened compared to the parallel direction in the anisotropic background, which is consistent with the top-down model.


title: Learning to Augment for Data-Scarce Domain BERT Knowledge Distillation

abstract: Despite pre-trained language models such as BERT have achieved appealing performance in a wide range of natural language processing tasks, they are computationally expensive to be deployed in real-time applications. A typical method is to adopt knowledge distillation to compress these large pre-trained models (teacher models) to small student models. However, for a target domain with scarce training data, the teacher can hardly pass useful knowledge to the student, which yields performance degradation for the student models. To tackle this problem, we propose a method to learn to augment for data-scarce domain BERT knowledge distillation, by learning a cross-domain manipulation scheme that automatically augments the target with the help of resource-rich source domains. Specifically, the proposed method generates samples acquired from a stationary distribution near the target data and adopts a reinforced selector to automatically refine the augmentation strategy according to the performance of the student. Extensive experiments demonstrate that the proposed method significantly outperforms state-of-the-art baselines on four different tasks, and for the data-scarce domains, the compressed student models even perform better than the original large teacher model, with much fewer parameters (only ${\sim}13.3\%$) when only a few labeled examples available.


abstract: We classify simple Whittaker modules for classical Lie superalgebras in terms of their parabolic decompositions. We establish a type of Mili\v{c}i\'c-Soergel equivalence of a category of Whittaker modules and a category of Harish-Chandra bimodules. For classical Lie superalgebras of type I, we reduce the problem of composition factors of standard Whittaker modules to that of Verma modules in their BGG categories $\mathcal O$. As a consequence, the composition series of standard Whittaker modules over the general linear Lie superalgebras $\mathfrak{gl}(m|n)$ and the ortho-symplectic Lie superalgebras $\mathfrak{osp}(2|2n)$ can be computed via the Kazhdan-Lusztig combinatorics.


title: Describing limits of integrable functions as grid functions of nonstandard analysis

abstract: In functional analysis, there are different notions of limit for a bounded sequence of $L^1$ functions. Besides the pointwise limit, that does not always exist, the behaviour of a bounded sequence of $L^1$ functions can be described in terms of its weak-$\star$ limit or by introducing a measure-valued notion of limit in the sense of Young measures. Working in Robinson's framework of analysis with infinitesimals, we show that for every bounded sequence $\{z_n\}_{n \in \mathbb{N}}$ of $L^1$ functions there exists a function of a hyperfinite domain (i.e.\ a grid function) that represents both the weak-$\star$ and the Young measure limits of the sequence. This result has relevant applications to the study of nonlinear PDEs. We discuss the example of an ill-posed forward-backward parabolic equation.


title: Joint quasiprobability distribution on the measurement outcomes of MUB-driven operators

abstract: We propose a method to define quasiprobability distributions for general spin-$j$ systems of dimension $n=2j+1$, where $n$ is a prime or power of prime. The method is based on a complete set of orthonormal commuting operators related to Mutually Unbiased Bases which enable (i) a parameterisation of the density matrix and (ii) construction of measurement operators that can be physically realised. As a result we geometrically characterise the set of states for which the quasiprobability distribution is non-negative, and can be viewed as a joint distribution of classical random variables assuming values in a finite set of outcomes. The set is an $(n^2-1)$-dimensional convex polytope with $n+1$ vertices as the only pure states, $n^{n+1}$ number of higher dimensional faces, and $n^3(n+1)/2$ edges.


title: Polynomial interpolation and residue currents

abstract: We show that a global holomorphic section of $\mathscr{O}(d)$ restricted to a closed complex subspace $X \subset \mathbb{P}^n$ has an interpolant if and only if it satisfies a set of moment conditions that involves a residue current associated with a locally free resolution of $\mathscr{O}_X$. When $X$ is a finite set of points in $\mathbb{C}^n \subset \mathbb{P}^n$ this can be interpreted as a set of linear conditions that a function on $X$ has to satisfy in order to have a polynomial interpolant of degree at most $d$.


title: Effective connectivity determines the critical dynamics of biochemical networks

abstract: Living systems operate in a critical dynamical regime -- between order and chaos -- where they are both resilient to perturbation, and flexible enough to evolve. To characterize such critical dynamics, the established 'structural theory' of criticality uses automata network connectivity and node bias (to be on or off) as tuning parameters. This parsimony in the number of parameters needed sometimes leads to uncertain predictions about the dynamical regime of both random and systems biology models of biochemical regulation. We derive a more accurate theory of criticality by accounting for canalization, the existence of redundancy that buffers automata response to inputs -- a known mechanism that buffers the expression of traits, keeping them close to optimal states despite genetic and environmental perturbations. The new 'canalization theory' of criticality is based on a measure of effective connectivity. It contributes to resolving the problem of finding precise ways to design or control network models of biochemical regulation for desired dynamical behavior. Our analyses reveal that effective connectivity significantly improves the prediction of critical behavior in random automata network ensembles. We also show that the average effective connectivity of a large battery of systems biology models is much lower than the connectivity of their original interaction structure. This suggests that canalization has been selected to dynamically reduce and homogenize the seemingly heterogeneous connectivity of biochemical networks.


title: Fundamental measure theory of inhomogeneous two-body correlation functions

abstract: For the three-dimensional hard-sphere model we investigate the inhomogeneous two-body correlations predicted by Rosenfeld's fundamental measure theory. For the special cases in which the density has either planar or spherical symmetry we provide analytic formulae for the Hankel and Legendre transforms, respectively, of the inhomogeneous two-body direct correlation function as explicit functionals of the density. When combined with the Ornstein-Zernike relation our analytical results allow for rapid calculation of inhomogeneous hard-sphere density correlations in real-space. These provide not only information about the packing structures of the hard-sphere system, but also form an essential building-block for constructing perturbation theories of more realistic models.


title: A variational formula for large deviations in First-passage percolation under tail estimates

abstract: Consider first passage percolation with identical and independent weight distributions and first passage time ${\rm T}$. In this paper, we study the upper tail large deviations $\mathbb{P}({\rm T}(0,nx)>n(\mu+\xi))$, for $\xi>0$ and $x\neq 0$ with a time constant $\mu$ and a dimension $d$, for weights that satisfy a tail assumption $ \beta_1\exp{(-\alpha t^r)}\leq \mathbb P(\tau_e>t)\leq \beta_2\exp{(-\alpha t^r)}.$ When $r\leq 1$ (this includes the well-known Eden growth model), we show that the upper tail large deviation decays as $\exp{(-(2d\xi +o(1))n)}$. When $1< r\leq d$, we find that the rate function can be naturally described by a variational formula, called the discrete p-Capacity, and we study its asymptotics. For $r<d$, we show that the large deviation event ${\rm T}(0,nx)>n(\mu+\xi)$ is described by a localization of high weights around the origin. The picture changes for $r\geq d$ where the configuration is not anymore localized.