title: Non-Extensive Behaviour of the QCD Strong Coupling Constant

abstract: The compatibility of theoretically calculated values for the strong coupling constant,through the renormalization group approach with experimental data is studied. There exist considerable divergence in-between theoretical and experimental results at low energies, which can not be explained by thermal field theory and considering chemical potential. Such great deviation can be treated successfully by considering the q-generalized statistical effects through adding a q-nonextensive parameter in the fitting of theoretical results on experimental data.


title: A semi-implicit high-order space-time scheme on staggered meshes for the 2D incompressible Navier-Stokes equations

abstract: A new high order accurate semi-implicit space-time Discontinuous Galerkin method on staggered grids, for the simulation of viscous incompressible flows on two-dimensional domains is presented. The designed scheme is of the Arbitrary Lagrangian Eulerian type, which is suitable to work on fixed as well as on moving meshes. In our space-time formulation, by expressing the numerical solution in terms of piecewise space-time polynomials, an arbitrary high order of accuracy in time is achieved through a simple and efficient method of Picard iterations. For the dual mesh, the basis functions consist in the union of continuous piecewise polynomials on the two subtriangles within the quadrilaterals: this allows the construction of a quadrature-free scheme, resulting in a very efficient algorithm. Some numerical examples confirm that the proposed method outperforms existing ones.


title: Active matter, microreversibility, and thermodynamics

abstract: Active matter, comprising many active agents interacting and moving in fluids or more complex environments, is a commonly occurring state of matter in biological and physical systems. By its very nature active matter systems exist in nonequilibrium states. In this paper the active agents are small Janus colloidal particles that use chemical energy provided by chemical reactions occurring on their surfaces for propulsion through a diffusiophoretic mechanism. As a result of interactions among these colloids, either directly or through fluid velocity and concentration fields, they may act collectively to form structures such as dynamic clusters. A general nonequilibrium thermodynamics framework for the description of such systems is presented that accounts for both self-diffusiophoresis and diffusiophoresis due to external concentration gradients, and is consistent with microreversibility. It predicts the existence of a reciprocal effect of diffusiophoresis back onto the reaction rate for the entire collection of colloids in the system, as well as the existence of a clustering instability that leads to nonequilibrium inhomogeneous system states.

categories: cond-mat.soft cond-mat.stat-mech nlin.PS physics.chem-ph


title: ADW: Blockchain-enabled Small-scale Farm Digitization

abstract: Farm records hold the static, temporal, and longitudinal details of the farms. For small-scale farming, the ability to accurately capture these records plays a critical role in formalizing and digitizing the agriculture industry. Reliable exchange of these record through a trusted platform could unlock critical and valuable insights to different stakeholders across the value chain in agriculture eco-system. Lately, there has been increasing attention on digitization of small scale farming with the objective of providing farm-level transparency, accountability, visibility, access to farm loans, etc. using these farm records. However, most solutions proposed so far have the shortcoming of providing detailed, reliable and trusted small-scale farm digitization information in real time. To address these challenges, we present a system, called Agribusiness Digital Wallet (ADW), which leverages blockchain to formalize the interactions and enable seamless data flow in small-scale farming ecosystem. Utilizing instrumentation of farm tractors, we demonstrate the ability to utilize farm activities to create trusted electronic field records (EFR) with automated valuable insights. Using ADW, we processed several thousands of small-scale farm-level activity events for which we also performed automated farm boundary detection of a number of farms in different geographies.


title: Crowd control: Reducing individual estimation bias by sharing biased social information

abstract: Cognitive biases are widespread in humans and animals alike, and can sometimes be reinforced by social interactions. One prime bias in judgment and decision-making is the human tendency to underestimate large quantities. Previous research on social influence in estimation tasks has generally focused on the impact of single estimates on individual and collective accuracy, showing that randomly sharing estimates does not reduce the underestimation bias. Here, we test a method of social information sharing that exploits the known relationship between the true value and the level of underestimation, and study if it can counteract the underestimation bias. We performed estimation experiments in which participants had to estimate a series of quantities twice, before and after receiving estimates from one or several group members. Our purpose was threefold: to study (i) whether restructuring the sharing of social information can reduce the underestimation bias, (ii) how the number of estimates received affects the sensitivity to social influence and estimation accuracy, and (iii) the mechanisms underlying the integration of multiple estimates. Our restructuring of social interactions successfully countered the underestimation bias. Moreover, we find that sharing more than one estimate also reduces the underestimation bias. Underlying our results are a human tendency to herd, to trust larger estimates than one's own more than smaller estimates, and to follow disparate social information less. Using a computational modelling approach, we demonstrate that these effects are indeed key to explain the experimental results. Overall, our results show that existing knowledge on biases can be used to dampen their negative effects and boost judgment accuracy, paving the way for combating other cognitive biases threatening collective systems.


title: Expected mean width of the randomized integer convex hull

abstract: Let $K \in \R^d$ be a convex body, and assume that $L$ is a randomly rotated and shifted integer lattice. Let $K_L$ be the convex hull of the (random) points $K \cap L$. The mean width $W(K_L)$ of $K_L$ is investigated. The asymptotic order of the mean width difference $W(\l K)-W((\l K)_L)$ is maximized by the order obtained by polytopes and minimized by the order for smooth convex sets as $\l \to \infty$.


title: Gabriel-Zisman Cohomology and spectral sequences

abstract: Extending constructions by Gabriel and Zisman, we develop a functorial framework for the cohomology and homology of simplicial sets with very general coefficient systems given by functors on simplex categories into abelian categories. Furthermore we construct Leray type spectral sequences for any map of simplicial sets. We also show that these constructions generalise and unify the various existing versions of cohomology and homology of small categories and as a bonus provide new insight into their functoriality.


title: Orlicz mixed chord integrals

abstract: We introduce a affine geometric quantity and call it Orlicz mixed chord integral, which generalize the chord integrals to Orlicz space. Minkoswki and Brunn-Minkowski inequalities for the Orlicz mixed chord integrals are establish. These new inequalities in special cases yield some isoperimetric inequalities for the usual chord integrals. The related concepts and inequalities of Lp-mixed chord integrals are also derived.


title: Bounds for exit times of Brownian motion and the first Dirichlet eigenvalue for the Laplacian

abstract: For domains in $\mathbb{R}^d$, $d\geq 2$, we prove universal upper and lower bounds on the product of the bottom of the spectrum for the Laplacian to the power $p>0$ and the supremum over all starting points of the $p$-moments of the exit time of Brownian motion. It is shown that the lower bound is sharp for integer values of $p$ and that for $p \geq 1$, the upper bound is asymptotically sharp as $d\to\infty$. For all $p>0$, we prove the existence of an extremal domain among the class of domains that are convex and symmetric with respect to all coordinate axes. For this class of domains we conjecture that the cube is extremal.


title: Causality-based Explanation of Classification Outcomes

abstract: We propose a simple definition of an explanation for the outcome of a classifier based on concepts from causality. We compare it with previously proposed notions of explanation, and study their complexity. We conduct an experimental evaluation with two real datasets from the financial domain.


title: $L_p$ optimal prediction of the last zero of a spectrally negative L\'evy process

abstract: Given a spectrally negative L\'evy process $X$ drifting to infinity, (inspired on the early ideas of Shiryaev (2002)) we are interested in finding a stopping time that minimises the $L^p$ distance ($p>1$) with $g$, the last time $X$ is negative. The solution is substantially more difficult compared to the case $p=1$, for which it was shown in Baurdoux and Pedraza (2020b) that it is optimal to stop as soon as $X$ exceeds a constant barrier. In the case of $p>1$ treated here, we prove that solving this optimal prediction problem is equivalent to solving an optimal stopping problem in terms of a two-dimensional strong Markov process that incorporates the length of the current excursion away from $0$. We show that an optimal stopping time is now given by the first time that $X$ exceeds a non-increasing and non-negative curve depending on the length of the current excursion away from $0$. We further characterise the optimal boundary and the value function as the unique solution of a non-linear system of integral equations within a subclass of functions. As examples, the case of a Brownian motion with drift and a Brownian motion with drift perturbed by a Poisson process with exponential jumps are considered.


title: HD 145263: Spectral Observations of Silica Debris Disk Formation via Extreme Space Weathering?

abstract: We report here time domain infrared spectroscopy and optical photometry of the HD145263 silica-rich circumstellar disk system taken from 2003 through 2014. We find an F4V host star surrounded by a stable, massive 1e22 - 1e23 kg (M_Moon to M_Mars) dust disk. No disk gas was detected, and the primary star was seen rotating with a rapid ~1.75 day period. After resolving a problem with previously reported observations, we find the silica, Mg-olivine, and Fe-pyroxene mineralogy of the dust disk to be stable throughout, and very unusual compared to the ferromagnesian silicates typically found in primordial and debris disks. By comparison with mid-infrared spectral features of primitive solar system dust, we explore the possibility that HD 145263's circumstellar dust mineralogy occurred with preferential destruction of Fe-bearing olivines, metal sulfides, and water ice in an initially comet-like mineral mix and their replacement by Fe-bearing pyroxenes, amorphous pyroxene, and silica. We reject models based on vaporizing optical stellar megaflares, aqueous alteration, or giant hypervelocity impacts as unable to produce the observed mineralogy. Scenarios involving unusually high Si abundances are at odds with the normal stellar absorption near-infrared feature strengths for Mg, Fe, and Si. Models involving intense space weathering of a thin surface patina via moderate (T < 1300 K) heating and energetic ion sputtering due to a stellar superflare from the F4V primary are consistent with the observations. The space weathered patina should be reddened, contain copious amounts of nanophase Fe, and should be transient on timescales of decades unless replenished.


title: On the last zero process with applications in corporate bankruptcy

abstract: For a spectrally negative L\'evy process $X$, consider $g_t$, the last time $X$ is below the level zero before time $t\geq 0$. We use a perturbation method for L\'evy processes to derive an It\^o formula for the three-dimensional process $\{(g_t,t, X_t), t\geq 0 \}$ and its infinitesimal generator. Moreover, with $U_t:=t-g_t$, the length of a current positive excursion, we derive a general formula that allows us to calculate a functional of the whole path of $ (U, X)=\{(U_t, X_t),t\geq 0\}$ in terms of the positive and negative excursions of the process $X$. As a corollary, we find the joint Laplace transform of $(U_{\mathbf{e}_q}, X_{\mathbf{e}_q})$, where $\mathbf{e}_q$ is an independent exponential time, and the q-potential measure of the process $(U, X)$. Furthermore, using the results mentioned above, we find a solution to a general optimal stopping problem depending on $(U, X)$ with some applications for corporate bankruptcy. Lastly, we establish a link between the optimal prediction of $g_{\infty}$ and optimal stopping problems in terms of $(U, X)$ as per Baurdoux and Pedraza (2020).


title: Multistage Curvilinear Coordinate Transform Based Document Image Dewarping using a Novel Quality Estimator

abstract: The present work demonstrates a fast and improved technique for dewarping nonlinearly warped document images. The images are first dewarped at the page-level by estimating optimum inverse projections using curvilinear homography. The quality of the process is then estimated by evaluating a set of metrics related to the characteristics of the text lines and rectilinear objects for measuring parallelism, orthogonality, etc. These are designed specifically to estimate the quality of the dewarping process without the need of any ground truth. If the quality is estimated to be unsatisfactory, the page-level dewarping process is repeated with finer approximations. This is followed by a line-level dewarping process that makes granular corrections to the warps in individual text-lines. The methodology has been tested on the CBDAR 2007 / IUPR 2011 document image dewarping dataset and is seen to yield the best OCR accuracy in the shortest amount of time, till date. The usefulness of the methodology has also been evaluated on the DocUNet 2018 dataset with some minor tweaks, and is seen to produce comparable results.


title: Square root of a multivector of Clifford algebras in 3D: A game with signs

abstract: An algorithm to extract the square root in radicals from a multivector (MV) in real Clifford algebras Cl(p,q) for n=p+q <=3 is presented. We show that in the algebras Cl(3,0), Cl(1,2) and Cl(0,3) there are up to four isolated roots in a case of the most general (generic) MV. The algebra Cl(2,1) makes up an exception and the MV here can have up to 16 isolated roots. In addition to isolated roots, a continuum of roots can appear in all algebras except p+q=1. A number of examples are provided to illustrate properties of various roots that may appear in n=3 Clifford algebras.


title: A dark matter telescope probing the 6 to 60 GHz band

abstract: In this article we present the Dark-photons$\&$Axion-Like particles Interferometer (DALI), a novel experiment designed for the detection of photon-mixing cold dark matter in the microwave band between 6 and 60 GHz. DALI is a haloscope for the simultaneous search for axions, ALPs and dark photons, with a number of novelties that make it unique. First, it is a dark matter telescope, with a capacity for pointing, tracking and rastering objects and areas in the sky. This potentially allows one to detect relativistic dark matter particles, substructures and flows, without compromising the simultaneous scanning for dark matter relic particles present in the laboratory. Second, it has been designed using commercial technology. This will allow feasible manufacture at a reasonable cost, thereby mitigating the need for R\&D and facilitating maintenance. Finally, it benefits from a high sensitivity over a broad band of frequencies with only minimal reconfiguration.

categories: physics.ins-det astro-ph.HE astro-ph.IM hep-ex hep-ph


title: Recommending Deployment Strategies for Collaborative Tasks

abstract: Our work contributes to aiding requesters in deploying collaborative tasks in crowdsourcing. We initiate the study of recommending deployment strategies for collaborative tasks to requesters that are consistent with deployment parameters they desire: a lower-bound on the quality of the crowd contribution, an upper-bound on the latency of task completion, and an upper-bound on the cost incurred by paying workers. A deployment strategy is a choice of value for three dimensions: Structure (whether to solicit the workforce sequentially or simultaneously), Organization (to organize it collaboratively or independently), and Style (to rely solely on the crowd or to combine it with machine algorithms). We propose StratRec, an optimization-driven middle layer that recommends deployment strategies and alternative deployment parameters to requesters by accounting for worker availability. Our solutions are grounded in discrete optimization and computational geometry techniques that produce results with theoretical guarantees. We present extensive experiments on Amazon Mechanical Turk and conduct synthetic experiments to validate the qualitative and scalability aspects of StratRec.


title: Convolution invariant linear functionals and applications to summability methods

abstract: We study topologically invariant means on $L^{\infty}(\mathbb{R})$, the set of all essentially bounded functions on the real line, and prove that invariance with respect to a single convolution operator is sufficient for a mean to be topologically invariant. We also consider some applications of this result to summability methods. In particular, the notion of almost convergence is introduced for a function in $L^{\infty}(\mathbb{R})$, and a Tauberian theorem concerning almost convergence and a summability method defined by a Wiener kernel is obtained. Further, for the $C_{\infty}$ summability method, which is defined by the limit of H\"{o}lder summability methods, we provide a necessary and sufficient condition for a given function to be $C_{\infty}$ summable.


title: Guidance and Evaluation: Semantic-Aware Image Inpainting for Mixed Scenes

abstract: Completing a corrupted image with correct structures and reasonable textures for a mixed scene remains an elusive challenge. Since the missing hole in a mixed scene of a corrupted image often contains various semantic information, conventional two-stage approaches utilizing structural information often lead to the problem of unreliable structural prediction and ambiguous image texture generation. In this paper, we propose a Semantic Guidance and Evaluation Network (SGE-Net) to iteratively update the structural priors and the inpainted image in an interplay framework of semantics extraction and image inpainting. It utilizes semantic segmentation map as guidance in each scale of inpainting, under which location-dependent inferences are re-evaluated, and, accordingly, poorly-inferred regions are refined in subsequent scales. Extensive experiments on real-world images of mixed scenes demonstrated the superiority of our proposed method over state-of-the-art approaches, in terms of clear boundaries and photo-realistic textures.


title: Diversity can be Transferred: Output Diversification for White- and Black-box Attacks

abstract: Adversarial attacks often involve random perturbations of the inputs drawn from uniform or Gaussian distributions, e.g., to initialize optimization-based white-box attacks or generate update directions in black-box attacks. These simple perturbations, however, could be sub-optimal as they are agnostic to the model being attacked. To improve the efficiency of these attacks, we propose Output Diversified Sampling (ODS), a novel sampling strategy that attempts to maximize diversity in the target model's outputs among the generated samples. While ODS is a gradient-based strategy, the diversity offered by ODS is transferable and can be helpful for both white-box and black-box attacks via surrogate models. Empirically, we demonstrate that ODS significantly improves the performance of existing white-box and black-box attacks. In particular, ODS reduces the number of queries needed for state-of-the-art black-box attacks on ImageNet by a factor of two.


title: Selecting Voting Locations for Fun and Profit

abstract: While manipulative attacks on elections have been well-studied, only recently has attention turned to attacks that account for geographic information, which are extremely common in the real world. The most well known in the media is gerrymandering, in which district border-lines are changed to increase a party's chance to win, but a different geographical manipulation involves influencing the election by selecting the location of polling places, as many people are not willing to go to any distance to vote. In this paper we initiate the study of this manipulation. We find that while it is easy to manipulate the selection of polling places on the line, it becomes difficult already on the plane or in the case of more than two candidates. Moreover, we show that for more than two candidates the problem is inapproximable. However, we find a few restricted cases on the plane where some algorithms perform well. Finally, we discuss how existing results for standard control actions hold in the geographic setting, consider additional control actions in the geographic setting, and suggest directions for future study.


title: Grammars for Document Spanners

abstract: We propose a new grammar-based language for defining information-extractors from documents (text) that is built upon the well-studied framework of document spanners for extracting structured data from text. While previously studied formalisms for document spanners are mainly based on regular expressions, we use an extension of context-free grammars, called {extraction grammars}, to define the new class of context-free spanners. Extraction grammars are simply context-free grammars extended with variables that capture interval positions of the document, namely spans. While regular expressions are efficient for tokenizing and tagging, context-free grammars are also efficient for capturing structural properties. Indeed, we show that context-free spanners are strictly more expressive than their regular counterparts. We reason about the expressive power of our new class and present a pushdown-automata model that captures it. We show that extraction grammars can be evaluated with polynomial data complexity. Nevertheless, as the degree of the polynomial depends on the query, we present an enumeration algorithm for unambiguous extraction grammars that, after quintic preprocessing, outputs the results sequentially, without repetitions, with a constant delay between every two consecutive ones.


title: Transformative Science from the Lunar Farside: Observations of the Dark Ages and Exoplanetary Systems at Low Radio Frequencies

abstract: The farside of the Moon is a pristine, quiet platform to conduct low radio frequency observations of the early Universe's Dark Ages, as well as space weather and magnetospheres associated with habitable exoplanets. In this paper, the astrophysics associated with NASA-funded concept studies will be described including a lunar-orbiting spacecraft, DAPPER, that will measure the 21 cm global spectrum at redshifts 40-80, and an array of low frequency dipoles on the lunar farside surface, FARSIDE, that would detect exoplanet magnetic fields. DAPPER observations (17-38 MHz), using a single cross-dipole antenna, will measure the amplitude of the 21 cm spectrum to the level required to distinguish the standard {\Lambda}CDM cosmological model from those of additional cooling models possibly produced by exotic physics such as dark matter interactions. FARSIDE has a notional architecture consisting of 128 dipole antennas deployed across a 10 km area by a rover. FARSIDE would image the entire sky each minute in 1400 channels over 0.1-40 MHz. This would enable monitoring of the nearest stellar systems for the radio signatures of coronal mass ejections and energetic particle events, and would also detect the magnetospheres of the nearest candidate habitable exoplanets. In addition, FARSIDE would determine the Dark Ages global 21 cm signal at yet lower frequencies and provide a pathfinder for power spectrum measurements.


title: Fate of priority pharmaceuticals and their main metabolites and transformation products in microalgae-based wastewater treatment systems

abstract: The present study evaluates the removal capacity of two high rate algae ponds (HRAPs) to eliminate 12 pharmaceuticals (PhACs) and 26 of their corresponding main metabolites and transformation products. The efficiency of these ponds, operating with and without primary treatment, was compared in order to study their capacity under the best performance conditions (highest solar irradiance). Concentrations of all the target compounds were determined in both water and biomass samples. Removal rates ranged from moderate (40-60%) to high (>60%) for most of them, with the exception of the psychiatric drugs carbamazepine, the beta-blocking agent metoprolol and its metabolite, metoprolol acid. O-desmethylvenlafaxine, despite its very low biodegradability in conventional wastewater treatment plants, was removed to certain extent (13-39%) Biomass concentrations suggested that bioadsorption/bioaccumulation to microalgae biomass was decisive regarding the elimination of some non-biodegradable compounds such as venlafaxine and its main metabolites. HRAP treatment with and without primary treatment did not yield significant differences in terms of PhACs removal efficiency. The implementation of HRAPs as secondary treatment is a viable alternative to CAS in terms of overall wastewater treatment (including organic micropollutants), with generally higher removal performances and implying a green, low-cost and more sustainable technology.


title: Night-time Scene Parsing with a Large Real Dataset

abstract: Although huge progress has been made on scene analysis in recent years, most existing works assume the input images to be in day-time with good lighting conditions. In this work, we aim to address the night-time scene parsing (NTSP) problem, which has two main challenges: 1) labeled night-time data are scarce, and 2) over- and under-exposures may co-occur in the input night-time images and are not explicitly modeled in existing pipelines. To tackle the scarcity of night-time data, we collect a novel labeled dataset, named {\it NightCity}, of 4,297 real night-time images with ground truth pixel-level semantic annotations. To our knowledge, NightCity is the largest dataset for NTSP. In addition, we also propose an exposure-aware framework to address the NTSP problem through augmenting the segmentation process with explicitly learned exposure features. Extensive experiments show that training on NightCity can significantly improve NTSP performances and that our exposure-aware model outperforms the state-of-the-art methods, yielding top performances on our dataset as well as existing datasets.


title: Collaborative Spectrum Sensing in Tactical Wireless Networks

abstract: In this paper, we propose an algorithm for channel sensing, collaboration, and transmission for networks of Tactical Communication Systems that are facing intrusions from hostile Jammers. Members of the network begin by scanning the spectrum for Jammers, sharing this information with neighboring nodes, and merging their respective sets of observation data into a decision vector containing the believed occupancy of each channel. A user can then use this vector to find a vacant channel for its transmission. We introduce the concept of nodes sharing these vectors with their peers, who then merge them into super-decision vectors, allowing each node to better identify and select transmission channels. We consider fading scenarios that substantially limit the reliability of the users' observations, unless they cooperate to increase the trustworthiness of their sensing data. We propose a pseudo-random channel selection algorithm that strikes a balance between sensing reliability with the number of channels being sensed. Simulation results show that the proposed system improves the network's overall knowledge of the spectrum and the rate of Jammer-free transmissions while limiting added computational complexity to the nodes of the network, despite the Jammers' unpredictable nature.


title: Evaluation of Rounding Functions in Nearest-Neighbor Interpolation

abstract: A novel evaluation study of the most appropriate round function for nearest-neighbor (NN) image interpolation is presented. Evaluated rounding functions are selected among the five rounding rules defined by the Institute of Electrical and Electronics Engineers (IEEE) 754-2008 standard. Both full- and no-reference image quality assessment (IQA) metrics are used to study and evaluate the influence of rounding functions on NN interpolation image quality. The concept of achieved occurrences over targeted occurrences is used to determine the percentage of achieved occurrences based on the number of test images used. Inferential statistical analysis is applied to deduce from a small number of images and draw a conclusion of the behavior of each rounding function on a bigger number of images. Under the normal distribution and at the level of confidence equals to 95%, the maximum and minimum achievable occurrences by each evaluated rounding function are both provided based on the inferential analysis-based experiments.


title: Hamiltonian Simulation Algorithms for Near-Term Quantum Hardware

abstract: The quantum circuit model is the de-facto way of designing quantum algorithms. Yet any level of abstraction away from the underlying hardware incurs overhead. In the era of near-term, noisy, intermediate-scale quantum (NISQ) hardware with severely restricted resources, this overhead may be unjustifiable. In this work, we develop quantum algorithms for Hamiltonian simulation "one level below" the circuit model, exploiting the underlying control over qubit interactions available in principle in most quantum hardware implementations. We then analyse the impact of these techniques under the standard error model where errors occur per gate, and an error model with a constant error rate per unit time. To quantify the benefits of this approach, we apply it to a canonical example: time-dynamics simulation of the 2D spin Fermi-Hubbard model. We derive analytic circuit identities for efficiently synthesising multi-qubit evolutions from two-qubit interactions. Combined with new error bounds for Trotter product formulas tailored to the non-asymptotic regime and a careful analysis of error propagation under the aforementioned per-gate and per-time error models, we improve upon the previous best methods for Hamiltonian simulation by multiple orders of magnitude. By our calculations, for a 5$\mathbf\times$5 Fermi-Hubbard lattice we reduce the circuit depth from 800,160 to 1460 in the per-gate error model, or the circuit-depth-equivalent to 440 in the per-time error model. This brings Hamiltonian simulation, previously beyond reach of current hardware for non-trivial examples, significantly closer to being feasible in the NISQ era.


title: How to Improve AI Tools (by Adding in SE Knowledge): Experiments with the TimeLIME Defect Reduction Tool

abstract: AI algorithms are being used with increased frequency in SE research and practice. Such algorithms are usually commissioned and certified using data from outside the SE domain. Can we assume that such algorithms can be used ''off-the-shelf'' (i.e. with no modifications)? To say that another way, are there special features of SE problems that suggest a different and better way to use AI tools? To answer these questions, this paper reports experiments with TimeLIME, a variant of the LIME explanation algorithm from KDD'16. LIME can offer recommendations on how to change static code attributes in order to reduce the number of defects in the next software release. That version of LIME used an internal weighting tool to decide what attributes to include/exclude in those recommendations. TimeLIME improves on that weighting scheme using the following SE knowledge: software comes in releases; an implausible change to software is something that has never been changed in prior releases; so it is better to use plausible changes, i.e. changes with some precedent in the prior releases. By restricting recommendations to just the frequently changed attributes, TimeLIME can produce (a)~dramatically better explanations of what causes defects and (b)~much better recommendations on how to fix buggy code. Apart from these specific results about defect reduction and TimeLIME, the more general point of this paper is that our community should be more careful about using off-the-shelf AI tools, without first applying SE knowledge. As shown here, it may not be a complex matter to apply that knowledge. Further, once that SE knowledge is applied, this can result in dramatically better systems.


title: Exploring nuclear exotica at the limits

abstract: The study of nuclear limits has been performed and new physical mechanisms and exotic shapes allowing the extension of nuclear landscape beyond the commonly accepted boundaries have been established. The transition from ellipsoidal to toroidal shapes plays a critical role in potential extension of nuclear landscape to hyperheavy nuclei. Rotational excitations leading to the birth of particle (proton or neutron) bound rotational bands provide a mechanism for an extension of nuclear landscape beyond spin zero proton and neutron drip lines.


title: A Deep Exposure in High Resolution X-Rays Reveals the Hottest Plasma in the $\zeta\,$Puppis Wind

abstract: We have obtained a very deep exposure (813 ks) of $\zeta\,$Puppis (O4 supergiant) with the Chandra/HETG Spectrometer. Here we report on analysis of the 1-9 \r{A} region, especially well suited for Chandra, which has a significant contribution from continuum emission between well separated emission lines from high-ionization species. These data allow us to study the hottest plasma present through the continuum shape and emission line strengths. Assuming a powerlaw emission measure distribution which has a high-temperature cut-off, we find that the emission is consistent with a thermal spectrum having a maximum temperature of 12 MK. This implies an effective wind shock velocity of $900\,\mathrm{km\,s^{-1}}$, well below the wind terminal speed of $2250\,\mathrm{km\,s^{-1}}$. For X-ray emission which forms close to the star, the speed and X-ray flux are larger than can be easily reconciled with strictly self-excited line-deshadowing-instability models, suggesting a need for a fraction of the wind to be accelerated extremely rapidly right from the base. This is not so much a dynamical instability as a nonlinear response to changing boundary conditions.


title: Fourier coefficients of $\operatorname{Sp}(4)$ Eisenstein Series

abstract: We compute explicit formulae for the constant terms and Fourier coefficients for Eisenstein series on $\operatorname{Sp}(4,\mathbb{R})$, in terms of zeta functions and Whittaker functions. We also develop a generalisation of Ramanujan sums to $\operatorname{Sp}(4,\mathbb{Z})$, which appears as coefficients in the Fourier expansion for the minimal Eisenstein series.


title: A novel staggered semi-implicit space-time discontinuous Galerkin method for the incompressible Navier-Stokes equations

abstract: A new high order accurate staggered semi-implicit space-time discontinuous Galerkin (DG) method is presented for the simulation of viscous incompressible flows on unstructured triangular grids in two space dimensions. The staggered DG scheme defines the discrete pressure on the primal triangular mesh, while the discrete velocity is defined on a staggered edge-based dual quadrilateral mesh. In this paper, a new pair of equal-order-interpolation velocity-pressure finite elements is proposed. On the primary triangular mesh (the pressure elements) the basis functions are piecewise polynomials of degree $N$ and are allowed to jump on the boundaries of each triangle. On the dual mesh instead (the velocity elements), the basis functions consist in the union of piecewise polynomials of degree $N$ on the two subtriangles that compose each quadrilateral and are allowed to jump only on the dual element boundaries, while they are continuous inside. In other words, the basis functions on the dual mesh are built by continuous finite elements on the subtriangles. This choice allows the construction of an efficient, quadrature-free and memory saving algorithm. In our coupled space-time pressure correction formulation for the incompressible Navier-Stokes equations, arbitrary high order of accuracy in time is achieved through the use of time-dependent test and basis functions, in combination with simple and efficient Picard iterations. Several numerical tests on classical benchmarks confirm that the proposed method outperforms existing staggered semi-implicit space-time DG schemes, not only from a computer memory point of view, but also concerning the computational time.


title: A motivic Segal theorem for open pairs of smooth schemes over an infinite perfect field

abstract: V. Voevodsyky laid the groundwork of delooping motivic spaces in order to provide a new, more computation-friendly, construction of the stable motivic category $SH(k)$, G. Garkusha and I. Panin made that project a reality, while collaborating with A. Ananievsky, A. Neshitov and A. Druzhinin. In particular, G. Garkusha and I. Panin proved that for an infinite perfect field $k$ and any $k$-smooth scheme $X$ the canonical morphism of motivic spaces $C_*Fr(X)\to \Omega^{\infty}_{\mathbb{P}^1} \Sigma^{\infty}_{\mathbb{P}^1} (X_+)$ is Nisnevich-locally a group-completion. In the present work, a generalisation of that theorem to the case of smooth open pairs $(X,U),$ where $X$ is a $k$-smooth scheme, $U$ is its open subscheme intersecting each component of $X$ in a nonempty subscheme. We claim that in this case the motivic space $C_*Fr((X,U))$ is Nisnevich-locally connected, and the motivic space morphism $C_*Fr((X,U))\to \Omega^{\infty}_{\mathbb{P}^1} \Sigma^{\infty}_{\mathbb{P}^1} (X/U)$ is Nisnevich-locally a weak equivalence. Moreover, we show that if the codimension of $S=X-U$ in each component of $X$ is greater than $r \geq 0,$ the simplicial sheaf $C_*Fr((X,U))$ is locally $r$-connected.


title: BARR-C:2018 and MISRA C:2012: Synergy Between the Two Most Widely Used C Coding Standards

abstract: The Barr Group's Embedded C Coding Standard (BARR-C:2018, which originates from the 2009 Netrino's Embedded C Coding Standard) is, for coding standards used by the embedded system industry, second only in popularity to MISRA C. However, the choice between MISRA C:2012 and BARR-C:2018 needs not be a hard decision since they are complementary in two quite different ways. On the one hand, BARR-C:2018 has removed all the incompatibilities with respect to MISRA C:2012 that were present in the previous edition (BARR-C:2013). As a result, disregarding programming style, BARR-C:2018 defines a subset of C that, while preventing a significant number of programming errors, is larger than the one defined by MISRA C:2012. On the other hand, concerning programming style, whereas MISRA C leaves this to individual organizations, BARR-C:2018 defines a programming style aimed primarily at minimizing programming errors. As a result, BARR-C:2018 can be seen as a first, dramatically useful step to C language subsetting that is suitable for all kinds of projects; critical projects can then evolve toward MISRA C:2012 compliance smoothly while maintaining the BARR-C programming style. In this paper, we introduce BARR-C:2018, we describe its relationship with MISRA C:2012, and we discuss the parallel and serial adoption of the two coding standards.


title: Exploring Gaussian mixture model framework for speaker adaptation of deep neural network acoustic models

abstract: In this paper we investigate the GMM-derived (GMMD) features for adaptation of deep neural network (DNN) acoustic models. The adaptation of the DNN trained on GMMD features is done through the maximum a posteriori (MAP) adaptation of the auxiliary GMM model used for GMMD feature extraction. We explore fusion of the adapted GMMD features with conventional features, such as bottleneck and MFCC features, in two different neural network architectures: DNN and time-delay neural network (TDNN). We analyze and compare different types of adaptation techniques such as i-vectors and feature-space adaptation techniques based on maximum likelihood linear regression (fMLLR) with the proposed adaptation approach, and explore their complementarity using various types of fusion such as feature level, posterior level, lattice level and others in order to discover the best possible way of combination. Experimental results on the TED-LIUM corpus show that the proposed adaptation technique can be effectively integrated into DNN and TDNN setups at different levels and provide additional gain in recognition performance: up to 6% of relative word error rate reduction (WERR) over the strong feature-space adaptation techniques based on maximum likelihood linear regression (fMLLR) speaker adapted DNN baseline, and up to 18% of relative WERR in comparison with a speaker independent (SI) DNN baseline model, trained on conventional features. For TDNN models the proposed approach achieves up to 26% of relative WERR in comparison with a SI baseline, and up 13% in comparison with the model adapted by using i-vectors. The analysis of the adapted GMMD features from various points of view demonstrates their effectiveness at different levels.


title: Measuring the tangle of three-qubit states

abstract: We present a quantum circuit that transforms an unknown three-qubit state into its canonical form, up to relative phases, given many copies of the original state. The circuit is made of three single-qubit parametrized quantum gates, and the optimal values for the parameters are learned in a variational fashion. Once this transformation is achieved, direct measurement of outcome probabilities in the computational basis provides an estimate of the tangle, which quantifies genuine tripartite entanglement. We perform simulations on a set of random states under different noise conditions to asses the validity of the method.


title: MACHO 311.37557.169: A VY Scl star

abstract: Optical surveys, such as the MACHO project, often uncover variable stars whose classification requires followup observations by other instruments. We performed X-ray spectroscopy and photometry of the unusual variable star MACHO 311.37557.169 with \XMM\ in April 2018, supplemented by archival X-ray and optical spectrographic data. The star has a bolometric X-ray luminosity of about $1\times 10^{32}$ erg s$^{-1}$ cm$^{-2}$ and a heavily absorbed two-temperature plasma spectrum. The shape of its light curve, its overall brightness, its X-ray spectrum, and the emission lines in its optical spectrum suggest that it is most likely a VY~Scl cataclysmic variable.


title: Relative Singularity categories and singular equivalences

abstract: Let $R$ be a right notherian ring. We introduce the concept of relative singularity category $\Delta_{\mathcal{X}}(R)$ of $R$ with respect to a contravariantly finite subcategory $\mathcal{X}$ of $\rm{mod}\mbox{-}R.$ Along with some finiteness conditions on $\mathcal{X}$, we prove that $\Delta_{\mathcal{X}}(R)$ is triangle equivalent to a subcategory of the homotopy category $\mathbb{K}_{\rm{ac}}(\mathcal{X})$ of exact complexes over $\mathcal{X}$. As an application, a new description of the classical singularity category $\mathbb{D}_{\rm{sg}}(R)$ is given. The relative singularity categories are applied to lift a stable equivalence between two suitable subcategories of the module categories of two given right notherian ring to get a singular equivalence between the rings. In different types of rings, including path rings, triangular matrix rings, trivial extension rings and tensor rings, we provide some consequences for their singularity categories.


title: Provably Efficient Exploration for Reinforcement Learning Using Unsupervised Learning

abstract: Motivated by the prevailing paradigm of using unsupervised learning for efficient exploration in reinforcement learning (RL) problems [tang2017exploration,bellemare2016unifying], we investigate when this paradigm is provably efficient. We study episodic Markov decision processes with rich observations generated from a small number of latent states. We present a general algorithmic framework that is built upon two components: an unsupervised learning algorithm and a no-regret tabular RL algorithm. Theoretically, we prove that as long as the unsupervised learning algorithm enjoys a polynomial sample complexity guarantee, we can find a near-optimal policy with sample complexity polynomial in the number of latent states, which is significantly smaller than the number of observations. Empirically, we instantiate our framework on a class of hard exploration problems to demonstrate the practicality of our theory.


title: Adversarial Encoder-Multi-Task-Decoder for Multi-Stage Processes

abstract: In multi-stage processes, decisions occur in an ordered sequence of stages. Early stages usually have more observations with general information (easier/cheaper to collect), while later stages have fewer observations but more specific data. This situation can be represented by a dual funnel structure, in which the sample size decreases from one stage to the other while the information increases. Training classifiers in this scenario is challenging since information in the early stages may not contain distinct patterns to learn (underfitting). In contrast, the small sample size in later stages can cause overfitting. We address both cases by introducing a framework that combines adversarial autoencoders (AAE), multi-task learning (MTL), and multi-label semi-supervised learning (MLSSL). We improve the decoder of the AAE with an MTL component so it can jointly reconstruct the original input and use feature nets to predict the features for the next stages. We also introduce a sequence constraint in the output of an MLSSL classifier to guarantee the sequential pattern in the predictions. Using real-world data from different domains (selection process, medical diagnosis), we show that our approach outperforms other state-of-the-art methods.


title: Vacancy diffusion in multi-principal element alloys: the role of chemical disorder in the ordered lattice

abstract: Many of the purported virtues of Multi-Principal Element Alloys (MPEAs), such as corrosion, high-temperature oxidation and irradiation resistance, are highly sensitive to vacancy diffusivity. Similarly, solute interdiffusion is governed by vacancy diffusion -- it is often unclear whether MPEAs are truly stable, or effectively stabilized by slow interdiffusion. The considerable composition space afforded to these alloys makes optimizing for desired properties a daunting task; theoretical and computational tools are necessary to guide alloy development. For diffusion, such tools depend on both a knowledge of the vacancy migration barriers within a given alloy and an understanding of how these barriers influence vacancy diffusivity. We present a generalized theory of vacancy diffusion in rugged energy landscapes, paired with Kinetic Monte Carlo simulations of MPEA vacancy diffusion. The barrier energy statistics are informed by nudged elastic band calculations in the equiatomic CoNiCrFeMn alloy. Theory and simulations show that vacancy diffusion in solid-solution MPEAs is not necessarily sluggish, but can potentially be tuned, and that trap models are an insufficient explanation for sluggish diffusion in the CoNiCrFeMn HEA. These results also show that any model that endeavors to faithfully represent diffusion-related phenomena must account for the full nature of the energy landscape, not just the migration barriers.


title: Solving constrained optimization problems without Lagrange multipliers

abstract: Constrained optimization problems exist in many domains of science, such as thermodynamics, mechanics, economics, etc. These problems are classically solved with the help of the Lagrange multipliers and the Lagrangian function. However, the disadvantage of this approach is that it artificially increases the dimensionality of the problem. Here, we show that the determinant of the Jacobian of the problem (function to optimize and constraints) is null. This extra equation transforms any equality-constrained optimization problem into a solving problem of same dimension. We also introduced the constraint matrices as the largest square submatrices of the Jacobian of the constraints. The boundaries of the constraint domain are given by the nullity of their determinants. The constraint matrices also permit to write the function to be optimized as a Taylor series of any of its variable, with its coefficients algebraically determined by an iterative process of partial derivation.


title: GENIEx: A Generalized Approach to Emulating Non-Ideality in Memristive Xbars using Neural Networks

abstract: The analog nature of computing in Memristive crossbars poses significant issues due to various non-idealities such as: parasitic resistances, non-linear I-V characteristics of the device etc. The non-idealities can have a detrimental impact on the functionality i.e. computational accuracy of crossbars. Past works have explored modeling the non-idealities using analytical techniques. However, several non-idealities have data dependent behavior. This can not be captured using analytical (non data-dependent) models thereby, limiting their suitability in predicting application accuracy. To address this, we propose a Generalized Approach to Emulating Non-Ideality in Memristive Crossbars using Neural Networks (GENIEx), which accurately captures the data-dependent nature of non-idealities. We perform extensive HSPICE simulations of crossbars with different voltage and conductance combinations. Following that, we train a neural network to learn the transfer characteristics of the non-ideal crossbar. Next, we build a functional simulator which includes key architectural facets such as \textit{tiling}, and \textit{bit-slicing} to analyze the impact of non-idealities on the classification accuracy of large-scale neural networks. We show that GENIEx achieves \textit{low} root mean square errors (RMSE) of $0.25$ and $0.7$ for low and high voltages, respectively, compared to HSPICE. Additionally, the GENIEx errors are $7\times$ and $12.8\times$ better than an analytical model which can only capture the linear non-idealities. Further, using the functional simulator and GENIEx, we demonstrate that an analytical model can overestimate the degradation in classification accuracy by $\ge 10\%$ on CIFAR-100 and $3.7\%$ on ImageNet datasets compared to GENIEx.


title: Old Problems, Classical Methods, New Solutions

abstract: We use a powerful extension of the classical method of heat potentials, recently developed by the present author and his collaborators, to solve several significant problems of financial mathematics. We consider the following problems in detail: (A) calibrating the default boundary in the structural default framework to a constant default intensity; (B) calculating default probability for a representative bank in the mean-field framework; (C) finding the hitting time probability density of an Ornstein-Uhlenbeck process. Several other problems, including pricing American put options and finding optimal mean-reverting trading strategies, are mentioned in passing. Besides, two non-financial applications -- the supercooled Stefan problem and the integrate-and-fire neuroscience problem -- are briefly discussed as well.


title: Probing the interactions between interstitial hydrogen atoms in niobium through density functional theory calculations

abstract: Past experiments about hydrogen absorption in niobium have revealed specific properties about interactions between interstitial hydrogen atoms. It has been reported that there are long-range attractive and short-range repulsive interactions between interstitial hydrogen atoms in niobium. It has also been reported that these interactions are of many-body nature. While previous understanding of these interactions is based on experimental inferences from past experiments, through these calculations, for the first time, we can understand the nature of the interactions at a fundamental level. In this work, we use Density Functional Theory calculations to study the interactions of interstitial hydrogen atoms in niobium. We report here that these interactions are a combination of an attractive, indirect image interaction and a repulsive, direct interaction. Through our calculations, we also infer here that these interactions indeed have many-body characteristics.


title: Bosonization based on Clifford algebras and its gauge theoretic interpretation

abstract: We study the properties of a bosonization procedure based on Clifford algebra valued degrees of freedom, valid for spaces of any dimension. We present its interpretation in terms of fermions in presence of $\mathbb{Z}_2$ gauge fields satisfying a modified Gauss' law, resembling Chern-Simons-like theories. Our bosonization prescription involves constraints, which are interpreted as a flatness condition for the gauge field. Solution of the constraints is presented for toroidal geometries of dimension two. Duality between our model and $(d-1)$-form $\mathbb{Z}_2$ gauge theory is derived, which elucidates the relation between the approach taken here with another bosonization map proposed recently.

categories: math-ph cond-mat.str-el hep-lat hep-th math.MP


title: Model-based Reinforcement Learning for Decentralized Multiagent Rendezvous

abstract: Collaboration requires agents to align their goals on the fly. Underlying the human ability to align goals with other agents is their ability to predict the intentions of others and actively update their own plans. We propose hierarchical predictive planning (HPP), a model-based reinforcement learning method for decentralized multiagent rendezvous. Starting with pretrained, single-agent point to point navigation policies and using noisy, high-dimensional sensor inputs like lidar, we first learn via self-supervision motion predictions of all agents on the team. Next, HPP uses the prediction models to propose and evaluate navigation subgoals for completing the rendezvous task without explicit communication among agents. We evaluate HPP in a suite of unseen environments, with increasing complexity and numbers of obstacles. We show that HPP outperforms alternative reinforcement learning, path planning, and heuristic-based baselines on challenging, unseen environments. Experiments in the real world demonstrate successful transfer of the prediction models from sim to real world without any additional fine-tuning. Altogether, HPP removes the need for a centralized operator in multiagent systems by combining model-based RL and inference methods, enabling agents to dynamically align plans.


title: Homogeneous metric ANR-compacta

abstract: This is a survey of most important results and unsolved problems about homogeneous finite-dimensional metric $ANR$-compacta. We also discuss some partial results and possible ways of solutions.


title: Will the real Hardy-Ramanujan formula please stand up?

abstract: The Hardy-Ramanujan formula for the number of integer partitions of $n$ is one of the most popular results in partition theory. While the unabridged final formula has been celebrated as reflecting the genius of its authors, it has become all too common to attribute either some simplified version of the formula which is not as ingenious, or an alternative more elegant version which was expanded on afterwards by other authors. We attempt to provide a clear and compelling justification for distinguishing between the various formulas and simplifications, with a summarizing list of key take-aways in the final section.


title: Resummation of QED radiative corrections in a strong constant crossed field

abstract: By considering radiative corrections of up to 3rd-loop order, Ritus and Narozhny conjectured that the proper expansion parameter for QED in a strong constant crossed field is $g=\alpha\chi^{2/3}$, where the dynamical quantum parameter $\chi=e\sqrt{-(Fp)^2}/m^3$ combines the particle momentum $p$ with the external field strength tensor $F$. Here we present and discuss the first non-perturbative result in this context, the resummed bubble-type polarization corrections to the electron self-energy in a constant crossed field. Our analysis confirms the relevance of the scaling parameter $g$ to the enhancement of bubble-type radiative corrections. This parameter actually represents the characteristic value of the ratio of the 1-loop polarization bubble to the photon virtuality. After an all-order resummation we identify and discuss two contributions to the self-energy with different formation regions and asymptotic behavior for $g\gg1$. Whereas the breakdown of perturbation theory occurs already for $g\gtrsim1$, the leading-order result remains dominant until the asymptotic regime $g\gg 1$ is reached. However, the latter is specific to processes like elastic scattering or photon emission and does not have to remain true for general higher-order QED processes.


title: Numerical analysis for a system coupling curve evolution attached orthogonally to a fixed boundary, to a reaction-diffusion equation on the curve

abstract: We consider a semi-discrete finite element approximation for a system consisting of the evolution of a planar curve evolving by forced curve shortening flow inside a given bounded domain $\Omega \subset \mathbb{R}^2$, such that the curve meets the boundary $\partial \Omega$ orthogonally, and the forcing is a function of the solution of a reaction-diffusion equation that holds on the evolving curve. We prove optimal error bounds for the resulting approximation and present numerical experiments.


title: Assessing the reliability and validity ranges of magnetic characterization methods

abstract: Evaluating the interaction fields of magnetic nanowires (MNWs) is of utmost importance for advancing their functionality in diverse applications including spintronic devices and nanomedicine. In recent years, several quantitative methods have been proposed and become inevitable tools to quantify the interaction fields and decouple their effects from the coercivity. However, the uncertainty of the attained results arose countless open questions leading to discrepancies among the literature. Here, we employ our novel experimental method, named the projection method, to resolve these discrepancies. Using a comparative analysis of the four most commonly used methods (hysteresis loops method, remanence curves method, first-order reversal curve method, and projection method), we unambiguously explicate the reliability and validity limits of these methods to elucidate the origin of the discrepancies. We show that the remanence curves method must solely be used for quantifying the interaction fields if they are considerably weaker compared to the coercivity. Furthermore, we show that both the remanence curves method and first-order reversal curve method fail to fully decouple the interaction fields effects from the coercivity, similar to the hysteresis loops method.


title: On unsteady flows of pore pressure-activated granular materials

abstract: We investigate mathematical properties of the system of nonlinear partial differential equations that describe, under certain simplifying assumptions, evolutionary processes in water-saturated granular materials. The unconsolidated solid matrix behaves as an ideal plastic material before the activation takes place and then it starts to flow as a Newtonian or a generalized Newtonian fluid. The plastic yield stress is non-constant and depends on the difference between the given lithostatic pressure and the pressure of the fluid in a pore space. We study unsteady three-dimensional flows in an impermeable container, subject to stick-slip boundary conditions. Under realistic assumptions on the data, we establish long-time and large-data existence theory.


title: Existence, Uniqueness and Asymptotic Behavior of Regular Time-Periodic Viscous Flow around a Moving Body

abstract: We show existence and uniqueness of regular time-periodic solutions to the Navier-Stokes problem in the exterior of a rigid body, $\mathscr B$, that moves by arbitrary (sufficiently smooth) time-periodic translational motion of the same period, provided the size of the data is suitably restricted. Moreover, we characterize the spatial asymptotic behavior of such solutions and prove, in particular, that if $\mathscr B$ has a nonzero net motion identified by a constant velocity $\bar{\mathbf{\xi}}$ (say), then the solution exhibit a wake-like behavior in the direction $-\bar{\mathbf{\xi}}$ entirely analogous to that of a steady-state flow around a body that moves with velocity $\bar{\mathbf{\xi}}$.


title: The inverse problem for Hamilton-Jacobi equations and semiconcave envelopes

abstract: We study the inverse problem, or inverse design problem, for a time-evolution Hamilton-Jacobi equation. More precisely, given a target function $u_T$ and a time horizon $T>0$, we aim to construct all the initial conditions for which the viscosity solution coincides with $u_T$ at time $T$. As it is common in this kind of nonlinear equations, the target might not be reachable. We first study the existence of at least one initial condition leading the system to the given target. The natural candidate, which indeed allows determining the reachability of $u_T$, is the one obtained by reversing the direction of time in the equation, considering $u_T$ as terminal condition. In this case, we use the notion of backward viscosity solution, that provides existence and uniqueness for the terminal-value problem. We also give an equivalent reachability condition based on a differential inequality, that relates the reachability of the target with its semiconcavity properties. Then, for the case when $u_T$ is reachable, we construct the set of all initial conditions for which the solution coincides with $u_T$ at time $T$. Note that in general, such initial conditions are not unique. Finally, for the case when the target $u_T$ is not necessarily reachable, we study the projection of $u_T$ on the set of reachable targets, obtained by solving the problem backward and then forward in time. This projection is then identified with the solution of a fully nonlinear obstacle problem, and can be interpreted as the semiconcave envelope of $u_T$, i.e. the smallest reachable target bounded from below by $u_T$.


title: Finite-Element Formulation for Advection-Reaction Equations with Change of Variable and Discontinuity Capturing

abstract: We propose a change of variable approach and discontinuity capturing methods to ensure physical constraints for advection-reaction equations discretized by the finite element method. This change of variable confines the concentration below an upper bound in a very natural way. For the non-negativity constraint, we propose to use a discontinuity capturing method defined on the reference element that is combined with an anisotropic crosswind-dissipation operator. This discontinuity capturing cannot completely eliminate negative values but effectively minimizes their occurrence. The proposed methods are applied to different biophysical models and show a good agreement with experimental results for the FDA benchmark blood pump for a physiological red blood cell pore formation model.

categories: math.NA cs.NA physics.bio-ph physics.flu-dyn


title: Complementary symmetric Rote sequences: the critical exponent and the recurrence function

abstract: We determine the critical exponent and the recurrence function of complementary symmetric Rote sequences. The formulae are expressed in terms of the continued fraction expansions associated with the S-adic representations of the corresponding standard Sturmian sequences. The results are based on a thorough study of return words to bispecial factors of Sturmian sequences. Using the formula for the critical exponent, we describe all complementary symmetric Rote sequences with the critical exponent less than or equal to 3, and we show that there are uncountably many complementary symmetric Rote sequences with the critical exponent less than the critical exponent of the Fibonacci sequence. Our study is motivated by a~conjecture on sequences rich in palindromes formulated by Baranwal and Shallit. Its recent solution by Curie, Mol, and Rampersad uses two particular complementary symmetric Rote sequences.


title: End-to-End Velocity Estimation For Autonomous Racing

abstract: Velocity estimation plays a central role in driverless vehicles, but standard and affordable methods struggle to cope with extreme scenarios like aggressive maneuvers due to the presence of high sideslip. To solve this, autonomous race cars are usually equipped with expensive external velocity sensors. In this paper, we present an end-to-end recurrent neural network that takes available raw sensors as input (IMU, wheel odometry, and motor currents) and outputs velocity estimates. The results are compared to two state-of-the-art Kalman filters, which respectively include and exclude expensive velocity sensors. All methods have been extensively tested on a formula student driverless race car with very high sideslip (10{\deg} at the rear axle) and slip ratio (~20%), operating close to the limits of handling. The proposed network is able to estimate lateral velocity up to 15x better than the Kalman filter with the equivalent sensor input and matches (0.06 m/s RMSE) the Kalman filter with the expensive velocity sensor setup.


title: Squeeze free space with nonlocal flat optics

abstract: There has been substantial interest in miniaturizing optical systems by flat optics. However, one essential optical component, free space, fundamentally cannot be substituted with conventional local flat optics with space-dependent transfer functions, since the transfer function of free space is momentum-dependent instead. Overcoming this difficulty is important to achieve utmost miniaturization of optical systems. In this work, we show that free space can be substituted with nonlocal flat optics operating directly in the momentum domain. We derive the general criteria for an optical device to replace free space, and provide a concrete design of a photonic crystal slab device. Such a device can substitute much thicker free space with a compression ratio as high as 144. Our work paves the way for utmost miniaturization of optical systems using combination of local and nonlocal flat optics.


title: Operators on positive semidefinite inner product spaces

abstract: We give canonical forms of selfadjoint and isometric operators on a complex vector space $U$ with scalar product given by a positive semidefinite Hermitian form, and of Hermitian forms on $U$. For an arbitrary system of semiunitary spaces and linear mappings on/between them, we give an algorithm that reduces their matrices to canonical form.


title: Getting Fairness Right: Towards a Toolbox for Practitioners

abstract: The potential risk of AI systems unintentionally embedding and reproducing bias has attracted the attention of machine learning practitioners and society at large. As policy makers are willing to set the standards of algorithms and AI techniques, the issue on how to refine existing regulation, in order to enforce that decisions made by automated systems are fair and non-discriminatory, is again critical. Meanwhile, researchers have demonstrated that the various existing metrics for fairness are statistically mutually exclusive and the right choice mostly depends on the use case and the definition of fairness. Recognizing that the solutions for implementing fair AI are not purely mathematical but require the commitments of the stakeholders to define the desired nature of fairness, this paper proposes to draft a toolbox which helps practitioners to ensure fair AI practices. Based on the nature of the application and the available training data, but also on legal requirements and ethical, philosophical and cultural dimensions, the toolbox aims to identify the most appropriate fairness objective. This approach attempts to structure the complex landscape of fairness metrics and, therefore, makes the different available options more accessible to non-technical people. In the proven absence of a silver bullet solution for fair AI, this toolbox intends to produce the fairest AI systems possible with respect to their local context.


title: Projected Gross-Pitaevskii equation for ring-shaped Bose-Einstein condensates

abstract: We propose an alternative implementation of the Projected Gross-Pitaevskki equation adapted for numerical modeling of the atomic Bose-Einstein condensate trapped in a toroidally-shaped potential. We present an accurate and efficient scheme to evaluate the required matrix elements and calculate time evolution of the matter wave field. We analyze the stability and accuracy of the developed method for equilibrium and nonequilibrium solutions in a ring-shaped trap with additional barrier potential corresponding to recent experimental realizations.


title: Willmore spheres in the 3-sphere revisited

abstract: Bryant \cite{Bryant84} classified all Willmore spheres in $3$-space to be given by minimal surfaces in $\mathbb R^3$ with embedded planar ends. This note provides new explicit formulas for genus 0 minimal surfaces in $\mathbb R^3$ with $2k+1$ embedded planar ends for all $k\geq4.$ Peng and Xiao claimed these examples to exist in \cite{PengXiao2000}, but in the same paper they also claimed the existence of a minimal surface with 7 embedded planar ends, which was falsified by Bryant \cite{Bryant88}.


title: RCNet: Incorporating Structural Information into Deep RNN for MIMO-OFDM Symbol Detection with Limited Training

abstract: In this paper, we investigate learning-based MIMO-OFDM symbol detection strategies focusing on a special recurrent neural network (RNN) -- reservoir computing (RC). We first introduce the Time-Frequency RC to take advantage of the structural information inherent in OFDM signals. Using the time domain RC and the time-frequency RC as the building blocks, we provide two extensions of the shallow RC to RCNet: 1) Stacking multiple time domain RCs; 2) Stacking multiple time-frequency RCs into a deep structure. The combination of RNN dynamics, the time-frequency structure of MIMO-OFDM signals, and the deep network enables RCNet to handle the interference and nonlinear distortion of MIMO-OFDM signals to outperform existing methods. Unlike most existing NN-based detection strategies, RCNet is also shown to provide a good generalization performance even with a limited training set (i.e, similar amount of reference signals/training as standard model-based approaches). Numerical experiments demonstrate that the introduced RCNet can offer a faster learning convergence and as much as 20% gain in bit error rate over a shallow RC structure by compensating for the nonlinear distortion of the MIMO-OFDM signal, such as due to power amplifier compression in the transmitter or due to finite quantization resolution in the receiver.


title: On the spatial and temporal shift in the archetypal seasonal temperature cycle as driven by annual and semi-annual harmonics

abstract: Statistical methods are required to evaluate and quantify the uncertainty in environmental processes, such as land and sea surface temperature, in a changing climate. Typically, annual harmonics are used to characterize the variation in the seasonal temperature cycle. However, an often overlooked feature of the climate seasonal cycle is the semi-annual harmonic, which can account for a significant portion of the variance of the seasonal cycle and varies in amplitude and phase across space. Together, the spatial variation in the annual and semi-annual harmonics can play an important role in driving processes that are tied to seasonality (e.g., ecological and agricultural processes). We propose a multivariate spatio-temporal model to quantify the spatial and temporal change in minimum and maximum temperature seasonal cycles as a function of the annual and semi-annual harmonics. Our approach captures spatial dependence, temporal dynamics, and multivariate dependence of these harmonics through spatially and temporally-varying coefficients. We apply the model to minimum and maximum temperature over North American for the years 1979 to 2018. Formal model inference within the Bayesian paradigm enables the identification of regions experiencing significant changes in minimum and maximum temperature seasonal cycles due to the relative effects of changes in the two harmonics.


title: A comparison of UV and optical metallicities in star-forming galaxies

abstract: Our ability to study the properties of the interstellar medium (ISM) in the earliest galaxies will rely on emission line diagnostics at rest-frame ultraviolet (UV) wavelengths. In this work, we identify metallicity-sensitive diagnostics using UV emission lines. We compare UV-derived metallicities with standard, well-established optical metallicities using a sample of galaxies with rest-frame UV and optical spectroscopy. We find that the He2-O3C3 diagnostic (He II 1640 / C III 1906,1909 vs. O III 1666 / C III 1906,1909) is a reliable metallicity tracer, particularly at low metallicity (12+log(O/H) < 8), where stellar contributions are minimal. We find that the Si3-O3C3 diagnostic (Si III 1883 / C III 1906,1909 vs. O III 1666 / C III 1906,1909) is a reliable metallicity tracer, though with large scatter (0.2-0.3 dex), which we suggest is driven by variations in gas-phase abundances. We find that the C4-O3C3 diagnostic (C IV 1548,1550 / O III 1666 vs. O III 1666 / C III 1906,1909) correlates poorly with optically-derived metallicities. We discuss possible explanations for these discrepant metallicity determinations, including the hardness of the ionizing spectrum, contribution from stellar wind emission, and non-solar-scaled gas-phase abundances. Finally, we provide two new UV oxygen abundance diagnostics, calculated from polynomial fits to the model grid surface in the He2-O3C3 and Si3-O3C3 diagrams.


title: Stochastic gradient descent with random learning rate

abstract: We propose to optimize neural networks with a uniformly-distributed random learning rate. The associated stochastic gradient descent algorithm can be approximated by continuous stochastic equations and analyzed within the Fokker-Planck formalism. In the small learning rate regime, the training process is characterized by an effective temperature which depends on the average learning rate, the mini-batch size and the momentum of the optimization algorithm. By comparing the random learning rate protocol with cyclic and constant protocols, we suggest that the random choice is generically the best strategy in the small learning rate regime, yielding better regularization without extra computational cost. We provide supporting evidence through experiments on both shallow, fully-connected and deep, convolutional neural networks for image classification on the MNIST and CIFAR10 datasets.


title: Ultrasound imaging with three dimensional full-wave nonlinear acoustic simulations. Part 2: sources of image degradation in intercostal imaging

abstract: Fullwave simulations are applied to an intercostal imaging scenario to determine the sources of fundamental and harmonic image degradation with respect to aberration and reverberation. These simulations are based on Part I of this two part paper, which established the Fullwave simulation methods to generate realistic ultrasound images based directly on the first principles of wave propagation in the human body. The ultasound images are generated based on the first principles of propagation and reflection and they describe interplay between distributed aberration and reverberation clutter. Three imaging scenarios that would not be realizable in vivo are investigated in silico. First, the ribs were completely removed and replaced with fat. Then, the ribs were maintained in their anatomically correct configuration to yield a reference image. Finally the ribs were placed closer together in elevation. The propagation based B-mode images show that of these three scenarios the second, anatomically correct configuration, has the best contrast-to-noise ratio. This is due to two competing effects. First the ribs effectively apodize the fundamental and harmonic beams by 3-5 dB. This effect alone would predict an improvement in image quality. However, the B-mode image quality, measured by the contrast-to-noise ratio degrades by 8%. To explain these changes, it is shown that a second effect, multiple reverberation, must be taken into account. A point spread function analysis shows that when the ribs are placed closer together they generate significantly more reverberation clutter (by 2.4 to 2.9 dB), degrading the image quality even though the beamplot has lower sidelobes. In this intercostal imaging scenario the effects of the ribs on beam shape and reverberation are therefore in competition in terms of image quality and there is an optimal acoustic window that balances them out.


title: Low-dimensional approximations of high-dimensional asset price models

abstract: We consider high-dimensional asset price models that are reduced in their dimension in order to reduce the complexity of the problem or the effect of the curse of dimensionality in the context of option pricing. We apply model order reduction (MOR) to obtain a reduced system. MOR has been previously studied for asymptotically stable controlled stochastic systems with zero initial conditions. However, stochastic differential equations modeling price processes are uncontrolled, have non-zero initial states and are often unstable. Therefore, we extend MOR schemes and combine ideas of techniques known for deterministic systems. This leads to a method providing a good pathwise approximation. After explaining the reduction procedure, the error of the approximation is analyzed and the performance of the algorithm is shown conducting several numerical experiments. Within the numerics section, the benefit of the algorithm in the context of option pricing is pointed out.


title: Asymptotic behaviour of Kac polynomials

abstract: We conjecture a formula supported by computations for the valuation of Kac polynomials of a quiver, which only depends on the number of loops at each vertex. We prove a convergence property of renormalized Kac polynomials of quivers when increasing the number of arrows: they converge in the ring of power series, with a linear rate of convergence. Then, we propose a conjecture concerning the global behaviour of the coefficients of Kac polynomials. All computations were made using SageMath.


title: Upconversion mid-infrared dual-comb spectroscopy

abstract: We devise a new detection technique for mid-infrared multi-heterodyne spectroscopy. As an experimental implementation, mid-infrared light interrogates a gas sample in the 3-$\mu$m region of the fundamental CH, OH, NH stretch in molecules and detection is performed in the near-infrared telecommunication region. Spectra showing 18000 resolved comb lines of 100-MHz spacing are recorded at a sensitivity close to the shot-noise limit. The demonstration adds to the unique advantages of dual-comb spectroscopy, such as potential very high resolution, negligible instrumental line width, and direct calibration of the frequency scale with an atomic clock. We also show that direct mid-infrared detection with a differential photo-detector module is beneficial in the limit of low mid-infrared average power.


title: On continuity of local epsilon factors of $\ell$-adic sheaves

abstract: Let $S$ be a noetherian scheme and $f\colon X\to S$ be a smooth morphism of relative dimension 1. For a locally constant sheaf on the complement of a divisor in $X$ at over $S$, Deligne and Laumon proved that the universal local acyclicity is equivalent to the local constancy of Swan conductors. In this article, assuming the universal local acyclicity, we show an analogous result of the continuity of local epsilon factors. We also give a generalization of this result to a family of isolated singularities.


title: Self-Constructing Graph Convolutional Networks for Semantic Labeling

abstract: Graph Neural Networks (GNNs) have received increasing attention in many fields. However, due to the lack of prior graphs, their use for semantic labeling has been limited. Here, we propose a novel architecture called the Self-Constructing Graph (SCG), which makes use of learnable latent variables to generate embeddings and to self-construct the underlying graphs directly from the input features without relying on manually built prior knowledge graphs. SCG can automatically obtain optimized non-local context graphs from complex-shaped objects in aerial imagery. We optimize SCG via an adaptive diagonal enhancement method and a variational lower bound that consists of a customized graph reconstruction term and a Kullback-Leibler divergence regularization term. We demonstrate the effectiveness and flexibility of the proposed SCG on the publicly available ISPRS Vaihingen dataset and our model SCG-Net achieves competitive results in terms of F1-score with much fewer parameters and at a lower computational cost compared to related pure-CNN based work. Our code will be made public soon.


title: Understanding Epidemic Data and Statistics: A case study of COVID-19

abstract: The 2019-Novel-Coronavirus (COVID-19) has affected 116 countries (By March 12) and out of more than 118,000 confirmed cases. Understanding the transmission dynamics of the infection in each country which affected on a daily basis and evaluating the effectiveness of control policies is critical for our further actions. To date, the statistics of COVID-19 reported cases show more than 80 percent of infected had a mild case of disease, while around 14 percent of infected experienced a severe one and about 5 percent are categorized as critical disease victims. Today's report (2020-03-12; daily updates in the prepared website) shows the confirmed cases of COVID-19 in China, South Korea, Italy, and Iran are 80932, 7869, 12462 and 10075; respectively. Calculating the total Case Fatality Rate (CFR) of Italy (2020-03-04), about 7.9% of confirmed cases passed away. Compared to South Korea's rate of 0.76% (10 times lower than Italy) and China's 3.8% (50% lower than Italy), the CFR of Italy is too high. There are some effective policies that yield significant changes in the trend of cases. The lockdown policy in China and Italy (the effect observed after 11 days), Shutdown of all nonessential companies in Hubei (the effect observed after 5 days), combined policy in South Korea and reducing working hours in Iran.


title: Ultrasound imaging of the human body with three dimensional full-wave nonlinear acoustics. Part 1: simulations methods

abstract: Simulations of three dimensional ultrasound propagation in heterogeneous media are computationally intensive due to the constraints arising from the large size of the domain, which is on the order of hundreds of wavelengths, and the small size of scatterers, which are much smaller than a wavelength. Consequently, three dimensional ultrasound imaging simulations are currently based on models that simplify the propagation physics. Here the full three dimensional wave physics is simulated with finite differences to generate ultrasound images of the human body based directly on the first principles of propagation and backscattering. The Visible Human project, a 3D data set of the human body that was generated with photographs of 0.33 mm cryosections, is converted into 3D acoustical maps. A full-wave nonlinear acoustic simulation tool is used to propagate ultrasound into the liver with a 2D intercostal ultrasound array in a $93 \times 39 \times 22$ mm domain with $6\times10^8$ points. Imaging metrics, based on the beamplots, root-mean-square phase aberration, spatial coherence lengths, and contrast-to-noise ratio are used to characterize the image quality. It is shown that the harmonic image quality is better than the fundamental image quality due, in part, to a narrower beam profile. The root-mean-square estimate of aberration after propagation through the simulated body wall is shown to be low (23.4 ns), consistently with previous reports of aberration measured experimentally in a human body wall. The spatial coherence measured at the transducer surface indicates that a transducer array element size of $<0.81 \lambda$ would be required to fully sample the acoustic field. These simulated three dimensional ultrasound images based directly on propagation physics provide a platform to investigate the sources of image degradation in three dimensions (included in Part II).


title: Control of chaos with minimal information transfer

abstract: This paper studies set-invariance and stabilization of hyperbolic sets over rate-limited channels for discrete-time control systems. We first investigate structural and control-theoretic properties of hyperbolic sets, in particular such that arise by adding small control terms to uncontrolled systems admitting (classical) hyperbolic sets. Then we derive a lower bound on the invariance entropy of a hyperbolic set in terms of the difference between the unstable volume growth rate and the measure-theoretic fiber entropy of associated random dynamical systems. We also prove that our lower bound is tight in two extreme cases. Furthermore, we apply our techniques to the problem of local uniform stabilization to a hyperbolic set. Finally, we discuss an example built on the H\'enon horseshoe.


title: Approximation algorithm for the Multicovering Problem

abstract: Let $\mathcal{H}=(V,\mathcal{E})$ be a hypergraph with maximum edge size $\ell$ and maximum degree $\Delta$. For given numbers $b_v\in \mathbb{N}_{\geq 2}$, $v\in V$, a set multicover in $\mathcal{H}$ is a set of edges $C \subseteq \mathcal{E}$ such that every vertex $v$ in $V$ belongs to at least $b_v$ edges in $C$. Set Multicover is the problem of finding a minimum-cardinality set multicover. Peleg, Schechtman and Wool conjectured that for any fixed $\Delta$ and $b:=\min_{v\in V}b_{v}$, the problem of \sbmultcov is not approximable within a ratio less than $\delta:=\Delta-b+1$, unless $\mathcal{P} =\mathcal{NP}$. Hence it's a challenge to explore for which classes of hypergraph the conjecture doesn't hold. We present a polynomial time algorithm for the Set Multicover problem which combines a deterministic threshold algorithm with conditioned randomized rounding steps. Our algorithm yields an approximation ratio of $ \max\left\{ \frac{148}{149}\delta, \left(1- \frac{ (b-1)e^{\frac{\delta}{4}}}{94\ell} \right)\delta \right\}$. Our result not only improves over the approximation ratio presented by Srivastav et al (Algorithmica 2016) but it's more general since we set no restriction on the parameter $\ell$. Moreover we present a further polynomial time algorithm with an approximation ratio of $\frac{5}{6}\delta$ for hypergraphs with $\ell\leq (1+\epsilon)\bar{\ell}$ for any fixed $\epsilon \in [0,\frac{1}{2}]$, where $\bar{\ell}$ is the average edge size. The analysis of this algorithm relies on matching/covering duality due to Ray-Chaudhuri (1960), which we convert into an approximative form. The second performance disprove the conjecture of peleg et al for a large subclass of hypergraphs.


title: A continuation approach to computing phase resetting curves

abstract: Phase resetting is a common experimental approach to investigating the behaviour of oscillating neurons. Assuming repeated spiking or bursting, a phase reset amounts to a brief perturbation that causes a shift in the phase of this periodic motion. The observed effects not only depend on the strength of the perturbation, but also on the phase at which it is applied. The relationship between the change in phase after the perturbation and the unperturbed old phase, the so-called phase resetting curve, provides information about the type of neuronal behaviour, although not all effects of the nature of the perturbation are well understood. In this chapter, we present a numerical method based on the continuation of a multi-segment boundary value problem that computes phase resetting curves in ODE models. Our method is able to deal effectively with phase sensitivity of a system, meaning that it is able to handle extreme variations in the phase resetting curve, including resets that are seemingly discontinuous. We illustrate the algorithm with two examples of planar systems, where we also demonstrate how qualitative changes of a phase resetting curve can be characterised and understood. A seven-dimensional example emphasises that our method is not restricted to planar systems, and illustrates how we can also deal with non-instantaneous, time-varying perturbations.


title: Increasing the Replicability for Linear Models via Adaptive Significance Levels

abstract: We put forward an adaptive alpha (Type I Error) that decreases as the information grows, for hypothesis tests in which nested linear models are compared. A less elaborate adaptation was already presented in \citet{PP2014} for comparing general i.i.d. models. In this article we present refined versions to compare nested linear models. This calibration may be interpreted as a Bayes-non-Bayes compromise, of a simple translations of a Bayes Factor on frequentist terms that leads to statistical consistency, and most importantly, it is a step towards statistics that promotes replicable scientific findings.


title: A Compact Fermion to Qubit Mapping

abstract: Mappings between fermions and qubits are valuable constructions in physics. To date only a handful exist. In addition to revealing dualities between fermionic and spin systems, such mappings are indispensable in any quantum simulation of fermionic physics on quantum computers. The number of qubits required per fermionic mode, and the locality of mapped fermionic operators strongly impact the cost of such simulations. We present a novel fermion to qubit mapping which outperforms all previous local mappings in both the qubit to mode ratio, and the locality of mapped operators.


title: Transient quantum beats, Rabi-oscillations and delay-time of modulated matter-waves

abstract: Transient phenomena of phase modulated cut-off wavepackets are explored by deriving an exact general solution to Schr\"odinger's equation for finite range potentials involving arbitrary initial quantum states. We show that the dynamical features of the probability density are governed by a virtual \textit{self-induced two-level system} with energies $E_{+}$, and $E_{-}$, due to the phase modulation of the initial state. The asymptotic probability density exhibits Rabi-oscillations characterized by a frequency $\Omega=(E_{+}-E_{-})/\hbar$, which are independent of the potential profile. It is also found that for a system with a bound state, the interplay between the virtual levels with the latter causes a \textit{quantum beat} effect with a beating frequency, $\Omega$. We also find a regime characterized by a \textit{time-diffraction} phenomenon that allows to measure unambiguously the delay-time, which can be described by an exact analytical formula. It is found that the delay-time agrees with the phase-time only for the case of strictly monochromatic waves.


title: Systematic studies and optimization of super sensitivity gaseous detectors of sparks, open flames and smoke

abstract: A review of the current progress in developments and tests of supersensitive gaseous detectors of open flame, sparks and smoke is given. A focus on the latest developments is on flat panel type sensors. This design, after further modification, will allow building not only a high efficiency detector but also offers UV imaging capability for flame visualization and several other applications. These studies were done in the framework of ATTRACT-SMART project.


title: Non-conventional Dynamical Bose Condensation

abstract: The paper presents a review of results concerning the non-conventional dynamical condensation versus conventional Bose-Einstein condensation, including the case of generalised van den Berg-Lewis-Pul\'e condensate. The review is based on detailed discussion of two models: a simple toy model and the Bogoliubov Weakly Imperfect Bose-Gas model, which was invented for explanation of superflfluidity of liquid 4He, but which is also instructive for analysis of non-conventional condensation regarding some recent interpretations of experimental data.


title: Observational signatures of strongly naked singularities: image of the thin accretion disk

abstract: We study the optical appearance of a thin accretion disk around the strongly naked static Janis-Newman-Winicour singularity. The solution does not possess a photon sphere, which results in the formation of a complex structure of bright rings in the central region of the disk image. Such structure is absent in the case of the Schwarzschild black hole with a thin accretion disk, where instead of the image we observe the black hole shadow. Some of the rings emit with the maximal observable radiation flux from the accretion disk, and should be experimentally detectable. Thus, this qualitatively new feature can be used to distinguish observationally black holes from naked singularities. We elucidate the appearance of the ring structure by revealing the physical mechanism of its formation, and explaining the nature of each of the ring images. We make the conjecture that a similar structure would also appear for other solutions without a photon sphere and it can serve as a general observational signature for distinguishing compact objects possessing no photon sphere from black holes.


title: Hyperspectral-Multispectral Image Fusion with Weighted LASSO

abstract: Spectral imaging enables spatially-resolved identification of materials in remote sensing, biomedicine, and astronomy. However, acquisition times require balancing spectral and spatial resolution with signal-to-noise. Hyperspectral imaging provides superior material specificity, while multispectral images are faster to collect at greater fidelity. We propose an approach for fusing hyperspectral and multispectral images to provide high-quality hyperspectral output. The proposed optimization leverages the least absolute shrinkage and selection operator (LASSO) to perform variable selection and regularization. Computational time is reduced by applying the alternating direction method of multipliers (ADMM), as well as initializing the fusion image by estimating it using maximum a posteriori (MAP) based on Hardie's method. We demonstrate that the proposed sparse fusion and reconstruction provides quantitatively superior results when compared to existing methods on publicly available images. Finally, we show how the proposed method can be practically applied in biomedical infrared spectroscopic microscopy.


title: Scene Completeness-Aware Lidar Depth Completion for Driving Scenario

abstract: This paper introduces Scene Completeness-Aware Depth Completion (SCADC) to complete raw lidar scans into dense depth maps with fine and complete scene structures. Recent sparse depth completion for lidars only focuses on the lower scenes and produces irregular estimations on the upper because existing datasets, such as KITTI, do not provide groundtruth for upper areas. These areas are considered less important since they are usually sky or trees of less scene understanding interest. However, we argue that in several driving scenarios such as large trucks or cars with loads, objects could extend to the upper parts of scenes. Thus depth maps with structured upper scene estimation are important for RGBD algorithms. SCADC adopts stereo images that produce disparities with better scene completeness but are generally less precise than lidars, to help sparse lidar depth completion. To our knowledge, we are the first to focus on scene completeness of sparse depth completion. We validate our SCADC on both depth estimate precision and scene-completeness on KITTI. Moreover, we experiment on less-explored outdoor RGBD semantic segmentation with scene completeness-aware D-input to validate our method.


title: Data-driven contact structures: from homogeneous mixing to multilayer networks

abstract: The modeling of the spreading of communicable diseases has experienced significant advances in the last two decades or so. This has been possible due to the proliferation of data and the development of new methods to gather, mine and analyze it. A key role has also been played by the latest advances in new disciplines like network science. Nonetheless, current models still lack a faithful representation of all possible heterogeneities and features that can be extracted from data. Here, we bridge a current gap in the mathematical modeling of infectious diseases and develop a framework that allows to account simultaneously for both the connectivity of individuals and the age-structure of the population. We compare different scenarios, namely, i) the homogeneous mixing setting, ii) one in which only the social mixing is taken into account, iii) a setting that considers the connectivity of individuals alone, and finally, iv) a multilayer representation in which both the social mixing and the number of contacts are included in the model. We analytically show that the thresholds obtained for these four scenarios are different. In addition, we conduct extensive numerical simulations and conclude that heterogeneities in the contact network are important for a proper determination of the epidemic threshold, whereas the age-structure plays a bigger role beyond the onset of the outbreak. Altogether, when it comes to evaluate interventions such as vaccination, both sources of individual heterogeneity are important and should be concurrently considered. Our results also provide an indication of the errors incurred in situations in which one cannot access all needed information in terms of connectivity and age of the population.