title: SCARLET-NAS: Bridging the Gap between Stability and Scalability in Weight-sharing Neural Architecture Search

abstract: To discover powerful yet compact models is an important goal of neural architecture search. Previous two-stage one-shot approaches are limited by search space with a fixed depth. It seems handy to include an additional skip connection in the search space to make depths variable. However, it creates a large range of perturbation during supernet training and it has difficulty giving a confident ranking for subnetworks. In this paper, we discover that skip connections bring about significant feature inconsistency compared with other operations, which potentially degrades the supernet performance. Based on this observation, we tackle the problem by imposing an equivariant learnable stabilizer to homogenize such disparities. Experiments show that our proposed stabilizer helps to improve the supernet's convergence as well as ranking performance. With an evolutionary search backend that incorporates the stabilized supernet as an evaluator, we derive a family of state-of-the-art architectures, the SCARLET series of several depths, especially SCARLET-A obtains 76.9% top-1 accuracy on ImageNet. Code is available at https://github.com/xiaomi-automl/ScarletNAS.


title: The physical origin of the cosmological constant: continuum limit and the analogy with the Casimir effect

abstract: In this paper we continue the investigations in \cite{1} concerning the origin of the cosmological constant. First of all, we generalise the results in \cite{1,2} by considering a continuum approximation for a radiation field in a cosmological background. In this way, we clearly show that the bare cosmological constant is obtained with wanishing temperature $T$ and that the specific heat $C$ is zero at the decoherence scale $L_D$. Moreover, we address the issue to fix the parameters present in our model. In particular, we push forward the analogy between our expression for ${\Lambda}_L$ at a given proper scale $L$ and the one extrapolated by the Casimir effect. As a consequence, we can fix the decoherence scale $L_D$ at which we have the crossover to classicality to be of the order of $\sim 10^{-5}$ meters. This implies that the actual observed value of the cosmological constant is fixed (frozen) at this new physical scale.


title: Tackling Online Abuse: A Survey of Automated Abuse Detection Methods

abstract: Abuse on the Internet represents an important societal problem of our time. Millions of Internet users face harassment, racism, personal attacks, and other types of abuse on online platforms. The psychological effects of such abuse on individuals can be profound and lasting. Consequently, over the past few years, there has been a substantial research effort towards automated abuse detection in the field of natural language processing (NLP). In this paper, we present a comprehensive survey of the methods that have been proposed to date, thus providing a platform for further development of this area. We describe the existing datasets and review the computational approaches to abuse detection, analyzing their strengths and limitations. We discuss the main trends that emerge, highlight the challenges that remain, outline possible solutions, and propose guidelines for ethics and explainability


title: Universality of scaling of correlations across probability distributions

abstract: Scale invariance and the resulting power law behaviours are seen in diverse systems. In this work we consider translation, rotational and scale invariant systems defined on a lattice, such that the variables defining the state at every lattice site take on the same range of finite values, with these values collectively picked up from probability distribution that can be arbitrary. We show that the exponent that describes the scaling of the two point correlation function in these systems will match the scaling exponent of a equilibrium statistical mechanical model described by a Boltzmannian distribution at criticality. This work therefore extends the concept of universality in statistical mechanics to probability distributions that do not have a Boltzmannian form.


title: Entire functions with prescribed singular values

abstract: We introduce a new class of entire functions $\mathcal{E}$ which consists of all $F_0\in\mathcal{O}(\mathbb{C})$ for which there exists a sequence $(F_n)\in \mathcal{O}(\mathbb{C})$ and a sequence $(\lambda_n)\in\mathbb{C}$ satisfying $F_n(z)=\lambda_{n+1}e^{F_{n+1}(z)}$ for all $n\geq 0$. This new class is closed under the composition and its is dense in the space of all non-vanishing entire functions. We prove that every closed set $V\subset \mathbb{C}$ containing the origin and at least one more point is the set of singular values of some locally univalent function in $\mathcal{E}$, hence this new class has non-trivial intersection with both the Speiser class and the Eremenko-Lyubich class of entire functions. As a consequence we provide a new proof of an old result by Heins which states that every closed set $V\subset\mathbb{C}$ is the set of singular values of some locally univalent entire function. The novelty of our construction is that these functions are obtained as a uniform limit of a sequence of entire functions, the process under which the set of singular values is not stable. Finally we show that the class $\mathcal{E}$ contains functions with an empty Fatou set and also functions whose Fatou set is non-empty.


title: Dispersion in an array of buildings in stable and convective atmospheric conditions

abstract: Wind tunnel experiments were conducted to study the impact of atmospheric stratification on flow and dispersion within and over a regular array of rectangular buildings. Three stable and two convective incoming boundary layers were tested with a Richardson number ranging from $-$1.5 to 0.29. Dispersion measurements were carried using a fast response flame ionisation detector. The results show that the stratification effect on the plume width is significantly lower than the effect on the vertical profiles. Stable stratification did not affect the plume central axis inside the canopy, but in the unstable case the axis appeared to deviate from the neutral case direction. Above the canopy both stratification types caused an increase in the plume deflection angle compared to the neutral case. Measured concentrations in stable stratification were up to two times larger in the canopy compared to the neutral case, while in convective conditions they were to three times smaller.The proportionality between the vertical turbulent fluxes and the vertical mean concentration gradient was also confirmed in the stratified cases. The high-quality experimental data produced during this work may help developing new mathematical models and parametrisation for non-neutral stratified conditions, as well as validating existing and future numerical simulations.


title: Slices of Parameter Space for Meromorphic Maps with Two Asymptotic Values

abstract: This paper is part of a program to understand the parameter spaces of dynamical systems generated by meromorphic functions with finitely many singular values. We give a full description of the parameter space for a specific family based on the exponential function that has precisely two finite asymptotic values and one attracting fixed point. It represents a step beyond the previous work in [GK] on degree 2 rational functions with analogous constraints: two critical values and an attracting fixed point. What is interesting and promising for pushing the general program even further, is that, despite the presence of the essential singularity, our new functions exhibit a dynamic structure as similar as one could hope to the rational case, and that the philosophy of the techniques used in the rational case could be adapted.


title: Pearson Distance is not a Distance

abstract: The Pearson distance between a pair of random variables $X,Y$ with correlation $\rho_{xy}$, namely, 1-$\rho_{xy}$, has gained widespread use, particularly for clustering, in areas such as gene expression analysis, brain imaging and cyber security. In all these applications it is implicitly assumed/required that the distance measures be metrics, thus satisfying the triangle inequality. We show however, that Pearson distance is not a metric. We go on to show that this can be repaired by recalling the result, (well known in other literature) that $\sqrt{1-\rho_{xy}}$ is a metric. We similarly show that a related measure of interest, $1-|\rho_{xy}|$, which is invariant to the sign of $\rho_{xy}$, is not a metric but that $\sqrt{1-\rho_{xy}^2}$ is. We also give generalizations of these results.


title: Signatures of noncommutativity in bar detectors of gravitational waves

abstract: The comparison between the noncommutative length scale $\sqrt{\theta}$ and the length variation $\delta L=h L$, detected in the GW detectors indicate that there is a strong possibility to detect the noncommutative structure of space in the GW detector set up. We therefore explore how the response of a bar detector gets affected due to the presence of noncommutative structure of space keeping terms upto second order in the gravitational wave perturbation ($h$) in the Hamiltonian. Interestingly, the second order term in $h$ shows a transition between the ground state and one of the perturbed second excited states that was absent when the calculation was restricted only to first order in $h$.


title: Gaia Parallax of Milky Way Globular Clusters -- A Solution of Mixture Model

abstract: We have established a mixture model approach to derive the parallax of the Milky Way globular clusters. It avoids the problem of cluster membership determination and provides a completely independent astrometrical solution by purely using the parallax data. This method is validated with simulated clusters of \cite{2017MNRAS.467..412P}. We have resolved 120 real globular clusters by the mixture model using parallaxes of the second data release of \gaia. They construct the largest direct parallax sample up to now. In comparing with other direct parallax results based on cluster members, including 75 clusters of \cite{2018A&A...616A..12G}, our method presents its accuracy, especially for some particular clusters. A systematic offset of $-27.6\pm1.7$ $\mu$as, together with a scatter of $22.8\pm1.3$ $\mu$as is found in comparing with other indirect parallax measurements. They are consistent with the global value and the variation of the zero-point of current \gaia parallaxes. Distances of several specific nearby globular clusters are discussed while the closest ones can reach high precisions, even taking the systematic error into account.


title: A case study of bilayered spin-$1/2$ square lattice compound [VO(HCOO)$_2\cdot$(H$_2$O)]

abstract: We present the synthesis and a detail investigation of structural and magnetic properties of polycrystalline [VO(HCOO)$_2\cdot$(H$_2$O)] by means of x-ray diffraction, magnetic susceptibility, high-field magnetization, heat capacity, and electron spin resonance measurements. It crystallizes in a orthorhombic structure with space group $Pcca$. It features distorted VO$_6$ octahedra connected via HCOO linker (formate anions) forming a two-dimensional square lattice network with a bilayered structure. Analysis of magnetic susceptibility, high field magnetization, and heat capacity data in terms of the frustrated square lattice model unambiguously establish quasi-two-dimensional nature of the compound with nearest neighbour interaction $J_1/k_{\rm B} \simeq 11.7$~K and next-nearest-neighbour interaction $J_2/k_{\rm B} \simeq 0.02$~K. It undergoes a N\'eel antiferromagnetic ordering at $T_{\rm N} \simeq 1.1$~K. The ratio $\theta_{\rm CW}/T_{\rm N} \simeq 10.9$ reflects excellent two-dimensionality of the spin-lattice in the compound. A strong in-plane anisotropy is inferred from the linear increase of $T_{\rm N}$ with magnetic field, consistent with the structural data.


title: Multi-waveband detection of quasi-periodic pulsations in a stellar flare on EK Draconis observed by XMM-Newton

abstract: Context. Quasi-periodic pulsations (QPPs) are time variations in the energy emission during a flare that are observed on both the Sun and other stars and thus have the potential to link the physics of solar and stellar flares. Aims. To characterise the QPPs detected in an X-ray flare on the solar analogue, EK Draconis, which was observed by XMM-Newton. Methods. We use wavelet and autocorrelation techniques to identify the QPPs in a detrended version of the flare. We also fit a model to the flare based on an exponential decay combined with a decaying sinusoid. The flare is examined in multiple energy bands. Results. A statistically significant QPP is observed in the X-ray energy band of 0.2-12.0 keV with a periodicity of 76+/-2 min. When this energy band is split, a statistically significant QPP is observed in the low-energy band (0.2-1.0 keV) with a periodicity of 73+/-2 min and in the high-energy band (1.0-12.0 keV) with a periodicity of 82+/-2 min. When fitting a model to the time series the phases of the signals are also found to be significantly different in the two energy bands (with a difference of 1.8+/-0.2 rad) and the high-energy band is found to lead the low-energy band. Furthermore, the first peak in the cross-correlation between the detrended residuals of the low- and high-energy bands is offset from zero by more than 3{\sigma} (4.1+/-1.3 min). Both energy bands produce statistically significant regions in the wavelet spectrum, whose periods are consistent with those listed above. However, the peaks are broad in both the wavelet and global power spectra, with the wavelet showing evidence for a drift in period with time, and the difference in period obtained is not significant. etc...


title: Slicing conditions for axisymmetric gravitational collapse of Brill waves

abstract: In numerical relativity, spacetimes involving compact strongly gravitating objects are constructed as numerical solutions of Einstein's equations. Success of such a process strongly depends on the availability of appropriate coordinates, which are typically constructed dynamically. A very robust coordinate choice is a so-called moving puncture gauge, commonly used for numerical simulations of black hole spacetimes. Nevertheless it is known to fail for evolving near-critical Brill wave data. We construct a new `quasi-maximal' slicing condition and demonstrate that it exhibits better behavior for such data. This condition is based on the 1+log slicing with an additional source term derived from maximal slicing. It is relatively simple to implement in existing moving puncture codes and computationally inexpensive. We also illustrate the properties of constructed spacetimes based on gauge-independent quantities in compactified spacetime diagrams. These invariants are also used to show how created black holes settle down to a Schwarzschild black hole.


title: PHY Research Is Sick but Curable: An Empirical Analysis

abstract: The controversy and argument on the usefulness of the physical layer (PHY) academic research for wireless communications are long-standing since the cellular communication paradigm gets to its maturity. In particular, researchers suspect that the performance improvement in cellular communications is primarily attributable to the increases in telecommunication infrastructure and radio spectrum instead of the PHY academic research, whereas concrete evidence is lacking. To respond to this controversy from an objective perspective, we employ econometric approaches to quantify the contributions of the PHY academic research and other performance determinants. Through empirical analysis and the quantitative evidence obtained, albeit preliminary, we shed light on the following issues: 1) what determines the cross-national differences in cellular network performance; 2) to what extent the PHY academic research and other factors affect cellular network performance; 3) what suggestions we can obtain from the data analysis for the stakeholders of the PHY research. To the best of our knowledge, this article is the first `empirical telecommunication research,' and the first effort to involve econometric methodologies to evaluate the usefulness of the PHY academic research.


title: A proposal for nonabelian (0,2) mirrors

abstract: In this paper we give a proposal for mirrors to (0,2) supersymmetric gauged linear sigma models (GLSMs), for those (0,2) GLSMs which are deformations of (2,2) GLSMs. Specifically, we propose a construction of (0,2) mirrors for (0,2) GLSMs with E terms that are linear and diagonal, reducing to both the Hori-Vafa prescription as well as a recent (2,2) nonabelian mirrors proposal on the (2,2) locus. For the special case of abelian (0,2) GLSMs, two of the authors have previously proposed a systematic construction, which is both simplified and generalized by the proposal here.


title: Needles in Haystacks: On Classifying Tiny Objects in Large Images

abstract: In some important computer vision domains, such as medical or hyperspectral imaging, we care about the classification of tiny objects in large images. However, most Convolutional Neural Networks (CNNs) for image classification were developed using biased datasets that contain large objects, in mostly central image positions. To assess whether classical CNN architectures work well for tiny object classification we build a comprehensive testbed containing two datasets: one derived from MNIST digits and one from histopathology images. This testbed allows controlled experiments to stress-test CNN architectures with a broad spectrum of signal-to-noise ratios. Our observations indicate that: (1) There exists a limit to signal-to-noise below which CNNs fail to generalize and that this limit is affected by dataset size - more data leading to better performances; however, the amount of training data required for the model to generalize scales rapidly with the inverse of the object-to-image ratio (2) in general, higher capacity models exhibit better generalization; (3) when knowing the approximate object sizes, adapting receptive field is beneficial; and (4) for very small signal-to-noise ratio the choice of global pooling operation affects optimization, whereas for relatively large signal-to-noise values, all tested global pooling operations exhibit similar performance.


title: Some hemivariational inequalities in the Euclidean space

abstract: The purpose of this paper is to study the existence of weak solutions for some classes of hemivariational problems in the Euclidean space $\mathbb{R}^d$ ($d\geq 3$). These hemivariational inequalities have a variational structure and, thanks to this, we are able to find a non-trivial weak solution for them by using variational methods and a non-smooth version of the Palais principle of symmetric criticality for locally Lipschitz continuous functionals, due to Krawcewicz and Marzantowicz. The main tools in our approach are based on appropriate theoretical arguments on suitable subgroups of the orthogonal group $O(d)$ and their actions on the Sobolev space $H^1(\mathbb{R}^d)$. Moreover, under an additional hypotheses on the dimension $d$ and in the presence of symmetry on the nonlinear datum, the existence of multiple pairs of sign-changing solutions with different symmetries structure has been proved. In connection to classical Schr\"{o}dinger equations a concrete and meaningful example of an application is presented.


title: Few-shot Text Classification with Distributional Signatures

abstract: In this paper, we explore meta-learning for few-shot text classification. Meta-learning has shown strong performance in computer vision, where low-level patterns are transferable across learning tasks. However, directly applying this approach to text is challenging--lexical features highly informative for one task may be insignificant for another. Thus, rather than learning solely from words, our model also leverages their distributional signatures, which encode pertinent word occurrence patterns. Our model is trained within a meta-learning framework to map these signatures into attention scores, which are then used to weight the lexical representations of words. We demonstrate that our model consistently outperforms prototypical networks learned on lexical knowledge (Snell et al., 2017) in both few-shot text classification and relation classification by a significant margin across six benchmark datasets (20.0% on average in 1-shot classification).


title: Performing Deep Recurrent Double Q-Learning for Atari Games

abstract: Currently, many applications in Machine Learning are based on define new models to extract more information about data, In this case Deep Reinforcement Learning with the most common application in video games like Atari, Mario, and others causes an impact in how to computers can learning by himself with only information called rewards obtained from any action. There is a lot of algorithms modeled and implemented based on Deep Recurrent Q-Learning proposed by DeepMind used in AlphaZero and Go. In this document, We proposed Deep Recurrent Double Q-Learning that is an implementation of Deep Reinforcement Learning using Double Q-Learning algorithms and Recurrent Networks like LSTM and DRQN.


title: On the eigenvalues of the Robin Laplacian with a complex parameter

abstract: We study the spectrum of the Robin Laplacian with a complex Robin parameter $\alpha$ on a bounded Lipschitz domain $\Omega$. We start by establishing a number of properties of the corresponding operator, such as generation properties, local analytic dependence of the eigenvalues and eigenspaces on $\alpha \in \mathbb C$, and basis properties of the eigenfunctions. Our focus, however, is on bounds and asymptotics for the eigenvalues as functions of $\alpha$: we start by providing estimates on the numerical range of the associated operator, which lead to new eigenvalue bounds even in the case $\alpha \in \mathbb R$. For the asymptotics of the eigenvalues as $\alpha \to \infty$ in $\mathbb C$, in place of the min-max characterisation of the eigenvalues and Dirichlet-Neumann bracketing techniques commonly used in the real case, we exploit the duality between the eigenvalues of the Robin Laplacian and the eigenvalues of the Dirichlet-to-Neumann map. We use this to show that every Robin eigenvalue either diverges to $\infty$ in $\mathbb C$ or converges to a point in the spectrum of the Dirichlet Laplacian, and also to give a comprehensive treatment of the special cases where $\Omega$ is an interval, a hyperrectangle or a ball. This leads to the conjecture that on a general smooth domain in dimension $d\geq 2$ all eigenvalues converge to the Dirichlet spectrum if ${\rm Re}\, \alpha$ remains bounded from below as $\alpha \to \infty$, while if ${\rm Re}\, \alpha \to -\infty$, then there is a family of divergent eigenvalue curves, each of which behaves asymptotically like $-\alpha^2$.


title: Thermomagnetic Ettingshausen-Nernst effect in tachocline, magnetic reconnection phenomenon in lower layers, axion mechanism of solar luminosity variations, coronal heating problem solution and mechanism of ADM variations around BH

abstract: It is shown that the holographic principle of quantum gravity (in the hologram of the Universe, and therefore in our Galaxy, and of course on the Sun!), in which the conflict between the theory of gravitation and quantum mechanics disappears, gives rise to the Babcock-Leighton holographic mechanism. Unlike the solar dynamo models, it generates a strong toroidal magnetic field by means of the thermomagnetic Ettingshausen-Nernst (EN) effect in the tachocline. Hence, it can be shown that with the help of the thermomagnetic EN effect, a simple estimate of the magnetic pressure of an ideal gas in the tachocline of e.g. the Sun can indirectly prove that by using the holographic principle of quantum gravity, the repulsive toroidal magnetic field of the tachocline ($B_{tacho}^{Sun} = 4.1 \cdot 10^7 ~G = - B_{core}^{Sun}$) precisely "neutralizes" the magnetic field in the Sun core, since the projections of the magnetic fields in the tachocline and the core have equal values but opposite directions. The basic problem is a generalized problem of the antidynamo model of magnetic flux tubes (MFTs), where the nature of both holographic effects (the thermomagnetic EN~effect and Babcock-Leighton holographic mechanism), including magnetic cycles, manifests itself in the modulation of asymmetric dark matter (ADM) and, consequently, the solar axion in the Sun interior.


title: A Shift Selection Strategy for Parallel Shift-Invert Spectrum Slicing in Symmetric Self-Consistent Eigenvalue Computation

abstract: The central importance of large scale eigenvalue problems in scientific computation necessitates the development of massively parallel algorithms for their solution. Recent advances in dense numerical linear algebra have enabled the routine treatment of eigenvalue problems with dimensions on the order of hundreds of thousands on the world's largest supercomputers. In cases where dense treatments are not feasible, Krylov subspace methods offer an attractive alternative due to the fact that they do not require storage of the problem matrices. However, demonstration of scalability of either of these classes of eigenvalue algorithms on computing architectures capable of expressing massive parallelism is non-trivial due to communication requirements and serial bottlenecks, respectively. In this work, we introduce the SISLICE method: a parallel shift-invert algorithm for the solution of the symmetric self-consistent field (SCF) eigenvalue problem. The SISLICE method drastically reduces the communication requirement of current parallel shift-invert eigenvalue algorithms through various shift selection and migration techniques based on density of states estimation and k-means clustering, respectively. This work demonstrates the robustness and parallel performance of the SISLICE method on a representative set of SCF eigenvalue problems and outlines research directions which will be explored in future work.


title: A finitary structure theorem for vertex-transitive graphs of polynomial growth

abstract: We prove a quantitative, finitary version of Trofimov's result that a connected, locally finite vertex-transitive graph G of polynomial growth admits a quotient with finite fibres on which the action of Aut(G) is virtually nilpotent with finite vertex stabilisers. We also present some applications. We show that a finite, connected vertex-transitive graph G of large diameter admits a quotient with fibres of small diameter on which the action of Aut(G) is virtually abelian with vertex stabilisers of bounded size. We also show that G has moderate growth in the sense of Diaconis and Saloff-Coste, which is known to imply that the mixing and relaxation times of the lazy random walk on G are quadratic in the diameter. These results extend results of Breuillard and the second author for finite Cayley graphs of large diameter. Finally, given a connected, locally finite vertex-transitive graph G exhibiting polynomial growth at a single, sufficiently large scale, we describe its growth at subsequent scales, extending a result of Tao and an earlier result of our own for Cayley graphs. In forthcoming work we will give further applications.


title: Neutrino, Electroweak and Nuclear Physics from COHERENT Elastic Neutrino-Nucleus Scattering with a New Quenching Factor

abstract: We present an updated analysis of the coherent neutrino-nucleus elastic scattering data of the COHERENT experiment taking into account the new quenching factor published recently in Phys. Rev. D100, 033003 (2019). Through a fit of the COHERENT time-integrated energy spectrum, we show that the new quenching factor leads to a better determination of the average rms radius of the neutron distributions of $^{133}\text{Cs}$ and $^{127}\text{I}$, while in combination with the atomic parity violation (APV) experimental results it allows to determine a data-driven APV measurement of the low-energy weak mixing angle in very good agreement with the Standard Model prediction. We also find a $3.7\sigma$ evidence of the suppression of coherence due to the nuclear structure. Neutrino properties are better constrained by considering the COHERENT time-dependent spectral data, that allow us to improve the bounds on the neutrino charge radii and magnetic moments. We also present for the first time constraints on the neutrino charges obtained with coherent neutrino-nucleus elastic scattering data. In particular, we obtain the first laboratory constraints on the diagonal charge of $\nu_{\mu}$ and the $\nu_{\mu}$-$\nu_{\tau}$ transition charge.


title: A Projection-based Augmented Reality for Elderly People with Dementia

abstract: As aging societies grow, researchers are actively studying care systems concerning the life and diseases of the elderly. Among these diseases, dementia makes it difficult to maintain daily life due to the degradation of cognitive functioning, memory, and reasoning, as well as the ability to perform actions. Moreover, dementia does not have a perfect cure, though therapy and care can slow its onset and provide patients with physical and mental support. In this paper, we developed a projection-based augmented reality system robot that can cover 360 degrees of space. We also propose an application that supports continuous monitoring of dementia patients to address the difficulties they face in daily life. The system is also designed to provide therapy applications, such as entertainment and spatial art, to provide mental care aids for the patients.


title: Robust Principal Component Analysis for Background Estimation of Particle Image Velocimetry Data

abstract: Particle Image Velocimetry (PIV) data processing procedures are adversely affected by light reflections and backgrounds as well as defects in the models and sticky particles that occlude the inner walls of the boundaries. In this paper, a novel approach is proposed for decomposition of the PIV data into background/foreground components, greatly reducing the effects of such artifacts. This is achieved by utilizing Robust Principal Component Analysis (RPCA) applied to the data matrix, generated by aggregating the vectorized PIV frames. It is assumed that the data matrix can be decomposed into two statistically different components, a low-rank component depicting the still background and a sparse component representing the moving particles within the imaged geometry. Formulating the assumptions as an optimization problem, Augmented Lagrange Multiplier (ALM) method is used for decomposing the data matrix into the low-rank and sparse components. Experiments and comparisons with the state-of-the-art using several PIV image sequences reveal the superiority of the proposed approach for background removal of PIV data.


title: Statistical analysis on X-ray flares from the nucleus and HST-1 knot in the M87 jet

abstract: The statistical properties of X-ray flares from two separate locations (nucleus and HST-1) in the M87 jet are investigated to reveal the physical origin of the flares. We analyse the archival \textit{Chandra} data for M87, and identify 14 flares in the nucleus and 9 flares in HST-1. The peak intensity ($I_{\rm{P}}$) and the flaring duration time ($T_{\rm{fl}}$) for each flare are obtained. It is found that the distributions of both $I_{\rm{P}}$ and $T_{\rm{fl}}$ for the nucleus obey a power-law form with a similar index. A similar result is also obtained for HST-1, and no significant inconsistency between the nucleus and HST-1 is found for the indices. Similar to solar X-ray flares, the power-law distributions of the flare event parameters can be well explained by a self-organized criticality (SOC) system, which are triggered by magnetic reconnection. Our results suggest that the flares from nucleus and HST-1 are possibly triggered by magnetic reconnection process. The consistent indices for the distributions of $I_{\rm{P}}$ and $T_{\rm{fl}}$ in the CORE and HST-1 indicate that the dimensions of the energy dissipation of the magnetic reconnection is identical in the two regions. A strong correlation between the flares in the two regions also suggests a similar physical origin for the flares.


title: ZeroER: Entity Resolution using Zero Labeled Examples

abstract: Entity resolution (ER) refers to the problem of matching records in one or more relations that refer to the same real-world entity. While supervised machine learning (ML) approaches achieve the state-of-the-art results, they require a large amount of labeled examples that are expensive to obtain and often times infeasible. We investigate an important problem that vexes practitioners: is it possible to design an effective algorithm for ER that requires Zero labeled examples, yet can achieve performance comparable to supervised approaches? In this paper, we answer in the affirmative through our proposed approach dubbed ZeroER. Our approach is based on a simple observation -- the similarity vectors for matches should look different from that of unmatches. Operationalizing this insight requires a number of technical innovations. First, we propose a simple yet powerful generative model based on Gaussian Mixture Models for learning the match and unmatch distributions. Second, we propose an adaptive regularization technique customized for ER that ameliorates the issue of feature overfitting. Finally, we incorporate the transitivity property into the generative model in a novel way resulting in improved accuracy. On five benchmark ER datasets, we show that ZeroER greatly outperforms existing unsupervised approaches and achieves comparable performance to supervised approaches.


title: Cosmological Inference using Gravitational Wave Standard Sirens: A Mock Data Challenge

abstract: The observation of binary neutron star merger GW170817, along with its optical counterpart, provided the first constraint on the Hubble constant $H_0$ using gravitational wave standard sirens. When no counterpart is identified, a galaxy catalog can be used to provide the necessary redshift information. However, the true host might not be contained in a catalog which is not complete out to the limit of gravitational-wave detectability. These electromagnetic and gravitational-wave selection effects must be accounted for. We describe and implement a method to estimate $H_0$ using both the counterpart and the galaxy catalog standard siren methods. We perform a series of mock data analyses using binary neutron star mergers to confirm our ability to recover an unbiased estimate of $H_0$. Our simulations used a simplified universe with no redshift uncertainties or galaxy clustering, but with different magnitude-limited catalogs and assumed host galaxy properties, to test our treatment of both selection effects. We explore how the incompleteness of catalogs affects the final measurement of $H_0$, as well as the effect of weighting each galaxy's likelihood of being a host by its luminosity. In our most realistic simulation, where the simulated catalog is about three times denser than the density of galaxies in the local universe, we find that a 4.4\% measurement precision can be reached using galaxy catalogs with 50\% completeness and $\sim 250$ binary neutron star detections with sensitivity similar to that of Advanced LIGO's second observing run.


title: Minimum Coprime Labelings of Generalized Petersen and Prism Graphs

abstract: A coprime labeling of a graph of order $n$ is an assignment of distinct positive integer labels in which adjacent vertices have relatively prime labels. Restricting labels to only the set $1$ to $n$ results in a prime labeling. In this paper, we consider families of graphs in which a prime labeling cannot exist with the goal being to minimize the largest value of the labeling set, resulting in a minimum coprime labeling. In particular, prism graphs, generalized Petersen graphs with $k=2$, and stacked prism graphs are investigated for minimum coprime labelings.


title: Recover and Identify: A Generative Dual Model for Cross-Resolution Person Re-Identification

abstract: Person re-identification (re-ID) aims at matching images of the same identity across camera views. Due to varying distances between cameras and persons of interest, resolution mismatch can be expected, which would degrade person re-ID performance in real-world scenarios. To overcome this problem, we propose a novel generative adversarial network to address cross-resolution person re-ID, allowing query images with varying resolutions. By advancing adversarial learning techniques, our proposed model learns resolution-invariant image representations while being able to recover the missing details in low-resolution input images. The resulting features can be jointly applied for improving person re-ID performance due to preserving resolution invariance and recovering re-ID oriented discriminative details. Our experiments on five benchmark datasets confirm the effectiveness of our approach and its superiority over the state-of-the-art methods, especially when the input resolutions are unseen during training.


title: Braidoids

abstract: Braidoids generalize the classical braids and form a counterpart theory to the theory of planar knotoids, just as the theory of braids does for the theory of knots. In this paper, we introduce basic notions of braidoids, a closure operation for braidoids, we prove an analogue of the Alexander theorem, that is, an algorithm that turns a knotoid into a braidoid, and we formulate and prove a geometric analogue of the Markov theorem for braidoids using the $L$-moves.


title: Angular momentum of the electron: One-loop studies

abstract: We combine bare perturbation theory with the imaginary time evolution technique to study one-loop radiative corrections to various components of angular momentum of the electron. Our investigations are based on the canonical decomposition of angular momentum, where spin and orbital components, associated with fermionic and electromagnetic degrees of freedom, are individually approached. We use for this purpose quantum electrodynamics in the general covariant gauge and develop a formalism, based on the repeated use of the Sochocki-Plemelj formula, for proper enforcement of the imaginary time limit. It is then shown that careful implementation of imaginary time evolutions is crucial for getting a correct result for total angular momentum of the electron in the bare perturbative expansion. We also analyze applicability of the Pauli-Villars regularization to our problem, developing a variant of this technique based on modifications of studied observables by subtraction of their ghost operator counterparts. It is then shown that such an approach leads to the consistent regularization of all angular momenta that we compute.


title: Itinerant ferromagnetism and intrinsic anomalous Hall effect in amorphous iron-germanium

abstract: The amorphous iron-germanium system ($a$-Fe$_x$Ge$_{1-x}$) lacks long-range structural order and hence lacks a meaningful Brillouin zone. The magnetization of \aFeGe is well explained by the Stoner model for Fe concentrations $x$ above the onset of magnetic order around $x=0.4$, indicating that the local order of the amorphous structure preserves the spin-split density of states of the Fe-$3d$ states sufficiently to polarize the electronic structure despite $\mathbf{k}$ being a bad quantum number. Measurements reveal an enhanced anomalous Hall resistivity $\rho_{xy}^{\mathrm{AH}}$ relative to crystalline FeGe; this $\rho_{xy}^{\mathrm{AH}}$ is compared to density functional theory calculations of the anomalous Hall conductivity to resolve its underlying mechanisms. The intrinsic mechanism, typically understood as the Berry curvature integrated over occupied $\mathbf{k}$-states but shown here to be equivalent to the density of curvature integrated over occupied energies in aperiodic materials, dominates the anomalous Hall conductivity of $a$-Fe$_x$Ge$_{1-x}$ ($0.38 \leq x \leq 0.61$). The density of curvature is the sum of spin-orbit correlations of local orbital states and can hence be calculated with no reference to $\mathbf{k}$-space. This result and the accompanying Stoner-like model for the intrinsic anomalous Hall conductivity establish a unified understanding of the underlying physics of the anomalous Hall effect in both crystalline and disordered systems.


title: Astrometric planet search around southern ultracool dwarfs IV. Relative motion of the FORS2/VLT CCD chips

abstract: We present an investigation of the stability of the two chips in the FORS2 camera CCD mosaic on the basis of astrometric observations of stars in 20 sky fields, some of which were monitored for four to seven years. We detected a smooth relative shear motion of the chips along their dividing line that is well approximated by a cubic function of time with an amplitude that reaches ~0.3 pixels (px) or ~38 mas over seven years. In a single case, we detected a step change of ~0.06 px that occurred within four days. In the orthogonal direction that corresponds to the separation between the chips, the motion is a factor of 5-10 smaller. This chip instability in the camera significantly reduces the astrometric precision when the reduction uses reference stars located in both chips, and the effect is not accounted for explicitly. We found that the instability introduces a bias in stellar positions with an amplitude that increases with the observation time span. When our reduction methods and FORS2 images are used, it affects stellar positions like an excess random noise with an RMS of ~0.5 mas for a time span of three to seven years when left uncorrected. We demonstrate that an additional calibration step can adequately mitigate this and restore an astrometric accuracy of 0.12 mas, which is essential to achieve the goals of our planet-search program. These results indicate that similar instabilities could critically affect the astrometric performance of other large ground-based telescopes and extremely large telescopes that are equipped with large-format multi-chip detectors if no precautions are taken.


title: Generalized group-based epidemic model for spreading processes on networks: GgroupEM

abstract: We develop a generalized group-based epidemic model (GgroupEM) framework for any compartmental epidemic model (for example; susceptible-infected-susceptible, susceptible-infected-recovered, susceptible-exposed-infected-recovered). Here, a group consists of a collection of individual nodes. This model can be used to understand the important dynamic characteristics of a stochastic epidemic spreading over very large complex networks, being informative about the state of groups. Aggregating nodes by groups, the state space becomes smaller than the individual-based approach at the cost of aggregation error, which is strongly bounded by the isoperimetric inequality. We also develop a mean-field approximation of this framework to further reduce the state-space size. Finally, we extend the GgroupEM to multilayer networks. Since the group-based framework is computationally less expensive and faster than an individual-based framework, then this framework is useful when the simulation time is important.


title: Lower bounds in the polynomial Szemer\'edi theorem

abstract: We construct large subsets of the first $N$ positive integers which avoid certain arithmetic configurations. In particular, we construct a set of order $N^{0.7685}$ lacking the configuration $\{x,x+y,x+y^2\},$ surpassing the $N^{3/4}$ limit of Ruzsa's construction for sets lacking a square difference. We also extend Ruzsa's construction to sets lacking polynomial differences for a wide class of univariate polynomials. Finally, we turn to multivariate differences, constructing a set of order $N^{1/2}$ lacking a difference equal to a sum of two squares. This is in contrast to the analogous problem of sets lacking a difference equal to a prime minus one, where the current record is of order $N^{o(1)}.$


title: Low-rank approximation in the Frobenius norm by column and row subset selection

abstract: A CUR approximation of a matrix $A$ is a particular type of low-rank approximation $A \approx C U R$, where $C$ and $R$ consist of columns and rows of $A$, respectively. One way to obtain such an approximation is to apply column subset selection to $A$ and $A^T$. In this work, we describe a numerically robust and much faster variant of the column subset selection algorithm proposed by Deshpande and Rademacher, which guarantees an error close to the best approximation error in the Frobenius norm. For cross approximation, in which $U$ is required to be the inverse of a submatrix of $A$ described by the intersection of $C$ and $R$, we obtain a new algorithm with an error bound that stays within a factor $k + 1$ of the best rank-$k$ approximation error in the Frobenius norm. To the best of our knowledge, this is the first deterministic polynomial-time algorithm for which this factor is bounded by a polynomial in $k$. Our derivation and analysis of the algorithm is based on derandomizing a recent existence result by Zamarashkin and Osinsky. To illustrate the versatility of our new column subset selection algorithm, an extension to low multilinear rank approximations of tensors is provided as well.


title: A gravitational-wave measurement of the Hubble constant following the second observing run of Advanced LIGO and Virgo

abstract: This paper presents the gravitational-wave measurement of the Hubble constant ($H_0$) using the detections from the first and second observing runs of the Advanced LIGO and Virgo detector network. The presence of the transient electromagnetic counterpart of the binary neutron star GW170817 led to the first standard-siren measurement of $H_0$. Here we additionally use binary black hole detections in conjunction with galaxy catalogs and report a joint measurement. Our updated measurement is $H_0 = 68.7^{+17.0}_{-7.8}$ km/s/Mpc (68.3\% of the highest density posterior interval with a flat-in-log prior) which is an improvement by a factor of 1.04 (about 4\%) over the GW170817-only value of $68.7^{+17.5}_{-8.3}$ km/s/Mpc. A significant additional contribution currently comes from GW170814, a loud and well-localized detection from a part of the sky thoroughly covered by the Dark Energy Survey. With numerous detections anticipated over the upcoming years, an exhaustive understanding of other systematic effects are also going to become increasingly important. These results establish the path to cosmology using gravitational-wave observations with and without transient electromagnetic counterparts.


title: Electromagnetic angular momentum of the electron: One-loop studies

abstract: We study angular momentum of the electron stored in its electric and magnetic fields. We use for this purpose quantum electrodynamics in the covariant gauge. We show that a finite one-loop result for such angular momentum can be obtained without invoking any renormalization procedure. We compare it to the classical estimation relying on a short-distance cutoff.


title: Adversarial shape perturbations on 3D point clouds

abstract: The importance of training robust neural network grows as 3D data is increasingly utilized in deep learning for vision tasks in robotics, drone control, and autonomous driving. One commonly used 3D data type is 3D point clouds, which describe shape information. We examine the problem of creating robust models from the perspective of the attacker, which is necessary in understanding how 3D neural networks can be exploited. We explore two categories of attacks: distributional attacks that involve imperceptible perturbations to the distribution of points, and shape attacks that involve deforming the shape represented by a point cloud. We explore three possible shape attacks for attacking 3D point cloud classification and show that some of them are able to be effective even against preprocessing steps, like the previously proposed point-removal defenses.


title: Improvements on "Secure multi-party quantum summation based on quantum Fourier transform"

abstract: Recently, a quantum multi-party summation protocol based on the quantum Fourier transform has been proposed [Quantum Inf Process 17: 129, 2018]. The protocol claims to be secure against both outside and participant attacks. However, a closer look reveals that the player in charge of generating the required multi-partite entangled states can launch two kinds of attacks to learn about other parties' private integer strings without being caught. In this paper, we present these attacks, and propose countermeasures to make the protocol secure again. The improved protocol not only can resist these attacks but also remove the need for the quantum Fourier transform and encoding quantum operations by participants.


title: Graviton Propagator in a 2-Parameter Family of de Sitter Breaking Gauges

abstract: We formulate the graviton propagator on de Sitter background in a 2-parameter family of simple gauges which break de Sitter invariance. Explicit results are derived for the first order perturbations in each parameter. These results should be useful in computations to check for gauge dependence of graviton loop corrections.


title: On Convex Duality in Linear Inverse Problems

abstract: In this article we dwell into the class of so called ill posed Linear Inverse Problems (LIP) in machine learning, which has become almost a classic in recent times. The fundamental task in an LIP is to recover the entire signal / data from its relatively few random linear measurements. Such problems arise in variety of settings with applications ranging from medical image processing, recommender systems etc. We provide an exposition to the convex duality of the linear inverse problems, and obtain a novel and equivalent convex-concave min-max reformulation that gives rise to simple ascend-descent type algorithms to solve an LIP. Moreover, such a reformulation is crucial in developing methods to solve the dictionary learning problem with almost sure recovery constraints.


title: Unicoder-VL: A Universal Encoder for Vision and Language by Cross-modal Pre-training

abstract: We propose Unicoder-VL, a universal encoder that aims to learn joint representations of vision and language in a pre-training manner. Borrow ideas from cross-lingual pre-trained models, such as XLM and Unicoder, both visual and linguistic contents are fed into a multi-layer Transformer for the cross-modal pre-training, where three pre-trained tasks are employed, including Masked Language Modeling (MLM), Masked Object Classification (MOC) and Visual-linguistic Matching (VLM). The first two tasks learn context-aware representations for input tokens based on linguistic and visual contents jointly. The last task tries to predict whether an image and a text describe each other. After pretraining on large-scale image-caption pairs, we transfer Unicoder-VL to caption-based image-text retrieval and visual commonsense reasoning, with just one additional output layer. We achieve state-of-the-art or comparable results on both two tasks and show the powerful ability of the cross-modal pre-training.


title: A Practical and Consistent Parametrization of Dark Matter Self-Interactions

abstract: Self-interacting dark matter has been proposed to explain the apparent mass deficit in astrophysical small-scale halos, while observations from galaxy clusters suggest that the corresponding cross section depends on the velocity. Accounting for this is often believed to be highly model-dependent with studies mostly focusing on scenarios with light mediators. Based on the effective-range formalism, in this work we point out a model-independent approach which accurately approximates the velocity dependence of the self-interaction cross section with only two parameters. We illustrate how this parameterization can be simultaneously interpreted in various well-motivated scenarios, including self-interactions induced by Yukawa forces, Breit-Wigner resonances and bound states. We investigate the astrophysical implications and discuss how the approximation can be improved in certain special regimes where it works poorly.


title: Covert Computation in Self-Assembled Circuits

abstract: Traditionally, computation within self-assembly models is hard to conceal because the self-assembly process generates a crystalline assembly whose computational history is inherently part of the structure itself. With no way to remove information from the computation, this computational model offers a unique problem: how can computational input and computation be hidden while still computing and reporting the final output? Designing such systems is inherently motivated by privacy concerns in biomedical computing and applications in cryptography. In this paper we propose the problem of performing ``covert computation'' within tile self-assembly that seeks to design self-assembly systems that ``conceal'' both the input and computational history of performed computations. We achieve these results within the growth-only restricted abstract Tile Assembly Model (aTAM) with positive and negative interactions. We show that general-case covert computation is possible by implementing a set of basic covert logic gates capable of simulating any circuit (functionally complete). To further motivate the study of covert computation, we apply our new framework to resolve an outstanding complexity question; we use our covert circuitry to show that the unique assembly verification problem within the growth-only aTAM with negative interactions is coNP-complete.


title: A Survey on Computational Politics

abstract: Computational Politics is the study of computational methods to analyze and moderate users' behaviors related to political activities such as election campaign persuasion, political affiliation, and opinion mining. With the rapid development and ease of access to the Internet, Information Communication Technologies (ICT) have given rise to massive numbers of users joining online communities and the digitization of analogous data such as political debates. These communities and digitized data contain both explicit and latent information about users and their behaviors related to politics. For researchers, it is essential to utilize data from these sources to develop and design systems that not only provide solutions to computational politics but also help other businesses, such as marketers to increase users, participation and interactions. In this survey, we attempt to categorize main areas in computational politics and summarize the prominent studies in one place to better understand computational politics across different and multidimensional platforms. e.g., online social networks, online forums, and political debates. We then conclude this study by highlighting future research directions, opportunities, and challenges.


title: Optimal scheduling strategy for networked estimation with energy harvesting

abstract: Joint optimization of scheduling and estimation policies is considered for a system with two sensors and two non-collocated estimators. Each sensor produces an independent and identically distributed sequence of random variables, and each estimator forms estimates of the corresponding sequence with respect to the mean-squared error sense. The data generated by the sensors is transmitted to the corresponding estimators, over a bandwidth-constrained wireless network that can support a single packet per time slot. The access to the limited communication resources is determined by a scheduler who decides which sensor measurement to transmit based on both observations. The scheduler has an energy-harvesting battery of limited capacity, which couples the decision-making problem in time. Despite the overall lack of convexity of the team decision problem, it is shown that this system admits globally optimal scheduling and estimation strategies under the assumption that the distributions of the random variables at the sensors are symmetric and unimodal. Additionally, the optimal scheduling policy has a structure characterized by a threshold function that depends on the time index and energy level. A recursive algorithm for threshold computation is provided.


title: QCD Axion on Hilltop by a Phase Shift of $\pi$

abstract: We show that the initial misalignment angle of the QCD axion (or axion-like particles) can be set very close to $\pi$, if the QCD axion has a mixing with another heavy axion which induces the phase shift $\approx \pi$ after inflation. In the simplest case, the heavy axion plays the role of the inflaton, and we call such inflation as "$\pi$nflation." The basic idea was first proposed by Daido and the present authors in Ref. [1702.03284] in 2017 and more recently discussed in Ref. [1903.00462]. We show that the QCD axion with a decay constant $f_a \gtrsim 3 \times 10^9\,$ GeV can explain dark matter by the $\pi$nflation mechanism. A large fraction of the parameter region has an overlap with the projected sensitivity of ORGAN, MADMAX, TOORAD and IAXO. We also study implications for the effective neutrino species and isocurvature perturbations. The $\pi$nflation can provide an initial condition for the hilltop inflation in the axion landscape, and in a certain set-up, a chain of the hilltop inflation may take place.


title: Antimagic orientations of graphs with large maximum degree

abstract: Given a digraph $D$ with $m $ arcs, a bijection $\tau: A(D)\rightarrow \{1, 2, \ldots, m\}$ is an antimagic labeling of $D$ if no two vertices in $D$ have the same vertex-sum, where the vertex-sum of a vertex $u $ in $D$ under $\tau$ is the sum of labels of all arcs entering $u$ minus the sum of labels of all arcs leaving $u$. We say $(D, \tau)$ is an antimagic orientation of a graph $G$ if $D$ is an orientation of $G$ and $\tau$ is an antimagic labeling of $D$. Motivated by the conjecture of Hartsfield and Ringel from 1990 on antimagic labelings of graphs, Hefetz, M\"{u}tze, and Schwartz in 2010 initiated the study of antimagic orientations of graphs, and conjectured that every connected graph admits an antimagic orientation. This conjecture seems hard, and few related results are known. However, it has been verified to be true for regular graphs and biregular bipartite graphs. In this paper, we prove that every connected graph $G$ on $n\ge9$ vertices with maximum degree at least $n-5$ admits an antimagic orientation.


title: Studying Cosmic-ray Interactions in Giant Molecular Clouds with the HAWC Gamma-ray Observatory

abstract: The cosmic-ray flux in the Galaxy can be characterized by combining the knowledge of the distribution of gas in the Galaxy and the observation of gamma rays. We analyze the data from the HAWC Observatory to look for gamma rays in three galactic giant molecular clouds, that are outside the galactic plane ($|b|>5^{\circ}$). We can then test the paradigm that the measured local cosmic-ray flux is the same as the "sea" of Galactic cosmic rays. Due to its large field of view, and high duty cycle, HAWC is suitable to search for gamma rays from large structures in the TeV gamma-ray regime. We present here preliminary results from measurements of the Aquila Rift, Hercules and Taurus molecular clouds.


title: Configurational entropy and instability of tachyonic braneworld

abstract: We consider tachyonic braneworld with a bulk cosmological constant and investigate a configurational entropy of various magnitudes of scale factor. It is found that for a bulk negative/zero cosmological constant, the configurational entropy has a global minimum when the magnitude of scale factor reaches the critical value. This result seems to have intriguing implications such that an accelerated rate of the universe and cosmological inflation rate for radiation/matter domination are able to be determined by such critical value. We also find that the configurational entropy almost monotonically decreases for a bulk positive cosmological constant as the magnitude of scale factor grows up. We find an exact solution of tachyonic braneworld in a bulk de Sitter space. It is shown that such system under scalar perturbations is stable for some constraint relation. Furthermore, we also find that tachyonic braneworld model with a bulk negative/zero cosmological constant is always stable under scalar perturbations.


title: An Exploratory Analysis of the Latent Structure of Process Data via Action Sequence Autoencoder

abstract: Computer simulations have become a popular tool of assessing complex skills such as problem-solving skills. Log files of computer-based items record the entire human-computer interactive processes for each respondent. The response processes are very diverse, noisy, and of nonstandard formats. Few generic methods have been developed for exploiting the information contained in process data. In this article, we propose a method to extract latent variables from process data. The method utilizes a sequence-to-sequence autoencoder to compress response processes into standard numerical vectors. It does not require prior knowledge of the specific items and human-computers interaction patterns. The proposed method is applied to both simulated and real process data to demonstrate that the resulting latent variables extract useful information from the response processes.


title: Number-Theoretic Characterizations of Some Restricted Clifford+T Circuits

abstract: Kliuchnikov, Maslov, and Mosca proved in 2012 that a $2\times 2$ unitary matrix $V$ can be exactly represented by a single-qubit Clifford+$T$ circuit if and only if the entries of $V$ belong to the ring $\mathbb{Z}[1/\sqrt{2},i]$. Later that year, Giles and Selinger showed that the same restriction applies to matrices that can be exactly represented by a multi-qubit Clifford+$T$ circuit. These number-theoretic characterizations shed new light upon the structure of Clifford+$T$ circuits and led to remarkable developments in the field of quantum compiling. In the present paper, we provide number-theoretic characterizations for certain restricted Clifford+$T$ circuits by considering unitary matrices over subrings of $\mathbb{Z}[1/\sqrt{2},i]$. We focus on the subrings $\mathbb{Z}[1/2]$, $\mathbb{Z}[1/\sqrt{2}]$, $\mathbb{Z}[1/i\sqrt{2}]$, and $\mathbb{Z}[1/2,i]$, and we prove that unitary matrices with entries in these rings correspond to circuits over well-known universal gate sets. In each case, the desired gate set is obtained by extending the set of classical reversible gates $\{X, CX, CCX\}$ with an analogue of the Hadamard gate and an optional phase gate.


title: NUQSGD: Provably Communication-efficient Data-parallel SGD via Nonuniform Quantization

abstract: As the size and complexity of models and datasets grow, so does the need for communication-efficient variants of stochastic gradient descent that can be deployed to perform parallel model training. One popular communication-compression method for data-parallel SGD is QSGD (Alistarh et al., 2017), which quantizes and encodes gradients to reduce communication costs. The baseline variant of QSGD provides strong theoretical guarantees, however, for practical purposes, the authors proposed a heuristic variant which we call QSGDinf, which demonstrated impressive empirical gains for distributed training of large neural networks. In this paper, we build on this work to propose a new gradient quantization scheme, and show that it has both stronger theoretical guarantees than QSGD, and matches and exceeds the empirical performance of the QSGDinf heuristic and of other compression methods.


title: Topotaxis of active Brownian particles

abstract: Recent experimental studies have demonstrated that cellular motion can be directed by topographical gradients, such as those resulting from spatial variations in the features of a micropatterned substrate. This phenomenon, known as topotaxis, is especially prominent among cells persistently crawling within a spatially varying distribution of cell-sized obstacles. In this article we introduce a toy model of topotaxis based on active Brownian particles constrained to move in a lattice of obstacles, with space-dependent lattice spacing. Using numerical simulations and analytical arguments, we demonstrate that topographical gradients introduce a spatial modulation of the particles' persistence, leading to directed motion toward regions of higher persistence. Our results demonstrate that persistent motion alone is sufficient to drive topotaxis and could serve as a starting point for more detailed studies on self-propelled particles and cells.


title: Task-Assisted Domain Adaptation with Anchor Tasks

abstract: Some tasks, such as surface normals or single-view depth estimation, require per-pixel ground truth that is difficult to obtain on real images but easy to obtain on synthetic. However, models learned on synthetic images often do not generalize well to real images due to the domain shift. Our key idea to improve domain adaptation is to introduce a separate anchor task (such as facial landmarks) whose annotations can be obtained at no cost or are already available on both synthetic and real datasets. To further leverage the implicit relationship between the anchor and main tasks, we apply our \freeze technique that learns the cross-task guidance on the source domain with the final network layers, and use it on the target domain. We evaluate our methods on surface normal estimation on two pairs of datasets (indoor scenes and faces) with two kinds of anchor tasks (semantic segmentation and facial landmarks). We show that blindly applying domain adaptation or training the auxiliary task on only one domain may hurt performance, while using anchor tasks on both domains is better behaved. Our \freeze technique outperforms competing approaches, reaching performance in facial images on par with a recently popular surface normal estimation method using shape from shading domain knowledge.


title: A young and obscured AGN embedded in the giant radio galaxy Mrk 1498

abstract: Mrk 1498 is part of a sample of galaxies with extended emission line regions (extended outwards up to a distance of $\sim$7 kpc) suggested to be photo-ionized by an AGN that has faded away or that is still active but heavily absorbed. Interestingly, the nucleus of Mrk 1498 is at the center of two giant radio lobes with a projected linear size of 1.1 Mpc. Our multi-wavelength analysis reveals a complex nuclear structure, with a young radio source (Giga-hertz Peaked Spectrum) surrounded by a strong X-ray nuclear absorption, a mid-infrared spectrum that is dominated by the torus emission, plus a circum-nuclear extended emission in the [OIII] image (with radius of $\sim$ 1 kpc), most likely related to the ionization of the AGN, aligned with the small and large scale radio jet and extended also at X-rays. In addition a large-scale extended emission (up to $\sim$ 10 kpc) is only visible in [OIII]. These data show conclusive evidence of a heavily absorbed nucleus and has recently restarted its nuclear activity. To explain its complexity, we propose that Mrk 1498 is the result of a merging event or secular processes, such as a minor interaction, that has triggered the nuclear activity and produced tidal streams. The large-scale extended emission that gives place to the actual morphology could either be explained by star formation or outflowing material from the AGN.


title: Analyzing the Fine Structure of Distributions

abstract: One aim of data mining is the identification of interesting structures in data. For better analytical results, the basic properties of an empirical distribution, such as skewness and eventual clipping, i.e. hard limits in value ranges, need to be assessed. Of particular interest is the question of whether the data originate from one process or contain subsets related to different states of the data producing process. Data visualization tools should deliver a clear picture of the univariate probability density distribution (PDF) for each feature. Visualization tools for PDFs typically use kernel density estimates and include both the classical histogram, as well as the modern tools like ridgeline plots, bean plots and violin plots. If density estimation parameters remain in a default setting, conventional methods pose several problems when visualizing the PDF of uniform, multimodal, skewed distributions and distributions with clipped data, For that reason, a new visualization tool called the mirrored density plot (MD plot), which is specifically designed to discover interesting structures in continuous features, is proposed. The MD plot does not require adjusting any parameters of density estimation, which is what may make the use of this plot compelling particularly to non-experts. The visualization tools in question are evaluated against statistical tests with regard to typical challenges of explorative distribution analysis. The results of the evaluation are presented using bimodal Gaussian, skewed distributions and several features with already published PDFs. In an exploratory data analysis of 12 features describing quarterly financial statements, when statistical testing poses a great difficulty, only the MD plots can identify the structure of their PDFs. In sum, the MD plot outperforms the above mentioned methods.


title: Shallow Domain Adaptive Embeddings for Sentiment Analysis

abstract: This paper proposes a way to improve the performance of existing algorithms for text classification in domains with strong language semantics. We propose a domain adaptation layer learns weights to combine a generic and a domain specific (DS) word embedding into a domain adapted (DA) embedding. The DA word embeddings are then used as inputs to a generic encoder + classifier framework to perform a downstream task such as classification. This adaptation layer is particularly suited to datasets that are modest in size, and which are, therefore, not ideal candidates for (re)training a deep neural network architecture. Results on binary and multi-class classification tasks using popular encoder architectures, including current state-of-the-art methods (with and without the shallow adaptation layer) show the effectiveness of the proposed approach.


title: Build it Break it Fix it for Dialogue Safety: Robustness from Adversarial Human Attack

abstract: The detection of offensive language in the context of a dialogue has become an increasingly important application of natural language processing. The detection of trolls in public forums (Gal\'an-Garc\'ia et al., 2016), and the deployment of chatbots in the public domain (Wolf et al., 2017) are two examples that show the necessity of guarding against adversarially offensive behavior on the part of humans. In this work, we develop a training scheme for a model to become robust to such human attacks by an iterative build it, break it, fix it strategy with humans and models in the loop. In detailed experiments we show this approach is considerably more robust than previous systems. Further, we show that offensive language used within a conversation critically depends on the dialogue context, and cannot be viewed as a single sentence offensive detection task as in most previous work. Our newly collected tasks and methods will be made open source and publicly available.


title: Polygamy Inequalities for Qubit Systems

abstract: Entanglement polygamy, like entanglement monogamy, is a fundamental property of multipartite quantum states. We investigate the polygamy relations related to the concurrence $C$ and the entanglement of formation $E$ for general $n$-qubit states. We extend the results in [Phys. Rev. A 90, 024304 (2014)] from the parameter region $\alpha\leq0$ to $\alpha\leq\alpha_0$, where $0<\alpha_0\leq2$ for $C$, and $0<\alpha_0\leq\sqrt{2}$ for $E$.


title: Arrow diagrams on spherical curves and computations

abstract: We give a definition of an integer-valued function $\sum_i \alpha_i x ^*_i$ derived from arrow diagrams for the ambient isotopy classes of oriented spherical curves. Then, we introduce certain elements of the free $\mathbb{Z}$-module generated by the arrow diagrams with at most $l $ arrows, called relators of Type~($\check{\rm{I}}$) (($\check{\rm{SI\!I} }$), ($\check{\rm{WI\!I}}$), ($\check{\rm{SI\!I\!I}}$), or ($\check{\rm{ WI\!I\!I}}$), resp.), and introduce another function $\sum_i \alpha_i \tilde{x}^*_i$ to obtain $\sum_i \alpha_i x^*_i$. One of the main results shows that if $\sum_i \alpha_i \tilde{x}^*_i$ vanishes on finitely many relators of Type~($\check{\rm{I}}$) (($\check{\rm{SI\!I}}$) , ($\check{\rm{WI\!I}}$), ($\check{\rm{SI\!I\!I}}$), or ($\check{\rm{WI\! I\!I}}$), resp.), then $\sum_i \alpha_i \tilde{x}$ is invariant under the deformation of type $\rm{RI}$ (strong$\rm{RI\!I}$, weak$\rm{RI\!I}$, strong$\rm{RI\!I\!I}$, or weak$\rm{RI\!I\!I}$, resp.). The other main result is that we obtain functions of arrow diagrams with up to six arrows. This computation is done with the aid of computers.


title: A Reliable IoT-Based Embedded Health Care System for Diabetic Patients

abstract: This paper introduces a reliable health care system for diabetic patients based on the Internet of Things technology. A diabetic health care system with a hardware implementation is presented. The proposed work employs Alaris 8100 infusion pump, Keil LPC-1768 board, and IoT-cloud to monitor the diabetic patients. The security of diabetic data over the cloud and the communication channel between health care system components are considered as part of the main contributions of this work. Moreover, an easy way to control and monitor the diabetic insulin pump is implemented. The \mbox{patient\textquotesingle s} records are stored in the cloud using the Keil board that is connected to the infusion pump. The reliability of the proposed scheme is accomplished by testing the system for five performance characteristics (availability, confidentiality, integrity, authentication, and authorization). The Kiel board is embedded with Ethernet port and Cortex-M3 micro-controller that controls the insulin infusion pump. The secure hash algorithm and secure socket shell are employed to achieve the reliability components of the proposed scheme. The results show that the proposed design is reliable, secure and authentic according to different test experiments and a case study of the Markov model. Moreover, a 99.3\% availability probability has been achieved after analyzing the case study.


title: 3D Rigid Motion Segmentation with Mixed and Unknown Number of Models

abstract: Many real-world video sequences cannot be conveniently categorized as general or degenerate; in such cases, imposing a false dichotomy in using the fundamental matrix or homography model for motion segmentation on video sequences would lead to difficulty. Even when we are confronted with a general scene-motion, the fundamental matrix approach as a model for motion segmentation still suffers from several defects, which we discuss in this paper. The full potential of the fundamental matrix approach could only be realized if we judiciously harness information from the simpler homography model. From these considerations, we propose a multi-model spectral clustering framework that synergistically combines multiple models (homography and fundamental matrix) together. We show that the performance can be substantially improved in this way. For general motion segmentation tasks, the number of independently moving objects is often unknown a priori and needs to be estimated from the observations. This is referred to as model selection and it is essentially still an open research problem. In this work, we propose a set of model selection criteria balancing data fidelity and model complexity. We perform extensive testing on existing motion segmentation datasets with both segmentation and model selection tasks, achieving state-of-the-art performance on all of them; we also put forth a more realistic and challenging dataset adapted from the KITTI benchmark, containing real-world effects such as strong perspectives and strong forward translations not seen in the traditional datasets.


title: Matrix Lie Maps and Neural Networks for Solving Differential Equations

abstract: The coincidence between polynomial neural networks and matrix Lie maps is discussed in the article. The matrix form of Lie transform is an approximation of the general solution of the nonlinear system of ordinary differential equations. It can be used for solving systems of differential equations more efficiently than traditional step-by-step numerical methods. Implementation of the Lie map as a polynomial neural network provides a tool for both simulation and data-driven identification of dynamical systems. If the differential equation is provided, training a neural network is unnecessary. The weights of the network can be directly calculated from the equation. On the other hand, for data-driven system learning, the weights can be fitted without any assumptions in view of differential equations. The proposed technique is discussed in the examples of both ordinary and partial differential equations. The building of a polynomial neural network that simulates the Van der Pol oscillator is discussed. For this example, we consider learning the dynamics from a single solution of the system. We also demonstrate the building of the neural network that describes the solution of Burgers' equation that is a fundamental partial differential equation.

categories: cs.NE cs.NA math.DS math.NA physics.comp-ph


title: Batch 1: Definition of several Weather & Climate Dwarfs

abstract: This document is one of the deliverable reports created for the ESCAPE project. ESCAPE stands for Energy-efficient Scalable Algorithms for Weather Prediction at Exascale. The project develops world-class, extreme-scale computing capabilities for European operational numerical weather prediction and future climate models. This is done by identifying weather & climate dwarfs which are key patterns in terms of computation and communication (in the spirit of the Berkeley dwarfs). These dwarfs are then optimised for different hardware architectures (single and multi-node) and alternative algorithms are explored. Performance portability is addressed through the use of domain specific languages. This deliverable contains the description of the characteristics of the weather & climate dwarfs that form key functional components of prediction models in terms of the science that they encapsulate and in terms of computational cost they impose on the forecast production. The ESCAPE work flow between work packages centres on these dwarfs and hence their selection, their performance assessment, code adaptation and optimization is crucial for the success of the project. At this stage of ESCAPE, a selection of established and new dwarfs has been made, their documentation been compiled and the software been made available on the software exchange platform. The selection of dwarfs will be extended throughout the course of the project (see Deliverable D1.2). The current selection includes the spectral transforms, the cloud microphysics scheme, two and three-dimensional elliptic solvers, a bi-Fourier spectral transform, an interpolation needed for the semi-Lagrangian advection scheme and a first version of the semi-Lagrangian advection scheme itself. This deliverable includes their scientific description and the guidance for installation, execution and testing.


title: Optimal Paired Comparison Experiments for Second-Order Interactions

abstract: In real life situations often paired comparisons involving alternatives of either full or partial profiles to mitigate cognitive burden are presented. For this situation the problem of finding optimal designs is considered in the presence of second-order interactions when all attributes have general common number of levels.


title: Development of Atlas, a flexible data structure framework

abstract: This document is one of the deliverable reports created for the ESCAPE project. ESCAPE stands for Energy-efficient Scalable Algorithms for Weather Prediction at Exascale. The project develops world-class, extreme-scale computing capabilities for European operational numerical weather prediction and future climate models. This is done by identifying Weather & Climate dwarfs which are key patterns in terms of computation and communication (in the spirit of the Berkeley dwarfs). These dwarfs are then optimised for different hardware architectures (single and multi-node) and alternative algorithms are explored. Performance portability is addressed through the use of domain specific languages. In this deliverable report, we present Atlas, a new software library that is currently being developed at the European Centre for Medium-Range Weather Forecasts (ECMWF), with the scope of handling data structures required for NWP applications in a flexible and massively parallel way. Atlas provides a versatile framework for the future development of efficient NWP and climate applications on emerging HPC architectures. The applications range from full Earth system models, to specific tools required for post-processing weather forecast products. Atlas provides data structures for building various numerical strategies to solve equations on the sphere or limited area's on the sphere. These data structures may contain a distribution of points (grid) and, possibly, a composition of elements (mesh), required to implement the numerical operations required. Atlas can also represent a given field within a specific spatial projection. Atlas is capable of mapping fields between different grids as part of pre- and post-processing stages or as part of coupling processes whose respective fields are discretised on different grids or meshes.


title: Optimal $2^K$ Paired Comparison Designs for Third-Order Interactions

abstract: In psychological research often paired comparisons are used in which either full or partial profiles of the alternatives described by a common set of two-level attributes are presented. For this situation the problem of finding optimal designs is considered in the presence of third-order interactions.


title: Additional key features required for different directives based porting approaches

abstract: This document is one of the deliverable reports created for the ESCAPE project. ESCAPE stands for Energy-efficient Scalable Algorithms for Weather Prediction at Exascale. The project develops world-class, extreme-scale computing capabilities for European operational numerical weather prediction and future climate models. This is done by identifying Weather & Climate dwarfs which are key patterns in terms of computation and communication (in the spirit of the Berkeley dwarfs). These dwarfs are then optimised for different hardware architectures (single and multi-node) and alternative algorithms are explored. Performance portability is addressed through the use of domain specific languages. This report summarizes key features required for OpenMP and OpenACC directives based on experience in the ESCAPE project. For OpenMP, the latest public draft standard, 5.0, contains the deep copy and multi-level memory features desired; for OpenACC, Technical Report 16 summarizes ongoing discussions beyond the Standard version 2.6. This document includes a summary of our recommendations on this approach. Additional work is also desirable in coordinating the runtime and debugging when both OpenACC and OpenMP directives are used; in particular the interoperability of the new OPDT debugging interface for OpenMP and its semantics in a mixed-directive program.


title: Report on the performance portability demonstrated for the relevant Weather & Climate Dwarfs

abstract: This document is one of the deliverable reports created for the ESCAPE project. ESCAPE stands for Energy-efficient Scalable Algorithms for Weather Prediction at Exascale. The project develops world-class, extreme-scale computing capabilities for European operational numerical weather prediction and future climate models. This is done by identifying Weather & Climate dwarfs which are key patterns in terms of computation and communication (in the spirit of the Berkeley dwarfs). These dwarfs are then optimised for different hardware architectures (single and multi-node) and alternative algorithms are explored. Performance portability is addressed through the use of domain specific languages. This deliverable provides an evaluation of the work performed within ESCAPE to port different dwarfs to accelerators, using different programming models. A key metric of the evaluation is the performance portability of the resulting porting efforts. Portability means that a single source code containing the numerical operators can be compiled and run in multiple architectures, while performance portability additionally requires that the single source code runs efficiently in all the different architectures. As results of other deliverables like D2.1, D2.4 ESCAPE provides a collection of dwarfs ported to different computing architectures like traditional CPUs, Intel XeonPhi and NVIDIA GPUs. Additionally D3.3 went through an optimization process to obtain efficient and energy efficient dwarfs. In this deliverable we present a review of the different programming models employed and their use to port various dwarfs of ESCAPE. A final evaluation of the different approaches based on different metrics like performance portability, readability of the numerical methods, efforts to port a dwarf and efficiency of the implementation obtained is reported.


title: Recommendations and specifications for data scope analysis tools

abstract: This document is one of the deliverable reports created for the ESCAPE project. ESCAPE stands for Energy-efficient Scalable Algorithms for Weather Prediction at Exascale. The project develops world-class, extreme-scale computing capabilities for European operational numerical weather prediction and future climate models. This is done by identifying Weather & Climate dwarfs which are key patterns in terms of computation and communication (in the spirit of the Berkeley dwarfs). These dwarfs are then optimised for different hardware architectures (single and multi-node) and alternative algorithms are explored. Performance portability is addressed through the use of domain specific languages. In today's computer architectures, moving data is considerably more time- and energy consuming than computing on this data. One of the key performance optimizations for any application is therefore to minimize data motion and maximize data reuse. Especially on modern supercomputers with very complex and deep memory hierarchies, it is mandatory to take data locality into account. Especially when targeting accelerators with directive systems like OpenACC or OpenMP, identifying data scope, access type and data reuse are critical to minimize the data transfers from and to the accelerator. Unfortunately, manually identifying data locality information in complex code bases can be a time consuming task and tool support is therefore desirable. In this report we summarize the results of a survey of currently available tools that support software developers and performance engineers with data locality information in complex code bases like numerical weather prediction (NWP) or climate simulation applications. Based on the survey results we then recommend a tool and specify some extensions for a tool to solve the problems encountered in an NWP application.


title: Performance report and optimized implementation of Weather & Climate Dwarfs on GPU, MIC and Optalysys Optical Processor

abstract: This document is one of the deliverable reports created for the ESCAPE project. ESCAPE stands for Energy-efficient Scalable Algorithms for Weather Prediction at Exascale. The project develops world-class, extreme-scale computing capabilities for European operational numerical weather prediction and future climate models. This is done by identifying Weather & Climate dwarfs which are key patterns in terms of computation and communication (in the spirit of the Berkeley dwarfs). These dwarfs are then optimised for different hardware architectures (single and multi-node) and alternative algorithms are explored. Performance portability is addressed through the use of domain specific languages. Here we summarize the work performed on optimizations of the dwarfs on CPUs, Xeon Phi, GPUs and on the Optalysys optical processor. We limit ourselves to a subset of the dwarf configurations and to problem sizes small enough to execute on a single node. Also, we use time-to-solution as the main performance metric. Multi-node optimizations of the dwarfs and energy-specific optimizations are beyond the scope of this report and will be described in Deliverable D3.4. To cover the important algorithmic motifs we picked dwarfs related to the dynamical core as well as column physics. Specifically, we focused on the formulation relevant to spectral codes like ECMWF's IFS code. The main findings of this report are: (a) Acceleration of 1.1x - 2.5x of the dwarfs on CPU based systems using compiler directives, (b) order of magnitude acceleration of the dwarfs on GPUs (23x for spectral transform, 9x for MPDATA) using data locality optimizations and (c) demonstrated feasibility of a spectral transform in a purely optical fashion.


title: Performance report and optimized implementations of Weather & Climate dwarfs on multi-node systems

abstract: This document is one of the deliverable reports created for the ESCAPE project. ESCAPE stands for Energy-efficient Scalable Algorithms for Weather Prediction at Exascale. The project develops world-class, extreme-scale computing capabilities for European operational numerical weather prediction and future climate models. This is done by identifying Weather & Climate dwarfs which are key patterns in terms of computation and communication (in the spirit of the Berkeley dwarfs). These dwarfs are then optimised for different hardware architectures (single and multi-node) and alternative algorithms are explored. Performance portability is addressed through the use of domain specific languages. Here we summarize the work performed on optimizations of the dwarfs focusing on CPU multi-nodes and multi-GPUs. We limit ourselves to a subset of the dwarf configurations chosen by the consortium. Intra-node optimizations of the dwarfs and energy-specific optimizations have been described in Deliverable D3.3. To cover the important algorithmic motifs we picked dwarfs related to the dynamical core as well as column physics. Specifically, we focused on the formulation relevant to spectral codes like ECMWF's IFS code. The main findings of this report are: (a) Up-to 30% performance gain with CPU based multi-node systems compared to optimized version of dwarfs from task 3.3 (see D3.3), (b) up to 10X performance gain on multiple GPUs from optimizations to keep data resident on the GPU and enable fast inter-GPU communication mechanisms, and (c) multi-GPU systems which feature a high-bandwidth all-to-all interconnect topology with NVLink/NVSwitch hardware are particularly well suited to the algorithms.


title: Projections of achievable performance for Weather & Climate Dwarfs, and for entire NWP applications, on hybrid architectures

abstract: This document is one of the deliverable reports created for the ESCAPE project. ESCAPE stands for Energy-efficient Scalable Algorithms for Weather Prediction at Exascale. The project develops world-class, extreme-scale computing capabilities for European operational numerical weather prediction and future climate models. This is done by identifying Weather & Climate dwarfs which are key patterns in terms of computation and communication (in the spirit of the Berkeley dwarfs). These dwarfs are then optimised for different hardware architectures (single and multi-node) and alternative algorithms are explored. Performance portability is addressed through the use of domain specific languages. This deliverable contains the description of the performance and energy models for the selected Weather & Climate dwarfs for different hardware architectures, multinode with GPU accelerators in particular. Presented performance models are extension to model provided in Deliverable 3.2. With some further enhancements, they are incorporated in the DCworms simulator. In particular, extended models allow to predict computational and energy performance on different architectures: single and multinodes, equipped with CPUs and GPUs accelerators. This allows to provide feasible performance projection at system scale.


title: Low frequency X-ray timing with Gaussian processes and reverberation in the radio-loud AGN 3C 120

abstract: A framework is developed to perform Fourier-domain timing analysis on X-ray light curves with gaps, employing Gaussian processes to model the probability distribution underlying the observed time series from which continuous samples can be drawn. A technique is developed to measure X-ray reverberation from the inner regions of accretion discs around black holes in the low frequency components of the variability, on timescales longer than can be probed employing standard Fourier techniques. This enables X-ray reverberation experiments to be performed using data from satellites in low-Earth orbit such as NICER, NuSTAR and the proposed X-ray timing mission STROBE-X, and enables long timescale reverberation around higher mass AGN to be measured by combining multiple observations. Gaussian processes are applied to observations of the broad line radio galaxy 3C120 spanning two orbits with XMM-Newton to measure the relative time lags of successive X-ray energy bands. The lag-energy spectrum between 5E-6 and 3E-5Hz, estimated using Gaussian processes, reveals X-ray reverberation from the inner accretion disc for the first time in this radio-loud AGN. Time lags in the relativistically broadened iron K line are significantly detected. The core of the line lags behind the continuum by (3800 +/- 1500)s, suggesting a scale height of the corona of (13 +/- 8)rg above the disc. The ability to compare the structure of coronae in radio loud AGN to their radio quiet counterparts will yield important insight into the mechanisms by which black holes are able to launch jets.


title: Gravitational waves, CMB polarization, and the Hubble tension

abstract: The discrepancy between the Hubble parameter inferred from local measurements and that from the cosmic microwave background (CMB) has motivated careful scrutiny of the assumptions that enter both analyses. Here we point out that the location of the recombination peak in the CMB B-mode power spectrum is determined by the light horizon at the surface of last scatter and thus provides an alternative early-Universe standard ruler. It can thus be used as a cross-check for the standard ruler inferred from the acoustic peaks in the CMB temperature power spectrum and to test various explanations for the Hubble tension. The measurement can potentially be carried out with a precision of $\lesssim2\%$ with stage-IV B-mode experiments. The measurement can also be used to measure the propagation speed of gravitational waves in the early Universe.


title: Parallel implementation of high-fidelity multi-qubit gates with neutral atoms

abstract: We report the implementation of universal two- and three-qubit entangling gates on neutral atom qubits encoded in long-lived hyperfine ground states. The gates are mediated by excitation to strongly interacting Rydberg states, and are implemented in parallel on several clusters of atoms in a one-dimensional array of optical tweezers. Specifically, we realize the controlled-phase gate, enacted by a novel, fast protocol involving only global coupling of two qubits to Rydberg states. We benchmark this operation by preparing Bell states with fidelity $\mathcal{F} \ge 95.0(2)\%$, and extract gate fidelity $\ge 97.4(3)\%,$ averaged across five atom pairs. In addition, we report a proof-of-principle implementation of the three-qubit Toffoli gate, in which two control atoms simultaneously constrain the behavior of one target atom. These experiments demonstrate key ingredients for high-fidelity quantum information processing in a scalable neutral atom platform.


title: The impact of stars stripped in binaries on the integrated spectra of stellar populations

abstract: Stars stripped of their envelopes from interaction with a binary companion emit a significant fraction of their radiation as ionizing photons. They are potentially important stellar sources of ionizing radiation, however, they are still often neglected in spectral synthesis simulations or simulations of stellar feedback. We modeled the radiative contribution from stripped stars by using detailed evolutionary and spectral models. We estimated their impact on the integrated spectra and specifically on the emission rates of HI-, HeI-, and HeII-ionizing photons from stellar populations. We find that stripped stars have the largest impact on the ionizing spectrum of a population in which star formation halted several Myr ago. In such stellar populations, stripped stars dominate the emission of ionizing photons, mimicking a younger stellar population in which massive stars are still present. Our models also suggest that stripped stars have harder ionizing spectra than massive stars. The additional ionizing radiation affects observable properties that are related to the ionizing emission from stellar populations. In co-eval stellar populations, the ionizing radiation from stripped stars increases the ionization parameter and the production efficiency of HI-ionizing photons. They also cause high values for these parameters for about ten times longer than what is predicted for massive stars. The hard ionizing radiation from stripped stars likely introduces a characteristic ionization structure of the nebula, which leads to the emission of highly ionized elements such as O$^{2+}$ and C$^{3+}$. We, therefore, expect that the presence of stripped stars affects the location in the BPT diagram and the diagnostic ratio of OIII to OII nebular emission lines. Our models are publicly available through CDS database and on the Starburst99 website.