title: Weighted Monte Carlo with least squares and randomized extended Kaczmarz for option pricing

abstract: We propose a methodology for computing single and multi-asset European option prices, and more generally expectations of scalar functions of (multivariate) random variables. This new approach combines the ability of Monte Carlo simulation to handle high-dimensional problems with the efficiency of function approximation. Specifically, we first generalize the recently developed method for multivariate integration in [arXiv:1806.05492] to integration with respect to probability measures. The method is based on the principle "approximate and integrate" in three steps i) sample the integrand at points in the integration domain, ii) approximate the integrand by solving a least-squares problem, iii) integrate the approximate function. In high-dimensional applications we face memory limitations due to large storage requirements in step ii). Combining weighted sampling and the randomized extended Kaczmarz algorithm we obtain a new efficient approach to solve large-scale least-squares problems. Our convergence and cost analysis along with numerical experiments show the effectiveness of the method in both low and high dimensions, and under the assumption of a limited number of available simulations.


title: Direct Observation of Plasmon Band Formation and Delocalization in Quasi-Infinite Nanoparticle Chains

abstract: Chains of metallic nanoparticles sustain strongly confined surface plasmons with relatively low dielectric losses. To exploit these properties in applications,such as waveguides, the fabrication of long chains of low disorder and a thorough understanding of the plasmon-mode properties, such as dispersion relations, are indispensable. Here, we use a wrinkled template for directed self-assembly to assemble chains of gold nanoparticles. With this up-scalable method, chain lengths from two particles (140 nm) to 20 particles (1500 nm) and beyond can be fabricated. Electron energy-loss spectroscopy supported by boundary element simulations, finite-difference time-domain, and a simplified dipole coupling model reveal the evolution of a band of plasmonic waveguide modes from degenerated single-particle modes in detail. In striking difference from plasmonic rod-like structures, the plasmon band is confined in excitation energy, which allows light manipulations below the diffraction limit. The non-degenerated surface plasmon modes show suppressed radiative losses for efficient energy propagation over a distance of 1500 nm.


title: Modular Riesz bases versus Riesz bases in Hilbert $C^*$-Modules

abstract: In this paper, we give new characterizations of modular Riesz bases in Hilbert $C^*$-modules. We prove that modular Riesz bases share many properties with Riesz bases in Hilbert spaces. Moreover we show that there are also important differences; for example, there exist exact frames that are not modular Riesz bases.


title: A new INARMA(1, 1) model with Poisson marginals

abstract: We suggest an INARMA(1, 1) model with Poisson marginals which extends the INAR(1) in a similar way as the INGARCH(1, 1) does for the INARCH(1) model. The new model is equivalent to a binomially thinned INAR(1) process. This allows us to obtain some of its stochastic properties and use inference methods for hidden Markov models. The model is compared to various other models in two case studies.


title: A characterization of the weighted weak type Coifman-Fefferman and Fefferman-Stein inequalities

abstract: We introduce a variant of the $C_p$ condition (denoted by $SC_p$), and show that it characterizes weighted weak type versions of the classical Coifman-Fefferman and Fefferman-Stein inequalities.


title: Covert Wireless Communication in Presence of a Multi-Antenna Adversary and Delay Constraints

abstract: Covert communication hides the transmission of a message from a watchful adversary while ensuring reliable information decoding at the receiver, providing enhanced security in wireless communications. In this letter, covert communication in the presence of a multi-antenna adversary and under delay constraints is considered. Under the assumption of quasi-static wireless fading channels, we analyze the effect of increasing the number of antennas employed at the adversary on the achievable throughput of covert communication. It is shown that in contrast to a single-antenna adversary, a slight increase in the number of adversary's antennas drastically reduces the covert throughput, even for relaxed covertness requirements.


title: A Notion of Harmonic Clustering in Simplicial Complexes

abstract: We outline a novel clustering scheme for simplicial complexes that produces clusters of simplices in a way that is sensitive to the homology of the complex. The method is inspired by, and can be seen as a higher-dimensional version of, graph spectral clustering. The algorithm involves only sparse eigenproblems, and is therefore computationally efficient. We believe that it has broad application as a way to extract features from simplicial complexes that often arise in topological data analysis.


title: Type-I superconductivity in PdTe$_2$ probed by $\mu$SR

abstract: The Dirac semimetal PdTe$_2$ was recently reported to be a type-I superconductor with $T_c = 1.64$ K and a critical field $\mu_0H_c = 13.6$ mT. Since type-I superconductivity is unexpected for binary compounds, we have conducted muon spin rotation experiments to probe the superconducting phase on the microscopic scale via its intermediate state. For crystals with a finite demagnetization factor, $N$, the intermediate state forms in applied fields $(1-N)H_c < H_a < H_c$. We have carried out transverse field muon spin rotation measurements on a thin disk-like crystal with the field perpendicular to ($N_{\perp}=0.86$) and in the plane ($N_{\parallel}=0.08$) of the disk. By analysing the $\mu$SR signal we find that the volume fraction of the normal domains grows quasi-linearly with applied field at the expense of the Meissner domain fraction. This then provides solid evidence for the intermediate state and type-I superconductivity in the bulk of PdTe$_2$.


title: Stumbling Through the Research Wilderness, Standard Methods to Shine Light on Electrically Conductive Nanocomposites for Future Health-Care Monitoring

abstract: Electrically conductive nanocomposites are an exciting ever expanding area of research that has yielded many new technologies for wearable health devices. Acting as strain sensing materials, they have paved the way towards real time medical diagnostic tools that may very well lead to a golden age of healthcare. Currently, the goal in research is to create a material that simultaneously has both a large gauge factor G and sensing range. However, a weakness in the area of electromechanical research is the lack of standardisation in the reporting of the figure of merit, i.e. G, and the need for new metrics to give researchers a more complete view of the research landscape of resistive type sensors. A paradigm shift in the way in which data is reported is required, to push research in the right direction and to facilitate achieving research goals. Here, we report a standardised method for reporting strain sensing performance and the introduction of the working factor W and the Young's modulus Y of a material as two new material criteria. Using this new method, we can now for the first time define the benchmarks for an optimum sensing material, G > 7, W > 1, Y < 300 kPa, using limits set by standard commercial materials and the human body. Using extrapolated data from 200 publications normalised to this standard method, we can review what composite types meet these benchmark limits, what governs composite performances, the literary trends in composites and individual nanomaterial performance and the future prospects of research.


title: Sharp bound on the radial derivatives of the Zernike circle polynomials (disk polynomials)

abstract: We sharpen the bound $n^{2k}$ on the maximum modulus of the $k^{{\rm th}}$ radial derivative of the Zernike circle polynomials (disk polynomials) of degree $n$ to $n^2(n^2-1^2)\cdot ... \cdot(n^2-(k-1)^2)/2^k(1/2)_k$. This bound is obtained from a result of Koornwinder on the non-negativity of connection coefficients of the radial parts of the circle polynomials when expanded into a series of Chebyshev polynomials of the first kind. The new bound is shown to be sharp for, for instance, Zernike circle polynomials of degree $n$ and azimuthal order $m$ when $m=O(\sqrt{n})$ by using an explicit expression for the connection coefficients in terms of squares of Jacobi polynomials evaluated at 0. Keywords: Zernike circle polynomial, disk polynomial, radial derivative, Chebyshev polynomial, connection coefficient, Gegenbauer polynomial.


title: Unfolding of sparse high-energy $\gamma$-ray spectra from LaBr$_{3}$:Ce detectors

abstract: Here we report on the characterization of one of the large-volume LaBr$_{3}$:Ce detectors for the ELIGANT project at ELI-NP. The main focus of this work is the response function for high-energy $\gamma$ rays of such detectors. In particular, we compare a selection of unfolding methods to resolve small structures in $\gamma$-ray spectra with high-energies. Three methods have been compared using $\gamma$-ray spectra with energies up to 12 MeV obtained in an experiment at the 3 MV Tandetron\texttrademark\ facility at IFIN-HH. The results show that the iterative unfolding approach gives the best qualitative reproduction of the emitted $\gamma$-ray spectrum. Furthermore, the correlation fluctuations in high-energy regime from the iterative method are two orders of magnitude smaller than when using the matrix inversion approach with second derivative regularization. In addition, the iterative method is computationally faster as it does not contain large matrix inversions. The matrix inversion method does, however, give more consistent results over the full energy range and in the low-statistics limit. Our conclusion is that the performance of the iterative approach makes it well suitable for semi-online analysis of experimental data. These results will be important, both for experiments with the ELIGANT setup, and for on-line diagnostics of the energy spread of the $\gamma$-ray beam which is under implementation at ELI-NP.


title: Distortion matrix concept for deep optical imaging in scattering media

abstract: In optical imaging, light propagation is affected by the inhomogeneities of the medium. Sample-induced aberrations and multiple scattering can strongly degrade the image resolution and contrast. Based on a dynamic correction of the incident and/or reflected wave-fronts, adaptive optics has been employed to compensate for those aberrations. However, it only applies to spatially-invariant aberrations or to thin aberrating layers. Here, we propose a global and non-invasive approach based on the distortion matrix concept. This matrix basically connects any focusing point of the image with the distorted part of its wave-front in reflection. A singular value decomposition of the distortion matrix allows to correct for high-order aberrations and forward multiple scattering over multiple isoplanatic modes. Proof-of-concept experiments are performed through biological tissues including a turbid cornea. We demonstrate a Strehl ratio enhancement up to 2500 and recover a diffraction-limited resolution until a depth of ten scattering mean free paths.


title: Guan-Li type mean curvature flow for free boundary hypersurfaces in a ball

abstract: In this paper we introduce a Guan-Li type volume preserving mean curvature flow for free boundary hypersurfaces in a ball. We give a concept of star-shaped free boundary hypersurfaces in a ball and show that the Guan-Li type mean curvature flow has long time existence and converges to a free boundary spherical cap, provided the initial data is star-shaped.


title: Audio-Conditioned U-Net for Position Estimation in Full Sheet Images

abstract: The goal of score following is to track a musical performance, usually in the form of audio, in a corresponding score representation. Established methods mainly rely on computer-readable scores in the form of MIDI or MusicXML and achieve robust and reliable tracking results. Recently, multimodal deep learning methods have been used to follow along musical performances in raw sheet images. Among the current limits of these systems is that they require a non trivial amount of preprocessing steps that unravel the raw sheet image into a single long system of staves. The current work is an attempt at removing this particular limitation. We propose an architecture capable of estimating matching score positions directly within entire unprocessed sheet images. We argue that this is a necessary first step towards a fully integrated score following system that does not rely on any preprocessing steps such as optical music recognition.


title: Derived deformation theory of algebraic structures

abstract: The main purpose of this article is to develop an explicit derived deformation theory of algebraic structures at a high level of generality, encompassing in a common framework various kinds of algebras (associative, commutative, Poisson...) or bialgebras (associative and coassociative, Lie, Frobenius...), that is algebraic structures parametrized by props. A central aspect is that we define and study moduli spaces of deformations of algebraic structures up to quasi-isomorphisms (and not just isotopies or isomorphisms). To do so, we implement methods coming from derived algebraic geometry, by encapsulating these deformation theories as classifying (pre)stacks with good infinitesimal properties and derived formal geometry, by means of derived formal moduli problems and derived formal groups. In particular, we prove that the Lie algebra describing the deformation theory of an object in a given $\infty$-category of dg algebras can be obtained equivalently as the tangent complex of loops on a derived quotient of this moduli space by the homotopy automorphims of this object. Moreover, we provide explicit formulae for such derived deformation problems of algebraic structures up to quasi-isomorphisms and relate them in a precise way to other standard deformation problems of algebraic structures. This relation is given by a fiber sequence of the associated dg-Lie algebras of their deformation complexes. Our results provide simultaneously a vast generalization of standard deformation theory of algebraic structures which is suitable (and needed) to set up algebraic deformation theory both at the $\infty$-categorical level and at a higher level of generality than algebras over operads. In addition, we study a general criterion to compare formal moduli problems of different algebraic structures and apply our formalism to $E_n$-algebras and bialgebras.


title: Sparsity of Hawking Radiation in $D+1$ Space-Time Dimensions Including Particle Masses

abstract: Hawking radiation from an evaporating black hole has often been compared to black body radiation. However, this comparison misses an important feature of Hawking radiation: Its low density of states. This can be captured in an easy to calculate, heuristic, and semi-analytic measure called "sparsity". In this letter we shall present both the concept of sparsities and its application to $D+1$-dimensional Tangherlini black holes and their evaporation. In particular, we shall also publish for the first time sparsity expressions taking into account in closed form effects of non-zero particle mass. We will also see how this comparatively simple method reproduces results of (massless) Hawking radiation in higher dimensions and how different spins contribute to the total radiation in this context.


title: Analysis of silicon nitride partial Euler waveguide bends

abstract: In this work we present a detailed analysis of individual loss mechanisms in silicon nitride partial Euler bends at a wavelength of 850 nm. This structure optimizes the transmission through small radii optical waveguide bends. The partial Euler bend geometry balances losses arising from the transition from the straight to the bend waveguide mode, and radiative losses of the bend waveguide mode. Numerical analyses are presented for 45-degree bends commonly employed in S-bend configurations to create lateral offsets, as well as 90- and 180-degree bends. Additionally, 90-degree partial Euler bends were fabricated on a silicon nitride photonic platform to experimentally complement the theoretical findings. The optimized waveguide bends allow for a reduced effective radius without increasing the total bend loss and, thus, enable a higher component density in photonic integrated circuits.


title: The Carbon State in Dilute Germanium Carbides

abstract: Conduction and valence band states for the highly mismatched alloy (HMA) Ge:C are projected onto Ge crystal states, Ge vacancy states, and Ge/C atomic orbitals, revealing that substitutional carbon not only creates a direct bandgap, but the new conduction band is optically active. Overlap integrals of the new Ge:C conduction band with bands of pure Ge shows the new band has almost no Ge band character. C sites structurally mimic uncharged vacancies in the Ge lattice, similar to Hjalmarson's model for other HMAs. C perturbs the entire Ge band structure even at the deepest crystal core energy levels. Projection onto atomic sites shows relatively weak localization compared with other HMAs, but does show a strong anisotropy in probability distribution. L-valley conduction band states in Ge are ruled out as major contributors to the carbon state in Ge:C, both by weak inner products between these states and by a negligible effect on optical transition strength when adding C.


title: Unified description of sound velocities in strongly coupled Yukawa systems of different spatial dimensionality

abstract: Sound velocities in classical single-component fluids with Yukawa (screened Coulomb) interactions are systematically evaluated and analyzed in one-, two-, and three spatial dimensions (${\mathcal D}=1,2,3$). In the strongly coupled regime the convenient sound velocity scale is given by $\sqrt{Q^2/\Delta m}$, where $Q$ is the particle charge, $m$ is the particle mass, $n$ is the particle density, and $\Delta=n^{-1/{\mathcal D}}$ is the unified interparticle distance. The sound velocity can be expressed as a product of this scaling factor and a dimension-dependent function of the screening parameter, $\kappa=\Delta/\lambda$, where $\lambda$ is the screening length. A unified approach is used to derive explicit expressions for these dimension-dependent functions in the weakly screened regime ($\kappa\lesssim 3$). It is also demonstrated that for stronger screening ($\kappa\gtrsim 3$), the effect of spatial dimensionality virtually disappears, the longitudinal sound velocities approach a common asymptote, and a one-dimensional nearest-neighbor approximation provides a relatively good estimate for this asymptote. This result is not specific to the Yukawa potential, but equally applies to other classical systems with steep repulsive interactions. An emerging relation to a popular simple freezing indicator is briefly discussed. Overall, the results can be useful when Yukawa interactions are relevant, in particular, in the context of complex (dusty) plasmas and colloidal suspensions.


title: Global Existence for some Cross Diffusion Systems with Equal Cross Diffusion/Reaction Rates

abstract: We consider some cross diffusion systems which is inspired by models in mathematical biology/ecology, in particular the Shigesada-Kawasaki-Teramoto (SKT) model in population biology. We establish the global existence of strong solutions to systems for multiple species having equal either diffusion or reaction rates. The systems are given on bounded domains of arbitrary dimension.


title: Electronic structure of semiconductor nanostructures: A modified localization landscape theory

abstract: In this paper we present a modified localization landscape theory to calculate localized/confined electron and hole states and the corresponding energy eigenvalues without solving a (large) eigenvalue problem. We motivate and demonstrate the benefit of solving $\hat{H}^2u=1$ in the modified localization landscape theory in comparison to $\hat{H}u=1$, solved in the localization landscape theory. We detail the advantages by fully analytic considerations before targeting the numerical calculation of electron and hole states and energies in III-N heterostructures. We further discuss how the solution of $\hat{H}^2u=1$ is used to extract an effective potential $W$ that is comparable to the effective potential obtained from $\hat{H}u=1$, ensuring that it can for instance be used to introduce quantum corrections to drift-diffusion transport calculations. Overall, we show that the proposed modified localization landscape theory keeps all the benefits of the recently introduced localization landscape theory but further improves factors such as convergence of the calculated energies and the robustness of the method against the chosen integration region for $u$ to obtain the corresponding energies. We find that this becomes especially important for here studied $c$-plane InGaN/GaN quantum wells with higher In contents. All these features make the proposed approach very attractive for calculation of localized states in highly disordered systems, where partitioning the systems into different subregions can be difficult.


title: Quasimaps and some examples of stacks for everybody

abstract: This note introduces the theory of quasimaps to GIT quotients with intuition and concrete examples, with the goal of explaining a closed formula for the quasimap $I$-function. Along the way, it emphasizes aspects of this story that illustrate general stacky concepts.


title: Bimodular continuous attractor neural networks with static and moving stimuli

abstract: We investigated the dynamical behaviors of bimodular continuous attractor neural networks, each processing a modality of sensory input and interacting with each other. We found that when bumps coexist in both modules, the position of each bump is shifted towards the other input when the intermodular couplings are excitatory and is shifted away when inhibitory. When one intermodular coupling is excitatory while another is moderately inhibitory, temporally modulated population spikes can be generated. On further increase of the inhibitory coupling, momentary spikes will emerge. In the regime of bump coexistence, bump heights are primarily strengthened by excitatory intermodular couplings, but there is a lesser weakening effect due to a bump being displaced from the direct input. When bimodular networks serve as decoders of multisensory integration, we extend the Bayesian framework to show that excitatory and inhibitory couplings encode attractive and repulsive priors, respectively. At low disparity, the bump positions decode the posterior means in the Bayesian framework, whereas at high disparity, multiple steady states exist. In the regime of multiple steady states, the less stable state can be accessed if the input causing the more stable state arrives after a sufficiently long delay. When one input is moving, the bump in the corresponding module is pinned when the moving stimulus is weak, unpinned at intermediate stimulus strength, and tracks the input at strong stimulus strength, and the stimulus strengths for these transitions increase with the velocity of the moving stimulus. These results are important to understanding multisensory integration of static and dynamic stimuli.


title: Perturbed Euler top and bifurcation of limit cycles on invariant Casimir surfaces

abstract: Analytical perturbations of the Euler top are considered. The perturbations are based on the Poisson structure for such a dynamical system, in such a way that the Casimir invariants of the system remain invariant for the perturbed flow. By means of the Poincar\'{e}-Pontryagin theory, the existence of limit cycles on the invariant Casimir surfaces for the perturbed system is investigated up to first order of perturbation, providing sharp bounds for their number. Examples are given.

categories: math-ph math.CA math.DS math.MP nlin.SI physics.class-ph


title: Optimising Individual-Treatment-Effect Using Bandits

abstract: Applying causal inference models in areas such as economics, healthcare and marketing receives great interest from the machine learning community. In particular, estimating the individual-treatment-effect (ITE) in settings such as precision medicine and targeted advertising has peaked in application. Optimising this ITE under the strong-ignorability-assumption -- meaning all confounders expressing influence on the outcome of a treatment are registered in the data -- is often referred to as uplift modeling (UM). While these techniques have proven useful in many settings, they suffer vividly in a dynamic environment due to concept drift. Take for example the negative influence on a marketing campaign when a competitor product is released. To counter this, we propose the uplifted contextual multi-armed bandit (U-CMAB), a novel approach to optimise the ITE by drawing upon bandit literature. Experiments on real and simulated data indicate that our proposed approach compares favourably against the state-of-the-art. All our code can be found online at https://github.com/vub-dl/u-cmab.


title: A new method for flow-based network intrusion detection using the inverse Potts model

abstract: Network Intrusion Detection Systems (NIDS) play an important role as tools for identifying potential network threats. In the context of ever-increasing traffic volume on computer networks, flow-based NIDS arise as good solutions for real-time traffic classification. In recent years, different flow-based classifiers have been proposed using Machine Learning (ML) algorithms. Nevertheless, classical ML-based classifiers have some limitations. For instance, they require large amounts of labeled data for training, which might be difficult to obtain. Additionally, most ML-based classifiers are not capable of domain adaptation, i.e. after being trained on an specific data distribution, they are not general enough to be applied to other related data distributions. And, finally, many of the models inferred by these algorithms are black boxes, which do not provide explainable results. To overcome these limitations, we propose a new algorithm, called Energy-based Flow Classifier (EFC). This anomaly-based classifier uses inverse statistics to infer a statistical model based on labeled benign examples. We show that EFC is capable of accurately performing binary flow classification and is more adaptable to different data distributions than classical ML-based classifiers. Given the positive results obtained on three different datasets (CIDDS-001, CICIDS17 and CICDDoS19), we consider EFC to be a promising algorithm to perform robust flow-based traffic classification.


title: Long Optimal or Small-Defect LRC Codes with Unbounded Minimum Distances

abstract: A code over a finite field is called locally recoverable code (LRC) if every coordinate symbol can be determined by a small number (at most r, this parameter is called locality) of other coordinate symbols. For a linear code with length n, dimension k and locality r, its minimum distance d satisfies a Singleton-like bound. A code attaining this bound is called optimal. Many families of optimal locally recoverable codes have been constructed by using different techniques in finite fields or algebraic curves. However no optimal LRC code over a general finite field of q elements with the length n around the square of q, the locality r larger than or equal to 24 and the minimum distance d larger than or equal to 9 has been constructed. In this paper for any given finite field of q elements, any given r between 1 and q-1 and given d in certain range, we give an optimal LRC code with length n around the square of q, locality r and minimum distance d. This is the only known family of optimal LRC codes with lengths n around the square of q and unbounded localities and minimum distances d larger than or equal to 9. We also gives an asymptotic bound for q-ary r locality LRC codes better than the previously konwn bound. Many long r-locality LRC codes with small defects are also constructed.


title: Exploring the fitness landscape of a realistic turbofan rotor blade optimization

abstract: Aerodynamic shape optimization has established itself as a valuable tool in the engineering design process to achieve highly efficient results. A central aspect for such approaches is the mapping from the design parameters which encode the geometry of the shape to be improved to the quality criteria which describe its performance. The choices to be made in the setup of the optimization process strongly influence this mapping and thus are expected to have a profound influence on the achievable result. In this work we explore the influence of such choices on the effects on the shape optimization of a turbofan rotor blade as it can be realized within an aircraft engine design process. The blade quality is assessed by realistic three dimensional computational fluid dynamics (CFD) simulations. We investigate the outcomes of several optimization runs which differ in various configuration options, such as optimization algorithm, initialization, number of degrees of freedom for the parametrization. For all such variations, we generally find that the achievable improvement of the blade quality is comparable for most settings and thus rather insensitive to the details of the setup. On the other hand, even supposedly minor changes in the settings, such as using a different random seed for the initialization of the optimizer algorithm, lead to very different shapes. Optimized shapes which show comparable performance usually differ quite strongly in their geometries over the complete blade. Our analyses indicate that the fitness landscape for such a realistic turbofan rotor blade optimization is highly multi-modal with many local optima, where very different shapes show similar performance.


title: "Get a Free Item Pack with Every Activation!" -- Do Incentives Increase the Adoption Rates of Two-Factor Authentication?

abstract: Account security is an ongoing issue in practice. Two-Factor Authentication (2FA) is a mechanism which could help mitigate this problem, however adoption is not very high in most domains. Online gaming has adopted an interesting approach to drive adoption: Games offer small rewards such as visual modifications to the player's avatar's appearance, if players utilize 2FA. In this paper, we evaluate the effectiveness of these incentives and investigate how they can be applied to non-gaming contexts. We conducted two surveys, one recruiting gamers and one recruiting from a general population. In addition, we conducted three focus group interviews to evaluate various incentive designs for both, the gaming context and the non-gaming context. We found that visual modifications, which are the most popular type of gaming-related incentives, are not as popular in non-gaming contexts. However, our design explorations indicate that well-chosen incentives have the potential to lead to more users adopting 2FA, even outside of the gaming context.


title: Thermal emission from the amorphous dust: An alternative possibility of the origin of the anomalous microwave emission

abstract: Complete studies of the radiative processes of thermal emission from the amorphous dust from microwave through far infrared wavebands are presented by taking into account, self-consistently for the first time, the standard two-level systems (TLS) model of amorphous materials. The observed spectral energy distributions (SEDs) for the Perseus molecular cloud (MC) and W43 from microwave through far infrared are fitted with the SEDs calculated with the TLS model of amorphous silicate. We have found that the model SEDs well reproduce the observed properties of the anomalous microwave emission (AME). The present result suggests an alternative interpretation for the AME being carried by the resonance emission of the TLS of amorphous materials without introducing new species. Simultaneous fitting of the intensity and polarization SEDs for the Perseus MC and W43 are also performed. The amorphous model reproduces the overall observed feature of the intensity and polarization SEDs of the Perseus MC and W43. However, the model's predicted polarization fraction of the AME is slightly higher than the QUIJOTE upper limits in several frequency bands. A possible improvement of our model to resolve this problem is proposed. Our model predicts that interstellar dust is amorphous materials having very different physical characteristics compared with terrestrial amorphous materials.


title: Representation of Polytopes as Polynomial Zonotopes

abstract: We prove that each bounded polytope can be represented as a polynomial zonotope, which we refer to as the Z-representation of polytopes. Previous representations are the vertex representation (V-representation) and the halfspace representation (H-representation). Depending on the polytope, the Z-representation can be more compact than the V-representation and the H-representation. In addition, the Z-representation enables the computation of linear maps, Minkowski addition, and convex hull with a computational complexity that is polynomial in the representation size. The usefulness of the new representation is demonstrated by range bounding within polytopes.


title: Nonlocal gauge equivalence: Hirota versus extended continuous Heisenberg and Landau-Lifschitz equation

abstract: We exploit the gauge equivalence between the Hirota equation and the extended continuous Heisenberg equation to investigate how nonlocality properties of one system are inherited by the other. We provide closed generic expressions for nonlocal multi-soliton solutions for both systems. By demonstrating that a specific auto-gauge transformation for the extended continuous Heisenberg equation becomes equivalent to a Darboux transformation, we use the latter to construct the nonlocal multi-soliton solutions from which the corresponding nonlocal solutions to the Hirota equation can be computed directly. We discuss properties and solutions of a nonlocal version of the nonlocal extended Landau-Lifschitz equation obtained from the nonlocal extended continuous Heisenberg equation or directly from the nonlocal solutions of the Hirota equation.


title: Sailing over three problems of Koszmider

abstract: We discuss three problems of Koszmider on the structure of the spaces of continuous functions on the Stone compact $K_{\mathcal A}$ generated by an almost disjoint family $\mathcal A$ of infinite subsets of $\omega$ -- we present a solution to two problems and develop a previous results of Marciszewski and Pol answering the third one. We will show, in particular, that assuming Martin's axiom the space $C(K_{\mathcal A})$ is uniquely determined up to isomorphism by the cardinality of $\mathcal A$ whenever $|{\mathcal A}|<{\mathfrak c}$, while there are $2^{\mathfrak c}$ nonisomorphic spaces $C(K_{\mathcal A})$ with $|{\mathcal A}|= {\mathfrak c}$. We also investigate Koszmider's problems in the context of the class of separable Rosenthal compacta and indicate the meaning of our results in the language of twisted sums of $c_0$ and some $C(K)$ spaces.


title: Universal Bounds for Size and Energy of Codes of Given Minimum and Maximum Distances

abstract: We employ signed measures that are positive definite up to certain degrees to establish Levenshtein-type upper bounds on the cardinality of codes with given minimum and maximum distances, and universal lower bounds on the potential energy (for absolutely monotone interactions) for codes with given maximum distance and cardinality. The distance distributions of codes that attain the bounds are found in terms of the parameters of Levenshtein-type quadrature formulas. Necessary and sufficient conditions for the optimality of our bounds are derived. Further, we obtain upper bounds on the energy of codes of fixed minimum and maximum distances and cardinality.


title: Algebraic Semantics for Relative Truth, Awareness, and Possibility

abstract: This paper puts forth a class of algebraic structures, relativized Boolean algebras (RBAs), that provide semantics for propositional logic in which truth/validity is only defined relative to a local domain. In particular, the join of an event and its complement need not be the top element. Nonetheless, behavior is locally governed by the laws of propositional logic. By further endowing these structures with operators -- akin to the theory of modal Algebras -- RBAs serve as models of modal logics in which truth is relative. In particular, modal RBAs provide semantics for various well known awareness logics and an alternative view of possibility semantics.


title: WASP-South detection of HD219666b transits provides an accurate ephemeris

abstract: The hot-Neptune HD219666b (= TOI-118.01) was an early discovery from Sector 1 of the TESS transit survey. Being a rare "Neptune desert" planet, and transiting a bright, V = 9.9 star, it is a prime target for further study including atmospheric characterisation. Targeting exoplanet transits with major facilities such as HST and the imminent JWST depends on an accurate ephemeris. Here we report a recovery of the transit in WASP-South data, enabling an ephemeris from data spanning 8 years that has a period uncertainty which is smaller by a factor 60 compared to the discovery paper.


title: CYGNO: Triple-GEM Optical Readout for Directional Dark Matter Search

abstract: CYGNO is a project realising a cubic meter demonstrator to study the scalability of the performance of the optical approach for the readout of large-volume, GEM-equipped TPC. This is part of the CYGNUS proto-collaboration which aims at constructing a network of underground observatories for directional Dark Matter search. The combined use of high-granularity sCMOS and fast sensors for reading out the light produced in GEM channels during the multiplication processes was shown to allow on one hand to reconstruct 3D direction of the tracks, offering accurate energy measurements and sensitivity to the source directionality and, on the other hand, a high particle identification capability very useful to distinguish nuclear recoils. Results of the performed R&D and future steps toward a 30-100 cubic meter experiment will be presented.


title: On the Relation between Weak Completion Semantics and Answer Set Semantics

abstract: The Weak Completion Semantics (WCS) is a computational cognitive theory that has shown to be successful in modeling episodes of human reasoning. As the WCS is a recently developed logic programming approach, this paper investigates the correspondence of the WCS with respect to the well-established Answer Set Semantics (ASP). The underlying three-valued logic of both semantics is different and their models are evaluated with respect to different program transformations. We first illustrate these differences by the formal representation of some examples of a well-known psychological experiment, the suppression task. After that, we will provide a translation from logic programs understood under the WCS into logic programs understood under the ASP. In particular, we will show that logic programs under the WCS can be represented as logic programs under the ASP by means of a definition completion, where all defined atoms in a program must be false when their definitions are false.


title: Zero temperature limits of equilibrium states for subadditive potentials and approximation of the maximal Lyapunov exponent

abstract: In this paper we study ergodic optimization problems for subadditive sequences of functions on a topological dynamical system. We prove that for $t\rightarrow \infty$ any accumulation point of a family of equilibrium states is a maximizing measure. We show that the Lyapunov exponent and entropy of equilibrium states converge in the limit $t\rightarrow \infty$ to the maximum Lyapunov exponent and entropy of maximizing measures. In the particular case of matrix cocycles we prove that the maximal Lyapunov exponent can be approximated by Lyapunov exponents of periodic trajectories under certain assumptions.


title: LIPSS on thin metallic films: New insights from multiplicity of laser-excited electromagnetic modes and efficiency of metal oxidation

abstract: Thin Cr films 28-nm thick deposited on glass substrates were processed by scanning low-intensity femtosecond laser pulses with energy well below single-pulse damage threshold. Two types of laser-induced periodic surface structures (LIPSS) were produced, depending on the scanning velocity, (1) parallel to laser light polarization with periodicity somewhat smaller than laser wavelength and (2) perpendicular to polarization with spatial period much smaller than wavelength. All structures are formed as protrusions above the initial film surface and exhibit a high degree of oxidation. To explain formation of the LIPSS and their conversion from one to another type, a rigorous numerical approach for modeling surface electromagnetic waves in thin-film geometry has been developed, which takes into account the change of optical properties of material due to laser-induced oxidation and porosity. The approach addresses the multiplicity of electromagnetic modes allowed for thin films. It has been found that the low spatial frequency LIPSS parallel to laser polarization, which are formed at low scanning velocities, are well described by the Sipe theory for surfaces of low roughness. The SEW mode responsible for high spatial frequency LIPSS formation at high scanning velocities has been identified. The mechanisms of optical feedback and transformation between types of LIPSS with scanning velocity have been proposed.


title: The diameter and radius of radially maximal graphs

abstract: A graph is called radially maximal if it is not complete and the addition of any new edge decreases its radius. In 1976 Harary and Thomassen proved that the radius $r$ and diameter $d$ of any radially maximal graph satisfy $r\le d\le 2r-2.$ Dutton, Medidi and Brigham rediscovered this result with a different proof in 1995 and they posed the conjecture that the converse is true, that is, if $r$ and $d$ are positive integers satisfying $r\le d\le 2r-2,$ then there exists a radially maximal graph with radius $r$ and diameter $d.$ We prove this conjecture and a little more.


title: Renormalization group approach to power-law modeling of complex metabolic networks

abstract: In the modeling of complex biological systems, the use of power-law models (such as S-systems and GMA systems) often provides a remarkable accuracy over several orders of magnitude in concentrations, an unusually broad range not fully understood at present. In order to provide additional insight in this sense, this article is devoted to the renormalization group analysis of reactions in fractal or self-similar media. In particular, the renormalization group methodology is applied to the investigation of how rate-laws describing such reactions are transformed when the geometric scale is changed. The precise purpose of such analysis is to investigate whether or not power-law rate-laws present some remarkable features accounting for the successes of power-law modeling. As we shall see, according to the renormalization group point of view the answer is positive, as far as power-laws are the critical solutions of the renormalization group transformation, namely power-law rate-laws are the renormalization group invariant solutions. Moreover, it is shown that these results also imply invariance under the group of concentration scalings, thus accounting for the reported power-law model accuracy over several orders of magnitude in metabolite concentrations.

categories: physics.bio-ph math.DS nlin.AO physics.chem-ph q-bio.MN


title: FISHDBC: Flexible, Incremental, Scalable, Hierarchical Density-Based Clustering for Arbitrary Data and Distance

abstract: FISHDBC is a flexible, incremental, scalable, and hierarchical density-based clustering algorithm. It is flexible because it empowers users to work on arbitrary data, skipping the feature extraction step that usually transforms raw data in numeric arrays letting users define an arbitrary distance function instead. It is incremental and scalable: it avoids the $\mathcal O(n^2)$ performance of other approaches in non-metric spaces and requires only lightweight computation to update the clustering when few items are added. It is hierarchical: it produces a "flat" clustering which can be expanded to a tree structure, so that users can group and/or divide clusters in sub- or super-clusters when data exploration requires so. It is density-based and approximates HDBSCAN*, an evolution of DBSCAN.


title: Overview of CMOS sensors for future tracking detectors

abstract: The depleted CMOS sensors are emerging as one of the main candidate technologies for future tracking detectors in high luminosity colliders. Its capability of integrating the sensing diode into the CMOS wafer hosting the front-end electronics allows for reduced noise and higher signal sensitivity. They are suitable for high radiation environments due to the possibility of applying high depletion voltage and the availability of relatively high resistivity substrates. The use of a CMOS commercial fabrication process leads to their cost reduction and allows faster construction of large area detectors. A general perspective of the state of the art of these devices will be given in this contribution as well as a summary of the main developments carried out with regard to these devices in the framework of the CERN RD50 collaboration.


title: The rotation angle distribution underlying magnetic field fluctuations in the 1/f range of solar wind turbulent spectra

abstract: We discuss properties of large amplitude magnetic field fluctuations during fast Alfvenic solar wind streams, focussing on the statistics of the rotation angle between consecutive magnetic field vector measurements for different scales in the plasma. Since in the fast solar wind fluctuations preserve the modulus of the magnetic field to a good approximation, the tip of the magnetic field vector is observed to move on a sphere of approximately constant radius |B|. We then compare statistics of solar wind measurements with that of a simple model of a random walk bounded on a spherical surface. The analogy consists in the fact that in both systems the geometrical constraint imposes a limiting amplitude at large separations and thus introduces a break scale in the power spectrum of the fluctuations, leading to a shallower slope for scales where the fluctuations amplitude becomes scale-independent. However, while in the case of the random walk the saturation of the fluctuations occurs when the pattern becomes uniform on the sphere (flat distribution of the cosine of the rotation angle), transitioning then to a white noise regime, in the solar wind magnetic field fluctuations saturate in amplitude maintaining a preferential direction. We suggest that this behaviour is due to the presence of the background interplanetary magnetic field, which keeps some long-range memory in the system also when the fluctuations becomes independent of the scale. This long-range correlation is a necessary ingredient in order to produce the 1/f spectrum observed at large scales in the solar wind.


title: A study on the dynamic spectral indices for SEP events on 2000 July 14 and 2005 January 20

abstract: We have studied the dynamic proton spectra for the two solar energetic particle (SEP) events on 2000 July 14 (hereafter GLE59) and 2005 January 20 (hereafter GLE69). The source locations of GLE59 and GLE69 are N22W07 and N12W58 respectively. Proton fluxes >30 MeV have been used to compute the dynamic spectral indices of the two SEP events. The results show that spectral indices of the two SEP events increased more swiftly at early times, suggesting that the proton fluxes >30 MeV might be accelerated particularly by the concurrent flares at early times for the two SEP events. For the GLE69 with source location at N12W58, both flare site and shock nose are well connected with the Earth at the earliest time. However, only the particles accelerated by the shock driven by eastern flank of the CME can propagate along the interplanetary magnetic field line to the Earth after the flare. For the GLE59 with source location at N22W07, only the particles accelerated by the shock driven by western flank of the associated CME can reach the Earth after the flare. Results show that there was slightly more than one hour during which the proton spectra for GLE69 are softer than that for GLE59 after the flares, suggesting that the shock driven by eastern flank of the CME associated with GLE69 is weaker than the shock driven by the western flank of the CME associated with GLE59. The results support that quasi-perpendicular shock has stronger potential in accelerating particles than the quasi-parallel shock. The results also suggest that only a small part of the shock driven by western flank of the CME associated with the GLE59 is quasi-perpendicular.


title: Continuous-Domain Assignment Flows

abstract: Assignment flows denote a class of dynamical models for contextual data labeling (classification) on graphs. We derive a novel parametrization of assignment flows that reveals how the underlying information geometry induces two processes for assignment regularization and for gradually enforcing unambiguous decisions, respectively, that seamlessly interact when solving for the flow. Our result enables to characterize the dominant part of the assignment flow as a Riemannian gradient flow with respect to the underlying information geometry. We consider a continuous-domain formulation of the corresponding potential and develop a novel algorithm in terms of solving a sequence of linear elliptic PDEs subject to a simple convex constraint. Our result provides a basis for addressing learning problems by controlling such PDEs in future work.


title: Existence and smoothness of the density of the solution to fractional stochastic integral Volterra equations

abstract: We consider stochastic Volterra integral equations driven by a fractional Brownian motion with Hurst parameter H > 1/2 . We first derive supremum norm estimates for the solution and its Malliavin derivative. We then show existence and smoothness of the density under suitable nondegeneracy conditions. This extends the results in Hu-Nualart and Nualart-Saussereau where stochastic differential equations driven by fractional Brownian motion are considered. The proof uses a priori estimates for deterministic differential equations driven by a function in a suitable Sobolev space.


title: Small-scale clustering of nano-dust grains in supersonic turbulence

abstract: We investigate the clustering and dynamics of nano-sized particles (nano-dust) in high-resolution ($1024^3$) simulations of compressible isothermal hydrodynamic turbulence. It is well-established that large grains will decouple from a turbulent gas flow, while small grains will tend to trace the motion of the gas. We demonstrate that nano-sized grains may cluster in a turbulent flow (fractal small-scale clustering), which increases the local grain density by at least a factor of a few. In combination with the fact that nano-dust grains may be abundant in general, and the increased interaction rate due to turbulent motions, aggregation involving nano dust may have a rather high probability. Small-scale clustering will also affect extinction properties. As an example we present an extinction model based on silicates, graphite and metallic iron, assuming strong clustering of grain sizes in the nanometre range, could explain the extreme and rapidly varying ultraviolet extinction in the host of GRB 140506A.


title: Understanding the origin of the magnetic field morphology in the wide-binary protostellar system BHR 71

abstract: We present 1.3 mm ALMA observations of polarized dust emission toward the wide-binary protostellar system BHR 71 IRS1 and IRS2. IRS1 features what appears to be a natal, hourglass-shaped magnetic field. In contrast, IRS2 exhibits a magnetic field that has been affected by its bipolar outflow. Toward IRS2, the polarization is confined mainly to the outflow cavity walls. Along the northern edge of the redshifted outflow cavity of IRS2, the polarized emission is sandwiched between the outflow and a filament of cold, dense gas traced by N$_2$D$^+$, toward which no dust polarization is detected. This suggests that the origin of the enhanced polarization in IRS2 is the irradiation of the outflow cavity walls, which enables the alignment of dust grains with respect to the magnetic field -- but only to a depth of ~300 au, beyond which the dust is cold and unpolarized. However, in order to align grains deep enough in the cavity walls, and to produce the high polarization fraction seen in IRS2, the aligning photons are likely to be in the mid- to far-infrared range, which suggests a degree of grain growth beyond what is typically expected in very young, Class 0 sources. Finally, toward IRS1 we see a narrow, linear feature with a high (10-20%) polarization fraction and a well ordered magnetic field that is not associated with the bipolar outflow cavity. We speculate that this feature may be a magnetized accretion streamer; however, this has yet to be confirmed by kinematic observations of dense-gas tracers.


title: Newton vs the machine: solving the chaotic three-body problem using deep neural networks

abstract: Since its formulation by Sir Isaac Newton, the problem of solving the equations of motion for three bodies under their own gravitational force has remained practically unsolved. Currently, the solution for a given initialization can only be found by performing laborious iterative calculations that have unpredictable and potentially infinite computational cost, due to the system's chaotic nature. We show that an ensemble of solutions obtained using an arbitrarily precise numerical integrator can be used to train a deep artificial neural network (ANN) that, over a bounded time interval, provides accurate solutions at fixed computational cost and up to 100 million times faster than a state-of-the-art solver. Our results provide evidence that, for computationally challenging regions of phase-space, a trained ANN can replace existing numerical solvers, enabling fast and scalable simulations of many-body systems to shed light on outstanding phenomena such as the formation of black-hole binary systems or the origin of the core collapse in dense star clusters.

categories: astro-ph.GA astro-ph.SR cs.LG physics.comp-ph


title: A manifestly covariant framework for causal set dynamics

abstract: We propose a manifestly covariant framework for causal set dynamics. The framework is based on a structure, dubbed covtree, which is a partial order on certain sets of finite, unlabeled causal sets. We show that every infinite path in covtree corresponds to at least one infinite, unlabeled causal set. We show that transition probabilities for a classical random walk on covtree induce a classical measure on the $\sigma$-algebra generated by the stem sets.


title: Magnetic Nanoparticle Chains in Gelatin Ferrogels: Bioinspiration from Magnetotactic Bacteria

abstract: Inspired by chains of ferrimagnetic nanocrystals (NCs) in magnetotactic bacteria (MTB), the synthesis and detailed characterization of ferrimagnetic magnetite NC chain-like assemblies is reported. An easy green synthesis route in a thermoreversible gelatin hydrogel matrix is used. The structure of these magnetite chains prepared with and without gelatin is characterized by means of transmission electron microscopy, including electron tomography (ET). These structures indeed bear resemblance to the magnetite assemblies found in MTB, known for their mechanical flexibility and outstanding magnetic properties and known to crystallographically align their magnetite NCs along the strongest <111> magnetization easy axis. Using electron holography (EH) and angular dependent magnetic measurements, the magnetic interaction between the NCs and the generation of a magnetically anisotropic material can be shown. The electro- and magnetostatic modeling demonstrates that in order to precisely determine the magnetization (by means of EH) inside chain-like NCs assemblies, their exact shape, arrangement and stray-fields have to be considered (ideally obtained using ET).


title: Reinforcement Learning for Robotic Manipulation using Simulated Locomotion Demonstrations

abstract: Mastering robotic manipulation skills through reinforcement learning (RL) typically requires the design of shaped reward functions. Recent developments in this area have demonstrated that using sparse rewards, i.e. rewarding the agent only when the task has been successfully completed, can lead to better policies. However, state-action space exploration is more difficult in this case. Recent RL approaches to learning with sparse rewards have leveraged high-quality human demonstrations for the task, but these can be costly, time consuming or even impossible to obtain. In this paper, we propose a novel and effective approach that does not require human demonstrations. We observe that every robotic manipulation task could be seen as involving a locomotion task from the perspective of the object being manipulated, i.e. the object could learn how to reach a target state on its own. In order to exploit this idea, we introduce a framework whereby an object locomotion policy is initially obtained using a realistic physics simulator. This policy is then used to generate auxiliary rewards, called simulated locomotion demonstration rewards (SLDRs), which enable us to learn the robot manipulation policy. The proposed approach has been evaluated on 13 tasks of increasing complexity, and can achieve higher success rate and faster learning rates compared to alternative algorithms. SLDRs are especially beneficial for tasks like multi-object stacking and non-rigid object manipulation.


title: Towards Resolving Propensity Contradiction in Offline Recommender Learning

abstract: We study offline recommender learning from explicit rating feedback in the presence of selection bias. A current promising solution for the bias is the inverse propensity score (IPS) estimation. However, the performance of existing propensity-based methods can suffer significantly from the propensity estimation bias. In fact, most of the previous IPS-based methods require some amount of missing-completely-at-random (MCAR) data to accurately estimate the propensity. This leads to a critical self-contradiction; IPS is ineffective without MCAR data, even though it originally aims to learn recommenders from only missing-not-at-random feedback. To resolve this propensity contradiction, we derive a propensity-independent generalization error bound and propose a novel algorithm to minimize the theoretical bound via adversarial learning. Our theory and algorithm do not require a propensity estimation procedure, thereby leading to a well-performing rating predictor without the true propensity information. Extensive experiments demonstrate that the proposed approach is superior to a range of existing methods both in rating prediction and ranking metrics in practical settings without MCAR data.


title: Engel elements in weakly branch groups

abstract: We study properties of Engel elements in weakly branch groups, lying in the group of automorphisms of a spherically homogeneous rooted tree. More precisely, we prove that the set of bounded left Engel elements is always trivial in weakly branch groups. In the case of branch groups, the existence of non-trivial left Engel elements implies that these are all $p$-elements and that the group is virtually a $p$-group (and so periodic) for some prime $p$. We also show that the set of right Engel elements of a weakly branch group is trivial under a relatively mild condition. Also, we apply these results to well-known families of weakly branch groups, like the multi-GGS groups.


title: Waves of space-time from a collapsing compact object

abstract: We study the partial time dependent collapse of a spherically symmetric compact object with initial mass $M_1+M_2$ and final mass $M_2$ and the waves of space-time emitted during the collapse via back-reaction effects. We obtain exact analytical solutions for the waves of space-time in an example in which $M_1=M_2=(M_1+M_2)/2$. The wavelengths of the space-time emitted waves during the collapse have the cut (we use natural units $c=\hbar=1$): $\lambda < (2/b)$, $(1/b)$-being the time scale that describes the decay of the compact object.


title: Intelligence in Strategic Games

abstract: The article considers strategies of coalitions that are based on intelligence information about moves of some of the other agents. The main technical result is a sound and complete logical system that describes the interplay between coalition power modality with intelligence and distributed knowledge modality in games with imperfect information.


title: On the complexity of acyclic modules in automata networks

abstract: Modules were introduced as an extension of Boolean automata networks. They have inputs which are used in the computation said modules perform, and can be used to wire modules with each other. In the present paper we extend this new formalism and study the specific case of acyclic modules. These modules prove to be well described in their limit behavior by functions called output functions. We provide other results that offer an upper bound on the number of attractors in an acyclic module when wired recursively into an automata network, alongside a diversity of complexity results around the difficulty of deciding the existence of cycles depending on the number of inputs and the size of said cycle.


title: An interferometric method for determining the losses of spatially multi-mode nonlinear waveguides based on second harmonic generation

abstract: The characterisation of loss in optical waveguides is essential in understanding the performance of these devices and their limitations. Whilst interferometric-based methods generally provide the best results for low-loss waveguides, they are almost exclusively used to provide characterization in cases where the waveguide is spatially single-mode. Here, we introduce a Fabry-P\'erot-based scheme to estimate the losses of a nonlinear (birefringent or quasi-phase matched) waveguide at a wavelength where it is multi-mode. The method involves measuring the generated second harmonic power as the pump wavelength is scanned over the phasematching region. Furthermore, it is shown that this method allows one to infer the losses of different second harmonic spatial modes by scanning the pump field over the separated phase matching spectra. By fitting the measured phasematching spectra from a titanium indiffused lithium niobate waveguide sample to the model presented in this paper, it is shown that one can estimate the second harmonic losses of a single spatial-mode, at wavelengths where the waveguides are spatially multi-mode.


title: Gevrey regularity for a system coupling the Navier-Stokes system with a beam: the non-flat case

abstract: We consider a bi-dimensional viscous incompressible fluid in interaction with a beam located at its boundary. We show the existence of strong solutions for this fluid-structure interaction system, extending a previous result where we supposed that the initial deformation of the beam was small. The main point of the proof consists in the study of the linearized system and in particular in proving that the corresponding semigroup is of Gevrey class.


title: Progress on the ISIS synchrotron digital low level RF system upgrade

abstract: The ISIS synchrotron at the Rutherford Appleton Laboratory in the UK now routinely uses a dual harmonic RF system to accelerate beam currents in excess of 230uA to run two target stations simultaneously. The acceleration in the ISIS synchrotron is provided by six fundamental frequency (1RF) and four second harmonic (2RF) RF cavities. The 1RF systems are required to sweep from 1.3MHz to 3.1MHz during the 10ms acceleration period, repeated at 50Hz, with the 2RF systems sweeping from 2.6MHz to 6.3MHz. The existing analogue LLRF control system has been in service for over 30 years and is now showing some signs of old age and spare parts are becoming difficult to source. In order to overcome this and to give more stable control of the phase of the RF voltage at each of the cavities, changes have been made to the LLRF control system. A new FPGA based combined frequency law generator . master oscillator has been implemented using 'off the shelf' National Instruments PXI-express based FlexRIO modules. This initial design has been successfully used during the ISIS operational cycles for over three years. This paper reports on the commissioning of the FlexRIO system, the implementation and recent testing of the cavity control loops and plans for the gradual replacement of remaining parts of the LLRF system.


title: Filter List Generation for Underserved Regions

abstract: Filter lists play a large and growing role in protecting and assisting web users. The vast majority of popular filter lists are crowd-sourced, where a large number of people manually label resources related to undesirable web resources (e.g., ads, trackers, paywall libraries), so that they can be blocked by browsers and extensions. Because only a small percentage of web users participate in the generation of filter lists, a crowd-sourcing strategy works well for blocking either uncommon resources that appear on "popular" websites, or resources that appear on a large number of "unpopular" websites. A crowd-sourcing strategy will perform poorly for parts of the web with small "crowds", such as regions of the web serving languages with (relatively) few speakers. This work addresses this problem through the combination of two novel techniques: (i) deep browser instrumentation that allows for the accurate generation of request chains, in a way that is robust in situations that confuse existing measurement techniques, and (ii) an ad classifier that uniquely combines perceptual and page-context features to remain accurate across multiple languages. We apply our unique two-step filter list generation pipeline to three regions of the web that currently have poorly maintained filter lists: Sri Lanka, Hungary, and Albania. We generate new filter lists that complement existing filter lists. Our complementary lists block an additional 3,349 of ad and ad-related resources (1,771 unique) when applied to 6,475 pages targeting these three regions. We hope that this work can be part of an increased effort at ensuring that the security, privacy, and performance benefits of web resource blocking can be shared with all users, and not only those in dominant linguistic or economic regions.


title: Stabilization of a rigid body moving in a compressible viscous fluid

abstract: We consider the stabilizability of a fluid-structure interaction system where the fluid is viscous and compressible and the structure is a rigid ball. The feedback control of the system acts on the ball and corresponds to a force that would be produced by a spring and a damper connecting the center of the ball to a fixed point $h_1$. We prove the global-in-time existence of strong solutions for the corresponding system under a smallness condition on the initial velocities and on the distance between the initial position of the center of the ball and $h_1$. Then, we show with our feedback law, that the fluid and the structure velocities go to 0 and that the center of the ball goes to $h_1$ as $t\to \infty$.


title: Algorithms and Hardness Results for the Maximum Balanced Connected Subgraph Problem

abstract: The Balanced Connected Subgraph problem (BCS) was recently introduced by Bhore et al. (CALDAM 2019). In this problem, we are given a graph $G$ whose vertices are colored by red or blue. The goal is to find a maximum connected subgraph of $G$ having the same number of blue vertices and red vertices. They showed that this problem is NP-hard even on planar graphs, bipartite graphs, and chordal graphs. They also gave some positive results: BCS can be solved in $O(n^3)$ time for trees and $O(n + m)$ time for split graphs and properly colored bipartite graphs, where $n$ is the number of vertices and $m$ is the number of edges. In this paper, we show that BCS can be solved in $O(n^2)$ time for trees and $O(n^3)$ time for interval graphs. The former result can be extended to bounded treewidth graphs. We also consider a weighted version of BCS (WBCS). We prove that this variant is weakly NP-hard even on star graphs and strongly NP-hard even on split graphs and properly colored bipartite graphs, whereas the unweighted counterpart is tractable on those graph classes. Finally, we consider an exact exponential-time algorithm for general graphs. We show that BCS can be solved in $2^{n/2}n^{O(1)}$ time. This algorithm is based on a variant of Dreyfus-Wagner algorithm for the Steiner tree problem.


title: Linking scalar elastodynamics and non-Hermitian quantum mechanics

abstract: Recent years have seen a fascinating pollination of ideas from quantum theories to elastodynamics---a theory that phenomenologically describes the time-dependent macroscopic response of materials. Here, we open route to transfer additional tools from non-Hermitian quantum mechanics. We begin by identifying the differences and similarities between the one-dimensional elastodynamics equation and the time-independent Schrodinger equation, and finding the condition under which the two are equivalent. Subsequently, we demonstrate the application of the non-Hermitian perturbation theory to determine the response of elastic systems; calculation of leaky modes and energy decay rate in heterogenous solids with open boundaries using a quantum mechanics approach; and construction of degeneracies in the spectrum of these assemblies. The latter result is of technological importance, as it introduces an approach to harness extraordinary wave phenomena associated with non-Hermitian degeneracies for practical devices, by designing them in simple elastic systems. As an example of such application, we demonstrate how an assembly of elastic slabs that is designed with two degenerate shear states according to our scheme, can be used for mass sensing with enhanced sensitivity by exploiting the unique topology near the exceptional point of degeneracy.


title: Force balance controls the relaxation time of the gradient descent algorithm in the satisfiable phase

abstract: We numerically study the relaxation dynamics of the single layer perceptron with the spherical constraint. This is the simplest model of neural networks and serves a prototypical mean-field model of both convex and non-convex optimization problems. The relaxation time of the gradient descent algorithm rapidly increases near the SAT-UNSAT transition point. We numerically confirm that the first non-zero eigenvalue of the Hessian controls the relaxation time. This first eigenvalue vanishes much faster upon approaching the SAT-UNSAT transition point than the prediction of Marchenko-Pastur law in random matrix theory derived under the assumption that the set of unsatisfied constraints are uncorrelated. This leads to a non-trivial critical exponent of the relaxation time in the SAT phase. Using a simple scaling analysis, we show that the isolation of this first eigenvalue from the bulk of spectrum is attributed to the force balance at the SAT-UNSAT transition point. Finally, we show that the estimated critical exponent of the relaxation time in the non-convex region agrees very well with that of frictionless spherical particles, which have been studied in the context of the jamming transition of granular materials.


title: Positivity of chromatic symmetric functions associated with Hessenberg functions of bounce number 3

abstract: We give a proof of Stanley-Stembridge conjecture on chromatic symmetric functions for the class of all unit interval graphs with independence number 3. That is, we show that the chromatic symmetric function of the incomparability graph of a unit interval order in which the length of a chain is at most 3 is positively expanded as a linear sum of elementary symmetric functions.


title: Stable quantum dots in an InSb two-dimensional electron gas

abstract: Indium antimonide (InSb) two-dimensional electron gases (2DEGs) have a unique combination of material properties: high electron mobility, strong spin-orbit interaction, large Land\'{e} g-factor, and small effective mass. This makes them an attractive platform to explore a variety of mesoscopic phenomena ranging from spintronics to topological superconductivity. However, there exist limited studies of quantum confined systems in these 2DEGs, often attributed to charge instabilities and gate drifts. We overcome this by removing the $\delta$-doping layer from the heterostructure, and induce carriers electrostatically. This allows us to perform the first detailed study of stable gate-defined quantum dots in InSb 2DEGs. We demonstrate two distinct strategies for carrier confinement and study the charge stability of the dots. The small effective mass results in a relatively large single particle spacing, allowing for the observation of an even-odd variation in the addition energy. By tracking the Coulomb oscillations in a parallel magnetic field we determine the ground state spin configuration and show that the large g-factor ($\sim$30) results in a singlet-triplet transition at magnetic fields as low as 0.3 T.


title: Linear chirp instability analysis for ultrafast pulse metrology

abstract: Pulse train instabilities have often given rise to confusion in misinterpretation in ultrafast pulse characterization measurements. Most prominently known as the coherent artifact, a partially mode-locked laser with non-periodic waveform may still produce an autocorrelation that has often been misinterpreted as indication for a coherent pulse train. Some modern pulse characterization methods easily miss the presence of a coherent artifact, too. Here we address the particularly difficult situation of a pulse train with chirp-only instability. This instability is shown to be virtually invisible to autocorrelation measurements, but can be detected with FROG, SPIDER, and dispersion scan. Our findings clearly show that great care is necessary to rule out a chirp instability in lasers with unclear mode-locking mechanism and in compression experiments in the single-cycle regime. Among all dynamical pulse train instabilities analyzed so far, this instability appears to be the best hidden incoherence and is most difficult to detect.


title: Toward a robust inference method for the likelihood of low-luminosity gamma-ray bursts to be progenitors of ultrahigh-energy cosmic rays correlating with starburst galaxies

abstract: Very recently, the Pierre Auger Collaboration reported a $4.5\sigma$ correlation between the arrival directions of the highest energy cosmic rays and nearby starburst galaxies. The cosmic rays producing the anisotropy signal have been proposed to originate in low-luminosity gamma-ray bursts (llGRBs). On the basis of the well-justified assumption that at redshift $z<0.3$ the host metallicity is a good indicator of the llGRB production rate, we show that the association of llGRBs and the starbursts correlating with Auger data is excluded at the 90% confidence level.


title: Throughput and Delay Analysis of Slotted Aloha with Batch Service

abstract: In this paper, we study the throughput and delay performances of the slotted Aloha with batch service, which has wide applications in random access networks. Different from the classical slotted Aloha, each node in the slotted Aloha with batch service can transmit up to M packets once it succeeds in channel competition. The throughput is substantially improved because up to M packets jointly undertake the overhead due to contention. In an innovative vacation model developed in this paper, we consider each batch of data transmission as a busy period of each node, and the process between two successive busy periods as a vacation period. We then formulate the number of arrivals during a vacation period in a renewal-type equation, which characterizes the dependency between busy periods and vacation periods. Based on this formulation, we derive the mean waiting time of a packet and the bounded delay region for the slotted Aloha with batch service. Our results indicate the throughput and delay performances are substantially improved with the increase of batch sizeM, and the bounded delay region is enlarged accordingly. As M goes to infinity, we find the saturated throughput can approach 100% of channel capacity, and the system remains stable irrespective of the population size and transmission probability.


title: Determination of the second virial coefficient of interacting Bosons using the method of Wigner distribution function

abstract: In a previous article \cite{kn:anirban1} a method has been introduced to derive the all order Bose-Einstein distribution of the non interacting Bosons as the solution of the Wigner equation. The process was a perturbative one where the Bose-Einstein distribution was taken as the unperturbed solution. In this article it is shown that the same formalism is also applicable in the case of interacting Bosons. The formalism has been applied to calculate the quantum second virial coefficient of the Bosons interacting pairwise via Lenard-Jones potential and compared with the previous result.


title: The NEXT experiment for neutrinoless double beta decay searches

abstract: The Neutrino Experiment with a Xenon TPC (NEXT) seeks to discover the neutrinoless double beta ($\beta\beta 0\nu$) decay of $^{136}$Xe using a high-pressure xenon gas time projection chamber with electroluminescent amplification. The observation of $\beta\beta 0\nu$ decay would prove that neutrinos are Majorana particles and that lepton number is not conserved. The NEXT detector concept offers several features of great value for $\beta\beta 0\nu$-decay searches, including excellent energy resolution, tracking for the active suppression of backgrounds and scalability to large source masses. The initial phase of the NEXT project (2009-2014) was devoted to R&D with two prototypes (DEMO and DBDM) of approximately 1 kg of active xenon mass that demonstrated the performance of the detector concept. During the second phase of the project (2015-2019), the NEXT Collaboration has operated underground at the Laboratorio Subterr\'aneo de Canfranc, in Spain, a radio-pure detector of about 5 kg of xenon mass. The goal of the current phase is the construction, commissioning and operation of the NEXT-100 detector, with a predicted sensitivity to the $\beta\beta0\nu$ decay half-life of $6\times10^{25}$ years (90% CL) after a run of 3 years. The Collaboration is planning as well a future tonne-scale phase to explore $\beta\beta 0\nu$-decay half-lives beyond $10^{26}$ years.


title: A space-time hybridizable discontinuous Galerkin method for linear free-surface waves

abstract: We present and analyze a novel space-time hybridizable discontinuous Galerkin (HDG) method for the linear free-surface problem on prismatic space-time meshes. We consider a mixed formulation which immediately allows us to compute the velocity of the fluid. In order to show well-posedness, our space-time HDG formulation makes use of weighted inner products. We perform an a priori error analysis in which the dependence on the time step and spatial mesh size is explicit. We provide two numerical examples: one that verifies our analysis and a wave maker simulation.


title: Evaporating Planets in SNe Ia

abstract: I suggest that in rare cases type Ia supernovae (SNe Ia) evaporate Jupiter-like planets that reside within a=50Ro from the SN, accounting for the presence of ~0.001Mo of hydrogen in the ejecta. It is not easy to form such planets. The common envelope evolution (CEE) might end with the formation of a circumbinary disk with a typical size of R(disk)~10-50Ro. However, the bright central star might prevent the formation of planets. Alternatively, the planets might be first-generation planets that survived the entire evolution. Interestingly, all five SN Ia binary scenarios (core degenerate; double degenerate; double detonation; single degenerate; WD-WD collision) might account for the presence of planets.


title: Machine-learning computation of distance modulus for local galaxies

abstract: Quickly growing computing facilities and an increasing number of extragalactic observations encourage the application of data-driven approaches to uncover hidden relations from astronomical data. In this work we raise the problem of distance reconstruction for a large number of galaxies from available extensive observations. We propose a new data-driven approach for computing distance moduli for local galaxies based on the machine-learning regression as an alternative to physically oriented methods. We use key observable parameters for a large number of galaxies as input explanatory variables for training: magnitudes in U, B, I, and K bands, corresponding colour indices, surface brightness, angular size, radial velocity, and coordinates. We performed detailed tests of the five machine-learning regression techniques for inference of $m-M$: linear, polynomial, k-nearest neighbours, gradient boosting, and artificial neural network regression. As a test set we selected 91 760 galaxies at $z<0.2$ from the NASA/IPAC extragalactic database with distance moduli measured by different independent redshift methods. We find that the most effective and precise is the neural network regression model with two hidden layers. The obtained root-mean-square error of 0.35 mag, which corresponds to a relative error of 16\%, does not depend on the distance to galaxy and is comparable with methods based on the Tully-Fisher and Fundamental Plane relations. The proposed model shows a 0.44 mag (20\%) error in the case of spectroscopic redshift absence and is complementary to existing photometric redshift methodologies. Our approach has great potential for obtaining distance moduli for around 250 000 galaxies at $z<0.2$ for which the above-mentioned parameters are already observed.


title: Scattering on square lattice from crack with damage zone

abstract: A semi-infinite crack in infinite square lattice is subjected to a wave coming from infinity, thereby leading to its scattering by the crack surfaces. A partially damaged zone ahead of the crack-tip is modeled by an arbitrarily distributed stiffness of the damaged links. While the open crack, with an atomically sharp crack-tip, in the lattice has been solved in closed form with help of {the} scalar Wiener-Hopf formulation (SIAM Journal on Applied Mathematics, 75, 1171--1192; 1915--1940), the problem considered here becomes very intricate depending on the nature of damaged links. For instance, in the case of partially bridged finite zone it involves a $2\times2$ matrix kernel of formidable class. But using an original technique, the problem, including the general case of arbitrarily damaged links, is reduced to a scalar one with the exception that it involves solving an auxiliary linear system of $N \times N$ equations where $N$ defines the length of the damage zone. The proposed method does allow, effectively, the construction of an exact solution. Numerical examples and the asymptotic approximation of the scattered field far away from the crack-tip are also presented.

categories: physics.class-ph cond-mat.mes-hall math-ph math.MP


title: Quantitative level lowering for Galois representations

abstract: We use Galois cohomology methods to produce optimal mod $p^d$ level lowering congruences to a $p$-adic Galois representation that we construct as a well chosen lift of a given residual mod $p$ representation. Using our explicit Galois cohomology methods, we construct for a reductive group $G$ and a given residual representation $\bar{\rho}: \Gamma_F \to G(k)$, ramified at a finite set of primes $S$, in favorable conditions that we identify, a finite set of lifts $\rho$, $\{\rho^q\}$ of $\bar{\rho}$ to $G(W(k))$ with the following properties: $\rho: \Gamma_F \to G(W(k))$ is ramified precisely at $S \cup Q$, with $Q$ a finite set of primes disjoint from $S$. For $q \in Q$, $\rho^q:G_F \to G(W(k))$ is unramified outside $S \cup Q \backslash \{q\}$ and $\rho$ and $\rho^q$ are congruent mod $p^d$ if $\rho$ mod $p^d$ is unramified at $q$. Furthermore, the Galois representations $\{\rho^q\}$ are "independent".


title: The Blessings of Multiple Causes: A Reply to Ogburn et al. (2019)

abstract: Ogburn et al. (2019, arXiv:1910.05438) discuss "The Blessings of Multiple Causes" (Wang and Blei, 2018, arXiv:1805.06826). Many of their remarks are interesting. But they also claim that the paper has "foundational errors" and that its "premise is...incorrect." These claims are not substantiated. There are no foundational errors; the premise is correct.


title: On $t$-relaxed 2-distant circular coloring of graphs

abstract: Let $k$ be an positive integer. For any two integers $i$ and $j$ in $\{0,1,\dots,k-1\}$, let $|i-j|_k=\min\{|i-j|,k-|i-j|\}$ be the circular distance between $i$ and $j$. Let $t$ be a nonnegative integer. Suppose $f$ is a mapping from $V(G)$ to $\{0,1,\dots,k-1\}$. If adjacent vertices receive different integers, and for each vertex $u$ of $G$, the number of neighbors $v$ of $u$ with $|f(u)-f(v)|_k=1$ is at most $t$, then $f$ is called a $t$-relaxed 2-distant circular $k$-coloring, or simply a $(\frac{k}{2},t)^*$-coloring of $G$. If $G$ has a $(\frac{k}{2},t)^*$-coloring, then $G$ is called $(\frac{k}{2},t)^*$-colorable. In this paper, we prove that, for any two fixed integers $k$ and $t$ with $k\geq2$ and $t\geq1$, deciding whether $G$ is $(\frac{k}{2},t)^*$-colorable is NP-complete expect the case $k=2$ and the case $k=3$ and $t\leq3$, which are polynomially solvable. For any outerplanar graph $G$, e show that all outerplanar graphs are $(\frac{5}{2},4)^*$-colorable, we prove that there is no fixed positive integer $t$ such that all outerplanar graphs are $(\frac{4}{2},t)^*$-colorable.


title: An Adaptive Broadcasting Strategy for Efficient Dynamic Mapping in Vehicular Networks

abstract: In this work, we face the issue of achieving an efficient dynamic mapping in vehicular networking scenarios, i.e., to obtain an accurate estimate of the positions and trajectories of connected vehicles in a certain area. State of the art solutions are based on the periodic broadcasting of the position information of the network nodes, with an inter-transmission period set by a congestion control scheme. However, the movements and maneuvers of vehicles can often be erratic, making transmitted data inaccurate or downright misleading. To address this problem, we propose to adopt a dynamic transmission scheme based on the actual positioning error, sending new data when the estimate passes a preset error threshold. Furthermore, the proposed method adapts the error threshold to the operational context according to a congestion control algorithm that limits the collision probability among broadcast packet transmissions. This threshold-based strategy can reduce the network load by avoiding the transmission of redundant messages, and is shown to improve the overall positioning accuracy by more than 20% in realistic urban scenarios.


title: Lead2Gold: Towards exploiting the full potential of noisy transcriptions for speech recognition

abstract: The transcriptions used to train an Automatic Speech Recognition (ASR) system may contain errors. Usually, either a quality control stage discards transcriptions with too many errors, or the noisy transcriptions are used as is. We introduce Lead2Gold, a method to train an ASR system that exploits the full potential of noisy transcriptions. Based on a noise model of transcription errors, Lead2Gold searches for better transcriptions of the training data with a beam search that takes this noise model into account. The beam search is differentiable and does not require a forced alignment step, thus the whole system is trained end-to-end. Lead2Gold can be viewed as a new loss function that can be used on top of any sequence-to-sequence deep neural network. We conduct proof-of-concept experiments on noisy transcriptions generated from letter corruptions with different noise levels. We show that Lead2Gold obtains a better ASR accuracy than a competitive baseline which does not account for the (artificially-introduced) transcription noise.


title: Bose-like few-fermion systems

abstract: Dealing with a few-fermion system in the canonical ensemble, rather than in the grand canonical ensemble, shows that a few-fermion system with odd number fermions behaves differently from a few-fermion system with even number fermions. An even-number-fermion system behaves like a Bose system rather than a Fermi system.


title: Multivariate Forecasting Evaluation: On Sensitive and Strictly Proper Scoring Rules

abstract: In recent years, probabilistic forecasting is an emerging topic, which is why there is a growing need of suitable methods for the evaluation of multivariate predictions. We analyze the sensitivity of the most common scoring rules, especially regarding quality of the forecasted dependency structures. Additionally, we propose scoring rules based on the copula, which uniquely describes the dependency structure for every probability distribution with continuous marginal distributions. Efficient estimation of the considered scoring rules and evaluation methods such as the Diebold-Mariano test are discussed. In detailed simulation studies, we compare the performance of the renowned scoring rules and the ones we propose. Besides extended synthetic studies based on recently published results we also consider a real data example. We find that the energy score, which is probably the most widely used multivariate scoring rule, performs comparably well in detecting forecast errors, also regarding dependencies. This contradicts other studies. The results also show that a proposed copula score provides very strong distinction between models with correct and incorrect dependency structure. We close with a comprehensive discussion on the proposed methodology.


title: Real Lines on Random Cubic Surfaces

abstract: We give an explicit formula for the expectation of the number of real lines on a random invariant cubic surface, i.e. a surface $Z\subset \mathbb{R}P^3$ defined by a random gaussian polynomial whose probability distribution is invariant under the action of the orthogonal group $O(4)$ by change of variables. Such invariant distributions are completely described by one parameter $\lambda\in [0,1]$ and as a function of this parameter the expected number of real lines equals: \begin{equation} E_\lambda=\frac{9(8\lambda^2+(1-\lambda)^2)}{2\lambda^2+(1-\lambda)^2}\left(\frac{2\lambda^2}{8\lambda^2+(1-\lambda)^2}-\frac{1}{3}+\frac{2}{3}\sqrt{\frac{8\lambda^2+(1-\lambda)^2}{20\lambda^2+(1-\lambda)^2}}\right). \end{equation} This result generalizes previous results by Basu, Lerario, Lundberg and Peterson for the case of a Kostlan polynomial, which corresponds to $\lambda=\frac{1}{3}$ and for which $E_{\frac{1}{3}}=6\sqrt{2}-3.$ Moreover, we show that the expectation of the number of real lines is maximized by random purely harmonic cubic polynomials, which corresponds to the case $\lambda=1$ and for which $E_1=24\sqrt{\frac{2}{5}}-3$.


title: Blade Products and the Angle Bivector of Subspaces

abstract: Principal angles are used to define an angle bivector of subspaces, which fully describes their relative inclination. Its exponential is related to the Clifford geometric product of blades, gives rotors connecting subspaces via minimal geodesics in Grassmannians, and decomposes giving Pl\"ucker coordinates, projection factors and angles with various subspaces. This leads to new geometric interpretations for this product and its properties, and to formulas relating other blade products (scalar, inner, outer, etc., including those of Grassmann algebra) to angles between subspaces. Contractions are linked to an asymmetric angle, while commutators and anticommutators involve hyperbolic functions of the angle bivector, shedding new light on their properties.


title: Segmentation Criteria in the Problem of Porosity Determination based on CT Scans

abstract: Porous materials are widely used in different applications, in particular they are used to create various filters. Their quality depends on parameters that characterize the internal structure such as porosity, permeability and so on. Computed tomography (CT) allows one to see the internal structure of a porous object without destroying it. The result of tomography is a gray image. To evaluate the desired parameters, the image should be segmented. Traditional intensity threshold approaches did not reliably produce correct results due to limitations with CT images quality. Errors in the evaluation of characteristics of porous materials based on segmented images can lead to the incorrect estimation of their quality and consequently to the impossibility of exploitation, financial losses and even to accidents. It is difficult to perform correctly segmentation due to the strong difference in voxel intensities of the reconstructed object and the presence of noise. Image filtering as a preprocessing procedure is used to improve the quality of segmentation. Nevertheless, there is a problem of choosing an optimal filter. In this work, a method for selecting an optimal filter based on attributive indicator of porous objects (should be free from 'levitating stones' inside of pores) is proposed. In this paper, we use real data where beam hardening artifacts are removed, which allows us to focus on the noise reduction process


title: On a Hybrid Version of the Vinogradov Mean Value Theorem

abstract: Given a family $\varphi= (\varphi_1, \ldots, \varphi_d)\in \mathbb{Z}[T]^d$ of $d$ distinct nonconstant polynomials, a positive integer $k\le d$ and a real positive parameter $\rho$, we consider the mean value $$ M_{k, \rho} (\varphi, N) = \int_{\mathbf{x} \in [0,1]^k} \sup_{\mathbf{y} \in [0,1]^{d-k}} \left| S_{\varphi}(\mathbf{x}, \mathbf{y}; N) \right|^\rho d\mathbf{x} $$ of exponential sums $$ S_{\varphi}( \mathbf{x}, \mathbf{y}; N) = \sum_{n=1}^{N} \exp\left(2 \pi i \left(\sum_{j=1}^k x_j \varphi_j(n)+ \sum_{j=1}^{d-k}y_j\varphi_{k+j}(n)\right)\right), $$ where $\mathbf{x} = (x_1, \ldots, x_k)$ and $\mathbf{y} =(y_1, \ldots, y_{d-k})$. The case of polynomials $\varphi_i(T) = T^i$, $i =1, \ldots, d$ and $k=d$ corresponds to the classical Vinaogradov mean value theorem. Here motivated by recent works of Wooley (2015) and the authors (2019) on bounds on $\sup_{\mathbf{y} \in [0,1]^{d-k}} \left| S_{\varphi}( \mathbf{x}, \mathbf{y}; N) \right|$ for almost all $\mathbf{x} \in [0,1]^k$, we obtain nontrivial bounds on $M_{k, \rho} (\varphi, N)$.


title: Linear Hyperelliptic Hodge Integrals

abstract: We provide a closed form expression for linear Hodge integrals on the hyperelliptic locus. Specifically, we find a succinct combinatorial formula for all intersection numbers on the hyperelliptic locus with one $\lambda$-class, and powers of a $\psi$-class pulled back along the branch map. This is achieved by using Atiyah-Bott localization on a stack of stable maps into the orbifold $[\mathbb{P}^1/\mathbb{Z}_2]$.


title: A Generalized and Robust Method Towards Practical Gaze Estimation on Smart Phone

abstract: Gaze estimation for ordinary smart phone, e.g. estimating where the user is looking at on the phone screen, can be applied in various applications. However, the widely used appearance-based CNN methods still have two issues for practical adoption. First, due to the limited dataset, gaze estimation is very likely to suffer from over-fitting, leading to poor accuracy at run time. Second, the current methods are usually not robust, i.e. their prediction results having notable jitters even when the user is performing gaze fixation, which degrades user experience greatly. For the first issue, we propose a new tolerant and talented (TAT) training scheme, which is an iterative random knowledge distillation framework enhanced with cosine similarity pruning and aligned orthogonal initialization. The knowledge distillation is a tolerant teaching process providing diverse and informative supervision. The enhanced pruning and initialization is a talented learning process prompting the network to escape from the local minima and re-born from a better start. For the second issue, we define a new metric to measure the robustness of gaze estimator, and propose an adversarial training based Disturbance with Ordinal loss (DwO) method to improve it. The experimental results show that our TAT method achieves state-of-the-art performance on GazeCapture dataset, and that our DwO method improves the robustness while keeping comparable accuracy.


title: What Did Your Adversary Believe? Optimal Filtering and Smoothing in Counter-Adversarial Autonomous Systems

abstract: We consider fixed-interval smoothing problems for counter-adversarial autonomous systems. An adversary deploys an autonomous filtering and control system that i) measures our current state via a noisy sensor, ii) computes a posterior estimate (belief) and iii) takes an action that we can observe. Based on such observed actions and our knowledge of our state sequence, we aim to estimate the adversary's past and current beliefs -- this forms a foundation for predicting, and counteracting against, future actions. We derive the optimal smoother for the adversary's beliefs (we treat the problem in a Bayesian framework). Moreover, we demonstrate how the smoother can be computed for discrete systems even though the corresponding backward variables do not admit a finite-dimensional characterization. Finally, we illustrate our results in numerical simulations.