title: Introduction to Gauge Theory of Gravitation

abstract: The fundamental interactions of nature, the electroweak and the quantum chromodynamics, are described in the Standard Model by the Gauge Theory under internal symmetries that maintain the invariance of the functional action. The fundamental interaction of gravitation is very well described by Einstein's General Relativity in a Riemannian spacetime metric, but General Relativity has been over time a gravitational field theory apart from the Standard Model. The theory of Gauge allows under symmetries of the group of Poincar\'e to impose invariances in the functional of the action of the spinor field that result in the gravitational interaction with the fermions. In this approach the gravitational field, besides being described by the equation similar to General Relativity, also brings a spin-gravitational interaction in a Riemann-Cartan spacetime.


title: Zero-Shot Knowledge Distillation in Deep Networks

abstract: Knowledge distillation deals with the problem of training a smaller model (Student) from a high capacity source model (Teacher) so as to retain most of its performance. Existing approaches use either the training data or meta-data extracted from it in order to train the Student. However, accessing the dataset on which the Teacher has been trained may not always be feasible if the dataset is very large or it poses privacy or safety concerns (e.g., bio-metric or medical data). Hence, in this paper, we propose a novel data-free method to train the Student from the Teacher. Without even using any meta-data, we synthesize the Data Impressions from the complex Teacher model and utilize these as surrogates for the original training data samples to transfer its learning to Student via knowledge distillation. We, therefore, dub our method "Zero-Shot Knowledge Distillation" and demonstrate that our framework results in competitive generalization performance as achieved by distillation using the actual training data samples on multiple benchmark datasets.


title: Framework for a novel mixed analytical/numerical approach for the computation of two-loop $N$-point Feynman diagrams

abstract: A framework to represent and compute two-loop $N$-point Feynman diagrams as double-integrals is discussed. The integrands are 'generalised one-loop type" multi-point functions multiplied by simple weighting factors. The final integrations over these two variables are to be performed numerically, whereas the ingredients involved in the integrands, in particular the "generalised one-loop type" functions, are computed analytically. The idea is illustrated on a few examples of scalar three- and four-point functions.


title: Computer Assisted Proofs of Attracting Invariant Tori for ODEs

abstract: This work studies existence and regularity questions for attracting invariant tori in three dimensional dissipative systems of ordinary differential equations. Our main result is a constructive method of computer assisted proof which applies to explicit problems in non-perturbative regimes. We obtain verifiable lower bounds on the regularity of the attractor in terms of the ratio of the expansion rate on the torus with the contraction rate near the torus. We consider separately two important cases of rotational and resonant tori. In the rotational case we obtain $C^k$ lower bounds on the regularity of the embedding. In the resonant case we verify the existence of tori which are only $C^0$ and neither star-shaped nor Lipschitz


title: The syzygy theorem for B\'ezout rings

abstract: We provide constructive versions of Hilbert's syzygy theorem for Z and Z/nZ following Schreyer's method. Moreover, we extend these results to arbitrary coherent strict B\'ezout rings with a divisibility test for the case of finitely generated modules whose module of leading terms is finitely generated.


title: The correspondence formula of Dolbeault complex on pair deformation

abstract: Given a holomorphic family of pairs $\{(X_t,E_t)\}$, where each $E_t$ is holomorphic vector bundle over compact complex manifold $X_t$. For small enough $t$, we get a correspondence between the Dolbeault complex of $E_t$-valued $(p,q)$-forms on $X_t$ and the one of $E_0$-valued $(p,q)$-forms on $X_0$.


title: Continual Learning in Deep Neural Network by Using a Kalman Optimiser

abstract: Learning and adapting to new distributions or learning new tasks sequentially without forgetting the previously learned knowledge is a challenging phenomenon in continual learning models. Most of the conventional deep learning models are not capable of learning new tasks sequentially in one model without forgetting the previously learned ones. We address this issue by using a Kalman Optimiser. The Kalman Optimiser divides the neural network into two parts: the long-term and short-term memory units. The long-term memory unit is used to remember the learned tasks and the short-term memory unit is to adapt to the new task. We have evaluated our method on MNIST, CIFAR10, CIFAR100 datasets and compare our results with state-of-the-art baseline models. The results show that our approach enables the model to continually learn and adapt to the new changes without forgetting the previously learned tasks.


title: A combinatorial approach for the state complexity of the Shuffle product

abstract: We investigate the state complexity of the shuffle operation on regular languages initiated by Campeanu et al. and studied subsequently by Brzozowski et al. We shift the problem into the combinatorics domain by turning the problem of state accessibility into a problem of intersection of partitions. This allows us to develop new tools and to reformulate the conjecture of Brzozowski et al. about the above-mentionned state complexity.


title: Quasilinear elliptic equations with sub-natural growth terms and nonlinear potential theory

abstract: We discuss recent advances in the theory of quasilinear equations of the type $ -\Delta_{p} u = \sigma u^{q} \; \; \text{in} \;\; \mathbb{R}^n, $ in the case $0<q< p-1$, where $\sigma$ is a nonnegative measurable function, or measure, for the $p$-Laplacian $\Delta_{p}u= \text{div}(|\nabla u|^{p-2}\nabla u)$, as well as more general quasilinear, fractional Laplacian, and Hessian operators. Within this context, we obtain some new results, in particular, necessary and sufficient conditions for the existence of solutions $u \in \text{BMO}(\mathbb{R}^n)$, $u \in L^r_{{\rm loc}}(\mathbb{R}^n)$, etc., and prove an enhanced version of Wolff's inequality for intrinsic nonlinear potentials associated with such problems.


title: Non-Parametric Estimation of Spot Covariance Matrix with High-Frequency Data

abstract: Estimating spot covariance is an important issue to study, especially with the increasing availability of high-frequency financial data. We study the estimation of spot covariance using a kernel method for high-frequency data. In particular, we consider first the kernel weighted version of realized covariance estimator for the price process governed by a continuous multivariate semimartingale. Next, we extend it to the threshold kernel estimator of the spot covariances when the underlying price process is a discontinuous multivariate semimartingale with finite activity jumps. We derive the asymptotic distribution of the estimators for both fixed and shrinking bandwidth. The estimator in a setting with jumps has the same rate of convergence as the estimator for diffusion processes without jumps. A simulation study examines the finite sample properties of the estimators. In addition, we study an application of the estimator in the context of covariance forecasting. We discover that the forecasting model with our estimator outperforms a benchmark model in the literature.


title: Sharp results concerning disjoint cross-intersecting families

abstract: For an $n$-element set $X$ let $\binom{X}{k}$ be the collection of all its $k$-subsets. Two families of sets $\mathcal A$ and $\mathcal B$ are called cross-intersecting if $A\cap B \neq \emptyset$ holds for all $A\in\mathcal A$, $B\in\mathcal B$. Let $f(n,k)$ denote the maximum of $\min\{|\mathcal A|, |\mathcal B|\}$ where the maximum is taken over all pairs of {\em disjoint}, cross-intersecting families $\mathcal A, \mathcal B\subset\binom{[n]}{k}$. Let $c=\log_2e$. We prove that $f(n,k)=\left\lfloor\frac12\binom{n-1}{k-1}\right\rfloor$ essentially iff $n>ck^2$ (cf. Theorem~1.4 for the exact statement). Let $f^*(n,k)$ denote the same maximum under the additional restriction that the intersection of all members of both $\mathcal A$ and $\mathcal B$ are empty. For $k\ge5$ and $n\ge k^3$ we show that $f^*(n,k)=\left\lfloor\frac12\left(\binom{n-1}{k-1}-\binom{n-2k}{k-1}\right)\right\rfloor+1$ and the restriction on $n$ is essentially sharp (cf. Theorem~5.4).


title: A blow-up criterion for the strong solutions to the nonhomogeneous Navier-Stokes-Korteweg equations in dimension three

abstract: This paper proves a Serrin's type blow-up criterion for the 3D density-dependent Navier-Stokes-Korteweg equations with vacuum. It is shown that if the density and velocity field satisfy some Serrin's type condition, then the strong solutions to the density-dependent Navier-Stokes-Korteweg equations can exist globally.


title: The transverse-momentum subtraction method at N$^3$LO applied to Higgs boson production at the LHC

abstract: We consider the extension of the transverse-momentum ($q_T$) subtraction method at next-to-next-to-next-to-leading order (N$^3$LO) in perturbative QCD. While all the $q_T$-subtraction ingredients at $q_T \neq 0$ are known in analytical form, the third-order collinear functions and helicity-flip functions, which contribute only at $q_T=0$, are approximated using a prescription which uses the known result for the total Higgs boson cross section at this order. As a first application of the third-order $q_T$-subtraction method, we present the N$^3$LO rapidity distribution of the Higgs boson at the LHC.


title: Can Bio-Inspired Swarm Algorithms Scale to Modern Societal Problems

abstract: Taking inspiration from nature for meta-heuristics has proven popular and relatively successful. Many are inspired by the collective intelligence exhibited by insects, fish and birds. However, there is a question over their scalability to the types of complex problems experienced in the modern world. Natural systems evolved to solve simpler problems effectively, replicating these processes for complex problems may suffer from inefficiencies. Several causal factors can impact scalability; computational complexity, memory requirements or pure problem intractability. Supporting evidence is provided using a case study in Ant Colony Optimisation (ACO) regards tackling increasingly complex real-world fleet optimisation problems. This paper hypothesizes that contrary to common intuition, bio-inspired collective intelligence techniques by their very nature exhibit poor scalability in cases of high dimensionality when large degrees of decision making are required. Facilitating scaling of bio-inspired algorithms necessitates reducing this decision making. To support this hypothesis, an enhanced Partial-ACO technique is presented which effectively reduces ant decision making. Reducing the decision making required by ants by up to 90% results in markedly improved effectiveness and reduced runtimes for increasingly complex fleet optimisation problems. Reductions in traversal timings of 40-50% are achieved for problems with up to 45 vehicles and 437 jobs.


title: Subcubic Equivalences Between Graph Centrality Measures and Complementary Problems

abstract: Despite persistent efforts, there is no known technique for obtaining unconditional super-linear lower bounds for the computational complexity of the problems in P. Vassilevska Williams and Williams introduce a fruitful approach to advance a better understanding of the computational complexity of the problems in P. In particular, they consider All Pairs Shortest Paths (APSP) and other fundamental problems such as checking whether a matrix defines a metric, verifying the correctness of a matrix product, and detecting a negative triangle in a graph. Abboud, Grandoni, and Vassilevska Williams study well-known graph centrality problems such as Radius, Median, etc., and make a connection between their computational complexity to that of two fundamental problems, namely APSP and Diameter. They show any algorithm with subcubic running time for these centrality problems, implies a subcubic algorithm for either APSP or Diameter. In this paper, we define vertex versions for these centrality problems and based on that we introduce new complementary problems. The main open problem of Abboud et al. is whether or not APSP and Diameter are equivalent under subcubic reduction. One of the results of this paper is APSP and CoDiameter, which is the complementary version of Diameter, are equivalent. Moreover, for some of the problems in this set, we show that they are equivalent to their complementary versions. Considering the slight difference between a problem and its complementary version, these equivalences give us the impression that every problem has such a property, and thus APSP and Diameter are equivalent. This paper is a step forward in showing a subcubic equivalence between APSP and Diameter, and we hope that the approach introduced in our paper can be helpful to make this breakthrough happen.


title: Fragile Topology Protected by Inversion Symmetry: Diagnosis, Bulk-Boundary Correspondence, and Wilson Loop

abstract: We study the bulk and boundary properties of fragile topological insulators (TIs) protected by inversion symmetry, mostly focusing on the class A of the Altland-Zirnbauer classification. First, we propose an efficient method for diagnosing fragile band topology by using the symmetry data in momentum space. Using this method, we show that among all the possible parity configurations of inversion-symmetric insulators, at least 17 percent of them have fragile topology in 2D while fragile TIs are less than 3 percent in 3D. Second, we study the bulk-boundary correspondence of fragile TIs protected by inversion symmetry. In particular, we generalize the notion of $d$-dimensional ($d$D) $k$th-order TIs, which is normally defined for $0<k\le d$, to the cases with $k>d$, and show that they all have fragile topology. In terms of the Dirac Hamiltonian, a $d$D $k$th-order TI has $(k-1)$ boundary mass terms. We show that a minimal fragile TI with the filling anomaly can be considered as the $d$D $(d+1)$th-order TI, and all the other $d$D $k$th-order TIs with $k>(d+1)$ can be constructed by stacking $d$D $(d+1)$th-order TIs. Although $d$D $(d+1)$th-order TIs have no in-gap states, the boundary mass terms carry an odd winding number along the boundary, which induces localized charges on the boundary at the positions where the boundary mass terms change abruptly. In the cases with $k>(d+1)$, we show that the net parity of the system with boundaries can distinguish topological insulators and trivial insulators. Also, by studying the (nested) Wilson loop spectra, we determine the minimal number of bands to resolve the Wannier obstruction, which is consistent with the prediction from our diagnosis method of fragile topology. Finally, we show that a $(d+1)$D $(k-1)$th-order TI can be obtained by an adiabatic pumping of $d$D $k$th-order TI, which generalizes the previous study of the 2D 3rd-order TI.


title: On the existence of Hopf bifurcations in the sequential and distributive double phosphorylation cycle

abstract: Protein phosphorylation cycles are important mechanisms of the post translational modification of a protein and as such an integral part of intracellular signaling and control. We consider the sequential phosphorylation and dephosphorylation of a protein at two binding sites. While it is known that proteins where phosphorylation is processive and dephosphorylation is distributive admit oscillations (for some value of the rate constants and total concentrations) it is not known whether or not this is the case if both phosphorylation and dephosphorylation are distributive. We study four simplified mass action models of sequential and distributive phosphorylation and show that for each of those there do not exist rate constants and total concentrations where a Hopf bifurcation occurs. To arrive at this result we use convex parameters to parameterize the steady state and Hurwitz matrices.


title: Toward Operator-to-Waveform 5G Radio Access Network Slicing

abstract: Radio access network (RAN) slicing realizes a vision where physical network resources that belong to a specific infrastructure provider can be shared among multiple mobile network operators (MNOs). Existing work in this area has addressed RAN slicing at different levels of network abstractions, but has often neglected the multitude of tightly intertwined inter-level operations involved in real-world slicing systems. For this reason, this article discusses a novel framework for operator-to-waveform 5G RAN slicing. In the proposed framework, slicing operations are treated holistically, including MNO's selection of base stations (BSs) and maximum number of users, down to the waveform-level scheduling of resource blocks. Experimental results show that the proposed framework provides up to 150% improvement in terms of number of resource blocks that can be used to enable 5G transmission technologies that require coordination and synchronization among BSs.


title: Longest common substring for random subshifts of finite type

abstract: In this paper, we study the behaviour of the longest common substring for random subshifts of finite type (for dynamicists) or of the longest common substring for random sequences in random environments (for probabilists). We prove that, under some exponential mixing assumptions, this behaviour is linked to the R\'enyi entropy of the stationary measure. We emphasize that what we establish is a quenched result.


title: Exoplanetary Monte Carlo Radiative Transfer with Correlated-k I. Benchmarking Transit and Emission Observables

abstract: Current observational data of exoplanets are providing increasing detail of their 3D atmospheric structures. As characterisation efforts expand in scope, the need to develop consistent 3D radiative-transfer methods becomes more pertinent as the complex atmospheric properties of exoplanets are required to be modelled together consistently. We aim to compare the transmission and emission spectra results of a 3D Monte Carlo Radiative Transfer (MCRT) model to contemporary radiative-transfer suites. We perform several benchmarking tests of a MCRT code, Cloudy Monte Carlo Radiative Transfer (CMCRT), to transmission and emission spectra model output. We add flexibility to the model through the use of k-distribution tables as input opacities. We present a hybrid MCRT and ray tracing methodology for the calculation of transmission spectra with a multiple scattering component. CMCRT compares well to the transmission spectra benchmarks at the 10s of ppm level. Emission spectra benchmarks are consistent to within 10% of the 1D models. We suggest that differences in the benchmark results are likely caused by geometric effects between plane-parallel and spherical models. In a practical application, we post-process a cloudy 3DHD 189733b GCM model and compare to available observational data. Our results suggest the core methodology and algorithms of CMCRT produce consistent results to contemporary radiative transfer suites. 3D MCRT methods are highly suitable for detailed post-processing of cloudy and non-cloudy 1D and 3D exoplanet atmosphere simulations in instances where atmospheric inhomogeneities, significant limb effects/geometry or multiple scattering components are important considerations.


title: The spiral pattern rotation speed of the Galaxy and the corotation radius with GAIA DR2

abstract: In this work we revisit the issue of the rotation speed of the spiral arms and the location of the corotation radius of our Galaxy. This research was performed using homogeneous data set of young open clusters (age < 50 Myr) determined from Gaia DR2 data. The stellar astrometric membership were determined using proper motions and parallaxes, taking into account the full covariance matrix. The distance, age, reddening and metallicity of the clusters were determined by our non subjective multidimensional global optimization tool to fit theoretical isochrones to Gaia DR2 photometric data. The rotation speed of the arms is obtained from the relation between age and angular distance of the birthplace of the clusters to the present-day position of the arms. Using the clusters belonging to the Sagittarius-Carina, Local and Perseus arms, and adopting the Galactic parameters $R_0$ = 8.3 kpc and $V_0$ = 240 km\,s$^{-1}$, we determine a pattern speed of $28.2 \pm 2.1$ km\,s$^{-1}$\,kpc$^{-1}$, with no difference between the arms. This implies that the corotation radius is $R_c = 8.51 \pm 0.64$ kpc, close to the solar Galactic orbit ($R_c/R_0 = 1.02\pm0.07$).


title: Micromagnetic Modelling of the Heat Assisted Switching Process in high Anisotropy FePt Granular Thin Films

abstract: The dynamic process of assisted magnetic switchins has been simulated to investigate the associated physics. The model uses a Voronoi construction to determine the physical structure of the nano granular thin film recording media; and the Landau-Lifshitz-Bloch (LLB) equation is solved to evolve the magnetic system in time. The reduction of the magnetization is determined over a range of peak system temperatures and for a number of anisotropy values. The results show that the HAMR process is not simply magnetization reversal over a thermally reduced energy barrier. To achieve full magnetization reversal (for all anisotropies investigated) an applied field strength of at least 6kOe is required and the peak system temperature must reach at least the Curie point (Tc). When heated to Tc the magnetization associated with each grain is destroyed, which invokes the non-precessional linear reversal mode. Reversing the magnetization through this linear reversal mode is favourable, as the reversal time is two orders of magnitude smaller than that associated with precession. Under these conditions, as the temperature decreases to ambient, the magnetization recovers in the direction of the applied field, completing the reversal process. Also the model produces results which are consistent with the concept of thermal writability; when heating the media to Tc, the smaller grains require a larger field strength to reverse the magnetization.


title: On the regularity of D-modules generated by relative characters

abstract: Following the ideas of Ginzburg, for a subgroup $K$ of a connected reductive $\mathbb{R}$-group $G$ we introduce the notion of $K$-admissible $D$-modules on a homogeneous $G$-variety $Z$. We show that $K$-admissible $D$-modules are regular holonomic when $K$ and $Z$ are absolutely spherical. This framework includes: (i) the relative characters attached to two spherical subgroups $H_1$ and $H_2$, provided that the twisting character $\chi_i$ factors through the maximal reductive quotient of $H_i$, for $i = 1, 2$; (ii) localization on $Z$ of Harish-Chandra modules; (iii) the generalized matrix coefficients when $K(\mathbb{R})$ is maximal compact. This complements the holonomicity proven by Aizenbud--Gourevitch--Minchenko. The use of regularity is illustrated by a crude estimate on the growth of $K$-admissible distributions which based on tools from subanalytic geometry.


title: Characteristic polynomials for random band matrices near the threshold

abstract: The paper continues previous works which study the behavior of second correlation function of characteristic polynomials of the special case of $n\times n$ one-dimensional Gaussian Hermitian random band matrices, when the covariance of the elements is determined by the matrix $J=(-W^2\triangle+1)^{-1}$. Applying the transfer matrix approach, we study the case when the bandwidth $W$ is proportional to the threshold $\sqrt{n}$


title: Derivation of the core mass -- halo mass relation of fermionic and bosonic dark matter halos from an effective thermodynamical model

abstract: We consider the possibility that dark matter halos are made of quantum particles such as fermions or bosons in the form of Bose-Einstein condensates. In that case, they generically have a "core-halo" structure with a quantum core that depends on the type of particle considered and a halo that is relatively independent of the dark matter particle and that is similar to the NFW profile of cold dark matter. We model the halo by an isothermal gas with an effective temperature $T$. We then derive the core mass -- halo mass relation $M_c(M_v)$ of dark matter halos from an effective thermodynamical model by extremizing the free energy $F(M_c)$ with respect to the core mass $M_c$. We obtain a general relation that is equivalent to the "velocity dispersion tracing" relation according to which the velocity dispersion in the core $v_c^2\sim GM_c/R_c$ is of the same order as the velocity dispersion in the halo $v_v^2\sim GM_v/r_v$. We provide therefore a justification of this relation from thermodynamical arguments. In the case of fermions, we obtain a relation $M_c\propto M_v^{1/2}$ that agrees with the relation found numerically by Ruffini {\it et al.}. In the case of noninteracting bosons, we obtain a relation $M_c\propto M_v^{1/3}$ that agrees with the relation found numerically by Schive {\it et al.}. In the case of bosons with a repulsive self-interaction in the Thomas-Fermi limit, we predict a relation $M_c\propto M_v^{2/3}$ that still has to be confirmed numerically. We also obtain a general approximate core mass -- halo mass relation that is valid for bosons with arbitrary repulsive or attractive self-interaction. For an attractive self-interaction, we determine the maximum halo mass that can harbor a stable quantum core (dilute axion "star").


abstract: During the last decades, we have witnessed a surge of interests of learning a low-dimensional space with discriminative information from one single view. Even though most of them can achieve satisfactory performance in some certain situations, they fail to fully consider the information from multiple views which are highly relevant but sometimes look different from each other. Besides, correlations between features from multiple views always vary greatly, which challenges multi-view subspace learning. Therefore, how to learn an appropriate subspace which can maintain valuable information from multi-view features is of vital importance but challenging. To tackle this problem, this paper proposes a novel multi-view dimension reduction method named Multi-view Locality Low-rank Embedding for Dimension Reduction (MvL2E). MvL2E makes full use of correlations between multi-view features by adopting low-rank representations. Meanwhile, it aims to maintain the correlations and construct a suitable manifold space to capture the low-dimensional embedding for multi-view features. A centroid based scheme is designed to force multiple views to learn from each other. And an iterative alternating strategy is developed to obtain the optimal solution of MvL2E. The proposed method is evaluated on 5 benchmark datasets. Comprehensive experiments show that our proposed MvL2E can achieve comparable performance with previous approaches proposed in recent literatures.


title: Self-Supervised Similarity Learning for Digital Pathology

abstract: Using features extracted from networks pretrained on ImageNet is a common practice in applications of deep learning for digital pathology. However it presents the downside of missing domain specific image information. In digital pathology, supervised training data is expensive and difficult to collect. We propose a self-supervised method for feature extraction by similarity learning on whole slide images (WSI) that is simple to implement and allows creation of robust and compact image descriptors. We train a siamese network, exploiting image spatial continuity and assuming spatially adjacent tiles in the image are more similar to each other than distant tiles. Our network outputs feature vectors of length 128, which allows dramatically lower memory storage and faster processing than networks pretrained on ImageNet. We apply the method on digital pathology WSIs from the Camelyon16 train set and assess and compare our method by measuring image retrieval of tumor tiles and descriptor pair distance ratio for distant/near tiles in the Camelyon16 test set. We show that our method yields better retrieval task results than existing ImageNet based and generic self-supervised feature extraction methods. To the best of our knowledge, this is also the first published method for self-supervised learning tailored for digital pathology.


title: On tau-functions for the Toda lattice hierarchy

abstract: We extend a recent result of [13] for the KdV hierarchy to the Toda lattice hierarchy. Namely, for an arbitrary solution to the Toda lattice hierarchy, we define a pair of wave functions, and use them to give explicit formulae for the generating series of $k$-point correlation functions of the solution. Applications to computing GUE correlators and Gromov--Witten invariants of the Riemann sphere are under consideration.


title: The s-process enriched star HD 55496: origin from a globular cluster or from the tidal disruption of a dwarf galaxy?

abstract: We present a new abundance analysis of HD 55496, previously known as a metal-poor barium star. We found that HD 55496 has a metallicity [Fe/H] $=-1.55$ and is s-process enriched. We find that HD 55496 presents four chemical peculiarities: (i) a Na-O abundance anti-correlation; (ii) it is aluminum rich; (iii) it is carbon poor for a s-process enriched star and (iv) the heavy 2$^{nd}$ s-process peak elements, such as Ba, La, Ce, and Nd, present smaller abundances than the lighter s-process elements, such as Sr, Y and Zr, which is not usually observed among the chemically-peculiar binary stars at this metallicity. The heavy-element abundance pattern suggests that the main source of the neutrons is the $^{22}$Ne($\alpha$,n)$^{25}$Mg reaction. Taken all these abundance evidence together into consideration, this strongly suggests that HD 55496 is a "second generation of globular cluster star" formed from gas already strongly enriched in s-process elements and now is a field halo object. Our dynamical analysis, however, indicates that the past encounter probabilities with the known globular clusters is very small ($\leq 6\%$). This evidence, together with the fact of having a retrograde motion, points to a halo intruder possibly originated from the tidal disruption of a dwarf galaxy.


title: Improved convergence analysis of Lasserre's measure-based upper bounds for polynomial minimization on compact sets

abstract: We consider the problem of computing the minimum value $f_{\min,K}$ of a polynomial $f$ over a compact set $K \subseteq \mathbb{R}^n$, which can be reformulated as finding a probability measure $\nu$ on $K$ minimizing $\int_K f d\nu$. Lasserre showed that it suffices to consider such measures of the form $\nu = q\mu$, where $q$ is a sum-of-squares polynomial and $\mu$ is a given Borel measure supported on $K$. By bounding the degree of $q$ by $2r$ one gets a converging hierarchy of upper bounds $f^{(r)}$ for $f_{\min,K}$. When $K$ is the hypercube $[-1, 1]^n$, equipped with the Chebyshev measure, the parameters $f^{(r)}$ are known to converge to $f_{\min,K}$ at a rate in $O(1/r^2)$. We extend this error estimate to a wider class of convex bodies, while also allowing for a broader class of reference measures, including the Lebesgue measure. Our analysis applies to simplices, balls and convex bodies that locally look like a ball. In addition, we show an error estimate in $O(\log r / r)$ when $K$ satisfies a minor geometrical condition, and in $O(\log^2 r / r^2)$ when $K$ is a convex body, equipped with the Lebesgue measure. This improves upon the currently best known error estimates in $O(1 / \sqrt{r})$ and $O(1/r)$ for these two respective cases.


title: Strong-field photoelectron holography beyond the electric dipole approximation: a semiclassical analysis

abstract: Strong-field photoelectron holography denotes the interference of various electron paths in laser-induced ionization of atoms, leading to interference patterns in the final momentum distribution. For a quantitative description of holography beyond the electric dipole approximation and in the presence of the Coulomb potential, we develop a semiclassical model in which the initial conditions of outgoing electrons are set according to the beyond-dipole strong-field approximation for the ionization step. The phases associated with trajectories are evaluated following the prescription for semiclassical propagators. The validity of the method is confirmed by comparison to the numerical solution of the time-dependent Schr\"odinger equation in two spatial dimensions. The semiclassical model reproduces correctly the backward and forward shifts of the photoelectron momenta along the laser propagation axis that arise from beyond-dipole dynamics. The position of the central holographic interference fringe can be estimated already from a simplified Coulomb-free interference model that provides closed-form expressions for the beyond-dipole shifts. In three dimensions, Coulomb focusing causes a breakdown of the semiclassical model for final momenta with directions close to the polarization axis. We implement a beyond-dipole regularization procedure based on the concept of glory scattering, which was recently used to describe Coulomb focusing in the dipole approximation. While the position of the central maximum and higher-order fringes can already be obtained approximately by simpler semiclassical modeling, this glory model is able to predict the shape of the distribution in the close vicinity of the central maximum. The violation of the dipole approximation in holography should be observable with mid-infrared fields, for which the forward/backward shifts can be comparable with the fringe spacing.


title: Incongruent equipartitions of the plane

abstract: R. Nandakumar asked whether there is a tiling of the plane by pairwise incongruent triangles of equal area and equal perimeter. Recently a negative answer was given by Kupavskii, Pach and Tardos. Still one may ask for weaker versions of the problem, or for the analogue of this problem for quadrangles, pentagons, or hexagons. Several answers were given by the first author in a previous paper. Here we solve three further cases. In particular, our main result shows that there are vertex-to-vertex tilings by pairwise incongruent triangles of unit area and bounded perimeter.


title: Almost inner derivations of Lie algebras II

abstract: We continue the algebraic study of almost inner derivations of Lie algebras over a field of characteristic zero and determine these derivations for free nilpotent Lie algebras, for almost abelian Lie algebras, for Lie algebras whose solvable radical is abelian and for several classes of filiform nilpotent Lie algebras. We find a family of $n$-dimensional characteristically nilpotent filiform Lie algebras $\mathfrak{f}_n$, for all $n\ge 13$, all of whose derivations are almost inner. Finally we compare the almost inner derivations of Lie algebras considered over two different fields $K\supseteq k$ for a finite-dimensional field extension.


title: In-medium effect on the thermodynamics and transport coefficients in van der Waals hadron resonance gas

abstract: An extension of the van der Waals hadron resonance gas (VDWHRG) model which includes in-medium thermal modification of hadron masses, the TVDWHRG model, is considered in this paper. Based on the 2+1 flavor Polyakov Linear Sigma Model(PLSM) and the scaling mass rule for hadrons we obtain the temperature behavior of all hadron masses for different fixed baryon chemical potentials $\mu_{B}$. We calculate various thermodynamic observables at $\mu_{B}=0$ GeV in TVDWHRG model. An improved agreement with the lattice data by TVDWHRG model in the crossover region ($T\sim 0.16-0.19$ GeV) is observed as compared to those by VDWHRG and Ideal HRG (IHRG) models. We further discuss the effects of in-medium modification of hadron masses and VDW interactions on the transport coefficients such as shear viscosity ($\eta$), scaled thermal ($\lambda/T^{2}$) and electrical ($\sigma_{el}/T$ conductivities in IHRG model at different $\mu_{B}$, by utilizing quasi-particle kinetic theory with relaxation time approximation.


title: Statistical limit laws for hyperbolic groups

abstract: Using techniques from ergodic theory and symbolic dynamics, we derive statistical limit laws for real valued functions on hyperbolic groups. In particular, our results apply to convex cocompact group actions on $\text{CAT}(-1)$ spaces, and provide a precise statistical comparison between word length and displacement. After generalising our methods to the multidimensional setting, we prove that the abelianisation map satisfies a non-degenerate multidimensional central limit theorem. We also obtain local limit theorems for group homomorphisms and for the displacement function associated to certain actions.


title: OGLE-2018-BLG-0532Lb: Cold Neptune With Possible Jovian Sibling

abstract: We report the discovery of the planet OGLE-2018-BLG-0532Lb, with very obvious signatures in the light curve that lead to an estimate of the planet-host mass ratio $q=M_{\rm planet}/M_{\rm host}\simeq 1\times10^{-4}$. Although there are no obvious systematic residuals to this double-lens/single-source (2L1S) fit, we find that $\chi^2$ can be significantly improved by adding either a third lens (3L1S, $\Delta\chi^2=81$) or second source (2L2S, $\Delta\chi^2=65$) to the lens-source geometry. After thorough investigation, we conclude that we cannot decisively distinguish between these two scenarios and therefore focus on the robustly-detected planet. However, given the possible presence of a second planet, we investigate to what degree and with what probability such additional planets may affect seemingly single-planet light curves. Our best estimates for the properties of the lens star and the secure planet are: a host mass $M\sim 0.25\,M_\odot$, system distance $D_L\sim 1\,$kpc and planet mass $m_{p,1}= 8\,M_\oplus$ with projected separation $a_{1,\perp}=1.4\,$au. However, there is a relatively bright $I=18.6$ (and also relatively blue) star projected within $<50\,$mas of the lens, and if future high-resolution images show that this is coincident with the lens, then it is possible that it is the lens, in which case, the lens would be both more massive and more distant than the best-estimated values above.


title: Dependent subsets of embedded projective varieties

abstract: Let $X\subset \mathbb {P}^r$ be an integral and non-degenerate variety. Set $n:= \dim (X)$. Let $\rho (X)''$ be the maximal integer such that every zero-dimensional scheme $Z\subset X$ smoothable in $X$ is linearly independent. We prove that $X$ is linearly normal if $\rho (X)''\ge \lceil (r+2)/2\rceil$ and that $\rho (X)'' < 2\lceil (r+1)/(n+1)\rceil$, unless either $n=r$ or $X$ is a rational normal curve.


title: Pretty Modular Symmetric Encryption (PMSE), compact algorithm for ''embedded cryptography'' with quite low computational cost

abstract: Nowadays, the dataflux shared between IOT systems must be secured from 8-bits to 64-bits processors systems. Several symmetric cryptographic algorithm already exist such as AES (Advanced Encryption Standard), RC4, Blowfish, etc. In this work, we propose an 8-bits encryption algorithm, combining several ideas of standard symmetric encryption algorithms. The aim is to provide a efficient, modular and compact algorithm able to tend to truly random one-time-pad encryption quality even for large data flux. The algorithm combines the implementation of a divergent polynomial with variable coefficients for pseudo-random keys generation, variable bit swapping and bitwise operations on data, and the use of one or two passwords. The encryption has been evaluated statistically, tested for image encryption and compared with One-Time-Pad encryption on the same data. Three implementations have been tested respectively in C, Javascript and GNU Octave. Encryption time performances are compared with AES on a 8-bits architecture: the Arduino Uno (ATmega328) microcontroller. A new concept of self-decryptionable web encrypted object, called Blocksnet is also presented and tested as a potential application of the compact and embedded PMSE algorithm.


title: An $L^p$-comparison, $p\in (1,\infty)$, on the finite differences of a discrete harmonic function at the boundary of a discrete box

abstract: It is well-known that for a harmonic function $u$ defined on the unit ball of the $d$-dimensional Euclidean space, $d\geq 2$, the tangential and normal component of the gradient $\nabla u$ on the sphere are comparable by means of the $L^p$-norms, $p\in(1,\infty)$, up to multiplicative constants that depend only on $d,p$. This paper formulates and proves a discrete analogue of this result for discrete harmonic functions defined on a discrete box on the $d$-dimensional lattice with multiplicative constants that do not depend on the size of the box.


title: Stochastic Variance Reduction for Deep Q-learning

abstract: Recent advances in deep reinforcement learning have achieved human-level performance on a variety of real-world applications. However, the current algorithms still suffer from poor gradient estimation with excessive variance, resulting in unstable training and poor sample efficiency. In our paper, we proposed an innovative optimization strategy by utilizing stochastic variance reduced gradient (SVRG) techniques. With extensive experiments on Atari domain, our method outperforms the deep q-learning baselines on 18 out of 20 games.


title: Quadratic formulas for split quaternions

abstract: Unlike the Hamilton quaternion algebra, the split-quaternions contain nontrivial zero divisors. In general speaking, it is hard to find the solutions of equations in algebras containing zero divisor. In this paper, we manage to derive explicit formulas for computing the roots of $x^{2}+bx+c=0$ in split quaternion algebra.


title: Towards the Limits of Existence of Nuclear Structure: Observation and First spectroscopy of the Isotope $^{31}$K by measuring its three-proton Decay

abstract: The most-remote from stability isotope $^{31}$K, which is located four atomic mass units beyond the proton drip line, has been observed. It is unbound in respect to three-proton (${3p}$) emission, and its decays have been detected in flight by measuring trajectories of all decay products using micro-strip detectors. The $3p$-emission processes have been studied by means of angular correlations $^{28}$S+3$p$ and the respective decay vertexes. The energies of the previously-unknown ground and excited states of $^{31}$K have been determined. This provides its $3p$ separation-energy value $S_{3p}$ of $-4.6(2)$ MeV. Upper half-life limits of 10 ps of the observed $^{31}$K states have been derived from distributions of the measured decay vertexes.


title: Analysis of numerical methods for spectral fractional elliptic equations based on the best uniform rational approximation

abstract: Here we study theoretically and compare experimentally an efficient method for solving systems of algebraic equations, where the matrix comes from the discretization of a fractional diffusion operator. More specifically, we focus on matrices obtained from finite difference or finite element approximation of second order elliptic problems in $\mathbb R^d$, $d=1,2,3$. The proposed methods are based on the best uniform rational approximation (BURA) $r_{\alpha,k}(t)$ of $t^{\alpha}$ on $[0,1]$. Here $r_{\alpha,k}$ is a rational function of $t$ involving numerator and denominator polynomials of degree at most $k$. We show that the proposed method is exponentially convergent with respect to $k$ and has some attractive properties. First, it reduces the solution of the nonlocal system to solution of $k$ systems with matrix $(A +c_j I)$ and $c_j>0$, $j=1,2,\ldots,k$. Thus, good computational complexity can be achieved if fast solvers are available for such systems. Second, the original problem and its rational approximation in the finite difference case are positivity preserving. In the finite element case, positivity preserving results when mass lumping is employed under some mild conditions on the mesh. Further, we prove that in the mass lumping case, the scheme still leads to the expected rate of convergence, at times assuming additional regularity on the right hand side. Finally, we present comprehensive numerical experiments on a number of model problems for various $\alpha$ in one and two spatial dimensions. These illustrate the computational behavior of the proposed method and compare its accuracy and efficiency with that of other methods developed by Harizanov et. al. and Bonito and Pasciak.


title: Charged Lovelock black holes in the presence of dark fluid with a nonlinear equation of state

abstract: We consider a model that charged static spherically-symmetric black hole is surrounded by dark fluid with nonlinear equation of state $p_d=-B/\rho_d$. We find that the energy density of the dark fluid can be characterized by two parameters. The derivation of metric solution, as well as the calculation of black hole thermodynamical quantities as functions of horizon radius, are performed. Specially, in $D$-dimensional Einstein gravity and Gauss-Bonnet gravity cases, we plot the metric functions and corresponding thermodynamical quantities, such as mass, Hawking temperature and heat capacity, by varying the values of spacetime dimensions and dark fluid parameters. The effects of the dark fluid parameters on black hole solutions as well as on thermodynamical stability of black holes are discussed.


title: On the Relation between the Extended Supporting Hyperplane Algorithm and Kelley's Cutting Plane Algorithm

abstract: Recently, Kronqvist et al.~\cite{KronqvistLundellWesterlund2016} rediscovered the supporting hyperplane algorithm of Veinott~\cite{Veinott1967} and demonstrated its computational benefits for solving convex mixed-integer nonlinear programs. In this paper we derive the algorithm from a geometric point of view. This enables us to show that the supporting hyperplane algorithm is equivalent to Kelley's cutting plane algorithm~\cite{J.E.Kelley1960} applied to a particular reformulation of the problem. As a result, we extend the applicability of the supporting hyperplane algorithm to convex problems represented by general, not necessarily convex, differentiable functions that satisfy a mild condition.


title: Dual-frequency spin resonance spectroscopy of diamond nitrogen-vacancy centers in zero magnetic field

abstract: The methods for controlling spin states of negatively charged nitrogen-vacancy (NV) centers using microwave (MW) or radiofrequency (RF) excitation fields for electron spin and nuclear spin transitions are effective in strong magnetic fields where a level anti-crossing (LAC) occurs. A LAC can also occur at zero field in the presence of transverse strain or electric fields in the diamond crystal, leading to mixing of the spin states. In this paper, we investigate zero-field LAC of NV centers using dual-frequency excitation spectroscopy. Under RF modulation of the spin states, we observe sideband transitions and Autler-Townes splitting in the optically detected magnetic resonance (ODMR) spectra. Numerical simulations show that the splitting originates from Landau-Zener transition between electron spin |$\pm$1> states, which potentially provides a new way of manipulating NV center spin states in zero or weak magnetic field.


title: On the management fourth-order Schr\"{o}dinger-Hartree equation

abstract: We consider the Cauchy problem associated to the fourth-order nonlinear Schr\"{o}dinger-Hartree equation with variable dispersion coefficients. The variable dispersion coefficients are assumed to be continuous or periodic and piecewise constant in time functions. We prove local and global well-posedness results for initial data in $H^s$-spaces. We also analyze the scaling limit of fast dispersion management and the convergence to a model with averaged dispersions.


title: Interpretable Neural Predictions with Differentiable Binary Variables

abstract: The success of neural networks comes hand in hand with a desire for more interpretability. We focus on text classifiers and make them more interpretable by having them provide a justification, a rationale, for their predictions. We approach this problem by jointly training two neural network models: a latent model that selects a rationale (i.e. a short and informative part of the input text), and a classifier that learns from the words in the rationale alone. Previous work proposed to assign binary latent masks to input positions and to promote short selections via sparsity-inducing penalties such as L0 regularisation. We propose a latent model that mixes discrete and continuous behaviour allowing at the same time for binary selections and gradient-based training without REINFORCE. In our formulation, we can tractably compute the expected value of penalties such as L0, which allows us to directly optimise the model towards a pre-specified text selection rate. We show that our approach is competitive with previous work on rationale extraction, and explore further uses in attention mechanisms.


title: Superconvergence of ultra-weak discontinuous Galerkin methods for the linear Schr\"odinger equation in one dimension

abstract: We analyze the superconvergence properties of ultra-weak discontinuous Galerkin (UWDG) methods with various choices of flux parameters for one-dimensional linear Schr\"odinger equation. In our previous work [10], stability and optimal convergence rate are established for a large class of flux parameters. Depending on the flux choices and if the polynomial degree $k$ is even or odd, in this paper, we prove $2k$ or $(2k-1)$-th order superconvergence rate for cell averages and numerical flux of the function, as well as $(2k-1)$ or $(2k-2)$-th order for numerical flux of the derivative. In addition, we prove superconvergence of $(k+2)$ or $(k+3)$-th order of the DG solution towards a special projection. At a class of special points, the function values and the first and second order derivatives of the DG solution are superconvergent with order $k+2, k+1, k$, respectively. The proof relies on the correction function techniques initiated in [8], and applied to [6] for direct DG (DDG) methods for diffusion problems. Compared with [6], Schr\"odinger equation poses unique challenges for superconvergence proof because of the lack of the dissipation mechanism from the equation. One major highlight of our proof is that we introduce specially chosen test functions in the error equation and show the superconvergence of the second derivative and jump across the cell interfaces of the difference between numerical solution and projected exact solution. This technique was originally proposed in [12] and is essential to elevate the convergence order for our analysis. Finally, by negative norm estimates, we apply the post-processing technique and show that the accuracy of our scheme can be enhanced to order $2k.$ Theoretical results are verified by numerical experiments.


title: Equivalent symmetric kernels of determinantal point processes

abstract: Determinantal point processes are point processes whose correlation functions are given by determinants of matrices. The entries of these matrices are given by one fixed function of two variables, which is called the kernel of the point process. It is well-known that there are different kernels that induce the same correlation functions. We classify all the possible transformations of a kernel that leave the induced correlation functions invariant, restricting to the case of symmetric kernels.


title: A Characterization of Semiglobal, Practical, Asymptotic Stability for Gain-Parametrized Systems

abstract: We consider a general class of nonlinear, constrained, discrete-time systems whose dynamics are parametrized by a set of gains. We define the semiglobal, practical, asymptotic stability (SPAS) of compact sets for this class of systems, and we provide a Lyapunov characterization of such sets. A set A that is SPAS with respect to a given system need not be an attractor for that system. Relative to existing characterizations of similar qualitative behaviors, our SPAS theorem does not require the existence of an asymptotically stable attractor associated to a nominal counterpart of the given dynamics.


title: simTPM: User-centric TPM for Mobile Devices (Technical Report)

abstract: Trusted Platform Modules are valuable building blocks for security solutions and have also been recognized as beneficial for security on mobile platforms, like smartphones and tablets. However, strict space, cost, and power constraints of mobile devices prohibit an implementation as dedicated on-board chip and the incumbent implementations are software TPMs protected by Trusted Execution Environments. In this paper, we present simTPM, an alternative implementation of a mobile TPM based on the SIM card available in mobile platforms. We solve the technical challenge of implementing a TPM2.0 in the resource-constrained SIM card environment and integrate our simTPM into the secure boot chain of the ARM Trusted Firmware on a HiKey960 reference board. Most notably, we address the challenge of how a removable TPM can be bound to the host device's root of trust for measurement. As such, our solution not only provides a mobile TPM that avoids additional hardware while using a dedicated, strongly protected environment, but also offers promising synergies with co-existing TEE-based TPMs. In particular, simTPM offers a user-centric trusted module. Using performance benchmarks, we show that our simTPM has competitive speed with a reported TEE-based TPM and a hardware-based TPM.


title: Gradient Ascent for Active Exploration in Bandit Problems

abstract: We present a new algorithm based on an gradient ascent for a general Active Exploration bandit problem in the fixed confidence setting. This problem encompasses several well studied problems such that the Best Arm Identification or Thresholding Bandits. It consists of a new sampling rule based on an online lazy mirror ascent. We prove that this algorithm is asymptotically optimal and, most importantly, computationally efficient.


title: Time and momentum-resolved tunneling spectroscopy of pump-driven non-thermal excitations in Mott insulators

abstract: We present a computational technique to calculate time and momentum resolved non-equilibrium spectral density of correlated systems using a tunneling approach akin scanning tunneling spectroscopy. The important difference is that our probe is extended, basically a copy of the sample, allowing one to extract the momentum information of the excitations. We illustrate the method by measuring the spectrum of a Mott-insulating extended Hubbard chain after a sudden quench with the aid of time-dependent density matrix renormalization group (tDMRG) calculations. We demonstrate that the system realizes a non-thermal state that is an admixture of spin and charge density wave states, with corresponding signatures that are recognizable as in-gap sub-bands. In particular, we identify a band of excitons and one of stable anti-bound states at high energies that gains enhanced visibility after the pump. We do not appreciate noticeable relaxation within the time-scales considered, which is attributed to the lack of decay channels due to spin-charge separation. These ideas can be readily applied to study transient dynamics and spectral signatures of correlation-driven non-equilibrium processes.


title: On the Fractional Riemann-Liouville Integral of Gauss-Markov processes and applications

abstract: We investigate the stochastic processes obtained as the fractional Riemann-Liouville integral of order $\alpha \in (0,1)$ of Gauss-Markov processes. The general expressions of the mean, variance and covariance functions are given. Due to the central rule, for the fractional integral of standard Brownian motion and of the non-stationary/stationary Ornstein-Uhlenbeck processes, the covariance functions are carried out in closed-form. In order to clarify how the fractional order parameter $\alpha$ affects these functions, their numerical evaluations are shown and compared also with those of the corresponding processes obtained by ordinary Riemann integral. The results are useful for fractional neuronal models with long range memory dynamics and involving correlated input processes. The simulation of these fractional integrated processes can be performed starting from the obtained covariance functions. A suitable neuronal model is proposed. Graphical comparisons are provided and discussed.


title: Classical Solutions for a Class of Burgers Equation

abstract: In this paper we consider a class of Burgers equation. We propose a new method of investigation for existence of classical solutions.


title: ATAC: A Tool for Automating Timed Automata Construction

abstract: In this paper, we focus on the design and verification of timed automata (TA). We introduce a new method for assisting construction and verification of TA models along with a tool implementing the proposed method, i.e., ATAC: Automated Timed Automata Construction. Our method provides two main functionalities, i.e., construction of TA models from descriptions and generation of temporal logic queries from specifications. Both description and specification sentences shall follow our well-defined structured natural language definition. TA models constructed from descriptions and temporal logic queries generated from specifications can be imported to UPPAAL, a verification tool for TA models. The goal is to accelerate the design phase for real-time systems by assisting the construction and verification of a formal model. We believe ATAC can be useful especially during the initial phases of the design process and help designers to avoid erroneous models.


title: DARC: Differentiable ARchitecture Compression

abstract: In many learning situations, resources at inference time are significantly more constrained than resources at training time. This paper studies a general paradigm, called Differentiable ARchitecture Compression (DARC), that combines model compression and architecture search to learn models that are resource-efficient at inference time. Given a resource-intensive base architecture, DARC utilizes the training data to learn which sub-components can be replaced by cheaper alternatives. The high-level technique can be applied to any neural architecture, and we report experiments on state-of-the-art convolutional neural networks for image classification. For a WideResNet with $97.2\%$ accuracy on CIFAR-10, we improve single-sample inference speed by $2.28\times$ and memory footprint by $5.64\times$, with no accuracy loss. For a ResNet with $79.15\%$ Top1 accuracy on ImageNet, we improve batch inference speed by $1.29\times$ and memory footprint by $3.57\times$ with $1\%$ accuracy loss. We also give theoretical Rademacher complexity bounds in simplified cases, showing how DARC avoids overfitting despite over-parameterization.


title: Semi-Supervised Learning by Augmented Distribution Alignment

abstract: In this work, we propose a simple yet effective semi-supervised learning approach called Augmented Distribution Alignment. We reveal that an essential sampling bias exists in semi-supervised learning due to the limited number of labeled samples, which often leads to a considerable empirical distribution mismatch between labeled data and unlabeled data. To this end, we propose to align the empirical distributions of labeled and unlabeled data to alleviate the bias. On one hand, we adopt an adversarial training strategy to minimize the distribution distance between labeled and unlabeled data as inspired by domain adaptation works. On the other hand, to deal with the small sample size issue of labeled data, we also propose a simple interpolation strategy to generate pseudo training samples. Those two strategies can be easily implemented into existing deep neural networks. We demonstrate the effectiveness of our proposed approach on the benchmark SVHN and CIFAR10 datasets. Our code is available at \url{https://github.com/qinenergy/adanet}.


title: Large enhancement of the thermoelectric power factor in disordered materials through resonant scattering

abstract: In the search for more efficient thermoelectric materials, scientists have placed high hopes in the possibility of enhancing the power factor using resonant states. In this study, we investigate theoretically the effects of randomly distributed resonant impurities on the power factor. Using the Chebyshev Polynomial Green's Function method, we compute the electron transport properties for very large systems (10 million atoms) with an exact treatment of disorder. The introduction of resonant defects can lead to a large enhancement of the power factor together with a sign inversion in the Seebeck coefficient. This boost depends crucially on the position of the resonant peak, and on the interplay between elastic impurity scattering and inelastic processes. Strong electron-phonon or electron-electron scattering are found detrimental. Finally, the robustness of our results is examined in the case of anisotropic orbitals and two-dimensional confinement. Our findings are promising for the prospect of thermoelectric power generation.


title: On Lipschitz-like continuity of a class of set-valued mappings

abstract: We study set-valued mappings defined by solution sets of parametric systems of equalities and inequalities. We prove Lipschitz-like continuity of these mappings under relaxed constant rank constraint qualification.


title: Generalized orbital varieties for Mirkovic-Vybornov slices as affinizations of Mirkovic-Vilonen cycles

abstract: We show that generalized orbital varieties for Mirkovic-Vybornov slices can be indexed by semi-standard Young tableaux. We also check that the Mirkovic-Vybornov isomorphism sends generalized orbital varieties to (dense subsets of) Mirkovic-Vilonen cycles, such that the (combinatorial) Lusztig datum of a generalized orbital variety, which it inherits from its tableau, is equal to the (geometric) Lusztig datum of its MV cycle.


title: Interacting topological edge channels

abstract: Electrical currents in a quantum spin Hall insulator are confined to the boundary of the system. The charge carriers can be described as massless relativistic particles, whose spin and momentum are coupled to each other. While the helical character of those states is by now well established experimentally, it is a fundamental open question how those edge states interact with each other when brought in spatial proximity. We employ a topological quantum point contact to guide edge channels from opposite sides into a quasi-one-dimensional constriction, based on inverted HgTe quantum wells. Apart from the expected quantization in integer steps of $2 e^2/h$, we find a surprising additional plateau at $e^2/h$. We explain our observation by combining band structure calculations and repulsive electron-electron interaction effects captured within the Tomonaga-Luttinger liquid model. The present results may have direct implications for the study of one-dimensional helical electron quantum optics, Majorana- and potentially para-fermions.


title: Estimation of Physiological Motion Using Highly Accelerated Continuous 2D MRI

abstract: Patient motion is well-known for degrading image quality during medical imaging. Especially positron emission tomography (PET) is susceptible to motion due to its usually long scan times. In hybrid PET/MRI (magnetic resonance imaging), simultaneously acquired dynamic MRI data can be used to correct for motion. Usually, MRI model-based motion correction approaches are applied to the PET data. However, these approaches may fail for non-predictable, irregular motion. We propose a novel approach for the continuous and real-time tracking of motion using highly accelerated, dynamic MRI for an accurate motion estimation. For this purpose, a TurboFLASH sequence is utilized in single-shot mode with additional exploiting GRAPPA acceleration. Sampling frequency for one slice is up to 26 ms and 520 ms for one 3D volume of 20 coronal slices. Principal component analysis and a phase-sensitive resorting of slices is performed to restore temporal consistency of the volumes. Motion is estimated from these volumes using hyper-elastic registration. The approach is validated with the help of a dynamic thorax phantom as well as with eleven healthy volunteers. Phantom ground-truth data demonstrates that the approach produces an accurate motion estimation. Volunteer validation proves that the approach is also valid for different respiratory amplitudes including highly irregular breathing. The approach could be proved to be promising for a continuous PET motion correction.


title: An issue with the classification of the scalar-tensor theories of gravity

abstract: In the bibliography a certain confusion arises in what regards to the classification of the gravitational theories into scalar-tensor theories and general relativity with a scalar field either minimally or non-minimally coupled to matter. Higher-derivatives Horndeski and beyond Horndeski theories that at first sight do not look like scalar-tensor theories only add to the confusion. To further complicate things, the discussion on the physical equivalence of the different conformal frames in which a given scalar-tensor theory may be formulated, makes even harder to achieve a correct classification. In this paper we propose a specific criterion for an unambiguous identification of scalar-tensor theories and discuss its impact on the conformal transformations issue. The present discussion carries not only pedagogical but also scientific interest since an incorrect classification of a given theory as a scalar-tensor theory of gravity may lead to conceptual issues and to the consequent misunderstanding of its physical implications.


title: Partial Redundancy Elimination using Lazy Code Motion

abstract: Partial Redundancy Elimination (PRE) is a compiler optimization that eliminates expressions that are redundant on some but not necessarily all paths through a program. In this project, we implemented a PRE optimization pass in LLVM and measured results on a variety of applications. We chose PRE because it is a powerful technique that subsumes Common Subexpression Elimination (CSE) and Loop Invariant Code Motion (LICM), and hence has the potential to greatly improve performance.


title: Does the Adam-Gibbs relation hold in simulated supercooled liquids?

abstract: We perform stringent tests of thermodynamic theories of the glass transition over the experimentally relevant temperature regime for several simulated glass-formers. The swap Monte Carlo algorithm is used to estimate the configurational entropy and static point-to-set lengthscale, and careful extrapolations are used for the relaxation times. We first quantify the relation between configurational entropy and the point-to-set lengthscale in two and three dimensions. We then show that the Adam-Gibbs relation is generally violated in simulated models for the experimentally relevant time window. Collecting experimental data for several supercooled molecular liquids, we show that the same trends are observed experimentally. Deviations from the Adam-Gibbs relation remain compatible with random first order transition theory, and may account for the reported discrepancies between Kauzmann and Vogel-Fulcher-Tammann temperatures. Alternatively, they may also indicate that even near $T_g$ thermodynamics is not the only driving force for slow dynamics.


title: The Largest Crab Glitch and the Vortex Creep Model

abstract: The Crab pulsar displayed its largest glitch on 2017 November. An extended initial spin-up phase of this largest glitch was resolved, for the first time with high cadence of observations both in radio and X-rays on a time-scale of 2 days. A combination of crustquake and vortex unpinning models is invoked to account for the extended spin-up, magnitude and post-glitch relaxation characteristics of this glitch. We evaluate the extended spin-up followed by the familiar spin-down as due to the creep response to the initial induced inward motion of some vortex lines pinned to broken crustal plates moving inward towards the rotation axis, together with the common and familiar post-glitch creep response to the sudden outward motion of vortices unpinned at the glitch. Our analysis confirms that the number of unpinned vortices participating in glitches are similar in all Crab glitches, and within an order of magnitude in all glitches from all pulsars. This typical number of unpinned vortices is related to the broken plate size in quakes as triggers for vortex unpinning avalanches. The physical determinant of this universal broken plate size is in turn the critical strain angle in the neutron star crust. Occurrence of this largest Crab glitch after a relatively long inactive period is consistent with accumulation of the pinned vorticity to be tapped.


title: A Neural, Interactive-predictive System for Multimodal Sequence to Sequence Tasks

abstract: We present a demonstration of a neural interactive-predictive system for tackling multimodal sequence to sequence tasks. The system generates text predictions to different sequence to sequence tasks: machine translation, image and video captioning. These predictions are revised by a human agent, who introduces corrections in the form of characters. The system reacts to each correction, providing alternative hypotheses, compelling with the feedback provided by the user. The final objective is to reduce the human effort required during this correction process. This system is implemented following a client-server architecture. For accessing the system, we developed a website, which communicates with the neural model, hosted in a local server. From this website, the different tasks can be tackled following the interactive-predictive framework. We open-source all the code developed for building this system. The demonstration in hosted in http://casmacat.prhlt.upv.es/interactive-seq2seq.


title: Some properties of geodesic $(\alpha,E)$-preinvex functions on a Riemannian manifold

abstract: In this article, we have introduced the concept of \textit{geodesic $(\alpha,E)$-invex set} and by using this concept the notion of \textit{geodesic $(\alpha,E)$-preinvex functions} and \textit{geodesic $(\alpha,E)$-invex functions} are developed on a Riemannian manifold. Moreover, several properties and results are deduced within aforesaid functions. An example is also constructed to illustrate the definition of geodesic $(\alpha,E)$-invex set. We have also established an important relation between geodesic $(\alpha,E)$-preinvex function and geodesic $(\alpha,E)$-invex function in a complete Riemannian manifold.


title: Injection locking at fractional frequencies of magnetic tunnel junction (MTJ)-based read sensors' ferromagnetic resonance modes

abstract: Being nonlinear dynamic systems, magnetic read sensors should respond to an excitation signal of a frequency considerably different from their natural ferromagnetic resonance (FMR) frequencies. Because of the magnetization dynamics' inherent nonlinear nature, the sensors' response should be measured at the DC, excitation frequency, and its multiples (harmonics). In this paper, we present results of such measurements, accomplished using a one-port nonlinear vector network analyzer (NVNA), which show distinct resonances at fractional frequencies of the free layer (FL) FMR mode. Identification of these resonances, resulting from the nonlinear nature of the spin-torque (ST)-induced magnetization dynamics, was performed using micromagnetic modeling. In particular, we show that the measured DC response at the above-mentioned fractional frequencies can be explained by a low-order nonlinearity and strong magnetodipolar feedback between magnetic layers adjacent to an MgO barrier. Additionally, we determined that the simulated harmonic response is strongly enhanced by the mutual ST effect between these layers. Finally, we demonstrate that the read sensors' nonlinear magnetization dynamics and, by extension, their harmonic response are highly sensitive to various magnetic and ST parameters. Thus, this study shows that using NVNA measurements in conjunction with micromagnetic modeling can clarify the uncertainty in the definition of these parameters.


title: Entanglement and non-locality between disparate solid-state quantum memories mediated by photons

abstract: Entangling quantum systems with different characteristics through the exchange of photons is a prerequisite for building future quantum networks. Proving the presence of entanglement between quantum memories for light working at different wavelengths furthers this goal. Here, we report on a series of experiments with a thulium-doped crystal, serving as a quantum memory for 794 nm photons, an erbium-doped fibre, serving as a quantum memory for telecommunication-wavelength photons at 1535 nm, and a source of photon pairs created via spontaneous parametric down-conversion. Characterizing the photons after re-emission from the two memories, we find non-classical correlations with a cross-correlation coefficient of $g^{(2)}_{12} = 53\pm8$; entanglement preserving storage with input-output fidelity of $\mathcal{F}_{IO}\approx93\pm2\%$; and non-locality featuring a violation of the Clauser-Horne-Shimony-Holt Bell-inequality with $S= 2.6\pm0.2$. Our proof-of-principle experiment shows that entanglement persists while propagating through different solid-state quantum memories operating at different wavelengths.


title: Localized Surface Plasmon Resonances of Simple Tunable Plasmonic Nanostructures

abstract: We derive and present systematic relationships between the analytical formulas for calculation of the localized surface plasmon resonances (LSPR) of some plasmonic nano-structures which we have categorized as simple. These relationships, including some new formulas, are summarized in a tree diagram which highlights the core-shell plasmons as the generators of solid and cavity plasmons. In addition, we show that the LSPR of complex structures can be reduced to that of simpler ones, using the LSPR of a nanorice as a case study, in the dipole limit. All the formulas were derived using a combination of the Drude model, the Rayleigh approximation, and the Fr\"ohlich condition. The formulas are handy and they are in good agreement with the results of the plasmon hybridization theory. The formulas also account for dielectric effects, which provide versatility in the tuning of the LSPR of the nanostructures. A simplified model of plasmon hybridization is presented, allowing us to investigate the weak-coupling regimes of solid and cavity plasmons in the core-shell nanostructures we have studied.


title: Load Characterization and Power Conditioner Synthesis Using Higher-Order Elements

abstract: It is shown that virtually all nonlinear and/or time-varying loads that generate harmonic current distortion can be characterized in terms of so-called higher-order circuit elements. The most relevant higher-order elements exploited in this paper are the memristor, meminductor, and memcapacitor. Such elements naturally arise by introducing constitutive relationships in terms of higher-order voltage and current differentials and integrals. Consequently, the power conditioner necessary to compensate for the load current distortions is synthesized similarly. The new characterization and compensation synthesis is applied to the half-wave rectifier and the controlled bridge converter.


title: The free field: realization via unbounded operators and Atiyah property

abstract: Let $X_1,\dots,X_n$ be operators in a finite von Neumann algebra and consider their division closure in the affiliated unbounded operators. We address the question when this division closure is a skew field (aka division ring) and when it is the free skew field. We show that the first property is equivalent to the strong Atiyah property and that the second property can be characterized in terms of the non-commutative distribution of $X_1,\dots,X_n$. More precisely, $X_1,\dots,X_n$ generate the free skew field if and only if there exist no non-zero finite rank operators $T_1,\dots,T_n$ such that $\sum_i[T_i,X_i]=0$. Sufficient conditions for this are the maximality of the free entropy dimension or the existence of a dual system of $X_1,\dots,X_n$. Our general theory is not restricted to selfadjoint operators and thus does also include and recover the result of Linnell that the generators of the free group give the free skew field. We give also consequences of our result for the question of atoms in the distribution of rational functions in free variables or in the asymptotic eigenvalue distribution of matrices over polynomials in asymptotically free random matrices. This solves in particular a conjecture of Charlesworth and Shlyakhtenko.


title: A Bayesian Approach to Robust Reinforcement Learning

abstract: Robust Markov Decision Processes (RMDPs) intend to ensure robustness with respect to changing or adversarial system behavior. In this framework, transitions are modeled as arbitrary elements of a known and properly structured uncertainty set and a robust optimal policy can be derived under the worst-case scenario. In this study, we address the issue of learning in RMDPs using a Bayesian approach. We introduce the Uncertainty Robust Bellman Equation (URBE) which encourages safe exploration for adapting the uncertainty set to new observations while preserving robustness. We propose a URBE-based algorithm, DQN-URBE, that scales this method to higher dimensional domains. Our experiments show that the derived URBE-based strategy leads to a better trade-off between less conservative solutions and robustness in the presence of model misspecification. In addition, we show that the DQN-URBE algorithm can adapt significantly faster to changing dynamics online compared to existing robust techniques with fixed uncertainty sets.


title: Third harmonic generation in gapped bilayer graphene

abstract: With the help of numerical simulations in microscopic nonlinear quantum theory of coherent electromagnetic radiation interaction with a gapped bilayer graphene, we find out the optimal values of pump wave intensity, graphene temperature, and energy gap induced by a constant electric field for practically significant third order harmonic coherent emission. The Liouville-von Neumann equation is treated numerically for the third harmonic generation in multiphoton excitation regime near the Dirac points of the Brillouin zone. We examine the rates of the third harmonic at the particle-hole annihilation in the field of a strong pump wave of linear polarization for practically real/optimal parameters of a considering system. The obtained results show that by choosing the optimal values of the main characteristic parameters, a gapped bilayer graphene can serve as an effective medium for generation of the third harmonic at room temperatures in the terahertz and far infrared domains.


title: Electrical conductivity of hot Abelian plasma with scalar charge carriers

abstract: We study the electrical conductivity of hot Abelian plasma containing scalar charge carriers in the leading logarithmic order in coupling constant $\alpha$ using the Boltzmann kinetic equation. The leading contribution to the collision integral is due to the M{\o}ller and Bhabha scattering of scalar particles with a singular cross section in the region of small momentum transfer. Regularizing this singularity by taking into account the hard thermal loop corrections to the propagators of intermediate particles, we derive the second order differential equation which determines the kinetic function. We solve this equation numerically and also use a variational approach in order to find a simple analytical formula for the conductivity. It has the standard parametric dependence on the coupling constant $\sigma\approx 2.38\, T/(\alpha \log\alpha^{-1})$ with the prefactor taking a somewhat lower value compared to the fermionic case. Finally, we consider the general case of hot Abelian plasma with an arbitrary number of scalar and fermionic particle species and derive the simple analytical formula for its conductivity.


title: Measuring Mountains on the Moon

abstract: Following a technique of Galileo we compute the height on a mountain on the Moon. It is based on a simple observation that precisely on a half Moon day, the Earth, the Moon, and the Sun form the vertices of a right triangle with the Moon at the 90 degree vertex.


title: Secure Extensibility for System State Extraction via Plugin Sandboxing

abstract: We introduce a new mechanism to securely extend systems data collection software with potentially untrusted third-party code. Unlike existing tools which run extension modules or plugins directly inside the monitored endpoint (the guest), we run plugins inside a specially crafted sandbox, so as to protect the guest as well as the software core. To get the right mix of accessibility and constraints required for systems data extraction, we create our sandbox by combining multiple features exported by an unmodified kernel. We have tested its applicability by successfully sandboxing plugins of an opensourced data collection software for containerized guest systems. We have also verified its security posture in terms of successful containment of several exploits, which would have otherwise directly impacted a guest, if shipped inside third-party plugins.


title: Exploring helical dynamos with machine learning

abstract: We use ensemble machine learning algorithms to study the evolution of magnetic fields in magnetohydrodynamic (MHD) turbulence that is helically forced. We perform direct numerical simulations of helically forced turbulence using mean field formalism, with electromotive force (EMF) modeled both as a linear and non-linear function of the mean magnetic field and current density. The form of the EMF is determined using regularized linear regression and random forests. We also compare various analytical models to the data using Bayesian inference with Markov Chain Monte Carlo (MCMC) sampling. Our results demonstrate that linear regression is largely successful at predicting the EMF and the use of more sophisticated algorithms (random forests, MCMC) do not lead to significant improvement in the fits. We conclude that the data we are looking at is effectively low dimensional and essentially linear. Finally, to encourage further exploration by the community, we provide all of our simulation data and analysis scripts as open source IPython notebooks.


title: Exploring the damping of Alfv\'en waves along a long off-limb coronal loop, up to 1.4 R$_\odot$

abstract: The Alfv\'en wave energy flux in the corona can be explored using the electron density and velocity amplitude of the waves. The velocity amplitude of Alfv\'en waves can be obtained from the non-thermal velocity of the spectral line profiles. Previous calculations of the Alfv\'en wave energy flux with height in active regions and polar coronal holes have provided evidence for the damping of Alfv\'en waves with height. We present off-limb Hinode EUV imaging spectrometer (EIS) observations of a long coronal loop up to 1.4~R$_\odot$. We have obtained the electron density along the loop and found the loop to be almost in hydrostatic equilibrium. We obtained the temperature using the EM-loci method and found the loop to be isothermal across, as well as along, the loop with a temperature of about 1.37 MK. We significantly improve the estimate of non-thermal velocities over previous studies by using the estimated ion (equal to electron) temperature. Estimates of electron densities are improved using the significant updates of the CHIANTI v.8 atomic data. More accurate measurements of propagating Alfv\'en wave energy along the coronal loop and its damping are presented up to distances of 1.4 R$_\odot$, further than have been previously explored. The Alfv\'en wave energy flux obtained could contribute to a significant part of the coronal losses due to radiation along the loop.


title: The role of zero-clusters in exchange-driven growth with and without input

abstract: The exchange-driven growth model describes the mean field kinetics of a population of composite particles (clusters) subject to pairwise exchange interactions. Exchange in this context means that upon interaction of two clusters, one loses a constituent unit (monomer) and the other gains this unit. Two variants of the exchange-driven growth model appear in applications. They differ in whether clusters of zero size are considered active or passive. In the active case, clusters of size zero can acquire a monomer from clusters of positive size. In the passive case they cannot, meaning that clusters reaching size zero are effectively removed from the system. The large time behaviour is very different for the two variants of the model. We first consider an isolated system. In the passive case, the cluster size distribution tends towards a self-similar evolution and the typical cluster size grows as a power of time. In the active case, we identify a broad class of kernels for which the the cluster size distribution tends to a non-trivial time-independent equilibrium in which the typical cluster size is finite. We next consider a non-isolated system in which monomers are input at a constant rate. In the passive case, the cluster size distribution again attains a self-similar profile in which the typical cluster size grows as a power of time. In the active case, a surprising new behavior is found: the cluster size distribution asymptotes to the same equilibrium profile found in the isolated case but with an amplitude that grows linearly in time.


title: Optimisation of Overparametrized Sum-Product Networks

abstract: It seems to be a pearl of conventional wisdom that parameter learning in deep sum-product networks is surprisingly fast compared to shallow mixture models. This paper examines the effects of overparameterization in sum-product networks on the speed of parameter optimisation. Using theoretical analysis and empirical experiments, we show that deep sum-product networks exhibit an implicit acceleration compared to their shallow counterpart. In fact, gradient-based optimisation in deep tree-structured sum-product networks is equal to gradient ascend with adaptive and time-varying learning rates and additional momentum terms.


title: Coding for Deletion Channels with Multiple Traces

abstract: Motivated by the sequence reconstruction problem from traces in DNA-based storage, we consider the problem of designing codes for the deletion channel when multiple observations (or traces) are available to the decoder. We propose simple binary and non-binary codes based on Varshamov-Tenengolts (VT) codes. The proposed codes split the codeword in blocks and employ a VT code in each block. The availability of multiple traces helps the decoder to identify deletion-free copies of a block, and to avoid mis-synchronization while decoding. The encoding complexity of the proposed scheme is linear in the codeword length; the decoding complexity is linear in the codeword length, and quadratic in the number of deletions and the number of traces. The proposed scheme offers an explicit low-complexity technique for correcting deletions using multiple traces.


title: [OIII] Emission Line Properties in a New Sample of Heavily Reddened Quasars at $z>2$

abstract: We present VLT-SINFONI near infra-red spectra of 26 new heavily reddened quasar candidates selected from the UKIDSS-LAS, VISTA VHS and VIKING imaging surveys. This new sample extends our reddened quasar search to both brighter and fainter luminosities. 25 of the 26 candidates are confirmed to be broad line quasars with redshifts $0.7<z<2.6$ and dust extinctions $0.5<E(B-V)<3.0$. Combining with previously identified samples, we study the H$\alpha$, H$\beta$ and [OIII] emission line properties in 22 heavily reddened quasars with $L_{bol}\approx 10^{47}$ erg/s and $z>2$. We present the first comparison of the [OIII] line properties in high luminosity obscured quasars to a large sample of 111 unobscured quasars in the same luminosity and redshift range. Broad wings extending to velocities of 2500 km/s are seen in the [OIII] emission line profiles of our reddened quasars, suggesting that strong outflows are affecting the ionised gas kinematics. However, we find no significant difference between the kinematics of the [OIII] emission in reddened and unobscured quasars when the two samples are matched in luminosity and redshift. Our results are consistent with a model where quasar-driven outflows persist for some time after the obscuring dust has been cleared from along the line of sight. Assuming the amount of ionised gas in obscured and unobscured quasars is similar, we use the equivalent width distribution of the [OIII] emission to constrain the location of the obscuring dust in our reddened quasars. We find that the dust is most likely to be located on galactic scales, outside the [OIII] emitting region.


title: Increasing the Security of Weak Passwords: the SPARTAN Interface

abstract: Password authentication suffers from the well-known tradeoff between security and usability. Secure passwords are difficult for users to remember, and memorable passwords are often easy to guess. SPARse Two-dimensional AuthenticatioN (SPARTAN) allows users to input their textual passwords in a two-dimensional grid instead of a linear textbox. This interface enables relatively short passwords to have a higher calculated level of security due to the need for an attacker to determine both the text of the password and the location of each character in the grid. We created a SPARTAN prototype and conducted a preliminary user study to evaluate the actual usability and security of the SPARTAN interface compared to the linear password entry interface. We find that while user-created SPARTAN passwords tend to be shorter than their linear counterparts, the calculated security of user-created SPARTAN passwords is higher than that of user-created linear passwords. We also asked participants to complete a survey on the usability of the SPARTAN interface and identified some areas of improvement, while prototype interaction provided evidence of users becoming more familiar with SPARTAN over time. Finally, we performed an investigation into password-cracking tools, and assert that SPARTAN passwords require more resources to crack than their linear counterparts. These findings suggest that SPARTAN is a promising alternative to linear passwords from a security standpoint. Usability of the interface and memorability of SPARTAN passwords is an interesting research question and should be further investigated in future work.


title: Principal component analysis of the primordial tensor power spectrum

abstract: We study how the shape of the spectrum of primordial gravitational waves can be constrained by future experiments looking at the B-mode of the Cosmic Microwave Background (CMB) polarization. We implement a Principal Component Analysis (PCA) including the effects of diffuse foreground residuals, following component separation, in the uncertainty of CMB angular power spectra, and taking into account the gravitational lensing by Large Scale Structure. We perform our study by considering the capabilities of future B-mode CMB experiments such as LiteBIRD, the Simons Observatory (SO) and Stage-IV (CMB-S4), in particular in detecting deviations of the primordial tensor spectrum from the scale-invariant behavior. We find that diffuse foreground residuals impact substantially both the derivation of the PCA basis and the corresponding constraining power, in all cases. In particular, depending on which experimental specifications and which value $r$ of tensor-to-scalar ratio for cosmological perturbations are considered, adding foregrounds residuals can determine an increase as large as a factor $\sim 4$ both on the uncertainty on $r$ and on the recovery of the PCA modes. We study the limitations of the methodology, including the effect of physicality priors on the PCA, which we quantify via a Monte Carlo Markov chain (MCMC) analysis of the combined cosmological and PCA power spectrum parameter space.


title: On the Emergence of Thousands of Absorption Lines in the Quasar PG1411+442: A Clumpy High-Column Density Outflow from the Broad Emission-Line Region?

abstract: Quasar outflows are fundamental components of quasar environments that might play an important role in feedback to galaxy evolution. We report on the emergence of a remarkable new outflow absorption-line system in the quasar PG1411+442 (redshift ~0.089) detected in the UV and visible with the Hubble Space Telescope Cosmic Origins Spectrograph and the Gemini Multi-Object Spectrograph, respectively. This new "transient'' system contains thousands of lines, including FeII and FeII* from excited states up to 3.89 eV, HI* Balmer lines, NaI D 5890,5896, and the first detection of HeI* 5876 in a quasar. The transient absorber is spatially inhomogeneous and compact, with sizes ~<0.003 pc, based on covering fractions on the quasar continuum source ranging from ~0.45 in strong UV lines to ~0.04 in NaI D. Cloudy photoionization simulations show that large total column densities log N_H(cm^-2) >~ 23.4 and an intense radiation field ~<0.4~pc from the quasar are needed to produce the observed lines in thick zones of both fully-ionised and partially-ionised gas. The densities are conservatively log n_H(cm-3) >~ 7 based on FeII*, HI*, and HeI* but they might reach log n_H(cm^-3) >~ 10 based on NaI D. The transient lines appear at roughly the same velocity shift, v ~ -1900 km/s, as a "mini-BAL'' outflow detected previously, but with narrower Doppler widths, b ~ 100 km/s, and larger column densities in more compact outflow structures. We propose that the transient lines identify a clumpy outflow from the broad emission-line region that, at its current speed and location, is still gravitationally bound to the central black hole.


title: Realizing realizability results with classical constructions

abstract: J.L. Krivine developed a new method based on realizability to construct models of set theory where the axiom of choice fails. We attempt to recreate his results in classical settings, i.e. symmetric extensions. We also provide a new condition for preserving well-ordered, and other particular type of choice, in the general settings of symmetric extensions.


title: On the sharp stability of critical points of the Sobolev inequality

abstract: Given $n\geq 3$, consider the critical elliptic equation $\Delta u + u^{2^*-1}=0$ in $\mathbb R^n$ with $u > 0$. This equation corresponds to the Euler-Lagrange equation induced by the Sobolev embedding $H^1(\mathbb R^n)\hookrightarrow L^{2^*}(\mathbb R^n)$, and it is well-known that the solutions are uniquely characterized and are given by the so-called ``Talenti bubbles''. In addition, thanks to a fundamental result by Struwe, this statement is ``stable up to bubbling'': if $u:\mathbb R^n\to(0,\infty)$ almost solves $\Delta u + u^{2^*-1}=0$ then $u$ is (nonquantitatively) close in the $H^1(\mathbb R^n)$-norm to a sum of weakly-interacting Talenti bubbles. More precisely, if $\delta(u)$ denotes the $H^1(\mathbb R^n)$-distance of $u$ from the manifold of sums of Talenti bubbles, Struwe proved that $\delta(u)\to 0$ as $\lVert\Delta u + u^{2^*-1}\rVert_{H^{-1}}\to 0$. In this paper we investigate the validity of a sharp quantitative version of the stability for critical points: more precisely, we ask whether under a bound on the energy $\lVert\nabla u\rVert_{L^2}$ (that controls the number of bubbles) it holds $\delta(u) \lesssim \lVert\Delta u + u^{2^*-1}\rVert_{H^{-1}}$. A recent paper by the first author together with Ciraolo and Maggi shows that the above result is true if $u$ is close to only one bubble. Here we prove, to our surprise, that whenever there are at least two bubbles then the estimate above is true for $3\le n\le 5$ while it is false for $n\ge 6$. To our knowledge, this is the first situation where quantitative stability estimates depend so strikingly on the dimension of the space, changing completely behavior for some particular value of the dimension $n$.


title: Custom Execution Environments with Containers in Pegasus-enabled Scientific Workflows

abstract: Science reproducibility is a cornerstone feature in scientific workflows. In most cases, this has been implemented as a way to exactly reproduce the computational steps taken to reach the final results. While these steps are often completely described, including the input parameters, datasets, and codes, the environment in which these steps are executed is only described at a higher level with endpoints and operating system name and versions. Though this may be sufficient for reproducibility in the short term, systems evolve and are replaced over time, breaking the underlying workflow reproducibility. A natural solution to this problem is containers, as they are well defined, have a lifetime independent of the underlying system, and can be user-controlled so that they can provide custom environments if needed. This paper highlights some unique challenges that may arise when using containers in distributed scientific workflows. Further, this paper explores how the Pegasus Workflow Management System implements container support to address such challenges.