title: Spin-diffusion model for micromagnetics in the limit of long times

abstract: In this paper, we consider spin-diffusion Landau-Lifshitz-Gilbert equations (SDLLG), which consist of the time-dependent Landau-Lifshitz-Gilbert (LLG) equation coupled with a time-dependent diffusion equation for the electron spin accumulation. The model takes into account the diffusion process of the spin accumulation in the magnetization dynamics of ferromagnetic multilayers. We prove that in the limit of long times, the system reduces to simpler equations in which the LLG equation is coupled to a nonlinear and nonlocal steady-state equation, referred to as SLLG. As a by-product, the existence of global weak solutions to the SLLG equation is obtained. Moreover, we prove weak-strong uniqueness of solutions of SLLG, i.e., all weak solutions coincide with the (unique) strong solution as long as the latter exists in time. The results provide a solid mathematical ground to the qualitative behavior originally predicted by Zhang, Levy, and Fert in [Physical Review Letters 88 (2002)] in ferromagnetic multilayers.


title: Turbulent Details Simulation for SPH Fluids via Vorticity Refinement

abstract: A major issue in Smoothed Particle Hydrodynamics (SPH) approaches is the numerical dissipation during the projection process, especially under coarse discretizations. High-frequency details, such as turbulence and vortices, are smoothed out, leading to unrealistic results. To address this issue, we introduce a Vorticity Refinement (VR) solver for SPH fluids with negligible computational overhead. In this method, the numerical dissipation of the vorticity field is recovered by the difference between the theoretical and the actual vorticity, so as to enhance turbulence details. Instead of solving the Biot-Savart integrals, a stream function, which is easier and more efficient to solve, is used to relate the vorticity field to the velocity field. We obtain turbulence effects of different intensity levels by changing an adjustable parameter. Since the vorticity field is enhanced according to the curl field, our method can not only amplify existing vortices, but also capture additional turbulence. Our VR solver is straightforward to implement and can be easily integrated into existing SPH methods.


title: Connected components in the invariably generating graph of a finite group

abstract: We prove that the invariably generating graph of a finite group can have an arbitrarily large number of connected components with at least two vertices.


title: Observation and measurement of forward proton scattering in association with lepton pairs produced via the photon fusion mechanism at ATLAS

abstract: The observation of forward proton scattering in association with lepton pairs ($e^+e^-+p$ or $\mu^+\mu^-+p$) produced via photon fusion is presented. The scattered proton is detected by the ATLAS Forward Proton spectrometer while the leptons are reconstructed by the central ATLAS detector. Proton-proton collision data recorded in 2017 at a center-of-mass energy of $\sqrt{s} = 13$ TeV are analyzed, corresponding to an integrated luminosity of 14.6 fb$^{-1}$. A total of 57 (123) candidates in the $ee+p$ ($\mu\mu+p$) final state are selected, allowing the background-only hypothesis to be rejected with a significance exceeding five standard deviations in each channel. Proton-tagging techniques are introduced for cross-section measurements in the fiducial detector acceptance, corresponding to $\sigma_{ee+p}$ = 11.0 $\pm$ 2.6 (stat.) $\pm$ 1.2 (syst.) $\pm$ 0.3 (lumi.) fb and $\sigma_{\mu\mu+p}$ = 7.2 $\pm$ 1.6 (stat.) $\pm$ 0.9 (syst.) $\pm$ 0.2 (lumi.) fb in the dielectron and dimuon channel, respectively.


title: A study of the decays of $S-$wave $\bar D^\ast K^\ast$ hadronic molecules: the scalar $X_0(2900)$ and its spin partners $X_{J(J=1,2)}$

abstract: In this work, we investigated the decays of the fully open-flavor tetraquark state $X_0(2900)$ which was observed by the LHCb Collaboration very recently. Here, the $X_0(2900)$ was assigned as a $S-$wave $\bar D^\ast K^\ast$ hadronic molecule with $I=0$, and the effective lagrangian approach was applied to estimate the partial decay widths. Moreover, we also predicted the decay behaviors of the other unobserved $X_{J(J=1,2)}$, which were the spin partners of the $X_0(2900)$ in the $S-$wave $\bar D^\ast K^\ast$ picture. It was pointed out that the $X_1$ state with $I=0$ was a broad state with the width more than one hundred MeV, while another $X_2$ state with $I=0$ was a narrow state with the width approaching half of that for the $X_0(2900)$. In addition, our results also showed that the $\bar D^\ast K$ mode was expected to be the dominant decay mode for both $X_1$ and $X_2$. Searching for those unobserved $X_{J(J=1,2)}$ in the future experiments might be helpful to understand the nature of $X_0(2900)$.


abstract: Most of the contemporary approaches for multi-hop Natural Language Inference (NLI) construct explanations considering each test case in isolation. However, this paradigm is known to suffer from semantic drift, a phenomenon that causes the construction of spurious explanations leading to wrong conclusions. In contrast, this paper proposes an abductive framework for multi-hop NLI exploring the retrieve-reuse-refine paradigm in Case-Based Reasoning (CBR). Specifically, we present Case-Based Abductive Natural Language Inference (CB-ANLI), a model that addresses unseen inference problems by analogical transfer of prior explanations from similar examples. We empirically evaluate the abductive framework on commonsense and scientific question answering tasks, demonstrating that CB-ANLI can be effectively integrated with sparse and dense pre-trained encoders to improve multi-hop inference, or adopted as an evidence retriever for Transformers. Moreover, an empirical analysis of semantic drift reveals that the CBR paradigm boosts the quality of the most challenging explanations, a feature that has a direct impact on robustness and accuracy in downstream inference tasks.


title: Comparison of SIA Defect Morphologies from Different Interatomic Potentials for Collision Cascades in W

abstract: The morphology of defects formed in collision cascades is an essential aspect of the subsequent evolution of the microstructure. The morphological composition of a defect decides its stability, interaction, and migration properties. We compare the defect morphologies in the primary radiation damage caused by high energy collision cascades simulated using three different interatomic potentials in W. An automated method to identify morphologies of defects is used. While most defects form 1/2\3 dislocation loops, other specific morphologies include \1 dislocation loops, multiple loops clustered together, rings corresponding to C15 configuration and its constituent structures, and a combination of rings and dislocations. The analysis quantifies the distribution of defects among different morphologies and the size distribution of each morphology. We show that the disagreement between predictions of the different potentials regarding defect morphology is much stronger than the differences in predicted defect numbers.


title: Deformed Shape Invariant Superpotentials in Quantum Mechanics and Expansions in Powers of $\hbar$

abstract: We show that the method developed by Gangopadhyaya, Mallow, and their coworkers to deal with (translationally) shape invariant potentials in supersymmetric quantum mechanics and consisting in replacing the shape invariance condition, which is a difference-differential equation, by an infinite set of partial differential equations can be generalized to deformed shape invariant potentials in deformed supersymmetric quantum mechanics. The extended method is illustrated by several examples, corresponding both to $\hbar$-independent superpotentials and to a superpotential explicitly depending on $\hbar$.


title: Weighted Tiling Systems for Graphs: Evaluation Complexity

abstract: We consider weighted tiling systems to represent functions from graphs to a commutative semiring such as the Natural semiring or the Tropical semiring. The system labels the nodes of a graph by its states, and checks if the neighbourhood of every node belongs to a set of permissible tiles, and assigns a weight accordingly. The weight of a labeling is the semiring-product of the weights assigned to the nodes, and the weight of the graph is the semiring-sum of the weights of labelings. We show that we can model interesting algorithmic questions using this formalism - like computing the clique number of a graph or computing the permanent of a matrix. The evaluation problem is, given a weighted tiling system and a graph, to compute the weight of the graph. We study the complexity of the evaluation problem and give tight upper and lower bounds for several commutative semirings. Further we provide an efficient evaluation algorithm if the input graph is of bounded tree-width.


title: Designing of Magnetic MAB Phases for Energy Applications

abstract: Based on high-throughput density functional theory calculations, we performed screening for stable magnetic MAB compounds and predicted potential strong magnets for permanent magnet and magnetocaloric applications. The thermodynamical, mechanical, and dynamical stabilities are systematically evaluated, resulting in 21 unreported compounds on the convex hull, and 434 materials being metastable considering convex hull tolerance to be 100 meV/atom. Analysis based on the Hume-Rothery rules revealed that the valence electron concentration and size factor difference are of significant importance in determining the stability, with good correspondence with the local atomic bonding. We found 71 compounds with the absolute value of magneto-crystalline anisotropy energy above 1.0 MJ/m$^3$ and 23 compounds with a uniaxial anisotropy greater than 0.4 MJ/m$^3$, which are potential gap magnets. Based on the magnetic deformation proxy, 99 compounds were identified as potential materials with interesting magnetocaloric performance.


title: Time-domain and Frequency-domain Signals and their Analysis

abstract: Depending on the application people use time-domain or frequency-domain signals in order to measure or describe processes. First we will look at the definition of these terms, produce some mathematical background and then apply the tools to measurements made in the accelerator domain. We will first look at signals produced by a single bunch passing once through a detector (transfer line, linac), then periodic single bunch passages (circular accelerator) and at the end multi-bunch passages in a circular accelerator.


title: Demographic Influences on Contemporary Art with Unsupervised Style Embeddings

abstract: Computational art analysis has, through its reliance on classification tasks, prioritised historical datasets in which the artworks are already well sorted with the necessary annotations. Art produced today, on the other hand, is numerous and easily accessible, through the internet and social networks that are used by professional and amateur artists alike to display their work. Although this art, yet unsorted in terms of style and genre, is less suited for supervised analysis, the data sources come with novel information that may help frame the visual content in equally novel ways. As a first step in this direction, we present contempArt, a multi-modal dataset of exclusively contemporary artworks. contempArt is a collection of paintings and drawings, a detailed graph network based on social connections on Instagram and additional socio-demographic information; all attached to 442 artists at the beginning of their career. We evaluate three methods suited for generating unsupervised style embeddings of images and correlate them with the remaining data. We find no connections between visual style on the one hand and social proximity, gender, and nationality on the other.


title: Fast reaction limits via $\Gamma$-convergence of the Flux Rate Functional

abstract: We study the convergence of a sequence of evolution equations for measures supported on the nodes of a graph. The evolution equations themselves can be interpreted as the forward Kolmogorov equations of Markov jump processes, or equivalently as the equations for the concentrations in a network of linear reactions. The jump rates or reaction rates are divided in two classes; `slow' rates are constant, and `fast' rates are scaled as~$1/\epsilon$, and we prove the convergence in the fast-reaction limit $\epsilon\to0$. We establish a $\Gamma$-convergence result for the rate functional in terms of both the concentration at each node and the flux over each edge (the level-2.5 rate function). The limiting system is again described by a functional, and characterizes both fast and slow fluxes in the system. This method of proof has three advantages. First, no condition of detailed balance is required. Secondly, the formulation in terms of concentration and flux leads to a short and simple proof of the $\Gamma$-convergence; the price to pay is a more involved compactness proof. Finally, the method of proof deals with approximate solutions, for which the functional is not zero but small, without any changes.


title: FAN: Frequency Aggregation Network for Real Image Super-resolution

abstract: Single image super-resolution (SISR) aims to recover the high-resolution (HR) image from its low-resolution (LR) input image. With the development of deep learning, SISR has achieved great progress. However, It is still a challenge to restore the real-world LR image with complicated authentic degradations. Therefore, we propose FAN, a frequency aggregation network, to address the real-world image super-resolu-tion problem. Specifically, we extract different frequencies of the LR image and pass them to a channel attention-grouped residual dense network (CA-GRDB) individually to output corresponding feature maps. And then aggregating these residual dense feature maps adaptively to recover the HR image with enhanced details and textures. We conduct extensive experiments quantitatively and qualitatively to verify that our FAN performs well on the real image super-resolution task of AIM 2020 challenge. According to the released final results, our team SR-IM achieves the fourth place on the X4 track with PSNR of 31.1735 and SSIM of 0.8728.


title: Modelling the Microstructure and Stress in Dense Suspensions Under Inhomogeneous Flow

abstract: Under inhomogeneous flow, dense suspensions exhibit complex behaviour that violates the conventional homogenous rheology. Specifically, one finds flowing regions with a macroscopic friction coefficient below the yielding criterion, and volume fraction above the jamming criterion. We demonstrate the underlying physics by incorporating shear rate fluctuations into a recently proposed tensor model for the microstructure and stress, and applying the model to an inhomogeneous flow problem. The model predictions agree qualitatively with particle-based simulations.


title: Optimal Control of Industrial Assembly Lines

abstract: This paper discusses the problem of assembly line control and introduces an optimal control formulation that can be used to improve the performance of the assembly line, in terms of cycle time minimization, resources' utilization, etc. A deterministic formulation of the problem is introduced, based on mixed-integer linear programming. A simple numerical simulation provides a first proof of the proposed concept.


title: Out-of-plane transport of 1T-TaS2/graphene-based van der Waals heterostructures

abstract: Due to their anisotropy, layered materials are excellent candidates for studying the interplay between the in-plane and out-of-plane entanglement in strongly correlated systems. A relevant example is provided by 1T-TaS2, which exhibits a multifaceted electronic and magnetic scenario due to the existence of several charge density wave (CDW) configurations. It includes quantum hidden phases, superconductivity and exotic quantum spin liquid (QSL) states, which are highly dependent on the out-of-plane stacking of the CDW. In this system, the interlayer stacking of the CDW is crucial for the interpretation of the underlying electronic and magnetic phase diagram. Here, thin-layers of 1T-TaS2 are integrated in vertical van der Waals heterostructures based on few-layer graphene (FLG) contacts and their electrical transport properties are measured. Different activation energies in the conductance and a gap at the Fermi level are clearly observed. Our experimental findings are supported by fully self-consistent DFT+U calculations, which evidence the presence of an energy gap in the few-layer limit, not necessarily coming from the formation of out-of-plane spin-paired bilayers at low temperatures, as previously proposed for the bulk. These results highlight dimensionality as a key effect for understanding the properties of 1T-TaS2 and opens the door to the possible experimental realization of low-dimensional QSLs.


title: Explainable Deep Reinforcement Learning for UAV Autonomous Navigation

abstract: Autonomous navigation in unknown complex environment is still a hard problem, especially for small Unmanned Aerial Vehicles (UAVs) with limited computation resources. In this paper, a neural network-based reactive controller is proposed for a quadrotor to fly autonomously in unknown outdoor environment. The navigation controller makes use of only current sensor data to generate the control signal without any optimization or configuration space searching, which reduces both memory and computation requirement. The navigation problem is modelled as a Markov Decision Process (MDP) and solved using deep reinforcement learning (DRL) method. Specifically, to get better understanding of the trained network, some model explanation methods are proposed. Based on the feature attribution, each decision making result during flight is explained using both visual and texture explanation. Moreover, some global analysis are also provided for experts to evaluate and improve the trained neural network. The simulation results illustrated the proposed method can make useful and reasonable explanation for the trained model, which is beneficial for both non-expert users and controller designer. Finally, the real world tests shown the proposed controller can navigate the quadrotor to goal position successfully and the reactive controller performs much faster than some conventional approach under the same computation resource.


title: Wasserstein Distributionally Robust Inverse Multiobjective Optimization

abstract: Inverse multiobjective optimization provides a general framework for the unsupervised learning task of inferring parameters of a multiobjective decision making problem (DMP), based on a set of observed decisions from the human expert. However, the performance of this framework relies critically on the availability of an accurate DMP, sufficient decisions of high quality, and a parameter space that contains enough information about the DMP. To hedge against the uncertainties in the hypothetical DMP, the data, and the parameter space, we investigate in this paper the distributionally robust approach for inverse multiobjective optimization. Specifically, we leverage the Wasserstein metric to construct a ball centered at the empirical distribution of these decisions. We then formulate a Wasserstein distributionally robust inverse multiobjective optimization problem (WRO-IMOP) that minimizes a worst-case expected loss function, where the worst case is taken over all distributions in the Wasserstein ball. We show that the excess risk of the WRO-IMOP estimator has a sub-linear convergence rate. Furthermore, we propose the semi-infinite reformulations of the WRO-IMOP and develop a cutting-plane algorithm that converges to an approximate solution in finite iterations. Finally, we demonstrate the effectiveness of our method on both a synthetic multiobjective quadratic program and a real world portfolio optimization problem.


title: Accounting for anomalous velocity curve in galaxies as a possible frame dragging effect

abstract: We investigate the possibility of accounting for the observed anomalous velocities of stars in galaxies to be a result of a dynamic prescription for having a galaxy generate and prescribe its own inertial frame. It is demonstrated that within observational uncertainties, an overall rotation of an inertial frame dragged by a galaxy is a comfortable fit to the observed anomalous velocity curves. A case for looking for a natural law for prescribing an inertial frame by a mass distribution over a finite (galaxy size) range is made out.


title: One Reflection Suffice

abstract: Orthogonal weight matrices are used in many areas of deep learning. Much previous work attempt to alleviate the additional computational resources it requires to constrain weight matrices to be orthogonal. One popular approach utilizes *many* Householder reflections. The only practical drawback is that many reflections cause low GPU utilization. We mitigate this final drawback by proving that *one* reflection is sufficient, if the reflection is computed by an auxiliary neural network.


title: A Diffuse Interface Model of Reactive-fluids and Solid-dynamics

abstract: This article presents a multi-physics methodology for the numerical simulation of physical systems that involve the non-linear interaction of multi-phase reactive fluids and elastoplastic solids, inducing high strain-rates and high deformations. Each state of matter is governed by a single system of non-linear, inhomogeneous partial differential equations, which are solved simultaneously on the same computational grid, and do not require special treatment of immersed boundaries. To this end, the governing equations for solid and reactive multiphase fluid mechanics are written in the same mathematical form and are discretised on a regular Cartesian mesh. All phase and material boundaries are treated as diffuse interfaces. An interface-steepening technique is employed at material boundaries to keep interfaces sharp whilst maintaining the conservation properties of the system. These algorithms are implemented in a highly-parallelised hierarchical adaptive mesh refinement platform, and are verified and validated using numerical and experimental benchmarks. Results indicate very good agreement with experiment and an improvement of numerical performance compared to certain existing Eulerian methods, without loss of conservation.


title: Multimillijoule terahertz radiation from laser interactions with microplasma-waveguides

abstract: When a relativistic, femtosecond laser pulse enters a waveguide, the pulse energy is coupled into waveguide optical modes. The longitudinal laser field effectively accelerates electrons along the axis of the channel, while the asymmetric transverse electromagnetic fields tend to expel fast electrons radially outwards. At the exit of the waveguide, the $\sim$${\rm nC}$, $\sim$$10\ {\rm MeV}$ electron bunch converts its energy to a $\sim$$10\ {\rm mJ}$ terahertz (THz) laser pulse through coherent diffraction radiation. In this paper, we present 3D particle-in-cell simulations and theoretical analyses of the aforementioned interaction process. We investigate the process of longitudinal acceleration and radial expulsion of fast electrons, as well as the dependence of the properties of the resulting THz radiation on laser and plasma parameters and the effects of the preplasma. The simulation results indicate that the conversion efficiency of energy can be over $5\%$ if the waveguide length is optimal and a high contrast pump laser is used. These results guide the design of more intense and powerful THz sources.


title: Varieties of Tropical Ideals are Balanced

abstract: Tropical ideals, introduced in arXiv:1609.03838, define subschemes of tropical toric varieties. We prove that the top-dimensional parts of their varieties are balanced polyhedral complexes of the same dimension as the ideal. This means that every subscheme of a tropical toric variety defined by a tropical ideal has an associated class in the Chow ring of the toric variety. A key tool in the proof is that specialization of variables in a tropical ideal yields another tropical ideal; this plays the role of hyperplane sections in the theory. We also show that elimination theory (projection of varieties) works for tropical ideals as in the classical case. The matroid condition that defines tropical ideals is crucial for these results.


title: Learning Object Detection from Captions via Textual Scene Attributes

abstract: Object detection is a fundamental task in computer vision, requiring large annotated datasets that are difficult to collect, as annotators need to label objects and their bounding boxes. Thus, it is a significant challenge to use cheaper forms of supervision effectively. Recent work has begun to explore image captions as a source for weak supervision, but to date, in the context of object detection, captions have only been used to infer the categories of the objects in the image. In this work, we argue that captions contain much richer information about the image, including attributes of objects and their relations. Namely, the text represents a scene of the image, as described recently in the literature. We present a method that uses the attributes in this "textual scene graph" to train object detectors. We empirically demonstrate that the resulting model achieves state-of-the-art results on several challenging object detection datasets, outperforming recent approaches.


title: Robust Utility Maximization in a Multivariate Financial Market with Stochastic Drift

abstract: We study a utility maximization problem in a financial market with a stochastic drift process, combining a worst-case approach with filtering techniques. Drift processes are difficult to estimate from asset prices, and at the same time optimal strategies in portfolio optimization problems depend crucially on the drift. We approach this problem by setting up a worst-case optimization problem with a time-dependent uncertainty set for the drift. Investors assume that the worst possible drift process with values in the uncertainty set will occur. This leads to local optimization problems, and the resulting optimal strategy needs to be updated continuously in time. We prove a minimax theorem for the local optimization problems and derive the optimal strategy. Further, we show how an ellipsoidal uncertainty set can be defined based on filtering techniques and demonstrate that investors need to choose a robust strategy to be able to profit from additional information.


title: Intersection theorems for triangles

abstract: Given a family of sets on the plane, we say that the family is intersecting if for any two sets from the family their interiors intersect. In this paper, we study intersecting families of triangles with vertices in a given set of points. In particular, we show that if a set $P$ of $n$ points is in convex position, then the largest intersecting family of triangles with vertices in $P$ contains at most $(\frac{1}{4}+o(1))\binom{n}{3}$ triangles.


title: Are cryptocurrencies becoming more interconnected?

abstract: This paper studies the dynamic market linkages among cryptocurrencies during August 2015 - July 2020 and finds a substantial increase in market linkages for both returns and volatilities. We use different methodologies to check the different aspects of market linkages. Financial and regulatory implications are discussed.


title: Friction on incommensurate substrates: Role of anharmonicity and defects

abstract: We present Molecular Dynamics simulations of one- and two-dimensional bead-spring models sliding on incommensurate substrates. We investigate how sliding friction is affected by interaction anharmonicity and structural defects. In their absence, we confirm earlier findings, namely, that at special resonance sliding velocities, friction is maximal. When sliding off-resonance, partially thermalized states are possible, whereby only a small number of vibrational modes becomes excited, but whose kinetic energies are already Maxwell-Boltzmann distributed. Anharmonicity and defects typically destroy partial thermalization, and instead lead to full thermalization, implying much higher friction. For sliders with periodic boundaries, thermalization begins with vibrational modes whose spatial modulation is compatible with the incommensurate lattice. For a disc-shaped slider, modes corresponding to modulations compatible with the slider radius are initially the most dominant. By tuning the mechanical properties of the slider's edge, this effect can be controlled, resulting in significant changes in the sliding distance covered.


title: Restoring Spatially-Heterogeneous Distortions using Mixture of Experts Network

abstract: In recent years, deep learning-based methods have been successfully applied to the image distortion restoration tasks. However, scenarios that assume a single distortion only may not be suitable for many real-world applications. To deal with such cases, some studies have proposed sequentially combined distortions datasets. Viewing in a different point of combining, we introduce a spatially-heterogeneous distortion dataset in which multiple corruptions are applied to the different locations of each image. In addition, we also propose a mixture of experts network to effectively restore a multi-distortion image. Motivated by the multi-task learning, we design our network to have multiple paths that learn both common and distortion-specific representations. Our model is effective for restoring real-world distortions and we experimentally verify that our method outperforms other models designed to manage both single distortion and multiple distortions.


title: Defining the spectral position of a Neumann domain

abstract: A Laplacian eigenfunction on a two-dimensional Riemannian manifold provides a natural partition into Neumann domains (a.k.a. a Morse--Smale complex). This partition is generated by gradient flow lines of the eigenfunction, which bound the so-called Neumann domains. We prove that the Neumann Laplacian defined on a Neumann domain is self-adjoint and has a purely discrete spectrum. In addition, we prove that the restriction of an eigenfunction to any one of its Neumann domains is an eigenfunction of the Neumann Laplacian. By comparison, similar statements about the Dirichlet Laplacian on a nodal domain of an eigenfunction are basic and well-known. The difficulty here is that the boundary of a Neumann domain may have cusps and cracks, so standard results about Sobolev spaces are not available. Another very useful common fact is that the restricted eigenfunction on a nodal domain is the first eigenfunction of the Dirichlet Laplacian. This is no longer true for a Neumann domain. Our results enable the investigation of the resulting spectral position problem for Neumann domains, which is much more involved than its nodal analogue.


title: Quantifying the Performance of a Hybrid Pixel Detector with GaAs:Cr Sensor for Transmission Electron Microscopy

abstract: Hybrid pixel detectors (HPDs) have been shown to be highly effective for diffraction-based and time-resolved studies in transmission electron microscopy, but their performance is limited by the fact that high-energy electrons scatter over long distances in their thick Si sensors. An advantage of HPDs compared to monolithic active pixel sensors (MAPS) is that their sensor does not need to be fabricated from Si. We have compared the performance of the Medipix3 HPD with a Si sensor and with a GaAs:Cr sensor using primary electrons in the energy range of 60 - 300keV. We describe the measurement and calculation of the detectors' modulation transfer function (MTF) and detective quantum efficiency (DQE), which show that the performance of the GaAs:Cr device is markedly superior to that of the Si device for high-energy electrons.


title: Gradient Estimates For $\Delta u + au^{p+1}=0$ And Liouville Theorems

abstract: In this short note, we use a unified method to consider the gradient estimates of the positive solution to the following nonlinear elliptic equation $\Delta u + au^{p+1}=0$ defined on a complete noncompact Riemannian manifold $(M, g)$ where $a > 0$ and $ p <\frac{4}{n}$ or $a < 0$ and $p >0$ are two constants. For the case $a>0$, this improves considerably the previous known results except for the cases $\dim(M)=4$ and supplements the results for the case $\dim(M)\leq 2$. For the case $a<0$ and $p>0$, we also improve considerably the previous related results. When the Ricci curvature of $(M,g)$ is nonnegative, we also obtain a Liouville-type theorem for the above equation.


title: A Bayesian ILC method for CMB B-mode posterior estimation and reconstruction of primordial gravity wave signal

abstract: The Cosmic Microwave Background (CMB) radiation B mode polarization signal contains the unique signature of primordial metric perturbations produced during the inflation. The separation of the weak CMB B-mode signal from strong foreground contamination in observed maps is a complex task, and proposed new generation low noise satellite missions compete with the weak signal level of this gravitational background. In this article, for the first time, we employ a foreground model-independent internal linear combination (ILC) method to reconstruct the CMB B mode signal using simulated observations over large angular scales of the sky of 6 frequency bands of future generation CMB mission Probe of Inflation and Cosmic Origins (PICO). We estimate the joint CMB B mode posterior density following the interleaving Gibbs steps of B mode angular power spectrum and cleaned map samples using the ILC method. We extend and improve the earlier reported Bayesian ILC method to analyze weak CMB B mode reconstruction by introducing noise bias corrections at two stages during the ILC weight estimation. By performing $200$ Monte Carlo simulations of the Bayesian ILC method, we find that our method can reconstruct the CMB signals and the joint posterior density accurately over large angular scales of the sky. We estimate Blackwell-Rao statistics of the marginal density of CMB B mode angular power spectrum and use them to estimate the joint density of scalar to tensor ratio $r$ and a lensing power spectrum amplitude $A^{\textrm{lens}}$. Using $200$ Monte Carlo simulations of the delensing approach, we find that our method can achieve an unbiased detection of the primordial gravitational wave signal $r$ with more than 8$\sigma$ significance for levels of $r \geqslant 0.01$.


title: A Note On Two Fiber Bundles and The Manifestations Of "Shtuka"

abstract: In this note we intend to look at the moduli stacks for global $G$-shtukas from a new perspective. We discuss a unifying interpretation of several moduli spaces (stacks) including moduli of global $G$-shtukas and (a variant of the) moduli of Higgs bundles. We view these spaces (stacks) as different fibers of a family over a scheme (stack) locally of finite type. We discuss (a relative version of) the local model theory for this family. We also consider the Hecke stacks over the moduli stack of $G$-shtukas and discuss the corresponding (motivic) Hecke operations.


title: Spin moir\'{e} engineering of topological magnetism and emergent electromagnetic fields

abstract: A superposition of spin helices can yield topological spin textures, such as skyrmion and hedgehog lattices. Based on the analogy with the moir\'e in optics, we study the magnetic and topological properties of such superpositions in a comprehensive way by modulating the interference pattern continuously. We find that the control of the angles between the superposed helices and the net magnetization yields successive topological transitions associated with pair annihilation of hedgehogs and antihedgehogs. Accordingly, emergent electromagnetic fields, magnetic monopoles and antimonopoles, and Dirac strings arising from the noncoplanar spin textures show systematic evolution. In addition, we also show how the system undergoes the magnetic transitions with dimensional reduction from the three-dimensional hedgehog lattice to a two-dimensional skyrmion lattice or a one-dimensional conical state. The results indicate that the concept of "spin moir\'{e}" provides an efficient way of engineering the emergent electromagnetism and topological nature in magnets.


title: A robustness measure for singular point and index estimation in discretized orientation and vector fields

abstract: The identification of singular points or topological defects in discretized vector fields occurs in diverse areas ranging from the polarization of the cosmic microwave background to liquid crystals to fingerprint recognition and bio-medical imaging. Due to their discrete nature, defects and their topological charge cannot depend continuously on each single vector, but they discontinuously change as soon as a vector changes by more than a threshold. Considering this threshold of admissible change at the level of vectors, we develop a robustness measure for discrete defect estimators. Here, we compare different template paths for defect estimation in discretized vector or orientation fields. Sampling prototypical vector field patterns around defects shows that the robustness increases with the length of template path, but less so in the presence of noise on the vectors. We therefore find an optimal trade-off between resolution and robustness against noise for relatively small templates, except for the "single pixel" defect analysis, which cannot exclude zero robustness. The presented robustness measure paves the way for uncertainty quantification of defects in discretized vector fields.


title: Unification of variational multiscale analysis and Nitsche's method, and a resulting boundary layer fine-scale model

abstract: We show that in the variational multiscale framework, the weak enforcement of essential boundary conditions via Nitsche's method corresponds directly to a particular choice of projection operator. The consistency, symmetry and penalty terms of Nitsche's method all originate from the fine-scale closure dictated by the corresponding scale decomposition. As a result of this formalism, we are able to determine the exact fine-scale contributions in Nitsche-type formulations. In the context of the advection-diffusion equation, we develop a residual-based model that incorporates the non-vanishing fine scales at the Dirichlet boundaries. This results in an additional boundary term with a new model parameter. We then propose a parameter estimation strategy for all parameters involved that is also consistent for higher-order basis functions. We illustrate with numerical experiments that our new augmented model mitigates the overly diffusive behavior that the classical residual-based fine-scale model exhibits in boundary layers at boundaries with weakly enforced essential conditions.


title: Uncovering Feature Interdependencies in High-Noise Environments with Stepwise Lookahead Decision Forests

abstract: Conventionally, random forests are built from "greedy" decision trees which each consider only one split at a time during their construction. The sub-optimality of greedy implementation has been well-known, yet mainstream adoption of more sophisticated tree building algorithms has been lacking. We examine under what circumstances an implementation of less greedy decision trees actually yields outperformance. To this end, a "stepwise lookahead" variation of the random forest algorithm is presented for its ability to better uncover binary feature interdependencies. In contrast to the greedy approach, the decision trees included in this random forest algorithm, each simultaneously consider three split nodes in tiers of depth two. It is demonstrated on synthetic data and financial price time series that the lookahead version significantly outperforms the greedy one when (a) certain non-linear relationships between feature-pairs are present and (b) if the signal-to-noise ratio is particularly low. A long-short trading strategy for copper futures is then backtested by training both greedy and stepwise lookahead random forests to predict the signs of daily price returns. The resulting superior performance of the lookahead algorithm is at least partially explained by the presence of "XOR-like" relationships between long-term and short-term technical indicators. More generally, across all examined datasets, when no such relationships between features are present, performance across random forests is similar. Given its enhanced ability to understand the feature-interdependencies present in complex systems, this lookahead variation is a useful extension to the toolkit of data scientists, in particular for financial machine learning, where conditions (a) and (b) are typically met.


title: Rain-Code Fusion : Code-to-code ConvLSTM Forecasting Spatiotemporal Precipitation

abstract: Recently, flood damage has become a social problem owing to unexperienced weather conditions arising from climate change. An immediate response to heavy rain is important for the mitigation of economic losses and also for rapid recovery. Spatiotemporal precipitation forecasts may enhance the accuracy of dam inflow prediction, more than 6 hours forward for flood damage mitigation. However, the ordinary ConvLSTM has the limitation of predictable range more than 3-timesteps in real-world precipitation forecasting owing to the irreducible bias between target prediction and ground-truth value. This paper proposes a rain-code approach for spatiotemporal precipitation code-to-code forecasting. We propose a novel rainy feature that represents a temporal rainy process using multi-frame fusion for the timestep reduction. We perform rain-code studies with various term ranges based on the standard ConvLSTM. We applied to a dam region within the Japanese rainy term hourly precipitation data, under 2006 to 2019 approximately 127 thousands hours, every year from May to October. We apply the radar analysis hourly data on the central broader region with an area of 136 x 148 km2 . Finally we have provided sensitivity studies between the rain-code size and hourly accuracy within the several forecasting range.


title: Sequential Algorithms and Independent Sets Discovering on Large Sparse Random Graphs

abstract: Computing the size of maximum independent sets is a NP-hard problem for fixed graphs. Characterizing and designing efficient algorithms to estimate this independence number for random graphs are notoriously difficult and still largely open issues. In a companion paper, we showed that a low complexity degree-greedy exploration is actually asymptotically optimal on a large class of sparse random graphs. Encouraged by this result, we present and study two variants of sequential exploration algorithms: static and dynamic degree-aware explorations. We derive hydrodynamic limits for both of them, which in turn allow us to compute the size of the resulting independent set. Whereas the former is simpler to compute, the latter may be used to arbitrarily approximate the degree-greedy algorithm. Both can be implemented in a distributed manner. The corresponding hydrodynamic limits constitute an efficient method to compute or bound the independence number for a large class of sparse random graphs. As an application, we then show how our method may be used to estimate the capacity of a large 802.11-based wireless network. We finally consider further indicators such as the fairness of the resulting configuration, and show how an unexpected trade-off between fairness and capacity can be achieved.


title: First-order Optimization for Superquantile-based Supervised Learning

abstract: Classical supervised learning via empirical risk (or negative log-likelihood) minimization hinges upon the assumption that the testing distribution coincides with the training distribution. This assumption can be challenged in modern applications of machine learning in which learning machines may operate at prediction time with testing data whose distribution departs from the one of the training data. We revisit the superquantile regression method by proposing a first-order optimization algorithm to minimize a superquantile-based learning objective. The proposed algorithm is based on smoothing the superquantile function by infimal convolution. Promising numerical results illustrate the interest of the approach towards safer supervised learning.


title: A String Diagrammatic Axiomatisation of Finite-State Automata

abstract: We develop a fully diagrammatic approach to the theory of finite-state automata, based on reinterpreting their usual state-transition graphical representation as a two-dimensional syntax of string diagrams. Moreover, we provide an equational theory that completely axiomatises language equivalence in this new setting. This theory has two notable features. First, the Kleene star is a derived concept, as it can be decomposed into more primitive algebraic blocks. Second, the proposed axiomatisation is finitary -- a result which is provably impossible to obtain for the one-dimensional syntax of regular expressions.


title: Variations in magnetotelluric response functions caused by Sq currents

abstract: The magnetotelluric (MT) responses of the Earth are biased by spatially heterogeneous source fields. Recently, such biases have been reported even in mid-latitude areas, where localized source currents rarely flow. This study focuses on shifts in the MT responses arising from variations in the focus latitude of the solar quiet (Sq) current. The MT responses at 60 s were calculated by changing the center of the Sq current. Slight variations in the focus latitude cause large shifts in the apparent resistivity and phase. During periods of quiet geomagnetic activity, the center varies within a wide range of 20-45 N, whereas this range narrows when there are disruptions. Even though this study considers only a limited case, the results demonstrate the instability of the impedances in periods of quiet geomagnetic activity and a correlation of the MT responses with the magnitude of the geomagnetic activity. As a consequence, even when a station is located at mid-latitudes, the impedances need to be treated carefully for time-lapse MT soundings, for example, by checking the ionospheric current conditions and their dependence on the geomagnetic activity.


title: Dilated Convolutional Attention Network for Medical Code Assignment from Clinical Text

abstract: Medical code assignment, which predicts medical codes from clinical texts, is a fundamental task of intelligent medical information systems. The emergence of deep models in natural language processing has boosted the development of automatic assignment methods. However, recent advanced neural architectures with flat convolutions or multi-channel feature concatenation ignore the sequential causal constraint within a text sequence and may not learn meaningful clinical text representations, especially for lengthy clinical notes with long-term sequential dependency. This paper proposes a Dilated Convolutional Attention Network (DCAN), integrating dilated convolutions, residual connections, and label attention, for medical code assignment. It adopts dilated convolutions to capture complex medical patterns with a receptive field which increases exponentially with dilation size. Experiments on a real-world clinical dataset empirically show that our model improves the state of the art.


title: A resolution of the Trans-Planckian problem in the R_h=ct universe

abstract: The recent measurement of a cutoff k_min in the fluctuation power spectrum P(k) of the cosmic microwave background may vitiate the possibility that slow-roll inflation can simultaneously solve the horizon problem and account for the formation of structure via the growth of quantum fluctuations in the inflaton field. Instead, we show that k_min may be interpreted more successfully in the R_h=ct cosmology, as the first mode exiting from the Planck scale into the semi-classical Universe shortly after the Big Bang. In so doing, we demonstrate that such a scenario completely avoids the well-known trans-Planckian problem plaguing standard inflationary cosmology.


title: Negative entropy production rates in Drude-Sommerfeld metals

abstract: It is a commonly accepted creed that in typical situations the rate of entropy production is non-negative. We show that this assertion is not entirely correct if a time-dependent, external perturbation is not compensated by a rapid enough decay of the response function. This is demonstrated for three variants of the Drude model to describe electrical conduction in noble metals, namely the classical free electron gas, the Drude-Sommerfeld model, and the Extended Drude-Sommerfeld model. The analysis is concluded with a discussion of potential experimental verifications and ramifications of negative entropy production rates.

categories: cond-mat.stat-mech cond-mat.quant-gas cond-mat.str-el quant-ph


title: Study of the spatial association between an active region jet and a nonthermal type~${\rm III}$ radio burst

abstract: We aim to investigate the spatial location of the source of an active region (AR) jet and its relation with associated nonthermal type~III radio emission. An emission measure (EM) method was used to study the thermodynamic nature of the AR jet. The nonthermal type~{\rm III} radio burst observed at meterwavelength was studied using the Murchison Widefield Array (MWA) radio imaging and spectroscopic data. The local configuration of the magnetic field and the connectivity of the source region of the jet with open magnetic field structures was studied using a nonlinear force-free field (NLFFF) extrapolation and potential field source surface (PFSS) extrapolation respectively. The plane-of-sky velocity of the AR jet was found to be $\sim$136~km/s. The EM analysis confirmed the presence of low temperature 2~MK plasma for the spire, whereas hot plasma, between 5-8 MK, was present at the footpoint region which also showed the presence of Fe~{\sc xviii} emission. A lower limit on the electron number density was found to be 1.4$\times$10$^{8}$ cm$^{-3}$ for the spire and 2.2$\times$10$^{8}$~cm$^{-3}$ for the footpoint. A temporal and spatial correlation between the AR jet and nonthermal type III burst confirmed the presence of open magnetic fields. An NLFFF extrapolation showed that the photospheric footpoints of the null point were anchored at the location of the source brightening of the jet. The spatial location of the radio sources suggests an association with the extrapolated closed and open magnetic fields although strong propagation effects are also present. The multi-scale analysis of the field at local, AR, and solar scales confirms the interlink between different flux bundles involved in the generation of the type III radio signal with flux transferred from a small coronal hole to the periphery of the sunspot via null point reconnection with an emerging structure.


title: Constraining Teleparallel Gravity through Gaussian Processes

abstract: We apply Gaussian processes (GP) in order to impose constraints on teleparallel gravity and its $f(T)$ extensions. We use available $H(z)$ observations from (i) cosmic chronometers data (CC); (ii) Supernova Type Ia (SN) data from the compressed Pantheon release together with the CANDELS and CLASH Multi-Cycle Treasury programs; and (iii) baryonic acoustic oscillation (BAO) datasets from the Sloan Digital Sky Survey. For the involved covariance functions, we consider four widely used choices, namely the square exponential, Cauchy, Mat\'{e}rn and rational quadratic kernels, which are consistent with one another within 1$\sigma$ confidence levels. Specifically, we use the GP approach to reconstruct a model-independent determination of the Hubble constant $H_0$, for each of these kernels and dataset combinations. These analyses are complemented with three recently announced literature values of $H_0$, namely (i) Riess $H_0^{\rm R} = 74.22 \pm 1.82 \,{\rm km\, s}^{-1} {\rm Mpc}^{-1}$; (ii) H0LiCOW Collaboration $H_0^{\rm HW} = 73.3^{+1.7}_{-1.8} \,{\rm km\, s}^{-1} {\rm Mpc}^{-1}$; and (iii) Carnegie-Chicago Hubble Program $H_0^{\rm TRGB} = 69.8 \pm 1.9 \,{\rm km\, s}^{-1} {\rm Mpc}^{-1}$. Additionally, we investigate the transition redshift between the decelerating and accelerating cosmological phases through the GP reconstructed deceleration parameter. Furthermore, we reconstruct the model-independent evolution of the dark energy equation of state, and finally reconstruct the allowed $f(T)$ functions. As a result, the $\Lambda$CDM model lies inside the allowed region at 1$\sigma$ in all the examined kernels and datasets, however a negative slope for $f(T)$ versus $T$ is slightly favored.

categories: gr-qc astro-ph.CO hep-th physics.data-an


title: Positional games on randomly perturbed graphs

abstract: Maker-Breaker games are played on a hypergraph $(X,\mathcal{F})$, where $\mathcal{F} \subseteq 2^X$ denotes the family of winning sets. Both players alternately claim a predefined amount of edges (called bias) from the board $X$, and Maker wins the game if she is able to occupy any winning set $F \in \mathcal{F}$. These games are well studied when played on the complete graph $K_n$ or on a random graph $G_{n,p}$. In this paper we consider Maker-Breaker games played on randomly perturbed graphs instead. These graphs consist of the union of a deterministic graph $G_\alpha$ with minimum degree at least $\alpha n$ and a binomial random graph $G_{n,p}$. Depending on $\alpha$ and Breaker's bias $b$ we determine the order of the threshold probability for winning the Hamiltonicity game and the $k$-connectivity game on $G_{\alpha}\cup G_{n,p}$, and we discuss the $H$-game when $b=1$. Furthermore, we give optimal results for the Waiter-Client versions of all mentioned games.


title: Progress on band structure engineering of twisted bilayer and two-dimensional moir\'e heterostructures

abstract: Artificially constructed van der Waals heterostructures (vdWHs) provide an ideal platform for realizing emerging quantum phenomena in condensed matter physics. Two methods for building vdWHs have been developed: stacking two-dimensional (2D) materials into a bilayer structure with different lattice constants, or with different orientations. The interlayer coupling stemming from commensurate or incommensurate superlattice pattern plays an important role in vdWHs for modulating the band structures and generating new electronic states. In this article, we review a series of novel quantum states discovered in two model vdWH systems -- graphene/hexagonal boron nitride (hBN) hetero-bilayer and twisted bilayer graphene (tBLG), and discuss how the electronic structures are modified by such stacking and twisting. We also provide perspectives for future studies on hetero-bilayer materials, from which an expansion of 2D material phase library is expected.


title: HI study of isolated and pair galaxies: the MIR SFR-M$\star$ sequence

abstract: Using mid-infrared star formation rate and stellar mass indicators in $\textit{WISE}$, we construct and contrast the relation between star formation rate and stellar mass for isolated and paired galaxies. Our samples comprise a selection of AMIGA (isolated galaxies) and pairs of ALFALFA galaxies with HI detections such that we can examine the relationship between HI content (gas fraction, HI deficiency) and galaxy location on the main sequence (MS) in these two contrasting environments. We derive for the first time an HI scaling relation for isolated galaxies using $\textit{WISE}$ stellar masses, and thereby establish a baseline predictor of HI content that can be used to assess the impact of environment on HI content when compared with samples of galaxies in different environments. We use this updated relation to determine the HI deficiency of both our paired and isolated galaxies. Across all the quantities examined as a function of environment in this work (MS location, gas fraction, and HI deficiency), the AMIGA sample of isolated galaxies is found to have the lower dispersion: $\sigma_{\rm{AMIGA}} = 0.37$ versus $\sigma_{\rm{PAIRS}} = 0.55$ on the MS, $\sigma_{\rm{AMIGA}} = 0.44$ versus $\sigma_{\rm{PAIRS}} = 0.54$ in gas fraction, and $\sigma_{\rm{AMIGA}} = 0.28$ versus $\sigma_{\rm{PAIRS}} = 0.34$ in HI deficiency. We also note fewer isolated quiescent galaxies, 3 (0.6$\%$), compared to 12 (2.3$\%$) quiescent pair members. Our results suggest the differences in scatter measured between our samples are environment driven. Galaxies in isolation behave relatively predictably, and galaxies in more densely populated environments adopt a more stochastic behaviour, across a broad range of quantities.


title: Classical Density Functional Theory applied to the solid state

abstract: The standard model of classical Density Functional Theory for pair potentials consists of a hard-sphere functional plus a mean-field term accounting for long ranged attraction. However, most implementations using sophisticated Fundamental Measure hard-sphere functionals suffer from potential numerical instabilities either due to possible instabilities in the functionals themselves or due to implementations that mix real- and Fourier-space components inconsistently. Here, we present a new implementation based on a demonstrably stable hard-sphere functional that is implemented in a completely consistent manner. The present work does not depend on approximate spherical integration schemes and so is much more robust than previous algorithms. The methods are illustrated by calculating phase diagrams for the solid state using the standard Lennard-Jones potential as well as a new class of potentials recently proposed by Wang et al (Phys. Chem. Chem. Phys. 22, 10624 (2020)). The latter span the range from potentials for small molecules to those appropriate to colloidal systems simply by varying a parameter. We verify that cDFT is able to semi-quantitatively reproduce the phase diagram in all cases. We also show that for these problems computationally cheap Gaussian approximations are nearly as good as full minimization based on finite differences.


title: Pointwise Universal Gysin formulae and Applications towards Griffiths' conjecture

abstract: Let $X$ be a complex manifold, $(E,h)\to X$ be a rank $r$ holomorphic hermitian vector bundle, and $\rho$ be a sequence of dimensions $0 = \rho_0 < \rho_1 < \cdots < \rho_m = r$. Let $Q_{\rho,j}$, $j=1,\dots,m$, be the tautological line bundles over the (possibly incomplete) flag bundle $\mathbb{F}_{\rho}(E) \to X$ associated to $\rho$, endowed with the natural metrics induced by that of $E$, with Chern curvatures $\Xi_{\rho,j}$. We show that the universal Gysin formula \textsl{\`{a} la} Darondeau--Pragacz for the push-forward of a homogeneous polynomial in the Chern classes of the $Q_{\rho,j}$'s also hold pointwise at the level of the Chern forms $\Xi_{\rho,j}$ in this hermitianized situation. As an application, we show the positivity of several polynomials in the Chern forms of a Griffiths (semi)positive vector bundle not previously known, thus giving some new evidences towards a conjecture by Griffiths, which in turn can be seen as a pointwise hermitianized version of the Fulton--Lazarsfeld Theorem on numerically positive polynomials for ample vector bundles.


title: EWS-GCN: Edge Weight-Shared Graph Convolutional Network for Transactional Banking Data

abstract: In this paper, we discuss how modern deep learning approaches can be applied to the credit scoring of bank clients. We show that information about connections between clients based on money transfers between them allows us to significantly improve the quality of credit scoring compared to the approaches using information about the target client solely. As a final solution, we develop a new graph neural network model EWS-GCN that combines ideas of graph convolutional and recurrent neural networks via attention mechanism. The resulting model allows for robust training and efficient processing of large-scale data. We also demonstrate that our model outperforms the state-of-the-art graph neural networks achieving excellent results


title: Hidden Markov Models for Pipeline Damage Detection Using Piezoelectric Transducers

abstract: Oil and gas pipeline leakages lead to not only enormous economic loss but also environmental disasters. How to detect the pipeline damages including leakages and cracks has attracted much research attention. One of the promising leakage detection method is to use lead zirconate titanate (PZT) transducers to detect the negative pressure wave when leakage occurs. PZT transducers can generate and detect guided stress waves for crack detection also. However, the negative pressure waves or guided stress waves may not be easily detected with environmental interference, e.g., the oil and gas pipelines in offshore environment. In this paper, a Gaussian mixture model based hidden Markov model (GMM-HMM) method is proposed to detect the pipeline leakage and crack depth in changing environment and time-varying operational conditions. Leakages in different sections or crack depths are considered as different states in hidden Markov models (HMM). Laboratory experiments show that the GMM-HMM method can recognize the crack depth and leakage of pipeline such as whether there is a leakage, where the leakage is.


title: Communication-Optimal Parallel Standard and Karatsuba Integer Multiplication in the Distributed Memory Model

abstract: We present COPSIM a parallel implementation of standard integer multiplication for the distributed memory setting, and COPK a parallel implementation of Karatsuba's fast integer multiplication algorithm for a distributed memory setting. When using $\mathcal{P}$ processors, each equipped with a local non-shared memory, to compute the product of tho $n$-digits integer numbers, under mild conditions, our algorithms achieve optimal speedup of the computational time. That is, $\mathcal{O}\left(n^2/\mathcal{P}\right)$ for COPSIM, and $\mathcal{O}\left(n^{\log_2 3}/\mathcal{P}\right)$ for COPK. The total amount of memory required across the processors is $\mathcal{O}\left(n\right)$, that is, within a constant factor of the minimum space required to store the input values. We rigorously analyze the Input/Output (I/O) cost of the proposed algorithms. We show that their bandwidth cost (i.e., the number of memory words sent or received by at least one processors) matches asymptotically corresponding known I/O lower bounds, and their latency (i.e., the number of messages sent or received in the algorithm's critical execution path) is asymptotically within a multiplicative factor $\mathcal{O}\left(\log^2_2 \mathcal{P}\right)$ of the corresponding known I/O lower bounds. Hence, our algorithms are asymptotically optimal with respect to the bandwidth cost and almost asymptotically optimal with respect to the latency cost.


title: Final state Sommerfeld effect on dark matter relic abundance

abstract: If the annihilation products of dark matter (DM) are non-relativistic and if there is some long-range force between them, there can be Sommerfeld effect for the final state particles. We study this effect on DM relic abundance in the thermal freeze-out scenario. As a proof of concept, we consider the case of a DM pair annihilation into a final state pair, assuming that the mutual interactions between the two final state particles give rise to a Coulomb-like potential, and that the masses of the initial and final state particles are similar, so that both the initial and final state particles are non-relativistic. The size of the final state Sommerfeld (FSS) effect depends on the strength of the potential, as well as on the mass ratio of the final and initial state particles. We find that the impact of the FSS effect on DM relic abundance can be significant, and an electroweak sized long-range interaction is large enough to make a correction well beyond the observational accuracy. Another feature of the FSS effect is that it could be suppressed when its time scale is longer than the lifetime of the final state particles. As a corollary, we also study in the DM coannihilation scenario where the initial state Sommerfeld effect between two coannihilators could be reduced due to their instability, which may need to be taken into account for an accurate calculation of the DM relic abundance.


title: Tunable critical field in Rashba superconductor thin-films

abstract: The upper critical field in type II superconductors is limited by the Pauli paramagnetic limit. In superconductors with strong Rashba spin-orbit coupling this limit can be overcome by forming a helical state. Here we quantitatively study the magnetic field-temperature phase diagram of finite-size superconductors with Rashba spin-orbit coupling. We discuss the effect of finite size and shape anisotropy. We demonstrate that the critical field is controllable by intrinsic parameters such as spin-orbit coupling strength and tunable parameters such as sample geometry and applied field direction. Our study opens new avenues for the design of superconducting spin-valves.


title: The Role of Isomorphism Classes in Multi-Relational Datasets

abstract: Multi-interaction systems abound in nature, from colloidal suspensions to gene regulatory circuits. These systems can produce complex dynamics and graph neural networks have been proposed as a method to extract underlying interactions and predict how systems will evolve. The current training and evaluation procedures for these models through the use of synthetic multi-relational datasets however are agnostic to interaction network isomorphism classes, which produce identical dynamics up to initial conditions. We extensively analyse how isomorphism class awareness affects these models, focusing on neural relational inference (NRI) models, which are unique in explicitly inferring interactions to predict dynamics in the unsupervised setting. Specifically, we demonstrate that isomorphism leakage overestimates performance in multi-relational inference and that sampling biases present in the multi-interaction network generation process can impair generalisation. To remedy this, we propose isomorphism-aware synthetic benchmarks for model evaluation. We use these benchmarks to test generalisation abilities and demonstrate the existence of a threshold sampling frequency of isomorphism classes for successful learning. In addition, we demonstrate that isomorphism classes can be utilised through a simple prioritisation scheme to improve model performance, stability during training and reduce training time.


title: Anomalous conductivity of two-dimensional Dirac electrons in organic conductor under pressures

abstract: The electric conductivity of Dirac electrons in the organic conductor $\alpha$-(BEDT-TTF)$_2$I$_3$ [BEDT-TTF = bis(ethylenedithio)tetrathiafulvalene] under pressure has been examined using a two-dimensional tight-binding (TB) model with both impurity and electron--phonon (e--p) scatterings. We study an anomalous temperature dependence of the conductivity, which shows a crossover from $\sigma_{x} < \sigma_{y}$ at low temperatures [region (I)] to $\sigma_{x} > \sigma_{y}$ at high temperatures [region (II)]. $\sigma_y$ and $\sigma_x$ are diagonal conductivities parallel and perpendicular to a stacking axis of molecules, respectively. The effect of Dirac cone tilting is dominant in region (I), whereas the anisotropy of the velocity of the Dirac cone is dominant in region (II). Such behavior is further examined by calculating the deviation of principal axes due to the off-diagonal conductivity $\sigma_{xy}$ and a nearly constant conductivity at high temperatures due to the e--p scattering, which is the extension of the previous result of the simple two-band model [Phys. Rev. B {\bf 98},161205 (2018)]. The relevance to experiments of organic conductors is discussed.


title: Philosophy of Natural Numbers

abstract: We discuss an extension of classical combinatorics theory to the case of spatially distributed objects.


title: Machine Learning and Computational Mathematics

abstract: Neural network-based machine learning is capable of approximating functions in very high dimension with unprecedented efficiency and accuracy. This has opened up many exciting new possibilities, not just in traditional areas of artificial intelligence, but also in scientific computing and computational science. At the same time, machine learning has also acquired the reputation of being a set of "black box" type of tricks, without fundamental principles. This has been a real obstacle for making further progress in machine learning. In this article, we try to address the following two very important questions: (1) How machine learning has already impacted and will further impact computational mathematics, scientific computing and computational science? (2) How computational mathematics, particularly numerical analysis, {can} impact machine learning? We describe some of the most important progress that has been made on these issues. Our hope is to put things into a perspective that will help to integrate machine learning with computational mathematics.


title: A multi-scale DNN algorithm for nonlinear elliptic equations with multiple scales

abstract: Algorithms based on deep neural networks (DNNs) have attracted increasing attention from the scientific computing community. DNN based algorithms are easy to implement, natural for nonlinear problems, and have shown great potential to overcome the curse of dimensionality. In this work, we utilize the multi-scale DNN-based algorithm (MscaleDNN) proposed by Liu, Cai and Xu (2020) to solve multi-scale elliptic problems with possible nonlinearity, for example, the p-Laplacian problem. We improve the MscaleDNN algorithm by a smooth and localized activation function. Several numerical examples of multi-scale elliptic problems with separable or non-separable scales in low-dimensional and high-dimensional Euclidean spaces are used to demonstrate the effectiveness and accuracy of the MscaleDNN numerical scheme.


title: $(t, n)$ Threshold $d$-level quantum secret sharing based on quantum Fourier transformation

abstract: Quantum secret sharing (QSS) is an important branch of secure multiparty quantum computation. Several schemes for $(n, n)$ threshold QSS based on quantum Fourier transformation (QFT) have been proposed. Inspired by the flexibility of $(t, n)$ threshold schemes, Song {\it et al.} (Scientific Reports, 2017) have proposed a $(t, n)$ threshold QSS utilizing $QFT$. Later, Kao and Hwang (arXiv:1803.00216) have identified a loophole in the scheme but have not suggested any remedy. In this present study, we have proposed a $(t, n)$ threshold QSS scheme to share a $d$ dimensional classical secret. This scheme can be implemented using local operations (such as $QFT$, generalized Pauli operators and local measurement) and classical communication. Security of the proposed scheme is described against outsider and participants' eavesdropping.


title: Effect of mantle oxidation state and escape upon the evolution of Earth's magma ocean atmosphere

abstract: The magma ocean period was a critical phase determining how Earth atmosphere developed into habitability. However there are major uncertainties in the role of key processes such as outgassing from the planetary interior and escape of species to space that play a major role in determining the atmosphere of early Earth. We investigate the influence of outgassing of various species and escape of H$_2$ for different mantle redox states upon the composition and evolution of the atmosphere for the magma ocean period. We include an important new atmosphere-interior coupling mechanism namely the redox evolution of the mantle which strongly affects the outgassing of species. We simulate the volatile outgassing and chemical speciation at the surface for various redox states of the mantle by employing a C-H-O based chemical speciation model combined with an interior outgassing model. We then apply a line-by-line radiative transfer model to study the remote appearance of the planet in terms of the infrared emission and transmission. Finally, we use a parameterized diffusion-limited and XUV energy-driven atmospheric escape model to calculate the loss of H$_2$ to space. We have simulated the thermal emission and transmission spectra for reduced or oxidized atmospheres present during the magma ocean period of Earth. Reduced or thin atmospheres consisting of H$_2$ in abundance emit more radiation to space and have larger effective height as compared to oxidized or thick atmospheres which are abundant in H$_2$O and CO$_2$. We obtain the outgassing rates of H2 from the mantle into the atmosphere to be a factor of ten times larger than the rates of diffusion-limited escape to space. Our work presents useful insight into the development of Earth atmosphere during the magma ocean period as well as input to guide future studies discussing exoplanetary interior compositions.


title: Accelerating Sparse Matrix-Matrix Multiplication with GPU Tensor Cores

abstract: Sparse general matrix-matrix multiplication (spGEMM) is an essential component in many scientific and data analytics applications. However, the sparsity pattern of the input matrices and the interaction of their patterns make spGEMM challenging. Modern GPUs include Tensor Core Units (TCUs), which specialize in dense matrix multiplication. Our aim is to re-purpose TCUs for sparse matrices. The key idea of our spGEMM algorithm, tSparse, is to multiply sparse rectangular blocks using the mixed precision mode of TCUs. tSparse partitions the input matrices into tiles and operates only on tiles which contain one or more elements. It creates a task list of the tiles, and performs matrix multiplication of these tiles using TCUs. To the best of our knowledge, this is the first time that TCUs are used in the context of spGEMM. We show that spGEMM, with our tiling approach, benefits from TCUs. Our approach significantly improves the performance of spGEMM in comparison to cuSPARSE, CUSP, RMerge2, Nsparse, AC-SpGEMM and spECK.


title: Electronic structure and finite temperature magnetism of yttrium iron garnet

abstract: Yttrium iron garnet is a complex ferrimagnetic insulator with 20 magnon modes which is used extensively in fundamental experimental studies of magnetisation dynamics. As a transition metal oxide with moderate gap (2.8 eV), yttrium iron garnet requires a careful treatment of electronic correlation. We have applied quasiparticle self-consistent GW to provide a fully ab initio description of the electronic structure and resulting magnetic properties, including the parameterisation of a Heisenberg model for magnetic exchange interactions. Subsequent spin dynamical modelling with quantum statistics extends our description to the magnon spectrum and thermodynamic properties such as the Curie temperature, finding favourable agreement with experimental measurements. This work provides a snapshot of the state-of-the art in modelling of complex magnetic insulators.


title: Covid-19 Modeling towards socioeconomic and health data from New South Wales (NSW) -- Australia: An approach via Geospatial Analysis and Geographically Weighted Poisson Regression (GWPR)

abstract: An integrated approach of spatial data analysis and Geographically Weighted Poisson Regression (GWPR) along with global regression techniques are used in this study. This approach aims to model relationships between dependent variable Covid-19 and independent variables from socioeconomic and pre-existing health conditions within the local government area (LGA) in New South Wales (NSW)-Australia. Based on geospatial data analysis and a step-by-step procedure in building both global and GWPR models, four (4) independent variables are finally selected to investigate relationships between dependent and independent variables at the local scale. The GWPR model's results with the Goodness-of-Fit (R2) range between 45-73% exhibit positive relationships between Covid-19 and the total population, the cancers, and the people with ages between 60 and 85 in most of the NSW state. Meanwhile, a negative relationship is observed between Covid-19 and the ischaemic heart disease; however, the estimated coefficients for this relationship are very low and close to zero; hence further investigation, including assessment from a different perspective, is necessary for validation. In conclusion, the model suggests that the relationships between the dependent variable and independent variables are nonstationary. Therefore, GWPR model calibration plays a vital role in geographic modelling at the local scale.


title: Superconducting Dirac point in proximetized graphene

abstract: Two-dimensional (2D) materials, composed of single atomic layers, have attracted vast research interest since the breakthrough discovery of graphene. One major benefit of such systems is the simple ability to tune the chemical potential by back-gating, in-principle enabling to vary the Fermi level through the charge neutrality point, thus tuning between electron and hole doping. For 2D Superconductors, this means that one may potentially achieve the strongly-coupled superconducting regime described by Bose Einstein Condensation physics of small bosonic tightly bound electron pairs. Furthermore, it should be possible to access both electron and hole based superconductivity in a single system. However, in most 2D materials, an insulating gap opens up around the charge neutrality point, thus preventing approach to this regime. Graphene is unique in this sense since it is a true semi-metal in which the un-gapped Dirac point is protected by the symmetries. In this work we show that single layer graphene, in which superconducting pairing is induced by proximity to regions of a low density superconductor, can be tuned from hole to electron superconductivity through the strong coupling regime. We study, both experimentally and theoretically, the vicinity of this "Superconducting Dirac point" and find an unusual situation where reflections at interfaces between normal and superconducting regions within the graphene, suppress the conductance and, at the same time, Andreev reflections maintain a large phase breaking length. In addition, the Fermi level can be adjusted so that the momentum in the normal and superconducting regimes perfectly match giving rise to ideal Andreev reflection processes.


title: Detection of a high-confidence quasi-periodic oscillation in radio light curve of the high redshift FSRQ PKS J0805-0111

abstract: In this work, we have searched quasi-periodic oscillations (QPOs) in the 15 GHz light curve of the FSRQ PKS J0805-0111 monitored by the Owens Valley Radio Observatory (OVRO) 40 m telescope during the period from January 9,2008 to May 9,2019, using the weighted wavelet Z-transform (WWZ) and the Lomb-Scargle Periodogram (LSP) techniques. This is the first time to search for periodic radio signal in the FSRQ PKS J0805-0111 by these two methods. All two methods consistently reveal a repeating signal with a periodicity of 3.38$\pm$ 0.8 years ($>$99.7\% confidence level). In order to determining the significance of the periods, the false alarm probability method was used, and a large number of Monte Carlo simulations were performed. As possible explanations, we discuss a number of scenarios including the thermal instability of thin disks scenario, the spiral jet scenario and the binary supermassive black hole scenario, we expected that the binary black hole scenario, where the QPO is caused by the precession of the binary black holes, is the most likely explanation. FSRQ PKS J0805-0111 thus could be a good binary black hole candidate. In the binary black hole scenario, the distance between the primary black hole and the secondary black hole is about $1.71\times 10^{16}$ cm.


title: Disruption in the Chinese E-Commerce During COVID-19

abstract: The recent outbreak of the novel coronavirus (COVID-19) has infected millions of citizens worldwide and claimed many lives. This paper examines its impact on the Chinese e-commerce market by analyzing behavioral changes seen from a large online shopping platform. We first conduct a time series analysis to identify product categories that faced the most extensive disruptions. The time-lagged analysis shows that behavioral patterns seen in shopping actions are highly responsive to epidemic development. Based on these findings, we present a consumer demand prediction method by encompassing the epidemic statistics and behavioral features for COVID-19 related products. Experiment results demonstrate that our predictions outperform existing baselines and further extend to the long-term and province-level forecasts. We discuss how our market analysis and prediction can help better prepare for future pandemics by gaining an extra time to launch preventive steps.


title: Improving Generalization of Deep Fault Detection Models in the Presence of Mislabeled Data

abstract: Mislabeled samples are ubiquitous in real-world datasets as rule-based or expert labeling is usually based on incorrect assumptions or subject to biased opinions. Neural networks can "memorize" these mislabeled samples and, as a result, exhibit poor generalization. This poses a critical issue in fault detection applications, where not only the training but also the validation datasets are prone to contain mislabeled samples. In this work, we propose a novel two-step framework for robust training with label noise. In the first step, we identify outliers (including the mislabeled samples) based on the update in the hypothesis space. In the second step, we propose different approaches to modifying the training data based on the identified outliers and a data augmentation technique. Contrary to previous approaches, we aim at finding a robust solution that is suitable for real-world applications, such as fault detection, where no clean, "noise-free" validation dataset is available. Under an approximate assumption about the upper limit of the label noise, we significantly improve the generalization ability of the model trained under massive label noise.


title: VMO spaces associated with Neumann Laplacian

abstract: In this paper, we establish several different characterizations of the vanishing mean oscillation space associated with Neumann Laplacian $\Delta_N$, written ${\rm VMO}_{\Delta_N}(\mathbb{R}^n)$. We first describe it with the classical ${\rm VMO}(\mathbb{R}^n)$ and certain ${\rm VMO}$ on the half-spaces. Then we demonstrate that ${\rm VMO}_{\Delta_N}(\mathbb{R}^n)$ is actually ${\rm BMO}_{\Delta_N}(\mathbb{R}^n)$-closure of the space of the smooth functions with compact supports. Beyond that, it can be characterized in terms of compact commutators of Riesz transforms and fractional integral operators associated to the Neumann Laplacian. Additionally, by means of the functional analysis, we obtain the duality between certain ${\rm VMO}$ and the corresponding Hardy spaces on the half-spaces. Finally, we present an useful approximation for ${\rm BMO}$ functions on the space of homogeneous type, which can be applied to our argument and otherwhere.


title: ATLASGAL -- Relationship between dense star forming clumps and interstellar masers

abstract: We have used catalogues from several Galactic plane surveys and dedicated observations to investigate the relationship between various maser species and Galactic star forming clumps, as identified by the ATLASGAL survey. The maser transitions of interest are the 6.7 & 12.2 GHz methanol masers, 22.2 GHz water masers, and the masers emitting in the four ground-state hyperfine structure transitions of hydroxyl. We find clump association rates for the water, hydroxyl and methanol masers to be 56, 39 and 82 per cent respectively, within the Galactic longitude range of 60{\deg} > $l$ > -60{\deg}. We investigate the differences in physical parameters between maser associated clumps and the full ATLASGAL sample, and find that clumps coincident with maser emission are more compact with increased densities and luminosities. However, we find the physical conditions within the clumps are similar for the different maser species. A volume density threshold of $n$(H$_{2}$) > 10$^{4.1}$ cm$^{-3}$ for the 6.7 GHz methanol maser found in our previous study is shown to be consistent across for all maser species investigated. We find limits that are required for the production of maser emission to be 500 L$_{\odot}$ and 6 M$_{\odot}$ respectively. The evolutionary phase of maser associated clumps is investigated using the L/M ratio of clumps coincident with maser emission, and these have similar L/M ranges (~10$^{0.2}$ - 10$^{2.7}$ L$_{\odot}$/M$_{\odot}$) regardless of the associated transitions. This implies that the conditions required for the production of maser emission only occur during a relatively narrow period during a star's evolution. Lower limits of the statistical lifetimes for each maser species are derived, ranging from ~0.4 - 2 x 10$^{4}$ yrs and are in good agreement with the "straw man" evolutionary model previously presented.


title: Magnetic (quasi-)modular forms

abstract: A (folklore?) conjecture states that no holomorphic modular form $F(\tau)=\sum_{n=1}^\infty a_nq^n\in q\mathbb Z[[q]]$ exists, where $q=e^{2\pi i\tau}$, such that its anti-derivative $\sum_{n=1}^\infty a_nq^n/n$ has integral coefficients in the $q$-expansion. A recent observation of Broadhurst and Zudilin, rigorously accomplished by Li and Neururer, led to examples of meromorphic modular forms possessing the integrality property. In this note we investigate the arithmetic phenomenon from a systematic perspective and discuss related transcendental extensions of the differentially closed ring of quasi-modular forms.

categories: math.NT hep-ph math-ph math.AC math.CO math.MP


title: Concurrent Neural Network : A model of competition between times series

abstract: Competition between times series often arises in sales prediction, when similar products are on sale on a marketplace. This article provides a model of the presence of cannibalization between times series. This model creates a "competitiveness" function that depends on external features such as price and margin. It also provides a theoretical guaranty on the error of the model under some reasonable conditions, and implement this model using a neural network to compute this competitiveness function. This implementation outperforms other traditional time series methods and classical neural networks for market share prediction on a real-world data set.


title: Real-Time Detection of Unmodelled Gravitational-Wave Transients Using Convolutional Neural Networks

abstract: Convolutional Neural Networks (CNNs) have demonstrated potential for the real-time analysis of data from gravitational-wave detector networks for the specific case of signals from coalescing compact-object binaries such as black-hole binaries. Unfortunately, training these CNNs requires a precise model of the target signal; they are therefore not applicable to a wide class of potential gravitational-wave sources, such as core-collapse supernovae and long gamma-ray bursts, where unknown physics or computational limitations prevent the development of comprehensive signal models. We demonstrate for the first time a CNN with the ability to detect generic signals -- those without a precise model -- with sensitivity across a wide parameter space. Our CNN has a novel structure that uses not only the network strain data but also the Pearson cross-correlation between detectors to distinguish correlated gravitational-wave signals from uncorrelated noise transients. We demonstrate the efficacy of our CNN using data from the second LIGO-Virgo observing run, and show that it has sensitivity comparable to that of the "gold-standard" transient searches currently used by LIGO-Virgo, at extremely low (order of 1 second) latency and using only a fraction of the computing power required by existing searches, allowing our models the possibility of true real-time detection of gravitational-wave transients associated with gamma-ray bursts, core-collapse supernovae, and other relativistic astrophysical phenomena.


title: Impact of dynamical dark energy on the neutron star equilibrium

abstract: We study the density distribution of the minimally-coupled scalar field dark energy inside a neutron star. The dark energy is considered in the hydrodynamical representation as a perfect fluid with three parameters (background density, equation of state, and effective sound speed). The neutron star matter is modeled with three unified equations of state, developed by the Brussels-Montreal group. With the calculated density distribution of the dark energy inside a neutron star (and its dependence on the dark energy parameters) we investigate how its presence impacts the macroscopic characteristics and the value of the mass limit for neutron stars. From this impact we derive the possible constrains on the effective speed of sound of dark energy with the help of maximal known masses of observed neutron stars. In this approach, we have found, that the squared effective speed of sound can not be smaller than $\sim 10^{-2}$ in units of squared speed of light.


title: A group-theorist's perspective on symmetry groups in physics

abstract: There are many Lie groups used in physics, including the Lorentz group of special relativity, the spin groups (relativistic and non-relativistic) and the gauge groups of quantum electrodynamics and the weak and strong nuclear forces. Various grand unified theories use larger Lie groups in different attempts to unify some of these groups into something more fundamental. There are also a number of finite symmetry groups that are related to the finite number of distinct elementary particle types. I offer a group-theorist's perspective on these groups, and suggest some ways in which a deeper use of group theory might in principle be useful. These suggestions include a number of options that seem not to be under active investigation at present. I leave open the question of whether they can be implemented in physical theories.


title: Detecting optical transients using artificial neural networks and reference images from different surveys

abstract: To search for optical counterparts to gravitational waves, it is crucial to develop an efficient follow-up method that allows for both a quick telescopic scan of the event localization region and search through the resulting image data for plausible optical transients. We present a method to detect these transients based on an artificial neural network. We describe the architecture of two networks capable of comparing images of the same part of the sky taken by different telescopes. One image corresponds to the epoch in which a potential transient could exist; the other is a reference image of an earlier epoch. We use data obtained by the Dr. Cristina V. Torres Memorial Astronomical Observatory and archival reference images from the Sloan Digital Sky Survey. We trained a convolutional neural network and a dense layer network on simulated source samples and tested the trained networks on samples created from real image data. Autonomous detection methods replace the standard process of detecting transients, which is normally achieved by source extraction of a difference image followed by human inspection of the detected candidates. Replacing the human inspection component with an entirely autonomous method would allow for a rapid and automatic follow-up of interesting targets of opportunity. The method will be further tested on telescopes participating in the Transient Optical Robotic Observatory of the South Collaboration.


title: Online Sparse Sliced Inverse Regression

abstract: Due to the demand for tackling the problem of streaming data with high dimensional covariates, we propose an online sparse sliced inverse regression (OSSIR) method for online sufficient dimension reduction. The existing online sufficient dimension reduction methods focus on the case when the dimension $p$ is small. In this article, we show that our method can achieve better statistical accuracy and computation speed when the dimension $p$ is large. There are two important steps in our method, one is to extend the online principal component analysis to iteratively obtain the eigenvalues and eigenvectors of the kernel matrix, the other is to use the truncated gradient to achieve online $L_{1}$ regularization. We also analyze the convergence of the extended Candid covariance-free incremental PCA(CCIPCA) and our method. By comparing several existing methods in the simulations and real data applications, we demonstrate the effectiveness and efficiency of our method.


title: The cooling-down central star of the planetary nebula SwSt\,1: a late thermal pulse in a massive post-AGB star?

abstract: SwSt 1 (PN G001.5-06.7) is a bright and compact planetary nebula containing a late [WC]-type central star. Previous studies suggested that the nebular and stellar lines are slowly changing with time. We studied new and archival optical and ultraviolet spectra of the object. The [OIII] 4959 and 5007 A to $\mathrm{H}\beta$ line flux ratios decreased between about 1976 and 1997/2015. The stellar spectrum also shows changes between these epochs. We modeled the stellar and nebular spectra observed at different epochs. The analyses indicate a drop of the stellar temperature from about 42 kK to 40.5 kK between 1976 and 1993. We do not detect significant changes between 1993 and 2015. The observations show that the star performed a loop in the H-R diagram. This is possible when a shell source is activated during its post-AGB evolution. We infer that a late thermal pulse (LTP) experienced by a massive post-AGB star can explain the evolution of the central star. Such a star does not expand significantly as the result of the LTP and does not became a born-again red giant. However, the released energy can remove the tiny H envelope of the star.


title: The pursuit of stability in halide perovskites: the monovalent cation and the key for surface and bulk self-repair

abstract: We find significant differences between degradation and healing at the surface or in the bulk for each of the different APbBr3 single crystals (A=CH3NH3+, methylammonium (MA); HC(NH2)2+, formamidinium (FA); and cesium, Cs+). Using 1- and 2-photon microscopy and photobleaching we conclude that kinetics dominate the surface, and thermodynamics the bulk stability. Fluorescence-lifetime imaging microscopy, as well as results from several other methods, relate the (damaged) state of the halide perovskite (HaP) after photobleaching to its modified optical and electronic properties. The A cation type strongly influences both the kinetics and the thermodynamics of recovery and degradation: FA heals best the bulk material with faster self-healing; Cs+ protects the surface best, being the least volatile of the A cations and possibly through O-passivation; MA passivates defects via methylamine from photo-dissociation, which binds to Pb2+. DFT simulations not only provide insight into the latter conclusion, but also show the importance and stability of the Br3- defect. These results rationalize the use of mixed A-cation materials for optimizing both solar cell stability and overall performance of HaP-based devices, and provide a basis for designing new HaP variants.


title: CuTe chains on Cu(111) by deposition of 1/3 ML Te: atomic and electronic structure

abstract: The surface atomic and electronic structure after deposition of 1/3 monolayer (ML) Te on Cu(111) was determined using a combination of low-energy electron diffraction (LEED), scanning tunneling microscopy and spectroscopy (STM/STS), angle-resolved single and two-photon photoelectron spectroscopy (ARPES /AR-2PPE) and density functional theory (DFT) calculations. Contrary to the current state in literature Te does not create a two-dimensional surface alloy but forms Cu$_2$Te$_2$ adsorbate chains in a $\left(2\sqrt{3} \times \sqrt{3}\right)\textrm{R30}^\circ$ superstructure. We establish this by a high-precision LEED-IV structural analysis with Pendry $R$ factor of $R = 0.099$ and corroborating DFT and STM results. The electronic structure of the surface phase is dominated by an anisotropic downward dispersing state at the Fermi energy $E_F$ and a more isotropic upward dispersing unoccupied state at $E-E_F = + 1.43\,\textrm{eV}$. Both states coexist with bulk states of the projected band structure and are therefore surface resonances.


title: X-ray Spectra and Light Curves of Cooling Novae and a Nova-Like

abstract: We present X-ray observations of novae V2491 Cyg and KT Eri about 9 years post-outburst, of the dwarf nova and post-nova candidate EY Cyg, and of a VY Scl variable. The first three objects were observed with XMM-Newton, KT Eri also with the Chandra ACIS-S camera, V794 Aql with the Chandra ACIS-S camera and High Energy Transmission Gratings. The two recent novae, similar in outburst amplitude and light curve, appear very different at quiescence. Assuming half of the gravitational energy is irradiated in X-rays, V2491 Cyg is accreting at $\dot{m}=1.4\times10^{-9}-10^{-8}M_\odot/yr$, while for KT Eri, $\dot{m}<2\times10^{-10}M_\odot/yr$. V2491 Cyg shows signatures of a magnetized WD, specifically of an intermediate polar. A periodicity of ~39 minutes, detected in outburst, was still measured and is likely due to WD rotation. EY Cyg is accreting at $\dot{m}\sim1.8\times10^{-11}M_\odot/yr$, one magnitude lower than KT Eri, consistently with its U Gem outburst behavior and its quiescent UV flux. The X-rays are modulated with the orbital period, despite the system's low inclination, probably due to the X-ray flux of the secondary. A period of ~81 minutes is also detected, suggesting that it may also be an intermediate polar. V794 Aql had low X-ray luminosity during an optically high state, about the same level as in a recent optically low state. Thus, we find no clear correlation between optical and X-ray luminosity: the accretion rate seems unstable and variable. The very hard X-ray spectrum indicates a massive WD.