title: Spectral Modeling of Charge Exchange in the Central Region of M51

abstract: Charge exchange (CX) emission reveals the significant interaction between neutral and ionized interstellar medium (ISM) components of the dense, multiphase, circumnuclear region of a galaxy. We use a model including a thermal and a CX components to describe the high-resolution XMM-Newton/RGS spectrum of the diffuse emission in the central region of M51. Representative signatures of CX emission -- especially the prominent OVII forbidden line and the excess emission in the OVIII Ly$\gamma$ lines -- can be well explained by the model. Combined with the Chandra images in the OVIII and the OVII bands, we find the soft X-ray emission is dominated by the jet-driven outflow and its interaction with the ambient neutral material. The jet-driven outflow itself is likely a thermal plasma of $\sim 0.59$ keV, with mostly sub-solar abundances. It runs into the ambient neutral gas, and produces significant CX emission that accounts for one-fifth of the diffuse X-ray emission in the 7--28 {\AA} band. The effective interface area in the CX process is one order of magnitude greater than the geometrical surface area of the jet-driven outflow. The tenuous outflow driven by the nuclear star formation may also contribute a small portion to both the diffuse thermal and CX emission. The photoionization by the active galactic nuclei (AGNs) and the resonance scattering by the hot gas itself are disfavored, though the effects from past AGN events may not be ruled out.


title: Geometry-aware Domain Adaptation for Unsupervised Alignment of Word Embeddings

abstract: We propose a novel manifold based geometric approach for learning unsupervised alignment of word embeddings between the source and the target languages. Our approach formulates the alignment learning problem as a domain adaptation problem over the manifold of doubly stochastic matrices. This viewpoint arises from the aim to align the second order information of the two language spaces. The rich geometry of the doubly stochastic manifold allows to employ efficient Riemannian conjugate gradient algorithm for the proposed formulation. Empirically, the proposed approach outperforms state-of-the-art optimal transport based approach on the bilingual lexicon induction task across several language pairs. The performance improvement is more significant for distant language pairs.


title: The 2-rank of the class group of some real cyclic quartic number fields

abstract: In this paper, we investigate the 2-rank of the class group of some real cyclic quartic number fields. Precisely, we consider the case where the quadratic subfield is Q(\sqrt{l}) with l congruent to 5 modulo 8 is a prime.


title: Quality Prediction on Deep Generative Images

abstract: In recent years, deep neural networks have been utilized in a wide variety of applications including image generation. In particular, generative adversarial networks (GANs) are able to produce highly realistic pictures as part of tasks such as image compression. As with standard compression, it is desirable to be able to automatically assess the perceptual quality of generative images to monitor and control the encode process. However, existing image quality algorithms are ineffective on GAN generated content, especially on textured regions and at high compressions. Here we propose a new naturalness-based image quality predictor for generative images. Our new GAN picture quality predictor is built using a multi-stage parallel boosting system based on structural similarity features and measurements of statistical similarity. To enable model development and testing, we also constructed a subjective GAN image quality database containing (distorted) GAN images and collected human opinions of them. Our experimental results indicate that our proposed GAN IQA model delivers superior quality predictions on the generative image datasets, as well as on traditional image quality datasets.


title: Res-CR-Net, a residual network with a novel architecture optimized for the semantic segmentation of microscopy images

abstract: Deep Neural Networks (DNN) have been widely used to carry out segmentation tasks in both electron and light microscopy. Most DNNs developed for this purpose are based on some variation of the encoder-decoder type U-Net architecture, in combination with residual blocks to increase ease of training and resilience to gradient degradation. Here we introduce Res-CR-Net, a type of DNN that features residual blocks with either a bundle of separable atrous convolutions with different dilation rates or a convolutional LSTM. The number of filters used in each residual block and the number of blocks are the only hyperparameters that need to be modified in order to optimize the network training for a variety of different microscopy images.


title: Real-time sparse-sampled Ptychographic imaging through deep neural networks

abstract: Ptychography has rapidly grown in the fields of X-ray and electron imaging for its unprecedented ability to achieve nano or atomic scale resolution while simultaneously retrieving chemical or magnetic information from a sample. A ptychographic reconstruction is achieved by means of solving a complex inverse problem that imposes constraints both on the acquisition and on the analysis of the data, which typically precludes real-time imaging due to computational cost involved in solving this inverse problem. In this work we propose PtychoNN, a novel approach to solve the ptychography reconstruction problem based on deep convolutional neural networks. We demonstrate how the proposed method can be used to predict real-space structure and phase at each scan point solely from the corresponding far-field diffraction data. The presented results demonstrate how PtychoNN can effectively be used on experimental data, being able to generate high quality reconstructions of a sample up to hundreds of times faster than state-of-the-art ptychography reconstruction solutions once trained. By surpassing the typical constraints of iterative model-based methods, we can significantly relax the data acquisition sampling conditions and produce equally satisfactory reconstructions. Besides drastically accelerating acquisition and analysis, this capability can enable new imaging scenarios that were not possible before, in cases of dose sensitive, dynamic and extremely voluminous samples.

categories: eess.IV cond-mat.mes-hall cs.LG physics.optics


title: Acoustical classification of different speech acts using nonlinear methods

abstract: A recitation is a way of combining the words together so that they have a sense of rhythm and thus an emotional content is imbibed within. In this study we envisaged to answer these questions in a scientific manner taking into consideration 5 (five) well known Bengali recitations of different poets conveying a variety of moods ranging from joy to sorrow. The clips were recited as well as read (in the form of flat speech without any rhythm) by the same person to avoid any perceptual difference arising out of timbre variation. Next, the emotional content from the 5 recitations were standardized with the help of listening test conducted on a pool of 50 participants. The recitations as well as the speech were analyzed with the help of a latest non linear technique called Detrended Fluctuation Analysis (DFA) that gives a scaling exponent {\alpha}, which is essentially the measure of long range correlations present in the signal. Similar pieces (the parts which have the exact lyrical content in speech as well as in the recital) were extracted from the complete signal and analyzed with the help of DFA technique. Our analysis shows that the scaling exponent for all parts of recitation were much higher in general as compared to their counterparts in speech. We have also established a critical value from our analysis, above which a mere speech may become a recitation. The case may be similar to the conventional phase transition, wherein the measurement of external condition at which the transformation occurs (generally temperature) is called phase transition. Further, we have also categorized the 5 recitations on the basis of their emotional content with the help of the same DFA technique. Analysis with a greater variety of recitations is being carried out to yield more interesting results.


title: Understanding the Difficulty of Training Transformers

abstract: Transformers have proved effective in many NLP tasks. However, their training requires non-trivial efforts regarding designing cutting-edge optimizers and learning rate schedulers carefully (e.g., conventional SGD fails to train Transformers effectively). Our objective here is to understand $\textit{what complicates Transformer training}$ from both empirical and theoretical perspectives. Our analysis reveals that unbalanced gradients are not the root cause of the instability of training. Instead, we identify an amplification effect that influences training substantially -- for each layer in a multi-layer Transformer model, heavy dependency on its residual branch makes training unstable, since it amplifies small parameter perturbations (e.g., parameter updates) and results in significant disturbances in the model output. Yet we observe that a light dependency limits the model potential and leads to inferior trained models. Inspired by our analysis, we propose Admin ($\textbf{Ad}$aptive $\textbf{m}$odel $\textbf{in}$itialization) to stabilize stabilize the early stage's training and unleash its full potential in the late stage. Extensive experiments show that Admin is more stable, converges faster, and leads to better performance. Implementations are released at: https://github.com/LiyuanLucasLiu/Transforemr-Clinic.


title: How to Teach DNNs to Pay Attention to the Visual Modality in Speech Recognition

abstract: Audio-Visual Speech Recognition (AVSR) seeks to model, and thereby exploit, the dynamic relationship between a human voice and the corresponding mouth movements. A recently proposed multimodal fusion strategy, AV Align, based on state-of-the-art sequence to sequence neural networks, attempts to model this relationship by explicitly aligning the acoustic and visual representations of speech. This study investigates the inner workings of AV Align and visualises the audio-visual alignment patterns. Our experiments are performed on two of the largest publicly available AVSR datasets, TCD-TIMIT and LRS2. We find that AV Align learns to align acoustic and visual representations of speech at the frame level on TCD-TIMIT in a generally monotonic pattern. We also determine the cause of initially seeing no improvement over audio-only speech recognition on the more challenging LRS2. We propose a regularisation method which involves predicting lip-related Action Units from visual representations. Our regularisation method leads to better exploitation of the visual modality, with performance improvements between 7% and 30% depending on the noise level. Furthermore, we show that the alternative Watch, Listen, Attend, and Spell network is affected by the same problem as AV Align, and that our proposed approach can effectively help it learn visual representations. Our findings validate the suitability of the regularisation method to AVSR and encourage researchers to rethink the multimodal convergence problem when having one dominant modality.


title: Sufficient and necessary conditions for hereditary of infinite category algebras

abstract: We describe necessary and sufficient conditions for the hereditarity of the category algebra of an infinite EI category satisfying certain combinatorial assumptions. More generally, we discuss conditions such that the left global dimension of a category algebra equals the maximal left global dimension of the endomorphism algebras of its objects, and classify its projective modules in this case. As applications, we completely classify transporter categories, orbit categories, and Quillen categories with left hereditary category algebras over a field.


title: Comparison between the non-crossing and the non-crossing on lines properties

abstract: In the recent paper [2], it was proved that the closure of the planar diffeomorphisms in the Sobolev norm consists of the functions which are non-crossing (NC), i.e., the functions which can be uniformly approximated by continuous one-to-one functions on the grids. A deep simplification of this property is to consider curves instead of grids, so considering functions which are non-crossing on lines (NCL). Since the NCL property is way easier to check, it would be extremely positive if they actually coincide, while it is only obvious that NC implies NCL. We show that in general NCL does not imply NC, but the implication becomes true with the additional assumption that $\det(Du)>0$ a.e., which is a very common assumption in nonlinear elasticity.


title: Hadron and lepton tensors in semileptonic decays including new physics

abstract: We extend our general framework for semileptonic decay, originally introduced in N. Penalva et al. [Phys. Rev. D100, 113007 (2019)], with the addition of new physics (NP) tensor terms. In this way, all the NP effective Hamiltonians that are considered in lepton flavour universality violation (LFUV) studies have now been included. Besides, we now also give general expressions that allow for complex Wilson coefficients. The scheme developed is totally general and it can be applied to any charged current semileptonic decay, involving any quark flavors or initial and final hadron states. We show that all the hadronic input, including NP effects, can be parametrized in terms of 16 Lorentz scalar structure functions. In the second part of this work, we use this formalism to obtain the complete NP effects in the $\Lambda_b\to \Lambda_c \tau\bar\nu_\tau$ semileptonic decay, where LFUV, if finally confirmed, is also expected to be seen. We stress the relevance of the center of mass (CM) $d^2\Gamma/(d\omega d\cos\theta_\ell)$ and laboratory (LAB) $d^2\Gamma/(d\omega dE_\ell)$ differential decay widths. While models with very different strengths in the NP terms give the same differential $d\Gamma/d\omega$ and total decay widths for this decay, they predict very different numerical results for some of the $\cos\theta_\ell$ and $E_\ell$ coefficient-functions that determine the above two distributions. Thus, the combined analysis of the CM $d^2\Gamma/(d\omega d\cos\theta_\ell)$ and LAB $d^2\Gamma/ (d\omega dE_\ell)$ differential decay widths will help clarifying what kind of NP is a better candidate in order to explain LFUV.


title: Event reconstruction for KM3NeT/ORCA using convolutional neural networks

abstract: The KM3NeT research infrastructure is currently under construction at two locations in the Mediterranean Sea. The KM3NeT/ORCA water-Cherenkov neutrino detector off the French coast will instrument several megatons of seawater with photosensors. Its main objective is the determination of the neutrino mass ordering. This work aims at demonstrating the general applicability of deep convolutional neural networks to neutrino telescopes, using simulated datasets for the KM3NeT/ORCA detector as an example. To this end, the networks are employed to achieve reconstruction and classification tasks that constitute an alternative to the analysis pipeline presented for KM3NeT/ORCA in the KM3NeT Letter of Intent. They are used to infer event reconstruction estimates for the energy, the direction, and the interaction point of incident neutrinos. The spatial distribution of Cherenkov light generated by charged particles induced in neutrino interactions is classified as shower- or track-like, and the main background processes associated with the detection of atmospheric neutrinos are recognized. Performance comparisons to machine-learning classification and maximum-likelihood reconstruction algorithms previously developed for KM3NeT/ORCA are provided. It is shown that this application of deep convolutional neural networks to simulated datasets for a large-volume neutrino telescope yields competitive reconstruction results and performance improvements with respect to classical approaches.


title: A Survey of Approximate Quantile Computation on Large-scale Data (Technical Report)

abstract: As data volume grows extensively, data profiling helps to extract metadata of large-scale data. However, one kind of metadata, order statistics, is difficult to be computed because they are not mergeable or incremental. Thus, the limitation of time and memory space does not support their computation on large-scale data. In this paper, we focus on an order statistic, quantiles, and present a comprehensive analysis of studies on approximate quantile computation. Both deterministic algorithms and randomized algorithms that compute approximate quantiles over streaming models or distributed models are covered. Then, multiple techniques for improving the efficiency and performance of approximate quantile algorithms in various scenarios, such as skewed data and high-speed data streams, are presented. Finally, we conclude with coverage of existing packages in different languages and with a brief discussion of the future direction in this area.


title: On the usage of randomized p-values in the Schweder-Spjotvoll estimator

abstract: We are concerned with multiple test problems with composite null hypotheses and the estimation of the proportion $\pi_{0}$ of true null hypotheses. The Schweder-Spj\o tvoll estimator $\hat{\pi}_0$ utilizes marginal $p$-values and only works properly if the $p$-values that correspond to the true null hypotheses are uniformly distributed on $[0,1]$ ($\mathrm{Uni}[0,1]$-distributed). In the case of composite null hypotheses, marginal $p$-values are usually computed under least favorable parameter configurations (LFCs). Thus, they are stochastically larger than $\mathrm{Uni}[0,1]$ under non-LFCs in the null hypotheses. When using these LFC-based $p$-values, $\hat{\pi}_0$ tends to overestimate $\pi_{0}$. We introduce a new way of randomizing $p$-values that depends on a tuning parameter $c\in[0,1]$, such that $c=0$ and $c=1$ lead to $\mathrm{Uni}[0,1]$-distributed $p$-values, which are independent of the data, and to the original LFC-based $p$-values, respectively. For a certain value $c=c^{\star}$ the bias of $\hat{\pi}_0$ is minimized when using our randomized $p$-values. This often also entails a smaller mean squared error of the estimator as compared to the usage of the LFC-based $p$-values. We analyze these points theoretically, and we demonstrate them numerically in computer simulations under various standard statistical models.


title: Duplication Detection in Knowledge Graphs: Literature and Tools

abstract: In recent years, an increasing amount of knowledge graphs (KGs) have been created as a means to store cross-domain knowledge and billion of facts, which are the basis of costumers' applications like search engines. However, KGs inevitably have inconsistencies such as duplicates that might generate conflicting property values. Duplication detection (DD) aims to identify duplicated entities and resolve their conflicting property values effectively and efficiently. In this paper, we perform a literature review on DD methods and tools, and an evaluation of them. Our main contributions are a performance evaluation of DD tools in KGs, improvement suggestions, and a DD workflow to support future development of DD tools, which are based on desirable features detected through this study.


title: Initial forms and a notion of basis for tropical differential equations

abstract: We show that solution sets of systems of tropical differential equations can be characterised in terms of monomial-freeness of an initial ideal. We discuss a candidate definition of tropical differential basis and give a nonexistence result for such bases in an example.


title: Majorization-Minimization-Based Levenberg--Marquardt Method for Constrained Nonlinear Least Squares

abstract: A new Levenberg--Marquardt (LM) method for solving nonlinear least squares problems with convex constraints is described. Various versions of the LM method have been proposed, their main differences being in the choice of a damping parameter. In this paper, we propose a new rule for updating the parameter so as to achieve both global and local convergence even under the presence of a convex constraint set. The key to our results is a new perspective of the LM method from majorization-minimization methods. Specifically, we show that if the damping parameter is set in a specific way, the objective function of the standard subproblem in LM methods becomes an upper bound on the original objective function under certain standard assumptions. Our method solves a sequence of the subproblems approximately using an (accelerated) projected gradient method. It finds an $\epsilon$-stationary point after $O(\epsilon^{-2})$ computation and achieves local quadratic convergence for zero-residual problems under a local error bound condition. Numerical results on compressed sensing and matrix factorization show that our method converges faster in many cases than existing methods.


title: Forecasting Multi-Dimensional Processes over Graphs

abstract: The forecasting of multi-variate time processes through graph-based techniques has recently been addressed under the graph signal processing framework. However, problems in the representation and the processing arise when each time series carries a vector of quantities rather than a scalar one. To tackle this issue, we devise a new framework and propose new methodologies based on the graph vector autoregressive model. More explicitly, we leverage product graphs to model the high-dimensional graph data and develop multi-dimensional graph-based vector autoregressive models to forecast future trends with a number of parameters that is independent of the number of time series and a linear computational complexity. Numerical results demonstrating the prediction of moving point clouds corroborate our findings.


title: The Bogomolov-Beauville-Yau Decomposition for Klt Projective Varieties with Trivial First Chern Class -Without Tears-

abstract: We give a simplified proof (in characteristic zero) of the decomposition theorem for complex projective varieties with klt singularities and numerically trivial canonical bundle. The proof rests in an essential way on most of the partial results of the previous proof obtained by many authors, but avoids those in positive characteristic by S. Druel. The single to some extent new contribution is an algebraicity and bimeromorphic splitting result for generically locally trivial fibrations with fibres without holomorphic vector fields. We give first the proof in the easier smooth case, following the same steps as in the general case, treated next.


title: Identification of heavy, energetic, hadronically decaying particles using machine-learning techniques

abstract: Machine-learning (ML) techniques are explored to identify and classify hadronic decays of highly Lorentz-boosted W/Z/Higgs bosons and top quarks. Techniques without ML have also been evaluated and are included for comparison. The identification performances of a variety of algorithms are characterized in simulated events and directly compared with data. The algorithms are validated using proton-proton collision data at $\sqrt{s} =$ 13 TeV, corresponding to an integrated luminosity of 35.9 fb$^{-1}$. Systematic uncertainties are assessed by comparing the results obtained using simulation and collision data. The new techniques studied in this paper provide significant performance improvements over non-ML techniques, reducing the background rate by up to an order of magnitude at the same signal efficiency.


title: Leveraging Mobility Flows from Location Technology Platforms to Test Crime Pattern Theory in Large Cities

abstract: Crime has been previously explained by social characteristics of the residential population and, as stipulated by crime pattern theory, might also be linked to human movements of non-residential visitors. Yet a full empirical validation of the latter is lacking. The prime reason is that prior studies are limited to aggregated statistics of human visitors rather than mobility flows and, because of that, neglect the temporal dynamics of individual human movements. As a remedy, we provide the first work which studies the ability of granular human mobility in describing and predicting crime concentrations at an hourly scale. For this purpose, we propose the use of data from location technology platforms. This type of data allows us to trace individual transitions and, therefore, we succeed in distinguishing different mobility flows that (i) are incoming or outgoing from a neighborhood, (ii) remain within it, or (iii) refer to transitions where people only pass through the neighborhood. Our evaluation infers mobility flows by leveraging an anonymized dataset from Foursquare that includes almost 14.8 million consecutive check-ins in three major U.S. cities. According to our empirical results, mobility flows are significantly and positively linked to crime. These findings advance our theoretical understanding, as they provide confirmatory evidence for crime pattern theory. Furthermore, our novel use of digital location services data proves to be an effective tool for crime forecasting. It also offers unprecedented granularity when studying the connection between human mobility and crime.


title: Machine learning reveals complex behaviours in optically trapped particles

abstract: Since their invention in the 1980s [1], optical tweezers have found a wide range of applications, from biophotonics and mechanobiology to microscopy and optomechanics [2, 3, 4, 5]. Simulations of the motion of microscopic particles held by optical tweezers are often required to explore complex phenomena and to interpret experimental data [6, 7, 8, 9]. For the sake of computational efficiency, these simulations usually model the optical tweezers as an harmonic potential [10, 11]. However, more physically-accurate optical-scattering models [12, 13, 14, 15] are required to accurately model more onerous systems; this is especially true for optical traps generated with complex fields [16, 17, 18, 19]. Although accurate, these models tend to be prohibitively slow for problems with more than one or two degrees of freedom (DoF) [20], which has limited their broad adoption. Here, we demonstrate that machine learning permits one to combine the speed of the harmonic model with the accuracy of optical-scattering models. Specifically, we show that a neural network can be trained to rapidly and accurately predict the optical forces acting on a microscopic particle. We demonstrate the utility of this approach on two phenomena that are prohibitively slow to accurately simulate otherwise: the escape dynamics of swelling microparticles in an optical trap, and the rotation rates of particles in a superposition of beams with opposite orbital angular momenta. Thanks to its high speed and accuracy, this method can greatly enhance the range of phenomena that can be efficiently simulated and studied.


title: Two-Pion production in the second resonance region in $\pi^- p$ collisions with HADES

abstract: Pion induced reactions provide unique opportunities for an unambiguous description of baryonic resonances and their coupling channels by means of a partial wave analysis. Using the secondary pion beam at SIS18, the two pion production in the second resonance region has been investigated to unravel the role of the $N(1520) \frac{3}{2}^-$ resonance in the intermediate $\rho$ production. Results on exclusive channels with one pion ($\pi^{-}p$) and two pions ($\pi^{+}\pi^{-}n$, $\pi^{0}\pi^{-}p$) in the final state measured in the $\pi^{-}-p$ reaction at four different pion beam momenta (0.650, 0.685, 0.733, and 0.786 GeV/c) are presented. The excitation function of the different partial waves and $\Delta\pi$, $N\sigma$ and $N\rho$ isobar configurations is obtained, using the Bonn-Gatchina partial wave analysis. The $N(1520) \frac{3}{2}^-$ resonance is found to dominate the $N\rho$ final state with the branching ratio $BR=12.2 \pm 1.9\%$.


title: The 2-rank of the class group of some real cyclic quartic number fields II

abstract: In this paper, we determine the 2-rank of the class group of certain classes of real cyclic quartic number fields. Precisely, we consider the case in which the quadratic subfield is Q(\sqrt{l}) with l=2 or a prime congruent to 1 mod 8.


title: To switch or not to switch -- a machine learning approach for ferroelectricity

abstract: With the advent of increasingly elaborate experimental techniques in physics, chemistry and materials sciences, measured data are becoming bigger and more complex. The observables are typically a function of several stimuli resulting in multidimensional data sets spanning a range of experimental parameters. As an example, a common approach to study ferroelectric switching is to observe effects of applied electric field, but switching can also be enacted by pressure and is influenced by strain fields, material composition, temperature, time, etc. Moreover, the parameters are usually interdependent, so that their decoupling toward univariate measurements or analysis may not be straightforward. On the other hand, both explicit and hidden parameters provide an opportunity to gain deeper insight into the measured properties, provided there exists a well-defined path to capture and analyze such data. Here, we introduce a new, two-dimensional approach to represent hysteretic response of a material system to applied electric field. Utilizing ferroelectric polarization as a model hysteretic property, we demonstrate how explicit consideration of electromechanical response to two rather than one control voltages enables significantly more transparent and robust interpretation of observed hysteresis, such as differentiating between charge trapping and ferroelectricity. Furthermore, we demonstrate how the new data representation readily fits into a variety of machinelearning methodologies, from unsupervised classification of the origins of hysteretic response via linear clustering algorithms to neural-network-based inference of the sample temperature based on the specific morphology of hysteresis.


title: Development of gamma insensitive silicon carbide diagnostics to qualify intense thermal and epithermal neutron fields

abstract: The e_LiBANS project aims at creating accelerator based compact neutron facilities for diverse interdisciplinary applications. After the successful setting up and characterization of a thermal neutron source based on a medical electron LINAC, a similar assembly for epithermal neutrons has been developed. The project is based on an Elekta 18 MV LINAC coupled with a photoconverter-moderator system which deploys the ({\gamma},n) photonuclear reaction to convert a bremsstrahlung photon beam into a neutron field. This communication describes the development of novel diagnostics to qualify the thermal and epithermal neutron fields that have been produced. In particular, a proof of concept for the use of silicon carbide photodiodes as a thermal neutron rate detector is presented.


title: Beat Detection and Automatic Annotation of the Music of Bharatanatyam Dance using Speech Recognition Techniques

abstract: Bharatanatyam, an Indian Classical Dance form, represents the rich cultural heritage of India. Analysis and recognition of such dance forms are critical for the preservation of cultural heritage. Like in most dance forms, a Bharatanatyam dancer performs in synchronization with structured rhythmic music, called Sollukattu, which comprises instrumental beats and vocalized utterances (bols) to create a rhythmic music structure. Computer analysis of Bharatanatyam, therefore, requires a structural analysis of Sollukattus. In this paper, we use speech processing techniques to recognize bols. Exploiting the predefined structures of Sollukattus and the detected bols, we recognize the Sollukattu. We estimate the tempo period by two methods. Finally, we generate a complete annotation of the audio signal by beat marking. For this, we also use the information of beats detected from the onset envelope of a Sollukattu signal. For training and test, we create a data set for Sollukattus and annotate them. We achieve 85% accuracy in bol recognition, 95% in Sollukattu recognition, 96% in tempo period estimation, and over 90% in beat marking. This is the maiden attempt to fully structurally analyze the music of an Indian Classical Dance form and the use of speech processing techniques for beat marking.


title: Weakly Supervised Geodesic Segmentation of Egyptian Mummy CT Scans

abstract: In this paper, we tackle the task of automatically analyzing 3D volumetric scans obtained from computed tomography (CT) devices. In particular, we address a particular task for which data is very limited: the segmentation of ancient Egyptian mummies CT scans. We aim at digitally unwrapping the mummy and identify different segments such as body, bandages and jewelry. The problem is complex because of the lack of annotated data for the different semantic regions to segment, thus discouraging the use of strongly supervised approaches. We, therefore, propose a weakly supervised and efficient interactive segmentation method to solve this challenging problem. After segmenting the wrapped mummy from its exterior region using histogram analysis and template matching, we first design a voxel distance measure to find an approximate solution for the body and bandage segments. Here, we use geodesic distances since voxel features as well as spatial relationship among voxels is incorporated in this measure. Next, we refine the solution using a GrabCut based segmentation together with a tracking method on the slices of the scan that assigns labels to different regions in the volume, using limited supervision in the form of scribbles drawn by the user. The efficiency of the proposed method is demonstrated using visualizations and validated through quantitative measures and qualitative unwrapping of the mummy.


title: Representations of fusion categories and their commutants

abstract: A bicommutant category is a higher categorical analog of a von Neumann algebra. We study the bicommutant categories which arise as the commutant $\mathcal{C}'$ of a fully faithful representation $\mathcal{C}\to\operatorname{Bim}(R)$ of a unitary fusion category $\mathcal{C}$. Using results of Izumi, Popa, and Tomatsu about existence and uniqueness of representations of unitary (multi)fusion categories, we prove that if $\mathcal{C}$ and $\mathcal{D}$ are Morita equivalent unitary fusion categories, then their commutant categories $\mathcal{C}'$ and $\mathcal{D}'$ are equivalent as bicommutant categories. In particular, they are equivalent as tensor categories: \[ \Big(\,\,\mathcal{C}\,\,\simeq_{\text{Morita}}\,\,\mathcal{D}\,\,\Big) \qquad\Longrightarrow\qquad \Big(\,\,\mathcal{C}'\,\,\simeq_{\text{tensor}}\,\,\mathcal{D}'\,\,\Big). \] This categorifies the well-known result according to which the commutants (in some representations) of Morita equivalent finite dimensional $\rm C^*$-algebras are isomorphic von Neumann algebras, provided the representations are `big enough'. We also introduce a notion of positivity for bi-involutive tensor categories. For dagger categories, positivity is a property (the property of being a $\rm C^*$-category). But for bi-involutive tensor categories, positivity is extra structure. We show that unitary fusion categories and $\operatorname{Bim}(R)$ admit distinguished positive structures, and that fully faithful representations $\mathcal{C}\to\operatorname{Bim}(R)$ automatically respect these positive structures.


title: Board Games for Quantum Computers

abstract: Scalable board games, including Five in a Row (or gomoku) and weiqi (or go), are generalized so that they can be played on or by quantum computers. We adopt three principles for the generalization: the first two are to ensure that the games are compatible with quantum computer and the third is to ensure that the standard classical games are the special cases. We demonstrate how to construct basic quantum moves and use them to set up quantum games. There are three different schemes to play the quantized games: one quantum computer with another quantum computer (QwQ), two classical computer playing with each other on one quantum computer (CQC), and one classical computer with another classical computer(CwC). We illustrate these results with the games of Five in a Row and weiqi.


title: 3-D Context Entropy Model for Improved Practical Image Compression

abstract: In this paper, we present our image compression framework designed for CLIC 2020 competition. Our method is based on Variational AutoEncoder (VAE) architecture which is strengthened with residual structures. In short, we make three noteworthy improvements here. First, we propose a 3-D context entropy model which can take advantage of known latent representation in current spatial locations for better entropy estimation. Second, a light-weighted residual structure is adopted for feature learning during entropy estimation. Finally, an effective training strategy is introduced for practical adaptation with different resolutions. Experiment results indicate our image compression method achieves 0.9775 MS-SSIM on CLIC validation set and 0.9809 MS-SSIM on test set.


title: Observation of symmetry protected zero modes in topolectrical circuits

abstract: Higher-order topological insulators are a new class of topological phases of matter, originally conceived for electrons in solids. It has been suggested that $\mathbb{Z}_N$ Berry phase (Berry phase quantized into $2\pi/N$) is a useful tool to characterize the symmetry protected topological states, while the experimental evidence is still elusive. Recently, topolectrical circuits have emerged as a simple yet very powerful platform for studying topological physics that are challenging to realize in condensed matter systems. Here, we present the first experimental observation of second-order corner states characterized by $\mathbb{Z}_3$ Berry phase in topolectrical circuits. We demonstrate theoretically and experimentally that the localized second-order topological states are protected by a generalized chiral symmetry of tripartite lattices, and they are pinned to "zero energy". By introducing extra capacitors within sublattices in the circuit, we are able to examine the robustness of the zero modes against both chiral-symmetry conserving and breaking disturbances. Our work paves the way for testing exotic topological band theory by electrical-circuit experiments.


title: Quasiconformal Gauss maps and the Bernstein problem for Weingarten multigraphs

abstract: We prove that any complete, uniformly elliptic Weingarten surface in Euclidean $3$-space whose Gauss map image omits an open hemisphere is a cylinder or a plane. This generalizes a classical theorem by Hoffman, Osserman and Schoen for constant mean curvature surfaces. In particular, this proves that planes are the only complete, uniformly elliptic Weingarten multigraphs. We also show that this result holds for a large class of non-uniformly elliptic Weingarten equations. In particular, this solves in the affirmative the Bernstein problem for entire graphs for that class of elliptic equations. To obtain these results, we prove that planes are the only complete multigraphs with quasiconformal Gauss map and bounded second fundamental form.


title: Longitudinal spin fluctuations in bcc and liquid Fe at high temperature and pressure calculated with a supercell approach

abstract: Investigation of magnetic materials at realistic conditions with first-principles methods is a challenging task due to the interplay of vibrational and magnetic degrees of freedom. The most difficult contribution to include in simulations is represented by the longitudinal magnetic degrees of freedom (LSF) due to their inherent many-body nature; nonetheless, schemes that enable to take into account this effect on a semiclassical level have been proposed and employed in the investigation of magnetic systems. However, assessment of the effect of vibrations on LSF is lacking in the literature. For this reason, in this work we develop a supercell approach within the framework of constrained density functional theory to calculate self-consistently the size of local-environment-dependent magnetic moments in the paramagnetic, high-temperature state in presence of lattice vibrations and for liquid Fe in different conditions. First, we consider the case of bcc Fe at the Curie temperature and ambient pressure. Then, we perform a similar analysis on bcc Fe at Earth's inner core conditions, and we find that LSF stabilize non-zero moments which affect atomic forces and electronic density of states of the system. Finally, we employ the present scheme on liquid Fe at the melting point at ambient pressure, and at Earth's outer core conditions ($p \approx 200$ GPa, $T \approx 6000$ K). In both cases, we obtain local magnetic moments of sizes comparable to the solid-state counterparts.


title: Learning Strategies for Radar Clutter Classification

abstract: In this paper, we address the problem of classifying clutter returns in order to partition them into statistically homogeneous subsets. The classification procedure relies on a model for the observables including latent variables that is solved by the expectation-maximization algorithm. The derivations are carried out by accounting for three different cases for the structure of the clutter covariance matrix. A preliminary performance analysis highlights that the proposed technique is a viable means to cluster clutter returns over the range.


title: The impact of CMEs on the critical frequency of F2-layer ionosphere (foF2)

abstract: The ionospheric critical frequency (foF2) from ionosonde measurements at geographic high, middle, and low latitudes is analyzed with the occurrence of coronal mass ejections (CMEs) in long term variability of the solar cycles. We observed trends of monthly maximum foF2 values and monthly averaged values of CME parameters such as speed, angular width, mass, and kinetic energy with respect to time. The impact of CMEs on foF2 is very high at high latitudes and low at low latitudes. The time series for monthly maximum foF2 and monthly-averaged CME speed is moderately correlated at high and middle latitudes.


title: Bypassing or flying above the obstacles? A novel multi-objective UAV path planning problem

abstract: This study proposes a novel multi-objective integer programming model for a collision-free discrete drone path planning problem. Considering the possibility of bypassing obstacles or flying above them, this study aims to minimize the path length, energy consumption, and maximum path risk simultaneously. The static environment is represented as 3D grid cells. Due to the NP-hardness nature of the problem, several state-of-theart evolutionary multi-objective optimization (EMO) algorithms with customized crossover and mutation operators are applied to find a set of non-dominated solutions. The results show the effectiveness of applied algorithms in solving several generated test cases.


title: In situ evidence of firehose instability in multiple reconnection

abstract: Energy conversion via reconnecting current sheets is common in space and astrophysical plasmas. Frequently, current sheets disrupt at multiple reconnection sites, leading to the formation of plasmoid structures between sites, which might affect energy conversion. We present in situ evidence of the firehose instability in multiple reconnection in the Earth's magnetotail. The observed proton beams accelerated in the direction parallel to magnetic field and ion-scale fluctuations of whistler type imply the development of firehose instability between two active reconnection sites. The linear wave dispersion relation, estimated for the measured plasma parameters, indicates a positive growth rate of firehose-related electromagnetic fluctuations. Simulations of temporal evolution of the observed multiple reconnection by using a 2.5D implicit particle-in-cell code show that, as the plasmoid formed between two reconnection sites evolves, the plasma at its edge becomes anisotropic and overcomes the firehose marginal stability threshold, leading to the generation of magnetic field fluctuations. The combined results of observations and simulations suggest that the firehose instability, operating between reconnection sites, converts plasma kinetic energy into energy of magnetic field fluctuations, counteracting the conversion of magnetic energy into plasma energy occurring at reconnection sites. This suggests that magnetic energy conversion in multiple reconnection can be less efficient than in the case of the single-site reconnection.

categories: physics.space-ph astro-ph.EP astro-ph.SR physics.plasm-ph


title: Perfect depolarization in single scattering of light from uncorrelated surface and volume disorder

abstract: We demonstrate that single scattering of p-polarized waves from uncorrelated surface and volume disorder can lead to perfect depolarization. The degree of polarization is shown to vanish in specific scattering directions that can be characterized based on simple geometric arguments. Depolarization results from the different polarization response of each source of disorder, providing a clear physical interpretation of the depolarization mechanism.


title: Distributed Lower Bounds for Ruling Sets

abstract: Given a graph $G = (V,E)$, an $(\alpha, \beta)$-ruling set is a subset $S \subseteq V$ such that the distance between any two vertices in $S$ is at least $\alpha$, and the distance between any vertex in $V$ and the closest vertex in $S$ is at most $\beta$. We present lower bounds for distributedly computing ruling sets. More precisely, for the problem of computing a $(2, \beta)$-ruling set in the LOCAL model, we show the following, where $n$ denotes the number of vertices, $\Delta$ the maximum degree, and $c$ is some universal constant independent of $n$ and $\Delta$. $\bullet$ Any deterministic algorithm requires $\Omega\left(\min \left\{ \frac{\log \Delta}{\beta \log \log \Delta} , \log_\Delta n \right\} \right)$ rounds, for all $\beta \le c \cdot \min\left\{ \sqrt{\frac{\log \Delta}{\log \log \Delta}} , \log_\Delta n \right\}$. By optimizing $\Delta$, this implies a deterministic lower bound of $\Omega\left(\sqrt{\frac{\log n}{\beta \log \log n}}\right)$ for all $\beta \le c \sqrt[3]{\frac{\log n}{\log \log n}}$. $\bullet$ Any randomized algorithm requires $\Omega\left(\min \left\{ \frac{\log \Delta}{\beta \log \log \Delta} , \log_\Delta \log n \right\} \right)$ rounds, for all $\beta \le c \cdot \min\left\{ \sqrt{\frac{\log \Delta}{\log \log \Delta}} , \log_\Delta \log n \right\}$. By optimizing $\Delta$, this implies a randomized lower bound of $\Omega\left(\sqrt{\frac{\log \log n}{\beta \log \log \log n}}\right)$ for all $\beta \le c \sqrt[3]{\frac{\log \log n}{\log \log \log n}}$. For $\beta > 1$, this improves on the previously best lower bound of $\Omega(\log^* n)$ rounds that follows from the 30-year-old bounds of Linial [FOCS'87] and Naor [J.Disc.Math.'91]. For $\beta = 1$, i.e., for the problem of computing a maximal independent set, our results improve on the previously best lower bound of $\Omega(\log^* n)$ on trees, as our bounds already hold on trees.


title: Learned Video Compression with Feature-level Residuals

abstract: In this paper, we present an end-to-end video compression network for P-frame challenge on CLIC. We focus on deep neural network (DNN) based video compression, and improve the current frameworks from three aspects. First, we notice that pixel space residuals is sensitive to the prediction errors of optical flow based motion compensation. To suppress the relative influence, we propose to compress the residuals of image feature rather than the residuals of image pixels. Furthermore, we combine the advantages of both pixel-level and feature-level residual compression methods by model ensembling. Finally, we propose a step-by-step training strategy to improve the training efficiency of the whole framework. Experiment results indicate that our proposed method achieves 0.9968 MS-SSIM on CLIC validation set and 0.9967 MS-SSIM on test set.


title: Time Series Data Cleaning with Regular and Irregular Time Intervals

abstract: Errors are prevalent in time series data, especially in the industrial field. Data with errors could not be stored in the database, which results in the loss of data assets. Handling the dirty data in time series is non-trivial, when given irregular time intervals. At present, to deal with these time series containing errors, besides keeping original erroneous data, discarding erroneous data and manually checking erroneous data, we can also use the cleaning algorithm widely used in the database to automatically clean the time series data. This survey provides a classification of time series data cleaning techniques and comprehensively reviews the state-of-the-art methods of each type. In particular, we have a special focus on the irregular time intervals. Besides we summarize data cleaning tools, systems and evaluation criteria from research and industry. Finally, we highlight possible directions time series data cleaning.


title: Learning Based Hybrid Beamforming Design for Full-Duplex Millimeter Wave Systems

abstract: Millimeter Wave (mmWave) communications with full-duplex (FD) have the potential of increasing the spectral efficiency, relative to those with half-duplex. However, the residual self-interference (SI) from FD and high pathloss inherent to mmWave signals may degrade the system performance. Meanwhile, hybrid beamforming (HBF) is an efficient technology to enhance the channel gain and mitigate interference with reasonable complexity. However, conventional HBF approaches for FD mmWave systems are based on optimization processes, which are either too complex or strongly rely on the quality of channel state information (CSI). We propose two learning schemes to design HBF for FD mmWave systems, i.e., extreme learning machine based HBF (ELM-HBF) and convolutional neural networks based HBF (CNN-HBF). Specifically, we first propose an alternating direction method of multipliers (ADMM) based algorithm to achieve SI cancellation beamforming, and then use a majorization-minimization (MM) based algorithm for joint transmitting and receiving HBF optimization. To train the learning networks, we simulate noisy channels as input, and select the hybrid beamformers calculated by proposed algorithms as targets. Results show that both learning based schemes can provide more robust HBF performance and achieve at least 22.1% higher spectral efficiency compared to orthogonal matching pursuit (OMP) algorithms. Besides, the online prediction time of proposed learning based schemes is almost 20 times faster than the OMP scheme. Furthermore, the training time of ELM-HBF is about 600 times faster than that of CNN-HBF with 64 transmitting and receiving antennas.


title: Greenhouse Gas Emission Prediction on Road Network using Deep Sequence Learning

abstract: Mitigating the substantial undesirable impact of transportation systems on the environment is paramount. Thus, predicting Greenhouse Gas (GHG) emissions is one of the profound topics, especially with the emergence of intelligent transportation systems (ITS). We develop a deep learning framework to predict link-level GHG emission rate (ER) (in CO2eq gram/second) based on the most representative predictors, such as speed, density, and the GHG ER of previous time steps. In particular, various specifications of the long-short term memory (LSTM) networks with exogenous variables are examined and compared with clustering and the autoregressive integrated moving average (ARIMA) model with exogenous variables. The downtown Toronto road network is used as the case study and highly detailed data are synthesized using a calibrated traffic microsimulation and MOVES. It is found that LSTM specification with speed, density, GHG ER, and in-links speed from three previous minutes performs the best while adopting 2 hidden layers and when the hyper-parameters are systematically tuned. Adopting a 30 second updating interval improves slightly the correlation between true and predicted GHG ERs, but contributes negatively to the prediction accuracy as reflected on the increased root mean square error (RMSE) value. Efficiently predicting GHG emissions at a higher frequency with lower data requirements will pave the way to non-myopic eco-routing on large-scale road networks {to alleviate the adverse impact on the global warming


title: Deep Neural Network for Respiratory Sound Classification in Wearable Devices Enabled by Patient Specific Model Tuning

abstract: The primary objective of this paper is to build classification models and strategies to identify breathing sound anomalies (wheeze, crackle) for automated diagnosis of respiratory and pulmonary diseases. In this work we propose a deep CNN-RNN model that classifies respiratory sounds based on Mel-spectrograms. We also implement a patient specific model tuning strategy that first screens respiratory patients and then builds patient specific classification models using limited patient data for reliable anomaly detection. Moreover, we devise a local log quantization strategy for model weights to reduce the memory footprint for deployment in memory constrained systems such as wearable devices. The proposed hybrid CNN-RNN model achieves a score of 66.31% on four-class classification of breathing cycles for ICBHI'17 scientific challenge respiratory sound database. When the model is re-trained with patient specific data, it produces a score of 71.81% for leave-one-out validation. The proposed weight quantization technique achieves ~4X reduction in total memory cost without loss of performance. The main contribution of the paper is as follows: Firstly, the proposed model is able to achieve state of the art score on the ICBHI'17 dataset. Secondly, deep learning models are shown to successfully learn domain specific knowledge when pre-trained with breathing data and produce significantly superior performance compared to generalized models. Finally, local log quantization of trained weights is shown to be able to reduce the memory requirement significantly. This type of patient-specific re-training strategy can be very useful in developing reliable long-term automated patient monitoring systems particularly in wearable healthcare solutions.


title: A data driven analysis and forecast of an SEIARD epidemic model for COVID-19 in Mexico

abstract: We propose an SEIARD mathematical model to investigate the current outbreak of coronavirus disease (COVID-19) in Mexico. We conduct a detailed analysis of this model and demonstrate its application using publicly reported data. We calculate the basic reproduction number ($R_0$) via the next-generation matrix method, and we estimate the per day infection, death and recovery rates. We calibrate the parameters of the SEIARD model to the reported data by minimizing the sum of squared errors and attempt to forecast the evolution of the outbreak until June 2020. Our results estimate that the peak of the epidemic in Mexico will be around May 2, 2020. Our model incorporates the importance of considering the aysmptomatic infected individuals, because they represent the majority of the infected population (with symptoms or not) and they could play a huge role in spreading the virus without any knowledge.


title: Disentangled Adversarial Transfer Learning for Physiological Biosignals

abstract: Recent developments in wearable sensors demonstrate promising results for monitoring physiological status in effective and comfortable ways. One major challenge of physiological status assessment is the problem of transfer learning caused by the domain inconsistency of biosignals across users or different recording sessions from the same user. We propose an adversarial inference approach for transfer learning to extract disentangled nuisance-robust representations from physiological biosignal data in stress status level assessment. We exploit the trade-off between task-related features and person-discriminative information by using both an adversary network and a nuisance network to jointly manipulate and disentangle the learned latent representations by the encoder, which are then input to a discriminative classifier. Results on cross-subjects transfer evaluations demonstrate the benefits of the proposed adversarial framework, and thus show its capabilities to adapt to a broader range of subjects. Finally we highlight that our proposed adversarial transfer learning approach is also applicable to other deep feature learning frameworks.


title: Empirical Study of Market Impact Conditional on Order-Flow Imbalance

abstract: In this research, we have empirically investigated the key drivers affecting liquidity in equity markets. We illustrated how theoretical models, such as Kyle's model, of agents' interplay in the financial markets, are aligned with the phenomena observed in publicly available trades and quotes data. Specifically, we confirmed that for small signed order-flows, the price impact grows linearly with increase in the order-flow imbalance. We have, further, implemented a machine learning algorithm to forecast market impact given a signed order-flow. Our findings suggest that machine learning models can be used in estimation of financial variables; and predictive accuracy of such learning algorithms can surpass the performance of traditional statistical approaches. Understanding the determinants of price impact is crucial for several reasons. From a theoretical stance, modelling the impact provides a statistical measure of liquidity. Practitioners adopt impact models as a pre-trade tool to estimate expected transaction costs and optimize the execution of their strategies. This further serves as a post-trade valuation benchmark as suboptimal execution can significantly deteriorate a portfolio performance. More broadly, the price impact reflects the balance of liquidity across markets. This is of central importance to regulators as it provides an all-encompassing explanation of the correlation between market design and systemic risk, enabling regulators to design more stable and efficient markets.


title: Longest cycles in 3-connected hypergraphs and bipartite graphs

abstract: In the language of hypergraphs, our main result is a Dirac-type bound: we prove that every $3$-connected hypergraph $H$ with $ \delta(H)\geq \max\{|V(H)|, \frac{|E(H)|+10}{4}\}$ has a hamiltonian Berge cycle. This is sharp and refines a conjecture by Jackson from 1981 (in the language of bipartite graphs). Our proofs are in the language of bipartite graphs, since the incidence graph of each hypergraph is bipartite.


title: Boson and neutron stars with increased density

abstract: We discuss boson stars and neutron stars, respectively, in a scalar-tensor gravity model with an explicitly time-dependent real scalar field. While the boson stars in our model -- in contrast to the neutron stars -- do not possess a hard core, we find that the qualitative effects of the formation of scalar hair are similar in both cases : the presence of the gravity scalar allows both type of stars to exist for larger central density as well as larger mass at given radius than their General Relativity counterparts. In particular, we find new types of neutron stars with scalar hair which have radii very close to the corresponding Schwarzschild radius and hence are comparable in density to black holes. This new branch of solutions is stable with respect to the decay into individual baryons.


title: Quarkyonic Matter Equation of State in Beta-Equilibrium

abstract: Quark matter may appear due to a hadronic-quark transition in the core of a hybrid star. Quarkyonic matter is an approach in which both quarks and nucleons appear as quasi-particles in a crossover transition, and provides an explicit realization of early ideas concerning quark matter (e.g., the MIT bag model). This description has recently been employed by McLerran and Reddy to model chargeless (pure neutron) matter with an approach that has the virtue that the speed of sound rises quickly at a neutron-quark transition so as to satisfy observational constraints on the neutron star maximum mass ($\gtrsim2M_\odot$) and the radius of a $1.4M_\odot$ star ($R_{1.4}\lesssim 13.5$ km). Traditional models involving first-order transitions result in softer pressure-energy density relations that have difficulty satisfying these constraints except with very narrow choices of parameters. We propose a variation of quarkyonic matter involving protons and leptons whose energy can be explicitly minimized to achieve both chemical and beta equilibrium, which cannot be done in the chargeless formulation. Quarkyonic stellar models are able to satisfy observed mass and radius constraints with a wide range of model parameters, avoiding the obligatory fine-tuning of conventional hybrid star models, including requiring the transition density to be very close to the nuclear saturation density. Our formulation fits experimental and theoretical properties of the nuclear symmetry energy and pure neutron matter, and contains as few as three free parameters. This makes it an ideal tool for the study of high-density matter that is an efficient alternative to piecewise polytrope or spectral decomposition methods.


title: Dimension of Restricted Classes of Interval Orders

abstract: Rabinovitch showed in 1978 that the interval orders having a representation consisting of only closed unit intervals have order dimension at most 3. This article shows that the same dimension bound applies to two other classes of posets: those having a representation consisting of unit intervals (but with a mixture of open and closed intervals allowed) and those having a representation consisting of closed intervals with lengths in $\{0,1\}$.


title: Self-Duality in the Context of the Skyrme Model

abstract: We study a recently proposed modification of the Skyrme model that possesses an exact self-dual sector leading to an infinity of exact Skyrmion solutions with arbitrary topological (baryon) charge. The self-dual sector is made possible by the introduction, in addition to the usual three SU(2) Skyrme fields, of six scalar fields assembled in a symmetric and invertible three dimensional matrix h. The action presents quadratic and quartic terms in derivatives of the Skyrme fields, but instead of the group indices being contracted by the SU(2) Killing form, they are contracted with the h-matrix in the quadratic term, and by its inverse on the quartic term. Due to these extra fields the static version of the model, as well as its self-duality equations, are conformally invariant on the three dimensional space R^3. We show that the static and self-dual sectors of such a theory are equivalent, and so the only non-self-dual solution must be time dependent. We also show that for any configuration of the Skyrme SU(2) fields, the h-fields adjust themselves to satisfy the self-duality equations, and so the theory has plenty of non-trivial topological solutions. We present explicit exact solutions using a holomorphic rational ansatz, as well as a toroidal ansatz based on the conformal symmetry. We point to possible extensions of the model that break the conformal symmetry as well as the self-dual sector, and that can perhaps lead to interesting physical applications.


title: Directional Anisotropy of the Vibrational Modes in 2D Layered Perovskites

abstract: The vibrational modes in organic/inorganic layered perovskites are of fundamental importance for their optoelectronic properties. The hierarchical architecture of the Ruddlesden-Popper phase of these materials allows for distinct directionality of the vibrational modes withrespect to the main axes of the pseudocubic lattice in the octahedral plane. Here, we study the directionality of the fundamental phonon modes in single exfoliated Ruddlesden-Popper perovskite flakes with polarized Raman spectroscopy at ultralow-frequencies. A wealth of Raman bands is distinguished in the range from 15-150 cm-1 (2-15 meV), whose features depend on the organic cation species, on temperature, and on the direction of the linear polarization of the incident light. By controlling the angle of the linear polarization of the excitation laser with respect to the in-plane axes of the octahedral layer, we gain detailed information on the symmetry of the vibrational modes. The choice of two different organic moieties, phenethylammonium (PEA) and butylammonium (BA) allows to discern the influence of the linker molecules, evidencing strong anisotropy of the vibrations for the (PEA)2PbBr4 samples. Temperature dependent Raman measurements reveal that the broad phonon bands observed at room temperature consist of a series of sharp modes, and that such mode splitting strongly differs for the different organic moieties and vibrational bands.


title: Towards data-driven stroke rehabilitation via wearable sensors and deep learning

abstract: Recovery after stroke is often incomplete, but rehabilitation training may potentiate recovery by engaging endogenous neuroplasticity. In preclinical models of stroke, high doses of rehabilitation training are required to restore functional movement to the affected limbs of animals. In humans, however, the necessary dose of training to potentiate recovery is not known. This ignorance stems from the lack of objective, pragmatic approaches for measuring training doses in rehabilitation activities. Here, to develop a measurement approach, we took the critical first step of automatically identifying functional primitives, the basic building block of activities. Forty-eight individuals with chronic stroke performed a variety of rehabilitation activities while wearing inertial measurement units (IMUs) to capture upper body motion. Primitives were identified by human labelers, who labeled and segmented the associated IMU data. We performed automatic classification of these primitives using machine learning. We designed a convolutional neural network model that outperformed existing methods. The model includes an initial module to compute separate embeddings of different physical quantities in the sensor data. In addition, it replaces batch normalization (which performs normalization based on statistics computed from the training data) with instance normalization (which uses statistics computed from the test data). This increases robustness to possible distributional shifts when applying the method to new patients. With this approach, we attained an average classification accuracy of 70%. Thus, using a combination of IMU-based motion capture and deep learning, we were able to identify primitives automatically. This approach builds towards objectively-measured rehabilitation training, enabling the identification and counting of functional primitives that accrues to a training dose.


title: IDDA: a large-scale multi-domain dataset for autonomous driving

abstract: Semantic segmentation is key in autonomous driving. Using deep visual learning architectures is not trivial in this context, because of the challenges in creating suitable large scale annotated datasets. This issue has been traditionally circumvented through the use of synthetic datasets, that have become a popular resource in this field. They have been released with the need to develop semantic segmentation algorithms able to close the visual domain shift between the training and test data. Although exacerbated by the use of artificial data, the problem is extremely relevant in this field even when training on real data. Indeed, weather conditions, viewpoint changes and variations in the city appearances can vary considerably from car to car, and even at test time for a single, specific vehicle. How to deal with domain adaptation in semantic segmentation, and how to leverage effectively several different data distributions (source domains) are important research questions in this field. To support work in this direction, this paper contributes a new large scale, synthetic dataset for semantic segmentation with more than 100 different source visual domains. The dataset has been created to explicitly address the challenges of domain shift between training and test data in various weather and view point conditions, in seven different city types. Extensive benchmark experiments assess the dataset, showcasing open challenges for the current state of the art. The dataset will be available at: https://idda-dataset.github.io/home/ .


title: DSTC8-AVSD: Multimodal Semantic Transformer Network with Retrieval Style Word Generator

abstract: Audio Visual Scene-aware Dialog (AVSD) is the task of generating a response for a question with a given scene, video, audio, and the history of previous turns in the dialog. Existing systems for this task employ the transformers or recurrent neural network-based architecture with the encoder-decoder framework. Even though these techniques show superior performance for this task, they have significant limitations: the model easily overfits only to memorize the grammatical patterns; the model follows the prior distribution of the vocabularies in a dataset. To alleviate the problems, we propose a Multimodal Semantic Transformer Network. It employs a transformer-based architecture with an attention-based word embedding layer that generates words by querying word embeddings. With this design, our model keeps considering the meaning of the words at the generation stage. The empirical results demonstrate the superiority of our proposed model that outperforms most of the previous works for the AVSD task.


title: Thermodynamic signatures of an antiferromagnetic quantum critical point inside a superconducting dome

abstract: Recent experiments in unconventional superconductors, and in particular iron-based materials, have reported evidence of an antiferromagnetic quantum critical point (AFM-QCP) emerging inside the superconducting dome of the phase diagram. Fluctuations associated with such an AFM-QCP are expected to promote unusual temperature dependencies of thermodynamic quantities. Here, we compute the $T$ dependence of the specific heat $C(T)$ deep inside a fully gapped $s^{+-}$ superconducting state as the AFM-QCP is approached. We find that, at the AFM-QCP, the specific heat $C(T)$ vanishes quadratically with temperature, as opposed to the typical exponential suppression seen in fully-gapped BCS superconductors. This robust result is due to a non-analytic contribution to the free-energy arising from the general form of the bosonic (AFM) propagator in the SC state. Away from the AFM-QCP, as temperature is lowered, $C(T)$ shows a crossover from a $T^2$ behavior to an exponential behavior, with the crossover temperature scale set by the value of the superconducting gap and the distance to the QCP. We argue that these features in the specific heat can be used to unambiguously determine the existence of AFM-QCPs inside the superconducting domes of iron-based and other fully gapped unconventional superconductors.


title: Belief Propagation for Maximum Coverage on Weighted Bipartite Graph and Application to Text Summarization

abstract: We study text summarization from the viewpoint of maximum coverage problem. In graph theory, the task of text summarization is regarded as maximum coverage problem on bipartite graph with weighted nodes. In recent study, belief-propagation based algorithm for maximum coverage on unweighted graph was proposed using the idea of statistical mechanics. We generalize it to weighted graph for text summarization. Then we apply our algorithm to weighted biregular random graph for verification of maximum coverage performance. We also apply it to bipartite graph representing real document in open text dataset, and check the performance of text summarization. As a result, our algorithm exhibits better performance than greedy-type algorithm in some setting of text summarization.


title: IMRPhenomTP: A phenomenological time domain model for dominant quadrupole gravitational wave signal of coalescing binary black holes

abstract: In this work we present IMRPhenomTP, a time domain phenomenological model for the dominant $l=2$, $m=|2|$ modes of coalescing black hole binary systems and its extension to describe general precessing systems within the "twisting up" approximation. The underlying non-precessing model is calibrated to the new release of Numerical Relativity simulations of the SXS Collaboration and its accuracy is comparable to the state-of-the-art non-precessing dominant mode models as IMRPhenomX and SEOBNRv4. The precessing extension allows for flexibility choosing the Euler angles of the time-dependent rotation between the co-precessing and the inertial reference systems, including the single spin NNLO and the double spin MSA PN descriptions present in other models, numerical integration of the orbit averaged spin evolution equations, different choices for the evolution of the orbital angular momentum norm and a simple approximation to the ringdown behaviour.


title: Canonical Maps of general Hypersurfaces in Abelian Varieties

abstract: The main theorem of this paper is that, for a general pair $(A,X)$ of an (ample) Hypersurface $X$ in an Abelian Variety $A$, the canonical map $\Phi_X$ of $X$ is birational onto its image if the polarization given by $X$ is not principal (i.e., its Pfaffian $d$ is not equal to $1$). We also show that, setting $g = dim (A)$, and letting $d$ be the Pfaffian of the polarization given by $X$, then if $X$ is smooth and $$\Phi_X : X \rightarrow \mathbb{P}^{N:=g+d-2}$$ is an embedding, then necessarily we have the inequality $ d \geq g + 1$, equivalent to $N : = g+d-2 \geq 2 \ dim(X) + 1.$ We also formulate the following interesting conjecture, motivated by work of the second author: if $ d \geq g + 1,$ then, for a general pair $(A,X)$, $\Phi_X$ is an embedding.


title: Travelling wave solutions on a non-zero background for the generalized Korteweg-de Vries equation

abstract: For the generalized $p$-power Korteweg-de Vries equation, all non-periodic travelling wave solutions with non-zero boundary conditions are explicitly classified for all integer powers $p\geq 1$. These solutions are shown to consist of: bright solitary waves and static humps on a non-zero background for odd $p$; dark solitary waves on a non-zero background and kink waves for even $p$ in the defocusing case; pairs of bright/dark solitary waves on a non-zero background, and also bright and dark heavy-tail waves (with power decay) on a non-zero background, for even $p$ in the focusing case. An explicit physical parameterization is given for each type of solutionin terms of the wave speed $c$, background size $b$, and wave height/depth $h$. The allowed kinematic region in $(c,b)$ as well as in $(h,b)$ for existence of the solutions is derived, and other main kinematic features are discussed. Explicit formulas are presented in the integrable cases $p=1,2$, and in the higher power cases $p=3,4$.


title: Symmetries of the Schroedinger-Pauli equation for neutral particles

abstract: With using the algebraic approach Lie symmetries of Schr\"odinger equations with matrix potentials are classified. Thirty three inequivalent equations of such type together with the related symmetry groups are specified, the admissible equivalence relations are clearly indicated. In particular the Boyer results concerning kinematical invariance groups for arbitrary potentials (C. P. Boyer, Helv. Phys. Acta, {\bf 47}, 450--605 (1974)) are clarified and corrected.


title: On first order amenability

abstract: We introduce the notion of first order amenability, as a property of a first order theory $T$: every complete type over $\emptyset$, in possibly infinitely many variables, extends to an automorphism-invariant global Keisler measure in the same variables. Amenability of $T$ follows from amenability of the (topological) group $Aut(M)$ for all sufficiently large $\aleph_{0}$-homogeneous countable models $M$ of $T$ (assuming $T$ to be countable), but is radically less restrictive. First, we study basic properties of amenable theories, giving many equivalent conditions. Then, applying a version of the stabilizer theorem from [Amenability, connected components, and definable actions; E. Hrushovski, K. Krupi\'{n}ski, A. Pillay], we prove that if $T$ is amenable, then $T$ is G-compact, namely Lascar strong types and Kim-Pillay strong types over $\emptyset$ coincide. This extends and essentially generalizes a similar result proved via different methods for $\omega$-categorical theories in [Amenability, definable groups, and automorphism groups; K. Krupi\'{n}ski, A. Pillay] . In the special case when amenability is witnessed by $\emptyset$-definable global Keisler measures (which is for example the case for amenable $\omega$-categorical theories), we also give a different proof, based on stability in continuous logic. Parallel (but easier) results hold for the notion of extreme amenability.


title: Fast self-testing Quantum Random Number Generator based on homodyne detection

abstract: Self-testing and Semi-Device Independent protocols are becoming the preferred choice for quantum technologies, being able to certify their quantum nature with few assumptions and simple experimental implementations. In particular for Quantum Random Number Generators the possibility of monitoring in real time the entropy of the source only by measuring the input/output statistics is a characteristic that no other classical system could provide. The cost of this new possibility is not necessarily increased complexity and reduced performance. Indeed, here we show that with a simple optical setup consisting of commercially available components, a high bit generation rate can be achieved. We manage to certify 145.5~MHz of quantum random bit generation rate.


title: Fast tests for probing the causal structure of quantum processes

abstract: The identification of causal relations is a cornerstone of the scientific method. Traditional approaches to this task are based on classical statistics. However, such classical approaches do not apply in the quantum domain, where a broader spectrum of causal relations becomes accessible. New approaches to quantum causal inference have been developed in recent years, and promising new features have been discovered. In this paper, we review and partly expand the framework and results of Ref. [1], which demonstrated quantum speedups in the identification of various types of causal relations induced by reversible processes.


title: Bayesian semiparametric long memory models for discretized event data

abstract: We introduce a new class of semiparametric latent variable models for long memory discretized event data. The proposed methodology is motivated by a study of bird vocalizations in the Amazon rain forest; the timings of vocalizations exhibit self-similarity and long range dependence ruling out models based on Poisson processes. The proposed class of FRActional Probit (FRAP) models is based on thresholding of a latent process consisting of an additive expansion of a smooth Gaussian process with a fractional Brownian motion. We develop a Bayesian approach to inference using Markov chain Monte Carlo, and show good performance in simulation studies. Applying the methods to the Amazon bird vocalization data, we find substantial evidence for self-similarity and non-Markovian/Poisson dynamics. To accommodate the bird vocalization data, in which there are many different species of birds exhibiting their own vocalization dynamics, a hierarchical expansion of FRAP is provided in Supplementary Materials.


title: Composition, Size, and Surface Functionalization dependent Optical Properties of Lead Bromide Perovskite Nanocrystals

abstract: The photoluminescence (PL), color purity, and stability of lead halide perovskite nanocrystals depend critically on the surface passivation. We present a study on the temperature dependent PL and PL decay dynamics of lead bromide perovskite nanocrystals characterized by different types of A cations, surface ligands, and nanocrystal sizes. Throughout, we observe a single emission peak from cryogenic to ambient temperature. The PL decay dynamics are dominated by the surface passivation, and a post-synthesis ligand exchange with a quaternary ammonium bromide (QAB) results in a more stable passivation over a larger temperature range. The PL intensity is highest from 50K-250K, which indicates that the ligand binding competes with the thermal energy at ambient temperature. Despite the favorable PL dynamics of nanocrystals passivated with QAB ligands (monoexponential PL decay over a large temperature range, increased PL intensity and stability), the surface passivation still needs improvement toward increased emission intensity in nanocrystal films.


title: Finding Small Proofs for Description Logic Entailments: Theory and Practice (Extended Technical Report)

abstract: Logic-based approaches to AI have the advantage that their behaviour can in principle be explained by providing their users with proofs for the derived consequences. However, if such proofs get very large, then it may be hard to understand a consequence even if the individual derivation steps are easy to comprehend. This motivates our interest in finding small proofs for Description Logic (DL) entailments. Instead of concentrating on a specific DL and proof calculus for this DL, we introduce a general framework in which proofs are represented as labeled, directed hypergraphs, where each hyperedge corresponds to a single sound derivation step. On the theoretical side, we investigate the complexity of deciding whether a certain consequence has a proof of size at most $n$ along the following orthogonal dimensions: (i) the underlying proof system is polynomial or exponential; (ii) proofs may or may not reuse already derived consequences; and (iii) the number $n$ is represented in unary or binary. We have determined the exact worst-case complexity of this decision problem for all but one of the possible combinations of these options. On the practical side, we have developed and implemented an approach for generating proofs for expressive DLs based on a non-standard reasoning task called forgetting. We have evaluated this approach on a set of realistic ontologies and compared the obtained proofs with proofs generated by the DL reasoner ELK, finding that forgetting-based proofs are often better w.r.t. different measures of proof complexity.


title: Identification of deregulated transcription factors involved in subtypes of cancers

abstract: We propose a methodology for the identification of transcription factors involved in the deregulation of genes in tumoral cells. This strategy is based on the inference of a reference gene regulatory network that connects transcription factors to their downstream targets using gene expression data. The behavior of genes in tumor samples is then carefully compared to this network of reference to detect deregulated target genes. A linear model is finally used to measure the ability of each transcription factor to explain those deregulations. We assess the performance of our method by numerical experiments on a breast cancer data set. We show that the information about deregulation is complementary to the expression data as the combination of the two improves the supervised classification performance of samples into cancer subtypes.


title: Traveling fronts in self-replicating persistent random walks with multiple internal states

abstract: Self-activation coupled to a transport mechanism results in traveling waves that describe polymerization reactions, forest fires, tumor growth, and even the spread of epidemics. Diffusion is a simple and commonly used model of particle transport. Many physical and biological systems are, however, better described by persistent random walks that switch between multiple states of ballistic motion. So far, traveling fronts in persistent random walk models have only been analyzed in special, simplified cases. Here, we formulate the general model of reaction-transport processes in such systems and show how to compute the expansion velocity for arbitrary number of states. For the two-state model, we obtain a closed-form expression for the velocity and report how it is affected by different transport and replication parameters. We also show that nonzero death rates result in a discontinuous transition from quiescence to propagation. We compare our results to a recent observation of a discontinuous onset of propagation in microtubule asters and comment on the universal nature of the underlying mechanism.

categories: cond-mat.stat-mech nlin.PS q-bio.PE q-bio.SC


title: Electromigration-guided composition patterns in thin alloy films: a computational study

abstract: Via computation of a continuum dynamical model of the diffusion and electromigration, this paper demonstrates the feasibility of guiding the formation of the stripe composition patterns in the thin surface layers of the crystal alloy films. By employing the systematic parametric computational analysis it is revealed how such properties of the pattern as the aerial number density of the stripes and the stripe in-plane orientation are influenced by the major physical factors that are not limited to the electric field strength and its direction angle in the plane, but also include a number of parameters that originate in the anisotropy of diffusion in the particular crystallographically-oriented surface layer. By following the insights from this analysis the real patterns hopefully can be created in a dedicated experiment.

categories: cond-mat.mtrl-sci cond-mat.mes-hall nlin.AO nlin.PS


title: Spanning the isogeny class of a power of an ordinary elliptic curve over a finite field. Application to the number of rational points of curves of genus $\leq 4$

abstract: Let $E$ be an ordinary elliptic curve over a finite field and $g$ be a positive integer. Under some technical assumptions, we give an algorithm to span the isomorphism classes of principally polarized abelian varieties in the isogeny class of $E^g$. The varieties are first described as hermitian lattices over (not necessarily maximal) quadratic orders and then geometrically in terms of their algebraic theta null point. We also show how to algebraically compute Siegel modular forms of even weight given as polynomials in the theta constants by a careful choice of an affine lift of the theta null point. We then use these results to give an algebraic computation of Serre's obstruction for principally polarized abelian threefolds isogenous to $E^3$ and of the Igusa modular form in dimension $4$. We illustrate our algorithms with examples of curves with many rational points over finite fields.


title: On Zeckendorf Related Partitions Using the Lucas Sequence

abstract: Zeckendorf proved that every positive integer has a unique partition as a sum of non-consecutive Fibonacci numbers. Similarly, every natural number can be partitioned into a sum of non-consecutive terms of the Lucas sequence, although such partitions need not be unique. In this paper, we prove that a natural number can have at most two distinct non-consecutive partitions in the Lucas sequence, find all positive integers with a fixed term in their partition, and calculate the limiting value of the proportion of natural numbers that are not uniquely partitioned into the sum of non-consecutive terms in the Lucas sequence.


title: Index Modulation-Based Flexible Non-Orthogonal Multiple Access

abstract: Non-orthogonal multiple access (NOMA) is envisioned as an efficient candidate for future communication systems. This letter proposes a novel orthogonal frequency division multiplexing (OFDM) with index modulation (IM)-based NOMA scheme, called OFDM-IM NOMA, for future multi-user communication systems. Inspired by IM and classical OFDM-NOMA, users utilize flexibility by adjusting power allocation factors and subcarrier activation ratios. Our new scheme allows different service users to share available resources as in classical NOMA, more efficiently. It is shown that OFDM-IM NOMA reliably supports a high and low data rate user at the same resources by adjusting their subcarrier activation ratios.


title: Causal Inference under Outcome-Based Sampling with Monotonicity Assumptions

abstract: We study causal inference under case-control and case-population sampling. Specifically, we focus on the binary-outcome and binary-treatment case, where the parameters of interest are causal relative and attributable risks defined via the potential outcome framework. It is shown that strong ignorability is not always as powerful as it is under random sampling and that certain monotonicity assumptions yield comparable results in terms of sharp identified intervals. Specifically, the usual odds ratio is shown to be a sharp identified upper bound on causal relative risk under the monotone treatment response and monotone treatment selection assumptions. We offer algorithms for inference on the causal parameters that are aggregated over the true population distribution of the covariates. We show the usefulness of our approach by studying three empirical examples: the benefit of attending private school for entering a prestigious university in Pakistan; the relationship between staying in school and getting involved with drug-trafficking gangs in Brazil; and the link between physicians' hours and size of the group practice in the United States.


title: Magnetic correlations in subsystems of the misfit [Ca$_2$CoO$_3$]$_{0.62}$[CoO$_2$] cobaltate

abstract: [Ca$_2$CoO$_3$]$_{0.62}$[CoO$_2$], a two dimensional misfit metallic compound, is famous for its rich phases accessed by temperature, $i.e.$ high temperature spin-state transition, metal-insulator transition (MIT) at intermediate temperature ($\sim$ 100 K) and low temperature spin density wave (SDW). It enters into SDW phase below T$_{MIT}$ which becomes long range at 27 K. Information on the independent role of misfit layers (rocksalt/Ca$_2$CoO$_3$ \& triangular/CoO$_2$) in these phases is scarce. By combining a set of complementary macroscopic (DC magnetization and resistivity) and microscopic (neutron diffraction and X-ray absorption fine structure spectroscopy) measurements on pure (CCO) and Tb substituted in the rocksalt layer of CCO (CCO1), magnetic correlations in both subsystems of this misfit compound are unraveled. CCO is found to exhibit glassiness, as well as exchange bias (EB) effects, while CCO1 does not exhibit glassiness, albeit it shows weaker EB effect. By combining local structure investigations from extended X-ray absorption fine structure (EXAFS) spectroscopy and neutron diffraction results on CCO, we confirm that the SDW arises in the CoO$_2$ layer. Our results show that the magnetocrystalline anisotropy associated with the rocksalt layer acts as a source of pinning, which is responsible for EB effect. Ferromagnetic clusters in the Ca$_2$CoO$_3$ affects SDW in CoO$_2$ and ultimately glassiness arises.


title: A Computational Model of Levodopa-Induced Toxicity in Substantia Nigra Pars Compacta in Parkinson's Disease

abstract: Parkinson's disease (PD) is caused by the progressive loss of dopaminergic cells in substantia nigra pars compacta (SNc). The root cause of this cell loss in PD is still not decisively elucidated. A recent line of thinking traces the cause of PD neurodegeneration to metabolic deficiency. Due to exceptionally high energy demand, SNc neurons exhibit a higher basal metabolic rate and higher oxygen consumption rate, which results in oxidative stress. Recently, we have suggested that the excitotoxic loss of SNc cells might be due to energy deficiency occurring at different levels of neural hierarchy. Levodopa (LDOPA), a precursor of dopamine, which is used as a symptom-relieving treatment for PD, leads to outcomes that are both positive and negative. Several researchers suggested that LDOPA might be harmful to SNc cells due to oxidative stress. The role of LDOPA in the course of PD pathogenesis is still debatable. We hypothesize that energy deficiency can lead to LDOPA-induced toxicity (LIT) in two ways: by promoting dopamine-induced oxidative stress and by exacerbating excitotoxicity in SNc. We present a multiscale computational model of SNc-striatum system, which will help us in understanding the mechanism behind neurodegeneration postulated above and provides insights for developing disease-modifying therapeutics. It was observed that SNc terminals are more vulnerable to energy deficiency than SNc somas. During LDOPA therapy, it was observed that higher LDOPA dosage results in increased loss of somas and terminals in SNc. It was also observed that co-administration of LDOPA and glutathione (antioxidant) evades LDOPA-induced toxicity in SNc neurons. We show that our proposed model was able to capture LDOPA-induced toxicity in SNc, caused by energy deficiency.


title: The local-filament pattern in the anomalous transparency of the Universe for energetic gamma rays

abstract: The propagation length of high-energy photons through the Universe is limited by $e^+e^-$ pair production on the extragalactic background radiation. Previous studies reported discrepancies between predicted and observed attenuation, suggesting explanations in terms of new physics. However, these effects are dominated by a limited number of observed sources, while many do not show any discrepancy. Here, we consider the distribution in the sky of these apparently anomalous objects, selected in two very different approaches: the study of unphysical hardenings at distance-dependent energies in deabsorbed spectra of TeV blazars, and the observation of ultra-high-energy air showers from the directions of BL Lac type objects. In both cases, directions to the anomalous sources follow the projected local distribution of galaxies: all the distant sources, contributing to the anomalies, are seen through the local filament. This matches the prediction of the proposed earlier explanation of the anomalies based on mixing of photons with axion-like particles in the filament's magnetic field. For ultra-high energies, this axion interpretation may be tested by the search of primary gamma rays.


title: Coping with dating errors in causality estimation

abstract: {We consider the problem of estimating causal influences between observed processes from time series possibly corrupted by errors in the time variable (dating errors) which are typical in palaeoclimatology, planetary science and astrophysics. "Causality ratio" based on the Wiener -- Granger causality is proposed and studied for a paradigmatic class of model systems to reveal conditions under which it correctly indicates directionality of unidirectional coupling. It is argued that in case of {\it a priori} known directionality, the causality ratio allows a characterization of dating errors and observational noise. Finally, we apply the developed approach to palaeoclimatic data and quantify the influence of solar activity on tropical Atlantic climate dynamics over the last two millennia. A stronger solar influence in the first millennium A.D. is inferred. The results also suggest a dating error of about 20 years in the solar proxy time series over the same period.


title: A polynomial-degree-robust a posteriori error estimator for N\'ed\'elec discretizations of magnetostatic problems

abstract: We present an equilibration-based a posteriori error estimator for N\'ed\'elec element discretizations of the magnetostatic problem. The estimator is obtained by adding a gradient correction to the estimator for N\'ed\'elec elements of arbitrary degree presented in [J. Gedicke, S. Geevers, and I. Perugia. An equilibrated a posteriori error estimator for arbitrary-order N\'ed\'elec elements for magnetostatic problems. Journal of Scientific Computing, 83:1-23, 2020]. This new estimator is proven to be reliable, with reliability constant 1, and efficient, with an efficiency constant that is independent of the polynomial degree of the approximation. These properties are demonstrated in a series of numerical experiments on three-dimensional test problems.


title: Hitting forbidden induced subgraphs on bounded treewidth graphs

abstract: For a fixed graph $H$, the $H$-IS-Deletion problem asks, given a graph $G$, for the minimum size of a set $S \subseteq V(G)$ such that $G\setminus S$ does not contain $H$ as an induced subgraph. Motivated by previous work about hitting (topological) minors and subgraphs on bounded treewidth graphs, we are interested in determining, for a fixed graph $H$, the smallest function $f_H(t)$ such that $H$-IS-Deletion can be solved in time $f_H(t) \cdot n^{O(1)}$ assuming the Exponential Time Hypothesis (ETH), where $t$ and $n$ denote the treewidth and the number of vertices of the input graph, respectively. We show that $f_H(t) = 2^{O(t^{h-2})}$ for every graph $H$ on $h \geq 3$ vertices, and that $f_H(t) = 2^{O(t)}$ if $H$ is a clique or an independent set. We present a number of lower bounds by generalizing a reduction of Cygan et al. [MFCS 2014] for the subgraph version. In particular, we show that when $H$ deviates slightly from a clique, the function $f_H(t)$ suffers a sharp jump: if $H$ is obtained from a clique of size $h$ by removing one edge, then $f_H(t) = 2^{\Theta(t^{h-2})}$. We also show that $f_H(t) = 2^{\Omega(t^{h})}$ when $H=K_{h,h}$, and this reduction answers an open question of Mi. Pilipczuk [MFCS 2011] about the function $f_{C_4}(t)$ for the subgraph version. Motivated by Cygan et al. [MFCS 2014], we also consider the colorful variant of the problem, where each vertex of $G$ is colored with some color from $V(H)$ and we require to hit only induced copies of $H$ with matching colors. In this case, we determine, under the ETH, the function $f_H(t)$ for every connected graph $H$ on $h$ vertices: if $h\leq 2$ the problem can be solved in polynomial time; if $h\geq 3$, $f_H(t) = 2^{\Theta(t)}$ if $H$ is a clique, and $f_H(t) = 2^{\Theta(t^{h-2})}$ otherwise.