title: Cosmic String Interpretation of NANOGrav Pulsar Timing Data

abstract: Pulsar timing data used to provide upper limits on a possible stochastic gravitational wave background (SGWB). However, the NANOGrav Collaboration has recently reported strong evidence for a stochastic common-spectrum process, which we interpret as a SGWB in the framework of cosmic strings. The possible NANOGrav signal would correspond to a string tension $G\mu \in (4 \times 10^{-11}, 10^{-10}) $ at the 68% confidence level, with a different frequency dependence from supermassive black hole mergers. The SGWB produced by cosmic strings with such values of $G\mu$ would be beyond the reach of LIGO, but could be measured by other planned and proposed detectors such as SKA, LISA, TianQin, AION-1km, AEDGE, Einstein Telescope and Cosmic Explorer.

categories: astro-ph.CO astro-ph.HE gr-qc hep-ph hep-th


title: Statistical mechanics of an integrable system

abstract: We provide here an explicit example of Khinchin's idea that the validity of equilibrium statistical mechanics in high dimensional systems does not depend on the details of the dynamics. This point of view is supported by extensive numerical simulation of the one-dimensional Toda chain, an integrable non-linear Hamiltonian system where all Lyapunov exponents are zero by definition. We study the relaxation to equilibrium starting from very atypical initial conditions and focusing on energy equipartion among Fourier modes, as done in the original Fermi-Pasta-Ulam-Tsingou numerical experiment. We find evidence that in the general case, i.e., not in the perturbative regime where Toda and Fourier modes are close to each other, there is a fast reaching of thermal equilibrium in terms of a single temperature. We also find that equilibrium fluctuations, in particular the behaviour of specific heat as function of temperature, are in agreement with analytic predictions drawn from the ordinary Gibbs ensemble, still having no conflict with the established validity of the Generalized Gibbs Ensemble for the Toda model. Our results suggest thus that even an integrable Hamiltonian system reaches thermalization on the constant energy hypersurface, provided that the considered observables do not strongly depend on one or few of the conserved quantities. This suggests that dynamical chaos is irrelevant for thermalization in the large-$N$ limit, where any macroscopic observable reads of as a collective variable with respect to the coordinate which diagonalize the Hamiltonian. The possibility for our results to be relevant for the problem of thermalization in generic quantum systems, i.e., non-integrable ones, is commented at the end.


title: Effective Federated Adaptive Gradient Methods with Non-IID Decentralized Data

abstract: Federated learning allows loads of edge computing devices to collaboratively learn a global model without data sharing. The analysis with partial device participation under non-IID and unbalanced data reflects more reality. In this work, we propose federated learning versions of adaptive gradient methods - Federated AGMs - which employ both the first-order and second-order momenta, to alleviate generalization performance deterioration caused by dissimilarity of data population among devices. To further improve the test performance, we compare several schemes of calibration for the adaptive learning rate, including the standard Adam calibrated by $\epsilon$, $p$-Adam, and one calibrated by an activation function. Our analysis provides the first set of theoretical results that the proposed (calibrated) Federated AGMs converge to a first-order stationary point under non-IID and unbalanced data settings for nonconvex optimization. We perform extensive experiments to compare these federated learning methods with the state-of-the-art FedAvg, FedMomentum and SCAFFOLD and to assess the different calibration schemes and the advantages of AGMs over the current federated learning methods.


title: Vector copulas

abstract: This paper introduces vector copulas associated with multivariate distributions with given multivariate marginals, based on the theory of measure transportation, and establishes a vector version of Sklar's theorem. The latter provides a theoretical justification for the use of vector copulas to characterize nonlinear or rank dependence between a finite number of random vectors (robust to within vector dependence), and to construct multivariate distributions with any given non overlapping multivariate marginals. We construct Elliptical and Kendall families of vector copulas, derive their densities, and present algorithms to generate data from them. The use of vector copulas is illustrated with a stylized analysis of international financial contagion.


title: A probabilistic approach to exhaustion in the infinite-genus case

abstract: We explore the use of Costa and Farber's model for random simplicial complexes to give probabilistic evidence for exhaustion via rigid expansions on random simplicial complexes which are analogous of curve complexes. This has potential applications to action rigidity following Ivanov's meta-conjecture.


title: Dual-Mandate Patrols: Multi-Armed Bandits for Green Security

abstract: Conservation efforts in green security domains to protect wildlife and forests are constrained by the limited availability of defenders (i.e., patrollers), who must patrol vast areas to protect from attackers (e.g., poachers or illegal loggers). Defenders must choose how much time to spend in each region of the protected area, balancing exploration of infrequently visited regions and exploitation of known hotspots. We formulate the problem as a stochastic multi-armed bandit, where each action represents a patrol strategy, enabling us to guarantee the rate of convergence of the patrolling policy. However, a naive bandit approach would compromise short-term performance for long-term optimality, resulting in animals poached and forests destroyed. To speed up performance, we leverage smoothness in the reward function and decomposability of actions. We show a synergy between Lipschitz-continuity and decomposition as each aids the convergence of the other. In doing so, we bridge the gap between combinatorial and Lipschitz bandits, presenting a no-regret approach that tightens existing guarantees while optimizing for short-term performance. We demonstrate that our algorithm, LIZARD, improves performance on real-world poaching data from Cambodia.


title: Critical field anisotropy in the antiferroelectric switching of PbZrO3 films

abstract: Antiferroelectrics have been recently sparking interest due to their potential use in energy storage and electrocaloric cooling. Their main distinctive feature is antiferroelectric switching, i.e. the possibility to induce a phase transition to a polar phase by an electric field. Here we investigate the switching behavior of the model antiferroelectric perovskite PbZrO3 using thin films processed by chemical solution deposition in different geometries and orientations. Both out-of-plane and in-plane switching configurations are investigated. The critical field is observed to be highly dependent on the direction of the electric field with respect to the film texture. We show that this behaviour is qualitatively consistent with a phase transition to a rhombohedral polar phase. We finally estimate the importance of crystallite orientation and film texturation in the variations observed in the literature.


title: Effective Proximal Methods for Non-convex Non-smooth Regularized Learning

abstract: Sparse learning is a very important tool for mining useful information and patterns from high dimensional data. Non-convex non-smooth regularized learning problems play essential roles in sparse learning, and have drawn extensive attentions recently. We design a family of stochastic proximal gradient methods by applying arbitrary sampling to solve the empirical risk minimization problem with a non-convex and non-smooth regularizer. These methods draw mini-batches of training examples according to an arbitrary probability distribution when computing stochastic gradients. A unified analytic approach is developed to examine the convergence and computational complexity of these methods, allowing us to compare the different sampling schemes. We show that the independent sampling scheme tends to improve performance over the commonly-used uniform sampling scheme. Our new analysis also derives a tighter bound on convergence speed for the uniform sampling than the best one available so far. Empirical evaluations demonstrate that the proposed algorithms converge faster than the state of the art.


title: Generating Functions for Some Families of the Generalized Al-Salam-Carlitz $q$-Polynomials

abstract: In this paper, by making use of the familiar $q$-difference operators $D_q$ and $D_{q^{-1}}$, we first introduce two homogeneous $q$-difference operators $\mathbb{T}({\bf a},{\bf b},cD_q)$ and $\mathbb{E}({\bf a},{\bf b}, cD_{q^{-1}})$, which turn out to be suitable for dealing with the families of the generalized Al-Salam-Carlitz $q$-polynomials $\phi_n^{({\bf a},{\bf b})}(x,y|q)$ and $\psi_n^{({\bf a},{\bf b})}(x,y|q)$. We then apply each of these two homogeneous $q$-difference operators in order to derive generating functions, Rogers type formulas, the extended Rogers type formulas and the Srivastava-Agarwal type linear as well as bilinear generating functions involving each of these families of the generalized Al-Salam-Carlitz $q$-polynomials. We also show how the various results presented here are related to those in many earlier works on the topics which we study in this paper.


title: Defects in the 3-dimensional toric code model form a braided fusion 2-category

abstract: It was well known that there are $e$-particles and $m$-strings in the 3-dimensional (spatial dimension) toric code model, which realizes the 3-dimensional $\mathbb{Z}_2$ topological order. Recent mathematical result, however, shows that there are additional string-like topological defects in the 3-dimensional $\mathbb{Z}_2$ topological order. In this work, we construct all topological defects of codimension 2 and higher, and show that they form a braided fusion 2-category satisfying a braiding non-degeneracy condition.


title: Classifying Intrinsically Linked Tournaments by Score Sequence

abstract: A tournament on 8 or more vertices may be intrinsically linked as a directed graph. We begin the classification of intrinsically linked tournaments by examining their score sequences. While many distinct tournaments may have the same score sequence, there exist score sequences $S$ such that any tournament with score sequence $S$ has an embedding with no nonsplit consistently oriented link. We call such score sequences $\textit{linkless}$, and we show that the vast majority of score sequences for 8 vertex tournaments are linkless. We also extend these results to $n$ vertex tournaments and are able to classify many longer score sequences as well. We show that for any $n$, there exist at least $O(n)$ linkless score sequences, but we conjecture that the fraction of score sequences of length $n$ that are linkless goes to 0 as $n$ becomes large.


title: Singular holomorphic foliations by curves. III: Zero Lelong numbers

abstract: Let $\mathcal{F}$ be a holomorphic foliation by curves defined in a neighborhood of $0$ in $\mathbb{C}^n$ ($n\geq 2$) having $0$ as a weakly hyperbolic singularity. Let $T$ be a positive harmonic current directed by $\mathcal{F}$ which does not give mass to any of the $n$ coordinate invariant hyperplanes $\{z_j=0\}$ for $1\leq j\leq n.$ Then we show that the Lelong number of $T$ at $0$ vanishes. Moreover, an application of this local result in the global context is given. We discuss also the relation between several basic notions such as directed positive harmonic currents, directed positive ddc-closed currents, Lelong numbers etc. in the framework of singular holomorphic foliations.


title: Conformal Defects from String Field Theory

abstract: Unlike conformal boundary conditions, conformal defects of Virasoro minimal models lack classification. Alternatively to the defect perturbation theory and the truncated conformal space approach, we employ open string field theory (OSFT) techniques to explore the space of conformal defects. We illustrate the method by an analysis of OSFT around the background associated to the $(1,2)$ topological defect in diagonal unitary minimal models. Numerical analysis of OSFT equations of motion leads to an identification of a nice family of solutions, recovering the picture of infrared fixed points due to Kormos, Runkel and Watts. In particular, we find a continuum of solutions in the Ising model case and 6 solutions for other minimal models. OSFT provides us with numerical estimates of the g-function and other coefficients of the boundary state.


title: The space of properly-convex structures

abstract: Suppose $G$ is finitely generated group and $\mathcal{C}(G)$ consists of all $\rho:G\to\operatorname{PGL}(n+1,\mathbb{R})$ for which there exists a properly convex set in $\mathbb{R}\mathbb{P}^n$ that is preserved by $\rho(G)$. Then the image of $\mathcal{C}(G)$ is closed in the character variety. Suppose $G$ does not contain an infinite, normal, abelian subgroup and $\mathcal{D}(G)\subset\mathcal{C}(G)$ is the subset of holonomies of properly-convex $n$-manifolds with fundamental group $G$. Then the image $\mathcal{D}(G)$ is closed in the character variety. If $M$ is the interior of a compact $n$-manifold and $G=\pi_1M$ is as above, and either $M$ is closed, or $\pi_1M$ contains a subgroup of infinite index isomorphic to $\mathbb{Z}^{n-1}$, then $\mathcal{D}(G)$ is closed. If, in addition, $M$ is the interior of a compact manifold $N$ such that every component of $\partial N$ is $\pi_1$-injective, and finitely covered by a torus, then every element of $\mathcal{D}(G)$ is the holonomy of a properly-convex structure on $M$, and $\mathcal{D}(G)$ is a union of connected components of a semi-algebraic set.


title: On properly convex real-projective manifolds with Generalized Cusp

abstract: Suppose $E$ is an end of an irreducible, properly convex, real-projective $n$-manifold $M$. If $\pi_1E$ contains a subgroup of finite index isomorphic to ${\mathbb Z}^{n-1}$, and $E\hookrightarrow M$ is $\pi_1$-injective, then $E$ is a generalized cusp. We list some consequences when all ends are of this type. Under certain hypotheses we prove the holonomy of a properly convex manifold is irreducible.


title: Spatial Differencing for Sample Selection Models with Unobserved Heterogeneity

abstract: This paper derives identification, estimation, and inference results using spatial differencing in sample selection models with unobserved heterogeneity. We show that under the assumption of smooth changes across space of the unobserved sub-location specific heterogeneities and inverse Mills ratio, key parameters of a sample selection model are identified. The smoothness of the sub-location specific heterogeneities implies a correlation in the outcomes. We assume that the correlation is restricted within a location or cluster and derive asymptotic results showing that as the number of independent clusters increases, the estimators are consistent and asymptotically normal. We also propose a formula for standard error estimation. A Monte-Carlo experiment illustrates the small sample properties of our estimator. The application of our procedure to estimate the determinants of the municipality tax rate in Finland shows the importance of accounting for unobserved heterogeneity.


title: Input Hessian Regularization of Neural Networks

abstract: Regularizing the input gradient has shown to be effective in promoting the robustness of neural networks. The regularization of the input's Hessian is therefore a natural next step. A key challenge here is the computational complexity. Computing the Hessian of inputs is computationally infeasible. In this paper we propose an efficient algorithm to train deep neural networks with Hessian operator-norm regularization. We analyze the approach theoretically and prove that the Hessian operator norm relates to the ability of a neural network to withstand an adversarial attack. We give a preliminary experimental evaluation on the MNIST and FMNIST datasets, which demonstrates that the new regularizer can, indeed, be feasible and, furthermore, that it increases the robustness of neural networks over input gradient regularization.


title: The optimal spectral gap for regular and disordered harmonic networks of oscillators

abstract: We consider one-dimensional chains and multi-dimensional networks of harmonic oscillators coupled to two Langevin heat reservoirs at different temperatures. Each particle interacts with its nearest neighbors by harmonic potentials and all individual particles are confined by harmonic potentials, too. In this article, we provide, for the first time, the sharp N dependence of the spectral gap of the associated generator under various physical assumptions and for different spatial dimensions. Our method of proof relies on a new approach to analyze a non self-adjoint eigenvalue problem involving low-rank non-hermitian perturbations of auxiliary discrete Schr\"odinger operators.


title: Themes Informed Audio-visual Correspondence Learning

abstract: The applications of short-term user-generated video (UGV), such as Snapchat, and Youtube short-term videos, booms recently, raising lots of multimodal machine learning tasks. Among them, learning the correspondence between audio and visual information from videos is a challenging one. Most previous work of the audio-visual correspondence(AVC) learning only investigated constrained videos or simple settings, which may not fit the application of UGV. In this paper, we proposed new principles for AVC and introduced a new framework to set sight of videos' themes to facilitate AVC learning. We also released the KWAI-AD-AudVis corpus which contained 85432 short advertisement videos (around 913 hours) made by users. We evaluated our proposed approach on this corpus, and it was able to outperform the baseline by 23.15% absolute difference.


title: Interactive Focus+Context Rendering for Hexahedral Mesh Inspection

abstract: The visual inspection of a hexahedral mesh with respect to element quality is difficult due to clutter and occlusions that are produced when rendering all element faces or their edges simultaneously. Current approaches overcome this problem by using focus on specific elements that are then rendered opaque, and carving away all elements occluding their view. In this work, we make use of advanced GPU shader functionality to generate a focus+context rendering that highlights the elements in a selected region and simultaneously conveys the global mesh structure in the surrounding. To achieve this, we propose a gradual transition from edge-based focus rendering to volumetric context rendering, by combining fragment shader-based edge and face rendering with per-pixel fragment lists. A fragment shader smoothly transitions between wireframe and face-based rendering, including focus-dependent rendering style and depth-dependent edge thickness and halos, and per-pixel fragment lists are used to blend fragments in correct visibility order. To maintain the global mesh structure in the context regions, we propose a new method to construct a sheet-based level-of-detail hierarchy and smoothly blend it with volumetric information. The user guides the exploration process by moving a lens-like hotspot. Since all operations are performed on the GPU, interactive frame rates are achieved even for large meshes.


title: Unobstructed deformation problems for GSp(4)

abstract: We prove, for many cuspidal automorphic representations for GSp(4), that the local obstructions to the deformation theory of the associated residual Galois representations generically vanish.


title: Disease control as an optimization problem

abstract: In the context of epidemiology, policies for disease control are often devised through a mixture of intuition and brute-force, whereby the set of logically conceivable policies is narrowed down to a small family described by a few parameters, following which linearization or grid search is used to identify the optimal policy within the set. This scheme runs the risk of leaving out more complex (and perhaps counter-intuitive) policies for disease control that could tackle the disease more efficiently. In this article, we use techniques from convex optimization theory and machine learning to conduct optimizations over disease policies described by hundreds of parameters. In contrast to past approaches for policy optimization based on control theory, our framework can deal with arbitrary uncertainties on the initial conditions and model parameters controlling the spread of the disease, and stochastic models. In addition, our methods allow for optimization over policies which remain constant over weekly periods, specified by either continuous or discrete (e.g.: lockdown on/off) government measures. We illustrate our approach by minimizing the total time required to eradicate COVID-19 within the Susceptible-Exposed-Infected-Recovered (SEIR) model proposed by Kissler \emph{et al.} (March, 2020).


title: Push-forward method for piecewise deterministic biochemical simulations

abstract: A biochemical network can be simulated by a set of ordinary differential equations (ODE) under well stirred reactor conditions, for large numbers of molecules, and frequent reactions. This is no longer a robust representation when some molecular species are in small numbers and reactions changing them are infrequent. In this case, discrete stochastic events trigger changes of the smooth deterministic dynamics of the biochemical network. Piecewise-deterministic Markov processes (PDMP) are well adapted for describing such situations. Although PDMP models are now well established in biology, these models remain computationally challenging. Previously we have introduced the push-forward method to compute how the probability measure is spread by the deterministic ODE flow of PDMPs, through the use of analytic expressions of the corresponding semigroup. In this paper we provide a more general simulation algorithm that works also for non-integrable systems. The method can be used for biochemical simulations with applications in fundamental biology, biotechnology and biocomputing.This work is an extended version of the work presented at the conference CMSB2019.


title: Braiding and all quantum operations with Majorana modes in 1D

abstract: We propose a scheme to perform braiding and all other unitary operations with Majorana modes in 1D that, in contrast to previous proposals, is solely based on resonant manipulation involving the first excited state extended over the modes. The detection of the population of the excited state also enables initialization and read-out. We provide an elaborated illustration of the scheme with a concrete device.


title: Reinforcement Learning for Dynamic Resource Optimization in 5G Radio Access Network Slicing

abstract: The paper presents a reinforcement learning solution to dynamic resource allocation for 5G radio access network slicing. Available communication resources (frequency-time blocks and transmit powers) and computational resources (processor usage) are allocated to stochastic arrivals of network slice requests. Each request arrives with priority (weight), throughput, computational resource, and latency (deadline) requirements, and if feasible, it is served with available communication and computational resources allocated over its requested duration. As each decision of resource allocation makes some of the resources temporarily unavailable for future, the myopic solution that can optimize only the current resource allocation becomes ineffective for network slicing. Therefore, a Q-learning solution is presented to maximize the network utility in terms of the total weight of granted network slicing requests over a time horizon subject to communication and computational constraints. Results show that reinforcement learning provides major improvements in the 5G network utility relative to myopic, random, and first come first served solutions. While reinforcement learning sustains scalable performance as the number of served users increases, it can also be effectively used to assign resources to network slices when 5G needs to share the spectrum with incumbent users that may dynamically occupy some of the frequency-time blocks.


title: Near-horizon aspects of acceleration radiation by free fall of an atom into a black hole

abstract: A two-level atom freely falling towards a Schwarzschild black hole was recently shown to detect radiation in the Boulware vacuum in an insightful paper [M. O. Scully et al., Proc. Natl. Acad. Sci. U.S.A. 115, 8131 (2018)]. The two-state atom acts as a dipole detector and its interaction with the field can be modeled using a quantum optics approach. The relative acceleration between the scalar field and the detector causes the atom to detect the radiation. In this paper, we show that this acceleration radiation is driven by the near-horizon physics of the black hole. This insight reinforces the relevance of near-horizon conformal quantum mechanics for all the physics associated with the thermodynamic properties of the black hole. We additionally highlight the conformal aspects of the radiation that is given by a Planck distribution with the Hawking temperature.


title: Twisted $L^2$-torsion on the character variety

abstract: We define a twisted $L^2$-torsion on the character variety of 3-manifold $M$ and study some of its properties. In the case where $M$ is hyperbolic of finite volume, we prove that the $L^2$-torsion is a real analytic function on a neighborhood of any lift of the holonomy representation.


title: The Space of Strictly-convex Real-projective structures on a closed manifold

abstract: This is an expository proof that, if $M$ is a compact $n$-manifold with no boundary, then the set of holonomies of strictly-convex real-projective structures on $M$ is a subset of $\operatorname{Hom}(\pi_1M,\operatorname{PGL}(n+1,\mathbb RR))$ that is both open and closed.


title: Rational values of powers of trigonometric functions

abstract: We extend the theorem by Olmsted (1945) and Carlitz-Thomas (1963) on rational values of trigonometric functions to powers of trigonometric functions.


title: Antiferromagnetic Dispersion Relations and Nature of Magnon Pressure

abstract: We derive higher-order corrections in the magnon dispersion relations for two- and three-dimensional antiferromagnets exposed to magnetic and staggered fields that are mutually aligned. "Dressing" the magnons is the prerequisite to separate the low-temperature representation of the pressure into a piece due to noninteracting magnons and a piece that corresponds to the magnon-magnon interaction. Both in two and three spatial dimensions, the interaction in the pressure turns out to be attractive. While concrete figures refer to the spin-$\frac{1}{2}$ square-lattice and the spin-$\frac{1}{2}$ simple cubic lattice antiferromagnet, our results are valid for arbitrary bipartite geometry.


title: Is GW190521 the merger of black holes from the first stellar generations?

abstract: GW190521 challenges our understanding of the late-stage evolution of massive stars and the effects of the pair-instability in particular. We discuss the possibility that stars at low or zero metallicity could retain most of their hydrogen envelope until the pre-supernova stage, avoid the pulsational pair-instability regime and produce a black hole with a mass in the mass gap by fallback. We present a series of new stellar evolution models at zero and low metallicity computed with the Geneva and MESA stellar evolution codes and compare to existing grids of models. Models with a metallicity in the range 0-0.0004 have three properties which favour higher BH masses as compared to higher metallicity models. These are (i) lower mass-loss rates during the post-MS phase, (ii) a more compact star disfavouring binary interaction and (iii) possible H-He shell interactions which lower the CO core mass. We conclude that it is possible that GW190521 may be the merger of black holes produced directly by massive stars from the first stellar generations. Our models indicate BH masses up to 70-75 Msun. Uncertainties related to convective mixing, mass loss, H-He shell interactions and pair-instability pulsations may increase this limit to ~85 Msun.


title: Zero-shot Synthesis with Group-Supervised Learning

abstract: Visual cognition of primates is superior to that of artificial neural networks in its ability to 'envision' a visual object, even a newly-introduced one, in different attributes including pose, position, color, texture, etc. To aid neural networks to envision objects with different attributes, we propose a family of objective functions, expressed on groups of examples, as a novel learning framework that we term Group-Supervised Learning (GSL). GSL allows us to decompose inputs into a disentangled representation with swappable components, that can be recombined to synthesize new samples. For instance, images of red boats & blue cars can be decomposed and recombined to synthesize novel images of red cars. We propose an implementation based on auto-encoder, termed group-supervised zero-shot synthesis network (GZS-Net) trained with our learning framework, that can produce a high-quality red car even if no such example is witnessed during training. We test our model and learning framework on existing benchmarks, in addition to anew dataset that we open-source. We qualitatively and quantitatively demonstrate that GZS-Net trained with GSL outperforms state-of-the-art methods.


title: Fast high-fidelity multi-qubit state transfer with long-range interactions

abstract: We describe an efficient protocol to perform quantum state transfer using Hamiltonian dynamics with long-range interactions. The time to transfer $n$ qubits a sufficiently large distance is proportional to $\sqrt{n}$. Even without error correction, the fidelity of this multi-qubit state transfer process remains finite for arbitrarily well-separated qubits in the presence of uncorrelated random errors in coupling constants.


title: Crack arrest through branching at curved weak interfaces: an experimental and numerical study

abstract: The phenomenon of arrest of an unstably-growing crack due to a curved weak interface is investigated. The weak interface can produce the deviation of the crack path, trapping the crack at the interface, leading to stable crack growth for certain interface geometries. This idea could be used as a technical solution for a new type of crack arrester, with a negligible impact on the global stiffness, strength and weight of the structure. In order to exploit this concept, an experimental campaign based on photo-elasticity and digital image correlation is carried out, showing the capability of curved weak interfaces to arrest cracks. The experiment is repeated for several geometrical configurations through the modification of the interface curvature radii. The phenomenon of crack deviation and subsequent arrest at the interface is also investigated with the assistance of a computational model based on the finite element method. The computational predictions provide the rationale for the interpretation of the experimental observations, and distinguish between the different behaviour of concave and convex interfaces. Consequently, as is shown in the present study, the curved interface concept fosters new routes for the attainment of structures with enhanced fracture resistance capacities, which are of paramount importance for materials and components used in extreme conditions.


title: Robust Deep Learning Ensemble against Deception

abstract: Deep neural network (DNN) models are known to be vulnerable to maliciously crafted adversarial examples and to out-of-distribution inputs drawn sufficiently far away from the training data. How to protect a machine learning model against deception of both types of destructive inputs remains an open challenge. This paper presents XEnsemble, a diversity ensemble verification methodology for enhancing the adversarial robustness of DNN models against deception caused by either adversarial examples or out-of-distribution inputs. XEnsemble by design has three unique capabilities. First, XEnsemble builds diverse input denoising verifiers by leveraging different data cleaning techniques. Second, XEnsemble develops a disagreement-diversity ensemble learning methodology for guarding the output of the prediction model against deception. Third, XEnsemble provides a suite of algorithms to combine input verification and output verification to protect the DNN prediction models from both adversarial examples and out of distribution inputs. Evaluated using eleven popular adversarial attacks and two representative out-of-distribution datasets, we show that XEnsemble achieves a high defense success rate against adversarial examples and a high detection success rate against out-of-distribution data inputs, and outperforms existing representative defense methods with respect to robustness and defensibility.


title: Gaussian phase sensitivity of boson-sampling-inspired strategies

abstract: In this work we study the phase sensitivity of generic linear interferometric schemes using Gaussian resources and measurements. Our formalism is based on the Fisher information. This allows us to separate the contributions of the measurement scheme, the experimental imperfections, and auxiliary systems. We demonstrate the strength of this formalism using a broad class of multimode Gaussian states that includes well-known results from single- and two-mode metrology scenarios. Using this, we prove that input coherent states or squeezing beat the non-classical states proposed in preceding boson-sampling-inspired phase-estimation schemes. We also develop a novel polychromatic interferometric protocol, demonstrating an enhanced sensitivity with respect to two-mode squeezed-vacuum states, for which the ideal homodyne detection is formally shown to be optimal.


title: Understanding Gesture and Speech Multimodal Interactions for Manipulation Tasks in Augmented Reality Using Unconstrained Elicitation

abstract: This research establishes a better understanding of the syntax choices in speech interactions and of how speech, gesture, and multimodal gesture and speech interactions are produced by users in unconstrained object manipulation environments using augmented reality. The work presents a multimodal elicitation study conducted with 24 participants. The canonical referents for translation, rotation, and scale were used along with some abstract referents (create, destroy, and select). In this study time windows for gesture and speech multimodal interactions are developed using the start and stop times of gestures and speech as well as the stoke times for gestures. While gestures commonly precede speech by 81 ms we find that the stroke of the gesture is commonly within 10 ms of the start of speech. Indicating that the information content of a gesture and its co-occurring speech are well aligned to each other. Lastly, the trends across the most common proposals for each modality are examined. Showing that the disagreement between proposals is often caused by a variation of hand posture or syntax. Allowing us to present aliasing recommendations to increase the percentage of users' natural interactions captured by future multimodal interactive systems.


title: Analogy-Making as a Core Primitive in the Software Engineering Toolbox

abstract: An analogy is an identification of structural similarities and correspondences between two objects. Computational models of analogy making have been studied extensively in the field of cognitive science to better understand high-level human cognition. For instance, Melanie Mitchell and Douglas Hofstadter sought to better understand high-level perception by developing the Copycat algorithm for completing analogies between letter sequences. In this paper, we argue that analogy making should be seen as a core primitive in software engineering. We motivate this argument by showing how complex software engineering problems such as program understanding and source-code transformation learning can be reduced to an instance of the analogy-making problem. We demonstrate this idea using Sifter, a new analogy-making algorithm suitable for software engineering applications that adapts and extends ideas from Copycat. In particular, Sifter reduces analogy-making to searching for a sequence of update rule applications. Sifter uses a novel representation for mathematical structures capable of effectively representing the wide variety of information embedded in software. We conclude by listing major areas of future work for Sifter and analogy-making in software engineering.


title: Phosphine Gas in the Cloud Decks of Venus

abstract: Measurements of trace-gases in planetary atmospheres help us explore chemical conditions different to those on Earth. Our nearest neighbor, Venus, has cloud decks that are temperate but hyper-acidic. We report the apparent presence of phosphine (PH3) gas in Venusian atmosphere, where any phosphorus should be in oxidized forms. Single-line millimeter-waveband spectral detections (quality up to ~15 sigma) from the JCMT and ALMA telescopes have no other plausible identification. Atmospheric PH3 at ~20 parts-per-billion abundance is inferred. The presence of phosphine is unexplained after exhaustive study of steady-state chemistry and photochemical pathways, with no currently-known abiotic production routes in Venusian atmosphere, clouds, surface and subsurface, or from lightning, volcanic or meteoritic delivery. Phosphine could originate from unknown photochemistry or geochemistry, or, by analogy with biological production of phosphine on Earth, from the presence of life. Other PH3 spectral features should be sought, while in-situ cloud and surface sampling could examine sources of this gas.


title: Memristive control of mutual SHNO synchronization for neuromorphic computing

abstract: Synchronization of large spin Hall nano-oscillators (SHNO) arrays is an appealing approach toward ultra-fast non-conventional computing based on nanoscale coupled oscillator networks. However, for large arrays, interfacing to the network, tuning its individual oscillators, their coupling, and providing built-in memory units for training purposes, remain substantial challenges. Here, we address all these challenges using memristive gating of W/CoFeB/MgO/AlOx based SHNOs. In its high resistance state (HRS), the memristor modulates the perpendicular magnetic anisotropy (PMA) at the CoFeB/MgO interface purely by the applied electric field. In its low resistance state (LRS), and depending on the voltage polarity, the memristor adds/subtracts current to/from the SHNO drive. The operation in both the HRS and LRS affects the SHNO auto-oscillation mode and frequency, which can be tuned up to 28 MHz/V. This tuning allows us to reversibly turn on/off mutual synchronization in chains of four SHNOs. We also demonstrate two individually controlled memristors to tailor both the coupling strength and the frequency of the synchronized state. Memristor gating is therefore an efficient approach to input, tune, and store the state of the SHNO array for any non-conventional computing paradigm, all in one platform.


title: Geometric properties of the Kazhdan-Lusztig Schubert basis

abstract: We study classes determined by the Kazhdan-Lusztig basis of the Hecke algebra in the $K$-theory and hyperbolic cohomology theory of flag varieties. We first show that, in $K$-theory, the two different choices of Kazhdan-Lusztig bases produce dual bases, one of which can be interpreted as characteristic classes of the intersection homology mixed Hodge modules. In equivariant hyperbolic cohomology, we show that if the Schubert variety is smooth, then the class it determines coincides with the class of the Kazhdan-Lusztig basis; this was known as the Smoothness Conjecture. For Grassmannians, we prove that the classes of the Kazhdan-Lusztig basis coincide with the classes determined by Zelevinsky's small resolutions. These properties of the so-called KL-Schubert basis show that it is the closest existing analogue to the Schubert basis for hyperbolic cohomology; the latter is a very useful testbed for more general elliptic cohomologies.


title: Sharp nonuniqueness for the Navier-Stokes equations

abstract: In this paper, we prove a sharp nonuniqueness result for the incompressible Navier-Stokes equations in the periodic setting. In any dimension $d \geq 2$ and given any $ p<2$, we show the nonuniqueness of weak solutions in the class $L^{p}_t L^\infty$, which is sharp in view of the classical Ladyzhenskaya-Prodi-Serrin criteria. The proof is based on the construction of a class of non-Leray-Hopf weak solutions. More specifically, for any $ p<2$, $q<\infty$, and $\varepsilon>0$, we construct non-Leray-Hopf weak solutions $ u \in L^{p}_t L^\infty \cap L^1_t W^{1,q}$ that are smooth outside a set of singular times with Hausdorff dimension less than $\varepsilon$. As a byproduct, examples of anomalous dissipation in the class $L^{ {3}/{2} - \varepsilon}_t C^{ {1}/{3}} $ are given in both the viscous and inviscid case.


title: On a nonlinear relation for computing the overpartition function

abstract: In 1939, H. S. Zuckerman provided a Hardy-Ramanujan-Rademacher-type convergent series that can be used to compute an isolated value of the overpartition function $\overline{p}(n)$. Computing $\overline{p}(n)$ by this method requires arithmetic with very high-precision approximate real numbers and it is complicated. In this paper, we provide a formula to compute the values of $\overline{p}(n)$ that requires only the values of $\overline{p}(k)$ with $k\leqslant n/2$. This formula is combined with a known linear homogeneous recurrence relation for the overpartition function $\overline{p}(n)$ to obtain a simple and fast computation of the value of $\overline{p}(n)$. This new method uses only (large) integer arithmetic and it is simpler to program.


title: Resolving the scales of the Yang-Mills theory by means of an extra dimension

abstract: The Yang-Mills theory is part of the Standard Model of particle physics. The lack of the mathematical understanding of the theory stands out in theoretical physics. In order to address this problem we observe that a recently proposed general model beyond the Standard Model resolves the energy scales of a lattice regularized Yang-Mills theory by means of an extra dimension. The extra dimension ensures that all intermediate length scales of the physical system are available by definition. In this paper we study the role of the extra dimension also in the case of the free boson field. We find that if the extra dimension size is large, the model describes the classical motion of the system. In the opposite limit we recover its standard quantum mechanical motion without loss of information. Therefore, the Hilbert space of states in the presence of the extra dimension describes physical phenomena at all energy scales. This observation allows us to raise the description of length scales by an extra dimension at the level of a principle for the theories beyond the Standard Model, the only modeling constraint being the correspondence principle. As it was shown recently, the fermion ground state energy of a gauge invariant Hamilton operator of Dirac fermions gives a particular lattice regularization of the Yang-Mills theory. Integration of gauge fields gives a pure fermion theory of color singlet fermion-antifermion pairs at each lattice site evolving along the extra dimension. Color confinement follows directly from this property. It allows for a saddle point solution in the limit of a large number of colors. In this paper we find that the glueball spectrum of the Yang-Mills regularized theory is of the Hagedron type and bounded below by a positive value. We show also that the color charge is screened and the quark-antiquark potential is constant.


title: Collaborative Attention Mechanism for Multi-View Action Recognition

abstract: Multi-view action recognition (MVAR) leverages complementary temporal information from different views to improve the learning performance. Obtaining informative view-specific representation plays an essential role in MVAR. Attention has been widely adopted as an effective strategy for discovering discriminative cues underlying temporal data. However, most existing MVAR methods only utilize attention to extract representation for each view individually, ignoring the potential to dig latent patterns based on mutual-support information in attention space. To this end, we propose a collaborative attention mechanism (CAM) for solving the MVAR problem in this paper. The proposed CAM detects the attention differences among multi-view, and adaptively integrates frame-level information to benefit each other. Specifically, we extend the long short-term memory (LSTM) to a Mutual-Aid RNN (MAR) to achieve the multi-view collaboration process. CAM takes advantages of view-specific attention pattern to guide another view and discover potential information which is hard to be explored by itself. It paves a novel way to leverage attention information and enhances the multi-view representation learning. Extensive experiments on four action datasets illustrate the proposed CAM achieves better results for each view and also boosts multi-view performance.


title: SNR-enhanced diffusion MRI with structure-preserving low-rank denoising in reproducing kernel Hilbert spaces

abstract: Purpose: To introduce, develop, and evaluate a novel denoising technique for diffusion MRI that leverages non-linear redundancy in the data to boost the SNR while preserving signal information. Methods: We exploit non-linear redundancy of the dMRI data by means of Kernel Principal Component Analysis (KPCA), a non-linear generalization of PCAto reproducing kernel Hilbert spaces. By mapping the signal to a high-dimensional space, better redundancy is achieved despite nonlinearities in the data thereby enabling better denoising than linear PCA. We implement KPCA with a Gaussian kernel, with parameters automatically selected from knowledge of the noise statistics, and validate it on realistic Monte-Carlo simulations as well as with in-vivo human brain submillimeter resolution dMRI data. We demonstrate KPCA denoising using multi-coil dMRI data also. Results: SNR improvements up to 2.7 X were obtained in real in-vivo datasets denoised with KPCA, in comparison to SNR gains of up to 1.8 X when using state-of-the-art PCA denoising, e.g., Marchenko- Pastur PCA (MPPCA). Compared to gold-standard dataset references created from averaged data, we showed that lower normalized root mean squared error (NRMSE) was achieved with KPCA compared to MPPCA. Statistical analysis of residuals shows that only noise is removed. Improvements in the estimation of diffusion model parameters such as fractional anisotropy, mean diffusivity, and fiber orientation distribution functions (fODFs)were demonstrated. Conclusion:Non-linear redundancy of the dMRI signal can be exploited with KPCA, which allows superior noise reduction/ SNR improvements than state-of-the-art PCA methods, without loss of signal information.


title: The Efficient Preparation of Normal Distributions in Quantum Registers

abstract: The efficient preparation of input distributions is an important problem in obtaining quantum advantage in a wide range of domains. We propose a novel quantum algorithm for the efficient preparation of arbitrary normal distributions in quantum registers. To the best of our knowledge, our work is the first to leverage the power of Mid-Circuit Measurement and Reuse (MCMR), in a way that is broadly applicable to a range of state-preparation problems. Specifically, our algorithm employs a repeat-until-success scheme, and only requires a constant-bounded number of repetitions in expectation. In the experiments presented, the use of MCMR enables up to a 862.6x reduction in required qubits. Furthermore, the algorithm is provably resistant to both phase-flip and bit-flip errors, leading to a first-of-its-kind empirical demonstration on real quantum hardware, the MCMR-enabled Honeywell System Models H0 and H1-2.


title: VacSIM: Learning Effective Strategies for COVID-19 Vaccine Distribution using Reinforcement Learning

abstract: A COVID-19 vaccine is our best bet for mitigating the ongoing onslaught of the pandemic. However, vaccine is also expected to be a limited resource. An optimal allocation strategy, especially in countries with access inequities and temporal separation of hot-spots, might be an effective way of halting the disease spread. We approach this problem by proposing a novel pipeline VacSIM that dovetails Deep Reinforcement Learning models into a Contextual Bandits approach for optimizing the distribution of COVID-19 vaccine. Whereas the Reinforcement Learning models suggest better actions and rewards, Contextual Bandits allow online modifications that may need to be implemented on a day-to-day basis in the real world scenario. We evaluate this framework against a naive allocation approach of distributing vaccine proportional to the incidence of COVID-19 cases in five different States across India (Assam, Delhi, Jharkhand, Maharashtra and Nagaland) and demonstrate up to 9039 potential infections prevented and a significant increase in the efficacy of limiting the spread over a period of 45 days through the VacSIM approach. Our models and the platform are extensible to all states of India and potentially across the globe. We also propose novel evaluation strategies including standard compartmental model-based projections and a causality-preserving evaluation of our model. Since all models carry assumptions that may need to be tested in various contexts, we open source our model VacSIM and contribute a new reinforcement learning environment compatible with OpenAI gym to make it extensible for real-world applications across the globe. (http://vacsim.tavlab.iiitd.edu.in:8000/).


title: Noise and charge discreteness as ultimate limit for the THz operation of ultra-small electronic devices

abstract: To manufacture faster electron devices, the industry has entered into the nanoscale dimensions and Terahertz (THz) working frequencies. The discrete nature of the few electrons present simultaneously in the active region of ultra-small devices generate unavoidable fluctuations of the current at THz frequencies. The consequences of this noise remain unnoticed in the scientific community because its accurate understanding requires dealing with consecutive multi-time quantum measurements. Here, a modeling of the quantum measurement of the current at THz frequencies is introduced in terms of quantum (Bohmian) trajectories. With this new understanding, we develop an analytic model for THz noise as a function of the electron transit time and the sampling integration time, which finally determine the maximum device working frequency. The model is confirmed by either semi-classical or full-quantum time-dependent Monte Carlo simulations. All these results show that intrinsic THz noise increases unlimitedly when the volume of the active region decreases. All attempts to minimize the low signal-to-noise ratio of these ultra-small devices to get effective THz working frequencies are incompatible with the basic elements of the scaling strategy. One can develop THz electron devices, but they cannot have ultra-small dimensions. Or, one can fabricate ultra-small electron devices, but they cannot be used for THz working frequencies.


title: GIA-Net: Global Information Aware Network for Low-light Imaging

abstract: It is extremely challenging to acquire perceptually plausible images under low-light conditions due to low SNR. Most recently, U-Nets have shown promising results for low-light imaging. However, vanilla U-Nets generate images with artifacts such as color inconsistency due to the lack of global color information. In this paper, we propose a global information aware (GIA) module, which is capable of extracting and integrating the global information into the network to improve the performance of low-light imaging. The GIA module can be inserted into a vanilla U-Net with negligible extra learnable parameters or computational cost. Moreover, a GIA-Net is constructed, trained and evaluated on a large scale real-world low-light imaging dataset. Experimental results show that the proposed GIA-Net outperforms the state-of-the-art methods in terms of four metrics, including deep metrics that measure perceptual similarities. Extensive ablation studies have been conducted to verify the effectiveness of the proposed GIA-Net for low-light imaging by utilizing global information.


title: Probing the limits of vortex mode generation and detection with spatial light modulators

abstract: Spatial light modulators (SLMs) are popular tools for generating structured light fields and have fostered numerous applications in optics and photonics. Here, we explore the limits of what fields these devices are capable of generating and detecting in the context of so-called vortex beams carrying orbital angular momentum (OAM). Our main contributions are to quantify (theoretically and experimentally) how the pixelation of the SLM screen affects the quality of the generated vortex mode and to offer useful heuristics on how to optimise the performance of the displayed digital hologram. In so doing, we successfully generate and detect a very high order optical vortex mode with topological charge $\ell = 600$, the highest achieved to date using SLMs. Since the OAM degree of freedom is frequently touted as offering a potentially unbounded state space, we hope that this work will inspire researchers to make more use of higher order vortex modes.


title: Adaptive KL-UCB based Bandit Algorithms for Markovian and i.i.d. Settings

abstract: In the regret-based formulation of Multi-armed Bandit (MAB) problems, except in rare instances, much of the literature focuses on arms with i.i.d. rewards. In this paper, we consider the problem of obtaining regret guarantees for MAB problems in which the rewards of each arm form a Markov chain which may not belong to a single parameter exponential family. To achieve a logarithmic regret in such problems is not difficult: a variation of standard Kullback-Leibler Upper Confidence Bound (KL-UCB) does the job. However, the constants obtained from such an analysis are poor for the following reason: i.i.d. rewards are a special case of Markov rewards and it is difficult to design an algorithm that works well independent of whether the underlying model is truly Markovian or i.i.d. To overcome this issue, we introduce a novel algorithm that identifies whether the rewards from each arm are truly Markovian or i.i.d. using a total variation distance-based test. Our algorithm then switches from using a standard KL-UCB to a specialized version of KL-UCB when it determines that the arm reward is Markovian, thus resulting in low regrets for both i.i.d. and Markovian settings.


title: Has NANOGrav found first evidence for cosmic strings?

abstract: The North American Nanohertz Observatory for Gravitational Waves (NANOGrav) has recently reported strong evidence for a stochastic common-spectrum process affecting the pulsar timing residuals in its 12.5-year data set. We demonstrate that this process admits an interpretation in terms of a stochastic gravitational-wave background emitted by a cosmic-string network in the early Universe. We study stable Nambu-Goto strings in dependence of their tension $G\mu$ and loop size $\alpha$ and show that the entire viable parameter space will be probed by an array of future experiments.


title: Bayesian decomposition of the Galactic multi-frequency sky using probabilistic autoencoders

abstract: All-sky observations of the Milky Way show both Galactic and non-Galactic diffuse emission, for example from interstellar matter or the cosmic microwave background (CMB). The different emitters are partly superimposed in the measurements, partly they obscure each other, and sometimes they dominate within a certain spectral range. The decomposition of the underlying radiative components from spectral data is a signal reconstruction problem and often associated with detailed physical modeling and substantial computational effort. We aim to build an effective and self-instructing algorithm detecting the essential spectral information contained Galactic all-sky data covering spectral bands from $\gamma$-ray to radio waves. Utilizing principles from information theory, we develop a state-of-the-art variational autoencoder specialized on the adaption to Gaussian noise statistics. We first derive a generic generative process that leads from a low-dimensional set of emission features to the observed high-dimensional data. We formulate a posterior distribution of these features using Bayesian methods and approximate this posterior with variational inference. The algorithm efficiently encodes the information of 35 Galactic emission data sets in ten latent feature maps. These contain the essential information required to reconstruct the initial data with high fidelity and are ranked by the algorithm according to their significance for data regeneration. The three most significant feature maps encode astrophysical components: (1) The dense interstellar medium (ISM), (2) the hot and dilute regions of the ISM and (3) the CMB. The machine-assisted and data-driven dimensionality reduction of spectral data is able to uncover the physical features encoding the input data. Our algorithm is able to extract the dense and dilute Galactic regions, as well as the CMB, from the sky brightness values only.


title: Self-dual codes over $GF(q)$ with symmetric generator matrices

abstract: We introduce a consistent and efficient method to construct self-dual codes over $GF(q)$ with symmetric generator matrices from a self-dual code over $GF(q)$ of smaller length where $q \equiv 1 \pmod 4$. Using this method, we improve the best-known minimum weights of self-dual codes, which have not significantly improved for almost two decades. We focus on a class of self-dual codes, including double circulant codes. Using our method, called a `symmetric building-up' construction, we obtain many new self-dual codes over $GF(13)$ and $GF(17)$ and improve the bounds of best-known minimum weights of self-dual codes of lengths up to 40. Besides, we compute the minimum weights of quadratic residue codes that were not known before. These are: a [20,10,10] QR self-dual code over $GF(23)$, two [24,12,12] QR self-dual codes over $GF(29)$ and $GF(41)$, and a [32,12,14] QR self-dual codes over $GF(19)$. They have the highest minimum weights so far.


title: Adaptive Text Recognition through Visual Matching

abstract: In this work, our objective is to address the problems of generalization and flexibility for text recognition in documents. We introduce a new model that exploits the repetitive nature of characters in languages, and decouples the visual representation learning and linguistic modelling stages. By doing this, we turn text recognition into a shape matching problem, and thereby achieve generalization in appearance and flexibility in classes. We evaluate the new model on both synthetic and real datasets across different alphabets and show that it can handle challenges that traditional architectures are not able to solve without expensive retraining, including: (i) it can generalize to unseen fonts without new exemplars from them; (ii) it can flexibly change the number of classes, simply by changing the exemplars provided; and (iii) it can generalize to new languages and new characters that it has not been trained for by providing a new glyph set. We show significant improvements over state-of-the-art models for all these cases.


title: Knowledge-Based Legal Document Assembly

abstract: This paper proposes a knowledge-based legal document assembly method that uses a machine-readable representation of knowledge of legal professionals. This knowledgebase has two components - the formal knowledge of legal norms represented as a rule-base and the tacit knowledge represented by a document template. A document assembly system is developed as a proof of concept. It collects input data in the form of an interactive interview, performs legal reasoning over input data, and generates the output document. The system also creates an argument graph as an explanation of the reasoning process providing the user with an interpretation of how the input data and the rule-base influence the content of the output document. The system also semantically marks up data in the output document, facilitating its further processing and providing support to the interoperability of information systems in the legal domain.


title: The powers of two as sums over partitions

abstract: In this paper, we investigate two methods to express the natural powers of $2$ as sums over integer partitions. First we consider a formula by N. J. Fine that allows us to express a binomial coefficient in terms of multinomial coefficients as a sum over partitions. The second method invokes the central binomial coefficients and the logarithmic differentiation of their generating function. Some experimental results suggest the existence of other methods of decomposing the power of $2$ as sums over partitions.


title: High-Resolution Deep Image Matting

abstract: Image matting is a key technique for image and video editing and composition. Conventionally, deep learning approaches take the whole input image and an associated trimap to infer the alpha matte using convolutional neural networks. Such approaches set state-of-the-arts in image matting; however, they may fail in real-world matting applications due to hardware limitations, since real-world input images for matting are mostly of very high resolution. In this paper, we propose HDMatt, a first deep learning based image matting approach for high-resolution inputs. More concretely, HDMatt runs matting in a patch-based crop-and-stitch manner for high-resolution inputs with a novel module design to address the contextual dependency and consistency issues between different patches. Compared with vanilla patch-based inference which computes each patch independently, we explicitly model the cross-patch contextual dependency with a newly-proposed Cross-Patch Contextual module (CPC) guided by the given trimap. Extensive experiments demonstrate the effectiveness of the proposed method and its necessity for high-resolution inputs. Our HDMatt approach also sets new state-of-the-art performance on Adobe Image Matting and AlphaMatting benchmarks and produce impressive visual results on more real-world high-resolution images.


title: The nonperturbative contribution to asymptotic masses

abstract: In gauge theories, charged particles obey modified dispersion due to medium interactions (forward scattering), leading at high energies $E \geq T$ to an asymptotic mass-squared $m^2_\infty$. We calculate the infrared part of this mass nonperturbatively for the theory of the strong interactions, QCD, through a lattice treatment of its low-energy effective description, Electrostatic QCD (EQCD). Incorporation of these results into a nonperturbative determination of the effective thermal mass will require a still-incomplete next-to-leading order perturbative matching of this quantity to full QCD.


title: Multilevel regression with poststratification for the national level Viber/Street poll on the 2020 presidential election in Belarus

abstract: Independent sociological polls are forbidden in Belarus. Online polls performed without sound scientific rigour do not yield representative results. Yet, both inside and outside Belarus it is of great importance to obtain precise estimates of the ratings of all candidates. These ratings could function as reliable proxies for the election's outcomes. We conduct an independent poll based on the combination of the data collected via Viber and on the streets of Belarus. The Viber and the street data samples consist of almost 45000 and 1150 unique observations respectively. Bayesian regressions with poststratification were build to estimate ratings of the candidates and rates of early voting turnout for the population as a whole and within various focus subgroups. We show that both the officially announced results of the election and early voting rates are highly improbable. With a probability of at least 95%, Sviatlana Tikhanouskaya's rating lies between 75% and 80%, whereas Aliaksandr Lukashenka's rating lies between 13% and 18% and early voting rate predicted by the method ranges from 9% to 13% of those who took part in the election. These results contradict the officially announced outcomes, which are 10.12%, 80.11%, and 49.54% respectively and lie far outside even the 99.9% credible intervals predicted by our model. The only marginal groups of people where the upper bounds of the 99.9% credible intervals of the rating of Lukashenka are above 50% are people older than 60 and uneducated people. For all other marginal subgroups, including rural residents, even the upper bounds of 99.9% credible intervals for Lukashenka are far below 50%. The same is true for the population as a whole. Thus, with a probability of at least 99.9% Lukashenka could not have had enough electoral support to win the 2020 presidential election in Belarus.


title: X-ray Spectroscopy of the Highly Magnetized Pulsar PSR J1846-0258, its Wind Nebula and Hosting Supernova Remnant Kes 75

abstract: We present broad-band X-ray spectroscopy of the energetic components that make up the supernova remnant (SNR) Kesteven 75 using concurrent 2017 Aug 17-20 XMM-Newton and NuSTAR observations, during which the pulsar PSR J1846-0258 is found to be in the quiescent state. The young remnant hosts a bright pulsar wind nebula powered by the highly-energetic (Edot = 8.1E36 erg/s) isolated, rotation-powered pulsar, with a spin-down age of only P/2Pdot ~ 728 yr. Its inferred magnetic field (Bs = 4.9E13 G) is the largest known for these objects, and is likely responsible for intervals of flare and burst activity, suggesting a transition between/to a magnetar state. The pulsed emission from PSR J1846-0258 is well-characterized in the 2-50 keV range by a power-law model with photon index Gamma_PSR = 1.24+/-0.09 and a 2-10 keV unabsorbed flux of (2.3+/-0.4)E-12 erg/s/cm^2). We find no evidence for an additional non-thermal component above 10 keV in the current state, as would be typical for a magnetar. Compared to the Chandra pulsar spectrum, the intrinsic pulsed fraction is 71+/-16% in 2-10 keV band. A power-law spectrum for the PWN yields Gamma_PWN = 2.03+/-0.03 in the 1-55 keV band, with no evidence of curvature in this range, and a 2-10 keV unabsorbed flux (2.13+/-0.02)E-11 erg/s/cm^2. The NuSTAR data reveal evidence for a hard X-ray component dominating the SNR spectrum above 10 keV which we attribute to a dust-scattered PWN component. We model the dynamical and radiative evolution of the Kes 75 system to estimate the birth properties of the neutron star, the energetics of its progenitor, and properties of the PWN. This suggests that the progenitor of Kes 75 was originally in a binary system which transferred most its mass to a companion before exploding.


title: Increasing the Accuracy of Binary Neutron Star Simulations with an improved Vacuum Treatment

abstract: Numerical relativity simulations are essential to study the last stages of the binary neutron star coalescence. Unfortunately, for stable simulations there is the need to add an artificial low-density atmosphere. Here we discuss a new framework in which we can effectively set the density surrounding the neutron stars to zero to ensure a more accurate simulation. We test our method with a number of single star test cases and for an equal mass binary neutron star simulation. While the bulk motion of the system is not influenced, and hence, there is no improvement with respect to the emitted gravitational-wave signal, we find that the new approach is superior with respect to mass conservation and it allows a much better tracking of outward moving material. This will allow a more accurate simulation of the ejected material and supports the interpretation of present and future multi-messenger observations with more accurate numerical relativity simulations.


title: Passive communication with ambient noise

abstract: Motivated by applications to wireless communications, this paper addresses the propagation of waves transmitted by ambient noise sources and interacting with metamaterials. We discuss a generalized Helmholtz-Kirchhoff identity that is valid in dispersive media and we characterize the statistical properties of the empirical cross spectral density of the wave field. We can then introduce and analyze an original communication scheme between two passive arrays that uses only ambient noise illumination. The passive transmitter array does not transmit anything but it is a tunable metamaterial surface that can modulate its dispersive properties and encode a message in the modulation. The passive receiver array made of two receivers that are half-a-wavelength apart from each other can decode the message from the empirical cross spectral density of the wave field.


title: Observation of multiplicity-dependent prompt $\chi_{c1}(3872)$ and $\psi(2S)$ production in $pp$ collisions

abstract: The production of $\chi_{c1}(3872)$ and $\psi(2S)$ hadrons is studied as a function of charged particle multiplicity in $pp$ collisions at a center-of-mass energy of 8 TeV, corresponding to an integrated luminosity of 2 fb$^{-1}$. For both states, the fraction that is produced promptly at the collision vertex is found to decrease as charged particle multiplicity increases. The ratio of $\chi_{c1}(3872)$ to $\psi(2S)$ cross-sections for promptly produced particles is also found to decrease with multiplicity, while no significant dependence on multiplicity is observed for the equivalent ratio of particles produced away from the collision vertex in $b$-hadron decays. This behavior is consistent with a calculation that models the $\chi_{c1}(3872)$ structure as a compact tetraquark. Comparisons with model calculations and implications for the binding energy of the $\chi_{c1}(3872)$ state are discussed.


title: Chiral current circulation and $\mathcal{PT}$ symmetry in a trimer of oscillators

abstract: We present a simple quantum theory of a bosonic trimer in a triangular configuration, subject to gain and loss in an open quantum systems approach. Importantly, the coupling constants between each oscillator are augmented by complex arguments, which give rise to various asymmetries. In particular, one may tune the complex phases to induce chiral currents, including the special case of completely unidirectional (or one-way) circulation when certain conditions are met regarding the coherent and incoherent couplings. When our general theory is recast into a specific non-Hermitian Hamiltonian, we find interesting features in the trimer population dynamics close to the exceptional points between phases of broken and unbroken $\mathcal{PT}$ symmetry. Our theoretical work provides perspectives for the experimental realization of chiral transport at the nanoscale in a variety of accessible nanophotonic and nanoplasmonic systems, and paves the way for the potential actualization of nonreciprocal devices.


title: The Frisch--Waugh--Lovell Theorem for Standard Errors

abstract: The Frisch--Waugh--Lovell Theorem states the equivalence of the coefficients from the full and partial regressions. I further show the equivalence between various standard errors. Applying the new result to stratified experiments reveals the discrepancy between model-based and design-based standard errors.


title: On the impact of galaxy bias uncertainties on primordial non-Gaussianity constraints

abstract: We study the impact that uncertainties on assumed relations between galaxy bias parameters have on constraints of the local PNG $f_{\rm NL}$ parameter. We focus on the relation between the linear density galaxy bias $b_1$ and local PNG bias $b_\phi$ in an idealized forecast setup with multitracer galaxy power spectrum and bispectrum data. We consider two parametrizations of galaxy bias: 1) one inspired by the universality relation where $b_\phi = 2\delta_c\left(b_1 - p\right)$ and $p$ is a free parameter; and 2) another in which the product of bias parameters and $f_{\rm NL}$, like $f_{\rm NL} b_\phi$, is directly fitted for. The constraints on the $f_{\rm NL}-p$ plane are markedly bimodal, and both the central value and width of marginalized constraints on $f_{\rm NL}$ depend sensitively on the priors on $p$. Assuming fixed $p=1$ in the constraints with a fiducial value of $p=0.55$ can bias the inferred $f_{\rm NL}$ by $0.5\sigma$ to $1\sigma$; priors $\Delta p \approx 0.5$ around this fiducial value are however sufficient in our setup to return unbiased constraints. In power spectrum analyses, parametrization 2, that makes no assumptions on $b_\phi$, can distinguish $f_{\rm NL} \neq 0$ with the same significance as parametrization 1 assuming perfect knowledge of $b_\phi$ (the value of $f_{\rm NL}$ is however left unknown). A drawback of parametrization 2 is that the addition of the bispectrum information is not as beneficial as in parametrization 1. Our results motivate strongly the incorporation of mitigation strategies for bias uncertainties in PNG constraint analyses, as well as further theoretical studies on the relations between bias parameters to better inform those strategies.


title: Suppression of Superconductivity and Nematic Order in Fe$_{1-y}$Se$_{1-x}$S$_x$ (0$\leq$$x$$\leq$1, $y$$\leq$0.1) Crystals by Anion Height Disorder

abstract: Connections between crystal chemistry and critical temperature $T_c$ have been in the focus of superconductivity, one of the most widely studied phenomena in physics, chemistry and materials science alike. In most Fe-based superconductors, materials chemistry and physics conspire so that $T_c$ correlates with the average anion height above the Fe plane, i. e. with the geometry of the FeAs4 or FeCh4 (Ch = Te, Se, or S) tetrahedron. By synthesizing Fe$_{1-y}$Se$_{1-x}$S$_x$ (0$\leq$$x$$\leq$1, $y$$\leq$0.1), we find that in alloyed crystals $T_c$ is not correlated with the anion height like it is for most other Fe superconductors. Instead, changes in $T_c$($x$) and tetragonal-to-orthorombic (nematic) transition $T_s$($x$) upon cooling are correlated with disorder in Fe vibrations in the direction orthogonal to Fe planes, along the crystallographic c-axis. The disorder stems from the random nature of S substitution, causing deformed Fe(Se,S)4 tetrahedra with different Fe-Se and Fe-S bond distances. Our results provide evidence of $T_c$ and $T_s$ suppression by disorder in anion height. The connection to local crystal chemistry may be exploited in computational prediction of new superconducting materials with Fe/S building blocks.


title: Beamed and unbeamed emission of $\gamma$-ray blazars

abstract: A two-component model of radio emission has been used to explain some radio observational properties of Active Galactic Nuclei (AGNs) and, in particular, of blazars. In this work, we extend the two-component idea to the $\gamma$-ray emission and assume that the total $\gamma$-ray output of blazars consists of relativistically beamed and unbeamed components. The basic idea leverages the correlation between the radio core-dominance parameter and the $\gamma$-ray beaming factor. To do so, we evaluate this correlation for a large sample of 584 blazars taken from the fourth source catalog of the Fermi Large Area Telescope (Fermi-LAT) and correlated their $\gamma$-ray core-dominance parameters with radio core-dominance parameters. The $\gamma$-ray beaming factor is then used to estimate the beamed and unbeamed components. Our analysis confirms that the $\gamma$-ray emission in blazars is mainly from the beamed component.


title: Revealing Secrets in SPARQL Session Level

abstract: Based on Semantic Web technologies, knowledge graphs help users to discover information of interest by using live SPARQL services. Answer-seekers often examine intermediate results iteratively and modify SPARQL queries repeatedly in a search session. In this context, understanding user behaviors is critical for effective intention prediction and query optimization. However, these behaviors have not yet been researched systematically at the SPARQL session level. This paper reveals secrets of session-level user search behaviors by conducting a comprehensive investigation over massive real-world SPARQL query logs. In particular, we thoroughly assess query changes made by users w.r.t. structural and data-driven features of SPARQL queries. To illustrate the potentiality of our findings, we employ an application example of how to use our findings, which might be valuable to devise efficient SPARQL caching, auto-completion, query suggestion, approximation, and relaxation techniques in the future.


title: Optimal Bounds on Nonlinear Partial Differential Equations in Model Certification, Validation, and Experimental Design

abstract: We demonstrate that the recently developed Optimal Uncertainty Quantification (OUQ) theory, combined with recent software enabling fast global solutions of constrained non-convex optimization problems, provides a methodology for rigorous model certification, validation, and optimal design under uncertainty. In particular, we show the utility of the OUQ approach to understanding the behavior of a system that is governed by a partial differential equation -- Burgers' equation. We solve the problem of predicting shock location when we only know bounds on viscosity and on the initial conditions. Through this example, we demonstrate the potential to apply OUQ to complex physical systems, such as systems governed by coupled partial differential equations. We compare our results to those obtained using a standard Monte Carlo approach, and show that OUQ provides more accurate bounds at a lower computational cost. We discuss briefly about how to extend this approach to more complex systems, and how to integrate our approach into a more ambitious program of optimal experimental design.


title: User Manual for the SU2 EQUiPS Module: Enabling Quantification of Uncertainty in Physics Simulations

abstract: This document serves as the manual for using the EQUiPS (Enabling Quantification of Uncertainty in Physics Simulations) module in SU2. The EQUiPS module uses the Eigenspace Perturbation methodology to provide interval bounds on Quantities of Interest (QoIs) that capture epistemic uncertainties arising from assumptions made in RANS turbulence models. This has been implemented and tested in SU2 for a variety of benchmark turbulence cases as well as flows of aerodynamic interest.


title: A fully-coupled framework for solving Cahn-Hilliard Navier-Stokes equations: Second-order, energy-stable numerical methods on adaptive octree based meshes

abstract: We present a fully-coupled, implicit-in-time framework for solving a thermodynamically-consistent Cahn-Hilliard Navier-Stokes system that models two-phase flows. In this work, we extend the block iterative method presented in Khanwale et al. [\textit{Simulating two-phase flows with thermodynamically consistent energy stable Cahn-Hilliard Navier-Stokes equations on parallel adaptive octree based meshes}, J. Comput. Phys. (2020)], to a fully-coupled, provably second-order accurate scheme in time, while maintaining energy-stability. The new method requires fewer matrix assemblies in each Newton iteration resulting in faster solution time. The method is based on a fully-implicit Crank-Nicolson scheme in time and a pressure stabilization for an equal order Galerkin formulation. That is, we use a conforming continuous Galerkin (cG) finite element method in space equipped with a residual-based variational multiscale (RBVMS) procedure to stabilize the pressure. We deploy this approach on a massively parallel numerical implementation using parallel octree-based adaptive meshes. We present comprehensive numerical experiments showing detailed comparisons with results from the literature for canonical cases, including the single bubble rise, Rayleigh-Taylor instability, and lid-driven cavity flow problems. We analyze in detail the scaling of our numerical implementation.


title: An explainable XGBoost-based approach towards assessing the risk of cardiovascular disease in patients with Type 2 Diabetes Mellitus

abstract: Cardiovascular Disease (CVD) is an important cause of disability and death among individuals with Diabetes Mellitus (DM). International clinical guidelines for the management of Type 2 DM (T2DM) are founded on primary and secondary prevention and favor the evaluation of CVD related risk factors towards appropriate treatment initiation. CVD risk prediction models can provide valuable tools for optimizing the frequency of medical visits and performing timely preventive and therapeutic interventions against CVD events. The integration of explainability modalities in these models can enhance human understanding on the reasoning process, maximize transparency and embellish trust towards the models' adoption in clinical practice. The aim of the present study is to develop and evaluate an explainable personalized risk prediction model for the fatal or non-fatal CVD incidence in T2DM individuals. An explainable approach based on the eXtreme Gradient Boosting (XGBoost) and the Tree SHAP (SHapley Additive exPlanations) method is deployed for the calculation of the 5-year CVD risk and the generation of individual explanations on the model's decisions. Data from the 5-year follow up of 560 patients with T2DM are used for development and evaluation purposes. The obtained results (AUC = 71.13%) indicate the potential of the proposed approach to handle the unbalanced nature of the used dataset, while providing clinically meaningful insights about the ensemble model's decision process.


title: Advanced Virus Monitoring and Analysis System

abstract: This research proposed an architecture and a system which able to monitor the virus behavior and classify them as a traditional or polymorphic virus. Preliminary research was conducted to get the current virus behavior and to find the certain parameters which usually used by virus to attack the computer target. Finally, test bed environment is used to test our system by releasing the virus in a real environment, and try to capture their behavior, and followed by generating the conclusion that the tested or monitored virus is classified as a traditional or polymorphic virus.


title: Large scale simulations of H and He reionization and heating driven by stars and more energetic sources

abstract: We present simulations of cosmic reionization and reheating from $z=18$ to $z=5$, investigating the role of stars (emitting soft UV-photons), nuclear black holes (BHs, with power-law spectra), X-ray binaries (XRBs, with hard X-ray dominated spectra), and the supernova-associated thermal bremsstrahlung of the diffuse interstellar medium (ISM, with soft X-ray spectra). We post-process the hydrodynamical simulation Massive-Black II (MBII) with multifrequency ionizing radiative transfer. The source properties are directly derived from the physical environment of MBII, and our only real free parameter is the ionizing escape fraction $f_{\rm esc}$. We find that, among the models explored here, the one with an escape fraction that decreases with decreasing redshift yields results most in line with observations, such as of the neutral hydrogen fraction and the Thomson scattering optical depth. Stars are the main driver of hydrogen reionization and consequently of the thermal history of the intergalactic medium (IGM). We obtain $\langle x_{\rm HII} \rangle = 0.99998$ at $z=6$ for all source types, with volume averaged temperatures $\langle T \rangle \sim 20,000~{\rm K}$. BHs are rare and negligible to hydrogen reionization, but conversely they are the only sources which can fully ionize helium, increasing local temperatures by $\sim 10^4~{\rm K}$. The thermal and ionization state of the neutral and lowly ionized hydrogen differs significantly with different source combinations, with ISM and (to a lesser extent) XRBs, playing a significant role and, as a consequence, determining the transition from absorption to emission of the 21 cm signal from neutral hydrogen.


title: Quantum smectic gauge theory

abstract: We present a gauge theory formulation of a two-dimensional quantum smectic and its relatives, motivated by their realizations in correlated quantum matter. The description gives a unified treatment of phonons and topological defects, respectively encoded in a pair of coupled gauge fields and corresponding charges. The charges exhibit subdimensional constrained quantum dynamics and anomalously slow highly anisotropic diffusion of disclinations inside a smectic. This approach gives a transparent description of a multi-stage quantum melting transition of a two-dimensional commensurate crystal (through an incommensurate crystal - a supersolid) into a quantum smectic, that subsequently melts into a quantum nematic and isotropic superfluids, all in terms of a sequence of Higgs transitions.


title: Socially competent robots: adaptation improves leadership performance in groups of live fish

abstract: Collective motion is commonly modeled with simple interaction rules between agents. Yet in nature, numerous observables vary within and between individuals and it remains largely unknown how animals respond to this variability, and how much of it may be the result of social responses. Here, we hypothesize that Guppies (\textit{Poecilia reticulata}) respond to avoidance behaviors of their shoal mates and that "socially competent" responses allow them to be more effective leaders. We test this hypothesis in an experimental setting in which a robotic Guppy, called RoboFish, is programmed to adapt to avoidance reactions of its live interaction partner. We compare the leadership performance between socially competent robots and two non-competent control behaviors and find that 1) behavioral variability itself appears attractive and that socially competent robots are better leaders that 2) require fewer approach attempts to 3) elicit longer average following behavior than non-competent agents. This work provides evidence that social responsiveness to avoidance reactions plays a role in the social dynamics of guppies. We showcase how social responsiveness can be modeled and tested directly embedded in a living animal model using adaptive, interactive robots.


title: Correlations between supermassive black holes and hot gas atmospheres in IllustrisTNG and X-ray observations

abstract: Recent X-ray observations have revealed remarkable correlations between the masses of central supermassive black holes (SMBHs) and the X-ray properties of the hot atmospheres permeating their host galaxies, thereby indicating the crucial role of the atmospheric gas in tracing SMBH growth in the high-mass regime. We examine this topic theoretically using the IllustrisTNG cosmological simulations and provide insights to the nature of this SMBH -- gaseous halo connection. By carrying out a mock X-ray analysis for a mass-selected sample of TNG100 simulated galaxies at $z=0$, we inspect the relationship between the masses of SMBHs and the hot gas temperatures and luminosities at various spatial and halo scales -- from galactic ($\sim R_{\rm e}$) to group/cluster scales ($\sim R_{\rm 500c}$). We find strong SMBH-X-ray correlations mostly in quenched galaxies and find that the correlations become stronger and tighter at larger radii. Critically, the X-ray temperature ($k_{\rm B}T_{\rm X}$) at large radii ($r\gtrsim5R_{\rm e}$) traces the SMBH mass with a remarkably small scatter ($\sim0.2$ dex). The relations emerging from IllustrisTNG are broadly consistent with those obtained from recent X-ray observations. Overall, our analysis suggests that, within the framework of IllustrisTNG, the present-time $M_{\rm BH}-k_{\rm B}T_{X}$ correlations at the high-mass end ($M_{\rm BH}\gtrsim 10^8 M_\odot$) are fundamentally a reflection of the SMBH mass -- halo mass relation, which at such high masses is set by the hierarchical assembly of structures. The exact form, locus, and scatter of those scaling relations are, however, sensitive to feedback processes such as those driven by star formation and SMBH activity.


title: PDFFlow: parton distribution functions on GPU

abstract: We present PDFFlow, a new software for fast evaluation of parton distribution functions (PDFs) designed for platforms with hardware accelerators. PDFs are essential for the calculation of particle physics observables through Monte Carlo simulation techniques. The evaluation of a generic set of PDFs for quarks and gluon at a given momentum fraction and energy scale requires the implementation of interpolation algorithms as introduced for the first time by the LHAPDF project. PDFFlow extends and implements these interpolation algorithms using Google's TensorFlow library providing the capabilities to perform PDF evaluations taking fully advantage of multi-threading CPU and GPU setups. We benchmark the performance of this library on multiple scenarios relevant for the particle physics community.


title: Eccentric Tidal Disruption Event Disks around Supermassive Black Holes: Dynamics and Thermal Emission

abstract: After the Tidal Disruption Event (TDE) of a star around a SuperMassive Black Hole (SMBH), if the stellar debris stream rapidly circularizes and forms a compact disk, the TDE emission is expected to peak in the soft X-ray or far Ultra-Violet (UV). The fact that many TDE candidates are observed to peak in the near UV and optical has challenged conventional TDE emission models. By idealizing a disk as a nested sequence of elliptical orbits which communicate adiabatically via pressure forces, and are heated by energy dissipated during the circularization of the nearly parabolic debris streams, we investigate the dynamics and thermal emission of highly eccentric TDE disks, including the effect of General-Relativistic apsidal precession from the SMBH. We calculate the properties of uniformly precessing, apsidally aligned, and highly eccentric TDE disks, and find highly eccentric disk solutions exist for realistic TDE properties (SMBH and stellar mass, periapsis distance, etc.). Taking into account compressional heating (cooling) near periapsis (apoapsis), we find our idealized eccentric disk model can produce emission consistent with the X-ray and UV/Optical luminosities of many optically bright TDE candidates. Our work attempts to quantify the thermal emission expected from the shock-heating model for TDE emission, and finds stream-stream collisions are a promising way to power optically bright TDEs.


title: Non-perturbative unitarity and fictitious ghosts in quantum gravity

abstract: We discuss aspects of non-perturbative unitarity in quantum field theory. The additional ghost degrees of freedom arising in "truncations" of an effective action at a finite order in derivatives could be fictitious degrees of freedom. Their contributions to the fully-dressed propagator -- the residues of the corresponding ghost-like poles -- vanish once all operators compatible with the symmetry of the theory are included in the effective action. These "fake ghosts" do not indicate a violation of unitarity.


title: An Increase in Small-planet Occurrence with Metallicity for Late-type Dwarf Stars in the Kepler Field and Its Implications for Planet Formation

abstract: While it is well established that giant-planet occurrence rises rapidly with host star metallicity, it is not yet clear if small-planet occurrence around late-type dwarf stars depends on host star metallicity. Using the Kepler Data Release 25 planet candidate list and its completeness data products, we explore planet occurrence as a function of metallicity in the Kepler field's late-type dwarf stellar population. We find that planet occurrence increases with metallicity for all planet radii $R_{\mathrm{p}}$ down to at least $R_{\mathrm{p}}\approx2~R_{\oplus}$ and that in the range $2~R_{\oplus}\lesssim R_{\mathrm{p}}\lesssim 5~R_{\oplus}$ planet occurrence scales linearly with metallicity $Z$. Extrapolating our results, we predict that short-period planets with $R_{\mathrm{p}}\lesssim 2~R_{\oplus}$ should be rare around early M dwarf stars with $\mathrm{[M/H]}\lesssim-0.5$ or late M dwarf stars with $\mathrm{[M/H]}\lesssim+0.0$. This dependence of planet occurrence on metallicity observed in the Kepler field emphasizes the need to control for metallicity in estimates of planet occurrence for late-type dwarf stars like those targeted by Kepler's K2 extension and the Transiting Exoplanet Survey Satellite (TESS). We confirm the theoretical expectation that the small planet occurrence--host star metallicity relation is stronger for low-mass stars than for solar-type stars. We establish that the expected solid mass in planets around late-type dwarfs in the Kepler field is comparable to the total amount of planet-making solids in their protoplanetary disks. We argue that this high efficiency of planet formation favors planetesimal accretion over pebble accretion as the origin of the small planets observed by Kepler around late-type dwarf stars.


title: Extracting the Subhalo Mass Function from Strong Lens Images with Image Segmentation

abstract: Detecting substructure within strongly lensed images is a promising route to shed light on the nature of dark matter. However, it is a challenging task, which traditionally requires detailed lens modeling and source reconstruction, taking weeks to analyze each system. We use machine-learning to circumvent the need for lens and source modeling and develop a neural network to both locate subhalos in an image as well as determine their mass using the technique of image segmentation. The network is trained on images with a single subhalo located near the Einstein ring across a wide range of apparent source magnitudes. The network is then able to resolve subhalos with masses $m\gtrsim 10^{8.5} M_{\odot}$. Training in this way allows the network to learn the gravitational lensing of light, and remarkably, it is then able to detect entire populations of substructure, even for locations further away from the Einstein ring than those used in training. Over a wide range of the apparent source magnitude, the false-positive rate is around three false subhalos per 100 images, coming mostly from the lightest detectable subhalo for that signal-to-noise ratio. With good accuracy and a low false-positive rate, counting the number of pixels assigned to each subhalo class over multiple images allows for a measurement of the subhalo mass function (SMF). When measured over three mass bins from $10^9M_{\odot}$--$10^{10} M_{\odot}$ the SMF slope is recovered with an error of 36% for 50 images, and this improves to 10% for 1000 images with Hubble Space Telescope-like noise.

categories: astro-ph.CO astro-ph.IM hep-ph physics.data-an stat.ML


title: What becomes of vortices when they grow giant

abstract: We discuss vortex solutions of the abelian Higgs model in the limit of large winding number $n$. We suggest a framework where a topological quantum number $n$ is associated with a ratio of dynamical scales and a systematic expansion in inverse powers of $n$ is then derived in the spirit of effective field theory. The general asymptotic form of ${\it giant}$ vortices is obtained. For critical coupling the axially symmetric vortices become ${\it integrable}$ in the large-$n$ limit and we present the corresponding analytic solution. The method provides simple asymptotic formulae for the vortex shape and parameters with accuracy that can be systematically improved, and can be applied to topological solitons of other models. After including the next-to-leading terms the approximation works remarkably well down to $n=1$.


title: Black Hole Glimmer Signatures of Mass, Spin, and Inclination

abstract: Gravitational lensing near a black hole is strong enough that light rays can circle the event horizon multiple times. Photons emitted in multiple directions at a single event, perhaps because of localized, impulsive heating of accreting plasma, take multiple paths to a distant observer. In the Kerr geometry, each path is associated with a distinct light travel time and a distinct arrival location in the image plane, producing black hole glimmer. This sequence of arrival times and locations uniquely encodes the mass and spin of the black hole and can be understood in terms of properties of bound photon orbits. We provide a geometrically motivated treatment of Kerr glimmer and evaluate it numerically for simple hotspot models to show that glimmer can be measured in a finite-resolution observation. We discuss potential measurement methods and implications for tests of the Kerr hypothesis.


title: OGLE-GAL-ACEP-091 -- The First Known Multi-Mode Anomalous Cepheid

abstract: Anomalous Cepheids (ACs) are metal-deficient, core-helium-burning pulsating stars with masses in the range 1.2-2.2 M_S. Until recently, all known ACs were pure single-mode pulsators. The first candidate for an AC pulsating in more than one radial mode - OGLE-GAL-ACEP-091 - was recently identified in the Milky Way based on the photometric database of the Optical Gravitational Lensing Experiment (OGLE) survey. We analyze this object showing that it is actually a triple-mode pulsator. Its position in the Petersen diagram, the light-curve morphology quantified by Fourier coefficients, and absolute magnitudes derived from the Gaia parallax are consistent with the assumption that OGLE-GAL-ACEP-091 is an AC. Our grid of linear pulsation models indicates that OGLE-GAL-ACEP-091 is a 1.8 M_S star with a metallicity of about [Fe/H] = -0.5 dex.